International Journal of Scientific & Engineering Research, Volume 6, Issue 4, April-2015 1672

ISSN 2229-5518

Effects of prestress on a hyperelastic, anisotropic and compressible tube

Dr Edouard Diouf

Laboratoire Mathématiques et Applications Université Assane Seck de Ziguinchor, BP 523 Ziguinchor, Sénégal

Abstract: We consider a hyperelastic, compressible, prestress and circular tube subjected to finite torsion and axial stretch. The mechanical response of the material is describes by an Ogden strain energy function with unidirectional reinforcement. The nonlinear equations governing the boundary problem are solved numerically by using the Runge-Kutta method. The effects of the prestress on the stress distributions are presented and permit to study a prototype vascular prosthesis of small diameter.

Keys-words: hyperelastic, compressible, anisotropic, prestress, small diameter vascular prosthesis

—————————— ——————————

1. Introduction

The following article describes valuable
techniques for fabricating devices deemed suitable for prosthetic applications in animals and humans. However, in vivo stability of the recommended constituent material is not completely demonstrated, and clinically important degenerative phenomena which may impact adversely on long-term safety are not addressed.
The design and fabrication of small diameter vascular prosthesis less than 6 mm which could assure a permanent functional service without the help of medicinal therapeutics, constitute a challenge subordinate to the necessity of well defining the fundamental characteristics of the ideal transplant [1]. The replacement of small-diameter blood vessels (≤ 6mm) is a clinical situation Frequently encountered by vascular and cardiothoracic surgeons [2].
Thus, in most cases, the use of the autologous saphenous vein, which is considered
to be excellent autologous graft material in
small diameter reconstructions [3], is the preferred solution. However, this vein is not available in one-third of patients, and the use of alternative
prostheses becomes necessary. Artificial materials are then proposed and are currently represented by commercially available polyester prostheses. Unfortunately, these biomaterials are not hemocompatible enough to allow a long- term
Patency rate (>5 years) and to be implanted
without the administration of
An adequate anticoagulotherapy [3].
The mimic the arterial materials [4], the elastic properties of biomaterials have to
be hyperelastic and anisotropic. Indeed, the mechanical properties of such biomaterials
are extremely important when selecting a material for the fabrication of vascular
The reader could consult the works of How
et al. [1] where the design requirement and the development of vascular grafts are investigated in details.
In this contribution, we propose a
compressible hyperelastic model to simulate the mechanical behavior of a tubular structure. This is a prototype of a small-diameter vascular graft composed of a silicone matrix embedded by fibers. Using experimental data [5], we show that the Ogden response can adequately represent the behavior of a silicone material. Furthermore, we show that such

IJSER © 2015

International Journal of Scientific & Engineering Research, Volume 6, Issue 4, April-2015 1673

ISSN 2229-5518

a fiber-reinforced silicone material offers a high mechanical similarity to natural arterial material. In the theoretical approach, we investigate the Ogden response augmented with unidirectional reinforcing characterized by a single additional constitutive parameter for strength of reinforcement. Then, the mechanical behavior of a prestressed tube subjected to a combined torsion and axial
For transversely isotropic compressible materials and in the particular case of four strain invariants, the Cauchy stress tensor σ is given by
(3) Where I d denotes the identity tensor and
stretch deformation is considered. The

Wk = ∂W

I k , (k = 1,2,3,4).
non-linear equations governing this
problem are established by the Runge-
Kutta method. The different deformations have been analyzed and are thought to result from certain local physiological conditions which should be supported by the stress distributions are then studied for different intensities of the prestress.

2. Method

Consider now an Ogden material reinforced in τ direction [6, 7]
where the scalar µo is the usual constant shear modulus for infinitesimal

2.1 Preliminaries

Consider a non-linearly elastic solid in its
undeformed state. In a cartesian coordinate system, let x = ( xi ) denote the coordinates of a particle in the undeformed state and
y = ( yi ) denote its coordinates in the
deformed state. We will note the components of the deformation gradient F, the left Cauchy-Green deformation tensor

