Three-dimensional mathematical model analysis of the patellofemoral joint

Age related changes in the Q angle of non-professional football players

Background

Football practice involves a great muscular demand, leading to the development of the lower limbs that, on occasions, can cause deviations from the normal anthropometric values. The quadriceps angle (Q angle) is a value often taken as a reference for the alignment of the lower limbs.

Objective

To observe the changes of the Q angle in young football players, because of muscular effort, analyzing the differences between four groups of different ages and to determine whether the playing position might influence these variations.

Methods

A cross sectional study was carried out with 104 male subjects divided into four groups according to age: under 8 years-old, 8-17 years-old, 17-21 years-old and over 21 years-old. A photograph was taken in standing position and the Q angle was plotted with KINOVEA® software. As for the reliability of the measurements, intraclass intra and interobserver coefficient were 0.958 and 0.860 respectively. The study was conducted in mid-season.

Results

Q angle value is greater in those under 8 years of age and decreases gradually and significantly (p < 0.005) until 17-21 years of age, where it stabilizes at values of 5.73° ± 2.78 for right Q angle and 5.88° ± 2.55 for left Q angle. Two way ANOVA demonstrated a significant group*position interaction for goalkeepers with a medium effect size in both angles (p < 0.001) with a medium effect (η² Right Q angle = 0.31; η² Left Q angle = 0.37). The values remain unchanged in subjects over 21 years of age (p > 0.005), except for goalkeepers, who suffered a difference in the evolution of the angle within their age category (p < 0.005) and with a high effect size with the other positions (value > 0.8) except forward (value < 0.5).

Conclusion

This study determines that the Q angle in football players decreases with growth, reaching values below 15° at the end of development. Playing positions only influence players over the age of 21, and the Q-angle of goalkeepers is greater than that of other players.

Computational biomechanics of human knee joint in stair ascent: Muscle-ligament-contact forces and comparison with level walking

About a third of knee joint disorders originate from the patellofemoral (PF) site that makes stair ascent a difficult activity for patients. A detailed finite element model of the knee joint is coupled to a lower extremity musculoskeletal model to simulate the stance phase of stair ascent. It is driven by the mean of measurements on the hip-knee-ankle moments-angles as well as ground reaction forces reported in healthy individuals. Predicted muscle activities compare well to the recorded electromyography data. Peak forces in quadriceps (3.87 BW, body weight, at 20% instance in our 607 N subject), medial hamstrings (0.77 BW at 20%), and gastrocnemii (1.21 BW at 80%) are estimated. Due to much greater flexion angles-moments in the first half of stance, large PF contact forces (peak of 3.1 BW at 20% stance) and stresses (peak of 4.83 MPa at 20% stance) are estimated that exceed their peaks in level walking by fourfold and twofold, respectively. Compared with level walking, ACL forces diminish in the first half of stance but substantially increase later in the second half (peak of 0.76 BW at 75% stance). Under nearly similar contact forces at 20% of stance, the contact stress on the tibiofemoral (TF) medial plateau reaches a peak (9.68 MPa) twice that on the PF joint suggesting the vulnerability of both joints. Compared with walking, stair ascent increases peak ACL force and both peak TF and PF contact stresses. Reductions in the knee flexion moment and/or angle appear as a viable strategy to mitigate internal loads and pain.

The Influence of Mathematical Definitions on Patellar Kinematics Representations

A correlation between patellar kinematics and anterior knee pain is widely accepted. However, there is no consensus on how they are connected or what profile of patellar kinematics would minimize anterior knee pain. Nevertheless, answering this question by merging existing studies is further complicated by the variety of ways to describe patellar kinematics. Therefore, this study describes the most frequently used conventions for defining patellar kinematics, focusing on the rotations. The similarities and differences between the Cardan sequences and angles calculated by projecting axes are analyzed. Additionally, a tool is provided to enable the conversion of kinematic data between definitions in different studies. The choice of convention has a considerable impact on the absolute values and the clinical characteristics of the patello-femoral angles. In fact, the angles that result from using different mathematical conventions to describe a given patello-femoral rotation from our analyses differ up to a Root Mean Squared Error of 111.49° for patellar flexion, 55.72° for patellar spin and 35.39° for patellar tilt. To compare clinical kinematic patello-femoral results, every dataset must follow the same convention. Furthermore, researchers should be aware of the used convention’s implications to ensure reproducibility when interpreting and comparing such data.

Impact of Power Output on Muscle Activation and 3D Kinematics During an Incremental Test to Exhaustion in Professional Cyclists

This study aimed to quantify the influence of an increase in power output (PO) on joint kinematics and electromyographic (EMG) activity during an incremental test to exhaustion for a population of professional cyclists. The hip flexion/extension and internal/external rotation as well as knee abduction/adduction ranges of motion were significantly decreased at 100% of the maximal aerobic power (MAP). EMG analysis revealed a significant increase in the root mean square (RMS) for all muscles from 70% of the MAP. Gastrocnemius muscles [lateralis gastrocnemius (GasL) and medialis gastrocnemius (GasM)] were the less affected by the increase of PO. Cross-correlation method showed a significant increase in the lag angle values for VM in the last stage compared to the first stage, meaning that the onset of the activation started earlier during the pedaling cycle. Statistical Parametric Mapping (SPM) demonstrated that from 70% MAP, biceps femoris (BF), tibialis anterior (TA), gluteus maximus (GM), and rectus femoris (RF) yielded larger ranges of the crank cycle on which the level of recruitment was significantly increased. This study revealed specific muscular and kinematic coordination for professional cyclists in response to PO increase.

Assessment of quadriceps angle in children aged between 2 and 8 years

Aim:

The quadriceps angle is the angle between the line drawn from the spina iliaca anterior superior to the midpoint of the patella, and the line drawn from the midpoint of the patella to the tuberositas tibiae. It is important for lower extremity posture. The aim of this study was to determine the normative quadriceps angle value by measurement, and to assess the probable effect of factors such as measurement position, age, sex, and presence of pes planus on these values.

Material and Methods:

A total of 599 children consisting of 296 (49.4%) girls and 303 (50.6%) boys aged between 2 and 8 years, were included in the study. The children were divided into three groups by age as 2–4 years, 4–6 years, and 6–8-years. After the children’s demographic data were collected, the quadriceps angle was measured using an electronic goniometer. Pes planus was assessed by drawing the Feiss line.

Results:

In bilateral measurement, it was found that the quadriceps angle decreased with age both in the supine and standing positions (p<0.05). It was observed that sex and presence of pes planus had no effect on the quadriceps angle independent from measurement positions (p>0.05). A low negative correlation was found between body mass index and the quadriceps angle in both measurement positions (p<0.05).

Conclusion:

It was found that positional changes and weight bearing on limbs did not cause any change in knee position in healthy children. We consider that the decrease in quadriceps angle in this age group is due to growth rate asymmetry between the femur shaft and pelvic diameter.

Increased Patellar Fracture Rate in Total Knee Arthroplasty With Preoperative Varus Greater Than 15°: A Case-Control Study

BACKGROUND

Management of severe varus deformity requires soft tissue balancing for implantation of low-constraint knee prosthesis. Patellar complications have been rarely studied in this specific group. Our hypothesis was that severe genu varum (>15°) would increase the rate of patellar complications.

METHODS

