On the development of an osseo-ligamentous finite element model of the human ankle joint

https://doi.org/10.1016/S0020-7683(00)00129-3Get rights and content

Abstract

A three-dimensional finite element model of the human ankle joint has been developed to study the mechanisms of impact injury to the major bones of the foot. The model is based on anatomically realistic bone geometry obtained from medical imaging and includes the major ligaments of the ankle. Careful consideration was given to model the soft tissues of the plantar surface of the foot. The model was used to simulate axial impulsive loading applied at the plantar surface of the foot. The stability of the ankle joint was achieved in the model strictly by the intrinsic anatomical geometry and the ligamenteous structure. The time history response of input/output accelerations and forces compared reasonably well with experimental data. Results indicate that the calcaneus experiences the highest stresses followed by the tibia and talus. This is in agreement with several experimental data on calcaneal fracture in axial dynamic loading. Also, the model gives stress localization in the lateral–collateral ligaments that agrees with injury observations for that region. The significance of the model lies in its potential uses as a research tool for understanding the mechanical response of the ankle related to injury and degenerative disease.

Introduction

Lower extremity trauma is common among athletes involved in high-energy sports, and among car crash occupants. According to Thomas et al. (1995), there are more than 100,000 ankle fractures that occur each year in the United States. Full recovery from these injuries is not always possible, resulting in great disruptions in the lives of the victims. They come at a great cost to society due to temporary and permanent impairment and disability (Mackenzie, 1986). It is estimated that the medical cost of ankle fracture repair is around $10,000 on average (Miller et al., 1995). The understanding of the mechanical behaviour of the ankle joint under injurious situations can enhance the delivery of treatment and rehabilitation services. It can also play a major role in the prevention of injury.

Epidemiological, clinical and biomechanical studies have implicated axial impact to plantar surface of the foot as a main cause of lower extremity trauma States, 1986, Pattimore et al., 1991, Kallina et al., 1995. Studies conducted by Radin et al. (1973) and Simon et al. (1972) have indicated a relationship between excessive levels of axial loads generated from gait and the damage of the articular cartilage in the lower extremity. Wosk and Voloshin (1981) and Voloshin et al. (1981) studied the capability of the human shock absorber system to attenuate input dynamic loading during gait. Using accelerometers mounted to specified points at the knee, pelvis, and head, the authors determined the percentage of peak acceleration attenuation offered by each body region. They observed a 30% reduction of the peak of the impulses at the knee and 70% attenuation at the head. Voloshin and Wosk (1981) also showed that healthy subjects possess more efficient shock-absorbing capability than do subjects with joint disorders. Recently, Yoganandan et al., 1995, Yoganandan et al., 1996 evaluated the biomechanics of the human-ankle complex under axial impact. The applied dynamic axial loads to the heel of cadaveric specimens and observed fracture. Values of 15.1 kN (SD=2.7) and 10.2 kN (SD=1.5) were determined for the mean dynamic forces at the plantar surface of the foot and the proximal end of the tibia, respectively. Both of these values were of sufficient magnitude to cause ankle fracture.

Over the past 20 years, the finite element method (FE) has evolved into a well-established computational tool in biomechanics Huiskes and Chao, 1981, King, 1984. Once validated, a FE model can be used as a tool for parametric investigation of stress distribution in biological structures. Recently, three-dimensional FE models for investigating the response of lower extremity to dynamic loads have been developed and reported Beaugonin et al., 1995, Tannous et al., 1996. The rationale for this research is that a FE model of the human lower leg can facilitate the understanding of lower leg injury mechanisms and can be used in a variety of applications in the medical and safety arenas.

Section snippets

Model geometry

The three-dimensional FE model was created using medical imaging data for the 26 bones of the foot along with the tibia and fibula. patran/p75 (msc/patran, 1997) was utilized to construct the model. The human bone was constructed by taking into account of the two distinct bone types with different material properties. These are referred to as compact (or cortical) bone and trabecular (or cancellous) bone. Compact bone always surrounds the trabecular bone and its thickness varies from one

Results

