Free Vibration Analysis of a Spinning Composite Laminated Truncated Conical Shell under Hygrothermal Environment

: This paper is concerned with free vibration characteristics of a spinning composite laminated truncated conical shell subjected to hygrothermal environment. Hygrothermal strains are introduced into the constitutive law of single-layer material, and ﬁber orientation lies symmetrically with respect to the midplane of the composite shell. Considering the spin-induced Coriolis and centrifugal forces, as well as initial hoop tension, the governing equations of free vibration of the composite conical shell with hygrothermal effects are derived on the basis of Love’s thin-shell theory and Hamilton’s principle. The solution of the equations is derived using the Galerkin approach. Then, a detailed parametric study on natural frequencies and critical spinning speeds is numerically performed. Results indicate that the Coriolis force induces an asymmetric inﬂuence on natural frequencies of forward and backward traveling waves, while the centrifugal force enhances the frequencies of both traveling waves symmetrically. Initial hoop tension plays a major role in the increase of critical spinning angular speed. Temperature, moisture concentration, and design parameters show the signiﬁcant inﬂuence on the free vibration characteristics of the conical shell, and thermal expansion deformation is nonnegligible in the free vibration analysis.


Introduction
Conical shell structures spinning about their symmetrical axis have been extensively utilized in numerous industrial applications, such as high-speed centrifugal separators, aircraft propellers, gas turbine engines, and so on. The Coriolis and centrifugal forces caused by spinning motion induce the traveling wave modes [1], which may result in complex vibration characteristics of the conical shell. Recently, composite materials have become of much interest to engineers and researchers, owing to their high strength and light weight [2]. However, composite shell structures are highly susceptible to a hygrothermal environment, which is a dominant influencing factor on vibration characteristics of the spinning conical shell made of composite. As a result, it is essential to know the dynamic behaviors of spinning composite conical shells under a hygrothermal environment.
Based on elastic theory, research on vibration characteristics of spinning shell structures has been performed for a long time. Pioneering work on spinning shell structures was reported in [3], which presented the phenomenon of traveling wave vibration of spinning cylindrical shells. After that, Srinivasan et al. [4] discussed the effects of centrifugal and Coriolis forces on traveling wave vibration of the spinning cylindrical shell, and found that the Coriolis force has an asymmetric effect on the traveling wave vibration, while the centrifugal force induces the symmetric effect on the vibration. Since then, abundant studies on vibration behavior of spinning conical shells have been carried out [5,6]. Various numerical and analytical methods have been developed to investigate the vibration characteristics of spinning conical shells, such as the generalized differential quadrature method [7], Galerkin method [8], Ritz method [9,10], finite element method [11], harmonic balance method [12], and multiple scales method [13]. Donnell's shell theory [14], Sanders' shell theory [15], Love's shell theory [16,17], first-order shear deformation theory [18], and higher-order shear deformation theory [19] are also widely utilized in dynamic modeling of the shell structures.
In recent years, the vibration behaviors of composite laminated beams [20], plates [21,22], and shells [23,24] have been of much interest to many engineers and researchers. Among them, composite laminated conical shells with spinning motion are normally paid much interest, especially in the following papers. Shakouri [25] analyzed the effect of asymmetrical distribution of functionally graded material on vibration characteristics of spinning conical shells, Civalek [26] investigated vibration behavior of spinning conical shells made of symmetrical laminated material, Heydarpour et al. [27] studied free vibration of the shells made of carbon nanotube-reinforced composite material, and Li et al. [28] discussed the influence of symmetrical and asymmetrical graphene dispersion patterns on dynamic stability of spinning functionally graded graphene reinforced conical shells. These works have improved our knowledge of spinning conical shells made of composites.
It is worth mentioning that most of the existing works focus on composite shell structures considering hygrothermal effects. In practical engineering, these structures are usually exposed in complex environments with strong hygrothermal effects, which can result in the variation of mechanical properties of the composites and then affect the vibration behaviors of the structures [29][30][31]. For example, some accidents of wind turbines happen in winter, which is because remarkable cooling tremendously changes mechanical properties of composite blades consisting of tapered cylindrical shells, and may result in resonance between blades and towers. In this case, study on hygrothermal effects of composite shell structures is of prime importance. For this purpose, Patel et al. [32] proposed a unified formulation on mechanical behaviors of composite laminated plates in different combinations of temperature and moisture concentration. In [33], an experimental and numerical study on vibration behaviors of composite shallow shells subjected to hygrothermal load was conducted, and the results indicated that temperature and moisture both impaired the frequencies of the shells. Li [34] analyzed the effects of hygrothermal conditions on nonlinear parametric resonances of laminated cylindrical shells. Qin et al. [35,36] investigated the effects of aerodynamic force and installation mode on hygrothermal vibration of rotating laminated composite beams. Bandyopadhyay et al. [37] investigated transient dynamic response of composite conical shells under hygrothermal load and low-velocity impact. Hajmohammad et al. [38] studied dynamic bucking of carbon nanotube-reinforced truncated conical shells, and found that dynamic instability region shifts to lower frequencies with increasing temperature and moisture. Parhi et al. [39] dealt with nonlinear vibration of composite conical shells in a hygrothermal environment, and found that hygrothermal load induce instability of the conical shell.
As mentioned above, significant works on vibration behaviors of spinning shell structures have been published in the literature. However, to the best of the authors' knowledge, no attention has been given to free vibration characteristics of spinning laminated composite truncated conical shells subjected to a hygrothermal environment. The novelty of this paper is the focus on the combined effect of hygrothermal conditions and spinning motion on free vibration of the composite conical shells. Some new and interesting findings are proposed by investigating the mechanism of how spin-induced Coriolis force, centrifugal force, initial hoop tension, hygrothermal effects, and design parameters affect free vibration behaviors, which can be utilized as a benchmark for optimization of design and vibration control of spinning composite truncated conical shells, considering hygrothermal effects in future research.

