*HEADING
Test of angular acceleration loading for C3D20E elements.
*RESTART,WRITE
*NODE,NSET=NALL
1,9.,-1.,-1.
2,11.,-1.,-1.
3,11,1.,-1.
4,9.,1.,-1.
5,9.,-1.,1.
6,11.,-1.,1.
7,11.,1.,1.
8,9.,1.,1.
9,10.,-1.,-1.
10,11.,0.,-1.
11,10.,1.,-1.
12,9.,0.,-1.
13,10.,-1.,1.
14,11.,0.,1.
15,10.,1.,1.
16,9.,0.,1.
17,9.,-1.,0.
18,11.,-1.,0.
19,11.,1.,0.
20,9.,1.,0.
*ELEMENT,TYPE=C3D20E,ELSET=EALL
1,1,2,3,4,5,6,7,8,9,10,
11,12,13,14,15,16,17,18,19,20
**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=STEEL,ELSET=EALL,ORIENT=RECT
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,.33
*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
*DENSITY
2000., 
*BOUNDARY
NALL,1,3
*STEP,PERTURBATION
*STATIC
*DLOAD
EALL,ROTA,10.,0.,0.,0.0,0,1
*NODE PRINT
RF, 
*EL PRINT,FREQ=0
*OUTPUT,FIELD
*NODE OUTPUT
RF, 
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT
*END STEP