B and the right Cauchy-Green deformation

tensor C, by
deformations from the undeformed state.
The constants (a1 , a2 , a3 ) are dimensionless parameters and the constant a4 and n represent respectively the density of
reinforcement and the fiber stiffness.
Also, from (3) and (4) yield the corresponding behavior law expressed as (5) Let us examine a state of plane strain

F = i , B = F
x j

F jk

, C = Fki


uniaxial stress parallel to the x1 − axis,
σ 11 = F ,σ 22 = 0, λ3 = 1.
Assuming that the material is transversely isotropic i.e. reinforced by fibers, the strain
When the fibers are parallel to the de
direction of loading
energy functional W will depend on four
(i.e.τ 1 = 1,τ 2 = 0,τ 3 = 0),
it may readily
strain invariants as follows
be shown that the applied stress F and the transverse stretch λ2 are related to the

where B* = (det(B))B −1 is the adjoint of B , J is the local volume change and τ = (τ i ) is the preferential direction vector in the undeformed state to the transverse
isotropic direction t = (ti ) in the deformed

principal stretch λ1

= Λ as follows

1 / 2



with λ2 =  2((a
+ 1) + (a
+ a )Λ2  .
state by t = F τ

I 4 .

1 1 2

IJSER © 2015

International Journal of Scientific & Engineering Research, Volume 6, Issue 4, April-2015 1674

ISSN 2229-5518

The engineering stress T is expressed as

T = Fλ2 .

Fig.1 Mechanical response (stress vs stretch) of silicone material in plane uniaxial stress

Figure 1 shows the stress-stretch relation for the isotropic case and the comparison with the experimental results obtained on homogeneous silicone specimen.
The fit parameters that have been found
(a1 = 0.25, a2 = −0,02, a3 = 3) and correspond to values which permit to approach measured data. We show clearly that, in tension, the isotropic model offers a close mechanical similarity to the silicone material.

2.2 Thick cylinder subjected to torsion and axial stretch

Consider a compressible and opened tube
defined by the angle Θ0 and made of a material des described by (4). Let us suppose that the tube undergoes two successive deformations. First including the closure and axial stretch of the tube which induce residual strains [8], and second including radial solicitations. The displacement is then prescribed by
where (R, Θ, Z ) and (r,θ , z) are respectively the reference and the deformed coordinates of a material particle in a cylindrical
system. Here Ψ is a twist angle per
unloaded length, Λ and λ are axial stretch ratios (respectively, for the first and the
second deformation). Let Ri and ri denote respectively the inner surfaces of the cylinder in the reference state and in the deformed state ( Ro and ro are the outer surfaces).
A routine calculation gives the physical
components of the F, B and C.

The four strain invariants I1 , I 2 , I 3 and I 4 are given by

where (τ R ,τ Θ ,τ Z ) are the components of the fibre direction in the undeformed configuration.
Using (5), it is easily follows that the stress
components with respect to cylindrical coordinates are

r = r (R), θ = 


 Θ + ΨΛZ ,
z = λΛZ ,
 Θ o

IJSER © 2015

International Journal of Scientific & Engineering Research, Volume 6, Issue 4, April-2015 1675

ISSN 2229-5518

Thus, the boundary-value problem consist of equations (11) with the boundary conditions, where the torsion angle Ψ and axial stretch Λ are given.
3. Results and discussion The boundary value problem is numerically solved by using the Runge- Kutta method of 4th order and a complementary iterative procedure, where
we take 100 points on the radial coordinate and check the circumferential stretch ration against the boundary condition until it is satisfied. Therefore, we consider the case where a tube made of silicone material is reinforced with axial fibers. The tube represents a prototype of small diameter vascular prosthesis [5] and is subjected to fixed twisting angle Ψ = 30°, axial stretch
Λ = 1,3, the density of reinforcement

a4 = 0,75

and different intensities of the
(10) From (5, 8) the equilibrium equations in
the absence of body forces are reduced to
prestress defined by the values of the parameter Θ 0 = {45°, 60°, 90°,120°}.

We note two step in Figure 2. First the circumferential stress decreases as

R ∈ R , R o Ri , then it grows rapidly
 
R ∈  R o Ri , R

as for Θ
= 45°.

dσ rr


+ 1 (σ

r rr

σ θθ

) = 0,

o o
 

Note also that σ θθ

is very large when R

dσ rθ


+ 2 σ = 0

r rθ

approaches R o .
dσ rz


+ σ rz


= 0.

This could be explained by the fact, the

residual stresses are much higher on the outer surface of the tube because of the opening angle.

Substituting (10) into (11.a), we obtain the
nonlinear differential equation

where f i , i = 1,2,3,4 are continuous for

R [Ri , Ro ].

In Figure 3, we note that the

circumferential stress increases with Θ o but all these circumferential stress are decreasing.

It should be noted the very large gap

between Figure 2 and Figure 3.

For an angle opening Θ o = 45°, σ θθ is

around 100 MPa while it is in the order of

Finally to determine r (R), we impose the
boundary conditions such as

2 MPa for Θ o

= 60°, 90°,120°.

r (Ri ) = Ri and r (Ro ) = Ro .

When R approach of

Ro for an angle

Θ o = 45°, the circumferential stress grows

IJSER © 2015

International Journal of Scientific & Engineering Research, Volume 6, Issue 4, April-2015 1676

ISSN 2229-5518

very rapidly, while it decreases for angles

Θ o = 60°, 90°,120°.

We can conclude in these cases, unless the opening angle is, the more σ θθ is great. There is more stress concentration for small angles than for high angles.

In Figure 4, we seek to establish the
influence of the parameter

n, the fiber

stiffness, in the circumferential stress distributions.

It appears that, for n = 2 or n = 3, and

Θ o = 45°, the circumferential stress decreases. Its maximum value is in the order of 1,8 MPa.
We note that for n = 4 in Figure 2, the stress σ θθ is very large so that it is slow for for n = 2 or n = 3. By comparing figures 2 and 4, we note that the fiber stiffness has a great influence in the circumferential stress distributions. In terms of maximum value, the function σ θθ is much larger than the
fiber stiffness is high.

Fig.2 Circumferential stress distributions vs radius for Θ 0 = 45°, n = 4

Fig.3 Circumferential stress distributions vs radius for different angles Θ 0 , n = 4
Fig.4 Circumferential stress distributions vs radius for Θ 0 = 45°, n = 2 and n = 3.


[1] How T.V,Guidoin R., young S.K.:
Engineering design of vascular prostheses. J.Eng.Med., Part H., p.61-71, (1992).
[2] Abbott, W.M.,Vignati,J.J.: Prosthetic
graft when are reasonable alternative?
Semin. Vasc.Surg.8, 235-245 (1995). [3] Greenwald S.E., Berry C.L.:
Improving vascular grafts: the important of the mechanical and

IJSER © 2015

International Journal of Scientific & Engineering Research, Volume 6, Issue 4, April-2015 1677

ISSN 2229-5518

haemodynamic properties. J. Pathol.
190,292-299, (2000).
[4] Fung Y.C.:Biomechanics,
Mechanical properties of living tissues, 2nd edition, Springer- Verlag.New-York, (1993).
[5] Cheref M.: Approche mécanique à la conception d’une prothèse vasculaire de petit diameter, PHD Thesis, université Paris 12 Val de Marne (1998).
[6] Diouf E., Exact solution of a problem of dynamic deformation and
nonlinear stability of a problem with a Blatz-Ko mateial, IJSER. Vol.3,
Isuue 2, (2012).
[7] Spencer AJ.M.: Continuum theory of
the mechanics of fibre-reinforced composites, Springer, New-York, (1984).
[8] Diouf E., Modélisation
hyperélastique, anisotrope, compressible et dynamique d’une structure tubulaire épaisse,PHD Thesis, université Paris 12 Val de Marne (2005).

IJSER © 2015