Specimen width affects vascular tissue integrity for in-vitro characterisation

The preparation of slender specimens for in-vitro tissue characterisation could potentially alter mechanical tissue properties. To investigate this factor, rectangular specimens were prepared from the wall of the porcine aorta for uniaxial tensile loading. Varying strip widths of 16 mm, 8 mm, and 4 mm were achieved by excising zero, one, and three cuts within the specimen along the loading direction, respectively. While specimens loaded along the vessel's circumferential direction acquired consistent tissue properties, the width of test specimens influenced the results of axially loaded tissue; vascular wall stiffness was reduced by approximately 40% in specimens with strips 4 mm wide. In addition, the cross-loading stretch was strongly influenced by specimen strip width, and fiber sliding contributed to the softening of slender tensile specimens, an outcome from finite element analysis of test specimens. We may, therefore, conclude that cutting orthogonal to the main direction of collagen fibers introduces mechanical trauma that weakens slender tensile specimens, compromising the determination of representative mechanical vessel wall properties.

Discrete network models of endothelial cells and their interactions with the substrate

Endothelial cell monolayers line the inner surfaces of blood and lymphatic vessels. They are continuously exposed to different mechanical loads, which may trigger mechanobiological signals and hence play a role in both physiological and pathological processes. Computer-based mechanical models of cells contribute to a better understanding of the relation between cell-scale loads and cues and the mechanical state of the hosting tissue. However, the confluency of the endothelial monolayer complicates these approaches since the intercellular cross-talk needs to be accounted for in addition to the cytoskeletal mechanics of the individual cells themselves. As a consequence, the computational approach must be able to efficiently model a large number of cells and their interaction. Here, we simulate cytoskeletal mechanics by means of molecular dynamics software, generally suitable to deal with large, locally interacting systems. Methods were developed to generate models of single cells and large monolayers with hundreds of cells. The single-cell model was considered for a comparison with experimental data. To this end, we simulated cell interactions with a continuous, deformable substrate, and computationally replicated multistep traction force microscopy experiments on endothelial cells. The results indicate that cell discrete network models are able to capture relevant features of the mechanical behaviour and are thus well-suited to investigate the mechanics of the large cytoskeletal network of individual cells and cell monolayers.

Investigation of hybrid ultrasonic machining process of Nomex honeycomb composite using a toothed disc cutter

Nomex honeycomb composite (NHC) has been increasingly used in the automotive, aerospace and defence applications due to its excellent thermal and mechanical properties. Its complex cellular hexagonal thin-walled configuration along with heterogeneous, soft and brittle nature pose substantial processing defects such as burr formation, tearing of walls, surface roughness, dimensional inaccuracy and low machining quality during conventional machining (CM). These surface defects have a substantial influence on the operating life and functional performance of its sandwiched structural members. Hybrid ultrasonic vibration assisted machining (HUSVAM) technology has been introduced to overcome such limitations. For an in depth research on HUSVAM of NHC using ultrasonic toothed disc cutter (UTDC), three dimensional finite element model was developed and experimental validation was carried out. The numerical simulation and experimental results were found to be in good agreement with one another. The influence of various machining parameters including the ultrasonic tool vibration amplitude (UTVA), feed rate (FR), depth of cut (DoC) and spindle angular speed (SAS) on NHC cutting forces, chip formation and surface quality was investigated using both HUSVAM and CM techniques. Furthermore, stresses, deformation and forces near tool-workpiece interaction along with the surface topography and morphologies were also analysed and compared for various operating conditions. A reduction in the cutting forces was found with the increase of SAS and UTVA (up to 65.47 % and 65.74 %, respectively). Although, the cutting forces were observed to increase by increasing the DoC and FR (up to 159.45 % and 126.33 %, respectively). DoC has a greater impact on the cutting forces among all machining parameters according to 4 levels, 4 factors (L16) orthogonal experiments based on Taguchi method. The results show that the chip formation and machining quality of NHC core can be improved with HUSVAM technique using UTDC. HUSVAM also reduced the cutting forces (up to 73 %) compared to CM. The ultrasonic toothed disc cutter generated large number of burr with very short length, no tearing defects and no uncut fibers as observed from scanning electron microscopy of NHC hexagonal cell structure, walls and triple points. A burr formation of less than 10 % was realized during HUSVAM of NHC cores for Fx≤3N, while it was found up to 20 % if Fx>3N, compared to at least 40 % during CM. The finite element model developed can be used to investigate the influence of HUSVAM on modern difficult to machine materials for enhanced surface quality.

Finite Element Model of the Shoulder with Active Rotator Cuff Muscles: Application to Wheelchair Propulsion

The rotator cuff is prone to injury, remarkably so for manual wheelchair users. To understand its pathomechanisms, finite element models incorporating three-dimensional activated muscles are needed to predict soft tissue strains during given tasks. This study aimed to develop such a model to understand pathomechanisms associated with wheelchair propulsion. We developed an active muscle model associating a passive fiber-reinforced isotropic matrix with an activation law linking calcium ion concentration to tissue tension. This model was first evaluated against known physiological muscle behavior; then used to activate the rotator cuff during a wheelchair propulsion cycle. Here, experimental kinematics and electromyography data was used to drive a shoulder finite element model. Finally, we evaluated the importance of muscle activation by comparing the results of activated and non-activated rotator cuff muscles during both propulsion and isometric contractions. Qualitatively, the muscle constitutive law reasonably reproduced the classical Hill model force-length curve and the behavior of a transversally loaded muscle. During wheelchair propulsion, the deformation and fiber stretch of the supraspinatus muscle-tendon unit pointed towards the possibility for this tendon to develop tendinosis due to the multiaxial loading imposed by the kinematics of propulsion. Finally, differences in local stretch and positions of the lines of action between activated and non-activated models were only observed at activation levels higher than 30%. Our novel finite element model with active muscles is a promising tool for understanding the pathomechanisms of the rotator cuff for various dynamic tasks, especially those with high muscle activation levels.

Biomechanical evaluation of a modified intramedullary nail for the treatment of unstable femoral trochanteric fractures

Background

The Proximal Femoral Nail Antirotation (PFNA) device is the most commonly used implant to fix unstable femoral trochanteric fractures (UFTFs), but it has a relatively high incidence of complications. Due to this factor, the modified intramedullary nail (MIN) was created to treat patients with UFTFs. The aim was to exhibit the MIN and make a comparison with PFNA and InterTAN using biomechanical methods.

Methods

An adult UFTF model was developed using Mimics software. The PFNA, InterTAN nail, and MIN models were drawn referring to the corresponding parameters and installed in the fracture models. Vertical, anteroposterior (AP) bending, and torsion loads of the femoral head were set in advance and loaded onto the fracture models. The value of maximal displacement and von Mises stress was evaluated via finite element analysis (FEA).

Results

The MIN model had smaller values for maximal displacement than that of the PFNA model, and the increase in displacement was less pronounced for the MIN compared to PFNA under increasing vertical loads. For the indicator of von Mises stress, the MIN model showed lower stress compared with the PFNA model in vertical loads ranging from 300 N to 2100 N. Except for the maximal stress at implants under AP bending loads, the MIN demonstrated the most superior biomechanical properties under AP bending and torsional loads.

Conclusion

The MIN offered obvious advantages in terms of mechanical stability and stress distribution among the three studied implants, providing a promising implant option for patients with UFTFs.

Failure mechanism and bearing force of CFRP strengthened square hollow section under compressive load

Carbon fibre-reinforced polymer (CFRP) plates can efficiently repair or enhance the mechanical properties of the square hollow section. However, the loading end of such a CFRP-strengthened member is prone to local bearing failure under compressive load. Given this limitation, an innovative CFRP-plate-strengthened square hollow section composite member (CFRP-SHSCM) was raised, and the thick-walled section was welded on both ends of the thin-walled steel column. The mechanical properties of CFRP-SHSCMs were investigated through parameter finite element (FE) analysis, focusing on the influence of the amount of CFRP layers (nc), the slenderness ratio (λ), the initial geometric imperfections (v0), the CFRP layouts (2S and 4S) and the length of the exposed steel column (Le). The load-displacement curves, the bearing force, and typical failure modes were also acquired. Results indicated that with increasing nc and v0, and decreasing λ, the conventional CFRP-SHSCMs were prone to local bearing failure with poor ductility, leading to the insufficient use of the CFRP plate, in contrast, the improved CFRP-SHSCMs primarily underwent overall buckling failure and exhibited better bearing force and ductility. Finally, the modified Perry-Robertson formula was put forward to predict the ultimate load of the CFRP-SHSCMs. The coefficients of variation between the FE simulation and the theoretical results were 0.00436 and 0.0292, respectively.