Using a prospective cohort of 4216 prostheses performed at a single center beginning in 1987, we analyzed 280 prostheses having preoperative varus greater than 15°, compared to 673 total knee arthroplasties (TKAs) with a preoperative hip-knee-ankle angle of 180° ± 2°. Preoperative and postoperative clinical and radiological characteristics were compared between the 2 groups, with particular attention paid to patellar complications.

RESULTS

Average follow-up was 40.2 months (24-239). The mean preoperative Knee Society Score (KSS) was statistically higher in the normal (hip-knee-ankle angle 180° ± 2) axis group (62.65 vs 37.47, P = .001). At the last follow-up, no significant difference was found between the 2 groups in terms of postoperative KSS (87.5 in the varus group vs 87.3 in the normal axis group, P = .87). The rate of satisfied patients was identical between the 2 groups (85.3% vs 88.8%, P = .49). However, at mid-term, there were more patellar fractures in the varus group (2.9% vs 0.9%, P = .005). A significantly lower patellar height in both the varus group and the group of patella fractures (P < .001) was also found.

CONCLUSION

TKA in severe varus knees produces a KSS equivalent to TKA in knees with a mechanical axis of 0 ± 2. The risk of patellar fracture could initiate a decline in patella resurfacing in patients with major varus deformation, especially in case of a preoperative patella baja.

Patellar instability can be classified into four types based on patellar movement with knee flexion: a three-dimensional computer model analysis

Objective

Patellar instability (PI) represents various underlying pathologies, including patellar malalignment. Continuous patellar alignment develops to patellar tracking and is regarded as the end product of combined predisposing factors. We quantitatively investigated the inhomogeneity of patellar tracking in PI.

Methods

Sixty knees of 56 patients with PI and 15 knees of 10 healthy volunteers (HVs) were studied. Three-dimensional (3D) computer models were created based on MRIs at 10° intervals over 0°–50° of flexion, and patellar tracking was quantitatively analysed using patellar 3D shift. Classification was performed according to the maximum 3D shift (max-shift), indicating the extent of lateral deviation, and the change of 3D shift from 0° to 50° (change0–50), indicating movement direction. First, the cut-off value (COV) of the max-shift was defined based on the data from HVs. When a value was greater than the COV, it was defined as a major subluxation, and when the value was smaller it was defined as a minor subluxation. Next, the two COVs of change0–50 were similarly defined. When a value was greater than the upper COV, it was defined as a major-lateral type, laterally moving the patella with flexion, and when smaller than the lower COV it was defined as a major-medial type, medially moving the patella with flexion. When a value fell between the two COVs, it was defined as a major-straight type.

Results

Fifty-three patellae (88%) with values larger than the COV of the max-shift (mean +1 SD of HV) were defined as major subluxations and seven (12%) showing smaller values as minor subluxations. Among the major subluxations, 25 (42%) showing a smaller value than the lower COV of change0–50 (mean –2 SD of HV) were defined as major-medial type, while 7 (12%) showing a larger value than the upper COV of change0–50 (mean +2 SD) were defined as major-lateral type. Twenty-one (35%) were defined as major-straight type. No further analysis was performed on the seven minor subluxations (the minor type).

Conclusion

PI was quantitatively classified into four types according to the extent of lateral deviation and the movement direction of the patellae with flexion, showing inhomogeneity of patellar tracking.

Athletes with Brain Injury: Pathophysiologic and Medical Challenges

Participation in elite sporting activities is becoming increasingly popular for individuals with brain injury. This article outlines the types of brain injury and the associated movement dysfunctions. In addition, specific pathophysiologic and medical challenges facing athletes with brain injury are discussed. Further research conducted using athletes with brain injury will add to the existing literature indicating the benefits of athletic training in this population. Increased scientific study within this area stands to further improve understanding of the complex interaction between neuromuscular impairment and athletic performance.

The coupled effects of crouch gait and patella alta on tibiofemoral and patellofemoral cartilage loading in children

BACKGROUND

Elevated tibiofemoral and patellofemoral loading in children who exhibit crouch gait may contribute to skeletal deformities, pain, and cessation of walking ability. Surgical procedures used to treat crouch frequently correct knee extensor insufficiency by advancing the patella. However, there is little quantitative understanding of how the magnitudes of crouch and patellofemoral correction affect cartilage loading in gait.

METHODS

We used a computational musculoskeletal model to simulate the gait of twenty typically developing children and fifteen cerebral palsy patients who exhibited mild, moderate, and severe crouch. For each walking posture, we assessed the influence of patella alta and baja on tibiofemoral and patellofemoral cartilage contact.

RESULTS

Tibiofemoral and patellofemoral contact pressures during the stance phase of normal gait averaged 2.2 and 1.0 MPa. Crouch gait increased pressure in both the tibofemoral (2.6-4.3 MPa) and patellofemoral (1.8-3.3 MPa) joints, while also shifting tibiofemoral contact to the posterior tibial plateau. For extended-knee postures, normal patellar positions (Insall-Salvatti ratio 0.8-1.2) concentrated contact on the middle third of the patellar cartilage. However, in flexed knee postures, both normal and baja patellar positions shifted pressure toward the superior edge of the patella. Moving the patella into alta restored pressure to the middle region of the patellar cartilage as crouch increased.

CONCLUSIONS

This work illustrates the potential to dramatically reduce tibiofemoral and patellofemoral cartilage loading by surgically correcting crouch gait, and highlights the interaction between patella position and knee posture in modulating the location of patellar contact during functional activities.

Biphasic Analysis of Cartilage Stresses in the Patellofemoral Joint

The objective of this study was to examine the state of stress within the solid matrix of articular cartilage in the patellofemoral joint, using anatomically faithful biphasic models of the articular layers, with the joint subjected to physiologic muscle force magnitudes. Finite element models of five joints were created from human cadaver knees. Biphasic sliding contact analyses were performed using FEBio software to analyze the response of the joint from 30 to 60 degrees of knee flexion. Results demonstrated that the collagen matrix always sustains tensile stresses, despite the fact that the articular layers are loaded in compression. The principal direction of maximum solid stresses was consistent with the known orientation of collagen fibrils in cartilage. The magnitudes of these tensile stresses under muscle forces representative of activities of daily living were well below tensile failure stresses reported in the prior literature. Results also hinted that solid matrix stresses were higher in the patellar versus femoral superficial zone. These anatomically correct finite element models predicted outcomes consistent with our understanding of structure-function relationships in articular cartilage, while also producing solid matrix stress estimates not observable from experiments alone, yet highly relevant to our understanding of tissue degeneration.

A descriptive comparison of sprint cycling performance and neuromuscular characteristics in able-bodied athletes and paralympic athletes with cerebral palsy

OBJECTIVE

This study investigated the sprint cycling performance and neuromuscular characteristics of Paralympic athletes with cerebral palsy (CP) during a fatiguing maximal cycling trial compared with those of able-bodied (AB) athletes.

DESIGN

Five elite athletes with CP and 16 AB age- and performance-matched controls performed a 30-sec Wingate cycle test. Power output (W/kg) and fatigue index (%) were calculated. Electromyography was measured in five bilateral muscles and expressed in mean amplitude (mV) and median frequency (Hz).

RESULTS

Power output was significantly higher in the AB group (10.4 [0.5] W/kg) than in the CP group (9.8 [0.5] W/kg) (P < 0.05). Fatigue index was statistically similar between the AB (27% [0.1%]) and CP (25% [0.1%]) groups. Electromyographic mean amplitude and frequency changed similarly in all muscle groups tested, in both affected and nonaffected sides, in the CP and AB groups (P < 0.05). Neuromuscular irregularities were identified in the CP group.

CONCLUSIONS

The similarity in fatigue between the CP and AB groups indicates that elite athletes with CP may have a different exercise response to others with CP. The authors propose that this may result from high-level training over many years. This has rehabilitative implications, as it indicates near-maximal adaptation of the CP body toward normal levels.

Correlation between Intrinsic Patellofemoral Pain Syndrome in Young Adults and Lower Extremity Biomechanics

[Purpose] The purpose of this study was to evaluate the correlation between intrinsic patellofemoral pain syndrome (PFPS) in young adults and lower extremity biomechanics. [Subjects] This experiment was carried out with sixty (24 men and 32 women), who are normal university students as subjects. [Methods] All subjects underwent 3 clinical evaluations. For distinguishing the intrinsic PFPS from controls, we used the Modified Functional Index Questionnaire (MFIQ), Clarke’s test and the Eccentric step test. Based on the results of the tests, subjects who were classified as positive for 2 more tests were allocated to the bilateral or unilateral intrinsic PFPS group (n=14), and the others were allocated to the control group (n=42). These two groups were tested for hamstring tightness, foot overpronation, and static Q-angle and dynamic Q-angle. These are the four lower extremity biomechanic, cited as risk factors of patellofemoral pain syndrome. [Results] The over pronation, static Q-angle and the dynamic Q-angle were not significantly different between the two groups. However, the hamstring tightness of the PFPS group was significantly greater than that of the controls. [Conclusion] We examined individuals for intrinsic patellofemoral pain syndrome in young adults and lower extremity biomechanics. We found a strong correlation between intrinsic PFPS and hamstring tightness.

A geometric approach to study the contact mechanisms in the patellofemoral joint of normal versus patellofemoral pain syndrome subjects

The biomechanics of the patellofemoral (PF) joint is complex in nature, and the aetiology of such manifestations of PF instability as patellofemoral pain syndrome (PFPS) is still unclear. At this point, the particular factors affecting PFPS have not yet been determined. This study has two objectives: (1) The first is to develop an alternative geometric method using a three-dimensional (3D) registration technique and linear mapping to investigate the PF joint contact stress using an indirect measure: the depth of virtual penetration (PD) of the patellar cartilage surface into the femoral cartilage surface. (2) The second is to develop 3D PF joint models using the finite element analysis (FEA) to quantify in vivo cartilage contact stress and to compare the peak contact stress location obtained from the FE models with the location of the maximum PD. Magnetic resonance images of healthy and PFPS subjects at knee flexion angles of 15°, 30° and 45° during isometric loading have been used to develop the geometric models. The results obtained from both approaches demonstrated that the subjects with PFPS show higher PD and contact stresses than the normal subjects. Maximum stress and PD increase with flexion angle, and occur on the lateral side in healthy and on the medial side in PFPS subjects. It has been concluded that the alternative geometric method is reliable in addition to being computationally efficient compared with FEA, and has the potential to assess the mechanics of PFPS with an accuracy similar to the FEA.

Recent advances in computational mechanics of the human knee joint

Computational mechanics has been advanced in every area of orthopedic biomechanics. The objective of this paper is to provide a general review of the computational models used in the analysis of the mechanical function of the knee joint in different loading and pathological conditions. Major review articles published in related areas are summarized first. The constitutive models for soft tissues of the knee are briefly discussed to facilitate understanding the joint modeling. A detailed review of the tibiofemoral joint models is presented thereafter. The geometry reconstruction procedures as well as some critical issues in finite element modeling are also discussed. Computational modeling can be a reliable and effective method for the study of mechanical behavior of the knee joint, if the model is constructed correctly. Single-phase material models have been used to predict the instantaneous load response for the healthy knees and repaired joints, such as total and partial meniscectomies, ACL and PCL reconstructions, and joint replacements. Recently, poromechanical models accounting for fluid pressurization in soft tissues have been proposed to study the viscoelastic response of the healthy and impaired knee joints. While the constitutive modeling has been considerably advanced at the tissue level, many challenges still exist in applying a good material model to three-dimensional joint simulations. A complete model validation at the joint level seems impossible presently, because only simple data can be obtained experimentally. Therefore, model validation may be concentrated on the constitutive laws using multiple mechanical tests of the tissues. Extensive model verifications at the joint level are still crucial for the accuracy of the modeling.

The effect of medial opening wedge proximal tibial osteotomy on patellofemoral contact

BACKGROUND

It has been suggested that patellofemoral contact pressures and contact forces may be altered secondary to an opening wedge high tibial osteotomy, yet few data are available that quantify the effect of varying degrees of medial opening wedge osteotomy on the patellofemoral joint contact characteristics.

HYPOTHESIS

Opening wedge medial proximal tibial osteotomy will increase patellofemoral contact force and pressure.

STUDY DESIGN

Controlled laboratory study.

METHODS

Nine human cadaver knees were used. Pressure-sensitive film was placed in the suprapatellar pouch, leaving the patellar tendon and medial and lateral retinacula intact. The quadriceps tendon was attached to a materials testing machine along the axis of the femur, whereby a pulley mechanism generated 950 N of force. Patellofemoral contact characteristics were measured with pressure-sensitive film at 30°, 60°, 90°, and 120° of flexion for the native knee and after subsequent 10-mm and 15-mm medial opening wedge proximal tibial osteotomies. The film was analyzed with imaging software.

RESULTS

There was a statistically significant increase (P < .05) in mean contact pressure at 30° and 120° between the 10-mm osteotomy and native knee and across all flexion angles between the 15-mm osteotomy and native knee. Furthermore, a significant difference was seen in peak pressures when native knees were compared with 10-mm and 15-mm opening wedge osteotomies at all flexion angles.

CONCLUSION

There was a significant increase in patellofemoral pressures at varying degrees of knee flexion after medial opening wedge proximal tibial osteotomies of only 10 mm; a larger osteotomy resulted in a greater increase.

CLINICAL RELEVANCE

When performing a medial opening wedge proximal tibial osteotomy, the surgeon should consider the negative effects of increased patellofemoral peak pressure.

A detailed and validated three dimensional dynamic model of the patellofemoral joint

A detailed 3D anatomical model of the patellofemoral joint was developed to study the tracking, force, contact and stability characteristics of the joint. The quadriceps was considered to include six components represented by 15 force vectors. The patellar tendon was modeled using four bundles of viscoelastic tensile elements. Each of the lateral and medial retinaculum was modeled by a three-bundle nonlinear spring. The femur and patella were considered as rigid bodies with their articular cartilage layers represented by an isotropic viscoelastic material. The geometrical and tracking data needed for model simulation, as well as validation of its results, were obtained from an in vivo experiment, involving MR imaging of a normal knee while performing isometric leg press against a constant 140 N force. The model was formulated within the framework of a rigid body spring model and solved using forth-order Runge-Kutta, for knee flexion angles between zero and 50 degrees. Results indicated a good agreement between the model predictions for patellar tracking and the experimental results with RMS deviations of about 2 mm for translations (less than 0.7 mm for patellar mediolateral shift), and 4 degrees for rotations (less than 3 degrees for patellar tilt). The contact pattern predicted by the model was also consistent with the results of the experiment and the literature. The joint contact force increased linearly with progressive knee flexion from 80 N to 210 N. The medial retinaculum experienced a peak force of 18 N at full extension that decreased with knee flexion and disappeared entirely at 20 degrees flexion. Analysis of the patellar time response to the quadriceps contraction suggested that the muscle activation most affected the patellar shift and tilt. These results are consistent with the recent observations in the literature concerning the significance of retinaculum and quadriceps in the patellar stability.

