Research Papers

Using a Statistically Calibrated Biphasic Finite Element Model of the Human Knee Joint to Identify Robust Designs for a Meniscal Substitute

[+] Author and Article Information
Erin R. Leatherman

Department of Statistics,
West Virginia University,
Morgantown, WV 26506
e-mail: erleatherman@mail.wvu.edu

Hongqiang Guo

Department of Biomechanics,
Hospital for Special Surgery,
New York, NY 10021
e-mail: GuoH@hss.edu

Susannah L. Gilbert

Department of Biomechanics,
Hospital for Special Surgery,
New York, NY 10021
e-mail: GilbertS@hss.edu

Ian D. Hutchinson

Department of Biomechanics,
Hospital for Special Surgery,
New York, NY 10021
e-mail: ihutchin@wakehealth.edu

Suzanne A. Maher

Department of Biomechanics,
Hospital for Special Surgery,
New York, NY 10021
e-mail: MaherS@hss.edu

Thomas J. Santner

Department of Statistics,
The Ohio State University,
Columbus, OH 43210
e-mail: tjs@stat.osu.edu

1Corresponding author.

Manuscript received October 11, 2013; final manuscript received March 31, 2014; accepted manuscript posted April 28, 2014; published online May 15, 2014. Assoc. Editor: Tammy L. Haut Donahue.

J Biomech Eng 136(7), 071007 (May 15, 2014) (8 pages) Paper No: BIO-13-1483; doi: 10.1115/1.4027510 History: Received October 11, 2013; Revised March 31, 2014; Accepted April 28, 2014

