Араметрическая идентификация тензоров инерции тел на сферических движениях с медленным собственным вращением Текст научной статьи по специальности «Физика»

Варжель Сергей Владимирович - Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант, vsv187@gmail.com

Куликов Андрей Владимирович - Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант, a.kulikov86@gmail.com

Брунов Вячеслав Сергеевич - Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, студент, brunov@oi.ifmo.ru

Асеев Владимир Анатольевич - Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, научный сотрудник, aseev@oi.ifmo.ru

УДК 681.5 + 531

ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ТЕНЗОРОВ ИНЕРЦИИ ТЕЛ НА СФЕРИЧЕСКИХ ДВИЖЕНИЯХ С МЕДЛЕННЫМ СОБСТВЕННЫМ ВРАЩЕНИЕМ С.Н. Шаховал, Г.И. Мельников

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

Элементы тензора инерции определяются либо на вращениях вокруг шести неподвижных осей, либо на прецессии, дополненной одним осевым вращением. При этом применяется реверсивно-симметричный характер движения [1, 2]. В одном из вариантов шесть осевых вращений связаны переходными процессами угловой выставки тела относительно неподвижной оси, что характеризуется снижением точности и дополнительными затратами времени. В связи с этим в данной работе предлагается вариант, в котором осуществляется двухосная полупрограммная реверсивно-симметричная прецессия тестируемого тела, при условии, что угловая скорость собственного вращения тела на порядок меньше угловой скорости прецессионного движения вокруг вертикальной оси (рисунок). При таком условии весьма медленное собственное вращение обеспечивает на протяжении полного собственного оборота совпадение с вертикальной осью прецессии всех линий вокруг наклонной оси собственного вращения кругового конуса, связанного условно с телом (подвижного аксоида), что заменяет пять переходных процессов смены осей. Собственная угловая скорость тела вносит лишь малые изменения в величину и направление (угол у на рисунке) вектора угловой скорости сферического движения тела. Вращения вокруг двух осей предполагаются синхронными реверсивно-симметричными, содержащими симметричное движение в обратном направлении.

Пусть ю1 = \\ к1 - переменная угловая скорость прецессии тела вокруг вертикальной оси 0хх, меняющая направление при обратном движении; е1 - угловое ускорение прецессии; ф - угол поворота тела вокруг наклонной оси собственного вращения тела; ю2 = ф к - угловая скорость вращения вокруг наклонной оси. Предполагаем ю1 >> ю2. Тогда в первом приближении угловую скорость ю2 вращения вокруг наклонной оси можно не учитывать. Нас интересует только изменение с течением времени углового положения тела относительно оси 02х. Пусть в процессе испытания тело совершило два оборота замедленного вращения вокруг оси 02 в положительном направлении и два оборота симметричного разгона в обратном направлении. Считаем, что в виду симметричности пары движений работы сил трения в обоих направлениях равны. Работы силы тяжести неуравновешенного тела на полном обороте по углу ф

равны нулю в виду вертикальности оси прецессии 021. Программные движения [3] осуществляются двумя электродвигателями (рисунок). Назначим на угловом интервале 0 < ф< 4л пять полных оборотов,

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2012, № 1 (77)

2л

отсчитываемых от следующих пяти угловых значений: ф0 = 0, ф1 = h , ..., ф4 = 4h ; h = ^ . На этих

оборотах имеем пять реверсивно-симметричных вращений (в одну и обратную стороны). Эти движения удовлетворяют парам уравнений энергии (в которых пренебрегаем угловой скоростью собственного вращения):

(J + Ik)(a>u+i -и?,k) = 2Ak + 2Fk, (J + Ik)(Ю12,к -и?,k+1) = 2A'k + 2Vk . (1)

