A strain energy function for lung parenchyma

The strain energy for the air-filled lung is calculated from a model of the parenchymal microstructure. The energy is the sum of the surface energy and the elastic energies of two tissue components. The first of these is the peripheral tissue system that provides the recoil pressure of the saline-filled lung, and the second is the system of line elements that form the free edges of the alveolar walls bordering the alveolar ducts. The computed strain energy is consistent with the observed linear elastic behavior of parenchyma and the data on large deformations around blood vessels.     

Parameter estimation and sensitivity analysis of a nonlinearly elastic static lung model

A model for the static pressure-volume behavior of the lung parenchyma based on a pseudo-elastic strain energy function was tested. Values of the model parameters and their variances were estimated by an optimal least-squares fit of the model-predicted pressures to the corresponding data from excised, saline-filled dog lungs. Although the model fit data from twelve lungs very well, the coefficients of variation for parameter values differed greatly. To analyze the sensitivity of the model output to its parameters, we examined an approximate Hessian, H, of the least-squares objective function. Based on the determinant and condition number of H, we were able to set formal criteria for choosing the most reliable estimates of parameter values and their variances. This in turn allowed us to specify a normal range of parameter values for these dog lungs. Thus the model not only describes static pressure-volume data, but also uses the data to estimate parameters from a fundamental constitutive equation. The optimal parameter estimation and sensitivity analysis developed here can be widely applied to other physiologic systems.     

The elastic behavior of the urinary bladder for large deformations

The purpose of this paper is to investigate the theoretical basis for the pressure-distension behavior of the urinary bladder. A finite strain theory is developed for hollow spherical structures and it is shown that the Treloar model is a good prototype only for rubber balloons. The pressure-extension ratio relationship is inverted to lead a general form of strain energy function, and fitted by an empirical relation involving one exponential. The following form of strain energy function is derived: W(lambda, lambda, lambda -2) = C1 (P(1), a) + P(1)C2 (a, lambda)ea(lambda -1). Where C1(P(1), a) is a constant (N m-2), P(1) is the initial pressure, a is the rate of pressure increase and C2 (a, lambda) a third degree polynomial relation. P(1) and a are experimentally determined through volumetric pressure-distension data. It is verified that this type of energy function is also valid for uniaxial loading experiments by testing strips coming from the same bladder for which P(1) and a were computed. There is a good agreement between the experimental points and the theoretical stress-strain relation. Finally, the strain energy function is plotted as a function of the first strain invariant and appears to be of an exponential nature.     

Strain energy density as a rupture criterion for the kidney: impact tests on porcine organs, finite element simulation, and a baseline comparison between human and porcine tissues

High-velocity (up to 25 m/s) impact tests were performed on pig kidneys to characterize failure behavior at deformation rates associated with traumatic injury. Cylindrical tissue samples (n = 45) and whole perfused organs (n = 34) were impacted using both falling weights and a high-velocity pneumatic projectile impactor. Impact energy was incrementally increased until visible rupture occurred. The strain energy density failure threshold fell between 25 and 60 kJ/m3 for excised porcine tissue samples, and between 15 and 30 kJ/m3 for whole, perfused organs. The relationship between localized failure in whole organ impacts and tissue level failure thresholds observed in cylindrical tissue samples was explored using a detailed finite element model of the human kidney. The model showed good correlation between experimentally observed injury patterns and predicted strain energy density distributions within the renal parenchyma. Finally, to facilitate interpretation of the porcine renal impact results with regard to human trauma, quasi-static compression test results of freshly excised human kidney cortex samples (n = 30) were compared against similar tests on pig kidneys. Human tissues failed at Lagrange strain levels similar to porcine tissue (63+/-6.3%), but at 52% lower Lagrange stress (116+/-28 kPa), and 35% lower strain energy density (17.1+/-4.4 kJ/m3). Thus conservative interpretation of porcine test results is recommended.     

Modeling of weak blast wave propagation in the lung

