Abstract

This paper carries out forced vibration analysis of graphene nanoplatelet-reinforced composite laminated shells in thermal environments by employing the finite element method (FEM). Material properties including elastic modulus, specific gravity, and Poisson’s ratio are determined according to the Halpin–Tsai model. The first-order shear deformation theory (FSDT), which is based on the 8-node isoparametric element to establish the oscillation equation of shell structure, is employed in this work. We then code the computing program in the MATLAB application and examine the verification of convergence rate and reliability of the program by comparing the data of present work with those of other exact solutions. The effects of both geometric parameters and mechanical properties of materials on the forced vibration of the structure are investigated.

1. Introduction

Nowadays, due to the development of science and technology, a variety of new materials have been applied widely to engineering applications such as composite materials, functionally graded materials (FGM), and piezoelectric materials (PZT). To enhance the strengths of structures, there are many common ways like adding stiffeners, using folded structures, or reinforcing the structures with smart materials. In these techniques, the carbon nanotube (CNT)-reinforced composite shell is one of the most modern structures [1]. However, in comparison with nanographene or graphene nanoplatelet (GPL) structures, carbon nanotube (CNT) structures are not only inflexible but also have a high cost [2]. In recent years, GPL reinforcement for composite materials (GPLRC) has been manufactured by Stankovich et al. [35]. The investigations found that only a small amount of GPL in the material could significantly improve contemporaneously its physical characteristics [6, 7]. In comparison with conventional composite materials containing more than 60% carbon fiber by the volume, GPLRC contains only a low graphene nanoplatelet (GPL) ratio (0.5%–20% weight) [810]. Therefore, GPLRC is concerned by many scientists all over the world. In civil engineering, cylindrical shells (CYL) and spherical shell (SPH) reinforced with GPL are considered with many different approaches.

In [11], nonlocal continuum theory and molecular dynamics simulations were used to discover the free vibration of single and double-layer graphene structures. More works dealt with nonlocal methods. Ghannadpour et al. [12] performed bending, buckling, and vibration analysis of Euler beams, and the same mechanical response analysis of their work [13] for Timoshenko beam was also carried out. Chandra et al. [14] investigated natural frequencies and oscillation mode shapes of GRP composite plates with several types of boundary supports using numerical approaches. The linear buckling and nonlinear buckling of multilayer beam structures reinforced by GPL with the Halpin–Tsai model were studied by Yang and his co-workers [15], where effects of material characteristics of each component on the vibration of the structure are examined. Feng and his colleagues [16] developed a relatively complete investigation of nonlinear static bending of multilayer beams reinforced with graphene. Then, he [17] extended his study to carry out the nonlinear free vibration problems. The material characteristics calculated by the Halpin–Tsai model and the von Kármán nonlinear deformation-displacement were used in both studies. Kitipornchai and colleagues [18] studied free vibration and stabilization of FGM beams with reinforced holes in the graphene. Lin et al. [19] carried out the investigation of mechanical characteristics of grapheme-reinforced components. In [20], Zhao and his co-workers performed the static bending and free vibration analysis of FGM plates reinforced with nanographene by using the newly proposed Halpin–Tsai model. Mirzaei and Kiani [21] examined the thermal buckling of composite structures (plate structures) reinforced with graphene based on FSDT where the volume of graphene in each component layer was different from one another.

Shen and his co-workers [22] carried out the nonlinear bending, buckling, and postbuckling analysis of composite beams reinforcing graphene on an elastic foundation or in a thermal environment. Shen et al. [23, 24] introduced the results of nonlinear bending and vibration analysis of FGM plate-reinforced graphene under the effect of thermal load. Then, they [25, 26] developed the linear and nonlinear thermostats of composite structures reinforcing graphene resting on an elastic foundation employing high-order shear deformation theory (HSDT). In [27], the author used both the Von Kármán-type nonlinearity and the third-order shear deformation theory (TSDT) to carry out the vibration analysis of graphene-reinforced cylinders resting on the elastic foundation (in the thermal environment) [28]. For buckling problems, Wang and colleagues [29] introduced the buckling analysis of cylindrical composite shells using FEM based on the Halpin–Tsai model and mixed rule in order to determine the mechanical characteristics of each component layer. Then, the buckling problem of composite CYL-reinforced graphene with holes and FGM CYL reinforcing graphene was presented in [30, 31]. They [32] also studied buckling of dynamic multilayer nanocomposite beams under the thermal load. Wu et al. [33] carried out the parameter buckling of the grapheme-reinforced multilayer nanocomposites under both thermal and mechanical loads. Kiani [34] examined the large deformation of the composite plate with reinforced graphene in which the masses of graphene in each component were divergent. Rout et al. [35] studied on nonlinear vibration of composite shell types (CYL, SPH, HYP, and so on) reinforced with graphene using higher-order shear deformation theory (HSDT). Recently, Duan et al. [36] introduced a novel formulation by employing the discrete singular convolution (DSC) for free vibration analysis of circular thin plates with uniform and stepped thickness. Civalek [37] investigated the free oscillation of CNT and FGM structures. Nguyen et al. [38, 39] computed dynamic responses of the composite shell with shear connectors using FEM.

