1 Introduction

The mechanical properties of the aorta are important to its physiological function. The concept of hemodynamic homeostasis enables large arteries such as the aorta to transform central pulsatile pressure and flow into continuous pressure and flow in the peripheral arterioles. In this transformation, central artery stiffness has been identified as a major independent risk factor for cardiovascular disease morbidity and overall mortality [1].

Arterial stiffness is attributed to, e.g., extracellular matrix (ECM) components, mainly elastin and collagen, vascular smooth muscle cell (VSMC) tone, VSMC stiffness and cell–ECM interactions [2]. The cell–ECM interaction affects and regulates arterial mechanical function and structural integrity [3]. Shear as well as circumferential and longitudinal stress is key mechanical determinants of arterial wall remodeling. Aortic mechanical properties are considered to evolve from an interdependency of circumferential–longitudinal coupling of stress and stretch where stiffness may play a part [4, 5].

Stress, stiffness and stretch can be measured in vitro, but especially stress and stiffness are difficult to measure in vivo. A mechanical model, usually computerized, may be used to simulate stress, stiffness and stretch in vivo. Various mechanical models describing cardiovascular growth and remodeling have been proposed [5, 6]. We have developed a model using in vivo data with nonlinear deformation and material behavior observed for arteries to compute stress and stiffness [7, 8]. The model has been validated [9] and produced findings relevant to age- and sex-dependent changes in vessel wall constituents [10].

In biological tissue, the stress–strain curve is nonlinear, implicating nonlinear stiffness. Our mechanical model enables computation of wall stress as well as direct and cross-coupled stiffnesses from the nonlinear stress–strain curve. A commonly used variable such as Young’s incremental elastic modulus calculated from stress–strain curves assumes a linear material property and is a measure of direct stiffness. In a clinical environment, aortic pulse wave velocity (PWV) is regarded as a reference parameter for central aortic stiffness [1] and mainly reflects circumferential direct stiffness. Nonlinear behavior as well as cross-coupled stiffness is much less investigated. Furthermore, cross-coupled stiffness links stretch in one direction to stress in another direction and might reveal information related to directional interdependency of stress and stretch. Therefore, the aim of this study is to investigate direct and cross-coupled stiffness and its possible implications on aortic remodeling.

2 Materials and Methods

2.1 Study Subjects

Thirty healthy, non-smoking male (n = 15) and female volunteers without any medications were included in the study and divided in three age-groups: young (23–30 years, n = 10), middle aged (41–54 years, n = 10) and elderly (67–72 years, n = 10). Exclusion criteria were a history of cardiopulmonary disease, diabetes or regular medication and an ankle brachial index < 1 suggestive of vascular obstruction. Oestrogen replacement therapy was not prescribed to anyone of the women. The acquisition of pressure and diameter is summarized below and has been described elsewhere [10].

3 Non-invasive Monitoring of Diameter Changes

Non-invasive monitoring of pulsatile diameter change in the distal abdominal aorta (AA) was carried out 3–4 cm proximal to the aortic bifurcation [10]. An electronic echo-tracking instrument (Diamove, Teltec AB, Lund, Sweden) was interfaced with a real-time ultrasound scanner (EUB-240, Hitachi, Tokyo, Japan) and fitted with a 3.5 Mhz linear array transducer. The instrument had dual echo-tracking loops; thus, two separate echoes from opposite vessel walls could be tracked simultaneously. The repetition frequency was 870 Hz, temporal resolution was 1.2 ms, and the smallest detectable displacement was 7.8 µm. For static (end diastolic and systolic) aortic diameter and for pulsatile diameter change, the coefficient of variation was 5% and 16%, respectively.

3.1 Invasive Blood Pressure Measurements

Abdominal aortic (AA) blood pressure was measured invasively at the midpoint between the renal arteries and the aortic bifurcation with a 3-F (SPC 330A) or 4-F (SPC 340) micromanometer tip catheter (Millar Instruments, Houston, TX) or with a fluid-filled catheter system (pressure monitoring kit DTX + with R.O.S.E, Viggo Spectramed, Oxnard, CA) depending on the availability. The frequency response of the Millar catheter (flat range to 10 kHz) was higher than in the fluid-filled system (flat range 35 Hz [3 dB]). However, curves from one cardiac cycle from each system were superimposed on each other using a Blood Systems Calibrator (Bio Tech Model 601 A, Old Mill Street, Burlington, VT) showing similar systolic blood pressures and pulsatile amplitude when compared.

The data acquisition system allowed for simultaneous monitoring of blood pressure and vessel diameter with a maximum registration duration of 11 s. The system contained a personal computer type 386 (Express, Tokyo, Japan) and a 12-bit analog-to-digital converter (Analogue Devices, Norwood, MA) allowing a sampling frequency of 290 Hz each for both signals. Example of acquired data is found in Fig. 1.

Fig. 1
figure 1

Example of acquired data. A and B illustrate measured data from the abdominal aorta and identification results from PIMMP. A Measured blood pressure (upper left panel) and inner radius (lower left panel) for a young male. Right panel shows the pressure-radius response and the used post-processed average signal (black). B Left panel shows wall stress vs. inner radius in the circumferential direction for a young male. Measured data with resulting fitted curve from PIMMP representing total stress is illustrated. Right panel shows circumferential and longitudinal direct as well as cross-coupled wall stiffness vs. inner radius for a young male

3.2 The Identification of Model Parameters and the Mechanical Model

The identification of model parameters and the mechanical model used have been described in previous publications [8, 10]. Also see Appendix for further information.

The parameter identification method for mechanical parameters (PIMMP) consists of two parts, a signal processing routine and a parameter identification routine including a nonlinear mechanical model. In the first part, measured blood pressure and diameter were processed in MATLAB. The data consisting of approximately 8–10 cycles (heartbeats) were lowpass filtered with a fourth-order Butterworth filter with a cutoff frequency of 15 Hz for noise reduction and automatically adjusted for time delays from the measurement setup. Furthermore, pressure and radius were averaged over cycles [8]. In part two, the model parameters were identified through a nonlinear curve fitting of the model response to the measured pressure–radius loop, according to the following iterative algorithm (Fig. 2A):

  1. 1.

    The stresses in the arterial wall were computed by the Laplace’s law (Eq. 8 in Appendix) using the pressure-radius loop together with an estimation of the cross-sectional area (A) of the aortic wall from Åstrand et al. [11].

  2. 2.

    A second set of stresses was computed using nonlinear continuum mechanics (Eqs. 11 and 12) [12]. These model stresses are dependent on six model parameters (explained below) describing the material characteristics and the in situ pre-stress of the aortic wall [13,14,15].

  3. 3.

    By comparing the first set of stresses from 1 to the stresses from 2, an estimate of an error (difference) was obtained (Eq. 9) [8]. If the difference between errors from two consecutive iterations was smaller than a pre-set tolerance (typically 10–5), the iteration was terminated, and the parameters were considered identified. If the error exceeded the tolerance, the parameters values were updated using standard identification techniques and steps 2 and 3 were repeated.

Fig. 2
figure 2

Identification algorithm and model parameters. A and B show the identification algorithm and definitions of some model parameters. A The identification algorithm used to compute the material parameters in the abdominal aorta, see text for description. B An overview of how model parameters β, R0 and λz are defined. The upper cylinder represents the unloaded state with unit length; the lower cylinder represents the pressurized state

Six model parameters were identified, describing the characteristic (material parameters: c, k1, k2, β) and geometrical (geometrical parameters: R0, and \({\lambda }_{z}\)) properties of the aortic wall [10]. The parameters are:

  • c (Pa)—relates to the stiffness of the isotropic constituents in the vascular wall, mainly elastin.

  • k1 (Pa)—relates to the stiffness of the anisotropic constituents in the vascular wall, mainly collagen.

  • k2 (dimensionless)—reflects the crimpling or folding, cross-linking and entanglement of collagen.

  • β (°)—the angle between the circumferential direction and the principal (mean) fiber direction in the unloaded configuration (Fig. 2B).

  • R0 (mm)—the radius in a stress and stretch free (ex situ) unloaded configuration (Fig. 2B)

  • \({\lambda }_{z}\) (dimensionless)—longitudinal stretch between the in situ configuration and the (ex situ) unloaded configuration (Fig. 2B).

Values for identified parameters are found in Table 4 in Appendix.

A single stiffness constant, e.g., Young’s modulus, cannot accurately predict arterial wall stress [16]. It must be computed from the observed deformation using a set of (nonlinear) equations and associated material parameters [8, 17]. The model used herein is based on a standard Holzapfel–Gasser–Ogden (HGO) nonlinear material model with a neo-Hookean matrix reinforced by a two-family fiber structure [12]. The stress–stretch curve in the circumferential direction can be described by:

$${\sigma }_{\theta }={\lambda }_{\theta }\frac{\partial \psi }{\partial {\lambda }_{\theta }}={\sigma }_{\theta }^{\mathrm{iso}}+{\sigma }_{\theta }^{\mathrm{aniso}}=2C\left[{\lambda }_{\theta }^{2}-\frac{1}{{\left({\lambda }_{\theta }{\lambda }_{z}\right)}^{2}}\right]+4{k}_{1}\left(I-1\right){\mathrm{e}}^{{k}_{2}{\left(I-1\right)}^{2}}{\lambda }_{\theta }{\mathrm{cos}}^{2}\beta$$
(1)

where c, k1, k2 > 0 guarantee energy dissipation and \(I={\lambda }_{\theta }^{2}{\mathrm{cos}}^{2}\beta +{\lambda }_{z}^{2}{\mathrm{sin}}^{2}\beta\) with \({\lambda }_{\theta }=\frac{{R}_{0}}{{r}_{0}}\frac{4\pi {r}_{0}^{2}+A}{4\pi {R}_{0}^{2}+{\lambda }_{z}A}\). Computed circumferential stress with the result from the identification routine is illustrated in Fig. 1B.

For our purpose, incremental stiffness for a nonlinear material such as biological tissue can be estimated from the slope of the nonlinear stress–stretch curve. This corresponds to the partial derivative of stress with respect to stretch. Hence, there will be two direct stiffnesses and two cross-coupled stiffnesses:

$${S}_{\theta \theta }=\frac{\partial {\sigma }_{\theta }}{\partial {\lambda }_{\theta }}, {S}_{\theta z}=\frac{\partial {\sigma }_{\theta }}{\partial {\lambda }_{z}}, {S}_{z\theta }=\frac{\partial {\sigma }_{z}}{\partial {\lambda }_{\theta }},{S}_{zz}=\frac{\partial {\sigma }_{z}}{\partial {\lambda }_{z}}$$
(2)

where \({S}_{\theta \theta }\) and \({S}_{\theta z}\) are circumferential direct and cross-coupled stiffnesses while \({S}_{zz}\) and \({S}_{z\theta }\) are longitudinal direct and cross-coupled stiffnesses. Cross-coupled stiffness carries information of how stress in one direction is affected when the artery is stretched in another direction. Thus, \({\mathrm{S}}_{\mathrm{\theta z}}\) represents how circumferential stress varies when the vessel wall is stretched in the longitudinal direction and \({\mathrm{S}}_{\mathrm{z\theta }}\) represents how longitudinal stress varies when the vessel wall is stretched in the circumferential direction.

3.3 Computed Variables

In the parameter identification routine, the material parameters were computed (c, k1, k2 and β) as well as the geometry for the unloaded state (R0, λz). Circumferential stretch (\({\lambda }_{\theta }\)) was approximated with \({\lambda }_{\theta }=r/{R}_{0}\) where r was measured radius and R0 was identified through the parameter identification process [8]. Furthermore, the circumferential (\({\sigma }_{\theta }\)) and longitudinal (\({\sigma }_{z}\)) stresses were computed. Stiffness was computed as the partial derivative of stress with respect to stretch according to Eq. 2; thus, there were four different stiffnesses computed: circumferential direct (\({S}_{\theta \theta }\)) and cross-coupled (\({S}_{\theta z}\)) as well as longitudinal direct (\({S}_{zz}\)) and cross-coupled (\({S}_{z\theta }\)) stiffnesses.

Mean arterial pressure (MAP) was calculated as 1/3 × (SBP − DBP) + DBP, unit Pa. To convert Pa to mmHg, 1 mmHg = 133.32 Pa, was used.

Pressure, radius, circumferential and longitudinal stretch, stress and stiffness were calculated at SBP, DBP and MAP. Notice that the identified parameters are constant during one heartbeat and are not subject to different values at different pressures.

3.4 Linearization

A small amplitude perturbation analysis has been carried out to determine how cross-coupled stiffness manifests through stress. It was done through a linearization close to a point of interest (systolic blood pressure). The linearization follows standard procedures and is a first order Taylor expansion close to the point of interest [18]. For the present study, the linearization will be two dimensional (2D):

$$\left(\begin{array}{c}{\sigma }_{\theta }^{\mathrm{L}}\left({\lambda }_{\theta },{\lambda }_{z}\right)\\ {\sigma }_{z}^{\mathrm{L}}\left({\lambda }_{\theta },{\lambda }_{z}\right)\end{array}\right)=\left(\begin{array}{c}{\sigma }_{\theta }^{0}\\ {\sigma }_{z}^{0}\end{array}\right)+\left(\begin{array}{cc}{S}_{\theta \theta }^{0}& {S}_{\theta z}^{0}\\ {S}_{z\theta }^{0}& {S}_{zz}^{0}\end{array}\right)\left(\begin{array}{c}{\lambda }_{\theta }-{\lambda }_{\theta }^{0}\\ {\lambda }_{z}-{\lambda }_{z}^{0}\end{array}\right)$$
(3)

Here \({\sigma }_{\theta }^{\mathrm{L}}\) and \({\upsigma }_{z}^{\mathrm{L}}\) are the linearized functions while \({\sigma }_{\theta }\) and \({\sigma }_{z}\) are the nonlinear functions. \({S}_{\theta \theta }\), \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) are the first partial derivatives of \({\sigma }_{\theta }\), and \({\sigma }_{z}\) with respect to \({\lambda }_{\theta }\) and \({\lambda }_{z}\). \({\lambda }_{\theta }^{0}\) and \({\lambda }_{z}^{0}\) denote the point of interest and superscript “0” serves as a marker for a function calculated at the point of interest.

The linearized functions were used to quantify the influence of stretches on stress through direct and cross-coupled stiffness. By changing the stretch a small amount, q (0 < q ≤ 0.03) and taking the ratio between the change (\({\sigma }_{i}^{\mathrm{L},\mathrm{j}q})\) and the original (\({\sigma }_{i}^{\mathrm{L}})\) values, the effect of the change in stretch can be analyzed. Choosing a too large value for q will end up in calculations outside the validity of the linearization. The higher the ratio, the more impact will stretch have on stress. Four ratios were calculated.

$${R}_{ij}=\frac{{\sigma }_{i}^{\mathrm{L},jq}-{\sigma }_{i}^{\mathrm{L}}}{{\sigma }_{i}^{\mathrm{L}}} \quad (i, j =\theta \; \mathrm{ or } \; z)$$
(4)

with

$${\sigma }_{i}^{\mathrm{L},jq}-{\sigma }_{i}^{\mathrm{L}}={S}_{ij}^{0}\times {\lambda }_{j}\times q \quad (i, j =\theta \; \mathrm{ or } \; z)$$
(5)

Here superscript “q” denotes stress where stretch has been changed a small amount; i and j denote θ or z. The ratios express:

  • \({R}_{\theta \theta }\): circumferential (θ) stretch impact on circumferential (θ) stress

  • \({R}_{\theta z}\): longitudinal (z) stretch impact on circumferential (θ) stress

  • \({R}_{z\theta }\): circumferential (θ) stretch impact on longitudinal (z) stress

  • \({R}_{zz}\): longitudinal (z) stretch impact on longitudinal (z) stress

The ratios from Eq. (4) were used to assess the influence of stretch on stress via direct stiffness or cross-coupled stiffness.

The impact of stiffness on stress without the effect of stretch was assessed through the stiffness matrix in linearized Eq. (3). Circumferential direct stiffness (\({S}_{\theta \theta }\)) was the largest of the four stiffnesses. Normalizing the elements with \({S}_{\theta \theta }\) will produce relative values which can be compared within the matrix as well as between different points of interest:

$${nS}_{ij}^{0}=\frac{{S}_{ij}^{0}}{{S}_{\theta \theta }^{0}} \quad (i, j=\theta \; \mathrm{ or } \; z)$$
(6)

The matrix in Eq. (3) can be rewritten as:

$${S}_{\theta \theta }^{0}\times \left(\begin{array}{cc}1& {nS}_{\theta z}^{0}\\ {nS}_{z\theta }^{0}& {nS}_{zz}^{0}\end{array}\right)$$
(7)

Here “n” denotes normalized values. The normalized values express stiffness in terms of \({S}_{\theta \theta }\); thus, the normalized matrix element for circumferential direct stiffness will always be 1. Since \({S}_{\theta \theta }\) is largest, the other three stiffnesses will always be less than one.

3.5 Statistics

Arithmetic mean and standard deviation (SD) were calculated for all variables and are expressed as mean ± SD, if not otherwise stated. Variable values calculated at SBP and DBP were regarded as maximum and minimum. A two-way ANOVA test with complementing general linear models was used to compare sex and age groups as suggested by Field [19]. Bonferroni correction was used when multiple comparisons were performed. All parameters were assessed for dependency of age within each sex with a linear regression analysis with Pearson correlation coefficient (R2). P < 0.05 was considered significant in the ANOVA and general linear model as well as in the linear regression analysis. Significance testing is used for a descriptive purpose.

3.6 Software

MATLAB (The Mathwork, Natick, MA, US) version 8.4 (R2014b) was used for computation. IBM SPSS Statistics Version 27 (IBM Corporation, Somers, NY, US) was used for statistical analysis.

4 Results

4.1 Aortic Blood Pressure and Vessel Diameter

Baseline clinical data for the study population are found in Table 1, while abdominal aortic (AA) diameters and blood pressures are shown in Table 2. In males, elderly compared with young had a higher systolic blood pressure (SBP) value, larger AA diameter, but smaller ΔD (P < 0.05). Elderly compared with young males appeared to have a higher mean arterial pressure (MAP) value (P = 0.07) as well as pulse pressure (PP) value (P = 0.08). In females, elderly compared with young had larger AA diameters, but smaller ΔD than (P < 0.05). Males compared to females had larger AA diameters (P < 0.05). Changes in the pulsatile diameter of the AA (ΔD) and blood pressure did not differ by sex.

Table 1 Characteristics of the studied population
Table 2 Abdominal aortic pressure and diameter

4.2 Linear Model

For the linearization, the groups of elder and middle aged were combined to one group (n = 20) of males (n = 10) and females. This unified group was investigated, and results are reported in Fig. 3.

Fig. 3
figure 3

Linearized stiffness matrix at SBP, elderly and middle-aged group. Results from linearization at SBP. Normalized linear model stiffness matrix (A) and ratio between the change and the original values (B) for males (left) and females (right) at SBP in the elderly and middle-aged group. Panel A reveals the impact of stiffness alone on stress, while panel B include stretch and thus show the combined effect of stiffness and stretch on stress. There were no differences between sexes. \({nS}_{\theta \theta }\) had a higher value compared to \({nS}_{\theta z}\), \({nS}_{z\theta }\) and \({nS}_{zz}\), which, on the other hand, were of the same magnitude. In males, \({R}_{\theta z}\) had a lower value compared to \({R}_{\theta \theta }\), \({R}_{z\theta }\) and \({R}_{zz}\) which, in turn, were of the same magnitude. In females, \({R}_{\theta z}\) was lower compared to \({R}_{\theta \theta }\), and appeared to be lower compared to \({R}_{z\theta }\) and \({R}_{zz}\) which, in turn were of the same magnitude. Error bars ± 1 SD.

4.2.1 Comparing Normalized Linear Model Stiffness Matrix Elements (\({{\varvec{n}}{\varvec{S}}}_{{\varvec{\theta}}{\varvec{\theta}}}\), \({{\varvec{n}}{\varvec{S}}}_{{\varvec{\theta}}{\varvec{z}}}\), \({{\varvec{n}}{\varvec{S}}}_{{\varvec{z}}{\varvec{\theta}}}\), \({{\varvec{n}}{\varvec{S}}}_{{\varvec{z}}{\varvec{z}}}\)) at SBP, DBP and MAP

At SBP in males and females, \({nS}_{\theta \theta }\) had a 70–100% higher value compared with \({nS}_{\theta z}\), \({nS}_{z\theta }\) and \({nS}_{zz}\) (males: 1.00 ± 0.00 vs. 0.59 ± 0.05 vs. 0.51 ± 0.03 vs. 0.56 ± 0.10, respectively, P < 0.01, and females: 1.00 ± 0.00 vs. 0.57 ± 0.06 vs. 0.51 ± 0.04 vs. 0.57 ± 0.12, respectively, P < 0.01). In males and females, there were no differences when comparing \({nS}_{\theta z}\), \({nS}_{z\theta }\) and \({nS}_{zz}\) (Fig. 3A).

4.2.2 Comparing Ratio of Stresses (\({{\varvec{R}}}_{{\varvec{\theta}}{\varvec{\theta}}}\), \({{\varvec{R}}}_{{\varvec{\theta}}{\varvec{z}}}\), \({{\varvec{R}}}_{{\varvec{\theta}}{\varvec{z}}}\), \({{\varvec{R}}}_{{\varvec{z}}{\varvec{z}}}\)) at SBP, i.e., Combined Effect of Stiffness and Stretch on Stress

At SBP in males, \({R}_{\theta \theta }\) had a higher value compared with \({R}_{\theta z}\) while \({R}_{\theta z}\) had a lower value compared to \({R}_{z\theta }\) and \({R}_{zz}\), (0.053 ± 0.027 vs. 0.027 ± 0.013 vs. 0.047 ± 0.023 vs. 0.048 ± 0.019, respectively, P < 0.05). In females, \({R}_{\theta \theta }\) had a higher value compared with \({R}_{\theta z}\) (0.045 ± 0.020 vs. 0.023 ± 0.010, P < 0.05), while \({R}_{\theta z}\) appeared to be lower than \({R}_{z\theta }\) and \({R}_{zz}\), (0.040 ± 0.017 vs. 0.040 ± 0.015, respectively, P < 0.1). At SBP in males and females, \({R}_{\theta \theta }\) was almost 100% higher compared to \({R}_{\theta z}\), and when comparing \({R}_{\theta \theta }\), \({R}_{z\theta }\) and \({R}_{zz}\), all three were of the same magnitude (Fig. 3B).

Figure 3A shows the impact of stiffness on stress, while Fig. 3B includes stretch and thus shows the combined effect of stiffness and stretch on stress. \({\lambda }_{\theta }\) had a higher value compared to longitudinal stretch \({\lambda }_{z}\) (Table 3). The effect of this difference is not seen in Fig. 3A but is considered in Fig. 3B. Comparing Fig. 3A and B expose the impact of stiffness and stretch on stress. First, \({S}_{\theta z}\) and \({S}_{z\theta }\) were of the same magnitude. Second, combining \({S}_{\theta z}\) and \({S}_{z\theta }\) with respective \({\lambda }_{\theta }\) and \({\lambda }_{z}\) showed that a small change in \({\lambda }_{z}\) will have a lesser effect on circumferential stress (second bar from left in Fig. 3B) compared to the effect a small change in \({\lambda }_{\theta }\) will have on longitudinal stress (third bar from left in Fig. 3B). Hence, it appears that the effect of a change in \({\lambda }_{\theta }\) will have a greater impact on longitudinal stress than a change in \({\lambda }_{z}\) will have on circumferential stress.

