This study follows an analytical approach to investigate the nonlinear dynamic response and vibration of eccentrically stiffened sandwich functionally graded material (FGM) cylindrical panels with metal–ceramic layers on elastic foundations in thermal environments. It is assumed that the FGM cylindrical panel is reinforced by the eccentrically longitudinal and transversal stiffeners and subjected to mechanical and thermal loads. The material properties are assumed to be temperature dependent and graded in the thickness direction according to a simple power law distribution. Based on the Reddy’s third-order shear deformation shell theory, the motion and compatibility equations are derived taking into account geometrical nonlinearity and Pasternak-type elastic foundations. The outstanding feature of this study is that both FGM cylindrical panel and stiffeners are assumed to be deformed in the presence of temperature. Explicit relation of deflection–time curves and frequencies of FGM cylindrical panel are determined by applying stress function, Galerkin method and fourth-order Runge-Kutta method. The influences of material and geometrical parameters, elastic foundations and stiffeners on the nonlinear dynamic and vibration of the sandwich FGM panels are discussed in detail. The obtained results are validated by comparing with other results in the literature.

Functionally graded materials (FGMs) are microscopically inhomogeneous made from a mixture of metal and ceramic, and their mechanical properties vary smoothly and continuously from one surface to the other. Functionally graded structures, such as cylindrical panels and cylindrical shells in recent years, play an important part in the modern industries. They are used as primary load-carrying parts in spacecraft and missile structures which are usually exposed to very severe loading conditions. Therefore, globally, researches on static and dynamic stability of these structures have received special attention. Song et al. [1] studied the active vibration control of carbon nanotube-reinforced functionally graded composite cylindrical shell using piezoelectric materials. Mehralian et al. [2] investigated buckling of anisotropic piezoelectric cylindrical shells subjected to axial compression and lateral pressure based on the new modified couple stress theory and using the shear deformation theory with the von Kármán geometrical nonlinearity. Mirzaei and Kiani [3] researched free vibration characteristics of composite plates reinforced with single-walled carbon nanotubes. Sofiyev et al. [4] dealt with the nonlinear vibration of orthotropic cylindrical shells on the nonlinear elastic foundations. Hosseini-Hashemi et al. [5] concentrated on presenting a reliable and accurate method for free vibration of functionally graded viscoelastic cylindrical panel under various boundary conditions. Based on 3D elasticity theory, Norouzi and Alibeigloo [6] carried out thermoviscoelastic analysis of FGM cylindrical panel under thermal and/or mechanical load, while Ahmed [7] investigated how the corrugation parameters and the Winkler foundation affect the buckling behavior of isotropic and orthotropic thin-elliptic cylindrical shells with cosine-shaped meridian subjected to radial loads. Quan et al. [8] focused on the nonlinear dynamic analysis and vibration of shear deformable imperfect eccentrically stiffened functionally graded thick cylindrical panels, taking into account the damping when subjected to mechanical loads. Further, Chen et al. [9] studied the structural–acoustic radiation problem of cylindrical shell structures with complex acoustic boundary conditions. Shen [10] presented modeling and analysis for the postbuckling of carbon nanotube-reinforced composite cylindrical panels resting on elastic foundations subjected to lateral pressure in thermal environments. Recently, Duc et al. [11] investigated the nonlinear buckling and postbuckling of eccentrically stiffened functionally graded thin-elliptical cylindrical shells surrounded on elastic foundations in thermal environment.

A thick structure is defined as a structure with a thickness which is large compared to its other dimensions, but in which deformations are still not large compared to thickness. When studying nonlinear analysis of thick structures, first-order shear deformation theory and higher order shear deformation theory are usually used, in which the effects of transverse shear and normal stresses in structures are taken into account. Despite their complexity, the mechanical behaviors of thick structures, such as bending, vibration, stability, buckling, etc., have attracted the attention of many researchers. Zhou and Zhu [12] utilized the third-order shear deformation plate theory to analyze the vibration and bending of the simply supported magneto-electro-elastic rectangular plates. Sayyaadi et al. [13] presented an analytical solution for power output from a piezoelectric shallow shell energy harvester using higher order shear deformation theory, and Selim et al. [14] studied the free vibration behavior of carbon nanotube-reinforced functionally graded composite plates in a thermal environment based on Reddy’s higher order shear deformation theory. Wattanasakulpong et al. [15] employed an improved third-order shear deformation theory to investigate free and forced vibration responses of FGM plates. Yang and Shen [16] and Huang and Shen [17] dealt with the nonlinear vibration and dynamic response of functionally graded material plates in thermal environments. Amirpour et al. [18] introduced the deformation solution of FGM plates with variation of material stiffness through their length using higher order shear deformation theory including stretching effects. Groh and Weaver [19] discussed the static inconsistencies that arise when modelling the flexural behavior of beams, plates, and shells with clamped boundary conditions using a certain class of axiomatic, higher order shear deformation theory. He et al. [20] presented a finite element formulation based on the classical laminated plate theory for the shape and vibration control of the FGM plates with integrated piezoelectric sensors and actuators. Quan and Duc [21] investigated the nonlinear static and dynamic stability of imperfect eccentrically stiffened FGM higher order shear deformable double-curved shallow shell on elastic foundations in thermal environments. Shao et al. [22] presented a unified formulation which is based on a general refined shear deformation beam theory to conduct free vibration analysis of composite laminated beams subjected to general boundary conditions.

