Response to Simply Supported Orthotropic Rectangular Plate Resting on a Variable Elastic Bi-Parametric Foundation under the Action of Moving Distributed Masses

§ Federal University of Technology

Send Message

To: Author

Response to Simply Supported Orthotropic Rectangular Plate Resting on a Variable Elastic Bi-Parametric Foundation under the Action of Moving Distributed Masses

Article Fingerprint

ReserarchID

L6JES

Response to Simply Supported Orthotropic Rectangular Plate Resting on a Variable Elastic Bi-Parametric Foundation under the Action of Moving Distributed Masses Banner

Key Research Insights

Synthesized scholarly intelligence & interactive research assistant
  • English
  • Afrikaans
  • Albanian
  • Amharic
  • Arabic
  • Armenian
  • Azerbaijani
  • Basque
  • Belarusian
  • Bengali
  • Bosnian
  • Bulgarian
  • Catalan
  • Cebuano
  • Chichewa
  • Chinese (Simplified)
  • Chinese (Traditional)
  • Corsican
  • Croatian
  • Czech
  • Danish
  • Dutch
  • Esperanto
  • Estonian
  • Filipino
  • Finnish
  • French
  • Frisian
  • Galician
  • Georgian
  • German
  • Greek
  • Gujarati
  • Haitian Creole
  • Hausa
  • Hawaiian
  • Hebrew
  • Hindi
  • Hmong
  • Hungarian
  • Icelandic
  • Igbo
  • Indonesian
  • Irish
  • Italian
  • Japanese
  • Javanese
  • Kannada
  • Kazakh
  • Khmer
  • Korean
  • Kurdish (Kurmanji)
  • Kyrgyz
  • Lao
  • Latin
  • Latvian
  • Lithuanian
  • Luxembourgish
  • Macedonian
  • Malagasy
  • Malay
  • Malayalam
  • Maltese
  • Maori
  • Marathi
  • Mongolian
  • Myanmar (Burmese)
  • Nepali
  • Norwegian
  • Pashto
  • Persian
  • Polish
  • Portuguese
  • Punjabi
  • Romanian
  • Russian
  • Samoan
  • Scots Gaelic
  • Serbian
  • Sesotho
  • Shona
  • Sindhi
  • Sinhala
  • Slovak
  • Slovenian
  • Somali
  • Spanish
  • Sundanese
  • Swahili
  • Swedish
  • Tajik
  • Tamil
  • Telugu
  • Thai
  • Turkish
  • Ukrainian
  • Urdu
  • Uzbek
  • Vietnamese
  • Welsh
  • Xhosa
  • Yiddish
  • Yoruba
  • Zulu
Reading Preferences
Font Size
Line Spacing
Background

I. INTRODUCTION

The problems connected with the analysis of thin structural bodies (rods, beams, plates, and shells) with other bodies have widespread application in various fields of science and technology. The physical phenomena involved in the impact event include structural responses, contact effects and wave propagation. These problems are topical issues of research in the field of applied mechanics. Since these problems belong to the problems of dynamic contact interaction, their solution is connected with severe mathematical and calculation difficulties. To this end, several researchers had worked and some are still working on the dynamic behavior of orthotropic rectangular plates. Analytical investigation of the low-velocity impact response of circular orthotropic and transversely isotropic plates possessing curvilinear anisotropy under compressive preloading has been carried out recently by Rossikhin and Shitikova in [1] and [2], respectively. The equations of plate motion take the rotatory inertia and transverse shear deformations into account. In the case of the orthotropic target [1], the changes in the geometrical dimensions of the contact domain have been ignored and the contact interaction is modeled by a linear spring, and a force arising in it is the linear approximation of Herts?z contact force. Ambartsumian [3] examined the five fundamental differential equations describing the equilibrium of an orthotropic plate with a cylindrical anisotropy for the case when all radial planes crossing the axis of anisotropy are the planes of elastic symmetry. Sveklo [4] suggested the contact theory for two anisotropic bodies under compression according to which the contact pressure is distributed over an elliptical contact region. The same structural effects are also true of the concrete slab in a composite girder bridge, but the steel orthotropic deck is considerably lighter, and therefore allows longer span bridges to be more efficiently designed. Awodola [5] studied the effect of plate parameters on the vibrations under moving masses of elastically supported plate resting on bi-parametric foundation with stiffness variation.Szekrenyes [6] investigated the interface fracture in orthotropic composite plates using second order shear deformation theory. Kadari [7] analyzed buckling in orthotropic nanoscale plates resting on elastic foundations. Yan [8] proposed elastic orthotropic models and used these in the nonlinear analysis of concrete structures subjected to monotonic or pseudo dynamic loading. Since these models can appropriately describe the strain softening behavior of concrete beyond the peak stress and show good agreement with the strength envelope obtained from experimental results Hu and Yao [9] studied the vibration solutions of rectangular orthotropic plates by symplectic geometry method. In the same vien, Alshaya, Hunt and Rowlands [10] investigated stresses and strains in thick perforated orthotropic plates. Gbadeyan and Dada [11] found the natural frequency of rectangular plates traversed by moving concentrated masses. Awodola and Adeoye [13] investigated the behavior of simply supported orthotropic rectangular plate by applying the technique of variable separable. Adeoye and Awodola [14] studied the dynamic behavior of orthotropic rectangular plate with clamped-clamped boundary conditions by making use of the technique of Shadnam Due to inability of researchers to solve orthotropic plates problems by analytical methods, this work aims at solving the governing equation by analytical solution and also considers the effect of the flexural rigidities in both x and y directions.

II. GOVERNING EQUATION

The dynamic transverse displacement W(x,y,t) of orthotropic rectangular plates when it is resting on a bi-parametric elastic foundation and traversed by distributed mass Mr moving with constant velocity cr along a straight line parallel to the x-axis issuing from point y=s on the y-axis with flexural rigidities Dx and Dy is governed by the fourth order partial differential equation given as

Ref

Dx∂4∂x4W(x,y,t)+2B∂4∂x2∂y2W(x,y,t)+Dy∂4∂y4W(x,y,t)+μ∂2∂t2W(x,y,t)−ρhR0
[∂4∂x2∂t2W(x,y,t)+∂4∂y2∂t2W(x,y,t)]+K0(4x−3x2+x3)W(x,y,t)+S0(−13+
(2.1)12x+3x2)∂∂xW(x,y,t)−S0(12−13x+6x2+x3)[∂2∂x2W(x,y,t)+∂2∂y2W(x,y,t)]
−∑r=1N[MrgH(x−ct)H(y−s)−Mr(∂2∂t2W(x,y,t)+2cr∂2∂x∂tW(x,y,t)+cr2
∂2∂x2W(x,y,t))H(x−crt)H(y−s)]=0

where Dx and Dy are the flexural rigidities of the plate along x and y axes respectively.

Notes

(2.2)Dx=Exh312(1−νxνy),Dy=Eyh312(1−νxνy),B=DxDy+Goh36

Ex and Ey are the Young's moduli along x and y axes respectively, Go is the rigidity modulus, νx and νy are Poisson's ratios for the material such that Exνy=Eyνx, ρ is the mass density per unit volume of the plate, h is the plate thickness, t is the time, x and y are the spatial coordinates in x and y directions respectively, Ro is the rotatory inertia correction factor, Ko is the foundation constant, So shear modulus and g is the acceleration due to gravity, H(.) is the Heaviside function.

Rewriting equation (2.1), one obtains

μ∂2∂t2W(x,y,t)+μωn2W(x,y,t)=ρhR0[∂4∂x2∂t2W(x,y,t)+∂4∂y2∂t2W(x,y,t)]−2B
∂4∂x2∂y2W(x,y,t)−Dy∂4∂y4W(x,y,t)−Dx∂4∂x4W(x,y,t)+μωn2W(x,y,t)−K0(4x−1)
(2.3)3x2+x3)W(x,y,t)+Go(−13+12x+3x2)∂∂xW(x,y,t)−G0(12−13x+6x2+x3)[
∂2∂x2W(x,y,t)+∂2∂y2W(x,y,t)]+∑r=1N[MrgH(x−crt)H(y−s)−Mr(∂2∂t2W(x,y,t)
+2cr∂2∂x∂tW(x,y,t)+cr2∂2∂x2W(x,y,t))H(x−crt)H(y−s)]

Simplifying equation (2.3) further, one obtains

