Validation of a personalized ligament-constraining discrete element framework for computing ankle joint contact mechanics

Optimal biomechanical choice of implant placement in various pilon fracture types: a finite element study

BACKGROUND

The research on the biomechanical characteristics of individual implant placement for pilon fractures based on the different initial direction of fracture displacement is still insufficient. This study's aim is to compare the stress distribution in bones and implants with various pilon fracture types.

METHODS

Varus, valgus, dorsiflexion, and plantarflexion type fractures were categorized as type I, II, III, and IV, respectively. The buttress plate was placed medially in subtypes IA and IIB, whereas it was placed anterolaterally in subtypes IB and IIA; The anterior or posterior buttress plate was utilized in subtypes IIIA and IVA, the lag screws were applied in subtypes IIIB and IVB. The maximum equivalent stress of tibia (TI-Smax) and implants (IF-Smax), stress of fracture fragments (Sfe), and axial displacement values of the fracture fragments (ADfe) in each subtype were analyzed when the ankle was in a neutral position, 15° of varus and valgus in types I and II, 15° of dorsiflexion and plantarflexion in types III and IV.

RESULTS

Under the same axial stress loading conditions, TI-Smax, Sfe, ADfe of subtypes IA and IIA were significantly lower than subtypes IB and IIB, while IF-Smax of subtypes IA and IIA were obviously larger than subtypes IB and IIB. Additionally, TI-Smax, Sfe, ADfe of subtypes IIIA and IVA were considerably lower as IF-Smax met expectations compared to subtypes IIIB and IVB.

CONCLUSION

Based on these results, when making decisions for open reduction and internal fixation in various pilon fractures, the choice and placement of the implant can be recommended as follows: the medial buttress plate for varus types; the anterolateral plate for valgus types; the anterior plate for dorsiflexion types; the posterior plate for plantarflexion types.

Congruent Weber-B ankle fractures do not alter tibiotalar contact mechanics

Current treatment strategy for managing Weber B ankle fractures is mainly governed by mortise congruency, malleolar alignment, deltoid ligament competence and fracture stability. While nonoperative treatment has yielded good functional outcomes in satisfactorily aligned stable injuries, a biomechanical rationale is not firmly established. Furthermore, current radiographic analysis is obscured by observer inaccuracy and beam rotation. This study aimed to utilize weightbearing CT and computational biomechanics to analyse 3D mortise displacement and contact mechanics in Weber-B ankle fractures and compare them with the uninjured contralateral side. 32 patients were included who sustained a Weber-B ankle fracture and underwent bilateral weightbearing CT imaging at injury. Segmentation into 3D models of bone was performed semi-automatically, and individualized cartilage layers were modeled based on a previously validated methodology. The 3D mortise congruency was evaluated by use of following parameters: alpha angle, fibular length, talocrural angle, medial gutter- and tibiofibular clear space- distance mapping. Contact mechanics were evaluated by the mean and maximum contact stress of the tibiotalar articulation. Statistical analysis revealed that there were no significant differences for all anatomical parameters. There was also no significant difference between the mean contact stress of the injured (2.10 ± 0.42 MPa) and uninjured side (2.10 ± 0.41 MPa), nor the maximum contact stress of the injured (7.67 ± 1.55 MPa) and uninjured (7.47 ± 1.67), respectively. Contact mechanics were unaffected in congruent Weber-B fractures. These findings support consideration of nonoperative management in such cases, given their presumably low risk for posttraumatic arthritic development in the long term. Level of Evidence: Level III; retrospective case-control study.

Hindfoot kinematics and kinetics - A combined in vivo and in silico analysis approach

BACKGROUND

The complex anatomical structure of the foot-ankle imposes challenges to accurately quantify detailed hindfoot kinematics and estimate musculoskeletal loading parameters. Most systems used to capture or estimate dynamic joint function oversimplify the anatomical structure by reducing its complexity.

RESEARCH QUESTION

Can four dimensional computed tomography (4D CT) imaging in combination with an innovative foot manipulator capture in vivo hindfoot kinematics during a simulated stance phase of walking and can talocrural and subtalar articular joint mechanics be estimated based on a detailed in silico musculoskeletal foot-ankle model.

METHODS

A foot manipulator imposed plantar/dorsiflexion and inversion/eversion representing a healthy stance phase of gait in 12 healthy participants while simultaneously acquiring 4D CT images. Participant-specific 3D hindfoot rotations and translations were calculated based on bone-specific anatomical coordinate systems. Articular cartilage contact area and contact pressure of the talocrural and subtalar joints were estimated using an extended foot-ankle model updated with an elastic foundation contact model upon prescribing the participant-specific rotations measured in the 4D CT measurement.

RESULTS

Plantar/dorsiflexion predominantly occurred at the talocrural joint (RoM 15.9±3.9°), while inversion/eversion (RoM 5.9±3.9°) occurred mostly at the subtalar joint, with the contact area being larger at the subtalar than at the talocrural joint. Contact pressure was evenly distributed between the talocrural and subtalar joint at the beginning of the simulated stance phase but was then redistributed from the talocrural to the subtalar joint with increasing dorsiflexion.

SIGNIFICANCE

In a clinical case study, the healthy participants were compared with four patients after surgically treaded intra-articular calcaneal fracture. The proposed workflow was able to detect small but meaningful differences in hindfoot kinematics and kinetics, indicative of remaining hindfoot pathomechanics that may influence the onset and progression of degenerative joint diseases.

A finite element model of human hindfoot and its application in supramalleolar osteotomy

BACKGROUND

The majority of the ankle osteoarthritis cases are posttraumatic and affect younger patients with a longer projected life span. Hence, joint-preserving surgery, such as supramalleolar osteotomy becomes popular among young patients, especially those with asymmetric arthritis due to alignment deformities. However, there is a lack of biomechanical studies on postoperative evaluation of stress at ankle joints. We aimed to construct a verifiable finite element model of the human hindfoot, and to explore the effect of different osteotomy parameters on the treatment of varus ankle arthritis.

METHODS

The bones of the hindfoot are reconstructed using normal CT tomography data from healthy volunteers, while the cartilages and ligaments are determined from the literature. The finite element calculation results are compared with the weight-bearing CT (WBCT) data to validate the model. By setting different model parameters, such as the osteotomy height (L) and the osteotomy distraction distance (h), the effects of different surgical parameters on the contact stress of the ankle joint surface are compared.

FINDINGS

The alignment and the deformation of hindfoot bones as determined by the finite element analysis aligns closely with the data obtained from WBCT. The maximum contact stress of the ankle joint surface calculated by this model increases with the increase of the varus angle. The maximum contact stresses as a function of the L and h of the ankle joint surface are determined.

INTERPRETATION

The relationship between surgical parameters and stress at the ankle joint in our study could further help guiding the planning of the supramalleolar osteotomy according to the varus/valgus alignment of the patients.

The Influence of Talar Displacement on Articular Contact Mechanics: A 3D Finite Element Analysis Study Using Weightbearing Computed Tomography

BACKGROUND

Talar displacement is considered the main predictive factor for poor outcomes and the development of post-traumatic osteoarthritis after ankle fractures. Isolated lateral talar translation, as previously studied by Ramsey and Hamilton using carbon powder imprinting, does not fully replicate the multidirectional joint subluxations seen in ankle fractures. The purpose of this study was to analyze the influence of multiple uniplanar talar displacements on tibiotalar contact mechanics utilizing weightbearing computed tomography (WBCT) and finite element analysis (FEA).

METHODS

Nineteen subjects (mean age = 37.6 years) with no history of ankle surgery or injury having undergone WBCT arthrogram (n = 1) and WBCT without arthrogram (n = 18) were included. Segmentation of the WBCT images into 3D simulated models of bone and cartilage was performed. Three-dimensional (3D) multiple uniplanar talar displacements were simulated to investigate the respective influence of various uniaxial displacements (including lateral translation, anteroposterior translation, varus-valgus angulation, and external rotation) on the tibiotalar contact mechanics using FEA. Tibiotalar peak contact stress and contact area were modeled for each displacement and its gradations.

RESULTS

Our modeling demonstrated that peak contact stress of the talus and tibia increased, whereas contact area decreased, with incremental displacement in all tested directions. Contact stress maps of the talus and tibia were computed for each displacement demonstrating unique patterns of pressure derangement. One millimeter of lateral translation resulted in 14% increase of peak talar contact pressure and a 3% decrease in contact area.

CONCLUSION

Our model predicted that with lateral talar translation, there is less noticeable change in tibiotalar contact area compared with prior studies whereas external rotation greater than 12 degrees had the largest effect on peak contact stress predictions.

LEVEL OF EVIDENCE

Level V, computational simulation study.

Quantifying walking speeds in relation to ankle biomechanics on a real-time interactive gait platform: a musculoskeletal modeling approach in healthy adults

Background: Given the inherent variability in walking speeds encountered in day-to-day activities, understanding the corresponding alterations in ankle biomechanics would provide valuable clinical insights. Therefore, the objective of this study was to examine the influence of different walking speeds on biomechanical parameters, utilizing gait analysis and musculoskeletal modelling. Methods: Twenty healthy volunteers without any lower limb medical history were included in this study. Treadmill-assisted gait-analysis with walking speeds of 0.8 m/s and 1.1 m/s was performed using the Gait Real-time Analysis Interactive Lab (GRAIL®). Collected kinematic data and ground reaction forces were processed via the AnyBody® modeling system to determine ankle kinetics and muscle forces of the lower leg. Data were statistically analyzed using statistical parametric mapping to reveal both spatiotemporal and magnitude significant differences. Results: Significant differences were found for both magnitude and spatiotemporal curves between 0.8 m/s and 1.1 m/s for the ankle flexion (p < 0.001), subtalar force (p < 0.001), ankle joint reaction force and muscles forces of the M. gastrocnemius, M. soleus and M. peroneus longus (α = 0.05). No significant spatiotemporal differences were found between 0.8 m/s and 1.1 m/s for the M. tibialis anterior and posterior. Discussion: A significant impact on ankle joint kinematics and kinetics was observed when comparing walking speeds of 0.8 m/s and 1.1 m/s. The findings of this study underscore the influence of walking speed on the biomechanics of the ankle. Such insights may provide a biomechanical rationale for several therapeutic and preventative strategies for ankle conditions.