Simulation of patella alta and the implications for in vitro patellar tracking in the ovine stifle joint

Patella alta is associated with adverse cartilage adaptations, patellofemoral pain, and instability. It is defined by a relatively long patellar tendon and patella positioned in a more proximal location within the patellar groove of the femur. This study used the ovine stifle joint model to investigate the effect of patellar tendon lengthening on the 3D passive kinematics of the patellofemoral and tibiofemoral joints. Eight patellar tendons were lengthened in 2 mm increments up to a maximum of 12 mm (20%) using a device placed in series with the transected patellar tendon. Three-dimensional kinematics were measured in the intact joint and at each increment of patellar tendon length (L(T)) during passively induced tibiofemoral flexion. Patellar flexion angle was linearly correlated with tibial flexion angle in the intact joint, and this correlation persisted after tendon lengthening (R = 0.897-0.965, p < 0.01). Patellofemoral kinematics expressed as a function of tibial flexion angle were significantly altered by L(T) increases >9%. In contrast, when patellofemoral kinematics were expressed as a function of patellar flexion angle they were not significantly altered by increases in L(T). Tibiofemoral kinematics were not affected by the L(T) increases. These results demonstrate that for a given tibial flexion angle, patellar tendon lengthening alters the patellar flexion angle. However, for a given patellar flexion angle, the orientation of the patella in the remaining five degrees of freedom is unchanged, implying a repeatable path of patellar motion.

The development of lower limb musculoskeletal models with clinical relevance is dependent upon the fidelity of the mathematical description of the lower limb. Part 2: patient-specific geometry

Musculoskeletal models have the potential to evolve into sensitive clinical tools that provide relevant therapeutic guidance. A key impediment to this is the lack of understanding as to the function of such models. In order to improve this it is useful to recognise that musculoskeletal modelling is the mathematical description of musculoskeletal movement – a process that involves the construction and solution of equations of motion. These equations are derived from standard mechanical considerations and the mathematical representation of anatomy. The fidelity of musculoskeletal models is highly dependent on the assumption that such representations also describe the function of the musculoskeletal geometry. In addition, it is important to understand the sensitivity of such representations to patient-specific variations in anatomy. The exploration of these twin considerations will be fundamental to the creation of musculoskeletal modelling tools with clinical relevance and a systematic enquiry of these key parameters is recommended.

MORPHOMETRICAL COMPARISON BETWEEN THE RESECTED SURFACES IN OSTEOARTHRITIC KNEES AND POROUS-COATED ANATOMIC KNEE PROSTHESIS

In total knee replacement, a good match of the prosthesis to bone is very important. Most knees that require total knee replacement are deformed. However, most of the design parameters of knee prosthesis were based on the normal knee. In this series, the dimensions of the resected surfaces in 77 osteoarthritic knees were measured intraoperatively and compared with the corresponding surfaces of the porous-coated anatomic (PCA) (Howmedica, Rutherford, NJ, USA) knee prosthesis. The results showed that the medial femoral condyle was wider than the lateral femoral condyle (p < 0.05) in the resected surfaces. The intercondylar notch of the resected femur was wider than that of the prosthesis (p < 0.05). In the resected tibial plateau, the ratio of the anteroposterior length to the mediolateral width was larger than that of the prosthesis (p < 0.05). The length and width of the resected patella were greater than those of the implant (p < 0.05). According to the difference in morphometrical parameters between the resected surfaces and the knee prosthesis, we suggest that the dimensions of the resected surfaces of the osteoarthritic knee should be important design parameters in total knee prosthesis.

Alloplastic reconstruction of the extensor mechanism after resection of tibial sarcoma

Reconstruction of the extensor mechanism is essential for good extremity function after endoprosthetic knee replacement following tumor resection. Only a few biological methods have been able to reliably restore a functional extensor mechanism, but they are often associated with significant complication rates. Reattachment of the patellar tendon to the prosthesis using an alloplastic patellar ligament (Trevira cord) can be an appropriate alternative. In vivo and in vitro studies have already shown that complete fibrous ingrowth in polyethylene chords can be seen after a period of six months. However, until now, no biomechanical study has shown the efficacy of an alloplastic cord and its fixation device in providing sufficient stability and endurance in daily life-activity until newly formed scar tissue can take over this function. In a special test bench developed for this study, different loading regimes were applied to simulate loads during everyday life. Failure loads and failure modes were evaluated. The properties of the cord were compared before and after physiological conditioning. It was shown that rubbing was the mode of failure under dynamic loading. Tensile forces up to 2558 N did not result in material failure. Thus, using an artificial cord together with this fixation device, temporary sufficient stable fixation can be expected.

Contact forces in several TKA designs during squatting: A numerical sensitivity analysis

Total knee arthroplasty (TKA) is a very successful procedure, but pain or difficulties during activities still persist in patients. Patient outcomes in TKA surgery can be affected by implant design, alignment or patient-related anatomical factors. This paper presents a numerical sensitivity analysis of several TKA types: a fixed bearing, posterior stabilized prosthesis, a high flexion fixed bearing guided motion prosthesis, a mobile bearing prosthesis and a hinge prosthesis. Each prosthesis was virtually implanted on the same cadaver leg model and it underwent a loaded squat, in 10s, between 0° and 120°, similar to several previous experimental tests performed on knee kinematics simulators. The aim of this examination was to investigate the sensitivity of the patello-femoral (PF) and tibio-femoral (TF) contact forces to patient-related anatomical factors, and component position in the different implant types. The following parameters were used for the sensitivity study: the proximo-distal patellar position, the patellar component tilting, the tibial component position and orientation, the locations of the medial and lateral collateral ligaments with respect to femur and tibia and the patellar tendon length. The sensitivity analysis showed that PF contact forces are mostly affected by patella height (increases up to 67% for one TKA type in patella-alta configuration), by an anterior tibial component translation (increases up to 30%), and by patellar component tilting (increases up to 29%); TF contact forces are mostly affected by the anterior displacement of the insertion points of the medial collateral ligament with respect to the reference position (increases up to 48%).

Quantification of Patellofemoral Joint Reaction Forces during Functional Activities Using a Subject-Specific Three-Dimensional Model

The purpose of this study was to describe an imaging based, subject specific model that was developed to quantify patellofemoral joint reaction forces (PFJRF’s). The secondary purpose was to test the model in a group of healthy individuals while performing various functional tasks. Twenty healthy subjects (10 males, 10 females) were recruited. All participants underwent two phases of data collection: 1) magnetic resonance imaging of the knee, patellofemoral joint, and thigh, and 2) kinematic, kinetic and EMG analysis during walking, running, stair ascent, and stair descent. Using data obtained from MRI, a subject specific representation of the extensor mechanism was created. Individual gait data were used to drive the model (via an optimization routine) and three-dimensional vasti muscle forces and subsequent three-dimensional PFJRF’s were computed. The average peak PFJRF was found to be highest during running (58.2 N/kg-bwt), followed by stair ascent (33.9 N/kg-bwt), stair descent (27.9 N/kg-bwt), and walking (10.1 N/kg-bwt). No differences were found between males and females. For all conditions, the direction of the PFJRF was always in the posterior, superior, and lateral directions. The posterior component of the PFJRF always had the greatest magnitude, followed by superior and lateral components. Our results indicate that estimates of the magnitude and direction of the PFJRF during functional tasks can be obtained using a 3D-imaging based model.

Is there a biomechanical explanation for anterior knee pain in patients with patella alta?: influence of patellar height on patellofemoral contact force, contact area and contact pressure

The purpose of this study was to test the hypothesis that patella alta leads to a less favourable situation in terms of patellofemoral contact force, contact area and contact pressure than the normal patellar position, and thereby gives rise to anterior knee pain. A dynamic knee simulator system based on the Oxford rig and allowing six degrees of freedom was adapted in order to simulate and record the dynamic loads during a knee squat from 30 degrees to 120 degrees flexion under physiological conditions. Five different configurations were studied, with variable predetermined patellar heights. The patellofemoral contact force increased with increasing knee flexion until contact occurred between the quadriceps tendon and the femoral trochlea, inducing load sharing. Patella alta caused a delay of this contact until deeper flexion. As a consequence, the maximal patellofemoral contact force and contact pressure increased significantly with increasing patellar height (p < 0.01). Patella alta was associated with the highest maximal patellofemoral contact force and contact pressure. When averaged across all flexion angles, a normal patellar position was associated with the lowest contact pressures. Our results indicate that there is a biomechanical reason for anterior knee pain in patients with patella alta.

Computational modeling to predict mechanical function of joints: application to the lower leg with simulation of two cadaver studies

Computational models of musculoskeletal joints and limbs can provide useful information about joint mechanics. Validated models can be used as predictive devices for understanding joint function and serve as clinical tools for predicting the outcome of surgical procedures. A new computational modeling approach was developed for simulating joint kinematics that are dictated by bone/joint anatomy, ligamentous constraints, and applied loading. Three-dimensional computational models of the lower leg were created to illustrate the application of this new approach. Model development began with generating three-dimensional surfaces of each bone from CT images and then importing into the three-dimensional solid modeling software SOLIDWORKS and motion simulation package COSMOSMOTION. Through SOLIDWORKS and COSMOSMOTION, each bone surface file was filled to create a solid object and positioned necessary components added, and simulations executed. Three-dimensional contacts were added to inhibit intersection of the bones during motion. Ligaments were represented as linear springs. Model predictions were then validated by comparison to two different cadaver studies, syndesmotic injury and repair and ankle inversion following ligament transection. The syndesmotic injury model was able to predict tibial rotation, fibular rotation, and anterior/posterior displacement. In the inversion simulation, calcaneofibular ligament extension and angles of inversion compared well. Some experimental data proved harder to simulate accurately, due to certain software limitations and lack of complete experimental data. Other parameters that could not be easily obtained experimentally can be predicted and analyzed by the computational simulations. In the syndesmotic injury study, the force generated in the tibionavicular and calcaneofibular ligaments reduced with the insertion of the staple, indicating how this repair technique changes joint function. After transection of the calcaneofibular ligament in the inversion stability study, a major increase in force was seen in several of the ligaments on the lateral aspect of the foot and ankle, indicating the recruitment of other structures to permit function after injury. Overall, the computational models were able to predict joint kinematics of the lower leg with particular focus on the ankle complex. This same approach can be taken to create models of other limb segments such as the elbow and wrist. Additional parameters can be calculated in the models that are not easily obtained experimentally such as ligament forces, force transmission across joints, and three-dimensional movement of all bones. Muscle activation can be incorporated in the model through the action of applied forces within the software for future studies.

The effects of trochlear groove geometry on patellofemoral joint stability-a computer model study

The effect of the variation in the femoral groove geometry on patellofemoral joint stability was studied using a two-dimensional transverse plane model with deformable articular surfaces. The femoral and patellar bony structures were modelled as rigid bodies with their profiles expressed by splines. The articular cartilage was discretized into compression springs, distributed along the femoral and patellar profiles, based on the rigid-body spring model. The medial and lateral retinacula were modelled as linear tensile springs, and the quadriceps muscles and patellar tendon as strings with known tension. The anatomical data were obtained from the transverse plane magnetic resonance images of a normal knee flexed at 20° and from the literature. A dynamic analysis approach was employed to solve the governing equations of the model, i.e. three static equilibrium equations of the patella and a constraint equation for each cartilage spring, explicitly. The results of the model suggest that alteration of the sulcus angle from 139° to 169° causes a lateral shift and tilt of less than 3 mm and 4°. This effect increased slightly with increasing total quadriceps force, however, to significantly more than 7 mm and 18° respectively when the medial retinaculum was released. It was suggested that this might be the combined effect of the medial retinaculum deficiency and trochlear dysplasia that is responsible for patellar subluxation and, particularly, dislocation disorders.

Electromyographic analysis of pedaling: a review

Although pedaling is constrained by the circular trajectory of the pedals, it is not a simple movement. This review attempts to provide an overview of the pedaling technique using an electromyographic (EMG) approach. Literature concerning the electromyographic analysis of pedaling is reviewed in an effort to make a synthesis of the available information, and to point out its relevance for researchers, clinicians and/or cycling/triathlon trainers. The first part of the review depicts methodological aspects of the EMG signal recording and processing. We show how the pattern of muscle activation during pedaling can be analyzed in terms of muscle activity level and muscle activation timing. Muscle activity level is generally quantified with root mean square or integrated EMG values. Muscle activation timing is studied by defining EMG signal onset and offset times that identify the duration of EMG bursts and, more recently, by the determination of a lag time maximizing the cross-correlation coefficient. In the second part of the review, we describe whether the patterns of the lower limb muscles activity are influenced by numerous factors affecting pedaling such as power output, pedaling rate, body position, shoe-pedal interface, training status and fatigue. Some research perspectives linked to pedaling performance are discussed throughout the manuscript and in the conclusion.

Surgical biomechanics of the patellofemoral joint

This review presents objective data, as far as possible, about the current understanding of the biomechanics of the patellofemoral joint as it pertains to the management of patellofemoral problems. When faced with a patellofemoral malfunction, it is important to check all the soft-tissue and articular geometry factors relating to the patella locally and not to neglect the overall lower limb alignment and function. It is important to remember that small alterations in alignment can result in significant alterations in patellofemoral joint stresses and that changes in the mechanics of the patellofemoral joint can also result in changes in the tibiofemoral compartments. Surgical intervention for patellofemoral problems needs to be planned carefully and take into account an individual's anatomy.

A rigid body spring model to investigate the lateral shift - restraining force behavior of the patella

Patellar lateral stability was studied using a 2D transverse plane model with deformable articular surfaces. Quadriceps muscles and patellar tendon were considered as strings with predefined forces and lateral and medial retinaculum as tensile springs. Deformation behavior of articular cartilage was modeled by a set of compression springs perpendicular to articular surfaces, based on rigid body spring model method (RBSM). Patellar lateral stability was investigated using restraining force method (the external force required to cause up to 10 mm lateral displacement on patella). The results were in good agreement with experimental reports for normal joint, vastus lateralis and vastus medialis relieved. Small changes in the femoral trochlear groove geometry provided significant variation in patellar stability. Simulation of different surgical treatments showed that the tibial tubercle medialization is the most effective procedure for patellar subluxation and dislocation disorders.

