Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 10 October 2022
Sec. Biomechanics
Volume 10 - 2022 | https://doi.org/10.3389/fbioe.2022.957435

Biomechanical effects of medial meniscus radial tears on the knee joint during gait: A concurrent finite element musculoskeletal framework investigation

www.frontiersin.orgSentong Wang1* www.frontiersin.orgKazunori Hase2* www.frontiersin.orgShunsuke Kita3,4 www.frontiersin.orgShinya Ogaya5
  • 1Human-Mechanical System Laboratory, Graduate School of Systems Design, Tokyo Metropolitan University, Hachioji, Japan
  • 2Human-Mechanical System Laboratory, Faculty of Systems Design, Tokyo Metropolitan University, Hachioji, Japan
  • 3Biomechanics of Exercise and Sports in Physical Therapy Laboratory, Graduate Course of Health and Social Services, Saitama Prefectural University, Koshigaya, Japan
  • 4Department of Rehabilitation, Soka Orthopedics Internal Medicine, Saitama, Japan
  • 5Biomechanics of Exercise and Sports in Physical Therapy Laboratory, Department of Physical Therapy, Saitama Prefectural University, Koshigaya, Japan

The biomechanical variation in the knee during walking that accompanies medial meniscal radial tears stemming from knee osteoarthritis (OA) has not been explored. This study introduced a finite element musculoskeletal model using concurrent lower limb musculoskeletal dynamics and knee joint finite element analysis in a single framework and expanded the models to include knees with medial meniscal radial tears and total medial meniscectomy. The radial tears involved three locations: anterior horn, midbody, and posterior horn with grades of 33%, 50%, and 83% of the meniscus width. The shear and hoop stresses of the tear meniscus and tibial cartilage contact load, accompanying tears, and postmeniscectomy were evaluated during the stance phase of the gait cycle using the models. In the 83% width midbody tear group, shear stress at the end of the tear was significantly greater than in the intact meniscus and other tear groups, and the maximum shear stress was increased by 310% compared to the intact meniscus. A medial meniscus radial tear has a much smaller effect on the tibial cartilage load (even though in the 83% width tear, the cartilage/total load ratio increased by only 9%). However, the contact force on the tibial cartilage with total postmeniscectomy was increased by 178.93% compared with a healthy intact meniscus, and the peak contact pressure after meniscectomy increased from 11.94 to 12.45 MPa to 17.64 and 13.76 MPa, at the maximum weight acceptance and push-off, respectively. Our study shows that radial tears with larger medial meniscus widths are prone to high stress concentrations at the end of the tears, leading to the potential risk of complete meniscal rupture. Furthermore, although the tears did not change the cartilage load distribution, they disrupted the circumferential stress-transmitting function of the meniscus, thus greatly increasing the likelihood of the onset of knee OA. The significant increase in the tibial cartilage load with total postmeniscectomy indicates a potential risk of OA flare-ups. This study contributes to a better understanding of meniscal tear-induced OA biomechanical changes during human activities and offers some potential directions for surgical guidance of meniscectomies and the prophylaxis and treatment of OA.

Introduction

A meniscus tear is the most common damage associated with the meniscus, and it more commonly occurs on the medial compartment than on the lateral compartment (Fox et al., 2012). Moreover, radial tears are not uncommon in young patients (<50 years), often caused from trauma or degenerative processes (Terzidis et al., 2006; Howell et al., 2014), and can occur anywhere on either meniscus, anterior horn, posterior horn, or midbody (Fox et al., 2012). A meniscus tear can significantly affect the functions of the meniscus, such as shock absorption and load transmission, resulting in abnormal kinematics and redistribution of load on the knee (Kurosawa et al., 1980; Fox et al., 2015). The occurrence of these abnormal kinematics and mechanics often leads to increased cartilage wear or knee osteoarthritis (OA) changes (Felson, 2013). A large proportion of radial tears of the meniscus are considered irreparable, because they occur mostly in the white zone that lacks a blood supply (Scotti et al., 2013), so meniscectomy is often used to relieve pain and instability in the knee (Rao et al., 2015). However, meniscectomies have been reported to induce an increased risk of the progress of knee OA (Roemer et al., 2016). Hence, the biomechanical impact of meniscus radial tears and resultant meniscectomies on the knee joint needs to be investigated. Because activity-related radial tears are the most common (Terzidis et al., 2006), it is crucial to study the effects of radial tears on the kinematics and mechanics of the knee during human activities such as walking.

Previously, some researchers have used implantable pressure-sensitive film to study the contact mechanics between the torn menisci and articular cartilages in cadaveric knees in a knee machine that simulates walking (Bedi et al., 2010; Bedi et al., 2012; Gilbert et al., 2014). These works aimed to analyze the mechanical effect of the meniscus tears on the knee joint during human activities using an in vitro knee simulating apparatus. However, the in vivo knee physiological motion and loading pattern during activities were hard to replicate, and the contraction activity of muscles that act as “motors” for joint motion cannot be reproduced in cadavers. Moreover, the use of these sensors can severely alter meniscus translation and deformation, resulting in errors in predictions of the joint kinematics and contact mechanics (Hu et al., 2019). Thus, the kinematics and mechanics effects of meniscus tears on the knee during human activity are difficult to study using an in vitro approach.

Currently, a computational approach, the finite element (FE) method, has been proven to provide detailed geometric representations of torn menisci and can explain the stress alterations of menisci and load redistribution on the knee that accompanies meniscus tears (Kedgley et al., 2019; Li et al., 2019; Zhang et al., 2019). However, these studies cannot investigate the effect of meniscus tears on the knee during natural human activities, such as walking, because knee biomechanics were predicted based on simplified assumed loads and boundary conditions, such as fixed compressive load and flexion angles (Kedgley et al., 2019; Li et al., 2019; Zhang et al., 2019). In addition, muscles that produce the force that causes joint movement were not considered in biomechanical analysis of the knee (Mononen et al., 2013). Therefore, the biomechanical effect of the meniscus tears in the actual physiological conditions on the knee during human activities is still unclear.

This study aimed to elucidate 1) the stress alterations on the meniscus resulting from radial tears of the medial meniscus and 2) the effect of the tears and total medial meniscectomy on the knee biomechanics during the stance phase of the gait cycle. In the above discussion, we pointed out the problems of previous gait analyses using FE models: the inability to maintain a realistic representation of the knee joint tissue deformation and motion under natural physiological conditions. Therefore, we applied our previously developed finite element musculoskeletal (FEMS) model using a single concurrent framework combining the entire lower limb musculoskeletal dynamics and knee joint FE analysis, which has been proven to describe natural knee biomechanics during the gait (Wang et al., 2021; Wang et al., 2022), to clarify the biomechanical effect of tears on the knee during gait.

Methods

Subject experiments