Stiffeners are secondary plates or sections which are attached to structures to stiffen them against out-of-plane deformations. Stiffeners can be critical to the performance of the structures since a properly placed stiffener can increase the capacity of the member it supports. Therefore, in recent years, many investigations have been carried out on the mechanical behaviors of structures which are reinforced by stiffeners. Yang et al. [23] tested eight stiffened square concrete-filled steel tubular stub columns with slender sections of encasing steel and two non-stiffened counterparts subjected to axial compressive load. Xu et al. [24] developed a new effective smeared stiffener method to compute the global buckling load of grid-stiffened composite panels. Turner and Vizzini [25] conducted a study to determine the effect of integral stiffeners on the damage growth and ultimate strength of sandwich panels with impact damage. Tong and Guo [26] investigated the elastic buckling behavior of steel trapezoidal corrugated shear walls with vertical stiffeners. Feng et al. [27] studied the effect of impact damage positions on the buckling and postbuckling behaviors of stiffened composite panels under axial compression. Duc et al. [28] presented an analytical approach to investigate the nonlinear stability analysis of eccentrically stiffened thin FGM cylindrical panels on elastic foundations subjected to mechanical loads, thermal loads and the combination of these loads. Farahani et al. [29] dealt with an analytical approach of the buckling behavior of an FGM circular cylindrical shell under axial pressure with external axial and circumferential stiffeners. Shen et al. [30] developed a theoretical model to predict sound transmission loss across periodically and orthogonally stiffened composite laminate sandwich structures using the first-order shear deformation theory. Recently, Duc [31] studied the nonlinear dynamic response of higher order shear deformable sandwich functionally graded circular cylindrical shells with outer surface-bonded piezoelectric actuator on elastic foundations subjected to thermo-electro-mechanical and damping loads.

New contribution of the paper is that this is the first investigation that successfully established the modeling and analytical formulations for the nonlinear vibration and dynamic response of eccentrically stiffened sandwich FGM thick cylindrical panels with metal–ceramic layers subjected to mechanical loads in thermal environments using the Reddy’s third-order shear deformation shell theory. The outstanding feature of this study is that both FGM cylindrical panel and stiffeners are assumed to be deformed in the presence of temperature. The nonlinear equations are solved by the Galerkin method and fourth-order Runge-Kutta method.

Consider an eccentrically stiffened FGM cylindrical panel with the radii of curvature, thickness, axial length and arc length of the panel are R, h, a and b, respectively. The FGM cylindrical panel is assumed to be rested on elastic foundations (Figure 1). A coordinate system (x,y,z) is established in which (x,y) plane on the middle surface of the panel and z on thickness direction (-h/2zh/2) as shown in Figure 1. The panel is reinforced by eccentrically longitudinal and transversal stiffeners (Figure 2). The width and thickness of longitudinal and transversal stiffeners are denoted by dx,hx and dy,hy, respectively; sx,sy are the spacing of the longitudinal and transversal stiffeners. The quantities Ax,Ay are the cross-sectional areas of stiffeners, and Ix,Iy,zx,zy are the second moments of cross-sectional areas and the eccentricities of stiffeners with respect to the middle surface of panel, respectively. E0 is Young’s modulus of stiffeners. In order to provide continuity between the panel and stiffeners, the stiffeners are made of full ceramic (E0=Ec).


                        figure

Figure 1. Geometry and coordinate system of eccentrically stiffened sandwich FGM cylindrical panels on elastic foundations.


                        figure

Figure 2. Configuration of longitudinal and transversal stiffeners of eccentrically stiffened sandwich FGM cylindrical panels.

By applying the power law distribution, the volume fractions of ceramic and metal of the FGM cylindrical panels are assumed as [21,28]

Vm(z)=(2z+h2h)N,Vc(z)=1-Vm(z)(1)
where N is volume fraction index (0N<), subscripts m and c stand for the metal and ceramic constituents, respectively.

It is assumed that the effective properties Preff(z,T) of FGM cylindrical panel, such as the elastic modulus E(z,T), the mass density ρ(z,T) and the thermal expansion coefficient α(z,T), vary in the thickness direction z and are determined by linear rule of mixture as

Preff(z,T)=Prc(T)Vc(z)+Prm(T)Vm(z)(2)

A material property Pr in equation (2) can be expressed as a nonlinear function of temperature [10,28,32]

Pr(T)=P0(P-1T-1+1+P1T+P2T2+P3T3)(3)
where T=T0+ΔT, ΔT is the temperature increment in the environment containing the cylindrical panel and T0=300K (room temperature), P0,P-1,P1,P2 and P3 are coefficients characterizing the constituent materials and determined by experiments.

The effective properties in equation (2) of the FGM cylindrical panel are obtained by substituting equation (1) into equation (2) as

[E(z,T),ρ(z,T),α(z,T)]=[Em(T),ρm(T),αm(T)]+[Ecm(T),ρcm(T),αcm(T)](2z+h2h)N(4)
where
Ecm(z,T)=Ec(z,T)-Em(z,T),ρcm(z,T)=ρc(z,T)-ρm(z,T),αcm(z,T)=αc(z,T)-αm(z,T)(5)
and the Poisson’s ratio is assumed to be constant ν(z)=v=const [33,34].

The cylindrical panel–foundation interaction of Pasternak model is given by

qe=k1w-k22w(6)
where 2=2/x2+2/y2, w is the deflection of the cylindrical panel, k1 is Winkler foundation modulus and k2 is the shear layer foundation stiffness of Pasternak model.

In this study, the Reddy’s third order-shear deformation shell theory is used to derive basic equations to investigate the nonlinear dynamic response and vibration of eccentrically stiffened FGM cylindrical panels on elastic foundations in thermal environments.

The normal strains ɛx,ɛy, in-plane shear strain γxy and transverse shear deformations γxz,γyz at the distance z from the mid-plane are given as

[ɛxɛyγxy]=[ɛx0ɛy0γxy0]+z[kx1ky1kxy1]+z3[kx3ky3kxy3],[γxzγyz]=[γxz0γyz0]+z2[kxz2kyz2](7)
where
[ɛx0ɛy0γxy0]=[ux+12(wx)2vy-wR+12(wy)2uy+vx+wxwy],[γxz0γyz0]=-3c1[wx+φxwy+φy],[kx1ky1kxy1]=[φxxφyyφxy+φyy],[kx3ky3kxy3]=-c1[φxx+2wx2φyy+2wy2φxy+φyy+22wxy],[kxz2kyz2]=-3c1[wx+φxwy+φy](8)
in which c1=4/(3h2) and geometrical nonlinearity is incorporated. Also, u,v are the displacement components along the x,y directions, respectively, and φx,φy are the rotations of normals to the mid-surface with respect to y and x axes, respectively.

