Applied Mathematics
Vol.4 No.5A(2013), Article ID:31913,7 pages DOI:10.4236/am.2013.45A011

The Effect of Ligament Modeling Technique on Knee Joint Kinematics: A Finite Element Study

Ata M. Kiapour1,2, Vikas Kaul1,2, Ali Kiapour1,2, Carmen E. Quatman3, Samuel C. Wordeman3, Timothy E. Hewett3, Constantine K. Demetropoulos1,2, Vijay K. Goel1,2*

1Engineering Center for Orthopaedic Research Excellence (ECORE), University of Toledo, Toledo, USA

2Departments of Orthopaedics and Bioengineering, The University of Toledo, Toledo, USA

3Sports Health and Performance Institute (SHPI), The Ohio State University, Columbus, USA

Email: *vijay.goel@utoledo.edu

Copyright © 2013 Ata M. Kiapour et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received March 1, 2013; revised April 2, 2013; accepted April 9, 2013

Keywords: Finite Element; Knee; Biomechanics; Constitutive Model

ABSTRACT

Finite element (FE) analysis has become an increasingly popular technique in the study of human joint biomechanics, as it allows for detailed analysis of the joint/tissue behavior under complex, clinically relevant loading conditions. A wide variety of modeling techniques have been utilized to model knee joint ligaments. However, the effect of a selected constitutive model to simulate the ligaments on knee kinematics remains unclear. The purpose of the current study was to determine the effect of two most common techniques utilized to model knee ligaments on joint kinematics under functional loading conditions. We hypothesized that anatomic representations of the knee ligaments with anisotropic hyperelastic properties will result in more realistic kinematics. A previously developed, extensively validated anatomic FE model of the knee developed from a healthy, young female athlete was used. FE models with 3D anatomic and simplified uniaxial representations of main knee ligaments were used to simulate four functional loading conditions. Model predictions of tibiofemoral joint kinematics were compared to experimental measures. Results demonstrated the ability of the anatomic representation of the knee ligaments (3D geometry along with anisotropic hyperelastic material) in more physiologic prediction of the human knee motion with strong correlation (r ≥ 0.9 for all comparisons) and minimum deviation (0.9˚ ≤ RMSE ≤ 2.29˚) from experimental findings. In contrast, non-physiologic uniaxial elastic representation of the ligaments resulted in lower correlations (r ≤ 0.6 for all comparisons) and substantially higher deviation (2.6˚ ≤ RMSE ≤ 4.2˚) from experimental results. Findings of the current study support our hypothesis and highlight the critical role of soft tissue modeling technique on the resultant FE predicted joint kinematics.

1. Introduction

The knee is the largest and one of the most complex joints within the human body, consisting of both patellofemoral and tibiofemoral articulations. Anatomical structures such as ligaments, menisci and articular cartilage provide stability across the knee joint during functional daily activities. However, abnormalities due to age, injury, disease and other factors can affect biomechanical function of the knee joint. Mechanistic computational models, if properly validated, can serve as effective tools in parametric analyses, as well as population-based clinical studies. In particular, the use of finite element (FE) analysis has became progressively popular in the study of joint biomechanics as it allows for detailed analysis of the joint/tissue behavior under complexclinically relevant loading conditions. FE methods have provided considerable insight into knee joint biomechanics, including ligament function, ligament reconstruction technique, and implant design. Due to inherent challenges associated with experiments (in vivo and ex vivo) and the associated high cost and time, FE analysis has long been recognized and trusted as a reliable alternative method in the study of human joints. Primary advantage of this numerical approach lies in precise control over boundary conditions, material properties and structural alterations in parametric studies. Moreover, the ligament forces/strains, contact forces/areas, and stress/strain distribution across soft and hard tissue structures are invaluable products of such a numerical approach, which are challenging, if not impossible, to obtain experimentally. The reliability of FE models strongly relies upon: a) appropriate representation of the geometry and assigned material properties, b) realistic simulation of interactions, constraints and boundary conditions, and finally c) thorough validation against experimental data.

