CC BY
9
https://doi.org/10.38013/2542-0542-2022-2-30-46 УДК 629.735.33.015.075
Повышение адекватности моделирования маневра летательного аппарата путем определения положений мгновенных осей его вращательного движения
А. А. Писковацкий
Публичное акционерное общество «Долгопрудненское научно-производственное предприятие», Долгопрудный, Московская область, Российская Федерация
Возможности разработки автоматического управления маневрами с критическими углами атаки и скольжения и автоматического вывода летательного аппарата (ЛА) из аварийных режимов полета ограничены наличием динамических погрешностей системы управления, превышение критического уровня которых приводит к срыву маневра, а при попадании в область неустойчивого полета чревато потерей ЛА. Одной из причин появления таких погрешностей являются методические погрешности, возникающие из-за того, что синтез законов управления (ЗУ) происходит на основе недостаточно точных математических моделей (ММ) динамики движения данных ЛА. Наличие методических погрешностей управления, определяемых принципом математического моделирования полета, не позволяет определять соответствующие данным процессам полета ЗУ и использовать их для выполнения таких маневров в автоматических режимах полета.
Цель исследования - повысить адекватность ММ маневрирования ЛА путем разработки новых принципов ее построения с использованием критерия с более жесткими требованиями к определению эквивалентности смоделированного и реального маневра ЛА (далее СМ и РМ ЛА).
В статье доказательством трех теорем сформированы условия принятия эквивалентными СМ и РМ ЛА и обоснована возможность снижения динамических ошибок системы управления в быстропротекающих процессах маневрирования ЛА путем использования нового принципа математического моделирования такого полета.
Несмотря на теоретико-исследовательскую направленность статьи, ожидаемые результаты от повышения адекватности ММ вследствие более точного прогнозирования параметров маневра ЛА и формирования сигналов управления данным маневрированием без методических погрешностей имеют большое знаем чение для практического применения.
0
pj Ключевые слова: адекватность модели, эквивалентность состояний, совпадение траекторий, мгновенная ц ось, распределение ускорения
J Для цитирования: Писковацкий А. А. Повышение адекватности моделирования маневра летательного £ аппарата путем определения положений мгновенных осей его вращательного движения // Вестник Кон-it церна ВКО «Алмаз - Антей». 2022. № 2. С. 30-46. https://doi.org/10.38013/2542-0542-2022-2-30-46
I
« For citation: Piskovatsky A. A. Improving the adequacy of aircraft manoeuvre imulation by determining the positions of the instantaneous axes of its rotational motion // Vestnik Koncerna VKO "Almaz - Antey". 2022. < No. 2. P. 30-46. https://doi.org/10.38013/2542-0542-2022-2-30-46
° Поступила 25.11.2020 Отрецензирована 26.02.2021 Одобрена 21.02.2022 Опубликована 17.05.2022
CQ та
1
о.
<D J I
0 ^
s
1 I-
o
<D CQ
CM ■Clio 9
CM ■Clio
CM
5
W
I! © Писковацкий А. А., 2022
Определения и сокращения
АТТ - абсолютно твердое тело вращательное движение ЛА - вращательное движение ЛА, рассматриваемое изолированно от его поступательного движения ДВС - действие внешних сил на ЛА, оказывающее влияние на динамику его движения действующее распределение сил (влияющее распределение массы) ЛА по координатам его ССК - распределение сил (массы) ЛА по координатам его ССК после линеаризации его значений
ЛА - летательный аппарат ММ - математическая модель маневрирования ЛА
МОДС - мгновенная ось действия момента внешних сил, влияющих на изменение параметров вращательного движения ЛА МОУС - мгновенная ось угловой скорости, определяемая как множество точек пространства с нулевыми значениями скорости, определенных для вращательного движения ЛА
МОУУ - мгновенная ось углового ускорения, определяемая как множество точек пространства с нулевыми значениями линейного ускорения ЛА, определенного для изолированного вращательного движения ЛА НС - неравномерно распределенные силы, формирующие угловое ускорение ЛА ПЛСК 1 типа - подвижная локальная система координат, связанная с МОУС ЛА ПЛСК 2 типа - подвижная локальная система координат, связанная с положением ССК ЛА в начальный момент шага моделирования равнодействующая распределенной силы -вектор, которым можно заменить действие данной распределенной силы, модуль которого равен объему этого распределения, а сам вектор делит площади проходящих через него сечений этого распределения пополам распределение сил (массы) ЛА по координатам ССК ЛА - распределение неравномерно распределенных нагрузок внешних сил (массы ЛА) по координатам ССК ЛА распределенные силы - нагрузки, действующие на определенную поверхность распределение скорости ЛА - распределение, в котором точкам ЛА (основания распределе-
ния), в соответствие векторам скоростей данных точек определено множество точек границы данного распределения РМ - реальный маневр ЛА РС - равномерно сбалансированные силы, формирующие линейное ускорение ЛА РУНС - распределение ускорения неравномерно сбалансированных сил РУРС - распределение ускорения равномерно распределенных сил СК - система координат ССК - связанная с ЛА СК с началом координат в носке ЛА
СМ - смоделированный маневр ЛА
ТПР - точка приложения равнодействующей
сил, приложенных к ЛА
ТПРВ - точка приложения результирующего
вектора
фигура скорости ЛА ФD для ЛА в положении D (нижний индекс) - неделимое множество точек, характеризующее границу объемного распределения скорости точек ЛА в момент времени t. Определяется потенциальным перемещением ЛА через единицу времени [4] фигура ускорения ЛА ФD для ЛА в положении D - неделимое множество точек, характеризующее границу распределения ускорения всех точек данного ЛА в момент времени ti [4] ЦМ - центр масс ЛА
эквивалентное преобразование - преобразование геометрического и аналитического представлений движения ЛА, выходные параметры которого позволяют его характеризовать как _ эквивалентное исходному движению ЛА i
ЭЧ - элементарные частицы ЛА g
характерные отрезки - отрезки ЛА, принад- g лежащие осям его ССК ь
[ ]ох ([ ]Атт) - знак отнесения действия к ха- « рактерному отрезку, принадлежащему оси Ox = (ко всему ЛА, определяемому АТТ) |
та ш о
Введение о
Стремительное развитие техники выполне- 8 ния сложных маневров (в т. ч. фигур высшего о пилотажа) современных летательных аппара- о тов открывает большие возможности их при- i менения по боевому назначению. Но такие | маневры выполняются на критических углах
см см о см
< I
со та
г
о со
о.
<и
о
о <и со
см ■ч-ю
с?
см ■ч-ю см
(П (П
атаки и скольжения, поэтому их выполнение под силу лишь опытным летчикам-профессионалам. Повсеместное выполнение таких маневров в автоматических режимах системами управления ЛА, а также беспилотными ЛА с целью достижения лучших условий их применения неприемлемо из-за существующих динамических погрешностей систем управления. Одна из причин таких погрешностей обусловлена использованием недостаточно точных ММ полета ЛА при синтезе их ЗУ. Для выполнения маневрирования ЛА в автоматическом режиме такие погрешности управления необходимо по возможности снизить до предельно малых вследствие невозвратных реакций ЛА, возникающих при превышении ключевых параметров маневра. Решение задачи отчасти позволило бы решить проблему выхода ЛА в автоматическом режиме из области неустойчивого полета в аварийных ситуациях.
Более точное понимание и, как следствие, прогнозирование происходящих процессов динамики ЛА позволит обеспечить компромисс устойчивости и управляемости, улучшить качество переходных процессов управления и стабилизации параметров полета ЛА.
Кроме того, использование существующих ММ полета ЛА, обладающих описанными погрешностями, при анализе маневров ЛА приводит к грубым погрешностям вычислительного эксперимента маневрирования ЛА, ошибкам в интерпретации полученных результатов и, в конечном счете, к возможности не всегда делать точные выводы об изменениях в протекающих процессах динамики полета данного ЛА.
Одной из причин динамических погрешностей управления являются погрешности, определяемые методом аналитического описания полета ЛА. В значительной степени снизить методические погрешности, возникающие из-за того, что синтез ЗУ происходит на основе недостаточно точных ММ динамики движения данных ЛА, возможно путем повышения адекватности ММ, описывающих маневрирование ЛА.
Повысить адекватность ММ полета ЛА для существующих условий принятия эквивалентными СМ и РМ ЛА достаточно проблемно.
Задача решается проще при разработке принципа математического моделирования с использованием более жесткого критерия принятия эквивалентными СМ и РМ ЛА.
Цель исследования - повысить адекватность ММ маневра ЛА путем обоснования принципов математического моделирования с использованием более жесткого критерия принятия эквивалентными СМ и РМ ЛА.
Условия и допущения исследования
1. ЛА принят абсолютно твердым телом (АТТ) с неравномерным распределением массы, определяемым относительно координат его связанной с ЛА системы координат, все элементарные частицы (ЭЧ) которого в той или иной степени подвержены действию распределенных внешних сил.
2. Анализ динамики движения производится для ЛА, определяемого свободным телом [3, с. 10], - АТТ, обладающего возможностью совершения эволюций вращательного движения с любым расположением мгновенной оси вращения, МОУС, МОУУ и МОДС.
3. Для каждого шага математического моделирования маневра ЛА принимается, что интервал времени, в который он происходит, настолько мал, что все параметры движения ЛА (кроме угловых и линейных перемещений) не изменяют значения, определенные на начальный момент этого шага моделирования. Изменения данных параметров происходят только в начальный момент следующего шага моделирования.
4. Анализ маневрирования ЛА производится по распределениям его скорости, угловой скорости, ускорения, углового ускорения, а также действующих на него сил, определенных относительно осей его связанной системы координат (ССК) [3, с. 58-61]. Начало координат ССК выбрано в носке ЛА, поскольку центр масс (ЦМ) ЛА при проведении описываемых исследований не играет такой роли, как при приведении параметров полета к этому центру.
Исключить погрешности измерения (вычисления) параметров вращательного движения можно только определением их относительно той мгновенной оси, относительно
которой происходит изменение данного параметра в РМ. Поэтому определение угловых скорости и ускорения ЛА, а также действие НС целесообразно производить относительно МОУС, МОУУ и МОДС соответственно.
Поэтому вращательное движение ЛА определяется в подвижной локальной (не связанной с положением других предметов) системе координат (ПЛСК 1 типа), определяемой для каждого шага моделирования в его начальный момент времени и в течение этого шага сохраняющей это положение. Началом данной ПЛСК выбран полюс О, определяемый как точка пересечения МОУС и плоскости, в которой происходит угловое перемещение носка ракеты. Положение осей ПЛСК 1 типа: оси Ох - по МОУС, с направлением, определяемым по правилу левой руки (для направления вращения ЛА); оси Оу - от центра О к носку ЛА, положение которого определено в начальный момент шага моделирования; ось Ог - дополняет оси Ох и Оу до правой СК.
При отсутствии угловой скорости ЛА параметры его поступательного движения исчисляются в ПЛСК 2 типа, определяемой положением ССК на начальный момент шага моделирования и остающейся в этом положении до конца этого шага.
5. Действие внешних сил, оказывающих влияние на динамику движения ЛА (ДВС), характеризуется описанными выше неравномерно действующими силами (НС) и равномерно сбалансированными силами (РС). Данные составляющие ДВС определяются из условия, что РС формируют линейное ускорение ЛА аЛА, а НС формируют угловое ускорение ЛА 8ла. Таким образом, для каждого шага моделирования характеристикой ДВС являются параметры равнодействующей и момента этих сил (действующего относительно МОДС).
Исключение замены вращательного движения ЛА его поступательным движением
Требованием к используемым СК и преобразованиям (определяемым принципом моделирования) при разработке ММ полета ЛА является не только возможность удобно представить действие сил на ЛА, но и обеспечить
максимальное соответствие их значениям в РМ. Синтез ЗУ на основе ММ, у которой выходные параметры отличаются от параметров РМ, приводит к тому, что управление ЛА будет осуществляться с динамическими погрешностями. По мнению автора одна из причин подобных погрешностей кроется в использовании принципа приведения сил к определенному центру, при котором происходит замена части вращательного движения ЛА его поступательным движением.
Исследуем адекватность такой замены.
Параметры вращательного движения точек ЛА (угловое ускорение и угловая скорость) [6] связаны с радиусами действий этих параметров. В то время как параметры поступательного движения (ускорение и скорость точек ЛА) не связаны с расстояниями до каких-либо центров движения ЛА.
Существенные различия траекторий поступательного и вращательного движений ЛА и параметров их движения позволяют сделать вывод, что поступательное и вращательное движения ЛА - принципиально разные движения ЛА. Следовательно, их действие необходимо определять без замены одного движения другим.
Различие данных движений ЛА, определяемого АТТ, по условиям возникновения и строгого перехода в зависимости от положения МОУС (у поступательного движения МОУС находится в бесконечности) определяет их как полную группу несовместных событий.
Принцип приведения сил, приложенных i
к АТТ, к определенному центру g
Основоположником существующего принципа g
анализа динамики движения АТТ является пэр ь
Франции Луи Пуансо. В 1803 г. им был сфор- £
мулирован принцип [1] преобразования векто- *
ров, действующих на АТТ, путем их приведе- |
ния к выбранному центру, который позже был g
им применен для решения задач динамики. о
Н. Е. Жуковским проводились исследо- о
вания движения АТТ с определением мгно- а>
венного центра вращения, мгновенного цен- о
тра скоростей (МЦС) и мгновенного центра i
ускорений (МЦУ) [2, с. 242-256]. Но возмож- g но, ввиду простоты принципа приведения
см см о см
< I
со та
г
о со
о.
<и
о
о <и со
см ■ч-ю
с?
см ■ч-ю см
(П (П
векторов к определенному центру [2] при решении задачи динамики полета ЛА им был использован именно этот принцип. В настоящее время принцип остается основным при определении параметров движения всех АТТ.
Описанный принцип Л. Пуансо определяет, что «любое число сил, приложенных к телу, можно привести к одной силе, проходящей через произвольно взятую точку, и одной паре... (описывающей действие НС. - Прим. автора)» [1].
Для определения влияния ДВС на маневрирующий ЛА принцип определяет:
- нахождение главного вектора путем определения равнодействующей и переноса точки ее приложения (ТПР) к центру приведения;
- нахождение главного момента путем определения действия пары, формируемой равнодействующей грс и вектором г РС, ей противоположно направленным. При этом вектор ГрС согласно правилу приведения приложен в ТПР, а вектор Гр* - в центре приведения.
Проанализируем данное преобразование в представлении распределенного действия сил на ЛА с равномерным распределением массы.
ДВС на маневрирующий ЛА проявляется в формировании РС и НС. На рисунке 1а действия данных сил показаны в виде распределений РС ГРС (х) и НС ГНС (х), определенных для значений координат по оси Ох ССК ЛА.
Проанализируем изменения распределения НС ГНС (х) после приведения сил к центру Оприв. После приведения сил к центру Оприв распределение НС ГНС (х) стало разбито на две составляющие, которые также описывают неравномерное ГНСив (х) и равномерно сбалансированное действие Гр£ив (х) данного распределения ГНС(х) (рис. 1б).
В результате такого приведения сил к центру Оприв действие НС ГНС определяется только составляющей НС ГН£ив. При этом вторая составляющая НС Гр£ив, определяющая действие НС на отрезке от ТПР ОТПР до центра приведения Оприв, рассматривается как распределение РС. Наблюдается замена части распределения НС распределением РС.
Кроме того, из-за приведения пространственного движения ЛА к центру Оприв разделение угловой и линейной скорости (а также ускорения и углового ускорения) также производится относительно центра приведения без учета положения МОУС (МОУУ). В этом случае погрешности приведения определяются расстояниями от Оприв до ТПР, от Оприв до точки разделения углового и линейного ускорения и от Оприв до точки разделения угловой и линейной скорости. Чем больше данные расстояния, тем погрешности определения действия НС, углового ускорения и угловой скорости в результате данного приведения будут больше. В классической механике при описании принципа приведения
рнТ (x)
аб Рис. 1. Замена части распределения НС распределением РС
утверждается, что вращательная часть от выбора полюса не зависит [3, с. 128].
Фактором, который компенсирует данное несоответствие, является выбор в качестве центра приведения ЦМ ЛА. Для устойчивого по углу атаки ЛА при небольших значениях управляющей силы ТПР перемещается на незначительном удалении от ее ЦМ. Поэтому в этом случае приведение сил к ЦМ в СМ не дает значительных расхождений с параметрами РМ ЛА.
В [4] и [5] автором статьи выдвинута гипотеза, альтернативная данному принципу Л. Пуансо. Согласно данной гипотезе, в любой момент времени движение ЛА обладает единственными мгновенными центрами вращательного движения, относительно которых угловые перемещение, скорость, ускорение, а также момент внешних сил, соответственно, определяются без погрешностей. В [5] также показано, что при анализе полета ЛА, движущегося с угловым ускорением, отказ от использования переноса пары и перемещения вектора силы вдоль линии его действия исключит связанные с этими преобразованиями методические погрешности, повысит точность прогнозирования параметров полета ЛА.
Выбор критерия эквивалентности преобразований движения ЛА
Для определения адекватности преобразований сил и параметров движения ЛА выберем критерий определения эквивалентности СМ и РМ ЛА. Такое утверждение не должно вызывать сомнений в истинности для его принятия без доказательства.
Распределение траекторий движения точек ЛА во времени является идеальным «интегральным» критерием для всех изменений скорости, ускорения и распределенных сил, действующих на данные точки ЛА. Это обусловлено тем, что для любой точки АТТ распределение ее траектории во времени является результатом двойного интегрирования ее ускорения и интегрирования ее скорости.
Поэтому за такой критерий принимается:
Движения летательных аппаратов являются эквивалентными, если данные ЛА перемещаются синхронно во времени
по одинаковым траекториям L1(x1, y1, z1, t) и L2(x2, y2, z2, t) при любом масштабировании:
L1(xb Уь zb t) = y^ ^ t). (1)
При доказательстве использовались свойства динамики движения АТТ. Из-за жесткой связи между соседними ЭЧ в АТТ данные ЭЧ не могут иметь резкого различия в скоростях и ускорениях их движения, поэтому:
- скорости и ускорения точек, находящихся на одной прямой, могут изменяться только пропорционально изменению их координат в ССК ЛА;
- при единичном повороте АТТ отрезки, соединяющие начальное и конечное положение точек, взятых на любом из радиусов этого поворота, образуют ряд параллельных прямых.
Используя доказательства [4] и [5], сформулируем основное утверждение статьи.
Для любого момента времени маневра ЛА разложение его угловой скорости, углового ускорения и момента внешних сил, действующих на ЛА, относительно любых осей (центров) приведения, которые не совпадают с осями действия этих параметров (МОУС, МОУУ и МОДС соответственно) в реальном маневре, не является эквивалентным преобразованием. Параметры вращательного движения маневрирующего ЛА определить достоверно можно только относительно осей, относительно которых они происходят.
Докажем данное утверждение статьи для каждой из перечисленных осей враща- _ тельного движения ЛА. i
г
ф о
Теорема 1 ь
Формулировка. В любой момент вре- ь мени движения ЛА его угловую скорость я определить достоверно можно только от- *
носительно МОУС - оси, относительно ко- |
„ я
торой в данный момент происходит враще- g
ние этого ЛА. Приведение вращательного о
движения ЛА относительно всех остальных о
центров (осей) не является эквивалентным а>
преобразованием. о
Постановка задачи. Для доказатель- i
s
ства изменения параметров вращательного g движения ЛА при использовании приведения —
см см о см
к определенному центру исследуем перемещение его отдельных точек. Исследование проводится для двух случаев вращательного движения ЛА в ПЛСК 1 типа из положения Б в положение Е за интервал времени At, настолько малый, что угловые скорости точек ЛА на нем можно принять неизменными (рис. 2а). В случае исходного движения ЛА совершает вращательное движение относительно центра 02* (анализируется в ПЛСК 1 типа, связанной с этим центром).
В случае приведенного движения ЛА совершает перемещение от тех же начальных к тем же конечным точкам ЛА, но с векторами скоростей его точек, приведенными к центру 01*. Для каждого из движений интервал времени Дt разбивается на т* равных участков (где т* > 10) и отслеживаются траектории (крос-сировки) перемещения этих точек 01 и 02*.
Доказательство 1. Сравним кроссиров-ки исходного и приведенного движений ЛА точек 01* и 02* на интервале моделирования Дt в увеличенном масштабе времени с Дtp = 0,Ш (рис. 2а).
Для движения, приведенного к центру 01*, приведение означает, что вращательное движение рассматривается в ПЛСК 1 типа, связанной с 01*, так, как будто оно происходит относительно этого центра, а сам центр 01' совершает только поступательное движение,
Траектория точки 02* в таком движении приведена на рисунке 2б (траектория показана в масштабе Дtp = 0,1 Д^. В случае исходного движения это траектория движения, состоящая из этой точки 02*.
Для точки 01* в плоскопараллельном движении, приведенном к данной точке, это поступательное перемещение на £ [м] к точке 01**. В случае исходного движения - это вращательное движение точки 0!* относительно центра 02*.
Вывод. Траектории исходного и приведенного движений ЛА не совпадают, при этом их несовпадение не связано с изменением дискретности, а отличается характером кросси-ровок отдельных точек ЛА. По критерию (1) на интервале времени, на котором параметры полета ЛА остаются постоянными, исследуемые движения от одинаковых исходных к одинаковым конечным точкам ЛА определяются как неэквивалентные.
Доказательство 2. Сравним параметры, приведенного на рисунке 2а плоскопараллельного движения ЛА с параметрами аналогичного движения, приведенного к точке 0{\
В случае исходного движения вращательное движение точки 0!* происходило с угловой скоростью х / Дt [рад/с] без линейной скорости.
В случае приведенного движения угловая скорость точки 01* относительно точки
< I
со та
г
о со
о.
<и
о
о <и со
см ■ч-ю
с?
см ■ч-ю см
(П (П
О'
о2*
аб Рис. 2. Кроссировки точек исходного и приведенного движения ЛА
01 равна нулю, а линейная скорость движения этой точки составляет ^ / Дt [м/с].
Различие линейных и угловых скоростей отдельных точек ЛА в двух рассмотренных случаях указывает на то, что в результате приведения исходного движения к новому центру полученное приведенное движение стало отличаться от исходного движения ЛА. Вследствие того что данные точки ЛА движутся со скоростями отличными от исходных угловых и линейных скоростей движения, исходное и преобразованное движения определяются как неэквивалентные.
Замечание. Отметим, что с точки зрения допущений существующей динамики описанные плоскопараллельные движения ЛА, рассматриваемые на одном интервале моделирования с разделением вращательного и поступательного движений относительно разных полюсов, принимаются как эквивалентные движения с разными центрами приведений.
Вывод. Допущение классического описания движения ЛА о возможности представления траекторий вращательного движения точек ЛА как траекторий поступательного движения данных точек для одного и того же рассматриваемого случая приводит
к изменению значений линейных и угловых скоростей точек ЛА. Отметим это как один из недостатков принципа приведения сил к определенному центру, поскольку зависимость параметров одного и того же движения ЛА от положения центра приведения наблюдаться не должна.
Таким образом, любой маневр ЛА может рассматриваться как вращательное движение относительно МОУС, с течением времени изменяющей свое положение в результате поступательного движения ЛА. В качестве примера на рисунке 3 показаны кроссировки перемещения МОУС для маневрирующего истребителя [4] (в данном случае МОУС показаны как мгновенные центры скоростей вращательного движения этого истребителя).
На рисунке 3 также показаны направления векторов мгновенной скорости ЦМ ЛА й ^ й10, определяемых как параметры мгновенных угловых скоростей движения этого центра в разные моменты времени выражением = Iй / г, (где / = 1,10).
На рисунке 3 видно, что в некоторых случаях МОУС уходит в бесконечность, ЛА переходит в состояние поступательного движения. На рисунке такие состояния определяются
Рис. 3. Кроссировки проекции МОУС ЛА на плоскость вращения при выполнении маневра [4]
е
о р
т с о т
тке
а р
а
ш
о Ч е л с с
к с е
у
и м с о К
см см о см
< I
со та
г
о со
точками пересечения прямых Ь, с, ё, е и траектории ЛА.
Качественное и количественное отличия распределений ускорений поступательного и вращательного движения ЛА определяют раздельное рассмотрение сил, которые их формируют. Поэтому далее распределение углового ускорения ЛА (ускорения вращательного движения ЛА) по координатам ССК данного ЛА будем определять как распределение ускорения неравномерно действующих на ЛА сил (РУНС). Соответственно распределение ускорения поступательного движения по координатам его ССК - как распределение ускорения ЛА равномерно сбалансированных сил, действующих на ЛА (РУРС).
Для более наглядного представления распределений скоростей (ускорений) точек ЛА помимо их границ в виде фигур скоростей (ускорений), например [Ф ^]Атт, будем показывать границы распределений скоростей (ускорений) для характерного отрезка ЛА [Ф^ох, принадлежащего оси 0х ССК ЛА.
Для определения алгоритма получения РУНС рассмотрим получение подобного распределения на более простом примере -на примере получения в ПЛСК 2 типа границы РУРС [Фал]0х ЛА, совершающего только поступательное движение (рис. 4а) в момент времени По исходным данным в моменты времени t1, ^ ЛА находился в положении фигур F, G, H соответственно, где ^ - t1 = t1 - t0 = Дt.
Исходя из того, что ускорение точки ЛА определяется из условия
, ч г Уу2(х)-Уу1(х)
■ 1ш1
ДГ->0
РУРС его характерного отрезка в момент времени ^ при его нахождении методом Эйлера (для чего значения скорости должны использоваться в масштабе - 1 / Д^ для ЛА в положении
F определяется выражением (рис. 4б): ]_
М
Теорема 2
Формулировка. Для любого момента времени движения ЛА разложение углового ускорения ЛА относительно любых центров (осей) приведения, которые не совпадают с МОУУ в РМ, не является эквивалентным преобразованием. Угловое ускорение маневрирующего ЛА определить достоверно можно только относительно МОУУ.
Постановка задачи. По трем положениям ЛА (ракеты) (фигуры Л, B и C, рис. 5а), зафиксированным в моменты времени t1 и для момента времени ^ определяется распределение его ускорения по координатам ССК ЛА в ПЛСК 1 типа.
Для однозначного определения вращательного движения ЛА согласно теореме 1 методом [2, с. 240] для нахождения положения ПЛСК 1 типа для ^ и t1 определены положения МОУС ЛА, проходящие перпендикулярно плоскостям вращения а и в, с их пересечением данных плоскостей в точках 0Х и 02 соответственно. Характерные отрезки фигур Л, В и С на рисунке 5а показаны как отрезки КК, К1Ы1, К2Ы2.
Необходимо доказать, что неизменное на интервале времени Дt = t1 - ^ угловое ускорение ЛА без погрешности может определяться только относительно единственной МОУУ.
.
<и
о
о <и со
см ■ч-ю
с?
см ■ч-ю см
(П (П
аб Рис. 4. Представление действия линейного ускорения ЛА в виде РУРС
1
L
а б
Рис. 5. Представление действия углового ускорения ЛА в виде РУНС
Доказательство. Поскольку потенциальное перемещение точек ЛА на рисунке 5 происходит как вращение относительно определенного центра, рассматриваемые распределения скорости нельзя считать прямолинейными. В этом случае больше подошел бы термин «секторные распределения вращательного движения ЛА», которые далее по тексту будем выделять верхним индексом сект, например
гжшН сект [ФШЮх .
Поэтому, по аналогии с понятиями распределений скорости и ускорения, фигур скорости и ускорения ЛА [4], введем понятия распределения угловой скорости и углового ускорения ЛА, фигур угловой скорости и углового ускорения ЛА.
Распределение угловой скорости ЛА -секторное распределение параметра углового перемещения ЛА, в соответствии с которым точкам текущего положения ЛА (основания) согласно их потенциальному угловому перемещению через единицу времени определено пространственное множество точек границы этого распределения.
Распределение углового ускорения ЛА -секторное распределение параметра углового перемещения ЛА, в соответствии с которым точкам текущего положения ЛА (основания) согласно прогнозируемому приращению угловой скорости точек ЛА через единицу времени определено пространственное множество точек границы этого распределения.
Фигура угловой скорости Ф™ с основанием Н (нижний индекс символа ФН) в момент времени - граница распределения угловой скорости ЛА как неделимое множество точек, которое характеризует положение этого ЛА через единицу времени после прогнозируемого углового перемещения.
Фигура углового ускорения ФБ с основанием Б в текущий момент времени Х - граница распределения углового ускорения ЛА как неделимое множество точек, которое определяет прогнозируемое приращение угловой скорости точек ЛА через единицу времени.
Для исключения замены параметров вращательного движения ЛА параметрами поступательного движения его анализ разложением вектора ускорения точки ЛА на касательное ат и нормальное ап ускорения данной точки [3, с. 110] не используется. При данном анализе используются описанные секторные распределения угловых скоростей и ускорений ЛА.
По разности распределений угловых скоростей [ФЙ^сГ и [Ф^]оХкт, взятой в масштабе 1 / ДХ (где ДХ - интервал времени между определением данных угловых скоростей), определим аналитическое выражение для РУНС [Ф1исг (рис. 5б):
r>fvCü2 псекг г>т*со1 псект L4^ SOx ~ L4^ iOx
Дt
где [Ф|УОХкт - РУНС с отрезком Ш в основании (с характерным отрезком по оси Ох ЛА,
е
о р
т с о т
тке
а р
а
ш
о ч е л с с
к с е
у
и м с о К
см см о см
< I
со та
г
о со
о.
<и
о
о <и со
см ■ч-ю
с?
см ■ч-ю см
(П (П
находящегося в момент времени ^ = 1 в положении А).
Определим положение РУНС [ФКм]«^ в пространстве.
1. В плоскости а в масштабе значений угловой скорости 1 / Дt определим распределение -[Ф^У сеК* / Д t как обратное распределению угловой скорости [Ф1^](9хКТ (на рисунке 5б показано синей стрелкой). Полученную границу распределения обозначим отрезком К3Д3.
2. Поворотом по часовой стрелке приведем распределение угловой скорости [ФК2N]оТт, определенное в масштабе 1 / Дt, до совмещения его основания с отрезком К3М3 в плоскости у, при у || в (на рисунке 5 показано зеленым цветом). Получим секторное распре-
деление [ФКм ]«Г / Д^
'2 "1 сект зМэ]
3. Суммированием секторных распределений угловой скорости -[ФК#]аГ / Дt и [ФК2эЛг3](0Хкт / Дt в плоскости х получим распределение [ФК^ОГ по значениям координат ССК ЛА с основанием КМ (на рисунке 5 сектор распределения выделен красным цветом).
4. Используя метод, изложенный в теореме [2, с. 240], для t0 определим положение МОУУ РУНС [ФАЧлА^ с основанием А. Для данного случая движения ЛА (рис. 5 а) рассматриваемая ось, проходящая через точку Ое1, располагается перпендикулярно плоскости X (рис. 5б).
Для доказательства от противного предположим, что [ФАЧлТ" может быть разложено на распределения линейного [Фа1]*лА и углового ускорения ЛА [Фе1]ЛсАкт* относительно произвольного центра:
[ФАЧлТ = [ФЕ1]лсАкт* + [фо1]*ла, (3)
Однако, исходя из условия, что [ФА1]сект
лА
может быть разложено на составляющие
Гфб1 псект
глш2 псекг псекг
^ ¡Ох ~ АОх
Дг
сект
его составляющие [Ф<А1]лсАкт и [Ф'2]лАкт тоже можно разложить на суммы аналогичные (3):
Г(Т»®1п сект _ глйИ сект* i г(ЪУП [ФА ]лА [фА ]лА + [ФА ]
тт ГгТ»®2п сект _ ГгТ»®2 сект* I г(Г\У2~\ и [ФА ]лА = [ФА ]лА + [ФА ]
лА
лА
(4)
и [Ф'2]лАкт на составляющие поступательного и вращательного движения приведет к тому, что траектории преобразованного движения (4) не будут совпадать с исходными траекториями. Следовательно, выражение (3) определяет движение лА с изменением исходной траектории движения, что противоречит критерию (1) и определяет доказательство теоремы 2 от противного.
По аналогии со следствием 2 для теоремы 2 можно показать изменение параметров исходного и приведенного маневров лА.
Приведем произвольно взятое РУНС (например, [ФА^лл") маневрирующего лА к центру приведения, не принадлежащего МОУУ этого РУНС (например, согласно выражению (3)). В результате получим распределения линейного и углового ускорения по координатам данного лА. Полученное и приведенное состояния движения лА, определяемые относительно разных МОУУ, имеют разные значения угловых ускорений точек лА (елтАк ~ объему [ФАЧлА^ и в;
■^прив
Однако, как показано в теореме 1, разложение угловых скоростей лА [ФАЧлА^
объему
[Ф|1]лс.Акт*). При угловом перемещении от одинаковых исходных положений к одинаковым конечным положениям лА это определяет, что характеризуемые ими состояния движения данного лА нельзя считать эквивалентными. Приведение РУНС к центру приведения, не принадлежащему МОУУ этого РУНС, приводит к изменению состояния движения лА, что подтверждает доказательство теоремы 2.
Таким образом, любое приведение РУНС относительно выбранного центра приведения, не принадлежащего оси, относительно которой в данный момент действует угловое ускорение точек этого лА в РМ (т.е. относительно текущей МОУС), приводит к неэквивалентности разложения этого РУНС.
Поскольку полученное РУНС, приведенное на рисунке 5б (РУРС, приведенное на рис. 4б), получено единичным вращением (единичным параллельным переносом), можно сделать вывод о том, что для любого момента времени границы этих распределений по отношению к основанию (фигуре самого лА) по соответствию размеров и, как следствие, и по форме не изменяются.
Сохранение размеров фигур линейных и угловых скоростей, а также линейных и угловых ускорений дает основание утверждать, что все значения распределений скоростей и ускорения для точек ЛА, лежащих на одной прямой, описываются пропорциональной зависимостью от координат основания. В частности, границы РУРС ЛА, определенные для характерных отрезков ЛА, можно характеризовать положением прямых аРс(х), арс(у), архс^), содержащих эти границы:
«РДх) = ар/(у) = = (5)
где арО, аРО, аРо - значения линейного ускорения границы распределения во всех точках данного распределения.
По аналогии, границы РУНС, определяемые для характерных отрезков ЛА, можно характеризовать положением прямых а1ус(х), аЧс(у), аИх(2) содержащих эти отрезки:
^(Хк = кн;(х)(х - хо)+анУ0;
ан/(у), = ¿^(У - Уо) + ан/о; (6)
а%с(г), = & - ¿о) + анх0,
где ¿1у<(х), ¿нСУ), ¿нхС(г) - коэффициенты пропорциональности прямых, содержащих границы РУНС характерных отрезков; а1ус0, а|0, а^хО - значения границ РУНС в точках ЛА с нулевыми значениями абсцисс ССК; х0, у0, - значения сдвига прямой от положения границы распределения над его основанием.
Перед тем как изложить следующую лемму, сделаем несколько пояснений по используемым в ней терминам и обозначениям.
В лемме под линеаризацией значений распределения параметра по координатам ССК ЛА понимается преобразование распределения, определяющее сглаживание их значений по значениям до линейных функций, определяемых от координат ССК ЛА, с сохранением параметров воздействия исходного распределения на динамику ЛА.
Под сохранением параметров воздействия на динамику ЛА понимается сохранение направления и модуля равнодействующей для равномерно сбалансированного распределения параметра по осям ССК ЛА и модуля, плоскости и оси действия для секторного распределения параметра по данным осям.
В интересах дальнейшего исследования РС и НС на ЛА представим второй закон Ньютона в размерности [м/с2] для равномерного и неравномерного действий внешних сил на ЛА и Е§х(0), определяемых в 1-й мо-
мент времени:
рс /_ ^со. нс _ F^d)
(7)
где аЛд('), аЛА.(0 - ускорения ЛА от действия РС и НС в 1-й момент времени, тЛА(1) - масса ЛА, зафиксированная в 1-й момент времени.
Отметим, что пространственными распределениями можно представить не только РУРС (2) и РУНС (4), но и распределения РС, НС и массы ЛА, с тем отличием, что границы последних распределений не всегда можно представить в виде сплошных неделимых фигур.
Это дает возможность для ЛА, находящегося в 1-й момент времени в положении фигуры Т, представить (7) в виде РУРС [Фт1]0х и РУНС [Фансг]0хкт с характерным отрезком в основании (в размерности [м/с2]):
[ФПаМФ ru
[Ф По,- IV
; [фтТс
• (8)
где [ФрЪ*, [Фнс ЧОГ, [Ф ¥']ох - распределения РС, НС и массы ЛА с характерным отрезком в основании, определенные в 1-й момент времени.
Лемма
Для любого произвольного момента времени 1 распределения НС [ФЕнсЧлГ и РС [ФЕг]ЛА по координатам ЛА, а также распределение массы ЛА [Ф™г]ЛА по этим координатам ЛА при влиянии на динамику маневра ЛА, определяемого АТТ, линеаризуются по значениям функций этих распределений.
Доказательство. Распределения углового и линейного ускорений по координатам ССК ЛА [Ф?нс 1 ]0ехКТ и [Фа%х с границами, характеризуемыми линейными функциями (6) и (5), согласно (8), являются частными, что определяет пропорциональность распределений делимого и делителя, т.е. пропорциональность распределений НС и РС и распределения массы ЛА.
е
о р
т с о т
тке
а р
а
ш
о ч е л с с
к с е
у
и м с о К
Парадокс отношений правых частей выражений (8) в очевидности, что неравномерное распределение РС [Фр]лА и неравномерное распределение НС [Фу нс г]лАкт во все моменты времени маневрирования не могут быть пропорциональны лА неравномерному распределению массы [Ф т' ]лА.
Объяснение данного парадокса лежит в понимании того, как для неравномерного ДВС по конструкции маневрирующего лА границы его РУНС и РУРС для любого момента времени оказываются без изменения формы и размеров по отношению к их основанию.
логическим объяснением этого несоответствия являются свойства лА как АТТ. В данном случае это свойство линеаризации действия на лА НС и РС так же, как и линеаризации инертного влияния неравномерно распределенной массы лА по координатам его ССК с сохранением свойств исходных распределений относительно этих координат. Свойство распространяется на НС и РС во все моменты времени рассматриваемого движения лА, что и требовалось доказать по условиям леммы.
Распределения НС и РС, определенные с учетом их линеаризации по координатам ^ ССК лА, будем называть действующими распределениями РС и НС этого лА, а распре-^ деление массы лА с учетом этого процесса -— влияющим распределением массы данного лА. '«и Свойство линеаризации (выравнива-
ния) значений распределения сил по коорди-Д натам ССК лА можно объяснить взаимны-5 ми компенсациями и перераспределениями ^ неравномерно действующих внешних сил
§ по конструкции лА. аз
щее распределение его массы
Решение задачи динамики для существующих в РМ неравномерных распределений НС, РС и массы сложно. Использование свойства линеаризации данных распределений делает определение РУРС и РУНС лА при решении данной задачи более простым и удобным.
На примере движения лА, происходящего в плоскости Оху, дадим полученным результатам геометрическую интерпретацию (рис. 6, 7).
Обозначения на рисунках 6 и 7: к • т(х), к • т*(х) - влияющие распределения массы лА, представленные для разных случаев. Для наглядности показаны в масштабе по значениям функции, где к - масштаб функций т(х) и т*(х); Унс(х)сект и У*(х) - секторное и линейное распределения действующих сил, с характерным отрезком лА в основании.
Условием возникновения только линейного ускорения лА с равномерно влияющим по его координатам распределением массы т(х) является равномерное действующее по этим координатам распределение сил У (х) (рис. 6). При неравномерно влияющем по координатам лА распределении массы т*(х) таким условием является пропорциональное по значениям распределения т*(х) действующее распределение сил У (х) (рис. 6). Отметим, в данном случае У *(х) - линейное, а не секторное распределение, поскольку созданное им ускорение линейное.
Необходимым условием перехода от движения с линейным ускорением к движению с угловым ускорением лА с равномерно влияющим распределением по его координатам массы является действие любой НС. Например, для лА с равномерно влияющим распределением по оси Ох ССК массы т (х), представленного на рисунке 6, это может быть действие НС с действующим распределением
р нс(х)сект (рис. 7).
Одним из путей перехода от движения с угловым ускорением лА а(х)сект к движению с линейным ускорением при действии на лА этого же распределения силы У нс(х)сект является изменение распределения массы лА до неравномерно влияющего по его координатам т*(х) (рис. 6). При этом значения влияющего распределения массы лА должны быть пропорциональны соответствующим значениям
не-
действующего распределения сил Г *(х). Таким образом, движение ЛА с угловым ускорением может выродиться, т.е. перейти к частному случаю полета с МОУУ, расположенной в бесконечности. Такое секторное распределение будем называть вырожденным секторным распределением параметра вращательного движения ЛА.
Использование распределений скорости (угловой скорости) ЛА позволяет наглядно продемонстрировать сложные переходы динамики полета с изменением данных скоростей. На рисунке 8 показан результат действия углового ускорения ?2 (относительно центра Ое2) как переход от движения с угловой скоростью (относительно исходного центра Ою2) к движению ЛА с вырожденным распределением угловой скорости ЛА ю2, т.е. к движению ЛА с линейной скоростью У^3.
Если не брать в расчет сферичность распределения НС и учитывать только положение границы этого распределения, создается ошибочное представление о том, что можно выделить в этом распределении постоянную составляющую. Это и определяет предпосылку классической динамики [3, 6] о возможности представления действия НС как одновременного действия сил, формирующих линейное и угловое ускорения ЛА (т.е. НС и РС). Однако свойства секторного распределения значительно отличаются от линейного распределения, что исключает такую возможность.
Отметим, что разделение ДВС на действие НС и РС возможно только для известного распределения массы ЛА на текущий момент времени. С изменением распределения массы ЛА НС могут перейти в РС, и, наоборот.
Теорема 3
Формулировка. Действие НС на маневрирующий ЛА с текущим распределением массы m (x, y, z) по осям его ССК в любой момент времени происходит только относительно одной МОДС.
Постановка задачи. Для движения ЛА с неравномерным распределением массы m (x, y, z) по осям его ССК и известным действующим распределением НС в текущий момент времени необходимо доказать единственность положения МОДС данных сил.
Для доказательства можно использовать свойства рассмотренного в теореме 2 РУНС или использовать свойства сил, действующих на АТТ.
Доказательство 1. Доказать единственность существования МОДС для каждого момента времени маневрирования ЛА можно, используя свойство действия его РУНС относительно единственной МОУС для этого момента времени. Согласно выражению (8) единственным распределением, которое может изменить положение МОДС, является распределение массы этого ЛА.
Рассматривая разные варианты, несложно показать, что для равномерного влияющего распределения массы ЛА МОДС совпадает с МОУУ. Пропорциональное изменение длин дуг действующего распределения НС и РУНС, соединяющих точки основания и границы этих
Рис. 7. Действующее распределение НС ЛА с равномерным влияющим распределением массы
Рис. 8. Переход ЛА к движению с линейной скоростью (показано в размерности [рад/сек])
е
о р
т с о т
тке
а р
а
ш
о ч е л с с
к с е
у
и м с о К
секторных распределений, приводит к тому, что положения МОУС и МОДС совпадают (рис. 7). Это можно показать аналитически, определяя пропорциональное убывание прямой, содержащей границу действующего распределения НС, и прямой, содержащей границу РУНС характерного отрезка лА, до общей для этих прямых точки с ординатой равной 0.
В случае неравномерного влияющего распределения массы полученное распределение сил (8) определяется не только «растянутым» по оси ординат распределением ускорения (значениям углового ускорения), но и пропорциональным смещением границы этого изменения по оси абсцисс. Поэтому для этого распределения сохраняется параллельность отрезков, соединяющих точки основания и границы распределения, что, как следствие, определяет существование единственной МОДС q этого распределения, не совпадающей с МОУС q. На рисунке 9 для распределения заданного ускорения НС лА а(х)сект, приведенного на рисунке 7, показано секторное распределение НС лА Унс*(х)сект, определенное для лА с неравномерным распределением влияющей массы к • т' (х).
Доказательство 2. Единственность существования МОДС можно доказать, поэтапно ^ определяя процессы формирования НС и РС ™ из суммарного ДВС на маневрирующий лА. а При определении происходящих про-
— цессов будем исходить из того, что лА яв-о ляется совокупностью ЭЧ, настолько малых, ^ что плотность каждой из них можно считать Д однородной. При этом все исходное неравно-5 мерное ДВС на лА, формируемое тремя видами разнонаправленных сил разной природы,
номерным влияющим распределением массы
можно разложить на совокупность векторов, каждый из которых приложен к центру ЭЧ лА и имеет модуль, характеризующий суммарное действие описанных сил на эту ЭЧ.
Поскольку для каждого момента времени в разных направлениях существует результирующее действие сил всех ЭЧ лА, логично предположить, что по одному из направлений их результирующее действие станет больше, чем в остальных, и его преобладающее действие станет формирующим для результирующего действия описываемых сил.
логично также предположить, что в результате суммарное действие сил по другим направлениям скорректирует параметры формирующего вектора до параметров результирующего вектора УРС всей системы сил, действующих на данный лА.
Отметим, что не все векторы сил, действующие на ЭЧ лА, участвуют в создании результирующего вектора УРС, который сформирует линейное ускорение лА в данный момент времени.
Поэтому выражение
п —«.
где - сумма всех векторов сил действу-1
ющих на ЭЧ лА, не может быть использовано для определения его модуля.
Принимая, что нам известны параметры результирующего вектора УРС, т.е. исходя от обратного, оценим, какие из векторов сил ЭЧ участвуют в его создании.
В первую очередь в создании результирующего вектора УРС принимают участие все векторы, с направлением содержащих их прямых, пересекающих точку приложения результирующего вектора УРС (ТПРВ).
Кроме того, к таким векторам необходимо отнести часть совокупности векторов сил с направлением содержащих их прямых, которые не пересекают ТПРВ. Такие векторы в формировании результирующего вектора УРС участвуют только в случае, если для каждого из них существует вектор равного модуля, действующий симметрично текущей ТПРВ ОТПРВ (например, векторы, приложенные к точкам
Ь и Ь1 на рисунке 10), или вектор, имеющий противоположно направленный относительно 0ТПРВ момент.
Другая часть этой совокупности векторов с направлением содержащих их прямых, не пересекающих ТПРВ, не участвует в формировании равнодействующей, если их момент не компенсируется другими векторами (т.е. которые являются НС). На рисунке 10 в качестве такого вектора показан вектор силы, приложенный к точке Q, не участвующий в формировании результирующего вектора РРС. Данный вектор, кроме того что направление содержащей его прямой не пересекает ТПРВ 0ТПР, не имеет парного вектора (из векторов, представленных на рисунке 10) относительно центра ОТПР.
Отметим, что существует множество векторов описанного типа, которые участвуют одновременно и в формировании неравномерного, и в формировании равномерно сбалансированного действия сил на ЛА.
Поскольку из известного распределения ДВС при одновременном действии НС определить парность всех векторов, образующих действие РС, сложно, вычислить положение ТПРВ сразу не представляется возможным.
В этом плане проще сначала определить положение МОДС всего ДВС (исходной совокупности разнонаправленных сил), как и момент, создаваемый относительно этой оси. После чего определять распределение РС (и, как следствие, параметры результирующего вектора сил РРС ЛА), вычитая из общего распределения ДВС действующее распределение НС.
Докажем, что все векторы НС формируют суммарный момент ЛА именно относительно ТПРВ. Исходя из факта, что каждый из векторов НС можно перевести в вектор РС, дополнительно приложив к ЛА ему парный вектор с противоположно направленным моментом относительно ТПРВ, следует, что все векторы НС действуют относительно ТПРВ. Например, вектор НС на рисунке 10, приложенный к точке Т, станет вектором РС, если к ЛА приложить дополнительный парный вектор, приложенный к точке Т*.
Продолжая описание процесса, отметим, что сформированный суммарный момент НС не может действовать только относительно точки (ТПРВ). В любой рассматриваемый момент времени в одной из плоскостей действие НС будет иметь больший момент, чем в других плоскостях. Это приведет к тому, что его действие, скорректированное действием НС в других плоскостях, создаст результирующий момент всех НС. Одним из параметров данного результирующего момента является положение плоскости его действия, которая содержит ТПРВ и которая определяет положение МОДС - единственной оси действия данного момента.
Рассмотренные свойства НС и РС указывают на то, что разная природа их возникновения, разные свойства возникающих от их действия углового и линейного ускорения ЛА исключают возможность замены действия части НС действием части РС при приведении к оси (центру) приведения (не совпадающему с МОДС). Получаемый при таком приведении момент внешних сил будет действовать относительно новой МОДС с модулем не равным исходному, что указывает на его неэквивалентность исходному моменту внешних сил, действующих на ЛА. Таким образом, МОДС является единственной осью действия НС для данных условий. Теорема 3 доказана.
Результатом выводов теорем 1-3 является доказательство основного утверждения статьи: СМ и РМ ЛА могут быть приняты эквивалентными по критерию (1), а их параметры одинаково изменяющимися во времени только при совпадении их МОУС, МОУУ, и МОДС во все моменты времени движения ЛА.
Принцип моделирования динамики полета ЛА, отвечающий данному утверждению,
Рис. 10. Формирование НС и РС
е
о р
т с о т
тке
а р
а
ш
о ч е л с с
к с е
у
и м с о К
см см о см
< I
со та
г
о со
о.
<и
о
о <и со
см ■ч-ю
с?
см ■ч-ю см
(П (П
заключается в том, что моделирование маневра лА для каждого момента времени производится не относительно ЦМ лА, а с учетом положений МОДС, МОУУ и МОУС движения данного лА.
Алгоритмически определение положения лА в пространстве через шаг моделирования при известных исходных данных (о его положении, положении его МОУС, распределении внешних сил и массы по координатам его ССК) означает:
- определение МОДС и момента НС лА относительно этой оси;
- определение действующих распределений НС и РС, влияющего распределения массы лА;
- определение параметров РУРС и РУНС относительно своей МОУУ;
- определение распределений приращений линейной и угловой скоростей лА за единицу моделирования (как следствие, определение положения МОУС по отношению к найденному распределению угловой скорости лА);
- определение положения МОУС в пространстве движения лА путем перемещения МОУС в предыдущий момент времени на значение приращения линейной скорости лА за шаг текущего моделирования;
- определение положения лА в пространстве его движения путем:
а) наложения распределения приращения угловой скорости за шаг моделирования на положение лА в пространстве его движения с совмещением:
- найденной МОУС в пространстве с МОУС распределения;
- точки носа лА в пространстве в предыдущий момент времени с точкой носа лА в распределении;
б) поворот лА от положения описанного совмещения на распределение приращения угловой скорости за шаг моделирования.
Принципиальное отличие выходных данных предлагаемого моделирования в том, что полученные распределения действующих на лА сил по координатам его ССК, а, следовательно, аналогичные распределения скорости и ускорения по координатам данного лА не зависят от положения центра приведения сил, действующих на лА.
Список литературы
1. Пуансо л. Начала статики. Перев. с фр. А. С. Долговой. Под ред. инж. А. Н. Долгова. Изд-во: Научно-техн. отдел ВСНХ. Научно-техн. изд-во. Петроградское отделение. Петроград, 1920. С. 11, 47-49.
2. Жуковский Н. Е. Полное собрание сочинений. лекции. Вып. 3. Теоретическая механика. Курс высшего технического училища. Под ред. проф. В. П. Ветчинкина М.-л.: Гос. изд. оборонной промышленности. 1939. С. 240-256.
3. Тарг С. М. Краткий курс теоретической механики: учеб. для втузов. 11-е изд., испр. М.: Высш. шк., 1995. С. 10, 38-39, 58-61, 108-110, 128, 153.
4. Писковацкий А. А. Метод определения положений центров вращательного движения при исследовании полета маневрирующего летательного аппарата // Информационно-измерительные и управляющие системы. 2020. Т. 18. № 4. С. 13-22. DOI: 10.18127/]20700814-202004-02
5. Писковацкий А. А. Одна из причин методических погрешностей решения задачи динамики полета летательного аппарата, маневрирующего с угловым ускорением // Информационно-измерительные и управляющие системы. 2020. Т. 18. № 6. С. 37-50. DOI: 10.18127/]20700814-202006-05
6. Голдстейн Г., Пул Ч., Сафко Дж. Классическая механика. М. - Ижевск. «Регулярная и хаотическая динамика». Ижевский институт компьютерных исследований. 2012. С. 174-186.
Об авторе
Писковацкий Андрей Анатольевич - кандидат технических наук, инженер-конструктор 1-й категории Публичного акционерного общества «Долгопрудненское научно-производственное предприятие», Долгопрудный, Московская область, Российская Федерация.
Область научных интересов: динамика полета, системы управления маневренными летательными аппаратами.
Improving the adequacy of aircraft manoeuvre simulation by determining the positions of the instantaneous axes of its rotational motion
Piskovatsky A. A.
PJSC "Dolgoprudny Research and Production Enterprise" Dolgoprudny, Moscow Region, Russian Federation
The potential of developing a system for automatic control over manoeuvres at critical stall and sideslip angles, and aircraft (AC) automatic recovery from emergency flight modes, is limited by the presence of the control system's dynamic errors, which, in case their critical level is exceeded, lead to a failed manoeuvre and may result in the loss of AC in an unstable flight situation. One of the causes of such errors are the measurement-method errors occurring due to the fact that synthesis of the control laws (CL) takes place on the basis of mathematical models (MM) of those AC motion dynamics that are not accurate enough. The presence of the measurement-method control errors, as determined by the principle of flight mathematical simulation, will not allow to determine CL corresponding to these flight processes and use them for implementation of such manoeuvres in the automatic flight modes.
The objective of the study is to improve the adequacy of the MM of AC manoeuvring through development of new principles of its configuration using a criterion with more stringent requirements to determining the equivalence of a simulated and a real manoeuvre (hereinafter SМ and RМ) of the AC.
By way of proving three theorems, the paper formulates the conditions for assuming AC's SМ and RМ as equivalent, and substantiates a possibility to reduce control system's dynamic errors in high-speed processes of AC manoeuvring through the use of a novel principle of mathematical modelling of such a flight. Despite the theoretical and investigative orientation of the paper, the expected outcomes from the improvement of MM adequacy due to a more accurate prediction of the AC manoeuvre parameters and generation of control signals under such manoeuvring that are free from the measurement-method errors will be of great importance for practical application.
Keywords: model adequacy, equivalence of states, trajectory coincidence, instantaneous axis, acceleration distribution
Information about the author
Piskovatsky Andrey Anatolyevich - Cand. Sci. (Engineering), Design Engineer of the 1st category, PJSC "Dolgoprudny Research and Production Enterprise", Dolgoprudny, Moscow Region, Russian Federation. Science research interests: flight dynamics, manoeuvrable aircraft control systems.
