Техническая механика часть 1
.pdf
МFtОBFOBcos.
пспсn псn
Подставив соотношения (8.48) и (8.49) в формулу (8.46) для к.п.д., получим связь между моментами Мдв и Мпс
Мпс = Мдв u1n , |
(8.50) |
где u1n = 1 / n – передаточное отношение механизма. Соотношение (8.50) также носит общий характер, в силу чего его называют универсальным законом передачи моментов.
8.9. Интегрирование дифференциального уравнения движения механизма
Как уже отмечалось, одной из задач анализа движения механизмов является определение истинного закона движения его звеньев, в том числе выходного звена. Для механизмов с одной степенью подвижности достаточно знать закон движения звена, которое выбрано в качестве звена приведения. В самом деле, в силу требований эквивалентности реального механизма и его динамической модели в любой момент времени кинетическая энергия звена приведения равна кинетической энергии всего механизма, а его угловая скорость совпадает с угловой скоростью начального звена. Как следует из выражений (8.7) и (8.14), изменение кинетической энергии может происходить как за счет изменения угловой или линейной скорости, так и за счет изменения приведенного момента инерции или приведенной массы. Те и другие величины зависят от различных факторов. Если первые определяются характером технологического процесса, то вторые зависят от передаточных отношений в механизме и, следовательно, являются функциями обобщенной координаты. Поэтому «синхронное» изменение сил и приведенного момента инерции или приведенной массы, приводящее к постоянству скорости звена приведения в пределах одного цикла работы механизма, невозможно. Таким образом, даже при установившемся режиме движение звеньев механизма, в том числе и начального, не является равномерным.
Закон движения звена приведения (а значит, и начального звена) может быть установлен в результате решения уравнения (8.24) или уравнения (8.26) в зависимости от того, какая динамическая модель принята. Рассмотрим возможность решения этих уравнений на примере уравнения (8.24):
191
M п |
|
2 |
dI |
п |
I |
п |
2 |
dI |
n |
In |
d |
(8.51) |
2 |
|
|
2 |
|
|
dt |
||||||
|
|
d |
|
d |
|
|||||||
Неизвестной функцией, подлежащей определению, является угловая скорость звена приведения ω(t) как функция времени. В самом деле, если явный вид функции ω(t) известен, то закон движения φ(t) может быть найден из соотношения, связывающего угловую скорость и угол поворота звена приведения φ:
(t) |
d |
|
dt . |
(8.52) |
Приведенный момент инерции Iп и приведенный момент сил Мп, входящие в уравнение движения (8.51), определяются с помощью соотношений (8.9) и (8.12) соответственно. В общем случае они могут быть функциями обобщенной координаты φ, а приведенный момент сил может зависеть также и от времени.
Дифференциальное уравнение (8.51) является нелинейным относительно неизвестной функции. Это обстоятельство сильно затрудняет его интегрирование с помощью аналитических методов, которые пригодны только для некоторых частных случаев. Так, если приведенный момент Мп зависит только от обобщенной координаты, то уравнение (8.51) может быть записано в виде:
|
|
|
2 |
|
d |
|
d |
|
In |
2 |
|
||
M п |
|
|
dIn |
In |
|
|
|
|
(8.53) |
||||
|
|
|
|
|
|
|
|||||||
|
|
d |
d |
|
|
|
|
. |
|||||
|
|
2 |
|
|
d |
2 |
|
|
|
||||
Здесь использовано несложное преобразование:
d |
|
d d |
|
d |
|||
|
|
|
|
|
. |
||
dt |
d |
dt |
d |
||||
Проинтегрируем обе части уравнения (8.53) от некоторого начального значения обобщенной координаты φ0 до произвольного значения . В результате получим:
|
|
In ( ) |
2 |
2 |
|
||
|
Мп ( )d |
|
|
|
In ( 0 ) 0 |
, |
|
|
|
|
|
||||
0 |
2 |
|
2 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
192 |
|
||
где 0 – значение угловой скорости звена приведения в том его положении, которому соответствует значение обобщенной координаты равное φ0. Из последнего равенства может быть найдено значение угловой скорости при произвольном значении :
|
2 |
|
( ) |
||
|
||
|
||
In ( ) |
||
|
|
|
|
In ( 0 ) |
0.5 |
|
||
Mn ( )d |
02 . |
(8.54) |
|||
In ( ) |
|||||
|
0 |
|
|
||
|
|
|
|
||
Полученное выражение задает угловую скорость звена приведения как функцию обобщенной координаты. Чтобы найти зависимость угловой скорости от времени воспользуемся соотношением (8.52). Из него следует, что дифференциалы времени и обобщенной координаты связаны между собой следующим образом:
|
d |
|
t |
|
d |
|
|
dt |
или |
dt |
. |
||||
|
|
||||||
|
( ) |
||||||
( ) |
|||||||
|
|
t0 |
0 |
|
|||
|
|
|
|
|
|||
Здесь t0 и t – моменты времени, в которые положение звена приведения определяется значениями обобщенной координаты φ0 и соответственно. Интегрируя обе части последнего равенства, нетрудно получить время движения механизма как функцию угла поворота звена приведения: t = t( ). Исключая из зависимостей ( ) и t( ) угол поворота , в итоге получим явный вид функции (t).
Изложенная процедура интегрирования уравнения движения (8.24) предполагает, что зависимости Мп( ) и Iп( ) известны. На практике значения приведенного момента сил и приведенного момента инерции определяются для нескольких значений обобщенной координаты. Полученные значения либо аппроксимируются некоторой функциональной зависимостью (как правило, степенной), либо представляются в виде таблицы. В первом случае все интегралы, входящие в приведенные соотношения, точно вычисляются, а во втором – они вычисляются приближенно с помощью численных методов.
Уравнение движения (8.24) может быть проинтегрировано в квадратурах еще в некоторых частных случаях. Однако чаще всего оно может быть решено только численными методами. Зная явный вид зависимости угловой скорости звена приведения от времени можно оценить степень неравномерности начального звена в пределах одного цикла работы механизма. Она определяется коэффициентом неравномерности движения :
193
|
max min |
, |
(8.55) |
|
ср |
||||
|
||||
|
|
|
где max и min – соответственно максимальное и минимальное значения угловой скорости начального звена за цикл, а ср – средняя за цикл его скорость. Она приближенно определяется как среднее арифметическое максимального и минимального значений:
ср |
max min |
. |
(8.56) |
|
2 |
||||
|
|
|
Неравномерность движения звеньев механизма сильно сказывается на условиях его работы. При высокой степени неравномерности большие ускорения вызывают появление значительных сил инерции, которые приводят к возрастанию динамических нагрузок, как на элементы кинематических пар, так и на опоры машины или технологического аппарата. Как следствие, могут возникнуть нежелательные переменные воздействия и вибрации на фундамент. Для каждого типа машин практикой устанавливаются допускаемые значения коэффициента неравномерности δ. Например, для грохотов и конвейеров = 0,04 0,1; для электрогенераторов = 0,003 0,01; для металлорежущего оборудования = 0,01 0,02.
Основным способом снижения коэффициента неравномерности движения начального звена механизма является увеличение инерционности начального звена, что приводит к увеличению сопротивления изменению его скорости. Практически это достигается добавлением к этому звену значительной постоянной массы, выполняемой в виде вращающегося вместе со звеном массивного маховика. Чем больше масса маховика (а, следовательно, и приведенный момент инерции механизма), тем равномернее движется начальное звено.
8.10. Подбор маховика и определение угловой скорости ведущего звена механизма с маховиком
Конструктивно маховик представляет собой, как правило, сплошной диск со ступицей или шкив с тяжелым массивным ободом, спицами и ступицей. Для маховика-диска зависимость между массой маховика mМ и моментом инерции IМ выражается формулой
194
|
m D2 |
|
|
|
IM |
M |
, |
(8.57) |
|
8 |
||||
|
|
|
а для маховика-шкива эта связь определяется соотношением:
|
|
m (D2 |
D2 ) |
|
|
IM |
M н |
в |
, |
(8.58) |
|
8 |
|
||||
|
|
|
|
|
где mМ – масса диска или обода; D – диаметр диска; Dн и Dв наружный и внутренний диаметры обода (кольца). Обычно Dн=1,1Dср, Dв=0,9Dср (Dср– средний диаметр). Моментами инерции ступицы и спиц можно пренебречь, т.к. они составляют менее 1% от суммарного момента инерции маховика.
Момент инерции маховика, обеспечивающего заданное значение коэффициента неравномерности хода δ, часто определяется с помощью графоаналитического метода энергомасс. Метод энергомасс сравнительно прост и в то же время позволяет обеспечить достаточную для инженерных расчетов точность. Его существо состоит в построении и дальнейшем анализе диаграммы энергомасс. Последняя представляет собой график зависимости приращения кинетической энергии Т от приведенного момента инерции механизма IП. В соответствии с соотношением (8.4) приращение кинетической энергии, очевидно, равно разности работ движущих сил АПдв и сил полезного сопротивления АПс, выполненных к данному моменту времени в пределах цикла установившегося движения. Различие в законах изменения приведенных работ движущих сил АПдв(φ) и сил сопротивления АПс(φ), а также сложная нелинейная зависимость приведенного момента инерции от обобщенной координаты внутри цикла движения механизма являются непосредственными причинами неравномерности хода.
Диаграмма энергомасс Т = f(IП ) при установившемся режиме движении представляет собой замкнутую петлеобразную кривую (рис. 8.8), построенную по ряду точек (обычно 8 – 12). Исходными данными для построения диаграммы являются законы изменения (в пределах одного цикла (φ=0 ÷ 2π)), приведенного момента сил сопротивления МПс = f(φ) (8.12) и приведенного момента движущих сил МПдв = f(φ), который зависит от типа двигателя и должен быть задан. Интегрирование данных зависимостей по обобщенной координате позволяет определить функции работ указанных моментов: АПс = f(φ), АПдв = f(φ). При получении зависимости АПдв=f(φ) следует учитывать, что за целое число циклов избыточная работа Аизб = АПдв – АПс = 0, т.е.
195
ΔT, T T, T
|
A i |
|
|
n i |
|
OΔT |
|
A |
|
|
|
|
T 0 |
I П, I |
|
|
|
|
β i |
|
OT |
m i |
|
|
I П, I |
Рисунок 8.8
в начале и в конце цикла значения АПдв и АПс совпадают. Вычисляя разность между АПдв и АПс для нескольких φ определяют Аизб = Т =
f(φ). Зависимости |
Т = f(φ) и IП = f(φ) позволяют построить диаграмму |
энергомасс, т.е. |
Т = f(IП ). |
Рассмотрим |
свойства диаграммы энергомасс (рис. 8.8). Если |
предположить, что для начального положения механизма угловая скорость ω0 ведущего звена известна, то нетрудно определить и величину полной кинетической энергии в данный момент времени:
T |
I |
2 |
|
|
|
П0 0 |
. |
(8.59) |
|
|
|
|||
0 |
|
2 |
||
|
|
|
|
Величину Т0 в соответствующем масштабе отложим от начала координат (точка О Т на рис. 8.8) вниз и получим точку ОТ. Из этой точки проведем луч в произвольную точку диаграммы Аi. Тангенс угла наклона этого луча к оси IП равен отношению длин отрезков:
A m
tg i Oi mi . (8.60)
T i
Из рис. 8.8 видно, что текущие значения кинетической энергии и приведенного момента инерции могут быть выражены длинами соответствующих отрезков графика:
196
T |
I |
Пi |
2 |
|
|
|
|
|
|
||
|
|
|
|
||||||||
|
i |
|
A m |
T ; |
(8.61) |
||||||
|
|
|
|
||||||||
i |
|
2 |
|
|
|
i i |
|
||||
|
|
|
|
|
|
|
|
|
|||
IПi |
|
|
OT mi |
|
I , |
|
(8.62) |
||||
|
|
|
|||||||||
где μТ – масштабный коэффициент по оси Т (размерность коэффициента μТ кгм2с-2/ мм = Нм/мм); μI - масштабный коэффициент по оси IП (его размерность кгм2/мм). После подстановки двух последних соотношений в формулу (8.60) для tg βi получим:
tg i |
|
T |
I |
|
|
I |
|
I |
Пi |
|
2 |
|
|
I |
i2 . |
|
i |
|
|
i |
|
(8.63) |
|||||||||||
T I Пi |
T |
|
|
|
2 T |
|||||||||||
|
|
|
|
2I Пi |
|
|
|
|
||||||||
Таким образом, тангенс угла наклона луча, проведенного на диаграмме Т = f(IП) из начала отсчета ОТ в произвольную точку Аi, оказывается пропорциональным квадрату угловой скорости ωi
ведущего звена механизма. Это основное свойство диаграммы энергомасс.
Из формулы (8.63) вытекает выражение для угловой скорости ведущего звена:
i |
|
|
2 T tg i |
|
. |
(8.64) |
|
||||||
|
|
|
I |
|
||
Максимальная угловая скорость ωmax будет соответствовать углу βmax (луч, проведенный из точки ОТ, касается верхней точки диаграммы энергомасс), а минимальная скорость ωmin - углу βmin (луч, проведенный из точки ОТ, касается нижней точки диаграммы). Подсчет величины δ по формуле (8.55) и сопоставление ее с заданным значением коэффициента неравномерности чаще всего указывает на необходимость использования в машине дополнительной (маховой) массы. В этом случае начало системы координат (Т, IП) помещают в точке ОМ, которая смещена вдоль оси IП влево на величину IМ (рис. 8.9). Точка ОМ находится на пересечении лучей, проведенных под углами βmin и βmax и охватывающих диаграмму энергомасс. Углы βmax и βmin рассчитываются с помощью соотношения (8.63):
197
tg max |
|
|
I |
|
max2 |
|
|
I |
|
ср2 (1 ) , |
||
|
2 T |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
2 T |
||||
tg min |
|
|
I |
min2 |
|
|
I |
ср2 (1 ) . |
||||
2 T |
2 T |
|||||||||||
|
|
|
|
|
|
|
||||||
при подстановке в него значений ωmax и ωmin , найденных в свою очередь по формулам:
max |
|
|
|
|
min |
|
|
|
cр 1 |
|
и |
cр 1 |
. |
||||
|
|
|
2 |
|
|
|
|
2 |
ΔT, T
|
A i |
|
q |
A p |
|
OΔT |
|
|
k |
I П, I |
|
T, T |
||
|
T 0
ℓβ i
βmax |
βmin |
|
|
|
|
OМ |
|
|
I М |
m |
I П, I |
|
|
Рисунок 8.9 |
Таким образом, точка ОМ является началом координат диаграммы энергомасс, при котором коэффициент неравномерности движения механизма и средняя угловая скорость имеют заданные значения.
Расстояние ОМm от нового начала координат до прежней оси ординат Т определяет искомое значение момента инерции маховика:
IM = |OM m| μI . |
(8.65) |
198
Может оказаться, что точка ОМ будет лежать на оси Т. Это означает, что при заданном значении средней угловой скорости ведущего звена механизм в состоянии обеспечить заданное значение коэффициента неравномерности движения без маховика. При графическом определении положения точки ОМ могут встретиться два крайних случая:
а) точка ОМ оказывается за пределами чертежа, а углы βmax и βmin близки к нулю. В этом случае ее положение находят аналитически
следующим образом. Из треугольников ОМmk |
и OMml (рис. 8.9) |
следует: |
|
km OM m tg max , |
(8.66) |
lm OM m tg min . |
(8.67) |
Отсюда
kl = km – lm = OM m ( tg βmax – tg βmin ).
Для величины отрезка OM m получим:
OM m |
|
kl |
|
|
|
|
. |
(8.68) |
|
tg max |
|
|||
|
tg min |
|
||
Этой формулой можно пользоваться во всех случаях, когда возможно графическое определение длины отрезка kl.
б) точка ОМ лежит за пределами чертежа, а углы βmax и βmin близки к 90 0. В этом случае охватывающие лучи могут не пересекать ось ординат Т в пределах чертежа. Но тогда они пересекают ось абсцисс IП в точках p и q (рис. 8.9). Из треугольников O Т kq и O Т lp следует:
O Т k = O Т q tg βmax ,
O Т l = O Т p tg βmin .
Следовательно,
kl = O Т l – O Т k = O Т p tg βmin – O Т q tg βmax .
Подставив полученное соотношение в формулу (4.23), окончательно получим:
199
OM m |
O Т ptg min |
O Т qtg max |
. |
(8.69) |
tg max |
|
|||
|
tg min |
|
||
Длину отрезка OM m, выраженную в мм, называют чертежным
значением момента инерции маховика, и обозначать через I M . Значение момента инерции в его обычной размерности (кг м 2) получается путем умножения чертежного значения на величину масштабного коэффициента μI. Чертежное значение среднего запаса
кинетической энергии T0 может быть получено по формуле (см. рис. 8.9):
|
|
|
|
0 O T k km O T k |
|
M tg max |
|
|||||
|
T |
I |
(8.70а) |
|||||||||
или по формуле |
|
|||||||||||
|
|
|
|
|
|
|
|
|
||||
|
T0 |
(O T q I |
M )tg max . |
(8.70б) |
||||||||
|
|
|
|
|
||||||||
Значение Т |
в Дж определяется путем умножения T |
на масштабный |
||||||||||
0 |
0 |
|
||||||||||
коэффициент μТ.
При известных координатах точки ОМ нетрудно найти углы βi (номера положений: i = 0, 1, …, n), характеризующие угловую скорость ведущего звена в соответствующих положениях:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
tg i |
|
T0 Ti |
|
. |
(8.71) |
||||||
|
|
|
|
|
|
|
|
|
|||||
I M I |
Пi |
||||||||||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|||||||||
Здесь Ti |
– чертежное значение |
приращения |
запаса кинетической |
||||||||||
энергии, взятое с учетом знака (знак «+», если точка Ai лежит выше оси IП, знак «-» – в противном случае); I Пi – чертежное значение
приведенного момента инерции механизма без учета маховой массы в соответствующий момент времени.
Далее по формуле (8.64) находятся истинные значения угловой скорости ведущего звена в промежуточных точках внутри отдельного цикла установившегося движения. Закон движения ведущего звена (изменение углового ускорения в форме ε=f(φ )) может быть получен с помощью дифференцирования зависимости ω=f(φ ) и последующего
200