Supramalleolar Osteotomies in Cavovarus Foot Deformity: Why Patient-Specific Instruments Make a Difference

Supramalleolar osteotomy enables correction of the ankle varus deformity and is associated with improvement of pain and function in the short term and long term. Despite these beneficial results, the amount of surgical correction is challenging to titrate and the procedure remains technically demanding. Most supramalleolar osteotomies are currently planned preoperatively on 2-dimensional weight-bearing radiographs and executed peroperatively using free-hand techniques. This article encompasses 3-dimensional planning and printing techniques based on weight-bearing computed tomography images and patient-specific instruments to correct ankle varus deformities.

Application of external torque enhances the detection of subtle syndesmotic ankle instability in a weight-bearing CT

PURPOSE

Acute syndesmotic ankle injuries continue to impose a diagnostic dilemma and it remains unclear whether weightbearing and/or external rotation should be added during the imaging process. Therefore, the aim of this study was to assess if combined weightbearing and external rotation increases the diagnostic sensitivity of syndesmotic ankle instability using weightbearing CT (WBCT) imaging, compared to isolated weightbearing.

METHODS

In this retrospective study, patients with an acute syndesmotic ankle injury were analysed using a WBCT (N = 21; Age = 31.6 ± 14.1 years old). Inclusion criteria were an MRI confirmed syndesmotic ligament injury imaged by a WBCT of the ankle during weightbearing and combined weightbearing-external rotation. Exclusion criteria consisted of fracture associated syndesmotic injuries. Three-dimensional (3D) models were generated from the CT slices. Tibiofibular displacement and talar rotation were quantified using automated 3D measurements (anterior tibiofibular distance (ATFD), Alpha angle, posterior Tibiofibular distance (PTFD) and Talar rotation (TR) angle in comparison to the contralateral non-injured ankle.

RESULTS

The difference in neutral-stressed Alpha angle and ATFD showed a significant difference between patients with a syndesmotic ankle lesion and contralateral control (P = 0.046 and P = 0.039, respectively). The difference in neutral-stressed PTFD and TR angle did not show a significant difference between patients with a syndesmotic ankle lesion and healthy ankles (n.s.).

CONCLUSION

Application of combined weightbearing-external rotation reveals an increased ATFD in patients with syndesmotic ligament injuries. This study provides the first insights based on 3D measurements to support the potential relevance of applying external rotation during WBCT imaging. In clinical practice, this could enhance the current diagnostic accuracy of subtle syndesmotic instability in a non-invasive manner. However, to what extent certain displacement patterns require operative treatment strategies has yet to be determined in future studies.

LEVEL OF EVIDENCE

Level III.

Association between the distal tibiofibular syndesmosis morphology classification and ankle osteoarthritis: a retrospective study

BACKGROUND

Syndesmosis injury is proposed to contribute to ankle stability and osteoarthritis (OA). However, whether distal tibiofibular syndesmosis structure is closely related to ankle OA is unclear. We hypothesized that different DTS morphology classifications would affect the biomechanics properties in ankle OA. The study aimed to determine the association between the distal tibiofibular syndesmosis (DTS) morphological classification and ankle OA.

METHODS

This is a retrospective study examining imaging data of 147 patients (87 males and 60 females) with ankle OA. Magnetic resonance imaging was used to access the DTS morphological classification, according to measuring various parameters. Joint space narrowing and osteophytes were measured using ankle weight-bearing radiography. The classification and parameters were analyzed to determine the relationship between the syndesmosis classification and the abnormality of ankle OA.

RESULTS

Five morphological classifications of the DTS, including Chevron (19.6%), Widow's peak (16.2%), Flat (22.3%), Trapezoid (32.0%), and Crescent (19.6%), were shown. There were statistical differences between DTS classification and tibial angle surface angle (TAS) (P = .009) and talar tilt angle (TTA) (P = .014). The TAS (degree) of the Crescent (86.47 ± 3.21) was less than Chevron (88.75 ± 2.72) (P = .006), Widow's peak (89.26 ± 3.15) (P = .001), Flat (88.83 ± 3.62) (P = .003) and Trapezoid (88.11 ± 2.62) (P = .041), respectively. The TTA (degree) of Crescent (86.83 ± 5.30) was less than Chevron (89.28 ± 2.46) and Widow's peak (89.82 ± 3.41). The men were greater than women for TAS (P = .008) and angle (P = .003), which are consistent with osteophyte (P = .019) and the modified Kellgren-Lawrence grades (P = .041) between gender.

CONCLUSIONS

DTS morphological classification might affect the biomechanics properties in TAS and TTA in ankle OA. In clinical practice, surgeons should pay attention to the effects of DTS on ankle OA.

LEVEL OF EVIDENCE

Level III, retrospective study.