A healthy male participant (age 29 years, height 175 cm, weight 80 kg) participated in the gait measurement. The participant was thoroughly informed about the purpose, methods, and caveats of the experiment. The experiment was approved by the research ethics committee of the Tokyo Metropolitan University. A straight, approximately 10-m-long walkway was prepared for the gait experiment. The participant was verbally instructed to “walk at a comfortable pace” with his preferred gait during the experiment. The marker-based motion trajectories (100-Hz sampling frequency, Vicon Nexus, Oxford Metrics Ltd., Oxford, United Kingdom) and ground reaction force (1000-Hz sampling frequency, TF-4060-D force plate, Tec Gihan Co., Ltd., Kyoto, Japan) were synchronously collected during a single stance phase.

Finite element musculoskeletal model

A FEMS model of the right lower limb, including a healthy knee with FE analysis previously developed in ABAQUS/Explicit (SIMULIA, Providence, RI, United States) (Figure 1), was used in the study. The modeling approach of the FEMS framework (Figure 2) of the lower limb has been discussed and validated previously (Wang et al., 2021; Wang et al., 2022) but is summarized below. The FEMS model included the geometry of all right lower limb bones and the articular cartilages and meniscus of the knee, which were segmented manually and reconstructed from computed tomography and magnetic resonance imaging scans (Figure 1). The femoral, tibial, and patellar bones were meshed using rigid triangular shell elements, and the cartilage and meniscus were defined using elastic eight-node hexahedral elements. A spherical joint with three degrees of freedom (DOF) was used to represent the hip, joints with six DOFs represented the tibiofemoral and patellofemoral joints, and a hinge joint at the ankle with one DOF were included in the FEMS model of the right lower limb, as shown in Figure 1. The anterior cruciate ligament, posterior cruciate ligament, medial collateral ligament, posteromedial capsule, lateral collateral ligament, anterolateral structure, and oblique popliteal ligament cross tibiofemoral joint were modeled as nonlinear spring bundles with a force-strain relationship (Abdel-Rahman and Hefzy, 1998):

fi={0εi0ki1(lili0)20<εi2εlki2[li(1+εl)li0]2εl<εi(1)

where fi is the force sustained by the ith ligament, li is the current length of the ith ligament, li0 is the slack length of the ith ligament, εl is the strain assumed to be constant at 0.03, and ki1 and ki2 are the stiffness coefficients of the spring elements representing the ith ligament for the nonlinear and the linear regions, respectively. The values of the material properties of the nonlinear spring elements are listed in Table 1 (Abdel-Rahman and Hefzy, 1998; Yu et al., 2001).

FIGURE 1
www.frontiersin.org

FIGURE 1. Right lower limb finite element musculoskeletal models containing knees with an intact healthy meniscus (A); a medial meniscus with radial tears in the anterior horn involving three grades with widths of (B) 33%, (C) 50%, and (D) 83%; midbody with widths of (E) 33%, (F) 50%, and (G) 83%; and posterior horns with widths of (H) 33%, (I) 50%, and (J) 83%. A model without medial meniscus, equivalent to total meniscectomy, was also prepared (not shown). The lower limb finite element musculoskeletal model includes a 12-DOF knee joint, along with 20 muscles, ligaments (ACL, anterior cruciate ligament; PCL, posterior cruciate ligament; MCL, medial collateral ligament; PMC, posteromedial capsule; LCL, lateral collateral ligament; ALS, anterolateral structure; OPL, oblique popliteal ligament; transverse ligament), cartilage, the meniscus, and meniscus horn attachments.

FIGURE 2
www.frontiersin.org

FIGURE 2. Overview of the dynamic gait simulation. Marker positions were input into the inverse kinematics to predict the joint primary angles (lumbar flexion–extension, hip joint 3-DOF rotations, knee joint flexion–extension, and ankle dorsiflexion–plantarflexion). The joint primary angles on each joint and ground reaction force measured from the force plate were input into the inverse dynamics to predict the muscle lengths, muscle moment arms, and right lower limb joint moments. The static optimization of the muscle force was performed with the joint moments, the muscle lengths, and the muscle moment arms as inputs to predict the muscle force. The knee contact mechanics were computed from the muscle forces and the ground reaction force to predict the knee joint secondary kinematics and contact mechanics. The predicted knee secondary kinematics were fed back to a new inverse dynamics analysis, and the muscle force optimization and knee contact mechanics were again computed.

TABLE 1
www.frontiersin.org

TABLE 1. Stiffness parameters of ligaments.

The FEMS model included 20 representative one-dimensional linear element muscles on the lower limb: the gluteus maximus (3 units), iliopsoas, pectineus, rectus femoris, vastus medialis, vastus intermedius, vastus lateralis, semimembranosus, semitendinosus, biceps femoris (short and long heads), gastrocnemius (medial and lateral heads), soleus, tibialis (anterior and posterior), extensor halluces longus, and extensor digitorum longus. The wrapping between the muscles and femoral and tibial bone was considered in the model to represent the muscle paths (Ali et al., 2016). The mechanical properties of the muscles were represented using a Hill-type dynamic model containing a contractile element (active fiber force-length property) in parallel with a passive elastic element (passive fiber force-length property) (Zajac, 1989).

Modeling of the meniscus with radial tears and meniscectomy

Modeling with FEMS, including knees with radial tears of the medial meniscus in the anterior horn, posterior horn, and midbody segment and without the medial meniscus (describing total meniscectomy), was an extension of the previously developed FEMS model with intact menisci. The medial meniscal model had radial tears in the anterior horn, midbody, and posterior horn involving three grades: 1) one-third width (33%), 2) one-half width (50%), and 3) five-sixths width (83%). These were created by partially removing the mesh from an intact medial meniscus (Figure 1). The intact and torn menisci were defined to be transversely isotropic with radial and axial moduli of 20 MPa and a circumferential modulus of 140 MPa (Donahue et al., 2002). The in-plane and out-of-plane Poisson’s ratios were 0.2 and 0.3, respectively (Donahue et al., 2002; Yao et al., 2006). The articular cartilage was considered as a linear elastic isotropic material with an elastic modulus of 5 MPa and a Poisson’s ratio of 0.46 (Li et al., 2001). It should be noted that to exclude the biomechanical effect of biological changes of tissue properties in the knee, only changes in the geometry of the meniscus tears were considered; in other words, the material property changes of tissues after tearing were not considered. The contact behaviors were defined with a coefficient of friction of 0.04 (McCann et al., 2009) for all contact pairs in the intact knee and knees with meniscus tear models: the femoral cartilage–medial and lateral meniscus, femoral cartilage–medial and lateral tibial cartilages, femoral cartilage–patellar cartilage, and between the meniscus and tibial cartilages on the medial and lateral sides. In the total meniscectomy model, the medial meniscus was resected, and the contact behavior of cartilage and the medial meniscus was not included.

The anterior transverse ligament connecting the anterior convex margin of the lateral meniscus to the anterior end of the medial meniscus was modeled as multiple linear spring elements with a stiffness of 12.5 N/mm (Donahue et al., 2003). The four meniscal horn attachments were assumed as multiple linear spring elements firmly connecting the meniscal horn faces and the tibial bone (Figure 1). The spring constants were calculated from Young’s modulus reported for the horn attachments (Hauch et al., 2009):

