Acessibilidade / Reportar erro

Viscoelastic Substrates Effects on the Elimination or Reduction of the Sandwich Structures Oscillations Based on the Kelvin-Voigt Model

Abstract

Effects of viscoelastic substrates on the sandwich structures oscillations are examined in this paper. In this regard, dynamic response of sandwich annular panels with FG polar orthotropic face sheets resting on viscoelastic substrate is presented. Young’s modulus in the radial and circumferential direction, shear modulus and density of each face sheet may be varied continuously in the radial direction. The viscoelastic substrate is modeled as Kelvin-Voigt foundation. To describe more accurately response of sandwich structures, the governing dynamical equations are derived based on the layerwise theory and five systems of second order coupled partial differential equations are obtained. The effects of the stiffness and damping coefficients of the foundation on the dynamic behavior of sandwich plate are investigated for various transient loads and boundary condition. Since no available results may be found in literature to demonstrate the efficiency and accuracy of the obtained results, the obtained results are verified by comparison with finite element results based on the three dimensional theory of elasticity for some special cases.

Keywords:
Dynamic response; FG polar orthotropic; viscoelasticfoundation; transient loads, Kelvin-Voigt

1 INTRODUCTION

Sandwich structures are widely used in the industries. In practical applications, sandwich structures are subjected to various transient dynamic loads where reduce vibration of these structures underdynamic loads is important in many engineering fields and viscoelastic substrate can be applied to reducethe vibration of structures. In addition, structures resting on viscous, elastic or viscoelastic foundationare extensively used in many engineering fields. However, most of the performed studies arelimited to the static and free vibration analyses of sandwich plates and rare researches are available onthe dynamic analysis especially for plates under viscoelastic foundation.

The dynamic response of a plate of infinite extent on a viscous Winkler foundation subjected to movingtandem-axle loads with amplitude variation was investigated by Kim and McCullough (2003Kim, S.M., McCullough, B.F., (2003). Dynamic response of plate on viscous Winkler foundation to moving loads of varyingamplitude, Engineering Structures 25: 1179-1188.), based onthe classical plate theory. Liang et al. (2014Liang, X.u., Wang, Z., Wang, L., Liu, G., (2014). Semi-analytical solution for three-dimensional transient response of functionallygraded annular plate on a two parameter viscoelastic foundation, Journal of Sound and Vibration 333: 2649-2663.) studied the transient response of functionally graded annularplate under a two parameter viscoelastic foundation based on the three-dimensional theory ofelasticity and using the differential quadrature method, state space method and Laplace transform.Vibration reduction of composite plates by piezoelectric patches was investigated by Bargh and Sadra (2014), based on the classical plate theory and using a modified artificial bee colony algorithm. Basedon the sinusoidal shear deformation theory and by using the Navier's and meshless methods, free vibrationanalyses of viscoelastic double-bonded polymeric nanocomposite plates are investigated byMohammadimehr et al. (2015Mohammadimehr, M., Navi, B.R., Arani, A.G., (2015). Free vibration of viscoelastic double-bonded polymeric nanocompositeplates reinforced by FG-SWCNTs using MSGT, sinusoidal shear deformation theory and meshless method, CompositeStructures 13: 654-671.). Forced vibration analysis of a single layer viscoelastic graphene sheetresting on viscoelastic medium was performed by Hosseini Hashemi et al. (2015Hosseini Hashemi, S.h., Mehrabani, H., Ahmadi-Savadkoohi, A., (2015). Forced vibration of nanoplate on viscoelastic substrate with consideration of structural damping: An analytical solution, Composite Structures 133: 8-15.). The governing equationis derived using Hamilton’s principle based on the classical plate theory and viscoelastic KelvinVoigt model. Luong-Van et al. (2014Luong-Van, H., Nguyen-Thoi, T., Liu, G.R., Phung-Van, P., (2014). A cell-based smoothed finite element method usingthree-node shear-locking free Mindlin plate element (CS-FEM-MIN3) for dynamic response of laminated compositeplates on viscoelastic foundation, Engineering Analysis with Boundary Elements 42: 8-19.) investigated the dynamic responses of sandwich and laminatedcomposite plates resting on viscoelastic foundation based on the first-order shear deformation theoryand finite element method. Arani et al. carried out the free vibration analysis of the coupled system ofdouble-layered rectangular (2012) and annular (2014) graphene sheets embedded in a visco-Pasternak foundation was by. The motion equations were derived based on the Hamilton’s principleand first-order shear deformation theory. The differential quadrature method was applied to obtain thefrequency ratio for various boundary conditions. Lepoittevin and Kress (2011Lepoittevin, G., Kress, G., (2011). Finite element model updating of vibrating structures under free-free boundary conditionsfor modal damping prediction, International Journal of Mechanical Sciences 25: 2203-2218.) proposed a method topredict resonance frequencies and modal loss factors of bare and damped samples, using constrainedlayer damping treatment, under free-free boundary condition. Karim and Chen (2012Karim, K.R., Chen, G.D., (2012). Surface damping effect of anchored constrained viscoelastic layers on the flexural response of simply supported structures, Mechanical Systems and Signal Processing 27: 419-432.) investigated thesurface damping effects of anchored constrained viscoelastic layers on the flexural response of simplysupported Euler beams or plate strips under base excitations. In order to reduce the structural vibrationsof a mechanical system used in aerospace, a new sketch of eddy current damper (ECD) was proposedby Pan et al. (2016Pan, Q., Tian He, T., Xiao, D., Liu, X., (2016). Design and Damping Analysis of a New Eddy Current Damper for AerospaceApplications, Latin American Journal of Solids and Structures 13: 1997-2011). Zenkour (2016Zenkour, A.M., (2016). Nonlocal transient thermal analysis of a single-layered graphene sheet embedded in viscoelastic medium, Physica E 79: 87-97.) applied the classical plate theory and Fourier transformmethod for transient thermal analysis of a single layered graphene sheet embedded in viscoelastic medium.The governing dynamical equation was obtained and solved for simply-supported grapheme sheet. By using the Rayleigh-Ritz method, Plattenburg et al. (2016Plattenburg, J., Dreyer, J. Singh, T. R., (2016). A new analytical model for vibration of a cylindrical shell and cardboardliner with focus on interfacial distributed damping. Mechanical Systems and Signal Processing 75: 176-195.) proposed a new analytical modelfor a thin cylindrical shell that utilizes a homogeneous cardboard liner to increase modal damping. Theproposed theory, incorporated material structural damping along with frequency-dependent viscousand Coulomb interfacial damping formulations for the shell-liner interaction. By using the Galerkinweighted residual method and the classical plate theory, Kiasat et al. (2014Kiasat, M.S., Zaman, H.A., Aghdam, M.M., (2014). On the transient response of viscoelastic beams and plates on viscoelasticmedium, International Journal of Mechanical Sciences 83: 133-145.) analyzed the transientresponse of viscoelastic beams and plates on viscoelastic medium. Based on Kelvin-Voigt model, freedamped vibration analysis of plates with hybrid material foundation viscoelasticity was performed by Zamani et al. (2017Zamani H.A., Aghdam, M.M., Salehi, M., (2017). Free damped vibration analysis of Mindlin plates with hybrid material foundation viscoelasticity. International Journal of Mechanical Sciences, 121: 33-43.). Free vibration analysis of a simply supported viscoelastic orthotropic nanoplates resting on viscoelastic medium was studied by Pouresmaeeli et al. (2013Pouresmaeeli, S., Ghavanloo, E., Fazelzadeh, S.A., (2013). Vibration analysis of viscoelastic orthotropic nanoplates restingon viscoelastic medium. Composite structures 96: 405-410.), using the classical plate theoryand the Navier method. Zhang et al. (2016Zhang, D.P., Lei, Y., Shen, Z.B., (2016). Free transverse vibration of double-walled carbon nanotubes embedded in viscoelastic médium, Acta Mechanica 227(12): 3657-3670.) analyzed the free transverse vibration of double walledcarbon nanotubes embedded in viscoelastic medium by using the Euler-Bernoulli beam theory. Luo etal. (2016Luo, W.L., Xia, Y., Zhou, X.Q., (2016). A closed-form solution to a viscoelastically supported Timoshenko beam under harmonicline load, Journal of Sound and Vibration 369: 109-118.) presented a closed-form solution to a viscoelastically supported Timoshenko beam under aharmonic line load. The differential governing equations were converted into algebraic equations byassuming the deflection and rotation of the beam in harmonic forms with respect to time and space.Some researchers analyzed beams resting on viscoelastic foundations subjected to moving loads. Thevibration instability analysis of an oscillator moving along a Timoshenko beam was performed byMetrikine and Verichev (2001Metrikine, A.V., Verichev, S.N., (2001). Instability of vibrations of a moving two-mass oscillator on a flexibly supported Timoshenko beam, Archives of Civil and Mechanical Engineering 71: 613-624.) and Mazilu et al. (2012Mazilu, T., Dumitriu, M., Tudorache, C., (2012). Instability of an oscillator moving along a Timoshenko beam on viscoelasticfoundation, Nonlinear Dynamics 67: 1273-1293.). The inHluence of a nonlinear foundation on thedynamic response of a periodically supported beam was investigated by Hoang et al. (2016Hoang, T., Duhamel, D., Foret, G., Yin, H.P., Cumunel, G., (2016). Response of a periodically supported beam on a nonlinearfoundation subjected to moving loads, Nonlinear Dynamics 86(2): 953-961.), based on the Euler-Bernoulli beam theory. Ding et al. (2013Ding, H., Shi, K.L., Chen, L.Q., Yang, S.P. (2013). Dynamic response of an infinite Timoshenko beam on a nonlinear viscoelastic foundation to a moving load, Nonlinear Dynamics 73: 285-298.) investigated the dynamic response of infinitebeams supported by nonlinear viscoelastic foundations. The differential equations were obtained byemploying the Timoshenko beam theory and were solved based on the Adomian decomposition methodand a perturbation method in conjunction with complex Fourier transformation. The dynamic response of finite Timoshenko beamsresting on a six parameter foundation was studied by Yang et al. (2013Yang, Y., Ding, H., Chen, L.Q., (2013). Dynamic response to a moving load of a Timoshenko beam resting on a nonlinearviscoelastic foundation, Acta Mechanica Sinica 29(5): 718-727.).

Most of the existing researches were performed based on the equivalent single layer theories, whereusing the equivalent single-layer theories for analysis of sandwich plates may be inaccurate or erroneousin most circumstances. Various theories are presented for analysis of sandwich structure Carrera and Brischetto (2008Carrera, E., Brischetto, S. (2008). A survey with numerical assessment of classical and refined theories for the analysis of sandwich plates. Applied Mechanics Reviews 62(1): 010803-010819.).Many researchers investigated the laminated composite and sandwich plates based on the sandwichand multilayered structures theories. Thermoelastic bending of functionally graded sandwich plates were analyzed by Houari et al (2013Houari, M. S. A., Tounsi, A., Bég O. A., (2013). Thermoelastic bending analysis of functionally graded sandwich plates using a new higher order shear and normal deformation theory, International Journal of Mechanical Sciences 76: 102-111) based on a new higher order shear and normal deformation theory and by Tounsi et al. (2013Tounsi, A., Houari, M. S. A.,Benyoucef, S.,Bedia, E. A. A., (2013). A refined trigonometric shear deformation theory for thermoelastic bending of functionally graded sandwich plates, Aerospace Science and Technology24 (1): 209-220.) based on A refined trigonometric shear deformation theory.Based on layerwise formulation, Alipour (2016aAlipour, M.M., (2016a). A novel economical analytical method for bending and stress analysis of functionally gradedsandwich circular plates with general elastic edge conditions, subjected to various loads, Composite Part B: Engineering95: 48-63.) presented anovel economical analytical method for bending and stress analysis of functionally graded sandwichcircular plates subjected to various loads with general elastic edge conditions. Alipour (2016bAlipour, M.M., (2016b). Effects of elastically restrained edges on FG sandwich annular plates by using a novel solutionprocedure based on layerwise formulation, Archives of Civil and Mechanical Engineering 16: 678-694.) analyzed the effects of elastically restrained edges on FG sandwich annular plates by using a novel solutionprocedure and layerwise method. Alipour (2018Alipour, M.M. (2018). Transient forced vibration response analysis of heterogeneous sandwich circular plates under viscoelastic boundary support, Archives of Civil and Mechanical Engineering 18: 12-31.) examined the dynamic response of sandwich plate with viscoelastic boundary support. Bending analysis of laminated compositeplates was performed based on a predictor-corrector approach and the zig-zag theory by Lee andCao (1996Lee, K.H., Cao, L., (1996). A predictor-corrector zig-zag model for the bending of laminated composite plates, InternationalJournal of Solids Structures 33(6): 879-897.). Alipour and Shariyat (2015Alipour, M.M., Shariyat, M. (2015). Analytical zigzag formulation with 3D elasticity corrections for bending and stress analysis of circular/annular composite sandwich plates with auxetic cores, Composite Structures 132: 175-197.) employed a zigzag-elasticity plate theory for bending andstress analysis of circular/annular sandwich plates with orthotropic composite face sheets and auxeticcores. Free vibration analysis of circular and annular composite sandwich plates with auxetic cores wasperformed by Shariyat and Alipour (2017aShariyat, M., Alipour, M.M., (2017a). Analytical layerwise free vibration analysis of circular/annular composite sandwichplates with auxetic cores, International Journal of Mechanics and Materials in Design 13, 125-157.).

