Ввод типов (т.е. их идентификаторов) элементов связи и их параметров
===> 1. ИДЕНТИФИКАТОРЫ ТИПА (MAX=4) ДЛЯ СВЯЗИ 5 (колесо-рельс:левая)
-1 = ПЕРЕЧИСЛИТЬ ВСЕ ВОЗМОЖНЫЕ
>> 21 0 0 0
===> 4. ОПЦИЯ ЗАДАНИЯ ОРИЕНТАЦИИ ВСПОМОГАТЕЛЬНОЙ СИСТЕМЫ КООРДИНАТ
ДЛЯ ЭЛЕМЕНТА СВЯЗИ: KENO = 0/1/2/3/4
5. 0<KENO<4: УГОЛ (KENO),
KENO=4 : БАЗИСНЫЕ ВЕКТОРЫ EPS1(1,2,3), EPS2(1,2,3)
ДЛЯ ТЕЛА 3 И СВЯЗИ 4 (колесо-рельс лев.)
>> 4 1.000 0.000 0.000 0.000 -1.000 0.0
===> 4. ОПЦИЯ ЗАДАНИЯ ОРИЕНТАЦИИ ВСПОМОГАТЕЛЬНОЙ СИСТЕМЫ КООРДИНАТ
ДЛЯ ЭЛЕМЕНТА СВЯЗИ: KENO = 0/1/2/3/4
5. 0<KENO<4: УГОЛ (KENO),
KENO=4 : БАЗИСНЫЕ ВЕКТОРЫ EPS1(1,2,3), EPS2(1,2,3)
ДЛЯ ТЕЛА 3 И СВЯЗИ 4 (колесо-рельс прав.)
>> 4 -1.000 0.000 0.000 0.000 1.000 0.0
· Задание движения отсчетной системы координат: прямолинейное со скоростью 35 м/с. » REFSYS
· Задание симметричных роторов, чье вращение влияет на динамику системы. » ROTOR
· Вычисление сил взаимодействия в связях. » NOMKRA
· Для дальнейших вычислений необходимо задать дополнительные степени свободы рамы. » KONFIG
===>1. ОПЦИЯ, KSKK(1:6), СТЕПЕНИ СВОБОДЫ ТВЕРДОГО ТЕЛА (ТРАНСЛ./ПОВОРОТ.)=0/1
2. ОПЦИЯ, KENHIB, ВСПОМОГАТЕЛЬНАЯ СИСТЕМА ОТСЧЕТА = 0/1
3. ИДЕНТИЧНОСТЬ С ДРУГИМИ ТЕЛАМИ = 0/>0
ТЕЛО 1 (РАМА )
>> 0 1 1 1 0 0 0 0
ТЕЛО 2 (КОЛЕСН ПАРА )
>> 1 1 1 1 1 0 0 0
Исследование устойчивости на прямом пути
Прежде, чем вычислять собственные числа системы уравнений, описывающей движение многомассовой системы, необходимо сгенерировать ее матрицу в равновесной конфигурации при скорости 35 м/с.
В отличие от модуля NOMKRA, который вычисляет силы, необходимые для приведения системы в состояние статического равновесия, основной модуль NEUKON используется для отыскания нового положения статического равновесия (соответствующих статическим силам перемещений) с учетом сил крипа, действующих между движущейся колесной парой и элементом пути.
· Отыскание новой конфигурации многомассовой системы. » NEUKON
· Генерация уравнений движения многомассовой системы, т.е. вычисление матриц системы и вектора обобщенных сил. » GEKSYS
· Вычисление собственных чисел и собственных векторов. »EIGEN
СОБСТВЕННЫЕ ЧИСЛА С ПОЛОЖИТЕЛЬНОЙ ВЕЩЕСТВЕННОЙ ЧАСТЬЮ (НЕУСТОЙЧИВОСТЬ)
НОМЕР RE ЧАСТЬ IM ЧАСТЬ FN(Гц) F0(Гц) К. ДЕМПФ.
WR WI WI/2PI WI/|WR,WI| WR/|WR,WI|
-------------------------------------------------------------------------
D-07 0. 0. 0. 1.
2/ 3 1.0599D-15 +/- 2.2206D-07 0.00 0.00 0.0000
4/ 5 1.5566D-16 +/- 3.6524D-07 0.00 0.00 0.0000
СОБСТВЕННЫЕ ЧИСЛА С НЕПОЛОЖИТЕЛЬНОЙ ВЕЩЕСТВЕННОЙ ЧАСТЬЮ (УСТОЙЧИВОСТЬ)
НОМЕР RE ЧАСТЬ IM ЧАСТЬ FN(Гц) F0(Гц) К. ДЕМПФ.
WR WI WI/2PI WI/|WR,WI| WR/|WR,WI|
-------------------------------------------------------------------------
6 0.0000D+00 0. 0. 0. 1.
7 0.0000D+00 0. 0. 0. 1.
D-07 0. 0. 0. 1.
9 -5.1669D+01 0. 0. 0. 1.
10 -3.4962D+02 0. 0. 0. 1.
11 -3.8063D+02 0. 0. 0. 1.
12 -1.0194D+03 0. 0. 0. 1.
13 -1.3370D+03 0. 0. 0. 1.
14 -1.5483D+03 0. 0. 0. 1.
15/ 16 -3.1250D+00 +/- 4.4612D+01 7.10 7.12 0.0699
17/ 18 -4.2083D+01 +/- 3.7705D+01 6.00 8.99 0.7448
19/ 20 -4.2487D+01 +/- 3.7656D+01 5.99 9.04 0.7484
21/ 22 -3.6696D+00 +/- 7.6500D+01 12.18 12.19 0.0479
Нулевые собственные числа системы соответствуют вращению колесной пары вокруг своей оси.
Поиск квазистатического положения равновесия экипажа в кривой
· Задание движения отсчетной системы координат: кривая постоянного радиуса со скоростью 20 м/с. » REFSYS
· Отыскание новой конфигурации многомассовой системы. » NEUKON
РАДИУС КРИВОЙ = 4.00D+02 M
ПОДЪЕМ ВИРАЖА = 0.200 М
СКОРОСТЬ = 20.00 М/С**2
ПОПЕРЕЧНОЕ УСКОРЕНИЕ = -0.32 М/С
ПРИРАЩЕНИЕ УСИЛИЯ = 139.25 Н
ЧИСЛО ИТЕРАЦИЙ = 2
ТЕОРИЯ КРИПА : FASTSIM ЛИНИЙ СЕТКИ: 10 11
ПОЛНАЯ МОДЕЛЬ
--------------------------------------------------------------------------
ТЕЛО X<MM> Y<MM> Z<MM> FX<MRAD> FY<MRAD> FZ<MRAD>
1 0.000 -4.617 0.156 0.435 0.000 0.000
2 -0.024 -4.297 0.059 0.499 0.000 0.000
3 0.000 -0.461 0.057 -0.022 0.000 0.000
Исследование устойчивости системы в кривой
· Выбор конфигурации, для которой производится исследование устойчивости.
===> ОПЦИИ ЧТЕНИЯ ДАННЫХ
БУДУТ ЧИТАТЬСЯ ИСХОДНЫЕ ДАННЫЕ
БУДУТ ЧИТАТЬСЯ ДАННЫЕ, СОЗДАННЫЕ NEUKON
>> 2 * ЧИТАЕМ ИЗ NEUKON
· Генерация уравнений движения многомассовой системы, т.е. вычисление матриц системы и вектора обобщенных сил. » GEKSYS
· Вычисление собственных чисел и собственных векторов. » EIGEN
СОБСТВЕННЫЕ ЧИСЛА С ПОЛОЖИТЕЛЬНОЙ ВЕЩЕСТВЕННОЙ ЧАСТЬЮ (НЕУСТОЙЧИВОСТЬ)
НОМЕР RE ЧАСТЬ IM ЧАСТЬ FN(Гц) F0(Гц) К. ДЕМПФ.
WR WI WI/2PI WI/|WR,WI| WR/|WR,WI|
-------------------------------------------------------------------------
D-02 0.0
D-02 0.0
СОБСТВЕННЫЕ ЧИСЛА С НЕПОЛОЖИТЕЛЬНОЙ ВЕЩЕСТВЕННОЙ ЧАСТЬЮ (УСТОЙЧИВОСТЬ)
НОМЕР RE ЧАСТЬ IM ЧАСТЬ FN(Гц) F0(Гц) К. ДЕМПФ.
WR WI WI/2PI WI/|WR,WI| WR/|WR,WI|
-------------------------------------------------------------------------
D-07 0. 0. 0. 1.
D-07 0. 0. 0. 1.
D-02 0. 0. 0. 1.
D-02 0. 0. 0. 1.
7 -5.1669D+01 0. 0. 0. 1.
8 -3.4152D+02 0. 0. 0. 1.
9 -3.8063D+02 0. 0. 0. 1.
10 -1.0194D+03 0. 0. 0. 1.
11 -1.3493D+03 0. 0. 0. 1.
12 -1.5483D+03 0. 0. 0. 1.
13/ 14 -1.2493D-07 +/- 8.9429D-03 0.00 0.00 0.0000
15/ 16 -3.1250D+00 +/- 4.4612D+01 7.10 7.12 0.0699
17/ 18 -4.2082D+01 +/- 3.7705D+01 6.00 8.99 0.7448
19/ 20 -4.2487D+01 +/- 3.7656D+01 5.99 9.04 0.7484
21/ 22 -3.6697D+00 +/- 7.6500D+01 12.18 12.19 0.0479
3.3 Модель с квазилинейным элементом «колесо – рельс»
Задание
1. Исследовать устойчивость движения со скоростью 35 м/с.
2. Исследовать частотные характеристики колесной пары в диапазоне от 0,1 до 100 Гц под действием поперечного гармонического возмущения левого и правого рельсов.
3. Провести стохастический анализ системы и определить спектральную плотность мощности выходных переменных в диапазоне частот от 0,1 до 100 Гц. Вычислить коэффициенты плавности хода по системе ORE в вертикальном и горизонтальном направлениях.
4. Вычислить усилия в контакте колесо – рельс при детерминированных неровностях пути, заданных как возвышение одной рельсовой нити над другой в виде полуволны синуса.
3.3.2 Генерация модели: описание конфигурации и параметров
· Запуск программы MEDYNA и инициализация файла данных для модели 1000.
· Копирование базы данных модели 1002 в файл данных модели 1000.
· Название модели. » MODELL
· Подготовка таблицы геометрии контакта для квазилинейного элемента связи колесо – рельс (тип 22). » RSGEO
· Задание конфигурации тел многомассовой системы. » KONFIG
Модель 1000 аналогична модели 1002. Степени свободы рамы в ней также должны быть ограничены вертикальными перемещениями для вычисления усилий в исходной конфигурации:
· Задание связей в системе. » VERBIN
В этом случае необходимо принять во внимание, что вспомогательная система координат для квазилинейного элемента связи 22 имеет ориентацию, отличную от нелинейного элемента связи 21.