Прямая задача динамики машины: определение закона движения при неустановившемся (переходном) режиме.
В отличие от установившегося режима движения режимы разгона и торможения называются неустановившимися. К этому режиму относят и режим движения "пуск-останов". Прямая задача динамики: определение закона движения машины при заданных внешних силовых воздействиях ( как сил и моментов сопротивления, так и движущих или управляющих сил ). Эта задача относится к задачам анализа, при которых параметры механизмов заданы, либо могут быть определены на предварительных этапах расчета. Для простоты и наглядности рассмотрим алгоритм решения этой задачи на примере конкретного механизма гидроподъемника. По условиям функционирования гидроподъемник за цикл движения должен переместить платформу 1 (рис. 7.6) на угол ∆φ1 и зафиксировать ее в конечном положении. При этом силы сопротивления определяются силами веса платформы и звеньев гидроцилиндра, движущие силы - давлением жидкости в цилиндре.
Алгоритм решения прямой задачи динамики при неустановившемся режиме.
Постановка задачи .
Дано: Кинематическая схема механизма и его размеры
lAB = 1 м, lBS1 = 2 м, lBD = 0.7м, lAC = 1.45м,
lBS2 = 0.35м, lBS3 = 0.4 м;
массы и моменты инерции звеньев m1 = 1000 кг,
IS1 = 800 кг * м 2, m2 = 50 кг, IS2 = 2 кг * м 2, m3 = 100 кг,
IS3 = 5 кг * м 2; w 1нач = 0, ∆φ1 = 30° , φ1нач = 0.
____________________________________________
Определить: w1 = f(φ1 ), t = f(φ1 ), w1 = f( t ), e1 = f(φ1 ).
1. Выбор динамической модели и определение ее параметров.
![]() |
Рис.7.7 |
В качестве динамической модели принимаем звено 1, совершающее вращательное движение вокруг точки А с круговой частотой w1 , положение которого определяется обобщенной координатой φ1 . Параметры динамической модели: суммарный приведенный момент инерции звеньев механизма Iпр∑ и суммарный приведенный момент, действующих на него внешних сил, Mпр∑определяются в следующей последовательности:
1.1. Определение кинематических передаточных функций для звеньев механизма u21 = u31 , центров масс VqS1 , VqS2 и VqS3 и точки приложения движущей силы VqD . Для определения этих функций воспользуемся методом проекций векторного контура механизма .
![]() |
Рис.7.8 |
Рассмотрим следующие векторные контуры:
l AB = l AC + l CB;
l AD = l AB + l BD;
l AS2 = l AC + l CS2;
l AS3 = l AC + l CS3;
l AS1 = xS1 + yS1 .
Для первого векторного контура l AB = l AC + l CB проекции на оси координат
Производные от этих выражений по φ1
позволяют определить первые передаточные функции
Для второго векторного контура l AD = l AB + l BD проекции на оси координат
Производные от этих выражений по φ1
позволяют определить первую передаточную функцию
![]() |
Для третьего векторного контура l AS2 = l AB + l BS2 проекции на оси координат
Производные от этих выражений
позволяют определить первую передаточную функцию
![]() |
Для четвертого векторного контура l AS3 = l AС + l С S3 проекции на оси координат
Производные от этих выражений
позволяют определить первую передаточную функцию
![]() |
Для последнего пятого векторного контура l AS1 = xS1 + yS1 проекции на оси координат
Производные от этих выражений по φ1
позволяют определить первую передаточную функцию
![]() |
Построим графики передаточных функций и передаточных отношений, которые необходимы для определения параметров динамической модели в нашем примере.
![]() |
Рис.7.9 |
1.2. Определение движущей силы по условиям в начале и в конце цикла.
Расчет проведем для закона изменения движущей силы, который изображен на рис.7.5. Величина движущей силы в начальном положении механизма рассчитывается по формуле
![]() |
Принимаем k=1.1 и получаем
В конечном положении величина движущей силы рассчитывается по формуле:
![]() |
Значение движущей силы в интервале ( β - α )* HD определим по формуле:
Примем α = 0.32 и β = 0.65 и рассчитаем перемещения центров масс
подставим полученные значения в формулу и получим
![]() |
![]() |
1.3. Определение приведенного суммарного момента .
![]() |
- 2. определение приведенного суммарного момента сил сопротивления
В нашем примере силами сопротивления являются силы веса звеньев механизма, поэтому расчет суммарного приведенного момента сил сопротивления проводим по формуле
![]() |
- определение приведенного момента движущей силы
В нашем примере только одна движущая сила, создаваемая давлением жидкости в гидроцилиндре. Приведенный момент от этой силы
![]() |
На рис. 7.13 приведены диаграммы приведенных моментов: сопротивления Мпр∑ с , движущего Мпр Fд i и суммарного Мпр∑ с = Мпр∑ + Мпр Fд i .
1.4. Определение суммарного приведенного момента инерции
В рассматриваемом механизме приведенный момент инерции суммируется из масс и моментов инерции звеньев и может быть рассчитан по следующей зависимости
![]() |
![]() |
Рис.7.14 |
![]() |
Рис.7.15 |
Графики переменной части суммарного приведенного момента инерции даны на рис. 7.13 и 7.14. Кроме того, имеется и постоянная часть Iпр∑ c, определяемая массой и моментом инерции звена 1
Суммарный приведенный момент инерции и равен сумме постоянной и переменной частей
.
2. Определение суммарной работы внешних сил.
Суммарную работу внешних сил получим интегрированием суммарного приведенного момента Мпр∑ по обобщенной координате dφ1
![]() |
Интегрирование можно проводить различными методами. Воспользуемся методом графического интегрирования. При этом методе участок изменения обобщенной координаты, на котором проводится интегрирование, разбирается на несколько малых частей (в нашем примере 6). В пределах каждого i -го участка кривая Мпр∑ = f (φ 1) заменяется прямой, соответствующей среднеинтегральному значению Мпр∑ i на этом участке. На продолжении оси абсцисс, влево от начала координат откладываем отрезок интегрирования k1 . Ординаты среднеинтегральных значений Мпр∑ i проецируем на ось ординат. Точки пересечения проецирующих линий с осью ординат соединяем прямыми с концом отрезка интегрирования. На диаграмме работы из начала первого участка (и до его конца) под углом ψ1 к оси абсцисс проводим прямую. Для второго участка аналогичная прямая проводится под углом ψ2.Ее начало выбирается в точке пересечения предыдущего отрезка прямой с вертикалью проходящей начало второго участка. Проведя построения для всего интервала интегрирования, получим график работы. Масштаб этого графика определим из подобия треугольников
Графики, иллюстрирующие построение диаграммы работы, приведены на рис.7.1 6 и 7.1 7
3. Определение угловой скорости звена приведения
Определение закона движения звена приведения в виде диаграммы изменения угловой скорости в функции обобщенной координаты w1= f(φ1) проводится по формуле
![]() |
![]() |
![]() |
![]() |
![]() |
Рис.7.18 |
Диаграмма w 1 = f (φ1 ) приведена на рис. 7.18.
4. Определение времени цикла.
Время цикла определяется по диаграмме t= f (φ1). Для построения этой диаграммы проведем интегрирование диаграммы угловой скорости
![]() |
Воспользуемся методом графического интегрирования обратной величины. При этом участок изменения обобщенной координаты, на котором проводится интегрирование, разбивается на несколько малых участков. В пределах каждого i -го участка кривая w1 = f (φ1) заменяется прямой, соответствующей среднеинтегральному значению w1ср i на этом участке. На оси ординат, откладываем отрезок интегрирования k2 (рис.7.19). Ординаты среднеинтегральных значений w1ср i проецируем на ось ординат. Точки пересечения проецирующих линий с осью ординат переносим по дугам окружности на продолжение оси абсцисс. Полученные на оси абсцисс точки, соединяем прямыми линиями с концом отрезка интегрирования. Из начала первого участка (на диаграмме времени) и до его конца под углом ψ1 к оси абсцисс проводим прямую линию. Для второго участка аналогичная прямая проводится под углом ψ2.Ее начало выбирается в точке пересечения предыдущего отрезка прямой с вертикалью проходящей начало второго участка. Проведя построения для всего интервала интегрирования, получим график времени. Масштаб этого графика определим из подобия треугольников
5. Построение диаграммы угловой скорости в функции времени
Диаграмма угловой скорости w1 = f ( t ) в функции времени строится по диаграммам w1 = f (φ1 ) и t= f (φ1 ), исключением переменной φ1 .
![]() |
![]() |
![]() |
6. Определение углового ускорения звена приведения
Для расчета углового ускорения звена приведения ε1 = f(φ1) можно воспользоваться двумя различными зависимостями:
Применение первой формулы приводит к большим погрешностям, так как она основывается на использовании одной из конечных зависимостей расчета w1 = f (φ1 ). Кроме того, в точках с нулевыми значениями w 1расчет по этой формуле дает неверный результат ε1 = 0. Поэтому проведем расчет зависимости ε1 = f(φ1) по второй формуле . Диаграмма функции ε1 = f(φ1) приведена на рис. 7.22.
![]() |
Дата добавления: 2020-02-05; просмотров: 660;