Although various studies have been presented for dynamic analysis of single layer, laminated compositeand sandwich plates, most of the existing studies were performed for plates subjected to load caseswith constant or specific amplitude. In this study dynamic response of sandwich annular plates withfunctionally graded polar orthotropic face sheets subjected to various dynamic loads are presented, asfirst time. Based on the combination of the Hinite Taylor’s transform (Alipour (2016cAlipour, M.M. (2016c). An analytical approach for bending and stress analysis of cross/angle-ply laminated composite plates under arbitrary non-uniform loads and elastic foundations, Archives of Civil and Mechanical Engineering 16(2): 193-210.), Shariyat and Alipour (2015Shariyat, M., Alipour, M.M., (2015). Novel Layerwise Shear Correction Factors for Zigzag Theories of Circular SandwichPlates with Functionally Graded Layers, Latin American Journal of Solids and Structures 12: 1362-1396., 2017bShariyat, M., Alipour, M.M., (2017b). Analytical bending and stress analysis of variable thickness FGM auxetic conical/cylindrical shells with general tractions, Latin American Journal of Solids and Structures 14: 805-843.)) and the fourth-order Runge-Kutta procedure, a semi-analytical solution procedureis developed for solution of the relatively complicated second order coupled partial differentialequations. It is worth to be mentioned that modern structures are designed based on the use of compositematerials which are actually anisotropic and inhomogeneous (Peng and Li (2012Peng, X.L., Li, X.F., (2012). Elastic analysis of rotating functionally graded polar orthotropic disks, International Journal ofMechanical Sciences 60: 84-91.) and Nie andBatra (2010Nie, G.J., Batra, R.C., (2010). Static deformations of functionally graded polar-orthotropic cylinders with elliptical innerand circular outer surfaces, Composites Science and Technology 70: 450-457.)). In the presented analysis, Young’s modulus in the radial and circumferential direction, shear modulus and density of each face sheet may be varied continuously in the radial direction. Thesandwich plate may be resting on viscous, elastic or viscoelastic foundation. To describe more accuratelyresponse of sandwich plates, the governing differential equations of motion are derived based on theminimum total potential energy principle by using the layerwise theory. The dynamic responses ofsandwich plate are examined for various stiffness and damping coefficients of the foundation, transientloads and boundary conditions. Since no existing work has been performed on the dynamic analysis ofFG polar orthotropic sandwich plates, accuracy and efficiency of the presented analysis are verified bycomparing the obtained results with results of the three-dimensional theory of elasticity extracted fromthe ABAQUS software for some special cases. The comparisons show that there is a very good agreementbetween present results and results of the three-dimensional theory of elasticity.

2 DESCRIPTION OF THE MATERIAL PROPERTIES, DISPLACEMENT, STRAIN AND STRESS FIELDS

As shown in Figure 1, a three-layer sandwich annular plate with functionally graded polar orthotropic face sheets resting on the viscoelastic substrate which is modeled as Kelvin-Voigt foundation is considered.

Figure 1:
Schematic of the sandwich annular plate on viscoelastic substrate.

Viscoelastic substrate is modeled as a continuously distributed medium with stiffness K w and damping coefficient C t .

Face sheets may be fabricated from functionally graded polar orthotropic materials. Young’s modulus in the radial and circumferential direction, shear modulus and density of each face sheet may be varied continuously in the radial direction according to a power-law fraction, as follows:

{ E r ( k ) ( r ) E θ ( k ) ( r ) G r z ( k ) ( r ) } = { E ¯ r ( k ) E ¯ θ ( k ) G ¯ r z ( k ) } [ 1 + α ( k ) ( r r o ) β ( k ) ] , ρ ( k ) ( r ) = ρ ¯ ( k ) [ 1 + γ ( k ) ( r r o ) η ( k ) ] k = 1,3 (1)

v r θ E r = v θ r E θ (2)

where α (k) , β (k) ,γ (k) and η (k) (k = 1, 3) are inhomogeneity parameters for top (k=1) and bottom (k=3) face sheets.

The three transverse local coordinates ξ (1), ξ (2) and ξ (3) are represented for top, core, and bottom layers, respectively. The local coordinates are measured from the mid plane of the corresponding layer and are positive upward.

Based on the layerwise theory with the piecewise linear local components, after some manipulations and imposing the continuity conditions of the displacement components at the interfaces between face sheets and core, the displacement field of the layers may be written as:

u 1 = u 0 + ( ξ ( 1 ) + h 1 2 ) ψ r ( 1 ) + h 2 2 ψ r ( 2 ) u 2 = u 0 + ξ ( 2 ) ψ r ( 2 ) u 3 = u 0 + ( ξ ( 3 ) h 3 2 ) ψ r ( 3 ) h 2 2 ψ r ( 2 ) , h i 2 ξ ( i ) h i 2 i = 1, 2, 3 (3)

Where u 0 is the radial displacement component of the mid plane of the core and ψr(i) are the local rotation of the layers of the plate.

Cauchy’s strain-displacement relations are:

ε r = u , r ε θ = u r ε r z = u , z + w , r (4)

where the symbol “,” stands for the partial derivative. Based on Hooke's generalized stress-strain law:

σ r ( k ) = E r ( k ) ( r ) 1 v r θ ( k ) v θ r ( k ) ε r ( k ) + v θ r ( k ) E r ( k ) ( r ) 1 v r θ ( k ) v θ r ( k ) ε θ ( k ) σ θ ( k ) = v θ r ( k ) E r ( k ) ( r ) 1 v r θ ( k ) v θ r ( k ) ε r ( k ) + E θ ( k ) ( r ) 1 v r θ ( k ) v θ r ( k ) ε θ ( k ) τ r z ( k ) = G r z ( k ) ( r ) γ r z ( k ) (5)

where the E, G and υ symbols denote Young’s modulus, shear modulus, and Poisson’s ratio, respectively.

Based on Eqs.(1) to (5), the stress-displacement relations of face sheets and core may be written as:

σ r ( 1 ) = E ¯ r ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] 1 v r θ ( 1 ) v θ r ( 1 ) [ u 0, r + v θ r ( 1 ) r u 0 + ( ξ ( 1 ) + h 1 2 ) ( ψ r , r ( 1 ) + v θ r ( 1 ) r ψ r ( 1 ) ) + h 2 2 ( ψ r , r ( 2 ) + v θ r ( 1 ) r ψ r ( 2 ) ) ] σ θ ( 1 ) = E ¯ θ ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] 1 v r θ ( 1 ) v θ r ( 1 ) [ v r θ ( 1 ) u 0, r + 1 r u 0 + ( ξ ( 1 ) + h 1 2 ) ( v r θ ( 1 ) ψ r , r ( 1 ) + 1 r ψ r ( 1 ) ) + h 2 2 ( v r θ ( 1 ) ψ r , r ( 2 ) + 1 r ψ r ( 2 ) ) ] τ r z ( 1 ) = G ¯ r z ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] ( ψ r ( 1 ) + w , r ) σ r ( 2 ) = E ( 2 ) 1 v ( 2 ) 2 [ u 0, r + v ( 2 ) r u 0 + ξ ( 2 ) ( ψ r , r ( 2 ) + v ( 2 ) r ψ r ( 2 ) ) ] σ θ ( 2 ) = E ( 2 ) 1 v ( 2 ) 2 [ v ( 2 ) u 0, r + 1 r u 0 + ξ ( 2 ) ( v ( 2 ) ψ r , r ( 2 ) + 1 r ψ r ( 2 ) ) ] τ r z ( 2 ) = G ( 2 ) ( ψ r ( 2 ) + w , r ) σ r ( 3 ) = E ¯ r ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] 1 v r θ ( 3 ) v θ r ( 3 ) [ u 0, r + v θ r ( 3 ) r u 0 + ( ξ ( 3 ) h 3 2 ) ( ψ r , r ( 3 ) + v θ r ( 3 ) r ψ r ( 1 ) ) h 2 2 ( ψ r , r ( 2 ) + v θ r ( 3 ) r ψ r ( 2 ) ) ] σ θ ( 3 ) = E ¯ θ ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] 1 v r θ ( 1 ) v θ r ( 1 ) [ v r θ ( 3 ) u 0, r + 1 r u 0 + ( ξ ( 3 ) h 3 2 ) ( v r θ ( 3 ) ψ r , r ( 3 ) + 1 r ψ r ( 3 ) ) h 2 2 ( v r θ ( 3 ) ψ r , r ( 2 ) + 1 r ψ r ( 2 ) ) ] τ r z ( 3 ) = G ¯ r z ( 1 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ( ψ r ( 3 ) + w , r ) (6)

3 THE GOVERNING EQUATIONS OF MOTION OF THE SANDWICH PLATE WITH FG POLAR ORTHOTROPIC FACE SHEETS

The governing equations of motion of the sandwich annular plate with functionally graded polar orthotropic face sheets on viscoelastic substrate are derived based on using the minimum total potential energy principle:

δ Π = δ U + δ K δ W = 0, (7)

where δU, δK and δW are increments of the strain energy, kinetic energy and work done by external applied loads, respectively:

δ U = V ( σ r δ ε r + σ θ δ ε θ + τ r z δ γ r z ) d V δ K = V ρ ( u ¨ δ u + w ¨ δ w ) d V δ W = A P δ w d A (8)

For sandwich plate subjected to transient load (q(t)) resting on viscoelastic substrate based on the Kelvin-Voigt model, the distributed transverse load can be defined as follows.

P = q ( t ) K w w ( r , t ) C t w ˙ ( r , t ) (9)

In which q(t) is the external dynamic loads, K w and C t are the spring and damper constant of viscoelastic foundation.

Substituting Eq. (8) into Eq. (7) and by using Eqs. (3) and (6), the following five systems of second order coupled partial differential equations are resulted after some manipulations that are not included here for brevity:

k = 0 3 ( N r ( k ) N θ ( k ) r + N r , r ( k ) ) = ( I 0 ( 1 ) + I 0 ( 2 ) + I 0 ( 3 ) ) u ¨ 0 + h 1 2 I 0 ( 1 ) ψ ¨ r ( 1 ) + h 2 2 ( I 0 ( 1 ) I 0 ( 3 ) ) ψ ¨ r ( 2 ) h 3 2 I 0 ( 3 ) ψ ¨ r ( 3 ) (10)

( 1 r + r ) ( h 1 2 N r ( 1 ) + M r ( 1 ) ) 1 r ( h 1 2 N θ ( 1 ) + M θ ( 1 ) ) Q r z ( 1 ) = h 1 2 I 0 ( 1 ) ( u ¨ 0 + h 2 2 ψ ¨ r ( 2 ) ) + ( I 2 ( 1 ) ψ ¨ r ( 1 ) + h 1 2 4 I 0 ( 1 ) ψ ¨ r ( 1 ) ) (11)

( 1 r + r ) ( h 2 2 N r ( 1 ) + M r ( 2 ) h 2 2 N r , r ( 3 ) ) 1 r h 2 2 ( N θ ( 1 ) N θ ( 3 ) ) Q r z ( 2 ) = h 2 2 I 0 ( 1 ) ( u ¨ 0 + h 1 2 ψ ¨ r ( 1 ) ) h 2 2 I 0 ( 3 ) ( u ¨ 0 h 3 2 ψ ¨ r ( 3 ) ) + ( h 2 2 4 I 0 ( 1 ) + I 2 ( 2 ) + h 2 2 4 I 0 ( 3 ) ) ψ ¨ r ( 2 ) (12)

( 1 r + r ) ( h 3 2 N r ( 3 ) + M r ( 3 ) ) Q r z ( 3 ) = h 3 2 I 0 ( 3 ) ( u ¨ 0 h 3 2 ψ ¨ r ( 3 ) ) + h 2 h 3 4 I 0 ( 3 ) ψ ¨ r ( 2 ) + I 2 ( 3 ) ψ ¨ r ( 3 ) (13)

k = 0 3 ( 1 r + r ) Q r ( k ) = q K w w C t w ˙ + ( I 0 ( 1 ) + I 0 ( 2 ) + I 0 ( 3 ) ) w ¨ (14)

where the stress resultants M, N, Q and the higher-order inertias are defined as:

{ N i ( k ) M i ( k ) } = h k 2 h k 2 σ i ( k ) { 1 ξ ( k ) } d ξ ( k ) , Q r ( k ) = h k 2 h k 2 τ r z ( k ) d ξ ( k ) , k = 1, 2, 3 i = r , θ (15)

I j ( k ) = h k 2 h k 2 ρ ( k ) ( r ) ξ ( k ) j d ξ ( k ) , I j ( k ) = [ 1 + γ ( k ) ( r r o ) η ( k ) ] I ¯ j ( k ) , k = 1,3 I ¯ j ( i ) = h i 2 h i 2 ρ ¯ ( i ) ξ ( i ) j d ξ ( i ) i = 1,2,3 j = 0,2 (16)

Based on Eq. (6), the stress resultants (Eq. (15)) may be rewritten in the following form:

{ N i ( 1 ) = A r i ( 1 ) ( u 0, r + h 1 2 ψ r , r ( 1 ) + h 2 2 ψ r , r ( 2 ) ) + 1 r A i θ ( 1 ) ( u 0 + h 1 2 ψ r ( 1 ) + h 2 2 ψ r ( 2 ) ) , N i ( 2 ) = A r i ( 2 ) u 0, r + 1 r A i θ ( 2 ) u 0 , N i ( 3 ) = A r i ( 3 ) ( u 0, r h 3 2 ψ r , r ( 3 ) h 2 2 ψ r , r ( 2 ) ) + 1 r A i θ ( 3 ) ( u 0 h 3 2 ψ r ( 3 ) h 2 2 ψ r ( 2 ) ) , i = r , θ (17)

M i ( k ) = D r i ( k ) ψ r , r ( k ) + 1 r D i θ ( k ) ψ r ( k ) , k = 1, 2, 3, i = r , θ (18)

Q r ( k ) = A r z ( k ) ( ψ r ( k ) + w , r ) , k = 1,2,3 (19)

where

{ A i i ( k ) D i i ( k ) } = h k 2 h k 2 E i ( k ) ( r ) 1 v r θ ( k ) v θ r ( k ) { 1 ξ ( k ) 2 } d ξ ( k ) , { A r θ ( k ) D r θ ( k ) } = h k 2 h k 2 v θ r ( k ) E r ( k ) ( r ) 1 v r θ ( k ) v θ r ( k ) { 1 ξ ( k ) 2 } d ξ ( k ) , A r z ( k ) = h k 2 h k 2 G r z ( k ) ( r ) d ξ ( k ) , k = 1,2,3 i = r , θ (20)

{ A i j ( k ) D i j ( k ) } = { A ¯ i j ( k ) D ¯ i j ( k ) } [ 1 + α ( k ) ( r r o ) β ( k ) ] , A r z ( k ) = A ¯ r z ( k ) [ 1 + α ( k ) ( r r o ) β ( k ) ] , i , j = r , θ k = 1,3 (21)

{ A ¯ i i ( k ) D ¯ i i ( k ) } = h k 2 h k 2 E ¯ i ( k ) 1 v r θ ( k ) v θ r ( k ) { 1 ξ ( k ) 2 } d ξ ( k ) , { A ¯ r θ ( k ) D ¯ r θ ( k ) } = h k 2 h k 2 v θ r ( k ) E ¯ r ( k ) 1 v r θ ( k ) v θ r ( k ) { 1 ξ ( k ) 2 } d ξ ( k ) , A ¯ r z ( k ) = h k 2 h k 2 G ¯ r z ( k ) d ξ ( k ) , i = r , θ k = 1, 3 (22)

Based on Eqs. (15) to (22), the governing equations (10) to (14) may be rewritten as:

[ A ¯ r r ( 1 ) + A r r ( 2 ) + A ¯ r r ( 3 ) + A ¯ r r ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) + A ¯ r r ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] ( u 0, r r + u 0, r r ) [ A ¯ θ θ ( 1 ) + A θ θ ( 2 ) + A ¯ θ θ ( 3 ) + A ¯ θ θ ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) + A ¯ θ θ ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] u 0 r 2 + h 1 2 [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] [ A ¯ r r ( 1 ) ( ψ r , r r ( 1 ) + ψ r , r ( 1 ) r ) A ¯ θ θ ( 1 ) ψ r ( 1 ) r 2 ] + h 2 2 [ A ¯ r r ( 1 ) A ¯ r r ( 3 ) + A ¯ r r ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) A ¯ r r ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] ( ψ r , r r ( 2 ) + ψ r , r ( 2 ) r ) h 2 2 [ A ¯ θ θ ( 1 ) A ¯ θ θ ( 3 ) + A ¯ θ θ ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) A ¯ θ θ ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] ψ r ( 2 ) r 2 h 3 2 [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] [ A ¯ r r ( 3 ) ( ψ r , r r ( 3 ) + ψ r , r ( 3 ) r ) A ¯ θ θ ( 3 ) ψ r ( 3 ) r 2 ] + α ( 1 ) β ( 1 ) ( r r o ) β ( 1 ) 1 [ ( A ¯ r r ( 1 ) r + A ¯ r θ ( 1 ) r ) ( u 0 + h 1 2 ψ r ( 1 ) + h 2 2 ψ r ( 2 ) ) ] + α ( 3 ) β ( 3 ) ( r r o ) β ( 3 ) 1 [ ( A ¯ r r ( 3 ) r + A ¯ r θ ( 3 ) r ) ( u 0 h 3 2 ψ r ( 3 ) h 2 2 ψ r ( 2 ) ) ] = ( I ¯ 0 ( 1 ) + I ¯ 0 ( 2 ) + I ¯ 0 ( 3 ) ) u ¨ 0 + γ ( 1 ) ( r r o ) η ( 1 ) I ¯ 0 ( 1 ) ( u ¨ 0 + h 1 2 ψ ¨ r ( 1 ) + h 2 2 ψ ¨ r ( 2 ) ) + h 1 2 I ¯ 0 ( 1 ) ψ ¨ r ( 1 ) + γ ( 3 ) ( r r o ) η ( 3 ) I ¯ 0 ( 3 ) ( u ¨ 0 h 2 2 ψ ¨ r ( 2 ) h 3 2 ψ ¨ r ( 3 ) ) + h 2 2 ( I ¯ 0 ( 1 ) I ¯ 0 ( 3 ) ) ψ ¨ r ( 2 ) h 3 2 I ¯ 0 ( 3 ) ψ ¨ r ( 3 ) (23)

h 1 2 [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] [ A ¯ r r ( 1 ) ( u 0, r r + u 0, r r ) A ¯ θ θ ( 1 ) u 0 r 2 ] + ( 1 + α ( 1 ) ( r r o ) β ( 1 ) ) [ ( h 1 2 4 A ¯ r r ( 1 ) + D ¯ r r ( 1 ) ) ( ψ r , r r ( 1 ) + ψ r , r ( 1 ) r ) ( h 1 2 4 A ¯ θ θ ( 1 ) + D ¯ θ θ ( 1 ) ) ψ r ( 1 ) r 2 ] + h 1 h 2 4 [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] [ A ¯ r r ( 1 ) ( ψ r , r r ( 2 ) + ψ r , r ( 2 ) r ) A ¯ θ θ ( 1 ) ψ r ( 2 ) r 2 ] A ¯ r z ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] ( ψ r ( 1 ) + w , r ) + α ( 1 ) β ( 1 ) ( r r o ) β ( 1 ) 1 ( D ¯ r r ( 1 ) ψ r , r ( 1 ) + D ¯ r θ ( 1 ) r ψ r ( 1 ) ) + h 1 2 α ( 1 ) β ( 1 ) ( r r o ) β ( 1 ) 1 [ A ¯ r r ( 1 ) ( u 0, r + h 1 2 ψ r , r ( 1 ) + h 2 2 ψ r , r ( 2 ) ) + A ¯ r θ ( 1 ) r ( u 0 + h 1 2 ψ r ( 1 ) + h 2 2 ψ r ( 2 ) ) ] = h 1 2 I ¯ 0 ( 1 ) ( u ¨ 0 + h 1 2 ψ ¨ r ( 1 ) + h 2 2 ψ ¨ r ( 2 ) ) + h 1 2 I ¯ 0 ( 1 ) γ ( 1 ) ( r r o ) η ( 1 ) ( u ¨ 0 + h 1 2 ψ ¨ r ( 1 ) + h 2 2 ψ ¨ r ( 2 ) ) + I ¯ 2 ( 1 ) ψ ¨ r ( 1 ) + I ¯ 2 ( 1 ) γ ( 1 ) ( r r o ) η ( 1 ) ψ ¨ r ( 1 ) (24)

h 2 2 [ A ¯ r r ( 1 ) A ¯ r r ( 3 ) + A ¯ r r ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) A ¯ r r ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] ( u 0, r r + u 0, r r ) h 2 2 [ A ¯ θ θ ( 1 ) A ¯ θ θ ( 3 ) + A ¯ θ θ ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) A ¯ θ θ ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] u 0 r 2 + h 1 h 2 4 A ¯ r r ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] ( ψ r , r r ( 1 ) + ψ r , r ( 1 ) r ) h 1 h 2 4 A ¯ θ θ ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] ψ r ( 1 ) r 2 + [ h 2 2 4 A ¯ r r ( 1 ) + D r r ( 2 ) + h 2 2 4 A ¯ r r ( 3 ) + h 2 2 4 A ¯ r r ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) + h 2 2 4 A ¯ r r ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] ( ψ r , r r ( 2 ) + ψ r , r ( 2 ) r ) [ h 2 2 4 A ¯ θ θ ( 1 ) + D θ θ ( 2 ) + h 2 2 4 A ¯ θ θ ( 3 ) + h 2 2 4 A ¯ θ θ ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) + h 2 2 4 A ¯ θ θ ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] ψ r ( 2 ) r 2 + h 2 h 3 4 A ¯ r r ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ( ψ r , r r ( 3 ) + ψ r , r ( 3 ) r ) h 2 h 3 4 A ¯ θ θ ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ψ r ( 3 ) r 2 A r z ( 2 ) ( ψ r ( 2 ) + w , r ) + h 2 2 α ( 1 ) β ( 1 ) ( r r o ) β ( 1 ) 1 [ A ¯ r r ( 1 ) ( u 0, r + h 1 2 ψ r , r ( 1 ) + h 2 2 ψ r , r ( 2 ) ) + A ¯ r θ ( 1 ) r ( u 0 + h 1 2 ψ r ( 1 ) + h 2 2 ψ r ( 2 ) ) ] h 2 2 α ( 3 ) β ( 3 ) ( r r o ) β ( 3 ) 1 [ A ¯ r r ( 3 ) ( u 0, r h 3 2 ψ r , r ( 3 ) h 2 2 ψ r , r ( 2 ) ) + A ¯ r θ ( 3 ) r ( u 0 h 3 2 ψ r ( 3 ) h 2 2 ψ r ( 2 ) ) ] = h 2 2 I ¯ 0 ( 1 ) ( u ¨ 0 + h 1 2 ψ ¨ r ( 1 ) + h 2 2 ψ ¨ r ( 2 ) ) + I ¯ 2 ( 2 ) ψ ¨ r ( 2 ) h 2 2 I ¯ 0 ( 3 ) ( u ¨ 0 h 2 2 ψ ¨ r ( 2 ) h 3 2 ψ ¨ r ( 3 ) ) + h 2 2 I ¯ 0 ( 1 ) γ ( 1 ) ( r r o ) η ( 1 ) ( u ¨ 0 + h 1 2 ψ ¨ r ( 1 ) + h 2 2 ψ ¨ r ( 2 ) ) h 2 2 I ¯ 0 ( 3 ) γ ( 3 ) ( r r o ) η ( 3 ) ( u ¨ 0 h 2 2 ψ ¨ r ( 2 ) h 3 2 ψ ¨ r ( 3 ) ) (25)

h 3 2 A ¯ r r ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ( u 0, r r + u 0, r r ) + h 3 2 A ¯ θ θ ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] u 0 r 2 + h 2 h 3 4 [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] [ A ¯ r r ( 3 ) ( ψ r , r r ( 2 ) + ψ r , r ( 2 ) r ) A ¯ θ θ ( 3 ) ψ r ( 2 ) r 2 ] + ( D ¯ r r ( 3 ) + h 3 2 4 A ¯ r r ( 3 ) ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ( ψ r , r r ( 3 ) + ψ r , r ( 3 ) r ) ( D ¯ θ θ ( 3 ) + h 3 2 4 A ¯ θ θ ( 3 ) ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ψ r ( 3 ) r 2 A ¯ r z ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ( ψ r ( 3 ) + w , r ) + α ( 3 ) β ( 3 ) ( r r o ) β ( 3 ) 1 ( D ¯ r r ( 3 ) ψ r , r ( 3 ) + D ¯ r θ ( 3 ) r ψ r ( 3 ) ) h 3 2 α ( 3 ) β ( 3 ) ( r r o ) β ( 3 ) 1 [ A ¯ r r ( 3 ) ( u 0, r h 3 2 ψ r , r ( 3 ) h 2 2 ψ r , r ( 2 ) ) + A ¯ r θ ( 3 ) r ( u 0 h 3 2 ψ r ( 3 ) h 2 2 ψ r ( 2 ) ) ] = h 3 2 I ¯ 0 ( 3 ) ( u ¨ 0 h 2 2 ψ ¨ r ( 2 ) h 3 2 ψ ¨ r ( 3 ) ) h 3 2 I ¯ 0 ( 3 ) γ ( 3 ) ( r r o ) η ( 3 ) ( u ¨ 0 h 2 2 ψ ¨ r ( 2 ) h 3 2 ψ ¨ r ( 3 ) ) + I ¯ 2 ( 3 ) ψ ¨ r ( 3 ) + I ¯ 2 ( 3 ) γ ( 3 ) ( r r o ) η ( 3 ) ψ ¨ r ( 3 ) (26)

[ A ¯ r z ( 1 ) + A r z ( 2 ) + A ¯ r z ( 3 ) + A ¯ r z ( 1 ) α ( 1 ) ( r r o ) β ( 1 ) + A ¯ r z ( 3 ) α ( 3 ) ( r r o ) β ( 3 ) ] ( w , r r + w , r r ) + A ¯ r z ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] ( ψ r , r ( 1 ) + ψ r ( 1 ) r ) + A 44 ( 2 ) ( ψ r , r ( 2 ) + ψ r ( 2 ) r ) + A ¯ r z ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ( ψ r , r ( 3 ) + ψ r ( 3 ) r ) + A ¯ r z ( 1 ) α ( 1 ) β ( 1 ) ( r r o ) β ( 1 ) 1 ( ψ r ( 1 ) + w , r ) + A ¯ r z ( 3 ) α ( 3 ) β ( 3 ) ( r r o ) β ( 3 ) 1 ( ψ r ( 3 ) + w , r ) = q + K w w + C t w ˙ + ( I ¯ 0 ( 1 ) + I ¯ 0 ( 2 ) + I ¯ 0 ( 3 ) ) w ¨ + γ ( 1 ) ( r r o ) η ( 1 ) I ¯ 0 ( 1 ) w ¨ + γ ( 3 ) ( r r o ) η ( 3 ) I ¯ 0 ( 3 ) w ¨ (27)

For solution of the governing equations (23) to (27), various combinations of the edge conditions may be employed.

  • I. Clamped immovable edge:

{ u 0 = 0 ψ r ( 1 ) = 0 ψ r ( 2 ) = 0 ψ r ( 3 ) = 0 w = 0 (28)

  • II. Simply-supported movable edge:

{ N r ( 1 ) + N r ( 2 ) + N r ( 3 ) = 0 h 1 2 N r ( 1 ) + M r ( 1 ) = 0 h 2 2 N r ( 1 ) + M r ( 2 ) h 2 2 N r ( 3 ) = 0 h 3 2 N r ( 3 + M r ( 3 ) = 0 w = 0 (29)

  • III. Free edge:

{ N r ( 1 ) + N r ( 2 ) + N r ( 3 ) = 0 h 1 2 N r ( 1 ) + M r ( 1 ) = 0 h 2 2 N r ( 1 ) + M r ( 2 ) h 2 2 N r ( 3 ) = 0 h 3 2 N r ( 3 + M r ( 3 ) = 0 Q r ( 1 ) + Q r ( 2 ) + Q r ( 3 ) = 0 (30)

Also, initial conditions of the sandwich plate are:

u ( r ,0 ) = u ˙ ( r ,0 ) = 0, ψ r ( 1 ) ( r ,0 ) = ψ ˙ r ( 1 ) ( r ,0 ) = 0, ψ r ( 2 ) ( r ,0 ) = ψ ˙ r ( 2 ) ( r ,0 ) = 0, ψ r ( 3 ) ( r ,0 ) = ψ ˙ r ( 3 ) ( r ,0 ) = 0, w ( r ,0 ) = w ˙ ( r ,0 ) = 0. (31)

4 SEMI-ANALYTICAL SOLUTION FOR DYNAMIC ANALYSIS OF FG POLAR ORTHOTROPIC SANDWICH PLATE

For solution of the relatively complicated five systems of second order coupled partial differential equations, the semi-analytical solution is developed based on the finite Taylor’s transform and the fourth-order Runge-Kutta procedure.

Based on Taylor’s expansion, the unknown displacement functions can be expressed by the following power series:

{ u 0 ( r , t = y Δ t ) u ˙ 0 ( r , t = y Δ t ) u ¨ 0 ( r , t = y Δ t ) } = x = 0 { U x , y U ¯ x , y U ˜ x , y } ( r r o ) x , { ψ r ( 1 ) ( r , t = y Δ t ) ψ ˙ r ( 1 ) ( r , t = y Δ t ) ψ ¨ r ( 1 ) ( r , t = y Δ t ) } = x = 0 { Φ x , y ( 1 ) Φ ¯ x , y ( 1 ) Φ ˜ x , y ( 1 ) } ( r r o ) x , { ψ r ( 2 ) ( r , t = y Δ t ) ψ ˙ r ( 2 ) ( r , t = y Δ t ) ψ ¨ r ( 2 ) ( r , t = y Δ t ) } = x = 0 { Φ x , y ( 2 ) Φ ¯ x , y ( 2 ) Φ ˜ x , y ( 2 ) } ( r r o ) x , { ψ r ( 3 ) ( r , t = y Δ t ) ψ ˙ r ( 3 ) ( r , t = y Δ t ) ψ ¨ r ( 3 ) ( r , t = y Δ t ) } = x = 0 { Φ x , y ( 3 ) Φ ¯ x , y ( 3 ) Φ ˜ x , y ( 3 ) } ( r r o ) x , { w ( r , t = y Δ t ) w ˙ ( r , t = y Δ t ) w ¨ ( r , t = y Δ t ) } = x = 0 { W x , y W ¯ x , y W ˜ x , y } ( r r o ) x (32)

Where ΔT is the time step, y is the time step counter and Ux,y, U¯x,y, U˜x,y, Φx,y(1), Φ¯x,y(1), Φ˜x,y(1), Φx,y(2), Φ¯x,y(2), Φ˜x,y(2), Φx,y(3), Φ¯x,y(3), Φ˜x,y(3), Wx,y, W¯x,y and W˜x,y are the coefficients of series in each time step.

On the other hands, Taylor transform of the functions 1/r and 1/r 2 may be expressed by the following power series whose center is located at r=ro.

1 r = s = 0 ( 1 r o ) s + 1 ( r r o ) s , 1 r 2 = s = 0 ( s + 1 ) ( 1 r o ) s + 2 ( r r o ) s (33)

In practical applications, the transformed form of functions must be expressed by means of finite series.

The transformed form of the governing equations may be obtained by substituting Eqs. (32) and (33) into the governing Eqs. (23) to (27) and performing some manipulations.

x = 0 X [ ( A ¯ r r ( 1 ) + A r r ( 2 ) + A ¯ r r ( 3 ) ) [ ( x + 2 ) ( x + 1 ) U x + 2, y s = 0 x χ s + 1 ( x s + 1 ) U x s + 1, y ] A ¯ r r ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) U x β ( 1 ) + 2, y s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) U x β ( 1 ) s + 1, y ] + A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) U x β ( 3 ) + 2, y s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) U x β ( 3 ) s + 1, y ] + ( A ¯ θ θ ( 1 ) + A θ θ ( 2 ) + A ¯ θ θ ( 3 ) ) s = 0 x ( s + 1 ) χ s + 2 U x s , y A ¯ θ θ ( 1 ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 U x β ( 1 ) s , y A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 3 ) ( s + 1 ) χ s + 2 U x β ( 3 ) s , y + h 1 2 A ¯ r r ( 1 ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 1 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 1 ) ] + h 1 2 A ¯ r r ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) Φ x β ( 1 ) + 2, y ( 1 ) s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) Φ x β ( 1 ) s + 1, y ( 1 ) ] h 1 2 A ¯ θ θ ( 1 ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 1 ) h 1 2 A ¯ θ θ ( 1 ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 Φ x β ( 1 ) s , y ( 1 ) + h 2 2 ( A ¯ r r ( 1 ) A ¯ r r ( 3 ) ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 2 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 2 ) ] + h 2 2 A ¯ r r ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) Φ x β ( 1 ) + 2, y ( 2 ) s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) Φ x β ( 1 ) s + 1, y ( 2 ) ] h 2 2 A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) Φ x β ( 3 ) + 2, y ( 2 ) s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) Φ x β ( 3 ) s + 1, y ( 2 ) ] h 2 2 ( A ¯ θ θ ( 1 ) A ¯ θ θ ( 3 ) ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 2 ) h 2 2 A ¯ θ θ ( 1 ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 Φ x β ( 1 ) s , y ( 2 ) + h 2 2 A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 3 ) ( s + 1 ) χ s + 2 Φ x β ( 3 ) s , y ( 2 ) h 3 2 A ¯ r r ( 3 ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 3 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 3 ) ] h 3 2 A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) Φ x β ( 3 ) + 2, y ( 3 ) s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) Φ x β ( 3 ) s + 1, y ( 3 ) ] + h 3 2 A ¯ θ θ ( 3 ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 3 ) + α ( 1 ) β ( 1 ) A ¯ r r ( 1 ) ( x β ( 1 ) + 2 ) ( U x β ( 1 ) + 2, y + h 1 2 Φ x β ( 1 ) + 2 ( 1 ) + h 2 2 Φ x β ( 1 ) + 2, y ( 2 ) ) + h 3 2 A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 3 ) ( s + 1 ) χ s + 2 Φ x β ( 3 ) s , y ( 3 ) α ( 1 ) β ( 1 ) A ¯ r θ ( 1 ) s = 0 x β ( 1 ) + 1 χ s + 1 ( U x s β ( 1 ) + 1 + h 1 2 Φ x s β ( 1 ) + 1, y ( 1 ) + h 2 2 Φ x s β ( 1 ) + 1, y ( 2 ) ) + α ( 3 ) β ( 3 ) A ¯ r r ( 3 ) ( x β ( 3 ) + 2 ) ( U x β ( 3 ) + 2, y h 3 2 Φ x β ( 3 ) + 2, y ( 3 ) h 2 2 Φ x β ( 3 ) + 2, y ( 2 ) ) α ( 3 ) β ( 3 ) A ¯ r θ ( 3 ) s = 0 x β ( 3 ) + 1 χ s + 1 ( U x s β ( 3 ) + 1, y h 3 2 Φ x s β ( 3 ) + 1, y ( 3 ) h 2 2 Φ x s β ( 3 ) + 1, y ( 2 ) ) ( I ¯ 0 ( 1 ) + I ¯ 0 ( 2 ) + I ¯ 0 ( 3 ) ) U ˜ x , y γ ( 1 ) I ¯ 0 ( 1 ) U ˜ x η ( 1 ) , y γ ( 2 ) I ¯ 0 ( 2 ) U ˜ x η ( 2 ) , y γ ( 3 ) I ¯ 0 ( 3 ) U ˜ x η ( 3 ) , y h 1 2 I ¯ 0 ( 1 ) Φ ˜ x , y ( 1 ) h 2 2 ( I ¯ 0 ( 1 ) I ¯ 0 ( 3 ) ) Φ ˜ x , y ( 2 ) + h 3 2 I ¯ 0 ( 3 ) Φ ˜ x , y ( 3 ) γ ( 1 ) I ¯ 0 ( 1 ) ( h 1 2 Φ ˜ x η ( 1 ) , y ( 1 ) + h 2 2 Φ ˜ x η ( 1 ) , y ( 2 ) ) + γ ( 3 ) I ¯ 0 ( 3 ) ( h 2 2 Φ ˜ x η ( 3 ) , y ( 2 ) + h 3 2 Φ ˜ x η ( 3 ) , y ( 3 ) ) ] ( r r o ) x = 0 (34)