Table 3 Stiffness and stretch at different blood pressures for sexes, divided into age groups

4.3 Aortic Wall Stiffness and Stretch

Stiffness and stretch are reported in Table 3.

4.3.1 Comparing Stiffness (\({{\varvec{S}}}_{{\varvec{\theta}}{\varvec{\theta}}}\), \({{\varvec{S}}}_{{\varvec{\theta}}{\varvec{z}}}\), \({{\varvec{S}}}_{{\varvec{z}}{\varvec{\theta}}}\), \({{\varvec{S}}}_{{\varvec{z}}{\varvec{z}}}\)) Between Age Groups Within a Sex at SBP, DBP and MAP

(From left to right in Table 3).

At SBP, DBP and MAP, in males, elderly and middle aged compared with young had a higher \({S}_{\theta \theta }\), \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) (P < 0.05, respectively). In females, elderly compared with young had a higher \({S}_{\theta \theta }\), \({S}_{\theta z}\) and \({S}_{z\theta }\) (P < 0.05, respectively); there appeared to be a higher \({S}_{zz}\) for elderly when compared to young (P < 0.1).

At SBP, DBP and MAP, in males, \({S}_{\theta \theta }\), \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) correlated positively with age (P < 0.05, respectively). This holds true for females as well (P < 0.05, respectively), except for \({S}_{zz}\), which appeared to correlate positively (P < 0.1). Altogether, our findings suggest that \({S}_{\theta \theta }\), \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) increased with age.

Furthermore, elderly compared with young males as well as females showed a higher value at SBP and DBP for \({S}_{\theta \theta }\), \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) (SBP: male 400%, 400%, 350%, 270%, respectively; female 200%, 200%, 150%, 100%, respectively. DBP: male 300%, 300%, 300%, 300%, respectively; female 100%, 100%, 100%, 100%, respectively). Thus, it seems that the increase in stiffness from the young to the elderly was constant and approximately two times higher in males compared with females both at SBP and DBP, i.e., for the physiological pressure range.

For elderly males and females at SBP, \({S}_{\theta \theta }\) compared with \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) had a higher value (90%, 100% and 95–100%, respectively). Hence, independent of sex, the elderly had values where \({S}_{\theta \theta }\) was approximately twofold higher than \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) which, in turn, were of the same magnitude. Independent of sex and pressure, it appears that for the elderly, \({S}_{\theta z}\) had approximately the same value as \({S}_{z\theta }\).

4.3.2 Comparing Stiffness Between Sexes Within an Age Group at SBP, DBP and MAP

(From left to right in Table 3).

At DBP, for the elderly, males compared to females appeared to have higher \({S}_{\theta \theta }\), \({S}_{\theta z}\), \({S}_{z\theta }\) and \({S}_{zz}\) (P < 0.1).

Males compared to females had a lower stiffness as young but a higher stiffness as elderly (Table 3). This suggests two things: first, for the young group, the female abdominal aortic wall may be stiffer compared to that of the male; second: the male abdominal aortic wall grows stiffer at a higher rate with age.

4.3.3 Comparing Stretch (\({{\varvec{\lambda}}}_{{\varvec{\theta}}}\), \({{\varvec{\lambda}}}_{{\varvec{z}}}\)) Between Age Groups Within a Sex at SBP, DBP and MAP

(From left to right in Table 3).

At SBP, DBP and MAP in both sexes, elderly and middle aged compared with young had a lower \({\lambda }_{\theta }\) and \({\lambda }_{z}\) (P < 0.05, respectively). In the model, \({\lambda }_{z}\) was assumed to be constant over blood pressure, [8].

At SBP, DBP and MAP in both males and females, \({\lambda }_{\theta }\) and \({\lambda }_{z}\) correlated negatively with age (P < 0.05, respectively) which suggests that \({\lambda }_{\theta }\) and \({\lambda }_{z}\) decreased with age.

From Table 3, it appears that the difference in \({\lambda }_{\theta }\) between young and elderly at SBP was approximately the same in males and females although at DBP males expressed a higher difference compared to females.

4.3.4 Comparing Stretch Between Sexes Within an Age Group at SBP, DBP and MAP

(From left to right in Table 3).

Males compared with females showed no statistically significant difference in \({\lambda }_{\theta }\) and \({\lambda }_{z}\).

5 Discussion

The main findings of this study were as follows:

  1. 1.

    The effect of a change in circumferential stretch appeared to have a greater impact on longitudinal stress than a change in longitudinal stretch might have on circumferential stress.

  2. 2.

    Circumferential and longitudinal direct and cross-coupled stiffnesses increased while circumferential and longitudinal stretches decreased with age, independent of sex.

  3. 3.

    For the young group, the female abdominal aortic wall is stiffer compared to that of the male while for the elderly group the male abdominal aortic wall seemed to be stiffer than the female wall.

The mechanical forces maintain a balance in the vascular wall through remodeling. A stretch of the vessel wall will cause an increased wall stress. The magnitude of the stress will be determined by the wall stiffness. Biological material such as vascular tissue will have a nonlinear stress–stretch relationship, meaning that the stiffness is nonlinear. Not only shear, circumferential and longitudinal stresses are fundamental mechanical properties in arterial wall remodeling, but the interdependency of circumferential–longitudinal coupling of stresses and stretches is considered to be of importance as well [4, 5, 20].

Our results suggest that cross-coupled stiffness participates in the interdependency between circumferential and longitudinal coupling of stresses and stretches. In independent processes, we propose that an increase in stretch in the circumferential direction affects longitudinal stress through cross-coupled circumferential stiffness and the axial force with longitudinal stretch affects circumferential stress through longitudinal cross-coupled stiffness. Interdependency between circumferential and longitudinal remodeling has been shown where a longitudinally stretched vessel grew into its new length recovering a prestretched stress value, while wall thickness increased but not circumferential stress [21]. Increased cell proliferation of smooth muscle and endothelial cells as well as increased internal elastic laminae fenestrae size, but not density was reported [21]. Increase in wall thickness is a reversible process [22]. An increase in vessel length is not reversible since matrix metalloproteinase (MMP) activity degrades proteins and prevents longitudinal strain from returning to normal values [23].

The effect of a change in circumferential stretch appears to have a greater impact on longitudinal stress than a change in longitudinal stretch might have on circumferential stress. The linearized model showed that an increase in longitudinal stress was composed by approximately 50% of stress depending on cross-coupled stiffness while an increase in circumferential stress was composed by approximately 33% of stress derived from cross-coupled stiffness. This suggests that circumferential stretch may have a significant impact on longitudinal remodeling while longitudinal stretch probably has a minor impact on circumferential remodeling.

The blood pressure-induced alterations of circumferential stretch may, through the proposed mechanism of cross-coupled stiffness, contribute significantly to increased longitudinal stress and stretch, resulting in the aorta growing into its new length. Considering that decreased longitudinal strain might exacerbate a lengthening and create tortuosity [23], there are two different mechanisms working together in an unfavorable way with respect to optimal aortic length. This may, at least in part, offer an explanation to why the aorta is prone to tortuosity.

It is well accepted that a stiffer arterial wall is accompanied by a higher blood pressure, e.g., hypertension [24]. In this context, stiffness is commonly associated with circumferential stiffness. The role of longitudinal stiffness in hypertension is much less investigated as well as cross-coupled stiffness. Our findings suggest that with increasing age, circumferential stiffness gets proportionally higher compared to longitudinal and cross-coupled stiffness while longitudinal and cross-coupled stiffnesses are of the same magnitude. This indicates that the composition of the vessel wall is changed in a stiffer direction by way of loadbearing structures in circumferential and longitudinal directions as well as the interdependency of the circumferential–longitudinal coupling.

Abdominal aortic aneurysms can have complicated geometries and aneurysm diameter is mainly governed by shear stress. Aneurysms rupture when wall stress exceeds wall strength [25]. Under pathological circumstances, direct and cross-coupled stiffness may adopt quite different values compared to ordinary conditions. However, our findings suggest that in an adverse situation, longitudinal stretch could contribute to circumferential stress exceeding wall strength resulting in wall rupture, even though the aneurysmatic diameter is small and/or the presence of a small circumferential stretch.

Our finding of an age-dependent increase in circumferential and longitudinal stiffness as well as decrease in stretch is in accordance with the literature, although we cannot show a difference between sexes in stretch [14, 24]. Furthermore, this age-dependent increase in circumferential stiffness agrees with previous reported findings [10].

For the young group, the female abdominal aortic wall may be stiffer compared to that of the male while for the elderly group the male abdominal aortic wall seems to be stiffer than in the female, indicating that the male abdominal aortic wall may grow stiffer at a higher rate with age. The gain in male abdominal aortic wall stiffness with age is in accordance with observations of a more accelerated ageing process in the aorta in males compared to females [26]. Female sex hormones seem to protect against aortic wall elastolysis, which in turn will contribute to vessel wall stiffening [27]. The effect of postmenopausal hormone replacement therapy in females seems to reduce arterial stiffness [28] further suggesting that females are better protected against aortic wall stiffening than males.

Different mechanisms contribute to the age-related alterations of the vascular wall. Increased wall stress has been suggested to induce synthesis of collagen in the fiber direction and thus a maintained stretch of the individual fibers [11, 29, 30]. The orientation of collagen fibers may also be important since it has been found to be more helical in the media and more longitudinal in the adventitia [31, 32]. Additionally, increased glycation of elastin and collagen as well as changed isoforms of collagen in the aortic wall might contribute to the age-related increase in stiffness [33].

It has been pointed out that constitutive equations describing biological material such as the arterial vessel wall suggest a circumferential–longitudinal coupling of stresses and stretches [5, 17, 20]. Our model expresses such an interdependency with a circumferential–longitudinal coupling where both circumferential and longitudinal stresses depend on stretch in circumferential and longitudinal directions (\({\sigma }_{\theta }({\lambda }_{\theta },{\lambda }_{z})\) and \({\sigma }_{z}({\lambda }_{\theta },{\lambda }_{z})\)) [8]. An attempt to quantify this interdependency was made through a linearization of the model equations. To the authors knowledge, there has been no in vivo quantification of the remodeling characteristics in one direction caused by stretching in cross-directions. Biaxial tests ex vivo provide some information regarding the interconnection of forces. Nevertheless, as a tool for assessing in vivo mechanisms, PIMMP may be used to explore mechanical properties of the vessel wall and may serve as a potential tool for risk assessment for the development of cardiovascular disease.

Based on our results, we propose that cross-coupled stiffness may play a part in the age-related remodeling of the aortic wall, where an increase in circumferential cross-coupled stiffness induced longitudinal stress constitute 50% of the change in total longitudinal stress and therefore may contribute to aortic tortuosity. Furthermore, with increasing age it seems that the interdependency of circumferential–longitudinal coupling increases possibly due to a change in composition of the aortic wall which may be a contributing factor in hypertension. Additionally, under adverse conditions our findings indicate that, e.g., a longitudinal stretch could contribute to circumferential stress increasing the likelihood of aneurysm rupture.

6 Limitations

Since there are differences between different parts of the vascular system and the abdominal aorta as well as between different segments of the aorta, with respect to histology and pulse-wave velocity, it must be emphasized that our findings should be extrapolated with caution to other segments of the aorta. Furthermore, assumptions in the model can affect the use of PIMMP.

The parameter identification algorithm with the attendant mechanical model was validated against finite element (FE) models of an artery, a procedure which is also known as “in silico” validation. The validation showed good agreement between PIMMP and the FE models [9]. In particular, the longitudinal response in the simulation is sensitive to the longitudinal prestretch (λz), and thus, this parameter is a natural candidate for the validation procedure. The PIMMP result shows a longitudinal prestretch within the range which has been reported from autopsy study of the human abdominal aorta [14]. In general, nonlinear arterial models are difficult to validate experimentally since each model parameter needs to be changed independently of the other parameters. Therefore, our model has not been validated against measurements in vivo.

The strain–energy function used in this study is based on Holzapfel et al. [12] (Eq. 10). Compared to Schulze-Bauer and Holzapfel [7], it allows for a better fit to young subjects, particularly in the low-pressure region where the mechanical behavior is primarily determined by isotropic material components such as elastin and the collagen recruitment is small [8].

The mechanical model is based on an assumption of the cylindrical wall being a membrane, i.e., wall thickness should be negligible when compared with the radius. However, wall thickness-to-radius ratio is ~ 0.1–0.2 for the abdominal aorta; therefore, the validity of the assumption might be argued [11]. Furthermore, the aortic wall consists of three distinct layers which all have different mechanical properties [34]. An exact model should consider this. However, such a high resolution in the model might introduce dependencies among the parameters during parameter identification [15]. Regarding the parameters from the membrane model as averages, it might be thought of as describing the global response of the three layers of the aortic wall.

Our model concerns the passive mechanical properties of the aorta and assumes that the outer boundary is traction free, i.e., the artery is a freestanding tube with neglected periadventitial support. It has been shown that the difference in circumferential stress between the freestanding and the tethered states is small, approximately 10% [35]. Considering the limited effect of tethering and that such mechanism would require the identification of the mechanical properties of surrounding tissue, we decided that a freestanding approach would be sufficient for this study.

7 Conclusion

This paper presents quantitative estimates for circumferential direct and cross-coupled as well as longitudinal direct and cross-coupled stiffnesses for the human abdominal aorta stratified for age and sex, based on in vivo and in situ measured radius and pressure. These stiffnesses have so far only been quantified through measurements ex vivo, ex situ in the laboratory. Our findings of an age-related circumferential remodeling acting through cross-coupled stiffness might have a major impact on longitudinal stress and remodeling and possibly aortic tortuosity. The findings suggest an age-related change in wall constituents. A potential explanation, although not studied here, could be related to chemical alterations of wall constituents and changes in the physical distribution of fibers.