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
ediouf@univ-zig.sn
—————————— ——————————
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
prosthesis.
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 http://www.ijser.org
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.
Consider now an Ogden material reinforced in τ direction [6, 7]
where the scalar µo is the usual constant shear modulus for infinitesimal
(4)
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
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
∂y F = i , B = F
∂x j
, C = Fki
(1)
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
(2)
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
− a
1 / 2
(6)
3
with λ2 = 2((a
+ 1) + (a
+ a )Λ2 .
state by t = F τ
I 4 .
1 1 2
IJSER © 2015 http://www.ijser.org
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 .
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.
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
(9)
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
(7)
IJSER © 2015 http://www.ijser.org
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
i
R ∈ R o − Ri , R
as for Θ
= 45°.
dσ rr
dr
+ 1 (σ
r rr
− σ θθ
) = 0,
o o
Note also that σ θθ
is very large when R
dσ rθ
dr
+ 2 σ = 0
r rθ
approaches R o .
(11)
dσ rz
dr
+ σ rz
r
= 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 http://www.ijser.org
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 http://www.ijser.org
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 http://www.ijser.org