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

Authors: Bo Li; Kai Hu; Guoguang Jin; Yanyan Song; Gen Ge
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

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. This 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. The effectiveness of the proposed method was verified by numerical simulation.


---

## Body

## 1. 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–7] and experimental researches [8, 9] on the harmonic excited cantilevers are carried out. These 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 multiple-degree-of-freedom externally excited vibrating function. Feng et al. [11–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–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–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 x¨+ω2x+gx+cx,x˙=ηt, where gx is the nonlinear stiffness term (geometric nonlinearity), cx,x˙is the damping term, and ηtis the random excitation, but the present methods cannot deal with the oscillator in the form of 1+αx2x¨+ω2x+gx+αxx˙2+cx,x˙=ηt, where the 1+αx2is the so-called inertia nonlinearity (also called as coordinate-dependent mass in some literature). Thus, 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–14] or Gaussian white noise [15]. The 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. The correlation function form is as follows:(1)Rijτ=λiDie−λiiτ,i=j,0,i≠j,where λ denotes the correlation time parameter and ξt denotes a Gaussian white noise with 0 mean and intensity of 2D. The 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.Figure 1
The correlation function of the colored 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. The 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.


## 2. Model Assumptions

A uniform inextensible cantilever was characterized by the mass per unit lengthm, the bending stiffness EI, length L, linear viscose damping coefficient c, and transversal displacement of the pedestal ht. The 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 ht 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. Therefore, a well-known dynamic model governing a slender pedestal transversely excited in-extendable cantilever beam is presented as follows [15]:(2)ρA∫0Lϕi2dsq¨i+c∫0Lϕi2dsq˙i+EI∫0Lϕi″2dsqi+2EI∫0Lϕi″2ϕi′2dsqi3+ρA∫0L∫0sϕi′2ds2dsqiq˙i2+qi2q¨i=−ρA∫0Lϕidsh¨t,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 wis,t=qitϕis, qit denotes the variable separated displacement, and ϕis denotes the i-th order mode functions:(3)ϕis=coshλis−cosλis+cosλi+coshλisinλi+sinhλisinλis−sinhλis.The 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–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. The fourth term 2EI∫0Lϕi″2ϕi′2dsqi3 is a static-type cubic nonlinearity associated with the potential energy stored in bending, and the fifth term ρA∫0L∫0sϕi′2ds2dsqiq˙i2+qi2q¨i represents longitudinal inertia nonlinearity.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 thei-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. Therefore, 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,x=q/L,η=h¨/L,τ=t/T,T=1/λ12ρAL4/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:(4)1+α2x2x¨+ω2x+α1x3+μx˙+α2x˙2x=γηt,where ηt is a Gaussian colored noise with a self-correlation function in the form of equation (1).The coefficients of equation (4) are as follows:(5)ω2=1λ14∫01ϕ1″2dξ=1,α1=2λ14∫01ϕ1″2ϕ1′2dξ=3.2713,α2=∫01∫0ξϕ1′2dξ2dξ=4.5968,γ=−∫01ϕ1dξ=0.7830,μ=cLλ12ρAEI∫01φ2dξ.At present, single dimension stochastic oscillator with only geometrical nonlinearity termα1x3 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.


## 3. Methodology

### 3.1. Definitions and Assumptions

Firstly, we defined the coordinate-dependent massMx=1+α2x2, stiffness termgx=ω2x+α1x3, and the damping term cx,x˙=μx˙+α2x˙2x. Crucial presumptions are made here. Mx and gx are strongly nonlinear terms, and the damping term cx,x˙ and the excitations ηt are weak.The kinetic energy in equation (4) is(6)T=12Mxx˙2=121+α2x2x˙2.The potential energy of the conservative part of equation (4) is(7)Ux=∫gxdx=12ω2x2+14α1x4.Thus, the Hamiltonian function is expressed as(8)Hx,x˙=T+U=121+α2x2x˙2+12ω2x2+14α1x4.Letting(9)x=Acosθ,y=−Aνsinθ,θ=Φt+τt,where A is the stochastic amplitude and τt is the stochastic phase. The stochastic frequency is(10)νA,θ=dΦdt.Substituting equation (9) into equation (8), when θ=0, the Hamiltonian function is(11)H1=12ω2A2+14α1A4.Whenθ≠0, the Hamiltonian function is(12)H2=121+α2A2cos2θA2ν2sin2θ+12ω2A2cos2θ+14α1A4cos4θ.Obviously, at any time, the conservative Hamiltonian function remains constant,H1=H2, which yields(13)νA,θsinθ=2UA−UAcosθA2MAcosθ.Equation (13) is crucial in the following deduction.


