Life estimation of the beam with normal distribution parameters and subjected to cyclic load

This research aims at estimating the life of the beam with normal distribution parameters and subjected to cyclic load. It is tested by Monte-Carlo simulation that the generalized displacement and velocity are normally distributed when the coefficient of variation ( ) of the random parameter is small (generally ≤ 0.01). The random perturbation method is employed to estimate the mean and variance of the generalized displacement and velocity. The random dynamic stress and its derivative with respect to the time of the beam is formulated according to the shape function of beam element and the stress equation in a Euler-Bernoulli beam. Their mean, variance and correlation coefficient are given using the first-order approximation in a Taylor series. Based on Palmgren-Miner rule, the expected cumulative damage equation is given and is used to estimate the life where the random dynamic stress is non-stationary and follows the normal distribution with the nonzero mean at any time . The presented method could also estimate the life of other structure or component which has several normal random parameters, is subjected to cyclic load and obeys the linear dynamics and elastic theory when the random parameter’s is small.


Introduction
The beam is the common structural element of the civil and mechanical engineering.It could withstand the bending moment and shear force primarily.If the beam failure occurs during service, it should lead to the huge economic lost or the casualties.The failure due to fatigue might be the common one.In mechanical engineering, there are a lot of factors which result in the beam with random parameters.The factors should include manufacturing error of machine tool, ambient temperature, random vibration of the foundation of machine tool, material uniformity and so on.Moreover, many beams used in mechanical engineering commonly carry the cyclic load, such as the loads due to all kinds of faults of the rotating machine.Therefore, it is important to estimate the life of the beam with random parameters and subjected to the cyclic load.
The random parameters of the beam result in the random stress under the deterministic cyclic load.The known methods of estimating the life of metallic structures and components under random stress can be classified into two groups: time-domain method based on a cycle counting method and frequency-domain method using power spectral density (PSD) of stochastic processes [1].One of time-domain method, the combination of rainflow counting and Palmgren-Miner rule, is recognized as the "gold standard" for estimating fatigue damage from variable amplitude stress or for comparison to other estimate techniques [2,3].However, it has been proved to spend a lot of time in processing data when time-domain method is employed to estimate the life under random stress [4,5].For the reason, frequency-domain method has been developed.It can produce (exact or approximated) analytical expressions by which we can estimate cycle distribution and fatigue life under a given counting procedure, without, however, requiring knowledge of each stress or strain time history [6,7].The life estimation of the beam under random stress has reported by some documents.The focus of the reported works is to predict the life of the beam with the deterministic parameters and withstanding the random loading.For example, the fatigue life prediction method of the dual cantilever beam subjected to the narrow band random load was introduced by Rajnikant Bhogilal Thakkar [8] and was based on the fatigue life under the variable sinusoidal load.The fatigue life of the Aluminum 6061 T6 cantilever beam, where the input was the white noise acceleration, were predicted using seven types of the frequency methods and Rainflow counting method respectively and they were compared with these obtained by experiments in [9].The effect of the damping ratios of the test structure on the fatigue life was also analyzed by the four types of tests with the different damping ratios.An analytic methodology was developed by Paulus et al. [10] to predict the fatigue life of the cantilever beam experiencing stationary and Gaussian random vibration excitation.In addition, Paulus et al. [11] proposed a semi-empirical life method of the cantilevered beam subjected to the random excitation with a wide range of PSD (power spectral density) shapes and amplitudes.Two methods were employed by Lin et al. [12,13] to predict the fatigue lives of RC beams strengthened with CFL (carbon fiber laminate) subjected to the exciting force due to the random traffic flow.They were based on the residual life model and the residual strength model respectively.Moreover, the fatigue life of the beam with the uncertain parameters was discussed by a few works.For example, the effect of the uncertainties on the fatigue life of a higher order sandwich beam was discussed by Wang and Sun [14].In their work, the uncertainties included the density of foam core materials and the thicknesses of the core and the top bottom skins, two independent white noise loadings with the same spectral density function were applied to two different points on the top skin.
It has not been reported to estimate the fatigue life of the beam which has several random parameters and is subjected to cyclic load.Therefore, it should be discussed in this work.Here, these random parameters are assumed to be normally distributed.The remainder of this paper is organized as follows: In Section 2, it should be tested by Monte-Carlo simulation that the generalized displacement and velocity follow normal distributions with the periodic mean and variance estimated by the random perturbation method, and the mean, variance and correlation coefficient of the random dynamic stress and its time derivative are formulated using the first-order approximation in a Taylor series.In Section 3, the expected cumulative damage based on Palmgren-Miner rule is modeled where the random dynamic stress is non-stationary and follows the normal distribution with the nonzero mean at any time .The conclusions are drawn finally.

