ДНГ. Исследование полученной модели 49
![]()
|
3.7. Вывод уравнений движения динамически настраиваемого гироскопа с учётом угловой податливости скоростной опоры.Для вывода уравнений математической модели ДНГ, учитывающей упругую угловую податливость скоростной опоры, воспользуемся II методом Лагранжа. Выбор данного метода обусловлен тем, что он несложно реализуется с помощью программных средств компьютера. Уравнение Лагранжа второго рода записывается следующим образом:
где T – кинетическая энергия системы, qi – обобщённая координата, Qi – обобщённая сила. При выводе уравнений примем ряд допущений:
В качестве обобщённых примем координаты, изменяющиеся в осях, не участвующих в быстром собственном вращении ротора, так как именно в этих осях происходит съём информации у управление движением маховика. Обобщёнными координатами будем считать углы наклона ротора относительно вала α1, β1 и углы наклона вала относительно корпуса (основания) α2, β2. Для вывода уравнений движения напишем программу в среде компьютерной алгебры Wolfram Mathematica. Алгоритм работы программы представлен на рис. 3.9. Исходный код представлен в приложении А. ![]() Рис. 3.9. Алгоритм работы программы вывода уравнений математической модели ДНГ Рассмотрим подробнее процесс получения уравнений движения динамически настраиваемого гироскопа с учётом угловой податливости скоростной опоры методом Лагранжа II рода. Определение кинетической энергии системы Из уравнения Лагранжа II рода (15) следует, что нам необходимо продифференцировать кинетическую энергию системы. Найдём её по формуле:
где Ji – момент инерции i-го тела (элемента механической системы); ωi – его абсолютная угловая скорость. Тензоры инерции тел, входящих в систему (4)–(7) описаны в механической модели ДНГ (см. п. 3.6). Необходимо определить абсолютную угловую скорость каждого элемента в осях, связанных с этим элементом. Основание движется относительно инерциального пространства с угловой скоростью ω0:
Абсолютную угловую скорость вала можно найти по формуле:
Абсолютная угловая скорость 1-го кольца:
Абсолютная угловая скорость 2-го кольца:
Абсолютная угловая скорость движения ротора:
Для получения уравнений движения ДНГ, записанных в невращающихся осях, необходимо в выражениях для угловых скоростей кардановых колец (19)–(20) заменить переменные αи β на обобщённые координаты α1 и β1. Для определения связи между α, β и α1, β1 найдём угловую скорость ротора относительно вала во вращающихся и в невращающихся осях, приравняем полученные выражения и из получившейся системы уравнений найдём выражения для ![]() ![]() Угловую скорость движения ротора относительно вала в осях x’y’z’ можно найти следующим образом:
Угловая скорость движения ротора относительно вала в осяхxyz:
При нахождении связи между координатами пренебрегаем величинами третьего и выше порядков малости, сохраняя второй [2]. Также ввиду большой величины скорости собственного вращения маховика считаем, что:
В итоге получаем следующие выражения, связывающие между собой углы α, β и α1, β1:
Подставив выражения для абсолютных угловых скоростей и тензоров инерции элементов системы в выражение (16) и произведя замену переменных согласно (25), получим кинетическую энергию системы. Определение обобщённых сил Правую часть уравнения Лагранжа II рода можно представить в следующем виде [2]:
где U – силовая функция системы; R – диссипативная функция Релея; x – координата. Силовая функция системы учитывает потенциальную энергию упругой деформации в перемычках подвеса и в скоростной опоре и имеет вид:
где Кα,β – жёсткость упругой перемычки подвеса вокруг оси α, β; Кα2,β2 – жёсткость скоростной опоры вокруг оси α2, β2. Диссипативная функция Релея учитывает вязкое трение в материале подвеса, воздушное сопротивление движению маховика в увлечённой массе газа и относительно корпуса:
где μподв – коэффициент вязкого сопротивления движению маховика в увлечённой массе газа; μнп – коэффициент вязкого сопротивления движению маховика относительно корпуса; μВ – коэффициент вязкого сопротивления движению вала относительно корпуса. Уравнения движения Подставив полученные выражения для кинетической энергии и обобщённой силы в уравнение Лагранжа II рода, а также проведя замену переменных согласно (25), получаем искомые уравнения движения. В полученных уравнениях пренебрегаем величинами второго и выше порядков малости, а так же составляющими, изменяющимися с частотами, кратными двойной частоте собственного вращения ротора. Уравнения движения маховика динамически настраиваемого гироскопа, учитывающие угловую податливость скоростной опоры, выглядят следующим образом:
где α1, β1 – углы наклона маховика относительно вала; α2, β2 – углы наклона вала относительно корпуса; ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Первые два из полученных уравнений описывают движение маховика относительно вала в осях, не участвующих в собственном вращении, другие два описывают движение вала относительно корпуса в осях, связанных с корпусом. Анализируя систему (29), можно сделать вывод, что уравнения движения вала и уравнения движения маховика связаны между собой через инерционные и гироскопические составляющие, обусловленные инерционными свойствами кардановых колец. |