Yuan-Chiao Lu

Material characterization of liver parenchyma using specimen-specific finite element models

The liver is one of the most frequently injured abdominal organs during motor vehicle crashes. Realistic car crash simulations require incorporating strain-rate dependent mechanical properties of soft tissue in finite element (FE) material models. This study presents a total of 30 tension tests performed on fresh bovine liver parenchyma at various loading rates in order to characterize the biomechanical and failure properties of liver parenchyma. Each specimen, cut in a standard dog-bone shape, was tested until failure at one of three loading rates (0.01 s(-1), 0.1s(-1), 1 s(-1)) using a tensile testing setup. Load and acceleration recorded from each specimen grip were employed to calculate the time history of force at specimen ends. The shapes of all specimens were reconstructed from laser scans recorded prior to each test and then used to develop specimen-specific FE models. A first-order Ogden material model and the time histories of specimen end displacement were assigned to each specimen FE model. The failure Green-Lagrangian strain showed averages around 50% and no significant dependence on loading rates, but the failure 2nd Piola-Kirchhoff stress showed rate-dependence with average values ranging from 33 kPa to 94 kPa. The FE models with material model parameters identified using a simulation-based optimization replicated well the time history of load recorded during the test. The FE simulations with model parameters identified using an analytical approach or based on the displacement of optical markers showed a significantly stiffer response and lower failure stress/strain than the FE specimen-specific models. This study provides novel biomechanical and failure data which can be easily implemented in FE models and used to assess injury risk in automobile collisions.

Statistical shape analysis of clavicular cortical bone with applications to the development of mean and boundary shape models

During car collisions, the shoulder belt exposes the occupants clavicle to large loading conditions which often leads to a bone fracture. To better understand the geometric variability of clavicular cortical bone which may influence its injury tolerance, twenty human clavicles were evaluated using statistical shape analysis. The interior and exterior clavicular cortical bone surfaces were reconstructed from CT-scan images. Registration between one selected template and the remaining 19 clavicle models was conducted to remove translation and rotation differences. The correspondences of landmarks between the models were then established using coordinates and surface normals. Three registration methods were compared: the LM-ICP method; the global method; and the SHREC method. The LM-ICP registration method showed better performance than the global and SHREC registration methods, in terms of compactness, generalization, and specificity. The first four principal components obtained by using the LM-ICP registration method account for 61% and 67% of the overall anatomical variation for the exterior and interior cortical bone shapes, respectively. The length was found to be the most significant variation mode of the human clavicle. The mean and two boundary shape models were created using the four most significant principal components to investigate the size and shape variation of clavicular cortical bone. In the future, boundary shape models could be used to develop probabilistic finite element models which may help to better understand the variability in biomechanical responses and injuries to the clavicle.

Effect of seat belt pretensioners on human abdomen and thorax: Biomechanical response and risk of injuries.

BACKGROUND: A better coupling of the occupant to the car seat in the early phase of a frontal or far side impacts using pretensioner systems may reduce the likelihood of the submarining effect or increases the likelihood of seat belt engaging the shoulder, respectively. However, the high belt forces may also increase the risk of upper body injuries to the vehicle occupant (especially in abdominal region). It was hypothesized that human body characteristics, such as body mass index (BMI) influence the biomechanical response and injury outcome to the abdominal regions during static pretensioning deployment tests. METHODS: Four postmortem human specimens (PMHS), in a BMI range from 15.6 to 31.2, were positioned in production seats in a normal passenger position and were restrained using a standard three-point belt system. The pretension forces in the belts were generated at two points (retractor and right anchorage) or at all three locations (retractor, left anchorage, and right anchorage). An optical motion capture system and acceleration cubes mounted to the lumbar spine were used to measure the abdomen deformation during testing. RESULTS: The normalized deflections of the thorax recorded at the level of fourth rib were under 10% (noninjury level). Two different patterns were observed in the time histories of abdominal penetration rate in the four PMHSs associated with lower and higher BMI. Abdominal injuries (spleen lacerations) were observed only in the two PMHS with highest BMI. CONCLUSION: Based on data from this study and similar data from the literature, belt velocity and FmaxCmax were shown to be the best injury predictors for injury risk analysis for Abbreviated Injury Scale 2+ and for Abbreviated Injury Scale 3+ injuries, respectively.

Effect of storage on tensile material properties of bovine liver.

Cadaveric tissue models play an important role in the assessment and optimization of novel restraint systems for reducing abdominal injuries. However, the effect of tissue preservation by means of freezing on the material properties of abdominal tissues remains unknown. The goal of this study was to investigate the influence of frozen storage time on the material responses of the liver parenchyma in tensile loading. Specimens from ten bovine livers were equally divided into three groups: fresh, 30-day frozen storage, and 60-day frozen storage. All preserved specimens were stored at -12°C. Dog-bone specimens from each preservation group were randomly assigned to one of three strain rates (0.01s(-1), 0.1s(-1), and 1.0s(-1)) and tested to failure in tensile loading. The local material response recorded at the tear location and the global material response of the whole specimen of the liver parenchyma specimens were investigated based on the experimental data and optimized analytical material models. The local and global failure strains decreased significantly between fresh specimens and specimens preserved for 30 days (p<0.05), and between fresh specimens and specimens preserved for 60 days (p<0.05) for all three loading rates. Changes on the material model parameters were also observed between fresh and preserved specimens. Preservation by means of frozen storage was found to affect both the material and failure response of bovine liver parenchyma in tensile loading. The stiffness of the tissue increased with increased preservation time and increased strain rate. In summary, significant changes (p<0.05) between the failure strain of previously frozen liver parenchyma samples and fresh samples were demonstrated at both global and local levels in this study. In addition, nonlinear and viscoelastic characteristics of the liver parenchyma were observed in tension for both fresh and preserved samples.

Modeling the biomechanical and injury response of human liver parenchyma under tensile loading

The rapid advancement in computational power has made human finite element (FE) models one of the most efficient tools for assessing the risk of abdominal injuries in a crash event. In this study, specimen-specific FE models were employed to quantify material and failure properties of human liver parenchyma using a FE optimization approach. Uniaxial tensile tests were performed on 34 parenchyma coupon specimens prepared from two fresh human livers. Each specimen was tested to failure at one of four loading rates (0.01s(-1), 0.1s(-1), 1s(-1), and 10s(-1)) to investigate the effects of rate dependency on the biomechanical and failure response of liver parenchyma. Each test was simulated by prescribing the end displacements of specimen-specific FE models based on the corresponding test data. The parameters of a first-order Ogden material model were identified for each specimen by a FE optimization approach while simulating the pre-tear loading region. The mean material model parameters were then determined for each loading rate from the characteristic averages of the stress-strain curves, and a stochastic optimization approach was utilized to determine the standard deviations of the material model parameters. A hyperelastic material model using a tabulated formulation for rate effects showed good predictions in terms of tensile material properties of human liver parenchyma. Furthermore, the tissue tearing was numerically simulated using a cohesive zone modeling (CZM) approach. A layer of cohesive elements was added at the failure location, and the CZM parameters were identified by fitting the post-tear force-time history recorded in each test. The results show that the proposed approach is able to capture both the biomechanical and failure response, and accurately model the overall force-deflection response of liver parenchyma over a large range of tensile loadings rates.

Effect of storage methods on indentation-based material properties of abdominal organs

To investigate the possible changes in material properties of cadaveric abdominal organs due to the preservation methods, the indentation data obtained from porcine abdominal organs (kidney, liver, and spleen) preserved by cooling and freezing are analyzed statistically in this study. Indentation tests were first conducted on fresh specimens. One half of the specimens of each organ were then frozen (preserved at −12 °C), and the other half of the specimens were cooled (preserved at 4 °C). All preserved specimens were retested after 20 days. Force and displacement data recorded during indentation were analyzed using a quasi-linear viscoelastic model. The results show that both cooling and freezing storage increased the kidney stiffness. In contrast, both storage methods decreased the stiffness of the spleen specimens. While cooling increased the liver stiffness, no significant changes of the instantaneous elastic response were observed in the liver specimens preserved by freezing. The liver and spleen’s reduced relaxation responses and the liver’s instantaneous elastic response were significantly different when comparing between cooling and freezing effects after 20 days of preservation. This study showed that both cooling and freezing storage methods significantly changed the material properties of abdominal organs, especially the instantaneous elastic response. More research is needed in investigating the effect of preservation on failure properties and mechanical properties under large deformation.

A statistical geometrical description of the human liver for probabilistic occupant models

Realistic numerical assessments of liver injury risk for the entire occupant population require incorporating inter-subject variations into numerical models. Statistical shape models of the abdominal organs have been shown to be useful tools for the investigation of the organ variations and could be applied to the development of statistical computational models. The main objective of this study was to establish a standard procedure to quantify the shape variations of a human liver in a seated posture, and construct three-dimensional (3D) statistical shape boundary models. Statistical shape analysis was applied to construct shape models of 15 adult human livers. Principal component analysis (PCA) was then utilized to obtain the modes of variation, the mean model, and a set of statistical boundary shape models, which were constructed using the q-hyper-ellipsoid approach. The first five modes of a human liver accounted for the major anatomical variations. The modes were highly correlated to the height, thickness, width, and curvature of the liver, and the concavity of the right lobe. The mean model and the principal components were utilized to construct four boundary models of human liver. The statistical boundary model approach presented in this study could be used to develop probabilistic finite element (FE) models. In the future, the probabilistic liver models could be used in FE simulations to better understand the variability in biomechanical responses and abdominal injuries under impact loading.

Statistical shape analysis of the human spleen geometry for probabilistic occupant models.

Statistical shape models are an effective way to create computational models of human organs that can incorporate inter-subject geometrical variation. The main objective of this study was to create statistical mean and boundary models of the human spleen in an occupant posture. Principal component analysis was applied to fifteen human spleens in order to find the statistical modes of variation, mean shape, and boundary models. A landmark sliding approach was utilized to refine the landmarks to obtain a better shape correspondence and create a better representation of the underlying shape contour. The first mode of variation was found to be the overall volume, and it accounted for 69% of the total variation. The mean model and boundary models could be used to develop probabilistic finite element (FE) models which may identify the risk of spleen injury during vehicle collisions and consequently help to improve automobile safety systems.

A bootstrap approach for lower injury levels of the risk curves

Survival analysis is widely applied to develop injury risk curves from biomechanical data. To obtain more accurate estimation of confidence intervals of parameters, bootstrap method was evaluated by a designed simulation process. Four censoring schemes and various sample sizes were considered to investigate failure time parameters corresponding to low-level injury probabilities. In the numerical simulations, the confidence interval ranges developed by bootstrapping were about two-third of the corresponding ranges calculated by asymptotical normal approximation and showed highest reduction for censored datasets with smaller sample size (≤ 40). In analysis of two experimental datasets with reduced sample sizes and mixed censored data, it was shown that the bootstrapping reduce significantly the confidence intervals as well. The results presented in this study recommend using bootstrapping in development of more accurate confidence intervals for risk curves in injury biomechanics, which consequently will lead to better regulations and safer vehicle designs.