A model for aortic growth based on fluid shear and fiber stresses.

Stress-modulated growth in the aorta is studied using a theoretical model. The model is a thick-walled tube composed of two pseudoelastic, orthotropic layers representing the intima/media and the adventitia. Both layers are assumed to follow a growth law in which the time rates of change of the growth stretch ratios depend linearly on the local smooth muscle fiber stress and on the shear stress due to blood flow on the endothelium. Using finite elasticity theory modified to include volumetric growth, we computed temporal changes in stress, geometry, and opening angle (residual strain) during development and following the onset of sudden hypertension. For appropriate values of the coefficients in the growth law, the model yields results in reasonable agreement with published data for global and local growth of the rat aorta.     

Measuring,reversing, and modeling the mechanical changes due to the absence of Fibulin-4 in mouse arteries

Mice with a smooth muscle cell (SMC)-specific deletion of Fibulin-4 (SMKO) show decreased expression of SMC contractile genes, decreased circumferential compliance, and develop aneurysms in the ascending aorta. Neonatal administration of drugs that inhibit the angiotensin II pathway encourages the expression of contractile genes and prevents aneurysm development, but does not increase compliance in SMKO aorta. We hypothesized that multidimensional mechanical changes in the aorta and/or other elastic arteries may contribute to aneurysm pathophysiology. We found that the SMKO ascending aorta and carotid artery showed mechanical changes in the axial direction. These changes were not reversed by angiotensin II inhibitors, hence reversing the axial changes is not required for aneurysm prevention. Mechanical changes in the circumferential direction were specific to the ascending aorta; therefore, mechanical changes in the carotid do not contribute to aortic aneurysm development. We also hypothesized that a published model of postnatal aortic growth and remodeling could be used to investigate mechanisms behind the changes in SMKO aorta and aneurysm development over time. Dimensions and mechanical behavior of adult SMKO aorta were reproduced by the model after modifying the initial component material constants and the aortic dilation with each postnatal time step. The model links biological observations to specific mechanical responses in aneurysm development and treatment.     

腹主动脉瘤生长过程的蠕变力学模型

基于增强对腹主动脉瘤生长过程的理解、为腹主动脉瘤临床手术提供参考的目的,本文根据腹主动脉瘤的生长物理机制,提出了以蠕变力学为基础模拟腹主动脉瘤生长过程的模型.建立了腹主动脉瘤简化模型,利用有限元方法进行模拟计算.结果显示蠕变模型能够有效模拟腹主动脉瘤生长过程中的形态变化.参数优化模型模拟结果符合临床统计数据所示的腹主动脉瘤生长过程.本文还讨论了腹主动脉材料力学参数对模型的影响.     

A model of the structural and functional development of the normal human fetal left ventricle based on a global growth law

The purpose of this research is to study the growth of the normal human left ventricle (LV) during the fetal period from 14 to 40 weeks of gestation. A new constitutive law for the active myocardium describing the mechanical properties of the active muscle during the whole cardiac cycle has been proposed. The LV model is a thick-walled, incompressible, hyperelastic cylinder, with families of helicoidal fibers running on cylindrical surfaces [1]. Based on the works of Lin and Taber [2] done on the embryonic chick heart, we use for the human fetal heart a growth law in which the growth rate depends on the wall stresses. The parameters of the growth law are adapted to agree with sizes and volumes inferred from two dimensional ultrasound measurements performed on 18 human fetuses.Then calculations are performed to extrapolate the cardiac performance during normal growth of the fetal LV. The results presented support the idea that a growth law in which the growth rate depends linearly on the mean wall stresses averaged through the space and during whole cardiac cycle, is adapted to the normal human fetal LV development.     

A theory for water and macromolecular transport in the pulmonary artery wall with a detailed comparison to the aorta

