Patient-specific Finite Element Models Research Articles

Effects of bone mineral density at the lateral sclerotic boundary on the femoral head collapse onset in osteonecrosis of the femoral head: A preliminary study

BackgroundIn the natural course of osteonecrosis of the femoral head, sclerotic changes at the boundary of necrotic lesion gradually occur until femoral head collapse. This study aims to examine the effects of bone mineral density at the lateral boundary of necrotic lesion on a subsequent femoral head collapse. MethodsWe developed patient-specific finite element models of 9 hips with subsequent collapse and 10 hips without subsequent collapse. Cubic regions of interest were selected at both subchondral areas of the lateral boundary and the adjacent necrotic lesion. Bone mineral density values of the regions of interest were quantitatively measured, and a ratio of bone mineral density values (lateral boundary/necrotic lesion) was calculated. Stress values at the lateral boundary were also evaluated. FindingsThe ratio of bone mineral density values was significantly higher in hips with subsequent collapse than that without subsequent collapse (p = 0.0016). The median equivalent stress and shear stress were significantly higher in hips with subsequent collapse than that without subsequent collapse (p = 0.0071, and p = 0.0143, respectively). The ratio of bone mineral density values showed a promising value in predicting the occurrence of subsequent femoral head collapse (AUC = 0.97). InterpretationOur results indicated that bone mineral density value at the lateral boundary of necrotic lesion may be associated with the occurrence of subsequent femoral head collapse in pre-collapse stage osteonecrosis of the femoral head.

A novel method for assigning bone material properties to a comprehensive patient-specific pelvic finite element model using biplanar multi-energy radiographs

The increasing prevalence of adult spinal deformity requires long spino-pelvic instrumentation, but pelvic fixation faces challenges due to distal forces and reduced bone quality. Bi-planar multi-energy X-rays (BMEX) were used to develop a patient-specific finite element model (FEM) for evaluating pelvic fixation. Calibration involved 10 patients, and an 81-year-old female test case was used for FEM customization and pullout simulation validation. Calibration yielded a root mean square error of 74.7 mg/cm3 for HU. The simulation accurately replicated the experimental pullout test with a force of 565 N, highlighting the method’s potential for optimizing biomechanical performance for pelvic fixation.

Assessment of Intramural Segment Compression in Anomalous Coronary Arteries through Patient-Specific Finite Element Modeling

Background: Anomalous Aortic Origin of a Coronary Artery (AAOCA) is a congenital condition that can lead to ischemia and sudden cardiac death. Current diagnostic tools are unable to fully quantify the pathological behavior that occurs mainly with physical effort. Methods: Patients’ computed tomography scans and centerline-based geometric quantities were used to develop three-dimensional computer-aided design models of the main anatomical variants of AAOCA. Blood pressure ranging from rest to extreme effort was simulated through structural finite element analyses, and the pressurized geometries were analyzed to evaluate coronary lumen cross-sectional areas and variations at the different loading conditions. Results: We simulated 39 subjects, demonstrating the ability to reproduce accurately the patient-specific anatomy of different AAOCA variants and capture pathological behaviors. AAOCAs with intramural courses showed compression along the proximal segment with a caliber reduction ranging from 0.14% to 18.87% at different pressure levels. The percentage of proximal narrowing relative to the distal segment was greater than any other type of anomalous course and exceeded 50% with simulated exertion. Conclusions: The present study proposes a computational pipeline to investigate conditions not reproducible in clinical practice, providing information to support decision-making in the management of AAOCA patients.

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.

A novel computer modeling and simulation technique for bronchi motion tracking in human lungs under respiration

