Patient-specific Finite Element Models Research Articles

Immediate Correction of Idiopathic Scoliosis with Nighttime Braces Created by A fully-automated Generative Design Algorithm: A Randomized Controlled Crossover Trial.

Single center, double blinded, prospective crossover randomized controlled trial. To clinically validate the efficacy of nighttime braces designed automatically by a generative design algorithm to treat idiopathic scoliosis (IS). The tested hypothesis was the clinical equivalence of immediate in-brace correction for the new automatically generated brace design versus a standard Providence-type brace. Documented efficacy of brace treatment varies between centers, and depends on the empirical expertise of the treating orthotist. Our group previously developed a fully-automated generative brace design algorithm that leverages a patient-specific finite element model (FEM) to optimize brace geometry and correction before its fabrication. Fifty-eight skeletally immature patients diagnosed with IS, aged between 10 to 16 years were recruited. All patients received both a nighttime brace automatically generated by the algorithm (Test) and a Providence-type brace designed by an expert orthotist (Control). Radiographs were taken for each patient with both braces in a randomized crossover approach to evaluate immediate in-brace correction. The targeted fifty-five patients (48 females, 7 males) completed the study. The immediate Cobb angle correction was 57% ± 19 (Test) versus 58% ± 21 (Control) for the main thoracic (MT) curve, while it was 89% ± 25 (Test) versus 87% ± 28 (Control) for the thoraco-lumbar/lumbar (TLL) spine. The immediate correction with the Test brace was non-inferior to that of the Control brace (P < 0.001). The order in which the braces were tested did not have a residual effect on the immediate correction. The fully-automated generative brace design algorithm proves to be clinically relevant, allowing for immediate in-brace correction equivalent to that of braces designed by expert orthotists. Patient 2 years follow-up will continue. This method's integration could help design and rationalize the design of braces for the treatment of IS. Level 2.

From Simulations to Inference: Using Machine Learning to Tune Patient-Specific Finite-Element Models of the Middle Ear Towards Objective Diagnosis.

Computational models, particularly finite-element (FE) models, are essential for interpreting experimental data and predicting system behavior, especially when direct measurements are limited. A major challenge in tuning these models is the large number of parameters involved. Traditional methods, such as one-by-one sensitivity analyses, are time-consuming, subjective, and often return only a single set of parameter values, focusing on reproducing averaged data rather than capturing the full variability of experimental measurements. In this study, we applied simulation-based inference (SBI) using neural posterior estimation (NPE) to tune an FE model of the human middle ear. The training dataset consisted of 10,000 FE simulations of stapes velocity, ear-canal (EC) input impedance, and absorbance, paired with seven FE parameter values randomly sampled within plausible ranges. The neural network learned the association between parameters and simulation outcomes, returning the probability distribution of parameter values that can reproduce experimental data. Our approach successfully identified parameter sets that reproduced three experimental datasets simultaneously. By accounting for experimental noise and variability during training, the method provided a probability distribution of parameters, representing all valid combinations that could fit the data, rather than tuning to averaged values. The network demonstrated robustness to noise and exhibited an efficient learning curve due to the large training dataset. SBI offers an objective alternative to laborious sensitivity analyses, providing probability distributions for each parameter and uncovering interactions between them. This method can be applied to any biological FE model, and we demonstrated its effectiveness using a middle-ear model. Importantly, it holds promise for objective differential diagnosis of conductive hearing loss by providing insight into the mechanical properties of the middle ear.

A Computer Modeling-Based Target Zone for Transposition Osteotomy of the Acetabulum in Patients with Hip Dysplasia.

