# A Novel Probabilistic Fatigue Life Prediction Method for Welded Structures Based on gPC

Authors: Huiying Gao; Xiaoqiang Zhang; Xiaoqiang Yang; Bo Zheng
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The traditional fatigue life prediction methods based on the S-N curve all believe that the parameters in the model are deterministic constants and can be categorized to the deterministic life prediction. However, in practice, it is difficult to carry out a large number of experiments due to the limitation of time or the possible shortage of funds. In addition, the specimens used in the experiments are not exactly the same, and the test operations and data reading depend on the accuracy of the test equipment as well as the subjective judgment of the testers, which result to the uncertainty of the S-N curve. Therefore, the uncertainty should be considered in order to improve the accuracy of the fatigue life prediction. In this paper, the uncertain factors affecting the fatigue life of welded joints are summarized, and the generalized polynomial chaos (gPC) is introduced into fatigue life prediction. A novel probabilistic fatigue life prediction method combined with the nonlinear cumulative damage model considering the uncertainty of the S-N curve is constructed. An illustrative example is presented to demonstrate the advantages of the proposed approach.


---

## Body

## 1. Introduction

The S-N curve is mainly obtained through a large number of fatigue tests and fitting analysis. Many fatigue life prediction methods based on the S-N curve believe that the parameters in the life prediction models are deterministic, that is, in a given loading environment, each parameter is a constant. Since this kind of methods considers the parameters in the models as constants, it is known as deterministic life prediction methods. In practice, these assumptions are not always true. Due to the influence of uncertain factors such as test samples, conditions, operations, and data reading, the fatigue life under constant amplitude load is often dispersed in different degrees, which leads to the uncertainty of the S-N curve. For welded joints, the S-N curve is affected not only by external uncertainties such as load and environment but also by internal factors such as the welding process, joint form, residual stress, and welding defects. The internal factors, such as welding defects, are the important causes of stress concentration. Under a cyclic loading, stress concentration usually becomes the birthplace of fatigue cracks, and it will affect the fatigue life of the welded structure. Generally, the shape, size, and stress distribution of the welded structure are uncertain, resulting in the uncertainty and dispersion of fatigue life. Even for the same type of welded joints and the same welding materials, different S-N curves will be obtained while different welding methods are used. Xu and Zhang [1] studied the fatigue properties of a TC4 thin plate of T-welded joint and found that the weld toe or weld surface was the most frequent location of fracture initiation and different welding processes had different characteristics. As a result, the welded joints with the same welding materials and joints have different S-N curves. Jonsson et al. [2] studied the fatigue performance of the welded bogie joint under different welding process parameters, that is, by changing the current and voltage, it was found that the static load strength and the shape of the joint were not sensitive to the change of welding parameters, but the fatigue limit decreased with the increase of current. There are many kinds of welded joints used in complex mechanical structures including butt joint, lap joint, and cross joint. The notch forms for each type of welded joints are different, and the degree of nonuniformity of stress distribution is usually different. For example, the stress distribution of the butt joint is relatively uniform while it is nonuniform for the cross joint because of the sharp transition from the weld to the base metal. For lap joints, the nonuniformity of stress distribution is reflected not only in the weld but also in the lap plate [3]. The nonuniform stress distribution will lead to stress concentration, and the stress concentration factors are quite different for each form of welded joints.Up to now, scholars have proposed many probabilistic fatigue life prediction models for various uncertainties. The main theoretical bases of these models include probability and statistics, entropy theory, Bayes theory, fuzzy theory, and rough set theory.The method of fatigue life prediction based on probability and statistics randomizes the relevant parameters in the models, including load, material parameters [4, 5], and structural size [6], and describes the objective uncertainty as a variable that obeys a random distribution. There are two ways to deal with this kind of methods. One is to randomize the fatigue life or damage under constant amplitude load. Rao et al. [7] considered the uncertainty of stress range and material properties and adopted the probabilistic S-N curve and a fatigue cumulative damage model to obtain the distribution of fatigue life. Shen et al. [8] considered the uncertainty of the loading process and the fatigue strength of the material, randomized the damage caused by a single cycle, and found the relationship between the distribution function and the stress amplitude. Liu and Mahadevan [9] regarded the fatigue life under different constant amplitude loads as random processes, and Karhunen–Loeve expansion was used to describe the variability of the corresponding random processes. Ni and Mahadevan [10] combined the probabilistic Miner rule and random stress-strain curve family to propose an energy-based probabilistic fatigue life prediction method under multiaxial variable amplitude loading. The other way is to randomize the parameters such as load, material properties and structural dimensions in the life prediction model. Based on the Kitagawa–Takahashi diagram, Liu and Mahadevan [11] proposed that the equivalent initial flaw size (EIFS) was only related to the fatigue limit and the threshold value of stress intensity factor. Maljaars et al. [12] predicted probabilistic fatigue crack growth of welded joints in civil engineering structures by randomizing the parameters in the model and considering the relative importance of random variables. Luo et al. [13] considered the effects of crack nucleation life and crack growth resistance variability and proposed a framework for fatigue life prediction based on crack nucleation and surface defects.Entropy theory is widely used in medicine, environment, management, statistics, physics, and so on. The main applications for entropy in probabilistic fatigue life prediction are maximum entropy and thermodynamic entropy. Yi et al. [14] proposed a method to determine the type of probability distribution of fatigue life and describe the dispersion of fatigue life. Guan et al. [15] proposed the overall framework of probabilistic fatigue damage prediction based on the maximum correlation entropy method and found that this method can give satisfactory prediction results. In addition, the maximum correlation entropy theory and Bayes theory are combined to provide support for health management decision making [16]. Thermodynamic entropy theory is also used to predict fatigue life. Zhu et al. [17] proposed a probabilistic fatigue life prediction model combining with thermodynamic entropy theory and Bayes theory. Naderi and Khonsari [18] used the theory of thermodynamic entropy to realize the real-time monitoring of fatigue life. Temfack and Basaran [19], Wang and Yao [20], and Kim et al. [21] also applied entropy theory to life prediction.The probabilistic fatigue life prediction method based on Bayes theory is another common method to deal with uncertainties. Guan et al. [22] adopted Bayes update to reduce the uncertainty in fatigue damage prediction and then used MCMC (Markov Chain Monte Carlo) simulation to select and update models and predict fatigue damage with mean probability. An et al. [23] introduced field failure data into prior information and adopted Bayes theory combined with MCMC simulation to obtain the parameters of fatigue life distribution. Chookah et al. [24] proposed a probabilistic failure physical model considering pitting fatigue degradation mechanism. The probability density function was used to represent the uncertainty of the parameters of the proposed model, and the Bayes theory combined with experimental data was used to evaluate the model. Zárate et al. [25] proposed a theoretical framework for updating and predicting crack length based on Bayes theory. The theoretical framework includes two main contents: the model update part, which is used to identify the probability density function of fracture mechanics parameters, and the prediction part, which is used to predict the crack length of the specimen and regards it as a function related to the number of cycles. Huang et al. [26], Mustafa et al. [27], Zhu et al. [28], and so on have been working on fatigue life prediction based on Bayes theory and have achieved many results.The probabilistic fatigue life prediction method based on fuzzy theory introduces fuzzy sets and membership functions to deal with the uncertain phenomena in fatigue life prediction. Since it was put forward, many models and methods have emerged. Muc [29] established the uncertainty model of fatigue life by using the fuzzy set method combined with preselection of material parameters. Nopiah et al. [30] identified low fatigue damage by using typical statistical characteristics based on fuzzy theory. Bhalla et al. [31] established a fuzzy probabilistic damage model to predict the residual life of bolted joints. Zhu et al. [32] proposed a fatigue life prediction model by the fuzzy set method.In addition to the abovementioned theories, there are other theories and methods that can be applied to probabilistic fatigue life prediction, such as rough set theory [33, 34] and grey theory [35].Among the abovementioned methods, the fatigue life prediction method based on probability and statistics is obtained on the basis of a large number of experiments and has high credibility. The life prediction method based on entropy theory, rough set theory, and grey theory has a short history and needs to be further discussed. Compared with the abovementioned methods, the development course of fatigue life prediction based on Bayes theory and fuzzy mathematics is longer while the maturity of them is not enough, and further research is needed.Although a lot of research works have been carried out, some theories and methods which have strong ability to deal with uncertain problems, such as polynomial chaos theory [36], evidence theory [37], and concept lattice theory [38], are rarely reported in the field of fatigue life prediction. In this paper, a novel probabilistic fatigue life prediction method combined with the nonlinear cumulative damage model based on gPC is proposed.The remainder of this paper is organized as follows. Section2 introduces the generalized polynomial chaos. Section 3 discusses the uncertainties and proposes a novel probabilistic fatigue life prediction method. Section 4 presents an illustrative example. Conclusions are finally summarized in Section 5.


