*HEADING
 THERMALLY COUPLED ANALYSIS OF A SYMMETRIC DISC BRAKE
*NODE
*NODE, NSET=CENTER
*NGEN, LINE=C, NSET=SUPPORT
*NGEN, LINE=C, NSET=INNER1
*NGEN, LINE=C, NSET=INNER2
*NSET, NSET=INNER
INNER1, INNER2
*NGEN, LINE=C, NSET=OUTER
*NFILL, NSET=DISK
INNER, OUTER, 3, 100
*NSET, NSET=SYMM, GENERATE
*ELEMENT, TYPE=C3D8T, ELSET=SUPPORT
*ELEMENT, TYPE=C3D8T, ELSET=DISK
*ELGEN, ELSET=SUPPORT
*ELGEN, ELSET=DISK
*ELSET, ELSET=FULLDISK
SUPPORT, DISK
*ELSET, ELSET=DISKSURF, GENERATE
*SOLID SECTION, ELSET=FULLDISK, MATERIAL=STEEL, ORIENTATION=CYL
*MATERIAL, NAME=STEEL
*ELASTIC
210.E9, .3
*EXPANSION
11.E-6,
*DENSITY
7800.,
*SPECIFIC HEAT
452.,
*CONDUCTIVITY
48.,
*ORIENTATION, SYSTEM=CYLINDRICAL, NAME=CYL
*MPC
BEAM, SUPPORT, 1000
*NODE
*NGEN, LINE=C, NSET=INPAD
*NGEN, LINE=C, NSET=OUTPAD
*NFILL, NSET=NPAD
INPAD, OUTPAD, 3, 100
*NSET,NSET=QA_TEST_NPAD
*NSET, NSET=NPADBACK, GENERATE
*NSET, NSET=ALL
DISK, NPAD, SUPPORT
*ELEMENT, TYPE=C3D8T
*ELGEN, ELSET=EPAD
*ELSET, ELSET=EPADBACK, GENERATE
*SOLID SECTION, ELSET=EPAD, MATERIAL=PADMAT, ORIENTATION=CYL
*MATERIAL, NAME=PADMAT
*ELASTIC
2200.E6, .25,  20.0
1300.E6, .25, 100.0
 530.E6, .25, 200.0
 320.E6, .25, 300.0
*EXPANSION
10.E-6, 20.0
30.E-6,200.0
*DENSITY
1550.,
*SPECIFIC HEAT
1200.,
*CONDUCTIVITY
0.9,
**
** DEFINE SURFACES AND CONTACT PAIR
**
*surface, NAME=MASTER
DISKSURF, S2
*surface, NAME=SLAVE
10001, S1
10002, S1
10003, S1
10004, S1
10101, S1
10102, S1
10103, S1
10104, S1
10201, S1
10202, S1
10203, S1
10204, S1
*CONTACT PAIR, INTERACTION=INTERF
SLAVE,MASTER
*SURFACE INTERACTION, NAME=INTERF
*FRICTION
0.37,,,20.0
0.38,,,100.0
0.41,,,200.0
0.39,,,300.0
0.24,,,400.0
*GAP CONDUCTANCE
*GAP HEAT GENERATION
*INITIAL CONDITIONS, TYPE=TEMPERATURE
ALL,300.
*RESTART, WRITE, FREQUENCY=100
*****
*STEP, NLGEOM, AMPLITUDE=RAMP, UNSYMM=YES
 PRESS THE PAD AGAINST THE DISK
*COUPLED TEMPERATURE-DISPLACEMENT,DELTMX=100.
0.0005, 0.001
*BOUNDARY
SYMM, 3, 3
CENTER, 1, 6
NPADBACK, 1, 2
*DLOAD
EPADBACK, P2, 1.7E6
*FILM
DISKSURF, F2, 20., 100.
**
*PRINT, CONTACT=YES
*CONTACT PRINT, FREQUENCY=5
*CONTACT FILE, FREQUENCY=5
*EL PRINT, FREQUENCY=100
S,
E,
*NODE PRINT, NSET=CENTER, FREQUENCY=1
UR3,RF
*NODE PRINT, NSET=NPAD, FREQUENCY=100
U1, U2, U3, NT11, RF1, RF2, RF3
*NODE FILE, NSET=CENTER, FREQUENCY=1
U, RF
*NODE FILE, NSET=NPAD, FREQUENCY=100
U, NT, RF
*OUTPUT,FIELD,Frequency=5
*CONTACT OUTPUT
CSTRESS,CDISP,PPRESS
*OUTPUT,FIELD,Frequency=5
*CONTACT OUTPUT, NSET= QA_TEST_NPAD
CSTRESS,CDISP,PPRESS
*OUTPUT,FIELD,Frequency=1
*NODE OUTPUT,NSET=CENTER
U, RF
*OUTPUT,FIELD,Frequency=100
*NODE OUTPUT,NSET=NPAD
U, NT, RF
*END STEP
*****
*STEP,NLGEOM,INC=100,AMPLITUDE=RAMP,UNSYMM=YES,EXTRAPOLATION=PARABOLIC
 ROTATE THE DISK BY 60 DEGREES
*COUPLED TEMPERATURE-DISPLACEMENT, DELTMX=100.
0.001, 0.015, , .004
*BOUNDARY
CENTER, 6, 6, 1.04717
*END STEP