x = 0 X { h 1 2 A ¯ r r ( 1 ) [ ( x + 2 ) ( x + 1 ) U x + 2, y s = 0 x χ s + 1 ( x s + 1 ) U x s + 1, y ] h 1 2 A ¯ θ θ ( 1 ) s = 0 x ( s + 1 ) χ s + 2 U x s , y + h 1 2 A ¯ r r ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) U x β ( 1 ) + 2, y s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) U x β ( 1 ) s + 1, y ] h 1 2 A ¯ θ θ ( 1 ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 U x β ( 1 ) s , y + ( h 1 2 4 A ¯ r r ( 1 ) + D ¯ r r ( 1 ) ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 1 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 1 ) ] + ( h 1 2 4 A ¯ r r ( 1 ) + D ¯ r r ( 1 ) ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) Φ x β ( 1 ) + 2, y ( 1 ) s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) Φ x β ( 1 ) s + 1, y ( 1 ) ] ( h 1 2 4 A ¯ θ θ ( 1 ) + D ¯ θ θ ( 1 ) ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 1 ) ( h 1 2 4 A ¯ θ θ ( 1 ) + D ¯ θ θ ( 1 ) ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 Φ x β ( 1 ) s , y ( 1 ) + h 1 h 2 4 A ¯ r r ( 1 ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 2 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 2 ) ] α ( 1 ) β ( 1 ) D ¯ r θ ( 1 ) s = 0 x β ( 1 ) + 1 χ s + 1 Φ x s β ( 1 ) + 1, y ( 1 ) + h 1 h 2 4 α ( 1 ) A ¯ r r ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) Φ x β ( 1 ) + 2, y ( 2 ) s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) Φ x β ( 1 ) s + 1, y ( 2 ) ] h 1 h 2 4 A ¯ θ θ ( 1 ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 2 ) h 1 h 2 4 α ( 1 ) A ¯ θ θ ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 Φ x β ( 1 ) s , y ( 2 ) A ¯ r z ( 1 ) [ Φ x , y ( 1 ) + ( x + 1 ) W x + 1, y ] A ¯ r z ( 1 ) α ( 1 ) [ Φ x β ( 1 ) , y ( 1 ) + ( x β ( 1 ) + 1 ) W x β ( 1 ) + 1, y ] + α ( 1 ) β ( 1 ) D ¯ r r ( 1 ) ( x β ( 1 ) + 2 ) Φ x β ( 1 ) + 2, y ( 1 ) + h 1 2 α ( 1 ) β ( 1 ) A ¯ r r ( 1 ) ( x β ( 1 ) + 2 ) ( U x β ( 1 ) + 2, y + h 1 2 Φ x β ( 1 ) + 2, y ( 1 ) + h 2 2 Φ x β ( 1 ) + 2, y ( 2 ) ) h 1 2 α ( 1 ) β ( 1 ) A ¯ r θ ( 1 ) s = 0 x β ( 1 ) + 1 χ s + 1 ( U x s β ( 1 ) + 1, y + h 1 2 Φ x s β ( 1 ) + 1, y ( 1 ) + h 2 2 Φ x s β ( 1 ) + 1, y ( 2 ) ) h 1 2 I ¯ 0 ( 1 ) ( U ˜ x , y + h 1 2 Φ ˜ x , y ( 1 ) + h 2 2 Φ ˜ x , y ( 2 ) ) h 1 2 I ¯ 0 ( 1 ) γ ( 1 ) ( U ˜ x η ( 1 ) , y + h 1 2 Φ ˜ x η ( 1 ) , y ( 1 ) + h 2 2 Φ ˜ x η ( 1 ) , y ( 2 ) ) I ¯ 2 ( 1 ) Φ ˜ x , y ( 1 ) I ¯ 2 ( 1 ) γ ( 1 ) Φ ˜ x η ( 1 ) , y ( 1 ) } ( r r 0 ) x = 0 (35)

x = 0 X { h 2 2 ( A ¯ r r ( 1 ) A ¯ r r ( 3 ) ) [ ( x + 2 ) ( x + 1 ) U x + 2, y s = 0 x χ s + 1 ( x s + 1 ) U x s + 1, y ] + h 2 2 A ¯ r r ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) U x β ( 1 ) + 2, y s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) U x β ( 1 ) s + 1, y ] h 2 2 A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) U x β ( 3 ) + 2, y s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) ) U x β ( 3 ) s + 1, y ] h 2 2 ( A ¯ θ θ ( 1 ) A ¯ θ θ ( 3 ) ) s = 0 x ( s + 1 ) χ s + 2 U x s , y h 2 2 A ¯ θ θ ( 1 ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 U x β ( 1 ) s , y + h 2 2 A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 U x β ( 3 ) s , y + h 1 h 2 4 A ¯ r r ( 1 ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 1 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 1 ) ] + h 1 h 2 4 A ¯ r r ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) Φ x β ( 1 ) + 2, y ( 1 ) s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) Φ x β ( 1 ) s + 1, y ( 1 ) ] h 1 h 2 4 A ¯ θ θ ( 1 ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 1 ) h 1 h 2 4 A ¯ θ θ ( 1 ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 Φ x β ( 1 ) s , y ( 1 ) + ( h 2 2 4 A ¯ r r ( 1 ) + D r r ( 2 ) + h 2 2 4 A ¯ r r ( 3 ) ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 2 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 2 ) ] + h 2 2 4 A ¯ r r ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) Φ x β ( 1 ) + 2, y ( 2 ) s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) Φ x β ( 1 ) s + 1, y ( 2 ) ] + h 2 2 4 A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) Φ x β ( 3 ) + 2, y ( 2 ) s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) Φ x β ( 3 ) s + 1, y ( 2 ) ] ( h 2 2 4 A ¯ θ θ ( 1 ) + D θ θ ( 2 ) + h 2 2 4 A ¯ θ θ ( 3 ) ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 2 ) h 2 2 4 A ¯ θ θ ( 1 ) α ( 1 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 Φ x β ( 1 ) s , y ( 2 ) h 2 2 4 A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 3 ) ( s + 1 ) χ s + 2 Φ x β ( 3 ) s , y ( 2 ) + h 2 h 3 4 A ¯ r r ( 3 ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 3 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 3 ) ] + h 2 h 3 4 A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) Φ x β ( 3 ) + 2, y ( 3 ) s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) Φ x β ( 3 ) s + 1, y ( 3 ) ] h 2 h 3 4 A ¯ θ θ ( 3 ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 3 ) h 2 h 3 4 A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 3 ) ( s + 1 ) χ s + 2 Φ x β ( 3 ) s , y ( 3 ) A r z ( 2 ) [ Φ x ( 2 ) + ( x + 1 ) W x + 1, y ] + h 2 2 α ( 1 ) β ( 1 ) A ¯ r r ( 1 ) ( x β ( 1 ) + 2 ) ( U x β ( 1 ) + 2, y + h 1 2 Φ x β ( 1 ) + 2, y ( 1 ) + h 2 2 Φ x β ( 1 ) + 2, y ( 2 ) ) h 2 2 α ( 1 ) β ( 1 ) A ¯ r θ ( 1 ) s = 0 χ s + 1 ( U x s β ( 1 ) + 1, y + h 1 2 Φ x s β ( 1 ) + 1, y ( 1 ) + h 2 2 Φ x s β ( 1 ) + 1, y ( 2 ) ) h 2 2 α ( 3 ) β ( 3 ) A ¯ r r ( 3 ) ( x β ( 3 ) + 2 ) ( U x β ( 3 ) + 2, y h 3 2 Φ x β ( 3 ) + 2, y ( 3 ) h 2 2 Φ x β ( 3 ) + 2, y ( 2 ) ) + h 2 2 α ( 3 ) β ( 3 ) A ¯ r θ ( 3 ) s = 0 χ s + 1 ( U x s β ( 3 ) + 1, y h 3 2 Φ x s β ( 3 ) + 1, y ( 3 ) h 2 2 Φ x s β ( 3 ) + 1, y ( 2 ) ) h 2 2 I ¯ 0 ( 1 ) ( U ˜ x , y + h 1 2 Φ ˜ x , y ( 1 ) + h 2 2 Φ ˜ x , y ( 2 ) ) I ¯ 2 ( 2 ) Φ ˜ x , y ( 2 ) + h 2 2 I ¯ 0 ( 3 ) ( U ˜ x , y h 2 2 Φ ˜ x , y ( 2 ) h 3 2 Φ ˜ x , y ( 3 ) ) h 2 2 I ¯ 0 ( 1 ) γ ( 1 ) ( U ˜ x η ( 1 ) , y + h 1 2 Φ ˜ x η ( 1 ) , y ( 1 ) + h 2 2 Φ ˜ x η ( 1 ) , y ( 2 ) ) I ¯ 2 ( 2 ) γ ( 2 ) Φ ˜ x η ( 2 ) , y ( 2 ) + h 2 2 I ¯ 0 ( 3 ) γ ( 3 ) ( U ˜ x η ( 3 ) , y h 2 2 Φ ˜ x η ( 3 ) , y ( 2 ) h 3 2 Φ ˜ x η ( 3 ) , y ( 3 ) ) } ( r r 0 ) x = 0 (36)

x = 0 X { h 3 2 A ¯ r r ( 3 ) [ ( x + 2 ) ( x + 1 ) U x + 2, y s = 0 x χ s + 1 ( x s + 1 ) U x s + 1, y ] h 3 2 A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) U x β ( 3 ) + 2, y s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) U x β ( 3 ) s + 1, y ] + h 3 2 A ¯ θ θ ( 3 ) s = 0 x ( s + 1 ) χ s + 2 U x s , y + h 3 2 A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 1 ) ( s + 1 ) χ s + 2 U x β ( 3 ) s , y + h 2 h 3 4 A ¯ r r ( 3 ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 2 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 2 ) ] + h 2 h 3 4 A ¯ r r ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) Φ x β ( 3 ) + 2, y ( 2 ) s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) Φ x β ( 3 ) s + 1, y ( 2 ) ] h 2 h 3 4 A ¯ θ θ ( 3 ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 2 ) h 2 h 3 4 A ¯ θ θ ( 3 ) α ( 3 ) s = 0 x β ( 3 ) ( s + 1 ) χ s + 2 Φ x β ( 3 ) s , y ( 2 ) + ( D ¯ r r ( 3 ) + h 3 2 4 A ¯ r r ( 3 ) ) [ ( x + 2 ) ( x + 1 ) Φ x + 2, y ( 3 ) s = 0 x χ s + 1 ( x s + 1 ) Φ x s + 1, y ( 3 ) ] + ( D ¯ r r ( 3 ) + h 3 2 4 A ¯ r r ( 3 ) ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) Φ x β ( 3 ) + 2, y ( 3 ) s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) Φ x β ( 3 ) s + 1, y ( 3 ) ] ( D ¯ θ θ ( 3 ) + h 3 2 4 A ¯ θ θ ( 3 ) ) s = 0 x ( s + 1 ) χ s + 2 Φ x s , y ( 3 ) ( D ¯ θ θ ( 3 ) + h 3 2 4 A ¯ θ θ ( 3 ) ) α ( 3 ) s = 0 x β ( 3 ) ( s + 1 ) χ s + 2 Φ x β ( 3 ) s , y ( 3 ) A ¯ r z ( 3 ) ( Φ x ( 3 ) + ( x + 1 ) W x + 1 ) A ¯ r z ( 3 ) α ( 3 ) [ Φ x β ( 3 ) ( 3 ) + ( x β ( 3 ) + 1 ) W x β ( 3 ) + 1 ] + α ( 3 ) β ( 3 ) D ¯ r r ( 3 ) ( x β ( 3 ) + 2 ) Φ x β ( 3 ) + 2, y ( 3 ) α ( 3 ) β ( 3 ) D ¯ r θ ( 3 ) s = 0 χ s + 1 Φ x s β ( 3 ) + 1, y ( 3 ) h 3 2 α ( 3 ) β ( 3 ) A ¯ r r ( 3 ) ( x β ( 3 ) + 2 ) ( U x β ( 3 ) + 2, y h 3 2 Φ x β ( 3 ) + 2, y ( 3 ) h 2 2 Φ x β ( 3 ) + 2, y ( 2 ) ) + h 3 2 α ( 3 ) β ( 3 ) A ¯ r θ ( 3 ) s = 0 χ s + 1 ( U x s β ( 3 ) + 1, y h 3 2 Φ x s β ( 3 ) + 1, y ( 3 ) h 2 2 Φ x s β ( 3 ) + 1, y ( 2 ) ) I ¯ 2 ( 3 ) γ ( 3 ) Φ ˜ x η ( 3 ) , y ( 3 ) I ¯ 2 ( 3 ) Φ ˜ x , y ( 3 ) + h 3 2 I ¯ 0 ( 3 ) ( U ˜ x , y h 3 2 Φ ˜ x , y ( 3 ) h 2 2 Φ ˜ x , y ( 2 ) ) + h 3 2 I ¯ 0 ( 3 ) γ ( 3 ) ( U ˜ x , y h 3 2 Φ ˜ x η ( 3 ) , y ( 3 ) h 2 2 Φ ˜ x η ( 3 ) , y ( 2 ) ) } ( r r o ) x = 0 (37)

x = 0 X { ( A ¯ r z ( 1 ) + A r z ( 2 ) + A ¯ r z ( 3 ) ) [ ( x + 2 ) ( x + 1 ) W x + 2, y s = 0 x χ s + 1 ( x s + 1 ) W x s + 1, y ] + A ¯ r z ( 1 ) α ( 1 ) [ ( x β ( 1 ) + 2 ) ( x β ( 1 ) + 1 ) W x β ( 1 ) + 2, y s = 0 x β ( 1 ) χ s + 1 ( x β ( 1 ) s + 1 ) W x β ( 1 ) s + 1, y ] + A ¯ r z ( 3 ) α ( 3 ) [ ( x β ( 3 ) + 2 ) ( x β ( 3 ) + 1 ) W x β ( 3 ) + 2, y s = 0 x β ( 3 ) χ s + 1 ( x β ( 3 ) s + 1 ) W x β ( 3 ) s + 1, y ] + A ¯ r z ( 1 ) ( x + 1 ) Φ x + 1, y ( 1 ) + A ¯ r z ( 1 ) α ( 1 ) ( x β ( 1 ) + 1 ) Φ x β ( 1 ) + 1, y ( 1 ) A ¯ r z ( 1 ) s = 0 χ s + 1 Φ x s , y ( 1 ) A ¯ r z ( 1 ) α ( 1 ) s = 0 χ s + 1 Φ x s β ( 1 ) , y ( 1 ) + A 44 ( 2 ) [ ( x + 1 ) Φ x + 1, y ( 2 ) s = 0 χ s + 1 Φ x s , y ( 2 ) ] + A ¯ r z ( 3 ) ( x + 1 ) Φ x + 1, y ( 3 ) + A ¯ r z ( 3 ) α ( 3 ) ( x β ( 3 ) + 1 ) Φ x β ( 3 ) + 1, y ( 3 ) A ¯ r z ( 3 ) s = 0 χ s + 1 Φ x s , y ( 3 ) A ¯ r z ( 3 ) α ( 3 ) s = 0 χ s + 1 Φ x β ( 3 ) s , y ( 3 ) + A ¯ r z ( 1 ) α ( 1 ) β ( 1 ) Φ x β ( 1 ) + 1, y ( 1 ) + q δ ( x ) K w W x , y C t W ¯ x , y + A ¯ r z ( 1 ) α ( 1 ) β ( 1 ) ( x β ( 1 ) + 2 ) W x β ( 1 ) + 2, y + A ¯ r z ( 3 ) α ( 3 ) β ( 3 ) Φ x β ( 3 ) + 1, y ( 3 ) + A ¯ r z ( 3 ) α ( 3 ) β ( 3 ) ( x β ( 3 ) + 2 ) W x β ( 3 ) + 2, y ( I ¯ 0 ( 1 ) + I ¯ 0 ( 2 ) + I ¯ 0 ( 3 ) ) W ˜ x , y γ ( 1 ) I ¯ 0 ( 1 ) W ˜ x η ( 1 ) , y γ ( 2 ) I ¯ 0 ( 2 ) W ˜ x η ( 2 ) , y γ ( 3 ) I ¯ 0 ( 3 ) W ˜ x η ( 3 ) , y } ( r r o ) x = 0 (38)

Where δ(x) is Kronecker’s delta function.

The transformed forms of the various boundary conditions at the inner (r=r i ) and outer (r=r o ) edges can be obtained by substituting Eq. (32) into the Eqs. (28-30).

u = 0 x = 0 X U x , y ( r r o ) x = 0 ψ r ( 1 ) = 0 x = 0 X Φ x , y ( 1 ) ( r r o ) x = 0, ψ r ( 2 ) = 0 x = 0 X Φ x , y ( 2 ) ( r r o ) x = 0, ψ r ( 3 ) = 0 x = 0 X Φ x , y ( 3 ) ( r r o ) x = 0, w = 0 x = 0 X W x , y ( r r o ) x = 0 (39)

N r ( 1 ) + N r ( 2 ) + N r ( 3 ) = 0 x = 0 X { [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] [ A ¯ r r ( 1 ) ( x + 1 ) ( U x + 1, y + h 1 2 Φ x + 1, y ( 1 ) + h 2 2 Φ x + 1, y ( 2 ) ) A ¯ r θ ( 1 ) r ( U x , y + h 1 2 Φ x , y ( 1 ) + h 2 2 Φ x , y ( 2 ) ) ] + [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] [ A ¯ r r ( 3 ) ( x + 1 ) ( U x + 1, y h 2 2 Φ x + 1, y ( 2 ) h 3 2 Φ x + 1, y ( 3 ) ) A ¯ r θ ( 3 ) r ( U x , y h 2 2 Φ x , y ( 2 ) h 3 2 Φ x , y ( 3 ) ) ] A r r ( 2 ) ( x + 1 ) U x + 1, y A ¯ r θ ( 2 ) U x , y r } ( r r o ) x = 0 h 1 2 N r ( 1 ) + M r ( 1 ) = 0 x = 0 X { h 1 2 [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] [ A ¯ r r ( 1 ) ( x + 1 ) ( U x + 1, y + h 1 2 Φ x + 1, y ( 1 ) + h 2 2 Φ x + 1, y ( 2 ) ) A ¯ r θ ( 1 ) r ( U x , y + h 1 2 Φ x , y ( 1 ) + h 2 2 Φ x , y ( 2 ) ) ] + [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] ( D ¯ r r ( 1 ) ψ r , r ( 1 ) + 1 r D ¯ r θ ( 1 ) ψ r ( 1 ) ) } ( r r o ) x = 0 h 2 2 N r ( 1 ) + M r ( 2 ) h 2 2 N r ( 3 ) = 0 x = 0 X { h 2 2 [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] [ A ¯ r r ( 1 ) ( x + 1 ) ( U x + 1, y + h 1 2 Φ x + 1, y ( 1 ) + h 2 2 Φ x + 1, y ( 2 ) ) A ¯ r θ ( 1 ) r ( U x , y + h 1 2 Φ x , y ( 1 ) + h 2 2 Φ x , y ( 2 ) ) ] h 2 2 [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] [ A ¯ r r ( 3 ) ( x + 1 ) ( U x + 1, y h 2 2 Φ x + 1, y ( 2 ) h 3 2 Φ x + 1, y ( 3 ) ) A ¯ r θ ( 3 ) r ( U x , y h 2 2 Φ x , y ( 2 ) h 3 2 Φ x , y ( 3 ) ) ] D r r ( 2 ) ψ r , r ( 2 ) + 1 r D r θ ( 2 ) ψ r ( 2 ) } ( r r o ) x = 0 h 3 2 N r ( 3 + M r ( 3 ) = 0 x = 0 X { h 3 2 [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] [ A ¯ r r ( 3 ) ( x + 1 ) ( U x + 1, y h 2 2 Φ x + 1, y ( 2 ) h 3 2 Φ x + 1, y ( 3 ) ) A ¯ r θ ( 3 ) r ( U x , y h 2 2 Φ x , y ( 2 ) h 3 2 Φ x , y ( 3 ) ) ] + [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] ( D ¯ r r ( 3 ) ψ r , r ( 3 ) + 1 r D ¯ r θ ( 3 ) ψ r ( 3 ) ) } ( r r o ) x = 0 Q r ( 1 ) + Q r ( 2 ) + Q r ( 3 ) = 0 x = 0 X { A ¯ r z ( 1 ) [ 1 + α ( 1 ) ( r r o ) β ( 1 ) ] [ Φ x , y ( 1 ) + ( x + 1 ) W x + 1, y ] + A r z ( 2 ) [ Φ x , y ( 2 ) + ( x + 1 ) W x + 1, y ] + A ¯ r z ( 3 ) [ 1 + α ( 3 ) ( r r o ) β ( 3 ) ] [ Φ x , y ( 3 ) + ( x + 1 ) W x + 1, y ] } ( r r o ) x = 0 (40)

The resulting governing differential equations can be solved in terms of time by the fourth-order Runge-Kutta method, numerically.

The unknown displacement, velocity and acceleration parameters can be expressed as following vectors:

F ( r , t ) = u ( r , t ) , ψ r ( 1 ) ( r , t ) , ψ r ( 2 ) ( r , t ) , ψ r ( 3 ) ( r , t ) , w ( r , t ) F ˙ ( r , t ) = u ˙ ( r , t ) , ψ ˙ r ( 1 ) ( r , t ) , ψ ˙ r ( 2 ) ( r , t ) , ψ ˙ r ( 3 ) ( r , t ) , w ˙ ( r , t ) F ¨ ( r , t ) = u ¨ ( r , t ) , ψ ¨ r ( 1 ) ( r , t ) , ψ ¨ r ( 2 ) ( r , t ) , ψ ¨ r ( 3 ) ( r , t ) , w ¨ ( r , t ) (41)

Also, the transformed forms of the unknown parameters may be expressed as following vectors:

F x , y = U x , y , Φ x , y ( 1 ) , Φ x , y ( 2 ) , Φ x , y ( 3 ) , W x , y F ¯ x , y = U ¯ x , y , Φ ¯ x , y ( 1 ) , Φ ¯ x , y ( 2 ) , Φ ¯ x , y ( 3 ) , W ¯ x , y F ˜ x , y = U ˜ x , y , Φ ˜ x , y ( 1 ) , Φ ˜ x , y ( 2 ) , Φ ˜ x , y ( 3 ) , W ˜ x , y (42)

Based on the fourth-order Runge-Kutta method, when the acceleration and velocity parameters at t = y∆t are known as F¯x,y and F˜x,y the displacement parameters at t = y∆t and t = (y + 1) ∆t are known as F x,y and F x,y+1 , the acceleration and velocity parameters at t = (y + 1) ∆t can be found by using the following formulations:

F ˜ x , y + 1 = 4 Δ t 2 ( F x , y + 1 F x , y Δ t F ¯ x , y 1 4 Δ t 2 F ˜ x , y ) (43)

F ¯ x , y + 1 = F ¯ x , y + Δ t 2 ( F ˜ x , y F ˜ x , y + 1 ) (44)

By using Eqs. (31) and (32), the initial conditions can be obtained as following transformed forms:

U x ,0 = U ¯ x ,0 = 0, Φ x ,0 ( 1 ) = Φ ¯ x ,0 ( 1 ) = 0, Φ x ,0 ( 2 ) = Φ ¯ x ,0 ( 2 ) = 0, Φ x ,0 ( 3 ) = Φ ¯ x ,0 ( 3 ) = 0, W x ,0 = W ¯ x ,0 = 0, (45)

The unknown functions can be obtained by using the transformed forms of the governing equations Eqs. (34)-(38), boundary conditions Eqs. (39, 40) and initial conditions Eq. (45) and the relations between the acceleration, velocity and displacement parameters Eqs. (43) and (44). From the Eqs. (39) and (40), five equations for each edge must be applied (based on the edge condition at the inner and outer edges).

The unknown functions may be obtained according to the following steps.

  1. Initial values of the acceleration parameters

    (F˜x,0)

    are determined based on the Eqs. (34)-(38) and Eq. (45).
  2. From Eqs. (42) and (43), the acceleration and velocity parameters at the end of each time step

    (F˜x,yandF¯x,y, y=1, 2, 3, ...)

    may be expressed based on the obtained parameters at the previous step

    (Fx,y1, F˜x,y1andF¯x,y1, y=1, 2, 3, ...)

    and unknown displacement parameters at the end of each time step (F x,y , y = 1,2,3...).
  3. The displacements parameters at the end of each time step (F x,y ) are determined based on the obtained results of the acceleration and velocity parameters from step (II), Eqs. (34)-(38) and the edge conditions (Eqs. (39, 40)) corresponding to the end of each time step.

  4. Repeating steps (II) and (III) till the final time instant is reached(t = y ∆T).

5 RESULTS AND DISCUSSIONS

In this section, effects of the viscoelastic foundation on the transient response of the functionally graded polar orthotropic sandwich annular plates are investigated. The presented results cover various dynamic loads, edge conditions and foundation parameters. Dynamic responses of the sandwich structures are extracted for various dynamic loads as q=q^ f(t), in which q^ is the maximum amplitude and f(t) is an arbitrary unit time-dependency function.

In the presented results, the following dynamic load cases are examined (as shown in Figure 2):

Figure 2:
The various kinds of dynamic loads in numerical examples.

  • Case I) Step loading: f(t) = 1

  • Case II) Pulse loading:

    f(t)={1 for 0 tτ0 for τ t

  • Case III) Sinusoidal loading: f(t) = sin(ωt)

  • Case IV) N-shaped pulse loading:

    f(t)={12tτ for 0 t τ0 for τ t

  • Case V) Sinusoidal loading acting in the same direction:

    f(t)=|sin(π τt)|

  • Case VI) A triangular-shaped pulse loading:

    f(t)={1tτ for 0 t τ0 for τ t

  • Case VII) Exponential loading: f(t) = e-ζ t

  • Case VIII) Repeated rectangular-shaped loading:

f ( t ) = { 1 for 2 i τ t ( 2 i + 1 ) τ 0 for ( 2 i + 1 ) τ t 2 ( i + 1 ) τ i = 0,1,2,...

  • Case IX) Repeated triangular-shaped loading:

f ( t ) = t i τ τ for 2 i τ t ( 2 i + 1 ) τ i = 0,1,2,...

Sandwich plates under dynamic loads with the maximum amplitude of q^ =1MPa are considered throughout the examples. Also, the following geometric information and material properties for the homogenous core are chosen.

r i = 0.2, r o = 1, h 1 = 0.1, h 2 = 0.2, h 3 = 0.1, E ( 2 ) = 20 GPa , ρ ( 2 ) = 950 / g N/m 3 , v ( 2 ) = 0.25, g = 9.81 m / s 2

5.1 Transient Response of Heterogeneous Sandwich Plate

In this study, dynamic response of sandwich circular plate with functionally graded orthotropic face sheets are presented as first time and no available results may be found in literature to demonstrate the efficiency and accuracy of the obtained results. For this reason, as a verification example, the obtained results are verified by comparison with finite element results for some special cases. The finite element results are extracted from the ABAQUS software based on the three dimensional theory of elasticity. In this example, sandwich plate with heterogeneous face sheets and asymmetric layups is considered and effects of the damping coefficient on the transient response of heterogeneous sandwich plate are examined. The face sheets have the same stiffness in the radial and circumferential directions as:

{ E r ( 1 ) ( r ) = E θ ( 1 ) ( r ) = 310 [ 1 ( r r o ) 2 ] GPa ρ ( 1 ) ( r ) = 1613 [ 1 ( r r o ) 2 ] / g N/m 3 v r θ ( 1 ) = v θ r ( 1 ) = 0.26 , { E r ( 3 ) ( r ) = E θ ( 3 ) ( r ) = 310 [ 1 + 0.5 ( r r o ) ] GPa ρ ( 3 ) ( r ) = 1613 [ 1 + 0.5 ( r r o ) ] / g N/m 3 v r θ ( 3 ) = v θ r ( 3 ) = 0.23 , G ( k ) ( r ) = E r ( k ) ( r ) 2 ( 1 + v r θ ( k ) )

Dynamic responses of the sandwich annular plate with clamped-free boundary condition (clamped inner edge and free outer edge) are shown in Figures 3-5 for the load cases (I), (II) and (III), respectively. In these figures, time variations of the lateral deflection at the outer edge of sandwich plate are plotted for k w =0 and various values of the damping coefficients of viscoelastic foundation C t . Also, the obtained results are compared with results of ABAQUS finite element software when there is no foundation (k w = C t =0). It can be seen that there is an excellent agreement between present results based on the analytical solution (AS) and results of ABAQUS software which are extracted based on the finite element method (FEM). It is shown that when the damping coefficient increases, vibration amplitude becomes smaller, as expected due to the dissipation of system energy. Figure 3 shows that the plate oscillates about the static deflection (W Static ) and for C t =0 (when there is no damping in the foundation) the maximum deflection will be as:W Dynamic-Maximum =2W Static . Also for plate subjected to step loading rested on viscous foundation, as time goes on, the dynamic lateral deflection tends to static response. Time variations of the lateral deflection of the plate under load case (II) are plotted in Figure 4 for time duration τ = 0.5 and 1.5 ms. In this type of loading, the dynamic responses of the plate are separated into two stages as follow:

  • a) Forced vibration (plate under the step loading) t < τ

The plate oscillates about the static deflection of the plate under the step loading.

  • b) Free vibration (the load is removed) t > τ

The plate oscillates about the zero.

Figure 3:
Dynamic response of clamped-free sandwich annular plate subjected to an abrupt uniformly distributed transverse load (case (I)).

Figure 4:
Dynamic response of clamped-free sandwich annular plate subjected to an abruptly-imposed uniform transverse load (Case (II)) with time duration: (a) t 0 = 0.5 ms and (b) t 0 = 2 ms.

Figure 5:
Dynamic response of clamped-free sandwich annular plate subjected to a uniform transverse load whose intensity varies according to a sinusoidal function (Case (III)): (a) ω = Ω, (b) ω = Ω/4 (c) ω = Ω/6 (d) ω = Ω/8.

As expected for plate subjected to a pulse loading rested on viscous foundation, the amplitude decreases and approaches zero as time increases.

Dynamic responses of sandwich plate under sinusoidal loading (case (III)) are illustrated in Figure 5 for various frequency of the harmonic excitation (ω). Value of the parameter ω is chosen as a ratio of the fundamental natural frequency (Ω = 3892.5) of the plate. As it is expected for plate under sinusoidal load with harmonic excitation ω = Ω, the oscillation amplitude becomes unbounded in undamped system (Figure 5(a)). However, there is not an obvious increase in the amplitude of oscillation when the damping coefficient increases up to 3E5 or 8E5.

It can be seen from Figures 5(b), (c) and (d) that for excitation frequencies ω = Ω/ 4, Ω/ 6 and Ω/ 8, influence of other vibration modes are superimposed on the main response in undamped system. On the other hands the dynamic response of plates resting on the viscous foundation (damped system) is smoother and the local oscillations are eliminated. It is observed that the peak deflection slightly decreases as excitation frequency decreases. These figures also show that increasing the damping coefficient of viscous foundation reduces the amplitude and frequency of vibration.

5.2 Transient Response of Functionally Graded Polar Orthotropic Sandwich Plate

In the present section, dynamic responses of circular sandwich plate fabricated from functionally graded polar orthotropic face sheets are investigated for various edge conditions. Also, effects of stiffness and damping coefficients of viscoelastic foundation on the transient response of sandwich plate are examined for various dynamic loads.

Variation of the Young’s modulus (in the radial and circumferential direction), shear modulus and density of each face sheet in the radial direction are chosen as:

  • a) Top face sheet:

{ E r ( 1 ) ( r ) E θ ( 1 ) ( r ) G r z ( 1 ) ( r ) ρ ( 1 ) ( r ) } = { 170 GPa 68 GPa 34 GPa 1413 / g N/m 3 } [ 1 + 0.4 ( r r o ) ]

  • b) Bottom face sheet:

{ E r ( 3 ) ( r ) E θ ( 3 ) ( r ) G r z ( 3 ) ( r ) ρ ( 3 ) ( r ) } = { 68 GPa 170 GPa 18 GPa 1413 / g N/m 3 } [ 1 + 0.7 ( r r o ) ]

Also, the Poisson’s ratios are assumed as: vrθ(1)=0.35, vrθ(3)=0.14 where vrθ Erθ=vθrEθr

Dynamic responses of clamped-simply supported sandwich plate subjected to load cases (I), (IV) and (V) are plotted in Figures 6-8, respectively. In these figures, time variations of the lateral deflection are presented for r=0.5. It is evident that the effects of elastic foundation (Winkler-type foundation) can be examined when damping coefficient C t is set to zero. Figure 6 shows that vibration amplitude and frequency of plate under elastic foundation (C t =0 and k w =2E10) decreases and increases, respectively, comparing with those of the plate without foundation. The plate oscillates about the static value of deflection when there is no damping in the foundation (C t =0), in which the vibration amplitude is constant. Also, the vibration amplitude of plate under viscous foundation decreases with time and tends to the static response of plate for C t =0.

Figure 6:
Dynamic response of clamped-simply supported sandwich annular plate subjected to an abrupt uniformly distributed transverse load (case (I)).

Figure 7:
Dynamic response of clamped- simply supported sandwich annular plate subjected to an N-shaped pulse loading (case (IV)).

Figure 8:
Dynamic response of clamped-simply supported sandwich annular plate subjected to a sinusoidal loading acting in the same direction (case (V)) with time duration: (a) τ =2 ms and (b) τ = 4 ms.

It can be seen from Figure 7 that for sandwich plate subjected to an N-shaped pulse loading, the center of oscillation decreases with time and becomes negative when τ < t (forced vibration), when the load is removed (τ > t free vibration) the plate oscillates about the zero. As expected, vibration amplitude of plate without damping (C t =0) is constant in each stage, but for plate subjected to a pulse loading rested on viscous foundation, the amplitude decreases with time and the lateral deflection tends to zero. Also for plate without foundation (C t =k w =0), vibration amplitude in the first stage (τ < t) is higher than that of the second stage (τ > t), but for plate under elastic foundation with C t =0 and k w =2E10, vibration amplitude in the second stage is higher. Indeed, in the first stage, vibration amplitude of plate without foundation is higher than that in the other cases and in the second stage, vibration amplitude of plate with C t =0 and k w =2E10 is higher.

Figure 8 shows that for the plate subjected to load case (V), the response to the loading with τ = 2ms is higher than the response to the loading with τ = 4ms. In the other words, the loading with the higher frequency of excitation leads to higher response.

This figure also shows that influence of other vibration modes which are superimposed on the main response of plate in undamped system, are eliminated for plates resting on the viscous foundation.

Effects of viscoelastic foundation on the dynamic response of clamped-clamped sandwich plate subjected to load cases (I), (VI) and (VII) are presented in Figures 9-11, respectively. In these figures, the results are reported for r=0.5. It can be seen that in comparison with the forgoing example, there is a same trend for the dynamic response of sandwich plate under step loading with clamped-clamped and clamped-simply supported boundary conditions (Figures 6 and 9). However, the vibration amplitude and frequency of clamped-clamped plate decreases and increases, respectively, comparing with those of the clamped-simply supported plate due to increasing the stiffness of plate with clamped-clamped edges.

Figure 9:
Dynamic response of clamped-clamped sandwich annular plate subjected to an abrupt uniformly distributed transverse load (case (I)).

