ProductsAbaqus/Standard Problem descriptionFor 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 discussionResults for the two-dimensional and three-dimensional analyses are discussed in the following sections. Two-dimensional resultsAbaqus 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 MPa. Results for the J-integral in Pa are presented in Table 1. Three-dimensional resultsThe 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
Input filesThe 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.
Tables
FiguresFigure 1. Crack model.
|