MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 828407 10.1155/2013/828407 828407 Research Article Reliability Model of Mechanical Components with Dependent Failure Modes Wang Shijie Efe M. Onder 1 School of Mechanical Engineering Shenyang University of Technology Shenyang Liaoning 110870 China sut.edu.cn 2013 6 2 2013 2013 27 07 2012 15 12 2012 17 12 2012 2013 Copyright © 2013 Shijie Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

A mechanical component may fail in many modes that are usually not independent. There is generally not a joint probability density function to describe these correlated failure modes. Thus, it is difficult to compute the reliability when considering the correlations between the failure modes. It is supposed that three or more failure modes arise synchronously to be a very small probability event. The relationship between ultimate state functions in different failure modes is established by utilizing linear regression method. A double integration model for reliability of mechanical components with dependent failure modes is built according to stress-strength interference model. In case of square, cube, or exponential relationship between two ultimate state functions, a linear transformation is made. An example of pin that may fail in shear fracture, bruise, or both is discussed. The reliability is compared with that obtained by using Monte Carlo method, which represents that the reliability model with dependent failure modes proposed is correct.

1. Introduction

Mechanical components may fail in many modes such as fatigue, creep, wear, resonance, plastic deformation, and instability which lead to several ultimate state functions. The ultimate state functions are often functions of the same random variables; therefore, they are not independent, and there may be some complex or even unknown relations between them. There is generally not a joint probability density function to describe these correlated failure modes; thus, it is a difficult and hot issue to compute the reliability with dependent failure modes in reliability study . In , each failure mode is assumed to be independent when calculating the reliability with many failure modes; then, the series model of system reliability theory is introduced. This method is thought to be less rigorous in theory because the relevance of ultimate state functions among different failure modes is ignored, and the calculation error caused by it is hard to estimate. According to the relevance between the two failure modes, the bound method of reliability calculation is proposed in , which distinguishes the two failure modes into primary failure mode and secondary failure mode; the two failure modes are assumed to be of complete correlation and complete dissociation, and the actual condition of the structure falls in between the two modes, and thus the interval estimation value of relevant reliability is obtained. This method is easy to compute and widely used in the reliability design of civilian industry product [9, 10]; however, it belongs to qualitative analytical method, and only the interval estimation value can be acquired. In [11, 12], the correlation coefficient is applied to characterize the relationship among different failure modes as the nuclear power stations do for the reliability design of critical parts in many failure modes; if they are positively correlated and the correlation is desirable, a satisfying result can be obtained; if they are of negative correlation or the correlation is poor, larger error is obtained. Actually, the correlation coefficient can just illustrate the two failure modes to the extent of linear and cannot represent the real relation between them [13, 14]; for example, X2+Y2=1, where X and Y are strongly correlated, but the correlation coefficient between them is 0. Errors will be introduced when using correlation coefficient to deal with the problems of reliability calculation in many failure modes. In [15, 16], the Copula function is used to describe the correlation of the random variables, and the Copula-Reliability method with dependent failure modes is obtained in which the multiple integral operation is replaced by multiple differential operation to simplify the reliability calculation process; but the approximate correlated relation between each failure modes must be identified firstly to reasonably choose Copula model during the implementation process of this method. Moreover, more sample values and some parameter estimation skills are needed to confirm the relevant parameter. The Monte Carlo method is used to calculate the reliability model of components and systems with dependent failure modes in [17, 18]. It needs no assumption or simplicity of the relation among failure modes, basic random variables, and ultimate state functions, and thus it has a wide application scope; however, when the reliability is high or the number of random variables is large, it needs higher sampling amount. Therefore, the Monte Carlo method is usually used to verify the theory method .

It is supposed in this paper that three or more failure modes existing simultaneously is a minimum probability event, the functional relation between ultimate state functions in different failure modes is obtained through linear regression method, double integration model for reliability is established by imitating the stress-strength interference model, and the Monte Carlo method is used to verify the new reliability model.

2. Reliability Model with Multiple Failure Modes

Suppose that a mechanical component is controlled by the vector of random variables X, (1)X=[X1,X2,,XH].Xh (h=1,2,,H) are the basic random variables. Assume that there are m probable failure modes of the mechanical components, and the relative ultimate state functions are (2)g1(X)=0,g2(X)=0,,gm(X)=0, where gi(X) (i=1,2,,m) are supposed to obey normal distribution for simplicity. Take the event that the ith failure mode occurs (gi(X)<0) as Ei, and the failure probability F of the components can be presented as (3)F=i=1mP(Ei)-1i<jmP(EiEj)++(-1)m-1P(E1E2Em).

