Research Article  Open Access
Vu H. Nam, Nguyen H. Nam, Pham V. Vinh, Dang N. Khoa, Do V. Thom, Phung V. Minh, "A New Efficient Modified FirstOrder Shear Model for Static Bending and Vibration Behaviors of TwoLayer Composite Plate", Advances in Civil Engineering, vol. 2019, Article ID 6814367, 17 pages, 2019. https://doi.org/10.1155/2019/6814367
A New Efficient Modified FirstOrder Shear Model for Static Bending and Vibration Behaviors of TwoLayer Composite Plate
Abstract
A twolayer (connected by stubs) partial composite plate is a structure with outstanding advantages which can be widely applied in many fields of engineering such as construction, transportation, and mechanical. However, studies are scarce in the past to investigate this type of structure. This paper is based on the new modified firstorder shear deformation plate theory and finite element method to develop a new fournode plate element with nine degrees of freedom per node for static bending and vibration analysis of the twolayer composite plate. The numerical results are compared to published data for some special cases. The effects of some parameters such as the boundary condition, stiffness of the connector stub, heighttowidth ratio, thicknesstothickness ratio between two layers, and aspect ratio are also performed to investigate new numerical results of static bending and free vibration responses of this structure.
1. Introduction
Twolayer beam and plate are among the most commonly used structures in many fields of engineering, such as construction, transportation, and mechanical, because of their advantages in comparison with onelayer structures; we can see these structures in practice such as steelconcrete composite beams and plates, layered wooden beams and plates, and woodconcrete and timbersteel floor structures. Some advantages of this type of composite structure can be considered such as easy manufacturing, taking advantage of the material properties of two components which made the structure. In fact, there are many different grafting methods to manufacture twolayer beams or plates; for example, two components of a structure are sandwiched together at the edges, using glue or connector stub to bond two parts at the contact surface. Because of these advantages, many researchers have focused on the mechanical analysis of these structures.
Twolayer beam has been investigated by Foraboschi [1, 2] using an analytical method. Kryžanowski et al. [3] developed an analytical method to evaluate exact critical forces of twolayer composite columns. Grognec et al. [4] used the Timoshenko beam model for buckling analysis of twolayer composite beams with partial interaction. He and Yang [5, 6] presented the finite element method and higher order beam theory for buckling and dynamic problem of the twolayer composite beam. He et al. [7] carried out an analytical solution for free vibration and buckling of two and threelayer composite beams based on a higher order beam theory. Hozan et al. [8] analyzed a twolayer composite planar beam with exact geometric and material nonlinearities as well as the finite slip between the layers. Wu and coworkers [9, 10] used the 2D elasticity solution for analysis of twolayer beams with the viscoelastic interlayer and arbitrarily shaped interface. Large deformation analysis of twolayer composite beams was studied [11]. Based on the higher order beam theory and analytical model, Wen et al. [12] investigated the flexural response of the twolayer composite beam. Do et al. [13] presented static bending analysis of threelayer beams made of functionally graded materials in the thermal environment based on FEM and Mindlin theory. Cas et al. [14] used an analytical method to investigate the mechanical behavior of twolayer composite beams with interlayer slips. Ke et al. [15] used experimental and numerical FE methods to study local web buckling of singlecoped beam connections with the slender web.
Latham et al. [16] developed a new finite element based on a new plate theory due to the Reissner plate theory for static bending of twolayer plates. Foraboschi and colleagues [17, 18] used an analytical model and experimental method for static bending of laminated glass plates. Layered plate with discontinuous connection was investigated by Foraboschi [19] using the exact mathematical model replacing discontinuous connections with a fictitious continuous medium. Foraboschi [20, 21] also detailed the mechanical behavior of threelayer plates using the exact mathematical model and elasticity solution. Wu et al. [22] applied the 3D exact solution to explore the twolayer plate bonded by a viscoelastic interlayer. Vidal et al. [23] presented a new approach to study the deflection and stress of the multilayer composite plate. Alinia et al. [24] applied the numerical nonlinear finite element method and theoretical pRitz energy method to analyze inelastic buckling and postbuckling behaviors of stocky plates under combined shear and inplane bending stresses.
This paper aims to use the finite element method (FEM) based on a new modified firstorder shear deformation plate theory (FSDT) to analyze static bending and free vibration of the twolayer composite plate. The proposed method shows the simple formulations and computational efficiency. The accuracy and reliability of the proposed method are validated with other published results. Several numerical examples and influence of some parameters on static bending and free vibration of the twolayer composite plate are also investigated; these new results have significant impacts on the use of this structure in practice.
The organization of this paper is as follows: Section 2 presents the finite formulation of static bending and free vibration of the twolayer composite plate based on a new modified firstorder shear deformation plate theory. New numerical results for bending and free vibration analysis of these plates are computed and discussed in Section 3. Some major conclusions are given in Section 4.
2. Geometry and Theoretical Formulation
2.1. Geometry and Assumptions
The geometry model of the problem is a twolayer plate including two isotropic plates which bonded together at the contact surface using connector stubs, as shown in Figure 1. Two parts are of the same size in the plane with thickness for the bottom layer and for the top layer. To avoid the repetition of similar equations for the bottom layer and top layer, we use the subscript for the bottom layer, subscript for the top layer, and subscript for the connector stub.
The assumptions of the twolayer composite plate include the following: the materials of each layer are linear, elastic, and isotropic; the displacement and rotation of the plate are small; there is no delamination phenomenon between two layers; the deflection of the stub is in the contact interface of two components; and the mass of the stub is much smaller than the mass of the plate, so we assume that it is neglected.
2.2. Discrete Finite Element Equations
2.2.1. New Modified FirstOrder Shear Deformation for Plate Elements
The basic equations of the Mindlin plate theory are [25]wherein which k = 5/6 is the shear correction.
Equilibrium of moments about  and axis and transverse force leads to
Substituting equation (1) into equation (3), we get the following:where is the Laplace differential operator.
Assume that the total deflection consists of two parts which are bending deflection and transverse shear, while angles of the plate crosssectional slope are a result of rotation of pure bending and shear angles:where are rotations due to pure bending and are the shear angles.
Substituting equation (5) into equation (4), we get two differential equations:
Equations (6) and (7) can be rewritten in the following forms:and their integrals per and read according to the structure of in equations (6) and (7), the reasonable solution is that both and must be set to zero; as a consequence of that consideration, the relation between transverse shear and bending is
The total deflection is
Using the new modified firstorder shear deformation theory, the displacement field of the plate can be rewritten in the following form [26]:
The straindisplacement relations may be written as follows [26]:
Equation (13) can be rewritten in the matrix form aswhereand
2.2.2. Finite Element Formulations
The plate is discretized using the quadrilateral fournode element, as shown in Figure 2. Each node has nine degrees of freedom including four axial, one transversal, and four rotational displacements.
(a)
(b)
The element nodal displacement vector can be given aswhere
The nodal displacement vector of each component iswhere
In equations (18) and (20), are, respectively, the axis displacements and transversal displacement and denote the bending rotations of the axis, respectively.
The static bending deflection is assumed in a polynomial form. Then, the transverse shear and bending rotations are calculated using the above equations, while the axial displacements are approximated using Lagrangian shape functions.
The displacements may be approximated aswhere
The bending deflection is approximated using a polynomial form aswhere is a vector of polynomial coefficients andwhere According to equation (10), the transverse shear is given bywherein which .
The total static deflection of the plate is
The rotations of cross sections on the neutral plane are
Taking coordinate values for each node into equations (27)–(29), we getwhere and matrix is shown in Appendix.
Taking equation (30) into account, we get
Substituting equations (21) and (31) into equation (14), we obtainwhere
By substituting equation (25) into equation (16), we getwhere
According to Hooke’s law, the relationship between stresses and strains can be given aswhere
The calculation model for the connector stub is shown in Figure 3.
By using Lagrangian approximation for both axial displacement and rotation of two components to determine the axial deflection of the connector stub, the interfacial slips of two components are given bywherein which u_{st} and are the deflections of the connector stub.
The strain energy of two components and connector stub can be given by
The kinetic energy isor
Equation (46) can be rewritten in the shorter form as follows:where .
The work done by external distribution normal force acting on the top surface of layer iswhere q is the uniform load.
The Lagrangian function of the plate element is defined by
Using Lagrange’s equations, the governing equations of motion of the plate element are given by
Substituting equations (44), (46), and (48) into equation (50) considering equation (49), the equation of motion of the plate can be obtained in the matrix form asin which , and are, respectively, the element mass matrix, the element stiffness matrix, and element force vector, which are defined as
Note that we need not any selective reduced integration or reduced integration scheme to calculate the matrices and vectors in equations (52)–(55). When k_{st} gets a high value, displacements of the top and bottom layers at the contact surface are the same; from equation (41), the stiffness matrix of the stub in equation (54) gets the value 0: it means that the structure has no stubs and then the displacement is continuous at the interfaces between two layers. After assembling the element matrix and the element vector into the global matrix and global vector, the governing equation of motion of the plate is obtained aswhere are, respectively, the global mass matrix, global stiffness matrix, and global force vector and is the vector of unknown nodal displacements.
2.3. Static Bending and Free Vibration Problem
For the static bending problem, ignoring , equation (56) becomes
Solving equation (57), we get the vector of unknown nodal displacements .
For free vibration analysis, neglecting the effect of external force vector and assuming in equation (56), we get the following eigenvalue problem:
Solving equation (58), we get , which denotes the natural frequencies of the plate; for each value of , we have an eigenvector.
3. Numerical Results and Discussion
We now focus on numerically studying the static bending and free vibration responses of the twolayer composite plate. Each edge of the plate can be under the simply supported boundary condition or clamped boundary condition. For the simply supported edge,
For the clamped edge,
3.1. Static Bending Analysis of TwoLayer Composite Plates
3.1.1. Accuracy Studies
To verify the proposed method, we compare the deflections in two cases of materials: homogeneous and composite materials. Firstly, a square homogeneous plate under contribution load is investigated, and the length and the width are . The plate is simply supported at all edges (SSSS) or fully clamped (CCCC). In the published works [27, 28], the plate has one layer made of homogeneous material with , Poisson’s ratio , and the thickness h. But in this work, we assume that the plate has two layers connected by stubs with two cases—case 1: h_{n} = 0.999 h, h_{m} = 0.001 h, and k_{st} = 0 and case 2: h_{n} = h_{m} = 0.5 h and k_{st} = 10^{10}E. In both cases, material properties of two layers are set to be , thicknesses of the top and bottom layer are h_{n} and h_{m}, respectively, and h_{n} + h_{m} = h. The nondimensional deflection is defined aswhere D is given in equation (2).
Table 1 shows a comparison of the nondimensional deflections of SSSS and CCCC plates obtained by this method and other approaches (FEM [28] and analytical solution [28]). We meet a good agreement among three solutions.
We see that when the thickness of one layer is much higher than that of the other layer and the stiffness of the connector stub is very small, the behavior of the twolayer plate is the same as the behavior of the onelayer plate with the same thickness. In addition, when two layers have the same thickness and material, the stiffness of the connector stub gets a much higher value k_{st} = 10_{10}E, and then the twolayer plate is considered as the onelayer plate with the thickness equal to the sum of thickness of the two components.
In order to further confirm the accuracy of this method, we compare the nondimensional deflection of a square SSSS composite plate ( = 1 and the total thickness h = ) with two layers (0/90°) which have the same thickness . In this comparison, the plate is modeled as h_{n} = h_{m} = 0.5 h, k_{st} = 10^{10}E_{2}, and the fiber directions of the top and bottom layers are 0° and 90°, respectively. Table 2 shows the comparison of the nondimensional deflection of this composite plate subjected to uniformly distributed load and central concentrated load. The nondimensional deflection of the plate is calculated by the following formula:
It is shown plainly that the present results have a good agreement with the results given in [29] (analytical method). According to the above two comparisons, the present results have a good agreement with other published ones.
3.1.2. Numerical Results for Static Bending Analysis of TwoLayer Composite Plates
In this section, a rectangular twolayer composite plate subjected to the uniform load is investigated. The boundary conditions of the plate are SSSS (fully simply supported), CCCC (fully clamped), CSCS (two opposite edges are simply supported and two other edges are clamped), and CCSS (two adjacent edges are clamped and two other edges are simply supported). The sides, thickness, and material properties of the plate are set to be , , , , , , and . The nondimensional displacement of this plate is normalized by
(1) Influence of the LengthtoWidth Ratio a/b. In order for studying the effect of the aspect ratio on the mechanical bending behavior of the twolayer composite plate, this ratio is set to vary from 0.5 to 4, while , , and the ratio of stiffness of the connector stub are considered. The computed results are shown in Table 3 and Figures 4(a) and 4(b). According to Table 3 and Figure 4, it is seen that when the aspect ratio increases, the nondimensional deflection decreases. In addition, when the ratio , the nondimensional deflections of both boundary conditions CCCC and CSCS of plates are very close to each other.

