CC BY
20
УДК 629.4.015:51.001.57
С. В. МЯМЛИН, А. Н. ПШИНЬКО, Е. А. ПИСЬМЕННЫЙ (ДИИТ), В. И. ПРИХОДЬКО (ОАО «Крюковский вагоностроительный завод»)
ОБЪЕКТНО-ОРИЕНТИРОВАННАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОСТРАНСТВЕННЫХ КОЛЕБАНИЙ РЕЛЬСОВЫХ ЭКИПАЖЕЙ
Докладно розглянуто методику складання диференщальних pÎBHHHb коливань рейкового ек1пажа, що до-зволяе визначити принципи «автоматично!» побудови цих piвнянь. Описано класи, що мiстять функцiï, якi виконують усi необхiднi опеpацiï над самими об'ектами i зв'язками, що ïx з'еднують, i з яких може бути по-будована модель ешпажа.
Подробно рассмотрена методика составления дифференциальных уравнений колебания рельсового экипажа, что позволило определить принципы «автоматического» построения этих уравнений. Описаны классы, которые содержат функции, выполняющие все необходимые операции над самими объектами и связями, их соединяющими, из которых может быть построена модель экипажа.
In article the technique of drawing up of the differential equations of fluctuation of rail crew that has allowed to define principles of «automatic» construction of these equations is in detail considered. Classes which contain the functions which are carrying out all necessary operations above objects and communications, their connecting from which the model of crew can be constructed are described.
В настоящее время существует множество математических моделей, описывающих пространственные колебания и динамическую нагруженность различных типов подвижного состава [1-5], однако большинство из них предназначено для моделирования какого-то одного типа подвижного состава. В настоящей работе предложен метод, позволяющий «собирать» модель любого типа подвижного состава. Рельсовый экипаж органично разделяется на объекты, в качестве которых могут быть выбраны элементы его конструкции. Например, для полувагона объектами динамической системы могут быть колесные пары, боковые рамы тележек, надрессорные балки и кузов, а для вагона с двойным рессорным подвешиванием - те же тела, только вместо боковых рам объектами являются рамы тележек.
Вначале рассмотрим методику составления уравнений пространственных колебаний объекта, на который действуют силы со стороны других объектов [6-9].
На рис. 1 показан объект и силы, действующие на него со стороны двух других объектов. Составим выражения для перемещений точек А и В объекта:
• для точки А
ZA = LZ - LX Ф + LY6,
Ya = ly + LX V - LZ6,
Xa = lx + lz Ф - Ly V ;
(1)
для точки B
ZB = Lz - Lx Ф - ly 6 > yb = ly + Lx V - Lz6, XB = Lx + Lz Ф + ly V,
(2)
где X А , ХВ, УА , УВ , XА , ХВ - координаты точек А и В по осям X, У, X соответственно; ф, 9, у - углы поворота объекта вокруг осей X, У и X соответственно.
Как видно из уравнений (1) и (2), формулы определения перемещений во всех плоскостях для точек А и В отличаются только знаками перед Ьу 9 и Ьу у, так как проекции этих точек на ось У находятся по разные стороны от центра масс объекта. Для других точек выражения перемещений запишутся аналогично с учетом знаков перед ЬХ, Ьу и . Например, уравнения для точки В будут иметь вид
ZD = LZ
lx Ф - ly 6,
yd = ly - Lx V + Lz6 > XD = LX - Lz Ф + LY V.
(3)
Таким образом, если под величинами Lx, LY и Lz понимать не абсолютные значения координаты точки на оси X, Y и z, а учитывать знаки координаты, то для произвольной точки объек-
та уравнение перемещений можно записать в общем виде (4), где у величин Ьх, ЬУ и Ь2 необходимо учесть знаки.
2 = Ь2 - Ьх ф + ЬУ 0, У = Ьу + Ьх у - ь2 0, X = Ьх + Ь2 ф - Ьу у.
Бхь
Рассмотрим методику составления уравнений динамики для объекта, изображенного на рис. 1. При этом будем считать, что силы
( 52л , 5 2в , 52э , 5ул , 5ув , 5уо , 5 хл , 5 хв , 5хо ),
действующие на объект со стороны объекта, который расположен «ниже» по оси 2 («правее» по оси У или «дальше» по оси X), положительные, а силы, действующие со стороны объекта, который расположен «выше» по оси 2 («левее» по оси У или «ближе» по оси X), отрицательные. Тогда
Рис. 1. Объект динамической системы
^X0 + 52Л (ЬУЛ - Ь2Л0) + 52В (ЬУВ - Ь2В0) + +(ЬУО - Ь2О0) - 5УЛ (Ь2Л + ЬУЛ0) -
-5УВ (Ь2В + ЬУВ0) - 5УО (Ь2О + ЬУО0) = 0,
где т, Ж - соответственно масса и вес объекта; 2, У , X - координаты центра масс объекта; Зу, Ь2, Jх - моменты инерции объекта
+ Ь7Аф)
т2 тУ тХ
5
2А
^В - $2О = Ж,
8уВ = О,
- 5x0 = 0.
5УЛ +
5ХЛ + 5ХВ
^ф - 52А (ЬХА + Ь2Афф - 52В (ЬХВ + Ь2Вфф --52О (ЬХО + Ь2О ф) + 5ХЛ (Ь2Л - ЬХАф) + + 5ХВ (Ь2В - ЬХВф) + 5ХО (Ь2О - ЬХОф) = 0, J2 у + 5УЛ (ЬХЛ 5УВ (ЬХВ -ЬУВ у) +
+5УО (ЬХО - ЬУОу) - 5ХЛ (ЬУА + ЬХАу) --5ХВ (ЬУВ + ЬХВу) - 5ХО (ЬУО + Ьхоу) =
О, (5)
относительно соответствующих осей. В формуле (5) выражения типа Б2А
представляют собой моменты соответствующих сил, а сумма в скобках - изменяющиеся в зависимости от угла поворота плечи сил.
Для использования полученных выше уравнений динамики для составления программы вычислений их целесообразно переписать в следующем виде:
т2 = Ж + ^ 52 тУ = £ 5У,
тХ = Е 5х, J2 V = Е М2,
^ ф = Х МУ,
Jх 00 = Е Мх.
Теперь в правых частях уравнений стоят суммы сил или суммы моментов сил. Для решения полученных уравнений динамики необходимо определить основные параметры объектов динамической системы. В качестве объектов рассматриваются кузов, элементы тележек, колесные пары. Перечень объектов и обозначения масс даны в табл. 1.
Таблица 1
Обозначения масс объектов рельсового экипажа
Объект модели Условное обозначение массы
Экипаж тэ
Кузов т
Тележка тт
Рама тележки тр
Колесная пара тк
S=кЛ + рЛ ,
(7)
Для моделирования пространственных колебаний реальных рельсовых экипажей необходимо использовать и другие возможные связи между объектами [6; 7; 9]:
1. Билинейная упруговязкая связь (рис. 2).
h/
Помимо этого объектами модели «экипаж-путь», общими для всех типов рельсовых экипажей, являются приведенные массы пути - в вертикальной тр и горизонтальной тр плоскостях.
Для удобства расчетов кроме перечисленных основных параметров объектов введены некоторые дополнительные. Это, прежде всего, параметры, необходимые для расчетов движения экипажа на криволинейных участках пути, название объекта, указатель на экипаж, к которому относится данный объект, и некоторые другие.
Параметры объектов рельсовых экипажей приводятся в технической документации.
Теперь для составления уравнений динамики группы объектов необходимы уравнения связи, которые определят зависимость величин сил в связях от величины относительных перемещений объектов.
Для линейной упруговязкой связи эта зависимость определяется выражением
Рис. 2. Характеристика билинейной упруговязкой связи:
А,, и Аг,8г - координаты точек, в которых происходит изменение жесткости связи; к1, к2, к3 - жесткость связи на соответствующих участках
Усилие в связи вычисляется по формуле (8)
Бы = + к1 (А-А,) + рА при А< А,, ^ = + к2 (А-А, ) + рА при А, < А < А г,[ (8) 8Ы = 8Г + к3(А-Аг) + рА при А> Аг.
Задавая различные значения параметров А,, 81, Аг, 8 и к1, к2, к3, можно получить практически любую желаемую характеристику связи. 2. Связь типа «сухое трение» (рис. 3).
h j T j
i- 1 - T
где к - жесткость связи в одном из трех направлений; А = £ - £ (- перемещение в одном из трех линейных или угловых направлений «нижнего» объекта, £>2 - перемещение в одном из трех линейных или угловых направлений «верхнего» объекта); Р - вязкость связи в одном
из трех направлений; А = - £ (£1, £2 - скорости перемещений соответствующих объектов).
Рис. 3. Характеристика связи типа «сухое трение»: T - сила сухого трения; к1, к2 - жесткости соответствующих участков характеристики связи
Усилие в связи вычисляется по формуле (9)
T = fN,
Sd = к2Л + рЛ + T sign (Л) при Л1 < Л < Л Sd = Sr + к (Л - Лr) при Л > Лr, Sd = Si + к (Л - Лг) при Л<Лг,
(9)
где / - коэффициент трения; N - сила, действующая по нормали к поверхности трения.
3. Связь между колесом и рельсом (крип) (рис. 4).
Рис. 4. Характеристика связи между колесом и рельсом (крип):
8 - относительная скорость проскальзывания колеса по рельсу; - полная сила крипа
Проскальзывания и силы крипа, согласно теории Картера, определяются следующими выражениями [6; 10; 11]:
-
b\i/ Ar — + —
V r
s„ =
V (кр - гё )'
s = Л'8х + sy •
f = 235P - 2,4P2 + 0,01P3,
h =
= f
F =
= f
1 + hV'
(10)
Sk = -F s,
Skx = Fsx,
Sky = Fsy ,
Таким образом, подробное рассмотрение методики составления дифференциальных уравнений колебания рельсового экипажа позволило определить принципы «автоматического» построения этих уравнений. Разработаны классы для описания объектов и связей рельсового экипажа. Описанные выше классы содержат функции, выполняющие все необходимые операции над самими объектами и связями, их соединяющими. То есть классы, описывающие объекты и связи рельсового экипажа как механической системы, являются теми «кирпичиками», из которых может быть построена («собрана») модель экипажа.
1.
2.
где Ь - половина расстояния между точками контакта колес одной колесной пары с рельсами; у -виляние колесной пары; - поперечное смещение колесной пары относительно рельса; Лг - отклонение от среднего значения радиуса колеса; 0 - боковая качка колесной пары; /т - коэффициент трения колеса о рельс; 5кх, 5ку - силы крипа в горизонтальном продольном и горизонтальном поперечном направлениях соответственно.
Перечисленные выше типы динамических связей используются в рельсовых экипажах. Анализ выражений (7)-(10) и рис. 1 позволяют определить основные параметры связей. Во-первых, это направления связи (вертикальная, горизонтальная продольная и горизонтальная поперечная), во-вторых, объекты, к которым приложена каждая связь, в-третьих, координаты точек приложения связи к объектам (на рис. 1 обозначены Ьх, ЬУ , Ь2), в-четвертых, деформация и скорость деформации связи - для расчета сил в связи. Полный текст программы вычислений характеристик связей между телами рельсовых экипажей приведен в [6].
3.
4.
5.
6.
7.
9.
10.
БИБЛИОГРАФИЧЕСКИМ СПИСОК
Блохин Е. П. Математическая модель пространственных колебаний четырехосного рельсового экипажа / Е. П. Блохин, В. Д. Данович, Н. И. Морозов; Днепропетровский институт инженеров железнодорожного транспорта. - Д., 1986. - 39 с. - Рус.-Деп. в ЦНИИТЭИ МПС 29.09.86, № 7252 ж. д. Коротенко М. Л. Дифференциальные уравнения пространственных колебаний четырехосного вагона с учетом конечной жесткости кузова и инерционных свойств основания / М. Л. Коро-тенко, В. Д. Данович // Проблемы механики наземного транспорта: Межвуз. сб. науч. тр. - Д.: ДИИТ, 1973. - Вып. 199/25. - С. 3-13. Математическое моделирование колебаний рельсовых транспортных средств / Под ред. В. Ф. Ушкалова. - К.: Наук. думка, 1989. - 240 с. М. Ф. Вериго, Взаимодействие пути и подвижного состава / М. Ф. Вериго, А. Я. Коган - М.: Транспорт, 1986. - 560 с.
Гарг В. К. Динамика подвижного состава / В. К. Гарг, Р. В. Дуккипати. - М.: Транспорт, 1988. - 391 с.
Мямлин С. В. Моделирование динамики рельсовых экипажей. - Д.: Новая идеология, 2002. - 240 с. Мямлин С. В. Методика моделирования пространственных колебаний железнодорожного экипажа // Залiзничний транспорт Украни. - 2001. -№ 2. - С. 2-5.
Pshin'ko O. Simulation of vehicle dynamic loading by means of object-oriented programming / O. Pshin'ko, Y. Blokhin, S. Myamlin // Archives of transport. Polish academy of sciences. Committee of transport. - 2002. -Vol.14, issue 1. - Р. 67-75. Blokhin Y. P. Mathematical modeling of body braces in a railway vehicle / Y. P. Blokhin, S. V. Myamlin, O. M. Pshin'ko // Problemy eksploatacji: Kwartalnik. -Radom: Wyd. Instytutu Technologii Eksploatacji. -2002. - № 1. - P. 269-276.
Carter F. W. On the Action of a Locomotive Driving Wheel // Proc. Royal Soc. - Vol. 112, ser. A. -1926. - P. 151-157. 11. Carter F. W. On the Stability of Running of Locomotives // Proc. Royal Soc. - Vol. 121, ser. A. -1928. - P. 585-611.
Поступила в редколлегию 25.02.2005.
