Ввод типов (Т. Е. Их идентификаторов) элементов связи и их параметров

===> 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

· Вычисление собственных чисел и собственных векторов. » EIGEN

СОБСТВЕННЫЕ ЧИСЛА С ПОЛОЖИТЕЛЬНОЙ ВЕЩЕСТВЕННОЙ ЧАСТЬЮ (НЕУСТОЙЧИВОСТЬ)

НОМЕР RE ЧАСТЬ IM ЧАСТЬ FN(Гц) F0(Гц) К. ДЕМПФ.

WR WI WI/2PI WI/|WR,WI| WR/|WR,WI|

-------------------------------------------------------------------------

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-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.