Comparison of side-cutting maneuvers versus low impact baseball swing on knee ligament loading in adolescent populations

BACKGROUND

High impact sports are associated with an increased incidence rate for knee ligament injuries, specifically pertaining to the anterior cruciate ligament and medial collateral ligament. What is less clear is (i) the extent to which high impact activities preferentially load the anterior cruciate ligament versus the medial collateral ligament, and (ii) whether both ligaments experience similar stretch ratios during high loading scenarios. Therefore, the goal of this project was to assess how different loading conditions experienced through more at-risk sporting maneuvers influence the relative displacements of the anterior cruciate ligament and medial collateral ligament. The focus of the study was on adolescent patients - a group that has largely been overlooked when studying knee ligament biomechanics.

METHODS

Through kinetic knee data obtained through motion capture experimentation, two different loading conditions (high vs low impact) were applied to 22 specimen-specific adolescent finite element knee models to investigate the biomechanical impact various sporting maneuvers place on the knee ligaments.

FINDINGS

The high impact side cutting maneuver resulted in 102% and 47% increases in ligament displacement compared to the low impact baseball swing (p < 0.05) for both the anterior cruciate ligament and medial collateral ligament.

INTERPRETATION

Quantifying biomechanical risks that sporting activities place on adolescent subjects provides physicians with insight into knee ligament vulnerability. More specifically, knowing the risks that various sports place on ligaments helps guide the selection of sports for at-risk patients (especially those who have undergone knee ligament surgery).

Sex disparities in tibia-fibula geometry and density are associated with elevated bone strain in females: A cross-validation study

Females are up to four times more likely to sustain a stress fracture than males. Our previous work, using statistical appearance modeling in combination with the finite element method, suggested that sex-related differences in tibial geometry may increase bone strain in females. The purpose of this study was to cross-validate these findings, by quantifying sex-related differences in tibia-fibula bone geometry, density, and finite element-predicted bone strain in a new cohort of young physically active adults. CT scans of the lower leg were collected for fifteen males (23.3 ± 4.3 years, 1.77 ± 0.09 m, 75.6 ± 10.0 kg) and fifteen females (22.9 ± 3.0 years, 1.67 ± 0.07 m, 60.9 ± 6.7 kg). A statistical appearance model was fit to each participant's tibia and fibula. The average female and male tibia-fibula complex, controlled for isotropic scaling, were then calculated. Bone geometry, density, and finite element-predicted bone strains in running were compared between the average female and male. The new cohort illustrated the same patterns as the cohort from the previous study: the tibial diaphysis of the average female was narrower and had greater cortical bone density. Peak strain and the volume of bone experiencing ≥4000 με were 10 % and 80 % greater, respectively, in the average female when compared to the average male, which was driven by a narrower diaphysis. The sex-related disparities in tibial geometry, density, and bone strain described by our previous model were also observed in this entirely new cohort. Disparities in tibial diaphysis geometry likely contribute to the elevated stress fracture risk observed in females.

Machine-learning defined precision tDCS for improving cognitive function

BACKGROUND

Transcranial direct current stimulation (tDCS) paired with cognitive training (CT) is widely investigated as a therapeutic tool to enhance cognitive function in older adults with and without neurodegenerative disease. Prior research demonstrates that the level of benefit from tDCS paired with CT varies from person to person, likely due to individual differences in neuroanatomical structure.

OBJECTIVE

The current study aims to develop a method to objectively optimize and personalize current dosage to maximize the functional gains of non-invasive brain stimulation.

METHODS

A support vector machine (SVM) model was trained to predict treatment response based on computational models of current density in a sample dataset (n = 14). Feature weights of the deployed SVM were used in a weighted Gaussian Mixture Model (GMM) to maximize the likelihood of converting tDCS non-responders to responders by finding the most optimum electrode montage and applied current intensity (optimized models).

RESULTS

Current distributions optimized by the proposed SVM-GMM model demonstrated 93% voxel-wise coherence within target brain regions between the originally non-responders and responders. The optimized current distribution in original non-responders was 3.38 standard deviations closer to the current dose of responders compared to the pre-optimized models. Optimized models also achieved an average treatment response likelihood and normalized mutual information of 99.993% and 91.21%, respectively. Following tDCS dose optimization, the SVM model successfully predicted all tDCS non-responders with optimized doses as responders.

