Динамика машин и механизмов в установившемся режиме движения
.pdf
Как отмечено ранее, приведенный момент сил MП представляется в виде алгебраической суммы
M П M ПД M ПС .
Определение M ПС выполняется из условия равенства мгновенных мощностей
M СП 1 Fi Vi M i i .
Откуда
|
|
M СП F X x i/ FY y i/ M i i i1 sign 1 |
, (1.12) |
|
где F |
X |
и F |
– проекции силы F на оси координат; |
|
|
Y |
i |
|
|
x i/ |
и y i/ – проекции аналога скорости точки |
приложения |
||
силы Fi ;
ii 1 – передаточная функция от i-го звена, к которому приложен момент Mi , к звену 1;
sign 1 1 при направлении вращения звена 1 против часо-
вой стрелки;
sign 1 1 - при направлении вращения звена 1 по часовой
стрелке.
В формуле M ПС (1.12) силы FХ , FY и моменты Mi берутся со
знаками, соответствующими правой системе координат (положительное направление вращения – против часовой стрелки).
Так, для горизонтального механизма (рис. 1.8, а) M ПС определяется из равенства
M СП 1 F ПС VB G 2 V S 2 G 3 VB ,
20
откуда
M СП F ПС x /B G 2 y /S 2 G 3 y /B sign 1 .
Рис. 1.8
Учитывая, что xB/ i31 , yB/ 0, sign 1 1, получим
M ПС FПС i 31 G 2 y /S 2 .
В рассматриваемом положении силаFПC имеет отрицательное
значение, так как она направлена против положительного направления оси X.
Для вертикального механизма (рис. 1.8, б) аналогичным образом можно получить
M СП FПС yB/ G2 yS/ 2 G3 yB/ sign 1 FПС i31 G2 yS/ 2 G3 i31 .
Сила FПC в изображенном случае положительна.
Приведенный момент движущих сил M ПД определяется из усло-
вия, что при установившемся режиме движения изменение кинетической энергии машины за цикл равно нулю, т.е.
21
Т AДЦ АСЦ 0,
откуда за цикл АД Ц АС Ц.
Работа сил сопротивления вычисляется как
i
АС МПС d 1.
0
Интегрирование выполняется численным методом по правилу трапеций:
|
AСi |
AС i 1 |
M ПСi M ПС |
i 1 |
|
|
1 |
|
, |
|
|||
|
|
|
|
||||||||||
|
2 |
|
|
|
|
|
|||||||
где 1 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
– шаг интегрирования в радианах. |
1 |
|
, |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
n |
||
|
|
|
|
|
|
|
|
|
|
|
|
||
n – число интервалов деления одного оборота кривошипа.
Сучетом AДЦ МПД Ц при Ц 2
МД АСЦ .
П2
1.2.3.Определение переменной составляющей приведенного момента инерции IП//
Переменная составляющая I П// определяется из условия равен-
ства кинетических энергий, т.е. кинетическая энергия звена приведения, имеющего момент инерции I П// , равна сумме кинетических
энергий звеньев, характеризуемых переменными передаточными функциями:
I // |
2 |
|
I S i i2 |
|
m iVS2i |
.. |
П |
1 |
|
|
|
|
|
|
2 |
2 |
2 |
|||
|
|
|
|
22
Разделив это выражение на 2 , с учетом того, что |
V 2 x |
2 |
y |
2 |
, |
||
1 |
|
|
S i |
S i |
|
S i |
|
получим |
|
|
|
|
|
|
|
I П// I S i i i21 m i x S/ |
i 2 y S/ |
i |
2 . |
|
|
|
|
Для звеньев 2, 3 кривошипно-ползунного механизма (рис. 1.8) |
|
|
|||||
I П// m 2 x S/ 2 2 y S/ 2 2 I S 2i 212 m 3 i 312 . |
|
|
|
|
|||
Производная dIП , |
необходимая в последующем для определе- |
|||||
|
|
d 1 |
|
|
|
|
ния закона движения звена приведения, имеет вид |
|
|||||
d I П |
2 m2 xS/ |
2 x S// |
2 yS/ |
2 y S// |
2 I S 2 i21 i21/ m3 i31 i31/ |
sign 1 . |
d |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
1.2.4.Определение постоянной составляющей приведенного момента инерции IП/
имомента инерции маховика IМ
Воснову расчета положен метод Н.И Мерцалова 4 . Для определения изменения кинетической энергии машины Т предвари-
тельно определяем работу движущих сил АД . Для i-го положения
|
АД i M ПД 1i , |
||||
где |
1i |
|
1 |
|
i 1 . |
|
|
||||
Тогда |
T i A Д i AСi . |
||||
|
|||||
23
Изменение кинетической энергии T/ звеньев с постоянным приведенным моментом инерции I П/ равно
T/ i Ti T// i ,
где T// i – кинетическая энергия звеньев, создающих переменную
составляющую I П// . По методу Н.И. Мерцалова, T// i определяется приближенно по средней угловой скорости 1cр :
|
|
I // |
2 |
|
T |
|
Пi |
1с р |
. |
|
|
|||
// i |
|
|
2 |
|
|
|
|
|
|
Далее из полученного за цикл массива значений находим максимальную T/ a и минимальную T/ b
пользуя которые, вычисляем максимальный перепад энергии:
T/ ab T/ a T/ b .
T/ (рис. 1.9)
величины, искинетической
Рис. 1.9
24
Тогда необходимая величина I П/ , при которой имеет место вра-
щение звена приведения с заданным коэффициентом неравномерности , равна
I / |
|
T |
/ a b |
, |
где 1ср |
n1 |
. |
(1.13) |
|
30 |
|||||||
|
|
|
||||||
П |
|
12cр |
|
|
|
|||
Момент инерции маховика определяется как
I М I П/ I П0 ,
где IП0 – приведенный момент инерции всех вращающихся масс машины (ротора двигателя, зубчатых колес, кривошипа).
Иногда величина IП0 может оказаться больше полученного значения I П/ . Это означает, что не требуется установки маховика. Ре-
альный коэффициент неравномерности вращения в этом случае из
(1.13) равен
|
T/ ab |
. |
||
I 0 |
2 |
|||
|
|
|||
|
П |
1cр |
|
|
1.2.5. Определение закона вращения звена приведения
С помощью зависимости T/ ( 1 ) , используемой при определе-
нии постоянной составляющей приведенного момента инерции I П/
по методу Мерцалова, можно получить зависимость угловой скорости звена приведения 1( 1 ) .
Из рис. 1.9 видно, что для любого положения кинетическая энергия звеньев, обладающих постоянным приведенным моментом
инерции I П/ , равна
T/ i T/ ср Т/ ср Т/ i ,
25
где T/ ср |
I П/ |
12ср |
, |
Т/ ср |
|
Т/ а |
Т/ b |
. |
||||
|
2 |
|
|
|
2 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
I / |
|
2 |
|
|
|
|
|
|
|
Так как T/ i |
П |
|
1i |
, то текущее значение угловой скорости |
||||||||
|
|
2 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1i |
|
|
2T/i |
sign 1 . |
|
|
|
|
|
|
|
|
|
/ |
||||
|
|
|
|
|
|
|
|
|
|
IП |
|
|
Угловое ускорение 1 определяется из дифференциального уравнения движения звена приведения:
|
M ПД МПСi |
21i d IПi |
|
|||
|
|
|
|
|
||
1i |
2 d 1 |
sign 1 . |
||||
|
||||||
IП/ IП// i |
||||||
|
|
|||||
1.2.6. Схема алгоритма программы исследования динамической нагруженности
машинного агрегата
Рассмотренные в предыдущих параграфах материалы позволяют разработать программу исследования динамической нагруженности машинного агрегата. В качестве объекта исследования взята технологическая машина, в которой основным исполнительным механизмом является кривошипно-ползунный механизм (например, го- ризонтально-ковочная машина). Примерная схема алгоритма такой программы приведена на рис. 1.10.
Осуществляется ввод исходных данных (блок 1). Пример подготовки исходных данных показан в табл. 1.1. Следует обратить внимание на соответствие направления вращения кривошипа 1cр , зна-
ка FПС по отношению к положительному направлению соот-
ветствующей оси координат, а также на знак величины эксцентриситета е.
26
|
Начало |
|
|
||
Исходные данные |
|
|
|||
1 , xB MAX ( yB MAX ) |
|||||
|
|
1 |
0 |
|
|
|
AC1 0 |
|
|
||
|
|
i 1, n |
|
|
|
Кинематические харак- |
|||||
теристики |
|
|
|||
M C |
, |
I // |
, dI Пi ,T |
|
|
Пi |
|
Пi |
d |
// i |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
i 1 |
|
нет |
|
|
|
|
|
||
да |
|
|
|
||
|
|
ACi |
|
|
|
1 |
1 1 |
|
|||
|
|
|
A |
|
|
Рис. 1.10 |
|
|
|||
27
A
Окончание рис. 1.10
M Д
П
i 1, n
AД i , Ti , T/ i
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Поиск |
|
|
|
T |
|
, T |
|
|
|
|
|
||||||
|
|
|
/ a |
/ b |
|
|
|
|
максималь- |
|||||||
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ного и мини- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
мального |
T |
, I / |
, I |
M |
T |
/ cр |
, T |
/ ср |
элементов |
||||||||
|
||||||||||||||||
|
/ ab |
П |
|
|
|
|
|
|
||||||||
i 1, n |
|
||
T |
, |
, |
1i |
/ i |
1i |
|
|
Печать
результатов
Конец
Окончание рис. 1.10
28
|
|
|
|
|
Таблица 1.1 |
|
|
|
|
|
|
|
|
№ |
Параметр |
Условное |
Единица |
|
Величина |
|
пп |
|
обозначение |
измерений |
|
|
|
1 |
2 |
|
3 |
4 |
|
5 |
1 |
Схема кривошипно- |
|
- |
- |
|
|
|
ползунного механизма |
|
|
|
|
|
2 |
Размеры звеньев |
l1 lOA |
м |
|
0,0742 |
|
|
|
|
||||
|
|
l3 |
l AS 2 |
м |
|
0,0741 |
|
|
l2 |
lAB |
м |
|
0,2225 |
|
|
|
e |
м |
|
0,01335 |
3 |
Начальная обобщенная |
|
о |
град |
|
2,58 |
|
координата |
|
|
|
|
|
4 |
Массы и моменты |
|
|
|
|
|
|
инерции звеньев |
|
m2 |
кг |
|
400 |
|
|
|
|
|||
|
|
|
m3 |
кг |
|
500 |
5 |
Сила полезного |
|
I S 2 |
кгм2 |
|
8,35 |
|
сопротивления FПС |
FПС1 |
H |
|
0 |
|
|
|
|
||||
|
|
FПС2 |
H |
|
0 |
|
|
|
FПС3 |
H |
|
0 |
|
|
|
FПС4 |
H |
|
0 |
|
|
|
FПС5 |
H |
|
0 |
|
|
|
FПС6 |
H |
|
0 |
|
|
|
FПС7 |
H |
|
0 |
|
|
|
FПС8 |
H |
|
0 |
|
|
|
FПС9 |
H |
|
0 |
|
|
|
FПС10 |
H |
|
0 |
|
|
|
FПС11 |
H |
|
-12321 |
|
|
|
FПС12 |
H |
|
-27142 |
|
|
|
FПС13 |
H |
|
-125000 |
|
|
|
|
|
|
|
|
29
