*HEADING
: DLOAD, DECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS [C3D8E]
*NODE,NSET=ALL
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=C3D8E, ELSET=EALL
 1,  1,2,3,4,5,6,7,8
**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=A1, ELSET=EALL,ORIENT=RECT
*MATERIAL,NAME=A1
*ELASTIC
 3.E6,0.3
*PIEZOELECTRIC,TYPE=S
 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,-10.
*BOUNDARY
 NS1,1,1
 NS1,3,3
 NS1,9,9
 NS2,2,2
 NS2,9,9
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES2,3000.
 1,ES4,3000.
*EL FILE
 S, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*NODE FILE
 U,EPOT,RF
*OUTPUT,FIELD
*NODE OUTPUT
U,EPOT,RF
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
 ALL,40.
*EL FILE
 S, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*NODE FILE
 U,EPOT,RF
*OUTPUT,FIELD
*NODE OUTPUT
U,EPOT,RF
*EL PRINT
 S, 
*NODE PRINT
 U,EPOT
*END STEP