Vibration and Reliability Analysis of Non-Uniform Composite Beam under Random Load

: Non-uniform structures and composite materials have advantages in engineering applications, such as light weight design, multi-functionality, and better buckling/ﬂutter load capacity. For composite structures under dynamic loading conditions, reliability is a key problem to be analyzed during practical operations. However, there is little research work on non-uniform composite structural reliability analysis under random load. The forced vibration response of non-uniform composite beam under random load is ﬁrstly solved by the Adomian Decomposition Method (ADM) and iterative process for reliability analysis. Different variation laws of the cross-section rigidity and mass distribution along the length of the non-uniform composite beam structures are analyzed. Various angular frequency and amplitude of random base motion acceleration following Gaussian distribution are considered. Inﬂuences of different random excitations and structural design on vibration responses and reliability are studied. The larger mean and variance of excitation frequency leads to the smaller amplitude and strain of the beam, while greater mean and variance of the base motion excitation amplitude will induce the higher maximum amplitude and strain values and lower reliability. The inﬂuences of structural design on reliability are studied. The reliability increases with the increment of taper ratios of the host beam and composite layer. The iteration mathematical model and numerical solutions proposed in this paper can be used to solve and analyze vibration responses and reliability of general non-uniform composite beam structures under arbitrary excitation during a certain period of time.


