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

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

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

Р— сила, с которой стержень выдергивается;

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.

ПО

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