*HEADING
EC3FEFS2: PRESSURE AND HYDROSTATIC PRESSURE FOR 3-D SOLIDS  [C3D15E]
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*RESTART,WRITE
*NODE,NSET=NALL
  1, 0.0, 0.0, 0.0
  2, 2.0, 0.0, 0.0
  3, 2.0, 2.0, 0.0
  4, 0.0, 2.0, 0.0
  5, 0.0, 0.0, 1.0
  6, 2.0, 0.0, 1.0
  7, 2.0, 2.0, 1.0
  8, 0.0, 2.0, 1.0
  9, 1.0, 0.0, 0.0
 10, 2.0, 1.0, 0.0
 11, 1.0, 2.0, 0.0
 12, 0.0, 1.0, 0.0
 13, 1.0, 0.0, 1.0
 14, 2.0, 1.0, 1.0
 15, 1.0, 2.0, 1.0
 16, 0.0, 1.0, 1.0
 17, 0.0, 0.0, 0.5
 18, 2.0, 0.0, 0.5
 19, 2.0, 2.0, 0.5
 20, 0.0, 2.0, 0.5
 21, 1.0, 1.0, 0.0
 22, 1.0, 1.0, 1.0
*ELEMENT,TYPE=C3D15E, ELSET=EALL
 1, 1,2,4,5,6,8,9,21,12,13,22,16,17,18,20
 2, 2,3,4,6,7,8,10,11,21,14,15,22,18,19,20
*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
0.,0.,0.,0.,0.,0.,0.,0.
0.,0.,0.,0.,0.,0.,0.,0.
0.,0.
*DIELECTRIC,TYPE=ISO
1.0E-3, 
*BOUNDARY
 1,1,3
 2,2
 4,3
 5,1
 1,9
 2,9
 3,9
 4,9
 9,9
10,9
11,9
12,9
21,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
  6,1,  -166.6667
  6,2,   166.6667
  7,1,   166.6667
  7,2,   166.6667
  7,3,   333.3333
  8,1,   166.6667
  8,2,  -166.6667
  9,1,   1333.3333
  9,2,   666.6667
  9,3,   666.6667
 10,1,   666.6667
 10,3,  -666.6667
 11,2,   666.6667
 11,3,  -666.6667
 12,1,   666.6667
 12,2,   1333.3333
 12,3,   666.6667
 13,2,  -666.6667
 13,3,   666.6667
 14,1,  -666.6667
 14,2,  -1333.3333
 14,3,  -666.6667
 15,1,  -1333.3333
 15,2,  -666.6667
 15,3,  -666.6667
 16,1,  -666.6667
 16,3,   666.6667
 17,1,   666.6667
 17,2,   666.6667
 17,3,   1333.3333
 18,1,   666.6667
 18,2,  -666.6667
 19,1,  -666.6667
 19,2,  -666.6667
 19,3,  -1333.3333
 20,1,  -666.6667
 20,2,   666.6667
 21,1,   1333.3333
 21,2,   1333.3333
 22,1,  -1333.3333
 22,2,  -1333.3333
*CECHARGE
1,,333.333
2,,333.333
3,,333.333
4,,333.333
5,,333.333
6,,333.333
7,,333.333
8,,333.333
9,,-1333.333
10,,-1333.333
11,,-1333.333
12,,-1333.333
13,,-666.667
14,,-666.667
15,,-666.667
16,,-666.667
17,,-1333.333
18,,-1333.333
19,,-1333.333
20,,-1333.333
21,,-1333.333
*EL PRINT
NFORC,
*SECTION PRINT,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL
 ,1.,1,1.
 ,0.,1.,1., ,1.,0.,1
SOAREA,SOF
SOCF,
SOE, 
*SECTION FILE,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL
 ,1.,1,1.
 ,0.,1.,1., ,1.,0.,1.
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