*HEADING
: *DSLOAD AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [S3R]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
1, 0., 0.
3, 7., 0.
7, 0., 7.
9, 7., 7.
*NSET,NSET=LEFT
1, 7
*NSET,NSET=BOT
1, 3
*ELEMENT,ELSET=ELEMS,TYPE=S3R
1, 1,3,9
2, 1,9,7
*surface, NAME=TOP1
 1,SPOS
*surface, NAME=TOP2
 2,SPOS
*SHELL GENERAL SECTION,MATERIAL=A1, ELSET=ELEMS
2.,  3
*MATERIAL,NAME=A1
*ELASTIC
3.E6      , 0.3
*EXPANSION
0.0001, 
*DENSITY
 1., 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL  , -10.
*BOUNDARY
ALL  , 1    , 6
*STEP,PERTURBATION
*STATIC
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*DSLOAD
 TOP1,P,1000.
 TOP2,P,1000.
*EL PRINT
 LOADS,
*OUTPUT,FIELD
*ELEMENT OUTPUT
LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 TOP1,HP,3000.,7.,0.
 TOP2,HP,3000.,7.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
 ALL,40.
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
U,
*EL PRINT
 S, 
*NODE PRINT
 RF,
 U,
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
 ALL,3,6
 LEFT,1
 BOT,2
*TEMPERATURE,OP=MOD
 ALL,40.
*END STEP