Prediction of Progressive Failure Behaviour of Composite Skewed Hypar Shells Using Finite Element Method

Progressive failure behaviour of laminated composite skewed hypar shells is investigated in this paper through finite element approach. For composite materials which are weak in transverse shear, failure may initiate at any inner lamina or interface and may remain undetected and unattended to. Such latent damages may progress gradually and lead to a sudden total ply failure. An iterative method is adopted using the various failure theories to predict the first ply failure load. After the first ply failure, the stiffness of the failed element is totally discarded from the laminate and the remaining laminate is considered for further analysis. The developed finite element code is validated through solution of a benchmark problem. Numerical experiments are carried out to obtain the first ply and progressive failure of simply supported hypar shells under transverse distributed load. While obtaining the failure loads, the shell is considered to be under plane stress condition and in linear elastic range. Results are studied meticulously to extract a set of conclusions of practical significance regarding the failure characteristics of composite hypar shells.


Introduction
In modern engineering and technology, a composite is defined as a material composed of stacked layers with fiber reinforcement held in proper location with the rigid binder.A number of laminae are bonded together to obtain a "laminate" which stands for a structural unit of a beam, plate, or shell.Nowadays researchers are emphasizing more on laminated composite shells realizing the strength and stiffness potentials of this advanced material.Shell structures are preferred in the industry to cover large column free open spaces as are found in airports, shopping malls, auditoriums, and car parking lots.A skewed hypar shell is aesthetically appealing and being doubly ruled is easy to cast.Moreover, this configuration can allow entry of north light and due to this advantage it finds use as roofing units in practical civil engineering applications.Researchers like Sahoo and Chakravorty [1,2] and Nayak and Bandyopadhyay [3] studied bending and vibration characteristics of graphite-epoxy hypar shells.
It is important to note that failure study of the laminated hypar shell is needed for their confident application in the industry.First and ultimate ply failure of composite plates have received attention from researchers like Singh and Kumar [4], Akhras and Li [5], and Ganesan and Liu [6] who reported that failure of the laminated composites is progressive in nature and the load at which failure initiates (first ply failure load) in the laminates is much lower than its ultimate load carrying capacity.First ply failure of laminated plates and shells also reported by Reddy and Pandey [7], Kam and Jan [8], Kam et al. [9] Prusty et al. [10], and Satish Kumar and Srivastava [11].The progressive failure analysis of laminated composite plates under transverse static loading was carried out in linear and elastic range by Pal and Ray [12].Philippidis and Antoniou [13] computed a progressive failure analysis model for glass/epoxy composite giving an extensive comparison between numerically calculated stress-strain response up to failure and experimental data.Very recently, Ellul et al. [14] studied progressive failure analysis of fiber reinforced composite plates subjected to out-of-plane bending.Tolson and Zabaras [15] developed two-dimensional finite element analysis for determining progressive failure in composite plates using seven degrees of freedom.A first ply and progressive failure analysis of laminated unstiffened and stiffened composite panels under static loadings were investigated by Prusty [16] using the geometric and material linearity.
It is noted from the literature review that the failure aspects of industrially important hypar shell is totally missing.To fulfill this lacuna the present investigation intends to study the failure progress of hypar shell by a linear finite element approach from practical engineering point of view.

