Published: 15 May 2015

The moment stability of response on a rectangular thin plate under in-plane narrow band random excitation

Gen Ge1
Zheng Chao Yun2
Wen Di Zhang3
1, 2Department of Mechanical Engineering, Tianjin Polytechnic University, Tianjin, China
3Department of Mechanical Engineering, Tianjin University, Tianjin, China
Corresponding Author:
Gen Ge
Views 12
Reads 6
Downloads 1125

Abstract

A dynamical model of a thin rectangular plate under in-plate parametrical narrow band noise excitation was proposed based on the elastic theory and the Galerkin’s method. Firstly, the model was perturbed applying the multiple scale method, subsequently the averaged equation in Ito form was obtained. Subsequently, the stochastic moment stability of the steady state responses was studied by Floquet theory and the moment method. Finally, the second order moment of the system was calculated. The moments can be used to estimate the statistical characters such as mean and variance of the responses. The numerical simulation was carried out and the results are consistent with the theatrical analysis.

1. Introduction

Much attention has been drawn to the nonlinear oscillation of thin plates in the past three decades. Numerous nonlinear phenomena such as bifurcations and chaos have been observed. Among these studies, the response of simply supported thin rectangular laminated plates subjected to harmonic excitation is analyzed [1]; the bifurcations and chaos of a rectangular thin plate with 1:1 internal resonance have been investigated [2]; both the local and global bifurcation of rectangular thin plate excited parametrically and externally have been calculated respectively [3, 4]; the bifurcations of a thin plate due to flow-induced vibrations are considered and a finite-dimensional analysis were given [5-7]. However, these studies about the vibration of the thin plate induced by stochastic excitations are few. The first-mode approximate response of thin rectangular plate subject to the Gauss white noise excitation with frictional boundaries has been analyzed [8]. The response of shape memory alloy thin plate subjected to Gauss white noise is detailed in reference [9]. Anyway, the studies on the thin plates subjected to narrow band noise excitation are few. As we know, the narrow band noise is a kind of more reasonable mathematical model than harmonic excitations to simulate the real condition. Thus, there remains a need for exploring the responses of the thin plate under narrow band noise excitation.

In this paper, we focused on a rectangular thin plate under narrow band noise excitation. In Section 2, the von Kármán equation and the Galerkin method were used to obtain the ordinary differential equation governing the motions of the frictionally-bounded rectangular thin plate. Then the multiple scale method was applied to get the stationary state response. Furthermore, the Ito-type stochastic differential equation was obtained by linearization, and the responses and stability of the noise included system were studied applying the moment method. Finally, the numerical simulation was carried out to verify the theoretical analysis. The results quite agree with the analysis.

2. Modeling of the plate

The length and width of the rectangular thin plate are a and b, the thickness is h. The plate is simply supported and x= 0 and x=a. The coordinate system is establish, as shown in Fig. 1. It is assumed that u, v, w represent the displacement of a middle plane point of the thin plate in the x, y and z directions. The in-plane narrow band noise excitation is imposed at y= 0 edge which has the expression that: p(t)=p0+p1cos[Ωt+δW(t)], where p0 denotes the steady pressure, Ω means the fundamental frequency of the excitation, δ denotes the strength of the noise, and W(t) is unit Wiener process. The friction force [10] is imposed at y=b edge, which is in the form of F=kt120bwy2dy, where k is the coefficient of boundary friction.

Fig. 1Model of a rectangular thin plate and the coordinate system

Model of a rectangular thin plate and the coordinate system

The transverse vibration function can be given by the von Kármán equation for the thin plate:

1
D4w+ρh2wt2+μwt=2ϕy22wx2-22ϕxy2wxy+2ϕx22wy2,
2
4ϕ=Eh2wxy2-2wx22wy2,

where 4=4/x4+24/x2y2+4/y4 is the multiple harmonic operators, ρ is the density of thin plate, μ is the damping coefficient, d is the bending rigidity, D=Eh3/12(1-ν2), E is Young’s modulus, ν is Possion’s ratio, and ϕ is the stress function:

