Static and Dynamic Stability of Carbon Fiber Reinforced Polymer Cylindrical Shell Subject to Non-Normal Boundary Condition with One Generatrix Clamped

Abstract

In this paper, static and dynamic stability analyses taking axial excitation into account are presented for a laminated carbon fiber reinforced polymer (CFRP) cylindrical shell under a non-normal boundary condition. The non-normal boundary condition is put forward to signify that both ends of the cylindrical shell are free and one generatrix of the shell is clamped. The partial differential motion governing the equations of the laminated CFRP cylindrical shell with a non-normal boundary condition is derived using the Hamilton principle, nonlinear von-Karman relationships and first-order deformation shell theory. Then, nonlinear, two-freedom, ordinary differential equations on the radial displacement of the cylindrical shell are obtained utilizing Galerkin method. The Newton-Raphson method is applied to numerically solve the equilibrium point. The stability of the equilibrium point is determined by analyzing the eigenvalue of the Jacobian matrix. The solution of the Mathieu equation describes the dynamic unstable behavior of the CFRP laminated cylindrical shells. The unstable regions are determined using the Bolotin method. The influences of the radial line load, the ratio of radius to thickness, the ratio of length to thickness, the number of layers and the temperature field of the laminated CFRP cylindrical shell on static and dynamic stability are investigated.

1. Introduction

Carbon fiber reinforced polymer (CFRP) laminates are widely used in many engineering fields, such as the ship, vehicle, and aerospace industries, because of their high strength, excellent material performance, light weight, high heat resistance and anti-corrosion properties. In recent years, scholars have carried out research on the mechanical properties of carbon fiber composite materials in order to expand their application range and maintain their reliability [1,2,3]. The stability characteristics of structures made of CFRP ensure their security and reliability. In unstable conditions, the vibration amplitude of the structure is unbounded and increases exponentially with time. Since the resulting vibration may completely destroy the structural members, leading to structural mutations, predictions of structural stability are of the utmost importance from the point of view of both design and optimization [4]. Cylindrical shells are among the most widely used structures in many engineering fields, such as rocket and aircraft propulsion systems and large deployable space annular antenna [5,6]. Hence, it is necessary to understand and predict the nonlinear stability characteristics of CFRP laminated cylindrical shells.

Numerous investigations of the stability characteristics of beam, plate and shell structures have been published to date. Kiral et al. [7] described the dynamic stability of a composite cantilever beam under periodic axial load delamination at predetermined positions. Ke et al. [8] studied the dynamic stability of functionally graded microbeams. In that report, the effects of gradient index, length scale parameters, the slenderness ratio and end supports on static buckling, free vibrations and the dynamic stability of FGM microbeams are discussed in detail. Couto et al. [9] studied the influence of non-uniform bending on transverse torsional buckling of slender steel beams at high temperature. Talebitooti [10] studied the buckling of laminated conical shells made of composite materials under uniformly distributed external loads according to first-order shear deformation theory. Maali et al. [11] studied the buckling behavior of thin defective conical plates under basic supported conditions. Bich et al. [12] studied the linear buckling behavior of functionally graded tapered plates under axial and external pressures. Gajdzicki et al. [13] carried out research on the stability of bi-directionally corrugated plates under compression and shear. Zeng et al. [14] studied the stability and vibrations of rectangular plates with side cracks. Dey et al. [15] studied the dynamic instability and post-buckling behavior of a composite, supported cylindrical shell plate under dynamic local edge load and transverse patch load. Finally, Han et al. [16] studied the dynamic stability of cylindrical shells under periodic axial loads with varying rotational speeds.

However, while there are numerous studies on the dynamic response of carbon fiber composites, few have examined their stability characteristics. Kolanua et al. [17] investigated the stability behavior and failure characteristics of carbon fiber reinforced polymer (CFRP) composite panels with a secondary bonded blade stiffener under compression. The suitability of a CFRP plate subjected to low-velocity impacts for the estimation of the critical load of delamination onset and the approximation of the load-displacement curve are investigated by Salvetti et al. [18]. Cui et al. [19] studied the failure process of CFRP electromagnetic riveting joints under high-speed loading. The deformation and stress capacity of CFRP was studied by Zhang et al. [20]. Juntanalikit et al. [21] studied the cyclic performance of reinforced concrete columns with non-ductile CFRP jackets by experimental and numerical methods. Reuter et al. [22] studied the shear strength of GFRP tubular structures using novel simulation methods. Time et al. [23] studied the fire stability of a CFRP shell structure with a medium-sized test device. Zhang and Zhao [24,25,26] studied the nonlinear response of a laminated CFRP cantilever plate under the action of moment excitation, in-plane airflow and supersonic airflow.

Cylindrical shells are often used as structural units. Hwu et al. [27], Viswanathan et al. [28] and Sarkheil et al. [29] respectively studied the free vibrations of a composite sandwich plate and cylindrical shell, an anti-symmetric cylindrical shell and a cylinder-conical shell. The nonlinear vibrations of water-filled cylindrical shells were studied by Amabili et al. [30]. Song et al. [31] studied the vibration behavior of carbon nanotube-reinforced, composite, closed cylindrical shells using Reddy’s high-order shear deformation theory. Zhang et al. [32] studied the nonlinear dynamics of a clamped, functional gradient material cylindrical shell under complex combined loads. Du et al. [33] discussed the internal resonance behavior of FGM cylindrical shells under a thermal environment. Sun et al. [34] studied the multi-pulse chaotic motion of a circular grid antenna and the nonlinear dynamics of an equivalent cylindrical shell. Liu et al. [35] studied the nonlinear vibrations of composite cylindrical shells with radial prestretched films at the ends. Wang [36] studied the nonlinear vibrations of rotating, composite laminated cylindrical shells with large amplitudes near the lowest resonance under radial harmonic excitation. Hao et al. [37] studied the aerodynamic and thermoelastic flutter characteristics of ceramic-metal gradient truncated conical shells. Wang et al. [38,39] studied the nonlinear dynamic response of rotating cylindrical shells under spectral neighborhood harmonic excitation using numerical methods and approximate analytical solutions. Shen et al. [40,41] studied large amplitude, nonlinear vibrations of shear deformed FGM cylindrical shells surrounded by elastic media.

Non-normal boundary conditions, i.e., when both ends are free and one generatrix of the shell is clamped, often occur in cylindrical shells, e.g., large annular antenna structures. However, few researchers have studied the stability of cylindrical shells under non-positive boundary conditions. In the present research, nonlinear static and dynamic stability analyses of CFRP laminated cylindrical shells with non-normal boundary conditions are carried out. Based on von-Karman-type nonlinear relationships, FSDT and the Hamilton principle, the nonlinear dynamic equation of CFRP laminated cylindrical shells was established using the Galerkin method and expressed as an ordinary differential equation describing radial displacement. The newton-Raphson method is used to numerically analyze the equilibrium point, and local stability is determined by the eigenvalues of the Jacobian matrix. The solution of the Mathieu equation describes the dynamic unstable behavior of a CFRP laminated cylindrical shell. The correctness of the results in this paper is verified by comparisons with the existing results. The influence of radial line load, the ratio of radius to thickness, the ratio of length to thickness, the number of layers and the temperature field on the static and dynamic stability of a CFRP laminated cylindrical shell is studied by parameterization.

2. Equations of Motion

A mechanical model of carbon fiber-reinforced, polymer laminated, cylindrical shells with length $L$, middle surface radius $R$ and uniform thickness $h$, as shown in Figure 1, is considered. There are $N_s$ layers with a ply stacking sequence of (45/−45)s. The curvilinear coordinate system $\left(x,\theta ,z\right)$ is located in the mid-surface of the CFRP laminated cylindrical shell along the axial direction, the circumferential direction and the radial direction, respectively. Displacement components $u$, $v$ and $w$ represent the displacements of an arbitrary point in directions $x$, $\theta$ and $z$, respectively. Non-normal boundary of cylindrical shells which are free at both ends and clamped at $\theta =0$, i.e., one of the longitudinal sections, are considered, as shown in Figure 1a. Figure 1b presents the sections of $x=L$ and $x=0$. The temperatures of the cylindrical shell surface are $T_o$ and $T_{ref}$, respectively. Axial excitation P is loaded at both ends ($x=0$, $x=L$) of the CFRP laminated cylindrical shell.

(1)$P=p_0+p_1\cos \left(\Omega t\right)$

where $p_0$ and $p_1\cos \left({\Omega }_{2}t\right)$ are static and dynamic harmonic excitation, respectively.

According to first-order shear deformation theory [42], it is assumed that the displacement field of CFRP laminated cylindrical shells is

(2)$u\left(x,\theta ,z\right)=u_0\left(x,\theta \right)+z{\phi }_x\left(x,\theta \right)$

(3)$v\left(x,\theta ,z\right)=v_0\left(x,\theta \right)+z{\phi }_{\theta }\left(x,\theta \right)$

(4)$w\left(x,\theta ,z\right)=w_0\left(x,\theta \right)$

where $u_0$, $v_0$ and $w_0$ represent the mid-plane displacements in directions $x$, $\theta$ and $z$, respectively. ${\phi }_x$ and ${\phi }_{\theta }$ denote radial rotations in the $\theta$ and $x$ directions, respectively.

Displacement field Equations (2)–(4) is substituted into the von Karman geometric nonlinear strain-displacement relation [43], and the nonlinear strain is determined as:

(5)$\left\{\begin{array}{l}{\epsilon }_x\\ {\epsilon }_{\theta }\\ {\gamma }_{\theta z}\end{array}\right\}=\left\{\begin{array}{l}{\epsilon }_x^{\left(0\right)}\\ {\epsilon }_{\theta }^{\left(0\right)}\\ {\gamma }_{x\theta }^{\left(0\right)}\end{array}\right\}+z\left\{\begin{array}{l}{\epsilon }_x^{\left(1\right)}\\ {\epsilon }_{\theta }^{\left(1\right)}\\ {\gamma }_{x\theta }^{\left(1\right)}\end{array}\right\}, \left\{\begin{array}{l}{\gamma }_{\theta z}\\ {\gamma }_{xz}\end{array}\right\}=\left\{\begin{array}{c}{\phi }_{\theta }+\frac{1}{R}\frac{\partial w_0}{\partial \theta }-\frac{1}{R}{v}_0\\ \frac{\partial w_0}{\partial x}+{\phi }_x\end{array}\right\}$

where

(6)$\left\{\begin{array}{l}{\epsilon }_x^{\left(0\right)}\\ {\epsilon }_{\theta }^{\left(0\right)}\\ {\gamma }_{x\theta }^{\left(0\right)}\end{array}\right\}=\left\{\begin{array}{c}\frac{\partial u_0}{\partial x}+\frac{1}{2}{\left(\frac{\partial w_0}{\partial x}\right)}^2\\ \frac{1}{R}\frac{\partial v_0}{\partial \theta }+\frac{1}{R}{w}_0+\frac{1}{2R^2}{\left(\frac{\partial w_0}{\partial \theta }\right)}^2\\ \frac{1}{R}\frac{\partial u_0}{\partial \theta }+\frac{\partial v_0}{\partial x}+\frac{1}{R}\frac{\partial w_0}{\partial x}\frac{\partial w_0}{\partial \theta }\end{array}\right\}, \left\{\begin{array}{l}{\epsilon }_x^{\left(1\right)}\\ {\epsilon }_{\theta }^{\left(1\right)}\\ {\gamma }_{x\theta }^{\left(1\right)}\end{array}\right\}=\left\{\begin{array}{c}\frac{\partial {\phi }_x}{\partial x}\\ \frac{1}{R}\frac{\partial {\phi }_{\theta }}{\partial \theta }\\ \frac{1}{R}\frac{\partial {\phi }_x}{\partial \theta }+\frac{\partial {\phi }_{\theta }}{\partial x}\end{array}\right\}$

where ${\epsilon }_x$ and ${\epsilon }_{\theta }$ are the principal strains, and ${\gamma }_{x\theta }$, ${\gamma }_{\theta z}$, and ${\gamma }_{xz}$ denote the shear strains.

The constitutive relationship of laminated CFRP cylindrical shell, considering thermal stress, may be written as

(7)${\left\{\begin{array}{l}{\sigma }_x\\ {\sigma }_{\theta }\\ {\sigma }_{x\theta }\\ {\sigma }_{\theta z}\\ {\sigma }_{xz}\end{array}\right\}}^{\left(k\right)}={\left[\begin{array}{ccccc}{\overline{Q}}_{11}& {\overline{Q}}_{12}& 0& 0& 0\\ {\overline{Q}}_{12}& {\overline{Q}}_{22}& 0& 0& 0\\ 0& 0& {\overline{Q}}_{66}& 0& 0\\ 0& 0& 0& {\overline{Q}}_{44}& 0\\ 0& 0& 0& 0& {\overline{Q}}_{55}\end{array}\right]}^{\left(k\right)}{\left\{\left\{\begin{array}{l}{\epsilon }_x\\ {\epsilon }_{\theta }\\ {\gamma }_{x\theta }\\ {\gamma }_{\theta z}\\ {\gamma }_{xz}\end{array}\right\}-\left\{\begin{array}{c}{\alpha }_x\\ {\alpha }_{\theta }\\ 2{\alpha }_{x\theta }\\ 0\\ 0\end{array}\right\}\Delta T\left(z\right)\right\}}^{\left(k\right)}$

