Patient-specific Finite Element Research Articles

Patient-specific finite element analysis of four different fixation methods for transversely unstable radial head fractures

Plate fixation is a common treatment option for radial head fractures (RHFs). Due to the benefits of less invasiveness and fewer complications of internal fixation, the application of small-diameter headless compression screws (HCSs) to treat RHFs has become a new trend. This study aimed to compare the mechanical stability of four distinct internal fixation protocols for transversely unstable RHFs via finite element analysis. Using computed tomography data from 10 patients, we developed 40 patient-specific FE models of transversely unstable RHFs fixed by parallel, crossed, and tripod HCSs and mini-T plate (MTP). Under simulated physiological loading of the elbow joint, the construct stiffness, displacement, and von Mises stresses were evaluated and verified by a biomechanical experiment. Under shear loading, the MTP group exhibited lower construct stiffness, larger displacement, and higher Von Mises stress than the HCSs group. The stiffness of tripod HCSs was greater than parallel and crossed screw fixation techniques. There was a strong relationship between apparent bone density and construct stiffness (R = 0.98 to 0.99). In the treatment of transversely unstable RHFs, HCSs have superior biomechanical stability than MTP. The tripod technique was also more stable than parallel and crossed fixation.

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.

A biomechanical investigation of miniscrews placed on a mandible for temporary anchorage; a patient-specific finite element analysis

Miniscrews are temporary skeletal anchorage devices that are widely used in orthodontic treatment, and their success depends on the placement area, angle, technique, and screw dimensions. This study aimed to investigate the effects of miniscrew lengths, insertion angles, and force directions on a mandible model consisting of teeth, cortical and cancellous bones. One Dental Volumetric Tomography (DVT) scan from a patient who had miniscrews were used for mandibular bone modeling to perform finite element analysis. The model variables included miniscrew lengths (6, 8, and 10 mm), insertion angles (−15°, 45°, 60°, and 90°), and force directions (30°, 45°, and 60°). The minimum and maximum stresses were calculated as 18.61 and 37.11 MPa at 6 mm and 10 mm, respectively. According to the insertion angles, the lowest stress was observed at 60°, while the highest stress was found at 15° in the ventral direction. At force directions, the lowest stress was at 60°, and the highest stress was at 45°. However, there were no significant differences in insertion angles and force directions. A statistically significant difference was determined in miniscrew length. As a result, the best result was calculated to be 6 mm inserted at a 60° angle, which could induce the lowest stress. Increasing the miniscrew length will increase the stress on the mandible. In addition, because of the higher force direction, stress decreases with shorter power arms.

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.

Validation of a patient-specific finite element analysis framework for identification of growing rod-failure regions in early onset scoliosis patients

PurposeGrowing rods are the gold-standard for treatment of early onset scoliosis (EOS). However, these implanted rods experience frequent fractures, requiring additional surgery. A recent study by the U.S. Food and Drug Administration (FDA) identified four common rod fracture locations. Leveraging this data, Agarwal et al. were able to correlate these fractures to high-stress regions using a novel finite element analysis (FEA) framework for one patient. The current study aims to further validate this framework through FEA modeling extended to multiple patients.MethodsThree patient-specific FEA models were developed to match the pre-operative patient data taken from both registry and biplanar radiographs. The surgical procedure was then simulated to match the post-operative deformity. Body weight and flexion bending (1 Nm) loads were then applied and the output stress data on the rods were analyzed.ResultsRadiographic data showed fracture locations at the mid-construct, adjacent to the distal and tandem connector across the patients. Stress analysis from the FEA showed these failure locations matched local high-stress regions for all fractures observed. These results qualitatively validate the efficacy of the FEA framework by showing a decent correlation between localized high-stress regions and the actual fracture sites in the patients.ConclusionsThis patient-specific, in-silico framework has huge potential to be used as a surgical tool to predict sites prone to fracture in growing rod implants. This prospective information would therefore be vital for surgical planning, besides helping optimize implant design for reducing rod failures.

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.

Improved Durability of a Modular Axial Fixator for Stable and Unstable Proximal Femoral Fractures: A Patient-Specific Finite Element Analysis

Abstract Femoral neck fractures, comprising 8–10% of all bodily fractures in the elderly, often necessitate alternatives to extensive surgical interventions. Despite limited research, external fixators are considered promising. This study evaluates the design and durability of a novel modular axial fixator (MAF) for stable and unstable proximal femoral fractures, using numerical method-based engineering analysis. Employing patient-specific CT scan data, three-dimensional (3D) solid modeling, and finite element analysis (FEA), the MAF-bone fixation is examined in eight simulation scenarios under static loading conditions. FEA results show a peak femur head displacement of 7.429 mm in FEA 001, with Schanz screw no. 2 reaching the maximum equivalent stress at 431.060 MPa in FEA-006. Notably, the 7.429-mm displacement improves stability compared to previous studies, yet interfragmentary movement surpasses the 100–200 μm reference range for primary fracture healing, posing challenges to direct healing despite enhanced stability. This study validates the durability of the innovative MAF for femoral neck fractures through engineering simulations. It contributes to understanding MAF durability issues, with implications for improving medical implant design in the industry. Simulation results offer opportunities for optimizing structure and production, enhancing the MAF's design, and ultimately benefiting medical implant manufacturing.

Impact of osteoporosis and Cement-Augmented fusion on adjacent spinal levels Post-Fusion Surgery: Patient-Specific finite element analysis

Cement-augmentation is a technique commonly used during posterior lumbar instrumented fusion (PLIF) to reinforce compromised osteoporotic vertebral bone, minimize the risk of loosening screws, enhance stability, and improve overall surgical outcomes. In this study, we introduce a novel segmented vertebral body regional modeling approach to investigate the effects of osteoporosis and cement-augmented lumbar fusion on disc biomechanics at spinal levels adjacent to the fused vertebrae. Using our previously validated personalized-poroelastic-osteoligamentous FE model of the spine, fusion was simulated at L4-L5, and the biomechanics of adjacent levels were studied for 30 patients (non-osteoporotic patients (N = 15), osteoporotic patients (N = 15)). PLIF models, with and without cement-augmentation, were developed and compared after an 8 h-rest period (200 N), following a 16 h-cyclic compressive loading of 500–1000 N (40 and 20 min, respectively). Movement in different directions (flexion/ extension/ lateral bending/ axial rotation) was simulated using 10Nm moment before and after cyclic loading. The material mapping algorithm was validated by comparing the results of voxel-based and parametric models. The FE cement-augmented models, subject to daily activity loading, demonstrated significant differences in disc height loss and fluid loss as compared to non-cemented models. The calculated axial stress and fiber strain values were also significantly higher for these models. This work demonstrates that although osteoporosis does not significantly alter the time-dependent characteristics of adjacent IVDs post-surgery, cement-augmentation increases the risk of adjacent segment disease (ASD) incidence. A holistic understanding of the trade-offs and long-term complex interplay between structural reinforcement modalities, including cement augmentation, and altered biomechanics warrants further investigation.

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.

FEA-Based Stress-Strain Barometers as Forecasters for Corneal Refractive Power Change in Orthokeratology.

To improve the effectivity of patient-specific finite element analysis (FEA) to predict refractive power change (RPC) in rigid Ortho-K contact lens fitting. Novel eyelid boundary detection is introduced to the FEA model to better model the effects of the lid on lens performance, and stress and strain outcomes are investigated to identify the most effective FEA components to use in modelling. The current study utilises fully anonymised records of 249 eyes, 132 right eyes, and 117 left eyes from subjects aged 14.1 ± 4.0 years on average (range 9 to 38 years), which were selected for secondary analysis processing. A set of custom-built MATLAB codes was built to automate the process from reading Medmont E300 height and distance files to processing and displaying FEA stress and strain outcomes. Measurements from before and after contact lens wear were handled to obtain the corneal surface change in shape and power. Tangential refractive power maps were constructed from which changes in refractive power pre- and post-Ortho-K wear were determined as the refractive power change (RPC). A total of 249 patient-specific FEA with innovative eyelid boundary detection and 3D construction analyses were automatically built and run for every anterior eye and lens combination while the lens was located in its clinically detected position. Maps of four stress components: contact pressure, Mises stress, pressure, and maximum principal stress were created in addition to maximum principal logarithmic strain maps. Stress and strain components were compared to the clinical RPC maps using the two-dimensional (2D) normalised cross-correlation and structural similarity (SSIM) index measure. On the one hand, the maximum principal logarithmic strain recorded the highest moderate 2D cross-correlation area of 8.6 ± 10.3%, and contact pressure recorded the lowest area of 6.6 ± 9%. Mises stress recorded the second highest moderate 2D cross-correlation area with 8.3 ± 10.4%. On the other hand, when the SSIM index was used to compare the areas that were most similar to the clinical RPC, maximum principal stress was the most similar, with an average strong similarity percentage area of 26.5 ± 3.3%, and contact pressure was the least strong similarity area of 10.3 ± 7.3%. Regarding the moderate similarity areas, all components were recorded at around 34.4% similarity area except the contact pressure, which was down to 32.7 ± 5.8%. FEA is an increasingly effective tool in being able to predict the refractive outcome of Ortho-K treatment. Its accuracy depends on identifying which clinical and modelling metrics contribute to the most accurate prediction of RPC with minimal ocular complications. In terms of clinical metrics, age, Intra-ocular pressure (IOP), central corneal thickness (CCT), surface topography, lens decentration and the 3D eyelid effect are all important for effective modelling. In terms of FEA components, maximum principal stress was found to be the best FEA barometer that can be used to predict the performance of Ortho-K lenses. In contrast, contact pressure provided the worst stress performance. In terms of strain, the maximum principal logarithmic strain was an effective strain barometer.

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.

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.