A mathematical dynamic model of the two-dimensional representation of the knee joint is presented. The profiles of the joint surfaces are determined from X-ray films and they are represented by polynomials. The joint ligaments are modelled as nonlinear elastic springs of realistic stiffness properties. Nonlinear equations of motion coupled with nonlinear constraint conditions are solved numerically. Time derivatives are approximated by Newmark difference formulae and the resulting nonlinear algebraic equations are solved employing the Newton-Raphson iteration scheme. Several dynamic loads are applied to the center of mass of the tibia and the ensuing motion is investigated. Numerical results on ligament forces, contact point locations between femur and tibia, and the orientation of tibia relative to femur are presented. The results are shown to be consistent with the anatomy of the knee joint.