- Research
- Open access
- Published:
Contact area and pressure changes of patellofemoral joint during stair ascent and stair descent
BMC Musculoskeletal Disorders volume 24, Article number: 767 (2023)
Abstract
Purpose
To investigate the differences of patellofemoral joint pressure and contact area between the process of stair ascent and stair descent.
Methods
The finite element models of 9 volunteers without disorders of knee (9 males) to estimate patellar cartilage pressure during the stair ascent and the stair descent. Simulations took into account cartilage morphology from magnetic resonance imaging, joint posture from weight-bearing magnetic resonance imaging, and ligament model. The three-dimension models of the patella, femur and tibia were developed with the medical image processing software, Mimics 11.1. The ligament was established by truss element of the non-linear FE solver. The equivalent gravity direction (-z direction) load was applied to the whole end of femur (femoral head) according to the body weight of the volunteers, and the force of patella was observed. A paired-samples t-test or Wilcoxon rank sum test to make comparisons between stair ascent and stair descent. Statistical analyses were performed using SPSS 22.0 using a P value of 0.05 to indicate significance.
Results
During the stair descent (knee flexion at 30°), the contact pressure of the patella was 2.59 ± 0.06Mpa. The contact pressure of femoral trochlea cartilage was 2.57 ± 0.06Mpa. During the stair ascent (knee flexion at 60°), the contact pressure with patellar cartilage was 2.82 ± 0.08Mpa. The contact pressure of the femoral trochlea cartilage was 3.03 ± 0.11Mpa. The contact area between patellar cartilage and femoral trochlea cartilage was 249.27 ± 1.35mm2 during the stair descent, which was less than 434.32 ± 1.70mm2 during the stair ascent. The area of high pressure was located in the lateral area of patella during stair descent and the area of high pressure was scattered during stair ascent.
Conclusion
There are small change in the cartilage contact pressure between stair ascent and stair descent, indicating that the joint adjusts the contact pressure by increasing the contact area.
Introduction
The principal functions of patella are to increase the moment arm of the quadriceps mechanisms and to transmit the tensile forces of quadriceps muscle to the patellar tendon. The stair ascent and stair descent are common activities of daily life that is biomechanically and physiologically more challenging than level walking. Analysis of the biomechanical and pathomechanics of knee requirements during the stair ascent and stair descent can add to our understanding of the diverse demands of this common activity in human motion. Comparing with walking, only a small number of studies have investigated stair ascent and stair descent of normal person [1,2,3,4].
The distribution of forces across the patellar articular surface during knee flexion involves the complex and dynamic interplay between soft tissue restraints and the bony geometry, so determination of the pressure distribution in vivo studies remains a challenging [5,6,7]. In order to solve this problem, the computational models of patellofemoral joint mechanics have been developed to understand patellofemoral function [8,9,10]. In particular, the finite element (FE) method offers the ability to predict spatial and temporal variations in stress, strain, and contact area/forces [11, 12]. Therefore, the finite element method is used to study the biomechanical characteristics of patellofemoral joint.
Earlier in vitro studies have demonstrated that, in weight bearing, contact pressures within the patellofemoral joint increase as the knee flexes from 0° to 90° [13, 14]. Comparison to descent stair, larger ranges of knee flexion angle and knee flexion moment are required during ascent stair. Reilly and Martens identified that the patellofemoral (PF)joint reaction force can be 0.8 times body weight during level walking [15]. The patellofemoral joint pressure was three times higher when climbing stairs than when walking horizontally [16]. The person extends the leg to the next step and contacts the stair at a higher flexion angle during ascent stair than that during stair descent [15, 17]. So, the patellofemoral joint contact area and pressure varies with the knee flexion angle in daily activities.
The goal of this study was to investigate the contact area and pressure changes of the patellofemoral joint during stair ascent and stair descent. It was assumed that the contact area and pressure of patellofemoral joint on the stair ascend than that on the stair descend.
Materials and methods
Subjects.
Experimental data were collected from 9 volunteers (9 males) without knee disorders or other neuromuscular disease that could affect the results of the experiment, and history of previous surgery. For experimental precision, these volunteers with height of 175 to 185Â cm and weight of 70 to 80Â kg were selected. Basic information was summarized in Table 1. For the accuracy of the experimental data, the gender, height and weight of the enrolled volunteers were restricted. Previous work estimating weightbearing patellofemoral joint contact area showed that males have 34% greater contact area compared with females [6]. All participants provided written informed consent and the experimental protocol was approved by the Ethics Committee of the Third Hospital of Hebei Medical University.
Two instantaneous actions were selected during the ascent stair and descent stair (knee flexion 60° on the stair ascent and 30° on the stair descent). (Fig. 1)There postures were chosen because they are positions at which peak knee extension moments are produced during stair ascent and stair descent [15].
The flow chart of obtaining experimental data and establishing the model was shown in Fig. 2. The geometrical data were obtained from computed tomography (CT) and magnetic resonance imaging (MRI) image acquisition from these healthy young male volunteers. Thin-slice CT (Spiral CT: Siemens 64-slice spiral CT, scan thickness 0.625 mm) was performed with 60° of knee flexion (upstairs) and 30° of knee flexion (downstairs), respectively. To segment the geometry of cartilage and the ligament of the knee, MR images of the knee was acquired with 3.0-T Siemens MR scanning system (Verio, Siemens, Erlangen, Germany) using fat-suppressed TIRM sequence (repetition time, 3500 ms; echo time, 32 ms; flip angle, 120°; matrix,320 × 224; field of view, 16 × 16 cm; slice thickness, 5 mm; scan time, 00:42 min). A 25% body weight resistance was provided through the pulley system to the footplate. Sagittal plane images were obtained at 30° and 60° of knee flexion.
Strategies for finite element models
The three-dimension models of the patella, femur and tibia were developed with the medical image processing software, Mimics 11.1 (Materialise, Inc., Belgium). The three-dimension models were shown in Fig. 3.
The CT and MRI scan data were transferred into the computer and the geometric model of the bone (femur, patella and tibia) the cartilage and meniscus were constructed by the Mimics11.1 software (Materialise, Inc., Belgium). Then the model was input into Pro/Engineer 4.0 software, and the advanced surface was established, meanwhile the geometry was corrected for establishing a complete knee joint model. To create a more perfect model, anatomical features were combined with digital models. The specific parameters are shown in Table 2.
The Meniscus was modeled with an elastic modulus of 59 MPa and a Poisson ratio of 0.475 [18, 19].The cartilage was modeled using eight-noded linear elastic solid elements with an elastic modulus of 15 MPa and a Poisson ratio of 0.3 [19, 20]. The cortical bone was modeled with an elastic modulus of 12,000 MPa and a Poisson ratio of 0.3 [21]. The Cancellous bone was modeled with an elastic modulus of 100 MPa and a Poisson ratio of 0.2 [21]. This linear elastic model is a reasonable assumption for dynamic activities such as a stair climb task. The patella ligament and quadriceps tendon were modeled using one-dimensional tension-only connector elements. The patella ligament was modeled distributed evenly across the attachment areas of the tibia and femur. The mesh models of knee with 30° and 60° flexion were established respectively by the ANSA software, utilizing the subdivided tetrahedral elements, truss elements (ligaments) and the data structure to fast search. The bulk element and the truss element of ligament are locally encrypted to ensure computational accuracy and speed. The tetrahedral elements were adopted C3D4 element and T3D2 element. As shown in Fig. 4
The ligament was established by truss element of the non-linear FE solver (ABAQUS Explicit, SIMULIA, Providence, RI). The nonlinear truss element was only strained but not compressed, which was consistent with the characteristics of knee ligament. Considering the lubrication and frictional force in the joint, a friction coefficient was modeled as 0.02 using surface-to-surface contact pairs.
In order to facilitate loading and reduce loading stress concentration in the model, reference points were established in the connection area between suprapatellar muscle and patella, the infrapatellar ligament and patella, and the infrapatellar ligament and tibia. The reference was shown as Fig. 5. According to the computational characteristics of the model, cortical bone, cancellous bone and cartilage were constrained by common nodes. The quadriceps femoris was simulated with a spring, which was used to simulate the interaction between the muscle and the bone.
Joint loading and muscle force model
Considering the gait of ascending and descending stair, this simulation only targeted a certain quasi-static angle in the ascending and descending process, so fixed constraints were applied to the bottom of the tibia. On account of the stair ascent and descent were carried by one leg, the equivalent gravity direction (-z direction) load was applied to the whole end of femur (femoral head) according to the body weight of the volunteers, and the force of patella was observed.
The length from the knee joint to the femoral head was measured by X-ray of the full length of the lower limb. Considering the torque caused by the length of femur and the position of the center of gravity. An equivalent gravity direction (-Z direction) load is applied to the entire femur end (femoral head) according to the body weight of the volunteer. The BEAM extension was used to determine the center of gravity position after measurement, and simulation was performed. (Fig. 6).
Data analysis and statistical methods
Statistical analyses were performed using SPSS 22.0 (IBM Co, Armonk, New York, USA) using a P value of 0.05 to indicate significance. The effects of variations were quantified in magnitude of the force and contact area of PF. Assessment of the pressure of patellofemoral cartilage focused on the maximum pressure of patellofemoral cartilage. Assessment of the contact pressure distribution focused on the maximum contact pressure. The data conform the normal distribution and were analyzed by paired t-test. When the data does not conform to the normal distribution, the non-parametric test (Wilcoxon rank sum test) is used for analysis. normal distributions Descriptive statistics were calculated as the Mean ± SD.
Results
The detailed results are in Table 3.
Magnitude of the force
During the stair descent (knee flexion at 30°), the pressure of patellar cartilage (Von Mises stress, the resultant force of normal stress and shear stress) of the patellar cartilage was 1.70 ± 0.04Mpa, the contact pressure of the patella cartilage was 2.59 ± 0.06Mpa, the pressure of femoral trochlea cartilage was 1.59 ± 0.03Mpa, the contact pressure of the patella cartilage was 2.57 ± 0.06Mpa.
During the stair ascent (knee flexion at 60°), the pressure of patellar cartilage was 1.79 ± 0.07 Mpa, the contact pressure of patellar cartilage was 2.82 ± 0.08Mpa. The pressure of femoral trochlea cartilage was 1.81 ± 0.07 Mpa, the contact pressure of femoral trochlea cartilage was 3.03 ± 0.11Mpa.
Characteristics of the contact area and position
The contact area between patellar cartilage and femoral trochlea cartilage was 249.27 ± 1.35mm2 during the stair descent, which was less than 434.32 ± 1.70mm2 during the stair ascent. Not only the total contact area but also the area of high stress was larger during the ascend stairs.
The contact position was located in the middle and lower part of the patella. The left side of this part was the high stress area during ascend stairs, demonstrating the force was relatively unbalanced. The contact position of patellofemoral joint was the middle and upper part of the patella during descent stairs, and the high stress area was scattered, which was relatively balance at this state. Pressure distribution characteristics of patellofemoral articular cartilage were shown in the Figs. 7  and 8.
Discussion
In this study, a finite element knee model with 3D was developed, it provided stresses and contact areas, thus allowing for more systematic biomechanical evaluation of the patellofemoral joint during stair ascent and stair descent. This study confirms that the patellofemoral joint contact area was greater during stair ascent than that during stair descent. The patellofemoral cartilage pressure increases but fluctuates less.
Wallace et al. [22] demonstrated that patellofemoral stress increased from 30° to 90°, peaking at 90° for muscle contractions. However, more recent study has shown that peak patellofemoral contact pressures are also observed at lower knee flexion angles [23,24,25]. It is proved, indicating the peak knee moment is reported to occur at about 60° during stair climbing [16]. Huang [26] hold the opinion that patellofemoral joint showed the highest peak contact pressure at 1.2 MPa at the initial flexion and the patellofemoral joint pressure increases slightly at the knee flexion reaches 60°. This result was also verified in this study, the patellofemoral pressure was greater during the stair ascent (60°) than that during stair descent(30°), it has increased 0.23Mpa.
The stair ascent is sometimes noted to be more demanding on the knee than descent stair. However, some scholars believe that the two activities are quite similar in terms of reaction load in the patellofemoral joint [15]. There is little difference in the contact pressure of patellofemoral cartilage between stair ascent and stair descent in this study. Why does the patellofemoral pressure increase with the knee flexion angle increase, but the maximum pressure difference of the patellofemoral joint increase less? With the knee flexes, the patella becomes engaged within the trochlear groove and contact area increases [27]. Using cadaver limbs and pressure sensitive film, Powers [25] reported a 68% increase in contact area between 15° and 60° knee flexion, whereas D’Agata [28]reported an 81% increase between 20° and 60°. Between knee flexion 20° and 60°, an 80% increase in total contact area was observed in the Gretchen’s research [29]. An 74% increase in total contact area was observed in our research with knee flexion from 30° to 60°. We believe that the reason for the small change in pressure may be the increase in the contact area as to the findings of other scholars [27, 30, 31]. The patellofemoral joint reduces the contact pressure of articular cartilage by increasing the contact area, which is a protective mechanism for the joint.
Akbarshahi [32] measured patellofemoral pressure in 4 healthy individuals during stair ascent and they found that the contact force and stress were greater on the lateral patellar facet compared with the medial facet. The peak patellofemoral contact pressure ranged from 3.7 to 6.1 times body weight on the lateral facet compared with 0.7 to 1.3 times body weight on the medial facet [32]. Other scholars have also identified that the distribution of pressure across the articular surface was remarkably uniform [13, 33]. In the knee flexion of 30°, with the patella engaged in the trochlea, the contact pressure on the lateral patellofemoral joint surface increased in this study. But, the medial–lateral patellofemoral joint compartment of the force during the 60° of knee flexion did not differ. Our results are similar to those of the cadaver specimens conducted by Huang et al. [26]. It can be considered that the change of pressure distribution is related to the characteristics of patellar tracking during knee flexion. Some studies revealed that the patella inclined medially first and then laterally during the extension to flexion [34]. Wilson [35] believed the patella inclined medially to -1.8° within the 0°-45° range, and inclined laterally to 2.5° within the 45°-90° range.
The discrepancies between these studies are most likely related to inherent differences in study design, such as detection methods or force simulation methods [2, 24,25,26]. There are some limitations in the current study. (1) The reconstruction and registration of bone models can induce errors. The errors could be minimized by improving the reconstruction and registration algorithm. (2) Only quadriceps loading was applied in the anatomical directions. (3) The effects of individual anatomical differences did not investigate. These limitations might have influenced our outcomes.
Conclusion
The patellofemoral joint contact area was greater during stair ascent (60° of knee flexion) than that during stair descent (30° of knee flexion). But there is a small change in patellofemoral cartilage contact pressure. It indicates that the joint adjusts the contact pressure by increasing the contact area.
Availability of data and materials
All data generated during this study are included in this published article and its supplementary information files.
References
Asay JL, Mündermann A, Andriacchi TP. Adaptive patterns of movement during stair climbing in patients with knee osteoarthritis. J Orthop Res. 2009;27:325–9.
Goudakos IG, König C, Schöttle PB, Taylor WR, Singh NB, Roberts I, Streitparth F, Duda GN, Heller MO. Stair climbing results in more challenging patellofemoral contact mechanics and kinematics than walking at early knee flexion under physiological-like quadriceps loading. J Biomech. 2009;42:2590–6.
Waiteman MC, Briani RV, Pazzinatto MF, Ferreira AS, Ferrari D, de Oliveira SD, de Azevedo FM. Relationship between knee abduction moment with patellofemoral joint reaction force, stress and self-reported pain during stair descent in women with patellofemoral pain. Clin Biomech (Bristol, Avon). 2018;59:110–6.
Thomeer LT, Guan S, Gray HA, Pandy MG. Articular contact motion at the knee during daily activities. J Orthop Res. 2022;40:1756–69.
Heino Brechter J, Powers CM. Patellofemoral stress during walking in persons with and without patellofemoral pain. Med Sci Sports Exerc. 2002;34:1582–93.
Besier TF, Draper CE, Gold GE, Beaupré GS, Delp SL. Patellofemoral joint contact area increases with knee flexion and weight-bearing. J Orthop Res. 2005;23:345–50.
Zhao D, Banks SA, D’Lima DD, Colwell CW Jr, Fregly BJ. In vivo medial and lateral tibial loads during dynamic and high flexion activities. J Orthop Res. 2007;25:593–602.
Wilharm A, Hurschler C, Dermitas T, Bohnsack M. Use of Tekscan K-scan sensors for retropatellar pressure measurement avoiding errors during implantation and the effects of shear forces on the measurement precision. Biomed Res Int. 2013;2013:829171.
Brimacombe JM, Wilson DR, Hodgson AJ, Ho KC, Anglin C. Effect of calibration method on Tekscan sensor accuracy. J Biomech Eng. 2009;131:034503.
Drewniak EI, Crisco JJ, Spenciner DB, Fleming BC. Accuracy of circular contact area measurements with thin-film pressure sensors. J Biomech. 2007;40:2569–72.
Li J, Lu Y, Miller SC, Jin Z, Hua X. Development of a finite element musculoskeletal model with the ability to predict contractions of three-dimensional muscles. J Biomech. 2019;94:230–4.
Besier TF, Pal S, Draper CE, Fredericson M, Gold GE, Delp SL, Beaupré GS. The role of cartilage stress in patellofemoral pain. Med Sci Sports Exerc. 2015;47:2416–22.
Huberti HH, Hayes WC. Contact pressures in chondromalacia patellae and the effects of capsular reconstructive procedures. J Orthop Res. 1988;6:499–508.
Hehne HJ. Biomechanics of the patellofemoral joint and its clinical relevance. Clin Orthop Relat Res. 1990;258:73–85.
Mason JJ, Leszko F, Johnson T, Komistek RD. Patellofemoral joint forces. J Biomech. 2008;41:2337–48.
Tanaka MJ, Voss A, Fulkerson JP. The anatomic midpoint of the attachment of the medial patellofemoral complex. J Bone Joint Surg Am. 2016;98:1199–205.
Taylor WR, Heller MO, Bergmann G, Duda GN. Tibio-femoral loading during human gait and stair climbing. J Orthop Res. 2004;22:625–32.
Harris MD, Anderson AE, Henak CR, Ellis BJ, Peters CL, Weiss JA. Finite element prediction of cartilage contact stresses in normal human hips. J Orthop Res. 2012;30:1133–9.
Kwon OR, Kang KT, Son J, Kwon SK, Jo SB, Suh DS, Choi YJ, Kim HJ, Koh YG. Biomechanical comparison of fixed- and mobile-bearing for unicomparmental knee arthroplasty using finite element analysis. J Orthop Res. 2014;32:338–45.
Kang KT, Kwon OR, Son J, Suh DS, Kwon SK, Koh YG. Effect of joint line preservation on mobile-type bearing unicompartmental knee arthroplasty: finite element analysis. Australas Phys Eng Sci Med. 2018;41:201–8.
Wang JP, Guo D, Wang SH, Yang YQ, Li G. Structural stability of a polyetheretherketone femoral component-A 3D finite element simulation. Clin Biomech (Bristol, Avon). 2019;70:153–7.
Wallace DA, Salem GJ, Salinas R, Powers CM. Patellofemoral joint kinetics while squatting with and without an external load. J Orthop Sports Phys Ther. 2002;32:141–8.
Csintalan RP, Schulz MM, Woo J, McMahon PJ, Lee TQ. Gender differences in patellofemoral joint biomechanics. Clin Orthop Relat Res. 2002;402:260–9.
Lee TQ, Yang BY, Sandusky MD, McMahon PJ. The effects of tibial rotation on the patellofemoral joint: assessment of the changes in in situ strain in the peripatellar retinaculum and the patellofemoral contact pressures and areas. J Rehabil Res Dev. 2001;38:463–9.
Powers CM, Lilley JC, Lee TQ. The effects of axial and multi-plane loading of the extensor mechanism on the patellofemoral joint. Clin Biomech (Bristol, Avon). 1998;13:616–24.
Huang W, Zeng X, Man GC, Yang L, Zhang Y. Simultaneous measurement of patellofemoral joint kinematics and contact mechanics in intact knees: a cadaveric study. Orthop Surg. 2022;14:2317–29.
Andriacchi TP, Yoder D, Conley A, Rosenberg A, Sum J, Galante JO. Patellofemoral design influences function following total knee arthroplasty. J Arthroplasty. 1997;12:243–9.
D’Agata SD. Pearsall AWt, Reider B, Draganich LF: An in vitro analysis of patellofemoral contact areas and pressures following procurement of the central one-third patellar tendon. Am J Sports Med. 1993;21:212–9.
Salsich GB, Ward SR, Terk MR, Powers CM. In vivo assessment of patellofemoral joint contact area in individuals who are pain free. Clin Orthop Relat Res. 2003;417:277–84.
O’Donnell P, Johnstone C, Watson M, McNally E, Ostlere S. Evaluation of patellar tracking in symptomatic and asymptomatic individuals by magnetic resonance imaging. Skeletal Radiol. 2005;34:130–5.
Frosch KH, Schmeling A. A new classification system of patellar instability and patellar maltracking. Arch Orthop Trauma Surg. 2016;136:485–97.
Akbarshahi M, Fernandez JW, Schache AG, Pandy MG. Subject-specific evaluation of patellofemoral joint biomechanics during functional activity. Med Eng Phys. 2014;36:1122–33.
Lee TQ, Anzel SH, Bennett KA, Pang D, Kim WC. The influence of fixed rotational deformities of the femur on the patellofemoral contact pressures in human cadaver knees. Clin Orthop Relat Res. 1994;302:69–74.
Xue Z, Pei Z, Zhang H, Tang C, Jia J, Zhang K, Zhang K, Guan Z. Development of an innovative measurement method for patellar tracking disorder. Aging (Albany NY). 2020;13:516–24.
Wilson NA, Press JM, Koh JL, Hendrix RW, Zhang LQ. In vivo noninvasive evaluation of abnormal patellar tracking during squatting in patients with patellofemoral pain. J Bone Joint Surg Am. 2009;91:558–66.
Acknowledgements
This work was supported by the Natural Science Foundation of Hebei Province [Grants number: H2021206505].
Funding
Natural Science Foundation of Hebei Province (Grants number: H2021206505).
Author information
Authors and Affiliations
Contributions
Xiaomeng Wang and Huixin Liu wrote the main manuscript text and Zhenyue Dong. Xiaobo Chen and Chenyue Xu collected the experimental data. Gang Ji and Huijun Kang prepared Figs. 1– 5. Wang Fei guided the experimental design. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All participants provided written informed consent and the experimental protocol was approved by the Ethics Committee of the Third Hospital of Hebei Medical University (S2021-032–1). All the methods were carried out in accordance with relevant guidelines and regulations.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, X., Liu, H., Dong, Z. et al. Contact area and pressure changes of patellofemoral joint during stair ascent and stair descent. BMC Musculoskelet Disord 24, 767 (2023). https://doi.org/10.1186/s12891-023-06882-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12891-023-06882-0