An instrument for methodological quality assessment of single-subject finite element analysis used in computational orthopaedics

Effect of forefoot transverse arch stiffness on foot biomechanical response--based on finite element method

Background

The plantar vault, comprising the transverse and longitudinal arches of the human foot, is essential for impact absorption, elastic energy storage, and propulsion. Recent research underscores the importance of the transverse arch, contributing over 40% to midfoot stiffness. This study aimed to quantify biomechanical responses in the ankle-foot complex by varying the stiffness of the deep metatarsal transverse ligament (DTML).

Methods

Using CT image reconstruction, we constructed a complex three-dimensional finite element model of the foot and ankle joint complex, accounting for geometric complexity and nonlinear characteristics. The focus of our study was to evaluate the effect of different forefoot transverse arch stiffness, that is, different Young's modulus values of DTML (from 135 MPa to 405 MPa), on different biomechanical aspects of the foot and ankle complex. Notably, we analyzed their effects on plantar pressure distribution, metatarsal stress patterns, navicular subsidence, and plantar fascial strain.

Results

Increasing the stiffness of the DTML has significant effects on foot biomechanics. Specifically, higher DTML stiffness leads to elevate von Mises stress in the 1st, 2nd, and 3rd metatarsals, while concurrently reducing plantar pressure by 14.2% when the Young's modulus is doubled. This stiffening also impedes navicular bone subsidence and foot lengthening. Notably, a 100% increase in the Young's modulus of DTML results in a 54.1% decrease in scaphoid subsidence and a 2.5% decrease in foot lengthening, which collectively contribute to a 33.1% enhancement in foot longitudinal stiffness. Additionally, doubling the Young's modulus of DTML can reduce the strain stretch of the plantar fascia by 38.5%.

Conclusion

Preserving DTML integrity sustains the transverse arch, enhancing foot longitudinal stiffness and elastic responsiveness. These findings have implications for treating arch dysfunction and provide insights for shoe developers seeking to enhance propulsion.

Biomechanical characteristics of Sanders type II and III calcaneal fractures fixed by open reduction and internal fixation and percutaneous minimally invasive fixation

BACKGROUND

This work investigated the differences in the biomechanical properties of open reduction and internal fixation (ORIF) and percutaneous minimally invasive fixation (PMIF) for the fixation of calcaneal fractures (Sanders type II and III calcaneal fractures as examples) through finite element analysis.

METHODS

Based on CT images of the human foot and ankle, according to the principle of three-point fixation, namely the sustentaculum tali, the anterior process and the calcaneal tuberosity were fixed. Three-dimensional finite element models of Sanders type II and III calcaneal fractures fixed by ORIF and PMIF were established. The proximal surfaces of the tibia, fibula and soft tissue were constrained, and ground reaction force and Achilles tendon force loads were added to simulate balanced standing.

RESULTS

The maximum stress was 80.54, 211.59 and 113.88 MPa for the calcaneus, screws and plates in the ORIF group and 70.02 and 209.46 MPa for the calcaneus and screws in the PMIF group, respectively; the maximum displacement was 0.26, 0.21 and 0.12 mm for the calcaneus, screws and plates in the ORIF group and 0.20 and 0.14 mm for the calcaneus and screws in the PMIF group, respectively. The values obtained from the simulation were within the permissible stress and elastic deformation range of the materials used in the model, and there was no significant stress concentration. The maximum stress and displacement of the calcaneus and implants were slightly lower in the PMIF group than in the ORIF group when fixing Sanders type II and III calcaneal fractures.

CONCLUSIONS

This study may provide a reference for optimising the design of implants, the development of individualised preoperative plans and the choice of clinical surgical approach.

Effects of plantar fascia stiffness on the internal mechanics of idiopathic pes cavus by finite element analysis: implications for metatarsalgia

Metatarsalgia occurring in individuals with pes cavus is typically associated with abnormal loading patterns in the forefoot resulting from structural alterations. Simultaneously, the frequent overstress of the plantar fascia (PF) caused by the persistence of this foot deformity may further exacerbate the chronic pain induced by metatarsal overload. We aimed to investigate and quantify the effects of PF stiffness on the internal biomechanics of pes cavus using a computational modelling approach. A patient-specific finite element model of the foot-ankle complex using the actual three-dimensional geometry of idiopathic pes cavus bones and soft tissues was reconstructed. A sensitivity study was conducted to evaluate the effects of varying elastic modulus (0-700 MPa) of the PF on the metatarsal stress distribution, and force transmission through the metatarsophalangeal (MTP) and tarsometatarsal (TMT) joints in the pes cavus. The results indicated that variations in PF stiffness led to stress redistribution in the metatarsal region. Peak stress gradually reduced with decreasing stiffness until the PF was released, eventually resulting in a reduction of 22.39% compared to the reference value of 350 MPa. Furthermore, adjusting the PF stiffness to twice the reference value (700 MPa) increased the contact forces through the TMT and MTP joints by up to 23% and 116%, respectively. The reduction of PF stiffness alleviated focal metatarsal loading, and therefore, surgical fascia release can be considered to alleviate metatarsalgia in patients with pes cavus.

Computational Biomechanics of Sleep: A Systematic Mapping Review

Biomechanical studies play an important role in understanding the pathophysiology of sleep disorders and providing insights to maintain sleep health. Computational methods facilitate a versatile platform to analyze various biomechanical factors in silico, which would otherwise be difficult through in vivo experiments. The objective of this review is to examine and map the applications of computational biomechanics to sleep-related research topics, including sleep medicine and sleep ergonomics. A systematic search was conducted on PubMed, Scopus, and Web of Science. Research gaps were identified through data synthesis on variants, outcomes, and highlighted features, as well as evidence maps on basic modeling considerations and modeling components of the eligible studies. Twenty-seven studies (n = 27) were categorized into sleep ergonomics (n = 2 on pillow; n = 3 on mattress), sleep-related breathing disorders (n = 19 on obstructive sleep apnea), and sleep-related movement disorders (n = 3 on sleep bruxism). The effects of pillow height and mattress stiffness on spinal curvature were explored. Stress on the temporomandibular joint, and therefore its disorder, was the primary focus of investigations on sleep bruxism. Using finite element morphometry and fluid-structure interaction, studies on obstructive sleep apnea investigated the effects of anatomical variations, muscle activation of the tongue and soft palate, and gravitational direction on the collapse and blockade of the upper airway, in addition to the airflow pressure distribution. Model validation has been one of the greatest hurdles, while single-subject design and surrogate techniques have led to concerns about external validity. Future research might endeavor to reconstruct patient-specific models with patient-specific loading profiles in a larger cohort. Studies on sleep ergonomics research may pave the way for determining ideal spine curvature, in addition to simulating side-lying sleep postures. Sleep bruxism studies may analyze the accumulated dental damage and wear. Research on OSA treatments using computational approaches warrants further investigation.

The effect of thickness and elastic modulus of the anterior talofibular ligament on anterior ankle joint stiffness: A subject-specific finite element study

Ankle sprain is a frequent type of sports injury leading to lateral ligament injury. The anterior talofibular ligament (ATFL) is a primary ligamentous stabilizer of the ankle joint and typically the most vulnerable ligament injured in a lateral ankle sprain (LAS). This study aimed to quantitively investigate the effect of the thickness and elastic modulus of ATFL on anterior ankle joint stiffness (AAJS) by developing nine subject-specific finite element (FE) models under acute injury, chronic injury, and control conditions of ATFL. A 120 N forward force was applied at the posterior calcaneus leading to an anterior translation of the calcaneus and talus to simulate the anterior drawer test (ADT). In the results, the ratio of the forward force to the talar displacement was used to assess the AAJS, which increased by 5.85% in the acute group and decreased by 19.78% in the chronic group, compared to those of the control group. An empirical equation described the relationship between AAJS, thickness, and elastic modulus (R-square 0.98). The equation proposed in this study provided an approach to quantify AAJS and revealed the effect of the thickness and the elastic modulus of ATFL on ankle stability, which may shed light on the potential diagnosis of lateral ligament injury.

Attachments for the Orthodontic Aligner Treatment-State of the Art-A Comprehensive Systematic Review

BACKGROUND

In recent years the burden of aligner treatment has been growing. However, the sole use of aligners is characterized by limitations; thus attachments are bonded to the teeth to improve aligner retention and tooth movement. Nevertheless, it is often still a challenge to clinically achieve the planned movement. Thus, the aim of this study is to discuss the evidence of the shape, placement and bonding of composite attachments.

METHODS

