Open access peer-reviewed chapter

# Appearance of Catastrophes and Plasticity in Porous and Cracked Media

By Boris Sibiryakov

Submitted: February 7th 2019Reviewed: May 23rd 2019Published: August 14th 2019

DOI: 10.5772/intechopen.87014

## Abstract

This chapter is devoted to study the properties of structured continuum, with specific surface and characteristic size of structure. This linear dimension means the absence of automatic transforming difference relations into differential equations. It is impossible to apply conservation laws at any point of the real structural body, because any closed points in vicinity of inner surface can represent both solid and liquid (gas) phases. We need use some representative minimal volume, which characterized the complicate body at hole. This approach leads to differential equations of motion of the infinite order. Solutions of them, along usual P and S waves, contain many waves with abnormally low velocities, which are not bounded below. It is shown that in such media, weak perturbations can increase or decrease without limit. The reason of the infinite order of differential equations is many degrees of freedom in such media. Catastrophes correspond to unstable solutions equations of motion. Plasticity begins in elastic state like continuous phenomenon, and there is a finite distance between the sliding lines on the contrary with classic plasticity, where distances between sliding lines are infinitely small.

### Keywords

• structure of pore space
• porous and cracked media
• instability
• plasticity

## 1. Introduction

The main idea of continuous mechanics is that any volume is the representative one. It means that the integral of loadings, which concentrates on the surface and bounds mentioned volume, is equal to zero in statics or to inertial forces in dynamics. The evident disagreement that the surface forces and inertial ones apply to different points (inertial forces apply to center of gravity of volume) overcomes due to an assumption about infinite small sizes of the mentioned volume. This assumption gives us a possibility to equal the volume forces (divergence of the stress tensor), which was created by the internal stresses, and the inertial forces, according to the second Newton law. Mathematical technique is based on the Gauss theorem about relation between the field flux across surface and divergence of this field in the volume, which is bounded by closed surface. However, in the structured bodies, there is a fundamentally different situation. The representative volume must contain some set of elementary structures. Otherwise, a small volume will contain only one of the phases, for example, liquid in the pores or the solid skeleton without liquid, and will not characterize the properties of the structured body. The characteristic size of the structure leads to fact that the average distance is between one of the cracks to another and one pore to another given by the specific surface of the sample. It is necessary to connect the integral geometric properties of a medium with physical processes of such bodies deforming. On the contrary with a classic continuum of Cauchy and Poisson, the new continuum for structured or blocked media must contain many degrees of freedom. It is evident because elementary blocks may translate the motion by contact interactions, by rotations, and by group of particle’s motion. It means that the energy contents not in first derivatives (strains) only. The potential energy contents in the second derivatives (curvatures) and other orders of ones. It means that the equation of motion of a blocked medium should contain many derivatives; in other words, the equation of motion may have been very high, probably, the infinite order. The static and dynamic processes in the classic continuum are divided by the Great Wall of China from each other. The equation of equilibrium never will pass in the equation of motion. However, it is evident that the dynamic processes often arise very slow and are quasi-static motions. It would be nice to destroy this mentioned wall by a newly structured continuum. It would be a good idea to destroy the abovementioned wall by means of justification of the newly structured continuum. The seismic emission, which causes due to static loading, maybe not a bad example of such phenomena, which are existed between statics and dynamics.

## 2. Equations of motion for structured media

In Figure 1, an element of the volume of structured body is shown, in which l0 is the average distance between one pore and another. Earlier presented was the result about the relation between the specific surface and the average length between cracks and pores. There is a theorem of integral geometry, which relates the specific surface σ0 and l0, namely Figure 1.Representative element of structured body for granular medium (left) and average distance l0 from one crack to another (right). On surface C, the equation of equilibrium is complied, and on surface D, it is not satisfied.
σ0l0=41fE1

where f is the porosity. Hence, if there is a specific surface of sample, there is automatically the average range of microstructure l0.

The distinction between classic and structured continuums is clear, see Figure 1. In the volume, which is inside into surface C, there is equation of equilibrium, because all forces delete to each other. In the volume, which is inside into surface D, there is equilibrium, because forces do not compensate to each other (on the one side of grain, we have forces; and on the other part of boundary surface D, we have no forces).

The idea of creation of the new model of space is as follows: consider some finite volume of the body (a sphere on a figure with radius l0). Surface forces act on a sphere of radius l0, while inertial forces applied at the center of the structure. There is no way for the volume element to tend to zero and to match the points of application of surface forces and inertial forces, as in the classical continuum. Therefore, since we must consider the representative finite volume, we have a problem of different positions of surface and inertial forces.

We need to translate the surface forces to the center of the structure by a special operator, and after this, it is possible to apply the law of conservation for some structural image continuum and to act as in a typical classical model of space. The main feature of this approach is to fill all the space, including the pores and cracks by field force. Because of it, we have a continuous image of a very complex media and a possibility to apply the physical laws into an image of the media.

The one-dimensional operator of field translation from point x into point x ± l0 is given by the symbolic formula 

ux±l0=expl0DxE2

The operator is Dx=x. The difference operator 1xis a difference between two translation operators

1=1l0ux+l02uxl02
=uxl0expl02Dxexpl02Dx=uxsinhl02Dxl02E3

This is a first difference for finite distance between two points. The second difference may be represented as quadrate of the first difference,

2=uxsinhl02Dx2l022E4

The formally expansion in Taylor’s series gives a finite increment of field. This expansion contains the infinite number of derivatives with different powers of l0. The factor l0 relates with the specific surface of the sample. The three-dimensional operator of field’s translation for some cube with length of l0 may be constructed as follows:

Pux=ux6coshl02Dx+coshl02Dy+coshl02DzE5

The analogous operator of translation for some spheres is given by expression

Pl0Dxl0Dyl0Dz=14π02π0πexpl0Dxsinθcosφ+Dysinθsinφ+DzcosθsinθdθdφE6

Because there is a Poisson formula 

02π0πfαsinθcosφ+βsinθsinφ+γcosθsinθdθdφ=2π0πfRcospsinpdp=2π11fRtdtE7

In the formula (7), parameters α,β,andγare some quantities. However, in Eqs. (6) and (7), parameters play the role of differential operators. The relation between quantities and operators is established by Maslov . Hence, P operator maybe rewritten as follows 

Pl0Dxl0Dyl0Dz=1211expl0tdt=01coshl0tdt=sinhl0l0=E+l023!+l045!∆∆+E8

In the classic continuum, we apply the impulse conservation law to any element of the medium. In this situation, we need to fill all pores over space by a force field. Instead of real stresses, which are changing very fast from one point to another, we can construct the continual image of real stresses. Namely, we use a continuous field, which is constructed by the application of the operator P to the real complicated force field. For this continuous image of real stress, Pσik,we can apply the impulse conservation law. In the classic continuous model, this operation is made by nature itself. This model of a continuum requires some mathematical operations in order to create the continuum medium. Using operator P, we can write the equation of motion of micro-inhomogeneous body, because for an average stresses in structure, the law of impulse conservation takes the usual form, namely 

xkPσik=ρu¨iE9

In a more detailed form Eq. (9) can be rewritten as follows

xkE+l023!+l045!∆∆+σik=ρu¨iE10

No wonder that Eq. (9) contains derivatives of the infinite order. This circumstance is due to many degrees of freedom for structured bodies. At l00, we have the usual equations of motion for classic continuous model of space.

## 3. Fundamental solutions

We can pass to the image space, following Hooke’s law and applying the Fourier transform along three coordinates, as 

uixyz=12π3expinxx+nyy+nzzUinxnynzdnE11

where n2=nx2+ny2+nz2;dn=dnxdnydnz. The operator P leads to

Puixyz=12π3sinhl0nl0nexpinxx+nyy+nzzUinxnynzdnE12

This allows us to calculate the Fourier transform for the fundamental solution of the system Eq. (9):

Gij=1μn2ρω2l0nsinl0nδijλ+μninjλ+2μn2ρω2l0nsinl0nl0nsinl0nE13