Constitutive Equations
In Figure 1, a composite truncated conical shell that spins about its horizontal and symmetrical axis at a constant angular speed Ω is presented, where L, h, and α denote the length, thickness, and semivertex angle of the shell, respectively. An orthogonal coordinate system (x, θ, z) is fixed at the middle surface of the conical shell. The displacement components along the x, θ, and z directions are denoted by u(x, θ, z, t), v(x, θ, z, t), and w(x, θ, z, t), respectively. R 1 and R 2 in the figure refer to the mean radii at both ends, and the cone radius of any point in the x direction is expressed as R(x) = R 1 + xsinα. timization of design and vibration control of spinning composite truncated conical shells, considering hygrothermal effects in future research.

Constitutive Equations
In Figure 1, a composite truncated conical shell that spins about its horizontal and symmetrical axis at a constant angular speed Ω is presented, where L, h, and α denote the length, thickness, and semivertex angle of the shell, respectively. An orthogonal coordinate system (x, θ, z) is fixed at the middle surface of the conical shell. The displacement components along the x, θ, and z directions are denoted by u(x, θ, z, t), v(x, θ, z, t), and w(x, θ, z, t), respectively. R1 and R2 in the figure refer to the mean radii at both ends, and the cone radius of any point in the x direction is expressed as R(x) = R1 + xsinα. According to Love's thin-shell theory [40], the strains of the conical shell in terms of the displacement components u(x, θ, z, t), v(x, θ, z, t) and w(x, θ, z, t) are expressed as follows: where ε1, ε2, and γ represent the middle surface strains of the shell, and κ1, κ2, and χ are the middle surface curvatures. They are given as follows: where R = R(x) is the cone radius of the middle surface. Considering hygrothermal expansion deformation, the constitutive equation of a single-layer composite conical shell is given by: where ε = (εx, εθ, εxθ) t and σ = (σx, σθ, σxθ) t represent total strain and stress vectors of the single-layer composite, the superscript t stands for transposition of the vectors. ε T and ε H According to Love's thin-shell theory [40], the strains of the conical shell in terms of the displacement components u(x, θ, z, t), v(x, θ, z, t) and w(x, θ, z, t) are expressed as follows: where ε 1 , ε 2 , and γ represent the middle surface strains of the shell, and κ 1 , κ 2 , and χ are the middle surface curvatures. They are given as follows: where R = R(x) is the cone radius of the middle surface. Considering hygrothermal expansion deformation, the constitutive equation of a single-layer composite conical shell is given by: where ε = (ε x , ε θ , ε xθ ) t and σ = (σ x , σ θ , σ xθ ) t represent total strain and stress vectors of the single-layer composite, the superscript t stands for transposition of the vectors. ε T and ε H refer to the thermal and humid strains, respectively, and Q is the stiffness matrix with the following elements [41]: where θ f represents fiber orientation angle defined as the angle from x-direction to the fiber orientation. In this study, fiber orientation lies symmetrically with respect to the midplane of the composite shell.
, and Q 66 = G 12 , in which E 1 , E 2 are elastic moduli in principal orientation, G 12 shear modulus, and v 12 , v 21 Poisson ratios. The thermal strain in Equation (3) generated by thermal expansion deformation is expressed as: t . α 12 = (α 1 , α 2 , 0) t , and α 1 , α 2 are the thermal expansion coefficients along the principal orientation. ∆T is temperature variation, and T t the transformation matrix given by the following form [41]: The humid strain in Equation (3) caused by hygroscopic expansion deformation is given by: where ε H = (ε H x , ε H θ , ε H xθ ) t , β 12 = (β 1 , β 2 , 0) t , and β 1 , β 2 are the hygroscopic expansion coefficients. ∆C represents the moisture concentration, written as: where Based on mechanics of composite materials [41], the in-surface force and moment resultants of the composite laminated truncated conical shell can be expressed as: Substituting Equation (3) into Equation (9), the hygrothermal loads due to the hygrothermal expansion deformation can be obtained:

Governing Equations
The strain energy of spinning composite laminated truncated conical shells consists of two contributions, which correspond to strain energy U D generated by initial strains and strain energy U h caused by spinning centrifugal force. They are formulated as follows: In which N 0 θ denotes initial hoop tension due to the spinning centrifugal force, defined as: where ρ is the density of the single-layer composite. Considering spinning motion, the displacement vector r and the speed vector v at any point of the spinning conical shell take the form of: In which i, j and k are the unit vectors along the x, θ, and z directions and Ω = −Ω cos αi + Ω sin αk.
Then, the kinetic energy of the spinning truncated conical shell in arbitrary time is expressed as: where ρ t is the density per unit area, described as: For the sake of deriving governing equations of the spinning composite laminated truncated conical shell in a hygrothermal environment, Hamilton's principle is utilized with the following form: where W is the work performed by external loads and δW = 0 in this analysis. Inserting Equations (11), (12) and (15) into Equation (17), the governing equations of the conical shell can be obtained: where L ij (i, j = 1, 2, 3) are differential operators listed in Appendix A.

Solution Procedure
In this subsection, the Galerkin approach is utilized to derive the partial differential Equation (18). The displacement functions of the conical shells along the x, θ, and z directions are assumed as [42,43]: where q mn (t), p mn (t) represent generalized displacement, and α mn and β mn are mode shape ratios [44]. ϕ im (x)(I = u, v, w) represent the modal functions, m is the number of axial half waves, and n is the number of circumferential waves.
Considering a simply supported truncated conical shell, the boundary condition at both small and large edges (Ss-Sl) are expressed as: Based on Equation (20), the corresponding modal functions of simply supported conical shells are defined as Introducing Equation (19) into Equation (18) and multiplying mode shapes of Equation (19), respectively, integrating over the x-θ plane yields: where I = 1, 2, . . . , mn. M i represents mass coefficient of the conical shell for the i-th mode (m, n), C i is damping coefficients generated by the Coriolis force owing to spinning motion, K i is stiffness coefficient of the conical shell, and K ri is stiffness coefficient caused by spinning motion. They are listed in Appendix B. For brevity, Equation (21) can be expressed as the following standard eigenvalue problem: It is obvious that λ represents the characteristic root of Equation (22), which includes real and imaginary parts. The real part is three coupled roots corresponding to each combination of m and n of a specific truncated conical shell, related to natural frequencies of u, v, and w in each mode. Among them, the frequencies related to w are normally paid more attention because other frequencies are much higher than the ones for w. When the conical shell undergoes the spinning motion, the Coriolis force appears and induces the natural mode of the shell bifurcating into forward and backward traveling wave modes. Usually, the positive root related to w refers to natural frequency of a backward traveling wave, while the negative root represents the frequency of a forward traveling wave.

