Biomechanical response of human spleen in tensile loading

Investigating the Impact of Blunt Force Trauma: A Probabilistic Study of Behind Armor Blunt Trauma Risk

Evaluating Behind Armor Blunt Trauma (BABT) is a critical step in preventing non-penetrating injuries in military personnel, which can result from the transfer of kinetic energy from projectiles impacting body armor. While the current NIJ Standard-0101.06 standard focuses on preventing excessive armor backface deformation, this standard does not account for the variability in impact location, thorax organ and tissue material properties, and injury thresholds in order to assess potential injury. To address this gap, Finite Element (FE) human body models (HBMs) have been employed to investigate variability in BABT impact conditions by recreating specific cases from survivor databases and generating injury risk curves. However, these deterministic analyses predominantly use models representing the 50th percentile male and do not investigate the uncertainty and variability inherent within the system, thus limiting the generalizability of investigating injury risk over a diverse military population. The DoD-funded I-PREDICT Future Naval Capability (FNC) introduces a probabilistic HBM, which considers uncertainty and variability in tissue material and failure properties, anthropometry, and external loading conditions. This study utilizes the I-PREDICT HBM for BABT simulations for three thoracic impact locations-liver, heart, and lower abdomen. A probabilistic analysis of tissue-level strains resulting from a BABT event is used to determine the probability of achieving a Military Combat Incapacitation Scale (MCIS) for organ-level injuries and the New Injury Severity Score (NISS) is employed for whole-body injury risk evaluations. Organ-level MCIS metrics show that impact at the heart can cause severe injuries to the heart and spleen, whereas impact to the liver can cause rib fractures and major lacerations in the liver. Impact at the lower abdomen can cause lacerations in the spleen. Simulation results indicate that, under current protection standards, the whole-body risk of injury varies between 6 and 98% based on impact location, with the impact at the heart being the most severe, followed by impact at the liver and the lower abdomen. These results suggest that the current body armor protection standards might result in severe injuries in specific locations, but no injuries in others.

Review of non-penetrating ballistic testing techniques for protection assessment: From biological data to numerical and physical surrogates

Human surrogates have long been employed to simulate human behaviour, beginning in the automotive industry and now widely used throughout the safety framework to estimate human injury during and after accidents and impacts. In the specific context of blunt ballistics, various methods have been developed to investigate wound injuries, including tissue simulants such as clays or gelatine ballistic, physical dummies and numerical models. However, all of these surrogate entities must be biofidelic, meaning they must accurately represent the biological properties of the human body. This paper provides an overview of physical and numerical surrogates developed specifically for blunt ballistic impacts, including their properties, use and applications. The focus is on their ability to accurately represent the human body in the context of blunt ballistic impact.

Relationship Between Spleen Pathologic Changes and Spleen Stiffness in Portal Hypertension Rat Model

OBJECTIVE

The aim of the study described here was to explore the influence of splenic pathology and hemodynamic parameters on spleen stiffness in portal hypertension (PH).

METHODS

A Sprague‒Dawley rat model of PH (n = 34) induced by CCl4 was established, and 9 normal rats were used as controls. All animals underwent a routine ultrasound examination, spleen stiffness measurement (SSM), liver stiffness measurement (LSM), portal vein pressure (PVP) measurement and histopathologic assessment. The diagnostic performance of SSM and LSM in PH was evaluated. SSMs were compared among the groups at different pathologic and hemodynamic levels. Multiple linear regression was used to analyze the factors affecting SSM.

RESULTS

SSM had excellent diagnostic efficacy for PH (area under the receiver operating characteristic curve [AUC] = 0.900) and was superior to LSM (AUC = 0.794). In a rat model of PH, pathologic changes such as splenic sinus widening, thickening of the splenic capsule and an increase in collagen fibers were observed in the spleen. There were significant differences in SSM at different splenic capsule thicknesses and splenic sinus widths (all p values <0.05), but there were no significant differences in the SSM at different levels of the splenic collagen fiber area and red pulp area (all p values >0.05). In addition, there were significant differences in SSM at different levels of portal vein diameter, blood flow and congestion index (all p values <0.05). Multiple linear regression analysis revealed that PVP, portal vein congestion index and splenic capsule thickness were significantly associated with SSM.

