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

книги / Механика деформируемого твердого тела.-1

.pdf
Скачиваний:
8
Добавлен:
20.11.2023
Размер:
15.15 Mб
Скачать

ниях как вследствие трения внутри материала, так и из-за конструкционного трения, возникающего, например^ между валом и ступицей диска, в цанговых зажимах и т. д.

Имея в виду изучить здесь влияние сил трения обо­ их типов, воспользуемся припятой в § 30 вращающейся координатной системой (см. рис. 30.4) ,и составим диф­ ференциальные уравнения относительного движения ро­ тора. Силу внешнего трения по-прежнему примем -про­

порциональной абсолютной

скорости центра масс v =

= ve + vr (ve — переносная

и vr — относительная скоро­

сти точки С) :

 

 

—b\—b(\e + vr) = —Ь[(—сот] + "|)i +

(со| + rj)j].

Учтем также силу внутреннего трения,

которая связана

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

-Ь 0\г = - b 0(li+ rjj),

где

bo — коэффициент вязкого

внутреннего трения.

 

С

учетом

приведенных

выражений

получим

вместо*

(30.6)

 

 

 

 

 

 

 

 

 

 

 

Ï — 2соц + (к2— со2) £ + 2 (п + По) g — 2тгсот] =

к2е,

 

 

 

 

 

 

 

 

 

 

(32.1)

 

 

г] + 2cog + (к2 — со2) г] +

2(п +

тг0)т] + 2дсо£ =

0,

 

 

где

п9= b j (2т) , п== Ъ/ (2т) .

покоя

соответствует

ре­

Состоянию

относительного

шение

 

 

 

 

 

 

 

 

^

 

к2е (к2— со2)

 

 

2к2епю

 

/п0

ё*

=

(к2- го2) 2 + 4 n V ’

 

(л2 - ( 0 2)2 + 4й2ш2 ’

(

>

которое отличается от решения* (30.7)

членами,

выража­

ющими влияние внешнего трения. Параметр щ сюда но входит/ так как силы внутреннего трения при относи­ тельном покое отсутствуют.

Рассматривая возмущенное движение I = £* + б£ = щ + 6г], придем к однородной системе дифференци­ альных уравнений для б| и бг), а затем — к характери­ стическому уравнению

Я4 + BsV + В2Х2+ В,К + Во = 0,

(32.3)

г д е

В3 = 4 (га + щ), i?2 = 4 (га + га0)2 + 2 (к2 + со2),

= 4[(&2— со2) (га + га0) + 8ю2га],

(32.4)

В0 = (А:2 — ю2)2 + 4га2ш2.

Для устойчивости рассматриваемой системы необходимо и достаточно, чтобы вещественные части всех кор­ ней .уравнения {32.3) были отрицательными, что соот­ ветствует выполнению условий Рауса — Гурвица, кото­ рые имеют вид (см. выше (18.12) в § 48)

В0> 0 ,

5 t >

0 ,

В2> О,

В3> О,

 

а д я

з -

В\ -

В0Ъ1 >

о..

(32.5)

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

случае

(п Ф0,

п0= 0) независимо от значений

ш и

к

выполняются все условия

Рауса — Гурвица, ‘ т.

е.

вал

устойчив (даже асимптотически устойчив!) при любых

угловых скоростях вращения. Во втором случае

(п = 0,

п0Ф 0)

характеристическое уравнение принимает

вид

 

 

(К2 + 2п 0К + к 2 -

ш 2) 2 + ( 2 о Ж ) 2 =

0 .

 