CONCLUSIONS

The results of this study serve as a foundation for a custom dose optimization strategy towards precision medicine in tDCS to improve outcomes in cognitive decline remediation for older adults.

A robust integrated control design for weight-on-bit fluctuation suppression in drilling process subject to fractured formations

Fractured formations lead to insufficient contact or loss of contact between the drill bit and rocks, which subsequently causes notable fluctuation in weight-on-bit due to flexible drill strings. This paper presents a robust integrated control design to suppress weight-on-bit fluctuation. First, a finite element drill-string longitudinal model is employed as the basis of control design, which has been verified using actual run data. The integrated controller contains a proportional-integral (PI) controller and a dynamic output feedback controller, dealing with the system's first-order and high-order dynamics. The dynamic output feedback controller only needs to focus on the system's high-order dynamics since the controller synthesis is based on the pre-designed PI controller. The controller decreases the flexible mode amplitude from disturbance to the output channel through H-infinity loop shaping, making the closed-loop system less flexible. The numerical results demonstrate that the developed approach can effectively suppress weight-on-bit fluctuation due to drill-string flexibility when drilling in fractured formations.

Finite element analysis of conductive hearing loss caused by fixation and detachment of ligament and tendon in the middle ear

BACKGROUND AND OBJECTIVE

The fixation of ligament and tendon of the middle ear often occurs after chronic otitis media surgery. However, there are relatively few studies on the effect of ligament and tendon on sound transmission in the human middle ear. Here, the finite element model and lumped parameter model are used to study the effect of ligament and tendon fixation and detachment on sound transmission in human ear.

METHODS

In this paper, the finite element model including the external auditory canal, middle ear and simplified inner ear is used to calculate and compare the middle ear frequency response of the normal and tympanosclerosis under pure tone stimulation. In addition, the lumped parametric model is taken into account to illustrate the effect of ligament and tendon stiffness on the human ear transmission system.

RESULTS

The results indicate that the motion of the tympanic membrane and stapes is reduced by ligament and tendon fixation. Although ligament and tendon detachment have a limited effect in the piston-motion direction, the stability of motion in the plane perpendicular to the piston-motion direction is significantly reduced. Most significantly, the ligament and tendon fixation cause a hearing effect of about 18 dB, which is greater in the plane perpendicular to the piston-motion direction after ligament and tendon detachment than in the piston-motion direction.

CONCLUSIONS

In this study, the calculation accuracy of the lumped parameter and the finite element model is studied, and the effect of ligament and tendon on hearing loss is further explored through the finite element model with high calculation accuracy, which is helpful to understand the role of ligament and tendon in the sound transmission mechanism of the human middle ear. The study of ligament and tendon on conductive hearing loss provides a reference for clinical treatment of tympanosclerosis.

Adaptation of Fibril-Reinforced Poroviscoelastic Properties in Rabbit Collateral Ligaments 8 Weeks After Anterior Cruciate Ligament Transection

Ligaments of the knee provide stability and prevent excessive motions of the joint. Rupture of the anterior cruciate ligament (ACL), a common sports injury, results in an altered loading environment for other tissues in the joint, likely leading to their mechanical adaptation. In the collateral ligaments, the patterns and mechanisms of biomechanical adaptation following ACL transection (ACLT) remain unknown. We aimed to characterize the adaptation of elastic and viscoelastic properties of the lateral and medial collateral ligaments eight weeks after ACLT. Unilateral ACLT was performed in six rabbits, and collateral ligaments were harvested from transected and contralateral knee joints after eight weeks, and from an intact control group (eight knees from four animals). The cross-sectional areas were measured with micro-computed tomography. Stepwise tensile stress-relaxation testing was conducted up to 6% final strain, and the elastic and viscoelastic properties were characterized with a fibril-reinforced poroviscoelastic material model. We found that the cross-sectional area of the collateral ligaments in the ACL transected knees increased, the nonlinear elastic collagen network modulus of the LCL decreased, and the amount of fast relaxation in the MCL decreased. Our results indicate that rupture of the ACL leads to an early adaptation of the elastic and viscoelastic properties of the collagen fibrillar network in the collateral ligaments. These adaptations may be important to consider when evaluating whole knee joint mechanics after ACL rupture, and the results aid in understanding the consequences of ACL rupture on other tissues.

