Материал: 2196

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам
h 0,585

СТРОИТЕЛЬСТВО. СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

Установим зависимости, характеризующие влияние на продолжительность промерзания

tпр таких теплофизических и температурных

параметров, как:

C и – объемная теплоемкость и тепло-

проводность материала теплового диода; в Tз

и Tзв – температура замерзания материала

телового диода и среднезимняя температура воздуха, соответственно.

Иначе говоря, выявим следующие закономерности:

tпр

f ( );

(1)

tпр

(C);

(2)

tпр

(Tз);

(3)

tпр

1(Tзв );

(4)

tпр

F g ск W Wн ,

(5)

где W и Wн

– объемные влажности (льдисто-

сти) замерзшей и незамерзшей воды материала теплового диода; ск – объемная мас-

са скелета грунта; – удельная теплота плавления льда; g – ускорение свободного падения.

Для этого воспользуемся трансцендентным уравнением [5]

2

Tзв

tпр

 

 

 

 

 

 

 

 

 

 

ск W Wн

0,5 С

 

Tзв

 

 

 

 

 

 

 

 

 

 

 

3,75

 

T0

Tзв

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a W Wн ск

 

 

 

 

 

 

 

 

 

h,

 

 

tл tпр

tл

(6)

где T0

– значение температуры с нулевой го-

довой амплитудой; а – температуропроводность теплового диода; tл – продолжитель-

ность летнего периода; h – высота теплового диода.

Решения tпр будем находить численными

методами, в частности, методом подбора. Для облегчения расчета воспользуемся

следующим приближенным соотношением [5]:

tл tпр tл 13,1 10 3 1,67 10 6 tл tпр (7)

 

 

 

 

1800 tл 3500 час ;

 

 

 

 

 

 

 

 

 

при

 

tпр 2400 час .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда уравнение (6) можно записать в бо-

лее простой для вычислений форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

Tзв

 

tпр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3,75

 

T0

Tзв

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ск W Wн 0,5 С

 

Tзв

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a W Wн ск

 

13,1 10 3

 

 

1,67 10 6tл tпр h.

 

 

 

 

 

(8)

 

Решение

 

 

 

этого

 

уравнения

 

относительно

tпр можно представить как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

2hN M

 

2hN

M

 

h

,

 

 

 

tпр

 

 

 

 

 

(9)

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2N

 

 

 

 

 

 

 

 

2N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

где

 

3,75

 

T0

Tзв

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

13,1 10 3

1,67 10 6tл ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a W Wн ск

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

2

 

Tзв

 

 

tпр

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ск W Wн

0,5 С

 

Tзв

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для проведения вычислительного эксперимента выберем следующие значения параметров. Неварьируемые значения параметров для всех экспериментов: м;

T0 1,206 град.

Для установления зависимостей (1), (2) и

(5) неварьируемыми значениями параметров

будут Tз 0,5

град; Tзв 25 град;

tл 2190 час, а вариация других параметров

будет производиться по следующему алгоритму:

0,25 1,25

 

ккал

 

; шаг 0,25

 

 

 

 

 

 

м час град

(5 шагов);

 

ккал

 

 

 

С 200 600

 

; шаг 100

м3 град

 

 

 

(5 шагов);

g ск W Wн 6500 8500 ккалм3 ;

шаг 500 (5 шагов).

Значения теплофизических параметров материала теплового диода соответствует торфу при различных значениях плотности и влажности.

36

Вестник СибАДИ, выпуск 3 (21), 2011

СТРОИТЕЛЬСТВО. СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

Для определения зависимостей (3) и (4)

значения:

для

зависимости

(3):

неварьируемыми

 

значениями

параметров

Tз ( 0,1) ( 0,5)

град, шаг ( 0,1) (5

ша-

будут:

 

 

 

 

ккал

гов); для зависимости (4): Tзв ( 20) ( 40)

tл 2190

час; 0,75

 

 

 

 

 

;

град, шаг ( 5) (5 шагов).

 

 

 

 

 

 

 

 

 

 

 

м час град

Анализ результатов вычислительного экс-

 

ккал

 

 

 

 

 

 

 

С 400

 

 

 

 

 

 

 

перимента позволил выявить характер зави-

м3 град ;

 

 

 

 

 

 

 

 

 

 

симости

времени

промерзания теплового

g ск W Wн 7500

ккал

 

 

диода от

теплофизических параметров

его

 

м3 ;

 

 

состояния (рисунки 2, 3, 4).

 

 

 

 

 

 

 

 

 

 

 

для зависимости (3): Tзв

25 град; для за-

 

 

 

 

 

висимости (4): Tз

0,7

 

град,

варьируемые

 

 

 

 

Рис. 2. Характер зависимости времени промерзания теплового диода от величин теплопроводности и объемной теплоемкости: 1 – tпр f ( ); 2 – tпр (C).

Рис. 3. Характер зависимости времени промерзания теплового диода от температуры замерзания материала и среднезимней температуры воздуха: 1 – tпр (Tз ); 2 – tпр 1(Tзв ).

Вестник СибАДИ, выпуск 3 (21), 2011

37

СТРОИТЕЛЬСТВО. СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

Рис. 4. Характер зависимости времени промерзания теплового диода от величины, характеризующей интенсивность переходных процессов при промерзании: tпр F g ск W Wн .

Вывод Качественный анализ исследуемых зависимостей показал возможность оптимизации величин и C (рисунок 2) и минимизации величины F (рисунок 4).

Библиографический список

1.Завьялов А. М. Аппарат математического моделирования процессов промерзанияпротаивания грунтов / А. М. Завьялов, Е. А. Бедрин, М. А. Завьялов // Омский научный вестник. – № 3 (93). – 2010. – С. 8-10.

2.Завьялов А. М. Моделирование температурного поля массива многолетнемерзлых грунтов / А. М. Завьялов, Е. А. Бедрин, М. А. Завьялов, В. Н. Лонский // Вестник СибАДИ. – № 3 (17). – Омск: Изд-во СибАДИ, 2010. – С. 49-52.

3.Завьялов А. М. Математическая модель деятельного слоя грунта, функционирующего как тепловой диод / А. М. Завьялов, М. А. Завьялов, Е. А. Бедрин // Омский научный вестник. – № 2 (100). – 2011. – С. 9-13.

4.Завьялов А. М. Влияние теплового диода на мощность подстилающих вечномерзлых грунтов / А. М. Завьялов, М. А. Завьялов, Е. А. Бедрин // Вестник СибАДИ. – № 1 (19). – Омск: Изд-во СибАДИ, 2011. – С. 25-28.

5.Фельдман, Г. М. Методы расчета температурного режима мерзлых грунтов / Г. М. Фельдман.

М.: Наука, 1973. – 254 с.

REVEALING THE DEPENDENCES OF THERMAL DIODE FREEZING TIME FROM THERMOPHYSICAL PARAMETERS

OF ITS CONDITION

E. A. Bedrin, M. A. Zavyalov, A. M. Zavyalov

Qualitative characteristics of thermal diode freezing time dependences from thermophysical and temperature parameters of its condition such, as are established: heat conductivity, volume thermal capacity, intensity of transients. Both temperatures of thermal diode material freezing and average winter air are received.

Завьялов Александр Михайлович – д-р техн. наук, проф., академик РАЕН, проректор по научной работе Сибирской государственной автомо- бильно-дорожной академии. Основное направление научных исследований – динамика рабочих процессов строительных и дорожных машин. Имеет 255 опубликованных работ, в том числе 3

монографии. E-mail: nis@sibadi.org

Завьялов Михаил Александрович – д-р техн. наук, профессор кафедры «Высшая математика» Сибирской государственной автомобильнодорожной академии. Основное направление научных исследований – транспортное строительство. Имеет более 50 опубликованных работ.

E-mail: zavyalov_ma@sibadi.org.

Бедрин Евгений Андреевич – канд. техн. наук, доцент кафедры «Экономика и управление дорожным хозяйством» Сибирской государственной автомобильно-дорожной академии. Основное направление научных исследований – разработка ресурсосберегающих технологий в дорожном строительстве. Имеет более 35 опубликованных работ. E-mail: Bedrin-ea@yandex.ru

38

Вестник СибАДИ, выпуск 3 (21), 2011

СТРОИТЕЛЬСТВО. СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

УДК 625.7

МОДЕЛИРОВАНИЕ НАПРЯЖЕНОГО СОСТОЯНИЯ АРМИРОВАННЫХ ЛЕДОВЫХ ОБРАЗЦОВ-БАЛОК

О.В. Якименко, С.А. Матвеев

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

Ключевые слова: ледовые образцы-балки, армирование, геосинтетические материалы, разрушение, плоское напряженное состояние, функция напряжений, напря- женно-деформированное состояние

Введение

В течение ряда лет в СибАДИ ведутся исследования по широкому кругу вопросов использования различных геосинтетических материалов в дорожных конструкциях [1-3]. При этом некоторые конструктивнотехнологические решения показали перспективные результаты, например – армирование геосетками (георешетками) ледовых переправ.

Основная часть

Для оценки эффективности применения геосинтетических материалов в качестве арматуры для льда был выполнен расчет на прочность ледовых образцов балок. На рисунках 1 и 2 представлены схема действия нагрузок на образец-балку и поперечное сечение балки – соответственно.

Площадь поперечного сечения образца (A) м2 определяется по формуле

A bh,