∂2∂t2W(x,y,t)+ωn2W(x,y,t)=∑r=1N[R0(∂4∂x2∂t2W(x,y,t)+∂4∂y2∂t2W(x,y,t))−2Bμ
∂4∂x2∂y2W(x,y,t)−Dyμ∂4∂y4W(x,y,t)−Dxμ∂4∂x4W(x,y,t)+ωn2W(x,y,t)−K0μ(4x−3x2)
+x3)W(x,y,t)+Goμ(−13+12x+3x2)∂∂xW(x,y,t)−G0μ(12−13x+6x2+x3)(∂2∂x2
W(x,y,t)=0=∂∂tW(x,y,t)

III. ANALYTICAL APPROXIMATE SOLUTION

In order to solve equation (2.4), one applies technique of Shadnam et al which requires that the deflection of the plates be in series form as

(3.1)W(x,y,t)=∑n=1NΨn(x,y)Qn(t)

where

Ψn(x,y)=Ψjm(x)Ψhm(y)
Ψjmx=sin⁡ζjmx+Ajmcos⁡ζjmx+Bjmsinh⁡ζjmx+Cjmcosh⁡ζjmx
Ψhmy=sin⁡φhmy+Ahmcos⁡φhmy+Bhmsinh⁡φhmy+Chmcosh⁡φhmy
ζjm=ϕjmLx,φhm=ϕhmLy

The right hand side of equation (2.4), taken into account equation (3.1), written in the form of series takes the form

∑n=1∞[R0(∂2∂x2Ψn(x,y)Q¨n(t)+∂4∂y2Ψn(x,y)Q¨n(t))−2Bμ∂4∂x2∂y2Ψn(x,y)Qn(t)−Dyμ∂4∂y4
Ψn(x,y)Qn(t)−Dxμ∂4∂x4Ψn(x,y)Qn(t)+ωn2Ψn(x,y)Qn(t)−K0μ(4x−3x2+x3)Ψn(x,y)
Qn(t)+Goμ(−13+12x+3x2)∂∂xΨn(x,y)Qn(t)−G0μ(12−13x+6x2+x3)(∂2∂x2Ψn(x,y)
Θn(t)=1θ∗∑n=1∞∫A[R0(∂2∂x2Ψn(x,y)Q¨n(t)+∂4∂y2Ψn(x,y)Q¨n(t))−2Bμ∂4∂x2∂y2Ψn(x,y)
Qn(t)−Dyμ∂4∂y4Ψn(x,y)Qn(t)−Dxμ∂4∂x4Ψn(x,y)Qn(t)+ωn2Ψn(x,y)Qn(t)−K0μ(4x−3x2+x3)Ψn(x,y)Qn(t)+G0μ(−13+12x+3x2)∂∂xΨn(x,y)Qn(t)−G0μ(12−13x+3x2)∂∂xΨn(x,y)Qn(t).
3x2+x3)Ψn(x,y)Qn(t)+Goμ(−13+12x+3x2)∂∂xΨn(x,y)Qn(t)−G0μ(12−13x+3x2)∂∂xΨn(x,y)Qn(t).
+6x2+x3)(∂2∂x2Ψn(x,y)Qn(t)+∂2∂y2Ψn(x,y)Qn(t))+∑r=1N(MrμgH(x−crt)H(y−s)
−Mrμ(Ψn(x,y)Q¨n(t)+2c∂∂xΨn(x,y)Q˙n(t)+cr2∂2∂x2Ψn(x,y)Qn(t))H(x−crt)H(y−s)
(3.3))]Ψm(x,y)dA

and zero when n≠m

where

(3.4)θ∗=∫AΨn2(x,y)dA

Making use of equation (3.3) and taking into account equation (3.2), equation (2.4) can be written as

Ψn(x,y)[ωn2Qn(t)+Q¨n(t)]=Ψn(x,y)θ∗∑q=1∞∫A[R0∂2Ψq(x,y)∂x2Ψm(x,y)Q¨q(t)+∂2Ψq(x,y)∂y2
Ψm(x,y)Q¨q(t))−2Bμ∂2Ψq(x,y)∂x2∂y2Ψm(x,y)Qq(t)−Dyμ∂4Ψq(x,y)∂y4Ψm(x,y)Qq(t)+Dxμ∂4Ψq(x,y)∂x4
Ψm(x,y)Qq(t)−K0μ(4x−3x2+x3)Ψq(x,y)Ψm(x,y)Qq(t)+Goμ(−13+12x+3x2)∂Ψq(x,y)∂x
Ψm(x,y)Qq(t)−Goμ(12−13x+6x2+x3)(∂2Ψq(x,y)∂x2Ψm(x,y)Qq(t)+∂2Ψq(x,y)∂y2Ψm(x,y)
Qq(t))+∑r=1N(MrμgΨm(x,y)H(x−crt)H(y−s)−Mrμ(Ψq(x,y)Ψm(x,y)Q¨q(t)+2crt)).

Notes

(3.5)∂Φq(x,y)∂xΦm(x,y)Q˙q(t)+cr2∂2Φq(x,y)∂x2Φm(x,y)Qq(t))H(x−crt)H(y−s)]dA

On further simplification of equation (3.5), one obtains

