Validation of computational models in biomechanics

Evaluating the stability of external fixators following pelvic injury: A systematic review of biomechanical testing methods

INTRODUCTION

This systematic review aims to identify previously used techniques in biomechanics to assess pelvic instability following pelvic injury, focusing on external fixation constructs.

METHODS

A systematic literature search was conducted to include biomechanical studies and to exclude clinical trials.

RESULTS

Of an initial 4666 studies found, 38 met the inclusion criteria. 84% of the included studies were retrieved from PubMed, Scopus, and Web of Science. The studies analysed 106 postmortem specimens, 154 synthetic bones, and 103 computational models. Most specimens were male (97% synthetic, 70% postmortem specimens). Both the type of injury and the classification system employed varied across studies. About 82% of the injuries assessed were of type C. Two different fixators were tested for FFPII and type A injury, five for type B injury, and fifteen for type C injury. Large variability was observed for external fixation constructs concerning device type and configuration, pin size, and geometry. Biomechanical studies deployed various methods to assess injury displacement, deformation, stiffness, and motion. Thereby, loading protocols differed and inconsistent definitions of failure were determined. Measurement techniques applied in biomechanical test setups included strain gauges, force transducers, and motion tracking techniques.

DISCUSSION AND CONCLUSION

An ideal fixation method should be safe, stable, non-obstructive, and have low complication rates. Although biomechanical testing should ensure that the load applied during testing is representative of a physiological load, a high degree of variability was found in the current literature in both the loading and measurement equipment. The lack of a standardised test design for fixation constructs in pelvic injuries across the studies challenges comparisons between them. When interpreting the results of biomechanical studies, it seems crucial to consider the limitations in cross-study comparability, with implications on their applicability to the clinical setting.

Validation and Estimation of Obesity-Induced Intervertebral Disc Degeneration through Subject-Specific Finite Element Modelling of Functional Spinal Units

(1) Background: Intervertebral disc degeneration has been linked to obesity; its potential mechanical effects on the intervertebral disc remain unknown. This study aimed to develop and validate a patient-specific model of L3-L4 vertebrae and then use the model to estimate the impact of increasing body weight on disc degeneration. (2) Methods: A three-dimensional model of the functional spinal unit of L3-L4 vertebrae and its components were developed and validated. Validation was achieved by comparing the range of motions (RoM) and intradiscal pressures with the previous literature. Subsequently, the validated model was loaded according to the body mass index and estimated stress, deformation, and RoM to assess disc degeneration. (3) Results: During validation, L3-L4 RoM and intradiscal pressures: flexion 5.17° and 1.04 MPa, extension 1.54° and 0.22 MPa, lateral bending 3.36° and 0.54 MPa, axial rotation 1.14° and 0.52 MPa, respectively. When investigating the impact of weight on disc degeneration, escalating from normal weight to obesity reveals an increased RoM, by 3.44% during flexion, 22.7% during extension, 29.71% during lateral bending, and 33.2% during axial rotation, respectively. Also, stress and disc deformation elevated with increasing weight across all RoM. (4) Conclusions: The predicted mechanical responses of the developed model closely matched the validation dataset. The validated model predicts disc degeneration under increased weight and could lay the foundation for future recommendations aimed at identifying predictors of lower back pain due to disc degeneration.

Custom orthotic design by integrating 3D scanning and subject-specific FE modelling workflow

The finite element (FE) foot model can help estimate pathomechanics and improve the customized foot orthoses design. However, the procedure of developing FE models can be time-consuming and costly. This study aimed to develop a subject-specific scaled foot modelling workflow for the foot orthoses design based on the scanned foot surface data. Six participants (twelve feet) were collected for the foot finite element modelling. The subject-specific surface-based finite element model (SFEM) was established by incorporating the scanned foot surface and scaled foot bone geometries. The geometric deviations between the scaled and the scanned foot surfaces were calculated. The SFEM model was adopted to predict barefoot and foot-orthosis interface pressures. The averaged distances between the scaled and scanned foot surfaces were 0.23 ± 0.09 mm. There was no significant difference for the hallux, medial forefoot, middle forefoot, midfoot, medial hindfoot, and lateral hindfoot, except for the lateral forefoot region (p = 0.045). The SFEM model evaluated slightly higher foot-orthoses interface pressure values than measured, with a maximum deviation of 7.1%. These results indicated that the SFEM technique could predict the barefoot and foot-orthoses interface pressure, which has the potential to expedite the process of orthotic design and optimization.

Model Properties and Clinical Application in the Finite Element Analysis of Knee Joint: A Review

The knee is the most complex joint in the human body, including bony structures like the femur, tibia, fibula, and patella, and soft tissues like menisci, ligaments, muscles, and tendons. Complex anatomical structures of the knee joint make it difficult to conduct precise biomechanical research and explore the mechanism of movement and injury. The finite element model (FEM), as an important engineering analysis technique, has been widely used in many fields of bioengineering research. The FEM has advantages in the biomechanical analysis of objects with complex structures. Researchers can use this technology to construct a human knee joint model and perform biomechanical analysis on it. At the same time, finite element analysis can effectively evaluate variables such as stress, strain, displacement, and rotation, helping to predict injury mechanisms and optimize surgical techniques, which make up for the shortcomings of traditional biomechanics experimental research. However, few papers introduce what material properties should be selected for each anatomic structure of knee FEM to meet different research purposes. Based on previous finite element studies of the knee joint, this paper summarizes various modeling strategies and applications, serving as a reference for constructing knee joint models and research design.

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.

Endoscopic versus laparoscopic bariatric procedures: A computational biomechanical study through a patient-specific approach

BACKGROUND AND OBJECTIVES

Within the framework of computational biomechanics, finite element models of the gastric district could be seen as a potential clinical tool not only to study the effects apported by bariatric surgery, but also to compare different surgical techniques such as the new emerging Endoscopic Sleeve Gastroplasty (ESG) with respect to well-established ones (such as the Laparoscopic Sleeve Gastrectomy, LSG).

METHODS

This work realized a fully computational comparison between the outcomes obtained from 10 patient-specific stomach models, which were used to simulate ESG, and the complementary results obtained from models representing the post-LSG of the same subjects. Specifically, once the ESG was simulated, a mechanical stimulus was applied by increasing an intragastric pressure up to a maximum of 5 kPa, in order to replicate the process of food intake, as well as for post-LSG models.

RESULTS

Results revealed non negligible differences between the techniques also within the same subject. In particular, not only LSG could lead to a greater reduction in the stomach volume (about 77 % at baseline, which is strictly linked to weight loss), but also influence the gastric distension (12 % less than pre-operative models). On the contrary, if ESG would be performed, a more similar pre-operative mechanical stimulation of the gastric walls may be seen (difference of about 1 %), thus preserving the mechanosensation, but the detriment of the volume reduction (about 56 % at baseline, and even decreases with increasing pressure). Moreover, since results suggested ESG may be more influenced by the pre-operative gastric cavity than LSG, a predictive model was proposed to support the surgical planning and the estimation of the volume reduction after ESG.

CONCLUSIONS

ESG and LSG have substantial differences in their protocols and post-surgical effects. This work pointed out that variations between the two procedures may be observed also from a computational point of view, especially when including patient-specific geometries. These insights support gastric modelling as a valuable tool to evaluate, design and critically compare emerging bariatric surgical procedures, not only from empirical aspects and clinical outcomes, but also from a mechanical point of view.

Reproducibility in modeling and simulation of the knee: Academic, industry, and regulatory perspectives

