Research Article  Open Access
Estimation of Local Delamination Buckling in Orthotropic Composite Plates Using Kirchhoff Plate Finite Elements
Abstract
We analyse the buckling process of composite plates with throughthewidth delamination and straight crack front applying uniaxial compression. We are focusing on the mixed mode buckling case, where the nonuniform distribution of the inplane forces controls the occurence of the buckling of the delaminated layers. For the analysis, semidiscrete finite elements will be derived based on the Lèvytype method. The method of harmonic balance is used for taking into account the force distribution that is generally non uniform inplane.
1. Introduction
Nowadays, fiberreinforced composite materials are applied everywhere in the industrial products, such as cars, air planes, ships, space vehicles, and sport equipment [1, 2]. The heterogeneous behaviour of these materials makes their damage analysis very important [3–6]. One of the main failure modes of layered composites is the delamination fracture, which can occur because of machining errors, or low velocity impact during the service life [7–16]. Buckling caused by uniaxial compression is one of the critical failures for these materials [17–24].
Many researchers have studied the field of delaminated composite plate buckling experimentally, analytically using finite element method (FEM). Chai et al. [25] have developed an analytical onedimensional model for simulating buckling of delaminated beam or plate specimens. The delamination was assumed to be close to the surface (thinfilm approach), and mixed buckling modes have been determined. On a beamplate model it was shown that the transverse shear effect reduces the critical buckling loads [26]. The effect of the boundary conditions was investigated on a 1D model, taking into account the large deformations, and the results were compared with the thinfilm approximation and experiments [27]. Anastasiadis and Simitses [28] presented a modified 1D model for improving the critical buckling loads. They used artificial springs along the crack line for preventing crack opening, which would result in inadmissible mode shapes [29]. Ovesy et al. [30] developed a layerwise theory based on the FirstOrder Shear Deformation Theory (FSDT) for analysing the postbuckling behaviour of multiple delaminated laminates. In order to prevent the inadmissible mode shapes they used contact constraints on the delaminated area. For the double delaminated beamplate model an analytical solution was introduced by Shu using a constrained model for the global buckling calculation [31]. Kim and Kedward [32] used the Classical Laminate Plate Theory (CLPT) for modelling a delaminated rectangular plate and analysed its behaviour with respect to the buckling using the Navier method. For the global stability analysis the delaminated region was treated as a reduced stiffness zone. The local stability analysis was carried out on a clamped plate assuming the same load along the local part as for the global model. The delamination buckling and delamination growth of crossply laminates were examined, induced by low velocity impact in [33]. The results of the optical measurement were compared with Xray nondestructive evaluation results. Wang and Lu [34] examined the failure mechanism of near the surface embedded delaminations under compression load using the energy method and experimental testing. Throughthewidth delaminated composite laminates subjected to compressive loads were investigated, using the RayleighRitz method and CLPT theory by Kharazi and Ovesy [35]. Experiments on multiple delaminated composite plates and a FEM study using ANSYS were carried out by Aslan and Şahin [36]. Kharazi developed a layerwise theory based on the FSDT where the prevention of the penetration of the delaminated layers was solved by using the penalty method and artificial springs [37]. The mixed mode buckling was captured with a contact model in the paper of Hwang and Liu [38]. A layerwise FEM solution using Abaqus and VICONOPT is given for composite plates with embedded rectangular delamination in Damghani et al. [39]. Marjanović and Vuksanović [40] used the Generalized Layerwise Plate Theory of Reddy and built a finite element (FE) model which was capable of handling more embedded delaminations in the composite structure. A 3D FEM model was developed for analysing the buckling behaviour of embedded and throughthewidth delaminated plates in [41].
For the global model an analytical solution was presented in Juhász and Szekrényes [42]. In that paper it was shown that the inplane force distribution greatly influences the buckling behaviour of the plate. In this paper a FE model is introduced, which concerns a plate with two simply supported opposite edges so the plate problem cannot be solved using beam theories. The other two boundary conditions (BCs) can be arbitrarily chosen. The plate is loaded with uniaxial compression. The critical inplane loads for the global buckling are calculated using a constrained model [29, 31]. Using this approach there is no need for using constraints along the delamination, which results in a very efficient calculation method contrary to the widely used contact models. The continuity between the delaminated and undelaminated parts is maintained using special transitional elements. This model is capable of predicting the global critical buckling loads and determines the corresponding inplane force distributions. Using the solution of the global model the local stability is analysed separately using a local FE model for each delaminated part. Because of the nonuniformity in the distribution of axial forces the stability analysis was carried out using the method of harmonic balance [43]. It was shown that these effects greatly influence the local buckling loads. In this paper the effects of different type of boundary conditions are analysed in a numerical example where the length of the delamination is varied, and the local and global critical amplitudes are calculated, which results in a stability map for the different BCs. It is shown that the different BCs greatly influence the buckling loads and the corresponding mixed mode buckling shapes. Using our method and the resulting stability maps the failure process of a compressed delaminated plate can be estimated, and it can be determined at which point will the plate buckle globally or locally, or in mixed mode case.
2. Model Creation
Let us consider a layered plate with several orthotropic plies and a closed throughthewidth delamination (see Figure 1). The plate length is , the length of the delamination is , and the thickness of the plate is . The plate is composed of two equivalent single layers (ESL) with thicknesses and , respectively.
The edges parallel to the axis are simply supported, and the constraint of the other two can be arbitrarily chosen. The plate is subjected to uniaxial compression in the midplane. We are considering thin structures; therefore, the solution is based on CLPT [44–47]:where and are constant parts in the inplane displacement functions of the undelaminated portions, which arise because of the kinematic coupling on the interface as Figure 2 shows. and are constant through the thickness in the top and bottom parts, where takes either “” or “,” respectively (see Figure 2). Moreover, is the deflection of the plate parts, respectively. Along the undelaminated parts we assume perfect adherence and no crack propagation at the crack tips. Therefore, the continuity of the undelaminated portions is ensured using the system of exact kinematic conditions (SEKC) [48]:And the inplane displacements in the global reference plane areFrom the literature it is well known that a free model allows the intersection of the delaminated portions into each other, which results in kinematically inadmissible mode shapes and wrong critical loads [29, 37, 49–51]. Avoiding this a constrained model is used, where the deflection of the delaminated top and bottom portions is common. Using these assumptions the strain fields of the undelaminated portions can be given asAnd for the delaminated portions the strain field will bewhere means the constant and means the linear part of the strain field in terms of the coordinate [52, 53]. For expressing the stresses of the composite laminate with the strain fields we have to use the constitutive equation of composite laminates (CEL) [52–55]:where the matrices , and can be calculated based on the literature.
Using these the strain energy density of the delaminated and undelaminated portions can be given as [17, 52, 56]where the inplane force and moment resultants are depending on the strain fields based on (6).
3. Finite Element Discretization
The finite element discretization concept can be seen in Figure 3.
As we have written in the previous section the inplane displacement fields are continuous on the undelaminated portions and independent of each other in the delaminated top and bottom portions (see Figure 3). Therefore, between the delaminated and undelaminated portions special transition elements were used on sections and which capture the crack tips. These elements ensure the continuity of the displacements [57]. Because of the opposite simply supported edges the Lèvytype method is applicable and the displacements fields in the direction can be given by the terms of Fourier series, HosseiniHashemi et al. [58], Bodaghi and Saidi [55], Thai and Kim [59], Szekrényes [48, 60], Thai and Kim [61], and Nguyen et al. [62]:where and , , are the amplitudes in the direction. As the displacement in the direction is given it is enough to discretize the model along , so we can write the nodal displacements aswhere is the nodal displacement vector and is a diagonal matrix containing the corresponding trigonometric terms from (8). This results in a semidiscrete finite element, which is capable of modelling platelike buckling. For the analysis we need to derive the material and the geometric stiffness matrices of the different elements of Figure 3.
Integrating the strain energy density the material stiffness matrix of an element can be derived based on the Hamilton principle [63–66]:where is the material or general stiffness matrix and is the element domain. Based on the literature the geometric stiffness matrix with respect to the axial compression can be given as [63, 67]As the inplane forces along the delaminated region will not be zero, for the local stability analysis we also need the geometric stiffness matrix with respect to the load:In ((11)(12)) the is the applied load on the element, is the vector of interpolation functions of the element, and is the derivative of the vector of interpolation functions.
3.1. Element of the Nondelaminated Parts
For the undelaminated parts an 8DoF element is used:where is the rotation of the crosssection in the node. The matrix of the trigonometric coefficients can be given asFor the inplane displacements linear interpolation functions were used, whereas for the transverse deflection a thirdorder function was applied:where , , and are constant coefficients and varies between and . Their value can be calculated using the CLPT from the following equations:where the derivation is carried out with respect to the dimensionless coordinate. By solving ((15)–(21)) the vector interpolation functions can be obtained [63]:where can be , , or .
Substituting the discretized displacement fields into (10) the strain energy can be given asCarrying out the integration over the element material stiffness matrix can be obtained.
Using the interpolation functions from (11) the geometric stiffness matrix can be derived.
3.2. Element of the Delaminated Parts
Because of the constrained model the transverse deflection is common, but the inplane displacements are independent in the delaminated portion, which results in a 12DoF element:The matrix can be composed based on the nodal displacement vector. For the inplane displacement the same linear interpolation functions were used, given by ((15)(16)), and for the transverse deflection the thirdorder function was applied in accordance with (17).
As the inplane displacements are independent, the potential energy has to be evaluated for both the top and bottom parts and the material stiffness matrix can be derived from the sum of the potential energies based on (23).
3.3. The Transition Elements
In accordance with Figure 3, the elements of sections and ensure the kinematic continuity between the delaminated and undelaminated portions. The vector of displacements of the element isAnd for the element denoted by we have
Based on the kinematic continuity the following 4 equations can be written based on the applied plate theory for the section:The equations take similar form for the section:Using the equations above and (20) and (21) the vector of interpolation functions can be obtained. The stiffness matrices can be calculated on the same way as it was shown before.
4. Stability Analysis
Based on Section 3 the structural matrices of the global model can be obtained. After applying the selected BCs on the and edges the critical loads and the corresponding eigenvectors can be calculated asThe corresponding global mode shapes and the resultant inplane force distributions can be obtained using the vector of interpolation functions by (22) and the CEL given by (6). Because of the inplane resultant forces at the crack tips the plate is able to buckle locally along the delamination. The local stability is analysed individually for the top and bottom delaminated plate portions, assuming plates with builtin end BCs along the crack tips. The local stability is affected by the distribution of the inplane forces (see Figure 4). For the local FE model we derived the elements of the individual top and bottom layers, using the same method as for the elements of the global model. The nonuniform resultant inplane forces of the global model are evaluated for every element at the middle. These values were normed with the value at the crack tip, and the element geometric stiffness matrices were multiplied with these values, taking into consideration the distribution along . Because of the simply supported edges the plate will have a half wave shape along the width, which results in the fact that the load along the crack tip will not be uniform (see Figure 4). Taking this aspect of the problem into consideration we applied the method of harmonic balance and wrote the Fourier series of the nodal displacements [43]:where are constant coefficients and is the vector of displacement values. Taking this back into (29) and applying some trigonometric identities we can obtain a system of equations in matrix form:The critical values and the corresponding mode shapes can be calculated from (31) and (29).
(a)
(b)
For validation purposes the model was solved using Abaqus. The plate is made by carbon/epoxy material using the following layup order: . Engineering constants of the layers are detailed in Table 1. The series expansion in (30) was carried out for two terms. Along the direction the plate was discretized using 14 elements, to capture the higher order mode shapes. The obtained critical values from (31) are ~40% higher than the loads of the problem with constant distribution along . The top ESL of the example in Section 6 was checked assuming constant force distribution along . The width of the plate was 100 mm, and the length of the plate was 105 mm. The S4R shell element was used for the analysis with 1 mm element size. The results show good agreement with the present calculations (see Table 2).


