Abstract
A new collocation methodology is presented to predict failure and progressive damage behavior of composite plates in this paper. The present work deals with composite plates containing initial geometric imperfections and different boundary conditions under uniaxial inplane compressive load. In the present study, the domain is discretized with LegendreGaussLobatto nodes and the approximation of displacement fields is performed by Legendre Basis Functions (LBFs). The onset of damage and damage evolution are predicted by Hashin’s failure criteria and by proposed material degradation models. Three geometric degradation models are also assumed to estimate the degradation zone around the failure location which are named complete, region and node degradation models.
Keywords
Progressive damage; Hashin failure criterion; Collocation; Composite plate; LegendreGaussLobatto nodes; Initial imperfection
1 INTRODUCTION
Over the last decade thinwalled structures have extensively used in various lightweight structural components. If these structures be subjected to inplane compressive loads, they are prone to buckling. Composite plates are one of these structures that can be optimized to carry inplane compressive loads beyond their buckling resistance. They combine high strength with low weight, which makes them ideal for use in many industries such as aerospace engineering. Several laminate theories have been exploited to explain composite laminated plates' behavior. The Classical Laminated Plate Theory (CLPT) can be considered as an extension of classical plate theory based on the Kirchhoff hypothesis for isotropic plates and can be applied, if thickness of the laminate is small and the effects of transverse shear stresses are neglected ( Dong et al. 1962 Dong, SB, KS Pister, and RL Taylor, (1962), On the theory of laminated anisotrophic plates and shells. Journal of Aeronautical Science 29(29):969975. ; Yang et al.1966 Yang, P Constance, Charles H Norris, and Yehuda Stavsky, (1966), Elastic wave propagation in heterogeneous plates. International Journal of solids and structures 2(4):665684. ; Ambart͡sumi͡an 1970 Ambart͡sumi͡an, Sergeĭ Aleksandrovich. (1970), Theory of anisotropic plates: strength, stability, vibration: Technomic Publishing Company (Lancaster). ). An improvement on the CLPT leads to creation of the First order Shear Deformation Theory (FSDT) which accounts for the transverse shear effects ( Whitney 1969 Whitney, JM, (1969), The effect of transverse shear deformation on the bending of laminated plates. Journal of Composite Materials 3(3):534547. ; Pagano 1970 Pagano, NJ, (1970), Shear Deformation in Heterogeneous Anisotropic Plates1. J. appl. Mech 37:10311036. ) in relatively thick laminates.
For moderately thick composite plates with geometric imperfection, linear, nonlinear and postbuckling analyses without taking account of damage and also material degradations are very conservative analyses and give uncertain estimates of critical loads. Hence, damage analyses and ultimate strength of such structures have been of considerable research interest. In the field of buckling and postbuckling analyses, many researchers have extensively investigated these behaviors for composite beams and plates without considering the damage effects. Turvey and Marshall (1995) Turvey, Geoffrey J, and Ian H Marshall. (1995), Buckling and postbuckling of composite plates: Springer Science & Business Media (Berlin). and Argyris and Tenek (1997) Argyris, John, and Lazarus Tenek, (1997), Recent advances in computational thermostructural analysis of composite plates and shells with strong nonlinearities. Applied Mechanics Reviews 50:285306. published excellent reviews on past studies of buckling and postbuckling of structures using different methods and recently, Pagani and Carrera (2017) Pagani, Alfonso, and Erasmo Carrera, (2017), Largedeflection and postbuckling analyses of laminated composite beams by Carrera Unified Formulation. Composite Structures 170:4052. have analyzed the large deflection and postbuckling behaviors of laminated beams by Carrera unified formulation.
Finite Element Method (FEM) is one the most common methods for analyzing the buckling and postbuckling of composite laminated plates. This method is based on dividing the plate’s domain into a finite number of simple subdomains called elements. The behavior of each element is described by some specified functions. Finally, the solution of the whole plate as an assembly of its elements can be obtained according to the procedures applicable to the standard discrete problems such as those described by Zienkiewicz (1977) Zienkiewicz, Olgierd Cecil, and Robert Leroy Taylor.(1977), The finite element method. Volume 3: McGrawhill (London). and Bathe (2006) Bathe, KlausJürgen. (2006), Finite element procedures: KlausJurgen Bathe. . Although FEM is an intensely powerful technique, it has some disadvantages. It needs large number of elements to reach acceptable stress results and in aspect of computational costs, it takes much longer computing time and requires larger amount of core storage.
Finite Strip Method (FSM) is another universally applicable method for buckling and postbuckling analyses of plates and plate structures. It can be considered as a particular kind of simplified finite element method in which a special element called strip is used. Finite strip method is based on discretization of the domain into longitudinal strips and interpolates the behavior in the longitudinal direction by different functions and in the transverse direction by polynomial functions. Cheung (1976) Cheung, YK, (1976), Finite strip method in structural analysis. Pergamon press (Oxford):26. may be considered as the pioneer who first proposed the concept of FSM. Cheung established FSM for the analysis of simply supported plates. The studies presented by Smith and Sridharan (1978) Smith, TR Graves, and S Sridharan, (1978), A finite strip method for the buckling of plate structures under arbitrary loading. International Journal of Mechanical Sciences 20(10):685693. proposed FSM for buckling of isotropic plate under edge loading. Recently, Ovesy, Ghannadpour and their coworkers ( Ovesy et al. 2005 Ovesy, HR, SAM Ghannadpour, and G Morada, (2005), Geometric nonlinear analysis of composite laminated plates with initial imperfection under end shortening, using two versions of finite strip method. Composite structures 71(3):307314. ) have made a contribution by introducing two different versions of finite strip methods, namely the fullenergy and semienergy finite strip approaches. Two other different versions of finite strip method, namely spline and semianalytical methods are also developed by them for predicting the response of rectangular laminates with nonsymmetric and symmetric forms of initial imperfection. They used both formulations to predict the nonlinear response of channel sections when subjected to uniform endshortening in their plane ( Ovesy et al. 2006 Ovesy, HR, J Loughlan, and SAM Ghannadpour, (2006), Geometric nonlinear analysis of channel sections under end shortening, using different versions of the finite strip method. Computers & structures 84(13):855872. ). An exact finite strip is introduced by Ghannadpour and Ovesy (2008) Ghannadpour, SAM, and HR Ovesy, (2008), An exact finite strip for the calculation of relative postbuckling stiffness of Isection struts. International Journal of Mechanical Sciences 50(9):13541364. to investigate the exact relative postbuckling stiffness of Isection struts. To extend their works, they developed a high accuracy finite strip for the buckling and postbuckling analyses of moderately thick symmetric crossply composite plates based on FSDT ( Ovesy et al. 2016 Ovesy, HR, E ZiaDehkordi, and SAM Ghannadpour, (2016), High accuracy postbuckling analysis of moderately thick composite plates using an exact finite strip. Computers & Structures 174:104112. ).
Furthermore, many other methods employed to investigate composite plate’s behavior subjected to compressive loads such as new class of numerical methods called meshless methods ( Liu, 2009 Liu, GuiRong. (2009), Meshfree methods: moving beyond the finite element method: Taylor & Francis (London). ). These methods have attracted the attention of many researchers in recent years due to the fact that meshless methods do not require a mesh to discretize the problem domain as in the finite element method and require only a scattered set of nodes to model the domain of interest. Meshless or meshfree methods have been proposed in multiple varieties such as Generalized Finite Difference method (GFD) ( Liszka, 1984 Liszka, Tadeusz, (1984), An interpolation method for an irregular net of nodes. International Journal for Numerical Methods in Engineering 20(9):15991612. ), which partial differential equations have been solved by applying irregular grids or clouds of points, the Smooth Particle Hydrodynamics (SPH) ( Monaghan, 1988 Monaghan, Joseph J, (1988), An introduction to SPH. Computer physics communications 48(1):8996. ), a fully Lagrangian meshless method which at first was used for simulating astrophysical phenomena, the Diffuse Element Method (DEM) ( Nayroles et al., 1992 Nayroles, B, G Touzot, and P Villon, (1992), Generalizing the finite element method: diffuse approximation and diffuse elements. Computational mechanics 10(5):307318. ), which is the first meshless method employed moving leastsquares approximation, the ElementFree Galerkin method (EFG) ( Belytschko et al., 1996 Belytschko, Ted, et al., (1996), Meshless methods: an overview and recent developments. Computer methods in applied mechanics and engineering 139(14):347. ) and Reproducing Kernel Particle Method (RKPM) ( Liu et al., 1995 Liu, Wing Kam, Sukky Jun, and Yi Fei Zhang, (1995), Reproducing kernel particle methods. International journal for numerical methods in fluids 20(8‐9):10811106. ). A corrected collocation method to apply essential boundary conditions in the meshless methods have been developed by Wagner and Liu (2000) Wagner, G. J., & Liu, W. K. (2000). Application of essential boundary conditions in mesh‐free methods: a corrected collocation method. International Journal for Numerical Methods in Engineering, 47(8), 13671379. .
In the field of studies that investigate the plates behavior using meshfree methods, it can be pointed out to the research conducted by Krysl and Belytschko (1995) Krysl, Petr, and Ted Belytschko, (1995), Analysis of thin plates by the elementfree Galerkin method. Computational Mechanics 17(1):2635. . They have implemented EFG to analyze the thin plates bending behavior. Lin and Jen (2005) Lin, C. H., & Jen, M. H. R. (2005). Analysis of a laminated anisotropic plate by Chebyshev collocation method. Composites Part B: Engineering, 36(2), 155169. proposed Chebyshev collocation method to solve the governing differential equations of a laminated anisotropic plate. A high order collocation method have been developed by Ferreira et al. (2009) Ferreira, A. J. M., Castro, L. M., & Bertoluzza, S. (2009). A high order collocation method for the static and vibration analysis of composite plates using a firstorder theory. Composite Structures, 89(3), 424432. for the static and vibration analyses of composite plates. Buckling of laminated composite plates subjected to various mechanical and thermal loads using meshless collocations have been studies by Singh et al. (2013) Singh, S., Singh, J., & Shukla, K. K. (2013). Buckling of laminated composite plates subjected to mechanical and thermal loads using meshless collocations. Journal of Mechanical Science and Technology, 27(2): 327336. . Free vibrations of beams were investigated by Carrera et al. (2013) Carrera, Erasmo, et al., (2013), Free vibrations of beams accounting for refined theories and radial basis functions. In 9th international symposium on vibrations of continuous systems. Italy. using Radial Basis Functions (RBFs). They have analyzed the free vibration characteristics of beams on the basis of unified formulation. Liew and Huang (2003) Liew, KM, and YQ Huang, (2003), Bending and buckling of thick symmetric rectangular laminates using the moving leastsquares differential quadrature method. International Journal of Mechanical Sciences 45(1):95114. and Liew et al. (2006) Liew, KM, et al., (2006), Postbuckling analysis of laminated composite plates using the meshfree kpRitz method. Computer methods in applied mechanics and engineering 195(7):551570. studied buckling and postbuckling of laminates using the moving leastsquares differential quadrature and meshfree kpRitz method. Also, Liew et al. (2011c) Liew, Kim Meow, Xin Zhao, and Antonio JM Ferreira, (2011c), A review of meshless methods for laminated and functionally graded plates and shells. Composite Structures 93(8):20312041. published a review paper for meshless methods and their applications in the buckling analysis of laminated or functionally graded plates. Recently, Ghannadpour and Barekati (2016) Ghannadpour, SAM, and M Barekati, (2016), Initial imperfection effects on postbuckling response of laminated plates under endshortening strain using Chebyshev techniques. ThinWalled Structures 106:484494. investigated the postbuckling behavior of composite plates using Chebyshev techniques by considering initial imperfection effects. More recently, Ghannadpour et al. (2017) Ghannadpour, S. A. M., Kiani, P., & Reddy, J. N. (2017). Pseudo spectral method in nonlinear analysis of relatively thick imperfect laminated plates under endshortening strain. Composite Structures, 182, 694710. have developed a high accuracy meshless analysis as a collocation method with Legendre Basis Functions (LBF) for nonlinear analysis of thin and moderately thick composite plates.
As emphasized earlier, all above mentioned research works, investigate the linear and nonlinear behaviors of plates without considering the effects of damage or failure and therefore the obtained results may be unrealistic. On the other hand, as was observed, all mentioned studies have been developed by different numerical or semianalytical methods. However, strength analyses of metal plates and damage analyses of composite plates are usually carried out using FEM in which the storage required is extremely large, and the computational time is too lengthy. Only a limited number of such analyses have been implemented by other numerical or semianalytical methods.
Several simplified semianalytical methods have been purposed by Brubak and Hellesland (2007a Brubak, Lars, and Jostein Hellesland, (2007a), Approximate buckling strength analysis of arbitrarily stiffened, stepped plates. Engineering Structures 29(9):23212333. , b Brubak, Lars, and Jostein Hellesland, (2007b), Semianalytical postbuckling and strength analysis of arbitrarily stiffened plates in local and global bending. ThinWalled Structures 45(6):620633. , ^{2008} Brubak, Lars, and Jostein Hellesland, (2008), Strength criteria in semianalytical, large deflection analysis of stiffened plates in local and global bending. ThinWalled Structures 46(12):13821390. , ^{2011} Brubak, Lars, and Jostein Hellesland, (2011), Semianalytical postbuckling analysis of stiffened imperfect plates with a free or stiffened edge. Computers & Structures 89(17):15741585. ) and Brubak et al. (2007c) Brubak, Lars, Jostein Hellesland, and Eivind Steen, (2007c), Semianalytical buckling strength analysis of plates with arbitrary stiffener arrangements. Journal of Constructional Steel Research 63(4):532543. for investigating the strength of stiffened and unstiffened metal plates with imperfection and cutout. Accelerated analysis techniques were proposed by Orifici et al. (2008) Orifici, AC, et al., (2008), Degradation investigation in a postbuckling composite stiffened fuselage panel. Composite Structures 82(2):217224. using FEM tools such as ABAQUS to model the degradation and fracture mechanics of composite stiffened fuselage panel in a postbuckling regime. Some experimental works with comparison by FEM analysis can be found on ultimate strength of composite plates such as study investigated by Hayman et al. (2011) Hayman, B, et al., (2011), Studies of the buckling of composite plates in compression. Ships and Offshore Structures 6(12):8192. . These studies are included a parametric study of ultimate strength analysis of rectangular simply supported composite plates with different geometrical imperfections. Yang and Hayman (2015a Yang, Qiao Jie, and Brian Hayman, (2015a), Prediction of postbuckling and ultimate compressive strength of composite plates by semianalytical methods. Engineering Structures 84:4253. , b Yang, Qiao Jie, and Brian Hayman, (2015b), Simplified ultimate strength analysis of compressed composite plates with linear material degradation. Composites Part B: Engineering 69:1321. ) and Yang et al. (2013) Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352. recently carried out a valuable studies established simplified methods to predict ultimate strength of composite plates with different initial imperfection values under inplane compressive loads.
In the present paper, a simplified progressive damage methodology is developed to predict the ultimate strength of imperfect composite plates with different boundary conditions under compression. The simplified methodology is referred to the assumptions of small deflection theory by which the computational time is greatly reduced. A new methodology is presented based on collocation method in which the interested domain is discretized with LegendreGaussLobatto nodes and therefore this approach do not require a mesh to discretize the problem domain. However, with less computational efforts, an acceptable field of stress can be obtained which is necessary for damage analyses. The formulations are based on the concept of the principle of minimum potential energy and the approximation of displacement fields is performed by Legendre Basis Functions (LBFs). In this new collocation approach which is based on the total potential energy, there is no need to enforce the natural boundary conditions, and the essential boundary conditions are satisfied very easily. Therefore, it is easy to apply different boundary conditions. The structural model is based on the first order shear deformation theory. The onset of damage and damage evolution are predicted by Hashin’s failure criteria ( Hashin and Rotem, 1973 Hashin, Zvi, and Assa Rotem, (1973), A fatigue failure criterion for fiber reinforced materials. Journal of composite materials 7(4):448464. ) and by proposed material degradation models. Three types of degradation models are assumed to estimate the degradation zone around the failure location which are named complete, region and node degradation models. With these assumed models, the accurate results can be extracted by reducing the area of failure around a node, and this is done very easily and quickly due to the proposed formulations. The properties of damaged materials instantaneously reduced to a 1% of primary value. Some examples involving various boundary conditions, initial imperfections and different thickness to width ratios are investigated to demonstrate the validity and capability of the proposed method. The accuracy of the present work is examined by comparing the numerical results with the previous studies.
2 Theoretical formulations
A typical rectangular laminated plate of dimensions of
To described the deformation of a laminated plate, displacement vector
Whose components of displacement vector are given by equation (2) .
Where
Due to moderately large displacements assumption, the vector of Green's strains
Where
where
The mechanical state is described by the stress vector
On the assumption that the plate is in a state of plane stress
The plate internal actions are including, namely, the inplane stress resultants
In which
Where
Where
By using equations (4) , (9) and (10) the total potential energy
3 Displacement fields and boundary conditions
In order to approximate the displacement functions of the problem, it is necessary to describe the specified boundary conditions. The laminates under consideration have two different types of boundary conditions as represented in Figure 2 . As shown in this figure, for both types of boundary conditions, the inplane displacement in xdirection at
However, as it can be seen, the labelling schemes are included in these figures to assign the related boundary conditions. The letter S refers to simply supported boundary condition and the letter C refers to clamped boundary condition on the specified edge. Therefore, the two aforementioned boundary conditions are different from the viewpoint of outofplane boundary conditions. The details of edges conditions for both Type A and Type B are expressed in Table 1 .
It is noted that the words “Free” and “Held” denote whether the edge is free to move or restricted against any movement, respectively. Also, the word “straight” denotes the edge can move longitudinally (for
Since the Rayleigh–Ritz technique is used in this study, approximated displacement fields have to satisfy only the aforementioned essential boundary conditions. The approximation of displacement fields is performed by Legendre basis functions
Also Legendre polynomials satisfy the threeterm recursion expressed by equation (13) .
Where
Where
The socalled boundary function
Where
Specifically, as noted earlier, to the aim of the present work Legendre orthogonal polynomials have been chosen. It is worth mentioning that different types of polynomials can be used, whose characteristic features have been discussed in the literature.
As can be seen in righthand side of equation (14) , each displacement field
4 Equilibrium equations and solution procedure
As mentioned before, the equilibrium equations are obtained based on the concept of the principle of minimum potential energy. Therefore, the elements in equation (11) should be written as described in previous section. For this purpose, the vectors
where
Accordingly, the strain vectors defined in equation (4) with the assumption of small deflection theory can also be rewritten as equations (20) and (21) .
Where
And finally by using equations (20) , the total potential energy
Since in the current study, the domain of the plates is discretized by a set of nodes therefore the above continuous integrals should be replaced by summations where they can be calculated over all nodes. For this purpose, LegendreGaussLobatto nodes are established here and can be obtained by solving equation (23).
Where the parameters m and n denote the number of nodes in both
In order to achieve better accuracy and also to avoid excessive number of nodes to reduce computational costs, an appropriate weight coefficient can be considered for each node. Calculation of the weight coefficients for nodes is performed by taking idea from GaussLobatto rules and therefore, the continuous integral of total potential energy (i.e. equation (22) ) after eliminating the constant factors is then converted to the following relation:
Where
To obtain the equilibrium equations of the problem using the principle of minimum potential energy, the discretized form of the total potential energy equation (24) should be minimized with respect to the unknown primary variables
5 Progressive damage model methodology
In this section, the methodology of progressive damage analysis including the failure criteria, material degradation model and ply geometric degradation models is described in details.
5.1 Failure criteria
In order to determine the failure load and the corresponding failure mode, a proper failure criterion must be established. The selected failure criteria applied in the present study are proposed by Hashin and Rotem (1973) Hashin, Zvi, and Assa Rotem, (1973), A fatigue failure criterion for fiber reinforced materials. Journal of composite materials 7(4):448464. .
The Hashin failure criteria used herein include four different damage functions which correspond to the different modes of failure namely fiber tension, fiber compression, matrix tension, and matrix compression. These criteria can be presented as equation (28) .
Fiber failure in tension
Fiber failure in compression
Matrix failure in tension
Matrix failure in compression
Here
5.2 Material degradation model
As it is known, when damage occurs in a composite structure, the effective material properties are reduced. This reduction can be modeled in this study by the following matrix.
Where
According to degradation model presented by Hayman et al. (2011) Hayman, B, et al., (2011), Studies of the buckling of composite plates in compression. Ships and Offshore Structures 6(12):8192. and also in order to better comparison with the results obtained by Yang et al. (2013) Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352. , the transverse shear stiffness matrix
In the current progressive damage model, the reduction of material properties is considered to happen instantaneously. When failure is detected in specific zone (complete, region or node), its properties are instantaneously reduced to %1 of their initial and undamaged values (i.e.
5.3 Ply geometric degradation models
In this study, three geometric degradation models are assumed to estimate the degradation zone around the failure location which are named complete, region and node degradation models. When failure occurs in a location, the material properties of its zone should be changed.

Complete Degradation Model (CDM) in which, the material properties of the entire ply are reduced. (as illustrated in Figure 4 (a))

Region Degradation Model (RDM) in which a plate is divided into 9 regions and material properties of the region where failure has occurred are reduced. (as illustrated in Figure 4 (b))

Node Degradation Model (NDM) in which, the material properties of the area around the failed node are reduced. (as illustrated in Figure 1 (c))
6 Numerical results and discussion
In this study, in order to implement the proposed formulations for analyzing the progressive damage of composite plates, a computer program is developed based on Fortran 77 software package. To reduce the execution time on multicore processors, a parallel programming technique is used using Open MultiProcessing (OpenMP) interface. The program is executed on a computer with TYAN FT48B8812 mainboard, 4 AMD CPU by 2.20 GHz frequency (4×16 cores) and 128.00 GB RAM.
The mentioned proposed formulations should be verified through a number of comparisons and this is done by comparing the results with those obtained by Hayman et al. (2011) Hayman, B, et al., (2011), Studies of the buckling of composite plates in compression. Ships and Offshore Structures 6(12):8192. and Yang et al. (2013) Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352. . For this purpose, the plates which are considered in this study are square plates (
The composite plates under consideration are modeled as
Since in this study, one of the geometric degradation models is region degradation model (RDM) in which a plate is divided into 9 regions, so specifying the size of each region is necessary. Therefore, in order to better comparison between the results, the size of the regions has been chosen as those considered by Yang et al. (2013) Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352. : regions dimensions 1, 3, 7 and 9 are 160×160mm. Regions 2 and 8 are 180×160mm each. Regions 4 and 6 are 160×180mm and region 5 is 180×180mm.
To do the convergence analysis, two composite plates with total thickness of 16 and 32 mm are selected. The boundary conditions are imposed as described for Type A and the plates have the maximum initial imperfection of 3% and 2%. Figure 5 shows the convergence study for both plates with the number of terms
convergence study of 16 mm and 32 mm composite plates with 3% and 2% initial imperfection respectively
As can be seen, for the plates under consideration the convergence studies with regard to the number of terms have revealed that 8 terms are sufficient to obtain converged results. However, the number of 9 terms is used to ensure accurate convergence in all analyses. Thus, the total number of unknown coefficients is 407. Also, similar analyzes were conducted with regards to the number of nodes for all three geometric degradation models and it was concluded that the total number of 169 nodes (m×n=13×13) are sufficient to obtain converged results.
6.1 Results for CDM.
To investigate the results for complete degradation model (CDM), two types of boundary conditions are considered as shown in Figure 2 . The results for composite plates with geometrical and mechanical properties mentioned in previous section and for CDM are given in Tables 3
and 4 for different maximum initial imperfections and total thicknesses. Other given data in these tables are first ply failure (FPF) stress, number and angle of that layer and the location failure, last ply failure (LPF) stress and also number and angle of last failed ply. Total number of failed plies is also given in these tables. Ultimate strength of composite plates with boundary condition Type A is tabulated in Table 3 and they are compared by results reported by Yang et al. (2013) Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352. and the results computed for boundary condition Type B are presented in Table 4 . Also, the loadcenter outofplane displacement history
Results for progressive damage of composite plates using CDM and for boundary condition Type A
Results for progressive damage of composite plates using CDM and for boundary condition Type B
As it can be seen, according to the results presented in Tables 3 and 4 for laminated plates with different boundary conditions, almost the ply fails first. However, for plates with boundary condition Type A, first failure usually occurs in the uppermost layer while for the plates having boundary condition Type B this mostly occurs in the lowermost ply. The only exceptions for both types are the two thicker plates with small initial imperfection (0.1%). For these special cases, the first ply failure occurred in ply with fiber orientation angle. Another important point is the location of first failure for both composite plates with different boundary conditions. For composite plates with all simplysupported edges (Type A), first failure usually occurs in the middle of the plates while for plates having one clamped edge (Type B) damage correspond to the first ply failure load mostly initiates from the midpoint of the clamped edge.
It is also seen that the ultimate strength of the composite plate is usually attained at the first incidence of fiber failure especially for plates having one clamped edge. Furthermore, the investigations show that for thin plates with small initial imperfections, there are little reserve strength beyond the first ply failure load while the ultimate strength of the thicker plates or the plates with larger initial imperfections have higher values than first ply failure loads. In these laminates, all plies usually have to experience matrix failure.
It is also seen from both tables and figures that by increasing the maximum value of initial imperfection of the plates, both FPF and LPF loads are decreased. Furthermore, as it can be observed from tables and figures, ultimate strength of composite plates having at least one clamped edge, is usually much greater than those plates with boundary condition Type A.
6.2 Results for RDM
In this section, the results for region degradation model (RDM) are presented and investigated. Tables 5 and 6 show the results for composite plates having boundary conditions Type A and Type B, respectively. Similar to the section of CDM results, first ply failure stresses and locations have been achieved. Also in both tables, last ply failure stresses and the number of matrix and fiber failed regions are included for each laminate thickness and initial imperfection. Furthermore, the values of ultimate strength obtained by CDM ( Tables 3 and 4 ) are also retabulated here for better comparison.
Results for progressive damage of composite plates using RDM and for boundary condition Type A
Results for progressive damage of composite plates using RDM and for boundary condition Type B
It can be observed that there is no difference between the first ply failure loads obtained by both CDM and RDM. However, the CDM model gives generally slightly lower ultimate loads than the RDM model. This difference is slightly higher for plates with boundary condition Type B.
In Figures 10 and 11 , the applied load is plotted against the central outofplane displacement and end shortening displacement respectively, for plates having total thickness of 32 mm and 2% initial imperfection. In order to better comparison, the results for both types of boundary conditions A and B are also represented in these figures.
Response of Load vs maximum outofplane displacement for 32 mm composite plates with 2% initial imperfection (Type A and Type B)
Response of Load vs end shortening for 32 mm composite plates with 2% initial imperfection (Type A and Type B)
It is seen that the results obtained for plates with boundary condition Type B have 1020% upper ultimate loads with respect to the results obtained for Type A. This is due to the fact that the boundary conditions have significant effects on the response of imperfect composite plates and as mentioned earlier the last ply failure loads of plates having one clamped edge are greater than the corresponding values in plates with all simplysupported edges.
However, regarding the first ply failure, for plates with small initial imperfection, Type A gives lower values than the results of Type B but by increasing amount of initial imperfection, FPF loads take higher values than Type B. In general, the laminates with boundary condition Type B can carry higher loads after their first ply fails.
Moreover, considering the location of first ply failure, plates with boundary conditions Type A fail in the center of the plate as mentioned in previous section, while laminates with boundary condition Type B fail in near the center of clamped edge.
6.3 Results of NDM.
The results for node degradation model (NDM) are presented and discussed here. Only the results correspond to the plates with boundary conditions Type A are presented in this section. These results are tabulated in Table 7 . As for previous models, first and last ply failure loads are reported in this table. In addition, the number of failed nodes in each fiber or matrix mode are also given. In order to better discussions, the results obtained by previous models ( Tables 3 and 5 ) and those reported by Hayman et al. (2011) Hayman, B, et al., (2011), Studies of the buckling of composite plates in compression. Ships and Offshore Structures 6(12):8192. are tabulated in the separate columns.
Results for progressive damage of composite plates using NDM and for boundary condition Type A
The behaviors of load against central outofplane displacement and load against endshortening displacement for plates with and 32 mm and for two different values of initial imperfection are depicted in Figures 12  15 . In these figures, in addition to the results obtained by NDM, the results taken from two other models, CDM and RDM, are also represented.
Response of load vs maximum outofplane displacement, compression between degradation model for 16 mm composite plate with %1 initial imperfection (Type A)
Response of load vs end shortening, compression between degradation model for 16 mm composite plate with %1 initial imperfection (Type A)
Response of load vs maximum outofplane displacement, compression between degradation model for 32 mm composite plate with %3 initial imperfection (Type A)
Response of load vs end shortening, compression between degradation model for 32 mm composite plate with %3 initial imperfection (Type A)
As it can be observed, there is no difference between the values of first ply failure loads calculated by three degradation models however, the complete model gives generally slightly lower values of ultimate load than the two others.
With comparison between three degradation models in Table 7 it is seen that the node degradation model gives upper last ply failure load and this can be due to the effects of dividing the plate (or ply) into smaller areas. Therefore, if one wants to obtain better and more accurate results, the node degradation model should be used in which more computation time is also needed but in the early stages of structural design, the results taken from region model can be accepted with slightly lower accuracy.
However, as it is seen, the ultimate loads predicted by all geometric degradation models in this study and those reported by Yang et al. (2013) Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352. are still smaller than those calculated by FEM analysis. This is mainly due to neglecting of the nonlinear terms in straindisplacement relations and therefore postbuckling effects, which are particularly considerable for thin plates. But for thick plates even by neglecting of those effects, the results have less deviation from FEM.
7 Conclusions
A new methodology has been introduced in this study to investigate the progressive damage analysis of composite plates with initial imperfection under inplane compressive load. The concept of the first order shear deformation theory and the assumption of small deflection have been established to drive the equilibrium equations. In presented method, the domain of the plate is discretized with LegendreGaussLobatto nodes. The onset of damage has been predicted by Hashin’s failure criteria and material properties of damaged zone have been degraded by instantaneous material degradation model. Three geometric degradation models were assumed to estimate the degradation zone around the failure location. The results were compared with linear analysis reported by Yang et al. (2013) Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352. and also nonlinear FE analysis presented by Hayman et al (2011) Hayman, B, et al., (2011), Studies of the buckling of composite plates in compression. Ships and Offshore Structures 6(12):8192. . By comparison, it has been concluded that the obtained results by increasing the thickness of the plates and decreasing the area around the location of failure can take better accuracy.
References
 Ambart͡sumi͡an, Sergeĭ Aleksandrovich. (1970), Theory of anisotropic plates: strength, stability, vibration: Technomic Publishing Company (Lancaster).
 Argyris, John, and Lazarus Tenek, (1997), Recent advances in computational thermostructural analysis of composite plates and shells with strong nonlinearities. Applied Mechanics Reviews 50:285306.
 Bathe, KlausJürgen. (2006), Finite element procedures: KlausJurgen Bathe.
 Belytschko, Ted, et al., (1996), Meshless methods: an overview and recent developments. Computer methods in applied mechanics and engineering 139(14):347.
 Brubak, Lars, and Jostein Hellesland, (2007a), Approximate buckling strength analysis of arbitrarily stiffened, stepped plates. Engineering Structures 29(9):23212333.
 Brubak, Lars, and Jostein Hellesland, (2007b), Semianalytical postbuckling and strength analysis of arbitrarily stiffened plates in local and global bending. ThinWalled Structures 45(6):620633.
 Brubak, Lars, and Jostein Hellesland, (2008), Strength criteria in semianalytical, large deflection analysis of stiffened plates in local and global bending. ThinWalled Structures 46(12):13821390.
 Brubak, Lars, and Jostein Hellesland, (2011), Semianalytical postbuckling analysis of stiffened imperfect plates with a free or stiffened edge. Computers & Structures 89(17):15741585.
 Brubak, Lars, Jostein Hellesland, and Eivind Steen, (2007c), Semianalytical buckling strength analysis of plates with arbitrary stiffener arrangements. Journal of Constructional Steel Research 63(4):532543.
 Carrera, Erasmo, et al., (2013), Free vibrations of beams accounting for refined theories and radial basis functions. In 9th international symposium on vibrations of continuous systems. Italy.
 Cheung, YK, (1976), Finite strip method in structural analysis. Pergamon press (Oxford):26.
 Dong, SB, KS Pister, and RL Taylor, (1962), On the theory of laminated anisotrophic plates and shells. Journal of Aeronautical Science 29(29):969975.
 Ferreira, A. J. M., Castro, L. M., & Bertoluzza, S. (2009). A high order collocation method for the static and vibration analysis of composite plates using a firstorder theory. Composite Structures, 89(3), 424432.
 Ghannadpour, SAM, and M Barekati, (2016), Initial imperfection effects on postbuckling response of laminated plates under endshortening strain using Chebyshev techniques. ThinWalled Structures 106:484494.
 Ghannadpour, SAM, and HR Ovesy, (2008), An exact finite strip for the calculation of relative postbuckling stiffness of Isection struts. International Journal of Mechanical Sciences 50(9):13541364.
 Ghannadpour, S. A. M., Kiani, P., & Reddy, J. N. (2017). Pseudo spectral method in nonlinear analysis of relatively thick imperfect laminated plates under endshortening strain. Composite Structures, 182, 694710.
 Hashin, Zvi, and Assa Rotem, (1973), A fatigue failure criterion for fiber reinforced materials. Journal of composite materials 7(4):448464.
 Hayman, B, et al., (2011), Studies of the buckling of composite plates in compression. Ships and Offshore Structures 6(12):8192.
 Krysl, Petr, and Ted Belytschko, (1995), Analysis of thin plates by the elementfree Galerkin method. Computational Mechanics 17(1):2635.
 Liew, Kim Meow, Xin Zhao, and Antonio JM Ferreira, (2011c), A review of meshless methods for laminated and functionally graded plates and shells. Composite Structures 93(8):20312041.
 Liew, KM, and YQ Huang, (2003), Bending and buckling of thick symmetric rectangular laminates using the moving leastsquares differential quadrature method. International Journal of Mechanical Sciences 45(1):95114.
 Liew, KM, et al., (2006), Postbuckling analysis of laminated composite plates using the meshfree kpRitz method. Computer methods in applied mechanics and engineering 195(7):551570.
 Liszka, Tadeusz, (1984), An interpolation method for an irregular net of nodes. International Journal for Numerical Methods in Engineering 20(9):15991612.
 Lin, C. H., & Jen, M. H. R. (2005). Analysis of a laminated anisotropic plate by Chebyshev collocation method. Composites Part B: Engineering, 36(2), 155169.
 Liu, GuiRong. (2009), Meshfree methods: moving beyond the finite element method: Taylor & Francis (London).
 Liu, Wing Kam, Sukky Jun, and Yi Fei Zhang, (1995), Reproducing kernel particle methods. International journal for numerical methods in fluids 20(8‐9):10811106.
 Monaghan, Joseph J, (1988), An introduction to SPH. Computer physics communications 48(1):8996.
 Nayroles, B, G Touzot, and P Villon, (1992), Generalizing the finite element method: diffuse approximation and diffuse elements. Computational mechanics 10(5):307318.
 Orifici, AC, et al., (2008), Degradation investigation in a postbuckling composite stiffened fuselage panel. Composite Structures 82(2):217224.
 Ovesy, HR, SAM Ghannadpour, and G Morada, (2005), Geometric nonlinear analysis of composite laminated plates with initial imperfection under end shortening, using two versions of finite strip method. Composite structures 71(3):307314.
 Ovesy, HR, J Loughlan, and SAM Ghannadpour, (2006), Geometric nonlinear analysis of channel sections under end shortening, using different versions of the finite strip method. Computers & structures 84(13):855872.
 Ovesy, HR, E ZiaDehkordi, and SAM Ghannadpour, (2016), High accuracy postbuckling analysis of moderately thick composite plates using an exact finite strip. Computers & Structures 174:104112.
 Pagani, Alfonso, and Erasmo Carrera, (2017), Largedeflection and postbuckling analyses of laminated composite beams by Carrera Unified Formulation. Composite Structures 170:4052.
 Pagano, NJ, (1970), Shear Deformation in Heterogeneous Anisotropic Plates1. J. appl. Mech 37:10311036.
 Reddy, J.N., (2004), Mechanics of laminated composite plates and shells: theory and analysis, CRC press (Boca Raton).
 Singh, S., Singh, J., & Shukla, K. K. (2013). Buckling of laminated composite plates subjected to mechanical and thermal loads using meshless collocations. Journal of Mechanical Science and Technology, 27(2): 327336.
 Smith, TR Graves, and S Sridharan, (1978), A finite strip method for the buckling of plate structures under arbitrary loading. International Journal of Mechanical Sciences 20(10):685693.
 Turvey, Geoffrey J, and Ian H Marshall. (1995), Buckling and postbuckling of composite plates: Springer Science & Business Media (Berlin).
 Wagner, G. J., & Liu, W. K. (2000). Application of essential boundary conditions in mesh‐free methods: a corrected collocation method. International Journal for Numerical Methods in Engineering, 47(8), 13671379.
 Whitney, JM, (1969), The effect of transverse shear deformation on the bending of laminated plates. Journal of Composite Materials 3(3):534547.
 Yang, P Constance, Charles H Norris, and Yehuda Stavsky, (1966), Elastic wave propagation in heterogeneous plates. International Journal of solids and structures 2(4):665684.
 Yang, Qiao Jie, and Brian Hayman, (2015a), Prediction of postbuckling and ultimate compressive strength of composite plates by semianalytical methods. Engineering Structures 84:4253.
 Yang, Qiao Jie, and Brian Hayman, (2015b), Simplified ultimate strength analysis of compressed composite plates with linear material degradation. Composites Part B: Engineering 69:1321.
 Yang, Qiao Jie, Brian Hayman, and Harald Osnes, (2013), Simplified buckling and ultimate strength analysis of composite plates in compression. Composites Part B: Engineering 54:343352.
 Zienkiewicz, Olgierd Cecil, and Robert Leroy Taylor.(1977), The finite element method. Volume 3: McGrawhill (London).
Publication Dates

Publication in this collection
04 June 2018 
Date of issue
2018
History

Received
13 July 2017 
Reviewed
29 Nov 2017 
Accepted
29 Nov 2017