where ${\overline{Q}}_{ij}$ $\left(i,j=1,2,4,5,6\right)$ are the stiffness coefficients, $\Delta T$ is the temperature increment, and ${\alpha }_x$, ${\alpha }_{\theta }$ and ${\alpha }_{x\theta }$ are the coefficients of thermal expansion, which are expressed by

(8)${\alpha }_x={\alpha }_1{\cos }^2\beta +{\alpha }_2{\sin }^2\beta$

(9)${\alpha }_{\theta }={\alpha }_1{\sin }^2\beta +{\alpha }_2{\cos }^2\beta$

(10)${\alpha }_{x\theta }=\left({\alpha }_1-{\alpha }_2\right)\sin \beta \cos \beta$

where ${\alpha }_1$ and ${\alpha }_2$ are the coefficients of the thermal expansion in the different material directions, respectively.

It is supposed that the laminated CFRP cylindrical shell is initially stress free at $T_{ref}$. Assuming that the temperature increment is linear, i.e.,

(11)$\Delta T=T_{ref}+\frac{z}{h}\left(T_0-T_{ref}\right)$

then the stiffness coefficients ${\overline{Q}}_{ij}$ are given by

(12)$\left\{\begin{array}{c}{\overline{Q}}_{11}\\ {\overline{Q}}_{12}\\ {\overline{Q}}_{22}\\ {\overline{Q}}_{16}\\ {\overline{Q}}_{26}\\ {\overline{Q}}_{66}\end{array}\right\}=\left\{\begin{array}{cccc}{C}^4& 2C^2{S}^2& S^4& 4C^2{S}^2\\ C^2{S}^2& C^4+S^4& C^2{S}^2& -4C^2{S}^2\\ S^4& 2C^2{S}^2& C^4& 4C^2{S}^2\\ C^{3}S& C{S}^3-C^{3}S& -C{S}^3& -2CS\left(C^2-S^2\right)\\ C{S}^3& C^{3}S-C{S}^3& -C^{3}S& 2CS\left(C^2-S^2\right)\\ C^2{S}^2& -2C^2{S}^2& C^2{S}^2& {\left(C^2-S^2\right)}^2\end{array}\right\}\left\{\begin{array}{c}{Q}_{11}\\ Q_{12}\\ Q_{22}\\ Q_{66}\end{array}\right\}$

(13)$\left\{\begin{array}{c}{\overline{Q}}_{44}\\ {\overline{Q}}_{45}\\ {\overline{Q}}_{55}\end{array}\right\}=\left\{\begin{array}{cc}{C}^2& S^2\\ -CS& CS\\ S^2& C^2\end{array}\right\}\left\{\begin{array}{c}{Q}_{44}\\ Q_{55}\end{array}\right\}, C=\cos \beta , S=\sin \beta$

where $\beta$ is the ply angle of the laminated CFRP cylindrical shell. The stiffness coefficients of material $Q_{ij}$ are

(14)$$\begin{gathered} Q_{11}=\frac{E_1}{1-{\nu }_{12}{\nu }_{21}}, Q_{12}=\frac{{\nu }_{12}{E}_2}{1-{\nu }_{12}{\nu }_{21}}, Q_{22}=\frac{E_2}{1-{\nu }_{12}{\nu }_{21}}, \\ {Q}_{66}=G_{12}, Q_{44}=G_{23}, Q_{55}=G_{13} \end{gathered}$$

where ${\nu }_{12}$ and ${\nu }_{21}$ are Poisson’s ratios, $E_1$ and $E_2$ are Young’s moduli, and $G_{12}$, $G_{23}$ and $G_{13}$ respectively are the shear modulus of the laminated CFRP cylindrical shell in different material directions.

Based on Hamilton’s principle, a set of nonlinear partial differential governing equations of motion for a CFRP laminated cylindrical shell are obtained, as follows:

(15)$N_{xx,x}+\frac{1}{R}{N}_{x\theta ,\theta }=I_0{\ddot{u}}_0+I_1{\ddot{\phi }}_x$

(16)$N_{x\theta ,x}+\frac{1}{R}{N}_{\theta \theta ,\theta }+\frac{1}{R}{Q}_{\theta }=I_0{\ddot{v}}_0+I_1{\ddot{\phi }}_{\theta }$

(17)$$\begin{gathered} N_{xx,x}\frac{\partial w_0}{\partial x}+N_{xx}\frac{{\partial }^2{w}_0}{\partial x^2}+\frac{1}{R}{N}_{x\theta ,\theta }\frac{\partial w_0}{\partial x}+\frac{2}{R^2}{N}_{x\theta ,\theta }\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+\frac{1}{R}{N}_{xy,x}\frac{\partial w_0}{\partial \theta }-\frac{1}{R}{N}_{\theta \theta }+\frac{1}{R^2}{N}_{\theta \theta ,\theta }\frac{\partial w_0}{\partial \theta } \\ +\frac{1}{R^2}{N}_{\theta \theta }\frac{{\partial }^2{w}_0}{\partial {\theta }^2}+Q_{x,x}+\frac{1}{R}{Q}_{\theta ,\theta }-P\frac{{\partial }^2{w}_0}{\partial x^2}-\gamma {\dot{w}}_0=I_0{\ddot{w}}_0 \end{gathered}$$

(18)$M_{xx,x}+\frac{1}{R}{M}_{x\theta ,\theta }-Q_x=I_1{\ddot{u}}_0+I_2{\ddot{\phi }}_x$

(19)$M_{x\theta ,x}+\frac{1}{R}{M}_{\theta \theta ,\theta }-Q_{\theta }=I_1{\ddot{v}}_0+I_2{\ddot{\phi }}_{\theta }$

where $\gamma$ is the damping coefficient and superscript dots represent the derivative with respect to time. The mass moments of inertia in Equations (15)–(19) are expressed as

(20)$I_{\eta }={\displaystyle \sum _{\eta =1}^N{\displaystyle {\int }_{z_{\eta }}^{z_{\eta +1}}\rho z^i}dz},\hspace{1em}\left(\eta =0,1,2\right)$

The resultant forces of stress and moment are calculated by

(21)$$\begin{gathered} \left\{\begin{array}{l}{N}_{xx}\\ N_{\theta \theta }\\ N_{x\theta }\end{array}\right\}=\left\{\left[A\right],\left[B\right]\right\}\left\{\begin{array}{l}{\epsilon }^{(0)}\\ {\epsilon }^{(1)}\end{array}\right\}-\left\{\begin{array}{c}{N}_{xx}^T\\ N_{\theta \theta }^T\\ N_{x\theta }^T\end{array}\right\}, \left\{\begin{array}{l}{M}_{xx}\\ M_{\theta \theta }\\ M_{x\theta }\end{array}\right\}=\left\{\left[B\right],\left[D\right]\right\}\left\{\begin{array}{l}{\epsilon }^{(0)}\\ {\epsilon }^{(1)}\end{array}\right\}-\left\{\begin{array}{c}{M}_{xx}^T\\ M_{\theta \theta }^T\\ M_{x\theta }^T\end{array}\right\}, \\ \left\{\begin{array}{l}{Q}_x\\ Q_{\theta }\end{array}\right\}=K\left[A\right]\left\{\begin{array}{l}{\gamma }_{xz}\\ {\gamma }_{\theta z}\end{array}\right\} \end{gathered}$$

where $K$ is the shear correction coefficient, given by Efraim as 5/6 [44]. The resulting thermal stress for the CFRP laminated cylindrical shell is defined as

(22)$\left(\left\{\begin{array}{c}{N}_{xx}^T\\ N_{\theta \theta }^T\\ N_{x\theta }^T\end{array}\right\},\left\{\begin{array}{c}{M}_{xx}^T\\ M_{\theta \theta }^T\\ M_{x\theta }^T\end{array}\right\}\right)={\displaystyle \sum _{k=1}^N{{\displaystyle {\int }_{z_k}^{z_{k+1}}\left[\begin{array}{ccc}{Q}_{11}& Q_{12}& 0\\ Q_{12}& Q_{22}& 0\\ 0& 0& Q_{66}\end{array}\right]}}^{\left(k\right)}{\left\{\begin{array}{l}{\alpha }_x\\ {\alpha }_{\theta }\\ {\alpha }_{x\theta }\end{array}\right\}}^{\left(k\right)}\left(\Delta T,\Delta Tz\right)dz}$

The tensile rigidity $A_{ij}$, bending-tensile coupling rigidity $B_{ij}$, and bending rigidity $D_{ij}$ of the laminated CFRP cylindrical shell determined as follows:

(23)$\left(A_{ij},B_{ij},D_{ij}\right)={\displaystyle \sum _{k=1}^N{\displaystyle {\int }_{z_k}^{z_{k+1}}{Q}_{ij}}\left(1,z,z^2\right)dz}, \left(i,j=1,2,6\right)$

(24)$A_{ij}={\displaystyle \sum _{k=1}^N{\displaystyle {\int }_{z_k}^{z_{k+1}}{Q}_{i,j}}\left(1,z,z^2\right)dz}, \left(i,j=4,5\right)$

According to Equations (20)–(24), the nonlinear motion equation can be expressed by the generalized displacement of laminated CFRP cylindrical shells, as follows:

