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

_________________________________________________________________

Определить: закон движения машины w ₁ = f(φ₁) и ε₁ = f(φ₁), момент инерции маховика I_доп , обеспечивающий заданную неравномерность вращения [δ].

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

Рассмотрим следующие векторные контуры, изображенные на рис. 8.4 рядом со схемой механизма:

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

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

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

Рис. 8.5

1.2. Определение приведенного момента движущих сил М^пр_д .

Индикаторную диаграмму (рис.8.3) строим по заданным значениям давления в цилиндре двигателя. Отрезок хода поршня Н_C* m _i делим на 10 интервалов. В каждой точке деления строим ординату диаграммы, задавшись (при p_i /p_max = 1) максимальной ординатой y_pmax . Тогда текущее значение ординаты

y_pi= y_pmax * ( p_i/p_max ),

где p_max= 4.4 МПа.

Масштаб индикаторной диаграммы

Площадь поршня

При построении графика силы, действующей на поршень, ординаты этого графика принимаем равными ординатам индикаторной диаграммы. Тогда масштаб силы

Для исследуемого механизма приведенный суммарной момент состоит из двух составляющих: движущей силы и момента сил сопротивления

.

Приведенный момент движущей силы определяется в текущем положении механизма по формуле

где F _дi - значение движущей силы,

где y_Fдi - ордината силы сопротивления,

μ _F - масштаб диаграммы сил.

V_qСi - значение передаточной функции в рассматриваемом положении механизма,

- угол между вектором силы и вектором скорости точки ее приложения.

Рис. 8.6

Масштаб диаграммы по оси абсцисс определяется по формуле

где b - база диаграммы ( отрезок оси абсцисс, который изображает цикл изменения обобщенной координаты).

1.3. Построение диаграммы приведенных моментов инерции I_v^пр = I _II^пр.

Инерционные характеристики звеньев механизма в его динамической модели представлены суммарным приведенным моментом инерции. При расчете эту характеристику динамической модели представляетсяв виде суммы двух составляющих переменной I_v^пр = I _II^пр и постоянной I_c^пр = I_I^пр. Первая определяется массами и моментами инерции звеньев, передаточные функции которых постоянны, вторые - массами и моментами инерции звеньев передаточные функции которых переменны.

Проведем расчет переменной части приведенного момента инерции I_v^пр = I _II^пр. Для рассматриваемого механизма во вторую группу звеньев входят звенья 2 и 3. Звено 3 совершает поступательное движение, звено 2 -плоское. Расчет переменной части приведенного момента проводится по следующим зависимостям:

где

Рис. 8.7

2. Построение диаграмм работы движущей силы, сил сопротивления и суммарной работы.

Диаграмму работы движущей силы получим интегрируя диаграмму ее приведенного момента

Интегрирование проведем графическим методом (рис.8.8), приняв при этом отрезок интегрирования равным k₁ . Тогда масштаб полученной диаграммы работы движущей силы будет равен

Рис. 8.8

Величина среднеинтегрального момента сил сопротивления определяется по формуле

3. Построение диаграмм кинетических энергий.

Диаграммы кинетических энергий для первой и второй групп звеньев получает на основании теоремы об изменении кинетической энергии системы

График кинетической энергии второй группы звеньев получим из зависимости

принимая, что w ₁ ≈ w _1ср . Тогда диаграмма приведенного момента инерции второй группы звеньев в масштабе рассчитанном по формуле

соответствует диаграмме кинетической энергии Т_II .

Рис. 8.9

График кинетической энергии первой группы звеньев приближенно строим по уравнению

В каждом положении механизма из ординат кривой A= f (φ₁) вычитаем ординаты y_TII и получаем ординаты искомой диаграммы T_I = f (φ₁). Для этого необходимо ординаты диаграммы T_II = f (φ₁) из масштаба μ_T перевести в масштаб μ_A^* по формуле

Диаграмма кинетической энергии первой группы звеньев представлена на рис. 8.9.

Рис. 8.10

4. Определение необходимого момента инерции маховых масс первой группы

Максимальное изменение кинетической энергии звеньев перD вой группы за цикл определяем по диаграмме

Тогда необходимый момент инерции маховых масс первой группы звеньев, обеспечивающий заданный коэффициент неравномерности, равен

4.1. Определение момента инерции дополнительной маховой массы.

В нашем случае момент инерции дополнительной маховой массы рассчитывается по следующей зависимости

,

где I₁₀ - момент инерции коленчатого вала .

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

Если считать, что w ₁ ≈ w _1ср , то

то есть диаграмма изменения кинетической энергии первой группы звеньев ∆T_I= f(φ₁) в другом масштабе соответствует диаграмме изменения угловой скорости ∆w_I= f(φ₁). Если считать что ординаты диаграмм равны, то

откуда

.

Ордината средней угловой скорости ( для определения положения начала координат на диаграмме угловой скорости )

После определения положения оси абсцисс на диаграмме угловой скорости можно определить начальное значение угловой скорости

а по ней кинетическую энергию механизма в начальном положении

6. Определение размеров маховика.

Принимаем конструктивное исполнение маховика - диск. Тогда его основные размеры и масса определятся по следующим зависимостям:

наружный диаметр

ширина b = y _b * D ,

масса m = 1230* D ³,

где r = 7.8 кг/дм³ - плотность материала маховика ,

ψ_b - коэффициент ширины .

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

Как отмечено ранее для расчета углового ускорения звена приведения ε₁ = f(φ₁ ) лучше пользоваться формулой :

Необходимые для расчета значения величин определяем по ранее построенным диаграммам. Диаграмма функции ε₁ = f(φ₁ ) приведена на рис. 8.11.

Рис. 8.11