Find the ply-by-ply failure loads for a [0/90]s graphite/epoxy laminate. Assume the thickness of each ply is 5 mm and use properties of unidirectional graphite/epoxy lamina from Table 2.1. The only load applied is a tensile normal load in the x-direction — that is, the direction parallel to the

Question

Find the ply-by-ply failure loads for a [latex][0/\overline{90}]_s[/latex] graphite/epoxy laminate. Assume the thickness of each ply is 5 mm and use properties of unidirectional graphite/epoxy lamina from Table 2.1. The only load applied is a tensile normal load in the x-direction — that is, the direction parallel to the fibers in the 0° ply.

TABLE 2.1
Typical Mechanical Properties of a Unidirectional Lamina (SI System of Units)

Property	Symbol	Units	Glass/ epoxy	Boron/ epoxy	Graphite/ epoxy
Fiber volume fraction	[latex]V_f[/latex]		0.45	0.50	0.70
Longitudinal elastic modulus	[latex]E_1[/latex]	GPa	38.6	204	181
Transverse elastic modulus	[latex]E_2[/latex]	GPa	8.27	18.50	10.30
Major Poisson’s ratio	[latex] u_{12}[/latex]		0.26	0.23	0.28
Shear modulus	[latex]G_{12}[/latex]	GPa	4.14	5.59	7.17
Ultimate longitudinal tensile strength	[latex](\sigma_1^T)_{ult}[/latex]	MPa	1062	1260	1500
Ultimate longitudinal compressive strength	[latex](\sigma_1^C)_{ult}[/latex]	MPa	610	2500	1500
Ultimate transverse tensile strength	[latex](\sigma_2^T)_{ult}[/latex]	MPa	31	61	40
Ultimate transverse compressive strength	[latex](\sigma_2^C)_{ult}[/latex]	MPa	118	202	246
Ultimate in-plane shear strength	[latex]( au_{12})_{ult}[/latex]	MPa	72	67	68
Longitudinal coefficient of thermal expansion	[latex]\alpha_1[/latex]	μm/m/°C	8.6	6.1	0.02
Transverse coefficient of thermal expansion	[latex]\alpha_2[/latex]	μm/m/°C	22.1	30.3	22.5
Longitudinal coefficient of moisture expansion	[latex]\beta_1[/latex]	m/m/kg/kg	0.00	0.00	0.00
Transverse coefficient of moisture expansion	[latex]\beta_2[/latex]	m/m/kg/kg	0.60	0.60	0.60

Accepted Answer

Because the laminate is symmetric and the load applied is a normal load, only the extensional stiffness matrix is required. From Example 4.4, the extensional compliance matrix is

[latex] \left[A^{*} ight]=\left[\begin{array}{ccc} 5.353 imes 10^{-10} & -2.297 imes 10^{-11} & 0 \ -2.297 imes 10^{-11} & 9.886 imes 10^{-10} & 0 \ 0 & 0 & 9.298 imes 10^{-9} \end{array} ight] \frac{1}{P a-m}, [/latex]

which, from Equation (5.1a), gives the midplane strains for symmetric laminates subjected to [latex]N_x = 1 N/m[/latex] as

[latex] \left[\begin{array}{c} \varepsilon_{x}^{0} \ \varepsilon_{y}^{0} \ \gamma_{x y}^{0} \end{array} ight]=\left[\begin{array}{c} 5.353 imes 10^{-10} \ -2.297 imes 10^{-11} \ 0 \end{array} ight] .[/latex]

The midplane curvatures are zero because the laminate is symmetric and no bending and no twisting loads are applied.
The global strains in the top 0° ply at the top surface can be found as follows using Equation (4.16),

[latex] \left[\begin{array}{l} \varepsilon_{x} \ \varepsilon_{y} \ \gamma_{x y} \end{array} ight]=\left[\begin{array}{c} 5.353 imes 10^{-10} \ -2.297 imes 10^{-11} \ 0 \end{array} ight]+(0.0075)\left[\begin{array}{l} 0 \ 0 \ 0 \end{array} ight] \ \space \ =\left[\begin{array}{c} 5.353 imes 10^{-10} \ -2.297 imes 10^{-11} \ 0 \end{array} ight] [/latex].

Using Equation (2.103), one can find the global stresses at the top surface of the top 0° ply as

