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

Authors: G. Ya. Popov; N. D. Vaysfel'd
Journal: Mathematical Problems in Engineering
 (2010)

---

## Abstract

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.


---

## Body

## 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 [1–5]. 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 [6–14]. 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 [15–17]. 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 cone0<r<∞,0≤θ≤ω,-π≤φ≤π (Poison’s coefficient is μ, the shear module is G,r,θ,φ is the spherical coordinate system)  at the vertex of which the concentrated force P is applied (Figure 1).Figure 1On the cone surface the stress is given:(2.1)τrθ(r,θ)|θ=ω=0,σθ(r,θ)|θ=ω=0.
The spherical crack is situated inside the cone at the distance from the vertex, its surface is described by the following relations:
(2.2)r=R,0≤θ≤ω0,-π≤φ≤π(ω0<ω).
The crack’s branches are free from the stress(2.3)τrθ(r,θ)|r=R±0=0,σr(r,θ)|r=R±0=0.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:(2.4)ur(r,θ)=ur0(r,θ)+ur1(r,θ),σθ(r,θ)=σθ0(r,θ)+σθ1(r,θ),uθ(r,θ)=uθ0(r,θ)+uθ1(r,θ),τrθ(r,θ)=τrθ0(r,θ)+τrθ1(r,θ).
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](2.5)(r2u')'-2u-μ**μ*(νsinθ)•sinθ+(sinθu•)•sinθ1μ*+r(ν′sinθ)•sinθμ0μ*=0,(r2ν')'+μ*((ν•sinθ)•sinθ-νsin2θ)+μ0ru'•+2μ*u•=0,μ*=1+μ0,μ**=μ0+2,μ0=(1-2μ)-1,u=u(r,θ)=ur1(r,θ),ν=ν(r,θ)=uθ1(r,θ).
(here the prime marks derivative with regard to the variable r, 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 χ(θ)=〈u(R,θ)〉,ψ(θ)=〈ν(R,θ)〉,〈f(R,θ)〉=f(R-0,θ)-f(R+0,θ).  To (2.5)  the  Mellin’s integral transformation is applied under the generalized scheme [2](2.6)fs(θ)=(∫0R-0+∫R+0∞)f(r,θ)rs-1dr,f(r,θ)=12πi∫γ-i∞γ+i∞fs(θ)r-sds,
(see Appendix A), and then the integral transformation with respect to the variable θ is used(2.7a)uks=∫0ω0us(θ)Pνk0(cosθ)sinθdθ,us(θ)=∑k=0∞uskPνk0(cosθ)∥Pνk0(cosθ)∥2,(2.7b)νks=∫0ω0νs(θ)Pνk1(cosθ)sinθdθ,νs(θ)=∑k=0∞νskPνk1(cosθ)∥Pνk1(cosθ)∥2νk are the roots of the transcendental equation Pνk1(cosω)=0,k=0,1,2,….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.(2.8)(s(s-1)-2)usk+1μ*νk(νk+1)usk+μ**μ*νsk+μ0μ*sνsk=Rs(s-1)χk-Rs+1〈uk′(R)〉-sinωμ*us•(ω)Pνk0(cosω)+μ**μ*νs(ω)sinωPνk0(cosω)+μ0μ*ssinωPνk0(cosω)νs(ω)+μ0μ**Rsψk,s(s-1)νsk-μ*νk(νk+1)νsk-μ*sνk(νk+1)usk+2μ*usk=Rs(s-1)ψk-Rs+1〈νk′(R)〉+μ*sinω(Pνk1(cosω))•νs(ω)-μ0Rsνk(νk+1)χk.
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) χk,ψk,νs(ω),us•(ω),〈uk′(R)〉,〈νk′(R)〉, we use conditions on the crack, having written them down in terms of the displacements:(2.9)12r[r2(νr)′+u•]|r=R±0=0,[μμ0θ+u']|r=R±0=0,
where
(2.10)θ=(r2u)'r2+(vsinθ)•rsinθ.
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 whenr=R±0, one must obtain the following relations: (2.11)〈uk′(R)〉=μR(1-μ)ψk-3μR(1-μ)χk,〈νk′(R)〉=ψkR-νk(νk+1)Rχk.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:(2.12)us•(ω)=(s+1)νs(ω).
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   usk and νsk:(2.13)usk=χkRsα(s,k)+ψkRsβ(s,k)+νs(ω)γ(s,k)Δsk,νsk=χkRsq(s,k)+ψkRsl(s,k)+νs(ω)p(s,k)Δsk
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:(1)
The roots of equationΔsk=0,s=sj,andj=1,4¯ are simple.(2)
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 rootss=s1 and s=s2)—the case r<R, or if on the right  (then one must take into consideration the simple roots s=s3,s4) the case r>R.After the calculations we obtain the transformations of the displacements existing in the cone because of the crack’s presence (see AppendixC). To get the displacements’ originals, the inverse transformations (2.7a), (2.7b) should be applied to the solutions (2.13) with the equalities: (2.14)χk=∫0ω0χ(η)Pνk0(cosη)sinηdη,ψk=∫0ω0ψ(η)Pνk1(cosη)sinηdη.
Finally, the discontinuous solutions of the equilibrium equations are obtained (2.15)u(r,θ)=∫0ω0χ(η)∑k=0∞{∑j=12[α(sj,k)Δj,k](Rr)sjPνk0(cosθ)∥Pνk0(cosθ)∥2Pνk0(cosη)}dη+∫0ω0ψ(η)∑k=0∞{∑j=12[β(sj,k)Δj,k](Rr)sjPνk0(cosθ)∥Pνk0(cosθ)∥Pνk1(cosη)}dη+∑k=0∞Fk(r)Pνk0(cosθ)∥Pνk0(cosθ)∥2,r<R,u(r,θ)=∫oω0χ(η)∑k=0∞{∑j=34[α(sj,k)Δj,k](Rr)sjPνk0(cosθ)∥Pνk0(cosθ)∥2Pνk0(cosη)}dη+∫0ω0ψ(η)∑k=0∞{∑j=34[β(sj,k)Δj,k](Rr)sjPνk0(cosθ)∥Pνk0(cosθ)∥2Pνk1(cosη)}dη+∑k=0∞Fk(r)Pνk0(cosθ)∥Pνk0(cosθ)∥2,r>R,ν(r,θ)=∫0ω0χ(η)∑k=0∞{∑j=12[g(sj,k)Δj,k](Rr)sjPνk1(cosθ)∥Pνk1(cosθ)∥2Pνk0(cosη)}dη+∫0ω0ψ(η)∑k=0∞{∑j=12[l(sj,k)Δj,k](Rr)sjPνk1(cosθ)∥Pνk1(cosθ)∥2Pνk1(cosη)}dη+∑k=0∞Gk(r)Pνk1(cosθ)∥Pνk1(cosθ)∥2,r<R,ν(r,θ)=∫0ω0χ(η)∑k=0∞{∑j=34[g(sj,k)Δj,k](Rr)sjPνk1(cosθ)∥Pνk1(cosθ)∥2Pνk0(cosη)}dη+∫0ω0ψ(η)∑k=0∞{∑j=34[l(sj,k)Δj,k](Rr)sjPνk1(cosθ)∥Pνk1(cosθ)∥2Pνk1(cosη)}dη+∑k=0∞Gk(r)Pνk1(cosθ)∥Pνk1(cosθ)∥2,r>R,Fk(r)=∫0∞ν(ψ,ω)gk(rξ)dξξ,Gk(r)=∫0∞ν(ξ,ω)yk(rξ)dξξ.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 branchr=R-0, will equal to zero. After substitution of the found discontinuous displacements in the boundary conditions when r=R-0 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: (3.1)∫0ω0χ(η)F1(θ,η)dη+∫0ω0ψ(η)F2(θ,η)+∫0ω0ν(ξ,ω)α1(ξ,θ)dξ+∫0∞ν'(ξ,ω)α2(ξ,θ)dξ=-τ∘(θ,R),∫0ω0χ(η)F3(θ,η)dη+∫0ω0ψ(η)F4(θ,η)dη+∫0∈ν(ξ,ω)α3(ξ,θ)dξ+∫0∞ν'(ξ,ω)α4(ξ,θ)dξ=-σ20(θ,R).
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, η=ω0 the stress has the singularity of order: −(1/2). 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 (1/2). For the estimation of function ν(r,ω) singularity let us use the Williams’s  method [22]. As it shown in Appendix E the searched order of the displacement’s singularity  is ν(r,ω)~rλ*-1.  The further researching of the integral equations’ kernels Fj(θ,η),j=1,4¯ consists in the obtaining of their asymptotical  expressions for k→∞. Therefore, they need to know the asymptotics of the functions α(sj,k),β(sj,k),l(sj,k),q(sj,k),yk(ξ),gk(ξ), and of the pairwise products, like  Pvk0(cosθ)Pvk0(cosη), Pvk0(cosθ)Pvk1(cosη),Pvk1(cosθ)Pvk1(cosη); 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: (3.2)νk~kπωwithk→∞,∥Pνk{q}(cosθ)∥2~k,q=0,1withk→∞.
After using of all these relations,  the following asymptotics of the pairwise products were derived and substituted in the kernelsFj(θ,η),j=1,4¯. The changing of the functions α(sj,k),β(sj,k),l(sj,k),q(sj,k),yk(ξ),gk(ξ),  was done with the asymptotical relations. It is necessary to sum the series in the kernels Fj(θ,η),j=1,4¯. The general scheme of this procedure is the following:  the series ∑k=1∞ak(θ,η) is divided on the sum∑k=1Nak(θ,η)+∑k=N+1∞ak(θ,η), in the second addend the general series term is changed by its asymptotical expression   ãk(θ,η) with the large values of  k. The next step is the adding and deduction of the sum ∑k=1Nãk(θ,η). The initial series is written as the sum of the two addends ∑k=1∞ak(θ,η)=∑k=1∞ãk(θ,η)+∑k=1N(ak(θ,η)-ãk(θ,η)). 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 kernelsFj(θ,η),j=1,4¯, 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  Fj(θ,η),j=1,4¯. F̂j(θ,η) 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 (3.3)d2dθ2[∫0ω0(χ(η)+ψ(η))ln1|η-θ|dη]+∫0ω0χ(η)F̃1(θ,η)dη+∫0ω0ψ(η)F̃2(θ,η)dη+∫0∞ν(ξ,ω)B1(ξ,θ)dξ=-τ∘(θ),d2dθ2[∫0ω0(χ(η)+ψ(η))ln1|η-θ|dη]+∫0ω0χ(η)F⌢3(θ,η)dη+∫0ω0ψ(η)F̃3(θ,η)dη+∫0∞ν(ξ,ω)B2(ξ,θ)dξ=-σr∘(θ).


## 4. The Solving of the Integro-Differential Equation System

One must realize the standard scheme of the orthogonal polynomial method [2]. The spectral relation (4.1)d2dx2∫-11ln1|x-y|1-y2Un(y)dy=-π(n+1)Un(x).
is needed for it (here Un(x)—Chebyshev’s polynomial of the second order). The variable changing is done for the passing to the interval (0,1)—x=2η-1,y=2ξ-1,(4.2)d2dη2∫01ln1|ξ-η|1-ξ2Un(2ξ-1)dξ=-π4(n+1)Un(2η-1).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:(4.3)χ(η)=∑k=0∞χkη-η2Un(2η-1),ψ(η)=∑k=0∞ψkη-η2Un(2η-1).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 functionBj(ξ,θ) decreasing  (4.4)∫0TV(ξ)Bj(ξ,θ)dξ=∑n=1NVnAnj-Bnj(θ),j=1,2,
where Anj are the quadrature Simpson formula coefficients, and Vn 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:(4.5)χl+ψl+∑k=1∞χkF̃kl1+∑k=1∞ψkF̃kl2=fl1+∑n=1NVnAn1Bnl1,χl+ψl+∑k=1∞χkF̃kl3+∑k=1∞ψkF̃kl4=fl2+∑n=1NVnAn2Bnl2.
(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 χk,ψk(k=1,∞¯) as the superposition of the N+1 unknown set of the constants:   (4.6)χk=∑l=1N+1χkl,ψk=∑l=1N+1ψkl,k=1,∞¯.Thus, it is necessary to solve theN+1 systems of the sort (4.5), differentiating one from another only by their right-hand parts:(4.7)(fl1,fl2),(B1l1A11,B1l2A12),…,(BNl1AN1,BNl2AN2).
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 Surfacev=(r,ω) and SIF Values

It is needed to use the second condition (2.1) to find the displacementv(r,ω). For the condition’s satisfaction, one must demand that(5.1)σθ(r,θ)|θ=ω=0whenr<R,σθ(r,θ)|θ=ω=0r>R.
These conditions are written with the displacements’ expression (5.2)ν(r,ω)-∫0∞ν(ξ,ω)G(ξ,r)dξ=-∫0σ0χ(η)φ1(η,r)dη-∫0ω0ψ(η)φ2(η,r)dη,r<R,ν(r,ω)-∫0∞ν(ξ,ω)G(ξ,r)dξ=-∫0ω0χ(η)φ3(η,r)dη-∫0ω0ψ(η)φ4(η,r)dη,r>R,G(ξ,r)=∑k=0∞[gk(rξ)+gk′(rξ)]Pνk1(cosω)∥Pνk1(cosω)∥2+∑k=0∞yk(rξ)(Pνk1(cosω))•∥Pνk1(cosω)∥2.
The functions φj(η,r) are defined by the discontinuous solutions’ kernels. Mellin’s transformation is applied to the relations (5.2) (5.3a)νs(ω)=-∫0ω0∫0∞φ1(η,r)rs-1dr1-G(s)χ(η)dηω0-∫0ω0ψ(η)∫0∞φ2(η,r)rs-1drdη,(5.3b)νs(ω)=-∫0ω0χ(η)∫0∞φ3(η,r)rs-1dr-∫0ω0ψ(η)∫0∞φ4(η,r)rs-1drdη1-G(s).One must show that the transformationνs(ω), 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 νs(ω) equality in the both parts gives the result, that the originals of this transformations are also equal on the intervals 0<r<R and R<r<+∞. That is why it is enough to solve the integral equation (5.2) on any of these intervals, for example, when r<R.  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  ν(r,ω) 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: (5.4)ν(r,ω)=χ̃(r)=∑n=0∞χ̅ne̅rr-λ*Ln(-λ*)(2r),
where Ln(-α)(r) is Chebyshev-Lager polynomials. The series (5.4) are substituted to the equation. The obtained expression is multiplied by r2-λLm(2-λ)(2r)e̅r(m=0,1,2,…), and each member of it is integrated on the interval  (0;+∞). As a result, the infinite system of linear equations is obtained: (5.5)Xm+∑n=0∞AmnXn=Cm(λ),m,n=0,1,2,…,
where(5.6)Amn=∫0∞r2-λe̅rLm(2-λ)(2r)K(r,λ)dr,K(r,λ)=12π∫-∞∞e̅ξLn(-λ)(2ξ)e̅ξG(r,ξ)dξ,Cm(λ)=∫0∞r2-λe-rLm(2-λ)(2r)F(r)dr.
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] KI+KII+KIII=C, where C is the material’s constant.In the stated problem, we haveKIII≡0, and KI, KII are the coefficients at the stress σr and  τrθ singularities correspondently: (5.7)KI=limθ→ω0+02πθ-ω0τrθ(R,θ),KII=limθ→σ0+02πθ-ω0σr(R,θ).
The stress in the formulas (5.7) is defined by the equalities(5.8)τrθ(R,θ)=d2dθ2∫0ω0χ(η)ln1|η-θ|dη+d2dθ2∫0ω0ψ(η)ln1|η-θ|dη+R1(θ)+τ0(θ,R),σr(R,θ)=(μμ0+1)d2dθ2∫0ω0ln1|η-θ|dη+(μμ0+1)∫0ω0ψ(η)ln1|η-θ|dη+R2(θ)+σr0(θ,R).For the limit calculation in (5.7), it is necessary to use the continuation of the spectral relation (4.2) on the interval |η|>1. The following equation is used [25](5.9)d2dx21π∫-11ln1|x-s|1-s2Vm(s)ds=(m+1)22m+2(x-1)m+2×[F(32+m,m+2;32;x+1x-1)-m+121-x-1-xΓ(32+m,m+1;12;x+1x-1)].Taking into consideration the variable change (4.2) and expressions (4.3), the result is obtained:(5.10)KI=π2(∑m=1∞(-1)m+1mXm+∑m=1∞(-1)m+1mψm),KII=π2(μμ0+1μ*)(∑m=1∞(-1)m+1mXm+∑m=1∞(-1)m+1mψm).


## 6. Numerical Results and Discussion

The dependence of the mode I SIF valuesKI and of the mode II SIF KII from the distance of a crack up to an edge is investigated at various values of a crack’s angle ω0.On Figure2 values of K̃I=(KI/PR),K̃II=(KII/PR)Figure 2, calculated for a steel cone, which angle isω=75°, are resulted.  The dotted curves correspond to the crack’s radius R=R1, and continuous to the radius R=2R1. 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 valuesKI∞,KII∞ 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 l=ω/ω0, 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 Figures3 and 4 the curves show the dependences of KI*=KI/KI∞ and  KII*=KII/KII∞on the ratiol=ω/ω0 correspondently.Figure 3Figure 4From 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 valuesl≥8 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 l≥6.


## 7. Conclusions

(1) 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.(2) 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.(3) 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.(4) 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.(5) 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.



---

*Source: 460852-2010-06-22.xml*