Additional Information on Load Sharing

For load sharing systems, one way to model the effect of load on life is to use a life distribution along with a life-stress relationship (as discussed in A Brief Introduction to Life-Stress Relationships) for each component. To illustrate this, consider a system of two units connected reliability-wise in parallel as shown below.

Assume that the units must supply an output of 8 volts and that if both units are operational, each unit is to supply 50% of the total output. If one of the units fails, then the surviving unit supplies 100%. Furthermore, assume that having to supply the entire load has a negative impact on the reliability characteristics of the surviving unit.

We will create the model starting from raw data. Assume that a total of 20 units were tested to failure at 7, 10 and 15 volts. The test data set is presented in the next table.

For this example, Units 1 and 2 are the same component. Therefore, only one set of data was collected. However, it is possible that the load sharing components in a system may not be the same. If that were the case, data would need to be collected for each component.

The data set was analyzed using ReliaSoft's ALTA software (as shown in the figure below) with the Inverse Power Law as the underlying life-stress relationship and Weibull as the life distribution.



The estimated model parameters, $$\beta \,\!$$, $$K\,\!$$ and $$n\,\!$$, are shown next.


 * $$\begin{align}

\beta = & 1.9239 \\ K= & 3.2387\times {{10}^{-7}} \\ n= & 3.4226 \end{align}\,\!$$

Or:


 * $$\begin{align}

{{R}_{1}}(t,{{S}_{1}})= & {{e}^{-{{\left( KS_{1}^{n}t \right)}^{\beta }}}} \\ = & {{e}^{-{{\left( 3.2387\times {{10}^{-7}}S_{1}^{3.4226}t \right)}^{1.9239}}}} \ \end{align}\,\!$$.


 * $${{f}_{1}}(t,{{S}_{1}})=\beta KS_{1}^{n}{{\left( KS_{1}^{n}t \right)}^{\beta -1}}{{e}^{-{{\left( KS_{1}^{n}t \right)}^{\beta }}}} \ \,\!$$

And for this case:


 * $$\begin{align}

{{R}_{1}}(t,{{S}_{1}})= & {{R}_{2}}(t,{{S}_{2}}) \\ {{f}_{1}}(t,{{S}_{1}})= & {{f}_{2}}(t,{{S}_{2}}) \end{align}\,\!$$

The figure below shows a plot of $${{R}_{1}}(t,{{S}_{1}})={{e}^{-{{\left( KS_{1}^{n}t \right)}^{\beta }}}} = {{e}^{-{{\left( 3.2387\times {{10}^{-7}}S_{1}^{3.4226}t \right)}^{1.9239}}}}\ \,\!$$.



Now that the failure properties have been determined using the test data, the reliability of the system at some time, $$t\,\!$$, can be calculated using the following equation:


 * $$\begin{align}

R(t,S)= & {{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}}) \\ & +\underset{o}{\overset{t}{\mathop \int }}\,{{f}_{1}}\left( x,{{S}_{1}} \right)\cdot {{R}_{2}}(x,{{S}_{2}})\cdot \left( \frac{{{R}_{2}}({{t}_{1e}}+(t-x),S)}{{{R}_{2}}({{t}_{1e}},S)} \right)dx \\ & +\underset{o}{\overset{t}{\mathop \int }}\,{{f}_{2}}\left( x,{{S}_{2}} \right)\cdot {{R}_{1}}(x,{{S}_{1}})\cdot \left( \frac{{{R}_{1}}({{t}_{2e}}+(t-x),S)}{{{R}_{1}}({{t}_{2e}},S)} \right)dx \end{align}\,\!$$

where:


 * $$\begin{align}

{{S}_{1}}= & {{P}_{1}}S \\ {{S}_{2}}= & {{P}_{2}}S \end{align}\,\!$$

and:


 * $$S\,\!$$ is the total load (or required output).


 * $${{P}_{1}}\,\!$$ and $${{P}_{2}}\,\!$$ are the portion of the total load that each unit supports when both units are operational. In this case, $${{P}_{1}}={{P}_{2}}=0.5=50%.\,\!$$


 * $${{S}_{1}}\,\!$$ and $${{S}_{2}}\,\!$$ represent the portions of the load that Unit 1 and Unit 2 must support when both units are operational.


 * $${{t}_}\,\!$$ is the equivalent operating time for Unit 1 if it had been operating at $$S\,\!$$ instead of $${{S}_{1}}\,\!$$. A graphical representation of the equivalent time is shown in the following figure, where the curve marked by L represents the low stress (load) and the curve marked by H represents the high stress (load).




 * $${{t}_{1e}}\,\!$$ can be calculated by:


 * $$\begin{align}

