*HEADING
: DLOADS AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [S4]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
1    , 0.        , 0.
3    , 7.        , 0.
7    , 0.        , 7.
9    , 7.        , 7.
*NODE, NSET=NNEW
11, 0.0, 0.0, 1.0
12, 1.0, 0.0, 1.0
13, 1.0, 1.0, 1.0
14, 0.0, 1.0, 1.0
15, 0.0, 0.0, 2.0
16, 1.0, 0.0, 2.0
17, 1.0, 1.0, 2.0
18, 0.0, 1.0, 2.0
*NSET,NSET=LEFT
1    , 7
*NSET,NSET=BOT
1    , 3
*ELEMENT,TYPE=S4, ELSET=EALL
1    , 1    , 3    , 9    , 7
*ELEMENT,TYPE=S4, ELSET=SSEC
2    , 11    , 12    , 13   , 14
*ELEMENT,TYPE=S4, ELSET=SGSEC
3    , 15    , 16    , 17   , 18
*SHELL GENERAL SECTION,MATERIAL=A1, ELSET=EALL
2. ,5
*SHELL SECTION,MATERIAL=MAT2, ELSET=SSEC, POISSON=-0.2
0.01,
*SHELL GENERAL SECTION,MATERIAL=MAT2, ELSET=SGSEC, POISSON=-0.2
0.01,
*MATERIAL,NAME=A1
*ELASTIC
3.E6      , 0.3
*EXPANSION
0.0001, 
*DENSITY
 1.0, 
*MATERIAL,NAME=MAT2
*ELASTIC
30.E6      , 0.3
*EXPANSION
0.000, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL, -10.
*BOUNDARY
ALL, 1,5
11, 1,6
14, 1,6
15, 1,6
18, 1,6
*STEP,PERTURB
*STATIC
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*DLOAD
1,BX,1000.
*EL PRINT, TOTALS=YES
1,2,3,4,5
 LOADS,
 IVOL,
 SVOL,
 EVOL,
*OUTPUT,FIELD
*ELEMENT OUTPUT
LOADS,
IVOL,
SVOL,
EVOL,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURB
*STATIC
*DLOAD
1,BY,2000.
*END STEP
*STEP,PERTURB
*STATIC
*DLOAD
 1,BZ,3000.
*END STEP
*STEP,PERTURB
*STATIC
*DLOAD
 1,GRAV,3000.,0,0,1.
*END STEP
*STEP,PERTURB
*STATIC
*DLOAD
1,CENT,3000.,0.,0.,0.,0.,1.,0.
*END STEP
*STEP,PERTURB
*STATIC
*DLOAD
1,HP,3000.,7.,0.
*END STEP
*STEP,PERTURB
*STATIC
*DLOAD
1,P,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
ALL,40.
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
MISESMAX,
*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
*STEP, NLGEOM
*STATIC
0.5,1.0
*MODEL CHANGE, REMOVE
EALL,
*BOUNDARY, OP=NEW
11, 1,6
14, 1,6
15, 1,6
18, 1,6
*CLOAD
12, 1, 10000.0
13, 1, 10000.0
16, 1, 10000.0
17, 1, 10000.0
*EL PRINT, ELSET=SSEC
S,
LE,
STH
*EL PRINT, ELSET=SGSEC
S,
E,
STH
*NODE PRINT, NSET=NNEW
U,
RF,
*OUTPUT,FIELD,OP=ADD
*ELEMENT OUTPUT
MISESMAX,
*END STEP