Analysis of stress and stabilization in adolescent with osteoporotic idiopathic scoliosis: finite element method

Imaging observation of intervertebral disc degeneration in patients with old thoracolumbar fracture-related kyphotic deformity

Old thoracolumbar fracture with kyphosis (OTLFK) often results in low back pain, with intervertebral disc degeneration being a significant contributor. We hypothesized that patients with OTLFK exhibit distinct patterns of disc degeneration compared to those with chronic low back pain without kyphotic deformity. This study aimed to investigate the characteristics of disc degeneration in OTLFK patients and explore its association with sagittal spinal parameters and endplate injury. A retrospective analysis was conducted on 52 patients with OTLFK (observation group, OG) and 104 age-, gender-, and BMI-matched patients with chronic low back pain (control group, CG) treated at our hospital between February 2017 and March 2023. Intervertebral disc degeneration from T11/12 to L5/S1 was assessed using MRI T2-weighted images and the Pfirrmann grading system. Sagittal spinal parameters-including lumbar lordosis (LL), thoracic kyphosis (TK), local kyphosis Cobb angle (LKCA), thoracolumbar kyphosis (TLK), pelvic tilt(PT), sacral slope(SS), and sagittal vertical axis (SVA)-and endplate injury grades were measured in the OG. Differences in disc degeneration between the two groups were compared, and correlations between disc degeneration, sagittal parameters, and endplate injury were analyzed. The OG exhibited significantly higher overall disc degeneration grades compared to the CG (p < 0.05), particularly at levels T11/12, T12/L1, L1/2, and L2/3. In the OG, grade IV and V degenerations were predominantly found from T11/12 to L2/3, whereas in the CG, they were mainly at L4/5 and L5/S1. Disc degeneration in the OG was significantly correlated with sagittal parameters and endplate injury grades (p < 0.05). Patients with OTLFK have higher grades of disc degeneration in the thoracolumbar region compared to those with chronic low back pain without kyphosis. Disc degeneration in OTLFK is associated with abnormal sagittal alignment and endplate injury, suggesting that kyphotic deformity and altered spinal biomechanics contribute to accelerated disc degeneration.

Characteristics and mechanical mechanisms of intervertebral disc degeneration in old thoracolumbar fractures with kyphosis: clinical observations and finite element analyses

BACKGROUND

Low back pain is a common complication in patients with old thoracolumbar fractures with kyphosis (OTLFK), and intervertebral disc degeneration (IDD) is a major contributor. Mechanical abnormalities are believed to play a key role in the development to IDD. This study aimed to investigate the characteristics of lumbar disc degeneration and underlying mechanical mechanisms in patients with OTLFK.

METHODS

A total of 52 patients with OTLFK were included from February 2017 to March 2023 as the observation group (OG). A control group (CG) of individuals with chronic low back pain were matched for age, body mass index, and gender. The disc degeneration grades and distribution in both groups were observed. Intact, 20°, 30°, and 40° kyphotic finite element (FE) models were established. Intervertebral disc pressures (IDPs) were calculated under standing, flexion, extension, lateral bending, and axial rotation conditions.

RESULTS

The overall IDD in the OG was higher than that in the CG. The grades from T11 to L3 were higher in the OG (p < 0.05), while there was no significant difference from L4 to S1. Degeneration levels IV and V were concentrated in the T11-L3 segment in the OG; whereas, in the CG, this was in L4/5 and L5/S1 (p < 0.05). The FE analysis results showed that, in the kyphotic model, the IDP was higher than the intact model in the standing position, flexion, lateral bending, and rotation, but lower in extension.

CONCLUSIONS

Patients with OTLFK exhibit higher-grade disc degeneration concentrated in the thoracolumbar segment. Abnormal mechanical stress may contribute to this degeneration, highlighting the importance of managing stress in kyphotic deformities.

Contribution of ankle motion pattern during landing to reduce the knee-related injury risk

BACKGROUND

Single-leg landing (SL) is an essential technique in sports such as basketball, soccer, and volleyball, which is often associated with a high risk of knee-related injury. The ankle motion pattern plays a crucial role in absorbing the load shocks during SL, but the effect on the knee joint is not yet clear. This work aims to explore the effects of different ankle plantarflexion angles during SL on the risk of knee-related injury.

METHODS