In the above studies, to our knowledge, most plate and shell structures are reinforced with graphenes, where a few of which dealt with graphene-reinforced composite shells. In addition, the forced vibration problems of graphene-reinforced composite shells with the different mass fractions of graphene of each component in thermal environments are still limitations. The main goals of this work are to analyze the forced vibration of the shell structures (Figure 1) with reinforced graphene performing in a thermal environment using FEM. Herein, we consider four models with different distributions of graphene including UD-type, X-type, V-type, and O-type to study in this paper. For UD structure, the dispensation rule of graphene shells is homogeneous in each component layer. On the other hand, for X structure, the graphene plates are found to be minimal in the middle surface and gradually increase to the top and bottom layers. For V-type structure, the volume of each part of graphene is detected as maximum in the top layer and minimum in the bottom layer. For O structure, the volume fraction of graphene is maximal in middle surface of the structure and minimum in the top layer and the bottom layer. The temperature-dependent material characteristics of the multilayer structure are evaluated by using the Halpin–Tsai extension method. In this paper, the authors investigated the influences of both geometrical and physical properties on the forced vibration with two cases of shells (cylindrical shell (CYL) and spherical shell (SPH)) in the thermal environment.

The organization of this paper is divided into 4 main sections. Section 2 presents briefly the equation of motion of the shell element by using finite element method (FEM). Numerical results and discussions are shown in Section 3. Some highlighted results are concluded in Section 4.

2. Finite Element Formulations

2.1. Equation of Motion of the Shell Element

Consider a curvilinear shell model reinforced graphene nanoplatelets with the following dimensions: length a, width b, and thickness h, as shown in Figure 1.

denote the radii of the shells:(i)CYL: (Case 1)(ii)SPH: (Case 2)

The displacement field is defined through FSDT as follows:where , and are the displacement components at any point of the shell; , and are the displacement components at the midsurface of the shell; and and are the rotation angles of the cross section of the structure around the y-axis and x-axis.

The stress-strain relationship is expressed as follows:

Equation (2) can be shortened in the matrix form as follows:in which

The stress and strain fields of layer k has a relationship as follows:where are the bending rigidity matrix and shear rigidity matrix of layer k:withwhere is the fiber angle of the k layer; are coefficients of thermal expansion of the k layer in the x and y directions. Herein, we employ the extended Halpin–Tsai method to estimate the effective elastic modulus as well as shear modulus of the structures which are reinforced by using graphene nanoplatelets. In the present considerations, we choose two options of reinforcement of graphene nanoplatelets, which are zigzag (θ = 0°) and armchair (θ = 90°). The operative material characteristics of composite structures reinforced with graphene are presented as [27]where (i = 1–3) represent graphene proficiency parameters calculated by comparing the results of the work in [27]. Em and Gm are, respectively, elastic modulus and shear modulus. The dimensions of the graphene sheet include length , width , and effective thickness . The volume component of graphene is VG. The coefficients are defined aswhere , , and are Young’s moduli and shear modulus of the graphene structures. The density of the mass and Poisson’s ratio of the composite structures are defined aswhere(i)Vm is the matrix volume fraction(ii) and are, respectively, Poisson’s ratio of the grapheme and matrix(iii) are are, respectively, the mass density of the graphene and the matrix

The parameters of the thermal development along the longitudinal and transverse directions of composite structures reinforced with graphene can be defined as follows:

In this work, we employ the 8-node isoparametric element, where each node consists of five degrees of freedom (DOF); the degree of freedom of node i is , and the total DOFs of the shell element is expressed as follows:in which Ni can be found in Ref. [39].

By substituting in the expression for verifying displacement of element, we havewhere are defined as follows:where , , and can be found in Appendix A.

