Perturbation of the stress-strain state of an elastic half-space by the spherical inhomogeneity of elastic properties under shear in horizontal plane taking account of gravitational forces

Abstract


The present paper is devoted to the solution of the problem related to the perturbation of the stress-strain state of an elastic isotropic half-space by the spherical inhomogeneity of elastic properties for the case of only gravity and the case of an additional superimposed pure shear on the half-space. The analytical solution is obtained for the displacement and volume strain components, expressed in terms of the potentials and quasipotentials of the sphere, for the case of only gravity forces acting on the half-spaces and the case of superposition of gravity and pure shear. Using the example of a specific seismic event, it was shown that the contribution of gravitational forces to the perturbation of the stress-strain state of a half-space by the spherical inhomogeneity of elastic properties (the distribution with respect to the volume deformation) is not small. Using the obtained analytical solution for the case of the action of gravity only, it is shown that the volume strain distribution is axisymmetric with respect to the axis . A spherical inhomogeneity of the shear modulus in an isotropic elastic half-space forms a zone of a relative compression located along the strike of the inhomogeneity, and zones of relative stretching above and below the inhomogeneity. A distinctive feature of the zones of relative stretching is their limiting size, while the size of the zone of a relative compression increases when the volume strain decreases. In this case, between these planes, the solution in the neighboring octants differs only in sign. When adding a term associated with the gravitational forces in the solution, the symmetry is not violated but the equality in magnitude of the volume deformations in neighboring octants is violated. Dimensions and geometry of isosurfaces of volumetric deformation varies significantly when gravity is taken into account. Depending on the direction of the axis of the maximum extension of the half-space with the inhomogeneity, the size of the zones of the relative compression increases and the corresponding reduction in the size of the zones of the relative stretching changes with the shape of both.

Full Text