Thirty healthy male subjects were recruited to perform SL biomechanics tests, and one standard subject was selected to develop the finite element model of foot-ankle-knee integration. The joint impact force was used to evaluate the impact loads on the knee at various landing angles. The internal load forces (musculoskeletal modeling) and stress (finite element analysis) around the knee joint were simulated and calculated to evaluate the risk of knee-related injury during SL. To more realistically revert and simulate the anterior cruciate ligament (ACL) injury mechanics, we developed a knee musculoskeletal model that reverts the ACL ligament to a nonlinear short-term viscoelastic mechanical mechanism (strain rate-dependent) generated by the dense connective tissue as a function of strain.

RESULTS

As the ankle plantarflexion angle increased during landing, both the peak knee vertical impact force (p = 0.001) and ACL force (p = 0.001) decreased significantly. The maximum von Mises stress of ACL, meniscus, and femoral cartilage decreased as the ankle plantarflexion angle increased. The overall range of variation in ACL stress was small and was mainly distributed in the femoral and tibial attachment regions, as well as in the mid-lateral region.

CONCLUSION

The current findings revealed that the use of larger ankle plantarflexion angles during landing may be an effective solution to reduce knee impact load and the risk of rupture of the medial femoral attachment area in the ACL. The findings of this study have the potential to offer novel perspectives in the optimized application of landing strategies, thus giving crucial theoretical backing for decreasing the risk of knee-related injury.

Influence of the biomechanical evaluation of rupture using two shapes of same intramedullary implant after proximal interphalangeal joint arthrodesis to correct the claw/hammer pathology: A finite element study

We used finite element analysis to study the mechanical stress distribution of a new intramedullary implant used for proximal interphalangeal joint (PIPJ) arthrodesis (PIPJA) to surgically correct the claw-hammer toe deformity that affects 20% of the population. After geometric reconstruction of the foot skeleton from claw toe images of a 36-year-old male patient, two implants were positioned, in the virtual model, one neutral implant (NI) and another one 10° angled (10°AI) within the PIPJ of the second through fourth HT during the toe-off phase of gait and results were compared to those derived for the non-surgical foot (NSF). A PIPJA was performed on the second toe using a NI reduced tensile stress at the proximal phalanx (PP) (45.83 MPa) compared to the NSF (59.44 MPa; p < 0.001). When using the 10°AI, the tensile stress was much higher at PP and middle phalanges (MP) of the same toe, measuring 147.58 and 160.58 MPa, respectively, versus 59.44 and 74.95 MPa at corresponding joints in the NSF (all p < 0.001). Similar results were found for compressive stresses. The NI reduced compressive stress at the second PP (-65.12 MPa) compared to the NSF (-113.23 MPa) and the 10°AI (-142 MPa) (all p < 0.001). The von Mises stresses within the implant were also significantly lower when using NI versus 10°AI (p < 0.001). Therefore, we do not recommend performing a PIPJA using the 10°AI due to the increase in stress concentration primarily at the second PP and MP, which could promote implant breakage.

Finite element analysis of precise puncture vertebral augmentation in the treatment of different types of osteoporotic vertebral compression fractures

BACKGROUND

Osteoporosis vertebral compression fracture (OVCF) secondary to osteoporosis is a common health problem in the elderly population. Vertebral augmentation (VA) has been widely used as a minimally invasive surgical method. The transpedicle approach is commonly used for VA puncture, but sometimes, it is limited by the anatomy of the vertebral body and can not achieve good surgical results. Therefore, we propose the treatment of OVCF with precise puncture vertebral augmentation (PPVA). This study used finite element analysis to explore the biomechanical properties of PPVA in the treatment of osteoporotic vertebral compression fractures (OVCFs) with wedge, biconcave, and collapse deformities.

METHOD

Three-dimensional finite element models of the fractured vertebral body and adjacent superior and inferior vertebral bodies were established using Computed Tomography (CT) data from patients with OVCF, both before and after surgery. Evaluate the stress changes of the wedged deformed vertebral body, biconcave deformed vertebral body, collapsed deformed vertebral body, and adjacent vertebral bodies before and after PPVA.

RESULT

In vertebral bodies with wedge deformity and collapsed deformity, PPVA can effectively reduce the stress on the vertebral body but increases the stress on the vertebral body with biconcave deformity. PPVA significantly decreases the stress on the adjacent vertebral bodies of the wedge deformed vertebral body, and decreases the stress on the adjacent superior vertebral body of biconcave deformity and collapsed deformed vertebral bodies, but increases the stress on the adjacent inferior vertebral bodies. PPVA improves the stress distribution of the vertebral body and prevents high-stress areas from being concentrated on one side of the vertebral body.

