*parameter
 neig=60
*Node,nset=retained
84, 50, 30, 0
87, 10, 30, 0
88, 10, -30, 0
89, 50, -30, 0
**87,       10.,             30. ,             0.
**84,   29.3649,             65. ,             0.
**89,   81.8649,         35.9526 ,             0.
**88,      62.5,    952.625E-03  ,            0.
*element,type=z101,file=box_ac2d4,elset=z1
10001,88,89,87,84
*substructure property,elset=z1
**0,0,0
**10, 30, 0,10, 30, 1,61.044975628
**------------------------------------------------
*amplitude,name=const
0,1.0,  1.0,1.0
*amplitude,name=constIm
0,1.0,  1.0,1.0
*STEP
*FREQUENCY,norm=mass
<neig>,,,-1
*boundary,op=new
87,1,2
84,1,2
*boundary,base name=excit1
88,1,1
89,1,1
*boundary,base name=excit2
88,2,2
89,2,2
*OUTPUT,FIELD
*NODE OUTPUT
U,
***SUBSTRUCTURE PATH, ENTER=10001
***OUTPUT,FIELD,freq=1
***NODE OUTPUT
**U,por
***SUBSTRUCTURE PATH, LEAVE
*END STEP
**------------------------------------------------
*Step, perturbation
*steadystatedynamics, subspaceprojection, interval=RANGE
0.2, 0.6, 50
*base motion,dof=1,loadcase=1,amplitude=const,base name=excit1,type=disp
*base motion,dof=2,loadcase=2,amplitude=constIm,base name=excit2,type=disp
*output, field, frequency=0
*output, history
*nodeoutput, nset=retained
U,
*endstep
*step
*steady state dynamics,direct,interval=range
0.2, 0.6, 50
*boundary,op=new
87,1,2
84,1,2
*boundary, op=new, amplitude=const,loadcase=1
88, 1, 1, 1.0
89, 1, 1, 1.0
*boundary, op=new, amplitude=constIm,loadcase=2
88, 2, 2, 1.0
89, 2, 2, 1.0
*output, field, frequency=0
*output, history
*nodeoutput, nset=retained
U,
*end step