Finite element analysis of a three-dimensional cervical spine model with muscles based on CT scan data.

BACKGROUND: The incidence of cervical spondylosis is increasing, gradually affecting people's normal lives. Establishing a finite element model of the cervical spine is one of the methods for studying cervical spondylosis. MRI (Magnetic Resonance Imaging) still has certain difficulties in transitioning from human imaging to establishing muscle models suitable for finite element analysis. Medical software provides specific morphologies and can generate muscle finite element models. Additionally, there is little research on the static analysis of cervical spine finite element models with solid muscle. PURPOSE: A new method is proposed for establishing a finite element model of the cervical spine based on CT (Computed Tomography) data and medical software, and the model's effectiveness is validated. Human movement characteristics based on the force distribution in various parts are analyzed and predicted. METHODS: The muscle model is reconstructed in medical software and a three-dimensional finite element model of the entire cervical spine (C0-C7) is established by combining muscle models with CT vertebral data models. 1.5 Nm of load is applied to the finite element model to simulate the cervical spine movement. RESULTS: The finite element model was successfully established, and effectiveness was verified. Stress variations in various parts under six movements were obtained. The effectiveness of the model was basically verified. CONCLUSION: The finite element model of the cervical spine for mechanical analysis can be successfully established by using medical software and CT data. In daily life, the C2-3, C3-4, C4-C5 intervertebral discs, rectus capitis posterior major, longus colli, and obliquus capitis inferior are more prone to injury.

Preliminary results of classifying otosclerosis and disarticulation using a convolutional neural network trained with simulated wideband acoustic immittance data.

Current noninvasive methods of clinical practice often do not identify the causes of conductive hearing loss due to pathologic changes in the middle ear with sufficient certainty. Wideband acoustic immittance (WAI) measurement is noninvasive, inexpensive and objective. It is very sensitive to pathologic changes in the middle ear and therefore promising for diagnosis. However, evaluation of the data is difficult because of large interindividual variations. Machine learning methods like Convolutional neural networks (CNN) which might be able to deal with this overlaying pattern require a large amount of labeled measurement data for training and validation. This is difficult to provide given the low prevalence of many middle-ear pathologies. Therefore, this study proposes an approach in which the WAI training data of the CNN are simulated with a finite-element ear model and the Monte-Carlo method. With this approach, virtual populations of normal, otosclerotic, and disarticulated ears were generated, consistent with the averaged data of measured populations and well representing the qualitative characteristics of individuals. The CNN trained with the virtual data achieved for otosclerosis an AUC of 91.1 %, a sensitivity of 85.7 %, and a specificity of 85.2 %. For disarticulation, an AUC of 99.5 %, sensitivity of 100 %, and specificity of 93.1 % was achieved. Furthermore, it was estimated that specificity could potentially be increased to about 99 % in both pathological cases if stapes reflex threshold measurements were used to confirm the diagnosis. Thus, the procedures' performance is comparable to classifiers from other studies trained with real measurement data, and therefore the procedure offers great potential for the diagnosis of rare pathologies or early-stages pathologies. The clinical potential of these preliminary results remains to be evaluated on more measurement data and additional pathologies.

[Oblique Lumbar Interbody Fusion Combined With 4-Screw Fixation for Treating Two-Level Degenerative Lumbar Diseases:A Finite Element Study].

Objective To demonstrate the feasibility of oblique lumbar interbody fusion (OLIF) combined with 4-screw fixation for treating two-level lumbar degenerative diseases.Methods An intact finite element model of L3-S1 (M0) was constructed and validated.Then,we constructed the M1 model by simulating OLIF surgery at L3/4 and L4/5 segments on the M0 model.By attachment of posterior 4-screw or 6-screw fixation to the M1 model,three 4-screw fixation models (M2-M4) and one 6-screw fixation model (M5) were established.The segmental and overall range of motion (ROM) and the peak von Mises stresses of superior endplate,cage,and posterior screw-rod were investigated under each implanted condition.Results Under the motion modes of forward flexion,backward extension,bilateral (left and right) flexion,and left and right rotation,the L3/4 ROM of M2 model and L4/5 ROM of M3 model increased,while the L3/4 and L4/5 ROM of M4 and M5 models significantly decreased compared with those of M1 model.Under all motion modes,the L4 superior endplate in M2 model and the L5 superior endplate in M3 model showed the maximum peak von Mises stress,and the peak von Mises stresses of L4 and L5 superior endplates in M4 and M5 models were close.The L3/4 cage in M2 model and the L4/5 cage in M3 model showcased the largest peak von Mises stress,and the peak von Mises stresses of cages in M4 and M5 models were close.The peak stresses of internal fixation in M2-M5 models were close.Conclusion Four-screw fixation can replace 6-screw fixation in the OLIF surgery for treating two-level degenerative lumbar diseases.

