Главная Учебники - Разные Лекции (разные) - часть 62
Факультет
Информатики
и систем управления.
Кафедра
ИУ-2.
Расчетно-пояснительная
записка
к
дипломной
работе на тему:
“Двухосный
индикаторный
гиростабилизатор
телекамеры
на ВОГ” Студент (Носов
Н.А.) Руководитель
проекта (Фатеев
В.В.)
1996 г.
Телевизионная
техника применяется
в различных
областях человеческой
деятельности
- экономике,
искусстве,
военном деле
и многих других.
Область ее
применения
постоянно
расширяется.
Это объясняет
активное развитие
в настоящее
время телевизионного
и кинотехнического
оборудования
и применение
новейших достижений
техники при
разработке
и производстве
устройств для
различных видов
теле и киносъемки.
К таким
устройствам
относят гироскопические
стабилизаторы
(ГС) углового
положения
телекамеры,
которые применяются
при съемках
с подвижных
объектов: вертолета,
судна, автомобиля,
а также с кранов.
Из требований,
предъявляемых
к этим устройствам,
наиболее важными
являются следующие:
— высокая
точность
стабилизации,
обусловленная
применением
объективов
с большим фокусным
расстоянием;
— дистанционное
управление
стабилизатором
и телекамерой,
что связано
с целесообразностью
или даже необходимостью
размещения
телекамеры
в месте, недоступном
для оператора,
например, на
кране;
— защита
от вибрационных
возмущений
как с целью
обеспечения
виброустойчивости
гиростабилизатора,
так и для повышения
качества отснятого
материала,
снижающегося
при поступательном
движении телекамеры;
— удобство
и простота в
обращении,
необходимость
“защиты от
дурака”;
— минимально
возможные
габариты и
масса;
— удобство
доступа к
стабилизируемой
камере;
— высокие
скорости (до
120 град./с.) и ускорения
(до 200 рад./с.) управления,
для получения
“динамичного”
отснятого
материала;
— большие
углы прокачки,
по некоторым
осям равные
360°;
— минимально
возможное
энергопотребление,
т.к. для работы
часто используются
автономные
источники
питания (аккумуляторы
и т.д.).
— большой
диапазон
балансировок,
необходимость
которого вызывается
использованием
различных типов
теле и кинокамер
с сильно различающимися
массогабаритными
показателями.
Из устройств,
в наибольшей
степени удовлетворяющих
перечисленным
требованиям,
в настоящее
время известны
и применяются
следующие. Это
силовой гироскопический
стабилизатор
"Wesscam" (Канада) и
трехосные
гиростабилизаторы
1ГСП (разработка
МВТУ и НИКФИ),
2ГСП, 4ГСП (разработка
МВТУ {МГТУ}).
Первый
представляет
собой платформу,
помещенную
в карданов
подвес, на которой
установлены
три двухстепенных
гироскопа.
Компенсация
внешнего момента
осуществляется
разгрузочным
двигателем
и маятниковым
устройством.
Сам карданов
подвес крепится
к основанию
с помощью упругой
связи, необходимой
для аммортизации
устройства
при поступательных
вибрациях
вертолета,
автомашины
и т.п. Гиростабилизатор
помещен в защитную
сферу имеющую
вырез напротив
объектива
киноаппарата.
Гиростабилизаторы
1ГСП, 2ГСП, 4ГСП
индикаторного
типа. Все они
в качестве
чувствительного
элемента используют
трехстепенный
гироскоп типа
МГТУ-05. Также
все эти три
гиростабилизатора
имеют дистанционное
управление
киноаппаратом
и пространственным
положением
стабилизированной
платформы.
Особенностью
гиростабилизатора
1ГСП является
использование
внутреннего
карданова
подвеса и
двухконтурной
системы стабилизации
с применением
маховиков
совместно с
двигателями
стабилизации.
К недостаткам
этого стабилизатора
относятся
низкая скорость
управления
(10 град/сек) и
отсутствие
защиты от
аэродинамических
воздействий.
В гиростабилизаторе
2ГСП применен
наружный карданов
подвес с подшипником
большого диаметра
и также используется
двухконтурная
система стабилизации
(маховик и двигатели
стабилизации).
Скорость управления
до 30 град/сек.
Для защиты от
аэродинамических
воздействий
введен следящий
частично прозрачный
обтекатель,
внутри которого
расположен
трехосный
аммортизатор
для предохранения
платформы с
киноаппаратом
от линейных
вибраций.
Гиростабилизатор
4ГСП также выполнен
с наружным
кардановым
подвесом. Он
имеет скорость
управления
до 60 град/сек.,
ограниченную,
в основном,
скоростью
управления
примененного
чувствительного
элемента. Защита
от аэродинамических
нагрузок
отсутствует.
В данной
работе продолжены
исследования
по возможности
построения
гироскопического
стабилизатора
углового положения
телекамеры,
в котором:
—скорость
управления
обеспечивается
не менее 100 град/сек;
—предусмотрено
дистанционное
управление
телекамерой
и самим гироcтабилизатором;
—в качестве
чувствительного
элемента использован
волоконно-оптический
гироскоп.
В частности,
рассматриваются
вопросы:
—обеспечения
устойчивости
канала стабилизации
при существенно
нежесткой
конструкции
ГС, нежестком
креплении
телекамеры
к ГС и расположении
чувствительного
элемента не
на стабилизируемом
объекте;
—проведено
исследование
инерционных
возмущающих
моментов, в том
числе моментов
возникающих
от несимметричности
конструкции
рам ГС;
— проводилась
доработка
усилителя
мощности с ШИМ;
— предложена
конструкция
датчика угла
фазового типа;
— разработана
конструкция
двухосного
ГС.
Одним
из основных
факторов,
определяющих
выбор принципиальной
схемы гиростабилизатора
телекамеры,
является тип
карданова
подвеса. В
гиростабилизаторе
телекамеры
может быть
использован
как внутренний,
так и наружный
карданов подвес.
Сравнительный
анализ и конструктивная
проработка
схем подвесов
[3] показывает,
что применение
наружного
карданова
подвеса для
стабилизатора
телекамеры
представляется
более целесообразным.
Это объясняется
следующими
причинами.
Использование
наружного
карданова
подвеса при
больших углах
прокачки позволяет
получить более
компактную
конструкцию
стабилизированной
платформы. В
этом случае
момент инерции
платформы
относительно
собственной
оси вращения
значительно
снижается, и
тогда величина
максимального
момента двигателя
стабилизации,
выбираемая
из условия
обеспечения
необходимого
углового ускорения
платформы при
управлении
может быть
уменьшена. Это
позволяет
повысить точность
стабилизации
за счет использования
двигателя
стабилизации
меньших габаритов,
имеющего меньший
момент сухого
трения вокруг
оси вращения
и меньший коэффициент
демпфирования.
Габаритные
размеры гиростабилизатора
телекамеры
с наружным
кардановым
подвесом оказываются
меньше, чем с
внутренним,
т.к. в последнем
случае для
получения
достаточных
рабочих углов
поворота платформы
необходимо
выполнение
подвеса по
гантельной
схеме, что приводит
к значительному
увеличению
одного габаритного
размера гиростабилизатора
по сравнению
с другим. Применение
наружного
карданова
подвеса позволяет
добиться минимального
различия между
габаритными
размерами
гиростабилизатора
по взаимноперпендикулярным
осям, что является
желательным.
В то же
время схема
гиростабилизатора
с наружным
кардановым
подвесом имеет
следующие
недостатки:
— увеличенный
возмущающий
инерционный
момент, действующий
вокруг оси
наружной рамы
подвеса, который
возникает при
переносных
поворотах
основания
гиростабилизатора;
— пониженная
жесткость рам
наружного
карданова
подвеса по
сравнению с
внутренним.
Возмущающий
инерционный
момент, действующий
по оси наружной
рамы карданова
подвеса, вызывает
появление
ошибок стабилизации,
а также создает
дополнительные
нагрузки на
приводы. Однако,
как показывают
исследования,
в реальных
условиях эксплуатации
гиростабилизатора
киноаппарата
на кране и на
вертолете
[8,9], величина
инерционного
момента при
симметричной
конструкции
рам оказывается
незначительной.
В связи с этим
первый из
перечисленных
недостатков
наружного
карданова
подвеса оказывается
несущественным.
В представленном
двухосном
гиростабилизаторе
телекамеры
применен наружный
карданов подвес.
Конструкция
подвеса гиростабилизатора
позволяет
получить угол
прокачки по
оси тангажа
+60...-80 град., а по оси
курса угол
вращения не
ограничен.
Описание
особенностей
конструкции
гиростабилизатора.
Особенностью
данного ГС
является выполнение
наружной рамы
в виде Г-образной
конструкции.
Это позволяет
уменьшить
габариты ГС
и упростить
доступ к телекамере.
Однако такая
конструкция
является существенно
несимметричной,
что вызывает
появление
дополнительных
возмущающих
моментов из-за
значительных
по величине
центробежных
моментов инерции
рам. Исследование
этих возмущающих
моментов проведено
в разделе “Анализ
инерционных
возмущающих
моментов”.
Кроме
того, требования
по минимизации
масс, моментов
инерции, требования
по динамике
управления
платформой
приводят к
тому, что наружная
рама оказывается
существенно
нежесткой. А
так как при
этом необходимо
учитывать
требования
по высокой
точности стабилизации
при значительных
возмущающих
моментах, то
возникает
необходимость
в проведении
специальных
исследований
по вопросам
обеспечения
устойчивости
канала стабилизации.
В частности,
исследований
по расположению
ЧЭ в конструкции
ГС. Исследование
устойчивости
канала стабилизации
приведено в
разделе “Исследование
влияния нежесткостей
элементов
гиростабилизатора
на его устойчивость”.
Особенностью
данной конструкции
ГС является
то, что стабилизация
положения
телекамеры
по курсу осуществляется
косвенным
образом, путем
стабилизации
положения
наружной рамы
карданова
подвеса. Эта
особенность
также учтена
в разделе
“Исследование
влияния нежесткостей
элементов ГС
на его устойчивость”.
Из требований
по минимизации
энергопотребления
вытекает
необходимость
увеличения
КПД канала
стабилизации.
Согласно этому
требованию,
а также с целью
получения
значительных
по величине
моментов привода,
в качестве
привода используется
редукторный
привод со встречным
включением
двигателей
стабилизации
серии ДПР. Кроме
того, использование
в качестве ЧЭ
волоконно-оптического
гироскопа
позволяет
снизить энергопотребление
собственно
ЧЭ до 3 Вт, как
в режиме измерения,
так и при управлении
положением
стабилизированной
платформы.
Наиболее
существенное
влияние на КПД
электронной
части канала
стабилизации
оказывает
коэффициент
полезного
действия УМ.
Поэтому УМ
выполнен импульсным,
с использованием
ШИМ модуляции
выходного
напряжения.
Это позволяет
примерно в два
раза увеличить
КПД УМ по сравнению
с линейными
схемами УМ.
Однако все
импульсные
УМ являются
мощными источниками
электромагнитных
помех, поэтому
в данной конструкции
ГС УМ располагается
на самом ГС, в
непосредственной
близости от
двигателей
стабилизации.
Кроме того,
непосредственно
на ГС расположены
схемы защиты
ВОГа.
Конструкция
крепления
телекамеры
позволяет
проводить
установку на
платформу
телекамер
отличающихся
по массогабаритным
параметрам
от базовой на
±
30 %. При этом
осуществляется
независимая
регулировка
положения
телекамеры
по трем взаимоперпендикулярным
осям.
Применение
в качестве ЧЭ
ВОГа вместо
механических
гироскопов
позволяет
практически
снять ограничения
по максимальным
скоростям
измерения и
управления,
накладываемых
на канал стабилизации
самим ЧЭ.
При
несимметричной
конструкции
рам гиростабилиза-тора
и значительных
угловых скоростях
движения основания
и управления
платформой
необходимо
учитывать
возмущающие
моменты, вызываемые
осевыми и
центробежными
моментами
инерции рам. В
данной работе
проводится
исследование
инерционных
возмущающих
моментов для
двухосного
гиростабилизатора,
с учетом влияния
центробежных
моментов инерции
рам и скоростей
управления
платформой.
Выражения
для инерционных
моментов получены
путем раскрытия
членов, зависящих
от параметров
движения основания
и платформы
входящих в
динамические
уравнения
Эйлера. Основные
математические
преобразования
выполнялись
с помощью программы
“DERIVE”. Системы
координат и
обозначения
используемые
далее. Рис.1.
X0,Y0,Z0
- система координат
связанная с
основанием. X1,Y1,Z1 -
система координат
связанная с
наружной
рамой.
X2,Y2,Z2 -
система координат
связанная с
платформой.
Qij -
момент количества
движения j-го
тела по i-й
оси.
wij -
угловая скорость
j-го тела по i-й
оси.
wij' -
угловое ускорение
j-го тела по i-й
оси.
Ji -
осевые моменты
инерции тела
относительно
i-й
оси.
Jij -
центробежные
моменты инерции.
Mij -
внешние возмущающие
моменты действующие
на j-е
тело по i-й оси.
a -
угол поворота
наружной рамы
по оси Y1.
a' -
угловая скорость
вращения наружной
рамы по
оси Y1.
a'' -
угловое ускорение
наружной рамы
по оси Y1.
b -
угол поворота
платформы по
оси Z2.
b' -
угловая скорость
вращ. платформы
по оси Z2.
b'' -
угловое ускорение
платформы по
оси Z2.
Динамические
уравнения
Эйлера для i-го
тела имеют вид: dQxi/dt
- QyiЧwzi
+ QziЧwyi
= Mxi dQyi/dt
- QziЧwxi
+ QxiЧwzi
= Myi dQyi/dt
- QziЧwxi
+ QxiЧwzi
= Myi
В случае
двухосного
гиростабилизатора
эти уравнения
преобразуются
в следующую
форму: а)
для наружной
рамы:
dQy1/dt
- Qz1Чwx1
+ Qx1Чwz1
= My1
б)
для платформы:
dQx2/dt
- Qy2Чwz2
+ Qz2Чwy2
= Mx2 dQy2/dt
- Qz2Чwx2
+ Qx2Чwz2
= My2
(1) dQz2/dt
- Qx2Чwy2
+ Qy2Чwx2
= Mz2
Полный
момент количества
движения наружной
рамы в проекциях
на оси X1,
Y1,
Z1
определяется
следующими
выражениями: Qx1
= Jx1Чwx1
- Jxy1Чwy1
- Jxz1Чwz1 Qy1
= Jy1Чwy1
- Jyx1Чwx1
- Jyz1Чwz1 (2) Qz1
= Jz1Чwz1
- Jzx1Чwx1
- Jzy1Чwy1
Полный
момент количества
движения платформы
в проекциях
на оси X2,
Y2,
Z2
определяется
следующими
выражениями: Qx2
= Jx2Чwx2
- Jxy2Чwy2
- Jxz2Чwz2 Qy2
= Jy2Чwy2
- Jyx2Чwx2
- Jyz2Чwz2 (3) Qz2
= Jz2Чwz2
- Jzx2Чwx2
- Jzy2Чwy2
Кинематические
уравнения
двухосного
гиростаби-лизатора,
для расположения
координатных
осей приве-денного
на рис.1, имеют
вид: а)
для наружной
рамы: wx1
= wx0Чcos(a)
- wz0Чsin(a) wy1
= wy0
+ a' (4*) wz1
= wx0Чsin(a)
+ wz0Чcos(a) wx1'
= wx0'Чcos(a)
- wz0'Чsin(a) wy1'
= wy0'
+ a'' (4*') wz1'
= wx0'Чsin(a)
+ wz0'Чcos(a) б)
для платформы:
wx2
= wx1Чcos(b)
+ wy1Чsin(b) wy2
= wy1Чcos(b)
- wx1Чsin(b) (5*) wz2
= wz1
+ b' wx2'
= wx1'Чcos(b)
+ wy1'Чsin(b) wy2'
= wy1'Чcos(b)
- wx1'Чsin(b) (5*')
wz2'
= wz1'
+ b'' Из
2-го уравнения
в (5*) следует, что: wy1=wx1Чtg(b)+wy2/cos(b) Из
2-го уравнения
в (5*') следует, что: wy1'=wx1'Чtg(b)+wy2'/cos(b)
Тогда,
учитывая, что
wy2,
wz2,
wy2',
wz2'
являются параметрами
движения
стабилизированного
объекта, т.е.
заданы, кинематические
уравнения можно
переписать
в следующем
виде: wx1
= wx0Чcos(a)
- wz0Чsin(a) wy1
= wx1Чtg(b)+wy2/cos(b) (4) wz1
= wx0Чsin(a)
+ wz0Чcos(a) wx1'
= wx0'Чcos(a)
- wz0'Чsin(a) wy1'
= wx1'Чtg(b)+wy2'/cos(b) (4') wz1'
= wx0'Чsin(a)
+ wz0'Чcos(a)
wx2
= wx1Чcos(b)
+ wy1Чsin(b) (5) wx2'
= wx1'Чcos(b)
+ wy1'Чsin(b) (5')
Подставляя
выражения для
полных моментов
количества
движения (2), (3) в
динамические
уравнения
Эйлера (1), получаем
следующий вид
уравнений
движения наружной
рамы и платформы:
Jy1Чwy1'
+ (Jx1-Jz1)Чwx1Чwz1
+ Jzx1Чwx12
- Jxz1Чwz12
+
+
Jzy1Чwx1Чwy1
- Jxy1Чwy1Чwz1
- Jyx1Чwx1'
- Jyz1Чwz1'
= My1 (6.1)
Jx2Чwx2'
+ (Jz2-Jy2)Чwy2Чwz2
- 2ЧJzyЧwy22
+ Jyz2Чwz22
+
+
Jyx2Чwx2Чwz2
- Jzx2Чwx2Чwy2
- Jxz2Чwz2'
- Jxy2Чwy2'
= Mx2 (6.2)
Jy2Чwy2'
+ (Jx2-Jz2)Чwx2Чwz2
+ Jzx2Чwx22
- Jxz2Чwz22
+
+
Jzy2Чwx2Чwy2
- Jxy2Чwy2Чwz2
- Jyx2Чwx2'
- Jyz2Чwz2'
= My2 (6.3)
Jz2Чwz2'
+ (Jy2-Jx2)Чwx2Чwy2
+ Jxy2Чwy22
- Jyx2Чwx22
+
+
Jxz2Чwy2Чwz2
- Jyz2Чwx2Чwz2
- Jzx2Чwx2'
- Jzy2Чwy2'
= Mz2 (6.4)
При
отсутствии
моментов внешних
сил правые
части уравнений
(6.2), (6.3), (6.4) обращаются
в нуль, а правая
часть (6.1) представляет
собой момент
реакции со
стороны платформы
на внешнюю раму
вокруг оси Y1.
Обозначив
левые части
уравнений
(6.1), (6.2), (6.3) буквами
A, B и C, соответственно,
получаем выражение
для полного
инерционного
момента относительно
оси внешней
рамы: My1ин
= A + B Ч
sin(b)
+ C Ч
cos(b) (7)
Раскрыв
в (7) сокращения
A, B и C и преобразовав
получаем выражение
для полного
инерционного
момента Мy1ин.
Мy1ин=Jxz1·{wx12-wz12}+
+Jxz2·cos(b)·wx22-Jyz2·sin(b)·wy22+
+{Jyz2·sin(b)-Jxz2·cos(b)}·wz22+
+{Jyz2·cos(b)-Jxz2·sin(b)}·wx2·wy2+
+{Jxy2·sin(b)+(Jx2-Jz2)·cos(b)}·wx2·wz2+
+{(Jz2-Jy2)·sin(b)-Jxy2·cos(b)}·wz2·wy2+ (8)
+{Jx2·sin(b)-Jxy2·cos(b)}·wx2'
+
+{Jy2·cos(b)-Jxy2·sin(b)}·wy2'-
-{Jxz2·sin(b)+Jyz2·cos(b)}·wz2'+
+Jyz1·wx1·wy1-
-Jxy1·wz1·wy1+
+(Jx1-Jz1)·wx1·wz1
-
-Jxy1·wx1'-
-Jyz1·wz1'+
+Jy1·wy1'
После
подстановки
в полученные
выражения для
инерционных
моментов Мy1ин,
Mz2ин
кинематических
уравнений (4),
(4'), (5), (5') и преобразования,
получим следующий
вид выражений
для Мy1ин,
Mz2ин: MZ2ИН={cos(2·b)-2}·cos(a)2·tg(b)2·Jxy2(·wx02+wz02)+
+{2·tg(b)2·sin(b)2-2·cos(b)2+4}·sin(a)·cos(a)·Jxy2·wx0·wz0+
+{(Jy2-Jx2)/cos(b)-2·Jxy2·sin(b)(1+tg(b)2)}·cos(a)·wx0·wy2+
+Jyz2·wz0·wz2·(sin(a)-cos(a))/cos(b)-
-Jxz2·wx0'·cos(a)/cos(b)+
+{2·Jxy2·(sin(b)·tg(b)2+sin(b))·sin(a)+(Jx2-Jy2)·sin(a)/cos(b)}·wy2·wz0+
+Jxz2·wz0'·sin(a)/cos(b)+
+{Jxz2-Jyz2}·wy2·wz2·tg(b)+
+{(Jy2-Jx2)·tg(b)+Jxy2·(1-tg(b)2)}·wy22-
-{Jxz2·tg(b)+Jyz2}·wy2'+
+Jz2·wz2' (9) My1ин={[Jxz2·(tg(b)4+2/cos(b)2-1)·cos(b)3+Jyz1·tg(b)+Jxz1]·cos(a)2+
+[[(Jx1-Jz1)-Jxy1·tg(b)]·cos(a)-Jxz1·sin(a)]·sin(a)}·wx02+
+{[[Jxy1·tg(b)+(Jz1-Jx1)]·sin(a)-Jxz1·cos(a)]·cos(a)+
+[Jxz2·cos(b)3·[2/cos(b)2+tg(b)4-1]+Jyz1·tg(b)+Jxz1]·sin(a)2}·wz02+
+{(Jx1-Jz1)·cos(2·a)+[1-tg(b)4-2/cos(b)2]·Jxz2·cos(b)3·sin(2·a)-
-[Jyz1·tg(b)+2·Jxz1]·2·sin(a)·cos(a)-
-Jxy1·tg(b)·cos(2·a)}·wx0·wz0+
+{[Jxy2·sin(b)·cos(b)(tg(b)2+1)+(Jx2-Jz2)]·cos(a)}·wx0·wz2+
+{[Jxz2·sin(b)·cos(b)+Jxz2·sin(b)3/cos(b)+Jyz2]·cos(a)+
+[Jyz1·cos(a)-Jxy1·sin(a)]/cos(b)}·wx0·wy2-
-{[Jxz2·sin(b)·cos(b)·(1+tg(b)2)+Jyz2]·sin(a)+
+[Jyz1·sin(a)+Jxy1·cos(a)]/cos(b)}·wz0·wy2+
+{-[tg(b)2+1]·sin(b)·cos(b)·Jxy2+(Jz2-Jx2)]·sin(a)}·wz0·wz2+
+{[Jx2·sin(b)·cos(b)·(1+tg(b)2)+Jy1·tg(b)-(Jxy1+
+Jxy2)]·cos(a)-Jyz1·sin(a)}·wx0'+
+{[-Jx2·sin(b)·cos(b)·(1+tg(b)2)+(Jxy1+Jxy2)-
-Jy1·tg(b)]·sin(a)-Jyz1·cos(a)}·wz0'+
+{Jyz2·sin(b)-Jxz2·cos(b)]·wz22-
-{Jxz2·sin(b)+Jyz2·cos(b)}·wz2'+
+{(Jx2-Jy2)·sin(b)+Jxy2·cos(b)·(tg(b)2-1)}·wz2·wy2+
+{Jx2·sin(b)2/cos(b)-2·Jxy2·sin(b)+Jy2·cos(b)+Jy1/cos(b)}·wy2'
Численный
анализ инерционных
возмущающих
моментов (9) провожу
для различных
режимов работы
ГС, типовая
конструкция
которого приведена
на рис 2.
Рис.2.
Пусть
ГС имеет следующие
инерционные
параметры
наружной рамы
и платформы:
Jx1
= -------//------ Jx2=
2000 гсмс2
= 0.2 кгм2
Jy1
= 1500 гсмс2
= 0.15 кгм2 Jy2=
9500 гсмс2
= 0.95 кгм2
Jz1
= -------//------ Jz2
= 10000 гсмс2
= 1 кгм2
Jxy1
= Jyx1
= 0 Jxy2 =
Jyx2
= 0.0085 кгм2
Jxz1
= Jzx1
= 0 Jxz2
= Jzx2
= 0.023 кгм2
Jzy1
= Jyz1 =1500
гсмс2
= 0.15 кгм2 Jzy2
= Jyz2
= 0.04 кгм2
Угловые
скорости и
ускорения
основания и
управления
платформой
принимаю равными
их типовым
значениям при
работе гиростабилизатора
на кране.
wx0
= ±1
рад/с wy2
= ±2
рад/с
wy0
= ±1
рад/с wz2
= ±2
рад/с
wz0
= ±1
рад/с wy2'
= ±3
рад/с2 (10)
wx0'=
±0,2
рад/с2 wz2'
= ±3
рад/с2
wy0'=
±0,2
рад/с2
wz0'=
±0,2
рад/с2
Углы
прокачки рам
изменяются
в диапазоне:
a
= ±
2 рад. »
±
120 град. (10)
b
= ±1
рад. »
±
60 град.
Исследование
величины численных
значений инерционных
возмущающих
моментов провожу
с помощью программы
“MOMIN” листинг
которой приведен
в “Приложении
1”.
Анализ
инерционных
возмущающих
моментов провожу
для следующих
случаев работы
гиро-стабилизатора:
1) Работа
на неподвижном
основании при
наличии скоростей
управления
платформой;
2) Работа
на подвижном
основании при
неподвижной
платформе;
3) Работа
на подвижном
основании при
управляемой
платформе;
т.е. при
условии:
wx0
= wy0
= wz0
= wx0'
= wy0'
= wz0'
= 0 (11)
a
№
0; b
№
0; wy2№
0; wz2
№
0; wy2'
№
0; wz2'
№
0
Тогда
подставляя
(11) в выражения
для инерционных
моментов (9),
получаем следующий
их вид:
MZ2ИН=+{Jxz2-Jyz2}·wy2·wz2·tg(b)+
+{(Jy2-Jx2)·tg(b)+Jxy2·(1-tg(b)2)}·wy22-
-{Jxz2·tg(b)+Jyz2}·wy2'+
+Jz2·wz2'
MY1ИН=+{Jyz2·sin(b)-Jxz2·cos(b)}·wz22-
-{Jxz2·sin(b)+Jyz2·cos(b)}·wz2'+
+{(Jx2-Jy2)·sin(b)+Jxy2·cos(b)·(tg(b)2-1)}·wz2·wy2+
+{Jx2·sin(b)2/cos(b)-
-2·Jxy2·sin(b)+Jy2·cos(b)+Jy1/cos(b)}·wy2'
Максимальные
значения инерционных
моментов,
полученные
при выполнении
условий (10), следующие:
а) ось
Y1: Мy1ин
= Мин
+ Мцб
= 5.68 + 0.14 = 5.82 НЧм.
при a
= 0.067 рад.
b
= 1 рад.
wy2
= -2.0 рад/с.
wy2'
= 3.0 рад/с2.
wz2
= 2 рад/с.
wz2'
= -3.0 рад/с2.
где Мин
- вклад в Мy1ин
возмущающих
моментов, связаных
с осевыми моментами
инерции наружной
рамы и платформы;
Мцб
- вклад в Мy1ин
возмущающих
моментов, связаных
с центробежными
моментами
инерции наружной
рамы и платформы;
Вклад
Мцб
в суммарный
возмущающий
момент составил:
Мцб
Мин
+ Мцб
б) ось
Z2:
Мz2ин
= Мин
+ Мцб
= 7.67 + 0.33 = 8.0 НЧм.
при a
= 0.067 рад.
b
= 1 рад.
wy2
= 2.0 рад/с.
wy2'
= -3.0 рад/с2.
wz2
= -2 рад/с.
wz2'
= 3.0 рад/с2.
Вклад
Мцб
в суммарный
возмущающий
момент составил:
Мцб
Мин
+ Мцб
т.е. при:
wy2=
wy2'=
wz2
= wz2'
= 0; a
№
0; b
№
0; (12)
wx0
№
0; wy0
№
0; wz0
№
0; wx0'
№
0; wy0'
№
0; wz0'
№
0
Тогда
подставляя
(12) в выражения
для инерционных
моментов (9)
получаем следующий
их вид:
MZ2ИН={cos(2·b)-2}·cos(a)2·tg(b)2·Jxy2(·wx02+wz02)+
+{2·tg(b)2·sin(b)2-2·cos(b)2+4}·sin(a)·cos(a)·Jxy2·wx0·wz0+
-Jxz2·wx0'·cos(a)/cos(b)+
+Jxz2·wz0'·sin(a)/cos(b)+
MY1ИН={[Jxz2·(tg(b)4+2/cos(b)2-1)·cos(b)3+Jyz1·tg(b)+
+Jxz1]·cos(a)2+
+[[(Jx1-Jz1)-Jxy1·tg(b)]·cos(a)-Jxz1·sin(a)]·sin(a)}·wx02+
+{[[Jxy1·tg(b)+(Jz1-Jx1)]·sin(a)-Jxz1·cos(a)]·cos(a)+
+[Jxz2·cos(b)3·[2/cos(b)2+tg(b)4-1]+Jyz1·tg(b)+
+Jxz1]·sin(a)2}·wz02+
+{(Jx1-Jz1)·cos(2·a)+[1-tg(b)4-2/cos(b)2]·Jxz2·cos(b)3·
·sin(2·a)-[Jyz1·tg(b)+2·Jxz1]·2·sin(a)·cos(a)-
-Jxy1·tg(b)·cos(2·a)}·wx0·wz0+
+{[Jx2·sin(b)·cos(b)·(1+tg(b)2)+Jy1·tg(b)-(Jxy1+Jxy2)]·cos(a)-
-Jyz1·sin(a)}·wx0'+
+{[-Jx2·sin(b)·cos(b)·(1+tg(b)2)+(Jxy1+Jxy2)-Jy1·tg(b)]·sin(a)-
-Jyz1·cos(a)}·wz0'+
При этом
получены следующие
максимальные
значения инерционных
возмущающих
моментов:
а) ось
Y1:
Мy1ин
= Мин
+ Мцб
= 0.154 + 0.551= 0.705 НЧм.
при a
= - 0.82 рад.
b
= 1 рад.
wx0
= wz0
= 1 рад/с.
wx0'
= wz0'
= 0.2 рад/с2.
wy0
= 0.167 рад/c.
wy0'
= 0.167 рад/с2.
Вклад
Мцб
в суммарный
возмущающий
момент составил:
Мцб
Мин
+ Мцб
б) ось
Z2:
Мz2ин
= Мин
+ Мцб
= 0 + 0.07= 0.07 НЧм.
при a
= - 0.785 рад.
b
= 1 рад.
wx0
= wz0
= 1 рад/с.
wx0'
= wz0'
= 0.2 рад/с2.
wy0
= 0.167 рад/с.
wy0'
= 0.167 рад/c2
Вклад
Мцб
в суммарный
возмущающий
момент составил:
Мцб
Мин
+ Мцб
При
подвижном
основании и
управляемой
платформе
инерционные
возмущающие
моменты определяются
выражениями
(9).
MZ2ИН={cos(2·b)-2}·cos(a)2·tg(b)2·Jxy2(·wx02+wz02)+
+{2·tg(b)2·sin(b)2-2·cos(b)2+4}·sin(a)·cos(a)·Jxy2·wx0·wz0+
+{(Jy2-Jx2)/cos(b)-2·Jxy2·sin(b)(1+tg(b)2)}·cos(a)·wx0·wy2+
+Jyz2·wz0·wz2·(sin(a)-cos(a))/cos(b)-
-Jxz2·wx0'·cos(a)/cos(b)+
+{2·Jxy2·(sin(b)·tg(b)2+sin(b))·sin(a)+(Jx2-
-Jy2)·sin(a)/cos(b)}·wy2·wz0+
+Jxz2·wz0'·sin(a)/cos(b)+
+{Jxz2-Jyz2}·wy2·wz2·tg(b)+
+{(Jy2-Jx2)·tg(b)+Jxy2·(1-tg(b)2)}·wy22-
-{Jxz2·tg(b)+Jyz2}·wy2'+
+Jz2·wz2'
MY1ИН={[Jxz2·(tg(b)4+2/cos(b)2-1)·cos(b)3+Jyz1·tg(b)+
+Jxz1]·cos(a)2+
+[[(Jx1-Jz1)-Jxy1·tg(b)]·cos(a)-Jxz1·sin(a)]·sin(a)}·wx02+
+{[[Jxy1·tg(b)+(Jz1-Jx1)]·sin(a)-Jxz1·cos(a)]·cos(a)+
+[Jxz2·cos(b)3·[2/cos(b)2+tg(b)4-1]+
+Jyz1·tg(b)+Jxz1]·sin(a)2}·wz02+
+{(Jx1-Jz1)·cos(2·a)+[1-tg(b)4-2/cos(b)2]·Jxz2·cos(b)3·
·sin(2·a)-[Jyz1·tg(b)+2·Jxz1]·2·sin(a)·cos(a)-
-Jxy1·tg(b)·cos(2·a)}·wx0·wz0+
+{[Jxy2·sin(b)·cos(b)(tg(b)2+1)+(Jx2-Jz2)]·cos(a)}·wx0·wz2+
+{[Jxz2·sin(b)·cos(b)+Jxz2·sin(b)3/cos(b)+Jyz2]·cos(a)+
+[Jyz1·cos(a)-Jxy1·sin(a)]/cos(b)}·wx0·wy2-
-{[Jxz2·sin(b)·cos(b)·(1+tg(b)2)+Jyz2]·sin(a)+
+[Jyz1·sin(a)+Jxy1·cos(a)]/cos(b)}·wz0·wy2+
+{-[tg(b)2+1]·sin(b)·cos(b)·Jxy2+(Jz2-Jx2)]·sin(a)}·wz0·wz2+
+{[Jx2·sin(b)·cos(b)·(1+tg(b)2)+Jy1·tg(b)-(Jxy1+Jxy2)]
·
·cos(a)-Jyz1·sin(a)}·wx0'+
+{[-Jx2·sin(b)·cos(b)·(1+tg(b)2)+(Jxy1+Jxy2)-Jy1·tg(b)]·
·sin(a)-Jyz1·cos(a)}·wz0'+
+{Jyz2·sin(b)-Jxz2·cos(b)}·wz22-
-{Jxz2·sin(b)+Jyz2·cos(b)}·wz2'+
+{(Jx2-Jy2)·sin(b)+Jxy2·cos(b)·(tg(b)2-1)}·wz2·wy2+
+{Jx2·sin(b)2/cos(b)-2·Jxy2·sin(b)+Jy2·cos(b)+
+Jy1/cos(b)}·wy2'
При этом
получены следующие
максимальные
значения инерционных
моментов.
а) ось
Y1:
Мy1ин
= Мин
+ Мцб
= 8.1 + 1.65 = 9.75 НЧм
при a
= 0.776 рад.
b
= 1.0 рад.
wy2
= -2 рад/с.
wy2'
= 3 рад/с2.
wz2
= 2 рад/с.
wz2'
= -3 рад/с2.
wx0
= wz0
= 1 рад/c.
wx0'
= 0.2 рад/c2.
wz0'
= - 0.2 рад/c2.
wy0
= 0.167 рад/c.
wy0'
= 0.167 рад/c2.
Вклад
Мцб
в суммарный
возмущающий
момент составил:
Мцб
Мy1ин+Мцб
б) ось
Z2:
Мz2ин
= Мин
+ Мцб
= 11.6 + 0.361 = 11.96 НЧм
при a
= -0.785 рад.
b
= 1.0 рад.
wy2
= 2 рад/с.
wy2'
= -3 рад/с2.
wz2
= -2 рад/с.
wz2'
= 3 рад/с2.
wx0
=
wz0
= 1 рад/c.
wx0'
= wz0'
= - 0.2 рад/c2.
wy0
= 0.167 рад/c.
wy0'
= 0.167 рад/c2.
Вклад
Мцб
в суммарный
возмущающий
момент составил:
Мцб
Мy1ин
+ Мцб
Анализ устойчивости
ГС с нежесткими
наружной рамой,
креплением
статора двигателя
стабилизации
к раме, с нежесткими
редуктором
и связью платформы
с объектом
стабилизации,
проводим основываясь
на следующей
физической
модели:
Рис. 1.
здесь Ji
- момент инерции
i-го элемента;
Ci,j - коэффициент
упругости;
Di,j - коэфф.
демпфирования
между i и j
элементами;
K - коэффициент
передачи цепи
обратной
связи.
Оценку влияния
каждого из
входящих в
модель элементов
(Ji,Ci,j,Di,j) выполняем
на основе анализа
поведения ЛАХ
разомкнутой
системы, при
вариациях
Ji,Ci,j,Di,j.
Уравнения
движения каждого
из элементов
модели в общем
виде могут быть
представлены
следующим
образом:
JiЧxi''+Di-1,iЧ(xi'-xi-1')-Di,i+1Ч(xi+1'-xi')+Ci-1,iЧ
(xi-xi-1)-Ci,i+1Ч(xi+1
- xi) = Мi (1)
где Мi -
внешний момент
действующий
на i-й элемент;
xi,xi', xi''- перемещение,
скорость и
ускорение i-го
элемента.
Расписав
уравнение (1)
для каждого
элемента, получим
следующюю
систему уравнений
движения модели:
J1Чx1''+D01Ч(x1'-x0')-D12Ч(x2'-x1')+C01Ч(x1-x0)-C12Ч(x2-x1)=
М1
J2Чx2''+D12Ч(x2'-x1')-D23Ч(x3'-x2')+C12Ч(x2-x1)-C23Ч(x3-x2)=
М2
J3Чx3''+D23Ч(x3'-x2')-D34Ч(x4'-x3')+C23Ч(x3-x2)-C34Ч(x4-x3)=
М3 (2)
J4Чx4''+D34Ч(x4'-x3')-D45Ч(x5'-x4')+C34Ч(x4-x3)-C45Ч
(x5-x4)= М4
J5Чx5''+D45Ч(x5'-x4')-D56Ч(x6'-x5')+C45Ч(x5-x4)-C56Ч(x6-x5)=
М5
Раскрыв в
(2) скобки и преобразовав
получаем следующий
вид уравнений
движения модели.
-D01Чx0'-C01Чx0+J1Чx1''+(D01+D12)Чx1'+(C01+C12)Чx1-D34Чx2'-
-C12Чx2=
М1
-D12Чx1'-C12Чx1+J2Чx2''+(D23+D12)Чx2'+(C12+C23)Чx2-D23Чx3'-
-C23Чx3=
М2
-D23Чx2'-C23Чx2+J3Чx3''+(D23+D34)Чx3'+(C23+C34)Чx3-D34Чx4'-
-C34Чx4=М3
(3)
-D34Чx3'-C34Чx3+J4Чx4''+(D34+D45)Чx4'+(C34+C45)Чx4-D45Чx5'-
-C45Чx5=
М4
-D45Чx4'-C45Чx4+J5Чx5''+(D45+D56)Чx5'+(C45+C56)Чx5-D56Чx6'-
-C56Чx6=
М5
Переписав
(3) в операторной
форме получаем
уравнения
движения модели
в следующем
виде.
-(D01Чs+C01)Чx0+(J1Чs2+(D01+D12)Чs+(C01+C12))Чx1
-
-(D12Чs+C12)Чx2=
М1
-(D12Чs+C12)Чx1+(J2Чs2+(D12+D23)Чs+(C12+C23))Чx2-
-(D23Чs+C23)Чx3=
КЧx4
-(D23Чs+C23)Чx2+(J3Чs2+(D23+D34)Чs+(C23+C34))Чx3-
-(D34Чs+C34)Чx4=-КЧx4
-(D34Чs+C34)Чx3+(J4Чs2+(D34+D45)Чs+(C34+C45))Чx4-
-(D45Чs+C45)Чx5=
М4 (4)
-(D45Чs+C45)Чx4+(J5Чs2+(D45+D56)Чs+(C45+C56))Чx5-
-(D56Чs+C56)Чx6=
М5
Для нахождения
передаточной
функции разомкнутой
системы по
управляющему
воздействию
Wp(s) составим два
определителя:
главный - D,
и характеризующий
входное воздействие
D1,
с учетом того,
что x0=0; D56=0; C56=0; C23=0.
a11 a12 0 0 0
a21 a22 a23 0 0
D= 0 a32 a33 a34 0 (5)
0 0 a43 a44 a45
0 0 0 a54 a55
где a11 =
J1Чs2+(D01+D12)Чs+C01+C12
a12 = -D12Чs-C12
a21 = a12
a22 = J2Чs2+(D12+D23)Чs+C12
a23 = -D23Чs
a32 = a23
a33 = J3Чs2+(D23+D34)Чs+C34
a34 = -D34Чs-C34
a43 = a34
a44 = J4Чs2+(D34+D45)Чs+C34+C45
a45 = -D45Чs-C45
a54 = a45
a55 = J5Чs2+D45Чs+C45
a21 a22 a23 -KЧx4 0
D1= 0 a32 a33
KЧx4 0 (6)
0 0 a43 0 a45
0 0 0 0 a55
Передаточная
функция разомкнутой
системы определяется
как:
Коэффициенты
ai, bi полиномов
числителя и
знаменателя
передаточной
функции Wp(s) выражаются
через параметры
элементов
модели следующим
образом:
(8)
a9=J1J2J3J4J5
a8=D01J2J3J4J5+D12J3J4J5(J1+J2)+J1(D23J4J5(J2+J3)+J2(D34J5(J3+J4)+D45J3(J4+J5)))
a7=C01J2J3J4J5+C12J3J4J5(J1+J2)+C34J1J2(J3J5+J4J5)+C45J1J2J3(J4+J5)+D01(D12J3J4J5+D23J4J5(J2+J3)+J2(D34J5(J3+J4)+D45J3(J4+J5)))+D12(D23J4J5(J1+J2+J3)+(J1+J2)(D34J5(J3+J4)+D45J3(J4+J5)))+J1(D23(D34J5(J2+J3+J4)+D45(J4+J5)(J2+J3))+D34D45J2(J3+J4+J5))
a6=C01(D12J3J4J5+D23J4J5(J2+J3)+J2(D34J5(J3+J4)+D45J3(J4+J5)))+C12(D01J3J4J5+D23J4J5(J1+J2+J3)+(J1+J2)(D34J5(J3+J4)+D45J3(J4+J5)))+C34(D01J2(J3J5+J4J5)+D12J5(J3+J4)(J1+J2)+J1(D23J5(J2+J3+J4)+D45J2(J3+J4+J5)))+C45(D01J2J3(J4+J5)+D12J3(J4+J5)(J1+J2)+J1(D23(J4+J5)(J2+J3)+D34J2(J3+J4+J5)))+D01(D12(D23J4J5+D34(J3J5+J4J5)+D45J3(J4+J5))+D23(D34(J2J5+J3J5+J4J5)+D45(J4+J5)(J2+J3))+D34D45J2(J3+J4+J5))+D12(D23(D34(J1J5+J2J5+J3J5+J4J5)+D45(J4+J5)(J1+J2+J3))+D34D45(J3+J4+J5)(J1+J2))+D23D34D45J1(J2+J3+J4+J5)
a5=C01(C12J3J4J5+C34J2(J3J5+J4J5)+C45J2J3(J4+J5)+D12(D23J4J5+D34(J3J5+J4J5)+D45J3(J4+J5))+D23(D34(J2J5+J3J5+J4J5)+D45(J4+J5)(J2+J3))+D34D45J2(J3+J4+J5))+C12(C34J5(J3+J4)(J1+J2)+C45J3(J4+J5)(J1+J2)+D01(D23J4J5+D34(J3J5+J4J5)+D45J3(J4+J5))+D23(D34(J1J5+J2J5+J3J5+J4J5)+D45(J4+J5)(J1+J2+J3))+D34D45(J3+J4+J5)(J1+J2))+C34(C45J1J2(J3+J4+J5)+D01(D12(J3J5+J4J5)+D23(J2J5+J3J5+J4J5)+D45J2(J3+J4+J5))+D12(D23(J1J5+J2J5+J3J5+J4J5)+D45(J3+J4+J5)(J1+J2))+D23D45J1(J2+J3+J4+J5))+C45(D01(D12J3(J4+J5)+D23(J2(J4+J5)+J3(J4+J5))+D34J2(J3+J4+J5))+D12(D23(J1(J4+J5)+J2(J4+J5)+J3(J4+J5))+D34(J3+J4+J5)(J1+J2))+D23D34J1(J2+J3+J4+J5))+D01(D12(D23(D34J5+D45(J4+J5))+D34D45(J3+J4+J5))+D23D34D45(J2+J3+J4+J5))+D12D23D34D45(J1+J2+J3+J4+J5)
a4=C01(C12(D23J4J5+D34(J3J5+J4J5)+D45J3(J4+J5))+C34(D12(J3J5+J4J5)+D23(J2J5+J3J5+J4J5)+D45J2(J3+J4+J5))+C45(D12J3(J4+J5)+D23(J2(J4+J5)+J3(J4+J5))+D34J2(J3+J4+J5))+D12(D23(D34J5+D45(J4+J5))+D34D45(J3+J4+J5))+D23D34D45(J2+J3+J4+J5))+C12(C34(D01(J3J5+J4J5)+D23(J1J5+J2J5+J3J5+J4J5)+D45(J3+J4+J5)(J1+J2))+C45(D01J3(J4+J5)+D23(J1(J4+J5)+J2(J4+J5)+J3(J4+J5))+D34(J3+J4+J5)(J1+J2))+D01(D23(D34J5+D45(J4+J5))+D34D45(J3+J4+J5))+D23D34D45(J1+J2+J3+J4+J5))+C34(C45(D01J2(J3+J4+J5)+D12(J1(J3+J4+J5)+J2(J3+J4+J5))+D23J1(J2+J3+J4+J5))+D01(D12(D23J5+D45(J3+J4+J5))+D23D45(J2+J3+J4+J5))+D12D23D45(J1+J2+J3+J4+J5))+C45(D01(D12(D23(J4+J5)+D34(J3+J4+J5))+D23D34(J2+J3+J4+J5))+D12D23D34(J1+J2+J3+J4+J5))+D01D12D23D34D45
a3=C01(C12(C34(J3J5+J4J5)+C45J3(J4+J5)+D23(D34J5+D45(J4+J5))+D34D45(J3+J4+J5))+C34(C45J2(J3+J4+J5)+D12(D23J5+D45(J3+J4+J5))+D23D45(J2+J3+J4+J5))+C45(D12(D23(J4+J5)+D34(J3+J4+J5))+D23D34(J2+J3+J4+J5))+D12D23D34D45)+C12(C34(C45(J1(J3+J4+J5)+J2(J3+J4+J5))+D01(D23J5+D45(J3+J4+J5))+D23D45(J1+J2+J3+J4+J5))+C45(D01(D23(J4+J5)+D34(J3+J4+J5))+D23D34(J1+J2+J3+J4+J5))+D01D23D34D45)+C34(C45(D01(D12(J3+J4+J5)+D23(J2+J3+J4+J5))+D12D23(J1+J2+J3+J4+J5))+D01D12D23D45)+C45D01D12D23D34
a2=C01(C12(C34(D23J5+D45(J3+J4+J5))+C45(D23(J4+J5)+D34(J3+J4+J5))+D23D34D45)+C34(C45(D12(J3+J4+J5)+D23(J2+J3+J4+J5))+D12D23D45)+C45D12D23D34)+C12(C34(C45(D01(J3+J4+J5)+D23(J1+J2+J3+J4+J5))+D01D23D45)+C45D01D23D34)+C34C45D01D12D23
a1=C01(C12(C34(C45(J3+J4+J5)+D23D45)+C45D23D34)+C34C45D12D23)+C12C34C45D01D23
a0=C01C12C34C45D23
b7=D34J1J2J5
b6=(C34J1J2J5+D34(D01J2J5+D12J5(J1+J2)+D45J1J2))
b5=(C01D34J2J5+C12D34J5(J1+J2)+C34(D01J2J5+D12J5(J1+J2)+D45J1J2)+C45D34J1J2+D34(D01(D12J5+D45J2)+D12D45(J1+J2)))
b4=(C01(C34J2J5+D12D34J5+D34D45J2)+C12(C34J5(J1+J2)+D01D34J5+D34D45(J1+J2))+C34(C45J1J2+D01(D12J5+D45J2)+D12D45(J1+J2))+C45D34(D01J2+D12(J1+J2))+D01D12D34D45)
b3=(C01(C12D34J5+C34(D12J5+D45J2)+C45D34J2+D12D34D45)+C12(C34(D01J5+D45(J1+J2))+C45D34(J1+J2)+D01D34D45)+C34(C45(D01J2+D12(J1+J2))+D01D12D45)+C45D01D12D34)
b2=(C01(C12(C34J5+D34D45)+C34(C45J2+D12D45)+C45D12D34)+C12(C34(C45(J1+J2)+D01D45)+C45D01D34)+C34C45D01D12)
b1=(C01(C12(C34D45+C45D34)+C34C45D12)+C12C34C45D01)
b0=C01C12C34C45
Представить
передаточную
функцию Wp(s) в виде
произведения
полиномов не
выше второго
порядка в числителе
и знаменателе
Wp(s) в аналитическом
виде не представляется
возможным даже
теоретически,
т.к. вид корней
характеристических
полиномов
ai,bi, а, следовательно,
и вид разложения
на полиномы
не выше второго
порядка, зависит
от численных
значений параметров
элементов
модели. Поэтому
исследование
влияния элементов
модели на
устойчивость
ГС проводилось
численно, путем
нахождения
корней характеристических
полиномов для
каждого частного
случая. Далее
по полученным
корням определялись
полиномы не
выше второго
порядка по
которым и строилась
ЛАХ разомкнутой
системы.
Все математические
операции проводилось
с использованием
пакета “MATHCAD”
с помощью которого
численно определялись
корни полиномов
в передаточной
функции разомкнутой
системы Wp(s), зная
которые можно
представить
Wp(s) в виде последовательного
соединения
элементарных
звеньев. Это
выполняется
следующим
образом. Пусть
полиномы числителя
и знаменателя
Wp(s) имеют корни
lai,
lbi
соответственно.
Эти корни могут
быть нулевыми,
действительными
и комплексно
сопряженными.
Каждый нулевой
корень знаменателя
lai=0
обеспечивает
появление в
составе Wp(s)
интегрирующего
звена Wi(s)= 1/s, соответственно
lbi=0
отвечает за
появление чисто
дифференцирующего
звена с Wi(s)= s. Каждый
из действительных
корней lai,
lbi
приносит в
числитель или
знаменатель
соответственно
выражение вида
(TiЧs+1)Ч(1/Ti),
где Ti=1/li
, что соответствует
появлению
апериодических
и дифференцирующих
звеньев в составе
Wp(s). Каждая пара
комплексно
сопряженных
корней li,
li*
в составе числителя
или знаменателя
передаточной
функции отвечает
за появление
в числителе
или знаменателе
соответственно
выражений вида
(Ti2 Ч
s2 +2ЧxiЧTiЧs
+1)Ч(1/Ti2),
где Ti=1 / |li|
, xi=Re(li)
/ |li|.
Таким образом,
зная корни
полиномов
числителя и
знаменателя
передаточной
функции можно
представить
её в виде:
П(si)ЧП(TgЧs+1)ЧП(
Tn2 Ч
s2 +2ЧxnЧTnЧs
+1)
П(sj)ЧП(TkЧs+1)ЧП(
Tm2 Ч
s2 +2ЧxmЧTmЧs
+1)
П(1/Ti) Ч
П(1/Ti2)
П(1/Ti) Ч
П(1/Ti2)
Для численных
расчетов примем
базовые параметры
модели характерными
для ГС данного
типа, которые
равны следующим
значениям:
J1 = 0.25 кгЧм2 C01
= 1Ч103
НЧм/рад. D01=0.001
НЧмЧс
J2 = 0.03 кгЧм2 C12
= 1Ч103
НЧм/рад. D12=0.001
НЧмЧс
J3 = 0.01 кгЧм2 C23
= 0 D23=0.1 НЧмЧс
J4 = 0.15 кгЧм2 C34
=1Ч104
НЧм/рад. D34=0.001
НЧмЧс
J5 = 1 кгЧм2 C45
=1Ч103
НЧм/рад. D45=0.01
НЧмЧс
К = 1000
Рассмотрим
следующие
варианты модели:
1) ГС с
“жесткими”
рамами и редуктором.
Начальные
параметры
модели принимают
следующие
знечения:
J1 = 0.25 кгЧм2 C01
= 1Ч1020
НЧм/рад. D01=
0.001 НЧмЧс
J2 = 0.03 кгЧм2 C12
= 1Ч1020
НЧм/рад. D12=
0.001 НЧмЧс
J3 = 0.01 кгЧм2 C23
= 0 D23=0.1 НЧмЧс
J4 = 0.15 кгЧм2 C34
=1Ч1020
НЧм/рад. D34=0.001
НЧмЧс
J5 = 1 кгЧм2 C45
=1Ч1020
НЧм/рад. D45=0.01
НЧмЧс
К = 1000
Варьируем
D23 = 0.01 ... 1 HЧмЧс
Передаточная
функция при
этом имеет вид:
k Ч
kw
s Ч
(TЧs+1)
Значения
постоянной
времени Т, w,
kw приведены в
Табл.1.
Табл.1. w=1/T
Т.о. ЛАХ модели
с бесконечно
жесткими пружинами
соответствует
ЛАХ идеализированного
индикаторного
ГС. Постоянная
времени Т
апериодического
звена апроксимируется
формулой:
J3 +J4 +J5
D23
2) ГС с
“нежестким”
редуктором.
Начальные
параметры
модели:
J1 = 0.25 кгЧм2 C01
= 1Ч1020
НЧм/рад. D01=
0.001 НЧмЧс
J2 = 0.03 кгЧм2 C12
= 1Ч1020
НЧм/рад. D12=
0.001 НЧмЧс
J3 = 0.01 кгЧм2 C23
= 0 D23=0.1 НЧмЧс
J4 = 0.15 кгЧм2 C34
=1Ч104
НЧм/рад. D34=0.001
НЧмЧс
J5 = 1 кгЧм2 C45
=1Ч1020
НЧм/рад. D45=0.01
НЧмЧс
К = 1000
Варьируем
нежесткость
редуктора
С34=103 ... 107 HЧм/рад.
Передаточная
функция при
этом имеет вид:
k Ч
kw
s Ч
(T1Чs+1)Ч(
T22 Ч
s2 +2Чx2ЧT2Чs
+1)
Значения
постоянных
времени Т1,
Т2, соответствующие
им частоты
“излома” ЛАХ
w1,
w1,
удельный коэффициент
демпфирования
x2
и коэффициент
передачи модели
kw приведены в
Табл.2. и Табл.3.
Табл.2. w1 w2 x2 3.1Ч10-4 1Ч10-4
Как видно
из Табл.2. нежесткость
редуктора
влияет только
на расположение
колебательного
звена на оси
частот (Т2,
w2)
и коэффициент
демпфирования
в этом звене
(x2).
Влияние
демпфирования
в редукторе
на поведение
ЛАХ определяем
варьируя D34=0.001
... 0.1 НЧмЧс
(при С34=104 = const.).
Табл.3. w1 w2 x2
Как видно
из Табл.3., изменение
демпфирования
в редукторе
влияет не только
на коэффициент
демпфирования
в колебательном
звене, но и на
расположение
на оси частот
апериодического
звена (Т1), и
на коэффициент
передачи модели.
3) ГС с
“нежесткой”
связью платформы
со стабилизируемым
объектом
(телекамерой).
Исходные
параметры
модели:
J1 = 0.25 кгЧм2 C01
= 1Ч1020
НЧм/рад. D01=
0.001 НЧмЧс
J2 = 0.03 кгЧм2 C12
= 1Ч1020
НЧм/рад. D12=
0.001 НЧмЧс
J3 = 0.01 кгЧм2 C23
= 0 D23=0.1 НЧмЧс
J4 = 0.15 кгЧм2 C34
=1Ч1020
НЧм/рад. D34=0.001
НЧмЧс
J5 = 1 кгЧм2 C45
=1Ч103
НЧм/рад. D45=0.01
НЧмЧс
К = 1000
Варьируем
С45 = 102 ... 106 HЧм/рад.
Передаточная
функция при
этом имеет вид:
k Ч
kwЧ(
T32 Ч
s2 +2Чx3ЧT3Чs
+1)
s Ч
(T1Чs+1)Ч(
T22 Ч
s2 +2Чx2ЧT2Чs
+1)
Влияние
жесткости
крепления
стабилизируемого
объекта к платформе
на передаточную
функцию Wp(s) приведено
в Табл.4.
Табл.4. T1
(w1) T2
(w2) x2 T3
(w3) x3 5Ч10-4 1.6Ч10-4 5Ч10-5 1.2Ч10-3(850) 3.2Ч10-3(313) 1.6Ч10-5
Влияние
демпфирования
в креплении
стабилизируемого
объекта к платформе
на передаточную
функцию Wp(s) приведено
в Табл.5. Коэффициент
демпфирования
изменяется
в пределах
D45=0.001 ... 0.1 НЧмЧс,
при постоянной
жесткости
крепления
объекта к платформе
равной C45=1000 HЧм/рад
=const.
Табл.5. T1
(w1) T2
(w2) x2 T3
(w3) x3 2.7Ч10-14 1.6Ч10-4 1.6Ч10-3
Как видно
из Табл.4. и 5.,
нежесткость
крепления
объекта к платформе
вызывает появление
в составе ЛАХ
двух звеньев:
колебательного
и антиколебательного,
причем антиколебательное
звено всегда
расположено
в области более
низких частот,
чем колебательное.
Это влечет
появление в
ЛАХ участка
с наклоном в
0 Дб/дек., который
в случае его
расположения
до частоты
среза, увеличивает
частоту среза,
что вызывает
трудности в
технической
реализации
такой системы
стабилизации.
Демпфирование
в креплении
объекта к платформе
влияет только
на удельные
коэффициенты
демпфирования
x2,
x3
в колебательном
и антиколебательном
звеньях, причем
особенно сильно
изменяется
x3.
4) ГС с
“нежестким”
креплением
статора двигателя
стабилизации
к наружной раме
(задняя нежесткость).
Параметры
модели:
J1 = 0.25 кгЧм2 C01
= 1Ч1020
НЧм/рад. D01=
0.001 НЧмЧс
J2 = 0.03 кгЧм2 C12
= 1Ч103
НЧм/рад. D12=
0.001 НЧмЧс
J3 = 0.01 кгЧм2 C23
= 0 D23=0.1 НЧмЧс
J4 = 0.15 кгЧм2 C34
=1Ч1020
НЧм/рад. D34=0.001
НЧмЧс
J5 = 1 кгЧм2 C45
=1Ч1020
НЧм/рад. D45=0.01
НЧмЧс
К = 1000
Варьируем
С12 = 102 ... 106 HЧм/рад.
Передаточная
функция при
этом имеет вид:
k Ч
kwЧ(
T32 Ч
s2 +2Чx3ЧT3Чs
+1)
s Ч
(T1Чs+1)Ч(
T22 Ч
s2 +2Чx2ЧT2Чs
+1)
Варьируем
С12 (при D12=0.001
НЧмЧс=const),
результаты
приведены в
Табл.6.
Табл.6. x2 x3 9.1Ч10-5 2.9Ч10-5 9.1Ч10-6
Варьируем
D12 (при С12=1000 HЧм/рад
= const.), результаты
приведены в
Табл.7.
Табл.7. T1/w1 T2
/ w2 x2 T3
/ w3 x3 8.3Ч10-14 9.1Ч10-5
Как видно
из Табл.6, нежесткость
крепления
статора двигателя
стабилизации
к основанию,
приводит к
появлению в
составе передаточной
функции Wp(s)
колебательного
и антиколебательного
звеньев с одинаковыми
постоянными
времени и различными
коэффициентами
демпфирования.
Т.к. постоянные
времени этих
звеньев одинаковы,
то наличие
“задней” нежесткости
никак не отражается
на виде ЛАХ,
однако различия
этих звеньев
в коэффициентах
демпфирования
влекут разную
скорость
“переключения”
фазы в каждом
звене, что вызывает
появление
незначительных
по амплитуде
выбросов на
фазо-частотной
характеристике.
5) ГС с
“нежесткой”
наружной рамкой.
Исходные
параметры
модели:
J1 = 0.25 кгЧм2 C01
= 1Ч103
НЧм/рад. D01=
0.001 НЧмЧс
J2 = 0.03 кгЧм2 C12
= 1Ч103
НЧм/рад. D12=
0.001 НЧмЧс
J3 = 0.01 кгЧм2 C23
= 0 D23=0.1 НЧмЧс
J4 = 0.15 кгЧм2 C34
=1Ч1020
НЧм/рад. D34=0.001
НЧмЧс
J5 = 1 кгЧм2 C45
=1Ч1020
НЧм/рад. D45=0.01
НЧмЧс
К = 1000
Варьируем
С01,С12 = 102 ... 106 HЧм/рад.
Передаточная
функция при
этом имеет вид:
k Ч
kwЧ(
T42 Ч
s2 +2Чx4ЧT4Чs
+1) Ч(
T52 Ч
s2 +2Чx5ЧT5Чs
+1)
s Ч
(T1Чs+1)Ч(
T22 Ч
s2 +2Чx2ЧT2Чs
+1) Ч(
T32 Ч
s2 +2Чx3ЧT3Чs
+1)
Вначале
варьируем С01,
при С12=const., результаты
приведены в
Табл.8.
Табл.8. x2 x3 x4 x5
1.8 10-13
9.7 10-5
1.76 10-13
5.7 10-5
4.8 10-5
9.4 10-6
8.9 10-5
4.37 10-13
Далее варьируем
С12, при С01=const.,
результаты
- в Табл.9.
Табл.9. x2 x3 x4 x5
1.5 10-4
1.8 10-4
3.1 10-5
1.8 10-13
0.9 10-5
1.8 10-13
Варьируя
последовательно
D01 и D12 выявляем
степень их
влияния на Ti,
при С01=С12=1000 HЧм/рад
= const. (Табл.10,11)
Табл.10. T1/w1 T2
/ w2 x2 T3
/ w3 x3 T4
/ w4 x4 T5
/ w5 x5 9.6Ч10-5
2Ч 10-14 9.7Ч10-5
17Ч 10-14 11Ч10-5 23Ч10-5
Табл.11. T1/w1 T2
/ w2 x2 T3
/ w3 x3 T4
/ w4
|