Научный журнал
Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955
ИФ РИНЦ = 0,580

ОЦЕНКА ВЛИЯНИЯ ПОЯСКА НА ДЕФОРМАЦИЮ РОТОРА ГИРОСКОПА В СЛУЧАЕ ЕГО ОСЕСИММЕТРИЧНОГО ВРАЩЕНИЯ

Карипбаев С.Ж. 1 Аринов Е. 2
1 АО «Академия гражданской авиации»
2 АО «Жезказганский университет им. О.А. Байконурова»
В данной работе изучается движение вокруг центра масс неоднородного шарового ротора электростатического гироскопа, взвешенного внутри вакуумированной полости. Неравенство моментов инерции обеспечивается за счет узкого кольцевого пояска в экваториальной области ротора, изготовленного из материала с более высокой плотностью. Для количественной оценки влияния пояска на деформацию ротора в случае его осесимметричного вращения, определяется уравнение радиальной деформации поверхности ротора без учета и с учетом кольцевого пояска. Определяется отношение амплитудных значений найденных деформации для реального ротора гироскопа.
шаровой ротор
электростатический гироскоп
момент инерции
неконтактный подвес
деформация
вызванная центробежными силами
1. Anfinoqenov А.S., Gusinsky V.Z., Parfenov 0.I. – Electrostatic gyro. // The second soviet – Chinese symposium on Inertial technology. 9-15 October 1992. S.-Petersburg, 1992. P. 71–80.
2. Vodicheva L.V., Lookin N.A. Strap down inertial navigation – system and special processors design problems, SINS development – experience // The IV Russian-Chinese symposium on inertial technology. September 27 – October 1, 1993. S. – Petersburg, 1993. – P. 67–76.
3. Duncan R.R. A strap down inertial navigation using miniature electrostatic gyro // Proceeding of the National Aerospace Meeting. Washington, – 13-14 March. – 1973. – P. 13.
4. Чарышев Ш.Ф., Брюшков В.Г. Базовые чувствительные элементы ИНС. Электростатический гироскоп. – М., 1988. – 148 с.
5. Мартыненко Ю.Г. Движение твердого тела в электрических и магнитных полях. – М.: Наука, 1988.
6. Булгаков Б.В. Прикладная теория гироскопов. – М.: Изд-во Моск. ун-та, 1976. – 401 с.
7. Карипбаев С.Ж., Ландау Б.Е., Мартыненко Ю.Г., Подалков В.В. Зависимость угловой скорости электростатического гироскопа от температуры окружающей среды // Изв. РАН. МТТ. -1993. – № 3. – С. 42–49.
8. Корецкий А.В. О коэффициентах несферичности формы ротора электростатического гироскопа // Тр. Моск. энерг. ин-та. 1982. – Вып. 573. – С. 27–31.
9. Александров А.М., Брюшков В.Г., Корецкий А.В., Мартыненко Ю.Г. Уход электростатического гироскопа, вызываемый упругими деформациями ротора // Межвузовск. сб. трудов. – М.: Моск. энерг. ин-т. – 1983. – № 14. –
С. 16–22.
10. «Разработка бескардановых гироскопов с шаровым ротором на электростатическом и шарикоподшипниковом подвесах» за 2012-2014 гг. Отчет о научно – исследовательской работе ГРНТИ 30.15.35, № госрегистрации: 0112РК02743, Инв: № 0212РК01519, Инв: № 0213РК01969

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

Основным достоинством неконтактного подвеса ротора является практическое полное отсутствие сил трения при его вращении. Это открывает принципиальную возможность повышения точности гироскопических приборов. Существенным преимуществом ЭГС является возможность его использования при неограниченных углах поворота летательного аппарата вокруг центра тяжести без каких либо дополнительных устройств типа карданова подвеса. В этом случае корпус гироскопа устанавливается на движущемся объекте, совершающем произвольное движение [1–3].

Активные исследования навигационных датчиков с электростатическими подвесами ведутся в США (Honeywell, Stanford University), Франции (Sagem), Китае (Tsinghua Jiaotong Universities) и в Росcии (ЦНИИ Электроприбор, МИЭА, РПКБ, Институт прикладной механики) [4].

В работах [5] проводится анализ проблем создания гироскопа с электрическим подвесом. Даются рекомендации по выбору материала. Проводятся данные о величинах деформации ротора при действии центробежных сил.

ЭСГ имеет ряд преимуществ по сравнению с другими датчиками ИНС: высокая точность (до 10-10I/с), длительная безотказная работа на выбеге ротора (до нескольких лет), малое энергопотребление1 (до нескольких ватт), небольшие габариты и масса. ЭСГ мало подвержен износу, вследствие чего надежность прибора в основном определяется надежностью и сроком службы электронных элементов. Опыт эксплуатации ЭСГ на морских объектах подтвердил высокую точность и достаточную надежность корабельных ИНС на ЭСГ [6].

Использование ЭСГ на космических аппаратах привлекательно по той причине, что в условиях космоса легче поддерживать необходимую степень вакуума в гироскопе, решать задачу поддержании ротора во взвешенном состоянии, снизить энергопотребление системы. Вместе с тем, использование ЭСГ в условиях космического пространства требует решении комплекса новых задач, как в области теории, так и в области высоких технологий.

Технологические проблемы обработки поверхности ротора с точностью до 0.1 мкм оказывается весьма сложными и требуют создания специального оборудование[5]. В силу того, что первая гармоника формы ротора описывает его дисбаланс, то далее термин не сферичность ротора будет рассматриваться в обобщенном смысле и включать в себя и понятие несбалансированности ротора.

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

Поддержание вакуума внутри камеры до величины 1,3*10-6Па осуществляется ионно-геттерным насосом. Разгон ротора и демпфирование его нутационных колебаний, возникающих во время первоначальной раскрутки ротора, производится специальными катушками, которые после достижения ротором номинальной скорости отличаются на все время функционирования прибора.

Рабочая (номинальная) скорость вращения роторов в созданных образцах гироскопов составляет 12, 30, 60, 150, 180 тыс. об/мин. [45, 47]. Средняя величина зазора между ротором и электродами выбирается в пределах (3-40)*10-6м. Не сферичность внутренней поверхности полости не превышает (1-5)*10-7м. [46]. Диаметр ротора колеблется в диапазоне (30-70)*10-9м. для обеспечения вращения вокруг определенной оси и улучшения характеристик начальной выставки, ротор имеет в экваториальной области узкое утолщение, изготовленное из материала с боле высокой плотностью. Материал ротора выбирается с учетом способности сохранять форму, иметь достаточную упругость и малый удельный вес. Таким требованием лучше других материалов удовлетворяет бериллий (удельный вес 1,85*103 кг/м3) и алюминий (удельный вес 2,7*103 кг/м3).

При наличии не сферичности поверхности ротора возможно появление уходов из-за не сферичности электродов, смещений центра масс ротора в подвесе, возникающих при перегрузках и вибрациях основания и при отсутствии нулевого электрода, заполняющим междуэлектродное пространство подвеса и т.д. [7, 8].

Деформации ротора ЭСГ, вызванные центробежными силами после его раскрутки до номинальной скорости, могут приводить при перегрузках в 1 g к весьма большим возмущающим моментам и уходом ЭСГ до десятых долей градуса в час и более. В работе [5] рассматривается ротор ЭСГ, предназначенного для БИНС. Ротор выполнен в виде сферической оболочки с кольцевым утолщением в экваториальной плоскости. Уходы такого ротора в ЭСГ с шести электродным подвесом при произвольных положения оси вращения относительно электродов подвеса могут доходить до 2 градусов в час, если не проводить дополнительной обработки поверхности ротора, придавая ему специальную форму с тем, чтобы после раскрутки поверхность ротора становилась сферической (предварительная «асферизация ротора»). Заметим, что в шести электродном подвесе основной вклад в уход ЭСГ вносят четвертая и в меньшей степени восьмая «гармоники» разложения поверхности ротора в ряд по полиномам Лежандра. В задаче, исследованной в [9], уводящие моменты при учете только второй гармоники имеют величину 0.20 градусов в час, при учете второй и четвертой гармоник – 1,5 градусов в час, при учете всех слагаемых до восьмой гармоники включительно – 1,6 градусов в час, учет последующих членов ряда может привести к уходам до 2 градусов в час.

Основная часть. Рассмотрим ротор электростатического гироскопа, в экваториальной области которого находится поясок, имеющий плотность, отличную от плотности шара. Предложим, что магнитное поле внутри кожуха отсутствует , вакуум идеальный и сам кожух гироскопа неподвижен. Центральный эллипсоид инерции ротора есть эллипсоид вращения [5]. Неравенство моментов инерции обеспечивается за счет узкого кольцевого пояска в экваториальной области ротора, изготовленного из материала с более высокой плотностью.

Тогда моменты инерции ротора относительно главных осей определяются формулами

kar01.wmf

kar02.wmf

где ρ* – плотность материала пояска, h, ф – толщина и угол, определяющий размеры пояска.

Формулы (1) были получены в предположении, что угол ф мал. Для реальных конструкции ф составляет 10 °÷20 °. Из соотношения (1) получаем

kar03.wmf.

Здесь

kar04.wmf. (2)

Если взять бериллиевый ротор с титановым пояском (для титана плотность ρ* = 5.2 г/см3) и угол ф = 15 °, a ϑ = 0.1 то для толщины титанового пояска имеем h = 0.1 R.

Рассмотрим деформированное состояние ротора без учета пояска.

Используя закон сохранения массы имеем,

kar05.wmf (3)

где ρ’, v’ – плотность и объем ротора после деформации, ρ, v – до деформации; После деформации ротора получаем для элементарного объема dv’ соотношение

kar06.wmf

kar07.wmf. (4)

Здесь kar08.wmf – приращения линейных размеров элементарного объема dv после деформации вдоль осей x соответственно. Подставляя (4) в (3), после необходимых преобразований найдем

kar09.wmf

где

kar10.wmf.

Момент инерции ротора относительно оси х3 после его деформации определяется по формуле

kar11.wmf

kar12.wmf.

Данное выражение получено при помощи соотношения (5) и преобразовано с учетом малости компонент вектора перемещений.

Переходим к сферической системе координат и после некоторых преобразований получим выражение для момента инерции

kar13.wmf (6)

kar14.wmf

kar15.wmf (7)

kar16.wmf

kar17.wmf

При расчетах обычно используется разности моментов инерции, поэтому обозначая
через kar18.wmf kar19.wmf получим

kar20.wmf

kar21.wmf (8)

kar22.wmf

kar23.wmf

В результате подстановки

kar24.wmf

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

kar25.wmf (9)

Для бериллиевого ротора радиусом R = 0.5 см, при ω = 3000 об⁄с, ϑ = π⁄2, имеем kar26.wmf. Этот числовой результат был получен после осреднения последней формулы по времени. Разность моментов инерции kar27.wmf до деформации этого ротора в реальных конструкциях имеет следующий порядок kar28.wmf. Сравнивая ∆I2 и kar29.wmf, можно сделать вывод, что изменением момента инерции шара, обусловленного наличием деформации можно пренебречь. Вычисляя интеграл ∆I1 можно убедиться, что kar30.wmf.

Теперь вычислим моменты инерции пояска относительно осей координат xi

kar31.wmf

kar32.wmf

kar33.wmf (10)

kar34.wmf

kar35.wmf.

Здесь

kar36.wmf;

kar37.wmf;

kar38.wmf.

Введем обозначение kar39.wmf, и, воспользовавшись (10) для разности моментов инерции, имеем

kar40.wmf

kar41.wmf. (11)

Подставляя kar42.wmf, kar43.wmf и kar44.wmf в последнее выражение (11) и интегрируя по углам α и β, окончательно получим

kar45.wmf. (12)

Здесь kar46.wmf

kar47.wmf

kar48.wmf

kar49.wmfkar50.wmf;

kar51.wmf.

Учитывая малость угла kar52.wmf, в данных соотношения можно положить kar53.wmf,
kar54.wmf и ограничиться линейными слагаемыми по kar55.wmf kar56.wmf kar57.wmf.

Для числовой оценки влияния деформации на разность моментов инерции пояска вычислим kar58.wmf по последней формуле (12) kar59.wmf, при рассмотренном выше случае

kar60.wmf.

Порядок разности моментов инерции до деформации kar61.wmf. Сравнивая kar62.wmf и kar63.wmf, можно сделать вывод, что изменением момента инерции пояска, обусловленного наличием деформации можно пренебречь.

Влияние пояска на деформацию ротора оценим в случае его осесимметричного вращения. Пренебрегая изменением упругих характеристик в зоне прикрепления пояска, заменяем его влияние на ротор инерционной нагрузкой

kar64.wmf

kar65.wmf, (13)

kar66.wmf.

На остальной поверхности ротора напряжения равны нулю.

Разложим правые части (13) граничных условий в ряд по полиномам Лежандра

kar67.wmf

kar68.wmf (14)

Коэффициенты этих разложений определяются по формулам [56]

kar69.wmf,

kar70.wmf.

Решение уравнений Ляме в осесимметричном случае известно и имеет вид

kar71.wmf

kar72.wmf (15)

kar73.wmf

kar74.wmf

где постоянные интегрирования An и Bn находится при удовлетворении граничных условии (13).

kar75.wmf (16)

kar76.wmf

При n = 2 для коэффициентов разложения в ряд по полиномам Лежандра имеем

kar77.wmf (17)

kar78.wmf

После подстановки (17) в (16) найдем коэффициенты системы уравнений при n = 2:

kar79.wmf

kar80.wmf. (18)

Воспользовавшись (13) и (18), находим уравнение радиальной деформации поверхности ротора, вызванное наличием кольцевого пояска в случае его осесимметричного вращения

kar81.wmf. (19)

Для количественной оценки влияния пояска на деформацию ротора в случае его осимметричного вращения определяем, воспользовавшись формулой

kar82.wmf

kar83.wmf

уравнение радиальной деформации поверхности ротора без учета кольцевого пояска

kar84.wmf

Затем вычислим отношение амплитудных значений kar85.wmf к ur для ротора, физические и геометрические характеристики были описаны выше

kar86.wmf (20)

Выводы

Из (20) видно, что наличие кольцевого пояска дает погрешность в вычислениях не более 10 процентов, следовательно, kar87.wmf, обусловленный пояском, несуществен. Поэтому при определении деформации ротора наличием пояска можно пренебречь.


Библиографическая ссылка

Карипбаев С.Ж., Аринов Е. ОЦЕНКА ВЛИЯНИЯ ПОЯСКА НА ДЕФОРМАЦИЮ РОТОРА ГИРОСКОПА В СЛУЧАЕ ЕГО ОСЕСИММЕТРИЧНОГО ВРАЩЕНИЯ // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 8-3. – С. 340-346;
URL: https://www.applied-research.ru/ru/article/view?id=10030 (дата обращения: 06.05.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074