## 2. The Generalized Polynomial Chaos (gPC)

The main idea of polynomial chaos theory is to approximately represent a random process by the sum of orthogonal polynomials corresponding to a specific distribution type. It was first proposed by Wiener [39] to establish a turbulence model. The polynomial chaos method belongs to the nonstatistical method. Compared with the traditional statistical method, the polynomial chaos method has the advantages of small amount of calculation and high accuracy. Because of these advantages, scholars applied this method to various engineering fields, such as finite deformation problems, transient heat transfer problems, and stochastic difference equations. At present, the polynomial chaos theory mainly includes generalized polynomial chaos (gPC), multielement generalized polynomial chaos (ME-gPC), and arbitrary polynomial chaos (aPC). Compared with gPC, ME-gPC can deal with the discontinuity of random space, and aPC can carry out polynomial chaos expansion while the uncertain parameters are distributed arbitrarily, but gPC has a longer development time and is more widely used and more mature.We suppose a model(1)Y=κz.Assuming thatY cannot be expressed analytically by z, but the probability space of z is smooth sufficiently, and the response Y corresponding to z can be obtained by numerical simulation, the relevant statistical properties of Y can be obtained by the surrogate model. For a system, if the input is affected by random parameters represented by the set ξii=1n, then, since the uncertainty of the output is only determined by the input, the output of the system can also be represented by the set ξii=1n. For the random field of Gaussian distribution, the series of Wiener–Askey polynomial chaos converges exponentially, so Wiener polynomial chaos can well deal with the problem of Gaussian distribution random field. For the random field with non-Gaussian distribution, Ghanem [40] proved that the series of Wiener polynomials chaos also converges, but the speed is relatively slow. In order to solve this problem, Xiu [41] and Sepahvand [42] extended Hermite polynomials and finally obtained a family of polynomial chaos for random variables with different distribution types. The polynomial chaos families corresponding to common distribution types are shown in Table 1.Table 1
The polynomial chaos corresponding to random variables of different distribution types.
Distribution typePolynomial chaos typeGaussianHermite chaosUniformLegendre chaosGammaLaguerre chaosBetaJacobi chaosPoissonCharlier chaosBinomialKrawtchouk chaosHypergeometricHahn chaosNonnegative binomialMeixner chaosAssuming that the input random variables of the system follow a Gaussian distribution, the polynomial chaos corresponding to the Gaussian distribution can be expanded as(2)Y=α0Γ0+∑i1∞αi1Γ1ξi1+∑i1=1∞∑i2=1i1αi1i2Γ2ξi1,ξi2+∑i1=1∞∑i2=1i1∑i3=1i2αi1i2i3Γ3ξi1,ξi2,ξi3+⋯,where αii=1,2,3,… are the undetermined constant coefficients, ξjj=1,2,3,… are independent variables and obey the standard normal distribution, and Γkξi1,ξi2,…,ξik is the kth multidimensional Hermite polynomial.For simplicity, equation (2) can be written in a compact form:(3)Y=∑j=0∞βjΠjξ,where βj is the coefficient of polynomial, Πjξ is the Wiener–Askey polynomial chaos, and ξ=ξ1,ξ2,…,ξn is the random variable vector.Assuming thatξ is a two-dimensional random variable vector, equation (2) can be rewritten as(4)Y=α0Γ0+α1Γ1ξ1+α2Γ1ξ2+α11Γ2ξ1,ξ1+α12Γ2ξ1,ξ2+α22Γ2ξ2,ξ2+α111Γ3ξ1,ξ1,ξ1+α112Γ3ξ1,ξ1,ξ2+α122Γ3ξ1,ξ2,ξ2+α222Γ3ξ2,ξ2,ξ2+⋯.Comparing equations (4) with (3), it can be seen that the function is changed only by the subscript. Therefore, equation (4) is equivalent to(5)Y=β0Π0+β1Π1+β2Π2+β3Π3+β4Π4+β5Π5+β6Π6+β7Π7+β8Π8+β9Π9+⋯.The expression of then-order Hermite polynomial is(6)Γnξi1,ξi2,ξi3,…,ξin=e1/2ξTξ−1n∂n∂ξi1,∂ξi2,…,∂ξine−1/2ξTξ.According to equation (6), it is easy to obtain a one-dimensional Hermite polynomial.(7)Γ0ξ=1,Γ1ξ=ξ,Γ2ξ=ξ2−1,Γ3ξ=ξ3−3ξ,…,Γnξ=ξΓn−1ξ−n−1Γn−2ξ.Tables2 and 3 show the one-dimensional and two-dimensional Hermite polynomial of each order, respectively.Table 2
The Hermite polynomial of each order for a one-dimensional random variable.
jThe orderχThejth Hermite polynomial Πj0χ=011χ=1ξ12χ=2ξ12−13χ=3ξ13−3ξ14χ=4ξ14−6ξ12+3Table 3
The Hermite polynomial of each order for two-dimensional random variables.
jThe orderχThejth Hermite polynomial Πj0χ=011χ=1ξ12ξ23χ=2ξ12−14ξ1ξ25ξ22−16χ=3ξ13−3ξ17ξ12ξ2−ξ28ξ1ξ22−ξ19ξ23−3ξ210χ=4ξ14−6ξ12+311ξ13ξ2−3ξ1ξ212ξ12ξ22−ξ12−ξ22+113ξ1ξ23−3ξ1ξ214ξ24−6ξ22+3If the random variables in Hermite polynomials follow the standard normal distribution and are independent from each other, then the polynomials are orthogonal, that is,(8)Πm,Πn=0,m≠n.Then, the Hermite polynomials form a set of complete orthogonal bases in 2-dimensional space, and it has(9)Πm,Πn=Πm2δmn,pξ,qξ=∫pξqξθξdξ,where ·,· is the inner product, δmn is the Kronecker function, θξ is the weight function of the Hermite polynomial, and it satisfies(10)θξ=12πe−1/2ξTξ.The coefficient of Hermite polynomialβj in equation (3) can be obtained by Galerkin projection method, that is,(11)βj=Y,ΠiΠi2=1Πi2∫YΠiθξdξ.This Galerkin projection method has some limitations; it requires thatY must have a clear analytical solution. Therefore, the random response surface method will be adopted in this paper.


