Response of Cantilever Model with Inertia Nonlinearity under Transverse Basal Gaussian Colored Noise Excitation

Considering the curvature nonlinearity and longitudinal inertia nonlinearity caused by geometrical deformations, a slender inextensible cantilever beam model under transverse pedestal motion in the form of Gaussian colored noise excitation was studied. Present stochastic averaging methods cannot solve the equations of random excited oscillators that included both inertia nonlinearity and curvature nonlinearity. In order to solve this kind of equations, a modified stochastic averaging method was proposed. )is method can simplify the equation to an Itô differential equation about amplitude and energy. Based on the Itô differential equation, the stationary probability density function (PDF) of the amplitude and energy and the joint PDF of the displacement and velocity were studied. )e effectiveness of the proposed method was verified by numerical simulation.


Introduction
Since A. H. Nayfeh [1,2] presented and analyzed a dynamical equation governing the vibration of slender uniform Euler-Bernoulli cantilever beams, the longitudinal inertia nonlinearity and the geometrical nonlinearity in the equation aroused much attention. Some basic theoretical analysis [3][4][5][6][7] and experimental researches [8,9] on the harmonic excited cantilevers are carried out. ese works have been extended into a wide range of fields such as nonlinear vibrations, bifurcations, and chaotic motions. Among the tremendous literatures, one essential and famous conclusion is that the geometrical nonlinearity will induce hardening effect, while the longitudinal inertia nonlinearity will lead to a softening effect.
During the enthusiasm of researching the dynamics of the cantilever excited by deterministic signals, some researchers explored the responses of randomly excited cantilever models. When Er et al. [10] investigated a cantilever model subject to lateral and axial excitations being Gaussian white noise, the model was treated as a linear multipledegree-of-freedom externally excited vibrating function. Feng et al. [11][12][13][14] studied responses of cantilever model with curvature nonlinearity and longitudinal inertia nonlinearity under narrow bounded noise excitation by the method of multiple scales. In Feng et al.'s literatures, curvature nonlinearity and longitudinal inertia nonlinearity were only treated as small perturbing terms rather than strongly nonlinear terms. Ge and Yan [15] studied a cantilever model with longitudinal inertia nonlinearity and the geometrical nonlinearity under basal white noise excitation. In this paper, the inertia nonlinearity term was simplified by Taylor series expansion, and thus, accuracy was reduced to some extent.
At present, the most effective tool to deal with wideband noise excited systems is the stochastic averaging method [16][17][18][19]. At first, the standard stochastic averaging method initialized by Landau and Stratonovich [16] and Khasminskii [17] is also called as the method of quasiconservative averaging. Later, Zhu and Lin [18,19] rederived an approach named the stochastic averaging of energy envelope by considering the conservative part of the Wong-Zakai correction terms of the system. Zhu and Huang [20,21] extended Xu and Chung's [22] method to handle strongly nonlinear oscillators with lightly linear and/or nonlinear damping subject to weakly external and/or parametric excitations of wideband random processes. Ge and Li [23] have simplified Zhu and Huang's method [20,21] with a cost of slightly reducing the prediction accuracy.
Although the improved energy envelope stochastic averaging methods [20][21][22][23] are effective on randomly excited oscillators with strong geometric nonlinearities, it cannot offer effective approach to solve randomly excited systems with inertia nonlinearity in the form of € is the nonlinear stiffness term (geometric nonlinearity), c(x, _ x)is the damping term, and η(t)is the random excitation, but the present methods cannot deal with the oscillator in the form of (1 + αx 2 is the so-called inertia nonlinearity (also called as coordinate-dependent mass in some literature). us, a new method is needed. We will investigate an improved stochastic averaging method to deal with the randomly excited differential equations with both inertia nonlinearity and the geometric nonlinearity. In other words, we will extend the present method into a new field.
Literatures on the cantilever model subject to noise are mostly concentrated on the narrow-banded noise [11][12][13][14] or Gaussian white noise [15]. e advantage of the white noise model over the others is that it is mathematically simple and idealistic. Compared to the white noise, Gaussian colored noise may be a more suitable mathematical model to simulate the environment behavior. As is well known [24,25], Gaussian colored noise can be obtained by a zero-mean Gaussian white noise passing through a linear filter _ η � − λη + λξ(t). e correlation function form is as follows: where λ denotes the correlation time parameter and ξ(t) denotes a Gaussian white noise with 0 mean and intensity of 2D. e expression of equation (1) can be also seen in [24,25]. If D � 1 and λ � 0.1, 0.2, Figure 1 shows the diagram of the correlation function R. As is well known, the autocorrelation function of the white noise isR(τ) � E[ξ(t)ξ(t + τ)] � Dδ(τ), and the δ(τ) is the Dirac delta function. If the parameter λ in equation (1) is infinitely large, the colored noise will decrease to be a Gaussian with noise. In order to study the cantilever model excited by Gaussian colored noise, a stochastic averaging method which can handle the equations with inertia nonlinearity and the geometrical nonlinearity is presented. To avoid being esoteric, the deduction procedure is elaborated in Methodology. After the averaging, an Itô type stochastic differential equation together with the draft coefficient and the diffusion coefficient is presented. Subsequently, the stationary responses of the cantilever model are given. e probability density function (PDF) of the equivalent amplitude, the vibrating energy, and the joint probability density function of the displacement and velocity are obtained. Numerical simulations on the randomly excited equations are employed to verify the theoretical deductions.

Model Assumptions
A uniform inextensible cantilever was characterized by the mass per unit length m, the bending stiffness EI, length L, linear viscose damping coefficient c, and transversal displacement of the pedestal h(t). e symbol s denotes the axial spatial Lagrange coordinate, (') denotes the partial derivative with respect to s, and (·) denotes the partial derivative with respect to time t. Only the planar motions of the beam are concerned, and the gravitational effects are neglected. If h(t) is in the form of harmonic signal and the frequency of which is near the i-th mode's frequency of the cantilever, the internal resonance between modes is negligible. erefore, a well-known dynamic model governing a slender pedestal transversely excited in-extendable cantilever beam is presented as follows [15]: where ρ is the density, A is the cross section area, EI is the bending moment, c is the damping coefficient, s is the coordinate along the axial, the transverse displacement of the cantilever is expressed as w i (s, t) � q i (t)ϕ i (s), q i (t) denotes the variable separated displacement, and ϕ i (s) denotes the ith order mode functions: e mode functions are so famous that one can find it in any textbook about vibration. Although the mode function is deduced from the linearly vibrating cantilever, it is available on nonlinear vibrations [3][4][5][6][7][8][9][10]. λ i is the roots of the transcendental frequency equation 1 + cos(λ)cosh(λ) � 0, and the first three of them are λ 1 � 1.8751, λ 2 � 4 .6941, and λ 3 � 7.8548. e fourth term 2EI i is a static-type cubic nonlinearity associated with the potential energy stored in bending, and the fifth term ρA One should notice that the randomly excited system is an utterly different form excited by deterministic signals. For a harmonic signal excited cantilever, if the exciting frequency is very near the i-th mode's frequency, the internal resonance between the modes can be ignored. As is well known, the frequency of the colored noise distributes in a very wide range. erefore, every mode may be under excitation when the basal motion is stochastic. However, the amplitude of the first mode predominates the amplitudes of other higher modes. In this paper, for simplicity, the contributions of higher modes to the total amplitude are ignored, and only the first vibration mode under colored noise excitation is studied.
By introducing the dimensionless transformation ξ � (s/L), ������� ρAL 4 /EI, the dimensionless equation excited by a Gaussian colored noise was obtained as equation (4). For convenience, the dimensionless time τ is still written as t in the following texture: where η(t) is a Gaussian colored noise with a self-correlation function in the form of equation (1). e coefficients of equation (4) are as follows: At present, single dimension stochastic oscillator with only geometrical nonlinearity term α 1 x 3 has been studied sufficiently, but there is no existing method which can do stochastic averaging on this type of stochastic differential equation with both inertia nonlinearity and geometrical nonlinearity. In the following section, a new method is proposed to solve this problem.

Definitions and Assumptions.
Firstly, we defined the coordinate-dependent mass M(x) � 1 + α 2 x 2 , stiffness termg(x) � ω 2 x + α 1 x 3 , and the damping term x 2 x. Crucial presumptions are made here. M(x) and g(x) are strongly nonlinear terms, and the damping term c(x, _ x) and the excitations η(t) are weak. e kinetic energy in equation (4) is e potential energy of the conservative part of equation (4) is us, the Hamiltonian function is expressed as where A is the stochastic amplitude and τ(t) is the stochastic phase. e stochastic frequency is Substituting equation (9) into equation (8), when θ � 0, the Hamiltonian function is When θ ≠ 0, the Hamiltonian function is Obviously, at any time, the conservative Hamiltonian function remains constant, H 1 � H 2 , which yields Equation (13) is crucial in the following deduction.

Transformations.
Next, by dividing the coordinate-dependent mass M(x) � 1 + α 2 x 2 on each term of equation (4), equation (4) is transformed to be Next, equation (14) can be transformed into a pair of equations about the stochastic transient equivalent amplitude A and the stochastic phase τ(t).
Taking derivative to x in equation (9) with respect to time t, Mathematical Problems in Engineering 3 which yields Taking derivative to y in equation (9) with respect to time t, Equations (16) and (17) can be simplified to be a pair of stochastic differential equations represented by equations (18) and (19). e detailed procedure is displayed in Appendix: 3.3. Averaging on the Frequency. Next, before we do stochastic averaging on (18) and (19), the stochastic instantaneous frequency ](A, θ) should be averaged to the averaged frequency ω(A). e averaged frequency can be obtained using various methods, such as He's energy balance method [26], variational approach [27], and He's homotopy perturbation method (HPM) [28]. Here, He's energy balance method [14] is chosen to get the averaged frequency ω(A).
Taking θ � (π/4) in equation (12), still considering H 1 � H 2 , the ω(A) can be solved out as follows: us, the stochastic instantaneous frequency ](A, θ) can be displaced to be the averaged frequency ω(A) in equations (18) and (19). is approach has been proved to be effective in dealing with single-degree-of-freedom strongly nonlinear stochastic equations by Ge [23].

