*HEADING
EB32QXD1: DLOADS AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [B31]
*PREPRINT,MODEL=YES
*NODE,NSET=ALL
1, 0.
2, 15.
*NODE, NSET=NEWNODES
3, 20.0
4, 21.0
5, 22.0
6, 23.0
*ELEMENT,TYPE=B31, ELSET=EALL
1, 1, 2
*ELEMENT,TYPE=B31, ELSET=BSECPOIS
2,3,4
*ELEMENT,TYPE=B31, ELSET=BGSECPOIS
3,5,6
*BEAM GENERAL SECTION,SECTION=RECT,DENSITY=0.16667, ELSET=EALL
2. , 3.
0. , 0. , -1.
3.E6, 1.1538E6, .0001
*BEAM SECTION,SECTION=RECT,MAT=MAT2,POISSON=-0.5 ,ELSET=BSECPOIS
0.1, 0.1
*BEAM GENERAL SECTION,SECTION=RECT,POISSON=-0.5, ELSET=BGSECPOIS, DEP=1
0.1, 0.1
0., 0., -1.
10.E6,11538462.0,0.0,,0.0
30.E6,11538462.0,0.0,,1.0
10.E6,11538462.0,0.0,,2.0
30.E6,11538462.0,0.0,,3.0
*SECTION POINTS
0.0,0.05, 0.0,-0.05
*MATERIAL, NAME=MAT2
*ELASTIC
30.E6,0.3
*EXPANSION
0.0,
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL, -10., 0.
*INITIAL CONDITIONS,TYPE=FIELD, VAR=1
 ALL, 1.0
NEWNODES,1.0
*BOUNDARY
ALL,1,6
NEWNODES,1,6
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, CENT, 3000., 7.5, 0., 0.,  0., 0., 1.
*EL PRINT
 LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
1,CENTRIF,3000.,7.5,0.,0.,0.,0.,1.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, GRAV, 1000., 1., 0., 0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, PX, 1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, PY, 2000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, PZ, 1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, P1, 1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, P2, 2000.
*END STEP
*STEP
*STATIC
*TEMPERATURE,OP=NEW
 ALL, 20., 6.6667
*EL PRINT
 SF, 
*NODE PRINT
 RF, 
*END STEP
*STEP
*STATIC
*BOUNDARY,OP=NEW
1,1,6
3,1,6
5,1,6
*TEMPERATURE,OP=MOD
*EL PRINT,FREQ=0
*NODE PRINT
 U, 
*END STEP
*STEP, NLGEOM
*STATIC
0.5,1.0
*MODEL CHANGE, REMOVE
EALL,
*BOUNDARY, OP=NEW
3,1,6
5,1,6
*CLOAD
 4,1, 1000.0
 6,1, 1000.0
*EL PRINT,FREQ=1,ELSET=BSECPOIS
S,
E,
*EL PRINT,FREQ=1,ELSET=BGSECPOIS
S,
E,
*NODE PRINT,FREQ=1, NSET=NEWNODES
 U, 
*END STEP