CC BY
1161
Melnikov Peter Nikolaevich, candidate of technical science, leading researcher, peter@olvs.miee.ru, Russia, Moscow, Zelenograd, National Research University of Electronic Technology,
Sazonov Aleksey Alexadrovich, candidate of technical science, head of department, als@olvs.miee.ru, Russia, Moscow, Zelenograd, National Research University of Electronic Technology
УДК 623.55.025
РАСЧЕТЫ ТРАЕКТОРИИ ПОЛЕТА АРТИЛЛЕРИЙСКОГО
СНАРЯДА
П.Н. Мельников, М.А. Терпигорев
Приведена методика перехода от табличного способа расчета траектории полета артиллерийского снаряда к расчету траектории полета снаряда способом численного интегрирования системы дифференциальных уравнений задачи внешней баллистики.
Ключевые слова: артиллерийская зенитная стрельба, вычислительная система.
Одна из основных задач, решаемых вычислительной системой артиллерийского зенитного комплекса - это задача встречи снаряда с целью. В свою очередь, ядро задачи встречи составляет задача прогнозирования траектории полета снаряда на произвольный отрезок времени в будущем. В большинстве зенитных комплексов, принятых в настоящее время на вооружение, прогнозирование траектории полета снаряда решается табличным способом. Для каждого разрабатываемого образца вооружения составляются оригинальные таблицы стрельбы, которые состоят из баллистических и поправочных таблиц. Ниже будет рассматриваться конкретный вариант наполнения таблиц, что не повлияет на общность результатов работы, но придаст материалу статьи большую стройность изложения.
В состав баллистических таблиц входят три таблицы:
г = ГТ(Б, Е); а = Та(Ц Е); р= ЩБ, Е), (1)
где т - двумерная таблица полетного времени снаряда в зависимости от полной дальности полета и угла места цели; а - двумерная таблица угла превышения в зависимости от полной дальности полета и угла места цели (вертикальная плоскость); в - двумерная таблица угла деривации в зависимости от полной дальности полета и угла места цели (горизонтальная плоскость).
Размер таблицы зависит от калибра артиллерийского орудия (для калибра 30 мм размер таблицы составляет около 200 значений). В связи с
компактностью таблиц их размещение в памяти вычислительной системы не вызывает трудностей. Выборка из двумерной таблицы, как правило, сопровождается линейной интерполяцией значений для произвольных аргументов дальности D и угла места E. Практическая реализация задачи встречи с использованием табличного прогнозирования траектории полета снаряда возможна, если производительность вычислительной системы превышает миллион операций в секунду. Для настоящего уровня развития вычислительной техники указанное требование по производительности вычислений является довольно скромным (например, микроконтроллер фирмы Atmel AT91R40807 производительностью 30 миллионов операций в секунду снят с производства как морально устаревший).
Альтернативой табличному способу расчета траектории полета снаряда является способ численного интегрирования системы дифференциальных уравнений, описывающих траекторию полета снаряда в атмосфере Земли. Задача построения системы дифференциальных уравнений, в основе которой лежит изучение влияния действия внешних сил и моментов сил на снаряд с момента прохождения его среза дула ствола орудия, называется основной задачей внешней баллистики. Этой задачей занимались довольно давно, теория вопроса достаточно проработана [1, 2]. В этой литературе представлена расчетная модель для вычисления кинематических элементов траектории полета снарядов, которая позволяет учесть форму и вращение Земли, возмущение атмосферы, принудительное вращение снаряда, колебания снаряда относительно центра масс. Практическая реализация задачи встречи с использованием метода численного интегрирования системы дифференциальных уравнений траектории полета снаряда возможна, если производительность вычислительной системы превышает 15 миллионов операций в секунду. Такая производительность вычислений вполне доступна на практике (например, микроконтроллер фирмы Atmel AT91R40008 производительностью 60 миллионов операций в секунду, AT91SAM9 - 200 миллионов операций в секунду, процессоры фирмы Texas Instruments TMS320C64x+ - 1000 миллионов операций в секунду). Таким образом, в настоящее время не существует объективных препятствий в практической реализации метода численного решения задачи внешней баллистики в вычислительной системе зенитного комплекса.
Хотя метод численного интегрирования по сравнению с табличным является более затратным, он превосходит последний в точности определения параметров траектории полета снаряда. Табличный метод представляет собой приближенное частное решение, в котором упрощенно учитывается взаимное влияние (иными словами, учитывается только главный член разложения) отдельных внешних факторов на результат решения задачи внешней баллистики снаряда. Это важно по той причине, что в табличном методе решение ищется в виде суммы: расчет траектории полета
снаряда в стандартных метеобаллистических условиях (баллистические таблицы) плюс изменения в траектории, вызванные отклонениями баллистических и метеорологических условий от принятых за стандартные (поправочные таблицы). Таким образом, точность табличного метода расчета ограничена учетом только линейных связей внешних факторов с результатом.
Система дифференциальных уравнений моделирует влияние физических процессов (имеющих место в атмосфере Земли) на формирование траектории полета снаряда. На вход задачи подаются не стандартные, а реальные метеобаллистические параметры, которые трансформируются в результат согласно алгоритму численного интегрирования, поэтому влияние внешних факторов на результат в этом методе учитывается наиболее полно. В дополнение к отмеченному, в системе дифференциальных уравнений учитывается влияние факторов, которых нет в табличном методе, например, влияние на параметры траектории снаряда вращения Земли вокруг своей оси. Фактор суточного вращения Земли приводит к учету изменения силы тяжести снаряда, учету зависимости силы тяжести от широты стояния орудия, к учету влияния ускорения Кориолиса на снаряд.
Опережающие темпы развития цифровой элементной базы по сравнению с темпами совершенствования материальной части артиллерийского вооружения создали условия повышения эффективности существующих зенитных комплексов за счет внедрения более совершенных алгоритмов управления, которые не удалось реализовать ранее из-за ограниченных возможностей вычислительной техники. В рамках предполагаемой модернизации предлагается заменить алгоритм табличного прогнозирования траектории полета снаряда на алгоритм численного интегрирования системы дифференциальных уравнений задачи внешней баллистики полета снаряда. Предметом настоящей работы являются исследование проблем, которые имеют место при реализации алгоритма интегрирования дифференциальных уравнений полета снаряда, при условии, что в качестве исходных данных принимаются известные баллистические таблицы артиллерийского орудия.
Исходная система дифференциальных уравнений траектории полета снаряда, в соответствии с [1, 2], представлена в фазовых координатах
{ V, О, ¥}, где V - производная земной скорости снаряда, О - производная
угла наклона траектории, ¥ - производная угла пути. Задача встречи снаряда с целью решается в земной прямоугольной системе координат ГОСТ 20058-80 [3](начало координат в точке орудия; ось ОбУб направлена вверх по местной вертикали; ось ОбХб параллельна линии топографической координатной сетки и направлена на север; ось Об/б дополняется до
получения правой системы координат). Для согласования систем координат предлагается преобразовать фазовые координаты { V, О, ХР} в фазовые координаты {\&сх,\&су,\&сг} стартовой прямоугольной системы координат ГОСТ 20058-80 (начало координат и ось ОбУс совпадают с началом и осью земной системы координат; ось ОбХс находится в плоскости стрельбы и направлена в сторону цели; ось Об/с дополняется до получения правой системы координат). Численное интегрирование системы дифференциальных уравнений траектории полета снаряда выполняется в стартовой прямоугольной системе координат. В результате интегрирования получают траекторию снаряда в фазовых координатах { Хс, Ус, Zc}. Для приведения результатов к земной прямоугольной системе координат { Хб, Уб, Zs} достаточно развернуть стартовую систему координат на величину азимута орудия. Так как баллистические таблицы представлены в стартовой системе координат, то и задачу согласования табличного и численного методов предлагается решать в стартовой системе координат.
Суть задачи согласования рассматриваемых методов заключается в том, чтобы рассчитать величины согласующих коэффициентов для метода численного интегрирования, при которых оба метода расчета дают одинаковые результаты с точностью до заданной погрешности. Задача согласования решается для неподвижной цели, неподвижного орудия в нормальных метеобаллистических условиях без учета фактора вращения Земли.
Входными величинами численного метода являются угол возвышения орудия О о и время полета снаряда т0. Входными величинами табличного метода являются угол места цели Ео и полная дальность полета снаряда Б о . Таблицы коэффициентов для согласования методов {1Х, ¡у, 12}
предлагается рассчитывать для каждой координаты { Хс, Ус, Zc }, при этом в качестве значений аргументов принимаются наборы величин О о и Б о .В качестве задающих значений дальностей Б о предлагается набор известных табличных значений (например, 2оо, 4оо, боо м, ...- баллистические таблицы представлены с шагом 200 м по дальности). В качестве задающих значений О о предлагается набор значений с фиксированным шагом (например, ©о = 0°,3°,6°,...). Значения коэффициентов согласования {1Х, 1у}
включаются в качестве множителей при расчете коэффициента сопротивления воздуха, а {} - коэффициента деривации.
Расчет коэффициентов согласования начинается с задания величин угла возвышения орудия О о и полной дальности полета снаряда Б о . Так
как угол возвышения О, угол места Е и угол превышения а связаны ме-
92
жду собой соотношением
0 = Е + а, (2)
то значение Е0 при известном значении 0 0 определяется из баллистической таблицы а = Та( £>0, Е). По рассчитанному значению Е0 и заданному значению В0 из баллистической таблицы т = Тт(£0, Е0) находится полетное время снаряда Т). По рассчитанному значению Е0 и заданному значению В0 из баллистической таблицы ¡ = Тр(£0,Е0) находится угол деривации ¡0. Таким образом, из данных баллистических таблиц следует, что снаряд, выпущенный под углом к горизонту 0 0, через время полета Т0 должен оказаться в точке { X0,У0, ^о}.
Х0 = В • Со5(Е0) • Со5(¡0); = В • Б1п(Е0);
70 = В0 • (Е0) • МД,)
В результате интегрирования системы дифференциальных уравнений полета снаряда, при тех же начальных условиях, получим точку с координатами { Хг, Уг, Т^г}. Значения коэффициентов согласования
{1Х, 1у, 1г} назначаются из множества допустимых значений [ 1п, 1т ], где 1п - минимальное допустимое значение коэффициента (например, 1 =0.5); 1т - максимальное допустимое значение коэффициента (например, 1 =1.9).
Расхождение результатов расчета траектории снаряда оценивается погрешностями
Ах = зЬб(Хг - X0); Ау = зЬб(Уг - У0); Az = зЬб(Zг - ). (4)
Фиксация величин коэффициентов согласования при условии минимума суммарной погрешности
А г = ^А* + Ау + А2) (5)
позволяет рассчитать их оптимальные значения {1°, 1°°, 1°}.
Контроль величин составляющих (4) и суммарной (5) погрешностей позволяет разработчику оценить степень соответствия табличного и численного способов расчета и обоснованно принять решение к дальнейшему использованию оптимальных коэффициентов согласования в решении задачи внешней баллистики способом численного интегрирования системы дифференциальных уравнений, описывающих траекторию полета артиллерийского снаряда в атмосфере Земли.
Список литературы
1. Дмитриевский А.А., Лысенко Л.Н. Внешняя баллистика: ечебник для студентов вузов. 4-е изд., перераб. и доп. М.: Машиностроение, 2005.
608 с.
2. Коновалов А. А., Николаев Ю.В. Внешняя баллистика. М.: ЦНИИ информации, 1979. 228 с.
3. ГОСТ 20058-80 Динамика летательных аппаратов в атмосфере. Термины, определения и обозначения.
Мельников Петр Николаевич, канд. техн. наук, вед. науч. сотр., peter@olvs.miee.ru, Россия, Москва, Зеленоград, Национальный исследовательский университет «МИЭТ»
Терпигорев Михаил Александрович, канд. техн. наук, ст. науч. сотр., terpm@olvs.miee.ги, Россия, Москва, Зеленоград, Национальный исследовательский университет «МИЭТ»
CALCULATIONS OF THE TRAJECTORY OF AN ARTILLERY PROJECTILE
P.N.Melnikov, M.A. Terpigorev
In the article the method of transition from the tabular method of calculating the trajectory of an artillery shell to the calculation of the trajectory of the projectile by the method of numerical integration of differential equations of tasks external ballistics.
Key words: artillery fire, calculating system.
Melnikov Peter Nikolaevich, candidate of technical science, leading researcher, peter@olvs.miee.ru, Russia, Moscow, Zelenograd, National Research University of Electronic Technology,
Terpigorev Micle Alexandrovich, candidate of technical science, senior researcher, terpm@olvs.miee.ru, Russia, Moscow, Zelenograd, National Research University of Electronic Technology
УДК 551.508.54
ОДНОПЛОСКОСТНОЙ ДАТЧИК СКОРОСТИ ВОЗДУШНОГО ПОТОКА
А.П. Панов, А.Н. Серов, И.Н. Мочегов
Рассмотрены основные методы измерения скорости ветра и предложена структурная схема одноплоскостного датчика скорости воздушного потока. Объяснены принципы функционирования датчика ветра. Приведены требования, предъявляемые к форме обтекателя. Представлены результаты моделирования и натурного испытания разработанного устройства.
Ключевые слова: дифференциальный датчик давления, анемометр, тензомо-
дуль.
На точность работы системы управления объектом специального назначения оказывает влияние множество факторов [1]. Одним из основ-