Fig. 5, Fig. 6, Fig. 7, Fig. 8, Fig. 9 show a comparison of the kinematic response between the numerical simulations and the cadaveric tests. The solid lines indicate the simulation results while the test results are indicated by the dotted lines. The shape and pulse duration of input/output acceleration and force time histories between the model and experiments are in good agreement for low impact velocity cases (7.3 and 11 ft/s) with the rise times and slopes being nearly identical. This

Discussion and conclusion

A three-dimensional FE model of the human ankle has been developed and subjected to preliminary verifications against a limited set of experimental data. The model was shown to have reasonable overall mass and stiffness characteristics and is capable of reproducing the experimental input/output responses for a range of impact velocities. At impact velocities greater than 14.6 ft/s (4.47 m/s), however, the model gave higher force and acceleration peak amplitudes than experiment. This discrepancy

References (40)

  • E Radin et al.

    Response of joints to impact loading-III

    Journal of Biomechanics

    (1973)
  • J Wosk et al.

    Wave attenuation in skeleton of young healthy persons

    Journal of Biomechanics

    (1981)
  • J.E Anderson

    Grants’s Atlas of Anatomy

    (1983)
  • D.E Attarian et al.

    Biomechanics characteristics of the human ankle ligaments

    Foot and Ankle

    (1985)
  • Beaugonin, M., Haug, E., Munvk, G., Cesari, D., 1995. A preliminary numerical model of the human ankle under impact...
  • Begeman, P., Aekbote, K., 1996. Axial load strength and some ligament properties of the ankle joint. Proceedings of the...
  • Benett, M.B., Kerr, R.F., 1990. The mechanical properties of the human subcalcaneal fat pad in compression. Journal of...
  • B.D Burdeaux

    Reduction of calcaneal fractures by the McReynolds medial approach technique and its experimental basis

    Clinical Orthopedics

    (1983)
  • D.R Carter et al.

    The behavior of bone as a two-phase porous structure

    Journal of Bone Joint Surgery

    (1977)
  • Cowin, C.S., 1991. The mechanical properties of cortical bone tissue. In: Cowin, C.S. (Ed.), Bone Mechanics, CRC,...
  • F Duck

    Physical Properties of Tissue: A Comprehensive Reference Book

    (1990)
  • R Huiskes et al.

    A survey of finite element analysis in orthopedics biomechanics: the first decade

    Journal of Biomechanics

    (1981)
  • Kallina, I., Scheunert, D., Scheerer, J., Breitner, R., Zeidler, R.B., 1995. Injuries to the lower leg: significance...
  • A.I King

    A review of biomechanical models

    Journal of Biomechanical Engineering

    (1984)
  • Kitagawa, Y., Ichikawa, H., King, A.I., Levine, R.S., 1998. A severe ankle and foot injury in frontal crashes and its...
  • Mackenzie, E.J., 1986. Public health impact of lower extremity trauma. Biomechanics and Medical Aspects of Lower Limb...
  • J.T Manter

    Movements of the subtalar and transverse tarsal joints

    Anatomical Records

    (1941)
  • Matejczyk, M.B., Greenwald, A.S., Black, J.D., Wilde, A.H., 1978. Ankle joint mechanics and implant evaluation....
  • Miller, T.R., Martin, P.G., Crandall, J.R., 1995. Costs of lower limb injuries in highway crashes. Proceedings of the...
  • Mow, V.C., Proctor, C.S., Kelly, M.A., 1989. Biomechanics of articular cartilage. In: Nordin, M., Frankel, V. (Eds.),...
  • Cited by (62)

    • Finite Element Modeling

      2022, Foot and Ankle Biomechanics
    • Finite element analysis of the effect of sagittal angle on ankle joint stability in posterior malleolus fracture: A cohort study

      2019, International Journal of Surgery
      Citation Excerpt :

      The geometric model was imported into Hypermesh 2017. Meshing and material property distributions (Table 1) were carried out [8]. Finally, the finite element model was obtained.

    • Bone cyst formation after ankle arthroplasty may be caused by stress shielding. A numerical simulation of the strain adaptive bone remodelling

      2017, Foot
      Citation Excerpt :

      Already several studies validated the FEM in order of prosthesis and bone interaction. FE models on the foot have mainly been used for the development of orthopedic footwear, insoles and orthoses [12–19]. Bouguecha et al. [20] validated the numerical simulation of the the ankle joint.

    View all citing articles on Scopus
    View full text