- Research
- Open access
- Published:
Towards a validated musculoskeletal knee model to estimate tibiofemoral kinematics and ligament strains: comparison of different anterolateral augmentation procedures combined with isolated ACL reconstructions
BioMedical Engineering OnLine volume 22, Article number: 31 (2023)
Abstract
Background
Isolated ACL reconstructions (ACLR) demonstrate limitations in restoring native knee kinematics. This study investigates the knee mechanics of ACLR plus various anterolateral augmentations using a patient-specific musculoskeletal knee model.
Materials and methods
A patient-specific knee model was developed in OpenSim using contact surfaces and ligament details derived from MRI and CT data. The contact geometry and ligament parameters were varied until the predicted knee angles for intact and ACL-sectioned models were validated against cadaveric test data for that same specimen. Musculoskeletal models of the ACLR combined with various anterolateral augmentations were then simulated. Knee angles were compared between these reconstruction models to determine which technique best matched the intact kinematics. Also, ligament strains calculated by the validated knee model were compared to those of the OpenSim model driven by experimental data. The accuracy of the results was assessed by calculating the normalised RMS error (NRMSE); an NRMSE < 30% was considered acceptable.
Results
All rotations and translations predicted by the knee model were acceptable when compared to the cadaveric data (NRMSE < 30%), except for the anterior/posterior translation (NRMSE > 60%). Similar errors were observed between ACL strain results (NRMSE > 60%). Other ligament comparisons were acceptable. All ACLR plus anterolateral augmentation models restored kinematics toward the intact state, with ACLR plus anterolateral ligament reconstruction (ACLR + ALLR) achieving the best match and the greatest strain reduction in ACL, PCL, MCL, and DMCL.
Conclusion
The intact and ACL-sectioned models were validated against cadaveric experimental results for all rotations. It is acknowledged that the validation criteria are very lenient; further refinement is required for improved validation. The results indicate that anterolateral augmentation moves the kinematics closer to the intact knee state; combined ACLR and ALLR provide the best outcome for this specimen.
Introduction
Despite continued advances in ACLR procedures, the failure rate is reported as high as 14% [1, 2]. Residual rotational instability is the reported reason for up to 25% of these ACLR failures [3,4,5]. This anterolateral instability is caused by damage to the anterolateral ligament complex at the time of ACL injury and typically presents as anterolateral rotational laxity. Techniques such as modifying ACL graft tunnel positions [6, 7] or using an anatomic double-bundle ACLR [8, 9] are not always sufficient to control this anterolateral rotational laxity [10].
An alternate approach is to combine ACLR with an anterolateral extra-articular reconstruction [11,12,13,14,15,16,17]. This concept pre-dates modern intra-articular reconstruction techniques, and several anterolateral extra-articular procedures have been developed [18,19,20,21], which have been shown to reduce anterolateral rotational instability [15, 22]. However, determining which method controls anterolateral rotational laxity without potentially over-constraining the joint remains challenging [23,24,25].
Neri et al. [25, 26] in a recent cadaveric study compared several of these anterolateral techniques. They found that adding either anterolateral ligament reconstruction or modified Ellison procedures returned the reconstructed knee closest to the native knee kinematics. A return to natural knee kinematics is thought to provide an ideal mechanical environment for the ACL graft during its integration [27]. However, they also reported that some anterolateral procedures, while improving overall kinematics, also appear to cause increased lateral tibiofemoral compartment pressures [25, 26]. Evidence suggests that this may accelerate knee osteoarthritis [28].
This study aimed to complement the cadaveric studies [25, 26] by developing a subject-specific musculoskeletal computational model based on one of the cadaveric specimens. Kinematic data for the intact and ACL-sectioned cadaveric knee were used for validation. This model was then modified to computationally assess the impact of the different reconstructions applied during the cadaveric study on the knee kinematics and ligament strains.
We hypothesised that our computational model would: accurately capture the rotations and translations of the tibiofemoral joint for passive knee flexion; approximate the strain present in the knee ligaments throughout these movements; and correlate well with the results of the previous cadaveric study [26]. Replicating these passive knee outcomes should provide an important step toward developing a computational model which can be used as a clinical tool for planning and optimising patient-specific augmented ACLR procedures.
Results
Model validation
For both the intact and ACL-sectioned knees, a kinematics comparison of inverse kinematics (IK) and forward dynamics (FD) data are presented in Figs. 1, 2, and Table 1. IK and FD data reflect cadaveric [25, 26] and musculoskeletal model predicted kinematics, respectively. The FD internal rotation (IR) pattern of the intact and ACL-sectioned knee models was similar to the IK results (Figs. 1, 2). This was reflected in NRMSE values below 30% for most DoFs. The only exception was the anterior–posterior translation, which was around 75% and 60% for the intact and ACL-sectioned knees, respectively (Table 1).
The FD ligament strains in the intact, and ACL-sectioned knees have been compared with the IK strains calculated using the experimental kinematics [26] (Fig. 3; Table 2). The FD ligament strains were within 13.9–63.8% NRMSE of the respective IK values across the entire flexion range in the intact knee. Across ligaments, NRMSE values were less than 34.6% in strains for PCL, MCL, DMCL, LCL, and POPL (Table 2) and substantially greater for ACL, with 63.8% of NRMSE in the intact knee. Considerable differences in strain profile were seen for ACL strain values in the whole flexion range. However, similarities in the patterns for the remaining ligament strain curves are evident.
Evaluation of reconstructed knee model predictions
Following validation of the intact and ACL-sectioned knee models, the tibiofemoral kinematics and ligament strains of the reconstructed FD knee models were measured and compared against the experimental IK values.
Predicted kinematics
The FD rotational pattern of the tibia relative to the femur (internal rotation (IR) and external rotation (ER)) across the different reconstructed knee models were predicted (Fig. 4) during the 0–100° of passive knee flexion. A comparison with IK kinematics under 5Nm IR torque [25, 26] is also presented in Table 3. The FD IR results reflect the trends in the cadaveric experiment [26], with the ACLR + ALLR procedure restoring the overall IR kinematics through the full flexion range to the intact knee. However, the FD IR values showed a smaller range of variation than the experimental results [26]. The 5Nm ER [25, 26] torque caused the models to externally rotate by around 2.5–10°, compared to the literature range of 0–30° [29, 30], across 0–100° of knee flexion.
Predicted ligament strains
The predicted FD ligament strains for the knee model with 5Nm IR applied torque [25, 26] are presented in Fig. 5. Compared with the intact model, the ACL-sectioned model showed a uniform increase in strain values across the PCL, MCL, POPL and DMCL. In contrast, the LCL experienced a uniform strain decrease in the ACL-sectioned model. The ACLR + ALLR produced the greatest overall reduction in ACL graft strain throughout the flexion cycle (Fig. 5A). ACLR + ALLR, ACLR + Mac, and ACLR + DL showed reduced ACL, MCL, DMCL, LCL, and POPL strains throughout the full flexion cycle compared to strains following isolated ACLR. The ACLR + Ell technique resulted in smaller strain values in MCL and DMCL throughout the full flexion cycle compared to isolated ACLR. In contrast, the strain patterns of ACL, PCL, LCL, and POPL changed slightly among the ACLR + Ell technique.
NRMS errors between the IK and FD strains within the ACL, PCL, MCL, DMCL, LCL and POPL bundles with 5 Nm IR applied torque [25, 26] have also been calculated and reported in Table 2. The maximum NRMS differences between the IK and FD strains within the ACL and PCL were seen in the intact knee with 63.8% and 34.6%, respectively. In comparison, the MCL (with NRMSE of 20.1%), and the DMCL (with NRMSE of 13.9%), showed the lowest NRMSE in the intact knee compared to the other knee states (Table 2). LCL and POPL had lower NRMSE in the intact (with NRMSE of 33.6% and 22.1%, respectively) and ACL-sectioned (with NRMSE of 33.1% and 20.7%, respectively) than the rest of the knee states. The NRMSEs for PCL, MCL, DMCL, LCL, and POPL ligament strains were on the order of 20–30%, while that relating to ACL was around 40–60% (Table 3).
Discussion
This musculoskeletal modelling study developed an OpenSim model of the intact knee consisting of subject-specific tibiofemoral contact surfaces and 19 ligament bundles. The predicted FD model kinematics were validated against experimental IK results for the matching cadaveric knee in intact and ACL injured states. Following validation, the reconstructed knee model was created by modifying the intact model to reflect the reconstruction grafts. FD was used to calculate the knee kinematics and ligament strains, which were compared to the cadaveric study data to determine which procedure best restored the knee to intact kinematics. The results indicate that anterolateral augmentation moves the kinematics closer to the intact knee state; combined ACLR and ALLR provides the best outcome for this specimen.
Model validation
The musculoskeletal model was validated for the intact and ACL injured states by comparing the predicted FD kinematics of these two models to corresponding experimental data. The ACL injured model was based on the intact model, with the ACL and ALL removed. It is acknowledged that the criterion for validation is very lenient (NRMSE < 30%). This was due to the ad hoc method of modifying the model to try to match the experimental kinematics.
The primary focus of this study was rotations, as these are typically the focus of both clinical and computational papers. However, tibiofemoral translations were also reported. FD rotational knee kinematics matched with the cadaveric experimental results for both intact and ACL-sectioned models (NRMSE < 30%) [31, 32]. The ranges of motion were also within the physiological range of motion [33,34,35]. Differences between the FD-predicted results and the cadaveric experimental data were more obvious in the translation results. However, the only anterior–posterior translation did not meet the acceptance criteria, with NRMSE > 30%. Differences between the FD and IK data are most likely attributed to the model's contact surfaces, ligament slack lengths, and insertion points, which should be further refined and optimised to produce a closer match [36,37,38,39,40].
The predicted ligament strains for the intact and ACL-sectioned models matched reasonably well with the experimental results calculated using IK in OpenSim [26]. The average NRMSE was less than 26%, except for the ACL strain in the intact case (NRMSE = 63.8%). The errors that can be seen highlight the results' sensitivity to slight variations in contact geometry and ligament placement.
The predicted ligament strain values agreed with other musculoskeletal modelling studies [41, 42], suggesting that applying a 5Nm IR torque [25, 26] increases MCL and DMCL strain at full extension and increases PCL, LCL, and POPL bundle strain in flexion. The published ACL strains are lower than those predicted by our FD model. The considerable difference in ACL strain values between the IK and FD results suggests that the model might not accurately predict the ACL strains. The authors hypothesise that perhaps the central location of the ACL makes it more sensitive to changes in contact geometry and the exclusion of other soft tissues, such as the menisci.
Finally, it is also worth noting that the maximum ligament strain values achieved during the passive flexion motion of the knee model of this study were 10% (Fig. 3) which is in close agreement with the limit proposed by Blankevoort et al. [43], adding to the validation of the defined ligament slack lengths.
Evaluation of reconstructed knee model predictions
Predicted kinematics: The predicted IR pattern of the reconstructed knee models generally reflected the cadaveric IR values [26]. The isolated ACLR model did not restore intact knee kinematics, consistent with the experimental results [26], supporting previous findings that the isolated ACLR can lead to residual rotational instability [4, 5].
Among the various anterolateral procedures investigated, this modelling study revealed that the ACLR + ALLR better restored the overall IR kinematics through the full flexion range toward the intact state, which correlated well with the experimental studies [25, 26, 44]. Neri et al. [25, 26], in a biomechanical study of 10 cadaveric knees that underwent ACLR combined with various lateral reconstructions, demonstrated that the ACLR + ALLR procedure provided additional rotational control whilst protecting the ACL graft without risking over-constraint of the joint. Sahanand et al. [44], in a clinical study of 25 patients who underwent ACLR + ALLR with an average of 31.5 months follow-up, reported a significant improvement in the patient outcomes compared to pre-surgery, with no graft failure and residual instability.
This modelling study showed that the ACLR + Ell procedure also provided rotational control whilst protecting the ACL graft without risking over-constraint for 0–40˚ of flexion cycle, compared to the isolated ACLR knee (Fig. 4). Our results confirm the findings of Neri et al. [26], who reported that the ACLR + Ell procedure could restore physiological kinematics. Also, our results support Devitt et al. [45], who suggested that ACLR + Ell can reduce anterolateral instability in the anterolateral capsule-injured knee and restore kinematics close to the intact state. While the ACLR + Mac and ACLR + DL models also provided rotational control for IR (Fig. 4), they over-constrained the knee kinematics throughout the flexion range (Fig. 5). These findings are supported by Neri et al. [25, 26], whose cadaveric results indicated the possibility of over-constraining the joint using these procedures. Also, Geeslin et al. [46] reported that the combination of ACLR + DL procedure resulted in both over-constraint of the joint and significant reductions in tibiofemoral motion at most knee flexion angles.
As defined by applying the 5 Nm internal and external torques [25, 26], the limits of tibial rotation are shown in Fig. 4 as a flexion function for all the knee models. These limits determine the freedom-of-motion range of the developed knee model, from full extension up to 100° of flexion, and are considered to comprise the rotation envelope of passive knee motion. Overall, the rotation envelope of passive knee flexion (IR and ER) presented in this developed model (Fig. 4) compares well with our experimental data (Table 3) [26] and similar literature [29, 30, 47,48,49] further supporting the validity of our knee models.
Predicted ligament strains: This study estimated the strains in the major knee ligaments based on musculoskeletal modelling techniques (Fig. 5). These findings indicate that ACLR combined with anterolateral procedures reduced the strain levels within the ACL graft compared to the isolated ACLR technique. Specifically, the ACL + ALLR technique had the greatest reduction in ACL graft strains, which supports evidence that this technique improves knee kinematics immediately post-surgery and may improve patient outcomes in the longer term [26, 50]. Specifically, all of the anterolateral augmented ACLR techniques saw a reduction of strain in the ACL graft to levels typical of a healthy knee, which is likely to be beneficial in protecting the ACL graft during healing, and in doing so, leading to improved outcomes [49].
Figure 4 shows that ACLR + DL and ACLR + Mac result in over-constrained internal tibia rotation compared to the other lateral extra-articular methods. In all the models, the axial rotation was restrained mainly by the MCL and DMCL bundles (with a 0–10.5% strain range) and less so with LCL and POPL (Fig. 5 and Table 2), which also agrees with the literature on ligament function [51,52,53].
This study also showed that as the ACL-sectioned knee model was passively flexed with 5Nm of IR torque [25, 26], LCL experienced a uniform reduction in strain values than the intact knee (Fig. 5E). These results are different from current thinking on LCL ligament elongation behaviour reported in the literature [54,55,56] and can likely be attributed to differences in the location of ligament attachment sites or ligament slack length values in the model.
Limitations
The main limitation of this model is related to the manual, ad hoc method of generating the contact surfaces. This method removed the possibility of achieving a close match between the cadaveric and predicted results and led to the adoption of a more lenient validation criteria for the musculoskeletal model. Further, the manual method of extracting the cartilage geometries means these results are not easily reproducible. Future work should develop a reproducible computational process for generating these contact surfaces and driving them towards an optimum which fulfils more stringent validation criteria. Developing an appropriate computational optimisation program was beyond the scope of this study. However, the agreement achieved between experimental and predicted results using our methodology highlights the potential for achieving a fully validated model in future studies. In addition, while the resulting surfaces created a reasonable approximation of the patient geometry, it must be acknowledged that critical soft tissue, such as the menisci, was not included. While there is precedent for this in several seminal papers on musculoskeletal modelling [36, 57,58,59], this omission has likely contributed to some of the higher NRMSE values found in this study. The menisci play an important role in stabilising the joint, providing cushioning, and restricting excessive movement by ensuring the femoral condyles are guided through motion [60, 61]. Also, it is noted that the torque applied during the experiment [26] was not continuously monitored. As the model applied a constant torque, any deviations from this during the experiment would lead to differences in the kinematics and strain outputs. Finally, the differences observed between our predicted FD knee kinematics and the experimental results [26] may be due to this study using a subject-specific knee model based on only one of the ten cadaveric models used in the experiment.
Conclusion
To our knowledge, this is the first study to provide a direct comparison between knee joint kinematics and ligament strains predicted by an OpenSim FD model and those measured experimentally for ACL reconstruction with different lateral augmentation procedures. The ability to approach validation of our musculoskeletal knee model by correlating it against an accepted cadaveric model is a crucial step in creating a non-invasive computational tool for assessing and optimising ACLR procedures for individual patients. A comparison of various currently used anterolateral augmentation procedures revealed that adding those procedures to the isolated ACLR restored intact knee kinematics in an ACL-sectioned knee. In comparison, combined ACLR and ALLR provide the best outcome for this specimen.
Materials and methods
Experimental data
This study utilised data from a previous controlled laboratory experiment that assessed the knee kinematics and lateral compartment pressures of ten lower extremities cadaveric specimens [25, 26]. For all knee states, knee kinematics and contact pressures were acquired for three cycles of passive knee flexion/extension combined with 5 Nm of internal rotation (IR) torque [25, 26] applied manually to the tibia. The test was repeated with 5Â Nm of external rotation (ER) [25, 26]. To track the knee kinematics, the limb was prepared with four bi-cortical pins (two for the femur and two for the tibia) with two retro-reflective markers fixed to each pin using a validated protocol [62]. Knee kinematics were acquired using a 3D optoelectronic motion capturing system consisting of five Bonita cameras on tripods (Vicon, LA, USA). A test with neutral rotation was also conducted, but these results were excluded from this study because the lack of applied torque led to an excessive amount of variation between flexion cycles due to knee laxity.
Knee model development
An intact, subject-specific musculoskeletal knee model consisting of a six-DoF tibiofemoral joint, a three-DoF patellofemoral joint, ligaments and capsular bundles was created in OpenSim (R3.3, Simtk, USA) [63] (Fig. 6). The intact knee model was initially adapted from Xu et al.'s [64] and Schmitz et al.'s [41] models and then further developed by incorporating subject-specific contact surfaces (Fig. 6).
Three-dimensional specimen-specific bone and cartilage geometries were reconstructed through the segmentation of CT and MR images. The contact surfaces were further amended using Fusion 360 software (R2020, Autodesk, USA) to guarantee medial and lateral contact over the entire natural knee rotation and match the experimentally measured kinematics for the intact and ACL-sectioned models [26]. Joint centres and bony landmarks, along with the precise location of drilled holes in the femur and tibia for reconstruction surgery, were also determined from segmented CT images (see Sects. 1.1–1.3 in Additional file 1).
Nineteen ligament bundles were included in the intact knee model (Table 4), with attachment points obtained from the CT and MRI images and the literature [65,66,67,68,69,70,71,72,73,74] (Fig. 6; Table 4). Two bundles represented the anterior cruciate ligament (ACL): anterior and posterior bundles (aACL, pACL). The posterior cruciate ligament (PCL) was represented by two bundles; anterior and posterior bundles of the posterior cruciate ligament (aPCL, pPCL). The medial collateral ligament (MCL) was divided into two groups of ligament bundles: superficial and deep. The superficial layer of MCL was further divided into three bundles: anterior, inferior, and posterior (aMCL, iMCL, pMCL). The deep layer of MCL was subdivided into anterior and posterior bundles (aDMCL, pDMCL). The lateral collateral ligament (LCL), the popliteofibular ligament (POPL), and the anterolateral ligament (ALL), were each represented with one bundle. The posterior capsule was represented by four bundles: the anterior, oblique, medial, and lateral bundles (CAPa, CAPo, CAPm, CAPl). The patellar tendon was modelled as three bundles: the central, medial, and lateral bundles (cPT, mPT, lPT).
Ligament bundles were modelled as non-linear elastic springs with linear damping [47, 70, 75,76,77], with predefined stiffness and slack length values (Table 4). To avoid penetration of the ligament bundles into the bones, wrapping surfaces were also included in the model [41] (Sect. 1.4 in Additional file 1).
Knee models
The intact OpenSim model was developed first using the process outlined above. The ACL-sectioned model was then created by removing the ACL ligament bundles from the intact model. Matching the kinematics of these two models to the corresponding cadaveric kinematic results provided validation for the OpenSim modelling process [26]. Following validation, new models (Fig. 7) were created for each reconstruction [26] by modifying the validated intact knee model, which was always used as the base model for consistency [26]. Seven of the eight knee groups explored in the Neri et al.'s cadaveric experiment [26] were modelled in this study (Fig. 7). Superficial Lemaire was excluded from the modelling as it required significant changes to the wrapping surface definitions of the validated base model used for all reconstruction models. Below are detailed descriptions of the seven modelled knee groups:
-
1)
The intact knee.
-
2)
The ACL-sectioned knee, in which the anterolateral ligament (ALL) was also sectioned to match the cadaveric experiment [26].
-
3)
ACLR, in which the ACL graft was modelled by one single ligament bundle fixed with the knee at 30° and with an 80 N tension [26]. The properties of a quadrupled hamstring autograft (semi-tendinosis) [78, 79] were applied to the graft as this was the graft used in the experiment [26].
-
4)
ACLR combined with the ALL-reconstructed knee (ACLR + ALLR) (Fig. 7A), in which the ALL graft was modelled as a single bundle ligament with properties similar to a typical gracilis graft [78, 79] passing under the ITB and over the LCL [80]. The tibial insertion point of the ALL graft was located equidistant between the centre of Gerdy's tubercle (GT) and the anterior margin of the fibular head and 10 mm distal to the joint line [81, 82]. The femoral insertion point was located 5 mm proximal and posterior to the LCL's femoral insertion [81,82,83].
-
5)
ACLR + DL was similar to the ACLR + ALLR, except that the tibial tunnel of the anterolateral reconstruction was located at the centre of Gerdy's tubercle (GT) (Fig. 7B).
-
6)
ACLR + Mac, in which the graft extended from the same tibial insertion point (GT) to the femoral insertion point located 70 mm proximal to the femoral epicondyle (Fig. 7C).
-
7)
ACLR + Ell technique, in which the graft was modelled as a strip of the iliotibial band (ITB) detached distally, passed underneath the LCL, and finally attached to the same tibial insertion site (GT) (Fig. 7D).
In all of the reconstructed models, the lateral graft was always fixed in the same condition: in neutral rotation at 30° of flexion and with 20 Nm of applied tension characterised by a non-linear tension bundle with an average stiffness as reported in the literature [78, 79].
Musculoskeletal simulation
The simulations were performed via OpenSim (R3.3, Simtk, USA) [63] using two main toolboxes: inverse kinematics (IK) to compute the knee joint kinematics and ligament strains from the cadaveric experiments; and forward dynamics (FD) to predict the knee kinematics and the ligament strains within the modelled knees. The pelvis and femur were fixed in all directions during the entire simulation for all simulations. For the remainder of the article, IK results will be referred to as experimental results and FD as predicted results.
For IK, the experimental marker trajectories were used to drive the tibiofemoral joint. For FD, only the flexion component of the tibia’s kinematics was prescribed, with other rotations and translations restricted only by the ligaments and contact surfaces. The FD analysis tool was used to predict the knee joint behaviour throughout the entire 100° flexion cycle while applying 5N of IR/ER torque to the tibia for all knee models [25, 26].
The FD kinematics for the intact and ACL-sectioned knee states were compared to their corresponding IK kinematics [25] to validate the OpenSim model. The kinematics measurements considered in the validation were: (1) transverse plane motion, including internal–external rotation; (2) sagittal plane motions, including anterior–posterior and proximal–distal translations; and (3) frontal plane motions, including lateral–medial translation and adduction–abduction rotation. In all three planes, the kinematics were defined as the motion of the tibia relative to the femur.
Once the intact and ACL-sectioned knee models were validated, then FD was used to predict the kinematics for all the reconstructed knee models. IR and ER torques of 5 Nm [25, 26] were applied to all the models during 0–100° of passive knee flexion. These predicted kinematics were compared to the average kinematics reported by Neri et al. [26].
The ligament strains were also compared for the cadaveric and modelled knees. For the experimental results, strains for ACL, PCL, MCL, DMCL, LCL, and POPL ligament bundles were output as part of the IK process over the full flexion range. The corresponding predicted ligament strains for the FD analysis were also output for all models. As the kinematics are validated, it is assumed that comparing the corresponding ligament strains is a valid approach.
Data analysis
For validation, the FD model outputs of a single cadaveric specimen were compared to the IK results for that same specimen. The model outputs were discretised using MATLAB into 1° increments for easier pointwise comparison.
Since the number of specimens involved in the validation process was n = 1, a normalised root mean square error (NRMSE) criterion was used to compare FD results with corresponding IK results to address our hypothesis that the model accurately represents passive knee motions and ligament strains. RMSE was normalised by the range value observed for each variable. NRMSE is given by the following equation [32]:
where \(o\) is the observed value of the parameter, \(p\) is the predicted value according to the model, and N is the total number of observations in the validation dataset.
Commonly, an NRMSE of less than 10% is considered acceptable [31, 32]. However, this criterion can vary depending on the application [84]. In this study, cadaveric data are compared with a computational model very sensitive to various modelling parameters, i.e. contact surfaces, ligament material properties, slack lengths, and insertion points. As these parameters are determined in an ad hoc manner in this study, we consider NRMSE values < 0.3 (30%) to indicate acceptable performance [85, 86]. Adopting this lenient criterion suggests the potential for a full validation under stricter criteria in future studies which will use computational methods to drive these sensitive parameters towards an optimal solution. This was beyond the scope of this study.
The ranges of FD tibiofemoral kinematics and ligament strains of our musculoskeletal models were compared to their corresponding IK results to assess the remaining reconstructed knee model outcomes. The internal–external rotation pattern of the tibia relative to the femur under 5 Nm of IR/ER torque [25, 26] across the 0–100° of passive knee flexion was measured. The strain borne by knee ligaments (ACL, PCL, MCL, DMCL, LCL, and POPL) during passive knee flexion from 0 to 100° with applied 5 Nm IR torque [25, 26] was also measured. NRMSE values were then calculated.
Availability of data and materials
The datasets generated and analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- ACL:
-
Anterior cruciate ligament
- ALL:
-
Anterolateral ligament
- ACLR:
-
Anterior cruciate ligament reconstruction
- ACLR + ALLR:
-
ACLR combined with the ALL-reconstructed knee
- ACLR + Mac:
-
ACLR combined with modified MacIntosh
- ACLR + Ell:
-
ACLR combined with modified Ellison
- ACLR + DL:
-
ACLR combined with deep-Lemaire
- NRMSE:
-
Normalised root mean square error
- RMSE:
-
Root mean square error
- IK:
-
Inverse kinematics
- FD:
-
Forward dynamics
- IR:
-
Internal rotation
- ER:
-
External rotation
- DoF:
-
Degrees of freedom
- aACL:
-
Anterior bundle of anterior cruciate ligament
- pACL:
-
Posterior bundle of anterior cruciate ligament
- PCL:
-
Posterior cruciate ligament
- aPCL:
-
Anterior bundle of posterior cruciate ligament
- pPCL:
-
Posterior bundle of posterior cruciate ligament
- MCL:
-
Medial collateral ligament
- aMCL:
-
Anterior bundle of medial collateral ligament
- iMCL:
-
Inferior bundle of medial collateral ligament
- pMCL:
-
Posterior bundle of medial collateral ligament
- DMCL:
-
Deep medial collateral ligament
- aDMCL:
-
Anterior bundle of deep medial collateral ligament
- pDMCL:
-
Posterior bundle of deep medial collateral ligament
- LCL:
-
Lateral collateral ligament
- POPL:
-
Popliteofibular ligament
- CAP:
-
Posterior capsule
- CAPa:
-
Anterior bundle of posterior capsule
- CAPo:
-
Oblique bundle of posterior capsule
- CAPm:
-
Medial bundle of posterior capsule
- CAPl:
-
Lateral bundle of posterior capsule
- PT:
-
Patellar tendon
- cPT:
-
Central bundle of patellar tendon
- mPT:
-
Medial bundle of patellar tendon
- lPT:
-
Lateral bundle of patellar tendon
- ITB:
-
Iliotibial band
- GT:
-
Gerdy's tubercle
References
Robinson JR, Haddad FS. ACL graft failure: surgical technique may affect outcomes. Br Editorial Soc Bone Jt Surg London. 2021;103:1439–41.
Van Eck CF, Schkrohowsky JG, Working ZM, Irrgang JJ, Fu FH. Prospective analysis of failure rate and predictors of failure after anatomic anterior cruciate ligament reconstruction with allograft. Am J Sports Med. 2012;40(4):800–7.
Kocher MS, Steadman JR, Briggs K, Zurakowski D, Sterett WI, Hawkins RJ. Determinants of patient satisfaction with outcome after anterior cruciate ligament reconstruction. JBJS. 2002;84(9):1560–72.
Amis AA, Bull AM, Lie DT. Biomechanics of rotational instability and anatomic anterior cruciate ligament reconstruction. Oper Tech Orthop. 2005;15(1):29–35.
Tashiro Y, Okazaki K, Miura H, Matsuda S, Yasunaga T, Hashizume M, et al. Quantitative assessment of rotatory instability after anterior cruciate ligament reconstruction. Am J Sports Med. 2009;37(5):909–16.
Lee MC, Seong SC, Lee S, Chang CB, Park YK, Jo H, et al. Vertical femoral tunnel placement results in rotational knee laxity after anterior cruciate ligament reconstruction. Arthrosc J Arthrosc Relat Surg. 2007;23(7):771–8.
Bedi A, Maak T, Musahl V, Citak M, O’Loughlin PF, Choi D, et al. Effect of tibial tunnel position on stability of the knee after anterior cruciate ligament reconstruction: is the tibial tunnel position most important? Am J Sports Med. 2011;39(2):366–73.
Kondo E, Merican AM, Yasuda K, Amis AA. Biomechanical comparison of anatomic double-bundle, anatomic single-bundle, and nonanatomic single-bundle anterior cruciate ligament reconstructions. Am J Sports Med. 2011;39(2):279–88.
Sastre S, Popescu D, Núñez M, Pomes J, Tomas X, Peidro L. Double-bundle versus single-bundle ACL reconstruction using the horizontal femoral position: a prospective, randomized study. Knee Surg Sports Traumatol Arthrosc. 2010;18(1):32–6.
Tashman S, Collon D, Anderson K, Kolowich P, Anderst W. Abnormal rotational knee motion during running after anterior cruciate ligament reconstruction. Am J Sports Med. 2004;32(4):975–83.
Vadalà AP, Iorio R, De Carli A, Bonifazi A, Iorio C, Gatti A, et al. An extra-articular procedure improves the clinical outcome in anterior cruciate ligament reconstruction with hamstrings in female athletes. Int Orthop. 2013;37(2):187–92.
Sonnery-Cottet B, Daggett M, Helito CP, Fayard J-M, Thaunat M. Combined anterior cruciate ligament and anterolateral ligament reconstruction. Arthrosc Tech. 2016;5(6):e1253–9.
Rezende FC, de Moraes VY, Martimbianco ALC, Luzo MV, da Silveira Franciozi CE, Belloti JC. Does combined intra-and extraarticular ACL reconstruction improve function and stability? A meta-analysis. Clin Orthopaedic Relat Res. 2015;473(8):2609–18.
DePhillipo NN, Cinque ME, Chahla J, Geeslin AG, LaPrade RF. Anterolateral ligament reconstruction techniques, biomechanics, and clinical outcomes: a systematic review. Arthrosc J Arthrosc Relat Surg. 2017;33(8):1575–83.
Sonnery-Cottet B, Thaunat M, Freychet B, Pupim BHB, Murphy CG, Claes S. Outcome of a combined anterior cruciate ligament and anterolateral ligament reconstruction technique with a minimum 2-year follow-up. Am J Sports Med. 2015;43(7):1598–605.
Rezende FC, de Moraes VY, Martimbianco ALC, Luzo MV, da Silveira Franciozi CE, Belloti JC. Does combined intra- and extraarticular ACL reconstruction improve function and stability? A meta-analysis. Clin Orthopaed Relat Res. 2015;473(8):2609–18.
Devitt BM, Bell SW, Ardern CL, Hartwig T, Porter TJ, Feller JA, et al. The role of lateral extra-articular tenodesis in primary anterior cruciate ligament reconstruction: a systematic review with meta-analysis and best-evidence synthesis. Orthop J Sports Med. 2017;5(10):2325967117731767.
Ellison A. Distal iliotibial-band transfer for anterolateral rotatory instability of the knee. J Bone Joint Surg Am. 1979;61(3):330–7.
Barrett D, Mackenney R. MacIntosh-Jones reconstruction for the unstable knee. Injury. 1991;22(4):282–6.
Lemaire M. Ruptures anciennes du ligament croise anterieur du genou. J Chir. 1967;93(3):311–20.
DePhillipo NN, Cinque ME, Chahla J, Geeslin AG, LaPrade RF. Anterolateral ligament reconstruction techniques, biomechanics, and clinical outcomes: a systematic review. Arthroscopy. 2017;33(8):1575–83.
Nitri M, Rasmussen MT, Williams BT, Moulton SG, Cruz RS, Dornan GJ, et al. An in vitro robotic assessment of the anterolateral ligament, part 2: anterolateral ligament reconstruction combined with anterior cruciate ligament reconstruction. Am J Sports Med. 2016;44(3):593–601.
Ferretti A. Extra-articular reconstruction in the anterior cruciate ligament deficient knee: a commentary. Joints. 2014;2(1):41.
Kang K-T, Kim S-H, Son J, Lee YH, Koh Y-G. Validation of a computational knee joint model using an alignment method for the knee laxity test and computed tomography. Bio-Med Mater Eng. 2017;28(4):417–29.
Neri T, Cadman J, Beach A, Grasso S, Dabirrahmani D, Putnis S, et al. Lateral tenodesis procedures increase lateral compartment pressures more than anterolateral ligament reconstruction, when performed in combination with ACL reconstruction: a pilot biomechanical study. J ISAKOS. 2021;6(2):66–73.
Neri T, Dabirrahmani D, Beach A, Putnis S, Oshima T, Cadman J, et al. A biomechanical comparison of the main anterolateral procedures used in combination with anterior cruciate ligament reconstruction. Orthopaed J Sports Med. 2020;8(2):2325967120S00002.
Abebe ES, Utturkar GM, Taylor DC, Spritzer CE, Kim JP, Moorman CT, et al. The effects of femoral graft placement on in vivo knee kinematics after anterior cruciate ligament reconstruction. J Biomech. 2011;44(5):924–9.
Andriacchi TP, Briant PL, Bevill SL, Koo S. Rotational changes at the knee after ACL injury cause cartilage thinning. Clin Orthopaed Relat Res. 2006;442:39–44.
Blankevoort L, Huiskes R, de Lange A. The envelope of passive knee joint motion. J Biomech. 1988;21(9):705–20.
Blache Y, Dumas R, de Guise J, Saithna A, Sonnery-Cottet B, Thaunat M. Technical considerations in lateral extra-articular reconstruction coupled with anterior cruciate ligament reconstruction: a simulation study evaluating the influence of surgical parameters on control of knee stability. Clin Biomech. 2019;61:136–43.
Schermelleh-Engel K, Moosbrugger H, Müller H. Evaluating the fit of structural equation models: tests of significance and descriptive goodness-of-fit measures. Methods Psychol Res Online. 2003;8:23–74.
Browne MW, Cudeck R. Alternative ways of assessing model fit. Sociolog Methods Res. 1992;21(2):230–58.
Wilson DR, Feikes JD, Zavatsky AB, O’Connor JJ. The components of passive knee movement are coupled to flexion angle. J Biomech. 2000;33(4):465–73.
Erdemir A. Open knee: open source modeling and simulation in knee biomechanics. J Knee Surg. 2016;29(02):107–16.
Moglo KE, Shirazi-Adl A. Cruciate coupling and screw-home mechanism in passive knee joint during extension–flexion. J Biomech. 2005;38(5):1075–83.
Guess TM, Liu H, Bhashyam S, Thiagarajan G. A multibody knee model with discrete cartilage prediction of tibio-femoral contact mechanics. Comput Methods Biomech Biomed Eng. 2013;16(3):256–70.
Machado M, Flores P, Ambrosio J, Completo A. Influence of the contact model on the dynamic response of the human knee joint. Proc Inst Mech Eng Part K J Multi-body Dyn. 2011;225(4):344–58.
Clouthier AL, Smith CR, Vignos MF, Thelen DG, Deluzio KJ, Rainbow MJ. The effect of articular geometry features identified using statistical shape modelling on knee biomechanics. Med Eng Phys. 2019;66:47–55.
Gerus P, Sartori M, Besier TF, Fregly BJ, Delp SL, Banks SA, et al. Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces. J Biomech. 2013;46(16):2778–86.
Giffin JR, Vogrin TM, Zantop T, Woo SL-Y, Harner CD. Effects of increasing tibial slope on the biomechanics of the knee. Am J Sports Med. 2004;32(2):376–82.
Schmitz A, Piovesan D. Development of an open-source, discrete element knee model. IEEE Trans Biomed Eng. 2016;63(10):2056–67.
Marieswaran M, Sikidar A, Goel A, Joshi D, Kalyanasundaram D. An extended OpenSim knee model for analysis of strains of connective tissues. Biomed Eng Online. 2018;17(1):42.
Blankevoort L, Huiskes R, de Lange A. Recruitment of knee joint ligaments. J Biomech Eng. 1991;113(1):94–103.
Sahanand S, Jose A, Kumar G, Rajan DV. Interlinked Hamstrings for combined anterolateral and anterior cruciate ligament reconstruction: a novel technique for ALL. Indian J Orthopaed. 2021;56:621–7.
Devitt BM, Lord BR, Williams A, Amis AA, Feller JA. Biomechanical assessment of a distally fixed lateral extra-articular augmentation procedure in the treatment of anterolateral rotational laxity of the knee. Am J Sports Med. 2019;47(9):2102–9.
Geeslin AG, Moatshe G, Chahla J, Kruckeberg BM, Muckenhirn KJ, Dornan GJ, et al. Anterolateral knee extra-articular stabilizers: a robotic study comparing anterolateral ligament reconstruction and modified lemaire lateral extra-articular tenodesis. Am J Sports Med. 2018;46(3):607–16.
Blankevoort L, Huiskes R. Ligament–bone interaction in a three-dimensional model of the knee. J Biomech Eng. 1991;113(3):263–9.
Kennedy A, Coughlin D, Metzger M, Tang R, Pearle A, Lotz J, et al. Biomechanical evaluation of pediatric anterior cruciate ligament reconstruction techniques. Am J Sports Med. 2011;39:964–71.
Reichl I, Auzinger W, Schmiedmayer H-B, Weinmuller E. Reconstructing the knee joint mechanism from kinematic data. Math Comput Model Dyn Syst. 2010;16:403–15.
Wang X, Bennell KL, Wang Y, Wrigley TV, Van Ginckel A, Fortin K, et al. Tibiofemoral joint structural change from 2.5 to 4.5 years following ACL reconstruction with and without combined meniscal pathology. BMC Musculoskelet Disord. 2019;20(1):1–11.
Zantop T, Schumacher T, Diermann N, Schanz S, Raschke MJ, Petersen W. Anterolateral rotational knee instability: role of posterolateral structures. Arch Orthop Trauma Surg. 2007;127(9):743–52.
Wilson WT, Deakin AH, Payne AP, Picard F, Wearing SC. Comparative analysis of the structural properties of the collateral ligaments of the human knee. J Orthopaed Sports Phys Ther. 2012;42(4):345–51.
Bowman KF Jr, Sekiya JK. Anatomy and biomechanics of the posterior cruciate ligament, medial and lateral sides of the knee. Sports Med Arthrosc Rev. 2010;18(4):222–9.
Moglo KE, Shirazi-Adl A. Biomechanics of passive knee joint in drawer: load transmission in intact and ACL-deficient joints. Knee. 2003;10(3):265–76.
Knapp A, Williams LN. Predicting the effect of localized ACL damage on neighbor ligament mechanics via finite element modeling. Bioengineering. 2022;9(2):54.
Van de Velde SK, DeFrate LE, Gill TJ, Moses JM, Papannagari R, Li G. The effect of anterior cruciate ligament deficiency on the in vivo elongation of the medial and lateral collateral ligaments. Am J Sports Med. 2007;35(2):294–300.
Hu J, Xin H, Chen Z, Zhang Q, Peng Y, Jin Z. The role of menisci in knee contact mechanics and secondary kinematics during human walking. Clin Biomech. 2019;61:58–63.
Guess TM, Thiagarajan G, Kia M, Mishra M. A subject specific multibody model of the knee with menisci. Med Eng Phys. 2010;32(5):505–15.
Farshidfar SS, Cadman J, Deng D, Appleyard R, Dabirrahmani D. The effect of modelling parameters in the development and validation of knee joint models on ligament mechanics: a systematic review. PLoS ONE. 2022;17(1): e0262684.
Aagaard H, Verdonk R. Function of the normal meniscus and consequences of meniscal resection. Scand J Med Sci Sports. 1999;9(3):134–40.
Rao AJ, Erickson BJ, Cvetanovich GL, Yanke AB, Bach BR Jr, Cole BJ. The meniscus-deficient knee: biomechanics, evaluation, and treatment options. Orthop J Sports Med. 2015;3(10):2325967115611386.
Neri T, Testa R, Laurendon L, Dehon M, Putnis S, Grasso S, et al. Determining the change in length of the anterolateral ligament during knee motion: a three-dimensional optoelectronic analysis. Clin Biomech (Bristol, Avon). 2019;62:86–92.
Delp SL, Anderson FC, Arnold AS, Loan P, Habib A, John CT, et al. OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans Biomed Eng. 2007;54(11):1940–50.
Xu H, Bloswick D, Merryweather A. An improved OpenSim gait model with multiple degrees of freedom knee joint and knee ligaments. Comput Methods Biomech Biomed Engin. 2015;18(11):1217–24.
Basso O, Johnson D, Amis A. The anatomy of the patellar tendon. Knee Surg Sports Traumatol Arthrosc. 2001;9(1):2–5.
LaPrade RF, Morgan PM, Wentorf FA, Johansen S, Engebretsen L. The anatomy of the posterior aspect of the knee: an anatomic study. JBJS. 2007;89(4):758–64.
Meister BR, Michael SP, Moyer RA, Kelly JD, Schneck CD. Anatomy and kinematics of the lateral collateral ligament of the knee. Am J Sports Med. 2000;28(6):869–78.
Davies H, Unwin A, Aichroth P. The posterolateral corner of the knee: anatomy, biomechanics and management of injuries. Injury. 2004;35(1):68–75.
Claes S, Vereecke E, Maes M, Victor J, Verdonk P, Bellemans J. Anatomy of the anterolateral ligament of the knee. J Anat. 2013;223(4):321–8.
Shelburne KB, Torry MR, Pandy MG. Contributions of muscles, ligaments, and the ground-reaction force to tibiofemoral joint loading during normal gait. J Orthop Res. 2006;24(10):1983–90.
Petersen W, Zantop T. Anatomy of the anterior cruciate ligament with regard to its two bundles. Clin Orthopaed Relat Res. 2007;454:35–47.
Kopf S, Musahl V, Tashman S, Szczodry M, Shen W, Fu FH. A systematic review of the femoral origin and tibial insertion morphology of the ACL. Knee Surg Sports Traumatol Arthrosc. 2009;17(3):213–9.
Girgis FG, Marshall JL, Monajem A. The cruciate ligaments of the knee joint. Anatomical, functional and experimental analysis. Clin Orthopaed Relat Res. 1975;106:216–31.
Edwards A, Bull AM, Amis AA. The attachments of the fiber bundles of the posterior cruciate ligament: an anatomic study. Arthrosc J Arthrosc Relat Surg. 2007;23(3):284–90.
Blankevoort L, Kuiper J, Huiskes R, Grootenboer H. Articular contact in a three-dimensional model of the knee. J Biomech. 1991;24(11):1019–31.
Shelburne KB, Kim HJ, Sterett WI, Pandy MG. Effect of posterior tibial slope on knee biomechanics during functional activity. J Orthop Res. 2011;29(2):223–31.
Wismans J, Veldpaus F, Janssen J, Huson A, Struben P. A three-dimensional mathematical model of the knee-joint. J Biomech. 1980;13(8):677–85.
Hamner DL, Brown CH, Steiner ME, Hecker AT, Hayes WC. Hamstring tendon grafts for reconstruction of the anterior cruciate ligament: biomechanical evaluation of the use of multiple strands and tensioning techniques. JBJS. 1999;81(4):549–57.
Benos L, Stanev D, Spyrou L, Moustakas K, Tsaopoulos DE. A review on finite element modeling and simulation of the anterior cruciate ligament reconstruction. Front Bioeng Biotechnol. 2020;8(967).
Chahla J, Menge TJ, Mitchell JJ, Dean CS, LaPrade RF. Anterolateral ligament reconstruction technique: an anatomic-based approach. Arthrosc Tech. 2016;5(3):e453–7.
Neri T, Palpacuer F, Testa R, Bergandi F, Boyer B, Farizon F, et al. The anterolateral ligament: anatomic implications for its reconstruction. Knee. 2017;24(5):1083–9.
Kennedy MI, Claes S, Fuso FAF, Williams BT, Goldsmith MT, Turnbull TL, et al. The anterolateral ligament: an anatomic, radiographic, and biomechanical analysis. Am J Sports Med. 2015;43(7):1606–15.
Daggett M, Ockuly AC, Cullen M, Busch K, Lutz C, Imbert P, et al. Femoral origin of the anterolateral ligament: an anatomic analysis. Arthroscopy. 2016;32(5):835–41.
Gu W, Pandy MG. Direct validation of human knee-joint contact mechanics derived from subject-specific finite-element models of the tibiofemoral and patellofemoral joints. J Biomech Eng. 2020;142(7).
Jamieson P, Porter J, Wilson D. A test of the computer simulation model ARCWHEAT1 on wheat crops grown in New Zealand. Field Crop Res. 1991;27(4):337–50.
Akumaga U, Tarhule A, Yusuf AA. Validation and testing of the FAO AquaCrop model under different levels of nitrogen fertilizer on rainfed maize in Nigeria, West Africa. Agric Forest Meteorol. 2017;232:225–34.
Acknowledgements
This study forms one part of SSF's PhD thesis. Therefore, the authors would like to acknowledge Macquarie University for funding this study under Macquarie University Research Excellence Scheme (MQRES) scholarship.
Funding
This study forms one part of SSF's PhD thesis. Therefore, the authors would like to acknowledge Macquarie University (www.mq.edu.au) for funding this study under Macquarie University Research Excellence Scheme (MQRES) scholarship. The Macquarie University had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author information
Authors and Affiliations
Contributions
Conceptualisation: SSF, DD, JC, TN. Data curation: SSF, DD, JC. Formal analysis: SSF, DD, JC. Investigation: SSF, DD, JC. Methodology: SSF, DD, JC, DP, TN. Project administration: DD, RA, DP. Resources: SSF, DD. Software: SSF, DD, JC. Supervision: DD, RA, DP. Validation: SSF, DD, JC. Visualisation: SSF, DD, JC. Writing—original draft preparation: SSF, DD, JC. Writing—review & editing: SSF, DD, JC, RA, DP, TN. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate.
Not applicable.
Consent for publication
All authors have given consent for publication.
Competing interests
The authors have no competing interest to declare.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
: Figure S1: Subject-specific intact knee model created in OpenSim [5], including one 6-DoF tibiofemoral joint and one 3-DoF patellofemoral joint. Figure S2: Representation of tibiofemoral contact surfaces through different developed tibial contact geometries (the tibial plateau): planar objects (A), curvature objects (B), and subject-specific objects(C). Figure S3: Representation of steps developing the subject-specific tibial contact surface (the tibial plateau), including, Boolean subtraction (A), meshing the 3D object (B), and cropping/smoothing (C). Figure S4: Representation of three wrapping objects included in the knee model placed at the medial epicondyle, lateral epicondyle, and patellofemoral joint. Table S1. Wrapping object parameters.
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
Farshidfar, S.S., Cadman, J., Neri, T. et al. Towards a validated musculoskeletal knee model to estimate tibiofemoral kinematics and ligament strains: comparison of different anterolateral augmentation procedures combined with isolated ACL reconstructions. BioMed Eng OnLine 22, 31 (2023). https://doi.org/10.1186/s12938-023-01094-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12938-023-01094-y