Material characterization of the pig kidney in relation with the biomechanical analysis of renal trauma

The objective of this study was an investigation of the material properties of the fresh pig kidney and parametric characterization of its elastic and inelastic material behavior. The material investigation included density measurements, uniaxial as well as three-dimensional compression tests, tensile tests. and shear tests on the samples extracted from the fresh pig kidney. For comparison, density measurements on a number of soft synthetic materials were also performed. Compression tests on the radial and the tangential specimens from the cortex tissue were performed at various loading rates. Three-axial compression tests were performed on the cortex tissues placed in a compression chamber. Shear tests were performed by punching a cylinder into a slice of the cortex. Tensile tests were carried out on the outer capsule. For characterization of the material behavior, a non-linear theoretical simulation based on a two parameter Blatz model was used. For characterization of the time-dependent behavior of the pig kidney cortex, a four-parameter linear viscoelastic model was employed. From the present experimental and theoretical studies, a number of conclusions were derived: (1) The general behavior of the pig kidney cortex samples under compression showed the general non-linear features typical of the soft tissues; the stress strain diagram was composed of a very flat part at very low stress level to about 30% relative deformation which was followed by a steeply rising stiffening leading to the radial rupture of samples marked by a maximum nominal rupture strain of about 50%. (2) The uniaxial compression tests on the radial and the tangential samples from the cortex tissue showed an increase of the rupture stress with the increase in the loading rate, but a decrease in the related rupture strain. (3) The long-term uniaxial compression tests on the cortex specimens under sustained constant load showed an instantaneous deformation followed by a creep response which eventually approached an asymptote. (4) Simulation of the non-linear material behavior of the cortex tissue under uniaxial compression by the Blatz model gave two pairs of material parameters for the cortex in the radial and the tangential directions. Furthermore, fitting of the assumed four-parameter linear viscoelastic model with the experimental data resulted in the viscoelastic material parameters.     

A constitutive formulation of vascular tissue mechanics including viscoelasticity and softening behaviour

Nearly all soft tissues, among which the vascular tissue is included, present a certain degree of viscoelastic response. This behaviour may be attributed in part to fluid transport within the solid matrix, and to the friction between its fluid and solid constituents. After being preconditioned, the tissue displays highly repetitive behaviour, so that it can be considered pseudo-elastic, that is, elastic but behaving differently in loading and unloading. Because of this reason, very few constitutive laws accounting for the viscoelastic behaviour of the tissue have been developed. Nevertheless, the consideration of this inelastic effect is of crucial importance in surgeries—like vascular angioplasty—where the mentioned preconditioning cannot be considered since non-physiological deformation is applied on the vessel which, in addition, can cause damage to the tissue. A new constitutive formulation considering the particular features of the vascular tissue, such as anisotropy, together with these two inelastic phenomena is presented here and used to fit experimental stress–stretch curves from simple tension loading–unloading tests and relaxation test on porcine and ovine vascular samples.     

A visco-hyperelastic constitutive approach for modeling polyvinyl alcohol sponge

This study proposes the quasi-linear viscoelastic (QLV) model to characterize the time dependent mechanical behavior of poly(vinyl alcohol) (PVA) sponges. The PVA sponges have implications in many viscoelastic soft tissues, including cartilage, liver, and kidney as an implant. However, a critical barrier to the use of the PVA sponge as tissue replacement material is a lack of sufficient study on its viscoelastic mechanical properties. In this study, the nonlinear mechanical behavior of a fabricated PVA sponge is investigated experimentally and computationally using relaxation and stress failure tests as well as finite element (FE) modeling. Hyperelastic strain energy density functions, such as Yeoh and Neo-Hookean, are used to capture the mechanical behavior of PVA sponge at ramp part, and viscoelastic model is used to describe the viscose behavior at hold part. Hyperelastic material constants are obtained and their general prediction ability is verified using FE simulations of PVA tensile experiments. The results of relaxation and stress failure tests revealed that Yeoh material model can define the mechanical behavior of PVA sponge properly compared with Neo-Hookean one. FE modeling results are also affirmed the appropriateness of Yeoh model to characterize the mechanical behavior of PVA sponge. Thus, the Yeoh model can be used in future biomechanical simulations of the spongy biomaterials. These results can be utilized to understand the viscoelastic behavior of PVA sponges and has implications for tissue engineering as scaffold.     

Nonlinear viscoelastic material property estimation of lower extremity residual limb tissues

Axisymmetric nonlinear finite-element analysis was used to simulate force-relaxation and creep data obtained during in vivo indentation of the residual limb soft tissues of six individuals with trans-tibial amputation [1]. The finite-element models facilitated estimation of an appropriate set of nonlinear viscoelastic material coefficients of extended James-Green-Simpson material formulation for bulk soft tissue at discrete, clinically relevant test locations. The results indicate that over 90% of the experimental data can be simulated using the two-term viscoelastic Prony series extension of James-Green-Simpson material formulation. This phenomenological material formulation could not, however, predict the creep response from relaxation experiments, nor the relaxation response from creep experiments [2-5]. The estimated material coefficients varied with test location and subject indicating that these coefficients cannot be readily extrapolated to other sites or individuals.     

A strain-hardening bi-power law for the nonlinear behaviour of biological soft tissues

Biological soft tissues exhibit a strongly nonlinear viscoelastic behaviour. Among parenchymous tissues, kidney and liver remain less studied than brain, and a first goal of this study is to report additional material properties of kidney and liver tissues in oscillatory shear and constant shear rate tests. Results show that the liver tissue is more compliant but more strain hardening than kidney. A wealth of multi-parameter mathematical models has been proposed for describing the mechanical behaviour of soft tissues. A second purpose of this work is to develop a new constitutive law capable of predicting our experimental data in the both linear and nonlinear viscoelastic regime with as few parameters as possible. We propose a nonlinear strain-hardening fractional derivative model in which six parameters allow fitting the viscoelastic behaviour of kidney and liver tissues for strains ranging from 0.01 to 1 and strain rates from 0.0151 s^?1 to 0.7 s^?1.     

Simulation of the aging face

A three-dimensional finite element program is described which attempts to simulate the nonlinear mechanical behavior of an aging human face with specific reference to progressive gravimetric soft tissue descent. A cross section of the facial structure is considered to consist of a multilayered composite of tissues with differing mechanical behavior. Relatively short time (elastic-viscoplastic) behavior is governed by equations previously developed which are consistent with mechanical tests. The long time response is controlled by the aging elastic components of the tissues. An aging function is introduced which, in a simplified manner, models the observed loss of stiffness of these aging elastic components due to the history of straining as well as other physiological and environmental influences. Calculations have been performed for 30 years of exposure to gravitational forces. The deformations and stress distributions in the layers of the soft tissues are described. Overall, the feasibility of using constitutive relations which reflect the highly nonlinear elastic-viscoplastic behavior of facial soft tissues in finite element based three-dimensional mechanical analyses of the human face is demonstrated. Further developments of the program are discussed in relation to possible clinical applications. Although the proposed aging function produces physically reasonable long-term response, experimental data are not yet available for more quantitative validation.     

A constituent-based model for the nonlinear viscoelastic behavior of ligaments

This paper presents a constitutive model for predicting the nonlinear viscoelastic behavior of soft biological tissues and in particular of ligaments. The constitutive law is a generalization of the well-known quasi-linear viscoelastic theory (QLV) in which the elastic response of the tissue and the time-dependent properties are independently modeled and combined into a convolution time integral. The elastic behavior, based on the definition of anisotropic strain energy function, is extended to the time-dependent regime by means of a suitably developed time discretization scheme. The time-dependent constitutive law is based on the postulate that a constituent-based relaxation behavior may be defined through two different stress relaxation functions: one for the isotropic matrix and one for the reinforcing (collagen) fibers. The constitutive parameters of the viscoelastic model have been estimated by curve fitting the stress relaxation experiments conducted on medial collateral ligaments (MCLs) taken from the literature, whereas the predictive capability of the model was assessed by simulating experimental tests different from those used for the parameter estimation. In particular, creep tests at different maximum stresses have been successfully simulated. The proposed nonlinear viscoelastic model is able to predict the time-dependent response of ligaments described in experimental works (Bonifasi-Lista et al., 2005, J. Orthopaed. Res., 23, pp. 67-76; Hingorani et al., 2004, Ann. Biomed. Eng., 32, pp. 306-312; Provenzano et al., 2001, Ann. Biomed. Eng., 29, pp. 908-214; Weiss et al., 2002, J. Biomech., 35, pp. 943-950). In particular, the nonlinear viscoelastic response which implies different relaxation rates for different applied strains, as well as different creep rates for different applied stresses and direction-dependent relaxation behavior, can be described.     

Time-dependent circumferential deformation of cortical bone upon internal radial loading

