МИНИСТЕРСТВО ТРАНСПОРТНОГО СТРОИТЕЛЬСТВА
ГОСУДАРСТВЕННЫЙ ВСЕСОЮЗНЫЙ ДОРОЖНЫЙ
НАУЧНО ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ
(СОЮЗДОРНИИ)
МЕТОДИЧЕСКИЕ РЕКОМЕНДАЦИИ
ПО РАСЧЕТУ ТЕМПЕРАТУРНЫХ ПОЛЕЙ,
НАПРЯЖЕНИЙ И ДЕФОРМАЦИЙ В
ЦЕМЕНТОБЕТОННЫХ ПОКРЫТИЯХ
Москва 1976
Одобрены Минтрансстроем
Изложены основные закономерности теплопередачи в дорожных цементобетонных покрытиях с учетом радиационного баланса и конвективного теплообмена, приведены формулы и примеры расчета по ним температурных полей, напряжений и деформаций в покрытиях.
Для практического использования подробно изложена методика применения гармонического анализа, позволяющего выразить аналитически любую функцию распределения температуры по толщине покрытия во времени. Учитывая пожелания проектных и строительных организаций, все формулы даны в доступной для практического применения форме.
Изменения суточного хода температуры воздуха, интенсивности солнечной радиации, конвективного и кондуктивного теплообмена вызывают на поверхности и по толщине дорожных покрытий температурные напряжения, величина которых в ряде случаев может достигать прочности бетона. Поэтому является обязательным учет как температурных напряжений и деформаций при расчете трещиностойкости покрытий и работы швов, так и температуры твердения бетона при разработке технологии строительства покрытий.
В «Методических рекомендациях по расчету температурных полей, напряжений и деформаций в цементобетонных покрытиях» обобщены отдельные результаты экспериментальных и аналитических исследований, проведенных в Союздорнии (В.А. Черниговым, Н.И. Мусориным, О.Б. Федотовой, В.А. Лапшиным, Е.И. Броницким, В.А. Зельмановичем, Б.Б. Самойленко, Г.С. Бабаяном), в других институтах (Л.И. Горецким, И.А. Медниковым, Б.С. Раевым-Богословским, Е.А. Палатниковым, А.В. Павловым, Э.Д. Бондаревой, В.Д. Садовым, В.А. Воейковым) и за рубежом. Настоящие «Методические рекомендации» предназначены для использования проектными и строительными организациями на стадии как проектирования, так и строительства бетонных покрытий.
«Методические рекомендации» разработал канд. техн. наук В.А. Чернигов.
Замечания и предложения просьба направлять по адресу: 143900 Балашиха-6 Московской обл., Союздорнии.
1. Температурный фактор определенным образом влияет на формирование физических свойств бетона и на напряженно-деформированное состояние дорожной одежды в период эксплуатации дороги. При известных функциях изменения температуры поверхности и по толщине цементобетонного покрытия можно рассчитать:
величину и повторяемость температурных напряжений в плитах разной длины;
продольную устойчивость покрытия;
температурные деформации плит и швов;
глубину промерзания и оттаивания грунта под дорожной одеждой.
2. Самостоятельное значение имеет теплотехническое обоснование технологических процессов, обеспечивающих:
трещиностойкость покрытия до нарезки поперечных швов в затвердевшем бетоне;
заданную скорость набора прочности бетоном;
требуемую температуру твердения бетона по условиям стойкости поверхностного слоя покрытия и формирования его оптимальной структуры;
выбор требуемого способа ухода за свежеуложенным бетоном;
бетонирование покрытий в зимних условиях;
бетонирование покрытий в условиях сухого и жаркого климата.
Основные закономерности теплопередачи в покрытии
3. Изменение температуры в поверхностном слое и по толщине покрытия происходит вследствие меняющегося в нем теплосодержания, обусловленного процессами радиационной, конвективной и кондуктивной теплопередач. Направление тепловых потоков, составляющих тепловой баланс на поверхности покрытия, схематически показано на рис. 1.
Рис. 1. Схема составляющих теплового баланса на поверхности покрытия днем (а) и ночью (б):
Qс - коротковолновая суммарная радиация (прямая и рассеянная); Qв - конвективная теплопередача; Qк - кондуктивная теплопередача; Jп - длинноволновое излучение покрытия; Jа - длинноволновое излучение атмосферы; Jэф = Jп - Ja - эффективное излучение покрытия
4. Радиационная теплопередача R выражается радиационным балансом и является мерой притока лучистой энергии к поверхности покрытия:
R = (Qп + Qр)(1 - A) - Jэф, (1)
где Qп, Qр - потоки тепла соответственно от прямой и рассеянной радиации;
A - величина, характеризующая отражение тепла в зависимости от цвета поверхности, называемая коэффициентом альбедо поверхности.
5. Конвективная теплопередача Qв между поверхностью покрытия и воздухом выражается законом Ньютона:
Qв = αк (Tп - Tв), (2)
где Tп, Tв - температура соответственно поверхности покрытия и воздуха*);
αк - коэффициент конвективной теплопередачи, зависящий от скорости ветра и перепада температуры ΔT = Tп - Tв.
*) Здесь и далее температура принята в °С.
6. Кондуктивная теплопередача от поверхности к подошве или от подошвы к поверхности покрытия определяется (в данном случае) приближенно законом Фурье:
где T(y,t) - функция распределения температуры по толщине покрытия во времени;
λ - коэффициент теплопроводности бетона.
7. С учетом направления тепловых потоков R, Qв и Qк (см. рис. 1) уравнения тепловых балансов на поверхности покрытия имеют вид:
днем: R = Qв + Qк или Qв = R - Qк; (4)
ночью: Qк = Qв + Jэф или Qв = Qк - Jэф. (5)
В уравнениях (4) и (5) не учтены затраты тепла на испарение влаги с поверхности покрытия, так как при сухой погоде эти затраты ничтожны.
8. Подставив в уравнения (4) и (5) значения Qв = α(Tп - Tв), найдем температуру поверхности покрытия Tп и коэффициент конвективной теплопередачи α:
днем: и (6)
ночью: и (7)
Формулы (6) и (7) дают ясное представление о влиянии каждой составляющей теплового баланса на температуру поверхности покрытия и будут далее использованы в расчетах.
9. Численные значения параметров радиационного баланса приведены в справочниках по климату или могут быть получены на метеорологических станциях, расположенных в районе строительства дороги. В СНиП II-А.6-72 приведены данные по радиации за июль и суточный ход освещенности горизонтальной поверхности по месяцам.
Известно, что суточный ход температуры воздуха является периодическим и приближается к простому гармоническому, например синусоидальному, виду с периодом 24 ч. Гармонической функцией, но более сложного вида, является также суточный ход радиационного баланса.
В зависимости от месяца длительность и величина прямой и рассеянной солнечной радиации в дневные часы суток различна. По этой причине происходит периодическое изменение температуры поверхности покрытия в течение любых суток года, которое всегда можно строго описать тригонометрическим рядом.
10. Тепловые потоки в эксплуатируемом покрытии практически всегда направлены перпендикулярно к поверхности покрытия. И только до устройства обочин, вблизи оголенных торцов покрытия, наблюдается двухмерное температурное поле.
11. Закономерности изменения температуры по толщине покрытия дороги можно описать качественно следующими известными тремя законами Фурье о распространении температурных волн в полупространстве:
амплитуда колебаний экспоненциально убывает с глубиной;
температурные колебания в покрытии происходят со сдвигом фазы, т.е. максимумы (минимумы) температуры, например, на подошве запаздывают по сравнению с максимумами (минимумами) температуры на поверхности;
глубина проникания температурной волны возрастает с увеличением периода колебаний температуры на поверхности покрытия.
12. При толщине покрытия 20 см в натурных условиях наблюдается отставание максимальной температуры на подошве покрытия в среднем на 4 - 5 ч по сравнению с максимальной температурой на поверхности.
13. При ясном небе (облачность 0 - 2 балла) максимальная температура поверхности покрытия наблюдается к 14 - 15 часам, а минимальная - в 4 - 6 часов, что соответствует наибольшей и наименьшей температуре воздуха в эти часы. В это время изменение температуры по толщине покрытия приближается к линейному виду.
14. В ясные зимние ночи температура поверхности покрытия может быть ниже температуры воздуха, что следует из формулы (7), при Jэф > Qк. Это явление наиболее вероятно для Средней Азии, Казахстана и Восточной Сибири.
15. Летом температура поверхности покрытия ночью, как правило, выше температуры воздуха на 2 - 6 °С и зависит от скорости ветра и облачности.
16. Через 2-3 ч после восхода солнца возникает почти безградиентное распределение температуры по толщине покрытия (12 - 24 см).
17. В период захода солнца температура поверхности близка к среднесуточной температуре поверхности и выравнивания температуры по толщине покрытия не наблюдается. С этого времени температура поверхности становится ниже температуры подошвы, и через 1 - 2 ч после захода солнца тепловые потоки полностью направлены к поверхности покрытия.
18. При установившихся суточных периодических колебаниях температуры поверхности покрытия колебания ее суточных амплитуд по толщине практически происходят относительно оси среднемесячных распределений температуры (в расчетном месяце) с периодом, равным году (рис. 2, а). Максимальная и минимальная среднемесячные температуры поверхности наблюдаются в конце июля и января.
Рис. 2. Схема колебаний суточных амплитуд A температуры по глубине y, относительно оси среднемесячной температуры в январе (Я), июле (И) и в первый месяц от начала промерзания относительно оси 1
19. В тех случаях, когда после ясных дней устанавливается пасмурная погода или происходит похолодание, среднемесячное распределение температуры по толщине до 2 - 3 м сильно искажается и не является осью, относительно которой происходят суточные колебания амплитуд температуры. В этих случаях при известном распределении температуры на поверхности покрытия для решения задачи распределения температуры и направления тепловых потоков по толщине покрытия может быть использован гармонический анализ.
20. В первый месяц от начала промерзания или оттаивания грунта основания и земляного полотна (примерно до 40 - 50 см) суточные амплитуды температуры колеблются относительно средней температуры по толщине покрытия, выражаемой линейным уравнением (рис. 2, б):
(8)
где Tc - среднесуточная температура на поверхности покрытия;
hб - толщина покрытия;
hпр - глубина промерзания-оттаивания грунта под подошвой покрытия;
y - ордината с отсчетом от поверхности покрытия.
21. Изложенные закономерности теплопередачи в бетонных покрытиях являются общими для любых климатических условий, влияние которых проявляется через величины и повторяемость суточных амплитуд температуры на поверхности покрытия в зависимости от суточных амплитуд температуры воздуха, интенсивности и длительности в течение дня солнечной радиации и числа ясных, полуясных и пасмурных дней в году.
22. Решение теплотехнических задач возможно при известных граничных и начальных условиях. Граничное условие задается функцией изменения температуры на поверхности покрытия и на определенной глубине, а начальное условие - функцией распределения температуры по толщине в начальное время.
23. Максимальная и минимальная температура поверхности покрытия в течение суток может стать в ряде задач граничным условием или оказаться вспомогательной предельной величиной температуры при задании граничного условия, например, гармонической функцией. При отсутствии кондуктивной теплопередачи в покрытии (после 2 - 3 ч от восхода солнца, см. п. 16) температура на поверхности покрытия может быть принята за начальное условие в задачах, например, распределения температуры по толщине.
24. Максимальную температуру поверхности покрытия Tпmax днем, соответствующую максимальной температуре воздуха Tв, найдем по формуле (6):
Тепловой поток Qк от поверхности покрытия толщиной h приближается в это время к стационарному и согласно формуле (3) выражается как
Qк = 5λ · град T1h, (10)
где град T1 - максимальный градиент температуры в покрытии, равный 60 °С/м в умеренном и 85 °С/м в континентальном климате;
λ - коэффициент теплопроводности бетона.
25. Минимальную температуру поверхности покрытия Tпmin ночью, соответствующую минимальной температуре воздуха Tв, найдем по формуле (7):
26. Тепловой поток Qк к поверхности покрытия в формуле (11) приближается в это время к стационарному и равен
Qк = 5 · λ · град T2 · h, (12)
где град T2 - максимальный градиент температуры в покрытии, равный 40 °С/м в умеренном и 57 °С/м в континентальном климате.
27. При безградиентном температурном поле в покрытии, наблюдаемом после 2 - 3 ч от восхода солнца, температуру поверхности Tп и по толщине T(y) определим по формуле (6), приняв Qк = 0:
28. Для практических расчетов суточное изменение температуры T(0,t) поверхности покрытия с периодом T0, равным, например, 24 ч, выразим простой гармонической функцией:
где Tc = Tп(max) - A0;
В данном случае начальная фаза выражает расчетное отклонение температуры в долях периода T0 от положения равновесия - средней температуры Tc.
29. Изменение суточной температуры поверхности T(0,t) покрытия с учетом колебаний температуры в годовом цикле (см. п. 18) найдем как сумму суточных и годовых колебаний:
T(0,t) = T1 + A1sinω1t1 + A0sin(ωt - φ0); (15)
где T1 - среднегодовая температура поверхности покрытия;
A1 - среднегодовая амплитуда колебаний температуры на поверхности покрытия;
Tп(max,июль) - максимальная температура поверхности покрытия в июле, определяемая по формуле (9);
Tп(min,январь) - минимальная температура поверхности покрытия в январе, определяемая по формуле (11);
A0(июль), A0(январь) - максимальные суточные амплитуды колебаний температуры на поверхности покрытия соответственно в июле и январе, определяемые по формуле (14);
T12 - период колебаний годовой температуры, равный 12 месяцам;
t1 - число месяцев, назначаемое от 1 до 12.
30. Расчет по формулам (9) - (15) приведен в приложении 1, а описание температуры поверхности тригонометрическим рядом - в приложении 2.
31. При стационарном кондуктивном теплообмене в покрытии распределение температуры T(y,t) по толщине выражается линейными функциями:
при T1 = Tп - Tпд и Tп > Tпд; (16)
при T2 = Tпд - Tп и Tп < Tпд; (17)
где Tпд - температура на подошве покрытия.
При этом в срединной плоскости плиты температура, заданная функциями (16), (17) и (18 - 19), соответственно равна:
32. Изменение температуры T(y,t) в любые часы суток по толщине покрытия при граничном условии на поверхности в виде гармонической функции может быть приближенно найдено на основе известного решения о распространении температурных волн в полупространстве:
0 ≤ y < ∞, Ty = 0 = A0cosωt,
где a - коэффициент температуропроводности бетона, м2/ч.
В формуле (20) экспонента дает коэффициент снижения амплитуды Aa по толщине покрытия, а - изменение амплитуды во времени со сдвигом фазы на величину (см. п. 11).
33. Используя формулу (20) и начальное условие о нулевом градиенте в покрытии после 2 - 3 ч от восхода солнца, получим расчетную формулу:
где φ0 - сдвиг фазы температуры поверхности покрытия от Tпmax до Tнач при t = 0,
выражаемый в градусах или радианах (например, при Tпmax в 14 часов, Tнач в 7 часов и To = 24 ч будем иметь или φ0 = 90 + 15 = 105° с отсчетом времени t = 0 в 7 часов).
Величину коэффициента b0 в формуле (21) найдем из равенства T(0,0) = T(y,0) или
откуда (22)
При наличии экспериментальных данных о суточном ходе температуры поверхности покрытия или расчетных данных о температуре поверхности через 2 - 3 ч после восхода, в 13 - 14 часов и во время захода солнца расчет распределения температур по толщине можно произвести приближенно на основе гармонического анализа путем подбора постоянных коэффициентов в ряде Фурье (см. приложение 2).
Расчетные формулы термонапряженного состояния покрытия
34. Температурные напряжения всегда возникают в покрытии при несостоявшихся температурных деформациях - продольных и коробления (изгиба). Поэтому, температурные напряжения в покрытии нельзя измерять датчиками деформаций; их можно определить расчетом, например на основе известной формулы С.П. Тимошенко, которая для цилиндрического изгиба плиты (края покрытия) имеет вид:
где σт - температурное напряжение в покрытии;
α - коэффициент линейного температурного расширения бетона;
E - модуль упругости бетона;
μ - коэффициент Пуассона, для бетона равный 0,15 - 0,20, что дает (1 - μ2) ≈ 1;
T(y) - расчетная функция распределения температуры по толщине покрытия.
В правой части формулы (23) первый член выражает напряжения σ1т, возникающие в покрытии при невозможности деформаций продольных и коробления; второй член - σ2т при невозможности только продольных деформаций; третий - σ3т при невозможности только деформаций коробления. Следовательно, для свободного края покрытия без поперечных швов получим
σ1т = EαT(y). (24)
Для свободного края середины плиты, которая под действием собственного веса не может коробиться и где беспрепятственно возникает продольная температурная деформация, получим
35. Если бы плита была невесомой и могла полностью коробиться и продольно смещаться, то при нелинейном распределении температуры по толщине в ней возникли бы напряжения, которые принято называть собственными напряжениями, определяемые по формуле (23). В этом случае напряжения σ1,2т найденные по формуле (25), означают сумму напряжений от невозможности коробления (изгиба) плиты и собственных напряжений, т.е., иначе, разность первых двух членов (σ1т - σ2т) не равна третьему члену σ3т в формуле (23) или
σ1,2т = σ1т - σ2т = σт + σ3т (26)
Если возможно частичное коробление плиты, то напряжения уменьшаются пропорционально коэффициенту Cx, определяемому по известному графику Бредбери. В этом случае формула (26) примет вид:
σ1,2т = σт + σ3т · Cx = σ1т - σ2т - σ3т + σ3т · Cx = σ1т - σ2т - σ3т(1 - Cx). (27)
При Cx = 1, что означает невозможность коробления плиты, получим формулу (26) или соответствующую ей (25). Вывод формулы (27) справедлив при наличии напряжений растяжения при изгибе на подошве плиты, обусловленных третьим членом в формуле (23).
36. С понижением или повышением температуры плита соответственно сокращается или удлиняется. При этом на ее подошве возникают силы трения-сцепления, в результате которых в плите появляются осевые напряжения σр растяжения или сжатия, определяемые на основе известной формулы Кулона N = hQ, или
N = σр · h · b; (28)
откуда
где h - коэффициент трения-сцепления между подошвой плиты и основанием;
γ - плотность бетона;
l - длина плиты между поперечными швами.
Более сложные формулы, по сравнению с выражением (28), учитывающие деформации плиты и функцию температуры, представляют интерес в исследовательских целях. Для инженерных расчетов вполне пригодна формула (28), так как напряжения при изменении коэффициента трения от 0,5 по 1,5 сравнительно малы.
При современных конструкциях дорожных одежд увеличение длины плиты на каждые 10 м приводит к незначительному увеличению напряжений: на 1 - 2 кгс/см2. Поэтому в коротких плитах длиной 4 - 8 м напряжения σр в большинстве случаев можно не учитывать.
Расчет температурных напряжений при линейном распределении температуры по толщине покрытия
37. При линейном распределении температуры в покрытии и невозможности продольных деформаций и коробления температурные напряжения по толщине определяют по формуле (24) с подстановкой соответствующей функции T(y) из выражений (16) - (19). Из формулы (24) следует, что взаимовлияние между напряжениями в различных плоскостях покрытия не возникает при повсюду одинаковой толщине покрытия и равных величинах E и α бетона по толщине и длине покрытия.
38. Если плита не имеет возможности коробиться (Cx = 1), а продольные деформации проявляются полностью, то расчет температурных напряжений σт производят по формуле (25) с учетом T(y) из (16) - (17), или
Из формул (29) - (30) видно, что по толщине плиты возникают только напряжения изгиба. Напряжения со знаком «плюс» соответствуют напряжениям сжатия при изгибе, со знаком «минус» - растяжения при изгибе. Продольные деформации и соответствующие им напряжения, выражаемые вторым членом правой части формулы (25), пропорциональны перепаду температуры в срединной плоскости плиты.
39. При частичном короблении плиты выпуклостью вверх (Cx < 1) температурные напряжения находят по формуле (27), что равносильно умножению правой части формулы (29) на коэффициент Cx < 1:
при (31)
Расчет температурных напряжений при гармоническом изменении температуры в покрытии
40. При отсутствии в покрытии продольных деформаций и коробления температурные напряжения определяют по формуле (24) с учетом гармонической функции температуры T(y,t), выраженной, например, формулой (21). Чтобы определить напряжения σ1т на подошве или поверхности покрытия по формуле (24), достаточно в формуле (21) принять соответственно y = h или y = 0.
41. При свободных продольных деформациях расчет температурных напряжений в плитах, не имеющих возможности коробиться (Cx = 1) или имеющих частичную возможность коробления (Cx < 1), производят по формуле (27), которую запишем в полном виде:
где T(y,t) - гармоническая функция температуры по толщине покрытия, заданная, например, формулой (21).
При Cx = 1 формула (32) соответствует формуле (25).
Для интегрирования функции (21) в формуле (32) следует выразить cos(α - β) = cosαcosβ + sinαsinβ и далее произвести интегрирование согласно известным табличным интегралам:
Конечная расчетная формула σ1,2т получается громоздкой и мало приемлемой для практики из-за вероятных ошибок при вычислениях. Поэтому целесообразное в каждом конкретном случае построить по формуле (21) график распределения температуры в плите, найти значения интегралов по правилу трапеций (приложение 3) или Симпсона и произвести вычисления σт по формулам (25) или (32).
Расчет температурных деформаций плит и поперечных швов в покрытии
42. В общем случае продольные температурные деформации Δlт или относительные температурные деформации εт плит равны:
Δlт = αl(Tн - Tк); (33)
(34)
где l - длина плиты;
Tн - начальная температура, соответствующая расчетному началу охлаждения или нагрева плиты;
Tк - конечная температура, соответствующая расчетному концу охлаждения или нагрева плиты.
43. Для определения расчетных температур Tн и Tк воспользуемся вторым членом правой части формулы (23):
(35)
Напряжения σ2т возникают в плите (см. п. 34) при невозможности только продольных деформаций. Поэтому средняя температура в покрытии, которую можно принять за расчетную начальную Tн или конечную Tк, будет равна
При линейном распределении температуры по толщине плиты расчетная температура Tc, как это следует из формулы (36), при будет равна температуре в срединной плоскости, или (см. п. 31).
При нелинейном распределении температуры по толщине покрытия средняя температура может значительно отличаться от температуры в срединной плоскости, что объясняется изменением волны температуры в плите по экспоненциальной зависимости и отставанием во времени максимальной амплитуды (см. п. 12). Находить среднюю температуру следует по формуле (36).
44. Продольную температурную деформацию Δlт плиты за любое время суток Δt = t1 - t2 определяем по формуле (33), вычислив предварительно Tн и Tк по формуле (36):
где T(y,t1), T(y,t2) - суточные функции температуры по толщине покрытия, выражаемые формулой (21), в моменты времени t1 и t2.
Интегрирование функции (21) в формуле (37) следует производить приближенно по правилу трапеций.
Продольные деформации Δlт плиты за месяц или год с достаточной точностью можно определить, найдя среднюю температуру в данный момент времени, равную температуре и срединной плоскости плиты при линейном перепаде по толщине.
Для этого вычисляют температуру поверхности плиты по формуле (15) и далее по формулам (16) или (17) находят расчетную температуру в срединной плоскости плиты, приняв градиенты температуры град T в зависимости от климата по пп. 24 - 26.
45. Ширина раскрытия поперечных швов между примыкающими плитами одинаковой длины, казалось бы, должна быть равна продольной деформации одной плиты и определяться по формулам (33) - (37). Однако в натурных условиях геометрическая и физическая середины плит не совпадают, что объясняется неравномерным распределением сил трения-сцепления по подошве плиты и сил от возможного смерзания подошвы или боковых граней плит с грунтом. Натурные измерений деформаций 10 - 12 последовательно расположенных швов дают обычно деформации отдельных швов, различающиеся в 1,3 - 1,5 раза. Сумма деформаций всех 10 - 12 швов обычно равна расчетной величине, а распределение деформаций соответствует закону Гаусса. Учитывая это обстоятельство, следует полученную по формуле (33) ширину раскрытия паза поперечного шва принимать за среднюю величину при среднеквадратическом отклонении σ = 0,15Δlт. Следовательно, максимальная и минимальная расчетные деформации паза шва будут равны:
Δl = αl(Tн - Tк)(1 ± 0,45).
Наряду с температурными деформациями в первые 1,5 - 2 года наблюдается усадка бетона плит, которая эквивалентна деформациям от понижения температуры на 10 - 12 °С. При резинобитумных мастиках, которыми периодически (через 2 - 3 года) заполняют швы, едва ли целесообразно учитывать в расчетах деформации паза от усадок бетона, так как усадка бетона приводит к увеличению ширины паза. Если применять герметики работающие в пазе более 4 лет, то при расчете деформации паза следует увеличивать расчетный перепад температуры (Tн - Tк) на 10 °С. Учет усадки бетона приведет к необходимости нарезки пазов шва большей ширины (приложение 4).
В период эксплуатации цементобетонных покрытий на уклонах дорог раскрытие швов между плитами меньше, чем на горизонтальных участках. Это явление обусловлено некоторым сползанием плит по направлению уклона покрытия. Поэтому на таких участках дороги для оценки предельной растяжимости мастик-герметиков нет надобности производить расчет раскрытия швов в летне-осенний период до начала промерзания основания.
46. У всех видов мастик-герметиков, применяемых в настоящее время для заполнения пазов швов, с понижением температуры резко снижаются вязкие и деформативные свойства. В зависимости от климатических условий расчетная минимальная температура, соответствующая предельной величине раскрытия шва, может изменяться от минус 5 °С в южных районах СССР до минус 50 °С и ниже в северных районах. Сама же мастика, находясь в пазе, будет иметь минимальную температуру зимой, практически равную температуре поверхности покрытия. Поэтому изложенные расчеты температурных полей и деформаций швов указывают на необходимость дифференцированного назначения ширины паза шва не только в зависимости от длины плит предельной продольной деформации, но и в зависимости от предельной растяжимости мастик в расчетном диапазоне температур заданного климатического района. Это означает, что при одинаковых свойствах мастики и геометрических размерах плит ширина паза шва в южных районах должна быть меньше, чем в северных.
Пример 1. Найти температуру поверхности покрытия в Московской обл. в 14 часов июля.
Дано: Tb = 25 °С; облачность 0 - 2 балла; Qп + Qр = 559 + 105 = 664 ; Jэф = 72 ; A = 0,3; град T = 0,6 °С/см = 60 °С/м; h = 20 см; α = 20 ; λ = 2 .
По формуле (10) найдем Qк = 5 · 2 · 60 · 0,2 = 120 ккал.
По формуле (9) получим искомую температуру поверхности
Пример 2. Найти температуру поверхности покрытия в Московской обл. в 7 часов июля.
Дано: Tb = 14 °С; облачность 0 - 2 балла; Qп + Qр = 270; Jэф = 72; град T = 0; Qк = 0; A = 0,3; αк = 20.
По формуле (13) найдем °С.
Пример 3. Найти минимальную температуру поверхности покрытия в Московской обл. в 4 часа июля.
Дано: Tb = 10 °С; Jэф = 70; λ= 2; Qп + Qр = 0; град T = 40 °С/м; h = 20 см; αк = 3.
По формуле (12) имеем Qк = 5 · 2 · 40 · 0,2 = 80 ккал.
По формуле (11) получим
Пример 4. Найти максимальную температуру поверхности покрытия в Московской обл. в 14 часов июля с учетом колебаний температуры в годовом цикле.
Дано: Tп(max,июль) = 38,6; Tп(min,июль) = 13,3, Tп(min,январь) = -30 °С.
Определим по формулам (14) и (15) величины А0, A1, T1:
A0 = (38,6 - 13,3) : 2 = 12,65 °С;
При t1 = 3 месяцам, φ0 = 0 и t = 6 ч по формуле (15) найдем
Расчет T(0,t) по формуле (15) дает лучшее приближение к результатам натурных измерений T(0,t) в дневные часы и худшее - в утренние, что объясняется отклонением в утренние часы хода температуры от синусоиды.
Пример 5. Найти температуру подошвы покрытия в Московской обл. в 14 и 19 часов июля.
Дано: Tср = 26 °С (см. пример 1 и 2); A0 = 12,65 °С; a = 0,003 ; h = 20 см; φ0 = 15 °C.
Найдем параметр b0 по формуле (22):
Температуру подошвы покрытия в 14 и 19 часов определим по формулe (21):
Максимальная температура подошвы покрытия 29,5 °C (так как cos0 = 1) наблюдается в 19 часов, а максимальная температура поверхности - в 14 часов.
Суточный ход температуры поверхности покрытия отличается от периодической функции, выражаемой, например, простой синусоидой. Если сложить отдельные синусоиды с частотами и периодами, кратными наименьшей из них, то получится периодическая функция более сложного вида:
Разложение каждого члена ряда (1) по формуле для синуса суммы приводит к ряду Фурье:
Для первых пяти гармоник коэффициенты a0, a1, a2, a3, a4, a5, a6, b1, b2, b3, b4, b5 находят по следующим формулам:
Вычисление величин u и v, входящих в формулы (3), производят, пользуясь табл. 1 для двенадцати ординат:
Таблица 1
Пример 1. В результате натурных замеров через каждые два часа в течение суток получены величины температуры поверхности покрытия, указанные в табл. 2 (крайние значения температуры на расстоянии одного периода должны находиться в одинаковых фазах).
Таблица 2
Время суток, ч |
7 |
9 |
11 |
13 |
15 |
17 |
19 |
21 |
23 |
1 |
3 |
5 |
7 |
Температура поверхности, °С |
17 |
23 |
29 |
33 |
34 |
32 |
27 |
22 |
20 |
17 |
15 |
14 |
17 |
Ординаты y |
y0 |
y1 |
y2 |
y3 |
y4 |
y5 |
y6 |
y7 |
y8 |
y9 |
y10 |
y11 |
y0 |
Требуется выразить аналитически суточный ход температуры T(0,t), т.е. определить коэффициенты a0, an, bn и число членов ряда (2), дающих в сумме достаточное приближение к измеренной температуре.
По табл. 1 найдем суммы u и разности v:
|
17 |
23 |
29 |
33 |
34 |
32 |
27 |
14 |
18 |
17 |
20 |
22 |
|
Суммы (u) |
44 |
37 |
44 |
50 |
54 |
54 |
Разности (v) |
-10 |
9 |
14 |
16 |
14 |
10 |
По формуле (3) определим коэффициенты:
Аналогично получены a2 = -1,6; a3 = 0; a4 = 0; a5 = 0; a6 = 0,08; b1 = 8,22; b2 = -0,146; b3 = 0,5; b4 = -0,146; b5 = 0,27.
Вычисления an и bn проверим по формулам:
y0 = a0 + a1 + a2 + a3 + a4 + a5 + a6 = 23,5 - 4,93 - 1,6 + 0 + 0 + 0 + 0,08 = 17,08;
Вычисления коэффициентов правильны, и ряд (2) примет вид:
При T0 = 24 ч и t = 0 получим T(0,0) = y0 = 17 °С в 7 часов. Если T0 = 24 ч, t = 6 ч и , то T(0,6) = 33,01 °С в 13 часов. Таким образом, имеем практически равенство расчетной температуры с измеренной. Причем, от первых двух членов ряда получаем T(0,6) = 32,88 °С.
Чтобы выразить T(0,t) в виде ряда (1), необходимо найти амплитуды An и начальные фазы φn каждой гармоники по формулам
и (4)
В зависимости от знаков an и bn, находят четверть, в которой лежит угол φn, по табл. 3.
Таблица 3
an |
bn |
φn |
φn = f(αn) |
+ |
+ |
I |
φn = αn |
+ |
- |
II |
φn = 180° - αn |
- |
- |
III |
φn = 180° + αn |
- |
+ |
IV |
φn = 360° - αn |
Тогда будем иметь угол φ1 в IV четверти, A0 = a0 и далее:
По числу -0,6 найдем угол, равный 31°, и обозначим его α1 = 31°. Согласно табл. 3 φ1 = 360 - 31 = 329°. Аналогично находят A2, A3, ..., An и φ2, φ3, ..., φn.
Таким образом, получим:
или
При известной функции температуры T(0,t) поверхности покрытия, выраженной рядом Фурье, расчет температуры по толщине покрытия можно приближенно произвести на основе решения о распространения температурных волн в полупространстве. Применительно к бетонным покрытиям такое решение использовал в 1938 г. К. Эберле в виде:
Формула (5) равнозначна формуле вида:
(6)
Используя третий закон Фурье и данные натурных измерений (через каждые 2 ч) температуры в покрытии за сутки, К. Эберле решал две задачи расчета:
- коэффициента теплопроводности бетона по известной зависимости отношений амплитуд:
где
- глубины проникновения суточных амплитуд по формуле (5).
Эберле пришел к выводу, что по формулам (5) и (7) нельзя найти коэффициент теплопроводности бетона. Вычисленный таким образом коэффициент λ почти в 2 раза больше измеренного. Данные расчета затухания по глубине амплитуд температуры хорошо согласуются с натурными измерениями, и на глубине 203 см от поверхности покрытия суточные амплитуды температуры практически равны нулю.
В 1974 г. болгарский инженер Б. Волчев опубликовал следующую формулу:
в которой задан нулевой градиент температуры при t = 0. Способа определения параметра z0 Б. Волчев не приводит, но обращает внимание на то, что результаты расчета T(y,t) по этой формуле почти совпадают с данными натурных измерений.
Пример 1. Найти температурные напряжения на поверхности и подошве плиты при невозможности коробления и линейном распределении температуры по толщине.
Дано: T1 = 38,6 - 27 = 11,6 °С; E = 350000 кгс/см2; α = 10-51/°С; h = 20 см; Cx = 1.
По формуле (29) найдем напряжения:
на поверхности покрытия
на подошве
На подошве возникнут напряжения растяжения при изгибе, а на поверхности - напряжения сжатия при изгибе.
Пример 2. Найти температурные напряжения в 14 часов на подошве плиты при невозможности коробления и волновом распределении температуры по толщине.
Дано: E = 350000 кгс/см2; α = 10-51/°С; h = 20 см; Cx = 1, T(у,t) по формуле (21); Tср = 26; A0 = 12,6; a = 0,003.
Согласно указаниям п. 41 вначале найдем распределение температуры по толщине плиты. В примере 5 приложения 1 были найдены T(20,14) = 27, T(0,14) = 38,6. По формуле (21) найдем температуру на глубине 5, 10 и 15 см от поверхности. На глубине 5 см будем иметь:
Аналогично найдем T(10,14) = 31,2 °С и T(15,14) = 28,5 °C.
Разобьем толщину плиты на 8 равных частей Δh = 2,5 см и по интерполяции найдем температуру в плите через 2,5 см. Значения температуры приведены в табл. 1.
Таблица 1
T(y,14) |
38,6 |
36,5 |
34,6 |
33 |
31,2 |
29,9 |
28,5 |
27,7 |
27 |
h, см |
0 |
2,5 |
5 |
7,5 |
10 |
12,5 |
15 |
17,5 |
20 |
y |
y0 |
y1 |
y2 |
y3 |
y4 |
y5 |
y6 |
y7 |
y8 |
По следующему правилу трапеций найдем значение интеграла в формуле (25):
Температурные напряжения на подошве плиты найдем по формуле (29):
Из примера 2 следует, что
при гармонической функции температуры (21) напряжения на подошве плиты почти на 4 кгс/см2 меньше, чем при линейном перепаде температуры (см. пример 1);
средняя температура покрытия, равная 31,7 = 635,5 : 20 и вызывающая продольные деформации плиты, больше температуры в срединной плоскости;
изменение температуры по толщине покрытия в различные часы суток вызывает соответствующие изменения температурных напряжений. Максимальные градиенты температуры в покрытии, приближающиеся к линейным, сохраняются в плите короткое время, которое в практических расчетах можно принять равным не более 1 ч.
Пример 3. Найти температурные напряжения на подошве плиты при частичной возможности коробления, волновом распределении температуры по толщине и расчетных параметрах примера 2.
Расчет температурных напряжений произведем по формуле (32) с коэффициентом меры коробления Cx = 0,6. Сумма первых двух членов правой части, в формуле (32) была найдена в примере 2 и равна -16,45 кгс/см2. Для нахождения интеграла третьего члена воспользуемся правилом трапеций. С этой целью разобьем толщину плиты на восемь равных частей Δh = 2,5 см, найдем ординаты каждой части подынтегральной функции и запишем значения ординат в табл. 2.
Таблица 2
y0 |
y1 |
y2 |
y3 |
y4 |
y5 |
y6 |
y7 |
y8 |
|
h, см |
0 |
2,5 |
5 |
7,5 |
10 |
12,5 |
15 |
17,5 |
20 |
|
386 |
273 |
173 |
83 |
0 |
-75 |
-142 |
-208 |
-270 |
Примечание. Значения функции T(y,14) взяты из табл. 1 примера 2. Интеграл третьего члена определим по правилу трапеций:
По формуле (32) получим:
Сравнение примеров 2 и 3 показало, что укорочение плиты (Cx = 0,6) уменьшает напряжения растяжения при изгибе на подошве плиты до 8 кгс/см2 вместо 16,45 кгс/см2 при Cx = 1.
Пример 1. Найти продольную температурную деформацию плиты при суточном изменении температуры в период от 14 часов до 4 часов следующих суток в июле.
Дано: T(y,t) выражена формулой (21); Tср = 26 °С; A0 = 12,6 °С; a = 0,003 ; h = 20 см; φ0 = 0; α = 10-51/°С; l = 6 м.
Расчетные параметры Tн и Tк определим по формулам (37).
В примере 2 приложения 3 была найдена Tн = 31,2 °С. Для определения Tк найдем распределение температуры по толщине (через 2,5 см) по формуле (21). Тогда температура поверхности в 4 часа при t = 14 ч равна:
Аналогично найдем температуру по толщине плиты, результаты расчета которой приведены в таблице.
y |
y0 |
y1 |
y2 |
y3 |
y4 |
y5 |
y6 |
y7 |
y8 |
h, см |
0 |
2,5 |
5 |
7,5 |
10 |
12,5 |
15 |
17,5 |
20 |
Т(y,4) |
15,2 |
16,2 |
17,2 |
18,3 |
19,6 |
20,7 |
21,9 |
22,6 |
23,3 |
По правилу трапеций вычислим интеграл и Tк:
Теперь найдем деформацию плиты по формуле (33):
Δlт = 10-5 · 600(31,2 - 19,5) ≈ 0,07 см.
Пример 2. Найти ширину паза шва сжатия в покрытии, расположенном в двух разных климатических районах. Для заполнения паза принята мастика одинакового состава.
Дано: α = 10-5 1/°С; l = 5 м. В первом районе Tн = 20 °С, Tк = -10 °С и предельная относительная деформация растяжения мастики равна ε1 = 0,3. Во втором районе Tн = 20 °С, Tк = -30 °С и ε2 = 0,20.
Очевидно, предельные относительные деформации ε1 и ε2, соответствующие условиям разрыва мастики или отрыва мастики от стенки паза, должны быть разные. Причем ε1 > ε2, так как с понижением температуры увеличивается модуль упругости мастики. В этом примере предполагается, что предельные относительные деформации мастики установлены в результате испытаний с учетом ее выносливости, старения и поперечного сечения мастики в пазе шва.
Относительные деформации мастики в шве εм и шва εт равны
где b - ширина шва;
к - коэффициент увеличения температурной деформации шва вследствие несовпадения физической и геометрической середин плит, наибольшее значение которого равно 1,45 (см. п. 45).
Для обеспечения предельного растяжения мастики в шве левые части выражений (1) должны быть равны, что позволяет из равенства правых частей найти ширину паза шва:
По формуле (2) получим ширину паза шва:
в первом районе
во втором районе:
Если швы раскрываются одинаково, то к =1,0; b1 = 5 мм и b2 = 12,5 мм.
СОДЕРЖАНИЕ