*HEADING
: *DSLOAD, *DSECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [CPS6E]
*NODE
1    , 0.        , 0.
3    , 7.        , 0.
7    , 0.        , 7.
9    , 7.        , 7.
*NGEN,NSET=BOTTOM
1    , 3
*NGEN,NSET=TOP
7    , 9
*NGEN,NSET=LEFT
1    , 7    , 3
*NGEN,NSET=RIGHT
3    , 9    , 3
*NGEN
4    , 6
*NSET,NSET=ALL,GENERATE
1    , 9
*NSET,NSET=SIDES
 1,2,3,4,5,6,7,8,9
*ELEMENT,TYPE=CPS6E,ELSET=ELEMS
1    , 1    , 3    , 9    , 2    , 6    , 5
2    , 1    , 9    , 7    , 5    , 8    , 4
*surface, NAME=SIDE_1
 1, S1
*surface, NAME=SIDE_2
 1, S2
*surface, NAME=SIDE_3
 2, S2
*surface, NAME=SIDE_4
 2, S3
*SOLID SECTION,MATERIAL=A1, ELSET=ELEMS
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
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=ISO
5.872E-9, 
*EXPANSION
0.0001, 
*DENSITY
5.E-5, 
*INITIAL CONDITION,TYPE=TEMPERATURE
ALL  , -10.
*BOUNDARY
SIDES,1,2
SIDES,9
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_1,HP,3000.,3.,0.
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*NODE FILE
 RF, 
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*EL PRINT
 LOADS,
*EL FILE
 LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_2,HP,3000.,3.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_3,HP,3000.,3.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_4,HP,3000.,3.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_1,P,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_2,P,2000.
*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
*DSECHARGE
SIDE_1,ES,1000.
*EL PRINT
 CHRGS,
*EL FILE
 CHRGS,
*NODE FILE
RCHG, 
*OUTPUT,FIELD
*NODE OUTPUT
RCHG,
*NODE PRINT
RCHG, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSECHARGE
SIDE_2,ES,2000.
*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
*TEMPERATURE
ALL,40.
*EL FILE
 S, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*NODE FILE
 RF,U
*OUTPUT,FIELD
*NODE OUTPUT
RF,U
*EL PRINT
 S, 
*NODE PRINT
 RF,U
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
LEFT,1
BOTTOM,2
ALL,9
*TEMPERATURE,OP=MOD
ALL,40.
*END STEP