Dynamic stress intensity factor of a finite crack based on a fractional differential model

Runtao Zhan1 , Zhaoxia Li2

1, 2Department of Engineering Mechanics, Southeast University, Jiangsu Key Laboratory of Engineering Mechanics, Sipailou 2, Nanjing, 210096, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 5, 2015, p. 128-133.
Accepted 24 July 2015; published 18 September 2015

Copyright © 2015 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

Fractional differential constitutive models are introduced for transient problem of Mode III finite length crack in a viscoelastic medium. The basic equations which govern the deformation behavior are converted to fractional wave-like equations. Integral transform method reduces the problem to Fredholm integral equation of second kind. Dynamics stress intensity factors of Mode III finite crack based on fractional differential constitutive are obtained by numerical solution of Fredholm integral equation.

Keywords: crack, fracture mechanics, dynamics, fractional calculus, wave.

1. Introduction

The increasing use of polymers and concrete demands an understanding of the response of cracked viscoelastic body to vibrational load. Rosakis et al. discovered cracks in a polymer, loaded under shear, propagate faster than the speed of shear waves [1]. Motivated by the experimental observations, it is desirable to make a study of cracks, loaded dynamically, in viscoelastic materials. Computation of the stress intensity factor (SIF) is one of the key tasks in fracture analysis. Dynamic stress intensity factors (DSIFs) of cracks imposed by dynamic load in viscoelastic materials have been explored by some researchers [2-6].

Previous work on DSIFs is mainly concerned with viscoelastic materials whose equation of motions are given in the form of wave equations, and constitutive equations are given in the form of various models of integer linear viscoelasticity, which are ‘‘integer differential” viscoelastic model and variations of them, however, they are not sufficiently accurate for more complex materials, i.e. some polymer. Recently the “fractional differential” models have been investigated for viscoelastic materials by a number of researchers [7, 8].

DSIFs of Mode III semi-infinite crack in a viscoelastic medium have been obtained by four parameters Bagley-Torvik fractional differential constitutive models [9]. This paper deal with the mode III finite length crack in a viscoelastic medium based on the fractional differential constitutive models.

2. Formulation of the problem

The basic equations of finite Mode III crack in viscoelastic medium consist of equation of motion as shown in Fig. 1.

Fig. 1. Finite Mode III crack in a viscoelastic medium

Finite Mode III crack in a viscoelastic medium
(1)
σ x z x + σ y z y = ρ 2 u z t 2 ,

Strain-displacement relations:

(2)
ε x z = 1 2 u z x ,       ε y z = 1 2 u z y ,

and constitutive equations defined by extending the one-dimension Bagley-Torvik fractional differential model to two-dimension [10, 11]:

(3)
1 + b 1 D α σ x z = 2 μ 1 + b 2 D β ε x z ,
(4)
1 + b 1 D α σ y z = 2 μ 1 + b 2 D β ε y z ,

where α, β (0 <α, β< 1) is the order of the fractional derivative and b1, b2 (b1, b2 0) are parameters. Obviously, when b1=b2= 0, Eqs. (3), (4) are equal to the linear elastic constitutive relationships.

There are various definitions of fractional integration and differentiation, such as Grunwald-Letnikov’s definition, Riemann-Liouville’s defintion, Caputo’s definition [12-14]. For the purpose of this paper, the Caputo’s definition of fractional integration and differentiation will be used, taking the advantage of Gaputo’s approach that the initial conditions for fractional differential equations with Caputo’s differentiation take on the same form as for integer-order differential equations [13, 14]:

(5)
D η f t = η f t η = 1 Γ n - η a t f n τ ( t - τ ) η + 1 - n d τ ,         n - 1 < R e η n ,       n N ,

where the parameter η is the order of fractional differentiation.

The Laplace transform of f(t) and its inverse formula are defined by:

(6)
L f t = f * p = 0 f t e - p t d t ,       f t = 1 2 π i c - i c + i f * p e p t d t .

The Laplace transform of the Caputo’s fractional derivative can be written as [14]:

(7)
L D η f t = p η L f t - k = 0 n - 1 f k 0 s η - k - 1 .

Taking the derivation of the Eqs. (3), (4) with respect to x, y, we have:

(8)
1 + b 1 D α σ x z x     = μ 1 + b 2 D β 2 u z x 2 ,
(9)
1 + b 1 D α σ y z y     = μ 1 + b 2 D β 2 u z y 2 .

The fractional wave-like equation can be obtained [9]:

