International Academic Institute for Science and Technology International Academic Journal of Science and Engineering Vol. 3, No. 5, 2016, pp. 57-67. ISSN 2454-3896 International Academic Journal of Science and Engineering www.iaiest.com Effects of Aging on the Mechanical Behavior of Human Arteries in Large Deformations Amin Safi Jahanshahi a,c, Alireza Saeedi b a PhD student, Shahid Bahonar Universiy, Kerman, Iran. b Professor, Shahid Bahonar University, Kerman, Iran. c Lecturer Sirjan University of Technology, Sirjan, Iran. Abstract The mechanical properties of living cells related to a lot of physical and Physiological factors. Aging is one important factor in changing physical and mechanical properties of living cells. The mechanical behavior of living cells and live tissues is nonlinear and their deformation is large. Finding a suitable mechanical model that can predict this behavior can have a lot of applications in prevention and treatment of various diseases. In this paper first based on the theory of nonlinear elasticity kinematic, equations has been extracted. Mooney-Rivlin nonlinear model considered on characteristic equations and using results of the report elastic stress test of biaxial Elastic constants have been calculated for the arteries. For simplicity, we considered artery into a long homogeneous cylinder. Elastic constants calculated for different age shows that arteries of people become coated by aging and their flexible capability decreases dramatically. Keywords: Nonlinear Elasticity, Large Deformations, Arteries, Biaxial Stress Test, Aging Introduction: Arteries and veins consisted of three layers (Figure 1) (Meyers et al.2008). Their Tunica intima is covered by a layer of endothelial cells. These cells are creating an inner membrane by mechanical and electrochemical connections for blood vessels. In fact, this membrane forms the boundary between the blood and the main lining of the vessels. Endothelial cells have consisted of collagen, fibronectin and laminin respectively (Meyers et al, 2008; Jafee 1987; Cox1978; Karsaj and Humphrey 2012). 57
The Tunica media or the main lining of the vessels includes the most vessels body and consisted of muscle cells and elastic by categories formation of collagen. Figure 1: Cross section of the vessel includes intima, media and adventitia layer. In arteries closer to the heart which have large-diameter, Elastic percentage is more than the muscle and in the arteries away from the heart that have the smaller diameter, this ratio is reversed. The outer layer of type 1 collagen, nerve cells of fibroblast (a type of cell that has the capability to become collagen fibers.) and formed elastic fibers (Meyers et al, 2008; Jafee 1987; Karsaj and Humphrey 2012). In macroscopic analysis, we can consider the entire arteries lining as isotropic or anisotropic, homogeneous and incompressible material (Carew et al. 1968; Vaishnav et al. 1972; Hudetz1980; Chuong and Fung1984; Papageorgiou and Jones1987). Deformation of large living tissues and their mechanical behavior is nonlinear. For mathematical modeling of arteries, we can consider them as long cylinder and double-layer or single-layer (Maltzahn 1981; Maltzahn 1982; Maltzahn et al. 1984; Hudetz 1980; Holzapfel and Weizsacker 1998). More modeling conducted, they have considered artery lining as homogenous with low thickness (Taber 1998; Rachev2000; Gleason et al. 2004; Tsamis et al. 2009; Valentin et al. 2009; Wan 2010) while to estimate actual behavior arteries, we need to provide a model with high thickness. In this paper, we consider a single-layer model with the actual thickness and calculate the Mooney-Rivlin constants based on biaxial stress test for presented model. Constants obtained from the biaxial stress test calculated and presented in different ages. Kinematic Suppose that we show the position of a point in the coordinate transformation to vector space X, after changing configuration, we illustrate this point with the position vector x screw. ( ) (1) Deformation gradient tensor expresses as follows: 58
In Cartesian coordinates, deformation gradient tensor components change as: that and place components of vectors x and X in the deformed and unreformed coordinate respectively. Polar decomposition theorem express as follows: That R is an arbitrary symmetric tensor and U and V are orthogonal and positive definite tensor, we can say them as right and left stretch tensor respectively and we write them as the shape of spectrum in terms of specific vectors as follow: (2) (3) ( ) ( ) ( ) ( ) (4) Where * + are the principal stretches, ( ) i.e.are the unit eigenvectors of, the so-called Lagrangian principal axes, ( ) the unit eigenvectors of, the Eulerian principal axes, and denotes the tensor product. We can write Jacobin of transformation based on the principal stretches as follows: (5) We can calculate Left Cauchy-Green deformation tensor by using the following formulas: and right Cauchy-Green deformation tensor (6) Green strain tensor is defined as follows: ( ) (7) That has equal invariables by the left and right Cauchy-Green strain tensors as follows: ( ), ( )- ( ) (8) That ( ) represents the trace of second order tensor. For an incompressible substance, equation (5) will be as follows: Stress and Equilibrium Equations According to Figure 2, consider an Element of the media in deformed and unreformed state. In the absence of body forces, we have: (9) (10) 59
Where is forming the outer surface of the object after deformation Figure 2: transfer an Element before and after deformation We can write the recent equation as follows: Now, that ( ) and also the relation between nominal and Cauchy stress, We can rewrite equation (10) as: (11) (12) That Forms the outer surface of the body before reshaping and is the first piola stress tensor. By applying the divergence theorem on, we have: The second Piola stress tensor is another stress tensor we show here with P: Now tensor showed another stress tensor called the Biot stress tensor with T and we express as following form: ( ) ( ) (15) That R is an arbitrary orthogonal tensor that we mentioned previously in the case of polar decomposition analysis. Characteristic Equations Strain energy function density W for a homogeneous material is only a function of the deformation gradient F. If the material was inhomogeneous, strain energy function in addition to deformation gradient (13) (14) 60
will also hinge to different parts to X. If there were no internal constraints such as incompressibility, nominal stress and deformation gradient will be co-producing that their multiplications produce work. So we have simply: We write recent equation in component form as. For an incompressible material equation (16) comes as following: That p is the Lagrange multiplier imported by adding an incompressible bond. We can write equation (16) according to Cauchy stress as follows: The mentioned equation for an incompressible material will be as follows: A logical consequence obtains from isotropy is that specific Cauchy stress tensor values, and also Baiot stress tensor and tensor U are equal. So we can write by similar process whatever we applied by writing the (4) equation: (16) (17) (18) (19) ( ) ( ) ( ) ( ) (20) That and are the principal stresses of Biot and Cauchy. Therefore: (21) And for an incompressible material, we have: (22) The equation (20) will be as following for Piola stress type two and nominal stress: ( ) ( ) ( ) ( ) (23) Which and. Biaxial Stress Test In this section, we review homogeneous deformation that in general we can consider it as follows: 61
(24) Deformation gradient tensor is: [ ] (25) And by considering the incompressibility conditions, we have: [ ] [ ] (26) We consider the following invariables based on principal stretches: (27) It will be as following by applying incompressibility condition: (28) Consider W potential energy subordinate as follows: ( ) ( ) (29) By eliminating p from the equation (22-b), we can write: (30) If we consider the biaxial test, we have and therefore the above equation will become: (31) If we consider: ( ) ( ) (32) Equations (31) conclude: 0 1 (33) And so it obtains a similar equation for that finally after simplification we have: ( )( ) ( )( ) (34) 62
By solving the mentioned equations, we have and : ( )( ) ( )( ) (35-a) ( )( ) ( )( ) In the above equations for simplicity in writing, we have used energy function for Mooney-Rivlin material is as follows: (35-b) equation. Helmholtz free ( ) ( ) (36) We can write Helmholtz free energy function for the Neo-Hockean by putting above equation as follows: equals zero in the ( ) (37) Comparing the equations (35) with (30) we have: This means that by doing a biaxial stress test, we can easily achieve Mooney-Rivlin constants. By using results of biaxial stress test that conducted by Mohan and Melvin (Mohan and Melvin 1983), we can easily calculate constants. Suppose we define two new parameters of μ and α: (38) (39) In this case we can write for the Mooney-Rivlin material: ( ) ( )( ) (40) In this case, it is clear that by placing equal amounts of α, we will have neo-hock in equation: ( ) (41) Mechanical analysis of arteries In this analysis, we consider whole lining arteries a homogeneous and isotropic substance. Deformation field for cylinder under tension and inside pressure is as follows (Batra &Bahrami 2009) ( ) (42) 63
That is stretch along the length of the cylinder. Given the right malformed deformation, we can write gradient tensor as follows: [ ] (43) That. By applying incompressibility condition, we have:. /. /. / (44) Integrating the last equation, we have: In equation (45), A is integration constant that we must identify it. By calculating left Cauchy deformation tensor and putting it in the equation (6), we have non-linear structure of Mooney-Rivlin equation: (45) (46) (47) In the recent equations If the constants p, D, A, distinguished, stress distribution is completely clear and solving the problem is accomplished. Balance equation in polar coordinates is expressed as follows: (48) ( ) (49) And we consider boundary conditions as follows: ( ) ( ) (50) (51) That and are internal and external artery pressure respectively and are axial force. By applying boundary conditions and a little math equations, we finally get the following equation. 64
( ) 0 ( ) ( )1 (52) *. /. /. / 0 ( ) ( )1 +( ) (53) So equations (52) and (53) are two equation based on constants A, D, which we can calculate in this two unknown coefficients by solving simultaneously and finally stress distribution is fully characterized. Results The experimental results conducted by Mohan and Melvin using the biaxial stress test (Mohan and Melvin 1983). Mooney-Rivlin material constants calculated for artery in different ages and it displayed in Table 1 according to the equations (35). Table1. Mooney-Rivlin equation constants for artery in different ages ( ) ( ) Age(yaer) -3147 10396 25-3421 11128 49-5082 19400 60-6175 22725 87 In order to show the results, in Table 2 hypertension environmental tensions and radius for 8 mmhg at different radii based on the theory of nonlinear elasticity Mooney-Rivlin modeling compared by resulting of numerical solution. Table 2 comparison between the calculated environmental stresses and radial on non-linear Mooney-Rivlin and Numerical solution Result r(mm) T θθ(mpa) -T rr(kpa) Mooney-Rivlin Numerical Mooney-Rivlin Numerical 5.42 2.1150 2.0814 20.2175 19.8194 5.52 0.5725 0.5635 6.5791 6.3497 5.62 0.3386 0.3333 4.3555 4.2475 5.72 0.2440 0.2402 3.1095 3.0382 5.82 0.1929 0.1898 2.2310 2.1771 5.92 0.1609 0.1583 1.5476 1.5040 6.02 0.1389 0.1367 0.9861 0.9495 6.12 0.1230 0.1211 0.5086 0.4769 65
T θθ (MPa) t t ( M P a ) International Academic Journal of Science and Engineering, In Figure 3 environmental stress distribution along the radius at different ages for internal pressure 16 mm Hg is drawn. This figure shows that with increasing age, increasing environmental stress in arteries. The reason for this is that with age, the flexibility of the arteries decreases. 3 2. 5 25 2 49 1. 5 60 1 87 0. 5 0 5.4 5.5 5. 6 5. 7 5.8 5. 9 6 6.1 r ( m m ) r (mm) Figure 3: Distribution of environmental tension in terms of the radius at different ages The graph shows that the inner layers in contrast to the outer layers of arteries play more role in stress tolerance.in This article arteries lining considered as single-layer that due to stress distribution obtained it is suggested that we should use a two-layer to investigate stress distribution References: A. G. Hudetz, (1980). Continuum mechanical methods and models in arterial biomechanics, Advances in Physiology, 8(6): 223-232. A. Rachev, (2000). A model of arterial adaptation to alterations in blood flow, Journal of Elasticity, 61(1): 83 111. A. Tsamis,N. Stergiopulos,A. Rachev, (2009). A structure-based model of arterial remodeling in response to sustained hypertension, Journal of Biomechanical Engineering, 131(10): 101-104. 66
A. Valentin,L. Cardamone,S. Baek,J. D. Humphrey, (2009).Complementary vasoactivity and matrix remodelling in arterial adaptations to altered flow and pressure, Journal of. The Royal. Society. Interface, 6(32): 293 306. C. J. Chuong, Y. C. Fung, (1984). Compressibility and constitutive equation of arterial wall in radial compression experiments, Journal of Biomechanics, 17(1): 35-40. D. Mohan, J. W. Melvin, (1983). Failure properties of passive human aortic tissue. II-Biaxial tension tests, Journal of Biomechanics, 16(1): 31-44. E. A. Jafee, (1987). Cell Biology of Endothelial Cells, Human Pathology, 18(3): 234-239, G. A. Holzapfel, H. W. Weizsacker, (1998). Biomechanical behavior of the arterial wall and its numerical characterization, Computers in Biology and Medicine, 28(4): 377-392. G.L. Papageorgiou, N.B. Jones, (1987). Physical modeling of the arterial wall. Part2: Simulation of the non-linear elasticity of the arterial wall, Journal of Biomedical Engineering, 9(3): 216-221. I. Karsaj, J. D. Humphrey, (2012). A multilayered wall model of arterial growth and remodeling, Mechanics of Materials, 44(1): 110-119. L. A. Taber, (1998). A model for aortic growth based on fluid shear and fiber stresses, Journal of Biomechanical Engineering, 120(3): 348-354. M. A. Meyers, P.Y. Chen, A.Y. M. Lin, Y. Seki, (2008). Biological materials: Structure and mechanical properties, Progress in Materials Science, 53(1): 1-206. R. C. Batra, A. Bahrami, (2009). Inflation and eversion of functionally graded non-linear elastic incompressible circular cylinders, International Journal of Non-Linear Mechanics, 44(3): 311-323. R. H. Cox, (1978). Regional variation of series elasticity in canine arterial smooth muscles, American Journal of Physiology, 234(5): 542-551. R. L. Gleason,L. A. Taber,J. D. Humphrey, (2004). A 2-d model of flow induced alterations in the geometry, structure, and properties of carotid arteries, Journal of Biomechanical. Engineering, 126(3): 371 381. R. N. Vaishnav, J. T. Young, J. S. Janicki, D. J. Patel, (1972). Nonlinear anisotropic elastic properties of the canine aorta, Biophysical Journal, 12(8): 1008-1027. T. E. Carew, R. N. Vaishnav, D. J. Patel, (1968). Compressibility of the arterial wall, Circulation Research, 23(1): 61-68. W. Wan, L. Hansen,R. L. Gleason, (2010). A 3-D constrained mixture model for mechanically mediated vascular growth and remodeling, Biomechanic Model. Mechanobiology, 9(4): 403 419. W. W. von Maltzahn, (1982). Stresses and strains in the cone-shaped carotid sinus and their effects on baroreceptor functions, Journal of Biomechanics, 15(10): 757-765. W. W. von Maltzahn, D. Besdo, W. Wiemer, (1981). Elastic properties of arteries: A non-linear two layer cylindrical model, Journal of Biomechanics, 14(6): 389-397. W. W. von Maltzahn, R. G. Warriyar, W. F. Keitzer (1984). Experimental measurments of elastic properties of media and adventitia of bovine caroted arteries, Journal of Biomechanics, 17(11): 839-847. 67