Residual Strength and Fatigue Lifetime of ... - Solid Mechanics
Residual Strength and Fatigue Lifetime of ... - Solid Mechanics Residual Strength and Fatigue Lifetime of ... - Solid Mechanics
B 2 1/ 2 ( S44 S55 S45) (4.17) where S44, S55 and S45 are compliance elements given by 1 S 44 G 23 1 S55 (4.18) G 13 S45 is zero for on-axis directions but appears for off-axis directions if G13G23. The total strain energy release rate is given by G G G (4.19) I II Mode II III x deflection at the crack tip Figure 4.13: Definition of x, y and z at the crack tip. The decomposition of the strain energy release rate to two components of GI+II and GIII is considered more useful for practical applications due to a present lack of experimental characterisation aimed at measuring the effect of GIII in terms of crack growth rate. In all recent studies the main focus has been on measuring the crack growth rate under pure mode I, II or mixed-mode loading at the crack tip. Due to difficulties associated with the mode III loading of the crack tip, this component has always been neglected. Thus, in the numerical routine presented, only the GI+II component of the energy release rate is used in the crack growth algorithm. This may introduce inaccuracy in the debond growth simulation if the mode III energy release rate contribution is large. These possible inaccuracies will be discussed later in this chapter. Debonded sandwich panels consisting of 2 mm thick plain-woven E-glass/polyester face sheets over 50 mm thick Divinycell H45 PVC foam are considered the simulation. Face sheet and core material properties are similar to those of the sandwich beam specimen analysed earlier, as listed 78 Mode I y deflection at the crack tip Mode III z deflection at the crack tip
in Table (4.1). The debonded panels are square with a side length of 310 mm. An elliptical face/core debond with a short radius (b) of 45 mm and a large radius (a) of 76.5 is created at the centre of the panel. 8-node isoparametric brick elements (SOLID45) are used in the finite element model. Due to the current lack of suitable experimental fatigue crack growth rate data, the crack growth rate vs. strain energy release rate is simply assumed to be constant for modemixity phase angles larger and smaller than -10 degrees and chosen arbitrarily as da dN da dN 2 0. 000005GI II for >-10 (4.20) 2 0. 000002GI II for -10 (4.21) where GI+II is the difference between maximum and minimum strain energy release rate in each cycle and da/dN is the crack growth rate. The simulation is conducted in load control with a maximum amplitude of 0.35 kN and loading ratio of R=Fmin/Fmax=0.1. To investigate the distribution of mode I, II and III components of strain energy release rate and mode-mixity phase angle along the debond front, radar diagrams from the analysis of the debonded panels exposed to maximum amplitude of the fatigue load are shown in the following figures. Debonded panels with a short radius of 45 mm and a ratio of large radius/short radius (a/b) of 1.7, 1.4 and 1.1 are analysed. In the diagrams 0 and 90 degrees correspond to the points on the debond front on the short and large radiuses of the ellipse. Figures 4.14 (a) and 4.14 (b) illustrate the distribution of mode I+II energy release rate (GI+II) and the related phase angle in the first cycle along the debond front. Maximum GI+II and mode-mixity phase angle occur at the short ellipse radius because of smaller crack length and decrease towards the larger radius. This can be attributed to the development of membrane forces in the face sheet at larger radiuses. As the radius of the ellipse increases the membrane forces become larger, and a subsequently larger part of the strain energy in the specimen should be used to stretch the debonded face sheet rather than create new crack surfaces, decreasing the energy release rate at the crack tip. As the ratio a/b decreases to one (circle) distribution of both GI+II and mode-mixity, the phase angle becomes more even as expected. The mode-mixity phase angle for all a/b ratios is between -5 and -10 degrees along the debond front, which indicates a mode I dominated loading at the crack tip. The mode III strain energy release rate along the debond front is shown in Figure 4.15. In the symmetry plane (0 and 90 degrees) - due to the symmetry effect and the boundary conditions - the out-of-plane deformation (crack plane) at the crack flanks is zero and consequently the mode III strain energy release rate is zero. The maximum GIII on the panels with an a/b ratio of 1.7 is almost 9% of the maximum GI+II , implying the importance of mode III loading at the crack tip in the elliptical debond case with a large a/b ratio. The mode III strain energy release rate is very small for the a/b ratio of 1.1 and is not shown in the diagrams. For debonds with a small a/b ratio the debond is close to a circle and the mode III effects are insignificant. Figure 4.15 reveals that the maximum mode III crack tip loading occurs close to the longer radius of the ellipse (around 79
- Page 50 and 51: previous observations of crack path
- Page 52 and 53: of contact elements (CONTACT173 and
- Page 54 and 55: the debond opening initially increa
- Page 56 and 57: Table 2.4: Instability loads determ
- Page 58 and 59: G (J/m2) 600 400 200 0 H100 IMP=0.1
- Page 60 and 61: This page is intentionally left bla
- Page 62 and 63: Hayman (2007) has described a damag
- Page 64 and 65: Table 3.1: Panel test specimens. La
- Page 66 and 67: sheet and the core thickness, respe
- Page 68 and 69: Load (N) 250 200 150 100 50 0 H130
- Page 70 and 71: Table 3.4: Parameters in the face/c
- Page 72 and 73: was used to monitor 3D surface disp
- Page 74 and 75: arising due to the junction between
- Page 76 and 77: (a) Debonded face sheet Figure 3.16
- Page 78 and 79: Figure 3.19 shows the determined en
- Page 80 and 81: H130 MMB H130 Panel H250 MMB Interf
- Page 82 and 83: 3.6 Conclusion In this chapter the
- Page 84 and 85: This page is intentionally left bla
- Page 86 and 87: composite laminates under thermal c
- Page 88 and 89: S S y( t ) y( t ) ( t ) 2 1 12 2
- Page 90 and 91: listed in Table 4.1. The length and
- Page 92 and 93: The strain energy release rate, G,
- Page 94 and 95: high growth rate of G (see Figure 4
- Page 96 and 97: 4.4 Face/Core Fatigue Crack Growth
- Page 98 and 99: Debond 310 mm Symmetry B. C. a Figu
- Page 102 and 103: 75), which illustrates possible ina
- Page 104 and 105: Figure 4.18 (a) presents the deflec
- Page 106 and 107: Debond radius (mm) 100 90 80 70 60
- Page 108 and 109: using the cycle jump method, more t
- Page 110 and 111: Chapter 5 Face/Core Interface Fatig
- Page 112 and 113: of this chapter, sandwich panels wi
- Page 114 and 115: H250 Specimen H100 Specimen H45 Spe
- Page 116 and 117: Figure 5.5: Test setup. Initially,
- Page 118 and 119: H100 Specimen Fibre bridging Figure
- Page 120 and 121: propagation, the crack continues to
- Page 122 and 123: Figure 5.14: Kinking of the crack i
- Page 124 and 125: Crack length [mm] Crack length [mm]
- Page 126 and 127: mode-mixity phase angle was chosen
- Page 128 and 129: should be taken into account for a
- Page 130 and 131: Figure 5.25: Kinking of the crack i
- Page 132 and 133: log (da/dN) (mm/cycle) 10 log G (J/
- Page 134 and 135: erroneous extrapolations in the tra
- Page 136 and 137: compared to the 65% efficiency obta
- Page 138 and 139: the case of uneven debond growth, t
- Page 140 and 141: Load (kN) 2 1.5 1 0.5 0 Panel 1 Pan
- Page 142 and 143: Figure 5.42: Zero and ninety degree
- Page 144 and 145: G(J/m 2 ) 180 150 210 120 150 100 5
- Page 146 and 147: 110 q=qG=0.4 Test #1 100 Test #2 Te
- Page 148 and 149: panels were determined for differen
B <br />
2 1/<br />
2<br />
( S44<br />
S55<br />
S45)<br />
(4.17)<br />
where S44, S55 <strong>and</strong> S45 are compliance elements given by<br />
1<br />
S 44 <br />
G<br />
23<br />
1<br />
S55 (4.18)<br />
G<br />
13<br />
S45 is zero for on-axis directions but appears for <strong>of</strong>f-axis directions if G13G23. The total strain<br />
energy release rate is given by<br />
G G G<br />
(4.19)<br />
I<br />
II<br />
Mode II<br />
III<br />
x deflection at the crack tip<br />
Figure 4.13: Definition <strong>of</strong> x, y <strong>and</strong> z at the crack tip.<br />
The decomposition <strong>of</strong> the strain energy release rate to two components <strong>of</strong> GI+II <strong>and</strong> GIII is<br />
considered more useful for practical applications due to a present lack <strong>of</strong> experimental<br />
characterisation aimed at measuring the effect <strong>of</strong> GIII in terms <strong>of</strong> crack growth rate. In all recent<br />
studies the main focus has been on measuring the crack growth rate under pure mode I, II or<br />
mixed-mode loading at the crack tip. Due to difficulties associated with the mode III loading <strong>of</strong><br />
the crack tip, this component has always been neglected. Thus, in the numerical routine<br />
presented, only the GI+II component <strong>of</strong> the energy release rate is used in the crack growth<br />
algorithm. This may introduce inaccuracy in the debond growth simulation if the mode III<br />
energy release rate contribution is large. These possible inaccuracies will be discussed later in<br />
this chapter.<br />
Debonded s<strong>and</strong>wich panels consisting <strong>of</strong> 2 mm thick plain-woven E-glass/polyester face sheets<br />
over 50 mm thick Divinycell H45 PVC foam are considered the simulation. Face sheet <strong>and</strong> core<br />
material properties are similar to those <strong>of</strong> the s<strong>and</strong>wich beam specimen analysed earlier, as listed<br />
78<br />
Mode I<br />
y deflection at the crack tip<br />
Mode III<br />
z deflection at the crack tip