Hooke’s law for an FGM cylindrical panel with temperature-dependent properties is defined as [8,11]

(σx,σy)=E(z,T)1-ν2[(ɛx,ɛy)+ν(ɛy,ɛx)-(1+ν)α(z,T)ΔT(1,1)](σxy,σxz,σyz)=E(z,T)2(1+ν)(γxy,γxz,γyz)(9)

The stress–strain relations of the stiffeners can be given as follows

(σxs,σys)=E0(z,T)(ɛx,ɛy)-E0(z,T)(1-2ν0)α0(z,T)ΔT(1,1)(10)

In this paper, we assumed that the force and moment resultants of an eccentrically stiffened FGM cylindrical panel under temperature can be constructed by the combination of the integrations of the stresses over the panel thickness and those of corresponding stiffener contributions (including the thermal stress in the stiffeners) as follows

(Ni,Mi,Pi)=-h/2h/2σi(1,z,z3)dz+h/2h/2+hiσis(1,z,z3)diTsiTdz,i=x,y(Nxy,Mxy,Pxy)=-h/2h/2σxy(1,z,z3)dz,(Qi,Ri)=-h/2h/2σiz(1,z2)dz,i=x,y(11)

It is assumed that both material properties and geometrical shapes of stiffeners depend on temperature, and they are deformed in the presence of temperature. After the thermal deformation process, the geometrical parameters of stiffeners can be determined as [8,11,32]

dxT=dx(1+α0ΔT),dyT=dy(1+α0ΔT),hxT=hx(1+α0ΔT),hyT=hy(1+α0ΔT)sxT=sx(1+α0ΔT),syT=sy(1+α0ΔT),zxT=zx(1+α0ΔT),zyT=zy(1+α0ΔT)(12)

Substitution of equation (7) into equations (9) and (10) and the result into equation (11) yields the constitutive relations as

Nx=B11ɛx0+B12ɛy0+B13kx1+B14ky1+B15kx3+B16ky3-B17Φ1-B18Φ1xs,Ny=B12ɛx0+B22ɛy0+B24ky1+B14kx1+B16kx3+B26ky3-B17Φ1-B28Φ1ys,Nxy=B31γxy0+B32kxy1+B33kxy3,Mx=B13ɛx0+B14ɛy0+B43kx1+B44ky1+B45kx3+B46ky3-B17Φ2-B18Φ2xs,My=B14ɛx0+B24ɛy0+B44kx1+B54ky1+B46kx3+B56ky3-B17Φ2-B28Φ2ys,Mxy=B32γxy0+B62kxy1+B63kxy3,Px=B71ɛx0+B16ɛy0+B73kx1+B46ky1+B75kx3+B76ky3-B17Φ4-B18Φ4xs,Py=B16ɛx0+B82ɛy0+B46kx1+B84ky1+B76kx3+B86ky3-B17Φ4-B28Φ4ys,Pxy=B33γxy0+B63kxy1+B93kxy3,Qx=B31γxz0+B62kxz2,Qy=B31γyz0+B62kyz2,Rx=B62γxz0+B63kxz2,Ry=B62γyz0+B63kyz2(13)
with the detail of coefficients B1i(i=1,8¯),B2j(j=2,4,6,8),B3k(k=1,3¯),B4l(l=3,6¯), B54,B56,B62,B63,B7m(m=1,3,5,6),B8n(n=2,4,6),B93 are given in Appendix 1.

The nonlinear motion equations of FGM cylindrical panel are [32]

Nxx+Nxyy=I1¯2ut2+I2¯2φxt2-I3¯3wt2x(14a)
Nxyx+Nyy=I1*¯2vt2+I2*¯2φyt2-I3*¯3wt2y(14b)
Qxx+Qyy-3c1(Rxx+Ryy)+c1(2Pxx2+22Pyxy+2Pyy2)+NyR+q+Nx2wx2+2Nxy2wxy+Ny2wy2-k1w+k22w=I12wt2+2ɛI1wt+I3¯3ut2x+I5¯3φxt2x+I3*¯3vt2y+I5*¯3φyt2y-c12I7(4wt2x2+4wt2y2)(14c)
Mxx+Mxyy-Qx+3c1Rx-c1(Pxx+Pxyy)=I2¯2ut2+I4¯2φxt2-I5¯3wt2x(14d)
Mxyx+Myy-Qy+3c1Ry-c1(Pxyx+Pyy)=I2*¯2vt2+I4*¯2φyt2-I5*¯3wt2y(14e)
in which q is an external pressure uniformly distributed on the surface of the FGM cylindrical panel, ɛ is the viscous damping coefficient and
I1¯=I1,I1*¯=I1+2I2R,I2¯=I2-c1I4,I2*¯=I2+I3R-c1I4-c1I5R,I3¯=c1I4,I3*¯=c1I4+c1I5R,I4¯=I4*¯=I3-2c1I5+c12I7,I5¯=I5*¯=c1I5-c12I7,(I1,I2,I3,I4,I5,I7)=-h/2h/2ρ(z)(1,z,z2,z3,z4,z6)dz(15)

The geometrical compatibility equation for an FGM cylindrical panel is written as [32]

2ɛx0y2+2ɛy0x2-2γxy0xy=2wxy2-2wx22wy2-1R2wx2(16)

From equation (13), one can write