(1)

где b – ширина поперечного сечения образца, м; h – высота поперечного сечения образца, м.

Площадь основного поперечного сечения

образца без учета арматуры (A1):

 

A1 A A2 ,

(2)

где A2 – площадь поперечного сечения арма-

туры, м2, определяемая по формуле

 

A2 2 2b2 ,

(3)

Здесь b2 – ширина поперечного сечения арматуры, м; 2 – высота поперечного сечения арматуры, м;

На рисунке 2 Y, Z – главные центральные оси площади А; Z1 – главная центральная ось площади А1; С – центр тяжести площади А1; E1 – модуль упругости льда; E2 – модуль упругости арматуры.

При действии нагрузки балка испытывает две стадии напряженно-деформированного состояния.

Рис. 1. Расчетная схема действия нагрузок: F

– сила, приложенная к образцу; l – расстояние между опорами.

Рис. 2. Поперечное сечение образца

Первая стадия развивается от момента нагружения до момента образования первой трещины, нормальной к его продольной оси. Первая трещина образуется в опасном сечении под нагрузкой при достижении растягивающими напряжениями значений, равных расчетному сопротивлению льда на растяжение ( 2 = Rизг). Рассмотрим данную стадию более подробно.

Вестник СибАДИ, выпуск 3 (21), 2011

39

СТРОИТЕЛЬСТВО. СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

Определим положение нейтральной оси Z0 поперечного сечения армированного образца.

Статический момент сечения льда (S1) и арматуры (S2) относительно нейтральной оси равны соответственно

S1 A1y1,

(4)

где y1 – расстояние от нейтральной оси Z0 до оси Z1, м.

S2 A2 y2 ,

(5)

где y2 – расстояние от нейтральной оси Z0 до центра тяжести поперечного сечения арматуры, м:

y

 

 

h

a y

 

y

 

,

(6)

2

 

1

0

 

2

 

 

 

 

здесь a – расстояние от центра тяжести поперечного сечения арматуры до нижней грани поперечного сечения образца, м;

Момент инерции J2 сечения арматуры относительно нейтральной оси Z0:

J

 

 

b 3

A y2

 

 

2

2 2

,

(11)

 

 

 

6

2

2

 

 

 

 

 

 

 

 

 

Момент инерции J1 площади А1 относительно нейтральной оси Z0:

J1 J Ay42 J2 ,

(12)

Модуль упругости и момент инерции приведенного сечения:

E

 

E

J1

E

 

 

J2

,

(13)

 

 

 

 

 

 

прив

 

1 J

 

 

2

 

J

 

 

 

 

J

прив

J

1

 

E2

 

J

2

,

(14)

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Eприв J E1Jприв .

 

 

 

A

h

 

 

 

y

0

 

2

 

 

a

,

(7)

A1

 

 

 

2

 

 

 

Если учесть, что

dA dA dA 0,

Aсж Aраст

тогда статический момент армированного сечения относительно нейтральной оси должен быть равен нулю:

Напряжения в основной части поперечного сечения, состоящей из льда ( h2 y h1):

 

1

 

E1

 

M

y ,

(15)

Eприв

 

 

 

 

J

 

Здесь М – изгибающий момент в опасном сечении:

M

Fl

,

(16)

 

4

 

 

SZ0

0,

 

 

 

 

 

 

(8)

где F – сила, действующая на обра-

 

 

 

 

 

 

 

 

зец, кН; l – расстояние между опорами, м.

Отметим, что в данном случае армирова-

Напряжения в арматуре:

 

 

 

 

 

 

 

 

 

 

 

ние делает рассматриваемое сечение неод-

 

 

 

 

 

E2

 

M

 

 

нородным. Для решения задачи приведем его

 

 

 

 

 

 

 

2

Eприв

 

J y2 ,

(17)

к однородному виду, и тогда нулю должен

 

быть равен приведенный статический момент

Определим

равнодействующие

внутренних

(Sприв):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E2

 

 

сил (рисунок 3), возникающих в растянутой и

S

 

 

A y A y

 

 

0,

 

сжатой зонах поперечного сечения (R2 и R1):

 

прив

 

 

 

1

1

 

2

 

2 E1

 

 

 

1

 

 

 

 

 

 

С учетом (7) и (8) получаем:

 

R1

 

2 1h1 b,

(18)

 

 

 

 

 

 

 

 

 

 

E

 

h

 

 

 

 

A

 

 

R2

1

2h2

 

b.

(19)

 

 

A

 

2

 

 

a

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1 2

 

 

 

 

A1

,

(9)

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A A

E2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Момент инерции J основной части сечения

 

 

 

 

 

 

 

 

 

 

без учета арматуры относительно оси Z:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

bh3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вестник СибАДИ, выпуск 3 (21), 2011