Научная статья на тему 'Математическая модель внутреннего тепло- и массообмена в бетонных плитах, подвергаемых тепловлажностной обработке'

Математическая модель внутреннего тепло- и массообмена в бетонных плитах, подвергаемых тепловлажностной обработке Текст научной статьи по специальности «Строительство и архитектура»

CC BY
152
40
i Надоели баннеры? Вы всегда можете отключить рекламу.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по строительству и архитектуре , автор научной работы — Аксенчик К. В., Шестаков Н. И.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Математическая модель внутреннего тепло- и массообмена в бетонных плитах, подвергаемых тепловлажностной обработке»

что цилиндр свободно расширяется вдоль оси z:

ст„ =-

1 г Iя

[—2 ¡а ■ t(r) -rdr + — Ja • t(r)rdr];

У Л п

1-ц r\J

ст = ———[—7 fa • /(г) • rdr + 1-М- г о

°экв = О/ »

где ст,

(CTi + (a2 -^з)2 + (°з -^i)2

интенсивность напряжения; ст,,ст2,ст3 - главные напряжения. Параметры контакта двух тел (цилиндра и плоскости) определяются по формуле

1 R

+ — ja ■ t(r) ■ rdr - a ■ t(r)~\ ;

¿ = 1,131,

If

1-m2 , i-rê

v 3

2 y

CT. = -

1-Ц

iV

Ja • ¿(r) • rdr - a • /(r)

л

Вычислить интеграл ja-i(r)-r<ir и опре-

0

делить напряжение можно, если известен закон изменения температуры t(r) по толщине цилиндра. Здесь сг,стф,ст2 - главные напряжения: радиальные, тангенциальные и осевые; a - коэффициент термического расширения; Е - модуль упругости; р. - коэффициент Пуассона, соответствующий средней температуре. По теории прочности точки с максимальным эквивалентным напряжением будут самыми опасными, и в них появятся первые признаки разрушения тела [3]:

где Ъ - полуширина полоски контакта; Р - сила сжатия двух тел; Я - радиус цилиндра; 1 -длина цилиндра; Ех и Ег - модули упругости материалов; ^ и ц2 - коэффициенты Пуассона для цилиндра и плоскости. Зная ширину полоски, можно определить угол контакта ср.

Список литературы

1. Блох А.Г., Журавлев Ю.А., Рыжков JI.H. Теплообмен: Справочник. - М.: Энергоатомиздат, 1991. - 432 с.

2. Кривандин В.А., Мастрюков B.C. и др. Металлургическая теплотехника: В 2 т. Т. 1: Теоретические основы: Учебник для вузов. - М.: Металлургия, 1986. - 424 с

3. Писаренко Г.С., Яковлев А.П., Матвеев В.В. Справочник по сопротивлению материалов. - Киев: Наука, 1988.-736 с.

УДК 666.97.035.55

К.В. Аксенчик, Н.И. Шестаков Череповецкий государственный университет

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВНУТРЕННЕГО ТЕПЛО- И МАССООБМЕНА В БЕТОННЫХ ПЛИТАХ, ПОДВЕРГАЕМЫХ ТЕПЛОВЛАЖНОСТНОЙ ОБРАБОТКЕ

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

обмена. В отличие от подобных моделей, в которых применяются аналитические решения [1], в частности для исследования температурных полей ограждающих одно- и многослойных конструкций [3], [7], в данной работе использован численный метод решения дифференциальных уравнений, описывающих процессы тепло- и массообмена в бетонах на шлаковом щебне.

Процессы тепломассообмена в капиллярно-пористых телах, таких как бетон, древесина и т.д., математически описываются системой

дифференциальных уравнений Лыкова [5], которая при отсутствии градиента давлений (УР = 0) имеет вид

^- = KUV2U + K1272T; (1) ох

~ = K2yU + K22V2T, (2)

где II - удельное влагосодержание; Т - температура.

Коэффициенты, входящие в уравнения (1) и (2), равны:

ки=а1=ап,8Т;

К2\ ~ат1Ж'>

с

^22=

С

где ат - коэффициент диффузии влаги; ат -

коэффициент термодиффузии влаги; 5Т - термоградиентный коэффициент; ат1 - коэффициент диффузии пара; гп - удельная теплота парообразования (конденсации); с - приведенная удельная теплоемкость; а - коэффициент температуропроводности; аТтХ - коэффициент

термодиффузии пара.

В пропарочных камерах тепловой обработке подвергают изделия различной формы. В качестве объекта моделирования были выбраны бетонные плиты размером 3200 х 1200 х х 100 мм. На основании того, что толщина плит значительно меньше двух других размеров (рис. 1), плиту можно рассматривать как неограниченную пластину толщиной Н= 25.

Система уравнений (1) и (2) для одномерного случая будет иметь вид:

(3)

(4)

dU дх d2U дх2 + Ки д2Т дх2

дТ дх = К21 д2и дх2 + К22 д2Т дх2

Рис. 1. Схема бетонного образца

Процессы переноса теплоты и влаги, согласно системе уравнений (3) и (4), являются связанными, т.е. температурные градиенты вызывают градиенты влагосодержания и наоборот. Если принять, что эффекты термодиффузии и внутренних фазовых превращений пренебрежимо малы [8], т.е. коэффициенты Ки = К2\ = 0, то система уравнений (3) и (4) распадается на два независимых уравнения -дифференциальное уравнение теплопроводности Фурье для твердого тела и уравнение диффузии, записанное через влагосодержание:

дТ д2Т L

■ = а дх дх

+ ■

ср

dU d2U

-~ат-Г

дх дх

(5)

(6)

где а - коэффициент температуропроводности; Iq - источник тепла, обусловленный процессами гидратации, протекающими в твердеющем бетоне; с - удельная теплоемкость; р - плотность; ат - коэффициент диффузии влаги.

Уравнения (5) и (6) позволяют моделировать сложные процессы тепло- и массообмена без учета их взаимного влияния друг на друга. Данное обстоятельство упрощает решение поставленной задачи.

Учитывая, что плиты направляются на тепловую обработку непосредственно после формования без предварительной выдержки, температура и влагосодержание приняты постоянными по всему объему плиты. Таким образом, получены начальные условия:

Т(х, 0) = Т0 = const при х = 0 и 0 < х < Я; (7) U(x,0) = U0= const прит = 0 и 0<х<Н. (8)

Граничные условия 3-го рода были получены из закона Ньютона - Рихмана:

9 = а[г(х,х)5-Гср(х)],

где ц - количество теплоты, воспринимаемое (отдаваемое) единицей поверхности тела в единицу времени; а - коэффициент конвективной теплоотдачи; Т (х., т) - температура на

поверхности материала; Гср (х) - температура

среды в зависимости от времени. Тогда

при х = 0 и 0<х<хц

а

при х = Н и 0 < х < хц

(9)

(10)

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

W =

T0 + wit-х, если 0<х<хпрог, 7Ш, если тпрог < х <хпрог + хю, тю ~ Щт ■ ec™ xmor + хиз < х < хи

(П)

где 7о - начальная температура в камере; м>\Т -скорость подъема температуры; х - время; м>1Т - скорость охлаждения; хпр0г, хиз, хц - длительность периода прогрева, изотермической выдержки и цикла пропаривания, соответственно. Скорость подъема температуры, скорость охлаждения, длительности периодов регламентируются инструкциями [6] и режимом пропаривания.

Значение д задается условиями внешнего тепло- и массообмена и изменяется по периодам тепловой обработки.

В период прогрева а - коэффициент теплоотдачи при конденсации насыщенного пара или пара из смеси с воздухом; при отсутствии пленки конденсата (в период изотермической выдержки) а - коэффициент теплоотдачи при испарении влаги с поверхности материала; в период охлаждения а - коэффициент теплоотдачи от поверхности материала к холодному воздуху.

Для влагосодержания были заданы симметричные граничные условия 1-го рода:

-при х = 0 и 0<х<хц : t/(0,x)s =ср(х), (12)

-при х = Н и 0<х<хц: U[H,x)s =cp(x), (13)

где ф(х) - функция вида ().

Таким образом, предлагаемая математическая модель внутреннего тепло- и массообмена состоит из двух подмоделей:

- подмодели внутреннего теплообмена, включающей дифференциальное уравнение

(5), начальное условие (7) и граничные условия (9)-(11);

- подмодели внутреннего массообмена, включающей дифференциальное уравнение

(6), начальное условие (8) и граничные условия (12), (13).

Уравнение (6) представляет собой модификацию уравнения (5), называемого в математической физике уравнением диффузии, поэтому метод решения уравнений (5) и (6) с начальными и граничными условиями один и тот же. Математическая модель решалась методом конечных разностей на ЭВМ в программе Mathcad по неявной схеме. Выбор схемы обусловлен ее абсолютной устойчивостью и уменьшением затрат времени счета на ЭВМ.

Для численного эксперимента были взяты следующие входные значения параметров модели: удельный расход цемента та = 324 кг/м3; теплоемкость, плотность, коэффициент теплопроводности бетона на шлаковом щебне с = 830 Дж/(кг • К), р = 1860 кг/м3, X = = 0,58 Вт/(м • К) [4]; начальная температура в плитах То = 20 °С; начальное влагосодержание плит U0 = 0,241 кг/кг; размеры плит 3,2 х 1,2 х х 0,1 м; режим пропаривания: прогрев от 20 до 85 °С - 4 ч, изотермическая выдержка при температуре 85 °С - 6 ч, охлаждение до 40 °С -4 ч; количество теплоты, воспринимаемое (отдаваемое) единицей поверхности тела в единицу времени, по периодам qa = 2500 Вт/м2, qm = 250 Вт/м2, #охл= 250 Вт/м2; среднее значение за весь цикл тепловлажностной обработки коэффициента теплоотдачи а = 508 Вт/(м ■ К)

(7); коэффициент диффузии влаги ат = 1,35 х х Ю"5 м2/с [5]; точность расчетов s = 0,01.

Объемный источник теплоты вследствие гидратации цемента определялся по формуле

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

q = Шц (о, 6t2 +16,5t - 3,5) • 10~3 Вт/м3,

где г - средняя температура периода, °С [2].

Было принято следующее изменение влаго-содержания поверхности материала (граничные условия):

ф(т) =

0,241 +1,465-10~5т при 0 < т < 4 ч,

0,452 при 4<х<10ч,

0,452-0,006т при 10<х<14ч.

Результаты моделирования представлены графически на рис. 2 и 3.

о

100

90

80

70

я 60 а.

£ 50

а.

<L>

s 40

<L>

н 30 20 10 0

ж -+4-+-Л

л*

0 0,02 0,04 0,06 0,08 0,1 Толщина бетонной плиты, м

поверхности и минимальна в центральной плоскости плиты (на глубине 0,05 м от поверхности). В период изотермической выдержки (следующие 6 ч цикла) градиент температур незначителен и составляет не более 5 °С. При этом максимум температур находится в центральной плоскости плиты. В конце периода максимум температур в центральной плоскости охлаждения сохраняется, но градиент температур увеличивается. Распределение влагосодержаний устанавливается к концу периода прогрева и сохраняется приблизительно постоянным в различные периоды тепловой обработки (минимум в центральной плоскости и максимум на поверхности), причем градиенты влагосодержания достаточно велики: АГ7 = = 0,35 кг/кг.

Таким образом, разработана математическая модель внутреннего тепло- и массообме-на в плитах из бетона на шлаковом щебне при тепловой обработке в пропарочных камерах ямного типа. После проверки адекватности созданной модели предполагается провести исследование процессов внутреннего тепло- и массообмена и термонапряженного состояния при различных условиях тепловой обработки: режимах пропаривания, составах бетона, несимметричных граничных условиях.

-0ч -7ч

2,1ч ■ -9,8 ч

-в-3,5 ч

-4-14 ч

Список литературы

Рис. 2. Поле температур в плите из бетона на шлаковом щебне

0,5

0,45I

иI

0,4

* 035

Я

5 0,3 I 0,25

I °'2

I 0,15 « 0,1 0,05 0

О 0,02 0,04 0,06 0,08 0,1 Толщина бетонной плиты, м

» 0ч х 3,5 ч —в—7ч —А—14ч

Рис. 3. Поле влагосодержаний в плите из бетона на шлаковом щебне

Как видно из рис. 2 и 3, в плитах наблюдается параболический характер распределения температур и влагосодержаний. Так, в период прогрева (4 ч) температура максимальна на

1. Гамаюнов Н.И., Испирян P.A., Клингер A.B. Расчет температурных полей в керамзитобетоне при его тепловой обработке // Инженерно-физический журнал. - 1977. -Т. 33, №2.-С. 360-361.

2. Заседателев И.Б., Петров-Денисов ВТ. Тепло- и массоперенос в бетоне специальных промышленных сооружений. — М.: Стройиздат, 1973. - 168 с.

3. Ибрагимов A.M. Нестационарный тепло- и массоперенос в строительных материалах и конструкциях при несимметричных граничных условиях. Ч. 1 // Строительные материалы. - 2006. - № 7. - С. 72-73.

4. Кучеренко A.A. Тепловые установки заводов сборного железобетона. Проектирование и примеры расчета. -Киев: В ища шк., 1977.-280 с.

5. Лыков A.B. Тепломассообмен: Справочник. - М.: Энергия, 1978. - 480 с.

6. Руководство по тепловой обработке бетона и железобетонных изделий. -М.: Стройиздат, 1974.-31 с.

7. Федосов C.B., Ибрагимов A.M., Гущин A.B. Влияние тепловлажностной обработки на прочность железобетонных ограждающих конструкций и изделий // Строительные материалы. - 2006. - № 9. - С. 7 - 8.

8. Федосов C.B., Ибрагимов A.M., Гущин A.B. Применение методов математической физики для моделирования массо- и энергопереноса в технологических процессах строительной индустрии // Строительные материалы. - 2008. - № 4. - С. 65 - 67.

i Надоели баннеры? Вы всегда можете отключить рекламу.