[latex] \left[\begin{array}{c} \sigma_{x} \ \sigma_{y} \ au_{x y} \end{array} ight]_{0^{\circ} ext {,top }}=\left[\begin{array}{ccc} 181.8 & 2.897 & 0 \ 2.897 & 10.35 & 0 \ 0 & 0 & 7.17 \end{array} ight]\left(10^{9} ight)\left[\begin{array}{c} 5.353 imes 10^{-10} \ -2.297 imes 10^{-11} \ 0 \end{array} ight] \ \space \ =\left[\begin{array}{c} 9.726 imes 10^{1} \ 1.313 \ 0 \end{array} ight] P a. [/latex]

Using the transformation Equation (2.94), the local stresses at the top surface of the top 0° ply are

[latex] \left[\begin{array}{l} \sigma_{1} \ \sigma_{2} \ au_{12} \end{array} ight]_{0^{\circ}, ext { top }}=\left[\begin{array}{lll} 1 & 0 & 0 \ 0 & 1 & 0 \ 0 & 0 & 1 \end{array} ight]\left[\begin{array}{c} 9.726 imes 10^{1} \ 1.313 imes 10^{0} \ 0 \end{array} ight] \ \space \ =\left[\begin{array}{c} 9.726 imes 10^{1} \ 1.313 \ 0 \end{array} ight] P a. [/latex]

All the local stresses and strains in the laminate are summarized in Table 5.1 and Table 5.2.

TABLE 5.1
Local Stresses (Pa) in Example 5.3

Ply no.	Position	[latex]σ_1[/latex]	[latex]σ_2[/latex]	[latex]τ_{12}[/latex]
1 (0°)	Top	9.726 × [latex]10^1[/latex]	1.313 × [latex]10^0[/latex]	0.0
	Middle	9.726 × [latex]10^1[/latex]	1.313 × [latex]10^0[/latex]	0.0
	Bottom	9.726 × [latex]10^1[/latex]	1.313 × [latex]10^0[/latex]	0.0
2 (90°)	Top	–2.626 × [latex]10^0[/latex]	5.472 × [latex]10^0[/latex]	0.0
	Middle	–2.626 × [latex]10^0[/latex]	5.472 × [latex]10^0[/latex]	0.0
	Bottom	–2.626 × [latex]10^0[/latex]	5.472 × [latex]10^0[/latex]	0.0
3 (0°)	Top	9.726 × [latex]10^1[/latex]	1.313 × [latex]10^0[/latex]	0.0
	Middle	9.726 × [latex]10^1[/latex]	1.313 × [latex]10^0[/latex]	0.0
	Bottom	9.726 × [latex]10^1[/latex]	1.313 × [latex]10^0[/latex]	0.0

TABLE 5.2
Local Strains in Example 5.3

Ply no.	Position	[latex]ε_1[/latex]	[latex]ε_2[/latex]	[latex]τ_{12}[/latex]
1 (0°)	Top	5.353 × [latex]10^{-10}[/latex]	–2.297 × [latex]10^{-11}[/latex]	0.0
	Middle	5.353 × [latex]10^{-10}[/latex]	–2.297 × [latex]10^{-11}[/latex]	0.0
	Bottom	5.353 × [latex]10^{-10}[/latex]	–2.297 × [latex]10^{-11}[/latex]	0.0
2 (90°)	Top	–2.297 × [latex]10^{-11}[/latex]	5.353 × [latex]10^{-10}[/latex]	0.0
	Middle	–2.297 × [latex]10^{-11}[/latex]	5.353 × [latex]10^{-10}[/latex]	0.0
	Bottom	–2.297 × [latex]10^{-11}[/latex]	5.353 × [latex]10^{-10}[/latex]	0.0
3 (0°)	Top	5.353 × [latex]10^{-10}[/latex]	–2.297 × [latex]10^{-11}[/latex]	0.0
	Middle	5.353 × [latex]10^{-10}[/latex]	–2.297 × [latex]10^{-11}[/latex]	0.0
	Bottom	5.353 × [latex]10^{-10}[/latex]	–2.297 × [latex]10^{-11}[/latex]	0.0

The Tsai–Wu failure theory applied to the top surface of the top 0° ply is applied as follows. The local stresses are

