2.1. Fracture-Parameter Calculation Process

Fracture parameters are evaluated during the solution phase of the analysis after a substep has converged.

For detailed information about the calculations for specific parameters, see Fracture-Parameter Calculation Types.

2.1.1. Step 1: Initiate a New Fracture-Parameter Calculation

Specify the CINT command's NEW option and provide a number to identify the input information for the calculation:

CINT,NEW,n

where n is the number identifying this fracture-parameter calculation.

Example 2.1: Initiating a New Fracture-Parameter Calculation

CINT,NEW,1   ! Initiate a new fracture-parameter calculation as # 1 

2.1.2. Step 2: Specify the Type of Fracture-Parameter Calculation

Request the calculation for the type of fracture parameter under consideration:

CINT,TYPE,FractureParameter

Example 2.2: Input to Evaluate J-integrals

CINT,TYPE,JINT

Example 2.3: Input to Evaluate Stress-Intensity Factors (SIFS)

CINT,TYPE,SIFS

2.1.3. Step 3: Define Crack Information

A fracture-parameter calculation requires the following input:

  • The crack-tip (front) nodes

  • The local coordinate systems at the crack tip (front)

Figure 2.1: Schematic Representation of the Crack-Tip Node(s) and Crack-Tip Local Coordinate System (2D and 3D)

Schematic Representation of the Crack-Tip Node(s) and Crack-Tip Local Coordinate System (2D and 3D)

The crack-extension direction (e1), crack-plane normal (e2) at the tip, and the crack-tangent (e3) at the tip together constitute an orthonormal local coordinate system at the crack tip.


Important:  The crack-tip nodes and the local coordinate system must be defined correctly. After specifying all parameters, issue the CINT,PLOT command to inspect the local coordinate system.


Two methods are available (via CINT) for specifying the crack-tip node component and one of the vectors of the local coordinate system. (The other vectors are constructed based on the local element and nodal connectivity information.)

  • Define the crack-tip node component and the crack-plane normal (CINT,CTNC).

    This method enables you to specify a crack-tip (front) node component and the normal to the crack plane. Use this method when the crack plane has the same crack normal at all the crack-tip (front) nodes. Generally, this approach applies to both plane 2D crack geometry and 3D flat crack surfaces.

  • Define the crack-extension node component and crack-extension direction (CINT,CENC).

    This method applies to problems where a unique normal may not exist. The method involves defining the crack-extension node component and the crack-extension direction at each crack-tip node location. Use it in any of the following cases:

    • The crack plane is not flat.

    • A set of nodes forms the crack tip, as in the case of a collapsed crack-tip mesh.

    • 2D axisymmetric problems exist.

    • The crack directions are already known.

2.1.3.1. Define the Crack-Tip Node Component and Crack-Plane Normal (CINT,CTNC)

2D crack geometry: Define a crack-tip node component (usually a node located at the crack tip). You can also define a group of nodes around the crack tip, including the node at the crack tip. Mechanical APDL uses this group of nodes as the starting nodes to form the necessary information for the contour integration.

3D flat crack geometry: Define a crack-tip node component that includes all nodes along the crack front. (Only one node can exist at each node location, however.) All nodes in the crack-tip node component must be connectable, and they must form a line (called the crack front) based on the element connectivity associated with it. Mechanical APDL uses the crack front to determine the elements needed for the contour integration. The procedure is similar to 2D crack geometry, and is done through all the nodes along the crack front.

The command syntax is:

CINT,CTNC,Par1

where Par1 is the component name.

After defining the crack-tip node component, use the CINT command's NORM option to define the normal of the crack plane:

CINT,NORM,Par1,Par2

where Par1 is the coordinate system number and Par2 is the axis of coordinate system defining the normal of the crack plane.

Mechanical APDL calculates crack-extension vector based on the information provided. The crack-extension vector is taken along the perpendicular direction to the plane formed by the crack-plane normal and the tangent direction of the crack-tip node, and is normalized to a unit vector.

Example 2.4: Defining the Crack-Tip Node Component and Crack-Plane Normal

! local coordinate system
LOCAL,11,0,,,,
! select nodes located along the crack front and 
!   define it as crack front/tip node component
NSEL,S,LOC,X,Xctip
NSEL,R,LOC,Y,Yctip
CM,CRACK_TIP_NODE_CM
! Define a new stress-intensity factors calculation
CINT,NEW,1
CINT,TYPE,SIFS
CINT,CTNC,CRACK_TIP_NODE_CM
CINT,NORM,11,2


