# Abstract

The paper describes the mathematical model for partially saturated porous media based on the Biot’s model with five basic functions to characterize wave process. The mathematical model of the boundary value problem for the three-dimension dynamic theory of poroelasticity is given in Laplace transform. On the basis of the operational calculus theorem about the original integration, the step method of numerical inversion of the Laplace transform is presented. The direct method of boundary integral equations is selected to solve value problems of the three-dimensional dynamic poroelasticity theory, the corresponding boundary integral equation is given. The corresponding matrices of fundamental and singular solutions of the three-dimensional dynamic poroelastic theory are given. A brief description of the boundary-element discretization is presented. Methodological assurance is based on the regularized boundary integral equation usage. The regularized boundary integral equations are written considering the problem of symmetry. The boundary surface of the investigated solid is divided by generalized eight-node quadrangular elements. The consistent elementwise approximation is used. Collocation solution points of the boundary integral equation coincide with the interpolation nodes of unknown boundary functions. To increase the integration accuracy for an element not containing the collocation point, the hierarchical integration algorithm and Gauss integration formulas are applied. Arising discrete analogs are solved by the Gauss-based step process to obtain the values of the boundary functions. Step process is determined by a step algorithm of the numerical Laplace transform inversion. The problem of a unit surface force shock on the free end of the prismatic partially saturated poroelastic solid is considered. Sandstone is selected as the porous material. An analytical solution of the corresponding one-dimensional problem is used for the boundary element model verification. The dependence of computational grid on convergence of the problem solution is investigated; and the impact of step scheme parameter on the solution is examined.

# Full Text

