Fracture, mathematical problems of

The part of mechanics which studies the stress and strain fields around various defects in solid bodies (cracks, inclusions, dislocations, etc.). Several fracture criteria [a1] are then used to define a dangerous situation which may lead to a fracture of a certain material.

In terms of the mathematical apparatus involved, one may specify two categories of problems: two-dimensional problems and three-dimensional ones. The theory of functions of a complex variable [a3] (cf. also Functions of a complex variable, theory of) and the Fourier transform method [a7] have been used for solving two-dimensional problems (cf. also Two-dimensional problems in fracture mechanics), while three-dimensional problems were treated by the integral-transform method [a4] and the Green's function method [a2] (cf. also Green function).

One may classify within each category problems related to elastic (linear and non-linear), elasto-plastic, visco-elastic bodies, etc. [a1]. The linear elastic two-dimensional theory is the most developed [a3], including the process of crack propagation. Much less is known about three-dimensional problems. A different mathematical treatment is required for infinite, semi-infinite and finite bodies. Exact solutions are known at the moment (1988) only for the case of a flat circular or elliptic crack in an elastic space [a4]. A complete closed form solution was obtained [a2] for a penny-shaped crack in an elastic space, subjected to arbitrary internal tractions.

As an example, consider an elastic space weakened by a flat crack of general shape, subjected to an arbitrary normal traction $ p $. The governing integro-differential equation takes the form

$$ \tag{a1 } p ( N _ {0} ) = - \frac{1}{4 \pi ^ {2} H } \Delta {\int\limits \int\limits } _ { S } \frac{w ( N ) d S }{R ( N _ {0} , N ) } , $$

where $ S $ denotes the crack area, $ \Delta $ is the two-dimensional Laplace operator, $ w $ denotes the crack opening displacement, $ H = ( 1 - \nu ^ {2} ) / ( \pi E ) $( $ E $ is the elastic modulus and $ \nu $ is the Poisson coefficient), and $ R ( N _ {0} , N ) $ stands for the distance between two points $ N _ {0} $ and $ N $, with $ ( N _ {0} , N ) \in S $. An exact solution of the equation (a1) is known only for a circle, and in polar coordinates it is

$$ \tag{a2 } w ( \rho , \phi ) = \frac{2} \pi H \int\limits _ { 0 } ^ { {2 } \pi } \int\limits _ { 0 } ^ { a } \frac{p ( \rho _ {0} , \phi _ {0} ) }{R} \mathop{\rm tan} ^ {-} 1 \left ( \frac \eta {R} \right ) \rho _ {0} d \rho _ {0} d \phi _ {0} , $$

where $ a $ is the radius of the circle and

$$ \tag{a3 } R = [ \rho ^ {2} + \rho _ {0} ^ {2} - 2 \rho \rho _ {0} \cos ( \phi - \phi _ {0} ) ] ^ {1/2} , $$

$$ \eta = \frac{( a ^ {2} - \rho ^ {2} ) ^ {1/2} ( a ^ {2} - \rho _ {0} ^ {2} ) ^ {1/2} }{a} . $$

An approximate solution of (a1) was obtained in [a2] for a crack of general shape subjected to a uniform pressure $ p $. Let the crack boundary be defined in polar coordinates as $ \rho = a ( \phi ) $. The crack opening displacement is approximated by

$$ \tag{a4 } w = \frac \delta {a ( \phi ) } [ a ^ {2} ( \phi ) - \rho ^ {2} ] ^ {1/2} , $$

where the maximum crack displacement $ \delta $ is related to the applied pressure $ p $ by

$$ \tag{a5 } p = \frac{\delta B }{8 \pi H } ,\ \ B = \int\limits _ { 0 } ^ { {2 } \pi } \frac{d \phi }{a ( \phi ) } . $$

The coordinate origin is located at the crack centre, which is defined as the point where the value of $ B $ attains its minimum. The solution (a4)–(a5) has proven to be fairly accurate for a vast variety of crack shapes [a2].