This study aimed to determine the acetabular position to optimize hip biomechanics after transposition osteotomy of the acetabulum (TOA), a specific form of periacetabular osteotomy, in patients with hip dysplasia. We created patient-specific finite-element models of 46 patients with hip dysplasia to simulate 12 virtual TOA scenarios: lateral rotation to achieve a lateral center-edge angle (LCEA) of 30°, 35°, and 40° combined with anterior rotation of 0°, 5°, 10°, and 15°. Joint contact pressure (CP) on the acetabular cartilage during a single-leg stance and simulated hip range of motion without osseous impingement were calculated. The optimal acetabular position was defined as satisfying both normal joint CP and the required range of motion for activities of daily living. Multivariable logistic regression analysis was used to identify preoperative morphological predictors of osseous impingement after virtual TOA with adequate acetabular correction. The prevalence of hips in the optimal position was highest (65.2%) at an LCEA of 30°, regardless of the amount of anterior rotation. While the acetabular position minimizing peak CP varied among patients, approximately 80% exhibited normalized peak CP at an LCEA of 30° and 35° with 15° of anterior rotation, which were the 2 most favorable configurations among the 12 simulated scenarios. In this context, the preoperative head-neck offset ratio (HNOR) at the 1:30 clock position (p = 0.018) was an independent predictor of postoperative osseous impingement within the required range of motion. Specifically, an HNOR of <0.14 at the 1:30 clock position predicted limitation of required range of motion after virtual TOA (sensitivity, 57%; specificity, 81%; and area under the receiver operating characteristic curve, 0.70). Acetabular reorientation to an LCEA of between 30° and 35° with an additional 15° of anterior rotation may serve as a biomechanics-based target zone for surgeons performing TOA in most patients with hip dysplasia. However, patients with a reduced HNOR at the 1:30 clock position may experience limited range of motion in activities of daily living postoperatively. This study provides a biomechanics-based target for refining acetabular reorientation strategies during TOA while considering morphological factors that may limit the required range of motion.

A Soft-Tissue Driven Bone Remodeling Algorithm for Mandibular Residual Ridge Resorption Based on Patient CT Image Data.

The role of the biomechanical stimulation generated from soft tissue has not been well quantified or separated from the self-regulated hard tissue remodeling governed by Wolff's Law. Prosthodontic overdentures, commonly used to restore masticatory functions, can cause localized ischemia and inflammation as they often compress patients' oral mucosa and impede local circulation. This biomechanical stimulus in mucosa is found to accelerate the self-regulated residual ridge resorption (RRR), posing ongoing clinical challenges. Based on the dedicated long-term clinical datasets, this work develops an in-silico framework with a combination of techniques, including advanced image post-processing, patient-specific finite element models and unsupervised machine learning Self-Organizing map algorithm, to identify the soft tissue induced RRR and quantitatively elucidate the governing relationship between the RRR and hydrostatic pressure in mucosa. The proposed governing equation has not only enabled a predictive simulation for RRR as showcased in this study, providing a biomechanical basis for optimizing prosthodontic treatments, but also extended the understanding of the mechanobiological responses in the soft-hard tissue interfaces and the role in bone remodeling.

Incorporating pathological gait into patient-specific finite element models of the haemophilic ankle

Haemarthrosis is an inherent clinical feature of haemophilia, a disease characterised by an absence or reduction in clotting proteins. Patients with severe haemophilia experience joint bleeding leading to blood-induced ankle arthropathy (haemarthropathy). Altered biomechanics of the ankle have been reported in people with haemophilia; however, the consequence of this on joint health is little understood. The aim of this study was to assess the changes in joint contact due to haemophilia disease-specific gait features using patient-specific modelling, to better understand the link between biomechanics and joint outcomes. Four, image-based, finite element models of haemophilic ankles were simulated through consecutive events in the stance phase of gait, using both patient-specific and healthy control group (n = 36) biomechanical inputs. One healthy control FE model was simulated through the healthy control stance phase of the gait cycle for a point of comparison. The method developed allowed cartilage contact mechanics to be assessed throughout the loading phase of the gait cycle. This showed areas of increased contact pressure in the medial and lateral regions of the talar dome, which may be linked to collapse in these regions. This method may allow the relationship between structure and function in the tibiotalar joint to be better understood.

Finite element modeling with patient-specific geometry to assess clinical risks of percutaneous pulmonary valve implantation.

Percutaneous pulmonary valve implantation (PPVI) is a non-surgical treatment for right ventricular outflow tract (RVOT) dysfunction. During PPVI, a stented valve, delivered via catheter, replaces the dysfunctional pulmonary valve. Stent oversizing allows valve anchoring within the RVOT, but overexpansion can intrude on the surrounding structures. Potentially dangerous outcomes include aortic valve insufficiency (AVI) from aortic root (AR) distortion and myocardial ischemia from coronary artery (CA) compression. Currently, risks are evaluated via balloon angioplasty/sizing before stent deployment. Patient-specific finite element (FE) analysis frameworks can improve pre-procedural risk assessment, but current methods require hundreds of hours of high-performance computation. We created a simplified method to simulate the procedure using patient-specific FE models for accurate, efficient pre-procedural PPVI (using balloon expandable valves) risk assessment. The methodology was tested by retrospectively evaluating the clinical outcome of 12 PPVI candidates. Of 12 patients (median age 14.5 years) with dysfunctional RVOT, 7 had native RVOT and 5 had RV-PA conduits. Seven patients had undergone successful RVOT stent/valve placement, three had significant AVI on balloon testing, one had left CA compression, and one had both AVI and left CA compression. A model-calculated change of more than 20% in lumen diameter of the AR or coronary arteries correctly predicted aortic valve sufficiency and/or CA compression in all the patients. Agreement between FE results and clinical outcomes is excellent. Additionally, these models run in 2-6 min on a desktop computer, demonstrating potential use of FE analysis for pre-procedural risk assessment of PPVI in a clinically relevant timeframe.

Patient-specific finite element analysis for assessing hip fracture risk in aging populations

The femur is one of the most important bone in the human body, as it supports the body’s weight and helps with movement. The aging global population presents a significant challenge, leading to an increasing demand for artificial joints, particularly in knee and hip replacements, which are among the most prevalent surgical procedures worldwide. This study focuses on hip fractures, a common consequence of osteoporotic fractures in the elderly population. To accurately predict individual bone properties and assess fracture risk, patient-specific finite element models (FEM) were developed using CT data from healthy male individuals. The study employed ANSYS 2023 R2 software to estimate fracture loads under simulated single stance loading conditions, considering strain-based failure criteria. The FEM bone models underwent meticulous reconstruction, incorporating geometrical and mechanical properties crucial for fracture risk assessment. Results revealed an underestimation of the ultimate bearing capacity of bones, indicating potential fractures even during routine activities. The study explored variations in bone density, failure loads, and density/load ratios among different specimens, emphasizing the complexity of bone strength determination. Discussion of findings highlighted discrepancies between simulation results and previous studies, suggesting the need for optimization in modelling approaches. The strain-based yield criterion proved accurate in predicting fracture initiation but required adjustments for better load predictions. The study underscores the importance of refining density-elasticity relationships, investigating boundary conditions, and optimizing models through in vitro testing for enhanced clinical applicability in assessing hip fracture risk. In conclusion, this research contributes valuable insights into developing patient-specific FEM bone models for clinical hip fracture risk assessment, emphasizing the need for further refinement and optimization for accurate predictions and enhanced clinical utility.

Optimization of S2-alar-iliac screw (S2AI) fixation in adult spine deformity using a comprehensive genetic algorithm and finite element model personalized to patient geometry and bone mechanical properties.

To optimize the biomechanical performance of S2AI screw fixation using a genetic algorithm (GA) and patient-specific finite element analysis integrating bone mechanical properties. Patient-specific pelvic finite element models (FEM), including one normal and one osteoporotic model, were created from bi-planar multi-energy X-rays (BMEXs). The genetic algorithm (GA) optimized screw parameters based on bone mass quality (BM method) while a comparative optimization method maximized the screw corridor radius (GEO method). Biomechanical performance was evaluated through simulations, comparing both methods using pullout and toggle tests. The optimal screw trajectory using the BM method was more lateral and caudal with insertion angles ranging from 49° to 66° (sagittal plane) and 29° to 35° (transverse plane). In comparison, the GEO method had ranges of 44° to 54° and 24° to 30° respectively. Pullout forces (PF) using the BM method ranged from 5 to 18.4kN, which were 2.4 times higher than the GEO method (2.1-7.7kN). Toggle loading generated failure forces between 0.8 and 10.1kN (BM method) and 0.9-2.9kN (GEO method). The bone mass surrounding the screw representing the fitness score and PF of the osteoporotic case were correlated (R2 > 0.8). Our study proposed a patient-specific FEM to optimize the S2AI screw size and trajectory using a robust BM approach with GA. This approach considers surgical constraints and consistently improves fixation performance.

Influence of build orientation and support structure on additive manufacturing of human knee replacements: a computational study.

Developing patient-specific implants has an increasing interest in the application of emerging additive manufacturing (AM) technologies. On the other hand, despite advances in total knee replacement (TKR), studies suggest that up to 20% of patients with elective TKR are dissatisfied with the outcome. By creating 3D objects from digital models, AM enables the production of patient-specific implants with complex geometries, such as those required for knee replacements. Previous studies have highlighted concerns regarding the risk of residual stresses and shape distortions in AM parts, which could lead to structural failure or other complications. This article presents a computational framework that uses CT images to create patient-specific finite element models for optimizing AM knee replacements. The workflow includes image processing in the open-source software 3DSlicer and MeshLab and AM process simulations in the commercial platform 3DEXPERIENCE. The approach is demonstrated on a distal femur replacement for a 50-year-old male patient from the open-access Natural Knee Data. The results show that build orientations have a significant impact on both shape distortions and residual stresses. Support structures have a marginal effect on residual stresses but strongly influence shape distortions, whereas conical support exhibits a maximum distortion of 18.5 mm. Future research can explore how these factors affect the functionality of AM knee replacements under in-service loading.

The effect of trapeziometacarpal joint passive stiffness on mechanical loadings of cartilages

Hypermobility of the trapeziometacarpal joint is commonly considered to be a potential risk factor for osteoarthritis. Nevertheless, the results remain controversial due to a lack of quantitative validation. The objective of this study was to evaluate the effect of joint laxity on the mechanical loadings of cartilage.A patient-specific finite element model of trapeziometacarpal joint passive stiffness was developed. The joint passive stiffness was modeled by creating linear springs all around the joint. The linear spring stiffness was determined by using an optimization process to fit force–displacement data measured during laxity tests performed on eight healthy volunteers. The estimated passive stiffness parameters were then included in a full thumb finite element simulation of a pinch grip task driven by muscle forces to evaluate the effect on trapeziometacarpal loading. The correlation between stiffness and the loading of cartilage in terms of joint contact pressure and maximum shear strain was analyzed.A significant negative correlation was found between the trapeziometacarpal joint passive stiffness and the contact pressure on trapezium cartilage during the simulated pinch grip task.These results therefore suggest that the hypermobility of the trapeziometacarpal joint could affect the contact pressure on trapezium cartilage and support the existence of an increased risk associated with hypermobility.

Development of patient-specific finite element model for study of composite dental implants

Traumatic dental injuries can occur due to various reasons such as accidents, sports injuries, fights, falls, and others. These injuries can affect the teeth, gums, and surrounding tissues, and can range from minor chips and cracks to severe fractures, dislocations, and avulsions (when the tooth is completely knocked out of the socket). The most common way to address this is by replacing affected teeth with dental implants. The purpose of this research is to evaluate the use of composite materials in dental implants and compare them with the traditionally used materials using a patient specific cone beam computed tomography (CBCT) based finite element model (FEM). To conduct this research, two different implant groups i.e., traditional implant and composite implant were designed using Titanium grade 4, zirconium oxide-reinforced lithium silicate (ZLS), and Zirconia (ZrO2). Six dental implants were designed namely Ti implant, ZLS implant, ZrO2 implant, Ti-ZrO2 composite, Ti-ZLS composite, and ZLS-ZrO2 composite using 3D modelling software. Detailed full-scale 3D models of patient specific dental implant were developed and traumatic loading conditions were applied to the enamel of central incisor teeth or crown of dental implant, and maxilla was constrained in all directions. It was found that the use of composite materials for dental implants can reduce the stresses over the surface of abutment and implant as compared to traditional implants. The detailed models developed as a part of this study can advance the research on dental implants, and with further experimental validation allow the use of composite materials for fabrication of more stable dental implants.

Automated design of nighttime braces for adolescent idiopathic scoliosis with global shape optimization using a patient-specific finite element model

Adolescent idiopathic scoliosis is a complex three-dimensional deformity of the spine, the moderate forms of which require treatment with an orthopedic brace. Existing brace design approaches rely mainly on empirical manual processes, vary considerably depending on the training and expertise of the orthotist, and do not always guarantee biomechanical effectiveness. To address these issues, we propose a new automated design method for creating bespoke nighttime braces requiring virtually no user input in the process. From standard biplanar radiographs and a surface topography torso scan, a personalized finite element model of the patient is created to simulate bracing and the resulting spine growth over the treatment period. Then, the topography of an automatically generated brace is modified and simulated over hundreds of iterations by a clinically driven optimization algorithm aiming to improve brace immediate and long-term effectiveness while respecting safety thresholds. This method was clinically tested on 17 patients prospectively recruited. The optimized braces showed a highly effective immediate correction of the thoracic and lumbar curves (70% and 90% respectively), with no modifications needed to fit the braces onto the patients. In addition, the simulated lumbar lordosis and thoracic apical rotation were improved by 5° ± 3° and 2° ± 3° respectively. Our approach distinguishes from traditional brace design as it relies solely on biomechanically validated models of the patient’s digital twin and a design strategy that is entirely abstracted from empirical knowledge. It provides clinicians with an efficient way to create effective braces without relying on lengthy manual processes and variable orthotist expertise to ensure a proper correction of scoliosis.

Biomechanical Effect of Distal Tibial Oblique Osteotomy: A Preliminary Finite-Element Analysis.

The biomechanical effect of distal tibial oblique osteotomy (DTOO) on osteoarthritic ankles has not been investigated. Using finite element (FE) models, we aimed to elucidate the effect of DTOO on the ankle contact pressure (CP) distribution. This study included two patients with ankle osteoarthritis who underwent DTOO and one asymptomatic control. Patient-specific FE models were reconstructed by matching standing radiographs with supine computed tomography scans. The joint contact area (CA) and maximumCPon the articular surface of the talus were calculated before and after DTOO and compared with those of the control. In the control, the CA was 584 mm2 and the maximum CP was 2.6 MPa. In case 1, the CA increased by 125% from 166 mm2 preoperatively to 375 mm2 postoperatively, accompanied by a 36% decrease in the maximum CP from 9.8 MPa to 6.3 MPa. Similarly, in case 2, the CA increased by 46% from 301 mm2 to 439 mm2, accompanied by a 27% decrease in the maximum CP from 6.7 MPa to 4.9 MPa. This study suggests DTOO improves the biomechanics of the ankle, but not sufficiently compared to the control. This analytical approach may enhance understanding of ankle pathophysiology and assist in the design of the ideal corrective osteotomy.

Optimized braces for the treatment of adolescent idiopathic scoliosis: A study protocol of a prospective randomised controlled trial.

Adolescent Idiopathic Scoliosis (AIS) is a 3D deformity of the spine that affects 3% of the adolescent population. Conservative treatments like bracing aim to halt the progression of the curve to the surgical threshold. Computer-aided design and manufacturing (CAD/CAM) methods for brace design and manufacturing are becoming increasingly used. Linked to CAD/CAM and 3D radiographic reconstruction techniques, we developed a finite element model (FEM) enabling to simulate the brace effectiveness before its fabrication, as well as a semi-automatic design processes. The objective of this randomized controlled trial is to compare and validate such FEM semi-automatic algorithm used to design nighttime Providence-type braces. Fifty-eight patients with AIS aged between 10 to 16-years and skeletally immature will be recruited. At the delivery stage, all patients will receive both a Providence-type brace optimized by the semi-automatic algorithm leveraging a patient-specific FEM (Test) and a conventional Providence-type brace (Control), both designed using CAD/CAM methods. Biplanar radiographs will be taken for each patient with both braces in a randomized crossover approach to evaluate immediate correction. Patients will then be randomized to keep either the Test or Control brace as prescribed with a renewal if necessary, and will be followed over two years. The primary outcome will be the change in Cobb angle of the main curve after two years. Secondary outcomes will be brace failure rate, quality of life (QoL) and immediate in-brace correction. This is a single-centre study, double-blinded (participant and outcome assessor) randomized controlled trial (RCT). ClinicalTrials.gov: NCT05001568.

Exploring uncertainty in hyper-viscoelastic properties of scalp skin through patient-specific finite element models for reconstructive surgery

Understanding skin responses to external forces is crucial for post-cutaneous flap wound healing. However, the in vivo viscoelastic behavior of scalp skin remains poorly understood. Personalized virtual surgery simulations offer a way to study tissue responses in relevant 3D geometries. Yet, anticipating wound risk remains challenging due to limited data on skin viscoelasticity, which hinders our ability to determine the interplay between wound size and stress levels. To bridge this gap, we reexamine three clinical cases involving scalp reconstruction using patient-specific geometric models and employ uncertainty quantification through a Monte Carlo simulation approach to study the effect of skin viscoelasticity on the final stress levels from reconstructive surgery. Utilizing the generalized Maxwell model via the Prony series, we can parameterize and efficiently sample a realistic range of viscoelastic response and thus shed light on the influence of viscoelastic material uncertainty in surgical scenarios. Our analysis identifies regions at risk of wound complications based on reported threshold stress values from the literature and highlights the significance of focusing on long-term responses rather than short-term ones.

Biomechanical effects of Evans versus Hintermann osteotomy for treating adult acquired flatfoot deformity: a patient-specific finite element investigation.

Evans and Hintermann lateral column lengthening (LCL) procedures are both widely used to correct adult acquired flatfoot deformity (AAFD), and have both shown good clinical results. The aim of this study was to compare these two procedures in terms of corrective ability and biomechanics influence on the Chopart and subtalar joints through finite element (FE) analysis. Twelve patient-specific FE models were established and validated. The Hintermann osteotomy was performed between the medial and posterior facets of the subtalar joint; while, the Evans osteotomy was performed on the anterior neck of the calcaneus around 10mm from the calcaneocuboid joint surface. In each procedure, a triangular wedge of varying size was inserted at the lateral edge. The two procedures were then compared based on the measured strains of superomedial calcaneonavicular ligaments and planter facia, the talus-first metatarsal angle, and the contact characteristics of talonavicular, calcaneocuboid and subtalar joints. The Hintermann procedure achieved a greater correction of the talus-first metatarsal angle than Evans when using grafts of the same size, indicating that Hintermann had stronger corrective ability. However, its distributions of von-Mises stress in the subtalar, talonavicular and calcaneocuboid joints were lesshomogeneous than those of Evans. In addition, the strains of superomedial calcaneonavicular ligaments and planter facia of Hintermann were also greater than those of Evans, but both generally within the safe range (less than 6%). This FE analysis study indicates that both Evans and Hintermann procedures have good corrective ability for AAFD. Compared to Evans, Hintermann procedure can provide a stronger corrective effect while causing greater disturbance to the biomechanics of Chopart joints, which may be an important mechanismof arthritis. Nevertheless, it yields a better protection to the subtalar joint than Evans osteotomy. Both Evans and Hintermann LCL surgeries have a considerable impact on adjacent joints and ligament tissues. Such effects alongside the overcorrection problem should be cautiously considered when choosing the specific surgical method. Level III, case-control study.

Patient-specific finite element modeling for fracture risk prediction in a canine model of osteosarcoma.

In cancer patients with bone tumors, pathological fractures are a major concern. Making treatment decision for these patients requires an evaluation of fracture risk, which is currently based on semi-qualitative criteria that lack patient-specificity. Because of this, there exists a need for quantitative fracture risk prediction tailored to the patient's individual bone geometry. To address this need, this study aims to develop and validate a finite element (FE) technique that can be used to create patient-specific models and more accurately identify fracture risk. Model validation was performed using canine radii. Radii were harvested from eight canines euthanized for reasons unrelated to the study. A semicircular osteotomy was made in the distal portion of each bone to simulate tumor lysis. Samples underwent computed tomography (CT) scanning and were randomly assigned to loading groups for destructive mechanical testing. Three samples were tested in torsion, three in cantilever bending, and two in compression. FE models were created for each bone from the corresponding CT scan to replicate patient-specific geometry. Material properties were based on equations relating scan properties to elastic modulus. Boundary conditions and loads were added to the models based on the sample's treatment group. Stiffness and strain data were collected from both the mechanical testing and FE simulation, and yield load predictions were made based on maximum principal strain. Experimental and computational results were compared using a linear regression. The FE models were most accurate in predicting stiffness, followed by strain, with yield load having the lowest accuracy. Linear regressions resulted in R2 values of 0.9335 for bending and compression and 0.8798 for torsion. The proposed FE technique is a valid method for predicting fracture in a canine model of osteosarcoma. This method could provide patient-specific, quantitative data to aid clinicians in decisions regarding surgical intervention for patients with bone tumors.

Study of the Influence of Boundary Conditions on Corneal Deformation Based on the Finite Element Method of a Corneal Biomechanics Model.

Implementing in silico corneal biomechanical models for surgery applications can be boosted by developing patient-specific finite element models adapted to clinical requirements and optimized to reduce computational times. This research proposes a novel corneal multizone-based finite element model with octants and circumferential zones of clinical interest for material definition. The proposed model was applied to four patient-specific physiological geometries of keratoconus-affected corneas. Free-stress geometries were calculated by two iterative methods, the displacements and prestress methods, and the influence of two boundary conditions: embedded and pivoting. The results showed that the displacements, stress and strain fields differed for the stress-free geometry but were similar and strongly depended on the boundary conditions for the estimated physiological geometry when considering both iterative methods. The comparison between the embedded and pivoting boundary conditions showed bigger differences in the posterior limbus zone, which remained closer in the central zone. The computational calculation times for the stress-free geometries were evaluated. The results revealed that the computational time was prolonged with disease severity, and the displacements method was faster in all the analyzed cases. Computational times can be reduced with multicore parallel calculation, which offers the possibility of applying patient-specific finite element models in clinical applications.

Single-plane osteotomy model is inaccurate for evaluating the optimal strategy in treating vertical femoral neck fractures: A finite element analysis

Background and objectivesThe conventional method for simulating vertical femoral neck fractures (vFNFs) is via a vertical single-plane osteotomy (SPO) across the entire femur. However, the accuracy of SPO for evaluating the optimal internal fixation strategy (IFS) and the appropriate assessment parameters is not clear. This study thus aimed to examine the accuracy of SPO in evaluating IFSs and to identify appropriate evaluation parameters using finite element analysis. MethodsEighty patient-specific finite element models were developed based on CT images from eight vFNF patients. The natural fracture model was built using structural features of the affected side, while the SPO was simulated on the healthy side. Five different IFSs were applied to both the natural fracture and SPO groups. Thirteen parameters, including stress, displacement, and stiffness, were subjected to a two-way repeated measures ANOVA to determine the effect of IFSs and fracture morphology on stability. A Pearson correlation analysis was performed on varied parameters with various IFSs to identify independent parameters. Based on these independent parameters, the entropy evaluation method (EEM) score was used to rank the performance of IFSs for each patient. ResultsEight of the thirteen parameters were significantly influenced by IFSs (p < 0.05), two by fracture morphology (p < 0.01), and none by the interaction between IFS and fracture morphology. In the natural fracture group, parameters including screw stress and displacement, bone cut rate (BCR), and compression effects varied independently with distinct IFSs. In the SPO group, trunk displacement, BCR, cut-out risk, and compression effects parameters changed independently. The BCR of the Alpha strategy was significantly higher than that of the Inverted strategy in the natural fracture group (p = 0.002), whereas the opposite was observed in the SPO group (p = 0.016). Regarding compression effects, two IFS pairings in the natural fracture group and seven IFS pairings in the SPO group exhibited significant differences. None of the five IFSs achieved the optimal EEM score for each patient. ConclusionsThe single-plane osteotomy model may have limitations in assessing IFSs, particularly when the bone cut rate and compression effects are the main influencing factors. Parameters of the screw stress and displacement, BCR, and compression effects appear to be relevant in evaluating IFSs for natural fracture models. It indicates that individualized natural fracture models could provide more comprehensive insights for determining the optimal IFS in treating vFNFs.

MINING OF BIOMECHANICAL AND GEOMETRY DATA OF INTERVERTEBRAL DISC FINITE ELEMENT SIMULATIONS

This study investigates the relationships between Intervertebral Disc (IVD) morphology and biomechanics using patient-specific (PS) finite element (FE) models and poromechanical simulations.169 3D lumbar IVD shapes from the European project MySpine (FP7-269909), spanning healthy to Pfirrmann grade 4 degeneration, were obtained from MRIs. A Bayesian Coherent Point Drift algorithm aligned meshes to a previously validated structural FE mesh of the IVD. After mesh quality analyses and Hausdorff distance measurements, mechanical simulations were performed: 8 and 16 hours of sleep and daytime, respectively, applying 0.11 and 0.54 MPa of pressure on the upper cartilage endplate (CEP). Simulation results were extracted from the anterior (ATZ) and posterior regions (PTZ) and the center of the nucleus pulposus (CNP). Data mining was performed using Linear Regression, Support Vector Machine, and eXtreme Gradient Boosting techniques. Mechanical variables of interest in DD, such as pore fluid velocity (FLVEL), water content, and swelling pressure, were examined. The morphological variables of the simulated discs were used as input features.Local morphological variables significantly impacted the local mechanical response. The local disc heights, respectively in the mid (mh), anterior (ah), and posterior (ph) regions, were key factors in general. Additionally, fluid transport, reflected by FLVEL, was greatly influenced (r2 0.69) by the shape of the upper and lower cartilage endplates (CEPs).This study suggests that disc morphology affects Mechanical variables of interest in DD. Attention should be paid to the antero-posterior distribution and local effects of disc heights. Surprisingly, the CEP morphology remotely affected the fluid transport in NP volumes around mid-height, and mechanobiological implications shall be explored. In conclusion, patient-specific IVD modeling has strong potential to unravel important correlations between IVD phenotypes and local tissue regulation.Acknowledgments: European Commission: Disc4All-MSCA-2020-ITN-ETN GA: 955735; O-Health-ERC-CoG-2021-101044828