The fidelity of a 3D model created using image segmentation must be precisely quantified and evaluated for the model to be trusted for use in subsequent biomechanical studies such as finite element analysis. The bones within the ankle joint vary significantly in size and shape. The purpose of this study was to test the hypothesis that the accuracy and reliability of a segmented bone geometry is independent of the particular bone being measured. Computed tomography (CT) scan data (slice thickness 1 mm, pixel size 808±7 µm) from three anonymous patients was used for the development of the ankle geometries (consisting of the tibia, fibula, talus, calcaneus, and navicular bones) using Simpleware Scan IP software (Synopsys, Exeter, UK). Each CT scan was segmented 4 times by an inexperienced undergraduate, resulting in a total of 12 geometry assemblies. An experienced researcher segmented each scan once, and this was used as the ‘gold standard’ to quantify the accuracy. The solid bone geometries were imported into CAD software (Inventor 2023, Autodesk, CA, USA) for measurement of the surface area and volume of each bone, and the distances between bones (tibia to talus, talus to navicular, talus to calcaneus, and tibia to fibula) were carried out. The intra-class coefficient (ICC) was used to assess intra-observer reliability. Bland Altman plots were employed as a statistical measure for criteria validity (accuracy) [1].Abstract
Objectives
Methods
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. 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.Abstract
Objectives
Methods
The aim of this study was to develop an open-source finite element model of the ankle for identification of the best clinical treatment to restore stability to the ankle after injury. The ankle geometry was defined from the Visible Human Project Female CT dataset available from the National Library of Medicine, and segmented using Dragonfly software (Object Research Systems, 2020). The finite element model was created with FEBio (University of Utah, 2021) using the dynamic nonlinear implicit solver. Linear isotropic material properties were assigned to the bones (E=7300MPa, ν=0.3, ρ=1730kg/m3) and cartilage (E=10MPa, ν=0.4, ρ=1100kg/m3). Spring elements were used to represent the ligaments and material properties were taken from Mondal et al. [1]. Lagrangian contact was defined between the cartilaginous surfaces with μ=0.003. A standing load case was modelled, assuming even distribution of load between the feet. A reaction force of 344.3N was applied to the base of the foot, a muscle force of 252.2N, and the proximal ends of the tibia and fibula were fully constrained.Abstract
Objectives
Methods
