Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2010 / Article

Research Article | Open Access

Volume 2010 |Article ID 460852 |

G. Ya. Popov, N. D. Vaysfel'd, "The Solution of Mitchell's Problem for the Elastic Infinite Cone with a Spherical Crack", Mathematical Problems in Engineering, vol. 2010, Article ID 460852, 20 pages, 2010.

The Solution of Mitchell's Problem for the Elastic Infinite Cone with a Spherical Crack

Academic Editor: Francesco Pellicano
Received10 Nov 2009
Accepted01 Apr 2010
Published22 Jun 2010


The new problem about the stress concentration around a spherical crack inside of an elastic cone is solved for the point tensile force enclosed to a cone's edge. The constructed discontinuous solutions of the equilibrium equations have allowed to express the displacements and stress in a cone through their jumps and the jumps of their normal derivatives across the crack's surface. The application of the integral transformation method under the generalized scheme has reduced the problem solving to the solving of the integrodifferential equation system with regard to the displacements' jumps. This system was solved approximately by the orthogonal polynomial method. The use of this method has allowed to take into consideration the order of the solution's singularities at the ends of an integral interval. The correlation between the crack's geometrical parameters, its distance from an edge, and the SIF values is established after the numerical analysis. The limit of the proposed method applicability is specified.

1. Introduction

The problems of the fracture mechanics and the non-destructive material testing demand the estimation of the stress intensity factor (SIF) around a crack located in an elastic body is important, as it is well known, because of it various applicatious in such engineering sciences as the fracture mechanics and the non-destructive material testing. The cracks’ researching in the unbounded elastic matrix can not reflect all the complicity of the crack’s phenomena in the real elastic body with boundaries. The mathematical complexity of the problems is caused by the necessity of the satisfaction to boundary conditions not only on the crack’s branches, but also on the boundaries of an elastic body. The topological form of a crack is also concerned with a number of factors complicating the problem’s solving. The big number of works, both in static and in dynamic statements, is devoted to the researching of the plane cracks with the different configurations of the contours [15]. The influence of the surface curvature, the variable curvature of a crack’s contour, the interference of the applied loading and geometrical parameters of a crack are investigated in these papers. In three-dimensional statement the nonplanar cracks in the unlimited elastic bodies were considered in [614]. In comparison with this, the number of the works investigating the nonplanar cracks, located in an elastic body with boundaries, is limited.

The problems on the cracks’ investigation in a finite elastic body are more often considered at the coincidence of a crack’s and bodies topology, because that allows to choose the same coordinate system for their description. So, the behavior of the penny-shaped cracks in the finite elastic cylinders is investigated in [1517]. The dynamic SIF around the spherical crack in a finite elastic shaft of the variable section is analyzed in [18]. The influence of an elastic cone boundaries on the SIF values around the spherical crack is shown in [19] in the assumption that at an cone’s edge the compressing force is applied. In the proposed paper the loading at the cone’s edge is the point tensile force, that essentially complicates the problem’s solving and allows to establish more general laws of the SIF correlation with the crack’s topology and the edge’s influence on its values.

2. Formulation of the Boundary Value Problem and the Discontinuous Solution Method

Let us consider the infinite elastic cone (Poison’s coefficient is , the shear module is , is the spherical coordinate system) at the vertex of which the concentrated force is applied (Figure 1).

On the cone surface the stress is given: The spherical crack is situated inside the cone at the distance from the vertex, its surface is described by the following relations: The crack’s branches are free from the stress

It is necessary to determine the SIF around the crack and to investigate the correlation between the SIF and the crack’s location and geometrical parameters.

The searched solution is constructed as the superposition of the continuous solution (in the assumption of the crack’s absence in the cone) and the discontinuous one (that one takes into consideration the existence of the crack). The first solution is marked by zero in the upper index and the second one by prim in the upper index: The continuous components are obtained in [20].

