*HEADING
: *DSLOAD AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [STRI65]
*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=STRI65, ELSET=EALL
    1,     1,     3,     9,     2,     6,     5
50001,     1,     9,     7,     5,     8,     4
*surface, NAME=TOP
 EALL,SPOS
*surface, NAME=BOT
 EALL,SNEG
*ELSET, ELSET=E1
1,
*SHELL GENERAL SECTION, ELSET=EALL, DENSITY=2.0
 6.59341E+06,1.97802E+06,6.59341E+06,0.0,0.0,2.30769E+06,0.0,0.0
 0.0,2.19780E+06,0.0,0.0,0.0,6.59341E+05,2.19780E+06,0.0
 0.0,0.0,0.0,0.0,7.69231E+05
 857.14,857.14
 1.0, 1.0,-10.0
*TRANSVERSE SHEAR STIFFNESS
 1.92308E+06,1.92308E+06
*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
TOP,HP,3000.,7.,0.
*EL PRINT, TOTALS=YES
1,2,3
 LOADS,
 IVOL,
 SVOL,
 EVOL,
*OUTPUT,FIELD
*ELEMENT OUTPUT
LOADS,
IVOL,
SVOL,
EVOL,
*NODE PRINT,TOTAL=YES
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
TOP,P,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
ALL,40.
*OUTPUT,FIELD
*ELEMENT OUTPUT
SF,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
U,
*EL PRINT
 SF, 
*NODE PRINT,TOTAL=YES
 RF,
 U,
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
BOUNDARY,3,6
LEFT,1
BOTTOM,2
*TEMPERATURE,OP=MOD
ALL,40.
*END STEP