Stakeholders in the modeling and simulation (M&S) community organized a workshop at the 2019 Annual Meeting of the Orthopaedic Research Society (ORS) entitled "Reproducibility in Modeling and Simulation of the Knee: Academic, Industry, and Regulatory Perspectives." The goal was to discuss efforts among these stakeholders to address irreproducibility in M&S focusing on the knee joint. An academic representative from a leading orthopedic hospital in the United States described a multi-institutional, open effort funded by the National Institutes of Health to assess model reproducibility in computational knee biomechanics. A regulatory representative from the United States Food and Drug Administration indicated the necessity of standards for reproducibility to increase utility of M&S in the regulatory setting. An industry representative from a major orthopedic implant company emphasized improving reproducibility by addressing indeterminacy in personalized modeling through sensitivity analyses, thereby enhancing preclinical evaluation of joint replacement technology. Thought leaders in the M&S community stressed the importance of data sharing to minimize duplication of efforts. A survey comprised 103 attendees revealed strong support for the workshop and for increasing emphasis on computational modeling at future ORS meetings. Nearly all survey respondents (97%) considered reproducibility to be an important issue. Almost half of respondents (45%) tried and failed to reproduce the work of others. Two-thirds of respondents (67%) declared that individual laboratories are most responsible for ensuring reproducible research whereas 44% thought that journals are most responsible. Thought leaders and survey respondents emphasized that computational models must be reproducible and credible to advance knee M&S.

The effect of different posterior inclinations of tibial component on tibiofemoral contact pressures after unicompartmental knee arthroplasty

BACKGROUND

Different posterior inclinations of tibial component after unicompartmental knee arthroplasty (UKA) may lead to different biomechanical characteristics of the knee joint. This finite element study was designed to investigate the tibiofemoral contact pressures after UKA with different posterior inclinations of tibial component.

METHODS

Finite element model of a healthy knee joint was constructed, and mobile-bearing (MB) UKA models with 5 different posterior inclinations (3°, 5°, 7°, 9° and 11°) of tibial components were simulated. The maximum contact pressures of tibial plateau cartilage in the lateral compartment and polyethylene insert in the medial compartment were calculated based on the ground reaction force and the angle of the knee flexion obtained by 3D motion capture system.

RESULTS

The loading ratio of medial and lateral compartments during standing stance (medial 54.49%, lateral 45.51%) and tibial anterior displacement (134 N, 3.89 mm) of healthy knee was basically consistent with previous experimental data. The maximum contact pressures of the medial meniscus and lateral tibial plateau cartilage of the healthy knee during standing stance were 2.14 MPa and 1.57 MPa, respectively. At the static standing phase, the maximum contact pressures of the polyethylene insert decreased from 17.90 to 17.29 Mpa, and the maximum contact pressures of the tibial plateau cartilage in the lateral compartment increased from 0.81 to 0.92 Mpa following an increase in the posterior inclination of the tibial component. At the first peak of ground reaction force, the maximum contact pressures of polyethylene insert increased from 22.37 to 25.16 MPa, and the maximum contact pressures of tibial plateau cartilage in the lateral compartment increased from 3.03 to 3.33 MPa, with the increase in the posterior inclination of the tibial component. At the second peak of ground reaction force, the maximum contact pressures of polyethylene insert decreased from 2.34 to 2.22 MPa with the increase in posterior inclination of tibial component.

CONCLUSION

The preoperative and postoperative finite element models of MB UKA were well established. The results showed that the maximum contact pressures of the polyethylene insert did not change significantly with the increase in the posterior inclination of the tibial prosthesis, while the maximum contact pressures of the tibial plateau cartilage of the lateral compartment increased when the posterior inclination of the tibial prosthesis was > 7°. Our results also show that the maximum contact pressures were greater with an excessive inclination angle (11°) of the tibial component, and the pressures of the tibial plateau cartilage in the lateral compartment were more concentrated on the posterior area. This study, therefore, proposes that excessive osteotomy should be avoided.

Modelling for design and evaluation of industrial exoskeletons: A systematic review

Industrial exoskeletons are developed to relieve workers' physical demands in the workplace and to alleviate ergonomic issues associated with work-related musculoskeletal disorders. As a safe and economical alternative to empirical/experimental methods, modelling is considered as a powerful tool for design and evaluation of industrial exoskeletons. This systematic review aims to provide a comprehensive understanding of the current literature on the design and evaluation of industrial exoskeletons through modelling. A systematic study was conducted by general keyword searches of five electronic databases over the last two decades (2003-2022). Out of the 701 records initially retrieved, 33 eligible articles were included and analyzed in the final review, presenting a variety of model inputs, model development, and model outputs used in the modelling. This systematic review study revealed that existing modelling methods can evaluate the biomechanical and physiological effects of industrial exoskeletons and provide some design parameters. However, the modelling method is currently unable to cover some of the main evaluation metrics supported by experimental assessments, such as task performance, user experience/discomfort, change in metabolic costs etc. Standard guidelines for model construction and implementation, as well as validation of human-exoskeleton interactions, remain to be established.

Definition and evaluation of a finite element model of the human heel for diabetic foot ulcer prevention under shearing loads

Diabetic foot ulcers are triggered by mechanical loadings applied to the surface of the plantar skin. Strain is considered to play a crucial role in relation to ulcer etiology and can be assessed by Finite Element (FE) modeling. A difficulty in the generation of these models is the choice of the soft tissue material properties. In the literature, many studies attempt to model the behavior of the heel soft tissues by implementing constitutive laws that can differ significantly in terms of mechanical response. Moreover, current FE models lack of proper evaluation techniques that could estimate their ability to simulate realistic strains. In this article, we propose and evaluate a FE model of the human heel for diabetic foot ulcer prevention. Soft tissue constitutive laws are defined through the fitting of experimental stretch-stress curves published in the literature. The model is then evaluated through Digital Volume Correlation (DVC) based on non-rigid 3D Magnetic Resonance Image Registration. The results from FE analysis and DVC show similar strain locations in the fat pad and strain intensities according to the type of applied loads. For additional comparisons, different sets of constitutive models published in the literature are applied into the proposed FE mesh and simulated with the same boundary conditions. In this case, the results in terms of strains show great diversity in locations and intensities, suggesting that more research should be developed to gain insight into the mechanical properties of these tissues.

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.

Best practices for exposure model peer review - A SciPinion advisory panel report

The extent and rigor of peer review that a model undergoes during and after development influences the confidence of users and managers in model predictions. A process for determining the breadth and depth of peer review of exposure models was developed with input from a panel of exposure-modeling experts. This included consideration of the tiers and types of models (e.g., screening, deterministic, probabilistic, etc.). The experts recommended specific criteria be considered when evaluating the degree to which a model has been peer reviewed, including quality of documentation and the model peer review process (e.g., internal review with a regulatory agency by subject matter experts, expert review reports, formal Scientific Advisory Panels, and journal peer review). In addition, because the determination of the confidence level for an exposure model's predictions is related to the degree of evaluation the model has undergone, irrespective of peer review, the experts recommended the approach include judging the degree of model rigor using a set of specific criteria: (1) nature and quality of input data, (2) model verification, (3) model corroboration, and (4) model evaluation. Other key areas considered by the experts included recommendations for addressing model uncertainty and sensitivity, defining the model domain of applicability, and flags for when a model is used outside its domain of applicability. The findings of this expert engagement will help developers as well as users of exposure models have greater confidence in their application and yield greater transparency in the evaluation and peer review of exposure models.

Biophysical models of early mammalian embryogenesis

Embryo development is a critical and fascinating stage in the life cycle of many organisms. Despite decades of research, the earliest stages of mammalian embryogenesis are still poorly understood, caused by a scarcity of high-resolution spatial and temporal data, the use of only a few model organisms, and a paucity of truly multidisciplinary approaches that combine biological research with biophysical modeling and computational simulation. Here, we explain the theoretical frameworks and biophysical processes that are best suited to modeling the early mammalian embryo, review a comprehensive list of previous models, and discuss the most promising avenues for future work.

