Об одном методе построения решений в динамических задачах для многослойного пороупругого полупространства
Аннотация
Дата поступления статьи: 14.12.2013Рассмотрено представление решений задачи об установившихся гармонических колебаниях для многослойного пороупругого полупространства в виде суперпозиции соответствующих решений для однородных пороупругих сред. Это позволяет с единых позиций исследовать влияние отдельных слоев структуры на формирование полей смещений и напряжений, а также дать анализ потоков энергии с разделением по отражающим границам и типам волн.
Ключевые слова: многослойное полупространство, установившиеся колебания, среда Био, принцип суперпозиции
При анализе динамического отклика геологической структуры на тарированное воздействие или сейсмические нагрузки одной из основных задач является создание адекватных математических моделей такой структуры. При этом на практике встречается достаточно большое разнообразие типов строения как оползневых поверхностных массивов, так и подстилающих сред. В последнем случае наиболее распространенной моделью является модель слоистого полупространства [1-3]. При этом наличие обводненности в целом или части слоев существенно сказывается на формировании спектрального состава отклика среды, регистрируемого на его поверхности или в оползневой зоне.
Рассматривается в декартовой системе координат 
-слойная пороупругая полуплоскость, свойства каждой составляющей которой определяются [4]: 
- соответственно плотностями твердой и жидкой фаз; 
- коэффициентом массовой связи; 
- упругими постоянными скелета; 
- объемной пористостью; 
- коэффициентом проницаемости; 
- динамической вязкостью.
Движение среды описывается амплитудными функциями смещения скелета (твердой фазы) - 
,
и жидкости в порах - 
,
удовлетворяющих уравнениям Био-Френкеля [5]:
(1)
На поверхности 
задан вектор напряжений
, (2)
здесь 
- напряжения в скелете, 
- сила, действующая на жидкость в порах.
На границах раздела слоев могут быть заданы следующие условия сопряжения сред (открытые поры):
- непрерывности векторов перемещений, раздельно для твердой и жидкой фракций;
- непрерывности векторов напряжений 
.
Отметим, что при отсутствии жидкости, в предположении, что 
, уравнения (1) переходят в уравнения Ламе для изотропных упругих сред.
При решении задачи для однородной полуплоскости с условиями на поверхности 
(3)
функции перемещений в скелете и жидкости можно представляются через волновые потенциалы [6]:
,
а решение краевой задачи отыскивается в пространстве преобразований Фурье [7]:
. (4)
Контур 
при интегрировании в плоскости комплексного параметра 
определяется принципом предельного поглощения [8]. Однако при ненулевой динамической вязкости жидкости в порах 
данный контур может совпадать с действительной осью.
В итоге решение может быть сформулировано в виде следующих матричных соотношений:
(5)
, (6)
где 
Перемещения для точек жидкой фракции осуществляется заменой в формулах для 
функций 
на 
и 
на 
.
Компоненты вектора 
определяются из граничного условия (3) в виде:
.
В итоге получим
.
При исследовании динамических задач для многослойной пороупругой среды требуется построение решений уравнений движения Био-Френкеля для слоя фиксированной толщины, а также выполнение условий сопряжения слоев между собой.
Рассмотрим слой с номером 
, занимающий в локальной системе координат область 
На гранях слоя считаем заданными векторы обобщенных напряжений
, 
. (7)
Для данного слоя решение будем отыскивать в виде [9]:
,
соответствующем решениям уравнений (1) для двух полуплоскостей, пересечение которых определяет требуемый слой.
В этом случае функции с индексами (1,2) можно получить из соотношений (5,6) путем замен:
; 
; 
; 
; 
- номер слоя.
Вектор-функции 
в данном представлении являются неизвестными и подлежат определению из условий (7), записанных в преобразованном по Фурье виде. В результате придем к системе линейных алгебраических уравнений, определяющих для каждого фиксированного 
значения трансформант Фурье 
.
, (8)
.
Отметим, что элементы полученной системы (8) 
, 
экспоненциально убывают при 
, что обеспечивает диагональное преобладание элементов и хорошую обусловленность матрицы системы [10].
После их определения восстановление напряженно-деформированного состояния среды осуществляется контурным интегрированием с использованием соотношений (4).
Литература:
- Ляпин, А.А. Механико-математические модели в задачах активной сейсмологии [Текст]: Монография / Ляпин А.А., Селезнев М.Г., Собисевич Л.Е., Собисевич А.Л. – М.: ГНИЦ ПГК, 1999. 294 с.
- Ewing W.M., Jardetzky W.S., Press F. Elastic waves in layered media. –New York etc.: Mc Graw-Hill Book Co., 1957. –380 p.
- Кадомцев М.И., Ляпин А.А., Тимофеев С.И. К вопросам построения эффективных алгоритмов расчета системы «сооружение-грунт» [Электронный ресурс] // «Инженерный вестник Дона», 2012, №1. – Режим доступа: http://ivdon.ru/magazine/archive/n1y2012/719 (доступ свободный) – Загл. с экрана. – Яз. рус.
- Николаевский, В.Н. Механика насыщенных пористых сред [Текст]: Монография / Николаевский В.Н., Басниев К.С., Горбунов А. Т., Зотов Г. А. - М.: Наука, 1970. -336 с.
- Biot M.A. Propagation of elastic waves in a cylindrical bore containing a fluid. / J. Appl. Phys. –1952. –V.23, 2. –P.997-1005.
- Поручиков, В.Б. Методы динамической теории упругости [Текст]: Монография / Поручиков В.Б. -М.: Наука, 1986. –328 с.
- Уфлянд, Я.С. Интегральные преобразования в задачах теории упругости [Текст]: Монография / Я.С. Уфлянд -Л., Наука. 1967. – 403 c.
- Ворович, И.И. Динамические смешанные задачи теории упругости для неклассических областей [Текст]: Монография / Ворович И.И., Бабешко В.А. - М.: Наука, 1989. –320 с.
- Кадыров Р.Р., Ляпин А.А. Особенности возбуждения слоистых сред внутренними источниками колебаний [Электронный ресурс] // «Инженерный вестник Дона», 2012, №3. – Режим доступа: http://ivdon.ru/magazine/archive/n1y2012/981 (доступ свободный) – Загл. с экрана. – Яз. рус.
- Малышев, А.Н. Введение в вычислительную линейную алгебру [Текст]: Монография / Малышев А.Н. –Новосибирск: Наука, 1991. –229 с.