Blast injuries of the lung are the most life-threatening after an explosion. The choice of physical parameters responsible for trauma is important to understand its mechanism. We developed a one-dimensional linear model of an elastic wave propagation in foam-like pulmonary parenchyma to identify the possible cause of edema due to the impact load. The model demonstrates different injury localizations for free and rigid boundary conditions. The following parameters were considered: strain, velocity, pressure in the medium and stresses in structural elements, energy dissipation, parameter of viscous criterion. Maximum underpressure is the most suitable wave parameter to be the criterion for edema formation in a rabbit lung. We supposed that observed scattering of experimental data on edema severity is induced by the physiological variety of rabbit lungs. The criterion and the model explain this scattering. The model outlines the demands for experimental data to make an unambiguous choice of physical parameters responsible for lung trauma due to impact load.     

A transversely isotropic constitutive model of excised guinea pig spinal cord white matter

Narrowing of the spinal canal generates an amalgamation of stresses within the spinal cord parenchyma. The tissue’s stress state cannot be quantified experimentally; it must be described using computational methods, such as finite element analysis. The objective of this research was to propose a compressible, transversely isotropic constitutive model, an augmentation of the isotropic Mooney–Rivlin hyperelastic strain energy function, to describe the guinea pig spinal cord white matter. Model parameters were derived from a combination of inverse finite element analysis on transverse compression experiments and least squared error analysis applied to quasi-static longitudinal tensile tests. A comparison of the residual errors between the predicted response and the experimental measurements indicated that the transversely isotropic constitutive law that incorporates an offset stretch reduced the error by a factor of four when compared to other commonly used models.     

The strain energy function in axial plant growth

A model for axial plant growth is formulated based on conservation of energy. The model derivation assumes that a strain energy function exists to describe the dissipation of potential energy associated with water uptake, mechanical deformation, and biosynthesis during growth. The derivation does not, however, make any further assumption on the mathematical form of this constitutive relation. The model is employed to investigate possible forms of the strain energy function as applied to steady root growth. Solutions of the nonlinear partial differential equations governing growth are given for cases when the third derivative of the strain energy function is >, <, or =0. These three cases encompass a multitude of mathematical forms of the strain energy function. The resulting solutions are compared with the realization of steady axial root growth. The results of this analysis indicate that a quadratic form of the strain energy function best described steady growth. This conclusion is consistent with previous assumptions on the form of constitutive relations for growth, and allows further interpretation on the water relations, mechanical, and biosynthetic energies associated with plant growth.Research support provided by state and federal funds appropriated to the OSU/OARDC. Journal article no. 12–88     

Nonlinear mechanical property of tracheal cartilage: a theoretical and experimental study

BACKGROUND: Despite being the stiffest airway of the bronchial tree, the trachea undergoes significant deformation due to intrathoracic pressure during breathing. The mechanical properties of the trachea affect the flow in the airway and may contribute to the biological function of the lung. METHOD: A Fung-type strain energy density function was used to investigate the nonlinear mechanical behavior of tracheal cartilage. A bending test on pig tracheal cartilage was performed and a mathematical model for analyzing the deformation of tracheal cartilage was developed. The constants included in the strain energy density function were determined by fitting the experimental data. RESULT: The experimental data show that tracheal cartilage is a nonlinear material displaying higher strength in compression than in tension. When the compression forces varied from -0.02 to -0.03N and from -0.03 to -0.04N, the deformation ratios were 11.03+/-2.18% and 7.27+/-1.59%, respectively. Both were much smaller than the deformation ratios (20.01+/-4.49%) under tension forces of 0.02 to 0.01N. The Fung-type strain energy density function can capture this nonlinear behavior very well, whilst the linear stress-strain relation cannot. It underestimates the stability of trachea by exaggerating the displacement in compression. This study may improve our understanding of the nonlinear behavior of tracheal cartilage and it may be useful for the future study on tracheal collapse behavior under physiological and pathological conditions.     

A model for the nonlinear elastic response of large arteries

The nonlinear elastic response of large arteries subjected to finite deformations due to action of biaxial principal stresses, is described by simple constitutive equations. Generalized measures of strain and stress are introduced to account for material nonlinearity. This also ensures the existence of a strain energy density function. The orthotropic elastic response is described via quasi-linear relations between strains and stresses. One nonlinear parameter which defines the measures of strain and stress, and three elastic moduli are assumed to be constants. The lateral strain parameters (equivalent to Poisson's ratios in infinitesimal deformations) are deformation dependent. This dependence is defined by empirical relations developed via the incompressibility condition, and by the introduction of a fifth material parameter. The resulting constitutive model compares well with biaxial experimental data of canine carotid arteries.     

Finite deformation theory for in vivo human skin

A finite deformation mathematical model of in vivo human skin has been developed for the normal physiological load range. Uniaxial load-deformation measurements were carried out with a non-invasive extensometer and utilized in formulating the model. The in vivo strain energy function was found to be a linear function of the first two strain invariants and a quadratic function of the third strain invariant. Only three independent constants were necessary to specify the strain energy function completely for the upper extremities of human volunteers.     

A model of non-uniform lung parenchyma distortion

A finite element model of mammalian lung parenchyma is used to study the effect of large non-uniform distortions on lung elastic behaviour. The non-uniform distortion is a uni-axial stretch from an initial state of uniform pressure expansion. For small distortions, the parenchymal properties are linearly isotropic and described by two elastic moduli. However, for large distortions, the parenchyma has anisotropic non-linear elastic properties described by five independent elastic moduli dependent on the degree of distortion; they are computed for a range of distortions and initial pressures. Ez, the Young's modulus in the direction of stretch, increases significantly with distortion (epsilon(z)) while Ex, the Young's modulus in the plane perpendicular to the stretch, is approximately constant. The greater the initial pressure, the bigger the difference between the two moduli at larger distortion strains. The shear modulus G(xz) is approximately independent of degree of distortion except at the highest initial pressure. The Poisson's ratio, nu(xz) is approximately constant with distortion strain for lower initial pressures, but increases significantly with epsilon(z) at higher pressures. Model predictions of the relation between G(xz) and initial uniform inflation pressure show a good correlation with reported experimental data for small distortion strains in a range of species. The model also exhibits similar behaviour to the experimentally measured uni-axial large deformations of a tri-axially pre-loaded block of parenchyma (Hoppin et al., 1975, Journal of Applied Physiology 39, 742-751).     

Model sensitivity and use of the comparative finite element method in mammalian jaw mechanics: mandible performance in the gray wolf

Finite Element Analysis (FEA) is a powerful tool gaining use in studies of biological form and function. This method is particularly conducive to studies of extinct and fossilized organisms, as models can be assigned properties that approximate living tissues. In disciplines where model validation is difficult or impossible, the choice of model parameters and their effects on the results become increasingly important, especially in comparing outputs to infer function. To evaluate the extent to which performance measures are affected by initial model input, we tested the sensitivity of bite force, strain energy, and stress to changes in seven parameters that are required in testing craniodental function with FEA. Simulations were performed on FE models of a Gray Wolf (Canis lupus) mandible. Results showed that unilateral bite force outputs are least affected by the relative ratios of the balancing and working muscles, but only ratios above 0.5 provided balancing-working side joint reaction force relationships that are consistent with experimental data. The constraints modeled at the bite point had the greatest effect on bite force output, but the most appropriate constraint may depend on the study question. Strain energy is least affected by variation in bite point constraint, but larger variations in strain energy values are observed in models with different number of tetrahedral elements, masticatory muscle ratios and muscle subgroups present, and number of material properties. These findings indicate that performance measures are differentially affected by variation in initial model parameters. In the absence of validated input values, FE models can nevertheless provide robust comparisons if these parameters are standardized within a given study to minimize variation that arise during the model-building process. Sensitivity tests incorporated into the study design not only aid in the interpretation of simulation results, but can also provide additional insights on form and function.     

Quantifying nonlinear anisotropic elastic material properties of biological tissue by use of membrane inflation

Determination of material parameters for soft tissue frequently involves regression of material parameters for nonlinear, anisotropic constitutive models against experimental data from heterogeneous tests. Here, parameter estimation based on membrane inflation is considered. A four parameter nonlinear, anisotropic hyperelastic strain energy function was used to model the material, in which the parameters are cast in terms of key response features. The experiment was simulated using finite element (FE) analysis in order to predict the experimental measurements of pressure versus profile strain. Material parameter regression was automated using inverse FE analysis; parameter values were updated by use of both local and global techniques, and the ability of these techniques to efficiently converge to a best case was examined. This approach provides a framework in which additional experimental data, including surface strain measurements or local structural information, may be incorporated in order to quantify heterogeneous nonlinear material properties.     

A bimodular polyconvex anisotropic strain energy function for articular cartilage

A strain energy function for finite deformations is developed that has the capability to describe the nonlinear, anisotropic, and asymmetric mechanical response that is typical of articular cartilage. In particular, the bimodular feature is employed by including strain energy terms that are only mechanically active when the corresponding fiber directions are in tension. Furthermore, the strain energy function is a polyconvex function of the deformation gradient tensor so that it meets material stability criteria. A novel feature of the model is the use of bimodular and polyconvex "strong interaction terms" for the strain invariants of orthotropic materials. Several regression analyses are performed using a hypothetical experimental dataset that captures the anisotropic and asymmetric behavior of articular cartilage. The results suggest that the main advantage of a model employing the strong interaction terms is to provide the capability for modeling anisotropic and asymmetric Poisson's ratios, as well as axial stress-axial strain responses, in tension and compression for finite deformations.     

An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment

This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated by means of a rheological model consisting of a spring coupled with a Maxwell element in parallel. In this context, the isochoric free energy function is decomposed into elastic equilibrium and viscous non-equilibrium parts. The baseline elastic response is modelled by the orthotropic model of Holzapfel and Ogden [Holzapfel GA, Ogden RW. 2009. Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans Roy Soc A Math Phys Eng Sci. 367:3445–3475]. The essential aspect of the proposed model is the account of distinct relaxation mechanisms for each orientation direction. To this end, the non-equilibrium response of the free energy function is constructed in the logarithmic strain space and additively decomposed into three anisotropic parts, denoting fibre, sheet and normal directions each accompanied by a distinct dissipation potential governing the evolution of viscous strains associated with each orientation direction. The evolution equations governing the viscous flow have an energy-activated nonlinear form. The energy storage in the Maxwell branches has a quadratic form leading to a linear stress–strain response in the logarithmic strain space. On the numerical side, the algorithmic aspects suitable for the implicit finite element method are discussed in a Lagrangian setting. The model shows excellent agreement compared to experimental data obtained from the literature. Furthermore, the finite element simulations of a heart cycle carried out with the proposed model show significant deviations in the strain field relative to the elastic solution.     

A membrane model from implicit elasticity theory: application to visceral pleura

A Fungean solid is derived for membranous materials as a body defined by isotropic response functions whose mathematical structure is that of a Hookean solid where the elastic constants are replaced by functions of state derived from an implicit, thermodynamic, internal energy function. The theory utilizes Biot’s (Lond Edinb Dublin Philos Mag J Sci 27:468–489, 1939) definitions for stress and strain that, in one-dimension, are the stress/strain measures adopted by Fung (Am J Physiol 28:1532–1544, 1967) when he postulated what is now known as Fung’s law. Our Fungean membrane model is parameterized against a biaxial data set acquired from a porcine pleural membrane subjected to three, sequential, proportional, planar extensions. These data support an isotropic/deviatoric split in the stress and strain-rate hypothesized by our theory. These data also demonstrate that the material response is highly nonlinear but, otherwise, mechanically isotropic. These data are described reasonably well by our otherwise simple, four-parameter, material model.     

A 3D active-passive numerical skeletal muscle model incorporating initial tissue strains. Validation with experimental results on rat tibialis anterior muscle

This paper presents a three-dimensional finite element model of skeletal muscle and its validation incorporating inital tissue strains. A constitutive relation was determined by using a convex free strain energy function (SEF) where active and passive response contributions were obtained fitting experimental data from the rat tibialis anterior (TA) muscle. The passive and active finite strains response was modelled within the framework of continuum mechanics by a quasi-incompressible transversely isotropic material formulation. Magnetic resonance images (MRI) were obtained to reconstruct the external geometry of the TA. This geometry includes initial strains also taken into account in the numerical model. The numerical results show excellent agreement with the experimental results when comparing reaction force-extension curves both in passive and active tests. The proposed constitutive model for the muscle is implemented in a subroutine in the commercial finite element software package ABAQUS.     

Determination of material parameters of the two-dimensional Holzapfel–Weizsäcker type model based on uniaxial extension data of arterial walls

Soft tissues are anisotropic materials yet a majority of mechanical property tests have been uniaxial, which often failed to recapitulate the tensile response in other directions. This paper aims to study the feasibility of determining material parameters of anisotropic tissues by uniaxial extension with a minimal loss of anisotropic information. We assumed that by preselecting a certain constitutive model, we could give the constitutive parameters based on uniaxial extension data from orthogonal strip samples. In our study, the Holzapfel–Weizsäcker type strain energy density function (H–W model) was used to determine the material parameters of arterial walls from two fresh donation bodies. The key points we applied were the relationships between strain components in uniaxial tensile tests and the methods of stochastic optimisation. Further numerical experiments were taken. The estimate–effect ratio, defined by the number of data with the precision of estimation less than 0.5% over whole size of data, was calculated to demonstrate the feasibility of our method. The material parameters for Chinese aorta and pulmonary artery were given with the maximum root mean square (RMS) errors 0.042, and the minimal estimate–effect ratio in numerical experiments was 90.79%. Our results suggest that the constitutive parameters of arterial walls can be determined from uniaxial extension data, given the passive mechanical behaviour governed by H–W model. This method may apply to other tissues using different constitutive models.     

Reassessment of brain elasticity for analysis of biomechanisms of hydrocephalus

This paper presents results from a finite element study of the biomechanics of hydrocephalus, with special emphasis on a reassessment of the parenchyma elastic modulus. A two-dimensional finite element model of the human brain/ventricular system is developed and analysed under hydrocephalic loading conditions. It is shown that the Young's modulus of the brain parenchyma used in previous studies (3000-10000 Pa) corresponds to strain rates much higher than those present in hydrocephalic brains. Consideration of the brain's viscoelasticity leads to the derivation of a considerably lower modulus value of approximately 584 Pa.     

Structure of the developing pea seed coat and the post-phloem transport pathway of nutrients

An important function of the seed coat is to deliver nutrients to the embryo. To relate this function to anatomical characteristics, the developing seed coat of pea (Pisum sativum L.) was examined by light- and cryo-scanning electron microscopy (cryo-SEM) from the late pre-storage phase until the end of seed filling. During this time the apparently undifferentiated seed coat tissues evolve into the epidermal macrosclereids, the hypodermal hourglass cells, chlorenchyma, ground parenchyma and branched parenchyma. Using the fluorescent symplast tracer 8-hydroxypyrene-1,3,6-trisulfonic acid, it could be demonstrated that solutes imported by the phloem move into the chlorenchyma and ground parenchyma, but not into the branched parenchyma. From a comparison with literature data of common bean (Phaseolus vulgaris L.) and broad bean (Vicia faba L.), it is concluded that in the three species different parenchyma layers, but not the branched parenchyma, may be involved in the post-phloem symplasmic transport of nutrients in the seed coat. In pea, the branched parenchyma dies during the storage phase, and its cell wall remnants then form the boundary layer between the living seed coat parenchyma cells and the cotyledons. Using cryo-SEM, clear images were obtained of this boundary layer which showed that many intracellular spaces in the seed coat parenchyma are filled with an aqueous solution. This is suggested to facilitate the diffusion of nutrients from the site of unloading towards the cotyledons.