SOLUTION
  # Model tests behavior under quasistatic loads associated with a constant sustained
  # Angluar velocity and angular acceleration
  statics
END

parameters
end

FILE
  # This defines the geometry to use in the analysis
  geometry_file hex20Beam40x.g
END

BOUNDARY
  sideset baseSurf
    fixed
END

LOADS
  # Angular acceleration of 200 radians/s^2 and angular velocity 500 rad/s.
  # These are defined in the with respect to the 'offset_Y' coordinate frame
  # defined below.  the 'offset_y' has directions the same as the global XYZ system
  # but with the origin point moved to '0 1 0'.  Thus these boundary conditions
  # define a spin around the a Z axis defined between the coordinates of
  # '0 1 0' and '0 1 1'
  #
  #  Angular acceleration creates an inertial force that is tangential around this
  #  axis while angular velocity creates a centripetal force around this axis.
  #  Even though angular acceleration implies the velocity velocity would be changing
  #  as a function of time this static analysis assumes both these forces are constant.
  #
  body
    angular_acceleration = 0 0 200
    coordinate offset_Y
  body
    angular_velocity = 0 0 500 
    coordinate offset_Y
END

OUTPUTS
  # Only output is deformed shape
  disp
END

ECHO
END	

BLOCK myTestBeam
  material example_mat
END

Material example_mat
  E 30.0E6         # Young's modulus (psi)
  nu 0.33          # Poisson ratio   (dimensionless)
  density 0.00074  # 'slinch'/in^3
END

begin rectangular coordinate system offset_Y
  # same coordinate directions as global XYZ, but with the origin at 0 1 0
  origin = 0 1 0
  z point = 0 1 5
  xz point = 1 1 0
end