Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Бричкин А.В. Направленное разрушение искусственных минеральных сред огнеструйными горелками

.pdf
Скачиваний:
11
Добавлен:
24.10.2023
Размер:
9.3 Mб
Скачать

 

 

 

 

 

 

Таблица 13

Д1, мм

i i < Д1"

At, сек.

m

R

R ^ 1 T

Rep.

 

 

'Id

 

 

 

 

 

5

10,8

10

4,50

0,50

2,00

1,2

0,10

5

10,8

20

5,76

1,76

0,57

1,2

0,16

8

28,8

15

5,12

1,12

0,89

1,2

0,16

10

45,0

30

6,00

2,00

0,50

1.2

0,20

Тепловые константы приняты равными [102, 103]

объемный вес бетона

 

у=2200

кг/м3;

 

теплоемкость

 

с=0,17

ккал/кг,

град;

коэффициент теплопроводности

Я= 1,2

ккал/м.

час.

град;

коэффициент теплоотдачи

с

поверх­

 

 

 

 

ности бетона к воздуху

(при ско­

 

 

 

 

рости ветра 5 м/сек)

 

а = 2 0

ккал/м2

час.

град.

Исходя из скорости резания, через каждые 15 сек. источник тепла продвигался на 24 мм, и температура на этом отрезке скачком поднималась до температуры плав­ ления бетона и поддерживалась в течение следующих 15 сек. Повышение температуры края плиты, вследствие наличия источника нагрева, за время At учитывается по­ вышением электрического напряжения в каждой точке схемы на соответствующую величину.

Результаты расчета представлены в виде изотерм, показывающих распределение температуры в плите в за­ висимости от времени при резании и при остывании (рис.23).

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

а

б

Рис. 23. Распределение температуры в бетонной плите:

a) D момент окончания резки; б) при ее остывании

(через 3 мни. после

окончания резки). Отрезанная часть плиты не удалена.

Направление реза сверху вниз.

 

мы высоких температур, и располагаются дальше в глубь плиты.

Из рис. 23 видно влияние на ход изотерм теплооб­ мена между плитой и окружающей средой. Изотермы на краях плиты загибаются к линии реза.

При остывании плиты температура соседних точек выравнивается. Это хорошо видно по ходу кривых изме­ нения во времени температур для различных точек плиты (рис. 24). Кривая 1 показывает изменение температуры точки на расстоянии 4 мм от края полости реза и 200 мм от его начала. Эта кривая через 7 мин. после начала рез-

м и н

Рис. 24. Изменение температуры точек плиты во времени на рас­ стоянии 200 мм от начала реза, Л=60 мм. Кривые 1—6 на рас­ стоянии от кромки реза, 4, 12, 20, 28, 36, 44 мм соответственно.

Пунктирная — кривая максимальных температур.

ки все более подходит к кривым 2, 3, 4, 5, 6-й, характери­ зующим температуры более удаленных точек (соответ­ ственно 12, 20, 28, 36, 44 мм от линии реза). Пунктирная кривая проведена через максимумы температур этих то­ чек. Область плиты слева от нее характеризуется нагре­ ванием, справа — остыванием. Очевидно, наклон кривой максимальных температур характеризует скорость рас­ пространения тепла в глубь плиты. При построении кри­ вых (рис. 24) использован второй случай охлаждения.

Для определения температурных полей в бетонных конструкциях других толщин (до 300 мм) при огневой резке и последующем остывании также был использован метод конечных разностей, но с применением ЭЦВМ [100]. Применение ЭЦВМ оправдывается, когда ведется боль­ шой арифметический счет и нетрудно составить програм­ му или когда даже при затрате большого труда на раз­ работку алгоритма и составление программы существует

I

і

Перевод числового материала из десятичной системы в двоич­ ную.

I

Засылка счетчика (число таблиц).

Засылка первоначальных данных в рабочую матрицу и выдача их на печать.

|

I

Вычисление температуры во внутренних точках области и вы­ числение Тгр.

Пересылка результатов промежуточной матрицы в первона­ чальную рабочую матрицу.

Пересылка результатов промежуточной матрицы из двоичной системы в десятичную и выдача на печать.

 

Переадресация засылки данных в рабочие ячейки.

И

Сравнение на конец цикла.

д а . !

 

 

нет j

 

останов.

Рис.

25. Блок-схема программы для резки бетонных плит.

114

необходимость в многократном использовании этой программы.

Для решения конечно-разностного уравнения тепло­ проводности (41) с граничными условиями (36) и (47) составлены три программы: для резки бетона и двух слу­ чаев остывания, которые отличаются одна от другой только засылкой первоначальных данных в рабочие ячейки (рабочую матрицу) и печатью таблиц (результа­ тов через интервал времени). Программы для различных толщин бетонных плит аналогичны, отличаются только засылкой температур на кромке реза (определяемых экс­ периментальной кривой, рис.,21), зависящих от скорости резания. Блок-схема программы приведена на рис. 25.

Задача для каждой из толщин на ЭЦВМ «Минск-1» решается за 25—30 мин. В результате расчета для каж­ дой из толщин получены изотермы температурного поля, имеющие вид, подобный изотермам, представленным на рис. 23. Анализ полученных результатов показывает, что с увеличением толщины плиты вид изотерм всех темпера­ тур изменяется, они располагаются дальше в глубь изде­ лия, так как факел горелки дольше находится в полости реза, а изотермы средних температур 200—500° С сильнее загибаются к линии реза — сказывается влияние тепло­ обмена между плитой и окружающей средой (рис. 26).

м м

Рис. 26. Распределение температуры в плите толщиной 300 мм, конец резки.

Изучая температурные поля, можно заметить, что возни­ кающая в начале резки область повышенных температур с течением времени увеличивается и, достигнув опреде­ ленных размеров (изотерма 300°С), прекращает даль­ нейшее продвижение в глубь плиты. Остывание плиты в присутствии отрезанной от нее части и в случае ее удале­ ния выражено изотермами температурного поля различ­ ного вида. В первом случае остывание плиты идет значи­ тельно дольше. Для остывания ее до изотермы макси­ мальной температуры 300° С в первом случае необходимо время в 2 раза большее, чем во втором.

В результате расчета получена подробная картина теплового состояния плиты толщиной от 25 до 300 мм при резке их линейным источником нагрева, а также для двух случаев остывания. Это служит основой для оценки прочности плит вблизи линии реза.

Сравнивая результаты расчета с данными изменения прочности бетона при нагревании от температуры [73— 80], установлено, что при резке бетонных плит различной толщины (от 25 до 300 мм) происходит нарушение проч­ ности максимум на глубину 23 мм от стенки полости ре-

60 .

4*7

Рис. 27. Зависимость зоны нарушения прочности от толщины изделия.

за (изотерма 300°С). На основании полученных данных построена кривая зависимости зоны поражения прочно­ сти от толщины разрезаемой плиты (рис. 27).

4. ТЕМПЕРАТУРНЫЕ ПОЛЯ

ПРИ РЕЗКЕ ЖЕЛЕЗОБЕТОННЫХ ИЗДЕЛИЙ

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

Постановка

задачи.

Пластина

конечного

размера

(0-<х<11 0

< у < : 12) пз армированного бетона (с постоян­

ным шагом

б)

находится

в начальный

момент

времени

при постоянной температуре То. Под воздействием

интен­

сивного теплового источника (факел

огнеструйной

горел­

ки протяженностью d), движущегося

вдоль оси у

с за­

данной скоростью w, плита разрезается. По мере продви­ жения факела на участке d поддерживается температу­ ра, равная температуре плавления бетона Т! : , а затем в каждой из точек температура изменяется по экспери­ ментальной зависимости ср (у, t ) . Закон изменения темпе­ ратуры в каждой из точек после прохождения реза будет одинаков, но смещен по времени. Требуется найти рас­ пределение температуры в теле для всего процесса резки. Очевидно, что полость реза будет осью симметрии, по­ этому решение ищется только для половины плиты. Вследствие низкой теплопроводности бетона и симмет­ ричности расположения арматуры можно ограничиться расчетом для отдельной ячейки, считая, что распределе­ ние температур будет подобным в остальных частях, но с некоторым запаздыванием, определяемым скоростью продвижения факела.

Математически задача сводится к интегрированию следующей системы дифференциальных уравнений (ин­

декс 1 относится

к арматуре, 2 — к бетону; рис. 28):