CONCLUSION

SSM is a good non-invasive way to assess PH. PVP, splenic capsule thickness and portal vein congestion index are responsible for spleen stiffness but not the proliferation of splenic fibrous tissue.

Quasi-Static Mechanical Properties and Continuum Constitutive Model of the Thyroid Gland

The purpose of this study is to obtain the digital twin parameters of the thyroid gland and to build a constitutional model of the thyroid gland based on continuum mechanics, which will lay the foundation for the establishment of a surgical training system for the thyroid surgery robot and the development of the digital twin of the thyroid gland. First, thyroid parenchyma was obtained from fresh porcine thyroid tissue and subjected to quasi-static unconfined uniaxial compression tests using a biomechanical test platform with two strain rates (0.005 s^-1 and 0.05 s^-1) and two loading orientations (perpendicular to the thyroid surface and parallel to the thyroid surface). Based on this, a tensile thyroid model was established to simulate the stretching process by using the finite element method. The thyroid stretching test was carried out under the same parameters to verify the validity of the hyperelastic constitutive model. The quasi-static mechanical property parameters of the thyroid tissue were obtained by a quasi-static unconstrained uniaxial compression test, and a constitutional model that can describe the quasi-static mechanical properties of thyroid tissue was proposed based on the principle of continuum media mechanics, which is of great value for the establishment of a surgical training system for the head and neck surgery robot and for the development of the thyroid digital twin.

Porcine spleen as a model organ for blunt injury impact tests: An experimental and histological study

The spleen is a large and highly vascularized secondary lymphatic organ. Spleen injuries are among the most frequent trauma-related injuries in the abdominal region. The aims of the study were to assess the volume fractions of the main splenic tissue components (red pulp, white pulp, trabeculae and reticular fibres) and to determine the severity of splenic injury due to the experimental impact test. Porcine spleens (n = 17) were compressed by 6.22 kg wooden plate using a drop tower technique from three impact heights (50, 100 and 150 mm corresponding to velocities 0.79, 1.24 and 1.58 m/s). The pressure was measured via catheters placed in the splenic vein. The impact velocity was measured using lasers. The severity of induced injuries was analysed on the macroscopic level. The volume fractions of splenic components were assessed microscopically using stereology. The volume fraction of the red pulp was 76.4%, white pulp 21.3% and trabeculae 2.7% respectively. All impact tests, even with the low impact velocities, led to injuries that occurred mostly in the dorsal extremity of the spleen, and were accompanied by bleeding, capsule rupture and parenchyma crushing. Higher impact height (impact velocity and impact energy) caused more severe injury. Porcine spleen had the same volume fraction of tissue components as human spleen, therefore we concluded that the porcine spleen was a suitable organ model for mechanical experiments. Based on our observations, regions around hilum and the diaphragmatic surface of the dorsal extremity, that contained fissures and notches, were the most prone to injury and required considerable attention during splenic examination after injury. The primary mechanical data are now available for the researchers focused on the splenic trauma modelling.

Mechanical properties of whole-body soft human tissues: a review

The mechanical properties of soft tissues play a key role in studying human injuries and their mitigation strategies. While such properties are indispensable for computational modelling of biological systems, they serve as important references in loading and failure experiments, and also for the development of tissue simulants. To date, experimental studies have measured the mechanical properties of peripheral tissues (e.g. skin)in-vivoand limited internal tissuesex-vivoin cadavers (e.g. brain and the heart). The lack of knowledge on a majority of human tissues inhibit their study for applications ranging from surgical planning, ballistic testing, implantable medical device development, and the assessment of traumatic injuries. The purpose of this work is to overcome such challenges through an extensive review of the literature reporting the mechanical properties of whole-body soft tissues from head to toe. Specifically, the available linear mechanical properties of all human tissues were compiled. Non-linear biomechanical models were also introduced, and the soft human tissues characterized using such models were summarized. The literature gaps identified from this work will help future biomechanical studies on soft human tissue characterization and the development of accurate medical models for the study and mitigation of injuries.

Biomechanical properties of abdominal organs under tension with special reference to increasing strain rate

Currently, abdominal finite element models overlook the organs such as gallbladder, bladder, and intestines, which instead are modeled as a simple bag that is not included in the analysis. Further characterization of the material properties is required for researchers to include these organs into models. This study characterized the mechanical properties of human and porcine gallbladder, bladder, and intestines using uniaxial tension loading from the rates of 25%/s to 500%/s. Small differences were observed between human and porcine gallbladder elastic modulus, failure stress, and failure strain. Strain rate was determined to be a significant factor for predicting porcine gallbladder elastic modulus and failure stress which were found to be 9.03 MPa and 1.83 MPa at 500%/s. Human bladder was observed to be slightly stiffer with a slightly lower failure stress than porcine specimens. Both hosts, however, demonstrated a strain rate dependency with the elastic modulus and failure stress increasing as the rate increased with the highest elastic modulus (2.16 MPa) and failure stress (0.65 MPa) occurring at 500%/s. Both human and porcine intestines were observed to be affected by the strain rate. Failure stress was found to be 1.6 MPa and 1.42 MPa at 500%/s for the human and porcine intestines respectively. For all properties found to be strain rate dependent, a numerical model was created to quantify the impact. These results will enable researchers to create more detailed finite element models that include the gallbladder, bladder, and intestines.

The differences in measured prostate material properties between probing and unconfined compression testing methods

Characterization of the mechanical properties of organs is important for determining their behavior under load and understanding and predicting their response. In order to appropriately understand behavior, including developing predictive models, the method used to measure the properties should match the application as different testing techniques can yield different results. One of the organs where little mechanical testing has been performed is the prostate. Therefore, the goal of this paper is to expand the knowledge of prostate gland mechanical behavior by using two different compressive testing methods under various loading rates. No differences were found between the elastic modulus measured using the compression and probing protocols for both human and porcine specimens. The elastic modulus ranged from 0.08 MPa at 1%/s to 0.24 MPa at 25%/s for human specimens and from 0.2 MPa at 1%/s to 0.4 MPa at 25%/s for porcine specimens. A strain rate dependency of the elastic modulus was observed for both testing methods. The dependency on strain rate started to saturate at higher rates and a material model was created to quantify this dependence as well as the stress-strain behavior. No strain rate dependency was observed for failure stress or failure strain. Overall, similar values of elastic modulus were found for both probing and compression protocols and the relationship developed between elastic modulus and strain rate could be implemented in models of the prostate gland to aid in understanding the response to dynamic loads.

Design framework for mechanically tunable soft biomaterial composites enhanced by modified horseshoe lattice structures

Soft biomaterials have a wide range of applications in many areas. However, one material can only cover a specific range of mechanical performance such as the elastic modulus and stretchability. In order to improve the mechanical performance of soft biomaterials, lattice structures are embedded to reinforce the biomaterials. In this paper, rectangular and triangular lattice structures formed by modified horseshoe microstructures are used because their mechanical properties are tunable and can be tailored precisely to match the desired properties by adjusting four geometrical parameters, the length L, radius R, width w and arc angle θ0. A theoretical design framework for the modified horseshoe lattice structures is developed to predict the dependence of the mechanical behaviors on geometrical parameters. Both experiments and finite element simulations on lattice structures are conducted to validate the theoretical models. Results show that a wide range of design space for the elastic modulus (a few kPa to hundreds of MPa), stretchability (strain up to 180%) and Poisson ratio (ranging from -0.5 to 1.2) can be achieved. Experiments on lattice-hydrogel composites are also conducted to verify the reinforcement effect of lattice structures on the hydrogel. This work provides a theoretical method to predict the mechanical behaviors of the lattice structures and aid the rational design of reinforced biomaterials, which has applications in tissue engineering, drug delivery and intraocular lenses.

The effect of steatosis and fibrosis on blunt force vulnerability of the liver

