Advances in Mathematical Physics

Volume 2016 (2016), Article ID 4792148, 9 pages

http://dx.doi.org/10.1155/2016/4792148

## Elastic Equilibrium of Porous Cosserat Media with Double Porosity

Ilia Vekua Institute of Applied Mathematics of Ivane Javakhishvili Tbilisi State University, 2 University Street, 0186 Tbilisi, Georgia

Received 11 May 2016; Accepted 30 June 2016

Academic Editor: John D. Clayton

Copyright © 2016 Roman Janjgava. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The static equilibrium of porous elastic materials with double porosity is considered in the case of an elastic Cosserat medium. The corresponding three-dimensional system of differential equations is derived. Detailed consideration is given to the case of plane deformation. A two-dimensional system of equations of plane deformation is written in the complex form and its general solution is represented by means of three analytic functions of a complex variable and two solutions of Helmholtz equations. The constructed general solution enables one to solve analytically a sufficiently wide class of plane boundary value problems of the elastic equilibrium of porous Cosserat media with double porosity. A concrete boundary value problem for a concentric ring is solved.

#### 1. Introduction

A model of elastic equilibrium of porous media with double porosity was constructed in the works [1–3]. The theory justified in these papers combines the previously proposed model of Barenblatt for media with double porosity [4] and that of Biot for media with ordinary porosity [5]. For a detailed account of the development of the theory of porous media and relevant references, see [6]. Various issues related to the elastic equilibrium of bodies with double porosities are treated in [7–15].

It should be noted that all the papers mentioned above dealt with a classical (symmetric) medium. But we do not know of any works where problems of double porous elasticity would have been considered for a nonsymmetric elastic Cosserat medium [16–34]. In our opinion, the investigation of such problems is interesting from both theoretical and practical standpoints. For this reason, we considered the elastic equilibrium of porous bodies with double porosity in the case of the nonsymmetric Cosserat theory.

In Section 2 we give the basic three-dimensional equations of the static equilibrium of porous elastic materials with double porosity in the case of an elastic Cosserat medium. In Section 3 we consider the case of a plane deformed state and write the corresponding two-dimensional system of equilibrium equations in the complex form. In Section 4 we construct the general solution of the abovementioned system of equations by means of analytic functions of a complex variable and solutions of Helmholtz equations. The obtained analogues of the Kolosov-Muskhelishvili formulas [35] make it possible to solve analytically plane boundary value problems of the elastic equilibrium of porous Cosserat media with double porosity. Finally, in Section 5 we solve a boundary value problem for a concentric circular ring.

#### 2. Basic Three-Dimensional Relations

Let an elastic body with double porosity occupy the domain . Denote by a point of the domain in the Cartesian coordinate system. Let the domain be filled with an elastic Cosserat medium having double porosity. The considered solid body is characterized by the displacement vector , rotation vector , and also the fluid pressures and occurring, respectively, in the pores and fissures of the porous medium.

Then a homogeneous system of static equilibrium equations is written in the form [25]where are stress tensor components, are moment stress tensor components, is the Levi-Civita symbol, and ; the summation over the recurring index is assumed to be done from 1 to 3.

Formulas that interrelate the stress and moment stress components, the displacement and rotation vector components, and the pressures , have the form [13, 25]where , are the Lamé parameters, , , , are the constants characterizing the microstructure of the considered elastic medium, and are the effective stress parameters, and is the Kronecker delta.

As is well known, for the internal energy to be positive it is necessary that the following conditions be fulfilled [28]:

In the stationary case, the values and satisfy the following system of [13]where , , , and , is fluid viscosity, and are the macroscopic intrinsic permeabilities associated with matrix and fissure porosity, and are the cross-coupling permeabilities for fluid flow at the interface between the matrix and fissure phases, is the internal transport coefficient and corresponds to fluid transfer rate with respect to the intensity of flow between the pore and fissures, and is the three-dimensional Laplace operator.

The three-dimensional system of (1), (2), and (4) describes the static equilibrium of a porous elastic Cosserat medium with double porosity. Substituting relations (2) into (1), we obtain equilibrium equations with respect to the components of the displacement and rotation vectors:

If we add boundary conditions on the boundary of the domain , to the system of equilibrium equations, then we can consider various classical boundary value problems.

