Knowledge of the muscle, ligament, and joint forces is important when planning orthopedic surgeries. Since these quantities cannot be measured in vivo under normal circumstances, the best alternative is to estimate them using musculoskeletal models. These models typically assume idealized joints, which are sufficient for general investigations but insufficient if the joint in focus is far from an idealized joint. The purpose of this study was to provide the mathematical details of a novel musculoskeletal modeling approach, called force-dependent kinematics (FDK), capable of simultaneously computing muscle, ligament, and joint forces as well as internal joint displacements governed by contact surfaces and ligament structures. The method was implemented into the anybody modeling system and used to develop a subject-specific mandible model, which was compared to a point-on-plane (POP) model and validated against joint kinematics measured with a custom-built brace during unloaded emulated chewing, open and close, and protrusion movements. Generally, both joint models estimated the joint kinematics well with the POP model performing slightly better (root-mean-square-deviation (RMSD) of less than 0.75 mm for the POP model and 1.7 mm for the FDK model). However, substantial differences were observed when comparing the estimated joint forces (RMSD up to 24.7 N), demonstrating the dependency on the joint model. Although the presented mandible model still contains room for improvements, this study shows the capabilities of the FDK methodology for creating joint models that take the geometry and joint elasticity into account.