Q¨n(t)+ωn2Qn(t)=1θ∗∑q=1∞∫A[R0(∂2Ψq(x,y)∂x2Ψm(x,y)Q¨q(t)+∂2Ψq(x,y)∂y2Ψm(x,y)Q¨q(t))
−2Bμ∂2Ψq(x,y)∂x2∂y2Ψm(x,y)Qq(t)−Dyμ∂4Ψq(x,y)∂y4Ψm(x,y)Qq(t)−Dxμ∂4Ψq(x,y)∂x4Ψm(x,y)
Qq(t)+ωq2Ψq(x,y)Ψm(x,y)Qn(t)−K0μ(4x−3x2+x3)Ψq(x,y)Ψm(x,y)Qq(t)+Goμ(−13
+12x+3x2)∂Ψq(x,y)∂xΨm(x,y)Qq(t)−Goμ(12−13x+6x2+x3)(∂2Ψq(x,y)∂x2Ψm(x,y)Qq(t)
+∂2Ψq(x,y)∂y2Ψm(x,y)Qq(t))+∑r=1N(MrμgΨm(x,y)H(x−crt)H(y−s)−Mrμ(Ψq(x,y)Ψm(x,y)
(3.6)Q¨q(t)+2cr∂Φq(x,y)∂xΦm(x,y)Q˙q(t)+cr2∂2Φq(x,y)∂x2Φm(x,y)Qq(t))H(x−crt)H(y−s)]dA

The system of equations in equation (3.6) is a set of coupled ordinary differential equations where H(x−crt) and H(y−s) are the Heaviside functions which are defined as

(3.7)H(x−crt)={1,for x≥crt0,for x<crt,H(y−s)={1,for y≥s0,for y<s

With the properties

(3.8)(i)ddx[H(x−crt)]=δ(x−crt),ddy[H(y−s)]=δ(y−s)
H(y−s)=14+1π∑r=1Nsin⁡(2n+1)π(y−s)2n+1,0<y<1

On putting equations (3.7) to (3.11) into equation (3.6) and simplifying, one obtains

Q¨n(t)+ωn2Qn(t)−1θ∗∑q=1∞[R0T1Q¨q(t)−2BμT2Qq(t)−DyμT3Qq(t)−DxμT4Qq(t)+(ωq2F4∗−ωq3)],
K0μF5∗)Qq(t)+G0μ(T6+T7)Qq(t)−∑r=1NMrμ((T8+1π2(∑j=1∞E1∗cos⁡(2j+1)πcrt2j+1−∑j=1∞E2∗cos⁡(2j+1)πcrt2j+1)
sin⁡(2j+1)πcrt2j+1)(∑k=1∞E3∗cos⁡(2k+1)πs2k+1−∑k=1∞E4∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E5∗cos⁡(2j+1)πcrt2j+1
−∑j=1∞E6∗sin⁡(2j+1)πcrt2j+1)+14π(∑k=1∞E7∗cos⁡(2k+1)πs2k+1−∑k=1∞E8∗sin⁡(2k+1)πs2k+1))Q¨q(t)+
2cr(T9+1π2(∑j=1∞E9∗cos⁡(2j+1)πcrt2j+1−∑j=1∞E10∗sin⁡(2j+1)πcrt2j+1)(∑k=1∞E11∗cos⁡(2k+1)πs2k+1
−∑k=1∞E12∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E13∗cos⁡(2j+1)πcrt2j+1−∑j=1∞E14∗sin⁡(2j+1)πcrt2j+1)
+14π(∑k=1∞E15∗cos⁡(2k+1)πs2k+1−∑k=1∞E16∗sin⁡(2k+1)πs2k+1)Q˙q(t)+cr2(T10+1π2(∑j=1∞E17∗))
cos⁡(2j+1)πcrt2j+1−∑j=1∞E18∗sin⁡(2j+1)πcrt2j+1)(∑k=1∞E19∗cos⁡(2k+1)πs2k+1−∑k=1∞E20∗sin⁡(2k+1)πs2k+1)
(3.13)T1=∫A[∂2∂x2Ψq(x,y)Ψm(x,y)+∂2∂y2Ψq(x,y)Ψm(x,y)]dA
(3.14)T2=∫A∂2∂x2[∂2∂x2Ψq(x,y)]Ψm(x,y)dA
(3.15)T3=∫A∂4∂y4[Ψq(x,y)]Ψm(x,y)dA
(3.16)T4=∫A∂4∂x4[Ψq(x,y)]Ψm(x,y)dA
(3.17)F4∗=∫AΨq(x,y)Ψm(x,y)dA
(3.18)T5=4U1−3U2+U3,T6=−13A1+12A2+3A3
(3.19)T7=12f1−13f2+6f3+f4+12f5−13f6+6f7+f8
(3.20)T8=116∫AΨq(x,y)Ψm(x,y)dA
(3.21)E1∗=∫AΨq(x,y)Ψm(x,y)sin⁡(2j+1)πxdA
(3.22)E2∗=∫AΨq(x,y)Ψm(x,y)cos⁡(2j+1)πxdA
(3.23)E3∗=∫AΨq(x,y)Ψm(x,y)sin⁡(2k+1)πydA
(3.24)E4∗=∫AΨq(x,y)Ψm(x,y)cos⁡(2k+1)πydA
(3.25)E5∗=E1∗,E6∗=E2∗,E7∗=E3∗,E8∗=E4∗
(3.26)T9=116∫A∂∂xΨq(x,y)Ψm(x,y)dA
(3.27)E9∗=∫A∂∂x(Ψq(x,y))Ψm(x,y)sin⁡(2j+1)πxdA

Notes

(3.28)E10∗=∫A∂∂x(Ψq(x,y))Ψm(x,y)cos⁡(2j+1)πxdA
(3.29)E11∗=∫A∂∂x(Ψq(x,y))Ψm(x,y)sin⁡(2k+1)πydA
(3.30)E12∗=∫A∂∂xΨq(x,y)Ψm(x,y)cos⁡(2k+1)πydA
(3.31)E13∗=E9∗,E14∗=E10∗,E15∗=E11∗,E16∗=E12∗
(3.32)T10=116∫A∂2∂x2(Ψq(x,y))Ψm(x,y)dA
(3.33)E17∗=∫A∂2∂x2(Ψq(x,y))Ψm(x,y)sin⁡(2j+1)πxdA
(3.34)E18∗=∫A∂2∂x2(Ψq(x,y))Ψm(x,y)cos⁡(2j+1)πxdA
(3.35)E19∗=∫AΨq(x,y)Ψm(x,y)sin⁡(2k+1)πydA
(3.36)E20∗=∫A∂2∂x2(Ψq(x,y))Ψm(x,y)cos⁡(2k+1)πydA
(3.37)E21∗=E17∗,E22∗=E18∗,E23∗=E19∗,E24∗=E20∗

Ψm(x,y) is assumed to be the products of functions Ψpm(x)Ψbm(y) which are the beam functions in the directions of x and y axes respectively. That is

(3.38)Ψm(x,y)=Ψpm(x)Ψbm(y)
(3.41)λpm=ξpmLx,λbm=ξbmLy

Considering a unit mass, equation (3.12) can be re-written as

Q¨n(t)+ωn2Qn(t)−1θ∗∑q=1∞[R0T1Q¨q(t)−2BμT2Qq(t)−DyμT3Qq(t)−DxμT4Qq(t)+
(ωq2F4∗−K0μT5)Qq(t)+G0μ(T6+T7)Qq(t)−αϱ((T8+1π2∑j=1∞E1∗cos⁡(2j+1)πct2j+1
−∑j=1∞E2∗sin⁡(2j+1)πct2j+1)(∑k=1∞E3∗cos⁡(2k+1)πs2k+1−∑k=1∞E4∗sin⁡(2k+1)πs2k+1)+14π)
∑j=1∞E5∗cos⁡(2j+1)πct2j+1−∑j=1∞E6∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E7∗cos⁡(2k+1)πs2k+1−
∑k=1∞E8∗sin⁡(2k+1)πs2k+1))Q¨q(t)+2c(T9+1π2(∑j=1∞E9∗cos⁡(2j+1)πct2j+1−∑j=1∞E10∗)),
sin⁡(2j+1)πct2j+1)(∑k=1∞E11∗cos⁡(2k+1)πs2k+1−∑k=1∞E12∗sin⁡(2k+1)πs2k+1)+14π∑j=1∞E13∗
cos⁡(2j+1)πct2j+1−∑j=1∞E14∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E15∗cos⁡(2k+1)πs2k+1−∑k=1∞E16∗)
sin⁡(2k+1)πs2k+1))Q˙q(t)+c2(T10+1π2(∑j=1∞E17∗cos⁡(2j+1)πct2j+1−∑j=1∞E18∗sin⁡(2j+1)πct2j+1

equation (3.42) is the fundamental equation of the problem. where

(3.43)α=Mμϱ,ϱ=LxLy
(3.44)Ψm(ct)=sin⁡χmt+Amcos⁡χmt+Bmsinh⁡χmt+Cmcosh⁡χmt
(3.45)Ψm(s)=sin⁡νm+Amcos⁡νm+Bmsinh⁡νm+Cmcosh⁡νm
(3.46)χm=ϕmcLx,νm=ϕmsLy

We shall consider the situation where the orthotropic rectangular plate is simply supported at all its edges. The boundary conditions for an orthotropic rectangular plate having simple supports at all its edges are given by

(3.47)W(0,y,t)=0=W(Lx,y,t)=0
(3.48)W(x,0,t)=0=W(x,Ly,t)
(3.49)∂2∂x2W(0,y,t)=0=∂2∂x2W(Lx,y,t)=0
(3.50)∂2∂y2W(0,y,t)=0=∂2∂y2W(x,Ly,t)=0
(3.51)Ψm(0)=Ψm(Lx)
(3.52)Ψm(0)=Ψm(Ly)
(3.53)∂2∂x2Ψm(0)=∂2∂x2Ψm(Lx)
(3.57)Am=0,Bm=0,Cm=0
(3.58)Γm=mπ

On putting equations (3.55) to (3.58) into equations (3.13) to (3.37), the integrals become

(3.59)T1=−[π2q2Lx2+π2q2Ly2]∫0Lxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.60)T2=π4q4Lx2Ly2∫0Lxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.61)T3=π4q4Ly4∫0Lxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.62)T4=π4q4Lx4∫0Lysin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
T5=4∫0Lxxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy−3∫0Lxx2sin⁡qπxLx
(3.63)sin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy+∫0Lxx3sin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLy
T6=−13πqLx∫0Lxcos⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy+12πqLx∫0Lxxcos⁡qπxLx
sin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy+6π2q2Lx2∫0Lxx2sin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLy
sin⁡mπyLydy+π2q2Lx2∫0Lxx3sin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy+12π2q2Ly2
∫0Lxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy−13π2q2Ly2∫0Lxxsin⁡qπxLxsin⁡mπxLxdx
∫0Lysin⁡qπyLysin⁡mπyLydy+6π2q2Ly2∫0Lxx2sin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.65)+π2q2Ly2∫0Lxx3sin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.66)T8=116∫0Lxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.67)E1∗=∫0Lxsin⁡qπxLxsin⁡mπxLxsin⁡(2j+1)πxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.68)E2∗=∫0Lxsin⁡qπxLxsin⁡mπxLxcos⁡(2j+1)πxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.69)E3∗=∫0Lxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLysin⁡(2k+1)πydy
(3.70)E4∗=∫0Lxsin⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLycos⁡(2k+1)πydy
(3.71)E5∗=E1∗,E6∗=E2∗,E7∗=E3∗,E8∗=E4∗
(3.76)E12∗=qπLx∫0Lxcos⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLycos⁡(2k+1)πydy
(3.77)E13∗=E9∗,E14∗=E10∗,E15∗=E11∗,E16∗=E12∗
(3.78)T10=−π2q216Lx2∫0Lxcos⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.79)E17∗=−q2π2Lx2∫0Lxcos⁡qπxLxsin⁡mπxLxsin⁡(2j+1)πxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.80)E18∗=−q2π2Lx2∫0Lxcos⁡qπxLxsin⁡mπxLxcos⁡(2j+1)πxdx∫0Lysin⁡qπyLysin⁡mπyLydy
(3.81)E19∗=q2π2Lx2∫0Lxcos⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLysin⁡(2k+1)πydy
(3.82)E20∗=−q2π2Lx2∫0Lxcos⁡qπxLxsin⁡mπxLxdx∫0Lysin⁡qπyLysin⁡mπyLycos⁡(2k+1)πydy
(3.83)E21∗=E17∗E22∗=E18∗E23∗=E19∗E24∗=E20∗