The pulmonary artery (PA) wall, which has much higher hydraulic conductivity and albumin void space and approximately one-sixth the normal transmural pressure of systemic arteries (e.g, aorta, carotid arteries), is rarely atherosclerotic, except under pulmonary hypertension. This study constructs a detailed, two-dimensional, wall-structure-based filtration and macromolecular transport model for the PA to investigate differences in prelesion transport processes between the disease-susceptible aorta and the relatively resistant PA. The PA and aorta models are similar in wall structure, but very different in parameter values, many of which have been measured (and therefore modified) since the original aorta model of Huang et al. (23). Both PA and aortic model simulations fit experimental data on transwall LDL concentration profiles and on the growth of isolated endothelial (horseradish peroxidase) tracer spots with circulation time very well. They reveal that lipid entering the aorta attains a much higher intima than media concentration but distributes better between these regions in the PA than aorta and that tracer in both regions contributes to observed tracer spots. Solutions show why both the overall transmural water flow and spot growth rates are similar in these vessels despite very different material transport parameters. Since early lipid accumulation occurs in the subendothelial intima and since (matrix binding) reaction kinetics depend on reactant concentrations, the lower intima lipid concentrations in the PA vs. aorta likely lead to slower accumulation of bound lipid in the PA. These findings may be relevant to understanding the different atherosusceptibilities of these vessels.     

A Model of the Structural and Functional Development of the Normal Human Fetal Left Ventricle Based on a Global Growth Law

The purpose of this research is to study the growth of the normal human left ventricle (LV) during the fetal period from 14 to 40 weeks of gestation. A new constitutive law for the active myocardium describing the mechanical properties of the active muscle during the whole cardiac cycle has been proposed. The LV model is a thick-walled, incompressible, hyperelastic cylinder, with families of helicoidal fibers running on cylindrical surfaces [1] . Based on the works of Lin and Taber [2] done on the embryonic chick heart, we use for the human fetal heart a growth law in which the growth rate depends on the wall stresses. The parameters of the growth law are adapted to agree with sizes and volumes inferred from two dimensional ultrasound measurements performed on 18 human fetuses. Then calculations are performed to extrapolate the cardiac performance during normal growth of the fetal LV. The results presented support the idea that a growth law in which the growth rate depends linearly on the mean wall stresses averaged through the space and during whole cardiac cycle, is adapted to the normal human fetal LV development.     

Ventricular mechanics in diastole: material parameter sensitivity

Models of ventricular mechanics have been developed over the last 20 years to include finite deformation theory, anisotropic and inhomogeneous material properties and an accurate representation of ventricular geometry. As computer performance and the computational efficiency of the models improve, clinical application of these heart mechanics models is becoming feasible. One such application is to estimate myocardial material properties by adjusting the constitutive parameters to match wall deformation from MRI or ultrasound measurements, together with a measurement (or estimate) of ventricular pressure. Pigs are now the principal large animal model for these studies and in this paper we present the development of a new three-dimensional finite element model of the heart based on measurements of the geometry and the fibre and sheet orientations of pig hearts. The end-diastolic deformation of the model is computed using the "pole-zero" constitutive law which we have previously used to model the mechanics of passive myocardial tissue specimens. The sensitivities of end-diastolic fibre-sheet material strains and heart shape to changes in the material parameters are computed for the parameters of the pole-zero law in order to assess the utility of the models for inverse material property determination.     

Numerical study of how creep and progressive stiffening affect the growth stress formation in trees

It is not fully understood how much growth stresses affect the final quality of solid timber products in terms of, e.g. shape stability. It is, for example, difficult to predict the internal growth stress field within the tree stem. Growth stresses are progressively generated during the tree growth and they are highly influenced by climate, biologic and material-related factors. To increase the knowledge of the stress formation, a finite element model was created to study how the growth stresses develop during the tree growth. The model is an axisymmetric general plane strain model where material for all new annual rings is progressively added to the tree during the analysis. The material model used is based on the theory of small strains (where strains refer to the undeformed configuration which is good approximation for strains less than 4%) where so-called biological maturation strains (growth-related strains that form in the wood fibres during their maturation) are used as a driver for the stress generation. It is formulated as an incremental material model that takes into account elastic strain, maturation strain, viscoelastic strain and progressive stiffening of the wood material. The results clearly show how the growth stresses are progressively generated during the tree growth. The inner core becomes more and more compressed, whereas the outer sapwood is subjected to slightly increased tension. The parametric study shows that the growth stresses are highly influenced by the creep behaviour and evolution of parameters such as modulus of elasticity, micro-fibril angle and maturation strain.     

