CC BY
43
Тиратурян А.Н., Углова Е.В., Ляпин А.А. Исследование распределения энергии динамического воздействия транспортных средств в слоях нежесткой дорожной конструкции // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 2. - С. 178-194. DOI: 10.15593/perm.mech/2017.2.10
Tiraturyan A.N., Uglova E.V., Lyapin A.A. Studying the energy distribution of the dynamic influences of road transport on the layers of nonrigid pavements. PNRPUMechanics Bulletin, 2017, no. 2, pp. 178-194. DOI: 10.15593/perm.mech/2017.2.10
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 2,2017 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/ini7
001 10.15593/регш.шесЬ/2017.2.10 УДК 531.3:625.7/8
ИССЛЕДОВАНИЕ РАСПРЕДЕЛЕНИЯ ЭНЕРГИИ ДИНАМИЧЕСКОГО ВОЗДЕЙСТВИЯ ТРАНСПОРТНЫХ СРЕДСТВ В СЛОЯХ НЕЖЕСТКОЙ ДОРОЖНОЙ КОНСТРУКЦИИ
А.Н. Тиратурян, Е.В. Углова, А.А. Ляпин
Донской государственный технический университет, Ростов-на-Дону, Россия
АННОТАЦИЯ
Рассмотрены вопросы изучения распределения рассеиваемой (диссипируе-мой) механической энергии, передаваемой на покрытие дорожной одежды при движении по ней автомобильного транспорта. Для численного моделирования процесса передачи энергии была усовершенствована аналитическая модель динамического напряженно-деформированного состояния многослойного полупространства путем введения подвижной системы координат. Изучение распределения энергии осуществлялось для трех разнопрочных дорожных конструкций, для каждой из которых были получены амплитудно-временные характеристики напряжений и деформаций на поверхности слоев покрытия, основания и грунта земляного полотна, на основе которых были построены динамические петли гистерезиса. Был осуществлен анализ площадей динамических петель гистерезиса на поверхности слоев покрытия, основания и грунта земляного полотна, что позволило выявить качественные и количественные зависимости распределения плотности рассеиваемой энергии в слоях дорожной одежды. Установлено, что с увеличением капитальности дорожной конструкции уменьшается плотность энергии, рассеиваемой на ее поверхности, при этом темпы затухания энергии волновых полей в дорожной конструкции, генерируемой воздействием расчетной нагрузки, в значительной степени отличаются в зависимости от капитальности дорожной одежды и от различных типов конструктивных слоев, из которых она состоит. Наибольшее различие при этом определяется свойствами материала, применяемого в качестве слоя основания (укрепленное, неукрепленное). На основе проведенных исследований предложен новый подход к оценке проектного ресурса нежестких дорожных конструкций с позиций энергии, передаваемой на ее поверхность в течение всего срока службы.
© ПНИПУ
О СТАТЬЕ
Получена: 23 февраля 2017 г. Принята: 05 июня 2017 г. Опубликована: 30 июня 2017 г.
Ключевые слова:
дорожная конструкция, динамический гистерезис, модуль упругости, декремент затухания, диссипация энергии, расчетная нагрузка, энергетический подход, многослойное полупространство, аналитическая модель, подвижная система координат.
© Тиратурян Артем Николаевич - кандидат технических наук, старший преподаватель, e-mail: [email protected]
Углова Евгения Владимировна - доктор технических наук, профессор, e-mail: [email protected] Ляпин Александр Александрович - доктор физико-математических наук, профессор, e-mail: [email protected]
Artem N. Tiraturyan - CSc in Technical Sciences, Senior Lecturer, e-mail: [email protected] Evgenia V. Uglova - Doctor of Technical Sciences, Professor, e-mail: [email protected] Alexander A. Lyapin - Doctor of Physical and Mathematical Sciences, Professor, e-mail: [email protected]
STUDYING THE ENERGY DISTRIBUTION OF THE DYNAMIC INFLUENCES OF ROAD TRANSPORT ON THE LAYERS OF NONRIGID PAVEMENTS
A.N. Tiraturyan, E.V. Uglova, A.A. Lyapin
Don State Technical University, Rostov-on-Don, Russian Federation
ABSTRACT
The paper deals with studying the distribution of the scattered (dissipating) mechanical energy transferred to the pavement cover when it is used by road transport. For the numerical simulation of the energy transfer process we improved the analytical model of the dynamic stress-strain state of the multilayered half-space by introducing the moving coordinate system. The energy distribution was studied for three road structures with different strengths. For each structure we obtained the amplitude-time characteristics of stress and strain on the surface of the coating layers, on the base and subgrade soil which have been used to build the dynamic hysteresis loops. We analyzed the areas of the dynamic hysteresis loops on the surface of the coating layers as well as the base and subgrade soil, which made it possible to reveal the qualitative and quantitative dependences of the density distribution of energy dissipation in the pavement layers. It is found that when the solidity of the road structure increases, the energy density which is dissipated on its surface decreases. The rate of the energy attenuation of wave fields in road structures which has been generated by the impact of the design load vary significantly depending both on the solidity of the pavement and various types of its structural layers. In this case the greatest difference is determined by the material properties which are used as the base layer (reinforced, unreinforced). Based on the studies, a new approach of evaluating the design service life of nonrigid road structures in terms of the energy transferred on its surface during its entire service life has been proposed.
© PNRPU
Современные автомобильные дороги представляют собой сложную систему, призванную обеспечить быстрое, комфортное и безопасное передвижение пользователя. Важнейшей частью этой системы является дорожная конструкция, непосредственно воспринимающая на себя нагрузки от транспортного потока и природно-климатических факторов. Основной задачей инженера-дорожника является обеспечение ее прочности и работоспособности, которые закладываются еще на стадии проектирования.
В практике Российской Федерации нормативным документом ОДН 218.046-01 регламентируется трехкритериальный расчет нежестких дорожных одежд, включающий:
- расчет конструкции в целом по допускаемому упругому прогибу;
- расчет по условию сдвигоустойчивости подстилающего грунта и малосвязных конструктивных слоев основания;
- расчет конструкции на сопротивление монолитных слоев усталостному разрушению от растяжения при изгибе.
Однако следует отметить, что основы используемого метода были разработаны в 5060-е гг. XX в. и базируются на приведении многослойных систем, которые и представляют собой дорожные конструкции, к системам, представляющим собой слой, лежащий на однородном полупространстве. В то же время современное развитие методов механико-математического моделирования позволяет решать задачи, касающиеся исследования напряженно-деформированного состояния многослойных сред без каких-либо упрощений.
Фундаментальные основы для решения данного класса задач были заложены в работах ученых И.И. Воровича, В.А. Бабешко, А.В. Белоконя, Е.В. Глушкова и др. [1-13]. Эффективный математический аппарат для решения прикладных задач анализа напряженно-деформированного состояния нежестких дорожных одежд был создан в работах
ARTICLE INFO
Received: 23 February 2017 Accepted: 05 June 2017 Published: 30 June 2017
Keywords:
pavement, dynamic hysteresis, modulus of elasticity, damping constant, dissipation of energy, design load, energy approach, multilayer half-space, analytical model, moving coordinate system.
С.К. Илиополова, М.Г Селезнева, Е.В. Угловой, А.А. Ляпина [14-27] и др. Использование в качестве базиса для проведения расчетных исследований, разработанных авторами математических моделей, позволяет решать широкий круг задач, в частности, касающихся исследования энергетических характеристик в дорожных конструкциях при воздействии движущегося транспорта.
Один из подходов к решению описанных выше задач был предложен в работах профессора А.В. Смирнова [28-30], посвященных разработке энергетического метода расчета нежестких дорожных одежд. Основным критерием разработанного им метода является условие
IП*
к,, > ктуп, (1)
лр
где П*. - суммарная потенциальная энергия разрушения слоев дорожной конструкции;
Кпу - проектное значение обобщенного критерия работоспособности; Лр - работа сил
веса движущегося транспорта; Ку1 - минимальное значение обособленного критерия
работоспособности (составляет 1,2, 0,65, 0,45, 0,25 для различных уровней транспортно-эксплуатационного состояния).
Выгодным отличием предлагаемых в рамках разработанного подхода решений являлся учет динамических упругих прогибов, формирующихся под действием движущейся нагрузки, и динамических модулей упругости материалов слоев нежесткой дорожной одежды, зависящих от температуры и времени приложения нагрузки. В то же время применение современных механико-математических моделей динамического НДС многослойного полупространства позволяет значительно развить энергетический метод с точки зрения анализа рассеивания энергии движущегося транспортного средства в структуре дорожной одежды.
В общем виде энергетический баланс в дорожной конструкции при воздействии нагрузки, движущейся по поверхности с постоянной скоростью, может быть представлен в виде зависимости
¡И
□ 0
ог]а = Еп0В + Ер,5+Е *
где левая часть приведенного выражения определяет энергию, закачиваемую в среду поверхностным источником по пятну контакта □; г^ - напряжения и деформации на поверхности дорожной конструкции; Епов - энергия, уносимая на бесконечность поверхностными волнами; ЕР 8 - энергия объемных волн; Е * - величина энергии, рассеиваемой в структуре
дорожной конструкции за счет проявления вязких свойств и внутренних аномалий.
При моделировании слоистой среды для рассматриваемых сочетаний физических и геометрических параметров слоев показано, что Епов = 0, а также приняты следующие гипотезы:
- все слои дорожной конструкции характеризуются вязкоупругими свойствами;
- все слои являются изотропными;
- слои не ограничены по ширине;
- грунтово-геологическии массив моделируется в виде упругого полупространства. Для изучения внутреннего распределения энергии в структуре дорожной конструкции при движении транспорта с постоянной скоростью V была разработана аналитическая модель динамического напряженно-деформированного состояния многослойного полупространства, базирующаяся на положениях теории упругости и вязкоупругости. Постановка задачи записывается следующим образом.
Пусть область Э представляет собой Ж-слойное упругое полупространство:
й
Э = Д и Э2 и... и Вм, = { у е ^Х г е (г^ ^ Ъг; (Л = О)-/-
г=1
слой 1=2,...,К) (рис. 1).
Упругие свойства сред в Б ., у = 0Д,...,ЛГ описываются плотностью р . и коэффициентами Ламе X/, ц/ или соответственно модулем упругости Е. и коэффициентом Пуассона V .:
ЗА,. +
А,..
v. =
Рис. 1. Многослойное полупространство «дорожная конструкция - грунт» Fig. 1. Multilayer half-space "pavement - soil"
В дальнейшем все параметры, имеющие размерность длины, будут отнесены к характерному линейному размеру а области распределения нагрузки О без изменения в их обозначениях.
Движение среды определяется решением системы уравнений в частных производных:
а v1) (г, г)
V-c( j )(r, t) = p.
dt2
(2)
где о(1 )(г, г) - тензор напряжений Коши; u( 1 )(г, г) - перемещения в точке наблюдения г = {х, у, г} в момент времени г, и( 1) = {ы^), иу1, и[]) } = {м[1), и21), ы( 1) }; индекс / соответствует номеру слоя (компоненты) слоистой конструкции (поиск решения уравнений осуществляется раздельно для каждого 1 с последующим сопряжением решений вдоль границ раздела однородных сред).
Условия стыковки слоев между собой, а также полупространства с вышележащим слоем могут приниматься как жесткими, так и однородными условиями скользящего контакта, по одной или нескольким границам раздела.
Для решения задачи с подвижными возмущениями будем опираться на решение задачи об установившихся гармонических колебаниях с частотой ш в подвижной, связанной с нагрузкой, системе координат. В предположении этого все соотношения для стационарной задачи далее выписаны в амплитудных функциях, в которых временной множитель ехр(-гш г) опущен.
На границе области считаем заданными векторы-функции напряжений:
t (N) = |T(N ) T(N) ст( N =
xz ' yz ' z
} = T (R),
z = z
N
R = (x, y) eQ.
(3)
Система (2) с использованием закона Гука для линейно-упругого материала, связывающего компоненты тензора напряжений а( 1) с компонентами тензора малых деформаций 8( 1)
а( 1) = 2(18( 1) + Л,Е ^8( 1), (4)
может быть сведена к эквивалентной системе уравнений Ламе относительно функций перемещений точек среды:
I2 02
УУ- и(1 )(г )--2г УхУх и( 1) (г ) + 0 ли (1) (г ) = 0, (5)
12
где 0 21=ш2а2/ Ур}-, 0 22 = ш2а2/ У^, - приведенные частоты колебаний,
УР]. = , У5]. = у^Цр! - скорости распространения продольных и попереч-
ных волн в 1-й среде.
Решение системы (5) строится с использованием интегрального преобразования Фурье по координатам х,у и применением принципа предельного поглощения [1, 3]:
и(1 )(г) = -1-- Ц К (а, Р, 2, ш) - Т(а, Р)ехр [-/ах - /'Ру ]ёхёу, (6)
Г1Г2
где Т(а, Р) - преобразование Фурье от нагрузки Т(И); Г . - контур интегрирования, определяемый принципом предельного поглощения.
Введя подвижную систему координат у = у - Vt, воспользуемся принципом соответствия [2, 3], по которому решение задачи с подвижным возмущением со скоростью V определяется численным интегрированием по (6) с заменой в подынтегральном выражении величины ш на ш* = ш + VР .
При наличии диссипации в среде, определяющейся вязкоупругими свойствами материала, воспользуемся принципом упругого - вязкоупругого соответствия, по которому коэффициенты Ламе становятся комплексными, что, в свою очередь, приводит к ком-плекснозначности приведенных частот. В этом случае диссипативные свойства материалов конструкции можно определять тангенсами углов потерь раздельно для продольных
♦ (1) ♦ (1) П ♦ (1) 1т(0л) ♦ (1) 1т(022) и поперечных волн ур', у; . При этом ур' =-—, у^ =-2—.
Яе(0;1) Яе(0; 2 )
Учитывая, что наибольшее проявление диссипации энергии устанавливается в верхних асфальтобетонных слоях дорожной конструкции, где данные величины также существенно зависят от температуры окружающей среды. Для численного определения тангенсов углов потерь примем 5-константную модель вязкоупругого материала [29-30], согласно которой
Е (ш) =- Н----комплексный температурно частотно-зависимый модуль
и 1
1 +-г +-Г"
(/шп)'1 (/шп)
упругости; Е - мгновенный модуль упругости, МПа; п - время релаксации материала; ё, к1, к2 - экспериментальные параметры модели.
Для асфальтобетона принято: Е = 20 000 МПа, ё = 1,7, к1 = 0,2, к2 = 0,5. Величина п зависит от температуры и имеет следующий вид:
( 1 1 V
n(T ) =
Y (T )
Y (T ) = exp
-X
1
T + 273
T0 + 273
где п0 - время релаксации при фиксированной температуре Г0, °С. В частности, можно принять: п0 = 0,1, Г0 = 15 °С, А = 25000.
По заданной величине Е(ю) для у'-го материала пересчитаем далее величины коэффи циентов Ламе и приведенные частоты, считая значения плотности материала и коэффи циента Пуассона постоянными:
Еу (ю)у у Еу (ю)
, Ну =
j
j (1 + v j )(1 - 2v j )
(8)
■у- -у 2(1 + )
На рис. 2 приведены характерные зависимости тангенса угла потерь продольных волн верхнего слоя асфальтобетона при различных температурах от частоты колебаний, Гц.
Путем дифференцирования по координатам выражения (6) и применения закона Гука аналогично виду (6) в форме двойного интеграла Фурье строится решение для напряжений (агу) и деформаций (гу) в структуре дорожной конструкции, при этом изменяется вид подынтегральной матрицы К. Вид полученных матриц для подынтегральных выражений отражен в [15].
Основным направлением теоретических исследований с привлечением данной модели явилось изучение закономерностей распределения рассеивания энергии от движения транспортного потока в структуре до-
tgy„
0,6
0,4
0,2
0
1
\ \ Лч
\ ч \
__
1
0
100
200
300
/, Гц
Рис. 2. Зависимость тангенса угла потерь продольных волн от частоты колебаний tgy^: 1 - при T = 0 °C; 2 - при T = 10 °C; 3 - при T = 20 °C; 4 - при T = 30 °C Fig. 2. The dependence of the loss tangent of the longitudinal waves on the frequency of oscillations Рожнои к°нструкции.
tgy : 1 is for T = 0 °C 2 is for T = 10 °C При проведении численного моделирования 3 is for T = 20 °C; 4 is for T = 30 °C нами рассматривалось динамическое воз-
действие расчетной нагрузки давлением 600 кПа, распределенной по круговому отпечатку диаметром 34 см и движущейся со скоростью 60 км/ч. Расчет плотности рассеиваемой энергии осуществлялся в соответствии с зависимостью
о
W = | Ojjd s jj, при i, j = 1,2,3.
(9)
Физический смысл зависимости (9) сводится к построению расчетной динамической петли гистерезиса на поверхности дорожной конструкции в координатах a — s (напряжение - деформация) для всех компонент напряжений. При этом суммарная плотность рас-
3
сеиваемой энергии W рассчитывается так: W = ^ Wij, где для W^ в формуле (9) суммиро-
i, j=1
вание по индексам не проводится.
о
В рамках численного эксперимента нами были рассмотрены три варианта дорожных конструкций. Первые два варианта характеризуются значительной прочностью, обеспечивающейся за счет применения в их конструкциях современных материалов и слоев значительной толщины. Как правило, эти варианты дорожных конструкций устраиваются на современных автомагистралях. Третий вариант представляет собой традиционную дорожную конструкцию, устраиваемую на дорогах II-III категорий. Варианты дорожных конструкций представлены в табл. 1-3:
Таблица 1
Конструкция дорожной одежды № 1
Table 1
Structure of pavement Nr. 1
Слой Толщина, см Модуль упругости, МПа
Щебеночно-мастичный асфальтобетон ЩМА-15 на ПБВ 60 4 5300
Асфальтобетон плотный тип А I марки на БНД 60/90 8 3200
Асфальтобетон пористый I марки на БНД 60/90 16 2000
Слой из асфальтогранулобетона, укрепленный комплексным вяжущим 20 700
Щебеночно-песчаная смесь, оптимального гранулометрического состава С4 35 240
Грунт земляного полотна-суглинок тяжелый пылеватый 41
Таблица 2
Конструкция дорожной одежды № 2 Table 2
Structure of pavement Nr. 2
Слой Толщина, см Модуль упругости, МПа
Щебеночно-мастичный асфальтобетон ЩМА-15 на ПБВ 60 4 5300
Асфальтобетон плотный тип А I марки на БНД 60/90 8 3200
Асфальтобетон пористый I марки на БНД 60/90 16 2000
Щебеночно-песчаная смесь, оптимального гранулометрического состава С4 55 240
Грунт земляного полотна-суглинок тяжелый пылеватый 41
Таблица 3
Конструкция дорожной одежды № 3 Table 3
Structure of pavement Nr. 3
Слой Толщина, см Модуль упругости, МПа
Щебеночно-мастичный асфальтобетон ЩМА-15 на ПБВ 60 5 5300
Асфальтобетон порситый тип А I марки на БНД 60/90 7 2000
Асфальтобетон пористый I марки на БНД 60/90 7 2000
Щебеночно-песчаная смесь, оптимального гранулометрического состава С4 35 240
Грунт земляного полотна-суглинок тяжелый пылеватый 41
Расчет плотности рассеянной энергии осуществлялся в критических точках в структуре дорожной конструкции:
- на поверхности покрытия;
- на нижней границе пакета асфальтобетонных слоев;
- на поверхности основания;
- на поверхности рабочего слоя грунта земляного полотна.
Пример результатов имитационного моделирования амплитудно-временных характеристик напряжений и деформаций для дорожной конструкции № 1 приведен на рис. 3-4. На рис. 5-8 представлены полученные динамические петли гистерезиса моделируемых нежестких дорожных конструкций в критических точках дорожной конструкции. Построение петель осуществлялось для всех нормальных напряжений и деформаций, действующих в точке приложения нагрузки.
Рис. 3. Нормальные напряжения на поверхности дорожной конструкции № 1 Fig. 3. Normal stresses on the surface of pavement Nr. 1
Время, с
Рис. 4. Нормальные деформации на поверхности дорожной конструкции № 1 Fig. 4. Normal strains on the surface of pavement Nr. 1
Поверхность дорожной конструкции № 1
Поверхность дорожной конструкции № 2
Микродеформация, е^ Микродеформация, eyy
Поверхность дорожной конструкции № 3
Микродеформация, ez.
я и н
е
н
же
я
Микродеформация, ех
Микродеформация, еу
Микродеформация, ez:
Рис. 5. Результаты моделирования эффекта динамического гистерезиса
на поверхности моделируемых дорожных конструкций Fig. 5. The results of simulating the dynamic hysteresis effect on the lower boundary of the package of asphalt-concrete layers
Анализ представленных на рис. 5-8 и в табл. 4 результатов позволяет сделать следующие выводы:
1. С увеличением капитальности дорожной конструкции уменьшается плотность энергии, рассеиваемой на ее поверхности. Так, при изменении значения эквивалентного модуля упругости с 340 МПа (конструкция № 3) до 620 МПа (конструкция № 1) суммарная плотность энергии (Фас, Wyy, Wzz) уменьшается с 223,16 Дж/м3 до 105 Дж/м3.
2. Темпы затухания энергии волновых полей в дорожной конструкции, генерируемой воздействием расчетной нагрузки, в значительной степени отличаются в зависимости от капитальности дорожной одежды и от различных типов конструктивных слоев, из которых она состоит.
3. Для дорожной конструкции № 1 характерно резкое затухание вертикальной составляющей рассеиваемой энергии, что обусловлено высокой жесткостью слоя укрепленного основания. При этом рассеиваемая энергия, направленная на формирование горизонтальных напряжений и деформаций на поверхности основания дорожной конструкции и грунта земляного полотна, превышает аналогичные значения для дорожной конструкции № 2, устроенной на неукрепленном основании.
4. Дорожная конструкция № 3 характеризуется наибольшими значениями плотности рассеиваемой энергии в вертикальном направлении на поверхности всех слоев дорожной конструкции и в горизонтальном направлении (на нижней границе пакета асфальтобетонных слоев), что напрямую связано с низкой капитальностью дорожной конструкции, обусловленной применением менее прочных материалов слоев, чем в конструкциях № 1 и № 2.
и
Нижняя граница пакета асфальтобетонных слоев дорожной конструкции № 1
Микродеформация, е^
1000
is
D
н
е
-566-
-400 -200
J,
200
-500
Микродеформация, еу
н
е
а На
450-
-т-
50
-200
t
О
200 400
-50
Микродеформация, ez:
Нижняя граница пакета асфальтобетонных слоев дорожной конструкции № 2
н
е
пр
а На
н
е
а На
я и н
е
р На
Микродеформация, gx
Микродеформация, gy
Микродеформация, ez:
я и н
е
Нижняя граница пакета асфальтобетонных слоев дорожной конструкции № 3
.<4 1ППП .<4
Микродеформация, gx
я и н
е
Микродеформация, gy
н
же
я
Микродеформация, ez:
Рис. 6. Результаты моделирования эффекта динамического гистерезиса
на нижней границе пакета асфальтобетонных слоев Fig. 6. The results of simulating the dynamic hysteresis effect on the surface of the simulated pavements
Установленные качественные и количественные закономерности распределения рассеиваемой на поверхности дорожной конструкции плотности энергии, генерируемой расчетной нагрузкой от колеса автотранспортного средства при движении, позволяют предположить возможность суммирования рассеиваемой энергии в структуре дорожной конструкции за весь ее расчетный срок службы (Фпр). Сопоставляя полученное значение с неким предельно допустимым значением плотности рассеиваемой энергии (Фдоп), можно сформулировать критерий работоспособности дорожной конструкции в виде
W < W .
пр доп
(10)
н
е
Поверхность основания дорожной конструкции № 1
Микродеформация, ех
¡й о
н
же
я
Я
Микродеформация, ey
Поверхность основания дорожной конструкции № 2
я и н
е
р
К
Микродеформация, ez
-400-
н
же
я
-50-
/
-500
-50
О
500 1000
Микродеформация, ех
Микродеформация, ey
Поверхность основания дорожной конструкции № 3
Микродеформация, ех
Микродеформация, ez
н
же
я
£ -500 К -50
1000
Микродеформация, е^
Микродеформация, ez:
Рис. 7. Результаты моделирования эффекта динамического гистерезиса
на поверхности основания моделируемых дорожных конструкций Fig. 7. The results of simulating the dynamic hysteresis effect on the surface of the base of the simulated pavements
Расчет дорожной конструкции на соответствие данному критерию может быть осуществлен в соответствии с подходом, представленным на рис. 9. Так, величину предельно допустимого значения плотности рассеиваемой энергии (Фдоп) от расчетной нагрузки предлагается определять путем имитационного моделирования проезда транс-
портных средств с расчетной нагрузкой на заднюю ось в течение расчетного срока службы дорожной одежды по поверхности полупространства, характеризующегося минимально требуемым общим модулем упругости, рассчитанным исходя из суммарного числа приложений расчетной нагрузки. После определения допустимого значения осуществляется конструирование проектной дорожной одежды с заданием перечня исходных параметров в соответствии с рис. 9. На заключительном этапе осуществляется проверка выполнимости критерия (10). Влияние температурных и влажностных факторов на механические параметры дорожной конструкции учитывается изменением модулей упругости материалов слоев асфальтобетона и грунта земляного полотна по временам года в процессе расчета.
н
е
а На
я и н
е
Микродеформация, g
Поверхность грунта земляного полотна дорожной конструкции № 1 а
н
же
я
Микродеформация, ew Поверхность грунта земляного полотна дорожной конструкции № 2
Микродеформация, ez:
Микродеформация, gx
Микродеформация, gy
Микродеформация, ez:
н
же
я
Поверхность грунта земляного полотна дорожной конструкции № 3 аа
Микродеформация, gx
н
е
а На
Микродеформация, gy
Микродеформация, ez:
Рис. 8. Результаты моделирования эффекта динамического гистерезиса на поверхности грунта земляного полотна моделируемых
дорожных конструкций Fig. 8. The results of simulating the dynamic hysteresis effect on the surface of the soil subgrade of the simulated pavements
Таблица 4
Результаты расчета рассеиваемой энергии динамического воздействия транспортных средств в структуре моделируемых дорожных конструкций
Computational results of the dissipated energy of dynamic impact of vehicles in the structure of the simulated road structures
Table 4
Критические точки в структуре дорожной конструкции Рассеиваемая энергия (W) в плоскостях x-x, y-y, z-z, Дж/м3
Конструкция № 1 Конструкция № 2 Конструкция № 3
Wxx W Wzz Wxx W Wzz Wxx W Wzz
Поверхность конструкции 43 25 37 61 34 32 129 81 13,16
Нижняя граница пакета а/б 13 9,4 3,84 28,8 21,2 4,45 85,9 61,5 13
Поверхность основания 0,9 0,58 3,84 0,4 0,13 4,45 1,2 0,33 13
Поверхность грунта 0,55 0,38 0,15 0,044 0,011 0,4 0,21 0,05 2,05
В случае если критерий не выполняется, то производится корректировка толщин и материалов слоев дорожной конструкции.
Рис. 9. Расчетный подход к оценке энергетического ресурса нежесткой дорожной конструкции Fig. 9 The computational approach to the assessment of energy resource of the nonrigid pavement
Выводы
1. Проведено имитационное моделирование динамического воздействия расчетной нагрузки (115 кН), движущейся по поверхности нежесткой дорожной конструкции со скоростью 60 км/ч, в пространственной постановке. На основе результатов моделирования установлено, что темпы затухания энергии волновых полей в дорожной конструкции, генерируемой воздействием расчетной нагрузки, в значительной степени отличаются в зависимости капитальности дорожной одежды и от различных типов конструктивных слоев, из которых она состоит.
2. Установлены качественные и количественные зависимости распределения плотности энергии от движущегося транспортного средства в структуре дорожной конструкции. Для дорожной конструкции № 1 характерно резкое затухание вертикальной составляющей рассеиваемой энергии, что обусловлено высокой жесткостью слоя укрепленного основания. При этом рассеиваемая энергия, направленная на формирование горизонтальных напряжений и деформаций на поверхности основания дорожной конструкции и грунта земляного полотна, превышает аналогичные значения для дорожной конструкции № 2, устроенной на неукрепленном основании.
3. Дорожная конструкция № 3 характеризуется наибольшими значениями плотности рассеиваемой энергии в вертикальном направлении на поверхности всех слоев дорожной конструкции и в горизонтальном направлении (на нижней границе пакета асфальтобетонных слоев), что напрямую связано с низкой капитальностью дорожной конструкции, обусловленной применением менее прочных материалов слоев, чем в конструкциях № 1 и № 2.
4. На основе проведенных исследований сформулирован расчетный подход к определению энергетического ресурса нежесткой дорожной конструкции.
Библиографический список
1. Ворович И.И., Селезнев М.Г., Бабешко В.А. Вибрация штампа на двуслойном основании // ПММ. - 1977. - Т. 41. - Вып. 1. - С. 166-173.
2. Белоконь А.В. К теории динамических задач с подвижными возмущениями для неоднородной упругой полосы // Докл. АН СССР. - 1981. - Т. 261, № 5. - С. 1079-1082.
3. Белоконь А.В. Колебания упругой неоднородной полосы, вызванные движущимися нагрузками // Прикл. математика и механика. - 1982. - Т. 46, № 2. - С. 296-302.
4. Бабешко В.А., Глушков Е.В., Зинченко Ж.В. Динамика неоднородных линейно-упругих сред. - М.: Наука, 1989. - 344 с.
5. Боев С.И., Румянцев А.Н., Селезнев М.Г. Решение задачи о возбуждении волн в упругом двухслойном полупространстве // Методы расширения частотного диапазона вибросейсмических колебаний / Ин-т геол. и геофиз. Сиб. отд-ия АН СССР. - Новосибирск, 1987. - 314 с.
6. Ворович И.И., Александров В.М., Бабешко В.А. Неклассические смешанные задачи теории упругости. - М.: Наука, 1974. - 322 с.
7. Гринченко, В.Т. Равновесие и установившиеся колебания упругих тел конечных размеров. -Киев: Наукова думка, 1978. - 264 с.
8. Приварников А.К., Мерзликин А.Е. Об использовании адаптивных программ интегрирования при решении прикладных задач теории многослойных оснований. Новое в проектировании конструкций дорожных одежд // Тр. СоюздорНИИ. - М., 1988. - С. 22-36.
9. Пряхина О.Д. Нестационарные колебания упругой балки на вязкоупругом основании // Изв. АН СССР. МТТ. - 1992. - № 1. - С. 164-169.
10. Слепян Л.И. Нестационарные упругие волны. - Л.: Судостроение, 1972. - 132 с.
11. Глушков Е.В. Расчет сооружений, заглубленных в грунт. - М.: Стройиздат, 1977. - 265 с.
12. Улитко А.Ф. Метод собственных векторных функций в пространственных задачах теории упругости. - Киев: Наукова думка, 1979. - 261 с.
13. Илиополов С.К., Селезнев М.Г., Углова Е.В. Необходимо разработать новые критерии расчета и конструирования дорожных одежд // Наука и техника в дорожной отрасли. - 2000. -№ 3. - С. 13.
14. Илиополов С.К., Селезнев М.Г., Углова Е.В. Динамика дорожных конструкций: моногр. -Ростов н/Д: Изд-во Рост. гос. строит. ун-та, 2002. - 258 с.
15. Исследование динамических характеристик дорожных конструкций / С.К. Илиополов, М.Г. Селезнев, А.Б. Суворов, Е.В. Углова // Наука и техника в дорожной отрасли. - 2002. - № 1. - С. 9.
16. Илиополов С.К., Селезнев М.Г. Разработка математических моделей и исследование на их основе энергетических характеристик воздействия автотранспорта на дорожную конструкцию и распределение колебаний в элементах системы «дорожная конструкция-грунт» // Дороги России XXI века. - 2003. - № 8. - С. 49.
17. Medvedeva T.A., Seleznev M.G., Sobisevich A.L. Analysis of wave fields induced in a layered half-space by a moving load // Doklady Physics. - 1999. - Vol. 44. - No. 7. - С. 450-453.
18. Боев С.И., Селезнев М.Г. Об одном подходе в нестационарных задачах теории упругости // Изв. вузов. Северо-Кавказский регион. Серия: Общественные науки. - 1989. - № 2. - С. 76.
19. Селезнев М.Г., Корабельников Г.Я. Некоторые особенности динамического взаимодействия геологической среды с поверхностными объектами // Изв. Рост. гос. строит. ун-та. - 2004. -№ 8. - С. 202.
20. Маринченко Е.В., Селезнев М.Г. Влияние дефектов на количественные и качественные характеристики колебаний элементов строительных конструкций // Изв. вузов. Северо-Кавказский регион. Серия: Естественные науки. - 2005. - № S1. - С. 17-21.
21. Селезнев М.Г. Аналитические методы в контактных задачах теории упругости // Изв. вузов. Северо-Кавказский регион. Серия: Естественные науки. - 2004. - № 4. - С. 97.
22. Углова Е.В., Илиополов С.К., Селезнев М.Г. Усталостная долговечность эксплуатируемых асфальтобетонных покрытий. - Ростов-н/Д: Изд-во Рост. гос. строит. ун-та, 2009. - 244 c.
23. Углова Е.В. Моделирование деформирования нежестких дорожных конструкций при воздействии движущегося транспорта / Изв. вузов. Строительство. - 2009. - № 3-4. - С. 87-94.
24. Тиратурян А.Н. Анализ программ для «обратного» расчета модулей упругости конструктивных элементов нежестких дорожных одежд на стадии эксплуатации // Строительство и реконструкция. - 2012. - № 4. - С. 72-78.
26. Смирнов A.B. Теоретические и экспериментальные исследования работоспособности нежёстких дорожных одежд: автореф. дис. ... д-ра техн. наук. - М.: Изд-во МАДИ, 1991. - 38 с.
27. Смирнов А.В., Илиополов С.К., Александров А.С. Динамическая устойчивость и расчет дорожных конструкций. - Омск: Изд-во СибАДИ, 2003. - 188 с.
28. Смирнов A.B. Прикладная механика дорожных и аэродромных конструкций: учеб. пособие. - Омск: Изд-во ОмГТУ, 1993. - 360 с.
29. Hürtgen H. Forschung Strassenbau und Strassenverkehrstechnik // Zum viskoelastischen und viskoplastischem Verhalten von Asphalt, 2010. - 210 р.
30. Beckedahl H. Forschung strassenbau und strassenverkehrstechnik. Bundesmin. fur Verkehr, Bau- u. Wohnungswesen. Abt. Strassenbau, Strassenverkehr, 1995, 250 p.
References
1. Vorovich, I.I. Vibratsiia shtampa na dvusloinom osnovanii [Vibration of the die on a double-ply substrate]. Prikladnaia matematika i mekhanika, 1977, vol. 41, iss. 1, pp. 166-173
2. Belokon' A.V. K teorii dinamicheskikh zadach s podvizhnymi vozmushcheniiami dlia neodnorodnoi uprugoi polosy [To the theory of dynamic problems with mobile perturbations for a non-homogeneous elastic strip]. DokladyAkademii naukSSSR, 1981, vol. 261, no 5, pp. 1079-1082.
3. Belokon' A.V. Kolebaniia uprugoi neodnorodnoi polosy, vyzvannye dvizhushchimisia nagruzkami [Oscillations of an elastic inhomogeneous band caused by moving loads]. Prikladnaia matematika i mekhanika. 1982, vol. 46, no. 2, pp. 296-302.
4. Babeshko V.A., Glushkov E.V., Zinchenko Zh.V. Dinamika neodnorodnykh lineino-uprugikh sred [Dynamics of inhomogeneous linearly elastic media]. Moscow, Nauka, 1989, 344 p.
5. Boev S.I., Rumiantsev A.N., Seleznev M.G. Reshenie zadachi o vozbuzhdenii voln v uprugom dvukhsloinom poluprostranstve [Solution of the problem of wave excitation in an elastic two-layer half-space]. Novosibirsk, Sbornik: "Metody rasshireniia chastotnogo diapazona vibroseismicheskikh kolebanii", Institut geologii i geofiziki Sibirskogo otdeleniia akademii nauk SSSR, 1987, 314 p.
6. Vorovich I.I., Aleksandrov V.M., Babeshko V.A. Neklassicheskie smeshannye zadachi teorii uprugosti [Nonclassical mixed problems of the theory of elasticity]. Moscow, Nauka, 1974, 322 p.
7. Grinchenko V.T. Ravnovesie i ustanovivshiesia kolebaniia uprugikh tel konechnykh razmerov [Equilibrium and steady-state oscillations of elastic bodies of finite dimensions]. Kiev, NaukovaDumka, 1978, 264 p.
8. Privarnikov A.K. Ob ispol'zovanii adaptivnykh programm integrirovaniia pri reshenii prikladnykh zadach teorii mnogosloinykh osnovanii. Novoe v proektirovanii konstruktsii dorozhnykh odezhd [On the use of adaptive integration programs in solving applied problems in the theory of multilayered bases. New in the design of pavement designs]. Moscow, Trudy Soiuzdornii, 1988, pp. 22-36
9. Priakhina O.D. Nestatsionarnye kolebaniia uprugoi balki na viazkouprugom osnovanii [Nonstationary vibrations of an elastic beam on a viscoelastic base].Moscow, Izvestiia Akademii nauk SSSR, Mekhanika tverdogo tela, 1992, no. 1, pp. 164-169.
10. Grinchenko V.T. Ravnovesie i ustanovivshiesia kolebaniia uprugikh tel konechnykh razmerov [Equilibrium and steady-state oscillations of elastic bodies of finite dimensions]. Kiev, Naukova Dumka, 1978, 264 p.
11. Slepian L.I. Nestatsionarnye uprugie volny [Nonstationary elastic waves]. Leningrad, Sudostroenie, 1972, 132 p.
12. Glushkov E.V. Raschet sooruzhenii zaglublennykh v grunt [Calculation of structures buried in the ground]. Moscow, Stroiizdat, 1977, 265 p.
13. Ulitko A.F. Metod sobstvennykh vektornykh funktsii v prostranstvennykh zadachakh teorii uprugosti [The method of eigenvector functions in the spatial problems of the theory of elasticity]. Kiev, Naukova dumka, 1979, 261 p.
14. Iliopolov S.K., Seleznev M.G., Uglova E.V. Neobkhodimo razrabotat' novye kriterii rascheta i konstruirovaniia dorozhnykh odezhd [It is necessary to develop new criteria for the calculation and design of road clothes]. Nauka i tekhnika v dorozhnoi otrasli, 2000, no. 3. pp. 5-18.
15. Iliopolov S.K., Seleznev M.G., Uglova E.V. Dinamika dorozhnykh konstruktsii: Monografiia [Dynamics of road structures: Monograph]. Rostov-na-Donu, Rostovskii gosudarstvennyi stroitel'nyi universitet, 2002, 258 p.
16. Iliopolov S.K., Seleznev M.G., Suvorov A.B., Uglova E.V. Issledovanie dinamicheskikh kharakteristik dorozhnykh konstruktsii [Investigation of dynamic characteristics of road structures]. Nauka i tekhnika v dorozhnoi otrasli, 2002, no. 1. pp. 11-20.
17. Iliopolov S.K., Seleznev M.G. Razrabotka matematicheskikh modelei i issledovanie na ikh osnove energeticheskikh kharakteristik vozdeistviia avtotransporta na dorozhnuiu konstruktsiiu i raspredelenie kolebanii v elementakh sistemy «dorozhnaia konstruktsiia-grunt» [Development of mathematical models and research based on them on the energy characteristics of the impact of vehicles on the road structure and the distribution of vibrations in the elements of the system "road construction-ground"]. Dorogi RossiiXXI veka, 2003, no 8. pp. 49-55.
18. Medvedeva T.A., Seleznev M.G., Sobisevich A.L. Analysis of wave fields induced in a layered half-space by a moving load. Doklady Physics, 1999, vol. 44, no. 7, pp. 450-453.
19. Boev S.I., Seleznev M.G. Ob odnom podkhode v nestatsionarnykh zadachakh teorii uprugosti [On an approach in the nonstationary problems of the theory of elasticity]. Izvestiia vysshikh uchebnykh zavedenii. Severo-Kavkazskii region. Seriia: Obshchestvennye nauki, 1989, no. 2, pp. 76-83.
20. Seleznev M.G., Korabel'nikov G.Ia. Nekotorye osobennosti dinamicheskogo vzaimodeistviia geologicheskoi sredy s poverkhnostnymi ob"ektami [Some features of the dynamic interaction of the geological environment with surface objects]. Izvestiia Rostovskogo gosudarstvennogo stroitel'nogo universiteta, 2004, no. 8, pp. 202-212.
21. Marinchenko E.V., Seleznev M.G. Vliianie defektov na kolichestvennye i kachestvennye kharakteristiki kolebanii elementov stroitel'nykh konstruktsii [Influence of defects on the quantitative and qualitative characteristics of vibrations of elements of building structures]. Izvestiia vysshikh uchebnykh zavedenii. Severo-Kavkazskii region. Seriia: Estestvennye nauki, 2005, no. S1, pp. 17-21
22. Seleznev M.G. Analiticheskie metody v kontaktnykh zadachakh teorii uprugosti [Analytical methods in contact problems of the theory of elasticity]. Izvestiia vysshikh uchebnykh zavedenii. Severo-Kavkazskii region. Seriia: Estestvennye nauki, 2004, no. 4, pp. 97-110.
23. Uglova E.V., Iliopolov S.K., Seleznev M.G. Ustalostnaia dolgovechnost' ekspluatiruemykh asfal'tobetonnykh pokiytii [Fatigue life of exploited asphalt-concrete coatings]. Rostov-na-Donu, Rostovskii gosudarstvennyi stroitel'nyi universitet, 2009, p. 244.
24. Uglova E.V. Modelirovanie deformirovaniia nezhestkikh dorozhnykh konstruktsii pri vozdeistvii dvizhushchegosia transporta [Modeling of deformation of nonrigid road structures under the influence of moving vehicles]. Izvestiia vysshikh uchebnykh zavedenii. Stroitel'stvo, 2009, no. 3-4, pp. 87-94.
25. Tiraturian A.N. Analiz programm dlia «obratnogo» rascheta modulei uprugosti konstruktivnykh elementov nezhestkikh dorozhnykh odezhd na stadii ekspluatatsii [Analysis of programs for the "inverse" calculation of the elastic moduli of structural elements of non-rigid pavements during the operation phase]. Stroitel'stvo i rekonstruktsiia, 2012, no. 4, pp. 72-78.
26. Smirnov A.B. Teoreticheskie i eksperimental'nye issledovaniia rabotosposobnosti nezhestkikh dorozhnykh odezhd [Theoretical and experimental studies of the performance of non-rigid pavements.]. Avtoreferat dissertatsii nasoiskanie uchenoi stepeni doktora tekhnicheskikh nauk. Moskva, Moskovskii avtomobil'no-dorozhnyi institut, 1991, 38 p.
27. Smirnov A.V., Iliopolov S.K., Aleksandrov A.S. Dinamicheskaia ustoichivost' i raschet dorozhnykh konstruktsii [Dynamic stability and calculation of road structures]. Omsk, Izdatel'stvo Sibirskoi gosudarstvennoi avtomobil'no-dorozhnoi akademii, 2003, 188 p.
28. Smirnov A.B. Prikladnaia mekhanika dorozhnykh i aerodromnykh konstruktsii: Ucheb. posobie [Applied mechanics of road and airfield structures: Proc. Allowance]. Omsk, Izdatel'stvo Omskogo gosudarstvennogo tekhnicheskogo universiteta, 1993, 360 p.
29. Hurtgen H. Forschung Strassenbau und Strassenverkehrstechnik. Zum viskoelastischen und viskoplastischem Verhalten von Asphalt, 2010 210 p.
30. Beckedahl H. Forschung strassenbau und strassenverkehrstechnik. Bundesmin. fur Verkehr, Bau- u. Wohnungswesen. Abt. Strassenbau, Strassenverkehr, 1995, 250 p.