khp=Enhlhpah(2)

where h indicates the meniscal horn; p indicates the linear spring, where khp is the pth spring stiffness for the hth meniscal horn; lhp is the pth spring length of the hth meniscal horn; E is Young’s modulus of the meniscal horn; ah is the hth horn face area; and nh is the number of spring elements for the hth meniscal horn. The spring length of each spring element was calculated from the insertions of the spring at the node of the horn face and the node of the tibial attachments.

Concurrent finite element musculoskeletal framework for gait analysis

The gait analysis considering lower limb motion achieved through inverse kinematics, inverse dynamics, muscle force optimization, and FE analysis on the knee, using inertial measurement unit sensors or motion capture system has been described previously (Wang et al., 2021; Wang et al., 2021; Wang et al., 2022) but is briefly presented below. In the study, considering that measurement accuracy was more important than measurement convenience, the gait analysis was performed using a motion capture system and force plate, generally regarded as the gold standard for motion analysis, instead of inertial measurement sensor units. An inverse kinematics analysis was used to derive the primary angles (hip joint flexion–extension, internal–external rotation, and abduction–adduction angles; knee joint flexion-extension angle; and ankle joint dorsiflexion–plantarflexion angle) of joints from marker trajectories measured in a human gait experiment. These primary joint angles and the measured ground reaction force were input to the inverse dynamics analysis to determine the joint moments, muscle lengths, and muscle moment arms. The muscle lengths and muscle moment arms were input to a static optimization algorithm to estimate the muscle forces by satisfying an equilibrium equation of the joint moment, and the muscle forces were optimized by minimizing the sum of the cube of the muscular activation sm as follows (Hase and Yamazaki, 2002):

Mj(θj,es)=m[rjm(θj,es)×fjm(sm)](3)
C=minmsm3(4)

where j indicates the joint, m indicates the muscle, s indicates the knee secondary kinematics, Mj(θj,es) is the moment of the jth joint, θj is the primary angle of the jth joint, es is the secondary kinematics of the knee joint (note that the hip and ankle joints only include primary angles) for the sth degree of freedom, and rjm(θj,es) is the moment arm vector for the jth joint. With each muscular activation considered as a variable, the force of the mth muscle can be expressed as a function fjm(sm) of the activation about the jth joint. In addition, the muscular activations were constrained to be in the range of zero to one (inclusive).

The static optimization method was performed by implementing subroutines written in MATLAB (R2018b, MathWorks, MA, United States) using the interior-point method, which combines concurrent FE contact analysis on the knee joint in ABAQUS/Explicit at 16 evenly distributed time points (0%, 6.67%, 13.33%, 20%, 26.67%, 33.33%, 40%, 46.67%, 53.33%, 60%, 66.67%, 73.33%, 80%, 86.67%, 93.33%, and 100% of the stance phase) during one stance phase. The secondary kinematics of the knee, which were determined entirely from the interaction of the joint contact mechanics, muscle forces, and ligament restraint, were fed back to a new inverse dynamics analysis updating muscle lengths and muscle moment arms until the stability of optimization was ensured with the convergence criteria for equilibrium in the joint moment. The FEMS models containing an intact knee joint, damaged knee joints with meniscal tears, and a joint without a meniscus were used to perform the optimization procedures, separately, under the same input condition (marker position and ground reaction force).

Results

The peak shear stress, occurred on the white zone of the medial meniscus in the intact meniscus during the maximum weight acceptance of the stance phase (1st peak of the vertical ground reaction force) (Figure 3). The shear stress distribution of the meniscus changed as the anterior horn tear widened, shear stress in the anterior horn tear increased significantly, and the peak shear stress occurred at the end (deepest portion of tear near the side of the red zone) of the 83% tear. The changes in shear stress distribution that occurred on the medial meniscus with the midbody tear were not observed at widths of 33% and 50%. Tears of the posterior horn did not alter the shear stress distribution on the meniscus during the maximum weight acceptance. During the push-off at the stance phase (2nd peak of the vertical ground reaction force), changes in the shear stress on the medial meniscus were barely observed at three locations with tears.

FIGURE 3
www.frontiersin.org

FIGURE 3. Results of shear stress distribution on the medial meniscus involving a healthy knee and a knee with radial tears in the medial meniscus during the maximum weight acceptance and push-off.

High shear stress was observed at the ends of meniscal anterior horn tears, including all three grades compared to the intact meniscus, except during the terminal stance (50%–70% stance phase) (Figure 4A). In the 83% midbody tear group, the shear stress located at the end of the tear was significantly more prominent than that in the intact meniscus and other tear groups, and the maximum shear stress at the 90% stance phase was increased by 310% compared to the intact meniscus (Figure 4B). Alterations in the ends of the meniscal posterior horn tear groups were not detected during the stance phase (Figure 4C).

FIGURE 4
www.frontiersin.org

FIGURE 4. Shear stress differences on an intact meniscus and a meniscus with a medial radial tear at the ends of the tears. The left subfigures show the magnitude of stress at the (A) anterior horn, (B) midbody, and (C) posterior horn, involving an intact meniscus and 33%, 50%, and 83% width tears during the stance phase of the gait cycle. The right subfigures indicated the tear locations (orange box) at the anterior horn, midbody, and posterior horn and enlarged views of the tear locations.

The tibial cartilage contact forces were hardly affected by radial tears in the medial meniscus (Figures 5A–C). The ratio of the medial meniscus contact force to the whole medial contact force was 0.61 on average in the intact meniscus; that is, the tibial cartilage shared 39% of the total medial load, which illustrates the essential function of the medial meniscus in load sharing. An insignificant increase in the load shared by the tibial cartilage to the total medial side was observed with tears (Figures 5E,F), consistent with previous experimental in vitro results (Walker et al., 2015). However, after meniscectomy, the maximum contact force on the medial tibial cartilage was 3.6 times the bodyweight, compared to 2.17 times in the intact healthy knee. The contact pressure distributions on the medial tibial cartilages in the intact, torn, and meniscectomy knees are shown in Figure 6 at the maximum weight acceptance and push-off. The peak contact pressure on the tibial cartilage increased after meniscectomy to 17.64 MPa and 13.76 MPa, compared with the healthy knee (11.94 MPa and 12.45 MPa) at the maximum weight acceptance and push-off, respectively.

FIGURE 5
www.frontiersin.org

FIGURE 5. Contact forces on medial tibial cartilages of the intact knee (black solid curves) compared with total meniscectomy (black dashed curves) and (A) anterior horn, (B) midbody, and (C) posterior horn tears with widths of 33% (blue dashed curve), 50% (green dotted curve), and 83% (red dashed curve) during the stance phase of the gait cycle. The lower panels show the ratios of the medial meniscal contact forces to total knee joint contact forces on the medial side with tears in the (D) anterior horn, (E) midbody, and (F) posterior horn.

FIGURE 6
www.frontiersin.org

