*HEADING
EC3FEFS2: PRESSURE AND HYDROSTATIC PRESSURE FOR 3-D SOLIDS  [C3D15E]
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*RESTART,WRITE
*NODE,NSET=NALL
*ELEMENT,TYPE=C3D15E, ELSET=EALL
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*surface,NAME=SIDE
1,S2
2,S2
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
 30.0E6,0.3
*density
 200
*PIEZOELECTRIC,TYPE=S
*DIELECTRIC,TYPE=ISO
*BOUNDARY
 1,1,3
 2,2
 4,3
 5,1
 1,9
 2,9
 3,9
 4,9
*STEP,PERTURBATION
*STATIC
*DLOAD
 1, P1, 1000.0
 1, P2, 1000.0
 1, P3, 1000.0
 1, P5, 1000.0
 2, P1, 1000.0
 2, P2, 1000.0
 2, P3, 1000.0
 2, P4, 1000.0
*DECHARGE
 1, ES1, 1000.0
 1, ES2, 1000.0
 1, ES3, 1000.0
 1, ES5, 1000.0
 2, ES1, 1000.0
 2, ES2, 1000.0
 2, ES3, 1000.0
 2, ES4, 1000.0
*CLOAD
  2,1,  -166.6667
  3,1,   166.6667
  3,2,   166.6667
  3,3,   333.3333
  4,1,   166.6667
  4,2,  -166.6667
  5,2,  -166.6667
  5,3,  -333.3333
*CECHARGE
*EL PRINT
NFORC,
*SECTION PRINT,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL
SOAREA,SOF
SOCF,
SOE,
*SECTION FILE,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL
SOAREA,SOF
SOCF,
SOE,
*ENERGY PRINT
*NODE PRINT
 U,COORD,EPOT
 RCHG,RF
*EL PRINT
 S,COORD,
 E,
 LOADS,
 CHRGS,
 EFLX,EPG
 ENER,
 ELEN,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S,COORD,E,
 EFLX,EPG
 ENER,
 ELEN,
*NODE OUTPUT
 U,COORD,EPOT,RCHG,RF
*OUTPUT,HISTORY
*ENERGY OUTPUT,VAR=ALL
*END STEP
*step,pert
*steady state dynamics, direct
50.,
*DLOAD
 1, P1, 2000.0
 1, P2, 2000.0
 1, P3, 2000.0
 1, P5, 2000.0
 2, P1, 2000.0
 2, P2, 2000.0
 2, P3, 2000.0
 2, P4, 2000.0
*DECHARGE
 1, ES1, 2000.0
 1, ES2, 2000.0
 1, ES3, 2000.0
 1, ES5, 2000.0
 2, ES1, 2000.0
 2, ES2, 2000.0
 2, ES3, 2000.0
 2, ES4, 2000.0
*node print
u,
rf,
*el print
s,
e,
*output, field, op=new
*node output
u,
rf,
*element output
s,
e,
*end step