Vibration Characteristics Analysis of Composite Laminated Annular/Circular Plate Using High-Order Shear Deformation Theory

Jin, Guodong; Ma, Lei; Zhang, Hong; Wang, Qingshan

doi:https://doi.org/10.1155/2021/5913860

Shock and Vibration

On this page

Abstract Introduction Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 5913860 | https://doi.org/10.1155/2021/5913860

Vibration Characteristics Analysis of Composite Laminated Annular/Circular Plate Using High-Order Shear Deformation Theory

Guodong Jin,¹Lei Ma,¹Hong Zhang,²and Qingshan Wang³

Academic Editor: Vasudevan Rajamohan

Received13 Apr 2021

Revised12 May 2021

Accepted21 Jul 2021

Published31 Jul 2021

Abstract

In this paper, the free vibration behaviors of composite laminated annular and circular plates under complex elastic boundary constraints are investigated. Firstly, Reddy’s high-order shear deformation theory (HSDT) and Jacobi polynomial method are effectively combined to establish the unified vibration analysis model of composite laminated annular and circular plates. Secondly, the simulation of complex elastic boundary and coupling boundary is realized by using artificial virtual spring technology. Then, the energy equation of the composite laminated plate is established by using Rayleigh–Ritz energy technology. Finally, the free vibration solution equation of the laminated plate is obtained through the Hamilton differential principle. The fast and uniform convergence of this method and the accuracy of the calculated results are verified by numerical examples and the model experimental method. On this basis, the parameterization study is conducted, and the effects of material parameters, geometric parameters, spring stiffness values, and lamination scheme on the vibration characteristics of the annular or circular plate are fully discussed, which can provide a theoretical basis for future research.

1. Introduction

The composite laminated annular or circular plates are regarded as the basic components of complex structures in various engineering and industrial application. Therefore, in-depth investigation of the free vibration characteristics of this kind of structure is quite meaningful for the vibration reduction, noise reduction, and impact resistance of the plate. Many scholars have been trying to find new methods or improve the existing methods to further study the vibration mechanical properties of the composite plate structure and have made fruitful research results. Among these research results, the commonly used theories are two-dimensional equivalent single-layer plate theory, zig-zag theory, layerwise theory, and three-dimensional elastic theory. The equivalent single-layer plate theory mainly includes classical thin plate theory (CPT), first-order shear deformation theory (FSDT), and high-order shear deformation theory (HSDT). Next, these theories are taken as the main line and briefly introduce the existing research results.

According to the Ritz method and CPT, Afsharmanesh et al. [1] solved the forced vibration problem of laminated circular plates resting on a Winkler-type foundation. Parametric research is mainly focused on the various classical edge constraints and fiber orientation. Amabili et al. [2] investigated the free vibrations of a bolted (or riveted) circular plate having elastic constraints variable by using the Rayleigh–Ritz method. Arshid and Khorshidvand [3] studied the free vibration analysis of saturated porous functionally graded (FG) circular plates based on the differential quadrature method (DQM) and CPT. Bahrami and Teimourian [4] presented a wave propagation approach to analyze free vibration of nonuniform annular and circular membranes. As the first theory to be put forward, CPT completely ignores the shear deformation of the plate, so it is only suitable for studying the low frequency vibration behavior of thin plates or shells.

Compared with CPT, the FSDT simply considers the effect of shear deformation on structural vibration by introducing the shear factor. Therefore, this theory can easily solve the vibration behaviors of the plate or shell with medium thickness. Mercan et al. [5] investigated important parametric results for annular plates and conical panels with isotropic, laminated, and FG materials by using FSDT and discrete singular convolution (DSC) method. On the basis of FSDT, Viswanathan et al. [6] expanded the Bickley-type spline function approximation to study the vibration problem of the symmetric angle-ply laminated annular and circular plate of variable thickness. Powmya and Narasimhan [7] discussed the free vibration analysis of orthotropic laminated circular and annular plates based on FSDT. They found that the transverse shear effects of orthotropic laminates are more prominent than that of isotropic laminates. Draiche et al. [8] presented an analytical model to predict the static analysis of laminated reinforced composite plates subjected to sinusoidal and uniform loads by using a simple first-order shear deformation theory. Although the FSDT further considers the effect of transverse shear deformation by introducing shear factor, this theory cannot reflect the warping deformation of the plate in the thickness direction more comprehensively.