Interaction of ligament bundles and articular contacts for the simulation of passive knee flexion

The purpose of this study is to investigate the effects of anterior bundle of ACL (aACL), anterior portion of PCL (aPCL), anterior and deep portions of MCL (aMCL, dMCL) and the tibio-femoral articular contacts on to passive knee motion. A three-dimensional simplistic anatomical dynamic model, based on the literature was used as a reference. This reference model attaches the bundles of the ligaments on medial and the lateral spherical condyles of the femur and tibial plateau giving us a representation close enough to a normal natural tibio-femoral joint, but does not allow to study abnormalities of the knee kinematics due to the assumptions of the femur shape. The proposed three-dimensional dynamic tibio-femoral model, however includes the isometric fascicles, aACL, aPCL, aMCL, dMCL, and irregularly shaped medial-lateral contact surfaces. The approach taken in this model is capable of ligament and bone surface modifications that will enable us to analyze bone shape and ligament related abnormalities of knee kinematics.

A new approach to C2 continuous piecewise bicubic representation of the articular surfaces of diarthrodial joints

Based on the force-deflection equation for a beam subjected to lateral point loads, a C2 continuous piecewise bicubic mathematical representation was proposed to model complicated geometrical surfaces, e.g. the articular surfaces of human joints. The method was then extended so that it could be used for mathematical modelling of incomplete nets of data points, as well as smoothing of noisy and/or filtering of erroneous data points. Mathematical techniques were also developed to calculate the required unknown parameters explicitly, with no need to solve the system of equations simultaneously. The performance of the proposed method was evaluated on a number of surface modelling problems, including two known analytical surfaces and the human femoral and patellar articular surfaces. The results indicate that the proposed method is precise, flexible, and easy to apply and has several advantages over the conventional smoothing methods, i.e. the B-spline approach.

The effect of modeling cartilage on predicted ligament and contact forces at the knee

This paper describes a two-dimensional dynamic model of the human knee joint by employing cartilage model. The model outputs of rigid and deformable models are compared. The predicted tibio-femoral contact and ligament forces remain consistent whether or not a model of cartilage was included in the calculations. For both rigid and deformable contact models, the maximum contact force decreases as the amplitude of external force is increased. The effect of cartilage on ligaments' responses is insignificant except for the behavior of the anterior portion of the anterior cruciate ligament at large flexion angles.

Biomechanics of the knee joint in flexion under various quadriceps forces

Bioemchanics of the entire knee joint including tibiofemoral and patellofemoral joints were investigated at different flexion angles (0 degrees to 90 degrees ) and quadriceps forces (3, 137, and 411 N). In particular, the effect of changes in location and magnitude of restraining force that counterbalances the isometric extensor moment on predictions was investigated. The model consisted of three bony structures and their articular cartilage layers, menisci, principal ligaments, patellar tendon, and quadriceps muscle. Quadriceps forces significantly increased the anterior cruciate ligament, patellar tendon, and contact forces/areas as well as the joint resistant moment. Joint flexion, however, substantially diminished them all with the exception of the patellofemoral contact force/area that markedly increased in flexion. When resisting extensor moment by a force applied on the tibia, the force in cruciate ligaments and tibial translation significantly altered as a function of magnitude and location of the restraining force. Quadriceps activation generated large ACL forces at full extension suggesting that post ACL reconstruction exercises should avoid large quadriceps exertions at near full extension angles. In isometric extension exercises against a force on the tibia, larger restraining force and its more proximal location to the joint substantially decreased forces in the anterior cruciate ligament at small flexion angles whereas they significantly increased forces in the posterior cruciate ligament at larger flexion angles.

The influence of patellofemoral joint contact geometry on the modeling of three dimensional patellofemoral joint forces

