Методы моделирования динамики железнодорожных колесных пар с учетом упругости

УДК 531.39:51-74

DOI: 10.30987/article 5cb58f50b38371.23941436

Г.В. Михеев, Д.Ю. Погорелов, АН. Родиков

МЕТОДЫ МОДЕЛИРОВАНИЯ ДИНАМИКИ ЖЕЛЕЗНОДОРОЖНЫХ КОЛЕСНЫХ ПАР С УЧЕТОМ УПРУГОСТИ

Предложены два подхода к моделированию динамики железнодорожных колесных пар с учетом упругости, в рамках которых рассмотрены ко-нечноэлементные модели с вращающейся и невра-щающейся сеткой. Уравнения движения колесных пар выведены с применением присоединенной системы координат и результатов модального анализа. Кинематика профиля колеса описана с учетом уп-

ругих перемещений узлов. Оба подхода реализованы в программе «Универсальный механизм». Приведены первые результаты моделирования, подтверждающие корректность предложенных методов.

Ключевые слова: динамика систем тел, железнодорожные колесные пары, упругость, моделирование.

G.V. Mikheev, D.Yu. Pogorelov, A.N. Rodikov

METHODS FOR WHEEL PAIR DYNAMICS MODELING TAKING INTO ACCOUNT ELASTICITY

The computer simulation of railway wheel pair (WP) dynamics taking into account elasticity allows analyzing WPs high-frequency oscillations and vi-broaccelerations of crew part units connected with them, modeling strain gage wheel pairs, estimating the WP deformation impact upon force distributions in gears of locomotive drives, and also solving other urgent problems.

In the paper there are offered two approaches to the analysis of the dynamics of elastic WPs within the limits of which their finite element models with the rotating and nonrotating grid are under consideration. The WP kinematics is presented as a sum of its motion as an absolute solid together with the jointed co-

Введение

Программы моделирования динамики систем тел предоставляют эффективный инструмент анализа механических систем, которые могут быть представлены абсолютно твердыми или упругими телами, соединенными шарнирами или силовыми элементами. Одной из основных задач, которая должна быть решена для применения данного подхода к исследованию динамики железнодорожных экипажей, является разработка моделей контакта между колесом и рельсом. Моделирование колесных пар абсолютно твердыми телами, взаимодействующими с невесомым рельсом, является вполне приемлемым подходом при анализе динамики экипажей в частотном диапазоне до 30 Герц. Однако многие исследования невозможны без уче-

ordinates (C) and small elastic removals regarding coordinates given which are calculated on the basis of the results of WP modal analysis. A wheel profile kinematics is described taking into account unit elastic motions. There are presented algorithms for the calculation of generalized forces of moving loads in the wheel-rail contact taking into account WP elasticity.

Both approaches are realized in the "Universal mechanism" program. The first results of modeling confirming the correctness of the methods offered are shown.

Key words: dynamics of body system, railway wheel pairs, elasticity, modeling.

та упругости колесных пар и без подробного моделирования пути с учетом его инерционных и упругих свойств. Примерами таких исследований являются анализ высокочастотных (вплоть до звуковых) вибраций колесных пар и рельсов, расчет упругодеформированного состояния и долговечности колесных пар, исследования износа колесных пар и рельсов и т.д.

В работе [1] представлена модель колесной пары, разработанная Морисом (МогуБ), для изучения процессов, приводящих к отклонениям от круглой формы колес высокоскоростного поезда 1СЕ-1. Она включает несколько абсолютно твердых тел, соединенных пружинами с изгиб-ной и крутильной жесткостью и демпферами. Основная проблема в применении

подобных подходов - идентификация параметров. Только несколько низших частот и соответствующих собственных форм могут учитываться корректно. В работе [2] представлена двухмерная модель, разработанная Штольцем (Б2о1с), для изучения совместных вибраций колесной пары и рельсов. Она включает ось, которая моделируется как сплошной стержень, имеющий изгибную и крутильную жесткость, но абсолютно твердый в продольном (осевом) направлении. Колеса и дисковые тормоза представляются жесткими кольцами, соединенными с осью невесомыми упругими изотропными мембранами.

В настоящее время модальный подход является наиболее распространенным и универсальным методом моделирования динамики упругих колесных пар. Модели создаются методом конечных элементов (МКЭ). Кинематика колесной пары представляется суммой её движения как абсолютно твердого тела вместе с присоединенной системой координат (СК) и малых упругих перемещений относительно данной СК. Подобный подход распространен в динамике систем тел в качестве общей методики моделирования упругих подсистем; он описан, например, в работе [11].

Основной трудностью его применения к моделированию динамики колесной пары является расчет подвижных контактных сил и обобщенных сил от них. Для решения данной проблемы предпринимаются попытки вывода уравнений движения с невращающейся конечноэлементной сеткой. Подобный подход обобщенно называется подходом Эйлера - в противоположность подходу Лагранжа, традиционно применяемому в динамике систем тел. Говоря простым языком, подход Лагранжа заключается в том, что исследователь наблюдает за точками объекта, движущегося в пространстве, а при подходе Эйлера внимание обращено на определенную область пространства, занятую сплошной средой. Основное преимущество подхода (координат) Эйлера заключается в малых смещениях точек контакта между колесом и рельсом относительно локальной системы координат колесной пары. Это позволяет прикладывать контактные силы к не-

скольким заранее выбранным узлам, или даже к одному узлу, и тем самым повысить эффективность моделирования.

Современные методы вывода уравнений движения упругой колесной пары, использующие подход Эйлера, предлагаются в недавних публикациях [3-9]. При этом принимают некоторые упрощения ввиду сложности получаемых выражений. В работах [3-6] уравнения описывают прямолинейное движение с постоянной угловой скоростью в инерциальной системе координат.

В работе [10] уравнения выводятся в неинерциальной системе отсчета, движущейся вместе с экипажем. Динамика упругой колесной пары изучается отдельно от экипажа в предположении, что её упругость не влияет на перемещения тележки. Для вычисления сил взаимодействия между колесной парой и тележкой выполняется предварительное моделирование в программе анализа динамики систем тел.

В данной работе рассматриваются методы моделирования упругих колесных пар, разработанные авторами и реализованные в программном комплексе «Универсальный механизм» (ПК «УМ») (www.universalmechani sm .com). Кинематика любой упругой подсистемы описывается в «УМ» с применением присоединенной системы координат и метода суперэлементов. Модальный анализ упругой подсистемы выполняется во внешних программах конечноэлементного анализа, таких как ANSYS, MSC.NASTRAN, NX NASTRAN, ABAQUS. Для вывода уравнений движения применяется подход Лагранжа. Для применения его к моделированию вращающегося железнодорожного колеса разработаны следующие алгоритмы: 1) расчет кинематики профиля колеса с использованием положения и скорости узлов обода, расположенных вблизи контактной площадки; 2) расчет обобщенных сил от сил взаимодействия колеса с рельсом в произвольных точках поверхности качения. Полученные уравнения включают выражения сил инерции без каких-либо упрощений. Таким образом, возможно моделирование динамики экипажа в кривых участках пути с ускорениями.

Другой подход, предлагаемый в данной работе, заключается в преобразовании уравнений движения, исключающем вра-

щение колесной пары на стадии интегрирования. Данное решение, по сути, является разновидностью подхода Эйлера.

Математическая модель упругого тела

В соответствии с подходом, лежащим в основе моделирования динамики упругих тел в ПК «УМ», колесная пара может произвольно перемещаться в пространстве как абсолютно твердое тело, при этом перемещения ее точек в результате деформаций полагаются малыми.

Кинематика упругой колесной пары. Положение произвольной точки К колесной пары относительно глобальной системы координат СК0 выражается суммой радиуса-вектора начала отсчета локальной СК1 и радиуса-вектора точки относительно СК1 (рис. 1):

гГ = Г1(0) + Ао1 Р1},

где А01 - матрица ориентации СК1 относительно СК0; векторы представлены в

системах координат, указанных верхним индексом в скобках.

Упругие перемещения точек колесной пары рассчитываются с применением метода конечных элементов и модального подхода. Перемещения узлов относительно локальной системы координат выражаются произведением модальной матрицы на матрицу-столбец модальных координат: н

х = 2 Ь ^ = Hw,

1=1

где х - N^1 матрица-столбец узловых координат; N - число степеней свободы ко-нечноэлементной модели; Ь- - 1-я форма упругой колесной пары; Wl ---я модальная координата; н - число используемых форм; Н - модальная матрица.

Рис. 1. Присоединенная система координат колесной Рис. 2. Локальные координаты точки К колесной пары пары

Если узел расположен в точке К, его упругие перемещения могут быть представлены произведением

d к = Н к w,

где Нк - часть модальной матрицы, соответствующая узлу К. Тогда выражение радиуса-вектора точки К может быть записано следующим образом (рис. 2):

гГ = Г(0) + Ао1(р^ + Нк w), где рР1 - постоянный в СК1 радиус-вектор точки К в недеформированном состоянии.

