Алгоритм решения прямой задачи динамики при неустановившемся режиме.

Постановка задачи .

Дано: Кинематическая схема механизма и его размеры

l_AB = 1 м, l_BS1 = 2 м, l_BD = 0.7м, l_AC = 1.45м,

l_BS2 = 0.35м, l_BS3 = 0.4 м;

массы и моменты инерции звеньев m₁ = 1000 кг,

I_S1 = 800 кг * м ², m₂ = 50 кг, I_S2 = 2 кг * м ², m₃ = 100 кг,

I_S3 = 5 кг * м ²; w _1нач = 0, D j ₁ = 30° , j _1нач = 0.

____________________________________________

Определить: w ₁ = f(j ₁ ), t = f(j ₁ ), w ₁ = f( t ), e ₁ = f(j ₁ ).

1. Выбор динамической модели и определение ее параметров.

рис. 7.7

В качестве динамической модели принимаем звено 1, совершающее вращательное движение вокруг точки А с круговой частотой w ₁ , положение которого определяется обобщенной координатой j ₁ . Параметры динамической модели: суммарный приведенный момент инерции звеньев механизма I^прå и суммарный приведенный момент, действующих на него внешних сил, M^прå определяются в следующей последовательности:

1.1. Определение кинематических передаточных функций для звеньев механизма u₂₁ = u₃₁ , центров масс V_qS1 , V_{qS2 и} V_qS3 и точки приложения движущей силы V_qD . Для определения этих функций воспользуемся методом проекций векторного контура механизма .

рис. 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 = x_S1 + y_S1 .

Для первого векторного контура l _AB = l _AC + l _CB проекции на оси координат

Производные от этих выражений по j ₁

позволяют определить первые передаточные функции

Для второго векторного контура l _AD = l _AB + l _BD проекции на оси координат

Производные от этих выражений по j ₁

позволяют определить первую передаточную функцию

Для третьего векторного контура l _AS2 = l _AB + l _BS2 проекции на оси координат

Производные от этих выражений

позволяют определить первую передаточную функцию

Для четвертого векторного контура l _AS3 = l _AС + l _{С S3} проекции на оси координат

Производные от этих выражений

позволяют определить первую передаточную функцию

Для последнего пятого векторного контура l _AS1 = x_S1 + y_S1 проекции на оси координат

Производные от этих выражений по j ₁

позволяют определить первую передаточную функцию

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

рис. 7.9

1.2. Определение движущей силы по условиям в начале и в конце цикла.

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

Принимаем k=1.1 и получаем

В конечном положении величина движущей силы рассчитывается по формуле:

Значение движущей силы в интервале ( b - a )* H_D определим по формуле:

F

Примем a = 0.32 и b = 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. Определение суммарной работы внешних сил.

Суммарную работу внешних сил получим интегрированием суммарного приведенного момента М_прå по обобщенной координате dj ₁

Интегрирование можно проводить различными методами. Воспользуемся методом графического интегрирования. При этом методе участок изменения обобщенной координаты, на котором проводится интегрирование, разбирается на несколько малых частей (в нашем примере 6). В пределах каждого i -го участка кривая М_прå = f (j ₁) заменяется прямой, соответствующей среднеинтегральному значению М_{прå i} на этом участке. На продолжении оси абсцисс, влево от начала координат откладываем отрезок интегрирования k₁ . Ординаты среднеинтегральных значений М_{прå i} проецируем на ось ординат. Точки пересечения проецирующих линий с осью ординат соединяем прямыми с концом отрезка интегрирования. На диаграмме работы из начала первого участка (и до его конца) под углом y ₁ к оси абсцисс проводим прямую. Для второго участка аналогичная прямая проводится под углом y ₂.Ее начало выбирается в точке пересечения предыдущего отрезка прямой с вертикалью проходящей начало второго участка. Проведя построения для всего интервала интегрирования, получим график работы. Масштаб этого графика определим из подобия треугольников

Графики, иллюстрирующие построение диаграммы работы, приведены на рис.7.1 6 и 7.1 7

3. Определение угловой скорости звена приведения

Определение закона движения звена приведения в виде диаграммы изменения угловой скорости в функции обобщенной координаты w ₁= f(j ₁) проводится по формуле

рис. 7.18

Диаграмма w ₁ = f (j ₁ ) приведена на рис. 7.18.

4. Определение времени цикла.

Время цикла определяется по диаграмме t= f (j ₁). Для построения этой диаграммы проведем интегрирование диаграммы угловой скорости

Воспользуемся методом графического интегрирования обратной величины. При этом участок изменения обобщенной координаты, на котором проводится интегрирование, разбивается на несколько малых участков. В пределах каждого i -го участка кривая w ₁ = f (j ₁) заменяется прямой, соответствующей среднеинтегральному значению w _{1ср i} на этом участке. На оси ординат, откладываем отрезок интегрирования k₂ (рис.7.19) . Ординаты среднеинтегральных значений w _{1ср i} проецируем на ось ординат. Точки пересечения проецирующих линий с осью ординат переносим по дугам окружности на продолжение оси абсцисс. Полученные на оси абсцисс точки, соединяем прямыми линиями с концом отрезка интегрирования. Из начала первого участка (на диаграмме времени) и до его конца под углом y ₁ к оси абсцисс проводим прямую линию. Для второго участка аналогичная прямая проводится под углом y ₂.Ее начало выбирается в точке пересечения предыдущего отрезка прямой с вертикалью проходящей начало второго участка. Проведя построения для всего интервала интегрирования, получим график времени. Масштаб этого графика определим из подобия треугольников

5. Построение диаграммы угловой скорости в функции времени

Диаграмма угловой скорости w ₁ = f ( t ) в функции времени строится по диаграммам w ₁ = f (j ₁ ) и t= f (j ₁ ), исключением переменной j ₁ .

6. Определение углового ускорения звена приведения

Для расчета углового ускорения звена приведения e ₁ = f(j ₁) можно воспользоваться двумя различными зависимостями:

Применение первой формулы приводит к большим погрешностям, так как она основывается на использовании одной из конечных зависимостей расчета w ₁ = f (j ₁ ). Кроме того, в точках с нулевыми значениями w ₁расчет по этой формуле дает неверный результат e ₁ = 0. Поэтому проведем расчет зависимости e ₁ = f(j ₁) по второй формуле . Диаграмма функции e ₁ = f(j ₁) приведена на рис. 7.22.