Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Сергеев В.И. Исследование динамики плоских механизмов с зазорами

.pdf
Скачиваний:
11
Добавлен:
23.10.2023
Размер:
4.26 Mб
Скачать

определяются моменты разрыва кинематической цепи и восстанов­ ления контакта, величины скоростей относительного движения элементов соударяющихся пар при восстановлении контакта. Подробно исследуется влияние зазора на динамику кривошипноползунного механизма при различных соотношениях параметров последнего.

Проводится исследование влияния трения в кинематической паре с зазором на величину реакции в паре и характер движения пальца в подшипнике. Показано, что наличие трения вносит значи­ тельные изменения в динамику рассматриваемого механизма с за­ зором. Так, кривая зависимости величины реакции в шатунном подшипнике с зазором от угла поворота кривошипа на участках безотрывного движения имеет вид затухающих колебаний, и с тече­ нием времени значение реакции оказывается равным соответствую­ щему значению в идеальном механизме, чего не наблюдается при решении поставленной задачи без учета сил трения. Резкие измене­ ния величины реакции, предшествующие моментам разрыва кине­ матической цепи, идентичны как в механизме без учета сил трения, так и при учете последних в паре с зазором.

Согласно проведенным расчетам характер изменения величины реакции в условиях безотрывного движения в качественном отноше­ нии находится в удовлетворительном соответствии с частными ре­ зультатами, полученными другими авторами при исследовании аналогичных по своей конструкции механизмов иными приближен­ ными расчетными и экспериментальными методами [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 и 23 от скорости относительного движения соответствующего шипа в подшипнике у имеет форму полинома

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.

механизм,

т. е.

механизм

Идеальный

кривошипно-ползунный

без зазоров, имеет одну степень свободы, положение всех звеньев такого механизма однозначно определяется заданием угла а. При наличии двух зазоров в соединениях 12 и 23 этот механизм получает дополнительно 2, 3 или 4 степени свободы в зависимости от того, сохраняется ли контакт в кинематических парах механиз­ мов или кинематическая цепь разомкнута. Возможны четыре вида движения рассматриваемого механизма:

I. С сохранением контакта во всех кинематических парах, в этом случае дополнительное движение механизма определяется

координатами Уі и Уз-

паре 23

и разрывом в

паре

II. С сохранением контакта в

12.

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

координатами у3 и хх, ух.

паре 12

и разрывом в

паре

III. С сохранением контакта в

23.

Тогда дополнительное движение будет

определяться

также

тремя

координатами ух и х3, у 3.

 

 

 

IV. С разрывом в парах 12 и 23. В этом случае дополни­

тельное движение будет определяться четырьмя координатами хх,

Уи х3 и у 3.

 

 

 

При переходе от одного вида движения к другому (от контакт­

ного к бесконтактному

или

наоборот) необходимо

вычислить

на­

чальные условия для

соответствующих уравнений

движения

и

одновременно осуществить

преобразование декартовых координат

в полярные или наоборот. Например, при выполнении условия на­ рушения контакта в паре 12, т. е. при обращении нормальной составляющей реакции Rx в нуль начальные условия для уравне­ ний свободного движения в поле зазора А* определяются по форму­

лам:

 

 

дс° = Дх cos TJ,

У\ = Лі sin T“,

(2), (3)

15

где Yi° — значение угла ух, при котором произошел разрыв кине­ матической цепи в паре 12. Дифференцируя равенства (2) и (3) по времени, получим два других начальных условия:

*і° = — ДіТ?sin TS.

Уі = Аіт! cos -г®,

(4), (5)

где fi° —значение угловой скорости дополнительного движения пальца шатуна по поверхности подшипника кривошипа в момент отрыва.

При переходе от свободного движения в поле зазора к контакт­ ному начальное значение угловой координаты у-, определяется по формуле

хт

 

 

 

arccos -А-

при у№ >

О

 

 

 

 

 

 

 

 

т і « =

 

 

1

, (ft)

 

 

 

 

 

 

 

 

 

 

(6)

 

 

 

2я — arccos-А-

при

 

О,

 

 

 

 

 

 

 

где

х\ ,

 

у\

— координаты

точки

обоймы

подшипника, в

кото­

рой произошло восстановление контакта в паре 12.

 

 

Начальное

значение

угловой

скорости

у ^

вычисляется

ис­

ходя

из

принятой модели

удара.

В общем

случае

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). За положительное направление нормальной составля­ ющей реакции примем направление к центру подшипника. Для нахождения нормальной составляющей реакции в паре 12 сог­ ласно принципу Даламбера приравняем нулю сумму моментов сил, действующих на шатун относительно точки В (см. рис. 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 и 23. Последние следует определить из решения дифференциальных уравнений движения механизма.

Аналогичным образом, приравнивая нулю сумму моментов сил,

действующих на шатун относительно его центра масс,

и принимая

во внимание принцип Даламбера, можно определить

нормальную

составляющую реакции R 3 в паре шатун — ползун:

 

R s = sin (Тз + ß) [ # i sin (Ті + ß) +

ß — Frp cos (Ti + ß) +

 

+ Frp cos (Тз + ß)]-

20

Соседние файлы в папке книги из ГПНТБ