Short and long duration tests were conducted on hollow femoral bone cylinders to study the circumferential (hoop) creep response of cortical bone subjected to an intramedullary radial load. It was hypothesized that there is a stress threshold above which nonlinear creep effects dominate the mechanical response and below which the response is primarily determined by linear viscoelastic material properties. The results indicate that a hoop stress threshold exists for cortical bone, where creep strain, creep strain rate and residual strain exhibited linear behavior at low hoop stress and nonlinear behavior above the hoop stress threshold. A power-law relationship was used to describe creep strain as a function of hoop stress and time and damage morphology was assessed.     

Large deformation shear properties of liver tissue

Characterization of the mechanical properties of soft biological tissues is important for establishing the mechanical tolerances of the tissues, and for input to computational models. In this work, the viscoelastic properties of bovine liver tissue in shear loading have been measured using relaxation and constant shear rate loading. The tissue is nonlinearly viscoelastic for strains greater than 0.2%, has a yield strain of approximately 10, and shows moderate strain-rate sensitivity. The response can be modelled using a nonlinear viscoelastic differential model previously developed for brain tissue.     

Inertia effects on characterization of dynamic response of brain tissue

Modeling and simulation of traumatic brain injury (TBI) resulted from collision or blast loading requires characterization of mechanical response over a wide range of loading rates under valid testing conditions. In this study, mechanical response of fresh bovine brain tissue was studied using the two modified Kolsky bar techniques. Radial deformation behavior of annular specimens, which are typically used to characterize the dynamic uniaxial compressive response of biological tissues, was examined using a modified Kolsky bar and a high speed camera to collect images while the specimen deforms at an axial strain rate of 2000s(-1). The high-speed images revealed inhomogeneous specimen deformation possibly brought about by radial inertia and causing a multi-axial stress state. To acquire valid stress-strain results that can be used to produce constitutive behavior of the soft materials, a novel torsion technique was developed to obtain pure shear response at dynamic loading rates. Experimental results show clear differences in the material response using the two methods. These results indicate that the previously demonstrated annular specimen geometry aimed at reducing inertia induced stress components for high rate soft materials uniaxial-compressive testing may still possess a significant component of radial inertia induced radial stress which consequently caused the observed inhomogeneous deformation in brain tissue test samples.     

Mechanical characterization of anisotropic planar biological soft tissues using large indentation: a computational feasibility study

Traditionally, the complex mechanical behavior of planar soft biological tissues is characterized by (multi)axial tensile testing. While uniaxial tests do not provide sufficient information for a full characterization of the material anisotropy, biaxial tensile tests are difficult to perform and tethering effects limit the analyses to a small central portion of the test sample. In both cases, determination of local mechanical properties is not trivial. Local mechanical characterization may be performed by indentation testing. Conventional indentation tests, however, often assume linear elastic and isotropic material properties, and therefore these tests are of limited use in characterizing the nonlinear, anisotropic material behavior typical for planar soft biological tissues. In this study, a spherical indentation experiment assuming large deformations is proposed. A finite element model of the aortic valve leaflet demonstrates that combining force and deformation gradient data, one single indentation test provides sufficient information to characterize the local material behavior. Parameter estimation is used to fit the computational model to simulated experimental data. The aortic valve leaflet is chosen as a typical example. However, the proposed method is expected to apply for the mechanical characterization of planar soft biological materials in general.     

Modeling of time-dependent force response of fingertip to dynamic loading

An extended exposure to repeated loading on fingertip has been associated to many vascular, sensorineural, and musculoskeletal disorders in the fingers, such as carpal tunnel syndrome, hand-arm vibration syndrome, and flexor tenosynovitis. A better understanding of the pathomechanics of these sensorineural and vascular diseases in fingers requires a formulation of a biomechanical model of the fingertips and analyses to predict the mechanical responses of the soft tissues to dynamic loading. In the present study, a model based on finite element techniques has been developed to simulate the mechanical responses of the fingertips to dynamic loading. The proposed model is two-dimensional and incorporates the essential anatomical structures of a finger: skin, subcutaneous tissue, bone, and nail. The skin tissue is assumed to be hyperelastic and viscoelastic. The subcutaneous tissue was considered to be a nonlinear, biphasic material composed of a hyperelastic solid and an invicid fluid, while its hydraulic permeability was considered to be deformation dependent. Two series of numerical tests were performed using the proposed finger tip model to: (a) simulate the responses of the fingertip to repeated loading, where the contact plate was assumed to be fixed, and the bone within the fingertip was subjected to a prescribed sinusoidal displacement in vertical direction; (b) simulate the force response of the fingertip in a single keystroke, where the keyboard was composed of a hard plastic keycap, a rigid support block, and a nonlinear spring. The time-dependent behavior of the fingertip under dynamic loading was derived. The model predictions of the time-histories of force response of the fingertip and the phenomenon of fingertip separation from the contacting plate during cyclic loading agree well with the reported experimental observations.     

Nonlinear elastic-viscoplastic constitutive equations for aging facial tissues

This paper reports on the initial stages of a project to simulate the nonlinear mechanical behavior of an aging human face. A cross-section of the facial structure is considered to consist of a multilayered composite of tissues with differing mechanical behavior. The constitutive properties of these tissues are incorporated into a finite element model of the three-dimensional facial geometry. Relatively short time (elastic-viscoplastic) behavior is governed by equations previously developed which are consistent with mechanical tests. The long time response is controlled by the aging elastic components of the tissues. An aging function is introduced which, in a simplified manner, captures the observed loss of stiffness of these aging elastic components due to the history of straining as well as other physiological and environmental influences. Calculations have been performed for 30 years of exposure to gravitational forces. Progressive gravimetric soft tissue descent is simulated, which is regarded as the main indication of facial aging. Results are presented for the deformations and stress distributions in the layers of the soft tissues.     

Identifying a minimal rheological configuration: a tool for effective and efficient constitutive modeling of soft tissues

We describe a modeling methodology intended as a preliminary step in the identification of appropriate constitutive frameworks for the time-dependent response of biological tissues. The modeling approach comprises a customizable rheological network of viscous and elastic elements governed by user-defined 1D constitutive relationships. The model parameters are identified by iterative nonlinear optimization, minimizing the error between experimental and model-predicted structural (load-displacement) tissue response under a specific mode of deformation. We demonstrate the use of this methodology by determining the minimal rheological arrangement, constitutive relationships, and model parameters for the structural response of various soft tissues, including ex vivo perfused porcine liver in indentation, ex vivo porcine brain cortical tissue in indentation, and ex vivo human cervical tissue in unconfined compression. Our results indicate that the identified rheological configurations provide good agreement with experimental data, including multiple constant strain rate load/unload tests and stress relaxation tests. Our experience suggests that the described modeling framework is an efficient tool for exploring a wide array of constitutive relationships and rheological arrangements, which can subsequently serve as a basis for 3D constitutive model development and finite-element implementations. The proposed approach can also be employed as a self-contained tool to obtain simplified 1D phenomenological models of the structural response of biological tissue to single-axis manipulations for applications in haptic technologies.     

Magnetic resonance elastography in nonlinear viscoelastic materials under load

Characterisation of soft tissue mechanical properties is a topic of increasing interest in translational and clinical research. Magnetic resonance elastography (MRE) has been used in this context to assess the mechanical properties of tissues in vivo noninvasively. Typically, these analyses rely on linear viscoelastic wave equations to assess material properties from measured wave dynamics. However, deformations that occur in some tissues (e.g. liver during respiration, heart during the cardiac cycle, or external compression during a breast exam) can yield loading bias, complicating the interpretation of tissue stiffness from MRE measurements. In this paper, it is shown how combined knowledge of a material’s rheology and loading state can be used to eliminate loading bias and enable interpretation of intrinsic (unloaded) stiffness properties. Equations are derived utilising perturbation theory and Cauchy’s equations of motion to demonstrate the impact of loading state on periodic steady-state wave behaviour in nonlinear viscoelastic materials. These equations demonstrate how loading bias yields apparent material stiffening, softening and anisotropy. MRE sensitivity to deformation is demonstrated in an experimental phantom, showing a loading bias of up to twofold. From an unbiased stiffness of \(4910.4 \pm 635.8\) Pa in unloaded state, the biased stiffness increases to 9767.5 \(\pm \,\)1949.9 Pa under a load of \(\approx \) 34% uniaxial compression. Integrating knowledge of phantom loading and rheology into a novel MRE reconstruction, it is shown that it is possible to characterise intrinsic material characteristics, eliminating the loading bias from MRE data. The framework introduced and demonstrated in phantoms illustrates a pathway that can be translated and applied to MRE in complex deforming tissues. This would contribute to a better assessment of material properties in soft tissues employing elastography.

     

A recruitment-based rheological model for mechanical behavior of soft tissues

When lung tissue is subjected to finite deformations, phenomena appear that can only be described using nonlinear models. This paper considers the lung as a material composed of two elements, a continuous phase that acts uninterruptedly and a second phase composed of fiber elements that are recruited progressively into the mechanical process. Each individual fiber participates in the mechanical response of the set only when the deformation is above a certain value. A nine-parameter model was designed adopting standard viscoelastic elements both for the matrix and for each of the fibers. The mechanical behavior of the lung can be reproduced by a fitting process with standard numerical procedures in both dynamic-mechanical measurements and stress relaxation processes. Mechanical stress relaxation tests and dynamic-mechanical measurements have been carried out on subpleural parenchymal strips from rat lung. The model permits the reproduction of lung behavior in both types of measurements. The results show a recruitment ratio that decreases with deformation and the nonparticipation of the parallel matrix fraction in the lung's mechanical response so that a uniaxial transmission of force in the lung occurs via the recruited elements and the matrix series.     

A model of fracture testing of soft viscoelastic tissues

Fracture, or tear, toughness of soft tissues can be computed from the work of fracture divided by the area of new crack surface. For soft tissues without significant plastic deformation, total work, which can be measured experimentally, is composed of the sum of fracture and viscoelastic work. In order to deduce fracture work, a method is needed to estimate viscoelastic work.Two different methods (Ph.D. Dissertation, University of Minnesota, 2000; J. Mater. Sci.: Mater. Med. 12 (2001) 327) have been proposed to estimate viscoelastic work in a fracture test of a soft tissue. The relative merits of these methods are unknown because the true viscoelastic work in an experiment is unknown. In order to characterize the accuracy of these methods, a theoretical model of crack propagation of viscoelastic soft tissue in a tensile test is presented, from which the exact viscoelastic work is calculated. The material is assumed to obey the standard linear solid model.The "exact" solution for the viscoelastic work during the fracture is computed from the model and compared with the work estimated by the two methods. It was found that both methods tend to underestimate the viscoelastic work done, and thus overestimate the fracture work and fracture toughness, although the errors were greater with the Fedewa method. It was further found that low displacement rates can give rise to a "snap" effect, where rapid crack growth can cause a disproportionate amount of viscoelastic energy to be dissipated during unloading. This modeling approach may be useful in evaluating other experimental methods of soft tissue fracture.     

Characterization of a hyper-viscoelastic phantom mimicking biological soft tissue using an abdominal pneumatic driver with magnetic resonance elastography (MRE)

The purpose of this study was to create a polymer phantom mimicking the mechanical properties of soft tissues using experimental tests and rheological models. Multifrequency Magnetic Resonance Elastography (MMRE) tests were performed on the present phantom with a pneumatic driver to characterize the viscoelastic (μ, η) properties using Voigt, Maxwell, Zener and Springpot models. To optimize the MMRE protocol, the driver behavior was analyzed with a vibrometer. Moreover, the hyperelastic properties of the phantom were determined using compressive tests and Mooney-Rivlin model. The range of frequency to be used with the round driver was found between 60 Hz and 100 Hz as it exhibits one type of vibration mode for the membrane. MRE analysis revealed an increase in the shear modulus with frequency reflecting the viscoelastic properties of the phantom showing similar characteristic of soft tissues. Rheological results demonstrated that Springpot model better revealed the viscoelastic properties (μ=3.45 kPa, η=6.17 Pas) of the phantom and the Mooney-Rivlin coefficients were C(10)=1.09.10(-2) MPa and C(01)=-8.96.10(-3) MPa corresponding to μ=3.95 kPa. These studies suggest that the phantom, mimicking soft tissue, could be used for preliminary MRE tests to identify the optimal parameters necessary for in vivo investigations. Further developments of the phantom may allow clinicians to more accurately mimic healthy and pathological soft tissues using MRE.     

Singular perturbation analysis of the nonlinear, flow-dependent compressive stress relaxation behavior of articular cartilage

The dominant mechanism giving rise to the viscoelastic response of articular cartilage during compression is the nonlinear diffusive interaction of the fluid and solid phases of the tissue as they flow relative to one another. The present study is concerned with the role of this interaction under uniaxial stress relaxation in compression. The model is a biphasic mixture of fluid and solid which incorporates the strain-dependent permeability found earlier from permeation experiments. When a ramp-displacement is imposed on the articular surface, simple, but accurate, asymptotic approximations are derived for the deformation and stress fields in the tissue for slow and moderately fast rates of compression. They are shown to agree very well with experiment and they provide a simple means for determining the material parameters. Moreover, they lead to important insights into the role of the flow-dependent viscoelastic nature of articular cartilage and other hydrated biological tissues.