At very small values, l0n, the sine and argument ratio approaches unity, and the Fourier transform becomes an ordinary equation for Green’s tensor in an elastic continuum. The inverse Fourier transform is obtained by integration of Eq. (13) which includes simple poles corresponding to P and S waves and a set of simple poles where the sine in the denominator of Eq. (13) becomes zero. The residuals are in the simple poles, n2=kS2l0n,where kS2is the wave number of both P and S waves. At very small l0,the ratio l0nsinl0n1and denominators in Eq. (13) become the classical equations that define the poles corresponding to compression and shear waves velocities (Figure 2). Assuming n/ks = m and ksl0 = ε, we obtain the equation for complex roots that describe waves from a focused source in porous and cracked solids as Figure 2.Wave number ratio as a function of dimensionless ratio ε = 2πl0/λs. Curves: 1—wave number ratio ks(ω)/ks(0), i.e., S-wave velocity decreasing with frequency; 2—γ = Vs/Vp increasing with frequency; and 3—wave number ratio ks(ω)/ks(0) of P waves.
msinεm=εE14

If m = x + iy is assumed to be a complex value, for the real and imaginary parts, we have the transcendental equations

xsinεxcoshεyysinhεycosεx=ε2
ysinεxcoshεy+xsinhεycosεx=0E15

We can rewrite Eq. (15) in a different form with x* = εx and y* = εy as new variables

tanxx=tanhyy;sin2x+sinh2y=ε4x2+y2E16

Equation (15) obviously has many real roots corresponding to у = 0. Indeed, at small ε, Eq. (15) gives the solution m = 1, which corresponds to the ordinary P- or S-wave velocity. At large values of m, Eq. (15) is satisfied only if εm approaches a value divisible by n, i.e., at near-zero sine that defines the characteristic anomalous velocity. The unbounded value of the wave number means that normal P and S waves coexist with arbitrarily small P and S velocity anomalies. The existence of these anomalies in a micro-heterogeneous medium has its physical explanation: energy is stored in strain (first derivatives of displacement) as well as in the curvature of higher derivatives. Therefore, there appear as velocities related to flexural and torsion waves and to numerous waves associated with oscillation of groups of particles (blocks) (Figure 3). Figure 3.The decreasing P-wave velocity (the upper line) and S-wave velocity (the middle curve) and the growth of their ratio γ = VS/VP (the lower line) due to increasing size of microstructure. The ratio γ = VS/VP more than 0.705 corresponds to the negative Poisson ratio . The vertical scale is the wave velocities (km/s) and a horizontal scale is the ratio between the size of the microstructure and wavelength.

The growing of ratio γ = VS/VP causes a very interesting phenomenon, namely an apparent negative Poisson value, for waves with the length not very small compared to size of a grain. The growing value of γ = VS/VP means that the Poisson ratio is decreasing up to negative values  (Figure 4). Figure 4.Gregory experimental data. Poison ratio (the vertical axis) versus pressure. Black color corresponds to water saturated porous shales and gray color corresponds to dry shales with the same porosity. In this case, negative Poisson ratios are possible.

At the same time, Eqs. (14) and (15) likewise have complex roots. The first Eq. (15) shows that complex roots arise only at some values of e, which are not so small, as they satisfy the inequality εx > π/2. Table 1 lists complex roots corresponding to some relatively small ε. Note that the parameter εcan be expressed via the linear size-to-wavelength ratio (l0/λs).

εxy
0.21472.02880.0548
0.25072.06450.5838
0.27712.10640.8880
0.32532.15601.1838
0.39182.21571.5122

### Table 1.

The value epsilon means dimensionless product of structure size into wavenumber of usual S waves in continuous medium. Value x means the real value of product of structural wavenumber into structure size. Value y is the imaginary part of it.

Complex roots can mean either damping or unlimited growth of wave amplitude, of course, in the presence of an energy-unbounded source. The minimum damping (growth) corresponds to (2.0288)−1 or about a half of the normal velocity. The same process can be expected to cause both excitation and damping in porous and cracked media, depending on the phase of stationary oscillations.

## 4. One-dimensional case: plane wave and instabilities

In one-dimensional case, the Eq. (10) takes more simple expression

