*preprint,model=yes
*heading
 Test of a circular plate
     
       ^                 ^
       |                 |
       o--o--o--X--o--o--o
 Arbitrary rotary inertia added to test anisotropic rotations
*preprint, echo=yes, model=yes
*parameter
**** Geometry of problem
 rota = 1.e-12
 rotaoff = rota*0.5
 time_step = 0.005
 time_interval = time_step*0.01
 ang45 = sqrt(2.0)/2.0
 d1 = 4.0
 d2 = 8.0
 d3 = 12.
 d4 = 16.
 d5 = 20.
 nd1 = -d1
 nd2 = -d2
 nd3 = -d3
 nd4 = -d4
 nd5 = -d5
 r1 = d1*ang45
 r2 = d2*ang45
 r3 = d3*ang45
 r4 = d4*ang45
 r5 = d5*ang45
 nr1 = -r1
 nr2 = -r2
 nr3 = -r3
 nr4 = -r4
 nr5 = -r5
**------------------------------------------------
*node, nset=nallUp
 1001, 0,0,0
  1, <d1>,0,0
  2, <d2>,0,0
  3, <d3>,0,0
  4, <d4>,0,0
  5, <d5>,0,0
 11, <r1>,0,<nr1>
 12, <r2>,0,<nr2>
 13, <r3>,0,<nr3>
 14, <r4>,0,<nr4>
 15, <r5>,0,<nr5>
 21, 0,0, <nd1>
 22, 0,0, <nd2>
 23, 0,0, <nd3>
 24, 0,0, <nd4>
 25, 0,0, <nd5>
 31, <nr1>,0,<nr1>
 32, <nr2>,0,<nr2>
 33, <nr3>,0,<nr3>
 34, <nr4>,0,<nr4>
 35, <nr5>,0,<nr5>
 41, <nd1>,0,0
 42, <nd2>,0,0
 43, <nd3>,0,0
 44, <nd4>,0,0
 45, <nd5>,0,0
 51, <nr1>,0,<r1>
 52, <nr2>,0,<r2>
 53, <nr3>,0,<r3>
 54, <nr4>,0,<r4>
 55, <nr5>,0,<r5>
 61, 0,0, <d1>
 62, 0,0, <d2>
 63, 0,0, <d3>
 64, 0,0, <d4>
 65, 0,0, <d5>
 71, <r1>,0,<r1>
 72, <r2>,0,<r2>
 73, <r3>,0,<r3>
 74, <r4>,0,<r4>
 75, <r5>,0,<r5>
*element, type=s3r, elset=shell3
 101, 1001,  1, 11
 111, 1001, 11, 21
 121, 1001, 21, 31
 131, 1001, 31, 41
 141, 1001, 41, 51
 151, 1001, 51, 61
 161, 1001, 61, 71
 171, 1001, 71,  1
*element, type=s4r, elset=shell4
  1,  1,  2, 12, 11
  2,  2,  3, 13, 12
  3,  3,  4, 14, 13
  4,  4,  5, 15, 14
 71, 71, 72,  2,  1
 72, 72, 73,  3,  2
 73, 73, 74,  4,  3
 74, 74, 75,  5,  4
*elgen, elset=shell4
 1, 7, 10, 10
 2, 7, 10, 10
 3, 7, 10, 10
 4, 7, 10, 10
*elset, elset=shellUp
 shell3, shell4
*shell section, elset=shellUp, material=mat
 3.0,
*element,type=rotaryI,elset=rotaryI
120041,41
120042,42
120052,52
**120051,51
**120032,32
**120031,31
**110041,10041
**110042,10042
**110052,10052
110051,10051
110032,10032
110031,10031
*rotary inertia,elset=rotaryI
<rota>,<rota>,<rota>,<rotaoff>,<rotaoff>,<rotaoff>
*boundary
 lipUp, 1,1
 lipUp, 3,3
*nset, nset=lipUp
 5, 15, 25, 35, 45, 55, 65, 75
**------------------------------------------------
*ncopy, old set=nallUp, new set=nallBot,change number=10000,shift
0,-3,0

*elcopy, old set =shellUp, new set=shellBot, element shift=10000, shift nodes=10000
*shell section, elset=shellBot, material=mat
 3.0,
*boundary
 lipBot, 1,1
 lipBot, 3,3
*nset, nset=lipBot
 10005, 10015, 10025, 10035, 10045, 10055, 10065, 10075
**------------------------------------------------
*elset,elset=up
31,41
*surface, name=up
 up,spos
*elset,elset=down
10031,10041
*surface, name=down
 down,spos
*node,nset=swpset
100355,-5.0, 0, 0
100356,-5.0, -3, 0
*MPC
BEAM,100356,100355
*COUPLING, CONSTRAINT NAME=up  , REF NODE=100355, SURFACE=up
*DISTRIBUTING,coupling=structural
*COUPLING, CONSTRAINT NAME=down, REF NODE=100356, SURFACE=down
*DISTRIBUTING,coupling=structural
**------------------------------------------------
*material, name=mat
*ELASTIC
  2.070000e+05,   3.000000e-01
*density
 7.850000e-09,
*PLASTIC
    292.0000,   0.0000e+00
    306.3399,   8.1205e-03
    325.2581,   2.0301e-02
    343.8255,   3.8573e-02
    364.3499,   6.5979e-02
    388.5583,   1.0709e-01
    416.8968,   1.6875e-01
    451.1431,   2.6125e-01
    494.2375,   4.0000e-01
*surfaceinteraction, name=INTPROP-1
**------------------------------------------------
*amplitude,name=ramp,def=smooth
0,0,<time_step>, 1.0
**------------------------------------------------
*step,nlgeom=yes
*dynamic, explicit
, <time_step>
*bulkviscosity
0.06, 1.2
*contact, op=NEW
*contactinclusions, allelementbased
*contactpropertyassignment
, , INTPROP-1
*boundary,amplitude=ramp
 lipUp,1,1,5.0
 lipBot,1,1,-5.0
*output,field,number interval=20
*node output
 u,
*contact output
Cforce,
**-----------------
*nset,nset=QA_TEST
100355,
*output, history, time interval=0.001e-2
*node output, nset=QA_TEST
u,
*node output,nset=LipUp
 rf1
*node output,nset=LipBot
 rf1
*energy output
*end step