A heating plate is heated up through a heat flux. Parts on top of this plate can exchange heat via contacts. This is shown in this example. A component is in direct contact and exchanges heat via heat conductivity. Another component is very close to the heat source and exchanges heat via a simplified heat convection and radiation formulation. Another component is located outside the contact range and does not exchange heat.
*BOUNDARY_FLUX_SET *CONTACT_SURFACE_TO_SURFACE_THERMAL_ID *CONTROL_SOLUTION *CONTROL_TERMINATION *CONTROL_THERMAL_SOLVER *CONTROL_THERMAL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_TPRINT *DEFINE_CURVE *END *INCLUDE *INITIAL_TEMPERATURE_SET *KEYWORD *MAT_THERMAL_ISOTROPIC *PART *SECTION_SOLID *SET_PART_LIST *TITLE
*KEYWORD
$
$=============================CONTROL DEFINITIONS ==============================
$
*TITLE
Thermal contact between heating plate and parts
*CONTROL_SOLUTION
$ thermal analysis
$ soln nlq isnan lcint lcacc ncdcf
1
*CONTROL_THERMAL_SOLVER
$ transient thermal analysis
$ atype ptype solver cgtol gpt eqheat fwork sbc
1 1 11
$ msglvl maxitr abstol reltol omega tsf
*CONTROL_THERMAL_TIMESTEP
$ ts tip its tmin tmax dtemp tscp lcts
1 0.1
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas nosol
10
$
$================================ OUTPUTDATA ===================================
$
*DATABASE_BINARY_D3PLOT
$ dt
0.1
*DATABASE_GLSTAT
$ dt
0.1
*DATABASE_MATSUM
$ dt
0.1
*DATABASE_TPRINT
$ dt
0.1
$
$============================== PART DEFINITIONS ===============================
$
*PART
$ title
heating_plate
$ pid secid mid eosid hgid grav adpopt tmid
1 1 0 0 0 0 0 2
*PART
$ title
cube_contact
$ pid secid mid eosid hgid grav adpopt tmid
2 1 0 0 0 0 0 1
*PART
$ title
cube_close
$ pid secid mid eosid hgid grav adpopt tmid
3 1 0 0 0 0 0 1
*PART
$ title
cube_far
$ pid secid mid eosid hgid grav adpopt tmid
4 1 0 0 0 0 0 1
$
$============================= SECTION PROPERTIES ==============================
$
*SECTION_SOLID
$ secid elform
1 -1
$
$======================== THERMAL MATERIAL PROPERTIES ==========================
$
*MAT_THERMAL_ISOTROPIC
$ zinc
$ tmid tro tgrlc tgmult tlat hlat
1 7.13e-9
$ hc tc
3.88e8 110.0
*MAT_THERMAL_ISOTROPIC
$ steel
$ tmid tro tgrlc tgmult tlat hlat
2 7.85e-9
$ hc tc
4.60e8 40.0
$
$=======================NODE/ELEMENT/SET/SEGMENT DEFINTIONS ====================
$
*INCLUDE
06_heating_plate_parts.k
$
$*SET_SEGMENT_TITLE segments_plate_bottom sid 1
$
*SET_PART_LIST
$ sid
5
$ pid1 pid2 pid3
2 3 4
$
$=========================== CONTACT DEFINITIONS ===============================
$
*CONTACT_SURFACE_TO_SURFACE_THERMAL_ID
1,contact
$ ssid msid sstyp mstyp sboxid mboxid spr mpr
5 1 2 3
$ fs fd dc vc vdc penchk bt dt
$ sfs sfm sst mst sfst sfmt fsf vsf
$ k frad h0 lmin lmax ftoslv
0.1 6.3E-12 50.0 0.01 2.0
$ frad based on emissivity e1=e2=0.2
$==================== THERMAL BOUNDARY and INITIAL CONDITIONS ==================
$
*INITIAL_TEMPERATURE_SET
$ nid/sid temp loc
0 293.15
*BOUNDARY_FLUX_SET
$ ssid
1
$ lcid mlc1 mlc2 mlc3 mlc4 loc nhisv
101 -1.0 -1.0 -1.0 -1.0
*DEFINE_CURVE
$ lcid sidr sfa sfo offa offo dattyp
101 0 1 1 0 0 0
$ time , heat flux
0.0,0.0
1.0,200.0
11.0,200.0
$
*END