uE+l023!+l045!∆∆++kS2u=0E17

This equation by substitution u=expikxgives us the dispersion equation for an unknown wave number k, or for unknown wave velocity, which depends on the size of structure l0 or specific surface of sample σ0:

sinkl0kl0=kS2k2E18

It is evident that by l0 → 0, the wave number k → kS, i.e., the wave velocity is equal to VP or VS, elastic wave velocity. However, if l0 is not a very small value, the wave velocity decreases up to zero by kl0 → , if m is the integer number. Hence, this model describes along with usual seismic waves many waves of very small velocities, which are not bound below.

This effect is more for P waves than for S ones. Eq. (14) shows that if the Poisson ratio is measured on samples by velocities VP and VS, their ratio VS/VP grows by growing l0, and this effect can produce abnormally small Poisson’s ratio, up to negative volume of it.

It is evident that at kl=, m is the integer number and the value k. It means that there are waves with arbitrary small velocities not bounded below. Beside it, Eq. (15) has complex roots too, because sinkl0may be negative, while the second term in Eq. (18) contains kS/k2. Eq. (18) means that the complex roots do not by small values of x, because the right-hand expression is a negative value. In order to be complex roots, an evident condition is necessary, i.e., tanx>π/2. The physical sense of it means that the complex roots are possible, if the wavelength is four times (or more than four times) more than the size of the structure. These complex roots mean that amplitude of oscillations may be increasing or decreasing up to infinity or, may be, to zero. These roots are responsible for catastrophe’s behavior of structured bodies.

Hence, if there is a source of sufficient energy, even some small oscillations can produce catastrophes. It is interesting that nonlinear deforming of samples decreases this effect, because a wave velocity for rocks is decreasing, by growing amplitude of wave. It means that the wave number is growing by the same frequency in the pure elastic process. In Figure 6, the real roots of dispersion, Eq. (18) are shown. The vertical axis shows a dimensionless frequency, namely ε, while horizontal axis shows us the real and imaginary parts of wave numbers. In Figures 5, 6, 7 , complex roots as a function of dimensionless frequency ε are shown. Every point is a position of some root, namely a real part, an imaginary one, and a dimensionless frequency. The more is the spreading of ε values, the greater is the number of complex roots. Figure 5.The position of complex roots depends on the value ε . The more the value ε , the more numbers of roots. The first value of ε corresponds to values of first row from Table 1. Figure 6.The position of complex roots depends on the value ε . The more the value ε , the more numbers of roots. The second value of ε corresponds to values of second row from Table 1. Figure 7.The position of complex roots depends on the value ε . The more the value ε , the more numbers of roots. The third value of ε corresponds to values of third row from Table 1.

## 5. Pointing vector and equation of equilibrium for blocked media

The equation of equilibrium for micro-structured media can be written from Eq. (9) as

xkPσik=Pσikxk=0;σikxk=P10E19

The inverse operator P10contains zero, but not zero only. It contains some periodic functions and the average value equal to zero. For example, such construction satisfies to Eq. (19)

PIm{expinπl0k1x+k2y+k3z=0E20

If n is integer number, k1=sinθcosφ;k2=sinθcosφ;k3=cosθ.Physical sense of it means that volume forces are equal to zero in average sense, not at any point. Using mentioned inverse operator, we can write the equilibrium equation for blocked media in the form

σikxk=ϕσ02σik0uk0n=1Imexpinπl0k1x+k2y+k3zE21

In Eq. (21) σ02is a quadrat of specific surface; σik0uk0=Ai0is the pointing vector of usual continuous body, andϕis a dimensionless constant, which must be obtained. These values we can put as constants in small structure volume. The integration with respect to spherical angles gives us a result that the imaginary part of exponent is zero in average sense, namely

14π02π0πexpπnl0ixsinθcosφ+ysinθsinφ+zcosθsinθdθdφ
0πexpirl0cospsinpdp=1211expirl0tdt=l0nπrsinrnπl0+i0E22

Partial solution of Eq. (22) is a convolution of Green tensor with right hand of Eq. (21), that is,

ui1x=ϕ1μσ02σmk0uk0xImΓmixyexpikmxmymdVyE23

Taking into account that the sizes of area much more, than sizes of structure, the area of integration is the infinite large one. In this case, integral Eq. (13) practically is the Fourier transform of fundamental solution of usual elastic equilibrium equations

ui1nx=ϕ1μσ2σmk0uk0xl02δmi1γ2kmkiexpiklxl)l0E24

