Characterization of mechanical stress in the occurrence of cortical opacification in age-related cataracts using three-dimensional finite element model of the human lens and RNA-seq

Cataract is the leading cause of blindness across the world. Age-related cataract (ARC) is the most common type of cataract, but its pathogenesis is not fully understood. Using three-dimensional finite element modeling combining experimental biotechnology, our study demonstrates that external forces during accommodation cause mechanical stress predominantly in lens cortex, basically matching the localization of opacities in cortical ARCs. We identified the cellular senescence and upregulation of PIEZO1 mRNA in HLECs under mechanical stretch. This mechano-induced senescence in HLECs might be mediated by PIEZO1-related pathways, portraying a potential biomechanical cause of cortical ARCs. Our study updates the fundamental insight towards cataractogenesis, paving the way for further exploration of ARCs pathogenesis and nonsurgical treatment.

An Improved Finite Element Model of Temporomandibular Joint in Maxillofacial System: Experimental Validation

Finite element (FE) analysis has become a popular method of exploring the biomechanical characteristics of temporomandibular joint (TMJ). However, the FE model should be improved and its reliability should be verified further. This study developed a complete maxillofacial model by cone-beam computed tomography (CBCT) and magnetic resonance imaging (MRI). The integrity and physiological environment of TMJ were considered. Then the FE model and corresponding 3D printed model were developed and loaded under the same conditions. The strains on the mandible and upper surface of the left articular disc were measured on the experimental model and compared with the FE model. The differences of the strains on the mandible were less than 6%. The strain distributions on the disc were also approximate between the experimental and simulated results. It indicated that the strains calculated from the improved FE model were reliable on the mandible and inside the TMJ.

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.

Traumatic axonal injury: Clinic, forensic and biomechanics perspectives

Identification of Traumatic axonal injury (TAI) is critical in clinical practice, particularly in terms of long-term prognosis, but also for medico-legal issues, to verify whether the death or the after-effects were attributable to trauma. Multidisciplinary approaches are an undeniable asset when it comes to solving these problems. The aim of this work is therefore to list the different techniques needed to identify axonal lesions and to understand the lesion mechanisms involved in their formation. Imaging can be used to assess the consequences of trauma, to identify indirect signs of TAI, to explain the patient's initial symptoms and even to assess the patient's prognosis. Three-dimensional reconstructions of the skull can highlight fractures suggestive of trauma. Microscopic and immunohistochemical techniques are currently considered as the most reliable tools for the early identification of TAI following trauma. Finite element models use mechanical equations to predict biomechanical parameters, such as tissue stresses and strains in the brain, when subjected to external forces, such as violent impacts to the head. These parameters, which are difficult to measure experimentally, are then used to predict the risk of injury. The integration of imaging data with finite element models allows researchers to create realistic and personalized computational models by incorporating actual geometry and properties obtained from imaging techniques. The personalization of these models makes their forensic approach particularly interesting.

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.

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.

The C2 isthmus screw provided sufficient biomechanical stability in the setting of atlantoaxial dislocation-a finite element study

BACKGROUND

The emerging of the C2 isthmus screw fixation technique is gaining popularity in the setting of atlantoaxial dislocation or other conditions requiring fixation of C2. However, the biomechanical stability of this fixation is poorly understood.

PURPOSE

To compare and elucidate the biomechanical stability of C2 pedicle screw (C2PS), C2 isthmus screw (C2IS) and C2 short isthmus screw (C2SIS) fixation techniques in atlantoaxial dislocation (AAD).

METHOD

A three-dimensional finite element model (FEM) from occiput to C3 was established and validated from a healthy male volunteer. Three FEMs, C1 pedicle screw (PS)-C2PS, C1PS-C2IS, C1PS-C2SIS were also constructed. The range of motion (ROM) and the maximum von Mises stress under flexion, extension, lateral bending and axial rotation loading were analyzed and compared. The pullout strength of the three fixations for C2 was also evaluated.

RESULT

C1PS-C2IS model showed the greatest decrease in ROM with flexion, extension, lateral bending and axial rotation. C1PS-C2PS model showed the least ROM reduction under all loading conditions than both C2IS and C2SIS. The C1PS-C2PS model had the largest von Mises stress on the screw under all directions followed by C1PS-C2SIS, and lastly the C1PS-C2IS. Under axial rotation and lateral bending loading, the three models showed the maximum and minimum von Mises stress on the screw respectively. The stress of the three models was mainly located in the connection of the screw and rod. Overall, the maximum screw pullout strength for C2PS, C2IS and C2SIS were 729.41N, 816.62N, 640.54N respectively.

CONCLUSION

In patients with atlantoaxial dislocations, the C2IS fixation provided comparable stability, with no significant stress concentration. Furthermore, the C2IS had sufficient pullout strength when compared with C2PS and C2SIS. C2 isthmus screw fixation may be a biomechanically favourable option in cases with AAD. However, future clinical trials are necessary for the evaluation of the clinical outcomes of this technique.

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.

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.

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.

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 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.