In order to overcome the limitations of the CPT and FSDT, the HSDT has been developed, which can better describe the vibration of the composite laminated plate. Amabili et al. [9] investigated forced nonlinear vibrations around the frequency of the fundamental mode by introducing a geometrically nonlinear damping model based on third-order shear deformation theory. In addition, he and his team have carried out a series of studies on the nonlinear mechanics of composite plate and shell structures, based on the third-order shear deformation theory [10–12]. Bisadi et al. [13] employed an exact closed-form solution for freely vibrating annular thick plates under traditional boundary constraints based on Reddy’s HSDT. Asadi and Jam [14] employed Mori–Tanaka method to calculate the buckling property of annular plates reinforced by carbon nanotubes on the basis of HSDT. The effects of CNTs’ volume fractions, orientation angles, boundary conditions, and geometric ratio of the plate had been analyzed. Vinyas et al. [15] studied the free vibration behavior of annular and circular magneto-electro-elastic plates through the combination of finite element method and HSDT. Hosseini-Hashemi et al. [16] proposed an exact analytical solution to discuss the freely vibrating piezoelectric coupled circular/annular thick plate based on Reddy’s HSDT. The boundary constraints of the plate are set to the combinations of free, simply support, or fixed support. Besides, he and his team [17] presented exact closed-form solutions to investigate free vibration behavior of circular thin, moderately thick, and thick plates with classical boundary constraints based on Reddy’s HSDT. Compared with the FSDT, the results obtained by HSDT get closer to the results obtained by 3D elastic theory. Najafizadeh and Heydari [18] considered mechanical buckling of FG circular plate with uniform radial compression based on HSDT. The study found that HSDT accurately predicts the behavior of the FG circular plate, whereas the CPT and FSDT overestimate buckling loads. In addition, some scholars have proposed a simplified higher-order shear deformation theory to further study the natural characteristics of composite plate or shell structures [19–23]. Katariya et al. [24, 25] studied the vibration characteristic of composite structures based on the HSDT and experimental test method. It is not difficult to find that there are still limitations in the study of vibration characteristics of composite laminated plate structures by using HSDT. Especially, for complex elastic boundary conditions of the laminated plate, the research results are relatively few.

According to the limitations of the above research results, this paper effectively combines Reddy’s HSDT and Jacobi polynomials method to establish the unified analytical model of the laminated annular plate and circular plate. This analysis model can quickly and accurately obtain the free vibration behaviors of laminated annular and circular plates under complex elastic boundary constraints. Through the parametric study of boundary spring stiffness, the range of elastic boundary conditions is determined. On this basis, several groups of different boundary settings are given. Besides, in order to ensure the continuity when the rotation angle is 360 degrees, the stiffness value of the connection spring is clearly given. Through the parametric study, the effects of geometric parameters, material parameters, ply angle, and other important parameters on the vibration behaviors of the plate structures are effectively revealed. These research results can provide a theoretical basis for future research.

2. Theoretical Formulations

2.1. Model Description of the Laminated Annular/Circular Plate

The theoretical analysis model can be established because of the correlation between the geometry and elasticity of the annular plate and the circular plate. According to the relationships between geometry and elasticity of the annular plate and circular plate, the theoretical analysis model can be found in Figure 1. First of all, as shown in Figure 1(a), it is not difficult to find that the cylindrical coordinate system of the plate is built on the midplane. In this coordinate system, r-axis represents the radius direction, and the corresponding geometric parameters are internal radius R₀ and external radius R₁, respectively. R represents the difference between R₁ and R₀. In particular, the circular plate is set when the inner radius is zero. θ-axis represents the rotation direction of the laminated plate, which ranges from 0 to 2π. Z-axis donates the thickness direction of the plate, which ranges from 0 to h. Then, the overall displacement functions of this kind of plate correspond to r-axis, θ-axis, and z-axis, respectively.

(a)

(b)

In this paper, simulations of various complex elastic boundary constraints are achieved by applying artificial virtual spring technology. All in all, spring components are placed evenly at the boundaries of the annular or circular plate. Among them, the number and type of spring components should be determined according to the displacement function on the middle plane. In this paper, the spring components are shown in Figure 1(b). Taking the outer radius of the laminated circular plate as an example, three linear springs (, , and ) and four rotary springs (, , , and ) form a boundary spring component. When the study model is changed from the annular plate to the circular plate, not only the inner radius is set to 0 but also the stiffness of these boundary springs at this position is set to 0. Through this series of settings, the simulation of different boundaries can be easily realized.

2.2. Kinematic Relations and Stress-Strain Relations

In the displacement field, the displacements U and can be extended to the cubic functions of the z-axis. Besides, transverse deflection is assumed to be constant along the thickness direction of the plate. Therefore, the relationship between the global displacement function and the midplane displacement function can be expressed as follows [26, 27]:where are the displacements of the midplane. Then, and denote the rotation foundations of normal to midplane about the r-axis and θ-axis, respectively. Besides, , , , and can be solved by using the vanish of transverse shear stresses at the top and bottom of the plate [26, 27]. Next, a simple solution process will be given.

According to the three-dimensional elasticity theory of the annular/circular plate, the strain-displacement relations in the cylindrical coordinate system can be obtained:

Based on the simple high-order theory put forward by Reddy [26, 27], , , , and can be obtained by setting . So, their expressions can be written as

In this way, combining equations (1) and (3), the displacement function can be rewritten as follows:in which and . It is not difficult to find from equation (4) that the unknown number of displacement functions involved is the same as that of the FSDT. But the present HSDT explains the parabolic distribution of transverse shear strains along the thickness of the composite laminated annular or circular plate.

Substituting equation (4) into (2), it can be rewritten aswhere the strain‐displacement relationships on the mid‐plane are given in Table 1.

The unidirectional fiber reinforced layers are the basic components of composite laminated structures. Therefore, the stress-strain relationship of the fiber reinforced layer of the laminated plate is focused. It assumes that the proposed plate has N layers. Then, the relationship of stress-strain for the kth layer can be expressed asin which (i = r and θ) and (j = rθ, θr, rz, and θz) denote normal stresses and shear stresses, respectively. Besides, the stiffness coefficients (p and l = 1, 2, 4, 5, and 6) are expressed asin which the angle between the fiber direction and the radius direction is defined as . is the material coefficient and is given aswhere and are Poisson’s ratios. Young’s modulus is expressed as and . The relationship between Poisson’s ratio and Young’s modulus is . Besides, the shear modulus is denoted as , , and . When and , the isotropic plate can be simulated.