ɛx0=I11*2fy2-I12*2fx2+I13*φxx+I14*φyy-c1I15*(2wx2+φxx)-c1I16*(2wy2+φyy)+I17*Φ1+I18*Φ1xs+I19*Φ1ys,ɛy0=I21*2fx2-I12*2fy2+I23*φxx+I24*φyy-c1I25*(2wx2+φxx)-c1I26*(2wy2+φyy)+I27*Φ1+I29*Φ1xs+I29*Φ1ys,γxy0=-I31*2fxy+I32*(φxy+φyx)-c1I33*(22wxy+φxy+φyx)(17)
in which
Δ=I11I22-I122,I11*=I22Δ,I12*=I12Δ,I13*=I14I12-I13I22Δ,I14*=I24I12-I14I22Δ,I15*=I16I12-I15I22Δ,I16*=I26I12-I16I22Δ,I17*=I22I17-I12I17Δ,I18*=I18I22Δ,I19*=-I28I12Δ,I21*=I11Δ,I23*=I13I12-I11I14Δ,I24*=I14I12-I11I24Δ,I25*=I15I12-I11I16Δ,I26*=I16I12-I11I26Δ,I27*=I11I17-I12I17Δ,I28*=-I18I12Δ,I29*=I11I28Δ,I31*=1I31,I32*=-I32I31,I33*=-I33I31(18)
and the stress function f(x,y,t) is defined as
Nx=2fy2,Ny=2fx2,Nxy=-2fxy(19)

Imposing equation (19) into equations (14a) and (14b) yields

2ut2=-I2¯I1¯2φxt2+I3¯I1¯3wt2x(20a)
2vt2=-I2*¯I1*¯2φyt2+I3*¯I1*¯3wt2y(20b)

Substituting equations (8), (13), (20a) and (20b) into equations (7c) to (7e) leads to

L11(w)+L12(φx)+L13(φy)+S(w,f)+q=I12wt2+2ɛI1wt+I5¯¯3φxt2x+I5*¯¯3φyt2y+I7¯¯4wt2x2+I7*¯¯4wt2y2(21a)
L21(w)+L22(φx)+L23(φy)+L24(f)=I3¯¯2φxt2-I5¯¯3wt2x(21b)
L31(w)+L32(φx)+L33(φy)+L34(f)=I3*¯¯2φyt2-I5*¯¯3wt2y(21c)
in which
L11(w)=D112wx2+D112wy2+D124wx4+D134wx2y2+D144wy4-k1w+k2(2wx2+2wy2),L12(φx)=D11φxx+D153φxx3+D163φxxy2,L13(φy)=D12φyy+D173φyy3+D183φxx2y,S(w,f)=2fy22wx2-22fxy2wy2+2fx22wy2+1R2fy2,L21(w)=D21wx+D223wx3+D233wxy2,L22(φx)=D21φx+D242φxx2+D252φxy2,L23(φy)=D262φyxy,L24(f)=D273fx3+D283fxy2,L31(w)=D21wy+D313wx2y+D323wx3,L32(φx)=D332φxxy,L33(φy)=D21φy+D342φyx2+D352φyy2,L34(f)=D363fx2y+B373fy3(22)
and the detail of coefficients D1i(i=1,8¯),D2j(j=1,8¯),D3k(k=1,7¯) may be found in Appendix 2.

Introducing equation (17) into equation (16) gives the compatibility equation of the eccentrically stiffened FGM cylindrical panel as

I21*4fx4+I11*4fy4+J14fx2y2+J23φxx3+J33φxxy2+J43φyy3+J53φyyx2-c1I25*4wx4-c1I16*4wy4+J64wx2y2-(2wxy2-2wx22wy2-1R2wy2)=0(23)
where
J1=I31*-2I12*,J2=I23*-c1I25*,J3=I13*-c1I15*-I32*+c1I33*,J4=I14*-c1I16*,J5=I24*-c1I26*-I32*+c1I33*,J6=-c1I15*-c1I26*+2c1I33*(24)

Equations (21) and (23) are nonlinear equations in terms of variables w,φx,φy and f and used to investigate the nonlinear vibration and dynamic response of FGM thick eccentrically stiffened cylindrical panels in thermal environments using the Reddy’s third-order shear deformation shell theory.

In the present study, four edges of the FGM cylindrical panel are assumed to be simply supported and freely movable. The associated boundary conditions are

w=Nxy=φy=Mx=Px=0,Nx=Nx0atx=0,aw=Nxy=φx=My=Py=0,Ny=Ny0aty=0,b(25)
where Nx0,Ny0 are in-plane compressive loads at movable edges.

The mentioned condition (25) can be satisfied identically if the approximate solutions are represented by [8,21,32]

[w(x,y,t)φx(x,y,t)φy(x,y,t)]=[W(t)sinλmxsinδnyΦx(t)cosλmxsinδnyΦy(t)sinλmxcosδny](26)
where λm=mπ/a, δn=nπ/b; m,n are odd natural numbers representing the number of half waves in the x and y directions, respectively, and W(t),Φx,Φy are the time-dependent amplitudes.

By introducing equation (26) into the compatibility equation (23), we define the stress function as

f=A1(t)cos2λmx+A2(t)cos2δny+A3(t)sinλmxsinδny+12Nx0y2+12Ny0x2(27)
with
A1=δn232I21*λm2W(W+2μh),A2=λm232I11*δn2W(W+2μh),A3=Q1W+Q2Φx+Q3Φy(28)
and
Q1=λm2R+c1I25*λm4+c1I16*δn4-J6λm2δn2I21*λm4+J1λm2δn2+I11*δn4,Q2=-(J2λm3+J3λmδn2)I21*λm4+J1λm2δn2+I11*δn4,Q3=-(J4δn3+J5λm2δn)I21*λm4+J1λm2δn2+I11*δn4(29)

Replacing equations (26) and (27) into equation (21) and then applying Galerkin method to the resulting equations yields

