*HEADING
: DLOADS AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [B21]
*NODE,NSET=ALL
*NODE, NSET=NEWNODES
*ELEMENT,TYPE=B21, ELSET=EALL
*ELEMENT,TYPE=B21, ELSET=BSECPOIS
*ELEMENT,TYPE=B21, ELSET=BGSECPOIS
*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
0.1, 0.1
0., 0., -1.
30.E6,11538462.0,0.0
*SECTION POINTS
*MATERIAL, NAME=MAT2
*ELASTIC
30.E6,0.3
*EXPANSION
0.0,
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL, -10., 0.
*BOUNDARY
ALL,1,2
ALL,6
NEWNODES,1,2
NEWNODES,6
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, CENT, 3000., 7.5, 0., 0.,  0., 0., 1.
*EL PRINT,FREQ=0
*NODE PRINT,NSET=ALL
 RF,
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, CENTRIF, 3000., 7.5, 0., 0.,  0., 0., 1.
*EL PRINT,FREQ=0
*NODE PRINT,NSET=ALL
 RF,
*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, P2, 2000.
*END STEP
*STEP
*STATIC
*TEMPERATURE,OP=NEW
 ALL, 20., 6.6667
*EL PRINT, ELSET=EALL
 SF,
*NODE PRINT, NSET=ALL
 RF,
*END STEP
*STEP
*STATIC
*BOUNDARY,OP=NEW
 1,1,2
 1,6
*TEMPERATURE,OP=MOD
*EL PRINT,FREQ=0
*NODE PRINT, NSET=ALL
 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