SOLUTION
  title 'test fixture model transient simulation'
  transient
    # Transient solution, nsteps is number of time steps
    time_step 1.0e-4
    nsteps    100
End

PARAMETERS
  # Weight to mass conversion factor, inverse of gravity accel
  wtmass=0.00259
End

FILE
  geometry_file 'fixture.exo'
End

LOADS
  # force on each and every node of nodeset 1
  nodeset 1
    # Total force is in Y direction, and multiples by both scale
    # and the value of function 8 at the current time
    force 0.0 1.0 0.0
    scale 1.0e7
    function 8
End

HISTORY
  nodeset 33, 148, 270
  displacement
  acceleration
End

OUTPUTS
End

Function 8 // Haversine pulse
  type analytic
  # The first three lines define variables, the value computed
  # by the final line is the actual value returned by the function
  evaluate expression = "amp = 1.5e3; 
                         period= 3.6e-4; 
                         omega = pi/period; 
                         (t>period)?(0.0):(amp*sin(omega*t)^2)"
End

ECHO
  # block-by-block mass output for debugging
  mass block
End

Block 1 // fixture
  material 1
End

Block 2
  # Rigid bar
  rbar
End

Block 3
  # Concentrated mass and inertia element 
  ConMass
    Mass 1.0e7
    Ixx  1.0e8
    Iyy  1.0e8
    Izz  1.0e8
    Offset= 0.0 0.0 0.0
End

MATERIAL 1 // fixture - Ti
  density=0.16
  E=1.6e+07     # Young's Modulus
  nu=0.3        # Poisson's ratio
End

GDSW
  # Tighter than default solver tolerance
  solver_tol 1.0e-8
End