УДК 539.422
РАЗРУШЕНИЕ ГОРНОЙ ПОРОДЫ. НЕЛОКАЛЬНАЯ ТЕОРИЯ УПРУГОСТИ
Виталий Николаевич Доровский
Новосибирский технологический центр компании «Бейкер Хьюз», 630090, Россия, г. Новосибирск, ул. Академика Кутателадзе, 4а, доктор физико-математических наук, технический советник, тел. (383)332-94-43, e-mail: [email protected]
Егор Владимирович Вторушин
Новосибирский технологический центр компании «Бейкер Хьюз», 630090, Россия, г. Новосибирск, ул. Академика Кутателадзе, 4а, кандидат физико-математических наук, научный сотрудник, тел. (383)332-94-43, e-mail: [email protected]
Евгений Игоревич Роменский
Институт математики им. С. Л. Соболева, 630090, Россия, г. Новосибирск, пр. Академика Коптюга, 4, доктор физико-математических наук, главный научный сотрудник, тел. (383)329-76-57, e-mail: [email protected]
Рассматривается разрушение слоя горной породы в рамках нелокальной теории упругости при пластическом течении. «Нелокальность» обусловлена зависимостью энергии деформируемой среды от производных метрического тензора деформаций. Показано: при продольном нагружении выделенного слоя пространственное распределение «поперечного» напряжения испытывает периодическое изменение. Максимальные значения напряжения достигаются в нестационарном режиме. Последующие разрушения горной породы ассоциируются с зонами максимального поперечного нагружения.
Ключевые слова: пластическое течение, разрушение горной породы, нелокальная теория упругости.
ROCK FAILURE AND NONLOCAL ELASTICITY THEORY
Vitaly N. Dorovsky
Novosibirsk Technology Center «Baker Hughes», 630090, Russia, Novosibirsk, 4a Kutateladze St., D. Sc., Technical Advisor, tel. (383)332-94-43, e-mail: [email protected]
Egor V. Vtorushin
Novosibirsk Technology Center «Baker Hughes», 630090, Russia, Novosibirsk, 4a Kutateladze St., Ph. D., Scientist, tel. (383)332-94-43, e-mail: [email protected]
Evgeniy I. Romenski
Sobolev Institute of Mathematics, 630090, Russia, Novosibirsk, 4 Koptyug Prospect, D. Sc., Chief Researcher, tel. (383)329-76-57, e-mail: [email protected]
The paper considers failure of a stratum within the framework of a nonlocal elasticity theory for a plastic flow. Here, the nonlocality is determined by dependence between deformation energy and metric deformation tensor derivatives. Our study has demonstrated that in case of a stratum under lateral stress, the spatial distribution of shear stress changes periodically. The maximum stress values are obtained in nonsteady regime. The following rock failures are associated with the zones of maximum shear loading.
Key words: plastic flow, rock fracture, nonlocal elastic theory.
In the classical elasticity theory, deformations are characterized by a metric tensor of elastic deformations glk. Plastic deformations also require considering a curvature tensor. IN their turn, nonelastic deformations take the deformation field out from Euclidean space. In [1] curvature is considered as a parameter of state having effect on deformation energy that makes the theory nonlocal. Unfortunately, the authors have provided no solid arguments for the introduction of additional degrees of freedom into the first law of thermodynamics, but still there is a whole class of problems in the failure theory, where the considered approach can be regarded as a reasonable one.
Let's consider a thin rod bended in a way its curvature radius at the local bend point is equal to 1/R. If one draws a line along the rod, the rod's upper part will be extended by 5/1, while the lower part will be compressed by S/2 relative to the rod. The first law of thermodynamics that unambiguously determines the process's dynamics, for a one -dimensional rod can be written as dE0 = TdS + f1d/1 + f2d/2, where E0 denotes the inner energy, T, S - temperature and entropy, f d/x, f2d/2 - the forces affecting the rod while extension and compression. The parameters S/x, S/2 can be altered by the length of 'average' line and the curvature according to 28/ = 8/ +S/2, 2tc-(-SR/R2) = S/-S/2, so the first law will be expressed as:
dE0 = TdS + (f + f2)d/ + f - f2)dR / R2.
Hence, the assumption that the rod is bended requires a degree of freedom related to the curvature to be introduced. At the same time, there are indications that the classical elasticity theory is not complete in case rock failures in the vicinity of mines are interpreted [2]. Solving a problem of stress distribution around a radial plane under a given radial centripetal stress gives us the monotonic reduction of the stress tensor's angular component. The periodic zone of rock failure observed in mines may be the evidence of periodic stress distribution before the moment when a failure occurs. The classic solution does not describe such failures.
In [1] the authors suggest the equations describing nonelastic deformations with account for inner- energy dependence on curvature. In the paper, the curvature has been considered as a destruction parameter in a steady-state problem of medium deformation around a plane. It has been shown that the stress state is periodical in the radial direction, in other words, the spatial pattern observed around cylindrical mines has been described [2]. However, the authors have not studied the properties of the nonsteady solutions. In the current paper demonstrates that it is the nonsteady process that determines failure zone evolution in time.
Equations of motion for small infinitesimal strains with a destruction parameter are expressed as [1]:
№ = dkaik, zik = (dkv{ /3= -d^sik,
aik = + 2№k + p R ■ Sik, (1) = C■Kv -4Ay), % = S{K-awSik/3)+Y-AySik/3)], y = P s + a R.
Here X, q, £ denote the elastic and plastic moduli, R = glkRik - the curvature determined by the Ricci tensor, and a, P - the parameters of the state equation
EqI p = const + r0(s -so) + Xsvv2 / 2po + iSkSk1R) + pR' sw1 ft) + aR212P0.
Alternative equations of destruction with plasticity are written as:
1 R 4
f -2 >
+ 2 3
a- P
A2R + 8 B—1--AR -
3 X + 2|
X + 2|
v
X + 2|u/3 f \ X + 2^/3^, 4 p
-u----R+--2— +--— Q----—AQ,
X + 2^i { + X + 2^i J 3X + 2^i (z)
ft) + ^ ^
i?
Changing of the density Q = -p under plastic conditions causes the failure R that satisfies the parabolic equation. As for the failure rate, it becomes a source in a wave equation determining the density change. The equations model density wave propagation and failure parameter diffusion as well as their mutual effect. The correctness of the parabolic equation is determined by the sign of the coefficient a - P21 (X + 2| > 0.
The character of the nonsteady stress-strain state under external strength loading has been analyzed using a 1-D problem. A plate with the thickness
z e
Zj, Z2 was considered. Its left boundary was under the compressive stress azz, while its right boundary remained steady (vz = 0). It was assumed that at the both boundaries y = dy / dz = 0. The system (1) had additional boundary conditions azz =-P0, y = dy I dz = 0 at z = Z and vz = 0, y=dy I dz = 0 at z = Z2. The considered medium had the following material constants: p = 2.2 kg/m3, X=12.63 GPa, I =4.95 GPa, a = 3.1110-4 kg-m3/s, P =-5.7-10-4 kg-m3/s, £=10-4 m-s/kg, ^ =10-4 m^s/kg. The thickness of the plate was 5 cm. The behavior of R can be seen in Fig. 1. In the process of reaching a steady state, the heterogeneities of maximumamplitude destruction have been concentrated next to the boundary. In the figure its spatial distribution is given for three moments of time.
The maximum rate of amplitude growth has been achieved at the first maximum value. Given the considered boundary conditions, an external force application point formally corresponds to the first destruction boundary. However, it is the value of the second of destruction extremum that was of our immediate interest (Fig. 2). This domain corresponded to the maximum shear stress value associated with the deformation part of strain a* = Xsw + 2psn = an +|P| R (Fig. 2).
When an external force affects a rock at z=0, the stress a* determined by increasing destruction rate R, starts growing until the failure criterion is met (the Ish-linsky criterion -max|2 a* |) in the zone where R reaches its maximum.
■40
Fig. 1. Destruction process at the time moments of 0.6, 3 and 300 ^s
Fig. 2. The stress a* the time moments of 0.6, 3 and 300 ^s
cc
<N
-0,5 -
Fig. 3. The destruction R and the stress a* the time moments of 0.6, 3 and 300 ^s
In the presented theory, rock failure is preceded by formation of periodic destruction zones in a medium. The maximum destruction value corresponds to the maximum destructing stress related to the deforming part of elastic stresses. The theory enables one to determine the characteristic scale of failure as a domain between two extreme values of destruction. The said failure scale can be achieved at the non-steady state of external strength loading put on a system.
REFERENCES
1. Guzev M.A., Paroshin A.A. Non-Euclidean Model of Zonal Rock Failure around Mines // Applied Mechanics and Technical Physics. - 2001. - Vol. 1. - P. 147-156.
2. Shemyakin E.I., Fisenko G.L., Kurlenya M.V. Zonal Rock Failure around Mines. Part 1. Field Study Data // Physical and Technical Problems of Mining. - 1986. - Vol. 3. - P. 3-15.
© B. H. ffopoecKUU, E. B. BmopywuH, E. H. POMCHCKUU, 2017