Число осей Масса, т База, мм Конструктивная
скорость, км/ч Расстояние от уровня головок рельсов до опорной поверхности
подпятника, м
|
2 4,8 1850 120 0,806
|
1.2
Анализ диапазона частот и амплитуд собственных колебаний объекта исследования
Данный полувагон модели 12-759 оснащен усиленными тележками модели
18-100. Тележка данной модели имеет максимальный статический прогиб и коэффициент динамичности .
Максимальная амплитуда при подпрыгивании вычисляется по формуле
(1)
где - статический прогиб, м;
- коэффициент динамичности.
Масса кузова рассчитывается по следующей формуле
, (2)
где -масса тары вагона, т;
- масса тележки 18-100, т.
Рассчитаем массу кузова в
гружёном состоянии
Частота колебаний кузова вагона при подпрыгивании вычисляется по
формуле
, (3)
где C-жесткость рессорного подвешивания тележки, Н/м;
- масса кузова вагона, кг.
По ТЭП тележки 18-100 жесткость одной пружины рессорного
подвешивания равна 7,5·105 (-0,3; +0,0). [4]
Определим жесткость рессорного
подвешивания тележки по формуле (4):
, (4)
где с - жесткость одной пружины, Н/м;
n -
количество пружин в рессорном комплекте;
m -
количество рессорных подвешиваний тележки.
Определим минимальную жесткость по
формуле для пяти пружин в рессорном комплекте
(5)
Определим максимальную жесткость по формуле
(6)
Определим минимальную жесткость по
формуле (5) для шести пружин в рессорном комплекте
Определим максимальную жесткость по формуле (6)
Частота колебаний для порожнего вагона составит.
Частота при жесткости C =
5250000 Н/м и пятью пружинами в рессорном комплекте
Частота при жесткости C = 7500000Н/м и пятью
пружинами в рессорном комплекте
Частота при жесткости C = 6300000 Н/м и шестью
пружинами в рессорном комплекте
Частота при жесткости C = 9000000 Н/м и шестью
пружинами в рессорном комплекте
Вычислим частоту колебаний кузова в груженном состоянии.
Частота при жесткости C = 5250000 Н/м и пятью пружинами
в рессорном комплекте
Частота при жесткости C =
7500000 Н/м и пятью пружинами в рессорном комплекте
Частота при жесткости C = 6300000 Н/м и шестью
пружинами в рессорном комплекте
Частота при жесткости C = 9000000 Н/м и шестью
пружинами в рессорном комплекте
Для вагона модели 12 - 759 пределы частоты колебаний должны
находятся в интервале от.
2. Разработка математической модели собственных
колебаний кузова вагона на рессорном подвешивании
2.1
Выбор и обоснование расчетной схемы
Расчетную схему для исследования собственных колебаний
подпрыгивания кузова вагона на пружинах рессорного комплекта представим виде
твердого тела, опирающегося на упругие элементы, жесткость которых
соответствует суммарной жесткости рессорного комплекта тележки. При колебаниях
подпрыгивания каждая точка твёрдого тела совершает плоскопараллельное движение.
Масса твердого тела равна массе кузова и приложена в геометрическом центре масс.
В качестве обобщенной координаты, характеризующий
колебательный процесс примем характеристику q, соответствующую
перемещению каждой точки твёрдого тела.
При составлении расчетной схемы примем следующие допущения
- Кузов считается абсолютно твердым телом;
- Масса кузова расположена в центре масс;
- Зазоры, износы в ходовых частях в местах
опирания кузова на тележку не учитываем;
- Зазоры в узле пятник - подпятник
отсутствуют;
- Считаем, что кузов постоянно связан с
ходовыми частями через упругие элементы;
- Считаем, что надрессорная балка
присоединена к кузову, поскольку расположена в схеме передачи нагрузки на
упругих элементах;
- Тормозное оборудование равномерно
распределено по длине кузова - кузов симметричное тело;
- Учитывая, что жесткость пути намного
больше жесткости рессорного подвешивания;
- Аэродинамическое сопротивление вагона при
движении не учитывается;
- Сила трения в фрикционных гасителях
колебаний равна 0. То есть работу фрикционных гасителей не учитываем, процесс
не затухающий;
- Суммарная жесткость рессорного
подвешивания двух тележек одинакова;
Согласно принципу Д'Аламбера, вырежем твёрдое тело (кузов
вагона), действие отброшенных связей заменим реакциями упругих элементов,
приложим к нему все внешние силы и силы инерции в центре масс. Сила инерции
направлена в сторону, противоположную направлению движения q.
Исходная схема для исследования собственных колебаний
подпрыгивания кузова вагона на рессорах будет иметь вид (рисунок 2.1).
Рисунок 2.1 - Исходная схема исследуемого объекта
Расчетная схема примет вид, изображенный на рисунке 2.2.
Рисунок 2.2 - Расчетная схема исследуемого объекта
2.2
Вывод уравнений математической модели
Математическая модель в динамике твердых тел представляет
собой систему, состоящую из уравнения движения и начальных условий.
Для получения математической модели используется принцип
Даламбера.
Алгоритм решения математической модели:
1 Непрерывный процесс колебаний делим на равные промежутки и
рассматриваем поведение систем именно в этих точках;
2 Поскольку тело находится в движении всегда
присутствует сила инерции, равная массе на ускорение, и приложенное к центру
масс в сторону противоположную движению;
Для каждого временного интервала составляем
уравнение равновесия;
Вводим ось движения, направленную вверх;
Проецируем все силы, действующие на кузов на ось
движения и получим уравнение движения
, (9)
где = - суммарная масса кузова
и груза;
- ускорение;
- сила реакции пружин.
Математическая модель в заданном движении твердых тел
включает в себя:
1 Уравнения движения;
2 Начальные условия.
Поскольку после установки груза в кузов, тело находится в
состоянии покоя (отсутствуют какие-либо перемещения, то есть , отсюда
(10)
Выразим и обозначим
Математическая модель собственных колебаний подпрыгивания
будет иметь вид
(12)
3. Выбор метода решения обыкновенных
дифференциальных уравнений
3.1
Анализ методов решения ОДУ
Существует множество методов решения ОДУ: Рунге-Кутта,
Эйлера-Коши, Милна, Адамса, Ньюмарка и разностные методы.
Самыми распространенными методами являются методы
Рунге-Кутта, Эйлера-Коши и разностные методы.
1 Метод Рунге-Кутта
Формулы этого метода предназначены для интегрирования
дифференциальных уравнений первого порядка вида
(8)
с начальными условиями
(9)
Существуют формулы Рунге - Кутта, предназначенные для
интегрирования дифференциального уравнения второго и третьего порядков, но они
очень громоздки и на практике используются очень редко.
Методы Рунге-Кутта обладают следующими отличительными свойствами
Они согласуются с рядом Тейлора вплоть до членов порядка , где степень р различна для
различных методов и называется порядком метода;
Эти методы являются одноступенчатыми: чтобы найти , нужна информация только о предыдущей
точке , ;
3 Не требуют вычисления производных от f (x, y), а требуют только
вычисления самой функции.
Именно благодаря третьему свойству методы Рунге-Кутта более
удобны для практических вычислений. Однако для вычисления одной последующей
точки решения приходится вычислять функцию f (x, y) несколько раз при
различных значениях х и у. Этот метод требует большой
квалификации и времени на отладки.
3 Итерационные
методы Эйлера-Коши
Итерационные методы Эйлера-Коши применимы для интегрирования
дифференциальных уравнений любого порядка.
Если имеем дифференциальное уравнение (10)
(10)
с начальными условиями (11)
(11)
то формула для интегрирования получается из следующих соображений.
Находят значения функции в виде трех членов разложения ряда Тейлора:
(12)
Затем вторая производная представляется в виде разностного выражения:
, (13)
Подставив выражения (12) в уравнение (13), получим
(14)
В выражении (12) значение неизвестно, поэтому, чтобы воспользоваться формулой (12),
строится итерационный процесс.
Таким образом, циклическая процедура вычислений по формуле (12)
представляет итерационный метод Эйлера-Коши для дифференциального уравнения
первого порядка типа (8).
Применение в программе ЭВМ таких, логических переходов позволяет
вести вычисления с автоматическим выбором шага интегрирования h в зависимости от заданной точности.
4 Разностные
методы интегрирования ОДУ
Эти методы основаны на замене производной в дифференциальных
уравнениях их приближенными разностными аналогами.
В итоге система дифференциальных уравнений преобразуется к
системе алгебраических уравнений. Алгоритм решения задачи прост, а программа
машинной реализации требует минимума времени при отладке. Однако разностные
методы являются неустойчивыми, т.е. решение зависит от величины шага
интегрирования. На практике в начале производят расчеты с целью определения шага
интегрирования. Критерием выбора является условие, когда разность значений двух
решений не превышает 5%.
Учитывая сказанное, в качестве метода решения полученной
математической модели примем разностный метод.
3.2
Описание алгоритма выбранного метода решения ОДУ
Для интегрирования полученного обыкновенного
дифференциального уравнения колебания кузова вагона на рессорах был выбран
разностный метод интегрирования.
Алгоритм решения разностным методом заключается в следующем:
1 Определяем наивысшую производную степень ОДУ. В нашем
случае - это 2 степень;
Определяем разностный аналог первой производной по
формуле (15):
, (15)
где - перемещение при текущем (i) значении времени;
- перемещение при предыдущем (i-1) значении времени;
h - шаг интегрирования или шаг разностной
аппроксимации по времени.
Определяем разностный аналог второй производной по формуле
(16)
(16)
где - перемещение при последующем (i+1) значении времени;
- перемещение при текущем значении времени;
- перемещение при предыдущем значении времени;
- шаг интегрирования или шаг разностной аппроксимации по времени.
(17)
Определяем начальные условия. Для этого выражаем из
уравнения (17) значение
(18)
(19)
Уравнение (19) является алгебраическим аналогом дифференциального
уравнения движения в системе (4).
С учетом выражения (19) и выражения (7) математическая модель
собственных колебаний подпрыгивания в разностной форме примет вид
(20)
Для Мк
Диапазон исследований для пяти
пружин в рессорном комплекте
Диапазон исследований для шести
пружин в рессорном комплекте
4. Разработка программы расчета собственных
колебаний кузова на рессорном подвешивании
Рисунок 3 - Блок-схема алгоритма решения задачи
Составляем таблицу идентификаторов (таблица 1).
Таблица 1 - Список идентификаторов
Обозначение в блок-схеме
|
Обозначение в программе
|
Наименование
|
Мк
|
Мк
|
Масса кузова
|
Мг
|
Mg
|
Масса груза
|
С
|
С
|
Жесткость рессорного комплекта
|
h
|
h
|
Шаг разностной Аппроксимации
|
g
|
g
|
Ускорение свободного падения
|
Tmax
|
Tmax
|
Время интегрирования
|
q0
|
q0
|
Начальные условия
|
qi-1
|
q1
|
Значение в предыдущий момент времени
|
qi+1
|
q2
|
Значение в следующий момент времени
|
.2 Исходный текст программы
Programvariant;
usescrt;
Varq0, qi, q1, q2, h, t, Tmax, C, g, m, Mk,
Mg:real;:text;;(‘BBOD
Mk’);(Mk);:=9.8;:=13000;:=0;:=Mk+Mg;:=(M*g)/(2*C);:=Mk;:=q0;:=q0;:= 2*qi - q1 +
sqr(h)*(g - (2*C*qi)/M);(q2:15:5, t:15:3);(F, q2:0:5);:=qi;:=q2;:=t +
h;t>=Tmax;;
End.
Пояснения к программе:
Программа - это упорядоченная последовательность команд на
выбранном языке программирования, соответствующая разработанному алгоритму.
Для выполнения программы предназначен оператор ввода Readln (читать). Как только при
выполнении программы встречается этот оператор, программа останавливается и
ожидает ввода числовых значений. После введения числовых значений выполнение
программы продолжается. Этот оператор сначала вводит значения, а затем
осуществляет переход на новую строку. В операторах в скобках указываются имена
переменных, которым будет принадлежать вводимые с клавиатуры значения.
Ввод данных:
Writeln - для ввода данных из памяти на экран, после
происходит переход на новую линию;
Var - для описания программа, с какими файлами она
взаимодействует;
Assign - для связывания ранее объявленной переменной с
именем создаваемого файла;
Rewrite - открыть файл для записи;
Sqr - квадрат числа.
Список литературы
колебание вагон дифференциальный нагруженность
1
Вагоны: Учебник для вузов ж.-д. трансп./ Под ред. Шадура. - 3-е изд., перераб.
и доп. - М.: Транспорт, 1980. - 439 с.
Вершинский
С.В., Данилов В.Н., Хусидов В.Д. Динамика вагона: Учебник для вузов ж.-д.
трансп. / Под ред. Вершинского. - 3-е изд., перераб. и доп.-М.: Транспорт,
1991-360 с.
Моделирование
собственных колебаний кузова вагона на рессорном подвешивании: метод. указания
/ В.Ф. Лапшин, Ю.Ю. Архипова. - Екатеринбург: УрГУПС, 2013. - 18.
Письменный
Д.П. Конспект лекций по высшей математике: полный курс - 4-е изд. - М.:
Айрис-пресс, 2006. - 608 с.: ил. - (Высшее образование).
Сенаторов
С.А. Прогнозирование нагруженности, износа и динамики подвижного состава. Часть
1-Динамические системы подвижного состава и методы их исследования.
Вагоны
и вагонное хозяйство: Методическое руководство к дипломному проектированию/
В.Ф. Лапшин, М.В. Орлов, А.Г. Пяткова и др.; Под общ. ред. проф. М.В. Орлова.
2-е изд., доп. и испр. Екатеринбург: Изд-во УрГУПС, 2005. - 120 с.
8
Конспект лекций по дисциплине «Математические модели вагонов и процессов».