To obtain the dynamic equation, we employ the weak form for each element, and then we obtain

By substituting equations (1) and (13) into equation (15), we have the dynamic equation of the shell element as follows:withwherewhere can be found in Appendix B.in which

For the case of taking into account the structural damping ratio, we obtain the forced oscillation connection of the structure element as follows:in which and are Rayleigh drag coefficients, which are calculated as in the literature [40].

2.2. The Differential Equation of Vibration

From the differential equation of forced vibration of the shell element (21), we obtain the differential equation of forced vibration of entire composite shell structure as follows:in which , , , , and are, respectively, the global mass matrix, the global structural damping matrix, the global stiffness matrix, the global load vector, and the global thermal load vector. These matrices and vectors are assembled from the element matrices and vectors, correspondingly. These are linear differential equations, where the right-hand side of the equation varies over time. To analyze these equations, we employ the Newmark-beta approach [40]. The program is coded in the MATLAB (MathWorks, Natick, MA, USA) environment with the algorithm flowchart of Newmark as shown in Figure 2.

When and , from equation (24), natural frequencies can be obtained by figuring out the roots of the following equation:

By assuming that , we obtainwhere is the natural frequency.

Now, we have the flowchart of the Newmark-beta method with symbols “n” is the total number of integral steps and “j” is the number of integral step.

Step 1. Determine the initial conditionsFrom the initial conditions, we obtain

Step 2. By approximating by , we havewherein which are defined by the assumption that the acceleration varies in each calculating step; the authors select the linear law for the varying acceleration:The condition to stabilize the roots is expressed as follows:

Step 3. Compute the efficiency stiffness matrix and the efficiency nodal force vector:

Step 4. Determine the nodal displacement vector :and then repeat the loop until the time runs out.

3. Numerical Examples and Discussion

Based on the finite element equations, the authors then code the calculation program in the MATLAB environment to examine the influence of some geometric parameters and mechanical properties on the forced vibration of the shell. In this study, we consider four different cases of graphene reinforcement: UD-type, X-type, O-type, and V-type with ten composite layers, which depend on the mass fraction of graphene. For the whole analysis, the volume component of graphene in the all considered cases of ten composite layers is expressed as follows [35]:(i)Case 1: UD-type: [0.07/0.07/0.07/0.07/0.07]s(ii)Case 2: X-type: [0.11/0.09/0.07/0.05/0.03]s(iii)Case 3: O-type: [0.03/0.05/0.07/0.09/0.11]s(iv)Case 4: V-type: [(0.11)2/(0.09)2/(0.07)2/(0.05)2/(0.03)2]

Herein, we choose polymethyl methacrylate (PMMA) as the matrix, which has temperature-dependent material properties as Em = (3.52–0.0034 T) GPa, , , and , in which T0 is the temperature inside the room, and it is equal to 300 K. Temperature-dependent material characteristics of monolayer graphene with  = 14.76 nm,  = 14.77 nm,  = 0.188 nm, and are shown in Table 1. PMMA (C5O2H8)n is a transparent thermoplastic, which is often used in the plate and shell structures due to its advantages as a light weight, shatter-resistant alternative to glass. The same material can be used as inks, coatings, a casting resin, and so on. Although PMMA is not considered as a type of familiar silica-based glass, the substance, similar to many thermoplastics, is often technically categorized as a type of glass (in that it is called as a no-crystalline vitreous substance); therefore, it is named as acrylic glass. PMMA is an economical alternative to polycarbonate (PC) with high tensile and flexural strength, transparency, polishability, and chemical and heat resistance; especially, it can be protected from UV tolerance. Nonmodified PMMA behaves in a brittle manner under load, typically when suffering impact load, and is more prone to scratching than other common inorganic glass, but modified PMMA sometimes could gain high ability in this property.

In Table 2, we illustrate the efficiency parameters which are related to the volume component of graphene.

3.1. Verification Examples