(a)
(b)
(2) Influence of the ThicknesstoThickness Ratio h_{n}/h_{m}. Next, the research on the variation of deflection affected by the relationship between the bottom layer thickness and the top layer thickness is now investigated. The ratio increases from 1 to 4, , and the lengthtowidth ratio is set to be 1 and 2. The nondimensional deflection of the plate with different boundary conditions is shown in Table 4 and Figure 5. The deflection of the plate increases when increasing ratio ; this can be explained by the fact that Young’s modulus of the layer n is less than that of layer m, and the stiffness of the twolayer plate becomes smaller when the ratio increases.

(3) Influence of the Stiffness of the Connector Stub. Table 5 and Figure 6 show the nondimensional deflections as a function of the ratio of a square twolayer composite plate with and . The stiffness of the stub varies from 0 to 10^{6}E_{m}; when it changes from 0 to E_{m}, this plate is stiffer and the deflection decreases; but when the stiffness of the stub > E_{m}, then the effect of the stub goes to zero so that the twolayer plate connected by the stub becomes a twolayer one with no stub and the displacements at the contact place are the same; it means that there is no slip between two layers at the contact place.

(a)
(b)
(a)
(b)
To show clearly the slip between the top and bottom layers at the contact place, we depicted the ratio and in Figures 7 and 8, respectively, in which and are defined as
We can see again in Figure 7 that when the stiffness of the stub increases, the slip between two layers decreases and tends to zero. By looking at Figure 8 (the ratio = 0.1 is employed), it is very interesting to see that the minimum of slip between two layers at the contact place appears at the central plate with symmetric boundary conditions (CCCC, SSSS, and CSCS). The slip gets a maximum value at the simply supported edge, and this slip strongly depends on the boundary condition.
(4) Influence of Boundary Condition. We further study the deflection affected by the boundary condition. A square plate with = 12 and = 01 is considered. Figure 9 shows the displacement at x = and along the y direction in four boundary conditions. The boundary condition has strong effects on the deflection magnitude and deformation shape of the plate. It is shown that the plate with the simple support at all edges has the highest deflection and the one with the fully clamped support has the smallest deflection. The deformation shapes of the SSSS plate, CCCC plate, and CSCS plate are symmetric, while the deformation of the CCSS plate is asymmetric.
(a)
(b)
(c)
(d)
3.2. Free Vibration Analysis
3.2.1. Accuracy Studies
To confirm the accuracy of this work, we first present a comparison of frequencies of homogeneous CCCC and SSSS plates, with a = b = 1 m, h = and , E = 10920 MPa, and Poisson’s ratio ; similar to the bending analysis, h_{n} = h_{m} = 0.5 h and k_{st} = 10^{10}E are again used for this comparison. The nondimensional frequencies are defined by
The results are listed in Tables 6 and 7, in which Ferreira theory [28] was based on FEM and Mindlin theory [28] was based on the analytical solution. As expected, the results obtained by this work match well with those derived from other solutions.
Next, we compare the nondimensional frequencies of an SSSS composite plate 0/90°. The parameters of this plate are set to be , , , , , and . The nondimensional frequencies of the plate are normalized by
The comparison is shown in Table 8; note that the results in [29] are obtained from the exact solution. We can see that the present results and reference results are close to each other. From two comparisons, we conclude that the accuracy of this approach is accepted.
3.2.2. Numerical Results for Free Vibration Analysis of TwoLayer Composite Plates
In this section, we research the natural frequency of the twolayer composite plate utilizing the proposed method. The twolayer composite plate (length a, width b, thickness of the top layer h_{n}, and thickness of the bottom layer h_{m}) is studied herein. The nondimensional fundamental frequency is defined as(1) Influence of LengthtoWidth Ratio. Table 9 and Figure 10 show the dimensionless fundamental frequencies as a function of aspect ratio with and of twolayer composite plates altered by the boundary condition. It is evident that the highest and the smallest natural frequencies are found for CCCC and SSSS plates, respectively. When the aspect ratio increases, the nondimensional fundamental frequency increases, and the fundamental frequencies of CCCC and CSCS plates are close when the aspect ratio > 2.6; this is a similar behavior to the static bending problem.

(2) Influence of ThicknesstoThickness Ratio . To study the influence of thicknesstothickness ratio , a twolayer square composite plate with is used for this analysis. The nondimensional fundamental frequencies as a function of the ratio are tabulated in Table 10 and shown in Figures 11 and 12. The nondimensional fundamental frequency varies nonlinearly with the increasing ratio ; the minimum value of fundamental frequency will be obtained when the ratio is in the range of 1.8 to 2.6.

(3) Influence of Stiffness of the Connector Stub. In this exploration, the effect of stiffness of the connector stub on the nondimensional fundamental frequency of the twolayer composite plate is studied. A square twolayer plate with and varying from 0 to is investigated. Table 11 and Figure 13 report the computed results of this plate with = 1; we can see that increasing the ratio from 0 to 1 leads to the increase in the frequencies. When the ratio , the nondimensional fundamental frequency is almost unchanged.

Figures 14 and 15 plot nondimensional fundamental frequencies of the plate as a function of ratios , , and . In addition, the first six mode shapes of plates with boundary conditions such as SSSS, CCCC, CSCS, and CCSS are shown in Figures 16–19 with = 1, = 1, and = 10^{−1}. It is seen that the boundary condition has an effect on mode shapes as well as on the value of frequencies.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
4. Conclusions
In this work, new numerical results of mechanical behaviors and free vibration responses of the twolayer composite plate are explored, in which two layers are connected by stubs. We used the finite element method combined with the new modified firstorder shear deformation theory, which has the following advantages:(i)Simple formulations for theoretical representation(ii)No need for reduced integration or selective reduced integration scheme for the proposed method
The computed results of static bending and free vibration obtained by this approach are also compared to other solutions, showing a good agreement. We gained insight into the responses of deflections and natural frequencies. Besides, some geometrical and physical properties of this structure are also examined. Finally, from new numerical results, several conclusions may be achieved as follows:(i)The stiffness of the stub has a strong effect on static bending deflections and natural frequencies of the twolayer composite plate when < 1. When > 1, there is no slip between the top and the bottom layers at the contact place, and the static bending and free vibration responses of this plate are similar to those of the onelayer composite plate without stubs.(ii)Likewise, the boundary condition has great effect on the slip between two layers, the minimum slip appears at the center of the plates with the symmetric boundary condition, and the slip has the greatest value at the simply supported edge.(iii)When the ratio increases, the static bending deflection increases, but the fundamental frequency varies as a nonlinear function of this ratio.(iv)The higher the ratio , the smaller the deflection and the higher the fundamental frequency.
The new results of this work are useful for calculation, design, and testing, as well as for giving the optimal solution for the twolayer plate and shell in engineering and technologies. This study suggests some further works on buckling, dynamic response, and heat transfer problems of the twolayer composite plate using different plate theories.
Appendix
The matrix in equation (30) is defined aswhere for layer and for layer .
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.022018.30.
References
 P. Foraboschi, “Behavior and failure strength of laminated glass beams,” Journal of Engineering Mechanics, vol. 133, no. 12, pp. 1290–1301, 2007. View at: Publisher Site  Google Scholar
 P. Foraboschi, “Analytical solution of twolayer beam taking into account nonlinear interlayer slip,” Journal of Engineering Mechanics, vol. 135, no. 10, pp. 1129–1146, 2009. View at: Publisher Site  Google Scholar
 A. Kryžanowski, S. Schnabl, G. Turk, and I. Planinc, “Exact slipbuckling analysis of twolayer composite columns,” International Journal of Solids and Structures, vol. 46, no. 1415, pp. 2929–2938, 2009. View at: Publisher Site  Google Scholar
 P. L. Grognec, Q. H. Nguyen, and M. Hjiaj, “Exact buckling solution for twolayer Timoshenko beams with interlayer slip,” International Journal of Solids and Structures, vol. 49, no. 1, pp. 143–150, 2012. View at: Publisher Site  Google Scholar
 G. He and X. Yang, “Finite element analysis for buckling of twolayer composite beams using Reddy's higher order beam theory,” Finite Elements in Analysis and Design, vol. 83, pp. 49–57, 2014. View at: Publisher Site  Google Scholar
 G. He and X. Yang, “Dynamic analysis of twolayer composite beams with partial interaction using a higher order beam theory,” International Journal of Mechanical Sciences, vol. 90, pp. 102–112, 2015. View at: Publisher Site  Google Scholar
 G. He, D. Wang, and X. Yang, “Analytical solutions for free vibration and buckling of composite beams using a higher order beam theory,” Acta Mechanica Solida Sinica, vol. 29, no. 3, pp. 300–315, 2016. View at: Publisher Site  Google Scholar
 T. Hozjan, M. Saje, S. Srpčič, and I. Planinc, “Geometrically and materially nonlinear analysis of planar composite structures with an interlayer slip,” Computers & Structures, vol. 114115, pp. 1–17, 2013. View at: Publisher Site  Google Scholar
 P. Wu, D. Zhou, W. Liu, L. Wan, and D. Liu, “Elasticity solution of twolayer beam with a viscoelastic interlayer considering memory effect,” International Journal of Solids and Structures, vol. 9495, pp. 76–86, 2016. View at: Publisher Site  Google Scholar
 P. Wu, D. Zhou, and W. Liu, “2D elasticity solutions of twolayer composite beams with an arbitrarily shaped interface,” Applied Mathematical Modelling, vol. 40, no. 2, pp. 1477–1493, 2016. View at: Publisher Site  Google Scholar
 M. A. Uddin, A. H. Sheikh, T. Bennett, and B. Uy, “Large deformation analysis of two layered composite beams with partial shear interaction using a higher order beam theory,” International Journal of Mechanical Sciences, vol. 122, pp. 331–340, 2017. View at: Publisher Site  Google Scholar
 J. Wen, A. H. Sheikh, M. A. Uddin, and B. Uy, “Analytical model for flexural response of twolayered composite beams with interfacial shear slip using a higher order beam theory,” Composite Structures, vol. 184, pp. 789–799, 2018. View at: Publisher Site  Google Scholar
 T. V. Do, T. Q. Bui, T. T. Yu, D. T. Pham, and C. T. Nguyen, “Role of material combination and new results of mechanical behaviour for FG sandwich plates in thermal environment,” Journal of Computational Science, vol. 21, pp. 164–181, 2017. View at: Publisher Site  Google Scholar
 B. Cas, I. Planinc, and S. Schnabl, “Analytical solution of threedimensional twolayer composite beam with interlayer slips,” Engineering Structures, vol. 173, pp. 269–282, 2018. View at: Publisher Site  Google Scholar
 K. Ke, M. C. H. Yam, A. C. C. Lam, and K. F. Chung, “Local web buckling of singlecoped beam connections with slender web,” Journal of Constructional Steel Research, vol. 150, pp. 543–555, 2018. View at: Publisher Site  Google Scholar
 C. T. Latham, A. Toledano, H. Murakami, and F. Seible, “A sheardeformable twolayer plate element with interlayer slip,” International Journal for Numerical Methods in Engineering, vol. 26, no. 8, pp. 1769–1789, 1988. View at: Publisher Site  Google Scholar
 P. Foraboschi, “Analytical model for laminatedglass plate,” Composites Part B: Engineering, vol. 43, no. 5, pp. 2094–2106, 2012. View at: Publisher Site  Google Scholar
 D. Baraldi, A. Cecchi, and P. Foraboschi, “Broken tempered laminated glass: nonlinear discrete element modeling,” Composite Structures, vol. 140, pp. 278–295, 2016. View at: Publisher Site  Google Scholar
 P. Foraboschi, “Layered plate with discontinuous connection: exact mathematical model,” Composites Part B: Engineering, vol. 47, pp. 365–378, 2013. View at: Publisher Site  Google Scholar
 P. Foraboschi, “Threelayered sandwich plate: exact mathematical model,” Composites Part B: Engineering, vol. 45, no. 1, pp. 1601–1612, 2013. View at: Publisher Site  Google Scholar
 P. Foraboschi, “Threelayered plate: elasticity solution,” Composites Part B: Engineering, vol. 60, pp. 764–776, 2014. View at: Publisher Site  Google Scholar
 P. Wu, D. Zhou, W. Liu, W. Lu, and H. Fang, “3D exact solution of twolayer plate bonded by a viscoelastic interlayer with memory effect,” Composite Structures, vol. 164, pp. 291–303, 2017. View at: Publisher Site  Google Scholar
 P. Vidal, L. Gallimard, and O. Polit, “Robust layerwise C 0 finite element approach based on a variable separation method for the modeling of composite and sandwich plates,” Finite Elements in Analysis and Design, vol. 139, pp. 1–13, 2018. View at: Publisher Site  Google Scholar
 M. M. Alinia, G. Soltanieh, and M. Amani, “Inelastic buckling behavior of stocky plates under interactive shear and inplane bending,” ThinWalled Structures, vol. 55, pp. 76–84, 2012. View at: Publisher Site  Google Scholar
 I. Senjanović, N. Vladimir, and N. Hadžić, “Modified Mindlin plate theory and shear lockingfree finite element formulation,” Mechanics Research Communications, vol. 55, pp. 95–104, 2014. View at: Publisher Site  Google Scholar
 H.T. Thai and D.H. Choi, “A simple firstorder shear deformation theory for laminated composite plates,” Composite Structures, vol. 106, pp. 754–763, 2013. View at: Publisher Site  Google Scholar
 K. M. Liew, J. Wang, T. Y. Ng, and M. J. Tan, “Free vibration and buckling analyses of sheardeformable plates based on FSDT meshfree method,” Journal of Sound and Vibration, vol. 276, no. 3–5, pp. 997–1017, 2004. View at: Publisher Site  Google Scholar
 A. J. M. Ferreira, MATLAB Codes for Finite Element Analysis: Solids and Structures, Springer, Dordrecht, Netherlands, 2009, ISBN 9789400789555.
 J. N. Reddy, Mechanics of Laminated Composite Plates and Shells: Theory and Analysis, CRC Press, Boca Raton, FL, USA, 2003, ISBN 9780849315923.
Copyright
Copyright © 2019 Vu H. Nam 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.