On solving equations (3.59) to (3.66),(3.72) and (3.78), and substituting into equation (47), one obtains

Q¨n(t)+ωn2Qn(t)−1θ∗∑q=1∞[−R0LxLy4π2q2Lx2+π2q2Ly2)Q¨q(t)−2Bπ4q44μLxLyQq(t)−Dyπ4q4Lx4μLy3]
Qq(t)−Dxπ4q4Ly4μLx3Qq(t)+(LxLy4ωq2−K0μT5)Qq(t)+G0μ(T6+T7)Qq(t)−ϖϱ∗((
cos⁡(2j+1)πct2j+1−∑j=1∞E10∗sin⁡(2j+1)πct2j+1)(∑k=1∞E11∗cos⁡(2k+1)πs2k+1−∑k=1∞E12∗sin⁡(2k+1)πs2k+1)
)+14π(∑j=1∞E13∗cos⁡(2j+1)πct2j+1−∑j=1∞E14∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E15∗cos⁡(2k+1)πs2k+1
−∑k=1∞E16∗sin⁡(2k+1)πs2k+1))Q˙q(t)+c2(−π2q2Ly64Lx+1π2(∑j=1∞E17∗cos⁡(2j+1)πct2j+1−∑j=1∞E18∗)).
E18∗sin⁡(2j+1)πct2j+1)(∑k=1∞E19∗cos⁡(2k+1)πs2k+1−∑k=1∞E20∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E21∗)
cos⁡(2j+1)πct2j+1−∑j=1∞E22∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E23∗cos⁡(2k+1)πs2k+1−∑k=1∞E24∗
(3.84)sin⁡(2k+1)πs2k+1))Qq(t))]=∑q=1∞Mgμθ∗Φm(ct)Φm(s)

