*HEADING
: BOX FRAME SECTION, THERMAL LOADING  [FRAME2D]
*REST,W 
*PREPRINT,MODEL=YES
*NODE
 1,  0.0, 0.0
11, 75., 0.0
*NGEN,NSET=ALL
 1,11,2
*ELEMENT,TYPE=FRAME2D, ELSET=EALL
1,1,3
2,3,5
3,5,7
4,7,9
5,9,11
*ELSET, ELSET=FOUR
4, 
*FRAME SECTION,SECTION=BOX, ELSET=EALL, DENSITY=8000.0E-4, ZERO=-10., DEPENDENCIES=1
4., 2., .2, .1, .2, .1
0., 0., -1.
3.E6, 1.1538E6, .001E-4, -10., 0.
0.5E6, 0.5769E6, .003E-4, 90., 0.
3.E6, 1.1538E6, .001E-4, -10., 1.
1.5E6, 0.5769E6, .002E-4, 90., 1.
*INITIAL CONDITIONS, TYPE=TEMPERATURE
ALL, -10. 
*INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1
ALL, 1.
*BOUNDARY
1, 1, 2
1, 6
*********************** STEP 1
*STEP
*STATIC
*CLOAD
11,2,25.E-4
*EL PRINT
SF,
SEE,
*EL FILE, ELSET=FOUR
SF, 
SEE,
*NODE FILE
 U, COORD
 RF, 
*NODE PRINT
 U,
 RF, 
*ENERGY FILE
*ENERGY PRINT
*OUTPUT, FIELD, FREQ=99, VAR=ALL 
*OUTPUT, HISTORY, FREQ=99, VAR=PRESELECT 
*ENERGY OUTPUT, VAR=PRESELECT  
*END STEP
*********************** STEP 2
*STEP
*STATIC
*DLOAD
EALL, GRAV,10.0E-4, 0., -1., 0.
*END STEP
*********************** STEP 3
*STEP,PERTURBATION
*STATIC
*DLOAD
EALL, PY, -12.16E-4
*EL PRINT
SF,
SEE,
*EL FILE, ELSET=FOUR
SF, 
SEE,
*NODE FILE
 U, COORD
 RF , 
*NODE PRINT
 U,
 RF, 
*ENERGY FILE
*ENERGY PRINT
*OUTPUT, FIELD, FREQ=99, VAR=PRESELECT
*OUTPUT, HISTORY, FREQ=99, VAR=ALL
*END STEP
*********************** STEP 4
*STEP
*STATIC
*DLOAD
EALL, P2, -12.16E-4
*END STEP
*********************** STEP 5
*STEP
thermal loading. 
 the cantilever tip is also fixed.
*STATIC
*CLOAD, OP=NEW
*DLOAD, OP=NEW
*TEMPERATURE, OP=NEW 
ALL, 90.
*FIELD
ALL, 1.
*BOUNDARY
11, 1, 2
11, 6
*EL PRINT
SF,
SEE,
TEMP, 
THE, 
*EL FILE, ELSET=FOUR
SF, 
SEE,
THE, 
*NODE FILE
 U, COORD,
 RF, 
*NODE PRINT
 U,
 RF,
*ENERGY FILE
*ENERGY PRINT 
*OUTPUT, FIELD, FREQ=99, OP=ADD
*ELEMENT OUTPUT, ELSET=EALL
SF, SEE, THE
*NODE OUTPUT, NSET=ALL
U, RF, 
*OUTPUT, HISTORY, FREQ=99
*ELEMENT OUTPUT, ELSET=EALL
SF, SEE, 
*NODE OUTPUT, NSET=ALL
U, RF, 
*ENERGY OUTPUT, VAR=ALL  
*END STEP
*********************** STEP 6
*STEP, PERTURBATION
 thermal loading. (2)
 the cantilever tip is also fixed.
 Note that the derivative of the thermal strain with respect
 to the temperature at 90 degree temperature is used as
 the thermal expansion coefficient for this perturbation step.
 Refer to the equation on page 12.1.2-1 of v5.7 /Standard
 User's manual for computing this value.  It turns out to be
 0.003E-4 and not 0.002E-4.
*STATIC
*TEMPERATURE
ALL, 100.
*EL PRINT
SF,
SEE,
THE, 
*EL FILE, ELSET=FOUR
SF,
SEE,
THE, 
*NODE FILE
 U, COORD,
 RF, 
*NODE PRINT
 U,
 RF,
 NT, 
*ENERGY FILE
*ENERGY PRINT
*OUTPUT, FIELD, FREQ=99, OP=ADD
*ELEMENT OUTPUT, ELSET=EALL
SF, SEE, 
*NODE OUTPUT, NSET=ALL
U, RF, 
*OUTPUT, HISTORY, FREQ=99, OP=REPLACE
*ELEMENT OUTPUT, ELSET=EALL
SF, SEE, THE
*NODE OUTPUT, NSET=ALL
U, RF, 
*ENERGY OUTPUT
ALLWK, ALLPD, ALLSE, ETOTAL
*END STEP