*heading
Coriolis forces for C3D20
Analytical results:   Disp  RF = 39.97  PRF = 170.96
                  : Coriol  RF = 12.57  PRF = +/-90
*node,nset=nall
1,0,0,0
2,1,0,0
3,1,1,0
4,0,1,0
5,0,0,1
6,1,0,1
7,1,1,1
8,0,1,1
9,0.5,0,0
10,1,0.5,0
11,0.5,1,0
12,0,0.5,0
13,0.5,0,1
14,1,0.5,1
15,0.5,1,1
16,0,0.5,1
17,0,0,0.5
18,1,0,0.5
19,1,1,0.5
20,0,1,0.5
*element,elset=one,type=c3d20
1,1,2,3,4,5,6,7,8,
9,10,11,12,13,14,15,16,
17,18,19,20
*solid section,elset=one,material=elas
*material,name=elas
*elastic
1000,.3
*density
1,
*damping,alpha=1.
*step,perturbation,unsymm=yes
*steady state dynamics,direct
1, 
*load case,name=case-1
*boundary,op=new
nall,1,1,1
nall,2,2
nall,3,3
*end load case
*load case,name=case-2
*dload,op=new
one,corio,1,.5,-10,0,0,0,1
*boundary,op=new
nall,1,1,1
nall,2,2
nall,3,3
*end load case
*load case,name=case-3
*boundary,op=new
nall,1,1
nall,2,2,1
nall,3,3
*end load case
*load case,name=case-4
*dload,op=new
one,corio,1,.5,-10,0,0,0,1
*boundary,op=new
nall,1,1
nall,2,2,1
nall,3,3
*end load case
*node print
u,pu
v, 
*node print,tot=yes
rf,prf
*el print,f=0
*output,field
rf,prf
*end step