The purpose of this study was to determine the influence of patellofemoral joint contact geometry on the modeling of three-dimensional patellofemoral joint forces. To achieve this goal, patellofemoral joint reaction forces (PFJRFs) that were measured from an in-vitro cadaveric set-up were compared to PFJRFs estimated from a computer model that did not consider patellofemoral joint contact geometry. Ten cadaver knees were used in this study. Each was mounted on a custom jig that was fixed to an Instron frame. Quadriceps muscle loads were accomplished using a pulley system and weights. The force in the patellar ligament was obtained using a buckle transducer. To quantify the magnitude and direction of the PFJRF, a six-axis load cell was incorporated into the femoral fixation system so that a rigid body assumption could be made. PFJRF data were obtained at 0 degrees , 20 degrees , 40 degrees and 60 degrees of knee flexion. Following in vitro testing, SIMM modeling software was used to develop computational models based on the three-dimensional coordinates (Microscribe digitizer) of individual muscle and patellar ligament force vectors obtained from the cadaver knees. The overall magnitude of the PFJRF estimated from the computer generated models closely matched the direct measurements from the in vitro set-up (Pearson's correlation coefficient, R(2)=0.91, p<0.001). Although the computational model accurately estimated the posteriorly directed forces acting on the joint, some discrepancies were noted in the forces acting in the superior and lateral directions. These differences however, were relatively small when expressed as a total of the overall PFJRF magnitude.

Relationship of Q angle and joint hypermobility and Q angle values in different positions

Patellar malalignment is the most important reason for anterior knee pain. Patellar alignment is commonly determined by the measurement of the quadriceps (Q) angle. In this study, our primary aim was to investigate the Q angle difference between healthy individuals with and without joint hypermobility. The other objectives were to compare the Q angle values, which were measured in supine and upright positions, to determine hypermobility frequency among healthy males in a certain population, and to investigate the correlation between the existent skeletal deformities and Beighton score levels. Two hundred fifty-three healthy male college students were examined for hypermobility and skeletal deformities. According to the Beighton scoring system, three groups (n=20) were formed, and Q angle measurements were performed on the 60 individuals in both supine and upright positions. In the comparison of the groups, the mean Q angle values in healthy hypermobile individuals were found to be significantly higher than that of the nonhypermobile ones (p<0.05). No statistical difference was found between Q angle values in supine and upright positions (p>0.05). The frequency of joint hypermobility (Beighton score 4 or more) was found to be 29.25% in this population. No correlation was determined between existent skeletal deformities and Beighton score values (p>0.05). In conclusion, the Q angle evaluation among healthy hypermobile individuals may have a prognostic value for probable knee pathologies that may appear in the future.

Computational Modeling of a Dynamic Knee Simulator for Reproduction of Knee Loading

As a first step towards reproducing desired three-dimensional joint loading and motion on a dynamic knee simulator, the goal of this study was to develop and verify a three-dimensional computational model that generated control profiles for the simulator using desired knee loading and motion as model inputs. The developed model was verified by predicting tibio-femoral loading on an instrumented analog knee for given actuator forces and the ability to generate simulator control profiles was demonstrated using a three-dimensional walking profile. The model predicted axial tibia loading for a sagittal-plane dual-limb squat within 1% of measured peak loading. Adding out-of-sagittal-plane forces decreased the accuracy of load prediction. The model generated control profiles to the simulator that produced axial tibia loading within 16% of desired for walking. Discrepancies in predicted and measured quadriceps forces influenced the accuracy of the generated control profiles. Future work will replace the analog knee in both the model and machine with a prosthetic knee.

An anatomically based patient-specific finite element model of patella articulation: towards a diagnostic tool

A 3D anatomically based patient-specific finite element (FE) model of patello-femoral (PF) articulation is presented to analyse the main features of patella biomechanics, namely, patella tracking (kinematics), quadriceps extensor forces, surface contact and internal patella stresses. The generic geometries are a subset from the model database of the International Union of Physiological Sciences (IUPS) (http://www.physiome.org.nz) Physiome Project with soft tissue derived from the widely used visible human dataset, and the bones digitised from an anatomically accurate physical model with muscle attachment information. The models are customised to patient magnetic resonance images using a variant of free-form deformation, called 'host-mesh' fitting. The continuum was solved using the governing equation of finite elasticity, with the multibody problem coupled through contact mechanics. Additional constraints such as tissue incompressibility are also imposed. Passive material properties are taken from the literature and implemented for deformable tissue with a non-linear micro-structurally based constitutive law. Bone and cartilage are implemented using a 'St-Venant Kirchoff' model suitable for rigid body rotations. The surface fibre directions have been estimated from anatomy images of cadaver muscle dissections and active muscle contraction was based on a steady-state calcium-tension relation. The 3D continuum model of muscle, tendon and bone is compared with experimental results from the literature, and surgical simulations performed to illustrate its clinical assessment capabilities (a Maquet procedure for reducing patella stresses and a vastus lateralis release for a bipartite patella). Finally, the model limitations, issues and future improvements are discussed.

Computational modelling of a total knee prosthetic loaded in a dynamic knee simulator

Dynamic knee simulators attempt to reproduce the estimated forces, moments, and motions of both the patello-femoral and tibio-femoral joints during ambulatory activities. As a continuation of work designed to reproduce desired three-dimensional joint loading and motion on a dynamic knee simulator, the goal of this study was to develop a computational model of a prosthetic knee placed within an existing computational model of a dynamic knee simulator. The resulting model was then used to produce inputs to the controllable axes of the simulator for reproduction of desired knee loading and motion. Previously, a three-dimensional computational model of the simulator was developed and verified using a simplified and instrumented analog knee. The work presented here replaced the simplified and constrained geometries of the analog knee with structures representing a prosthetic knee. Three-dimensional geometries were determined based on digitized surface points of a right total knee replacement. Deformable contacts between the articulating surfaces of the tibio-femoral and patello-femoral joints were then defined and model sensitivities were identified. Predicted results from the computational model were compared to experimental results for force profiles applied at the simulator's controllable axes. Within identified limitations, the model was then used to generate inputs to the simulator to reproduce desired patellar tendon load during a squat and desired out-of-sagittal-plane motion during a squat.

In vivo and noninvasive load sharing among the vasti in patellar malalignment

PURPOSE

It is not clear how the knee extension torque is distributed quantitatively among the lateral and medial vasti in patellofemoral pain (PFP) patients with patellar malalignment, which was investigated in vivo and noninvasively in ten PFP patients and eleven controls. We hypothesized that the vastus medialis oblique (VMO) and vastus medialis longus (VML) of PFP patients contribute less to knee extension than that in controls.

METHODS

Electrical stimulation was used to activate each vastus component selectively. The relationship between the knee extension torque generated by each individual vastus component and the corresponding compound muscle action potential (M-wave) was established over different contraction levels, which was used to calibrate the corresponding voluntary EMG signal and determine torque ratios of VMO/VL (vastus lateralis), VMO/VML, VML/VL and (VMO+VML)NL during voluntary isometric knee extension.

RESULTS

The VMO and VML of PFP patients contributed significantly less to knee extension than their counterparts in controls. The combination of VMO and VML generated comparable amount of extension torque as the VL in the controls, while it produced significantly lower extension torque than that of the VL in the PFP patients. In addition, the VMO/VL was lower than VMO/VML and VML/VL in both PFP and control groups.

CONCLUSIONS

Compared to controls, the VMO and VML in the PFP patients contributed significantly less to the knee extension torque. The approach can be used to investigate load sharing among quadriceps muscles in vivo and noninvasively, in both healthy subjects and patients with patellofemoral disorder and patellar malalignment.

Simultaneous In Vitro Measurement of Patellofemoral Kinematics and Forces

This study involved the development and testing of a system for the simultaneous in vitro measurement of tibiofemoral kinematics and patellofemoral kinematics and forces. Knee motion was tracked using a Vicon 370, and patellofemoral force was measured using a six degree-of-freedom transducer based on the design of Singerman et al. Using this system, twelve knee specimens were tested in supine leg extension under a simulated quadriceps force. The comprehensive set of results corresponds well to the individual results of previous studies. The measurement system will be of value in assessing the effects of total knee arthroplasty on patellar function.

Evaluation of a computational model used to predict the patellofemoral contact pressure distribution

One possible cause of patellofemoral pain syndrome is excessive lateral force acting on the patella. Although several treatment methods focus on decreasing the lateral force acting on the patella, the relationship between the lateral force and the patellofemoral contact pressure distribution is unclear. A computational model has been developed to determine how loading variations alter the patellofemoral force and pressure distributions for individual knees. The model allows variation in the quadriceps and patella tendon forces, and calculates the predicted contact pressure distribution using the discrete element analysis technique. To characterize the accuracy of the model, four cadaver knees were flexed on a knee simulator with three initial Q-angles, while recording the force and pressure distributions with a pressure sensor. A model of each knee was created from CT data. Using the external force applied to the knee, the geometry of the knee, and the quadriceps origin as input, the pressure distribution was calculated during flexion. Similar trends were noted for the computational and experimental results. The percentage of the total force applied to the lateral cartilage increased with the Q-angle. The maximum contact pressure increased during flexion. The maximum lateral contact pressure increased with the Q-angle for three knees. For the other knee, increasing the Q-angle decreased the maximum lateral pressure. The maximum medial contact pressure decreased as the Q-angle increased. By characterizing the influence of patellofemoral loading on the force and pressure distributions, the computational model could be used to evaluate treatment methods prescribed for patellofemoral pain.

3-D Anatomically Based Dynamic Modeling of the Human Knee to Include Tibio-Femoral and Patello-Femoral Joints

An anatomical dynamic model consisting of three body segments, femur, tibia and patella, has been developed in order to determine the three-dimensional dynamic response of the human knee. Deformable contact was allowed at all articular surfaces, which were mathematically represented using Coons’ bicubic surface patches. Nonlinear elastic springs were used to model all ligamentous structures. Two joint coordinate systems were employed to describe the six-degrees-of-freedom tibio-femoral (TF) and patello-femoral (PF) joint motions using twelve kinematic parameters. Two versions of the model were developed to account for wrapping and nonwrapping of the quadriceps tendon around the femur. Model equations consist of twelve nonlinear second-order ordinary differential equations coupled with nonlinear algebraic constraint equations resulting in a Differential-Algebraic Equations (DAE) system that was solved using the D_ifferential/A_lgebraic S_ystem S_ol_ver (DASSL) developed at Lawrence Livermore National Laboratory. Model calculations were performed to simulate the knee extension exercise by applying non-linear forcing functions to the quadriceps tendon. Under the conditions tested, both “screw home mechanism” and patellar flexion lagging were predicted. Throughout the entire range of motion, the medial component of the TF contact force was found to be larger than the lateral one while the lateral component of the PF contact force was found to be larger than the medial one. The anterior and posterior fibers of both anterior and posterior cruciate ligaments, ACL and PCL, respectively, had opposite force patterns: the posterior fibers were most taut at full extension while the anterior fibers were most taut near 90° of flexion. The ACL was found to carry a larger total force than the PCL at full extension, while the PCL carried a larger total force than the ACL in the range of 75° to 90° of flexion.

Computer simulations of patellofemoral joint surgery. Patient-specific models for tuberosity transfer

BACKGROUND

Variable clinical outcomes of tibial tuberosity transfer surgery have been reported.

HYPOTHESES

The biomechanical outcome of surgery is patient-specific; no single procedure produces superior results for all patients. Use of patient-specific computer models can optimize choice of procedure.

STUDY DESIGN

Computer simulation study using clinical data.

METHODS

We used patient-specific multibody models of the patellofemoral joints of 20 patients with a diagnosis of patellar subluxation and osteoarthritis. Four tibial tuberosity transfer procedures (two anterior and two anteromedial) were simulated for each patient and compared with their preoperative model.

RESULTS

When results for all patients were averaged, all simulated operations produced a statistically significant decrease in surface-wide mean contact stress, although no significant difference was found among them.

CONCLUSIONS

The simulated surgical outcomes were patient-specific: no single procedure was consistently superior at decreasing peak or mean stress and each procedure produced a potentially detrimental outcome, an increase in either mean stress or peak stress, in at least one patient.

CLINICAL RELEVANCE

Computer simulation may serve as a valuable tool for tailoring procedures to specific patients.

The effect of vastus medialis forces on patello-femoral contact: a model-based study

A mathematical model of the patello-femoral joint was introduced to investigate the impact of the vastus medialis (longus, obliquus) forces on the lateral contact force levels. In the model, the quadriceps were represented as five separate forces: vastus lateralis, vastus intermedius, rectus femoris, vastus medialis longus (VML), and obliquus (VMO). By varying the relative force generation ratios of the quadriceps heads, the patello-femoral contact forces were estimated. We sought to analytically determine the range of forces in the VMO and VML that cause a reduction or an increase of lateral contact forces, often the cause of patello-femoral pain. Our results indicated that increased contact forces are more dependent on combinations of muscle forces than solely VMO weakness. Moreover, our simulation data showed that the contact force levels are also highly dependent on the knee flexion angle. These findings suggest that training targeted to reduce contact forces through certain joint angles could actually result in a significant increase of the contact forces through other joint angles.

Patellofemoral stress during walking in persons with and without patellofemoral pain

OBJECTIVE

To determine whether individuals with patellofemoral pain (PFP) demonstrate elevated patellofemoral joint (PFJ) stress compared with pain-free controls during free and fast walking.

DESIGN

A cross-sectional study utilizing an experimental and a control group.

BACKGROUND

Although the cause of PFJ pathology is believed to be related to elevated joint stress (force per unit area), this hypothesis has not been adequately tested and causative mechanisms have not been clearly defined.

METHODS

Ten subjects with a diagnosis of PFP and 10 subjects without pain participated. All subjects completed two phases of data collection: 1) magnetic resonance imaging (MRI) assessment to determine PFJ contact area and 2) comprehensive gait analysis during self-selected free and fast walking velocities. Data obtained from both phases were required as input variables into a biomechanical model to quantify PFJ stress.