The solutions to equation (3.84) shall be obtained by considering two cases:

a) Simply Supported Orthotropic Rectangular Plate Tranversed by Moving Force

For moving force problem, one sets ϖ=0 in equation (3.84) which becomes

Q¨n(t)+(1−LxLy4θ∗)ωn2Qn(t)−1μθ∗(−μR0LxLy4(π2n2Lx2+π2n2Ly2)Q¨n(t)−2Bπ4n44LxLy)
Qn(t)−Dyπ4n4Lx4Ly3Qn(t)−Dxπ4n4Ly4Lx3Qn(t)−K0T5Qn(t)+G0μ(T6+T7)Qn(t))
−1μθ∗∑q=1,q≠=n∞(−μR0LxLy4π2q2Lx2+π2q2Ly2)Q¨q(t)−2Bπ4q44LxLyZq(t)−Dyπ4q4Lx4Ly3
Zq(t)−Dxπ4q4Ly4Lx3Qq(t)+(LxLy4μωq2−K0T5)Qq(t)−G0(T6+T7)Qq(t)
=∑q=1∞Mgμθ∗sin⁡mπsLysin⁡mπctLx

Notes

(3.85) On further simplification and re-arrangement, one obtains

[1+λμR0LxLy4π2n2Lx2+π2n2Ly2)]Q¨n(t)+(ϵn2−λ(−2Bπ4n44LxLy−Dyπ4n4Lx4Ly3−Dxπ4q4Ly4Lx3
−K0T5Qn(t)+G0(T6+T7)Qn(t)))Qn(t)−λ∑q=1,q≠=n∞(−μR0LxLy4(π2q2Lx2+π2q2Ly2)
Q¨q(t)−2Bπ4q44LxLyQq(t)−Dyπ4q4Lx4Ly3Qq(t)−Dxπ4q4Ly4Lx3Qq(t)+(LxLy4μωq2−K0T5)
(3.86)Qq(t)+G0(T6+T7)Qq(t))=∑q=1∞Mgλsin⁡mπsLysin⁡mπctLx

