УДК 625.072:531.8
Е. В. Вайнштейн, В. М. Вайнштейн, П. А. Нехорошков
ИССЛЕДОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ДОРОЖНОЙ ОДЕЖДЫ ЛЕСОВОЗНОЙ АВТОДОРОГИ ПО КАСАТЕЛЬНЫМ НАПРЯЖЕНИЯМ МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ
Приведены результаты исследований напряженно-деформированного состояния конструкции дорожной одежды лесовозной автомобильной дороги по касательным напряжениям методом конечных элементов с использованием компьютерного моделирования по программному комплексу «РLAXIS». Параметры напряженно-деформированного состояния конструкции дорожной одежды лесовозной автомобильной дороги определены при размещении двух лесовозных автопоездов по полосам движения. Данные расчетов позволяют сравнить их с данными, полученными по отраслевым нормам дорожного хозяйства расчета нежестких дорожных одежд.
Ключевые слова: щебеночно-мастичный асфальтобетон, лесовозный автопоезд, конструкция дорожной одежды.
Введение. Для автомобильного лесотранспорта, как и для всего лесопромышленного комплекса в целом, основной задачей развития является дальнейшее повышение эффективности производства. Дорожные условия оказывают значительное влияние на технико-эксплуатационные показатели работы автомобильного лесотранспорта. Дорожная одежда лесовозной автомобильной дороги является сложной многослойной инженерной конструкцией, на которую воздействуют нагрузки от лесовозных автопоездов и других транспортных средств [1].
Целью работы является определение параметров напряженно-деформированного состояния конструкции дорожной одежды лесовозной автомобильной дороги при размещении двух лесовозных автопоездов по полосам движения при компьютерном моделировании методом конечных элементов по программному комплексу «PLAXIS», а также установление коэффициента прочности на сдвиг в слоях земляного полотна и сравнение его с расчетным по отраслевым дорожным нормам (ОДН 218.046 - 2001) [2].
Решаемые задачи, направленные на достижение цели, заключаются в определении коэффициента прочности на сдвиг, в слоях земляного полотна путем моделирования и определения параметров напряженно-деформированного состояния дорожной одежды по касательному напряжению.
Основные уравнения имеют следующий вид.
Рассмотрим силу F (рис.1.), действующую на упругое полупространство, которое имеет модуль упругости Е и коэффициент Пуассона ц. Выберем произвольную систему координат, в которой ось Z направлена вниз по направлению действия силы, а две другие оси расположены в плоскости действия силы. Координаты точки приложения силы обозначим х, у, z, причём z = 0.
© Вайнштейн Е. В., Вайнштейн В. М., Нехорошков П. А., 2012.
Рис. 1. Расчетная схема действия сосредоточенной силы на полупространство и напряжения в точке М
Пусть имеется некоторая точка М с координатами хт, ут, zm. Расстояние от точки М до точки приложения силы обозначим R и найдём по формуле
Напряжения в точке М вычисляются по формулам Буссинеска:
12 ■ -Z 2 m
(1)
3F z3
>5 '
а„ = -
2л R
3F\(Xm - X)2 Zm , 1 - 2Ц
- + -
а = — y 2л
т„ = —
2л [ R5 3
3F f(*m — У)" Zm , 1 — 2М
[ R5
3F (Xm — x)zl
1
(2R + Zm)(Xm — X)2 , Z,
+
3
R(R + Zm) 1
(R + Zm)' R3
+
R3
(2R + Zn)(yn — y)2 , Z,
R(R + Zm) (R + Zm ) R
+
R3
(2)
t = —-Zy 2л
= — 3F TXy = — 2Л
2л Rs
3F (Уп — y)Z,
R
(Xm — X )(Уm — y)Zm 1 — 2И (2R + Zm)(Xm — X)(ym — У)
R5
3
(R + Zm )) R3
От колеса лесовозного автомобиля на полупространство действует распределённая по некоторой области ю нагрузка с интенсивностью q(x,y) (рис.2). Выделим в этой области элементарную площадку da = dxdy.
Действующая на неё сила будет равна
dF = q • dю = q • dxdy. (3)
Одновременно с напряжением от колес лесовозного автомобиля на грунт земляного полотна действует распределенная нагрузка от веса дорожной одежды. На выделенную элементарную площадку da=dxdy действует сила
dP = Ё у" Zi •da = Ё y1Z1 ■ dXdy>
(4)
где у1 - плотность (Г) слоя дорожной одежды, кг/м ; zi - толщина слоя, м.
а = —
Z
Рис. 2. Расчетная схема действия распределенной нагрузки от подвижного состава и веса дорожной
одежды на полупространство (В)
Интенсивность нагрузки от дорожной одежды определяется от каждого слоя в зависимости от плотности его материала и толщины. На поверхности грунта земляного полотна она равна
£ = У ■ ^ +72 ■ ^ + Уз ■ z3- (5)
где у1 - плотность асфальтобетона, кг/м3; z1 - толщина слоя асфальтобетона, м; у2 -плотность укрепленных цементом отсевов дробления известняков, кг/м3; Z2 - толщина слоя укрепленных цементом отсевов дробления известняков, м; у3 - плотность песка, кг/м3; z3 - толщина песчаного слоя, м.
Используя принцип суперпозиции сил и интегрируя по области ю, найдём полные значения напряжений [3]. 3z^ С^^у
о, = —
2л
3
/Р
о = —-
2л 3
о =—2Л11 9
Я<
а
Я<
Д5 '
(хт — х)2 гт 1 — 2м
1
(2Д + Zm)(Xm — X)2 ,
Д5 3
(Ут — У)2 ^ , 1 —
д
- + -
3
Д(Д + Zm) (Д + Zm)2Д3 Д
1 — (2Д + Zm)(Уm — У)" + ^
Д(Д + Zm) (Д + Zm)2Д3 Д
dxdy dxdy
(6)
= —
3z
_m_
2л
flq
— х)
dxdy;
3z 2
,(Уm — У)
dxdy;
=—2Л119
Я<
(^ — x)(Уm — УК 1 — 2^( 2Д + — x)(Уm — У)
3
(Д + z )2 Д3
Аналогично определяем полные значения напряжений от веса дорожной одежды
3
= -
3zm ^ CCYtZidxdy
2п i=1 3 ^
Z ff
2n
= -
r„. = -
Z 1 R5
а ^
' V ffv Z (Xm - X)
R5 '
frx„ - x/Zm . 1 -
- + •
3
1
(2R + zm)(xm -x)2 + z^ R3
R(R + Zm) (R + zm)2R3
> dxdy;
3Zm Z £
dxdy;
(7)
3Z
m
Z Я™ ^dxdy;
2ж t
=-£ Z ih
(xm - x)(ym - У)Zm 1 - 2V (2R + Zm )(xm - x)(ym - У)
R5
3
i2 r»3
(R + Zm;2 R
dxdy
Под действием заданных нагрузок в теле появляются напряжения. При анализе напряжений в окрестности рассматриваемой точки М (рис.3) выделим бесконечно малый элемент в форме параллелепипеда, ребра которого параллельны координатным осям, а их длина равна ёх, ёу, dz.
На гранях этого параллелепипеда действуют три нормальных и шесть касательных напряжений, совокупность которых образует тензор напряжения. Одноименные и параллельные напряжения, действующие на параллельных гранях бесконечно малого параллелепипеда, отличаются на бесконечно малую величину и потому их можно считать одинаковыми. Полученные зависимости позволяют найти наибольшие нормальные и касательные напряжения для рассматриваемой точки.
Среди бесчисленного множества площадок, которые можно провести через исследуемую точку, имеются три взаимно перпендикулярные площадки, касательные напряжения на которых отсутствуют. Эти площадки и действующие по ним нормальные напряжения будут главными [4]. В общем случае пространственного напряженного состояния связь между главными напряжениями и напряжениями по любой наклонной площадке, проходящей через данную точку, описывается кубическим уравнением [5]
а3 - 11а2 + 12а - 1з = 0, (8)
где 11 ,12 ,13 - соответственно первый, второй, третий инварианты тензора напряжений.
II = ах + ау + аz ; (9)
2 2 2
12 = ахау + ау аz + аzах - т ху - т уz - т zх; (10)
222
13 = ах ау аz +2тху туzтzх - ах т уz - аут zх - аzт ху. (11)
Решив это уравнение, найдем главные напряжения аь а2; а3. Кубическое уравнение решаем методом половинного деления.
Рис. 3 . Схема действия напряжений на гранях параллелепипеда в исследуемой точке элемента
т
Все напряжения отличны от нуля, и напряженное состояние является объемным, трехосным.
По площадкам, наклоненным под углом 45° к главным, возникают наибольшие касательные напряжения. Значения этих напряжений определяются по формуле
" " ■ (12)
i max 2
При а1>а2 > а3 Tmax является максимальным касательным напряжением для данной точки тела, где площадка параллельна о2.
Используя условия третьей теории прочности, найдены максимальные касательные напряжения. Условие прочности будет иметь вид
т = < [Тпр], (13)
t max 2
где [тпр] - предельная величина активного напряжения сдвига (в той же точке), превышение которой вызывает нарушение прочности на сдвиг.
Предельную величину активного напряжения сдвига в грунте земляного полотна определяется по теории прочности Мора-Кулона.
т <
i max
пр
или < С + Gmtgф . (14)
В случае напряженно-деформированного состояния оценка прочности в данной точке конструкции производится путем непосредственного сопоставления возникающего в ней рабочего напряжения с предельным. Коэффициент запаса прочности равен отношению предельного напряжения к максимальному:
(15)
Кпр
Т тах
Рассмотрим влияние параметров лесовозных автопоездов и их расстановку на автодороге при определении максимальных касательных напряжений в покрытии, основании автодороги, рабочем слое земляного полотна и в теле насыпи на глубине 84200 см от поверхности покрытия. Расстановку лесовозных автопоездов в сечении автомобильной дороги произведем таким образом, когда две самые нагруженные вторые оси тягача будут находиться в одном сечении по оси Z. Схема расстановки двух лесовозных автопоездов приведена на рис.4. Нагрузка на вторую ось принята с коэффициентом перегруза, равным Кп = 1,4.
Вариант расстановки двух лесовозных автопоездов следующий: правая сторона проезжей части автодороги - лесовозный автопоезд МАЗ 509А+ГКБ-9383-011, груженый с нагрузкой на вторую ось 178 кН [2]. Левая сторона проезжей части автодороги -лесовозный автопоезд МАЗ 509А+ГКБ-9383-011, груженый с нагрузкой на вторую ось 178 кН.
Выберем произвольную систему координат, у которой ось Z направлена вдоль трассы по ходу движения по правой полосе, ось X перпендикулярно ей, ось Y вниз, плоскость XOZ на поверхности покрытия.
Автомобили двигаются по правой и левой полосам. Введём местную систему координат относительно автомобиля. На рис. 4 показана схема определения координат расстановки лесовозов. Параметры и координаты отпечатков колес лесовозных автопоездов по вариантам их расстановки приведены в таблице.
47,50м
50,00м
52,50м
Обочина
1 ось 2 ось
7,м
3 ось 4 ось
Х=47,50м Х=47,55м
Z46,6te Е^Х=47,50м^=6040^ ^Х4^
I I 7=50,00- , I I Z=61,75м
"Х=49,5См | | Х=4Я45м ~ (
' Z=50,00м Z=61,75м
Х=50,55м Х=50,50м
c^aX=50,55мZ=5000M I I Z=39,60м
Х=52,45м | | _ Х=52,50м
7=38,255^^ ^Х=52,45м7=50,00м _7=39,60м
"-Т-Т-м
Х=50,50м 7=53,40м_
Х=52,50м 2=53Г40м
Обочина
^ Х, м 40,00м 50,00м 60,00м
Рис.4. Схема расстановки двух лесовозных автопоездов Параметры и координаты отпечатков колес лесовозных автопоездов
Номер лесовоза Номер оси Диаметр отпеч. колеса,™ Нагрузка на колеса, кН Координаты колес автопоезда, м
Левая полоса движения Правая полоса движения
Правые колеса Левые колеса Левые колеса Правые колеса
X 7 X 7 X 7 X Z
1 2 3 4 5 6 7 8 9 10 11 12
Вариант 1
1, 2 1 26 32,5 47,50 46,60 49,50 46,60 50,50 53,40 52,50 53,40
2 43 89,18 47,50 50,00 49,50 50,00 50,50 50,00 52,50 50,00
3 36 61,75 47,55 60,40 49,45 60,40 50,55 39,60 52,45 39,60
4 36 61,75 47,55 61,75 49,45 61,75 50,55 38,25 52,45 38,25
Исследования параметров напряженно-деформированного состояния конструкции дорожной одежды лесовозной автомобильной дороги после расстановки лесовозных автопоездов и определения его параметров произведены методом конечных элементов.
В настоящее время существует ряд программ, реализующих данный метод. В работе был использован наиболее адаптированный к дорожной специфике программный комплекс «PLAXIS». Расчетный комплекс «PLAXIS» представляет собой пакет прикладных вычислительных программ для конечно-элементного анализа напряженно-деформированного состояния системы основание-фундамент-сооружение в условиях плоской и осесимметричной задач [6].
Исходными для расчета послужили следующие данные:
1 слой - асфальтобетонное покрытие из ЩМА, толщиной 0,06 м. Модуль упругости асфальтобетона при +20о С равен Еупр = 1477 МПа [6]; коэффициент Пуассона ц=0,1; плотность р = 25,10 кН/м3;
2 слой - асфальтобетонное покрытие из пористого асфальтобетона М-1 толщиной 0,06 м. Модуль упругости асфальтобетона при +20о С равен Еупр = 800 МПа [2]; коэффициент Пуассона ц=0,1; плотность р = 24,30 кН/м3;
3 слой - щебеночное основание толщиной 0,38 м. Модуль упругости щебня равен Еупр = 250 МПа [2]; коэффициент Пуассона ц=0,15; плотность р = 20,0 кН/м3;
4 слой - песчаный слой толщиной 0,34 м. Модуль упругости равен Еупр = 100 МПа; коэффициент Пуассона ц=0,25; плотность р = 16,5 кН/м3, сцепление С = 5 кПа (5 кН/м2), угол внутреннего трения ф = 420 [2];
5 слой - грунт суглинок легкий. Модуль упругости равен Еупр = 43 МПа; коэффициент Пуассона ц=0,35; плотность р = 18,0 кН/м3, сцепление С = 4,7 кПа (4,7 кН/м2), угол внутреннего трения ф = 260.
На основании проведенных расчетов при расстановке лесовозов получены следующие данные по определению максимальных касательных напряжений.
На рис.5. показана схема расстановки двух лесовозов и нагрузки, действующей от колес на покрытие.
Схема деформирования дорожной одежды приведена на рис. 6.
Рис. 5. Расчетная схема метода конечных элементов: а - большими стрелками показано положение наиболее нагруженной задней оси лесовоза; б - маленькими стрелками показано положение всех остальных, менее нагруженных осей
Рис. 6. Схема деформирования конструкции дорожной одежды
Схема деформирования дорожной одежды в опасном сечении приведена на рис. 7. Изолинии касательных напряжений в конструкции представлены на рис.8.
Рис. 7. Схема деформирования дорожной одежды в опасном сечении
Рис. 8. Изолинии касательных напряжений в конструкции
На основании построенных изолиний касательные напряжения точек полупространства от нагрузки двух лесовозных автопоездов приведены на рис. 9.
14
<N
I '2 S 1°
И
I 8
& r
5 6
45 46 47 48 49 50 51 52 53 54 Координаты поперечника дорожного полотна
• асфальтобетонное покрытие рабочий слой грунта земполотна
"щебеночное основание грунт земляного полотна
Рис. 9. График максимальных касательных напряжений от нагрузки на вторую ось
лесовозных автопоездов
Анализ рис. 9. показывает, что наибольшие касательные напряжения имеют место под колесами второй оси лесовозного автопоезда. Касательные напряжения от действия нагрузки под колесами второй оси тягача на левой и правой полосе движения в асфальтобетонном покрытии составляют Тф1 = 0,013 МПа, на кромках покрытия и между колесами - Тф1 = 0,009 МПа. Соответственно касательные напряжения на кромке покрытия в 1,4 раза меньше, чем под колесами перегруженного тягача. По глубине конструкции дорожной одежды касательные напряжения уменьшаются. Так, в щебеночном основании на глубине 0,12 м от поверхности покрытия касательные напряжения составляют Тф2= 0,011 МПа под колесами второй оси лесовозного автопоезда. На уровне рабочего слоя грунта земляного полотна на глубине 0,84 м от поверхности покрытия касательные напряжения активного напряжения сдвига составляют Тф3= 0,008 МПа.
По толщине конструкции дорожной одежды с увеличением глубины от поверхности покрытия касательные напряжения уменьшаются. Касательные напряжения на уровне щебеночного основания уменьшаются в 1,2 раза, на уровне верха рабочего слоя земляного полотна - в 1,6 раза, в подошве земляного полотна - в 4,3 раза.
По данным расчета конструкции дорожной одежды по сдвигу в грунте земляного полотна по [2], расчетное активное напряжение сдвига равно Тф.расч= 0,0123 МПа, а пре-
дельная величина активного напряжения сдвига т
пр.расч
= 0,0148 МПа.
Коэффициент прочности по сдвигу в грунте земляного полотна расчетным методом
составляет:к = -
.пред
0,0148
= 1,20.
пр тф 0,0123
ф.расч. '
Фактические значения касательных напряжений на уровне рабочего слоя грунта земляного полотна, рассчитанные по расчетному комплексу «PLAXIS», составляют Тф.эксп= 0,008 МПа, что меньше, чем расчетные по [2] Тф.расч= 0,0123 МПа. Коэффициент прочности по сдвигу в грунте земляного полотна, полученный в сравнении по расчетному комплексу «PLAXIS» и отраслевым нормам, составляет:
0,0123
К =
пр J
ф.расч
ф.эксп.
0,008
= 1,53.
4
2
0
Относительно предельного состояния сдвига коэффициент прочности составляет Кпр = 1,85, что показывает большой запас прочности по сдвигу в грунте земляного полотна в сравнении с касательными напряжениями, возникающими от построенной модели, и расчетами по отраслевым нормам [2].
Выводы
1. Анализ расчета касательных напряжений по расчетному комплексу «PLAXIS» показал: а) максимальные касательные напряжения возникают под колесами второй оси автопоезда; б) влияние соседних колес автопоездов, расположенных на разных полосах движения, не выявлено; в) по толщине конструкции дорожной одежды с увеличением глубины от поверхности покрытия касательные напряжения уменьшаются. На уровне щебеночного основания, глубиной 12 см от покрытия - в 1,2 раза, в рабочем слое земляного полотна, на глубине 84 см от покрытия - в 1,6 раза, в подошве земляного полотна на глубине 144 см - в 4,3 раза.
2. Касательные напряжения, полученные по расчетному комплексу «PLAXIS», для неблагоприятной расстановки лесовозных автопоездов оказались в 1,53 раза меньше, чем рассчитанные по отраслевым нормам [2].
3. Касательные напряжения, полученные по расчетному комплексу «PLAXIS», в сравнении с предельным значением активного напряжения сдвига, оказались в 1,85 раза меньше, что показывает большой запас прочности по сдвигу в грунте земляного полотна.
Список литературы
1. Радовский, Б. С. Проектирование дорожных одежд для движения большегрузных автомобилей / Б. С.Радовский, А. С. Супрун, И. И. Козаков. - Киев: Будивэльнык, 1989. - 168 с.
2. Инструкция по проектированию нежестких дорожных одежд ОДН 218.046-01. - М.: Гос. служба дор. хоз-ва м-ва транспорта РФ, 2001. - 145 с.
3. Вайнштейн, В.М. Разработка конструкции и технологии строительства лесовозных дорог с использованием отсевов дробления известняков (для условий РМЭ): автореф. дисс. ... канд. техн. наук / В. М. Вайнштейн. - Йошкар-Ола, 2002.-22 с.
4. Ицкович, Г.М. Сопротивление материалов / Г.М.Ицкович. - М.: Высш. шк., 1960. - 530 с.
5. Александров, А.В. Основы теории упругости и пластичности: Учебник для строит. спец. вузов /А.В.Александров, В.Д.Потапов. - М.: Высш.шк., 1990. - 400 с.: ил.
6. Вайнштейн, Е. В. Технология строительства лесовозных дорог из щебеночно-мастичных асфальтобетонов с отсевами дробления известняков: автореф. дисс. ... канд. техн. наук / Е. В. Вайнштейн. -Йошкар-Ола, 2010.-16 с.
Статья поступила в редакцию 04.06.12.
Работы выполнены в рамках государственного задания Министерства образования и науки РФ ФГБОУ ВПО «Марийский государственный технический университет» на 2012 год (шифр проекта 5.3516.2011, № 26).
E. V. Vainshtein, V. M. Vainshtein, P. A. Nekhoroshkov
RESEARCH OF STRAIN-STRESS STATE OF ROAD STRUCTURE OF FOREST MOTOR ROAD BY TANGENTIAL STRESS WITH THE USE OF FINITE-ELEMENT
METHOD
Research results of strain-stress state of road structure of forest motor road by tangential stress are presented. A finite-element method with the use of computer simulation in software solution «PLAXIS» was used in the course of research. Characteristics of strain-stress state of road structure offorest motor road are determined with the help of two timber trucks placed at the traffic lanes. Analytical data make it possible to compare the results with the road facilities trade standards data of flexible pavement calculation.
Key words: rubble and mastic asphalt concrete, timber truck, road structure.
ВАЙНШТЕЙН Виктор Мейлехович - кандидат технических наук, доцент кафедры автомобильных дорог Поволжского государственного технологического университета (Россия, Йошкар-Ола). Область научных интересов - рациональное лесопользование, дорожно-строительные материалы и технологии. Автор 36 публикаций.
E-mail: [email protected]
НЕХОРОШКОВ Петр Аркадьевич - кандидат технических наук, доцент кафедры автомобильных дорог Поволжского государственного технологического университета (Россия, Йошкар-Ола). Область научных интересов - рациональное лесопользование, транспорт леса, оптимизация и моделирование транспорта леса, геоинформационные и логистические системы. Автор 38 публикаций.
E-mail: [email protected]
ВАЙНШТЕЙН Евгений Викторович - кандидат технических наук, доцент кафедры автомобильных дорог Поволжского государственного технологического университета (Россия, Йошкар-Ола). Область научных интересов - рациональное лесопользование, дорожно-строительные материалы и технологии. Автор 16 публикаций.
E-mail: [email protected]