So far, the relations between generalized force and strain can be written as the matrix form:in which A_pl, B_pl, C_pl, D_pl, E_pl, F_pl, and are marked as the stiffness coefficients. Besides, Z_k is the thickness coordinate value of the under surface of the kth layer plate. Z_k+1 is the thickness coordinate value of the upper surface of the kth layer plate. So, the thickness of this layer plate is z_k = Z_k+1 − Z_k.

2.3. Energy Equation of Laminated Annular/Circular Plate

In order to decrease the requirement of displacement function selection, a kind of multisegment segmentation technique is extended which can effectively relax the boundary constraints of the laminated annular or circular plate. Specifically, the laminate structure is subdivided into N_r segments along the r-axis direction. Then, a penalty parameter is introduced to ensure the continuity and boundary constraints between segments. Firstly, the general expression of the Lagrange energy equation L for composite laminated annular or circular plates is given, which fully considers the effect of boundary constraints:where the kinetic energy and potential energy of the qth segment plate are expressed by the symbols and . When the plate is segmented, the coupling potential energy is introduced, which is represented by symbol . For the whole plate, the total kinetic energy and total potential energy should be the superposition of N_r segments and , while the total coupling potential energy should be the superposition of (N_r − 1) segment . According to the boundary spring simulation technique, the spring potential stored at the plate boundary can be represented by symbol . Next, according to Reddy’s HSDT, their specific expressions are given.

For the qth segment plate, the kinetic energy is expressed asin which

The specific expression of stored in the qth segment laminated plate is

As mentioned above, when the plate structure is segmented along the radius direction, the coupling potential energy between adjacent segments is needed to ensure the continuity of displacements. Similarly, the coupling boundary is simulated by virtual spring technology. In short, the coupling springs (, , , , , , and ) are added between adjacent segments. Then, the reasonable spring stiffness is set, which will be further discussed in the following numerical analysis. For the qth segment and (q + 1)th segment plate, the coupling potential energy can be written as

There are a lot of research results on the elastic boundary conditions of plates and shells by using artificial virtual spring technology [28–31]. In Figure 1(b), the setting and types of elastic boundary springs have been given. The expressions of coupling potential energy of the boundary springs are given directly below:

So far, the energy equation of the elastic laminate structure has been clearly expressed. Next, we will focus on the setting of the displacement function on the reference midsurface and the specific solving process of the natural frequencies and mode shapes.

2.4. Solution Process of Laminated Annular/Circular Plate

In the previous section, the global energy equation of the laminated annular or circular plate has been established, so this section will solve it by using the Hamilton differential principle. Firstly, it is necessary to determine the concrete displacement function expression of the middle plane. This paper applies the multisegment segmentation technique to investigate the free vibration behaviors of composite annular or circular plates, which can reduce the requirements of choosing displacement functions. Because compared with directly building the displacement function of the whole structure, the displacement function only needs the low order polynomials for the segment structure. Here, the Jacobi orthogonal polynomials are applied to construct displacement functions [32–36], which can be expressed aswhere the Jacobi coefficients are expressed as , , , , and for the qth segment plate. The Jacobi polynomial of the mth in the radius direction is set to , and α and β denote the type of polynomials. Jacobi polynomials are orthogonal polynomials, so we can get Legendre polynomials, Chebyshev polynomials, and so on by setting the values of α and β. For example, if the values of α and β are 0, the Legendre polynomial will be obtained. When α = β = −0.5, it is the first kind of Chebyshev polynomial. When α = β = 0.5, it is the second kind of Chebyshev polynomial. This also shows the universality of the Jacobi polynomials. In addition, n represents the number of half wave in the rotation direction. Then, j, ω, and t are pure imaginary number, circular frequency, and time variable.

According to the Hamilton differential principle, the partial derivative of L to the unknown Jacobi coefficient matrix E is zero, which can be written as

By introducing equations (11)–(15) into (21), a vibration characteristic solution equation can be obtained:in which K, K_cp, and K_BC are stiffness matrix, coupling matrix between segments, and boundary spring matrix, separately. M denotes the mass matrix. From this equation, it is easy to find that the circular frequencies can be obtained by finding the generalized eigenvector. In addition, the corresponding generalized eigenvectors are introduced in equations (19)–(23), and the modal shapes of the laminated plate could be easily gained.

3. Numerical Analysis and Discussion

By effectively combining Reddy’s HSDT and Jacobi polynomials, the equations for solving vibration behaviors of laminated annular or circular plates are established. Next, some numerical examples are applied to study the free vibration characteristics of the laminated plate structures. This section mainly studies the convergence and accuracy, new results, and experimental tests. In addition, it should be pointed out that different boundary constraints are represented by a combination of corresponding first characters. Specifically, since there is only one boundary in the direction of the outer radius for a circular plate, there is only one character to represent it. Among them, the character sign corresponding to the fixed support is C, the sign corresponding to the simple support is S, the sign corresponding to the free boundary is F, and the sign corresponding to the elastic support is E. For annular plates with fixed support in the inner radius direction and simply supported in the outer radius direction, the boundary constraint can be expressed as C-S.

