*parameter
neig = 75
**set neig=150 to obtain the results in the manuals
*amplitude,name=const
0,1.0,  1.0,1.0
*Node,nset=retained
273, 30, 0, 5
274, 0, 0, 5
275, 0, 0, 0
276, 30, 0, 0
277, 25, 5, 5
285, 0, 40, 5
286, 0, 40, 0
287, 30, 40, 0
288, 30, 40, 5
*element,type=z101,file=box_ac3d4,elset=z1
10001,273,274,275,276,285,286,287,288
*substructure property,elset=z1
*STEP
*FREQUENCY,norm=mass
<neig>,,,-1
*boundary,op=new
285,1,3
286,1,3
287,1,3
288,1,3
275,1,3
276,1,3
273,2,3
274,2,3
*boundary,op=new,base name=excit1
273,1,1
274,1,1
*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
1.1, 1.9, 50
*base motion,dof=1,loadcase=1,amplitude=const,base name=excit1,type=disp
*output, field, frequency=0
*output, history
*nodeoutput, nset=retained
U,
*endstep
*step
*steady state dynamics,direct,interval=range
1.1, 1.9, 50
*boundary,op=new
285,1,3
286,1,3
287,1,3
288,1,3
275,1,3
276,1,3
273,2,3
274,2,3
*boundary, op=new, amplitude=const,loadcase=1
273, 1, 1, 1.0
274, 1, 1, 1.0
*output, field, frequency=0
*output, history
*nodeoutput, nset=retained
U,
*end step