SLIDE 1 Energy Minimizing Multi-Crack Growth in Linear Elastic Fracture Using The Extended Finite Element Method
Danas Sutula
- Prof. Stéphane Bordas
- Dr. Pierre Kerfriden
01/04/2016
SLIDE 2
Content
1. Motivation 2. Problem statement 3. Crack growth 4. Discretization by XFEM 5. Implementation 6. Verification 7. Results 8. Summary
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM 2
SLIDE 3 Problem statement
- Consider a cracked linear-elastic isotropic solid subject to an
external load whose quasistatic behavior can be described by the following total Lagrangian form:
- The solution for u(a) and a(t) are obtained by satisfying the
stationarity of L(u,a) during the evolution of t, subject to Δai≥ 0:
3 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 4 Problem statement
- The solution procedure at time tk consists of
1. solving the variational form for u(ak): 2. advancing the fracture fronts, such that Π(u,ak) → Π(u,ak+1) follows the path of steepest descent while satisfying Griffith’s energy balance
4
Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 5 Crack growth
maximum hoop stress
- Post processing of solution to evaluate SIF [Yau, 1980]
- Crack growth direction [Erdogan & Shi, 1963]
- Growth criterion [Irwin, 1957; Hayashi & Nemat-Nasser, 1981]
5 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 6
- Energy release rate w.r.t. crack increment direction, θi:
- The rates of energy release rates:
- Updated directions (using Newton):
Crack growth
energy minimization
6 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 7
- The discrete potential energy is given by:
- Energy release rate w.r.t. crack increment direction θi :
- The rates of the energy release rate:
Crack growth
energy minimization
7
, where: , where:
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 8 Discretization
XFEM
- Approximation function [Belytschko et al., 2001]
8
singular tip enrichment discontinuous enrichment standard part
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 9
9 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
Implementation
how to compute 𝜀K ?
SLIDE 10
10 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
Implementation
how to compute 𝜀K ?
SLIDE 11 Verification
rotational energy release rates
11
Π vs. θ F F
Test case: square plate with an edge crack with a small kink loaded in vertical tension
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 12 Verification
rotational energy release rates
12
G vs. θ Test case: square plate with an edge crack with a small kink loaded in vertical tension F F
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 13 Verification
rotational energy release rates
13
(topological enr.) dG/dθ vs. θ Test case: square plate with an edge crack with a small kink loaded in vertical tension F F
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 14 Verification
rotational energy release rates
14
(geometrical enr.) dG/dθ vs. θ Test case: square plate with an edge crack with a small kink loaded in vertical tension F F
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 15
Verification
energy min. VS. max-hoop
15
F F Test case: square plate with an inclined center crack in vertical tension
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 16
Verification
energy min. VS. max-hoop
16
Gmin(Π)/Ghoop vs. θ
F F θ Test case: square plate with an inclined center crack in vertical tension
Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 17
Results
10 crack problem
20 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 18
Results
10 crack problem
21 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 19
Results
10 crack problem
22 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 20
Results
10 crack problem
23 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 21
Results
10 crack problem
24 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 22
Results
10 crack problem
25 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 23
Results
10 crack problem
26 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 24
Results
10 crack problem
27 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 25
Results
10 crack problem
29 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 26
Results
double cantilever problem
30 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 27
Results
double cantilever problem
34 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 28
Results
2 edge cracks; internal pressure loading
36 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 29
Results
3 cracks; center crack pressure loading
37 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 30
Results
Edge crack in a PMMA beam with 3 holes
38 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 31
39 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
Results
Edge crack in a PMMA beam with 3 holes
SLIDE 32
40 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
Results
Edge crack in a PMMA beam with 3 holes
SLIDE 33
41 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
Results
Edge crack in a PMMA beam with 3 holes
SLIDE 34
42 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
Results
Edge crack in a PMMA beam with 3 holes
SLIDE 35
Results
2 edge cracks and 2 holes (Khoeil et al. 2008)
43 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 36
Results
2 edge cracks and 2 holes (Khoeil et al. 2008)
24 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 37
Results
2 edge cracks and 2 holes (Khoeil et al. 2008)
45 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM
SLIDE 38
- A robust approach to determining multiple crack growth
based on the principle of minimum energy within XFEM;
- Limitations undermining the max. hoop-stress criterion are
- vercome, e.g. assumptions about geometry and loading;
- The energy minimization approach is characterized by mode-I
field dominance at the crack tip (post-increment);
- Both criteria lead to fracture paths solutions that are in close
agreement (strong correlation with local symmetry, i.e. KII=0);
- Better accuracy and faster convergence of fracture path
solutions can be obtained by taking a bi-section of the interval that is bounded by the respective criteria.
Summary
46 Danas Sutula Energy minimizing multi-crack growth in linear elastic fracture using XFEM