Упругие формы Ь- рассчитываются в соответствии с методом Крэйга - Бэмптона

[12; 13]. Вначале выбираются интерфейсные узлы в шарнирных точках и точках присоединения силовых элементов. Затем последовательно рассчитываются статические формы от единичных смещений и единичных поворотов по всем степеням свободы в интерфейсных узлах и собственные формы при закрепленных интерфейсных узлах. Число используемых собственных форм выбирается исследователем в зависимости от требуемого частотного диапазона.

Статическая форма от единичного поворота вокруг оси X

£ %

Собственная форма при закрепленных интерфейсных

узлах

Рис. 3. Примеры статической и собственной форм колесной пары

Примеры статической и собственной форм представлены на рис. 3. Выбраны два интерфейсных узла на торцах оси.

Чтобы исключить твердотельное движение относительно локальной системы координат, набор статических и собственных форм преобразуется на основе решения обобщенной проблемы собственных значений с редуцированными матрицами масс и жесткости колесной пары [20].

Уравнения движения. Упругая колесная пара представляется в ПК «УМ» упругой подсистемой, уравнения движения которой выводятся на основе единого подхода, независимо от области использования моделируемого объекта. В рамках данного подхода применяются уравнения Лагранжа второго рода с упрощенным выражением кинетической энергии, которое записывается с учетом допущения, что масса тела распределена по узлам конеч-ноэлементной сетки. Полный вывод уравнений приводится в работе [20].

Модель контактных сил. Моделирование железнодорожных экипажей требует точного вычисления контактных сил

Рис. 4. Пример результата расчета многоточечной контактной модели

между колесом и рельсом. Многоточечная

модель неэллиптического контакта применяется для расчета контактных сил в ПК «УМ» (рис. 4). Для решения нормальной задачи используется так называемая полу-герцевская контактная модель Пиотровского и Кика (1 рюйг^бы, К1к) [14]. Расчет касательных сил крипа выполняется на основе алгоритма РАЗТБГМ, разработанного Калкером (Ка1кег) [15].

Кинематика профиля колеса. Одной из основных проблем моделирования динамики упругой колесной пары является расчет контактных сил, которые следует рассматривать как подвижную нагрузку. То есть точки их приложения изменяются относительно локальной системы координат и в общем случае не совпадают ни с одним узлом. Кроме того, необходимо предложить способ описания кинематики профиля колеса с учетом упругости колесной пары.

Профиль колеса абсолютно твердой колесной пары в ПК «УМ» задается в системе координат СКР0 У0Ъ0 с началом отсчета в центральной точке на поверхности рельса (рис. 5). Ось Ъ0 параллельна оси Ъ СК пути с учетом возвышения рельса. Система координат профиля рельса УгЪг повернута на угол аг0 относительно СКР0 в соответствии с подуклонкой рельса. Положение и ориентация системы координат профиля колеса СКК относительно

СКР0 задаются сдвигами АЪ и АУ начала отсчета и углом Аа. Углы аг0 и Аа полагаются малыми.

Профиль колеса недеформируемый; положение любой точки профиля относи-

тельно глобальной СК0 определяется положением и ориентацией локальной системы координат колесной пары СК1 и положением начала отсчета СКК относительно СК1.

Рис. 5. Относительное положение и ориентация профилей колеса и рельса

Профиль упругой колесной пары, так же как и для абсолютно твердой пары, полагается недеформируемым. Однако его положение и ориентация рассчитываются методом наименьших квадратов с учетом упругих перемещений узлов на поверхности качения колеса.

Рассмотрим кратко данный подход. Трехмерная конечноэлементая модель колесной пары создается вращением половины её сечения, разбитого на плоские элементы, вокруг оси с некоторым угловым шагом Аа. Эта плоская сетка должна содержать узлы, расположенные точно на профиле колеса. Таким образом, КЭ модель включает узлы, принадлежащие профилю колеса, в каждом у'-м сечении, по-

вернутом на угол у'Аа. Положение и ориентация профиля колеса, находящегося в контакте с рельсом, рассчитываются с учетом упругих перемещений точек пересечения профиля с линиями КЭ сетки между двумя сечениями (рис. 6). Перемещения точек вычисляются путем интерполяции соответствующих значений в соседних узлах.

Введем следующие обозначения: р/0 - положение 1-й точки недеформированно-го профиля (рис. 7); А гг- - упругие перемещения точки; АХ, (А Я), Аа - неизвестные смещение и поворот профиля в результате деформаций.