Spinal Cord Boundary Conditions Affect Brain Tissue Strains in Impact Simulations

Brain and spinal cord injuries have devastating consequences on quality of life but are challenging to assess experimentally due to the traumatic nature of such injuries. Finite element human body models (HBM) have been developed to investigate injury but are limited by a lack of biofidelic spinal cord implementation. In many HBM, brain models terminate with a fixed boundary condition at the brain stem. The goals of this study were to implement a comprehensive representation of the spinal cord into a contemporary head and neck HBM, and quantify the effect of the spinal cord on brain deformation during simulated impacts. Spinal cord tissue geometries were developed, based on 3D medical imaging and literature data, meshed, and implemented into the GHBMC 50th percentile male model. The model was evaluated in frontal, lateral, rear, and oblique impact conditions, and the resulting maximum principal strains in the brain tissue were compared, with and without the spinal cord. A new cumulative strain curve metric was proposed to quantify brain strain distribution. Presence of the spinal cord increased brain tissue strains in all simulated cases, owing to a more compliant boundary condition, highlighting the importance of the spinal cord to assess brain response during impact.

Development, validation and a case study: The female finite element head model (FeFEHM)

BACKGROUND AND OBJECTIVE

Traumatic brain injuries are one of the leading causes of death and disability in the world. To better understand the interactions and forces applied in different constituents of the human head, several finite element head models have been developed throughout the years, for offering a good cost-effective and ethical approach compared to experimental tests. Once validated, the female finite element head model (FeFEHM) will allow a better understanding of injury mechanisms resulting in neuronal damage, which can later evolve into neurodegenerative diseases.

METHODS

This work encompasses the approached methodology starting from medical images and finite element modelling until the validation process using novel experimental data of brain displacements conducted on human cadavers. The material modelling of the brain is performed using an age-specific characterization of the brain using microindentation at dynamic rates and under large deformation, with a similar age to the patient used to model the FeFEHM.

RESULTS

The numerical displacement curves are in good accordance with the experimental data, displaying similar peak times and values, in all three anatomical planes. The case study result shows a similarity between the pressure fields of the FeFEHM compared to another model, highlighting the future potential of the model.

CONCLUSIONS

The initial objective was met, and a new female finite element head model has been developed with biofidelic brain motion. This model will be used for the assessment of repetitive impact scenarios and its repercussions on the female brain.

Numerical model characterization of the sound transmission mechanism in the tympanic membrane from a high-speed digital holographic experiment in transient regime

A methodology for the development of a finite element numerical model of the tympanic membrane (TM) based on experiments carried out in the time domain on a cadaveric human temporal bone is presented. Using a high-speed digital holographic (HDH) system, acoustically-induced transient displacements of the TM surface are obtained. The procedure is capable to generate and validate the finite element model of the TM by numerical and experimental data correlation. Reverse engineering approach is used to identify key material parameters that define the mechanical response of the TM. Finally, modal numerical simulations of the specimen are performed. Results show the feasibility of the methodology to obtain an accurate model of a specific specimen and to help interpret its behaviour with additional numerical simulations. STATEMENT OF SIGNIFICANCE: Improving knowledge of the dynamic behavior of the tympanic membrane is key to understanding the sound transmission system in human hearing and advance in the treatment of its pathologies. Recently we acquired a new tool to carry out experiments in transient regime by means of digital laser holography, capable of providing a large amount of information in a controlled transient test. In this work, these data are used to develop a methodology that generates a numerical model of the tympanic membrane based on numerical-experimental correlations. It is important to be able to develop models that fit specific patients. In this work, additional modal simulations are also presented that, in addition to validating the results, provide more information on the specimen.

Numerical investigation of the biomechanical effects of orthodontic intermaxillary elastics on the temporomandibular joint