Figure 10:
Dynamic response of clamped-clamped sandwich annular plate subjected to a triangular-shaped pulse loading (case (VI)).

Figure 11:
Dynamic response of clamped-clamped sandwich annular plate subjected to an exponential loading (case (VII)): (a) ζ = 100 (b) ζ = 500.

It is observed that for dynamic response of plate under the triangular-shaped pulse and exponential loadings (Figures 10 and 11), the center of oscillation decreases and approaches zero as time increases. These figures show that the vibration amplitude of plate without foundation (C t =k w =0) is higher than that in the other cases.

Time variations of the lateral deflection at the outer edge of the sandwich plate with clamped-free boundary condition subjected to load cases (VIII) and (IX) are plotted in Figure 12 and 13, respectively. In Figure 12, results of sandwich plate under repeated rectangular-shaped loading (load case (VII)) with time duration of 2ms and 4ms are illustrated.

Figure 12:
Dynamic response of clamped-free sandwich annular plate subjected to a repeated rectangular-shaped loading (Case (VIII)) with time duration: (a) τ = 2 ms and (b) τ = 4 ms.

Figure 13:
Dynamic response of clamped-free sandwich annular plate subjected to a repeated triangular-shaped loading (Case (IX)) with time duration: (a) τ = 2 ms and (b) τ = 5 ms.

It can be seen that in the first duration of loading (t < τ), the peak deflection of plate without foundation is higher than the other cases and then (t > τ) the peak deflection of plate under elastic foundation with K w =5E9 in Figure 12(a) and K w =2E9 in Figure 12(b) is higher. Increasing elastic foundation stiffness value up from 5E9 to 9E9 in Figure 12(a) and 2E9 to 2E10 in Figure 12(b), causes an obvious decrease of vibration amplitude.

It is interesting to note that the viscoelastic foundation causes a reduction in the vibration amplitude of plate in the loading time durations (2 iτ< t < ( 2i + 1)τ i = 0,1,2,...), but in the time durations when the load is removed ((2i+1)τ < t < 2(i + 1)τ i = 0,1,2,...), the vibration amplitude of plate without foundation is smaller and its response is smoother.

Figure 13 shows that the oscillation amplitude of plate subjected to the repeated triangular-shaped loading is unbounded in undamped system. It is seen that the repeated triangular-shaped loading with lower time durations (τ ) increases the amplitude rapidly. Also, vibration amplitude of plate under viscoelastic foundation is bounded in some cases of viscoelastic foundation, where the stiffness and viscosity coefficients of foundation become higher than specific values.

6 CONCLUSIONS

By using the layerwise theory, effects of the viscoelastic foundation on the dynamic response of sandwich annular plates are studied. The viscoelastic substrate is modeled as Kelvin-Voigt foundation. Sandwich plates may be fabricated from functionally graded polar orthotropic face sheets, where Young’s modulus in the radial and circumferential direction, shear modulus and density of each face sheet may be varied continuously in the radial direction. The relatively complicated second order coupled partial differential equations are solved based on a developed semi-analytical solution procedure by using a combination of the finite Taylor’s transform procedure in the space domain and the fourth-order Runge-Kutta method in the time domain. The results are obtained for various dynamic loads and edge conditions. For some special cases, the obtained results based on the proposed solution procedure are compared with finite element results based on the three dimensional theory of elasticity. The comparisons show that there is a very good agreement between present results and results of the three-dimensional theory of elasticity. Results reveal that FG polar orthotropic sandwich plates under various transient loads, edge conditions and foundation parameters may be analyzed by using the proposed solution procedure. Results also show that influence of other vibration modes which are superimposed on the main response of plate in undamped system, are eliminated for plates resting on the viscous foundation.

References

  • Alipour, M.M., (2016a). A novel economical analytical method for bending and stress analysis of functionally gradedsandwich circular plates with general elastic edge conditions, subjected to various loads, Composite Part B: Engineering95: 48-63.
  • Alipour, M.M., (2016b). Effects of elastically restrained edges on FG sandwich annular plates by using a novel solutionprocedure based on layerwise formulation, Archives of Civil and Mechanical Engineering 16: 678-694.
  • Alipour, M.M. (2016c). An analytical approach for bending and stress analysis of cross/angle-ply laminated composite plates under arbitrary non-uniform loads and elastic foundations, Archives of Civil and Mechanical Engineering 16(2): 193-210.
  • Alipour, M.M. (2018). Transient forced vibration response analysis of heterogeneous sandwich circular plates under viscoelastic boundary support, Archives of Civil and Mechanical Engineering 18: 12-31.
  • Alipour, M.M., Shariyat, M. (2015). Analytical zigzag formulation with 3D elasticity corrections for bending and stress analysis of circular/annular composite sandwich plates with auxetic cores, Composite Structures 132: 175-197.
  • Carrera, E., Brischetto, S. (2008). A survey with numerical assessment of classical and refined theories for the analysis of sandwich plates. Applied Mechanics Reviews 62(1): 010803-010819.
  • Ding, H., Shi, K.L., Chen, L.Q., Yang, S.P. (2013). Dynamic response of an infinite Timoshenko beam on a nonlinear viscoelastic foundation to a moving load, Nonlinear Dynamics 73: 285-298.
  • Ghashochi-Bargh, H., Sadr, M.H., (2014). Vibration reduction of composite plates by piezoelectric patches using a modifiedartificial bee colony algorithm, Latin American Journal of Solids and Structures 11: 1846-1863.
  • Ghorbanpour Arani, A., Shiravand, A., Rahi, M., Kolahchi, R., (2012). Nonlocal vibration of coupled DLGS systems embedded on Visco-Pasternak foundation, Physica B. 407: 4123-4131.
  • Ghorbanpour Arani, A., Maboudi, M.J., Kolahchi, R. (2014), Nonlinear vibration analysis of viscoelastically coupled DLAGS-systems. European Journal of Mechanics - A/Solids 45: 185-197.
  • Hoang, T., Duhamel, D., Foret, G., Yin, H.P., Cumunel, G., (2016). Response of a periodically supported beam on a nonlinearfoundation subjected to moving loads, Nonlinear Dynamics 86(2): 953-961.
  • Hosseini Hashemi, S.h., Mehrabani, H., Ahmadi-Savadkoohi, A., (2015). Forced vibration of nanoplate on viscoelastic substrate with consideration of structural damping: An analytical solution, Composite Structures 133: 8-15.
  • Houari, M. S. A., Tounsi, A., Bég O. A., (2013). Thermoelastic bending analysis of functionally graded sandwich plates using a new higher order shear and normal deformation theory, International Journal of Mechanical Sciences 76: 102-111
  • Karim, K.R., Chen, G.D., (2012). Surface damping effect of anchored constrained viscoelastic layers on the flexural response of simply supported structures, Mechanical Systems and Signal Processing 27: 419-432.
  • Kiasat, M.S., Zaman, H.A., Aghdam, M.M., (2014). On the transient response of viscoelastic beams and plates on viscoelasticmedium, International Journal of Mechanical Sciences 83: 133-145.
  • Kim, S.M., McCullough, B.F., (2003). Dynamic response of plate on viscous Winkler foundation to moving loads of varyingamplitude, Engineering Structures 25: 1179-1188.
  • Lee, K.H., Cao, L., (1996). A predictor-corrector zig-zag model for the bending of laminated composite plates, InternationalJournal of Solids Structures 33(6): 879-897.
  • Lepoittevin, G., Kress, G., (2011). Finite element model updating of vibrating structures under free-free boundary conditionsfor modal damping prediction, International Journal of Mechanical Sciences 25: 2203-2218.
  • Liang, X.u., Wang, Z., Wang, L., Liu, G., (2014). Semi-analytical solution for three-dimensional transient response of functionallygraded annular plate on a two parameter viscoelastic foundation, Journal of Sound and Vibration 333: 2649-2663.
  • Luo, W.L., Xia, Y., Zhou, X.Q., (2016). A closed-form solution to a viscoelastically supported Timoshenko beam under harmonicline load, Journal of Sound and Vibration 369: 109-118.
  • Luong-Van, H., Nguyen-Thoi, T., Liu, G.R., Phung-Van, P., (2014). A cell-based smoothed finite element method usingthree-node shear-locking free Mindlin plate element (CS-FEM-MIN3) for dynamic response of laminated compositeplates on viscoelastic foundation, Engineering Analysis with Boundary Elements 42: 8-19.
  • Mazilu, T., Dumitriu, M., Tudorache, C., (2012). Instability of an oscillator moving along a Timoshenko beam on viscoelasticfoundation, Nonlinear Dynamics 67: 1273-1293.
  • Metrikine, A.V., Verichev, S.N., (2001). Instability of vibrations of a moving two-mass oscillator on a flexibly supported Timoshenko beam, Archives of Civil and Mechanical Engineering 71: 613-624.
  • Mohammadimehr, M., Navi, B.R., Arani, A.G., (2015). Free vibration of viscoelastic double-bonded polymeric nanocompositeplates reinforced by FG-SWCNTs using MSGT, sinusoidal shear deformation theory and meshless method, CompositeStructures 13: 654-671.
  • Nie, G.J., Batra, R.C., (2010). Static deformations of functionally graded polar-orthotropic cylinders with elliptical innerand circular outer surfaces, Composites Science and Technology 70: 450-457.
  • Pan, Q., Tian He, T., Xiao, D., Liu, X., (2016). Design and Damping Analysis of a New Eddy Current Damper for AerospaceApplications, Latin American Journal of Solids and Structures 13: 1997-2011
  • Peng, X.L., Li, X.F., (2012). Elastic analysis of rotating functionally graded polar orthotropic disks, International Journal ofMechanical Sciences 60: 84-91.
  • Plattenburg, J., Dreyer, J. Singh, T. R., (2016). A new analytical model for vibration of a cylindrical shell and cardboardliner with focus on interfacial distributed damping. Mechanical Systems and Signal Processing 75: 176-195.
  • Pouresmaeeli, S., Ghavanloo, E., Fazelzadeh, S.A., (2013). Vibration analysis of viscoelastic orthotropic nanoplates restingon viscoelastic medium. Composite structures 96: 405-410.
  • Reddy, J. (2006). An introduction to the Hinite element method. 2ndEdition, New York: Wiley.
  • Shariyat, M., Alipour, M.M., (2015). Novel Layerwise Shear Correction Factors for Zigzag Theories of Circular SandwichPlates with Functionally Graded Layers, Latin American Journal of Solids and Structures 12: 1362-1396.
  • Shariyat, M., Alipour, M.M., (2017a). Analytical layerwise free vibration analysis of circular/annular composite sandwichplates with auxetic cores, International Journal of Mechanics and Materials in Design 13, 125-157.
  • Shariyat, M., Alipour, M.M., (2017b). Analytical bending and stress analysis of variable thickness FGM auxetic conical/cylindrical shells with general tractions, Latin American Journal of Solids and Structures 14: 805-843.
  • Tounsi, A., Houari, M. S. A.,Benyoucef, S.,Bedia, E. A. A., (2013). A refined trigonometric shear deformation theory for thermoelastic bending of functionally graded sandwich plates, Aerospace Science and Technology24 (1): 209-220.
  • Yang, Y., Ding, H., Chen, L.Q., (2013). Dynamic response to a moving load of a Timoshenko beam resting on a nonlinearviscoelastic foundation, Acta Mechanica Sinica 29(5): 718-727.
  • Zamani H.A., Aghdam, M.M., Salehi, M., (2017). Free damped vibration analysis of Mindlin plates with hybrid material foundation viscoelasticity. International Journal of Mechanical Sciences, 121: 33-43.
  • Zenkour, A.M., (2016). Nonlocal transient thermal analysis of a single-layered graphene sheet embedded in viscoelastic medium, Physica E 79: 87-97.
  • Zhang, D.P., Lei, Y., Shen, Z.B., (2016). Free transverse vibration of double-walled carbon nanotubes embedded in viscoelastic médium, Acta Mechanica 227(12): 3657-3670.

Publication Dates

  • Publication in this collection
    Dec 2017

History

  • Received
    07 June 2017
  • Reviewed
    29 Aug 2017
  • Accepted
    04 Sept 2017
Individual owner www.lajss.org - São Paulo - SP - Brazil
E-mail: lajsssecretary@gmsie.usp.br