Dynamic finite element analyses to compare the influences of customised total talar replacement and total ankle arthroplasty on foot biomechanics during gait

Objective, Total talar replacement (TTR) using a customised talus prosthesis is an emerging surgical alternative to conventional total ankle arthroplasty (TAA) for treating ankle problems. Upon satisfying clinical reports in the literature, this study explored the advantages of TTR in restoring foot biomechanics during walking compared with TAA through computational simulations.Methods, A dynamic finite element foot model was built from the MRIs of a healthy participant and modified into two implanted counterparts (TTR and TAA) by incorporating the corresponding prosthetic components into the ankle joint. Twenty bony parts, thirty-nine ligament/tendon units, nine muscle contractors, and bulk soft tissue were included in the intact foot model. The TTR prosthesis was reconstructed from the mirror image data of the participant's contralateral talus and the TAA prosthesis was modelled by reproducing the Scandinavian ankle replacement procedure in the model assembly. The model was meshed with explicit deformable elements and validated against existing experimental studies that have assessed specific walking scenarios. Simulations were performed using the boundary conditions (time-variant matrix of muscle forces, segment orientation, and ground reaction forces) derived from motion capture analyses and musculoskeletal modelling of the participant's walking gait. Outcome variables, including foot kinematics, joint loading, and plantar pressure were reported and compared among the three model conditions.

RESULTS

Linear regression indicated a better agreement between the TTR model and intact foot model in plots of joint motions and foot segment movements during walking (R² ​= ​0.721-0.993) than between the TAA and intact foot (R² ​= ​0.623-0.990). TAA reduced talocrural excursion by 21.36%-31.92% and increased (MTP) dorsiflexion by 3.03%. Compared with the intact foot, TTR and TAA increased the midtarsal joint contact force by 17.92% and 10.73% respectively. The proximal-to-distal force transmission within the midfoot was shifted to the lateral column in TTR (94.52% or 210.54 ​N higher) while concentrated on the medial column in TAA (41.58% or 27.55 ​N higher). The TTR produced a plantar pressure map similar to that of the intact foot. TAA caused the plantar pressure centre to drift medially and increased the peak forefoot pressure by 7.36% in the late stance.

CONCLUSION

The TTR better reproduced the foot joint motions, segment movements, and plantar pressure map of an intact foot during walking. TAA reduced ankle mobility while increasing movement of the adjacent joints and forefoot plantar pressure. Both implant methods changed force transmission within the midfoot during gait progression.The translational potential of this article Our work is one of the few to report foot segment movements and the internal loading status of implanted ankles during a dynamic locomotion task. These outcomes partially support the conjecture that TTR is a prospective surgical alternative for pathological ankles from a biomechanical perspective. This study paves the way for further clinical investigations and systematic statistics to confirm the effects of TTR on functional joint recovery.

Kinematic and biomechanical responses of the spine to distraction surgery in children with early onset scoliosis: A 3-D finite element analysis

Periodical and consecutive distraction is an effective treatment for severe early onset scoliosis (EOS), which enables the spinal coronal and sagittal plane deformity correction. However, the rate of rod fractures and postoperative complications was still high mainly related to the distraction process. Previous studies have primarily investigated the maximum safe distraction force without a rod broken, neglecting the spinal re-imbalance and distraction energy consumption, which is equally vital to evaluate the operative value. This study aimed to reveal the kinematic and biomechanical responses occurring after spinal distraction surgery, which were affected by traditional bilateral fixation. The spinal models (C6-S1) before four distractions were reconstructed based on CT images and the growing rods were applied with the upward displacement load of 0–25 mm at an interval of 5 mm. Relationships between the distraction distance, the distraction force and the thoracic and lumbar Cobb angle were revealed, and the spinal displacement and rotation in three-dimensional directions were measured. The spinal overall imbalance would also happen during the distraction process even under the safe force, which was characterized by unexpected cervical lordosis and lateral displacement. Additionally, the law of diminishing return has been confirmed by comparing the distraction energy consumption in different distraction distances, which suggests that more attention paid to the spinal kinematic and biomechanical changes is better than to the distraction force. Notably, the selection of fixed segments significantly impacts the distraction force at the same distraction distance. Accordingly, some results could provide a better understanding of spinal distraction surgery.

Musculoskeletal Modeling of the Wrist via a Multi Body Simulation

In this study, three different musculoskeletal modeling approaches were compared to each other. The objective was to show the possibilities in the case of a simple mechanical model of the wrist, using a simple multi-body-simulation (MBS) model, and using a more complex and patient-specific adaptable wrist joint MBS model. Musculoskeletal modeling could be a useful alternative, which can be practiced as a non-invasive approach to investigate body motion and internal loads in a wide range of conditions. The goal of this study was the introduction of computer-based modelling of the physiological wrist with (MBS-) models focused on the muscle and joint forces acting on the wrist.

Reproducing Morphological Features Of Intervertebral Disc Using Finite Element Modeling To Predict The Course Of Cervical Spine Dorsopathy

Study objective — To evaluate how morphological features of intervertebral disc would affect the outcomes of finite element modeling of axial load in the cervical spine, C3-C5, in order to predict the risk of occurrence and course of dorsopathies. Material and Methods — Three-dimensional models of the cervical spine vertebrae were generated from the computed tomography data of a volunteer (24 years old male without detected pathology of his neck). Intervertebral disc models were developed in two configurations. For each model, we performed a finite element investigation of the stress-strain state with the same loading conditions. The load-displacement curves were compared with the experimental data generated from the results of previously conducted in vitro experiments. Results — The maximum and mean displacement values for the isotropic model were 1.15 mm and 0.73 ± 0.45 mm, respectively. For anisotropic model, maximum and mean displacement values were 0.86 mm and 0.47 ± 0.24 mm, correspondingly. Predicted displacement values for both models matched the experimental data fairly well. Stress profiles of intervertebral discs and stress diagrams of facet joints were calculated. Conclusion — The proposed geometric and constitutive configurations of the intervertebral disc take into account specific morphological features at low computational costs, thereby facilitating the modeling of degenerative disc changes.

Finite Element Modelling of Planus and Rectus Foot Types for the Study of First Metatarsophalangeal and First Metatarsocuneiform Joint Contact Mechanics

Evaluating the contact mechanics of human joints is an important element in understanding the pathomechanics of orthopaedic diseases. Although physical testing is essential in the evaluation process, reliable computational models can augment these experiments by non-invasive predictions of biomechanical or surgical variables. The objective of this study was to perform verification of a framework for developing a medial forefoot finite element. Verification was conducted by comparing computational predictions to experimental measurements of first metatarsophalangeal and first metatarsocuneiform joint contact mechanics. A custom-built force-controlled cadaveric test-rig was used to derive measurements of contact pressure, force, and area. A quasi-static finite element was developed and driven under the same boundary and loading conditions. Calibration of cartilage moduli and mesh sensitivity analyses were performed. Mean errors in contact pressures, forces, and areas were 24%, 4%, and 40% at the first metatarsophalangeal joint and 23%, 12%, and 19% at the first Metatarsocuneiform joint, respectively. Verification of a medial forefoot finite element model development framework was presented and found to be within 30% for contact pressure and contact force of both joints. This study presents a method to verify and simulate realistic physiological loading to investigate orthopaedic diseases of the medial forefoot.

Mathematical Modeling of Proliferative Immune Response Initiated by Interactions Between Classical Antigen-Presenting Cells Under Joint Antagonistic IL-2 and IL-4 Signaling

