Электронный научный журнал
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,791

ГЕОМЕХАНИЧЕСКАЯ МОДЕЛЬ БЛОКОВОГО ДВИЖЕНИЯ ЗЕМНОЙ КОРЫ

Трофименко С.В. 1 Гриб Н.Н. 1
1 Технический институт (филиал) ФГАОУ ВПО "Северо-Восточный федеральный университет имени М.К. Аммосова", г. Нерюнгри, Республика Саха (Якутия)
Рассмотрена задача движения блока земной коры для случая периодической силы трения. В системе блоков с упругим взаимодействием и наличием одного дезинтегрированного контакта (разлома) задача взаимодействия системы блоков сводится к одномерному уравнению движения блока в периодическом потенциале внешних сил. Приведение уравнения движения блока к уравнению в обобщенных безразмерных координатах позволило сопоставить задачу о движении блока с известной задачей о колебании нелинейного математического маятника. Сравнение уравнения математического маятника и полученное уравнение движения блока показало, что движение блока происходит в виде затухающих колебаний, когда затухание пропорционально первой степени скорости. Переход колебательного движения блоков в поступательное (в режим «stick-slip») возможен при уменьшении сил трения вследствие изменения физических свойств контакта. Фактически это означает, что длительные слабые периодические воздействия на массив горных пород могут инициировать землетрясение либо модулировать реальный сейсмический процесс.
землетрясение
нелинейная сила трения
периодическая нагрузка
блоковое строение
земная кора
1. Астафуров С. В., Шилько Е. В., Димаки А. В., Ружич В. В., Лопатин В. В., Попов В. Л., Псахье С. Г. Изучение особенностей отклика границ раздела в разломно-блоковых средах на изменение их состояния и динамические воздействия // Известия Томского политехнического университета. – 2005. – Т. 308, № 5. – С. 25-32.
2. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. – М.: Наука, 1976. – 576с.
3. Садовский М. А., Болховитинов Л. Г., Писаренко В. Ф. Деформирование геофизической среды и сейсмический процесс. – М.: Наука, 1987. – 100 с.
4. Трофименко С. В. Тектоническая интерпретация статистической модели распределений азимутов аномалий гравимагнитных полей Алданского щита // Тихоокеанская геология. – 2010. – Т. 29, №3. – С. 64-77.
5. Tomlinson G. A. A molecular theory of friction // Phil. Mag. Series. – 1929. – P. 935-939.
Введение

По результатам статистических исследований аномалий геофизических полей магнитной и гравитационной природы в пределах изучаемой территории установлена блоковая структура земной коры [4] в виде ромбических структур, что отражает общие закономерности строения литосферы [3].

При изучении закономерностей сейсмического процесса Олекмо - Становой зоны было установлено, что в динамике сейсмичности содержатся периодические компоненты. Этот  результат свидетельствует о том, что динамическая функция состояния массива горных пород (энергонасыщенность) должна содержать периодические компоненты вида

ff,

(1)

где g, f - энергетический и периодический параметры сейсмической активности, формирующиеся под действием внешней силовой нагрузки f; f - количество разнородных силовых нагрузок.

Статистический анализ распределений количества сейсмических событий на различных периодах реализаций  показал, что неоднородность сейсмической активности наиболее контрастно проявляется на суточных и недельных периодах. Фактически это означает,  энергонасыщенность массива горных пород (1) можно  аппроксимировать функциями вида

f,

(2)

где k1 = 1,2,..., 24 - дискретное время (часы суток) и аналогично для недельной статистики

f

(3)

где k2 = 1,2,..., 7 - дискретное время (дни недели).

Изменение энергонасыщенности массива горных пород в виде функций распределения (2) - (3) связано с действием внешней силовой нагрузки и, как следствие, с деформацией геофизической среды [3].

Источником суточных аномалий (2) могут быть силы инерции, связанные с неравномерностью суточного вращения Земли, и как следствие, с изменением кинематики взаимодействия блоков земной коры. Для недельных аномалий (3) - причиной может быть семидневная цикличность техногенной нагрузки на геологическую среду посредством массовых взрывов.

Фактически это означает, что необходимо рассмотреть задачи о движении блоков под действием внешней периодической силовой нагрузки в потенциальном поле силы трения. Данные задачи определили цель настоящего исследования, направленную на построение моделей движения блоков земной коры под действием слабых возмущений.

Результаты моделирования

Поставленная задача о движении (относительном) блоков может быть рассмотрена в одномерном случае, если из всего набора шовных межблоковых зон одна отличается по своим физическим свойствам (верхний фрагмент, рис.1). В этом случае, выбирая направление оси x вдоль данного контакта  и проецируя все векторные величины на данное направление, получим одномерную модель.