5. Boundary and Continuity Conditions
In this paper the process of loss of stability is determined by using a displacement controlled model based on Section 2. For solution, the Lèvytype method is used with the statespace approach [52]. From (7) using Hamilton’s principle, the governing PDEs of each section can be derived [52]. Applying (8) the obtained ODEs can be rearranged into the statespace model [52, 68]:where is the state vector. The general solution of (32) is [52, 68, 69]where is the vector of constants . At the crack tips we have to define 1010 continuity conditions (CCs) between the plate portions ①, ② and ②, ③. Because of the closed delamination (see Figure 1) the socalled Mujumdar conditions have to be used for fitting the moment and the Kirchhoff equivalent shear force [29]:where can take either, or respectively, and depends on and their derivatives. can take 1, 2, or 3 depending on the sections, which will be fit. In the BCs an axial displacement at has to be prescribed, and there is no other load. Substituting the solution of the statespace model into the BCs and CCs a system of inhomogeneous equations can be obtained:which can be solved for the constants. Using (33) we can get the displacement functions, and the inplane forces can be calculated using (6).
Using this model we calculated the arising forces at the edge of the plate and at the crack tips with respect to the axial displacement . The critical values of the global and local stability analysis were compared with these results.
5.1. Criterion of Constant Arc Length
All of the mode shapes were calculated with a maximum amplitude of 1 mm and scaled to fit the physical requirements. The amplitudes of the global and local modes were controlled using an arc length criterion [57]. This means that the arc length of the superimposed eigenshapes minus the axial displacement has to be equal to the length of the plate or the delamination:where is the scale factor for the mode shapes and is the buckled shape of the th buckling mode. For global mode shapes . For local mode shapes it is the signed sum of the axial displacements at the left and right crack tips. In case of mixed mode buckling iswhere is the axial displacement of the global model for the th mode scaled with the amplitude and is the static axial displacement.
5.2. Superposition of the Mode Shapes
Based on the linear relationship between the inplane normal force and the axial displacement of the displacement controlled model, critical axial displacements can be calculated for the critical loads. From these values the amplitudes are rising linearily:
We assume that the dominant inplane force distribution is determined by the first global mode. Therefore, local buckling in mixed mode case occurs only if the critical values of the local modes, which belong to the first global mode, are reached.
6. Numerical Example
In this section we adopt the method on a carbon/epoxy layered plate. The ply order of the plate is . The plate consists of 9 layers. The corresponding material data can be found in Table 1. The plate is symmetrically delaminated, and its geometric data is presented in Table 3. The stiffness matrices of each single layer were determined based on the elastic properties given by Table 1. The analysis was carried out with condition in (8). The order of the matrix in (31) was set to 2. The plate was discretized using 12 elements in all sections and 11 additional transitional elements were used at the crack tips. The position of the delamination was set above the 5th layer. At the edges and the same simply supported (SS) or builtin (BB) BCs were used. The length of the delamination was varied from 10 mm to 100 mm. The global critical forces with respect to the delamination length can be seen in Figure 5.

It can be seen that the obtained critical loads of the builtin plate are higher, but as the length of the delamination increases the effect of BCs gets less significant. The critical amplitudes of the local top and bottom delaminated portions can be seen in Figures 67.
As it can be seen the local critical values are higher in the simply supported cases. This is because different eigenshape belongs to the different BCs which results in different inplane force distribution. Again as the delamination length increases the effect of the BCs gets less significant. Using the displacement controlled model the critical axial displacements can be calculated for each critical amplitude. Based on this calculation stability diagrams can be obtained with respect to the axial displacement and the delamination length (see Figures 89).
On both pictures below the blue line the plate is stable. In the orange region the plate buckles globally, in the green region it buckles globally and the crack opens as the local top plate loses its stability, and above the green line the delaminated bottom portion buckles too. It has to be remarked that in the BB case the bottom part buckles only at higher axial compression; therefore, the green line is outside the range shown in Figure 9. The maximal critical amplitude was set to 2 mm. It can be seen that the builtin end plate is more stable, and its bottom part does not lose its stability up to the maximal axial displacement, whereas the simply supported plate loses its stability on smaller amplitudes. It can be noticed that as the delamination length increases the point of the global and local stability loss of the top plate gets close to each other. The presented critical loads are the first critical amplitudes. But if the plate is weak against uniaxial compression, higher order mode shapes are also feasible. These mode shapes can be superimposed using the arc length criterion. In the following we will show the process of stability loss of the simply supported and builtin end plates with 100 mm delamination length. The global critical amplitudes for the two types of BCs are listed in Table 4. For these values the critical axial compressions can be determined based on the displacement controlled model. The resulting forces with respect to the axial displacement for the simply supported case are shown in Figure 10. On the same way the critical axial displacements of the builtin end plate can be determined. Whereas the critical loads are higher than in case of simply supported BCs the critical axial displacements of the first 2 modes are smaller and only the third mode appears at higher displacement: 0.14 mm, 0.15 mm, and 0.33 mm. The maximal axial compression was chosen in both cases for the 120% of the third mode. The critical values of the delaminated portions were calculated for the local buckling case where the nonuniform distribution of the inplane forces does not count, but the calculated critical axial displacements were higher than the critical axial displacement of the first global mode; therefore, the plate loses its stability first globally.

From Figures 11 and 12 it can be seen that because of the different BCs different mode shapes appear. For the mixed mode buckling the local critical values were calculated for both cases using the nonuniform force distribution of the global modes. Here we present only the critical loads of the realizing local modes (see Table 5). As it can be seen only the first two local modes appear in both cases. The third mode would only appear at higher axial compression. At the builtin end case the local modes, calculated with the force distribution of the second global mode, are not present during the stability loss, because the critical values of these modes are much higher. The plate was also examined for the forces, but according to the results no stability loss occurs with respect to the forces at the crack tip. In accordance with Figures 8 and 9 the delaminated bottom part does not lose its stability at the selected maximal axial displacement.

The shapes of these modes were calculated with the inplane force distribution resulting from the corresponding global modes and were superimposed using the arc length criterion. Figures 13 and 14 show the buckled shapes at , , , and . On the superimposed shapes it can be seen that the dominant part of the solution is always the global and local first modes, but the higher order modes influence the shape slightly.
7. Conclusion
In this paper the buckling process of a delaminated layered plate was investigated. The formulation of the problem is based on the system of exact kinematic conditions (SEKC) by cutting the plate in the plane of the delamination and forming the continuity conditions. The problem was solved using FEM with selfdeveloped semidiscrete finite elements. The model contains special transitional elements which ensure the kinematic continuity between the delaminated and undelaminated portions. The delaminated region was modelled as a constrained section in the global model; therefore, there is no need for using contact along the delaminated area, which results in a calculation efficient and simple method for the estimation of the global critical buckling loads and the corresponding shapes. The local behaviour of the delaminated portion was analysed by a separate FE model. For the consideration of the nonuniform inplane force distribution the method of harmonic balance was used. On a numerical example the effects of the simply supported and builtin end BCs were determined with respect to the delamination length. It was shown that the BCs are influencing not only the critical loads but also the corresponding global mode shapes. Because of the different global mode shapes the local behaviour of the delaminated portions is different, as the inplane force distributions differ significantly. This results in the fact that whereas the simply supported plate buckles globally at lower values this configuration is more stable locally than the builtin end configuration. It was also shown that this effect is more significant if the delamination length is small. Stability diagrams with respect to the axial displacement and the delamination length were given where the global and mixed mode stability loss cases were shown. At one delamination length the process of stability loss was presented for both BCs. Here the effect of the BCs and the nonuniform inplane force distribution can be seen. This nonuniform distribution was not observed with respect to the different type of BCs in the literature, and we can state that it greatly alerts the buckled shape of the delaminated layers.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the Hungarian National Scientific Research Fund (OTKA) under Grant no. 4461506615 (108414).
References
 O. Gohardani and D. W. Hammond, “Ice adhesion to pristine and eroded polymer matrix composites reinforced with carbon nanotubes for potential usage on future aircraft,” Cold Regions Science and Technology, vol. 96, pp. 8–16, 2013. View at: Publisher Site  Google Scholar
 S. Giannis and K. Hansen, “Investigation on the joining of CFRPtoCFRP and CFRPtoaluminium for a small aircraft structural application,” in Proceedings—American Society for Composites: 25th Technical Conference of the American Society for Composites and 14th USJapan Conference on Composite Materials, 2022 September 2010, Dayton, Ohio, USA, J. B. Lantz, Ed., vol. 1, pp. 333–346, Curran Associates, Red Hook, NY, USA, 2011. View at: Google Scholar
 D. A. Hills, P. A. Kelly, D. N. Dai, and A. M. Korsunsky, Solution of Crack Problems: The Distributed Dislocation Technique, Kluwer Academic, Dordrecht, The Netherlands, 1996. View at: Publisher Site  MathSciNet
 D. F. Adams, L. A. Carlsson, and R. B. Pipes, Experimental Characterization of Advanced Composite Materials, CRC Press, Boca Raton, Fla, USA, 3rd edition, 2000.
 B. D. Davidson, F. O. Sediles, and K. D. Humphreys, “A sheartorsionbending test for mixedmode IIIIII delamination toughness determination,” in Proceedings of the 25th Technical Conference of the American Society for Composites and 14th USJapan Conference on Composite Materials, pp. 1001–1020, Dayton, Ohio, USA, September 2010. View at: Google Scholar
 M. F. S. F. De Moura, R. M. Guedes, and L. Nicolais, “Fracture: interlaminar,” in Wiley Encyclopedia of Composites, pp. 60–78, John Wiley & Sons, 2011. View at: Google Scholar
 L. N. Phillips, Ed., Design with Advanced Composite Materials, Springer, The Design Council, Berlin, Germany, 1989.
 V. Rizov, A. Shipsha, and D. Zenkert, “Indentation study of foam core sandwich composite panels,” Composite Structures, vol. 69, no. 1, pp. 95–102, 2005. View at: Publisher Site  Google Scholar
 V. I. Rizov, “Nonlinear indentation behavior of foam core sandwich composite materials—a 2D approach,” Computational Materials Science, vol. 35, no. 2, pp. 107–115, 2006. View at: Publisher Site  Google Scholar
 A. D. Zammit, S. Feih, and A. C. Orifici, “2D numerical investigation of pretension on low velocity impact damage of sandwich structures,” in Proceedings of the 18th International Conference on Composite Materials (ICCM18 '11), pp. 1–6, Jeju International Convention Center, Jeju, Republic of Korea, August 2011. View at: Google Scholar
 R. A. Chaudhuri and K. Balaraman, “A novel method for fabrication of fiber reinforced plastic laminated plates,” Composite Structures, vol. 77, no. 2, pp. 160–170, 2007. View at: Publisher Site  Google Scholar
 N. Carrere, T. Vandellos, and E. Martin, “Multilevel analysis of delamination initiated near the edges of composite structures,” in Proceedings of the 17th International Conference on Composite Materials (ICCM '09), pp. 1–10, Edinburgh, UK, July 2009. View at: Google Scholar
 V. N. Burlayenko and T. Sadowski, “A numerical study of the dynamic response of sandwich plates initially damaged by lowvelocity impact,” Computational Materials Science, vol. 52, no. 1, pp. 212–216, 2012. View at: Publisher Site  Google Scholar
 J. Rhymer, H. Kim, and D. Roach, “The damage resistance of quasiisotropic carbon/epoxy composite tape laminates impacted by high velocity ice,” Composites Part A: Applied Science and Manufacturing, vol. 43, no. 7, pp. 1134–1144, 2012. View at: Publisher Site  Google Scholar
 G. Goodmiller and S. TerMaath, “Investigation of composite patch performance under lowvelocity impact loading,” in Proceedings of the 55th AIAA/ASME/ASCE/AHS/SC Structures, Structural Dynamics, and Materials Conference, National Harbor, Md, USA, 2014. View at: Google Scholar
 C. Elanchezhian, B. V. Ramnath, and J. Hemalatha, “Mechanical behaviour of glass and carbon fibre reinforced composites at varying strain rates and temperatures,” Procedia Materials Science, vol. 6, pp. 1405–1418, 2014, Proceedings of the 3rd International Conference on Materials Processing and Characterisation (ICMPC '14). View at: Publisher Site  Google Scholar
 R. Guo and A. Chattopadhyay, “Development of a finiteelementbased design sensitivity analysis for buckling and postbuckling of composite plates,” Mathematical Problems in Engineering, vol. 1, no. 3, pp. 255–274, 1995. View at: Publisher Site  Google Scholar
 L. P. Kollár, “Buckling of rectangular composite plates with restrained edges subjected to axial loads,” Journal of Reinforced Plastics and Composites, vol. 33, no. 23, pp. 2174–2182, 2014. View at: Publisher Site  Google Scholar
 G. Tarján, Á. Sapkás, and L. P. Kollár, “Stability analysis of long composite plates with restrained edges subjected to shear and linearly varying loads,” Journal of Reinforced Plastics and Composites, vol. 29, no. 9, pp. 1386–1398, 2010. View at: Publisher Site  Google Scholar
 H.T. Thai and D.H. Choi, “Analytical solutions of refined plate theory for bending, buckling and vibration analyses of thick plates,” Applied Mathematical Modelling, vol. 37, no. 1819, pp. 8310–8323, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 H.T. Thai, M. Park, and D.H. Choi, “A simple refined theory for bending, buckling, and vibration of thick plates resting on elastic foundation,” International Journal of Mechanical Sciences, vol. 73, pp. 40–52, 2013. View at: Publisher Site  Google Scholar
 C. Klobedanz, A study of the effect of delamination size on the critical sublaminate buckling load in a composite plate using the Ritz method [Ph.D. thesis], Rensselaer Polytechnic Institute, Troy, NY, USA, 2014.
 S. A. M. Ghannadpour, H. R. Ovesy, and E. ZiaDehkordi, “Buckling and postbuckling behaviour of moderately thick plates using an exact finite strip,” Computers & Structures, vol. 147, pp. 172–180, 2015. View at: Publisher Site  Google Scholar
 H. R. Ovesy, A. Totounferoush, and S. A. M. Ghannadpour, “Dynamic buckling analysis of delaminated composite plates using semianalytical finite strip method,” Journal of Sound and Vibration, vol. 343, pp. 131–143, 2015. View at: Publisher Site  Google Scholar
 H. Chai, C. D. Babcock, and W. G. Knauss, “One dimensional modelling of failure in laminated plates by delamination buckling,” International Journal of Solids and Structures, vol. 17, no. 11, pp. 1069–1083, 1981. View at: Publisher Site  Google Scholar
 G. A. Kardomateas and D. W. Schmueser, “Buckling and postbuckling of delaminated composites under compressive loads including transverse shear effects,” AIAA Journal, vol. 26, no. 3, pp. 337–343, 1988. View at: Publisher Site  Google Scholar
 G. A. Kardomateas, “Large deformation effects in the postbuckling behavior of composites with thin delaminations,” AIAA Journal, vol. 27, no. 5, pp. 624–631, 1989. View at: Publisher Site  Google Scholar
 J. S. Anastasiadis and G. J. Simitses, “Spring simulated delamination of axiallyloaded flat laminates,” Composite Structures, vol. 17, no. 1, pp. 67–85, 1991. View at: Publisher Site  Google Scholar
 P. M. Mujumdar and S. Suryanarayan, “Flexural vibrations of beams with delaminations,” Journal of Sound and Vibration, vol. 125, no. 3, pp. 441–461, 1988. View at: Publisher Site  Google Scholar
 H. R. Ovesy, M. A. Mooneghi, and M. Kharazi, “Postbuckling analysis of delaminated composite laminates with multiple throughthewidth delaminations using a novel layerwise theory,” ThinWalled Structures, vol. 94, pp. 98–106, 2015. View at: Publisher Site  Google Scholar
 D. Shu, “Buckling of multiple delaminated beams,” International Journal of Solids and Structures, vol. 35, no. 13, pp. 1451–1465, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. Kim and K. T. Kedward, “A method for modeling the local and global buckling of delaminated composite plates,” Composite Structures, vol. 44, no. 1, pp. 43–53, 1999. View at: Publisher Site  Google Scholar
 J. T. Ruan, F. Aymerich, J. W. Tong, and Z. Y. Wang, “Optical evaluation on delamination buckling of composite laminate with impact damage,” Advances in Materials Science and Engineering, vol. 2014, Article ID 390965, 9 pages, 2014. View at: Publisher Site  Google Scholar
 X. Wang and G. Lu, “Local buckling of composite laminar plates with various delaminated shapes,” ThinWalled Structures, vol. 41, no. 6, pp. 493–506, 2003. View at: Publisher Site  Google Scholar
 M. Kharazi and H. R. Ovesy, “Postbuckling behavior of composite plates with throughthewidth delaminations,” ThinWalled Structures, vol. 46, no. 7–9, pp. 939–946, 2008. View at: Publisher Site  Google Scholar
 Z. Aslan and M. Şahin, “Buckling behavior and compressive failure of composite laminates containing multiple large delaminations,” Composite Structures, vol. 89, no. 3, pp. 382–390, 2009. View at: Publisher Site  Google Scholar
 M. Kharazi, H. R. Ovesy, and M. Asghari Mooneghi, “Buckling analysis of delaminated composite plates using a novel layerwise theory,” ThinWalled Structures, vol. 74, pp. 246–254, 2014. View at: Publisher Site  Google Scholar
 S.F. Hwang and G.H. Liu, “Buckling behavior of composite laminates with multiple delaminations under uniaxial compression,” Composite Structures, vol. 53, no. 2, pp. 235–243, 2001. View at: Publisher Site  Google Scholar
 M. Damghani, D. Kennedy, and C. Featherston, “Global buckling of composite plates containing rectangular delaminations using exact stiffness analysis and smearing method,” Computers & Structures, vol. 134, pp. 32–47, 2014. View at: Publisher Site  Google Scholar
 M. Marjanović and D. Vuksanović, “Layerwise solution of free vibrations and buckling of laminated composite and sandwich plates with embedded delaminations,” Composite Structures, vol. 108, no. 1, pp. 9–20, 2014. View at: Publisher Site  Google Scholar
 J. D. Whitcomb, “Mechanics of instabilityrelated delamination growth,” in Composite Materials: Testing and Design, vol. 9, pp. 215–230, ASTM, 1990. View at: Google Scholar
 Z. Juhász and A. Szekrényes, “Progressive buckling of a simply supported delaminated orthotropic rectangular composite plate,” International Journal of Solids and Structures, 2015. View at: Publisher Site  Google Scholar
 W. W. Bolotin, Kinetische Stabilität Elastischer Systeme, VEB Deutscher Verlag der Wissenschaften, Berlin, Germany, 1961.
 A. Szekrényes, “Analysis of classical and firstorder shear deformable cracked orthotropic plates,” Journal of Composite Materials, vol. 48, no. 12, pp. 1441–1457, 2014. View at: Publisher Site  Google Scholar
 L. S. Ma and T. J. Wang, “Relationships between axisymmetric bending and buckling solutions of FGM circular plates based on thirdorder plate theory and classical plate theory,” International Journal of Solids and Structures, vol. 41, no. 1, pp. 85–101, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Amabili and S. Farhadi, “Shear deformable versus classical theories for nonlinear vibrations of rectangular isotropic and laminated composite plates,” Journal of Sound and Vibration, vol. 320, no. 3, pp. 649–667, 2009. View at: Publisher Site  Google Scholar
 A. M. Zenkour, “Exact mixedclassical solutions for the bending analysis of shear deformable rectangular plates,” Applied Mathematical Modelling, vol. 27, no. 7, pp. 515–534, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Szekrényes, “The system of exact kinematic conditions and application to delaminated firstorder shear deformable composite plates,” International Journal of Mechanical Sciences, vol. 77, pp. 17–29, 2013. View at: Publisher Site  Google Scholar
 C. N. Della and D. Shu, “Vibration of delaminated multilayer beams,” Composites Part B: Engineering, vol. 37, no. 23, pp. 227–236, 2006. View at: Publisher Site  Google Scholar
 Y. Guo, M. Ruess, and Z. Gürdal, “A contact extended isogeometric layerwise approach for the buckling analysis of delaminated composites,” Composite Structures, vol. 116, pp. 55–66, 2014. View at: Publisher Site  Google Scholar
 J. Wang and L. Tong, “A study of the vibration of delaminated beams using a nonlinear antiinterpenetration constraint model,” Composite Structures, vol. 57, no. 1–4, pp. 483–488, 2002. View at: Publisher Site  Google Scholar
 J. N. Reddy, Mechanics of Laminated Composite Plates and Shells—Theory and Analysis, CRC Press, Boca Raton, Fla, USA, 2004.
 L. Kollár and G. Springer, Mechanics of Composite Structures, Cambridge University Press, Cambridge, UK, 2002.
 J. Ye, Laminated Composite Plates and Shells—3D modelling, Springer, London, UK, 2003.
 M. Bodaghi and A. R. Saidi, “Lévytype solution for buckling analysis of thick functionally graded rectangular plates based on the higherorder shear deformation plate theory,” Applied Mathematical Modelling, vol. 34, no. 11, pp. 3659–3673, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 S. W. Tsai, Theory of Composites Design, Think Composites, Dayton, Ohio, USA, 1992.
 A. Szekrényes, “A special case of parametrically excited systems: free vibration of delaminated composite beams,” European Journal of Mechanics—A/Solids, vol. 49, pp. 82–105, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 S. HosseiniHashemi, M. Fadaee, and H. Rokni Damavandi Taher, “Exact solutions for free flexural vibration of Lévytype rectangular thick plates via thirdorder shear deformation,” Applied Mathematical Modelling, vol. 35, no. 2, pp. 708–727, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 H.T. Thai and S.E. Kim, “Levytype solution for buckling analysis of orthotropic plates based on two variable refined plate theory,” Composite Structures, vol. 93, no. 7, pp. 1738–1746, 2011. View at: Publisher Site  Google Scholar
 A. Szekrényes, “Application of Reddy's thirdorder theory to delaminated orthotropic composite plates,” European Journal of Mechanics A: Solids, vol. 43, pp. 9–24, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 H.T. Thai and S.E. Kim, “Lévytype solution for free vibration analysis of orthotropic plates based on two variable refined plate theory,” Applied Mathematical Modelling, vol. 36, no. 8, pp. 3870–3882, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Q.H. Nguyen, E. Martinelli, and M. Hjiaj, “Derivation of the exact stiffness matrix for a twolayer Timoshenko beam element with partial interaction,” Engineering Structures, vol. 33, no. 2, pp. 298–307, 2011. View at: Publisher Site  Google Scholar
 K.J. Bathe, Finite Element Procedures, Prentice Hall, Upper Saddle River, NJ, USA, 1996.
 M. Petyt, Introduction to Finite Element Vibration Analysis, Cambridge University Press, Cambridge, UK, 2nd edition, 2010. View at: Publisher Site  MathSciNet
 E. Ventsel and T. Krauthammer, Thin Plates and Shells—Theory, Analysis and Applications, Marcel Dekker, New York, NY, USA, 2001. View at: Publisher Site
 T. Özben and N. Arslan, “FEM analysis of laminated composite plate with rectangular hole and various elastic modulus under transverse loads,” Applied Mathematical Modelling, vol. 34, no. 7, pp. 1746–1762, 2010. View at: Publisher Site  Google Scholar
 R. Szilard, Theories and Applications of Plate Analysis, John Wiley & Sons, Hoboken, NJ, USA, 2004. View at: Publisher Site
 W. Q. Chen, Y. F. Wu, and R. Q. Xu, “State space formulation for composite beamcolumns with partial interaction,” Composites Science and Technology, vol. 67, no. 1112, pp. 2500–2512, 2007. View at: Publisher Site  Google Scholar
 K. Xu, A. K. Noor, and Y. Y. Tang, “Threedimensional solutions for coupled thermoelectroelastic response of multilayered plates,” Computer Methods in Applied Mechanics and Engineering, vol. 126, no. 34, pp. 355–371, 1995. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Zoltán Juhász and András Szekrényes. 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.