Example 1. Let us consider a (0°/90°/0°/90°/0°)s composite shell with geometrical parameters: a = b, thickness h = a/10, radii Rx = Ry = 50 h (for elliptic paraboloid-EPS), Rx = 50 h, Ry = ∞ (for cylindrical-CLY). Material properties of the shell are presented in Tables 1 and 2. The comparative data of the first nondimensional free vibration with those of Rout et al. [35] are shown in Table 3 ( and are mass density and elastic modulus of shell at the room temperature T = T0 = 300 K, respectively).
In ref. [35], Rout et al. employed the finite element method (FEM) based on high-order shear deformation theory (HSDT) with the 8-node isoparametric element, where each node consists of 7 degrees of freedom (DOF). For shell structures with medium thicknesses (h/a = 1/100 ÷ 1/10) and the radius (R/a > 50), using the 8-node isoparametric element with 5 DOF per one element to analyze free vibration of the shell is very convenient. From Table 3, we can see clearly that the computed results between this work and the analytical method [35] are in good agreement; therefore, it demonstrates that our proposed theory and program are verified for free vibration problems and the convergence rate is guaranteed with the mesh 8 × 8. However, in order to improve accuracy, we used 10 × 10 mesh.

Example 2. Consider a fully clamped square plate with parameters that can be found in [41]: a = b = 1 m and h/a = 10. Material characteristics are Young’s modulus E = 30 GPa, Poisson’s ratio  = 0.3, and the density  = 2800 kg/m3. The structure is subjected to distribution blast load p0 = 104 Pa. The nondimensional displacement is calculated by the formula . The calculation program is used with the mesh 10 × 10 (let the radii of the shell be very large , the elastic modulus be E1 = E2 = E, and the fiber angle of the k layer be  = 0°). The deflection of the center point of the structure and those of [41] are shown in Figure 3 (integral time is 5 ms, and acting time of load is 2 ms).
We can see from Figure 3 that the results are similar to each other (both shape and value). This proves that our program is verified.

3.2. Effects of Some Parameters on Forced Vibration of the Shell

Next, we study the effects of some parameters on forced vibration of the shell structure; we consider a fully clamped composite shell consisting of ten composite layers with the following geometrical parameters: the length a = 1 m, the width b, the thickness h, and radii of the shell Rx = Ry = R; material parameters are shown in Tables 1 and 2. All cases of graphene reinforcement are UD-type, X-type, O-type, and V-type as mentioned. The external uniform load acts on the upper surface of the shell in the z-direction.

Nondimensional displacement and velocity of the center point over time are given as follows:where and are the deflection and velocity of the centroid of the shell, respectively.

The nondimensional frequencies are calculated as (where is value of at T = 300 K). For fully clamped shell, the first five nondimensional frequencies are presented in Table 4. From numerical results, it can be seen that with the same boundary condition and the temperature in 4 cases of graphene reinforcement, X-type gives the largest nondimensional frequencies, followed by UD-type, V-type, and O-type, respectively. This demonstrates that the composite shell reinforced with graphene X-type is the strongest one and the smallest one is O-type. For all distributions of graphene, the temperature decreases the nondimensional fundamental frequencies. Therefore, “stiffness” of the shell would be reduced.

First four vibration mode shapes of the shell structure are presented in Figure 4. For the case of clamping all edges of the shell, the second and the third vibration mode shapes are similar to each other (the second fundamental frequency is approximated the third fundamental frequency). This phenomenon is suitable to the actual symmetrical shell structures under the same supporting conditions.

3.3. Influence of Four Different Types of Graphene-Reinforced Structures

In this section, we investigate the effect of four different types of graphene reinforcement. Let us consider a fully clamped (CCCC) composite shell with geometrical parameters a = b = 1 m, h = a/25, and R/a = 10. Nondimensional deflection, velocity, and the stress (when the deflection obtains the maximum value) of the centroid of the shell at the room temperature T = T0 = 300 K are presented in Figure 5. In the case of T = 400 K, nondimensional results (including velocity, stress and deflection) of the center point of the shell are presented in Figure 6 and the maximum value is introduced in Table 5.

From Figures 5 and 6 and Table 5, we can see that in four cases of graphene reinforcement, X-type gives the smallest displacement, velocity, and stress responses of the centroid of shell structure, followed by the UD-type, V-type, and O-type, respectively. This demonstrates that with the same geometrical parameters, the composite shell reinforced with graphene X-type is the strongest one and the smallest one is O-type. Due to the amount of graphene nanoplatelets added in the original material, the structure will become stiffer. The distribution law of graphene nanoplatelets also affects strongly on the stiffness of structures. For example, X-type has more nanographene particles focused on the surface of the shell while O-type has more nanographene particles focused on the neutral plane of the shell; therefore, the stiffness of X-type is higher than that of O-type. Similarly, in the cases of nanographene particles distributing uniformly in layers of the structure (UD-type) and reducing smoothly from the top surface to the bottom surface, the stiffeners are smaller than those of X-type. It can be confirmed that nanographene particles are distributed more on the surface of the shell and less gradually into the core (X-type); we can obtain the structure, where the surface layer is stiffer than the core layer. This characteristic will help the structure can absorb energy better. One more thing is that in the thermal environment, all types of reinforced shell structures give the smaller displacement, velocity, and stress responses than those of the structures in the room temperature. This phenomenon shows that the temperature reduces the stiffness of the structure, in which the thermal stress appears; after the acting process, the displacement response tends to move away from position 0.