(10)
1 + b 2 D β 2 u z = 1 c 2 2 u z t 2 + b 1 2 + α u z t 2 + a ,

where c=μ/ρ.

For the finite length crack, the boundary conditions are:

(11)
u z x , 0 , t = 0 ,       x > a ,       t > 0 ,
(12)
σ y z x , 0 , t = - τ H t ,       x < a ,       t > 0 ,
(13)
σ i j x , y , t = 0 ,         x 2 + y 2 ,         t > 0 ,

where τ is the amplitude of the anti-plane shear load, H(t) is Heaviside unit step function and crack length is 2a.

The initial conditions are:

(14)
u z ( x , y , 0 ) = 0 ,         u z x , y , t t t = 0 = 0 .

3. Dynamic stress intensity factors

Taking the Laplace transform of the Eq. (10) with respect to t, we have [9]:

(15)
2 u z * = ξ p c 2 u z * ,

where ξ2=1+b1pα/1+b2pβ,   ξ>0.

The Fourier transform and its inverse transform is:

(16)
f ¯ s = - f x e i s x d x ,       f x = 1 2 π - f ¯ s e - i s x d s .

Taking the Fourier transform of the Eq. (15) with respect to x, we have:

(17)
d 2 u ¯ z * ( s , y , p ) d y 2 = ξ p c 2 + s 2 u ¯ z * s , y , p .

The solution to the Eq. (17) by applying the boundary condition Eq. (13) is:

(18)
u ¯ z * s , y , p = C s , p e - γ y ,

where γ2=ξp/c2+s2, (γ> 0).

The Laplace transform to Eqs. (11) and (12) gives by using the symmetry:

(19)
u * z x , 0 , p = 0 ,       x > a ,
(20)
σ * y z x , 0 , p = - τ p ,       0 < x < a .

Taking the inverse Fourier cosine transform of the Eq. (18), we have [15]:

(21)
u z * x , y , p = 2 π 0 C s , p e - γ y cos s x d s .

Eqs. (19) and (21) can be applied to lead to:

(22)
0 C ( s , p ) c o s ( s x ) d s = 0 .

The Laplace transform to Eq. (9) with applying the Eq. (2) leads to:

(23)
1 + b 1 p α σ y z * = 1 + b 2 p β u z * y .

From Eq. (20) and (23), we find:

(24)
0 γ C ( s , p ) c o s s x d s = π τ 1 + b 1 p α 2 μ p 1 + b 2 p β .

Eqs. (20) and (23) form a pair of dual integral equations. The function C(s,p) is defined by:

(25)
C s , p = π τ 1 + b 1 p α a 2 2 μ p 1 + b 2 p β 0 1 ζ Φ 3 * ζ , p J 0 s a ζ d ζ .

The dual integral equations can be solved by a Fredholm integral equation of the second kind [15]:

(26)
Φ 3 * ζ , p + 0 1 K ζ , η , p Φ 3 * η , p d η = ζ ,

where:

(27)
ζ η 1 2 0 s 2 + ξ p a c 2 1 2 - s J 0 s ζ J 0 s η d s = K ζ , η , p ,

where J0 is a zero order Bessel function. Fredholm integral equation Eq. (26) is treated by the composite Simpson of rules to obtain Φ3*(1,p), the numerical solution to Eq. (26) can be found in Appendix Φ3*(1,p) is function of c/ξpa which can be seen in Fig. 2.

Fig. 2. Solution of Fredholm equation for finite Mode III crack

Solution of Fredholm equation for finite Mode III crack

The DSIF of cracked viscoelastic body can be found [15]:

(28)
K 3 t = L - 1 Φ 3 * 1 , p ξ p τ π a = L - 1 1 + b 2 p β 1 2 1 + b 1 p α 1 2 Φ 3 * 1 , p p τ π a .

From Eq. (1), the viscoelastic DSIF is equal to elastic DSIF when 1+b1pα=1+b2pβ.

The Miller and Guy numerical inverse Laplace transform is used for Eq. (1).

4. Numerical results

Let K3#(t)=K3/τπa. Numerical results for the dynamic SIF were obtained for 4 cases:

– Case a: b1=b2=8, β= 0.5, α= 0.25, 0.5, 0.75.

– Case b: b1=b2= 8, α= 0.5, β= 0.25, 0.5, 0.75.

– Case c: α=β= 0.5, b1= 8, b2= 4, 8, 12.

– Case d: α=β= 0.5, b2= 8, b1= 4, 8, 12.

