*HEADING
: *DSLOAD AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [S9R5]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
1    , 0.        , 0.
3    , 7.        , 0.
7    , 0.        , 7.
9    , 7.        , 7.
5, 3.5, 3.5
*NGEN,NSET=BOTTOM
1    , 3
*NGEN,NSET=TOP
7    , 9
*NGEN,NSET=LEFT
1    , 7    , 3
*NGEN,NSET=RIGHT
3    , 9    , 3
*NSET,NSET=ALL,GENERATE
1    , 9
*NSET,NSET=BOUNDARY,GENERATE
1,4
6,9
*ELEMENT,TYPE=S9R5, ELSET=EALL
1, 1, 3, 9, 7, 2, 6, 8, 4, 5
*surface, NAME=TOP
 1,SPOS
*surface, NAME=BOT
 1,SNEG
*SHELL SECTION,MATERIAL=A1, ELSET=EALL
2.        , 3
*MATERIAL,NAME=A1
*ELASTIC
3.E6      , 0.3
*EXPANSION
0.0001, 
*DENSITY
 1.0, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL, -10.
*BOUNDARY
ALL,1,6
*STEP,PERTURBATION
*STATIC
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*OUTPUT,FIELD
*ELEMENT OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*DSLOAD
BOT,HP,-3000.,7.,0.
*EL PRINT
 LOADS,
*OUTPUT,FIELD
*ELEMENT OUTPUT
LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
TOP,P,1000.
*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
BOTTOM,2
*TEMPERATURE,OP=MOD
ALL,40.
*END STEP