During an adaptive immune response from pathogen invasion, multiple cytokines are produced by various immune cells interacting jointly at the cellular level to mediate several processes. For example, studies have shown that regulation of interleukin-4 (IL-4) correlates with interleukin-2 (IL-2) induced lymphocyte proliferation. This motivates the need to better understand and model the mechanisms driving the dynamic interplay of proliferation of lymphocytes with the complex interaction effects of cytokines during an immune response. To address this challenge, we adopt a hybrid computational approach comprising of continuous, discrete and stochastic non-linear model formulations to predict a system-level immune response as a function of multiple dependent signals and interacting agents including cytokines and targeted immune cells. We propose a hybrid ordinary differential equation-based (ODE) multicellular model system with a stochastic component of antigen microscopic states denoted as Multiscale Multicellular Quantitative Evaluator (MMQE) implemented using MATLAB. MMQE combines well-defined immune response network-based rules and ODE models to capture the complex dynamic interactions between the proliferation levels of different types of communicating lymphocyte agents mediated by joint regulation of IL-2 and IL-4 to predict the emergent global behavior of the system during an immune response. We model the activation of the immune system in terms of different activation protocols of helper T cells by the interplay of independent biological agents of classic antigen-presenting cells (APCs) and their joint activation which is confounded by the exposure time to external pathogens. MMQE quantifies the dynamics of lymphocyte proliferation during pathogen invasion as bivariate distributions of IL-2 and IL-4 concentration levels. Specifically, by varying activation agents such as dendritic cells (DC), B cells and their joint mechanism of activation, we quantify how lymphocyte activation and differentiation protocols boost the immune response against pathogen invasion mediated by a joint downregulation of IL-4 and upregulation of IL-2. We further compare our in-silico results to in-vivo and in-vitro experimental studies for validation. In general, MMQE combines intracellular and extracellular effects from multiple interacting systems into simpler dynamic behaviors for better interpretability. It can be used to aid engineering of anti-infection drugs or optimizing drug combination therapies against several diseases.

Posterior spinal stabilization: A biomechanical comparison of Laminar Hook Fusion to a Pedicle Screw System

BACKGROUND

Several spine instrumentation techniques have been introduced to correct inter-segmental alignment, or provide long-term stability. Whilst pedicle screws are considered the intervention of reference, we hypothesize that the week hold of osteoporotic bone, might be a clinical indicator for an alternative surgical approach.

METHODS

To put this to the test, a non-linear Finite Element model, of a ligamentous lumbosacral spine, was employed to examine a stabilization spanning over L3-L5. Two different immobilization techniques (a Pedicle Screw System and Laminar Hook Fusion) are compared as to their biomechanical response during 7.5 Nm flexion, lateral flexion and torsion, while considering a 280 N follower load. Fifteen analyses performed in total, simulating patients of healthy and osteoporotic Bone Mineral Density.

FINDINGS

Range of Motion was significantly reduced after instrumentation for both implant systems. This trend was more pronounced in the Pedicle Screw models, which were stressed to a higher degree. To evaluate implant loosening risk, we introduce the consideration of strain energy patterns around the screw tract. The notably higher intensity of these, for the osteoporotic model, taken into consideration with the weaker strength of the tissue and inconsistencies in the stress allocation between implant and bone, affirmed an increased risk for loosening of the Pedicle Screws in osteoporotic patients.

INTERPRETATION

The analysis provided refined insight as to the treatment of osteoporotic patients as well as to their postoperative care, as restriction of specific movements (e.g. through bracing), could significantly restrict the stress values in the bone-implant interface and thus, reduce implant failure.

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.

Biomechanical analysis of the cervical spine segment as a method for studying the functional and dynamic anatomy of the human neck

BACKGROUND

Traditionally, dynamic and functional anatomy, in particular the dynamic anatomy of the neck, is studied on cadaveric material. However, the development of in vivo visualization technologies and in silico modeling has made it possible to expand these possibilities. Despite significant progress in the study of dynamic and functional anatomy of the neck by means of in silico methods, the issues of validating the developed models and taking into account the pronounced nonlinearity of soft tissues as well as local anisotropy remain open. The aim of this study was to develop a virtual dynamic anatomical model of the human neck and reproduce the dynamic processes in the cervical spine from this model using the finite element method.

MATERIALS AND METHODS

Reverse engineering was used to generate a dynamic anatomical model of the neck from CT data (both male, 24 and 22 years old). Two segments of the cervical spine (C3-C5, C2-T1) were isolated from the resulting model for finite element analysis. Finite element mesh generation and contact interactions were performed using the HyperMesh software (Altair Engineering Inc, Troy, Michigan, USA). The anisotropic hyperelastic Holzapfel-Gasser-Ogden model was used to describe the material behavior of the fibrous rings of the disc. Material modeling and finite element analysis were performed using Abaqus CAE 6.14 software (Simulia, Johnston, Rhode Island, USA).

RESULTS

A technique for creating a virtual dynamic anatomical model of the neck was elaborated and implemented. The model includes 79 major anatomical structures of the neck segmented from radiological data. A finite element analysis of the cervical spine was performed. The results of finite element analysis of the C3-C5 segment under axial load were compared with in vitro data. The proposed model shows nonlinear deformation of the disc under static loading; the model predicted displacement values agree well with the experimental ones. The displacement of the С3-С5 central vertebra with an axial load of 800 N reaches a value of 0.65 mm. For the segment C2-T1, data on intradiscal pressure, stress plots and displacements during flexion were obtained. The maximum stress value of 10.036 MPa is observed in the C3-C4 disc.

CONCLUSION

Simulation results using the proposed methodology are in good agreement with experimental data. The generated biomechanical models allow describing dynamic phenomena in the cervical spine and obtaining a wide range of quantitative properties of anatomical objects, which are otherwise inaccessible to classical methods for studying dynamic and functional anatomy.

Biomechanical comparison among five mid/hindfoot arthrodeses procedures in treating flatfoot using a musculoskeletal multibody driven finite element model

BACKGROUND AND OBJECTIVE

Mid/hindfoot arthrodesis could modify the misalignment of adult-acquired flatfoot and attenuate pain. However, the long-term biomechanical effects of these surgical procedures remain unclear, and the quantitative evidence is scarce. Therefore, we aimed to investigate and quantify the influences of five mid/hindfoot arthrodeses on the internal foot biomechanics during walking stance.

METHODS

A young participant with flexible flatfoot was recruited for this study. We reconstructed a subject-specific musculoskeletal multibody driven-finite element (FE) foot model based on the foot magnetic resonance imaging. The severe flatfoot model was developed from the flexible flatfoot through the attenuation of ligaments and the unloading of the posterior tibial muscle. The five mid/hindfoot arthrodeses simulations (subtalar, talonavicular, calcaneocuboid, double, and triple arthrodeses) and a control condition (no arthrodesis) were performed simultaneously in the detailed foot multibody dynamics model and FE model. Muscle forces calculated by a detailed multi-segment foot model and ground reaction force were used to drive the foot FE model. The internal foot loadings were compared among control and these arthrodeses conditions at the first and second vertical ground reaction force (VGRF) peak and VGRF valley instants.

RESULTS

The results indicated that the navicular heights in double and triple arthrodeses were higher than other surgical procedures, while the subtalar arthrodesis had the smallest values. Five mid/hindfoot arthrodeses reduced the peak plantar fascia stress compared to control. However, double and triple arthrodeses increased the peak medial cuneo-navicular joint contact pressures and peak foot pressures as well as the metatarsal bones stresses.

CONCLUSION

Although mid/hindfoot arthrodesis generally reduced the collapse of medial longitudinal arch and plantar fascia loading during the stance phase, the increased loading in the adjacent unfused joint and metatarsal bones for double and triple arthrodeses should be noted. These findings could account for some symptoms experienced by flatfoot patients after surgery, which may facilitate the optimization of surgical protocols.

Assessing the biomechanical properties of nitinol staples in normal, osteopenic and osteoporotic bone models: A finite element analysis

OBJECTIVE