Numerical Study on Mechanical Behavior and Electromechanical Properties of Solder-Jointed REBCO-Coated Conductors.

As the second-generation high-temperature superconducting conductors, rare earth-barium-copper-oxide (REBCO) coated conductor (CC) tapes have good potential as high-field and high-energy superconductors. In superconducting applications, several joints are required for conjugating comparatively short REBCO CC tapes. Soldering lap joints are the simplest and most commonly applied REBCO CC joints. In addition to joint resistance, the mechanical behavior and electromechanical properties are also crucial for superconducting applications. In this paper, the electromechanical properties and mechanical behaviors of soldering lap joints at 77 K under a self-field were studied. The mechanical behavior was addressed by using a full three-dimensional multilayer elastic-plastic finite element model (FEM) with REBCO CC tape main layers and solder connecting layers. Then, the electromechanical properties were analyzed by using Gao's strain-Ic degradation general model on the basis of the FEM results. Both the mechanical behavior and electromechanical properties were verified by experimental results. The effects of soldering lap conditions including lap length, soldering thickness and lap style on the electromechanical properties and mechanical behaviors were discussed. The results indicate that shorter overlap lengths and a thinner solder can reduce the premature degradation of Ic due to stress concentrations nearby the joint edges; moreover, the irreversible critical strain is significantly higher in the back-to-back joint approach compared to the widely used face-to-face joint approach.

Mitigating Measurement Inaccuracies in Digital Twins of Construction Machinery through Multi-Objective Optimization.

The advent of digital twins facilitates the generation of high-fidelity replicas of actual systems or assets, thereby enhancing the design's performance and feasibility. When developing digital twins, precise measurement data is essential to ensure alignment between the actual and digital models. However, inherent uncertainties in sensors and models lead to disparities between observed and predicted (simulated) behaviors. To mitigate these uncertainties, this study originally proposes a multi-objective optimization strategy utilizing a Gaussian process regression surrogate model, which integrates various uncertain parameters, such as load angle, bucket cylinder stroke, arm cylinder stroke, and boom cylinder stroke. This optimization employs a genetic algorithm to indicate the Pareto frontiers regarding the pressure exerted on the boom, arm, and bucket cylinders. Subsequently, TOPSIS is applied to ascertain the optimal candidate among the identified Pareto optima. The findings reveal a substantial congruence between the experimental and numerical outcomes of the devised virtual model, in conjunction with the TOPSIS-derived optimal parameter configuration.

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.

A Finite Element Model for Simulating Stress Responses of Permeable Road Pavement.

Permeable road pavements, due to their open-graded design, suffer from low structural strength, restricting their use in areas with light traffic volume and low bearing capacity. To expand application of permeable road pavements, accurate simulation of stress parameters used in pavement design is essential. A 3D finite element (3D FE) model was developed using ABAQUS/CAE 2021 to simulate pavement stress responses. Utilizing a 53 cm thick permeable road pavement and a 315/80 R22.5 wheel as prototypes, the model was calibrated and validated, with its accuracy confirmed through t-test statistical analysis. Simulations of wheel speeds at 11, 15, and 22 m/s revealed significant impact on pavement depths of 3 cm and 8 cm, while minimal effects were observed at depths of 13 cm and 33 cm. Notably, stress values at a depth of 3 cm with 15 m/s speed in the open-graded asphalt concrete (OGFC) surface layer exceeded those at the speed of 11 m/s, while at a depth of 8 cm in the porous asphalt concrete (PAC) base layer, an opposite performance was observed. This may be attributed to the higher elastic modulus of the OGFC surface layer, which results in different response trends to velocity changes. Overall, lower speeds increase stress responses and prolong action times for both layers, negatively affecting pavement performance. Increasing the moduli of layers is recommended for new permeable road pavements for low-speed traffic. Furthermore, considering the effects of heavy loads and changes in wheel speed, the recommended design depth for permeable road pavement is 30 cm. These conclusions provide a reference for the design of permeable road pavements to address climate change and improve performance.

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.

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.

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.