3
2ϕy2=Nx, 2ϕx2=Ny, 2ϕxy=-Nxy,

where Nx, Ny, Nxy are the components of the in-plane stress resultant. The boundary conditions are assumed to be:

4
w=2wx2=0u=0a1E2ϕy2-ν2ϕx2-122wx2dx=0x=0, x=a.

The boundary condition of y=b edge can be expressed as:

5
0a(2ϕ/x2)dx=-p0-p1cosΩt+δWta.

The transverse deflection w of the first order modular is expressed as:

6
wt,x,y=qtsinπxasinπyb.

Substituting Eq. (6) into Eq. (2), and considering the boundary condition Eqs. (4), (5), the stress function can be solved, then applying Eq. (3) the internal force can be obtained:

7
Nx=π2f2Eh8a21-cos2πyb-vp0-p1cosΩt+δWt,
Ny=-π2f2Eh8b2cos2πxa-p0-p1cosΩt+δWt+F,
Nxy=0.

Substituting Eq. (6) into Eq. (2), and by applying the Galerkin approach, the nonlinear vibration equation is given:

8
q¨+μmq˙+Rπ4m1a2+1b22q+4π2abmq0a0bNxa2+Nyb2sin2πxasin2πyb
-Nxy2absin2πxasin2πybdydx=0.

Considering the Eq. (7), the Eq. (9) can be rewritten as:

9
q¨+ω2q+2c+βq2q˙+αq3=fqcosΩt+δWt,

where:

c=μ2m, ω2=Dπ2m1a2+1b22-p0π2mνa2+1b2,
β=3kπ28a2bm, α=π4Eh4m3a4+1b4, f=p1π2mνa2+1b2.

Considering, t*=ωt, x=u/ab, δ*=ωδ, α*=αab/ω, μ*=μ/(mω), β*=βab/ω, f*=f/ω, Ω*=Ω/ω, then substitute the above no dimensional variables into Eq. (9) and dropping the over star yields the dimensionless function as following:

10
x¨+x+ε2μ+βx2x˙+εαx3=εfxcosΩt+δWt,

where the ε denotes the mark for small perturbation terms.

3. Modeling of the plate

The uniformly approximate solution of Eq. (10) is expressed in the following form:

11
xt,ε=x0T0,T1+εx1T0,T1+,

where T0=t, T1=εt are fast and slow time scale respectively.

In this study only the first-order uniform expansion of the solution x0(T0,T1) is calculated. Denoting D0=/T0, D1=/T1, the ordinary-time derivatives can be transformed into partial derivatives as:

12
ddt=D0+εD1+d2dt2=D02+ε2D0D1+.

Substituting Eqs. (11), (12) into Eq. (10), and comparing coefficients of ε with equal powers, one obtains:

13
D02x0+x0=0,
D02x1+x1=-2D0D1x0-αx03-2μ+βx02D0x0+fx0cosΩT0+δ-Wt.

The general solution of Eq. (13) can be expressed as:

14
x0T0,T1=AT1expiT0+cc,

where the cc is the complex conjugate of its preceding terms, and A(T1) is the slowly varying amplitude of the response. Substituting Eq. (12) and Eq. (14) into Eq. (13), yields:

15
D02x1+x1=-2D1AieiT0-αA3ei3T0+3A2A-eiT0-2μiAeiT0
-βiA3ei3T0+iA2A-eiT0+12feiT0+δWT1+cc,

where the over bar stands for the complex conjugate, δ=δ-/ε. For Wiener process W(t), E[W(t)]=0, E[W2(t)]=t, one has:

16
δ-Wt=δ-εWεt=δWT1.

We focus on the primary resonances of the system Eq. (1). Introducing the detuning parameter σ as follows, Ω=2+εσ. One has:

17
Ωt=2T0+σT1.

Using the above Eq. (17) and eliminating the secular terms of Eq. (15) yields:

18
12A-eiWTδ+iσT1f-3A2A-α-2iD1A-iA2A-β-2iAμ=0.

Expressing A(T1) into the polar form:

19
AT1=aT1expiφT1.