3.1. Model Validation

Convergence and accuracy analysis are important indexes for evaluating an algorithm. Therefore, this section focuses on the convergence of this method and the accuracy of the analysis model in the study of vibration characteristics. By parameterizing the stiffness of the boundary spring, different boundary constraints are defined. In addition, the parameterization of the stiffness of the coupling spring between segments is examined. From the expression of the displacement function, it is not difficult to find that the truncated value M of Jacobi polynomial and segment number N_r directly represent the convergence of the proposed method. In Figure 2, the convergence of frequency parameters for the laminated annular or circular plate with fixed boundary constraints is investigated. For the laminated plate, the material parameters are set to E₁ = 20E₂, E₂ = 10 GPa, G₁₂ = G₁₃ = 0.6E₂, G₂₃ = 0.5E₂, μ₁₂ = 0.25, and ρ = 1600 kg/m³. The lamination scheme is [0°/90°], and the geometric parameters are R₁ = 1 m and h = 0.1 m. In addition, the internal radius R₀ is 0.3 m. Figure 2 shows the curves of the frequency parameters of the annular plate and the circular plate with the increase of M or N_r. It shows intuitively that the frequency parameters of the two types of plate structures converge when M is 8 or N_r is 2. The proposed method shows good and consistent convergence in solving the natural frequency of this structure. In order to ensure the accuracy of the numerical example, the truncated value in the following examples is defined as M = 10 and N_r = 2.

(a)

(b)

(c)

(d)

As previously described, different types of polynomials could be obtained by setting the values of Jacobi parameters α and β. Table 2 shows error analysis of frequency parameter for the clamped composite laminated annular plate with various Jacobi parameters α and β. The material parameters are as follows: E₁ = 15E₂, E₂ = 10 GPa, G₁₂ = G₁₃ = 0.6E₂, G₂₃ = 0.5E₂, μ₁₂ = 0.25, and ρ = 1500 kg/m³. The lamination scheme is [0°/90°0°/90°], and the geometric parameters are R₀ = 1 m, R₁ = 3 m, and h_p = 0.1 m. It is not difficult to find out from Table 2 that the Jacobi parameters α and β have little effect on the calculation results of the frequency parameter Ω under the accuracy of five effective digits. In addition, the results obtained in this paper are in excellent agreement with those in [37]. If there is no special description, the Jacobi parameters are chosen as α = β = 0 in other calculation examples in this paper.

As mentioned above, the complex boundary conditions are simulated by artificial virtual spring technology. Next, Figure 3 displays the influence of boundary spring values on natural frequency of the laminated annular plate and a circular plate. Through numerical analysis, the variation range of spring stiffness for the elastic boundary can be determined. Except for the boundary constraints, the structural and material parameters of Figure 3 are exactly the same as those used in Table 2. The range of values of the linear boundary spring (defined as k) and rotary boundary spring (defined as K) is both 10⁻⁴ to 10¹⁶. It can be seen from Figure 3 that the waterfall plots of the frequencies to the change of boundary spring stiffness values have the same trends. In brief, when the linear boundary spring stiffness value (k) and the rotary boundary spring stiffness value (K) are less than 10⁴, the frequency value is 0 corresponding to the free boundary condition. When k and K are higher than 10¹⁰, the frequency reaches the peak corresponding to the clamped boundary condition. Besides, the elastic boundary conditions can be determined when k and K appear in the following four combinations: 10⁴<k < 10¹⁰, K < 10⁴; 10⁴ < k < 10¹⁰, K > 10¹⁰; 10⁴ < K < 10¹⁰, k < 10⁴; 10⁴ < K < 10¹⁰, k > 10¹⁰. In light of the conclusions above, the values k and K corresponding to different boundary conditions are given in Table 3.

(a)

(b)

(c)

(d)

Now, let us focus on the connective spring. Figure 4 displays the influence of the connective spring stiffness value on the natural frequency of laminated annular and circular plates. In Figure 4, the annular plate is simply supported and the circular plate is clamped. Other parameters are exactly the same with Figure 3. Stiffness values of the connective spring vary from 10⁻² to 10¹⁶. The results show that the effect of connective spring stiffness on the natural frequencies is the same whether it is an annular plate under the SS boundary constraint or a circular plate under the C boundary constraint. The frequencies are almost unaffected when the value of the connective spring stiffness is less than 10⁴. Thereafter, the frequencies increase with the increase of the connective spring stiffness until the value of connective spring stiffness is greater than 10¹⁰. To ensure continuity, the value of connective spring stiffness is chosen as 10¹⁰.

(a)

(b)

So far, a numerical analysis model for free vibration behaviors of the composite laminated annular or circular plate has been established. To further validate the correctness of the proposed method, Table 4 gives the comparisons of the frequency parameter for the laminated annular plate between the results obtained by HSDT and existing research [37]. Obviously, these two methods have a great agreement, which proves the correctness of the present vibration analysis model. In addition, it is not difficult to see from this table that the stronger the boundary constraint, the greater the natural frequency of the plate.