Bone staples are internal fixation devices that are frequently used in the foot, ankle, and hand to provide stabilization. Fixation stability is vital after fusion or fracture surgeries to ensure proper bone healing. Patients undergoing surgeries that require fixation to keep bones aligned and stable may present with diminishing bone mechanical properties, and this may compromise the ability of the fixation hardware to maintain a stable construct. The purpose of this study was to investigate the mechanical performance of shape memory and superelastic nitinol bone staples with different bridge geometries in normal, osteopenic, and osteoporotic bone models. Contact forces and maximum principal stress and strain in the bone were recorded.

METHODS

Finite element simulations of a bone staple fixation procedure were performed to examine the initial and post-surgery contact force, as well as the maximum principal stress and strain of 15 mm bridge and 20 mm bridge staple-bone constructs.

RESULTS

Shape memory nitinol staples exhibited higher contact forces compared to superelastic nitinol staples. Nitinol bone staples with 20 mm bridge lengths displayed higher contact forces and lower stresses in all bone types, as well as lower strains in osteoporotic bone models compared to nitinol staples with a 15 mm bridge length.

CONCLUSION

Nitinol bone staple constructs with 20 mm bridge length staples provide higher contact forces and display lower stresses in the bone than 15 mm bridge staple-bone constructs, which may be beneficial in bone with diminishing mechanical properties. Both superelastic and shape memory effect nitinol staples provide adequate compression and stress relief. However, if osteopenia is present, shape memory effect nitinol staples with a 20 mm bridge length may provide more stress relief and compression, if the bone anatomy allows.

Computational modelling of the scoliotic spine: A literature review

Scoliosis is a deformity of the spine that in severe cases requires surgical treatment. There is still disagreement among clinicians as to what the aim of such treatment is as well as the optimal surgical technique. Numerical models can aid clinical decision-making by estimating the outcome of a given surgical intervention. This paper provided some background information on the modelling of the healthy spine and a review of the literature on scoliotic spine models, their validation, and their application. An overview of the methods and techniques used to construct scoliotic finite element and multibody models was given as well as the boundary conditions used in the simulations. The current limitations of the models were discussed as well as how such limitations are addressed in non-scoliotic spine models. Finally, future directions for the numerical modelling of scoliosis were addressed.

Finite Element Analysis of Fracture Fixation

PURPOSE OF REVIEW

Fracture fixation aims to provide stability and promote healing, but remains challenging in unstable and osteoporotic fractures with increased risk of construct failure and nonunion. The first part of this article reviews the clinical motivation behind finite element analysis of fracture fixation, its strengths and weaknesses, how models are developed and validated, and how outputs are typically interpreted. The second part reviews recent modeling studies of the femur and proximal humerus, areas with particular relevance to fragility fractures.

RECENT FINDINGS

There is some consensus in the literature around how certain modeling aspects are pragmatically formulated, including bone and implant geometries, meshing, material properties, interactions, and loads and boundary conditions. Studies most often focus on predicted implant stress, bone strain surrounding screws, or interfragmentary displacements. However, most models are not rigorously validated. With refined modeling methods, improved validation efforts, and large-scale systematic analyses, finite element analysis is poised to advance the understanding of fracture fixation failure, enable optimization of implant designs, and improve surgical guidance.

Finite element modelling and experimental validation of a total implanted shoulder joint

Background Replicating a total shoulder arthroplasty in laboratory is a difficult task due to complex geometry of the structures and degrees of freedom of the joint. Implanted joint shoulders have been investigated using numerical tools, but models developed lack of experimental validation. The objective of this study was to develop a finite element model that replicated correctly an experimental simulator of an implanted joint shoulder based on the comparison of measured and calculated strains. The methods used include a non-cemented Anatomical Comprehensive© Total Shoulder System that was implanted in 4^th generation composite bones. The finite element model designed replicates adequately the experimental model. Both models included the most important muscles of shoulder abduction and the same boundary conditions (loads, fixation, and interface conditions). Strain gauge rosettes were used to measure strain responses on the shoulder in 90° abduction. The results of linear regression analysis between numerical and experimental results present a high correlation coefficient of 0.945 and a root-mean-square-error of 35 µε, suggesting adequate agreement between the experimental and the numerical models. Small strains were obtained and changes in load distribution from posterior to anterior region were observed. As conclusion we can say that the experiments allowed good replication of the finite element model, and the use of strain gauges is suitable for numerical-experimental validation of bone joints.

Measuring and modeling in vivo human gracilis muscle-tendon unit length

Musculoskeletal models rely heavily on the use of an anatomical dataset and clearly defined assumptions to accurately model the subject being studied. Therefore, it is important to understand the limitations of using musculoskeletal models to study individuals. This paper describes a method of measuring in vivo gracilis muscle-tendon unit length and presents a comparison of experimental data versus predictions from four musculoskeletal models in OpenSim. The largest errors occurred when the knee was fully extended. At this position, the absolute average muscle-tendon unit length error was 7% and the absolute average fiber length error was between 15% and 32%. However, the variability of these errors was significant. Manual linear scaling based on an anthropometric database did not capture the variability observed in subjects. The fiber length errors observed are predicted to have a significant impact on muscle force production that may not represent true subject specific force-length relationship of the gracilis.

A review of foot finite element modelling for pressure ulcer prevention in bedrest: Current perspectives and future recommendations

Pressure ulcers (PUs) are a major public health challenge, having a significant impact on healthcare service and patient quality of life. Computational biomechanical modelling has enhanced PU research by facilitating the investigation of pressure responses in subcutaneous tissue and skeletal muscle. Extensive work has been undertaken on PUs on patients in the seated posture, but research into heel ulcers has been relatively neglected. The aim of this review was to address the key challenges that exist in developing an effective FE foot model for PU prevention and the confusion surrounding the wide range of outputs reported. Nine FE foot studies investigating heel ulcers in bedrest were identified and reviewed. Six studies modelled the posterior part of the heel, two included the calf and foot, and one modelled the whole body. Due to the complexity of the foot anatomy, all studies involved simplification or assumptions regarding parts of the foot structure, boundary conditions and material parameters. Simulations aimed to understand better the stresses and strains exhibited in the heel soft tissues of the healthy foot. The biomechanical properties of soft tissue derived from experimental measurements are critical for developing a realistic model and consequently guiding clinical decisions. Yet, little to no validation was reported in each of the studies. If FE models are to address future research questions and clinical applications, then sound verification and validation of these models is required to ensure accurate conclusions and prediction of patient outcomes. Recommendations and considerations for future FE studies are therefore proposed.

Biomechanical assessment of different fixation methods in mandibular high sagittal oblique osteotomy using a three-dimensional finite element analysis model

With modern-day technical advances, high sagittal oblique osteotomy (HSOO) of the mandible was recently described as an alternative to bilateral sagittal split osteotomy for the correction of mandibular skeletal deformities. However, neither in vitro nor numerical biomechanical assessments have evaluated the performance of fixation methods in HSOO. The aim of this study was to compare the biomechanical characteristics and stress distribution in bone and osteosynthesis fixations when using different designs and placing configurations, in order to determine a favourable plating method. We established two finite element models of HSOO with advancement (T1) and set-back (T2) movements of the mandible. Six different configurations of fixation of the ramus, progressively loaded by a constant force, were assessed for each model. The von Mises stress distribution in fixations and in bone, and bony segment displacement, were analysed. The lowest mechanical stresses and minimal gradient of displacement between the proximal and distal bony segments were detected in the combined one-third anterior- and posterior-positioned double mini-plate T1 and T2 models. This suggests that the appropriate method to correct mandibular deformities in HSOO surgery is with use of double mini-plates positioned in the anterior one-third and posterior one-third between the bony segments of the ramus.

Predictive Simulations of Musculoskeletal Function and Jumping Performance in a Generalized Bird