3.3.1. Influence of a/h Ratio

Next, to investigate the influence of a/h ratio, we consider the shell structure with geometrical parameters b = a (while a is fixed), the value of a/h ratio gets 10, 25, 50 and 75, R/a = 10. The nondimensional deflection, velocity, and the stress (when the deflection reaches the maximum value) of the centroid of the shell are presented in Figures 7 and 8, and the maximum value is listed in Table 6.

We can see in Figures 7 and 8 and Table 6, for all cases of graphene-reinforced shells (UD-type, X-type, O-type, and V-type), when the a/h ratio increases, nondimensional deflection, velocity, and stress (change over time) of the center point also increase rapidly. For the cases of nanographene distribution, X-type gives the smallest displacement, velocity, and stress responses, followed by UD-type and V-type; the largest one is O-type. From Figures 8(c) and 8(d), we can see that the stress (by the thickness direction) of the center point of the shell (when the displacement obtains the maximum value) is suitable to the nanographene distribution. It is appropriate in this situation because with the same geometric dimensions, materials, boundary conditions, and thinner shells have smaller stiffness. And the temperature is also the reason of increasing displacement, velocity, and stress of the centroid point.

3.3.2. Influence of a/b Ratio

Next, we examine the influence of the length-to-width ratio on the nondimensional deflection and velocity of the center point of the shell with a = 1 m; a/b gets the values of 0.5, 1, 1.5, and 2. Geometrical parameters are h = a/25 and R/a = 10. The numerical results of the nondimensional deflection, velocity, and stress (when the deflection reaches the maximum value) of the centroid of the shell are shown in Figures 9 and 10, and the maximum values are listed in Table 7.

Now, we can see in Figures 9 and 10 and Table 7 that when increasing a/b ratio, the nondimensional deflection and velocity of the center point overtime decrease for all cases of graphene-reinforced shells (UD-type, X-type, O-type, and V-type). This demonstrates that the stiffness of the shell gets larger, especially when a/b ratio equals to 4. The stress (by the thickness direction) of the center point of the shell is suitable to the nanographene distribution (see Figures 10(c) and 10(d)). This can be understood obviously that when the shape of structure gets smaller, while boundary condition and other parameters are not changed, the structure will become stronger.

4. Conclusions

In this paper, the numerical results of mechanical responses of laminated composite shells reinforced with graphene nanoplatelets are explored. The authors use FEM based on FSDT that has the following advantages:(i)Simple formulations for a theoretical representation(ii)The computed results of free vibration and forced vibration obtained by this approach compared to other solutions show a good agreement

From new computed results, several conclusions may be drawn as follows:(i)With the same geometric parameters, the X-type has the largest hardness and the smallest one is the O-type structure. In the temperature environment, the shells with all reinforcement types are greater in the displacement response, velocity, and stress than those of the shells in the room temperature. This demonstrates that heat reduces the stiffness of the shell, and thermal stresses appear in the structure. Thus, due to the effect of the external load, the displacement response tends to move away from position 0.(ii)When increasing a/b ratio, the stiffness of the shell gets larger, especially when a/b ratio equals 4. This can be understood obviously that when the shape of structure gets smaller, while boundary conditions and other parameters are not changed, the structure will become stronger.(iii)We also understand that in all cases, the maximum displacement of the centroid of the structures performing in the temperature environment T = 400 K is approximately 10–20% larger than that in case of T = 300 K. Hence, we would like to express that the temperature has a small influence on the deflections of the structure under the impact loads in the thermal environment. These remarkable results are a reliable suggestion for engineers when using these laminated composite shells reinforced with graphene nanoplatelets in the thermal environment.

Finally, based on achieved numerical results, the proposed program is able to analyze the static bending, dynamic response, nonlinear problems, etc., with complicated structures, which is not convenient to solve by analytical and other methods.

Appendix

A. The Matrices in Equation (14)

where .

B. The Matrices in Equation (18)

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

DVT gratefully acknowledges the support of Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant no. 107.02-2018.30.