Template:Gll relationship: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
 
(11 intermediate revisions by 2 users not shown)
Line 1: Line 1:
==General Log-Linear Relationship==
#REDIRECT [[Multivariable_Relationships:_General_Log-Linear_and_Proportional_Hazards]]
<br>
When a test involves multiple accelerating stresses or requires the inclusion of an engineering variable, a general multivariable relationship is needed. Such a relationship is the general log-linear relationship, which describes a life characteristic as a function of a vector of  <math>n</math>  stresses, or  <math>\underline{X}=({{X}_{1}},{{X}_{2}}...{{X}_{n}}).</math>  ALTA PRO includes this relationship and allows up to eight stresses. Mathematically the relationship is given by:
<br>
::<math>L(\underline{X})={{e}^{{{\alpha }_{0}}+\underset{j=1}{\overset{n}{\mathop{\sum }}}\,{{\alpha }_{j}}{{X}_{j}}}}</math>
 
<br>
:where:
<br>
• <math>{{\alpha }_{0}}</math>  and  <math>{{\alpha }_{j}}</math>  are model parameters.
 
• <math>X</math>  is a vector of  <math>n</math>  stresses.
 
This relationship can be further modified through the use of transformations and can be reduced to the relationships discussed previously, if so desired. As an example, consider a single stress application of this relationship and an inverse transformation on  <math>X,</math>  such that  <math>V=1/X</math>  or:
 
<br>
::<math>\begin{align}
  & L(V)= & {{e}^{{{\alpha }_{0}}+\tfrac{{{\alpha }_{1}}}{V}}} \\
& = & {{e}^{{{\alpha }_{0}}}}{{e}^{\tfrac{{{\alpha }_{1}}}{V}}} 
\end{align}</math>
 
<br>
It can be easily seen that the generalized log-linear relationship with a single stress and an inverse transformation, Eqn. (argll), has been reduced to the Arrhenius relationship discussed in Chapter 6, where:
<br>
::<math>\begin{align}
  & C= & {{e}^{{{\alpha }_{0}}}} \\
& B= & {{\alpha }_{1}} 
\end{align}</math>
 
<br>
:or:
<br>
 
::<math>L(V)=C{{e}^{\tfrac{B}{V}}}</math>
 
 
<br>
Similarly, when one chooses to apply a logarithmic transformation on  <math>X</math>  such that  <math>V=\ln (X)</math> , the relationship would reduce to the inverse power relation discussed in Chapter 8. Furthermore, if more than one stress is present, one could choose to apply a different transformation to each stress to create combination relationships similar to the ones discussed in Chapters 9 and 10. ALTA PRO has three built-in transformation options, namely:
 
<br>
{|align="center"
|-
|None|| X=V|| Exponential LSR
|-
|Reciprocal|| <math>V=1/X</math>|| Arrhenius LSR
|-
|Logarithmic|| <math>V=\ln (X)</math>|| Power LSR
|}
 
<br>
The power of the relationship and this formulation becomes evident once one realizes that 6,651 unique life-stress relationships are possible (when allowing a maximum of eight stresses). When combined with the life distributions available in ALTA PRO, almost 20,000 models can be created.
<br>
===Using the GLL Model===
<br>
Like the previous relationships, the general log-linear relationship can be combined with any of the available life distributions by expressing a life characteristic from that distribution with the GLL relationship. A brief overview of the GLL-distribution models available in ALTA PRO is presented next.
<br>
====GLL Exponential====
<br>
The GLL-exponential model can be derived by setting  <math>m=L(\underline{X})</math>  in Eqn. (GLL1), yielding the following GLL-exponential  <math>pdf</math> :
 
<br>
::<math>f(t,\underline{X})={{e}^{-\left( {{\alpha }_{0}}+\underset{j=1}{\overset{n}{\mathop{\sum }}}\,{{\alpha }_{j}}{{X}_{j}} \right)}}{{e}^{-\left( {{\alpha }_{0}}+\underset{j=1}{\overset{n}{\mathop{\sum }}}\,{{\alpha }_{j}}{{X}_{j}} \right)\cdot t}}</math>
 
<br>
The total number of unknowns to solve for in this model is  <math>n+1</math>  (i.e.  <math>{{a}_{0}},{{a}_{1}},...{{a}_{n}}).</math>
<br>
{{gll weibull}}
====GLL Lognormal====
<br>
The GLL-lognormal model can be derived by setting  <math>\breve{T}=L(\underline{X})</math> 
in Eqn. (GLL1), yielding the following GLL-lognormal <math>pdf</math> :
 
