Демонтаж бетона: rezkabetona.su

Главная  Переработка нефти и газа 

Скачать эту книгу

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 [ 72 ] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139

z = z

(5.21)

или эквивалентно [114]

P = r-V

(5.22)

Можно воспользоваться балансом масс (5.9) на фронте волны:

и р ~е + 1

(5.23)

Если объемная деформация на фронте постоянна (например, e{R) = е.), то на фронте (г = К)

р = Я"-Пхг) = Po{J + е,);

(5.24)

X, = R"- / {1 + е.),

и расчет упрощается.

Теперь можно определить функцию f{x) как

2-л л

f{x) = PR-" = Ро{1 + е,){х{1 + е,)у-". (5.25)

Кинематические интегралы (5.17) и (5.25) совместно с пластическим условием (5.6) и балансом импульса на фронте (5.10) позволяют свести динамическое уравнение (5.2) к обыкновенному дифференциальному уравнению, которое может быть проинтегрировано численно.

5.1.3. ЧИСЛЕННЫЙ РАСЧЕТ КАМУФЛЕТНОГО ВЗРЫВА

Рассмотрим вариант построения решения, при котором численный расчет применяется с самого начала. Такой подход позволяет учитывать изменения скорости дилатансии (вплоть



до смены ее знака) одновременно с изменениями коэффициента сухого трения а и сцепления chs = аН. Более того, удается также учитывать упругую составляющую деформаций в пластической зоне II (равно как и существование внешней упругой зоны III).

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

5v дг

дг г J

(5.26)

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

В пластической зоне II упругая часть деформаций исключается предварительно с помощью дилатансионного соотношения (5.5):

дг г J

+ - К

да да

+ V-

= -г=Л

dv дг

(5.27)

В обеих зонах выполнены уравнения баланса (5.1) и (5.2), но в зоне II выполнены упругопластические определяющие соотношения при пределе пластичности (5.6).

Необходимо отметить, что непосредственно перед R (т.е. при /- = R- 0) между зонами II и III выполнено то же самое пластическое условие (5.6).

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

хотя ударный фронт может и не совпадать с этой границей.)

Поскольку вычислительный метод включает в себя численную вязкость, фронтальные скачки заменяются на подвижные зоны высоких градиентов. Некоторые типичные результаты расчета [200] подземного камуфлетного взрыва даны на рис. 5.2-5.4.

Здесь использованы приведенная координата г / s, где S - энергия взрыва. Профили взрывной волны (рис. 5.2) приве-



дены для трех моментов времени (t = 28 t; 56 t; 84 t, где 0 ~ / P скорость продольной волны).

Сплошные линии соответствуют дилатансионному разрых-

У, Р, кеар

СМ/Г


•Z/gt м/кг*-

Рис. 5.2. Профили волн, иалучаемых динамически растущей полостью с учетом дилатансии массива (сплощные линии) и без нее (пунктир)

V,„, м/с

, Т„ КБАР


Г / е"-, м7кгз

Рис. 5.3. Затухание сейсмических волн в окрестности подземного взрыва (сплощные линии соответствуют учету дилатансии, энергия дана в весе эквивалентного заряда ТИТ; 1 м/кг1/=100 м/кт/)




0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 [ 72 ] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139



Яндекс.Метрика