*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
  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
*element, type=s4r, elset=shell4
*elgen, elset=shell4
*elset, elset=shellUp
 shell3, shell4
*shell section, elset=shellUp, material=mat
 3.0,
*element,type=rotaryI,elset=rotaryI
**120051,51
**120032,32
**120031,31
**110041,10041
**110042,10042
**110052,10052
*rotary inertia,elset=rotaryI
<rota>,<rota>,<rota>,<rotaoff>,<rotaoff>,<rotaoff>
*boundary
 lipUp, 1,1
 lipUp, 3,3
*nset, nset=lipUp
**------------------------------------------------
*ncopy, old set=nallUp, new set=nallBot,change number=10000,shiftnallup
*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
**------------------------------------------------
*elset,elset=up
*surface, name=up
 up,spos
*elset,elset=down
*surface, name=down
 down,spos
*node,nset=swpset
*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
*surfaceinteraction, name=INTPROP-1
**------------------------------------------------
*amplitude,name=ramp,def=smooth
0,0,<time_step>, 1.0
**------------------------------------------------
*step,nlgeom=yes
*dynamic, explicit
, <time_step>
*bulkviscosity
*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
*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