l11W+l12Φx+l13Φy+l14(W+μh)Φx+l15(W+μh)Φy+[n1-(Nx0λm2+Ny0δn2)](W+μh)+n2W(W+μh)+n3W(W+2μh)+n4W(W+μh)(W+2μh)+n5q+n5Ny0R=I02Wt2+2ɛI1Wt-λmI5¯¯2Φxt2-δnI5*¯¯2Φyt2(30a)
l21W+l22Φx+l23Φy+n6(W+μh)+n7W(W+2μh)=I3¯¯2Φxt2-λmI5¯¯2Wt2(30b)
l31W+l32Φx+l33Φy+n8(W+μh)+n9W(W+2μh)=I3¯¯2Φyt2-δnI5¯¯2Wt2(30c)
in which the detail of coefficients l1i(i=1,3¯),ljk(j=2,3¯,k=1,2¯),nm(m=1,9¯) may be found in Appendix 3.

This is the basic equation to determine the nonlinear vibration of an eccentrically stiffened thick FGM cylindrical panel on elastic foundations in thermal environments.

Nonlinear dynamic response

Consider an eccentrically stiffened FGM cylindrical panel acted on by an uniformly distributed transverse load q=QsinΩt (Q is the amplitude of uniformly excited load, Ω is the frequency of the load). The system equation (30) has the form as

l11W+l12Φx+l13Φy+l14(W+μh)Φx+l15(W+μh)Φy+[n1-(Nx0λm2+Ny0δn2)](W+μh)+n2W(W+μh)+n3W(W+2μh)+n4W(W+μh)(W+2μh)+n5QsinΩt+n5Ny0R=I02Wt2+2ɛI1Wt-λmI5¯¯2Φxt2-δnI5*¯¯2Φyt2(31a)
l21W+l22Φx+l23Φy+n6(W+μh)+n7W(W+2μh)=I3¯¯2Φxt2-λmI5¯¯2Wt2(31b)
l31W+l32Φx+l33Φy+n8(W+μh)+n9W(W+2μh)=I3¯¯2Φyt2-δnI5¯¯2Wt2(31c)

The nonlinear dynamic response of eccentrically stiffened FGM cylindrical panel can be obtained by solving this equation combined with initial conditions to be assumed as W(0)=0,dWdt(0)=0 by using the fourth-order Runge–Kutta method.

Natural frequencies

In the case of q=0, the natural frequencies of the eccentrically stiffened FGM cylindrical panel can be determined by solving the following equation

|l11+n1+I0ω2l12-λmI5¯¯ω2l13-δnI5*¯¯ω2l21+n6-λmI5¯¯ω2l22+I3¯¯ω2l23l31+n8-δnI5¯¯ω2l32l33++I3¯¯ω2|=0(32)

Three angular frequencies of the FGM cylindrical panel in the axial, circumferential and radial directions are determined by solving equation (30) and the smallest one is being considered.

Validation

Firstly, Table 1 shows the comparison of the natural frequencies ω(Hz) for simply supported Al2O3/Ti-6Al-4 V square FGM plate (a=b=0.4m,h=0.005m) for the two special cases of isotropy in this paper with the results presented by Yang and Shen [16], Huang and Shen [17] based on Reddy’s higher order shear deformation plate theory and He et al. [20] using the classical laminated plate theory. The material properties are given as: Em=105.7GPa,νm=0.2981,ρm=4429kg/m3 for Ti-6Al-4 V and Ec=320.24GPa,νc=0.26,ρc=3750kg/m3 for Al2O3. From Table 1, good agreements are observed for the present solution and the numerical results of three publications.

Table

Table 1. Comparison of natural frequencies for isotropic square FGM plates (a=b=0.4m,h=0.005m).

Table 1. Comparison of natural frequencies for isotropic square FGM plates (a=b=0.4m,h=0.005m).

Table

Table 2. Material properties of the constituent materials of the considered FGM panel.

Table 2. Material properties of the constituent materials of the considered FGM panel.

Secondly, Figure 3 compares the dynamic response of simply supported eccentrically stiffened S-FGM cylindrical panel with metal–ceramic–metal layers under uniform external pressure using the results of Quan et al. [8] in the case of temperature-independent properties. Again, this comparison study also shows that the present result agrees very well the existing result.


                        figure

Figure 3. Comparisons of nonlinear dynamic response of eccentrically stiffened S-FGM cylindrical panel under uniform external pressure with the results of Quan et al. [8].

Next, we will investigate the effects of the volume fraction index, the geometrical dimensions, elastic foundations, imperfections and stiffeners on the nonlinear response of the eccentrically stiffened FGM cylindrical panel.

The effective material properties with dependent temperature in equation (5) are listed in Table 2, and the Poisson’s ratio is v=0.3 [21,28,3234].

The parameters for the stiffeners are [29]

s1=s2=0.4m,z1=z2=0.0225m,h1=h2=0.003m,d1=d2=0.004m(33)

Natural frequencies

Table 3 shows the influences of elastic foundations with two coefficients k1 and k2, geometrical parameter b/h, volume fraction index N and stiffeners on the natural frequencies of simply supported FGM cylindrical panel in thermal environments. Obviously, the value of the natural oscillation frequency increases when the values k1 and k2 increase. Furthermore, the Pasternak elastic foundation influences on the natural oscillation frequency larger than the Winkler foundation. The results of Table 3 also shows that decrease of volume fraction index N and ratio b/h leads to increase of the natural oscillation frequencies of FGM cylindrical panel. Moreover, the natural frequency of eccentrically stiffened FGM cylindrical panel is higher than one of FGM cylindrical panels without stiffeners with the same geometrical and material parameters.

Table

Table 3. Effects of elastic foundations, ratio b/h, volume fraction index N and stiffeners on natural frequencies (s-1) of FGM cylindrical panel in thermal environments.

Table 3. Effects of elastic foundations, ratio b/h, volume fraction index N and stiffeners on natural frequencies (s-1) of FGM cylindrical panel in thermal environments.

Effects of volume fraction index

