Processing of the CT files and three-dimensional solid reconstruction
An image of the right ankle joint of an adult man in a neutral position was obtained by CT tomography (64 slices, SIEMENS, US), and the image was input into the 3D reconstruction software Mimics in the Dicom format to obtain a clear skeleton outline. After mask processing, the image was read in Geomagic (Geomagic,US) in the STL format, reverse engineering reconstruction was completed, and the 3D graphics in the IGES file format were generated (Fig. 1).
Construction of the working condition model
After the foot skeleton and foot contour were built, all ligaments were connected with lines in the physiological position, and a complete foot model was generated. According to the anatomical data of the joint surface, cartilage boundaries were established, and cartilage joints were built with Geomagic with an offset thickness of 1 mm. According to the requirements of the analysis, only the fibula and tibia of the calcaneal talus, as well as the related cartilages and ligaments, were required to be retained; thus, a relatively complete three-dimensional finite element model of the ankle joint of the normal adults was built. On the basis of the normal model, the talus cartilage was divided into 9 regions by the nine-grid partition method. Studies have shown that area 4 is the most common area for talus cartilage injuries [13, 14]. In this study, defects in the cartilage and subchondral bone in region 4 of the talus were simulated. Because the existing literature does not study the area size of finite element and talus injury, when we conduct finite element analysis, the experimental measurement depth 1mm, and the area defect size is 2 mm * 2 mm, 4 mm * 4 mm, 6 mm * 6 mm, 8 mm *8 mm, 10 mm *10 mm, 12 mm * 12 mm (Fig. 2).
Meshing
The assembled solid model was imported into Ansys Workbench (Ansys, USA), a Boolean operation was carried out, material parameters were assigned, contact was defined, and then, the grid division process was completed. The solid unit comprised Solid 187 and Solid 95, the ligament was a Link180 unit, and its nonlinear characteristics were set under tension without pressure (Fig. 3).
Material parameters and Contact Settings
All kinds of tissue materials involved in this model were simplified into isotropic homogeneous elastic materials, and the material parameters are listed in Table 1 and Table 2.
Table 1
Properties of the bone and cartilage materials
Material | Modulus of Elasticity (MPa) | Poisson's ratio |
Bone | 7300 | 0.3 |
Cartilage | 12 | 0.42 |
Table 2
Material properties of the ligaments
Ligament | Modulus of Elasticity (MPa) | Poisson's ratio | Sectional area (mm2) |
AtiF | 260 | 0.4 | 18.4 |
PtiF | 260 | 0.4 | 18.4 |
AtaFi | 255.5 | 0.4 | 12.9 |
PtaFi | 216.5 | 0.4 | 21.9 |
CaTi | 512 | 0.4 | 9.7 |
AtiTa | 184.5 | 0.4 | 13.5 |
PtiTa | 99.5 | 0.4 | 22.6 |
TiCa | 512 | 0.4 | 9.7 |
TiNa | 320.7 | 0.4 | 7.1 |
AtiF = Anterior tibiofibular ligament,PtiF = Posterior tibiofibular ligament༌AtaFi = Anterior talofibular ligament༌PtaFi = Posterior talofibular ligament༌CaTi = Calcaneofibular ligament༌AtiTa = Anterior tibial ligament༌PtiTa = Posterior tibial talus ligament༌TiCa = Tibiocalcanean ligament༌TiNa = Tibionavicular ligament |
The settings of the contact between the components were set according to the actual condition. The cartilage was bound to the corresponding bones, and the friction coefficient between the articular surfaces of the cartilage was 0.01.
Applying loads and constraints
The grid direction of the corresponding sites of the calcaneus and scaphoid was constrained so that the degree of freedom was 0. Three gait patterns were selected for analysis according previous studies, as shown below, and it was assumed that the body weight was 600 N and the foot length was 25.4 cm (Fig. 4). After the model was established, it was verified that it was close to those in previous studies. [8–9].
Experimental groups and data acquisition process
After the above model was established, 7 groups were established for the experiment: the normal talus osteochondral group and groups the area defect size is 2 mm * 2 mm, 4 mm * 4 mm, 6 mm * 6 mm, 8 mm *8 mm, 10 mm *10 mm, 12 mm * 12 mm. In each group, the finite element method and the above model were used to simulate the stress on the ankle joint when it was in the push-off phase, midstance phase and heel-strike phase to determine the contact pressure on the joint surface, the equivalent stress of the cartilage of the proximal talus and distal tibia in each phase, and the displacement of the talus. The stress, contact state and displacement of each component of the ankle joint in the different groups were observed to determine its maximum value and location. The maximum pressure was recorded as the experimental data and analyzed to obtain the column diagram, and the changes in pressure were discussed.