FIGURE 6. Results of the contact pressure distribution and resultant peak values (white words) occurring on the medial tibial cartilage obtained by a healthy intact knee, total meniscectomy knee, and radial tears during the maximum weight acceptance and push-off.

Discussion

In this study, models of the lower limb with a medial meniscus having radial tears and without a medial meniscus were developed based on a previously developed FEMS model (Wang et al., 2021; Wang et al., 2022), which has been proven to describe realistic knee biomechanics during the gait cycle. The predicted knee biomechanics results of the model containing an intact knee were compared with experimental in vivo and in vitro results (Kozanek et al., 2009; Bergmann et al., 2014; Gilbert et al., 2014; Clément et al., 2018) and proven to be reliable. We evaluated the knee biomechanical changes after meniscus tears and total meniscectomy using the models.

The anterior horn with intact meniscus bore stable shear stress with an average value of 12.27 MPa (Figure 4A) from the knee during the stance phase of the gait, compared with 6.72 MPa (Zhang et al., 2019) and 15.34 MPa (Zhang et al., 2019) predicted by other finite element simulation during slight flexion. However, some investigators also reported smaller results of less than 4 MPa (Li et al., 2019; Li et al., 2020); these differences might be due to the definition of the geometry, material properties, and contact behavior of the knee joint. In the medial meniscus, the mobility of the anterior horn was most prominent (Thompson et al., 1991; Vedi et al., 1999), and there was basically no obvious fluctuation in the shear stress of the anterior horn in the intact knee during the stance phase (Figure 4A), which may suggest that the anterior horn plays a role in stabilizing load transmission, but only during knee extension. As shown in Figure 7A, circumferential tensile stress occurred in the anterior horn of the intact meniscus, which showed a uniform decreasing trend from the inner circumference to the outer circumference. The shear stresses that accompanied the tear increased significantly at the ends of the anterior horn tears (Figure 4A); this was due to the destruction of the circumferential collagen fibers on the inner circumference, which results in the transfer of the load to the tear and produces more significant circumferential stress (Figure 7A). When the 83% width tear occurs, the changing trend of the circumferential tensile stress from the inner to the outer circumference changes from a uniform decrease to a gradual increase. Especially, the stress at the end of the tear increases extremely, and all the tear areas were stretched, which may lead to further tearing or complete rupture. On the intact meniscus, the average shear stress experienced by the midbody was 8.77 MPa, which is less than 12.27 MPa for the anterior horn and 11.54 MPa for the posterior horn. The anterior and posterior horns carry the anteroposterior shear forces applied to the femur, while the midbody primarily provides stability against medial subluxation by the midbody of the meniscus through hoop stress transmission (Hwang et al., 2012). The midbody of the intact meniscus exhibited a uniform hoop tensile stress trend that increases from the inner circumference to the outer circumference, as shown in Figure 7B. Therefore, the higher hoop stress from the outer circumference can make the meniscus firmly lock onto the medial femoral condyle (Figure 7A). The tears with widths of 33% and 50% hardly affected changes in the midbody shear stress distribution (Figures 3, 4B). However, substantial stress concentrations were detected at the end of the 83% midbody tear (Figure 4B), which may mean that longer tears at the midbody severely disrupt more circumferential fibers and hoop stress transmission, causing an increase in stress. Stresses reached the experimentally measured failure stress range of 20–170 MPa (Kohn et al., 1992). We observed that the hoop tensile stress in the tear region was replaced by compressive stress as the tear increased, resulting in the need for extremely high tensile stress at the tear to maintain circumferential stress transmission, as shown in Figure 7B. The observed compressive stress was favorable for compressing the tear surfaces together, resulting in a stable tear and possibly more favorable healing conditions. However, a study reported that five of the six radial tears had no evidence of healing and one had become longer (Weiss et al., 1989). Henning et al. (1988) suggested that only shallow radial tears do not require treatment because they heal spontaneously or remain asymptomatic. Our observations contradict these conclusions from clinical reports. This may be because the defined isotropic material models cannot accurately calculate the stresses within a fibrillar tissue such as the meniscus, and a fibril-reinforced material model, which can represent fluid pressure and site-specific collagen fiber orientation, should be used for accurate simulations (Klets et al., 2016). Even so, this finding could indicate that alongside vascularity, biomechanics may play an essential role in healing for nonoperative treatment or postoperative meniscal repair. The translation is significant in the anterior horn and small in the posterior horn of the medial meniscus. Therefore, stress readily accumulated at the junction of the middle and the posterior portions of the medial meniscus (Figure 3), consistent with the results reported in other studies (Bin et al., 2004). The posterior horn of the intact meniscus carried a high percentage of the load overall, with shear stress reaching a peak of 25.16 MPa at higher flexion angles during preswing, compared with 17.08 MPa for the anterior horn and 10.98 for the midbody during the same gait phase. The low translation of the posterior horn of the medial meniscus is a potential mechanism for meniscal tears, with a resultant “trapping” of the fibrocartilage between the femoral condyle and the tibial plateau during high flexion (Fox et al., 2012), possibly contributing to the high circumferential compressive stress distribution. As shown in Figure 7C, the tears did not change the hoop compressive stress distribution, and the stress that occurred at the tear surfaces caused them to be compressed together. Apart from vascularity, this may explain why radial tears in this region are most amenable to healing (Belzer and Cannon, 1993). In the addition, the slackness of the meniscus horn attachment also affects the loading of the meniscus horn. Guess et al. (2016) showed that increased meniscal extrusion occurs with laxity of the attachment, and increased laxity results in loss of meniscal function. In contrast, simulating an overly stiff meniscus horn attachment might lead to an overloading of the meniscal horns, and the study provides important guidance for related research on meniscus modeling.

FIGURE 7
www.frontiersin.org

FIGURE 7. Tensor plots of circumferential stress for the (A) anterior horn, (B) midbody, and (C) posterior horn with 83% width tears. The color changes from blue to red represent the stress variation from small (negative, compressive stress) to large (positive, tensile stress) on the tensor plot. Dashed arrows indicated circumferential stress in tensile or compressive directions. Thicker longer arrows represent higher magnitudes.