Figure 4 considers the influences of volume fraction index N on the nonlinear dynamic response of eccentrically stiffened FGM cylindrical panels with movable edges in thermal environments. It is clear that the fluctuation amplitude of the eccentrically stiffened FGM cylindrical panel increases when the volume fraction index N decreases. This is reasonable because from equation (1), it is easy to see that when the volume fraction index N decreases, the volume fraction of ceramic decreases (the volume fraction of metal increases). Furthermore, the elastic module of ceramic is higher than metal (Ec>Em). Therefore, a decrease of volume fraction index N leads to a decrease of the elastic module of FGM, which is the reason for an increase of the fluctuation amplitude of the eccentrically stiffened FGM cylindrical panel.


                        figure

Figure 4. Effect of volume fraction index N on the nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments.

The ceramic volume fraction is decreased; moreover, elastic module of ceramic is higher than metal (Ec>Em).

Effects of temperature increment

Effect of temperature increment, ΔT, on the nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments is shown in Figure 5. The result from this figure shows that the fluctuation amplitude of eccentrically stiffened FGM cylindrical panel increases when temperature increment ΔT increases.


                        figure

Figure 5. Effects of temperature increment on the nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments.

Effects of elastic foundations

Figures 6 and 7 indicate the effect of elastic foundations on the nonlinear dynamic response of eccentrically stiffened FGM cylindrical panel in thermal environments with a/b=1,a/h=20,N=1,Px=0,Py=0. As can be seen, elastic foundations have beneficial effects on the nonlinear dynamic response of eccentrically stiffened FGM cylindrical panels. Especially, the panel fluctuation amplitude becomes considerably lower due to the support of elastic foundations. In addition, the beneficial effect of the Pasternak foundation with the module k2 on the dynamic response of the eccentrically stiffened FGM cylindrical panels is better than the Winkler one with the parameter k1.


                        figure

Figure 6. Effect of the linear Winkler foundation on nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments.


                        figure

Figure 7. Effect of the Pasternak foundation on nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments.

Effects of stiffeners

Figure 8 illustrates the effects of stiffeners on the nonlinear dynamic response of the FGM cylindrical panels with movable edges in thermal environments. Obviously, the fluctuation amplitude of the eccentrically stiffened FGM cylindrical panel is lower than the FGM cylindrical panel without stiffeners. In other words, the stiffeners strongly decrease the amplitude of the FGM cylindrical panels.


                        figure

Figure 8. Effects of stiffeners on the nonlinear dynamic response of FGM cylindrical panels in thermal environments.

Effects of force amplitude

Figure 9 indicates the effect of exciting force amplitude Q on nonlinear dynamic response of FGM cylindrical panels in thermal environments. Three values Q=25,500N/m2,Q=27,500N/m2 and Q=30,500N/m2 are used. It is easy to see that the nonlinear dynamic response amplitude of the FGM cylindrical panel increases when the value of the exciting force amplitude Q increases.


                        figure

Figure 9. Effects of exciting force amplitude on the nonlinear dynamic response of FGM cylindrical panels in thermal environments.

Effects of geometrical parameters

Figures 10 to 12 show the effects of geometrical parameters on nonlinear dynamic response of eccentrically stiffened FGM cylindrical panel in thermal environments with N=1,k1=0.3GPa/m,k2=0.01GPa.m,ΔT=300K,ɛ=0.1. Specifically, Figures 10 and 11 show the influences of ratios b/a and b/h on the nonlinear dynamic response of eccentrically stiffened FGM panels, which indicates that an increase of ratio b/a or b/h leads to an increase of the panel fluctuation amplitude. The nonlinear dynamic response of the eccentrically stiffened FGM panels with various values of ratio R/h is illustrated in Figure 12. As can be observed, the amplitude of the FGM shell increases when decreasing the ratio R/h.


                        figure

Figure 10. Effects of ratio b/a on the nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments.


                        figure

Figure 11. Effects of ratio b/h on the nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments.


                        figure

Figure 12. Effects of ratio R/h on the nonlinear dynamic response of the eccentrically stiffened FGM cylindrical panels in thermal environments.

This work studied the nonlinear vibration and dynamic response of eccentrically stiffened sandwich FGM thick cylindrical panels with metal–ceramic layers on Pasternalk elastic foundations in thermal environments. The sandwich FGM panel is reinforced by the eccentrically longitudinal and transversal stiffeners and subjected to mechanical load and temperature. The sandwich panel’s properties with metal–ceramic layers depend on temperature and change according to the nonlinear function of the thickness. Governing equations are derived using the Reddy’s third-order shear deformation shell theory and solved by the Galerkin method and fourth-order Runge-Kutta method. The present numerical results are compared with other results obtained from the literature. Some conclusions can be obtained from the numerical analysis in this study:

  1. The stiffener system strongly reduces the amplitude of the eccentrically stiffened sandwich FGM cylindrical panels.

  2. Temperature increment has a significant negative influence on the nonlinear vibration of the eccentrically stiffened sandwich FGM cylindrical panels.

  3. Elastic foundations have beneficial effect on the nonlinear dynamic response and natural frequencies of eccentrically stiffened sandwich FGM cylindrical panels. Moreover, the influence of nonlinear Pasternalk foundation is stronger than the linear Winkler foundation.

  4. The excitation force strongly influences the dynamic response of the eccentrically stiffened sandwich FGM cylindrical panels.

  5. The geometrical parameters have strong influences on the nonlinear vibration of the eccentrically stiffened sandwich FGM cylindrical panels.

Nguyen D Duc is also affiliated with National Research Laboratory, Department of Civil and Environmental Engineering, Sejong University, Korea.

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Grant in Mechanics code 107.02-2015.03 of the National Foundation for Science and Technology Development of Vietnam – NAFOSTED. The authors are grateful for this support.