where

(3.87)ϵn2=(1−LxLy4θ∗)ωn2λ=1μθ∗

Consider a parameter λ∗<1 for any arbitrary mass ratio λ, defined as

(3.88)λ=λ∗1+λ∗

It can be shown that

(3.89)λ=λ∗−o(λ∗2)

Retaining only o(λ∗), one obtains

(3.90)λ=λ∗

On putting equation (3.89) into equation (3.86), rewriting and simplifying further, one obtains

)+o(λ∗2)+…))∑q=1,q≠=n∞(−μR0LxLy4π2q2Lx2+π2q2Ly2)Q¨q(t)−(2Bπ4q44LxLy+Dxπ4q4Ly4Lx3)
+Dyπ4q4Lx4Ly3+K0Ly4φ1(x)−G0Ly4(πqLxφ2(x)+π2q2Lx2+π2q2Ly2)φ3(x))))Qq(t)=
(3.91)Mgλ∗sin⁡mπsLysin⁡mπctLx

Expanding equation (3.91), and retaining only o(λ∗), one obtains

Q¨n(t)+[ωn2(1−η∗μR0LxLy4π2n2Lx2+π2n2Ly2))−η∗(2Bπ4n44LxLy−Dxπ4n4Ly4Lx3−Dyπ4n4Lx4Ly3).
+(μωn2−K0Ly4φ1(x))+G0Ly4(πnLxφ2(x)+(π2n2Lx2+π2n2Ly2)φ3(x)))]Qn(t)−η∗(1−π2n2Lx2+π2n2Ly2).
−η∗μR0LxLy4(π2n2Lx2+π2n2Ly2)∑q=1,q≠=n∞(−μR0LxLy4(π2q2Lx2+π2q2Ly2)Q¨q(t)−(2Bπ4q44LxLy
+Dxπ4q4Ly4Lx3+Dyπ4q4Lx4Ly3+K0Ly4φ1(x)+G0Ly4(πqLxφ2(x)+(π2q2Lx2+π2q2Ly2)φ3(x)))
(3.92)Qq(t)=Mgη∗sin⁡mπsLysin⁡mπctLx

Using Struble's technique, equation (3.92) can be rewritten as

(3.93)Q¨n(t)+νn2Qn(t)=0

Hence, the entire equation (3.92) becomes

(3.94)Q¨n(t)+νn2Qn(t)=Mgλ∗sin⁡mπsLysin⁡mπctLx
(3.95)νn=ϵn−λ∗2ϵn[ϵn2R0LxLy4π2n2Lx2+π2n2Ly2)+2Bπ4n44LxLy−Dyπ4n4Lx4Ly3−Dxπ4n4Ly4Lx3−Dyπ4n4Lx4Ly3].
K0LxLy4+G0LxLy4π2n2Lx2+π2n2Ly2)]

Notes is the modified frequency for simply supported orthotropic rectangular plate traversed by moving force.

On solving equation (3.94) by Laplace transformation techniques one obtains

(3.96)Qn(t)=Mgλ∗νnsin⁡mπsLy×[νnsin⁡mπcLxt−(mπcLx)sin⁡νnt(mπcLX)2−νn2]

which on inversion becomes

(3.97)W(x,y,t)=Mgλ∗νnsin⁡mπsLy×[νnsin⁡mπcLxt−(mπcLx)sin⁡νnt(mπcLX)2−νn2]×sin⁡mπxLxsin⁡mπyLy

is the transverse displacement response to a moving force of a simply supported orthotropic rectangular plate.

b) Simply Supported Orthotropic Rectangular Plate Transversed by Moving Mass

Here, one seeks solution to the entire equation (3.84). To solve this problem, one makes use of the modified asymptotic method of Struble. The equation becomes,

Q¨n(t)+νn2Qn(t)+αϱ∗∑q=1∞[(LxLy64+1π2(∑j=1∞E1∗cos⁡(2j+1)πct2j+1−∑j=1∞E2∗))],
sin⁡(2j+1)πct2j+1)(∑k=1∞E3∗cos⁡(2k+1)πs2k+1−∑k=1∞E4∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E5∗cos⁡(2k+1)πs2k+1)
cos⁡(2j+1)πct2j+1−∑j=1∞E6∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E7∗cos⁡(2k+1)πs2k+1−∑k=1∞E8∗
sin⁡(2k+1)πs2k+1))Q¨q(t)+2c(−4mLy64(q2−m2)π+1π2(∑j=1∞E9∗cos⁡(2j+1)πct2j+1−∑j=1∞E10∗
sin⁡(2j+1)πct2j+1(∑k=1∞E11∗cos⁡(2k+1)πs2k+1−∑k=1∞E12∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E13∗cos⁡(2j+1)πct2j+1−∑j=1∞E14∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E15∗cos⁡(2k+1)πs2k+1−∑k=1∞E16∗sin⁡(2k+1)πs2k+1))Q˙q(t)+c2(−π2q2Ly64Lx+1π2(∑j=1∞E17∗cos⁡(2j+1)πct2j+1−∑j=1∞E18∗sin⁡(2j+1)πct2j+1)(∑k=1∞E19∗cos⁡(2k+1)πs2k+1−∑k=1∞E20∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E21∗−∑j=1∞E22∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E23∗cos⁡(2k+1)πs2k+1−∑k=1∞E24∗sin⁡(2k+1)πs2k+1))Qq(t)]=gαθ∗mπsLysin⁡mπctLx

