Considerations for reporting finite element analysis studies in biomechanics

Topology Optimization Driven Bone-Remodeling Simulation for Lumbar Interbody Fusion

This study proposes a numerical approach for simulating bone remodeling in lumbar interbody fusion (LIF). It employs a topology optimization method to drive the remodeling process and uses a pixel function to describe the structural topology and bone density distribution. Unlike traditional approaches based on strain energy density or compliance, this study adopts von Mises stress to guide the remodeling of LIF. A novel pixel interpolation scheme associated with stress criteria is applied to the physical properties of the bone, directly addressing the stress shielding effect caused by the implanted cage, which significantly influences the bone remodeling outcome in LIF. Additionally, a boundary inverse approach is utilized to reconstruct a simplified analysis model. To reduce computational cost while maintaining high structural resolution and accuracy, the scaled boundary finite element method (SBFEM) is introduced. The proposed numerical approach successfully generates results that closely resemble human lumbar interbody fusion.

Finite element analysis in the Dental Sciences: A Bibliometric and a Visual Study

INTRODUCTION AND AIMS

Finite element analysis (FEA) is an incrementally practical and precise tool for the prediction of stress effects on different tissue structures and has therefore interested dental researchers for decades. This bibliometric and visualized study was aimed to assess the research progress related to FEA in the dental sciences in terms of research trends and frontiers.

METHODS

The articles about FEA studies in this field during 1999 to 2024 were obtained from Web of Science Core Collection. Then, these results were analysed and plotted using Microsoft Excel, VOSviewer, and CiteSpace in order to find out the historical evolution, current hotspots, and future directions.

RESULTS

Total 2838 literature records related to the topic were retrieved from Web of Science Core Collection. The most active country and institution were USA (538 documents) and Universidade Estadual Paulista (140 documents), respectively. Baggi et al from University of Naples Federico II was the author with the most highly cited article (352 citations), which was published on the Journal of Prosthetic Dentistry in 2008. Dental Materials ranked first (231 documents) among the 10 journals with the greatest numbers of relevant publications. The top three trending keywords were 'dental implant', 'stress distribution', and 'fracture'. The endocrown, clear aligner, and posterior edentulism were scientific frontiers in this field.

CONCLUSION

The present study provides a comprehensive bibliometric analysis of research in the dental science by FEA approaches, which will identify active hotspots of scientific interest to guide further research endeavours.

Kinematic-kinetic compliant acetabular cup positioning based on preoperative motion tracking and musculoskeletal modeling for total hip arthroplasty

The invention of the surgical robot enabled accurate component implantation during total hip arthroplasty (THA). However, a preoperative surgical planning methodology is still lacking to determine the acetabular cup alignment considering the patient-specific hip functions during daily activities such as walking. To simultaneously avoid implant edgeloading and impingement, this study established a kinematic-kinetic compliant (KKC) acetabular cup positioning method based on preoperative gait kinematics measurement and musculoskeletal modeling. Computed tomography images around the hip joint and their biomechanical data during gait, including motion tracking and foot-ground reaction forces, were collected. Using the reconstructed pelvic and femur geometries, the patient-specific hip muscle insertions were located in the lower limb musculoskeletal model via point cloud registration. The designed cup orientation has to be within the patient-specific safe zone to prevent implant impingement, and the optimized value selected based on the time-dependent hip joint reaction force to minimize the risk of edgeloading. As a validation of the proposed musculoskeletal model, the predicted lower limb muscle activations for seven patients were correlated with their surface electromyographic measurements, and the computed hip contact force was also in quantitative agreement with data from the literature. However, the designed cup orientations were not always within the well-known Lewinnek safe zone, highlighting the importance of KKC surgical planning based on patient-specific biomechanical evaluations.

Comparative biomechanical study of different screw fixation methods for minimally invasive hallux valgus surgery: A finite element analysis

BACKGROUND

There are different screw configurations utilised for minimally invasive hallux valgus (HV) deformity despite limited biomechanical data assessing the stability and strength of each construct. We aimed to compare the strength of various screw configurations for minimally invasive HV surgery using finite element analysis (FEA).

METHODS

A FEA model was developed from a CT of a female with moderate HV deformity. Five screw configurations utilizing one or two bicortical or intramedullary screws were tested. Stress analysis considered osteotomy displacement, maximum and minimum principal stresses, and von Mises stress for both implants and bone for each screw configuration.

RESULTS

Fixation with two screws (one bicortical and one intramedullary) demonstrated the lowest values for osteotomy displacement, minimum and maximum total stress, and equivalent von Mises stress on the bone and screws in both loading conditions.

CONCLUSION

The optimal configuration when performing minimally invasive surgery for moderate HV is one bicortical and one intramedullary screw.

LEVEL OF EVIDENCE

Level III.

Simulation studies of social systems: telling the story based on provenance patterns

Social simulation studies are complex. They typically combine various data sources and hypotheses about the system's mechanisms that are integrated by intertwined processes of model building, simulation experiment execution and analysis. Various documentation approaches exist to increase the transparency and traceability of complex social simulation studies. Provenance standards enable the formalization of information on sources and activities, which contribute to the generation of an entity, in a queryable and computationally accessible manner. Provenance patterns can be defined as constraints on the relationships between specific types of activities and entities of a simulation study. In this paper, we refine the provenance pattern-based approach to address specific challenges of social agent-based simulation studies. Specifically, we focus on the activities and entities involved in collecting and analysing primary data about human decisions, and the collection and quality assessment of secondary data. We illustrate the potential of this approach by applying it to central activities and results of an agent-based simulation project and by presenting its implementation in a web-based tool.

Advancements in diabetic foot insoles: a comprehensive review of design, manufacturing, and performance evaluation

The escalating prevalence of diabetes has accentuated the significance of addressing the associated diabetic foot problem as a major public health concern. Effectively offloading plantar pressure stands out as a crucial factor in preventing diabetic foot complications. This review comprehensively examines the design, manufacturing, and evaluation strategies employed in the development of diabetic foot insoles. Furthermore, it offers innovative insights and guidance for enhancing their performance and facilitating clinical applications. Insoles designed with total contact customization, utilizing softer and highly absorbent materials, as well as incorporating elliptical porous structures or triply periodic minimal surface structures, prove to be more adept at preventing diabetic foot complications. Fused Deposition Modeling is commonly employed for manufacturing; however, due to limitations in printing complex structures, Selective Laser Sintering is recommended for intricate insole designs. Preceding clinical implementation, in silico and in vitro testing methodologies play a crucial role in thoroughly evaluating the pressure-offloading efficacy of these insoles. Future research directions include advancing inverse design through machine learning, exploring topology optimization for lightweight solutions, integrating flexible sensor configurations, and innovating new skin-like materials tailored for diabetic foot insoles. These endeavors aim to further propel the development and effectiveness of diabetic foot management strategies. Future research avenues should explore inverse design methodologies based on machine learning, topology optimization for lightweight structures, the integration of flexible sensors, and the development of novel skin-like materials specifically tailored for diabetic foot insoles. Advancements in these areas hold promise for further enhancing the effectiveness and applicability of diabetic foot prevention measures.

Developing Advanced Patient-Specific In Silico Models: A New Era in Biomechanical Analysis of Tooth Autotransplantation

INTRODUCTION

As personalized medicine advances, there is an escalating need for sophisticated tools to understand complex biomechanical phenomena in clinical research. Recognizing a significant gap, this study pioneers the development of patient-specific in silico models for tooth autotransplantation (TAT), setting a new standard for predictive accuracy and reliability in evaluating TAT outcomes.

METHODS

Development of the models relied on 6 consecutive cases of young patients (mean age 11.66 years ± 0.79), all undergoing TAT procedures. The development process involved creating detailed in silico replicas of patient oral structures, focusing on transplanting upper premolars to central incisors. These models underpinned finite element analysis simulations, testing various masticatory and traumatic scenarios.

RESULTS

The models highlighted critical biomechanical insights. The finite element models indicated homogeneous stress distribution in control teeth, contrasted by shape-dependent stress patterns in transplanted teeth. The surface deviation in the postoperative year for the transplanted elements showed a mean deviation of 0.33 mm (±0.28), significantly higher than their contralateral counterparts at 0.05 mm (±0.04).

CONCLUSIONS

By developing advanced patient-specific in silico models, we are ushering in a transformative era in TAT research and practice. These models are not just analytical tools; they are predictive instruments capturing patient uniqueness, including anatomical, masticatory, and tissue variables, essential for understanding biomechanical responses in TAT. This foundational work paves the way for future studies, where applying these models to larger cohorts will further validate their predictive capabilities and influence on TAT success parameters.

A genetic algorithm optimization framework for the characterization of hyper-viscoelastic materials: application to human articular cartilage

This study aims to develop an automated framework for the characterization of materials which are both hyper-elastic and viscoelastic. This has been evaluated using human articular cartilage (AC). AC (26 tissue samples from 5 femoral heads) underwent dynamic mechanical analysis with a frequency sweep from 1 to 90 Hz. The conversion from a frequency- to time-domain hyper-viscoelastic material model was approximated using a modular framework design where finite element analysis was automated, and a genetic algorithm and interior point technique were employed to solve and optimize the material approximations. Three orders of approximation for the Prony series were evaluated at N = 1, 3 and 5 for 20 and 50 iterations of a genetic cycle. This was repeated for 30 simulations of six combinations of the above all with randomly generated initialization points. There was a difference between N = 1 and N = 3/5 of approximately ~5% in terms of the error estimated. During unloading the opposite was seen with a 10% error difference between N = 5 and 1. A reduction of ~1% parameter error was found when the number of generations increased from 20 to 50. In conclusion, the framework has proved effective in characterizing human AC.

Automated 2D and 3D finite element overclosure adjustment and mesh morphing using generalized regression neural networks

Computer representations of three-dimensional (3D) geometries are crucial for simulating systems and processes in engineering and science. In medicine, and more specifically, biomechanics and orthopaedics, obtaining and using 3D geometries is critical to many workflows. However, while many tools exist to obtain 3D geometries of organic structures, little has been done to make them usable for their intended medical purposes. Furthermore, many of the proposed tools are proprietary, limiting their use. This work introduces two novel algorithms based on Generalized Regression Neural Networks (GRNN) and 4 processes to perform mesh morphing and overclosure adjustment. These algorithms were implemented, and test cases were used to validate them against existing algorithms to demonstrate improved performance. The resulting algorithms demonstrate improvements to existing techniques based on Radial Basis Function (RBF) networks by converting to GRNN-based implementations. Implementations in MATLAB of these algorithms and the source code are publicly available at the following locations: https://github.com/thor-andreassen/femors; https://simtk.org/projects/femors-rbf; https://www.mathworks.com/matlabcentral/fileexchange/120353-finite-element-morphing-overclosure-reduction-and-slicing.

Computational models and their applications in biomechanical analysis of mandibular reconstruction surgery

Advanced head and neck cancers involving the mandible often require surgical removal of the diseased parts and replacement with donor bone or prosthesis to recreate the form and function of the premorbid mandible. The degree to which this reconstruction successfully replicates key geometric features of the original bone critically affects the cosmetic and functional outcomes of speaking, chewing, and breathing. With advancements in computational power, biomechanical modeling has emerged as a prevalent tool for predicting the functional outcomes of the masticatory system and evaluating the effectiveness of reconstruction procedures in patients undergoing mandibular reconstruction surgery. These models offer cost-effective and patient-specific treatment tailored to the needs of individuals. To underscore the significance of biomechanical modeling, we conducted a review of 66 studies that utilized computational models in the biomechanical analysis of mandibular reconstruction surgery. The majority of these studies employed finite element method (FEM) in their approach; therefore, a detailed investigation of FEM has also been provided. Additionally, we categorized these studies based on the main components analyzed, including bone flaps, plates/screws, and prostheses, as well as their design and material composition.

Finite element analysis in the optimization of posterior-stabilized total knee arthroplasty

Posterior-stabilized total knee arthroplasty (PS-TKA) is associated with high rates of satisfaction and functional recovery. This is notably attributed to implant optimization in terms of design, choice of materials, positioning and understanding of biomechanics. Finite elements analysis (FEA) is an assessment technique that contributed to this optimization by ensuring mechanical results based on numerical simulation. By close teamwork between surgeons, researchers and engineers, FEA enabled testing of certain clinical impressions. However, the methodological features of the technique led to wide variations in the presentation and interpretation of results, requiring a certain understanding of numerical and biomechanical fields by the orthopedic community. The present study provides an up-to-date review, aiming to address the following questions: what are the principles of FEA? What is the role of FEA in studying PS design in TKA? What are the key elements in the literature for understanding the role of FEA in PS-TKA? What is the contribution of FEA for understanding of tibiofemoral and patellofemoral biomechanical behavior? What are the limitations and perspectives of digital simulation and FEA in routine practice, with a particular emphasis on the "digital twin" concept? LEVEL OF EVIDENCE: V, expert opinion.

A systematic comparison between FEBio and PolyFEM for biomechanical systems

BACKGROUND AND OBJECTIVES

Finite element simulations are widely employed as a non-invasive and cost-effective approach for predicting outcomes in biomechanical simulations. However, traditional finite element software, primarily designed for engineering materials, often encountered limitations in contact detection and enforcement, leading to simulation failure when dealing with complex biomechanical configurations. Currently, a lot of model tuning is required to get physically accurate finite element simulations without failures. This adds significant human interaction to each iteration of a biomechanical model. This study addressed these issues by introducing PolyFEM, a novel finite element solver that guarantees inversion- and intersection-free solutions with completely automatic collision detection. The objective of this research is to validate PolyFEM's capabilities by comparing its results with those obtained from a well-established finite element solver, FEBio.

METHODS

To achieve this goal, five comparison scenarios were formulated to assess and validate PolyFEM's performance. The simulations were reproduced using both PolyFEM and FEBio, and the final results were compared. The five comparison scenarios included: (1) reproducing simulations from the FEBio test suite, consisting of static, dynamic, and contact-driven simulations; (2) replicating simulations from the verification paper published alongside the original release of FEBio; (3) a biomechanically based contact problem; (4) creating a custom simulation involving high-energy collisions between soft materials to highlight the difference in collision methods between the two solvers; and (5) performing biomechanical simulations of biting and quasi-stance.

RESULTS

We found that PolyFEM was capable of replicating all simulations previously conducted in FEBio. Particularly noteworthy is PolyFEM's superiority in high-energy contact simulations, where FEBio fell short, unable to complete over half of the simulations in Scenario 4. Although some of the simulations required significantly more simulation time in PolyFEM compared to FEBio, it is important to highlight that PolyFEM achieved these results without the need for any additional model tuning or contact declaration.

DISCUSSION

Despite being in the early stages of development, PolyFEM currently provides verified solutions for hyperelastic materials that are consistent with FEBio, both in previously published workflows and novel finite element scenarios. PolyFEM exhibited the ability to tackle challenging biomechanical problems where other solvers fell short, thus offering the potential to enhance the accuracy and realism of future finite element analyses.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Assessing Model Calibration Workflows and Outcomes

Model reproducibility is a point of emphasis for the National Institutes of Health (NIH) and in science, broadly. As the use of computational modeling in biomechanics and orthopedics grows, so does the need to assess the reproducibility of modeling workflows and simulation predictions. The long-term goal of the KneeHub project is to understand the influence of potentially subjective decisions, thus the modeler's "art", on the reproducibility and predictive uncertainty of computational knee joint models. In this paper, we report on the model calibration phase of this project, during which five teams calibrated computational knee joint models of the same specimens from the same specimen-specific joint mechanics dataset. We investigated model calibration approaches and decisions, and compared calibration workflows and model outcomes among the teams. The selection of the calibration targets used in the calibration workflow differed greatly between the teams and was influenced by modeling decisions related to the representation of structures, and considerations for computational cost and implementation of optimization. While calibration improved model performance, differences in the postcalibration ligament properties and predicted kinematics were quantified and discussed in the context of modeling decisions. Even for teams with demonstrated expertise, model calibration is difficult to foresee and plan in detail, and the results of this study underscore the importance of identification and standardization of best practices for data sharing and calibration.

An analytical model of lateral condylar plate working length

BACKGROUND

The locking plate is a common device to treat distal femur fractures. Healing is affected by construct stiffness, thus many surgeon-controlled variables such as working length have been examined for their effects on strain at the fracture. No convenient analytical model which aids surgeons in determining working length has yet been described. We propose an analytical model and compare it to finite element analysis and cadaveric biomechanical testing.

METHODS

First, an analytical model based on a cantilever beam equation was derived. Next, a finite element model was developed based on a CT scan of a "fresh-frozen" cadaveric femur. Third, biomechanical testing in single-leg stance loading was performed on the cadaver. In all methods, strain at the fracture was recorded. An ANCOVA test was conducted to compare the strains.

FINDINGS

In all models, as the working length increased so did strain. For strain at the fracture, the shortest working length (35 mm) had a strain of 8% in the analytical model, 9% in the finite element model, and 7% for the cadaver. The longest working length (140 mm) demonstrated strain of 15% in the analytical model, and the finite element and biomechanical tests both demonstrated strain of 14%.

INTERPRETATION

The strain predicted by the analytical model was consistent with the strain observed in both the finite element and biomechanical models. As demonstrated in existing literature, increasing the working length increases strain at the fracture site. Additional work is required to refine and establish validity and reliability of the analytical model.

Validation and evaluation of subject-specific finite element models of the pediatric knee

Finite element (FE) models have been widely used to investigate knee joint biomechanics. Most of these models have been developed to study adult knees, neglecting pediatric populations. In this study, an atlas-based approach was employed to develop subject-specific FE models of the knee for eight typically developing pediatric individuals. Initially, validation simulations were performed at four passive tibiofemoral joint (TFJ) flexion angles, and the resulting TFJ and patellofemoral joint (PFJ) kinematics were compared to corresponding patient-matched measurements derived from magnetic resonance imaging (MRI). A neuromusculoskeletal-(NMSK)-FE pipeline was then used to simulate knee biomechanics during stance phase of walking gait for each participant to evaluate model simulation of a common motor task. Validation simulations demonstrated minimal error and strong correlations between FE-predicted and MRI-measured TFJ and PFJ kinematics (ensemble average of root mean square errors < 5 mm for translations and < 4.1° for rotations). The FE-predicted kinematics were strongly correlated with published reports (ensemble average of Pearson's correlation coefficients (ρ) > 0.9 for translations and ρ > 0.8 for rotations), except for TFJ mediolateral translation and abduction/adduction rotation. For walking gait, NMSK-FE model-predicted knee kinematics, contact areas, and contact pressures were consistent with experimental reports from literature. The strong agreement between model predictions and experimental reports underscores the capability of sequentially linked NMSK-FE models to accurately predict pediatric knee kinematics, as well as complex contact pressure distributions across the TFJ articulations. These models hold promise as effective tools for parametric analyses, population-based clinical studies, and enhancing our understanding of various pediatric knee injury mechanisms. They also support intervention design and prediction of surgical outcomes in pediatric populations.

The relative influence of model parameters on finite element analysis simulations of intervertebral body fusion device static compression performance

Non-clinical mechanical performance testing is a critical aspect of intervertebral body fusion device (IBFD) development and regulatory evaluation. Recently, stakeholders have begun leveraging computational modeling and simulations such as finite element analysis (FEA) in addition to traditional bench testing. FEA offers advantages such as reduced experiment time, lower costs associated with elimination of bench testing (e.g. specimen manufacture and test execution), and elucidating quantities of interest that traditional testing cannot provide (e.g. stress and strain distributions). However, best practices for FEA of IBFDs are not well defined, and modeler decision making can significantly influence simulation setup and results. Therefore, the goal of this study was to determine the relative influence of modeling parameters when using FEA to assess non-clinical mechanical performance of IBFDs. FEA was used to conduct a series of IBFD static uniaxial compression simulations. Several parameters relating to implant geometry, loading/boundary conditions, and material properties were carefully controlled to assess their relative influence on two output variables (IBFD stiffness and yield load). Results were most influenced by device geometry, while the effects of boundary conditions and material properties were more significant within IBFDs of identical or similar geometries. These results will aid stakeholders in the development of standardized best practices for using FEA to assess non-clinical mechanical performance of IBFDs.

Validation of a novel finite-element model for evaluating patellofemoral forces and stress during squatting after posterior-stabilized total knee arthroplasty

INTRODUCTION

Several studies have documented the relationship between patellofemoral pain and patient dissatisfaction after total knee arthroplasty (TKA). However, few computer simulations have been designed to evaluate the patellofemoral joint during flexion. The aim of this study was to validate a new computational simulation, driven by forces and moments, and to analyze patellofemoral reaction forces and stress under squat loading conditions after TKA implantation.

HYPOTHESIS

This computational simulation of a squat using a model driven by forces and moments is comparable to in vitro and in silico data from the literature.

MATERIAL AND METHODS

We developed a finite element model of the lower limb after implantation of a fixed-bearing posterior-stabilized TKA. To simulate squat loading conditions when standing on both legs, an initial load of 130N was applied to the center of the femoral head. Quadriceps force, patellofemoral contact force and Von Mises stress on the patellar implant, tibiofemoral contact forces and pressure on the tibial insert, and post-cam contact force were evaluated from 0° to 100° of knee flexion.

RESULTS

Quadriceps force increased during flexion, up to 6 times the applied load. Von Mises stress on patellar implant increased up to 16MPa at 100° flexion. Tibiofemoral contact forces increased up to 415 N medially and 339 N laterally, with 64% distributed medially on the tibial insert. Post-cam contact started slightly before 70° of flexion.

DISCUSSION

In this simulation, tibiofemoral, patellofemoral and post-cam contact forces, and pressure distribution on the tibial insert were consistent with various published studies. This agreement suggests that computational simulation driven by forces and moments can reproduce squat loading conditions during knee flexion after TKA, without experimental kinematic data used to drive the simulation.

CONCLUSION

This study represents an initial step towards validating tibiofemoral and patellofemoral mechanical behavior under squat conditions, from this computational simulation driven by forces and moments. This model will help us better understand the influence of various implantation techniques on patellofemoral forces and stress during flexion.

LEVEL OF EVIDENCE

IV, biomechanical computational study.

Elucidation of the Mechanism of Occasional Anterior Longitudinal Ligament Rupture with Posterior Correction Procedure for Adult Spinal Deformity Using LLIF-Finite Element Analysis of the Impact of the Lordotic Angle of Intervertebral LLIF Cage

Background and Objectives: There are several advantages of using lateral lumbar interbody fusion (LLIF) for correction surgeries for adult spinal deformity (ASD); however, we currently have unresolved new issues, including occasional anterior longitudinal ligament (ALL) rupture during the posterior correction procedure. When LLIF was initially introduced, only less lordotic cages were available and ALL rupture was more frequently experienced compared with later periods when more lordotic cages were available. We performed finite element analysis (FEA) regarding the mechanism of ALL rupture during a posterior correction procedure. Methods: A spring (which mimics ALL) was introduced at the location of ALL in the FEA and an LLIF cage with two different lordotic angles, 6 and 12 degrees (6DC/12DC), was employed. To assess the extent of burden on the ALL, the extension length of the spring during the correction procedure was measured and the location of the rotation center was examined. Results: We observed a significantly higher degree of length extension of the spring during the correction procedure in the FEA model with 6DC compared with that of 12DC. We also observed that the location of the rotation center was shifted posteriorly in the FEA model with 6DC compared with that of 12DC. Conclusions: It is considered that the posterior and rostral edge of the less lordotic angle cage became a hinge, and the longer lever arm increased the burden on ALL as the principle of leverage. It is important to use an LLIF cage with a sufficient lordotic angle, that is compatible with the degree of posterior osteotomy in ASD correction.

Advanced Mechanical Testing Technologies at the Cellular Level: The Mechanisms and Application in Tissue Engineering

Mechanics, as a key physical factor which affects cell function and tissue regeneration, is attracting the attention of researchers in the fields of biomaterials, biomechanics, and tissue engineering. The macroscopic mechanical properties of tissue engineering scaffolds have been studied and optimized based on different applications. However, the mechanical properties of the overall scaffold materials are not enough to reveal the mechanical mechanism of the cell-matrix interaction. Hence, the mechanical detection of cell mechanics and cellular-scale microenvironments has become crucial for unraveling the mechanisms which underly cell activities and which are affected by physical factors. This review mainly focuses on the advanced technologies and applications of cell-scale mechanical detection. It summarizes the techniques used in micromechanical performance analysis, including atomic force microscope (AFM), optical tweezer (OT), magnetic tweezer (MT), and traction force microscope (TFM), and analyzes their testing mechanisms. In addition, the application of mechanical testing techniques to cell mechanics and tissue engineering scaffolds, such as hydrogels and porous scaffolds, is summarized and discussed. Finally, it highlights the challenges and prospects of this field. This review is believed to provide valuable insights into micromechanics in tissue engineering.

Anatomical Tissue Engineering of the Anterior Cruciate Ligament Entheses

The firm integration of anterior cruciate ligament (ACL) grafts into bones remains the most demanding challenge in ACL reconstruction, since graft loosening means graft failure. For a functional-tissue-engineered ACL substitute to be realized in future, robust bone attachment sites (entheses) have to be re-established. The latter comprise four tissue compartments (ligament, non-calcified and calcified fibrocartilage, separated by the tidemark, bone) forming a histological and biomechanical gradient at the attachment interface between the ACL and bone. The ACL enthesis is surrounded by the synovium and exposed to the intra-articular micromilieu. This review will picture and explain the peculiarities of these synovioentheseal complexes at the femoral and tibial attachment sites based on published data. Using this, emerging tissue engineering (TE) strategies addressing them will be discussed. Several material composites (e.g., polycaprolactone and silk fibroin) and manufacturing techniques (e.g., three-dimensional-/bio-printing, electrospinning, braiding and embroidering) have been applied to create zonal cell carriers (bi- or triphasic scaffolds) mimicking the ACL enthesis tissue gradients with appropriate topological parameters for zones. Functionalized or bioactive materials (e.g., collagen, tricalcium phosphate, hydroxyapatite and bioactive glass (BG)) or growth factors (e.g., bone morphogenetic proteins [BMP]-2) have been integrated to achieve the zone-dependent differentiation of precursor cells. However, the ACL entheses comprise individual (loading history) asymmetric and polar histoarchitectures. They result from the unique biomechanical microenvironment of overlapping tensile, compressive and shear forces involved in enthesis formation, maturation and maintenance. This review should provide a road map of key parameters to be considered in future in ACL interface TE approaches.

Positive medial cortical support versus anatomical reduction for trochanteric hip fractures: Finite element analysis and biomechanical testing

BACKGROUND AND OBJECTIVES

The anatomical reduction (AR) is usually considered the best option for fractures. Nevertheless, in unstable trochanteric hip fractures (UTHF), previous clinical reports found that the positive medial cortical support (PMCS, an over-reduction technique) attained higher mechanical stability, but this challenging clinical finding still needs experimental validation.

METHODS

This study constructed in-silico and biomechanical PMCS and AR models, with the use of the most clinically-representative geometry design of fracture models, the multi-directional design in FE analysis, and the subject-specific (osteoporotic) bone material properties, to make the models better mimic the actual condition in clinical settings. Then multiple performance variables (von-Mises stress, strain, integral axial stiffness, displacement, structural changes, etc.) were assessed to uncover details of integral and regional stability.

RESULTS

Among in-silico comparison, PMCS models showed significantly lower maximum displacement than AR models, and the maximum von Mises stress of implants (MVMS-I) was significantly lower in PMCS models than in AR models (highest MVMS-I in -30°-A3-AR of 1055.80 ± 93.37 MPa). Besides, PMCS models had significantly lower maximum von Mises stress along fracture surfaces (MVMS-F) (highest MVMS-F in 30°-A2-AR of 416.40 ± 38.01 MPa). Among biomechanical testing comparison, PMCS models showed significantly lower axial displacement. Significantly lower change of neck-shaft angle (CNSA) was observed in A2-PMCS models. A fair amount of AR models converted into the obvious negative medial cortical support (NMCS) condition, whereas all PMCS models kept the PMCS condition. The results were also validated through comparison to previous clinical data.

CONCLUSIONS

The PMCS is superior to the AR in the UTHF surgery. The current study opens up the second thought of the role of over-reduction technique in bone surgery.

Clinically useful finite element models of the natural ankle - A review

BACKGROUND

Biomechanical simulation of the foot and ankle complex is a growing research area but compared to simulation of joints such as hip and knee, it has been under investigated and lacks consistency in research methodology. The methodology is variable, data is heterogenous and there are no clear output criteria. Therefore, it is very difficult to correlate clinically and draw meaningful inferences.

METHODS

The focus of this review is finite element simulation of the native ankle joint and we will explore: the different research questions asked, the model designs used, ways the model rigour has been ensured, the different output parameters of interest and the clinical impact and relevance of these studies.

FINDINGS

The 72 published studies explored in this review demonstrate wide variability in approach. Many studies demonstrated a preference for simplicity when representing different tissues, with the majority using linear isotropic material properties to represent the bone, cartilage and ligaments; this allows the models to be complex in another way such as to include more bones or complex loading. Most studies were validated against experimental or in vivo data, but a large proportion (40%) of studies were not validated at all, which is an area of concern.

INTERPRETATION

Finite element simulation of the ankle shows promise as a clinical tool for improving outcomes. Standardisation of model creation and standardisation of reporting would increase trust, and enable independent validation, through which successful clinical application of the research could be realised.

A Validated Three-Dimensional, Heterogenous Finite Element Model of the Rotator Cuff and The Effects of Collagen Orientation

Continuum mechanics-based finite element models of the shoulder aim to quantify the mechanical environment of the joint to aid in clinical decision-making for rotator cuff injury and disease. These models allow for the evaluation of the internal loading of the shoulder, which cannot be measured in-vivo. This study uses human cadaveric rotator cuff samples with surface tendon strain estimates, to validate a heterogeneous finite element model of the supraspinatus-infraspinatus complex during various load configurations. The computational model was considered validated when the absolute difference in average maximum principal strain for the articular and bursal sides for each load condition estimated by the model was no greater than 3% compared to that measured in the biomechanical study. The model can predict the strains for varying infraspinatus loads allowing for the study of load sharing between these two tightly coordinated tendons. The future goal is to use the modularity of this validated model to study the initiation and propagation of rotator cuff tear and other rotator cuff pathologies to ultimately improve care for rotator cuff tear patients.

Comparison of 3 different bone-borne type expansion appliances used in surgically-assisted rapid palatal expansion: A finite element analysis

INTRODUCTION

This study aimed to compare the effects of 3 different bone-borne type expansion appliances used in the surgically-assisted rapid palatal expansion (SARPE) by finite element analysis.

METHODS

Three different miniscrew-supported palatal expansion appliances were modeled. Median and lateral osteotomies were performed without pterygomaxillary suture separation. Model I consisted of a palatal expander with 2 miniscrews placed 4 mm far from the midpalatal suture. In model II, 2 miniscrews were located at the alveolar ridge between the first molar and the second premolar. In model III, 4 miniscrews were placed as a combination of the first and second models. Stress distributions and amount of displacements were evaluated with Ansys software (version 19.2; Ansys, Canonsburg, Pa) for 5-mm expansion in a symmetrical finite element analysis model to reflect the clinical situation.

RESULTS

SARPE simulation using miniscrew-assisted maxillary expanders for all models showed a rotation and tipping of the maxilla. The largest displacement was found for the anterior part of the palate in model II and the posterior part in model III. Although a wedge-shaped expansion pattern was observed in all models, this form was more prominent in model II. The highest stress value (0.91 MPa) was measured in model I, and the lowest value (0.004 MPa) was measured in model II for the anterior nasal spine region. The highest stress value (0.51 MPa) was measured in model III, and the lowest value (0.12 MPa) was measured in model II for the posterior nasal spine region. The lowest stress values were measured in model II for all the craniofacial and maxillofacial structures.

CONCLUSIONS

Among the models, the lowest stress distribution conditions for craniofacial and maxillofacial structures were found in model II. The largest displacement was found at the incisors and anterior part of the maxilla for model II. The greatest displacement was found at the posterior region for model III.

A review of the characterizations of soft tissues used in human body modeling: Scope, limitations, and the path forward

Soft tissue material properties are vital to human body models that evaluate interactions between the human body and its environment. Such models evaluate internal stress/strain responses in soft tissues to investigate issues like pressure injuries. Numerous constitutive models and parameters have been used to represent mechanical behavior of soft tissues in biomechanical models under quasi-static loading. However, researchers reported that generic material properties cannot accurately represent specific target populations due to large inter-individual variability. Two challenges that exist are experimental mechanical characterization and constitutive modeling of biological soft tissues and personalization of constitutive parameters using non-invasive, non-destructive bedside testing methods. It is imperative to understand the scope and appropriate applications for reported material properties. Thus, the goal of this paper was to compile studies from which soft tissue material properties were obtained and categorize them by source of tissue samples, methods used to quantify deformation, and material models used to describe tissues. The collected studies displayed wide ranges of material properties, and factors that affected the properties included whether tissue samples were in vivo or ex vivo, from humans or animals, the body region tested, body position during in vivo studies, deformation measurements, and material models used to describe tissues. Because of the factors that affected reported material properties, it is clear that much progress has been made in understanding soft tissue responses to loading, yet there is a need to broaden the scope of reported soft tissue material properties and better match reported properties to appropriate human body models.

Assessment of reporting practices and reproducibility potential of a cohort of published studies in computational knee biomechanics

Reproducible research serves as a pillar of the scientific method and is a foundation for scientific advancement. However, estimates for irreproducibility of preclinical science range from 75% to 90%. The importance of reproducible science has not been assessed in the context of mechanics-based modeling of human joints such as the knee, despite this being an area that has seen dramatic growth. Framed in the context of five experienced teams currently documenting knee modeling procedures, the aim of this study was to evaluate reporting and the perceived potential for reproducibility across studies the teams viewed as important contributions to the literature. A cohort of studies was selected by polling, which resulted in an assessment of nine studies as opposed to a broader analysis across the literature. Using a published checklist for reporting of modeling features, the cohort was evaluated for both "reporting" and their potential to be "reproduced," which was delineated into six major modeling categories and three subcategories. Logistic regression analysis revealed that for individual modeling categories, the proportion of "reported" occurrences ranged from 0.31, 95% confidence interval (CI) [0.23, 0.41] to 0.77, 95% CI: [0.68, 0.86]. The proportion of whether a category was perceived as "reproducible" ranged from 0.22, 95% CI: [0.15, 0.31] to 0.44, 95% CI: [0.35, 0.55]. The relatively low ratios highlight an opportunity to improve reporting and reproducibility of knee modeling studies. Ongoing efforts, including our findings, contribute to a dialogue that facilitates adoption of practices that provide both credibility and translation possibilities.

A Review of Image-Based Simulation Applications in High-Value Manufacturing

Image-Based Simulation (IBSim) is the process by which a digital representation of a real geometry is generated from image data for the purpose of performing a simulation with greater accuracy than with idealised Computer Aided Design (CAD) based simulations. Whilst IBSim originates in the biomedical field, the wider adoption of imaging for non-destructive testing and evaluation (NDT/NDE) within the High-Value Manufacturing (HVM) sector has allowed wider use of IBSim in recent years. IBSim is invaluable in scenarios where there exists a non-negligible variation between the 'as designed' and 'as manufactured' state of parts. It has also been used for characterisation of geometries too complex to accurately draw with CAD. IBSim simulations are unique to the geometry being imaged, therefore it is possible to perform part-specific virtual testing within batches of manufactured parts. This novel review presents the applications of IBSim within HVM, whereby HVM is the value provided by a manufactured part (or conversely the potential cost should the part fail) rather than the actual cost of manufacturing the part itself. Examples include fibre and aggregate composite materials, additive manufacturing, foams, and interface bonding such as welding. This review is divided into the following sections: Material Characterisation; Characterisation of Manufacturing Techniques; Impact of Deviations from Idealised Design Geometry on Product Design and Performance; Customisation and Personalisation of Products; IBSim in Biomimicry. Finally, conclusions are drawn, and observations made on future trends based on the current state of the literature.

Enhanced bone formation in locally-optimised, low-stiffness additive manufactured titanium implants: An in silico and in vivo tibial advancement study

A tibial tuberosity advancement (TTA), used to treat lameness in the canine stifle, provides a framework to investigate implant performance within an uneven loading environment due to the dominating patellar tendon. The purpose of this study was to reassess how we design orthopaedic implants in a load-bearing model to investigate potential for improved osseointegration capacity of fully-scaffolded mechanically-matched additive manufactured (AM) implants. While the mechanobiological nature of bone is well known, we have identified a lower limit in the literature where investigation into exceedingly soft scaffolds relative to trabecular bone ceases due to the trade-off in mechanical strength. We developed a finite element model of the sheep stifle to assess the stresses and strains of homogeneous and locally-optimised TTA implant designs. Using additive manufacturing, we printed three different low-stiffness Ti-6Al-4 V TTA implants: 0.8 GPa (Ti1), 0.6 GPa (Ti2) and an optimised design with a 0.3 GPa cortex and 0.1 GPa centre (Ti3), for implantation in a 12-week in vivo ovine pilot study. Static histomorphometry demonstrated uniform bone ingrowth in optimised low-modulus Ti3 samples compared to homogeneous designs (Ti1 and Ti2), and greater bone-implant contact. Mineralising surfaces were apparent in all implants, though mineral apposition rate was only consistent throughout Ti3. The greatest bone formation scores were seen in Ti3, followed by Ti2 and Ti1. Results from our study suggest lower stiffnesses and higher strain ranges improve early bone formation, and that by accounting for loading environments through rational design, implants can be optimised to improve uniform osseointegration. STATEMENT OF SIGNIFICANCE: The effect of different strain ranges on bone healing has been traditionally investigated and characterised through computational models, with much of the literature suggesting higher strain ranges being favourable. However, little has been done to incorporate strain-optimisation into porous orthopaedic implants due to the trade-off in mechanical strength required to induce these microenvironments. In this study, we used finite element analysis to optimise the design of additive manufactured (AM) titanium orthopaedic implants for different strain ranges, using a clinically-relevant surgical model. Our research suggests that there is potential for locally-optimised AM scaffolds in the use of orthopaedic devices to induce higher strains, which in turn encourages de novo bone formation and uniform osseointegration.

Open Knee(s): A Free and Open Source Library of Specimen-Specific Models and Related Digital Assets for Finite Element Analysis of the Knee Joint

There is a growing interest in the use of virtual representations of the knee for musculoskeletal research and clinical decision making, and to generate digital evidence for design and regulation of implants. Accessibility to previously developed models and related digital assets can dramatically reduce barriers to entry to conduct simulation-based studies of the knee joint and therefore help accelerate scientific discovery and clinical innovations. Development of models for finite element analysis is a demanding process that is both time consuming and resource intensive. It necessitates expertise to transform raw data to reliable virtual representations. Modeling and simulation workflow has many processes such as image segmentation, surface geometry generation, mesh generation and finally, creation of a finite element representation with relevant loading and boundary conditions. The outcome of the workflow is not only the end-point knee model but also many other digital by-products. When all of these data, derivate assets, and tools are freely and openly accessible, researchers can bypass some or all the steps required to build models and focus on using them to address their research goals. With provenance to specimen-specific anatomical and mechanical data and traceability of digital assets throughout the whole lifecycle of the model, reproducibility and credibility of the modeling practice can be established. The objective of this study is to disseminate Open Knee(s), a cohort of eight knee models (and relevant digital assets) for finite element analysis, that are based on comprehensive specimen-specific imaging data. In addition, the models and by-products of modeling workflows are described along with model development strategies and tools. Passive flexion served as a test simulation case, demonstrating an end-user application. Potential roadmaps for reuse of Open Knee(s) are also discussed.

The comparison of stress and strain between custom-designed bone plates (CDBP) and locking compression plate (LCP) for distal femur fracture

BACKGROUND

Distal femur fracture is considered one of the most common fractures due to high-energy traumas such as car accidents or low-energy traumas such as osteoporosis. Locking plates are orthopedic implants used for stabilized femur fracture. Thus, designing a bone plate fitted exactly with the patient's bone and correctly fixing bone segments are required for better fracture healing.

OBJECTIVES

This study aims to design a bone plate based on anthropometric characteristics of patients' femurs and compare performing custom-designed bone plates (CDBP) with the locking compression plate (LCP) by finite element method.

MATERIALS AND METHODS

In this analytical study, a 3D model of four patients' femur and CDBP were firstly designed in MIMICS 19.0 based on the patient's femur anatomy. After designing the bone plate, the CDBPs and LCP were fixed on the bone and analyzed by finite element method (FEM) in ANSYS, and stress and strain of bone plates were also compared.

RESULTS

The maximum principal stress for all 3D models of patients' fracture femur by CDBPs was stabilized better than LCP with a decrease by 39.79, 12.54, 9.49, and 20.29% in 4 models, respectively. Also, in all models, the strain of CDBPs is less than LCP. Among the different thicknesses considered, the bone plate with 5 mm thickness showed better stress and strain distribution than other thicknesses.

CONCLUSION

Customized bone plate designed based on patient's femur anatomical morphology shows better bone-matching plate, resulting in increasing the quality of the fracture healing and fails to any need for additional shaping.

TRIAL REGISTRATION NUMBER

Design and analysis of an implant were investigated in this study. There was no intervention in the diagnosis and treatment of patients and the study was not a clinical trial.

Effect of screw tunnels on proximal femur strength after screw removal: A finite element analysis

BACKGROUND

The presence of screw tunnels in the femoral neck is a problem for patients with proximal femoral fractures after removal of internal fixation. The question of how much does the existence of the screw tunnels affect the strength of the femur and whether the patient needs to be protected with an adjunctive device has been controversial. The objective of this finite element analysis was to determine (1) whether the screw tunnels affects normal weight bearing after removal of internal fixation of a proximal femur fracture, (2) which screw tunnels parameters affect the weight bearing capacity of the entire femur.

HYPOTHESIS

The presence of the screw tunnels reduces the load-bearing capacity of the femur, and the arrangement, diameter and wall thickness of the screw tunnels affect the load-bearing capacity of the femur.

MATERIALS AND METHODS

Twenty patients who underwent surgical treatment for proximal femur fracture at our hospital were included in the study. Computed tomography (CT) values of the screw tunnel wall in the femur after removal of internal fixations were analysed. Mimics v16.0 and Hypermesh v13.0 software programs were used to generate 3-dimensional (3D) tetrahedral finite element models of the proximal femur with different screw tunnel numbers, diameters, thicknesses, and arrangements. An acetabulum exerting a vertical pressure load of 600N on the femoral head was simulated and the force on various parts of the femur in each model was calculated.

RESULTS

There was no difference in the Hounsfield Units of the tunnel walls and cortical bone of the proximal femur (893.48±61.28 vs. 926.34±58.43; p=0.091). In each of the 3D models, the cancellous bone was the first structure to reach maximal stress. The compressive strength of the femur decreased with increasing thickness of the screw tunnel wall and decreased with increasing tunnel diameter. The femoral neck model with the inverted triangle screw tunnel arrangement had the highest compressive strength.

DISCUSSION

The femoral neck with screw tunnels can withstand day-to-day stress without special intervention. For femoral neck fractures fixed with cannulated screws, inverted triangle screws are recommended; For a single screw tunnel in the femoral neck, the larger the diameter of the femoral neck internal screw channel, the weaker the load-bearing capacity of the femur.

LEVEL OF EVIDENCE

III; well-designed computational non-experimental study.

Personalizing Revision Tibial Baseplate Position and Stem Trajectory With Custom Implants Using 3D Modeling to Optimize Press-fit Stem Placement

Background

A common tibial construct for revision total knee arthroplasty includes a long diaphyseal engaging press-fit stem. Due to tibial canal bowing, compromises are often necessary to match patient anatomy when choosing stemmed implants. The objective of this study is to determine through 3-D modeling whether current implant press-fit options appropriately fit patient anatomy, or whether an alternative angle between the stem and baseplate could increase the cortical engagement of long press-fit tibial stems.

Methods

Preoperative computerized tomography scans from 100 patients undergoing TKA were imported into an image-processing software program. Three-dimensional models were created with tibial stems placed at a fixed perpendicular angle and a custom angle to the revision tibial baseplate. Stem diameter, depth, offset, and contact surface area were measured and analyzed between the 2 groups.

Results

Significantly more cortical contact, larger stem diameter, and smaller offset of the custom keel from the center of the baseplate were associated with free custom tibial stem placement vs a fixed perpendicular baseplate-stem interface (P < .001). Statistically significant differences were also found between different patient demographics.

Conclusions

Custom free-angle stem placement allows for increased stem diameter and cortical contact of press-fit tibial stems compared to existing constructs that must interface with the baseplate at a 90-degree angle. Current revision tibia implants limit fixation of tibial press-fit stems and often mismatch with patient anatomy. Alternative ways to fit patient anatomy may be beneficial for patients with extreme mismatch. In the future, custom keel angles may help to resolve this problem.

EXPERIMENTAL SUBSTANTIATION OF THE BIOMECHANICAL EFFICIENCY OF THE ANTERIOR PLATE COMBINED WITH TWO SPONGIOUS SCREWS FOR ANKLE ARTHRODESIS

Актуальность. Несмотря на существующее значительное количество разнообразных методик, применяемых для артродеза голеностопного сустава, ряд авторов указывают на определенные технические сложности указанных операций, потерю коррекции заданного положения таранной и большеберцовой костей в процессе анкилозирования, несращения. Проблема совершенствования методики фиксации при артродезе голеностопного сустава требует новых решений. Цель провести сравнительный биомеханический анализ стабильности фиксации при артродезе голеностопного сустава тремя спонгиозными винтами и разработанной пластиной, комбинируемой с двумя спонгиозными винтами, методом конечных элементов. Материалы и методы. Методом конечных элементов выполнена оценка биомеханических характеристик трех вариантов систем фиксации голеностопного сустава при артродезе: три спонгиозных винта, разработанная пластина, комбинируемая с двумя спонгиозными винтами, проксимальный винт в пластине кортикальный, а также разработанная пластина, комбинируемая с двумя спонгиозными винтами, проксимальный винт в пластине с угловой стабильностью. Изучены напряжения и деформации при приложении различных видов нагрузок. Результаты. В модели фиксации голеностопного сустава передней пластиной, комбинируемой с двумя спонгиозными винтами и проксимальным кортикальным винтом, имплантаты и таранная кость испытывали наименьшие напряжения по сравнению с двумя другими моделями. Так, максимальное эквивалентное напряжение в имплантатах при втором варианте составило 68-124 МПа, при первом варианте 92-147 МПа, при третьем 130-331 МПа. Эквивалентное напряжение в таранной кости во втором варианте фиксации составило от 20 до 46 МПа, в первом и третьем вариантах 28-58 МПа и 47-65 МПа, соответственно. Показатели максимального контактного давления на границе большеберцовой и таранной костей оказались наибольшими в первом варианте по сравнению с двумя другими моделями (34 МПа, 31 МПа и 31 МПа соответственно). Заключение. Среди изученных систем фиксации голеностопного сустава для артродеза, применение передней пластины, комбинируемой с двумя спонгиозными винтами и проксимальным кортикальным винтом, является наиболее предпочтительным с точки зрения биомеханики.

Precision medicine using patient-specific modelling: state of the art and perspectives in dental practice

The dental practice has largely evolved in the last 50 years following a better understanding of the biomechanical behaviour of teeth and its supporting structures, as well as developments in the fields of imaging and biomaterials. However, many patients still encounter treatment failures; this is related to the complex nature of evaluating the biomechanical aspects of each clinical situation due to the numerous patient-specific parameters, such as occlusion and root anatomy. In parallel, the advent of cone beam computed tomography enabled researchers in the field of odontology as well as clinicians to gather and model patient data with sufficient accuracy using image processing and finite element technologies. These developments gave rise to a new precision medicine concept that proposes to individually assess anatomical and biomechanical characteristics and adapt treatment options accordingly. While this approach is already applied in maxillofacial surgery, its implementation in dentistry is still restricted. However, recent advancements in artificial intelligence make it possible to automate several parts of the laborious modelling task, bringing such user-assisted decision-support tools closer to both clinicians and researchers. Therefore, the present narrative review aimed to present and discuss the current literature investigating patient-specific modelling in dentistry, its state-of-the-art applications, and research perspectives.

Mechanism of Staple Line Leak After Sleeve Gastrectomy via Isobaric Pressurisation Concentrating Stress Forces at the Proximal Staple Line

Purpose

Staple line leak following sleeve gastrectomy is a significant problem and has been hypothesised to be related to hyperpressurisation in the proximal stomach. There is, however, little objective evidence demonstrating how these forces could be transmitted to the luminal wall. We aimed to define conditions in the proximal stomach and simulate the transmission of stress forces in the post-operative stomach using a finite element analysis (FEA).

Materials and Methods

The manometry of fourteen patients post sleeve gastrectomy was compared to ten controls. Manometry, boundary conditions, and volumetric CT were integrated to develop six models. These models delineated luminal wall stress in the proximal stomach. Key features were then varied to establish the influence of each factor.

Results

The sleeve gastrectomy cohort had a significantly higher peak intragastric isobaric pressures 31.58 ± 2.1 vs. 13.49 ± 1.3 mmHg (p = 0.0002). Regions of stress were clustered at the staple line near the GOJ, and peak stress was observed there in 67% of models. A uniform greater curvature did not fail or concentrate stress under maximal pressurisation. Geometric variation demonstrated that a larger triangulated apex increased stress by 17% (255 kPa versus 218 kPa), with a 37% increase at the GOJ (203kPA versus 148kPA). A wider incisura reduced stress at the GOJ by 9.9% (128 kPa versus 142 kPa).

Conclusion

High pressure events can occur in the proximal stomach after sleeve gastrectomy. Simulations suggest that these events preferentially concentrate stress forces near the GOJ. This study simulates how high-pressure events could translate stress to the luminal wall and precipitate leak.

Graphical Abstract

Effect of Intramedullary Nailing Patterns on Interfragmentary Strain in a Mouse Femur Fracture: A Parametric Finite Element Analysis

Delayed long bone fracture healing and nonunion continue to be a significant socioeconomic burden. While mechanical stimulation is known to be an important determinant of the bone repair process, understanding how the magnitude, mode, and commencement of interfragmentary strain (IFS) affect fracture healing can guide new therapeutic strategies to prevent delayed healing or nonunion. Mouse models provide a means to investigate the molecular and cellular aspects of fracture repair, yet there is only one commercially available, clinically-relevant, locking intramedullary nail (IMN) currently available for studying long bone fractures in rodents. Having access to alternative IMNs would allow a variety of mechanical environments at the fracture site to be evaluated, and the purpose of this proof-of-concept finite element analysis study is to identify which IMN design parameters have the largest impact on IFS in a murine transverse femoral osteotomy model. Using the dimensions of the clinically relevant IMN as a guide, the nail material, distance between interlocking screws, and clearance between the nail and endosteal surface were varied between simulations. Of these parameters, changing the nail material from stainless steel (SS) to polyetheretherketone (PEEK) had the largest impact on IFS. Reducing the distance between the proximal and distal interlocking screws substantially affected IFS only when nail modulus was low. Therefore, IMNs with low modulus (e.g., PEEK) can be used alongside commercially available SS nails to investigate the effect of initial IFS or stability on fracture healing with respect to different biological conditions of repair in rodents.

Reporting guidelines for in-silico studies using finite element analysis in medicine (RIFEM)

BACKGROUND

To the best of our knowledge, there are no reporting guidelines for design, conduct and reporting of Finite Element studies in health sciences. We intend to propose specific and detailed guidelines for reporting these studies.

METHOD

After recognizing the need to have uniform guidelines for reporting of finite element analysis in medicine and dentistry, a group of 5 researchers working on FEA as their research area met in the summer of 2020 and drafted the methodology for the development of such guidelines. Each researcher individually made a list of major headings required for reporting these studies and met again in September 2020 to finalize the domains. Subsequently, sub headings and details were charted. The draft list of items for reporting the guidelines were presented to a larger team of 15 experts and some changes were further made based on their inputs.

RESULTS

The guidelines entail seven major domains and their sub-domains, including parameters for model structure, segmentation, mesh structure, force application and model validation, etc. This checklist aims to improvise the reporting and consistency of FEA studies.

CONCLUSION

We hope that the usage and adoption of these guidelines by the scientific community would result in more thoughtful and uniform documentation. Also, the confidence in the results would be enhanced through model reproducibility, reusability and accountability. The proposed guidelines were named as 'Reporting of in-silico studies using finite element analysis in medicine' and the term 'RIFEM' was used as acronym.

FINITE ELEMENT SPINE MODELS AND SPINAL INSTRUMENTS: A REVIEW

There is considerable biomechanics literature on finite element modeling and analysis of the spine. To accurately mimic the biomechanical behavior of the vertebral column, a generated computational model has to include anatomical structures that are consistent with physiological reality. In this review article, we focused on the finite element spine models that have been developed by various approaches in the literature. Firstly, the anatomical features of the spine and the spinal components have been briefly explained. We then focused on the modeling stages of vertebrae, ligaments, facet joints, intervertebral discs, and spinal instruments. With this paper, we expect to provide a comprehensive resource regarding the modeling preferences used in spine modeling.

Finite element analysis of stress distribution in autotransplanted molars

OBJECTIVE

The biomechanical response of an autotransplanted tooth and surrounding bone to occlusal loads is not well-known. The aim of the present study was to investigate the effect of root form and occlusal morphology on stress distribution in autotransplanted teeth and surrounding bone by using finite element analysis (FEA).

METHODS

Seven FEA models representing different autotransplanted tooth situations were generated: (a) first molar, (b) third molar, (c) root canal-treated third molar, (d) root canal-treated, ankylosed, third molar, (e) crowned third molar, (f) crowned and root canal-treated third molar, (g) root canal-treated, ankylosed, and crowned third molar. Load (200 N) was applied on the occlusal surface, parallel to the long axis of the tooth. Maximum von Mises stress values on dentin and surrounding bone were calculated for each situation.

RESULTS

Differences in stress distribution were observed among models. In ankylosed model, stress was primarily observed at the coronal region of the tooth. The stress was observed more at the coronal region of the tooth in crowned models compared with the non-crowned models. The stress distribution was homogeneous with root canal-treated and crowned autotransplanted tooth.

CONCLUSIONS

The occlusal morphology and root form of the autotransplanted tooth affected the stress in surrounding bone at the transfer site and the biomechanical response of the tooth. The stress was more homogeneous in crowned tooth and primarily observed at the coronal region, which may decrease the risk for root resorption.

CLINICAL SIGNIFICANCE

Root configuration, occlusal form and root canal treatment induce significant changes on the stress distribution on teeth and bone, including characteristic stress concentration and increased stress values. Clinicians can consider crowning autotransplanted teeth for improved stress distribution within the tooth structure.

Sealing mechanism study of laryngeal mask airways via 3D modelling and finite element analysis

Proper sealing of laryngeal mask airways (LMAs) is critical for airway management in clinical use. Understanding the sealing mechanism can significantly help front-line anaesthetists to reduce the incidence of adverse events. However, anaesthetists, who may not have the most substantial engineering backgrounds, lack intuitive ways to develop an understanding of the LMA sealing mechanism effectively. The paper aims to study the LMA-pharynx sealing mechanisms from the perspective of front-line anaesthetists. We use a computer-aided 3D modelling technique to visualise the LMA—pharynx interactions, which helps anaesthetists identify the critical areas of complications. Furthermore, we conduct a quantitative pressure distribution analysis of the LMA-pharynx contacting surface using the finite element analysis technique, which helps further understand the sealing mechanics in those areas. We present two cases studies based on one male volunteer, aged 50, inserted with a ProSeal LMA. In the first case, a relatively low cuff pressure (CP) was applied to simulate the clinical circumstances in which complications related to air leakage are most likely to happen; in the second case, we increase the CP to a relatively high value to simulate the scenarios with an increased risk of complications related to high mucosal pressure. The experiments suggest the follows: (1) Sore throat complications related to high mucosal pressure is most likely to occur in the hypopharynx with a high CP setting, particularly in the areas where the cricoid cartilage presses the mucosa. (2) The narrow hyoid bone super horn width likely causes LMA insertion difficulties. (3) Insufficient CP will significantly increase the risk of air leakage in the oropharynx. A complete sealing pressure line in the contacting surface will be formed with sufficient CP, thereby preventing the air leakage into the oral.

Biomechanical finite element analysis of short-implant-supported, 3-unit, fixed CAD/CAM prostheses in the posterior mandible

Objective

To assess the biomechanical effects of different prosthetic/implant configurations and load directions on 3-unit fixed prostheses supported by short dental implants in the posterior mandible using validated 3-D finite element (FE) models.

Methods

Models represented an atrophic mandible, missing the 2nd premolar, 1st and 2nd molars, and rehabilitated with either two short implants (implant length-IL = 8 mm and 4 mm) supporting a 3-unit dental bridge or three short implants (IL = 8 mm, 6 mm and 4 mm) supporting zirconia prosthesis in splinted or single crowns design. Load simulations were performed in ABAQUS (Dassault Systèmes, France) under axial and oblique (30°) force of 100 N to assess the global stiffness and forces within the implant prosthesis. Local stresses within implant/prosthesis system and strain energy density (SED) within surrounding bone were determined and compared between configurations.

Results

The global stiffness was around 1.5 times higher in splinted configurations vs. single crowns, whereby off-axis loading lead to a decrease of 39%. Splinted prostheses exhibited a better stress distribution than single crowns. Local stresses were larger and distributed over a larger area under oblique loads compared to axial load direction. The forces on each implant in the 2-implant-splinted configurations increased by 25% compared to splinted crowns on 3 implants. Loading of un-splinted configurations resulted in increased local SED magnitude.

Conclusion

Splinting of adjacent short implants in posterior mandible by the prosthetic restoration has a profound effect on the magnitude and distribution of the local stress peaks in peri-implant regions. Replacing each missing tooth with an implant is recommended, whenever bone supply and costs permit.

Experimental validation of a micro-CT finite element model of a human cadaveric mandible rehabilitated with short-implant-supported partial dentures

PURPOSE

This study aimed to address the predictive value of a micro-computed tomography (μCT)-based finite element (μFE) model of a human cadaveric edentulous posterior mandible, rehabilitated by short dental implants. Hereby, three different prosthetic/implant configurations of fixed partial dentures ("Sp"-3 splinted crowns on 3 implants, "Br" - Bridge: 3 splinted crowns on 2 implants, and "Si"- 3 single crowns) were analysed by comparing the computational predictions of the global stiffness with experimental data.

METHODS

Experimental displacement of the bone/implant/prosthesis system was measured under axial and oblique loads of 100 N using an optical deformation system (GOM Aramis) and the overall movement of the testing machine (Zwick Z030). Together with the measured machine force, an "Aramis" (optical markers) and "Zwick" (test machine) stiffness were calculated. FE models were created based on μCT-scans of the cadaveric mandible sample (n = 1) before and after implantation and using stl-files of the crowns. The same load tests and boundary conditions were simulated on the models and the μFE-results were compared to experimental data using linear regression analysis.

RESULTS

The regression line through a plot of pooled stiffness values (N/mm) for the optical displacement recording (true local displacement) and the test machine (machine compliance included) had a slope of 0.57 and a correlation coefficient R² of 0.82. The average pooled correlation of global stiffness between the experiment and FE-analysis (FEA) showed a R² of 0.80, but the FEA-stiffness was 7.2 times higher. The factor was highly dependent on the test configuration. Sp-configuration showed the largest stiffness followed by Br-configuration (17% difference in experiment and 21% in FEA).

CONCLUSIONS

The current study showed good qualitative agreement between the experimental and predicted global stiffness of different short implant configurations. It could be deduced that 1:1 splinting of the short implants by the crowns is most favorable for the stiffness of the implant/prosthesis system. However, in the clinical context, the absolute in silico readings must be interpreted cautiously, as the FEA showed a considerable overestimation of the values.

Bone remodelling in implanted proximal femur using topology optimization and parameterized cellular model

The clinical relevance of bone remodelling predictions calls for accurate finite element (FE) modelling of implant-bone structure and musculoskeletal loading conditions. However, simplifications in muscle loading, material properties, has often been used in FE simulations. Bone adaptation induces changes in bone apparent density and its microstructure. Multiscale simulations, involving optimization methods and biomimetic microstructural models, have proven to be promising for predicting changes in bone morphology. The objective of the study is to develop a novel computational framework to predict bone remodelling around an uncemented femoral implant, using multiscale topology optimization and a parameterized cellular model. The efficacy of the scheme was evaluated by comparing the remodelling predictions with those of isotropic strain energy density (SED) and orthotropy based formulations. The characteristic functional groups and low-density regions of Ward's triangle, predicted by the optimization scheme, were comparable to micro-CT images of the proximal femur. Although the optimization scheme predicted well comparable material distribution in the 2D femur models, the obscured material orientations in some planes of the 3D model indicate the need for a more robust modelling of the boundary conditions. Regression analysis revealed a higher correlation (0.6472) between the topology optimization and SED models than the orthotropic predictions (0.4219). Despite higher bone apposition of 10-20% around the distal tip of the implant, the bone density distributions were well comparable to clinical observations towards the proximal femur. The proposed computational scheme appears to be a viable method for including bone anisotropy in the remodelling formulation.

Effect of Biologically Oriented Preparation Technique on the Stress Concentration of Endodontically Treated Upper Central Incisor Restored with Zirconia Crown: 3D-FEA

The aim of this study was to evaluate the effect of biologically oriented preparation technique on the stress concentration of endodontically treated upper central incisors restored with zirconia crown (yttria-stabilized zirconia polycrystalline ceramic) through finite element analysis (FEA). Four models of maxillary central incisors containing enamel, dentin, periodontal ligament, cortical and medullary bone were created in CAD. Each model received a polymeric core-build up with nanofilled dental resin composite. The evaluated models were SM-preparation in shoulder 90°; CM-chamfer preparation; BOPT-biologically oriented preparation technique and BOPTB-BOPT preparation 1 mm below the cement-enamel junction. All models received zirconia crowns (5Y-TZP), fiberglass post and 1 mm ferrule. The models were imported into the analysis software with parameters for mechanical structural testing using the maximum principal stress and the tensile strength as the analysis criteria. Then, load of 150 N was applied at the cingulum with 45° slope to the long axis of the tooth, with the fixed base for each model. The type of marginal preparation affected the stresses concentration in endodontically treated teeth and in the zirconia crown margin. Considering the stress magnitude only, BOPT is a viable option for anterior monolithic zirconia crowns; however, with the highest stress magnitude at the restoration margin.

In Silico Finite Element Analysis of the Foot Ankle Complex Biomechanics: A Literature Review

Computational approaches, especially finite element analysis (FEA), have been rapidly growing in both academia and industry during the last few decades. FEA serves as a powerful and efficient approach for simulating real-life experiments, including industrial product development, machine design, and biomedical research, particularly in biomechanics and biomaterials. Accordingly, FEA has been a "go-to" high biofidelic software tool to simulate and quantify the biomechanics of the foot-ankle complex, as well as to predict the risk of foot and ankle injuries, which are one of the most common musculoskeletal injuries among physically active individuals. This paper provides a review of the in silico FEA of the foot-ankle complex. First, a brief history of computational modeling methods and finite element (FE) simulations for foot-ankle models is introduced. Second, a general approach to build an FE foot and ankle model is presented, including a detailed procedure to accurately construct, calibrate, verify, and validate an FE model in its appropriate simulation environment. Third, current applications, as well as future improvements of the foot and ankle FE models, especially in the biomedical field, are discussed. Finally, a conclusion is made on the efficiency and development of FEA as a computational approach in investigating the biomechanics of the foot-ankle complex. Overall, this review integrates insightful information for biomedical engineers, medical professionals, and researchers to conduct more accurate research on the foot-ankle FE models in the future.

Characterization and comparison of hyper-viscoelastic properties of normal and osteoporotic bone using stress-relaxation experiment

Bone tissue behavior under various loads is nonlinear elastic due to irreversible energy absorption. Also, viscoelasticity is one of the most important properties of bone which is very important in dynamic analyses and helps a lot in making artificial bone. In this study, rat tibia bone specimens were subjected to compression stress-relaxation test for normal (n = 5) and osteoporotic (n = 5) groups in order to characterize their mechanical properties using finite element modeling coupled with an optimization algorithm. Using this method, the structural equation parameters for the Neo-Hookean model and the Prony series coefficients were used to describe the hyper-elastic and the viscoelastic behavior of specimens, respectively; moreover, the properties of materials including the bulk, shear and Young's moduli for both groups were obtained and compared. The shear modulus was also gained as a function of time. In addition, the percentage of stress reduction and its relation to the initial stress were investigated for specimens. Finally, the effect of changes in each of the parameters of the hyper-viscoelastic structural equation on the force response was determined. Results of this study can be used in predicting the transient response and dynamic analysis of the bone.

Relating simulation studies by provenance-Developing a family of Wnt signaling models

For many biological systems, a variety of simulation models exist. A new simulation model is rarely developed from scratch, but rather revises and extends an existing one. A key challenge, however, is to decide which model might be an appropriate starting point for a particular problem and why. To answer this question, we need to identify entities and activities that contributed to the development of a simulation model. Therefore, we exploit the provenance data model, PROV-DM, of the World Wide Web Consortium and, building on previous work, continue developing a PROV ontology for simulation studies. Based on a case study of 19 Wnt/β-catenin signaling models, we identify crucial entities and activities as well as useful metadata to both capture the provenance information from individual simulation studies and relate these forming a family of models. The approach is implemented in WebProv, a web application for inserting and querying provenance information. Our specialization of PROV-DM contains the entities Research Question, Assumption, Requirement, Qualitative Model, Simulation Model, Simulation Experiment, Simulation Data, and Wet-lab Data as well as activities referring to building, calibrating, validating, and analyzing a simulation model. We show that most Wnt simulation models are connected to other Wnt models by using (parts of) these models. However, the overlap, especially regarding the Wet-lab Data used for calibration or validation of the models is small. Making these aspects of developing a model explicit and queryable is an important step for assessing and reusing simulation models more effectively. Exposing this information helps to integrate a new simulation model within a family of existing ones and may lead to the development of more robust and valid simulation models. We hope that our approach becomes part of a standardization effort and that modelers adopt the benefits of provenance when considering or creating simulation models.

We revise a provenance ontology for simulation studies of cellular biochemical models. Provenance information is useful for understanding the creation of a simulation model because it not only contains information about the entities and activities that have led to a simulation model but also their relations, all of which can be visualized. It provides additional structure by explicitly recording research questions, assumptions, and requirements and relating them along with data, qualitative models, simulation models, and simulation experiments through a small set of predefined but extensible activities.

We have applied our concept to a family of 19 Wnt signaling models and implemented a web-based tool (WebProv) to store the provenance information from these studies. The resulting provenance graph visualizes the story line of simulation studies and demonstrates the creation and calibration of simulation models, the successive attempts of validation and extension, and shows, beyond an individual simulation study, how the Wnt models are related. Thereby, the steps and sources that contributed to a simulation model are made explicit.

Our approach complements other approaches aimed at facilitating the reuse and assessment of simulation products in systems biology such as model repositories as well as annotation and documentation guidelines.

Assessment of paradoxical anterior translation in a CR total knee prosthesis coupling dynamic RSA and FE techniques

Purpose

The study aims were to assess the kinematic data, Internal-External (IE) rotation, and Antero-Posterior (AP) translation of the contact points between the femoral condyles and polyethylene insert and to develop a combined dynamic RSA-FE (Radiostereometric – Finite Element) model that gives results congruent with the literature.

Methods

A cohort of 15 patients who underwent cemented cruciate-retaining highly congruent mobile-bearing total knee arthroplasty were analyzed during a sit-to-stand motor task. The kinematical data from Dynamic RSA were used as input for a patient-specific FE model to calculate condylar contact points between the femoral component and polyethylene insert.

Results

The femoral component showed an overall range about 4 mm of AP translation during the whole motor task, and the majority of the movement was after 40° of flexion. Concerning the IE rotation, the femoral component started from an externally rotate position (− 6.7 ± 10°) at 80° of flexion and performed an internal rotation during the entire motor task. The overall range of the IE rotation was 8.2°.

Conclusions

During the sit to stand, a slight anterior translation from 40° to 0° of flexion of the femoral component with respect to polyethylene insert, which could represent a paradoxical anterior translation. Despite a paradoxical anterior femoral translation was detected, the implants were found to be stable. Dynamic RSA and FE combined technique could provide information about prosthetic component’s stress and strain distribution and the influence of the different designs during the movement.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Variations in Model Development

The use of computational modeling to investigate knee joint biomechanics has increased exponentially over the last few decades. Developing computational models is a creative process where decisions have to be made, subject to the modelers' knowledge and previous experiences, resulting in the "art" of modeling. The long-term goal of the KneeHub project is to understand the influence of subjective decisions on the final outcomes and the reproducibility of computational knee joint models. In this paper, we report on the model development phase of this project, investigating model development decisions and deviations from initial modeling plans. Five teams developed computational knee joint models from the same dataset, and we compared each teams' initial uncalibrated models and their model development workflows. Variations in the software tools and modeling approaches were found, resulting in differences such as the representation of the anatomical knee joint structures in the model. The teams consistently defined the boundary conditions and used the same anatomical coordinate system convention. However, deviations in the anatomical landmarks used to define the coordinate systems were present, resulting in a large spread in the kinematic outputs of the uncalibrated models. The reported differences and similarities in model development and simulation presented here illustrate the importance of the "art" of modeling and how subjective decision-making can lead to variation in model outputs. All teams deviated from their initial modeling plans, indicating that model development is a flexible process and difficult to plan in advance, even for experienced teams.

Reducing the risk of rostral bending failure in Curculio Linnaeus, 1758

With over 300 species worldwide, the genus Curculio Linnaeus, 1758 is a widespread, morphologically diverse lineage of weevils that mainly parasitize nuts. Females use the rostrum, an elongate cuticular extension of the head, to excavate oviposition sites. This process causes extreme bending and deformation of the rostrum, without apparent harm to the structure. The cuticle of the rostral apex exhibits substantial modifications to its composite structure that enhance the elasticity and resiliency of this structure. Here we develop finite element models of the head and rostrum for three Curculio species representing disparate North American clades and rostral morphotypes. The models were subjected to varying apical loads and to prescribed dislocation of the head capsule, with and without representing the cuticular modifications of the rostral apex. We found that the altered layer thicknesses and macrofiber orientation angles of the rostral apex fully explain the observed elasticity of the rostrum. These modifications have a synergistic effect that greatly enhances the flexibility of the rostral apex. Consequently, the cuticle composite profile of the rostral apex substantially mitigates the risk of fracture in dorso-apical flexion. Removing the cuticular modifications, in turn, causes a negative margin of safety for rostral bending, implying strong risk of catastrophic structural failure. The occipital sulci were identified as an important source of biomechanical constraint upon the elasticity of the rostrum, and exhibit the greatest risk of failure within this structure. The apical cuticle profile greatly reduced the maximum stresses and strain energy accumulated in the rostrum, thereby resulting in a positive margin of safety and reducing the risk of fracture. Our findings imply that the primary selective pressure influencing the evolution of the rostral cuticle was most likely negative selection of structural failure caused by bending. STATEMENT OF SIGNIFICANCE: Weevils are among the most diverse and evolutionarily successful animal lineages on Earth. Their success is driven in part by a structure called the rostrum, which gives weevil heads a characteristic "snout-like" appearance. Nut weevils in the genus Curculio use the rostrum to drill holes into developing fruits and nuts, into which they deposit their eggs. During oviposition this exceedingly slender structure is bent into a straightened configuration - in some species up to 90^∘ - but does not suffer any damage during this process. Using finite element models of the rostra of three morphologically distinct species, we show that the Curculio rostrum is only able to withstand repeated, extreme bending because of modifications to the composite structure of the cuticle in the rostral apex. These modifications were shown previously to enhance the intrinsic toughness of the cuticle; in this study, we demonstrate that modification of the rostral cuticle also results in more evenly distributed bending stresses, further reducing the risk of fracture. This is the first time that the laminate profile, orthotropic behavior, and functional gradation of the cuticle have been incorporated into a three-dimensional finite element model of an insect cuticular structure. Our models highlight the significance of biomechanical constraint - i.e., avoidance of catastrophic structural failure - on the evolution of insect morphology.