Stress-Strength Analysis

Stress-strength analysis has been used in mechanical component design. The probability of failure is based on the probability of stress exceeding strength. The following equation is used to calculate the expected probability of failure, F:


 * $$F=P[Stress\ge Strength]=\int_{0}^{\infty }{{{f}_{Strength}}(x)\cdot {{R}_{Stress}}(x)}dx$$

The expected probability of success or the expected reliability, R, is calculated as:


 * $$R=P[Stress\le Strength]=\int_{0}^{\infty }{{{f}_{Stress}}(x)\cdot {{R}_{Strength}}(x)}dx$$

The above calculation assumes both stress and strength are in the positive domain. For general cases, the expected reliability can be calculated using the following equation:


 * $$R=P[{{X}_{1}}\le {{X}_{2}}]=\frac{1}{{{F}_{1}}(U)-{{F}_{1}}(L)}\int_{L}^{U}{{{f}_{1}}(x)\cdot {{R}_{2}}(x)}dx$$

where:
 * $$L\le {{X}_{1}}\le U$$
 * $$\begin{align}

& {{X}_{1}}:\text{ Stress } \\ & {{X}_{2}}:\text{ Strength } \\ \end{align}$$ When U = infinite and L = 0, the above two equations are the same.

Confidence Intervals on the Probability
Both the stress and strength distributions can be estimated from actual data or specified by engineers based on engineering knowledge or existing references. Based on the source of the distribution, there are two types of variation associated with the calculated probability: variation in the model parameters and variation in the probability values. Both are described next.

Variation in Model Parameters
If both the stress and strength distributions are estimated from data sets, then there are uncertainties associated with the estimated distribution parameters. These uncertainties will cause some degree of variation of the probability calculated from the stress-strength analysis. Therefore, we can use these uncertainties to estimate the confidence intervals on the calculated probability. To get the confidence intervals, we first calculate the variance of the reliability, $$R$$, using:


 * $$Var\left[ R \right]=\int_{0}^{\infty }{Var\left[ {{f}_{1}}(x) \right]}{{\left[ {{R}_{2}}(x) \right]}^{2}}dx+{{\int_{0}^{\infty }{\left[ {{f}_{1}}(x) \right]}}^{2}}Var\left[ {{R}_{2}}(x) \right]dx$$

Variance of $${{f}_{1}}(x)$$ and $${{R}_{2}}(x)$$ can be estimated from the Fisher Information Matrix. For details, please see Confidence Bounds.

Once the variance of the expected reliability is obtained, the two-sided confidence intervals can be calculated using:


 * $$[\frac{R}{R+(1-R)w},\frac{R}{R+(1-R)/w}]$$

where:
 * CL is the confidence level
 * $$\alpha$$ is 1-CL
 * $$w=\exp \{{{z}_{1-\alpha /2}}\sqrt{Var(R)}/[R(1-R)]\}$$
 * $${{Z}_{1-\alpha /2}}$$ is the $$1-\alpha/2$$ percentile of a standard normal distribution.

If the upper bound (U) and lower bound (L) are not infinite and 0, respectively, then the calculated variance of $$ R $$ is adjusted by $${{\left[ {1}/{\left( {{F}_{1}}(U)-{{F}_{1}}(L) \right)}\; \right]}^{2}}$$.

Variation in Probability Values
Assume the distributions for stress and strength are known. From the stress-strength equation


 * $$R=P[Stress\le Strength]=\int_{0}^{\infty }{{{f}_{Stress}}(x)\cdot {{R}_{Strength}}(x)}dx$$

we know that the calculated reliability is the expected value of the probability that a strength value is larger than a stress value. Since both strength and stress are random variables from their distributions, the reliability is also a random variable. This can be explained using the following example. Let's first assume that stress is a fixed value of 567. The reliability then is:


 * $$R(567)=\Pr (Strength>567)={{R}_{2}}(567)$$

This is the reliability when the stress value is 567 and when the strength distribution is given. If stress is not a fixed value (i.e., it follows a distribution instead), then it can take values other than 567. For instance, if stress takes a value of 700, then we get another reliability value of $${R}(700)$$. Since stress is a random variable, for any stress value $${x}_{i}$$, there is a reliability value of $$R({{x}_{i}})$$ calculated from the strength distribution. We will end up with many $$R({{x}_{i}})$$s or $$R_{2}({{x}_{i}})$$s. From these $$R({{x}_{i}})$$s, we can get the mean and variance of the reliability. In fact, its mean is the result from:


 * $$R=P[Stress\le Strength]=\int_{0}^{\infty }{{{f}_{Stress}}(x)\cdot {{R}_{Strength}}(x)}dx$$

and its variance is:


 * $$\begin{align}

& Var\left[ R \right]=Var\left[ {{R}_{2}}({{X}_{1}}) \right]=E\left[ {{R}_{2}}{{\left( {{X}_{1}} \right)}^{2}} \right]-{{\left( E\left[ {{R}_{2}}\left( {{X}_{1}} \right) \right] \right)}^{2}} \\ & =\int_{0}^{\infty }{{{f}_{1}}(x){{\left[ {{R}_{2}}(x) \right]}^{2}}dx}-{{\left( E\left[ {{R}_{2}}\left( {{X}_{1}} \right) \right] \right)}^{2}} \\ & =\int_{0}^{\infty }{{{f}_{1}}(x){{\left[ {{R}_{2}}(x) \right]}^{2}}dx}-{{\left( R \right)}^{2}} \\ \end{align}$$

where R is the expected value of the reliability.

Once the variance of the expected reliability is obtained, the two-sided confidence intervals can be calculated using:


 * $$[\frac{R}{R+(1-R)w},\frac{R}{R+(1-R)/w}]$$

where:
 * CL is the confidence level
 * $$\alpha$$ is 1-CL
 * $$w=\exp \{{{z}_{1-\alpha /2}}\sqrt{Var(R)}/[R(1-R)]\}$$
 * $${{Z}_{1-\alpha /2}}$$ is the $$1-\alpha/2$$ percentile of a standard normal distribution

If the upper bound (U) and lower bound (L) are not infinite and 0, the above calculated variance of R is adjusted by $${{\left[ {1}/{\left( {{F}_{1}}(U)-{{F}_{1}}(L) \right)}\; \right]}^{2}}$$.

Stress-Strength Analysis in Design for Reliability
As we know, the expected reliability is called from the following stress-strength calculation:


 * $$R=P[Stress\le Strength]=\int_{0}^{\infty }{{{f}_{Stress}}(x)\cdot {{R}_{Strength}}(x)}dx$$

The stress distribution is usually estimated from customer usage data, such as the mileage per year of a passenger car or the load distribution for a beam. The strength distribution, on the other hand, is affected by the material used in the component, the geometric dimensions and the manufacturing process.

Because the stress distribution can be estimated from customer usage data, we can assume that $${f}_{Stress} $$ is known. Therefore, for a given reliability goal, the strength distribution $$ {R}_{Strength}$$ is the only unknown in the given equation. The factors that affect the strength distribution can be adjusted to obtain a strength distribution that meets the reliability goal. The following example shows how to use the Target Reliability Parameter Estimator tool in the stress-strength folio to obtain the parameters for a strength distribution that will meet a specified reliability goal.