In this work, we proposed a novel computer modeling and simulation technique for motion tracking of lung bronchi (or tumors) under respiration using 9 cases of computed tomography (CT)-based patient-specific finite element (FE) models and Ogden’s hyperelastic model. In the fabrication of patient-specific FE models for the respiratory system, various organs such as the mediastinum, diaphragm, and thorax that could affect the lung motions during breathing were considered. To describe the nonlinear material behavior of lung parenchyma, the comparative simulation for biaxial tension-compression of lung parenchyma was carried out using several hyperelastic models in ABAQUS, and then, Ogden’s model was adopted as an optimal model. Based on the aforementioned FE models and Ogden’s material model, the 9 cases of respiration simulation were carried out from exhalation to inhalation, and the motion of lung bronchi (or tumors) was tracked. In addition, the changes in lung volume, lung cross-sectional area on the axial plane during breathing were calculated. Finally, the simulation results were quantitatively compared to the inhalation/exhalation CT images of 9 subjects to validate the proposed technique. Through the simulation, it was confirmed that the average relative errors of simulation to clinical data regarding to the displacement of 258 landmarks in the lung bronchi branches of total subjects were 1.10%~2.67%. In addition, the average relative errors of those with respect to the lung cross-sectional area changes and the volume changes in the superior-inferior direction were 0.20%~5.00% and 1.29 ~ 9.23%, respectively. Hence, it was considered that the simulation results were coincided well with the clinical data. The novelty of the present study is as follows: (1) The framework from fabrication of the human respiratory system to validation of the bronchi motion tracking is provided step by step. (2) The comparative simulation study for nonlinear material behavior of lung parenchyma was carried out to describe the realistic lung motion. (3) Various organs surrounding the lung parenchyma and restricting its motion were considered in respiration simulation. (4) The simulation results such as landmark displacement, lung cross-sectional area/volume changes were quantitatively compared to the clinical data of 9 subjects.

Determinants of spinal cord stress and strain in degenerative cervical myelopathy: a patient-specific finite element study.

Degenerative cervical myelopathy (DCM) is the commonest cause of spinal cord dysfunction in older adults and is characterized by chronic cervical spinal cord compression. Spinal cord stress and strain during neck motion are also known contributors to the pathophysiology of DCM, yet these factors are not routinely assessed for surgical planning. The aim of this study was to measure spinal cord stress/strain in DCM using patient-specific 3D finite element models (FEMs) and determine whether spinal cord compression is the primary determinant of spinal cord stress/strain. Three-dimensional patient-specific FEMs were created for six DCM patients (mild [n = 2], moderate [n = 2] and severe [n = 2]). Flexion and extension of the cervical spine were simulated with a pure moment load of 2 Nm. Segmental spinal cord von Mises stress and maximum principal strain were measured. Measures of spinal cord compression and segmental range of motion (ROM) were included in a regression analysis to determine associations with spinal cord stress and strain. Segmental ROM in flexion-extension and axial rotation was independently associated with spinal cord stress (p < 0.001) and strain (p < 0.001), respectively. This relationship was not seen for lateral bending. Segmental ROM had a stronger association with spinal stress and strain as compared to spinal cord compression. Compared to the severity of spinal cord compression, segmental ROM is a stronger determinant spinal cord stress and strain. Surgical procedures that address segmental ROM in addition to cord compression may best optimize spinal cord biomechanics in DCM.

Biomechanical modeling and assessment of lumbar vertebral body tethering configurations.

Vertebral body tethering (VBT) is a fusionless spinal growth modulation technique, which shows promise for pediatric idiopathic scoliosis (IS) curve correction. This technique, mainly used for thoracic curves, is increasingly being used to treat lumbar curves in order to preserve spine flexibility. It remains necessary to adequately define the cord tension to be applied during the operation and the instrumented levels to biomechanically predict correction over time for the lumbar spine. Twelve pediatric patients with lumbar IS, treated with lumbar-only or lumbar and thoracic VBT, were selected for this study. Three independent variables were tested alternately using a patient-specific finite element model (FEM), which includes an algorithm modeling vertebra growth and spinecurvechanges due to growth modulation for 24months post-operatively according to the Hueter-Volkmann principle. Parameters included cable tensioning (150N/250N), upper instrumented level (actual UIV, UIV-1) and lower instrumented level (actual LIV, LIV + 1). Each FEM was personalized using 3D radiographic reconstruction and flexibility supine radiographs. An increase in cord tension (from 150 to 250N) had significant effects on main thoracic and thoraco-lumbar/lumbar Cobb angles, as well as on lumbar lordosis, after surgery (supplementary average correction of 3° and 8°, and increase of 1.4°, respectively) and after 24months (4°, 10° and 1.1°) (p < 0.05). Adding a level to the actual UIV or LIV did not improve correction. This parametric study showed that cord tension is the most important biomechanical parameter on the simulated immediate and 2-year increase in lumbar curve correction. Our preliminary model suggests that it is not advantageous to add additional instrumented levels. This computational study uses a retrospective validation cohort (level of evidence 3).

