*HEADING
ECE3EFDM: *DSLOAD, *DSECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [CPE3E]
*NODE,NSET=ALL
1    , 0.        , 0.
3    , 7.        , 0.
7    , 0.        , 7.
9    , 7.        , 7.
*NSET,NSET=LEFT
1    , 7
*NSET,NSET=BOT
1    , 3
*ELEMENT,TYPE=CPE3E,ELSET=ELEMS
1    , 1    , 3    , 9
2    , 1    , 9    , 7
*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
ALL  , 1    , 2
ALL  , 9
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_1,HP,3000.,7.,0.
*ENERGY FILE
*NODE FILE
 RF, 
*EL PRINT
 LOADS,
*EL FILE
 LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_2,HP,3000.,7.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_3,HP,3000.,7.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_4,HP,3000.,7.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DSLOAD
SIDE_1,P,3000.
*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
*DSECHARGE
SIDE_1,ES,3000.
*EL PRINT
 CHRGS,
*EL FILE
 CHRGS,
*NODE FILE
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
*TEMPERATURE
ALL,40.
*EL FILE
 S, 
*NODE FILE
 RF,U
*EL PRINT
 S, 
*NODE PRINT
 RF,U
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
LEFT,1
BOT,2
ALL,9
*TEMPERATURE,OP=MOD
ALL,40.
*END STEP