*parameter
neig = 50
** 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
285, 0, 40, 5
286, 0, 40, 0
287, 30, 40, 0
288, 30, 40, 5
*element,type=z101,file=box_ac3d4_moddyn,elset=z1
10001,273,274,275,276,285,286,287,288
*substructure property,elset=z1
*boundary
retained,1,3
*STEP
*FREQUENCY,norm=mass
<neig>,,,-1
*OUTPUT,FIELD
*NODE OUTPUT
U,
*END STEP
**------------------------------------------------
*STEP
*DYNAMIC,alpha=0,initial=no
0.001,0.1
*sload
10001,load401,1.0
*output,history
*node output,nset=retained
rf,
*SUBSTRUCTURE PATH, ENTER=10001
*OUTPUT,HISTORY
*NODE OUTPUT,nset=cload
u,
*SUBSTRUCTURE PATH, LEAVE
*end step