Ligaments are soft connective tissues with a composite structure that connect bones together. As the main contributor to the overall joint stability, the mechanical function of these connective tissue structures is to guide normal joint motion and restrict abnormal joint movement. This is assisted by the topology of the articulating surfaces, muscle forces and other soft tissue constraints such as joint capsule. Physiologic characteristic of soft tissue material composition has always challenged the accuracy of the simplified numerical models of anatomical joints, specifically the knee joint which has been the scope of numerous studies due to its critical role in stability of human body during various physiological activities. A wide variety of modeling techniques have been utilized to model knee joint ligaments [1-12]. In majority of earlier FE studies of the knee joint, uniaxial discrete line elements (truss or spring) with simplified material properties were used to model ligaments [1-4,7]. Such an assumption of soft tissue geometry is associated with shortcomings such as inability to predict non-uniform 3-dimensional (3D) stresses and strains across the tissue [6, 8,11]. Using Image processing techniques, ligaments were modeled with a 3D reconstructed geometry coupled with isotropic hyperelastic constitutive material models [9,10,12]. More recently, transversely isotropic hyperelastic constitutive models were developed and used to study knee ligaments [5,6,8,11].

Considering the critical role of ligaments in providing joint stability and associated assumptions with each modeling techniques in the characterization of tissue material properties, joint kinematics are expected to differ. However, the effect of selected constitutive model to simulate the ligaments on knee kinematics remains unclear. Hence, this study was designed to investigate alterations in knee joint kinematics under functional loading resulting from two different ligament modeling techniques: 1) uniaxial representation with isotropic nonlinear elastic properties and 2) anatomic 3D representation with anisotropic hyperelastic properties. We hypothesized that anatomic representations of the knee ligaments with anisotropic hyperelastic material property will result in more realistic kinematics.

2. Methods

2.1. Model Development

Following IRB approval, computerized tomography (CT) and magnetic resonance imaging (MRI) scans of a young adult female athlete’s lower limb (Age: 25 years, Height: 170 cm, Weight: 64.4 Kg) were used to capture bony and soft tissue geometry, respectively. Scans were obtained while the subject was supine with the leg in an unloaded neutral position. CT and MRI scans were co-registered for bony and soft tissue alignment. 3D geometry of the pelvis, leg (upper and lower) and foot segments were reconstructed from high resolution CT images in all three anatomical planes. Sagittal, coronal and axial MR images of the left knee were used to generate the 3D geometry of the knee articular cartilage, menisci, and knee cruciate and collateral ligaments. These geometries were then converted into solid 8-node hexahedral elements and subsequently imported into the ABAQUS FE package v6.11 (SIMULIA, Providence, RI, USA) to generate the FE model (Figure 1). While cruciate and collateral ligaments, articular cartilage and menisci were modeled as 3D structures, the rest of the simulated knee ligaments, joint capsule and muscle tendons were modeled as uniaxial truss elements (Figure 2).

To optimize computational expense, pelvis, proximal femur (from 10 cm above the joint line), distal tibia (from 10 cm below the joint line), fibula and foot were

Figure 1. FE model development steps.

Figure 2. Developed FE model of lower extremity. (ACL: anterior cruciate ligament; PCL: posterior cruciate ligament; LCL: lateral collateral ligament; sMCL, dMCL and oMCL: superficial, deep and oblique bundles of medial collateral ligament; CAPm, CAPl, CAPo and CAPa: medial, lateral, oblique popliteal and arcuate popliteal bundles of posterior capsule; ALS: anterolateral structure; PFL: popliteofibular ligament; MPFL: medial patellofemoral ligament; LPFL: lateral patellofemoral ligament; PT: patellar tendon)

modeled as rigid bodies, while the remaining structures were considered deformable. Following assembly, proper material properties taken from literature were assigned to each segment [1,2,5,13-20]. Bones were modeled as linear elastic [21-25] with different moduli assigned to cortical and trabecular regions consistent with earlier FE studies of the human knee joint [2,13]. Tibiofemoral and patellofemoral articular cartilage were modeled as isotropic linear elastic [16]. Menisci were modeled as transversely isotropic linear elastic with different mechanical properties in circumferential, axial and radial directions [14,15,26]. Horn-meniscus attachment was simulated with multiple linear elastic truss elements [13].

