SOLUTION
  # Compute modal basis of 2-modes and use that to create
  # frequency response functions over the covered frequency
  # range
  case eig
    eigen
    nmodes = 20
  case test2
     modalfrf
END

FILE
   geometry_file = 'beam_frf.e'
END

LOADS
  # Applied force and location which the frequency response
  # functions are output with respect to
  nodeset 500
  # Defines force in the Z direction
  force = 0.0 0.0 1.0
  scale = 1
  function = 1
END

FUNCTION 1
  # Force defined in frequency domain.  Constant white
  # noise spectrum 
  type LINEAR
  name "white noise"
  data 0.0 1.0
  data 200. 1.0
END

DAMPING
  # Mass proportional damping coefficient
  alpha = 5
END

BLOCK 1
   material = 1 // rubber linear
END

BLOCK 90
  # Set of rigid bar elements
  rbar
END

BLOCK 91
  # Specified per node mass and inertia
  conmass
  mass = 1e-3
  Ixx = 1e-3
  Iyy = 1e-3
  Izz = 1e-3
END

MATERIAL 1 // linear
   isotropic
   density 0.0343
   E 218
   nu = .499
END

PARAMETERS
  # Change the unit system for mass-inputs to weight
  wtmass = 0.002588
END

OUTPUTS
   disp
   stress
END

FREQUENCY
   freq_min = .1
   freq_step = .1
   freq_max = 50
   acceleration
   disp
   nodeset 2
END

ECHO
   mass=block
END