( 3 2

Первая

группа

условий

(32.5)

выполняется при лю­

бых значениях ш и ку и для устойчивости необходимо

выполнение лишь последнего условия (32.5), которое сво­

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

 

 

(32.7)

 

 

 

 

(ù<k.

 

 

Таким образок, при отсутствии внешнего трения си­

стема

устойчива

только

в

д о к р и т и ч е с к о й

о б л а ­

сти,

а в закритической

области

внутреннее

трение вы­

зывает

неустойчивость.

Если ш >

/с, то характеристиче­

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

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

работу; работа

этих сил, конечно,'

отрицательна, но

они — именно

они — создают условия

для перекачки

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

который обеспечивает предположенное выше постоянство угловой скорости; эта работа расходуется на увеличение энергии ротора и на потери, связанные с действием сил внутреннего трения. На эти обстоятельства обратил вни­ мание В. В. Болотин в книге [13].

Напомним, что условие устойчивости в виде (32.7) соответствует отсутствию внешнего трения. При одно­ временном учете внешнего трения порог возникновения неустойчивости отодвигается, и с помощью (32.5) можно получить условие устойчивости в выразительной форме

(32.8)

Иначе обсуждается вопрос о влиянии внутреннего трения на устойчивость вращающихся роторов в получив­ шей всемирное признание книге С. П. Тимошенко [65], автор которой славится мастерством простого и ясного анализа самых сложпых механических явлений. Верный своему принципу, С. П. Тимошенко и в рассматриваемом

случае

предельно

упростил зада­

 

чу. Он принял, что сила внутрен­

 

него трения

просто

пропорцио­

 

нальна прогибу вала F аг и на­

 

правлена

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

пло­

 

скости

изгиба

(рис. 32.1),

так что

- с г

проекции

этой

силы

на

оси

не­

подвижной

координатной

системы

 

соответственно

равны

Fx = —ау,

 

Fy= ах

— коэффициент

про­

Рис. 32.1. Схема Тимо­

порциональности) .

Дифференци­

шенко

альные

уравнения

возмущенного

 

движения центра вполне уравновешенного диска С. П. Ти­ мошенко записал в виде

тх + сх + ау = 0, ту + су ах = О

и. из простого анализа этих уравнений получил, что ро­ тор всегда неустойчив. Сразу настораживает, по-видимо­ му, незамеченное С. П. Тимошенко простое обстоятель­ ство: к тому же выводу о неустойчивости можно прийти, изменив направление силы внутреннего трения па обрат­ ное, когда дифференциальные уравнения записываются в виде

тх + сх — ау = 0, ту + су + ах = 0.

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

вости ротора не только в закритической области, но так­ же в докритической области, и даже при отсутствии вра­ щения. В данном случае стремление к упрощению по­ становки задачи оказалось чрезмерным, а границы ра­ зумной схематизации были явно нарушены*). Изящное высказывание Пуанкаре: «Наши задачи мы должны ре­ шать сГо всей возможной простотой, но, не проще того» следует относить не только к процедуре решения, но и к самой постановке задач.

§ 33. Влияние внутренней и внешней упругой анизотропии

Чаще всего понятие «анизотропия» связывается со* свойствами материалов. В прошлом обычным примером анизотропного материала служило дерево, механические свойства которого вдоль и поперек волокон существенно* неодинаковы; широко применяемые в современной тех*- нике композитные материалы чаще всего также анизо­ тропны. Но помимо «физической» анизотропии иногда приходится иметь дело с системами, которые естественно называть «конструктивно» анизотропными.

Рис. 33.1. Анизотропия: а) внутренняя; б) внешняя

Так, например, всюду выше было принято,, что изгибные жесткости вала одинаковы во всех направлениях. То же было приписано упругим опорам в варианте мо­ дели Феппля, показанном на рис. 29.3. Однако в ряде' случаев конструкция не обладает этим свойством. Та­ ковы, например, роторы электрических машин, имеющие1

продольные

прорези

для размещения обмоток

(рис.

33.1, а),

или

жесткие

роторы, установленные

на

анизо­

тропные

упругие опоры, имеющие, различные

жесткости

*) Ошибка С. П. Тимошенко была отмечена в работе Э. Л. Позняка [47].

т вертикальном и горизонтальном направлениях (рис. 33.1,6). Можно сказать, что обе эти конструкции анизотропны, но в первом случае главные оси жесткости в р а щ а ю т с я вместе с ротором, а во втором случае они остаются н е п о д в и ж н ы м и относительно основания системы.

Свойства систем этих двух типов существенно неоди­ наковы, и поэтому следует различать внутреннюю и внешнюю анизотропию. (Иногда пользуются также стран­ но звучащими терминами «подвижная анизотропия» и «неподвижная анизотропия», которые нельзя назвать удачными; пока еще не поздно и они не успели проч­ но внедриться в лексикон, от них следовало бы поско­ рее отказаться). Встречаются такжр смешанные случаи, когда внешняя анизотропия сочетается с внутренней.

Впервые динамику внутренне анизотропных роторов исследовал Прандтль *) в 1918 г. В работе Прандтля рассматривается модель Феппля, но с усложнением — считается, что вал обладает различными главными жест­ костями при изгибе Ci и с2. Движение изучается во

вращающейся координатной системе, причем

ось

на­

правляется параллельно первой

главной

оси

жесткости,

а ось Ог\ — параллельно второй

главной

оси.

Для

упро­

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

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

1-2(4) + (& ? -2H = 0,

г) + 2 (4 + (kl — о 2) г) = 0

*) Людвиг Прапдтль (1875—1953)— профессор механики в Ганновере и Геттингене. Известен прежде всего как один из осно­ вателей современной, аэродинамики, но Л. Прандтлто принадле­ жит ряд результатов также в механике твердого деформируемого чела (исследование устойчивости плоской формы изгиба, разработ­ ка мембранной аналогии в теории кручения упругих стержней

и др.)«

(величины к\ = cjrrt и к\ = с21пь представляют собой соб­ ственные частоты невращающегося вала при его колеба­ ниях в направлениях главных осей жесткости). Далее*

полагая

1, = Ауеи

и ц = А 2ех\ приходим к

однородной

системе

алгебрайческих уравнений для A t

и А2,, а за­

тем — из

условия

существования ненулевых

решений —

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

Я4 + Я2 (к\ + к\ + 2со2) -f- (kl — о)2) (к\ г— а>2) == 0.

Решив это уравнение относительно Я2, получим

X2=4 [ - (*ï + *5 + 2®2)± V (*5 - **)“ + 8ш2(*1 + А|)].

(33.1>

Так как подкоренное выражение положительно, та оба найденных значения Я2 вещественны. Следовательно* соответствующие значения Я либо чисто мнимые, либо вещественные; комплексных значений Я в этой задача нет. Мнимым значениям Я соответствуют колебания с по­ стоянными амплитудами, а вещественным отрицательным значениям Я —. неколебательный затухающий процесс. В обоих этих случаях вал может быть признан устой­ чивым.

Однако, если какой-либо из в е щ е с т в е н н ы х кор­ ней Я оказывается п о л о ж и т е л ь н ы м , то с течением времени отклонения | и ц будут неограниченно возра­ стать, что свидетельствует о неустойчивости вала. Для положительности одного из корней Я необходимо, чтобы

выполнялось

неравенство

 

 

 

{к\ -

kl) + 8(Ù2 (ftï +

kl) > {k\ + k\ +

2co2)2,

или эквивалентное неравенство

 

 

 

(ю2 -

kf) (to2 -

kl) < 0.

(33.2)

Отсюда видно, что

если

ki <

ю < к2% то

система н е-

у с т о й ч и в а

(теорема Прандтля) .

 

Таким образом, различие главных жесткостей вызы­ вает появление конечной зоны неустойчивости, соответ­ ствующей интервалу частот [ки к2]; при ki = к2 эта зона вырождается в единственную опасную точку ю = к.

Примечательный случай соответствует бесконечна большой жесткости в одном из главных направлений.. Если, например, с2-*• «>, То зона неустойчивости начина­ ется при о) == kt и простирается до бесконечности. Суще­

ство этого результата легко понять, обратив внимание на то, что при с2->• °° вал может изгибаться только в одной из главных плоскостей, т. е. исчезает одна степень сво­ боды и ротор приобретает свойства системы, показанной на рис. 30.3 (как мы писали в § 30, при со > к она не­ устойчива).

При учете внешнего трения вместо неравенства (33.2) получается более сильное неравенство

(со3 - kl) (ю3 - Al) + {2ruùf < 0,

(33.3)

которое определяет и более узкую зопу неустойчивости.

При достаточно большом

трении, когда

 

п >

2\к%— /cj,

(33.4)

неравенство (33.3) вообще не может

быть выполнено,

т. е. неустойчивость, связанная ; с внутренней анизотро­ пией, полностью подавляется внешним трением.

Случай в н е ш н е й анизотропии естественно и удобно анализировать в неподвижной координатной системе. Есци пренебречь трением, то дифференциальные урав­ нения движения можно сразу получить из (29.2), введя, конечно, различие между коэффициентами жесткости ci и с2 в двух главные направлениях, параллельно которым ориентируем координатные оси. После этого вместо (29.3) получится при п = 0:

х + к\х = к\е cos ю£,

(33.5)

У + к\у = к\е sin оо£,

где к\ = cjm, к\ = с2/пг. Стационарный режим описыва­ ется частными решениями уравнений (33.5)

х — АI cos Ü)£, ï/ = 4 2sin(ü£,

(33.6)

в которых

(33.7)

Очевидно, что в данном случае критические состояния возникают при со = ki и со = к2.

Если из (33.6) исключить время, то мы найдем траекторию точки С — эллипс, по которому точка С

обращается вокруг центра

О:

 

 

 

(33.8)

Отметим, что хотя п е р и о д

обращения центра масс

вокруг центра О равен

периоду

собственного в р а щ е ­

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

На рис. 33.2 в качестве иллюстраций показаны тра­ ектории и направления движения центра масс для трех

Ж

Рис. 33.2. Примеры движения в условиях внешней анизотропии; в случаях а) и в) — прямая прецессия, б) — обратная прецессия

