- •1. Теоретические основы математического моделирования
- •Проблемы в управлении дорожным движением
- •1.2.Методы моделирования транспортных систем
- •1.3. Классификация моделей
- •1.4. Классификация математических моделей
- •1.5. Этапы разработки математических моделей
- •1.6. Вопросы для самоконтроля
- •2. Основные характеристики движения транспортных средств
- •2.1. Параметры транспортного потока
- •2.2. Макроскопические модели транспортного потока
- •2.3. Микроскопические модели транспортного потока
- •2.4. Вопросы для самоконтроля
- •3. Движение транспортных средств на перекрестках
- •3.1. Типы перекрестков. Анализ конфликтных точек
- •3.2. Определение сложности и опасности перекрёстка
- •3.3. Круговое движение на пересечениях
- •3.4. Управляющие параметры светофорной сигнализации
- •3.5. Пропускная способность регулируемого перекрестка
- •3.6. Задержки автомобилей при регулярном транспортном потоке
- •4.7. Вопросы для самоконтроля
- •4 250 . Транспортные заторы
- •4.1. Одиночные пробки
- •4.2. Серия пробок
- •4.3. Классификация фаз потока
- •4.4. Причина затора – геометрические особенности дороги
- •4.5. Задержки транспортных средств
- •4.6. Множественные устойчивые состояния и хаос
- •4.7. Вопросы для самоконтроля
- •5. Моделирование тягово-скоростных свойств транспортных средств
- •5.1. Силы, действующие на автомобиль при прямолинейном движении
- •5.2. Уравнение прямолинейного движения автомобиля
- •5.3.Динамика автомобиля, транспортирующего емкости с жидкостью
- •5.3.1 Динамика торможения автомобиля при движении по прямой
- •5.3.2. Динамики равномерного движения автомобиля по повороту
- •5.3.3. Расчет дополнительных усилий, вызванных движением жидкости в баках
- •5.4. Контрольные вопросы
- •6. Методы моделирования транспортных средств на микроуровне
- •6.1. Основы построения математических моделей на микроуровне
- •6.2. Моделирование тепловых систем
- •6.3. Технология вычислительного эксперимента
- •6.3.1. Основы метода сеток
- •6.3.2. Схемы аппроксимации уравнения теплопроводности
- •6.3.3. Решение сеточных уравнений
- •Основы теории подобия и моделирования
- •6.4.1. Подобие физических явлений
- •6.4.1.Числа подобия при моделировании обтекания автомобиля
- •Вопросы для самоконтроля
- •7. Основные понятия стохастического моделирования
- •7.1. Моделирование в условия неопределенности
- •7.2. Функция и плотность распределения случайной величины
- •7.3.Меры положения и рассеяния кривой распределения
- •7.4.Теоретические законы распределения
- •7.4.1. Закон нормального распределения (закон Гаусса)
- •Кривая Гаусса имеет следующие особенности.
- •7.4.2. Экспоненциальное распределение
- •7.5. Представление распределения скоростей автомобилей в транспортном потоке
- •7.6. Основы корреляционного и регрессионного анализа
- •7.6.1. Метод наименьших квадратов
- •7.6.2. Выборочный коэффициент корреляции
- •7.7. Вопросы для самоконтроля
- •8. Материалы для самостоятельной работы
- •8.1. Курсовая работа
- •8.2. Методические рекомендации по выполнению курсовой работы
- •Список рекомендуемой литературы
- •8.3. Пример выполнения курсовой работы
- •8.4. О приближенных вычислениях
- •8.5. Тест для проверки уровня обученности
- •4) Сдвиг фаз относительно соседних перекрестков.
- •Пример титульного листа курсовой работы
2.2. Макроскопические модели транспортного потока
Транспортный поток можно рассматривать в соответствии с гидродинамической моделью как поток одномерной сжимаемой жидкости, допуская, что поток сохраняется и существует однозначная зависимость между скоростью и плотностью транспортного потока.
Первое допущение выражается уравнением неразрывности. Второе – функциональной зависимостью между скоростью и плотностью для учета уменьшения скорости движения автомобилей с ростом плотности потока. Это интуитивно верное допущение теоретически может привести к отрицательной величине плотности или скорости. Поскольку одному значению плотности может соответствовать несколько значений скорости, для второго допущения средняя скорость потока в каждый момент времени должна соответствовать равновесному значению при данной плотности автомобилей на дороге. Равновесная ситуация – чисто теоретическое допущение и может наблюдаться только на однополосных участках дорог без пересечений и при движении автомобилей без обгонов.
Закон сохранения транспортного потока
При равновесном движении с плотностью автомобилей k(x,t) в направлении x число автомобилей в интервале длины дороги (x1, x2) в момент времени t равно
(2.13)
Пусть v(x, t) – скорость автомобилей в точке x в момент времени t. Число проходящих через x (единицу длины) автомобилей в момент t равно k(x,t)v(x, t). Найдем уравнение изменения плотности. Число автомобилей в интервале (x1, x2) за время t изменяется в соответствии с числом въезжающих и выезжающих машин:
(2.14)
Интегрируя по времени и полагая, что k и v – непрерывные функции, получим
(2.15)
Поскольку x1, x2 и t1, t2 > 0 произвольны,
(2.16)
Это уравнение дополняется начальным условием
(2.17)
Найдем уравнение для скорости v. Положим, что v зависит только от плотности k. Если дорога пуста (k =0), автомобили едут с максимальной скоростью v=vmax. При наполнении дороги скорость падает вплоть до полной остановки (v=0), когда машины расположены "бампер-к-бамперу" (k=k max). Эта простейшая модель выражается следующим линейным соотношением
(2.18)
Тогда уравнение (2.16) принимает вид
(2.19)
и является законом сохранения количества автомобилей. Действительно, интегрируя (2.19) по x, получим
(2.20)
и, следовательно, количество автомобилей на отрезке дороги (x1, x2) постоянно для любых значений t ≥ 0.
Уравнение (2.18), впервые полученное Гриншилдсом, носит его имя. Линейная аппроксимация Гриншилдса представлена на рис. 2.1.
Можно построить макроскопическую модель, в которой уравнение Гриншилдса является частным случаем. Рассмотрим связь между скоростью v и плотностью k автомобилей на дороге. В общем случае, когда плотность k повышается, водители снижают скорость и наоборот, поэтому
(2.21)
где x(t) – координата движения элемента потока.
Проследим изменение скорости для некоторого передвигающегося элемента потока во времени, которое определяется как полная производная по времени
(2.22)
Из уравнения (2.16) следует соотношение
(2.23)
которое после подстановки в (2.23) принимает вид
(2.24)
Поскольку в соответствии с уравнениям (2.21)
(2.25)
из соотношения (2.24) получаем
(2.26)
где
,
а отрицательный коэффициент
пропорциональности
можно
интерпретировать как вязкость в жидкости.
Для классической сжимаемой жидкости
уравнение (2.26) называется уравнением
Эйлера, в этом случае:
(2.27)
где С – неотрицательная константа с размерностью скорости.
Принято рассматривать более общий класс моделей, в которых
(2.28)
Уравнение
(2.27) соответствует случаю n≠–1,
следовательно, из уравнений (2.26) и (2.27)
.
Решением этого уравнения будет
(2.29)
при n=–1 и
(2.30)
при n≠–1.
Формула (2.29) впервые получена Гринбергом и носит его имя. Обозначив v0 – скорость при k=0, для значений n≤0 можно записать
(2.31)
Следует отметить, что полученное ранее уравнение Гриншилдса (2.18) является частным случаем уравнения Гринберга (2.31) при n=1.
Кинематические волны в транспортном потоке
Зависимость интенсивности движения от плотности называют фундаментальной диаграммой транспортного потока. Эта диаграмма может быть получена по модели Лайтхилла-Уизема при следующих допущениях:
транспортный поток непрерывен, его плотность k(x, t) есть число машин, занимающих единицу длины дороги;
интенсивность движения q(x, t) равна числу машин, пересекающих черту x за единицу времени, определяется локальной плотностью k:
(2.32)
скорость потока
т.е. средняя скорость является функцией
плотности
на участке дороги без съездов-въездов количество машин сохраняется.
Уравнения (2.32) и (2.19) образуют полную систему. После подстановки получим
(2.33)
где
–
скорость распространения возмущений.
Соотношение q(k)=kvв(k) играет важную роль в теории транспортных потоков и называется фундаментальной диаграммой (рис. 2.2). В модели Лайтхилла-Уизема эта зависимость непрерывна, следовательно, предельная пропускная способность участка дороги определяется плотностью потока.
Общий вид решения нелинейного уравнения (2.33):
(2.34)
где F –
произвольная функция. Соотношение
(2.34)
описывает бегущую волну, рассматриваемую
как волну
уплотнения в среде.
Волны типа (2.34) называют кинематическими
волнами,
что подчеркивает их кинематическое
происхождение в противоположность
динамической природе акустических и
упругих волн.
Ударные волны в транспортном потоке
Зависимость скорости от интенсивности движения приводит к образованию ударной волны в транспортном потоке. Пусть скорость автомобилей в потоке лежит в пределах
(2.35)
так, что [(dq)/(dv)] > 0. Если по какой-либо причине скорость некоторой части потока понизится на Δv, интенсивность движения понизится на kj(1– 2v/v0)Δv. Плотность этой части потока k повысится, и скорость будет далее снижаться. Возмущение скорости является незатухающим, что и демонстрирует неустойчивость поведения транспортного потока. В этих случаях автомобили в потоке вынуждены неоднократно трогаться с места и останавливаться. Это явление носит название ударной волны.
