The Effect of Strain Hardening on the Dynamic Response of Human Artery Segments

Background: When subjected to time-dependent blood pressure, human arteries undergo large deformations, exhibiting mainly nonlinear hyperelastic type of response. The mechanical response of arteries depends on the health of tissues that comprise the artery walls. Typically, healthy arteries exhibit convex strain hardening under tensile loads, atherosclerotic parts exhibit stiffer response, and aneurysmatic parts exhibit softening response. In reality, arterial dynamics is the dynamics of a propagating pulse, originating in heart ventricle, propagating along aorta, bifurcating, etc. Artery as a whole cannot be simulated as a lump ring, however its cross section can be simulated as a vibrating ring having a phase lag with respect to the other sections, creating a running pressure wave. A full mathematical model would require fluid-solid interaction modeling continuity of blood flow in a compliant vessel and a momentum equation. On the other hand, laboratory testing often uses small-length arteries, the response of which is covered by the present work. In this way, material properties that change along the artery length can be investigated. Objective: The effect of strain hardening on the local dynamic response of human arteries (excluding the full fluid-structure interaction) is examined through appropriate hyperelastic models related to the health condition of the blood vessel. Furthermore, this work aims at constituting a basis for further investigation of the dynamic response of arteries accounting for viscosity. Method: The governing equation of motion is formulated for three different hyperelastic material behaviors, based on the constitutive law proposed by Skalak et al., Hariton, and Mooney-Rivlin, associated with the hardening behavior of healthy, atherosclerotic, and aneurysmatic arteries, respectively. The differences between these modelling implementations are caused by physiology, since aneurysmatic arteries are softer and often sclerotic arteries are stiffer than healthy arteries. The response is investigated by proper normalization of the involved material parameters of the arterial walls, geometry of the arteries, load histories, time effects, and pre-stressing. The effect of each problem parameter on the arterial response has been studied. The peak response of the artery segment is calculated in terms of radial displacements, principal elongations, principal stresses, and strain-energy density. The validity of the proposed analytical models is demonstrated through comparison with previous studies that investigate the dynamic response of arterial models. Results: Important metrics that can be useful to vascular surgery are the radial deformation and the maximum strain-energy density along with the radial resonance frequencies. These metrics are found to be influenced heavily by the nonlinear strain-hardening characteristics of the model and the longitudinal pre-stressing. Conclusion: The proposed formulation permits a systematic and generalizable investigation, which, together with the low computational cost of analysis, makes it a valuable tool for calculating the response of healthy, atherosclerotic, and aneurysmatic arteries. The radial resonance frequencies can explain certain murmures developed in stenotic arteries.