Fig. 3. Viscoelastic DSIF time history for different cases

Viscoelastic DSIF time history for different cases

a) Case a

Viscoelastic DSIF time history for different cases

b) Case b

Viscoelastic DSIF time history for different cases

c) Case c

Viscoelastic DSIF time history for different cases

d) Case d

As shown in Fig. 3(a)-Fig. 3(d), DSIFs are proportional to the parameters β, b2, and inversely proportional to the parameters α, b1.

As shown in Fig. 3(a)-Fig. 3(b), The fractional differential order have strong effects on the viscoelastic DSIF after DSIF curves reach the peak of time history curve. When fractional order α>β, the curve is negative gradient, while when α<β, the curve is negative gradient after the peak of DSIF time history.

As shown in Fig. 3(c)-Fig. 3(d), parameters b1, b2 have little effect on the gradient of the DSIF curves after peak of DSIF curve.

5. Conclusion

Bagley-Torvik fractional differential model are introduced to obtain a fractional wave-like equation for finite Mode III crack. Bifurcation phenomenon of viscoelastic DSIF can be observed by choosing appropriate fractional differential order α, β. Parameters b1, b2 have some effect on amplitude of viscoelastic DSIF curve but little effect on the gradient of viscoelastic DSIF after the peak of the DSIF curve.

References

  1. Rosakis A., Samudrala O., Coker D. Cracks faster than the shear wave speed. Science, Vol. 284, 1999, p. 1337-1340 [Search CrossRef]
  2. Willis J. Crack propagation in viscoelastic media. Journal of the Mechanics and Physics of Solids, Vol. 15, 1967, p. 229-240. [Search CrossRef]
  3. Walton J. R. Dynamic energy release rate for a steadily propagating antiplane shear crack in a linearly viscoelastic body. Journal of Applied Mechanics, Transactions ASME, Vol. 54, 1987, p. 635-641. [Search CrossRef]
  4. Walton J. R. Dynamic energy release rate for a steadily propagating mode I crack in an infinite, linearly viscoelastic body. Journal of Applied Mechanics, Transactions ASME, Vol. 57, 1990, p. 343-353. [Search CrossRef]
  5. Walton J. R. On a dynamically accelerating dugdale-zone in elastic and viscoelastic material. Journal of the Mechanics and Physics of Solids, Vol. 44, 1996, p. 1353-1370. [Search CrossRef]
  6. Khalmanova D., Walton J. R. Dynamic steady-state crack propagation in a transversely isotropic viscoelastic body. International Journal of Fracture, Vol. 131, 2005, p. 107-129. [Search CrossRef]
  7. Rossikhin Y., Shitikova M. Applications of fractional calculus to dynamic problems of linear and nonlinear hereditary mechanics of solids. Applied Mechanics Review, Vol. 50, 1997, p. 15-67. [Search CrossRef]
  8. Rossikhin Y., Shitikova M. Application of fractional calculus for dynamic problems of solid mechanics: novel trends and recent results. Applied Mechanics Review, Vol. 63, 2010, p. 1-52. [Search CrossRef]
  9. Zhan R. T., Li Z. X., Wang L. A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials. Acta Mechanica Sinica, Vol. 30, 2014, p. 403-409. [Search CrossRef]
  10. Bagley R. L., Torvik P. J. Fractional calculus – a different approach to the analysis of viscoelastically damped structures. AIAA Journal, Vol. 21, 1983, p. 741-748. [Search CrossRef]
  11. Bagley R. L., Torvik P. J. Fractional calculus in the transient analysis of viscoelastically damped structures. AIAA Journal, Vol. 23, 1985, p. 918-925. [Search CrossRef]
  12. Oldham K. B., Spanier J. Fractional Calculus: Theory and Applications, Differentiation and Integration to Arbitrary Order. Academic Press, New York, London, 1974. [Search CrossRef]
  13. Caputo M. Linear models of dissipation whose Q is almost frequency independent – II. Geophysical Journal International, Vol. 13, 1967, p. 529-539. [Search CrossRef]
  14. Podlubny I. Fractional Differential Equation. Academic Press, New York, 1999. [Search CrossRef]
  15. Sih G. C. Elastodynamic Crack Problems. Leyden, Noordhoff, 1977. [Search CrossRef]
  16. Miller M. K., Guy J. W. T. Numerical inversion of the Laplace transform by use of Jacobi polynomials. SIAM Journal on Numerical Analysis, Vol. 3, 1966, p. 624-635. [Search CrossRef]