Introduction
Non-uniform composite beams, including functionally graded materials (FGM) beams, have important applications in engineering because they can change geometric and material properties, while most of the excitations applied on mechanical structures in nature and engineering practice are random forces. In recent years, the problems of mechanical structure damage caused by vibration have achieved more and more attention, which make the vibration reliability analysis of non-uniform composite beam structures under random load of great significance.
Before the study of non-uniform composite beams, researchers have done a lot of work on uniform composite beams dynamics. In references [1][2][3][4][5][6][7][8][9], the free vibration of a composite beam structure is solved by the finite element method. In addition to the finite element method, there are some other methods to solve the free vibration responses of composite beams. By using a global higher order beam theory, Matsunaga analyzed natural frequencies and buckling stresses of simply supported multilayered composite beams through the method of power series expansion of displacement components [10]. Atlihan et al. used differential quadrature method (DQM) to analyze the free vibration of laminated composite beams [11]. The natural frequencies of longitudinal and transverse vibration modes have been obtained by using the method of power series [12]. Except for the free vibration analysis, Raja, Prathap and Sinha developed a control scheme based on the linear quadratic regulator/independent modal space control (LQR/IMSC) method and used this to do the active vibration control of composite sandwich beams [13]. Ramanamurthy, Chandrasekaran and Nishant used finite element method for damage detection [14]. Xu et al. proved the feasibility of the ultrasonic non-destructive testing with the combination of a strain measurement for fatigue crack details detection of headed shear studs in composite beams [15]. Tao et al. obtained the fitting formula for calculating the equivalent flexural stiffness of composite beams through extensive parameter analysis, and they studied and analyzed the equivalent flexural stiffness of composite frame beams to evaluate their vertical deflection [16]. Kim et al. studied the free and forced vibration of cracked laminated composite beam through the Jacobi-Ritz method and the first-order shear deformation theory (FSDT) [17,18].
In addition to research and understanding of the uniform composite beam, some researchers have also used various methods to solve and analyze vibration responses of the non-uniform beams. For the geometrically non-uniform beams, the vibration responses of non-uniform beams were solved and analyzed by the amatrix transfer method in [19][20][21][22]. In addition, Wu investigated the free vibration characteristics of a non-uniform cantilever beam carrying multiple two degree-of-freedom spring-damper-mass systems by means of two finite element methods, FEM1 and FEM2 [23]. Ho and Chen used Timoshenko beam theory, Hamilton principle and differential transformation method (DTM) to study the vibration of axially loaded non-uniform spinning twisted Timoshenko beam [24]. Martinez Castro et al. gave the semi-analytic solution of moving load problem, which is used to analyze the stress of multi-span uniform beam and non-uniform beam under a moving load [25]. Mazanoglu and Sabuncu studied the bending vibration of non-uniform Rayleigh beams with single-edge and double-edge cracks [26]. Chen et al. presented the design, test and analysis of a non-uniform thickness piezoelectric beam for impact vibration energy harvesting. [27]. Clementi et al. determined the frequency response curve of a non-uniform beam under nonlinear vibration by the multiple time scale method analytically [28]. Based on the nonlocal elastic theory, Chakraverty and Behera studied the free vibration of non-uniform Euler-Bernoulli nanobeams by using the Rayleigh-Ritz method [29]. Celik studied the free vibration of non-uniform Euler-Bernoulli beams under different support conditions through the Chebyshev wavelet collocation method [30]. For beam structures with anisotropic material properties (FGM), the finite element method was used to solve and analysis the free vibration of FGM beam in references [31][32][33]. Differential transformation method and transfer matrix method (TMM) were also used to study the free vibration of an FGM beam in references [34,35]. Abdelrahman et al. developed an analytical solution method using Navier's procedure to study the dynamic behavior of carbon nanotube reinforced functional gradient (FG) beams on two parameter elastic foundation under a moving load [36]. Esen et al. studied and analyzed the free vibration and dynamic response of an S-shaped FG Timoshenko beam model under a moving load based on Navier's method [37]. There is also some research work on the non-uniform FGM beam. Rajasekaran used the differential transformation (DT) based dynamic stiffness approach to study the free vibration of axially functionally graded non-uniform beams under different boundary conditions [38]. Heshmati and Daneshmand described the effect of different profile variations on vibrational properties of non-uniform beams made of graded porous materials [39]. Under the framework of nonlocal strain gradient theory, Rajasekaran and Khaniki used the Lagrange interpolation method, Gaussian quadrature method and Wilson multiplier method to study the bending, buckling and vibration behavior of depth functionally graded in depth direction non-uniform nanobeams [40]. Adomian Decomposition Method (ADM) has been used to analyze the free and forced vibration of general non-uniform cone beams with axial functionally graded material by Keshmiri et al. [41][42][43].
Based on dynamic modeling and analysis of different structures, their reliability analysis under random load excitation can be conducted. Siddiqui and Ahmad carried out nonlinear dynamic analysis of the platform for response calculation. The response history obtained has been used for the fatigue reliability analysis of Tension leg platform (TLP) tethers under long crested random sea and associated wind [44]. Jensen et al. used a standard gradient-based algorithm with a line search to solve the reliability optimization problem of a structural system under random load [45]. Liu and Zhang deduced the calculation expression of the mean value of structural dynamic reliability under stationary random excitation, and conveniently obtained the dynamic reliability analysis results of random structures under random excitation [46]. Singh et al. determined the fatigue reliability of an automotive crankshaft by stochastically inducing random loads [47]. In summary, there have been many studies on the dynamic analysis of composite structures and non-uniform beams, and the vibration and reliability analysis of uniform structures and mechanical parts with certain design under random loads is also a common and wellstudied problem in engineering. However, to the best of the authors' knowledge, there is no research on the vibration reliability analysis of a non-uniform composite beam under a random load. Considering the various mechanical properties as well as unknown response and strain/stress distribution of the non-uniform structures made from multiple composite layers, the reliability of such structures can be different from traditional uniform structures as novel research [48,49], especially under different random dynamic loads. There is no explicit expression for the random load whose excitation frequency and amplitude are changing, so the analytical solution of forced vibration at each time point needs to be solved by iterative method. To deal with this challenge, a model combining ADM and iteration numerical process is proposed to solve the vibration response and dynamic stressdistribution of the non-uniform composite beams for reliability analysis in this research.
In this paper, the vibration responses and reliability of a general non-uniform beam with tapered composite layers and host beam under random load are calculated and analyzed. The cross-section rigidity and mass distribution can be adjusted by changing the volume proportion of different materials of the non-uniform laminated composite beam. The angular frequency and amplitude of random base motion acceleration accord with Gaussian distribution. The vibration mode shapes and nature frequencies are calculated, and the vibration solution under random base motion excitation is solved by the iteration numerical method. In the parameter studies, the influences of the base motion excitation frequencies variation on amplitude and strain of beam are studied. Then, the effects of base motion excitation amplitudes on vibration responses and reliability are analyzed. In addition, the vibration responses and reliability in different periods of time are calculated. Finally, the reliability of non-uniform beams with different taper ratios of composite layers and host beam is solved and analyzed.

Theoretical Model
In this section, ADM is used to solve the natures (nature frequencies and mode shapes) of the non-uniform composite beam. Then, an iterative model is used to describe and solve the vibration response and reliability of the non-uniform composite beam under random load. This random excitation is caused by base motion, and its angular frequency and amplitude of acceleration conform to Gaussian distribution with specified mean and variance. Finally, the definition and calculation method of vibration reliability are given.

Natures (Nature Frequencies and Mode Shapes) of Non-Uniform Composite Beam
In this study, a sample non-uniform rectangular composite cantilever beam structure with varied cross section is analyzed and shown in Figure 1. In the middle of the composite beam is a host beam made of aluminum, and composite layers of piezoelectric ceramic material are pasted above and below the host beam. At the left fixed end, the length, width and thickness of the composite beam and the thickness of the host beam are L, b 0 , 2H 0 , and 2h 0 , respectively. The thicknesses of the host beam and composite layer at x position Appl. Sci. 2022, 12, 2700 4 of 20 are 2h 1 (x) and 2h 2 (x). x is the position variable along the beam length. The width and thickness of the composite beam at x position are b(x) and 2H(x) = 2h 1 (x) + 2h 2 (x). F(x, t) is the random excitation.
( , ) is the random excitation. The mass distribution and cross-section rigidity at position of the composite beam are ( ) ( ) = 1 1 ( ) + 2 2 ( ), where ( ) = 2 ( ) ( ), 1 , 1 ( ), 1 , 2 , 2 ( ) and 2 are mass density, cross-sectional area and modulus of elasticity of host beam and composite layer.  The equation of free vibration of a non-uniform composite beam without considering damping is given below, where ( , ) is the deflection function and is the time.
The boundary conditions are defined as: By using the mode superposition method, ( , ) in Equation (2) can be decomposed into two parts, The mass distribution and cross-section rigidity at x position of the composite beam are where and E 2 are mass density, cross-sectional area and modulus of elasticity of host beam and composite layer. The equation of free vibration of a non-uniform composite beam without considering damping is given below, where w(x, t) is the deflection function and t is the time.
The boundary conditions are defined as: By using the mode superposition method, w(x, t) in Equation (2) can be decomposed into two parts, where W i (x) is the i-th mode shape and q i (t) is the i-th corresponding generalized coordinate of the free vibration response or external force. For free vibration analysis solving the natural frequencies and mode shape functions, considering the orthogonality of a linear vibration problem, substituting Equation (6) into Equation (3), one ordinary differential equation corresponding to the i-th mode is obtained as, Appl. Sci. 2022, 12, 2700

of 20
where ω i is the i-th natural frequency. In order to solve the above equation, ADM is applied. The operator form of the Equation (7) is rewritten as, where L x is the fourth order differential operator. L x −1 is applied on the both sides of Equation (8) at the same time, where L x −1 is the fourth-order integral operator.
where C 1 to C 4 are constants that can be determined by boundary conditions, the detailed progress defining C 1 -C 4 can be found in [36][37][38].
where k is the number of terms in series form. The larger the value of k, the more accurate the solution. A precise solution is often obtained with very small values of k [50]. Substituting Equation (10) into Equation (9), we have, For each term in the series, we can have, The natural frequencies can be obtained by introducing the mode functions into the boundary conditions and solving the eigenvalue problem. The detailed solving process of the mode shape functions and nature frequencies of the non-uniform composite beam can be found in references [41][42][43] and is hence not provided here.

Vibration Model of Non-Uniform Composite Beam under Random Load
After the natures of the beam structure are solved by free vibration analysis and ADM, the forced vibration response can be solved using the modal analysis and Duhamel integral. The forced vibration governing equation of non-uniform composite beam considering damping under action of random load F(x, t) is given below, Appl. Sci. 2022, 12, 2700 is the distribution force generated by base motion. The displacement and acceleration of base motion are y = Y sin(ϕt) and a = −ϕ 2 Y sin(ϕt), and A a = ϕ 2 Y is the amplitude of acceleration. Y and ϕ are amplitude and angular frequency of base motion. Substituting Equation (6) into Equation (14), the following equation can be obtained, Based on the understanding from the free vibration governing equation, for the i-th mode of free vibration, we have, With the understanding from the relationship given in the above equation, both sides of Equation (15) multiplied by W j (x) and integrated from 0 to L in space domain leads to, Since in general, the damping function does not possess the orthogonality property, it is assumed that the structural damping is in the form of According to the orthogonality of normal vibration modes, Equation (18) for the i-th mode can be obtained, while αω 2 i = 2ξ i ω i is defined, where ξ i is the modal damping ratio of the corresponding i-th order natural mode. Equation (18) can be solved by the Duhamel integral, and the final time domain solution for i-th mode is, where ω d i is the damped frequency corresponding to the i-th vibration mode, is the force coefficient corresponding to the i-th vibration mode. According to Equations (12), (14) and (19), the response of the non-uniform composite beam structure is,

Iteration Process Considering Random Load
The premise of the iterative method is that at t = 0 s, the initial condition of vibration is known, and there is no force applied to the beam. At this time, the beam is considered to be at rest, and the deflection and velocity at any position of the beam are, Because some variables, like deflection, corresponding generalized coordinate, excitation, in each iteration step are different in the process of iteration of the system, they are defined by iteration step subscripts. The subscript 1 represents these variables in the first time period/iteration step in period t 1 -t 2 . The subscript 2 represents these variables in the second time period/iteration step in period t 2 -t 3 and so on. The subscript n represents these variables in the n-th time period/iteration step in period t n -t n+1 .
According to the engineering experience, the step length of each iteration is determined by base motion excitation frequency, ∆t=t n+1 − t n = 1 10× f (n) , 1 ≤ n ≤ ∞, and n is the number of iteration steps, f (n) is the base motion excitation frequency during time period t n -t n+1 .
The initial condition for time t 1 (t 1 = 0 s) is, When the time is in t 1 -t 2 , there is no free vibration and the response of the forced vibration is, To clarify the iteration progress, we start the derivation from the second time step, t 2 -t 3 , while the first step vibration solutions is only from the Duhamel integral assuming the structure is at rest before excitation as described in Equation (23). The free vibration response from the initial condition at t 2 is, A 2 i and B 2 i can be obtained by the above formula, During the period, t 2 -t 3 , the total vibration response is, where For general vibration solution during the iteration, when the time is t n -t n+1 , the total vibration response is, hence, where q f ree n i (t) is i-th corresponding generalized coordinate of free vibration during t n -t n+1 , q f orce n i (t) is i-th corresponding generalized coordinate of the external force during t n -t n+1 .
A n i and B n i are given as follows with known q n−1 i (t n ) from the previous iteration step,

Vibration Reliability Analysis
In the process of vibration, the excessive vibration response (stress or amplitude) of the structure will lead to the damage and failure of the structure or system. According to the reliability theory, the state function of general system with vibration is defined as, where ε j is the deflection of non-uniform composite beam under a random load at a certain time point, ε i is the reference safety value of the strain of the Bernoulli-Euler beam with the parameters researched in this work, γ is a specific value and is generally 10% times of ε i . The vibration reliability of the non-uniform composite beam under a random load during a period of time is, where N G ij (ε j ,ε i )≤0 is the number of iteration points satisfying the state function greater than 0, and N total is the total number of iteration points in the iterative calculation results within a period of time.