[latex]σ_1 = 9.726 × 10^1 Pa \ \space \ σ_2 = 1.313 Pa \ \space \ τ_{12} = 0[/latex]

Using the parameters [latex]H_1, H_2, H_6, H_{11}, H_{22}, H_{66}, and H_{12}[/latex] from Example 2.19, the Tsai–Wu failure theory Equation (2.152) gives the strength ratio as

[latex](0) (9.726 × 10^1) SR + (2.093 × 10^{–8}) (1.313) SR+ (0 × 0) + (4.4444 × 10^{–19}) (9.726 × 10^1)^2(SR)^2 + (1.0162 × 10^{–16}) (1.313)^2(SR)^2 + (2.1626 × 10^{–16}) (0)^2 + 2(–3.360 × 10^{–18}) (9.726 × 10^1) (1.313)(SR)^2=1 \ SR = 1.339 × 10^7.[/latex]

The maximum strain failure theory can also be applied to the top surface of the top 0° ply as follows. The local strains are

[latex] \left[\begin{array}{l} \varepsilon_{1} \ \varepsilon_{2} \ \gamma_{12} \end{array} ight]=\left[\begin{array}{c} 5.353 imes 10^{-10} \ -2.297 imes 10^{-11} \ 0.000 \end{array} ight] .[/latex]

Then, according to maximum strain failure theory (Equation 2.143), the strength ratio is given by

[latex]SR = \min \{[(1500 × 10^6)/(181 × 10^9)]/(5.353 × 10^{–10}),[(246 × 10^6)/(10.3 × 10^9)]/(2.297 × 10^{–11})\} = 1.548 × 10^7.[/latex]

The strength ratios for all the plies in the laminate are summarized in Table 5.3 using the maximum strain and Tsai–Wu failure theories. The symbols in the parentheses in the maximum strain failure theory column denote the mode of failure and are explained at the bottom of Table 2.3.

TABLE 5.3
Strength Ratios in Example 5.3

Ply no.	Position	Maximum strain	Tsai–Wu
1 (0°)	Top	1.548 × [latex]10^{7}[/latex] (1T)	1.339 × [latex]10^{7}[/latex]
	Middle	1.548 × [latex]10^{7}[/latex] (1T)	1.339 × [latex]10^{7}[/latex]
	Bottom	1.548 × [latex]10^{7}[/latex] (1T)	1.339 × [latex]10^{7}[/latex]
2 (90°)	Top	7.254 × [latex]10^{6}[/latex] (2T)	7.277 × [latex]10^{6}[/latex]
	Middle	7.254 × [latex]10^{6}[/latex] (2T)	7.277 × [latex]10^{6}[/latex]
	Bottom	7.254 × [latex]10^{6}[/latex] (2T)	7.277 × [latex]10^{6}[/latex]
3 (0°)	Top	1.548 × [latex]10^{7}[/latex] (1T)	1.339 × [latex]10^{7}[/latex]
	Middle	1.548 × [latex]10^{7}[/latex] (1T)	1.339 × [latex]10^{7}[/latex]
	Bottom	1.548 × [latex]10^{7}[/latex] (1T)	1.339 × [latex]10^{7}[/latex]

From Table 5.3 and using the Tsai–Wu theory, the minimum strength ratio is found for the 90° ply. This strength ratio gives the maximum value of the allowable normal load as

[latex] N_{x}=7.277 imes 10^{6} \frac{ N }{ m } [/latex]

and the maximum value of the allowable normal stress as

[latex] \frac{N_{x}}{h}=\frac{7.277 imes 10^{6}}{0.015}, \ \space \ =0.4851 imes10^9 Pa[/latex]

where h = thickness of the laminate.
The normal strain in the x-direction at this load is

[latex] \left(\varepsilon_{x}^{0} ight)_{ ext {first ply failure }}\left(5.353 imes 10^{-10} ight)\left(7.277 imes 10^{6} ight) \ \space \ =3.895 imes 10^{-3} [/latex]

Now, degrading the 90° ply completely involves assuming zero stiffnesses and strengths of the 90° lamina. Complete degradation of a ply does not allow further failure of that ply. For the undamaged plies, the [latex][0/\overline{90}]_s[/latex] laminate has two reduced stiffness matrices as

