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

книги / Основы прикладной теории упругих колебаний

..pdf
Скачиваний:
4
Добавлен:
13.11.2023
Размер:
23.13 Mб
Скачать

и находят 6 (v o ). Найденное значение 6 определяет абсциссу центра окружности (рис. 154).

Теперь можно провести малую дугу окружности из начальной точки фазовой траектории (г/о. So) по ходу часовой стрелки; та­ ким образом определяется первый элемент фазовой траектории. С чертежа можно снять новое значение фазовой координаты vi, вновь подставить ее в выражение (496) и найти 6(vi), т. е. определить на оси абсцисс положение нового центра окружности. При помощи это­ го центра строят второй элемент фа­

зовой траектории и т. д.

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

Пример 31. Построить фазовую кривую для процесса установления фрик­ ционных автоколебаний, описываемых дифференциальным уравнением (480), если сила трения дана в виде (486). Параметры системы: пг = 0,102 кГ сек/см,

с = 1000 кГ/см, R *

= 1 0

кг, v0 = 9,5 см/сек, у* = 10 см/сек.

По выражению

(486)

находим

ЗЯ* _

U

6Я*У0

з

vi

г г г

6я*

* о - =

V* -

и выражение (493) принимает вид

/ (0l) =

+

2 ,87v -f 2 ,80i>2 — 0 ,10y3 .

 

 

 

 

Определим частоту автоколебаний:

 

 

 

 

 

 

V

Т

100 сек

 

 

 

 

 

 

-

 

 

 

 

По выражению (496) находим

 

 

 

 

 

 

6 (v) =

— 0 ,029v — 2,80v2 - f 10v8

 

 

 

(497)

и приступаем к построению

фазовой

траектории, начиная с

любой

(заданной

начальными условиями) точки

фазовой плоскости. Пусть

х0=

0,

v0 =

0,09.

Тогда по формуле (497) 6 = — 0,018.

 

 

фазовой

тра­

На рис. 155, а показано

начало построения. Первый отрезок

ектории имеет центр в точке

с координатами 0,018; 0. Из этого центра прове­

дена дуга окружности из исходной точки Со до

точки Сь в которой v = 0,08.

По формуле (497) получится

б = — 0,015. Центр

второй дуги находится в точ­

302

ке с координатами 0,015; 0; из этого центра проведена вторая дуга до точки С2, в которой v = 0,07, и т. д.

'Фазовая траектория в целом изображена на рис. 155, б и обозначена циф­ рой II. Она представляет собой ceep-i

тывающуюся спираль. Другая фазовая траектория, начинающаяся в точке 0; 0,045, является развертывающейся спи­ ралью; она обозначена цифрой I. Фазо­ вые траектории типа I и II неограничен­ но приближаются к замкнутой траекто­ рии А, являющейся предельным циклом.

Кривая А несимметрична, причем особенно значительно нарушение симмет­ рии относительно вертикальной оси. Мак­ симальное и минимальное отклонения системы при ее движении по предельно­ му циклу равны соответственно 0,06 и 0,05 см. Таким образом, центр колебаний песколько смещен в направлении оси х

иполуразмах колебаний составляет

0,055

см.

 

Наибольшее

значение

v =

= 0,055

см

и

максимальная

скорость

Omax =

vp =

100 X 0,055

= 5.5

см/сек.

Эти

результаты

удовлетворительно

со­

гласуются

с

решением

(484),

согласно

которому

амплитуда автоколебаний

рав­

на а = 0,064

см и максимальная скорость

о та1 = ар =

6,4 см/сек.

В данном случае более точными сле­ дует считать результаты графо-аналити­ ческого решения при помощи дельта-ме­ тода; во всяком случае, оно свободно от

произвольного предположения о

гармо­

ническом характере процесса,

которое

было принято в аналитическом решении

энергетическим методом.

Рис. 155

30. РЕЛАКСАЦИОННЫЕ АВТОКОЛЕБАНИЯ

Автоколебания безмассовой системы

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

Ниже рассмотрен приближенный метод, относящийся к слу­ чаю весьма большой силы трения. Вернемся к схеме на рис. 149

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

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

имеет чисто статический вид

cx = R(v),

(498)

где сх — сила упругости пружины;

308

R{v) — сила трения, являющаяся функцией относительной скорости.

Нужно помнить, что соотношением (498) установлена про­ порциональность координаты х силе трения R\ поскольку коор­ дината х не может во времени изменяться скачком, то и сила трения должна оказаться непрерывной функцией времени t.

На рис. 156 точке О соответствует состояние равновесия гру­ за, при этом сила трения R0 уравновешивает силу упругости сх0

( * 0 — статическое отклонение груза). Убедимся, что это

состоя­

ние равновесия неустойчиво.

 

Пусть груз оттянут дополнительно на величину Ах некоторой

добавочной силой АР и удерживается в этом положении;

ввиду

неподвижности груза относительная скорость движения равна о,-» и состояние системы характеризуется той же точкой О.

 

 

Пусть

в

некоторое

 

 

мгновение,

которое будем

 

 

считать

начальным,

сила

 

 

АР исчезает

и груз

пре­

 

 

доставляется

самому

се­

 

 

бе. Так как начальная си­

 

 

ла

упругости

с(х0+Ах) =

 

 

— RQ -}- АР больше,

 

чем

 

 

Ro,

то

равенство

(498)

 

 

потребует

появления

си­

 

 

лы

трения R\,

большей

 

 

чем

 

Как видно из рис.

ние трения

 

156, мгновенное возраста­

от значения R0 до значения

Ry

возможно

только

вследствие

скачка относительной

скорости — уменьшения

ее до

величины Vy.

 

 

 

 

 

 

 

 

Началу колебаний соответствует точка 1. Конечно, мгновен­

ное изменение скорости является

лишь удобным

приблизитель­

ным описанием весьма быстрого изменения скорости. Уменьшение относительной скорости до значения Vy означает

появление положительной

абсолютной скорости

груза

х =

= v0Vy. Следовательно,

как только

оттянутый

вправо

груз

будет предоставлен самому

себе, он

немедленно

приобретает

направленную также вправо скорость х. Это вызовет дальнейшее перемещение груза вправо от положения равновесия. Такому процессу соответствует участок характеристики 1—2; с кинема­ тической стороны он характеризуется непрерывным ростом как

координаты х, так и скорости х, а со статической стороны — од­ новременным (и одинаковым) возрастанием как силы упругости пружины, так и силы трения. В точке 2 относительная скорость

обращается в нуль, т. е. абсолютная скорость х становится рав­ ной скорости ленты v0 и груз перестает скользить по ленте.

304

Продолжение движения вправо означало бы рост силы упру­ гости сх, и соответственно формуле (498) потребовалось бы уве­ личение силы трения R.

Однако сила трения возрастать далее не может; поэтому, как только изображающая точка достигает положения 2, движение груза вправо мгновенно прекращается.

Для определения нового значения скорости х нужно иметь в виду, что сила трения R должна остаться неизменной. Следова­ тельно, новое значение скорости отвечает изображающей точ­ ке 3.

Так как о3 > о0, то скорость груза х меняет знак. Следова­ тельно, после мгновенного изменения скорости от значения х2до

значения хз начнется движение груза влево. При этом движении деформация пружины будет убывать и, следовательно, будут убывать как сила упругости, так и сила трения R. Убывание си­ лы R возможно лишь вследствие уменьшения относительной ско­ рости v и соответствует участку 34; на этом участке уменьше­ ние относительной скорости будет непрерывным до точки 4.

Однако дальнейшее движение не может описываться ветвью 40— / —2, так как этой ветви соответствует рост силы трения R. Это противоречит уравнению (498), показывающему, что при движении груза влево сила упругости должна убывать, а не воз­ растать.

Поэтому в точке 4 скорость вновь должна мгновенно изме­

ниться— изображающая точка скачком переходит в

положе­

ние 5. Как и на участке 23, скачок 45 обозначает

разрыв

скорости при неизменном значении R. Затем груз будет двигать­ ся вправо вместе с лентой без скольжения, и сила упругости вновь начнет увеличиваться; при этом будет увеличиваться и си­ ла трения R. Это увеличение соответствует участку 52.

В точке 2 вновь произойдет разрыв скорости, изображающая точка переместится в положение 3, и т. д. Таким образом устано­ вится периодическое движение по циклу 234—5—23... без возврата на ветвь 40—12.

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

Весь автоколебательный цикл состоит из двух этапов (описа­ ние начинаем с точки 5).

П е р в ы й э тап. Движение груза с постоянной скоростью v0 вправо совместно с лентой (участок 52). Смещение х опре­ деляется через силу трения по формуле (498) и имеет следующие значения: в начале этапа (точка 5) х5 = R5 : с и в конце этапа (точка 2) х2 = R2 : с.

305

В т о р о й

э т а п .

Движение груза с убывающей скоростью

злево

(участок 34).

В начале этапа (точка 3) скорость Хз =

= — (оз — v0)

и смещение Хз = R$: с, причем х3 = х2. В конце

этапа

(точка 4) скорость х4 = — {v4— о0) и смещение х4 = R4: с

(■^4 =

 

 

На рис. 157, а и б представлены графики движения х = x(t)

и скорости x = x(t). График движения имеет «пилообразный» характер и резко отличается от закона гармонического движе­ ния. Наибольшие отклонения от значения Хо неодинаковы в обе стороны; поэтому удоб­ нее говорить не об ам­ плитуде, а о размахе автоколебаний; он оп­

ределяется суммой

I Х 2 |+ |Х ь |— /?2Rb

(499)

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

Длительность пер­ вого этапа движения вычисляется по зако­ нам равномерного дви­ жения

ti —

RzRb

(500)

 

cv0

 

Несколько сложнее Рис. 157 вычисление длительно­ сти второго этапа дви­

жения.

Рассмотрим производную

dR

dx

взятую по скорости движения х. Из уравнения (498) получим dR — cdx — cxdt.

Следовательно,

dR • dt

— - = cx — г- ; dx dx

отсюда находим длительность второго этапа движения

I *?.

R'.

(501)

U =

— dx.

J1

-X

 

Х%

306

Нижний предел интегрирования соответствует изображающей

точке 3, где хъ =

и0 — vz\верхний предел соответствует точке

для которой v4 =

Vov4.

составляет

 

 

Период автоколебаний

 

 

 

 

 

 

 

Т = /1

—[“to.

 

 

Пример 32. Как и в примере 31, характеристика трения имеет вид

Я =

ЗЯ* ( 1 — —

+

 

= 3Я„

1

v0— x

(Vb —x),3 1

 

 

3vd

 

 

 

 

3vc

 

 

 

и vо = 0,5

о*. Найти период и размах автоколебаний.

 

По характеристике трения находим

 

 

 

 

 

=

V 3

~

3,468у0;

Яд R o = 3Я*,'

 

 

 

 

= 2v0;

R4 = R6==Я*.

 

Пределы интегрирования в формуле (501)

 

 

 

*з = vo — v3 ~ — 2,468u0;

 

x4=

v0— o4 = — v0.

Длительность первого этапа движения по формуле

(500)

 

 

 

h=-3Я* — Я*

2Я*

 

 

 

 

 

 

 

cvn

 

cvо

 

 

Производная R по скорости х

 

 

 

 

 

 

 

 

 

 

\т +

 

X

Xа \

 

 

 

*

V*

 

°*

/

 

 

 

 

 

 

 

 

 

Длительность второго

этапа движения

по формуле

(501)

 

to= ЗЯ*

—»«

3

Л'

 

 

djc

1 ,023Я,

 

 

 

 

 

 

 

 

(4

 

 

 

 

 

 

— 2.468и0«

~1_ и*

„ 2

 

cvn

оказывается почти вдвое меньшей длительности первого этапа. Период автоколебаний

T= h + h~- 3,02Я*

cvn

Размах колебаний по формуле (499)

ЗЯ$ R^ 2Я»

сс

Автоколебания при упрощенной характеристике трения

Одновременный учет массы колеблющегося объекта и значи­ тельных сил трения оказывается в общем случае затруднитель­ ным. Но в некоторых случаях указанный учет необходим. Так, например, отмечено, что не всегда обеспечивается плавное дви-

307

<*N i
Рис. .158

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

Приближенное решение подобных вопросов возможно при помощи упрощенной характеристики трения, приведенной на рис. 158. Эта характеристика дает два значения силы трения:

максимальную силу трения по­ коя Ri и постоянную силу трения движения R2.

Рассмотрим движение мас­

тсивного груза 1 (рис. 150,6), свя­

занного пружиной 2 с ведущим звеном 3. Скорость движения по- Г следнего будем считать постоян­

ной и равной Do-

Пусть движение груза 1 и звена 3 совершается с общей ско­ ростью D0Сила упругости пружины Р равна силе трения /?2-

Если скорость Do весьма мала, то малое случайное препятствие может оказаться достаточным для остановки груза. Рассмот­ рим, что произойдет после этого.

Ведущее звено, продолжая движение вправо, будет сжимать пружину до тех пор, пока сила сжатия Р не сравняется с мак­ симальной силой трения покоя Ri. Так как после этого дальней­ ший рост силы трения невозможен, то произойдет «срыв» груза 1. При этом сила трения Ri мгновенно уменьшится до значения Rz, тогда как сила сжатия пружины Р = Ri мгновенно не может из­ мениться и в первое мгновение начавшегося движения будет попрежнему равна R\. Мгновение срыва примем за начало отсчета

времени (t — 0); при этом равны нулю

как смещение х, так и

скорость груза х, т. е.

 

х = 0; х = 0*.

(502)

Рассмотрим последующее движение груза. К текущему мгно­ вению длина пружины изменится на величину х — Dot и сила упругости пружины уменьшится до величины

Р (t) = Rt с(х v0t).

(503)

* В данном случае нельзя предполагать мгновенного скачка скорости. Такому скачку соответствует бесконечно большое ускорение, а следовательно, и бесконечно большие силы; здесь на груз в первое мгновение движения дей­ ствует конечная сила Ri R2.

308

Дифференциальное уравнение движения груза запишется в виде

Ri — с{х — v0t) — Rz = тх,

или

х-\-р*х = P W + - —

 

 

 

т

 

Решение этого уравнения, отвечающее начальным услови­

ям (502), имет вид

 

 

 

х = v0t — — sin pt +

 

(1 — cos pt)\

(504)

P

c

 

 

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

Скорость движения груза меняется по закону

 

х = v0v0cos pt + -p^ ~ — sin pt.

(505)

c

 

Следующая остановка груза произойдет в мгновение, когда х

вновь обратится в нуль. Условие остановки х = 0 приводит со­ гласно соотношению (505) к уравнению

v0— »о cos рк +

^ sin ptx = 0,

 

С

где t\— время от момента срыва до новой остановки. Введем безразмерный параметр

а =

= _JL A f

(506)

cv0

pv0

 

где Af — разность коэффициентов трения покоя и движения. Условие остановки принимает вид

a sin ptx = cos ptx— 1.

(507)

Решая это трансцендентное уравнение, находим

sin ptx= —

cos ptx= 1 —аа

(508)

 

1 + а2 ;

1 + а3"

 

Получив отсюда значение t\, можно по формуле (504) опре­ делить абсциссу Х\ груза в момент новой остановки, т. е. путь, пройденный грузом за время tx:

хх = vQtx— — sinptx+ --1 ~ ~

(1 — cosptx) = v0tx +

p

c

p

С учетом выражений (508) найдем по формуле (503) силу сжатия пружины в момент остановки:

P(tx) = 2Rz- R x.

309

Так как # 2 < то P(t\) < R\; следовательно, после останов­ ки груз некоторое время будет оставаться на месте, пока сила упругости пружины вновь не достигнет значения предельной силы трения покоя Ri.

За время, в течение которого груз покоится, сила сжатия пру­ жины постепенно возрастает на величину

bP = R1- P ( t 1) = 2(R1- R z),

и соответствующее укорочение пружины составит д / = А Р = 2(RX- R 2)

сс

Этой же величине равен путь, который проходит ведущее зве­ но за время, пока груз стоит на месте. Следовательно, длитель­ ность состояния покоя груза равна

А/ _

2 ( —R2) _

2g

(509)

Щ

CV0

р-

 

Тот же результат можно найти из условия t'o - г к ) =

выражающего равенство перемещений груза и ведущего звена за период.

Таким образом, период автоколебаний груза определяется формулой

Т — к + к,

для пользования которой нужно сначала найти t\ из выражения (508), а затем к — из формулы (509). В момент t = Т происхо­ дит следующий срыв груза и начинается новый цикл авто­ колебаний.

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

sin pty -s- 0; t1 — — .

P

При учете выражения (509) период автоколебаний прибли­ женно равняется

Т = я + 2а

Р

Значение второго слагаемого в числителе увеличивается с уменьшением скорости iv, характер движения при двух различ­ ных малых значениях v0 показан на рис. 159. С уменьшением

310

скорости vQрастет период автоколебаний Т и одновременно уве­ личивается значение Х \ .

о

t o

t

Рис. 159

Пример 33. Определить, при каком значении параметра а длительность остановки груза равна длительности его движения. Найти соответствующую* скорость v0, если р — 50 сек-1 и Д/ = 0,05.

Согласно условию t2 = U', поэтому, обозначив

pk = Рк = Р,

имеем из выражений (508) и (509)

sin 23 = —

*

 

 

 

1-Нра

 

 

Наименьший отличный от нуля корень этого трансцендентного уравнения

Р = 2,03.

 

соответствующую

среднюю

Теперь из формулы (506) можно получить

скорость движения

 

 

 

 

v0= —— Д/ = 0,48 см/сек.

 

 

ар

 

 

 

 

Если средняя скорость движения и0 <

0,48 см/сек, то интервалы

остановок

будут короче интервалов движения.

 

 

собственная часто­

Как видно, чем жестче система (т. е. чем больше ее

та р), тем меньшей оказывается скорость

о0, отвечающая

условию

задачи.