Measuring absorbed energy in the human auditory system using finite element models: A comparison with experimental results

BACKGROUND

There are different ways to analyze energy absorbance (EA) in the human auditory system. In previous research, we developed a complete finite element model (FEM) of the human auditory system.

OBJECTIVE

In this current work, the external auditory canal (EAC), middle ear, and inner ear (spiral cochlea, vestibule, and semi-circular canals) were modelled based on human temporal bone histological sections.

METHODS

Multiple acoustic, structure, and fluid-coupled analyses were conducted using the FEM to perform harmonic analyses in the 0.1-10 kHz range. Once the FEM had been validated with published experimental data, its numerical results were used to calculate the EA or energy reflected (ER) by the tympanic membrane. This EA was also measured in clinical audiology tests which were used as a diagnostic parameter.

RESULTS

A mathematical approach was developed to calculate the EA and ER, with numerical and experimental results showing adequate correlation up to 1 kHz. Another published FEM had adapted its boundary conditions to replicate experimental results. Here, we recalculated those numerical results by applying the natural boundary conditions of human hearing and found that the results almost totally agreed with our FEM.

CONCLUSION

This boundary problem is frequent and problematic in experimental hearing test protocols: the more invasive they are, the more the results are affected. One of the main objectives of using FEMs is to explore how the experimental test conditions influence the results. Further work will still be required to uncover the relationship between middle ear structures and EA to clarify how to best use FEMs. Moreover, the FEM boundary conditions must be more representative in future work to ensure their adequate interpretation.

Ex vivo, in vivo and in silico studies of corneal biomechanics: a systematic review

Healthy cornea guarantees the refractive power of the eye and the protection of the inner components, but injury, trauma or pathology may impair the tissue shape and/or structural organization and therefore its material properties, compromising its functionality in the ocular visual process. It turns out that biomechanical research assumes an essential role in analysing the morphology and biomechanical response of the cornea, preventing pathology occurrence, and improving/optimising treatments. In this review, ex vivo, in vivo and in silico methods for the corneal mechanical characterization are reported. Experimental techniques are distinct in testing mode (e.g., tensile, inflation tests), samples' species (human or animal), shape and condition (e.g., healthy, treated), preservation methods, setup and test protocol (e.g., preconditioning, strain rate). The meaningful results reported in the pertinent literature are discussed, analysing differences, key features and weaknesses of the methodologies adopted. In addition, numerical techniques based on the finite element method are reported, incorporating the essential steps for the development of corneal models, such as geometry, material characterization and boundary conditions, and their application in the research field to extend the experimental results by including further relevant aspects and in the clinical field for diagnostic procedure, treatment and planning surgery. This review aims to analyse the state-of-art of the bioengineering techniques developed over the years to study the corneal biomechanics, highlighting their potentiality to improve diagnosis, treatment and healing process of the corneal tissue, and, at the same, pointing out the current limits in the experimental equipment and numerical tools that are not able to fully characterize in vivo corneal tissues non-invasively and discourage the use of finite element models in daily clinical practice for surgical planning.

Regional shear wave speeds track regional axial stress in nonuniformly loaded fibrous soft tissues

Ligaments and tendons undergo nonuniform deformation during movement. While deformations can be imaged, it remains challenging to use such information to infer regional tissue loading. Shear wave tensiometry is a promising noninvasive technique to gauge axial stress and is premised on a tensioned beam model. However, it is unknown whether tensiometry can predict regional stress in a nonuniformly loaded structure. The objectives of this study were to (1) determine whether regional shear wave speed tracks regional axial stress in nonuniformly loaded fibrous soft tissues, and (2) determine the sensitivity of regional axial stress and shear wave speed to nonuniform load distribution and fiber alignment. We created a representative set of 12,000 dynamic finite element models of a fibrous soft tissue with probabilistic variations in fiber alignment, stiffness, and aspect ratio. In each model, we applied a randomly selected nonuniform load distribution, and then excited a shear wave and tracked its regional propagation. We found that regional shear wave speed was an excellent predictor of the regional axial stress (RMSE = 0.57 MPa) and that the nature of the regional shear wave speed-stress relationship was consistent with a tensioned beam model (R2 = 0.99). Variations in nonuniform load distribution and fiber alignment did not substantially alter the wave speed-stress relationship, particularly at higher loads. Thus, these findings suggests that shear wave tensiometry could provide a quantitative estimate of regional tissue stress in ligaments and tendons.

Finite element model of ocular adduction with unconstrained globe translation

Details of the anatomy and behavior of the structures responsible for human eye movements have been extensively elaborated since the first modern biomechanical models were introduced. Based on these findings, a finite element model of human ocular adduction is developed based on connective anatomy and measured optic nerve (ON) properties, as well as active contractility of bilaminar extraocular muscles (EOMs), but incorporating the novel feature that globe translation is not otherwise constrained so that realistic kinematics can be simulated. Anatomy of the hemisymmetric model is defined by magnetic resonance imaging. The globe is modeled as suspended by anatomically realistic connective tissues, orbital fat, and contiguous ON. The model incorporates a material subroutine that implements active EOM contraction based on fiber twitch characteristics. Starting from the initial condition of 26° adduction, the medial rectus (MR) muscle was commanded to contract as the lateral rectus (LR) relaxed. We alternatively modeled absence or presence of orbital fat. During pursuit-like adduction from 26 to 32°, the globe translated 0.52 mm posteriorly and 0.1 mm medially with orbital fat present, but 1.2 mm posteriorly and 0.1 mm medially without fat. Maximum principal strains in the optic disk and peripapillary reached 0.05-0.06, and von-Mises stress 96 kPa. Tension in the MR orbital layer was ~ 24 g-force after 6° adduction, but only ~ 3 gm-f in the whole LR. This physiologically plausible simulation of EOM activation in an anatomically realistic globe suspensory system demonstrates that orbital connective tissues and fat are integral to the biomechanics of adduction, including loading by the ON.

Investigation of stair ascending and descending activities on the lifespan of hip implants

Total hip arthroplasty (THA) surgeries among young patients are on the increase, so it is crucial to predict the lifespan of hip implants correctly and produce solutions to improve longevity. Current implants are designed and tested against walking conditions to predict the wear rates. However, it would be reasonable to include the additional effects of other daily life activities on wear rates to predict convergent results to clinical outputs. In this study, 14 participants are recruited to perform stair ascending (AS), descending (DS), and walking activities to obtain kinematic and kinetic data for each cycle using marker based Qualisys motion capture (MOCAP) system. AnyBody Modeling System using the Calibrated Anatomical System Technique (CAST) full body marker set are performed Multibody simulations. The 3D generic musculoskeletal model used in this study is a marker-based full-body motion capture model (AMMR,2.3.1 MoCapModel) consisting of the upper extremity and the Twente Lower Extremity Model (TLEM2). The dynamic wear prediction model detailing the intermittent and overall wear rates for CoCr-on-XLPE bearing couple is developed to investigate the wear mechanism under 3D loading for AS, DS, and walking activities over 5 million cycles (Mc) by using finite element modelling technique. The volumetric wear rates of XLPE liner under AS, DS, and walking activities over 5-Mc are predicted as 27.43, 23.22, and 18.84 mm3/Mc respectively. Additionally, the wear rate was predicted by combining stair activities and gait cycles based on the walk-to-stair ratio. By adding the effect of stair activities, the volumetric wear rate of XLPE is predicted as 22.02 mm3/Mc which is equivalent to 19.41% of walking. In conclusion, in this study, the effect of including other daily life activities is demonstrated and evidence is provided by matching them to the clinical data as opposed to simulator test results of implants under ISO 14242 boundary conditions.

Investigation of Head Kinematics and Brain Strain Response During Soccer Heading Using a Custom-Fit Instrumented Mouthguard

