Nonuniform crack-face loading and J-integrals

This example verifies nonuniform crack-face loading in two-dimensional and three-dimensional analyses.

The following topics are discussed:

ProductsAbaqus/Standard

Problem description

For the two-dimensional case an edge crack of length 1 m is modeled in a linear elastic specimen. The results are effectively for an infinitely long plate. The geometry is symmetric about the crack line, so only the top half is modeled. The geometry is meshed using CPE8R elements. The crack faces are loaded in five steps. In the first step a load of constant magnitude 1 MPa is applied. In all subsequent steps the load is zero at the surface of the specimen and has magnitude 1 MPa at the crack tip. The load varies linearly in Step 2, quadratically in Step 3, cubically in Step 4, and quartically in Step 5.

For the three-dimensional case the model from 3DDoubleEdgedNotchC3D20_model.py in Contour integral evaluation: two-dimensional case, is modified to apply a uniform crack-face loading via user subroutine DLOAD.

Results and discussion

Results for the two-dimensional and three-dimensional analyses are discussed in the following sections.

Two-dimensional results

Abaqus results are compared with the results taken from page 8.8 of The Stress Analysis of Cracks Handbook by H. Tada, P. C. Paris, and G. R. Irwin. The crack-face loading is given by p=xγ MPa. Results for the J-integral in Pa are presented in Table 1.

Three-dimensional results

The results should be the same as those shown in Contour integral evaluation: two-dimensional case. The results are within a difference of 0.1%.

Python scripts

Nonuniform crack face loading on a two-dimensional model

Run the 2DEdgeCrackCPE8R_model.py script to create the model. Then run the 2DEdgeCrackCPE8R_job.py script to analyze the model. User subroutine DLOAD (2DEdgeCrackCPE8R.for) is used to apply the nonuniform loading on the crack face.

Uniform crack face loading on a three-dimensional model via user subroutine DLOAD

Run the 3DCrackC3D20_model.py script to create the model. Then run the 3DCrackC3D20_job.py script to analyze the model. User subroutine DLOAD (3DCrackC3D20.for) is used to apply the uniform loading on the crack face.

Input files

The input files listed below are provided for users who prefer to use the Abaqus keyword interface instead of Abaqus/CAE. The meshes created in these input files are different from those created by using the Python scripts; however, the results are of the same accuracy.

pjinnu2d.inp

Checks the nonuniform loads applied to plane strain elements via user subroutine DLOAD.

pjinnu2d.f

User subroutine DLOAD used in pjinnu2d.inp.

pjinnu3d.inp

Uses subroutine DLOAD and, therefore, “nonuniform” load types, to apply a uniform load to the crack faces.

pjinnu3d.f

User subroutine DLOAD used in pjinnu3d.inp.

Tables

Table 1. J-integral results in Pa.
γ 1 2 3 4
Tada et al. 17.98 6.67 3.94 2.78 2.27
Abaqus 18.57 6.81 4.01 2.81 2.16

Figures

Figure 1. Crack model.