1. Introduction
The problem of singularities in the theory of elasticity and in fracture mechanics is widely discussed in the related scientific literature [
1,
2,
3]. The singularity of solutions for stresses at the crack tip in the linear theory of elasticity excludes the use of traditional criteria for the strength of bodies with stress concentration. Moreover, the formally obtained singular solutions contradict not only the physical meaning, but also the postulates of the theory of elasticity. Such paradoxes still require an explanation [
3]. Gradient elasticity allows one to describe size effects [
4], and provides regularization of singular solutions of differential equations of elasticity theory [
5,
6]. In this regard, it would be quite natural to expect the development of gradient fracture mechanics. Nevertheless, so far only non-singular solutions have been constructed in gradient fracture mechanics for test problems corresponding to cracks of Mode III [
7].
For the general theory of gradient deformation, the variational model is described by the density of potential energy:
where
is the displacement vector,
εij the strain tensor,
Cijkm and
are classical and gradient moduli of elasticity respectively, having different dimensions. The force model is determined by the Cauchy stress tensor
and by the tensor of double stresses of the third rank
. Constitutive equations are obviously given by the following formulas:
In the general case, the physical properties of an isotropic gradient theory are described using seven physical constants, two of which are Lamé constants, and five others determine the tensor of gradient moduli of the sixth rank.
As a rule, one-parameter models are used to solve applied problems. One of the most popular is a model of the continuum in which the tensors of stresses and double stresses are determined by the formulas [
8,
9]:
and accordingly, the elastic modules have the form:
Here, , are the Lame coefficients, is the scale parameter, is the dilation, and is the Kronecker delta.
The equilibrium equations in such a model are written through the tensor of “total” stresses (see [
9]):
where
. Constitutive equations for such applied gradient models have the following form (see [
9]):
where
is the Laplace operator.
Note that static boundary conditions for the considered gradient model are specified through a linear combination of “classical” stresses and derivatives of double stresses . As a result, the boundary conditions cannot be written only through the tensor of stresses , which makes it difficult to obtain solutions. This probably explains why gradient solutions are constructed only for harmonic problems of cracks of Mode III.
Let us note that the problem of the regularization of the classical singular crack fields was still considered in [
10,
11,
12] where the attempts were made to construct non-singular solutions. Furthermore, in [
13] these studies were discussed because the obtained solutions did not satisfy the strain compatibility equations. Indeed, in the indicated works, it was possible to construct only approximate non-singular solutions, which is associated with incomplete consideration of the biharmonic components on the right side of the Helmholtz equations. In any case, in the indicated papers there are not the common analytical solutions to problems for cracks of Mode I. In the interesting work [
14] the search of solution for the theory of cracks comes down to integral equations.
The complexity of the obtained expressions and the method of constructing approximate solutions using the Chebyshev series did not allow the find of explicit analytical expressions that would show the regularity of the solution. The authors do not provide analytical solutions in an accurate study conducted in [
15] using micropolar elasticity, as well as in the fundamental work [
16] since, during the construction of an approximate solution, the authors are not able to ensure fully the continuation of the solution directly to the crack tip. In the present work, we are going to construct an exact non-singular solution for Mode I cracks using a new representation of the classical solution for generalized stresses through a complex-valued scalar potential.
The modified variant of the generalized elasticity theory was proposed by the authors of [
17,
18], where the defining relations are written for the generalized non-local stress tensor and the generalized non-local strain tensor. Note that the idea of the non-local theory of elasticity was proposed in the fundamental works [
19,
20]. Connection of non-local theory with gradient theories was established in the work [
21]. Theoretical aspects of non-local theories and their generalized variants and also their applications in the area of mechanics of nano-structures and nano devices discussed in the interesting works [
22,
23,
24].
We use the applied variant of model when a generalized non-local tensor function is introduced using the averaging operation of a local function on a given finite fragment and then the local function is expanded in a Taylor series on local coordinates on the fragment under consideration with retention of a finite number of terms in the series. As a result, after integrating the resulting finite Taylor series with respect to local coordinates over the considered fragment, the non-local generalized functions are constructed. In the proposed procedure, generalized functions are determined not only through the values of the local function at some point of the fragment, but also in aggregate through its derivatives of the second, fourth, etc. orders, depending on the number of terms held in the Taylor series. Here we can see a difference from the traditional differential calculus where the function or its derivatives are determined by their values at the point. Generalized functions in particular are determined by the corresponding local functions through the Helmholtz operator.
The construction of the solution in the model [
17,
18] is carried out in two stages. At the first stage, the traditional boundary value problem of the theory of elasticity is considered and non-local stresses and displacements are found. At the second stage, the obtained solutions are substituted on the right-hand side of the Helmholtz equation, which is solved with respect to local functions. Thus, to construct the local solutions it is first proposed to solve the classical elasticity problem for generalized stresses and displacements, and then find local stresses
and displacements
, by solving the corresponding Helmholtz equations:
We can see that in using this model for the mechanics of cracks the generalized non-local solutions formally coincide with the classical solution of the static boundary conditions and the local stress field must be found as common solutions of the non-homogeneous Helmholtz equation (2). The solution for local stresses and displacements has the form of a boundary layer.
We note that the procedure for constructing local solutions is similar to the construction of gradient solutions considered in [
25]. On the other hand, in the discussed model [
18] the constitutive equations are written on non-local generalized stresses and deformations
, which is quite justified, physically. As is evident, the constitutive equations of the discussed model are different from the constitutive equation (1). For the mechanics of cracks, the generalized non-local solutions formally coincide with the classical solution of the static boundary conditions, and the local stress field must be found as a common solution of the non-homogeneous Helmholtz equations. We assume that the local stresses are important for assessing the strength and that local stresses are involved in the strength criteria. These stresses are found after the definition of generalized stresses (“classical”)
from the equation
.
The purpose of this work is to show that, for the model under consideration, the local stress field in the vicinity of the crack tip is regular, i.e., it does not have a singularity. Non-singular solutions could be used as the basis for constructing and testing a new concept of stress concentration in crack mechanics proposed in [
26]. We belive that for brittle cracks the non-singular cracks dependences allow us to determine the important role for the scale parameter as a characteristic of fracture for brittle cracks along with tensile strength. In fact, a simular algorithm was implemented in [
26] for the case of finite cracks of normal separation in a plane strip. It was shown that the parameter
for a particular brittle material is a constant and can be considered as a critical parameter of fracture, giving with high accuracy a strength forecast for brittle materials.
We propose a new representation of the solution of the plane problem of the theory of elasticity through the complex-valued harmonic potentials included in the Papkovich–Neuber representation in a symmetric form, convenient for applications. A condition is determined under which the solution of plane problems of elasticity can be written through one harmonic potential. This condition is also fulfilled for problems of crack mechanics. As a result, we found the form of solution through one harmonic potential for displacement, and stresses that allow us to construct the non-singular solutions of mechanics of cracks.
2. New Representation of Solution through Harmonic Potential in Mechanics of Cracks
Let us consider the Papkovich–Neuber representation as a convenient analytical tool for describing the stress-strain state of an elastic body [
27]. This representation allows us to represent the displacements and the corresponding stresses, through two auxiliary harmonic potentials:
is the vector potential and another
is scalar,
:
here
is the displacement vector,
is the tensor of stresses,
, where
.
We propose a new representation of the solution of the plane problem of the theory of elasticity through the complex-valued harmonic potentials included in the Papkovich–Neuber representation in a symmetric form, convenient for applications. We show that these representations obtained directly from the representation of the general Papkovich–Neuber solution can be reduced to the well-known representations of Kolosov–Muskhelishvili for the plane problem of the theory of elasticity through two complex potentials. However, then we prove that under certain conditions a solution to a plane problem can be represented through one complex-valued harmonic potential. Such a representation is new and is convenient in constructing generalized solutions in crack mechanics.
Lemma 1. Papkovich–Neuber’s representation allows us to present a solution to the plane problem of the theory of elasticity through three complex potentials in the following form: Papkovich–Neuber’s representation leads to the following form of the common solution for complex potentials,,of the plane problem of the theory of elasticity through two analytic functionsand:where.
Proof. We consider the plane problem of the theory of elasticity and write the displacements and stresses in complex form, introducing in (3) complex potentials
,
instead of the corresponding vector expressions
,
and complex variables
instead variables
,
. Also, we use differentiation by complex coordinates
,
:
Let us define the first invariant of the stress tensor
and the complex potential
. Then the components of the displacement vector
,
and the stress tensor
,
,
can be expressed explicitly through the introduced complex potentials in the form of Equation (4). The Equation (3) determines the link between
,
,
and potentials
and
. Using representation (3) and taking into account the equalities
,
we can get
It is easy to establish also the following symmetric expressions for the stresses:
We take into account that the potential
is harmonic:
. Then Equation (7) leads to the following equalities for complex potentials
,
:
Furthermore, we take into account that for a plane statement the general complex solution of the Laplace equation
is expressed across two analytical functions
. On the other hand, the scalar potential
φ in (3) is real and can be written through one analytical function
. The following sequence of equalities shows the possibility of writing Equations (6) and (8) through potentials
F1,
F2 and
Φ1
here
and
, and
is the antiderivative function
. We took into account that the derivative of the analytic function with respect to the conjugate complex variable is equal to zero.
The stresses potentials are converted similarly:
Thus, Equations (6) and (8) are reduced to the form of Equation (5), which are defined with two analytical functions and . Note that for problems of plane stress state in Equations (3), (5), (6) and (8) the value should be replaced by the value . The lemma is proved. □
Theorem 1. Equation (5), obtained solely on the basis of the Papkovich–Neuber representation, fully corresponds to the Kolosov–Muskhelishvili representation (see [28]), in which stresses and displacements are expressed through two analytical functions using the formulas:where, χ are Kolosov–Muskhelishvili’s potentials
There is a relation between the functionsand under which the Papkovich–Neuber relation can be written through one harmonic potential, assuming that .
For problems of crack mechanics, the Pakovich–Neuber representation can always be written through one harmonic vector potential.
Proof. Comparing the right and left parts of Equations (5) and (9) it is easy to see that they are completely consistent with each other, if accept that
,
, or
Substitution of the functions and from Equation (11) into Equation (5) for complex displacements shows that the resulting equation exactly coincides with the expression for displacements. The first part of the theorem is proved.
To prove the second part of the theorem, we consider again the representation (3) and show that there is somewhat overdetermined, since between the functions
and
a relationship can be established, which does not change the stress–strain state. Indeed, we can introduce the following harmonic scalar function satisfying the equality:
As a result, if Equality (12) holds, then for any harmonic function it can be resolved with respect to the function and the Papkovich–Neuber representation (3) can be rewritten by redefining the vector potential with condition .
Let us briefly consider the cases when relation (12) is not satisfied. We represent the function in the form , where and are the angular coordinates, which are independent from the radial coordinate r. Equality (12) does not satisfied if . For negative values of the exponent, Equation (12) always has a solution, and without loss of generality the relation (3) can be written assuming condition . Therefore, for problems of crack mechanics, the Papkovich–Neuber representation can always be written in terms of a single harmonic vector potential. The theorem is proved. □
Thus, the general solution of problems of the plane theory of elasticity can be expressed not only using two analytical functions in the form of Equation (5), but also in equivalent form using one harmonic function
(i.e., through Papkovich–Neuber complex potential):
Let us consider the canonical singular problem of the theory of elasticity, which defines the problem of fracture mechanics for the plane problem:
and leads to the determination of the correct solution in stresses, tending to zero when
as
,
,
.
It is known that the correct solutions to this problem for a plane problem for normal separation crack and transversal shear, respectively are given through two Kolosov–Muskhelishvili potentials. This solution for a crack open displacement is written as follows (see [
1]):
Here , , , , are components of the displacement vector and stress tensor, respectively, are polar coordinates in the vicinity of the top of the crack is stress intensity factors.
Consider again the Papkovich–Neuber representation (see [
27]), which allows us to present a general solution to the problem of the theory of elasticity; that is, both displacement vector
and their corresponding stresses
on a surface element with a normal vector
, through only one harmonic vector potential
(in accordance with Theorem 1):
Here, the differential combination means the expression , , where is the stress tensor. Let us prove that the solution of the classical singular problem (for the example of a crack open displacement (15)–(17)) can be represented also through one harmonic complex-valued potential.
Theorem 2. An explicit expression exists for the harmonic potential for which the general representation of the solution of the plane Papkovich–Neuber problem (13), (14) is completely equivalent to the classical representation for the singular solution for a crack Mode I given by Equations (15)–(17).
Proof. Let us find an explicit expression of the harmonic potential
for a crack Mode I. We first consider the classical solution (13) for the displacement components
,
and transform them, distinguishing the coordinate
. We obtain:
Using the written equalities, we can get a complex-valued value of displacements
:
In the resulting expression, we select the derivative with respect to the coordinate
:
Assuming in (19) that
, we rewrite the complex value of displacements
in the form:
In order to reduce Equation (20) to the form (13) we must show that the last term in Equation (20) can always be compensated by the corresponding choice of the harmonic function. Instead of the harmonic function in Equation (20), we propose to consider another harmonic function obtained as the sum of the function and the harmonic function with some unknown coefficient , so far .
It is easy to verify that one can find a constant
such that Expression (20), in which instead of the function
stands
, coincides completely in appearance with Formula (13) for complex displacements, found in accordance with the Papkovich–Neuber representation. Indeed, after some transformations we find
. Now the statement of the theorem is verified by direct substitution. Thus, for a crack open displacement in the classical formulation, the field of elastic displacements and stresses in the vicinity of the crack tip is described by the Papkovich–Neuber representation (13) and (14) in complex form with a potential
having the following expression through complex variables
and
:
We have thus proved that the mathematical description of the crack is based on harmonic functions with a half-integer value of the exponent for the complex variable . The form obtained for representing the solution is based on one complex potential, which represents the solution in a form convenient for construction of the gradient solutions. The theorem is proved. □
The statement proved above is also valid for crack of mode II and crack of mode III. For example, for the classical singular solution mode II:
the following harmonic function of complex variables
can be found:
This harmonic function establishes a correspondence of the classical singular solution (22)–(24) with the solution given by the Papkovich–Neuber representation (13) and (14).
Finally, we consider the longitudinal shear crack (Mode III), which in the spatial theory of elasticity is described by the following relations for displacements and stresses:
It is easy to check that the stress–strain state for this case (26) is also completely described by the spatial version (3) of the Papkovich–Neuber representation with non-zero plane harmonic potential
and its compensating potential
, where:
Thus, strain–stress state for Mode I, II, and III singular cracks can be described on the basis of the Papkovich–Neuber formulas with the aid of one complex potential with a fractional degree. These forms are convenient for constructing a generalized gradient regular solution.
Remark 1. Note that harmonic functions establishing a correspondence between singular classical solutions of crack mechanics and solutions constructed using the Papkovich–Neuber representation are not analytic because, as is easily verifiable, the Cauchy relations are not satisfied for them.
In the future, our goal is to construct the local stresses and displacements associated with classical stresses and displacements through the Helmholtz equation using generalized elasticity. Consequently, the field of local stresses will describe the non-singular crack solutions.
3. Regular Gradient Solutions in Crack Mechanics
Let us construct a regular solution for crack Mode I, using an applied version of the non-local theory of elasticity. Using the algorithm discussed earlier in the Introduction, we can determine the local displacement and stress fields as general solutions of the inhomogeneous Helmholtz Equation (2), on the right-hand sides of which the well-known classical solutions are taken. We rewrite relation (2), introducing on the right-hand sides of the Helmholtz equations the classical solutions written in terms of complex displacements
and complex potentials
,
(see Equation (4):
,
. The generalized theory of elasticity allows finding the local stresses and displacements associated with classical stresses and displacements through the Helmholtz equation. The fields of local stresses describe the regular crack solutions as:
Here , , are the potentials of displacements and stresses expressed through the Papkovich–Neuber potential (21) by the formulas (13) and (14).
The general solution of Equation (28) is constructed as the sum of the general solution of the homogeneous Helmholtz equation and the particular solution of the inhomogeneous Helmholtz equation. It holds the following lemma, indicating the structure of a particular solution of the Helmholtz equation.
Lemma 2. Assume that function is the right-hand side of the Helmholtz equation, written with respect to a function. Assume also that function is a biharmonic function:
where . Then, a particular solution of the inhomogeneous Helmholtz equation has the form:
Proof. The lemma is proved by directly substituting Expression (30) into Equation (29). □
The regularity of the generalized solution of Equation (29) is ensured by the structure of a homogeneous solution and the special choice of arbitrary constants in this solution. Since the functions , , in Equation (28) and their particular solutions (see Equation (30)) contain singularities with only a half-integer value of the exponent in the degree of the complex variable , among the solutions of the homogeneous Helmholtz equation we are interested only in singular functions with the same value of the exponent in the degree of the complex variable .
It was shown in [
29,
30] that a system of such singular solutions can be constructed in an explicit analytical form using radial multipliers satisfying a special recurrence relation:
Thus, solutions of Equation (28) in the gradient theory of cracks can be constructed by explicitly taking into account Expressions (31) and (32).
A complete representation of the generalized solution (classical solutions) for a crack Mode I for generalized displacements
and generalized complex stress potentials
,
is given by Formulas (13), (14) and (21). As a result, taking into account the statement of Lemma 2 and Relations (31) and (32), we find:
Here , , are unknown coefficients that are selected from the compensation condition for singular terms in gradient potentials , , .
The specific structure of the solutions in the last terms of the above-stated equations, which are general solutions of the homogeneous Helmholtz equations, is found taking into account the order of the singularities in the corresponding particular solutions. So, in Expression (33), the first three terms determine a particular solution of the Helmholtz equation for generalized displacement . They include derivatives with respect to the argument up to the second order, which give rise to singularities of the form . Therefore, in order to compensate for the singularities in (33), we should use .
Furthermore, in Equality (34), the first term, which is a particular solution of the Helmholtz equation for the harmonic potential of stresses P, contains only the first derivative with respect to variables . Therefore, in the last term we should put , , .
Finally, since in Equation (35) a particular solution includes derivatives with respect to the argument up to the third order and contains singularities , and , then in the term representing the general solution of the homogeneous equation, we should put .
The constants
in the solution obtained in this way are explicitly found in such a way as to compensate for the singularities in the general solutions (33)–(35). We get:
Thus, Equalities (33)–(36) allow us to construct regular generalized solutions for complex displacements
and complex stress potentials P and
. The components of generalized displacements
,
, and generalized stresses
,
,
in the gradient theory of elasticity are calculated through the real and imaginary parts of complex potentials (33)–(35) according to the following formulas:
As a result, we obtain the following regular generalized solutions describing the non-singular stress–strain state of a normal separation crack in a gradient formulation:
where:
When deriving (38)–(40), we take into account the relations
,
in Expressions (39) and (40) written down taking into account (21). We also take into account that functions
,
and
are real functions. The coefficients
for the radial multipliers
,
and
are chosen so that in the expressions (38)–(40) there are no singular components. We can verify this, if we take into account the specific form of functions
and the asymptotic behavior of the exponential in functions
,
and
in Equations (38)–(40):
Using the polar coordinates
in Expressions (38)–(40), we obtain an explicit regular everywhere generalized gradient solution for the normal separation crack, written for the components of displacement:
and stress components:
Note that as the scale parameter tends to zero , the solutions (42)–(46) go into the classical singular solution.
The asymptotics of Expressions (44)–(46) up to terms of the order of
is calculated on the basis of Relation (41), where
. To use (41), we transform the expression in parentheses in (44)–(46) as follows:
As a result, we obtain the asymptotic behavior for stresses, which shows the absence of singularities at the origin:
Using the radial multipliers method and the Papkovich–Neuber representations for classical solutions of cracks Mode II and Mode III through the one harmonic potential, the local solutions for displacement and stresses of these modes of cracks were established:
and
Analytical solutions for displacements and stresses (42)–(46) and (47)–(52) supplement the results presented in the work [
29].
Non-singular solutions (44)–(46) are the basis for constructing and testing a new concept of stress concentration in crack mechanics proposed in [
26]. Indeed, we can see that for the problem considered above, the normal stresses
along the crack extension can be represented as
,
, where
is some parameter of the length. It is easy to verify, that the real value
is a function of
and is a positive bounded continuous function of the parameter
over the entire determination interval
that becomes zero when
. The maximum of this function is realized in some neighborhood of the point
.
Figure 1 shows the distribution of normalized stresses
in the vicinity of the crack tip alone the normalized distance from the crack tip (
). Obviously, the distribution of this function over the parameter
is a typical picture for stress concentration.
Therefore, to evaluate the strength, it is proposed to use the following procedure. First, the dependence of the maximums of relative stresses on the parameter is found. Then, using the known explicit Equation (45), the dependence of the maximum points of function on the parameter is constructed. This dependence is essentially a dependence of the stress concentration coefficient on a parameter . In the case when the ultimate stress level characterizing failure is determined by the ultimate stresses for brittle cracks, the constructed dependence allows one to determine the relative scale parameter , which along with tensile strength is characteristic of the fracture for brittle cracks.