
книги из ГПНТБ / Сергеев В.И. Исследование динамики плоских механизмов с зазорами
.pdfопределяются моменты разрыва кинематической цепи и восстанов ления контакта, величины скоростей относительного движения элементов соударяющихся пар при восстановлении контакта. Подробно исследуется влияние зазора на динамику кривошипноползунного механизма при различных соотношениях параметров последнего.
Проводится исследование влияния трения в кинематической паре с зазором на величину реакции в паре и характер движения пальца в подшипнике. Показано, что наличие трения вносит значи тельные изменения в динамику рассматриваемого механизма с за зором. Так, кривая зависимости величины реакции в шатунном подшипнике с зазором от угла поворота кривошипа на участках безотрывного движения имеет вид затухающих колебаний, и с тече нием времени значение реакции оказывается равным соответствую щему значению в идеальном механизме, чего не наблюдается при решении поставленной задачи без учета сил трения. Резкие измене ния величины реакции, предшествующие моментам разрыва кине матической цепи, идентичны как в механизме без учета сил трения, так и при учете последних в паре с зазором.
Согласно проведенным расчетам характер изменения величины реакции в условиях безотрывного движения в качественном отноше нии находится в удовлетворительном соответствии с частными ре зультатами, полученными другими авторами при исследовании аналогичных по своей конструкции механизмов иными приближен ными расчетными и экспериментальными методами [33, 50, 54, 59, 60, 61, 64, 66, 72].
Проведен анализ точности решения уравнений движения. Выяв лено, что при решении рассматриваемой задачи следует задавать точность интегрирования уравнений движения порядка ІО-4. Мень шая точность приводит к значительным искажениям траекторий движения и графиков изменения величины реакции в шатунном подшипнике с зазором. В ряде случаев при необходимости получить лишь качественную картину движения механизмов с зазорами с целью сокращения непроизводительных затрат машинного времени можно решать поставленную задачу с точностью порядка 10-2.
Далее в работе проводится исследование влияния на динамику рассматриваемого механизма величины зазора, внешней нагрузки, отношения длины кривошипа к длине шатуна и соотношений масс шатуна и ползуна, причем параметры механизма изменяются в пре делах, соответствующих широкому диапазону реально существую щих механизмов. Расчеты проводились как с учетом различного вида сил трения в паре кривошип — шатун, так и при отсутствии указанной силы трения. Результаты расчетов при разных значениях параметров приведены в виде таблиц.
Следует отметить, что предложенный подход к исследованию динамических моделей механизмов с зазорами не накладывает никаких ограничений на величины зазоров и геометрические раз меры звеньев механизма.
11
f л а в а !
УРАВНЕНИЯ ДВИЖЕНИЯ И ОБЩИЙ МОДЕЛИРУЮЩИЙ АЛГОРИТМ
ИССЛЕДОВАНИЯ МЕХАНИЗМОВ С ЗАЗОРАМИ
Влияние зазоров в кинематических парах на динамику механизмов
Будем называть «идеальным» механизмом такой, в кинематических парах которого отсутствуют зазоры. Если при этом механизм об ладает одной степенью свободы, то положение любой точки его звеньев однозначно определяется заданием одной координаты ве дущего звена. При заданных параметрах механизма и внешних си лах, действующих на него, движение ведущего звена будет опи сываться одним обыкновенным дифференциальным уравнением второго порядка относительно выбранной координаты. Подобные уравнения могут быть получены различными методами на основа нии общих законов механики [1, 4, 7, 25, 30]. Зная закон движения ведущего звена, его скорость и ускорение в каждый момент вре мени, можно, исходя из кинематических соотношений для рассмат риваемого механизма, получить выражения для скорости и ускоре ния любой точки звеньев механизма. Величины реакций в его ки нематических парах могут быть определены из последующего кинетостатического анализа.
Однако «реальный» механизмАтличается от идеального как не точностью изготовления отдельных его звеньев, так и наличием зазоров в кинематических парах. Условия собираемости механизма и возможности относительного движения его звеньев приводят к возникновению зазоров в кинематических парах. Величины зазо ров выбираются в соответствии с предъявляемыми к механизму технологическими и точностными требованиями и ограничиваются принятыми допусками. По назначению зазоров механизмы можно разделить на две группы:
1) механизмы, в которых зазоры служат только для обеспече ния собираемости механизмов и относительного движения звеньев, их по возможности стремятся уменьшить;
2) виброударного действия, принцип работы которых основан на использовании эффектов соударения смежных звеньев в зазоре.
Настоящая работа посвящена исследованию влияния зазоров на динамику механизмов первой группы.
Наличие зазоров в кинематических парах реальных механизмов намного усложняет анализ движения последних, а также определе-
12
йиё |
реакций |
в кинематических |
парах |
с зазорами [10, 40]. Так, |
|
при |
наличии |
зазора только |
в одной |
низшей |
кинематической |
паре появляются одна или две |
дополнительные |
степени свободы. |
На определенных участках движения такого механизма будет осу ществляться контакт в кинематической паре с зазором, на других участках может происходить независимое движение одного звена относительно другого, при этом будет иметь место разрыв кине матической цепи. В течение одного цикла движения ведущего зве на контактное и бесконтактное движения могут многократно чере доваться. Движение механизма с контактом и без контакта описы вается разными уравнениями, решения которых должны сопря гаться в точках перехода от одного вида движения к другому. Сле довательно, в рассматриваемом случае приходится проводить ис следование динамики механизмов с переменной структурой.
При чередовании замкнутого и разомкнутого состояний кине матической цепи механизма с зазорами будет происходить много кратное соударение элементов кинематической пары, при этом уве личивается износ трущихся поверхностей по сравнению с предпо лагаемым износом в идеальном механизме [19]. Кроме того, как показывают расчеты, даже в условиях сохранения контакта в кине матической паре скорость относительного движения элементов па ры с зазором на отдельных участках намного превышает скорость относительного движения этих же элементов в идеальном механиз ме, а величина реакции в этих соединениях также превосходит соответствующие значения для идеального механизма. Все это при водит к ускоренному износу элементов кинематических пар и поте ре работоспособности механизма вследствие возникновения в его рабочем цикле «постепенных» отказов [15, 16].
Характерной особенностью задачи исследования динамики ме ханизмов с зазорами является несоизмеримая малость величин последних по сравнению с номинальными размерами отдельных звеньев механизма. Однако наличие этих зазоров может оказать существенное влияние как на законы движения его звеньев, так и на работоспособность механизма в целом.
Динамическая модель механизма с зазорами
Не ограничивая общности исследования механизмов с зазорами (с точки зрения разработки общих методов подхода к исследованию подобных механизмов), рассмотрим аксиальный кривошипно-пол зунный механизм с зазорами в соединениях кривошип—шатун и шатун — ползун. Схематически этот механизм показан на рис. 1 (зазоры изображены в увеличенном масштабе). Величина зазора равна разности радиусов подшипника и шипа (А = гх —г2). Подоб ная схема механизма используется в двигателях внутреннего сго рания, различного вида насосах, механизмах плоско-печатных стан ков и во многих других машинах, служащих для преобразования вращательного движения в поступательное или наоборот.
13
Будем считать заданным движение кривошипа, т. е. зададим |
|
со = а (і) |
в виде известной функции времени и положим, в частнос |
ти, со = |
const. Для простоты рассуждений рассмотрим шатун в ви |
де прямолинейного тонкого стержня с равномерно распределенной массой. Это позволяет при выводе уравнений движения считать массу шатуна сосредоточенной в его геометрическом центре (в точ ке s). Массу ползуна будем считать сосредоточенной в его центре масс (в точке 03). Звенья механизма будем предполагать абсолют но жесткими, зазор между ползуном и направляющими, а также трение в этой паре учитывать не будем.
Предполагается, что зависимость сил трения в парах 1 — 2 и 2—3 от скорости относительного движения соответствующего шипа в подшипнике у имеет форму полинома
FTP = — R (/гтрі sign т + &тр,Т + |
&тр,Т2), |
|
|
(1) |
|||||
где |
R — нормальная |
составляющая |
реакции |
в данном |
соедине |
||||
нии; |
Ьр. |
(і = 1,2,3) — соответственно |
коэффициенты |
сухого, |
|||||
жидкостного и квадратичного трения; |
у |
— угловая скорость дви |
|||||||
жения точки контакта по поверхности подшипника, |
которая на |
||||||||
ходится из |
решения |
уравнений |
движения |
звеньев |
механизма. |
Вдальнейшем при численных расчетах на ЭЦВМ будем задавать
висходной информации различные значения коэффициентов тре
14
ния, в частности равные нулю, исследуя таким образом динамику рассматриваемого механизма в условиях, когда трение в парах с зазорами отсутствует либо имеется только сухое, жидкостное или квадратичное трение, либо их комбинация.
Выберем в качестве неподвижной системы координат прямо угольные декартовы координаты ХОУ, начало отсчета поместим
в точку О, относительно которой происходит вращение кривошипа, |
|
а ось X направим вдоль прямой, по которой перемещается ползун. |
|
Введем в рассмотрение две подвижные системы координат: х10 1у1 |
|
с началом в центре |
подшипника кривошипа и x 3Osy3 с началом в |
центре подшипника |
ползуна (см. рис. 1). Оси подвижных систем |
координат |
направим параллельно осям неподвижной системы |
|
( О л I Оз-ѵ'зі! ОХ, О^ЦОзУзЦ OY). |
При выводе уравнений движения |
|
с контактом |
вместо декартовых |
систем координат х101у1 и х 30 3у3 |
воспользуемся полярными (р^ ух) и (р3, у3), что позволяет записать
уравнения движения |
в более простом |
виде, при |
этом |
рх = |
= |
|
= const И Рз = |
Д3 = |
const. |
механизм, |
т. е. |
механизм |
|
Идеальный |
кривошипно-ползунный |
без зазоров, имеет одну степень свободы, положение всех звеньев такого механизма однозначно определяется заданием угла а. При наличии двух зазоров в соединениях 1—2 и 2—3 этот механизм получает дополнительно 2, 3 или 4 степени свободы в зависимости от того, сохраняется ли контакт в кинематических парах механиз мов или кинематическая цепь разомкнута. Возможны четыре вида движения рассматриваемого механизма:
I. С сохранением контакта во всех кинематических парах, в этом случае дополнительное движение механизма определяется
координатами Уі и Уз- |
паре 2—3 |
и разрывом в |
паре |
|
II. С сохранением контакта в |
||||
1— 2. |
В этом случае дополнительное движение определяется тремя |
|||
координатами у3 и хх, ух. |
паре 1—2 |
и разрывом в |
паре |
|
III. С сохранением контакта в |
||||
2— 3. |
Тогда дополнительное движение будет |
определяться |
также |
|
тремя |
координатами ух и х3, у 3. |
|
|
|
IV. С разрывом в парах 1—2 и 2—3. В этом случае дополни |
тельное движение будет определяться четырьмя координатами хх,
Уи х3 и у 3. |
• |
|
|
|
При переходе от одного вида движения к другому (от контакт |
||||
ного к бесконтактному |
или |
наоборот) необходимо |
вычислить |
на |
чальные условия для |
соответствующих уравнений |
движения |
и |
|
одновременно осуществить |
преобразование декартовых координат |
в полярные или наоборот. Например, при выполнении условия на рушения контакта в паре 1—2, т. е. при обращении нормальной составляющей реакции Rx в нуль начальные условия для уравне ний свободного движения в поле зазора А* определяются по форму
лам: |
|
|
дс° = Дх cos TJ, |
У\ = Лі sin T“, |
(2), (3) |
15
где Yi° — значение угла ух, при котором произошел разрыв кине матической цепи в паре 1—2. Дифференцируя равенства (2) и (3) по времени, получим два других начальных условия:
*і° = — ДіТ?sin TS. |
Уі = Аіт! cos -г®, |
(4), (5) |
где fi° —значение угловой скорости дополнительного движения пальца шатуна по поверхности подшипника кривошипа в момент отрыва.
При переходе от свободного движения в поле зазора к контакт ному начальное значение угловой координаты у-, определяется по формуле
хт
|
|
|
arccos -А- |
при у№ > |
О |
|
|
|
|
|
|
|
|
||||
т і « = |
|
|
1 |
, (ft) |
|
|
|
|
|
|
|
|
|
|
(6) |
||
|
|
|
2я — arccos-А- |
при |
|
О, |
|
|
|
|
|
|
|
||||
где |
х\ , |
|
у\ |
— координаты |
точки |
обоймы |
подшипника, в |
кото |
|||||||||
рой произошло восстановление контакта в паре 1—2. |
|
|
|||||||||||||||
Начальное |
значение |
угловой |
скорости |
у ^ |
вычисляется |
ис |
|||||||||||
ходя |
из |
принятой модели |
удара. |
В общем |
случае |
ti° = f |
(4É)> |
||||||||||
y[k), |
x[k\ |
y[k)). Воспользуемся формулами (2) |
и |
(3) |
перехода от |
||||||||||||
полярных |
координат |
(рд, уг) |
к декартовым |
(лу, гу) |
для движения |
||||||||||||
внутри зазора |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
* 1 — Pi cos Ті, |
Уі = Pi sin Yi- |
|
|
|
|
|
|
|
|
||||||||
Дифференцируя их по времени, получим |
|
|
|
|
|
|
|||||||||||
М = |
рх cos Ті — PiTi sin Ті, |
|
Уі = |
Pi sin Yx + |
pxTi cos ту. |
|
|||||||||||
Решая |
эту |
систему |
алгебраических |
уравнений |
относительно |
||||||||||||
f 1 и |
рх и учитывая, что sin ух = |
гу/Рі. cos ух = |
ху/Рі> |
получим |
|||||||||||||
Ті = |
(УЛ — ХіУі)ІрІ, |
|
|
PI = |
(ВД. + |
УіУіУрі, |
|
|
|
|
где рх — скорость движения пальца шатуна по нормали к окруж ности.
В настоящей работе принята модель неупругого удара, т. е. коэффициент восстановления нормальной составляющей скорости соударения принимается равным нулю, следовательно,
ТІ4) = ( 4 М Й,- Э Д / А?- |
(7) |
'S. |
уравнении |
■'Аналощчно вычисляются начальные условия для |
свободной) движения и движения с контактом в зазоре А3. Момент перехода от свободного движения -в поле соответствующего зазора
к |
контактному определяется нарушением неравенств х\ + у\ < |
< |
А? или ХІ-+УІ < Ад. |
}6
Вывод уравнений движения кривошипно-ползунного механизма с зазорами
Поскольку в рассматриваемой модели механизма с зазорами движение ведущего звена считается заданным, то для получения уравнений добавочного движения механизма достаточно записать уравнения относительно координат xlt уІУх 3, у 3или ух, УзВосполь зуемся уравнениями Лагранжа второго рода:
где Т — кинетическая энергия ползуна и шатуна, выраженная че
рез обобщенные координаты q£\ Q,- — обобщенные силы, |
соответ |
||
ствующие выбранным обобщенным координатам. |
|
|
|
в |
В качестве обобщенных координат при сохранении |
контакта |
|
кинематических парах выберем соответственно углы |
ух и |
у3, |
|
а |
при разрыве кинематической цепи — декартовы координаты |
хъ |
|
Уи |
и Уз- |
|
|
Определение кинетической энергии шатуна и ползуна. Кине тическая энергия шатуна в неподвижной системе координат XOY выражается формулой
Тш = тгѵ\і2 -f / sß72, |
(9) |
где m2 — масса шатуна; vs — скорость поступательного движения его центра масс в неподвижной системе координат XOY\ Is — мо мент инерции шатуна относительно его центра масс; ß —угловая
скорость |
вращения шатуна. |
|
в |
виде |
|||
Квадрат скорости |
представим |
||||||
ѵі = Xi + Y l |
|
|
|
|
|
|
|
где X s и Ys — проекции скорости |
vs на оси неподвижной системы |
||||||
координат XOY. Для того |
чтобы их определить, запишем коорди |
||||||
наты центра масс шатуна |
в системе XOY следующим образом: |
||||||
X s = |
г cos а + |
хг + Y cos ß, |
|
(10) |
|||
Ys = |
-isin ß + |
г/g. |
|
|
|
(11) |
|
Дифференцируя равенства (10) и (11) по времени, получим |
|||||||
X s= — reosin а -1- хх — у |
ß sin ß, |
(12) |
|||||
Ys = ~ ßcosß + |
y3. |
|
|
|
(13) |
||
Определим теперь |
ß. Для |
этого |
выразим длину отрезка, А К |
(ем. рис. 1) через геометрические параметрьГмехан^змшнякоориина-
I иаучно-трхн.-1'."-»**»аг! Г
б библио '..-к« ',У. С Р I J-
L |
ЭКЗЕМПЛЯР |
I |
ЧИТАЛЬНОГО ЗАЛА |
ты точек А и В в подвижных системах координат. Тогда |
|
|
|||
AK = г sin а -|- уа = / sin ß -|- у3. |
|
(14) |
|||
Дифференцируя последнее равенство по времени, можно записать |
|||||
reocos а + уг = |
/ß cos ß + y3. |
|
(15) |
||
Откуда |
|
|
|
|
|
ß = (reocos а + |
ух — уз)//cos ß. |
|
(16) |
||
Для упрощения дальнейших выкладок введем обозначения |
|
|
|||
U = |
— reosin а + хъ |
V = reo cos а + уг — у3. |
(17), |
(18) |
|
Соответственно в полярных координатах |
|
|
|||
U = — reo sin а — AxXxSin |
|
(19) |
|||
V = |
reocos а -f |
ДіУх cos ух — Д3у3cos у3. |
|
(20) |
|
Тогда вместо (16) |
имеем |
|
|
|
|
ß = |
17//cosß |
|
|
|
(21) |
и |
|
|
|
|
|
Xs = |
f/ -■ £ tgß, |
Y , = ^ + y3. |
|
|
|
Возводя в квадрат обе части этих равенств и складывая полу |
|||||
ченные выражения, будем иметь |
|
|
|||
о* = |
XI + У| = |
Д 2 - |
UV tg ß + V4t (tg2ß + 1) + y3 (V -f y3)- |
|
Поскольку в рассматриваемой модели механизма шатун пред полагается прямолинейным тонким стержнем, его момент инерции
относительно |
центра |
масс, расположенного в геометрическом цен |
|||
тре шатуна, |
равен |
/ 5 = т 2/2/ 12. |
|
|
|
Подставляя найденные значения ß,u2 и Is в формулу для кине |
|||||
тической энергии шатуна (9), получим |
|
|
|||
т — ИІІ |
|
V2 |
ys (V + у,) + |
V2 - |
|
1ш — |
2 |
и * + U V tg ß + т (tg2ß + 1) + |
1 2 cos2 ß. * |
Ползун совершает только возвратно-поступательное движение, следовательно, его кинетическая энергия может быть записана в виде Тп = т3и3/2, где т3 — масса ползуна, ѵ3 — скорость его цен тра масс в неподвижной системе координат XOY.
Продифференцируем по времени следующее выражение, при помощи которого определяется координата центра масс ползуна
внеподвижной системе XOY:
Х3 — г cos а -|- д'х -|- I cos ß — х3.
18
Врезультате получим
Х3 =: ѵ3 — — л о sin а -|- х1 — д'з — 1$sin ß.
|
С |
учетом |
обозначений |
(17), |
(18) и |
(21) можно |
записать |
ѵ3 == |
|||
= |
U — V tg ß — л'3. Тогда |
выражение |
для |
кинетической энергии |
|||||||
ползуна принимает вид |
|
|
|
|
|
|
|
||||
|
Т п = f |
( t / - y t g ß - x 3)2. |
|
|
|
|
|
||||
|
Общая кинетическая энергия звеньев в 2 и 3 определяется сле |
||||||||||
дующим образом: |
|
|
|
|
|
|
|
||||
|
т = Тш+ Тп = f Гт/ 2 — UV tg ß + ^ ( t g 2ß + !) + |
|
|||||||||
|
|
|
|
|
|
1/2 |
|
V i g V - x 3)\ |
(22) |
||
|
|
|
+ |
У&(У + Уз) |
|
12 cos2 ß . + % ( u |
|||||
|
Для записи кинетической энергии в полярной системе коорди |
||||||||||
нат |
следует |
воспользоваться |
преобразованиями, |
аналогичными |
|||||||
(4) |
и |
(5) |
|
|
|
|
|
|
|
|
|
|
*з = — ДзТзsin Тз. |
Уз = |
ДзТз cos т3. |
|
|
|
|||||
|
Тогда формулу (22) можно переписать в виде |
|
|
||||||||
|
Т = f |
[£/» - UV tg ß + |
£ |
(tg2 ß + |
1) + |
ДзТз cos тз (V + |
|
||||
|
|
+ |
ДзГз С05Гз) + 1 2 |
^ |
р + у |
(u — V tg ß + Д 3Тз sin Тз)2- |
Определение величины реакции в кинематических парах с зазо рами. Контакт внутренней поверхности подшипника с поверхностью шипа осуществляется по их общей образующей. Силу реакции в кинематической паре с зазором можно представить как сумму двух составляющих: нормальной реакции, совпадающей по направ лению с общей нормалью к соприкасающимся окружностям, и силы трения, направленной по касательной к указанным окружностям (см. рис. 1). За положительное направление нормальной составля ющей реакции примем направление к центру подшипника. Для нахождения нормальной составляющей реакции в паре 1—2 сог ласно принципу Даламбера приравняем нулю сумму моментов сил, действующих на шатун относительно точки В (см. рис. 1) с учетом сил инерции шатуна. За положительное направление момента сил примем направление возрастания угла ß. Тогда в соответствии с рис. 1 будем иметь
= — Rxl sin %+ |
F'^1 sin Ti + |
Mn + |
РИХу sin ß + |
|
|
|
+ Рад у |
cos ß — P2y co sß = |
0, |
(23) |
|
где M „ — момент сил |
инерции шатуна |
относительно |
точки |
В\ |
Р„X и Рщ — проекции сил инерции шатуна на оси X и Y, прило-
19
&eüHi>te к его Центру масс; Р 2 — сила веса Шатуна. Согласно прин ципу Даламбера
= Р’ Рн.ѵ' — m2x s, Рни — пцУsi
где ß — угловое ускорение шатуна; X s и |
Ks — проекции ускоре |
ния центра масс шатуна соответственно на |
оси X и У. |
Воспользовавшись формулой (1) для силы трения и решая (23) относительно Rlt можно записать выражение для нормальной сос
тавляющей силы реакции в паре 1—2 в виде |
|
|||
Ri = |
т-2 |
X |
|
|
( k Tpi sign Ti + |
*TpJ i + b TpJ |
|
||
2 [sin X + |
2) sin T-J] |
|
||
|
Ф |
-Xssinß + |
Kscosß + geos ß ) , |
(24) |
|
X (-fT + |
|||
где g — ускорение |
силы тяжести, а углы т и т* определяются |
из |
геометрических соотношений для рассматриваемого механизма.
Согласно |
рис. |
1, т = |
+ ß; |
= (л/2) — т. |
|
Вторые производные ß,X's, ^определяются следующим образом. |
|||||
Дифференцируя равенства (12) и (13) по времени |
и учитывая, что |
||||
d = © = |
const, |
получим |
|
|
|
Xs = |
г©2cos а + Xi — Y (ß s>n ß + ß2cos ß), |
(25) |
|||
Ys = |
4- (ß cos ß — ß2sinß)-Ly3. |
|
(26) |
Для нахождения ß продифференцируем равенство (15) по вре мени:
—г©2sin а + |
г/і = /ß cos ß — /ß'2sin ß -f- y3. |
|
Откуда |
|
|
- r a r sin а + (/) — г/з + /ß2 sin ß |
(27) |
|
ß = |
I cos ß |
|
Таким образом, оказывается, что нормальная составляющая реакции Rx в паре кривошип — шатун зависит от геометрических параметров механизма и от первых и вторых производных по вре мени от добавочных координат в кинематических парах 1—2 и 2—3. Последние следует определить из решения дифференциальных уравнений движения механизма.
Аналогичным образом, приравнивая нулю сумму моментов сил,
действующих на шатун относительно его центра масс, |
и принимая |
|
во внимание принцип Даламбера, можно определить |
нормальную |
|
составляющую реакции R 3 в паре шатун — ползун: |
|
|
R s = sin (Тз + ß) [ # i sin (Ті + ß) + |
ß — Frp cos (Ti + ß) + |
|
|
+ Frp cos (Тз + ß)]- |
20