Considering that three or more failure modes occur synchronously is an event of very small probability, (3) can be expressed as (4)Fi=1mP(Ei)-1i<jmP(EiEj).

Failure probability for single failure mode of the components can be obtained as follows: (5)P(Ei)=-0fi(gi)dgi, where fi(gi) are the ultimate state functions, and gi are the probability density functions.

There are some inevitable functional relationships between generalized strength and generalized stress for the same components in each failure mode; therefore, the probability that Ei and Ej occur simultaneously cannot be calculated directly. Suppose that the relationships of ultimate state functions of each failure mode are approximately linear, as is shown in Figure 1. Even if there are other functional relationships between them, it can be translated into linear relationship through linear conversion.

Relationship of ultimate state functions of different failure modes.

The linear function can be presented as (6)gi=Ai-jgj+Bi-j+εi-j, where Ai-j, Bi-j are constants; εi-j obey normal distribution N(0,σi-j), and the intercept of the line on gj axis is recorded as ai-j.

In order to obtain the distribution function of ai-j, linear regression analysis is introduced on the N sets of observed value (gjn,gin)(n=1,2,,N) of gj, gi, and the estimation value of the unknown parameter in the straight line equation is acquired. Consider (7)A^i-j=n=1Ngjngin-(1/N)(n=1Ngjnn=1Ngin)n=1N(gjn)2-(1/N)(n=1Ngjn)2,B^i-j=n=1Ngjnn=1N(gin)2-n=1Ngjnginn=1NgjnN[n=1N(gjn)2-(1/N)(n=1Ngjn)2],σ^i-j2=1Nn=1N(gin-1Nn=1Ngin)2-A^i-j2Nn=1N(gjn-1Nn=1Ngjn)2.

Notice that the intercept ai-j of the fitting line on gj axis is (8)a^i-j=-B^i-j+εi-jA^i-j. It is thus clear that a^i-j obey normal distribution N(-B^i-j/A^i-j,σ^i-j/A^i-j).

It can be concluded from Figure 1 that

if Ai-j>0; (9)P(EiEj)=P[(gj<ai-j0)(gj<0ai-j>0)]=P(gj<ai-j0)+P(gj<0ai-j>0);

if Ai-j<0; (10)P(EiEj)=P(ai-j<gj0).

The distribution curves of gj, ai-j are presented in Figure 2, in which fi-j(ai-j) and fi(gi) denote the probability density function of gj and ai-j, respectively.

Relationship of gj and ai-j distribution functions.

It can be concluded from Figure 2 that (11)P(gj<ai-j<0)=-0fi-j(ai-j)[-ai-jfi(gj)dgj]dai-j,(12)P(gj<0ai-j>0)=-0fi(gj)dgj0+fi-j(ai-j)dai-j,(13)P(ai-j<gj0)=-0fi(gj)[-gifi-j(ai-j)dai-j]dgj. Equations (11)–(13) can be solved through numerical integration method.

3. Validation of Reliability Model with Dependent Failure Modes Utilizing Monte Carlo Method

The Monte Carlo method is also called stochastic simulation . It is a numerical computation method that is on the basis of probability statistics theory, in terms of the laws of large numbers, with the utilization of numerical simulation technology and through statistical sampling test or stochastic simulation of relevant random variable to seek for approximate solution of complex engineering and technology problems. The basic idea of the method is first to build the probability model of the problem to be solved, then produce a large number of sample values of random variable according to the statistical property of random variable of probability model, and last through sampling calculation of probability model to get the approximate solution of the problem to be solved. Mathematical difficulties in structure reliability analysis can be avoided in Monte Carlo simulation method, and no matter if the state functions are nonlinear or not, the random variables are normal distribution or not, and the limit state functions are correlated or not, as long as the number of simulation is large enough, accurate failure probability and reliability index can be obtained. Thus, the Monte Carlo simulation method is used to assess the correctness of the new method in this paper.

The hypothesis that three or more failure modes exist simultaneously is a minimum probability event which is unnecessary when using the Monte Carlo method to calculate the reliability with dependent failure modes, and sampling calculation of (14) can be directly used. Consider (14)R=P(i=1mE-i).

Concrete steps are listed as follows:

initialize loop control variable ii, sampling amount r and safe times Rg, in which r is a larger value and the others are zero;

generate N dimensions vector of random number X(ii) according to the basic random variables Xh (h=1,2,,N);

calculate g1(ii),g2(ii),,gm(ii);

test if i=1mgi(ii)>0 (i=1,2,,m) is satisfied, and if it is the case, (15)Rg(ii+1)=Rg(ii)+1;

