Nonlinear vibration of porous funcationally graded cylindrical panel using reddy’s high order shear deformation

The study results for dynamic response of PFGCP present the effect of geometrical ratio, elastic foundations: Winkler foundation and Paskternak foundation; loads: mechanical load and thermal load; and the material properties and distribution type of porosity. The results are shown numerically and are determined by using Galerkin methods and Fourth-order Runge-Kutta method.

VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 Original Article Nonlinear Vibration of Porous Funcationally Graded Cylindrical Panel Using Reddy’s High Order Shear Deformation Vu Minh Anh1, Nguyen Dinh Duc1,2,* Department of Construction and Transportation Engineering, VNU University of Engineering and Technology, 144 Xuan Thuy, Cau Giay, Hanoi, Vietnam Infrastructure Engineering Program, VNU Vietnam-Japan University (VJU), My Dinh 1, Tu Liem, Hanoi, Vietnam Received 02 December 2019 Accepted 06 December 2019 Abstract: The nonlinear dynamic response and vibration of porous functionally graded cylindrical panel (PFGCP) subjected to the thermal load, mechanical load and resting on elastic foundations are determined by an analytical approach as the Reddy’s third order shear deformation theory, Ahry’s function… The study results for dynamic response of PFGCP present the effect of geometrical ratio, elastic foundations: Winkler foundation and Paskternak foundation; loads: mechanical load and thermal load; and the material properties and distribution type of porosity The results are shown numerically and are determined by using Galerkin methods and Fourth-order Runge-Kutta method Keywords: Nonlinear dynamic response, porous functionally graded cylindrical panel, the high order shear deformation theory, mechanical load, thermal load, nonlinear vibration Introduction With the requirements of working ability such as bearing, high temperature in the harsh environment of some key industries such as defense industries, aircraft, space vehicles, reactor vessels and other engineering structures, in the world, many advanced materials have appeared, including Functionally Graded Materials (FGMs) FGMs is a composites material and is made by a combination of two main materials: metal and ceramic Therefore, the material properties of FGMs will have all the Corresponding author Email address: ducnd@vnu.edu.vn https//doi.org/ 10.25073/2588-1124/vnumap.4444 V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 outstanding properties of the two-component materials and vary with the thickness of the structure With these outstanding properties, FGMs have attracted the attention of scientists in the world In recent years, a lot of research carried out on FGMs, dynamic FGMs Contact mechanics of two elastic spheres reinforced by functionally graded materials (FGM) thin coatings are studied by Chen and Yue [1] Li et al [2] determined nonlinear structural stability performance of pressurized thin-walled FGM arches under temperature variation field Dastjerdi and Akgöz [3] presented new static and dynamic analyses of macro and nano FGM plates using exact three-dimensional elasticity in thermal environment Wang and Zu [4] researched about nonlinear dynamic thermoelastic response of rectangular FGM plates with longitudinal velocity Dynamic response of an FGM cylindrical shell under moving loads is investigated by Sofiyev [5] Wang and Shen [6] published nonlinear dynamic response of sandwich plates with FGM face sheets resting on elastic foundations in thermal environments Shariyat [7] gave vibration and dynamic buckling control of imperfect hybrid FGM plates with temperature-dependent material properties subjected to thermo-electro-mechanical loading conditions Nonlinear dynamic analysis of sandwich S- FGM plate resting on pasternak foundation under thermal environment are reseached by Singh and Harsha [8] Nonlinear dynamic buckling of the imperfect orthotropic E- FGM circular cylindrical shells subjected to the longitudinal constant velocity are studied by Gao [9] Reddy and Chin [10] studied thermomechanical analysis of functionally graded cylinders and plates Babaei et al [11] investigated thermal buckling and post-buckling analysis of geometrically imperfect FGM clamped tubes on nonlinear elastic foundation 3D graphical dynamic responses of FGM plates on Pasternak elastic foundation based on quasi-3D shear and normal deformation theory is presented by Han et al [12] Ghiasian et al [13] researched about dynamic buckling of suddenly heated or compressed FGM beams resting on nonlinear elastic foundation Geometrically nonlinear rapid surface heating of temperature-dependent FGM arches are gave by Javani et al [14] Shariyat [15] presented dynamic thermal buckling of suddenly heated temperature-dependent FGM cylindrical shells, under combined axial compression and external pressure Babaei et al [16] published large amplitude free vibration analysis of shear deformable FGM shallow arches on nonlinear elastic foundation In fact, porous materials appear around us and play in many areas of life such as fluid filtration, insulation, vibration dampening, and sound absorption In addition, porous materials have high rigidity leading to the ability to work well in harsh environments In recent years, porous materials have been researched and applied along with other materials to create new materials with the preeminent properties of component materials Porous functionally graded (PFG) is one of the outstanding materials among them and has research works such as: Gao et al [17] researched about dynamic characteristics of functionally graded porous beams with interval material properties Dual-functional porous copper films modulated via dynamic hydrogen bubble template for in situ SERS monitoring electrocatalytic reaction are proposed by Yang et al [18] Foroutan et al [19] investigated nonlinear static and dynamic hygrothermal buckling analysis of imperfect functionally graded porous cylindrical shells Li [20] presented nonlinear vibration and dynamic buckling analyses of sandwich functionally graded porous plate with graphene platelet reinforcement resting on Winkler–Pasternak elastic foundation Li [21] carried out experimental research on dynamic mechanical properties of metal tailings porous concrete Transient [22] published response of porous FG nanoplates subjected to various pulse loads based on nonlocal stress-strain gradient theory Esmaeilzadeh and Kadkhodayan [23] researched about dynamic analysis of stiffened bi-directional functionally graded plates with porosities under a moving load by dynamic relaxation method with kinetic damping Vibration analysis of magneto-electro-elastic heterogeneous porous material plates resting on elastic foundations are studied by Ebrahimi et al [24] Arshid and Khorshidvand [25] presented free vibration analysis of saturated porous FG circular plates integrated with piezoelectric actuators via differential quadrature V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 method Micro temperature-dependent FG porous plate: Free vibration and thermal buckling analysis using modified couple stress theory with CPT and FSDT are researched by Shojaeefard [26] Demirhan and Taskin [27] investigated bending and free vibration analysis of Levy-type porous functionally graded plate using state space approach The proposed method used in this study is the third-order shear deformation theory and the effect of the thermal load has been applied in a number of case studies such as Zhang [28] published nonlinear bending analysis of FGM beams based on physical neutral surface and high order shear deformation theory Reddy and Chin [29] investigated thermo-mechanical analysis of functionally graded cylinders and plates Modeling and analysis of FGM rectangular plates based on physical neutral surface and high order shear deformation theory is presented by Zang [30] In [31], Thon and Bélanger presented EMAT design for minimum remnant thickness gauging using high order shear horizontal modes Dynamic analysis of composite sandwich plates with damping modelled using high order shear deformation theory are carried out by Meunier and Shenoi [32] Cong et al.[33] investigated nonlinear dynamic response of eccentrically stiffened FGM plate using Reddy’s TSDT in thermal environment Stability of variable thickness shear deformable plates—first order and high order analyses are gave by Shufrin and M Eisenberger [34] In [35], Stojanović et al showed exact closed-form solutions for the natural frequencies and stability of elastically connected multiple beam system using Timoshenko and high order shear deformation theory A high order shear element for nonlinear vibration analysis of composite layered plates and shells are proposal by Attia and El-Zafrany [36] Khoa et al [37] observed nonlinear dynamic response and vibration of functionally graded nanocomposite cylindrical panel reinforced by carbon nanotubes in thermal environment In [38], Allahyari and Asgari found that thermo-mechanical vibration of double-layer graphene nanosheets in elastic medium considering surface effects; developing a nonlocal third order shear deformation theory Wang and Shi [39] presented a refined laminated plate theory accounting for the third order shear deformation and interlaminar transverse stress continuity From the literature review, the authors often used the high order shear deformation theory to investigate the nonlinear static, nonlinear dynamic or nonlinear vibration of FGM or plate porous funcationally graded (PFG) For the nonlinear dynamic and vibration of PFGCP has not carried out Therefore, in order to observe the nonlinear dynamic and vibration of PFGCP under mechanical load and thermal load, using the Reddy’s high order shear deformation theory and Ahry’ function are proposaled in this paper The natural frequency of PFGCP is obtained by using cylindrical panel fourthorder Runge-Kutta method Besides, the effect of geometrical ratio, elastic foundations: Winkler foundation and Paskternak foundation, the material properties and distribution type of porous on the modeling will be shown Theoretical formulation Fig.1 show a PFGCP resting on elastic foundations included Winkler foundation and Pasternak foundation in a Cartesian coordinate system x, y, z , with xy - the midplane of the panel z - the thickness coordinator, h /  z  h / a - the length b - the width h - the thickness of the panel R - the radius of the cylinderical panel V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 Fig Geometry of the PFGCP on elastic foundation The volume fractions of metal and ceramic, Vm and Vc , are assumed as [25] with using a simple power-law distribution:  2z  h  Vc ( z )    ;Vm ( z )   Vc ( z ),  2h  N (1) in which the volume fraction index N , subscripts m and c stand for the metal and ceramic constituents, respectively The material properties of PFGCP with prosity distribution are given as Porosity – I:  E ( z ), ( z),  ( z)   Em , m , m    2z  h    Ecm , cm , cm      Em  Ec , m   c ,  m  c  ,  2h  N (2) Porosity – II:  E ( z ), ( z),  ( z)   Em , m , m   | z |  2z  h    Ecm , cm , cm     1    Em  Ec , m   c ,  m  c  , 2 h   2h  N (3) with Ecm  Ec  Em , cm  c   m , cm  c  m , Kcm  Kc  Km , (4) and in this paper, the Poisson ratio   z  can be considered constant  ( z )  v b Porosity - II a Porosity - I Fig The porosity distribution type V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 2.1 Governing equations By using the Reddy’s high order shear deformation theory (HSDT), the governing equations and the dynamic analysis of PFGCP are determined The train-displacement relations taking into account the Von Karman nonlinear terms as [10,29,33]:  k x3    x    x0   k x1            xz    xz   k xz      z k  z k ;   z   ,  y         y  y  y yz  yz     k yz       k   k xy   xy   xy   xy    (5) With  u w  w 2   x         x R  x     x   x     kx    y    v  w       y    y   y   ;  k y    y      xy0     k xy        u  v  w w   x  y  y x x y  x  y       ;       x  w    x x    k x3     y  w  3  ,  k y   c1  y y    k xy3        w  x  y   x xy   y    ;    3c1  ,   xz     w   k xz     w      x x   k   x x   yz     yz       w     w   y y   y y      in which c1  / 3h2 ,  x ,  y  x ,  y - normal strains,  xy - the in-plane shear strain,  xz ,  yz - the transverse shear deformations u, v, w - the displacement components along the x, y, z directions, respectively x , y - the slope rotations in the  x, z  and  y, z  planes, respectively The strains are related in the compatibility equation [27, 28]: (6) V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 2  2 x0   y   xy   w   w  w  w       y x xy  xy  x y R x 2 (7) The effect of temperature will be descirbed in Hooke's law for a PFGCP as:  E [( x ,  y )    y ,  x   (1   )T (1,1)],   xy , xz , yz   2(1E )  xy ,  xz ,  yz  , x , y   (8) in which T is temperature rise The force and moment resultants of the PFGCP can be obtained with equations of stress components along with thickness of PFGCP as: h/2  Ni , M i , Pi     i 1, z, z  dz; i  x, y, xy, h/ (9) h/  Qi , Ri     j (1, z )dz; i  x, y; j  xz , yz h/ Replacing Eqs (3), (5) and (8) into Eq (9) gives:  N x , M x , Px     E1 , E2 , E4    x0   y0    E2 , E3 , E5   k x1   k 1y    ,    E4 , E5 , E7   k x3   k y3   (1   )   a ,  b ,  c     0 1  E1 , E2 , E4    y   x    E2 , E3 , E5   k y   k x   ,  N y , M y , Py     3   E , E , E k   k  (1   )  ,  ,        y x a b c    N xy , M xy , Pxy   2(11  )  E1 , E2 , E4   xy0   E2 , E3 , E5  k xy1   E4 , E5 , E7  k xy3  ,  E1 , E3   xz0   E3 , E5  k xz2  ,  Qx , Rx   2(1   )  Qy , Ry   2(11  )  E1 , E3   yz0   E3 , E5  k yz2  , (10) with h /2  E1 , E2 , E3 , E4 , E5 , E7    1, z, z , z , z , z  E ( z)dz,  h /2 h /2 ( a ,  b ,  c )   (1, z, z ) E ( z ) ( z )T ( z )dz,  h /2 and the coefficients Ei (i   5,7) are give in Appendix From Eq (10), The inverse expression are obtained: (11) V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21  x0   xy 1  N x   N y  E2 k x1  E4 k x3   a  ;  y0   N y   N x  E2 k 1y  E4 k y3   a  , E1 E1 1   2(1   ) N xy  E2 k xy  E4 k xy3  E1 (12) The equations of motion are [26] based on HSDT: N x N xy  2  2u 3 w   I  J1 2x  c1 I , x y t t xt  y  2v 3w   I  J1  c1 I , x y t t yt N xy N y  R Ry  Qx Qy 2w 2w 2w   c2  x   N  N  N  x xy y x y y  x xy y  x  2 P  Pxy  Py   2w 2w  N y  c1  2x     k1w  k2     q  x   x  y  y  x  y R      4w  w w 4w   I  2 I  c12 I  2  2  t t y t   x t    3u  3v  c1  I   yt   xt   3x  3 y   J    4 yt   xt (13)    ,    P Pxy  M x M xy  2x  2u 3w   c1  x   Q  c R  J  K  c1 J ,  x x 2 x y y  t t xt  x M xy M y  2 y  Pxy Py   2v 3 w   c1    Q  c R  J  K  c J ,  y y x y y  t t yt  x where h /2 Ii     z  z dz, i (i  0,1,2,3,4,6), (14)  h /2 J i  Ii  c1Ii  , K2  I  2c1I  c I , c2  3c1 , and the coefficients Ii (i   4,6) are noted in Appendix, and k1 - Winkler foundation modulus k - Pasternak foundation model q - an external pressure uniformly distributed on the surface of the plate  - damping coefficient The stress function f  x, y, t  is introduced as V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 Nx  2 f 2 f 2 f , N  , N   y xy y x xy (15) Replacing Eq (15) into the two first Eqs (13) yields J1  2x c1 I  w  2u    , t I t I xt (16a) J1   y c1 I  w  2v    t I t I yt (16b) By substituting Eqs (16) into the Eqs (13) and Eqs (13) can be rewritten  R Ry  Qx Qy 2w 2w 2w   c2  x   Ny   N x  N xy x y y  x xy y  x 2  2 P  Pxy  Py   2w 2w  N y  c1  2x     k1 w  k2     q  x xy y  y  R  x    w 2w w  c12 I 32 4w   I  2 I   c12 I  2  2   t t  I y t    x t (17a)  J1 I 3c1   3x  J1 I 3c1    y J c   J c  ,     I  xt  I  yt    3 w  P Pxy  M x M xy J12   2x  c1 I J1   c1  x   Q  c R  K    c J ,     x x 4 x y y  I  t  I  x   xt   3 w  Pxy Py  J12    y  c1 I J1   c1    c1 J    Qy  c2 Ry   K     x y y  I  t  x   I0  yt M xy M y (17b) (17c) Substituting Eqs (12) and (15) into the equation (7), we have:  4 f 4 f 4 f   2 w  2 w 2 w       E1  x x y y   xy  x y (18) The system of motion Eqs (17) is rewritten as Eq.(19) by replacing Eq (6) into Eq (10) and then into Eqs (17): 2 w w  2 I t t 2 4 c I   w J I c   x  J1 I 3c1    y  w      c12 I   2  2    J c1    J c  ,   y t   I  xt  I  yt  I0   x t L11  w   L12 x   L13  y   P  w, f   q  I   3 w J   2  c I J L21  w   L22 x   L23  y    K   2x    c1 J  , I  t   I0  xt (19a) (19b) V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21   3 w J   y  c I J L31  w   L32 x   L33  y    K      c1 J  , I  t   I0  yt (19c) and the Lij  i   3, j   3 and P are noted in Appendix In next section, in order to determined nonlinear dynamical analysis of PFGCP using the third order shear deformation theory, the Eqs (18)-(19) are used along with boundary conditions and initial conditions 2.2 Boundary conditions In order to determined nonlinear dynamical analysis of PFGCP, the edges of the PFGCP are simply supported and immovable (IM) and boundary conditions are [27, 28]: w  u   y  M x  Px  0, N x  N x at x  0, a, (20) w  v = x  M y  Px  0, N y  N y at y  0, b, in which N x , N y are forces along the x and y axis Wiht the boundary condition (20) and In order to solve Eqs (18) and (19), the approximate solutions can be written as [26]: w  x, y, t   W  t  sin  x sin  y, x  x, y, t    x  t  cos  x sin  y, (21)  y  x, y, t    y  t  sin  xcos y, with m n  ,  and m, n  1,2, are the numbers of half waves in the direction x, y , W ,  x ,  y a b the amplitudes which are functions dependent on time The stress function is defined as: 1 f  x, y, t   A1  t  cos 2 x  A2  t  cos 2 y  N x y  N y x , (22a) 2 with E 2 E2 A1  W ; A2  W (22b) 32 32 2.3 Nonlinear dynamical analysis In ordet to obtain Eqs (23), Eqs (21) and (22) are replaced into the Eq (19) after that applying Galerkin method: 16 2 W W l11   N x 0 +N y    W  l12  x  l13 y  n1W  q  n  2 I 2   mn t t m   x n   y  2  2 , a t b t  2 x m  W l21W  l22  x  l23 y  1   , t a t (23a) (23b) V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 10 l31W  l32  x  l33 y  1 Where 1  K   2 y t  2 n  W , b t cI J J12 ; 2   c1 J , I0 I0   m 2  n   l11    E1  3c1 E3  c2  E3  3c1 E5        a   b   2(1   )    2   m   n    m  c12   m   n    E42      k1  k2    +    E7         ,    a   b    E1   a   b    a  R  c  m   E2 E4 E42   m  E1  3c1 E3  c2  E3  3c1 E5     E   c E    1  2(1   ) a   a   E1 E1     E2 E4  E42   m  n       c1   E   c  E   1     ,  E1     E1   a  b        v   l12    c  n   E2 E4 E42   n  E1  3c1 E3  c2  E3  3c1 E5    E   c E   1   2(1   ) b   b   E1 E1     E2 E4  E42    m  n v    v c1   E   ,   c1      E7    E1     E1    a  b          l13   E1   m   n      16   a   b  n1     ,  2  I    m   n   n2  I  c12  I         , I    a   b    2 c1 m   m   n    c1 E42 EE  l21    E5      +    c1 E7   a   a   b    E1 E1  m   E1  3c1 E3   c2  E3  3c1 E5   , 2(1   ) a    m   E22 E2 E4   E42    n    l22     E   c E   c E    1  1     a  2(1   )  b    E1 E1  E1         E1  3c1 E3   c2  E3  3c1 E5   , 2(1   )  l23     E22 EE  E  mn   2c1  E5    c12  E7    ,  E3  1    ab  E1 E1  E1     (23c) V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 l31   11 2 c1 n  E42 E2 E4    m   n   n c E  c  E    E1  3c1 E3  c2  E3  3c1 E5  , 1       b  E1 E1    a   b   2(1   ) b    E22 E2 E4   E42   E   c E   c E    1  , 1 E1 E1  E1        n   E22 E2 E4   E42    m    l33     E   c E   c E    1  1     b  2(1   )  a    E1 E1  E1         E1  3c1 E3  c2  E3  3c1 E5   , 2(1   )  l32   mn 2 1  v  ab In order to investigate the effect of thermal load on the PFGCP The modeling are simply supported and immovable in all edges In other words, boundary conditions of case stud is u  (on x  0, a ) and v  (on y  0, b ), and and is satisfied equations as [10, 29]: u v dxdy  0,  dxdy   x y 0 0 With Eqs (6) and (12) along with (15), Ep (25) is determined as: b a a b (24) E cE u w  w   a   f, yy  vf, xx   x , x  x , x  w, xx     ,   x E1 E1 E1 R  x  E1 (25) E2 c1 E4 v 1  w   a   f , xx  vf , yy    y , y  y, y  w, yy    y   E y E1 E1 E1   In order to obtain Eqs (26), Eqs (21) and (22) are placed into Eq (25) then the obtained result are replaced into Eq (24) we have:   E2  c1 E4   x  v y   c1 E4   v W  N x0   a  2    v mn 1  v  (26a) E1 2    v   W , 1  v  N y0    a  E2  c1 E4   v x   y   c1 E4  v   W   2    v mn 1  v  E1 v   W  1  v  (26b) Putting Eqs (26) into the Eq (23) We obtain: é æ æ mp ö æ np ö ö F ự l11 + ỗ ỗ + ÷ a úW + l12 F x + l13F y + l14F x W + l15F y W + l16W + l17W ỗố ố a ữứ ỗố b ữứ ữứ 1- v ú ê ë û 2 16 ¶ W ¶W mp ¶2 F x np ¶ F y + q = n + e I + r + r , 2 ¶t a ¶t b ¶t mnp ¶t l21W  l22  x  l23 y  1  2 x m  W   , t a t (27a) (27b) V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 12 l31W  l32  x  l33 y  1  2 y t  2 n  W , b t (27c) with   m 2  n    E2  c1 E4     v     a   b   m  l14  , a mn 1  v  l15    m 2  n    E2  c1 E4    v       a   b   n  l16   mn 1  v  4c1 E4   2v     mn 1  v  b , , 2  m 4 E1  m   n   n   l17  n1     2v       , 1  v   a   a   b   b   The Eqs (27) are used to dertermine dynamic analysis of PFGCP with suported by elastic foundations along with the effect of temperature with the appearance of the coefficient  a In order to obtain the natural frequencies   of PFGCP, solving determinant of matrix (28) with linear coefficients taken from Eqs (27) along with q  The obtained results are three natural frequencies, the smallest value will be selected   m   n    a l11      n2   a   b    v   m l21    a n l31    b l12   m  a l13   n  b l22  1 l23 l32 l33  1  (28) Numerical results and discussion In this paper, consider PFGCP is under the influence of a uniformly distributed load on the surface with equations in the harmonic function q  Q sin t with Q is the amplitude of uniformly excited load Besides, the fourth-order Runge-Kutta method is used to solve equations (27) and the material properties of the component material are given as [40]: Ec  380  109 N / m2 , c  3800kg / m3 ,  c  7.4 10 6 o C 1 , K c  10.4W / mK , Em  70  109 N / m2 , m  2702kg / m3 ,  m  23 10 6 o C 1 , K m  204W / mK , v  0.3, V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 13 3.1 Natural frequency In order to verify the reliability of the used method, the numerical results of this paper will be compared with the published results of Duc et al [40] In the Duc author's paper, the used method is high order shear deformation for FGM plate In order to obtain the highest accuracy of the comparison results, the modeling of this paper was brought to the same format as the modeling in Duc’s paper by R = infinity The comparison results of the free natural frequency are shown in Table It can be seen that the results of both papers are not much different It verifies that the used method is reliable   Table Comparison of natural frequency s 1 of PFGCP with other paper with same conditions as a / b  1, a / h  20, N   k1; k2   GPa / m; GPa.m  0;0   0.25;0  0.25;0.02   T  o C  Natural frequency s 1 Duc et al [40] Present study 50 105 50 105 50 105 1941.73 1570.25 1015.19 2303.00 1999.79 1601.28 2779.46 2533.92 2232.86 1922 1522 1082 2341 1998 1687 2790 2576 2293 The natural frequencies  s 1  of PFGCP with the influence of index N , modes  m, n  along with o a / b  1, a / h  20, k1  0, k2  , T  100 C and using type of porosity distribution We can see that the natural frequencies will raise when modes  m, n  increase in which index N is constant In contrast, when index N increase in which modes  m, n  is constant, the natural frequencies increase Form equation (1), it can be explained that index N = corresponding to isotropic uniform panel is made from ceramic materials, N = is the case when the ceramic and metal components are distributed linearly over the thickness of the structure wall and when N increase, the volume ratio of the metal component in the structure increases   Table The natural frequencies s 1 of PFGCP with the influnece of index N , modes  m, n  along with o a / b  1, a / h  20, k1  0, k2  0, R / h  400 and T  100 C 1 N 2 3 4 5  m  1, n  1  m  1, n  3  m  1, n  5  m  3, n  5  m  5, n  5 0.5 2702 11175 27051 34564 48871 2434 9766 24060 30838 43767 2251 8817 21037 26808 37785 2172 8505 20218 25719 36137 12 2080 8174 19369 24594 34440 ∞ 1629 6195 14636 18610 26149 V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 14 Table shows the natural frequencies  s 1  of PFGCP with the influnece of ratio a / b , elastic a / h  20, N  1, T  100o C, distribution type along with (m, n)  (1,1), R / h  400 It can easily see that the natural frequencies increase when ratio a / b raises as well as when modeling is supported by elastic foundations Furthermore, the comparison of the effect between type of porosity distribution and type of porosity distribution are shown in table The effect of type of porosity distribution on the natural frequencies is smaller than the effect type of porosity distribution foundation, porosity   Table The natural frequencies s 1 of PFGCP with the influnece of ratio a / b , elastic foundation (m, n)  (1,1), R / h  400 a / h  20, N  1, T  100o C  k1; k2    Natural frequency s 1 a/b  GPa / m; GPa.m 1.5 1.5 1.5  0;0   0.2;0  0.2;0.04 Type Typ 2525 3757 5448 2757 3917 5559 3528 4816 6554 2434 3393 5068 2449 3569 5187 3292 4538 6242 3.2 Nonlinear dynamic response 10-6 Porosity - II Porosity - I W(m) -2 -4 (m,n)=(1,1), a/b=1, a/h=20, N = 1, k 1=0.1 GPa/m, -6 k2=0.02 GPa/m, q=1500sin(600t), -8 0.005 0.01 0.015 0.02 0.025 0.03 T = 100oC 0.035 0.04 0.045 0.05 t(s) Fig Comparision between two type of porosity distribution in PFGCP Figure shows the effect of porosity distribution on nonlinear dynamic response of the PFGCP with a / b  1, a / h  20, T  0,   0.1 It can easily see that the amplitude of deflection for PFGCP in case Porosity I bigger more than the amplitude of deflection for PFGCP in Porosity II In other words, the porosity distribition Porosity II will enhance the loading carrying capacity of PFGCP more than the Porosity I Thus, other results will use the Porosity II to investigate the effect of other factors on the nonlinear dynamic response of the panel V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 15 Figures and give the nonlinear dynamic response of PFGCP with the influence of ratio a / b , ratio a / h and supported by elastic foundations The amplitude of deflection of the PFGCP will raise when increasing ratio a / h or decrease a / b In figure 3, the amplitude of deflection has the highest value with a / b  When increasing a / b from to 2, the amplitude of deflection decreases about times In figure 4, the amplitude of deflection has the smallest value with a / h  40 When increasing a / h from 20 to 40, the amplitude of deflection raise about times 10-6 a/b = a/b = 1.5 a/b = W(m) -1 -2 -3 -4 (m,n) = (1,1), a/h=20, N=1, k =0.1 GPa/m; k = 0.02 GPa.m, -5 q = 1500sin(600t), 0.005 0.01 0.015 0.02 o T = 100 C 0.025 0.03 0.035 0.04 0.045 0.05 t(s) Fig The nonlinear dynamical response of the PFGCP with influence of effect of ratio a / b The influence of index N on the nonlinear dynamic response of the PFGCP is investigated in figure From figure can see that the amplitude vibration of the PFGCP will raise when the index N increase and change about 2.3 times This behavior can be explained that the index N increase results in a volume ratio of ceramic more bigger than a volume ratio of metal It leads to bigger change in Elastic modulus, in other words, it leads to decrease load carrying capacity of the PFGCP 10 -5 a/h = 20 a/h = 30 a/h = 40 W(m) -1 -2 -3 Porosity - II (m,n) = (1,1), a/b = 1, N=1, k =0.1 GPa/m; k2 = 0.02 GPa.m, q = 1500sin(600t), -4 0.005 0.01 0.015 0.02 0.025 0.03 o T = 100 C 0.035 0.04 0.045 0.05 t(s) Fig The nonlinear dynamical response of the PFGCP with influence of effect of ratio a / h 16 V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 10-5 N=1 N=2 N=3 0.8 0.6 0.4 W(m) 0.2 -0.2 -0.4 -0.6 Porosity - II (m,n)=(1,1), a/b=1, a/h=20, k =0.1 GPa/m, -0.8 k2 =0.02 GPa/m,q=1500sin(600t), -1 0.005 0.01 0.015 0.02 0.025 0.03 o T = 100 C 0.035 0.04 0.045 0.05 t(s) Fig Effect of power law index N on nonlinear dynamical response of the PFGCP The influence of elastic foundations on the nonlinear dynamic response of the PFGCP with a / b  1, a / h  20, N  1, T  100o C are shown in figures and It clear that the elastic foundations: Winkler foundation (figure 6) and Pasternak foundations (figure 7) has a positive influence on the amplitude vibrations When increasing the Winkler foundation, the Pasternak foundations, the amplitude vibrations will decrease But the influence of the Pasternak foundation better than Winkler foundations, the amplitude vibrations will enhance about 200% with the change of the Winkler foundation in which the amplitude vibrations will enhance about 250% In order to enhance the load carrying capacity of modeling, the elastic foundations have an important role 10-6 k 1=0.1 GPa/m k 1=0.3 GPa/m k 1=0.7 GPa/m W(m) -2 -4 Porosity - II: (m,n) = (1,1), a/b = 1, N=1, a/h = 20; k = 0.02 GPa.m, q = 1500sin(600t), T = 100 oC -6 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 t(s) Fig Effect of the linear Winkler foundation on nonlinear dynamical response of the PFGCP V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 10 17 -6 k 2=0.02 GPa.m k 2=0.04 GPa.m k 2=0.07 GPa.m W(m) -1 -2 -3 Porosity - II (m,n) = (1,1), a/b = 1, N=1, a/h = 20, -4 k1 =0.1 GPa/m; q = 1500sin(600t), -5 0.005 0.01 0.015 0.02 0.025 0.03 o T = 100 C 0.035 0.04 0.045 0.05 t(s) Fig Effect of the Pasternak foundation on nonlinear dynamical response of the PFGCP The influence of exciting force amplitude on nonlinear dynamic response of the PFGCP with Q  1500,1800,2500 N / m2  and T  100o C are proposed in figure The amplitude of vibration will raise when the excited force amplitude increase In other words, the exciting force amplitude decrease, the amplitude of vibration will decrease 10 -6 Q = 1500 N/m Q = 1800 N/m Q = 2500 N/m W(m) -2 -4 -6 Porosity - II (m,n) = (1,1), a/b = 1, N=1, a/h=20, k =0.1 GPa/m, -8 k = 0.02 GPa.m, q = Qsin(600t), 0.005 0.01 0.015 0.02 0.025 0.03 o = 100 C 0.035 0.04 0.045 0.05 t(s) Fig Nonlinear dynamic responses of the PFGCP with different loads The nonlinear response of the PFGCP with the effect of uniform temperature is presented in figure From figure 9, it can easily seen that the uniform temperature has a negative influence on the amplitude deflection of the cylindrical panel PFG The amplitude deflection of the PFGCP will increase when the uniform temperature rise It can be explained that the uniform temperature changes the stiffness matrix 18 V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 10-6 10 o T = 100 C Porosity - II (m,n) = (1,1), a/b = 1, N=1, a/h=20, k =0.1 GPa/m, k = 0.02 GPa.m, q = 1500sin(600t) o T = 200 C T = 300o C W(m) -2 -4 -6 -8 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 t(s) Fig Effect of uniform temperature rise on nonlinear response of the PFGCP Conclusions This is the first paper uses the high order shear deformation theory and Ahry function to investigate the nonlinear dynamic response and vibration of panel FGM with porosity distribution Besides, the PFGCP is supported by elastic foundations and under mechanical load, thermal load By using the Galerkin method as well as the four-order Runge-Kutta method, the numerical results for nonlinear dynamic response of modelings are described by figure and table with the influences of the materials properties, thermal load, mechanical load, elastic foundations, geometrical parameters, the porosity distribution type This paper obtained some remarkable results as: - The porosity distribution type will enhance the loading carrying capacity of PFGCP more than the porosity distribution type - The amplitude of deflection of the PFGCP will raise when increasing ratio a / h or decrease a / b - The amplitude vibration of the PFGCP will raise when the index N increase and changes about 2.3 times - When increasing the Winkler foundation, the Pasternak foundations, the amplitude vibrations will decrease The influence of the Pasternak foundation better than Winkler foundations - The amplitude of vibration will raise when the exciting force amplitude increase - The uniform temperature has a negative influence on the amplitude deflection of the cylindrical panel PFG The amplitude deflection of the PFGCP will increase when the uniform temperature rise Acknowledgements This work has been supported by VNU Hanoi - University of Engineering and Technology under the project number CN.19.04 The authors are grateful to this support V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 19 References [1] X.W Chen, Z.Q Yue, Contact mechanics of two elastic spheres reinforced by functionally graded materials (FGM) thin coatings, Engineering Analysis with Boundary Elements 109 (2019) 57-69 [2] Z Li, J Zheng, Q Sun, H He, Nonlinear structural stability performance of pressurized thin-walled FGM arches under temperature variation field, International Journal of Non-Linear Mechanics 113 (2019) 86 – 102 [3] S Dastjerdi, B Akgöz, New static and dynamic analyses of macro and nano FGM plates using exact threedimensional elasticity in thermal environment, Composite Structures 192 (2018) 626-641 [4] Y.Q Wang, W.Z Jean, Nonlinear dynamic thermoelastic response of rectangular FGM plates with longitudinal velocity, Composites Part B: Engineering 117 (2017) 74-88 [5] A.H Sofiyev, Dynamic response of an FGM cylindrical shell under moving loads, Composite Structures 93 (2010) 58-66 [6] Z.X Wang, H.S Shen, Nonlinear dynamic response of sandwich plates with FGM face sheets resting on elastic foundations in thermal environments, Ocean Engineering 57 (2013) 99-110 [7] M Shariyat, Vibration and dynamic buckling control of imperfect hybrid FGM plates with temperature-dependent material properties subjected to thermo-electro-mechanical loading conditions, Composite Structures 88 (2009) 240-252 [8] S.J Singh, S.P Harsha, Nonlinear dynamic analysis of sandwich S- FGM plate resting on pasternak foundation under thermal environment, European Journal of Mechanics - A/Solids 76 (2019) 155-179 [9] K Gao, W Gao, D Wu, C Song, Nonlinear dynamic buckling of the imperfect orthotropic E- FGM circular cylindrical shells subjected to the longitudinal constant velocity, International Journal of Mechanical Sciences 138– 139 (2018) 199-209 [10] J N Reddy, C.D Chin, Thermo-mechanical analysis of functionally graded cylinders and plates, J Thermal Stress; 21 (1998) 593-626 [11] H Babaei, Y Kiani, M.R Eslami, Thermal buckling and post-buckling analysis of geometrically imperfect FGM clamped tubes on nonlinear elastic foundation, Applied Mathematical Modelling 71 (2019) 12-30 [12] S.C Han, W.T Park, W.Y Jung, 3D graphical dynamic responses of FGM plates on Pasternak elastic foundation based on quasi-3D shear and normal deformation theory, Composites Part B: Engineering 95 (2016) 324-334 [13] S.E Ghiasian, Y Kiani, M.R Eslami, Dynamic buckling of suddenly heated or compressed FGM beams resting on nonlinear elastic foundation, Composite Structures 106 (2013) 225-234 [14] M Javani, Y Kiani and M.R Eslami, Geometrically nonlinear rapid surface heating of temperaturedependent FGM arches, Aerospace Science and Technology 90 (2019) 264-274 [15] M Shariyat, Dynamic thermal buckling of suddenly heated temperature-dependent FGM cylindrical shells, under combined axial compression and external pressure, International Journal of Solids and Structures 45 (2008) 25982612 [16] H Babaei, Y Kiani, M.Z Eslami, Large amplitude free vibration analysis of shear deformable FGM shallow arches on nonlinear elastic foundation, Thin-Walled Structures 144 (2019) 106237 [17] K Gao, R Li and J Yang, Dynamic characteristics of functionally graded porous beams with interval material properties, Engineering Structures 197 (2019) 109441 [18] H Yang, X Hao, J Tang, W Jin, J Hu, Dual-functional porous copper films modulated via dynamic hydrogen bubble template for in situ SERS monitoring electrocatalytic reaction Applied Surface Science 494 (2019) 731739 [19] K Foroutan, A Shaterzadeh, H Ahmadi, Nonlinear static and dynamic hygrothermal buckling analysis of imperfect functionally graded porous cylindrical shells, Applied Mathematical Modelling January 77 (2020) 539553 [20] Q Li, D Wu, X Chen, L Liu, W Gao, Nonlinear vibration and dynamic buckling analyses of sandwich functionally graded porous plate with graphene platelet reinforcement resting on Winkler–Pasternak elastic foundation, International Journal of Mechanical Sciences 148 (2018) 596-610 [21] H.N Li, P.F Liu, C Li, G Li, H Zhang, Experimental research on dynamic mechanical properties of metal tailings porous concrete, Construction and Building Materials 213 (2019) 20-31 20 V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 [22] S.S Mirjavadi, B M Afshari, M.R Barati, A.M.S Hamouda, Transient response of porous FG nanoplates subjected to various pulse loads based on nonlocal stress-strain gradient theory, European Journal of Mechanics A/Solids March–April 74 (2019) 210-220 [23] M Esmaeilzadeh, M Kadkhodayan, Dynamic analysis of stiffened bi-directional functionally graded plates with porosities under a moving load by dynamic relaxation method with kinetic damping, Aerospace Science and Technology 93 (2019) 105333 [24] F Ebrahimi, A Jafari, M.R Barati, Vibration analysis of magneto-electro-elastic heterogeneous porous material plates resting on elastic foundations, Thin-Walled Structures 119 (2017) 33-46 [25] E Arshid, A.R Khorshidvand, Free vibration analysis of saturated porous FG circular plates integrated with piezoelectric actuators via differential quadrature method, Thin-Walled Structures 125 (2018) 220-233 [26] M.H Shojaeefard, H.S Googarchin, M Ghadiri, M Mahinzare, Micro temperature-dependent FG porous plate: Free vibration and thermal buckling analysis using modified couple stress theory with CPT and FSDT, Applied Mathematical Modelling 50 (2017) 633-655 [27] P.A Demirhan, V Taskin, Bending and free vibration analysis of Levy-type porous functionally graded plate using state space approach, Composites Part B: Engineering 160 (2019) 661-676 [28] D.G Zhang, Nonlinear bending analysis of FGM beams based on physical neutral surface and high order shear deformation theory, Composite Structures 100 (2013) 121-126 [29] J.N Reddy, C.D Chin, Thermo-mechanical analysis of functionally graded cylinders and plates, J Thermal Stress 21 (1998) 593-626 [30] D.G Zhang, Modeling and analysis of FGM rectangular plates based on physical neutral surface and high order shear deformation theory International Journal of Mechanical Science 68 (2013) 92-104 [31] A Thon, P Bélanger, EMAT design for minimum remnant thickness gauging using high order shear horizontal modes Ultrasonics 95 (2019) 70-78 [32] M Meunier, R.A Shenoi, Dynamic analysis of composite sandwich plates with damping modelled using high order shear deformation theory Composite Structures 54 (2001) 243-254 [33] P.H Cong, V.M Anh, N.D Duc, Nonlinear dynamic response of eccentrically stiffened FGM plate using Reddy’s TSDT in thermal environment, Journal of Thermal Stresses 40 (2017) 704-732 [34] Shufrin, M Eisenberger, Stability of variable thickness shear deformable plates—first order and high order analyses, Thin-Walled Structures 43 (2005) 189-207 [35] V Stojanović, P Kozić, G Janevski, Exact closed-form solutions for the natural frequencies and stability of elastically connected multiple beam system using Timoshenko and high order shear deformation theory, Journal of Sound and Vibration 332 (2013) 563-576 [36] O Attia, A El-Zafrany, A high order shear element for nonlinear vibration analysis of composite layered plates and shells, International Journal of Mechanical Sciences 41 (1999) 461-486 [37] N.D Khoa, V.M Anh, N.D Duc, Nonlinear dynamic response and vibration of functionally graded nanocomposite cylindrical panel reinforced by carbon nanotubes in thermal environment, Journal of Sandwich Structures & Materials, 2019 https://doi.org/10.1177/1099636219847191 [38] E Allahyari, M Asgari, Thermo-mechanical vibration of double-layer graphene nanosheets in elastic medium considering surface effects; developing a nonlocal third order shear deformation theory, European Journal of Mechanics - A/Solids 75 (2019) 307-321 [39] X Wang, G Shi, A refined laminated plate theory accounting for the third order shear deformation and interlaminar transverse stress continuity, Applied Mathematical Modelling 39 (2015) 5659-5680 [40] N.D Duc, D.H Bich, P.H Cong, Nonlieanr thermal dynamic response of shear deformable FGM plates on elastic foundations Journal of Thermal Stresses 29 (2016) 278 – 297 V.M Anh, N.D Duc / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-21 Appendix L11  w    2 w 2 w   E1  3c1 E3  c2  E3  3c1 E5      2(1   ) y   x   2 w 2 w    k1 w  k2    , y   x    c  EE E2  E1  3c1 E3  c2  E3  3c1 E5   x   E5   c1  E7  L12 x   2(1   ) x    E1 E1   c12   w E2 4 w 4 w    2 +   E7   E1    x x y y      3x    x  E2 E4  E42    3x      c1   E   c  E  ,        2 E1     E1   xy       v    y  c1  E2 E4 E42     y  E1  3c1 E3  c2  E3  3c1 E5   L13  y    c1  E7   E5   2(1   )  y    E1 E1   y   E2 E4  E42     y v    v c1   E   c  E  ,    1   E1     E1   x y          P  x, f   L21  w   2 f 2 w 2 f 2 w 2 f 2 w 2 f    , xy xy x y R x y x c1  c1 E42 E2 E4   w  w  w c E   E     E1  3c1 E3   c2  E3  3c1 E5  ,   2   E1 E1   xy x  2(1   ) x   2x   2x   E22 EE L22 x     E   2c1  E5   2  2(1   ) y   E1 E1    x    E1  3c1 E3   c2  E3  3c1 E5   x , 2(1   )  L23  y    2 E42    c E   1  E1        E22 EE  E   y  2c1  E5    c12  E7    ,  E3  1     E1 E1  E1   xy    E42 E2 E4    w  w  w c E  c  E   E1  3c1 E3  c2  E3  3c1 E5   , 1     E1 E1   y x y  2(1   ) y    E22 EE  E    2x  L32 x    2c1  E5    c12  E7    ,  E3  1  v   E1 E1  E1   xy     2 y  2 y    E22 E2 E4   E42   L33  y     E   c E   c E           y 2(1   ) x  E1 E1  E1         E1  3c1 E3  c2  E3  3c1 E5    y , 2(1   )  L31  w   c1  21 ... and stability of elastically connected multiple beam system using Timoshenko and high order shear deformation theory A high order shear element for nonlinear vibration analysis of composite layered... static, nonlinear dynamic or nonlinear vibration of FGM or plate porous funcationally graded (PFG) For the nonlinear dynamic and vibration of PFGCP has not carried out Therefore, in order to... minimum remnant thickness gauging using high order shear horizontal modes Dynamic analysis of composite sandwich plates with damping modelled using high order shear deformation theory are carried