## 3. The Proposed Probabilistic Fatigue Life Prediction Method

### 3.1. The Uncertainty Analysis and Probabilistic Fatigue Life Prediction of the S-N Curve

For the life prediction method based on the S-N curve and fatigue cumulative damage theory, the S-N curve is crucial, which affects the prediction accuracy. In the case of a given stress ratio, the S-N curve is tested under different levels of constant amplitude load, and the corresponding fatigue life of all levels of load is recorded and then fitted in the later stage. In the deterministic fatigue life prediction method, the “N” in the S-N curve usually refers to the average fatigue life under all levels of stress, without considering its dispersion. However, for welded structures and parts in engineering, this treatment method often cannot achieve satisfactory results and may even lead to catastrophic accidents.The expression of the S-N curve is more commonly used in the form of power function, that is, the relationship between fatigue life and stress level is expressed by(12)NfΔσm=C,where Nf is the fatigue life cycle, Δσ is the stress, and C and m are the material parameters.Taking logarithms on both sides of the abovementioned equation, it can be transformed into(13)lnNf+mlnΔσ=lnC.From equation (12), it is found that the fatigue life cycle Nf will be affected by material parameters C and m under the action of steel. The study shows that, for a given material, the parameters C and m are not fixed constants, and there is a certain degree of dispersion and uncertainty, which is due to the uncertainty of inclusions, pores, and microcracks in the specimen. For welded joints, the internal welding defects are inevitable, including more pores, inclusions, and microcracks, than nonwelded structures, and the distribution is more complex than that of nonwelded structures. Therefore, when the method based on the S-N curve and fatigue cumulative damage theory is used to predict the life of welded joints, it is necessary to consider the uncertainty of the S-N curve which originates from the uncertainty of material parameters C and m. In this paper, the gPC is proposed to consider the uncertainty of the S-N curve and combine with the improved Manson–Halford model and the improved Corten–Dolan model to predict the probabilistic fatigue life.The fatigue life under constant amplitude load is obtained on the basis of the S-N curve. Generally speaking, it is considered that the fatigue lifeNx follows a lognormal distribution or a Weibull distribution. x=x1,x2,…,xn is the uncertain parameter vector of the fatigue life. For example, Nx follows a lognormal distribution, while lnNx follows a normal distribution. In order to establish the polynomial chaos expansion model of fatigue life, it is assumed that the uncertain parameters x1,x2,…,xn are independent of each other. Therefore, the fatigue life corresponding to all levels of stress can be approximately calculated by(14)lnNx=∑j=0∞βj∏jx,x=x1,x2,…,xn.According to the properties of the polynomial chaos expansion, the finite term can be used to approximately estimate the fatigue life. Assuming that the number of terms iss, the fatigue life prediction model can be simplified as(15)lnNx=∑j=0s−1βj∏jx,x=x1,x2,…,xn,where s=n+χ!/n!χ! and χ is the order of expansion for polynomial chaos.In order to improve the speed of calculation,x1,x2,…,xn can be normalized into random variables ζ1,ζ2,…,ζn.(16)ζi=xi−ExiVarxi,where Exi and Varxi denote the expectation and variance of the random variable xi, respectively. Since the random variables are independent from each other, after standardization, the random variables in ζ=ζ1,ζ2,…,ζn are independent from each other. The multidimensional orthogonal polynomial ∏jζ can be derived as(17)∏jζ=∏i=1nφiliζi,where j=0,1,2,…,s−1, li=0,1,2,…,χ, ∑i=1nli≤χ. φiliζi is the li-order orthogonal polynomial, whose parameters are random variables ζii=1,2,⋯,n. According to the distribution type of ζi, the corresponding relations in Table 1, the base of one-dimensional orthogonal polynomial φiliζi can be determined.It is assumed that the fatigue life of the material under constant amplitude loadΔσ is expressed by Nx. It is known that the material parameters C and m are random variables and x=lnC,m. Then, the polynomial chaos of fatigue life under constant amplitude load can be obtained by substituting x=lnC,m into the abovementioned model and solving the coefficients of the polynomial chaos. The probabilistic fatigue life prediction model based on gPC is finally obtained. The whole calculation process is shown in Figure 1.Figure 1
The flow chart of probabilistic fatigue life prediction based on gPC.From Figure1, it is known that the key point for the uncertainty analysis of the S-N curve and probabilistic fatigue life prediction based on gPC is solving the coefficients of the polynomial chaos. In this paper, the stochastic response surface method is used.