Consideration of a general flat crack subjected to a shear loading $ \tau = \tau _ {zx} + i \tau _ {yz} $ leads to the governing integro-differential equation [a2]

$$ \tag{a6 } \tau ( N _ {0} ) = $$

$$ = \ - \frac{1}{2 \pi ^ {2} ( G _ {1} ^ {2} - G _ {2} ^ {2} ) } \left [ G _ {1} \Delta {\int\limits \int\limits } _ { S } \frac{u ( N) }{R ( N , N _ {0} ) } d S _ {N\right} . + $$

$$ + \left . G _ {2} \Lambda ^ {2} {\int\limits \int\limits } _ { S } \frac{\overline{u}\; ( N) }{R ( N , N _ {0} ) } d S _ {N} \right ] , $$

where the complex tangential displacement $ u = u _ {x} + i u _ {y} $ was introduced, and

$$ \tag{a7 } \Lambda = \frac \partial {\partial x } + i \frac \partial {\partial y } , $$

$$ G _ {1} = \frac{( 2 - \nu ) ( 1 + \nu ) }{\pi E } ,\ G _ {2} = \frac{\nu ( 1 + \nu ) }{\pi E } . $$

An exact closed form solution of (a6) was obtained recently for a circle, and is [a2]

$$ \tag{a8 } u ( \rho , \phi ) = $$

$$ \frac{G _ {1} } \pi \int\limits _ { 0 } ^ { {2 } \pi } \int\limits _ { 0 } ^ { a } \left [ \frac{1}{R} \mathop{\rm tan} ^ {-} 1 \frac \eta {R} - \frac{G _ {2} ^ {2} }{G _ {1} ^ {2} } \frac{( 3 - t bar ) \eta }{a ^ {2} ( 1 - \overline{t}\; ) ^ {2} } \right ] \tau ( \rho _ {0} , \phi _ {0} ) \rho _ {0} d \rho _ {0} d \phi _ {0} + $$

$$ + \frac{G _ {2} } \pi \int\limits _ { 0 } ^ { {2 } \pi } \int\limits _ { 0 } ^ { a } \left [ \frac{q}{R \overline{q}\; } \mathop{\rm tan} ^ {-} 1 \frac \eta {R} + \frac{\eta \left [ \left ( \frac{q}{\overline{q}\; } \right ) - t e ^ {2 i \phi _ {0} } \right ] }{a ^ {2} ( 1 - t ) ( 1 - \overline{t}\; ) } \right ] \tau bar ( \rho _ {0} , \phi _ {0} ) \rho _ {0} d \rho _ {0} d \phi _ {0} , $$

where $ R $ and $ \eta $ are defined by (a3), an overbar indicates the complex conjugate value, and

$$ t = \frac{\rho \rho _ {0} }{a ^ {2} } e ^ {i ( \phi - \phi _ {0} ) } ,\ \ q = \rho e ^ {i \phi } - \rho _ {0} e ^ {i \phi _ {0} } . $$

An approximate solution can be derived for a general flat crack subjected to a polynomial shear loading [a2]. In the simplest case of a uniform shear, the crack face displacements are approximated by the expression

$$ \tag{a9 } u ( \rho , \phi ) = u _ {0} \frac{[ a ^ {2} ( \phi ) - \rho ^ {2} ] ^ {1/2} }{a ( \phi ) } , $$

where $ u _ {0} $ is related to the applied loading $ \tau $ by

$$ \tag{a10 } u _ {0} = \ \frac{4 \pi ( G _ {1} ^ {2} - G _ {2} ^ {2} ) }{G _ {1} ^ {2} B ^ {2} - 9 G _ {2} ^ {2} B _ {2} ^ {2} } [ G _ {1} B \tau - 3 G _ {2} B _ {2} \overline \tau \; ] , $$

with $ B $ defined by (a5) and

$$ B _ {2} = \int\limits _ { 0 } ^ { {2 } \pi } \frac{e ^ {2 i \phi } d \phi }{a ( \phi ) } , $$

Formula (a10) is exact for an ellipse and is expected to perform fairly accurately for various crack shapes.

The stress intensity factors (SIF) are very important parameters in fracture mechanics. They are defined as follows:

$$ k _ {1} = \lim\limits _ {\rho \rightarrow a } \ \{ ( \rho - a ) ^ {1/2} \sigma _ {z} \} $$

for the stress intensity factor of mode I, and

$$ k _ {2} + i k _ {3} = \lim\limits _ {\rho \rightarrow a } \ \{ ( \rho - a ) ^ {1/2} \tau e ^ {- i \phi } \} $$

for the modes II and III loadings. In the case of a penny-shaped crack, one can write for an arbitrary distributed pressure [a1]

$$ k _ {1} = \frac{1}{\pi ^ {2} ( 2 a ) ^ {1/2} } \int\limits _ { 0 } ^ { {2 } \pi } \int\limits _ { 0 } ^ { a } \frac{( a ^ {2} - \rho _ {0} ^ {2} ) p ( \rho _ {0} , \phi _ {0} ) \rho _ {0} d \rho _ {0} d \phi _ {0} }{a ^ {2} + \rho _ {0} ^ {2} - 2 a \rho _ {0} \ \cos ( \phi - \phi _ {0} ) } . $$

In the general case of arbitrarily distributed shear loading, the stress intensity factors take the form [a2]

$$ k _ {2} + i k _ {3} = \ \int\limits _ { 0 } ^ { {2 } \pi } \int\limits _ { 0 } ^ { a } \frac{( a ^ {2} - \rho _ {0} ^ {2} ) ^ {1/2} }{\pi ^ {2} ( 2 a ) ^ {1/2} } \left [ \frac{e ^ {- i \phi } \tau ( \rho _ {0} , \phi _ {0} ) }{\rho _ {0} ^ {2} + a ^ {2} - 2 a \rho \cos ( \phi - \phi _ {0} ) } \right . + $$

$$ + \left . \frac{G _ {2} }{G _ {1} } \frac{( 3 a - \rho _ {0} e ^ {i ( \phi - \phi _ {0} ) } ) e ^ {i \phi } \overline \tau \; ( \rho _ {0} , \phi _ {0} ) }{a ( a - \rho _ {0} e ^ {i ( \phi - \phi _ {0} ) } ) ^ {2} } \right ] \rho _ {0} d \rho _ {0} d \phi _ {0} . $$

There is no general expression for the stress intensity factor for any shape other than circular. Some particular cases of loading of an elliptical crack were considered in [a4]. Various three-dimensional problems of thin inclusions are solved in [a5]. The mathematical description of crack propagation is still in the process of development.

References

[a1]	G.P. Cherepanov, "Mechanics of brittle fracture" , McGraw-Hill (1979) (Translated from Russian)
[a2]	V.I. Fabrikant, "Applications of potential theory in mechanics. Selection of new results" , Kluwer (1989)
[a3]	N.I. Muskhelishvili, "Some basic problems of the mathematical theory of elasticity" , Noordhoff (1953) (Translated from Russian)
[a4]	M.K. Kassir, G. Sih, "Three-dimensional crack problems" , Noordhoff (1975)
[a5]	V.V. Panasyuk, M.M. Stadnik, V.P. Silovanyuk, "Stress concentration in three-dimensional bodies with thin inclusions" , Kiev (1986) (In Russian)
[a6]	G.C. Sih, H. Liebowitz, "Mathematical theories of brittle fracture" , 2. Fracture , Acad. Press (1968)
[a7]	I.N. Sneddon, M. Lowengrub, "Crack problems in the classical theory of elasticity" , Wiley (1969)