RESULTS

On the average, PFJ stress was significantly greater in subjects with PFP compared with control subjects during level walking. The observed increase in PFJ stress in the PFP group was attributed to a significant reduction in PFJ contact area, as the PFJ reaction forces were similar between groups.

CONCLUSION

Our results are consistent with the hypothesis that increased patellofemoral joint stress may be a predisposing factor with respect to development of PFP. Clinically, these findings indicate that treatments designed to increase the area of contact between the patella and the femur may be beneficial in reducing the PFJ stress during functional activities.

RELEVANCE

Patellofemoral pain affects about 25% of the population, yet its etiology is unknown. Knowledge of the biomechanical factors contributing to patellofemoral joint pain may improve treatment techniques and guide development of prevention strategies.

Dynamic magnetic resonance imaging

Among the devices helping with an accurate diagnosis, neither MRI nor arthroscopy is perfect; both delineate pathology in the knee joint with reasonable sensitivity and specificity. MRI, as a noninvasive and nonionizing modality, has made a significant contribution to the understanding of musculoskeletal disturbances. Static images through the patellofemoral joint in different degrees of flexion reveal only the degree of patellar tilt or subluxation, parameters that can be measured also on the axial view of conventional radiography. The accuracy of patellar position on static axial MRI is limited by the absence of muscle contraction, movement, and loading. Dynamic axial images of patellofemoral articulation can demonstrate the degree of flexion where patellar malalignment is maximal and assess whether or not it reduces. Arthroscopy, aside from its diagnostic values, provides the opportunity for treatment of intra-articular changes contributing to knee joint disturbances, but it is an invasive technique with potential risks of complications. The performed cost-effectiveness analysis of MRI is based mainly on estimation of intra-articular pathology of the acutely-injured knee [49,52,56]. There are scarce data on the cost-effectiveness of MRI of patellofemoral alignment in patellofemoral pain knees. Total examination time for active movement dynamic MRI procedure is approximately 8 to 10 minutes, thus it can be performed during routine MRI examination of the knee. In cases of suspected patellofemoral malalignment with symptoms that mimic other types of internal derangement of the knee joint, dynamic MRI can be a procedure of choice for detection of transient patellar dislocation, whereas a single clinical examination cannot differentiate from other internal knee pathologies. Dynamic MRI, although in an experimental phase, gives us a new perspective for dynamic study of the patellofemoral joint.

Q-angle influences tibiofemoral and patellofemoral kinematics

Numerous surgical procedures have been developed to correct patellar tracking and improve patellofemoral symptoms by altering the Q-angle (the angle between the quadriceps load vector and the patellar tendon load vector). The influence of the Q-angle on knee kinematics has yet to be specifically quantified, however. In vitro knee simulation was performed to relate the Q-angle to tibiofemoral and patellofemoral kinematics. Six cadaver knees were tested by applying simulated hamstrings, quadriceps and hip loads to induce knee flexion. The knees were tested with a normal alignment, after increasing the Q-angle and after decreasing the Q-angle. Increasing the Q-angle significantly shifted the patella laterally from 20 degrees to 60 degrees of knee flexion, tilted the patella medially from 20 degrees to 80 degrees of flexion, and rotated the patella medially from 20 degrees to 50 degrees of flexion. Decreasing the Q-angle significantly tilted the patella laterally at 20 degrees and from 50 degrees to 80 degrees of flexion, rotated the tibia externally from 30 degrees to 60 degrees of flexion, and increased the tibiofemoral varus orientation from 40 degrees to 90 degrees of flexion. The results show that an increase in the Q-angle could lead to lateral patellar dislocation or increased lateral patellofemoral contact pressures. A Q-angle decrease may not shift the patella medially, but could increase the medial tibiofemoral contact pressure by increasing the varus orientation.