Association football, also known as soccer in some regions, is unique in encouraging its participants to intentionally use their head to gain a competitive advantage, including scoring a goal. Repetitive head impacts are now being increasingly linked to an inflated risk of developing long-term neurodegenerative disease. This study investigated the effect of heading passes from different distances, using head acceleration data and finite element modelling to estimate brain injury risk. Seven university-level participants wore a custom-fitted instrumented mouthguard to capture linear and angular acceleration-time data. They performed 10 headers within a laboratory environment, from a combination of short, medium, and long passes. Kinematic data was then used to calculate peak linear acceleration, peak angular velocity, and peak angular acceleration as well as two brain injury metrics: head injury criterion and rotational injury criterion. Six degrees of freedom acceleration-time data were also inputted into a widely accepted finite element brain model to estimate strain-response using mean peak strain and cumulative strain damage measure values. Five headers were considered to have a 25% concussion risk. Mean peak linear acceleration equalled 26 ± 7.9 g, mean peak angular velocity 7.20 ± 2.18 rad/s, mean peak angular acceleration 1730 ± 611 rad/s2, and 95th percentile mean peak strain 0.0962 ± 0.252. Some of these data were similar to brain injury metrics reported from American football, which supports the need for further investigation into soccer heading.

Computational Analysis of the Effects of Fiber Deformation on the Microstructure and Permeability of Blood Oxygenator Bundles

Mechanical loads on the polymeric fibers of oxygenating bundles are commonly present due to bundle press-fitting during device assembly and blood pressure load. However, computational fluid dynamics (CFD) simulations for fiber bundle optimization neglect possible changes in microstructure due to such deformations. The aim of this study is to investigate the impact of fiber deformability on bundle microstructure and fluid dynamics mainly in terms of permeability. Fibers from commercial mats typically used for blood oxygenators were mechanically tested and based on these experimental data, a material model was developed to simulate the structural deformations the fibers undergo under press-fitting and blood pressure loads. Then, CFD simulations were performed on deformed bundle repetitive units to investigate permeability under varying loading conditions. The effects of different bundle geometric parameters on the variation of bundle permeability due to press-fitting were evaluated. Bundle press-fitting results in significant changes in microstructure that are reflected in a bundle permeability more than halved for a 15% press-fitting. This impact on permeability is present in all the simulated fiber bundles and becomes more pronounced as the pitch between fibers and thus bundle porosity decreases. Instead, the analyses on pressurized bundle show only small deformations caused by pressure load, with permeability changes below 1%. While blood pressure effects could be neglected, bundle press-fitting turns out to have a significant impact on bundle microstructure and permeability. Neglecting such microstructure variations during CFD simulations could also lead to incorrect assessment of the local fluid dynamics within the bundle.

Development and validation of a semi-automated and unsupervised method for femur segmentation from CT

Quantitative computed tomography (QCT)-based in silico models have demonstrated improved accuracy in predicting hip fractures with respect to the current gold standard, the areal bone mineral density. These models require that the femur bone is segmented as a first step. This task can be challenging, and in fact, it is often almost fully manual, which is time-consuming, operator-dependent, and hard to reproduce. This work proposes a semi-automated procedure for femur bone segmentation from CT images. The proposed procedure is based on the bone and joint enhancement filter and graph-cut algorithms. The semi-automated procedure performances were assessed on 10 subjects through comparison with the standard manual segmentation. Metrics based on the femur geometries and the risk of fracture assessed in silico resulting from the two segmentation procedures were considered. The average Hausdorff distance (0.03 ± 0.01 mm) and the difference union ratio (0.06 ± 0.02) metrics computed between the manual and semi-automated segmentations were significantly higher than those computed within the manual segmentations (0.01 ± 0.01 mm and 0.03 ± 0.02). Besides, a blind qualitative evaluation revealed that the semi-automated procedure was significantly superior (p < 0.001) to the manual one in terms of fidelity to the CT. As for the hip fracture risk assessed in silico starting from both segmentations, no significant difference emerged between the two (R2 = 0.99). The proposed semi-automated segmentation procedure overcomes the manual one, shortening the segmentation time and providing a better segmentation. The method could be employed within CT-based in silico methodologies and to segment large volumes of images to train and test fully automated and supervised segmentation methods.

Optimizing 3D printed ankle-foot orthoses for patients with stroke: Importance of effective elastic modulus and finite element simulation

Patients with stroke often use ankle-foot orthoses (AFOs) for gait improvement. 3D printing technology has become a popular tool in recent years for the production of AFOs due to its strengths on customization and rapid manufacturing. However, the porosity of the 3D printed materials affects the kinetic features of these orthoses, leading to its lower-strength than solid ones. The effective elastic modulus of 3D printed material was measured following standard test method to obtain the kinetic features precisely in a finite element simulation. This study demonstrated that the porosity of 3D printed samples using 100% fill density was 11% for PLA and 16% for Nylon. As a result, their effective elastic modulus was reduced to 1/3 and 1/12 of fully solid objects, respectively, leading to a lower stiffness of 3D printed orthoses. A fatigue testing platform was built to verify our finite element model, and the findings of the fatigue test were consistent with the analysis of the finite element model. Further, our AFO has been proven to have a lifespan exceeding 200 thousand steps. Our study highlights the significance of determining the actual porosity of 3D printed samples by calculating the effective elastic modulus, which leads to a more precise finite element simulation and enables reliable prediction of the kinetic features of the AFO. Overall, this study provides valuable insights into the production and optimization of 3D printed AFOs for patients with stroke.

Analysis of Extracellular ATP Distribution in the Intervertebral Disc

Progressive loss of proteoglycans (PGs) is the major biochemical change during intervertebral disc (IVD) degeneration. Adenosine triphosphate (ATP) as the primary energy source is not only critical for cell survival but also serves as a building block in PG synthesis. Extracellular ATP can mediate a variety of physiological functions and was shown to promote extracellular matrix (ECM) production in the IVD. Therefore, the objective of this study was to develop a 3D finite element model to predict extracellular ATP distribution in the IVD and evaluate the impact of degeneration on extracellular ATP distribution. A novel 3D finite element model of the IVD was developed by incorporating experimental measurements of ATP metabolism and ATP-PG binding kinetics into the mechano-electrochemical mixture theory. The new model was validated by experimental data of porcine IVD, and then used to analyze the extracellular distribution of ATP in human IVDs. Extracellular ATP was shown to bind specifically with PGs in IVD ECM. It was found that annulus fibrosus cells hydrolyze ATP faster than that of nucleus pulposus (NP) cells whereas NP cells exhibited a higher ATP release. The distribution of extracellular ATP in a porcine model was consistent with experimental data in our previous study. The predictions from a human IVD model showed a high accumulation of extracellular ATP in the NP region, whereas the extracellular ATP level was reduced with tissue degeneration. This study provides an understanding of extracellular ATP metabolism and its potential biological influences on the IVD via purinergic signaling.

Smoothed particle hydrodynamics implementation to enhance vertebral fracture finite element model in a cervical spine segment under compression

Spinal cord injuries (SCIs) can arise from compression loading when a vertebra fractures and bone fragments are pushed into the spinal canal. Experimental studies have demonstrated the importance of both fracture initiation and post-fracture response in the investigation of vertebral fractures and spinal canal occlusion resulting from compression. Finite element models, such as the Global Human Body Models Consortium (GHBMC) model, focused on predicting the initiation location of fractures using element erosion to model hard tissue fracture. However, the element erosion method resulted in a loss of material and structural support during compression, which limited the ability of the model to predict the post-fracture response. The current study aimed to improve the post-fracture response by combining strain-based element erosion with smoothed particle hydrodynamics (SPH) to preserve the volume of the trabecular bone during compression fracture. The proposed implementation was evaluated using a model comprising two functional spinal units (FSUs) (C5-C6-C7) extracted from the GHBMC 50th percentile male model, and loaded under central compression. The original and enhanced models were compared to experimental force-displacement data and measured occlusion of the spinal canal. The enhanced model with SPH improved the shape and magnitude of the force-displacement response to be in good agreement with the experimental data. In contrast to the original model, the enhanced SPH model demonstrated occlusion on the same order of magnitude as reported in the experiments. The SPH implementation improved the post-fracture response by representing the damaged material post-fracture, providing structural support throughout compression loading and material flow leading to occlusion.

Head injury: Importance of the deep brain nuclei in force transmission to the brain

Finite element modeling provides a digital representation of the human body. It is currently the most pertinent method to study the mechanisms of head injury, and is becoming a scientific reference in forensic expert reports. Improved biofidelity is a recurrent aim of research studies in biomechanics in order to improve earlier models whose mechanical properties conformed to simplified elastic behavior and mechanic laws. We aimed to study force transmission to the brain following impacts to the head, using a finite element head model with increased biofidelity. To the model developed by the Laboratory of Applied Biomechanics of Marseille, we added new brain structures (thalamus, central gray nuclei and ventricular systems) as well as three tracts involved in the symptoms of head injury: the corpus callosum, uncinate tracts and corticospinal tracts. Three head impact scenarios were simulated: an uppercut with the prior model and an uppercut with the improved model in order to compare the two models, and a lateral impact with an impact velocity of 6.5 m/s in the improved model. In these conditions, in uppercuts the maximum stress values did not exceed the injury risk threshold. On the other hand, the deep gray matter (thalamus and central gray nuclei) was the region at highest risk of injury during lateral impacts. Even if injury to the deep gray matter is not immediately life-threatening, it could explain the chronic disabling symptoms of even low-intensity head injury.

An optimal de Quervain's tenosynovitis splint with ergonomic thumb support and evenly distributed pressure

Splinting is a conventional treatment for de Quervain's tenosynovitis (dQt). However, existing splints have problems such as excessive thermal discomfort and poor fit, which have been pointed out in previous studies. This study proposes a new functional splint consisting of both hard and soft materials with the aim of providing wear comfort with a good fit and sufficient stability of the injured hand. Thumb support of the splint is an important component that controls and protects the affected thumb. To develop an ergonomically shaped thumb support, 16 participants with dQt were recruited for three-dimensional (3D) scanning of their hands. The angles of the wrist and the curvature of the thumb were measured using computer software, and the results were used as a reference for the design of the prototype supports. Excessive pressure on particular regions, such as bony areas, may cause discomfort or pain. To ensure the wear comfort of the proposed splint, a finite element model (FEM) was built to simulate the wear process of the splint and hence to predict the pressure distribution exerted from the splint onto the hand of the wearer. The simulated results show that the pressure is evenly distributed over the hand, indicating that patients are likely to wear the proposed splint comfortably during their treatment period.

Investigation on vehicle occupant dummy applicability for under-foot impact loading conditions

PURPOSE

Under-foot impact loadings can cause serious lower limb injuries in many activities, such as automobile collisions and underbody explosions to military vehicles. The present study aims to compare the biomechanical responses of the mainstream vehicle occupant dummies with the human body lower limb model and analyze their robustness and applicability for assessing lower limb injury risk in under-foot impact loading environments.

METHODS

The Hybrid III model, the test device for human occupant restraint (THOR) model, and a hybrid human body model with the human active lower limb model were adopted for under-foot impact analysis regarding different impact velocities and initial lower limb postures.

RESULTS

The results show that the 2 dummy models have larger peak tibial axial force and higher sensitivity to the impact velocities and initial postures than the human lower limb model. In particular, the Hybrid III dummy model presented extremely larger peak tibial axial forces than the human lower limb model. In the case of minimal difference in tibial axial force, Hybrid III's tibial axial force (7.5 kN) is still 312.5% that of human active lower limb's (2.4 kN). Even with closer peak tibial axial force values, the biomechanical response curve shapes of the THOR model show significant differences from the human lower limb model.

CONCLUSION

Based on the present results, the Hybrid III dummy cannot be used to evaluate the lower limb injury risk in under-foot loading environments. In contrast, potential improvement in ankle biofidelity and related soft tissues of the THOR dummy can be implemented in the future for better applicability.

Finite element biomechanical analysis of 3D printed intervertebral fusion cage in osteoporotic population

OBJECTIVE

To study the biomechanical characteristics of each tissue structure when using different 3D printing Cage in osteoporotic patients undergoing interbody fusion.

METHODS

A finite element model of the lumbar spine was reconstructed and validated with regarding a range of motion and intervertebral disc pressure from previous in vitro studies. Cage and pedicle screws were implanted and part of the lamina, spinous process, and facet joints were removed in the L4/5 segment of the validated mode to simulate interbody fusion. A 280 N follower load and 7.5 N·m moment were applied to different postoperative models and intact osteoporotic model to simulate lumbar motion. The biomechanical characteristics of different models were evaluated by calculating and analyzing the range of motion of the fixed and cephalic adjacent segment, the stress of the screw-rod system, the stress at the interface between cage and L5 endplate, and intervertebral disc pressure of the adjacent segment.

RESULTS

After rigid fixation, the range of motion of the fixed segment of model A-C decreased significantly, which was much smaller than that of the osteoporotic model. And with the increase of the axial area of the interbody fusion cages, the fixed segment of model A-C tended to be more stable. The range of motion and intradiscal pressure of the spinal models with different interbody fusion cages were higher than those of the complete osteoporosis model, but there was no significant difference between the postoperative models. On the other hand, the L5 upper endplate stress and screw-rod system stress of model A-C show a decreasing trend in different directions of motion. The stress of the endplate is the highest during flexion, which can reach 40.5 MPa (model A). The difference in endplate stress between models A-C was the largest during lateral bending. The endplate stress of models A and B was 150.5% and 140.9% of that of model C, respectively. The stress of the screw-rod system was the highest during lateral bending (model A, 102.0 MPa), which was 108.4%, 102.4%, 110.4%, 114.2% of model B and 158.5%, 110.1%, 115.8%, 125.4% of model C in flexion, extension, lateral bending, and rotation, respectively.

CONCLUSIONS

For people with osteoporosis, no matter what type of cage is used, good immediate stability can be achieved after surgery. Larger cage sizes provide better fixation without significantly increasing ROM and IDP in adjacent segments, which may contribute to the development of ASD. In addition, larger cage sizes can disperse endplate stress and reduce stress concentration, which is of positive significance in preventing cage subsidence after operation. The cage and screw rod system establish a stress conduction pathway on the spine, and a larger cage greatly enhances the stress-bearing capacity of the front column, which can better distribute the stress of the posterior spine structure and the stress borne by the posterior screw rod system, reduce the stress concentration phenomenon of the nail rod system, and avoid exceeding the yield strength of the material, resulting in the risk of future instrument failure.

Biomechanical evaluation of the novel assembled internal fixed system in C2-C3 anterior cervical discectomy and fusion: a finite element analysis

BACKGROUND

To analyze and compare the biomechanical characteristics of the new combined cervical fusion device (NCCFD) and the traditional cage-plate construct (CPC) to ascertain its effectiveness in anterior cervical discectomy and fusion (ACDF) using finite element analysis.

METHODS

A finite element model of the cervical spine, inclusive of the occipital bone was created and validated. In the ACDF model, either CPC or NCCFD was implanted at the C2-C3 segment of the model. A pure moment of 1.0 Nm combined with a follower load of 50 N was directed onto the superior surfaces of the occipital bone to determine flexion, extension, lateral bending (left and right), and axial rotation (left and right). The range of motion (ROM), stress distribution at the bone-implant interface, and facet joint forces were investigated and compared between CPC and NCCFD systems.

RESULT

The results showed that the ROMs of the fused levels in both models were nearly zero, and the motions of the unfused segments were similar. In addition, the maximum displacement exhibited nearly identical values for both models. The maximum stress of NCCFD screws in lateral bending and rotational conditions is significantly higher than that of the CPC, while the NCCFD model's maximum stress remains within an acceptable range. Comparing the maximum fusion stress, it was found that the CPC experiences much lower fusion stress in anterior flexion and extension than the NCCFD, with no significant difference between the two in lateral bending and rotational states. Stress on the cage was mainly concentrated on both sides of the wings. Comparing the maximum IDP in the CPC and NCCFD, it was observed that maximum stresses rise in extension and lateral bending for both models. Lastly, stress distributions of the facet joints were generally similar across the two devices.

CONCLUSION

NCCFD not only provides the same level of biomechanical stability as CPC but also avoids postoperative complications associated with uneven force damage to the implant. The device offers a novel surgical alternative for ACDF in C2-C3 level.

Modeling of single-lap joints with auxetic adhesive utilizing homogeneous and heterogeneous bondline microstructures

Producing lightweight structures can be effectively achieved using adhesive bonding that incorporates a high-performance adhesive, eliminating the mechanical joints (rivet, bolts). Here, the proof-of-concept of tailoring the high-performance adhesive using microstructures exhibiting a negative Poisson's ratio (auxetic) is investigated using two single lap joint (SLJ) models developed in ABAQUS/Standard. The SLJ models consist of two rigid adherends that are bonded using either homogeneous or heterogeneous adhesive that produces auxetic response. In homogeneous adhesive, a negative value of Poisson's ratio is defined in the adhesive part of the models. In heterogeneous adhesive, the negative Poisson's ratio is obtained by explicitly building orthogonally-arranged elliptical voids in the adhesive part of the models. In addition, the adherend-adhesive interface is represented by cohesive elements with bi-linear traction-separation rule. The effects of using two different adhesive models and thickness on peel and shear stresses, as well as failure mechanism and joint strength, are evaluated. We found that the effect of auxetic adhesive on SLJ response is more profound in the model using heterogeneous adhesive rather than the homogeneous one, where the joint strength enhancement could achieve 45 % in reference to the baseline model. The heterogeneous adhesive embedded between adherends is able to activate ligaments, the mechanism that could not be obtained using merely homogeneous adhesive. Our proposed modeling strategy can be a starting point to further the modeling approach of bonded joints utilizing auxetic adhesive.

Human body numerical simulation: An accurate model for a thigh subjected to a cold treatment

This article presents the development of a digital twin model of a thigh portion subjected to various thermal treatments. Two scenarios are investigated: cold water immersion (CWI) and whole body cryotherapy (WBC), for which the comparison of numerical results with experimental measurements validates the consistency of the developed model. The use of real geometry on a first subject demonstrates the high heterogeneity of the temperature field and the need for accurate geometry. A second subject with thicker adipose tissue highlights the impact of the subject's actual morphology on the validity of the treatment and the necessity to work with real geometry in order to optimize cold modalities and develop personalized treatments.

Evaluation of bone mineral density in adolescent idiopathic scoliosis using a three-dimensional finite element model: a retrospective study

BACKGROUND

Adolescent idiopathic scoliosis (AIS) is often accompanied by osteopenia and osteoporosis, which can cause serious complications. The aim of this study was to determine the specific bone mineral density (BMD) of each vertebral body in patients with AIS using biomechanical finite element modeling based on three-dimensional (3D) reconstruction.

METHODS

This retrospective study involved 56 patients with AIS. Computed tomography (CT) and radiography were performed. Spinal vertebrae were segmented from the spinal CT images of patients with AIS to reconstruct 3D vertebral models. The vertebral models were meshed into tetrahedral finite elements to assess the BMD.

RESULTS

The mean main curve Cobb angle was 88.6 ± 36.7°, and the mean kyphosis angle was 36.8 ± 31.5°. The mean BMD of the global spine was 0.83 ± 0.15 g/cm2. The highest BMD was measured on the concave side of the apex (0.98 ± 0.16 g/cm2). Apical vertebral BMD was negatively correlated with age and height (r = - 0.490, p = 0.009 and r =  - 0.478, p = 0.043, respectively). There were no significant differences in BMD values between the concave and convex sides (p > 0.05).

CONCLUSIONS

The 3D finite element modeling of BMD in patients with AIS is a reliable and accurate BMD measurement method. Using this method, the overall BMD of patients with AIS was shown to gradually decrease from the top to the bottom of the spine. Our findings provide valuable insights for surgical planning, choice of screw trajectories, and additional biomechanical analyzes using finite element models in the context of scoliosis.

The major influence of anterior and equatorial zonular fibres on the far-to-near accommodation revealed by a 3D pre-stressed model of the anterior eye

PURPOSE

To explore the synergistic function of the ligaments in eye, the zonular fibres, that mediate change in eye lens shape to allow for focussing over different distances.

METHODS

A set of 3D Finite Element models of the anterior eye together with a custom developed pre-stress modelling approach was proposed to simulate vision for distant objects (the unaccommodated state) to vision for near objects (accommodation). One of the five zonular groups was cut off in sequence creating five models with different zonular arrangements, the contribution of each zonular group was analysed by comparing results of each specific zonular-cut model with those from the all-zonules model in terms of lens shape and zonular tensions.

RESULTS

In the all-zonular model, the anterior and equatorial zonules carry the highest tensions. In the anterior zonular-cut model, the equatorial zonular tension increases while the posterior zonular tension decreases, resulting in an increase in the change in Central Optical Power (COP). In the equatorial zonular-cut model, both the anterior and posterior zonular tensions increase, causing a decreasing change in COP. The change in COP decreases only slightly in the other models. For vitreous zonular-cut models, little change was seen in either the zonular tension or the change in COP.

CONCLUSIONS

The anterior and the equatorial zonular fibres have the major influence on the change in lens optical power, with the anterior zonules having a negative effect and the equatorial zonules contributing a positive effect. The contribution to variations in optical power by the equatorial zonules is much larger than by the posterior zonules.

Numerical investigation of intra-abdominal pressure and spinal load-sharing upon the application of an abdominal belt

Chronic low back pain patients may experience spinal instability. Abdominal belts (ABs) have been shown to improve spine stability, trunk stiffness, and resiliency to spinal perturbations. However, research on the contributing mechanisms is inconclusive. ABs may increase intra-abdominal pressure (IAP) and reduce paraspinal soft tissue contribution to spine stability without increasing spinal compressive loads. A finite element model (FEM) of the spine inclusive of the T1-S1 vertebrae, intervertebral discs (IVDs), ribcage, pelvis, soft tissues, and abdominal cavity, without active muscle forces was developed. An identical FEM with an AB was developed. Both FEMs underwent trunk flexion. Following validation, the models' intervertebral rotation (IVR), IAP, IVD pressure, and tensile stress in the multifidus (MF), erector spinae (ES), and thoracolumbar fascia (TLF) were compared. The inclusion of an AB resulted in a 3.8 kPa IAP increase, but a decreased average soft tissue tensile stress of 0.28 kPa. The TLF withstood the majority of tension being transferred across the paraspinal soft tissues (>70 %). The average IVR in the AB model decreased by 10 %, with the lumbar spine experiencing the largest reduction. The lumbar IVDs of the AB model likewise showed a 31 % reduction in average IVD pressure. Using an AB improved trunk bending stiffness, primarily in the lumbar spine. Wearing an AB had minimal effect on reducing tensile stress in theES. The skewed stress distribution towards the TLF suggests its large contribution to spine stability and the potential advantage in unloading the structure when wearing an AB, measured herein at8 %.

Biomechanical effects of typical lower limb movements of Chen-style Tai Chi on knee joint

The load and stress distribution on cartilage and meniscus of the knee joint in typical lower limb movements of Chen-style Tai Chi (TC) and deep squat (DS) were analyzed using finite element (FE) analysis. The loadings for this analysis consisted of muscle forces and ground reaction force (GRF), which were calculated through the inverse dynamic approach based on kinematics and force plate measurements obtained from motion capture experiments. Thirteen experienced practitioners performed four typical TC movements, namely, single whip (SW), brush knee and twist step (BKTS), stretch down (SD), and part the wild horse's mane (PWHM), which exhibit lower posture and greater lower limb force compared to other TC styles. The results indicated that TC required greater lower limb muscle strength than DS, resulting in greater knee joint forces. The stress on the medial cartilage in SW and BKTS fell within a range conductive to maintaining the balance between anabolism and catabolism of cartilage matrix. This was due to the fact that SW and BKTS reduce the medial to total tibiofemoral contact force ratios through knee abduction, which may effectively alleviate mild medial knee osteoarthritis (KOA). However, the greater medial contact force ratios observed in SD and PWHM resulted in great contact stresses that may aggravate the pain of patients with KOA. To mitigate these effects, practitioners should consider elevating their postures appropriately to reduce knee flexion angles, especially during the single-leg support phase. This adjustment can decrease the required muscle strength, load and stress on the knee joint.

QCT-based computational bone strength assessment updated with MRI-derived 'hidden' microporosity

Microdamage accumulated through sustained periods of cyclic loading or single overloading events contributes to bone fragility through a reduction in stiffness and strength. Monitoring microdamage in vivo remains unattainable by clinical imaging modalities. As such, there are no established computational methods for clinical fracture risk assessment that account for microdamage that exists in vivo at any specific timepoint. We propose a method that combines multiple clinical imaging modalities to identify an indicative surrogate, which we term 'hidden porosity', that incorporates pre-existing bone microdamage in vivo. To do so, we use the third metacarpal bone of the equine athlete as an exemplary model for fatigue induced microdamage, which coalesces in the subchondral bone. N = 10 metacarpals were scanned by clinical quantitative computed tomography (QCT) and magnetic resonance imaging (MRI). We used a patch-based similarity method to quantify the signal intensity of a fluid sensitive MRI sequence in bone regions where microdamage coalesces. The method generated MRI-derived pseudoCT images which were then used to determine a pre-existing damage (Dpex) variable to quantify the proposed surrogate and which we incorporate into a nonlinear constitutive model for bone tissue. The minimum, median, and maximum detected Dpex of 0.059, 0.209, and 0.353 reduced material stiffness by 5.9%, 20.9%, and 35.3% as well as yield stress by 5.9%, 20.3%, and 35.3%. Limb-specific voxel-based finite element meshes were equipped with the updated material model. Lateral and medial condyles of each metacarpal were loaded to simulate physiological joint loading during gallop. The degree of detected Dpex correlated with a relative reduction in both condylar stiffness (p = 0.001, R2 > 0.74) and strength (p < 0.001, R2 > 0.80). Our results illustrate the complementary value of looking beyond clinical CT, which neglects the inclusion of microdamage due to partial volume effects. As we use clinically available imaging techniques, our results may aid research beyond the equine model on fracture risk assessment in human diseases such as osteoarthritis, bone cancer, or osteoporosis.

Changes in load distribution after unilateral condylar fracture: A finite element model study

OBJECTIVE

Premature dental contact on the fractured side and a contralateral open bite are signs of a unilaterally fractured condyle of the temporomandibular joint (TMJ). The lateral pterygoid muscle pulls the condyle inwards, causing angulation of the fractured part and shortening of the ramus. This imbalance after fracture might change the load in both TMJs and consequently induce remodeling. The present study aimed to calculate this change in load. It is hypothesized to decrease on the fractured side and increase on the non-fractured side.

DESIGN

For these calculations, a finite element model (FEM) was used. In the FEM, shortening of the ramus varied from 2 mm to 16 mm; angulation, from 6.25° to 50°.

RESULTS

After fracture, load on the non-fractured side increased, but only at maximal mouth opening (MMO). Simultaneously, load on the fractured side decreased, at both timepoints, i.e., MMO and closed mouth. When comparing all simulations at those time points, i.e., from 2 mm and 6.25° to 16 mm and 50°, the load in the fractured condyle declines steadily. However, for both timepoints, a threshold stands out around 6 mm shortening and 18.75° angulation: visualization of the fractured condyle showed, apart from load on the condylar head, a second point of load more medial in the TMJ which was most evident in the 6 mm - 18.75° simulation.

CONCLUSIONS

These findings could implicate that the balance between both TMJs is more difficult to restore after a fracture with more than 6 mm shortening and more than 18.75° angulation.

Adjacent tissues modulate shear wave propagation in axially loaded tendons

Shear wave tensiometry is a noninvasive approach for gauging tendon loads based on shear wave speed. Transient shear waves are induced and tracked via sensors secured to the skin overlying a superficial tendon. Wave speeds measured in vivo via tensiometry modulate with tendon load but are lower than that predicted by a tensioned beam model of an isolated tendon, which may be due to the added inertia of adjacent tissues. The objective of this study was to investigate the effects of adjacent fat tissue on shear wave propagation measurements in axially loaded tendons. We created a layered, dynamic finite element model of an elliptical tendon surrounded by subcutaneous fat. Transient shear waves were generated via an impulsive excitation delivered across the tendon or through the subcutaneous fat. The layered models demonstrated dispersive behavior with phase velocity increasing with frequency. Group shear wave speed could be ascertained via dispersion analysis or time-to-peak measures at sequential spatial locations. Simulated wave speeds in the tendon and adjacent fat were similar and modulated with tendon loading. However, wave speed magnitudes were consistently lower in the layered models than in an isolated tendon. For all models, the wave speed-stress relationship was well described by a tensioned beam model after accounting for the added inertia of the adjacent tissues. These results support the premise that externally excited shear waves are measurable in subcutaneous fat and modulate with axial loading in the underlying tendon. The model suggests that adjacent tissues add inertia to the system, which in turn lowers shear wave speeds. This information must be considered when using tensiometry as a clinical or research tool to infer absolute tendon loading.

Contact stress distribution after unilateral condylar fracture with angulation of the fractured part: A finite element model study

After a fracture of the condyle, the head of the condyle is often pulled inwards, which causes the fractured part to angulate medially. This change can cause a disbalance in the masticatory system. The disbalance could lead to contact stress differences within the temporomandibular joints (TMJs) which might induce remodelling within the TMJ to restore the balance. The contact stress in the fractured condyle during open and closing movements is expected to decrease, while the contact stress in the non-fractured condyle will increase. In a clinical situation this is hard to investigate. Therefore, a finite element model (FEM) was used. In the FEM a fractured right condyle with an angulation was induced, which was placed at different degrees, varying from 5° to 50° in steps of 5°. This study shows only minor differences in amount of contact stress between the fractured and the non-fractured condyle. The amount of contact stress in the condyles does not increase with a higher degree of angulation. However, with larger angulations, the contact stress within the fractured condyle is more centralized. Clinically, this more centralized area could be associated with complaints, such as pain. In conclusion, due to the more centralized contact stress in the fractured condyle, one would expect some minor remodelling on the fractured side with more angulation.

Comparison of different surgical treatments for stage II progressive collapsing foot deformity: a finite element analysis

BACKGROUND

This study analyzed the advantages and disadvantages of different procedures for stage IIA progressive collapsing foot deformity (PCFD) through three-dimensional finite element models.

METHODS

A previous validated stage IIA PCFD FEA model was established consisting of 16 bones, 56 ligaments, 5 muscles and soft tissues. The ligament properties of the spring, deltoid, short plantar and long plantar ligaments, and plantar fascia were attenuated according to a previous publication. Medial column fusion (MCF), medializing calcaneal osteotomy (MCO), lateral column lengthening (LCL), and subtalar joint arthroereisis (SJA) operations were simulated in this model. The indexes of plantar stress distribution, maximum von Mises of the medial and lateral columns, strain of the medial ligaments and plantar fascia that supported the medial longitudinal arch, arch height, talo-first metatarsal angle, calcaneus pitch angle, and talonavicular coverage angle were all compared before and after simulated single-foot weight loading.

RESULTS

The maximum plantar stress of PCFD decreased with MCO and SJA but increased with MCF and LCL. MCF and LCL failed to significantly reduce the stress on the medial column fragments, thereby increasing their stress. Both MCO and SJA relieved medial plantar stress. MCF had no significant effect on stress relief of the medial ligament. MCO, LCL, and SJA were all shown to reduce the pressure on the medial plantar ligament, with LCL having the most obvious effect. All four procedures corrected the arch deformity; however, MCF was not as effective as the other methods. SJA is the best method for restoring arch height and correcting arch deformities. For stage IIA PCFD, isolated MCF failed to reduce pressure on the medial column; however, isolated MCO significantly reduced the pressure on the medial plantar and ligamentous soft tissues while restoring the foot's arch and correcting the hindfoot valgus.

CONCLUSION

SJA with type II sinus tarsi implant effectively transferred pressure from the medial plantar tract to the lateral side and restored the arch. Isolated LCL was not found suitable for stage IIA PCFD.

Effect of muscle activation on dynamic responses of neck of pilot during emergency ejection: a finite element study

