VM23

VM23
Thermal-structural Contact of Two Bodies

Overview

Reference: Any Basic Mechanics Text
Analysis Type(s): Static Analysis (ANTYPE = 0)
Element Type(s):
2D Coupled-Field Solid Element (PLANE13)
2D Coupled-Field Solid Element (PLANE223)
2D/3D Node-to-Surface Contact Element (CONTA175)
2D Target Segment Element (TARGE169)
Input Listing: vm23.dat

Test Case

Two bodies, A and B, are initially at a temperature of 100°C. A temperature of 500°C is then imposed at the left edge, EA1, of A. Further, the right edge, EB2, of B is heated to attain a temperature of 850°C and is subsequently cooled back to 100°C. Compute the interface temperature (right edge) of A, EA2, and the amount of heat flow through the interface when the right edge of EB2 is at 600°C and 850°C, respectively. Also, compute the heat flow through the interface when B is subsequently cooled to 100°C.

Figure 29: Contact Problem Sketch

Contact Problem Sketch

Material PropertiesGeometric PropertiesLoading
E = 10 x 106 N/m2
K = 250 W/m°C
α = 12 x 10-6 m/m°C
Contact Conductance = 100 W/°C (per contact element)
h = 0.1 m
1 = 0.4 m
2 = 0.5 m
δ = 0.0035 m
Tref = 100°C
Load Step 1 T@E1 = 500°C
Load Step 2 T@E2 = 850°C
Load Step 3 T@E2 = 100°C

Analysis Assumptions and Modeling Notes

A coupled-field analysis is performed to solve this thermal/structural contact problem using PLANE13 elements and PLANE223 elements with weak coupling. The interface is modeled using CONTA175 and TARGE169 elements with the heat conduction option. For PLANE13 elements under large deflection (NLGEOM,ON), the thermal results are not updated during solution.

Results Comparison

 TargetMechanical APDLRatio
Using PLANE13 Elements
T @ EB2 = 600°CT @ EA2 (°C)539.0539.01.000
Q W2439.02439.01.000
T @ EB2 = 850°CT @ EA2 (°C)636.6636.61.000
Q W8536.68536.61.000
T @ EB2 = 100°CQ W0.00.0--
Using PLANE223 Elements
T @ EB2 = 600°CT @ EA2 (°C)539.0539.01.000
Q W2450.62450.61.000
T @ EB2 = 850°CT @ EA2 (°C)636.6636.61.000
Q W8588.78588.71.000
T @ EB2 = 100°CQ W0.00.0--