{{R}_{1}}(t)= & {{R}_{1}}({{t}_{1e}}) \\ {{e}^{-{{(tKS_{1}^{n})}^{\beta }}}}= & {{e}^{-{{({{t}_{1e}}K{{S}^{n}})}^{\beta }}}} \\ tS_{1}^{n}= & {{t}_{1e}}{{S}^{n}} \\ {{t}_{1e}}= & t{{\left( \frac{{{S}_{1}}}{S} \right)}^{n}},\text{    }{{S}_{1}}={{P}_{1}}S \\ \therefore & {{t}_{1e}}=tP_{1}^{n} \end{align}\,\!$$


 * $${{t}_{2e}}\,\!$$ can be calculated the same way, or:


 * $$\begin{align}

{{R}_{2}}(t)= & {{R}_{2}}({{t}_{2e}}) \\ \therefore & {{t}_{2e}}=tP_{2}^{n} \end{align}\,\!$$

In this example, the reliability equations for Unit 1 and Unit 2 are the same since they are the same type of component and demonstrate the same failure properties. In addition, the total output is divided equally between the two units (when both units are operating), so $${{t}_{1e}}\,\!$$ and $${{t}_{2e}}\,\!$$ will also be the same.

The next step is to determine the reliability of the system after 8,760 hours, $$R(t=8,760)\,\!$$. Using
 * $$\begin{align}

R(t,S)= & {{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}}) \\ & +\underset{o}{\overset{t}{\mathop \int }}\,{{f}_{1}}\left( x,{{S}_{1}} \right)\cdot {{R}_{2}}(x,{{S}_{2}})\cdot \left( \frac{{{R}_{2}}({{t}_{1e}}+(t-x),S)}{{{R}_{2}}({{t}_{1e}},S)} \right)dx \\ & +\underset{o}{\overset{t}{\mathop \int }}\,{{f}_{2}}\left( x,{{S}_{2}} \right)\cdot {{R}_{1}}(x,{{S}_{1}})\cdot \left( \frac{{{R}_{1}}({{t}_{2e}}+(t-x),S)}{{{R}_{1}}({{t}_{2e}},S)} \right)dx \end{align}\,\!$$

the system reliability is found to be:


 * $$\begin{align}

R(t=8760)= & 0.8567 \\ = & 85.67% \end{align}\,\!$$

=Load Sharing= As presented in earlier chapters, a reliability block diagram (RBD) allows you to graphically represent how the components within a system are reliability-wise connected. In most cases, independence is assumed across the components within the system. For example, the failure of component A does not affect the failure of component B. However, if a system consists of components that are sharing a load, then the assumption of independence no longer holds true.

If one component fails, then the component(s) that are still operating will have to assume the failed unit's portion of the load. Therefore, the reliabilities of the surviving unit(s) will change. Calculating the system reliability is no longer an easy proposition. In the case of load sharing components, the change of the failure distributions of the surviving components must be known in order to determine the system's reliability.

To illustrate this, consider a system of two units connected reliability-wise in parallel as shown below. Assume that the units must supply an output of 8 volts and that if both units are operational, each unit is to supply 50% of the total output. If one of the units fails, then the surviving unit supplies 100%. Furthermore, assume that having to supply the entire load has a negative impact on the reliability characteristics of the surviving unit.

Because the reliability characteristics of the unit change based on the load it is sharing, a method that can model the effect of the load on life should be used. One way to do this is to use a life distribution along with a life-stress relationship (as discussed in A Brief Introduction to Life-Stress Relationships) for each component. The detailed discussion for this method can be found at Additional Information on Load Sharing. Another simple way is to use the concept of acceleration factors and assume that the load has a linear effect on the failure time. If the load is doubled, then the life of the component will be shortened by half.

