Findings from basic science and clinical studies agree that arterial stiffness is fundamental to both the mechanobiology and the biomechanics that dictate vascular health and disease. There is, therefore, an appropriately growing literature on arterial stiffness. Perusal of the literature reveals, however, that many different methods and metrics are used to quantify arterial stiffness, and reported values often differ by orders of magnitude and have different meanings. Without clear definitions and an understanding of possible inter-relations therein, it is increasingly difficult to integrate results from the literature to glean true understanding. In this paper, we briefly review methods that are used to infer values of arterial stiffness that span studies on isolated cells, excised intact vessels, and clinical assessments. We highlight similarities and differences and identify a single theoretical approach that can be used across scales and applications and thus could help to unify future results. We conclude by emphasizing the need to move toward a synthesis of many disparate reports, for only in this way will we be able to move from our current fragmented understanding to a true appreciation of how vascular cells maintain, remodel, or repair the arteries that are fundamental to cardiovascular properties and function.

## Introduction

Mechanical factors have long been known to play vital roles in arterial health, disease, and treatment. Studies dating back to the late 19th century showed that altered blood flow and pressure typically result in changes in luminal radius and wall thickness [1,2], thus foreshadowing the need to quantify the mechanoregulation of arterial structure and function. It was not until the 1970s, however, that it was discovered how such changes occur. Experiments on cultured endothelial cells and vascular smooth muscle cells revealed flow- and stretch-induced changes in secreted proteins that resulted from altered gene expression [3,4]. Hence was born the field of vascular mechanobiology, reviews of the first three decades of which can be found elsewhere [5,6]. Importantly, careful studies of the aorta during this same defining period revealed that circumferential lamellar tension (∼2 N/m), and by inference circumferential lamellar stress (∼150 kPa), is similar within the medial layer across multiple mammalian species [7]. This finding, coupled with subsequent observations that aortic material stiffness at mean arterial pressure is also similar (∼500 kPa) across invertebrates and vertebrates [8], strongly suggests the existence of a “mechanical homeostasis”—that is, vascular cells seek to establish and then maintain particular mechanical quantities near target (homeostatic) values. Indeed, findings at subcellular, cellular, tissue, and organ levels suggest that such homeostatic processes exist across many spatial and temporal scales [9].

Two seminal contributions by Fung and his colleagues advanced our ability to quantify arterial wall stress and the associated stiffness, which are important both for assessing the mechanics of the wall and understanding its mechanobiology. First, Fung observed via uniaxial experiments on soft tissues that material stiffness (in this case, a change in the first Piola-Kirchhoff stress with respect to stretch) relates nearly linearly to the stress itself [10]. Importantly, this finding suggests directly that the first Piola-Kirchhoff stress increases nearly exponentially with stretch. Indeed, a similar observation had previously been reported for the overall structural stiffness of the pressurized eye based on its pressure–volume relation [11]. Motivated by these findings, Fung later hypothesized the existence of an exponential stored energy function $W$ that yields an exponential relationship between the second Piola-Kirchhoff stress $S$ and Green strain $E$, where $S=\u2202W/\u2202E$. Now known as Fung elasticity, this hyperelastic function can be written as $W=c(exp(Q)\u22121)$, where $c$ is a material parameter and the scalar function $Q$ depends quadratically on $E$ [12]. An associated metric of material stiffness, in referential form, is thus $C=\u2202S/\u2202E=\u22022W/\u2202E\u2202E$. Notwithstanding the importance of Fung's constitutive hypothesis in quantifying the mechanical behavior of many soft tissues and solving associated initial-boundary value problems, the observation that material stiffness relates linearly to stress implies that it should be challenging to delineate whether mechanobiological responses correlate better with stress or stiffness, consistent with the aforementioned observations that homeostatic values of aortic stress and material stiffness are similar across species [7,8]. Second, Fung showed that the existence of residual stress in arteries dramatically affects the calculated distribution of Cauchy stress $t$ across the arterial wall [13], which in combination with basal smooth muscle cell tone results in a nearly homogenized transmural distribution [12]. This finding supports the concept that vascular cells seek to establish and then maintain mechanical stress or stiffness near a homeostatic target, regardless of location within the wall. Importantly, it now seems that different homeostatic targets may exist for different cell types that populate different layers of the wall [14], thus reinforcing the concept of a cell-specific mechanical homeostasis that can manifest at the vessel level as well.

The importance of arterial stiffening in human disease was anticipated as early as the turn of the 20th century by Sir W^{m} Osler, one of four founding fathers of Johns Hopkins medicine, but confirmation had to await a seminal clinical study at Paris, France that revealed that an increased pulse wave velocity ($PWV$) is an initiator and indicator of diverse cardiovascular, neurovascular, and renovascular disease and thus all-cause mortality [15]. Carotid-to-femoral $PWV$ is now considered to be the gold standard for clinical assessment of arterial stiffness [16], which increases with diabetes, hypertension, natural aging, particular connective tissue disorders, and so forth [17–19]. This metric of stiffness is best measured by dividing the vascular centerline distance between two recording sites by the time that it takes the (foot of the) pressure pulse wave to travel between these sites (“foot-to-foot” $PWV$). Intuitive understanding of the $PWV\u2009$ often comes from the Moens–Korteweg equation, derived in the 1870s and commonly written as $PWV=Eh/2\rho a$, where $E$ is an intrinsic (linear) material stiffness, $h$ is the wall thickness, $a$ is the luminal radius, and $\rho $ is the mass density of the fluid. Understood in this way, we see that $PWV$ depends on both the intrinsic material stiffness and the geometry of the wall, hence rendering it an integrated (along the centerline distance) measure of structural, not material, stiffness.

Not surprisingly, there continues to be increasingly greater interest in measuring and understanding arterial stiffness in basic science studies and clinical assessments. Quantifying stiffness helps us to understand critical questions related to, among other issues, fundamental aspects of vascular mechanobiology, effects of a distensible wall on the hemodynamics, disease progression, and even the efficacy of particular clinical treatments. In this paper, we briefly review and contrast different methods for measuring and interpreting arterial stiffness both to emphasize that which has been learned and that which remains unclear. In particular, we emphasize that many different metrics of arterial stiffness reported in the literature have different theoretical underpinnings and thus different meanings. Caution is thus warranted when comparing values of stiffness across studies as well as when interpreting fundamental implications of a particular metric on either the cell biology or the systems physiology.

## Approaches

We first review the theoretical basis for some of the experimental methods that are commonly used to infer material properties of the arterial wall and its constituent parts. The presentation is organized by physical scale, not chronological introduction within the field. Regardless of scale, it is critical to delineate concepts of material versus structural stiffness, and similarly, it is essential to note that values of stiffness depend on the type (tensile, compressive, or shear) and magnitude (small or large) of the loading as well as the mechanical state in which the experimental loads are imposed, especially if the sample is otherwise unloaded during testing. Unfortunately, many papers often simply report a value for stiffness without emphasizing the precise definition or experimental conditions, hence obscuring the range of applicability. Here, we attempt to delineate some of these issues.

### Atomic Force Microscopy.

Invented in the early 1980s at IBM, atomic force microscopy (AFM) uses a laser to detect small deflections of a cantilever probe as it interacts with the surface of a specimen. The atomic force microscope can be used in various modes to assess the topography of a surface or to probe mechanical properties. With regard to the latter, one can perform precise indentation (compression) tests or functionalize the probe to enable tractions (tension) to be applied at the surface. As implied by its name, the most precise atomic force microscopes measure nanoscale forces and motions, thus it is not surprising that AFM has been used extensively in biology, biophysics, and bioengineering, often to assess aspects of cell stiffness, cell-matrix interactions, or local matrix stiffness. With larger probes, one can also assess bulk tissue stiffness. Associated data consist primarily of the applied force (inferred from knowledge of the cantilever properties and deflections) and the imposed motion of the probe, often the depth of penetration into the surface. Quantities such as stress and stiffness can be determined by solving the associated initial-boundary value problem.

