INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL

Volume 2, No 1, 2011

© Copyright 2010 All rights reserved Integrated Publishing Association

RESEARCHARTICLE ISSN - 0976-4259

Article title goes here

Author name1

1- Author Affiliation, do not use Dr. Mr etc.

Author email id

ABSTRACT

Paste your abstract here and format it so that it can match the template, do not edit the header or footer, also remove all section break, page break from the article by using View- Normal option in MS word. Give table name and number in top of table and figure name and number below the figure, use figure instead of fig.1. Reference citation is an important concept, give at least 10-15 references, mostly from recent journal articles and also cite all the references I the journal article content.Paste your abstract here and format it so that it can match the template, do not edit the header or footer, also remove all section break, page break from the article by using View- Normal option in MS word. Give table name and number in top of table and figure name and number below the figure, use figure instead of fig.1. Reference citation is an important concept, give at least 10-15 references, mostly from recent journal articles and also cite all the references I the journal article content.

Keywords:Keyword, keyword, keyword

Nomenclature

E1 = Young’s modulus of the lamina in the fiber direction

E2 = E3 = Young’s modulus of the lamina in the transverse direction of the fiber

G12 = G13 = Shear modulus in the longitudinal plane of the fiber

G23 = Shear modulus in the transverse plane of the fiber

ν12 = ν13 = Poisson’s ratio in the longitudinal plane of the fiber

ν23 = Poisson’s ratio in the transverse plane of the fiber

I-1 = Interface between 00 and 900

I-2 = Interface between 900 and 900

I-3 = Interface between 900 and 00

s = span-to-depth ratio = A/h

A/B = In-Plane aspect ratio

A = Span of the plate

B = Breadth of the plate

h = Thickness of the plate

2a = Length of the minor axis

2b = Length of the major axis

1. Introduction

The increasing use of fiber reinforced laminates in space vehicles, aircrafts, automobiles, ships and chemical vessels have necessitated the rational analysis of structures for their mechanical response. In addition, the anisotropy and non homogeneity and larger ratio of longitudinal to transverse modulii of these new materials demand improvement in the existing analytical tools. As a result, the analysis of laminated composite structures has attracted many research workers, and has been considerably improved to achieve realistic results. Depending upon the nature of application, these structural elements are acted upon by mechanical and thermal loads of varied nature. Usually, the anisotropy in laminated composite structures causes complicated responses under different loading conditions by creating complex couplings between extensions, bending, and shear deformation modes. To capture the full mechanical behavior, it must be described by three dimensional elasticity theories.

Give reference citation as below and not as [1] or (1) etc.

(Bjeletich, Crossman, and Warren, 1979) conducted experiments on symmetrically stacked “quasi isotropic” laminates by altering the stacking sequence of 0o, +45o, and 90o plies. They observed that laminates with certain stacking sequences had a higher tendency to delaminate. (Wang and Crossman, 1982) used a series of [+25o/90on]s graphite epoxy laminates to investigate the thickness effect of one ply on edge delamination. (K. C. Kairouz and F. L. Mathews, 1993) have manufactured and tested the adhesively bonded joint with cross ply adherends having 0° and 90° surface layer. The stresses in the joint were determined using a continuum method of analysis and layer displacement finite element analysis was also undertaken.

Finite element analyses for the prediction of interlaminar stress in symmetric laminates subjected to torsion loading have been conducted by several investigators. A review of published works is provided by (Sen and Fish, 2002). A quasi-three-dimensional finite element stress analysis was used by (Sen and Fish, 2002 and 1996) to determine the stress in glass-epoxy laminates under torsion and combined-torsion loads. The strains were assumed to be independent of the longitudinal coordinate. Thus, only the laminate cross section was modeled.

(S. Raju, R. Sistla and T. Krishnamurthy, 1996) have performed fracture mechanics analysis on two debond configurations; flange-ski strip and skin-stiffener three-dimensional finite element analysis were performed. The stress characterization around the hole in a composite plate under tensile stress had been studied by (Lei Li and Zhu Feng Yue, 2006) using the analytical method. The distribution of the stress components was presented. With the Hashin Failure Criterion, the failure characterization and different failure modes by the analytical method and finite element method were then studied. The possibility of the failure in different modes was compared. (S. K. Panigrahi and Dr. B. Pradhan, 2007) have performed a three-dimensional finite element analysis and computed the out-of-plane normal and shear stresses in an adhesively bonded single lap joint (SLJ) with laminated FRP composite plates which in comparison to other analytical methods for bonded joint analysis, is capable of handling more general situations related to initiation of damages and its growth.

2. Problem Statement

The objective of the present work is to estimate the failure index ’e’ at the interfaces of a four layered symmetric cross-ply laminate using three-dimensional finite element analysis. The finite element models created in ANSYS software are validated and extended to evaluate ‘e’ by varying the span-to-depth ratio (s) and in-plane aspect ratio (A/B) of the plate.

3. Modeling the Problem

3.1 Geometric modeling

A rectangular plate of length A and breadth B=100mm is considered for the present analysis as shown in the Fig. 1. The value of A is calculated from A/B ratio (0.5, 0.75 and 1.0). The thickness of the plate ‘h’ is selected from span-to-depth ratio (s=A/h) which is varied as 40, 60, 80 and 100. The aspect ratio of cut out is maintained as that of plate aspect ratio and the value of ‘a’ is obtained from a fixed value of ‘b’=5mm.

Figure 1: Geometry of laminate with elliptical cut-out

3.2 Finite element modeling

Due to symmetry in geometry, loading and boundary conditions, one-quarter portion of the plate is modeled as shown in Fig. 2. The finite element mesh is generated using SOLID 95 element of ANSYS software (Lei Li, Zhu Feng Yue, 2006). SOLID 95 is a second order brick element. It can tolerate irregular shapes without much loss of accuracy. SOLID 95 elements have compatible displacement shapes and are well suited to model curved boundaries. This element is defined by 20 nodes having three degrees of freedom per node: translations in the nodal x, y, and z directions. The element is suitable to model isotropic as well as orthotropic materials.

Figure 2:FE model on one-fourth portion of laminate

All images should be clear with at least 300 dpi and in layout it should be in line with text

3.3 Boundary Conditions and Loads

Clamped boundary conditions are applied along the side of the plate at y=A. Symmetric boundary conditions are applied on the surfaces at x=0 and y=0 of FE model. A uniform transverse pressure load of 1 MPa is applied on the top surface of the plate.

3.4 Material Properties

The following material properties are considered for the present analysis.

  1. Young’s Modulii, E1= 147GPa, E2 = 10.3GPa, E3= 10.3GPa
  2. Poisson’s Ratios, ν12= 0.27, ν23= 0.54, ν13 = 0.27
  3. Rigidity Modulii, G12= G13= 7GPa, G23= 3.7GPa
  4. Allowable tensile strengths in the three principal material directions:
  5. Xt = 2280MPa, Yt = 57MPa, Zt = 57MPa
  6. Allowable compressive strengths in the three principal material directions:
  7. Xc = 1725MPa, Yc = 228MPa, Zc = 228MPa and
  8. Shearing strengths of the orthotropic layer in various coupling modes:
  9. Sxy = 76MPa, Syz = 76MPa, Sxz = 76MPa
  10. Z = Inter-laminar normal strength = 0.49 MPa
  11. S = Inter-laminar shear strength = 0.96 MPa

3.5 Validation of FE model

The present finite element model is validated by computing the out of plane stresses at the free surface i.e. at the bottom surface of the plate and away from the boundaries. At this location, all the out-of-plane stresses computed are found to be close to zero (Table-1).

Table 1: Validation of the finite element model

Location / σz (MPa) / τyz (MPa) / τzx (MPa)
(Bottom surface and away from boundaries) / 0.093382 / 0.056243 / 0.092651

All table should be within the borders and not more than 5 columns should be used for tables.

4. Results and Discussions

Out-of-planestresses computed at the interfaces in case of a laminate (s=100, A/B=0.5) with an elliptical cut out are shown in Figs. 3 to 11. It is observed that the stress distribution is almost uniform except at the boundaries. At the first interface of the laminate the maximum σz andτyz are found at the mid span and on the hole boundary of the plate and the maximum value of τzx is at x=0 and y=2.5 as shown in Figs. 3 to 5 respectively. At the second interface of the laminate the maximum σz is at x=0 and y=2.5, the maximumτyz is found at the mid span and on the hole boundary of the plate and the maximum value of τzx is at x=4.379512 and y=1.206221 as shown in Figs. 6 to 8 respectively. At the third interface of the laminate the maximum σz andτyz are found at the mid span and on the hole boundary of the plate and the maximum value of τzx is at x=0 and y=2.5 as shown in Figs. 9 to 11 respectively.

After knowing the out-of-plane stresses from the stress analysis, which are responsible for the damages over the different interfaces, the next step is to use these data to compute the location of damage initiation. The failure is considered as a de-lamination damage which is mainly attributed to the inter-laminar stress effects, so only the inter-laminar shear stresses (τyz and τzx) and through-the-thickness normal stress (σz) are required to predict the damage initiation. Under the three dimensional stress states, the failure can be evaluated by the Tsai-Wu quadratic failure criterion equation given below (S.K.Panigrahi, B. Pradhan, 2007).

(σz/Z) 2 + (τyz/S) 2 + (τzx/S) 2 = e 2 if e < 1, no failure