### 3.2. gPC-Based Stochastic Response Surface Method

As mentioned above, how to calculate the coefficients of the polynomial chaosβj is the fundamental and critical procedure. In the space expanded by the base vector of ξii=1n, each set of determined sample value corresponds to a point named sampling point. For the model with multidimensional random variables, the common prescription is to accomplish polynomial chaos expansion and calculate the coefficients via choosing appropriate sampling points; then, the polynomial chaos of the model is obtained. The key of the stochastic response surface method is how to choose the appropriate sampling points. Taking the random field of Gaussian distribution as an example, it can be expanded into an approximate expression based on Hermite polynomials according to the Wiener polynomial chaos theory. It is assumed that the highest order in the expansion is χ, and then, the roots of the χ+1 order Hermite primary function can be taken as the collocation points.Letn=2,χ=2, that is, the model containing two-dimensional random variables is expanded to a second-order Hermite polynomial, then s=χ+n!/χ!n!=2+2!/2!2!=6, and the combinatorial number of the corresponding roots is χ+1n=2+12=9. The partial Hermite primary functions corresponding to the Gaussian distribution random field are shown in Table 4.Table 4
The Hermite primary functions.
Γ0x=1Γ1x=xΓ2x=1/2x2−1Γ3x=1/6x3−3xΓ4x=1/26x4−6x2+3Γ5x=1/230x5−10x3+15xΓ6x=1/125x6−15x4+45x2−15From Table4, the root of the third-order Hermite primary function can be obtained, that is, 0,3,−3; thus, the corresponding sampling points include 0,0, 0,3, 0,−3, 3,0, 3,3, 3,−3, −3,0, −3,3, and −3,−3. The selection of the sampling points should follow some principles, including that the sampling points should be chosen as close as possible to the origin and it should be tried to make them symmetrical about the origin. For the random field of Gaussian distribution, the collocation point should include the origin, and the number of sampling points is generally twice as the expansion terms of the polynomial chaos. Then, the unknown coefficients are obtained by the collocation approximation method, and the specific formula satisfies(18)Π0η0Π1η0…Πs−1η0Π0η1Π1η1…Πs−1η1⋮⋮⋮⋮Π0ηNΠ1ηN…Πs−1ηNβ0β1⋮βs−1=Yη0Yη1⋮YηN,where η0,η1,…,ηN are the sampling points and N is the number of the sampling points.The whole process of solving the expansion coefficients of polynomials by the random response surface method can be represented by the flow chart shown in Figure2.Figure 2
The implementation process of the stochastic response surface method.