Important:  Before proceeding, inspect the local coordinate systems at the crack tip to ensure that the crack-extension directions are consistent (CINT,PLOT). If they are inconsistent, you must provide additional information to obtain the correct crack-extension directions. See Inspect the Local Coordinate Systems at the Crack Tip.


2.1.3.1.1. Inspect the Local Coordinate Systems at the Crack Tip

The specification of the crack-normal and the crack-tangent (along the crack-front) is sometimes inadequate for evaluating the crack-extension directions.

After defining crack information for fracture-parameter calculation, plot the crack-tip extension directions and normal (CINT,PLOT) before proceeding with the analysis.


Important:  If the crack-extension directions are inconsistent, you must provide additional information to obtain the correct crack-extension directions.


You can use either of these methods to provide additional information to control crack-extension directions:

If neither method provides the crack-extension directions needed, use the crack-extension node-component approach (CINT,CENC) to define the local coordinate system.

2.1.3.1.1.1. Specify the Crack-Assist Extension Direction (CINT,EDIR)

The crack-extension directions are primarily calculated using the tangent to the crack-front and the crack-normal. The crack-assist node or the crack-assist extension direction is used only to flip the crack-extension in the appropriate direction.

Issue the CINT,EDIR command to provide an approximate crack-extension direction at any node on the crack front. You can either use the vector form or provide a coordinate system and an axis.

An approximate direction at the reference node provided (optional) is sufficient for Mechanical APDL to properly evaluate the crack-extension directions. If the crack-front extension directions at the crack tips change by more than 180 degrees, the reference node is required.

The command syntax is, respectively:

CINT,EDIR,CS,CSID,AXIS,,REFERENCENODE

CINT,EDIR,COMP,XVAL,YVAL,ZVAL,REFERENCENODE

2.1.3.1.1.2. Specify the Crack-Assist Extension Direction via a Crack-Assist Node (CINT,CTNC)

Alternatively, you can issue the CINT,CTNC command to use the assist node to provide an approximate crack-extension direction:

CINT,CTNC,Par1,Par2,Par3

where Par2 defines the crack-extension direction calculation-assist node (any node on the open side of the crack), and Par3 is the crack-front’s end-node crack-extension direction-override flag.

By default, Mechanical APDL uses the external surface to determine the crack-extension direction and normal when the crack-tip node meets the free surface; however, you can use Par3 to override the default behavior with the calculated coordinate system.

2.1.3.2. Define the Crack-Extension Node Component and Crack-Extension Direction (CINT,CENC)

Use this method when the crack plane is not flat, or when a set of nodes forms the crack tip, as in the case of a collapsed crack-tip mesh:

  1. Define a node component consisting of one or more nodes forming the crack tip. The node component can have one or more nodes.

    Example: CINT,CENC,ComponentName

  2. Identify the crack-tip node separately if the node component has more than one node. If a crack-tip node is not identified, the program considers the first node of the node component to be the first node.

    Example: CINT,CENC,ComponentName,Node1

  3. Define the crack-extension direction.

    Identify the local coordinate system associated with the crack under consideration.

    Identify the local axis for the given local coordinate system along which the crack should extend.

    Example: CINT,CENC,ComponentName,Node1,11,2

    Alternatively, define the crack-extension direction by directly specifying the global X Y and Z components of the crack-extension vector.

    Example: CINT,CENC,ComponentName,Node1,,,compx,compy,compz

Repeat for all node locations along the crack front.

Although Mechanical APDL calculates the local coordinate system at the crack tip to determine stress-intensity factors, it is generally a better practice to use the CINT command's NORM option to help align the calculated normals of the crack-tip nodes.

Example 2.5: Defining the Crack-Extension Node Component and Crack-Extension Direction via Local Coordinate Systems (CINT,CENC)

! Crack-extension node component and
!   crack-extension direction specification using local coordinate systems
local,11,0,,,,
local,12,0,,,,
! …
local,n,0,,,,

NSEL,S,LOC,X,Xctip1
NSEL,R,LOC,Y,Yctip1
NSEL,R,LOC,Z,Zctip1
CM,CRACK_FRONT_NODE_CM1
NSEL,S,LOC,X,Xctip2
NSEL,R,LOC,Y,Yctip2
NSEL,R,LOC,Z,Zctip2
CM,CRACK_FRONT_NODE_CM2
! …
NSEL,S,LOC,X,Xctipn
NSEL,R,LOC,Y,Yctipn
NSEL,R,LOC,Z,Zctipn
CM,CRACK_FRONT_NODE_CMn

CINT,NEW,1
CINT,TYPE,SIFS
CINT,CENC,CRACK_FRONT_NODE_CM1,,11,2
CINT,CENC,CRACK_FRONT_NODE_CM2,,12,2
! …
CINT,CENC, CRACK_FRONT_NODE_Cn,,n,2

Example 2.6: Defining the Crack-Extension Node Component and Crack-Extension Direction via Vectors (CINT,CENC)

! Crack-extension node component and
!   crack-extension direction specification using vectors
NSEL,S,LOC,X,Xctip1
NSEL,R,LOC,Y,Yctip1
NSEL,R,LOC,Z,Zctip1
CM,CRACK_FRONT_NODE_CM1
NSEL,S,LOC,X,Xctip2
NSEL,R,LOC,Y,Yctip2
NSEL,R,LOC,Z,Zctip2
CM,CRACK_FRONT_NODE_CM2
! …
NSEL,S,LOC,X,Xctipn
NSEL,R,LOC,Y,Yctipn
NSEL,R,LOC,Z,Zctipn
CM,CRACK_FRONT_NODE_CMn

CINT,NEW,1
CINT,TYPE,SIFS
CINT,CENC,CRACK_FRONT_NODE_CM1,,,,Vx1,Vy1,Vz1
CINT,CENC,CRACK_FRONT_NODE_CM2,,,,Vx2,Vy2,Vz2
! …
CINT,CENCOMP,CRACK_FRONT_NODE_Cn,,,,Vxn,Vyn,Vzn

2.1.4. Step 4: Specify the Number of Contours (if Needed)

The command syntax for specifying the number of contours is:

CINT,NCON,n

where n is the number of contours.

For 3D crack geometry, every node along the crack front has the same number of contours.

Example 2.7: Specifying the Number of Contours

CINT,NEW,1
CINT,TYPE,SIFS
CINT,CTNC,CRACK_FRONT_NODE_SET,NODE
CINT,NCON,6

2.1.5. Step 5: Define a Crack-Symmetry Condition (if Needed)

The command for defining a symmetric condition is:

CINT,SYMM,ON

Example 2.8: Defining a Crack-Symmetry Condition

CINT,NEW,1
CINT,TYPE,SIFS
CINT,SYMM,ON    ! crack #1 is a symmetrical crack

2.1.6. Step 6: Specify Crack-Surface Traction to be Derived from Initial Stress (if Needed)

This command specifies that mesh-independent initial-stress data (if any) are to be converted to crack-surface traction loading:

CINT,CSFL,INIS,STRE

To define the mesh-independent initial-stress data, see Mesh-Independent Initial-State Application in the Advanced Analysis Guide. To read the data, issue INISTATE,READ,Fname,Ext,Path,DOBJ.

For a SMART-based crack-growth analysis, the program applies the traction on the growing crack faces after each remeshing automatically.

For more information, see Table 2.2: Material and Load Support for Fracture-Parameter Calculation and Figure 2.2: Superposition Principle: Initial Stress vs. Crack-Face Traction.

2.1.7. Step 7: Specify Output Controls

Mechanical APDL calculates the fracture parameters during the solution phase of the analysis and stores it in the results file for postprocessing. The output uses all defaults from the OUTRES command.

The command OUTRES,ALL includes CINT command results. To control the output for specific fracture-parameter results only, issue OUTRES,CINT.

Example 2.9: Specifying Output Controls

CINT,NEW,1
CINT,TYPE,SIFS
CINT,CTNC,CRACK_TIP_NODE_CM
CINT,SYMM,ON
CINT,NCON,5
OUTRES,CINT,10   ! Output CINT results every 10 substeps

2.1.7.1. Postprocessing Fracture Parameters

To print or display various fracture parameters during postprocessing, issue PRCINT or PLCINT, respectively.