Research Papers

Finite Element Modeling of the Contact Geometry and Deformation in Biomechanics Applications1

[+] Author and Article Information
Ahmed A. Shabana

Department of Mechanical and
Industrial Engineering,
University of Illinois at Chicago,
Chicago, IL 60607

A version of this paper was presented at the 2012 ASME Design Engineering Technical Conferences and Computer and Information in Engineering Conference.

Contributed by the Design Engineering Division of ASME for publication in the Journal of Computational and Nonlinear Dynamics. Manuscript received August 21, 2012; final manuscript received April 26, 2013; published online June 10, 2013. Assoc. Editor: Aki Mikkola.

J. Comput. Nonlinear Dynam 8(4), 041013 (Jun 10, 2013) (11 pages) Paper No: CND-12-1127; doi: 10.1115/1.4024541 History: Received August 21, 2012; Revised April 26, 2013

The main contribution of this paper is to demonstrate the feasibility of using one computational environment for developing accurate geometry as well as performing the analysis of detailed biomechanics models. To this end, the finite element (FE) absolute nodal coordinate formulation (ANCF) and multibody system (MBS) algorithms are used in modeling both the contact geometry and ligaments deformations in biomechanics applications. Two ANCF approaches can be used to model the rigid contact surface geometry. In the first approach, fully parameterized ANCF volume elements are converted to surface geometry using parametric relationship that reduces the number of independent coordinate lines. This parametric relationship can be defined analytically or using a spline function representation. In the second approach, an ANCF surface that defines a gradient deficient thin plate element is used. This second approach does not require the use of parametric relations or spline function representations. These two geometric approaches shed light on the generality of and the flexibility offered by the ANCF geometry as compared to computational geometry (CG) methods such as B-splines and NURBS (Non-Uniform Rational B-Splines). Furthermore, because B-spline and NURBS representations employ a rigid recurrence structure, they are not suited as general analysis tools that capture different types of joint discontinuities. ANCF finite elements, on the other hand, lend themselves easily to geometric description and can additionally be used effectively in the analysis of ligaments, muscles, and soft tissues (LMST), as demonstrated in this paper using the knee joint as an example. In this study, ANCF finite elements are used to define the femur/tibia rigid body contact surface geometry. The same ANCF finite elements are also used to model the MCL and LCL ligament deformations. Two different contact formulations are used in this investigation to predict the femur/tibia contact forces; the elastic contact formulation which allows for penetrations and separations at the contact points, and the constraint contact formulation in which the nonconformal contact conditions are imposed as constraint equations, and as a consequence, no separations or penetrations at the contact points are allowed. For both formulations, the contact surfaces are described in a parametric form using surface parameters that enter into the ANCF finite element geometric description. A set of nonlinear algebraic equations that depend on the surface parameters is developed and used to determine the location of the contact points. These two contact formulations are implemented in a general MBS algorithm that allows for modeling rigid and flexible body dynamics.

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