Knee cruciate and collateral ligaments were modeled as incompressible anisotropic hyperelastic structures using the Holzapfel-Gasser-Ogden (HGO) material model [27]. HGO model is a hyperelastic, anisotropic material model that was developed to model the criss-crossed fibrous soft tissues like the illiac adventitia [16]. Briefly, isotropic non-collagenous ground matrix is modeled by the incompressible hyperelastic neo-Hookean component of the strain energy density (SED) function, whereas the transversely isotropic fibrous component is modeled by the following function developed by Gasser et al. [27]:

where and are the respective isotropic and anisotropic components of the SED, is the mean orientation of the fibers, is the structure tensor, and is the dispersion parameter for the fiber family. A statistical distribution function allows for a spatial distribution of the fiber orientation. Fibrous component of the SED supports tensile loads only and is defined as:

where is the first invariant of and is a generalized structure tensor defined as:

The non-collagenous ground substance is modeled using the following incompressible isotropic neo-Hookean model:

Cruciate ligaments were modeled using two fiber families each in order to simulate bundles within ACL and posterior cruciate ligament (PCL) [17,18]. Both MCL (superficial bundle) and LCL were modeled using one family of fibers. Given the microstructure of the MCL and LCL, the HGO model was modified to account for a single family of fibers:

FE simulation of experimental uniaxial tensile tests along the longitudinal direction as per Butler et al. [28] for the ACL and PCL, and Quapp and Weiss [29] for the MCL were used to derive a series of coefficients for the constitutive model using a curve fitting technique (Figure 3). Coefficients for the lateral collateral ligament (LCL) were assumed to be identical to those of the MCL [5,11]. All other simulated knee ligaments were modeled as non-linear elastic, tension-only materials using truss elements with theoretically defined cross-sectional area. Further, 13 uniaxial truss-connector elements were used to simulate trans-knee muscle forces (Figure 2).

All other simulated knee ligaments were modeled as non-linear elastic, tension-only materials using truss elements with theoretically defined cross-sectional area. Further, 13 uniaxial truss-connector elements were used to simulate trans-knee muscle forces (Figure 2).

A frictionless surface-to-surface tangential contact with non-linear finite sliding interaction was used to simulate articular surfaces [3,11,13]. Since the current FE model was developed to investigate phenomena associated with knee biomechanics and relevant injuries, key knee joint soft tissue structures have been incorporated into the model. Both tibiofemoral and patellofemoral

Figure 3. FE predictions vs. experimental data of the uniaxial tensile test for ACL, PCL (Top) and MCL (Bottom).

joints were simulated as six degree-of-freedom (DOF) joints with their motion defined by their surrounding soft tissue constraints and the topology of the articular surfaces. The hip and ankle joints were simplified as virtual ball-and-socket joints controlled by imported kinematic data, while optimizing for computational efficiency. The kinematics of the hip, knee and ankle joints were defined using the local coordinate systems proposed by Grood and Suntay [30]. Subsequently, the model was extensively validated against direct experimental measures of tibiofemoral kinematics, ACL and MCL strains and tibiofemoral articular cartilage pressure distribution under a wide range of quasi-static and dynamic loading conditions [31].

2.2. Loading Profile

Four quasi-static loading conditions were simulated in order to compare the predicted FE kinematics with experimental measurements from an ex vivo study of 19 fresh frozen cadaveric legs [32,33]:

1) 0 to 50 Nm of knee abduction (at 25˚ of flexion), 2) 0 to 50 Nm of knee abduction + 20 Nm of internal tibial rotation (at 25˚ of flexion), 3) baseline (no external load, 0˚ - 90˚ of flexion), 4) 15 Nm of internal tibial rotation (0˚ - 90˚ of flexion), all under simulated muscle loads (quadriceps: 400 N and hamstrings: 200 N).

In order to study the effects of soft tissue material models, 3D reconstructed cruciate and collateral ligaments were substituted with multiple uniaxial representations (truss elements) with isotropic non-linear elastic material properties [2], while maintaining the same origins, insertions and initial orientation as the 3D model. Finally, the quasi-static simulations were repeated using uniaxial ligaments.

3. Results

3.1. Frontal Plane Kinematics (Valgus Rotation)