Axial nonuniformity of geometric and mechanical properties of mouse aorta is increased during postnatal growth

The hemodynamic conditions of aorta are relatively uniform prenatally and become more heterogeneous postnatally. Our objective was to quantify the heterogeneity of geometry and mechanical properties during growth and development. To accomplish this objective, we obtained a systematic set of data on the geometry and mechanical properties along the length of mouse aorta during postnatal development. C57BL/6 mice of ages 1-33 days were studied. The ascending aorta was cannulated in situ and preconditioned with several cyclic changes in pressure. We investigated the axial variations of geometry (diameter and length) and mechanical properties (stress-stain relation, elastic modulus and compliance) of the mouse aorta from the aortic valve to the common iliac. Our results show that the arterial blood pressure of mice increased from approximately 30 to 80 mmHg during the first 2 wk of life. The stretch ratio, diameter, wall (intima-media) thickness, and total lumen volume of mouse aorta increased with age. The aorta was transformed from a cylindrical tube at birth to a tapered structure during growth. Furthermore, we found the mechanical properties were fairly uniform along the length of the aorta at birth and become more nonuniform with age. We conclude that the rapid change of blood pressure and blood flow after birth alter the geometric and mechanical properties differentially along the length of the aorta. Hence, the axial nonuniformity of the aorta increases as the organ becomes more specialized during growth and development.     

A model of growth and rupture of abdominal aortic aneurysm

We present here a coupled mathematical model of growth and failure of the abdominal aortic aneurysm (AAA). The failure portion of the model is based on the constitutive theory of softening hyperelasticity where the classical hyperelastic law is enhanced with a new constant indicating the maximum energy that an infinitesimal material volume can accumulate without failure. The new constant controls material failure and it can be interpreted as the average energy of molecular bonds from the microstructural standpoint. The constitutive model is compared to the data from uniaxial tension tests providing an excellent fit to the experiment. The AAA failure model is coupled with a phenomenological theory of soft tissue growth. The unified theory includes both momentum and mass balance laws coupled with the help of the constitutive equations. The microstructural alterations in the production of elastin and remodeling of collagen are reflected in the changing macroscopic parameters characterizing tissue stiffness, strength and density. The coupled theory is used to simulate growth and rupture of an idealized spherical AAA. The results of the simulation showing possible AAA ruptures in growth are reasonable qualitatively while the quantitative calibration of the model will require further clinical observations and in vitro tests. The presented model is the first where growth and rupture are coupled.     

Analysis of protein distribution in budding yeast

Flow cytometry is a fast and sensitive method that allows monitoring of different cellular parameters on large samples of a population. Protein distributons give relevant information on growth dynamics, since they are related to the age distribution and depend on the law of growth of the population and the law of protein accumulation during the cell cycle. We analyzed protein distributions to evaluate alternative growth models for the budding yeast Saccharomyces cerevisiae and to monitor the changes in population dynamics that result from environmental modifications; such an analysis could potentially give parameters useful in the control of biotechnological processes. Theoretical protein distributions (taking into account the unequal division of yeast cells and the exponential law of protein accumulation during a cell cycle) quantitatively fit experimental distributions, once appropriate variability sources are introduced. Best fits are obtained when the protein threshold required for bud emergence increases at each new generation of parent cells.     

A finite element model of stress-mediated vascular adaptation: application to abdominal aortic aneurysms

Despite rapid expansion of our knowledge of vascular adaptation, developing patient-specific models of diseased arteries is still an open problem. In this study, we extend existing finite element models of stress-mediated growth and remodelling of arteries to incorporate a medical image-based geometry of a healthy aorta and, then, simulate abdominal aortic aneurysm. Degradation of elastin initiates a local dilatation of the aorta while stress-mediated turnover of collagen and smooth muscle compensates the loss of elastin. Stress distributions and expansion rates during the aneurysm growth are studied for multiple spatial distribution functions of elastin degradation and kinetic parameters. Temporal variations of the degradation function are also investigated with either direct time-dependent degradation or stretch-induced degradation as possible biochemical and biomechanical mechanisms for elastin degradation. The results show that this computational model has the capability to capture the complexities of aneurysm progression due to variations of geometry, extent of damage and stress-mediated turnover as a step towards patient-specific modelling.     

A finite element model of stress-mediated vascular adaptation: application to abdominal aortic aneurysms

Despite rapid expansion of our knowledge of vascular adaptation, developing patient-specific models of diseased arteries is still an open problem. In this study, we extend existing finite element models of stress-mediated growth and remodelling of arteries to incorporate a medical image-based geometry of a healthy aorta and, then, simulate abdominal aortic aneurysm. Degradation of elastin initiates a local dilatation of the aorta while stress-mediated turnover of collagen and smooth muscle compensates the loss of elastin. Stress distributions and expansion rates during the aneurysm growth are studied for multiple spatial distribution functions of elastin degradation and kinetic parameters. Temporal variations of the degradation function are also investigated with either direct time-dependent degradation or stretch-induced degradation as possible biochemical and biomechanical mechanisms for elastin degradation. The results show that this computational model has the capability to capture the complexities of aneurysm progression due to variations of geometry, extent of damage and stress-mediated turnover as a step towards patient-specific modelling.     

A constrained mixture model for developing mouse aorta

Mechanical stresses influence the structure and function of adult and developing blood vessels. When these stresses are perturbed, the vessel wall remodels to return the stresses to homeostatic levels. Constrained mixture models have been used to predict remodeling of adult vessels in response to step changes in blood pressure, axial length and blood flow, but have not yet been applied to developing vessels. Models of developing blood vessels are complicated by continuous and simultaneous changes in the mechanical forces. Understanding developmental growth and remodeling is important for treating human diseases and designing tissue-engineered blood vessels. This study presents a constrained mixture model for postnatal development of mouse aorta with multiple step increases in pressure, length and flow. The baseline model assumes that smooth muscle cells (SMCs) in the vessel wall immediately constrict or dilate the inner radius after a perturbation to maintain the shear stress and then remodel the wall thickness to maintain the circumferential stress. The elastin, collagen and SMCs have homeostatic stretch ratios and passive material constants that do not change with developmental age. The baseline model does not predict previously published experimental data. To approximate the experimental data, it must be assumed that the SMCs dilate a constant amount, regardless of the step change in mechanical forces. It must also be assumed that the homeostatic stretch ratios and passive material constants change with age. With these alterations, the model approximates experimental data on the mechanical properties and dimensions of aorta from 3- to 30-day-old mice.     

Optimizing cardiac material parameters with a genetic algorithm

Determining the unknown material parameters of intact ventricular myocardium can be challenging due to highly nonlinear material behavior. Previous studies combining a gradient-search optimization procedure with finite element analysis (FEA) were limited to two-dimensional (2D) models or simplified three-dimensional (3D) geometries. Here we present a novel scheme to estimate unknown material parameters for ventricular myocardium by combining a genetic algorithm (GA) with nonlinear finite element analysis. This approach systematically explores the domain of the material parameters. The objective function to minimize was the error between simulated strain data and finite element model strains. The proposed scheme was validated for a 2D problem using a realistic material law for ventricular myocardium. Optimized material parameters were generally within 0.5% of the true values. To demonstrate the robustness of the new scheme, unknown material parameters were also determined for a realistic 3D heart model with an exponential hyperelastic material law. When using strains from two material points, the algorithm converged to parameters within 5% of the true values. We conclude that the proposed scheme is robust when estimating myocardial material parameters in 2D and 3D models.     

Towards in vivo aorta material identification and stress estimation