Results and Discussion
In this section, a detailed parametric study on free vibration characteristics of spinning composite laminated truncated conical shells in a hygrothermal environment is numerically performed. The graphite-epoxy composite laminated conical shell is adopted in the numerical simulation, and the variation in material properties of the composite shell with respect to temperature and moisture concentrations is presented in Table 1. A symmetric (θ f/ − θ f ) s lay up of laminates is considered, and the density of the single-layer composite is 1600 kg/m 3 . ∆T in Equation (5) is defined as T − T 0 , and reference temperature T 0 is assumed to be 25 • C. Unless otherwise specified, the fiber orientation angle θ f = 60 • , L/R 1 = 5, and h/R 1 = 0.01 obtain in the following. Table 1. Variation of elastic modulus of graphite-epoxy laminate [32].
It can be seen in Table 1 that elastic modulus E 1 is unaffected by temperature or moisture concentration in a certain range. The reason is that the elastic modulus of graphite fiber, which has the predominant contribution in elastic modulus E 1 of the conical shell, is much higher than that of the epoxy matrix, and the graphite fiber is negligibly affected by temperature and moisture. However, the epoxy matrix is the dominant influencing factor on elastic modulus E 2 of the shell. It is highly susceptible to temperature variation and moisture concentration, and the influence of temperature on the matrix is more prominent than that of moisture. Moreover, elastic modulus G 12 of the shell is related to both the graphite fiber and the epoxy matrix. Specifically, the graphite fiber and epoxy matrix are less affected by moisture, which demonstrates that G 12 is unaffected by the moisture. However, the matrix is greatly affected by temperature, which indicates that G 12 is reduced with the increase in temperature.
In what follows, the dimensionless values are introduced: where ω*, ω 1 *, ω 2 *, and Ω* stand for dimensionless natural frequencies of the conical shell, frequencies of the forward traveling wave, frequencies of the backward traveling wave, and spinning angular speed, respectively.

Validation
In order to validate the accuracy of the present theory, the dimensionless natural frequencies of truncated conical shells with Ss-Sl boundary condition are compared with results in the literature, as listed in Tables 2 and 3. Table 2 compares the frequency ω* of nonspinning conical shells for different semivertex angles. Table 3 compares the frequencies of forward and backward traveling waves of isotropic spinning conical shells, which correspond to ω 1 * and ω 2 *. It is observed that the present results agree well with those in the literature. For the sake of demonstrating the validity of the present solution considering hygrothermal effects, the natural frequencies of composite truncated conical shells are calculated and compared with the results of the finite element method (FEM). In the modeling process, the quadrilateral shell elements with four nodes are used in Nastran software. The comparison results are presented in Table 4, where the relative error R = |ω FEM − ω present |/ω FEM and spinning speed are ignored. It is shown that the difference between the FEM results and the present results is very slight, and the maximum error is less than 0.83%, which validates the validity of the present solution.