Design and biomechanical analysis of patientspecific porous tantalum prostheses for knee joint revision surgery.

Artificial joint revision surgery, as an increasingly common surgery in orthopedics, often requires patient-specific prostheses to repair the bone defect. Porous tantalum is a good candidate due to its excellent abrasion and corrosion resistance and good osteointegration. Combination of 3D printing technology and numerical simulation is a promising strategy to design and prepare patient-specific porous prostheses. However, clinical design cases have rarely been reported, especially from the viewpoint of biomechanical matching with the patient's weight and motion and specific bone tissue. This work reports a clinical case on the design and mechanical analysis of 3D-printed porous tantalum prostheses for the knee revision of an 84-year-old male patient. Particularly, standard cylinders of 3D-printed porous tantalum with different pore size and wire diameters were first fabricated and their compressive mechanical properties were measured for following numerical simulation. Subsequently, patientspecific finite element models for the knee prosthesis and the tibia were constructed from the patient's computed tomography data. The maximum von Mises stress and displacement of the prostheses and tibia and the maximum compressive strain of the tibia were numerically simulated under two loading conditions by using finite element analysis software ABAQUS. Finally, by comparing the simulated data to the biomechanical requirements for the prosthesis and the tibia, a patient-specific porous tantalum knee joint prosthesis with a pore diameter of 600 μm and a wire diameter of 900 μm was determined. The Young's modulus (5719.32 ± 100.61 MPa) and yield strength (172.71 ± 1.67 MPa) of the prosthesis can produce both sufficient mechanical support and biomechanical stimulation to the tibia. This work provides a useful guidance for designing and evaluating a patient-specific porous tantalum prosthesis.

RETHINKING THE 10% RULE: PATIENT-SPECIFIC MODELLING OF STRAIN IN DISTAL FEMUR FRACTURES

The objective of this study was to use patient-specific finite element modeling to measure the 3D interfragmentary strain environment in clinically realistic fractures. The hypothesis was that in the early post-operative period, the tissues in and around the fracture gap can tolerate a state of strain in excess of 10%, the classical limit proposed in the Perren strain theory.Eight patients (6 males, 2 females; ages 22–95 years) with distal femur fractures (OTA/AO 33-A/B/C) treated in a Level I trauma center were retrospectively identified. All were treated with lateral bridge plating. Preoperative computed tomography scans and post-operative X-rays were used to create the reduced fracture models. Patient-specific materials properties and loading conditions (20%, 60%, and 100% body weight (BW)) were applied following our published method.[1]Elements with von Mises strains &gt;10% are shown in the 100% BW loading condition. For all three loading scenarios, as the bridge span increased, so did the maximum von Mises strain within the strain visualization region. The average gap closing (Perren) strain (mean ± SD) for all patient-specific models at each body weight (20%, 60%, and 100%) was 8.6% ± 3.9%, 25.8% ± 33.9%, and 39.3% ± 33.9%, while the corresponding max von Mises strains were 42.0% ± 29%, 110.7% ± 32.7%, and 168.4% ± 31.9%. Strains in and around the fracture gap stayed in the 2–10% range only for the lowest load application level (20% BW).Moderate loading of 60% BW and above caused gap strains that far exceeded the upper limit of the classical strain rule (&lt;10% strain for bone healing). Since all of the included patients achieved successful unions, these findings suggest that healing of distal femur fractures may be robust to localized strains greater than 10%.

Differences in spinal cord biomechanics after laminectomy, laminoplasty, and laminectomy with fusion for degenerative cervical myelopathy.

Spinal cord stress/strain during neck motion contributes to spinal cord dysfunction in degenerative cervical myelopathy (DCM), yet the effect of surgery on spinal cord biomechanics is unknown. It is expected that motion-preserving and fusion surgeries for DCM will have distinct effects on spinal cord biomechanics. The aim of this study was to compare changes in spinal cord biomechanics after laminectomy with fusion, laminectomy, and laminoplasty using a patient-specific finite element model (FEM) for DCM. A patient-specific FEM of the cervical spine and spinal cord was created using MRI from a subject with mild DCM. Multilevel laminectomy with fusion, laminectomy, and laminoplasty were simulated for DCM using the patient-specific FEM. Spinal cord von Mises stress and maximum principal strain during neck flexion-extension, lateral bending, and axial rotation were recorded. Segmental range of motion, intradiscal pressure, and capsular ligament strain were also measured. FEM outputs were calculated as a change with respect to the preoperative values and compared between the three models. Across the surgical levels, spinal cord stress increased after laminectomy for neck flexion (+50%), neck extension (+37.8%), and axial rotation (+23%). Similarly, spinal cord strain increased in neck extension (+118.4%) and axial rotation (+75.1%) after laminectomy. Laminoplasty was associated with greater spinal cord stress in neck flexion (+57.4%) and increased strain in lateral bending (+56.7%) and axial rotation (+20.9%). Compared with laminectomy and laminoplasty, spinal cord biomechanics for laminectomy with fusion revealed significantly reduced median extension stress (13.7 kPa vs 9.7 kPa, p = 0.03), lateral bending strain (0.01 vs 0.007, p = 0.007), axial rotation stress (3.7 kPa vs 2.1 kPa, p = 0.04), and axial rotation strain (0.017 vs 0.009, p = 0.04). Spinal cord strain decreased in neck flexion in all three models, yet spinal cord stress increased with neck flexion for laminectomy and laminoplasty. Changes in spinal cord biomechanics for laminoplasty parallel those for laminectomy with fusion except during neck flexion, lateral bending, and axial rotation. Compared with motion-preserving approaches such as laminectomy and laminoplasty, laminectomy with fusion was associated with the lowest spinal cord stress and strain in flexion-extension, lateral bending, and axial rotation of the neck.

Biomechanical evaluation of a healed acetabulum with internal fixators: finite element analysis

BackgroundTreatment of complicated acetabular fracture with internal fixation usually has high risk of failure because of unbefitting fixation. However, evaluation of the biomechanical effect of internal fixation under physiological loading for fracture healing is still generally rarely performed. The purpose of this study is to analyze the biomechanical characteristics of a healed acetabulum with designed internal fixators under gait and to explore the biomechanical relationship between the healed bone and the internal fixator.MethodsA patient-specific finite element model of whole pelvis with designed internal fixators was constructed based on the tomographic digital images, in which the spring element was used to simulate the main ligaments of the pelvis. And the finite element analysis under both the combination loading of different phases and the individual loading of each phase during the gait cycle was carried out. The displacement, von Mises stress, and strain energy of both the healed bone and the fixation were calculated to evaluate the biomechanical characteristics of the healed pelvis.ResultsUnder the combination loading of gait, the maximum difference of displacement between the left hip bone with serious injury and the right hip bone with minor injury is 0.122 mm, and the maximum stress of the left and right hemi-pelvis is 115.5 MPa and 124.28 MPa, respectively. Moreover, the differences of average stress between the bone and internal fixators are in the range of 2.3–13.7 MPa. During the eight phases of gait, the stress distribution of the left and right hip bone is similar. Meanwhile, based on the acetabular three-column theory, the strain energy ratio of the central column is relatively large in stance phases, while the anterior column and posterior column of the acetabular three-column increase in swing phases.ConclusionsThe acetabular internal fixators designed by according to the anatomical feature of the acetabulum are integrated into the normal physiological stress conduction of the pelvis. The design and placement of the acetabular internal fixation conforming to the biomechanical characteristics of the bone is beneficial to the anatomical reduction and effective fixation of the fracture, especially for complex acetabular fracture.

Finite element simulation of growth modulation during brace treatment of adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a spine deformity whose progression during growth is affected by asymmetrical loads acting on the spine. The conservative brace treatment aims to limit the deformity's progression until the end of skeletal growth. This study's objective was to develop a patient-specific finite element model (FEM) simulating immediate in-brace correction and subsequent growth modulation over 2 years of treatment. Thirty-five retrospective AIS cases with documented correction over 2 years were analyzed. For each case, a patient-specific FEM was built, and in-brace correction was simulated. Vertebral growth and its modulation were modeled using simulated pressures on epiphyseal vertebral growth plates, including a compliance factor representing the recorded brace wear. The simulated Cobb angles, thoracic kyphosis, lumbar lordosis and apical vertebral rotation were compared with the actual measurements immediately in-brace and out-of-brace at the 2-year follow-up. Treatment outcomes according to simulated compliance scenarios of no brace-wear vs. full brace-wear were also computed. The average immediate in-brace difference between the simulated and actual Cobb angle was 4.9° (MT) and 3.7° (TL/L). 2-year out-of-brace, it was 5.6° (MT) and 5.4° (TL/L). The no brace-wear and full brace-wear compliance scenarios resulted respectively in 15/35 (43%) and 31/35 (89%) simulated spine deformities progressing by less than 5° over 2 years of treatment. Clinical significance: the FEM's ability to simulate the final correction with an accuracy on the order of the radiological measurements' interoperator reproducibility, combined with its sensitivity to brace-wear compliance, provides confidence in the model's predictions for a comparative context of use like improving a brace's design before its application. This article is protected by copyright. All rights reserved.

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.28mm, 19.99mm, 27.06mm, 43.99mm, and 50.1mm 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.3mm 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.

The role of muscle forces and gait cycle discretization when assessing acetabular cup primary stability: A finite element study

The aim of this study was to investigate the influence of the muscle force contribution and loading cycle discretization on the predicted micromotion and interfacial bone strains in the implanted acetabulum. To this end, a patient specific finite element model of the hemipelvis was developed, based on the CT-scan and gait analysis results, collected as part of the authors’ previous work. Outcomes of this study suggests that the acetabular cup micromotion and interfacial bone strains can be predicted just using the joint contact force. This helps to reduce the complexity of the finite element models by ignoring the contribution of muscle forces and the associated challenges of mapping these forces to the pelvis. However, the gait cycle needs to be adequately discretised to capture the micromotion at the bone-implant interface. Background and ObjectiveThe Dalstra load case, which includes muscle forces, has been widely adopted in the literature for studying the mechanical environment in the intact and implanted acetabulum. To simplify the modelling approach, some researchers ignore the contribution of muscle forces. The Dalstra load case is also divided into eight separate load steps (five in the stance phase and three in the swing phase), however, it is unclear whether this adequately captures the micromotions, for a cementless acetabular cup, during a simulated activity. The aim of this study was to investigate the influence of the muscle force contribution and loading cycle discretization on the predicted micromotion and interfacial bone strains. MethodsIn this work, a patient specific finite element model of the hemipelvis was developed, based on the CT-scan and gait analysis results, collected as part of the authors’ previous work. Finite element simulations were performed using the joint contact and muscle forces derived from two sources. The first approach was used the load case proposed by Dalstra et al. The second approach used joint contact and muscle forces predicted by a musculoskeletal model. Additionally, the musculoskeletal load case was discretised into 50 equal load steps and the results compared with the equivalent Dalstra load steps. ResultsThe results showed that the contribution of the muscle forces resulted in minor differences in both the magnitude and distribution of the predicted acetabular micromotion (up to 4.01% in the mean acetabular micromotion) and interfacial bone strains (up to 10.34% in the mean interfacial bone strains). The degree of gait cycle discretisation had a significant influence on the acetabular micromotion with a difference of 20.89% in the mean acetabular micromotion. ConclusionOutcomes of this study suggests that the acetabular cup micromotion and interfacial bone strains can be predicted just using the joint contact force. This helps to reduce the complexity of the finite element models by ignoring the contribution of muscle forces and the associated challenges of mapping these forces to the pelvis. However, the gait cycle needs to be adequately discretised to capture the micromotion at the bone-implant interface.

Patient-specific finite element modeling of scoliotic curve progression using region-specific stress-modulated vertebral growth

PurposeThis study describes the creation of patient-specific (PS) osteo-ligamentous finite element (FE) models of the spine, ribcage, and pelvis, simulation of up to three years of region-specific, stress-modulated growth, and validation of simulated curve progression with patient clinical angle measurements. Research Question: Does the inclusion of region-specific, stress-modulated vertebral growth, in addition to scaling based on age, weight, skeletal maturity, and spine flexibility allow for clinically accurate scoliotic curve progression prediction in patient-specific FE models of the spine, ribcage, and pelvis?MethodsFrontal, lateral, and lateral bending X-Rays of five AIS patients were obtained for approximately three-year timespans. PS-FE models were generated by morphing a normative template FE model with landmark points obtained from patient X-rays at the initial X-ray timepoint. Vertebral growth behavior and response to stress, as well as model material properties were made patient-specific based on several prognostic factors. Spine curvature angles from the PS–FE models were compared to the corresponding X-ray measurements.ResultsAverage FE model errors were 6.3 ± 4.6°, 12.2 ± 6.6°, 8.9 ± 7.7°, and 5.3 ± 3.4° for thoracic Cobb, lumbar Cobb, kyphosis, and lordosis angles, respectively. Average error in prediction of vertebral wedging at the apex and adjacent levels was 3.2 ± 2.2°. Vertebral column stress ranged from 0.11 MPa in tension to 0.79 MPa in compression.ConclusionIntegration of region-specific stress-modulated growth, as well as adjustment of growth and material properties based on patient-specific data yielded clinically useful prediction accuracy while maintaining physiological stress magnitudes. This framework can be further developed for PS surgical simulation.

Finite element analysis of screw fixation durability under multiple boundary and loading conditions for a custom pelvic implant.

Despite showing promising functional outcomes for pelvic reconstruction after sarcoma resection, custom-made pelvic implants continue to exhibit high complication rates due to fixation failures. Patient-specific finite element models have been utilized by researchers to evaluate implant durability. However, the effect of assumed boundary and loading conditions on failure analysis results of fixation screws remains unknown. In this study, the postoperative stress distributions in the fixation screws of a state-of-the-art custom-made pelvic implant were simulated, and the risk of failure was estimated under various combinations of two bone-implant interaction models (tied vs. frictional contact) and four load cases from level-ground walking and stair activities. The study found that the average weighted peak von Mises stress could increase by 22-fold when the bone-implant interactions were modeled with a frictional contact model instead of a tied model, and the likelihood of fatigue and pullout failure for each screw could change dramatically when different combinations of boundary and loading conditions were used. The inclusion of additional boundary and loading conditions led to a more reliable analysis of fixation durability. These findings demonstrated the importance of simulating multiple boundary conditions and load cases for comprehensive implant design evaluation using finite element analysis.

Spinal Cord Stress After Anterior Cervical Diskectomy and Fusion: Results from a Patient-Specific Finite Element Model.

Degenerative cervical myelopathy (DCM) is the commonest cause of cervical spinal cord dysfunction in older adults and is characterized by spinal cord compression and stress during neck motion. Although surgical decompression eliminates static spinal cord compression, cord stress resulting from flexion-extension motion of the spinal column has not been determined for single and multi-level surgical interventions. The effect of surgery on spinal cord stress is expected to change with the number of surgical levels as well as patient-specific anatomy. Using a MRI-derived patient-specific finite element model, we simulated 1-, 2- and 3-level anterior cervical diskectomy and fusion (ACDF) surgery for DCM. A substantial decrease in spinal cord stress at the level of spinal cord decompression was noted in all simulations. This was associated with a considerable increase in spinal cord stress rostral to the surgical level, and the magnitude of stress was higher in multi-level surgery. Increased spinal cord stress at the rostral adjacent segment correlated with increased segmental range of motion (r = 0.69, p = 0.002) and disk pressure (r = 0.57, p = 0.05). Together, these results indicate that ACDF for DCM is associated with adversespinal cord stress patterns adjacent to the fusion construct, and further research is needed to determine if the altered stress is associated with clinical outcomes after surgery for DCM.

Rethinking the 10% strain rule in fracture healing: A distal femur fracture case series.

Since the 1970s, the 2%-10% rule has been used to describe the range of interfragmentary gap closure strains that are conducive for secondary bone healing. Interpreting the available evidence for the association between strain and bone healing remains challenging because interfragmentary strain is impossible to directly measure in vivo. The question of how much strain occurs within and around the fracture gap is also difficult to resolve using bench tests with osteotomy models because these do not reflect the complexity of injury patterns seen in the clinic. To account for these challenges, we used finite element modeling to assess the three-dimensional interfragmentary strain in a case series of naturally occurring distal femur fractures treated with lateral plating under load conditions representative of the early postoperative period. Preoperative computed tomography scans were used to construct patient-specific finite element models and plate fixation constructs to match the operative management of each patient. The simulations showed that gap strains were within 2%-10% only for the lowest load application level, 20% static body weight (BW). Moderate loading of 60% static BW and above caused gap strains that far exceeded 10%, but in all cases, strains in the periosteal region external to the fracture line remained low. Comparing these findings with postoperative radiographs suggests that in vivo secondary healing of distal femur fractures may be robust to early gap strains much greater than 10% because formation of new bone is initiated outside the gap where strains are lower, followed by later consolidation within the gap.

Experimental validation of a subject-specific finite element model of lumbar spine segment using digital image correlation.

Pathologies such as cancer metastasis and osteoporosis strongly affect the mechanical properties of the vertebral bone and increase the risk of fragility fractures. The prediction of the fracture risk with a patient-specific model, directly generated from the diagnostic images of the patient, could help the clinician in the choice of the correct therapy to follow. But before such models can be used to support any clinical decision, their credibility must be demonstrated through verification, validation, and uncertainty quantification. In this study we describe a procedure for the generation of such patient-specific finite element models and present a first validation of the kinematics of the spine segment. Quantitative computed tomography images of a cadaveric lumbar spine segment presenting vertebral metastatic lesions were used to generate the model. The applied boundary conditions replicated a specific experimental test where the spine segment was loaded in compression-flexion. Model predictions in terms of vertebral surface displacements were compared against the full-field experimental displacements measured with Digital Image Correlation. A good agreement was obtained from the local comparison between experimental data and simulation results (R2 > 0.9 and RMSE% <8%). In conclusion, this work demonstrates the possibility to apply the developed modelling pipeline to predict the displacement field of human spine segment under physiological loading conditions, which is a first fundamental step in the credibility assessment of these clinical decision-support technology.

The relationship between subchondral bone cysts and cartilage health in the Tibiotalar joint: A finite element analysis

BackgroundSubchondral bone cysts are a common presentation in ankle haemarthropathy. The relationship with ankle joint health has however not previously been investigated. The aim of this study was to assess the influence of subchondral bone cysts of differing shapes, volumes and depths on joint health. MethodsChronologically sequential Magnetic Resonance imaging scans of four hemophilic ankles with subchondral bone cysts present (N = 18) were used to build patient specific finite element models under two cystic conditions to assess their influence on cartilage contact pressures. Variables such as location, volume and depth were considered individually, to investigate whether certain cystic conditions may be more detrimental to cartilage health. FindingsSignificant quantifiable contact redistribution was seen in the presence of subchondral bone cysts and this redistribution reflected the shape and size of the cysts, however, with the presence of cysts in both bones in 10 of the 18 cases a direct relationship to volume could not be correlated. InterpretationThis work demonstrated a redistribution of contact pressures in the presence of subchondral bone cysts. This alteration to loading history could be linked to cartilage degeneration due to the biological response to abnormal loading.