*Heading
*Node
      1,          -1.,           1.
      5,          -1.,          -1.
     21,           1.,           1.
     25,           1.,          -1.
*Element,type=z1,elset=one,file=lib
10001,1,5,21,25
*substructure property,elset=one
*node,nset=nref
201, 0.0, 2.0, 0.0
*element,type=mass,elset=mass
200001,201
*mass,elset=mass
40,
*mpc
beam, 1, 201
beam,21, 201
*initial conditions, type=velocity
201,1,0.1
201,2,0.05
** 
*Step, name=Stp1, inc=1000,nlgeom
Rigid Body Rotation
*Dynamic,alpha=-0.05,nohaf
0.1, 1.
*controls,param=field
1.e-8,
*Restart, write, frequency=1
*output, history
*node output,nset=nref
u,v
*substructure path,enter=10001
*output,field
*node output
u,
*element output
s,
*substructure path,leave
*End Step