SOLUTION
  # Adds the 6 canonical rigid body modes
  case out
     geometric_rigid_body_modes
  # Add the first 10 flexible modes building on from the
  # 6 geometry rigid body modes
  case flexible_modes
     eigen 
     nmodes 10
END

FILE
   geometry_file 'simpleTied.exo'
END

BOUNDARY
  # free-free, no boundary conditions allowed when using
  # geometry rigid body modes
END

PARAMETERS
  # Number of expected rigid body modes, typically 6,
  # three translation and three rotation.
  num_rigid_mode 6
  # Allowed tolerance of rigid body modes, higher number allow
  # less perfect modes, such as may be caused by constraints
  # that don't perfectly allow rigid body rotation.
  RbmTolerance 2.e-6
  # Conversion between 'weight' and 'mass'.  Allows specifying
  # input deck parameters in pounds instead of slinches
  wtmass=0.00259
END

OUTPUTS
   disp
END

ECHO
  # Write block-by-block mass to rslt file
  mass block
END

LOADS
  # Apply uniform traction to all faces in sideset 3 
  sideset 3
    traction 1 1 1
    scale = 1.0
END

DAMPING
  # Uniform damping per node
  beta 2.0e-6
END

TIED JOINT
  # Adds node on face constraints that allow tangential
  # motion at the joint, but no normal motion.
  normal definition = slip
  # Joint is defined between surfaces 1 and 2
  surface 1,2
  # Slip constraints are defined using this geometric search
  # tolerance
  search tolerance 1.0e-3
  # Sidesets are not rigid
  side = free
  # Spring-like connection between the two sides of the
  # joint
  connect to block 3
END

BLOCK 1
   material 1
   nonlinear=no
END

BLOCK 2
   material 1
   nonlinear=no
END

BLOCK 3
 # the tied joint X, Y, Z response is with respect to coordinate system 2
 coordinate 2
 joint2g
 # The X and Y behavior of the joint are defined by PROPERTY 1' below
 kx = Iwan 1
 ky = Iwan 1
 # The rotational Z behavior is an elastic torsional spring
 krz = elastic 1.0e9
END

MATERIAL 1
   density 0.3
   E = 3.0e7    # Young's Modulus 
   nu = 0.3     # Poisson's Ratio
END

PROPERTY 1
    # Slope of force dissipation curve
    chi -.82
    # Breakage force
    phi_max = 1.75e-4
    #Constant coefficient
    R = 5.5050e+6
    # Strength of singularity in break free condition
    S = 2.1097e+6
END

Begin rectangular coordinate system 2
    origin   = 0 -3.83232e-2 -5.96407
    z point  = 1.0 -3.83232e-2 -5.96407
    xz point = 1.0 0.4616768 -6.46407
end

GDSW
  # Allows handling larger constraints as 'type 1' which is
  # usually better for the solver
  max_numterm_C1 500
  # Print constraint matrix for debugging
  prt_constraint 1
END