*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=s4, 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=shellUp1, generate
   1, 71, 10
   2, 72, 10
   3, 73, 10
 101,171, 10
*elset, elset=shellBot1, generate
 10001, 10071, 10
 10002, 10072, 10
 10003, 10073, 10
 10101, 10171, 10
*elset, elset=shellUp
 shell3, shell4
*shell section, elset=shellUp, material=mat
 3.0,
*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,
*nset, nset=lipBot
 10005, 10015, 10025, 10035, 10045, 10055, 10065, 10075
*nset, nset=fix2
 10005,10035,10055
*boundary
 fix2, 2,2
 lipUp, 3,3
 lipBot, 3,3
**------------------------------------------------
*surface, name=up
 shellUp1,spos
*surface, name=down
 shellBot1,spos
*elset,elset=up1
31,41
*surface, name=up1
 up1,spos
*elset,elset=down1
10031,10041
*surface, name=down1
 down1,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=up1
*DISTRIBUTING,coupling=structural
*COUPLING, CONSTRAINT NAME=down, REF NODE=100356, SURFACE=down1
*DISTRIBUTING,coupling=structural
**------------------------------------------------
*material, name=mat
*ELASTIC
  2.070000e+05,   3.000000e-01
*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
*contact pair, interaction=smooth, type=surface to surface
 up, down
*surfaceinteraction, name=smooth
**------------------------------------------------
*step,nlgeom=yes
*static
 0.005, 1.
*boundary
 lipUp,1,1,5.0
 lipBot,1,1,-5.0
*output,field, freq=2
*node output
 u,
*contact output
 cstress,
**-----------------
*nset,nset=QA_TEST
100355,
*output, history
*node output, nset=QA_TEST
u,
*node output,nset=LipUp
 rf1
*node output,nset=LipBot
 rf1
*energy output
*end step