Вестник КРАУНЦ. Физ.-мат. науки. 2018. № 5(25). C. 74-82. ISSN 2079-6641
DOI: 10.18454/2079-6641-2018-25-5-74-82
УДК 519.25 +519.237+ 550.34.01
МОДЕЛЬ ПРОЦЕССОВ РЕЛАКСАЦИИ В РАЗЛИЧНЫХ РЕЖИМАХ ПЛАСТИЧЕСКИХ ДЕФОРМАЦИЙ
О.В. Шереметьева1, 2
1 Институт космофизических исследований и распространения радиоволн ДВО РАН, 684034, Камчатский край, с. Паратунка, ул. Мирная, 7
2 Камчатский государственный университет имени Витуса Беринга, 683032, г. Петропавловск-Камчатский край, ул.Пограничная, 4
E-mail: olga.v.sheremetyeva@gmail.com
Рассмотрена модель пяти режимов деформационной активности на основе составного процесса Пуассона и его фрактальных обобщений. Проведено сравнение с зависимостями экспериментально полученными авторами работ [1, 2].
Ключевые слова: функция Миттаг-Леффлёра, масштабные законы, стохастическая модель, пластические деформации, релаксационные процессы.
© Шереметьева О. В., 2018
MSC 60G99
MODEL RELAXATION PROCESSES IN THE DIFFERENT MODES OF
PLASTIC DEFORMATION
O. V. Sheremetyeva1, 2
1 Institute of Cosmophysical Researches and Radio Wave Propagation Far-Eastern Branch, Russian Academy of Sciences, 684034, Kamchatskiy Kray, Paratunka, Mirnaya st., 7, Russia
2 Vitus Bering Kamchatka State University, 683031, Petropavlovsk-Kamchatsky, Pogranichnaya st., 4, Russia
E-mail: olga.v.sheremetyeva@gmail.com
A model of five modes of deformation activity based on the compound Poisson process and its fractal generalizations is considered. Comparison of the analytical dependences with the experimental ones obtained by the authors Carbon and Buck, was carried out [1, 2].
Key words: Mittag-Leffler function, scale laws, stochastic model, plastic deformations, relaxation processes.
© Sheremetyeva O.V., 2018
Введение
В рамках работы [3] обсуждались пять видов деформационного процесса: нормальные (фоновые) пульсации, их замедление (замирания), затем активизация (разрядка), которая имеет три фазы: форшоков, главного удара и афтершоков, - это реализуются те события, которые не состоялись при замирании, и снова фоновый режим. Процесс рассматривался с точки зрения дислокационных изменений, вызванных работой внешних сил, которые поддерживают критический уровень упругих напряжений. Случайные дислокационные изменения в заданном объеме V и временном интервале Т рассматривались как однопараметрический процесс
N
^) = £ Ь(хк, гк), к=0
где хк, гк - координаты, а Ь(хк,¿к) - пространственный масштаб к-ой дислокации.
В первом приближении случайные величины хк, ¿к и Ь(хк,¿к) можно считать независимыми и процесс - пуассоновским. Такое описание неприменимо в случае роста плотности дислокаций, а следовательно, увеличения их связанности и образования дислокационного кластера, в котором плотность дислокаций велика настолько, что перекрываются пространственно-временные радиусы влияния, зависящие от Ь(хк,¿к) [4, 5]. Такие кластеры образуют последовательности форшоков и афтершоков (разреженная упаковка дислокаций), главный удар (плотная упаковка). С этой точки зрения процесс Пуассона может рассматриваться как фоновый режим, а отклонения от него, связанные с локальном упрочнением или ослаблением среды, описываются фрактальным процессом Пуассона.
Данное исследование является логическим продолжением работ [4, 5]. В нём будет рассмотрена математическая модель деформационной активности, построены зависимости, описывающие различные режимы деформационного процесса и проведено сравнение полученных аналитических закономерностей с имеющимися результатами обработки данных сейсмических каталогов.
Математическая модель режимов деформационного процесса
Процессы дислокационных изменений имеют одинаковую структуру на разных масштабных уровнях деформационного процесса, поэтому достаточно рассмотреть процесс пластических деформаций для одного из них.
Пусть зафиксирован пространственный масштаб Ь. Пять режимов деформационной активности можем рассматривать как пять состояний, в которых может находиться деформационный процесс: £1 - фоновый режим (фоновые или нормальные пульсации), £2 - замирания (замедление пульсаций), £3 - режим активизации в фазе форшоков, £4 - режим активизации в фазе главного удара, £5 -режим активизации в фазе афтершоков. Опишем каждое состояние и вероятность его сохранения или перехода в другое состояние используя понятие фрактального процесса Пуассона и функцию Миттаг-Леффлёра - фрактальное обобщение экспоненты
Е(х)= I ПГ+ук) ■
где 0 < V < 1, в случае V = 1 получаем обычную экспоненциальную зависимость, описывающую простейший поток событий (процесс Пуассона).
Релаксационной характеристикой деформационного процесса является фрактальная функция Миттаг-Леффлёра со степенным аргументом
F ( (а, ) v V HMt)V)k
Fv (-(M ' 1 )- E rT+Vk) ■
(1)
выражающая вероятностные характеристики этого процесса. Сама функция (х) даёт описание эффектов памяти (эридитарность процесса, при уменьшении память увеличивается) (рис. 1б), а степенной аргумент (дг)фрактальной функции Миттаг-Леффлёра отвечает за эффекты последействия (нелокальность процесса, процессы ветвления, которые проявляются как фрактальное замедление потока, чем ближе к нулю значение V, тем меньше интенивность потока)(рис. 1а). Вместе (дг)и (х) дают замедленную релаксацию [3].
Таким образом, при варьировании параметров V, V можно получить различные режимы деформационного процесса либо с преобладанием памяти (нелокальности во времени), либо процессов ветвления (пространственной нелокальности). В данном исследовании будем считать параметры равными друг другу V = V.
а) б)
Рис. 1. Изменение характера поведения функции Миттаг-Леффлёра при вариации параметров V, V: а) график I - график функции е-дг ( V = V = 1), график II
- график функции (1) при V = 1, V = 0.9, график III - график функции (1) при V = 1, V = 0.4, б) график I - график функции е-дг ( V = V = 1), график II
- график функции (1) при V = 0.9, V = 1, график III - график функции (1) при V = 0.4, V = 1
Состояние ¿1. В фоновом режиме деформационный процесс является пуассоновским (поток независимых событий) со средней интенсивностью д. Тогда вероятность перехода процесса в следующее состояние зависит от времени экспоненциально и определяется как
Р12 (г ) = 1 - е-дг = 1 - Е1(-д г), а вероятность сохранения процесса в исходном состоянии (рис.2)
Рц(г ) = е-дг = Е1 (-д г).
1
Рис. 2. I - график экспоненциального распределения Рц, II - график распределения вероятности Р22
Состояние £2. Отклонение от фонового режима может произойти, если возникает локальное упрочнение или ослабление среды (изменение реологии). При упрочнении среды возникает область, в которой пластический процесс замедляется, и поток событий ослабевает (область сейсмического затишья или сейсмической бреши [6], в которой формируется деформационная неоднородность). В результате возникновения неоднородности в объёме V исходная интенсивность д потока событий падает. Временной промежуток ожидания каждого последующего события увеличивается, что говорит о появлении памяти (эредитарности) у процесса, которая связана с поведением вероятности сохранения его в исходном состоянии [7, 8] и может быть описана с помощью фрактальной функции Миттаг-Леффлёра с параметрами дг > 1 и 0 < V < 1 (рис.2)
р22(г) = ^(-(дгГ) = £ ГЩ^.
Как видно, вероятность Р22(г) сохранения процесса в исходном состоянии в начальный момент времени убывает сильнее, чем вероятность Рц(г) для потока независимых событий, но с увеличением временного промежутка вероятность Р22(г) меняется медленнее, появляется затянутость процесса во времени (рис.2). Отклонение Р22(г) от обычной экспоненты указывает на наличие эффектов последействия и памяти у фрактального процесса, которые обусловлены тем, что несостоявшиеся события локализованы и образуют кластер связанных между собой событий.
Отметим, что параметр V отвечает за фрактальную размерность распределения событий на временном интервале [9]. Таким образом, упрочнение среды меняет фрактальность процесса (статистику событий), что проявляется как замедленная релаксация, обусловленная эффектами памяти, которые возникли в результате упрочнения.
Состояние £3. Несостоявшиеся в затишье события можно рассматривать как кластер, энергия которого приводит к росту упругих напряжений, в результате чего упрочнение среды преодолевается, а накопленная дополнительная упругая энергия высвобождается. Эту активизацию можно рассматривать как событие более высокого масштаба по отношению к тем событиям, которые не состоялись во время затишья, вызванного локальным упрочнением среды, и энергетика деформационного
возмущения должна соответствовать сумме энергий несостоявшихся событий. Результатом накопления упругой энергии в состоянии ¿2 является возникающее в некоторый момент времени г\ деформационное возмущение и переход в состояние ¿з - фазу форшоков. Для аналитического описания состояния ¿з воспользуемся вероятностью возникновения главного удара в момент времени г *, т. е. вероятностью Р34 перехода из состояния ¿з в состояние 54, которую зададим также с помощью возрастающей фрактальной функции Миттаг-Леффлёра следующим образом (рис.За):
Р34(г ) = Е* (-[-Д (г - г ')]") = £ (~ -к'/' •
где Д - средний поток событий в деформационном возмущении, 0 < V < 1. Тогда вероятностью сохранения состояния ¿з (рис. За), т.е. отсутствия главного удара, будет
Рзз(г) = 1 - Рз4(г) = 1 - Еу (-[-Д (г - г')]у).
р(1)
' ' -
V / ___
4
\
а)
б)
Рис. 3. а) I - график распределения вероятности Рз4, II - график распределения вероятности Р33, б) график распределения плотности потока форшоков (обратный закон Омори)
При дифференцировании выражения Р34 получается распределение нарастающего потока форшоков, которое можно интерпретировать как обратный закон Омори (рис.Зб)
dP34(t) dFv(-[-Mi(, -,*)]v) - (-1)k+1 Mvk(-j(, -,*))
* N Nvk— 1
d,
d,
I
k-0
Г(1 + vk)
Состояние ¿5. По достижении функцией Рз4(г) значения единица процесс переходит в состояние ¿4. Фаза афтершоков может быть определена аналогичным образом, используя выражение для вероятности сохранения главного удара (рис.4а), т. е. вероятности сохранения состояния ¿4,
P44(,)- Fv'(-[Mi(, -,*)]v)- I
k-0
(-[ü (, -, * )]v) Г(1 + v' k)
' \ k
где 0 < V7 < 1.
Дифференцируя это выражение, получим плотность убывающего потока афтершоков, т.е. закон Омори (рис.4б),
dP44(,) dFv/ (-[j(, -,*)]v)
d,
d,
I
k-0
(-1)k jü v'k(jü (, -, *)) Г(1 + v'k)
*\w'k- 1
В таком случае, вероятность не сохранения главного удара (рис.4а), т.е. перехода из состояния в состояние £5 к фазе афтершоков, будет
Р45(г ) = 1 - Р44(г).
а)
б)
Рис. 4. а) I - график распределения вероятности Р44, II - график распределения плотности вероятности Р45, б) график распределения убывающего потока афтершоков (закон Омори)
Сравнение с зависимостями, полученными по данным сейсмических каталогов
На основании построенной модели аппроксимируем зависимости, полученные авторами статей [1, 2]. В них исследуются полные каталоги сейсмических событий [1, 2] и афтершоковые последовательности [1]. В исследовании [1] показано, что масштабный закон для афтершоковых последовательностей отличается от аналогичного закона для последовательностей главных событий из общих каталогов. Как считают авторы, это связано с тем, что даже слабые корреляции между афтершоками сохраняются в течение длительного времени в одном и том же регионе, и это существенно отличается от того, что наблюдается для некоррелированных основных сейсмических событий (главных ударов) [2]. Полученные распределения плотности потока событий представлены в двойном логарифмическом масштабе, множества А и В взяты из полных каталогов главных событий, а множества С и В - это афтершоковые последовательности.
Для аппроксимации использовалась плотность вероятности Р45 не сохранения главного удара (2).
На масштабный закон, полученный авторами [1] для афтершоковых последовательностей, достаточно хорошо ложится кривая с параметрами д = 0.75 и V = 0.5 (рис.5а). Отметим, что хорошее приближение получается и при помощи плотности вероятности, заданной функцией (1), и в этом случае наилучшее приближение получено для параметров V = 0.7 и V = 0.4. Масштабный закон, полученный на основе последовательностей из общих каталогов, аппроксимировать при одинаковых значениях V и V не удалось. График, близкий к экспериментальному закону, получен для значений парметров д = 2, v = 1.15 и V = 0.4 (рис.5б).
В работе [2] исследовались только последовательности событий из общего каталога землетрясений. Авторами были выбраны пять пространтсвенных масштабов разложения и три ограничения на энергии сейсмических событий.
10"
10
10"
10
10
10
v
• • у** * */ ^ ш ■ • V. 1
• Set] C^t t 3 ■ч i ... ■
о Set D i
10" 104 10" 10'
bT D
x = cTS L а)
10
10
10-
II 10
x
10
102 10
10
1 1 i
aBlSl^
V «»
: ■ SetA i \
• Set В
- и Set С w
о Set D i i 1
10"6 104 10 10"
bT D
x = cTS L
б)
10 10
Рис. 5. а) I - аппроксимация масштабного закона для афтершоковых последовательностей (график синего цвета), б) аппроксимация масштабного закона для последовательностей главных событий (график синего цвета, асимптотика выделена красным цветом)
Полученные экспериментальные законы масштабирования аппроксимированы с использованием фрактальной функциии Миттаг-Леффлёра (1) для значений парметров равных ¡1 = 100, V = 1.3 и V = 0.2 (рис.ба) и ¡1 = 500, V = 0.85 и V = 0.2 (рис.бб). Изменяя параметр V и масштабный коэффициент ¡1 могут быть получены графики для остальных случаев.
с Т 8"" ь"1 с Т 8"" I/'1
а) б)
Рис. 6. Аппроксимация масштабного закона для последовательностей главных событий [2] (графики красного и синего цветов)
Обсуждение и выводы
В работе деформационный процесс рассмотрен с вероятностной точки зрения как переход из одного сотояния (режима) в другое. Предложенная аппроксимация с использованием фрактальной функции Миттаг-Леффлёра учитывает историю процесса (свойства нелокальности и во времени, и в пространстве) и является логическим продолжением и расширением модели пуассоновского процесса с независимыми событиями. Кроме того, фрактальный параметр v определяется параметрами среды, что открывает возможности для более полного описания её свойств.
На данный момент хорошо исследованными являются закономерности в последовательностях афтершоков. Проведённое сравнение зависимостей, полученных авторами исследований [1, 2] в результате обработки данных сейсмических каталогов, с модельными показало, что варьирование фрактальных параметров и масштабного коэффициента позволяет достаточно хорошо аппроксимировать эти зависимости с помощью функции Миттаг-Леффлёра. Отметим, что в случае афтершоковых последовательностей аппроксимация функцией (2) достигалась при близких значениях параметров v и V. Для полных каталогов аппроксимация получена только при значительно отличающихся значениях фрактальных параметров.
Список литературы
[1] Carbone V., Sorriso-Valvo L., Harabaglia P., Guerra I., "Unified scaling law for waiting times between seismic events", Europhys. Lett., 6:71 (2005), 1036-104210.
[2] Bak P., Christensen K., Danon L., Scanlon T., "Unified Scaling Law for Earthquakes", Phys. Rev. Lett, 88:17. (2002), 178501-1-178501-4..
[3] Shevtsov B., Sheremetyeva O., "Fractional models of seismoacoustic and electromagnetic activity", E3S Web of Conferences, 20 (2017), 02013. DOI: 10.1051/e3sconf/20172002013.
[4] Шевцов Б. М., Сагитова Р. Н., "Статистический анализ сейсмических процессов на основе диффузионного подхода", ДАН, 426:2 (2009), 254-256. [Shevtsov B.M., Sagito-va R. N., "Statistical analysis of seismic processes on the basis of the diffusion approach", Doklady Earth Sciences, 426:1 (2009), 642-644].
[5] Шевцов Б.М., Сагитова Р. Н., "Диффузионный подход в статистическом анализе сейсмичности Камчатки", Вулканология и сейсмология, 2012, №2, 56-66. [Shevtsov B. M., Sagitova R. N., "A diffusion approach to the statistical analysis of Kamchatka Seis-micity", Journal of Vulcanology and Seismology, 6:2 (2012), 116-125].
[6] Федотов С. А., "О закономерностях распределения сильных землетрясений Камчатки, Курильских островов и северо-восточной Японии", Тр. ИФЗАН СССР, Наука, М., 1968, 121-150. [Fedotov S. A., "O zakonomernostyah raspredeleniya sil'nyh zemletryasenij Kam-chatki, Kuril'skih ostrovov i severo-vostochnoj Yaponii", Tr. IFZAN SSSR, Nauka, M., 1968, 121-150].
[7] Cahoy D.O., Uhaikin V. V., Woyczyski W. A., "Parameter estimation for fractional Poisson processes", Journal of Statistical Planning and Inference, 2010, № 140, 3106-3120.
[8] Laskin N., "Fractional Poisson processes", Communications in Nonlinear Science and Numerical Simulation, 2003, №8, 201-213.
[9] Учайкин В. В., "Автомодельная аномальная диффузия и устойчивые законы", УФН, 2003, № 173, 847-876. [Uchajkin V. V., "Avtomodel'naya anomal'naya diffuziya i usto-jchivye zakony", UFN, 2003, № 173, 847-876].
Список литературы (ГОСТ)
[1] Carbone V., Sorriso-Valvo L., Harabaglia P., Guerra I. Unified scaling law for waiting times between seismic events // Europhys. Lett. 2005. vol. 6. no 71. p. 1036-1042. DOI: 10.1209/epl/i2005-10185-0
[2] Bak P., Christensen K., Danon L., Scanlon T. Unified Scaling Law for Earthquakes // Phys. Rev. Lett. 2002. vol. 88. no 17. p. 178501-1-178501-4.
[3] Shevtsov B., Sheremetyeva O. Fractional models of seismoacoustic and electromagnetic activity // E3S Web of Conferences. 2017. vol. 20, 02013. DOI: 10.1051/e3sconf/20172002013
[4] Шевцов Б.М., Сагитова Р. Н. Статистический анализ сейсмических процессов на основе диффузионного подхода // ДАН. 2009. Т. 426. № 2. C. 254-256.
[5] Шевцов Б.М., Сагитова Р.Н. Диффузионный подход в статистическом анализе сейсмичности Камчатки // Вулканология и сейсмология. 2012. № 2. C. 56-66.
[6] Федотов С.А. О закономерностях распределения сильных землетрясений Камчатки, Курильских островов и северо-восточной Японии // Тр. ИФЗАН СССР. М.: Наука, 1968. С. 121-150.
[7] Cahoy D.O., Uhaikin V.V., Woyczyski W.A. Parameter estimation for fractional Poisson processes // Journal of Statistical Planning and Inference. 2010. no 140. pp. 3106-3120.
[8] Laskin N. Fractional Poisson processes // Communications in Nonlinear Science and Numerical Simulation. 2003. no. 8. pp. 201-213.
[9] Учайкин В.В. Автомодельная аномальная диффузия и устойчивые законы // УФН. 2003. № 173. С. 847-876.
Для цитирования: Шереметьева О. В. Модель процессов релаксации в различных режимах пластических деформаций // Вестник КРАУНЦ. Физ.-мат. науки. 2018. № 5(25). C. 74-82. DOI: 10.18454/2079-6641-2018-25-5-74-82
For citation: Sheremetyeva O. V. Model relaxation processes in the different modes of plastic deformation, Vestnik KRAUNC. Fiz.-mat. nauki. 2018, 25: 5, 74-82. DOI: 10.18454/20796641-2018-25-5-74-82
Поступила в редакцию / Original article submitted: 25.10.2018