Numerical Studies, Results and Discussion
A non-uniform rectangular cantilever composite beam with positive exponential crosssection variation function is chosen for numerical case studies. The length of the beam is fixed as L = 0.2 m and the thickness of the left end of the non-uniform composite beam is H 0 = 0.003 m. The dimensions at x position along the length of the beam are b(x) = b 0 e hx , h 1 (x) = h 0 e nx and h 2 (x) = (H 0 − h 0 )e mx . h is the taper ratio of the width variation of the composite beam. n and m are taper ratios for the variations of thickness of the host beam and the composite layers, respectively. The material properties of the composite beam are ρ 1 = 2.7 × 10 3 kg/m 3 , ρ 2 = 7.8 × 10 3 kg/m 3 , E 1 = 7 × 10 10 Pa, E 2 = 68.9 × 10 9 Pa, for the host beam and composite layers, respectively. The reference value of maximum strain is ε i = 0.0035.

Model Validation
In order to prove the correctness of the mathematical model proposed in this paper, the first two order dimensionless natural frequencies of the non-uniform composite beam structure shown in Figure 1 and the non-uniform beam structure in reference [43] are calculated and compared. The thickness of the composite layer in the current model is set to be a very small value of 0.00001 m to match the structure studied in reference [43]. The parameters are shown in Table 1. The calculation results are shown in Table 2. It can be seen from Table 2 that when the thickness of the composite layer is very small, the first two natural frequencies of the two beam structures are very close, which proves the correctness of the mathematical model in this paper.

Influences of Different Design Parameters on Natural Frequencies
For composite structures, different design parameters will change the inherent properties of the structure. By adjusting the structural design parameters, the damage caused by excessive structural amplitude (resonance) can be effectively avoided. The influences of different design parameters on the natural frequencies (first three orders, unit rad/s) of non-uniform composite beam are shown in Tables 3-7. Constant parameters for each study are shown in the table captions. The variation ranges of m and n are to study the influence trend of different parameters taper ratios of different composite layers and host beam on the natural frequency of non-uniform composite non-uniform beam structure. As for the core material in the composite beam in our paper with relatively less effect on stiffness, the variation of the host beam taper ratio, n, is set in a relatively larger range to see its effect on natural frequency variation more obviously. It can be seen from the results in Tables 3-7 that within the range of parameters studied, for the three layered non-uniform composite beam with parameter variation ranges in this paper, the nature frequencies increase with the increment of h 0 , n, m and decrease a bit with the increment of h. b 0 has little effect on nature frequencies. First, for the change of h and b 0 (change of width), both stiffness and mass at certain cross section of the beam changes linearly with the beam width variation leading to less effect on natural frequency variation. However, for the ranges studied of h 0 , n and m, their increment will increase the thickness and stiffness of the beam structure in cubic function but the mass linearly. Therefore, within a certain change range, the more significant increment of stiffness compared with the mass change will lead to the increase of natural frequency. On the other hand, for the first vibration mode as an example, if n and m are further increased, the mass can be more concentrated at the beam free end, where the large beam thickness will not affect the whole modal stiffness significantly, so the natural frequency can decrease. In summary, for different non-uniform composite beams, natural frequencies variations for different modes depend on the stiffness and mass distribution following different shapes of the cross-section variation along the beam length. It is almost impossible to show all natural frequency variations of different non-uniform composite beam. The studied parameter ranges were just chosen to show a clear natural frequency variation trend.

Vibration Responses under Random Base Motion Excitation
In order to study the influences of the random base motion excitation on the vibration characteristics of non-uniform composite beam, the vibration responses of non-uniform From Figure 2a,b, while the frequency and amplitude of acceleration of random base motion excitation change with Gaussian distribution, with the increment of time, the deflection at the free end of the non-uniform composite beam will not approach a steady state with variable amplitude. This fully shows the randomness characteristics of the random load. However, for a non-uniform composite beam excited by base motion with constant frequency and amplitude of acceleration which are the average of the ones of the random load, the vibration responses (deflection) at the free end experience a process from From Figure 2a,b, while the frequency and amplitude of acceleration of random base motion excitation change with Gaussian distribution, with the increment of time, the deflection at the free end of the non-uniform composite beam will not approach a steady state with variable amplitude. This fully shows the randomness characteristics of the random load. However, for a non-uniform composite beam excited by base motion with constant frequency and amplitude of acceleration which are the average of the ones of the random load, the vibration responses (deflection) at the free end experience a process from transient to steady state with the time increment. The amplitude at the free end of beam under a random load is larger than the one under base motion with constant frequency and amplitude of acceleration, and the variation trend of deflection is unpredictable. From Figure 2c,d, PSD of time domain signal under random frequency and amplitude has obvious 'volatility'. This kind of uncontrollability explains the significance of studying the vibration reliability under a random load.