Random dynamical response analysis
When the beam undergoes small random vibrations around a position of stable equilibrium and is assumed to be with viscous damping, the dynamical equation of the beam with random parameters is: where , and are the × dimension mass matrix, damping matrix and stiffness matrix respectively and ( ) is the dimension external load vector of which the element is the cyclic external load., and are the dimension generalized acceleration, velocity and displacement vector respectively.is the total number of the degree of freedom.One or multiple elements of , and are the random variables.All random parameters of the considered beam are written by one random vector .Generally, is composed of the total length , width , height ℎ, modulus of elasticity and density which are supposed to be normally distributed and mutually independent.The random perturbation method [15,16] is employed to analyze the random dynamical response of the dynamical Eq. ( 1).The mean of the generalized displacement vector is: where ( ) is the deterministic part of the generalized displacement vector and is solved by: where , and are the deterministic part of , and respectively.They are obtained by substituting the mean of all random parameters into , and .The covariance matrix of the generalized displacement vector is: Var (1) s y m m e t r y Cov ( 2), ( 1) where is the sensitivity matrix of ( ) to the random vector and is: Cov( ) is the covariance matrix of the random parameter vector and is: where is the total number of random parameters and , ( ) = ( ) ⁄ which is solved by: where is the deterministic part of .The skewness and kurtosis of the generalized displacement and velocity of node 1, 3 and 5 of the random cantilever beam at 0.005 and 0.012 second are shown in table 1.All random parameters' is the identical number in one Monte-Carlo simulation.The samples of the generalized displacement are generated by Monte-Carlo method and Newmark-method.The sample size is 5000.The considered cantilever beam with normal distribution parameters is shown in Fig. 1(b).It is divided into 5 elements.The means of the total length , width , height ℎ, density and elastic modulus are 0.25 m, 0.02 m, 0.03 m, 7850 kg/m 3 and 2.06×10 11 Pa respectively.The deterministic cyclic load ( ) = 3500sin(100 ) is applied to the free end node.From Table 1, it could be trusted that the random generalized displacement and velocity at any time are the normal distribution when all CVs of the random parameters are small (generally ≤ 0.01).
The displacement mean and standard deviation of node 1 of the cantilever beam with normal distribution parameters calculated by Monte-Carlo simulation and random perturbation method respectively with = 0.01 are shown in Fig. 2 and 3. Here, the values of the parameters except are identical with the numerical calculation of the skewness and kurtosis of the generalized displacement and velocity.The difference between Monte-Carlo simulation and random perturbation method is very little according to Fig. 2 and 3.In fact, the similar results are found in the mean and standard deviation of the generalized displacement and velocity of other nodes.Therefore, it might be one solution to estimating the mean and variance by random perturbation method when the random generalized displacement and velocity of the beam with normal distribution parameters are normally distributed at any time.