значений угловой скорости вала в характерных диапазо­

нах: о < &4, ki < со < к2 и

со > к2. Анизотропия вала ха­

рактеризуется

отношением

собственных частот

k2/kt = 2.

Рис. 33.2, а

относится

к

случаю,

когда

с6/кх=

0,5У

рис. 33.2, б — к случаю co/fti = 1,5

и

рис. 33L2, в — к

слу­

чаю o)/&i = 2,5.

 

второй

случай (рис. 33.2, б):

Наиболее

примечателен

направление

обращения

п р о т и в о п о л о ж н о направ­

лению вращения вала. Если вглядеться в соотношения (33.7), то станет ясным, что это явлепие (обратная пре­ цессия) связано с тем, что величины А х и Az имеют различные знаки, и возникает всегда, если угловая ско­ рость вращения располагается в промежутке между соб­ ственными частотами (Лг4< со < Лг2) . Отметим, что эти же неравенства нам встречались в задаче о роторе с внут­ ренней анизотропией, где, однако, они определяют не

условия обратной прецессии, а согласно теореме Прандтля — область неустойчивости *).

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

но

жесткими

направляющими, соз­

 

 

дающими

внешнюю

анизотропию

 

 

(рис_33.3). В этом случае система

 

 

имеет

только

одну степень свободы

 

 

и колебания возможны только в пло­

 

 

скости,

фиксированной

направляю­

 

 

щими. Если

сi и с2— коэффициен­

 

 

ты жесткости ротора в собственных

 

 

главных плоскостях, то коэффициент

 

 

жесткости системы в плоскости ко­

Рис. 33.3. Особый слу­

лебаний

определяется

выражением:

с =

Ci cos2 cot + съsin2 £,

т. e. оказы­

чай

сочетания внеш­

ней

и внутренней

вается явной функцией времени. Со­

 

амортизации

ответственно,

дифференциальное

 

 

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

У + т ( 1 + С' С* sin2(ûi) У= 0

п описывает параметрические колебания. Здесь также возможна неустойчивость (параметрический резонанс) , но границы областей неустойчивости определяются не­ сколько сложнее, чем в рассмотренных выше случаях. Решение этой задачи было дано В. В. Болотиным (см. его книгу [12J).

§ 34. Влияние нелинейной упругости

Нелипейность

упругих свойств

роторов возникает,

в частности, когда

в конструкцию

опор входят шарико-

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

P = Kyzfz.

(34.1)

*) Внешняя анизотропия — не единственная причина возник­ новения устойчивой обратной прецессии в роторных системах; то же явление возникает из-за гироскопического эффекта при пере­ косах плоскости диска.

Рис. 341. Упругий шарик между обоймами

В записанном здесь законе Герца *) коэффициент К оп­ ределяется геометрией контактирующих тел и упругими постоянными их материалов. Например, при сжатии уп­ ругого шара по схеме, изображенной на рис. 34.1, К =

= —

при одинаковых материалах), где Е и v —

3 (1 — v*)

 

соответственно модуль упругости и коэффициент Пуас­ сона. В данном случае у — сближение тел А и В (точнее говоря — сближение точек этих телг достаточно удаленных от контакт­

ной зоны).

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

Соотношение (34.1) получено Герцем в результате решения стати­ ческой задачи теории упругости и безупречно в той же степени, в ко­

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

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

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

*) Генрих Рудольф Герц (1857— 1894) — профессор физики в Карлсруэ и Бонне. Экспериментально доказал существование элект­ ромагнитных волн, открыл внешний фотоэффект и др. «Хотя тео­ рия упругости занимала в научных достижениях Герца сравни­ тельно скромное место, мы все же обязаны ему решением ряда трудных проблем, представлявших к тому же, и большое практи­ ческое значение» [64, с. 417].