Finite-Element Modeling of the Temperature Effect on Extended Avalanche Damage of Gas Main Pipelines.

The dynamic stress-strain state and fracture of a steel main gas pipe section between supports with a straight-through crack was analyzed with consideration of the temperature effect on changes in the mechanical properties of the pipe material. The numerical solution of the problem was implemented in the ANSYS-19.2/Explicit Dynamics software package. The process of fracture in a section of the gas pipeline "Beineu-Bozoy-Shymkent" with a linear crack in the temperature range of -40 °C to +50 °C at the operating pressure of 7.5 MPa and critical pressure equal to 9.8 MPa was considered. As a result, it was found that at the initial growth of the internal pressure from working pressure to critical pressure, the length of the crack doubled. At the same time, the process had a local characteristic. Further development of the crack had the nature of avalanche fracture and depended on the temperature of the steel pipeline. With increasing temperature, there was also an increase in the length of the crack at the avalanche fracture. Thus, at a temperature of 40 °C, the crack lengthened 67.75-fold; at a temperature of -10 °C, the crack lengthened 68-fold; at a temperature of +20 °C, the crack lengthened 68.25-fold; and at a temperature of +50 °C, the crack lengthened 68.5-fold. In this work, this difference was 75% of the initial crack length. This fact will be used for further development of the technique of strengthening damaged pipe sections using bandages.

Integrating Experimental and Computational Analyses for Mechanical Characterization of Titanium Carbide/Aluminum Metal Matrix Composites.

This study investigates the mechanical properties of titanium carbide/aluminum metal matrix composites (AMMCs) using both experimental and computational methods. Through accumulative roll bonding (ARB) and cryorolling (CR) processes, AA1050 alloy surfaces were reinforced with TiCp particles to create the Al-TiCp composite. The experimental analysis shows significant improvements in tensile strength, yield strength, elastic modulus, and hardness. The finite element analysis (FEA) simulations, particularly the microstructural modeling of RVE-1 (the experimental case model), align closely with the experimental results observed through scanning electron microscopy (SEM). This validation underscores the accuracy of the computational models in predicting the mechanical behavior under identical experimental conditions. The simulated elastic modulus deviates by 5.49% from the experimental value, while the tensile strength shows a 6.81% difference. Additionally, the simulated yield strength indicates a 2.85% deviation. The simulation data provide insights into the microstructural behavior, stress distribution, and particle-matrix interactions, facilitating the design optimization for enhanced performance. The study also explores the influence of particle shapes and sizes through Representative Volume Element (RVE) models, highlighting nuanced effects on stress-strain behavior. The microstructural evolution is examined via transmission electron microscopy (TEM), revealing insights regarding grain refinement. These findings demonstrate the potential of Al-TiCp composites for lightweight applications.

Improved Finite Element Model Updating of a Highway Viaduct Using Acceleration and Strain Data.

Most finite element model updating (FEMU) studies on bridges are acceleration-based due to their lower cost and ease of use compared to strain- or displacement-based methods, which entail costly experiments and traffic disruptions. This leads to a scarcity of comprehensive studies incorporating strain measurements. This study employed the strain- and acceleration-based FEMU analyses performed on a more than 50-year-old multi-span concrete highway viaduct. Mid-span strains under heavy vehicles were considered for the strain-based FEMU, and frequencies and mode shapes for the acceleration-based FEMU. The analyses were performed separately for up to three variables, representing Young's modulus adjustment factors for different groups of structural elements. FEMU studies considered residual minimisation and the error-domain model falsification (EDMF) methodology. The residual minimisation utilised four different single-objective optimisations focusing on strains, frequencies, and mode shapes. Strain- and frequency-based FEMU analyses resulted in an approximately 20% increase in the overall superstructure's design stiffness. This study shows the benefits of the intuitive EDMF over residual minimisation for FEMU, where information gained from the strain data, in addition to the acceleration data, manifests more sensible updated variables. EDMF finally resulted in a 25-50% overestimated design stiffness of internal main girders.

Voxel Design of Grayscale DLP 3D-Printed Soft Robots.