The method of the discontinuous solutions has been proved in Popov's works, and developed in the further in the papers [2, 12, 19]. The kernel of it lays in the construction of such solutions of Lame’s equations [21] (here the prime marks derivative with regard to the variable , the point represents derivative with regard to the variable ), which satisfy to these equations everywhere in the medium, except for the points of a defect. As a defect it can be understood both a crack, and an inclusion. At the transition across the defect’s surface the mechanical characteristics have the discontinuities of a continuity of the first kind. The jumps of the displacements and stress are assumed set. They are determined further in the problem’s statement and from the satisfaction of the boundary conditions. The constructed solutions allow to calculate the displacements and stress in any point of the medium with taking into consideration the discontinuity inside it.

We construct such solutions of (2.5) for a case of the crack defect of the spherical form. We will set jumps of the displacements and stress . To (2.5) the Mellin’s integral transformation is applied under the generalized scheme [2] (see Appendix A), and then the integral transformation with respect to the variable is used are the roots of the transcendental equation .

The equation system (2.5) is reduced by all of these transformations to the system of the linear algebraic equations in regard to the transformations of the displacements’ jumps. Let us resolve (2.8) and receive the transformations of the equilibrium equations’ discontinuous solutions which allow to express the displacements in any point of a medium through the jumps of the displacements and their normal derivatives across the crack’s surface.

To reduce the quantity of the unknown functions in the right parts (2.8) , we use conditions on the crack, having written them down in terms of the displacements: where One must satisfy to each condition on each crack’s branch. The integral transformation (2.7b) is applied to the first condition, and the integral transformation (2.7a) to the second one. After deducting the equalities when, one must obtain the following relations:

For the further shortening of the unknown functions in the right-hand parts of the equation system, one must satisfy to the first condition on the conical surface (2.1). They must write them in the displacements and apply Mellin’s integral transformation to it. After all these conversions the boundary condition will be as follows: Now one must substitute the equalities (2.11), (2.12) to the equation system (2.8), and solve it in regard to the unknown displacements’ transforms and : All coefficients are given in Appendix B. For obtaining of the discontinuous solution originals we apply the inverse Mellin’s transformation to (2.13).

One must use the residue theorem for the integral calculation with the following notes:()The roots of equation are simple. ()For the Jordan lemma satisfaction it is necessary to close the contour or if on the left (then one must take into consideration the simple roots and )—the case , or if on the right (then one must take into consideration the simple roots ) the case .

After the calculations we obtain the transformations of the displacements existing in the cone because of the crack’s presence (see Appendix C). To get the displacements’ originals, the inverse transformations (2.7a), (2.7b) should be applied to the solutions (2.13) with the equalities: Finally, the discontinuous solutions of the equilibrium equations are obtained

Let us note, that the specified procedure allows to receive the discontinuous solutions as for a crack case, so for a case of an inclusion, and under various conditions on a defect’s surface.

3. The Problem Reducing to the System of the Integral-Differential Equations

With the aim to satisfy to the conditions on the crack, one must finally demand the stress on the any of crack’s branches, for example, on the branch , will equal to zero. After substitution of the found discontinuous displacements in the boundary conditions when with the expression of solution in the form (2.4), the system of the integral equations with regard to the unknown displacement jumps is obtained: All taken notifications are written in Appendix D.

We must estimate the singularity order of the unknown functions in (3.1). As it is known, on the ends of the integration interval, the stress has the singularity of order: −. In the integral equations (3.1) the unknown functions , are the displacements’ jumps, and hence, with the formulas of the displacements and stress correlation, one can make the conclusion that these functions have on the ends of the integration interval the singularity of order . For the estimation of function singularity let us use the Williams’s method [22]. As it shown in Appendix E the searched order of the displacement’s singularity is . The further researching of the integral equations’ kernels consists in the obtaining of their asymptotical expressions for . Therefore, they need to know the asymptotics of the functions , and of the pairwise products, like , the last is possible because of well known formula [23], describing the asymptotic of Legandre’s functions with the large values of the order. Also one needs the formulas that were obtained in [24] for the asymptotics of the eigenvalues and of the functions’ norms: After using of all these relations, the following asymptotics of the pairwise products were derived and substituted in the kernels. The changing of the functions , was done with the asymptotical relations. It is necessary to sum the series in the kernels . The general scheme of this procedure is the following: the series is divided on the sum, in the second addend the general series term is changed by its asymptotical expression with the large values of . The next step is the adding and deduction of the sum . The initial series is written as the sum of the two addends . The series in this formula is the well-known one and could be found at the tables of the series.

After all these operations with the kernels , the table series where obtained [25]. The system of the integral equations (3.1) is reduced to the system of the two integro-differential equations (the derivation operator is exported from the integral sign with the aim of avoiding the divergent integrals in the kernels . are the regular kernels that were obtained by the scheme described earlier). In the last addends of both equations, the integration in parts was done. Finally the expression of the equation system (3.1) is

4. The Solving of the Integro-Differential Equation System

One must realize the standard scheme of the orthogonal polynomial method [2]. The spectral relation is needed for it (here —Chebyshev’s polynomial of the second order). The variable changing is done for the passing to the interval ,

In accordance with the singularity orders of the searched function , and the spectral relation (4.2), the unknown functions will be searched as the following expansions:

The infinite integral in the both integro-differential equations is changed by the finite one, and then the quadrature Sympson formula is applied with regard to the exponential character of function decreasing where are the quadrature Simpson formula coefficients, and are the unknown coefficients of the expansion.

The realization of the orthogonal polynomial method standard scheme leads to the system of the two-linear algebraic infinite equation system with regard to the following expansion coefficients: (In Appendix F one could see the linear algebraic equation system coefficients.) Taking into consideration the linearity of the SLAE solution, one must perform the unknown coefficients as the superposition of the unknown set of the constants:

Thus, it is necessary to solve the systems of the sort (4.5), differentiating one from another only by their right-hand parts: Each of these systems is solved by the reduction method. The argumentation of its availability could be done by the method proposed in [26].

After the solving of the equation system (4.5), the coefficients of the expansion (4.3) are obtained, and this would be the final step of the displacement jumps searching. For the estimation of the cone’s stress state, all that is left is to define the displacements along the conical surface.

5. The Calculation of the Displacements on the Conical Surface and SIF Values

It is needed to use the second condition (2.1) to find the displacement. For the condition’s satisfaction, one must demand that These conditions are written with the displacements’ expression The functions are defined by the discontinuous solutions’ kernels. Mellin’s transformation is applied to the relations (5.2)

One must show that the transformation , which is searched by the formula (5.3a) and by the formula (5.3b) is the same one. Really after the deduction from the right-hand part of (5.3a) of the right-hand part of (5.3b), zero will be the answer, so the coinciding of the left-hand parts is also proved. The fact of Mellin’s transformation equality in the both parts gives the result, that the originals of this transformations are also equal on the intervals and . That is why it is enough to solve the integral equation (5.2) on any of these intervals, for example, when . The solving is done by the method which was first used in [21]. It is important to use the fact that the singularity of the function in the vicinity of zero is equal to , as it is shown earlier. On the base of this fact the solution of the equation is constructed as the following expansion: where is Chebyshev-Lager polynomials. The series (5.4) are substituted to the equation. The obtained expression is multiplied by , and each member of it is integrated on the interval . As a result, the infinite system of linear equations is obtained: where This system is solved approximately by the reduction method. For the proof of the method’s convergence, the method proposed in [2] could be used.

The destruction criterion for the space case is Cherepanov’s formula [1] , where is the material’s constant.

In the stated problem, we have , and , are the coefficients at the stress and singularities correspondently: The stress in the formulas (5.7) is defined by the equalities

For the limit calculation in (5.7), it is necessary to use the continuation of the spectral relation (4.2) on the interval . The following equation is used [25]

Taking into consideration the variable change (4.2) and expressions (4.3), the result is obtained:

6. Numerical Results and Discussion

