*HEADING
: DLOAD, DECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [CPS8RE], MULTIPLE *SUBMODEL OPTIONS
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*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
*NSET,NSET=ALL,GENERATE
1    , 9
*ELEMENT,TYPE=CPS8RE, ELSET=EALL
1    , 1    , 3    , 9    , 7    , 2    , 6    , 8    , 4
*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.,
12.710,0.
*DIELECTRIC,TYPE=ISO
5.872E-9, 
*EXPANSION
0.0001, 
*DENSITY
5.E-5, 
*BOUNDARY
LEFT,1,1
LEFT,9,9
BOTTOM,2,2
BOTTOM,9,9
*STEP,PERT
*STATIC
*DECHARGE
1,ES2,2000.
1,ES3,3000.
*CECHARGE
9,,-150.0
*NODE PRINT
U,EPOT
*NODE FILE
U,EPOT,RCHG
*OUTPUT,FIELD
*NODE OUTPUT
U,EPOT,RCHG
*END STEP
*STEP,PERT
*FREQUENCY, eigensolver=lanczos
5, , 
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*END STEP
*STEP,PERT
*STEADY STATE DYNAMICS
50., 
*DLOAD
1,P2,200.
1,P3,300.
*NODE PRINT
U,EPOT,PU,PHPOT
*NODE FILE
U,EPOT,RCHG,PU,PHPOT
*OUTPUT,FIELD
*NODE OUTPUT
U,EPOT,RCHG
*END STEP