Stochastic Averaging Procedure.
With the transformations accomplished, equations (18) and (19) can be written as where Equation (21) is the standard form of stochastic averaging. A(t) has been proved to converge weakly to a diffusive Markov process even in a semi-infinite time interval [29,30].

Mathematical Problems in Engineering
If one take the velocity y � 0 in equation (31), the stationary PDF of displacement x is obtained.

Numerical Simulations
One precondition of applying the stochastic averaging procedure in Section 3 is that the damping and the noise strength are all small enough compared with the conservative energy so that damping coefficient and the densities of the noise are chosen as weak ones in order to ensure that the energy dissipated by damping can balance the input energy induced by the noise. e densities of white noise are set as D � 0.05 and D � 0.03, and the dimensionless damping coefficient is set as μ � 0.1. Another precondition of applying the stochastic averaging procedure is that the noise must be wideband noise. As seen from Figure 1 and equation (1), the larger the correlation time parameter is, and the narrower the noise band becomes. In this article, the correlation time parameter was set as λ � 0.1 to ensure that the excitation is wideband. For each density, totally 2000 sets of colored noise  (6)). are imposed to the system (6). Each set noise includes 40000 numbers. e dimensionless time step is set as Δt � 0.005. e numerical solutions of the oscillator equation (6) could be obtained by an order-2 stochastic Runge-Kutta algorithm [24]. e last 10000 steps are kept as the steady-state responses. For each step, the transient Hamiltonian function H value is calculated out; simultaneously, the transient am- Finally, doing statistic on total stationary n � 2000 × 10000 values of H and A, the stationary probability density function (PDF) of the amplitude is shown in Figure 2(a), and the stationary PDF of the energy is shown in Figure 2(b). Subsequently, the joint PDF of displacement x and velocity y is numerically simulated. We take 40×40 grids in a range − 0.2 ≤ x ≤ 0.2, − 0.2 ≤ y ≤ 0.2 with a space gap Δx � Δy � 0.01. en, we count the numbers n ij (i, j � 1 ∼ 40) of the n � 2000×10000 stationary (x, y) dots which fall within each grid. Finally, the joint PDF P(x, y) value of each grid can be calculated out byP(x i , y i ) � (n ij /(n × Δx × Δy)). e joint PDFs are illustrated in Figure 3. From Figures 2 and 3, we found the theoretical prediction by the new stochastic averaging method deduced in this paper coincides with the digital Monte Carlo method. Figure 4 is the stationary PDF of the displacement x.

Conclusion
A new stochastic averaging method for solving the cantilever beam model with both inertia nonlinearity and geometric nonlinearity which is excited by colored noise is proposed in this paper. e old method is extended to a new filed. e key issue of this method is balancing the input energy and diffusion energy of the Hamiltonian function.
Expressions for the drift coefficient and diffusion coefficient in the Itô function are obtained by the averaging procedure.
e probability of stationary responses of the amplitude, the Hamiltonian function, the stationary PDF of displacement x, and the joint probability of displacement and velocity are studied. e analytical approximations are consistent with the Monte Carlo numerical simulation. erefore, the presented approach is effective to solve the equations for single-degree-of-freedom (SDOF) oscillators with inertia and curvature coexisting nonlinearities.
Appendix is appendix is a detailed deduction procedure to transform equations (16) and (17) to be equations (18) and (19). It is obvious that equation (17) contains two important expressions (z(] sin θ)/zA) and (z(] sin θ)/zθ) which can be gained by taking partial derivative to equation (13).