1. Introduction
It is unanimously accepted that the non-classical theories of elasticity have as their main goal to eliminate experimental aspects that cannot be described by elasticity based on the classical theory. If a body’s intimate structure is considered for behavior, classical elasticity theory does not offer a good description, for example, in polymers, human and animal bones, and graphite. This is even more pronounced in materials with specific porous or granular internal structures. For more information, we recommend [
1,
2].
The first researcher who approached microstructure theories in a large number of studies was Eringen, who is considered a true pioneer in this field. We suggest [
3,
4]. After the way was opened by Eringen, many other researchers tackled a large number of aspects of microstructure, such as [
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21]. In [
7], Ciarletta investigated the bending of microstretch elastic material plates with no thermal effects using the main ideas of Eringen. In [
5], Iesan and Pompei find a exact solution of the Boussinesq–Somigliana–Galerkin type for a boundary-value problem in the same context of an elastic material with microstretch structure. In [
12,
14], the authors present some aspects regarding the propagation of waves through bodies with different microstructures. One of the first and most important microstructures is the Cosserat type. The Cosserat continuum, proposed by the French brothers Cosserat, introduces the rotational kinematic degree of freedom in addition to the translation components of the Cauchy model. A large number of influential scientists have granted great importance to this type of structure. A long time ago, the study of crack propagation was raised, at first, obviously, in the case of elastic or thermoelastic media. The first studies dedicated to this matter include Eringen [
22], Freund [
23], Gurtin [
24], Gurtin and Yatomi [
25], and so on. For instance, Freund [
18] proved that energy flux depends on the elastic field near the tip. Further, Gurtin [
24] discussed some thermodynamically cohesive models of mechanical fractures.
Gurtin and Yatomi [
25] provided some expressions for the energy release rates in the presence of a propagating crack for a body that is hyperelastic and is subject to finite deformation.
The authors of [
26] address the physical meaning of path-independent integrals for elastodynamically propagating cracks, the relation of these integrals to the energy release rates for propagating cracks, and the relation between these integrals and time-dependent stress-intensity factors in general mixed-mode dynamic crack propagation. Recently, some studies have presented new results regarding the energy release rate for crack propagation in bodies with microstructures. For examples, Zhang and Wang [
27] theoretically investigated magneto–thermo–electro–elastic crack branching in these materials. A mathematical model was presented by Tian and Rajapakse [
28] to find the parameters of fracture in the case of a finite impermeable crack in a magneto–electro–elastic plane actuated by magneto–electro–mechanical loading.
A study of a semi-infinite crack that is constantly propagated in an elastic material was conduected by Gourgiotis and Piccolroaz [
29]. The body is subjected to some plane-strain shear loadings.
Our present paper tackles the problem of the propagation of a straight crack in a Cosserat elastic material.
2. Preliminaries
An open region D of the Euclidean three-dimensional space is considered, and an elastic Cosserat material fills this domain. As usual, the closure of D is denoted by , and its boundary is assumed to be a regular surface so as to allow application of the divergence theorem. We adopt the Cartesian vector and tensor notation and use two types of variables: the points of D have the spatial variables , and we use the temporal variable t, with . When the meaning is unambiguous, we do not specify the dependence of a function on its spatial and/or temporal variables. We use the known rule of summation in the case of repeated indices.
Regarding the differentiation of functions, we use the following two rules: (1) a superposed dot on a function denotes the partial differentiation of the function with respect to t; and (2) the partial differentiation with respect to spatial variables will be designated with the respective subscript preceded by a comma: .
Notations for the mechanical quantities:
- -
the components of the displacement vector field;
- -
the components of the microrotation vector;
- -
the components of the strain tensor;
- -
the components of the couple strain tensor;
- -
the permutation symbol (the Ricci’s tensor);
- -
the components of the stress tensor;
- -
the components of the couple stress tensor;
- -
the reference mass density;
- -
the components of the microinertia;
- -
the components of the body force;
- -
the components of the couple body force;
- -
are the director cosines of the unit normal, outward to the boundary ;
- -
the components of the surface traction vector;
- -
the components of the couple surface traction vector;
- -
the dynamic energy release rate;
The deformation of a Cosserat elastic body will be described using the following internal constitutive variables, as functions depending on the time and on point position :
- -
the component elements of the vector of displacement;
- -
the component elements of the tensor of Cosserat displacement.
The strain tensors have the component elements
and
, and these are introduced with help of the components of displacement by means of the next kinematic relations:
Because our subsequent considerations are made only in the context of a linear theory, it is normal to consider that the density of the internal energy
E is a quadratic form regarding all its constitutive variables. Thus, by assuming that our body in its reference state is free from stress and has zero intrinsic body couples,
E has the following expression:
Taking into account the expression of the density of internal energy using Clausius–Duhem inequality, we can obtain the constitutive relations that give the expressions for the stress tensors
and
in terms of the strain tensors:
obtained in the case of a homogenous Cosserat elastic body.
Above tensors
,
, and
are elastic constants that characterize the material properties and obey these relations of symmetry:
By using the procedure of Eringen [
4], we can deduce the equations of motion. If it is assumed that mass forces and the couples of forces are missing, then the motion equations have the following form:
in which
is the density of the body in its reference state, and
are the components of the tensor of inertia, with
.
Further, is the permutation symbol (the Ricci’s tensor).
We assume that Equations (1), (3), and (5) take place for all .
Further, we suppose that our body D in moment t has a straight crack, denoted by , and the tip of this crack at moment t is denoted by . Assume that is a function of class regarding the time variable t, and its velocity is .
Let us denote by the domain covered by for all .
In what follows, we consider certain applications of the form
, depending on
and on
. According to a suggestion given by Gurtin and Yatomi in [
20], we introduce the next definition.
Definition 1. We say that a function f is called a fracture field for a positive integer m if the following two conditions are met:
(1) f admits derivatives of any order i, , except for the points in the crack;
(2) all derivatives , , are continuous applications except for the points in the crack, and functions are continuous at points approaching the crack from both sides except in points at the tip of the crack.
Let us specify some theoretical issues that we use below.
(i) The notation means that we are in the case for any ;
(ii) Let us introduce the family
of smooth surfaces
of
D, having the property:
as
. If
, we have
(iii) If the limit from (6) takes place uniformly with regards to , we consider that belongs uniformly in time.
(iv) Consider that
is the vector of the position from the tip. If
f is a function that depends on
t in the form
, then the derivative of
f with respect to
t, denoted by
, is computed as follows:
holding
fixed. Here, we use the notation
(v) In the case of three scalar functions
u,
v, and
w with properties
, and
, we use the rule:
in which
is a normal versor oriented outward from
.
(vi) If
is a
fracture field so that
, then the divergence theorem gives:
called the “flow-divergence formula”.
Statements (v) and (vi) above are easy to prove. However, for more details, we recommend [
20].
3. Basic Results
We assume that the displacement vector of components and the couple displacement vectors of components are fracture fields in the sense of the above definition. Further, we suppose that the tensors of stress and are fracture fields.
We can now define the surface traction vector of components
and the couple surface traction vector of components
with the help of the next equations:
where
are the director cosines of the unit normal outward to the border
.
By definition, the application
defined in any
t of the interval
is called
a rate of the dynamic energy release if
has the following expression:
We will consider
S a surface that includes its inside at time
at the tip of the crack, where
and so
and
intersects the crack
only once.
It can be shown that the rate of dynamic energy release does not depend on the choice of surface S.
Now, we prove a property of the crack tip.
Proposition 1. The rate of dynamic energy release is intrinsic to the tip of the crack.
Proof. First of all, observing the expression of the internal density, we can write:
Then, considering the theorem of divergence and Equation (
8), we can write:
Further, if we consider the motion Equations (5) and (10), we deduce:
After these calculations, dynamic energy release rate receives the form:
which is even (9), and thus we conclude the proof of this proposition. □
In our next result, we obtain another expression for the rate of dynamic energy release.
Theorem 1. For any , the rate of dynamic energy release can be written in the following form:where S is a surface that includes the inside of the tip of the crack. Further, represent the elements of the speed of the tip , and appears in (7). Proof. If we consider Rule (7), we obtain
If we take into account the kinematic Equation (
1), we are led to
where the following notations are used:
Applying the second derivative, according to the rule in (7), we deduce:
We obtain a relation similar to the one in (11) using (7) again:
and by using (16), we also deduce that
Now, if we take into account Equation (
2) and use the relation (7), the following equality is obtained
so that from previous equations we deduce:
where
represent the elements of the speed of the tip, and
is the outward normal to the border
.
By taking into account Equations (11), (12), (17) and (20), the following relation is obtained:
Finally, we take into account that surface
S includes the tip of its crack in its inside at time
t. Based on Proposition 1, we know that the rate of the dynamic energy release
is intrinsic relative to the tip of the crack, as such,
can be written in the form:
In this last relation, we take into account that the area of
S is tending to zero, and so we obtain (13), which concludes the demonstration of Theorem 1. □
In Theorem 2 that follows, we obtain another expression of the rate of dynamic energy release at an arbitrary fixed time .
Theorem 2. For a fixed time , the rate of dynamic energy release can be written in the form: Proof. We start with the observation that:
Then, we have
Now, it is easy to verify that we have:
According to these last two equations, we deduce that:
Now, we consider the motion Equation (
5), the kinematic relations (1) and the tractions on the surface from (8), so that we obtain the next identity:
According to the Equation (
24), one can write:
Further, it is not difficult to verify that:
Now, by integrating over
D for both sides of Equation (
26), we deduce:
By using some similar calculus, we can obtain an analogous relation:
Now, we consider Equation (
10) in order to obtain:
Ultimately, we consider the relations (27)–(29) so that from (25), we arrive at identity (22), and the proof of the theorem ends. □
Our results are placed within a linear theory of Cosserat elastic media. After formulating the basic equations and conditions in this context, we approached the rate of dynamic energy release due to the propagation of a straight crack in this type of material. We insisted on highlighting the effect of microinertia and the couples stresses on this rate, these being contributions of the Cosserat structure. Now, we emphasize the fact that compared to classical elasticity:
a. There was only one balance equation, namely ; now there is another one, namely , in which microinertia appears.
b. There was only one stress tensor; now there is another one that appears in both balance equations.
c. There was only one strain tensor; now there is another one that appears in both constitutive equations.
d. There was only one constitutive equation; now there is another one.
It is not difficult to observe that these contributions are not so obvious as to cause spectacular and meaningful changes in the propagation of a straight crack. In our basic results, concretized in estimations (13) and (22), we deduced certain alternative formulations for the rate of dynamic energy release. These refer to some boundary concrete integrals calculated on a surface that includes the crack in its inside, and this surface shrinks to the tip of the crack.