- 1Medical Research Institute, Pusan National University, Busan, South Korea
- 2Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
- 3Department of Convergent Biosystems Engineering, Sunchon National University, Suncheon, South Korea
- 4Institute of Human Harmonized Robotics, Sunchon National University, Suncheon, South Korea
The simulation of the cardiovascular system and in silico clinical trials have garnered attention in the biomedical engineering field. Physics-based modeling is essential to associate with physical and clinical features. In physics-based constitutive modeling, the identification of the parameters and estimation of their ranges based on appropriate experiments are required. Uniaxial tests are commonly used in the field of vascular mechanics, but they have limitations in fully characterizing the regional mechanical behavior of the aorta. Therefore, this study is aimed at identifying a method to integrate constitutive models with experimental data to elucidate regional aortic behavior. To create a virtual two-dimensional dataset, a pair of uniaxial experimental datasets in the longitudinal and circumferential directions was combined using a one-to-many correspondence method such as bootstrap aggregation. The proposed approach is subsequently applied to three constitutive models, i.e., the Fung model, Holzapfel model, and constrained mixture model, to estimate the material parameters based on the four test regions of the porcine thoracic aorta. Finally, the regional difference in the mechanical behavior of the aorta, the correlation between the experimental characteristics and model parameters, and the inter-correlation of the material parameters are confirmed. This integrative approach will enhance the prediction capability of the model with respect to the regions of the aorta.
Introduction
Recently, the simulation of the cardiovascular system and in silico clinical trials have garnered attention in biomedical engineering. Virtual patient cohorts and digital twins are emerging as promising strategies in precision medicine, although they are not yet fully established for clinical applications (Niederer et al., 2020; Chakshu et al., 2021). In physics-based modeling, a virtual model comprises two layers, where one layer includes a set of clinical features, and the other includes a set of model parameters. It is important to associate the model with physical or clinical features. Inadequate anatomical, physiological, and functional parameters of the heart and vasculature hinder the further development of cardiovascular devices and novel treatments (Hose et al., 2019).
The constitutive model represents the behavior of materials in the form of relationship between strain and stress. To establish physics-based modeling, an essential step is to identify the constitutive parameters and estimate their ranges based on appropriate experiments such as the tensile test of flat segments or the extension-inflation test of cylindrical segments. The mechanical properties of the blood vessel are generally determined by using a specific constitutive model. Blood vessels exhibit hyperelastic and anisotropic properties. It is well known that elastin dominates the initial linear behavior in the low-stiffness regime, whereas collagen is recruited in the high-stiffness regime (Harkness et al., 1957). To identify the anisotropic behavior of blood vessels, pairs of experimental datasets in two loading directions for the same sample are required. Biaxial testing is more reliable in determining the material parameters in the constitutive models, while uniaxial testing is advantageous in investigating the regional variation in the heterogeneous properties of the blood vessel. We previously demonstrated that the posterior side of the porcine thoracic aorta is significantly stiffer than the anterior side in the extension-inflation test, and that the longitudinal difference in the aortic mechanical properties is dependent on the circumferential region (Kim and Baek, 2011; Kim et al., 2013). Therefore, the uniaxial tensile test must be revisited in terms of the spatial description of arterial heterogeneity and its potential application to in silico studies of the virtual aorta.
To elucidate the regional variation of blood vessels, this study was conducted to identify a new approach to integrate a constitutive model and uniaxial test data. To this end, uniaxial tensile tests were performed on the eight sets of aortic segments in consideration of the test regions and loading directions. For an integrated prediction tool based on experimental data, we propose a three-step approach to describe the regional and nonlinear anisotropic mechanical behavior, as well as the relationship between the physical characteristics (arterial stiffness) and constitutive model parameters. First, each virtual two-dimensional (2D) dataset is constructed using a pair of uniaxial tests in the longitudinal and circumferential loading directions for the same sample, followed by the estimation of a set of anisotropic constitutive model parameters. Consequently, a number of virtual datasets are generated via bootstrap aggregation, which can be used to improve the predictive performance of the model. Second, the data populated using the models allows the analysis of their regional variabilities and differences in the mechanical behavior of the aorta through various statistical tools. Finally, the correlations between the arterial physical characteristics and model parameters, as well as the correlations among the model parameters are determined. The proposed approach is evaluated using three constitutive models: the Fung model, Holzapfel model, and constrained mixture model (CMM).
Methods
To investigate the regional variations in the mechanical behavior and material parameters of the aorta, the thoracic aorta was segmented into eight sets according to the circumferential direction (i.e., anterior vs. posterior), longitudinal direction (i.e., proximal vs. distal), and loading directions (i.e., circumferential vs. longitudinal) in the experiment (Figure 1), as follows:
1) circumferential specimen of the anterior side of proximal descending thoracic aorta (PAC)
2) longitudinal specimen of the anterior side of proximal descending thoracic aorta (PAL)
3) circumferential specimen of the anterior side of distal descending thoracic aorta (DAC)
4) longitudinal specimen of the anterior side of distal descending thoracic aorta (DAL)
5) circumferential specimen of the posterior side of proximal descending thoracic aorta (PPC)
6) longitudinal specimen of the posterior side of proximal descending thoracic aorta (PPL)
7) circumferential specimen of the posterior side of distal descending thoracic aorta (DPC)
8) longitudinal specimen of the posterior side of distal descending thoracic aorta (DPL).
FIGURE 1. Test regions of the thoracic aorta in (A) the longitudinal direction and (B) the circumferential direction, and (C) pairs of datasets in two loading directions.
Sample and uniaxial tensile test
Eight fresh descending thoracic aortas from approximately 6-month-old pigs were obtained from a local slaughterhouse. They were immediately used after their loose connective tissues were removed from the adventitia. The aorta was cut in the order of circumferential and longitudinal directions with a width of 6 mm along the aortic tree. The total number of aortic segments used in this study was 35, 29, 26, 28, 25, 25, 28, and 28 for the PAC, PAL, DAC, DAL, PPC, PPL, DPC, and DPL, respectively. The thickness of the sample was measured twice at three different points using a Vernier caliper. The ends of the sample were attached to sandpaper with cyanoacrylate to avoid slipping during the experiments. Two microspheres with a 10 mm gap were attached to the center of the sample. Subsequently, the sample was mounted on the experimental device at a length of 20 mm between the clamps of the device, which was originally developed for the extension-inflation test (Kim and Baek, 2011) but was modified for the uniaxial tensile test. Before commencing the tensile test, the distance between two clamps was adjusted to flatten the sample in the absence of bending, which was confirmed by the microspheres attached to the sample. The sample was maintained in the moist state and then stretched at a speed of 2.9 mm/s until rupture. During the test, the axial force exerted on the sample was measured continuously at a sampling rate of 100 Hz with a 50 N load cell. The images of the sample were also taken using a CCD camera, and the position of the microspheres were tracked through the post imaging process with custom-written Matlab codes and were used to determine the stretch.
Calculation of slope in stress–stretch curve
The experimental data up to fiber breakage points in which the stress–stretch curve deviated from the high-stiffness regime was used for analysis. Two distinct regimes were observed in the stress–stretch curve: A moderate gradient slope (low-stiffness) regime and a steep gradient slope (high-stiffness) regime. To quantitatively characterize the nonlinear behavior of the aortic tissues, the slope of the tangent to the stress–stretch curve was determined. The slope was determined in two steps for each regime to represent the stiffness of the aortic tissue. Since the stress-stretch curve can be approximated as a piecewise linear function, the optimal range of each linear fit before and after transition was first determined based on the coefficient of determination (R2). The first point of the stress-stretch curve served as a starting point for the optimal range in the low-stiffness regime, and the last point of the curve was used as an endpoint in the high-stiffness regime. Next, the fitting equation to describe the stress-stretch curve was defined and the derivative of the fitting equation was determined at one point within the optimal range for each regime.
Constitutive models for hyperelasticity
Stress tensor and incompressibility
A majority of constitutive models consider the elastic behavior of blood vessels to be hyperelastic. Aortic tissues are typically characterized by the hyperelastic strain energy function. The aortic wall is assumed to be an incompressible material, i.e., J2 = det C = 1, where J is the Jacobian, and C is the right Cauchy–Green deformation tensor. The Cauchy stress tensor, T, which describes the true stress state in the material configuration, can be defined by the hyperelastic strain energy as follows:
where p is the Lagrange multiplier that can be determined by the boundary conditions, F the deformation gradient, and W the strain energy function. For the uniaxial test, the deformation gradient of the incompressible material is expressed as
assuming that stress is applied along the loading direction (Martins et al., 2006; Annaidh et al., 2012; de Rooij and Kuhl, 2016), where λ is the stretch ratio.
Fung model
Chuong and Fung (Chuong and Fung, 1983; Zulliger et al., 2004) proposed a strain energy function to describe the mechanical behavior of an artery under internal pressure or stretching. They considered the anisotropic property and nonlinearity of the artery, and the mechanical behavior is described through the strain energy function W.
where
where E is the Green–Lagrange strain tensor, c1 is an elastic constant, and b1–b6 are parameters describing the contribution of the principal strains. The uniaxial stress–stretch relationship for the Fung model, based on Eqs 3, 5, is expressed as (by solving for the Lagrange multiplier p using the free boundary condition in each uniaxial loading direction):
where
where
Holzapfel model
Holzapfel et al. (Holzapfel et al., 2000) proposed a strain energy function that includes two families of collagen fibers. They assumed that each layer of the arterial wall demonstrated similar mechanical characteristics (Schroeder et al., 2018). In this model, the strain energy function is divided into two: Wiso associated with isotropic deformations and Waniso associated with anisotropic deformations. The neo-Hookean model is used in the isotropic response component, as follows (Ogden, 1997):
where c > 0 is a stress-like material parameter. The strain energy function associated with anisotropic deformation is proposed to be an exponential function, as follows:
where k1 > 0 is a stress-like material parameter, and k2 > 0 is a dimensionless parameter. The fiber direction can be described as a model parameter using I4 and I6. The uniaxial stress–stretch relationship in each loading direction for the Holzapfel model based on Eq. 12 and Eq. 13 is expressed as follows:
where
where
In Eqs 14–17,
CMM
The CMM proposed by Humphrey and Rajagopal (2002) presents general concepts for the modeling of biological tissues. The basic principle of the model is that living tissues have a preferred state. When the preferred state is changed, the rates of production and removal and the natural configurations of multiple constituents can be changed to restore them to the preferred state. Typically, the CMM is utilized in the computational simulation of vascular growth and remodeling; however, in this study, the CMM is presented as a constitutive model that describes the arterial mechanical behavior. In addition, in the CMM, the tissue is considered to be a homogenized mixture and the deformation of the mixture as a whole is constrained by incompressibility. In the model, the stored energy functions are given by
where
for elastin and collagen fiber families, respectively. Here,
where
where
where
The uniaxial stress–stretch relationship in each loading direction is expressed as
where
Calculation of material parameters
The biaxial tensile or inflation test can be used to investigate the anisotropic behavior of the aorta; however, it is limited to the characterization of spatially heterogeneous material parameters in local regions. The uniaxial tensile test allows the intrinsic material characteristics in the local regions of the blood vessel to be estimated. However, it is difficult to estimate all parameters simultaneously, particularly when numerous parameters are involved. Thus, a novel data-oriented approach, i.e., the bootstrap approach, was employed in this study. More specifically, cross-mapping involving one circumferential-to-many longitudinal mapping (or vice versa) was performed to generate a virtual biaxial dataset using custom-written Matlab codes. For a virtual dataset, a set of material parameters in the region of interest was estimated individually for each loading direction. They were then integrated to determine the common range of the material parameters for the biaxial behavior in a given region. However, a single set of material parameters for the biaxial behavior may result in a large error, if the ranges of the material parameters for each loading direction are too wide or the common ranges are too narrow such as in transition. To minimize the error, the error tolerance for the iterative solution in the Matlab codes was adjusted based on the goodness of the curve fit.
Statistical analyses
In this study, two separate statistical analyses were performed. Analysis of variance was used to identify regional differences in the experimental results, followed by the Bonferroni correction as a post-hoc test. Additionally, the Student’s t-test was performed to examine the significance of the regional differences as well as the relationship between the material parameters and experimental results.
To determine the appropriate method for performing the Student’s t-test and correlation study, the following procedures were performed: 1) The normal distribution of the calculated material parameters and stiffness as well as the experimental data were examined using the Shapiro–Wilk test, as shown in Figure 2. 2) In the normality test, a significant difference was determined based on the p-value of 0.05. 3) The Student’s t-test was conducted to examine the differences between the material parameters based on the test regions. 4) If the normality test of the data fails, then the data can be transformed to the normally distributed data via Box–Cox transformations with an optimal Lambda using MATLAB codes, followed by the Student’s t-test. In addition, Pearson and Spearman rank correlation studies were conducted between the estimated material parameters and experimental data, as well as a correlation study among material parameters based on the test regions. Once the correlation coefficients (r) were calculated, the relationship strengthened in the positive and negative directions as the correlation coefficients approached +1 and −1, respectively. In addition, no correlation was observed when the correlation coefficient was 0 or the p-value was 0.05 or more. The SPSS software was used for the statistical analysis.
Results
Stress–stretch response
The overall stress–stretch responses were nonlinear with pronounced regional differences (Figures 3, 4). In terms of the anisotropy of the blood vessel, the circumferential specimens were stiffer than the longitudinal specimens. In addition, the transition of the circumferential specimens from the low- to high-stiffness regimes was greater than that of the longitudinal specimens.
FIGURE 3. Stress–stretch response for each test region. Left column shows circumferential specimens [(A) PAC, (C) PPC, (E) DAC, and (G) DPC], and right column shows longitudinal specimens [(B) PAL, (D) PPL, (F) DAL, and (H) DPL].
FIGURE 4. Mean and standard error of mean for slopes of stress–stretch curves of aortic tissues in (A) low-stiffness regimes and (B) high-stiffness regimes for circumferential and longitudinal specimens, respectively.
The slopes of the stress-stretch curves for aortic tissues represent the stiffness of the aortas, which depend on the test region, the loading direction, and the regime of nonlinear behavior. The slope of the distal posterior region (DPC and DPL) was the greatest in both regimes for each loading direction. Significant difference between the anterior and posterior sides was found in the high-stiffness regime (DPC vs. DAC, and PPL vs. PAL). Significant difference between the proximal and distal regions was found at the posterior side, but only in the longitudinal specimens of both regimes (PPL vs. DPL). In addition, the interaction between the circumferential and longitudinal regions was significant for the circumferential specimens (PAC vs. DPC) in both regimes and for the longitudinal specimens (PPL vs. DPL) in the high-stiffness regime.
Material parameters
The representative fitting results for each constitutive model along with the experimental data are shown in Figure 5. The results indicate that the fitting method used in this study is reliable.
FIGURE 5. Comparison between experimental data and constitutive models based on regional aortic specimens. (A) DAC–DAL, (B) DPC–DPL, (C) PAC–PAL, and (D) PPC–PPL. C: circumferential specimen; L: longitudinal specimen.
Figures 6–8 show the Box–Whisker plots for the distribution of the estimated material parameters for three constitutive models based on the experimental data. Material parameters were calculated by sequentially applying pairs of experimental datasets for the circumferential and longitudinal specimens and were estimated without constraints. Subsequently, curve fitting was performed for 9,456 cases [{PAC–PAL (35 × 29 = 1,015 cases), PPC–PPL (25 × 25 = 625 cases), DAC–DAL (26 × 28 = 728 cases), DPC–DPL (28 × 28 = 784 cases)} × three constitutive models]. The minimum, maximum, first quartile (1Q), third quartile (3Q), and median of the estimated material parameters for all three models are listed in the Supplementary Appendix S1.
FIGURE 6. Random value distribution of material parameters estimated using Fung model for PA, PP, DA, and DP specimens. (A) c1, (B) b1, (C) b2, (D) b3, (E) b4, (F) b5, and (G) b6.
FIGURE 7. Random value distribution of material parameters estimated using Holzapfel model for PA, PP, DA, and DP specimens. (A) c, (B) k1, (C) k2, and (D) γ.
FIGURE 8. Random value distribution of material parameters estimated using CMM for PA, PP, DA, and DP specimens. (A)
The values of the material parameters, c1, b1, b2, b3, b4, b5, and b6, were obtained using the Fung model, as shown in Figure 6. The material parameters of c1 on the DA specimens indicated the widest interquartile ranges (IQRs) and the highest median value of 0.3012. The values of b1 and b2 for the DP specimens were the most widely distributed in the IQR, with the highest median values of 2.4510 and 1.1120, respectively. Additionally, material parameters b3 and b5 indicated a wide IQR, with slight differences in the median over all test regions. The means ± standard deviations of all median values for parameters b3 and b5 were 1.2538 ± 0.1706 and 0.9191 ± 0.1705, respectively. The median values of b4 and b6 were relatively low compared with those of the other parameters. In particular, the material parameters of b4 on the DP specimens and b6 in the PP region were approximately zero.
The values of material parameters c, k1, k2, and γ obtained using the Holzapfel model are shown in Figure 7. In the regional comparison, the material parameter c of the PP specimens exhibited the highest median value of 0.0581. Meanwhile, the material parameter k2 of the DP specimens indicated the widest IQR, with the highest median value of 3.3020. Material parameters k1 and γ indicated only slight differences over all the test regions.
The values of material parameters
Significant difference of the material parameters
The Student’s t-test for the estimated material parameters was performed to determine the difference among all test regions. It was discovered that material parameters c1 and b6 in the Fung model, k2 and γ in the Holzapfel model, and
Statistical correlations
The estimated material parameters were used sequentially for the correlation study. Spearman rank correlation coefficients (r) were computed to quantify the correlation degree between the regional experimental data (i.e., slopes of stress-stretch curves) and material parameters or among the material parameters. Only cases exhibiting a correlation in all test regions are reported, as shown in Figures 9–11. To understand the general trend from the correlation study, the average values of the correlation coefficients over the four test regions (PA, PP, DA, and DP) were obtained.
FIGURE 9. Spearman rank correlation coefficients (A) between material parameters of Fung model and experimental results and (B) among the material parameters.
FIGURE 10. Spearman rank correlation coefficients (A) between material parameters of Holzapfel model and experimental results and (B) among the material parameters.
FIGURE 11. Spearman rank correlation coefficients (A) between material parameters of CMM and experimental results and (B) among the material parameters.
In the correlation study between the parameters of Fung model and experimental results, b1 (r = 0.665 on average) indicated the strongest correlation with the slope in the high-stiffness regime, hereinafter known as the high slope, in the circumferential specimens, followed by b5 (0.408) (Figure 9A). Meanwhile, b2 showed the strong correlation with both high slope (0.617) and low slope (0.310) in the longitudinal specimens.
In the correlation study among the material parameters of the Fung model, c1 indicated the strongest correlation with b1 (−0.578), followed by b2 (−0.492) (Figure 9B). Next, b1 and b2 indicated a strong correlation with b5 (0.458) and b6 (0.364), respectively.
For the correlation study between the parameters of Holzapfel model and experimental results, γ indicated the strongest correlation with the high slope (0.583) and low slope (0.380) in the longitudinal specimens (Figure 10A). Meanwhile, k2 (0.485) showed the strong correlation with the high slope in the circumferential specimens.
In the correlation study among the material parameters of the Holzapfel model, c indicated the strongest correlation with k1 (−0.759), followed by k2 (0.459) (Figure 10B). In addition, k1 was negatively correlated with k2 (-0.539).
In the correlation study between the parameters of CMM and experimental results,
In the correlation study among the material parameters of the CMM,
Discussion
In this study, first of all, a virtual 2D regional dataset is generated by experimental data from uniaxial tensile tests to overcome the limitations of conventional experimental setups to characterize the anisotropic behavior of the artery. We use cross-mapping to determine the material parameters by applying experimental data in both loading directions. The relationship between the physical features (arterial stiffness) and constitutive parameters is determined in different regions. Subsequently, we present the ranges of material parameters and propose an integrative method to characterize the region-dependent material parameters for the nonlinear anisotropic materials.
Kim and Baek (Kim and Baek, 2011) and Kim et al. (Kim et al., 2013) report the spatial variations in the stiffness of the aorta for both circumferential and longitudinal regions via the extension-inflation test. In this study, uniaxial tensile tests present consistent significant differences. The circumferential specimen is stiffer than the longitudinal specimen in both the low- and high-stiffness regimes, and the circumferential and longitudinal differences are found among several regions. In particular, the DP exhibits the stiffest slope for each of the circumferential and longitudinal loading directions in both low- and high-stiffness regimes. This is likely due to the non-uniform fiber orientations on the different test regions of the aorta. Various histological and mechanical analyses indicate diverse fiber orientations of the intimal, medial, and adventitial strips (Huang et al., 2016; Niestrawska et al., 2016; Yu et al., 2018; Babici et al., 2020; Concannon et al., 2020; Jadidi et al., 2020; Díaz et al., 2021). Because the DP is the thinnest region of the aorta, the transmural variation in the fiber orientation and the contribution of fibers to the stress–stretch relationship in the DP may differ from those of the other test regions.
Although the material behavior of the aorta is described by a constitutive model, it is not uniquely determined for the best representation. Complex anisotropic materials, such as biological tissues, typically have a single isotropic plane, so that they can be idealized as transversely isotropic in hyperelastic models (O’Shea et al., 2019). In this study, the aorta is assumed to be a transversely isotropic material to simulate the uniaxial tensile test. Although this approach has limitations in ensuring the positive definiteness of the fourth-order stiffness tensor, it provides a stable solution to elasticity problems if a proposed model is calibrated based on experimental data (O’Shea et al., 2019). In addition, in the uniaxial test, the transversely isotropic model to consider only I4 invariant shows responses similar to that of the method considering both I4 and I5 invariants (Feng et al., 2016). In shear deformation, however, both should be considered, as the results depend on whether I5 invariant is considered.
In this study, three constitutive models, i.e., the Fung model, Holzapfel model, and CMM, are considered with this methodology to describe the material responses in the uniaxial test. All three models present good fits with the stress–stretch curves obtained from the uniaxial test. The Fung model has the advantage of fitting curves with a large variation, whereas the Holzapfel model and CMM have the significant advantage of considering vascular constituents such as elastin and collagen. In the CMM, the highest number of material parameters is used, including the pre-stretch value in the homeostatic state. Nevertheless, it presents a substantial advantage in that it can accurately describe nonlinear behavior of the aorta, particularly in the high-stiffness regime. To improve the computational efficiency of the CMM composed of many parameters, the estimation method is classified into two steps in a previous study (Seyedsalehi et al., 2015). However, unlike the previous study, the method used in this study involves integrating material parameters for the two loading directions into the virtual biaxial behavior. No significant difference in the computational efficiency is observed between the two methods regardless of integration.
This study shows that the material parameters of each constitutive model can express the mechanical behavior of the regional aortic tissue. There are significant differences in the material parameters with respect to the regions. For example, the material parameter γ in the Holzapfel model and the material parameter
In addition to regional differences, strong correlations are presented between regional aortic behavior and material parameters or among material parameters for the same test region. As shown in Figures 9–11, a higher absolute value of the correlation coefficient indicates a stronger correlation. In the Fung model, the material parameters c1, b1, b2, b3, b5, and b6 correlate with the slope of the high-stiffness regime. Although the material parameters of the Fung model provide a description of the stress–stretch relationship phenomenologically, they are distinct with respect to the regions and regimes when describing the mechanical behavior with the physical characteristics and slopes. In particular, the material parameters b1 and b2 have a direct effect on the simulating the stiffness of the materials in each loading direction, because they are multiplied by the stretch in each loading direction. The ratio of b1/b2 can be used to characterize the anisotropic behavior of the materials. In this study, this ratio may give a better way to represent the behavior of materials in the correlation study between the ratio and stiffness regime. It ranges from −0.473 to −0.171 among the slopes in the low-stiffness regime for longitudinal specimens, and from 0.298 to 0.503/from −0.645 to −0.473 among the slopes in the high-stiffness regime for circumferential/longitudinal specimens, respectively. In the Holzapfel model, all parameters correlate with the slope in the high-stiffness regime, while the parameters k1, k2, and γ correlate with the slope in the low-stiffness regime for the longitudinal specimens. In addition, the parameter c is inter-correlated with all other parameters, and particularly has the strongest negative correlation with k1. Finally, in the CMM, the parameters
There are limitations of this study. For the experimental protocol, the deformation of aortic segments is assumed to be uniform, whereas the sample might not deform uniformly since the blood vessel is a composite material. However, the specimens are small, and their spatial variations are generally neglected in the uniaxial tensile test. Especially, the stretch is defined by the change in distance between two microspheres because it gives an accurate estimation of local deformation at the central region, instead of measuring the grip movement in the uniaxial test. Nonetheless, the specimens are not necessarily ruptured in the central region between the two microspheres so the uniformity assumption cannot apply to the rupture behavior. To further investigate the nonuniform deformation, adding more markers to the sample or utilizing the digital image correlation technique may be an alternative.
The statistical analysis performed in this study confirms that the statistical significance of key parameters is important in differentiating the regional mechanical behavior of the aorta. However, the constitutive models used in this study are nonlinear and, as shown in the correlation studies, the physical characteristics (e.g., stiffness in high-stiffness regime) are highly correlated with the material parameters of each model. Analyzing such multidimensional data (regions, loading directions, low- and high-stiffness regimes, constitutive models and material parameters, and enhanced experimental measurements) requires further studies using a machine learning approach to predict the regional mechanical behavior of the blood vessel for the virtual aorta.
Due to the complexity between the clinical features and model parameters, the number of model parameters is reduced to estimate the physical features within the model. A decrease in the number of independent parameters may be more advantageous, based on overall nonlinear sensitivity analysis and model uncertainty quantification. For example, a multidimensional virtual dataset can quantify model uncertainties through correlation studies. A linear regression model with Bayesian inference approaches can be used to describe significant material parameters using backward elimination (Akkoyun et al., 2020). Alternatively, the Gaussian mixture model is based on an unsupervised machine learning approach involving an uncertainty measure (probability) (Paalanen et al., 2006), which determines how much the material parameters are associated with a specific region by considering the mean and covariance of each cluster and mixing probability. Therefore, the combination of material parameters can improve the prediction capability of the model with respect to the different regions and the selection of a constitutive model, appropriately characterizing the regional mechanical behavior of the blood vessel based on a machine learning approach.
In summary, this study presents the region-dependent mechanical characterization of the porcine thoracic aorta by generating virtual 2D datasets from the uniaxial test using three constitutive models. Physics-based modeling is essential to associate with physical and clinical features. The development of computational models can further promote the era of “virtual human” or “virtual aorta.” However, the development of specific models is hampered by the scarcity of material test data for individual patients, which are relevant to clinical setups. In particular, tissue mechanical tests of young human aortic specimens are rarely available, even if they exist. The virtual biaxial data presented herein may serve as an alternative to material behavior in young human aortas.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
Ethical review and approval were not required in this study because the porcine aortas used in this study were provided by a local slaughterhouse.
Author contributions
DR, SB, and JK contributed to conception and design of the study. JK performed the experiments and data analysis. DR performed computation and statistical analysis. SB developed the theoretical and computational framework and verified the analytical methods. JK and DR wrote the first draft of the manuscript. All authors discussed and contributed to manuscript revision and approved the submitted version.
Funding
This study was supported partially by the National Institute of Health (R01HL115185), the National Research Foundation of Korea (NRF) grant funded by the Korea government (No. 2021R1C1C2004703), and the Basic Science Research Program through the NRF funded by the Ministry of Education (No. 2021R1I1A3059740).
Acknowledgments
The authors acknowledge that Dr. Emrah Akkoyun provided data-driven learning expertise and assisted in the editing of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2022.937326/full#supplementary-material
References
Akkoyun, E., Kwon, S. T., Acar, A. C., Lee, W., and Baek, S. (2020). Predicting abdominal aortic aneurysm growth using patient-oriented growth models with two-step Bayesian inference. Comput. Biol. Med. 117, 103620. doi:10.1016/j.compbiomed.2020.103620
Annaidh, A. N., Bruyère, K., Destrade, M., Gilchrist, M. D., Maurini, C., Otténio, M., et al. (2012). Automated estimation of collagen fibre dispersion in the dermis and its contribution to the anisotropic behaviour of skin. Ann. Biomed. Eng. 40 (8), 1666–1678. doi:10.1007/s10439-012-0542-3
Babici, D., Kudej, R. K., McNulty, T., Zhang, J., Oydanich, M., Berkman, T., et al. (2020). Mechanisms of increased vascular stiffness down the aortic tree in aging, premenopausal female monkeys. Am. J. Physiology-Heart Circulatory Physiology 319 (1), H222–H234. doi:10.1152/ajpheart.00153.2020
Chakshu, N. K., Sazonov, I., and Nithiarasu, P. (2021). Towards enabling a cardiovascular digital twin for human systemic circulation using inverse analysis. Biomech. Model. Mechanobiol. 20 (2), 449–465. doi:10.1007/s10237-020-01393-6
Chuong, C. J., and Fung, Y. C. (1983). Three-dimensional stress distribution in arteries. J. Biomech. Eng. 105 (3), 268–274. doi:10.1115/1.3138417
Concannon, J., Dockery, P., Black, A., Sultan, S., Hynes, N., McHugh, P. E., et al. (2020). Quantification of the regional bioarchitecture in the human aorta. J. Anat. 236 (1), 142–155. doi:10.1111/joa.13076
de Rooij, R., and Kuhl, E. (2016). Constitutive modeling of brain tissue: Current perspectives. Appl. Mech. Rev. 68 (1), 010801. doi:10.1115/1.4032436
Díaz, C., Peña, J. A., Martínez, M. A., and Peña, E. (2021). Unraveling the multilayer mechanical response of aorta using layer-specific residual stresses and experimental properties. J. Mech. Behav. Biomed. Mat. 113, 104070. doi:10.1016/j.jmbbm.2020.104070
Feng, Y., Okamoto, R. J., Genin, G. M., and Bayly, P. V. (2016). On the accuracy and fitting of transversely isotropic material models. J. Mech. Behav. Biomed. Mat. 61, 554–566. doi:10.1016/j.jmbbm.2016.04.024
Harkness, M. L. R., Harkness, R. D., and McDonald, D. A. (1957). The collagen and elastin content of the arterial wall in the dog. Proc. R. Soc. Lond. B Biol. Sci. 146 (925), 541–551. doi:10.1098/rspb.1957.0029
Holzapfel, G. A., Gasser, T. C., and Ogden, R. W. (2000). A new constitutive framework for arterial wall mechanics and a comparative study of material models. J. Elast. 61, 1–48. doi:10.1023/A:1010835316564
Hose, D. R., Lawford, P. V., Huberts, W., Hellevik, L. R., Omholt, S. W., and van de Vosse, F. N. (2019). Cardiovascular models for personalised medicine: Where now and where next? Med. Eng. Phys. 72, 38–48. doi:10.1016/j.medengphy.2019.08.007
Huang, A. H., Balestrini, J. L., Udelsman, B. V., Zhou, K. C., Zhao, L., Ferruzzi, J., et al. (2016). Biaxial stretch improves elastic fiber maturation, collagen arrangement and mechanical properties in engineered arteries. Tissue Eng. Part C. Methods 22 (6), 524–533. doi:10.1089/ten.TEC.2015.0309
Humphrey, J. D., and Rajagopal, K. R. (2002). A constrained mixture model for growth and remodeling of soft tissues. Math. Models Methods Appl. Sci. 12 (3), 407–430. doi:10.1142/S0218202502001714
Jadidi, M., Habibnezhad, M., Anttila, E., Maleckis, K., Desyatova, A., Mactaggart, J., et al. (2020). Mechanical and structural changes in human thoracic aortas with age. Acta Biomater. 103, 172–188. doi:10.1016/j.actbio.2019.12.024
Kim, J., and Baek, S. (2011). Circumferential variations of mechanical behavior of the porcine thoracic aorta during the inflation test. J. Biomech. 44 (10), 1941–1947. doi:10.1016/j.jbiomech.2011.04.022
Kim, J., Hong, J.-W., and Baek, S. (2013). Longitudinal differences in the mechanical properties of the thoracic aorta depend on circumferential regions. J. Biomed. Mater. Res. Part A 101A, 1525–1529. doi:10.1002/jbm.a.34445
Martins, P. A. L. S., Natal Jorge, R. M. N., and Ferreira, A. J. M. (2006). A comparative study of several material models for prediction of hyperelastic properties: Application to silicone-rubber and soft tissues. Strain 42 (3), 135–147. doi:10.1111/j.1475-1305.2006.00257.x
Niederer, S. A., Aboelkassem, Y., Cantwell, C. D., Corrado, C., Coveney, S., Cherry, E. M., et al. (2020). Creation and application of virtual patient cohorts of heart models. Phil. Trans. R. Soc. A 378, 20190558. doi:10.1098/rsta.2019.0558
Niestrawska, J. A., Viertler, C., Regitnig, P., Cohnert, T. U., Sommer, G., and Holzapfel, G. A. (2016). Microstructure and mechanics of healthy and aneurysmatic abdominal aortas: Experimental analysis and modelling. J. R. Soc. Interface 13 (124), 20160620. doi:10.1098/rsif.2016.0620
O’Shea, D. J., Attard, M. M., and Kellermann, D. C. (2019). Hyperelastic constitutive modelling for transversely isotropic composites and orthotropic biological tissues. Int. J. Solids Struct. 169, 1–20. doi:10.1016/j.ijsolstr.2018.07.013
Paalanen, P., Kamarainen, J. K., Ilonen, J., and Kälviäinen, H. (2006). Feature representation and discrimination based on Gaussian mixture model probability densities—practices and algorithms. Pattern Recognit. DAGM. 39 (7), 1346–1358. doi:10.1016/j.patcog.2006.01.005
Schroeder, F., Polzer, S., Slažanský, M., Man, V., and Skácel, P. (2018). Predictive capabilities of various constitutive models for arterial tissue. J. Mech. Behav. Biomed. Mat. 78, 369–380. doi:10.1016/j.jmbbm.2017.11.035
Seyedsalehi, S., Zhang, L., Choi, J., and Baek, S. (2015). Prior distributions of material parameters for Bayesian calibration of growth and remodeling computational model of abdominal aortic wall. J. Biomech. Eng. 137 (10), 101001. doi:10.1115/1.4031116
Yu, X., Wang, Y., and Zhang, Y. (2018). Transmural variation in elastin fiber orientation distribution in the arterial wall. J. Mech. Behav. Biomed. Mat. 77, 745–753. doi:10.1016/j.jmbbm.2017.08.002
Zeinali-Davarani, S., Sheidaei, A., and Baek, S. (2011). A finite element model of stress-mediated vascular adaptation: Application to abdominal aortic aneurysms. Comput. Methods Biomechanics Biomed. Eng. 9, 803–817. doi:10.1080/10255842.2010.495344
Keywords: virtual dataset, regional variation, arterial stiffness, Fung model, holzapfel model, constrained mixture model, statistical correlations, bootstrapping
Citation: Ryu D, Baek S and Kim J (2022) Region-dependent mechanical characterization of porcine thoracic aorta with a one-to-many correspondence method to create virtual datasets using uniaxial tensile tests. Front. Bioeng. Biotechnol. 10:937326. doi: 10.3389/fbioe.2022.937326
Received: 06 May 2022; Accepted: 29 August 2022;
Published: 11 October 2022.
Edited by:
Emmanuel A. Audenaert, Ghent University, BelgiumReviewed by:
Liang Wang, Southeast University, ChinaHiroshi Yamada, Kyushu Institute of Technology, Japan
Copyright © 2022 Ryu, Baek and Kim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jungsil Kim, jungsilkim@scnu.ac.kr