книги / Основы прикладной теории упругих колебаний
..pdfТогда
и
i
j Ndx
_о_____
i
J mx2dx
о
Если продольные силы N появляются в результате вращения стержня с постоянным распределением массы (рис. 26), то
q = т&2х',
Л /= | ? & = |
ПИй" |
(/2 — х2); |
|
|||
|
2 |
|
|
|||
j* Ndx = |
ffl(D2/3 |
|
||||
О |
|
|
|
|
з |
|
/ |
|
|
|
|
|
|
Imx2dx = |
|
ml3 |
|
|||
Следовательно, pf = со2, |
это приводит к формуле Лэмба — |
|||||
Саисвелла |
|
|
|
|
|
(52) |
р2 = |
|
Р? + |
|
“ 2, |
т. е. квадрат собственной частоты вращающегося стержня при
ближенно (с преуменьшением) |
равен сумме квадратов собствен |
||||
|
|
ной частоты, |
вычисленной без |
||
о сj |
|
учета вращения, и угловой ско- |
|||
|
рости вращения. |
|
|||
|
т |
Эта |
формула |
справедлива |
|
I- |
Л |
и в случае, показанном на рис. |
|||
- L -------- |
f |
28, так |
как |
|
|
|
|
|
|
||
Рис. 28 |
|
|
|
|
|
но N = та>21, что опять приводит к формуле (52). |
|
||||
Пример 6. Сравнить |
результаты, получаемые |
для |
схемы, |
приведенной |
на рис. 28, по формуле Лэмба— Саусвелла, энергетическим способом и при по мощи формулы для коэффициента жесткости (см. А® 10, табл. 1).
52
По формуле Лэмба— Саусвелла находим
3EJ
Р2 = — -г + ю2. ml8
Воспользуемся энергетическим способом, приняв для f(x) кривую прогибов от равномерно распределенной нагрузки:
Тогда получим
|
* |
|
m g |
|
f £/(/')*<** = |
i |
L |
J |
N (f')*dx = mefll J ( f ) 4 x = Y |
о |
о |
По формуле (45) находим (с завышением)
= |
Г3^о |
, |
JL |
2,2]. «2 |
3£У |
1,20со2. |
/3 |
+ |
5 |
/o j *т /о = |
mJ3 |
|
Истинное значение собственной частоты располагается в довольно узком интервале
3EJ |
-t- со-4< |
< |
3£У |
1,2ю>. |
ml3 |
тР |
Точное решение задачи может быть получено при помощи формулы для
коэффициента жесткости
a?EJ
с— аIch al — shat
вкоторую следует подставить
т Г |
N |
, Г |
ml |
o=V |
-E T - * V I T - |
||
После этого собственная частота |
определяется |
по формуле ро = — . |
|
|
|
|
т |
5. ВЛИЯНИЕ СИЛ НЕУПРУГОГО СОПРОТИВЛЕНИЯ НА СВОБОДНЫЕ
КОЛЕБАНИЯ ЛИНЕЙНОЙ СИСТЕМЫ С ОДНОЙ СТЕПЕНЬЮ СВОБОДЫ
Общие сведения
Выше считалось, что рассеяние энергии при колебаниях не происходит, и был установлен незатухающий характер процесса свободных колебаний. Опыт, однако, показывает, что колебания упругой системы, вызванные однократным возмущением, посте пенно затухают. Причина затухания состоит в том, что при ко лебаниях, кроме упругих сил, развиваются силы неупругого со противления. На преодоление неупругих сопротивлений непре-
53
рывно в необратимой форме расходуется энергия, вследствие че го постепенно убывает общий запас энергии и уменьшаются ам плитуды колебаний.
В зависимости от природы сил неупругого сопротивления для их описания 'пользуются следующими упрощенными представ лениями.
Вязкое сопротивление. Так называют сопротивление, пропор циональное скорости. Если, как обычно, обозначить скорость че рез v, то сила вязкого сопротивления описывается выражением
R — kv, |
(53) |
где k — коэффициент пропорциональности.
Такова, например, система с гидравлическим амортизатором (рис. 29), который создает сопротивление движению поршня, за
висящее не от перемещения |
(как это свойственно |
упругим |
свя- |
|
w ////////////////// |
зям), а от скорости |
и пропорцио |
||
НЕ |
нальное |
ее первой |
степени. |
По |
добные |
устройства |
применяются, |
шнапример, в конструкциях авто мобильной подвески. Гидравли
Рис. 29 |
ческий |
амортизатор |
состоит из |
||
|
одного |
или нескольких цилинд |
|||
ров с поршнями или из камеры, |
в которой |
может |
вращаться |
||
крыльчатка. Цилиндры и камера |
наполнены |
амортизационной |
жидкостью. При движении поршней или крыльчатки эта жидкость продавливается через калиброванные отверстия; этим создается сопротивление, по характеру близкое к вязкому.
В формуле (53) под R подразумевается сила, действующая на амортизатор; вязкая реакция амортизатора на колеблющееся тело имеет противоположное направление.
Сопротивление, пропорциональное п-й степени скорости. В ря де технических устройств, в частности в гидравлических аморти заторах, возможна более сильная зависимость сил неупругого сопротивления от скорости. Такие силы часто записываются в виде
R = kv\v\n~mly |
(54) |
где k — постоянная.
Кулоново трение. В этом случае силу неупругого сопротивле ния принимают постоянной по величине, но направленной проти воположно скорости:
R = ± k — kv\v\~{.
На рис. 30 даны зависимости сил R от скорости v для рассмот ренных случаев: а — вязкое сопротивление; б — сопротивление, пропорциональное квадрату скорости; в — сухое трение.
54
Вследствие неупругих сопротивлений при циклическом, нагру жении механических систем кривая сила — перемещение откло няется от закона Гука и приобретает вид петли гистерезиса
Рис. 30
(рис. 31). Если перемещению х по закону Гука соответствует сила упругого сопротивления Р, то полная сила будет больше или меньше Р в зависимости от направления движения (знака ско рости). При нагружении сила будет равна Р + R (верхняя ветвь петли), а при разгрузке она будет
равна Р — R (нижняя ветвь петли); на рис. 31 стрелками показаны на гружение и последующее разгружение.
В реальных механических си стемах причиной гистерезисных яв лений служит внутреннее трение в материале и, главным образом, тре ние в неподвижных соединениях (прессовых, болтовых, заклепочных и т. п.), возникающее при малых
проскальзываниях по контактным поверхностям. Интенсивность этого вида неупругих сопротивлений оценивается площадью пет ли гистерезиса, которая определяет энергию, рассеиваемую в конструкции за один цикл нагружения — разгружения.
На основе анализа многочисленных опытных данных Н. Н. Давиденков показал, что для ряда конструкционных материалов величина скорости не влияет на площадь Ч*1петли гистерезиса, и эту ллощадь можно представить в виде
¥ = kan+1, |
(55) |
где k и п — постоянные, зависящие от типа конструкции и мате риала;
а — амплитуда колебаний.
Последовательно рассмотрим процессы затухающих свобод ных колебаний при различных типах неупругих сопротивлений.
55
Свободные колебания при вязком сопротивлении
Рассмотрим схему, приведенную на рис. 29. При учете сил вяз кого сопротивления дифференциальное уравнение движения гру
за (рис. 32) получит вид |
*• |
|
|
• |
|
|
— сх— kx = тх. |
|
k |
С |
|
Обозначим — |
= п\ — = р2, здесь коэффициент п характери- |
|
2т |
т |
|
зует вязкость системы (его не следует смешивать с коэффициен том п на стр. 54—55). Теперь запишем дифференциальное урав нение в стандартной форме:
|
|
х + |
2пх + |
pH — 0. |
(56) |
|
Обычно |
выполняется неравенство |
|||
Рис. 32 |
р2> п2, иобщее решение уравнения |
(56) |
|||
может быть представлено в виде |
|
||||
х = ae~nt sin (|/ р2 — пН+ |
ос). |
(57) |
|||
Постоянные а и а определяются из начальных условий. Вме |
|||||
сто выражений (8а) получим |
|
|
|
|
|
/ |
2 - |
(о0 + ях0)а |
|
|
|
* > + |
|
■ |
|
||
а = arctg |
Хр У р а — na |
|
|
||
|
|
о0 + пх0 |
|
|
|
Решение (57) может быть также -переписано в виде, подобном |
|||||
выражению (8), |
|
|
|
|
|
х — e~nt ( хпcos У д 2 — пЧ 4- |
Р° + п*° sin У р 2 — пч). |
(57а) |
|||
V |
|
Ур2~Па |
) |
|
Кривая колебаний представлена на рис. 33, где отчетливо ви ден затухающий характер процесса.
Как видно из выражения (57), частота колебаний опреде ляется формулой
Р» = V Р2— ла |
(576) |
и обычно мало отличается от частоты незатухающих колебаний той же системы, но лишенной демпфирования.
Рассмотрим |
последовательные |
амплитуды, соответствующие |
|
тем моментам |
времени, |
когда sin (У р 2 — n2t + а) обращается |
|
в единицу: |
|
|
|
а±= ae~ntl\а2= |
аё~п(*1+Г); |
а3 = аё~п(<*+2Г)Г •••7 |
56
где ti — время, соответствующее первому наибольшему |
откло |
нению; Т — 2я : У р2— п2— период колебаний (следует |
иметь |
в виду условность этого термина, так как, строго говоря, движе ние не является периодическим).
Отношение двух последовательных амплитуд (декремент коле баний) остается все время постоянным:
Ч _ |
а2 _ |
-лг |
| |
'■■ |
-- |
•■•““ С |
|
CZo |
Q$ |
|
|
т. е. последовательность амплитуд образует геометрическую про грессию. Следовательно, при любом значении i справедливо ра венство
б =пТ = In |
(58) |
|
at+i |
Выражение (58) определяет темп затухания и называется ло гарифмическим декрементом колебаний.
Рис. зз
Пример 7. При колебаниях упругой системы обнаружено, что за один пе риод амплитуда уменьшается вдвое. Определить логарифмический декремент и изменение собственной частоты вследствие затухания.
По формуле (58) |
находим |
|
|
|
откуда |
|
б = пТ = In 2 = |
0,693, |
|
|
0,693 |
0,693 |
|
|
|
п |
Vp*-nK |
||
|
Т |
2я |
||
Решая это уравнение, находим, что п2 весьма мало сравнительно с р2: |
||||
|
|
я2 = |
0, 012р3. |
|
Собственная частота колебаний |
|
|
||
Р* = |
У Р2—п2 = V р2— 0, 012р3 = 0, 994р |
|||
отличается от частоты |
соответствующих незатухающих колебаний всего на |
|||
0,6% . |
|
|
|
|
57
Влияние небольшого вязкого сопротивления на частоту весь ма незначительно; вместе с тем даже малое сопротивление ин тенсивно гасит свободные колебания. Это позволяет, с одной сто роны, при вычислении частоты не считаться с наличием вязкого сопротивления, а, с другой стороны, считать свободные колебания практически исчез нувшими по истечении доста точно большого промежутка времени. В рассмотренном при мере после 10 периодов ампли
туда составляет всего (0.5) 10« «0,001 от начальной амплиту ды.
Чтобы построить фазбвый портрет системы с затуханием, нужно рассматривать совокуп ность уравнений
х = ae~nisin (p j + |
а); |
х = ae~nt [/?* cos (p*t + |
а) — |
— nsin(p*f + а)] |
|
Рис. 34
как уравнение фазовой траек тории в параметрический фор ме.
Типичная фазовая траектория изображена на рис. 34. Она представляет собой спираль, накручивающуюся на начало коор динат. Фазовый портрет образуется системой таких спиралей, окружающих особую точку О, О; последняя в этом случае назы вается устойчивым фокусом.
Свободные колебания при сопротивлении, пропорциональном л-й степени скорости
В этом случае вместо уравнения (56) получим нелинейное дифференциальное уравнение движения в виде
mx + kx п — I -j- сх — 0;
это уравнение не может быть решено в замкнутой форме. Для практических целей можно воспользоваться приближенным спо собом энергетического баланса.
Идея этого способа состоит в приравнивании энергии, теряе мой системой, работе, совершаемой силой иеупругого сопротив ления за один цикл колебаний; при этом обе приравниваемые величины выражаются приближенно (для точного их определе ния необходимо было бы знать закон движения).
58
Прежде всего выразим убывание энергии за один период ко лебаний (рис. 35, а). Потенциальная энергия в двух последова тельных состояниях наибольшего отклонения (см. точки in i + 1)
П-i — ~~ccii, llj+i = Сй^i,
где с — коэффициент жесткости системы. Уменьшение энергии за рассматриваемый период
ДП, = -j- (а? — а?+1).
При умеренном затухании последнее выражение может быть
записано проще: |
|
ДП, = у {cii + af+1) (at — a*+i) ^ caAah |
(59) |
где Adi = di — di+\ представляет собой уменьшение амплитуды за один период колебаний. Согласно рис. 35, б можно принять
Да = - Г - | - , |
(60) |
где Т — период колебаний, который допустимо считать совпадающим с периодом свободных колебаний систе мы без затухания. Оконча тельно имеем, опуская ин дексы,
Д П = — cdT-Ц -. |
(61) |
|
||
dt |
v |
1 |
*) |
|
Определим теперь |
рабо |
|||
Рис. 35 |
ту, совершаемую силой не упругого сопротивления R=
kx\x\n~l также за один цикл колебаний; при этом предположим, что в течение этого цикла колебания приближенно описываются синусоидальной зависимостью
х — dcos pt |
(62) |
(здесь принимается, что частота затухающих колебаний не отли чается от собственной частоты колебаний системы без сопротив
ления) .
Выражение указанной работы имеет вид
Y = j Rxdt = |
j kx%|x I"-1 *# = |
k j |x |n+ldt. |
0 |
0 |
0 |
Подставляя сюда согласно зависимости (62)
х = —ар sin pt,
59
получим
|
|
|
Ч1*= kan+1pn+1 j |
I sinpt\n+ldt, |
|
|
|||
или |
|
|
|
0 |
|
|
|
||
|
4 |
= kan+1pnS, |
|
|
|
||||
где интеграл |
|
|
(63) |
||||||
|
|
|
|
|
|
|
|||
|
|
|
S = |
J| sinY |rt+1 dy |
|
|
(64) |
||
|
|
|
|
о |
|
|
|
|
|
зависит только от показателя степени п: |
|
|
|
||||||
п |
................. |
О |
0,5 |
1,0 |
|
1,5 |
2 |
2,5 |
3 |
S |
.................' |
4,000 |
3,500 |
3,142 |
|
2,874 |
2,666 |
2,498 |
2,356 |
|
Приравнивая выражения Ля и ¥ (формулы |
(61) и |
(63)], по |
||||||
лучим дифференциальное уравнение для верхней огибающей |
|||||||||
|
|
|
— Ьап = |
da |
|
|
(65) |
||
где коэффициент |
|
|
dt ’ |
|
|
|
|||
|
kpn+lS |
|
|
|
|||||
|
|
|
|
Ь= |
|
|
(66) |
||
|
|
|
|
|
2пс |
|
|
||
зависит только от параметров системы. |
|
|
|
||||||
|
Рассмотрим |
случай, когда п — 1, т. е. случай вязкого сопро |
|||||||
тивления, Вместо уравнения (65) имеем |
|
|
|
||||||
|
|
|
— Ьа — |
da |
|
|
(67) |
||
|
|
|
|
|
|
И Г ’ |
|
|
|
причем выражение (66) в данном случае принимает вид
Ь = k
2т ‘
Разделяя в формуле (67) переменные
—bdt = —
иинтегрируя в пределах от t — 0 до t, находим
In— = — fc, Оо
т. е.
_ л Л—М |
_kt_ |
„ „ 2т |
|
а = а0е |
— а0е |
60
Таким образом, вновь |
получено прежнее выражение |
для |
|
огибающей кривой затухающих колебаний. |
|
||
Обратимся к случаю, |
когда п Ф 1. Разделяя переменные |
||
в уравнении (65), имеем |
|
|
|
— belt = da |
|
|
|
Интегрируя в пределах (0,/), получаем |
|
||
а1—П |
1—п |
|
|
ао |
— Ы, |
|
|
1 — п |
1 — п |
|
|
|
|
||
откуда |
|
|
(68) |
а = п—1 |
|
Этим выражением и определяется закон затухания свободных колебаний при любых значениях п, отличных от единицы. Конеч но, это решение не может дать Q
подробностей процесса затухаю щих колебаний x(t) но такие 0 подробности малосущественны.
Остановимся на двух частных случаях.
Случай п = |
0. |
В этом случае |
0 |
|
из выражения |
(68) получаем |
|
||
а = а0— Ы, |
(69) |
Рис. 36 |
||
т. е. огибающая |
имеет вид пря |
|
мой линии. Согласно формуле (66) при п = 0
(70)
этим выражением определяется темп затухания.
Случай п = 2 (гидродинамическое или турбулентное демпфи рование). Из основного выражения (68) имеем уравнение оги бающей
1 -f- ba0t
т. е. огибающая имеет характер гиперболы.
На рис. 36 схематически представлены огибающие для случа ев п = 0, 1, 2.
Свободные колебания при кулоновом трении
В случае сухого трения |
уравнение |
движения записывается |
в виде |
|
|
т х + |
сх + k = |
0, |
61