Введение Пористые материалы моделируются упругим скелетом, или матрицей, с неким распределением закрытых и открытых пор, содержащих жидкий или газообразный наполнитель. Примерами пористых материалов могут служить водонасыщенные грунты, горные породы, биологические ткани, пенометаллы и т.д. Механика пористых сред имеет большое значение во многих областях науки и техники, таких как геотехника, геомеханика, инженерная геология, биомеханика, машиностроение и материаловедение. В сейсмостойком строительстве важную роль играет прогнозирование потери прочности грунта в результате скачка порового давления жидкости, вызванного сейсмическим воздействием. Началом исследований волновых процессов в насыщенных пористых средах послужила работа Я.И. Френкеля [1]. Л.Я. Косачевский [2] показал, что теория М. Био [3,4], как и подход Я.И. Френкеля, опирается на те же соотношения между напряжениями и деформациями, но отличается большей общностью. Естественным продолжением этих работ стали попытки различных авторов обобщить теорию на случай пористых сред, насыщенных двумя несмешивающимися текучими наполнителями. Принципиальным отличием в данном случае является необходимость учитывать поверхностные натяжения и капиллярные эффекты, возникающие в результате взаимодействия наполнителей. Brutsaert [5] первым расширил теорию Био на случай трехфазной среды, однако исключив из рассмотрения капиллярные эффекты. Общая система уравнений многофазной пористой среды была предложена Николаевским [6]. Работы Garg and Nayfeh [7], Santos и др. [8, 9], Tuncay и Corapcioglu [10, 11], Wei и Muraleetharan [12], Lo и др. [13], Lu и др. [14] по теории пористых сред с двумя жидкостями принадлежат концептуально к общему подходу, отличаясь в деталях реализации. Метод в целом состоял в разработке соотношений между напряжениями и деформациями для составляющих фаз с последующим использованием их в уравнении баланса в рамках стандартной механики сплошных сред. Berryman [15] вывел уравнения для частично насыщенных сред с использованием предположения о незначительном изменении капиллярного давления при прохождении через них акустического сигнала. Smeulders и др. [16] изучали распространение волн сжатия в пористой среде, заполненной жидкостью с небольшой объемной долей газа. В работах Thomas [17, 18] предложена конечно-элементная модель для анализа волновых полей напряжений, перемещений и порового давления в ненасыщенном морском дне. Мейснер и Беккер [19] установили влияние капиллярных эффектов на динамическое поведение частично насыщенного песка. Обзор доступных аналитических решений задач о распространении волн в насыщенных пористых средах можно найти в статье Schanz [20]. Однако для ненасыщенных моделей известно немного существующих аналитических решений. Можно отметить решение, полученное в работе Li и Schanz [21], для одномерного частично насыщенного пороупругого стержня. Поскольку получение аналитического решения часто возможно лишь в некоторых специальных случаях и при определенных видах граничных условий, решение задач динамики частично насыщенных пористых сред требует привлечения развитых численных методов, таких как МКЭ и МГЭ. В случаях бесконечных и полубесконечных тел и сред метод МГЭ является наиболее предпочтительным универсальным подходом к исследованию волн. Для рассмотрения задачи о распространении волны необходимы динамические формулировки исходной системы дифференциальных уравнений и соответствующие ей ГИУ, фундаментальные и сингулярные решения. В работах Gatmiri и Jabbari [22, 23] получены фундаментальные решения для статических и квазистатических задач теории ненасыщенных грунтов в двух- и трехмерной постановке. Соответствующие термо-упругие фундаментальные решения двух и трехмерных статических и квазистатических задач для частично насыщенных сред представлены в работах Jabbari и Gatmiri [24], Gatmiri и др. [25], Maghoul и др. [26]. Динамические формулировки ГИУ и фундаментальные решения для двухмерных задач представлены Maghoul и др. в [27], для трехмерных Ashayeri и др. в [28], Li в [29]. Авторами рассматривается трехкомпонентная пороупругая среда Био с пятью базовыми функциями описания волнового процесса - перемещениями скелета и поровым давлением жидкости и газа. ГИУ получены на основе метода взвешенных невязок, фундаментальные решения представлены в терминах преобразования Лапласа по переменной времени, при построении шаговой гранично-элементной схемы используется метод квадратур свёрток. В настоящей работе представлена гранично-элементная методика на основе шагового метода численного обращения интегрального преобразования Лапласа для решения начально-краевых трехмерных задач динамики частично насыщенных пороупругих тел. Шаговый метод численного обращения преобразования Лапласа опирается на квадратурные формулы для интеграла, получаемого из теоремы операционного исчисления об интегрировании оригинала. Методика исследований основана на граничных интегральных уравнениях прямого подхода трехмерной пороупругости в преобразованиях по Лапласу, на компьютерном моделировании искомых решений методом граничных элементов в сочетании с методом коллокации, локальной поэлементной аппроксимацией на основе согласованной модели интерполирования Гольдштейна. Математическая модель Рассматривается частично насыщенный пористый материал, который представляет собой упругий скелет с порами, заполненными жидкостью и газом, при этом исключается испарение жидкости и растворение газа [29]. Кроме того, предполагается, что все три фазы являются сжимаемыми, имеют постоянную температуру, а любые изменения температуры не учитываются. Для частично насыщенного пористого материала пористость обозначается как , (1) где - объем взаимосвязанных пор в образце объема . Закрытые поры рассматриваются как часть упругого скелета. Степень насыщения материала жидкостью и газом определяется как отношение занимаемого наполнителем объема и соответственно к общему объему пор, т.е. . (2) Существенным отличием частично насыщенных пористых сред от частично ненасыщенных является необходимость учитывать взаимодействие наполнителей. В результате искривления поверхности жидкости, вызванного действием сил поверхностного натяжения на границе раздела сред, возникает разность давлений в жидкости и газе , которая называется капиллярным давлением . Капиллярное давление может быть представлено как функция степени насыщения . (3) Величина определяет давление газа, необходимое для вытеснения жидкости из пор. Коэффициент , характеризующий распределение размера зерен скелета, обычно принимает значения от 0,2 до 3. - эффективная водонасыщенность определяется следующим образом [29]: (4) где - остаточная водонасыщенность; - остаточная газонасыщенность. Определяющие соотношения для суммарного напряжения вслед за [30] могут быть записаны в виде (5) Здесь обозначено эффективное напряжение и - коэффициент эффективных напряжений, определяющий зависимость сжимаемости упругого скелета от его объемного модуля и объемного модуля его зерен . Скелет считается изотропным и однородным с двумя константами упругого материала: объемный модуль и модуль сдвига . Исходя из этих предположений получена формулировка для напряжения: (6) Кроме того, в допущении о малых градиентах деформации полагаем, что компоненты тензора деформаций твердого тела и смещения связаны соотношением (7) Уравнения баланса массы для твердого тела и наполнителей имеют вид , (8) , (9) (10) где , , - плотности твердотельного скелета, жидкости и газа соответственно, и - компоненты смещений наполнителей, относительных по отношению к скелету. Формулировка предполагает, что растворения газообразного наполнителя в жидкости не происходит, а градиенты плотностей, насыщенности и пористости пренебрежимо малы. Для трехфазного материала можно записать , (11) , (12) . (13) Уравнение (11) записано для скелета, уравнения (12), (13) - для жидкости и газа соответственно. Введена объемная сила и усредненная плотность . Перенос наполнителей определяется законом Дарси. Соответствующие фазовые проницаемости определяются как и , где и - относительные фазовые проницаемости; - внутренняя проницаемость материала; и - вязкости наполнителей. Относительная фазовая проницаемость определяется как отношение фазовой проницаемости к абсолютной. Используем вслед за [31] следующие экспериментально установленные зависимости для относительных проницаемостей в случае газоводяной смеси в порах: , . (14) Система уравнений (8)-(13) полностью описывает поведение частично насыщенного пороупругого континуума. Выберем в качестве базовых функций давление в порах , и перемещение скелета . Для дальнейшего рассмотрения применяется интегральное преобразование Лапласа. Пусть - функция по , = 0 для . Функция преобразования по Лапласу и ее обращение определено следующим образом: , (15) (формула Римана-Меллина), (16) где s - комплексный параметр преобразования; a - вещественное число большее, чем вещественные части всех особенностей . Применяя преобразование Лапласа к уравнениям (12), (13) и исключая относительные смещения и , получаем , (17) где , - потоки наполнителей. В (17) введены следующие обозначения: . (18) Подстановка (17) в уравнения (8)-(11) позволяет получить формулировку , (19) , (20) , (21) где , , , . Уравнения (19)-(21) могут быть записаны в операторной форме: (22) с несамосопряженным оператором , , , , , , , , , , . Применение интегрального преобразования позволяет строить метод граничных интегральных уравнений. 2. Граничные интегральные уравнения и фундаментальные решения Фундаментальное решение находится из следующего операторного соотношения: , (23) где - единичная матрица; - дельта-функция Дирака; . Итоговая матрица фундаментальных решений (24) содержит девять компонент - компоненты фундаментальных решений перемещений и давлений , , а также шесть спаренных компонент. Все компоненты матрицы приведены в [29]. Гранично-интегральное уравнение для динамической пороупругости в области Лапласа можно получить на основе традиционных формулировок: теоремы взаимности или метода взвешенных невязок. Пороупругий оператор является несамосопряженным, а значит, для этих методов требуются разные фундаментальные решения, но в конечном счете оба метода дают одно и то же интегральное уравнение. Интегральное уравнение динамической пороупругости, построенное на основе метода взвешенных невязок, имеет вид , (25) где - матрица сингулярных решений; - поверхностная сила; - поток жидкости; - поток газа; - вектор внешней нормали. Построенные граничные интегральные уравнения являются сингулярными и существуют в смысле Коши. Для нахождения решения граничных интегральных уравнений требуется соответствующая методика их дискретного представления. В работе используется гранично-элементный подход. 3. Гранично-элементная (ГЭ) дискретизация Чтобы ввести ГЭ-дискретизацию, рассмотрим регуляризованное уравнение [32]: где - матрицы фундаментальных и сингулярных решений; - матрица компонент, содержащих особенности; - вектор обобщенных усилий; - вектор обобщенных перемещений. Базовый процесс ГЭ-дискретизации состоит в разбиении поверхности на граничных элементов совокупностью четырехугольных и треугольных восьмиузловых биквадратичных элементов. При этом треугольные элементы рассматриваются как вырожденные четырехугольные элементы, каждый из которых отображается на контрольный элемент (каждый - это либо квадрат , либо треугольник ). Неизвестные граничные поля определяются через узловые значения и в интерполяционных узлах. При этом для расчетного значения параметра s будем иметь следующие выражения граничных перемещений и поверхностных сил внутри элемента Sk: Здесь - функции формы для линейного четырехугольного элемента. Для получения дискретного аналога ГИУ применим метод коллокации. В качестве узлов коллокации будем выбирать узлы аппроксимации исходных граничных функций. В итоге формируется система линейных алгебраических уравнений . Уравнения записаны соответственно в узлах аппроксимации обобщенных граничных перемещений и обобщенных поверхностных сил. , 4. Шаговая схема численного обращения преобразования Лапласа Рассмотрим метод, опирающийся на теорему об интегрировании оригинала, - шаговый метод численного обращения преобразования Лапласа. Рассмотрим следующий интеграл: . (26) Интеграл (26) порождает задачу Коши для обыкновенного дифференциального уравнения Интеграл (26) заменяется квадратурной суммой, весовые множители которой определяются с помощью изображения по Лапласу и линейного многошагового метода [33-36]. Дальнейшее изложение идет с учетом результатов этих работ. Традиционный шаговый метод интегрирования оригинала состоит в том, что интеграл (26) считается по следующему соотношению: (27) где . Единственная аппроксимация, которая используется при выводе формулы (27), заключается в применении линейного многошагового метода для решения задачи Коши дифференциального уравнения первого порядка. Все остальные вычисления проводятся непосредственно. Многошаговый метод должен быть порядка точности , являясь строго ноль-устойчивым или A-устойчивым, а функция ограничена в правой полуплоскости относительно прямой , то есть при Если функция аналитична и ограничена в области , где , критерий устойчивости может быть ослаблен до A(a)-устойчивости. Выражаясь в терминах характеристической функции : функция не должна содержать нулей и полюсов внутри замкнутого единичного круга за исключением однократного нуля при где при , где при К соответствующим примерам многошаговых методов относятся методы дифференцирования назад порядка : для A-устойчивого метода дифференцирования назад второго порядка можем записать следующую характеристическую функцию: Следующая аппроксимация, которая используется в методе, - численное интегрирование, при получении значений аппроксимируем с равными шагами : . (28) При условии, что функция в уравнении (28) вычисляется с некоторой погрешностью выбор и допускает погрешность вычисления порядка . 5. Численные эксперименты Рассматривается задача о действии торцевой силы на призматическое пороупругое тело длиной 3 м, высотой и шириной 1 м, с жестко закрепленным концом. Краевая задача представлена на рис. 1. Отклики перемещения и поверхностные усилия , вызванные силой = , наблюдаются соответственно на нагруженном и на закрепленном торцах. Параметры частично насыщенного пористого материала соответствуют параметрам песчаника: = , = , = , = , = , = , = , = , = , = , = , = , = 0,9, = 0, = 1, = 1,5. На временном интервале от 0 с до 0,035 с шаг по времени выбирается из условия L = N, параметр R принят равным 0,997. Задача решается в трёхмерной постановке. Заданные граничные условия позволяют сравнивать решение исходной задачи с решением соответствующей одномерной задачи, представленным в [21]. При решении использовались три гранично-элементные сетки с различной степенью пространственной дискретизации (рис. 2). Наличие двух плоскостей симметрии дает возможность рассматривать только часть сетки, при этом четверть сетки a содержит 14 элементов, сетки б - 56 элементов, сетки в - 126 элементов. Рис. 1. Геометрическая модель и граничные условия задачи о частично насыщенной пороупругой консоли Fig. 1. Geometry and boundary conditions of a partially saturated poroelastic column а б в Рис. 2. Гранично-элементные сетки: a - 56 элементов; б - 224 элемента; в - 504 элемента Fig. 2. Boundary-element meshes: a - 56 elements; b - 224 elements; c - 504 elements Рис. 3. Гранично-элементные решения, полученные на трех различных сетках: a - перемещения; б - усилия Fig. 3. Boundary-element solutions on the three different meshes for: a - displacements; b - tractions Рис. 4. Гранично-элементные решения, полученные при трех различных значениях параметра шагового метода: a - перемещения; б - усилия Fig. 4. Boundary-element solutions for the three different values of the time-step method parameter for: a - displacements; b - tractions Результаты исследования сходимости по расчетной сетке показывают, что с ростом числа элементов гранично-элементные решения для перемещения (рис. 3, а) и напряжения (рис. 3, б) принимают значения, близкие к рассчитанным по аналитическим формулам. Варьирование параметром шаговой схемы приводит к тому, что при снижении его до значения 0,00002185c в расчетах на сетке c наблюдается незначительное уточнение амплитуды решения (рис. 4, а, б), однако при дальнейшем его уменьшении решение расходится. Заключение Дано описание математической модели частично насыщенной пористой среды, основанной на модели Био, с пятью базовыми функциями. В изображениях по Лапласу приведена математическая модель краевой задачи трехмерной динамической теории пороупругости. Представлен шаговый метод численного обращения преобразования Лапласа. Дано описание методики гранично-элементного моделирования. Решена задача о действии нагрузки на торец призматического частично насыщенного пороупругого тела. Полученные гранично-элементные и рассчитанные аналитически значения искомых величин графически неразличимы.

### L A Igumnov

Research Institute of Mechanics, National Research Lobachevsky State University of Nizhny Novgorod

Email: igumnov@mech.unn.ru

### A N Petrov

Don State Technical University

# References

1. Frenkel J. On the theory of seismic and seismoelectric phenomena in a moist soil // Journal of Physics. - 1944. - Vol. 8. - P. 230-241.
2. Косачевский Л.Я. О распространении упругих волн в двухкомпонентных средах // Прикладная математика и механика. - 1959. - Т. 23, № 6. - С. 1115-1123.
3. Biot M.A. Theory of propagation of elastic waves in a fluid-saturatedporous solid. I. Low-frequency range // Journal of the Acoustical Society of America. - 1956. - Vol. 28. - P. 168-178. doi: 10.1121/1.1908239
4. Biot M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher-frequency range // Journal of the Acoustical Society of America. - 1956. - Vol. 28. - P. 179-191. doi: 10.1121/1.1908241
5. Brutsaert W. The propagation of elastic waves in unconsolidated unsaturated granular mediums // Journal of Geophysical Research. - 1964. - Vol. 69. - P. 243-257. doi: 10.1029/JZ069i002p00243
6. Nikolaevskij V.N. Mechanics of porous and fractured media. Singapore, World Scientific, 1990. - 472 p.
7. Garg S.K., Nayfeh A.H. Compressional wave propagation in liquid and/or gas saturated elastic porous media // Journal of Applied Physics. - 1986. - Vol. 60. - P. 3045-3055. doi: 10.1063/1.337760
8. Santos J.E., Corberó J.M., Douglas J.Jr. Static and dynamic behavior of a porous solid saturated by a two-phase fluid // Journal of the Acoustical Society of America. - 1990. - Vol. 87. - P. 1428-1438. doi: 10.1121/1.399439
9. A model for wave propagation in a porous medium saturated by a two-phase fluid / J.E. Santos, J.Jr. Douglas, J. Corberó, O.A. Lovera // Journal of the Acoustical Society of America. - 1990. - Vol. 87. - P. 1439-1448. doi: 10.1121/1.399440
10. Tuncay K., Corapcioglu M.Y. Body waves in poroelastic media saturated by two immiscible fluids // Journal of Geophysical Research. - 1996. - Vol. 101. - P. 25149-25159. doi: 10.1029/96JB02297
11. Tuncay K., Corapcioglu M.Y. Wave propagation in poroelastic media saturated by two fluids // Journal of Applied Mechanics. - 1997. - Vol. 64. - P. 313-320. doi: 10.1115/1.2787309
12. Wei C., Muraleetharan K.K. A continuum theory of porous media saturated by multiple immiscible fluids: I. Linear poroelasticity // International Journal of Engineering Science. - 2002. - Vol. 40. - P. 1807-1833. doi: 10.1016/S0020-7225(02)00068-X
13. Lo W.-C., Sposito G., Majer E. Wave propagation through elastic porous media containing two immiscible fluids // Water Resources Research. - 2005. - Vol. 41. - P. 140-160. doi: 10.1029/2004WR003162
14. Lu J.-F., Hanyga A., Jeng D.-S. A mixture-theory-based dynamic model for a porous medium saturated by two immiscible fluids // Journal of Applied Geophysics. - 2007. - Vol. 62. - P. 89-106. doi: 10.1016/j.jappgeo.2006.08.002
15. Berryman J.G. Waves in partially saturated porous media / Eds. In W.E. Fitzgibbon, M.F. Wheeler // Wave Propagation and Inversion (Proceedings of SIAM Conference on Mathematical and Computational Issues in Geophysical Fluid and Solid Mechanics) SIAM, Philadelphia, 1992. - P. 1-25.
16. Smeulders D., De La Rosette J., Van Dongen M. Waves in partially saturated porous media // Transp. Porous. Med. - 1992. - Vol. 9. - P. 25-37.
17. Thomas S.D. A finite element model for the analysis of wave induced stresses, displacements and pore pressures in an unsaturated seabed I // Theory. Comput. Geotech. - 1989. - Vol. 8. - P. 1-38.
18. Thomas S.D. A finite element model for the analysis of wave induced stresses, displacements and pore pressures in an unsaturated seabed II // Model verification. Comput. Geotech. - 1995. - Vol. 17. - P. 107-132.
19. Meissner H., Becker A. Dynamic behaviour of partially saturated sand // T. Built Environ. - 1995. - Vol. 14. - P. 45-55.
20. Schanz M. Poroelastodynamics: Linear models, analytical solutions, and numerical methods // Appl. Mech. Re. - 2009. - Vol. 62 (3). - 15 p.
21. Li P., Schanz M. Wave propagation in a 1-d partially saturated poroelastic column // Geophys. J. Int. - 2011. - Vol. 184. - P. 1341-1353.
22. Gatmiri B., Jabbari E. Time-domain green’s function for unsaturated soils. Part I: two-dimensional solution // International Journal of Solids and Structures. - 2005. - Vol. 42. - P. 5971-5990.
23. Gatmiri B., Jabbari E. Time-domain Green’s function for unsaturated soils. Part II: three-dimensional solution // International Journal of Solids and Structures. - 2005. - Vol. 42. -P. 5991-6002.
24. Jabbari E., Gatmiri B. Thermo-Poro-Elastostatic Green Functions for Unsaturated Soils // Comput. Model. Eng. Sci. - 2008. - Vol. 18 (1). - P. 31-43.
25. Gatmiri B., Maghoul P., Duhamel, D. Two-dimensional transient thermo-hydro-mechanical fundamental solutions of multiphase porous media in frequency and time domains // Int. J. Sol. Struct. - 2010. - Vol. 47. - P. 595-610. doi: 10.1016/j.ijsolstr.2009.10.022
26. Maghoul P., Gatmiri B., Duhamel D. Three-dimensional transient thermo-hydro-mechanical of unsaturated soils // Int. J. Numer. Anal. Methods Geomech. - 2009. - Vol. 34. - P. 297-329. doi: 10.1002/nag.820
27. Maghoul P., Gatmiri B., Duhamel D. Boundary integral formulation and two-dimensional fundamental solutions for dynamic behaviour analysis of unsaturated soils // International Journal of Soil Dynamics and Earthquake Engineering. - 2011. - Vol. 31 (11). - P. 1480-1495.
28. Analytical 3d transient elastodynamic fundamental solution of unsaturated soils / I. Ashayeri, M. Kamalian, M.K. Jafari, B. Gatmiri // Int. J. Numer. Anal. Methods Geomech. - 2010. - Vol. 35. - P. 1801-1829.
29. Li P. Boundary element method for wave propagation in partially saturated poroelastic continua. Computation in engineering and science. Graz: Verlag der Technischen Universität, 2012.
30. Bishop A.W., Blight G.E. Some aspects of effective stress in saturated and partly saturated soils // Geotechnique. - 1963. - Vol. 13. - P. 177-197.
31. Brooks R.H., Corey A.T. Hydraulic properties of porous media. - In Hydraulic Papers. Colorado State Univ., 1964.
32. Баженов В.Г., Игумнов Л.А. Методы граничных интегральных уравнений и граничных элементов в решении задач трехмерной динамической теории упругости с сопряженными полями. - М.: Физматлит, 2008. - 352 с.
33. Lubich C. Convolution quadrature and discretized operational calculus // Numer. Math. I. - 1988. - Vol. 52 (2). - P. 129-145.
34. Lubich C. Convolution quadrature and discretized operational calculus // Numer. Math. II. - 1988. - Vol. 52 (4). - P. 413-425.
35. Schanz M. Wave Propagation in Viscoelastic and Poroelastic Continua. - Berlin Springer, 2001. - 170 p.
36. Литвинчук С.Ю., Петров А.Н., Сабаева Т.А., Решение задачи о действии давления внутри полости в пороупругом полупространстве // Проблемы прочности и пластичности. - 2014. - № 76 (2). - C. 97-105.

# Statistics

#### Views

Abstract - 287

PDF (Russian) - 77

#### PlumX

Copyright (c) 2016 Igumnov L.A., Petrov A.N.