Nonlinear dynamic response of nano-composite sandwich annular spherical shells

In this research, the nonlinear dynamic response of functionally graded carbon nanotube reinforced composite (FG-CNTRC) sandwich annular spherical shells supported by Pasternak’ foundation is considered by using the analytical approach. Unlike existing works, the structure has three layers: FG-CNTRC layer – homogeneous core – FG-CNTRC layer.

VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 Original Article Nonlinear Dynamic Response of Nano-composite Sandwich Annular Spherical Shells Nguyen Dinh Duc1,2,*, Vu Thi Thuy Anh1, Vu Thi Huong1, Vu Dinh Quang1,2, Pham Dinh Nguyen1 Advanced Materials and Structures Laboratory, VNU Hanoi, University of Engineering and Technology (UET), 144 Xuan Thuy, Cau Giay, Hanoi, Vietnam Infrastructure Engineering Program, VNU, Vietnam-Japan University (VJU), My Dinh 1, Tu Liem, Hanoi, Vietnam Received 10 July 2019 Accepted 05 August 2019 Abstract: In this research, the nonlinear dynamic response of functionally graded carbon nanotube reinforced composite (FG-CNTRC) sandwich annular spherical shells supported by Pasternak’ foundation is considered by using the analytical approach Unlike existing works, the structure has three layers: FG-CNTRC layer – homogeneous core – FG-CNTRC layer Several examples are considered to analyse the behaviour of this sandwich-structured composite The classical shell theory (CST) is used to derive theoretical formulation delineating nonlinear dynamic response of FG-CNTRC sandwich annular spherical shells The numerical results explain the effect of material, geometrical parameters, and elastic foundations on the nonlinear dynamic response of the annular spherical shell Keywords: Nonlinear dynamic response, annular spherical shell, sandwich structures, FG-CNTRC Introduction Carbon nanotubes (CNTs) was discovered in 1991 by Sumio Iijima [1] with many precious advantages such as high electrical and thermal conductivity, very high tensile strength, highly flexible and elastic, high aspect ratio, goof field emission Therefore, CNTs are applied to the advanced materials as reinforcement components of composites in order to enhance structures’ mechanical properties There Corresponding author Email address: ducnd@vnu.edu.vn https//doi.org/ 10.25073/2588-1124/vnumap.4360 N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 are many publications that studied the static and dynamic behavior for CNTs reinforced composite plates and shells Shen reported that the load-bending moment curves can be increased by CNT reinforcement in [2] Farzam and Hassani [3] researched the thermal buckling and mechanical buckling of FG-CNTRC plate using isogeometric analysis Trang and Tung [4] investigated the static behaviour of CNTRC cylindrical panels using the analytical approach Zhong et al [5] analysed the vibration of FG-CNTRC plate that are circular, annular and sector with arbitrary boundary conditions The behaviour of FGCNTRC plate and shell structures with forced vibration was studied by Frikha [6] The nonlinear dynamic response and vibration of FG-CNTRC plates using third-order shear deformation was analysed by Thanh et al [7] The sandwich-structures are more attention to researchers and engineers due to their outstanding properties such as incredible bending rigidity, lightweight, excellent vibration characteristics These sandwich-structures are made up of a moderately thick homogeneous core connected to thin FGM or nano-composite top and bottom layers The sandwich configuration takes advantage of precious mechanical properties of top and bottom layers that are composed of super-stiffer materials Dastjerdi and Aghadavoudi [8] analysed static behaviour of nano-composite sandwich plates which is reinforced by three types of defected carbon nanotubes using first order shear deformation theory Wang and Shen [9] presented the examination of CNT’s volume fraction and thermal loads on the nonlinear free vibration and bending behavior analysis of a sandwich plate placed on elastic foundation with top and bottom layers reinforced CNTs Then, Shen et al [10] developed this problem to analyze the nonlinear vibration of sandwich plates using the HSDT and von Kármán’s assumption Jalali and Heshmati [11] studied the buckling behavior of circular sandwich plates with three layers (tapered cores and nanocomposite top and bottom layers) The annular spherical shells are widely applied in many engineering fields such as civil, mechanical, aerospace engineering Alwar and Narasimhan [12] studied the nonlinear static behaviour of laminated orthotropic annular spherical shells Dumir et al [13] investigated the axisymmetric buckling behavior of laminated shallow annular spherical cap with various boundary conditions The results of fundamental frequencies of the annular spherical shells are presented by Li et al [14] Duc et al [15-16] presented nonlinear static response analysis of FGM annular spherical shells based on CST In the literature, there are many publications that studied the behaviour of many structures which made from functionally graded carbon nano-tube reinforced composite, but there isn’t publication that research nonlinear dynamic response of FG-CNTRCs sandwich annular spherical shell on elastic foundation subjected to thermal and mechanical loads This paper presents the nonlinear dynamic response of the nano-composite sandwich annular spherical shell supported by elastic foundation The advantages of this work are applied the sandwich-structured configuration with three layers (top and bottom layers FG-CNTRC and a homogeneous core) to improve the mechanical properties of the annular spherical shell A form of stress function for CNTs reinforced nano-composite annular spherical shells are presented The governing equations of the CNTRC sandwich annular spherical shell are obtained by using the CST and von Kármán’s geometrical nonlinearity Material properties and formulation 2.1 Sandwich annular spherical shells Consider an annular spherical shell resting on elastic foundations with the radius of curvature R, base radii r1 , r0 and thickness h The shell is subjected to external pressure q uniformly distributed on the outer surface as shown in Fig.1 The shell is defined in the coordinate system ( , ,z) , where  and N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12  are in the meridional and circumferential direction of the shells, respectively and z is perpendicular to the middle surface positive inwards For a thin annular spherical shell, it is convenient to introduce a variable r , referred to as the radius of parallel circle with the base of shell and defined by r  R sin  Moreover, due to the shallowness of the shell, it is approximately assumed that cos  1, Rd  dr Fig Configuration of a sandwich nano-composite annular spherical shell In this present study, the sandwich annular spherical shells composed of three layers, where the top and bottom layers are made of nano-composite reinforced CNTs and a core layer is homogeneous materials Two different models of CNT reinforcements are namely uniformly distributed (UD-UD) and functionally graded ( X-X ) The annular spherical shell is symmetric and the symmetry line is in midplane of the shell The thickness of the homogeneous core and the face sheets are respectively presented as hc , h f , the total thickness of sandwich annular spherical shells is h  hc  2h f (Fig 1) The volume fractions of the two distribution types, which are considered here, are expressed as follows: * VCNT  VCNT   z * VCNT  k VCNT h  UD (1) FG-X wCNT , with subscript CNT , m are indicators which wCNT   CNT / m    CNT / m  wCNT used to represent CNTs and matrix, w and  are mass fraction and density, respectively The *  where VCNT N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 superscript k  k  1, 2,3 represents the top layer  hc h h h  h f  z   c , the core layer  c  z  c , and 2 2 hc h  z  c  hf 2 The elastic modulus ( E ) and shear modulus ( G ) of the nano-composite layers can be achieved as follows [2]: the bottom layer E11  EmVm  1 E11CNT VCNT , 2 E22 3 G12  Vm VCNT  CNT , Em E22  Vm VCNT  Gm G12CNT (2) in which, V is volume fractions The corrective parameters i (i   3) are shown below Table The corrective parameters for CNTs-polymer composites [7] * VCNT 1 2 3 0.12 0.137 1.022 0.715 0.17 0.142 1.626 1.138 0.28 0.141 1.585 1.110 2.2 Basic equations The strain fields of the shell talking von Karman’s geometrical nonlinearity into account [16]:       z  ,   u w  w        r R  r     r   2    v u w  w      2       , r   R R   r        r   v u v w w      r r  r r r       2w   r   r      w  w         r  r r        r    w w      r r  r   (3) The relationship between strain and stress:  k  Qk  k where Qijk are elastic coefficients Q11k  (4) k k E11k v21 E11k E22 k k k , Q  , Q  , Q66  G12k 12 22 k k k k k k  v12 v21  v12 v21  v12 v21 The force and moment fields of the shell can be defined as: N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12  Ni    M i   hc  h  c h f 1   iT   dz  z hc  h  c 1   iC   dz  z hc hf  1  z  iB   dz i  r , , r  hc (5) Specifically, substituting the results of Eq (4) to Eq (5) leads to:  N    A     M    B  B       D     (6) All the detailed coefficients are shown in Appendix The nonlinear equilibrium equations of the annular spherical shell are given as [16]: N r N r N r N     0, r r  r r N N r N r    0, r  r r   M r  M r M r M r   M  M       r r  r    r   r r  R  N r  N   r r r     w w    w N w    N r N r   q  qE   rN r   r r  r   r   r r   (7)  w w  w ; w is the deflection of the shell, k1 ( N / m3 )   r r r r  is the spring layer foundation stiffness and k2 ( N / m) is the shear layer foundation stiffness of Pasternak foundation The stress function method is applied in this paper, so the selected stress function must be satisfied: F  F 2 F 2 F F Nr   , N  , N     r 2 r r r  r r  r  r A transformation has been used for the FGM annular spherical shell in [16], reused in this study to simplify the calculation process such as : where qE  k1w  k2 w , w  w  w( ), F  F0 ( )e2 with r  r0 e ;   ln r r  ;  r0e ;  r0  r r The annular spherical shell’s compatibility equation is shown:  2 r0  r0   0   w    r r r0     r2   r  ,   2 r r r r  r  r r r r    (8) Substituting the results of Eq (6) into Eqs (7) and (8), after a complex series of calculations yields: T1*  T2*  T3*  T4*  T5*  T6*  q  qE  1 (9a) 2 w t (9b) N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 where the Ti*  i    coefficients are shown as:  *  F0 F0 * * * *  F0 * *  F0 A  A  A  A  A  A  A  11 11 12 66 22 11 22       4   F0 A*  A*  A*  F0   F0 A*  A*  2 F0 A*  A* 22 11 22 12 66 12 66           T1*  2 , r0 e     *  w w  * * *   w * * * *    2 B21  B22  B11  B12  B12  B12  B22  B66  r0 e            w  w  w   * * * * * * * *   B21 T2*  4  5B21  B22  B11  B12  r0 e   4 B21  B22  B11 ,     r0 e        4w 3w * * * * * *   B11   B11   B22  B66   B66  3B22 2                 T5*  T6*   r04 e2 r04 e4 r04 e 2           T3* T4*    4 r0 e    w w 2  w  w  w  w       ,                     F0  F0 F C   C  C  C   C11  C12  C22  C21      21  31 21 22        F0  F0  F0  ,   C11  4C12  C22    5C12   C11  C22   C21   C12         F0  F0   C11  4C31  3C22    C  C  C   11 31 22        4w 2w w C  C  C  C   C14  C23  C13  C24      24  14 23 24        w  w  w    5C13   C14  C23   C24    4C13  C14  C23   C13   ,         w  w   3C14  4C32  C23    C  C  C   14 32 23        F0   F0 F0 F0  F0  w  F0  F0   w     F                        2 2  2   F0  F0   w    F0  F0  F   w                      F0  F0 1  F0  F    R r02     w    ,   w w  w  w  w    2   2 2 r r r  r r  r  r0 e  2w 2w          w           N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 The Eqs (9a, 9b) with the new expressions Ti* ,(i 1, ,6) are basic equations to analyze the dynamic response of annular spherical shells resting on elastic foundations 2.3 Analytical solution The shell is assumed to be simply supported along the periphery and subjected to mechanical loads uniformly distributed on the outer surface and the base edges of the shell Depending on the in-plane behavior at the edge of boundary conditions will be considered in case the edges are simply supported and immovable For this case, the boundary conditions are  w w u  0, w  0,   0, N r  0, N r  0, with   (i.e at r  r0 ) (10)   with N is the fictitious compressive load rendering the immovable edges The solution of the annual spherical shell satisfied the boundary condition can be chosen as: (11) w  We sin(1 )sin(n ), r in which 1  m , a  ln , and m, n are buckling mode parameters in the r and  directions, a r0 respectively By using the supper position method, the form of stress function F0 for the nano-composite annular spherical shells is given as: F0  f1 sin  1   f 21  e f 22 sin  n  sin  1   f3 cos  n  sin  1   f cos  n  (12)  f5 cos  1   f 61  f 62 e sin  n  cos  1   f cos  n  cos  1  ,     in which, fi  i    are obtained by solving the geometrical compatibility equation talking stress function and approximate solution into account All the detailed coefficients are shown in Appendix Applying Galerkin method for the nano-composite sandwich annular spherical shell: ln r1 r0 2   e  sin( 1 )sin(n )d  d  0, (13) where  is the Eq (9b) after introducing Eq.(12) and Eq.(11) The differential equation is obtained to analysis the nonlinear dynamic response as follow: d 2W  A1  dt  r0 B4    A  A3 R  A2 r0  B k  B2 k2   W   8B3  q A W   25   1 W     B4 Rr0 B4  B4    Rr0 B4 r0 B4  (14) Numerical results and discussion The sandwich-structured composite shells composed of three layers The core layer is poly (methyl methacrylate) (PMMA) with material properties shown as [7]: Em  GPa    3.52  0.0034T  , vm  0.34 Two face layers are made of composite reinforced CNTs with PMMA is considered as the matrix The material properties of SWCNTs are shown as [17]: CNT E11 TPa   6.3998  4.338417 103T  7.43 106 T  4.458333 109 T , N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 CNT E22 TPa   8.02155  5.420375 103T  9.275 106 T  5.5625 109 T , CNT G12 TPa   1.40755+3.476208 103 T  6.965 106 T + 4.479167 109 T , CNT v12  0.175,  CNT  1400Kg / m3 3.1 Natural frequencies * Table The effect of R / h ration and VCNT on the natural frequencies  1 / s  of nano-composite sandwich annular spherical shells R h 200 250 300 500 * VCNT  0.12 X X 1034.10 827.34 691.48 425.84 UD  UD * VCNT  0.28 X X UD  UD 505.95 410.65 349.04 233.61 1284.80 1025.90 855.51 520.65 724.16 582.37 489.61 310.80 From the table, it’s clear that the natural frequencies are significantly changed with different types * of R / h ratio, VCNT and CNT reinforcements The table indicates that the natural frequencies increases * when the ratio R / h decreases and VCNT increases In addition, the natural frequencies of the X-X model have higher than the UD-UD mode 3.2 Dynamic response analysis The reaction of the sandwich-structured shell with two models of CNT reinforcements (UD-UD, XX) is shown in figure In this figure, the frequencies of the loading are approximately similar to the natural frequencies of the nano-composite sandwich annular spherical shell (shown in table 2) The frequency of loading is  force  750 1/ s  and  force  350 1/ s  for the X-X and UD-UD models, respectively The effect of the core to face sheet thickness ratio hc / h f on the deflection-time curve of the nanocomposite sandwich shell is considered in both X-X and UD-UD models in figure For the X-X type, when the hc / h f ratio is increased, the amplitude of vibration decreases On the contrary, for the UDUD model, the amplitude of fluctuation increases as the ratio increases * The influence of CNTs’ volume fraction VCNT on the dynamic response of nano-composite sandwich annular spherical shell in figure It shows that the fluctuates of the nano-composite sandwich annular * spherical shell increase when the volume fraction VCNT decrease In addition, it is possible to see the deflection amplitude of the X-X model is smaller than the UD-UD model In other words, the X-X model has better load-bearing capacity than UD-UD model Figure presents the influence of R / h ratio on the deflection-time curve of the nano-composite sandwich annular spherical shell reinforced by CNTs It’s clear that the strength of the sandwichstructure shell is better with lower R / h ratio N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 x 10 -4 x 10 -4 1.5 FG X-X UD-UD UD-UD hc=hf hc=2hf 0.5 W(m) W(m) 0 -0.5 -2 -1 -4 X-X (m,n)=(1,1), V*CNT =0.12, R/r0=25, R/r1=2 -1.5 -6 (m,n)=(1,1),V*CNT =0.12, hc=2hf, R/r0=20 R/r1=2, R/h=250, k1=2MPa/m, k2=1MPa.m -8 0.05 0.1 0.15 0.2 0.25 -2 R/h=250, k1=2MPa/m, k2=1MPa.m, q=5000sin(150t) 0.05 0.1 0.15 0.2 0.25 t(s) 0.3 t(s) Fig.2 Effect of CNT models on the nonlinear dynamical response x 10 Fig Effect of core to face sheet thickness ratio hc / h f on the nonlinear dynamical response -4 V*CNT =0.12 UD-UD V*CNT =0.28 x 10 -4 R/h=200 R/h=250 R/h=300 W(m) W(m) 0 -1 -1 X-X -2 (m,n)=(1,1), UD-UD Type, hc=2hf, R/r0=20 (m,n)=(1,1), hc=2hf, R/r0=20, R/r1=2, R/h=250 -2 k1=2MPa/m, k2=1MPa.m, q=5000sin(150t) 0.05 0.1 0.15 0.2 -3 0.25 R/r1=2, k1=2MPa/m, k2=1MPa.m, q=5000sin(100t) 0.05 0.1 Fig Effect of CNTs volume fraction on the deflection-time curve x 10 -5 x 10 R/r1=3.0 R/r1=4.0 0.25 0.3 0.35 -5 R/r0=20 R/r0=25 R/r0=30 W(m) W(m) 0.2 Fig Effect of geometric ratio R / h on the deflection-time curve R/r1=2.0 0.15 t(s) t(s) -2 -1 -4 -2 -6 (m,n)=(1,1), X-X Type, hc=2hf, R/h=250, R/r0=20, k1=2MPa/m, k2=1MPa.m, q=5000sin(150t) -8 0.05 0.1 0.15 0.2 0.25 t(s) -3 -4 (m,n)=(1,1), X-X Type, hc=2hf, R/h=250, R/r1=2, k1=2MPa/m, k2=1MPa.m, q=5000sin(100t) 0.05 0.1 0.15 0.2 0.25 0.3 t(s) Fig Effect of R / r1 ratio on the deflection-time curve Fig Effect of R / r0 ratio on the deflection-time curve 10 N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 Figures and indicate the influence of curvature to radii ratios  R / r1 , R / r0  on the dynamic response of the nano-composite sandwich annular spherical shell Clearly, the fluctuation of the shell is affected by the curvature to radius ratio Conclusions In this study, by using the supper position method, the form of stress function has found to solve the problem of the nanocomposite annular spherical shells The conclusions are obtained: The deflection-time curve and natural frequencies of the nano-composite sandwich shells are considerably changed by geometrical parameters Among the CNT reinforcement models, the X-X model is the best for strength of the sandwich annular spherical shell CNTs’ volume fraction does wonder for the load capacity of the annular spherical shell Acknowledgment This research is funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 107.02-2018.314 The authors are grateful for this support References [1] S Iijima, Helical microtubules of graphitic carbon, Nature 354 (6348) (1991) 56 [2] H.S Shen, Nonlinear bending of functionally graded carbon nanotube-reinforced composite plates in thermal environments, Compos Struct., 91 (1) (2009) 9-19 [3] A Farzam, B Hassani, Thermal and mechanical buckling analysis of FG carbon nanotube reinforced composite plates using modified couple stress theory and isogeometric approach, Compos Struct 206 (2018) 774-790 [4] L.T.N Trang, H Tung, Buckling and postbuckling of carbon nanotube-reinforced composite cylindrical panels subjected to axial compression in thermal environments, Vietnam Journal of Mechanics 40 (1) (2018) 47-61 [5] R Zhong, Q Wang, J Tang, C Shuai, B Qin, Vibration analysis of functionally graded carbon nanotube reinforced composites (FG-CNTRC) circular, annular and sector plates, Compos Struct 194 (2018) 49-67 [6] A Frikha, S Zghal, F Dammak, Dynamic analysis of functionally graded carbon nanotubes-reinforced plate and shell structures using a double directors finite shell element, Aerosp Sci Technol 78 (2018) 438-451 [7] N.V Thanh, N.D Khoa, N.D Tuan, P Tran, N.D Duc, Nonlinear dynamic response and vibration of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) shear deformable plates with temperature-dependent material properties and surrounded on elastic foundations, J Therm Stress 40 (10) (2017) 1254-1274 [8] R Moradi-Dastjerdi, F Aghadavoudi, Static analysis of functionally graded nano-composite sandwich plates reinforced by defected CNT, Compos Struct 200 (2018) 839-848 [9] Z.X Wang, H.S Shen, Nonlinear vibration and bending of sandwich plates with nanotube-reinforced composite face sheets, Compos Part B-Eng 43 (2012) 411- 421 [10] H.S Shen, H Wang, D.Q Yang, Vibration of thermally postbuckled sandwich plates with nanotube-reinforced composite face sheets resting on elastic foundations, Int J Mech SCI 124-125 (2018) 253- 262 [11] S.K Jalali, M Heshmati, Buckling analysis of circular sandwich plates with tapered cores and functionally graded carbon nanotubes-reinforced composite face sheets, Thin-walled struct 100 (2016) 14-24 [12] R.S Alwar, M.C Narasimhan, Axisymmetric non-linear analysis of laminated orthotropic annular spherical shells, Int J Nonlinear Mech 27 (1992) 611-622 N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 11 [13] P.C Dumir, G.P Dube, Mallick A Axisymmetric buckling of laminated thick annular spherical cap Commun Nonlinear Sci Numer Simul 10 (2015) 191-204 [14] H Li, F Pang, H Chen A semi-analytical approach to analyze vibration characteristics of uniform and stepped annular-spherical shells with general boundary conditions Eur J Mech A Solids 74 (2019) 48-65 [15] N.D Duc, D.H Bich, V.T Thuy Anh On the nonlinear stability of eccentrically stiffened functionally graded annular spherical segment shells Thin-Walled Struct., 106 (2016) 258-267 [16] V.T Thuy Anh, D.H Bich, N.D Duc Nonlinear stability analysis of thin FGM annular spherical shells on elastic foundations under external pressure and thermal loads Eur J Mech A Solids 50 (2015) 28-38 [17] M Mirzaei, Y Kiani Thermal buckling of temperature dependent FG-CNT reinforced composite conical shells Aerosp Sci Tech 47 (2015) 42-53 Appendix f1  W a11W ; f3   16 a12 f 61   a42 a43  a41a44  ; f W  a22 a24  a21a84  a22  a212 a42  a412 ; f7    a62W a W W a51 ; f5   ; f 62   25 ; a61 16 a12 a23 W  a41a43  a44 a42  a42  a412 ; f 21   W  a21a24  a84 a22  a22  a212 ; f 22   a85W ; a23 a11  12 A1112 n  A1114  A11n  A22 n  A1112  A22 12  A11  A22  1 ; 16 16 A112  24 A112 14  814 A22 A11  A112 12  A22 A1112  12 A22  a12    ;  A  2A A  A 22 11 22  11  a21   2n A12 1  n A66 1  A1113  1 A11  1 A22  ;  n A12 12  / A1112  1/ A22 12  A22 n  n A12  1/ 2n A66  a22    ;  1/ A n  1/ A   1/ 2n A  22 11 66   a23  A22 n  n A12 12  0.5 A22 12  A1112  A22 n A66 12  A1114  A22 n4  ; a24  1r0 , a25  1 B21  /  1 B11  13 B22  n 1 B22  1 B22  n 1 B11  13 B11  ; a84  /  r0  n r0  12 r0  ;   B12 12  B12  n B11 12  B21 14  a85    B11  B21  n B12  /  ;   B n  n B   B   n B     22 21 66    12  a41   16 n 1 A12  n 1 A66  16 A11 13  1 A11  1 A22   A22 n  16 n A12 12  n A66 12  A11 14  n A12  A22 n  a42     2 n A  10 A   A  66 11 22   a43   1 / 13  / n 1  / 1  , a44   1 / 12    A n  A11 12 n  12 A22 n  A11 n  A22 n  A11 12  a51   112   A   A  A 11 22  22  a61   16 A22 n  n A12  A22 n  n A66  , a62   1 / n 12  N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 12 A1    a62 2 n m  m 2  n a  a  4a61a   a51 2 n m  m 2  n a  a 2 2 64a12 a     ma11 m 2  3n a  a 64a12 a     m 2  a  n a  a41a43  a44 a42   m m 2  a  n a  a41a44  a42 a43      a42 a     m   a na61   2   a412 8 am n  1  e a a62 A2        m 2  a m 2  a   2ana12 m   a a  na 41a44 41 2   m  e a a11   a42 a43   a42       9m   a   na 12  8 m 3m 4  3a m 2 n  a m 2  n a  a  e a    a42 a  ma51 3m 2  a  e a 16 m m 2  n a  a  e a    a412 a  3a m 2  a m 2  a n a412  a42 41a43  a44 a42  ;  ; N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12  m a51  e 3 a A3  12 m C12  a  7C11  C12  5C21  C22    1  2     a m 17C22  17 C11  51C12  3C21    2na  m   a  m   a  a12 4 m  e 2  a   am3 a85  9n  10    a ma85  5n      1    3 a m a25  3n    a a25 1  n   3 m a25    3a23 a  m  a  9 m  a  n   a  3 m n  2 m  a n   a21 a24  a22 a84     8   a m  2n  1  a   m3  a a  a a   21 84 22 24    2 9m  a21  a22  a n    16 amn  e  a  1 n C21  C21  C22  C31   m  a     m  6C11  40 C12  C22    2   2   a m  17 C  31 C  C  17 C     11 12 21 22 2a  m   a  m   a  na12     a  C11  C12  C21  C22    4 12 a m n  C11  C22  C31    4    a m  3C21  17C22  17C11  51C12     6 4 8 m  e  a  1  a41 a43  a44 a42   12 m C12  12 a m n C21  a n C21   2 3na  m  a  9 m  a  a41  a42    a m2 n2  46C31  12C11  6C21  34C22      a m  7C  C  5C  C   11 12 21 22     2a n  C  C  C   21 22 31   4 4  m  6C11  40 C12  C22   16 a n C21   2 2     a m n  22C11  34 C22  56 C31   2 a 8 m  e  1  a41 a44  a42 a43   2     a m 17C11  31C12  C21  17 C22   ; 2 2 2 2 3na  m  a  9 m  a  a41  a42      a n  2C11  8C21  14 C22  16 C31   a4 C  C  C  C   11 12 21 22   m   e  a  1 a11 A4     m2  n2 a  a 4a ; 13 N.D Duc et al / VNU Journal of Science: Mathematics – Physics, Vol 35, No (2019) 1-12 14          a85 m 2  n a  a m 2 m 2  n a  a e a   a21a24  a84 a22    2  4aa23 a m 2  a a21  a22  A5    m m2  n a  a e a   a21 a84  a22 a24     ma25  2  2a23 2a m 2  a a21  a22            ;       m 4C13   a m n  C14  C23  2C32   a n 4C24    m  1  e   2 A6  2   a m  C13  C14  C23  C24    8a   m  a     a n  C14  3C23  2C24  2C32   a  C23  C24  C14     4 4 2 2  2 m C12  2a n C21  2 a m n  C11  C22  2C31    m 2  1  e  a   a21a24  a84 a22   2    a m  7C22  7C11  14C12  2C21   2 2 2a  4 m  a  a21  a22     a n  8C31  C11  4C21  7C22   2a  C11  C12  C21  C22     4 4 2 2   m C12  a n C21   a m n  C11  C22  2C31     m 2 1  e 2 a    a m  C12  C21  C22  C11   a85    a n  C11  2C21  C22   a  C21  C22  C11      8a  m  a  a23 2 a 4 2 2   m 1  e 2 a  a25   m  C11  C12  C22   2 a m n  C22  C31    8a  m  a  a23   a 4C21  n  1   a m  C11  C12  C21  C22   2a n 2C21    m  2C11  9C12  2C22   a m n  3C11  7C22  10C31      m  1  e  a   a21a84  a22 a24    a m  7C22  7C11  9C12  5C21     2a n  C31  C21  C22   a n 4C21  ;  2 2 2a  4 m  a  a21  a22   A7   ; 4 am e a  m  2  a n  ; B   a  n  1  m   m   n a  B  8 a  m   4ar  a  m   2 am  e 1    m   e  1 B  ;B  a m  a  m  n   3m2 e2 a  1 2 2 2 a 2 2a 2 2 2 2 2 2 ... results of fundamental frequencies of the annular spherical shells are presented by Li et al [14] Duc et al [15-16] presented nonlinear static response analysis of FGM annular spherical shells. .. fraction VCNT on the dynamic response of nano-composite sandwich annular spherical shell in figure It shows that the fluctuates of the nano-composite sandwich annular * spherical shell increase when... spherical shell A form of stress function for CNTs reinforced nano-composite annular spherical shells are presented The governing equations of the CNTRC sandwich annular spherical shell are obtained