The aim of our study was to examine the possible effect of steatosis and fibrosis on the blunt force vulnerability of human liver tissue. 3.5 × 3.5 × 2-cm-sized liver tissue blocks were removed from 135 cadavers. All specimens underwent microscopical analysis. The tissue samples were put into a test stand, and a metal rod with a square-shaped head was pushed against the capsular surface. The force (Pmax) causing liver rupture was measured and registered with a Mecmesin AFG-500 force gauge. Six groups were formed according to the histological appearance of the liver tissue: intact (group 1), mild steatosis (group 2), moderate steatosis (group 3), severe steatosis (group 4), fibrosis (group 5), and cirrhosis (group 6). The average Pmax value was 34.1 N in intact liver samples (range from 18.1 to 60.8 N, SD ± 8.7), 45.1 N in mild steatosis (range from 24.2 to 79.8 N SD ± 12.6), 55.4 N in moderate steatosis (range from 28.9 to 92.5 N, SD ± 16.0), 57.6 N in severe steatosis (range from 39.8 to 71.5 N, SD ± 11.9), 63.7 N in fibrosis (range from 37.8 to 112.2 N, SD ± 19.5), and 87.1 N in the case of definite cirrhosis (range from 52.7 to 162.7 N, 30.3). The Pmax values were significantly higher in samples with visible structural change than in intact liver sample (p = 0.023, 0.001, 0.009, 0.0001, 0.0001 between group 1 and groups 2 to 6 respectively). Significant difference was found between mild steatosis (group 2) and cirrhosis (group 6) (p = 0.0001), but the difference between mild, moderate, and severe steatosis (groups 2, 3, and 4) was not significant. Our study demonstrated that contrary to what is expected as received wisdom dictates, the diseases of the parenchyma (steatosis and presence of fibrosis) positively correlate with the blunt force resistance of the liver tissue.

Changes in splenic capsule with aging; beliefs and reality

BACKGROUND

Research describing the splenic capsule and its effect on non-operative management of splenic injuries is limited. The aim of this study is to identify the current beliefs about the splenic capsule thickness and investigate changes in the splenic capsule with age.

METHODS

Trauma Medical Directors were surveyed on their beliefs regarding splenic capsule thickness changes with age. Thicknesses of cadaveric splenic capsule samples were measured.

RESULTS

The majority of trauma medical directors (59%) believe the capsule thickness decreases with age. There were 94 splenic specimens obtained. The splenic capsules of infants were thin and had a uniform layer of elastin fibers. With aging, the capsule becomes thick and develops a collagen layer.

CONCLUSION

Most trauma directors believe the splenic capsule thickness decreases with age. However, our results demonstrate that the splenic capsule thickness increases during childhood but remains constant in adulthood.

Monolithic 3D printing of embeddable and highly stretchable strain sensors using conductive ionogels

Medical training simulations that utilize 3D-printed, patient-specific tissue models improve practitioner and patient understanding of individualized procedures and capacitate pre-operative, patient-specific rehearsals. The impact of these novel constructs in medical training and pre-procedure rehearsals has been limited, however, by the lack of effectively embedded sensors that detect the location, direction, and amplitude of strains applied by the practitioner on the simulated structures. The monolithic fabrication of strain sensors embedded into lifelike tissue models with customizable orientation and placement could address this limitation. The demonstration of 3D printing of an ionogel as a stretchable, piezoresistive strain sensor embedded in an elastomer is presented as a proof-of-concept of this integrated fabrication for the first time. The significant hysteresis and drift inherent to solid-phase piezoresistive composites and the dimensional instability of low-hysteresis piezoresistive liquids inspired the adoption of a 3D-printable piezoresistive ionogel composed of reduced graphene oxide and an ionic liquid. The shear-thinning rheology of the ionogel obviates the need to fabricate additional structures that define or contain the geometry of the sensing channel. Sensors are printed on and subsequently encapsulated in polydimethylsiloxane (PDMS), a thermoset elastomer commonly used for analog tissue models, to demonstrate seamless fabrication. Strain sensors demonstrate geometry- and strain-dependent gauge factors of 0.54-2.41, a high dynamic strain range of 350% that surpasses the failure strain of most dermal and viscus tissue, low hysteresis (<3.5% degree of hysteresis up to 300% strain) and baseline drift, a single-value response, and excellent fatigue stability (5000 stretching cycles). In addition, we fabricate sensors with stencil-printed silver/PDMS electrodes in place of wires to highlight the potential of seamless integration with printed electrodes. The compositional tunability of ionic liquid/graphene-based composites and the shear-thinning rheology of this class of conductive gels endows an expansive combination of customized sensor geometry and performance that can be tailored to patient-specific, high-fidelity, monolithically fabricated tissue models.

