This paper is concerned with a mathematical model analysis of the patellofemoral joint in the human knee, taking into account the articular surface geometry and mechanical properties of the ligament. It was made by the application of a computer-aided design theory (previously studied) and it was possible to express the articular surface geometries in a mathematical formulation and hence elucidate the joint movement mechanics.

This method was then applied to a three-dimensional geometrical model of the patellofemoral joint. For the modelling of tendofemoral contact at large angles of knee flexion, the geodestic line theory was adopted. Applying the Newton-Raphson method and the Runge-Kutta Gil method to the model, variables such as patellar attitudes, patellofemoral contact force and tensile force of the patellar ligament for various knee flexion angles were computed. Applying the Hertzian elastic theory, contact stress was also computed. These results showed good agreement with the previously reported experimental results.

As an application for the model, some parameter analyses were performed in terms of the contact stress variations and compared with those of the normal knee. The simulation results indicated that both the Q-angle increase and decrease increased contact stress, the patella alta showed undulating variations of stress while the patella infera showed little change of stress, and the tibial tuberositas elevation showed 20–30% reduction of stress.