Then, Table 5 shows the comparison of the frequency parameter for the isotropic circular plate with different thickness ratio h/R₁. The parameters of the clamped circular plate are E₁ = E₂ = 185 GPa, G₁₂ = G₁₃ = G₂₃ = 7.12 GPa, μ₁₂ = 0.3, ρ = 1600 kg/m³, and R₁ = 1 m. The results of FEM in the table are obtained by ABAQUS software. It is not difficult to find that the results obtained by the present method are in good agreement with those obtained by FEM when h/R₁ is 0.2. In addition, the material used throughout this numerical example is actually an isotropic material. Therefore, the accuracy and universality of the present method in solving isotropic plate structures are also verified by this example. From this table, it can be found that the natural frequency of the plate decreases with the increase of the ratio of h/R₁.

3.2. Study on Free Vibration Characteristics

In Section 3.1, the vibration analysis model for the composite annular and circular plate based on Reddy’s HSDT and Jacobi polynomials method has been established. The correctness and universality of the present method have been validated, too. In this part, we will give some new results and conclusions.

Firstly, Table 6 shows the frequency parameter for the composite annular plate with various radius ratios and boundary constraints. The material parameters are E₁ = 30E₂, E₂ = 10 GPa, G₁₂ = G₁₃ = 0.6E₂, G₂₃ = 0.5E₂, μ₁₂ = 0.27, and ρ = 1600 kg/m³. And, the lamination scheme is [60°/−30°/60°]. The invariant geometric parameters are R₁ = 1 m and h = 0.15 m. It can be found from Table 6 that radius ratios and boundary constraints have a profound impact on the frequency parameter. On the one hand, the frequency parameter increases when the boundary spring stiffness value increases. On the other hand, the influence of the radius ratio on the frequency parameter is more complex and does not show monotonicity. Then, Table 7 gives the frequency parameter for the laminated annular plate with various boundary constraints and thicknesses. The material parameters involved in Table 7 are E₁ = 15E₂, E₂ = 10 GPa, G₁₂ = G₁₃ = 0.5E₂, G₂₃ = 0.2E₂, μ₁₂ = 0.25, and ρ = 1500 kg/m³. The lamination scheme is [30°/−30°], and R₁ = 1 m. It is noticeable that the frequency parameter increases when the boundary spring stiffness and plate thickness increase.

Next, Table 8 gives the frequency parameter for the laminated annular plate and a circular plate with diverse anisotropy degrees and lamination schemes. The material parameters implicated in Table 8 are E₂ = 10 GPa, G₁₂ = G₁₃ = 0.5E₂, G₂₃ = 0.6E₂, μ₁₂ = 0.25, and ρ = 1500 kg/m³. The geometrical parameters are R₁ = 1 m, h = 0.08 m, and R₀ = 0.3 m (for the annular plate). The boundary constraints of the annular plate and circular plate are E2-E2 and E2, respectively. It is not difficult to find that the frequency parameter increases with the increase of anisotropy degrees. It is noteworthy that the lamination scheme has an important influence on the frequency parameters. However, due to data limitations, it is difficult to draw general conclusions.

For the sake of fully revealing the effect of material parameters on vibration behaviors of the laminated annular plate and circular plate, the parametric research is conducted. Figure 5 studies the effect of the layer number on natural frequencies of the laminated annular plate and circular plate with fixed boundary constraints. The material parameters are E₁ = 5E₂, E₂ = 10 GPa, G₁₂ = G₁₃ = 0.5E₂, G₂₃ = 0.6E₂, μ₁₂ = 0.25, and ρ = 1500 kg/m³. The geometrical parameters of the clamped plates are R₁ = 1 m, h = 0.08 m, and R₀ = 0.3 m (for an annular plate). It should be noted that p represents the layer number. And, [0/−30°]₁ represents [0/−30°], while [0/−30°]₂ represents [0/−30°/0/−30°], and so on. As can be seen from Figure 5, the frequencies increase when the layer number increases and tend to be stable after the layer number exceeds 4. Moreover, the frequency changes most significantly when the layer number changes from 1 to 2.

(a)

(b)

(c)

(d)

In addition, for the annular plate, its frequency decreases monotonously with the increase of the ply angle, but not for the circular plate. In order to reveal the effect of the ply angle on the vibration behaviors of the laminated annular and circular plate with fixed boundary constraints, Figure 6 shows the variation of frequency with various ply angles. The lamination scheme used in Figure 6 is chosen as [0/°], and α ranges from 10° to 170°. From Figure 6, it is easy to find that the frequency variation is symmetrical with respect to [0/90°]. On the contrary, when α is less than 90°, the frequency of the annular plate decreases monotonously with the increase of . For the circular plate, the frequency is not significantly affected by the change of the ply angle and does not have monotonicity.

(a)

(b)

3.3. Model Experimental Study

In Section 3.2, free vibration behaviors of the laminated annular plate and circular plate are explored in depth based on the high-order shear deformation theory proposed by Reddy and Jacobi polynomials method. However, the previous studies only consider numerical validation. There is a lack of experimental test data to support the qualitative conclusion of the correctness of the above analysis model. In order to more fully verify the correctness of this study, this part applies the model experimental method. At present, equipment and principle of the structural model test can give reliable and more accurate experimental data. In addition, the fixed or free boundary conditions are insensitive to the external environment and easy to get realized. So, the model experiments of annular and circular plates are completed under the fixed boundary condition or the combination of free and fixed boundary conditions, which can provide more accurate experimental results. The experimental details on circular and annular plates can be found in [38, 39].

In this paper, the force measurement method is used for the model experiment. The measurement method is a single-point vibration pickup method, that is, the position of the acceleration sensor remains unchanged all the time, and the hammer strikes all measuring points in turn. This test method requires that the vibration pick-up point avoid the joint position of the test piece and tries to select the position with large deformation. During the experiment, the 32-channel DH5922N dynamic signal test and analysis system and the corresponding modal analysis software are used, which were produced by Jiangsu Donghua Testing Technology Co., Ltd. In addition, the general piezoelectric acceleration sensor (1A116E) and the charge amplifier (DH5857-1) are also needed. Table 9 shows the experimental instruments and their related information used in the model experiment.

In this experiment, according to the experimental schematic diagram in Figure 7(a), the above experimental equipment and workpiece are installed and arranged. Figures 7(b) and 7(c) show the connection diagram of experimental equipment and test object of the annular plate and circular plate. In this experiment, the 30 cm thick foundation frame, 10 cm thick batten, and annular plate or circular plate are connected by 60 bolts with the same diameter 8 mm to simulate the fixed boundary conditions. Table 10 gives specifications of test plates and arrangement. It should be noted that the geometric dimensions given in this table are all effective dimensions. The first test step is to select a point on the plate as the vibration pick-up point and place the acceleration sensor on it. In the second step, the hammer is used to excite each grid intersection on the plate, and the excitation force signal is measured by the force sensor on the hammer. Then, the acceleration sensor at the pick-up point of the plate collects the vibration acceleration signal and amplifies the signal through the charge adapter. Finally, using the force measurement method modal analysis software in the computer, modal parameters identification and frequency response function analysis are carried out on the experimental data. Through frequency response curve fitting, steady-state diagram calculation, and mode calculation, the vibration mode and natural frequency of the plate are finally obtained.

(a)

(b)

(c)

Figures 8 and 9 show the difference of experimental test results, FEM results, and present results. Among them, the maximum deviation of the experimental test results and present results is 7.57% for the natural frequency of the annular plate and 2.57% for the natural frequency of the circular plate. Besides, the maximum deviation of the experimental test results and FEM results is 9.83% for the natural frequency of the annular plate and 1.90% for the natural frequency of the circular plate. For the prediction results of the model shape, the two methods are highly consistent. The deviation of the above experimental results is within the acceptable range, which fully proves the correctness of the analysis model for the vibration behaviors of annular and circular plates. The experimental error is caused by many reasons. First of all, the fixed boundary conditions of the plate cannot be fully simulated by the way of clamping the foundation frame and the battens. Then, the material parameters used in the numerical calculation of the plate deviate from the actual material parameters of the work piece, and the work piece cannot be completely ideal isotropic material. In addition, the accuracy deviation of the force sensor and acceleration sensor and the human error of the experimenter in the process of hammering will cause the error of experimental data.

(a)

(b)

(c)

(a)

(b)

(c)

4. Conclusions

This paper applies the Jacobi polynomial method to establish a unified analytical model for the vibration characteristics of laminated annular and circular plates under complex elastic boundary constraints based on Reddy’s HSDT. Due to the adoption of multisegment segmentation technique, coupling springs are applied between adjacent segments to ensure continuity. In addition, various elastic boundary constraints are simulated by setting boundary springs at the edges of the plate. The Lagrangian energy equation of composite laminated annular or circular plates under elastic boundary constraints is established. The solution equation of free vibration behavior for the laminated plate is easily obtained by the Hamilton differential principle. On the basis of verifying the fast convergence and great accuracy of this method, further parametric study is carried out. Some important conclusions are revealed:(1)Through the analysis of the stiffness value of the boundary springs, the value range of the elastic boundary condition is determined which are 10⁴ < k < 10¹⁰ and K < 10⁴, 10⁴ < k < 10¹⁰ and K > 10¹⁰, 10⁴ < K < 10¹⁰ and k < 10⁴, and 10⁴ < K < 10¹⁰ and k > 10¹⁰. When θ = 2π, the value of connective spring stiffness should be chosen as 10¹⁰ to ensure continuity.(2)Material parameters, geometric parameters, and lamination schemes have an important effect on the vibration characteristics of the annular and circular plates. Specifically, the natural frequency will increase with the increase of plate thickness and anisotropy degrees.(3)The frequency increases with the increase of the layer number and tends to be stable after the layer number exceeds 4. Moreover, the frequency changes most significantly when the layer number changes from 1 to 2.(4)The deviation of the model experimental results is within the acceptable range, which fully proves the correctness of the unified analysis model for the vibration behaviors of annular and circular plates.

Data Availability

The (TIF, OPTU, and DOCX) data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This project was supported by the National Key R&D Program of China (Grant no. 2018YFB2000100). The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (Grant no. 52005255), Natural Science Foundation of Jiangsu Province of China (Grant no. BK20200430), and Fundamental Research Funds for the Central Universities (Grant no. NS2020033).

References

