Choosing Penalty vs. Augmented Lagrange in ANSYS Mechanical Contact Modeling

Desired Outcome

Accurately enforce bonded or frictional interfaces in a finite element analysis while maintaining numerical stability and reasonable memory usage.

Prerequisites

• ANSYS Mechanical (2025 R1 or later) installed on a machine with sufficient RAM for the model size.
• Basic familiarity with creating geometry, meshing, and applying loads in the ANSYS Workbench environment.
• Access to the Mechanical project and the ability to edit contact settings.

Procedure

1. Define the Contact Pair
   • In the Model tree, right‑click Contact → New Contact Pair.
   • Choose the Surface to Surface type for two distinct parts.
   • Assign the contacting surfaces (e.g., BlockA.Top and BlockB.Bottom).
   • Set the Contact Formulation to either Penalty or Augmented Lagrange based on the decision table below.
2. Configure Penalty Parameters (if selected)
   • In the Contact Pair properties, open the Penalty tab.
   • Specify a Contact Stiffness (e.g., 1e6 for a typical steel interface). The value should be high enough to prevent penetration but not so high that it causes ill‑conditioning.
   • Optionally, enable Automatic Penalty Stiffness for a quick estimate, but review the suggested value.
3. Configure Augmented Lagrange Parameters (if selected)
   • Open the Augmented Lagrange tab in the Contact Pair properties.
   • Define the Multiplier Field (e.g., Multiplier).
   • Set the Penalty Coefficient (e.g., 1e5) to control the conditioning of the multiplier equations.
   • Choose the Multiplier Initialization method (default is zero).
4. Mesh Refinement Near Interfaces
   • Right‑click the mesh of each part → Insert Mesh Control.
   • Apply a Refinement to the contacting surfaces with a factor of 2–3, ensuring element size ≤ 1/10 of the expected contact gap.
   • Generate the mesh.
5. Apply Loads and Boundary Conditions
   • Apply a known load (e.g., a compressive force of 100 kN on BlockA top).
   • Fix the opposite side of BlockB to prevent rigid body motion.
6. Run the Analysis
   • Set the solver to Static, General or Nonlinear as appropriate.
   • Click Solve and monitor the Solution Status panel.

Expected Checks

• Contact Gap: In the post‑processing Result tree, view Contact Gap. Acceptable penetration is ≤ 1e-6 m for most engineering cases.
• Reaction Forces: Verify that the sum of reaction forces equals the applied load within a tolerance of 0.5 %.
• Solver Residuals: Residuals should drop below 1e-6 for the last iteration. If not, consider adjusting the penalty or multiplier values.
• Iteration Count: Excessive iterations (> 20 for a linear problem) may indicate ill‑conditioning.

Recovery Options

1. Adjust Penalty Stiffness
   • If penetration is too large, increase the penalty stiffness by a factor of 10.
   • If residuals are high, reduce it gradually.
2. Refine Mesh Near Contact
   • Increase the refinement factor or add an additional Mesh Control element layer.
   • Re‑run the analysis to observe changes in contact gap and residuals.
3. Switch Formulation
   • If the penalty method leads to poor conditioning, switch to Augmented Lagrange and re‑run.
   • Conversely, if memory usage is high with Augmented Lagrange, revert to penalty with a moderate stiffness.
4. Check Contact Surface Quality
   • Ensure surfaces are planar and aligned; mis‑aligned meshes can create artificial gaps.
   • Use the Contact Check tool to visualize potential penetration.

Practical Example: Two‑Block Compression Test

Below is a concise APDL script that demonstrates the key steps for a two‑block compression test using both formulations. Replace placeholder names with your actual part IDs.

! Create parts
/Preproc
ET, 1, SOLID185
! Define geometry for BlockA and BlockB
! (Geometry creation omitted for brevity)
!
! Mesh the parts with refinement at contact surfaces
MESH, ALL
!
! Define contact pair
CP, 1, 1, 2
! For Penalty formulation
CP, 1, PENALTY, 1e6
!
! For Augmented Lagrange formulation
! CP, 1, AUGLAG, 1e5
!
! Apply loads and boundary conditions
D, 1, ALL, 0
S, 1, FORCE, 100e3
!
! Solve
SOLVE
FINISH

After running, use the post‑processing commands PRINT, GAP, 1 and PRINT, REAC, 1 to extract the contact gap and reaction forces for verification.

Decision Table: Penalty vs. Augmented Lagrange

Aspect	Penalty	Augmented Lagrange
Constraint Enforcement	Approximate; depends on stiffness	Exact (within solver tolerance)
Numerical Conditioning	Sensitive to stiffness; high values cause ill‑conditioning	Better conditioning due to multiplier field
Memory Footprint	Lower (no additional DOFs)	Higher (extra multiplier DOFs)
Implementation Simplicity	Very simple	Requires setting multiplier field
Best For	Large models where memory is a concern and contact gaps are small	High‑accuracy contact problems, especially where penetration must be strictly zero

Summary

Choosing between penalty and Augmented Lagrange in ANSYS Mechanical hinges on the trade‑off between enforcement accuracy and computational resources. Start with a modest penalty stiffness and a coarse mesh to gauge convergence. If penetration persists or residuals remain high, switch to Augmented Lagrange or refine the mesh. Always validate the contact gap and reaction forces against expected physics before finalizing the model.