Глава третья
КРУТИЛЬНЫЕ КОЛЕБАНИЯ СИСТЕМЫ С ОДНОЙ СТЕПЕНЬЮ СВОБОДЫ
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) |
|
где |
|
яг* |
|
|
|
|
1р |
|
|
|
|
|
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; |
|
|
О |
|
|
2я |
|
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) |
8е |
(3.43) |
|
|
Подставляя выражения (3.42) и (3.43) соответственно в уравнения (3.37), будем иметь
|
“ЗГ -еЛ ^ м ) |
2л |
(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 — |
|||
2Л |
|
|
|
2К |
|
о |
|
2л |
|
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П |
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 |
|
|
|
2Я |
|
|
|
] |
+ |
Wl) cosxd% = 0» |
|
так как |
|
|
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лисо2 |
|
|
|
|
|
|
|
|
|
Исходя из выражения (3.62), можно |
найти косинус сдви |
|||||
га фаз: |
|
|
|
|
|
|
|
cos ф — ± |
|
|
|
(3.66) |
|
Подставляя уравнения (3.66) в (3.65), получаем |
|
|||||
|
^еФ (ср) cosxdx ± V |
(е^)2я2 — |
|
|
||
р |
_ |
|
|
2яисо2 |
|
(3.67) |
ю |
= 1+ - |
|
|
|
|
|
Учитывая согласно выражениям (3.40) |
и (3.41), |
что |
||||
|
|
2я |
|
|
|
|
j еФ (<р) sin x d x = са^ и ;
о
2«
£ 6 Ф (<р) COSxdx = З ж и ^ п ,
о
формулы (3.66) и (3.67) можно окончательно представить в виде
|
cos ф = ± |
л Г (е?)ал2 — ш28,и |
(3.68) |
|
|
L_ -------- |
а - |
||
|
Зя |
щп |
|
|
|
У(е?)2 я2 |
— (o28s« |
||
р_ |
|
|||
1 + |
2л«ш2 |
|
(3.69) |
|
ш |
|
Формулы (3.68) и (3.69) и являются теми окончатель ными формулами, пользуясь которыми можно построить резонансную амплитудно-частотную кривую вынужденных крутильных колебаний системы с одной степенью свободы.