Random dynamical stress analysis
The coordinate system of the considered beam is shown in Fig. 1(a).The normal stress of a cross-section of the Euler-Bernoulli beam for the application of bending moment and transverse force is [17]: where is the axial direction coordinate of the cross-section, ( , , ) is the bending moment of the cross-section at the time , ( , ) is area moment of inertia and is the distance to neutral axis of the cross-section.By substituting the approximate differential equation of the deflection curve of the beam into Eq.( 8), the normal stress is: where , ( , , ) are the modulus of elasticity and the deflection of the cross-section at the time respectively.The maximal normal stress of the cross-section at the time is: In the finite element method for the Euler-Bernoulli beam, d ( , , ) d ⁄ between any two adjacent node and ( , = 1, 2, 3,…, < and ≤ 2 ⁄ ) is: where the -coordinate of the node is more than that of the node , , , and are the lateral displacement and rotation of the node and respectively and is the element length.According to Eq. ( 11), the maximal normal stress of the cross-section at the time is normally distributed if , , and are normally distributed.Then, by expanding the maximal normal stress of the cross-section to first-order approximation in a Taylor series at a point = , the mean and variance are respectively calculated by [18]: where: is × 1 vector and: According to the random perturbation method and Eq. ( 10) and ( 11), the partial derivative of the maximal normal stress of the cross-section with respect to the random variable ( = 1, 2,…, ) between any two adjacent node and : where 2) and ⁄ = , (2 ).By combining Eqs. ( 10) and ( 11), the time derivative of the maximal normal stress of the cross-section is: Then, its approximate mean and variance could be calculated by Eqs. ( . According to Eq. ( 16), the time derivative of the maximal normal stress of the cross-section at the time is normally distributed if , , and are normally distributed.By expanding the maximal normal stress of the cross-section and its time derivative to first-order approximation in a Taylor series at a point = , ( , , ) − ( , , ) and ( , , ) − ( , , ) respectively are: where , and , are the mean of the random variables and respectively.Then, the correlation coefficient between ( , , ) and ( , , ) is: The mean and standard deviation of the maximal fixed-end stress of the cantilever beam calculated using the presented method are compared with these by Monte-Carlo simulation in Fig. 4 and 5 respectively.The correlation coefficient between the maximal fixed-end stress and its time derivative calculated by Monte-Carlo simulation and Eq. ( 19) is shown in Fig. 6.Here, = 0.005 and other parameters are the same as Section 2.1.It could be seen that the mean, standard deviation and correlation coefficient estimated by the presented method are very close to these calculated by Monte-Carlo simulation from Figs. 4, 5 and 6.In fact, there is only a small difference between the mean, standard deviation and correlation coefficient estimated by the presented method and Monte-Carlo simulation when of all random parameters is less than about 0.01.

Life estimation of a beam with normal distribution parameters and subjected to cyclic load
According to Palmgren-Miner rule, the expected cumulative damage during the time interval to is [19]: where C and are constants depending on the material, ( , , ) and ( , , ) are replaced by and respectively for simplicity, ̅ ( ) is the expected number of crossings (at level , with positive slope) and is calculated by [20]: If and are normally distributed, their joint probability density function is: where = ( , ) is estimated by Eq. ( 19), = ( , , ) and = ( , , ).Then, the expected number of crossings is: where Erfc( ) = 1 − erf( ), erf(z) is the error function defined by:  Then, the partial derivative of the expected number of crossings with respect to is: where: Then, the expected cumulative damage of the beam with normal distribution parameters and subjected to cyclic load during the time interval t 1 to t 2 could be calculated by substituting Eq. ( 25) and ( 26) into Eq.( 20) and using the numerical integration method.Using the proposed method, the expected instantaneous damages ( ∂ ̅ ( ) ∂ ⁄ ) of the cantilever beam with normal distribution parameters are shown in Fig. 7 and 8 when = 0.005 and 0.010 respectively.Other parameters of the considered beam are the same as Section 2.1.The expected instantaneous damage of the discussed random cantilever beam subjected to cyclic load is the periodic function of the time and its amplitude with = 0.010 is more than one with = 0.005 from Fig. 7 and 8.The expected cumulative damage with = 0.005 and 0.010 during one second are 4.945×10 -4 and 1.742×10 -3 and their corresponding lives (numbers of cycles to failure) are 1.072×10 5 and 3.043×10 4 respectively.Li Changyou offers the idea and the most of the theoretical innovation of this article.The random dynamical response was analyzed and the most of the revision works of the manuscript were done by Wang Xuchu.The random dynamical stress was analyzed by Wang Wei.The prominent assistance in the theory analysis was offered by Zhang Yimin.The simulation of Section 3 was done by Guo Song.

Conclusions
One method for estimating the life of the beam with normal distribution parameters and subjected to cyclic load was presented in this work.The random generalized displacement and velocity of the beam with random parameters at any time are the normal distributions when CVs of all random parameters are small (generally ≤ 0.01).Their means and variances could be accurately estimated by the random perturbation method.By compared with Monte-Carlo simulation, the means, variances and their correlation coefficient of the random dynamic stress and its time derivative could be estimated using the first-order approximation method in a Taylor series.The presented method could estimate the life of the beam with normal distribution parameters and subjected to cyclic load where the random dynamic stress is non-stationary and follows the normal distribution with the nonzero mean at any time .

Fig. 2 .Fig. 3 .
Fig. 2. Displacement mean of node 1 of the considered cantilever beam with = 0.01 2028.LIFE ESTIMATION OF THE BEAM WITH NORMAL DISTRIBUTION PARAMETERS AND SUBJECTED TO CYCLIC LOAD.CHANGYOU LI, XUCHU WANG, WEI WANG, YIMIN ZHANG, SONG GUO

Fig. 4 . 1 Fig. 5 . 1 Fig. 6 .
Fig. 4. Mean of the maximal fixed-end stress of the considered cantilever beam when = 0.005 and other parameters are the same as Section 2.1 Fig. 5. Standard deviation of the maximal fixed-end stress of the considered cantilever beam when = 0.005 and other parameters are the same as Section 2.1

Fig. 7 . 1 Fig. 8 .
Fig. 7. Expected instantaneous damage of the cantilever beam with normal distribution parameters and subjected to cyclic load when = 0.005 and other parameters are the same as Section 2.1 Fig. 8. Expected instantaneous damage of the cantilever beam with normal distribution parameters and subjected to cyclic load when = 0.01 and other parameters are the same as Section 2.1
2028.LIFE ESTIMATION OF THE BEAM WITH NORMAL DISTRIBUTION PARAMETERS AND SUBJECTED TO CYCLIC LOAD.CHANGYOU LI, XUCHU WANG, WEI WANG, YIMIN ZHANG, SONG GUO 2028.LIFE ESTIMATION OF THE BEAM WITH NORMAL DISTRIBUTION PARAMETERS AND SUBJECTED TO CYCLIC LOAD.CHANGYOU LI, XUCHU WANG, WEI WANG, YIMIN ZHANG, SONG GUO