This paper addresses the problem of constructing a mechanical model for the abdominal aorta and calibrating its parameters to in vivo measurable data. The aorta is modeled as a pseudoelastic, thick-walled, orthotropic, residually stressed cylindrical tube, subjected to an internal pressure. The model parameters are determined by stating a minimization problem for the model pressure and computing the optimal solution by a minimization algorithm. The data used in this study is in vivo pressure–diameter data for the abdominal aorta of a 24-year-old man. The results show that the axial, circumferential and radial stresses have magnitudes in the span 0 to 180 kPa. Furthermore, the results show that it is possible to determine model parameters directly from in vivo measurable data. In particular, the parameters describing the residual stress distribution can be obtained without interventional procedures.     

Evolution of mechanoregulation of bone growth will lead to non-optimal bone phenotypes

Mechanical forces acting on the bones during growth affect their final shape and strength. Mechanoregulation of bone growth may be recognized in embryogenesis, and also in the adaptation of the adult skeleton to changes in mechanical loading. Mechanoregulatory responses for tissues have arisen during evolution, but does evolution give rise to responses that produce optimal skeletal phenotypes? In this paper, we investigate the emergence of an optimal mechanoregulation response in a population. By combining equations describing long bone growth with a genetic algorithm to describe evolutionary change, we created a computational model to simulate the evolution of mechanoregulation in bone growth. A population of individuals is created where each individual is assigned a diploid gene set which controls the growth and remodelling of the bone. At maturity, each bone is assessed and its 'fitness' calculated; fitness is quantified as bone strength relative to bone mass. The simulation continues for many generations, and includes mutations and a varying environment. The genes present in the population are tracked and the evolution of parameters governing mechanoregulation is calculated. The results indicate that a population may converge to one bone growth algorithm but, more usually, a range of mechanoregulation algorithms for different individuals will persist after many generations. Even if the population converges to one mechanoregulation law, convergence to the 'optimum' bone was never found. Although many researchers propose that natural selection has pushed skeletal structure towards an optimum, our computational model suggests that this is unlikely to be the case.     

The orthotropic viscoelastic behavior of aortic elastin

In this paper, we studied the viscoelastic behaviors of isolated aortic elastin using combined modeling and experimental approaches. Biaxial stress relaxation and creep experiments were performed to study the time-dependent behavior of elastin. Experimental results reveal that stress relaxation preconditioning is necessary in order to obtain repeatable stress relaxation responses. Elastin exhibits less stress relaxation than intact or decellularized aorta. The rate of stress relaxation of intact and decellularized aorta is linearly dependent on the initial stress levels. The rate of stress relaxation for elastin increases linearly at stress levels below about 60 kPa; however, the rate changes very slightly at higher initial stress levels. Experimental results also show that creep response is negligible for elastin, and the intact or decellularized aorta. A quasi-linear viscoelasticity model was incorporated into a statistical mechanics based eight-chain microstructural model at the fiber level to simulate the orthotropic viscoelastic behavior of elastin. A user material subroutine was developed for finite element analysis. Results demonstrate that this model is suitable to capture both the orthotropic hyperelasticity and viscoelasticity of elastin.     

A calcium-driven mechanochemical model for prediction of force generation in smooth muscle

A new model for the mechanochemical response of smooth muscle is presented. The focus is on the res- ponse of the actin–myosin complex and on the related generation of force (or stress). The chemical (kinetic) model describes the cross-bridge interactions with the thin filament in which the calcium-dependent myosin phosphorylation is the only regulatory mechanism. The new mechanical model is based on Hill’s three-component model and it includes one internal state variable that describes the contraction/relaxation of the contractile units. It is characterized by a strain-energy function and an evolution law incorporating only a few material parameters with clear physical meaning. The proposed model satisfies the second law of thermodynamics. The results of the combined coupled model are broadly consistent with isometric and isotonic experiments on smooth muscle tissue. The simulations suggest that the matrix in which the actin–myosin complex is embedded does have a viscous property. It is straightforward for implementation into a finite element program in order to solve more complex boundary-value problems such as the control of short-term changes in lumen diameter of arteries due to mechanochemical signals.