# Contact-based Thermal Analysis

In OptiStruct, structural models involving contact are solved by using Small Displacement Nonlinear Analysis.

The analysis involves finding the contact status, such as contact clearance and pressure. Contact clearance spans the distance between the main and secondary, while contact pressure is developed between two surfaces in contact.

Thermal Contact Analysis via PCONTHT and PGAPHT is supported for Linear Steady-State Heat Transfer, Linear Transient Heat Transfer, Nonlinear Steady-State Heat Transfer, and Nonlinear Transient Heat Transfer analyses. Thermal Contact is also supported for One-Step Thermal Transient Stress Analysis (OSTTS).

In Figure 1, you can see that a change in contact status does not affect the thermal problem. This may lead to inaccurate solutions if thermal conductance depends on the contact status. In Figure 2, the contact clearance and/or pressure changes during the course of the quasi-static nonlinear analysis, the corresponding change in the thermal conductance will affect the solution of the thermal problem.

## Thermal Contact Solver

Thermal analysis is performed first using initial contact status.

An iterative solution process is developed to solve fully coupled nonlinear thermal structural problem, as shown in Figure 3. Nonlinear structural analysis is employed to find contact status. Thermal conductance at the contact interface is calculated based on contact clearance or pressure, or based on user-defined values. Coupling is essential because the contact status is used to determine thermal conductance. Temperature results from thermal analysis are used as convergence criteria.

- For thermal contact problems with
CGAP/CGAPG,
PGAPHT is required. The PGAPHT
entry should have the same
`PID`as PGAP. For problems with CONTACT and PCONT, the PCONTHT entry should be used and it requires the same`PID`as PCONT. - Thermal conductivity based on the AUTO option
(
`KCHTC`/`KCHT`fields on PCONTHT/PGAPHT entries) can be used in thermal analysis to allow OptiStruct to automatically determine the conductivity values based on the conductivity of surrounding elements. AUTO conductivity is chosen in such a way that it works as a perfect conductor for closed gap/contact and insulator for open gap/contact. - For problems without PCONT, PCONTHT is not required. Thermal conductivity is internally calculated by the AUTO method.
- High conductance at the interface is automatically enforced for FREEZE contact.
- In Heat transfer analysis with Contact (except for
FREEZE contact), the secondary set ID
(
`SSID`) should not be defined as a GRID set.

Theoretically, while higher conductance values enforce a perfect conductor, excessively high values may cause poor conditioning of the conductivity matrix. If such effects are observed, it may be beneficial to reduce the value of conductance, or use conductance based contact clearance and pressure.

### Clearance based Thermal Conductance (TCID on PCONTHT/PGAPHT, via TABLED#)

### Pressure based Thermal Conductance (TPID on PCONTHT, via TABLED#)

### Clearance and Pressure based Thermal Conductance (TCID and TPID on PCONTHT, via TABLED#)

### Thermal Contact with FREEZE Status

For thermal contact with FREEZE status, the actual contact status (open or closed) based on geometry will be used in heat transfer analysis.

## Thermal Contact without Static Analysis

Pure heat transfer analysis with thermal contact is solved based on initial contact status. Contact clearance and area are calculated based on geometry.

`KAHT`, `KBHT` and `TCID` on
PGAPHT Bulk Data Entry and
`KCHTC`, `KOHTC` and
`TCID` on PCONTHT Bulk Data
Entry can be used. Contact pressure is not available without static
analysis. Therefore, `TPID` can not be used in such a
scenerio.