Human neck musculature plays an important structural role during tension and bending modes of neck loading. Computational models of the head and neck require that the muscles follow the curvature of the neck during bending to maintain anatomically correct lines of action. The method used to model muscle wrapping affects head kinematics, neck loads, model stability and computational runtimes. This study presents a review of existing methods to model muscle wrapping. One additional method that takes advantage of a new contact interaction available within LS-DYNA (LSTC, Livermore, CA) is also presented. A comparison is made between methods and their effects on head kinematics, neck loads, model stability and computational runtimes during simulated frontal impact and airbag loading. This comparison shows that the new contact interaction within LS-DYNA has advantages over the other methods.