In Eq. (24) the imaginary part of the exponent is used. Hence, the additional value in average sense is equal to zero. Using relation Eq. (1) σ0l0=41f, we get a partial solution, which depends on porosity only

ui1nx=ϕ1μσmk0uk0x41f2δmi1γ2kmkiexpiklxll0E25

If these indexes coincide, i=m,we get

<ui1nx>=ϕμUi041f211γ23expiklxl)l0E26

Take into account that the average value of a quadrat of cosine is <kmki>=δkm3. There is a summation with respect to n,and Ui0is a Pointing vector for usual continuous model of the medium. This value is a small one of the second order compared to usual displacement, because a Pointing vector, divided on the shear module is order to strain, multiplied to size of structure l0.

Strains. By differentiating of an integral Eq. (23) take into account that the main part of the field contains in fast changing exponent, not in Green tensor itself, i.е.,

ui,j1xϕμσ2σmk0uk0xikjImΓmixyexpiklxlyldVyE27
eij1n=i2ui,j1+uj,i1=ϕ2μσmk0uk0xl04π41f2{δmi1γ2kmkikj
+δmj1γ2kmkjki}1nexpik1x+k2y+k3zl0=ϕ8γ23πμ1l0U0ikj+U0jki1f21nexpiklxll0E28

According to Eq. (9) the additional dilatation is

θn=ϕ16γ23πμ1l0U0n1f21nexpiklxll0E29

Let us integrate the normal component of the Pointing vector on the small sphere with radius r. This integral must be equal to density of potential energy E(divergence of Pointing vector) namely,

<U0n>=Eσ0=41fEl0E30

The average value of fast-changing exponent in Eqs. (28) and (29) on spherical angles is

<exp1niklxll0>=41fπn20sinxxdx=41fπn2SiE31

The additional dilatation due to randomly oriented volume forces (an average value of these forces is zero) may be written as

θ=ϕ16γ23πμ41f3Eπn=1Sin2E32

In Eq. (32) the symbol Simeans an integral sine of argument . The left hand in Eq. (32) is an additional expansion or compression, so called as dilatancy. It depends on the potential energy of the continuous body E, which may contain shear energy only, but it produces additional expansion or compression. It is a quadrat effect too.

More strong effect is related with product of high-changing volume force (equal to zero in average) into displacement. This product in not equal to zero in average, because it contains a quadrat of high-changing sine, which is equal to number one third in three dimension space.

En=ϕ1μσmk0uk0x41f2δmi1γ2kmki13ϕσ02U0iE33

If indexes coincide, m=i, we get the additional potential energy, due to fluctuations

En=ϕ23λ+2μU012+U022+U032σ0241f2E34

The summation with respect to index nfrom unit up to infinity gives

E=8ϕ29λ+2μσ02U012+U022+U0321f2E35

In spite of a fact that the Pointing vector is the small value of more high order, than stresses, the high value σ02(quadrat of specific surface) in Eq. (35) can produce not small common effect. The indefinite factor ϕdepends on the real structure of pore space and macro-stress-strain state. However, in some simple situations, it can calculate elementary. For example, at rigid pressing of globe by spherical force (radial displacements are constants), the stress-strain state is a hydrostatic state in average, but not such state at any point. The compressional energy is proportional to compress module of skeleton and its volume plus the incompressibility of fluid and its volume, namely [8, 9, 10, 11]

E=λ+2μ3θ1221f+ρc2θ022fE36

Indexes unit and zero in Eq. (36) mean solid and liquid parameters. The dilatation of two-phase body gives by the formula

θ=1fθ1+fθ0;θ0=θ1E37

If we have uniform random distribution of phases, the average energy is

E=E11f+E0fE38

