Chapter 7: Interface Delamination Analysis of Double Cantilever Beam

This tutorial demonstrates the use of Interface Delamination feature available in Mechanical by examining the displacement of two 2D parts on a double cantilever beam. This same problem is demonstrated in VM248 in the Ansys Mechanical APDL Verification Manual. The following example is provided to demonstrate the steps to set up and analyze the same model using Mechanical.

Analysis TypeInterface Delamination
Features DemonstratedEngineering data/materials, static structural analysis, match control, interface delamination
Licenses RequiredAnsys Mechanical Enterprise/Enterprise Solver/Enterprise PrepPost
Help Resources VM248: Delamination Analysis of Double Cantilever Beam
Tutorial Files 2D_Fracture_Geom.agdb

7.1. Problem Description

As illustrated below, a two dimensional beam has a length of 100 mm and an initial crack of length of 30 mm at the free end that is subjected to a maximum vertical displacement (Umax) at the top and bottom of the free end nodes. Two vertical displacements, one positive and one negative, are applied to determine the vertical reaction at the end point. The point of fracture is at the vertex of the crack and the interface edges.

This image illustrates the dimension of the model.

This tutorial also examines how to prepare the necessary materials and mesh controls that work in cooperation with the Interface Delamination feature.

7.2. Create Static Structural Analysis

  1. Open Ansys Workbench.

  2. On the Workbench Project page, drag a Static Structural system from the Toolbox to the Project Schematic. The Project Schematic should appear as follows. The properties window does not display unless you have made the required selection; right-click a cell and select Properties.


    Note:  The Interface Delamination feature is only available for Static Structural and Transient Structural analyses.


7.3. Assign Materials

This analysis requires creating proper materials using the Engineering Data feature of Workbench. We will define a new Orthotropic Elastic material for the model as well as a Cohesive Zone Bilinear material for the Interface Delamination feature.

  1. In the Static Structural schematic, right-click the Engineering Data cell and choose Edit. The Engineering Data tab opens and displays Structural Steel as the default material.

  2. Click the box labeled "Click here to add new material" and enter the name "Interface Body Material".

  3. Expand the Linear Elastic option in the Toolbox and right-click Orthotropic Elasticity. Select Include Property. The required properties for the material are highlighted in yellow.

  4. Define the new material by entering the following property values and units of measure into the corresponding fields.

    PropertyValueUnit
    Young’s Modulus X Direction1.353E+05MPa
    Young’s Modulus Y Direction9000MPa
    Young’s Modulus Z Direction9000MPa
    Poisson’s Ratio XY0.24na
    Poisson’s Ratio YZ0.46na
    Poisson’s Ratio XZ0.24na
    Shear Modulus XY5200MPa
    Shear Modulus YZ0.0001MPa
    Shear Modulus XZ0.0001MPa

    The properties for the material should appear as follows:

  5. Click the box labeled "Click here to add new material" and enter the name “CZM Material”. This material will specify the formulation used to introduce the fracture mechanism (Cohesive Zone Material method).

  6. Expand the Cohesive Zone option in the Toolbox and right-click Exponential for Interface Delamination. Select Include Property. The required properties for the material are highlighted in yellow.

  7. Define the new material by entering the following property values and units of measure into the corresponding fields.

    PropertyValueUnit
    Maximum Normal Traction2.5E+07Pa
    Normal Separation Across the Interface4E-06m
    Shear Separation at Maximum Shear Traction1m

    The properties for the material should appear as follows.

7.4. Attach Geometry

  1. In the Static Structural schematic, right-click the Geometry cell and select Import Geometry>Browse.

  2. Browse to the proper location and open the file 2D_Fracture_Geom.agdb. This file is available here on the Ansys customer site.

  3. Right-click the Geometry cell and select Properties. In the Properties window, set the Analysis Type property to 2D.

    The Project Schematic should appear as follows:

  1. Launch Mechanical. Right-click the Model cell and then choose Edit. (Tip: You can also double-click the cell to launch Mechanical).

  2. Define unit system. From the Tools group on the Home tab, open the Units drop-down menu and select Metric (mm, kg, N, s, mV, mA).

7.5. Define 2D Behavior

  1. Define 2D behavior.

    1. Highlight the Geometry object.

    2. In the Details pane, specify the 2D Behavior property as Plane Strain. This constrains all of the UZ degrees of freedom. See the 2D Analyses section of the Mechanical User's Guide for additional information about this property.

  2. Apply Material: Expand the Geometry folder and select the Part 2 folder. Set the Assignment property to "Interface Body Material". Selecting the Part 2 folder allows you to assign the material to both parts at the same time.

  3. Suppress Contact.


    Caution:  Contact cannot be present for this analysis.


    1. Expand the Connections folder and then expand the Contacts folder.

    2. Right-click the Contact Region object and select Suppress.

7.6. Define Coordinate Systems