A query was carried out in six databases on 10 December 2022 using the search string ("orthodontics" OR "malocclusion" OR "Tooth movement techniques AND ("aligner*" OR "thermoformed splints" OR "invisible splint*" AND ("attachment*" OR "accessor*" OR "auxill*" AND "position*").

RESULTS

There were 209 potential articles identified. Finally, twenty-six articles were included. Four referred to attachment bonding, and twenty-two comprised the influence of composite attachment on movement efficacy. Quality assessment tools were used according to the study type.

CONCLUSIONS

The use of attachments significantly improves the expression of orthodontic movement and aligner retention. It is possible to indicate sites on the teeth where attachments have a better effect on tooth movement and to assess which attachments facilitate movement. The research received no external funding. The PROSPERO database number is CRD42022383276.

Biomechanical evaluation of two modified intramedullary fixation system for treating unstable femoral neck fractures: A finite element analysis

Purpose: The existing implants for fixation of femoral neck fractures have poor biomechanical stability, so the failure rate is high. We designed two modified intramedullary implants for treating unstable femoral neck fractures (UFNFs). We tried to improve the biomechanical stability of fixation by shortening the moment and reducing stress concentration. Each modified intramedullary implant was compared with cannulated screws (CSs) through finite element analysis (FEA). Methods: Five different models were included: three cannulated screws (CSs, Model 1) in an inverted triangle configuration, the dynamic hip screw with an anti-rotation screw (DHS + AS, Model 2), the femoral neck system (FNS, Model 3), the modified intramedullary femoral neck system (IFNS, Model 4), and the modified intramedullary interlocking system (IIS, Model 5). Three-dimensional (3D) models of femur and implants were constructed by using 3D modelling software. Three load cases were simulated to assess the maximal displacement of models and fracture surface. The maximal stress at the bone and implants was also evaluated. Results: FEA data showed that Model 5 had the best performance in terms of maximum displacement while Model 1 had the worst performance for this index under axial load of 2100 N. With respect to Maximum stress, Model 4 had the best performance while Model 2 had the worst performance under axial load. The general trends under bending and torsion load were consistent with that under axial load. Our data demonstrated that the two modified intramedullary implants exhibited the best biomechanical stability, followed by FNS and DHS + AS, and then three cannulated screws in axial, bending, and torsion load cases. Conclusion: The two modified intramedullary designs showed the best biomechanical performance among the five implants included in this study. Therefore, this might provide some new options for trauma surgeons to deal with unstable femoral neck fractures.

Design feature combinations effects of running shoe on plantar pressure during heel landing: A finite element analysis with Taguchi optimization approach

Large and repeated impacts on the heel during running are among the primary reasons behind runners’ injuries. Reducing plantar pressure can be conducive to reducing running injury and improving running performance and is primarily achieved by modifying the design parameters of running shoes. This study examines the effect of design parameters of running shoes (i.e., heel-cup, insole material, midsole material, and insole thickness) on landing peak plantar pressure and determines the combination of different parameters that optimize cushion effects by employing the Taguchi method. We developed the foot–shoe finite element (FE) model through reverse engineering. Model assembly with different design parameters was generated in accordance with the Taguchi method orthogonal table. The effectiveness of the model was verified using the static standing model in Ansys. The significance and contribution of different design parameters, and the optimal design to reduce plantar pressure during landing, were determined using the Taguchi method. In the descending order of percentage contribution was a conforming heel-cup (53.18%), insole material (25.89%), midsole material (7.81%), and insole thickness (2.69%). The more conforming heel-cup (p < 0.001) and softer insole (p = 0.001) reduced the heel pressure during landing impact. The optimal design of running shoe in this study was achieved with a latex insole, a 6 mm insole thickness, an Asker C-45 hardness midsole, and a 100% conforming heel-cup. The conforming heel-cup and the insole material significantly affected the peak plantar pressure during heel landing. The implementation of a custom conforming heel-cup is imperative for relieving high plantar pressure for long-distance heel-strike runners.

How Does Orthodontic Mini-Implant Thread Minidesign Influence the Stability?—Systematic Review with Meta-Analysis

Background: Clinical guidelines are lacking for the use of orthodontic mini-implants (OMIs) in terms of scientific evidence referring to the choice of proper mini-design. Thus, the present study aimed to investigate to what extent orthodontic mini-implant thread design influences its stability. Methods: Search was conducted in five search engines on 10 May. Quality assessment was performed using study type specific scales. Whenever possible, meta-analysis was performed. Results: The search strategy identified 118 potential articles. Twenty papers were subjected to qualitative analysis and data from 8 papers—to meta-analysis. Studies included were characterized by high or medium quality. Four studies were considered as low quality. No clinical studies considering the number of threads, threads depth, or TSF have been found in the literature. Conclusions: Minidesign of OMIs seems to influence their stability in the bone. Thread pitch seems to be of special importance for OMIs retention—the more dense thread—the better stability. Thread depth seems to be of low importance for OMIs stability. There is no clear scientific evidence for optimal thread shape factor. Studies present in the literature vary greatly in study design and results reporting. Research received no external funding. Study protocol number in PROSPERO database: CRD42022340970.

Tibio-Femoral Contact Force Distribution of Knee Before and After Total Knee Arthroplasty: Combined Finite Element and Gait Analysis

Objective

To assess the tibio‐femoral contact forces before and after total knee arthroplasty (TKA) in patients with knee osteoarthritis (KOA) by three‐dimensional (3D) finite element analysis (FEA) models and gait analysis.

Methods

Two hospitalized patients with Kellgren–Lawrence grade IV varus KOA and two healthy subjects were enrolled in this study. Both patients underwent unilateral TKA. FEA models were established based on CT and MR images of the knees of the patients with KOA and healthy subjects. Gait analysis was performed using a three‐dimensional motion capture system with a force plate. Three direction forces at the ankle joints were calculated by inverse dynamic analysis, which provided the load for the FEA models. The total contact forces of the knee joints were also calculated by inverse dynamic analysis to enable comparisons with the results from the FEA models. The total knee contact forces, maximum von Mises stress, and stress distribution of the medial plateau were compared between the patients and healthy subjects. The distributions of the medial plateau force at 2 and 6 months postoperatively were compared with the distributions of the forces preoperatively and those in the healthy subjects.

Results

During static standing, the medial plateau bore the most of the total contact forces in the knees with varus KOA (90.78% for patient 1 and 93.53% for patient 2) compared with 64.75 ± 3.34% of the total force in the healthy knees. At the first and second peaks of the ground reaction force during the stance phase of a gait cycle, the medial plateau bore a much higher percentage of contact forces in patients with KOA (74.78% and 86.48%, respectively, for patient 1; 70.68% and 83.56%, respectively, for patient 2) than healthy subjects (61.06% ± 3.43% at the first peak and 72.09% ± 1.83% at the second peak). Two months after TKA, the percentages of contact forces on the medial tibial plateau were 79.65%–85.19% at the first and second peaks of ground reaction forces during the stance phase of a gait cycle, and the percentages decreased to 53.99% – 68.13% 6 months after TKA.

Conclusion

FEA showed that TKA effectively restored the distribution of tibio‐femoral contact forces during static standing and walking, especially 6 months after the surgery. The changes in the gait were consistent with the changes in the contact force distribution calculated by the FEA model.

Finite element analysis of shank and ankle with different boot collar heights in parachuting landing on inversion ground surface

This study aimed to the biomechanics of the foot-ankle-shank complex with different boots collar heights in parachuting landing on inversion ground surface. A finite element model including tibia, fibula, ankle, foot and parachuting boot was developed. Three collar heights (low, medium, high) of the parachuting boot were simulated. Von-Mises stress, ankle inversion angle, ligament force and bone displacement were analyzed. Compared with that of the high and low collar heights, boots with medium collar height produced the lowest peak stress on the tibia and the articular cartilage of the subtalar joint. In addition, the medium collar height can better control the ankle inversion and minimize the tensile forces on the lateral ankle ligaments.

Different Design Feature Combinations of Flatfoot Orthosis on Plantar Fascia Strain and Plantar Pressure: A Muscle-Driven Finite Element Analysis With Taguchi Method

Customized foot orthosis is commonly used to modify foot posture and relieve foot pain for adult acquired flexible flatfoot. However, systematic investigation of the influence of foot orthotic design parameter combination on the internal foot mechanics remains scarce. This study aimed to investigate the biomechanical effects of different combinations of foot orthoses design features through a muscle-driven flatfoot finite element model. A flatfoot-orthosis finite element model was constructed by considering the three-dimensional geometry of plantar fascia. The plantar fascia model accounted for the interaction with the bulk soft tissue. The Taguchi approach was adopted to analyze the significance of four design factors combination (arch support height, medial posting inclination, heel cup height, and material stiffness). Predicted plantar pressure and plantar fascia strains in different design combinations at the midstance instant were reported. The results indicated that the foot orthosis with higher arch support (45.7%) and medial inclination angle (25.5%) effectively reduced peak plantar pressure. For the proximal plantar fascia strain, arch support (41.8%) and material stiffness (37%) were strong influencing factors. Specifically, higher arch support and softer material decreased the peak plantar fascia strain. The plantar pressure and plantar fascia loading were sensitive to the arch support feature. The proposed statistics-based finite element flatfoot model could assist the insole optimization and evaluation for individuals with flatfoot.

Biomechanical Analysis of a Novel Double-Point Fixation Method for Displaced Intra-Articular Calcaneal Fractures

The development of minimally invasive procedures and implant materials has improved the fixation strength of implants and is less traumatic in surgery. The purpose of this study was to propose a novel “double-point fixation” for calcaneal fractures and compare its biomechanical stability with the traditional “three-point fixation.” A three-dimensional finite element foot model with a Sanders type IIIAB calcaneal fracture was developed based on clinical images comprising bones, plantar fascia, ligaments, and encapsulated soft tissue. Double-point and three-point fixation resembled the surgical procedure with a volar distal radius plate and calcaneal locking plate, respectively. The stress distribution, fracture displacement, and change of the Böhler angle and Gissane’s angle were estimated by a walking simulation using the model, and the predictions between the double-point and three-point fixation were compared at heel-strike, midstance, and push-off instants. Double-point fixation demonstrated lower bone stress (103.3 vs. 199.4 MPa), but higher implant stress (1,084.0 vs. 577.9 MPa). The model displacement of double-point fixation was higher than that of three-point fixation (3.68 vs. 2.53 mm). The displacement of the posterior joint facet (0.127 vs. 0.150 mm) and the changes of the Böhler angle (0.9° vs. 1.4°) and Gissane’s angle (0.7° vs. 0.9°) in double-point fixation were comparably lower. Double-point fixation by volar distal radius plates demonstrated sufficient and favorable fixation stability and a lower risk of postoperative stress fracture, which may potentially serve as a new fixation modality for the treatment of displaced intra-articular calcaneal fractures.

Risk of proximal femoral nail antirotation (PFNA) implant failure upon different lateral femoral wall thickness in intertrochanteric fracture: a finite element analysis

Proximal Femoral Nail Antirotation (PFNA) has been commonly used to treat intertrochanteric fractures, despite the risk of implant failure. The integrity of the femur could influence the risk of implant failure. This study evaluated the influence of lateral femoral wall thickness on the potential of implant failure. A finite element model of the hip was reconstructed from the Computed Tomography of a female patient. Five intertrochanteric fracture models at different lateral femoral wall thickness (T1 = 27.6 mm, T2 = 25.4 mm, T3 = 23.4 mm, T4 = 21.4 mm, and T5 = 19.3 mm) were then created and fixed with PFNA. We simulated a critical loading condition by a high loading case during walking. Elastoplastic material models with yield stress and failure strain were applied to the bone and implant in which breakage can be simulated using the element deletion function. In addition, the stress and displacement of the implant and femur were analysed. Implant breakage occurred at the sides of the proximal nail canal in cases of T4 and T5 which was further supported by the higher maximum von Mises stress and nail displacement. The increased stress and displacement of the implant may implicate a reduction of stability and risk of implant failure. We suggested that precaution shall be taken when the wall thickness was less than 21.4 mm.

Finite element analysis of lumbar spine with different backpack positions in parachuting landing

The purpose of this study was to investigate the lumbar spine stress with different backpack positions in parachuting landing using a finite element model of lumbar vertebra 1-5. The backpack gravity center was set at three positions (posterior-high (case PH), posterior-low (case PL), and anterior-low (case AL)) respectively. In results, the peak Von-Mises stresses of the matrix, nucleus, fibers, endplate and ligament in case AL were 2.765 MPa, 0.534 MPa, 6.561 MPa, 4.045 MPa and 1.790 MPa respectively, lower than those in case PL (6.913 MPa, 1.316 MPa, 20.716 MPa, 10.917 MPa and 5.147 MPa respectively) and case PH (7.328 MPa, 1.394 MPa, 22.147 MPa, 11.617 MPa and 5.464 MPa respectively). In conclusion, setting the gravity center of backpack at anterior-low position would reduce lumbar spine stress and reduce lumbar spine injuries.