SOLUTION
   # Frequency response function calculation
   case test2
      directfrf
END

FILE
   geometry_file = 'beam_frf.e'
END

GDSW
  # Option which effects GPU solution
  tacho_alg_variant 2
END

LOADS
  # Define frequency domain load function that drives the structure.
  # Same force applied to every node of nodeset 500
  nodeset 500
   # force in Z direction
   force = 0.0 0.0 1.0
   scale = 1
   function = 1
END

FUNCTION 1
   # This defines a spectrum of constant force of the frequency range
   type LINEAR
   name "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
  # All elements in block 90 are rigid bars.  Any two connected
  # rigid bars meld into a single rigid body.
  rbar
END

BLOCK 91
  # Concentrated mass element with specified mass and rotational
  # inertia
  conmass
  mass = 1e-3
  Ixx = 1e-3
  Iyy = 1e-3
  Izz = 1e-3
END

MATERIAL 1
   isotropic
   density 0.0343
   E 218             # Young's Modulus
   nu = .499         # Poisson's ratio.  A poisson ratio
                     # would be incompressible.  Specified Poisson
                     # ratio must be at least slightly below 0.5
END

PARAMETERS
   # Conversion factor between weight and mass.
   # With this parameter input deck terms like density
   # should be specified in terms of weight rather than
   # mass
   wtmass = 0.002588
END

OUTPUTS
   disp
   stress
END

FREQUENCY
  # Frequency domain outputs are calculated on
  # nodeset 2, this includes the acceleration
  # and displacement magnitudes due to the
  # input forces
  freq_min = .1
  freq_step = .1
  freq_max = 50
  acceleration
  disp
  nodeset 2
END

ECHO
  # Output block-by-block masses for debugging purposes
  mass=block
END