VM262

VM262
2D Fracture Problem Under Thermal Loading

Overview

Reference:Wilson, W.K. et al. "The Use of the J-integral in Thermal Stress Crack Problems." International Journal of Fracture. (1979): 377-387.
Analysis Type(s):Static (ANTYPE = 0)
Element Type(s):2D Structural Solid Elements (PLANE182)
Input Listing:

vm262.dat

VM262 requires a supplemental .cdb input file which is too long to include full input listings. This file must be downloaded and placed in your working directory for the test case to run properly. Additionally, the geometry and mesh should be regenerated.

Download link: MAPDL Test Case Files for 2024 R2

vm262.cdb

Test Case

The problem deals with an edge cracked strip with its ends constrained. The strip is subject to a linear temperature gradient through the thickness starting with zero at mid thickness and reaching its final value T0 at the right end edge. Stress intensity factor for the cracked strip is calculated and compared against analytical value.

Figure 452: 2D Fracture Problem Sketch

2D Fracture Problem Sketch

Material Properties

Geometric Properties

Loading

E = 1 x 105 Pa

µ = 0.3

Thermal expansion, α = 1 x 10-4

Crack length = 1 m

L = 4 m

W = 2 m

T0 = 10 °

Analysis Assumptions and Modeling Notes

The problem is solved using 2D PLANE182 element with plain strain element behavior. A half plate is modeled and symmetry boundary conditions are considered. The crack tip nodes and the number of paths surrounding the crack tip nodes are defined using CINT command. The plate is subjected to linear temperature loading and the J integral values are computed for the crack tip node. From the J integral values the stress intensity factor KI is calculated.

Results Comparison

 Target Mechanical APDLRatio
Stress intensity, KI126.604128.6331.016