## 3.1. The Uncertainty Analysis and Probabilistic Fatigue Life Prediction of the S-N Curve

For the life prediction method based on the S-N curve and fatigue cumulative damage theory, the S-N curve is crucial, which affects the prediction accuracy. In the case of a given stress ratio, the S-N curve is tested under different levels of constant amplitude load, and the corresponding fatigue life of all levels of load is recorded and then fitted in the later stage. In the deterministic fatigue life prediction method, the “N” in the S-N curve usually refers to the average fatigue life under all levels of stress, without considering its dispersion. However, for welded structures and parts in engineering, this treatment method often cannot achieve satisfactory results and may even lead to catastrophic accidents.The expression of the S-N curve is more commonly used in the form of power function, that is, the relationship between fatigue life and stress level is expressed by(12)NfΔσm=C,where Nf is the fatigue life cycle, Δσ is the stress, and C and m are the material parameters.Taking logarithms on both sides of the abovementioned equation, it can be transformed into(13)lnNf+mlnΔσ=lnC.From equation (12), it is found that the fatigue life cycle Nf will be affected by material parameters C and m under the action of steel. The study shows that, for a given material, the parameters C and m are not fixed constants, and there is a certain degree of dispersion and uncertainty, which is due to the uncertainty of inclusions, pores, and microcracks in the specimen. For welded joints, the internal welding defects are inevitable, including more pores, inclusions, and microcracks, than nonwelded structures, and the distribution is more complex than that of nonwelded structures. Therefore, when the method based on the S-N curve and fatigue cumulative damage theory is used to predict the life of welded joints, it is necessary to consider the uncertainty of the S-N curve which originates from the uncertainty of material parameters C and m. In this paper, the gPC is proposed to consider the uncertainty of the S-N curve and combine with the improved Manson–Halford model and the improved Corten–Dolan model to predict the probabilistic fatigue life.The fatigue life under constant amplitude load is obtained on the basis of the S-N curve. Generally speaking, it is considered that the fatigue lifeNx follows a lognormal distribution or a Weibull distribution. x=x1,x2,…,xn is the uncertain parameter vector of the fatigue life. For example, Nx follows a lognormal distribution, while lnNx follows a normal distribution. In order to establish the polynomial chaos expansion model of fatigue life, it is assumed that the uncertain parameters x1,x2,…,xn are independent of each other. Therefore, the fatigue life corresponding to all levels of stress can be approximately calculated by(14)lnNx=∑j=0∞βj∏jx,x=x1,x2,…,xn.According to the properties of the polynomial chaos expansion, the finite term can be used to approximately estimate the fatigue life. Assuming that the number of terms iss, the fatigue life prediction model can be simplified as(15)lnNx=∑j=0s−1βj∏jx,x=x1,x2,…,xn,where s=n+χ!/n!χ! and χ is the order of expansion for polynomial chaos.In order to improve the speed of calculation,x1,x2,…,xn can be normalized into random variables ζ1,ζ2,…,ζn.(16)ζi=xi−ExiVarxi,where Exi and Varxi denote the expectation and variance of the random variable xi, respectively. Since the random variables are independent from each other, after standardization, the random variables in ζ=ζ1,ζ2,…,ζn are independent from each other. The multidimensional orthogonal polynomial ∏jζ can be derived as(17)∏jζ=∏i=1nφiliζi,where j=0,1,2,…,s−1, li=0,1,2,…,χ, ∑i=1nli≤χ. φiliζi is the li-order orthogonal polynomial, whose parameters are random variables ζii=1,2,⋯,n. According to the distribution type of ζi, the corresponding relations in Table 1, the base of one-dimensional orthogonal polynomial φiliζi can be determined.It is assumed that the fatigue life of the material under constant amplitude loadΔσ is expressed by Nx. It is known that the material parameters C and m are random variables and x=lnC,m. Then, the polynomial chaos of fatigue life under constant amplitude load can be obtained by substituting x=lnC,m into the abovementioned model and solving the coefficients of the polynomial chaos. The probabilistic fatigue life prediction model based on gPC is finally obtained. The whole calculation process is shown in Figure 1.Figure 1
The flow chart of probabilistic fatigue life prediction based on gPC.From Figure1, it is known that the key point for the uncertainty analysis of the S-N curve and probabilistic fatigue life prediction based on gPC is solving the coefficients of the polynomial chaos. In this paper, the stochastic response surface method is used.


