Strojniški vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 © 2015 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2014.2061 Original Scientific Paper Received for review: 2014-07-15 Received revised form: 2014-10-28 Accepted for publication: 2014-11-03 Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State Janez Urevc1* - Milan Brumen2. 3 - Vojko Flis4 - Boris Štok1 iUniversity of Ljubljana, Faculty of Mechanical Engineering, Slovenia 2University of Maribor, Faculty of Medicine, Slovenia 3 Jožef Stefan Institute, Slovenia 4University Medical Centre Maribor, Department of Vascular Surgery, Slovenia The prediction of the arterial zero-stress state is as of yet an unresolved problem in the field of modelling the mechanical response of patient-specific arteries. This is because the configuration associated with arterial zero-stress state is impossible to obtain experimentally. However, the zero-stress configuration (or, equivalently, the residual stresses related to this configuration) represents the crucial data of any numerical analysis. In this study, the mechanical response of a pre-stressed, pressurised, hyperelastic tube (representing the artery) is determined without knowing the initial zero-stress configuration of the artery. Instead, to predict the arterial residual (bending and stretching) stresses, a corresponding thermomechanical analogy is used. As shown in the paper, the arterial residual stress state can equally be obtained by thermally loading a properly defined closed tube. Thus, based on the loaded state of the corresponding thermomechanical model, the arterial residual stress sate is constructed, from where the arterial loaded state can be obtained. The initial configuration of the thermomechanical model is defined on the basis of the arterial loaded configuration. The methodology is validated by predicting the zero-stress state of the artery. The predicted equilibrium state of the artery, when cut longitudinally and transversally, has the form of an opened-up tube with a relatively low stress state in comparison to the arterial residual stresses. The results thus demonstrate that arterial residual stresses that are predicted with the corresponding thermomechanical model exhibited the bending distribution, which proves the methodology to be adequate. Keywords: residual stresses, zero-stress state of arteries, finite element method, thermomechanics Highlights • Predicting the arterial mechanical response based on the in vivo arterial data. • Proposed methodology for predicting the arterial mechanical response. • Using thermomechanical analogy for predicting the arterial residual stress state. • Predicting the arterial zero-stress state based on the in vivo arterial data. • Determination of bending stresses by means of thermomechanical analogy. 0 INTRODUCTION The progress of medical imaging has enabled the increasingly accurate acquisition of patient-specific arterial geometries. This, in turn, allows for the study of arterial response under various conditions, giving additional insight into human physiology that cannot be experimentally obtained. Medical imaging is most often performed in vivo, thus providing the arterial geometry that corresponds to the arterial loaded (i.e., physiological) state. To numerically model the arterial mechanical response, however, either the arterial zero-stress configuration or the arterial residual stress state is needed. This study addresses the problem of modelling the arterial response in cases when the arterial zero-stress configuration is not given whereas the arterial loaded configuration is known. The artery is treated as a hyperelastic, residually pre-stressed and pressurised tube. The arterial unloaded state represents a part of the artery removed from the surrounding tissue, freed of intraluminal pressure, and removed from the body. In this state, the artery is free of external loads but is still subject to some amount of stresses, referred to as the residual stresses. These stresses are relieved by cutting the artery longitudinally. The artery then springs open and assumes a sector-like shape in which the size of this opening quantifies the magnitude of the relieved residual stresses. Because of their transmural distribution, with compressive stresses on the arterial inside and tensile stresses on the outside [1], the residual stresses are normally associated with pure bending stresses [2]. Although some residual stresses exist in the arterial longitudinal direction [3] and [4], we will assume that only one longitudinal cut is sufficient to release almost all residual stresses [5] and [6]. In this work, the artery is treated as single-layered structure. Even though this is a rough approximation of the complex structure of an artery, it is able to provide a reasonable assessment of the arterial physiological response [7]. *Corr. Author's Address: University of Ljubljana, Faculty of Mechanical Engineering, Aškerčeva 6, 1000 Ljubljana, Slovenia, janez.urevc@fs.uni-lj.si 5 Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 Residual stress produces a favourable (physiological) environment for cell functions [8], which is created by lowering the gradient of the transmural stresses [9] and [10]. According to the so-called uniform strain hypothesis [11], the distribution and magnitude of the residual stresses cause the physiological intraluminal distribution of the circumferential strain X^(r) to be uniform (X^(r) = const), which implies a low gradient of transmural stresses. Furthermore, it has been experimentally observed that the distribution of the residual stresses in the vessel wall is three dimensional, meaning that they depend on both the radial position and on the axial position [5], [12] and [13]. A crucial datum needed for any numerical prediction of the arterial response, as already mentioned, is the zero-stress configuration of the artery. To the knowledge of the authors, the only methodology that allows the inclusion of residual stresses in a patient-specific arterial geometry is presented in the study performed by Alastrué et al. [2]. There, the residual stresses are accounted for by applying a simplified initial strain field that is obtained from a computation of the residual stress field in a circular arterial segment. The drawback of this approach is that the zero-stress configuration of the circular segment, which must be adopted from the literature, does not correspond with a patient-specific artery. Furthermore, due to the complexity of the problem, various simplifications have been made in the literature in order to obtain the patient-specific arterial zero-stress configuration. The most general approach is to take the loaded configuration of the artery for the arterial initial stress-free state. This approach leads to unphysically large arterial deformations [14] where, subsequently, the predicted arterial configuration obtained by pressurising the artery with intraluminal pressure strongly differs from the one observed in vivo. Another approach is to neglect the arterial residual stresses and assume the arterial unloaded configuration for the arterial initial zero-stress configuration. For example, in the study conducted by Gee et al. [14], two such methods are presented that enable a prediction of the pre-stress state of a patient-specific geometries under in vivo loads. The first is the inverse design (ID) analysis. In this approach, which is also used in [15] and [16], the zero-stress configuration is computed and then used to compute a (non-unique) strain and stress state that satisfied equilibrium with respect to the known spatial configuration. The second method is based on a modified updated Lagrangian formulation (MULF), in which a multiplicative split of the deformation gradient is used to produce a displacement-free pre-strained/pre-stressed state. Despite the widespread use of the first method in particular, the drawback of such an approach is that it neglects the residual stresses. On the other hand, progress has been made in obtaining the arterial unloaded state starting from the opened-up configuration (see for instance [12], [17] and [18]). In this case, however, the arterial opened-up configuration needed to be provided (i.e., given). In this paper, a new methodology is presented that predicts the mechanical response of a hyperelastic, residually pre-stressed and pressurised tube (representing an artery). The loaded configuration of the artery is assumed known whereas the arterial opened-up zero-stress configuration is unknown. The methodology for predicting the arterial residual (i.e., bending and stretching) stresses is based on an equivalent approach where a corresponding thermomechanical problem of an associated closed tube is defined. In Section 1, a case study is presented and the assumptions of the considered mechanical problem are specified. Section 2 gives a theoretical background on (i) the characterisation of the circular arterial segment, (ii) the thermomechanical analogy for obtaining the arterial residual stresses, and (iii) the use of the thermomechanical analogy to predict the mechanical response of a realistic arterial segment. The methodology presented for the arterial segment is then generalised to treat a real artery in Section 3. The analysis results obtained by applying the developed methodology to the case study are given in Section 4 together with a corresponding validation of the methodology. Finally, in Section 5, the conclusions are given along with some guidelines for future research. 1 PROBLEM DEFINITION The methodology for predicting the arterial residual stress state is demonstrated using a case study analysis considering an analytically defined, hyperelastic tube. The tube, here referred to as the arterial model (AM), is defined so as to resemble the in vivo observed artery. The inner, i.e., the luminal surface of the AM (rA , Fig. 1a) is defined as a straight tube of an elliptical cross section having respective principal directions coincident with the x and y axis along the entire AM. The assumed variation of the elliptical semi-minor and semi-major axis (a(z) and b(z)), specifying the surface r M, is shown in Fig. 1b. Based on this variation, the AM luminal cross-sectional area A(z) variation is computed as show in Fig. 1b. Furthermore, the thickness t of the AM is assumed to be constant along the whole section, its 6 Urevc, J. - Brumen, M. - Flis, V. - Stok, B. Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 z [mm] " "" z [mm] Fig. 1. Presentation of the case-study: a) Three-dimensional representation of the arterial model (AM); b) the inner, i.e., the luminal surface of the AM (rAM) represents a straight tube of an elliptical cross section which is defined by the ellipse semi-axes a(z) and b(z) value being t = 0.5 mm. The length of the model is 160 mm. Let us assume that the AM geometry is obtained from medical imaging and that the configuration presented in Fig. 1a corresponds to the AM diastolic state with an intraluminal pressure of P = 70 mmHg. For the material behaviour of the AM, an isotropic hyperelastic material model is assumed, described by a strain energy density function W [17]: cb exp -1 (1) where I1 is the first invariant of the strain tensor and ca and cb are the constitutive parameters. In [17], their values were experimentally determined to be ca = 44.2 kPa and cb =16.7. Remark: The artery is treated as a single-layered tube with isotropic properties. This is an oversimplification because it is well known that arteries are composed of three distinct layers (intima, media and adventitia), which have, in general, anisotropic properties [19]. Nevertheless, as discussed in the introduction, the considered AM is still able to provide a reasonable approximation of the arterial response. In this regard, when taking into consideration that the primary goal of this paper is to give a presentation of the new developed methodology, the adopted model approximation may be considered acceptable and justified for the given purpose. Additionally, the effect of the tissue surrounding the artery, which acts as a pressure load on the arterial outer surface, is neglected in this study. 1.1 Arterial Configurations The arterial configurations are denoted by the symbol Q. The subscript supplementing Q specifies the corresponding arterial loading state whereas the superscript of Q indicates the corresponding model. The superscript "AM" will refer to the arterial model and "C" to the configurations of the circular arterial segments. The following arterial configurations are considered (Fig. 2): Loaded state (Q2): In its loaded state the artery is subject to intraluminal pressure P and longitudinal (axial) as well as circumferential pre-stretch. Stretched state (Q1x): The arterial stretched state is obtained by removing the pressure P from Q2. (Note that the term arterial stretched state is not common in the literature; it is defined here because it offers a good targeting state for the methodology presented in the following sections.) Unloaded state (Qj): The unloaded state is obtained by releasing the axial pre-stretch from Q1X. The associated axial contraction is described by the axial stretch X. Since the in vivo data is insufficient to determine X [20], we have adopted here the data of Delfino et al. [17], where X is set to 10% (X = 1.1). Zero-stress state (Q0): Once the artery in its unloaded state is cut open, it is assumed that the stresses in the artery are close to zero. The obtained state is taken as the arterial zero-stress state. The axial contraction A, that is associated with the longitudinal cut is normally neglected, therefore A = 1 [17]. Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State 11 Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 2 MECHANICAL RESPONSE OF A CIRCULAR ARTERIAL SEGMENT The subsequent subsections begin with some theoretical background on the characterisation of the circular arterial segment, and are followed by a presentation on an analogical approach aimed to obtain the stretched state of a circular arterial segment. The latter is achieved by defining a corresponding thermomechanical problem of a closed circular tube. Finally, the application of this analogy to predict the stretched state of a realistic arterial segment is schematically presented. 2.1 Stress Analysis of a Circular Arterial Segment The mechanical response of an arterial segment (i.e., a slice of the artery) can be characterised by simplifying the segment geometry to a straight circular tube [20] to [23]. The zero-stress state configuration of the open circular segment Q^ (Fig. 2) is thus defined with the radii R, and Ro and the opening angle &0. The unloaded state configuration Q^ is then obtained by bending (closing) the segment from its zero-stress state configuration Q^. By further applying the axial stretch X and the inner pressure P to the unloaded state configuration Qf, the stretched state configuration Qu and the loaded state configuration q22 of the segment are obtained, respectively. The deformation gradient tensor F associated with the loaded state configuration q2 accounts for all three successive motions of material particles, displayed in Fig. 2. A material particle (r, < r < ro) can be represented in terms of the principal stretches as: F = diag dr nr . t —,-,XA dR ©0 R - diagX^,XZ ], (2) where Xr, X^, Xz are, respectively, the principal stretches in the radial, circumferential, and axial directions. The right Cauchy-Green tensor is defined as C = FtF , where FT is the transpose of F. The reference radius R (r ) of the material particle in the zero-stress state configuration Q^ is determined by imposing the incompressibility condition J = det F = = 1 and considering the functional expressions for the principal stretches given in Eq. (2): ( 0 , , ,V/2 R = 2 nXAl 2 2 Ro--—(ro - r &t 0 (3) In addition, the radial stretch is written as: ^=i Wz=ntxA> (4) which also comes from the incompressibility condition. Considering the strain energy density function W the Cauchy stress tensor a is obtained as: dW T a=-pl +& = -pI + 2F——Ft , dC (5) where p is the Langrange multiplier, which enforces the incompressibility constraint det F = 1, and I is the identity tensor. Using Eqs. (1) and (5) the non-zero components of the stress tensor are written as: °kk =~P + °kk =~P + hca (6) • exp ^ [Û +1} +1} - 3 ; k = r,},z. Because the considered problem is a radial one, all the quantities in the above equation depend only on the radial direction. In the absence of the body forces the equilibrium equations reduce to: Fig. 2. Configurations of the circular arterial segment: loaded Q2, stretched Q1X , unloaded Q1 and zero-stress state Q0; c the symbols B, S, and P represent bending, axial stretching and pressurising; r represents the inner (luminal) surface of the loaded state configuration Q2 of the segment 8 Urevc, J. - Brumen, M. - Flis, V. - Stok, B. Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 dorr dr = 0. (7) The natural boundary conditions, Eq. (8), are applied both to the inner and outer surface of the segment. The inner surface is subject to the intraluminal pressure P, whereas the outer surface is considered traction-free: z — n atjz, srz — n arz, £rj — n arj, (14) Gp GT GT where Sj and a^ are the strain and the stress tensor. Next, the components of the strain tensor are expressed in terms of the radial and axial displacements u(r) and w(r) [26]: du u £rr = ~t , £(b(b = , dr rr r dw ezz = — = e = const, zz dz sr(=srz = s(z = (15) Using Eqs. (14) and (15), the following relationships between the stress tensor components and radial displacement u(r) are obtained: du u arr = Ap — + Apt — + Ajp e - ßpT, dr r du u °J ■ /aT jrJ Fig. 5. Methodology that predicts the stretched state (Q ^) of a realistic arterial segment (AS) is known, the stretched state of the circular arterial segment Q^ can be computed. Next, the thermomechanical analogy to predict the residual stresses is defined on the basis of the obtained Qii state. This is done by identifying the fictitious circular segment corresponding to the characterised circular arterial segment. The initial geometry of is defined based on the surface C r (see Fig. 5) whereas the thickness and material properties of should be identified so that the loaded state of the fictitious circular segment ^^ resembles the stretched state of the circular arterial segment Q^. So far in this subsection, the stretched state of the AS has been quantified on the basis of a simplified AS geometry assumed to have the form of a circular tube. Based on this simplification, the parameters that characterise the fictitious circular segment have been obtained. Next, the fictitious AS model needs to be defined. This is done by applying the identified parameters of the fictitious circular segment (represented as FIC in Fig. 5) to the luminal surface of the AS, i.e., the surface rAS. With thus defined state, the thermomechanical response of the fictitious AS (T^f) can be computed. Note that in general this computation, i.e., solving the thermomechanical problem of the fictitious AS model, cannot be performed analytically because of the irregular geometry of the AS model. An approximate method, such as FEM, needs to be considered. Finally, based on the results of the Tstate the AS stretched state Q AiS is approximated - note that AS 1X O O AS O1X is just defined based on the state TJS and thus still includes the nonlinearities of the arterial AS mechanical response. The geometry of the O a state Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State 11 Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 equals the geometry of the state, and the stresses of the Q ^ state equal the stresses of the state. 3 APPLICATION OF THE THERMOMECHANICAL ANALOGY TO THE ARTERIAL MODEL The methodology for determining the residual stresses in the arterial segment, which were discussed thoroughly in the previous section, is generalised in this section to enable a prediction of the arterial response based on the in vivo arterial data. The individual steps of the methodology are demonstrated considering the AM, defined in Section 1 (Fig. 1). To give an overview of the methodology, let us consider the flowchart presented in Fig. 6: (i) the starting data, needed for a prediction of the AM stretched state Qp^, are the AM loaded geometry, the corresponding intraluminal pressure, the in vivo stretch X, and the material properties of the AM. (ii) The stretched state of the AM is quantified based on characterising the mechanical response of two idealised (circular) arterial segments, corresponding to the arterial smallest (m) and largest (M) luminal cross-sectional area (see Fig. 1). (iii) The fictitious arterial model (FAM), i.e., the fictitious model of the AM, is defined next. This is done by defining the fictitious circular segments for the aforementioned circular arterial segments. The identified parameters of the fictitious circular segments are then interpolated/ extrapolated over the AM luminal cross-sectional area A(z) and applied to the AM luminal surface r^4. In this way, the FAM in its zero-stress state is obtained. Its loaded state Tj^ is then predicted by solving the corresponding steady-state thermomechanical problem. (iv) Finally, the AM stretched state Q is defined based on the obtained FAM T^M state (i.e., based on the T^M geometry and stress state). Further on, the loaded state of the AM Q , as well as the AM zero-stress state Q , are predicted. The individual steps of the methodology are presented in more detail in the following subsections. 3.1 Quantification of the AM Stretched State The quantification of the AM stretched state is performed based on the observation of two idealised (circular) arterial segments that refer to the arterial smallest (m) and largest (M) luminal cross-sectional area, see Fig. 1b. The respective segments' geometry is obtained from the AM cross-sectional area A(z) on the area equivalence basis: A(z) = na(z)b(z) = n(ri (z))2, (20) Fig. 6. The methodology that predicts the mechanical response of the artery based on the arterial in vivo geometry (I.e., without the arterial zero-stress state being known); the symbols P, T, CT, and CL represent pressurising, thermal loading, transversal cut, and longitudinal cut, respectively 14 Urevc, J. - Brumen, M. - Flis, V. - Stok, B. Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 from where the equivalent radius rt (z ), see ri Fig. 2, is computed, and in particular, k e {m, M} for segment (m) and (M), respectively. Next, the zero-stress state configurations of both segments, (m) and (M), are separately identified (see Section 2.1). The identification results are tabulated in Table 1, where R(otc} and 00k} ; k e{m,M} are respectively the outer radius and the opening angle for the segment (m) and (m). With the zero-stress state of both segments determined, their stretched states, designated as Q^ and QjM for the segment (m) and (M), respectively, can be computed. Table 1. The identified properties of the circular arterial segments Segment: k (m) (M) (m ) RO] [mm] 4.45 5.34 4.9 00k > [°] 132.89 129.94 131.34 3.2 Fictitious Arterial Model As discussed in Section 2.2, the AM stretched state ^^ can be predicted by defining an associated fictitious arterial model, designated as FAM and solving the corresponding thermomechanical problem. The advantage of this approach is that the arterial response is predicted without the arterial zero-stress state being known. How to define the FAM is discussed in detail in the following subsections. 3.2.1 Thermomechanical Characterisation of the Fictitious Circular Segments The thermomechanical analysis of the FAM is basically defined by the characterisation of the fictitious circular segments corresponding to the circular arterial segments (m) and (M), as presented in step 2 in Fig. 5. Note that in this way, essentially, the properties of the FAM for its location A(m) and A(M), see Fig. 1, are defined. The variables that need to be identified in order to characterise both fictitious segments (see Section 2.3) are: the material properties, the thicknesses dfk, the thermal loadings ATfk}, and AT0(i}; k e{m,M} for segment (m) and (M), respectively. During the identification of the unknowns, the material properties of both segments were initially assumed to be equal because they represent the material properties of the FAM. However, what was found is that in order for the stresses of both fictitious segments to correspond to the stresses of the associated arterial segments, the coefficient of thermal expansion aT has to be identified separately for each segment. Because of a relatively large number of unknowns, the identification procedure was performed in the following way. First, the identification of the parameters defining the fictitious segment (m) (see Section 2.2.1.) is performed. The unknowns to be determined: thickness 5(m), material constants EP, aP, Et and aTm), and thermal loadings A7^m) and A7^m), are identified by targeting the Qstate of the circular arterial segment (m). Next, the identification of the fictitious segment (M) is performed by adopting the material properties identified for the segment (m). The remaining variables to be determined: AT(M), A7;(m) , and a(M, are obtained by targeting the QM state of the circular arterial segment (M). The result of the identification is given in Table 2. Table 2. The identified properties of the fictitious circular segments Segment: k (m) (M) AT-k -1 [°] 1.0059 1.069 AT{0k k [°] 1.255 1.295 f > [mm] 0.6931 0.71 4 > [1/K] -0.05 -0.048 EP [kPa] 299.29 Et [kPa] 299.29 aP [1/K] -0.1388 3.2.2 Thermomechanical Characterisation of the Fictitious Arterial Model As mentioned, the parameters of the FAM are defined based on the parameters identified for the fictitious circular segments (m) and (M). The parameters that are common to both fictitious segments are also assumed constant in the FAM. These are: EP, aP, ET, GP, Gt, vp, vpt, k, c and p (see Section 2.2). However, some parameters of the fictitious segments, namely their thickness cf, thermal loadings ATi and ATo , and the coefficient of thermal expansion aT, are not the same in both fictitious segments. We will refer to these variables as the varying variables because they depend on the cross-sectional area of the fictitious circular segments. Note that the cross-sectional area of the FAM also varies - from A(m) to A(M). Therefore, Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State 11 Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 in order to apply the varying variables to the FAM, their values need to be generalised (interpolated) over the cross-sectional area A(z) of the FAM. In this paper, a linear interpolation of the varying variables across the A(z), i.e., from A(m) to A(M), is adopted. For example, a linear interpolation of the variable aTAM (where the superscript "AM" refers to the interpolated variable) as a function of A(z) can be written as: arM (A) = aF> Ym (A) + + 4M) Ym (A) = a£M (z), (21) where ((m)) and and the shape functions Ym (A) and YM (A) are A _ A(m) defined as: Ym (a) = 1 - A(M) _ A(m) =Ym (z) A - A(m) and YM (A) = a(M) - A(m) = YM (Z) . Note that a^M(A) (as well as Ym(A) and YM(A)) is essentially a function of the coordinate z, thus aTAM(z) because A = A(z) (see Fig. 1). The same interpolation approach is used to define the remaining varying variables: SfM (z), ATiAM (z) and AToAM (z). Their variations as a function of the coordinate z are plotted in Fig. 7. Remark: For the purpose of defining the properties of the FAM, it would be possible to characterise more than two circular arterial segments. For these segments, the corresponding fictitious segments could be obtained and characterised, thus allowing for a higher-order interpolation used in the FAM characterisation. However, as shown in Section 4.1.1, the linear interpolation, which is used to define the FAM properties, is sufficient. (A) (M) 0 50 100 150 i z [mm] Fig. 7. Variations of the FAM properties: a^(z ), (z ), atAM (z), and AToAM (z) (see Eq. (21)), along the coordinate z 3.2.3 Construction of the Fictitious Arterial Model (FAM) Considering in vivo Measured Data Based on the characterisation of the FAM parameters that were determined in previous subsections under the assumption of idealised circular arterial segments, it is now possible to construct a computational model of the FAM taking in vivo measured data into account. Accordingly, the geometry of the FAM is constructed based on the inner surface of the AM r^M (Fig. 1), as shown in Fig. 8, considering its variation rAM along the z direction. As discussed, we define the inner surface of the FAM to be the same as the inner (luminal) surface of the associated AM, whereas the outer surface of the FAM is obtained by applying the thickness 5AM (z) (Fig. 7) to the surface rfM . Solving the steady-state thermomechanical problem Fig. 8. Constructing the fictitious arterial model 16 Urevc, J. - Brumen, M. - Flis, V. - Stok, B. Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 Note that the thickness of the FAM is assumed to be constant in the circumferential direction of the model while changing in the z direction. As assumed, the material behaviour of the FAM is transversely isotropic with respect to the z axis. Some of the material properties are constant (see Table 2) whereas aATM(z) is given in Fig. 7. The remaining material parameters that are not presented in Table 2 are given in Section 2.2.1. 3.2.4 Loaded State of the Fictitious Arterial Model In order to evaluate the residual stresses of the AM, the above-defined FAM is thermally loaded by the temperature changes A7}am (z) and AToAM (z) applied respectively on the inner and outer surfaces of the model (Fig. 8). Note that both ATiAM (z) and ATo (z ) (Fig. 7) are assumed to be constant in the circumferential direction of the model, while changing in the z direction. The stress distribution in the loaded state of the FAM T^ , obtained by solving the given thermomechanical problem with associated mechanical and thermal boundary conditions (see Section 2.2), yields an approximation of the residual stresses distribution in the AM and its configuration. 3.3 Modelling the Mechanical Response of the Arterial Model with Residual Stresses Considered The geometry and the stresses of the loaded state TjM of the FAM represent an approximation to the geometry and stresses of the arterial model in its stretched state Q^. The computational model of the AM stretched state Q^ is thus specified by adopting the geometry of the Tj^4 configuration and the given arterial material behaviour. For the initial stress state of the thus-defined model, representing, actually, the residual stresses in the artery, the stresses from T are adopted whereas the stretches corresponding to the Q^ configuration with the initial stresses considered are obtained by taking into account the arterial material model, Eq. (1). From the arterial stretched state all the remaining arterial states can be predicted (see the lower grey box in Fig. 6). The loaded state Q^M of the AM is obtained by inflating the AM from its stretched state. The AM unloaded state Q is predicted by releasing the axial stretch A from the AM stretched state, whereas the approximation of the AM zero-stress state is obtained by further releasing the constraints on one of the AM symmetry surfaces (Fig. 2). Remark: The reason why the arterial stretched state Q^ is predicted with the fictitious model, and not, for instance, by considering the arterial unloaded state QAM , is the following: normally, the goal of the arterial characterisation is to predict the arterial loaded state. Therefore it is more convenient to predict the arterial loaded state from the arterial stretched state given that in between both states there is only the inflation process (see Fig. 2). For example, by predicting the arterial loaded state from the arterial unloaded state, the artery first needs to be stretched for A in order to obtain its stretched state, and can only later be inflated. 4 RESULTS, VALIDATION OF THE METHODOLOGY AND DISCUSSION In this section, the developed methodology of determining the arterial response based on the in vivo measured arterial data and taking the residual stresses into account is applied in the analysis of the given case study (Fig. 1). In the analysis, regardless of whether the computation applies to the real or the fictitious model, the computations that relate to the circular models are performed using the finite difference method while the FEM is used, because of the complexity of respective model geometries, for the computations related to the arterial models. First, after giving a brief description of the numerical model in Section 4.1, the results for the circular arterial and fictitious segments are presented in Section 4.2. The results regarding the arterial response are presented in Sections 4.3.1 and 4.3.2 for the AM stretched and loaded state, respectively. Finally, at the end of the section, the methodology is validated based on the prediction of the AM zero-stress state. 4.1 FEM Modelling Taking into consideration that the geometry of the FAM loaded state Tj^ configuration may be considered to be an approximation to the geometry of the AM stretched state Qp^ configuration, only one single FEM discretisation is needed for all of the computations regardless of whether the computation applies to the real or the fictitious model. The FEM model is constructed based on the geometry of the FAM in its stress-free configuration ¥AM (Fig. 8) taking into account the data specified in Section 3.2.3. In principle, because the considered case study (Fig. 1) is symmetric with respect to two perpendicular planes (note that the cross section of the Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State 11 Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 AM is elliptical), only one quarter of the model could be considered for computations. However, in order for the computation of the AM zero-stress state to also be performed (see Fig. 2), half of the model is considered for the FEM discretisation with the x - z plane being the plane of symmetry. The FEM model is built of 17640 eight-node brick elements with displacement and temperature degrees of freedom. The conceived FEM model is used first for the computation of the FAM loaded state • For this purpose, a linearly elastic, transversely isotropic behaviour of the material is assumed taking the data specified in Section 3.2.3 into account. The thermomechanical analysis is performed using the commercial FEM software Abaqus/Standard [27] under prescribed boundary and loading conditions. The front and rear surfaces of the model are fixed in the z direction whereas they are free to move in the x -y plane. On the same surfaces, the adiabatic thermal conditions are applied. Moreover, proper mechanical and thermal conditions are applied on the symmetry surfaces (see Section 3.2.4). The stresses of the FAM loaded state represent an approximation of the residual stresses in the artery. The computational model of the AM stretched state Qu^ is thus specified by adopting the geometry of the configuration and related stresses as initial stresses, and, importantly, taking the given arterial material behaviour into account, Eq. (1). Clearly, the stretches corresponding to the configuration are readily computed considering the constitutive relations and the initial stresses. When the thus-defined computational model Q ^ is inflated with the intraluminal pressure P, the loaded state Qam is computed. If, on the contrary, the constraints on the axial displacement are removed and the axial stretch A is released, the unloaded state can be computed. Finally, by further releasing the constraints on one of the AM symmetry surfaces, the approximation of the zero-stress state QAM can be attained. 4.2 Arterial and Fictitious Circular Segment The stretched states of the circular arterial segments (m) QC(m) and (M) qc(m) are presented in Fig. 9. Both states are computed as presented in Section 2.1 based on the identified data of their zero-stress states (Table 1). As seen in Fig. 9, the stresses predicted for both segments are in the same order of magnitude. The circumferential stresses are somewhat bigger for the segment (m) whereas the axial stresses are higher in the segment (M). The same graphs also represent the stresses predicted for the corresponding fictitious segments ^C(m) and as well as for the fictitious circular segments (m) and (M), respectively. Despite a minimal deviation between the axial stresses (see Fig. 9), an excellent agreement of the stress distributions and geometries between the circular arterial segments () and the fictitious circular segments (is obtained. The results presented demonstrate that the thermomechanical state, obtained by thermally loading a fixed, closed (fictitious) tube, is equivalent to the stretched state that is obtained by bending and stretching an open arterial section. As discussed, the fictitious circular segments have essentially been defined in order to characterise the thermomechanical procedure and to apply the determined parameters of fictitious segments to the AM. So far, we have demonstrated that the stretched state of the arterial segments (m) and (M) can be a) 19 14 ra Q a. y 3 2 -6 vpC(rn) MX - axial stress circumferential stress - . radial stress 3.2 r [mm] r [mm] Fig. 9. Stresses computed for the ^C(m) and states of the circular arterial segments (full line), and their comparison with the stresses obtained for the fictitious circular segments ^C(m) and (broken line); the results for the segment (m) are given in a) whereas b) shows the results for the segment (M) 18 Urevc, J. - Brumen, M. - Flis, V. - Stok, B. Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 obtained with the fictitious segments (m) and (M), respectively. In the next subsection, we shall assess the used linear interpolation that was applied to interpolate the properties of the fictitious segments with the purpose of defining the FAM (see Section 3.2.2). 4.2.1 The Adequacy of the Considered Interpolation The variables cc£M(z), (z), ATiAM (z) and (z) were obtained by linearly interpolating the corresponding variables that have been determined for the fictitious segments (m) and (M) (see Section 3.2.2). The adequacy of this approximation is assessed by characterizing an arterial segment, denoted by (m), with the equivalent radius r m _ _ Cr- em) + ~ (M) i i )/2 - C(m ) R 1S For the segment (m), its stretched state Q identified and then compared to the stretched state, obtained through defining an associated fictitious segment (state). The parameters of the fictitious segment (m) are obtained by considering the data in Fig. 7 and Table 2. The loaded state of the fictitious segment (^j^111"1) is computed as presented in Section 2.2. To determine the stretched state of the arterial segment (m ), its zero-stress configuration needs to be identified as presented in Section 2.1. The identified values of the unknowns R^ and are given in Table 1. The computed stresses of the arterial segment ( ) and the fictitious segment ( ) are presented in Fig. 10. ,C(m) r (mm) Fig. 10. The calculated stresses of the QC(m'1 state of the circular artery segment (m) (full line), compared to stresses of the state, computed with the associated fictitious circular segment (broken line) a) [kPa] 6.4 _.4 -.2 0 -2 -3 ■--6.15 (% view) (Mx) fit r-hm 1 il ■ rn IIBc —1 1 1 ..... 11 r-4 (MS I The transmural stress distribution over the lines (mx), (m ), (Ms) and (M ) is shown in b). b) n^ (mx) (my) ---axial stresses -circumferential stresses 19 14 ■ 9 : 0.60 fiF/M> (Mx) (My) ---axial stresses -circumferential stresses 0.30 0.40 r [mm] Fig. 11. The predicted AM loaded state (Q ): a) 3D view and b) the transmural distribution of the circumferential and axial stresses at the AM positions (mx), (my), (Mx), and (My); for their location see Fig. 11a Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State 11 Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 Despite a slight discrepancy between the axial stresses (Fig. 10), the agreement between the obtained arterial and fictitious state of the segment (m) is excellent. Based on these results, it is possible to conclude that the first-order interpolation, which is used to interpolate the properties of the fictitious segments, is justified. 4.3 The Predicted Mechanical Response of AM 4.3.1 AM Stretched State The circumferential stresses o^ of the predicted AM stretched state (Qij^) are presented in Fig. 11a. The transmural variation of the o^ shows the bending distribution (see the zoomed-in detail along the wall thickness), as was also observed for the circular arterial segments. In contrast to the straight circular arterial segment, the stresses in the AM vary along the AM length because of the varying AM cross-sectional area. The transmural distribution of the stresses at the AM positions (m) and (M) are shown in Fig. 11b (red and black curves). In the figure, the stresses are presented at the position of the AM cross-sectional x and y axis. Note that the stresses at position (m) are higher on the (mx) axis compared to the values on the (my) axis. At position (M), on the other hand, the opposite is observed. This is interesting because the x axis represents the cross-sectional semi-minor axis at both position (m) as well as at position (M) (see Fig. 1). The reason for this can be probably attributed to the variable cross-sectional area of the AM. Next, let us compare the transmural distribution of the stresses at the AM locations (m) and (M) with the stresses predicted for the circular arterial segments and . The stress distribution in Fig. 11b that corresponds to the circular arterial segment is shown in blue. At position (m), the stresses predicted for the AM are somewhat higher in comparison to the stresses of , whereas the opposite is true for the position (M). Such an error between the stresses predicted for the AM and for the circular arterial segments is a result of the AM variable cross-sectional area. The behaviour of the individual AM cross sections or segments is limited by the behaviour of the whole AM. Therefore, the individual AM segments cannot extend freely as is the case of the circular arterial segment. The error between the maximal circumferential stresses predicted for the circular arterial segments and those predicted for the AM is around 30%. 4.3.2 AM Loaded State The predicted AM loaded state (Q^) is presented in Fig. 12a. The stresses, as seen in the figure, are somewhat higher in the areas where the arterial cross-sectional area is larger. The maximal circumferential stress o^ is 73 kPa, which appears on the outer surface of the AM wall at the AM location (My) (see the graphs in Fig. 12b). The range of stresses o^, as well as their transmural distribution (Fig. 12b), agree with the corresponding data from literature [24] and [28]. Next, let us see the stresses that would be predicted for the AM sections (m) and (M) if they were treated as straight circular segments. The loaded state of the circular arterial segments (m) and (M), presented as (m) and (M), respectively, are shown in blue in Fig. 12b. The loaded states of both circular arterial segments (m) and (M) are computed as presented in Section 2.1. Similar to what has been observed in the previous subsection, the stresses predicted for the AM at the position (m) are somewhat higher in comparison to QC(m), whereas the opposite is true for the position (M). The reason for this lies in the AM variable cross section which prevents the individual cross sections of the AM to extend freely. The error between the stresses predicted for the circular arterial segments and for those of the AM can be as high as 15%. Given the size of the error, the suitability of using the circular arterial segments in order to determine the stress states of arteries, as is common in literature, is questionable. 4.3.3 Predictions of the AM Zero-Stress State A way to partially verify the adequacy of the presented method is to simulate the opening angle experiment. If the stresses predicted for the arterial opened-up configuration are almost null, the method can be considered adequate. The same approach, for instance, can be seen in [2]. The zero-stress state of the AM (^oM) is predicted based on the AM stretched state () in two steps. First, the axial constraint that is applied to the state is released. The AM thereby contracts axially and occupies its unloaded state QjAM. In addition, the AM zero-stress state is obtained with a longitudinal cut which is performed in a numerical analysis by releasing the constraint on one of the AM symmetry surfaces (located on the x - z plane, see Fig. 13a). The predicted equilibrium state of the AM opened-up configuration, obtained with a single longitudinal cut is designated as and 20 Urevc, J. - Brumen, M. - Flis, V. - Stok, B. Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 un Q2M (Va view) (Mx) I The transmural stress distribution over the lines I I (mx), (m ), (MJ and (M ) is shown in (b). (b) > 80 n 75 70 ; ■-V 65 : M te 60 55 (mj (,u i ■ axial stresses ■ circumferential stresses 50 0.00 0.10 0.20 0.30 0.40 0.50 r (mm) axial stresses circumferential stresses 0.00 0.10 0.20 0.30 0.40 0.50 r (mm) Fig. 12. The predicted AM stretched state ( Qj^ ): a) 3D view and b) the transmural distribution of the circumferential and axial stresses at positions (mx), (my), (Mx), and (My): their location is presented in Fig. 12a (a) aeq (kPa) 1.2 1 — 0.8 — 0.6 0.4 É-0.2 0 (b) Q AM #1 Q AM #2 Fig. 13. Predictions of the AM zero-stress state: a) zero-stress state Q^^1, obtained by longitudinal cutting of the AM unloaded state, and b) zero-stress state , obtained by transversely cutting the predicted AM Q^^1 state into smaller segments is presented in Fig. 13a. The maximal von Mises equivalent stresses in QAM#1 are around 1.2 kPa (Fig. 13a). These stresses occur locally in the areas where the AM cross-sectional area is larger. As seen in Fig. 13a, the predicted AM QAM#1 state is not stressfree. The reason for this is that the AM variable cross- sectional area limits the AM's ability to extend freely, and consequently, it prevents the AM from occupying its zero-stress configuration. In order to obtain a better estimation of the AM zero-stress state, the predicted AM QitM#1 state would have to be cut into smaller segments. This procedure is presented next. Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State 11 Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 The predicted AM zero-stress state Q^M#2, which is obtained by cutting the Q^^1 state into 10 separate sections, is presented in Fig. 13b. As seen in the figure, the maximal stresses obtained for the individually cut segments ( Q^M#2 state) are less than 1 kPa. These stresses occur in the second, third, and fourth segments (from both ends) whereas the maximal stresses range up to 0.8 kPa in the remaining ones (the first and the fifth segment from both ends). For these segments, the stresses are lower by approximately 30% compared to the AM QAM#1 state. Higher stresses that appear in the second, third, and fourth segments are the result of the shape of the segments, which is, if compared to the first segment, somewhat curved in the z direction. Based on the results presented, however, and given that the stresses of the AM Q^M#2 state are relatively low, it is possible to conclude that the use of the thermomechanical analogy in order to enforce the residual (bending and stretching) stresses into arterial model is justified. 5 CONCLUSION In this paper, a method is presented which determines the arterial mechanical response without knowing the arterial opened-up zero-stress configuration. The arterial configuration on which the method is based is the arterial loaded configuration. We approach the problem of predicting the arterial residual (bending and stretching) stresses by considering a thermomechanical analogy. These stresses are determined by thermally loading a clamped, closed tube (called a fictitious model). The initial geometry of the fictitious model is constructed based on the luminal surface of the observed artery. The consequence of the applied thermal loading is that the fictitious model shrinks and the thermal (residual) stresses occur. In this paper, we have demonstrated that the loaded state of the fictitious model (i.e., its geometry and its stress state) is equivalent to the stretched state of the arterial model (see Section 4.1). However, in order for the presented analogical procedure to be adequate, the parameters of the fictitious model need to be properly defined. The adequacy of the method has been verified by computing the zero-stress state of the arterial model. As presented in Fig. 13, the equilibrium configurations of the segments that are obtained following a longitudinal and transversal cut of the arterial model have the form of an opened-up tube. Furthermore, the residual stresses remaining in the segments are fairly mild compared to the stresses of the arterial stretched state. The results obtained show that the residual stresses predicted for the arterial stretched state with the thermomechanical analogy exhibit what is known as the bending distribution. This has proved the method to be adequate. An interesting aspect that aroused in the analysis of the results is suitability of using the straight circular arterial segment for the characterisation of the arteries. This is a common approach in the literature [19], [28] and [29]. The segment is considered straight and stand-alone, meaning that the effect of the remaining artery in the segment's axial direction is not taken into account. Based on our results (see Section 4.2.2), the error between the stresses that are predicted for a straight circular arterial segment, and the stresses predicted for the arterial model, can be as high as 15%. Additionally, the gradient of transmural distribution of the stresses for both cases differs (see Fig. 12). Nevertheless, it has to be said that a realistic artery does not exhibit such a variable cross-sectional area as was considered in this article. In a realistic artery, in fact, the influence of neglecting the effect of the remaining artery is smaller when observing one straight circular segment as a separate unit. The results of the methodology presented in the paper show a distinction between the mechanical state predicted for a circular segment that is observed as a separate unit and the mechanical state that is predicted when the arterial variable cross-sectional area is taken into account. This topic therefore warrants further research, and presents one of the goals of our further work. Moreover, some of the arterial properties that have been neglected in this work, such as their multi-layered structure or associated longitudinal residual stress, are planned to be implemented further in the methodology. As for the clinical applicability of the presented methodology, the methodology could - supported by the medical imaging techniques - serve as a useful tool for predicting the patient-specific arterial mechanical response. 6 REFERENCES [1] Holzapfel, G.A., Sommer, G., Auer, M., Regitnig, P., Ogden, R.W. (2007). Layer-Specific 3D residual deformations of human aortas with non-atherosclerotic intima thickening. Annals of Biomedical Engineering, vol. 35, no. 4, p. 530-545, DOI:10.1007/s10439-006-9252-z. [2] Alastrue, V., Pena, E., Martinez, M.A., Doblare, M. (2007). Assessing the use of the "opening angle method" to enforce residual stresses in patient-specific arteries. Annals of Biomedical Engineering, vol. 35, no. 10, p. 1821-1837, DOI:10.1007/s10439-007-9352-4. 22 Urevc, J. - Brumen, M. - Flis, V. - Stok, B. Strojniski vestnik - Journal of Mechanical Engineering 61(2015)1, 5-23 [3] Vossoughi, J. (1992). Longitudinal residual strains in arteries. Proceedings of the 11th Southern Biomedical Engineering Conference, Memphis. [4] Wang, R., Gleason, R.L. (2010). A mechanical analysis of conduit arteries accounting for longitudinal residual strains. Annals of Biomedical Engineering, vol. 38, no. 4, p. 13771387, DOI:10.1007/s10439-010-9916-6. [5] Fung, Y.C., Liu, S.Q. (1989). Change of residual strains in arteries due to hypertrophy caused by aortic constriction. Circulation Research, vol. 65, no. 5, p. 1340-1349, D0I:10.1161/01.RES.65.5.1340. [6] Humphrey, J.D. (2002). Cardiovascular Solid Mechanics: Cells, Tissues, and Organs. Springer-Verlag, New York, D0I:10.1007/978-0-387-21576-1. [7] Bazilevs, Y., Hsu, M.C., Zhang, Y., Wang, W., Kvamsdal, T., Hentschel, S. Isaksen, J.G. (2010). Computational vascular fluid-structure interaction: methodology and application to cerebral aneurysms. Biomechanics and Modeling in Mechanobiology, vol. 9, no. 4, p. 481-498, D0I:10.1007/ s10237-010-0189-7. [8] Fung, Y.C. (1991). What are the residual stresses doing in our blood vessels? Annals of Biomedical Engineering, vol. 19, p. 237-249, D0I:10.1007/BF02584301. [9] Chuong, C.J., Fung, Y.C. (1983). Three-dimensional stress distribution in arteries. Journal of Biomechanical Engineering, vol. 105, no. 3, p. 268-274, D0I:10.1115/1.3138417. [10] Chuong, C.J., Fung, Y.C. (1986). On residual stresses in arteries. Journal of Biomechanical Engineering, vol. 108, no. 2, p. 189-192, D0I:10.1115/1.3138600. [11] Takamizawa, K., Hayashi, K. (1987). Strain energy density function and uniform strain hypothesis for arterial mechanics. Journal of Biomechanics, vol. 20, no. 1, p. 7-17, D0I:10.1016/0021-9290(87)90262-4. [12] Bustamante, R., Holzapfel, G.A. (2010). Methods to compute 3D residual stress distributions in hyperelastic tubes with application to arterial walls. International Journal of Engineering Science, vol. 48, p. 1066-1082, D0I:10.1016/j. ijengsci.2010.06.005. [13] Greenwald, M., Moore, J., Rachev, A., Kang, T., Meister, J.J. (1997). Experimental investigation of the distribution of residual strains in the artery wall. Journal of Biomechanical Engineering, vol. 119, no. 4, p. 438-444, D0I:10.1115/1.2798291. [14] Gee, M.W., Förster, Ch., Wall, W.A. (2010). A computational strategy for prestressing patient-specific biomechanical problems under finite deformation. International Journal for Numerical Methods in Engineering, vol. 26, no. 1, p. 52-72, D0I:10.1002/cnm.1236. [15] Govindjee, S., Mihalic, P. (1996). Computational methods for inverse finite elastostatics. Computer Methods in Applied Mechanics and Engineering, vol. 136, no. 1-2, p. 47-57, D0I:10.1016/0045-7825(96)01045-6. [16] Lu, J., Zhou, X., Raghavan, M.L. (2008). Inverse method of stress analysis for cerebral aneurysms. Biomechanics and Modeling in Mechanobiology, vol. 7, no. 6, p. 477-486, D0l:10.1007/s10237-007-0110-1. [17] Delfino, A., Stergiopulos, N., Moore, J.E.Jr., Meister, J.J. (1997). Residual strain effects on the stress field in a thick wall finite element model of the human carotid bifurcation. Journal of Biomechanics, vol. 30, no. 8, p. 777-786, D0I:10.1016/ S0021-9290(97)00025-0. [18] Raghavan, M.L., Trivedi, S., Nagaraj, A., McPherson, D.D., Chandran, K.B. (2004). Three-dimensional finite element analysis of residual stress in arteries. Annals of Biomedical Engineering, vol. 32, no. 2, p. 257-263. [19] Holzapfel, G.A., Ogden, R.W. (2010). Modelling the layer-specific three-dimensional residual stresses in arteries, with an application to the human aorta. Journal of the Royal Society Interface, vol. 7, no. 46, p. 787-99, D0I:10.1098/ rsif.2009.0357. [20] Schulze-Bauer, C.A.J., Holzapfel, G.A. (2003). Determination of constitutive equations for human arteries from clinical data. Journal of Biomechanics, vol. 36, no. 2, p. 165-169, D0I:10.1016/S0021-9290(02)00367-6. [21] Masson, I., Boutouyrie, P., Laurent, S., Humphrey, J.D., Zidi, M. (2008). Characterization of arterial wall mechanical behavior and stresses from human clinical data. Journal of Biomechanics, vol. 41, no. 12, p. 2618-2627, D0I:10.1016/j. jbiomech.2008.06.022. [22] Stalhand, J., Klarbing, A., Karlsson, M. (2004). Towards in vivo aorta material identification and stress estimation. Biomechanics and Modeling in Mechanobiology, vol. 2, no. 3, p. 169-186, D0I:10.1007/s10237-003-0038-z. [23] Stalhand, J. (2009). Determination of human arterial wall parameters from clinical data. Biomechanics and Modeling in Mechanobiology, vol. 8, no. 2, p. 141-148, D0I: 10.1007/ s10237-008-0124-3. [24] Cardamone, L., Valentin, A., Eberth, J.F., Humphrey, J.D. (2009). Origin of axial prestretch and residual stress in arteries. Biomechanics and Modeling in Mechanobiology, vol. 8, no. 6, p. 431-446, D0I:10.1007/s10237-008-0146-x. [25] Abaqus 6.7.1 (2008). Users Manual. Simulia, Providence. [26] Birger, B.I. (1972). Temperature stresses in an anisotropic cylinder, NASA, Washington. [27] Fonck, E., Prod'hom, G., Ray, S., Augsburger, L., Rufenacht, D.A., Stergiopulos, N. (2007). Effect of elastin degradation on carotid wall mechanics as assessed by a constituent-based biomechanical model. American Journal of Physiology - Heart and Circulatory Physiology, vol. 292, no. 6, p. H2754-H2763, D0I:10.1152/ajpheart.01108.2006. [28] Stalhand, J., Klarbing, A. (2005). Aorta in vivo parameter identification using an axial force constraint. Biomechanics and Modeling in Mechanobiology, vol. 3, no. 4, p. 191-199. [29] Masson, I., Beaussier, H., Boutouyrie, P., Laurent, S., Humphrey, J.D., Zidi, M. (2011). Carotid artery mechanical properties and stresses quantified using in vivo data from normotensive and hypertensive humans. Biomechanics and Modeling in Mechanobiology, vol. 10, no. 6, p. 867-882, D0I:10.1007/s10237-010-0279-6. Applying Thermomechanical Analogy to Predict the Arterial Residual Stress State 11