This paper describes a methodology for selecting a set of biomechanical engineering design variables to optimize the performance of an engineered meniscal substitute when implanted in a population of subjects whose characteristics can be specified stochastically. For the meniscal design problem where engineering variables include aspects of meniscal geometry and meniscal material properties, this method shows that meniscal designs having simultaneously large radial modulus and large circumferential modulus provide both low mean peak contact stress and small variability in peak contact stress when used in the specified subject population. The method also shows that the mean peak contact stress is relatively insensitive to meniscal permeability, so the permeability used in the manufacture of a meniscal substitute can be selected on the basis of manufacturing ease or cost. This is a multiple objective problem with the mean peak contact stress over the population of subjects and its variability both desired to be small. The problem is solved by using a predictor of the mean peak contact stress across the tibial plateau that was developed from experimentally measured peak contact stresses from two modalities. The first experimental modality provided computed peak contact stresses using a finite element computational simulator of the dynamic tibial contact stress during axial dynamic loading. A small number of meniscal designs with specified subject environmental inputs were selected to make computational runs and to provide training data for the predictor developed below. The second experimental modality consisted of measured peak contact stress from a set of cadaver knees. The cadaver measurements were used to bias-correct and calibrate the simulator output. Because the finite element simulator is expensive to evaluate, a rapidly computable (calibrated) Kriging predictor was used to explore extensively the contact stresses for a wide range of meniscal engineering inputs and subject variables. The predicted values were used to determine the Pareto optimal set of engineering inputs to minimize peak contact stresses in the targeted population of subjects.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Makris, E. A., Hadidi, P., and Athanasiou, K. A., 2011, “The Knee Meniscus: Structure-Function, Pathophysiology, Current Repair Techniques, and Prospects for Regeneration,” Biomaterials, 32, pp. 7411–7431. [CrossRef] [PubMed]
Rodkey, W. G., 2000, “Basic Biology of the Meniscus and Response to Injury,” Instrum. Course Lect., 49, pp. 189–193.
Fairbank, T. J., 1948, “Knee Joint Changes After Meniscectomy,” J. Bone Jt. Surg. Br. Vol., 30, pp. 664–670.
Allen, P. R., Denham, R. A., and Swan, A. V., 1984, “Late Degenerative Changes After Meniscectomy,” J. Bone Jt. Surg. Br. Vol., 66, pp. 666–671.
Bloecker, K., Englund, M., Wirth, W., Hudelmaier, M., Burgkart, R., Frobell, R., and Eckstein, F., 2011, “Size and Position of the Healthy Meniscus, and Its Correlation With Sex, Height, Weight, and Bone Area—A Cross-Sectional Study,” BMC Musculoskeletal Disord., 12(1), p. 248. [CrossRef]
Haut Donahue, T. L., Hull, M. L., Rashid, M. M., and Jacobs, C. R., 2004, “The Sensitivity of Tibiofemoral Contact Pressure to the Size and Shape of the Lateral and Medial Menisci,” J. Orthop. Res., 22(4), pp. 807–814. [CrossRef] [PubMed]
Haut, T. L., Hull, M. L., and Howell, S. M., 2000, “Use of Roentgenography and Magnetic Resonance Imaging to Predict Meniscal Geometry Determined With a Three-Dimensional Coordinate Digitizing System,” J. Orthop. Res., 18(2), pp. 228–237. [CrossRef] [PubMed]
Tissakht, M., and Ahmed, A. M., 1995, “Tensile Stress-Strain Characteristics of the Human Meniscal Material,” J. Biomech., 28(4), pp. 411–422. [CrossRef] [PubMed]
Fithian, D. C., Kelly, M. A., and Mow, V. C., 1990, “Material Properties and Structure-Function Relationships in the Menisci,” Clin. Orthop. Relat. Res., 252, pp. 19–31. [CrossRef] [PubMed]
Sweigart, M. A., Zhu, C. F., Burt, D. M., DeHoll, P. D., Agrawal, C. M., Clanton, T. O., and Athanasiou, K. A., 2004, “Intraspecies and Interspecies Comparison of the Compressive Properties of the Medial Meniscus,” Ann. Biomed. Eng., 32(11), pp. 1569–1579. [CrossRef] [PubMed]
Li, G., Park, S. E., DeFrate, L. E., Schutzer, M. E., Ji, L., Gill, T. J., and Rubash, H. E., 2005, “The Cartilage Thickness Distribution in the Tibiofemoral Joint and Its Correlation With Cartilage-to-Cartilage Contact,” Clin. Biomech.20, pp. 736–744. [CrossRef]
Eckstein, F., Yang, M., Guermazi, A., Roemer, F., Hudelmaier, M., Picha, K., Baribaud, F., Wirth, W., and Felson, D., 2010, “Reference Values and Z-Scores for Subregional Femorotibial Cartilage Thickness Results From a Large Population-Based Sample (Framingham) and Comparison With the Non-Exposed Osteoarthritis Initiative Reference Cohort,” Osteoarthritis Cartilage, 18(10), pp. 1275–1283. [CrossRef] [PubMed]
Akizuki, S., Mow, V. C., Müller, F., Pita, J. C., Howell, D. S., and Manicourt, D. H., 1986, “Tensile Properties of Human Knee Joint Cartilage: I. Influence of Ionic Conditions, Weight Bearing, and Fibrillation on the Tensile Modulus,” J. Orthop. Res., 4(4), pp. 379–392. [CrossRef] [PubMed]
Mow, V. C., Gu, W. Y., and Chen, F. H., 2005, “Structure and Function of Articular Cartilage and Meniscus,” Basic Orthopaedic Biomechanics and Mechano-Biology, V. C.Mow and R.Huiskes, eds., Lippincott, Williams, and Wilkins, Philadelphia, pp. 181–258.
Armstrong, C. G., and Mow, V. C., 1982, “Variations in the Intrinsic Mechanical Properties of Human Articular Cartilage With Age, Degeneration, and Water Content,” J. Bone Joint Surg. Am., 64(1), pp. 88–94. [PubMed]
Guo, H., Maher, S. A., and Spilker, R. L., 2013, “Biphasic Finite Element Contact Analysis of the Knee Joint Using an Augmented Lagrangian Method,” Med. Eng. Phys., 35(9), pp. 1313–1320. [CrossRef] [PubMed]
Donzelli, P. S., 1995, “A Mixed-Penalty Contact Finite Element Formulation for Biphasic Soft Tissue,” Ph.D. thesis, Rensselaer Polytechnic Institute, Troy, NY.
International Organization for Standardization, 2009, “Implants for Surgery—Wear of Total Knee-Joint Prostheses—Loading and Displacement Parameters for Wear-Testing Machines With Load Control and Corresponding Environmental Conditions for Test,” International Organization for Standardization, Geneva, ISO 14243-1:2002(E).
Guo, H., and Spilker, R. L., 2011, “Biphasic Finite Element Modeling of Hydrated Soft Tissue Contact Using an Augmented Lagrangian Method,” ASME J. Biomech. Eng., 133(11), p. 111001. [CrossRef]
Guo, H., Nickel, J. C., Iwasaki, L. R., and Spilker, R. L., 2012, “An Augmented Lagrangian Method for Sliding Contact of Soft Tissue,” ASME J. Biomech. Eng., 134(8), p. 084503. [CrossRef]
Guo, H., and Spilker, R. L., 2014, “An Augmented Lagrangian Finite Element Formulation for 3D Contact of Biphasic Tissues,” Comput. Methods Biomech. Biomed. Eng., 17(11), pp. 1206–1216. [CrossRef]
Guo, H., Shah, M., and Spilker, R. L., 2014, “A Finite Element Implementation for Biphasic Contact of Hydrated Porous Media Under Finite Deformation and Sliding,” Proc. Inst. Mech. Eng., Part H, 228(3), pp. 225–236. [CrossRef]
Higdon, D., Kennedy, M., Cavendish, J., Cafeo, J., and Ryne, R., 2004, “Combining Field Data and Computer Simulations for Calibration and Prediction,” SIAM J. Sci. Comput., 26, pp. 448–466. [CrossRef]
Saltelli, A., Chan, K., and Scott, E., 2000, Sensitivity Analysis, John Wiley & Sons, Chichester.
Santner, T. J., Williams, B. J., and Notz, W. I., 2003, The Design and Analysis of Computer Experiments, Springer Verlag, New York.
Gattiker, J. R., 2008, “Gaussian Process Models for Simulation Analysis (GPM/SA) Command, Function, and Data Structure Reference,” Los Alamos National Laboratory, Technical Report LA-UR-08-08057.
Higdon, D., Gattiker, J., Williams, B., and Rightley, M., 2008, “Computer Model Calibration Using High Dimensional Output,” J. Am. Stat. Assoc., 103, pp. 570–583. [CrossRef]
Draguljić, D., Santner, T. J., and Dean, A. M., 2012, “Non-Collapsing Spacing-Filling Designs for Bounded Polygonal Regions,” Technometrics, 54, pp. 169–178. [CrossRef]
Efron, B., and Tibshirani, R., 1997, “Improvements on Cross-Validation: The .632+ Bootstrap Method,” J. Am. Stat. Assoc., 92, pp. 548–560. [CrossRef]
Coello Coello, C. A., Lamont, G. B., and Van Veldhuizen, D. A., 2006, Evolutionary Algorithms for Solving Multi-Objective Problems (Genetic and Evolutionary Computation), Springer-Verlag, New York, Inc., Secaucus, NJ.
Forrester, A., Sóbester, A., and Keane, A., 2007, “Multi-Fidelity Optimization Via Surrogate Modeling,” Proc. R. Soc. A, 463(2088), pp. 3251–3269. [CrossRef]
Audet, C., Savard, G., and Zghal, W., 2010, “A Mesh Adaptive Direct Search Algorithm for Multiobjective Optimization,” Eur. J. Oper. Res., 204, pp. 545–556. [CrossRef]
Meakin, J. R., Shrive, N. G., Frank, C. B., and Hart, D. A., 2003, “Finite Element Analysis of the Meniscus: The Influence of Geometry and Material Properties on Its Behaviour,” Knee, 10(1), pp. 33–41. [CrossRef] [PubMed]
Haut Donahue, T. L., 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(1), pp. 19–34. [CrossRef] [PubMed]
Haut Donahue, T. L., Hull, M. L., and Howell, S. M., 2006, “New Algorithm for Selecting Meniscal Allografts that Best Match the Size and Shape of the Damaged Meniscus,” J. Orthop. Res., 24(7), pp. 1535–1543. [CrossRef] [PubMed]
Mow, V. C., Kuei, S. C., Lai, W. M., and Armstrong, C. G., 1980, “Biphasic Creep and Stress Relaxation of Articular Cartilage in Compression, Theory, and Experiments,” ASME J. Biomech. Eng., 102(1), pp. 73–84. [CrossRef]
Kelly, B. T., Robertson, W., Potter, H. G., Deng, X. H., Turner, A. S., Lyman, S., Warren, R. F., and Rodeo, S. A., 2007, “Hydrogel Meniscal Replacement in the Sheep Knee: Preliminary Evaluation of Chondroprotective Effects,” Am. J. Sports Med., 35, pp. 43–52. [CrossRef] [PubMed]
Messner, K., and Gillquist, J., 1993, “Prosthetic Replacement of the Rabbit Medial Meniscus,” J. Biomed. Mater. Res., 27, pp. 1165–1173. [CrossRef] [PubMed]
Setton, L. A., Guilak, F., Hsu, E. W., and Vail, T. P., 1999, “Biomechanical Factors in Tissue Engineered Meniscal Repair,” Clin. Orthop. Relat. Res., 367S, pp. 254–272. [CrossRef]
Hutchinson, I. D., Moran, C. J., Potter, H. G., Warren, R. F., and Rodeo, S. A., 2014, “Restoration of the Meniscus: Form and Function,” Amer. J. Sports Med., 42(4), pp. 987–998. [CrossRef]
Alhalki, M. M., Howell, S. M., and Hull, M. L., 1999, “How Three Methods for Fixing a Medial Meniscal Autograft Affect Tibial Contact Mechanics,” Am. J. Sports Med., 27, pp. 320–328. [PubMed]
Dienst, M., Greis, P. E., Ellis, B. J., Bachus, K. N., and Burks, R. T., 2007, “Effect of Lateral Meniscal Allograft Sizing on Contact Mechanics of the Lateral Tibial Plateau: An Experimental Study in Human Cadaveric Knee Joints,” Am. J. Sports Med., 35(1), pp. 34–42. [CrossRef] [PubMed]
Rodeo, S. A., 2001, “Meniscal Allografts—Where Do We Stand,” Am. J. Sports Med., 29(2), pp. 246–261. [PubMed]
Karvonen, R. L., Negendank, W. G., Teitge, R. A., Reed, A. H., Miller, P. R., and Fernandez-Madrid, F., 1994, “Factors Affecting Articular Cartilage Thickness in Osteoarthritis and Aging,” J. Rheumatol., 21(7), pp. 1310–1318. [PubMed]
Holmes, M., and Mow, V., 1990, “The Nonlinear Characteristics of Soft Gels and Hydrated Connective Tissues in Ultrafiltration,” J. Biomech., 23, pp. 1145–1156. [CrossRef] [PubMed]
Gilbert, S., Chen, T., Hutchinson, I. D., Choi, D., Voigt, C., Warren, R. F., and Maher, S. A., 2013, “Dynamic Contact Mechanics on the Tibial Plateau of the Human Knee During Activities of Daily Living,” J. Biomech., (in press). [CrossRef]


Grahic Jump Location
Fig. 1

A schematic diagram of the 2D axisymmetric knee joint

Grahic Jump Location
Fig. 2

ISO 2009 dynamic axial force during gait. This profile was used in runs of the bFEM; the profile, scaled to 10 s, was used in the cadaver model.

Grahic Jump Location
Fig. 3

The smoothed surface connecting the discrete set of bFEM-determined contact stresses across radial position and time when run at meniscal design inputs (hm, hc, Erm, Ecm, km) = (6.128, 2.980, 10.376, 108.800, 2.202) and patient inputs (ht, hf, Ec, kc) = (2.073, 2.044, 6.709, 5.289). The radial position of 0 mm corresponds to the center of the lateral compartment of the knee.

Grahic Jump Location
Fig. 4

(a) The tibial component with attached sensor; note the box depicting the line of interest. (b) The experimental setup: during loading the tibia was aligned with the femoral component which was fixed to the crosshead of the MTS. The load profile in Fig. 2, scaled to 10 s, was applied to represent the axial force profile that occurs during gait.

Grahic Jump Location
Fig. 5

(a) Contact stress data from a representative, experimentally tested, specimen. A line of interest was identified and the region covered by the meniscus was identified. (Each square represents a 2 × 2 mm sensing element.) (b) The contact stress, as a function of radial distance from the edge of the tibial spine, was plotted for both sensors and for the computational model; the data shown is from a representative knee at 14% of the loading cycle. The vertical line near 12 mm in (b) marks the inner edge of the meniscus.