Influences of Different Ranges of Base Motion Excitation Angular Frequency on Vibration Responses
During the vibration process under random base motion excitation, the influences of exaction amplitude and frequency with different means and variances on vibration characteristics are different. Firstly, the influences of base motion excitation frequency on vibration responses were studied. The amplitude of acceleration of base motion accord with Gaussian distribution with A a ∼ N 50, 10 2 . The amplitude of base motion is Y = A a ϕ 2 . Four groups of base motion excitation angular frequencies with Gaussian distribution and standard deviation of 1/15 of the average value were selected for analysis. The specific base motion excitation frequency variation ranges are ϕ 1 ∼ N 200, (200/15) 2 rad/s, ϕ 2 ∼ N 440, (440/15) 2 rad/s, ϕ 3 ∼ N 3600, (3600/15) 2 rad/s, and ϕ 4 ∼ N 11000, (11000/15) 2 rad/s. The first three natural frequencies of non-uniform composite beam with b 0 = 0.005 m, h 0 = 0.001 m, n = 2, m = 2.5 and h = 1.5 are 439.2 rad/s, 3555.4 rad/s and 10832.3 rad/s. The frequency ranges of ϕ 2 , ϕ 3 and ϕ 4 can cover the first, second and third beam natural frequency, respectively. In addition, another four groups of base motion excitation angular frequencies with Gaussian distribution and different means and variances were selected. The specific value ranges are ϕ 5 ∼ N 200, 20 2 rad/s, ϕ 6 ∼ N 440, 44 2 rad/s, ϕ 7 ∼ N 2810, 518 2 rad/s and ϕ 8 ∼ N 8360, 1628 2 rad/s. The values of ϕ 6 , ϕ 7 and ϕ 8 can ensure the following relationships 1/2ω 1 ≤ ϕ 6 ≤ 3/2ω 1 , 1/2ω 1 ≤ ϕ 7 ≤ 3/2ω 2 , 1/2ω 1 ≤ ϕ 8 ≤ 3/2ω 3 and cover the first one, two and three beam natural frequencies, respectively. The maximum vibration amplitude and strain at all locations along the length of the beam with different variation ranges of base motion excitation frequency are shown in Figure 3.
groups of base motion excitation angular frequencies with Gaussian distribution and standard deviation of 1/15 of the average value were selected for analysis. The specific base motion excitation frequency variation ranges are 1~( 200, (200/15) 2 ) rad/s, 2~( 440, (440/15) 2 ) rad/s, 3~( 3600, (3600/15) 2 ) rad/s, and 4~( 11000, (11000/15) 2 ) rad/s. The first three natural frequencies of non-uniform composite beam with 0 = 0.005 m, ℎ 0 = 0.001 m, = 2, = 2.5 and ℎ = 1.5 are 439.2 rad/s, 3555.4 rad/s and 10832.3 rad/s. The frequency ranges of 2 , 3 and 4 can cover the first, second and third beam natural frequency, respectively. In addition, another four groups of base motion excitation angular frequencies with Gaussian distribution and different means and variances were selected. The specific value ranges are 5~( 200, 20 2 ) rad/s, 6~( 440, 44 2 ) rad/s, 7~( 2810, 518 2 ) rad/s and 8~( 8360, 1628 2 ) rad/s. The values of 6 , 7 and 8 can ensure the following relationships 1/2 1 ≤ 6 ≤ 3/2 1 , 1/2 1 ≤ 7 ≤ 3/2 2 , 1/2 1 ≤ 8 ≤ 3/2 3 and cover the first one, two and three beam natural frequencies, respectively. The maximum vibration amplitude and strain at all locations along the length of the beam with different variation ranges of base motion excitation frequency are shown in Figure 3.  As can be seen from Figure 3, when the angular frequencies of several groups of random loads cover the first, second and third and the first one, two and three order natural frequencies, the larger the mean and variance of base motion excitation frequency, the smaller the maximum amplitude and strain of beam in the same period of time can be excited. In engineering practice, the vibration amplitude of the structure can be controlled by changing the base motion excitation frequency. From the iteration numbers with different excitation frequency, we can know that when the base motion excitation frequency only covers the third order natural frequency, the high-frequency component accounts for As can be seen from Figure 3, when the angular frequencies of several groups of random loads cover the first, second and third and the first one, two and three order natural frequencies, the larger the mean and variance of base motion excitation frequency, the smaller the maximum amplitude and strain of beam in the same period of time can be excited. In engineering practice, the vibration amplitude of the structure can be controlled by changing the base motion excitation frequency. From the iteration numbers with different excitation frequency, we can know that when the base motion excitation frequency only covers the third order natural frequency, the high-frequency component accounts for the largest proportion in base motion excitation frequency within the same period of time.