For the above load sharing system, the reliability of each component is a function of time and load. For example, for Unit 1, the reliability and the probability density function are:


 * $${{R}_{1}}(t,{{S}_{1}})\,\!\,\!$$ and $${{f}_{1}}(t,{{S}_{1}})\,\!\,\!$$

where $${{S}_{1}}\,\!\,\!$$ is the load shared by Unit 1 at time t and the total load of the system is $$S={{S}_{1}}+{{S}_{2}}\,\!\,\!$$. At the beginning, both units are working. Assume that Unit 1 fails at time x and Unit 2 takes over the entire load. The reliability for Unit 2 at time x is:


 * $$\begin{align}

& {{R}_{2}}(x,{{S}_{2}})={{R}_{2}}({{t}_{2e}},S) \\ & {{t}_{2e}}=\frac{S}x \\ \end{align}\,\!\,\!$$

$${{t}_{2e}}\,\!\,\!$$ is the equivalent time for Unit 2 at time x if it is operated with load S. The equivalent time concept is illustrated in the following plot. The system reliability at time t is:

& R(t)={{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}}) \\ & +\int_{0}^{t}{{{f}_{1}}(x,{{S}_{1}})\cdot {{R}_{2}}({{t}_{2e}}+(t-x),S)dx} \\ & +\int_{0}^{t}{{{f}_{2}}(x,{{S}_{2}})\cdot {{R}_{1}}({{t}_{1e}}+(t-x),S)dx} \end{align}\,\!\,\!$$
 * $$\begin{align}

In BlockSim, the failure time distribution for each component is defined at the load of S. The reliability function for a component at a given load is calculated as:


 * $${{R}_{i}}(t,{{S}_{i}})={{R}_{i}}(t\times \frac{S},S)\,\!\,\!$$
 * $${{f}_{i}}(t,{{S}_{i}})=\frac{S}{{f}_{i}}\left( t\times \frac{S},S \right)\,\!\,\!$$

From the above equation, it can be seen that the concept used in the calculation for load sharing is the same as the concept used in the calculation for duty cycle.

Example
In the following load sharing system, Block 1 follows a Weibull failure distribution with$${{\beta }_{1}}=1.5\,\!$$, and$${{\eta }_{1}}=1,000\,\!$$ . Block 2 follows a Weibull failure distribution with$${{\beta }_{2}}=2\,\!$$, and$${{\eta }_{2}}=2,000\,\!$$. The load for Block 1 is 1 unit, and for Block 2 it is 3 units. Calculate the system reliability at time 1,500.

Block 1 shares 25% (P1) of the entire load, and Block 2 shares 75% (P2) of it. Therefore, we have the following equations for calculating the system reliability: $${{R}_{i}}(t,{{S}_{i}})={{R}_{i}}({{P}_{i}}\times t)\,\!$$, 	$${{f}_{i}}(x,{{S}_{i}})={{P}_{i}}{{f}_{i}}({{P}_{i}}\times x)\,\!$$ and $${{R}_{1}}({{t}_{1e}}+(t-x))={{R}_{1}}({{P}_{1}}x+t-x)={{R}_{1}}(t-{{P}_{2}}x)\,\!$$ $${{R}_{2}}({{t}_{2e}}+(t-x))={{R}_{2}}({{P}_{2}}x+t-x)={{R}_{2}}(t-{{P}_{1}}x)\,\!$$ Using the above equations in the system reliability function, we get: $$\begin{align} & R(t)={{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}})+\int_{0}^{t}{{{f}_{1}}(x,{{S}_{1}})\cdot {{R}_{2}}({{t}_{2e}}+(t-x),S)dx} \\ & +\int_{0}^{t}{{{f}_{2}}(x,{{S}_{2}})\cdot {{R}_{1}}({{t}_{1e}}+(t-x),S)dx} \\ & ={{R}_{1}}({{P}_{1}}\times t)\cdot {{R}_{2}}({{P}_{2}}\times t)+\int_{0}^{t}{{{P}_{1}}{{f}_{1}}({{P}_{1}}x)\cdot {{R}_{2}}(t-{{P}_{1}}x)dx} \\ & +\int_{0}^{t}{{{P}_{2}}{{f}_{2}}({{P}_{2}}x)\cdot {{R}_{1}}(t-{{P}_{2}}x)dx} \end{align}\,\!$$ The calculated system reliability at time 1,500 is 0.8569, as given below.