XFEM
XFEM
Konstantinos Agathos Dipl. Civil Eng. Aristotle University of Thessaloniki Prof. Dr. Eleni Chatzi, Chair of Structural Mechanics, IBK, D-BAUG
Introduction
The Extended Finite Element Method (XFEM) is a numerical method, based on the Finite Element Method (FEM), that is especially designed for treating discontinuities. Discontinuities are generally divided in strong and weak discontinuities.
Displacement jump
Cracked Bar
Strain jump
Bimaterial Bar
Institute of Structural Engineering Method of Finite Elements II 4
Fracture Mechanics
The biggest part of this presentation will be dealing with the modeling of strong discontinuities and more specificaly with cracks. All formulations will be derived for the 2D cracked domain case and in the end the corresponding formulations for weak discontinuities will be given. So some basic concepts of fracture mechanics will be briefly mentioned
Fracture Mechanics
Problem Statement Determine the stress, strain and displacement distribution in structures in the presence of flaws such as cracks and small holes.
Small Hole
Crack
Fracture Mechanics
Problem geometry (cracked domain case)
Y
X
Crack tip Crack plane
Crack front
x z
Crack tip
2D Crack
3D Crack
Fracture Mechanics
Crack opening modes (i.e. how can the two crack surfaces deform)
Mode I
Mode II
Mode III
Fracture Mechanics
Fracture Mechanics
10
Fracture Mechanics
11
Fracture Mechanics
Westergaard solution parameters
Y
12
Fracture Mechanics
FEM Solution
In order to model the crack with FEM, the geometry has to be explicitly represented by the mesh, i.e. nodes have to be placed across the crack and on the crack tip. Example:
FEM Solution
Remarks: Mesh refinement is usually necessary near the crack tips in order to represent the asymptotic fields asociated with the crack tips. As the crack propagates remeshing is needed which is computationally expensive especially in complex geometries and 3D domains. In some cases when remeshing, results need to be projected from one mesh to the other which further increases the computational cost.
15
Partition of Unity
16
Partition of Unity
17
Partition of Unity
In the above equation two factors have to be determined: 1. The type of enrichment functions used (next section). 2. The parts of the approximation that are going to be enriched. In the case of cracks, the nature of the discontinuity is local since stress, strain and displacement fields are discontinuous or singular only near the crack tips or along the crack, so enrichment should be local too, i.e. only nodes near the crack are enriched. This matter is going to be addressed in more detail in a following section.
Institute of Structural Engineering Method of Finite Elements II 18
19
20
Jump Enrichment
A simple example is considered first in order to demostrate the concept, the results will the be generalized to more complex cases: The objective is to represent mesh 1 using mesh 2 plus some enrichment terms
1 2 Y 9 10 4 X 5 11 3 1 2 Y 4 X 5 3
Mesh 1
Institute of Structural Engineering
Mesh 2
Method of Finite Elements II 21
Jump Enrichment
22
Jump Enrichment
23
Jump Enrichment
24
Jump Enrichment
Bimaterial interface
n X
X
d
26
27
28
30
32
33
34
2.5
2.5
2
2
1.5
1.5
1
1
0.5
0.5
Non-regular discontinuity
Institute of Structural Engineering
FE approximation
X r
Crack tip
2D crack
Institute of Structural Engineering
3D crack
Method of Finite Elements II 37
38
39
Numerical Integration
Gauss quadrature is not apropriate for the numerical integration of the discontinuous enrichment functions, so usually one of the following approaches is employed:
Post processing
Apart from stresses strains and displacements, one quantity of interest when post processing XFEM results is the stress intensity factors. Their calculation is based on the evaluation of an integral (interaction integral) over an area around the crack tip. The procedure is similar to the one for the FE case. Stress intensity factors are necessary for the calculation of the stress fields around the crack tip as well as for determining the direction of crack propagation.
41
Extentions
The method can be extended in a very straightforward manner to more general and complex problems such as: Crack propagation Branched and intersecting cracks Plastic enrichment Nonlinear finite elements Dynamic problems
42
References
References/recomended reading: Extended Finite Element Method for Fracture Analysis of Structures by Soheil Mohammadi, Blackwell Publishing, 2008 Extended Finite Element Method for Crack Propagation by Sylvie Pommier, John Wiley & Sons, 2011 Elastic Crack Growth in Finite Elements by T. Belytschko and T. Black, International Journal for Numerical Methods in Engineering, 1999 A Finite Element Method for Crack Growth without Remeshing by N. Mos, J. Dolbow and T. Belytschko, International Journal for Numerical Methods in Engineering, 1999
Institute of Structural Engineering Method of Finite Elements II 43