*HEADING
: *DSLOAD AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [B21]
*NODE,NSET=ALL
1, 0.
2, 15.
*ELEMENT,TYPE=B21, ELSET=EALL
1,1,2
*surface, NAME=TOP
 1,SPOS
*surface, NAME=BOT
 1,SNEG
*BEAM GENERAL SECTION,SECTION=RECT,DENSITY=0.16667, ELSET=EALL
2., 3.
0., 0., -1.
3.E6,1.1538E6, .0001
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL, -10., 0.
*BOUNDARY
ALL,1,2
ALL,6
*STEP,PERTURBATION
*STATIC
*DSLOAD
 TOP, P, 2000.
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*EL PRINT,FREQ=0
*NODE PRINT
 RF, 
*END STEP
*STEP
*STATIC
*TEMPERATURE,OP=NEW
 ALL, 20., 6.6667
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*OUTPUT,FIELD
*ELEMENT OUTPUT
 SF, 
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*EL PRINT
 SF, 
*NODE PRINT
 RF, 
*END STEP
*STEP
*STATIC
*BOUNDARY,OP=NEW
 1,1,2
 1,6
*TEMPERATURE,OP=MOD
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
 U, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
 SE, 
*EL PRINT,FREQ=0
*NODE PRINT
 U, 
*END STEP