## 3.2. gPC-Based Stochastic Response Surface Method

As mentioned above, how to calculate the coefficients of the polynomial chaosβj is the fundamental and critical procedure. In the space expanded by the base vector of ξii=1n, each set of determined sample value corresponds to a point named sampling point. For the model with multidimensional random variables, the common prescription is to accomplish polynomial chaos expansion and calculate the coefficients via choosing appropriate sampling points; then, the polynomial chaos of the model is obtained. The key of the stochastic response surface method is how to choose the appropriate sampling points. Taking the random field of Gaussian distribution as an example, it can be expanded into an approximate expression based on Hermite polynomials according to the Wiener polynomial chaos theory. It is assumed that the highest order in the expansion is χ, and then, the roots of the χ+1 order Hermite primary function can be taken as the collocation points.Letn=2,χ=2, that is, the model containing two-dimensional random variables is expanded to a second-order Hermite polynomial, then s=χ+n!/χ!n!=2+2!/2!2!=6, and the combinatorial number of the corresponding roots is χ+1n=2+12=9. The partial Hermite primary functions corresponding to the Gaussian distribution random field are shown in Table 4.Table 4
The Hermite primary functions.
Γ0x=1Γ1x=xΓ2x=1/2x2−1Γ3x=1/6x3−3xΓ4x=1/26x4−6x2+3Γ5x=1/230x5−10x3+15xΓ6x=1/125x6−15x4+45x2−15From Table4, the root of the third-order Hermite primary function can be obtained, that is, 0,3,−3; thus, the corresponding sampling points include 0,0, 0,3, 0,−3, 3,0, 3,3, 3,−3, −3,0, −3,3, and −3,−3. The selection of the sampling points should follow some principles, including that the sampling points should be chosen as close as possible to the origin and it should be tried to make them symmetrical about the origin. For the random field of Gaussian distribution, the collocation point should include the origin, and the number of sampling points is generally twice as the expansion terms of the polynomial chaos. Then, the unknown coefficients are obtained by the collocation approximation method, and the specific formula satisfies(18)Π0η0Π1η0…Πs−1η0Π0η1Π1η1…Πs−1η1⋮⋮⋮⋮Π0ηNΠ1ηN…Πs−1ηNβ0β1⋮βs−1=Yη0Yη1⋮YηN,where η0,η1,…,ηN are the sampling points and N is the number of the sampling points.The whole process of solving the expansion coefficients of polynomials by the random response surface method can be represented by the flow chart shown in Figure2.Figure 2
The implementation process of the stochastic response surface method.


