CC BY
51
Моделирование многомерного случайного возмущения в задачах динамики рельсовых экипажей
А. Н. Савоськин,
д-р техн. наук, профессор кафедры «Электропоезда и локомотивы» Института транспортной техники и систем управления Московского государственного университета путей сообщения (ИТТСУМИИТ)
Ю. С.Ромен,
д-р техн. наук, профессор, главный научный сотрудник отделения «Перспективное развитие сети железных дорог и целевые программы» ОАО «ВНИИЖТ»
А. А. Акишин,
аспирант кафедры «Электропоезда и локомотивы» ИТТСУ МИИТ
Предложена методика выбора параметров генерации многомерного случайного процесса возмущения, которая обеспечивает реализацию заданного частотного диапазона, необходимого для моделирования колебаний исследуемых рельсовых экипажей при движении со скоростью от 10 до 150 м/с.
Генерация многомерного случайного процесса возмущения ||у (£)| на ЭВМ во временной области проводится по его известным характеристикам - матрице авто- и взаимных корреляционных функций (т)| или матрице спектральных и взаимных спектральных плотностей |Ф,Ш|(//) [2]. При этом не коррелирующие между собой временные последовательности x¡[nT\, получаемые с помощью типовых программ генераторов белого шума, необходимо преобразовать в многомерный гауссовский стационарный случайный процесс || т]ш (пАх - VпТ)| с заданной матрицей спектральных плотностей ЦСДЯЩ размером ixu = mxm, где i = 1, 2, 3, ..., m и u = 1, 2, 3, ..., m, с использованием матриц частотных характеристик, получаемых на основе решения задачи факторизации.
Для решения этой задачи предложен новый алгоритм, основанный на использовании временных характеристик формирующих фильтров. Процессы на выходе ГШ, дискретные, и вначале, как и в [2], необходимо получить матрицу дискретных спектральных плотностей |С„, (г) выходного процесса п[пТ], выполнив дискретное преобразование Лапласа
у(2)=г{х(лг)} = !>[«>-
' п=0
{Т. = еЧ=еР^ (1)
матрицы корреляционных функций
||с^2)|| = г{||^лГ)||).
(2)
Затем необходимо представить матрицу ЦС^Дг)! в виде:
1М2)НК(2)1х НМ41М2)||,
№1 -
(3)
где | гг1и [¿.ц - матрица дискретных частотных характеристик, эрмитовски сопряженная с матрицей частотных характеристик формирующих фильтров; ||СИ.Ц (г)|| - матрица единичных спектральных плотностей некоррелированных сигналов на выходе генераторов белого шума размером шхш:
100. .. 0
010 ..0
001 ..0
000 ..1
(4)
Задача решается во временной области, поэтому в отличие от [2] матрицу
(г)| необходимо преобразовать в матрицу импульсных характеристик формирующих фильтров с помощью обратного дискретного преобразования Лапласа:
-¿Ткм^«-
(5)
Зная входной процесс x¡[nT\ и импульсные характеристики, определяющие динамические свойства формирующего фильтра во временной области, можно найти выходной процесс п,„(пТ) с помощью интеграла свертки:
Реализация этого алгоритма может быть выполнена по новой схеме формирующего механизма (рис. 1), который
Рис. 1. Структурная схема формирующего механизма во временной области
содержит импульсные характеристики кии[пТ] и интеграторы, выполняющие вычисление интеграла свертки (6).
В качестве примера рассмотрим генерацию четырехмерного гауссовского стационарного случайного процесса возмущений в виде вертикальных и горизонтальных неровностей (х = V) левого и правого рельсов, вызывающего колебания рельсовых экипажей по аналитическому выражению корреляционной функции, приведенному в [1].
Выполним прямое дискретное преобразование Лапласа аналитического выражения для корреляционной функции в программном пакете МаНАВ с
+<ехр
[-<Лд*)2| [-"UA*)2]}
+2 jexp [-<sl (А*)2 ] j • [cos^„s,Ax] z +1
(7)
жение для элементов матрицы передаточных функций:
к
jexp [■(cos/? «•USIAX)Z+
1 jexp[-
+|ехр -«Мд*)2
+2 jexp [■(cos/?, a.„s,Ax)z + l
(8)
Выполнив обратное г-преобразо-вание для №ы(г) с помощью команды перейдем к импульсным характеристикам:
К* [пД*] = ^ХАм^К-Д*)2]* >[(«Иш1Д*)2|
(9)
ехр
или окончательно
= SiuslIXusl jexp
помощью команды <^гата» и получим выражение для дискретных спектральных плотностей в виде
к
ехр[-а^, (А*)2]} • [со80Иш,Дх>2 +
xcos[/3
kiusl
X (10)
Факторизовав (7) в программном пакете MatLab с помощью оператора «factor», найдем аналитическое выра-
Было выполнено сравнение графиков импульсных характеристик, построенных по этому выражению, и корреляционных функций со значениями параметров, приведенными в [1]. Для примера на рис. 2 приведено такое сравнение для автокорреляционных функций и импульсных характеристик горизонтальной неровности правого рельса и взаимной корреляционной функции, а также импульсной характеристики вертикальной неровности левого и горизонтальной неровности правого рельсов. Сравнение графиков R¡us|[nAx]
Rur, W -R.Tsravw1
Рис. 2. Графики импульсной характеристики (сплошные линии) и корреляционной функции (пунктир) горизонтальной неровности правого рельса (а); взаимной корреляционной функции и импульсной характеристики вертикальной неровности левого рельса и горизонтальной неровности правого рельса (б)
и кш[пАх\ свидетельствует о том, что эти процессы затухают практически синхронно. Это означает, что к„ш1[пАх\ достаточно адекватно описывает динамические свойства формирующего фильтра во временной области.
Выражение (10) для импульсной характеристики формирующего механизма соответствует дифференцируемому случайному процессу и может быть использовано не только при решении задачи генерации многомерного случайного процесса геометрических
неровностей пути, но и при решении аналогичных задач в различных приложениях.
По найденным импульсным характеристикам с использованием интеграла свертки была выполнена генерация многомерного случайного процесса во временной области по схеме формирующего механизма, приведенной на рис. 1. Результаты генерации одной из совокупностей реализаций четырехмерного случайного процесса геометрических неровностей пути (рис. 3), выполнен-
ные по рассмотренному алгоритму, показывают, что сгенерированные случайные процессы в среднем близки к исходным.
Для более подробного сравнения вычислены корреляционные функции и спектральные плотности для сгенерированных случайных процессов по [1], выполнено их сравнение с графиками заданных корреляционных функций и спектральных плотностей реального рельсового пути по выражениям, приведенным в [1].
Для примера на рис. 3 приведены некоторые графики для авто- и взаимных корреляционных функций, а на рис. 4, 5 - соответствующие им спектральные и взаимные спектральные плотности. Эти графики, построенные по сгенерированным реализациям, имеют удовлетворительную сходимость с заданными. Некоторое расхождение между ними объясняется особенностями функционирования генераторов белого шума, а также погрешностями спектрального анализа и не превышает расхождений, имеющих место при аппроксимации характеристик, полученных по экспериментальным данным [1].
Таким образом, решена задача по разработке алгоритма генерации многомерной геометрической неровности рельсового пути. Используя этот алгоритм, можно генерировать такую неровность при различных значениях скорости движения с любыми длиной реализации и шагом дискретизации. Это позволит численно решать во временной области многомерные задачи динамики подвижного состава с нелинейными характеристиками рессорного подвешивания.
Однако вследствие ограниченной длины реализации случайных возмущений, равной =1 км, максимум низкочастотной части спектра был ограничен высшей длиной волны неровности Ьв = = 50 м, а максимум высокочастотной составляющей соответствовал Ьн = 6,25 м. При решении задач динамики частотный диапазон возмущения должен быть в пределах от /н = 0,2 Гц до /в = 10 Гц.
Эти низшая и высшая частоты связаны с соответствующими длинами волн соотношениями
Рис. 3. Реализации четырехмерного случайного процесса геометрической неровности рельсов: 1 - сгенерированные, 2 - исходные; горизонтальные левого (а) и правого (б) рельсов; вертикальные левого (в) и правого (г) рельсов
где V - скорость движения, м/с.
(11) (12)
Л.тв.ч«2 Дт-.-чч2
МГ.....;......•......;.....1......;......:......у...........:......; £Г)10 ......;.........
Рис. 4. Автокорреляционные функции неровностей пути: левого рельса вертикальная (а) и правого рельса горизонтальная (б); взаимные корреляционные функции между правой вертикальной и горизонтальной неровностью рельсов (в); правой вертикальной и левой горизонтальной неровностью рельсов (г): здесь и на рис. 5, 6 сгенерированные - сплошные линии; по аналитическому выражению - пунктир
Рис. 5. Спектральные плотности неровностей: левого рельса вертикальных (а) и правого рельса горизонтальных (б)
В связи с этим при увеличении скорости движения частотный диапазон возмущения смещается в сторону высоких частот и может выйти за пределы диапазона частот колебаний системы. В этом случае возмущение практически перестанет возбуждать колебания и полученные результаты окажутся некорректными.
Определим требования к длине волн неровностей, которые должны быть в спектре возмущения, для решении задач динамики при значениях скорости
движения от 10 до 150 м/с, т. е. от 36 до 540 км/ч. Необходимую длину реализации определим по эмпирической формуле
5 + 10
Т =
(з+ю)/. •
(13)
(14)
лизации составит = ^^ = 50 с. Такой
величине { соответствует длина реализации (метры), равная
Ь = vt = 50у м.
р р
(15)
При принятом значении низшей частоты процесса /н = 0,2 Гц длина реа-
Для выделения заданного значения высшей частоты процесса /в = 10 Гц шаг дискретизации по времени примем
равным Т = ~гг= 0,01 с, что соответст-
' в
вует шагу по пути
Рис. 6. Составляющие взаимной спектральной плотности между правой вертикальной и правой горизонтальной неровностями рельсов: а - вещественная; б - мнимая, и составляющие взаимной спектральной плотности между правой вертикальной и левой горизонтальной неровностями рельсов: в - вещественная; г - мнимая
Дх = уТ = 0,01у м.
(16)
При заданных значениях низшей и высшей частот количество точек реализации составит
К к
ЛТ„= — = — = 5 000 р Т Ах 5
(17)
В соответствии с этим при скорости 150 м/с длина реализации процесса возмущения Ь , необходимая для корреляционно-спектрального анализа, должна
быть не менее Ь = 50*150 = 7500 м. За-
р
писи случайных процессов возмущения такой протяженности должны удовлетворять условиям стационарности, т. е. должны быть выполнены на площадке с однородной конструкцией пути при движении по прямому пути (без укло-
нов и кривых). Выполнить эти условия на пути такой протяженности практически невозможно, поэтому для получения достаточно продолжительной реализации восьмикратно повторим с соответствующими переходами использованную в [1] реализацию длиной Ь = 999 м, каждый раз изменяя ее направление, чтобы избежать совпадений. При получении этой реализации в [1] было сделано N = 5400 замеров неровностей с шагом Дх = 0,185 м. Таким образом, общая длина реализации составит Ьро6 = 8Ьр = 7992 м с шагом дискретизации Дх = 0,185 м и общим числом замеров Nо6 = 43 200. Спектральный анализ реализации четырехмерного случайного процесса неровности рельсов с такими параметрами, отнесенной к скорости движения у = 1 м/с, был выполнен по
методике, изложенной в [1]. Спектральная плотность (рис. 7) была аппроксимирована аналитическим выражением с восемнадцатью слагаемыми (к = 18). Некоторые параметры слагаемых приведены в таблице.
В составе спектральной плотности имеется целый ряд максимумов, приходящихся на длины волн неровностей от 6,5 до 750 м, что во много раз превышает диапазон длин волн 12,5-50 м, полученный в [1] при анализе исходных экспериментальных данных. Такой расширенный диапазон длин волн обеспечивает возможность моделирования колебаний рельсовых экипажей при скорости движения вплоть до 150 м/с, так как в соответствии с (11) низшей частоте процесса /н = 0,2 Гц соответствует длина волны 750 м.
Параметры аналитического выражения корреляционной функции горизонтальной неровности левого рельса
К 1, м в' в, м-1 а, м-1 а к 1, м в в, м-1 а, м-1 а
1 6,25 1610-2 32-10-3 234-10-4 10 350 28-10-4 810-5 1610-4
2 12,5 800-10-4 22-10-3 1953 10-4 11 400 25-10-4 14-10-5 14-10-4
3 25,0 400-10-4 1010-3 5078-10-4 12 450 22-10-4 1110-5 65-10-5
4 50,0 200-10-4 28-10-4 1823-10-4 13 500 20-10-4 25-10-3 39-10-5
5 100 10010-4 12-10-4 521-10-4 14 550 18-10-4 45-10-6 36-10-5
6 150 67-10-4 610-4 208-10-4 15 600 16-10-4 41-10-6 20-10-5
7 200 60-10-4 35-10-5 99-10-4 16 650 1510-4 61-10-6 41-10-5
8 250 40-10-4 45-10-5 26-10-4 17 700 14-10-4 21-10-6 2-10-5
9 300 33-10-4 3110-5 1810-4 18 750 13-10-4 2-10-6 1,4-10-5
Рис. 7. Графики спектральной плотности горизонтальной неровности левого рельса: 1 - экспериментальная; 2 - по аналитическому выражению
а) ||..ии
Рис. 8. Графики сгенерированных реализаций случайного процесса горизонтальной неровности левого рельса, построенные для скорости движения 20 м/с (а), 60 м/с (б), 120 м/с (в)
Рис. 9. График спектральных плотностей горизонтальной неровности левого рельса для решения задач динамики при скорости движения: 20 м/с (а), 60 м/с (б), 120 м/с (в)
Аналогичные процедуры были выполнены и для других неровностей, что дало возможность получить полную матрицу спектральных и взаимных спектральных плотностей, а также авто- и взаимных корреляционных функций.
Для примера были сгенерированы возмущения для скорости V = 20 м/с, 60 м/с и 120 м/с. Из графиков сгене-
рированных реализаций случайного процесса горизонтальной неровности рельса (рис. 8) следует, что с ростом скорости движения, как и следовало ожидать, увеличивается частота возмущения.
Графики спектральных плотностей возмущения для исследуемых значений скорости движения показывают, что вся энергия процессов приходит-
ся на выбранный диапазон частот 0,2-10 Гц (рис. 9). С ростом скорости вид спектральных плотностей изменяется: увеличивается доля высокочастотных составляющих. Если при скорости движения V = 20 м/с (рис. 9а) максимум Оф приходится на частоты 0,8-1,0 Гц, то при V =60 м/с (рис. 9б) - на диапазон 1,5-2,7 Гц, а при v = 120 м/с (рис. 9в) - на диапазон 2,2-4,8 Гц. При этом изменяются и ординаты максимумов. Кроме того, с ростом скорости движения уменьшаются отклонения спектральных плотностей сгенерированных и заданных случайных процессов.
Таким образом, разработанная методика генерации случайных процессов обеспечивает получение возмущений в необходимом диапазоне частот при заданных значениях скорости движения от 10 до 150 м/с. При этом реализацию возмущения необходимо генерировать заново для каждого исследуемого значения скорости, выбирая длину реализации и шаг дискретизации на основе предложенных соотношений.
Таким образом, разработанный алгоритм генерации многомерного случайного процесса на основе импульсных характеристик и интеграла свертки позволяет решать эту задачу на ЭВМ во временной области. Предложенная методика выбора параметров генерации многомерного случайного процесса возмущения обеспечивает реализацию заданного частотного диапазона, необходимого для моделирования колебаний исследуемых рельсовых экипажей. Приведенные примеры генерации многомерного случайного процесса возмущения при различных значениях скорости движения свидетельствуют, что с ростом скорости увеличивается высокочастотная составляющая процесса, а спектральный анализ сгенерированных реализаций возмущения соответствует заданному частотному диапазону колебаний системы от 0,2 до 10 Гц. и
Литература
1. Ромен Ю. С., Савоськин А. Н., Аки-шин А. А. Анализ случайных процессов геометрических неровностей рельсовых нитей // Изв. ПГУПС. 2014. № 1 (38). С. 22-32.
2. Матыаш И., Шилханек Я. Генератор случайных процессов с заданной матрицей спектральных плотностей // Автоматика и телемеханика. 1960. Т. 21, № 1. С. 29-34.
