Файл: Кушнарев Д.М. Использование энергии взрыва в строительстве.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 126

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

жается непосредственно через внутреннюю энергию в идеаль­ ном газе:

Е =

РоѴо

 

 

• 1

 

Если отнести внутреннюю

энергию к единице веса

заряда

Q, т. е. е = — , то

 

 

р 0 = е Р ; ( ѵ * - 1 ) ,

(ХІ.7)

где у*— некоторый эффективный показатель адиабаты (можно считать у*=1,25). Подставляя выражение (XI.7) в (XI.3), по­ лучим:

/ ' \

1/1.25

 

 

1/ = (0,25)1 / 1 , Я 5 К0 ( ^ - j

.

 

(ХІ.8)

Массовую плотность энергии обычно определяют по справоч­

нику. По физическому смыслу она совпадает с энергией

хими­

ческого разложения взрывчатого вещества и по порядку

величи­

ны ее можно считать равной е = 1 ккал/кг

= 4,18-1013

эрг/кг (для

аммонита).

 

 

 

Соотношения (XI.3) — (XI.6) и (XI.8)

в первом

приближении

описывают количественным образом явление расширения по­ лости газовой камеры под действием взрыва. Существенным не­

достатком указанных формул является их незамкнутость,

так

как в них входит

одна неизвестная величина

р* — упругая

кон­

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

физико-

механическим и пластическим свойствам данного грунта.

 

Рассмотрим

вначале

общую

картину,

возникающую

при

взрыве в пластическом грунте, в соответствии с выбранной

нами

моделью грунта.

 

 

 

 

 

 

 

Напишем граничные условия

для упругой области:

 

 

<т„ =

- а і ;

" U . = 0 ,

 

 

(XI.9)

где в ы — тензор напряжений;

и — смещение;

о\—-радиальное

напряжение на границе упругой зоны с пластической.

 

Известно, что в общем

случае, когда деформация

вызывает­

ся не объемными силами, а силами, приложенными к поверх­ ности,

(1 — 2а) Au-fgraddivu = О,

(Х.10)

где а — коэффициент Пуассона.

Введем цилиндрические координаты с осью z по оси цилинд­ ра. При однородных вдоль оси условиях деформация представ­

ляет собой чисто радиальное смещение iir=u(r).

Воспользовав­

шись известной формулой векторного анализа, получим:

grad div г/ = Au

(XI. 11

207


и из уравнения (XI.10) с учетом (XI.11) получим окончательное уравнение равновесия в упругой области:

 

 

graddivu =

0,

(XI. 12

которое надо решать при граничных

условиях

(XI.9).

Находим первый интеграл:

 

 

div и

J_

d(ru) =

const = 2C.

 

 

r

dr

 

 

Интегрируя последнее соотношение, получим:

и — сг-\

.

(XI. 13)

r

Известно, что в цилиндрических координатах r, ср, z компо­ ненты тензора деформаций выражаются через смещение сле­ дующим образом:

- _ L ^ Ф _ I

}h_.

_

д±

дг •

w

г

дц>

г

 

 

02

ф2

г

оф

ôz

 

 

dz

dr

'

(XI. 14)

 

 

 

 

 

 

dur

 

 

 

 

 

 

 

 

 

 

Оф

 

 

 

 

Отличные от нуля компоненты тензора

деформаций

(XI. 13) и

(XI.14) равны:

du

D

 

 

 

 

D

 

 

 

 

 

 

 

 

(XI. 15)

 

 

 

 

 

 

 

 

 

 

Из условия ti[«, = 0 находим

С = 0 и w> r

= —

D

U

= D

 

 

 

 

 

 

 

 

ФФ

 

r*

Постоянная определяется из условия на границе с пластиче­

ской зоной:

 

 

 

 

 

 

 

 

 

 

Предварительно запишем закон Гука в цилиндрических ко­

ординатах:

Е

 

[(l—2o)ulk

 

+

OUetÔik\;

 

 

 

 

 

 

 

( 1 + а ) ( 1 - 2а)

 

 

 

 

Е

 

[(1 -

а) ы

+

сшф ф ];

[

(XI. 16)

(1

+ а ) ( 1 - 2о)

 

 

 

 

 

 

оФФ =

Е

 

К 1 ст) и ф Ф

+ °Чг1.

 

 

+ а ) ( 1 -- 2а)

 

 

(1

 

 

 

 

 

 

 

 

где Е — модуль Юнга.

 

 

 

 

 

 

 

 

 

Затем находим D:

 

 

 

 

 

 

 

 

 

 

« , r | r = * = - a 1 -

R2(l

 

+ о )

'

 

 

 

 

 

 

 

 

 

208


ß = Pl С + 0 ) ^2

 

E

 

Поэтому тензоры деформаций и напряжений

описываются

формулами:

 

^ = _ M l _ h £ L ( A j 2 ;

(VI. 17

ФФß

М тг Р

" « • " " • ( - г ) •

< Х , Л 8 )

 

 

Зная тензор деформаций, легко вычислить перемещение гра­ ниц упругой зоны:

ô = \ arrdr = - g*

0 ) R* j ' - ^ - = 0 1

G ) Д. (XI.19

Продольные напряжения в данном случае, очевидно, равны пулю:

а = (а 4- er ) а = О

Во второй зоне главные напряжения связаны условием пла­ стичности:

 

f a

— а

а

=

а

,

(XI.20)

 

<т2 _|_ гт22

п

п

г, 2

 

 

 

где От —

' m m

 

(сю

т ''

V

'

параметр пластичности.

 

 

 

 

 

 

Для определения главных напряжений требуется решить

уравнение равновесия —— о . а

= 0.

 

 

 

 

 

 

Чтобы записать последнее уравнение в цилиндрических пе­ ременных, следует воспользоваться формулами перехода от де­ картовых переменных к полярным:

 

А- = г cos ф;

у = г sin ф; z = г.

Однако

практически бывает удобно непосредственно выве­

сти условия

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

основе статических соображений

в интересующей нас криволинейной ортогональной системе ко­ ординат.

Тензор напряжений и тензор деформаций для плоской зада­ чи в полярной системе координат могут быть записаны в виде:

,<т„ 0

0 a

ФФ

На рис. 87 представлены система напряжений для элементар­ ного объема, выделенного в полярных координатах, и картина его смещения.

14—50

209



В соответствии с принятыми обозначениями уравнение рав­ новесия, записанное в виде суммы проекций всех усилий на направление радиуса, запишется в виде:

a,rcd + arr

- I -

ab — 2a add(p

0.

 

dr

ФФ

 

Производя замену cd=rdQ; ab = (r-j-dr)dip и отбрасывая беско­ нечно малые высшего порядка, получим уравнение

° " - ° w

+ r - ^ - = 0. (XI.21)

 

dr

Исключив из уравнения рав­ новесия (XI.21) и условия плас­ тичности (XI .20) нормальное азимутальное напряжение от , получаем дифференциальное уравнение

 

 

 

 

 

 

 

 

 

dOrr = = _

_рѴг

 

Рис. 87. Система

напряжений

для

 

 

 

dr

2

'

элементарного

объема

и картина

 

 

 

 

 

 

его

смещения

 

 

 

- Г

От

1 — 0,75 ' °-

 

 

 

 

 

 

 

 

 

 

1

От )

Для решения данного уравнения воспользуемся методом раз­

деления переменных. Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

=

Г - ^ - f - C ;

2стт

= X;

 

 

 

 

 

 

dx

 

 

Г С.

 

 

 

 

 

 

V 1 — З*2

x

 

 

 

 

 

 

 

 

 

Вычислим входящий сюда интеграл:

 

 

 

V Зх

=

Sin СС;

У 1

Зх-

cos а-

,

cos a

j

 

dx = ——

da;

 

 

 

 

 

 

 

 

 

 

Уз

 

 

l/3

=

t g 6 0 ° ;

 

 

 

 

 

 

 

 

J

 

 

dx

 

 

 

 

cos

a da

 

(XI.22)

 

 

 

 

 

 

 

 

sin a

 

 

 

Vi

Зх2

—.

1/3

cos

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

J1

cos a da

 

 

 

 

 

 

 

 

2

sin (60° —

a ) '

 

 

 

210