B. Afsharmanesh, A. Ghaheri, and F. Taheri-Behrooz, “Buckling and vibration of laminated composite circular plate on Winkler-type foundation,” Steel and Composite Structures, vol. 17, no. 1, pp. 1–19, 2014.
View at: Publisher Site | Google Scholar
M. Amabili, R. Pierandrei, and G. Frosali, “Analysis of vibrating circular plates having non-uniform constraints using the modal properties of free-edge plates: application to bolted plates,” Journal of Sound and Vibration, vol. 206, no. 1, pp. 23–38, 1997.
View at: Publisher Site | Google Scholar
E. Arshid and A. R. Khorshidvand, “Free vibration analysis of saturated porous FG circular plates integrated with piezoelectric actuators via differential quadrature method,” Thin-Walled Structures, vol. 125, pp. 220–233, 2018.
View at: Publisher Site | Google Scholar
A. Bahrami and A. Teimourian, “Free vibration analysis of composite, circular annular membranes using wave propagation approach,” Applied Mathematical Modelling, vol. 39, no. 16, pp. 4781–4796, 2015.
View at: Publisher Site | Google Scholar
K. Mercan, A. K. Baltacıoglu, and Ö. Civalek, “Free vibration of laminated and FGM/CNT composites annular thick plates with shear deformation by discrete singular convolution method,” Composite Structures, vol. 186, pp. 139–153, 2018.
View at: Publisher Site | Google Scholar
K. K. Viswanathan, S. Javed, Z. A. Aziz, and K. Prabakar, “Free vibration of symmetric angle-ply laminated annular circular plate of variable thickness under shear deformation theory,” Meccanica, vol. 50, no. 12, pp. 3013–3027, 2015.
View at: Publisher Site | Google Scholar
A. Powmya and M. C. Narasimhan, “Free vibration analysis of axisymmetric laminated composite circular and annular plates using Chebyshev collocation,” International Journal of Advanced Structural Engineering, vol. 7, no. 2, pp. 129–141, 2015.
View at: Publisher Site | Google Scholar
K. Draiche, A. A. Bousahla, A. Tounsi, A. S. Alwabli, A. Tounsi, and S. Mahmoud, “Static analysis of laminated reinforced composite plates using a simple first-order shear deformation theory,” Computers and Concrete, vol. 24, pp. 369–378, 2019.
View at: Google Scholar
M. Amabili, P. Balasubramanian, R. Garziera, and G. Royer-Carfagni, “Blast loads and nonlinear vibrations of laminated glass plates in an enhanced shear deformation theory,” Composite Structures, vol. 252, Article ID 112720, 2020.
View at: Publisher Site | Google Scholar
M. Amabili, “A new third-order shear deformation theory with non-linearities in shear for static and dynamic analysis of laminated doubly curved shells,” Composite Structures, vol. 128, pp. 260–273, 2015.
View at: Publisher Site | Google Scholar
M. Amabili, Nonlinear Mechanics of Shells and Plates in Composite, Soft and Biological Materials, Cambridge University Press, Cambridge, UK, 2018.
M. Amabili and J. N. Reddy, “The nonlinear, third-order thickness and shear deformation theory for statics and dynamics of laminated composite shells,” Composite Structures, vol. 244, Article ID 112265, 2020.
View at: Publisher Site | Google Scholar
H. Bisadi, M. Es’haghi, H. Rokni, and M. Ilkhani, “Benchmark solution for transverse vibration of annular Reddy plates,” International Journal of Mechanical Sciences, vol. 56, no. 1, pp. 35–49, 2012.
View at: Publisher Site | Google Scholar
E. Asadi and J. Jam, “Analytical and numerical buckling analysis of carbon nanotube reinforced annular composite plates,” International Journal of Advanced Design and Manufacturing Technology, vol. 7, p. 35, 2014.
View at: Google Scholar
M. Vinyas, A. Sandeep, T. Nguyen-Thoi, F. Ebrahimi, and D. Duc, “A finite element-based assessment of free vibration behaviour of circular and annular magneto-electro-elastic plates using higher order shear deformation theory,” Journal of Intelligent Material Systems and Structures, vol. 30, no. 16, pp. 2478–2501, 2019.
View at: Publisher Site | Google Scholar
S. Hosseini-Hashemi, M. Es’haghi, and H. Rokni Damavandi Taher, “An exact analytical solution for freely vibrating piezoelectric coupled circular/annular thick plates using Reddy plate theory,” Composite Structures, vol. 92, no. 6, pp. 1333–1351, 2010.
View at: Publisher Site | Google Scholar
S. Hosseini-Hashemi, M. Es’haghi, H. Rokni Damavandi Taher, and M. Fadaie, “Exact closed-form frequency equations for thick circular plates using a third-order shear deformation theory,” Journal of Sound and Vibration, vol. 329, no. 16, pp. 3382–3396, 2010.
View at: Publisher Site | Google Scholar
M. M. Najafizadeh and H. R. Heydari, “An exact solution for buckling of functionally graded circular plates based on higher order shear deformation plate theory under uniform radial compression,” International Journal of Mechanical Sciences, vol. 50, no. 3, pp. 603–612, 2008.
View at: Publisher Site | Google Scholar
M. Abualnour, A. Chikh, H. Hebali et al., “Thermomechanical analysis of antisymmetric laminated reinforced composite plates using a new four variable trigonometric refined plate theory,” Computers and Concrete, vol. 24, pp. 489–498, 2019.
View at: Google Scholar
O. Allam, K. Draiche, A. A. Bousahla et al., “A generalized 4-unknown refined theory for bending and free vibration analysis of laminated composite and sandwich plates and shells,” Computers and Concrete, vol. 26, pp. 185–201, 2020.
View at: Google Scholar
N. Belbachir, M. Bourada, K. Draiche et al., “Thermal flexural analysis of anti-symmetric cross-ply laminated plates using a four variable refined theory,” Smart Structures and Systems, vol. 25, pp. 409–422, 2020.
View at: Google Scholar
N. Belbachir, K. Draich, A. A. Bousahla, M. Bourada, A. Tounsi, and M. Mohammadimehr, “Bending analysis of anti-symmetric cross-ply laminated plates under nonlinear thermal and mechanical loadings,” Steel and Composite Structures, vol. 33, pp. 81–92, 2019.
View at: Google Scholar
M. Sahla, H. Saidi, K. Draiche, A. A. Bousahla, F. Bourada, and A. Tounsi, “Free vibration analysis of angle-ply laminated composite and soft core sandwich plates,” Steel and Composite Structures, vol. 33, pp. 663–679, 2019.
View at: Google Scholar
P. V. Katariya, S. K. Panda, and K. Mehar, “Theoretical modelling and experimental verification of modal responses of skewed laminated sandwich structure with epoxy-filled softcore,” Engineering Structures, vol. 228, Article ID 111509, 2021.
View at: Publisher Site | Google Scholar
K. Mehar, S. K. Panda, and T. R. Mahapatra, “Theoretical and experimental investigation of vibration characteristic of carbon nanotube reinforced polymer composite structure,” International Journal of Mechanical Sciences, vol. 133, pp. 319–329, 2017.
View at: Publisher Site | Google Scholar
J. N. Reddy, “A refined nonlinear theory of plates with transverse shear deformation,” International Journal of Solids and Structures, vol. 20, no. 9-10, pp. 881–896, 1984.
View at: Publisher Site | Google Scholar
J. N. Reddy, “A simple higher-order theory for laminated composite plates,” Journal of Applied Mechanics, vol. 51, no. 4, pp. 745–752, 1984.
View at: Publisher Site | Google Scholar
T. Ye, G. Jin, and Y. Zhang, “Vibrations of composite laminated doubly-curved shells of revolution with elastic restraints including shear deformation, rotary inertia and initial curvature,” Composite Structures, vol. 133, pp. 202–225, 2015.
View at: Publisher Site | Google Scholar
C. Zhang, G. Jin, T. Ye, and Y. Zhang, “Harmonic response analysis of coupled plate structures using the dynamic stiffness method,” Thin-Walled Structures, vol. 127, pp. 402–415, 2018.
View at: Publisher Site | Google Scholar
Z. Su, G. Jin, and T. Ye, “Electro-mechanical vibration characteristics of functionally graded piezoelectric plates with general boundary conditions,” International Journal of Mechanical Sciences, vol. 138-139, pp. 42–53, 2018.
View at: Publisher Site | Google Scholar
Z. Su, G. Jin, and T. Ye, “Vibration analysis of multiple-stepped functionally graded beams with general boundary conditions,” Composite Structures, vol. 186, pp. 315–323, 2018.
View at: Publisher Site | Google Scholar
B. Qin, R. Zhong, Q. Wu, T. Wang, and Q. Wang, “A unified formulation for free vibration of laminated plate through Jacobi-Ritz method,” Thin-Walled Structures, vol. 144, Article ID 106354, 2019.
View at: Publisher Site | Google Scholar
B. Qin, K. Choe, T. Wang, and Q. Wang, “A unified Jacobi-Ritz formulation for vibration analysis of the stepped coupled structures of doubly-curved shell,” Composite Structures, vol. 220, pp. 717–735, 2019.
View at: Publisher Site | Google Scholar
Q. Wang, K. Choe, J. Tang, C. Shuai, and A. Wang, “Vibration analyses of general thin and moderately thick laminated composite curved beams with variable curvatures and general boundary conditions,” Mechanics of Advanced Materials and Structures, vol. 27, pp. 991–1005, 2019.
View at: Publisher Site | Google Scholar
Y. Kurylov and M. Amabili, “Polynomial versus trigonometric expansions for nonlinear vibrations of circular cylindrical shells with different boundary conditions,” Journal of Sound and Vibration, vol. 329, no. 9, pp. 1435–1449, 2010.
View at: Publisher Site | Google Scholar
Y. Kurylov and M. Amabili, “Nonlinear vibrations of clamped-free circular cylindrical shells,” Journal of sound and vibration, vol. 330, no. 22, pp. 5363–5381, 2011.
View at: Publisher Site | Google Scholar
G. Jin, T. Ye, X. Jia, and S. Gao, “A general Fourier solution for the vibration analysis of composite laminated structure elements of revolution with general elastic restraints,” Composite Structures, vol. 109, pp. 150–168, 2014.
View at: Publisher Site | Google Scholar
E. Askari, K.-H. Jeong, and M. Amabili, “Hydroelastic vibration of circular plates immersed in a liquid-filled container with free surface,” Journal of sound and vibration, vol. 332, no. 12, pp. 3064–3085, 2013.
View at: Publisher Site | Google Scholar
M. Kwak and M. Amabili, “Hydroelastic vibration of free-edge annular plates,” Journal of Vibration and Acoustics, vol. 121, pp. 26–32, 1999.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Guodong Jin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies