Научная статья на тему 'Моделирование движений гироскопа методом электроаналогий'

Моделирование движений гироскопа методом электроаналогий Текст научной статьи по специальности «Математика»

CC BY
565
215
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ СИСТЕМ / МЕТОД ЭЛЕКТРОМЕХАНИЧЕСКИХ АНАЛОГИЙ / ГИРОСКОП / A GYROSCOPE. / MATHEMATICAL MODELLING OF SYSTEMS / A METHOD OF ELECTROMECHANICAL ANALOGIES

Аннотация научной статьи по математике, автор научной работы — Чикуров Н. Г.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Modelling of movements of a gyroscope by a method of electro-mechanical analogies

On an example of a gyroscope with three degrees of freedom the technique of mathematical modelling of mechanical systems is shown by a method of electro-mechanical analogies. The dynamic model of a gyroscope in the form of the differentsialno-algebraic equations is constructed.

Текст научной работы на тему «Моделирование движений гироскопа методом электроаналогий»

Уфа : УГАТУ, 2013 УДК 531.3.07

'Вютньк, ^ГДОиОу

Т. 17, № 2 (55). С. 53-63

Моделирование движений гироскопа методом электроаналогий

Н. Г. Чикуров

[email protected]

ФГБОУ ВПО «Уфимский государственный авиационный технический университет» (УГАТУ)

Поступило в редакцию 22.12.2012

Аннотация. На примере гироскопа с тремя степенями свободы показана методика математического моделирования механических систем методом электроаналогий. Построена динамическая модель гироскопа в виде дифференциально-алгебраических уравнений.

Ключевые слова. Математическое моделирование систем; метод электромеханических аналогий; гироскоп

ВВЕДЕНИЕ

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

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

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

важное следствие, вытекающее из первого закона Кирхгофа: сумма производных токов (применительно к механике - ускорений) в узле также равна нулю. Данное утверждение в терминах механики не столь очевидно, как в электротехнике.

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

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

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

ОПРЕДЕЛЕНИЕ ИНЕРЦИОННЫХ ПАРАМЕТРОВ ГИРОСКОПА

Из-за ограниченного объема статьи рассмотрим простейший гироскоп с тремя степенями свободы в виде твердого тела, одна из точек которого неподвижна (рис. 1). Положение такого тела можно определить в неподвижной (базовой) системе координат путем задания трех направляющих косинусов главной центральной оси инерции (оси симметрии) тела.

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

Выведем уравнения движения рассматриваемого тела в неподвижной системе координат Охуг. Рассмотрим вектор угловой

скорости тела ю = |шх, шу,ш*} и вектор главной оси симметрии г = |гх, гу,г} . Тогда вектор

линейной скорости конца вектора г (точки с ) будет равен [3]

У=юXг = |ШуГ2 - Гуш, ш2Гх - г Шх, ШхГу - ГхШу }.

Рис. 1

Полагая V = ^, приходим к системе

дифференциальных уравнений

где

Лгх

—^ = г ш - г ш = V ,

Л * у у * х’

Лгу

— = гт - г тх = Vv,

Ж х * * х у’

Лг* т,

— = Г т - гт = V;,

Лі у х х у ^

V = Г Шу - ГуШ , Vy = ГхШ - ГШх , Vz = ГуШх - ГхШу.

(1)

(2)

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

М - М" = 0,

(3)

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

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

Лт Л|/||

-----+ —!!~^ю = М"

Лі Л

(4)

где р|| - тензор инерции тела.

Тензор инерции представляет собой

матрицу, вычисленную системы координат

в заданной точке

(

3 - 3 - 3

X ху х

- 3 3 - 3

ух у у*

- 3у 3*

Таким образом, дальнейшее решение задачи сводится к определению входящих в тензор инерции осевых моментов инерции , Jy, Jz и

центробежных моментов инерции Jxy, Jyz, Jx ,

а также их производных.

Определим осевые моменты инерции тела относительно координатных осей базовой системы координат Охуг. Будем считать моменты

инерции тела J'x, J'y, J' определенными в подвижной системе координат Оху 2Г, ось г' кото-

рой совмещена с осью симметрии тела (рис. 2). Поскольку тело симметрично относительно этой оси, то J'x = J' .

* 7

Сначала найдем момент инерции тела Jx относительно оси x базовой системы координат. Координатные оси Oxyz являются главными осями инерции тела. Поэтому центробежные моменты инерции раны нулю и момент инерции Jx относительно орта i, определяющего направление оси x, и составляющего с осями xуz углы соответственно ax, px, gx, определяется по формуле [3]

Jx = J'x cos2 ax + jy cos2 bx + J'z cos2 gx. (5)

Сумма квадратов направляющих косинусов cos2 ax + cos2 px + cos2 gx = 1.

Откуда

cos2 px + cos2 gx = 1 - cos2 ax.

С учетом последнего выражения и равенства J'x = Jy преобразуем уравнение (5).

Jx = J'x (cos2 ax + cos2 px) + J'z cos2 gx =

= J (1 - cos2 gx ) + JZ cos2 gx = ( JZ - J ) cos2 gx + Jx .

где cos gx = —.

r

Определяя аналогично моменты инерции относительно двух других осей у и z , получаем

J =

J- - Jx . 2

гг г*

Тf Тf

rl + JX,

J - J

j=ry+j: ,

j. =

J' - J'

r2 + JX.

Обозначим константу

J' - J' C = J-^-J.

(6)

Тогда последняя система уравнений приобретает вид

J, = сТ + гх,

JУ=о;+^, (7)

Л = + J'x.

Принимая во внимание, что моменты инерции !', J,, J'z в подвижной системе координат

величины постоянные, продифференцируем уравнения (7).

dJx = 2Crx drx_

dt dt

dJy = 2Cry dr,

dt y dt

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

dJz = 2Cr_ drz

dt ~dt ‘

Подставляя в последние уравнения выражения для производных из уравнений (1), получаем

dt

= 2CrxVx,

dJ

—y = 2CrV , dt y y

(8)

J

dt

Определим центробежные моменты инерции 3х , 3 *, 3*х в точке О неподвижной системы

координат. Воспользуемся формулами [3]

,1_у, = шгуу + а,^, (J0 - Jx0) + а12а22 (JZ - J0),

./у: = тгХ + а21а31 (^ ) + а22а32 (^ - !),

Jx = тт>; +а,1аи (J: - J:) + а33ав (J.0 - J°y), где тсх, тус, тсг - координаты центра масс тела.

В силу симметрии тела относительно центральной оси !х = J0y и система уравнений принимает следующий вид

Jxy = тТ>С +(а11а21 +а12а22 )( ^ - ^ ) ,

Jу: = тТу Т + (а21а31 + а22а32 )(г - Jх ) ,

Jzx = тТгТх +(а31а11 +а33а13 )(Jz - Jx ) .

Исходя из ортогональности векторов х0, у0 и г0, имеем

а11а 21 + а12 а 22 + а13а 23 = 0,

а 21 а31 +а 22 а32 +а 23а33 = 0 а31а11 +а32 а12 +а33а13 = 0.

r

Поэтому можем записать

г г

х у

^11^21 +^!2^22 = -^з^з _ 2 >

Г

, _ _ ГуГ

^21^31 +^22^32 ^23^33 2 ,

Г

ГГ

«31«11 + а32а12 _ а33а13 _ 2 •

Г2

Подставим эти выражения в уравнения для центробежных моментов инерции.

гг

■у = тг>; - ру (■ - ■),

гг

■, = тГуГ, - ^ ( ■ - ■ ).

■х = тгХ - ^ (■ - ■).

(9)

Поскольку центральная ось симметрии ъ0 проходит через точку О неподвижной системы координат, то - тг2, = ■у - тг2,

. Подставляя эти соотношения в уравнения (9) и полагая координаты центра масс тела равными г. = г,, г. = гу, г. = гг , получаем компактную форму уравнений для центробежных моментов инерции.

■ху = -Сг,гу , ■уъ = -Сгуг, , = -Сгъгх,

(10)

где коэффициент С определяется с помощью формулы (6).

Теперь найдем производные центробежных моментов. Для этого продифференцируем уравнения (10).

( Лг Лгу

ху -С х гу + У_

Лі ч Лі у Лі

( Лгу Лгъ

уг -С у г2 +

Лі ч Лі Лі

■ -С ( Лг Лг

ъх 2 гх +

Лі V Лі Лі

Подставим выражения для производных из уравнений (1) в последние уравнения. В результате получим

■ __х

Лі ■ ___у

Лі

Лі

= -С(V г + V г )

\ х у у х )

= -С (V г + V г )

\ у 2 2 у )

= -С (У,г, + Кп)

(11)

Вернемся к векторному уравнению (4). Запишем его подробно

(Лю ^

Лі

Лі Л ю.

- 3ху -

- 3ух ■у - ^

- -

Лі

■х ■ху

Лі Лі Лі

■ух Л1у ■уъ

Лі Лі Лі

■х ■у

Лі Лі Лі

ю

(М; Л и; ми

V 2у

Или в координатной форме

М; = ■ ю - ■ ю - ■ ю + ■ ю - ■ ю - ■ ю ,

х хх ху у х, ъ хх ху у X, ъ ^

М; =-■ ю+ ■ ю - ■ ю - ■ ю+ ■ ю - ■ ю

у ух х у у уъ ъ ух х у у у ъ

М; =-■ ю -■ ю + ■ ю -■ ю -■ ю + ■ ю .

ъ ъх х ъу у г г ъх х ъу у г г

(12)

РЕШЕНИЕ ЗАДАЧИ МЕТОДОМ ЭЛЕКТРОАНАЛОГИЙ

Чтобы построить математическую модель гироскопа с применением метода электроаналогий, необходимо составить уравнения кинетостатики для главных векторов и для главных моментов активных сил, реакций связей и сил инерции, действующих на гироскоп. Все моменты должны быть вычислены относительно полюса, совмещенного с центром масс гироскопа. Тогда результирующее движение гироскопа, происходящее под действием указанных выше сил и моментов, можно рассматривать как сумму поступательного движения гироскопа с его вращательным движением вокруг центра масс.

Все векторные уравнения будем проектировать на оси подвижной системы координат Сх1 у^, начало которой совмещено с центром масс гироскопа, т. е. с точкой С , а оси параллельны соответствующим осям системы Охуг (рис. 3).

Силы реакции Кх, Ку, Кг в опорной точке О

гироскопа теперь следует рассматривать как внешние силы. Моменты, создаваемые этими силами относительно координатных осей х, у и г, зависят от величин проекций гх, гу и г2 вектора оси г гироскопа.

Осевые моменты инерции гироскопа относительно осей подвижной системы координат

г

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

г

г

у

но а13 = —, а23 =—, а33 = —

г

(

+

V

/

г

х

\

г

у

г

г

Сх1 у1 ^ будем рассчитывать с помощью формул, которые выводятся аналогично формулам (7).

4х = С,г; + /0,

/ у = Сг2 + /I, (13)

/„ = сг + /0,

где /°, /0, / - моменты инерции относительно главных центральных осей х0, у0, ,0; Й1 - постоянный коэффициент.

J0 - J0

C = J-^-J

(14)

Рис. 3

Уравнения для производных моментов инерции (8), уравнения для центробежных моментов инерции (10) и уравнения для производных центробежных моментов (11) в системе координат Cx1 y zx, отличаются лишь постоянной C,.

Рассмотрим движение гироскопа в каждой из трех плоскостей, образованных системой координат Cxyz.

Движение в плоскости Cxy

Рассмотрим движение, которое совершает гироскоп в плоскости Cxy (рис. 4). В соответствии с данной схемой запишем уравнения равновесия сил и моментов. Уравнения равновесия сил для оси, соединяющей две смежные координатные плоскости, будут одинаковы для этих плоскостей. Чтобы исключить дублирование уравнений, каждое из уравнений сил и моментов отметим соответствующей литерой, номер которой свяжем с индексом координатной оси.

Рис. 4

Пусть оси x соответствует номер 1 литеры, оси y - номер 2 литеры, а оси z - номер 3 литеры. Буквенно-цифровые символы литеров будем записывать с левой стороны уравнений. В результате уравнения сил и моментов для плоскости Cxy примут следующий вид:

C1) Rx - Fxu = 0, с 2) Ry - f; = 0,

M3) - MRy + MR - MI = 0, где MR = Rxry, MR = Ryrx.

Далее запишем компонентные уравнения, номера литеров для которых свяжем с индексами осей

K1) Fl = mVx,

K2) F; = mVy.

В качестве компонентных уравнений для инерционных моментов используем уравнения (12). Перепишем третье из них в сжатом виде. Для данных уравнений связанный с индексом оси номер литеры надо увеличить на 3

K 6) MU = -j1zx wx - j1; Wy + j1z wz - j1zx Wx - j1; Wy + j1z Wz.

Входящие в уравнение осевой момент инерции J1z и его производную j1z вычислим с помощью формул (13) и (8). Центробежные моменты инерции j1zx , j1zy и их производные j 1zx

и j 1zy необходимо вычислить для центра масс

гироскопа, т. е. в точке с, с помощью формул (10) и (11). Принять в формулах (8), (10), (11) коэффициент C = с .

На основании записанных уравнений равновесия и компонентных уравнений строим эквивалентную электрическую схему замещения (рис. 5) [2]. Электрические величины на схеме будем обозначать с помощью соответствующих символов, принятых в механике. Эта схема строится в два этапа. Вначале следует вычертить нижние два контура, в которых напряжения и ЭДС. должны соответствовать уравнениям равновесия сил C1 и C2 . В соответствии со вторым законом Кирхгофа сумма напряжений и

r

ЭДС. в каждом контуре должна быть равна нулю.

На втором этапе чертим верхний контур, отвечающий уравнению равновесия моментов М3 . Чтобы выполнить содержащиеся в уравнении М 3 операции умножения, в электрическую схему введены два идеальных трансформатора с коэффициентами трансформации, соответствующими плечам действующих сил.

М‘

Рис. 5

Запишем уравнения идеальных трансформаторов, входящих в рассматриваемую электрическую схему, причем соотношения токов продифференцируем

71.1) ЯхГу = Му,

71.2) ю,Гу = V?,

71.3) юг +ю7Гу =“7гу + “7Уу = У“2,

7 У 7 у 7 у 7 у ху

Т2.1) Я/х = му,

Т2.2) о>л = У™,

Т2.3) ю7Гх +ю7Гх =ю7Гх +ю7Кх = V“7.

/ 7 х 7 х 7 х 7 х у

Движение в плоскости Схг

Рассмотрим движение, которое совершает гироскоп в плоскости Схг (рис. 6).

Рис. 6

Используя установленные выше правила для присвоения номеров литерным уравнениям, запишем уравнения равновесия сил и моментов

С1) у -Б = 0,

С3) - Р + у - ¥и2 = 0,

М 2) - Му + Му - М; = 0,

где Мух = Яхг2 , М^ = Я7гх , в свою очередь, гх , г7 -

проекции вектора г на координатные оси. Запишем компонентные уравнения.

К1) Б^и = тУх,

К 3) Б? = шУг,

К5)М; = /“х + /1у “у - /1 у7 “7 - /1ух “х + /1уЮу - /1у7“7.

На основании записанных уравнений равновесия и компонентных уравнений строим эквивалентную электрическую схему замещения (рис. 7).

. К .

Рис. 7

Уравнения идеальных трансформаторов Т3.1) Яхг2 = М ух,

Т3.2) шугг = Г“у,

Т3.3) “угг +ШуГ7 =“уГ7 +Шу V = У“,

Т4.1) Яг = Му,

Т4.2) ШуГх = V“,

Т4.3) “г + шуг =“г + = У7щу.

х х х х 7

Движение в плоскости Суг

Рассмотрим движение гироскопа в плоскости Суг (рис. 8).

V, А

Рис. 8

Запишем уравнения сил и моментов.

с2) яу - б; = 0,

С3) - Р + - Б7и = 0,

М1) - му + му - М; = 0, где МУ = ЯуГ7, МУ = Яг,

Т5.1) яуг2 = му,

Т5.2) Шхг2 = Ку“,

т 5.3) юхгъ + юхгъ =юхгъ + ЩЛ = ^

Т6.1) Кг = му

г К - х ,

в свою очередь, гу ,гъ

у ъ ’ х ъ у

проекции вектора г на координатные оси. Компонентные уравнения

К 2)Б; = mV у,

К 3)Б; = тУ2,

К 4)М = Лх “х - /1ху “у - ^хЛ + А^х“х - /1ху “у - А “7 -

На основании записанных уравнений равновесия и компонентных уравнений строим электрическую схему замещения (рис. 9).

Рис. 9

Уравнения идеальных трансформаторов

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Т6.2) юхгу = Vr,

Т6.3) юхгу + Шх^ = юхгу + ЮxVy = VWX.

Общая электрическая схема замещения

Объединяя полученные выше схемы замещения, построим общую электрическую схему замещения гироскопа (рис. 10).

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

Рассмотренный переход от эквивалентных электрических схем замещения к натуральным электрическим схемам имеет лишь теоретический интерес.

При построении математических моделей такой переход не требуется, поскольку из тео-

