книги из ГПНТБ / Давыдов, В. В. Динамические расчеты прочности судовых конструкций учебник
.pdfт. е. если отклонение ф (t) выражается функцией косинуса, то откло
нение ф ( / + — ) — функцией синуса. Но из (а)
ф ( ^ ) = —сш sin (со^+ Р), |
(в) |
следовательно [сопоставляя (б) и (в)], получим, что
Поэтому силу гистерезисного сопротивления можно считать не только пропорциональной отклонению, которое будет иметь место через четверть периода, но и скорости в данный момент
xNФ ^ |
ф(0- |
(г) |
В случае, если постоянная сила трения Q, как это обычно бывает, невелика, ее приближенно также можно считать зависящей от скоро сти
Qsi ky(t),
что позволяет составить уравнение |
для всего процесса колебаний, |
а не только для одного полупериода, |
как приходится делать при рас- |
смотрении свободных колебаний при наличии сухого трения (§ 2, с. 27). Коэффициент k следует подобрать так, чтобы работа действитель
ной силы трения за полный цикл одного колебания
Wx= — 4Qa
равнялась работе условной силы ky, за такое же время.
При равенстве этих работ будет одинаковым и поглощение энергии, обусловленное этим видом сопротивления при обоих выражениях для
силы трения. |
|
Условная сила трения |
на элементарном перемещении |
|
|
|
dy = фdt |
|
|
совершит отрицательную работу |
|
|
|||
|
|
dW2= —kq4t, |
|
|
|
и если ф выражается формулой (а), то |
|
|
|||
|
W2= |
Т |
Т |
|
|
|
—J &ф2dt = |
— ka?со2 J sin2 (со/ -f Р) dt. |
|
|
|
|
|
о |
о |
|
|
Последний |
интеграл, как |
известно, равен половине |
своего |
верх- |
|
него предела |
Т |
и, следовательно, учитывая также, |
(ОТ |
—п, |
|
— , |
что — |
получим
W2= — ka2n<d.
30
Приравнивая Wx = W%, найдем
ал о
и приближенно выражение, учитывающее силу сухого трения в урав нениях (1.21) и (1.22), примет вид
±<2 = - ^ ф ( 0 . |
(Д) |
ала» |
|
Произведя замену в уравнении (1.21) выражений для гистерезис ного сопротивления и постоянной силы трения их значениями (г) и (д), получим
Л4ф (i) -f Ry (t) + |
m (t) + |
ф (t) + Ncp (0 = F cos сot |
|
40 |
|
ажв |
|
или |
|
|
|
M v + ( R + ^ + - ^ \ < f + N(f = Fcosa>t. |
(1.23) |
Поскольку все функции теперь относятся к одному и тому же мо менту времени, подчеркивание их зависимости от времени можно опу стить.
Поделив последнее уравнение на М и введя обозначение
2р = 2г |
у.Х2 |
4Q |
|
СО |
аМла» |
||
|
получим
Ф + 2рф -j- ^2ф = f cos a t.
Здесь по-прежнему
г, R |
\2 N |
t F |
, ч |
2г = — , |
№ = — , /■= — , а |
р —по выражению (е). |
|
Л4 |
М |
М |
|
(е)
(1.24)
Интеграл уравнения (1.24) состоит из общего интеграла такого же уравнения, но без правой части
Ф + 2рф + ^2ф = 0, (ж)
который отображает собственные затухающие колебания, и из част ного интеграла уравнения (1.24). Общий интеграл зависит как от на чального отклонения массы и ее скорости в момент t = О, так и от на чальной фазы возмущающей силы. Частный интеграл определяет вы нужденные колебания.
Поскольку собственные колебания рано или поздно затухнут, нас интересует чаще всего лишь вынужденные колебания, т. е. только частный интеграл уравнения (1.24).
31
Будем его искать в форме
<р = A cos (ot-j-B sin at, |
(з) |
где А, В — некоторые постоянные, подбираемые так, чтобы решение удовлетворяло уравнению (1.24).
Подставив (з) в (1.24), получим
—Л со2 cos соt —Boo2 sin соt —2рсоЛ sin at +
+2 р а В cos соt + Я2Л cos a t + Я2В sin соt = f cos сot.
Это уравнение должно быть справедливо в любой момент времени. Подбирая моменты времени так, чтобы либо sin со^, либо cos со / равня лись нулю, получим два уравнения для определения А и В
(Я2— со2) А + 2рсоВ = /; |
|
|
|
|||
—2рсоЛ + (X2— со2) В = О, |
|
|
|
|||
откуда |
|
|
|
|
|
|
(Я2 — to3)f |
Q ________2ры/_____ |
(и) |
||||
(Я2 — со2)2 + (2рш)2 |
’ |
(Я2 — со2)2 + (2рц>)2 |
||||
|
||||||
и, следовательно, |
|
|
|
|
|
|
(Я2 — со2) f |
|
2р<а/ |
|
|
|
|
ф = |
|
(Я2 — со2)2 + |
(2рш)2 |
■sin at. |
|
|
(Я2 — и2)2 + (2р<о)2 - cos соt - |
|
|||||
Синусоиду и косинусоиду одного |
аргумента |
можно объединить |
в косинусоиду с амплитудой, равной геометрической сумме амплитуд объединяемых зависимостей
Ф= |
| / Л 2 + |
|
f cos (at— у) |
||
В2 cos(co/—у) = |
|
(2pw)2 |
|||
|
|
|
V (Я2 — со2)2 + |
||
Поделив числитель |
и знаменатель последнего |
выражения на Я2 |
|||
и заменив 2р его выражением (е), |
получим |
|
|
||
|
Ф = аСт |
cos (at — у) |
|
(1.25) |
|
|
|
2 _j_ (^г<л |
|
||
|
|
|
4Q \2 |
||
|
|
|
I Я2 |
nNa) |
|
/ |
F |
|
|
т. е. |
отклонение, ко- |
где а„ = -1—= |
------- статическое отклонение, |
||||
Я2 |
N |
|
|
|
|
торое масса имела бы под действием силы F при статическом ее при |
|||||
ложении и жесткости системы N, т. е. без учета динамики явления. |
|||||
При отсутствии сопротивлений |
|
|
|
||
|
|
F |
cos at |
|
( 1.26) |
|
|
N ' |
j _<о2 |
|
|
|
|
|
|
||
|
|
|
Я2 |
|
|
32
Угол отставания вынужденных колебаний от фазы возмущающей силы определяется зависимостями:
cosy =
2рсо \ а
|
|
|
X*2 ) |
|
|
2 рсо |
|
(1.27) |
|
|
|
|
||
sm у- |
X2 |
|
|
|
\2 |
( |
2 рсо |
||
|
||||
|
/(■- X2 |
\ |
X2 |
либо первой из этих формул при условии, что
0 < у < 180°.
Решение уравнения (1.22) для случая, когда амплитуда силы воз
растает пропорционально |
квадрату частоты о, может быть получено |
|||
из выражения (1.25) путем замены |
р |
|||
в нем множителя аст= — мно- |
||||
|
\а |
|
|
|
С00 |
F 0a 2 |
|
|
|
жителем |
М ар ? |
|
|
|
N |
|
|
|
|
При этом вместо (1.25) получим |
|
|||
|
|
|
а* |
|
|
|
|
1? |
cos (to/ —у), (1.28) |
|
/ F ? |
а |
|
|
|
(2га |
4Q \ а |
||
|
+ |
\ X2 |
nNa ) |
|
где |
|
|
а°° |
(1.29) |
|
|
|
Ма^. |
— амплитуда колебания при неограниченном возрастании частоты возмущающей силы (о оо).*
Фаза отставания колебания от возмущающей силы определяется теми же зависимостями (1.27).
В случае, если необходимо получить не только установившиеся вынужденные колебания, но изучить и начальный период движения,
следует к частному интегралу уравнения |
(1.24), т. е. к выражению |
||
* |
Легко видеть, что при а -> оо выражение |
СО* |
стремится к еди- |
------ — |
|||
нице. |
|
X,2 / . . . |
|
|
|
|
|
2 В. |
В. Давыдов, Н. В. Маттес |
|
33 |
(з), прибавить общий интеграл уравнения (ж), имеющий вид (1.7). Выполнив это, будем иметь
Ф —e~pt (Dxc o s + |
D2 sinX]/)-f A cos соt-\-B sin at, |
|
||
где |
|
______ |
|
|
Найдя постоянные D x и D %, получим |
|
|
||
Ф =>=e~pt ^ф0 cos ^ + Р<Ро^~ фо sin j — |
|
|||
— e~pt ^4 cos X±t + |
Ap |
Ba sin |
+ A cos at + В sin at. |
(к) |
Первый член правой |
части соответствует свободным колебаниям, |
т. е. колебаниям при отсутствии возмущающей силы. При нулевых
начальных условиях ф0 = Фо = 0 эти колебания отсутствуют. Второй член отображает сопровождающие колебания. Они зависят
от амплитуды и частоты возмущающей силы и не зависят от начальных условий. Они происходят с собственной частотой и всегда в начале движения сопровождают вынужденные колебания. Так же, как и сво бодные колебания, они сравнительно скоро затухают. Последние два слагаемых соответствуют вынужденным колебаниям.
Качественный анализ явления. Вынужденные колебания происхо дят всегда с частотой возмущающей силы со. Это основное свойство вы нужденных колебаний позволяет при экспериментальном обследовании вибрации установить ее причину и разыскать ее источник. Вынуж денные колебания — незатухающие, их амплитуды неизменны в тече ние. всего времени действия гармонической возмущающей силы.
Амплитуду вынужденных колебаний, т. е. множитель в выражении (1.25) при функции косинуса при наличии силы сухого трения
/ |
_____аст |
|
(1.30) |
|
а1\2 /2га |
х-+ Щ ' |
|||
ХЧ + 1 |
Я2 |
|||
л Na |
||||
|
|
|
следует рассчитывать методом последовательных приближений. Сна чала нужно положить Q = 0 и вычислить по (1.30) амплитуду а; най денное значение следует подставить в правую часть (1.30) и найти но вое значение а (второе). Найденное (второе) значение снова подставить в правую часть, найти третье значение и продолжить этот процесс, пока вычисленное значение будет мало отличаться от подставленного.
Аналогично следует поступать и при вычислении амплитуды для колебаний в случае квадратичной зависимости возмущающей силы от частоты
а — |
( 1.31) |
4Q \а
лЫа }
34
Если сухое трение отсутствует Q = 0, то |
амплитуды колебаний |
по формулам (1.30) или (1.31) находятся сразу. |
Если, кроме того, от |
сутствует и гистерезисное сопротивление, то из выражения (1.30) по лучим
а = аг |
0)2 |
|
(1.32) |
|
\2 |
2га) |
|
V |
'~№ |
+ |
У? |
Первый множитель последнего выражения при неизменных F и N есть величина постоянная.
Второй множитель
1
(1.33)
носит название динамического множителя или коэффициента дина мичности и представляет собой поправочный коэффициент, на кото рый нужно помножить статическое отклонение, чтобы получить амплитуду вынужденных колебаний. Динамический множитель—функ ция трех величин: собственной и вынужденной частот Я. и со и коэффи циента г, учитывающего вязкое сопротивление.
При низких частотах возмущающей силы, когда со X, действие силы близко к статическому, т. е. амплитуда близка к аст, а коэффи циент динамичности — к единице.
При частотах возмущающей силы, близких к собственной частоте со ~ %, наблюдается резкое увеличение амплитуд вынужденных ко лебаний, так называемое явление резонанса, играющее важную роль в динамических расчетах.
Точное значение резонансной частоты, при которой и амплитуды, и динамический множитель достигают максимальных значений, мо
жет быть получено, если производную подкоренной величины по —
приравнять нулю. При этом можно получить |
/ V |
|
|
||
сор = Х |
|
(1.34) |
Из формулы (1.34) видно, что |
резонансная частота |
несколько |
меньше К. |
Я возмущающая сила |
конечной |
При очень высоких частотах со |
амплитуды настолько быстро изменяет свое направление, что система «не успевает» за ней следовать, и амплитуды колебаний оказываются
очень незначительными: при со -» |
оо и динамический множитель ц, |
|||
и амплитуды а стремятся к нулю. |
зависимости |
амплитуд колебаний |
||
На рис. 10, а показан график |
||||
от частоты возмущающей силы со |
для частного |
случая X = 39 с-1, |
||
а аст = 1 мм и разных |
значений сопротивления г. |
При выборе соот |
||
ветствующего масштаба |
кривые на |
рис. 10, а дают |
представление и |
об изменении динамического множителя.
2* |
35 |
Следует отметить, что при сопротивлении г > - р = кривые ампли-
У
туд и динамического множителя монотонно убывают с возрастанием о (см. нижнюю кривую при г = 32 с-1).
Явление резонанса наблюдается и в случае, когда амплитуда воз мущающей силы растет пропорционально квадрату частоты. В этом случае и при наличии лишь вязкого сопротивления (г =j=0, к = Q =0)
Рис. |
10. Зависимость амплитуд колебаний |
от час |
|
тоты |
возмущающей силы |
при %= 39 с—1 |
и раз |
|
личных сопротивлениях |
|
|
из формулы (1.31) следует |
о-1 |
|
|
|
|
|
|
|
’ |
А* |
(1.35) |
|
|
|
Характер зависимости динамических амплитуд от частоты возму щающей силы со для частного случая А, = 39 с-1, аот = 0,625 мм и различных сопротивлениях показан на рис. 10, б.
• В этом случае формула. (1.34) дает мнимую резонансную частоту.
36
Максимальных значений амплитуда колебаний достигает при ©pt несколько больших Я
(0„ |
(1.36) |
Итак, в случае, когда амплитуды возмущающих усилий пропорциональны квадрату возмущающей частоты, амплитуды колебаний возрастают от нуля до некоторого максимального значения (резонанс), а после резонанса уменьшаются, стремясь при неограниченном воз растании со к определенному пределу аоо.
Сказанное, разумеется, справедливо при относительно небольших
X
сопротивлениях г < -^-=, когда формула (1.36) дает вещественное зна
чение ©р. При более высоком сопротивлении г амплитуды колебаний монотонно возрастают по мере возрастания со и достигают максималь ного значения а при со = оЪ.
Как видно из формул (1.25) и (1.28), отклонения при колебаниях
несколько отстают по фазе от возмущающей силы. |
Угол отставания у |
определяется формулами (1.27). |
Я угол этот, как |
При небольших частотах возмущающей силы со |
видно из выражений (1.27), невелик, т. е. фазы возмущающей силы и отклонения почти совпадают: масса отклоняется в ту же сторону, в ко торую действует сила, и, когда сила достигает своего наибольшего значения, наибольшего значения достигает и отклонение.
При резонансе со ~ Я из (1.27) следует, что у ~ 90°. Нетрудно понять, что в этом случае с фазой возмущающей силы совпадает фаза скорости, т. е. в каждый момент времени направления возмущающей силы и скорости совпадают, сила все время стремится увеличить ско рость и, как результат этого, отклонения оказываются весьма боль шими.. Силы инерции и силы упругости находятся в противофазе и взаимно уравновешивают друг друга.
При очень большой возмущающей частоте м Я в фазе с силой находится ускорение, отклонение же отстает на полпериода (у ss 180°) и находится в противофазе с силой: масса отклоняется в сторону, противоположную направлению силы.
В промежуточных случаях опоздание по фазе зависит как от со противления, так и от близости к резонансу. Чем меньше сопротивле ние и чем дальше возмущающая частота © от собственной Я, тем ближе ■у к своим крайним пределам: 0 и 180°. При отсутствии сопротивления у — 0, если ©<Я, и у = 180°, если ©>>Я. Чем больше сопротивление и чем ближе © к Я, тем ближе у к резонансному опозданию, т. е. к 90°.
При гистерезисном сопротивлении амплитудно-частотные харак теристики (см. рис. 10) и явление прохождения резонанса будут та кими же. Значения амплитуд колебаний для этого случая можно по
лучить из (1.25), (1.28), |
положив в них г = |
Q — 0 |
а = |
а СТ |
( 1.37) |
37
|
(Г)2 |
|
а = а |
U |
(1.38) |
|
где ас1 и ат имеют прежние значения. Наибольших значений ампли
туда, вычисленная по формуле (1.37), т. е. при резонансе, имеет место, очевидно, при сор = к. Применительно к формуле (1.38)
|
|
(1.39) |
П р и м е р . |
Исследовать |
амплитуду вынужденных колебаний массы М = |
= 1000 кг = 1 |
кН -с2/м [0,102 |
тс-с2/м] на упругой невесомой балке (см. рис. 7) |
при различных частотах в предположении, что на указанный груз действует вертикальная сила, пропорциональная квадрату частоты, с амплитудой 1 кН
[102 кгс] при частоте |
40 |
с- 1 . Коэффициент сопротивления т принять равным |
|
4 с- 1 , а жесткость N == 1536 кН/м [157 тс/м]. |
|||
Собственная частота |
без |
сопротивления X — 39,2 с-1 была определена |
|
раньше (см. с. 22). Амплитуда |
колебаний при неопределенно большой частоте |
||
по формуле (1.29) |
|
|
|
я = |
M(Oq |
= 0,000625м = 0,625 мм. |
|
00 |
1- 403 |
Наибольшие амплитуды колебаний, согласно формуле (1.36), имеют место
при
39,2
:39,4с- 1 .
2-42
/39,22
Амплитуды колебаний для нескольких частот вычисляем в таблич ной форме по формуле (1.35) (табл. 1).
|
|
|
|
|
|
Таблица 1 |
|
Определение амплитуд колебаний |
по формуле (1.35) |
|
|||
|
|
|
|
|
|
Амплитуда |
|
ш2 |
и 2 |
2г(О |
|
|
(0а |
и , с- 1 |
|
|
|
|||
X2 |
1— V |
— =0.00527 |
У . . . |
a = flo o _ ^ _ |
||
|
|
|
||||
|
|
|
|
|
|
Y77. |
|
|
|
|
|
|
ММ |
0 |
0 |
1 |
0 |
|
1 |
0 |
20 |
0 ,2 6 3 |
0 ,7 3 7 |
0 ,1 0 5 |
|
0 ,7 4 5 |
0 ,2 2 |
30 |
0 ,5 9 2 |
0 ,4 0 8 |
0 ,1 5 8 |
|
0 ,4 3 8 |
0 ,8 4 |
’ 39 ' |
1 ' |
0 ' |
0 ,2 0 6 |
|
0 ,2 0 6 |
3 ,0 3 |
ОО |
— |
— |
— |
|
— |
0 ,6 2 5 |
Аналогично вычислены и остальные амплитуды, по которым по строена верхняя из кривых, изображенных на рис. 10, б.
38
Амортизация возмущающих усилий. Работающий поршневой дви гатель вследствие своей неуравновешенности и непостоянства крутя щего момента прилагает к своему фундаменту периодически изменяю щиеся усилия, которые могут вызвать нежелательные колебания кор пуса судна.
Для уменьшения этого вредного явления иногда применяется уста новка двигателей на амортизаторы — упругие связи, помещаемые между станиной двигателя и судовым фундаментом (рис. 11).
При учете лишь вязкого сопротивления для колебаний двигателя
формула (1.28) |
дает |
|
св* |
Ф |
cos (со?—у). |
1 ♦ t
77/7777777777777777Т/
Рис. 11. Установка двигателя на амортизаторе с вяз ким сопротивлением
Усилие в пружинах амортизатора будет пропорционально откло нению Л/'ср, а сила вязкого сопротивления — скорости /?ф. Суммарное усилие, передающееся фундаменту, составит, учитывая, что N = МК2
и R = 2гМ ,
Щ + R ^ F o ^ J |
c o s ( a t — у) |
|
(-£ № )’ |
||
Y |
||
2rcp |
||
- F a _© \* |
s in (со/ — у ) |
|
X* |
co0
Функции синуса и косинуса разных амплитуд, но одного аргу мента могут быть объединены в одну функцию синуса или коси нуса с амплитудой, равной геометрической сумме амплитуд состав ляющих гармоник
39