Natural Frequency and Critical Spinning Angular Speed
In this subsection, natural frequency and critical spinning speed of composite laminated spinning truncated conical shells in a hygrothermal environment are studied. Figure 2 displays the change in dimensionless natural frequencies against spinning angular speed Ω* for modes (m, n) = (1, 1-4). One can find a bifurcation of the frequencies when the speed Ω* exists, where ω 1 * and ω 2 * represent the frequencies of forward and backward traveling waves, respectively. For mode (1,1), ω 1 * descends monotonically with Ω* until the one intersects the abscissa, and then ascends again. The value of the corresponding abscissa is critical spinning angular speed Ω c *, which means a potential divergent instability. While ω 2 * ascends with Ω* before the instability occurs, after that ω 2 * descends until it coincides with ω 1 *. This phenomenon is attributed to the spin-induced Coriolis and centrifugal forces, which both affect the stiffness of the conical shells. In a physical sense, the Coriolis force induces the asymmetric influence on natural frequencies of forward and backward traveling waves, which weakens the frequency of the forward traveling wave and enlarges the frequency of the backward traveling wave. However, centrifugal force stiffens the shells and enhances the frequencies of both traveling waves symmetrically. In the lower spinning angular speed region (before critical spinning angular speed), the Coriolis force has the dominant contribution in the frequencies and leads to different trends of forward and backward traveling waves. After that, centrifugal force has the more significant effect on the frequencies in the higher angular speed region, which causes the increase in frequencies of both traveling waves.
inated spinning truncated conical shells in a hygrothermal environment are studied. Figure 2 displays the change in dimensionless natural frequencies against spinning angular speed Ω* for modes (m, n) = (1, 1-4). One can find a bifurcation of the frequencies when the speed Ω* exists, where ω1* and ω2* represent the frequencies of forward and backward traveling waves, respectively. For mode (1, 1), ω1* descends monotonically with Ω* until the one intersects the abscissa, and then ascends again. The value of the corresponding abscissa is critical spinning angular speed Ωc*, which means a potential divergent instability. While ω2* ascends with Ω* before the instability occurs, after that ω2* descends until it coincides with ω1*. This phenomenon is attributed to the spin-induced Coriolis and centrifugal forces, which both affect the stiffness of the conical shells. In a physical sense, the Coriolis force induces the asymmetric influence on natural frequencies of forward and backward traveling waves, which weakens the frequency of the forward traveling wave and enlarges the frequency of the backward traveling wave. However, centrifugal force stiffens the shells and enhances the frequencies of both traveling waves symmetrically. In the lower spinning angular speed region (before critical spinning angular speed), the Coriolis force has the dominant contribution in the frequencies and leads to different trends of forward and backward traveling waves. After that, centrifugal force has the more significant effect on the frequencies in the higher angular speed region, which causes the increase in frequencies of both traveling waves. For modes (1,(2)(3)(4), it is found that with the increase of Ω*, ω1* initially ascends and then descends until critical spinning speed Ωc* appears. It can be inferred that the initial hoop tension generated by the spinning centrifugal effect results in the difference, as shown in Equation (13). Figure 3 plots the variation of the frequencies against spinning speed for modes (1,(1)(2)(3)(4), where initial hoop tension is ignored. It is clearly seen that for modes (1,(2)(3)(4), the variation trend of ω1* against Ω* becomes similar with the mode ones (1,1). It is also found that the values of critical spinning speed are reduced obviously compared with the case in Figure 2b-d, considering initial hoop tension. However, the critical speed of mode (1, 1) only has a slight decrease compared with that in the presence of initial hoop tension. From Figures 2 and 3, one can obtain the conclusion that the initial hoop tension is the predominant factor resulting in the increase of critical spinning speed for higher modes. For modes (1,(2)(3)(4), it is found that with the increase of Ω*, ω 1 * initially ascends and then descends until critical spinning speed Ω c * appears. It can be inferred that the initial hoop tension generated by the spinning centrifugal effect results in the difference, as shown in Equation (13). Figure 3 plots the variation of the frequencies against spinning speed for modes (1,(1)(2)(3)(4), where initial hoop tension is ignored. It is clearly seen that for modes (1,(2)(3)(4), the variation trend of ω 1 * against Ω* becomes similar with the mode ones (1, 1). It is also found that the values of critical spinning speed are reduced obviously compared with the case in Figure 2b-d, considering initial hoop tension. However, the critical speed of mode (1, 1) only has a slight decrease compared with that in the presence of initial hoop tension. From Figures 2 and 3, one can obtain the conclusion that the initial hoop tension is the predominant factor resulting in the increase of critical spinning speed for higher modes.
(1, 1). It is also found that the values of critical spinning speed are reduced obviously compared with the case in Figure 2b-d, considering initial hoop tension. However, the critical speed of mode (1, 1) only has a slight decrease compared with that in the presence of initial hoop tension. From Figures 2 and 3, one can obtain the conclusion that the initial hoop tension is the predominant factor resulting in the increase of critical spinning speed for higher modes.  In Figure 4, the variation of dimensionless natural frequencies against number of circumferential waves for various spinning angular speeds is displayed. It is seen that the frequencies of forward and backward traveling waves firstly descend and then ascend with number of circumferential waves n, which illustrates that the frequency of mode (1, 1) is not the smallest among all vibration modes of the conical shells. It is also shown that the frequencies of both traveling waves with higher angular speed are more affected by circumferential wave number than those with lower speed.

Hygrothermal Analysis
In this section, hygrothermal analysis on vibration characteristics of the conical shells is carried out. Figure 5 presents the effects of temperature variation and moisture concentration on natural frequencies of forward and backward traveling waves, respec- In Figure 4, the variation of dimensionless natural frequencies against number of circumferential waves for various spinning angular speeds is displayed. It is seen that the frequencies of forward and backward traveling waves firstly descend and then ascend with number of circumferential waves n, which illustrates that the frequency of mode (1, 1) is not the smallest among all vibration modes of the conical shells. It is also shown that the frequencies of both traveling waves with higher angular speed are more affected by circumferential wave number than those with lower speed.

Hygrothermal Analysis
In this section, hygrothermal analysis on vibration characteristics of the conical shells is carried out. Figure 5 presents the effects of temperature variation and moisture concentration on natural frequencies of forward and backward traveling waves, respectively. It is manifest in Figure 5a that the frequencies of both traveling waves all descend with temperature variation ∆T, which is because the increase in temperature would weaken the stiffness of the shell. As a general reality, temperature variation can induce thermal expansion deformation and changes in material properties of the composite. In Figure 6, the effects of thermal expansion deformation and material property variation on the frequencies of both traveling waves are plotted. It can be observed that when the same temperature variation happens, the frequencies decrease with ∆T only considering thermal expansion deformation (as shown in Figure 6a), while the ones are little affected by ∆T considering solely material property variation (as shown in Figure 6b). It can be concluded that the thermal expansion deformation is nonnegligible in free vibration analysis of the spinning composite laminated truncated conical shells. circumferential waves for various spinning angular speeds is displayed. It is seen that the frequencies of forward and backward traveling waves firstly descend and then ascend with number of circumferential waves n, which illustrates that the frequency of mode (1, 1) is not the smallest among all vibration modes of the conical shells. It is also shown that the frequencies of both traveling waves with higher angular speed are more affected by circumferential wave number than those with lower speed.

Hygrothermal Analysis
In this section, hygrothermal analysis on vibration characteristics of the conical shells is carried out. Figure 5 presents the effects of temperature variation and moisture concentration on natural frequencies of forward and backward traveling waves, respectively. It is manifest in Figure 5a that the frequencies of both traveling waves all descend with temperature variation ΔT, which is because the increase in temperature would weaken the stiffness of the shell. As a general reality, temperature variation can induce thermal expansion deformation and changes in material properties of the composite. In Figure 6, the effects of thermal expansion deformation and material property variation on the frequencies of both traveling waves are plotted. It can be observed that when the same temperature variation happens, the frequencies decrease with ΔT only considering thermal expansion deformation (as shown in Figure 6a), while the ones are little affected by ΔT considering solely material property variation (as shown in Figure 6b). It can be From Figure 5b, it is indicated that the frequencies of both traveling waves decrease slightly with moisture concentration ΔC, which is attributed to hygroscopic expansion deformation and the reduction of elastic modulus E2. Figure 7 displays the variation in dimensionless natural frequencies against ΔC. The variation is similar with the variation trend of E2 in Table 1 against ΔC, which reveals that the elastic modulus E2 is the dominant factor in the frequencies of the shell, rather than the hygroscopic expansion deformation. Figure 8, which plots the effects of temperature variation and moisture concentration on critical spinning angular speed, respectively, demonstrates that temperature and moisture both impair the critical spinning speed of the shells, but temperature has more effect on the reduction of critical speed than moisture. Results in Figures 5-8 reveal From Figure 5b, it is indicated that the frequencies of both traveling waves decrease slightly with moisture concentration ΔC, which is attributed to hygroscopic expansion deformation and the reduction of elastic modulus E2. Figure 7 displays the variation in dimensionless natural frequencies against ΔC. The variation is similar with the variation trend of E2 in Table 1 against ΔC, which reveals that the elastic modulus E2 is the dominant factor in the frequencies of the shell, rather than the hygroscopic expansion deformation. Figure 8, which plots the effects of temperature variation and moisture concentration on critical spinning angular speed, respectively, demonstrates that temperature and moisture both impair the critical spinning speed of the shells, but temperature has more effect on the reduction of critical speed than moisture. Results in Figures 5-8 reveal From Figure 5b, it is indicated that the frequencies of both traveling waves decrease slightly with moisture concentration ∆C, which is attributed to hygroscopic expansion deformation and the reduction of elastic modulus E 2 . Figure 7 displays the variation in dimensionless natural frequencies against ∆C. The variation is similar with the variation trend of E 2 in Table 1 against ∆C, which reveals that the elastic modulus E 2 is the dominant factor in the frequencies of the shell, rather than the hygroscopic expansion deformation. Figure 8, which plots the effects of temperature variation and moisture concentration on critical spinning angular speed, respectively, demonstrates that temperature and moisture both impair the critical spinning speed of the shells, but temperature has more effect on the reduction of critical speed than moisture. Results in Figures 5-8 reveal that the influence of moisture concentration on the frequencies and critical spinning speed of the shell are less prominent than that of temperature.

Effects of Design Parameters on Vibration Characteristics
In this section, the effects of design parameters on vibration characteristics are presented, including semivertex angle and fiber orientation angle. Figure 9 plots the influence of semivertex angle on dimensionless natural frequencies and critical spinning speed. One can see that the frequencies of forward and backward traveling waves vary obviously with semivertex angle α, and the critical spinning speed is reduced with α. It is a fact that the influence of the semivertex angle on the frequencies of both traveling waves depends on the vibration modes. In Figure 10, the variation in frequencies against semivertex angle for various circumferential wave numbers is presented. It is seen that the frequencies of both traveling waves initially increase and then decrease with α for all vibration modes, but the frequencies are more sensitive to the lower-order modes than the higher ones.

Effects of Design Parameters on Vibration Characteristics
In this section, the effects of design parameters on vibration characteristics are presented, including semivertex angle and fiber orientation angle. Figure 9 plots the influence of semivertex angle on dimensionless natural frequencies and critical spinning speed. One can see that the frequencies of forward and backward traveling waves vary obviously with semivertex angle α, and the critical spinning speed is reduced with α. It is a fact that the influence of the semivertex angle on the frequencies of both traveling waves depends on the vibration modes. In Figure 10, the variation in frequencies against semivertex angle for various circumferential wave numbers is presented. It is seen that the frequencies of both traveling waves initially increase and then decrease with α for all vibration modes, but the frequencies are more sensitive to the lower-order modes than the higher ones.

Effects of Design Parameters on Vibration Characteristics
In this section, the effects of design parameters on vibration characteristics are presented, including semivertex angle and fiber orientation angle. Figure 9 plots the influence of semivertex angle on dimensionless natural frequencies and critical spinning speed. One can see that the frequencies of forward and backward traveling waves vary obviously with semivertex angle α, and the critical spinning speed is reduced with α. It is a fact that the influence of the semivertex angle on the frequencies of both traveling waves depends on the vibration modes. In Figure 10, the variation in frequencies against semivertex angle for various circumferential wave numbers is presented. It is seen that the frequencies of both traveling waves initially increase and then decrease with α for all vibration modes, but the frequencies are more sensitive to the lower-order modes than the higher ones. In Figure 11, the variation in dimensionless natural frequencies of forward and backward traveling waves against fiber orientation angle θf is presented, where θf represents (θf/−θf)s layup of laminate. It is seen that both frequencies firstly ascend and then descend with θf for all vibration modes, but the variation in frequencies of lower-order modes is more significant than the ones of higher modes. In a physical sense, lamination schemes of composites would affect the stiffness of the conical shells, which leads to the variation of the frequencies. In Figures 9-11, we can see that the dynamic behavior of spinning composite truncated conical shells can be distinctly controlled by employing suitable values of design parameters. In Figure 11, the variation in dimensionless natural frequencies of forward and backward traveling waves against fiber orientation angle θf is presented, where θf represents (θf/−θf)s layup of laminate. It is seen that both frequencies firstly ascend and then descend with θf for all vibration modes, but the variation in frequencies of lower-order modes is more significant than the ones of higher modes. In a physical sense, lamination schemes of composites would affect the stiffness of the conical shells, which leads to the variation of the frequencies. In Figures 9-11, we can see that the dynamic behavior of spinning composite truncated conical shells can be distinctly controlled by employing suitable values of design parameters. In Figure 11, the variation in dimensionless natural frequencies of forward and backward traveling waves against fiber orientation angle θ f is presented, where θ f represents (θ f / − θ f ) s layup of laminate. It is seen that both frequencies firstly ascend and then descend with θ f for all vibration modes, but the variation in frequencies of lower-order modes is more significant than the ones of higher modes. In a physical sense, lamination schemes of composites would affect the stiffness of the conical shells, which leads to the variation of the frequencies. In Figures 9-11, we can see that the dynamic behavior of spinning composite truncated conical shells can be distinctly controlled by employing suitable values of design parameters.

Conclusions
Free vibration characteristics of a spinning composite laminated truncated conical shell subjected to hygrothermal environment were examined in the present work. The governing equations of the free vibration of the conical shell were established based on Love's thin-shell theory and Hamilton's principle. The Galerkin approach was used to derive the free vibration characteristics of the conical shell. The combined effects of spinning angular speed, temperature, moisture concentration, semivertex angle, and fiber orientation angle on natural frequencies and critical spinning speeds were studied systematically. The conclusions are as follows: 1. The Coriolis force is the dominant factor resulting in the asymmetric variation of natural frequencies of forward and backward traveling waves. The centrifugal force would stiffen the conical shell and enhance the frequencies of both traveling waves symmetrically, and initial hoop tension plays a major role in the increase of critical spinning angular speed. 2. Temperature and moisture concentration both weaken natural frequencies and critical spinning speeds of the shell, and the influence of temperature on the frequencies and the speed is more prominent than that of moisture concentration. In addition, thermal expansion deformation is nonnegligible in free vibration analysis of spinning composite laminated truncated conical shells. 3. Natural frequencies of forward and backward traveling waves all initially increase and then decrease with semivertex angle and fiber orientation angle, which indicates that vibration behavior of the conical shell can be distinctly controlled by using suitable values of design parameters. Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available from the corresponding author upon reasonable request.

Conclusions
Free vibration characteristics of a spinning composite laminated truncated conical shell subjected to hygrothermal environment were examined in the present work. The governing equations of the free vibration of the conical shell were established based on Love's thin-shell theory and Hamilton's principle. The Galerkin approach was used to derive the free vibration characteristics of the conical shell. The combined effects of spinning angular speed, temperature, moisture concentration, semivertex angle, and fiber orientation angle on natural frequencies and critical spinning speeds were studied systematically. The conclusions are as follows:

1.
The Coriolis force is the dominant factor resulting in the asymmetric variation of natural frequencies of forward and backward traveling waves. The centrifugal force would stiffen the conical shell and enhance the frequencies of both traveling waves symmetrically, and initial hoop tension plays a major role in the increase of critical spinning angular speed.

2.
Temperature and moisture concentration both weaken natural frequencies and critical spinning speeds of the shell, and the influence of temperature on the frequencies and the speed is more prominent than that of moisture concentration. In addition, thermal expansion deformation is nonnegligible in free vibration analysis of spinning composite laminated truncated conical shells.

3.
Natural frequencies of forward and backward traveling waves all initially increase and then decrease with semivertex angle and fiber orientation angle, which indicates that vibration behavior of the conical shell can be distinctly controlled by using suitable values of design parameters.
Author Contributions: Conceptualization, X.L. and Z.Z.; methodology, X.L.; validation, X.L. and Z.Z.; data curation, X.L., X.Z. and Z.Z.; writing-original draft preparation, X.L.; writing-review and editing, X.L. and Z.Z.; supervision, Z.Z.; funding acquisition, X.L. All authors have read and agreed to the published version of the manuscript. Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
On behalf of all authors, the corresponding author states that there are no conflicts of interest.