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

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

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

Тогда

и

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, здесь коэффициент п характери-

т

 

зует вязкость системы (его не следует смешивать с коэффициен­ том п на стр. 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 =

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

Разделяя в формуле (67) переменные

bdt = —

иинтегрируя в пределах от t — 0 до t, находим

In— = — fc, Оо

т. е.

_ л Л—М

_kt_

„ „

а = а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