![](/user_photo/_userpic.png)
книги из ГПНТБ / Бричкин А.В. Направленное разрушение искусственных минеральных сред огнеструйными горелками
.pdf
|
|
|
|
|
|
Таблица 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). От подобных |
изотерм |
для бетона они отличаются появлением выступов в местах пересечения арматуры, которые указывают на большую скорость распространения в ней тепла, чем по бетону. В дальнейшем при остывании плиты выступы становятся пологими (что объясняется отдачей тепла от арматуры к бетону), изотермы высоких температур разрываются в местах расположения арматуры, образуя замкнутые кри вые с одной и другой стороны. При длительным остыва-