Journal of Structures

Journal of Structures / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 147578 | 8 pages | https://doi.org/10.1155/2014/147578

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

Academic Editor: Hanxing Zhu
Received19 May 2014
Accepted20 Aug 2014
Published16 Sep 2014

Abstract

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.

1. 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.

2. 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 .

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 and 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: Lamina stresses are obtained using the following constitutive relation of the lamina: where , , , and .

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 stress-strain 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 = Normal stress along 1st direction, that is, along the direction of lamina. = 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 stress in the 1-2 plane. = 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 and = Normal stress along 1st and 2nd direction, respectively. and = Strength of the material in failure along 1st and 2nd direction, respectively. = Shear stress in the 1-2 plane. = 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 and = Tensile strength along 1st and 2nd direction, respectively. and = Compressive strength along 1st and 2nd direction, respectively. and = Normal stress along 1st and 2nd direction, respectively. = Shear stress in the 1-2 plane. = Shear strength of the material.

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. and = Normal stress along 1st and 2nd direction, respectively. = 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.

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.


Failure criteriaSide/thicknessFirst ply failure loads (Kam et al. [9]) NewtonFirst ply failure loads (present formulation) Newton

Maximum stress105.26108.26112.14
Hoffman106.45104.40
Tsai-Wu112.77110.50
Tsai-Hill107.06104.40


LaminationReddy [17]Present approach

00/90016.98017.010
00/900/006.6976.701
00/900/900/006.8336.836

, , , , , ,  N/cm2, and .

Hypar shell dimensionsValues mm

Length ()1000
Width ()1000
Thickness ()10
Rise of hypar shell ()200


Lamina nameStacking orders

Antisymmetric Cross Ply (ASCP)00/900
Symmetric Cross Ply (SYCP)00/900/00
Antisymmetric Angle Ply (ASAP)450/−450
Symmetric Angle Ply (SYAP)450/−450/450


Material constantsStrengths

142.5GPa2193.5MPa
9.79GPa2457.0MPa
4.72GPa41.3MPa
1.192GPa206.8MPa
0.2761.28MPa
0.2578.78MPa


Type of lamina Failure criteriaFirst ply failure Ultimate Ply Failure Load (UPFL) MPa Load for 4 mm deflection where such load is less than FPFL MPa
Failure load (FPFL) MPaFailed plyFailed Gauss pointFailed element

ASCPMaximum stress 0.8238118180.05500.1993
Hoffman0.7714L118190.4308
Tsai-Hill0.8177118179.0895L
Tsai-Wu 0.7722218191.8979

SYCPMaximum stress 0.7980118265.7825 0.1870
Hoffman0.7244L218296.7450
Tsai-Hill 0.7922118265.3874L
Tsai-Wu 0.7249218300.7547

ASAPMaximum stress 17.79682159499.8536
Hoffman0.8798111495.2707
Tsai-Hill 0.99871349509.8078
Tsai-Wu 0.8122L111493.5450L

SYAPMaximum stress 25.64081157504.5590
Hoffman1.3737148356.7432
Tsai-Hill 1.51891157418.6368
Tsai-Wu 1.2917L148328.4164L

Note: “L” represents the least failure load obtained from the four failure criteria.

4. Results and Discussion

4.1. First Ply and Ultimate Ply Failure Load

The first ply failure load and ultimate ply failure load for composite skewed hypar shells are mentioned in Table 6. Both the FPF and UPF loads are evaluated using four different well-established failure criteria and the minimum value among those is picked up as the failure load. Symmetric angle ply laminate gives the maximum FPF load value while antisymmetric angle ply shows the highest UPF load value. In general the performances of the angle ply shells are better. It may be inferred from the results that among two laminates the one stiffer in terms of FPF may not be the stiffer one in terms of the UPF load. It is also interesting to note that the UPF load is by far higher than the FPF load. In practical civil engineering apart from strength, serviceability is also considered as an important parameter in design. In the present study the permissible deflection is assumed as span/250 (according to Indian Standard IS: 456-2000) which is 4 mm for the shell geometry considered here. It is very interesting to note that 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 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.

4.2. 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, and UPF load. The typical failure pattern of ASCP lamina is shown in Figures 5(a), 5(b), and 5(c).

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  cm and  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.

5. 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.

Notations

:Length and width of the shell
:Rise of hypar shell
:In-plane and transverse displacements
:Slopes along the plan directions
:Elastic moduli
Shear moduli
:Twist radius of curvature of the hypar 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
, :In-plane normal strains along 1 and 2 axes of a lamina, respectively
:In-plane shear strain in 1-2 plane of a lamina
In-plane lamina stresses
:Poisson’s ratio
, :Natural co-ordinates of isoparametric elements
, , :Curvatures of the shell due to load
:Interpolation functions
and :Normal stress strengths of a lamina along the fiber direction in tension and compression respectively
and :Normal stress strengths of the matrix along the perpendicular to the fiber direction in tension and compression respectively
:Shear stress strengths of a lamina
:Intensity of uniformly distributed load
:Transverse displacement in cm
:Non-dimensional transverse displacement of shell .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The first author gratefully acknowledges the financial assistance of Technical Education Quality Improvement Programme, Phase-II (a World Bank aided project of Government of India) of serial no. 1893 of Jadavpur University, Kolkata 700032, India.

References

  1. S. Sahoo and D. Chakravorty, “Finite element vibration characteristics of composite hypar shallow shells with various edge supports,” Journal of Vibration and Control, vol. 11, no. 10, pp. 1291–1309, 2005. View at: Publisher Site | Google Scholar
  2. S. Sahoo and D. Chakravorty, “Finite element bending behaviour of composite hyperbolic paraboloidal shells with various edge conditions,” Journal of Strain Analysis, vol. 39, no. 5, pp. 1–15, 2004. View at: Google Scholar
  3. A. N. Nayak and J. N. Bandyopadhyay, “On the free vibration of stiffened shallow shells,” Journal of Sound and Vibration, vol. 255, no. 2, pp. 357–382, 2003. View at: Publisher Site | Google Scholar
  4. S. B. Singh and A. Kumar, “Postbuckling response and failure of symmetric laminates under in-plane shear,” Composites Science and Technology, vol. 58, no. 12, pp. 1949–1960, 1998. View at: Publisher Site | Google Scholar
  5. G. Akhras and W. C. Li, “Progressive failure analysis of thick composite plates using the spline finite strip method,” Composite Structures, vol. 79, no. 1, pp. 34–43, 2007. View at: Publisher Site | Google Scholar
  6. R. Ganesan and D. Y. Liu, “Progressive failure and post-buckling response of tapered composite plates under uni-axial compression,” Composite Structures, vol. 82, no. 2, pp. 159–176, 2008. View at: Publisher Site | Google Scholar
  7. J. N. Reddy and A. K. Pandey, “A first-ply failure analysis of composite laminates,” Computers and Structures, vol. 25, no. 3, pp. 371–393, 1987. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  8. T. Y. Kam and T. B. Jan, “First-ply failure analysis of laminated composite plates based on the layerwise linear displacement theory,” Composite Structures, vol. 32, no. 1-4, pp. 583–591, 1995. View at: Publisher Site | Google Scholar
  9. T. Y. Kam, H. F. Sher, T. N. Chao, and R. R. Chang, “Predictions of deflection and first-ply failure load of thin laminated composite plates via the finite element approach,” International Journal of Solids and Structures, vol. 33, no. 3, pp. 375–398, 1996. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  10. B. G. Prusty, S. K. Satsangi, and C. Ray, “First ply failure analysis of stiffened panels—a finite element approach,” Composite Structures, vol. 51, no. 1, pp. 73–81, 2001. View at: Publisher Site | Google Scholar
  11. Y. V. Satish Kumar and A. Srivastava, “First ply failure analysis of laminated stiffened plates,” Composite Structures, vol. 60, no. 3, pp. 307–315, 2003. View at: Publisher Site | Google Scholar
  12. P. Pal and C. Ray, “Progressive failure analysis of laminated composite plates by finite element method,” Journal of Reinforced Plastics and Composites, vol. 21, no. 16, pp. 1505–1513, 2002. View at: Publisher Site | Google Scholar
  13. T. P. Philippidis and A. E. Antoniou, “A progressive damage FEA model for glass/epoxy shell structures,” Journal of Composite Materials, vol. 47, no. 5, pp. 623–637, 2013. View at: Publisher Site | Google Scholar
  14. B. Ellul, D. Camilleri, and J. C. Betts, “A progressive failure analysis applied to fiber-reinforced composite plates subject to out-of-plane bending,” Mechanics of Composite Materials, vol. 49, no. 6, pp. 605–620, 2014. View at: Publisher Site | Google Scholar
  15. S. Tolson and N. Zabaras, “Finite element analysis of progressive failure in laminated composite plates,” Computers and Structures, vol. 38, no. 3, pp. 361–376, 1991. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  16. B. G. Prusty, “Progressive failure analysis of laminated unstiffened and stiffened composite panels,” Journal of Reinforced Plastics and Composites, vol. 24, no. 6, pp. 633–642, 2005. View at: Publisher Site | Google Scholar
  17. J. N. Reddy, “Exact solutions of moderately thick laminated shells,” Journal of Engineering Mechanics, vol. 110, no. 5, pp. 794–809, 1984. View at: Publisher Site | Google Scholar

Copyright © 2014 Arghya Ghosh and Dipankar Chakravorty. 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.

565 Views | 210 Downloads | 2 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder