*HEADING
: SUBMODEL FROM PGC38EFM  [C3D8E]
*NODE,NSET=ALL
1,    3.5,      0.,        0.
2,    7.,       0.,        0.
3,    7.,       0.,       -7.
4,    3.5,      0.,       -7.
5,    3.5,      7.,        0.
6,    7.,       7.,        0.
7,    7.,       7.,       -7.
8,    3.5,      7.,       -7.
9,    3.5,      3.5,       0.
10,   7.,       3.5,       0.
11,   7.,       3.5,      -7.
12,   3.5,      3.5,      -7.
*NSET,NSET=NS1
1,4,5,8,9,12
*NSET,NSET=NS2
1,2,3,4
*NSET,NSET=NS3
1,2,5,6
*ELEMENT,TYPE=C3D8E, ELSET=EALL
 1,  1,2,3,4,9,10,11,12
 2,  9,10,11,12,5,6,7,8
*SUBMODEL
NS1, 
**local orientation matching global system 
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 1, 0.0
*SOLID SECTION,MATERIAL=A1, ELSET=EALL,ORIENT=RECT
*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=ANISO
 5.872E-9,0.,6.752E-9,0.,0.,6.752E-9
*EXPANSION
 .0001, 
*DENSITY
 10.0, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL,-10.
*BOUNDARY
 2,2,2
 2,9,9
 3,2,2
 3,9,9
*STEP,PERTURBATION
*STATIC
*BOUNDARY,SUBMODEL,STEP=1
NS1,1,3
NS1,9,9
*DECHARGE
 2,ES2,3000.
 1,ES4,3000.
 2,ES4,3000.
*EL FILE
 S, 
*NODE FILE
 U,EPOT,RF
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,SUBMODEL,STEP=2
NS1,1,3
NS1,9,9
*TEMPERATURE
 ALL,40.
*EL FILE
 S, 
*NODE FILE
 U,EPOT,RF
*EL PRINT
 S, 
*NODE PRINT
 U,EPOT
*END STEP