Grahic Jump Location
Fig. 6

ME plots for inputs hm (left), ht (center), and Ec (right) on the peak contact stress predictions at 14% of the axial loading profile

Grahic Jump Location
Fig. 7

The distribution of predicted peak contact stress at 14% (left) and 45% (right) of the axial loading profile for 10,000 different patient conditions, for the meniscal design (hm, hc, Erm, Ecm, km) = (4.513, 1.698, 9.352, 83.065, 1.750)

Grahic Jump Location
Fig. 8

Means and standard deviations of peak contact stress predictions at 14% of the axial loading profile for the 10,000 different patient conditions corresponding to 50 fixed meniscal designs. Desirable meniscal designs (those with low mean and low standard deviation) are marked with a filled circle.

Grahic Jump Location
Fig. 9

Distributions of the peak contact stresses at 14% of the axial loading profile over the 10,000 draws of patient cartilage values for the four optimal meniscal designs listed in Table 3. The filled circles mark the 95th-percentile of the peak contact stresses, and the open circles mark the 99th-percentile of these stresses.

Grahic Jump Location
Fig. 10

2D projections of 50 fixed meniscal designs corresponding to the designs used to construct Fig. 8. Desirable meniscal designs (those corresponding to low mean and low standard deviation) are marked with a filled circle.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In