*HEADING
ECA4EFDM: *DSLOAD, *DSECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [CAX4E]
*NODE,NSET=ALL
 1,1.,0.
 3,4.,0.
 7,1.,3.
 9,4.,3.
*ELEMENT,TYPE=CAX4E, ELSET=EALL
 1,1,3,9,7
*surface, NAME=SIDE_1
 1, S1
*surface, NAME=SIDE_2
 1, S2
*surface, NAME=SIDE_3
 1, S3
*surface, NAME=SIDE_4
 1, S4
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*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.,
0.,0.
*DIELECTRIC,TYPE=ISO
5.872E-9, 
*DENSITY
5.E-5, 
*BOUNDARY
 ALL,1,2
 ALL,9
*STEP,PERTURBATION
*STATIC
*DSLOAD
 SIDE_1,HP,3000.,3.,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.,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, 
*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