DEVELOPMENT OF AN OPEN-SOURCE FINITE ELEMENT ANKLE MODEL FOR IN-SILICO INVESTIGATION OF THE ANKLE JOINT
Authors: Mondal Subrata , Mangwani Jitendra , Brockett Claire , Gulati Aashish , Pegg Elise 
Published: November 2023

Objectives

This abstract provides an update on the Open Ankle Models being developed at the University of Bath. The goal of this project is to create three fully open-source finite element (FE) ankle models, including bones, ligaments, and cartilages, appropriate musculoskeletal loading and boundary conditions, and heterogeneous material property distribution for a standardised representation of ankle biomechanics and pre-clinical ankle joint analysis.

Methods

A computed tomography (CT) scan data (pixel size of 0.815 mm, and slice thickness of 1 mm) was used to develop the 3D geometry of the bones (tibia, talus, calcaneus, fibula, and navicular). Each bone was given the properties of a heterogeneous elastic material based on the CT greyscale. The density values for each bone element were calculated using a linear empirical relation, ρ= 0.0405 + (0.000918) HU and then power law equations were utilised to get the Young's Modulus value for each bone element [1]. At the bone junction, a thickness of cartilage ranging from 0.5–1 mm, and was modelled as a linear material (E=10 MPa, ν=0.4 [2]). All ligament insertions and positions were represented by four parallel spring elements, and the ligament stiffness and material attributes were applied in accordance with the published literature [2]. The ankle model was subjected to static loading (balance standing position). Four noded tetrahedral elements were used for the discretization of bones and cartilages. All degrees of freedom were restricted at the proximal ends of the tibia and fibula. The ground reaction forces were applied at the underneath of the calcaneus bone. The interaction between the cartilages and bones was modelled using an augmented contact algorithm with a sliding elastic contact between each cartilage. A tied elastic contact was used between the cartilages and the bone. FEbio 2.1.0 (University of Utah, USA) was used to construct the open-source ankle model.

Results

When the double-legged stance phase loading condition was taken into consideration, stress at the antero-medial tibial wall (ranged from 1 to 7 MPa) was found to be similar to the prior work [2], indicating bulk of the load transfer was through this region. The maximum principal strain was predicted at the different regions on bones around the ankle joint. The proximal surface of the talus, and tibial distal surface were shown to have the highest maximum principal strains followed by antero-medial walls of the tibia bone, at the proximal location.

Conclusions

The present open 3D FE model of the ankle will assist researchers in better understanding ankle biomechanics, precisely predicting load transfer, and examining the ankle to address unmet clinical needs for this joint. The results of the current investigation are realistic in terms of load transfer and stress-strain distribution across the ankle joint and well comparable to those reported in the literature [2]. However, sensitivity and ankle instability simulations will be performed in future work to investigate the model's reliability and robustness.

Declaration of Interest

(b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project.