Both FE models resulted in similar frontal plane quality of motion as the experimental measurements under both singleand multi-axial loading conditions (Figure 4). Models also replicated coupled motion as observed in cadaveric experiments shown by knee valgus rotation under an additional internal tibial rotation moment of 20 Nm (Figure 4). The anatomic 3D representation of ligaments resulted in strong correlations (Pure abduction:, Combined abduction and internal rotation:) with minimum deviation (Pure abduction:, Combined abduction and internal rotation:) between FE model predictions and experimental measures of tibiofemoral frontal plane kinematics. Moreover, model predictions were within the range of 95% confidence intervals of average experi-

Figure 4. FE predictions Vs. experimental data for tibiofemoral frontal plane kinematics (Shaded area represent experimental 95% confidence intervals).

mental measurements. In contrast, the uniaxial assumption coupled with simplified constitutive model of the knee ligaments resulted in substantially lower correlations (Pure abduction:, Combined abduction and internal rotation:) and higher deviation (Pure abduction:, Combined abduction and internal rotation:) from the average experimentally quantified tibiofemoral kinematics. In addition to lower correlation and higher deviation from average experimental data, model predictions of joint kinematics were demonstrated to be outside the range of 95% confidence intervals of average experimental measurements (Figure 4).

3.2. Axial Plane Kinematics (Internal Rotation)

Both FE models demonstrated similar trends as the experimental measurements under both singleand multiaxial loading conditions (Figure 5). Models also replicated knee joint screw-home mechanism [34] as observed in cadaveric experiments shown by internal tibial rotation during the early phase of flexion (Figure 5). The anatomic 3D representation of ligaments resulted in strong correlations (Baseline:, Internal rotation:) with minimum deviation (Baseline: RMSE = 1.1˚, Internal rotation: RMSE = 2.2˚) between FE model predictions and experimental measures of tibiofemoral axial plane kinematics. Moreover, model predictions were within the range of 95% confidence intervals of average experimental measurements. In contrast, the uniaxial assumption coupled with simplified constitutive model of the knee ligaments resulted in substantially lower correlations (Baseline:, Internal rotation:) and higher deviation (Baseline:, Internal rotation:) from the average experimentally quantified tibiofemoral kinematics. In addition to lower correlation and higher deviation from average experimental data, model predictions of joint kinematics were demonstrated to be outside the range of 95% confidence intervals of average experimental data (Figure 5).

4. Discussion

FE analysis is a powerful numerical technique that makes it feasible to investigate the biomechanical behavior of complex biological structures. During the past three decades, a large number of knee FE models with varying degrees of complexity, accuracy and functionality have been reported in the literature [1-12]. Simplified uniaxial representations of ligaments coupled with non-physiologic constitutive material models have been associated with the majority of these models [1-4,7]. More recent studies have used a 3D representation of knee ligaments with various degrees of anatomical and constitutive model complexity [5,6,8,9,11,12,28]. Song and colleagues developed a 3D FE model of the tibiofemoral joint which included 3D representation of the femur, tibia and ACL (with distinct AM and PL bundles) modeled as an isotropic hyperelastic material [9]. A similar model was developed by Gardiner and Weiss in order to study MCL biomechanics under functional loading [6]. They utilized a novel transversely isotropic, incompressible hyperelastic material model in order to simulate the MCL (superficial bundle) as a composite soft tissue structure [6]. Limbert et al. used a similar constitutive material model to study ACL biomechanics under passive tibial translation and flexion in a 3D FE model of an isolated ACL [8]. Others [5,11] have used similar constitutive modeling approaches with 3D simulations of key knee ligaments incorporated in 3D FE models of the knee joint. Despite substantial research efforts to develop soft tissue constitutive material models, little is known about the effects of

Figure 5. FE predictions Vs. experimental data for tibiofemoral axial plane kinematics (Shaded area represent experimental 95% confidence intervals).

such techniques on resultant joint function.

The purpose of the current study was to determine the effect of two most common techniques utilized to model knee ligaments on joint kinematics under functional loading conditions. A previously developed, extensively validated anatomic FE model of the knee developed from a healthy, young female athlete was used. FE models with 3D anatomic and simplified uniaxial representations of main knee ligaments (ACL, PCL, MCL and LCL) were recruited to simulate four quasi-static loading conditions as conducted in the cadaveric experiments.

