Flaws and defects in structures and components sometimes lead to disastrous results even though the stress level in a perfect structure may indicate a satisfactory design. The study of crack initiations and growth is a complicated topic which involves physics, chemistry, mechanics, and so on. The engineering field of fracture mechanics was established to develop a basic understanding of such crack propagation problems. For example, engineers usually want to know the conditions under which an existing crack will continue to grow.
The fracture analysis supports the following features:
- 2D and 3D brick models with elastic and elasto-plastic nonlinear materials. Elements can include midside nodes. Tetrahedron elements are supported for Linear Static and MES analyses.
- Analysis types are supported are the following:
- Static Stress with Linear Material Models
- Mechanical Event Simulation (MES)
- Static Stress with Nonlinear Materials
- MES Riks Analysis.
- Thermal loads.
- Dynamic loads and effects (although some of the calculation procedures are based on quasi-static theory).
- Crack growth direction and stress intensity factors based on mode separation.
Furthermore, the theory of fracture mechanics implemented are based on the following:
- The cracks do not propagate during the analysis.
- The surfaces of the crack do not have any loads applied to them.
Procedure
Fracture analysis is a post-processing function, meaning that the stress analysis is performed first, and the fracture analysis is performed on the existing results in the Results environment (post-processing). The basic steps to performing a fracture analysis are as follows:
- Create a mesh of the model that includes the crack or defect. A crack is a separation between elements. The tip of the crack should be a sharp point. See Figure 1.
- Perform a stress analysis.
- Define the parameters of the crack or cracks in the Results environment using the tree view.
- Run the fracture analysis from the Results environment.
- Review the results.
A part with a crack (exaggerated)
|
|
detail of mesh around crack.
|
|
Key: |
The crack modeled in the part. |
|
Nodes on the crack front (or crack tip) |
|
Approximate direction of the crack propagation. |
Figure 1: Crack Terminology
|
Define the Crack
The parameters associated with each crack are entered in the Results environment using the Fracture Analysis branch of the tree view. Right-click the branch and choose New to add a crack (or Edit to edit an existing crack definition). These actions open the Fracture Crack Definition dialog. The input is entered on the following tabs: (Refer also to Figure 1.)
Front tab
The Front tab is used to specify the nodes that belong to a crack (the crack front) and the approximate direction in which the crack will grow.
- Use the Selection command to select the nodes at the tip of the crack. If midside nodes are included in the part, those nodes must also be selected. A given node cannot be assigned to more than one crack.
- Click the Add Selected button to add the nodes to the list.
- Using the propagation direction X, Y, and Z columns, enter the virtual crack extension direction (the direction normal to the crack front). This direction affects the results, especially the calculations of KI and KII. In general, the direction normal to the plane of the crack, the tangential direction of the crack front and this virtual extension direction generate an orthogonal frame. See Figure 2.
|
A slice through the plane of the crack. Axis N is perpendicular to the plane of the crack. Axis T is tangent to the crack at each point. (The crack front is shown in red.) Axis V is the virtual crack extension direction and is perpendicular to both N and T.
Figure 2: Virtual Extension Direction
|
Tip: If the direction of the crack growth is the same for all nodes, enter the X, Y and Z directions for the first row. Then highlight the row, right-click, and choose Copy Direction to Others. If direction is the same for a series of sequential nodes but changes in groups, then add the first node to the list, enter the direction, then add the remaining nodes from the group to the list. The crack direction for newly added nodes is set to the same as the last node in the list. Repeat for the next group of nodes.
- If the crack front has no beginning and end, then activate the Crack Front is Closed option. For example, a crack around the outside of a cylinder is continuous; therefore this option needs to be activated. A crack across the thickness of a bar begins on one face and ends at the other face; this option does not need to be activated.
Highlighting a row in the Front list will highlight the corresponding node on the model. Use the Remove button as needed to remove the selected rows from the list.
Note:
- The order of the nodes in the list needs to be in the same order as the crack front, from one end of the crack to the other. If a node is out of order, highlight the row in the list, right-click, and choose Move Up or Move Down to move the selected row.
- The node list can also be sorted based on the X, Y, or Z coordinates of the nodes (not to be confused with the X, Y, Z crack propagation directions shown in the list). If only one row is highlighted when you right-click, the options Sort By - Coordinate will sort the entire list by the chosen direction. If multiple rows are selected when you right-click, the options Sort By - Coordinate will sort only the selected rows.
Region tab
Fracture analysis only needs to be performed on the elements surrounding the crack tip. It does not need to be performed on the entire part or model. Therefore, use the Region tab to indicate which elements are used in the analysis. There are two ways to indicate which elements to use:
- Rings of Elements When chosen, use this field to enter how many rings of elements around the crack tip should be selected. When you type a number and press Enter (or click in the Region list), the elements that make up the ring will replace the current entries in the list. See Figure 2. Generally, 3 or 4 rings of elements are sufficient depending on the mesh size. (The goal is to choose the elements around the crack tip with a stress higher than the nominal stress value.)
- Specified When chosen, selected elements are added to the list. Use the Selection command to select the elements around the crack. Then click the Add Selected button to add the elements to the list.
Highlighting a row in the Region list will highlight the corresponding element on the model. Use the Remove button as needed to remove the selected rows from the list.
The numbers correspond to the elements that would be selected when 1, 2, or 3 rings of elements are designated. (Two rings would include all elements shown with '1' and '2'; three rings would include all elements shown with '1', '2', and '3', and so on.)
Figure 2: Examples of Rings of Elements
Output tab
The Output tab controls what calculations are performed.
The choices in the Propagation Direction Criterion control which theory is used to calculate the direction of the crack growth. These criteria predict slightly different angles for the initial crack propagation, but they all have the implication that KII = 0 at the crack tip as the crack extends. The options are as follows:
The fracture analysis can be performed on the chosen load cases or time steps instead of the entire analysis. Enter the first and last load case/time step to be analyzed with the Start and End fields. Generally, only the load cases/time steps around the highest stresses need to be analyzed.
Tip: If a crack already has results and is modified, a confirmation dialog is given when the
OK button is clicked to close the Fracture Crack Definition. The confirmation dialog warns that the results will be deleted and gives the option to continue or not. The confirmation dialog also includes a
Do not warn me about this again check box. If the box is activated, results will be deleted automatically when the crack is modified. To activate the appearance of the confirmation dialog again, use the

Options
Results tab and activate the option
Display warning message when fracture analysis results become invalid due to changes to crack definitions.
Note: The input is saved between sessions. If the model is changed, then the previously entered node and element numbers may no longer be applicable to the new model.
Perform Fracture Analyses
To analyze all the defined cracks, right-click the Fracture Analysis branch of the tree view and choose Analyze All. To analyze individual cracks, select the cracks in the tree view, right-click, and choose Analyze. Analyzing selected cracks does not delete the results for other cracks. The fracture analysis on an MES model may be noticeably slower than on a static model since the fracture analysis computes the velocity and accelerations for each step.
The results for each crack are written to a separate file, and then all the crack results are combined together. The Results environment reads the results from the combined file. If a crack is deleted or modified (right-click the entry in the tree view, then Delete or Edit), the results for the other cracks need to be recombined. These actions could produce a noticeable delay.
View Results
After the fracture analysis is completed, the Results Contours
Other Results panel will contain the following items:
(a) Mode I: tension.
|
(b) Mode II: in-plane shear
|
(c) Mode III: out-of-plane shear
|
Figure 3: Deformation Modes
|
All the results are displayed as vector arrows at the nodes defined along the crack front. The magnitude and color of the arrows represent the magnitude of the selected result. The direction of the arrow represents the direction of the crack growth.
Tip: Select a node (
Selection
Select
Nodes) and inquire on the results (
Results Inquire
Inquire
Current Results) to get the crack propagation direction and the magnitude of the result. The components of the results are the direction components of the crack propagation, multiplied by the result. Thus, the direction of the arrows will change by 180 degrees when the displayed result is negative. The result itself is not a vector and does not have components.