1. Song, ZG, Zhang, LW, Liew, KM. Active vibration control of CNT-reinforced composite cylindrical shells via piezoelectric patches. Compos Struct 2016; 158: 92100.
Google Scholar | Crossref | ISI
2. Mehralian, F, Tadi Beni, Y, Ansari, R. On the size dependent buckling of anisotropic piezoelectric cylindrical shells under combined axial compression and lateral pressure. Int J Mech Sci 2016; 119: 155169.
Google Scholar | Crossref
3. Mirzaei, M, Kiani, Y. Free vibration of functionally graded carbon nanotube reinforced composite cylindrical panels. Compos Struct 2016; 142: 4556.
Google Scholar | Crossref | ISI
4. Sofiyev, AH, Karaca, Z, Zerin, Z. Non-linear vibration of composite orthotropic cylindrical shells on the non-linear elastic foundations within the shear deformation theory. Compos Struct 2017; 159: 5362.
Google Scholar | Crossref
5. Hosseini-Hashemi, S, Abaei, AR, Ilkhani, MR. Free vibrations of functionally graded viscoelastic cylindrical panel under various boundary conditions. Compos Struct 2015; 126: 115.
Google Scholar | Crossref
6. Norouzi, H, Alibeigloo, A. Three dimensional thermoviscoelastic analysis of a simply supported FGM cylindrical panel. Compos Struct 2016; 148: 181190.
Google Scholar | Crossref
7. Ahmed, MK . Buckling behavior of a radially loaded corrugated orthotropic thin-elliptic cylindrical shell on an elastic foundation. Thin-Walled Struct 2016; 107: 90100.
Google Scholar | Crossref
8. Quan, TQ, Tran, P, Tuan, ND, et al. Nonlinear dynamic analysis and vibration of shear deformable eccentrically stiffened S-FGM cylindrical panels with metal-ceramic-metal layers resting on elastic foundations. Compos Struct 2015; 126: 1633.
Google Scholar | Crossref | ISI
9. Chen, L, Liang, X, Yi, H. Vibro-acoustic characteristics of cylindrical shells with complex acoustic boundary conditions. Ocean Eng 2016; 126: 1221.
Google Scholar | Crossref
10. Shen, HS . Postbuckling of nanotube-reinforced composite cylindrical panels resting on elastic foundations subjected to lateral pressure in thermal environments. Eng Struct 2016; 122: 174183.
Google Scholar | Crossref | ISI
11. Duc, ND, Tuan, ND, Tran, P, et al. Nonlinear stability of eccentrically stiffened S-FGM elliptical cylindrical shells in thermal environment. Thin-Walled Struct 2016; 108: 280290.
Google Scholar | Crossref
12. Zhou, Y, Zhu, J. Vibration and bending analysis of multiferroic rectangular plates using third-order shear deformation theory. Compos Struct 2016; 153: 712723.
Google Scholar | Crossref
13. Sayyaadi, H, Rahnama, F, Askari Farsangi, MA. Energy harvesting via shallow cylindrical and spherical piezoelectric panels using higher order shear deformation theory. Compos Struct 2016; 147: 155167.
Google Scholar | Crossref
14. Selim, BA, Zhang, LW, Liew, KM. Vibration analysis of CNT reinforced functionally graded composite plates in a thermal environment based on Reddy’s higher-order shear deformation theory. Compos Struct 2015; 156: 276290.
Google Scholar | Crossref
15. Wattanasakulpong, N, Prusty, GB, Kelly, DW. Free and forced vibration analysis using improved third-order shear deformation theory for functionally graded plates under high temperature loading. J Sandw Struct Mater 2013; 15: 583606.
Google Scholar | SAGE Journals | ISI
16. Yang, J, Shen, HS. Vibration characteristics and transient response of shear-deformable functionally graded plates in thermal environments. J Sound Vibr 2002; 255: 579602.
Google Scholar | Crossref | ISI
17. Huang, XL, Shen, HS. Nonlinear vibration and dynamic response of functionally graded plates in thermal environments. Int J Solids Struct 2004; 41: 24032427.
Google Scholar | Crossref | ISI
18. Amirpour, M, Das, R, Saavedra Flores, EI. Analytical solutions for elastic deformation of functionally graded thick plates with in-plane stiffness variation using higher order shear deformation theory. Compos Part B: Eng 2016; 94: 109121.
Google Scholar | Crossref
19. Groh, RMJ, Weaver, PM. Static inconsistencies in certain axiomatic higher-order shear deformation theories for beams, plates and shells. Compos Struct 2015; 120: 231245.
Google Scholar | Crossref | ISI
20. He, XQ, Ng, TY, Sivashanker, S, et al. Active control of FGM plates with integrated piezoelectric sensors and actuators. Int J Solids Struct 2001; 38: 16411655.
Google Scholar | Crossref | ISI
21. Quan, TQ, Duc, ND. Nonlinear thermal stability of eccentrically stiffened FGM double curved shallow shells. J Therm Stress 2016; 5739: 126.
Google Scholar
22. Shao, D, Hu, S, Wang, Q, et al. Free vibration of refined higher-order shear deformation composite laminated beams with general boundary conditions. Compos Part B: Eng 2017; 108: 7590.
Google Scholar | Crossref
23. Yang, Y, Wang, Y, Fu, F. Effect of reinforcement stiffeners on square concrete-filled steel tubular columns subjected to axial compressive load. Thin-Walled Struct 2014; 82: 132144.
Google Scholar | Crossref
24. Xu, Y, Tong, Y, Liu, M, et al. A new effective smeared stiffener method for global buckling analysis of grid stiffened composite panels. Compos Struct 2016; 158: 8391.
Google Scholar | Crossref
25. Turner, KM, Vizzini, AJ. Response of impacted sandwich panels with integral stiffeners. J Sandw Struct Mater 2004; 6: 313326.
Google Scholar | SAGE Journals | ISI
26. Tong, JZ, Guo, YL. Elastic buckling behavior of steel trapezoidal corrugated shear walls with vertical stiffeners. Thin-Walled Struct 2015; 95: 3139.
Google Scholar | Crossref
27. Feng, Y, Zhang, H, Tan, X, et al. Effect of impact damage positions on the buckling and post-buckling behaviors of stiffened composite panel. Compos Struct 2016; 155: 184196.
Google Scholar | Crossref
28. Duc, ND, Tuan, ND, Quan, TQ, et al. Nonlinear mechanical, thermal and thermo-mechanical postbuckling of imperfect eccentrically stiffened thin FGM cylindrical panels on elastic foundations. Thin-Walled Struct 2015; 96: 155168.
Google Scholar | Crossref
29. Farahani, H, Azarafza, R, Barati, F. Mechanical buckling of a functionally graded cylindrical shell with axial and circumferential stiffeners using the third-order shear deformation theory. Comptes Rendus – Mecanique 2014; 342: 501512.
Google Scholar | Crossref
30. Shen, C, Xin, FX, Lu, TJ. Sound transmission across composite laminate sandwiches: influence of orthogonal stiffeners and laminate layup. Compos Struct 2016; 143: 310316.
Google Scholar | Crossref
31. Duc, ND . Nonlinear thermo-electro-mechanical dynamic response of shear deformable piezoelectric sigmoid functionally graded sandwich circular cylindrical shells on elastic foundations. J Sandw Struct Mater 2016. doi:10.1177/1099636216653266.
Google Scholar
32. Duc, ND . Nonlinear static and dynamic stability of functionally graded plates and shells, Hanoi, Vietnam: Vietnam National University Press, 2014.
Google Scholar
33. Duc, ND, Quan, TQ. Nonlinear dynamic analysis of imperfect FGM double curved thin shallow shells with temperature-dependent properties on elastic foundation. J Vibr Control 2015; 21: 13401362.
Google Scholar | SAGE Journals | ISI
34. Duc, ND, Quan, TQ. Transient responses of functionally graded double curved shallow shells with temperature-dependent material properties in thermal environment. Eur J Mech – A: Solids 2014; 47: 101123.
Google Scholar | Crossref

