begin sierra beam_sm_fft
  # This is a Sierra solid mechanics input deck that computes time domain
  # reponse to the input loads.  I can be used to compare to the Sierra/SD
  # frequency domain results after FFT transformation

  begin function prescribed_force
    type is piecewise analytic
    begin expressions
        0.0 "1e-4*sin(2*pi*t)"
    end expressions
  end

  begin material rubber
    density = {0.0343*0.002588}

    begin parameters for model elastic
      poissons ratio = 0.499
      youngs modulus = 218
    end parameters for model elastic
  end material rubber

  begin material rbar
    density = 0
    begin parameters for model elastic
      poissons ratio = 0
      youngs modulus = 1e-7
    end parameters for model elastic
  end material rbar  

  begin rigid body rbar
  end rigid body rbar 

  begin beam section rbar_sec
  rigid body = rbar
  section = bar
  width = 1e-7
  height = 1e-7
  t axis = 0 0 1
  end

  begin point mass section conmass
  mass = {1e-3*0.002588}
  end

  begin finite element model fft_run
    database name = beam_frf.e
    database type = exodusII
    
    # - Block id 1 had name bar
    begin parameters for block block_1
      material = rubber
      model = elastic
    end parameters for block block_1
        
    # - Block id 90 had name rbar
    begin parameters for block block_90
      material =rbar
      model = elastic
      section = rbar_sec
    end parameters for block block_90
    
    # - Block id 91 had name conmass
    begin parameters for block block_91
      section = conmass
    end parameters for block block_91

  end finite element model fft_run

  begin presto procedure beam_fft

    #
    # *** Time step control information
    begin time control

      begin time stepping block p1
        start time = 0.0
        begin parameters for presto region presto
        time step scale factor =  1.0
          step interval = 100
        end parameters for presto region
      end time stepping block p1

      termination time = 100

    end time control

    begin presto region presto
      
      begin viscous damping
         include all blocks
         mass damping coefficient = 5
      end viscous damping

      use finite element model fft_run
      ###  output description ###
      begin results output results
        start time = 0
        database name = beam_frf-out.e
        database type = exodusII
	At Time 0.0, Increment = 1.0e-1
	#At Time 0.0, Increment = 1.0e-5
        nodal Variables = displacement as displ
        nodal Variables = velocity as vel
        nodal Variables = acceleration as accel
      end results output results
      
      begin prescribed force
      node set = nodelist_500
      component = z
      function = prescribed_force
      scale factor = 1
      end prescribed force

    end presto region presto

  end presto procedure beam_fft

end sierra beam_sm_fft