Exact Computation of Contact Regions for Nonlinear Beams with Cubic Interpolation