On further simplifications and rearrangements of equation (3.98), one obtains

Q¨n(t)+2cαϱ∗1+αϱ∗δ(i,j))(−4mLy64(q2−m2)π+1π2(∑j=1∞E9∗cos⁡(2j+1)πct2j+1−∑j=1∞E10∗sin⁡(2j+1)πct2j+1)(∑k=1∞E11∗cos⁡(2k+1)πs2k+1−∑k=1∞E12∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E13∗cos⁡(2j+1)πct2j+1−∑j=1∞E14∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E15∗cos⁡(2k+1)πs2k+1−∑k=1∞E16∗sin⁡(2k+1)πs2k+1))Q˙n(t)+11+αϱ∗δ(i,j))[νn2+c2αϱ∗−π2q2Ly64Lx+1π2(∑j=1∞E17∗cos⁡(2j+1)πct2j+1−∑j=1∞E18∗sin⁡(2j+1)πct2j+1)(∑k=1∞E19∗cos⁡(2k+1)πs2k+1−∑k=1∞E20∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E21∗cos⁡(2j+1)πct2j+1−∑j=1∞E22∗sin⁡(2j+1)πct2j+1)+14π
sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E7∗cos⁡(2k+1)πs2k+1−∑k=1∞E8∗sin⁡(2k+1)πs2k+1))Q¨q(t)
+2c(−4mLy64(q2−m2)π+1π2(∑j=1∞E9∗cos⁡(2j+1)πct2j+1−∑j=1∞E10∗sin⁡(2j+1)πct2j+1)
(∑k=1∞E11∗cos⁡(2k+1)πs2k+1−∑k=1∞E12∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E13∗cos⁡(2j+1)πct2j+1
−∑j=1∞E14∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E15∗cos⁡(2k+1)πs2k+1−∑k=1∞E16∗sin⁡(2k+1)πs2k+1
))Q˙q(t)+c2(−π2q2Ly64Lx+1π2(∑j=1∞E17∗cos⁡(2j+1)πct2j+1−∑j=1∞E18∗sin⁡(2j+1)πct2j+1)
(∑k=1∞E19∗cos⁡(2k+1)πs2k+1−∑k=1∞E20∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E21∗cos⁡(2j+1)πct2j+1
(3.99)−∑j=1∞E22∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E23∗cos⁡(2k+1)πs2k+1−∑k=1∞E24∗sin⁡(2k+1)πs2k+1
))Qq(t)]=gα(θ∗1+αρ∗δ(i,j))sin⁡mπsLysin⁡mπctLx

We shall consider a parameter α∗<1 for any arbitrary mass ratio defined by

(3.100)α=α∗1+α∗

On putting equation (3.101) into equation (3.99), expanding and retaining only o(α∗2), one obtains

Q¨n(t)+2cα∗ϱ∗(−4mLy64(n2−m2)π+1π2(∑j=1∞E9∗cos⁡(2j+1)πct2j+1−∑j=1∞E10∗sin⁡(2j+1)πct2j+1)
(∑k=1∞E11∗cos⁡(2k+1)πs2k+1−∑k=1∞E12∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E13∗cos⁡(2j+1)πct2j+1
−∑j=1∞E14∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E15∗cos⁡(2k+1)πs2k+1−∑k=1∞E16∗sin⁡(2k+1)πs2k+1)
Q˙n(t)+[νn2(1−α∗ϱ∗(LxLy64+1π2(∑j=1∞E1∗cos⁡(2j+1)πct2j+1−∑j=1∞E2∗sin⁡(2j+1)πct2j+1
)(∑k=1∞E3∗cos⁡(2k+1)πs2k+1−∑k=1∞E4∗sin⁡(2k+1)πs2k+1)+14π(∑j=1∞E5∗cos⁡(2j+1)πct2j+1
−∑j=1∞E6∗sin⁡(2j+1)πct2j+1)+14π(∑k=1∞E7∗cos⁡(2k+1)πs2k+1−∑k=1∞E8∗sin⁡(2k+1)πs2k+1))
+c2α∗ϱ∗(−π2q2Ly64Lx+1π2(∑j=1∞E17∗cos⁡(2j+1)πct2j+1−∑j=1∞E18∗sin⁡(2j+1)πct2j+1)