## 4. An Illustrative Example

To verify the feasibility of the proposed probabilistic fatigue life prediction method, two different nonlinear fatigue cumulative damage models including the improved Manson–Halford model [43] and the improved Corten–Dolan model [44] are used to predict the probabilistic fatigue life, respectively. The prediction results are compared with those of the Monte Carlo method. The probabilistic fatigue life prediction of aluminum alloy welded joints is carried out considering the uncertainty of the S-N curve. Constant amplitude loading or two-stage loading mode is usually an ideal load-bearing state, but in engineering practice, the load of welded joints is more complex and changeable. Therefore, this section will predict the probabilistic fatigue life of welded joints under multistage variable amplitude loading on the basis of gPC, the improved Manson–Halford model, and the improved Corten–Dolan model. The data in this section are from the fatigue tests of 45 # steel welded joints under constant amplitude and program block loads. A total of five groups of 10 specimens were tested under constant amplitude loading, and the test fatigue life under applied stress amplitude and all levels of load is shown in Table 5. The program block load is loaded by multistage variable amplitude load, and the amplitude and cycle times of all levels of stress are shown in Figure 3. A total of 6 specimens are used to carry out the program block load test, and the fatigue life obtained is shown in Table 6.Table 5
The fatigue life test data of a 45 # steel welded joint under constant amplitude load.
Group IDExperiment lifeNf/105 (105 cycles)Δσ = 750 MPaΔσ = 650 MPaΔσ = 630 MPaΔσ = 590 MPaΔσ = 520 MPa10.20800.49400.70891.32191.967020.22150.81750.90411.41292.095130.24300.85900.92811.52693.236740.25500.96290.94301.56784.562550.27001.02000.96161.64214.613260.29001.08691.00301.70024.738170.34701.12691.19811.78485.233680.37801.18301.20201.92805.481590.46101.19291.53212.46386.5569100.61001.68501.89192.627212.9211Figure 3
The schematic diagram of multistage variable amplitude loading for a 45 # steel welded joint.Table 6
The fatigue life test data of a 45 # steel welded butt joint under multistage variable amplitude load.
Group IDExperimental value of fatigue life (load block)Predicted value of fatigue life (load block)Modified Manson–Halford modelModified Corten–Dolan model13.7179.391813.856625.18936.72947.45257.943612.151First of all, the statistical characteristics of the parametersC and m which affect the fatigue life are analyzed according to Table 5. There are 5 groups of fatigue tests carried out in reference [45]. Taking one from the 10 test data under each stress amplitude, a group of C and m can be obtained by fitting. There are 50 experimental data here, and there are a total of C1015=105 groups. It is not realistic to combine and fit all the forms. Therefore, 200 combinations are randomly selected for statistical analysis of C and m. The statistical results are shown in Figures 4 and 5, where lnC and m follow normal distributions. Therefore, the influence of the uncertainty of the S-N curve on the fatigue life can be analyzed by the Hermite polynomial chaos expansion, as shown in Figures 6 and 7.Figure 4
The statistical distribution results forC of the S-N curve (45 # steel welded butt joint).Figure 5
The statistical distribution results form of the S-N curve (45 # steel welded butt joint).Figure 6
The probabilistic fatigue life prediction results under the modified Manson–Halford model (45 # steel welded butt joints).Figure 7
The probabilistic fatigue life prediction results under the modified Corten–Dolan model (45 # steel welded butt joints).The probability of each predicted value can be judged from Figures6 and 7, that is, when the uncertainty of the S-N curve is considered, the fatigue life predicted by the improved Manson–Halford model and the improved Corten–Dolan model has a certain uncertainty, but not a definite value. This can provide a more reasonable explanation for the abovementioned phenomenon, that is, the welded joint designed according to the model with high prediction accuracy may still have fatigue failure during its safety period.The gPC method uses the addition and subtraction algorithm to express the dispersion and randomness of fatigue life caused by the uncertainty of material parametersC and m under constant amplitude load, so as to conside the uncertainty of the S-N curve. Compared with the power function form, the calculation is faster. However, compared with the Monte Carlo method, the polynomial chaos theory has a limitation, that is, the calculation process of the polynomial chaos expansion method is more complex and tedious when there are more parameters, and the calculation efficiency is higher when the parameters are less. Therefore, in the fatigue life prediction of welded joints, if there are more uncertain parameters input, the polynomial chaos theory still needs to be studied more deeply in order to simplify the calculation process and save time and cost. In the analysis of this chapter, the uncertainty of the S-N curve is due to the dispersion and randomness of the two material parameters. In the process of calculation, it is found that the fatigue life prediction under second-order and multistage loading by the gPC method is more efficient than that by the Monte Carlo method, and the prediction results are close. Therefore, the fatigue life prediction method based on gPC expansion and considering the uncertainty of the S-N curve proposed in this chapter can obtain satisfactory probabilistic life prediction results through a simple and efficient calculation process.By comparing the probabilistic fatigue life of a 45 steel welded butt joint obtained by the gPC expansion method and Monte Carlo simulation method under multistage variable amplitude block load, a conclusion similar to that of second-order loading mode can be obtained, that is, the probabilistic fatigue life prediction result of polynomial chaos expansion is quite close to that of the Monte Carlo method. This shows that the polynomial chaos theory is applicable in the field of fatigue life prediction and is an effective method to reflect and characterize the uncertainty of the S-N curve.


## 5. Summary and Conclusions

Based on the gPC, the uncertain parameters that affect the fatigue life under constant amplitude load are regarded as random variables to consider the uncertainty of the S-N curve of welded joints. Combined with the improved Manson–Halford model and the improved Corten–Dolan model, the probabilistic fatigue life of welded joints is predicted. By comparing the results of probabilistic fatigue life prediction under secondary and multistage load, it is found that the probabilistic fatigue life prediction method based on gPC is close to the Monte Carlo method when it reaches the second or third order. Therefore, this method can be applied to fatigue life prediction under secondary or more complex loading conditions. The polynomial chaotic expansion method expresses the fatigue life as the polynomial of the addition and subtraction rule, and the calculation efficiency is higher than that of the Monte Carlo method, which verifies the feasibility of gPC expansion in fatigue life prediction, broadens the engineering application field of polynomial chaos theory, and provides a more practical way to deal with the uncertainty in the process of life prediction.



---

*Source: 5534643-2021-04-26.xml*