МАТЕМАТИКА
DOI: 10.5862ZJPM.218.12 УДК 519.216.2
Т.А. Алексеева
Национальный исследовательский университет «Высшая школа экономики»,
Санкт-Петербургский филиал
математическое моделирование
случайных дорожных возмущений
при исследовании динамики мобильных машин
В работе представлена методика моделирования случайных дорожных возмущений на основе метода неканонического разложения случайных функций в виде детерминированных функций, зависящих только от трех случайных величин при любых законах распределения их вероятностей.
СЛУЧАЙНОЕ ВОЗМУЩЕНИЕ, НЕКАНОНИЧЕСКОЕ РАЗЛОЖЕНИЕ СЛУЧАЙНЫХ ФУНКЦИЙ, КОРРЕЛЯЦИОННАЯ ФУНКЦИЯ, СПЕКТРАЛЬНАЯ ПЛОТНОСТЬ, ИНТЕРПОЛЯЦИОННЫЙ МЕТОД.
Одной из основных задач при исследовании низкочастотных колебаний мобильных машин, в частности колесных тракторов, является учет случайных воздействий со стороны микрорельефа поля (дороги). Микрорельеф представляет собой поверхность поля или дороги со случайным расположением неровностей. Его принято рассматривать как случайную функцию, удовлетворяющую следующим требованиям: функция стационарна, ординаты микропрофиля подчиняются нормальному закону распределения, длины неровностей ограничены по верхнему и нижнему пределам, микропрофиль меняется случайным образом только в вертикальной продольной плоскости [8, 12 — 14].
Существуют различные методы представления случайных функций: с помощью обобщенных рядов Фурье, методами К. Ка-рунена, В.С. Пугачева, В.А. Котельникова и др. [1—5]. Эти методы представления позволяют построить для любой случайной функции бесконечный функциональный ряд, коэффициентами которого являются нормированные некоррелированные случайные величины. Однако все эти пред-
ставления случайной функции обладают рядом существенных недостатков. В частности, для получения приемлемой точности необходимо использовать достаточно большое число членов ряда. Кроме того, линейная форма этих представлений мало пригодна при исследовании нелинейных колебательных систем.
В данной работе в основу методики моделирования случайных дорожных возмущений положено неканоническое представление случайной функции, предложенное В.И. Чернецким [6]. Оно отличается от других тем, что дает точное совпадение случайной функции и ее представления в рамках корреляционной теории при использовании для стационарной случайной функции всего трех случайных величин. Такое представление случайной функции имеет широкий спектр применения, в частности используется для анализа электрических цепей [7], а также представления реализаций случайных процессов сложной структуры [8—11].
Согласно принятым допущениям, достаточными статистическими характеристиками случайной функции дорожного
4
Математика
возмущения являются корреляционная функция и спектральная плотность. В результате анализа применяемых на практике выражений для аппроксимации корреляционных функций была выбрана функция следующего вида [12, 13]:
(
R (т) = D
Л0в а°(1 + а0 |т|) +
I Лв
а.
(1)
cos В.т + —L sin В. т p. 11
/у
Соответствующая ей спектральная плотность имеет вид
Sh (ш) = Dh
С 4 Ло а0
(ш2 +а0)2
(2)
I
44а,и о.
1(ш2 -ш),) + 4а2га2
где А0, ^ — безразмерные коэффициенты, учитывающие величину спектральной плотности микрорельефа в точках максимума; Бк, см2 — дисперсия возмущения; т, с — время; N — число максимумов спектральной плотности; ю, с1 — частота.
Корреляционные коэффициенты а0, а;, р ; связаны со скоростью движения машины V соотношениями:
а„ = а V; а. = а. V; В. = В. V, (3)
0 0 ед. " I гед. " 'I <гед. "
где а0ед., а,.ед., р ,.ед. - коэффициенты, соответствующие единичной скорости.
Такой выбор корреляционной функции обусловлен тем, что, во-первых, она позволяет учитывать наличие периодической составляющей в случайной функции возмущения, во-вторых, на основе экспериментальных исследований [12, 13] были получены значения коэффициентов в выражении (1) для ряда основных типов микропрофилей. Это явилось основой для построения реализаций случайных функций возмущения в аналитической форме.
Из-за сложности проблемы построения реализаций случайных функций возмущений по их статистическим характеристикам, обычно в выражении для корреляционной функции (1) используют одно или два слагаемых.
В данной работе решена задача построения реализаций случайных дорожных возмущений при удержании в указанном выражении пяти и более членов, что позволяет учесть максимумы спектральных плотностей и получить более достоверную информацию о спектральном составе микрорельефа..
Эта методика используется для построения реализаций случайных функций возмущений ) и /¡2 (/), входящих в систему дифференциальных уравнений, моделирующих колебания мобильных машин. Эти функции имеют вид:
(5)
hí(t) = m^ + ^ sin ш t + A,2 cos ш t; (4)
h^(t) = m^ + ^ sin ш (t - t0) + + X2 cos ш (t - t0), где t0 = l / v (l — база мобильной машины, v — скорость движения); mh, m^ — математические ожидания случайных функций; Х1, ю — независимые случайные величины.
Необходимые для построения реализаций вероятностные характеристики случайных функций hí(t) и /¡2(tt заданы: их математические ожидания равны нулю, а корреляционная функция определяется выражением (1).
Законы распределения вероятностей случайных величин и Х2 произвольны, а плотность распределения случайной величины ю определяется по корреляционной функции (1):
f (ш) =
1 r R/ (т) -
2п D,
d т =
2 Л а0
п(а0 + ш )
N
+ I
2Л,. а,. шо,.
(6)
tí п(шо,. - ш ) + (2а,.ш)
2
При интерполяционном методе анализа необходимая выборка значений случайных величин Х1, ю рассчитывается при помощи соответствующего ^-преобразования и специальных таблиц стандартных узлов типа Чебышева для любых законов распределения вероятностей случайных величин.
Для случайных величин и А,2 были использованы таблицы узлов типа чебы-шева нормального распределения. Рабочие
узлы Х1к[ и Х2кг для всевозможных индексов кг =1,2, ..., = 1,2, ..., рассчитываются по формулам:
ч = + ^ хк X ^ = + стх2 х*2;
^ = = л/А; М (X,) = М (А 2) = 0,
где х^, х^ — стандартные узлы Чебышева для нормального распределения.
Для случайной величины ю по указанной плотности распределения вероятностей (6) построено ^-преобразование к равномерному закону распределения:
х +1
2A а0
N +i
_»1п(а2+®2)2 2 Д-а,.ю0,-
(7)
1 л(ш°; - ш ) + (2а (ш)
d ш.
Полученное уравнение не дает возможности выразить ю в явной форме в зависимости от х. Поэтому решение этого уравнения было сведено к численному интегрированию дифференциального уравнения
d ш dx
= (1/2)
2 Ao а;
N + 1
п(а2 + ш2)2
24а,-ш0,- Л
(8)
,=1 тс(ш°,. - ш ) + (2а ,.ш)
при начальных условиях х = 0, ю = 0.
В результате получены рабочие узлы соответствующие стандартным узлам
шк х.
(к3 = 1, 2, ..., q 3).
Аналогично строится ^-преобразование для случайной величины ю с плотностью распределения (6); при этом используются таблицы стандартных узлов типа чебышева для нормального и экспоненциального законов распределения.
При расчетах для случайных величин Х1
и Х2 выбрано по два узла, а для случайной величины ю — пять узлов интерполирования, т. е. д1 =2, #2 = 2, #3 =5; к1 = 1, 2; к2 = 1, 2; к3 = 1, 2, 3, 4, 5. В этом случае необходимо рассчитать по 20 значений X^, Хк2, н^, п = 4424з = 2 • 2 • 5 = 20.
Такой выбор числа узлов интерполирования обеспечивает сочетание требуемой точности и минимального объема вычислений в рамках интерполяционного метода.
При программном синтезе реализаций случайных дорожных возмущений по заданной корреляционной функции (1) корреляционные коэффициенты а и р; варьировались в зависимости от скорости движения машины в соответствии с выражением (3).
По разработанной методике были рассчитаны таблицы рабочих узлов типа чебы-
шева X
X
'V
ш
кз
случайных величин Х1, Х2, ю для восьми основных фонов дорог при различных скоростях движения, что является исходным материалом для исследования динамикимобильныхмашин.
Итан, разработанная методика математического моделирования случайных дорожных возмущений дает точное представление случайной функции возмущения в рамках корреляционной теории при использовании всего трех случайных величин, позволяет получить более достоверную информацию о спектральном составе микрорельефа, упростить исследование динамики мобильных машин и сократить объем вычислений по сравнению с другими методами. Компьютерная модель не критична к программному средству и может быть реализована в средах, допускающих стохастическое моделирование.
список литературы
[1] Пугачев В.С. Теория случайных функций и ее применение к задачам автоматического управления. М.: Физматгиз, 3-е изд., 1962. 883 с.
[2] Karhunen K. Uber lineare Methoden in der Wahrscheinlichkeitsrechnung // Ann. Acad. Sci. Fennicae. Ser. A. I. Math.-Phys. 1947. No. 37. Pp. 1-79.
[3] Свешников А.А. Прикладные методы теории случайных функций. М.: Наука, 1968. 463 с.
[4] Фурунжиев Р.и. Проектирование оптимальных виброзащитных систем. Монография. Минск: Вышейшая школа, 1971. 320 с.
[5] Ламперти Дж. Случайные процессы. Обзор математической теории. Киев: Вища школа, 1983. 227 с.
[6] Чернецкий В.и. Анализ точности нелинейных систем управления. М.: Машиностроение, 1968. 246 с.
со
2
Математика -►
[7] Мищенко т.Н., Михаличенко П.Е., Костин Н.А. Вероятностные характеристики случайной функции напряжения на токоприемнике первого украинского электровоза ДЭ-1 // Научные журналы НТУ «ХПИ». Электротехника и электромеханика. 2003. № 2. С. 43-46.
[8] Фурунжиев Р.и. Автоматизированное проектирование колебательных систем. Монография. Минск: Вышейшая школа, 1977. 452 с.
[9] Чернецкий В.и. Математическое моделирование стохастических систем. Петрозаводск: Петрозаводский гос. ун-т, 1994. 486 с.
[10] Перов С.Н., Скворцов Ю.В. Представление случайных процессов с помощью неканонического разложения // Вестник Самарского государственного аэрокосмического ун-та. 2008.
№ 1 (14). С. 226-235.
[11] Павлюк Ю.С., Сакулин В.Д. Нелинейная динамика механических систем при случайном воздействии // Динамика, прочность и износостойкость машин. 1995. № 1. С. 33-43.
[12] Росляков В.П. Аппроксимация корреляционных функций случайных процессов в задачах динамики сельхозмашин // Механизация и электрификация социалистического сельского хозяйства. 1969. № 8. С. 15-19.
[13] Росляков В.П. Статистические характеристики рельефов обрабатываемых полей // Труды Курского СХИ. Воронеж. 1969. Вып. 3.
[14] Хачатуров А.А. Динамика системы дорога-шина-автомобиль-водитель. М.: Машиностроение, 1976. 535 с.
СВЕДЕНИЯ ОБ АВТОРЕ
АЛЕКСЕЕВА татьяна Анатольевна - кандидат физико-математических наук, доцент департамента прикладной математики и бизнес-информатики Санкт-Петербургского филиала Национального исследовательского университета «Высшая школа экономики».
190008, Российская Федерация, г. Санкт-Петербург, ул. Союза Печатников, 16 tatyanalexeeva@gmail.com
Alexeeva T.A. SIMULATION OF OCCASIONAL ROAD DISTURBANCES FOR MOBILE VEHICLES DYNAMIC STUDIES.
This paper focuses on a simulation procedure of occasional road disturbances reasoning from non-canonical occasional functions analysis method in terms of determinate functions depending on only three casual quantities with any distribution law of their probabilities. A developed mathematical method for simulation of random road disturbances gives an accurate representation of the random function perturbations in the framework of the correlation theory using only three random variables. The proposed approach allows obtaining more reliable information about the spectral composition of the microrelief, simplifying studies of the mobile vehicles dynamics and reducing the computation amount as compared to other methods.
OCCASIONAL DISTURBANCE, NON-CANONICAL ANALYSIS OF OCCASIONAL FUNCTIONS, CORRELATIONAL FUNCTION, SPECTRAL DENSITY, INTERPOLATION METHOD.
REFERENCES
[1] V.S. Pugachev, Teoriya sluchajnykh funk-tsij i ee primenenie k zadacham avtomaticheskogo upravleniya [The theory of random functions and its application to problems of automatic control], Fizmatgiz, Moscow, 3rd ed., 1962.
[2] K. Karhunen, Uber lineare Methoden in der Wahrscheinlichkeitsrechnung, Ann. Acad. Sci. Fen-nicae, Ser. A. I. Math.-Phys. 37 (1947) 1-79.
[3] A.A. Sveshnikov, Prikladnye metody teorii sluchajnykh funktsij [Applied methods of the theory of random functions], Nauka, Moscow, 1968.
[4] R.I. Furunzhiev, Proektirovanie optimal'nykh vibrozashchitnykh system [Designing optimal vibration isolation systems], Vyshejshaya shkola, Minsk, 1971.
[5] J. Lamperti, Sluchajnye protsessy. Obzor matematicheskoj teorii [Random processes. Over-
view of the mathematical theory], Vishcha shkola, Kiev, 1983.
[6] V.I. Chernetskij, Analiz tochnosti neline-jnykh sistem upravleniya [Accuracy analysis of nonlinear control systems], Mashinostroenie, Moscow, 1968.
[7] T.N. Mishchenko, P.E. Mikhalichenko, N.A.
Kostin, Veroyatnostnye kharakteristiki sluchajnoj funktsii napryazheniya na tokopriemnike pervogo ukrainskogo elektrovoza DE-1 [Probabilistic characteristics of a random function of voltage in the current collector of the first Ukrainian electric locomotive DE-1], Nauchnye zhurnaly NTU «KhPI». Elektrotekhnika i elektromekhanika. 2 (2003) 43-46.
[8] R.I. Furunzhiev, Avtomatizirovannoe proek-tirovanie kolebatel'nykh sistem, [Computer-aided
design of vibrating systems], Minsk: Vyshejshaya shkola, 1977.
[9] V.I. Chernetskij, Matematicheskoe mod-elirovanie stokhasticheskikh system [Mathematical modeling of stochastic systems], Petrozavodskij gos. un-t, Petrozavodsk, 1994.
[10] S.N. Perov, Yu.V. Skvortsov, Predstavle-nie sluchajnykh protsessov s pomoshch'yu nekanon-icheskogo razlozheniya [Representation of stochastic processes using non-canonical decomposition], Vestnik Samarskogo gosudarstvennogo aerokosmi-cheskogo un-ta. 1 (14) (2008) 226-235.
[11] Yu.S. Pavlyuk, V.D. Sakulin, Nelinejnaya dinamika mekhanicheskikh sistem pri sluchajnom vozdejstvii [Nonlinear dynamics of mechanical systems at accidental exposure], Dinamika, prochnost'
i iznosostojkost' mashin. 1 (1995) 33—43.
[12] V.P. Roslyakov, Approksimatsiya kor-relyatsionnykh funktsij sluchajnykh protsessov v zadachakh dinamiki sel'khozmashin [Approximation of correlation functions of random processes in dynamics problems of agricultural machinery], Mekhanizatsiya i elektrifikatsiya sotsialisticheskogo sel'skogo khozyajstva. 8 (1969) 15-19.
[13] V.P. Roslyakov, Statisticheskie kharakteris-tiki rel'efov obrabatyvaemykh polej [Statistical characteristics of the reliefs of cultivated fields], Trudy Kurskogo SKhl, Voronezh. 3 (1969).
[14] A.A. Khachaturov, Dinamika sistemy dor-oga-shina-avtomobil'-voditel' [Dynamics of a road-tire-car-driver system], Mashinostroenie, Moscow, 1976.
THE AuTHOR
ALEXEEVA Tatiana A.
National Research University Higher School of Economics
16 Soyuza Pechatnikov St., St. Petersburg, 190008, Russian Federation
tatyanalexeeva@gmail.com
© Санкт-Петербургский политехнический университет Петра Великого, 2015