[latex] [Q]=\left[\begin{array}{ccc} 181.8 & 2.897 & 0 \ 2.897 & 10.35 & 0 \ 0 & 0 & 7.17 \end{array} ight] G P a [/latex]

and, for the damaged ply,

[latex] [Q]=\left[\begin{array}{lll} 0 & 0 & 0 \ 0 & 0 & 0 \ 0 & 0 & 0 \end{array} ight] G P a. [/latex]

Using Equation (4.28a), the extensional stiffness matrix

[latex] A_{i j}=\sum_{k=1}^{3}\left[\bar{Q}_{i j} ight]_{k}\left(h_{k}-h_{k-1} ight) \ \space \ [A]=\left[\begin{array}{ccc} 181.8 & 2.897 & 0 \ 2.897 & 10.35 & 0 \ 0 & 0 & 7.17 \end{array} ight]\left(10^{9} ight)(0.005) \ \space \ +\left[\begin{array}{lll} 0 & 0 & 0 \ 0 & 0 & 0 \ 0 & 0 & 0 \end{array} ight]\left(10^{9} ight)(0.005) \ \space \ +\left[\begin{array}{ccc} 181.8 & 2.897 & 0 \ 2.897 & 10.35 & 0 \ 0 & 0 & 7.17 \end{array} ight]\left(10^{9} ight)(0.005) \ \space \ [A]=\left[\begin{array}{ccc} 181.8 & 2.897 & 0 \ 2.897 & 10.35 & 0 \ 0 & 0 & 7.17 \end{array} ight]\left(10^{7} ight) Pa - m. [/latex]

Inverting the new extensional stiffness matrix [A], the new extensional compliance matrix is

[latex] \left[A^{*} ight]=\left[\begin{array}{ccc} 5.525 imes 10^{-10} & -1.547 imes 10^{-10} & 0 \ -1.547 imes 10^{-10} & 9.709 imes 10^{-9} & 0 \ 0 & 0 & 1.395 imes 10^{-8} \end{array} ight] \frac{1}{P a-m} ,[/latex]

which gives midplane strains subjected to [latex]N_x = 1 N/m[/latex] by Equation (5.1a) as

[latex] \left[\begin{array}{c} \varepsilon_{x}^{0} \ \varepsilon_{y}^{0} \ \gamma_{x y} \end{array} ight]=\left[\begin{array}{ccc} 5.525 imes 10^{-10} & -1.547 imes 10^{-10} & 0 \ -1.547 imes 10^{-10} & 9.709 imes 10^{-9} & 0 \ 0 & 0 & 1.395 imes 10^{-8} \end{array} ight]\left[\begin{array}{l} 1 \ 0 \ 0 \end{array} ight] [/latex]

TABLE 5.4
Local Stresses after First Ply Failure in Example 5.3

Ply no.	Position	[latex]σ_1[/latex]	[latex]σ_2[/latex]	[latex]τ_{12}[/latex]
1 (0°)	Top	1.0000 × [latex]10^{2}[/latex]	0.0	0.0
	Middle	1.0000 × [latex]10^{2}[/latex]	0.0	0.0
	Bottom	1.0000 × [latex]10^{2}[/latex]	0.0	0.0
2 (90°)	Top	—	—	—
	Middle	—	—	—
	Bottom	—	—	—
3 (0°)	Top	1.0000 × [latex]10^{2}[/latex]	0.0	0.0
	Middle	1.0000 × [latex]10^{2}[/latex]	0.0	0.0
	Bottom	1.0000 × [latex]10^{2}[/latex]	0.0	0.0

TABLE 5.5
Local Strains after First Ply Failure in Example 5.3

Ply no.	Position	[latex]ε_1[/latex]	[latex]ε_2[/latex]	[latex]γ_{12}[/latex]
1 (0°)	Top	5.25 × [latex]10^{-10}[/latex]	–1.547 × [latex]10^{-10}[/latex]	0.0
	Middle	5.525 × [latex]10^{-10}[/latex]	–1.547 × [latex]10^{-10}[/latex]	0.0
	Bottom	5.525 × [latex]10^{-10}[/latex]	–1.547 × [latex]10^{-10}[/latex]	0.0
2 (90°)	Top	—	—	—
	Middle	—	—	—
	Bottom	—	—	—