Machado, M., Flores, P., Pimenta Claro, J. C., Ambrosio, J., Silva, M., and Completo, A., 2010, “Development of a Planar Multibody Model of the Human Knee Joint,” Nonlinear Dyn., 60(3), pp. 459–478. [CrossRef]
Blankevoort, L., and Huiskes, R., 1996, “Validation of a Three-Dimensional Model of the Knee,” J. Biomech., 29, pp. 955–961. [CrossRef] [PubMed]
Bendjaballah, M. Z., Shirazi-adl, A., and Zukor, D. J., 1995, “Biomechanics of the Human Knee Joint in Compression: Reconstruction, Mesh Generation and Finite Element Analysis,” The Knee, 2(2), pp. 69–79. [CrossRef]
Donahue, T. L. H., Hull, M. L., Rashid, M. M., and Jacobs, C. R., 2002, “A Finite Element Model of the Human Knee Joint for the Study of Tibio-Femoral Contact,” ASME J. Biomech. Eng., 124, pp. 273–280. [CrossRef]
McLean, S. G., Su, A., and Van den Bogert, A. J., 2003, “Development and Validation of a 3-D Model to Predict Knee Joint Loading during Dynamic Movement,” ASME J. Biomech. Eng., 125, pp. 864–874. [CrossRef]
Bei, Y., and Fregly, B., 2004, “Multibody Dynamic Simulation of the Knee Contact Mechanics,” J. Med. Eng. Phys., 26, pp. 777–789. [CrossRef]
Blankevoort, L., Kuiper, J. H., Huiskes, R., and Grootenboer, H. J., 1991, “Articular Contact in a Three-Dimensional Model of the Knee,” J. Biomech., 24, pp. 1019–1031. [CrossRef] [PubMed]
Abdel-Rahman, E. M., and Hefzy, M. S., 1998, “Biomechanics of the Human Knee Joint in Compression: Reconstruction, Mesh Generation and Finite Element Analysis Three-Dimensional Dynamic Behavior of the Human Knee Joint under Impact Loading,” Med. Eng. Phys., 20, pp. 276–290. [CrossRef] [PubMed]
Elias, J. J., Wilson, D. R., Adamson, R., and Cosgarea, A. J., 2004, “Evaluation of a Computational Model Used to predict the Patella Femoral Contact Pressure Distribution,” J. Biomech., 37, pp. 295–302. [CrossRef] [PubMed]
Piazza, S. J., and Delp, S. L., 2001, “Three-Dimensional Simulation of Total Knee Replacement Motion During a Step-Up Task,” ASME J. Biomech. Eng., 123, pp. 599–606. [CrossRef]
Engel, K., Herpers, R., and Hartmann, U., 2011, Biomechanical Computer Models, Theoretical Biomechanics, Vaclav Klika, ed., In Tech, New York.
Guess, T. M., Thiagarajan, G., Kia, M., and Mishra, M., 2010, “A Subject Specific Multibody Model of the Knee With Menisci,” Med. Eng. Phys., 32(5), pp. 505–515. [CrossRef] [PubMed]
Gantoi, F. M., Brown, M. A., and Shabana, A. A., 2010, “ANCF Finite Element/Multibody System Formulation of the Ligament/Bone Insertion Site Constraints,” ASME J Comput. Nonlinear Dyn., 5(3), p. 031006. [CrossRef]
Bartel, D. L., Davy, D. T., and Keaveny, T. M., 2006, Orthopaedic Biomechanics: Mechanics and Design in Musculoskeletal Systems, Pearson, Delhi.
Shabana, A. A., Zaazaa, K. E., and Sugiyama, H., 2008, Railroad Vehicle Dynamics: A Computational Approach, CRC Press, Boca Raton, FL.
Cheng, R. C. K., Brown, T. D., and Andrews, J. G., 1990, “Non-Uniqueness of the Bicompartmental Contact Force Solution in a Lumped-Parameter Mathematical Model of the Knee,” J. Biomech., 23, pp. 353–355. [CrossRef] [PubMed]
Mohamed, A. A., Brown, M., and Shabana, A. A., 2010, “Study of the Ligament Tension and Cross-Section Deformation using Nonlinear Finite Element/Multibody System Algorithms,” Multibody Syst. Dyn., 23, pp. 227–248. [CrossRef]
Weed, D., Maqueda, L., Brown, M., and Shabana, A., 2008, “A Multibody/Finite Element Nonlinear Formulation of a Two-Ligament Knee Joint,” Proceedings of the 2008 ASME International Mechanical Engineering Congress and Exhibition, Boston, MA, Oct. 31–Nov. 6, pp. 407–416.
Ashraf, T., Beard, D. J., and Newman, J. H., 2003, “Symmetrical vs. Asymmetrical Total Knee Replacement: A Medium Term Comparative Analysis,” The Knee, 10(1), pp. 61–66. [CrossRef] [PubMed]
Shabana, A. A., 2012, Computational Continuum Mechanics, 2nd ed., Cambridge University Press, Cambridge, UK.
Sanborn, G. G., and Shabana, A. A., 2009, “On the Integration of the Computer Aided Design and Analysis using the Finite Element Absolute Nodal Coordinate Formulation,” J. Multibody Syst. Dyn., 22, pp. 181–197. [CrossRef]
Mikkola, A., and Shabana, A. A., 2012, “Comparison Between ANCF and B-Spline Surfaces,” Proceedings of the Second Joint International Conference on Multibody System Dynamics, Stuttgart, Germany, 29 May 29–June 1.
Nordin, M., and Frankel, V. H., 2001, Basic Biomechanics of the Musculoskeletal System, 3rd ed., Lippincott Williams & Willkins, Philadelphia, PA.
Kapandji, I. A., 1970, The Physiology of the Joints, Churchill Livingstone, Edinburgh, UK.
Leondes, C. T., 2007, Biomechanical Systems Technology, Vol. 1: Computational Methods, World Scientific Publishing Co. Pte. Ltd., New York, USA.
Mow, V. C., and Hayes, W. C., 1997, Basic Orthopedic Biomechanics, Lippincott-Raven Publishers, Philadelphia, PA.
Milner, C. E., 2008, Functional Anatomy for Sport and Exercise, Routledge, New York.
Kennedy, J. C., Weinberg, H. W., and Wilson, A. S., 1974, “The Anatomy and Function of the Anterior Cruciate Ligament,” J. Bone Jt. Surg., 56A(2), pp. 223–235.
Amis, A. A., Bull, A. M. J., Gupte, C. M., Hijazi, I., and Robinson, J. R., 2003, “Biomechanics of the PCL and Related Structures: Postero-Lateral, Postero-Medial and Menisco-Femoral Ligaments,” Knee Surg. Sports Traumatol. Arthrosc., 11, pp. 271–281. [CrossRef] [PubMed]
Wang, C. J., 2002, “Injuries to the Posterior Cruciate Ligament and Posterolateral Instabilities of the Knee,” Chang Gung Med. J., 25, pp. 288–297. [PubMed]
Peña, E., Calvo, B., Martínez, M. A., and Doblaré, M., 2006, “A Three-Dimensional Finite Element Analysis of the Combined Behavior of Ligaments and Menisci in the Healthy Human Knee Joint,” J. Biomech., 39(9), pp. 1686–1701. [CrossRef] [PubMed]
Vollebregt, E. A. H., 2008, “Survey of Programs on Contact Mechanics Developed by J. J. Kalker,” Veh. Syst. Dyn., 46(1), pp. 85–92. [CrossRef]
Fung, Y. C., 1993, Biomechanics. Mechanical Properties of Living Tissues, 2nd ed., Springer-Verlag, New York.
Shabana, A. A., 2005, Dynamics of Multibody Systems, 3rd ed., Cambridge University Press, Cambridge, UK.
İşcan, M. Y., 2005, “Forensic Anthropology of Sex and Body Size, Forensic Science International,” Forensic Sci. Int., 147(2–3), pp. 107–112. [CrossRef]
Mommersteeg, T. J. A., Blankevoort, L., Huiskes, R., Kooloos, J. G. M., and HendriksJ. C. M., 1994, “The Effect of Variable Relative Insertion Orientation of Human Knee Bone Ligament Bone Complexes on the Tensile Stiffness,” J. Biomech., 28(6), pp. 745–752. [CrossRef]
Song, Y., Debski, R. E., Musahl, V., Thomas, M., and Woo, S. L., 2004, “A Three-Dimensional Finite Element Model of the Human Anterior Cruciate Ligament: A Computational Analysis with Experimental Validation,” J. Biomech., 37, pp. 383–390. [CrossRef] [PubMed]
Cheng, C. K., 1988, “A Mathematical Model for Predicting Bony Contact Forces and Muscles Forces at the Knee During the Human Gait,” Ph.D. thesis, The University of Iowa, Iowa City, IA.
Sasaki, K., and Neptune, R. R., 2010, “Individual Muscle Contributions to the Axial Knee Joint Contact Force during Normal Walking,” J. Biomech., 43, pp. 2780–2784. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 2

Femur and tibia profile

Grahic Jump Location
Fig. 3

Surface parameterization

Grahic Jump Location
Fig. 4

ANCF 3D beam element

Grahic Jump Location
Fig. 5

X-displacement of the femur's center of mass ( ECF, ECCF)

Grahic Jump Location
Fig. 6

Z-displacement of the femur's center of mass ( ECF, ECCF)

Grahic Jump Location
Fig. 7

Normal contact force for the knee lateral side ( ECF, ECCF)

Grahic Jump Location
Fig. 8

Normal contact force for the knee medial side ( ECF, ECCF)

Grahic Jump Location
Fig. 9

Sum of normal contact forces (lateral side, medial side, sum)

Grahic Jump Location
Fig. 10

Deformation of the ACL ligament ( ECF, ECCF)

Grahic Jump Location
Fig. 11

Deformation of the PCL ligament ( ECF, ECCF)

Grahic Jump Location
Fig. 12

MCL cross section deformation at the midpoint ( ECF, ECCF)

Grahic Jump Location
Fig. 13

Axial strain at the midpoint for MCL ( ECF, ECCF)

Grahic Jump Location
Fig. 14

Axial strain at the midpoint for MCL (—— previous knee model, actual knee model)

Grahic Jump Location
Fig. 15

Normal contact force for the knee lateral side (—— ECF, ECCF)

Grahic Jump Location
Fig. 16

Normal contact force for the knee medial side (—— ECF, ECCF)



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