This analysis requires a mesh Match Control property to match the elements of the two parts. To properly define the Match Control property, you need to define coordinate systems for the element faces that will be matched with one another. In theory, for this model, one coordinate system could facilitate the specification of the Mesh Match Control because the coordinate systems you are about to create are virtually identical.

  1. Right-click the Coordinate Systems object in the tree and select Insert>Coordinate System.

  2. Right-click the new coordinate system object, select Rename, and name the object "High Coordinate System."

  3. In the Details pane of the newly-created Coordinate System object, select the Geometry property field Click to Change.

  4. Select the Edge selection filter (on the Graphics Toolbar) and highlight an edge in the center of the model.

    This tutorial employs the Depth Picking tool because of the close proximity of the two edges involved in the interface, as well as the crack. As illustrated here, the graphics window displays a stack of rectangles in the lower left corner. The rectangles are stacked in appearance, with the topmost rectangle representing the visible (selected) geometry and subsequent rectangles representing additional geometry selections. For this example, the topmost geometry is the "high" edge.

  5. Click Apply in the Geometry property. The "High Coordinate System" is defined.

  6. Right-click the Coordinate Systems object again and insert another Coordinate System object. Rename this object "Low Coordinate System."

  7. Select the Edge selection filter and highlight an edge in the center of the model. Using the Depth Picking tool, select the second rectangle in the stack, and then scope the edge as the geometry (Apply in the Geometry property). This scoping is illustrated below.

7.7. Define Mesh Options and Controls

  1. Select the Mesh object. Define the following Mesh object properties:

    • Set Element Order (Defaults category) to Quadratic (kept).

    • Enter an Element Size (Defaults category) of 0.75.

    • Set Use Adaptive Sizing (Sizing category) to No

  2. Right-click the Mesh object and select Insert>Match Control.

  3. Activate the High Geometry Selection property by selecting its field (that is highlighted in yellow). The Apply and Cancel buttons display. Select the Edge selection tool and highlight one of the edges in the center of the model. Use the Depth Picking tool to select the topmost geometry. Click the Apply button.

  4. Perform the same steps to specify the Low Geometry Selection property, as illustrated below.

  5. Change the Transformation property from Cyclic to Arbitrary and specify the High Coordinate System and Low Coordinate System properties using the coordinate systems created in the previous step of the tutorial. The object should appear as illustrated below.

  6. Select the Edge selection filter (on the Graphics Toolbar) and, holding the Ctrl key, select the four side edges.

  7. Right-click the Mesh object and select Insert>Sizing. This mesh sizing control should be scoped to the four side edges.

  8. In the Details view, enter 0.75 mm as the Element Size.

  9. Select the Edge selection filter (on the Graphics Toolbar) and highlight an edge in the center of the model. Use the Depth Picking tool and, holding the Ctrl key, select both rectangles in the lower left corner of the graphics window. Continue to hold the Ctrl key, and select an edge of the crack. Again, use the Depth Picking tool and select both rectangles in the lower left corner of the graphics window. Still holding the Ctrl key, select the top and bottom edges on the model.

  10. Right-click the Mesh object and select Insert>Sizing. This mesh sizing control should be scoped to six (top and bottom and the four interface edges) edges.

  11. In the Details view, enter 0.5 mm as the Element Size.

  12. Right-click the Mesh object and select Generate Mesh.

7.8. Define Interface Delamination Object

  1. Insert a Fracture folder into the tree by highlighting the Model object and selecting the Fracture option from the Define group on the Model Context tab.

  2. Select the Interface Delamination option from the Model group on the Fracture Context tab.

  3. In the Details pane, set the Method property to CZM.

  4. Set the Material property to CZM Material.

  5. Select the Match Control that was created earlier in the tutorial for the Match Control property.

    The Interface Delamination object is complete.

7.9. Configure the Analysis Settings

  1. Select the Analysis Settings object.

  2. Set the Auto Time Setting property to On and then enter 40 for the Initial Substeps, Minimum Substeps, and Maximum Substeps properties.

  3. In the Details pane, set the Large Deflection property to On to activate geometric nonlinearities.

7.10. Define Boundary Conditions

  1. Select the Edge selection filter and select the two edges on the side of the model that is opposite of the crack. Select one edge, press the Ctrl key, and then select the next edge.

  2. Highlight the Static Structural object and select the Fixed Support option from the Structural group on the Environment Context tab.

  3. Highlight the Static Structural object. With the Vertex selection filter active, select the vertex illustrated below, and select the Displacement option from the Structural group on the Environment Context tab.

  4. Highlight the Displacement object in the tree and enter 10 (mm in the positive Y direction) as the loading value for the Y Component property.

  5. Create another Displacement. With the Vertex selection filter active, select the bottom vertex and select the Displacement option from the Structural group on the Environment Context tab. Enter -10 (mm in the negative Y direction) as the loading value for the Y Component property.

7.11. Specify Result Objects and Solve

  1. Highlight the Solution object. From the Results group of the Solution Context tab, open the Deformation drop-down menu and select Directional.

  2. Under the Definition category in the Details view, set the Orientation property to Y Axis.

  3. From the Results group of the Solution Context tab, open the Probe drop-down menu and select Force Reaction.

  4. Select Displacement for the Boundary Condition property.

  5. Click the Solve button.

7.12. Review the Results

  1. Highlight the Directional Deformation and Force Reaction objects. Results appear as follows:

You may wish to validate results against those outlined in the verification test case (VM248 in the Ansys Mechanical APDL Verification Manual). This is most easily accomplished by creating User Defined Results using the Worksheet.

Congratulations! You have completed the tutorial.