Obtaining tibio-femoral (TF) contact forces, ligament deformations and loads during daily life motor tasks would be useful to better understand the aetiopathogenesis of knee joint diseases or the effects of ligament reconstruction and knee arthroplasty. However, methods to obtain this information are either too simplified or too computationally demanding to be used for clinical application. A multibody dynamic model of the lower limb reproducing knee joint contact surfaces and ligaments was developed on the basis of magnetic resonance imaging. Several clinically relevant conditions were simulated, including resistance to hyperextension, varus-valgus stability, anterior -posterior drawer, loaded squat movement. Quadriceps force, ligament deformations and loads, and TF contact forces were computed. During anterior drawer test the anterior cruciate ligament (ACL) was maximally loaded when the knee was extended (392 N) while the posterior cruciate ligament (PCL) was much more stressed during posterior drawer when the knee was flexed (319 N). The simulated loaded squat revealed that the anterior fibres of ACL become inactive after 608 of flexion in conjunction with PCL anterior bundle activation, while most components of the collateral ligaments exhibit limited length changes. Maximum quadriceps and TF forces achieved 3.2 and 4.2 body weight, respectively. The possibility to easily manage model parameters and the low computational cost of each simulation represent key points of the present project. The obtained results are consistent with in vivo measurements, suggesting that the model can be used to simulate complex and clinically relevant exercises.