The 3D anisotropic hyperelastic model resulted in a more physiologic prediction of the human knee motion under a range of singleand multi-planar functional loading conditions with strong correlation and minimal deviation from experimental data. In contrast, lower correlations in addition to notable deviations were observed using simplified uniaxial modeling technique. The current findings support our hypothesis and highlight the critical role of soft tissue modeling technique on resultant FE predicted joint kinematics. Anatomically accurate 3D representation of such structures coupled with structurally motivated constitutive models [27] facilitate implementation of realistic ligament mechanical properties such as finite deformation, anisotropy and non-linear incompressible fiber-reinforced structures. This approach also permits incorporation of realistic interactions between adjacent structures such as ligament-bone interaction that may also result in a more realistic simulation of lines of action as they vary with changes in joint orientation [8,11]. Moreover, anatomic representation of the ligament will also make it feasible to quantify local stress-strain distribution across the tissue, which is critical in study of ligament injury mechanisms.

5. Acknowledgements

The authors acknowledge funding support from the National Institutes of Health/National Institute of Arthritis and Musculoskeletal and Skin Diseases grants R01- AR049735 and R01-AR056259. The authors would also like to thank Dr. Fernando Cacho-Nerin for his assistance in anisotropic hyperelsatic material modeling.

REFERENCES

  1. E. M. Abdel-Rahman and M. S. Hefzy, “Three-Dimensional Dynamic Behaviour of the Human Knee Joint under Impact Loading,” Medical Engineering & Physics, Vol. 20, No. 4, 1998, pp. 276-290. doi:10.1016/S1350-4533(98)00010-1
  2. P. Beillas, G. Papaioannou, S. Tashman and K. H. Yang, “A New Method to Investigate in Vivo Knee Behavior Using a Finite Element Model of the Lower Limb,” Journal of Biomechanics, Vol. 37, No. 7, 2004, pp. 1019-1030. doi:10.1016/j.jbiomech.2003.11.022
  3. M. Z. Bendjaballah, A. Shirazi-Adl and D. J. Zukor, “Finite Element Analysis of Human Knee Joint in VarusValgus,” Clinical Biomechanics, Vol. 12, No. 3, 1997, pp. 139-148. doi:10.1016/S0268-0033(97)00072-7
  4. L. Blankevoort and R. Huiskes, “Validation of a ThreeDimensional Model of the Knee,” Journal of Biomechanics, Vol. 29, No. 7, 1996, pp. 955-961. doi:10.1016/0021-9290(95)00149-2
  5. Y. Y. Dhaher, T. H. Kwon and M. Barry, “The Effect of Connective Tissue Material Uncertainties on Knee Joint Mechanics under Isolated Loading Conditions,” Journal of Biomechanics, Vol. 43, No. 16, 2010, pp. 3118-3125. doi:10.1016/j.jbiomech.2010.08.005
  6. J. C. Gardiner, J. A. Weiss and T. D. Rosenberg, “Strain in the Human Medial Collateral Ligament during Valgus Loading of the Knee,” Clinical Orthopaedics and Related Research, 2001, pp. 266-274. doi:10.1097/00003086-200110000-00031
  7. G. Li, J. Gil, A. Kanamori and S. L. Woo, “A Validated Three-Dimensional Computational Model of a Human Knee Joint,” Journal of Biomechanical Engineering, Vo1. 21, No. 6, 1999, pp. 657-662. doi:10.1115/1.2800871
  8. G. Limbert, M. Taylor and J. Middleton, “Three-Dimensional Finite Element Modelling of the Human ACL: Simulation of Passive Knee Flexion with a Stressed and Stress-Free ACL,” Journal of Biomechanics, Vol. 37, No. 11, 2004, pp. 1723-1731. doi:10.1016/j.jbiomech.2004.01.030
  9. Y. Song, R. E. Debski, V. Musahl, M. Thomas and S. L. Woo, “A Three-Dimensional Finite Element Model of the Human Anterior Cruciate Ligament: A Computational Analysis with Experimental Validation,” Journal of Biomechanics, Vol. 37, No. 3, 2004, pp. 383-390. doi:10.1016/S0021-9290(03)00261-6
  10. N. A. Ramaniraka, A. Terrier, N. Theumann and O. Siegrist, “Effects of the Posterior Cruciate Ligament Reconstruction on the Biomechanics of the Knee Joint: A Finite Element Analysis,” Clinical Biomechanics, Vol. 20, No. 4, 2005, pp. 434-442. doi:10.1016/j.clinbiomech.2004.11.014
  11. E. Pena, B. Calvo, M. A. Martinez and M. Doblare, “A Three-Dimensional Finite Element Analysis of the Combined Behavior of Ligaments and Menisci in the Healthy Human Knee Joint,” Journal of Biomechanics, Vol. 39, No. 9, 2006, pp. 1686-1701. doi:10.1016/j.jbiomech.2005.04.030
  12. F. Xie, L. Yang, L. Guo, Z. J. Wang and G. Dai, “A Study on Construction Three-Dimensional Nonlinear Finite Element Model and Stress Distribution Analysis of Anterior Cruciate Ligament,” Journal of Biomechanical Engineering, Vol. 131, No. , 2009, Article ID: 121007. doi:10.1115/1.4000167
  13. T. L. Donahue, M. L. Hull, M. M. Rashid and C. R. Jacobs, “A Finite Element Model of the Human Knee Joint for the Study of Tibio-Femoral Contact,” Journal of Biomechanical Engineering, Vol. 124, No. 3, 2002, pp. 273- 280. doi:10.1115/1.1470171
  14. D. L. Skaggs, W. H. Warden and V. C. Mow, “Radial Tie Fibers Influence the Tensile Properties of the Bovine Medial Meniscus,” Journal of Orthopaedic Research, Vol. 12, No. 2, 1994, pp. 176-185. doi:10.1002/jor.1100120205
  15. M. Tissakht and A. M. Ahmed, “Tensile Stress-Strain Characteristics of the Human Meniscal Material,” Journal of Biomechanics, Vol. 28, No. 4, 1995, pp. 411-422. doi:10.1016/0021-9290(94)00081-E
  16. D. E. Shepherd and B. B. Seedhom, “The ‘Instantaneous’ Compressive Modulus of Human Articular Cartilage in Joints of the Lower Limb,” Rheumatology, Vol. 38, No. 2, 1999, pp. 124-132. doi:10.1093/rheumatology/38.2.124
  17. F. G. Girgis, J. L. Marshall and A. Monajem, “The Cruciate Ligaments of the Knee Joint. Anatomical, Functional and Experimental Analysis,” Clinical Orthopaedics and Related Research, 1975, pp. 216-231. doi:10.1097/00003086-197501000-00033
  18. S. L. Woo, A. Kanamori, J. Zeminski, M. Yagi, C. PapaGeorgiou and F. H. Fu, “The Effectiveness of Reconstruction of the Anterior Cruciate Ligament with Hamstrings and Patellar Tendon. A cadaveric Study Comparing Anterior Tibial and Rotational Loads,” The Journal of Bone & Joint Surgery, Vol. 84A, 2002, pp. 907-914.
  19. C. S. Shin, A. M. Chaudhari and T. P. Andriacchi, “The Influence of Deceleration Forces on ACL Strain during Single-Leg Landing: A Simulation Study,” Journal of Biomechanics, Vol. 40, No. 5, 2007, pp. 1145-1152. doi:10.1016/j.jbiomech.2006.05.004
  20. P. Atkinson, T. Atkinson, C. Huang and R. Doane, “A Comparison of the Mechanical and Dimensional Properties of the Human Medial and Lateral Patellofemoral Ligaments,” Proceedings of the 46th ORS Annual Meeting, Orlando, 12-15 March 2000. http://www.ors.org/Transactions/46/0776.pdf
  21. P. L. Mente and J. L. Lewis, “Elastic Modulus of Calcified Cartilage Is an Order of Magnitude Less than That of Subchondral Bone,” Journal of Orthopaedic Research, Vol. 12, No. 5, 1994, pp. 637-647. doi:10.1002/jor.1100120506
  22. F. Linde, “Elastic and Viscoelastic Properties of Trabecular Bone by a Compression Testing Approach,” Danish Medical Bulletin, Vol. 41, 1994, pp. 119-138.
  23. J. C. Lotz, T. N. Gerhart and W. C. Hayes, “Mechanical Properties of Metaphyseal Bone in the Proximal Femur,” Journal of Biomechanics, Vol. 24, No. 5, 1991, pp. 317- 329. doi:10.1016/0021-9290(91)90350-V
  24. J. L. Kuhn, S. A. Goldstein, M. J. Ciarelli and L. S. Matthews, “The Limitations of Canine Trabecular Bone as a Model for Human: A Biomechanical Study,” Journal of Biomechanics, Vol. 22, No. 2, 1989, pp. 95-107. doi:10.1016/0021-9290(89)90032-8
  25. S. A. Goldstein, “The Mechanical Properties of Trabecular Bone: Dependence on Anatomic Location and Function,” Journal of Biomechanics, Vol. 20, No. 11-12, 1987, pp. 1055-1061. doi:10.1016/0021-9290(87)90023-6
  26. J. Yao, P. D. Funkenbusch, J. Snibbe, M. Maloney and A. L. Lerner, “Sensitivities of Medial Meniscal Motion and Deformation to Material Properties of Articular Cartilage, Meniscus and Meniscal Attachments Using Design of Experiments Methods,” Journal of Biomechanical Engineering, Vol. 128, No. 3, 2006, pp. 399-408. doi:10.1115/1.2191077
  27. T. C. Gasser, R. W. Ogden and G. A. Holzapfel, “Hyperelastic Modelling of Arterial Layers with Distributed Collagen Fibre Orientations,” Journal of the Royal Society, Interface/the Royal Society, Vol. 3, No. 6, 2006, pp. 15-35. doi:10.1098/rsif.2005.0073
  28. D. L. Butler, M. Y. Sheh, D. C. Stouffer, V. A. Samaranayake and M. S. Levy, “Surface Strain Variation in Human Patellar Tendon and Knee Cruciate Ligaments,” Journal of Biomechanical Engineering, Vol. 112, No. 1, 1990, pp. 38-45. doi:10.1115/1.2891124
  29. K. M. Quapp and J. A. Weiss, “Material Characterization of Human Medial Collateral Ligament,” Journal of Biomechanical Engineering, Vol. 120, No. 6, 1998, pp. 757- 763. doi:10.1115/1.2834890
  30. E. S. Grood and W. J. Suntay, “A Joint Coordinate System for the Clinical Description of Three-Dimensional Motions: Application to the Knee,” Journal of Biomechanical Engineering, Vol. 105, No. 2, 1983, pp. 136-144. doi:10.1115/1.3138397
  31. A. Kiapour, A. M. Kiapour, V. Kaul, C. E. Quatman, R. C. Ditto, J. W. Levine, et al., “Finite Element Model of the Knee for Investigation of High Rate Injury Mechanisms: Development and Validation,” Proceedings of the 58th ORS Annual Meeting, San Francisco, 4-7 February 2012. http://online.ors.org/web/Transactions/58/0101.pdf
  32. A. M. Kiapour, C. E. Quatman, R. C. Ditto, J. W. Levine, S. C. Wordeman, T. E. Hewett, et al., “Influence of Axial Rotation Moments on ACL Strain: A Cadaveric Study of Singleand Multi-Axis Loading of the Knee,” Proceedings of the 35th ASB Annual Meeting, Long Beach, 10-13 August 2011. http://www.asbweb.org/conferences/2011/pdf/216.pdf
  33. A. M. Kiapour, C. E. Quatman, V. K. Goel, R. C. Ditto, S. C. Wordeman, J. W. Levine, et al., “Knee Articular Cartilage Pressure Distribution under Singleand Multi-Axis Loading Conditions: Implications for ACL Injury Mechanism,” Proceedings of the 36th ASB Annual Meeting, Geinsville, 15-18 August 2012. http://www.asbweb.org/conferences/2012/topics/ASB%202012%20oral.pdf
  34. M. A. Freeman and V. Pinskerova, “The Movement of the Normal Tibio-Femoral Joint,” Journal of Biomechanics, Vol. 38, No. 2, 2005, pp. 197-208. doi:10.1016/j.jbiomech.2004.02.006

NOTES

*Corresponding author.