Temporomandibular joint disorder (TMD) often coincides with malocclusion, and in some cases TMDs are reported after orthodontic treatment. Intermaxillary elastics (also known as orthodontic elastics, OE) are a common way to apply force during orthodontic treatment, and they might cause mechanical effects on the temporomandibular joint (TMJ), thereby lead to joint remodeling. It is still a controversial topic whether the adapted remodeling of the TMJ or of the alveolar bone is the main cause for the alteration of occlusion after treatment with OEs. It was the aim of this study to analyze whether variations of OEs would develop harmful effects on the healthy TMJ. A TMJ model with a masticatory system based on Hill-type muscle actuators was established. Mouth opening and closure with and without OEs were simulated, and maximum principal stresses in the disc and condylar cartilage as well as the displacement of the mandible were analyzed. We found no considerably difference in the mandibular movement without and with symmetrical OEs during mouth opening and closing. At full mouth opening, stresses in the disc and condylar cartilage of some models with OEs were much smaller than without OEs, but we did not find consistency in the results from the left and right sides of the same model (e.g. the lowest compressive stress on the left side of disc from the model with Class II OEs is much smaller than without OEs, -17.3 MPa compared with -28.2, while on the right side, there was no obvious difference). Hence, we could not conclude that OEs would develop deleterious effects on the healthy TMJ.

Capacitive stimulation-sensing system for instrumented bone implants: Finite element model to predict the electric stimuli delivered to the interface

BACKGROUND

Prevalence of orthopaedic replacements are increasing around the world. The main cause of revision remains associated to the interface loosening. In this work, a computational study using the Finite element method was developed to predict the electric field stimuli delivered to trabecular bone structures, as well as to predict the sensing ability to detect different bone-implant interface scenarios.

METHODS

Three finite element models were developed: two simplified models, including a Gyroid TMP structure, and a realistic model based on microCT scan of a trabecular bone from sheep vertebra. Simulations were performed using a co-surface capacitive technology for stimulating and sensing bone-implant interfaces. Different fixation scenarios were considered, namely by establishing bone-stimulator gap sizes up to 1 mm (from fixation to massive loosening scenario). Electrodes were excited with sinusoidal and square electric signals up to 10V voltage and 64kHz frequency.

RESULTS

Simplification of bone geometry resulted in significant electric stimuli differences compared to the realistic bone geometry. Realistic modelling allowed to observe that, in the fixation scenario, the electric field stimuli decreased 85% from the sensor interface to a parallel plane 2 mm apart from such interface. A significant influence of the bone-stimulator distance on the electric stimuli was found: the electric stimuli magnitudes varied in the range between 0.38 V/mm (fixation scenario) and 4.8 mV/mm (massive loosening scenario) for voltages up to 10V. Strong frequency-dependent behaviours were also observed in the electric stimuli: their magnitudes can reach 106-fold decreases when the excitation frequency is decreased from 32 kHz to 14 Hz CONCLUSION: This study points out the inability of our two simplified models to predict the electric stimulation provided to different bone-implant interface scenarios. Results highlight that co-surface stimulators can deliver osteogenic electric stimuli along trabecular bone structures, ensuring low electric power excitations. Moreover, realistic models strongly enhance the sensing predictability of the bone-implant fixation states. These new and significant evidences provide a strong support to integrate co-surface capacitive into bioelectronic implants for both therapeutic and sensing operations.

A deep dive into the static force transmission of the human masticatory system and its biomechanical effects on the temporomandibular joint

OBJECTIVE

This study aims to investigate the biomechanical behavior and reveal the force transmission patterns of the human masticatory system through advanced three-dimensional finite element (FE) models.

METHODS

The FE model was constructed according to the medical images of a healthy male adult. It contains full skull structures, detailed temporomandibular joints (TMJs) with discs, complete dentitions, masticatory muscles, and related ligaments. Several static bite scenarios were simulated to demonstrate the effects of bite positions and muscle force recruitments on the force transmission patterns.

RESULTS

Molar occlusal surfaces are the primary force transmission region for clenching. Sensitivity analysis demonstrated that the stiffness of the bite substance would not alter the force transmission patterns but could affect the maximum contact stresses on the discs and the occlusal surfaces. During the unilateral clenching tasks, the high-stress region on the discal surfaces shifted ipsilaterally. The presence or absence of the molar cushions would significantly affect the biomechanical response of the masticatory system.

SIGNIFICANCE

FE analysis is an effective way of investigating biomechanical responses involving complicated interactions. Enriching the static analysis of the masticatory system with a detailed model can help understand better how the forces were transmitted and the significance of TMJs during the clenching process.

Comparative Study of Miniplate vs Reconstruction Plate in the Management of Bilateral Parasymphysis Mandible Fracture: FEM Analysis

Aim and Objective

The present study compared the stability of fracture fragments in the management of bilateral parasymphysis mandible fracture with Miniplate fixation and Reconstruction plate fixation using finite element analysis.

Material and Method

3D FE Mandible model was created using CT scanner. Two bilateral parasymphysis mandible fracture models were created. Model 1 was fixed with Miniplates, and Model 2 was fixed with Reconstruction plate. Loading forces of 120 N at molar region and 62.5 N at incisor region were applied. These two models were imported to ANSY'S Workbench Software.

Result

Miniplate fixation model showed comparatively reduced gap between fragments than Reconstruction plate. But the gap values of both the models were within the physiologic limit of healing under this specific loading. Analytically Miniplates were superior to Reconstruction plate in the management of bilateral parasymphysis fracture.

Conclusion

Analytically Miniplates are superior to Reconstruction plate in the management of bilateral parasymphysis fracture. As the masticatory forces were reduced during fracture healing period, both fixations provide satisfactory healing. So both Miniplate and Reconstruction plate can be considered as fixation method for bilateral parasymphysis mandible fracture.

Biomechanical Effect of Osteoporosis on Adjacent Segments After Anterior Cervical Corpectomy and Fusion

BACKGROUND

Adjacent segmental degeneration (ASD) is one common long-term complication of anterior cervical corpectomy and fusion (ACCF), and osteoporosis is one basic disease in the elderly. After ACCF, patients may experience osteoporosis with age. However, the influence of osteoporosis on ASD remains unclear. The purpose of this study was to determine whether osteoporosis could affect the development of ASD following ACCF.

METHODS

Three finite element models of the cervical spine, including 1 normal model, 1 ACCF model, and 1 ACCF with osteoporosis model, were constructed. ACCF was simulated at the C4-C6 level. A 73.6 N follower load and a 1 Nm moment were imposed on the normal model, and the same follower load together with an adjusted moment was applied to the ACCF model and the ACCF with osteoporosis model, to simulate movement in each direction. The range of motion, intradiscal pressure, shear stress on anulus fibrosus, and facet joint stress at C3-C4 and C6-C7 levels of the models were calculated.

RESULTS

In this study, the normal model was well validated. In flexion, extension, right lateral bending, and right axial rotation, the overall range of motion was 8.92°, 19.7°, 15.37°, and 45.27° in the normal model, and the adjusted moment was 1.4 Nm, 2.7 Nm, 1.1 Nm, and 2.6 Nm in the ACCF model, and 1.3 Nm, 2.5 Nm, 1.1 Nm, and 2.4 Nm in the ACCF with osteoporosis model. Despite of a few exceptions, the maximum values of the outcome measurements were mostly found in the ACCF model, and the minimum values in the normal model. Compared with the ACCF model, most of the outcome measurements were decreased in the ACCF with osteoporosis model.

CONCLUSIONS

Osteoporosis can retard the adverse influence of ACCF on adjacent segments.

Biomechanical study of the male lower urinary tract: Simulation of internal and external sphincters dyssynergia

Urethral sphincter dysfunction is an important cause of stress urinary incontinence (SUI). The most effective treatment is the insertion of an artificial urethral sphincter (AUS), which relies to a large extent on the surgeon's experience. However, there is no quantitative standard for cuff tightness, resulting in frequent postoperative complications. This study aimed to investigate the effect of internal and external sphincter dyssynergia on urodynamic parameters in the lower urinary tract. A geometric model of male lower urinary tract tissue was constructed from collodion slices, accounting for the active behavior of the internal and external sphincters. Normal and dyssynergic internal and external sphincters (active sphincter behavior was individually injured by 25%, 50%, 75%, or 100%) were simulated with fluid-structure interactions and changes in urethral stress, displacement, and urine flow rate were detected. We found that when the internal sphincter was injured by 25%, 50%, 75%, and 100%, urethral stress near the internal sphincter decreased by 8.3%, 15.6%, 24.3%, and 35.7%, respectively. Additionally, when the external sphincter was injured by 25%, 50%, 75%, and 100%, urethral stress near the external sphincter was reduced by 13.3%, 24.3%, 38.6%, and 46.6%, respectively. Internal sphincter injury primarily affects positions near the internal sphincter and prostate, while external sphincter injury affects the area between the prostate and urethral outlet. These data could facilitate the standardized evaluation of internal and external sphincter dysfunction and lead to novel methods of preoperative assessment for AUS surgery.

Investigating the post-yield behavior of mineralized bone fibril arrays using a 3D non-linear finite element unit-cell model

In this study, we propose a 3D non-linear finite element (FE) unit-cell model to investigate the post-yield behavior of mineralized collagen fibril arrays (FAY). We then compare the predictions of the model with recent micro-tensile and micropillar compression tests in both axial and transverse directions. The unit cell consists of mineralized collagen fibrils (MCFs) embedded in an extrafibrillar matrix (EFM), and the FE mesh is equipped with cohesive interactions and a custom plasticity model. The simulation results confirm that MCF plays a dominant role in load bearing prior to yielding under axial tensile loading. Damage was initiated via debonding in shear and progressive sliding at the MCF/EFM interface, and resulted in MCF pull-out until brittle failure. In transverse tensile loading, EFM carried most of the load in pre-yield deformation, and then mixed normal/shear debonding between MCF and EFM began to form, which eventually produced brittle delamination of the two phases. The loading/unloading FE analysis in compression along both axial and transverse directions demonstrated perfect plasticity without any reduction in elastic modulus, i.e., damage due to the interfaces as seen in micropillar compression. Beyond the brittle and ductile nature of the stress-strain curves, in tensile and compressive loading, the simulated post-yield behavior and failure mechanism are in good quantitative agreement with the experimental observations. Our rather simple but efficient unit-cell FE model can reproduce qualitatively and quantitatively the mechanical behavior of bone ECM under tensile and compressive loading along the two main orientations. The model's integration into higher length scales may be useful in describing the macroscopic post-yield and failure behavior of trabecular and cortical bone in greater detail.

Using Finite Element Modeling in Bone Mechanoadaptation

PURPOSE OF THE REVIEW

Bone adapts structure and material properties in response to its mechanical environment, a process called mechanoadpatation. For the past 50 years, finite element modeling has been used to investigate the relationships between bone geometry, material properties, and mechanical loading conditions. This review examines how we use finite element modeling in the context of bone mechanoadpatation.

RECENT FINDINGS

Finite element models estimate complex mechanical stimuli at the tissue and cellular levels, help explain experimental results, and inform the design of loading protocols and prosthetics. FE modeling is a powerful tool to study bone adaptation as it complements experimental approaches. Before using FE models, researchers should determine whether simulation results will provide complementary information to experimental or clinical observations and should establish the level of complexity required. As imaging technics and computational capacity continue increasing, we expect FE models to help in designing treatments of bone pathologies that take advantage of mechanoadaptation of bone.

Failure analysis of a locking compression plate with asymmetric holes and polyaxial screws

Locking compression plates (LCP) with asymmetrical holes and polyaxial screws are effective for treating mid-femoral fractures, but are prone to failure in cases of bone nonunion. To understand the failure mechanism of the LCP, this study assessed the material composition, microhardness, metallography, fractography and biomechanical performance of a retrieved LCP used for treating a bone fracture of AO type 32-A1. For the biomechanical assessment, a finite element surgical model implanted with the intact fixation-plate system was constructed to understand the stresses and structural stiffness on the construct. Also, to avoid positioning screws around the bone fracture, different working lengths of the plate (the distance between the two innermost screws) and screw inclinations (±5°, ±10° and ±15°) were investigated. The fracture site of the retrieved LCP was divided into a narrow side and broad side due to the asymmetrical distribution of holes on the plate. The results indicated that the chemical composition and microhardness of the LCP complied with ASTM standards. The fatigue failure was found to originate on the narrow side of the hole, while the broad side showed overloading characteristics of crack growth. When the screws were inserted away from the region of the bone fracture by increasing the working length, the stress of the fixation-plate system decreased. Regardless of the screw insertion angle, the maximum stress on the LCP always appeared on the narrow side, and there was little change in the structural stiffness. However, angling the screws at -10° resulted in the most even stress distribution on the fixation-plate system. In conclusion, the LCP assessed in this study failed by fatigue fracture due to bone nonunion and stress concentration. The narrow side of the LCP was vulnerable to failure and needs to be strengthened. When treating an AO type 32-A1 fracture using an LCP with asymmetrical holes and polyaxial screws, inserting the screws at -10° may reduce the risk of implant failure and positing screws around the fractured area of the bone should be avoided.