CONCLUSION

PPVA has shown positive surgical outcomes in treating wedge deformed and collapsed deformed vertebral bodies. However, its effectiveness in treating biconcave vertebral body is limited. Furthermore, PPVA has demonstrated favorable results in addressing adjacent superior vertebral body in three types of fractures.

Can the Entire Function of the Foot Be Concentrated in the Forefoot Area during the Running Stance Phase? A Finite Element Study of Different Shoe Soles

The goal of this study was to use the finite element (FE) method to compare and study the differences between bionic shoes (BS) and normal shoes (NS) forefoot strike patterns when running. In addition, we separated the forefoot area when forefoot running as a way to create a small and independent area of instability. An adult male of Chinese descent was recruited for this investigation (age: 26 years old; body height: 185 cm; body mass: 82 kg) (forefoot strike patterns). We analyzed forefoot running under two different conditions through FE analysis, and used bone stress distribution feature classification and recognition for further analysis. The metatarsal stress values in forefoot strike patterns with BS were less than with NS. Additionally, the bone stress classification of features and the recognition accuracy rate of metatarsal (MT) 2, MT3 and MT5 were higher than other foot bones in the first 5%, 10%, 20% and 50% of nodes. BS forefoot running helped reduce the probability of occurrence of metatarsal stress fractures. In addition, the findings further revealed that BS may have important implications for the prevention of hallux valgus, which may be more effective in adolescent children. Finally, this study presents a post-processing method for FE results, which is of great significance for further understanding and exploration of FE results.

Asymmetric Load Transmission Induces Facet Joint Subchondral Sclerosis and Hypertrophy in Patients with Idiopathic Adolescent Scoliosis: Evaluation Using Finite Element Model and Surgical Specimen

Adolescent idiopathic scoliosis (AIS) with thoracic curvature primarily progresses from the thoracolumbar region, causing abnormal twisting and rotation of the spinal column. This results in unbalanced, asymmetric loads on each vertebrae and increased demands on the thoracic facet joints to withstand rotational stress from adjacent vertebrae. However, no studies have focused on the stress distribution on the facet joints of the thoracic spine in patients with AIS. This study aimed to investigate the mechanical loading and its distribution on the thoracic facet joints of AIS patients using finite element (FE) analysis and surgical specimens. FE models of the thoracic spine were created from a total of 13 female AIS patients (Lenke type 1, n = 4; Lenke type 2, n = 4; Lenke type 3, n = 5). A load of 200 N on the T3 vertebrae and 30 N each on the bilateral superior articular processes were applied vertically to quantify the contact force on the facet joints from T3 to T11. In addition, morphological and histological analyses were performed on the inferior articular processes obtained during surgery. FE analysis demonstrated that contact forces of the facet joint progressively increased from the mid to lower thoracic spine of the concave side, reaching a maximum around the apex. More than 91% of the load was transmitted by the facet joints at the concave side, resulting in facet joint subchondral sclerosis and hypertrophy. The apical facet joint in AIS helps counteract rotational stress between vertebrae and transfers most stress through the concave side. In conclusion, this study found that asymmetric load transfer in the facet joints leads to subchondral sclerosis and hypertrophy. These findings can enhance our understanding of the stress loading on facet joints and the resulting biological changes and help clarify the mechanisms involved in scoliosis progression. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Biomechanical effects of osteoporosis severity on the occurrence of proximal junctional kyphosis following long-segment posterior thoracolumbar fusion

BACKGROUND

Proximal junctional kyphosis is a common long-term complication in adult spinal deformity surgery that involves long-segment posterior spinal fusion. However, the underlying biomechanical mechanisms of the impact of osteoporosis on proximal junctional kyphosis remain unclear. The present study was to evaluate adjacent segment degeneration and spine mechanical instability in osteoporotic patients who underwent long-segment posterior thoracolumbar fusion.

METHODS

Finite element models of the thoracolumbar spine T1-L5 with posterior long-segment T8-L5 fusion under different degrees of osteoporosis were constructed to analyze intervertebral disc stress characterization, vertebrae mechanical transfer, and pedicle screw system loads during various motions.

FINDINGS

Compared with normal bone mass, the maximum von Mises stresses of T7 and T8 were increased by 20.32%, 22.38%, 44.69%, 4.49% and 29.48%, 17.84%, 40.95%, 3.20% during flexion, extension, lateral bending, and axial rotation in the mild osteoporosis model, and by 21.21%, 18.32%, 88.28%, 2.94% and 37.76%, 15.09%, 61.47%, -0.04% in severe osteoporosis model. The peak stresses among T6/T7, T7/T8, and T8/T9 discs were 14.77 MPa, 11.55 MPa, and 2.39 MPa under lateral bending conditions for the severe osteoporosis model, respectively. As the severity of osteoporosis increased, stress levels on SCR8 and SCR9 intensified during various movements.

INTERPRETATION

Osteoporosis had an adverse effect on proximal junctional kyphosis. The stress levels in cortical bone, intervertebral discs and screws were increased with bone mass loss, which can easily lead to intervertebral disc degeneration, bone destruction as well as screw pullout. These factors have significantly affected or accelerated the occurrence of proximal junctional kyphosis.

Oblique lateral interbody fusion with internal fixations in the treatment for cross-segment degenerative lumbar spine disease (L2-3 and L4-5) finite element analysis

Multi-segmental lumbar degenerative disease, including intersegmental disc degeneration, is found in clinical practice. Controversy still exists regarding the treatment for cross-segment degeneration. Oblique Lateral Interbody Fusion (OLIF) with several internal fixations was used to treat cross-segment lumbar degenerative disease. A whole lumbar spine model was extracted from CT images of the whole lumbar spine of patients with lumbar degeneration. The L2-3 and L4-5 intervertebral spaces were fused with OLIF using modeling software, the Pedicle screws were performed on L2-3 and L4-5, and different internal fixations were performed on L3-4 in Finite Element (FE) software. Among the six 10 Nm moments of different directions, the L3-4 no surgery (NS) group had the relatively largest Range of Motion (ROM) in the whole lumbar spine, while the L2-5 Long segmental fixation (LSF)group had the smallest ROM and the other groups had similar ROM. The ROM in the L1-2 and L5-S1 was relatively close in the six group models, and the articular cartilage stress and disc stress on the L1-2 and L5-S1 were relatively close. In contrast, the L3-4 ROM differed relatively greatly, with the LSF ROM the smallest and the NS ROM the largest, and the L3-4 Coflex (Coflex) group more active than the L3-4 Bacfuse (Bacfuse) group and the L3-4 translaminar facet screw fixation (TFSF) group. The stress on the articular cartilage and disc at L3-4 was relatively greater in the NS disc and articular cartilage, and greater in the Coflex group than in the Bacfuse and TFSF groups, with the greatest stress on the internal fixation in the TFSF group, followed by the Coflex group, and relatively similar stress in the Bacfuse, LSF, and NS groups. In the TFSF group, the stress on the internal fixation was greater than the yield strength among different directional moments of 10 Nm, which means it is unsuitable to be an internal fixation. The LSF group had the greatest overall ROM, which may lead to postoperative low back discomfort. The NS group has the greatest overall ROM, but its increased stress on the L3-4 disc and articular cartilage may lead to accelerated degeneration of the L3-4 disc and articular cartilage. The Coflex and Bacfuse groups had a reduced L3-4 ROM but a greater stress on disc compared to the LSF group, which may lead to disc degeneration in the long term. However, their stress on the articular cartilage was relatively low. Coflex and Bacfuse can still be considered better surgical options.

Comparative Efficacy of Vibration foam Rolling and Cold Water Immersion in Amateur Basketball Players after a Simulated Load of Basketball Game

To compare the efficacy of different recovery strategies (sitting; cold water immersion, CWI; vibration foam rolling, VFR) on the lower extremities of amateur basketball players after the simulated load of a basketball game, we assessed the power, agility, and dynamic balance before and after interventions. Ten amateur basketball players alternately underwent 12 min of sitting, 12 min of CWI at 5 °C, and 12 min of VFR. The power, agility, and dynamic balance were measured immediately post-warm-up, immediately post-game, immediately post-intervention, 1 h after interventions, and 24 h after interventions. To simulate the load of a basketball game, specific movements were designed and implemented. Jump height was measured using a Kistler force plate. Reaction time and dynamic balance score were assessed using the Pavigym agility response system and the Y balance test, respectively. The data were analyzed with a two-way repeated measures analysis of variance (ANOVA). The results showed that the vertical jump height significantly decreased after the CWI intervention compared to the CON and VFR groups (p < 0.001). At 1 h after the intervention, the vertical jump height in the CON group showed delayed recovery compared to the CWI and VFR groups (p = 0.007; p < 0.001). At 24 h after the intervention, the vertical jump height in the CWI group further increased and was significantly different from the CON and VFR groups (p < 0.001; p = 0.005). Additionally, reaction times significantly increased immediately after the CWI intervention (p = 0.004) but showed further recovery at 24 h compared to the CON group (p < 0.001). The dynamic balance score significantly rebounded after the CWI intervention compared to the CON group (p = 0.021), with further improvement at 24 h (p < 0.001). CWI initially showed negative effects, but over time, its recovery effect was superior and more long-lasting. VFR had the best immediate effect on lower limb recovery after the game.

Applications of Finite Element Modeling in Biomechanical Analysis of Foot Arch Deformation: A Scoping Review

Excessive foot arch deformation is associated with plantar tissue overload and ligamentous injury pathologies. Finite element (FE) analysis, as an effective tool for modeling and simulation, has been utilized clinically for providing insights into arch biomechanics. This systematic scoping review aimed to summarize the current state of computational modeling techniques utilized in arch biomechanics from 2000 onwards and outline the main challenges confronting the further development of accurate models in clinical conditions. English-language searches of the electronic databases were conducted in the Web of Science, PubMed, and Scopus until July 2022. Articles that investigated arch deformation mechanisms by FE modeling were included. The methodological quality was assessed utilizing the Methodological Quality Assessment of Subject-Specific Finite Element Analysis Used in Computational Orthopedics (MQSSFE). Seventeen articles were identified in this systematic scoping review, mostly focusing on constructing models for specific pathological conditions, such as progressive collapsing foot deformity, valgus foot, and posterior tibial tendon dysfunction. However, given the complexity of the arch problem, geometrical simplifications regarding the balance between accurate detail and computational cost and assumptions made in defining modeling parameters (material properties and loading and boundary conditions) may bring challenges to the accuracy and generalizability of models applied to clinical settings. Overall, advances in computational modeling techniques have contributed to reliable foot deformation simulation and analysis in modern personalized medicine.

Effect of pregnancy on female gait characteristics: a pilot study based on portable gait analyzer and induced acceleration analysis

Introduction: The changes in physical shape and center of mass during pregnancy may increase the risk of falls. However, there were few studies on the effects of maternal muscles on gait characteristics and no studies have attempted to investigate changes in induced acceleration during pregnancy. Further research in this area may help to reveal the causes of gait changes in women during pregnancy and provide ideas for the design of footwear and clothing for pregnant women. The purpose of this study is to compare gait characteristics and induced accelerations between non-pregnant and pregnant women using OpenSim musculoskeletal modeling techniques, and to analyze their impact on pregnancy gait. Methods: Forty healthy participants participated in this study, including 20 healthy non-pregnant and 20 pregnant women (32.25 ± 5.36 weeks). The portable gait analyzer was used to collect participants' conventional gait parameters. The adjusted OpenSim personalized musculoskeletal model analyzed the participants' kinematics, kinetics, and induced acceleration. Independent sample T-test and one-dimensional parameter statistical mapping analysis were used to compare the differences in gait characteristics between pregnant and non-pregnant women. Results: Compared to the control group, pregnancy had a 0.34 m reduction in mean walking speed (p < 0.01), a decrease in mean stride length of 0.19 m (p < 0.01), a decrease in mean stride frequency of 19.06 step/min (p < 0.01), a decrease in mean thigh acceleration of 0.14 m/s² (p < 0.01), a decrease in mean swing work of 0.23 g (p < 0.01), and a decrease in mean leg falling strength of 0.84 g (p < 0.01). Induced acceleration analysis showed that pregnancy muscle-induced acceleration decreased in late pregnancy (p < 0.01), and the contribution of the gastrocnemius muscle to the hip and joint increased (p < 0.01). Discussion: Compared with non-pregnant women, the gait characteristics, movement amplitude, and joint moment of pregnant women changed significantly. This study observed for the first time that the pregnant women relied more on gluteus than quadriceps to extend their knee joints during walking compared with the control group. This change may be due to an adaptive change in body shape and mass during pregnancy.

Simulation of Lower Limb Muscle Activation Using Running Shoes with Different Heel-to-Toe Drops Using Opensim

BACKGROUND

Although numerous studies have been conducted to investigate the acute effects of shoe drops on running kinematics and kinetic variables, their effects on muscle forces remain unknown. Thus, the primary aim of this study was to compare the muscle force, kinematics, and kinetic variables of habitually rearfoot runners with heel-to-toe drops of negative 8 mm shoes (minimalist shoes) and positive 9 mm shoes (normal shoes) during the running stance phase by using musculoskeletal modeling and simulation techniques.

METHODS

Experimental data of lower limb kinematics, ground reaction force, and muscle activation from 16 healthy runners with rearfoot strike patterns were collected and analyzed in OpenSim. Using Matlab, the statistical parameter mapping paired t-test was used to compare the joint angle, moment, and muscle force waveform.

RESULTS

The results revealed differences in the sagittal ankle and hip angles and sagittal knee moments between the different heel-to-toe drops of running shoes. Specifically, it showed that the negative 8 mm running shoes led to significantly smaller values than the positive 9 mm running shoes in terms of the angle of ankle dorsiflexion, ankle eversion, knee flexion, hip flexion, and hip internal and hip external rotation. The peak ankle dorsiflexion moment, ankle plantarflexion moment, ankle eversion moment, knee flexion moment, knee abduction moment, and knee internal rotation also decreased obviously with the minimalist running shoes, while the lateral gastrocnemius, Achilleas tendon, and extensor hallucis longus muscles were obviously greater in the minimalist shoes compared to normal shoes. The vastus medialis, vastus lateralis and extensor digitorum longus muscles force were smaller in the minimalist shoes.

CONCLUSIONS

Runners may shift to a midfoot strike pattern when wearing negative running shoes. High muscle forces in the gastrocnemius lateral, Achilleas tendon, and flexor hallucis longus muscles may also indicate an increased risk of Achilleas tendonitis and ankle flexor injuries.

Biomechanics of Topspin Forehand Loop in Table Tennis: An Application of OpenSim Musculoskeletal Modelling

Topspin is one of the most attacking strokes in table tennis, and topspin forehand loop is an effective way to score. The aim of this study was to investigate the kinematics of the lower extremities in topspin forehand loop between different levels via OpenSim Musculoskeletal Modelling. Ten elite athletes (NL1) and ten medium athletes (NL2) performed the topspin forehand loop without muscle and joint injuries. An eight-camera Vicon motion capture system was used to measure the kinematics data. During the topspin forehand loop, the forward phase (FP) and the entire phase (EP) of the NL1 were significantly shorter than that of the NL2. In the sagittal plane, NL1 significantly had greater hip and ankle flexion and extension at range of motion (ROM) but less hip flexion and knee flexion at FP and less ankle flexion at BP than NL2. In the frontal plane, NL1 displayed less ROM in the hip joint and significantly less hip abduction ROM at the backward phase (BP). In the transverse plane, NL1 had a significantly greater ROM in the hip joint and displayed significantly less hip ROM at the BP. The level differences presented in this study could help table tennis athletes to improve performance and coaches to develop technical training.

Effect of Interbody Implants on the Biomechanical Behavior of Lateral Lumbar Interbody Fusion: A Finite Element Study

Porous titanium interbody scaffolds are growing in popularity due to their appealing advantages for bone ingrowth. This study aimed to investigate the biomechanical effects of scaffold materials in both normal and osteoporotic lumbar spines using a finite element (FE) model. Four scaffold materials were compared: Ti6Al4V (Ti), PEEK, porous titanium of 65% porosity (P65), and porous titanium of 80% porosity (P80). In addition, the range of motion (ROM), endplate stress, scaffold stress, and pedicle screw stress were calculated and compared. The results showed that the ROM decreased by more than 96% after surgery, and the solid Ti scaffold provided the lowest ROM (1.2-3.4% of the intact case) at the surgical segment among all models. Compared to solid Ti, PEEK decreased the scaffold stress by 53-66 and the endplate stress by 0-33%, while porous Ti decreased the scaffold stress by 20-32% and the endplate stress by 0-32%. Further, compared with P65, P80 slightly increased the ROM (<0.03°) and pedicle screw stress (<4%) and decreased the endplate stress by 0-13% and scaffold stress by approximately 18%. Moreover, the osteoporotic lumbar spine provided higher ROMs, endplate stresses, scaffold stresses, and pedicle screw stresses in all motion modes. The porous Ti scaffolds may offer an alternative for lateral lumbar interbody fusion.

Three-Dimensional Computational Model Simulating the Initial Callus Growth during Fracture Healing in Long Bones: Application to Different Fracture Types

Bone fractures are among the most common and potentially serious injuries to the skeleton, femoral shaft fractures being especially severe. Thanks to recent advances in the area of in silico analysis, several approximations of the bone healing process have been achieved. In this context, the objective of this work was to simulate the initial phase of callus formation in long bones, without a pre-meshed domain in the 3D space. A finite element approach was computationally implemented to obtain the values of the cell concentrations along the whole domain and evaluate the areas where the biological quantities reached the thresholds necessary to trigger callus growth. A voxel model was used to obtain the 3D domain of the bone fragments and callus. A mesh growth algorithm controlled the addition of new elements to the domain at each step of the iterative procedure until complete callus formation. The implemented approach is able to reproduce the generation of the primary callus, which corresponds to the initial phase of fracture healing, independently of the fracture type and complexity, even in the case of several bone fragments. The proposed approach can be applied to the most complex bone fractures such as oblique, severely comminuted or spiral-type fractures, whose simulation remains hardly possible by means of the different existing approaches available to date.

The Lower Limb Movements of the Fetus in Uterus: A Narrative Review

The fetus movements play an important role in fetal well-being. With the continuous advancement of real-time scanning machines, it is feasible to observe the fetus movement in detail. The characteristics of fetal lower limb movements in prenatal examination have not been systematically investigated. This review proposes the patterns of fetal lower limb movements, the maternal influence on fetal lower limb movements, and the application of fetal lower limb movements for the diagnosis of prenatal diseases. A systematic search of literature on the lower limb movements of the fetus in uterus was performed in the databases, namely, Web of Science and Scopus. Thirty-four publications were selected. This review demonstrates that isolated fetal lower limb movements are rare and always accompanied with the movements of other body segments. Detection of the presence of fetal leg movements seems to be of no diagnostic value for fetuses with prenatal diseases. The isolated lower limb movement was statistically significant different between fetuses of low- and high-risk pregnant women. The coordinated movements of the fetal lower limbs and other parts should be considered when analyzing fetal movements in the future study.

Material sensitivity of patient-specific finite element models in the brace treatment of scoliosis

Objectives: To study the mechanical sensitivity of different intervertebral disc and bone material parameters and ligaments under different force configurations and magnitudes in the scoliosis model. Methods: The finite element model of a 21-year-old female is built using computed tomography. Local range of motion testing and global bending simulations are performed for the model verification. Subsequently, Five force of different directions and configurations were applied to the finite element model applying the brace pad position. The material parameters of the model were related to different spinal flexibilities and included different material parameters of cortical bone, cancellous bone, nucleus and annulus. The virtual X-ray technique measured Cobb angle, thoracic Lordosis, and lumbar Kyphosis. Results: The difference in peak displacement is 9.28 mm, 19.99 mm, 27.06 mm, 43.99 mm, and 50.1 mm under five force configurations. The maximum Cobb angle difference due to material parameters are 4.7° and 6.2°, which are converted to thoracic and lumbar in-brace correction difference of 18% and 15.5%. The maximum difference in Kyphosis and Lordosis angle is 4.4° and 5.8°. The average thoracic and lumbar Cobb angle variation difference in intervertebral disc control group is larger than that in bone control group, while the average Kyphosis and Lordosis angle is inverse. The displacement distribution of models with or without ligaments is similar, with a peak displacement difference of 1.3 mm in C5. The peak stress occurred at the junction of the cortical bone and ribs. Conclusion: Spinal flexibility largely influences the treatment effect of the brace. The intervertebral disc has a greater effect on the Cobb angle, the bone has a greater effect on the Kyphosis and Lordosis angles, and the rotation is affected by both. Patient-specific material is the key to increasing accuracy in the personalized finite element model. This study provides a scientific basis for using controllable brace treatment for scoliosis.

Automated recognition of asymmetric gait and fatigue gait using ground reaction force data

Introduction: The purpose of this study was to evaluate the effect of running-induced fatigue on the characteristic asymmetry of running gait and to identify non-linear differences in bilateral lower limbs and fatigued gait by building a machine learning model. Methods: Data on bilateral lower limb three-dimensional ground reaction forces were collected from 14 male amateur runners before and after a running-induced fatigue experiment. The symmetry function (SF) was used to assess the degree of symmetry of running gait. Statistical parameter mapping (Paired sample T-test) algorithm was used to examine bilateral lower limb differences and asymmetry changes pre- and post-fatigue of time series data. The support vector ma-chine (SVM) algorithm was used to recognize the gait characteristics of both lower limbs before and after fatigue and to build the optimal algorithm model by setting different kernel functions. Results: The results showed that the ground reaction forces were asymmetrical (SF > 0.5) both pre-and post-fatigue and mainly concentrated in the medial-lateral direction. The asymmetry of the medial-lateral direction increased significantly after fatigue (p < 0.05). In addition, we concluded that the polynomial kernel function could make the SVM model the most accurate in classifying left and right gait features (accuracy of 85.3%, 82.4%, and 82.4% in medial-lateral, anterior-posterior and vertical directions, respectively). Gaussian radial basis kernel function was the optimal kernel function of the SVM algorithm model for fatigue gait recognition in the medial-lateral and vertical directions (accuracy of 54.2% and 62.5%, respectively). Moreover, polynomial was the optimal kernel function of the anterior-posterior di-rection (accuracy = 54.2%). Discussion: We proved in this study that the SVM algorithm model depicted good performance in identifying asymmetric and fatigue gaits. These findings can provide implications for running injury prevention, movement monitoring, and gait assessment.

The Effect of Concave-Side Intertransverse Ligament Laxity on the Stress of AIS Lumbar Spine Based on Finite Element Method

(1) Background: Scoliosis has the mechanical characteristic of asymmetric stress distribution, which is one of the reasons for the aggravation of scoliosis. Bracing therapy is the best treatment for AIS, but it is difficult and costly to operate. Is it possible to reduce pressure in the concave side by relaxing the ITL in the concave side of scoliosis, so as to improve the abnormal stress distribution of scoliosis? In this paper, a finite element method was used to simulate the effect of the relaxation of concave-side ITL on the stress of a lumbar spine with scoliosis, which provides some guidance for the treatment of scoliosis. (2) Methods: Using CT images of a patient with scoliosis whose Cobb Angle was 43° and Lordosis Angle was 45, a scoliosis lumbar was established, and Young's modulus of the ITL of the concave-side lumbar spine was reduced by 95% to simulate ligament relaxation. By comparing the stress condition of the model vertebral body with no ligament relaxation, the effect of concave-side ITL relaxation on the mechanical characteristics of scoliosis lumbar spine was explored. (3) Results: An effective and complete model of the lumbar spine was established. The concave ITL relaxed, which only had a great impact on the bending loads. After the ligament was relaxed, the stability of the spine was reduced. Stress concentration on the concave side of vertebrae and the IVD was aggravated. Under loads on the convex side, the maximum stress on the vertebral body and the IVD increased significantly, making lumbar vertebrae more vulnerable to injury. (4) Conclusions: Laxity of the ITL on the concave side of the AIS lumbar only affects the bending load. Laxity of the concave-side ligament will reduce the stability of the lumbar, aggravate the uneven stress distribution of scoliotic lumbar vertebrae, increase the risk of IVD injury, and be unfavorable for the scoliotic lumbar spine. Relaxation of the concave ITL alone is not an effective way to treat scoliosis.

A new method proposed to explore the feline's paw bones of contributing most to landing pattern recognition when landed under different constraints

Felines are generally acknowledged to have natural athletic ability, especially in jumping and landing. The adage “felines have nine lives” seems applicable when we consider its ability to land safely from heights. Traditional post-processing of finite element analysis (FEA) is usually based on stress distribution trend and maximum stress values, which is often related to the smoothness and morphological characteristics of the finite element model and cannot be used to comprehensively and deeply explore the mechanical mechanism of the bone. Machine learning methods that focus on feature pattern variable analysis have been gradually applied in the field of biomechanics. Therefore, this study investigated the cat forelimb biomechanical characteristics when landing from different heights using FEA and feature engineering techniques for post-processing of FEA. The results suggested that the stress distribution feature of the second, fourth metacarpal, the second, third proximal phalanx are the features that contribute most to landing pattern recognition when cats landed under different constraints. With increments in landing altitude, the variations in landing pattern differences may be a response of the cat's forelimb by adjusting the musculoskeletal structure to reduce the risk of injury with a more optimal landing strategy. The combination of feature engineering techniques can effectively identify the bone's features that contribute most to pattern recognition under different constraints, which is conducive to the grasp of the optimal feature that can reveal intrinsic properties in the field of biomechanics.