Bergmann group has conducted measurements of the in vivo tibiofemoral joint contact loads, which are 2.61 body weight (Kutzner et al., 2010) and 3.98 body weight (Dreyer et al., 2022), respectively, using instrumented implants. These studies directly measured the knee joint load of multiple patients and performed statistical analysis. Although these measurement results are not natural knee data, and due to the material properties of the total knee arthroplasty prosthesis, prosthesis geometry, and interprosthesis contact, the represented behavior varies widely from natural or removed meniscal knees. However, they can be used as a quantitative reference to properly verify the accuracy of our prediction results. The peak contact forces of 4.42 BW (Navacchia et al., 2019) and 4.02 BW (Hume et al., 2019) during the gait cycle were predicted by the Shelburne group using a finite element knee model without a meniscus. As far as we know, the joint contact force is generally accepted as being between two and 4.5 BW (Kumar et al., 2013; Bergmann et al., 2014; Richards et al., 2018; Hume et al., 2019; Shu et al., 2022). In general, the results predicted by computational models are usually greater than those predicted by in vivo experiments (Kutzner et al., 2010; Esrafilian et al., 2020; Wang et al., 2022). The role of the meniscus in sharing knee joint loading on the medial side was most significant during the loading response (0%–20% stance phase) and preswing (80%–100% stance phase) (Figure 5E,F). The ratio of the medial meniscus contact force to the whole medial contact force was 0.57 and 0.78 on average on the intact meniscus during the loading response and preswing, respectively. The trends of increased cartilage load caused by tearing were not apparent, compared with the ratio of meniscus–medial total contact force on the intact meniscus, which increased by 2.6%, 5.1%, and 7.7% on average in the anterior horn, midbody, and posterior horn tear groups, respectively, during the stance phase. Of these, even though the presence of an 83% width tear of the posterior horn during the preswing resulted in the most significant increase in cartilage load, the cartilage/total load ratio increased by only 9%. The findings suggested that medial meniscus radial tears have lesser influence on the tibial cartilage loading (Figure 5), consistent with previous experimental results (Bedi et al., 2010; Bedi et al., 2012), however, some studies have also concluded that meniscal tears lead to increased cartilage loading (Zhang et al., 2015; Sukopp et al., 2021). Previously, investigators have pointed out that meniscal tears in a healthy knee may ultimately lead to the occurrence of knee OA due to the loss of much of the meniscal functionality, such as circumferential stress transmission (Englund et al., 2009). According to results from previous case-control studies, OA is more likely to occur in the knee with meniscus tears without arthropathy than in knees with intact menisci, suggesting that meniscal tears occur before any visible cartilage changes (Englund et al., 2009). The cartilage loads were not significantly increased by the meniscal tears, but OA may begin to occur as tears result in loss of most of the meniscus function, which requires extreme vigilance. Therefore, the evaluation of meniscal function in addition to mechanical evaluation of cartilage is critical.

The peak contact pressure was 12.45 MPa on a healthy knee in our model. The results might be acceptable, compared with other computational studies, which reported results of 7.9 (Park et al., 2019), 13.0 (Halonen et al., 2017), and 15 MPa (Kłodowski et al., 2016). Bedi et al. (2010) reported maximum contact pressures of 7.4 ± 0.6 MPa and 6.4 ± 1.1 MPa (Bedi et al., 2012) in vitro studies on a cadaver knee joint with film contact pressure sensors to simulate gait conditions. Our calculated maximum contact pressures were high compared to these in vitro experiment results. However, with in vitro experiments, the contraction activity of muscles cannot be reproduced in cadavers, which may have resulted in the measurement of relatively small contact loads. Moreover, the use of these sensors can severely alter meniscus translation and deformation, possibly resulting in prediction errors of the contact mechanics. The maximum peak contact pressure on the tibial cartilage increased to 17.64 MPa after meniscectomy, which indicated a risk of cartilage damage (Kempson et al., 1971). The contact force on the tibial cartilage after total meniscectomy was increased by 178.93% compared with the intact meniscus knee (Figure 5). The magnitude of the peak contact pressure on the tibial cartilage with total meniscectomy greatly increased, and the contact distribution was significantly altered compared to the intact and tear meniscus groups. Previous research reported that radial tears extending to the periphery result in a loss of hoop tension, which has been described as functionally equivalent to a total meniscectomy (Allaire et al., 2008). Therefore, total meniscectomy is the mainstream treatment for long radial tears. However, despite the benefits of short-term pain relief, meniscectomy resulted in a substantial cartilage load that may increase the risk of OA progression; in other words, the effect on the knee biomechanics can be much worse than the meniscus tear itself. Our finding demonstrated that although tears can disrupt the meniscal circumferential stress transfer function, even wider meniscal radial tears have little effect on cartilage loading, whereas meniscectomy has a significant effect on the cartilage load. Therefore, we emphasize the importance of meniscal preservation whenever feasible, including consideration of meniscus repair or the implantation of an artificial meniscus.

This study had some limitations. First, the FEMS lower limb model and gait data were not from the same subject. The most accurate computational analysis of knee mechanics requires subject-specific model geometry and human motion data. Second, only one subject participated in the gait measurement experiment. Due to individual differences, the generalizability of our findings may be affected by random individual traits. However, using the same walking pattern of the same subject as input data while only changing the meniscus tear model eliminates the influence of the walking pattern itself on changes in joint load; the effects of meniscal tears can be directly evaluated. Therefore, we decided to perform the analysis using the data of one trial with one subject without considering the walking pattern and the subject’s own gait variation. In addition, analysis using the same gait measurement data was also beneficial for eliminating the influence of different gait patterns themselves in meniscal tear patients on changes in joint load. Third, the ligaments were simplified as one-dimensional elastic elements to reduce the computational expense. However, because wrapping and realistic soft tissue contact are not included in the model, the effect of the medial collateral ligament and oblique popliteal ligament contacting the medial meniscus could not be represented, which may have implications for meniscal mechanical prediction. Fourth, to exclude the influence of mesh conditions, we did not reconstruct and remesh the geometries of the menisci with tears; instead, these were created by partially removing the mesh from an intact medial meniscus. However, this modelling was unable to replicate a more realistic meniscus with tears geometry and may lead in errors in predictions of the joint contact mechanics. In addition, the division of the mesh was based on experience in the study; however, to identify the optimum mesh size, a mesh convergence analysis should be conducted. A final limitation was that the articular cartilages were defined as linear elastic isotropic material, whereas a biphasic fibril-reinforced material might better approximate the representation of cartilage’s dynamic response (Brindle et al., 2001). However, studies have also shown that the properties of linear elastic materials can accurately mimic the overall behavior of cartilage (Bell et al., 2009).

Conclusion

In this study, FEMS models of the lower limb, including knees with radial tears of the medial meniscus and total medial meniscectomy, were developed and used to investigate the effects of radial meniscal tears and total meniscectomy on the meniscal stress changes and biomechanical redistribution on the knee. Our research demonstrated that wider radial tears of the medial meniscus are prone to high-stress concentrations at the end of the tears, leading to the potential risk of developing complete meniscal ruptures. In addition, although the tears did not cause changes in cartilage load distribution, they disrupted the circumferential stress transmitting function of the meniscus; thus, knee OA may begin to develop. Significantly increased load on the tibial cartilage with the postmeniscectomy model indicated potential risks for the onset of OA. Therefore, surgical procedures such as meniscectomy should be applied conservatively; that is, excessive removal of meniscal tissue should be avoided. The modeling may provide a potential clinical tool for surgical decisions for patients with meniscus and other soft tissue injures, including more accurate pathology analysis and treatment of OA. Overall, the FEMS is expected to provide doctors with a reference for diagnosis.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by The research ethics committee of the Tokyo Metropolitan University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

This paper was conceptualized by SW, KH, SK, and SO. The methodology was developed by SW, KH, and SO. Subject motion data acquisition and analysis were carried out by SW and SK. Resources were arranged by KH and SO. The original draft was prepared by SW. All authors have read and agreed to the published version of the manuscript.