Jumping is a common, but demanding, behavior that many animals employ during everyday activity. In contrast to jump-specialists such as anurans and some primates, jumping biomechanics and the factors that influence performance remains little studied for generalized species that lack marked adaptations for jumping. Computational biomechanical modeling approaches offer a way of addressing this in a rigorous, mechanistic fashion. Here, optimal control theory and musculoskeletal modeling are integrated to generate predictive simulations of maximal height jumping in a small ground-dwelling bird, a tinamou. A three-dimensional musculoskeletal model with 36 actuators per leg is used, and direct collocation is employed to formulate a rapidly solvable optimal control problem involving both liftoff and landing phases. The resulting simulation raises the whole-body center of mass to over double its standing height, and key aspects of the simulated behavior qualitatively replicate empirical observations for other jumping birds. However, quantitative performance is lower, with reduced ground forces, jump heights, and muscle–tendon power. A pronounced countermovement maneuver is used during launch. The use of a countermovement is demonstrated to be critical to the achievement of greater jump heights, and this phenomenon may only need to exploit physical principles alone to be successful; amplification of muscle performance may not necessarily be a proximate reason for the use of this maneuver. Increasing muscle strength or contractile velocity above nominal values greatly improves jump performance, and interestingly has the greatest effect on more distal limb extensor muscles (i.e., those of the ankle), suggesting that the distal limb may be a critical link for jumping behavior. These results warrant a re-evaluation of previous inferences of jumping ability in some extinct species with foreshortened distal limb segments, such as dromaeosaurid dinosaurs.