д 1 - я (

d3 Ti(x,y,t) • d'Trfx.y.t) \ .._

1 Г - М — ш — + — w — h 1 - 1 , 2 {50)

с начальными и граничными

условиями

 

 

 

 

 

T(x,y,t) =

T0 ,

t = 0

 

 

 

х

ЭТ,(хл,0

.=

« ( т п о в - Т о к р ) ,

у = 0

 

 

 

T(x,y,t) =

T(y,t),

х = 0

(51)

 

 

 

 

 

 

- І Ї ! ^ М -

= 0

 

y = b

 

 

 

 

 

 

 

dx

 

 

 

 

J

 

 

 

 

 

 

Здесь

Т п о в , Т о к р соответствен-

 

W'W&ffffli

и 0

т е м п

е Р а т

У Р ы поверхности и

n-t

^ f P ^ f i + i ^

окружающей

среды,

d — коэф­

 

 

 

 

 

 

фициент

 

температуропровод­

 

 

 

 

 

 

ности, X — коэффициент тепло-

БЄТ0Н

шжъпгзоЩ

проводиости,

а коэффициент

теплоотдачи,

х, у— простран­

Рис. 28. Схема

определения

ственные

координаты, t — вре­

мя.

 

 

 

 

 

 

температуры

в точке сопри­

 

Уравнение

(50) в

конечных

косновения

арматуры

с

бе­

 

тоном.

 

 

 

разностях

для случая:

At

 

At

 

1

 

 

 

 

М > 4

(условие

Дхз

 

Дуэ

 

• М

 

 

 

 

 

 

 

 

 

 

 

 

 

устойчивости)

получим

 

 

 

 

 

 

 

Тк, п, ш + 1 —

а , 2

( Т к + 1 .

п , ш

+ Т | с — 1 . n , m

+

Т к ,

m + Т

к , n — 1, ш ) +

 

 

 

 

 

 

 

 

 

М

 

 

 

 

 

 

 

+

( М — 4 а 1 2 ) Т к . п .

 

 

(52)

 

 

 

 

 

 

 

М

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В случае прямоугольной сетки (Дхт^Ау) системы раз­

ностных уравнений

(51)

сводятся,к

виду:

 

Тк.

 

а 12 { М 2

( Т к - 1 . n , m + Т к + 1 ,

п , m ) - f - Щ Т к , n - 1, m - f -

 

 

 

 

 

 

 

 

 

 

 

 

• Т к .

 

ш ))

|_ { М 1 М 3 - 2 а , 2 ( М 1

+ М 3

) }

 

_ДХ2

=

М , > 2 ,

 

 

_ДуЗ =

М 2 > 2 .

 

At

 

 

 

 

 

 

 

At

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения

(52) и (53)

справедливы в пределах

каж­

дой из составных

частей плиты

(полагая

i = l , ai 2

= a i 2 ,

получим соотношения, справедливые для арматуры;

i = 2,

а 1 2 = а 2 2 = 1

— для

бетона). На стыке бетона и арматуры

уравнения

(53) преобразуются к виду;

 

 

Tit, л, ПН-1

(Tk -I, n, m +

Tk+1, n, m )M3 +a1 3 'MiTk. n-1, in +

 

 

 

 

 

 

 

4-MrTk-, п+1. in

+

(Mj.VL, — 2Mo(l ai,,)M,}

Tk, n. m

(54)

 

 

 

 

 

 

 

 

 

 

 

 

 

граничное условие (51)

 

 

 

 

 

 

 

 

Tk+i +

а- Дх

-T o up

 

 

Tk+i

=

X

0,

 

 

 

 

 

 

 

 

 

1

+ _

T ~

 

 

 

где

 

 

 

 

 

 

 

 

При решении задачи распределения температуры для

случая резки

и последующего

остывания

железобетон­

ной плиты

толщиной

60 мм использован статический

электроинтегратор

[101].

 

 

 

 

 

Нахождение температурных полей в железобетонных изделиях толщиной до 300 мм при резке и последующем остывании для двух случаев осуществлялось на ЭЦВМ «Минск-1» по программе, блок-схема которой приведена на рис. 29 с соответствующими поправками для различ­ ных толщин [33].

Теплофизические характеристики материалов, исполь- ' зуемые для расчетов, представлены в таблице 14 [31, 102].

 

 

ч

 

 

грга

 

 

га

Материал

С

час.

 

я

 

а-.

 

 

 

2200

0.27

Арматура (сталь)

7300

0,17

га гр.час.

S

I ,1 25,0

 

Таблица

14

грг

 

О

га

 

ычас.

Я

я

«

о

 

20

0,0018

12

0.02QJ

Ьвоа программы и

получение /

Перевод числового материала

из Ю - 2

2 0 - р яч

Засылка счетчика

[число таблиц) * = 5

Засылка счетчика (черезCKQODKO гзьлиц печага1ыу:4

Заыллка данных в рабочую

матрицу и яоілача и* на печато

Засылка счетчика (зависит

от толщины Бетона г=<(для/) = б)

г=4(дляМ<()

г'//(для/1'30)

^с " r'wifa | И "8 0 столбцов)

 

 

 

Переадресация оегона,

 

 

 

арматуры,

стыка,

 

 

 

 

гран, точек

бетона

 

 

 

 

| Вычисление

 

арматуры

Ьычнсление

внутрен.точек

ічнеление

гран гоуек ветона |

Бетона

и

переадресация

 

 

 

 

дат

 

 

 

 

 

Вычисление стыка

Лереадресэция стыка

и

емугр. точек веюна

 

Рис. 29. Блок-схема программы

для резки

железобетонных

 

 

плит.

 

 

 

 

Результаты

расчета представлены в виде изотерм тем­

пературного

поля (рис. 30, 31). От подобных

изотерм

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

Соседние файлы в папке книги из ГПНТБ