e ≥ 1, failure

From the above equation the failure index values are obtained at three interfaces. It is observed that at every interface of the laminate, the maximum failure index is greaterdhan the critical failure index and the failure occures at every interface of the laminate. The maximum failure index values at all the three interfaces of the laminate are located at the same point i.e. at x=5 and y=0.

Figure 3: Contour plot ofσzat I-1 (s=100, A/B=0.5)

To reduce this failure effect on the laminate, critical pressure load corresponding to e=1 is calculated by solving the problem at 10 equal load steps. For the estimation of pressure load, the graph (Fig. 12) is drawn between the load and the failure index (e) at A/B=0.5 and s=40 of the first interface of the laminate. It is observed that the failure index varies linearly with respect to load. From the linear relationship between the load and the failure index, the critical pressure load is found to be 0.1718 MPa. At critical conditions also the plate undergo failure. To reduce this failure effect, the pressure load of 0.1 MPa is applied to study the effect of other parameters on interlaminar failure of the laminate.

Figure 12: Variation of ‘e’ at I-1 (A/B=0.5, s= 40)

Figures 13 to 15 show the variation of interlaminar failure index with respect to ‘s’ at three interfaces respectively. At the first interface, the laminate is safe up to (s=60, A/B=0.5), (s=40, A/B=0.75), and in remaining cases it undergoes failure as shown in Fig. 13. At the second interface, the laminate is safe up to (s=60 for A/B=0.75 and 1.0), (s=40, A/B=0.5), and in remaining cases it undergoes failure as shown in Fig. 14. At the third interface, the laminate is safe up to (s=80, A/B=0.75), (s=60, A/B=1.0) and in remaining cases it undergoes failure as shown in Fig. 15. From the overall observation of the graphs, the laminate is safe up to (s=40, A/B=0.5). These graphs reveal that, at the first interface of the laminate, the failure index is increasing with span-to-depth ratio (s) and in-plane aspect ratios (A/B), and is minimum at (s=40, A/B=0.5) and maximum at (s=100, A/B=1) as shown in Fig. 13. At the second interface of the laminate the failure index is increasing with ‘s’ and decreasing with (A/B), and is minimum at (s=40, A/B=1) and maximum at (s=100, A/B=0.5) as shown in Fig. 14. At the third interface of the laminate the failure index is increasing with ‘s’ and (A/B), and is minimum at (s=40, A/B=0.5) and maximum at (s=100, A/B=1) as shown in Figure 15.

Figure 13: Variation of ‘e’ with respect to ’s’ at I-1

Figure 14: Variation of ‘e’ with respect to ’s’ at I-2

Figure 15: Variation of ‘e’ with respect to ’s’ at I-3

5. Conclusions

Interlaminar failure analysis of a four layered cross-ply laminate with elliptical cutout under cylindrical bending subjected to transverse pressure is performed using three dimensional FEM. The following conclusions are identified from the above failure analysis

  1. The laminate is safe up to (s=40, A/B=0.5)
  2. Failure index is maximum for s=100 and A/B=0.5.
  3. It is identified that the failure index ‘e’ increasing with ‘s’ and (A/B) at first and third interfaces of the laminate, and is increasing with ‘s’ and decreasing with (A/B) at second interface of the laminate.

In any case there should be only 6 sections in the article and not more than that, reduce the sections to only 6 in total

6. References

  1. Kairouz K. C. and Mathews F. L., (1993), Strength and Failure Modes of Bonded Single Lap Joints between Cross-ply Adherends, Composite Structures , 24(6),pp 475-484.
  2. Sen J. K. and Fish J. c., Zhang C., Di S. and Zhang N., (2002), A new procedure for static analysis of thermo-electric laminated composite plates under cylindrical bending, Composite Structures, 56 , pp 131-40.
  3. Sen J. K. and Fish J. C., (1996), Failure prediction of composite laminates under Torsion, key engineering materials, Fracture of composites (Editor: E. Armanios Transtec publications) 121-122,pp 285-306.
  4. RajuI. S., Sistla R. and Krishnamurty T., (1996), Fracture Mechanics Analysis for Skin-Stiffener Debonding, Engineering Fracture Mechanics (1996); 54(I),pp 371-385.
  5. Lei Li, Zhu Feng Yue, (2006), Analysis of Interlaminar Stresses and Failure around HoleEdgefor Composite Laminates under In-Plane Loading, Key Engineering Materials 324 – 325,pp 1027-1030.
  6. S. K. Panigrahi and Dr. B. Pradhan, (2007), Three Dimensional Failure Analysis and Damage Propagation Behavior of Adhesively Bonded Single Lap Joints in Laminated FRP Composites, Journal of Reinforced Plastics and Composites 26,pp 183-201.

1