Simulations prédictives de la fonction musculo-squelettique et des performances de saut chez un oiseau généralisé Sauter est un comportement commun, mais exigeant, que de nombreux animaux utilisent au cours de leurs activités quotidiennes. Contrairement aux spécialistes du saut tels que les anoures et certains primates, la biomécanique du saut et les facteurs qui influencent la performance restent peu étudiés pour les espèces généralisées qui n’ont pas d’adaptations marquées pour le saut. Les approches de modélisation biomécanique computationnelle offrent un moyen d’aborder cette question de manière rigoureuse et mécaniste. Ici, la théorie du contrôle optimal et la modélisation musculo-squelettique sont intégrées pour générer des simulations prédictives du saut en hauteur maximal chez un petit oiseau terrestre, le tinamou. Un modèle musculo-squelettique tridimensionnel avec 36 actionneurs par patte est utilisé, et une méthode numérique nommée “direct collocation” est employée pour formuler un problème de contrôle optimal rapidement résoluble impliquant les phases de décollage et d’atterrissage. La simulation qui en résulte élève le centre de masse du corps entier à plus du double de sa hauteur debout, et les aspects clés du comportement simulé reproduisent qualitativement les observations empiriques d’autres oiseaux sauteurs. Cependant, les performances quantitatives sont moindres, avec une réduction des forces au sol, des hauteurs de saut et de la puissance musculo-tendineuse. Une manœuvre de contre-mouvement prononcée est utilisée pendant le lancement. Il a été démontré que l’utilisation d’un contre-mouvement est essentielle à l’obtention de hauteurs de saut plus importantes, et il se peut que ce phénomène doive exploiter uniquement des principes physiques pour réussir; l’amplification de la performance musculaire n’est pas nécessairement une raison immédiate de l’utilisation de cette manœuvre. L’augmentation de la force musculaire ou de la vitesse de contraction au-dessus des valeurs nominales améliore grandement la performance de saut et, fait intéressant, a le plus grand effet sur les muscles extenseurs des membres plus distaux (c'est-à-dire ceux de la cheville), ce qui suggère que le membre distal peut être un lien critique pour le comportement de saut. Ces résultats justifient une réévaluation des déductions précédentes de la capacité de sauter chez certaines espèces éteintes avec des segments de membres distaux raccourcis, comme les dinosaures droméosauridés.

Voorspellende simulaties van musculoskeletale functie en springprestaties bij een gegeneraliseerde vogel Springen is een veel voorkomend, maar veeleisend, gedrag dat veel dieren toepassen tijdens hun dagelijkse bezigheden. In tegenstelling tot de springspecialisten zoals de anura en sommige primaten, is de biomechanica van het springen en de factoren die de prestaties beïnvloeden nog weinig bestudeerd voor algemene soorten die geen uitgesproken adaptaties voor het springen hebben. Computationele biomechanische modelbenaderingen bieden een manier om dit op een rigoureuze, mechanistische manier aan te pakken. Hier worden optimale controle theorie en musculoskeletale modellering geïntegreerd om voorspellende simulaties te genereren van maximale hoogtesprong bij een kleine grondbewonende vogel, een tinamou. Een driedimensionaal musculoskeletaal model met 36 actuatoren per poot wordt gebruikt, en directe collocatie wordt toegepast om een snel oplosbaar optimaal controleprobleem te formuleren dat zowel de opstijg-als de landingsfase omvat. De resulterende simulatie verhoogt het lichaamszwaartepunt tot meer dan het dubbele van de stahoogte, en belangrijke aspecten van het gesimuleerde gedrag komen kwalitatief overeen met empirische waarnemingen voor andere springende vogels. De kwantitatieve prestaties zijn echter minder, met verminderde grondkrachten, spronghoogtes en spierpeeskracht. Tijdens de lancering wordt een uitgesproken tegenbewegingsmanoeuvre gebruikt. Aangetoond is dat het gebruik van een tegenbeweging van cruciaal belang is voor het bereiken van grotere spronghoogten, en dit fenomeen hoeft alleen op fysische principes te berusten om succesvol te zijn; versterking van de spierprestaties hoeft niet noodzakelijk een proximate reden te zijn voor het gebruik van deze manoeuvre. Het verhogen van de spierkracht of van de contractiesnelheid boven de nominale waarden verbetert de sprongprestatie aanzienlijk, en heeft interessant genoeg het grootste effect op de meer distale extensoren van de ledematen (d.w.z. die van de enkel), wat suggereert dat de distale ledematen een kritieke schakel kunnen zijn voor het springgedrag. Deze resultaten rechtvaardigen een herevaluatie van eerdere conclusies over springvermogen bij sommige uitgestorven soorten met voorgekorte distale ledematen, zoals dromaeosauride dinosauriërs.

Prädiktive Simulationen der muskuloskelettalen Funktion und Sprungleistung bei einem generalisierten Vogel Springen ist ein übliches jedoch anstrengendes Verhalten, das viele Tiere bei ihren täglichen Aktivitäten einsetzen. Im Gegensatz zu Springspezialisten, wie Fröschen und einigen Primaten, sind bei allgemeinen Arten, welche keine ausgeprägten Anpassung für Sprungverhalten aufweisen, die Biomechanik beim Springen und die Faktoren, welche die Leistungsfähigkeit beeinflussen, noch wenig untersucht. Computergestützte biomechanische Modellierungsverfahren bieten hier eine Möglichkeit, dies in einer gründlichen, mechanistischen Weise anzugehen. In dieser Arbeit werden die optimale Steuerungstheorie und Muskel-Skelett-Modellierung zusammen eingesetzt, um die maximale Sprunghöhe eines kleinen bodenlebenden Vogels, eines Perlsteisshuhns, zu simulieren und zu prognostizieren. Es wird ein dreidimensionales Muskel-Skelett-Modell mit 36 Aktuatoren pro Bein verwendet, und durch direkte Kollokation wird ein schnell lösbares optimales Steuerungsproblem formuliert, das sowohl die Abstoss- als auch die Landephase umfasst. Die daraus folgende Simulation bringt den Ganzkörperschwerpunkt auf mehr als das Doppelte seiner Standhöhe und entscheidende Aspekte des simulierten Verhaltens entsprechen qualitativ empirischen Beobachtungen für andere springende Vögel. Allerdings ist die quantitative Leistungsfähigkeit geringer, mit reduzierten Bodenkräften, Sprunghöhen und Muskel-Sehnen-Kräften. Beim Abstossen wird ein ausgeprägtes Gegenbewegungsmanöver durchgeführt. Die Durchführung einer Gegenbewegung ist nachweislich entscheidend für das Erreichen grösserer Sprunghöhen, wobei dieses Phänomen möglicherweise nur physikalische Prinzipien auszuschöpfen braucht, um erfolgreich zu sein. Die Verstärkung der Muskelleistung ist daher möglicherweise nicht zwingend ein unmittelbarer Grund für die Verwendung dieses Manövers. Eine Erhöhung der Muskelkraft oder der Kontraktionsgeschwindigkeit über die Nominalwerte hinaus führt zu einer erheblichen Zunahme der Sprungleistung und hat interessanterweise den grössten Effekt bei den weiter distal gelegenen Streckmuskeln der Beine (d.h. bei denjenigen des Sprunggelenks), was darauf hindeutet, dass die distale Gliedmasse ein entscheidendes Element für das Sprungverhalten sein könnte. Diese Ergebnisse geben Anlass zur Überprüfung früherer Schlussfolgerungen hinsichtlich der Sprungfähigkeit einiger ausgestorbener Arten mit verkürzten distalen Gliedmassen, wie beispielsweise bei dromaeosauriden Dinosauriern.

Reporting checklist for verification and validation of finite element analysis in orthopedic and trauma biomechanics

Finite element analysis (FEA) has become a fundamental tool for biomechanical investigations in the last decades. Despite several existing initiatives and guidelines for reporting on research methods and results, there are still numerous issues that arise when using computational models in biomechanical investigations. According to our knowledge, these problems and controversies lie mainly in the verification and validation (V&V) process as well as in the set-up and evaluation of FEA. This work aims to introduce a checklist including a report form defining recommendations for FEA in the field of Orthopedic and Trauma (O&T) biomechanics. Therefore, a checklist was elaborated which summarizes and explains the crucial methodologies for the V&V process. In addition, a report form has been developed which contains the most important steps for reporting future FEA. An example of the report form is shown, and a template is provided, which can be used as a uniform basis for future documentation. The future application of the presented report form will show whether serious errors in biomechanical investigations using FEA can be minimized by this checklist. Finally, the credibility of the FEA in the clinical area and the scientific exchange in the community regarding reproducibility and exchangeability can be improved.

Influence of arch support heights on the internal foot mechanics of flatfoot during walking: A muscle-driven finite element analysis

BACKGROUND

Different arch support heights of the customized foot orthosis could produce different effects on the internal biomechanics of the foot. However, quantitative evidence is scarce. Therefore, we aimed to investigate and quantify the influence of arch support heights on the internal foot biomechanics during walking stance.

METHODS

We reconstructed a foot finite element model from a volunteer with flexible flatfoot. The model enabled a three-dimensional representation of the plantar fascia and its interactions with surrounding osteotendinous structures. The volunteer walked in foot orthosis with different arch heights (low, neutral, and high). Muscle forces during gaits were calculated by a multibody model and used to drive a foot finite element model. The foot contact pressures and plantar fascia strains in different regions were compared among the insole conditions at the first and second vertical ground reaction force (VGRF) peak and VGRF valley instants.

RESULTS

The results indicated that peak foot pressures decreased in balanced standing and second VGRF as the arch support height increased. However, peak midfoot pressures increased during all simulated instants. Meanwhile, high arch support decreased the plantar fascia loading by 5%-15.4% in proximal regions but increased in the middle and distal regions.

CONCLUSION

Although arch support could generally decrease the plantar foot pressure and plantar fascia loading, the excessive arch height may induce high midfoot pressure and loadings at the central portion of the plantar fascia. The consideration of fascia-soft tissue interaction in modeling could improve the prediction of plantar fascia strains towards design optimization for orthoses.

Finite Element Simulations of the ID Venous System to Treat Venous Compression Disorders: From Model Validation to Realistic Implant Prediction

The ID Venous System is an innovative device proposed by ID NEST MEDICAL to treat venous compression disorders that involve bifurcations, such as the May-Thurner syndrome. The system consists of two components, ID Cav and ID Branch, combined through a specific connection that prevents the migration acting locally on the pathological region, thereby preserving the surrounding healthy tissues. Preliminary trials are required to ensure the safety and efficacy of the device, including numerical simulations. In-silico models are intended to corroborate experimental data, providing additional local information not acquirable by other means. The present work outlines the finite element model implementation and illustrates a sequential validation process, involving seven tests of increasing complexity to assess the impact of each numerical uncertainty separately. Following the standard ASME V&V40, the computational results were compared with experimental data in terms of force-displacement curves and deformed configurations, testing the model reliability for the intended context of use (differences < 10%). The deployment in a realistic geometry confirmed the feasibility of the implant procedure, without risk of rupture or plasticity of the components, highlighting the potential of the present technology.

Using inverse finite element analysis to identify spinal tissue behaviour in situ

In computational modelling of musculoskeletal applications, one of the critical aspects is ensuring that a model can capture intrinsic population variability and not only representative of a "mean" individual. Developing and calibrating models with this aspect in mind is key for the credibility of a modelling methodology. This often requires calibration of complex models with respect to 3D experiments and measurements on a range of specimens or patients. Most Finite Element (FE) software's do not have such a capacity embedded in their core tools. This paper presents a versatile interface between Finite Element (FE) software and optimisation tools, enabling calibration of a group of FE models on a range of experimental data. It is provided as a Python toolbox which has been fully tested and verified on Windows platforms. The toolbox is tested in three case studies involving in vitro testing of spinal tissues.

Biomechanical Investigation of the Stomach Following Different Bariatric Surgery Approaches

Background: The stomach is a hollow organ of the gastrointestinal tract, on which bariatric surgery (BS) is performed for the treatment of obesity. Even though BS is the most effective treatment for severe obesity, drawbacks and complications are still present because the intervention design is largely based on the surgeon’s expertise and intraoperative decisions. Bioengineering methods can be exploited to develop computational tools for more rational presurgical design and planning of the intervention. Methods: A computational mechanical model of the stomach was developed, considering the actual complexity of the biological structure, as the nonhomogeneous and multilayered configuration of the gastric wall. Mechanical behavior was characterized by means of an anisotropic visco-hyperelastic constitutive formulation of fiber-reinforced conformation, nonlinear elastic response, and time-dependent behavior, which assume the typical features of gastric wall mechanics. Model applications allowed for an analysis of the influence of BS techniques on stomach mechanical functionality through different computational analyses. Results: Computational results showed that laparoscopic sleeve gastrectomy and endoscopic sleeve gastroplasty drastically alter stomach capacity and stiffness, while laparoscopic adjustable gastric banding modestly affects stomach stiffness and capacity. Moreover, the mean elongation strain values, which are correlated to the mechanical stimulation of gastric receptors, were elevated in laparoscopic adjustable gastric banding compared to other procedures. Conclusions: The investigation of stomach mechanical response through computational models provides information on different topics such as stomach capacity and stiffness and the mechanical stimulation of gastric receptors, which interact with the brain to control satiety. These data can provide reliable support to surgeons in the presurgical decision-making process.

Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis

There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.

Understanding the Patterns of Deformity of Wrist Fractures Using Computer Analysis

Computer modeling of the wrist has followed other fields in the search for descriptive methods to understand the biomechanics of injury. Using patient-specific 3D computer models, we may better understand the biomechanics of wrist fractures in order to plan better care. We may better estimate fracture morphology and stability and evaluate surgical indications, design more adequate or effective surgical approaches and develop novel methods of therapy. The purpose of this review is to question the actual advances made in the understanding of wrist fractures using computer models.

Uncertainty quantification reveals the physical constraints on pumping by peristaltic hearts

Most biological functional systems are complex, and this complexity is a fundamental driver of diversity. Because input parameters interact in complex ways, a holistic understanding of functional systems is key to understanding how natural selection produces diversity. We present uncertainty quantification (UQ) as a quantitative analysis tool on computational models to study the interplay of complex systems and diversity. We investigate peristaltic pumping in a racetrack circulatory system using a computational model and analyse the impact of three input parameters (Womersley number, compression frequency, compression ratio) on flow and the energetic costs of circulation. We employed two models of peristalsis (one that allows elastic interactions between the heart tube and fluid and one that does not), to investigate the role of elastic interactions on model output. A computationally cheaper surrogate of the input parameter space was created with generalized polynomial chaos expansion to save computational resources. Sobol indices were then calculated based on the generalized polynomial chaos expansion and model output. We found that all flow metrics were highly sensitive to changes in compression ratio and insensitive to Womersley number and compression frequency, consistent across models of peristalsis. Elastic interactions changed the patterns of parameter sensitivity for energetic costs between the two models, revealing that elastic interactions are probably a key physical metric of peristalsis. The UQ analysis created two hypotheses regarding diversity: favouring high flow rates (where compression ratio is large and highly conserved) and minimizing energetic costs (which avoids combinations of high compression ratios, high frequencies and low Womersley numbers).

Finite Element Analysis of Generalized Ligament Laxity on the Deterioration of Hallux Valgus Deformity (Bunion)

Hallux valgus is a common foot problem affecting nearly one in every four adults. Generalized ligament laxity was proposed as the intrinsic cause or risk factor toward the development of the deformity which was difficult to be investigated by cohort clinical trials. Herein, we aimed to evaluate the isolated influence of generalized ligament laxity on the deterioration using computer simulation (finite element analysis). We reconstructed a computational foot model from a mild hallux valgus participant and conducted a gait analysis to drive the simulation of walking. Through parametric analysis, the stiffness of the ligaments was impoverished at different degrees to resemble different levels of generalized ligament laxity. Our simulation study reported that generalized ligament laxity deteriorated hallux valgus by impairing the load-bearing capacity of the first metatarsal, inducing higher deforming force, moment and malalignment at the first metatarsophalangeal joint. Besides, the deforming moment formed a deteriorating vicious cycle between hallux valgus and forefoot abduction and may result in secondary foot problems, such as flatfoot. However, the metatarsocuneiform joint did not show a worsening trend possibly due to the overriding forefoot abduction. Controlling the deforming load shall be prioritized over the correction of angles to mitigate deterioration or recurrence after surgery.

Computational Biomechanics: In-Silico Tools for the Investigation of Surgical Procedures and Devices

Biomechanical investigations of surgical procedures and devices are usually developed by means of human or animal models. The exploitation of computational methods and tools can reduce, refine, and replace (3R) the animal experimentations for scientific purposes and for pre-clinical research. The computational model of a biological structure characterizes both its geometrical conformation and the mechanical behavior of its building tissues. Model development requires coupled experimental and computational activities. Medical images and anthropometric information provide the geometrical definition of the computational model. Histological investigations and mechanical tests on tissue samples allow for characterizing biological tissues' mechanical response by means of constitutive models. The assessment of computational model reliability requires comparing model results and data from further experimentations. Computational methods allow for the in-silico analysis of surgical procedures and devices' functionality considering many different influencing variables, the experimental investigation of which should be extremely expensive and time consuming. Furthermore, computational methods provide information that experimental methods barely supply, as the strain and the stress fields that regulate important mechano-biological phenomena. In this work, general notes about the development of biomechanical tools are proposed, together with specific applications to different fields, as dental implantology and bariatric surgery.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Overall Strategy

Recent explorations of knee biomechanics have benefited from computational modeling, specifically leveraging advancements in finite element analysis and rigid body dynamics of joint and tissue mechanics. A large number of models have emerged with different levels of fidelity in anatomical and mechanical representation. Adapted modeling and simulation processes vary widely, based on justifiable choices in relation to anticipated use of the model. However, there are situations where modelers' decisions seem to be subjective, arbitrary, and difficult to rationalize. Regardless of the basis, these decisions form the "art" of modeling, which impact the conclusions of simulation-based studies on knee function. These decisions may also hinder the reproducibility of models and simulations, impeding their broader use in areas such as clinical decision making and personalized medicine. This document summarizes an ongoing project that aims to capture the modeling and simulation workflow in its entirety-operation procedures, deviations, models, by-products of modeling, simulation results, and comparative evaluations of case studies and applications. The ultimate goal of the project is to delineate the art of a cohort of knee modeling teams through a publicly accessible, transparent approach and begin to unravel the complex array of factors that may lead to a lack of reproducibility. This manuscript outlines our approach along with progress made so far. Potential implications on reproducibility, on science, engineering, and training of modeling and simulation, on modeling standards, and on regulatory affairs are also noted.

A Combined Geometric Morphometric and Discrete Element Modeling Approach for Hip Cartilage Contact Mechanics

Finite element analysis (FEA) provides the current reference standard for numerical simulation of hip cartilage contact mechanics. Unfortunately, the development of subject-specific FEA models is a laborious process. Owed to its simplicity, Discrete Element Analysis (DEA) provides an attractive alternative to FEA. Advancements in computational morphometrics, specifically statistical shape modeling (SSM), provide the opportunity to predict cartilage anatomy without image segmentation, which could be integrated with DEA to provide an efficient platform to predict cartilage contact stresses in large populations. The objective of this study was, first, to validate linear and non-linear DEA against a previously validated FEA model and, second, to present and evaluate the applicability of a novel population-averaged cartilage geometry prediction method against previously used methods to estimate cartilage anatomy. The population-averaged method is based on average cartilage thickness maps and therefore allows for a more accurate and individualized cartilage geometry estimation when combined with SSM. The root mean squared error of the population-averaged cartilage geometry predicted by SSM as compared to the manually segmented cartilage geometry was 0.31 ± 0.08 mm. Identical boundary and loading conditions were applied to the DEA and FEA models. Predicted DEA stress distribution patterns and magnitude of peak stresses were in better agreement with FEA for the novel cartilage anatomy prediction method as compared to commonly used parametric methods based on the estimation of acetabular and femoral head radius. Still, contact stress was overestimated and contact area was underestimated for all cartilage anatomy prediction methods. Linear and non-linear DEA methods differed mainly in peak stress results with the non-linear definition being more sensitive to detection of high peak stresses. In conclusion, DEA in combination with the novel population-averaged cartilage anatomy prediction method provided accurate predictions while offering an efficient platform to conduct population-wide analyses of hip contact mechanics.

Finite element analysis of the rotator cuff: A systematic review

BACKGROUND

Finite element modeling serves as a promising tool for investigating underlying rotator cuff biomechanics and pathology. However, there are currently no concrete guidelines for reporting in finite element model studies. This has compromised the reliability, validity, and reproducibility of literature due to omission of pertinent items within publications. Recently a Finite Element Model Grading Procedure has been proposed as a reporting guideline for model developers. The aim of this study was to conduct a systematic review of rotator cuff focused finite element models and characterize the reporting quality of those articles.

METHODS

A comprehensive literature search was performed in PubMed, Web of Science, and Embase to find relevant articles. Each article was graded and given a reporting quality ranking based on a score generated from the Finite Element Model Grading Procedure.

FINDINGS

We found that only 5/22 articles had scores of 75% or higher and fell within the "exceptional" reporting quality range. Most of the articles (16/22) fell within the "good" reporting quality range with scores between 50% and 75%. However, 9/16 articles within the "good" reporting quality range had scores below 60%.

INTERPRETATION

This study indicates that improved guidelines and standards for good reporting practices must be made in the field of finite element modeling. Furthermore, it supports the use of the Finite Element Model Grading Procedure as an objective method for evaluating the quality of finite element model reporting in the literature.