*HEADING
EC34EFDM: *DSLOAD, *DSECURRENT, AND TEMPERATURE FOR PIEZO. ELEMENTS  [C3D4E]
*NODE,NSET=ALL_NODES_IN_MODEL
1    , 0.        , 0.        , 0.
2    , 7.        , 0.        , 0.
3    , 7.        , 0.        , -7.
4    , 0.        , 0.        , -7.
5    , 0.        , 7.        , 0.
6    , 7.        , 7.        , 0.
7    , 7.        , 7.        , -7.
8    , 0.        , 7.        , -7.
*NSET,NSET=NS1
1,4,5,8
*NSET,NSET=NS2
1,2,3,4
*NSET,NSET=NS3
1,2,5,6
*ELEMENT,TYPE=C3D4E,ELSET=ALL_ELEMENTS
 1, 4,3,8,1
 2, 1,3,6,2
 3, 3,7,8,6
 4, 1,6,8,5
 5, 1,6,3,8
*surface, NAME=SIDE_1
 1, S1
 3, S1
*surface, NAME=SIDE_2
 2, S2
 2, S4
*surface, NAME=SIDE_3
 1, S2
 2, S2
*surface, NAME=SIDE_4
 2, S3
 3, S2
*surface, NAME=SIDE_5
 3, S3
 4, S3
*surface, NAME=SIDE_6
 1, S4
 4, S4
**local orientation matching global system 
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 1, 0.0
*SOLID SECTION,MATERIAL=ELMAT, ELSET=ALL_ELEMENTS,ORIENT=RECT
*MATERIAL,NAME=ELMAT
*ELASTIC
 3.E6,0.3
*PIEZOELECTRIC
 15.08,-5.207,-5.207,0.,0.,0.,0.,0.
 0.,12.710,0.,0.,0.,0.,0.,0.,
 12.710,0.
*DIELECTRIC,TYPE=ANISO
 5.872E-9,0.,6.752E-9,0.,0.,6.752E-9
*EXPANSION
 .0001, 
*DENSITY
 10.0, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL_NODES_IN_MODEL,-10.
*BOUNDARY
 ALL_NODES_IN_MODEL,1,3
 ALL_NODES_IN_MODEL,9
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_1,P,3000.
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*EL PRINT
 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
*DSECHARGE
 SIDE_1,ES,3000.
*EL PRINT
 CHRGS,
*OUTPUT,FIELD
*NODE OUTPUT
RCHG, 
*NODE PRINT
RCHG, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSECHARGE
 SIDE_2,ES,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSECHARGE
 SIDE_3,ES,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSECHARGE
 SIDE_4,ES,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSECHARGE
 SIDE_5,ES,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSECHARGE
 SIDE_6,ES,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
 ALL_NODES_IN_MODEL,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
 ALL_NODES_IN_MODEL,9
*TEMPERATURE,OP=MOD
 ALL_NODES_IN_MODEL,40.
*OUTPUT,FIELD
*NODE OUTPUT
 U, 
*EL PRINT,FREQ=0
*NODE PRINT
 U, 
*END STEP