Рис. 10. Общая электрическая схема замещения гироскопа

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

Схема замещения (рис.10) содержит узлы 1, 2 и 3. Запишем для этих узлов уравнения по первому закону Кирхгофа и продифференцируем эти уравнения.

Кр1.1) - Ух - У™ + У“ = 0,

Кр1.2) - Ух - У“ + У“ = 0,

Кр2.1) - уу - уу°х +у“ = 0,

= 0,

(15)

Кр2.2) - Уу - У™ + уу;

Кр3.1) - У2 -У“у + У“х = 0,

Кр3.2) - У г - У“ + У “х = 0.

Законы Кирхгофа позволяют составить полную систему линейных уравнений, из которой могут быть найдены все неизвестные переменные на входах интегрирующих звеньев, входящих в математическую модель. В данной системе используются девять интегрирующих звеньев, т. е. это система 9-го порядка. Запишем дифференциальные уравнения указанных интегрирующих звеньев.

Е1) р“х = “х , Е6) РГ = К ,

Е2) рщу =йу, Е7) рУх = Ух,

Е3) р“1=“1, Е8) рУу = У,, (16)

Е4) ргх = Уу, Е9) рУ2 = У 2.

Е5) ргу = У,,

Далее в соответствии с процедурой метода электроаналогий построим сеть связей (рис. 11) [2]. Каждая стрелка в сети связей отмечена литерой соответствующего уравнения, с помощью которого будет определяться переменная, на которую указывает стрелка. Сеть связей используется как основа для построения системы линейных алгебраических уравнений (СЛАУ), которая имеет следующий вид:

М и = В,

где и - вектор неизвестных переменных, В -вектор свободных членов, М - матрица коэффициентов системы.

Решение СЛАУ определяет правые части системы дифференциальных уравнений (16).

Сеть связей представляет собой прототип структурной схемы динамической системы. За-

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

Анализируя сеть связей, выделим из нее СЛАУ. Для упрощения записей введем обозначения переменных:

/1х

/1ху = -С1ГЛ >

/1У = сГ+/0

С1ГУГ7

/17х =-СгГх,

/ =-С1 (УхГу + УуГх),

= -С1 (УуГ7 + У/у),

= -С1 (КГх+УхГ).

Рис. 11. Сеть связей

Можно сократить число алгебраических уравнений гироскопа путем исключения внутренних переменных. Объединяя уравнения

Кр1.2, Кр2.2, Кр3.2, М1, М2, М3 с трансформаторными уравнениями, а уравнения С1, С2, С3 с компонентными уравнениями К1, К2, КЗ, придем к системе из 12 линейных уравнений.

Запишем эти уравнения в виде следующего реестра:

и4 = м; и5 = м; и 6 = м;

и7 = К* и8 = Ку

и10 = “х

и11 = “у

-V, - ю2Гу + ю/2 = ю2Гу - ЮуГ2,

-V у - “Л + юЛ = - ^,

-¥; - ЮуГ, + Ю,Гу = ЮуГ, - Ю,Гу ,

-кг + Дуг2 - м; = о,

-Дл + кЛ - м; = о,

-Дуг + К,Гу -м; = о,

к - А = 0,

Ку - ^у = 0

з,- А“у - -м; =

-/1х<Их + / “у + 31Л,

-/1ху “ + /1у“у - /1yz“ -м; = /у“ - / “у + / “ ,

/ “ - А “у + “ - м; =

-/1^ “ + / “у - А “ .

Соответственно вектор свободных членов принимает вид:

В = WVy ^у^ ,

В2 = ®Л -“V,

В =Юу^ -ю,Ку,

В = Р; ,

В10 = - Ах“х + Аху“у + ,

В11 = Аху “х - Ау “у + Ау “ ,

В12 = ■их “х + </1уг “у - ■ “ .

С учетом новых обозначений переменных перепишем систему дифференциальных уравнений (16)

р II р^1 = и 10..

р у II 8- у и11 II р

PЮz = Юz ь2 II р

II ь" II ¡хг р

II О Ь” II р

II С II р

II II р

II ь II р

II ь" II р

РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ

В целом мы получили совместную систему, которая включает в себя 12 линейных уравнений и 9 дифференциальных уравнений.