The dependence of the mode I SIF values and of the mode II SIF from the distance of a crack up to an edge is investigated at various values of a crack’s angle .

On Figure 2 values of

, calculated for a steel cone, which angle is , are resulted. The dotted curves correspond to the crack’s radius , and continuous to the radius . The analysis has shown that the angle of a crack, at which the mode I SIF reaches a maximum, almost twice is less than the value of the crack’s angle at which the mode II SIF reaches its one. The distance to the cone’s edge, as it is appreciable, influences the SIF absolute values, which are noticeably larger for the normal stress. The change of the crack’s distance up to an edge insignificantly influences the value of a crack’s angle at which maximum of SIF is reached.

The increase in loading essentially increases absolute SIF values though the cracks’ angle at which SIF reach the maxima vary insignificantly. The maximum of mode I SIF is reached by the smaller values of the crack’s angle than the values of the crack’s angle at which the mode II SIF reaches its peak.

Comparison has been lead and gave enough good concurrence with numerical results of SIF values calculation for the case of a spherical crack, located in the unlimited elastic medium at its stretching [13]. It is necessary to specify, that at ratio of the cone’s angle to the crack’s angle , smaller than 1, 2, calculations lost stability that testifies that the proposed approach to the problem solving in this case is inapplicable, and it is necessary to use, for example, a method of a small parameter.

On Figures 3 and 4 the curves show the dependences of and on the ratio correspondently.

From the analysis of the curves, notice that the crack’s distance from the cone’s surface essentially influences the SIF absolute values. Noticeably, already at values for the mode I SIF results coincide with a case of the infinite medium with a spherical crack, that is, the influence of the boundaries becomes insignificant. For the mode II SIF, the edge ceases to influence at .

7. Conclusions

() The approach proposed in the paper allows to solve the new axisymmetrical problem about stress concentration near the spherical crack located in an elastic cone at conditions of the first main elasticity problem on the cone’s surface at the point tensile force enclosed to edge.

() It is established that values of the mode I SIF are larger on absolute values than the mode II SIF. The maximal crack’s angle, at which the mode I SIF reaches the peak, almost is twice less than the crack’s angle at which reaches the maximum the mode II SIF.

() The crack’s distance to a cone’s surface renders more essential influence on size of absolute values of SIF than the distance on which the crack is located from an edge.

() The proposed approach to the problem solving is available when the ratio of a cone’s angle to the crack’s angle is not less than 1, 2.

() The method that was used in the paper allows to solve a similar problem for a defect of the inclusion type, and also to solve the more complicate problem for a compound elastic cone in which the crack settles down on a surface of an elastic constants changing, that is, interphase crack. Moreover, the proposed approach will allow to solve a problem for a case of the arbitrary oriented force enclosed to a cone’s edge for which it is necessary to construct the discontinuous solutions of the equilibrium equations for such case.


A. The Equilibrium Equations’ Representation in the Space of Mellin’s Integral Transformation Which was Applied by the Generalized Scheme

B. The Representation of the Transformations’ Coefficients of the Equilibrium Equations’ Discontinuous Solution

C. The Discontinuous Solutions in the Space of the Integral Transformation with Regard to the Variable

D. The Kernels and the Right-Hand Parts of the Integral Equations for the Unknown Jumps Searching

E. The Order of the Displacement’s Singularity

Guttmann’s representation of the equilibrium equation solutions [26] was used to obtain the order of the displacement’s singularity: where The operators and are defined by the equalities ,.

Correspondently to [22], the function is represented in the form , where is the searched order of the singularity. This representation is substituted in (E.2), and operator is applied to it. It leads to the solving of the differential equation: where are the unknown constants. The solution of this equation is

Taking into consideration the regularity of the solution of variable , one must demand that . It allows one to write: The final expression for the displacements will be acquired after the substitution of the solutions (E.6) into the formulas (E.1) when it follows that

The conditions of the problem (2.1) should be satisfied on the conical surface in order for the stress to be found from the known relations of the displacements and stress connections:

One must substitute the equalities (E.6) in the conditions (2.13) and pass to . With that, the homogenous system of equations with regard to the unknown constants is obtained. Its determinant should be equal to zero for its unique solution. It yields the transcendental equation for obtaining:

Equation (E.9) is solved numerically with MAPLE. By results of the roots’ analysis that root, which brings the strongest singularity in the solution, gets out. The searched value is, hence the searched order of the displacement’s singularity is .

F. The Coefficients of the Linear Algebraic Equation System with Regard to the Expansion Coefficients (4.5)


The research is supported by Ukrainian Department of Science and Education under project no. 0101U008297. The authors are grateful to Roman Bromblin for the help in the article’s text editing.


  1. H. S. Kit and M. V. Khaj, Method of Potentials in 3-D Thermoelasticity Problems for Solids with Cracks, Naukova Dumka, Kiev, Ukraine, 1989.
  2. G. Ya. Popov, Elastic Stress Concentration Around Dies, Cuts, Thin Inclusions and Supports, Nauka, Moscow, Russia, 1982.
  3. V. Z. Parton and V. G. Boriskovsky, Dynamic Fracture Mechanics: Stationary Cracks, vol. 1, Hemisphere, New York, NY, USA, 1989. View at: Zentralblatt MATH | MathSciNet
  4. J. Balas, J. Sladek, and V. Sladek, Stress Analysis by Boundary Element Method, Elsevier, Amsterdam, The Netherlands, 1989.
  5. Y. Shindo, “Axisymmetric elastodynamic response of a flat annular crack to normal impact waves,” Engineering Fracture Mechanics, vol. 19, no. 5, pp. 837–848, 1984. View at: Publisher Site | Google Scholar
  6. A. Boström and P. Olsson, “Scattering of elastic waves by non-planar cracks,” Wave Motion, vol. 9, no. 1, pp. 61–76, 1987. View at: Google Scholar
  7. B. Budiansky and J. R. Rice, “An integral equation for dynamic elastic response of an isolated 3-D crack,” Wave Motion, vol. 1, no. 3, pp. 187–192, 1979. View at: Google Scholar
  8. J. Dominguez and M. P. Ariza, “A direct traction BIE approach for three-dimensional crack problems,” Engineering Analysis with Boundary Elements, vol. 24, no. 10, pp. 727–738, 2000. View at: Publisher Site | Google Scholar
  9. R. V. Goldstein, “Three-dimensional elasticity problems related to cracks: exact solution by inversion transformation,” Theoretical and Applied Fracture Mechanics, vol. 5, no. 3, pp. 143–149, 1986. View at: Google Scholar
  10. V. V. Mikhas'kiv and I. O. Butrak, “Stress concentration around a spheroidal crack caused by a harmonic wave incident at an arbitrary angle,” International Applied Mechanics, vol. 42, no. 1, pp. 61–66, 2006. View at: Publisher Site | Google Scholar
  11. V. V. Mikhas'kiv, “Opening-function simulation of the three-dimensional nonstationary interaction of cracks in an elastic body,” International Applied Mechanics, vol. 37, no. 1, pp. 75–84, 2001. View at: Publisher Site | Google Scholar
  12. G. Ya. Popov, “Problems of stress concentration in the neighbourhood of a spherical defect,” Advances in Mechanics, vol. 15, no. 1-2, pp. 71–110, 1992. View at: Google Scholar | MathSciNet
  13. A. F. Ulitko, Vektor Expanding in the Space Elasticity Theory, Kiev, Ukraine, 2002.
  14. V. V. Panasuk, M. P. Savruk, and Z. T. Nazarchuk, Singular Integral Equations' Method in the Two Dimensional Elasticity Problems of Diffraction, Kiev, Ukraine, 1984.
  15. D. K. L. Tsang, S. O. Oyadiji, and A. Y. T. Leung, “Multiple penny-shaped cracks interaction in a finite body and their effect on stress intensity factor,” Engineering Fracture Mechanics, vol. 70, no. 15, pp. 2199–2214, 2003.