*HEADING
EC3RSFDI: *DSLOAD AND *TEMPERATURE FOR STRESS-DISPLACEMENT ELEMENTS  [C3D27]
*NODE
1    , 0.        , 0.        , 0.
3    , 7.        , 0.        , 0.
19   , 0.        , 7.        , 0.
21   , 7.        , 7.        , 0.
7    , 0.        , 0.        , -7.
9    , 7.        , 0.        , -7.
25   , 0.        , 7.        , -7.
27   , 7.        , 7.        , -7.
*NGEN,NSET=ALL
1    , 3
1    , 7    , 3
3    , 9    , 3
4,6
7    , 9
1    , 19   , 9
3    , 21   , 9
7    , 25   , 9
9    , 27   , 9
10,16,3
12,18,3
10,12
13,15
16,18
19   , 25   , 3
21   , 27   , 3
19   , 21
25   , 27
22,24
*NSET,NSET=NS1
1,4,7,16,25,22,19,10,13
*NSET,NSET=NS2
1,2,3,6,9,8,7,4,5
*NSET,NSET=NS3
1,2,3,12,21,20,19,10,11
*ELEMENT,TYPE=C3D27, ELSET=EALL
1,  1,3,9,7,19,21,27,25,2,6,8,4,20,24,26,
     22,10,12,18,16,14,5,23,11,15,17,13
*surface, NAME=SIDE_1
 1, S1
*surface, NAME=SIDE_2
 1, S2
*surface, NAME=SIDE_3
 1, S3
*surface, NAME=SIDE_4
 1, S4
*surface, NAME=SIDE_5
 1, S5
*surface, NAME=SIDE_6
 1, S6
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
 3.E6,0.3
*DENSITY
 10.0, 
*EXPANSION
 .0001, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL,-10.
*BOUNDARY
 ALL,1,3
*STEP,PERTURBATION
*STATIC
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*DSLOAD
 SIDE_1,P,3000.
*EL PRINT
 LOADS,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_2,P,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_3,P,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_4,P,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_5,P,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_6,P,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_1,HP,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_2,HP,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_3,HP,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_4,HP,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_5,HP,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_6,HP,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
 ALL,40.
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S, 
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*EL PRINT
 S, 
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
 NS1,1
 NS2,2
 NS3,3
*TEMPERATURE,OP=MOD
 ALL,40.
*OUTPUT,FIELD
*NODE OUTPUT
 U, 
*EL PRINT,FREQ=0
*NODE PRINT
 U, 
*END STEP