Grayscale digital light processing (DLP) printing is a simple yet effective way to realize the variation of material properties by tuning the grayscale value. However, there is a lack of available design methods for grayscale DLP 3D-printed structures due to the complexities arising from the voxel-level grayscale distribution, nonlinear material properties, and intricate structures. Inspired by the dexterous motions of natural organisms, a design and fabrication framework for grayscale DLP-printed soft robots is developed by combining a grayscale-dependent hyperelastic constitutive model and a voxel-based finite-element model. The constitutive model establishes the relationship between the projected grayscale value and the nonlinear mechanical properties, while the voxel-based finite-element model enables fast and efficient calculation of the mechanical performances with arbitrarily distributed material properties. A multiphysics modeling and experimental method is developed to validate the homogenization assumption of the degree of conversion (DoC) variation in a single voxel. The design framework is used to design structures with reduced stress concentration and programmable multimodal motions. This work paves the way for integrated design and fabrication of functional structures using grayscale DLP 3D printing.

Quantitative assessment of traumatic brain injury risk in diverse age groups of females: Insights from computational biomechanics.

Traumatic Brain Injury (TBI) stands as a multifaceted health concern, exhibiting varying influences across human population. This study delves into the biomechanical complexities of TBI within gender-specific contexts, focusing on females. Our primary objective is to investigate distinctive injury mechanisms and risks associated with females, emphasizing the imperative for tailored investigations within this cohort. By employing Fluid-Structure Interaction (FSI) Analysis, we conducted simulations to quantify biomechanical responses to traumatic forces across diverse age groups of females. The study utilized a scaling technique to create finite element models (FEMs). The young female FEM, based on anthropometric data, showcased a 15 % smaller head geometry compared to the young male FEM. Moreover, while the elderly female FEM closely mirrored the young female FEM in most structural aspects, it showed distinctive features such as brain atrophy and increased cerebrospinal fluid (CSF) layer thickness. Notably, the child female FEM (ages 7-11 years) replicated around 95 % of the young female FEM's geometry. These structural distinctions meticulously captured age-specific variations across our modeled female age groups. It's noteworthy that identical conditions, encompassing impact intensity, loading type, and boundary conditions, were maintained across all FEMs in this biomechanical finite element analysis, ensuring comparative results. The findings unveiled significant variations in frontal and occipital pressures among diverse age groups, highlighting potential age-related discrepancies in TBI susceptibility among females. These variations were primarily linked to differences in anatomical features, including brain volume, CSF thickness, and brain condition, as the same material properties were used in the FEMs. These results were approximately 4.70, 6.33 and 6.43 % in frontal area of brain in diverse age groups of females (young, elderly, and child) respectively compared to young male FEM. Comparing the FEM results between the young female and the elderly female, we observed a decrease in occipital brain pressure at the same point, reducing from 171,993 to 167,793 Pa, marking an approximate 2.5 % decrease. While typically the elderly exhibit greater brain vulnerability compared to the young, our findings showcase a reduction in brain pressure. Notably, upon assessing the relative movement between the brain and the skull at the point located in occipital area, we observed greater relative movement in the elderly (1.8 mm) compared to the young female (1.04 mm). Therefore, brain atrophy increases the range of motion of the brain within the cranial space. The study underscores the critical necessity for nuanced TBI risk assessment tailored to age and gender, emphasizing the importance of age-specific protective strategies in managing TBIs across diverse demographics. Future research employing individual modeling techniques and exploring a wider age spectrum holds promise in refining our understanding of TBI mechanisms and adopting targeted approaches to mitigate TBI in diverse groups.

Experimental and Numerical Study of Membrane Residual Stress in Q690 High-Strength Steel Welded Box Section Compressed Member.

High-strength steel (HSS) members with welded sections exhibit a notably lower residual compressive stress ratio compared with common mild steel (CMS) members. Despite this difference, current codes often generalize the findings from CMS members to HSS members, and the previous unified residual stress models are generally conservative. This study focuses on the membrane residual stress distribution in Q690 steel welded box sections. By leveraging experimental results, the influence of section sizes and welding parameters on membrane residual stress was delved into. A larger plate size correlates with a decrease in the residual compressive stress across the section, with a more pronounced reduction observed in adjacent plates. Additionally, augmenting the number of welding passes tends to diminish residual stresses across the section. Results showed that membrane residual stress adhered to the section's self-equilibrium, while the self-equilibrium in the plates was not a uniform pattern. A reliable residual stress simulation method for Q690 steel welded box sections was established using a three-dimensional thermal-elastic-plastic finite element model (3DTEFEM) grounded in experimental data. This method served as the cornerstone for parameter analysis in this study and set the stage for subsequent research. As a result, an accurate unified residual stress model for Q690 steel welded box sections was derived.

