Template:Stress-Strength Confidence Intervals: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
(Redirected page to Stress-Strength Analysis)
 
(15 intermediate revisions by 2 users not shown)
Line 1: Line 1:
===Confidence Intervals on the Probability===
#REDIRECT [[Stress-Strength_Analysis]]
Both the stress and strength distributions can be either 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 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, <math>R</math> using:
<center>
<math>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</math> </center>
 
Variance of <math>{{f}_{1}}(x)</math> and <math>{{R}_{2}}(x)</math> can be estimated from the Fisher Information Matrix. For detail, please see chapter [[Confidence Bounds]].
 
Once the variance of the expected reliability is obtained, the two-sided confidence interval of it can be calcualted using:
 
<center><math>[\frac{R}{R+(1-R)w},\frac{R}{R+(1-R)/w}]</math></center>
 
where:
:: CL is the confidence level;
:: <math>\alpha</math> = 1-CL;
:: <math>w=\exp \{{{z}_{1-\alpha /2}}\sqrt{Var(R)}/[R(1-R)]\}</math>;
:: <math>{{Z}_{1-\alpha /2}}</math> is the <math>1-\alpha/2</math> percential of a standard normal distribution.
 
If the upper bound (U) and lower bound (L) are not infinite and 0, the above calcualted variance of ''R'' is adjusted by <math>{{\left[ {1}/{\left( {{F}_{1}}(U)-{{F}_{1}}(L) \right)}\; \right]}^{2}}</math>.
 
 
'''Variation in Probability Values'''
 
Assume the distributions for stress and strength are known. From the stress-strength equation:
 
<center><math>R=P[Stress\le Strength]=\int_{0}^{\infty }{{{f}_{Stress}}(x)\cdot {{R}_{Strength}}(x)}dx</math></center>
 
we know, the calcualted 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 exampe. Let's first assume stress is a fixed value of 567. The reliability then is:
<center>
<math>R(567)=\Pr (Strength>567)={{R}_{2}}(567)</math>
</center>
This is the reliability when the stress value is 567 and when the strength distribution is given. If stress is not a fixed value, instead it follows a distribution, then it can take values other than 567. For instance, it can be a value of 700. Therefore, we get another reliability value of <math>{R}(700)</math>. Since stress is a random variable, for any stress value <math>{x}_{i}</math>, there is a reliability value of
<math>R({{x}_{i}})</math> calculated from the strength distribution. We will end up with many <math>R({{x}_{i}})</math>s or <math>R_{2}({{x}_{i}})</math>s. From these <math>R({{x}_{i}})</math>s, we can get the mean and variance of the reliability. In fact, its mean is the the result from:
 
<center><math>R=P[Stress\le Strength]=\int_{0}^{\infty }{{{f}_{Stress}}(x)\cdot {{R}_{Strength}}(x)}dx</math></center>
 
and its variacne is:
 
<center>
<math>\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}</math></center>
 
where ''R'' is the expected value of the reliability.
 
Once the variance of the expected reliability is obtained, the two-sided confidence interval of it can be calcualted using:
 
<center><math>[\frac{R}{R+(1-R)w},\frac{R}{R+(1-R)/w}]</math></center>
 
where:
:: CL is the confidence level;
:: <math>\alpha</math> = 1-CL;
:: <math>w=\exp \{{{z}_{1-\alpha /2}}\sqrt{Var(R)}/[R(1-R)]\}</math>;
:: <math>{{Z}_{1-\alpha /2}}</math> is the <math>1-\alpha/2</math> percential of a standard normal distribution.
 
If the upper bound (U) and lower bound (L) are not infinite and 0, the above calcualted variance of ''R'' is adjusted by <math>{{\left[ {1}/{\left( {{F}_{1}}(U)-{{F}_{1}}(L) \right)}\; \right]}^{2}}</math>.
 
 
'''Example 1:'''
{{Example: Stress-Strength Analysis with Parameter Uncertainty}}

Latest revision as of 06:24, 13 August 2012