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

Глава третья

КРУТИЛЬНЫЕ КОЛЕБАНИЯ СИСТЕМЫ С ОДНОЙ СТЕПЕНЬЮ СВОБОДЫ

1. Исходные уравнения

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

ми напряжениями т и деформацией сдвига у применитель­

но к симметричному циклу для восходящей и нисходящей ветвей петли гистерезиса, записанных по аналогии с не­ линейными зависимостями (2.32):

7 - o [ T + 4 » ( * - * , - £

) j ;

Г - о [т - 4 и (т. + 2 т - ^

- ) ] .

где G — модуль упругости при сдвиге; Ya — амплитудное

значение деформации сдвига; 26 — сумма декрементов ко­ лебаний, обусловленных различными источниками потерь энергии колебаний.

Деформация сдвига элемента пружины — материала стержня, находящегося на расстоянии р от центра сечения круглого стержня,— определяется известной формулой

7 = Р-$Г = РФ'>

(3-2)

а амплитудное значение сдвига в данной точке

р (-§-) = РФт>

(3-3)

\ил /шах

 

где <р — угол закручивания; х — координатная ось, совпа­

дающая с осью стержня.

Подставляя выражения (3.2) и (3.3) в (3.1), йолуч£%

 

т = о[Рф '± 4 б2 (рф; =Р W

- Р

(М ;

Крутящий момент, действующий в сечении стержня, Оп­

ределяется по формуле

 

 

 

 

 

 

М =

J трdF.

 

(3-5)

 

 

 

F

 

 

 

.Для стержня сплошного круглого се­

 

чения

наружным

радиусом

 

Pmax="r

 

формула (3.5) с учетом (3.4) мол^т

 

быть представлена так:

 

 

 

 

 

2яг*

Оф' ±

 

 

 

4

Рис. 8.

Схема кру­

Зя

 

 

 

± 4

Ф«Ч=

 

тильной

колебатель­

 

ной системы с одной

 

ф'а

 

 

степенью свободы.

 

>

 

 

 

 

 

 

 

Фт

 

ИЛИ

 

 

 

 

 

 

 

 

 

 

М = GIP

Ф* =F 2ф'

 

(3.6)

где

 

яг*

 

 

 

 

 

 

 

 

2 ’

 

 

 

Поскольку

 

 

 

 

 

O Irf~G Ip-%- = My

 

(3-7)

представляет собой «упругий» крутящий момент в любом

поперечном сечении стержня, то формулу (3.6)

сокращение!

можно переписать так:

 

М = Му + М„

(3.8)

-*•

 

где М3 — крутящий момент, обусловленный внутренним «тре-

нием». Согласно выражению (3.6)

4 - ± f e s0 /p (t'„ = F 2 < p '~ ^ .) .

(3.9)

Учитывая, что

4,' “ Tf-= T = t‘“lst'

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

ЛО

(ЗЛО)

<P=0^ - = qy

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

Ml

МУ1

,

Msl

(3.11)

ф - Glp -

GIp

+

GIp

 

Обозначим второе слагаемое

правой

части

выражения

(ЗЛ1) через <ps. С учетом формулы (3.9)

 

 

Ф ,« ± -§ -6 я (ф « = Р 2 ф --^ ).

(3.12)

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

ф= Фу Л- Ф«1

(3-13)

а уравнение свободных колебаний груза с моментом инер­ ции /, пользуясь принципом Даламбера, можно предста­ вить в виде

/ - ^ - + ^ = ° .

(З-М)

где с' — жесткость стержня при кручении, с'=

GI

. Учи­

тывая (3.13), уравнение (3.14) можно переписать в виде

J "ЩГ" + сг (Фу + Ф«) = 0-

(3.15)

Так как здесь присутствует составляющая угла <р5, значе-

—¥

4—

ния которого при восходящем (<р«)

и нисходящем (ф$) дви­

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

^ -^ Е-+ с'[ф у +«Г(<Р')]«0,

(3.16)

8/ (ф')= Ф«= ± -I-

fq>m=F 2ф-

(3.17)

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

Дифференциальное уравнение установившихся вынуж­ денных колебаний под действием на концевой груз внеш­ ней периодической силы eq cos pt может быть записано в виде

J -^ Г + С [фу + е/(ф)] = eqt cos pt.

(3.18)

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

2. Свободные колебания

 

 

 

 

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

уравнение

свободных

колебаний

(3.16) и вводя обозначение

 

 

 

 

®

/

II

 

(3.19)

 

 

 

j*

 

 

 

 

— (о2е/ (ф) =

еФ (ф),

 

(3.20)

уравнение (3.16) можно представить так:

 

 

dp

+ о2ф =

еФ (ф),

 

(3.21)

где

 

 

 

 

 

е| = q=

(Фто =F 2Ф -

.

(3.22)

Напомним, что член еФ(<р), входящий в уравнение (3.21), учитывает рассеяние энергии в материале и при вращении груза в разных направлениях его значение будет различ­ ным. При этом верхние знаки в уравнении (3.22) следует,

брать при направлении стрелки вправо (Ф), а нижние —

при направлении стрелки влево (Ф). В идеализированном

44

случае, если бы материал подчинялся строго закону Гука, то возмущение в уравнении (3.21) отсутствовало бы (е=0) и колебания системы были бы чисто гармоническими:

 

cp =

Mcos0

 

(3.23)

с постоянной амплитудой и равномерно

вращающимся фазо­

вым углом

 

 

 

 

4 г = о;

7 Г =

0;

0 =

(at + Ф,

т. е. амплитуда и и фаза колебания 0 будут постоянными во времени величинами, зависящими от начальных условий.

Наличие возмущения (ет^О) приводит к появлению в решении уравнения (3.21) обертонов, обусловленных зави-

симостыо мгновенной частоты d0 от амплитуды, и вызы­

вает в рассматриваемом случае постепенное уменьшение со временем амплитуды колебаний из-за рассеяния энер­

гии в колебательной системе.

 

М. Крылова,

В связи с этим, следуя методам Н.

Н. Н. Боголюбова

[1], общее решение указанного уравне­

ния (3.21) со слабым возмущением будем

искать в виде

разложения в ряд

по степени малого параметра е:

 

Ф = « cos0 +

(и, 0) + гЧг (и, 0) + ... +

&тит(и,0),

(3.24,

где щ{и, 0), и2(и,

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

функции 0

с пе­

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

= еЛ, (и) + еМ2(п) + ... + еМ т (и)\

(3.25)

= © + zBi (и) + &гВг («) + ... + &тВт (и).

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

_|_ (|)2ф =

е J — 2©A sin 0 — 2mBi cos 0 + ©2

_j_

+ e2 y[At

—uB\ — 2©«jB2j cos 0 — (2©Л2 -f 2Л1В1 +

+

u)sin9 + 2 » Л ,-Ц . +

+

+ о>2

д*и3

4- 8* ...

(3.26)

ае*

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

 

еф (ф) = еФ cos 0) 4

в2| M^

(«, cos 0) -f- Ai cos 0 +

 

4

mBt sin 0

+

Ф;, (и, cos 0)J + 83 ...

(3.27)

Приравнивая коэффициенты

при

одинаковых

степенях е

в правых частях выражений

(3.26) и (3.27) к членам т -го

порядка включительно, получаем

 

 

 

и2

4

=

F0 (и, 0) 4 2аА1sin 0 4

2<miBl cos 0;

(о2

4

«2j =

Ft (и, 0) 4

2о)Л2sin 0 4

2(ouBacos 0;

(О2 (ip r + um) = фт - 1 (и, 0) 4 2соЛт sin 0 4 2т Втcos 0

 

 

 

 

 

 

(3.28)

где

Ф0(«, 0) = Ф(исоз0);

 

 

 

 

 

 

 

 

Ф1(«, 0) =

щФ’и(иcos 0)

cos 0 — иВ±sin 0

4

 

+

ш

ф;. (и COS 6) +

A,

cos в +

^3'29'

+

(М А

+ Л, £ ») sin в -

2А, g f r

- Ы » ,

 

 

Интересуясь в дальнейшем решением в первом приближе­ нии, определим щ ( и , 0), A i { u ) , B i ( u ) . С этой целью пред­ ставим разложения Ф0(«, 0) и 0), входящие в первое уравнение (3.28), в виде рядов Фурье

ф0 (и. 0) =

ёо («) 4

2

[£п («) cos П0 4- h (и) sinпЩ;

 

 

Л=1

 

 

 

(3.30)

 

 

00

*

щ (и, 0) =

и0(и) 4 -

2

[*>» («) cos п0 4- wn (a) sin пО],

 

 

л=1

 

 

2 я

 

£Го (w) =

у <Р0 («»в) rf0;

 

 

О

 

 

 

8п(“) ~

| Ф0 (и, 0) cos п0<20;

(3.31)

 

о

 

1 2я

hn (и) — ( Ф„ (и, 0) sin nQdQ.

Подставляя правые части уравнений (3.30) в первое урав­ нение системы (3.28), получаем

 

00

 

 

 

ш2п0 (и) +

2

ш2(1 — «2) [vn (и)cos пв -f- wn (и) sin Л0] =

 

л = 2

 

 

= go (и) +

[£Ti (ti) + 2шВ1) cos 0

+ [hi (и) + 2соЛ1] sin 0 +

 

+

2 t£n (“) cos л0 +

К (и) sin л0].

(3.32)

 

 

л = 2

 

 

Приравнивая в последнем уравнении коэффициенты при одинаковых гармониках, будем иметь

gi (и) +

2cauBi = 0;

ht (и) + 2<оЛА=

 

0;

о л л

_

So (и) .

(

) _____Sn (“ )

.

(3.33)

V0 \U)

0)2

Vn W

0 ) 2 ( 1 — Л2)

wn(u)=

h

(и)

;

(« = 2»3,4,...),

 

0)2 (Г-П»)

 

откуда с учетом выражений (3.31) получим

Ai

М “) .

Bt = -

gi(“)

(3.34)

 

2© *

 

2

 

Так же определены все гармонические компоненты функ­ ции «i(«, 0), кроме первых t>i(и) и щ (и). Однако в силу того, что эти функции не содержат первой главной гар­ моники, будем иметь i>i(«)=0; Wi(u) = 0 . Тогда на осно­ вании выражений (3.30) найдем

„ /„ m _ go («) .

1 V U n W cos n0 + К (“>sin

(3.35)

- 5 Г +

-з г 2 а

 

 

л=2

 

С помощью изложенного метода можно определить «п(и, 0). Ап(и), Вп{и) (/1=1, 2, 3,...) до какого угодно значения п и тем самым построить приближенное реше­

ние,

удовлетворяющее

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

(3.21)

с точностью до

величины сколь угодно высокого

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

В силу обоснования, изложенного во второй главе, при­ менительно к рассматриваемой задаче в качестве первого приближения для функции углового перемещения следует

принять упрощенное выражение

 

(p = ucos0,

(3.36)

в котором функции времени и и 0 должны быть

опреде­

лены из дифференциальных уравнений

 

т г = в + вв,(и),

(3.37)

где согласно выражению (3.34) с учетом (3.31)

 

2п —►

 

4 М- - - - таг f «Ъ(«• в)staел;

(3.38)

= -------

3 5 5 Г $ е ф (“ - в>с м ш ■

 

О

Выражение еФо(«, 0) в явном виде на основании (3.22) и

(3.29) с учетом (3.39) можно записать так:

еФ (и, 0) =

q= -i-

(1 =F 2 cos 0 — cos20), (3.39)

а интегралы, входящие в выражения (3.48), будут

2я 2

 

Г2П

j еФ0 (и, 0) sin 0d0 =

-g-

zu I Г (1 — 2 cos 0 — co$20J0) X

о

 

Lo

= оуЧ^и;

(3.40?

\ еФ (и, 0) cos Ш =

3

г мя

 

-g- CO2S2U

J (1 2 cos 0 — cos2 0rf0) x

X cos Odd — ^ (1 +

2 cos 0 — cos2 0) cos BdQ

x

 

 

 

 

X [2я] =

Зл

(3.41)

 

~ o ) 26yu.

Подставляя выражения (3.40) и (3.41) в (3.38), будем иметь

Л ( « ) “ -

ш8su

(3.42)

2пе ’

Bt (tt)

(3.43)

 

 

Подставляя выражения (3.42) и (3.43) соответственно в уравнения (3.37), будем иметь

 

“ЗГ -еЛ ^ м )

(3.44)

 

 

 

40

со еВ (и) = to —

(3-45)

4/

 

 

 

В случае, когда декремент колебаний 6 j= /(« ),

урав­

нение (3.42) легко интегрируется и мы будем иметь

 

или

1пв= . т г + с

Р « )

©бя*

■с

 

 

 

 

f u r

(3.47)

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

(а)<=о = а0 = е°-

(3-48)

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

©б^

 

а = а0е 2п

(3-49)

В том случае, если декремент бг зависит от амплитуды де­ формации, как это имеет место при учете рассеяния энер­ гии в материале, интегрирование уравнения (3.45) будет

сложнее. Интегрируя уравнение (3.45), найдем выражение для фазы

0 =

0o + d>f— -f-ооб/

(3.50)

где 0о — начальная

фаза, определяемая из

условия

0о= (0)<=о.

 

 

3.Вынужденные колебания

врезонансной золе

Перепишем дифференциальное уравнение (3.18) вынуж­ денных крутильных колебаний (3.18) с учетом (3.19) в виде

~dt*

(ог<р = еq cos p t — еФ (ф),

(3.51)

 

 

где

 

 

 

еФ(ф) = co2ef (ф),

(3.52)

со — собственная круговая частота; q = - j- .

Используя принятый выше общий метод построения асимптотических решений, будем искать общее решение (3.51) в виде следующего разложения по степеням малого параметра:

Ф = « COST -f e«t (« ,т) + е2«2(«, т) +

(3.53)

где

г = p t -f- ар.

Амплитуду деформаций «, фазу т, а следовательно, и сдвиг фаз определим из дифференциальных уравнений

~ljf- = eAj (и) + еМ2(и) + ...;

(3.54)

■^- = (о + еВ1 (и) + егЯ2 (« )+ •••

или уравнения

= © _ р + еВt (и) + е2В2 («) +

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

~ р

«

еЛ, (и, ф) + еМа (и, ф) +

 

=

©— р +

eBj (м, ф) +

е25 2 (и, ф) +

(3.55)

Однако, поскольку

нас

интересует

резонансный

случай,

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

Взяв вторую производную от выражения угла закручи­ вания qp согласно выражению (3.53) с учетом (3.54) и под­ ставив в левую часть уравнения (3.51) с учетом (3.54) и затем в полученном выражении собрав члены при одина­ ковых степенях е и приравняв их правой части уравнения (3.6), содержащего малый параметр в первой степени,

получим

 

 

 

 

 

2(лА1sin т — 2m Bl cos т +

со2 ( +

и } = q cos pt Ф(ф).

 

 

 

 

1

(3.56)

Далее, умножая

уравнение

(3.66) на

smxdx

и cos xdx и

интегрируя от 0 до 2я, будем иметь

 

 

<о2 J

+ ui\ s*nт^т=

2nco/4i -j- q ^cos pt sin xdx

 

 

 

 

о

 

 

6

 

 

 

 

 

 

 

 

— jj Ф (ф) sin xdx]

(3.57)

2rt

 

о

 

2rt

 

 

 

 

 

to2 j (4 ^ —

Wjj cos xdx =

2ижоВ1 +

q ^ cos pt cos xdx —

6

 

 

 

о

 

 

 

— ^ Ф (ф) cos xdx.

 

(3.58)

 

 

6

 

 

 

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

точно в резонансе

= 0,

т. е. ф = const, обозначая 0 = pt,

имеем

 

 

dx = d0;

т =

0 +

ф;

2П

in

Гcos 0cos (0 if) dQ =

cos

j cos20d0 — sin я|> f cos 0 sin 0d0 =

J

n

 

n

 

=

ncosi|);

 

 

 

2lt

Г cos 0 sin (0 +

г|>) dO = cos iji Гcos 0 sin 0d0 -f

0

n

 

°

 

2rt

 

 

+ sin i|>j

cos2 0 = JX sin %

а также

 

 

 

j

+«iJsinT dT = 0;

0

 

 

 

 

 

 

]

+

Wl) cosxd% =

так как

 

 

2Jl

 

 

 

^-^1- sinxdx =

— j u t sinxdx = 0.

о

 

 

о

Выражения для определения искомых At и Вх найдем со­ ответственно из уравнений (3.57) и (3.58):

2Я-»

^ Ф ((р) sinxdx nq sin

 

Л = :

2гов

(3.59)

2Я_

 

 

J Ф (ф) cosxdx — nq cos i|>

 

Bi =

2unco

(3.60)

 

 

Подставляя полученные выражения Ai и Bi в уравнения

(3.54), будем иметь

 

 

du

t 2Я

^

 

 

о

®®(ф) s in t d r - -g-sinil);

(3.61)

dx

2я ^

 

1

 

Т Г *" Р +

2лыйГ 1 еФ (Ф)C0STrfT 2и&Г C0SФ

C3 '62)

или

 

0

 

 

 

 

dib

i

2n 2

 

7 r = ffl_ p + 5 _ J e ® ( ? ) c o S T * - - ^ r cosi|j. (3 .6 3 )

Ли

 

dt

=* 0;

I f

= 0 или

= /?.

 

Выражение для

синуса сдвига фаз найдем из уравнения (3.61):

 

 

 

t

2я ^

 

 

 

 

sin ^ =

 

f £ф М sin TdT*

(3-64)

 

 

 

 

6

 

 

а выражение для частоты — из уравнения (3.63):

 

 

 

 

 

2Я ^

 

 

 

 

 

cos ф — ^еФ (ср) costdr

 

 

 

 

_________ о____________

(3.65)

 

 

 

 

исо2

 

 

 

 

 

 

 

Исходя из выражения (3.62), можно

найти косинус сдви­

га фаз:

 

 

 

 

 

 

cos ф — ±

 

 

 

(3.66)

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

 

 

^еФ (ср) cosxdx ± V

(е^)2я2 —

 

 

р

_

 

 

2яисо2

 

(3.67)

ю

= 1+ -

 

 

 

 

Учитывая согласно выражениям (3.40)

и (3.41),

что

 

 

 

 

 

 

j еФ (<р) sin x d x = са^ и ;

о

£ 6 Ф (<р) COSxdx = З ж и ^ п ,

о

формулы (3.66) и (3.67) можно окончательно представить в виде

 

cos ф = ±

л Г (е?)ал2 — ш28,и

(3.68)

 

L_ --------

а -

 

Зя

щп

 

 

 

У(е?)2 я2

— (o28s«

р_

 

1 +

2л«ш2

 

(3.69)

ш

 

Формулы (3.68) и (3.69) и являются теми окончатель­ ными формулами, пользуясь которыми можно построить резонансную амплитудно-частотную кривую вынужденных крутильных колебаний системы с одной степенью свободы.