Положения точек с учетом их упругих перемещений могут быть записаны как сумма:

Р = Р 0 +Аг .

С другой стороны, они могут быть выражены через положение и ориентацию профиля следующим образом:

Я/ = Р 0 +А Я -Р 0 е х А« ,

где ех - единичный вектор вдоль оси X,

ортогональной плоскости рисунка (направление движения экипажа). Тогда неизвестные значения А Я и Аа могут быть найдены из условия минимизации невязки: е(6Я,8а) = 2 (Я Т-Рр )(Я г-р).

Рис. 6. Точки пересечения профиля, контактирующего с рельсом, с линиями КЭ сетки

Расчет контактных сил. Конечно-элементная модель колесной пары состоит из восьмиузловых шестигранников (кубиков). Рассмотрим основные идеи для вы-

Рис. 7. Положение /-й точки профиля в СКК

числения обобщенных сил, соответствующих силе, приложенной в произвольной точке грани конечного элемента. Другими словами, сила, действующая в произволь-

ной точке поверхности, должна быть преобразована к узловым силам, а затем обобщенные силы рассчитываются с применением общего алгоритма для упругих тел.

В «УМ» реализованы два алгоритма расчета узловых сил: 1) упрощенный алгоритм, рассматривающий только геометрию полигона, к которому приложена сила; 2) распределение приложенной силы с использованием функций форм конечного элемента.

Согласно упрощенному алгоритму, внешняя сила распределяется между узлами полигона в два этапа обратно пропорционально расстояниям от точки приложения до узлов (рис. 8).

Для описания алгоритма распределения силы с использованием функций форм введем систему безразмерных координат П, С с началом отсчета в центре шестигранника (рис. 9) [16].

Ь,2 Ь22 \Ек2

г

Ь21 Л, к21

Рис. 8. Алгоритм упрощенного расчета узловых сил

Значения безразмерных координат в узлах элемента равны ±1. Перемещения точек элемента вдоль осей х, у и г обозна-

Рис. 9. Распределение приложенной силы с использованием функций форм конечного элемента

чим соответственно и, V и н. Тогда поля перемещений определяются следующими выражениями:

где и е = [и!,..., Щ~]Т , V е = [V!,..., V8]Т , W в = ^]Т , N = N8],