How to characterize a nonlinear elastic material? A review on nonlinear constitutive parameters in isotropic finite elasticity

The mechanical response of a homogeneous isotropic linearly elastic material can be fully characterized by two physical constants, the Young’s modulus and the Poisson’s ratio, which can be derived by simple tensile experiments. Any other linear elastic parameter can be obtained from these two constants. By contrast, the physical responses of nonlinear elastic materials are generally described by parameters which are scalar functions of the deformation, and their particular choice is not always clear. Here, we review in a unified theoretical framework several nonlinear constitutive parameters, including the stretch modulus, the shear modulus and the Poisson function, that are defined for homogeneous isotropic hyperelastic materials and are measurable under axial or shear experimental tests. These parameters represent changes in the material properties as the deformation progresses, and can be identified with their linear equivalent when the deformations are small. Universal relations between certain of these parameters are further established, and then used to quantify nonlinear elastic responses in several hyperelastic models for rubber, soft tissue and foams. The general parameters identified here can also be viewed as a flexible basis for coupling elastic responses in multi-scale processes, where an open challenge is the transfer of meaningful information between scales.

Evaluation of kinematics and injuries to restrained occupants in far-side crashes using full-scale vehicle and human body models

OBJECTIVE

The objective of the current study was to perform a parametric study with different impact objects, impact locations, and impact speeds by analyzing occupant kinematics and injury estimations using a whole-vehicle and whole-body finite element-human body model (FE-HBM). To confirm the HBM responses, the biofidelity of the model was validated using data from postmortem human surrogate (PMHS) sled tests.

METHODS

The biofidelity of the model was validated using data from sled experiments and correlational analysis (CORA). Full-scale simulations were performed using a restrained Global Human Body Model Consortium (GHBMC) model seated on a 2001 Ford Taurus model using a far-side lateral impact condition. The driver seat was placed in the center position to represent a nominal initial impact condition. A 3-point seat belt with pretensioner and retractor was used to restrain the GHBMC model. A parametric study was performed using 12 simulations by varying impact locations, impacting object, and impact speed using the full-scale models. In all 12 simulations, the principal direction of force (PDOF) was selected as 90°. The impacting objects were a 10-in.-diameter rigid vertical pole and a movable deformable barrier. The impact location of the pole was at the C-pillar in the first case, at the B-pillar in the second case, and, finally, at the A-pillar in the third case. The vehicle and the GHBMC models were defined an initial velocity of 35 km/h (high speed) and 15 km/h (low speed). Excursion of the head center of gravity (CG), T6, and pelvis were measured from the simulations. In addition, injury risk estimations were performed on head, rib cage, lungs, kidneys, liver, spleen, and pelvis.

RESULTS

The average CORA rating was 0.7. The shoulder belt slipped in B- and C-pillar impacts but somewhat engaged in the A-pillar case. In the B-pillar case, the head contacted the intruding struck-side structures, indicating higher risk of injury. Occupant kinematics depended on interaction with restraints and internal structures-especially the passenger seat. Risk analysis indicated that the head had the highest risk of sustaining an injury in the B-pillar case compared to the other 2 cases. Higher lap belt load (3.4 kN) may correspond to the Abbreviated Injury Scale (AIS) 2 pelvic injury observed in the B-pillar case. Risk of injury to other soft anatomical structures varied with impact configuration and restraint interaction.

CONCLUSION

The average CORA rating was 0.7. In general, the results indicated that the high-speed impacts against the pole resulted in severe injuries, higher excursions followed by low-speed pole, high-speed moving deformable barrier (MDB), and low-speed MDB impacts. The vehicle and occupant kinematics varied with different impact setups and the latter kinematics were likely influenced by restraint effectiveness. Increased restraint engagement increased the injury risk to the corresponding anatomic structure, whereas ineffective restraint engagement increased the occupant excursion, resulting in a direct impact to the struck-side interior structures.

Development of a computationally efficient full human body finite element model

INTRODUCTION

A simplified and computationally efficient human body finite element model is presented. The model complements the Global Human Body Models Consortium (GHBMC) detailed 50th percentile occupant (M50-O) by providing kinematic and kinetic data with a significantly reduced run time using the same body habitus.

METHODS

The simplified occupant model (M50-OS) was developed using the same source geometry as the M50-O. Though some meshed components were preserved, the total element count was reduced by remeshing, homogenizing, or in some cases omitting structures that are explicitly contained in the M50-O. Bones are included as rigid bodies, with the exception of the ribs, which are deformable but were remeshed to a coarser element density than the M50-O. Material models for all deformable components were drawn from the biomechanics literature. Kinematic joints were implemented at major articulations (shoulder, elbow, wrist, hip, knee, and ankle) with moment vs. angle relationships from the literature included for the knee and ankle. The brain of the detailed model was inserted within the skull of the simplified model, and kinematics and strain patterns are compared.

RESULTS

The M50-OS model has 11 contacts and 354,000 elements; in contrast, the M50-O model has 447 contacts and 2.2 million elements. The model can be repositioned without requiring simulation. Thirteen validation and robustness simulations were completed. This included denuded rib compression at 7 discrete sites, 5 rigid body impacts, and one sled simulation. Denuded tests showed a good match to the experimental data of force vs. deflection slopes. The frontal rigid chest impact simulation produced a peak force and deflection within the corridor of 4.63 kN and 31.2%, respectively. Similar results vs. experimental data (peak forces of 5.19 and 8.71 kN) were found for an abdominal bar impact and lateral sled test, respectively. A lateral plate impact at 12 m/s exhibited a peak of roughly 20 kN (due to stiff foam used around the shoulder) but a more biofidelic response immediately afterward, plateauing at 9 kN at 12 ms. Results from a frontal sled simulation showed that reaction forces and kinematic trends matched experimental results well. The robustness test demonstrated that peak femur loads were nearly identical to the M50-O model. Use of the detailed model brain within the simplified model demonstrated a paradigm for using the M50-OS to leverage aspects of the M50-O. Strain patterns for the 2 models showed consistent patterns but greater strains in the detailed model, with deviations thought to be the result of slightly different kinematics between models. The M50-OS with the deformable skull and brain exhibited a run time 4.75 faster than the M50-O on the same hardware.

CONCLUSIONS

The simplified GHBMC model is intended to complement rather than replace the detailed M50-O model. It exhibited, on average, a 35-fold reduction in run time for a set of rigid impacts. The model can be used in a modular fashion with the M50-O and more broadly can be used as a platform for parametric studies or studies focused on specific body regions.

Isotropic incompressible hyperelastic models for modelling the mechanical behaviour of biological tissues: a review

Modelling the mechanical behaviour of biological tissues is of vital importance for clinical applications. It is necessary for surgery simulation, tissue engineering, finite element modelling of soft tissues, etc. The theory of linear elasticity is frequently used to characterise biological tissues; however, the theory of nonlinear elasticity using hyperelastic models, describes accurately the nonlinear tissue response under large strains. The aim of this study is to provide a review of constitutive equations based on the continuum mechanics approach for modelling the rate-independent mechanical behaviour of homogeneous, isotropic and incompressible biological materials. The hyperelastic approach postulates an existence of the strain energy function – a scalar function per unit reference volume, which relates the displacement of the tissue to their corresponding stress values. The most popular form of the strain energy functions as Neo-Hookean, Mooney-Rivlin, Ogden, Yeoh, Fung-Demiray, Veronda-Westmann, Arruda-Boyce, Gent and their modifications are described and discussed considering their ability to analytically characterise the mechanical behaviour of biological tissues. The review provides a complete and detailed analysis of the strain energy functions used for modelling the rate-independent mechanical behaviour of soft biological tissues such as liver, kidney, spleen, brain, breast, etc.

Statistical modeling of human liver incorporating the variations in shape, size, and material properties

