3D Fracture Mechanics with ADINA
The finite element analysis of components in fracture, including the
possible behavior of cracks, can be of utmost importance.
Specifically, power plants and airplanes need to be continuously
checked for cracks and whether cracks may grow. In the first instance,
linear elastic fracture mechanics analyses are peformed in which the Jintegral is calculated.
ADINA is used effectively for such analyses. Here we feature typical
analysis steps to calculate the Jintegral for a 3D crack with ADINA.
For 3D fracture mechanics calculations, it is customary to use
specialized meshes that are highly refined around the crack front
line. The ADINA User Interface, AUI, includes special features for
creating these meshes. These features are based on the ADINAM geometry modelers, either Parasolid based or Open Cascade based.
To demonstrate the features, we describe the modeling of a typical surface crack in 3D analysis. Symmetry is not used, so there is geometry on both sides of the crack. Figure 1 shows a schematic of this typical crack.
Figure 1 Schematic of typical crack considered
Geometry modeling
The uncracked geometry is modeled using the usual ADINAM modeling features.
The crack front is represented by a geometry edge. This edge is conveniently defined using a local coordinate system oriented with respect to the crack.
The geometry immediately adjacent to the crack front is modeled using two bodies, one body on either side of the crack, as shown in Figure 2. These bodies are called crack front bodies. They are topologically halfcylinders, with one edge of each body corresponding to the crack front, and one face of each body corresponding to the cracked area adjacent to the crack front.
Figure 2 Crack front bodies (shown in exploded view)
The geometry surrounding the crack front bodies is modeled using two additional bodies, one body on either side of the crack, see Figure 3. These bodies are called crack sleeve bodies. Notice that the crack sleeve bodies are large enough to include all of the cracked area not already included in the crack front bodies. Thus the entire crack is surrounded by the crack front bodies and the crack sleeve bodies.
Figure 3 Crack sleeve and crack front bodies (shown in exploded view)
In order to include the crack in the model, the crack front bodies and the crack sleeve bodies are subtracted from uncracked geometry. Facelinks are used to create compatible meshes between all of the bodies.
Meshing
The crack front bodies are meshed with a special mapped meshing command. This command creates a mapped mesh with wedge elements adjacent to the crack front line and elements that fan out radially from the wedge elements. This command automatically does not use nodal coincidence checking on the faces corresponding to the cracked area. There is also an option to move the midside nodes of higherorder elements to the quarterpoints (this option is applied only to the wedge elements adjacent to the crack front line). The mapped meshes considered are shown in Figure 4.
Figure 4 Mapped meshes used for crack front bodies (exploded view)
The crack sleeve bodies are meshed using a special freeform meshing command, see Figure 5. This freeform meshing command automatically does not use nodal coincidence checking on the faces corresponding to the cracked area.
The remaining bodies are meshed as usual using freeform meshing.
Figure 5 Meshes used for bodies (exploded view)
Virtual shifts
The virtual shifts required for fracture mechanics calculations are automatically generated according to the following table:
Virtual shift number 
Corner node number on crack front 
Ring number 
1 
1^{st} corner node 
0 
2 
2^{nd} corner node 
0 

... 

N 
N^{th} corner node 
0 
N + 1 
1^{st} corner node 
1 

... 

2 x N 
N^{th} corner node 
1 

... 

nring x N 
N^{th} corner node 
nring  1 
Here N is the number of corner nodes on the crack front and nring is the number of rings in the mapped mesh surrounding the crack front.
Ring 0 represents the motion of only the corner node at the crack front. Ring 1 represents the motion of the same node, and also the motion of the first ring of elements surrounding that node. This process is repeated for higher ring numbers. The virtual shifts considered are shown in Figure 6.
Figure 6 Virtual shifts illustrated
Due to the path independence of the Jintegral, the Jintegral values for virtual shifts j, j + N, j + 2 x N, ... should be identical. Thus pathindependence can be checked using the virtual shifts.
Example
The movie above shows the meshing used and stresses due to pressure
in a pressure vessel in which a surface crack is present.
The results for the Jintegral are shown in Figure 7.
Figure 7 Plots of the Jintegral for rings 0, 1 and 2
For the first and last corner nodes on the crack front, the virtual shifts for rings 1, 2, ... are invalid because these virtual shifts are on the boundary of the model, the model boundary is curved and the crack front is not normal to the boundary. Hence these virtual shifts perturb the geometric boundary.
For the other corner nodes on the crack front, as we should have, the
pathindependence of the Jintegral is observed.
The modeling features of the ADINA User Interface and calculations
performed using ADINA are clearly of much value in practical fracture
mechanics solutions.
Keywords:
Fracture mechanics, crack front, Jintegral, virtual shift, finite
elements, modeling cracks