Computational Evaluation of the Fracture Behavior of Porous Asphalt Concrete Exposed to Moisture and Salt Erosion.

Salt erosion has an adverse impact on the durability of asphalt pavements. Porous asphalt concrete is particularly susceptible to the influence of salt. In this study, a finite element model was developed to investigate the fracture behavior of PAC exposed to salt erosion. The 2D heterogeneous structure of PAC was generated with an image-aided approach to computationally study the fracture behavior of PAC. Laboratory SCB tests were conducted to validate the finite element model. The simulation results of the SCB tests indicate that the peak load of the PAC decreased by 21.8% in dry-wet cycles and 26.1% in freeze-thaw cycles compared to the control group. The salt solution accelerated the degradation of the durability of PAC under both dry-wet cycles and freeze-thaw cycle conditions, which is consistent with laboratory tests. After flushing treatment before the drying phase, the peak load of the PAC in salt environments increased by 5.3% compared to that of the samples with no flushing. Salt erosion also results in a higher average value of scalar stiffness degradation (SDEG), and the damaged elements were primarily the cohesive elements in the fracture of the PAC. Additionally, the influence of crucial factors including the void content, adhesion and cohesion, and loading rate on the fracture behavior of the PAC was analyzed. As the void content increases, the average SDEG value of the cohesive elements increases and surpasses the average SDEG value of the adhesive elements at a void content of approximately 9%. The performance of the fine aggregate matrix (FAM) has a much greater impact than the FAM-aggregate interface on the durability of the PAC. And there were more damaged CZM elements with the increase in the loading rate. Salt erosion results in higher SDEG values and a larger number of cohesive damaged elements at each loading rate.

Design of a simplified cranial substitute with a modal behavior close to that of a human skull.

Individuals exposed to the propagation of shock waves generated by the detonation of explosive charges may suffer Traumatic Brain Injury. The mechanism of cranial deflection is one of many hypotheses that could explain the observed brain damage. To investigate this physical phenomenon in a reproducible manner, a new simplified cranial substitute was designed with a mechanical response close to that of a human skull when subjected to this type of loading. As a first step, a Finite Element Model was employed to dimension the new substitute. The objective was indeed to obtain a vibratory behavior close to that of a dry human skull over a wide range of frequencies up to 10 kHz. As a second step, the Finite Element Model was used together with Experimental Modal Analyses to identify the vibration modes of the substitute. A shaker excited the structure via a metal rod, while a laser vibrometer recorded the induced vibrations at defined measurement points. The results showed that despite differences in material properties and geometry, the newly developed substitute has 10/13 natural frequencies in common with those of dry human skulls. When filled with a simulant of cerebral matter, it could therefore be used in future studies as an approximation to assess the mechanical response of a simplified skull substitute to a blast threat.

Semi-automatic algorithm to build finite element numerical models of the human hearing system from Micro-CT data.

Finite Element modeling has been an extended methodology to build numerical model to simulate the behavior of the hearing system. Due to the complexity of the system and the difficulties to reduce the uncertainties of the geometric data, they result in computationally expensive models, sometimes generic, representative of average geometries. It makes it difficult to validate the model with direct experimental data from the same specimen or to establish a patient-oriented modeling strategy. In the present paper, a first attempt to automatize the process of model building is made. The source information is geometrical information obtained from CT of the different elements that compose the system. Importing that data, we have designed the complete procedure to build a model including tympanic membrane, ossicular chain and cavities. The methodology includes the proper coupling of all the elements and the generation of the corresponding finite element model. The whole automatic procedure is not complete, as we need to make some human-assisted decisions; however, the model development time is reduced from 4 weeks to approximately 3 days. The goal of the modeling algorithm is to build a Finite Element Model with a limited computational cost. Several tasks as contour identification or model decimation are designed and integrated in order to follow a semi-automated process that allows generating a patient-oriented model.

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.