The liver is one of the most frequently injured abdominal organs during motor vehicle crashes. Realistic numerical assessments of liver injury risk for the entire occupant population require incorporating inter-subject variations into numerical models. The main objective of this study was to quantify the shape variations of human liver in a seated posture and the statistical distributions of its material properties. Statistical shape analysis was applied to construct shape models of the livers of 15 adult human subjects, recorded in a typical seated (occupant) posture. The principal component analysis was then utilized to obtain the modes of variation, the mean model, and 95% statistical boundary shape models. In addition, a total of 52 tensile tests were performed on the parenchyma of three fresh human livers at four loading rates (0.01, 0.1, 1, and 10 s^-1) to characterize the rate-dependent and failure properties of the human liver. A FE-based optimization approach was employed to identify the material parameters of an Ogden material model for each specimen. The mean material 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 parameters. Results showed that the first five modes of the human liver shape models account for more than 60% of the overall anatomical variations. The distributions of the material parameters combined with the mean and statistical boundary shape models could be used to develop probabilistic finite element (FE) models, which may help to better understand the variability in biomechanical responses and injuries to the abdominal organs under impact loading.

Effect of Strain Rate on the Material Properties of Human Liver Parenchyma in Unconfined Compression

The liver is one of the most frequently injured organs in abdominal trauma. Although motor vehicle collisions are the most common cause of liver injuries, current anthropomorphic test devices are not equipped to predict the risk of sustaining abdominal organ injuries. Consequently, researchers rely on finite element models to assess the potential risk of injury to abdominal organs such as the liver. These models must be validated based on appropriate biomechanical data in order to accurately assess injury risk. This study presents a total of 36 uniaxial unconfined compression tests performed on fresh human liver parenchyma within 48 h of death. Each specimen was tested once to failure at one of four loading rates (0.012, 0.106, 1.036, and 10.708 s−1) in order to investigate the effects of loading rate on the compressive failure properties of human liver parenchyma. The results of this study showed that the response of human liver parenchyma is both nonlinear and rate dependent. Specifically, failure stress significantly increased with increased loading rate, while failure strain significantly decreased with increased loading rate. The failure stress and failure strain for all liver parenchyma specimens ranged from −38.9 kPa to −145.9 kPa and from −0.48 strain to −1.15 strain, respectively. Overall, this study provides novel biomechanical data that can be used in the development of rate dependent material models and the identification of tissue-level tolerance values, which are critical to the validation of finite element models used to assess injury risk.

Shear mechanical properties of the porcine pancreas: experiment and analytical modelling

We provide the first account of the shear mechanical properties of porcine pancreas using a rheometer both in linear oscillatory tests and in constant strain-rate tests reaching the non-linear sub-failure regime. Our results show that pancreas has a low and weakly frequency-dependent dynamic modulus and experiences a noticeable strain-hardening beyond 20% strain. In both linear and non-linear regime, the viscoelastic behaviour of porcine pancreas follows a four-parameter bi-power model that has been validated on kidney, liver and spleen. Among the four solid organs of the abdomen, pancreas proves to be the most compliant and the most viscous one.

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.

Experimental mechanical characterization of abdominal organs: liver, kidney & spleen

Abdominal organs are the most vulnerable body parts during vehicle trauma, leading to high mortality rate due to acute injuries of liver, kidney, spleen and other abdominal organs. Accurate mechanical properties and FE models of these organs are required for simulating the traumas, so that better designing of the accident environment can be done and the organs can be protected from severe damage. Also from biomedical aspect, accurate mechanical properties of organs are required for better designing of surgical tools and virtual surgery environments. In this study porcine liver, kidney and spleen tissues are studied in vitro and hyper-elastic material laws are provided for each. 12 porcine kidneys are used to perform 40 elongation tests on renal capsule and 60 compression tests on renal cortex, 5 porcine livers are used to perform 45 static compression tests on liver parenchyma and 5 porcine spleens are used to carry out 20 compression tests. All the tests are carried out at a static speed of 0.05 mm/s. A comparative analysis of all the results is done with the literature and though the results are of same order of magnitude, a slight dissonance is observed for the renal capsule. It is also observed that the spleen is the least stiff organ in the abdomen whereas the kidney is the stiffest. The results of this study would be essential to develop the FE models of liver, kidney and spleen which can be further used for impact biomechanical and biomedical applications.