3 (0°)	Top	5.525 × [latex]10^{-10}[/latex]	–1.547 × [latex]10^{-10}[/latex]	0.0
	Middle	5.525 × [latex]10^{-10}[/latex]	–1.547 × [latex]10^{-10}[/latex]	0.0
	Bottom	5.525 × [latex]10^{-10}[/latex]	–1.547 × [latex]10^{-10}[/latex]	0.0

[latex]\left[\begin{array}{c} \varepsilon_{x}^{0} \ \varepsilon_{y}^{0} \ \gamma_{x y}^{0} \end{array} ight]=\left[\begin{array}{c} 5.525 imes 10^{-10} \ -1.547 imes 10^{-10} \ 0 \end{array} ight] .[/latex]

These strains are close to those obtained before the ply failure only because the 90° ply takes a small percentage of the load out of the normal load in the x-direction.
The local stresses in each layer are found using earlier techniques given in this example and are shown in Table 5.4. The strength ratios in each layer are also found using methods given in this example and are shown in Table 5.5.
From Table 5.6 and using Tsai–Wu failure theory, the minimum strength ratio is found in both the 0° plies. This strength ratio gives the maximum value of the normal load as

[latex] N_{x}=1.5 imes 10^{7} \frac{N}{m} [/latex]

TABLE 5.6
Strength Ratios after First Ply Failure in Example 5.3

Ply no.	Position	Maximum strain	Tsai–Wu
1 (0°)	Top	1.5000 × [latex]10^{7}[/latex] (1T)	1.5000 × [latex]10^{7}[/latex]
	Middle	1.5000 × [latex]10^{7}[/latex] (1T)	1.5000 × [latex]10^{7}[/latex]
	Bottom	1.5000 × [latex]10^{7}[/latex] (1T)	1.5000 × [latex]10^{7}[/latex]
2 (90°)	Top	—	—
	Middle	—	—
	Bottom	—	—
3 (0°)	Top	1.5000 × [latex]10^{7}[/latex] (1T)	1.5000 × [latex]10^{7}[/latex]
	Middle	1.5000 × [latex]10^{7}[/latex] (1T)	1.5000 × [latex]10^{7}[/latex]
	Bottom	1.5000 × [latex]10^{7}[/latex] (1T)	1.5000 × [latex]10^{7}[/latex]

and the maximum value of the allowable normal stress as

[latex] \frac{N_{x}}{h}=\frac{1.5 imes 10^{7}}{0.015} \ \space \ =1.0 imes 10^9 Pa[/latex]

where h is the thickness of the laminate.
The normal strain in the x-direction at this load is

[latex] \left(\varepsilon_{x}^{0} ight)_{ ext {last ply failure }}=\left(5.525 imes 10^{-10} ight)\left(1.5 imes 10^{7} ight) \ \space \ =8.288 imes 10^{-3}.[/latex]

The preceding load is also the last ply failure (LPF) because none of the layers is left undamaged. Plotting the stress vs. strain curve for the laminate until last ply failure shows that the curve will consist of two linear curves, each ending at each ply failure. The slope of the two lines will be the Young’s modulus in x direction for the undamaged laminate and for the FPF laminate — that is, using Equation (4.35),

[latex]E_{x}=\frac{1}{(0.015)\left(5.353 imes 10^{-10} ight)} \ \space \ =124.5 GPa[/latex]

until first ply failure, and

[latex] E_{x}=\frac{\left(N_{x} / h ight)_{ ext {last ply failure }}-\left(N_{x} / h ight)_{ ext {first play failure }}}{\left(\varepsilon_{x}^{o} ight)_{ ext {last play failure }}-\left(\varepsilon_{x}^{o} ight)_{ ext {first play failure }}} \ \space \ =\frac{0.1 imes 10^{10}-0.4851 imes 10^{9}}{8.288 imes 10^{-3}-3.895 imes 10^{-3}} \ \space \ = 117.2 GPa[/latex]

after first ply failure and until last ply failure (Figure 5.1).

____________________________________

(2.94): [latex] \left[\begin{array}{l} \sigma_{x} \ \sigma_{y} \ au_{x y} \end{array} ight]=[T]^{-1}\left[\begin{array}{l} \sigma_{1} \ \sigma_{2} \ au_{12} \end{array} ight] [/latex]