Mathematical Formulation
A laminated composite hypar shell (Figure 1) of uniform thickness ℎ and twist radius of curvature   is considered.Keeping the total thickness the same, the thickness may consist of any number of thin laminae each of which may be arbitrarily oriented at an angle  with reference to the  axis of the coordinate system.The surface equation of this shell is The radius of cross curvature may be evaluated by differentiating the surface equation of shell and for shallow shells which are taken up for the present study the same may be expressed as 1/  =  2 / = 4/.An eight-noded curved quardratic isoparametric element (Figure 2) having five degrees of freedom at each node which includes the in-plane, transverse displacements, and slopes along the plan directions.The governing equations and the systematic development of the stiffness matrix of the shell have been reported in Sahoo and Chakravorty [2].
The strain-displacement relation is given by where In-plane strain components for a lamina situated at a distance "" from the lamina mid-plane are evaluated in global axes as Lamina strains are transformed from the global axes of the shell to the local axes of the lamina using the following transformation matrix: where  = sin ,  = cos . ( Lamina stresses are obtained using the following constitutive relation of the lamina: where , and  66 =  12 .
The basic problem of static takes the following form: where [] is the overall stiffness matrix and {} and {} are generalized displacement and load vectors, respectively.The above equation is solved by Gauss elimination technique and from the global nodal displacement vector {} thus obtained, the element displacement vectors {  } are known.Using {  } in (( 2)) the strains can be evaluated at the Gauss points, which give the generalized force and moment resultants at the Gauss points.The stresses are obtained using the stressstrain constitutive relationship.Lamina stresses are used in well-accepted failure theories like maximum stress, Tsai-Hill, Tsai-Wu, and Hoffman failure criteria, shown in (( 8)) to (( 11)), to evaluate the first ply and ultimate ply failure loads of the composite hypar shells under present study.
Maximum Stress Failure Criteria.According to maximum stress theory, the failure initiates if at least one of the following criteria is satisfied: where  1 = Normal stress along 1st direction, that is, along the direction of lamina.
2 = Normal stress along 2nd direction, that is, along the perpendicular direction of lamina.
and  = Strength of the material in failure along 1st and 2nd direction, respectively.
= Shear strength of the material.
Tsai-Hill Criteria.According to this theory, the first ply failure occurs when the following criteria are satisfied in any one of the laminae: where  1 and  2 = Normal stress along 1st and 2nd direction, respectively.
and  = Strength of the material in failure along 1st and 2nd direction, respectively.
= Shear strength of the material.
Tsai-Wu Criteria.According to this theory the first ply failure occurs when the following criteria are satisfied in any one of the laminae: where Hoffman's Criteria.According to this theory, the lamina failure is assumed to occur when the following criteria are satisfied: where   and   = Tensile strength along 1st and 2nd direction, respectively.  and   = Compressive strength along 1st and 2nd direction, respectively. 1 and  2 = Normal stress along 1st and 2nd direction, respectively. 12 = Shear stress in the 1-2 plane.
= Shear strength of the material.
After the first ply failure (FPF), the stiffness of the failed element is totally discarded from the laminate and the remaining shell is considered for further analysis.The failed element is replaced with an equivalent element of the degraded material properties.Displacements and stresses are recalculated and the stresses for the remaining shell portions are verified against the selected failure criteria to compute the next failure load and failure point that would fail immediately under the increased share of stresses.The progressive failure study for the shell is carried out using this iterative process, which ultimately achieves the ultimate ply failure (UPF) load.This iterative process is repeated at each load step to investigate the failure of the elements one by one.The ultimate ply failure occurs when it is found that such iterations do not converge.The schematic algorithm of progressive failure is shown Figure 3.

Numerical Problems
First, ply failure loads obtained from the present formulation are compared with the linear failure loads reported by Kam et al. [9] in Table 1.The results match closely.Also, the results displayed in Table 2 show a good match with the results obtained by Reddy [17] using exact method.This confirms the correct incorporation of hypar shell curvature in the present formulation.For the authors' own problems geometric properties of the skewed hypar shells under simply supported condition and different stacking sequences are presented in Tables 3 and 4, respectively.The Gauss point locations for an element are clearly shown in Figure 4.The Q-1115 graphite-epoxy is considered as fabrication material for the shells under present study and the corresponding material properties are reported in Table 5.The shells are subjected to uniformly distributed transverse static loading.The failure loads, failure locations on the shell surface, and the failed lamina, counted from the top of the laminate downwards, correspond to the first ply failure and the ultimate ply failure loads of the laminated composite skewed hypar shells are reported in Table 6.    the deflections remain arrested below 4 mm even when the FPF load is exceeded.The geometry of a skewed hypar shell is such that the shell tends to transfer the loads in the diagonal directions.In the angle ply laminates that are selected here, the fiber directions are aligned along the diagonal directions too.This is why the angle ply shells are in general stiffer than the cross ply ones and deflect less even after getting partly damaged.

Failure Progress.
The failure criterion from which the FPF load is obtained is more important than that of UPF load.This is why the failure progress according to this criterion is studied pictorially in discrete steps and lamina wise.Typical antisymmetric and symmetric angle ply shells are taken up for such study where the failed elements are made dark corresponding to first to fifth ply failure loads and corresponding to 0.25 UPF load, 0.5 UPF load, 0.75 UPF load, Note: "L" represents the least failure load obtained from the four failure criteria.The lamina-wise failure progress mechanism with gradually increasing failure load and knowledge of failure propagation may be utilised for getting a warning in advance about the gross failure of the shell by appropriate instrumentation related to different nondestructive tests (NDT).If an engineer has a prior knowledge regarding the flow of damage then appropriate points on the shell surfaces may be tested only to detect the hidden flaws without going into an elaborate testing program.The maximum area of damage can be allowed under a practical situation may be correlated with the NDT results at predetermined points.
For example, in case of a 0 ∘ /90 ∘ shell, if one has an ultrasonic flaw detection point set up at  = 12.5 cm and  = 87.5 cm, then once he/she detects a flaw at that point he/she can be sure that area of damage is about 0.8%.In other words if he/she wants to restrict the shell within a maximum limit of 0.8% damage he/she should be vigilant about the ultrasonic flaw detection results only at the point as mentioned above without going into elaborate instrumentation.

Conclusion
The present finite element code is capable of successfully assessing the first ply and progressive failure behaviour of laminated composite skewed hypar shells with simply supported boundary condition.The present study indicates that angle ply shells should be preferred compared to cross ply ones from engineering point of view.The location of the first ply failure point is extremely important to be known to a practising engineer because any instrumentation needed for hidden flaw detection should start from that point.In other words, if it is found that the point prone to first ply failure damage is free of any hidden flaw, it can safely be concluded that the shell surface is free of any damage caused due to overloading.

𝑎, 𝑏:
Length and width of the shell : R i s eo fh y p a rs h e l l , V, : In-plane and transverse displacements , : Slopes along the plan directions  11 ,  22 ,  33 : Elastic moduli  12 ,  13 ,  23 : Shear moduli   : T w i s tr a d i u so fc u r v a t u r eo ft h eh y p a r shell 1, 2, and 3: Local coordinates of a lamina , , : Global coordinates of the laminate   ,   : In-plane normal strain components along and -axis, respectively   : In-plane shear strain in - plane  1 ,  2 : In-plane normal strains along 1 and 2 axes of a lamina, respectively
and   = Tensile strength along 1st and 2nd direction, respectively.  and   = Compressive strength along 1st and 2nd direction, respectively. 1 and  2 = Normal stress along 1st and 2nd direction, respectively. 12 = Shear stress in the 1-2 plane. = Shear strength of the material.

Figure 3 :
Figure 3: Schematic algorithm of progressive failure analysis.

Table 3 :
Geometric dimensions of the hypar shell.

Table 4 :
Different types of lamina.
for the cross ply shells, the FPF load corresponds to deflection greater than 4 mm and for the angle ply shells the deflection corresponding to FPF load is less than 4 mm.In case of angle ply shells when the failure progress is deeply examined it is found that over a wide range of load values one of the plies remains almost undamaged and this accounts for the fact that

Table 6 :
First ply failure load and ultimate ply failure load.