Development of a 10-year-old paediatric thorax finite element model validated against cardiopulmonary resuscitation data

A fast methodology for generating skeletal FEM with detailed human geometric features based on CPD and RBF algorithms

Due to the significant effects of the human anatomical characteristics on the injury mechanism of passenger in traffic accidents, it is necessary to develop human body FEM (Finite Element Model) with detailed anatomical characteristics. However, traditional development of a human body FEM is an extremely complicated process. In particular, the meshing of human body is a huge and time-consuming project. In this paper, a new fast methodology based on CPD (Coherent Point Drift) and RBF (Radial Basis Function) was proposed to achieve the rapid developing the FEM of human bone with detailed anatomical characteristics. In this methodology, the mesh morphing technology based the RBF was used to generate FEM mesh in the geometry extracted from the target CT (Computed Tomography) data. In order to further improve the accuracy and speed of mesh morphing, the target geometric feature points required in the mesh morphing process were realized via the rapid and automatic generation based on the point-cloud registration technology of the CPD algorithm. Finally, this new methodology was used to generate a 3-year-old ribcage FEM consisting of a total of 27,728 elements with mesh size 3-5 mm based on the THUMS (Total Human Model for Safety) adult model. In the entire process of generating this new ribcage model, it only took about 2.7 s. The average error between the new FEM and target geometries was only about 2.7 mm. This indicated that the new FEM well described the detailed anatomical characteristics of target geometry, thus importantly revealing that the mesh quality of the new FEM was basically similar to that of source FEM.

Identifying Vulnerable Impact Locations to Reduce the Occurrence of Deadly Commotio Cordis Events in Children's Baseball: A Computational Approach

Commotio cordis is the second leading cause of sudden cardiac death in young athletes. Currently available chest protectors on the market are ineffective in preventing cases of commotio cordis in young athletes who play baseball. This study focused on using contour maps to identify specific baseball impact locations to the chest that may result in instances of commotio cordis to children during baseball games. By identifying these vulnerable locations, we may design and develop chest protectors that can provide maximum protection to prevent commotio cordis in young athletes. Simulation cases were run using the validated CHARM-10 chest model, a detailed finite element model representing an average 10-year-old child's chest. A baseball model was developed in company with the chest model, and then used to impact the chest at different locations. A 7 × 8 impact location matrix was designed with 56 unique baseball impact simulations. Left ventricle strain and pressure, reaction force between the baseball and chest, and rib deformations were analyzed. Left ventricle strain was highest from baseball impacts directly over the left ventricle (0.34) as well as impacts slightly lateral and superior to the cardiac silhouette (0.34). Left ventricle pressure was highest with impacts directly over the left ventricle (82.94 kPa). We have identified the most dangerous impact locations resulting in high left ventricle strain and pressure. This novel study provided evidence of where to emphasize protective materials for establishing effective chest protectors that will minimize instances of commotio cordis in young athletes.

Development of a Finite Element Model of the Pediatric Thoracic and Lumbar Spine, Ribcage, and Pelvis with Orthotropic Region-Specific Vertebral Growth

Finite element (FE) modeling of the spine has increasingly been applied in orthopedic precision-medicine approaches. Previously published FE models of the pediatric spine growth have made simplifications in geometry of anatomical structures, material properties, and representation of vertebral growth. To address those limitations, a comprehensive FE model of a pediatric (10-year-old) osteo-ligamentous thoracic and lumbar spine (T1-L5 with intervertebral discs (IVDs) and ligaments), ribcage, and pelvis with age- and level-specific ligament properties and orthotropic region-specific vertebral growth was developed and validated. Range of motion (ROM) measures, namely lateral bending, flexion-extension, and axial rotation, of the current 10 YO FE model were generally within reported ranges of scaled in vitro adult ROM data. Changes in T1-L5 spine height, as well as kyphosis (T2-T12) and lordosis (L1-L5) angles in the current FE model for two years of growth (from ages 10 to 12 years) were within ranges reported from corresponding pediatric clinical data. The use of such comprehensive pediatric FE models can provide clinically relevant insights into normative and pathological biomechanical responses of the spine, and also contribute to the development and optimization of clinical interventions for spine deformities.

Developing commotio cordis injury metrics for baseball safety: unravelling the connection between chest force and rib deformation to left ventricle strain and pressure

Commotio cordis is a sudden death mechanism that occurs when the heart is impacted during the repolarization phase of the cardiac cycle. This study aimed to investigate commotio cordis injury metrics by correlating chest force and rib deformation to left ventricle strain and pressure. We simulated 128 chest impacts using a simulation matrix which included two initial velocities, 16 impact locations spread across the transverse and sagittal plane, and four baseball stiffness levels. Results showed that an initial velocity of 17.88 m/s and an impact location over the left ventricle was the most damaging setting across all possible settings, causing the most considerable left ventricle strain and pressure increases. The impact force metric did not correlate with left ventricle strain and pressure, while rib deformations located over the left ventricle were strongly correlated to left ventricle strain and pressure. These results lead us to the recommendation of exploring new injury metrics such as the rib deformations we have highlighted for future commotio cordis safety regulations.

[Mechanism of blunt spleen injury: a finite element analysis]

OBJECTIVE

To explore the biomechanical mechanism of blunt spleen injury based on finite element analysis.

OBJECTIVE

A fist finite element model was used to simulate the impact at 4-8 m/s in the spleen area of THUMS4.0 human body model from the front of the left costal area, the left anterior axillary line and the rear scapular line. The strain distribution and damage of the spleen under different conditions were observed. The simulation results were compared with the clinical cases of spleen rupture to analyze the mechanism of spleen injury.

OBJECTIVE

The damage location and strain distribution of the spleen could vary under different conditions. Due to the special anatomical location of the spleen, a blunt impact at the speed of 4-8 m/s on the front side did not easily cause spleen injury, and the strain was distributed mainly in the front of the spleen and the spleen hilum; a similar blunt impact on the left side was likely to cause spleen diaphragmatic surface injury, the splenic visceral surface could be injured by the compression of the medial tissue and organs and the traction of the splenic pedicle, and the strain was distributed in the spleen diaphragmatic and visceral surfaces; an impact on the back side was likely to cause injuries in the posterior portion and hilum of the spleen, and the strain was mainly concentrated in the injured area.

OBJECTIVE

Blunt spleen injuries caused by punches on the abdomen are mostly caused by direct impact on the ribs, the compression by the surrounding tissues and organs and the traction by the spleen pedicle.

Study on the Long Bone Failure Behaviors Under the Indenter Rigid-Contact by Experiment Analysis and Subject-Specific Simulation

The bending fracture behaviors of long bone have gained great attention due to the high bending fracture risk during sports events, traffic accidents, and falling incidents, etc. For evaluating bone bending behaviors, most of the previous studies used an indenter in three point bending experiments while the effect of its rigidity was never considered. In this work, using the porcine long bones, the three point bending tests were conducted to explore the bone fracture behaviors under a rigid indenter. In addition to collecting the force applied, the bone fracture dynamic process was recorded by high-speed photography, and the fracture surface profile in mesoscale was observed by the scanning electron microscope (SEM). Based on CT scanning of long bones, the cross section properties of test specimens were calculated by a homemade matlab script for correlating with their failure strengths. Also, a subject-specific finite element (FE) model was developed to identify the outcomes induced by a rigid indenter on simulation. Findings led to conclusions as follows: (1) The tension fracture came with fracture path deflection, which was caused by the bone indentation induced mesoscale crack-opening. Due to this damage before the whole bone fracture, a bone fracture moment correction was established to compensate experimental data. (2) The plastic indentation caused the force fluctuation as suggested by correlation analysis. (3) The bone failure moment correlated with the inertial moment of the bone cross section at the fracture location higher than the traditional cross section area. (4) In the subject-specific simulation, the indentation caused compression fracture under a much lower failure force. Removing the element erosion on the indenter-contacted area only during the validation was verified as a good option to solve this issue.

Validation of a finite element model with six-year-old child anatomical characteristics as specified in Euro NCAP Pedestrian Human Model Certification (TB024)

Accident statistics show that more than 80% of car-to-pedestrian collisions (CPC) occur when pedestrians cross the road. It is very important to establish a finite element model with natural walking posture to study the kinematics and injury mechanism of pedestrians. In this study, a finite element model of six-year-old child pedestrian is developed with detailed anatomical characteristics and posture parameters as specified in Euro NCAP Pedestrian Human Model Certification (TB024). The numerical human body model is validated in total twelve simulations in which the pedestrian is impacted against four generic vehicle models at speeds 30, 40, 50 km/h prescribed in TB024. The Head Impact Time (HIT), Contact Force and the Trajectories of HC, T12 and AC of all twelve simulations are compared with the reference corridors provided by Technical Bulletin 024. The results indicate that the numerical human body model of a six-year-old child can be used to demonstrate the suitability of the sensing system for the range of pedestrian sizes; the timing of system deployment, and the bonnet deflection due to body loading. Furthermore, the model could be a good tool for further research on pedestrian injury mechanism and the development of pedestrian protection devices.

Development and validation of a semi-automatic landmark extraction method for mesh morphing

Currently, landmark-based mesh morphing technology is widely used to rapidly obtain meshes with specific geometry, which is suitable to develop parametric human finite element (FE) models. However it takes too much time for landmark extraction to obtain high geometric accuracy. The purpose of this study is to develop and validate a semi-automatic landmark extraction method to reduce the time of manual selection of landmarks without sacrificing the accuracy of identifying landmarks in the process of mesh morphing. A few contour edge landmarks were extracted manually. Mathematical landmarks and pseudo-landmarks were extracted automatically by user-defined algorithm. The radial basis function (RBF) was used to morph the baseline FE model into the target geometry based on these landmarks. The cervical vertebra (C5), rib (R7) and femur were selected as the target geometries to verify the effectiveness of the method. The maximum mean geometric error of the three types of target geometries was less than 1 mm. The mesh quality of the morphed FE model was similar to that of the baseline FE model. Compared to the traditional manual method, 2/3 to 3/4 of the time for landmark extraction was saved by the semi-automatic method.

Morphing the feature-based multi-blocks of normative/healthy vertebral geometries to scoliosis vertebral geometries: development of personalized finite element models

Personalized Finite Element (FE) models and hexahedral elements are preferred for biomechanical investigations. Feature-based multi-block methods are used to develop anatomically accurate personalized FE models with hexahedral mesh. It is tedious to manually construct multi-blocks for large number of geometries on an individual basis to develop personalized FE models. Mesh-morphing method mitigates the aforementioned tediousness in meshing personalized geometries every time, but leads to element warping and loss of geometrical data. Such issues increase in magnitude when normative spine FE model is morphed to scoliosis-affected spinal geometry. The only way to bypass the issue of hex-mesh distortion or loss of geometry as a result of morphing is to rely on manually constructing the multi-blocks for scoliosis-affected spine geometry of each individual, which is time intensive. A method to semi-automate the construction of multi-blocks on the geometry of scoliosis vertebrae from the existing multi-blocks of normative vertebrae is demonstrated in this paper. High-quality hexahedral elements were generated on the scoliosis vertebrae from the morphed multi-blocks of normative vertebrae. Time taken was 3 months to construct the multi-blocks for normative spine and less than a day for scoliosis. Efforts taken to construct multi-blocks on personalized scoliosis spinal geometries are significantly reduced by morphing existing multi-blocks.

Numerical simulations of the 10-year-old head response in drop impacts and compression tests

BACKGROUND AND OBJECTIVE

Studies on traumatic injuries of children indicate that impact to the head is a major cause of severe injury and high mortality. However, regulatory and ethical concerns very much limit development and validation of computer models representing the pediatric head. The purpose of this study was to develop a child head finite element model with high-biofidelity to be used for studying pediatric head injury mechanisms.

METHODS

A newly developed 10-year-old (YO) pediatric finite element head model was limitedly validated for kinematic and kinetic responses against data from quasi-static compressions and drop tests obtained from an experimental study involving a child-cadaver specimen. The validated model was subsequently used for a fall accident reconstruction and associated injury analysis.

RESULTS

The model predicted the same shape of acceleration-time histories as was found in drop tests with the largest discrepancy of -8.2% in the peak acceleration at a drop height of 15 cm. Force-deflection responses predicted by the model for compression loading had a maximum discrepancy of 7.5% at a strain rate of 0.3 s(-1). The model-predicted maximum von Mises stress (σv) and principal strain (εp) in the skull, intracranial pressure (ICP), maximum σv and maximum εp in the brain, head injury criterion (HIC), brain injury criterion (BrIC), and head impact power (HIP) were used for analyzing risks of injury in the accident reconstruction.

CONCLUSIONS

Based on the results of the injury analyses, the following conclusions can be drawn: (1) ICP cannot be used to accurately predict the locations of brain injury, but it may reflect the overall energy level of the impact event. (2) The brain regions predicted by the model to have high σv coincide with the locations of subdural hematoma with transtentorial herniation and the impact position of an actual injury. (3) The brain regions with high εp predicted by the model coincide with locations commonly found where diffuse axonal injuries (DAI) due to blunt-impact and rapid acceleration have taken place.

Computational Modeling of Traffic Related Thoracic Injury of a 10-Year-Old Child Using Subject-Specific Modeling Technique

Traffic injuries have become a major health-related issue to school-aged children. To study this type of injury with numerical simulations, a finite element model was developed to represent the full body of a 10-year-old (YO) child. The model has been validated against test data at both body-part and full-body levels in previous studies. Representing only the average 10-YO child, this model did not include subject-specific attributes, such as the variations in size and shape among different children. In this paper, a new modeling approach was used to morph this baseline model to a subject-specific model, based on anthropometric data collected from pediatric subjects. This mesh-morphing method was then used to rapidly morph the baseline mesh into the subject-specific geometry while maintaining a good mesh quality. The morphed model was subsequently applied to simulate a real-world motor vehicle crash accident. A lung injury observed in the accident was well captured by the subject-specific model. The findings of this study demonstrate the feasibility of the proposed morphing approach to develop subject-specific human models, and confirm their capability in prediction of traffic injuries.