Прикладные задачи
^^^^^^^^^^»нелинейной теории колебаний и вслн
УДК 531.47+531.536+531.555.1
О ДИНАМИКЕ ВРАЩЕНИЯ ТВЕРДОГО ТЕЛА ВОКРУГ НЕПОДВИЖНОЙ ОСИ, ПРОХОДЯЩЕЙ ЧЕРЕЗ ЦЕНТР МАСС, ПРИ СУХОМ ТРЕНИИ В ОПОРАХ
В. В. Чистяков ФГБОУ ВПО Ярославский ГПУ им. К.Д. Ушинского
Рассмотрена динамика вращения твердого тела вокруг неподвижной оси Oz, проходящей через его центр масс C и не являющейся главной осью инерции. Возникающие инерционные пары Mx,y (Jxz, Jyz, ю, е) обуславливают сухое трение с коэффициентом S, делающее величину углового ускорения е и причиной, и следствием действия сил одновременно. Проинтегрированы аналитически уравнения свободного и вынужденного (Mz = const) движений с учетом квадратичного по угловой скорости ю аэродинамического сопротивления, а также численно - колебательного движения под действием гармонического и ангармонических («ступенька», «пила») моментов. Полученные предельные циклы сравниваются с аналогичными при линейных вращательных колебаниях.
Ключевые слова: Астатический ротатор, ось инерции, инерционные пары, сухое трение, квадратичное сопротивление, вращательные колебания.
Введение
Если ось вращения проходит через центр масс насаженного тела, но не является для него осью свободного вращения из-за отличных от нуля центробежных моментов инерции <1хг и ,1уг, то поперечные реакции N0, Кв в опорах (рис. 1) и обусловленный ими тормозящий момент сухого трения Mfr появляются лишь при вращении тела и зависят от закона вращения. Потому явление самоотторжения (заклинивания) оси, как это бывает, например, в случае маятника с трением [1], для астатического ротатора невозможно в принципе.
В качестве примера такого ротатора можно привести знаменитый кельтский камень [2], насаженный на поворотную ось 2-го порядка, или нечто, изготовленное по его образцу, например, просверленный вкось через центр однородный эллипсоид с неравными осями.
Для нахождения необходимых сил нормальных реакций в местах контактов тел в самом общем случае движения применяют изложенный в § 7.2 монографии А.И. Лурье [3] метод для определения реакций связи контактирующих
А
в т точках твердых тел через введение т избыточных обобщенных координат дп+1,Яи+2, ■ ■ ■ Чи+т- После составления уравнений Лагранжа эти координаты вкупе с обобщенными скоростями принимаются равными нулю.
со, е
Мх
О
с
Так, автор [1] использует данный прием для нахождения выражения для силы нормальной реакции N и самой силы трения Ф в оси колеблющегося маятника с целью последующего составления и решения уравнения нелинейных колебаний с возможностью уже упомянутого самоотторжения, когда маятник «зависает».
Рис. 1. Астатический ротатор с неглавной осью вращения (ЗО-кельтский камень)
Нетрудно показать, что применение данного метода через введение избыточных угловых переменных д2 и
ответственных за поворот вокруг направлений Сх и Су, [3] даёт то же уравнение динамики, что и использование принципа Д'Аламбера для нахождения поперечных реакций N0, Nв (см. рис. 1) в осях вращающегося тела, а с ними и момента М^.
Последний при этом будет зависеть не только от угловой скорости ю, но и от углового ускорения е (!), и при этом вкупе с прочими парами сил определяет само это ускорение. То есть формально величины М^ и е выступают одновременно и ведущими и ведомыми факторами.
Эта двусторонняя причинно-следственная связь не может не вызывать фундаментального интереса с точки зрения динамики нелинейной, так как в некоторых вырожденных случаях может приводить к несовместным в рамках динамики абсолютно твердого тела (а.т.т.) уравнениям, внося тем самым вклад в копилку парадоксов сухого трения [4].
Кроме того, подобная связь часто реализуется в социодинамике (эконофизике), динамике биологических процессов и явлений и т.д., и потому есть смысл исследовать возможный механический аналог столь сложных систем.
Равно как и само незаклинивающее сухое трение, не существующее вне движения и чем-то напоминающее трение жидкостное, представляет чисто практический интерес рассмотрение ситуаций, когда требуется чрезвычайно точное динамическое уравновешивание высокоскоростного вращательного устройства (турбины, ротора, вала станка и т.д.), а также в тех случаях, когда такое уравновешивание не предполагается. Последнее, как ни печально, может означать аварийную, а в некоторых случаях - катастрофическую ситуацию утраты гидротурбиной, ротором и пр. динамической уравновешенности, когда и вступает в силу вышеозначенное трение и не только оно (см. далее).
Настоящая работа посвящена изучению вопроса, как влияет такое трение в осях на динамику вращения ротатора, насколько она схожа или различна с динамикой, описываемой стандартными уравнениями с линейным либо квадратичным трением.
Проблема решается в рамках классической динамики твердого тела и рассматривает основные типические случаи вращения динамически неуравновешенного астатического ротатора: свободное вращение под действием резистивных факторов; вращение под действием постоянного тягового момента; свободные колебания под действием упругой пары и вращательные колебания под действием вынуждающей гармонической и ангармонической («ступенька», «пила») пар.
1. Тормозящий момент
Решение классической проблемы нахождения динамических реакций в осях вращающегося тела дает следующие выражения для возникающих усилий в подвижной (Б) и неподвижной (С) опорах (см., например, [5])
Д еJxz Ю -Зух ( ^ еJyz + Ю JZХ
ВУ = ^ 1 °Х =
y h ' ) — h
+ ю 2 Jzx | O Jx~ h ■ [иУ h
Bx = ^yz ! ^2jzx. | 0y = _eJxz ю Jyz , h = OB. ()
Эти несимметричные по отношению к обращению времени величины, тем не менее, складываясь своими проекциями по теореме Пифагора, дают пару поперечных сил, не зависящую от направления вращения
Jjxz + J2z) (£2 + ю4)
No = Nb = ^- У
h
и вызывающую силы трения с их тормозящим осевым моментом
26л¡ГЛ.
y/täz + J2yz) (е2 + ю4)
Mfr =--1---= _aV е2 + ю4. (2)
h
Размерность коэффициента а здесь совпадает с размерностью инвариантного при
повороте эффективного центробежного момента инерции у + JУjz.
Не имеет принципиального значения, какой вид сухого трения имеет место в опорах - качения или скольжения. И в том и в другом случае возникающий момент М^ пропорционален итоговой силе N0 + Nв нормального давления на опоры.
Отличительной чертой такого трения (назовем его инерционным) является необращение его величины в ноль в точке разворота вращения. В этот момент из-за смены направления Mfr следует ожидать скачка ускорения е и, следовательно, удара.
Что же касается сил аэродинамического сопротивления вращению, то они традиционно описываются квадратичным моментом Мв = —с |ю| ю, от ускорения не зависят и во времени изменяются плавным образом. Коэффициент пропорциональности с здесь также имеет смысл некоторого момента инерции.
2. Вращение по инерции
Пусть тело имеет угловую скорость юо, когда перестает действовать ведущий
момент, и далее оно движется лишь под действием тормозящих пар Мв и М^. Очевидно, что направление вращения не изменится, и динамика его будет описываться
уравнением
ф = -а\/ф4 + ф2 - сер2. (3)
Его разрешение относительно старшей производной возможно через перегруппировку и возведение в квадрат
т2 •• 2 , о т •• • 2 , 2-4 2 •• 2 , 2-4
З^ф + 23 ф ер + с ер = а ф + а ф ,
и оно неоднозначно
.. -сЗ:: ± аТГГс2 —О2 . 2
ф = -т 2 _ о-ф •
т:: а
Знак перед радикалом выбирается из очевидного условия
т .. . 2 Л + с2 — а2 — са
ф + сф2 = а - < 0. (4)
3::- а
Тогда при Jzz > а должен быть выбран «—»:
• • -2 cJzz + ал/ Jz2z + С2 — а2 -2 (5)
Ф = —Ф -J2 _а2-= —YCP • (5)
J zz а
Уравнение можно представить как
7 .. -2 -2 са2 + аJzz\| Jz2z + c2 — а2 Jzz Ф = — Сф2 — Ф2- 2_ 2-, (6)
J zz а
и момент трения, согласно А. И. Лурье [3], определится квадратичным выражением
M ( m) са2 + ^zz^Jz2z + с2 — а2 „2 (7)
Mfr (m) =--а2-m • (7)
J zz а
В самом общем случае параметров уравнение легко интегрируется и дает верное решение
Ф 11 / ч m0 / ч ln(1 + movi)
— ^ = Y ^---= Yt ^ m(t) =-0—, ф(í) = —-— • (8)
ф2 r m m0 r w 1 + Y
Здесь, в отличие от классического линейного по m сопротивления и от тормозящего постоянного момента, угловой путь до полной остановки логарифмически расходится со временем из-за медленно убывающей угловой скорости.
Таким образом, кельтский камень, закрепленный на проходящей через центр масс, то есть на центральной оси с трением в опорах, вращается по инерции точно так же, как и в отсутствие такового трения, но при наличии квадратичного момента сопротивления. Но множитель такого сопротивления помимо аэродинамических величин определяется также коэффициентом трения ô, центробежными моментами инерции и расстоянием между опорами
, = cJz2z + а Jzz J + с2 — а2 = + J2yz
С Jz2z — а2 ' а h (9)
и не обращается в нуль даже в безвоздушном пространстве (с = 0).
Последнее означает, что инерционное трение при значительном своем коэффициенте а поддерживает квадратичность диссипативной пары для астатического ротатора вплоть до очень малых оборотов, несмотря на то, что в аэродинамике на малых скоростях ньютоновский квадратичный закон переходит в линейный стоксовый.
Стоит отметить, что в частных (вырожденных) ситуациях, например Jzz = а = = с, уравнение (3) решения не имеет при начальном условии ю(0) = юо = 0. Также бессмысленно решение (8) с уходом угловой скорости в бесконечность при у < 0, то есть когда осевой момент инерции ротатора Jzz < а.
По мнению Пэнлеве [4], подобные ситуации свидетельствует о неэффективности динамики а.т.т. и законов Кулона-Амонтона в решении механических проблем, связанных с сухим трением. Эти и им подобные частные случаи требуют особого подхода (в печати), учитывающего связанное с конечностью скорости звука в материале ротатора запаздывание углового ускорения е(Ь — т) (и скорости ю(Ь — т)), определяющего динамические реакции, от ускорения е(Ь), задающего темп изменения осевого кинетического момента ротатора Kz. Другими словами, подхода, так или иначе уводящего проблему в область динамики реального деформируемого тела.
3. Вращение под действием постоянного разгоняющего момента
Пусть постоянный во времени осевой (Oz) момент М приложен к неподвижному телу в момент времени Ь = 0. Соответствующее уравнение разгона тела вокруг этой центральной оси будет отличаться от (5) постоянным слагаемым в правой части
Jzz ф = М — а\/ф4 + ф2 — сф2. ( 1 0)
Предельная скорость разгона определится из условия ф ^ 0
М — сф2 — аУф4 = 0 ^ юоо = л/(11)
V с + а
Аналогичным образом оно неоднозначно разрешится относительно старшей производной
ф2 ^ — а2) — 2фJzz (М — сф2) + ф4 (с2 — а2) + М2 — 2сф2М = 0 ^
.. -Угг (М—сф2) J2z (М—сф2)2 — ^—а2) (ф4 (с2—а2) — 2сф2М + М2)
Jzz — а2
Jzz (М — сф2) — аЛЩТ^+с2—а2)—2«ф2МТМ2
ф =
J2z а2
(12)
Отрицательный знак при Jzz > а выбирается аналогичным (4) образом. Анализ данного уравнения по методу А.И. Лурье
ф = (М — сср2) — (М — сср2) +
Jz2z (М — сф2) — аJzzУф^Л^+с2—О2)—2сф2Ш+М2
+ J2 а2 ^^
2 а (М — сф2) — ф4 ^ + с2 — а2) — 2сф2М + М2 Jzzф = М — сф2 + а—---J2 _ а2-
дает выражение для момента трения
а (М - сю2) - 322\/(М - сю2)2 + 32ю4 Мг (ю, М) = а^---"-— (13)
и его начального и предельного значений -
Мь (0, М) = Мь (юте,М) = - — ■ (14)
3 хх + а с + а
Само уравнение разгона (12) автономно, и оно решается стандартным методом перехода к переменной ю = ф (ф)
йю 3 гг (М - сю2) - а^ю4 ' + с2 - а2) - 2сю2М + а2М2 йф - а2
Последнее эквивалентно уравнению
йф - а2) ю
йю 3гх (М - сю2) - а^/ю4 ('2х + с2 - а2) - 2сю2М + а2М2 : решение которого в параметрической форме -
)
(15)
Г '^х- а2) ю'йю'
ф (ю) = -, =,
У 322 (М - сю'2) - а^ю'4 (+ с2 - а2) - 2сю'2М + а2М2
(16)
о
ю
- а2) йю'
г (ю) = -. 4 22 7 ■ (17)
0 .322(М - сю'2) - а^ю'4 ' + с2 - а2) - 2сю'2М + а2М2
Выражение (16) для угла поворота в принципе может быть проинтегрировано после цепочки стандартных замен как рациональная дробь со знаменателем 4-го порядка. Что касается интеграла (17), то он может быть сведен к эллиптическому.
Численное интегрирование дает в общем-то стандартную кинетику разгона ротатора (рис. 2). Однако при одинаковых тяговых моментах и одинаковых предельных скоростях темп разгона при наличии только квадратичного аэродинамического сопротивления представляется более быстрым, нежели когда таковое сопротивление дополняется еще инерционным трением.
Без учета квадратичной добавки интегрирование (12) дает фазовый портрет системы на плоскости (ф, ю)
. . V32х - а2 ,
ф (ю) = у - 1п
а
(
^ю4 (322 - а2) + М2 - М - ю2д/31- а2
^ю4 (322 - ах2) + М2 - М + ю2д/31- а2 4 (ю4 (32 - а2) + М2 - М)
(18)
-4 1п
'¿г + а
(
^ю4 ('2 - а2) + М2 - М - ю2 (3хх - а)
^ю4 (322 - а2) + М2 - М + ю2 (322 - а)
Примечательно, что требование астатичности ротатора в рассмотренной задаче не является обязательным. Даже при наличии статических составляющих в реакциях опор и вызванных ими дополнительных сил трения решение (16), (17) не изменится по существу. Просто это дополнительное постоянное по величине и направлению трение может быть включено со знаком минус в разгоняющий момент M.
Вообще же реально при высоких оборотах вращения и резких разгонах и торможениях наличие статических реакций в опорах и связанного с ними трения не является проблемой, и может быть учтено как малое возмущение.
Рис. 2. Зависимость угловой скорости от времени при постоянном тяговом моменте при наличии (1) и отсутствии (2) сухого трения
4. Вращательные колебания (гармонический момент)
При наличии возвращающего упругого момента Mei (ф) = —Кф такой ротатор способен совершать колебания как свободные, так и вынужденные под действием гармонически изменяющейся со временем пары M¡nd (t) = Ao sin (Qt). Первый тип без учета квадратичной добавки описывается уравнением
Jzzф = — sign (ф) ад/ф4 + ф2 — Кф. (19)
Разрешение его относительно второй производной дает
.. = —JzКф — sign (ф) а^/ф4 (J2 — а2) + К2ф2 ф= Jz2z — а2 . ( )
Для вынужденных колебаний, соответственно —Jz (Кф — Ao sin Qt)
ф =-J^o2--
Jzz u_
sign (ф) a^/—2кА0ф sin Qt + A0 sin2 Qt + ф4 (J2 — a2) + К2ф2 ( )
Jfz — a2
Простым преобразованием правой части можно выделить член, ответственный за трение в осях и величину поперечных усилий в опорах O и B:
.. — Кф + Ao sin Qt + Кф — Ao sin Qt + — Jz (кф — Ao sin Qt)
Jz _Jz_Jz2z a_
sign (ф) a^j—2кА0ф sin Qt + A^ sin2 Qt + ф4 (J2 — a2) + к2ф2
& Jz ф = — Кф + Ao sin Qt +
2 Jzz a2
a2 (Ao sin Qt — Кф)
a sign (ф) Jz \J—2кА0ф sin Qt + A^ sin2 Qt + ф4 (J| — a2) + к2ф2
Mfr (ф, ю, t) =
а2 (А0 sin Qt — Кф)
а sign (ю) Jzy^к2ф2 — 2кА0ф sin Qt + А2 sin2 Qt + ю4 (J| — а2)
(22)
= = f 2о
Достойно внимания то, что своей плавной, не зависящей от угловой скорости частью, момент Mfr может значительно усиливать консервативную возвращающую и вынуждающую гармоническую пары при примерном равенстве характеристик Jzz и а.
Что же касается разрывной части с множителем sign (ю), то она ответственна за удар, амплитуда которого
Im
2aJz |кфт — А0 sin Qt|
Jz2z - а2
помимо соотношения между Jzz и а зависит также от углового положения точки разворота фт (в установившемся режиме - просто угловой амплитуды) и от сдвига фаз между ф^) и А0 sin Qt. Как поперечная сила этот удар будет передаваться обеим опорам O и B с множителем 1 /(26) (22).
Уравнения (20) и (21) в принципе можно приближенно проинтегрировать одним из принятых в теории нелинейных колебаний способов (см. например [6]). В частности, при помощи теории возмущений, рассматривая в качестве такового всё или только разрывную часть Mfr.
Численный расчет (Maple 15) дает при высоких значениях коэффициента трения зависимости для угла и скорости, существенно отличающиеся по форме от свободных линейных колебаний (рис. 3). Так, для кривой w(t) характерны изломы в точках разворота, тем ярче выраженные, чем больше параметр а. Примечательно также, что затухание идет существенно быстрее, чем в линейных системах с аналогичными параметрами (штрихпунктир).
У вынужденных колебаний ангармоничность также сильнее выражена при высоких показателях трения (рис. 4).
Рис. 3. Кинетика свободных колебаний
Рис. 4. Вынужденные колебания
Рис. 5. Фазовый портрет системы при малых (а) и больших (б) параметрах трения а
Интересен фазовый портрет таких колебаний с характерной неправильностью предельного цикла (рис. 5), однако сохраняющего центральную симметричность. Примечательно, что максимальная угловая скорость на этом цикле достигается вдали от положения равновесия, скорее ближе к поворотным точкам.
При малых параметрах а и вид колебаний, и предельный цикл близки по форме к соответствующим кривым линейных вынужденных колебаний.
5. Колебания под действием ангармонического момента
5.1. «Ступенька». В ряде чисто технических проблем представляет интерес движение ротатора при чередующихся по направлению постоянных внешних моментах, то есть когда периодическое внешнее воздействие имеет в идеале ступенчатую форму. Последнюю можно описать сингулярной формулой
Mind (t) = sign (sin Qt) . (23)
Каким образом наличие инерционного трения влияет на форму вынужденных колебаний, на портрет системы на фазовой плоскости? Чтоб ответить на этот вопрос, необходимо исследовать динамическое уравнение, имеющее вид
.. — Jz (кф — A0 sign (sin Qt))
J2 - a2
° zz __(24)
sign (ф) ад/—2кА0ф sign (sin Qt) + AO + ф4 (J| — a2) + к2ф2 v '
J2z — a2
Его решение можно получить аналитически из (16) путем сшивки в особых точках, а также и численно для конкретных значений параметров.
Помимо двух, по закону положенных изломов в точках изменения полярности момента на квазипериоде кривой w(t) (рис. 6), есть еще два, соответствующие точке разворота. Изломы наблюдаются и на предельном цикле (рис. 7) в фазовой плоскости.
Рис. 6. Кинетика ф(£) и ю(£) при больших тре- Рис. 7. Фазовый портрет астатичного ротатора
ниях при действии ступенчатого момента
при больших а
J = l, 01=0.8, К=1 COA
Рис. 8. Фазовый портрет стандартной линейной системы с большим трением («ступенька»)
J= 1, а=0.4. к=8, £2=0.7, Q2=197989898
¡ Н о -1 \у И А á щ 'ц ¡¡i I 0.05М У ^ Ч i Ж / ü Щ 'i b / A /i A /¡A Ф''|П /!|5 /И /
Ж i у ч U' 1 ■ HiM щ ! Д! i ¡ / V i !/ i' г i i I/ V ¡/ У / V I' t /
Рис. 9. Кинетика ф и ю при подаче пилообразного вынуждающего момента
Однако аналогичные изломы имеют место и на предельном цикле линейной колебательной системы, диссипативный момент в которой пропорционален угловой скорости (рис. 8).
Что касается вынуждающей частоты, то для обеих систем она выбиралась равной резонансной для линейной системы Q2 = л/k/Jzz (1 — а2/J2z). Уже упомянутая сдвоенность изломов на кривой угловой скорости свидетельствует о вдвойне опасных ударных воздействиях на ось в случае наличия инерционного трения.
5.2. «Пила». Периодический пилообразный момент моделировался формулой
M(t) = M0 cos (arceos (Qt)) sign (sin (Qt)).
(25)
Частоты Q выбирались как близкие, так и далекие от «резонанса» Q2, начальные условия - по типу «отклонение» и «толчок» (импульс).
При различных значениях вынуждающей частоты Q установившаяся кинетика угла поворота и угловой скорости (рис. 9, 10) ротатора практически совпа-
а б
Рис. 10. Кинетика угла поворота (а) и угловой скорости (б) у астатичного ротатора и классического линейного осциллятора при Jz = 2.0, а = 0.9, к = 4, Я = 0.4Я2
дала с аналогичной кинетикой классической линейной системы с тем же моментом инерции , но диссипативным моментом Mfr(í) = —а • 1 • ю (множитель 1 с-1 - для нужной размерности), отличаясь только дополнительными изломами ю(Ь) в точках разворота.
Таким образом, можно заключить, что в исследованной многомерной области параметров как при гармоническом вынуждающем воздействии, так и при различных ангармонических исследуемый ротатор, в основном, ведет себя как линейная осциллирующая система, отличаясь лишь в деталях второго порядка.
Разумеется, проблема схожести предельного поведения астатичного колебательного ротатора и линейного осциллятора с инерционным трением требует более глубокой проработки и выходит за рамки настоящей статьи. Однако полученные предварительные результаты оправдывают повсеместное применение простой и удобной линейной модели для демонстрации принципиальных свойств колебательного устройства.
Выводы
Изучено как свободное, так и вынужденное движение важной динамической системы, весьма близкой по физическим свойствам к реальной: астатического ротатора с неглавной осью вращения. Примечательным свойством такой системы является возникновение инерционного трения в оси, а вместе с ним и двусторонней причинно-следственной связи между угловым ускорением и тормозящим моментом такого трения.
Показано: а) свободное вращение по инерции эквивалентно движению с квадратичным по угловой скорости сопротивлением; б) вынужденные колебания под действием как гармонического, так и различных ангармонических моментов по основным характеристикам схожи с аналогичными колебаниями простых линейных систем с диссипацией.
Различие с последними выражается в наличии дополнительных ударных воздействий на ось в поворотных точках, где скачком изменяется момент сухого трения М^.
Изучение таких ударов представляет практический интерес в плане анализа и моделирования хода и последствий аварийных ситуаций, связанных с развитием у ротатора (вращательного осциллятора) динамической неуравновешенности.
Кроме того, несмотря на схожесть предельного поведения с классической линейной системой, рассмотренный ротатор с часовой пружиной и трением может быть включен в библиотеку нелинейных колебательных моделей для существующих компьютерных имитационных программ типа Nonlinear Oscillations (Physics Academic Soft) и др.
Работа выполнена в рамках государсвенного задания Минобрнауки РФ по теме «Разработка компетентностно-ориентированной модели обучения физике в общеобразовательной средней школе», № государственной регистрации 01201261999.
Библиографический список
1. Смирнов Ю.П. Об уравнениях динамики систем с трением// Сборник научно-методических статей по теоретической механике. Вып. 11. М.: Высшая школа, 1981. С. 184.
2. Borisov A.V., Jalnin A.Yu., Kuznetsov S.P., Sataev I.R., Sedova J.V. Dynamical phenomena occurring due to phase volume compression in nonholonomic model of the rattleback // Regular and Chaotic Dynamics. 2012. Vol. 17, № 6. P. 512.
3. Лурье А.И. Аналитическая механика. М.: Наука, 1961.
4. Painleve P. Leçons sur le frottement. Paris: Hermann, 1895. (Пэнлеве П. Лекции о трении. М.: Гостехиздат, 1954.)
5. Маркеев А.П. Теоретическая механика: Учебник для университетов. М.: ЧеРо, 1999.
6. Бать М.И., Джанелидзе Г.Ю., Кельзон А.С. Теоретическая механика в примерах и задачах. М.: Наука, 1973.
Поступила в редакцию 20.01.2014 После доработки 3.02.2014
ABOUT ROTATIONAL DYNAMICS OF A RIGID BODY AROUND NON-PRINCIPAL AXIS PASSING THROUGH THE CENTER OF MASS UNDER DRY FRICTION ACTING
V. V. Chistyakov
Yaroslavl State Pedagogical University named after K.D. Ushinsky
The dynamics is studying for rigid body rotating around fixed axis Oz being central but not principal. Therefore the inertial torques Mx and My arose depending both on mass geometry Jxz, Jyz and on angular velocity rn and acceleration e. Dry friction acting on axis's supports with coefficient 6 leads to that the value of e serves as the reason and result of the motion simultaneously. There were integrated numerically and/or analytically the dynamical equations of free and forced motion including rotational harmonic and
inharmonic oscillations too. The results obtained are comparing with those following from the standard linear equations.
Keywords: Astatic rotator, axis of inertia, inertial torque, dry friction, quadratic drag, rotational oscillations.
Чистяков Виктор Владимирович - родился в городе Буй Костромской области (1957), окончил факультет общей и прикладной физики Московского физико-технического института (1981). После окончания МФТИ работал инженером-программистом в НИИ Вычислительных Комплексов (Москва). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук (1992, НИФХИ им. Л.Я. Карпова, Москва). Работал в учреждениях науки и высшего образования Ярославля. Автор более 60 научных и научно-методических трудов и изобретений. Последние 4 года работал в области динамики резистивного движения, опубликовал цикл работ по интегрированию динамических уравнений в условиях, приближенных к реальным. В настоящее время доцент Ярославского ГПУ им. К.Д. Ушинского, Ярославской ГСХА и ЯФ ЛГУ им А.С. Пушкина.
150000 Россия, Ярославль, ул. Республиканская, 108 ФГБОУ ВПО Ярославский ГПУ им. К.Д. Ушинского E-mail: v.chistyakov@yarcx.ru, chistiakov_v_v@rambler.ru