On the development of an osseo-ligamentous finite element model of the human ankle joint
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)
- et al.
Response of joints to impact loading-III
Journal of Biomechanics
(1973) - et al.
Wave attenuation in skeleton of young healthy persons
Journal of Biomechanics
(1981) Grants’s Atlas of Anatomy
(1983)- 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...
Reduction of calcaneal fractures by the McReynolds medial approach technique and its experimental basis
Clinical Orthopedics
(1983)- 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,...
Physical Properties of Tissue: A Comprehensive Reference Book
A survey of finite element analysis in orthopedics biomechanics: the first decade
Journal of Biomechanics
A review of biomechanical models
Journal of Biomechanical Engineering
Movements of the subtalar and transverse tarsal joints
Anatomical Records
Cited by (62)
Clinically useful finite element models of the natural ankle – A review
2023, Clinical BiomechanicsFinite Element Modeling
2022, Foot and Ankle BiomechanicsDevelopment of customized finite element models of medial column fixation using an intramedullary beam: A computational sensitivity analysis
2021, Medical Engineering and PhysicsFinite element analysis of the effect of sagittal angle on ankle joint stability in posterior malleolus fracture: A cohort study
2019, International Journal of SurgeryCitation 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, FootCitation 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.