### 3.2. Transformations

Next, by dividing the coordinate-dependent massMx=1+α2x2 on each term of equation (4), equation (4) is transformed to be(14)x˙=y,y˙=−gxMx−cx,x˙Mx+γMxηt.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 tox in equation (9) with respect to time t,(15)dAdtcosθ−Asinθv+dτdt=y=−Avsinθ,which yields(16)dAdtcosθ−Asinθdτdt=0.Taking derivative toy in equation (9) with respect to time t,(17)−dAdtνsinθ−A∂νsinθ∂A×dAdt−A∂νsinθ∂θν+dτdt=−gAcosθMAcosθ−cAcosθ,−AνsinθMAcosθ+γMAcosθηt.Equations (16) and (17) can be simplified to be a pair of stochastic differential equations represented by equations (18) and (19). The detailed procedure is displayed in Appendix:(18)dAdt=cAcosθ,−AνsinθνAsinθgA+A2ν3sin2θ2gA×∂MAcosθ∂θ−γ×ηtAνsinθgA,(19)dτdt=cAcosθ,−AνsinθνcosθgA+Aν3sinθcosθ2gA×∂MAcosθ∂θ−γηt×cosθνgA.


### 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. The 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 H1=H2, the ωA can be solved out as follows:(20)ωA=3A2α1+4ω24+2A2α2.Thus, the stochastic instantaneous frequencyνA,θ can be displaced to be the averaged frequency ωA in equations (18) and (19). This approach has been proved to be effective in dealing with single-degree-of-freedom strongly nonlinear stochastic equations by Ge [23].


### 3.4. Stochastic Averaging Procedure

With the transformations accomplished, equations (18) and (19) can be written as(21)dAdt=m1A,τ+b1A,τηt,dτdt=m2A,τ+b2A,τηt,where(22)m1=cAcosθ,−AωAsinθωAAsinθgA+A2ω3Asin2θ2gA×∂MAcosθ∂θ,m2=cAcosθ,−AωAsinθωAcosθgA+Aω3Asinθcosθ2gA×∂MAcosθ∂θ,b1=−γAωAsinθgA,b2=−γωAcosθgA.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]. Thus, the Itô equation of the limiting diffusion process is of the following form:(23)dA=mAdt+σAdWt,where(24)mA=m1+∫0+∞∂b1i∂A|tb1j|t+τ+∂b1i∂τ|tb2j|t+τRijτdτσ2A=∫−∞+∞b1i|tb1j|t+τRijτdτ,where W(t) is the standard unit Brownian motion, ⋅is the averaging with respect to time, and Rijτ is the correlation function shown as equation (1). The obtained drift coefficient mA and the diffusion coefficient σ2A are as follows:(25)mA=a10A10+a8A8+a6A6+a4A4+a2A2+a04A2+A2α2A2α1+ω233A2α1+4ω2+22+A2α2λ2,where(26)a10=−6α13α2μ−9α14μ,a8=−12α13μλ2−20α12α2μω2λ2−42α13λ2μω2,a6=−40α12μω2λ2−22α1α2μω4λ2−73α12μω4−6α12α2γ2Dλ2,a4=−44α1μω4λ2−8α2μω6λ2−56α1μω6−10α1α2γ2ω2Dλ2,a2=−16μω6−16λ2μω6+8α1γ2ω2Dλ2,a0=16γ2ω4Dλ2,σ2A=γ23A2α1+4ω2Dλ2A2α1+ω223A2α1+4ω2+22+A2α2λ2.The averaged Fokker–Planck–Kolmogorov (FPK) equation associated with Itô equation (23) is of the following form [30, 31]:(27)∂p∂t=−mA∂p∂A+12σ2A∂2p∂A2,where pA,t is the transition probability density of displacement amplitude. A nontrivial stationary solution exists due to nonvanishing external excitation. Under assumption of zero probability flow at the two boundaries A=0 and A=∞, the stationary solution of FPK equation (27) [31] for system (23) is of the following form:(28)pA=Nσ2Ae2∫mA/σ2AdA,where constant N is the normalization constant. Substituting equations (25) and (26) into equation (28),(29)PA=NAA2α1+ω23A2α1+4ω2+22+A2α2λ23A2α1+4ω22+A2α21/2−1+6α12−7α1α2ω2+2α22ω4/α23γ2Dλ2eA4−3α12α2μ−2α1α22μλ2+A212α12μ−14α1α2μω2−4α22μω2λ2/8α22γ2Dλ2.The stationary probability density of total energyH or energy envelope can be obtained as follows:(30)pH=pAdAdH=pAgA|A=g−1H,where A=g−1H=−ω2/α1+4Hα1+ω4/α1. It is solved from equation (11).The joint stationary probability density of displacementx and velocity y can be further obtained from p(H) as follows:(31)px,y=pHTH|H=1/2ω2x2+1/4α1x4+1/21+α2x2y2,where T(H) can be obtained from TH=2π/ωA|A=g−1H.If one take the velocityy = 0 in equation (31), the stationary PDF of displacement x is obtained.


## 3.1. Definitions and Assumptions

Firstly, we defined the coordinate-dependent massMx=1+α2x2, stiffness termgx=ω2x+α1x3, and the damping term cx,x˙=μx˙+α2x˙2x. Crucial presumptions are made here. Mx and gx are strongly nonlinear terms, and the damping term cx,x˙ and the excitations ηt are weak.The kinetic energy in equation (4) is(6)T=12Mxx˙2=121+α2x2x˙2.The potential energy of the conservative part of equation (4) is(7)Ux=∫gxdx=12ω2x2+14α1x4.Thus, the Hamiltonian function is expressed as(8)Hx,x˙=T+U=121+α2x2x˙2+12ω2x2+14α1x4.Letting(9)x=Acosθ,y=−Aνsinθ,θ=Φt+τt,where A is the stochastic amplitude and τt is the stochastic phase. The stochastic frequency is(10)νA,θ=dΦdt.Substituting equation (9) into equation (8), when θ=0, the Hamiltonian function is(11)H1=12ω2A2+14α1A4.Whenθ≠0, the Hamiltonian function is(12)H2=121+α2A2cos2θA2ν2sin2θ+12ω2A2cos2θ+14α1A4cos4θ.Obviously, at any time, the conservative Hamiltonian function remains constant,H1=H2, which yields(13)νA,θsinθ=2UA−UAcosθA2MAcosθ.Equation (13) is crucial in the following deduction.


## 3.2. Transformations

Next, by dividing the coordinate-dependent massMx=1+α2x2 on each term of equation (4), equation (4) is transformed to be(14)x˙=y,y˙=−gxMx−cx,x˙Mx+γMxηt.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 tox in equation (9) with respect to time t,(15)dAdtcosθ−Asinθv+dτdt=y=−Avsinθ,which yields(16)dAdtcosθ−Asinθdτdt=0.Taking derivative toy in equation (9) with respect to time t,(17)−dAdtνsinθ−A∂νsinθ∂A×dAdt−A∂νsinθ∂θν+dτdt=−gAcosθMAcosθ−cAcosθ,−AνsinθMAcosθ+γMAcosθηt.Equations (16) and (17) can be simplified to be a pair of stochastic differential equations represented by equations (18) and (19). The detailed procedure is displayed in Appendix:(18)dAdt=cAcosθ,−AνsinθνAsinθgA+A2ν3sin2θ2gA×∂MAcosθ∂θ−γ×ηtAνsinθgA,(19)dτdt=cAcosθ,−AνsinθνcosθgA+Aν3sinθcosθ2gA×∂MAcosθ∂θ−γηt×cosθνgA.


## 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. The 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 H1=H2, the ωA can be solved out as follows:(20)ωA=3A2α1+4ω24+2A2α2.Thus, the stochastic instantaneous frequencyνA,θ can be displaced to be the averaged frequency ωA in equations (18) and (19). This approach has been proved to be effective in dealing with single-degree-of-freedom strongly nonlinear stochastic equations by Ge [23].


## 3.4. Stochastic Averaging Procedure

With the transformations accomplished, equations (18) and (19) can be written as(21)dAdt=m1A,τ+b1A,τηt,dτdt=m2A,τ+b2A,τηt,where(22)m1=cAcosθ,−AωAsinθωAAsinθgA+A2ω3Asin2θ2gA×∂MAcosθ∂θ,m2=cAcosθ,−AωAsinθωAcosθgA+Aω3Asinθcosθ2gA×∂MAcosθ∂θ,b1=−γAωAsinθgA,b2=−γωAcosθgA.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]. Thus, the Itô equation of the limiting diffusion process is of the following form:(23)dA=mAdt+σAdWt,where(24)mA=m1+∫0+∞∂b1i∂A|tb1j|t+τ+∂b1i∂τ|tb2j|t+τRijτdτσ2A=∫−∞+∞b1i|tb1j|t+τRijτdτ,where W(t) is the standard unit Brownian motion, ⋅is the averaging with respect to time, and Rijτ is the correlation function shown as equation (1). The obtained drift coefficient mA and the diffusion coefficient σ2A are as follows:(25)mA=a10A10+a8A8+a6A6+a4A4+a2A2+a04A2+A2α2A2α1+ω233A2α1+4ω2+22+A2α2λ2,where(26)a10=−6α13α2μ−9α14μ,a8=−12α13μλ2−20α12α2μω2λ2−42α13λ2μω2,a6=−40α12μω2λ2−22α1α2μω4λ2−73α12μω4−6α12α2γ2Dλ2,a4=−44α1μω4λ2−8α2μω6λ2−56α1μω6−10α1α2γ2ω2Dλ2,a2=−16μω6−16λ2μω6+8α1γ2ω2Dλ2,a0=16γ2ω4Dλ2,σ2A=γ23A2α1+4ω2Dλ2A2α1+ω223A2α1+4ω2+22+A2α2λ2.The averaged Fokker–Planck–Kolmogorov (FPK) equation associated with Itô equation (23) is of the following form [30, 31]:(27)∂p∂t=−mA∂p∂A+12σ2A∂2p∂A2,where pA,t is the transition probability density of displacement amplitude. A nontrivial stationary solution exists due to nonvanishing external excitation. Under assumption of zero probability flow at the two boundaries A=0 and A=∞, the stationary solution of FPK equation (27) [31] for system (23) is of the following form:(28)pA=Nσ2Ae2∫mA/σ2AdA,where constant N is the normalization constant. Substituting equations (25) and (26) into equation (28),(29)PA=NAA2α1+ω23A2α1+4ω2+22+A2α2λ23A2α1+4ω22+A2α21/2−1+6α12−7α1α2ω2+2α22ω4/α23γ2Dλ2eA4−3α12α2μ−2α1α22μλ2+A212α12μ−14α1α2μω2−4α22μω2λ2/8α22γ2Dλ2.The stationary probability density of total energyH or energy envelope can be obtained as follows:(30)pH=pAdAdH=pAgA|A=g−1H,where A=g−1H=−ω2/α1+4Hα1+ω4/α1. It is solved from equation (11).The joint stationary probability density of displacementx and velocity y can be further obtained from p(H) as follows:(31)px,y=pHTH|H=1/2ω2x2+1/4α1x4+1/21+α2x2y2,where T(H) can be obtained from TH=2π/ωA|A=g−1H.If one take the velocityy = 0 in equation (31), the stationary PDF of displacement x is obtained.


## 4. Numerical Simulations

One precondition of applying the stochastic averaging procedure in Section3 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. The 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 are imposed to the system (6). Each set noise includes 40000 numbers. The dimensionless time step is set as Δt=0.005. The numerical solutions of the oscillator equation (6) could be obtained by an order-2 stochastic Runge–Kutta algorithm [24]. The 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 amplitude A=−ω2/α1+4Hα1+ω4/α1 is solved.Finally, doing statistic on total stationaryn = 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. Then, we count the numbers niji,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 byPxi,yi=nij/n×Δx×Δy. The 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.Figure 2
Stationary probability density of displacement amplitude and Hamilton: (a) PDF of amplitude A; (b) PDF of Hamiltonian function H (solid lines: stochastic averaging equations (29) and (30); dots: Monte Carlo simulation on equation (6)).
(a)(b)Figure 3
Stationary joint probability densityλ = 0.1: (a) D = 0.05, stochastic averaging equation (31); (b) D = 0.05, Monte Carlo simulation; (c) D = 0.03, stochastic averaging equation (31); (d) D = 0.03, Monte Carlo simulation.
(a)(b)(c)(d)Figure 4
Stationary probability density of displacementx (D = 0.03 and λ = 0.1) (solid lines: stochastic averaging; dots: Monte Carlo simulation on equation (6)).


## 5. 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. The old method is extended to a new filed. The 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. The probability of stationary responses of the amplitude, the Hamiltonian function, the stationary PDF of displacementx, and the joint probability of displacement and velocity are studied. The analytical approximations are consistent with the Monte Carlo numerical simulation. Therefore, the presented approach is effective to solve the equations for single-degree-of-freedom (SDOF) oscillators with inertia and curvature coexisting nonlinearities.



---

*Source: 1823596-2021-01-20.xml*