Введение Пионерские работы Эшелби о влиянии включений и неоднородностей упругих свойств на напряженно-деформированное состояние линейно-упругого пространства, опубликованные в середине прошлого века [1, 2], положили начало активному использованию предложенного подхода для решения различных задач механики структурно-неоднородных сред [3-7]. В настоящее время в ряде работ предложено обобщение задач о неоднородности и включении на случай анизотропных сред [8-11] и случай неупругих сред [12-17]. Отдельным направлением использования результатов решения задач о включении и неоднородности являются задачи расчета предвестниковых изменений геофизических полей различной природы при подготовке тектонических коровых (внутриплитных) землетрясений. В таких задачах за неоднородность упругих свойств принимается область геосреды, которая возмущает асейсмичное течение тектонических процессов и является непосредственной причиной возникновения и развития предвестников. Момент распада неоднородности при этом ассоциируется с моментом наступления землетрясения и появлением очага. Одним из возможных физических механизмов появления таких неоднородностей в геосреде служит зацепление двух блоков (отдельностей), которое в поле приложенных тектонических напряжений приводит к появлению консолидированной области (неоднородности), осредненные свойства которой отличаются от свойств среды до зацепления. В таком случае распаду неоднородности соответствует подвижка по границе этих блоков, приводящая к генерации сейсмических волн. Приведенное представление о процессе формирования очага является, безусловно, сильно упрощенным и модельным, однако позволяет проводить прямые расчеты эволюции деформационных, фильтрационных, электромагнитных полей, вызванных подготовкой конкретного сейсмического события, а также оценивать размеры зоны проявления предвестников различной природы. В работах [18-20], посвященных задачам расчета предвестниковых изменений при подготовке и реализации одиночного землетрясения, принимается гипотеза о существенно незначительном вкладе гравитационных сил в возмущение напряженно-деформированного состояния геосреды, вызванное неоднородностью упругих свойств (областью формирования будущего очага). Правомочность такой гипотезы обосновывается тем, что коровые (внутриплитные) землетрясения являются в основном неглубокофокусными (глубина расположения гипоцентра не более 50 км) и гравитационные напряжения на таких глубинах существенно меньше тектонических напряжений, вызванных движением плит. Настоящая работа посвящена оценке вклада гравитационных сил (силы тяжести) в возмущение напряженно-деформированного состояния упругого полупространства, вызванное шаровой неоднородностью упругих свойств, расположенной на некоторой глубине. Полученное аналитическое решение для определения возмущения напряженно-деформированного состояния (НДС) упругого полупространства c шаровой неоднородностью упругих свойств при действии только гравитационных сил сравнивается с решением, полученным в [20] для случая чистого сдвига без учета сил гравитации; решение [20] получено в развитие подхода, основанного на методе функций Грина для однородной упругой среды и использованного ранее в [21-23] для вычисления полей возмущений (пульсаций) деформационных полей в неограниченных композитных средах. 1. Математическая постановка задачи В работе задача о неоднородности упругих свойств формулируется и решается на основе подхода, предложенного в [18, 20]. Рассмотрим исходно однородное упругое полупространство и полупространство с неоднородностью упругих свойств, занимающую область центр которой расположен на некоторой глубине . Граница полупространства свободна от напряжений. Краевая задача теории упругости в перемещениях для полупространства с неоднородностью имеет вид (1) где - вектор перемещений; - тензорное поле упругих модулей, компоненты которого обладают следующими свойствами симметрии: , , - компоненты вектора удельной объемной силы; - линейный однородный оператор граничных условий, система координат - декартова с началом координат в центре неоднородности. На границе между неоднородностью и полупространством предполагается выполнение условия непрерывности векторов перемещений и напряжений. Здесь и далее используется общепринятая система индексов: по повторяющимся индексам проводится суммирование, запятая в нижнем индексе обозначает дифференцирование по соответствующей координате, все индексы принимают значения от единицы до трех. Наряду с задачей (1) рассмотрим задачу нахождения напряженно-деформированного состояния полупространства без неоднородности (2) Предполагается, что решение системы (2) для полупространства известно и известна ее функция Грина , которая удовлетворяет системе (3) где - символы Кронекера; - дельта-функция; - номер единичной силы, приложенной в точке и действующей вдоль оси . Будем искать возмущение поля перемещений, вызванное неоднородностью, как разницу между решениями (1) и (2), полагая [24], что (4) где - характеристическая функция неоднородности (внутри и на границе неоднородности , снаружи неоднородности ); - тензорное поле упругих модулей внутри неоднородности; - числовой параметр, отражающий степень неоднородности свойств, . Тогда из (1), (2) с учетом (4) получаем задачу для определения : (5) Рассмотрим случай, когда модули упругости в неоднородности не существенно отличаются от модулей окружающего полупространства. Тогда для решения задачи (5) используем метод малого параметра, изложенный в [24], и будем искать решение в виде ряда (6) Подставив (5) в (6) и раскрывая итоговое соотношения для каждой степени приходим к следующему набору систем при нулевых граничных условиях: (7) Решение системы (7) можно представить через функции Грина (упростив с использованием формулы Гаусса-Остроградского) в виде (8) где - область неоднородности. В результате рекуррентная последовательность (8) позволяет найти компоненты вектора перемещений для конкретных граничных условий, формы и степени неоднородности. 2. Шаровая неоднородность в упругом полупространстве в поле силы тяжести Для задач расчетов предвестниковых изменений геофизических полей различной природы при подготовке тектонического корового землетрясения наибольший интерес представляет дальняя зона возмущений НДС геосреды (расстояние от неоднородности превосходит несколько ее радиусов), в этом случае форма неоднородности слабо влияет на характер возмущений. В связи с этим далее будем рассматривать случай шаровой неоднородности упругих свойств. Из полевых наблюдений за вариациями отношения скоростей продольных и поперечных волн известно, что вариация модулей упругости в районе формирования очага землетрясения не превосходит 20 %, поэтому для нахождения решения задачи (5) ограничимся первым приближением в рекуррентной последовательности (8). Первое приближение для случая однородного изотропного полупространства с неоднородностью согласно (8) имеет вид (9) где - тензор Миндлина, являющийся функцией Грина для однородного изотропного полупространства [25, 26]. В свою очередь, компоненты тензора Миндлина задаются соотношениями (10) где ; - модуль сдвига среды; - коэффициент Пуассона среды. Вспомогательные функции задаются соотношениями (11) Предположим, что модули упругости внутри неоднородности одинаковы для всех точек объема, т.е. , тогда соотношение (9) может быть записано как (12) Случай, когда на полупространство действует только сила тяжести, может быть сведен к случаю однородного напряженно-деформированного состояния невозмущенного полупространства (согласно [26]) соответствующим заданием граничных условий. Тогда соотношение (12) примет вид (13) где (14) При интегрировании компонент тензора Миндлина решение задачи (13) будет содержать линейные комбинации потенциалов (гармонического и бигармонического ) и квазипотенциалов (бигармонического , гармонического , логарифмического , смешанного ) однородных тел, свойства которых в настоящее время хорошо изучены. Свойства, использованные для решения настоящей задачи, представлены в [32-33]. Для этих потенциалов и квазипотенциалов выполняются соотношения (15) где Для случая полупространства введенные потенциалы по отношению к неоднородности могут быть внешними (индекс «е») и внутренними (индекс «i»), тогда как квазипотенциалы могут быть только внешними. Далее рассмотрим случай однородной и изотропной неоднородности, свойства в которой отличаются от свойств окружающего пространства только модулем сдвига. Физические предпосылки такого предположения описаны в [20]. Тогда приращение модулей упругости в неоднородности можно записать в виде (16) где - модуль сдвига окружающей неоднородность среды; - степень отличия модуля сдвига неоднородности . Раскроем соотношение (13) с учетом (16) и запишем его в напряжениях, учтя, что на дневной поверхности касательные напряжения , , равны нулю ( = 0 во всей области в силу осевой симметрии задачи): (17) Вычисление компонент вектора перемещений с использованием соотношения (17) требует конкретизации граничных условий на границах полупространства. В случае действия на полупространство только сил тяжести компоненты тензора равны и ( - ускорение свободного падения, - коэффициент Пуассона). Видно, что для дневной поверхности автоматически выполняется условие отсутствия напряжений. Раскроем выражение (17) с учетом правила дифференцирования потенциалов и квазипотенциалов (27) и соотношений (15). В результате перемещения определяются формулами (18) где . Дифференцирование соотношений (18) по соответствующим координатам позволяет найти все компоненты тензора деформаций. В задачах расчета гидрогеологических и электрокинетических предвестников землетрясений определяющим фактором является изменение объемной деформации среды, влияющее на интенсивность и направление фильтрационных процессов. После дифференцирования и сложения соответствующих компонент объемная деформация геосреды с неоднородностью упругих свойств может быть найдена как (19) Преимуществом используемого подхода является запись итогового решения через объемные потенциалы и квазипотенциалы однородных тел. Использование соответствующих потенциалов для конкретной формы неоднородности позволяет сразу найти итоговое решение задачи о возмущении НДС упругого полупространства этой неоднородностью. В случае шаровой неоднородности, центр которой располагается в начале координат, потенциалы и квазипотенциалы определяются следующими выражениями: (20) где - объем неоднородности; - радиус неоднородности. Подставляя (20) в (19) с учетом правила переноса потенциалов и квазипотенциалов (26), получаем (21) где , , , выражения с индексом «i» имеют место только внутри неоднородности, а с индексом «e» только вне неоднородности. Для графического представления полученного решения воспользуемся параметрами сейсмического события и свойствами земной коры, указанными в [20] и [27, 28]. Сейсмическое событие магнитудой 5,2 произошло в октябре 2009 г. в районе Петропавловск-Камчатского геодинамического полигона, глубина залегания гипоцентра 10 км, радиус неоднородности, оцененный по магнитуде согласно [29], составил 2862 м. Принятые для расчетов свойства среды: плотность 2500 кг/м3, коэффициент Пуассона 0,25, модуль сдвига 20 ГПа, степень неоднородности На рис. 1 представлены изоповерхности объемной деформации для величин , , и . Положительная величина объемной деформации отражает режим относительного сжатия, тогда как отрицательная - режим относительного растяжения. Таким образом, из полученного решения видно, что распределение объемной деформации является осесимметричным относительно оси z (ввиду осевой симметрии задачи). Симметрия решения относительно плоскости , проходящей через центр неоднородности, отсутствует, что вызвано дневной поверхностью, свободной от напряжений. Вокруг шаровой неоднородности упругих свойств в изотропном упругом полупространстве формируются зона относительного сжатия (обозначена темно-серым цветом на рисунке), расположенная вокруг неоднородности, и зоны относительного растяжения (светло-серый цвет), расположенные над и под неоднородностью. Отличительной особенностью зон относительного растяжения является их предельный размер, тогда как размер зоны относительного сжатия растет при уменьшении величины объемной деформации. Наличие зоны относительного растяжения над гипоцентром сейсмического события было предсказано ранее в работе [30]. Необходимо отметить, что размеры возмущения НДС упругого полупространства шаровой неоднородностью упругих свойств, совпадающие с размерами области подготовки землетрясения или зоной проявления предвестников (которые согласно [18] определяются деформацией ), в случае действия только гравитационных сил достигают 190 км. Этот факт позволяет уже утверждать, что вклад гравитационных сил в общее возмущения НДС не является пренебрежительно малым. Рис. 1. Изоповерхности объемной деформации, вычисленной согласно (21), для величин деформации (а); (б); (в); (г) (положительная деформация (относительное сжатие) - тёмно-серый цвет, отрицательная деформация (относительное растяжение) - светло-серый цвет) Fig. 1. The isosurfaces of volumetric deformation, calculated according to (21), for the strain values (a); (b); (c); (d) (positive deformation (relative compression) - is shown in a dark gray color; negative deformation (relative tension) is shown in a light gray color Полученное решение (18) для случая отличных от нуля осевых компонент тензора является частным случаем решения задачи об упругом полупространстве с шаровой неоднородностью, полученным ранее в работе [20], здесь - тензорное поле, представляющее собой произведение вариаций упругих свойств в полупространстве с неоднородностью и частных производных по пространственным координатам ее функций Грина. 3. Шаровая неоднородность в упругом полупространстве при сдвиге в горизонтальной плоскости с учетом силы тяжести В случае действия на полупространство, помимо гравитационных сил, сил тектонической природы, вызванных движением и взаимодействием плит, решение задачи (13) сводится к дополнительному учету граничных условий для горизонтальных напряжений на бесконечности. Рассмотрим случай сдвига в горизонтальной плоскости, который для задач геодинамики и тектонофизики является широко распространенным. В случае сдвига в горизонтальной плоскости граничные условия могут быть записаны как (22) где - касательные напряжения на бесконечности. Тогда из соотношения (13) получаем (23) Видно, что решение в этом случае складывается из двух решений: отдельно для гравитационных сил (первое слагаемое) и касательных напряжений (второе слагаемое). Решение для случая сдвига в горизонтальной плоскости было рассмотрено ранее в работах [18, 20]. Согласно этому решению выражение для объемной деформации имеет вид (24) где , , , как и ранее, выражение с индексом «e» имеет место только вне неоднородности. Таким образом, решение задачи о возмущении напряженно-деформированного состояния упругого полупространства шаровой неоднородностью упругих свойств при сдвиге в горизонтальной плоскости с учетом сил тяжести задается суммой решений (21) и (24). Для оценки вклада гравитационных сил в общее возмущение НДС при сдвиге в горизонтальной плоскости на рис. 2 представлены изоповерхности объемной деформации в случае учета и неучета гравитации. В случае только сдвига в горизонтальной плоскости видно, что решение симметрично относительно двух биссекторных плоскостей октантов 5-7 и октантов 6-8. При этом между этими плоскостями решение в соседних октантах отличается только знаком. При добавлении в решении слагаемого, связанного с гравитационными силами, симметрия не нарушается, но нарушается равенство по модулю объемных деформаций в соседних октантах. Особенно необходимо подчеркнуть, что размеры и геометрия изоповерхностей объемной деформации при учете сил тяжести существенно меняется. В зависимости от направления оси максимального растяжения полупространства с неоднородностью происходит увеличение размеров зон относительного сжатия и соответствующее уменьшение размеров зон относительного растяжения с изменением формы обеих. Полученное решение для случая сдвига в горизонтальной плоскости с учетом гравитационных сил показало, что гипотеза о пренебрежении силами тяжести при прямом моделировании предвестниковых возмущений, когда зона формирующегося очага землетрясения представляется в виде неоднородности упругих свойств, несостоятельна. Рис. 2. Изоповерхности объемной деформации для величин деформации (а-б); (в-г); (д-е) при чистом сдвиге (а, в, д) и чистом сдвиге с учетом гравитационных сил (б, г, е) (положительная деформация (относительное сжатие) - тёмно-серый цвет, отрицательная деформация (относительное растяжение) - светло-серый цвет) Fig. 2. Surface deformation isosurfaces for strain values (а-b); (c-d); (e-f) under pure shear (а, c, e) and pure shear with the allowance for gravitational forces (b, d, f) (positive deformation (relative compression) is shown in a dark gray color, negative deformation (relative tension) is shown in a light gray color Выводы В работе рассмотрена задача о возмущении напряженно-деформированного состояния упругого полупространства шаровой неоднородностью упругих свойств при сдвиге в горизонтальной плоскости с учетом вклада в решение гравитационных сил. Получено аналитическое решение для компонент вектора перемещений отдельно для случая действия на полупространство только сил тяжести и для случая сдвига в горизонтальной плоскости с учетом сил тяжести. На основе этих решений получены выражения для нахождения объемной деформации полупространства в обоих случаях (соотношения (21) и сумма (21) и (24)). Показано, что в случае действия только сил тяжести объемная деформация среды имеет осевую симметрию (ось симметрии z) и разделяется на зону относительного сжатия и зоны относительного растяжения. Последние имеют предельные размеры и не увеличиваются при уменьшении величины объемной деформации. Для задачи сдвига в горизонтальной плоскости с учетом действия сил гравитации показано, что исходная симметрия распределения объемной деформации при сдвиге в горизонтальной плоскости в случае учета сил тяжести не меняется, существенно меняется форма и размеры зоны относительного сжатия/растяжения в зависимости от направления оси максимального растяжения. Полученные результаты указывают на то, что при моделировании предвестниковых возмущений полей различной природы, когда зона формирующегося очага землетрясения представляется в виде неоднородности упругих свойств, пренебрегать действием гравитационных сил нельзя ввиду их существенного влияния на размер и форму возмущения НДС среды.

About the authors

I A Panteleev

Institute of Continuous Media Mechanics

E V Poltavceva

Institute of Volcanology and Seismology Fat Eastern Brunch of Russian Academy of Science

V A Mubassarova

Research Station of Russian Academy of Science

V A Gavrilov

Institute of Volcanology and Seismology Fat Eastern Brunch of Russian Academy of Science

References

  1. Eshelby J.D. The Determination of the Elastic Field of an Ellipsoidal Inclusion, and Related Problems // Proc. R. Soc. London. - 1957. - A241. - P. 376-396.
  2. Эшелби Дж. Континуальная теория дислокаций. - М.: Иностр. лит., 1963. - 248 с.
  3. Tsuchida E., Mura T. The Stress Field in an Elastic Half Space Having a Spheroidal Inhomogeneity Under All-Around Tension Parallel to the Plane Boundary // J. Appl. Mech. - 1983. - Vol. 50. - No. 4a. - P. 807-816.
  4. Гольдштейн Р.В., Городцов В.А., Устинов К.Б. Влияние поверхностных остаточных напряжений и поверхностной упругости на деформирование шарообразных включений нанометровых размеров в упругой матрице // Физическая мезомеханика. - 2010. - Т. 13, № 5. - С. 127-138.
  5. Гольдштейн Р.В., Шифрин Е.И. Интегральные уравнения задачи об упругом включении. Полное аналитическое решение задачи об эллиптическом включении // Изв. Рос. акад. наук. Механика твердого тела. - 2004. - № 1. - С. 50-76.
  6. Chang-Hung Kuo. Stress disturbances caused by the inhomogeneity in an elastic half-space subjected to contact loading // International Journal of Solids and Structures. - 2007. - Vol. 44 - No. 3-4. - P. 860-873.
  7. Eshelby problem of an arbitrary polygonal inclusion in anisotropic piezoelectric media with quadratic eigenstrains / Y. M. Yue, K. Y. Xu, Q. D. Chen, E. Pan // Acta Mechanica. - 2015. - Vol. 226. - No. 7. - P. 2365-2378.
  8. Withers P.J. The Deformation of the Elastic Field of an Ellipsoidal Inclusion in Transversely Isotropic Medium, and Its Relevance to Composite Materials // Philosophical Magazine. - 1989. - A59. - P. 759-781.
  9. Jiang Dazhi. Viscous inclusions in anisotropic materials: Theoretical development and perspective applications // Tectonophysics. - 2016. - Vol. 693. - P. 116-142.
  10. Семенова Д.В., Устинов К.Б. Некоторые асимптотические разложения компонент тензора Эшелби в случае кубической и гексагональной анизотропии // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2010. - № 4. - С. 80-100.
  11. Березкин В.А. О построении тензора Эшелби для слабоанизотропной упругой среды // Вестн. Моск. ун-та. Сер. 1. Математика, Механика. - 2011. - № 1. - С. 73-76.
  12. Song Zhi-ping, Yin Xiang-chu, Mei Shi-rong. Theoretical analysis of the spatio-temporal evolution of the bulk-strain field based on a rheologic inclusion model // Acta Seismologica Sinica. - 2000. - Vol. 13. - No. 5. - P. 525-535.
  13. Song Zhi-ping, Yin Xiang-chu, Mei Shi-rong. 3-D rheologic model of earthquake preparation (I) - Displacement field // Acta Seismologica Sinica. - 2003. - Vol. 16. - No. 6. - P. 588-597.
  14. Song Zhi-ping, Yin Xiang-chu, Mei Shi-rong. 3-D rheologic model of earthquake preparation (II): Strain field and its applications // Acta Seismologica Sinica. - 2004. - Vol. 17. - No. 2. - P. 131-143.
  15. Song Zhi-ping, Yin Xiang-chu, Mei Shi-rong. 3-D rheologic model of earthquake preparation (III): Precursor field // Acta Seismologica Sinica. - 2006. - Vol. 19. - No. 1. - P. 20-29.
  16. The Displacement and Strain Field of Three-dimensional Rheologic Model of Earthquake Preparation / Zhiping Song, Xiangchu Yin, Shirong Mei, Yucang Wang, Can Yin, Huihui Zhang, Langping Zhang // Pure and Applied geophysics. - 2006. - Vol. 163. - No. 9. - P. 1991-2009.
  17. Цвелодуб И.Ю. О тензоре Эшелби // Прикладная математика и механика. - 2010. - Т. 74, № 2. - C. 346-351.
  18. Добровольский И.П. Распределение деформаций и напряжений при подготовке тектонического землетрясения // Физика Земли. - 2003. - № 10. - С. 37-40.
  19. Добровольский И.П. Гидродинамические явления // Физика Земли. - 2005. - № 1. - С. 89-96.
  20. Добровольский И.П. Математическая теория подготовки и прогноза тектонического землетрясения. - М.: Физматлит, 2009. - 240 с.
  21. Волков С.Д., Ставров В.П. Статистическая механика композитных материалов. - Минск.: Изд-во Белорус. гос. ун-та, 1978. - 208 с.
  22. Хорошун Л.П., Маслов Б.П., Лещенко П.В. Прогнозирование эффективных свойств пьезоактивных композитных материалов. - Киев: Наук. думка, 1989. - 208 с.
  23. Шеpмеpгоp Т.Д. Теоpия упpугости микpонеодноpодных сpед. - М.: Наука, 1976. - 399 с.
  24. Ломакин В.А. Статистические задачи механики твердых деформируемых тел. - М.: Наука, 1970. - 139 с.
  25. Лурье А.И. Теория упругости. - М.: Наука, 1970. - 940 с.
  26. Новацкий В. Теория упругости. - М.: Мир, 1975. - 872 с.
  27. Panteleev I.A., Gavrilov V.A. Implications of electrokinetic processes for the intensity of geoacoustic emission in the time vicinity of a tectonic earthquake: a theoretical study // Russian Journal of Earth Sciences. - 2015. - Vol. 15. - No. 4. - С. 1a-14.
  28. Пантелеев И.А., Полтавцева Е.В., Гаврилов В.А. Изменение НДС в районе измерительной скважины Г-1 на Камчатском геодинамическом полигоне при подготовке близких сильных землетрясений по данным скважинных измерений: численное моделирование // Проблемы геодинамики и геоэкологии внутриконтинентальных орогенов: тез. докл. VII Междунар. симпозиума. К 80-летию со дня рождения выдающегося ученого, основателя и первого директора НС РАН Ю.А. Трапезникова (18.10.1936 - 13.04.1999). - Бишкек, 2017. - С. 257.
  29. Ризниченко Ю.В. Размеры очага корового землетрясения и сейсмический момент // Исследования по физике очага землетрясения. - М.: Наука, 1976. - С. 9-18.
  30. Алексеев А.С., Белоносов А.С., Петренко В.Е. О концепции многодисциплинарного прогноза землетрясений с использованием интегрального предвестника // Проблемы динамики литосферы и сейсмичности: сб. науч. тр. ГЕОС. Вычислительная сейсмология. - 2001. - Вып. 32. - С. 81-97.
  31. Сретенский Л.Н. Теория ньютоновского потенциала. - М.-Л.: ОГИЗ-Гостехиздат, 1946. - 318 с.
  32. Гюнтер Н.М. Теория потенциала и ее применение к основным задачам математической физики. - М.: Гостехиздат, 1953. - 406 с.
  33. Купрадзе В.Д. Методы потенциала в теории упругости. - М.: Физматлит, 1963. - 472 с.

Statistics

Views

Abstract - 319

PDF (Russian) - 69

Cited-By


PlumX


Copyright (c) 2017 Panteleev I.A., Poltavceva E.V., Mubassarova V.A., Gavrilov V.A.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies