In biomechanical investigations, geometrically accurate computer models of anatomical structures can be created readily using computed-tomography scan images. However, representation of soft tissue structures is more challenging, relying on approximations to predict the muscle loading conditions that are essential in detailed functional analyses. Here, using a sophisticated multi-body computer model of a reptile skull (the rhynchocephalian Sphenodon), we assess the accuracy of muscle force predictions by comparing predicted bite forces against in vivo data. The model predicts a bite force almost three times lower than that measured experimentally. Peak muscle force estimates are highly sensitive to fibre length, muscle stress, and pennation where the angle is large, and variation in these parameters can generate substantial differences in predicted bite forces. A review of theoretical bite predictions amongst lizards reveals that bite forces are consistently underestimated, possibly because of high levels of muscle pennation in these animals. To generate realistic bites during theoretical analyses in Sphenodon, lizards, and related groups we suggest that standard muscle force calculations should be multiplied by a factor of up to three. We show that bite forces increase and joint forces decrease as the bite point shifts posteriorly within the jaw, with the most posterior bite location generating a bite force almost double that of the most anterior bite. Unilateral and bilateral bites produced similar total bite forces; however, the pressure exerted by the teeth is double during unilateral biting as the tooth contact area is reduced by half. © 2010.