Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
книги / Обобщенная нелинейная модель учета рассеяния энергии при колебаниях..pdf
Скачиваний:
0
Добавлен:
19.11.2023
Размер:
29.24 Mб
Скачать

Глава пятая

СВОБОДНЫЕ КОЛЕБАНИЯ СИСТЕМ

с р а с п р е д е л е н н ы м и п а р а м е т р а м и

1. Продольные колебания призматических стержней

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

Дифференциальное уравнение продольных колебаний стержня получим из рассмотрения условий динамического равновесия элемента стержня длиной dx (рис. 11) с ис­ пользованием принципа Даламбера:

 

d N = > p F d x - $ - ,

 

(5.1)

где и — перемещение поперечного сечения вдоль оси;

р —

плотность

материала стержня; F — площадь поперечного

сечения стержня, dN = N z—Nl — разница

внутренних

уси­

лий в сечениях 2 и 1, находящихся на расстоянии dx:

^! =

Г((Т + ав); Nz = F(o + oa)+ F

+ aa)dx

 

или

—^

 

 

 

 

 

dN — F -Цг (а + <тв) dx.

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

! = - § - ; | = | acos0; e - e r f + ф,

где (cos0),=o = соэфл* 1, можно записать

; О. = ± 4

^ 2 Cos0 - C°s20)‘

Обозначая

д

 

 

 

 