where 2$ro$ is the diameter (*d _{o}*) of a flat-ended cylindrical probe and $\lambda $ is the equibiaxial in-plane stretch. A Hertz-type solution can be recovered as $\lambda \u21921$, for which $\alpha \u219216ro\mu \u2009=\u20094doE/3$, with Young's modulus $E(=6\mu )$ having units of N/m

^{2}. Hence, $E$ can be computed easily given the diameter of the probe and the slope of the force-depth data ($f,\delta $) for an otherwise unloaded specimen. Of course, because of inherent uncertainty in all experimental data, the stiffness parameter(s) should be estimated using least-squares regressions of data via an over-determined system of equations. Moreover, in the case of a nonlinear behavior, the intrinsic material stiffness (e.g., shear modulus $\mu $ for a neo-Hookean response) should be inferred from force-depth data over a range of in-plane stretches that are relevant to the biological or physiological regime of interest in order to characterize well the overall material behavior. Noting that adherent cells spread out and develop significant cytoskeletal stress when cultured on any substrate, which is to say that they develop a finite nearly in-plane state of stress or stretch, it has been shown that a small-on-large framework can be used to assess cell stiffness from AFM for more general exponential stored energy functions [24].

Indeed, although a neo-Hookean relation admits a simple analytical solution that allows one to intuit effects of finite in-plane stretches on AFM measured stiffness (Eq. (1)), the Fung exponential model is more appropriate for most soft tissues and even cells. Appendix B summarizes a general Fung model and Appendix C provides results for a small-on-large solution for AFM based on an isotropic Fung exponential model. As expected, the numerical implementation is straightforward. Nevertheless, results from AFM for assessing vascular cell stiffness [25,26] and arterial properties [27–30] typically have been interpreted in terms of a single Young's modulus inferred from the Hertz solution, which tacitly ignores the important nonlinear dependence on in-plane prestretch and the nonlinear (often exponential) material behavior. Hence, the inferred values of stiffness are typically well below in vivo values and could be misleading regarding in vivo mechanobiology.

Although not discussed in detail here, additional methods are also used to assess cell stiffness, including magnetic twisting cytometry (MTC) and optical tweezers [31]. MTC is similar to AFM except that one either affixes onto or embeds within a living cell a ferromagnetic bead that can be cyclically twisted using a magnetic field. Associated data are often interpreted in terms of the so-called storage ($G\u2032$) and loss ($G\u2033$) shear moduli, which are basic descriptors of a one-dimensional (1D) linearly viscoelastic behavior over small strains and rotations [32]. For example, MTC was used to show that the stiffness of vascular smooth muscle cells (i.e., elastic storage and dissipative loss moduli) increases in aging [33]. Similar to the Hertz solution, however, this inference of material properties does not explicitly account for the underlying nonlinear material behaviors or the finite deformations that adherent cells experience in vitro when spreading on a surface or to which they are likely exposed in vivo. Again, a small-on-large approach for data analysis could be more appropriate.

### Biaxial Biomechanical Testing.

where $r\u2208[ri,ro]$, with $ri$ and $ro$ being the inner and the outer radius of the sample, respectively, $P$ is the distending pressure, and $L$ is the reduced axial load. If only the mean (i.e., radially averaged, $\u2329\u2026\u232a$) values of stress are of interest, then these two integral relations can be replaced with algebraic ones, $\u2329t\theta \theta \u232a\u2261\sigma \theta \theta =Pri/h$ and $\u2329tzz\u232a\u2261\sigma zz=(L+\pi ri2P)/\pi h(2ri+h)$, where $h=ro\u2212ri$ is the wall thickness. Mean values of stress are surprisingly useful because of the aforementioned effect of residual stress in homogenizing the transmural distribution of stress. Regardless, some investigators infer stiffness from plots of stress versus stretch or strain (yielding so-called tangent moduli), but stiffness depends on the full multiaxial deformation and is best computed from an appropriate constitutive relation as noted above for the referential stiffness ($\u2202S/\u2202E$). Given the complex microstructure of the vascular wall, most investigators now prefer the so-called fiber-family constitutive models [35–37] to characterize multiaxial data, though the Fung exponential provides good fits to data in many cases (recall that Appendix B reviews a general orthotropic Fung relation). Best-fit values of material parameters for any appropriate constitutive relation can be determined via nonlinear regression of pressure-diameter and axial force-length data, with data from multiple biaxial protocols typically combined to improve the parameter estimation, again via (nonlinear) least squares regression. In this regard, it is important to note that many studies nevertheless report only pressure-diameter data at a single value of axial stretch and often do not measure the associated axial force. Such data are essentially one-dimensional and not useful for calculating in vivo relevant biaxial stress or stiffness [38].

where $\delta ij$ is the Kronecker delta and a superscript $o$ denotes a deformation quantity that is associated with an original finite deformation about which the small superimposed deformation occurs. For example, the large deformation could be for an artery from a traction-free reference configuration to a finitely distended and extended in vivo configuration near mean arterial pressure about which relatively small motions occur over a cardiac cycle. Second, there is a pressing need to understand better how central arterial stiffness and resistance vessel function interact to affect both local and global hemodynamics and cardiovascular function [19,39]. Third, the effects of perivascular tethering can be as important as the structural stiffness (i.e., material stiffness and geometry) in affecting the hemodynamics though there has not been much quantification of this effect [40,41]. Fourth, constitutive relations that are inferred in vitro are generally established under quasi-static conditions, whereas in vivo loading is pulsatile and hence intrinsically dynamic. Associated differences in strain-rate may influence calculated values of measured stiffness [42,43]. Finally, although well suited for nearly cylindrical samples, standard distension–extension tests are not sufficient for more complex arterial geometries, particularly those manifesting in disease. Fortunately, new methods are emerging that enable local material properties to be inferred using full-field strain measurements and inverse methods for material characterization [44,45].

Although not discussed in detail here, there are also many reports of vascular function and properties based on either uniaxial loading tests (e.g., Ref. [46]) or ring myography (e.g., Ref. [25]). Originally conceived to study isometric contraction, ring tests are performed by placing a short ring-like sample of the vascular wall on two mounting fixtures that are separated by a finite distance. The sample thus stretches primarily in the circumferential direction as the fixtures are separated and the sample deforms from a circular to an oval to a more uniaxial geometry within the central region of measurement. Measurement of the force acting on (passive) or generated by (active) the sample thus allows one to estimate 1D stress-stretch information or to construct dose-response curves while holding a sample at a fixed separation distance (isometric). This approach is particularly valuable for high-throughput comparisons of different drugs or their doses, but the associated contraction of the smooth muscle is under nonphysiological loading [47]. Inference of passive or active stiffness from these ring tests is further compromised by the lack of an exact solution to the full boundary value problem (finite bending and uniaxial extension of an annulus) and the lack of biaxial loading [12,34]. In particular, both the radial and axial stretches reduce below one as circumferential stretch increases with the separation of the mounting fixtures. Axial stretch is typically much greater than unity in vivo and is important to both the mechanics and the mechanobiology [48]. Finally, in-plane biaxial tests can be performed on excised samples (e.g., Ref. [49]) and are generally very informative for they can largely mimic the in vivo state of stress, with the exception of radial compression due to a distending blood pressure. The primary caveat with in-plane biaxial tests is that large specimens can retain some of their natural cylindrical curvature in the unloaded state, thus requiring combined bending and extension to load the sample biaxially. Uniaxial tests performed on different samples with different gross orientations can avoid some issues of residual curvature, but associated protocols are limited and cannot explore directly the inherent coupling that is important in dictating the multiaxial mechanical behavior.

### Pulse Wave Velocity and Distensibility.

Of course, $PWV$ can also be estimated using computational methods given much more information: the spatially distributed geometry, material properties, and boundary conditions. Such calculations avoid implicit assumptions inherent in the Moens–Korteweg and related equations and they enable informative parametric studies, as, for example, how vascular taper or distal resistances affect the pulse wave. Unsteady 1D and three-dimensional models have been used to compute $PWV$ [57,58], but again additional intuitive insight can be gleaned from analytical solutions. Toward this end, the theory of small deformations superimposed on large has also been used to obtain analytical results for a simple cylindrical geometry [59]. This solution shows explicitly that both the diastolic distension and the axial prestretch are important contributors to the computed values of $PWV$; in the limit as the strains become small, this relation recovers the Moens–Korteweg equation. As in the case of the AFM, however, this analysis shows clearly that complexities due to nonlinear material behaviors and finite deformations play important roles in determining the precise value of stiffness that is computed.

## Comparisons

Fung and colleagues used nonlinear regression to determine best-fit values of the seven material parameters in an exponential stored energy function $W$ (Eqs. (A1) and (A2) in Appendix B) from both uniaxial (radial compression [60]) and biaxial (pressure-diameter-axial force-stretch [13]) data for a rabbit thoracic artery. Not surprisingly, the parameter values differed for the predominantly compressive versus tensile behaviors. Figure 1 compares these results together; note the anisotropy and strong nonlinearity in the tensile (circumferential and axial) behaviors and the asymmetry between the tensile and compressive behaviors even at modest strains. Results are not shown for greater compressive strains since the original tests used modest levels of compression and the associated analytical solution reveals potential bifurcations in the equilibrium solutions at higher compressive strains. Importantly, given that the slope of the stress-stretch curves reflects (but does not define) the stretch-dependent spatial material stiffness, note the expected tremendous differences between the low values of material stiffness near the unloaded configuration (low stress) and high values near the in vivo configuration (physiologic stress). Figure 2 shows calculated circumferential behaviors for different degrees of fixed axial stretch from one to the in vivo value (1.691 for this particular rabbit artery). This result reveals the strong biaxial coupling, with circumferential stress and stiffness affected dramatically by the value of axial stretch. The aforementioned ring myography and uniaxial tests disregard such coupling and thus can underestimate the actual stiffness dramatically.

Figure 3 shows calculated results for a simulated AFM indentation test with the material modeled with the same Fung-exponential type of constitutive relation and the indentation performed at different levels of fixed equibiaxial in-plane stretch. Since the analytical small-on-large solution holds for isotropic material behaviors [21], we first estimated new values of the material parameters in the Fung-exponential that yield an isotropic response similar to the mean of the anisotropic response shown in Fig. 1. Again, it can be seen that there is a strongly coupled response between different directions of loading, here out-of-plane versus in-plane. Most published works on arterial wall and vascular cell stiffness disregard this coupling effect and use AFM to test samples, in the absence of a pre-existing in-plane stress, to compute Young's modulus that strictly holds only for small strains. Such a situation is not physiological and again is expected to underestimate the actual in vivo stiffness dramatically.

To emphasize further both the multiaxial nature of stress and stiffness, and couplings therein, we show three-dimensional plots in Fig. 4 for simple cases of combined circumferential and axial extensions of Fung elastic materials. That is, assuming principal homogeneous deformations and incompressibility, the stored energy function and thus associated wall stress and stiffness depend on the in-plane principal stretches alone. The primary in-plane components of the mean Cauchy stress can thus be thought of conceptually in terms of components of the left stretch tensor $V=FFT$, namely $\sigma \theta \theta =\sigma \u0302\theta \theta (\lambda \theta ,\lambda z)$ and $\sigma zz=\sigma \u0302zz(\lambda \theta ,\lambda z)$. As it can be seen in the figure, for both isotropic (left columns; cf. Fig. 3) and anisotropic (right columns; cf. Figs. 1 and 2) behaviors, circumferential and axial deformations contribute similarly to the overall elastic energy storage and thus stress and material stiffness, here in a spatial rather than referential description since spatial quantities (defined relative to the current configuration) are most relevant in vivo.

## Discussion

It is now widely accepted that both cell- and matrix-level stiffness are fundamental to mechanobiological responses within the vasculature, including modulation of cell phenotype. Similarly, it is widely accepted that tissue-level stiffness is fundamental to the hemodynamics, particularly propagation of the arterial pressure wave that dictates many aspects of end organ health or disease. There is, therefore, an appropriately growing literature on arterial stiffness. Yet, the many different reported metrics have different meanings because of different types of loads that are imposed relative to different biomechanical states. We summarize in Table 1 some of these methods and metrics and, for purposes of illustration, we list in Table 2 some values of stiffness resulting from these different methods. As it can be seen, reported values of stiffness (some compressive, some shear, most tensile) differ by orders of magnitude, as should be expected for highly nonlinear material behaviors when assessed relative to different configurations, ranging from otherwise unloaded to in vivo relevant. Assuming that the associated calculations were performed correctly, each of these different values of stiffness should be viewed as reliable. The key question, which we do not attempt to answer here, is therefore: How can we extract from these disparate metrics, having different meanings, a unified understanding of vascular cell mechanobiology and biomechanics and their roles in dictating vascular health or disease progression? We submit that focused effort should be directed toward answering this question, for without such our ability to use basic science findings to inform clinical decisions will remain largely wanting since understanding will remain largely fragmented and incomplete.

In some ways, we have emphasized the obvious—the value of stiffness depends on its definition (material versus structural), the configuration to which it refers (referential-unloaded, versus spatial-current), and the conditions under which it is evaluated (compression versus tension versus shear in an otherwise unloaded state or not). Yet, we are not aware of a prior consistent discussion of methods used in vascular mechanics across scales from atomic force microscopy to biaxial tests on cylindrical segments to in vivo measurements. As we have noted, a number of metrics reported in the literature are based on solutions from classical elasticity because of the associated simplicity, not the theoretical relevance. Strictly speaking, therefore, most of these results are applicable only for linearly elastic isotropic responses under small strains and rotations when measured about an unloaded state—conditions not applicable to mechanobiologically, physiologically, or clinically relevant situations. Although we noted two analytical examples wherein a material stiffness applicable to a nonlinear behavior can be compared directly in the limit to a small strain Young's modulus [21,59], it is generally problematic and uninformative to compare findings based on Hertz, Moens–Korteweg, or similar equations with those based on nonlinear constitutive relations or appropriate linearizations thereof, including the theory of small deformations superimposed on large. Yet, values of Young's modulus continue to be reported and, in some cases, results are similarly presented based on concepts from linearized viscoelasticity, including storage and loss moduli. Of course, one could argue that such results can be insightful when one consistently compares values across studies using the same methods and metrics whether they strictly hold theoretically or not in configurations relevant to in vivo conditions or not. For example, a compressive Young's modulus inferred from AFM has been reported by multiple groups to increase for isolated smooth muscle cells in aging and hypertension relative to that in normalcy [25,69], hence these data are trying to tell us something—the question is, What? That is, although results may be reproducible and reliable, the key question should focus on their possible relevance in vivo. Similar methods have been used to show that the nuclear protein lamin-A scales with tissue stiffness [70], with tissues ranging from compliant (adipose or liver) to stiff (ligaments and bones). Again, the results are reproducible, reliable, and provocative, though not evaluated at in vivo values of stiffness. How such results should be interpreted or compared to values that are relevant to the in vivo condition remains an open question.

We submit that in vivo relevant conditions and metrics should be used when possible. Toward this end, the theory of small deformations superimposed on large deformations [71] can serve as a theoretically appropriate method to compute linearized values of stiffness while accounting for underlying nonlinear material behaviors and in vivo relevant finite deformations, with applicability including interpretations of AFM data [21], defining biaxial material stiffness in intact excised vessels for simulations of hemodynamics [35], and computing $PWV$ either analytically [59] or numerically [58]. With regard to the need to measure and report metrics that are relevant to the in vivo conditions, it would be prudent to remember the words of Y. C. Fung written ∼50 years ago. First, “The main difficulty [problem] lies in the customary use of infinitesimal theory of elasticity to the media which normally exhibit finite deformations” [10] and “the greatest need lies in the direction of collecting data in multiaxial loading conditions and formulating a [constitutive] theory for the general rheological behavior of living tissues…” [72].

In conclusion, it is critical to quantify arterial stiffness—material and structural—because deviations from normal values associate with both the phenotypic modulation of vascular cells and the clinical severity of disease or disease risk. That said, there is also a fundamental conceptual issue that must be considered carefully as we seek to advance our understanding of the underlying mechanobiology. Although vascular cells clearly attempt to establish and then maintain certain mechanical quantities near homeostatic values [9,73], the continuum quantities of stress and strain, and metrics such as material stiffness that are derived from them, are actually mathematical concepts, not physical realities [74]. Hence, even though it appears that mean wall stress and stiffness are normally regulated near homeostatic targets across mammalian species [7,8,75], we should not expect a cell to necessarily respond to a stress (i.e., a linear transformation, or tensor, that transforms an outward unit vector into a traction vector at a point). Rather, it is more likely that forces acting at the molecular level change the conformation of important biomolecules and thereby stimulate cell signaling and downstream gene products. There is, therefore, a pressing need to understand better the micromechanics of mechanosensing by cells and the associated mechanoregulation of matrix [76] and to associate such phenomena with convenient continuum metrics such as stress and stiffness. Nevertheless, until, and possibly after, such multiscale understanding is achieved, direct correlations of mechanobiological and (patho)physiological responses with wall stress and stiffness should continue to be identified. Toward this end, an increased use of concepts of nondimensionalization and allometric scaling [77] should also become a priority. There is, therefore, a need for continued development of new concepts and techniques in vascular mechanics and mechanobiology and we conclude with words of Fung in his foreword to the inaugural issue of the journal *Biomechanics and Modeling in Mechanobiology* in 2002 [78]—“let us enjoy the work.”

## Acknowledgment

JDH also gratefully acknowledges the Y.C. Fung Young Investigator Award from ASME that he received in 1990.

## Funding Data

US NIH (Funder ID: 10.13039/100000002; Grant Nos. R01 HL105297, P01 HL134605).

Netherlands Organization for Scientific Research (NWO) (Funder ID: 10.13039/501100003246; Rubicon Grant 452172006).

### Small on Large

The polymath A. Cauchy knew equations of nonlinear elasticity in the 1820s, but analytical solutions to this class of problems had to await the semi-inverse approach of R. Rivlin in the late 1940s. Because of the inherent complexities, analytical solutions to problems in nonlinear elasticity were yet possible for only a relatively small class of motions [79] and many employed the advances in finite element methods that soon arrived [80]. It is often difficult to develop intuition from numerical solutions of highly nonlinear problems, however, thus the concept of “small deformations superimposed on large deformations” became useful in extending the range of possible analytical solutions. As noted herein, small-on-large solutions have been found useful in interpreting experimental results associated with atomic force microscopy [21], quantifying mechanical properties from biaxial tests on excised arteries [35], and relating the structural stiffness that can be inferred in vivo from measurements of pulse wave velocity to the material stiffness inferred from in vitro tests [59]. It seems appropriate, therefore, to briefly outline steps of this approach; the interested reader is referred elsewhere for more details [71].

Briefly, let the location of a material particle in an original configuration be denoted by $X\u2009$ and in a finitely deformed configuration by $x$. In addition, let the location of this particle in a configuration that is close to the finitely deformed one be denoted by $y$. The total deformation gradient is thus $F=(\u2202y/\u2202x)(\u2202x/\u2202X)$, which can be written as $F=FsFo$, where subscript $o\u2009$ denotes the original finite deformation and subscript $s$ denotes the superimposed small deformation. It proves convenient to write $y=x+u$, where ** u** is a small displacement vector. Hence, $Fs=\u2202x/\u2202x+\u2202u/\u2202x$ which we write as $Fs=I+h$, where $h$ is a displacement gradient with $h=tr(hhT)\u226a1$. The right Cauchy–Green tensor can then be written as $C=FTF=[I+hFo]TI+hFo=FoTFo+FoThTFo+FoThFo+FoThThFo\u2245Co+FoThT+hFo=Co+FoT2\epsilon Fo$ if we neglect terms that are higher order in the (small) displacement gradient and we recognize the small strain tensor $\epsilon =(h+hT)/2$. Hence, the right Cauchy–Green tensor in the final configuration simply equals that in an intermediate (finitely deformed) configuration, denoted by subscript

*o*, plus a part that is linear in the small displacement gradient.

Just as the total deformation gradient $F=I+hFo=\u2009Fo+hFo$ consists of a contribution that is due to a finite deformation to an intermediate configuration plus an addition, one can also assume that the total second Piola-Kirchhoff stress can be computed as $S=So+S*$, where the first contribution is associated with the original finite deformation and the second with a superimposed small deformation, namely $S*\u2245(\u2202S/\u2202C)o:C*$ with $S=\u2202W/\u2202Co$ and $C*=C\u2212Co=FoT2\epsilon Fo$ from above. Importantly, it can be shown that we can write the total Cauchy stress $t$, for an incompressible response and to order $h$, as [35] $t=to+ht\u0302o+t\u0302ohT\u2212p*I+FoS\u0302*FoT$, where $t\u0302o$ is the extra (i.e., deformation-dependent) part of the Cauchy stress in the finitely deformed intermediate configuration and $p*$ is a Lagrange multiplier that enforces isochoric motions during the small superimposed deformation. Finally, using the result above for the second Piola-Kirchhoff stress in terms of the stored energy function $W$, and recognizing that a small displacement gradient can be written in terms of the infinitesimal strain and rotation (i.e., $h=h+hT/2+(h\u2212hT)/2$, or, $h=\epsilon +\Omega $), we arrive at our final result: $t+p*I=to+C\epsilon +D\Omega $, where the spatial material stiffness (fourth order) tensor $C$ can be written for an artery as given in component form in the third equation in the main text. Noting that the infinitesimal rotation tensor $\Omega $ vanishes with assumed principal deformations, this derivation reveals clearly that the spatial material stiffness associated with a superimposed small deformation depends directly and strongly on the prestress (or prestretch) and the nonlinear properties represented by $W$, which are amplified nonlinearly by the initial finite deformation when evaluated in the current (spatial) configuration [35].

### Fung Elasticity

##### Uniaxial Compression.

Simulations of the compression experiment were then performed using reported values of the constitutive parameters for a rabbit thoracic artery (sample 11, incompressible case) given in Ref. [60]: $c=43.12\u2009kPa$, $b1=0.8230$, $b2=0.9125$, $b3=1\xd710\u22127$, $b4=1.1237$, $b5=0.4125$, and $b6=0.3768$.

##### Biaxial Stretching.

Simulations of tensile loading were performed using values of the constitutive parameters, again for a rabbit thoracic artery [13]: $c=22.40\u2009kPa$,$\u2009b1=1.0672$, $b2=0.4775$, $b3=0.0499$, $b4=0.0903$, $b5=0.0585$, and $b6=0.0042$. Note the differences between the values for compression and tension despite the same functional form of $W\u2009$ and the same sample tested by the same investigators.

Increasing values of $ri$ thus yield desired increases in pressure (e.g., up to 120 mmHg).

### Simulations of Atomic Force Microscopy for Fung Elasticity

Given $W(Q)$ written in this form, the equations in Ref. [21], written as functions of $\u2202W/\u2202IC$ and $\u2202W/\u2202IIC$, can be used directly to simulate AFM indentations for an isotropic Fung elastic material, similar to the calculation discussed in the text for the simpler neo-Hookean material.