Определение полей напряжений по ориентациям напряжений

Содержание

Слайд 2

Ориентация осей современного максимального горизонтального сжатия SH,max в Западной Европе и

Ориентация осей современного максимального горизонтального сжатия SH,max в Западной Европе и

континентальной Австралии из мировой базы данных
Слайд 3

Ориентация осей современного максимального горизонтального сжатия SH,max в районе желоба Сунда

Ориентация осей современного максимального горизонтального сжатия SH,max в районе желоба Сунда

у о. Суматра из мировой базы данных
Слайд 4

Ориентация осей современного максимального горизонтального сжатия SH,max в районе конвергенции Аравийской

Ориентация осей современного максимального горизонтального сжатия SH,max в районе конвергенции Аравийской

и Евразийской плит из мировой базы данных
Слайд 5

Моделирование полей напряжений Векторы граничных напряжений ПОСТАНОВКА КЛАССИЧЕСКОЙ КРАЕВОЙ ЗАДАЧИ МСС

Моделирование полей напряжений

Векторы граничных напряжений

ПОСТАНОВКА КЛАССИЧЕСКОЙ КРАЕВОЙ ЗАДАЧИ МСС

Источник массовых сил

?

СИТУАЦИЯ

В ГЕОДИНАМИЧЕСКИХ ПРИЛОЖЕНИЯХ
Слайд 6

Поля напряжений моделируются на основе решения классических краевых задач, поставленных обычно

Поля напряжений моделируются на основе решения классических краевых задач, поставленных обычно

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

Традиционный подход к определению напряжений в геодинамике

Слайд 7

Величины полученных напряжений существенно различаются в разных решениях: Экспериментально определенные генерализованные

Величины полученных напряжений существенно различаются в разных решениях:

Экспериментально определенные генерализованные направления

оси максимального сжатия [Coblentz et al., 1995]

Разница в решениях при традиционном подходе к определению напряжений

Слайд 8

Прямой подход к определению полей напряжений ОСНОВНАЯ ЦЕЛЬ: объединить достоинства измерений

Прямой подход к определению полей напряжений

ОСНОВНАЯ ЦЕЛЬ: объединить достоинства измерений

напряжений in-situ (реальность характеристик получаемого напряженного состояния) и методов математического моделирования (получение сплошных полей полного тензора напряжений)

ОСНОВНЫЕ ИДЕИ ПОДХОДА:
использование экспериментальных данных об ориентации главных напряжений непосредственно в качестве входной информации и отказ от умозрительных гипотез традиционного математического моделирования

Слайд 9

Представление плоского тензора напряжений SH,max SH,min ϕ SH,max и SH,min –

Представление плоского тензора напряжений

SH,max

SH,min

ϕ

SH,max и SH,min – взаимно ортогональные главные значения

тензора напряжений (максимальное и минимальное горизонтальные сжимающие напряжения);
ϕ - угол наклона SH,max к некоторому выделенному направлению в плоскости)

τmax= ½(SH,max-SH,min) – максимальное касательное напряжение;
P= ½(SH,max+SH,min) – среднее напряжение

Определить поле тензора напряжений T(x)– значит найти
3 скалярных поля SH,max, SH,min, ϕ или, что эквивалентно, 3 скалярных поля τmax, P, ϕ

(SH, max=- T1, SH, min=- T2, T1< T2)

Слайд 10

Траектории главных напряжений Траектории главных напряжений (ТГН) – кривые, касательные к

Траектории главных напряжений

Траектории главных напряжений (ТГН) – кривые, касательные к которым

в каждой точке ориентированы по направлению SH,max или SH,min

SH,min

SH,max

Модельный пример построения плоской сетки взаимно ортогональных ТГН по заданным дискретным ориентациям максимального сжимающего напряжения

Слайд 11

2D уравнения равновесия T=T1m1⊗m1+ T2m2⊗m2 В декартовой с.к. x1, x2 с

2D уравнения равновесия

T=T1m1⊗m1+ T2m2⊗m2

В декартовой с.к. x1, x2 с ортами e1,

e2

В ортогональной криволинейной с.к., совпадающей с ТГН

(*)

bi=b∙mi, κi – кривизна i-той ТГН

Слайд 12

2D уравнения равновесия ϕ − угол наклона T1 к оси x;

2D уравнения равновесия

ϕ − угол наклона T1 к оси x; sj

− длина дуги вдоль траектории j-го семейства

− кривизна траектории j-го семейства

Для трех неизвестных функций напряжения T1, T2, ϕ имеется два уравнения равновесия

Если траектории главных напряжений (ТГН) заданы внутри исследуемой области Ω, то функция ϕ известна и уравнения равновесия представляют замкнутую систему уравнений. Величины T1, T2 определяются без предположений об определяющих соотношениях среды.

Слайд 13

2D уравнения равновесия Тензор T(x(s)) задан на линии x=x(s) и в

2D уравнения равновесия

Тензор T(x(s)) задан на линии x=x(s) и в окрестности

линии известно поле ТГН. Задача: определить T(x(s)) в окрестности линии из уравнений равновесия.

При заданном 2D поле ТГН уравнения равновесия являются системой уравнений гиперболического типа с характеристиками, совпадающими с ТГН

Слайд 14

Модельный пример реконструкции ТГН Постановки краевых задач

Модельный пример реконструкции ТГН

Постановки краевых задач

Слайд 15

Аналитические решения для простейших ТГН ТГН, совпадающие с линиями декартовой с.к.

Аналитические решения для простейших ТГН

ТГН, совпадающие с линиями декартовой с.к.

ТГН,

совпадающие с линиями полярной с.к.
Слайд 16

Пример пошаговой процедуры определения поля напряжений по данным сейсмического мониторинга I

Пример пошаговой процедуры определения поля напряжений по данным сейсмического мониторинга

I –

напряжения в блоке перед выемкой грунта

II – определение локальных ориентаций оси SH,max на основе сейсмического мониторинга

III – реконструкция поля ТГН по дискретным ориентациям SH,max

IV – решение уравнений равновесия для определения величин главных напряжений

ABCD – задача Гурса; CDE – задача Коши; CFGE – задача Гурса

CH – разлом, неизвестный до начала работ; CD – выявленная часть разлома

Слайд 17

0A – Срединно-Атлантический хребет (САХ) 0C – зона коллизии Африканской и

0A – Срединно-Атлантический хребет (САХ)

0C – зона коллизии Африканской и Евразийской

плит

tR – сила отталкивания от САХ

tС – сила коллизии

D′DBB′ – Западно-Европейская провинция напряжений (ЗЕПН)

Поле напряжений в литосфере Западно-Европейской платформы

Слайд 18

Поле напряжений в литосфере Западно-Европейской платформы

Поле напряжений в литосфере Западно-Европейской платформы

Слайд 19

Решение накладывает ограничение на напряжение коллизии tC. Если известна величина этого

Решение накладывает ограничение на напряжение коллизии tC. Если известна величина этого

напряжения, то определяется его ориентация, и наоборот. Выбор действительной ориентации tC основывается на сравнении теоретически и экспериментально определяемых режимах напряжений в литосфере ЗЕПН.

Для согласования с наблюдаемыми напряжениями растяжения и сдвиговым режимом напряжений на большей части ЗЕПН вектор tC должен быть близок по направлению к направлению SH,max и вращаться против часовой стрелки при движении на восток.

Поле напряжений в литосфере Западно-Европейской платформы

Слайд 20

Некоторые сведения из ТФКП Δf1=0, Δf2=0 Определить голоморфную функцию f в

Некоторые сведения из ТФКП

Δf1=0, Δf2=0

Определить голоморфную функцию f в круге по

заданной на границе действительной части f1.

Оператор Шварца

Слайд 21

Основные соотношения плоской задачи теории упругости P – вещественная функция среднего

Основные соотношения плоской задачи теории упругости

P – вещественная функция среднего напряжения,

D – комплекснозначная функция девиатора напряжений

Уравнения равновесия

Закон Гука

T11, T12, T22 – компоненты T в декартовых координатах x1, x2

Слайд 22

Основные соотношения плоской задачи теории упругости Уравнения равновесия Закон Гука

Основные соотношения плоской задачи теории упругости

Уравнения равновесия

Закон Гука

Слайд 23

Основные соотношения плоской задачи теории упругости Формулы Колосова-Мусхелишвили Φ, Ψ -

Основные соотношения плоской задачи теории упругости

Формулы Колосова-Мусхелишвили

Φ, Ψ - голоморфные функции,

D – биголоморфная функция

|D| - модуль D

α=arg(D)

α - аргумент D

Искомые функции

Слайд 24

Общий вид вещественной биголоморфной функции (Im c0 =Im c2 =0) Основные

Общий вид вещественной биголоморфной функции

(Im c0 =Im c2 =0)

Основные соотношения плоской задачи теории упругости

-

биголоморфная функция с тем же аргументом (если Dreal>0)
Слайд 25

Типы задач для упругой литосферы Определение поля напряжений по: Общие свойства

Типы задач для упругой литосферы

Определение поля напряжений по:

Общие свойства решений:

1) решение

для поля напряжений зависит от m≥2 произвольных констант

2) решение для ТГН и τmax не зависит от объемных сил b, если они имеют гармонический потенциал

Слайд 26

Особые точки первого рода поля траекторий главных напряжений z=0 – особая

Особые точки первого рода поля траекторий главных напряжений

z=0 – особая

точка первого порядка

T1=T2 при x=x* → D(x*)=0

Существует 3 типа асимптотического поведения ТГН, которые топологически неэквивалентны и устойчивы к возмущениям коэффициентов A, B:

ϕ→ϕ+π

ϕ→ϕ-π

c

b

a

При обходе особой точки угол ϕ касательной к траекториям одного семейства получает приращение +π при |A|>|B| (траектории выпуклы по отношению к особой точке) и приращение -π при |A|<|B| (траектории вогнуты)

Слайд 27

Особые точки поля ТГН Поля ТГН в окрестности особых точек в частных случаях

Особые точки поля ТГН

Поля ТГН в окрестности особых точек в частных

случаях
Слайд 28

Определение напряжений при заданных ТГН в упругости Не каждое поле траекторий

Определение напряжений при заданных ТГН в упругости

Не каждое поле траекторий допустимо

в упругой среде, т.к. при произвольно заданной функции α необходимо двумя функциями P и τmax удовлетворить трем уравнениям упругости: двум уравнениям равновесия и условию гармоничности P.

В частности, гармоническая функция α (Δα=0) допустима в 2D упругости, т.к. в этом случае α является аргументом некоторой голоморфной функции A(z) и, следовательно, аргументом би- голоморфной функции DrealA(z)

Слайд 29

Определение напряжений при заданных ТГН в упругости Определить би-голоморфную функцию D

Определение напряжений при заданных ТГН в упругости

Определить би-голоморфную функцию D по

заданному полю ее аргумента α. Затем функция P вычисляется простым интегрированием уравнений равновесия
Слайд 30

Частное решение Определение напряжений при заданных ТГН в упругости

Частное решение




Определение напряжений при заданных ТГН в упругости

Слайд 31

Частное решение (продолжение) Решая задачу для двух кругов r* и r**

Частное решение (продолжение)

Решая задачу для двух кругов r* и r**

Определение

напряжений при заданных ТГН в упругости

Если в рассчитанной функции D1 зависимость от r* остается, то заданная функция α не допустима в упругости

Слайд 32

Определение напряжений при заданных ТГН в упругости Общее решение

Определение напряжений при заданных ТГН в упругости

Общее решение


Слайд 33

Определение напряжений при заданных ТГН в упругости Общее решение на γ

Определение напряжений при заданных ТГН в упругости

Общее решение

на γ является граничным

значением голоморфной функции


Φ′(е)≠0, радиус R выбирается так, чтобы знаменатель внутри круга не обращался в 0.

Т.к. ImB(z)= 0 на γ, то B(z)=const=c3 >0.

Слайд 34

Определение напряжений при заданных ТГН в упругости Декартовые ТГН Полярные ТГН

Определение напряжений при заданных ТГН в упругости

Декартовые ТГН

Полярные ТГН

Слайд 35

Определение напряжений при заданных ТГН в упругости Параболические ТГН η=const → →

Определение напряжений при заданных ТГН в упругости

Параболические ТГН

η=const