дх [ ± т б2 - ^ - (1 т 2 C°S0 — cos*0)| = еф (м),

(5.2)

выражение для dN можно записать в виде

 

dN = E

F ^ +

Ф (U)jdx,

(5.3

а подставляя выражение (5.3) в (5.1), получаем

 

+

еФ(ы) =

I -

 

дх*

^

К )

Е dt*

 

ИЛИ

 

 

 

 

д2и

 

rfin

ч—

(5.4)

дх*

P2jjZ - = - e ® (u ),

где еФ(и) — член, учитывающий потери энергии в колеба­ тельной системе, вносящий определенную нелинейность, которую полагаем малой, о чем свидетельствует стоящий

при нем

множителем

малый

и I

 

ди

,

параметр е.

уравнения

(5.4)

 

Т

и + Ц

Г

^

Решение

 

 

 

 

 

для определения

перемещения

£ ■

 

 

 

 

любого

сечения

стержня

на

 

 

dx

 

 

расстоянии

х

от

 

выбранного

 

 

 

 

 

 

 

 

 

 

начала координат 0 будем ис­

 

 

 

 

 

кать в виде ряда

 

 

 

 

N

 

" 'З г *

и (х, t) =

Яф (х) cos 0 + е% (х, 0 )+

 

 

 

 

 

X

+

е2иа(*,0) +

...,

(5.5)

 

 

 

 

1

к

ди

,

 

 

 

 

 

 

 

 

 

 

где амплитуда колебания Я и

 

 

дх

 

 

фаза 0 могут быть определены

Рис. 11. Схема

к

продольным

из следующих дифференциаль­

колебаниям.

 

 

 

 

ных зависимостей:

 

 

 

 

 

 

 

 

 

 

dk =

eAL(Я) +

еМ2 (Я) +

 

 

(5.6)

 

 

- § -

=

 

м + е В 1(А.) +

е‘В2 (Ь) + ...

 

 

(5.7)

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

jp- и -^5-в уравнение (5.5), а также учитывая уравнен*1*

(5.6) и (5.7), собирая коэффициенты при одинаковых с'г^ пенях в полученном при этом выражении и приравнив^з их нулю, получаем следующую систему уравнений:

cos6 + ty2(o2<p(x)cos0 = 0;

(6-0

-■jfip- 4-

(х) sin 0 +

2%ргВ1щ (х) со s 0 —

 

- PV ^ . = T

4 - S^ - S

- ( 1=I= 2 co s9 - cos29)-

(**£

Выражение функции ф(х) найдем из уравнения (5.8), Кс

торое может быть представлено в виде

 

^ & М _ + *М *) = 0,

(5-1*

где

№ = р2©2.

 

 

 

Решением этого уравнения будет

 

Ф (х )

= Сх cos kx 4- С2 sin kx.

(5.1)

Постоянную интегрирования Сх

и значения параметра k опр-

делим

из граничных

условий

<р (0) =

0; <р' (/) =

0,

откуа

Сг = 0,

т. е. ср (х) =

C2sin kx.

Тогда

cosxl =

0;

kl =

-у- <

X (2п — 1) (п — 1,2,3 ...). Для

п — 1,

kl = Y

или ра>1 —

следовательно, © = , ср (х) = Сг sin -jj~x.

Так как <р(х) выражает форму колебаний, то моя^о положить Сг=1 с точностью до постоянного коэффицие^Д и записать выражение для <р(х) в виде

<р(х) = sin -^ -x .

(5.2)

Используя условие, что щ, и%,... не содержат главных гРмоник, можно определить Ai и Bt. Для этого умноя{М (5.9) сначала на ф(х) cos QdxdQ и проинтегрируем от (Но I и от 0 до 2я, затем — на ф (х) sin QdxdQ и также npdHтегрируем в тех же пределах. При этом будем иметь

12л

о о

~ * г^ ) ф(*)со5в‘ы в =

 

12П

=

^

 

|2©/)2/4i(p (х) sin 0 -f 2Ха)р2В1ф (х) cos0 ±

 

о о

 

^

 

 

 

4 Е Г

[ ~ § £ "

^

(1 =F 2cos0 — cos20)1j <р(л:) cosMxdB\

 

 

 

I 2Л

 

 

(5.13)

 

 

 

д2и

дЧ

 

 

 

 

 

 

 

-

ш

дх2

я аэа ) Ф (^) sin0djcrf0 =

 

*

 

 

 

 

о

о

 

 

 

 

12Л

 

 

 

 

=

j j

^(«рМхФ (л:) sin 0 4* 2Ыр2В1(р (х) cos 0 ±

 

о о

1

 

 

 

± ~3г[~1” бл*'

 

=F2cos0 — cos20)]W(*)sin0d*d0.

 

 

 

 

 

 

(5.14)

Имея в виду,

что в выражении (5.13)

 

 

 

$ -^г* Ф (х) dx = J

щфг,

 

 

 

Г

- cos 0rf0 = — j

ut cos OdO,

 

 

 

о

 

0

 

учитывая граничные условия w , _ „ = 0 ;

l * W U = 0 ; (-§ ■ )„ , = ° и интегрируя по частям, будем иметь

I

I

* У

m (х) =

% - * w

 

аха

ф w

дх

 

О

 

О о

 

 

 

 

I

 

 

 

| ' +

£ lW — udx

 

 

. ( ■ £ « « * * - J-

dx2 Uiax‘

Следовательно, левая часть уравнения (5.13) обратится в нуль:

12л

Оо

tin м , .

л

- f J' [_ *i£ L + /ftp(*)]«,cose<taffl = 0.

0 6

Тогда из уравнения (5.13) найдем

12П

 

 

 

t2n

ХВхф2(л:) coszQdxdQ =

f f Аг<p2 (<p)sin0 cos ddxdQ+

j J

oo

(2lt

 

 

oo

 

 

 

 

 

(1T 2 cos e - cos’6) X

- *

J } TF

( V - ^

)

* w

 

о 0

X cos QdxdQ.

(5.15)

 

 

 

2jt

 

 

Учитывая,

что f cos2 0d0 = я;

j* sin 6dQ = 0, а также выра-

 

o

 

о

 

 

жения (2.52) и (2.53), из уравнения (5.17) найдем

 

 

I

 

 

 

 

3J Чх (62Л"

 

 

ЯЛ* . ) - — *

 

 

(5.16)

 

 

8&р2а% J <ра (х) dx

 

 

 

о

 

 

Аналогично из уравнения (5.14) получим

 

 

 

 

 

dx

а

д -

 

 

 

(5.17)

 

 

2пее>р* J Ф (х) dx

Подставляя выражения (5.17) и (5.16) в (5.6) и (5.7) со* ответственно, получаем

(5.18)

dt

 

2mopaj

Фа (•*) dx

 

 

 

dx

dQ

_

(5.19)

Л_ss

Q X

&р*<й%J фа (x) dX

или, учитывая, что для рассматриваемого конкретного случая согласно формуле (5.12)

<P(*) = sin-JL*; JvgL<p(x) = -±-cos-TXi

 

l

l

 

 

f q * ( x ) d x = f s in * ~ xdx =

- j - ,

 

0

0

 

уравнения

(5.18) и (5.19)

можно переписать в виде

 

i

 

 

 

О

 

(5-20)

 

 

 

"ЗГ =

ю “ 4 /V b I Ч

Г [бая cos -J T

* \ sin ч г dx• (5-21)

 

о

 

J

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

литуды относительной деформации В том случае, когда

в « * /(Ь ),

(5.22)

из уравнения (5.20) будем иметь

 

 

d \

-

м

(5.23)

% ~

°2

4

 

Интегрируя уравнение (5.22), получаем

 

1пЛ =

Мр

+ С

(5.24)

 

или

 

 

 

 

 

 

(5.25)

Постоянная интегрирования С найдется из начальных

условий:

 

ес.

 

(Ч=о =

Ьо =

(5-26)

Тогда закон изменения амплитуды свободных продольных колебаний (5.23) можно представить формулой

л - V 4"

(5-27)

Интегрируя дифференциальное уравнение (5.19), полагая, что 6s=£f(h), будем иметь

(5.28)

о - * ' - - щ г ‘ + с »

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

(9)f=0 =

= О,-

С учетом (5.26) закон изменения фазы 0 во времени мо­ жет быть окончательно выражен формулой

(5.29)

В случае, когда бs=f(X), задача значительно усложня­ ется из-за возникающей существенной нелинейности и ее решение требует новых подходов, которые нашли отраже­ ние в работе [1].

2. Крутильные колебания стерж ня постоянного сечения

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

Учитывая, что относительная деформация сдвига эле­ мента скручиваемого стержня, находящегося на расстоя­ нии г от центра сечения, и ее максимальное значение че­ рез относительный угол закручивания <р' определяются формулами

выражение крутящего момента в сечении стержня радиу­ сом R, исходя из зависимости (3.1) с учетом формул (5.30) можно представить формулой (3.6), которую запишем в виде

= G j* r^'dF ± у 6г J г2 ^ =F 2ф' - X ~J dFj , (5.31)

где

 

 

 

Му =

G j r2(p'dF =

GIP<p';

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

M s =

± 4

G \62 (фт =F2<P' “

'ф^* )r W = GIP6® (Ф)*

(532

 

 

F

 

 

 

 

 

 

 

 

 

Для вывода дифференциального уравнения крутильных

колебаний

стержня воспользуемся

принципом

Даламбера

н составим уравнение равнове-

^

 

 

 

сия

стержня

длиной

dx

|

 

 

 

(рис.

12)

 

 

 

 

 

 

 

 

 

 

 

 

 

—►

—►

 

 

 

 

X

дх

 

 

 

— Af-f- М -f*

 

 

 

 

 

—>

 

 

 

 

 

 

 

 

 

 

 

+ Jsr< l* - U

r4F- w - d x ~ 0'

 

ft м

у

У

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

)

*

откуда

 

 

 

 

 

 

 

dx

Ал

 

 

 

 

 

 

 

 

 

ЭМ

_

г

д \

 

 

 

Рис. 12. Схема к крутильны»

 

 

 

 

колебаниям.

 

 

 

дх

~~ р1р

дР

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в последнее уравнение выражение момец

та согласно формуле

(5.31),

получаем

 

 

 

4 г [01> i t

+

0 /*еф « ]

= Р'р Ж

 

<5-33

Обозначая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ь2 —

Р

>

 

 

(5.34)

 

 

 

 

 

«

 

Q

 

 

дифференциальное уравнение (5.33) можно переписать в виде

(5.35)

где

е ф ( Ф) = ± т б2 f(pm =F 2Ф' -

Напомним, что член еФ(ф), стоящий в правой части уравнения (5.35), учитывает рассеяние энергии в колеба­ тельной системе за счет различных факторов.

Следуя принятой теории нелинейной механики, решение уравнения (5.35), обладающего слабой нелинейностью, бу­ дем искать в виде разложения функции угла закручивания ф(х, /) в ряд по степеням малого параметра е:

(х, t) = шр0 (JC) cos (ю/ + ф) + еф! (и, х, 0) + в2ф2 (и, х, 0) 4- ..., (5.36)

где ф о (х )— решение невозмущенного дифференциального уравнения (3.35) в случае е = 0 , т. е. уравнения

 

дх*

— bzJ^Ss.

= 0;

(5.37)

 

к дР

 

 

ф1(и, х, 0),

фг(и, х,

0),... — периодические

функции

0=шЦ-ф с периодом 2я; и и 0 — функции времени, опре­ деляемые из дифференциальных уравнений

-ЗГ = е4(«)-ЬеМ 2(м)-{-

(5.38)

- f -

= со + еВ, (и) + в*В2 («) + ...;

(5.39)

со — собственная

частота колебаний стержня

при отсут­

ствии рассеяния энергии, движение которого описывается уравнением (5.37).

Подбираем выражения для функций

 

Ф!(ц,х,0), фа(«, х, 0 ),...;

 

М 2(«),...; Вх(и),Вг (и) ,...

(5.30)

так, чтобы уравнение (5.36) с учетом (5.38) и (5.39) ока­ залось бы решением уравнения (5.35). Для обеспечения однозначности будем полагать, что функции фДм, х, 0), Фг(и, х, 0),... не содержат первых гармоник, т. е. должны выполняться условия

 

 

Ф! (и, х, 0) cos QdS =

0;

j

фх (и, х, 0) sin QdQ= 0;

 

ЗЯ

(5.41)

 

ф2 (и, х, 0) sin BdQ 0;

ф2 (и, х, 0) cos QdO=

0;

J

Далее, взяв вторые производные от уравнения (5.36) по х и по t и подставив их в левую часть уравнения (5.35),

получим

 

4 5 — *

= “ r S t +

0 +

 

+

№& 2(0Л («) Фо sin 0 — 2ШБ! (и) ф0 cos 0 + и2

+

+

-jr - 0 - ] +

*2е2 { [Л

(«)

" " иВ*^ “ 2(вдВа («)] X

X Фо cos 0 — [2о>4 («) +

2АХ(и) Вх (и) +

Ах- Bj a(a)-jф0 sin 0 +

+ 2о)Лх (м)

+

2 ( о В х

+

w2 op

+ Аа -^JrJ +

в3 ...

(5.42)

Разложим также в ряд Тейлора правую часть уравнения (5.35);

еФ (ф) = е

Ф' [и, q>' (х) cos 0] +

 

+ e2®" [мфц (х) cos 0] Фх + е3 ...

(5.43)

Приравняем коэффициенты при одинаковых степенях 8 в правых частях выражений (5.42) и (5.43) в соответст­ вии с уравнением (5.35) до членов m-го порядка включи­ тельно.

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

- ^ + Л

» = 0;

(5.44)

[— 2ооф0Л1 sin 0 — 2а>и<р0Вхcos 0] +

 

+

 

<5 -« )

где

k 4 2

(5.46)

р2=

Уравнения (5.44) — (5.45) и являются теми исходными уравнениями, с помощью которых рассматриваемая задача может быть решена в нулевом и в первом приближениях.

Решение уравнения (5.44) может быть записано в виде

Фо (х) = A cos рх + В sin рх,

(5.47)

где А и В — постоянные интегрирования, которые должны быть определены из граничных условий. Для стержня со

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

Фо(0) = 0; Фо (О = 0.

(5.4^

Тогда на основании уравнения (5.47) получим уравнен^

частоты sin pi = 0,

откуда л1 tit.

Учитывая,

что р2 =

^2со

получаем sin>feco/ =

0, следовательно,

k<al — in

(i =

1 ,2 ,3 ,...

Собственная частота основного тона колебаний о> =

-гг-

ил:

Таким образом, решение уравнения (5.47) окончательш может быть записано в виде

Ф0 (х) = В sin рх = В sin -у- х,

(5.49)

функция угла закручивания в нулевом приближении со­ гласно выражению (5.36) примет вид

Ф (*. 9) = Нф0 (х) cos 0 = Ви sin х cos 9;

обозначая Ви С, получаем

Ф (х, 0) = Сsin ~ х cos 0.

(5.50)

Постоянная С определится из условия, по которому посре­ дине стержня (x=lj2) в начальный момент при 0 = 0 угол закручивания будет максимальным:

Ф° ( у ) = Сsin

= С — и.

Следовательно, С — Ви = и; В — 1.

Таким образом, согласно

выражению (5.50) угол за­

кручивания в произвольном

поперечном сечении стержня

в любой момент времени может быть определен выраже­ нием

ф0 = зт -у -х ; ф (х, 0) = иsin -у- х cos 0.

(5.51)

Для решения задачи в первом приближении

следует

ввести в рассмотрение уравнение (5.45), из которого можно найти A i ( u ) и B i ( u ) . Умножим уравнение (5.25) на Фо cos QdxdQ, а затем на <р0 sin QdxdQ и полученные выраже­

ния проинтегрируем по длине стержня (от 0 до /) и по

циклу (от 0 до 2л):

t 2л

11

+

рг 4S^“) ф°м cos Qdxde=

 

о

о

 

 

 

I 2Л

 

 

 

 

{Ф' [и, Фо (х)>cos 0] + £2 (2соЛх (и) Фо (л:) sin 0 +

 

о о

 

 

 

 

+

2©ыф0 (х)

(и) cos ®]} Фо (х) cos 04x40;

(5.52)

 

 

j

f

фо W sin Qdxdd =

 

 

 

О 6

 

 

 

=

12л

 

 

 

j* J

{ф' [а, ф (x), cos0] + А2[2(оЛх (ы) ф0 (x) sin 0 +

 

° °

 

+ 2m<p0 (x)

(u) cos eO Фо (x) s>n Qdxdd.

(5.53)

Интегрированием по частям можно показать, что

 

 

 

 

/

*

 

 

 

 

 

1

о

 

 

 

 

 

ая

 

 

 

 

 

1- ^ 8 Ш 0 4 0 = - | фхsin 040,

 

 

 

 

 

 

л

 

 

в связи

с

чем на основании выражения (5.53)

получим

 

 

 

/

 

 

 

 

 

(]■ ( % р - + Ц > ^ - ) ч ,М И п ( к Ш

=

 

 

 

12л

 

 

 

 

=

j

J

- Я2Фо (*)] Ф1sin 04x40 =

0.

 

 

о

о

 

 

 

Можно также показать, что равна нулю и левая часть уравнения (5.52). Тогда из уравнений (5.52) и (5.53) соот­ ветственно получим

I 2Я-+

^ j

Ф' [н, <p' (х) cos0] фо (х)

cos0d*d0

О, («) = - —

-----------------,------------------------ ; (5.54)

 

2 u k 2a m ^ <pjf (X) d x

 

о

/2rt-.

^ J Ф' [и, ф' (ж) cos 0] фо (ж) sin Qdxdb

Ax(u)-------

LS-----------------

т-------------------------

, (5.55)

2й2шя £ фц (ж) rfjc

о

где согласно выражениям (5.24), (5.43) и (5.51)

Ф' [и, ф' (х), cos 0J = ± 4

б2ы

(! =F 2 cos 0 ~ cos2 6)1

 

 

 

 

 

(5.56)

Подставляя уравнение (5.56) в (5.54) и (5.55), получаем

12Я

З

d /

 

 

\

И±

зт

зт

- g j - 3 7

у

cos у

ж) (1 =F 2 cos 0 — cos2 0) X

о о

 

 

 

я

 

 

 

X

sin

 

 

 

- у - 0 cos 0dxd0

Л (#)-

 

 

 

 

 

 

 

2пй2(оя

xdx

/ 2Я

 

 

 

 

(5.57)

 

 

 

 

 

1 J ^

Ж

“dx” ( fis “ 7

cos у

ж| (1 q : 2 cos 0 — cos2 0) X

Оо

 

 

 

 

 

Ai(u) =

 

X

sin -j- х cos 0dxd0

 

 

 

 

 

2kaan щ)f - f xdx

о

(5.58) Заметим, что при интегрировании от я до 2я в послед­ них уравнениях следует брать верхние знаки, а при интег­ рировании от 0 до я — нижние; учитывая это, выражения (5.57) и (5.58) после интегрирования числителей по 6 и знаменателей по х и элементарных преобразований можно

представить в виде

<5-59>

О

A^ = - i m r \ 4 r [ s^ T

* * • (5-60)

Подставляя выражения для Ai(u) и Bt(u) согласно вы­ ражениям (5.60) и (5.59) соответственно в дифференци­ альные уравнения (5.38) и (5.39), ограничиваясь при этом первым приближением, получаем

i

du

 

l

C

d

&2и cos

-j- xj sin -j-xdx;

(5.61)

dt

/2A2©

J

dx

 

 

 

о

 

 

 

 

-ЗГ

=

“ + 1 ^

r

i 4

- { ^

T x) s'm T xdx- <5-62>

 

 

 

 

о

 

 

 

Имея значения декрементов, условно показанных сум­ мой бл как функции тех или иных факторов, полученных из эксперимента путем интегрирования дифференциальных уравнений (5.61) и (5.62), найдем выражения для ампли­ туды и и фазы колебаний 0 в функции времени:

и — fi (?) + Ср 0 = (at -f- /2 (0 +

Постоянные интегрирования Ci и Сг должны быть опре­ делены из начальных условий.

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