Если рас­сматривать границу раздела блоков как тривиальный трибологический контакт, то для анализа движения блока в режиме «stick-slip» по контакту можно использовать формализм обобщенной одномерной модели Томлинсона [5], разработанной для изучения процессов трения. В рамках этой мо­дели рассматривается одномерное движение тела с массой m под действием силы F в поле периодической потенциальной силы сопротивления c амплитудой N и периодом a. Уравнение движения тела имеет вид:

g,

(4)

где f - вязкость среды, x - координата тела. Если ввести новые безразмерные переменные r   и f, то f, f, и уравнение (4) может быть переписано как

f,

(5)

где f - безразмерная сила, f - безразмерный коэффициент затухания.

Теперь рассмотрим трибологическую систему (рис.1, нижний фрагмент), состоящую из двух блоков массой m и M, разделенных межблоковой шовной зоной, тонкой пленкой  между блоками (трибологический контакт). Блок m имеет возможность двигаться относительно блока M за счет внешней силы F1 в поле силы F2. Блок M соединен с основанием двумя пружинами с коэффициентом жесткости k. Пусть сила F1, вызывающая смещение блока m, изменяется по периодическому закону f. Проекции этой силы равны f на ось x и f на направление силы F2, где γ - угол между направлением действия силы F1 и осью x.  Координату его центра тяжести f относительно неподвижной системы координат и скорость движения блока f можно записать в виде f, f. В этом случае сила трения Fтр будет зависеть от нагрузки N( F2, F1) нелинейным образом (периодический потенциал).

Действительно, если упруго-вязкие свойства контакта не меняются за период действия силы F1, то силу трения можно представить в виде суммы  f, где f - упругая и f - пластическая составляющие силы трения. Пусть f - коэффициент трения в отсутствие периодической силы,  тогда силу df можно представить в виде

d

где f периодическая компонента коэффициента трения. Сила d пропорциональна скорости движения блока в первой степени f, где η - вязкость среды контакта.

pic 

Рис. 1. Блоковая (а) и трибологическая (б) модели движения блока. Обозначения в тексте

Если контакт не однороден, то дополнительно к силе f появится сила трения f, пропорциональная линейному размеру неоднородности b=L/2n, где L - линейный размер блока, n=1, 2, 3,... .

f,   f  .

(6)

Под действием силы f блок М начнет движение. Уравнение движения блока будет зависеть от силы горизонтального сжатия F2 (направление силы выбрано условно), которую за период действия силы F1 можно принять постоянной и сил упругости f при взаимодействии с соседними блоками.

Таким образом, уравне­ние движения (4) с учетом (6) можно записать в виде

d  .

(7)

Если блоки расположены горизонтально, то роль нагрузки играет нормальная проекция силы горизонтального сжатия . Отличие двух положений блоков заключается в том, что при вертикальном расположении Fтр N(m). В горизонтальном положении сила трения не зависит от массы подвижного блока, и сила трения может принимать любое значение (в зависимости от величины ).

Ограничением горизонтального смещения без проскальзывания служат максимумы f в точках ±f  (см. рис.1). При дальнейшем увеличении x движение переходит в режим «stick-slip», то есть блок из одной потенциальной ямы переходит в другую. При этом величина смещения f будет пропорциональна (± b/4), т.е. f, где k = 0,1,2,... Этот режим может наблюдаться и при изменении физических параметров контакта, как это показано в работе [1] при изучении движения блоков в режиме «stick-slip», в пределах f изменения силы трения.

Уравнение (7) обобщает уравнение (4) и эквивалентно ему при отсутствии сил дальнодействия d и периодической нагрузки d. Сравнивая уравнение математического маятника [2, стр. 487] и уравнение (5), можно заключить, что уравнение (5) описывает затухающие колебания, когда затухание пропорционально первой степени скорости. Переход колебательного движения блоков в поступательное (в режим «stick-slip»)  возможен при уменьшении сил трения (6) вследствие изменения физических свойств контакта.

Рецензенты:

Максимов Евгений Петрович, доктор геолого-минералогических наук, профессор по кафедре ТИТР, Технический институт (филиал) ФГАОУ ВПО "Северо-Восточный федеральный университет имени М.К. Аммосова", г. Нерюнгри.

Заровняев Борис Николаевич, доктор технических наук, профессор, декан Горного факультета Северо-Восточного федерального университета, г. Якутск.


Библиографическая ссылка

Трофименко С.В., Гриб Н.Н. ГЕОМЕХАНИЧЕСКАЯ МОДЕЛЬ БЛОКОВОГО ДВИЖЕНИЯ ЗЕМНОЙ КОРЫ // Современные проблемы науки и образования. – 2012. – № 3.;
URL: http://www.science-education.ru/ru/article/view?id=6462 (дата обращения: 24.09.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074