Ni = 1/80 + ££ )0 + П 1, )0 + , г=1. 8, например N3 = 1/8(1 - £)(1 + п)(1 - О .

Узловые силы от силы, приложенной в точке К, рассчитываются по следующим формулам:

Рех = N (£ , ПК ,Ск )Рх , К = NТ , Пк ,Ск )РУ , ^ = NТ (£ , П К ¿К Жг ,

где Гех = [Fx1,..., Рх8]Т , Реу = [Fy1,..., ^Г ,

У1'

у8 -

= [ fл,..., ^8]т; £к , Пк ,Ск - безразмерные координаты точки К; р1хх,¥1уу,¥1^2 -

проекции Гк. Координаты £К,пк,СК рассчитываются по известным значениям декартовых координат точки К с использованием итерационного алгоритма типа метода Ньютона.

Моделирование динамики колесной пары без вращения конечноэлементной сетки. Уравнения движения интегрируются неявным многошаговым методом типа метода Парка [17; 18]. Этот метод реализует алгоритм прогноза-коррекции.

Рассмотрим моделирование упругой колесной пары с невращающейся конечно-элементной сеткой. Данный прием имитирует подход Эйлера на этапе интегрирования уравнений движения.

Перед расчетом прогнозных значений обобщенных координат на новом шаге интегрирования узловые координаты поворачиваются на угол da. Используется прогнозное приращение угла поворота колесной пары вокруг оси вращения.

Пусть ёа = еАа , где А а - угловой шаг конечноэлементной сетки колесной тары (рис. ШХ ее [Мт^ етах =1 для линейной интерполяции, которая рассматривается ниже.

2

1

А г, Н1

А г = =

А ^ Н N

Рис. 10. К описанию подхода Эйлера

Интерполяция применятся к упругим перемещениям узлов сетки. Пусть А гк -упругое перемещение узла с номером к на предыдущем шаге интегрирования, А гк -

перемещение точки к' после вращения колеса на da (рис. 10). Тогда можно записать следующее выражение для интерполяции:

А гк = аА гк + ЬА А А гр

Здесь р(к) - предыдущий узел для узла с номером к, то есть узел, который займет место узла к после поворота колеса на Аа ; А А = Ау (Аа) - матрица поворота

вокруг оси колесной пары Аа . Для линейной интерполяции значения коэффициентов: а = 1 -е, Ь = е.

Упругие перемещения зависят от модальных координат следующим образом:

= .

По столбцам матрицы Н располагаются преобразованные статические и собственные формы, рассчитанные в соответствии с методом Крэйга - Бэмптона. Они взаимно ортогональны и нормализованы в М-норме, то есть НтМН = Е , где М - матрица масс конечноэлементной модели колесной пары. Тогда выражения (1) и (2) принимают следующий вид:

А Ч = fа Н к + ЬА АН р (к)) w, А г = 1 а Н+ Ь Н I w = Hw .

Поскольку модальные координаты колесной пары являются частью набора обобщенных координат всего экипажа, их значения W должны быть рассчитаны после поворота на da при условии минимизации невязки

ш^Шн^ -

в некоторой норме. Выбор М-нормы в данном случае кажется разумным, то есть

пип^ ТНТ - w1H1) M(HW - Hw) = ш_т(/ (W)).

д д

Чтобы найти решение, необходимо рассчитать экстремум функции /(W):

^ = 2НТМ(Н^ - Н^ = 2(W - HTмHw) = 0.

дд

Окончательно выражение преобразованных модальных координат принимает следующий вид:

ш = Н^МН№ = Н^М(аН+ ЬНp)w = ауя+ЬХр Матрица X должна быть рассчитана перед началом моделирования:

т N т

X р = Нт МН р = 1^ Нк А АН р

к=1

Модели пути. В ПК УМ реализованы три модели пути, которые отличаются уровнем детализации путевой структуры:

- невесомый рельс;

- инерционный рельс, моделируемый абсолютно твердым телом;

- подробная (упругая) модель пути.

V

Wheel-rail penetration contact model

~3*_i г 7L dj ' 4

5!

Рис. 11. Модель упругого пути: 1 - рельс;

2 - крепление рельса; 3 - твердая шпала;

4 - прокладка под шпалой; 5 - абсолютно твердое или упругое основание

Невесомый рельс применяется для исследований в частотном диапазоне 0...30 Гц. Инерционный рельс обеспечивает достоверное моделирование в диапазоне до 100 Гц, упругий путь - до 1000 Гц.

Упругая колесная пара может взаимодействовать с любой из упомянутых моделей; в данной работе приведены результаты, полученные при движении по упругому пути.

Подробная трехмерная модель пути включает упругие рельсы, прокладки, скрепления, шпалы и подшпальное основание (рис. 11). Рельсы моделируются балками Тимошенко. Шпалы могут представляться абсолютно твердыми телами. Скре-

пления и подшпальное основание моделируются с помощью нелинейных силовых элементов. Верификация упругой модели пути представлена в работе [19]. Результаты моделирования

В этом разделе рассмотрен тестовый пример. Представлены результаты моделирования автомотрисы АС4 с постоянной скоростью 80 км/ч в криволинейном участке пути.

Модель автомотрисы включает следующие элементы (рис. 12, 13):

- упругую переднюю колесную пару;

- абсолютно твердую заднюю колесную пару;

- кузов и четыре тела, моделирующие буксы;

- один шарнир, задающий шесть координат кузова, и четыре вращательных шарнира, соединяющих буксы с колесными парами;

- 12 биполярных силовых элементов: четыре наклонных и четыре поперечных гасителя колебаний, четыре продольные тяги;

- восемь силовых элементов, моделирующих пружины подвески.

Рис. 12. Общий вид модели автомотрисы

Рис. 13. Модель подвески автомотрисы: 1 - колесная пара; 2 - пружина подвески; 3 - наклонные фрикционные гасители; 4 - поперечный демпфер; 5 - букса; 6 - продольная тяга

Основные параметры модели упругой колесной пары представлены в таблице. Основные параметры модели упругой колесной пары

Таблица

Параметр Значение Параметр Значение

Число конечных элементов 46 442 Число упругих форм 35

Число узлов 56 561 Низшая частота 29,6 Hz

Число степеней свободы 339 366 Высшая частота 5456,6 Hz

Угловой шаг конечноэлементой сетки 2° Коэффициент критического демпфирования для каждой формы 0,01

Модель упругого пути включает 1070 балочных конечных элементов, моделирующих рельсы, и 1072 тела, моделирующих шпалы. Длина каждого балочного элемента 0,6 м соответствует межшпаль-ному расстоянию. Модель упругого пути имеет 9648 степеней свободы. Крепления представлены линейными силовыми элементами.

Геометрия пути показана на рис. 14. Он включает начальный прямой участок, первую переходную кривую, кривую постоянного радиуса, вторую переходную

примера

Контрольный узел

Рис. 16. Контрольный узел

Динамические параметры рассчитываются в точках, положение которых может быть задано двумя способами. Первый способ: координаты точки задаются относительно локальной СК колесной пары и, как правило, совпадают с координатами выбранного контрольного узла в недефор-мированном состоянии. Таким образом, точка вращается вместе с локальной СК.

кривую и конечный прямой участок. Параметры, указанные на рисунке, имеют следующие значения: Ь0 = 10 м, Р1 = 50 м, 5 = 200 м, Я = 300 м, Р2 = 50 м, Ь1 = 1000 м.

Движение автомотрисы исследуется с учетом неровностей пути. Они сгенерированы по спектральным плотностям мощности и соответствуют плохому состоянию пути согласно нормам ШС (Международное объединение железных дорог) (рис. 15).

35 2(1 15 111

0

Рис. 15. Неровности пути: 1 - левый рельс, вертикальные; 2 - правый рельс, вертикальные; 3 - оба рельса, горизонтальные

Согласно второму способу, параметры вычисляются в заданной точке пространства, связанного с колесной парой. Ее координаты не зависят от угла поворота колесной пары вокруг поперечной оси У. Например, если выбран контрольный узел, расположенный в начальный момент времени на вертикальном радиусе (рис. 16), то в соответствии со вторым способом параметры вычисляются в точке колеса, находящейся на вертикальном радиусе в данный момент времени. Первый способ будем называть «вращающаяся сетка», второй - «невра-щающаяся сетка».

Результаты моделирования представлены на рис. 17, 18, 19.

Рис. 17. Радиальные упругие перемещения в контрольном узле: 1 - вращающаяся сетка, подход Ла-гранжа; 2 - вращающаяся сетка, подход Эйлера; 3 -невращающаяся сетка, подход Лагранжа; 4 - невра-щающаяся сетка, подход Эйлера

Рис. 18. Эквивалентные напряжения фон Мизеса в контрольном узле: 1 - вращающаяся сетка, подход Лагранжа; 2 - вращающаяся сетка, подход Эйлера; 3 -невращающаяся сетка, подход Лагранжа; 4 - невра-щающаяся сетка, подход Эйлера

Заключение

Рис. 19. Спектральная плотность мощности ускорения в контрольном узле

Применение двух предложенных в работе методов моделирования динамики упругих колесных пар приводит практически к одинаковым результатам расчета упругих перемещений и напряжений в контрольном узле.

Спектр ускорений, рассчитанных в контрольном узле, лежит в частотном диапазоне от 0 до 1300 Гц (рис. 19). Подобные результаты не могут быть получены с использованием твердотельной колесной пары.

Один численный эксперимент продолжается от 72 до 75 минут при использовании параллельных вычислений на компьютере с процессором Intel Core i7 3,5 ГГц. Моделирование с применением подхода Лагранжа выполняется быстрее по сравнению с подходом Эйлера на 8.. .10 %. Повышение эффективности расчетов с использованием подхода Эйлера является предметом будущих работ.

Работа выполнена при поддержке Российского фонда фундаментальных исследований (РФФИ), грант 17-01-00815.

СПИСОК ЛИТЕРАТУРЫ

Morys, B. Enlargement of Out-of-Round Wheel Profiles on High Speed Trains / B. Morys // Journal of Sound and Vibration. - 1999. - Vol. 227 (5). - P. 965-978.

Szolc, T. Medium Frequency Dynamic Investigation of the Railway Wheelset-Track System Using a Discrete-Continuous Model / T. Szolc // Vehicle System Dynamics. - 1998. - Vol. 30. - P. 473-508.

3. Fayos, J. An Eulerian coordinate-based method for analyzing the structural vibrations of a solid of revolution rotating about its main axis / J. Fayos, L. Baeza, F.D. Denia, J.E. Tarancon // Journal of Sound and Vibration. - 2007. - Vol. 306. - P. 618635.

4. Baeza, L. High frequency railway vehicle-track dynamics through flexible rotating wheelsets /

L. Baeza, J. Fayos, A. Roda, R. Insa // Vehicle System Dynamics. - 2008. - Vol. 46 (7). - P. 647-659.

5. Vila, P. Simulation of the evolution of rail corrugation using a rotating flexible wheelset model / P. Vila, J. Fayos, L. Baeza // Vehicle System Dynamics. - 2001. - Vol. 49 (11). - P. 1749-1769.

6. Frigerio, M. Flexible wheelset-track model for numerical estimation of stresses in railway axles: PhD Thesis / M. Frigerio. - Politecnico di Milano, 2010.

7. Kaiser, I. Interaction of elastic wheelsets and elastic rails: modelling and simulation / I. Kaiser, K. Popp // Vehicle System Dynamic. - 2006. - Vol. 44. - P. 932-939.

8. Kaiser, I. On an ALE-approach for rotating elastic structures / I. Kaiser, A. Heckmann, F.v.d. Linden // Multibody dynamics: ECCOMAS Thematic Conference. - 2007.

9. Kaiser, I. Refining the modelling of vehicle-track interaction / I. Kaiser // Vehicle System Dynamics. - 2012. - Vol. 50. - P. 229-243.

10. Guiral, A. Vehicle-track interaction at high frequencies - Modelling of a flexible rotating wheelset in non-inertial reference frames / A. Guiral, A. Alonso, J.G. Giménez // Journal of Sound and Vibration. - 2015. - Vol. 355. - P. 284-304.

11. Shabana, A. Flexible Multibody Dynamics: Review of Past and Recent Developments / A.Shabana // Multibody System Dynamics. - 1997. - Vol. 1. - P. 189-222.

12. Craig, R.R.Jr. Coupling of substructures for dynamic analysis / R.R.Jr. Craig, M.C.C. Bampton // AIAA Journal. - 1997. - Vol. 6. - № 7. - P. 13131319.

13. Craig, R.R.Jr. Coupling of substructures for dynamic analysis: an overview / R.R.Jr. Craig // AIAA

1. Morys, B. Enlargement of Out-of-Round Wheel Profiles on High Speed Trains / B. Morys // Journal of Sound and Vibration. - 1999. - Vol. 227 (5). - P. 965-978.

2. Szolc, T. Medium Frequency Dynamic Investigation of the Railway Wheelset-Track System Using a Discrete-Continuous Model / T. Szolc // Vehicle System Dynamics. - 1998. - Vol. 30. - P. 473-508.

3. Fayos, J. An Eulerian coordinate-based method for analyzing the structural vibrations of a solid of revolution rotating about its main axis / J. Fayos, L. Baeza, F.D. Denia, J.E. Tarancon // Journal of Sound and Vibration. - 2007. - Vol. 306. - P. 618635.

4. Baeza, L. High frequency railway vehicle-track dynamics through flexible rotating wheelsets / L. Baeza, J. Fayos, A. Roda, R. Insa // Vehicle System Dynamics. - 2008. - Vol. 46 (7). - P. 647-659.

5. Vila, P. Simulation of the evolution of rail corrugation using a rotating flexible wheelset model / P. Vila, J. Fayos, L. Baeza // Vehicle System Dynamics. - 2001. - Vol. 49 (11). - P. 1749-1769.

6. Frigerio, M. Flexible wheelset-track model for numerical estimation of stresses in railway axles: PhD Thesis / M. Frigerio. - Politecnico di Milano, 2010.

Paper, No 2000-1573. - AIAA Dynamics Specialists Conference, Atlanta.

14. Piotrowski, J. A simplified model of wheel/rail contact mechanics for non-Hertzian problems and its application in rail vehicle dynamic simulations / J. Piotrowski, W. Kik // Vehicle System Dynamics. -2008.- Vol. 45. - Is. 1-2. - Р. 27-48.

15. Kalker, J.J. A Fast Algorithm for the Simplified Theory of Rolling Contact / J.J. Kalker // Vehicle System Dynamics. - 1982. - Vol. 11. - Is. 1. - Р. 113.

16. Zienkiewicz, O.C. The Finite Element Method. Vol. I. Basic Formulations and Linear Problems / O.C. Zienkiewicz, R.L. Taylor. - 4th ed. - London: McGraw-Hill, 1989.

17. Park, K.C. An improved stiffly stable method for direct integration of nonlinear structural dynamic equations / K.C. Park // Journal of Applied Mechanics. - 1975. - Vol. 42 (2). - Р. 464-470.

18. Pogorelov, D. Differential-algebraic equations in multibody system modeling / D. Pogorelov // Numerical Algorithms. -1998. - Vol. 19. - Р. 183194.

19. Rodikov, A. Computer simulation of train-trackbridge interaction / A. Rodikov, D. Pogorelov, G. Mikheev, R. Kovalev, Q. Lei, Y. Wang // CORE 2016: Proceedings of International Conference on Railway Excellence (May 16-18, 2016, Melbourne, Australia).

20. Михеев, Г.В. Компьютерное моделирование динамики систем абсолютно твердых и упругих тел, подверженных малым деформациям: дис. ... канд. техн. наук / Г.В. Михеев. - Брянск, 2004. -153 с.

7. Kaiser, I. Interaction of elastic wheelsets and elastic rails: modelling and simulation / I. Kaiser, K. Popp // Vehicle System Dynamic. - 2006. - Vol. 44. - P. 932-939.

8. Kaiser, I. On an ALE-approach for rotating elastic structures / I. Kaiser, A. Heckmann, F.v.d. Linden // Multibody dynamics: ECCOMAS Thematic Conference. - 2007.

9. Kaiser, I. Refining the modelling of vehicle-track interaction / I. Kaiser // Vehicle System Dynamics. - 2012. - Vol. 50. - P. 229-243.

10. Guiral, A. Vehicle-track interaction at high frequencies - Modelling of a flexible rotating wheelset in non-inertial reference frames / A. Guiral, A. Alonso, J.G. Giménez // Journal of Sound and Vibration. - 2015. - Vol. 355. - P. 284-304.

11. Shabana, A. Flexible Multibody Dynamics: Review of Past and Recent Developments / A.Shabana // Multibody System Dynamics. - 1997. - Vol. 1. - P. 189-222.

12. Craig, R.R.Jr. Coupling of substructures for dynamic analysis / R.R.Jr. Craig, M.C.C. Bampton // AIAA Journal. - 1997. - Vol. 6. - № 7. - P. 13131319.

13. Craig, R.R.Jr. Coupling of substructures for dynamic analysis: an overview / R.R.Jr. Crag // AIAA Paper, No 2000-1573. - AIAA Dynamics Specialists Conference, Atlanta.

14. Piotrowski, J. A simplified model of wheel/rail contact mechanics for non-Hertzian problems and its application in rail vehicle dynamic simulations / J. Piotrowski, W. Kik // Vehicle System Dynamics. -2008.- Vol. 45. - Is. 1-2. - P. 27-48.

15. Kalker, J.J. A Fast Algorithm for the Simplified Theory of Rolling Contact / J.J. Kalker // Vehicle System Dynamics. - 1982. - Vol. 11. - Is. 1. - P. 113.

16. Zienkiewicz, O.C. The Finite Element Method. Vol. I. Basic Formulations and Linear Problems / O.C. Zienkiewicz, R.L. Taylor. - 4th ed. - London: McGraw-Hill, 1989.

17. Park, K.C. An improved stiffly stable method for direct integration of nonlinear structural dynamic equations / K.C. Park // Journal of Applied Mechanics. - 1975. - Vol. 42 (2). - P. 464-470.

18. Pogorelov, D. Differential-algebraic equations in multibody system modeling / D. Pogorelov // Numerical Algorithms. -1998. - Vol. 19. - P. 183194.

19. Rodikov, A. Computer simulation of train-trackbridge interaction / A. Rodikov, D. Pogorelov, G. Mikheev, R. Kovalev, Q. Lei, Y. Wang // CORE 2016: Proceedings of International Conference on Railway Excellence (May 16-18, 2016, Melbourne, Australia).

20. Mikheev, G.V. Computer simulation of the dynamics of systems of absolutely solid and elastic bodies subject to small deformations: dis. ... Cand. tech. Sciences / G.V. Mikheev. - Bryansk, 2004. - 153 p.

Статья поступила в редакцию 21.11.18 Рецензент: д.т.н., профессор Брянского государственного

технического университета Сакало В.И.

Статья принята к публикации 22. 03. 19.

Сведения об авторах:

Михеев Геннадий Викторович, к.т.н., директор ООО «Вычислительная механика», е-mail: mikheev @umlab. ru.

Погорелов Дмитрий Юрьевич, д.физ.-мат.н., профессор кафедры «Прикладная механика»

Mikheev Gennady Victorovich, Can. Sc. Tech., Director of PC "Computer Mechanics", е-mail: mikheev@umlab. ru.

Pogorelov Dmitry Yurievich, Dr. Phys-Math., Prof. of the Dep. "Applied Mechanics", Bryansk State Technical University, е-mail: pogorelov@umlab.ru.

Брянского государственного технического университета, е-mail: po gorelov@umlab.ru. Родиков Александр Николаевич, инженер-программист ООО «Вычислительная механика», е-mail: rodikov@umlab.ru.

Rodikov Alexander Nikolaevich, Engineer-Programmer of PC "Computer Mechanics", е-mail: rodikov@umlab.ru.