Численное решение системы ДУ на каждом шаге интегрирования включает несколько стадий, каждая из которых предусматривает вычисление правых частей системы ДУ при различных значениях аргумента. Значения сигналов на входе каждого интегратора математической модели определяются с помощью СЛАУ, которая решается на каждом шаге интегрирования ДУ столько раз, сколько стадий содержит алгоритм численного решения системы ДУ.

Например, если использовать для решения системы ДУ численный метод Рунге-Кутты 4го порядка, то на каждом шаге интегрирования надо четыре раза решать систему линейных алгебраических уравнений.

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

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

т = 1.0, К = 0.1, г = 0.2, Н = 0.05, “ = 500,

0 = /0, / = 1 тК2, Ф = Р. ч12 4 ) у х 2 3

На рис. 12 представлены графики движений гироскопа. На этом рисунке показано: а, б, в -изменения проекций оси гироскопа, г, д, е - изменения реакций связи в точке N . Время интегрирования дифференциальных уравнений выбрано равным одному периоду прецессии.

Из теории гироскопов известно, что этот период определяется по формуле [1, 3]

2/“

Р • Г

Графики движения подтверждают, что период прецессии испытуемого гироскопа соответствует данной формуле.

Кроме основного движения прецессии в гироскопе происходят высокочастотные нутационные колебания.

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

/ = т I — Н2 +1 К2

Rz -тУ z = Pz

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

На основе данной модели можно построить модели других, более сложных гироскопов. Они накладываются на все основные движения гироскопа. На графиках линейных скоростей V и V

(рис. 13, а, б) видны огибающие скоростей, на которые наложены высокочастотные нутационные колебания. В увеличенном масштабе эти графики показаны на рис. 13, в, г. Нутационные колебания влияют также на величины проекций оси гироскопа (рис. 13, д, е).

Рис. 13. Нутационные колебания гироскопа

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ЗАКЛЮЧЕНИЕ

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

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Пельпор Д. С. Гироскопические системы. Теория гироскопов и гироскопических стабилизаторов: учеб. для вузов. 2-е изд. М.: Высш. шк., 1986. 423 с.

2. Чикуров Н. Г. Синтез математических моделей технических систем методом электроаналогий // Вестник УГАТУ. 2009. Т. 12, № 2 (31). С. 156-165.

3. Бутенин Н. В., Лунц Я. Л., Меркин Д. Р. Курс теоретической механики: в 2 т. СПб.: Лань, 2007. 736 с.

4. Введение в математическое моделирование: учеб. пособие / Под ред. П. В. Трусова. М.: Логос, 2004. 440 с

ОБ АВТОРЕ

Чикуров Николай Георгиевич, доц., каф. автоматизации технол. процессов. Дипл. инж. по станкостроен. и технол. машиностроен. (УАИ, 1965). Канд. техн. наук по электро-оборуд. металлореж. станков (МОССТАНКИН, 1971). Иссл. в обл. комп. систем ЧПУ класса ГСЫС.

METADATA

Title: Modelling of movements of a gyroscope by a method of electro-mechanical analogies.

Authors: N. G. Chikurov

Affiliation: Ufa State Aviation Technical University (UGATU), Russia.

Email: [email protected].

Language: Russian.

Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 17, no. 2 (55), pp. 53-63, 2013. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).

Abstract: On an example of a gyroscope with three degrees of freedom the technique of mathematical modeling of mechanical systems is shown by a method of electro-mechanical analogies. The dynamic model of a gyroscope in the form of the differential algebraic equations is constructed.

Key words: Mathematical modeling of systems; a method of electromechanical analogies; a gyroscope.

References (English Transliteration):

1. D. S. Pelpor, Gyroscopic Systems. The Theory of Gyroscopes and Gyroscopic Stabilizers, (in Russian). Moskow, Vysshaya shkola, 1986.

2. N. G. Tchikurov, “Mathematical simulation of technical systems by a method of electroanalogies,” (in Russian), Vestnik UGATU, vol. 12, no. 2 (31), pp. 156-165, 2009.

3. N. V. Butenin, Y. L. Lunz, D. R. Merkin, Course of Theoretical Mechanics, (in Russian). St.-Petersburg: Lan’, 2007.

4. Introduction in Mathematical Modelling, (in Russian), P. V. Trusov Ed. Moskow: Logos, 2004.

About author:

Chikurov, Nikolai Georgievich, Technological Processes Automation Dept., Cand. (PhD) Tech. Sci. Researches in the field of control systems on the basis of computers.

i Надоели баннеры? Вы всегда можете отключить рекламу.