To determine the effect of muscle activation on the dynamic responses of the neck of a pilot during simulated emergency ejections. A complete finite element model of the pilot's head and neck was developed and dynamically validated. Three muscle activation curves were designed to simulate different activation times and levels of muscles during pilot ejection: A is the unconscious activation curve of the neck muscles, B is the pre-activation curve, and C is the continuous activation curve. The acceleration-time curves obtained during ejection were applied to the model, and the influence of the muscles on the dynamic responses of the neck was investigated by analyzing both angles of rotation of the neck segments and disc stresses. Muscle pre-activation reduced fluctuations in the angle of rotation in each phase of the neck. Continuous muscle activation caused a 20% increase in the angle of rotation compared to pre-activation. Moreover, it resulted in a 35% increase in the load on the intervertebral disc. The maximum stress on the disc occurred in the C4-C5 phase. Continuous muscle activation increased both the axial load on the neck and the posterior extension angle of rotation of the neck. Muscle pre-activation during emergency ejection has a protective effect on the neck. However, continuous muscle activation increases the axial load and rotation angle of the neck. A complete finite element model of the pilot's head and neck was established and three neck muscle activation curves were designed to investigate the effects of muscle activation time and level on the dynamic response of the pilot's neck during ejection. This increased insights into the protection mechanism of neck muscles on the axial impact injury of the pilot's head and neck.

Effectiveness of Wearable Protection Equipment for Seated Pedestrians

This study used finite element models to investigate the efficacy of seated pedestrian protection equipment in vehicle impacts. The selected safety equipment, a lap belt, an airbag vest, and a bicycle helmet, were chosen to mitigate the underlying biomechanical causes of seated pedestrian injuries reported in the literature. The impact conditions were based on the three most dangerous impact scenarios from a previous seated pedestrian impact study. Serious injury (AIS 3+) risks were compared with and without protective equipment. A 50th percentile male occupant model and two generic vehicle models, the family car (FCR) and sports utility vehicle (SUV), were used to simulate vehicle collisions. Three impact conditions were run with every combination of protective equipment (n = 24). The helmet reduced head and brain injury risks from the vehicle-head and ground-head contacts. The airbag reduced the head injury risk in the FCR vehicle-head contact but increased the brain injury risks in the SUV impacts from increased whiplash. The lap belt increased head injury risks for both the FCR and the SUV impacts because it created a stronger FCR vehicle-head contact and SUV ground-head contact. When the belt and airbag were used together the head injury risks dramatically decreased because the pedestrian body impacted the ground arm or leg first and slowly rolled onto the ground which resulted in softer ground-head contacts and in two instances, no ground-head contact. Only the helmet proved effective in all impact conditions. Future testing must be completed before recommending the belt or airbag for seated pedestrians.

Patterns of stress distribution of endodontically treated molar under different types of loading using finite element models-the exploring of mechanism of vertical root fracture

INTRODUCTION

The purpose of this study was to explore the mechanism of vertical root fracture (VRF) using three-dimensional finite element models (FEMs).

METHODS

An endodontically treated mandibular first molar with a subtle VRF was collected and scanned with cone beam CT (CBCT). Three finite element analysis models were created: Model 1 had the actual endodontically treated root canal size; Model 2 had the same root canal size as the contralateral homonymous tooth; and Model 3 had the root canal size expanded by 1 mm based on Model 1. Different types of loading were performed on these 3 FEMs. The stress distribution on the cervical, middle, and apical planes was analyzed, and the maximum stress on the root canal wall was calculated and compared.

RESULTS

In Model 1, the maximum stress around the root canal wall occurred in the cervical part of the mesial root under vertical masticatory force and in the middle part of the mesial root under buccal and lingual lateral masticatory forces. Additionally, there was a stress change zone in a bucco-lingual direction that corresponded with the actual fracture line. In Model 2, the maximum stress around the root canal was in the cervical part of the mesial root under both vertical and buccal lateral masticatory forces. For Model 3, the stress distribution was similar to that of Model 1, but greater under buccal lateral masticatory force and occlusal trauma force. In all three models, the maximum stress around the root canal wall was in the middle part of the distal root under occlusal trauma force.

CONCLUSIONS

The uneven stress around the root canal in the middle part (presented as a stress change zone in a bucco-lingual direction) may be the cause of VRFs.

Can pelvic tilt cause cam morphology? A computational model of proximal femur development mechanobiology

Cam deformity of the proximal femur is a risk factor for early osteoarthritis. While cam morphology is related to mechanical force at a formative time in skeletal growth, the specific problematic forces contributing to the development of cam morphology remain unknown. Individuals with femoroacetabular impingement syndrome exhibit an increased anterior pelvic tilt during walking, which alters their hip joint forces. This study aims to investigate the influence of altered joint force caused by anterior pelvic tilt on proximal femur epiphyseal growth and the potential association between increased anterior pelvic tilt and the development of cam morphology. A computational model is utilized to simulate the endochondral ossification in the proximal femur and predict cam formation. Cartilage growth and ossification patterns for a gait cycle with and without anterior pelvic tilt were modeled. The simulated growth results indicated an increased alpha angle (53° for typically developing to 68° for anterior pelvic tilt) and aspherical femoral head in the model with anterior pelvic tilt. We conclude that anterior pelvic tilt may be sufficient to cause the formation of the cam morphology. Identifying the critical mechanical conditions that increase the risk of cam deformity could help prevent this condition by adjusting the physical activities before skeletal maturity.

Load distribution after unilateral condylar fracture with shortening of the ramus: a finite element model study

OBJECTIVES

After a fracture of the condyle, the fractured ramus is often shortened, which causes premature dental contact on the fractured side and a contralateral open bite. The imbalance could change the load in the temporomandibular joints (TMJs). This change could lead to remodelling of the TMJs to compensate for the imbalance in the masticatory system. The load in the non-fractured condyle is expected to increase, and the load in the fractured condyle to decrease.

MATERIALS AND METHODS

These changes cannot be measured in a clinical situation. Therefore a finite element model (FEM) of the masticatory system was used. In the FEM a fractured right condyle with shortening of the ramus was induced, which varied from 2 to 16 mm.

RESULTS

Results show that, with a larger shortening of the ramus, the load in the fractured condyle decreases and the load in the non-fractured condyle increases. In the fractured condyle during closed mouth a major descent in load, hence a cut-off point, was visible between a shortening of 6 mm and 8 mm.

CONCLUSIONS

In conclusion, the change of load could be associated with remodelling on both condyles due to shortening of the ramus.

CLINICAL RELEVANCE

The cut-off point implies that shortening over 6 mm could present more difficulty for the body to compensate.

Structure and mechanical properties of ladybird elytra as biological sandwich panels

The armour of the ladybird, elytra, protect the body from injury and are well-adapted to flight. However, experimental methods to decipher their mechanical performances had been challenging due to the small size, making it unclear how the elytra balance mass and strength. Here, we provide insights to the relationship between the microstructure and multifunctional properties of the elytra by means of structural characterization, mechanical analysis and finite element simulations. Micromorphology analysis on the elytron revealed the thickness ratio of the upper lamination, middle layer and lower lamination is approximately 51:139:7. The upper lamination had multiple cross fibre layers and the thickness of each fibre layer is not the same. In addition, the tensile strength, elastic modulus, fracture strain, bending stiffness and hardness of elytra were obtained through in-situ tensile and nanoindentation-bending under the influence of multiple loading conditions, which also serve as references for finite element models. The finite element model revealed that structural factors such as thickness of each layer, angle of fibre layer and trabeculae are key to affecting the mechanical properties, but the effect is different. When the thickness of upper, middle and lower layers is the same, the tensile strength provided by unit mass of the model is 52.78% lower than that provided by elytra. These findings broaden the relationship between the structural and mechanical properties of the ladybird elytra, and are expected to inspire the development of sandwich structures in biomedical engineering.

Physiological and degenerative loading of bovine intervertebral disc in a bioreactor: A finite element study of complex motions

Intervertebral disc (IVD) degeneration and regenerative therapies are commonly studied in organ-culture experiments with uniaxial compressive loading. Recently, in our laboratory, we established a bioreactor system capable of applying loads in six degrees-of-freedom (DOF) to bovine IVDs, which replicates more closely the complex multi-axial loading of the IVD in vivo. However, the magnitudes of loading that are physiological (able to maintain cell viability) or mechanically degenerative are unknown for load cases combining several DOFs. This study aimed to establish physiological and degenerative levels of maximum principal strains and stresses in the bovine IVD tissue and to investigate how they are achieved under complex load cases related to common daily activities. The physiological and degenerative levels of maximum principal strains and stresses were determined via finite element (FE) analysis of bovine IVD subjected to experimentally established physiological and degenerative compressive loading protocols. Then, complex load cases, such as a combination of compression + flexion + torsion, were applied on the FE-model with increasing magnitudes of loading to discover when physiological and degenerative tissue strains and stresses were reached. When applying 0.1 MPa of compression and ±2-3° of flexion and ±1-2° of torsion the investigated mechanical parameters remained at physiological levels, but with ±6-8° of flexion in combination with ±2-4° of torsion, the stresses in the outer annulus fibrosus (OAF) exceeded degenerative levels. In the case of compression + flexion + torsion, the mechanical degeneration likely initiates at the OAF when loading magnitudes are high enough. The physiological and degenerative magnitudes can be used as guidelines for bioreactor experiments with bovine IVDs.

Computational Seated Pedestrian Impact Design of Experiments with Ultralight Wheelchair

Pedestrians who use wheelchairs (seated pedestrians) report higher mortality rates than standing pedestrians in vehicle-to-pedestrian collisions but the cause of this mortality is poorly understood. This study investigated the cause of seated pedestrian serious injuries (AIS 3+) and the effect of various pre-collision variables using finite element (FE) simulations. An ultralight manual wheelchair model was developed and tested to meet ISO standards. The GHBMC 50th percentile male simplified occupant model and EuroNCAP family car (FCR) and sports utility vehicle (SUV) were used to simulate vehicle collisions. A full factorial design of experiments (n = 54) was run to explore the effect of pedestrian position relative to the vehicle bumper, pedestrian arm posture, and pedestrian orientation angle relative to the vehicle. The largest average injury risks were at the head (FCR: 0.48 SUV: 0.79) and brain (FCR: 0.42 SUV: 0.50). The abdomen (FCR: 0.20 SUV: 0.21), neck (FCR: 0.08 SUV: 0.14), and pelvis (FCR: 0.02 SUV: 0.02) reported smaller risks. 50/54 impacts reported no thorax injury risk, but 3 SUV impacts reported risks ≥ 0.99. Arm (gait) posture and pedestrian orientation angle had larger effects on most injury risks. The most dangerous arm posture examined was when the hand was off the wheelchair handrail after wheel propulsion and the two more dangerous orientations were when the pedestrian faced 90° and 110° away from the vehicle. Pedestrian position relative to the vehicle bumper played little role in injury outcomes. The findings of this study may inform future seated pedestrian safety testing procedures to narrow down the most concerning impact scenarios and design impact tests around them.

Experimental and numerical investigation of 3D-Printed bone plates under four-point bending load utilizing machine learning techniques

The fused deposition modeling (FDM) technique is widely used to produce components for various applications and has the potential to revolutionize orthopedic research through the production of custom-fit and readily available biomedical implants. The properties of FDM-produced implants are significantly influenced by processing parameters, with layer thickness being a crucial parameter. This study investigated the effect of layer thickness on the flexural properties of Polylactic Acid (PLA) bone plate implants produced by the FDM technique. Experimental results showed that the flexural strength is inversely proportional to the layer thickness due to the variation of voids in the specimens. A 3D finite element (FE) model was developed using Abaqus/Explicit software by incorporating the Gurson-Tvergaard (GT) porous plasticity model to predict the elastoplastic and damage behavior of specimens with different layer thicknesses. The characterization of the elastoplastic and GT parameters was done using a tensile test and by the calibration of a machine learning algorithm. It was shown that the FE model was able to predict the flexural behavior of 3D-printed solid plates with a maximum error of 6.13% in the maximum load. The optimal layer height was found to be 0.1 mm, providing both high flexural strength and adequate bending stiffness.

Dynamic computational wear model of PEEK-on-XLPE bearing couple in total hip replacements

Understanding wear mechanisms is a key factor to prevent primary failures causing revision surgery in total hip replacement (THR) applications. This study introduces a wear prediction model of (Polyetheretherketone) PEEK-on-XLPE (cross-linked polyethylene) bearing couple utilized to investigate the wear mechanism under 3D-gait cycle loading over 5 million cycles (Mc). A 32-mm PEEK femoral head and 4-mm thick XLPE bearing liner with a 3-mm PEEK shell are modeled in a 3D explicit finite element modeling (FEM) program. The volumetric and linear wear rates of XLPE liner per every million cycles were predicted as 1.965 mm³/Mc, and 0.0032 mm/Mc respectively. These results are consistent with the literature. PEEK-on-XLPE bearing couple exhibits a promising wear performance used in THR application. The wear pattern evolution of the model is similar to that of conventional polyethylene liners. Therefore, PEEK could be proposed as an alternative material to the CoCr head, especially used in XLPE-bearing couples. The wear prediction model could be utilized to improve the design parameters with the aim of prolonging the life span of hip implants.

Optimal hinge level in opening wedge high tibial osteotomy: Biomechanical analysis using finite element method

BACKGROUND

While the concept of a safe zone, which can minimize the hinge fracture when performing opening wedge high tibial osteotomy, has been introduced, there is a lack of understanding of the biomechanical environment at the lateral tibial cortex. This study aimed to evaluate the effect of the hinge level on the biomechanical environment at the lateral cortex of the tibia with heterogeneous finite element models.

METHODS

Finite element models of biplanar opening wedge high tibial osteotomy based on computed tomography images of a control subject and three patients with medial compartment knee osteoarthritis were created. In each model, three different hinge levels (proximal, middle, and distal) were set. The process of opening the gap during the operation was simulated, and the maximum von Mises stress values at the lateral tibial cortex were calculated for each hinge level and correction angle.

FINDINGS

The maximum von Mises stress value at the lateral tibial cortex was the lowest when the hinge was at the middle, while the value was the highest when the hinge was at the distal level. Furthermore, it was demonstrated that a higher correction angle yielded a higher probability of lateral tibial cortex fracture.

INTERPRETATION

The findings of this study demonstrate that the hinge at the point where the upper end of the articular cartilage of the proximal tibiofibular joint is located provides the least possibility of lateral tibial cortex fracture, as this is an anatomically independent position from the fibula.

Extended mechanical loads for the analysis of acetabular cages

To analyse the strength and mechanical behaviour of hip implants, it is essential to employ an appropriate loading model. Generating computational models supplemented with muscle forces is a complicated task, especially in the initial phase of implant development. This research aims to expand the possibilities of the simpler acetabular cage model based on joint loads without significantly increasing the demand for computing resources. A Python script covered and grouped the loads from daily activities. The ten calculated major loads were compared with the maximum of the walking and stair climbing loads through the finite element analyses of a custom-made acetabular cage. Sensitivity analyses were performed for the surrounding bones' elastic modulus and the pelvis boundary conditions. The major loads can geometrically cover the entire load spectrum of daily activities. The effect of many high-magnitude force vectors is uncertain in the approach that uses the most common maximum loads. Using these resultant major loads, a new stress concentration area could be detected on the acetabular cage, besides the stress concentration areas induced by the loads reported in the literature. The qualitative correctness of the results is also supported by a control computed tomography scan: a fracture occurred in an extensive, high-stress zone. The results are not sensitive to changes in the elastic modulus of the surrounding bone and the boundary conditions of the model. The presented load vectors and the algorithm make more extensive static analyses possible with little computational overhead. The proposed method can be used for checking the static strength of similar implants.

Equivalent loads from the life-cycle of acetabular cages in relation to bone-graft transformation

BACKGROUND AND OBJECTIVES

Bone grafts placed behind acetabular cages change their structure in response to mechanical stimuli. The full consideration of lifestyle loads is extremely resource-intensive, so a method using substitutive loads was proposed to reduce the calculation cost. The aim of the study is to present and prove this method.

METHODS

By means of mechanical equations and using the force vectors from the literature which have the same initial point and their relative frequency, while applying a linear model, the average strain energy density distribution for all load cases can be calculated, compiling a matrix from the external loads. From the elements of this matrix, three substitutive load vectors can be calculated, which can be proven to produce the same strain energy density distribution by averaging their effects. The feasibility of using this to model the transformation of bone grafts placed behind acetabular cages is demonstrated with a finite element model, along with a reference calculation.

RESULTS

The substitutive load vectors could be calculated in closed form and the simulations showed that they produced a similar density distribution to the reference model with a numerical calculation error range. Accordingly, the density distribution calculated from bone graft transformation is almost the same.

CONCLUSIONS

In addition to the aforementioned linearity and the same initial point limitations, the applied method is able to produce the substitutive load vectors with which the calculation of the strain energy density distribution and the bone graft's new density distributions can be carried out faster.

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.