Problem Description:
I am simulating a tool assembly (Extrusion Die and Head) made of dissimilar materials. The goal is to observe how differential thermal expansion causes the initial contact interface to separate (opening a gap), and how this separation affects the heat transfer between the two components.
Model Setup:
-
Physics: Coupled Thermomechanical (Static Structural + Heat Transfer).
-
Initial Condition: Components are in contact at t=0.
-
Loading: Internal pressure ramp and thermal loading to operating temperature.
-
Contact: Physical Contact with “Enable Heat Transfer” set to Yes.
Current Issues:
-
Singularity & Stability: As soon as the components begin to separate due to thermal expansion, the “floating” component (the Die) loses its primary constraint, leading to rigid body motion errors and solver crashes (Parallel Process/Matrix Singular errors).
-
Thermal Discontinuity: The solver struggles with the sudden drop in thermal conductance as the contact transitions from “Solid-to-Solid” to a “Gap” (Air/Vacuum).
-
Hardware Constraints: The high mesh density required to resolve the contact interface (approx. 0.01m on a 40cm part) combined with 2nd-order elements is causing “Out of Memory” errors on standard computing nodes.
Steps Taken So Far:
-
Implemented Automatic Time Stepping to capture the non-linear contact transition.
-
Applied Load Ramping via tables for both pressure and temperature to improve convergence.
-
Added Surface Refinement on the contact faces.
-
Ensured no conflicting “Bonded” contacts exist on the same faces.
Seeking Advice On:
-
Stabilization: Best practices for using Elastic Supports (artificial springs) to prevent rigid body motion without over-constraining the thermal expansion.
-
Contact Stiffness: Recommendations for scaling contact stiffness (Penalty Factor) to ease convergence in metal-to-metal separation.
-
Thermal Gap Conductance: How to realistically model the thermal transition when the gap opens to prevent numerical oscillation.