Substituting Eq. (19) into Eq. (18) and separating the real and imaginary parts of Eq. (18) yields:

20
dadT1=-μa-18βa3-14fasinη, adηdT1=-σa+34αa3-12fcosη-δaW'T1,

where η=2φ-σT1-δW'(T1), Eq. (11) is the first-order equation governing the modulation of the amplitude and phase.

The steady state response da/dT1=dη/dT1=0 is decided by:

21
-μa--18βa-3=14fa-sinη-, -σa-+34αa-3-δa-W'T=12fa-cosη-,

where a- and η- denote the steady amplitude and phase.

4. Stability

Different definitions of stochastic stability have been reported: the almost certain stability and the moment stability. Subsequently, we discuss the almost certain stability of the trivial solution a-=0 by the linearizing Eq. (20) at (0, 0):

22
dadT1=μa-14fasinη, adηdT1=-σa-12fcosη-δaW'T1.

Let v=lna, Eq. (22) can be written as the following Ito equations:

23
dv=-μ-f4sinηdT1, dη=-σ-f2cosηdT1-δW'T1.

The largest Lyapunov exponent λmax of the corresponding solution is given by [11]:

24
λmax=limT11T1lna(T1)a(T0)=limT11T10T1-μ-f4sinηdT1=-μ-f4Esinη.

The almost certain stability of the trivial solution can be determined by the largest Lyapunov exponent λmax, when λmax<0 the trivial solution is almost certainly stable. Otherwise the trivial one is almost certainly unstable. In the deterministic case when f= 0, as known that the Lyapunov exponent of system Eq. (10) is λmax=-μ. Such that the trivial solution of Eq. (10) is stable if and only if μ>0. When μ=0.3, δ= 0.1, the variation of λmax with σ and f are illustrated in Fig. 2(a), and the corresponding level lines are observed in Fig. 2(b). From Fig. 2(a) and (b), it can be seen that λmax is a decreasing function of parameter σ. On the other hand λmax is an increasing function of parameter f. If the exciting frequency is located near the parametrical resonance frequency Ω2: the Lyapunov exponent rises up. In short, trivial solution stability of Eq. (10) is as same as the linear system. Hence there remains a demand that the nontrivial solutions steady state response caused by the nonlinearity should be studied.

Fig. 2a) The Lyapunov exponent, b) the corresponding level lines

a) The Lyapunov exponent, b) the corresponding level lines

a)

a) The Lyapunov exponent, b) the corresponding level lines

b)

The Floquet theory is applied to investigate the stability of the nontrivial steady state response a-0 by linearization [12].

Assuming a=a-+Δaη=η-+Δη, where Δa,Δη are the small perturbation terms, substituting the above equations into Eq. (11) and ignoring the high order terms, the linearization is expressed as:

25
dΔadT1=μΔa-38βa-2Δa-14fsinη-Δa-14fcosη-Δη,
dΔηdT1=32αa--δW'tΔa+12fsinη-Δη.

Substituting Eq. (21) into Eq. (25) yields the Ito type stochastic equation:

26
dΔa=-14βa-2Δa+12σa--38αa-3ΔηdT1,
dΔη=3α2a-Δa+2μ-14βa-2ΔηdT1-δdWt,

where dW(t) is the standard unit Gauss white noise.

It is significant to explore the stable conditions for the steady state solutions. The moment method is an effective tool [12] to estimate the steady state solutions statistically. Firstly, we take expectations on both sides of Eq. (27), considering the expectation of the Gauss white noise dW(t) equals to zero, one obtains:

dE[Δa]dT1=dE[Δη]dT1=0,

which yields:

27
[-14βa-2E[Δa]+12σa--38αa-3E[Δη]=0,
3α2a-E[Δa]+(2μ-14βa-2)E[Δη]=0.

The Jacoby matrix is a frequently applied approach to judge the stability of the moment, which is expressed as follows:

28
det-14βa-2-λ12σa--38αa-33α2a--(2μ+14βa-2)-λ=0.

The characteristic equation is in the form:

29
λ2-Pλ+Q=0,

where:

30
P=-2μ-12βa-2Q=9α2+β216a-4-μβ2+34ασa-2.

The roots of the Eq. (29) are:

31
λ1,2=P±P2-4Q2.

Obviously, the stable condition of the system Eq. (11) without noise is P<0,Q>0 (condition 1), which demonstrates that the solution is realizable by numerical simulation. And the condition for existing complex roots is P2-4Q<0 (condition 2), which means there exists periodical vibration. We can observe the bifurcation from Fig. 3.

Then we consider the second order moment of the steady state solutions since the first two order moments are of the most importance. Similarly, one has:

32
dE[Δa2]dT1=dE[Δη2]dT1=dE[ΔaΔη]dT1=0.

Applying the Ito’s rule of stochastic differentiation [13], one obtains:

33
E[Δa2]=δ2(4σ-3αa-2)24(4μ+βa-2)[(9α2+β2)-(12ασ-8βμ)],E[ΔaΔη]=-βδ2a-(3αa-2-4σ)2(4μ+βa-2)[(9α2+β2)-(12ασ-8βμ)]EΔη2=δ2a-29α2+2β2-43ασ-2βμ4μ+βa-29α2+β2-12ασ-8βμ.,

Considering Eqs. (14) and (32), one obtains:

34
Ea=a-,Ea2=a- 2+EΔa2, Eη=η-,Eη2=η- 2+EΔη2.

Form Eq. (33) it is easy to see that, for E[Δa2]0,E[Δη2]0, there must exit 4μ-a-2β0 and (9α2+β2)-(12ασ+8βμ)0. One finds that the existing condition for second order moment is the same as the stability condition for the first order moment. Then we can get the approximate probability density function of the amplitude with the form of normal distribution N(a-,E[Δa2]) by Eq. (36), where a- denotes the mean and E[Δa2] denotes the variance.

Fig. 3Bifurcation of the nontrivial solution with detuning frequency

Bifurcation of the nontrivial solution with detuning frequency

Fig. 4The amplitude frequency responses of the system (the solid line: theoretical solution, red circle: responses of noise free system, black square: responses of δ= 0.05 noise, solid black dots: responses of δ= 0.1 noise)

The amplitude frequency responses of the system (the solid line: theoretical solution, red circle: responses of noise free system, black square: responses of δ= 0.05 noise, solid black dots: responses of δ= 0.1 noise)

5. Numerical simulation

Letting f= 4, μ= 0.3, β= 0.2, α= 0.4, ε= 0.1, and applying the method proposed in reference [14], the numerical results of the first and second moment of Eq. (10) with σ are given in Fig. 4. For comparison, the deterministic solutions are also given.

As shown in Fig. 4, it is clear that the numerical results verify the theoretical analysis quite well if the density of noise is small. Now, we study the effect of the noise.

The noise free system has a stable periodical solution which is shown as a limit cycle in the Fig. 5(a), while the limit cycle is diffused by the noise as shown in Fig. 5(b)-(c). The noise changes the periodical vibration into a quasi-periodical motion. The stronger the strength of noise is, the wider the width of the limit cycle will be.

Eventually, we discuss the probability density of the steady amplitude influenced by the strength of the noise. All 1000 sets of narrow band noise with same fundamental frequency and strength are imposed to the system. The statistical results and the approximate normal distribution are both given in Fig. 6. From the simulation results, it is clear that large noise strength results in large variance.

Fig. 5The phase plot of the system

The phase plot of the system

a)δ= 0

The phase plot of the system

b)δ= 0.05

The phase plot of the system

c) δ= 0.1

Fig. 6The probability density of the response (circles: numerical simulation, lines: approximate normal distribution)

The probability density of the response  (circles: numerical simulation, lines: approximate normal distribution)

a) δ= 0.05, Ea-2=0.042

The probability density of the response  (circles: numerical simulation, lines: approximate normal distribution)

b)δ= 0.1, Ea-2=0.092

6. Conclusions

In this present paper, a rectangular thin plate excited by in-plane narrow band noise is studied. Some results are obtained as follows:

1) The stochastic stability of the trivial solution is determined by the system’s linear damping coefficient and the exciting amplitude.

2) The noise will diffuse the limit cycle of the noise free system. The width of the limit cycle is determined by the noise density. The stronger the noise strength is, the wider it will be.

3) Appling the moment method, the approximate probability density of the steady amplitude is obtained. With small enough noise density, the response is in form of the Gauss distribution. The results of the digital solution agree with the analysis.

This paper enriches the studies on the thin palate vibration problems due to the new exciting form. The nonlinear vibration research has been extended to stochastic vibrating field which has many potential applications in manufacturing, aerospace engineering and civil engineering.

References

  • Abe A., Kobayashi Y., Yamada G. Two-mode response of simply supported, rectangular laminated plates. International Journal of Nonlinear Mechanics, Vol. 33, Issue 4, 1998, p. 675-690.
  • Chang S. L., Bajaj A. K., Krousgrill C. M. Nonlinear vibrations and chaos in harmonically excited rectangular plates with one-to-one internal resonance. Nonlinear Dynamics, Vol. 4 Issue 5, 1993, p. 433-460.
  • Zhang W., Liu Z. M. Global dynamics of a parametrically and externally excited thin plate. Nonlinear Dynamics, Vol. 24 Issue 3, 2001, p. 245-268.
  • Yu P., Zhang W., Bi Q. S. Vibration analysis on a thin plate with the aid of computation of normal forms. International Journal of Nonlinear Mechanics, Vol. 36 Issue 4, 2001, p. 597-627.
  • Holmes P. J. Bifurcations to divergence and flutter in flow-induced oscillations: a finite-dimensional analysis. Journal of Sound and Vibration, Vol. 53 Issue 4, 1977, p. 161-174.
  • Holmes P. J., Marsden J. E. Bifurcations to divergence and flutter in flow-induced oscillations: an infinite-dimensional analysis. Automatica, Vol. 14 Issue 4, 1978, p. 367-384.
  • Holmes P. J. Center manifolds, normal forms and bifurcations of vector fields with application to coupling between periodic and steady motion. Physica D, Vol. 2, 1981, p. 449-481.
  • Ge G., Wang H. L., Xu J. Stochastic bifurcation for a thin rectangular plate subject to in-plane stochastic parametrical excitation. Journal of Vibration and Shock, Vol. 28, Issue 9, 2009, p. 91-94, (in Chinese).
  • Zhu Z. Z., Zhang Q. X. Xu J. Nonlinear dynamic characteristics and optimal control of giant magnetostrictive laminated plate subjected to in-plane stochastic excitation. Journal of Superconductivity and Novel Magnetism, Vol. 26, Issue 4, 2013, p. 1343-1347.
  • Ye M., Xiao X. L. Analytical Mechanics. Tianjin University Publisher, Tianjin, 2001.
  • Li J. R., Xu W., Ren Z. Z., et al. Maximal Lyapunov exponent and almost-sure stability for stochastic Mathieu-Duffing systems. Journal of Sound and Vibration, Vol. 28, Issue 6, 2005, p. 395-402.
  • Rong H. W., Xu W., Tong F. Principal response of Duffing oscillator to combined deterministic and narrowband random parametric excitation. Journal of Sound and Vibration, Vol. 210 Issue 4, 1998, p. 483-515.
  • Zhu W. Q., Lu M. Q., Wu Q. T. Stochastic jump and bifurcation of a duffing oscillator under narrow-band excitation. Journal of Sound and Vibration, Vol. 165, Issue 2, 1993, p. 285-304.

About this article

Received
01 March 2015
Accepted
05 May 2015
Published
15 May 2015
SUBJECTS
Mechanical vibrations and applications
Keywords
thin rectangular plate
narrow band noise
moment method
Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant No. 11402168, 11272229 and Grant No. 11302144) and Natural Science Foundation of Tianjin (Grant No. 14JCQNJC05600, 13JCYBJC17900) and the Science Foundation of Tianjin Education Committee (Grant No. 20120902).