I11=E11-v2+E0AxTsxT,I12=vE11-v2,I13=E21-v2+E0AxTzxTsxT,I14=vE21-v2,I15=E41-v2+E0AxT(zxT)3sxT+dxT(hxT)3E0zxT4sxT,I16=vE41-v2,I17=11-ν,I18=11-2ν0dxTsxT,I22=E11-v2+E0AyTsyT,B24=E21-v2+E0A2Tz2Ts2T,I26=E41-v2+E0AyT(zyT)3syT+dyT(hyT)3E0zyT4syT,I28=11-2ν0dyTsyT,I31=E12(1+ν),I32=E22(1+ν),I33=E42(1+ν),I43=E31-v2+E0AxT(zxT)2sxT+E0(hxT)3dxT12sxT,I44=vE31-v2,I45=E51-v2+E0AxT(zxT)4sxT+E0dxT(hxT)3(zxT)22sxT+E0dxT(hxT)580sxT,I46=vE51-v2,I54=E31-v2+E0AyT(zyT)2syT+E0(hyT)3dyT12syT,I56=E51-v2+E0AyT(zyT)4syT+E0dyT(hyT)3(zyT)22syT+E0dyT(hyT)580syT,I62=E32(1+ν),I63=E52(1+ν),I71=E41-v2+E0AxT(zxT)3+E0zxTdxT(hxT)34,I73=E51-v2+E0AxT(zxT)4+E0dxT(zxT)2(hxT)32+E0dxT(hxT)580,I75=E71-v2+E0AxT(zxT)6+E0dxT(hxT)7448+15E0dxT(zxT)4(hxT)312+15E0dxT(zxT)2(hxT)580,I76=vE71-v2,I82=E41-v2+E0AyT(zyT)3+E0zyTdyT(hyT)34,I84=E51-v2+E0AyT(zyT)4+E0dyT(zyT)2(hyT)32+E0dyT(hyT)580,I86=E71-v2+E0AyT(zyT)6+E0dyT(hyT)7448+15E0dyT(zyT)4(hyT)312+15E0dyT(zyT)2(hyT)580,I93=E72(1+ν),(E1,E2,E3,E4,E5,E7)=-h/2h/2(1,z,z2,z3,z4,z6)E(z)dz,(Φ1,Φ2,Φ4)=11-v-h/2h/2(1,z,z3)E(z)α(z)ΔTdz,(Φ1is,Φ2is,Φ4is)=11-2ν0-h/2-h1-h/2(1,z,z3)E0α0ΔTdiTsiTdz,i=x,y

Appendix 2

Appendix 3

l11=-k1-k2(λm2+δn2)+D12λm4+D13λm2δn2+D14δn4+D19Q1λm4+D110Q1λm2δn2+D111Q1δn4-Q1λm2R,l12=-D11λm+D15λm3+D16λmδn2+D19Q2λm4+D110Q2λm2δn2+D111Q2δn4-Q2λm2R,l13=-D11δn+D17δn3+D18λm2δn+D19Q3λm4+D110Q3λm2δn2+D111Q3δn4-Q3λm2R,l14=32Q2λmδn3ab,l15=32Q3λmδn3ab,n1=-D11(λm2+δn2),n2=32Q1λmδn3ab,n3=2δn3abI21*λmR-8D19λmδn3abI21*-8D111λmδn3abI11*,n4=-λm416I11*-δn416I21*,n5=16mnπ2,l21=-λm3(D22+Q1D27)-λmδn2(D23+Q1D28),l22=D21-D24λm2-D25δn2-D27Q2λm3-D28Q2λmδn2,l23=-D26λmδn-D27Q3λm3-D28Q3λmδn2,n6=D21λm,n7=8D27δn3abI21*,l31=-δn3(D32+Q1D37)-λm2δn(D31+Q1D36),l32=-D33λmδn-D37Q2δn3-D36Q2λm2δn,l33=D21-D34λm2-D35δn2-D37Q3δn3-D36Q3λm2δn,n8=D21δn,n9=8D37λm3abI11*.

Cookies Notification

This site uses cookies. By continuing to browse the site you are agreeing to our use of cookies. Find out more.
Top