In Eq. (38) fis the porosity and E1andE0are the energies of solid and liquid. The dispersion of random value relates with random volume forces, i.e.,

E1E1(1f+E0f)21f+fE0E1(1f+E0f)2=E1E02f1fE39

Equation (39) gives the additional energy for very simple macro-hydrostatic state in average. This is the additional of interphase acting. It is equal to additional energy, which is given by Eq. (15). It is reasonable that at unit or zero porosity, an additional energy is equal to zero. The second result is, if the phase energy is equal, the mentioned additional one is equal to zero too. Hence, the indefinite factor ϕ2given by the simple equation is

1μE1E0f1f=8ϕ29λ+2μσ02U012+U022+U0321f2E40

## 6. The arriving of plasticity

In spite of that, the additional average strains is small, does not means, that these strains are small in the any point of the volume. Equations (28) and (29) show that on the planes k1x±k2y±k3z=2l0q(q is an integer number), the exponent is not a highly changed value, because it is equal to1or −1.

In plane situation, the role of these planes plays orthogonal lines k1x±k2y=2l0q.

The series Eqs. (28) and (29) with respect to n in vicinity of mentioned planes are divergent (harmonic) series. It means that the field is decreasing inside of quadrats, making planes, and concentrating in vicinity of planes. Mentioned planes are analogs of slipping lines (lines of Luders)  in classic plasticity of the compressible medium. In practical, the number n in Eqs. (28) and (29) is bounded by the elastic limit of the second strain invariant. The field of strains is growing into planes (lines) and decreasing inside of them. This process is called as localization of strains. This localization begins in elasticity, with contrary of classic plasticity and elasticity. The other specific feature of this process is the finite distance between planes (lines). This distance is equal to l0 (the inverse value of specific surface of sample), while in classic plasticity, this distance is infinitely small. The geological sense of it is interesting. In order to transform the matter from elasticity to plasticity, there is no necessary to have the plastic state at any point of the medium. Plasticity may concentrated near planes, and the other volume can be in elastic state. Rock may flow comparatively light, if they have pores and cracks. On the Figure 8 shown successive process of localization of strains due to decreasing field inside of quadrats, making by orthogonal lines and increasing them near lines itself. Figures 8.Successive process of strains localization due to decreasing strains inside of quadrats, making orthogonal lines and increasing them near lines itself.

## 7. Conclusions

1. The model of the structured continuum with specific surface of the blocked medium or average size of structure, gives us the differential equations of motion of the infinite order. This model includes collective properties of pore space like the porosity and specific surface and predicts besides usual elastic waves many unusual waves with very small velocities.

2. This model predicts the decreasing of the Poisson ratio (up to negative values) due to finite size of microstructure. The reason for this is the decreasing of wave velocity with finite specific surface of the rock.

3. The localization of stresses and strains in structured media begins in elastic state of deforming.

4. The small areas of a stress-strain concentration looks like usual orthogonal sliding lines in classic plasticity. However, they have a finite effective thickness, which depends on the average size of the structure and the elastic strain limit. Besides, there is a finite distance between analogs of sliding lines, which is equal to the average distance from one pore to another one, or between cracks.

## How to cite and reference

### Cite this chapter Copy to clipboard

Boris Sibiryakov (August 14th 2019). Appearance of Catastrophes and Plasticity in Porous and Cracked Media, Seismic Waves - Probing Earth System, Masaki Kanao and Genti Toyokuni, IntechOpen, DOI: 10.5772/intechopen.87014. Available from:

### Related Content

Next chapter

#### Efficient Simulation of Fluids

By Pierre Thuillier Le Gac, Emmanuelle Darles, Pierre-Yves Louis and Lilian Aveneau

First chapter

#### Seismic Wave Interactions Between the Atmosphere - Ocean - Cryosphere System and the Geosphere in Polar Regions

By Masaki Kanao, Alessia Maggi, Yoshiaki Ishihara, Masa-yuki Yamamoto, Kazunari Nawa, Akira Yamada, Terry Wilson, Tetsuto Himeno, Genchi Toyokuni, Seiji Tsuboi, Yoko Tono and Kent Anderson

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books