double arrow

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

Решение задачи регулирования хода машины по методу Н.И.Мерцалова.

При расчете маховика (или решении задачи регулирования хода машины) по методу Н.И.Мерцалова задача решается в следующей последовательности:

  • Определяются параметры динамической модели, например для ДВС Мпрд - приведенный суммарный момент движущих сил и IпрII - приведенный момент инерции второй группы звеньев.
  • Определяется работа движущих сил Ад интегрированием функции Мпрд = f(j 1) за цикл движения машины (допустим 2p);
  • Определяется работа движущих сил за цикл и приравнивается к работе сил сопротивления Адц = Асц. Из этого равенства определяется среднеинтегральное значение момента сил сопротивления

Мпрсср = Асц/ (2p);

и для него строится диаграмма работы Ас = f(j 1). Суммированием этой диаграммы и диаграммы Ад = f(j 1) получаем диаграмму А = f(j 1).

  • Делается допущение w 1» w 1ср, при котором TII » IпрII *w 1ср2/ 2 (первое допущение метода Мерцалова), и определяется TII = f(j 1).
  • Определяется кинетическая энергия первой группы звеньев

TI = А -TII + Tнач = А-TII + TIнач + TIIнач .

Так как начальные значения кинетической энергии неизвестны, то если учесть, что Tнач = TIнач + TIIнач, D TI = TI - TIнач, D TII = TII - TIIнач , получим

D TI = А- D TII ,

то есть, вычитая из суммарной работы приращение кинетической энергии второй группы, получим приращение кинетической энергии первой группы.

По функции D TI = f(j 1) определяется максимальное изменение кинетиской энергии за цикл D TImax. Второй раз делаем допущение w 1» w 1ср на основании которого, как показано выше, можно записать

IпрI = D TImax / (d *w 1ср2).

Из этого выражения, определив предварительно D TImax, можно решить две задачи:

  • задачу синтеза - при заданном [d ] определить необходимый для его обеспечения приведенный момент инерции IпрI нб ,
  • задачу анализа - при заданном IпрI определить обеспечиваемый им коэффициент неравномерности d.

Решение этой задачи рассмотрим на конкретном примере машинного агрегата привода буровой установки.

Дано: Кинематическая схема машины - lAB = 0.12м, lBC = 0.528м, lBS2 = 0.169м, средняя частота вращения кривошипа - w 1ср = 47.124 рад/с2, массы звеньев -

m2 = 24.2 кг, m3 = 36.2 кг, момент инерции - I 2S = 1.21 кг* м2, I 10 = 2.72 кг* м2, максимальное давление в цилиндре - pmax = 4.4 МПа, коэффициент неравномерности вращения [d ] = 1/80, индикаторная диаграмма (приведена на рис. 8.3).

_________________________________________________________________

Определить: закон движения машины w 1 = f(j 1) и e 1 = f(j 1), момент инерции маховика Iдоп , обеспечивающий заданную неравномерность вращения [d ].

  1. Определение параметров динамической модели: Мпрд - приведенного суммарного момента движущих сил и IпрII - приведенного момента инерции второй группы звеньев.
  2. Определение первых кинематических передаточных функций. Определение кинематических передаточных функций для звеньев механизма u21 = u31 , центров масс VqS1 , VqS2 и VqS3 и точки приложения движущей силы VqD . Для определения этих функций воспользуемся методом проекций векторного контура механизма.

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

l AB + l CB = l AC ; l AS2 = l AB + l BS2 .

Для первого векторного контура l AB+ l CB = l AC проекции на оси координат

lAB * cos j 1 + lCB * cos j 2 = xC = 0,

lAB * sin j 1 + lCB * sin j 2 = yC = SC,

j 2 = arccos (- lAB * cos j 1 / lBC ).

Рис. 8.4

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

- lAB * sin j 1 - lCB * u21* sin j 2 = 0,

lAB * cos j 1 + lCB * u21* cos j 2 = VqC ,

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

u21 = - lAB * sin j 1 / (lCB* sin j 2 ),

VqC = lAB * cos j 1 + lCB * u21* cos j 2 .

Для третьего векторного контура l AS2 = l AB + l BS2 проекции на оси координат

xS2 = lAB * cos j 1 + lBS2 * cos j 2 ,

yS2 = lAB * sin j 1 + lBS2 * sin j 2 .

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

VqS2x = - lAB * sin j 1 - lBS2 * u21* sin j 2 ,

VqS2y = lAB * cos j 1 + lBS2 * u21* cos j 2,

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

Рис. 8.5

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

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

ypi= ypmax * (pi/pmax ),

где pmax= 4.4 МПа.

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

m p = ypmax /pmax.

Площадь поршня Sп = p *dп2 /4.

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

m F = m p/Sп.

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

Mпр = Mпрд + Mпрс.

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

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

F дi = yFдi / m F,

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

m F - масштаб диаграммы сил.

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

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

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

m j = b / 2*p,

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

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

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

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

Ivпр = I IIпр = I2Впр + I2Ппр+ I3пр,

где

I2Ппр = m 2 * VqS22, I2Впр = IS2 * u212, I3пр= m3* VqС2,

Рис. 8.7

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

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

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

tg y 1 = yD / xD j 1 = yMпрд1/ k1 D Aд * m А/ (D j 1 * m j) = Mпр д1 * m М / k1

так как D Aд / D j 1 = Mпр д1, то m А / m j = m М / k1, откуда

m А = m М *m j / k1 .

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

Мпрсср = Асц/ (2p).

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

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

D Т = Т - Тнач, A = D Т I + D Т II .

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

Т II = III пр*w 1ср2 /2,

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

yI = yT IпрII * m I = (IпрII * w 1ср2 / 2) * m T, откуда

m T = 2* m I /w 1ср2,

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

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

ТI = Т - ТII.

В каждом положении механизма из ординат кривой A= f (j 1) вычитаем ординаты yTII и получаем ординаты искомой диаграммы TI = f (j 1). Для этого необходимо ординаты диаграммы TII = f (j 1) из масштаба m T перевести в масштаб m A* по формуле

yTII* = yTII * m A*/ m T.

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

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

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

D TImax = (y D TImax )/ m A.

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

IIпр = D TImax / (w 1ср2 * [d ] ).

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

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

Iдоп = IIпр - I10 ,

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

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

Если считать, что w 1» w 1ср, то

D TI = IIпр*w 1ср * D w 1,

то есть диаграмма изменения кинетической энергии первой группы звеньев D TI= f(j 1) в другом масштабе соответствует диаграмме изменения угловой скорости D w 1= f (j 1). Если считать что ординаты диаграмм равны, то

yD w 1 = yD TI m A* D TI = m w * D w 1 m A* IIпр*w 1ср * D w 1 = m w * D w 1 ,

откуда

m w = m A* IIпр *w 1ср .

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

yw 1ср = w 1ср *m w.

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

w 10 = yw 10 /m w,

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

TI нач = IIпр *w 1ср2 /2.

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

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

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

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

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

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

y b - коэффициент ширины.

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

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

e 1 = dw 1/dt = М пр/ Iпр- w 12/(2* Iпр) * (d Iпр /dj 1).

Необходимые для расчета значения величин определяем по ранее построенным диаграммам. Диаграмма функции e 1 = f(j 1 ) приведена на рис. 8.10.

Рис. 8.10


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями: