*heading
 BASIC VERIFICATION - MPC UNIVERSAL ON BEAM ELEMENTS
 FOR LINEAR DYNAMICS
*node
1,0,0
2,1,0
4,1,0
3,2,0
5,1,0,1
6,1,1,0
*element, elset=eall, type=b31
1,1,2
2,4,3
*beam section, elset=eall, section=rect, material=mat1
4.0,3.5
*MPC
PIN,2,4
UNIVERSAL,2,4,5,6
*material, name=mat1
*elastic
4.32e9, 
*density
1.0e-5, 
*boundary
1,1,6
3,1,6
*SPECTRUM,NAME=SPEC,TYPE=ACCELERATION
 0., 0., 0.
 1.E5, 1000., 0.
 1.E5, 10000., 0.
*restart,write,freq=999
********
*step
*FREQUENCY, eigensolver=lanczos
9, , 
*EL PRINT,FREQ=999
*EL FILE,FREQ=999
 S, 
 E, 
*NODE PRINT
 U, 
 RF, 
*NODE FILE
 U, 
 RF, 
*endstep
*STEP
*STEADY STATE DYNAMICS,DIRECT
0.01, 
0.1, 
1., 
*DLOAD
1,PY,1000.0
1,PZ,1000.0
*EL PRINT,FREQ=999
*EL FILE,FREQ=999
 S, 
 E, 
*NODE PRINT
 U, 
 RF, 
*NODE FILE
 U, 
 RF, 
*endstep
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
0.01, 
0.1, 
1., 
*DLOAD
1,PY,1000.0
1,PZ,1000.0
*EL PRINT,FREQ=999
*EL FILE,FREQ=999
 S, 
 E, 
*NODE PRINT
 U, 
 RF, 
*NODE FILE
 U, 
 RF, 
*endstep
*STEP
*MODAL DYNAMIC
.01,.03
*DLOAD
1,PY,1000.0
1,PZ,1000.0
*EL PRINT,FREQ=999
*EL FILE,FREQ=999
 S, 
 E, 
*NODE PRINT
 U, 
 RF, 
*NODE FILE
 U, 
 RF, 
*endstep
*STEP
*RESPONSE SPECTRUM,SUM=ABS,COMP=ALG
 SPEC,0.,0.,1.,386.4
*EL PRINT,FREQ=999
*EL FILE,FREQ=999
 S, 
 E, 
*NODE PRINT
 U, 
*NODE FILE
 U, 
 RF, 
*END STEP