Funding

Supported in part by the Japan Society for the Promotion of Science [Grant-in-Aid for Scientific Research (B) 19H03970] (No. 50452199; Funder ID: 10.13039/501100001691).

Acknowledgments

We would like to express our gratitude to the Nagoya Orthopedic Surgery and Artificial Joint Clinic for providing medical images of patients.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Abdel-Rahman, E. M., and Hefzy, M. S. (1998). Three-dimensional dynamic behaviour of the human knee joint under impact loading. Med. Eng. Phys. 20, 276–290. doi:10.1016/S1350-4533(98)00010-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Ali, A. A., Shalhoub, S. S., Cyr, A. J., Fitzpatrick, C. K., Maletsky, L. P., Rullkoetter, P. J., et al. (2016). Validation of predicted patellofemoral mechanics in a finite element model of the healthy and cruciate-deficient knee. J. Biomech. 49, 302–309. doi:10.1016/j.jbiomech.2015.12.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Allaire, R., Muriuki, M., Gilbertson, L., and Harner, C. D. (2008). Biomechanical consequences of a tear of the posterior root of the medial meniscus. Similar to total meniscectomy. J. Bone Jt. Surgery-American Volume 90, 1922–1931. doi:10.2106/JBJS.G.00748

PubMed Abstract | CrossRef Full Text | Google Scholar

Bedi, A., Kelly, N. H., Baad, M., Fox, A. J. S., Brophy, R. H., Warren, R. F., et al. (2010). Dynamic contact mechanics of the medial meniscus as a function of radial tear, repair, and partial meniscectomy. J. Bone Jt. Surgery-American Volume 92, 1398–1408. doi:10.2106/JBJS.I.00539

PubMed Abstract | CrossRef Full Text | Google Scholar

Bedi, A., Kelly, N. H., Sc, B., Baad, M., Fox, A. J. S., Sc, M., et al. (2012). Dynamic contact mechanics of radial tears of the lateral meniscus: Implications for treatment. Arthrosc. J. Arthrosc. Relat. Surg. 28, 372–381. doi:10.1016/j.arthro.2011.08.287

PubMed Abstract | CrossRef Full Text | Google Scholar

Bell, J. S., Winlove, C. P., Smith, C. W., and Dehghani, H. (2009). Modeling the steady-state deformation of the solid phase of articular cartilage. Biomaterials 30, 6394–6401. doi:10.1016/j.biomaterials.2009.08.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Belzer, J. P., and Cannon, W. D. (1993). Meniscus tears: Treatment in the stable and unstable knee. J. Am. Acad. Orthop. Surg. 1, 41–47. doi:10.5435/00124635-199309000-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergmann, G., Bender, A., Graichen, F., Dymke, J., Rohlmann, A., Trepczynski, A., et al. (2014). Standardized loads acting in knee implants. PLoS ONE 9, 86035. doi:10.1371/journal.pone.0086035

PubMed Abstract | CrossRef Full Text | Google Scholar

Bin, S., Kim, J., and Shin, S. (2004). Radial tears of the posterior horn of the medial meniscus. Arthrosc. J. Arthrosc. Relat. Surg. 20, 373–378. doi:10.1016/j.arthro.2004.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Brindle, T., Nyland, J., and Johnson, D. L. (2001). The meniscus: Review of basic principles with application to surgery and rehabilitation. J. Athl. Train. 36, 160–169.

PubMed Abstract | Google Scholar

Clément, J., Toliopoulos, P., Hagemeister, N., Desmeules, F., Fuentes, A., and Vendittoli, P. A. (2018). Healthy 3D knee kinematics during gait: Differences between women and men, and correlation with X-ray alignment. Gait Posture 64, 198–204. doi:10.1016/j.gaitpost.2018.06.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Donahue, T. L. H., Hull, M. L., Rashid, M. M., and Jacobs, C. R. (2002). A finite element model of the human knee joint for the study of tibio-femoral contact. J. Biomech. Eng. 124, 273–280. doi:10.1115/1.1470171

PubMed Abstract | CrossRef Full Text | Google Scholar

Donahue, T. L. H., Hull, M. L., Rashid, M. M., and Jacobs, C. R. (2003). How the stiffness of meniscal attachments and meniscal material properties affect tibio-femoral contact pressure computed using a validated finite element model of the human knee joint. J. Biomech. 36, 19–34. doi:10.1016/s0021-9290(02)00305-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Dreyer, M. J., Trepczynski, A., Nasab, S. H. H., Kutzner, I., Schütz, P., Weisse, B., et al. (2022). European Society of Biomechanics S.M. Perren Award 2022: Standardized tibio-femoral implant loads and kinematics. J. Biomech. 141, 111171. doi:10.1016/j.jbiomech.2022.111171

PubMed Abstract | CrossRef Full Text | Google Scholar

Englund, M., Guermazi, A., and Lohmander, S. L. (2009). The role of the meniscus in knee osteoarthritis: A cause or consequence? Radiol. Clin. North Am. 47, 703–712. doi:10.1016/j.rcl.2009.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Esrafilian, A., Stenroth, L., Mononen, M. E., Tanska, P., Avela, J., and Korhonen, R. K. (2020). EMG-assisted muscle force driven finite element model of the knee joint with fibril-reinforced poroelastic cartilages and menisci. Sci. Rep. 10, 3026. doi:10.1038/s41598-020-59602-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Felson, D. T. (2013). Osteoarthritis as a disease of mechanics. Osteoarthr. Cartil. 21, 10–15. doi:10.1016/j.joca.2012.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox, A. J. S., Bedi, A., and Rodeo, S. A. (2012). The basic science of human knee menisci structure, composition, and function. Sports Health. 4, 340–351. doi:10.1177/1941738111429419

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox, A. J. S., Wanivenhaus, F., Burge, A. J., Warren, R. F., and Rodeo, S. A. (2015). The human meniscus: A review of anatomy, function, injury, and advances in treatment. Clin. Anat. 28, 269–287. doi:10.1002/ca.22456

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilbert, S., Chen, T., Hutchinson, I. D., Choi, D., Voigt, C., Voigt, C., et al. (2014). Dynamic contact mechanics on the tibial plateau of the human knee during activities of daily living. J. Biomech. 47, 2006–2012. doi:10.1016/j.jbiomech.2013.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Guess, T. M., Razu, S., Jahandar, H., and Stylianou, A. (2016). Predicted loading on the menisci during gait: The effect of horn laxity. J. Biomech. 48, 1490–1498. doi:10.1016/j.jbiomech.2015.01.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Halonen, K. S., Dzialo, C. M., Mannisi, M., Venäläinen, M. S., de Zee, M., and Andersen, M. S. (2017). Workflow assessing the effect of gait alterations on stresses in the medial tibial cartilage - combined musculoskeletal modelling and finite element analysis. Sci. Rep. 7, 17396. doi:10.1038/s41598-017-17228-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hase, K., and Yamazaki, N. (2002). Computer simulation study of human locomotion with a three-dimensional entire-body neuro-musculo-skeletal model (I. Acquisition of normal walking). JSME Int. J. Ser. C 45, 1040–1050. doi:10.1299/jsmec.45.1040

CrossRef Full Text | Google Scholar

Hauch, K. N., Oyen, M. L., Odegard, G. M., and Donahue, T. L. H. (2009). Nanoindentation of the insertional zones of human meniscal attachments into underlying bone. J. Mech. Behav. Biomed. Mat. 2, 339–347. doi:10.1016/j.jmbbm.2008.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Henning, C. E., Clark, J. R., Lynch, M. A., Stallbaumer, R., Yearout, K. M., and Vequist, S. W. (1988). Arthroscopic meniscus repair with a posterior incision. Instr. Course Lect. 37, 209–221.

PubMed Abstract | Google Scholar

Howell, R., Kumar, N. S., Patel, N., and Tom, J. (2014). Degenerative meniscus: Pathogenesis, diagnosis, and treatment options. World J. Orthop. 5, 597–602. doi:10.5312/wjo.v5.i5.597

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Xin, H., Chen, Z., Zhang, Q., Peng, Y., and Jin, Z. (2019). The role of menisci in knee contact mechanics and secondary kinematics during human walking. Clin. Biomech. (Bristol, Avon. 61, 58–63. doi:10.1016/j.clinbiomech.2018.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Hume, D. R., Navacchia, A., Rullkoetter, P. J., and Shelburne, K. B. (2019). A lower extremity model for muscle-driven simulation of activity using explicit finite element modeling. J. Biomech. 84, 153–160. doi:10.1016/j.jbiomech.2018.12.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, S. H., Jung, K. A., Lee, W. J., Yang, K. H., Lee, D. W., Carter, A., et al. (2012). Morphological changes of the lateral meniscus in end-stage lateral compartment osteoarthritis of the knee. Osteoarthr. Cartil. 20, 110–116. doi:10.1016/j.joca.2011.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kedgley, A. E., Saw, T., Segal, N. A., Hansen, U. N., Bull, A. M. J., and Masouros, S. D. (2019). Predicting meniscal tear stability across knee-joint flexion using finite-element analysis. Knee Surg. Sports Traumatol. Arthrosc. 27, 206–214. doi:10.1007/s00167-018-5090-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Kempson, G. E., Spivey, C. J., Swanson, S. A. V., and Freeman, M. A. R. (1971). Patterns of cartilage stiffness on normal and degenerate human femoral heads. J. Biomech. 4, 597–608. doi:10.1016/0021-9290(71)90049-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Klets, O., Mononen, M. E., Tanska, P., Nieminen, M. T., Korhonen, R. K., and Saarakkala, S. (2016). Comparison of different material models of articular cartilage in 3D computational modeling of the knee: Data from the Osteoarthritis Initiative (OAI). J. Biomech. 16, 3891–3900. doi:10.1016/j.jbiomech.2016.10.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Kłodowski, A., Mononen, M. E., Kulmala, J. P., Valkeapää, A., Korhonen, Rk, Avela, J., et al. (2016). Merge of motion analysis, multibody dynamics and finite element method for the subject-specific analysis of cartilage loading patterns during gait: Differences between rotation and moment-driven models of human knee joint. Multibody Syst. Dyn. 37, 271–290. doi:10.1007/s11044-015-9470-y

CrossRef Full Text | Google Scholar

Kohn, D., Wirth, C. J., Plitz, G. R. W., Maschek, H., Erhardt, W., Wülker, N., et al. (1992). Medial meniscus replacement by a tendon autograft. experiments in sheep. J. Bone Jt. Surg. Br. volume 74, 910–917. doi:10.1302/0301-620X.74B6.1447257

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozanek, M., Hosseini, A., Liu, F., Van de Velde, S. K., Gill, T. J., Rubash, H. E., et al. (2009). Tibiofemoral kinematics and condylar motion during the stance phase of gait. J. Biomech. 43, 1877–1884. doi:10.1016/j.jbiomech.2009.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, D., Manal, K. T., and Rudolph, K. S. (2013). Knee joint loading during gait in healthy controls and individuals with knee osteoarthritis. Osteoarthr. Cartil. 21, 298–305. doi:10.1016/j.joca.2012.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurosawa, H., Fukubayashi, T., and Nakajima, H. (1980). Load-bearing mode of the knee joint physical behavior of the knee joint with or without menisci. Clin. Orthop. Relat. Res. 149, 283–290. doi:10.1097/00003086-198006000-00039

PubMed Abstract | CrossRef Full Text | Google Scholar

Kutzner, I., Heinlein, B., Graichen, F., Bender, A., Rohlmann, A., Halder, A., et al. (2010). Loading of the knee joint during activities of daily living measured in vivo in five subjects. J. Biomech. 43, 2164–2173. doi:10.1016/j.jbiomech.2010.03.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Lopez, O., and Rubash, H. (2001). Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis. J. Biomech. Eng. 123, 341–346. doi:10.1115/1.1385841

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Yang, L., Zhang, K., Zhu, L., Wang, X., and Jiang, Q. (2020). Three-dimensional finite-element analysis of aggravating medial meniscus tears on knee osteoarthritis. J. Orthop. Transl. 20, 47–55. doi:10.1016/j.jot.2019.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Yang, X., Yang, L., Zhang, K., Shi, J., Zhu, L., et al. (2019). Biomechanical analysis of the effect of medial meniscus degenerative and traumatic lesions on the knee joint. Am. J. Transl. Res. 11, 542–556. eCollection 2019.

PubMed Abstract | Google Scholar

McCann, L., Ingham, E., Jin, Z., and Fisher, J. (2009). Influence of the meniscus on friction and degradation of cartilage in the natural knee joint. Osteoarthr. Cartil. 17, 995–1000. doi:10.1016/j.joca.2009.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Mononen, M. E., Jurvelin, J. S., and Korhonen, R. K. (2013). Effects of radial tears and partial meniscectomy of lateral meniscus on the knee joint mechanics during the stance phase of the gait cycle—A 3D finite element study. J. Orthop. Res. 31, 1208–1217. doi:10.1002/jor.22358

PubMed Abstract | CrossRef Full Text | Google Scholar

Navacchia, A., Hume, D. R., Rullkoetter, P. J., and Shelburne, K. B. (2019). A computationally efficient strategy to estimate muscle forces in a finite element musculoskeletal model of the lower limb. J. Biomech. 84, 94–102. doi:10.1016/j.jbiomech.2018.12.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S., Lee, S., Yoon, J., and Chae, S. (2019). Finite element analysis of knee and ankle joint during gait based on motion analysis. Med. Eng. Phys. 63, 33–41. doi:10.1016/j.medengphy.2018.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Rao, A. J., Erickson, B. J., Cvetanovich, G. L., Yanke, A. B., and Cole, B. J. (2015). The meniscus-deficient knee. Orthop. J. Sports Med. 3, 232596711561138. doi:10.1177/2325967115611386

CrossRef Full Text | Google Scholar

Richards, D. P., Barber, F. A., and Herbert, M. A. (2008). Meniscal tear biomechanics: Loads across meniscal tears in human cadaveric knees. Orthopedics 31, 347–350. doi:10.3928/01477447-20080401-32

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, R. E., Andersen, M. S., Harlaar, J., and van den Noort, J. C. (2018). Relationship between knee joint contact forces and external knee joint moments in patients with medial knee osteoarthritis: Effects of gait modifications. Osteoarthr. Cartil. 26, 1203–1214. doi:10.1016/j.joca.2018.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Roemer, F. W., Kwoh, C. K., Hannon, M. J., Hunter, D. J., Eckstein, F., Grago, J., et al. (2016). Partial meniscectomy is associated with increased risk of incident radiographic osteoarthritis and worsening cartilage damage in the following year. Eur. Radiol. 27, 404–413. doi:10.1007/s00330-016-4361-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Scotti, C., Hirschmann, M. T., Antinolfi, P., Martin, I., and Peretti, G. M. (2013). Meniscus repair and regeneration: Review on current methods and research potential. Eur. Cell. Mat. 26, 150–170. doi:10.22203/ecm.v026a11

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, L., Yamamoto, K., Yoshizaki, R., Yao, J., Takashi, S., and Sugita, N. (2022). Multiscale finite element musculoskeletal model for intact knee dynamics. Comput. Biol. Med. 141, 105023. doi:10.1016/j.compbiomed.2021.105023

PubMed Abstract | CrossRef Full Text | Google Scholar

Sukopp, M., Schall, F., Hacker, S. P., Ignatius, A., Dürselen, L., and Seitz, A. M. (2021). Influence of menisci on tibiofemoral contact mechanics in human knees: A systematic review. Front. Bioeng. Biotechnol. 9, 765596. doi:10.3389/fbioe.2021.765596

PubMed Abstract | CrossRef Full Text | Google Scholar

Terzidis, L. P., Christodoulou, A., Ploumis, A., Givissis, P., Natsis, K., and Koimtzis, M. (2006). Meniscal tear characteristics in young athletes with a stable knee: Arthroscopic evaluation. Am. J. Sports Med. 34, 1170–1175. doi:10.1177/0363546506287939

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, W. O., Thaete, F. L., Fu, F. H., and Dye, S. F. (1991). Tibial meniscal dynamics using three-dimensional reconstruction of magnetic resonance images. Am. J. Sports Med. 19, 210–216. doi:10.1177/036354659101900302

PubMed Abstract | CrossRef Full Text | Google Scholar

Vedi, V., Williams, A., Tennant, S. J., Spouse, E., Hunt, D. M., and Gedroyc, W. M. (1999). Meniscal movement. An in-vivo study using dynamic MRI. J. Bone Jt. Surg. Br. volume 81, 37–41. doi:10.1302/0301-620x.81b1.0810037

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, P. S., Arno, S., Bell, C., Salvadore, G., Borukhov, I., and Oh, C. (2015). Function of the medial meniscus in force transmission and stability. J. Biomech. 48, 1383–1388. doi:10.1016/j.jbiomech.2015.02.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Cai, Y., Hase, K., Uchida, K., Kondo, D., Saitou, T., et al. (2021). Estimation of knee joint angle during gait cycle using inertial measurement unit sensors: A method of sensor-to-clinical bone calibration on the lower limb skeletal model. J. Biomech. Sci. Eng. 17, 21-00196–00196. doi:10.1299/jbse.21-00196

CrossRef Full Text | Google Scholar

Wang, S., Hase, K., and Ota, S. (2022). A computationally efficient lower limb finite element musculoskeletal framework directly driven solely by inertial measurement unit sensors. J. Biomech. Eng. 144, 051011. doi:10.1115/1.4053211

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Hase, K., and Ota, S. (2021). Development of a lower limb finite element musculoskeletal gait simulation framework driven solely by inertial measurement unit sensors. Biomechanics 1, 293–306. doi:10.3390/biomechanics1030025

CrossRef Full Text | Google Scholar

Wiess, C. B., Lundberg, M., Hamberg, P., DeHaven, K. E., and Gillquist, J. (1989). Non-operative treatment of meniscal tears. J. Bone Jt. Surg. 71, 811–822. doi:10.2106/00004623-198971060-00003

CrossRef Full Text | Google Scholar

Yao, J., Snibbe, J., Maloney, M., and Lerner, A. L. (2006). Stresses and strains in the medial meniscus of an ACL deficient knee under anterior loading: A finite element analysis with image-based experimental validation. J. Biomech. Eng. 128, 135–141. doi:10.1115/1.2132373

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, C., Walker, P. S., and Dewar, M. E. (2001). The effect of design variables of condylar total knees on the joint forces in step climbing based on a computer model. J. Biomech. 34, 1011–1021. doi:10.1016/s0021-9290(01)00060-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zajac, F. E. (1989). Muscle and tendon: Properties, models, scaling, and application to biomechanics and motor control. Crit. Rev. Biomed. Eng. 17, 359–411.

PubMed Abstract | Google Scholar

Zhang, A. L., Miller, S. L., Coughlin, D. G., Lotz, J. C., and Feeley, B. T. (2015). Tibiofemoral contact pressures in radial tears of the meniscus treated with all-inside repair, inside-out repair and partial meniscectomy. Knee 22, 400–404. doi:10.1016/j.knee.2015.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, K., Li, L., Yang, L., Shi, J., Zhu, L., Liang, H., et al. (2019). Effect of degenerative and radial tears of the meniscus and resultant meniscectomy on the knee joint: A finite element analysis. J. Orthop. Transl. 18, 20–31. doi:10.1016/j.jot.2018.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, K., Li, L., Yang, L., Shi, J., Zhu, L., Liang, H., et al. (2019). The biomechanical changes of load distribution with longitudinal tears of meniscal horns on knee joint: A finite element analysis. J. Orthop. Surg. Res. 14, 237. doi:10.1186/s13018-019-1255-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: meniscus tear, finite element, musculoskeletal model, gait, knee, osteoarthritis

Citation: Wang S, Hase K, Kita S and Ogaya S (2022) Biomechanical effects of medial meniscus radial tears on the knee joint during gait: A concurrent finite element musculoskeletal framework investigation. Front. Bioeng. Biotechnol. 10:957435. doi: 10.3389/fbioe.2022.957435

Received: 31 May 2022; Accepted: 26 September 2022;
Published: 10 October 2022.

Edited by:

Jun Pan, Chongqing University, China

Reviewed by:

Tammy Haut Donahue, University of Memphis, United States
Ginu Unnikrishnan, BlueHalo, United States
Andreas Martin Seitz, Ulm University Medical Center, Germany

Copyright © 2022 Wang, Hase, Kita and Ogaya. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sentong Wang, sentongwang@gmail.com; Kazunori Hase, kazunori.hase@tmu.ac.jp

Download