Influences of Mean and Variance of Base Motion Excitation Amplitude on Vibration Responses and Reliability
The same as base motion excitation frequency, the influences of different base motion excitation amplitudes on vibration responses and reliability were studied in this section. The value ranges of frequency and amplitude of acceleration of base motion are   As can be seen from Figure 4, when the frequency and amplitude of acceleration of base motion accord with Gaussian distribution, the greater the mean and variance of the base motion excitation amplitude, the greater the value of the maximum amplitude and strain and lower reliability of beam are induced. With a larger mean value of the excitation amplitude, the upper and lower limits of the interval become larger and the reliability hence becomes smaller. The excitation with larger variance but the same mean value leads As can be seen from Figure 4, when the frequency and amplitude of acceleration of base motion accord with Gaussian distribution, the greater the mean and variance of the base motion excitation amplitude, the greater the value of the maximum amplitude and strain and lower reliability of beam are induced. With a larger mean value of the excitation amplitude, the upper and lower limits of the interval become larger and the reliability hence becomes smaller. The excitation with larger variance but the same mean value leads to a higher possibility of large strain happening on the tested structure. When the safe and failure data points appear in the same number, the failure probability will hence increase, and the reliability will decrease. The range and magnitude of vibration response of non-uniform composite beam structure can be controlled and adjusted by changing the mean and variance of excitation amplitude of the random load. The calculation results are shown in Figure 5.
As can be seen from Figure 5, the maximum amplitude and strain of beam in different periods of time is different; because the reliability is not only related to the maximum strain in this time period but also to the number of loading cycles. The result shows that in the relatively short period of time (<120 s), the reliability may increase or decrease with the increment of the time period. However, when the period of time of analysis increases above a certain value (>120 s), the reliability shows a downward trend reaching to a converged value as a whole. The results of reliability under different time lengths are different, and the calculation time duration for reliability analysis should be according to the actual operation and measurement conditions. to different vibration responses and reliability results within different time periods. The maximum amplitude and strain of beam with 0 = 0.005 m, ℎ 0 = 0.001 m, = 2, = 2.5 and ℎ = 1.5 in the time periods of 0-20 s, 20-40 s, 40-60 s, 60-80 s, 80-100 s, 100-120 s, 120-140 s, 140-160 s and 160-180 s are calculated, and vibration reliability in the time points of 20 s, 40 s, 60 s, 80 s, 100 s, 120 s, 140 s, 160 s and 180 s under the same random load with (~(50, 10 2 ) m/s 2 and ~(200, (200/15) 2 ) rad/s) were analyzed, respectively. The calculation results are shown in Figure 5. As can be seen from Figure 5, the maximum amplitude and strain of beam in different periods of time is different; because the reliability is not only related to the maximum strain in this time period but also to the number of loading cycles. The result shows that in the relatively short period of time (<120 s), the reliability may increase or decrease with the increment of the time period. However, when the period of time of analysis increases above a certain value (>120 s), the reliability shows a downward trend reaching to a converged value as a whole. The results of reliability under different time lengths are different, and the calculation time duration for reliability analysis should be according to the actual operation and measurement conditions.