test if iir is satisfied, and if it is the case, the reliability is (16)R=Rgr. If it is not case, ii=ii+1, back to the second step;

output reliability R.

The higher is the reliability, the more samples are required in the Monte Carlo method, and the Monte Carlo method is only used to verify the correctness of the model for reliability of components with dependent failure modes mentioned earlier.

4. Examples

Compute the reliability of the pin connection component in Figure 3. Given that Fa~N(6000,15002)N, the diameter of pin d=10 mm, h1=h3=30 mm, and h2=8 mm. The shear strength [τ] and bearing strength of 20 pins are obtained by experiment and listed in Table 1.

The Statistics Data of Shear Strength and Bearing Strength of 20 Pins.

Serial number 1 2 3 4 5 6 7 8 9 10
[ τ ] /MPa 60.36 61.84 56.73 58.71 59.27 56.60 56.31 62.61 60.03 58.20
[ σ j y ] /MPa 142.40 143.30 138.45 146.37 139.34 132.82 134.63 150.87 159.65 154.58
Serial number 11 12 13 14 15 16 17 18 19 20
[ τ ] /MPa 62.26 60.65 57.23 63.75 60.81 64.14 63.19 58.08 56.38 59.80
[ σ j y ] /MPa 167.03 144.04 139.4 148.12 144.02 149.54 149.06 145.31 133.34 137.20

A pin connection component structure.

The solution of the problem is listed as follows.

Force analysis. The force analysis of the pin is illustrated in Figure 4(a). The pin is divided into three parts I, II, and III with assumed plane m-m and plane n-n. Because the force condition of part I is the same as part III, only part I and II are taken as research objects as shown in Figure 4(b). From equilibrium conditions, we can get Q1=Q2=Fa/2. Obviously, the bearing force of part II is bigger and the bearing area is smaller. Therefore, only part II is possible to fail in the example.

Failure mode analysis. The structure has two single failure modes. They are, respectively, E1: shear fracture of the pin and E2: bruise of the pin. Moreover, the structure is possible to have two failure modes synchronously. Therefore, the failure probability F is (17)F=P(E1)+P(E2)+P(E1E2).

The limit state functions of the two failure modes are (18)g1=[τ]-2Faπd2=[τ]-0.0064Fa,g2=[σjy]-Fadh2=[σjy]-0.0125Fa. The shear strength of the pin [τ]~N(60.5,32) MPa and the bearing strength of the pin [σjy]1~N(162.9,6.12) MPa are both obtained by statistics data. The relationship between [τ] and [σjy]1 is shown in Figure 4, and the relationships between the two limit state functions are (19)g1~N(21.93,9.90);g2~N(69.42,20.72).

According to the distribution pattern of the basic variables, the accounted functional relation between g1 and g2 is shown in Figure 5. The unknown parameters are obtained from (7) as (20)A^2-1=1.90;B^2-1=26.36;σ2-12=75.45, and intercept ai-j is obtained from (8) as (21)a1-2~N(-13.87,4.57).

Calculation of failure probability F and reliability R. According to the standard normal distribution list, the failure probabilities P(E1) and P(E2) are (22)P(E1)=Φ(0-21.939.90)=1.39×10-2,P(E2)=Φ(0-69.4220.72)=5.38×10-4. The failure probability P(E1E2) can be obtained from (9), (11), and (12) as (23)P(E1E2)=-014.572πexp[-12(a2-1+13.874.57)2]{-a2-119.902πexp[-12(g1-21.939.90)2]dg1}da2-1+-019.902πexp[-12(g1-21.939.90)2]dg10+14.572πexp[-12(a2-1+13.874.57)2]da2-1=3.04×10-4+1.61×10-5=3.20×10-4,F=P(E1)+P(E2)-P(E1E2)=1.41×10-2,R=1-F=0.9859.

According to (4), conducting 105 samples with the method of Monte Carlo, the reliability R=0.985613; when the sampling amount is 106, the reliability R=0.98618, and when the sampling amount is 1.5×107, the reliability R=0.985943.

Force analysis of the pin.

Functional relation between g1 and g2.

It is clear that the reliability of the pin obtained from the theoretical approach and Monte Carlo simulation shows good agreements; therefore, the double integration model for reliability of mechanical components with dependent failure modes established in this paper is correct.

5. Conclusion

It is supposed that three or more failure modes arise synchronously to be a minimum probability event. The relationship between different ultimate state functions with dependent failure modes is established with the method of linear regression; the double integration model for reliability of mechanical components with dependent failure modes is established by imitating the stress-strength interference model, and it is verified with the method of Monte Carlo.