<br>
::<math>f(t,\underline{X})=\frac{1}{t\text{ }{{\sigma }_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{{T}'-{{\alpha }_{0}}-\underset{j=1}{\overset{n}{\mathop{\sum }}}\,{{\alpha }_{j}}{{X}_{j}}}{{{\sigma }_{{{T}'}}}} \right)}^{2}}}}</math>
 
<br>
The total number of unknowns to solve for in this model is  <math>n+2</math>  (i.e.  <math>{{\sigma }_{{{T}'}}},{{a}_{0}},{{a}_{1}},...{{a}_{n}}).</math>
<br>
 
====GLL Likelihood Function====
<br>
The maximum likelihood estimation method can be used to determine the parameters for the GLL relationship and the selected life distribution. For each distribution, the likelihood function can be derived, and the parameters of model (the distribution parameters and the GLL parameters) can be obtained by maximizing the log-likelihood function. For example, the log-likelihood function for the Weibull distribution is given by:
 
<br>
::<math>\begin{align}
  & \ln (L)= & \Lambda =\underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,{{N}_{i}}\ln \left[ \beta \cdot T_{i}^{\beta -1}{{e}^{-T_{i}^{\beta }\cdot {{e}^{-\beta \left( {{\alpha }_{0}}+\mathop{}_{j=1}^{n}{{a}_{j}}{{x}_{i,j}} \right)}}}}{{e}^{-\beta \left( {{\alpha }_{0}}+\mathop{}_{j=1}^{n}{{a}_{j}}{{x}_{i,j}} \right)}} \right] \\
&  & -\underset{i=1}{\overset{S}{\mathop \sum }}\,N_{i}^{\prime }{{\left( T_{i}^{\prime } \right)}^{\beta }}{{e}^{-\beta \left( {{\alpha }_{0}}+\mathop{}_{j=1}^{n}{{a}_{j}}{{x}_{i,j}} \right)}}+\overset{FI}{\mathop{\underset{i=1}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\ln [R_{Li}^{\prime \prime }-R_{Ri}^{\prime \prime }] 
\end{align}</math>
 
<br>
:where:
 
<br>
::<math>\begin{align}
  & R_{Li}^{\prime \prime }= & {{e}^{-{{\left( T_{Li}^{\prime \prime }{{e}^{{{\alpha }_{0}}+\underset{j=1}{\mathop{\overset{n}{\mathop{\mathop{}_{}^{}}}\,}}\,{{\alpha }_{j}}{{x}_{j}}}} \right)}^{\beta }}}} \\
& R_{Ri}^{\prime \prime }= & {{e}^{-{{\left( T_{Ri}^{\prime \prime }{{e}^{{{\alpha }_{0}}+\underset{j=1}{\mathop{\overset{n}{\mathop{\mathop{}_{}^{}}}\,}}\,{{\alpha }_{j}}{{x}_{j}}}} \right)}^{\beta }}}} 
\end{align}</math>
 
<br>
:and:
<br>
• <math>{{F}_{e}}</math>  is the number of groups of exact times-to-failure data points.
<br>
• <math>{{N}_{i}}</math>  is the number of times-to-failure in the  <math>{{i}^{th}}</math>  time-to-failure data group.
<br>
• <math>\lambda </math>  is the failure rate parameter (unknown).
<br>
• <math>{{T}_{i}}</math>  is the exact failure time of the  <math>{{i}^{th}}</math>  group.
<br>
• <math>S</math>  is the number of groups of suspension data points.
<br>
• <math>N_{i}^{\prime }</math>  is the number of suspensions in the  <math>{{i}^{th}}</math>  group of suspension data points.
<br>
• <math>T_{i}^{\prime }</math>  is the running time of the  <math>{{i}^{th}}</math>  suspension data group.
<br>
• <math>FI</math>  is the number of interval data groups.
<br>
• <math>N_{i}^{\prime \prime }</math>  is the number of intervals in the  <math>{{i}^{th}}</math>  group of data intervals.
<br>
• <math>T_{Li}^{\prime \prime }</math>  is the beginning of the  <math>{{i}^{th}}</math>  interval.
<br>
• <math>T_{Ri}^{\prime \prime }</math>  is the ending of the  <math>{{i}^{th}}</math>  interval.
<br>
<br>
 
===Example===
<br>
Consider the data summarized in Table 1 and Table 2. These data illustrate a typical three stress type accelerated test.
<br>
<br>
<center>Table 1: Stress Profile Summary</center>
[[Image:ALTA11t1.gif|thumb|center|300px|]]
<br>
 
<center>Table 2: Failure Data</center>
[[Image:ALTA11t2.gif|thumb|center|300px|]]
 
<br>
The data of Table 2 is analyzed assuming a Weibull distribution, an Arrhenius life-stress relationship for temperature and an inverse power life-stress relationship for voltage. No transformation is performed on the operation type. The operation type variable is treated as an indicator variable,  taking the discrete values of 0 and 1, for on/off and continuous operation, respectively. The following figure, shows the stress types and their transformations as in ALTA.
<br>
[[Image:StressTransformation.gif|thumb|center|300px|]]
 
<br>
Eqn. (GLL1) then becomes:
 
<br>
::<math>\eta ={{e}^{{{\alpha }_{0}}+{{\alpha }_{1}}\tfrac{1}{{{V}_{1}}}+{{\alpha }_{2}}\ln ({{V}_{2}})+{{\alpha }_{3}}{{V}_{3}}}}</math>
 
<br>
The resulting relationship after performing these transformations is:
 
<br>
::<math>\begin{align}
  & \eta = & {{e}^{{{\alpha }_{0}}}}{{e}^{{{\alpha }_{1}}\tfrac{1}{{{V}_{1}}}}}{{e}^{{{\alpha }_{2}}\ln ({{V}_{2}})}}{{e}^{{{\alpha }_{3}}{{V}_{3}}}} \\
& = & {{e}^{{{\alpha }_{0}}}}{{e}^{{{\alpha }_{1}}\tfrac{1}{{{V}_{1}}}}}V_{2}^{{{\alpha }_{2}}}{{e}^{{{\alpha }_{3}}{{V}_{3}}}} 
\end{align}</math>
 
<br>
Therefore, the parameter  <math>B</math>  of the Arrhenius relationship is equal to the log-linear coefficient  <math>{{\alpha }_{1}}</math> , and the parameter  <math>n</math>  of the inverse power relationship is equal to ( <math>-{{\alpha }_{2}}</math> ). Therefore  <math>\eta </math>  can also be written as:
 
<br>
::<math>\eta ={{e}^{{{\alpha }_{0}}}}{{e}^{\tfrac{B}{{{V}_{1}}}}}V_{2}^{n}{{e}^{{{\alpha }_{3}}{{V}_{3}}}}</math>
 
<br>
The activation energy of the Arrhenius relationship can be calculated by multiplying B with Boltzmann's constant.
<br>
The best fit values for the parameters in this case are:
 
<br>
::<math>\begin{align}
  & \beta = & 3.7483;\text{ }{{\alpha }_{0}}=-6.0220;\text{ }{{\alpha }_{1}}=5776.9341 \\
& {{\alpha }_{2}}= & -1.4340;\text{ }{{\alpha }_{3}}=0.6242. 
\end{align}</math>
 
<br>
Once the parameters are estimated, further analysis on the data can be performed. First, using ALTA PRO, a Weibull probability plot of the data can be obtained, as shown in Fig. 1.
<br>
 
<br>
[[Image:ALTA11.1.gif|thumb|center|300px|Weibull probability plot for all covariates.]]
<br>
<br>
Several types of information about the model as well as the data can be obtained from a probability plot. For example, the choice of an underlying distribution and the assumption of a common slope (shape parameter) can be examined. In this example, the linearity of the data supports the use of the Weibull distribution. In addition, the data appear parallel on this plot, therefore reinforcing the assumption of a common beta. Further statistical analysis can and should be performed for these purposes as well.
The Life vs. Stress plot is a very common plot for the analysis of accelerated data. Life vs. Stress plots can be very useful in assessing the effect of each stress on a product's failure. In this case, since the life is a function of three stresses, three different plots can be created. Such plots are created by holding two of the stresses constant at the desired use level, and varying the remaining one. The use stress levels for this example are 328K for temperature and 10V for voltage. For the operation type, a decision has to be made by the engineers as to whether they implement on/off or continuous operation. Figs. 2 and 3 display the effects of temperature and voltage on the life of the product.
<br>
[[Image:ALTA11.2.gif|thumb|center|300px|Effects of temperature on life.]]
<br>
<br>
[[Image:ALTA11.3.gif|thumb|center|300px|Effects of voltage on life.]]
<br>
<br>
The effects of the two different operation types on life can be observed in Fig. 4. It can be seen that the on/off cycling has a greater effect on the life of the product in terms of accelerating failure than the continuous operation. In other words, a higher reliability can be achieved by running the product continuously.
 
<br>
[[Image:ALTA11.4.gif|thumb|center|300px|Effect of operation type on life.]]
<br>
 
{{PH Model}}

Latest revision as of 22:44, 15 August 2012