![](/user_photo/_userpic.png)
книги из ГПНТБ / Бричкин А.В. Направленное разрушение искусственных минеральных сред огнеструйными горелками
.pdfР— сила, с которой стержень выдергивается;
U— периметр арматуры;
1— длина заделки стержня в бетоне. Первоначальное увеличение сил сцепления гладкой
арматуры с бетоном при 100—250° С объясняется увели чением его прочности при таких температурах. Сцепле нию арматуры с бетоном способствует также увеличение сил трения за счет расширения арматуры при нагрева нии. При дальнейшем возрастании температуры, вслед ствие увеличения разности температурных деформаций
бетона и |
арматуры, трение уменьшается. Около стержня |
в бетоне |
возникают большие местные напряжения и об |
разуются |
остаточные деформации, из-за которых арма |
тура |
при охлаждении образца может отстать от |
бетона, |
|||||
и величина |
сцепления сильно снизится. Связь |
арматуры |
|||||
с бетоном |
нарушается. При 275—325° С сцепление арма |
||||||
туры с бетоном близко |
к величине их сцепления |
при |
|||||
обычных температурах, |
при 450° С —• снижается |
на |
60— |
||||
70%. |
При |
нагревании |
усадочные |
напряжения |
и |
напря |
|
жения от разности температурного |
расширения |
металла |
и бетона создают сложное напряженное состояние бето на. Возникают температурные контактные напряжения вследствие разницы (более чем в 2 раза) коэффициентов линейного расширения бетона и стали и в 30—35 раз коэффициента теплопроводности. В цементном камне и заполнителях происходит деформация, связанная с изме нением тепловлажностного состояния бетона. Одновре менно в нем происходят пластические деформации, при водящие к перераспределению напряжения между бето ном и арматурой и к увеличению напряжения на контакте между ними.
Сцепление бетона с арматурой периодического про филя значительно выше, Чем с гладкой. При 450° С сцеп ление арматуры периодического профиля с бетоном на портландцементе несколько ниже, чем в охлажденном состоянии.
Определение зоны нарушения прочности при огневой резке может быть выполнено двумя способами. Первый — сравнение результатов расчета температурного поля бе тонной или железобетонной плиты с прочностью бетона при нагревании. Второй — различные эксперименталь ные методы [109—127]. Ниже приводится описание этих способов и полученных результатов.
2. ТЕМПЕРАТУРНЫЕ ПОЛЯ
ПРИ ОГНЕВОЙ РЕЗКЕ БЕТОННЫХ ИЗДЕЛИЙ
Исходные положения. При огневой резке плиты стен ки образующейся полости на участках длиной 60—65 мм вдоль линии реза находятся в течение 35—40 сек. (ско рость резания около 6 мічас) при температуре плавления бетона. Остальная часть плиты находится в состоянии теплообмена с окружающей средой, температура кото рой постоянна.
При движении,огпеструйного прибора по линии реза влияние источника теплоты на начальные участки реза постепенно уменьшается. Установлено, что уже на рас
стоянии 250—300 мм от начала реза |
струя прекращает |
свое действие на' начальный участок, а тепло по нормали |
|
к направлению реза распространяется |
не более, чем на |
10 см. Вследствие этого при расчетах принята плита раз мером 60X300X300 мм. По .толщине ее температура принимается одинаковой и равной средней линии.
Так как распространение тепла происходит симмет рично в обе стороны от линии реза, достаточно рассмат ривать температурное поле только одной половины пли ты. Расчет проведен в два этапа. На первом —• определя ется распространение тепла в плите поперек канавки в процессе резки. После прекращения действия факела на чинается остывание плиты. При таком решении задачи учитывается теплообмен с окружающей средой.
Рассматриваются два случая выравнивания темпера туры в плите. Первый, когда отрезанная часть плиты не удаляется и влияет па ход остывания последующего участка плиты. Второй случай, когда отрезанная часть плиты удаляется, остывание краев происходит интенсив ней из-за лучшего теплообмена с окружающей средой.
Дифференциальные уравнения теплопроводности со ответственно для первого и второго этапов можно выра зить:
(33)
г д е . а = — |
—коэффициент |
температуропроводности, |
|
|
м21час; |
|
|
f(x s t) |
—член, характеризующий источник |
тепла, |
|
|
зависит от координаты и времени. |
|
|
Начальные условия |
|
|
|
|
при t = 0; |
Т = Т0 . |
(35) |
Более сложно установление граничных условий, ха рактеризующих теплообмен и температурный режим на внешнем контуре исследуемого тела [81].
Условия первого рода характеризуются тем, что на поверхности тела задана температура, могущая зависеть от координат и времени.
Trp. = T ( x , y , t ) . |
(36) |
При условиях второго рода на границе задается вели чина теплового потока, т. е. температура на границе мо жет иметь любую величину при заданном градиенте:
- g - ( 0 , t ) = Kt)- |
(37) |
Граничные условия третьего рода определяют тепло обмен на границе с окружающей средой:
X -|T -(0,t) = a[Tr p .(Q,t)-Tc p .(t)]. |
(38) |
Для случая резки на границе реза и при первом слу чае остывания приняты граничные условия первого рода. Для других краев плиты учитывается теплообмен с окру жающей средой, поэтому принимаются граничные усло вия третьего рода. Учет теплообмена верхней и нижней поверхностей делает задачу трехмерной и значительно усложняет ее решение. Как показала предварительная оценка, величиной поверхностного теплообмена можно пренебречь ввиду его малости (2—3%). Это объясняется тем, что поверхность плиты, особенно нагретого участка, способного отдать теплоту, не велика по сравнению с тол щиной. Время выравнивания температуры также незна чительно, так как рассматривается охлаждение только до 300° С. Пренебрегая поверхностным теплообменом, мы создаем для плиты худшие условия (остывание происхо-
днт дольше действительного), что не противоречит полу ченным ниже выводам.
Для задания граничных условий первого рода экспе риментально определено изменение температуры во вре мени в каждой точке плиты по линии перемещения факе
ла. Измерение проводилось с помощью термопар |
[82], |
||||
которые |
вводились в |
расплав вслед за резаком |
на |
рас |
|
стоянии |
50 мм одна |
от другой. Температура |
измерялась |
||
через каждые 10 сек. в начальный момент |
и через |
30— |
|||
60 сек. при плавном |
ее падении. Применялись |
платино- |
родий-платииовые термопары, выдерживающие длитель ное нагревание до 1300° С и кратковременное до 1600° С. Величина возникающей термоэлектродвижущей силы регистрировалась пирометрическими милливольтмет рами.
Градуировка вольтметров производилась при темпе ратуре свободногоконца термопары 0°С. В опытах температура свободного конца поддерживалась равной
20° С, что потребовало |
введения поправок |
[82], вычислен |
||||||
ных по формуле: |
|
|
|
|
|
|
|
|
|
Т„ст. = |
Ту к . + |
к ( Т 1 |
- Т 0 ) > |
|
|
(39) |
|
где Т„ст. |
—истинная |
температура; |
|
|
|
|||
Т'уК. |
— температура, |
показываемая |
прибором; |
|
||||
Т0 |
—• температура |
свободного конца, |
при |
кото |
||||
|
рой производилась градуировка; |
|
|
|||||
Ті |
— действительная |
температура |
свободного |
|||||
|
конца; |
|
|
|
|
|
|
|
1< |
— коэффициент, |
зависящий от |
типа |
термопа |
||||
|
ры и интервала измеряемой температуры. |
|||||||
Значения «к» взяты из таблиц [82]. Результаты изме |
||||||||
рений приведены в виде кривой (рис. 21) |
изменения тем |
|||||||
пературы |
во времени |
от момента |
плавления до полного |
|||||
остывания |
в данной точке. Опыты |
показали, |
что |
бетон |
марки 300 (па щебне из гранита) плавится при темпера туре 1300° С.
На рис. 21 даны значения температуры, измеренные тремя термопарами, находящимися на расстоянии 50, 100, 150 мм от начала реза. Как видно, разброс точек не значителен, кривая построена по средним значениям. Та ким образом, граничные условия вдоль линии реза за даны в виде: Тг р . =ср(х, I) .
5 |
10 |
15 |
20 |
25 |
ЗО |
Рис. 21. Изменение температуры оплавленного слоя бетона.
Исследованию и методам решения уравнения тепло проводности посвящено большое число работ [81, 83—86].
Аналитическое решение задач теплопроводности со сложными краевыми условиями и изменением термиче ских коэффициентов, часто встречающихся на практике, связано сбольшими техническими трудностями.
В настоящее время получили распространение раз личные методы упрощенного решения задач теплопровод ности. Среди численных методов интегрирования диффе ренциальных уравнений теплопроводности часто приме няется метод конечных разностей [87—91]. Этот метод не накладывает ограничений на температуру среды, на коэффициенты X, а, а и удобен для решения различных технических задач.
3. МЕТОД КОНЕЧНЫХ РАЗНОСТЕЙ
Сущность метода конечных разностей (метода сеток) заключается в том, что дифференциальное уравнение, начальные и граничные условия заменяются системой ко-
нечно-разностных |
уравнений |
[89—91]. При таком пред |
||
ставлении время |
также изменяется |
дискретно |
с ша |
|
гом At. |
|
|
|
|
Метод дает достаточно |
точные |
решения дифферен |
||
циальных уравнений даже |
при сравнительно |
больших |
интервалах между узлами сетки. Область задания функ ции разбивают в простейшем случае на равные участки, хотя это не обязательно [92].
Конечно-разностная аппроксимация уравнения (34) выразится в виде
Тк. п. ш |
— Тк.іі.пі |
/ Т к , п—1 пі - 2 Т к . п. in-!-Тк. п -!- 1. m і |
||||||
|
|
дГ |
- а |
\ |
Ш |
|
|
+ |
і |
Тк—1, n, i n — 2 Тк, п, m ~\~ Тк+1. п. т\ |
|
/лг\\ |
|||||
+ |
|
|
ду! |
|
)• |
|
(40) |
|
Если в (40) принять Дх = Ау = Д1и "^Г" = ^> т |
0 после пре |
|||||||
образования получим |
|
|
|
|
|
|||
гр |
Тк. п + 1 , т + Т к , п — l . m + T k + 1 , iLm+Tk — 1, n.m + (*) — 4)Тк,п га |
|||||||
ik, n,m+l — |
|
|
|
|
~ |
|
|
(41) |
|
|
|
|
|
|
|
|
|
Это уравнение дает возможность рассчитывать |
темпера |
|||||||
турное поле шаг за шагом, исходя из начального |
распре |
|||||||
деления температуры. |
|
|
|
|
|
|||
При выполнении условия [90, 93] |
|
|
|
|||||
|
|
|
|
|
|
|
|
(42) |
погрешность |
применения |
уравнения (41) будет |
не ниже |
|||||
порядка малости относительно Д1, и решение |
разностного |
|||||||
уравнения |
будет сходиться |
к решению |
дифференциаль |
|||||
ного уравнения |
(34), т. е. |
|
|
|
|
|||
|
|
|
lim ( Т * п |
t - |
Тк , „, t ) = О, |
|
|
|
|
|
|
Д1-»-0 |
|
|
|
|
|
где T*k,n,t |
— истинное значение температуры в точке п, к. |
|||||||
Выражение |
(42) связывает шаг At по времени с ша |
|||||||
гом по координате. |
|
|
|
|
|
|||
Для частного случая при т]=4 уравнение |
(41) перепи |
|||||||
шется в виде |
|
|
|
|
|
|
||
rp . |
Tk,n+1, ш И- Тк, п—І.пі + |
T k + l , n,m + |
Tic— 1,п, т |
, , „ \ |
||||
1п,к.т + 1 == |
|
|
|
|
т |
|
. |
(43) |
Температура элемента п, к в каждый следующий мо мент времени (m'-f-.l)At равна среднеарифметическому значению температуры смежных элементов.
Недостаток метода конечных разностей — значитель ная вычислительная работа, необходимая для получения значений температурного поля в различные моменты времени.
Применение различных вычислительных устройств значительно облегчает процесс счета.
Большое распространение получили различного рода
моделирующие устройства, |
действие которых |
основано |
|
на общности дифференциальных |
уравнений, |
описыва |
|
ющих различные процессы. |
Такая |
замена |
возможна |
вследствие аналогии между |
тепловыми, гидравлически |
ми, электрическими и другими явлениями.
Электрическое моделирование [94—96] обладает ря дом преимуществ, обусловивших его широкое примене ние: надежность, стабильность, высокая точность, прису щая электроизмерительным приборам, возможности из менения параметров модели и проч.
Ниже дается принцип действия электростатических интеграторов, созданных в лаборатории кафедры общей физики Казахского государственного университета им. С. М. Кирова [97—99]. На одном из них и была ре шена задача.
В основе действия статических электроинтеграторов (СЭИ) положено воспроизведение конечно-разностной (пространственной и временной) аппроксимации уравне ний в частных производных. Дискретность решения во времени определяет ряд преимуществ статических инте граторов: возможность сравнительно легкого учета источ ников и стоков, коэффициентов, зависящих от функции или координаты, времени и т. д., так как можно активно вмешиваться в процесс решения на любой его стадии.
Интегратор состоит из низкоомного потенциометра, коммутационной панели, счетно-решающего элемента и стрелочного нуль-гальванометра. Обмотка потенциометра выполнена нихромовым проводом диаметром 1,0 мм. От обмотки сделано 1000 отводов на коммутационную па нель.
Таким образом, деление приложенного к концам по тенциометра напряжения (около 1—2 вольт) обеспечи вается с точностью до третьего знака после запятой.- -
В качестве счетно-решающего элемента используется пятплучевая звезда из магазинов сопротивлений, на кон цы которых подаются соответствующие напряжения с коммутационной панели. Результат конечно-разностной аппроксимации уравнения теплопроводности находится путем сравнения напряжений в узле счетно-решающего элемента и на коммутационной панели с помощью щупа
инуль-прибора и заносится в таблицу. Принципиальная схема интегратора представлена на'
рис. 22.
а)
Рис. 22.
а) принципиальная схема статического электроинтегратора КП —
коммутационная панель, H U — нуль-гальванометр; |
Rj |
—сопротивления |
|
решающего элемента; С) |
решающий элемент |
для |
осуществления |
граничного |
условия третьего рода. |
|
Соответственно заданным сопротивлениям |
Rk - i,n , |
Rk, п - 1 , Rlt, n, Rk+l ,n, Rk.n+l С ПОМОЩЬЮ ШТЄККЄрОВ 11Э |
коммутационной панели набираются величины потенциа
лов V к—1,п, m i |
Vk, П - 1 . П 1 , |
Vk, n, i n , Vk, ц+1, п і , Vk+l,n, m |
имитирующих |
значения |
температур в соответствую |
щих точках к + 1 , и; к, n—1; к, п; к—1, п; к, п + 1 в мо
мент |
времени - mAt. |
Шестым |
поисковым |
штек- |
кером с помощью нуль-гальванометра |
на коммутацион |
|||
ной панели отыскивается |
значение потенциала Vk, п. m+i, |
|||
дающее |
температуру в |
точке k, п в момент |
времени |
( m + l ) A t .
На основании закона Кирхгофа можно записать:
Vk.n, m + l - Vk-l,n, m |
, |
Vk, n,m + 1 — Vk, n-l.m |
, |
|
||
Rk - l,n |
"t" |
Rk,n-1 |
"t~ |
|
||
і Vk, n.m + 1 — Vk, n, m |
, |
Vk, n. m + l — Vk, n+1, m |
+ |
|
||
Ri |
T |
|
Rk,n + i |
|
||
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
Vk, n. m + l |
— Vk + I. n, |
0 |
|
(44) |
|
|
d |
|
= |
|
||
|
Rk+1, 11 |
|
v |
' |
||
|
|
|
|
|
|
Если все сопротивления счетно-решающего элемента
ВЫбраТЬ ОДИНаКОВЫМИ Rk-I, n = Rk,n-l =Rk,n + l = Rk+l,n = R, то после соответствующих преобразований получим:
, 7 |
_ |
Rt(Vk, n+l,m + Vk, n-1. ni+Vk+1, п, ni+Vk-l,ii,ni)+Vk, n. m |
||
V k, n, m + l |
|
4 R T + R |
|
|
|
|
|
|
(45) |
Сравнивая |
(45) с (41), |
находим условие моделирования: |
||
|
|
R t = l ; |
R = 7 , - 4 . |
(46) |
|
Задавая значения At и А1 при известном а, |
получим ве |
||
личину її, а, следовательно, и значения R и R t . |
||||
|
При работе на интеграторе для удобства |
обычно со |
||
противления R принимают равными единице |
и устанав |
|||
ливают эту величину на четырех магазинах |
сопротивле |
|||
ний, а на пятом магазине задают R t —-ур |
|
Условия (46) устанавливают правила моделирования
для внутренних точек рассматриваемой области, причем предполагается, что'на границе функция Т задана.
Рассмотрим моделирование граничных условий. Гра ничные условия третьего рода (38) в конечных разностях
запишутся в виде:
-ду-Т^ + аТср. |
|
XTj + аДІТср. |
|
|
Т г Р - = |
і X |
= |
Х + <Ш |
( 4 7 ) |
а |
+ ~кг |
|
|
|
На интеграторе условие (47) выполняется следующим образом. При подходе к границе тела три сопротивления
из пяти, составляющих счетно-решающий элемент, от ключают. В результате получают ячейку из двух сопро тивлений RCp и R (рис. 22).
По закону Кирхгофа для такой схемы справедливо равенство:
Vrp.— Vcp. i
RcP . 1
Vrp.— Vj
Ri
Отсюда
|
V,RC p.+R,VcP . |
|
|
= |
R1 + ReP. |
• |
<4 8 > |
І Із уравнения (48) и (47) находим условие моделиро вания на границе бетон — среда
RcP. = X; |
R, = o.Al. |
(49) |
Следовательно, выбрав два магазина сопротивлений, устанавливают на них RCp. =Я и Ri = a - A l и штеккерами набирают значения потенциалов V c p . и V], численно рав ные (или пропорциональные) температурам Тс р . и Ті и определяют Vrp. поисковым штеккером и нуль-гальвано метром. Найденное значение потенциала пропорциональ но или равно температуре поверхности в момент време ни (m-f-l)At. Продолжая таким образом вычисления, получают распределение температуры в рассматриваемой области Для любого момента времени:
m - At(m = 0, 1, 2, 3....).
При решении задачи контур плиты был разбит.сеткой с шагом Д1==8 мм. На выбор размера ячеек влияют два противоположных фактора. Для более точного и подроб ного решения требуется уменьшать шаг сетки, а для со кращения времени решения необходимо укрупнение сет ки и увеличение временного интервала. При этом размер шагов сетки связан между собой условием устойчиво сти решения.
Значения величин, входящих в условие моделирова ния при выбранных А1 и At, сведены в табл. 13.
ПО