The following lemma is easy to prove.

Lemma 1. *If , , then the system of (4) is equivalent to two independent equations: to the Laplace equation with respect to the combination and to the Helmholtz equation with respect to the difference where *

*Proof. *Adding the first equation of system (4) to the second equation of this system, we immediately obtain (6).

Let us write system (4) in the matrix form: By assumption, the determinant of the matrix is positive. The left multiplication of all members of the latter equation by the matrix gives the systemIf from the first equation of the latter system we subtract the second equation of this system, then we obtain (7).

Since the transformation determinant is defined as the validity of the lemma is proved.

Corollary 2. *If on the boundary of the domain , , then throughout the body .*

This corollary follows from the fact that the homogeneous Helmholtz equation (7) with zero boundary conditions has only the trivial solution.

#### 3. The Plane Deformation Case

From the basic three-dimensional equations we obtain the basic equations for the case of plane deformation. Let be a sufficiently long cylindrical body with generatrix parallel to the -axis. Denote by the cross section of this cylindrical body; thus . In the case of plane deformation , , and , while the functions , , , , and do not depend on the coordinate (see, e.g., [29]).

As follows from formulas (2), in the case of plane deformation Therefore the system of equilibrium equations (1) takes the form Relations (2) are rewritten as where

Equations (6) and (7) take the form where is the Laplace operator in two dimensions.

If relations (13) are substituted into system (12), then we obtain the following system of equilibrium equations with respect to the functions , , and :

On the plane , we introduce the complex variable and the operators , , , and

To write system (12) in the complex form, the second equation of this system is multiplied by and summed up with the first equation [35–37] (see also [38, 39]):where, by formulas (13),We write (14) and (15) as If relations (18) are substituted into system (17), then system (16) is written in the complex form:

#### 4. The General Solution of System (19)-(20)

In this section, we construct the analogues of the Kolosov-Muskhelishvili formulas [35] (see also [36–39]) for system (19)-(20).

Equations (19) imply that where , is an arbitrary analytic function of a complex variable in the domain , and is an arbitrary solution of the Helmholtz equation From system (21) we easily obtain the expressions for the pressures and :

Theorem 3. *The general solution of the system of (20) is represented as follows:where , , and are arbitrary analytic functions of a complex variable in the domain , and is an arbitrary solution of the Helmholtz equation where *

*Proof. *We take the operator out of the brackets in the left-hand part of the first equation of system (20): Since (28) is a system of Cauchy-Riemann equations, we have where is an arbitrary analytic function of .

A conjugate equation to (29) has the form Summing up (29) and (30) and taking into account that we obtain If from (29) we subtract (30) and write the expression , then we haveThe second equation of system (20) is written as Substituting formula (33) into formula (34) we obtain the equation The general solution of (35) is written in the form where is a general solution of the Helmholtz equation The multiplier has been introduced for convenience in writing our subsequent formulas.

Substituting formulas (32) and (36) into (29) and taking into account that is a solution of (37), we obtain From formulas (23) we find the following expression for the combination :Substituting the latter formula into (38), integrating over , and using we obtain formula (24) which we are proving: Thus, if the solution of system (20) is sufficiently smooth, then it is represented in the form of (24) and (25). Conversely, if expressions (24) and (25) are substituted into (20), then this system will be satisfied.

Substituting expressions (24) and (25) into formulas (18), for combinations of stress tensor components we obtain the following formulas:Thus, the general solution of a two-dimensional system of differential equations that describes the static equilibrium of a porous elastic medium with double porosity is represented by means of three analytic functions of a complex variable and two solutions of the Helmholtz equation. By an appropriate choice of these functions we can satisfy five independent classical boundary conditions.

Let mutually perpendicular unit vectors and be such that where is the unit vector directed along the -axis. The vector forms the angle with the positive direction of the -axis. Then the displacement components , , as well as the stress and moment stress components acting on an area of arbitrary orientation, are expressed by the formulas

#### 5. A Problem for a Concentric Circular Ring

In this section, we solve a concrete boundary value problem for a concentric circular ring. On the boundary of the considered domain which is free from stresses and moment stresses, the values of pressures and are given.

Let a porous elastic body with double porosity occupy the domain which is bounded by the concentric circumferences and with radii and , respectively () (Figure 1).