Reliability Analysis of Non-Uniform Composite Beam with Different Design Parameters
In order to study the influences of different design parameters on the vibration reliability, the reliability of non-uniform beam with different taper ratios, and , were calculated under the same random load with (~(50, 10 2 ) m/s 2 and ~(200, (200/ 15) 2 ) rad/s). Other constant parameters are 0 = 0.005 m, ℎ 0 = 0.001 m and ℎ = 1.5. When studying the influence of the change of on vibration reliability, is a constant

Reliability Analysis of Non-Uniform Composite Beam with Different Design Parameters
In order to study the influences of different design parameters on the vibration reliability, the reliability of non-uniform beam with different taper ratios, n and m, were calculated under the same random load with ( A a ∼ N 50, 10 2 m/s 2   It can be seen from Figure 6 that within the range of parameters studied, the reliability of the non-uniform composite beam increases with the increment of and decreases with the increment of . It can be concluded that the taper design itself of the non-uniform composite beam will affect the structural reliability. It should be noted that the exact structural reliability variation trend will also depend on the excitation. The iteration mathematic model proposed in this paper can be used to solve and analyze the vibration responses and reliability of non-uniform composite beam with different design parameters. For the multiple layered composite beam, different design parameters lead to different It can be seen from Figure 6 that within the range of parameters studied, the reliability of the non-uniform composite beam increases with the increment of n and decreases with the increment of m. It can be concluded that the taper design itself of the nonuniform composite beam will affect the structural reliability. It should be noted that the exact structural reliability variation trend will also depend on the excitation. The iteration mathematic model proposed in this paper can be used to solve and analyze the vibration responses and reliability of non-uniform composite beam with different design parameters. For the multiple layered composite beam, different design parameters lead to different reliability results. Reasonable design under random loading condition with certain variation range of its frequency and amplitude can enhance structural reliability.

Conclusions
In this paper, an iterative method is proposed to solve the forced vibration responses and reliability of non-uniform composite beam under random load. Considering the change of frequency and acceleration amplitude of base motion excitation in the iteration process, the vibration response of beam in each short time period/iteration step is calculated by ADM and Duhamel integral considering both the transient and steady response. Through numerical studies, the results reveal the following conclusions: (1) ADM is used to be combined with iteration progress solving non-uniform composite beam structural vibration under random load. The correctness of the mathematical model is proved. Comparing solutions under random and harmonic base motion, the vibration response under random load shows the randomness and non-frequency characteristics, and its vibration amplitude is larger than the one under the same level harmonic excitation.
(2) Within the ranges of studied design parameters, the composite beam nature frequencies increase with the increment of h 0 , n, m and decrease with the increment of h. b 0 has little effect on nature frequencies.
(3) The higher mean value and variance of base motion excitation frequency, the smaller amplitude and maximum strain of the non-uniform composite beam can be excited.
(4) The greater mean and variance of the base motion excitation amplitude, the greater vibration amplitude and strain of the non-uniform composite beam, and the lower the vibration reliability in the same time period.
(5) The vibration reliability will increase or decrease with the increment of operation time considered, but the overall trend is downward if the operation time of the system is long enough.
(6) The reliability of the non-uniform composite beam increases with the increment of n and decreases with the increment of m. For the multiple layered composite beam, different design parameters lead to different reliability results. Reasonable design can enhance structural reliability.
(7) However, the research results obtained within the range of parameters studied occur when the operation time is long enough. When the calculation/analysis time duration is very short, the results may change. The statement that the operation time is long enough means that the analyzed period should be enough to ensure that the excitation amplitude and frequency distribution law can be illustrated and judged from the data points of the random loads in the calculation process.