В уравнениях (1) и1 k (k = 0,...,4 ) - узловые значения угловой скорости прецессии; J - приведенный момент инерции устройства; Ik - моменты инерции тела относительно оси Oz1 в угловой позиции k (относительно мгновенной оси, проведенной в теле при узловом значении ф = <k); Ak , A'k - полезные работы электродвигателей на движениях в одном и обратном направлениях; Vk и Vk = Vk - отрицательные работы диссипативных сил трения и аэродинамического сопротивления. Из (1) получаем расчетную формулу для пяти осевых моментов инерции тела:

Ik = ?Ak- Ak--J , k = 0, ..., 4. (2)

и1, k -и1, k+1

Эти оси расположены в теле на круговом конусе, описанном вокруг собственной оси Oz с углом а и 63,4° . Формулы (2) являются приближенными, так как в них не учитывается угловая скорость собственного вращения ю2. Точная формула, очевидно, имеет следующий вид:

Т At - Ak Т

Ik = 2k 2k - J , И -rok+1

где rak = w1 k + ю2 k ; =ro2k +ю2 k + 2и1 k и2 k cos а k ; угол ak - угол между ю1 k и ю2 k . При этом осевые моменты инерции тела Ik определяются относительно пяти мгновенных осей, направленных по векторам , на подвижном аксоиде. Дополнительно определяют в эксперименте момент инерции относительно собственной оси вращения Oz. Затем по шести моментам инерции находят матрицу тензора инерции в заданной точке O тела. Малые погрешности формулы (2) оцениваются с использованием способа [4]. Работа выполнена при финансовой поддержке гранта РФФИ проект № 10-08-01046.

1. Мельников В .Г. Энергетический метод параметрической идентификации тензоров инерции тел // Научно-технический вестник СПбГУ ИТМО. - 2010. - № 1 (65). - С. 59-63.

2. Мельников В.Г., Едачев А.С., Мельников Г.И., Шаховал С.Н. Метод определения тензора инерции на программных движениях // Изв. Самарского научного центра РАН. - 2010. - Т. 12 (33). - № 1 (2). -С. 445-448.

3. Мельников В .Г. Использование программных движений для идентификации тензора инерции и центра масс твердого тела // Изв. вузов. Приборостроение. - 2007. - Т. 50. - № 8. - С. 33-36.

4. Шаховал С.Н. Исследование матричных алгебраических уравнений, определяющих тензор инерции через осевые моменты инерции // Научно-технический вестник СПбГУ ИТМО. - 2008. - Вып. 47. -С. 196-201.

Шаховал Сергей Николаевич - Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант, shakhovalsergey@gmail.com

Мельников Геннадий Иванович- Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, доктор физ.-мат. наук, профессор, melnikov@mail.ifmo.ru

УДК 531

ИДЕНТИФИКАЦИЯ ТЕНЗОРА ИНЕРЦИИ ТЕЛА НА РЕВЕРСИВНО-СИММЕТРИЧНЫХ ПРЕЦЕССИЯХ В ОГРАНИЧЕННОМ УГЛОВОМ ИНТЕРВАЛЕ В.Г. Мельников, Р.Ю. Кравчук, Г.И. Мельников, С.Н. Шаховал

Решается задача определения тензора инерции на основании прецессионного двухосного реверсивно-симметричного сферического движения тела при ограничениях на два угловых интервала. Получены расчетные формулы для инерционных параметров.

Ключевые слова: матрица тензора инерции, угол прецессии, угол собственного вращения, реверсивно-симметричная прецессия.

Известен способ параметрической идентификации пяти осевых моментов инерции тела на реверсив-но-симметричном прецессионном двухосном движении тела вокруг неподвижной точки. Обобщим этот способ на случай отсутствия полных оборотов по углу собственного вращения. Пусть твердое тело вместе со сцепленной с ним системой Oxyz, собственной осью Oz c известным центром масс С, собственным углом поворота ф совершает сферическое движение вокруг неподвижной точки O. При этом за рассматриваемый конечный интервал времени мгновенная ось вращения тела описывает вокруг собственной оси Oz тела

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2012, № 1 (77)