Assessment of bone ingrowth around beaded coated tibial implant for total ankle replacement using mechanoregulatory algorithm.

The long-term performance of porous coated tibial implants for total ankle replacement (TAR) primarily depends on the extent of bone ingrowth at the bone-implant interface. Although attempts were made for primary fixation for immediate post-operative stability, no investigation was conducted on secondary fixation. The aim of this study is to assess bone ingrowth around the porous beaded coated tibial implant for TAR using a mechanoregulatory algorithm. A realistic macroscale finite element (FE) model of the implanted tibia was developed based on computer tomography (CT) data to assess implant-bone micromotions and coupled with microscale FE models of the implant-bone interface to predict bone ingrowth around tibial implant for TAR. The macroscale FE model was subjected to three near physiological loading conditions to evaluate the site-specific implant-bone micromotion, which were then incorporated into the corresponding microscale model to mimic the near physiological loading conditions. Results of the study demonstrated that the implant experienced tangential micromotion ranged from 0 to 71 µm with a mean of 3.871 µm. Tissue differentiation results revealed that bone ingrowth across the implant ranged from 44 to 96 %, with a mean of around 70 %. The average Young's modulus of the inter-bead tissue layer varied from 1444 to 4180 MPa around the different regions of the implant. The analysis postulates that when peak micromotion touches 30 µm around different regions of the implant, it leads to pronounced fibrous tissues on the implant surface. The highest amount of bone ingrowth was observed in the central regions, and poor bone ingrowth was seen in the anterior parts of the implant, which indicate improper osseointegration around this region. This macro-micro mechanical FE framework can be extended to improve the implant design to enhance the bone ingrowth and in future to develop porous lattice-structured implants to predict and enhance osseointegration around the implant.

Biomechanical assessment of different transforaminal lumbar interbody fusion constructs in normal and osteoporotic condition: a finite element analysis.

BACKGROUND CONTEXT: With the aging population, osteoporosis, which leads to poor fusion, has become a common challenge for lumbar surgery. In addition, most people with osteoporosis are elderly individuals with poor surgical tolerance, and poor bone quality can also weaken the stability of internal fixation. PURPOSE: This study compared the fixation strength of the bilateral traditional trajectory screw structure (TT-TT), the bilateral cortical bone trajectory screw structure (CBT-CBT), and the hybrid CBT-TT (CBT screws at the cranial level and TT screws at the caudal level) structure under different bone mineral density conditions. STUDY DESIGN: A finite element (FE) analysis study. METHODS: Above all, we established a healthy adult lumbar spine model. Second, under normal and osteoporotic conditions, three transforaminal lumbar interbody fusion (TLIF) models were established: bilateral traditional trajectory (TT-TT) screw fixation, bilateral cortical bone trajectory (CBT-CBT) screw fixation, and hybrid cortical bone trajectory screw and traditional trajectory screw (CBT-TT) fixation. Finally, a 500-N compression load with a torque of 10 N/m was applied to simulate flexion, extension, lateral bending, and axial rotation. We compared the range of motion (ROM), adjacent disc stress, cage stress, and posterior fixation stress of the different fusion models. RESULTS: Under different bone mineral density conditions, the range of motion of the fusion segment was significantly reduced. Compared to normal bone conditions, the ROM of the L4-L5 segment, the stress of the adjacent intervertebral disc, the surface stress of the cage, and the maximum stress of the posterior fixation system were all increased in osteoporosis. Under most loads, the ROM and surface stress of the cage and the maximum stress of the posterior fixation system of the TT-TT structure are the lowest under normal bone mineral density conditions. However, under osteoporotic conditions, the fixation strength of the CBT-CBT and CBT-TT structures are higher than that of the TT-TT structures under certain load conditions. At the same time, the surface stress of the intervertebral fusion cage and the maximum stress of the posterior fixation system for the two structures are lower than those of the TT-TT structure. CONCLUSION: Under normal bone mineral density conditions, transforaminal lumbar interbody fusion combined with TT-TT fixation provides the best biomechanictability. However, under osteoporotic conditions, CBT-CBT and CBT-TT structures have higher fixed strength compared to TT-TT structures. The hybrid CBT-TT structure exhibits advantages in minimal trauma and fixation strength. Therefore, this seems to be an alternative fixation method for patients with osteoporosis and degenerative spinal diseases. CLINICAL SIGNIFICANCE: This study provides biomechanical support for the clinical application of hybrid CBT-TT structure for osteoporotic patients undergoing TLIF surgery.

Fracture risk assessment of vascularized medial femoral condylar bone graft: A finite element analysis.

BACKGROUND: Vascularized medial femoral condyle (MFC) bone graft is useful for pseudarthrosis and osteonecrosis, but has the risk of fracture as a complication. This study aimed to create multiple three-dimensional (3D) finite element (FE) femur models to biomechanically evaluate the fracture risk in the donor site of a vascularized MFC bone graft. METHODS: Computer tomography scans of the femurs of nine patients (four males and five females) with no left femur disease were enrolled in the study. A 3D FE model of the left femur was generated based on the CT images taken from the patients. The descending genicular artery (DGA), the main nutrient vessel in vascularized MFC bone grafts, divides into the proximal transversal branch (TB) and the distal longitudinal branch (LB) before entering the periosteum. Thirty-six different bone defect models with different sizes and locations of the harvested bone were created. RESULTS: The highest stress was observed in the proximal medial and metaphyseal portions under axial and external rotation, respectively. In the bone defect model, the stress was most elevated in the extracted region's anterior or posterior superior part. Stress increased depending on proximal location and harvested bone size. CONCLUSION: Increasing the size of the bone graft proximally raises the stress at the site of bone extraction. For bone grafting to non-load-bearing areas, bone grafting distally using LB can reduce fracture risk. If TB necessitates a larger proximal bone extraction, it is advisable to avoid postoperative rotational loads.

Comparative assessment of orthodontic clear aligner versus fixed appliance for anterior retraction: a finite element study.

BACKGROUND: The aim of this study is to conduct a comparative evaluation of different designs of clear aligners and examine the disparities between clear aligners and fixed appliances. METHODS: 3D digital models were created, consisting of a maxillary dentition without first premolars, maxilla, periodontal ligaments, attachments, micro-implant, 3D printed lingual retractor, brackets, archwire and clear aligner. The study involved the creation of five design models for clear aligner maxillary anterior internal retraction and one design model for fixed appliance maxillary anterior internal retraction, which were subsequently subjected to finite element analysis. These design models included: (1) Model C0 Control, (2) Model C1 Posterior Micro-implant, (3) Model C2 Anterior Micro-implant, (4) Model C3 Palatal Plate, (5) Model C4 Lingual Retractor, and (6) Model F0 Fixed Appliance. RESULTS: In the clear aligner models, a consistent pattern of tooth movement was observed. Notably, among all tested models, the modified clear aligner Model C3 exhibited the smallest differences in sagittal displacement of the crown-root of the central incisor, vertical displacement of the central incisor, sagittal displacement of the second premolar and second molar, as well as vertical displacement of posterior teeth. However, distinct variations in tooth movement trends were observed between the clear aligner models and the fixed appliance model. Furthermore, compared to the fixed appliance model, significant increases in tooth displacement were achieved with the use of clear aligner models. CONCLUSIONS: In the clear aligner models, the movement trend of the teeth remained consistent, but there were variations in the amount of tooth displacement. Overall, the Model C3 exhibited better torque control and provided greater protection for posterior anchorage teeth compared to the other four clear aligner models. On the other hand, the fixed appliance model provides superior anterior torque control and better protection of the posterior anchorage teeth compared to clear aligner models. The clear aligner approach and the fixed appliance approach still exhibit a disparity; nevertheless, this study offers a developmental direction and establishes a theoretical foundation for future non-invasive, aesthetically pleasing, comfortable, and efficient modalities of clear aligner treatment.

Macro-indentation testing of soft biological materials and assessment of hyper-elastic material models from inverse finite element analysis.

Mechanical characterization of hydrogels and ultra-soft tissues is a challenging task both from an experimental and material parameter estimation perspective because they are much softer than many biological materials, ceramics, or polymers. The elastic modulus of such materials is within the 1 - 100 kPa range, behaving as a hyperelastic solid with strain hardening capability at large strains. In the current study, indentation experiments have been performed on agarose hydrogels, bovine liver, and bovine lymph node specimens. This work reports on the reliable determination of the elastic modulus by indentation experiments carried out at the macro-scale (mm) using a spherical indenter. However, parameter identification of the hyperelastic material properties usually requires an inverse finite element analysis due to the lack of an analytical contact model of the indentation test. Hence a comprehensive study on the spherical indentation of hyperelastic soft materials is carried out through robust computational analysis. Neo-Hookean and first-order Ogden hyperelastic material models were found to be most suitable. A case study on known anisotropic hyperelastic material showed the inability of the inverse finite element method to uniquely identify the whole material parameter set.

Assessment of biomechanical behavior of immature non-vital incisors with various treatment modalities by means of three-dimensional quasi-static finite element analysis.

The objectives of this study were to evaluate the stress distribution and risk of fracture of a non-vital immature maxillary central incisor subjected to various clinical procedures using finite element analysis (FEA). A three-dimensional model of an immature central incisor was developed, from which six main models were designed: untreated immature tooth (C), standard apical plug (AP), resin composite (RC), glass-fibre post (GFP), regeneration procedure (RET), and regeneration with induced root maturation (RRM). Mineral trioxide aggregate (MTA) or Biodentine® were used as an apical or coronal plug. All models simulated masticatory forces in a quasi-static approach with an oblique force of 240 Newton at a 120° to the longitudinal tooth axis. The maximum principal stress, maximum shear stress, risk of fracture, and the strengthening percentage were evaluated. The mean maximum principal stress values were highest in model C [90.3 MPa (SD = 4.4)] and lowest in the GFP models treated with either MTA and Biodentine®; 64.1 (SD = 1.7) and 64.0 (SD = 1.6) MPa, respectively. Regarding the shear stress values, the dentine tooth structure in model C [14.4 MPa (SD = 0.8)] and GFP models [15.4 MPa (SD = 1.1)] reported significantly higher maximum shear stress values compared to other tested models (p < 0.001), while no significant differences were reported between the other models (p > 0.05). No significant differences between MTA and Biodentine® regarding maximum principal stress and maximum shear stress values for each tested model (p > 0.05). A maximum strain value of 4.07E-03 and maximum displacement magnitude of 0.128 mm was recorded in model C. In terms of strengthening percentage, the GFP models were associated with the highest increase (22%). The use of a GFP improved the biomechanical performance and resulted in a lower risk of fracture of a non-vital immature maxillary central incisor in a FEA model.

Assessment of the elastic stiffness of human cardiac fibres after an apical infarction using finite element simulation.

Several research works in the literature have focused on understanding the post-infarction ventricular remodelling phenomenon, but few works have considered the evaluation of the elastic behaviour of the cardiac tissue after a myocardial infarction. This paper presents an investigation focused on predicting the elastic performance of the human heart after a left ventricular apical infarction. The aim is to understand the elastic alterations of the cardiac fibres at different periods after an apical infarct. For this purpose, a hybrid method based on pressure and volume measurements of the left ventricle (LV) at different periods of ventricular remodelling, and the Finite Element Method (FEM), is developed. In addition, several performance indexes are defined to evaluate the heart performance during the ventricular remodelling process. The results show that during the first 2 weeks after a heart infarction, the cardiac fibres must support a much higher structural overload than during normal conditions. This structural overload is proportional to the aneurysm size but diminishes with the time, together with a significant reduction of the ventricular pumping capacity.

Assessment of finite element human body and ATD models in estimating injury risk in far-side impacts using field-based injury risk.

The objective of this study was to assess the ability of finite element human body models (FEHBMs) and Anthropometric Test Device (ATD) models to estimate occupant injury risk by comparing it with field-based injury risk in far-side impacts. The study used the Global Human Body Models Consortium midsize male (M50-OS+B) and small female (F05-OS+B) simplified occupant models with a modular detailed brain, and the ES-2Re and SID-IIs ATD models in the simulated far-side crashes. A design of experiments (DOE) with a total of 252 simulations was conducted by varying lateral ΔV (10-50kph; 5kph increments), the principal direction of force (PDOF 50°, 60°, 65°, 70°, 75°, 80°, 90°), and occupant models. Models were gravity-settled and belted into a simplified vehicle model (SVM) modified for far-side impact simulations. Acceleration pulses and vehicle intrusion profiles used for the DOE were generated by impacting a 2012 Camry vehicle model with a mobile deformable barrier model across the 7 PDOFs and 9 lateral ΔV's in the DOE for a total of 63 additional simulations. Injury risks were estimated for the head, chest, lower extremity, pelvis (AIS 2+; AIS 3+), and abdomen (AIS 3+) using logistic regression models. Combined AIS 3+ injury risk for each occupant was calculated using AIS 3+ injury risk estimations for the head, chest, abdomen, and lower extremities. The injury risk calculated using computational models was compared with field-based injury risk derived from NASS-CDS by calculating their correlation coefficient. The field-based injury risk was calculated using risk curves that were created based on real-world crash data in a previous study (Hostetler et al., 2020). Occupant age (40 years), seatbelt use (belted occupant), collision deformation classification, lateral ΔV, and PDOF of the crash event were used in these curves to estimate field injury risk. Large differences in the kinematics were observed between HBM and ATD models. ATD models tended to overestimate risk in almost every case whereas HBMs yielded better risk estimates overall. Chest and lower extremity risks were the least correlated with field injury risk estimates. The overall risk of AIS 3+ injury risk was the strongest comparison to the field data-based risk curves. The HBMs were still not able to capture all the variance but future studies can be carried out that are focused on investigating their shortfalls and improving them to estimate injury risk closer to field injury risk in far-side crashes.

Computational assessment of leg response to extreme loadings using a detailed finite element model.

This study focuses on evaluating the response of the Total Human Model for Safety™ lower extremity finite element model under blast loading. Biofidelity of the lower extremity model was evaluated against experiments with impact loading equivalent to underbody blast. The model response was found to match well with the experimental data for the average impactor speeds of 7 and 9.3 m/s resulting in an overall correlation and analysis rating of 0.86 and 0.82, respectively. The model response was then used to investigate response for antipersonnel mine explosion where the numerical setup consists of a charge mass of 40 g trinitrotoluene placed at a depth of 50 mm below the heel. The explosion was modeled using Multi Material-Arbitrary Lagrangian Eulerian method. The model was subjected to the graded input in terms of variation in standoff distance and mass of explosive to investigate the sensitivity of the model. The model found sensitive to the threat definition and predicted an increase of 110% in peak fluid-structure interaction force with 20% reduction in its time to peak and 29% increase in peak calcaneus axial force with a reduction of 33% in its time to peak when explosive mass varied from 40 g to 100 g. The location of the explosive below the foot was discovered to have significant effect on the injury pattern in near-field explosion. A comparative study suggested that the model predicted similar response and damage pattern compared to experimental data.

Cortical and Trabecular Bone Stress Assessment during Periodontal Breakdown-A Comparative Finite Element Analysis of Multiple Failure Criteria.

Background and Objectives: This numerical analysis investigated the biomechanical behavior of the mandibular bone as a structure subjected to 0.5 N of orthodontic force during periodontal breakdown. Additionally, the suitability of the five most used failure criteria (Von Mises (VM), Tresca (T), maximum principal (S1), minimum principal (S3), and hydrostatic pressure (HP)) for the study of bone was assessed, and a single criterion was identified for the study of teeth and the surrounding periodontium (by performing correlations with other FEA studies). Materials and Methods: The finite element analysis (FEA) employed 405 simulations over eighty-one mandibular models with variable levels of bone loss (0-8 mm) and five orthodontic movements (intrusion, extrusion, tipping, rotation, and translation). For the numerical analysis of bone, the ductile failure criteria are suitable (T and VM are adequate for the study of bone), with Tresca being more suited. S1, S3, and HP criteria, due to their distinctive design dedicated to brittle materials and liquids/gas, only occasionally correctly described the bone stress distribution. Results: Only T and VM displayed a coherent and correlated gradual stress increase pattern for all five movements and levels of the periodontal breakdown. The quantitative values provided by T and VM were the highest (for each movement and level of bone loss) among all five criteria. The MHP (maximum physiological hydrostatic pressure) was exceeded in all simulations since the mandibular bone is anatomically less vascularized, and the ischemic risks are reduced. Only T and VM displayed a correlated (both qualitative and quantitative) stress increase for all five movements. Both T and VM displayed rotation and translation, closely followed by tipping, as stressful movements, while intrusion and extrusion were less stressful for the mandibular bone. Conclusions: Based on correlations with earlier numerical studies on the same models and boundary conditions, T seems better suited as a single unitary failure criterion for the study of teeth and the surrounding periodontium.

The impact of modeling choices on the assessment of Ni-Ti fatigue properties through surrogate specimens.

The implant of self-expandable Ni-Ti stents for the treatment of peripheral diseases has become an established medical practice. However, the reported failure in clinics highlights the open issue of the fatigue characterization of these devices. One of the most common approaches for calculating the Ni-Ti fatigue limit (commonly defined in terms of mean and alternate strain for a fixed number of cycles) consists of using surrogate specimens which replicate the strain distributions of the final device but in simplified geometries. The main drawback lies in the need for computational models to determine the local distribution and, hence, interpret the experimental results. This study aims at investigating the role of different choices in the model preparation, such as the mesh refinement and the element formulation, on the output of the fatigue analysis. The analyses show a strong dependency of the numerical results on modeling choices. The use of linear reduced elements enriched by a layer of membrane elements is successful to increase the accuracy of the results, especially when coarser meshes are used. Due to material nonlinearity and stent complex geometries, for the same loading conditions and element type, (i) different meshes result in different couples of mean and amplitude strains and (ii) for the same mesh, the position of the maximum mean strain is not coincident with the maximum amplitude, making difficult the selection of the limit values.

Computational assessment on the impact of collagen fiber orientation in cartilages on healthy and arthritic knee kinetics/kinematics.

BACKGROUND: The inhomogeneous distribution of collagen fiber in cartilage can substantially influence the knee kinematics. This becomes vital for understanding the mechanical response of soft tissues, and cartilage deterioration including osteoarthritis (OA). Though the conventional computational models consider geometrical heterogeneity along with fiber reinforcements in the cartilage model as material heterogeneity, the influence of fiber orientation on knee kinetics and kinematics is not fully explored. This work examines how the collagen fiber orientation in the cartilage affects the healthy (intact knee) and arthritic knee response over multiple gait activities like running and walking. METHODS: A 3D finite element knee joint model is used to compute the articular cartilage response during the gait cycle. A fiber-reinforced porous hyper elastic (FRPHE) material is used to model the soft tissue. A split-line pattern is used to implement the fiber orientation in femoral and tibial cartilage. Four distinct intact cartilage models and three OA models are simulated to assess the impact of the orientation of collagen fibers in a depth wise direction. The cartilage models with fibers oriented in parallel, perpendicular, and inclined to the articular surface are investigated for multiple knee kinematics and kinetics. FINDINGS: The comparison of models with fiber orientation parallel to articulating surface for walking and running gait has the highest elastic stress and fluid pressure compared with inclined and perpendicular fiber-oriented models. Also, the maximum contact pressure is observed to be higher in the case of intact models during the walking cycle than for OA models. In contrast, the maximum contact pressure is higher during running in OA models than in intact models. Additionally, parallel-oriented models produce higher maximum stresses and fluid pressure for walking and running gait than proximal-distal-oriented models. Interestingly, during the walking cycle, the maximum contact pressure with intact models is approximately three times higher than on OA models. In contrast, the OA models exhibit higher contact pressure during the running cycle. INTERPRETATION: Overall, the study indicates that collagen orientation is crucial for tissue responsiveness. This investigation provides insights into the development of tailored implants.

Assessment of Bone Healing: Opportunities to Improve the Standard of Care.

➤ Bone healing is commonly evaluated by clinical examination and serial radiographic evaluation. Physicians should be mindful that personal and cultural differences in pain perception may affect the clinical examination. Radiographic assessment, even with the Radiographic Union Score, is qualitative, with limited interobserver agreement.➤ Physicians may use serial clinical and radiographical examinations to assess bone healing in most patients, but in ambiguous and complicated cases, they may require other methods to provide assistance in decision-making.➤ In complicated instances, clinically available biomarkers, ultrasound, and magnetic resonance imaging may determine initial callus development. Quantitative computed tomography and finite element analysis can estimate bone strength in later callus consolidation phases.➤ As a future direction, quantitative rigidity assessments for bone healing may help patients to return to function earlier by increasing a clinician's confidence in successful progressive healing.

Hip Fracture Risk Assessment in Elderly and Diabetic Patients: Combining Autonomous Finite Element Analysis and Machine Learning.

Autonomous finite element analyses (AFE) based on CT scans predict the biomechanical response of femurs during stance and sidewise fall positions. We combine AFE with patient data via a machine learning (ML) algorithm to predict the risk of hip fracture. An opportunistic retrospective clinical study of CT scans is presented, aimed at developing a ML algorithm with AFE for hip fracture risk assessment in type 2 diabetic mellitus (T2DM) and non-T2DM patients. Abdominal/pelvis CT scans of patients who experienced a hip fracture within 2 years after an index CT scan were retrieved from a tertiary medical center database. A control group of patients without a known hip fracture for at least 5 years after an index CT scan was retrieved. Scans belonging to patients with/without T2DM were identified from coded diagnoses. All femurs underwent an AFE under three physiological loads. AFE results, patient's age, weight, and height were input to the ML algorithm (support vector machine [SVM]), trained by 80% of the known fracture outcomes, with cross-validation, and verified by the other 20%. In total, 45% of available abdominal/pelvic CT scans were appropriate for AFE (at least 1/4 of the proximal femur was visible in the scan). The AFE success rate in automatically analyzing CT scans was 91%: 836 femurs we successfully analyzed, and the results were processed by the SVM algorithm. A total of 282 T2DM femurs (118 intact and 164 fractured) and 554 non-T2DM (314 intact and 240 fractured) were identified. Among T2DM patients, the outcome was: Sensitivity 92%, Specificity 88% (cross-validation area under the curve [AUC] 0.92) and for the non-T2DM patients: Sensitivity 83%, Specificity 84% (cross-validation AUC 0.84). Combining AFE data with a ML algorithm provides an unprecedented prediction accuracy for the risk of hip fracture in T2DM and non-T2DM populations. The fully autonomous algorithm can be applied as an opportunistic process for hip fracture risk assessment. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Assessment of the Maximum Amount of Orthodontic Force for PDL in Intact and Reduced Periodontium (Part I).

This study examines 0.6 N and 1.2 N as the maximum orthodontic force for periodontal ligament (PDL) at multiple levels of periodontal breakdown, and the relationships with the ischemic, necrotic, and resorptive risks. Additionally, this study evaluates if Tresca failure criteria is more adequate for the PDL study. Eighty-one 3D models (from nine patients; nine models/patients) with the 2nd lower premolar and different degrees of bone loss (0-8 mm) where subjected to intrusion, extrusion, rotation, translation, and tipping movements. Tresca shear stress was assessed individually for each movement and bone loss level. Rotation and translation produced the highest PDL stresses, while intrusion and extrusion determined the lowest. Apical and middle third PDL stresses were lower than the cervical stress. In intact periodontium, the amount of shear stress produced by the two investigated forces was lower than the 16 KPa of the maximum physiological hydrostatic pressure (MHP). In reduced periodontium (1-8 mm tissue loss), the apical amount of PDL shear stress was lower than MHP for both applied forces, while cervically for rotation, translation and tipping movements exceeded 16 KPa. Additionally, 1.2 N could be used in intact periodontium (i.e., without risks) and for the reduced periodontium only in the apical and middle third of PDL up to 8 mm of bone loss. However, for avoiding any resorptive risks, in the cervical third of PDL, the rotation, translation, and tipping movements require less than 0.2-0.4 N of force after 4 mm of loss. Tresca seems to be more adequate for the study of PDL than other criteria.

Assessment of the Maximum Amount of Orthodontic Force for Dental Pulp and Apical Neuro-Vascular Bundle in Intact and Reduced Periodontium on Bicuspids (Part II).

This study examines 0.6 N-4.8 N as the maximum orthodontic force to be applied to dental pulp and apical NVB on intact and 1-8 mm reduced periodontal-ligament (PDL), in connection with movement and ischemic, necrotic and resorptive risk. In addition, it examines whether the Tresca finite-element-analysis (FEA) criterion is more adequate for the examination of dental pulp and its apical NVB. Eighty-one (nine patients, with nine models for each patient) anatomically correct models of the periodontium, with the second lower-premolar reconstructed with its apical NVB and dental pulp were assembled, based on X-ray CBCT (cone-beam-computed-tomography) examinations and subjected to 0.6 N, 1.2 N, 2.4 N and 4.8 N of intrusion, extrusion, translation, rotation, and tipping. The Tresca failure criterion was applied, and the shear stress was assessed. Forces of 0.6 N, 1.2 N, and 2.4 N had negligible effects on apical NVB and dental pulp up to 8 mm of periodontal breakdown. A force of 4.8 N was safely applied to apical NVB on the intact periodontium only. Rotation and tipping seemed to be the most invasive movements for the apical NVB. For the dental pulp, only the translation and rotation movements seemed to display a particular risk of ischemia, necrosis, and internal orthodontic-resorption for both coronal (0-8 mm of loss) and radicular pulp (4-8 mm of loss), despite the amount of stress being lower than the MHP. The Tresca failure criterion seems more suitable than other criteria for apical NVB and dental pulp.

Impact of secondary particles on the magnetic field generated by a proton pencil beam: a finite-element analysis based on Geant4-DNA simulations.

PURPOSE: To investigate the static magnetic field generated by a proton pencil beam as a candidate for range verification by means of Monte Carlo simulations, thereby improving upon existing analytical calculations. We focus on the impact of statistical current fluctuations and secondary protons and electrons. METHODS: We considered a pulsed beam (10 µ ${\umu}$ s pulse duration) during the duty cycle with a peak beam current of 0.2 µ $\umu$ A and an initial energy of 100 MeV. We ran Geant4-DNA Monte Carlo simulations of a proton pencil beam in water and extracted independent particle phase spaces. We calculated longitudinal and radial current density of protons and electrons, serving as an input for a magnetic field estimation based on a finite element analysis in a cylindrical geometry. We made sure to allow for non-solenoidal current densities as is the case of a stopping proton beam. RESULTS: The rising proton charge density toward the range is not perturbed by energy straggling and only lowered through nuclear reactions by up to 15%, leading to an approximately constant longitudinal current. Their relative low density however (at most 0.37 protons/mm3 for the 0.2  µ ${\umu}$ A current and a beam cross-section of 2.5 mm), gives rise to considerable current density fluctuations. The radial proton current resulting from lateral scattering and being two orders of magnitude weaker than the longitudinal current is subject to even stronger fluctuations. Secondary electrons with energies above 10 eV, that far outnumber the primary protons, reduce the primary proton current by only 10% due to their largely isotropic flow. A small fraction of electrons (<1%), undergoing head-on collisions, constitutes the relevant electron current. In the far-field, both contributions to the magnetic field strength (longitudinal and lateral) are independent of the beam spot size. We also find that the nuclear reaction-related losses cause a shift of 1.3 mm to the magnetic field profile relative to the actual range, which is further enlarged to 2.4 mm by the electron current (at a distance of ρ = 50 $\rho =50$  mm away from the central beam axis). For ρ > 45 $\rho >45$  mm, the shift increases linearly. While the current density variations cause significant magnetic field uncertainty close to the central beam axis with a relative standard deviation (RSD) close to 100%, they average out at a distance of 10 cm, where the RSD of the total magnetic field drops below 2%. CONCLUSIONS: With the small influence of the secondary electrons together with the low RSD, our analysis encourages an experimental detection of the magnetic field through sensitive instrumentation, such as optical magnetometry or SQUIDs.

Fracture risk assessment and evaluation of femoroplasty in metastatic proximal femurs. An in vivo CT-based finite element study.

The goal of this study was twofold. First, we aimed to evaluate the accuracy of a finite element (FE) model to predict bone fracture in cancer patients with proximal femoral bone metastases. Second, we evaluated whether femoroplasty could effectively reduce fracture risk. A total of 89 patients were included, with 101 proximal femurs affected with bone metastases. The accuracy of the model to predict fracture was evaluated by comparing the FE failure load, normalized for body weight, against the actual occurrence of fracture during a 6-month follow-up. Using a critical threshold, the model could identify whether femurs underwent fracture with a sensitivity of 92% and a specificity of 66%. A virtual treatment with femoroplasty was simulated in a subset of 34 out of the 101 femurs; only femurs with one or more well-defined lytic lesions were considered eligible for femoroplasty. We modeled their lesions, as well as the surrounding 4 mm of trabecular bone, to be augmented with bone cement. The simulation of femoroplasty increased the median failure load of the FE model by 57% for lesions located in the head/neck of the femur. At this lesion location, all high risk femurs that had fractured during follow-up effectively moved from a failure load below the critical threshold to a value above. For lesions located in the trochanteric region, no definite improvement in failure load was found. Although additional validation studies are required, our results suggest that femoroplasty can effectively reduce fracture risk for several osteolytic lesions in the femoral head/neck.

Biomechanical assessment of disease outcome in surgical interventions for medial meniscal posterior root tears: a finite element analysis.

BACKGROUND: The adverse consequences of medial meniscus posterior root tears have become increasingly familiar to surgeons, and treatment strategies have become increasingly abundant. In this paper, the finite element gait analysis method was used to explore the differences in the biomechanical characteristics of the knee joint under different conditions. METHODS: Based on CT computed tomography and MR images, (I) an intact knee (IK) model with bone, cartilage, meniscus and main ligaments was established. Based on this model, the posterior root of the medial meniscus was resected, and (ii) the partial tear (PT) model, (iii) the entire radial tear (ERT) model, and (iv) the entire oblique tear (EOT) model were established according to the scope and degree of resection. Then, the (v) meniscus repair (MR) model and (vi) partial meniscectomy (PM) model were developed according to the operation method. The differences in stress, displacement and contact area among different models were evaluated under ISO gait loading conditions. RESULTS: Under gait loading, there was no significant difference in the maximum stress of the medial and lateral tibiofemoral joints among the six models. Compared with the medial tibiofemoral joint stress of the IK model, the stress of the PM model increased by 8.3%, while that of the MR model decreased by 18.9%; at the same time, the contact stress of the medial tibiofemoral joint of the ERT and EOT models increased by 17.9 and 25.3%, respectively. The displacement of the medial meniscus in the ERT and EOT models was significantly larger than that in the IK model (P < 0.05), and the tibial and femoral contact areas of these two models were lower than those of the IK model (P < 0.05). CONCLUSIONS: The integrity of the posterior root of the medial meniscus plays an important role in maintaining normal tibial-femoral joint contact mechanics. Partial meniscectomy is not beneficial for improving the tibial-thigh contact situation. Meniscal repair has a positive effect on restoring the normal biomechanical properties of the medial meniscus.

Assessment of the Best FEA Failure Criteria (Part II): Investigation of the Biomechanical Behavior of Dental Pulp and Apical-Neuro-Vascular Bundle in Intact and Reduced Periodontium.

The aim of this study was to biomechanically assess the behavior of apical neuro-vascular bundles (NVB) and dental pulp employing Tresca, Von Mises, Pressure, S1 and S3 failure criterions in a gradual periodontal breakdown under orthodontic movements. Additionally, it was to assess the accuracy of failure criteria, correlation with the maximum hydrostatic pressure (MHP), and the amount of force safe for reduced periodontium. Based on cone-beam computed tomography, 81 3D models of the second lower premolar were subjected to 0.5 N of intrusion, extrusion, rotation, tipping, and translation. A Finite Elements Analysis (FEA) was performed. In intact and reduced periodontium apical NVB, stress (predominant in all criteria) was significantly higher than dental pulp stress, but lower than MHP. VM and Tresca displayed identical results, with added pulpal stress in translation and rotation. S1, S3 and Pressure showed stress in the apical NVB area. 0.5 N seems safe up to 8 mm periodontal breakdown. A clear difference between failure criteria for dental pulp and apical NVB cannot be proved based only on the correlation quantitative results-MHP. Tresca and VM (adequate for ductile materials) showed equivalent results with the lowest amounts of stress. The employed failure criteria must be selected based on the type of material to be analyzed.