A numerical model for fibril remodeling in articular cartilage

BACKGROUND

Collagen fibrils of articular cartilage have a distinct organization in mature human knee joints. It seems that a mechanobiological process drives the remodeling of newborn collagen fibrils with maturation. Therefore, the goal of the present study was to develop a collagen fibril remodeling algorithm that describes the unique collagen fibril organization in a 3D knee model.

METHOD

A fibril-reinforced, biphasic cartilage model was used with a cuboid and a 3D human knee joint geometries. An isotropic collagen fibril distribution was assigned to the cartilage at the start of the analysis. Each fibril was rotated towards the direction that resulted in a maximum stretch at each time increment of the loading cycle.

RESULTS

The resulting pattern for the collagen fibrils was compared with split line patterns of porcine knee joint cartilage and also data published in the literature. Fibrils on the articular surface had a radial pattern towards the geometrical centroid of the tibial and femoral cartilage. In the tibiofemoral contact regions of superficial zone, fibrils were oriented circumferentially and randomly. In the porcine samples, the split-line patterns were similar to those obtained theoretically. Depth-wise organization of fibril network was characterized by fibrils perpendicular to the subchondral bone in the deeper layers, and fibrils parallel to the surface of cartilage in the superficial zone.

CONCLUSIONS

The maximum stretch criterion, coupled with a biphasic constitutive model, successfully predicted the collagen fibril organization observed in the articular cartilage throughout the depth and on the articular surface.

Patient-specific finite element modeling of scoliotic curve progression using region-specific stress-modulated vertebral growth

PURPOSE

This study describes the creation of patient-specific (PS) osteo-ligamentous finite element (FE) models of the spine, ribcage, and pelvis, simulation of up to three years of region-specific, stress-modulated growth, and validation of simulated curve progression with patient clinical angle measurements.

RESEARCH QUESTION

Does the inclusion of region-specific, stress-modulated vertebral growth, in addition to scaling based on age, weight, skeletal maturity, and spine flexibility allow for clinically accurate scoliotic curve progression prediction in patient-specific FE models of the spine, ribcage, and pelvis?

METHODS

Frontal, lateral, and lateral bending X-Rays of five AIS patients were obtained for approximately three-year timespans. PS-FE models were generated by morphing a normative template FE model with landmark points obtained from patient X-rays at the initial X-ray timepoint. Vertebral growth behavior and response to stress, as well as model material properties were made patient-specific based on several prognostic factors. Spine curvature angles from the PS-FE models were compared to the corresponding X-ray measurements.

RESULTS

Average FE model errors were 6.3 ± 4.6°, 12.2 ± 6.6°, 8.9 ± 7.7°, and 5.3 ± 3.4° for thoracic Cobb, lumbar Cobb, kyphosis, and lordosis angles, respectively. Average error in prediction of vertebral wedging at the apex and adjacent levels was 3.2 ± 2.2°. Vertebral column stress ranged from 0.11 MPa in tension to 0.79 MPa in compression.

CONCLUSION

Integration of region-specific stress-modulated growth, as well as adjustment of growth and material properties based on patient-specific data yielded clinically useful prediction accuracy while maintaining physiological stress magnitudes. This framework can be further developed for PS surgical simulation.

Evaluation of Deep Brain Stimulation (DBS) Lead Biomechanical Interaction with Brain Tissue