(25)$$\begin{gathered} A_{11}\frac{{\partial }^2{u}_0}{\partial x^2}+A_{66}\frac{1}{R^2}\frac{{\partial }^2{u}_0}{\partial {\theta }^2}+\left(A_{12}+A_{66}\right)\frac{1}{R}\frac{{\partial }^2{v}_0}{\partial x\partial \theta }+B_{11}\frac{{\partial }^2{\phi }_x}{\partial x^2}+B_{66}\frac{1}{R^2}\frac{{\partial }^2{\phi }_x}{\partial {\theta }^2} \\ +\left(B_{12}+B_{66}\right)\frac{1}{R}\frac{{\partial }^2{\phi }_{\theta }}{\partial x\partial \theta }+A_{11}\frac{\partial w_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial x^2}+A_{66}\frac{1}{R^2}\frac{\partial w_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial {\theta }^2} \\ +\left(A_{12}+A_{66}\right)\frac{1}{R^2}\frac{\partial w_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+\frac{A_{12}}{R}\frac{\partial w_0}{\partial x}=I_0{\ddot{u}}_0+I_1{\ddot{\phi }}_x \end{gathered}$$

(26)$$\begin{gathered} A_{66}\frac{{\partial }^2{v}_0}{\partial x^2}+A_{22}\frac{1}{R^2}\frac{{\partial }^2{v}_0}{\partial {\theta }^2}+\left(A_{12}+A_{66}\right)\frac{1}{R}\frac{{\partial }^2{u}_0}{\partial x\partial \theta }+B_{66}\frac{{\partial }^2{\phi }_{\theta }}{\partial x^2}+B_{22}\frac{1}{R^2}\frac{{\partial }^2{\phi }_{\theta }}{\partial {\theta }^2} \\ +\left(B_{12}+B_{66}\right)\frac{1}{R}\frac{{\partial }^2{\phi }_x}{\partial x\partial \theta }+A_{66}\frac{1}{R}\frac{\partial w_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x^2}+A_{22}\frac{1}{R^3}\frac{\partial w_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial {\theta }^2}+\left(A_{12}+A_{66}\right)\frac{1}{R}\frac{\partial w_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial x\partial \theta } \\ +\left(\frac{A_{22}}{R}+\frac{K{A}_{44}}{R}\right)\frac{1}{R}\frac{\partial w_0}{\partial \theta }+\frac{K}{R}{A}_{44}\left({\phi }_{\theta }-\frac{v_0}{R}\right)=I_0{\ddot{v}}_0+I_1{\ddot{\phi }}_{\theta } \end{gathered}$$

(27)$$\begin{gathered} A_{11}\frac{\partial u_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial x^2}+A_{12}\frac{1}{R^2}\frac{\partial u_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial {\theta }^2}+2A_{66}\frac{1}{R^2}\frac{\partial u_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+\left(A_{12}+A_{66}\right)\frac{1}{R^2}\frac{{\partial }^2{u}_0}{\partial x\partial \theta }\frac{\partial w_0}{\partial \theta } \\ +A_{11}\frac{{\partial }^2{u}_0}{\partial x^2}\frac{\partial w_0}{\partial x}+A_{66}\frac{1}{R^2}\frac{{\partial }^2{u}_0}{\partial {\theta }^2}\frac{\partial w_0}{\partial x}-\frac{A_{12}}{R}\frac{\partial u_0}{\partial x}+\left(A_{12}+A_{66}\right)\frac{1}{R}\frac{{\partial }^2{v}_0}{\partial x\partial \theta }\frac{\partial w_0}{\partial x} \\ -\frac{\left(A_{22}+K{A}_{44}\right)}{R}\frac{1}{R}\frac{\partial v_0}{\partial \theta }+2A_{66}\frac{1}{R}\frac{\partial v_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+A_{22}\frac{1}{R^3}\frac{\partial v_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial {\theta }^2}+A_{12}\frac{1}{R}\frac{\partial v_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x^2} \\ +A_{22}\frac{1}{R^3}\frac{{\partial }^2{v}_0}{\partial {\theta }^2}\frac{\partial w_0}{\partial \theta }+A_{66}\frac{1}{R}\frac{{\partial }^2{v}_0}{\partial x^2}\frac{\partial w_0}{\partial \theta }+\left(K{A}_{55}-\frac{B_{21}}{R}\right)\frac{\partial {\phi }_x}{\partial x}+B_{12}\frac{1}{R^2}\frac{\partial {\phi }_x}{\partial x}\frac{{\partial }^2{w}_0}{\partial {\theta }^2} \\ +2B_{66}\frac{1}{R^2}\frac{\partial {\phi }_x}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+\left(B_{12}+B_{66}\right)\frac{1}{R^2}\frac{{\partial }^2{\phi }_x}{\partial x\partial \theta }\frac{\partial w_0}{\partial \theta }+B_{11}\frac{{\partial }^2{\phi }_x}{\partial x^2}\frac{\partial w_0}{\partial x}+B_{11}\frac{\partial {\phi }_x}{\partial x}\frac{{\partial }^2{w}_0}{\partial x^2} \\ +B_{66}\frac{1}{R^2}\frac{{\partial }^2{\phi }_x}{\partial {\theta }^2}\frac{\partial w_0}{\partial x}-\frac{A_{22}{w}_0}{R^2}+\left(K{A}_{44}-\frac{B_{22}}{R}\right)\frac{1}{R}\frac{\partial {\phi }_{\theta }}{\partial \theta }+B_{22}\frac{1}{R^3}\frac{\partial {\phi }_{\theta }}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial {\theta }^2} \\ +2B_{66}\frac{1}{R}\frac{\partial {\phi }_{\theta }}{\partial x}\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+B_{66}\frac{1}{R}\frac{{\partial }^2{\phi }_{\theta }}{\partial x^2}\frac{\partial w_0}{\partial \theta }+B_{12}\frac{1}{R}\frac{\partial {\phi }_{\theta }}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x^2}+\left(B_{12}+B_{66}\right)\frac{1}{R}\frac{{\partial }^2{\phi }_{\theta }}{\partial x\partial \theta }\frac{\partial w_0}{\partial x} \\ +B_{22}\frac{1}{R^3}\frac{{\partial }^2{\phi }_{\theta }}{\partial {\theta }^2}\frac{\partial w_0}{\partial \theta }+2\left(A_{12}+2A_{66}\right)\frac{1}{R}\frac{\partial w_0}{\partial x}\frac{\partial w_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x^2}+\frac{3}{2}{A}_{11}{\left(\frac{\partial w_0}{\partial x}\right)}^2\frac{{\partial }^2{w}_0}{\partial x^2} \\ +\frac{3}{2}{A}_{22}\frac{1}{R^4}{\left(\frac{\partial w_0}{\partial \theta }\right)}^2\frac{{\partial }^2{w}_0}{\partial {\theta }^2}+\left(\frac{A_{12}w}{R}+K{A}_{55}\right)\frac{{\partial }^2{w}_0}{\partial x^2}+\left(\frac{A_{22}w}{R}+K{A}_{44}\right)\frac{1}{R^2}\frac{{\partial }^2{w}_0}{\partial {\theta }^2} \\ +\frac{A_{12}}{2R}{\left(\frac{\partial w_0}{\partial x}\right)}^2+\frac{A_{22}}{2R}\frac{1}{R^2}{\left(\frac{\partial w_0}{\partial \theta }\right)}^2+\left(\frac{1}{2}{A}_{12}+A_{66}\right)\frac{1}{R^2}{\left(\frac{\partial w_0}{\partial x}\right)}^2\frac{{\partial }^2{w}_0}{\partial {\theta }^2} \\ +\left(\frac{1}{2}{A}_{12}+A_{66}\right)\frac{1}{R^2}{\left(\frac{\partial w_0}{\partial \theta }\right)}^2\frac{{\partial }^2{w}_0}{\partial x^2}+N_{xx}^T\frac{{\partial }^2{w}_0}{\partial x^2}+\frac{1}{R^2}{N}_{yy}^T\frac{{\partial }^2{w}_0}{\partial {\theta }^2}+2\frac{1}{R}{N}_{x\theta }^T\frac{{\partial }^2{w}_0}{\partial x\partial \theta }-\frac{N_{yy}^T}{R} \\ -P\frac{{\partial }^2{w}_0}{\partial x^2}-\gamma \frac{\partial w_0}{\partial t}=I_0{\ddot{w}}_0 \end{gathered}$$

(28)$$\begin{gathered} B_{11}\frac{{\partial }^2{u}_0}{\partial x^2}+B_{66}\frac{1}{R^2}\frac{{\partial }^2{u}_0}{\partial {\theta }^2}+\left(B_{12}+B_{66}\right)\frac{1}{R}\frac{{\partial }^2{v}_0}{\partial x\partial \theta }+D_{11}\frac{{\partial }^2{\phi }_x}{\partial x^2}+D_{66}\frac{1}{R^2}\frac{{\partial }^2{\phi }_x}{\partial {\theta }^2}+\left(D_{12}+D_{66}\right)\frac{1}{R}\frac{{\partial }^2{\phi }_{\theta }}{\partial x\partial \theta } \\ -K{A}_{55}{\phi }_x+B_{11}\frac{\partial w_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial x^2}+\left(B_{12}+B_{66}\right)\frac{1}{R^2}\frac{\partial w_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+B_{66}\frac{1}{R^2}\frac{\partial w_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial {\theta }^2} \\ +\left(\frac{B_{12}}{R}-K{A}_{55}\right)\frac{\partial w_0}{\partial x}=I_1{\ddot{u}}_0+I_2{\ddot{\phi }}_x \end{gathered}$$

(29)$$\begin{gathered} B_{66}\frac{{\partial }^2{v}_0}{\partial x^2}+B_{22}\frac{1}{R^2}\frac{{\partial }^2{v}_0}{\partial {\theta }^2}+\left(B_{12}+B_{66}\right)\frac{1}{R}\frac{{\partial }^2{u}_0}{\partial x\partial \theta }+D_{66}\frac{{\partial }^2{\phi }_{\theta }}{\partial x^2}+D_{22}\frac{1}{R^2}\frac{{\partial }^2{\phi }_{\theta }}{\partial {\theta }^2}+\left(D_{12}+D_{66}\right)\frac{1}{R}\frac{{\partial }^2{\phi }_x}{\partial x\partial \theta } \\ -K{A}_{44}\left({\phi }_{\theta }-\frac{v_0}{R}\right)+B_{66}\frac{1}{R}\frac{\partial w_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial x^2}+\left(B_{12}+B_{66}\right)\frac{1}{R}\frac{\partial w_0}{\partial x}\frac{{\partial }^2{w}_0}{\partial x\partial \theta }+B_{22}\frac{1}{R^3}\frac{\partial w_0}{\partial \theta }\frac{{\partial }^2{w}_0}{\partial {\theta }^2} \\ +\frac{B_{22}}{R^2}\frac{\partial w_0}{\partial \theta }-K{A}_{44}\frac{1}{R}\frac{\partial w}{\partial \theta }=I_1{\ddot{v}}_0+I_0{\ddot{\phi }}_{\theta } \end{gathered}$$

The surface of $\theta =0$ is clamped and both ends of the shell are free. This may be expressed by

(30)$u_0=v_0=w_0={\phi }_x={\phi }_y=0 \mathrm{at} \theta =0 \mathrm{and} \theta =2\pi$

(31)$N_{xx}=N_{x\theta }=M_{xx}=M_{x\theta }=Q_x=0 \mathrm{at} x=0 \mathrm{and} x=L$

(32)${\displaystyle {\int }_{-\frac{h}{2}}^{\frac{h}{2}}{N}_{xx}\left|{}_{x=0,L}\right.}Rd\theta ={\displaystyle {\int }_{-\frac{h}{2}}^{\frac{h}{2}}P}Rd\theta$

According to [5,42], displacements u₀, v₀, w₀, ${\phi }_x$ and ${\phi }_{\theta }$ of the shell, which satisfy the non-normal conditions, are written as

(33)$u_0={\displaystyle \sum _{n=1}^M}{\displaystyle \sum _{m=1}^N{u}_{mn}\left(t\right)}\cos \left(\frac{m\pi x}{L}\right)Y_n\left(\theta \right)$

(34)$v_0={\displaystyle \sum _{n=1}^M}{\displaystyle \sum _{m=1}^N{v}_{mn}\left(t\right)}{X}_m\left(x\right)\sin \left(n\theta \right)$

(35)$w_0={\displaystyle \sum _{n=1}^M}{\displaystyle \sum _{m=1}^N{w}_{mn}\left(t\right)}{X}_m\left(x\right)Y_n\left(\theta \right)$

(36)${\phi }_x={\displaystyle \sum _{n=1}^M}{\displaystyle \sum _{m=1}^N{\phi }_x{}_{mn}\left(t\right)}\cos \left(\frac{m\pi x}{L}\right)Y_n\left(\theta \right)$

(37)${\phi }_{\theta }={\displaystyle \sum _{n=1}^M}{\displaystyle \sum _{m=1}^N{\phi }_{\theta }{}_{mn}\left(t\right)}{X}_m\left(x\right)\sin \left(n\theta \right)$

where

(38)$X_i(x)=\sin \frac{{\lambda }_{i}x}{L}+\sinh \frac{{\lambda }_{i}x}{L}-{\alpha }_i(\cosh \frac{{\lambda }_{i}x}{L}+\cos \frac{{\lambda }_{i}x}{L})$

(39)$Y_j(\theta )=\sin \frac{{\mu }_j\theta }{2\pi }-\sinh \frac{{\mu }_j\theta }{2\pi }+{\beta }_j(\cosh \frac{{\mu }_j\theta }{2\pi }-\cos \frac{{\mu }_j\theta }{2\pi })$

(40)$\cos {\lambda }_{i}L\cdot \cosh {\lambda }_{i}L-1=0, \cos {\mu }_{j}2\pi \cdot \cosh {\mu }_{j}2\pi -1=0$

(41)${\alpha }_i=\frac{\sinh {\lambda }_{i}L+\sin {\lambda }_{i}L}{\cosh {\lambda }_{i}L+\cos {\lambda }_{i}L}, {\beta }_j=\frac{\sinh {\mu }_{j}2\pi +\sin {\mu }_{j}2\pi }{\cosh {\mu }_{j}2\pi +\cos {\mu }_{j}2\pi }$

According to Noseir and Bhimaraddi [45,46], the influence of the inertia terms of $u_0$, $v_0$, ${\phi }_x$ and ${\phi }_{\theta }$ in the rotation and in-plane on the nonlinear vibrations of the CFRP laminated cylindrical shell is very small compared to the radial inertia term given in Equation (15). Therefore, inertia terms $u_0$, $v_0$, ${\phi }_x$ and ${\phi }_{\theta }$ can be omitted. Thus, we now focus on the first two modes of transverse displacement $w$. Using Galerkin’s method, both the in-plane and rotational displacement can be expressed as functions of the radial displacement. On this basis, the second order, nonlinear, ordinary differential equation of radial motion of CFRP laminated cylindrical shells is established

(42)$$\begin{gathered} {\ddot{w}}_1+{\mu }_1{\dot{w}}_1+{\omega }_1^2{w}_1+m_2{w}_1^2+m_3{w}_1{w}_2+m_4{w}_2^2+m_5{w}_1^3+m_6{w}_1^2{w}_2^{}+m_7{w}_1{w}_2^2+m_8{w}_2^3 \\ +m_9{w}_1\left(p_1\cos \Omega t\right)=0 \end{gathered}$$

(43)$$\begin{gathered} {\ddot{w}}_2+{\mu }_2{\dot{w}}_2+{\omega }_2^2{w}_2+n_2{w}_1^2+n_3{w}_1{w}_2+n_4{w}_2^2+n_5{w}_1^3+n_6{w}_1^2{w}_2^{}+n_7{w}_1{w}_2^2+n_8{w}_2^3 \\ +n_9{w}_2\left(p_1\cos \Omega t\right)=0 \end{gathered}$$

where ${\omega }_1^2=m_1+m_9{p}_0$ and ${\omega }_2^2=n_1+n_9{p}_0$. All coefficients in Equation (19) can be found in Appendix A.

In order to obtain the dimensionless equation of laminated CFRP cylindrical shells, the following variables and parameters are introduced

(44)$$\begin{gathered} \tau ={\omega }_{1}t, w_1=q_{1}h, w_2=q_{2}h, \overline{\Omega }=\frac{\Omega }{{\omega }_1}, {\overline{\mu }}_1=\frac{{\mu }_1}{{\omega }_1}, {\overline{\mu }}_2=\frac{{\mu }_2}{{\omega }_1}, {\overline{\omega }}_1=\frac{{\omega }_1}{{\omega }_1}, \\ {\overline{\omega }}_2=\frac{{\omega }_2}{{\omega }_1}, {\overline{p}}_0=\frac{p_0}{{\omega }_1^2}, {\overline{p}}_1=\frac{p_1}{{\omega }_1^2}, {\overline{m}}_{\zeta }=\frac{m_{\zeta }h}{{\omega }_1^2}, {\overline{n}}_{\zeta }=\frac{n_{\zeta }h}{{\omega }_1^2}, \left(\zeta =2,3,4\right), \\ {\overline{m}}_{\varsigma }=\frac{m_{\varsigma }{h}^2}{{\omega }_1^2}, {\overline{n}}_{\varsigma }=\frac{n_{\varsigma }{h}^2}{{\omega }_1^2}, \left(\varsigma =5,6,7,8\right) \end{gathered}$$

Equation (19) can be rewritten in non-dimensional form:

(45)$$\begin{gathered} {\ddot{q}}_1+{\overline{\mu }}_1{\dot{q}}_1+{\overline{\omega }}_1^2{q}_1+{\overline{m}}_2{q}_1^2+{\overline{m}}_3{q}_1{q}_2+{\overline{m}}_4{q}_2^2+{\overline{m}}_5{q}_1^3+{\overline{m}}_6{q}_1^2{q}_2^{}+{\overline{m}}_7{q}_1{q}_2^2+{\overline{m}}_8{q}_2^3 \\ +{\overline{m}}_9{\overline{p}}_1{q}_1\cos \overline{\Omega }\tau =0 \end{gathered}$$

(46)$$\begin{gathered} {\ddot{q}}_2+{\overline{\mu }}_2{\dot{q}}_2+{\overline{\omega }}_2^2{q}_2+{\overline{n}}_2{q}_1^2+{\overline{n}}_3{q}_1{q}_2+{\overline{n}}_4{q}_2^2+{\overline{n}}_5{q}_1^3+{\overline{n}}_6{q}_1^2{q}_2^{}+{\overline{n}}_7{q}_1{q}_2^2+{\overline{n}}_8{q}_2^3 \\ +{\overline{n}}_9{\overline{p}}_1{q}_2\cos \overline{\Omega }\tau =0 \end{gathered}$$

where “$\dot{}$” signifies the derivative with respect to dimensionless time “$\tau$”.

3. Static Bifurcation and Stability

In this section, dynamic harmonic excitation is set to zero and static excitation is selected as the controlling parameter to analyze the bifurcation and stability of the CFRP laminated cylindrical shell. The Newton-Raphson method is applied to numerically analyze the equilibrium points. Then, by solving the eigenvalues of the Jacobian matrix, the stability of the equilibrium point is obtained.

Equations (45) and (46) can be rewritten as a first-order system as follows:

(47)${\dot{q}}_1=q_{01}$

(48)${\dot{q}}_{01}=-{\overline{\mu }}_1{q}_{01}-m_1{q}_1-m_9{p}_0{q}_1-{\overline{m}}_2{q}_1^2-{\overline{m}}_3{q}_1{q}_2-{\overline{m}}_4{q}_2^2-{\overline{m}}_5{q}_1^3-{\overline{m}}_6{q}_1^2{q}_2^{}-{\overline{m}}_7{q}_1{q}_2^2-{\overline{m}}_8{q}_2^3$

(49)${\dot{q}}_2=q_{02}$

(50)$$\begin{gathered} {\dot{q}}_{02}=-{\overline{\mu }}_2{q}_{02}-n_1{q}_2-n_9{p}_0{q}_2-{\overline{n}}_2{q}_1^2-{\overline{n}}_3{q}_1{q}_2-{\overline{n}}_4{q}_2^2-{\overline{n}}_5{q}_1^3-{\overline{n}}_6{q}_1^2{q}_2^{}-{\overline{n}}_7{q}_1{q}_2^2-{\overline{n}}_8{q}_2^3 \\ -{\overline{n}}_9{\overline{p}}_1{q}_2\cos \overline{\Omega }\tau \end{gathered}$$

Setting the left parts of Equations (47)–(50) to zero, the nonlinear algebraic equations are expressed as

(51)$q_{01}=0$

(52)$-{\overline{\mu }}_1{q}_{01}-{\overline{m}}_1{q}_1-{\overline{m}}_9{p}_0{q}_1-{\overline{m}}_2{q}_1^2-{\overline{m}}_3{q}_1{q}_2-{\overline{m}}_4{q}_2^2-{\overline{m}}_5{q}_1^3-{\overline{m}}_6{q}_1^2{q}_2^{}-{\overline{m}}_7{q}_1{q}_2^2-{\overline{m}}_8{q}_2^3=0$

(53)$q_{02}=0$

(54)$-{\overline{\mu }}_2{q}_{02}-{\overline{n}}_1{q}_2-{\overline{n}}_9{p}_0{q}_2-{\overline{n}}_2{q}_1^2-{\overline{n}}_3{q}_1{q}_2-{\overline{n}}_4{q}_2^2-{\overline{n}}_5{q}_1^3-{\overline{n}}_6{q}_1^2{q}_2^{}-{\overline{n}}_7{q}_1{q}_2^2-{\overline{n}}_8{q}_2^3=0$

The Jacobian matrix is indicated as

(55)$J=\left[\begin{array}{cccc}0& 1& 0& 0\\ \frac{\partial {\dot{q}}_{01}}{\partial q_1}& \frac{\partial {\dot{q}}_{01}}{\partial q_{01}}& \frac{\partial {\dot{q}}_{01}}{\partial q_2}& \frac{\partial {\dot{q}}_{01}}{\partial q_{02}}\\ 0& 0& 0& 1\\ \frac{\partial {\dot{q}}_{02}}{\partial q_1}& \frac{\partial {\dot{q}}_{02}}{\partial q_{01}}& \frac{\partial {\dot{q}}_{02}}{\partial q_2}& \frac{\partial {\dot{q}}_{02}}{\partial q_{02}}\end{array}\right]$

where

(56)$\frac{\partial {\dot{q}}_{01}}{\partial q_1}=-{\overline{m}}_1-{\overline{m}}_9{p}_0-{\overline{m}}_2{q}_1^{}-{\overline{m}}_3{q}_2-{\overline{m}}_5{q}_1^2-{\overline{m}}_6{q}_1^{}{q}_2^{}-{\overline{m}}_7{q}_2^2$

(57)$\frac{\partial {\dot{q}}_{01}}{\partial q_{01}}=-{\overline{\mu }}_1$

(58)$\frac{\partial {\dot{q}}_{01}}{\partial q_2}=-{\overline{m}}_3{q}_1-{\overline{m}}_4{q}_2^{}-{\overline{m}}_6{q}_1^2-{\overline{m}}_7{q}_1{q}_2^{}-{\overline{m}}_8{q}_2^2$

(59)$\frac{\partial {\dot{q}}_{01}}{\partial q_{02}}=0$

(60)$\frac{\partial {\dot{q}}_{02}}{\partial q_1}=-{\overline{n}}_2{q}_1^{}-{\overline{n}}_3{q}_2-{\overline{n}}_5{q}_1^2-{\overline{n}}_6{q}_1^{}{q}_2^{}-{\overline{n}}_7{q}_2^2$

(61)$\frac{\partial {\dot{q}}_{02}}{\partial q_{01}}=-{\overline{\mu }}_2$

(62)$\frac{\partial {\dot{q}}_{02}}{\partial q_2}=-{\overline{n}}_1-{\overline{n}}_9{p}_0-{\overline{n}}_3{q}_1-{\overline{n}}_4{q}_2^{}-{\overline{n}}_6{q}_1^2-{\overline{n}}_7{q}_1{q}_2^{}-{\overline{n}}_8{q}_2^2$

(63)$\frac{\partial {\dot{q}}_{02}}{\partial q_{02}}=0$

By calculating the equilibrium points of Equations (51)–(54), critical static in-plane load $p_{cr}$, which has nonzero equilibrium points, is found. The stability of the equilibrium point is determined by examining the maximum real part of the eigenvalue of the Jacobian matrix expressed in Equation (55).

In following analysis, a N_s-layer, antisymmetric angle-ply shell (45/−45)s with length $L=1 \mathrm{m}$ is considered, and the shell’s material properties $E_1=140\times {10}^3 \mathrm{MPa}$, $E_2=10\times {10}^3 \mathrm{MPa}$, $G_{12}=7\times {10}^3 \mathrm{MPa}$, $G_{13}=7\times {10}^3 \mathrm{MPa}$, $G_{23}=7\times {10}^3 \mathrm{MPa}$, ${\nu }_{12}=0.25$, ${\alpha }_1=-0.3\times {10}^{-6} \mathrm{m}/\mathrm{K}$ and ${\alpha }_2=28\times {10}^{-6} \mathrm{m}/\mathrm{K}$ are utilized. The temperatures of inner surfaces of the shell is 300 K.

In order to validate the present results, the dimensionless natural frequencies (${\Omega }_n={\omega }_{n}R\sqrt{(1-{\nu }_{}^2)\rho /E}$) are compared with the results of Zhang et al. [47] and Song et al. [48] in Table 1, taking into account simply supported isotropic cylindrical shells with $L/R=20$, $m=1$ and $\nu =0.3$. As shown, the calculated results are in good agreement with existing ones. In addition, the dimensionless axial static buckling load $P_{cr}{L}^2/\left(E_{02}{h}^3\right)$ of the simply supported orthotropic cylindrical shell is calculated and compared with the results of Lee et al. [49] and Gao et al. [50] in Table 2. The geometric parameters and material properties are as follows: $h=0.00254 \mathrm{m}$, $R/h=100$, $L/R=2$, $E_{01}=275.8\times {10}^9 \mathrm{Pa}$, $E_{02}=27.58\times {10}^9 \mathrm{Pa}$, $G_0=10.34\times {10}^9 \mathrm{Pa}$, ${\nu }_{12}=0.25$, ${\nu }_{21}=0.025$ and $\rho =1619.27 {\mathrm{kg}/\mathrm{m}}^3$. As can be seen, the results compare reasonably well.

Figure 2 shows the effects of the ratio of radius to thickness and temperature field on the critical static in-plane load. It is observed that with an increase of the ratio of radius to thickness, the critical static in-plane load decreases monotonically. On the other hand, with an increase of temperature difference between the inner and outer surface, the critical static in-plane load decreases. This is because increasing the ratio of radius to thickness and the temperature field can lead to a decrease in the stiffness of the system. The curves of critical in-plane load versus the ratio of length to thickness $L/h$ with different temperature fields are shown in Figure 3. One can find that the critical static in-plane load decreases while the ratio of length to thickness or the temperature field increases. Figure 4 shows that the curves for a critical in-plane load versus the number of layers $N_s$ when the outer surface temperature $T_o$ is set at 400, 500 and 600, respectively. As in Figure 2 and Figure 3, with the increase of temperature field, the critical static in-plane load decreases. In addition, we observe that as the number of layers increases, the stiffness of the system increases monotonously, as does the critical in-plane load.

Now, the nonlinear static bifurcations and the stabilities of equilibrium points will be investigated. In this regard, point $\left(L/2,\pi ,0\right)$ on the CFRP laminated cylindrical shell is the referential location. Figure 5 illustrates the solution curves of the transverse displacement of the CFRP laminated cylindrical shell with different temperature fields when $N_s=8$, $L/h=120$ and $R/h=30$. As noted in Figure 5a–c, the outer surface temperatures $T_o$ are set at 400, 500 and 600, respectively. Here, the solid line is the stable equilibrium solution and the dashed line is the unstable equilibrium solution. Three solutions, i.e., two stable nonzero solutions and one unstable zero solution, occur when the in-plane load is greater than the critical load. By contrasting Figure 5a–c, we see that with an increase in the outer surface temperature, the critical static load increases, as does the nonzero equilibrium displacement of the referential location. The solution curves for the transverse displacement of the CFRP laminated cylindrical shell versus the static in-plane load with different numbers of layers $N_s$ when $L/h=120$, $R/h=30$ and $T_o=400 \mathrm{K}$ are shown in Figure 6. As shown, the static bifurcation point are 31.6, 28.4 and 24.6 when the number of layers is set as $N_s=8$, $N_s=6$ and $N_s=4$, respectively. With an increase in the numbers of layers $N_s$, the critical static load increases. Figure 7 illustrates the effects of the static in-plane load and the ratio of radius to thickness on the nonlinear static bifurcations and the stabilities of the equilibrium points of the CFRP laminated cylindrical shell when $N_s=8$, $R/h=20$ and $T_o=400 \mathrm{K}$. As shown in Figure 7a–c, the values of the ratios of length to thickness are set at 80, 100 and 120, respectively. As the static in-plane load increases, nonlinear static bifurcation occurs in the system. Additionally, static bifurcation occurs earlier when the ratio of length to thickness is bigger. Figure 8 shows the solution curves for the transverse deflection of the CFRP laminated cylindrical shell versus the static in-plane load with different ratios of radius to thickness $R/h$ when $N_s=6$, $L/h=80$ and $T_o=400 \mathrm{K}$. Increasing the ratio of radius to thickness $R/h$ may cause nonlinear static bifurcation to occur sooner.

4. Dynamic Stability Analysis

In this section, the dynamic stabilities of the CFRP laminated cylindrical shell are investigated. As noted in Equations (45) and (46), dynamic harmonic excitation is selected as the controlling parameter to investigate the dynamic stability of the system. Based on the Liapunov principle and studies [51,52], the dynamic unstable region of the nonlinear dynamic system can be determined by its linear parts.

The Mathieu equations, obtained by omitting all the nonlinear terms in Equations (45) and (46), can be written in the following form:

(64)${\ddot{q}}_1+{\overline{m}}_1{q}_1+{\overline{m}}_9{\alpha }_0{p}_{cr}{q}_1+{\overline{m}}_9{\alpha }_1{p}_{cr}{q}_1\cos \overline{\Omega }\tau =0$

(65)${\ddot{q}}_2+{\overline{n}}_1{q}_2+{\overline{n}}_9{\alpha }_0{p}_{cr}{q}_2+{\overline{n}}_9{\alpha }_1{p}_{cr}{q}_2\cos \overline{\Omega }\tau =0$

where ${\alpha }_0$ and ${\alpha }_1$ are the static and dynamic in-plane load factors, respectively. The static and dynamic loads can be expressed as $p_0={\alpha }_0{p}_{cr}$ and $p_1={\alpha }_1{p}_{cr}$, respectively.

Using the Bolotin method, the approximated solutions with period $T=2\pi /\overline{\Omega }$ are assumed to be

(66)$q_1=a_1\sin \left(\frac{\overline{\Omega }\tau }{2}\right)+b_1\cos \left(\frac{\overline{\Omega }\tau }{2}\right)$

(67)$q_2=a_2\sin \left(\frac{\overline{\Omega }\tau }{2}\right)+b_2\cos \left(\frac{\overline{\Omega }\tau }{2}\right)$

Substituting Equation (27) into Equation (26), and combining the coefficients of the sine and cosine function, we obtain the following equations:

(68)$$\begin{gathered} \left(-\frac{1}{4}{\overline{\Omega }}^2+{\overline{m}}_1+{\overline{m}}_9{\alpha }_0{p}_{cr}-\frac{1}{2}{\overline{m}}_9{\alpha }_1{p}_{cr}\right)a_1\sin \left(\frac{\overline{\Omega }\tau }{2}\right)+\frac{1}{2}{\overline{m}}_9{\alpha }_1{p}_{cr}{a}_1\sin \left(\frac{3\overline{\Omega }\tau }{2}\right) \\ +\left(-\frac{1}{4}{\overline{\Omega }}^2+{\overline{m}}_1+{\overline{m}}_9{\alpha }_0{p}_{cr}+\frac{1}{2}{\overline{m}}_9{\alpha }_1{p}_{cr}\right)b_1\cos \left(\frac{\overline{\Omega }\tau }{2}\right)+\frac{1}{2}{\overline{m}}_9{\alpha }_1{p}_{cr}{b}_1\cos \left(\frac{3\overline{\Omega }\tau }{2}\right)=0 \end{gathered}$$

(69)$$\begin{gathered} \left(-\frac{1}{4}{\overline{\Omega }}^2+{\overline{n}}_1+{\overline{n}}_9{\alpha }_0{p}_{cr}-\frac{1}{2}{\overline{n}}_9{\alpha }_1{p}_{cr}\right)a_2\sin \left(\frac{\overline{\Omega }\tau }{2}\right)+\frac{1}{2}{\overline{n}}_9{\alpha }_1{p}_{cr}{a}_2\sin \left(\frac{3\overline{\Omega }\tau }{2}\right) \\ +\left(-\frac{1}{4}{\overline{\Omega }}^2+{\overline{n}}_1+{\overline{n}}_9{\alpha }_0{p}_{cr}+\frac{1}{2}{\overline{n}}_9{\alpha }_1{p}_{cr}\right)b_2\cos \left(\frac{\overline{\Omega }\tau }{2}\right)+\frac{1}{2}{\overline{n}}_9{\alpha }_1{p}_{cr}{b}_2\cos \left(\frac{3\overline{\Omega }\tau }{2}\right)=0 \end{gathered}$$

Setting the coefficients of $\sin \left(\frac{\overline{\Omega }\tau }{2}\right)$ and $\cos \left(\frac{\overline{\Omega }\tau }{2}\right)$ of Equation (28) to zero, we obtain a series of algebraic equations which can be written as

(70)$-\frac{1}{4}{\overline{\Omega }}^2+{\overline{m}}_1+{\overline{m}}_9{\alpha }_0{p}_{cr}-\frac{1}{2}{\overline{m}}_9{\alpha }_1{p}_{cr}=0$

(71)$-\frac{1}{4}{\overline{\Omega }}^2+{\overline{m}}_1+{\overline{m}}_9{\alpha }_0{p}_{cr}+\frac{1}{2}{\overline{m}}_9{\alpha }_1{p}_{cr}=0$

(72)$-\frac{1}{4}{\overline{\Omega }}^2+{\overline{n}}_1+{\overline{n}}_9{\alpha }_0{p}_{cr}-\frac{1}{2}{\overline{n}}_9{\alpha }_1{p}_{cr}=0$

(73)$-\frac{1}{4}{\overline{\Omega }}^2+{\overline{n}}_1+{\overline{n}}_9{\alpha }_0{p}_{cr}+\frac{1}{2}{\overline{n}}_9{\alpha }_1{p}_{cr}=0$

Based on Equation (29), the dynamic stability of the CFRP laminated cylindrical shell subjected to axial excitation may be analyzed numerically. The unstable regions are plotted by the dynamic load factor against the excitation frequency on the plane $\left({\alpha }_1,\overline{\Omega }\right)$ for the first two modes.

The present dynamic unstable regions of the laminated composite cylindrical shell ($L/R=1$, $R/h=100$) are compared with those of Ganapathi et al. [53] and Dey et al. [54] in Figure 9. In this regard, a laminated composite cylindrical shell with the following material properties is considered: $E_{11}/E_{22}=25$, $G_{23}=0.2E_{22}$, $G_{12}=G_{13}=0.5E_{22}$ and ${\nu }_{12}=0.25$. The present results agree well with the results reported by Ganapathi et al. [53] and Dey et al. [54].

Figure 10 illustrates the effect of the temperature fields on the dynamic unstable regions of the CFRP laminated cylindrical shell when $N_s=8$, $L/h=120$, $R/h=30$ and ${\alpha }_0=0$. Figure 10 show the dynamic unstable regions of the first and second modes, respectively. The part between the two lines is the dynamic unstable region. Inside the dynamic unstable region, the CFRP laminated cylindrical shell vibrates with unbounded amplitudes, and as such, unstable behavior occurs. Outside the dynamic unstable region, the amplitudes of the CFRP laminated cylindrical shell are bounded, i.e., the shell is stable. One can observe that all the dynamic unstable regions become wider with an increase of dynamic in-plane load factor ${\alpha }_1$. Furthermore, with an increase of the temperature field, the unstable regions of both modes are translated to the lower parametric excitation frequency.

Figure 11 shows the linear response for the CFRP laminated cylindrical shell with an excitation frequency in the unstable region. It may be observed that the linear response of the CFRP laminated cylindrical shell grows exponentially and the shell becomes unstable. Figure 12 shows the nonlinear response of the CFRP laminated cylindrical shell with the same parameters. The amplitude of the nonlinear response is the same as that of the assumed initial amplitude, so we may consider this shell to also be unstable. Figure 13 and Figure 14 show the linear and nonlinear responses for the CFRP laminated cylindrical shell with an excitation frequency in the stability region. The frequency of the nonlinear response is higher than that of the linear response.

The dynamic unstable regions are shown in Figure 15; these illustrate the effect of the ratio of length to thickness $L/h$ on the CFRP laminated cylindrical shell with $T_o=400 \mathrm{K}$, $N_s=8$, $R/h=30$ and ${\alpha }_0=0$. It is found that with an increase of $L/h$, the dynamic unstable regions in both modes shift downward. The dynamic unstable regions of the CFRP laminated cylindrical shell with $R/h=10$, $R/h=20$ and $R/h=30$ are depicted in Figure 16 when $T_o=400 \mathrm{K}$, $N_s=8$, $L/h=100$ and ${\alpha }_0=0$. As shown, the value of the unstable frequency decreases with an increase in the ratio of radius to thickness. Figure 17 shows the effect of number of layers $N_s$ on the dynamic unstable regions of the CFRP laminated cylindrical shell with $T_o=400 \mathrm{K}$, $L/h=100$, $R/h=30$ and ${\alpha }_0=0$. The number of layers was set at 8, 6, and 4, respectively. It is well known that when $N_s$ decreases, the unstable regions of the system begin at a lower frequency.

5. Conclusions

This paper presents static and dynamic stability analyses of a carbon fiber reinforced polymer (CFRP) laminated cylindrical shell under axial excitation. Non-normal boundary conditions were applied, i.e., both ends of the cylindrical shell were free and one generatrix of the shell was clamped. Based on von-Karman-type nonlinear relationships, first-order shear deformation theory and the Hamilton principle, the partial differential motion control equation of CFRP laminated cylindrical shells was derived. Using the Galerkin method, the nonlinear ordinary differential motion equation of the shell along the radial displacement was obtained. The newton-Raphson method was used to numerically analyze the equilibrium point, and the local stability was obtained by the eigenvalues of the Jacobian matrix. The Mathieu equation describes the dynamic unstable behavior of the CFRP laminated cylindrical shell. The correctness of the results in this paper was verified by comparisons with existing results. A parametric study was conducted to investigate the effects of the radial line load, the ratio of radius to thickness, the ratio of length to thickness, the number of layers and the temperature field on the static and dynamic stability of a CFRP laminated cylindrical shell. It can be concluded that:

• (1). Bifurcation phenomena might occur when the static in-plane load is greater than the critical load.

• (2). The ratio of radius to thickness, the ratio of length to thickness, the number of layers and the temperature field have significant effects on static bifurcation characteristics of a CFRP laminated cylindrical shell.

• (3). With an increase of ratio of radius to thickness, the ratio of length to thickness and the temperature field, the unstable regions in both modes are translated to a lower parametric excitation frequency.

• (4). When the number of layers decreases, the unstable regions of the system begin at a lower frequency.

Author Contributions

Conceptualization, S.Y. and W.Z.; methodology, Y.H.; software, L.L.; validation, S.Y. and W.M.; investigation, Y.H.; data curation, S.Y.; writing—original draft preparation, S.Y.; writing—review and editing, W.Z.; funding acquisition, S.Y., Y.H., W.Z. and W.M. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge the supports of National Natural Science Foundation of China (NNSFC) through grant Nos. 12002057, 11872127, 11832002 and 12102207, Scientific Research Project of Beijing Educational Committee No. KM202111232023, Qin Xin Talents Cultivation Program, Beijing Information Science & Technology University QXTCP C202102.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Model of a CFRP laminated cylindrical shell: (a) the mechanical model, (b) the sections of [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.].

Figure 2. Critical in-plane load versus the ratio of radius to thickness [Forumla omitted. See PDF.] with different temperature fields.

Figure 3. Critical in-plane load versus the ratio of length to thickness [Forumla omitted. See PDF.] with different temperature fields.

Figure 4. Critical in-plane load versus the number of layers [Forumla omitted. See PDF.] with different temperature fields.

Figure 5. The transverse deflection of the CFRP laminated cylindrical shell versus the static in-plane load with different temperature fields: (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.].

Figure 6. The transverse deflection of the CFRP laminated cylindrical shell versus the static in-plane load with different numbers of layers [Forumla omitted. See PDF.]: (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.].

Figure 7. The transverse deflection of the CFRP laminated cylindrical shell versus the static in-plane load with different ratios of length to thickness [Forumla omitted. See PDF.]: (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.].

Figure 8. The transverse deflection of the CFRP laminated cylindrical shell versus the static in-plane load with different ratios of radius to thickness [Forumla omitted. See PDF.]: (a) [Forumla omitted. See PDF.], (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.].

Figure 9. The dynamic unstable regions of a cross-ply laminated composite cylindrical shell ([Forumla omitted. See PDF.], [Forumla omitted. See PDF.]) subjected to uniform periodic in-plane loading.

Figure 10. The dynamic unstable regions of the CFRP laminated cylindrical shell with different temperature fields: (a) first mode, and (b) second mode.

Figure 11. The linear response for the CFRP laminated cylindrical shell with an excitation frequency in the unstable region: (a) the time history on the plane [Forumla omitted. See PDF.], and (b) the phase portrait on the plane [Forumla omitted. See PDF.].

Figure 12. The nonlinear response for the CFRP laminated cylindrical shell with an excitation frequency in the unstable region: (a) the time history on the plane [Forumla omitted. See PDF.], and (b) the phase portrait on the plane [Forumla omitted. See PDF.].

Figure 13. The linear response for the CFRP laminated cylindrical shell with an excitation frequency in the stable region: (a) the time history on the plane [Forumla omitted. See PDF.], and (b) the phase portrait on the plane [Forumla omitted. See PDF.].

Figure 14. The nonlinear response for the CFRP laminated cylindrical shell with an excitation frequency in the stable region: (a) the time history on the plane [Forumla omitted. See PDF.], and (b) the phase portrait on the plane [Forumla omitted. See PDF.].

Figure 15. The dynamic unstable regions of the CFRP laminated cylindrical shell with different [Forumla omitted. See PDF.]: (a) the first mode, and (b) the second mode.

Figure 16. The dynamic unstable regions of the CFRP laminated cylindrical shell with different [Forumla omitted. See PDF.]: (a) the first mode, and (b) the second mode.

Figure 17. The dynamic unstable regions of the CFRP laminated cylindrical shell with different [Forumla omitted. See PDF.]: (a) the first mode, and (b) the second mode.

Table 1

Comparison of the frequency parameters (${\Omega }_n={\omega }_{n}R\sqrt{(1-{\nu }_{}^2)\rho /E}$) for a simply supported isotropic cylindrical shell with $L/R=20$, $m=1$ and $\nu =0.3$.

$\mathit{h}/\mathit{R}$	$\mathit{n}$	Zhang et al. [47]	Song et al. [48]	Present Study
0.05	0	0.0929586	0.0929392	0.0929465
	1	0.0161065	0.0161299	0.0151185
	2	0.0393038	0.0393231	0.0393236
	3	0.1098113	0.1097653	0.1096523
	4	0.1098113	0.1097653	0.1098103
0.002	0	0.0929296	0.0929296	0.0929236
	1	0.0161011	0.0161011	0.01610032
	2	0.0054532	0.0054536	0.0054532
	3	0.0050418	0.0050424	0.0050423
	4	0.0085340	0.0085344	0.0085354

Table 2

Comparison of a dimensionless axial static buckling load $P_{cr}{L}^2/\left(E_{02}{h}^3\right)$ on a simply supported orthotropic cylindrical shell.

$(\mathit{m},\mathit{n})$	Lee et al. [49]	Gao et al. [50]	Present Study
(1, 1)	78,139.72	78,145.73	78,151.26
(1, 2)	29,556.79	29,580.83	29,578.52
(1, 3)	13,850.67	13,904.75	13,895.28
(2, 1)	32,341.27	32,347.28	32,343.52
(2, 2)	19,852.90	19,876.93	19,856.36
(2, 3)	12,046.73	12,100.82	12,089.65

Appendix A

The coefficients obtained in Equations (42) and (43) are presented as follows (A1)$$\begin{gathered} m_1=\frac{e_5}{e_4}, m_9=\frac{e_2}{e_4}, m_{10}=\frac{e_1}{e_4}, {\mu }_1=\frac{e_3}{e_4}, n_1=\frac{g_5}{g_4}, n_9=\frac{g_2}{g_4}, n_{10}=\frac{g_1}{g_4}, {\mu }_2=\frac{g_3}{g_4}, \\ {m}_2=-\frac{a_1{b}_3{e}_{10}}{e_4\Gamma }+\frac{a_3{b}_1{e}_{10}}{e_4\Gamma }+\frac{a_2{b}_3{e}_9}{e_4\Gamma }-\frac{a_3{b}_2{e}_{10}}{e_4\Gamma }+\frac{a_1{b}_2{d}_3{e}_{10}}{d_2{e}_4\Gamma }+\frac{a_1{b}_3{d}_1{e}_{12}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{d}_3{e}_{12}}{d_2{e}_4\Gamma }-\frac{a_3{b}_1{d}_1{e}_{12}}{d_2{e}_4\Gamma } \\ -\frac{a_1{b}_2{c}_3{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_3{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_2{b}_3{c}_1{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_3{b}_2{c}_1{e}_{11}}{c_2{e}_4\Gamma }+\frac{e_6}{e_4}, \\ {m}_3=+\frac{a_4{b}_1{e}_{10}}{e_4\Gamma }+\frac{a_2{b}_3{e}_{13}}{e_4\Gamma }-\frac{a_3{b}_2{e}_{13}}{e_4\Gamma }-\frac{a_1{b}_3{e}_{14}}{e_4\Gamma }+\frac{a_3{b}_1{e}_{14}}{e_4\Gamma }-\frac{a_4{b}_2{e}_9}{e_4\Gamma }+\frac{a_2{b}_1{c}_4{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_4{e}_{11}}{c_2{e}_4\Gamma } \\ -\frac{a_4{b}_1{d}_4{e}_{12}}{d_2{e}_4\Gamma }-\frac{a_1{b}_2{d}_3{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_1{b}_3{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{d}_3{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_3{b}_2{c}_1{e}_{15}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_3{e}_{15}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_3{e}_{15}}{c_2{e}_4\Gamma } \\ -\frac{a_2{b}_3{c}_1{e}_{15}}{c_2{e}_4\Gamma }+\frac{a_4{b}_2{c}_1{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_3{b}_1{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{e_7}{e_4}, \\ {m}_4=-\frac{a_4{b}_2{e}_{13}}{e_4\Gamma }+\frac{a_4{b}_1{e}_{14}}{e_4\Gamma }-\frac{a_4{b}_4{d}_4{e}_{16}}{d_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_4{e}_{15}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_4{e}_{15}}{c_2{e}_4\Gamma }+\frac{a_4{b}_2{c}_1{e}_{15}}{c_2{e}_4\Gamma }+\frac{e_8}{e_4}, \\ {m}_5=-\frac{a_1{b}_4{e}_{10}}{e_4\Gamma }+\frac{a_5{b}_1{e}_{10}}{e_4\Gamma }+\frac{a_2{b}_4{e}_9}{e_4\Gamma }-\frac{a_5{b}_2{e}_9}{e_4\Gamma }-\frac{a_1{b}_2{c}_5{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_5{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_2{b}_4{c}_1{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_5{b}_2{c}_1{e}_{11}}{c_2{e}_4\Gamma } \\ -\frac{a_1{b}_2{d}_4{e}_{12}}{d_2{e}_4\Gamma }+\frac{a_1{b}_4{d}_1{e}_{12}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{d}_4{e}_{12}}{d_2{e}_4\Gamma }-\frac{a_5{b}_1{d}_1{e}_{12}}{d_2{e}_4\Gamma }+\frac{e_{11}}{e_4}, \\ {m}_6=-\frac{a_5{b}_1{e}_{14}}{e_4\Gamma }+\frac{a_2{b}_5{e}_9}{e_4\Gamma }-\frac{a_6{b}_2{e}_9}{e_4\Gamma }-\frac{a_1{b}_5{e}_{10}}{e_4\Gamma }+\frac{a_6{b}_1{e}_{10}}{e_4\Gamma }+\frac{a_2{b}_4{e}_{13}}{e_4\Gamma }-\frac{a_5{b}_2{e}_{13}}{e_4\Gamma }-\frac{a_1{b}_4{e}_{14}}{e_4\Gamma } \\ +\frac{a_2{b}_1{d}_5{e}_{12}}{d_2{e}_4\Gamma }-\frac{a_6{b}_1{d}_1{e}_{12}}{d_2{e}_4\Gamma }-\frac{a_1{b}_2{d}_4{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_1{b}_4{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{d}_4{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_1{b}_4{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{d}_4{e}_{16}}{d_2{e}_4\Gamma } \\ -\frac{a_5{b}_1{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_6{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_2{b}_5{c}_1{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_6{b}_2{c}_1{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_5{e}_{15}}{c_2{e}_4\Gamma }-\frac{a_2{b}_4{c}_1{e}_{15}}{c_2{e}_4\Gamma } \\ +\frac{a_2{b}_1{c}_5{e}_{15}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_6{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_5{b}_2{c}_1{e}_{15}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_5{e}_{12}}{c_2{e}_4\Gamma }+\frac{a_1{b}_5{c}_1{e}_{12}}{c_2{e}_4\Gamma }+\frac{e_{18}}{e_4}, \\ {m}_7=-\frac{a_1{b}_6{e}_{10}}{e_4\Gamma }+\frac{a_7{b}_1{e}_{10}}{e_4\Gamma }+\frac{a_2{b}_5{e}_{13}}{e_4\Gamma }-\frac{a_6{b}_2{e}_{13}}{e_4\Gamma }-\frac{a_1{b}_5{e}_{14}}{e_4\Gamma }+\frac{a_6{b}_1{e}_{14}}{e_4\Gamma }+\frac{a_2{b}_6{e}_9}{e_4\Gamma }-\frac{a_7{b}_2{e}_9}{e_4\Gamma } \\ -\frac{a_1{b}_2{d}_5{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_1{b}_5{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{d}_5{e}_{16}}{d_2{e}_4\Gamma }-\frac{a_6{b}_1{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_6{b}_2{c}_1{e}_{15}}{c_2{e}_4\Gamma }+\frac{a_7{b}_2{c}_1{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_6{e}_{15}}{c_2{e}_4\Gamma } \\ -\frac{a_2{b}_5{c}_1{e}_{15}}{c_2{e}_4\Gamma }-\frac{a_2{b}_6{c}_1{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_7{e}_{11}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_7{e}_{11}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_6{e}_{15}}{c_2{e}_4\Gamma }-\frac{a_1{b}_2{c}_6{e}_{12}}{c_2{e}_4\Gamma } \\ +\frac{a_1{b}_6{c}_1{e}_{12}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_6{e}_{12}}{c_2{e}_4\Gamma }-\frac{a_7{b}_1{c}_1{e}_{12}}{c_2{e}_4\Gamma }+\frac{e_{19}}{e_4}, \\ {m}_8=-\frac{a_2{b}_6{e}_{13}}{e_4\Gamma }+\frac{a_7{b}_2{e}_{13}}{e_4\Gamma }+\frac{a_1{b}_6{e}_{14}}{e_4\Gamma }-\frac{a_7{b}_1{e}_{14}}{e_4\Gamma }-\frac{a_1{b}_2{d}_6{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_1{b}_6{d}_1{e}_{16}}{d_2{e}_4\Gamma }+\frac{a_2{b}_1{d}_6{e}_{16}}{d_2{e}_4\Gamma }-\frac{a_7{b}_1{d}_1{e}_{16}}{d_2{e}_4\Gamma } \\ -\frac{a_1{b}_2{c}_7{e}_{15}}{c_2{e}_4\Gamma }+\frac{a_2{b}_1{c}_7{e}_{15}}{c_2{e}_4\Gamma }-\frac{a_2{b}_6{c}_1{e}_{15}}{c_2{e}_4\Gamma }+\frac{a_7{b}_2{c}_1{e}_{15}}{c_2{e}_4\Gamma }+\frac{e_{20}}{e_4}, \\ {n}_2=+\frac{a_2{b}_3{g}_9}{g_4\Gamma }-\frac{a_3{b}_2{g}_9}{g_4\Gamma }-\frac{a_1{b}_3{g}_9}{g_4\Gamma }+\frac{a_3{b}_1{g}_{10}}{g_4\Gamma }+\frac{a_1{b}_3{d}_1{g}_{12}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_3{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_3{b}_1{d}_1{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_1{b}_2{c}_3{g}_{11}}{c_2{g}_4\Gamma } \\ +\frac{a_2{b}_1{c}_3{g}_{11}}{c_2{g}_4\Gamma }-\frac{a_2{b}_3{c}_1{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_3{b}_2{c}_1{g}_{11}}{c_2{g}_4\Gamma }-\frac{a_1{b}_2{d}_3{g}_{12}}{d_2{g}_4\Gamma }+\frac{g_6}{g_4}, \\ {n}_3=+\frac{a_2{b}_3{g}_{13}}{g_4\Gamma }-\frac{a_3{b}_2{g}_{13}}{g_4\Gamma }-\frac{a_1{b}_3{g}_{14}}{g_4\Gamma }+\frac{a_3{b}_1{g}_{14}}{g_4\Gamma }-\frac{a_4{b}_2{g}_9}{g_4\Gamma }+\frac{a_4{b}_1{g}_{10}}{g_4\Gamma }-\frac{a_1{b}_2{c}_4{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_1{b}_2{c}_3{g}_{15}}{c_2{g}_4\Gamma } \\ -\frac{a_4{b}_1{d}_1{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_1{b}_2{d}_3{g}_{16}}{d_2{g}_4\Gamma }+\frac{a_1{b}_3{d}_1{g}_{16}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_3{g}_{16}}{d_2{g}_4\Gamma }-\frac{a_3{b}_1{d}_1{g}_{16}}{d_2{g}_4\Gamma } \\ +\frac{a_2{b}_1{c}_3{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_2{b}_1{c}_4{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_4{b}_2{c}_1{g}_{11}}{c_2{g}_4\Gamma }-\frac{a_2{b}_3{c}_1{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_3{b}_2{c}_1{g}_{15}}{c_2{g}_4\Gamma }+\frac{g_7}{g_4}, \\ {n}_4=-\frac{a_4{b}_2{g}_{13}}{g_4\Gamma }+\frac{a_4{b}_1{g}_{14}}{g_4\Gamma }-\frac{a_1{b}_2{c}_4{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_2{b}_1{c}_4{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_4{b}_2{c}_1{g}_{15}}{c_2{g}_4\Gamma }-\frac{a_4{b}_1{d}_1{g}_{11}}{d_2{g}_4\Gamma }+\frac{g_8}{g_4}, \\ {n}_5=+\frac{a_2{b}_4{g}_9}{g_4\Gamma }-\frac{a_5{b}_2{g}_9}{g_4\Gamma }-\frac{a_1{b}_4{g}_{10}}{g_4\Gamma }+\frac{a_5{b}_1{g}_{10}}{g_4\Gamma }+\frac{a_1{b}_2{d}_4{g}_{12}}{d_2{g}_4\Gamma }+\frac{a_1{b}_4{d}_1{g}_{12}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_4{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_5{b}_1{d}_1{g}_{12}}{d_2{g}_4\Gamma } \\ +\frac{a_1{b}_2{c}_5{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_2{b}_1{c}_5{g}_{11}}{c_2{g}_4\Gamma }-\frac{a_2{b}_4{c}_1{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_5{b}_2{c}_1{g}_{11}}{c_2{g}_4\Gamma }+\frac{g_{17}}{g_4}, \\ {n}_6=+\frac{a_2{b}_4{g}_{13}}{g_4\Gamma }-\frac{a_5{b}_2{g}_{13}}{g_4\Gamma }-\frac{a_1{b}_4{g}_{14}}{g_4\Gamma }+\frac{a_5{b}_1{g}_{14}}{g_4\Gamma }+\frac{a_2{b}_5{g}_9}{g_4\Gamma }-\frac{a_6{b}_2{g}_9}{g_4\Gamma }-\frac{a_1{b}_5{g}_{10}}{g_4\Gamma }+\frac{a_6{b}_1{g}_{10}}{g_4\Gamma } \\ -\frac{a_1{b}_2{d}_4{g}_{16}}{d_2{g}_4\Gamma }+\frac{a_1{b}_4{d}_1{g}_{16}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_4{g}_{16}}{d_2{g}_4\Gamma }-\frac{a_5{b}_1{d}_1{g}_{16}}{d_2{g}_4\Gamma }-\frac{a_2{b}_4{c}_1{g}_{15}}{c_2{g}_4\Gamma }-\frac{a_2{b}_5{c}_1{g}_{11}}{c_2{g}_4\Gamma }-\frac{a_1{b}_2{c}_6{g}_{11}}{c_2{g}_4\Gamma } \\ +\frac{a_2{b}_1{c}_5{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_6{b}_2{c}_1{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_2{b}_1{c}_6{g}_{11}}{c_2{g}_4\Gamma }-\frac{a_1{b}_2{c}_5{g}_{15}}{c_2{g}_4\Gamma }-\frac{a_1{b}_2{d}_5{g}_{12}}{d_2{g}_4\Gamma }+\frac{a_5{b}_2{c}_1{g}_{15}}{c_2{g}_4\Gamma } \\ +\frac{a_1{b}_5{d}_1{g}_{12}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_5{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_6{b}_1{d}_1{g}_{12}}{d_2{g}_4\Gamma }+\frac{g_{18}}{g_4}, \\ {n}_7=+\frac{a_2{b}_6{g}_9}{g_4\Gamma }-\frac{a_7{b}_2{g}_9}{g_4\Gamma }-\frac{a_1{b}_6{g}_{10}}{g_4\Gamma }+\frac{a_7{b}_1{g}_{10}}{g_4\Gamma }+\frac{a_2{b}_5{g}_{13}}{g_4\Gamma }-\frac{a_6{b}_2{g}_{13}}{g_4\Gamma }-\frac{a_1{b}_5{g}_{14}}{g_4\Gamma }+\frac{a_6{b}_1{g}_{14}}{g_4\Gamma } \\ -\frac{a_2{b}_5{c}_1{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_2{b}_1{c}_7{g}_{11}}{c_2{g}_4\Gamma }-\frac{a_2{b}_6{c}_1{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_7{b}_2{c}_1{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_2{b}_1{c}_6{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_6{b}_2{c}_1{g}_{15}}{c_2{g}_4\Gamma }-\frac{a_1{b}_2{c}_6{g}_{15}}{c_2{g}_4\Gamma } \\ -\frac{a_1{b}_2{d}_6{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_1{b}_2{c}_7{g}_{11}}{c_2{g}_4\Gamma }+\frac{a_1{b}_6{d}_1{g}_{12}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_6{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_7{b}_1{d}_1{g}_{12}}{d_2{g}_4\Gamma }-\frac{a_1{b}_2{d}_5{g}_{16}}{d_2{g}_4\Gamma } \\ +\frac{a_1{b}_5{d}_1{g}_{16}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_5{g}_{16}}{d_2{g}_4\Gamma }-\frac{a_6{b}_1{d}_1{g}_{16}}{d_2{g}_4\Gamma }+\frac{g_{19}}{g_4}, \\ {n}_8=-\frac{a_1{b}_6{g}_{14}}{g_4\Gamma }+\frac{a_7{b}_1{g}_{14}}{g_4\Gamma }+\frac{a_2{b}_6{g}_{13}}{g_4\Gamma }-\frac{a_7{b}_2{g}_{13}}{g_4\Gamma }-\frac{a_1{b}_2{d}_6{g}_{16}}{d_2{g}_4\Gamma }+\frac{a_1{b}_6{d}_1{g}_{16}}{d_2{g}_4\Gamma }+\frac{a_2{b}_1{d}_6{g}_{16}}{d_2{g}_4\Gamma }-\frac{a_7{b}_1{d}_1{g}_{16}}{d_2{g}_4\Gamma } \\ -\frac{a_1{b}_2{c}_7{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_2{b}_1{c}_7{g}_{15}}{c_2{g}_4\Gamma }-\frac{a_2{b}_6{c}_1{g}_{15}}{c_2{g}_4\Gamma }+\frac{a_7{b}_2{c}_1{g}_{15}}{c_2{g}_4\Gamma }+\frac{g_{20}}{g_4}, \end{gathered}$$ where (A2)$$\begin{gathered} \Gamma =a_1{b}_2-a_2{b}_1, a_1=-31.00835726A_{11}/L-7.871121054L{A}_{66}/R^2, \\ {a}_2=-21.19550611A_{16}/L-0.01374976696L{A}_{26}/R^2, \\ {a}_3=-1.316898663L{A}_{26}/R^3, a_4=0.05435034699L{A}_{26}/R^3, \\ {a}_5=-629.3381897A_{11}/L^2+11.58524211A_{66}/R^2-11.58524380A_{12}/R^2, \\ {a}_6=-3130.142027A_{11}/L^2-56.62296173A_{66}/R^2+95.01717950A_{12}/R^2, \\ {a}_7=-1642.614350A_{11}/L^2+5.21998485A_{66}/R^2-47.79767774A_{12}/R^2, \\ {b}_1=-0.5428190419A_{16}/L-0.01374966901L{A}_{26}/R^2, \\ {b}_2=-144.4460509A_{66}/L-0.7841786683L{A}_{22}/R^2-3.136714672KL{A}_{44}/R^2, \\ {b}_3=-0.1641087045L{A}_{22}/R^3, b_4=58.03431294A_{26}/R^2, \\ {b}_5=-118.4772154A_{16}/L^2-49.61690248A_{26}/R^2, \\ {b}_6=2185.537826A_{16}/L^2+171.2842070A_{26}/R^2, \\ {c}_1=-31.00835726B_{11}/L-7.871121054L{B}_{66}/R^2, \\ {c}_2=-3.141803460L{A}_{55}K-31.00835727D_{11}/L-7.871121053L{D}_{66}/R^2, \\ {c}_3=-1.316898663L{B}_{26}/R^3, c_4=0.05435034702L{B}_{26}/R^3, \\ {c}_5=-629.3381897B_{11}/L^2+11.58524211B_{66}/R^2-11.58524380B_{12}/R^2, \\ {c}_6=-3130.142024B_{11}/L^2-56.62296171B_{66}/R^2+95.01717948B_{12}/R^2, \\ {c}_7=-1642.614350B_{11}/L^2+5.21998492B_{66}/R^2-47.79767769B_{12}{R}^2, \\ {d}_1=-144.4460509B_{66}/L-0.7841786683L{B}_{22}/R^2-3.136714672KL{B}_{44}/R^2, \\ {d}_2=-0.5428190419B_{16}/L-0.01374966901L{B}_{26}/R^2, \\ {d}_3=-0.1641087045L{B}_{22}/R^3, d_4=58.03431294B_{26}/R^2, \\ {d}_5=-118.4772154B_{16}/L^2-49.61690248B_{26}/R^2, \\ {d}_6=2185.537826B_{16}/L^2+171.2842070B_{26}/R^2, e_1=25.09540907, \\ {e}_2=288.9114871/L, e_3=-6.273850301L\gamma , e_4=-6.273850303L{I}_0, \\ {e}_5=-288.9114871A_{55}K/L-6.273850307L{A}_{22}/R^2-15.71779897L{A}_{44}K/R^2 \\ -288.9114871N_{xx}^T/L-15.71779897L{N}_{\theta \theta }^T/R^2, \\ {e}_6=-39.72963636A_{26}/R^3, e_7=-165.6765179A_{26}/R^3, e_8=-103.2852305A_{26}/R^3, \\ {e}_9=545.8105948A_{11}/L^2-6.64989922A_{66}/R^2+6.649896554A_{12}/R^2, \\ {e}_{10}=-58.03431074A_{26}/R^2, e_{12}=-58.03431074B_{26}/R^2, \\ {e}_{11}=545.8105948B_{11}/L^2-6.64989922B_{66}/R^2+6.649896554B_{12}/R^2, \\ {e}_{13}=-182.5445517A_{11}/L^2-143.4703132A_{66}/R^2-8.169858374A_{12}/R^2, \\ {e}_{14}=-640.3323995A_{16}/L^2+17.60816090A_{26}/R^2, \\ {e}_{15}=-182.5445517B_{11}/L^2-143.4703132B_{66}/R^2-8.169858374B_{12}/R^2, \\ {e}_{16}=-640.3323995B_{16}/L^2+17.60816090B_{26}/R^2, \\ {e}_{17}=-18183.70969A_{11}/L^3-18.0373938L{A}_{22}/R^4+\left(197.449892A_{12}+394.89978A_{66}\right)/L{R}^2, \\ {e}_{18}=1209.128345A_{11}/L^3+5.069883616L{A}_{22}/R^4-\left(242.5993347A_{12}+485.19867A_{66}\right)L{R}^2, \\ {e}_{19}=74011.22204A_{11}/L^3-59.75908863L{A}_{22}/R^4-\left(356.845848A_{12}+2841.38288A_{66}\right)/L{R}^2, \\ {e}_{20}=9504.768786A_{11}/L^3+1.592447035L{A}_{22}/R^4-\left(21.736372A_{12}+43.4727439A_{66}\right)/L{R}^2, \\ {g}_1=25.0954091, g_2=1078.441452/L, g_3=-6.285173437L\gamma , g_4=-6.285173437L{I}_0, \\ {g}_5=-1078.441451A_{55}K/-6.285173436L{A}_{22}/R^2-42.00198451L{A}_{44}K/R^2 \\ -1078.441451N_{xx}^T/L-42.00198451L{N}_{\theta \theta }^T/R^2, \\ {g}_6=-49.75797336A_{26}/R^3, g_7=-165.6765179A_{26}/R^3, g_8=-62.58299370A_{26}/R^3, \\ {g}_9=-182.5480497A_{11}/L^2-27.31625097A_{66}/R^2+8.169718595A_{12}/R^2, \\ {g}_{10}=-120.4374000A_{26}/R^2, g_{12}=-120.4374000B_{26}/R^2, \\ {g}_{11}=-182.5480497B_{11}/L^2-27.31625097B_{66}/R^2+8.169718595B_{12}/R^2, \\ {g}_{13}=1188.606609A_{11}/L^2-3.31779223A_{66}/R^2+30.35534344A_{12}/R^2, \\ {g}_{14}=-1536.656782A_{16}/L^2-120.4374000A_{26}/R^2, \\ {g}_{15}=1188.606609B_{11}/L^2-3.31779223B_{66}/R^2+30.35534344B_{12}/R^2, \\ {g}_{16}=-1536.656782B_{16}/L^2-120.4374000B_{26}/R^2, \\ {g}_{17}=-180.6408093A_{11}/L^3+1.690008865L{A}_{22}/R^4-\left(33.432213A_{12}+66.864426A_{66}\right)/L{R}^2, \\ {g}_{18}=-14810.2121A_{11}/L^3-59.7583037L{A}_{22}/R^4-\left(1241.088955A_{12}+2482.1779A_{66}\right)/L{R}^2, \\ {g}_{19}=1585.526437A_{11}/L^3+4.777494122L{A}_{22}/R^4-\left(291.12256A_{12}+581.5127135A_{66}\right)/L{R}^2, \\ {g}_{20}=-23540.79961A_{11}/L^3-119.1995L{A}_{22}/R^4-\left(119.3011A_{12}+238.60223A_{66}\right)/L{R}^2 \end{gathered}$$

References

1. Montesano, J.; Bougherara, H.; Fawaz, Z. Application of infrared thermography for the characterization of damage in braided carbon fiber reinforced polymer matrix composites. Compos. Part B Eng.; 2014; 60, pp. 137-143. [DOI: https://dx.doi.org/10.1016/j.compositesb.2013.12.053]

2. Zhu, J.H.; Wei, L.; Wang, Z.; Liang, C.K.; Fang, Y. Application of carbon fiber reinforced polymer anode in electrochemical chloride extraction of steel-reinforced concrete. Constr. Build. Mater.; 2016; 120, pp. 275-283. [DOI: https://dx.doi.org/10.1016/j.conbuildmat.2016.05.103]

3. Zhu, J.H.; Zhu, M.; Han, N.; Liu, W.; Xing, F. Electrical and Mechanical Performance of Carbon Fiber-Reinforced Polymer Used as the Impressed Current Anode Material. Materials; 2014; 7, pp. 5438-5453. [DOI: https://dx.doi.org/10.3390/ma7085438] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28788137]

4. Asadi, H.; Wang, Q. Dynamic stability analysis of a pressurized FG-CNTRC cylindrical shell interacting with supersonic airflow. Compos. Part B Eng.; 2017; 118, pp. 15-25. [DOI: https://dx.doi.org/10.1016/j.compositesb.2017.03.001]

5. Zhang, W.; Yang, S.W.; Mao, J.J. Nonlinear radial breathing vibrations of CFRP laminated cylindrical shell with non-normal boundary conditions subjected to axial pressure and radial line load at two ends. Compos. Struct.; 2018; 190, pp. 52-78. [DOI: https://dx.doi.org/10.1016/j.compstruct.2018.01.091]

6. Yang, S.W.; Zhang, W.; Mao, J.J. Nonlinear vibrations of carbon fiber reinforced polymer laminated cylindrical shell under non-normal boundary conditions with 1:2 internal resonance. Eur. J. Mech./A Solids; 2019; 74, pp. 317-336. [DOI: https://dx.doi.org/10.1016/j.euromechsol.2018.11.014]

7. Kiral, B.G.; Kiral, Z.; Ozturk, H. Stability analysis of delaminated composite beams. Compos. Struct.; 2015; 93, pp. 342-350.

8. Ke, L.L.; Wang, Y.S. Size effect on dynamic stability of functionally graded microbeams based on a modified couple stress theory. Compos. Struct.; 2011; 79, pp. 406-418. [DOI: https://dx.doi.org/10.1016/j.compstruct.2010.09.008]

9. Couto, C.; Maia, E.; Real, P.V.; Lopes, N. The effect of non-uniform bending on the lateral stability of steel beams with slender cross-section at elevated temperatures. Eng. Struct.; 2018; 163, pp. 153-166. [DOI: https://dx.doi.org/10.1016/j.engstruct.2018.02.033]

10. Talebitooti, M. Analytical and finite-element solutions for the buckling of composite sandwich conical shell with clamped ends under external pressure. Arch. Appl. Mech.; 2017; 87, pp. 59-73. [DOI: https://dx.doi.org/10.1007/s00419-016-1176-y]

11. Maali, M.; Showkati, H.; Fatemi, S.M. Investigation of the buckling behavior of conical shells under weld-induced imperfections. Thin-Walled Struct.; 2012; 57, pp. 13-24. [DOI: https://dx.doi.org/10.1016/j.tws.2012.04.003]

12. Bich, D.H.; Phuong, N.T.; Tung, H.V. Buckling of functionally graded conical panels under mechanical loads. Compos. Struct.; 2012; 94, pp. 1379-1384. [DOI: https://dx.doi.org/10.1016/j.compstruct.2011.11.029]

13. Gajdzicki, M.; Perlinski, W.; Michalak, B. Stability analysis of bi-directionally corrugated steel plates with orthotropic plate model. Eng. Struct.; 2018; 160, pp. 519-534. [DOI: https://dx.doi.org/10.1016/j.engstruct.2018.01.057]

14. Zeng, H.C.; Huang, C.S.; Leissa, A.W.; Chang, M.J. Vibrations and stability of a loaded side-cracked rectangular plate via the MLS-Ritz method. Thin-Walled Struct.; 2016; 106, pp. 459-470. [DOI: https://dx.doi.org/10.1016/j.tws.2016.05.013]

15. Dey, T.; Ramachandra, L.S. Static and dynamic unstable analysis of composite cylindrical shell panels subjected to partial edge loading. Int. J. Non-Linear Mech.; 2014; 64, pp. 46-56. [DOI: https://dx.doi.org/10.1016/j.ijnonlinmec.2014.03.014]

16. Han, Q.K.; Qin, Z.Y.; Lu, W.X.; Chu, F.L. Dynamic stability analysis of periodic axial loaded cylindrical shell with time-dependent rotating speeds. Nonlinear Dyn.; 2015; 81, pp. 1649-1664. [DOI: https://dx.doi.org/10.1007/s11071-015-2097-y]

17. Kolanua, N.R.; Rajub, G.; Ramji, M. Experimental and numerical studies on the buckling and post-buckling behavior of single blade-stiffened CFRP panels. Compos. Struct.; 2018; 196, pp. 135-154. [DOI: https://dx.doi.org/10.1016/j.compstruct.2018.05.015]

18. Salvetti, M.; Gilioli, A.; Sbarufatti, C.; Manes, A.; Giglio, M. Analytical model of the dynamic behaviour of CFRP plates subjected to low velocity impacts. Compos. Part B Eng.; 2018; 142, pp. 47-55. [DOI: https://dx.doi.org/10.1016/j.compositesb.2018.01.005]

19. Cui, J.; Dong, D.; Zhang, X.; Huang, X.; Lu, G. Influence of thickness of composite layers on failure behaviors of carbon fiber reinforced plastics/aluminum alloy electromagnetic riveted lap joints under high-speed loading. Int. J. Impact Eng.; 2018; 115, pp. 1-9. [DOI: https://dx.doi.org/10.1016/j.ijimpeng.2018.01.004]

20. Zhang, F.S.; Xu, J.Z.; Xu, G.D.; Zu, L. The buckling behavior of radome with different braided angles based on CFRP. Compos. Struct.; 2017; 176, pp. 597-607. [DOI: https://dx.doi.org/10.1016/j.compstruct.2017.05.067]

21. Juntanalikit, P.; Jirawattanasomkul, T.; Pimanmas, A. Experimental and numerical study of strengthening non-ductile RC columns with and without lap splice by Carbon Fiber Reinforced Polymer (CFRP) jacketing. Eng. Struct.; 2016; 125, pp. 400-418. [DOI: https://dx.doi.org/10.1016/j.engstruct.2016.07.019]

22. Reuter, C.; Sauerland, K.H.; Tröster, T. Experimental and numerical crushing analysis of circular CFRP tubes under axial impact loading. Compos. Struct.; 2017; 174, pp. 33-44. [DOI: https://dx.doi.org/10.1016/j.compstruct.2017.04.052]

23. Timme, S.; Trappe, V.; Korzen, M.; Schartel, B. Fire stability of carbon fiber reinforced polymer shells on the intermediate-scale. Compos. Struct.; 2017; 178, pp. 320-329. [DOI: https://dx.doi.org/10.1016/j.compstruct.2017.07.025]

24. Zhang, W.; Zhao, M.H. Nonlinear vibrations of a composite laminated cantilever rectangular plate with one-to-one internal resonance. Nonlinear Dyn.; 2012; 70, pp. 295-313. [DOI: https://dx.doi.org/10.1007/s11071-012-0455-6]

25. Zhao, M.H.; Zhang, W. Nonlinear dynamics of composite laminated cantilever rectangular plate subject to third-order piston aerodynamics. Acta Mech.; 2014; 225, pp. 1985-2004. [DOI: https://dx.doi.org/10.1007/s00707-013-1035-7]

26. Zhang, W.; Zhao, M.H.; Guo, X.Y. Nonlinear responses of a symmetric cross-ply composite laminated cantilever rectangular plate under in-plane and moment excitations. Compos. Struct.; 2013; 100, pp. 554-565. [DOI: https://dx.doi.org/10.1016/j.compstruct.2013.01.013]

27. Hwu, C.; Hsu, H.W.; Lin, Y.H. Free vibration of composite sandwich plates and cylindrical shells. Compos. Struct.; 2017; 171, pp. 528-537. [DOI: https://dx.doi.org/10.1016/j.compstruct.2017.03.042]

28. Viswanathan, K.; Javed, S. Free vibration of anti-symmetric angle-ply cylindrical shell walls using first-order shear deformation theory. J. Vib. Control; 2016; 22, pp. 1757-1768. [DOI: https://dx.doi.org/10.1177/1077546314544893]

29. Sarkheil, S.; Foumani, M.S. Free vibrational characteristics of rotating joined cylindrical-conical shells. Thin-Walled Struct.; 2016; 107, pp. 657-670. [DOI: https://dx.doi.org/10.1016/j.tws.2016.07.009]

30. Amabili, M.; Balasubramanian, P.; Ferrari, G. Travelling wave and non-stationary response in nonlinear vibrations of water-filled circular cylindrical shells: Experiments and simulations. J. Sound Vib.; 2016; 381, pp. 220-245. [DOI: https://dx.doi.org/10.1016/j.jsv.2016.06.026]

31. Song, Z.G.; Zhang, L.W.; Liew, K.M. Vibration analysis of CNT-reinforced functionally graded composite cylindrical shells in thermal environments. Int. J. Mech. Sci.; 2016; 115–116, pp. 339-347. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2016.06.020]

32. Zhang, W.; Hao, Y.X.; Yang, J. Nonlinear dynamics of FGM circular cylindrical shell with clamped-clamped edges. Compos. Struct.; 2012; 94, pp. 1075-1086. [DOI: https://dx.doi.org/10.1016/j.compstruct.2011.11.004]

33. Du, C.C.; Li, Y.H. Nonlinear resonance behavior of functionally graded cylindrical shells in thermal environments. Compos. Struct.; 2013; 102, pp. 164-174. [DOI: https://dx.doi.org/10.1016/j.compstruct.2013.02.028]

34. Sun, Y.; Zhang, W.; Yao, M.H. Multi-pulse chaotic dynamics of circular mesh antenna with 1:2 internal resonance. Int. J. Appl. Mech.; 2017; 9, 1750060. [DOI: https://dx.doi.org/10.1142/S1758825117500600]

35. Liu, T.; Zhang, W.; Wang, J.F. Nonlinear dynamics of composite laminated circular cylindrical shell clamped along a generatrix and with membranes at both ends. Nonlinear Dyn.; 2017; 90, pp. 1393-1417. [DOI: https://dx.doi.org/10.1007/s11071-017-3734-4]

36. Wang, Y.Q. Nonlinear vibration of a rotating laminated composite circular cylindrical shell: Traveling wave vibration. Nonlinear Dyn.; 2014; 77, pp. 1693-1707. [DOI: https://dx.doi.org/10.1007/s11071-014-1410-5]

37. Hao, Y.X.; Yang, S.W.; Zhang, W.; Yao, M.H.; Wang, A.W. Flutter of high-dimension nonlinear system for a FGM truncated conical shell. Mech. Adv. Mater. Struct.; 2018; 25, pp. 47-61. [DOI: https://dx.doi.org/10.1080/15376494.2016.1255815]

38. Wang, Y.Q.; Guo, X.H.; Chang, H.H.; Li, H.Y. Nonlinear dynamic response of rotating circular cylindrical shells with precession of vibrating shape, Part I: Numerical solution. Int. J. Mech. Sci.; 2010; 52, pp. 1217-1224. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2010.05.008]

39. Wang, Y.Q.; Guo, X.H.; Chang, H.H.; Li, H.Y. Nonlinear dynamic response of rotating circular cylindrical shells with precession of vibrating shape, Part II: Approximate analytical solution. Int. J. Mech. Sci.; 2010; 52, pp. 1208-1216. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2010.05.007]

40. Shen, H.S. Nonlinear vibration of shear deformable FGM cylindrical shells surrounded by an elastic medium. Compos. Struct.; 2012; 94, pp. 1144-1154. [DOI: https://dx.doi.org/10.1016/j.compstruct.2011.11.012]

41. Shen, H.S.; Wang, H. Nonlinear vibration of shear deformable FGM cylindrical panels resting on an elastic medium in thermal environments. Compos. Part B Eng.; 2014; 60, pp. 167-177. [DOI: https://dx.doi.org/10.1016/j.compositesb.2013.12.051]

42. Reddy, J.N. Mechanics of Laminated Composite Plates and Shells: Theory and Analysis; CRC Press: New York, NY, USA, 2004.

43. Hao, Y.X.; Zhang, W.; Yang, L.; Wang, J.H. Dynamic Response of Cantilever FGM Cylindrical Shell. Appl. Mech. Mater.; 2012; 130, pp. 3986-3993. [DOI: https://dx.doi.org/10.4028/www.scientific.net/AMM.130-134.3986]

44. Efraim, E.; Eisenberger, M. Exact vibration analysis of variable thickness thick annular isotropic and FGM plates. J. Sound Vib.; 2007; 299, pp. 720-738. [DOI: https://dx.doi.org/10.1016/j.jsv.2006.06.068]

45. Noseir, A.; Reddy, J.N. A study of non-linear dynamic equations of higher-order deformation plate theories. Int. J. Non-Linear Mech.; 1991; 26, pp. 233-249. [DOI: https://dx.doi.org/10.1016/0020-7462(91)90054-W]

46. Bhimaraddi, A. Large amplitude vibrations of imperfect antisymmetric angle-ply laminated plates. J. Sound Vib.; 1993; 162, pp. 457-470. [DOI: https://dx.doi.org/10.1006/jsvi.1993.1133]

47. Zhang, X.M.; Liu, G.R.; Lam, K.Y. Vibration analysis of thin cylindrical shells using wave propagation approach. J. Sound Vib.; 2001; 239, pp. 397-403. [DOI: https://dx.doi.org/10.1006/jsvi.2000.3139]

48. Song, Z.G.; Li, F.M. Aerothermoelastic analysis and active flutter control of supersonic composite laminated cylindrical shells. Compos. Struct.; 2013; 106, pp. 653-660. [DOI: https://dx.doi.org/10.1016/j.compstruct.2013.07.029]

49. Lee, D.S. Nonlinear dynamic buckling of orthotropic cylindrical shells subjected to rapidly applied loads. J. Eng. Math.; 2000; 38, pp. 141-154. [DOI: https://dx.doi.org/10.1023/A:1004729816574]

50. Gao, K.; Gao, W.; Wu, D.; Song, C. Nonlinear dynamic stability of the orthotropic functionally graded cylindrical shell surrounded by Winkler-Pasternak elastic foundation subjected to a linearly increasing load. J. Sound Vib.; 2018; 415, pp. 147-168. [DOI: https://dx.doi.org/10.1016/j.jsv.2017.11.038]

51. Darabi, M.; Ganesan, R. Non-linear dynamic unstable analysis of laminated composite cylindrical shells subjected to periodic axial loads. Compos. Struct.; 2016; 147, pp. 168-184. [DOI: https://dx.doi.org/10.1016/j.compstruct.2016.02.064]

52. Lei, Z.X.; Zhang, L.W.; Liew, K.M.; Yu, J.L. Dynamic stability analysis of carbon nanotubereinforced functionally graded cylindrical panels using the element-free kp-Ritz method. Compos. Struct.; 2014; 113, pp. 328-338. [DOI: https://dx.doi.org/10.1016/j.compstruct.2014.03.035]

53. Ganapathi, M.; Balamurugan, V. Dynamic unstable analysis of laminated composite circular cylindrical shell. Comput. Struct.; 1998; 69, pp. 181-189. [DOI: https://dx.doi.org/10.1016/S0045-7949(98)00095-9]

54. Dey, T.; Ramachandra, L.S. Dynamic stability of simply supported composite cylindrical shells under partial axial loading. J. Sound Vib.; 2015; 353, pp. 272-291. [DOI: https://dx.doi.org/10.1016/j.jsv.2015.05.021]