For the simplicity and ease of description of the problem, suppose that the relationship between ultimate state functions is linear. In case of square, cubic, or exponential relationship, linear transformation is used before utilizing this model.

Acknowledgments

This project is supported by the Ministry of Major Science and Technology of Liaoning (Grant no. 2012219020) and Shenyang Key Scientific and Technological Program (Grant no. F12-006-2-00).

Zhili S. Naisu Z. Reliability calculation of a machine component failed in different modes Journal of Northeastern University 2004 35 3 5 8 Zhang Y. Wang S. Liu Q. Wen B. Reliability analysis of multi-degree-of-freedom nonlinear random structure vibration systems with correlation failure modes Science in China E 2003 46 5 498 508 2-s2.0-22744459836 10.1360/02ye0392 Zhang Y. Zhang X. Huang X. Robust reliability-based design optimization of mechanical components with multiple failure modes China Mechanical Engineering 2009 20 2 142 149 2-s2.0-61349118601 Frost D. F. Poole K. F. A method for predicting VLSI-device reliability using series models for failure mechanisms IEEE Transactions on Reliability 1987 36 2 234 242 2-s2.0-0023365853 10.1109/TR.1987.5222353 Chung W. K. Reliability analysis of a series repairable system with multiple failures Microelectronics Reliability 1991 31 2-3 371 373 2-s2.0-0025757404 10.1016/0026-2714(91)90223-T Kodama M. Sawa I. Reliability and maintainability of a multicomponent series-parallel system with simultaneous failure and repair priorities Microelectronics Reliability 1984 24 1 147 164 2-s2.0-0021141990 10.1016/0026-2714(84)90644-9 Panko D. Vanhoucke T. Campos R. Hurkx G. A. M. Time-to-fail extraction model for the “mixed-mode” reliability of high-performance SiGE bipolar transistors Proceedings of the 44th Annual IEEE International Reliability Physics Symposium (IRPS '06) March 2006 512 515 2-s2.0-34250695472 10.1109/RELPHY.2006.251271 Ditlevsen O. Narrow reliability bounds for structural system Journal of Structural Mechanics 1979 7 4 453 472 2-s2.0-0018767337 10.1080/03601217908905329 Kozine I. Krymsky V. An interval-valued reliability model with bounded failure rates International Journal of General Systems 2012 41 8 760 773 10.1080/03081079.2012.721201 Boushaba M. Azouz Z. Reliability bounds of a 3-dimensional consecutive-k-out-of-n : F system International Journal of Reliability, Quality and Safety Engineering 2011 18 1 51 59 2-s2.0-79959777482 10.1142/S0218539311004081 Zhou Y. Yimin Z. Reliability sensitivity-based correlation coefficient calculation in structural reliability analysis Chinese Journal of Mechanical Engineering 2012 25 3 608 614 10.3901/CJME.2012.03.608 Xianzhen H. Probabilistic approach to system reliability of mechanism with correlated failure models Mathematical Problems in Engineering 2012 2012 11 465853 10.1155/2012/465853 Steiner F. Hajagos B. Reliability of the estimations of the correlation coefficient Acta Geodaetica et Geophysica Hungarica 2001 36 4 399 416 2-s2.0-0035714711 10.1556/AGeod.36.2001.4.3 Radim B. Martorell G. Reliability, Risk and Safety 2009 Boca Raton, Fla, USA CRC Press Singh A. Mourelatos Z. P. Time-dependent reliability estimation for dynamic systems using a random process approach SAE International Journal of Materials and Manufacturing 2010 3 1 339 355 2-s2.0-78649877654 Sun Y. Wang P. Cheng L. Liu H. Operational reliability assessment of power systems considering condition-dependent failure rate IET Generation, Transmission and Distribution 2010 4 1 60 72 2-s2.0-77953346144 10.1049/iet-gtd.2009.0006 Kuczera R. C. Mourelatos Z. P. Latcha M. A Monte Carlo reliability assessment for multiple failure region problems using approximate metamodels Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference (IDETC '07) September 2007 1135 1147 2-s2.0-44949204557 Solver T. Amelin M. State duration based Monte Carlo Simulation model with independent failures for distribution system reliability analysis Proceedings of the 9th International Conference on Probabilistic Methods Applied to Power Systems (PMAPS '06) June 2006 1 8 2-s2.0-46149108129 10.1049/iet-gtd.2009.0423 David P. Binder K. A Guide to Monte Carlo Simulations in Statistical Physics 2010 New York, NY, USA Cambridge University Press