(2.103): [latex] \left[\begin{array}{l} \sigma_{x} \ \sigma_{y} \ au_{x y} \end{array} ight]=\left[\begin{array}{lll} \bar{Q}_{11} & \bar{Q}_{12} & \bar{Q}_{16} \ \bar{Q}_{12} & \bar{Q}_{22} & \bar{Q}_{26} \ \bar{Q}_{16} & \bar{Q}_{26} & \bar{Q}_{66} \end{array} ight]\left[\begin{array}{c} \varepsilon_{x} \ \varepsilon_{y} \ \gamma_{x y} \end{array} ight], [/latex]

(2.143):

[latex] -\left(\varepsilon_{1}^{C} ight)_{u l t}<\varepsilon_{1}<\left(\varepsilon_{1}^{T} ight)_{u l t}, ext { or } \ \space \ -\left(\varepsilon_{2}^{C} ight)_{u l t}<\varepsilon_{2}<\left(\varepsilon_{2}^{T} ight)_{u l t}, ext { or } \ \space \ -\left(\gamma_{12} ight)_{u l t}<\gamma_{12}<\left(\gamma_{12} ight)_{u l t} [/latex]

(2.152): [latex] H_{1} \sigma_{1}+H_{2} \sigma_{2}+H_{6} au_{12}+H_{11} \sigma_{1}^{2}+H_{22} \sigma_{2}^{2}+H_{66} au_{12}^{2}+2 H_{12} \sigma_{1} \sigma_{2}<1 [/latex]

(4.16): [latex] \left\{\begin{array}{c} \varepsilon_{x} \ \varepsilon_{y} \ \gamma_{x y} \end{array} ight\}=\left\{ \begin{array}{c} \varepsilon_{x}^{0} \ \varepsilon_{y}^{0} \ \gamma_{x y}^{0} \end{array} ight\}+z\left\{\begin{array}{c} \kappa_{x} \ \kappa_{y} \ \kappa_{x y} \end{array} ight\} [/latex]

(4.28a): [latex] A_{i j}=\sum_{k=1}^{n}\left[\left(\bar{Q}_{i j} ight) ight]_{k}\left(h_{k}-h_{k-1} ight), \quad i=1,2,6 ; \quad j=1,2,6, [/latex]

(4.35): [latex] E_{x} \equiv \frac{\sigma_{x}}{\varepsilon_{x}^{0}}=\frac{N_{x} / h}{A_{11}^{*} N_{x}}=\frac{1}{h A_{11}^{*}} [/latex]

(5.1a): [latex] \left[\begin{array}{c} N_{x} \ N_{y} \ N_{x y} \end{array} ight]=\left[\begin{array}{lll} A_{11} & A_{12} & A_{16} \ A_{12} & A_{22} & A_{26} \ A_{16} & A_{26} & A_{66} \end{array} ight]\left[\begin{array}{c} \varepsilon_{x}^{0} \ \varepsilon_{y}^{0} \ \gamma_{x y}^{0} \end{array} ight] [/latex]

TABLE 2.3
Effect of Sign of Shear Stress as a Function of Angle of Lamina

Angle, Degrees	Positive [latex]τ_{xy}[/latex] MPa	Negative [latex]τ_{xy}[/latex] MPa	Shear strength MPa
0	68.00 (S)	68.00 (S)	68.00
15	78.52 (S)	78.52 (S)	78.52
30	136.0 (S)	46.19 (2T)	46.19
45	246.0 (2C)	40.00 (2T)	40.00
60	136.0 (S)	46.19 (2T)	46.19
75	78.52 (S)	78.52 (S)	78.52
90	68.00 (S)	68.00 (S)	68.00

Question 5.3: Find the ply-by-ply failure loads for a [0/90]s graphite/epo...

Related Answered Questions

Verified Answer:

Repeat Example 5.3 for the first ply failure and use Tsai–Wu failure theory now with an additional thermal load: a temperature change of –75°C. ...

Verified Answer:

Verified Answer:

Verified Answer:

Verified Answer:

Show that the extensional stiffness matrix for a general N-ply quasi-isotropic laminate is given by [A] = [U1 U4 0 U4 U1 0 0 0 U_1-U_4/2]h. where U1 and U4 are the stiffness invariants given by Equation (2.132) and h is the thickness of the laminate. Also, find the in-plane engineering stiffness ...

Verified Answer:

Verified Answer: