Rethinking the use of finite element simulations in comparative biomechanics research.
Ontology highlight
ABSTRACT:
SUBMITTER: Tseng ZJ
PROVIDER: S-EPMC8035905 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Simulation-based medicine and the development of complex computer models of biological structures is becoming ubiquitous for advancing biomedical engineering and clinical research. Finite element analysis (FEA) has been widely used in the last few decades to understand and predict biomechanical phenomena. Modeling and simulation approaches in biomechanics are highly interdisciplinary, involving novice and skilled developers in all areas of biomedical engineering and biology. While recent advances in model development and simulation platforms offer a wide range of tools to investigators, the decision making process during modeling and simulation has become more opaque. Hence, reliability of such models used for medical decision making and for driving multiscale analysis comes into question. Establishing guidelines for model development and dissemination is a daunting task, particularly with the complex and convoluted models used in FEA. Nonetheless, if better reporting can be established, researchers will have a better understanding of a model's value and the potential for reusability through sharing will be bolstered. Thus, the goal of this document is to identify resources and considerate reporting parameters for FEA studies in biomechanics. These entail various levels of reporting parameters for model identification, model structure, simulation structure, verification, validation, and availability. While we recognize that it may not be possible to provide and detail all of the reporting considerations presented, it is possible to establish a level of confidence with selective use of these parameters. More detailed reporting, however, can establish an explicit outline of the decision-making process in simulation-based analysis for enhanced reproducibility, reusability, and sharing.
| S-EPMC3278509 | biostudies-literature
Project description:Internal fixation of bone fractures using plates and screws involves many choices-implant type, material, sizes, and geometric configuration-made by the surgeon. These decisions can be important for providing adequate stability to promote healing and prevent implant mechanical failure. The purpose of this study was to develop mathematical models of the relationships between fracture fixation construct parameters and resulting 3D biomechanics, based on parametric computer simulations. Finite element models of hundreds of different locked plate fixation constructs for midshaft diaphyseal fractures were systematically assembled using custom algorithms, and axial, torsional, and bending loadings were simulated. Multivariate regression was used to fit response surface polynomial equations relating fixation design parameters to outputs including maximum implant stresses, axial and shear strain at the fracture site, and construct stiffness. Surrogate models with as little as three regressors showed good fitting (R 2 = 0.62-0.97). Inner working length was the strongest predictor of maximum plate and screw stresses, and a variety of quadratic and interaction terms influenced resulting biomechanics. The framework presented in this study can be applied to additional types of bone fractures to provide clinicians and implant designers with clinical insight, surgical optimization, and a comprehensive mathematical description of biomechanics.
| S-EPMC6358019 | biostudies-literature
Project description:Literature and field data from CIREN database have shown that lumbar spine injuries occur during car crashes. There are multiple hypotheses regarding how they occur; however, there is no biomechanical explanation for these injuries during collisions with road safety barriers (RSBs). Therefore, the objective of this study was to investigate the mechanics of vertebral fractures during car collisions with concrete RSBs. The finite element method was used for the numerical simulations. The global model of the car collision with the concrete RSB was created. The lumbar spine kinematics were extracted from the global simulation and then applied as boundary conditions to the detailed lumbar spine model. The results showed that during the collision, the occupant was elevated, and then dropped during the vehicle landing. This resulted in axial compression forces 2.6 kN with flexion bending moments 34.7 and 37.8 Nm in the L2 and L3 vertebrae. It was shown that the bending moment is the result of the longitudinal force on the eccentricity. The lumbar spine index for the L1-L5 section was 2.80, thus indicating a lumbar spine fracture. The minimum principal strain criterion of 7.4% and damage variable indicated L2 and L3 vertebrae and the inferior part of L1, as those potentially prone to fracture. This study found that lumbar spine fractures could occur as a consequence of vehicle landing during a collision with a concrete RSB mostly affecting the L1-L3 lumbar spine section. The fracture was caused by a combination of axial forces and flexion bending moments.
| S-EPMC8591065 | biostudies-literature
Project description:Computational modeling of tissue-scale cardiac electrophysiology requires numerically converged solutions to avoid spurious artifacts. The steep gradients inherent to cardiac action potential propagation necessitate fine spatial scales and therefore a substantial computational burden. The use of high-order interpolation methods has previously been proposed for these simulations due to their theoretical convergence advantage. In this study, we compare the convergence behavior of linear Lagrange, cubic Hermite, and the newly proposed cubic Hermite-style serendipity interpolation methods for finite element simulations of the cardiac monodomain equation. The high-order methods reach converged solutions with fewer degrees of freedom and longer element edge lengths than traditional linear elements. Additionally, we propose a dimensionless number, the cell Thiele modulus, as a more useful metric for determining solution convergence than element size alone. Finally, we use the cell Thiele modulus to examine convergence criteria for obtaining clinically useful activation patterns for applications such as patient-specific modeling where the total activation time is known a priori.
| S-EPMC4525671 | biostudies-literature
Project description:A finite element analysis was performed to evaluate the stresses around nails and cortical bones in subtrochanteric (ST) fracture models fixed using short cephalomedullary nails (CMNs). A total 96 finite element models (FEMs) were simulated on a transverse ST fracture at eight levels with three different fracture gaps and two different distal locking screw configurations in both normal and osteoporotic bone. All FEMs were fixed using CMNs 200 mm in length. Two distal locking screws showed a wider safe range than 1 distal screw in both normal and osteoporotic bone at fracture gaps ≤ 3 mm. In normal bone FEMs fixed even with two distal locking screws, peak von Mises stresses (PVMSs) in cortical bone and nail constructs reached or exceeded 90% of the yield strength at fracture levels 50 mm and 0 and 50 mm, respectively, at all fracture gaps. In osteoporotic bone FEMs, PVMSs in cortical bone and nail constructs reached or exceeded 90% of the yield strength at fracture levels 50 mm and 0 and 50 mm, respectively, at a 1-mm fracture gap. However, at fracture gaps ≥ 2 mm, PVMSs in cortical bone reached or exceeded 90% of the yield strength at fracture levels ≥ 35 mm. PVMSs in nail showed the same results as 1-mm fracture gaps. PVMSs increased and safe range reduced, as the fracture gap increased. Short CMNs (200 mm in length) with two distal screws may be considered suitable for the fixation of ST transverse fractures at fracture levels 10 to 40 mm below the lesser trochanter in normal bone and 10 to 30 mm in osteoporotic bone, respectively, under the assumptions of anatomical reduction at fracture gap ≤ 3 mm. However, the fracture gap should be shortened to the minimum to reduce the risk of refracture and fixation failure, especially in osteoporotic fractures.
| S-EPMC8248632 | biostudies-literature
Project description:Arthrodesis is a recommended treatment in advanced stages of degenerative disc disease. Despite dynamic fixations were designed to prevent abnormal motions with better physiological load transmission, improving lumbar pain and reducing stress on adjacent segments, contradictory results have been obtained. This study was designed to compare differences in the biomechanical behaviour between the healthy lumbar spine and the spine with DYNESYS and DIAM fixation, respectively, at L4-L5 level. Behaviour under flexion, extension, lateral bending and axial rotation are compared using healthy lumbar spine as reference. Three 3D finite element models of lumbar spine (healthy, DYNESYS and DIAM implemented, respectively) were developed, together a clinical follow-up of 58 patients operated on for degenerative disc disease. DYNESYS produced higher variations of motion with a maximum value for lateral bending, decreasing intradiscal pressure and facet joint forces at instrumented level, whereas screw insertion zones concentrated stress. DIAM increased movement during flexion, decreased it in another three movements, and produced stress concentration at the apophyses at instrumented level. Dynamic systems, used as single systems without vertebral fusion, could be a good alternative to degenerative disc disease for grade II and grade III of Pfirrmann.
| S-EPMC5706716 | biostudies-literature
Project description:Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.
| S-EPMC9166403 | biostudies-literature
Project description:Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.
| S-EPMC9792499 | biostudies-literature
Project description:Human infants are born neurologically immature, potentially owing to conflicting selection pressures between bipedal locomotion and encephalization as suggested by the obstetrical dilemma hypothesis. Australopithecines are ideal for investigating this trade-off, having a bipedally adapted pelvis, yet relatively small brains. Our finite-element birth simulations indicate that rotational birth cannot be inferred from bony morphology alone. Based on a range of pelvic reconstructions and fetal head sizes, our simulations further imply that australopithecines, like humans, gave birth to immature, secondary altricial newborns with head sizes smaller than those predicted for non-human primates of the same body size especially when soft tissue thickness is adequately approximated. We conclude that australopithecines required cooperative breeding to care for their secondary altricial infants. These prerequisites for advanced cognitive development therefore seem to have been corollary to skeletal adaptations for bipedal locomotion that preceded the appearance of the genus Homo and the increase in encephalization.
| S-EPMC9018746 | biostudies-literature
Project description:Masticatory overload on dental implants is one of the causes of marginal bone resorption. The implant-abutment connection (IAC) design plays a critical role in the quality of the stress distribution, and, over the years, different designs were proposed. This study aimed to assess the mechanical behavior of three different types of IAC using a finite element model (FEM) analysis. Three types of two-piece implants were designed: two internal conical connection designs (models A and B) and one internal flat-to-flat connection design (model C). This three-dimensional analysis evaluated the response to static forces on the three models. The strain map, stress analysis, and safety factor were assessed by means of the FEM examination. The FEM analysis indicated that forces are transmitted on the abutment and implant's neck in model B. In models A and C, forces were distributed along the internal screw, abutment areas, and implant's neck. The stress distribution in model B showed a more homogeneous pattern, such that the peak forces were reduced. The conical shape of the head of the internal screw in model B seems to have a keystone role in transferring the forces at the surrounding structures. Further experiments should be carried out in order to confirm the present suppositions.
| S-EPMC6803952 | biostudies-literature