The current study aims to examine the effect of material properties on implanted leads used for deep brain stimulation (DBS) using finite element (FE) analysis to investigate brain deformation around an implanted DBS lead in response to daily head accelerations. FE analysis was used to characterize the relative motion of the DBS lead in a suite of fifteen cases sampled from a previously derived kinematic envelope representative of everyday activities describing translational and rotational pulse shape, magnitude, and duration. Load curves were applied to the atlas-based brain model (ABM) with a scaled Haversine acceleration pulse in four directions of rotation: + X, - Y, + Y, and + Z. In addition to the fifteen sampled cases, six experimental cases taken from a previous literature review were also simulated for comparison. The current investigation found that there was very little difference in brain response for the DBS leads with two different material properties. In general, the brain and DBS lead experienced the greatest deformation during rotation about the Z axis for similar load cases. In conclusion, this study showed that there was no significant difference in implanted DBS lead deformation based on lead material properties.

The dependency of fetal left ventricular biomechanics function on myocardium helix angle configuration

The helix angle configuration of the myocardium is understood to contribute to the heart function, as finite element (FE) modeling of postnatal hearts showed that altered configurations affected cardiac function and biomechanics. However, similar investigations have not been done on the fetal heart. To address this, we performed image-based FE simulations of fetal left ventricles (LV) over a range of helix angle configurations, assuming a linear variation of helix angles from epicardium to endocardium. Results showed that helix angles have substantial influence on peak myofiber stress, cardiac stroke work, myocardial deformational burden, and spatial variability of myocardial strain. A good match between LV myocardial strains from FE simulations to those measured from 4D fetal echo images could only be obtained if the transmural variation of helix angle was generally between 110 and 130°, suggesting that this was the physiological range. Experimentally discovered helix angle configurations from the literature were found to produce high peak myofiber stress, high cardiac stroke work, and a low myocardial deformational burden, but did not coincide with configurations that would optimize these characteristics. This may suggest that the fetal development of myocyte orientations depends concurrently on several factors rather than a single factor. We further found that the shape, rather than the size of the LV, determined the manner at which helix angles influenced these characteristics, as this influence changed significantly when the LV shape was varied, but not when a heart was scaled from fetal to adult size while retaining the same shape. This may suggest that biomechanical optimality would be affected during diseases that altered the geometric shape of the LV.

Interfacing finite elements with deep neural operators for fast multiscale modeling of mechanics problems

Multiscale modeling is an effective approach for investigating multiphysics systems with largely disparate size features, where models with different resolutions or heterogeneous descriptions are coupled together for predicting the system's response. The solver with lower fidelity (coarse) is responsible for simulating domains with homogeneous features, whereas the expensive high-fidelity (fine) model describes microscopic features with refined discretization, often making the overall cost prohibitively high, especially for time-dependent problems. In this work, we explore the idea of multiscale modeling with machine learning and employ DeepONet, a neural operator, as an efficient surrogate of the expensive solver. DeepONet is trained offline using data acquired from the fine solver for learning the underlying and possibly unknown fine-scale dynamics. It is then coupled with standard PDE solvers for predicting the multiscale systems with new boundary/initial conditions in the coupling stage. The proposed framework significantly reduces the computational cost of multiscale simulations since the DeepONet inference cost is negligible, facilitating readily the incorporation of a plurality of interface conditions and coupling schemes. We present various benchmarks to assess the accuracy and efficiency, including static and time-dependent problems. We also demonstrate the feasibility of coupling of a continuum model (finite element methods, FEM) with a neural operator, serving as a surrogate of a particle system (Smoothed Particle Hydrodynamics, SPH), for predicting mechanical responses of anisotropic and hyperelastic materials. What makes this approach unique is that a well-trained over-parametrized DeepONet can generalize well and make predictions at a negligible cost.

Factors affecting the numerical response and fracture location of the GHBMC M50 rib in dynamic anterior-posterior loading

Rib fractures are common traumatic injuries, with links to increased morbidity and mortality. Finite element ribs from human body models have struggled to predict the force-displacement response, force and displacement at fracture, and the fracture location for isolated rib tests. In the current study, the sensitivity of a human body model rib with updated anisotropic and asymmetric material models to changes in boundary conditions, material properties, and geometry was investigated systematically to quantify contributions to response. The updated material models using uncalibrated average material properties from literature improved the force-displacement response of the model, whereas the cross-sectional geometry was the only parameter to effect fracture location. The resulting uncalibrated model with improved material models and cross-sectional geometry closely predicted experimental average force-displacement response and fracture location.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.