Понятие и типы баллистики

Именно для быстрого составления баллистических таблиц по заказу армии США в Лаборатории баллистических исследований. в 1946 г. была создана первая вычислительная машина «Эниак» (ENIAC - Electronic Number Integrator And Computer - Электронный цифровой интегратор и вычислитель). Вычисления на «Эниаке» велись в десятичной системе, а для изменения программы требовалось установить в определенное положение тысячи переключателей и соединить сотни кабелей, и поэтому в среднем на подготовку машины к вычислению одной таблицы уходило два дня кропотливой ручной работы .

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

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

Общая формула для силы сопротивления воздуха

На рисунке 1 показана схема сил, приложенных к снаряду на траектории: - сила тяжести; - равнодействующая аэродинамических сил, т.е. сила сопротивления воздуха. Она приложена в центре давления С, не совпадающем с центром массы O. Расстояние между этими точками определяется по формуле Гобара [2]. На рисунке д - угол атаки, т.е. угол между осью снаряда и касательной к траектории в данной точке (на касательной лежит вектор скорости снаряда ); - угол между вектором скорости и горизонтом. Если силу перенести в центр масс О и одновременно приложить к этой точке уравновешивающую силу (), то возникает пара сил, создающая опрокидывающий момент (его учитывают при исследовании движения снаряда как твердого тела). Силу , приложенную в центре масс, раскладывают на две составляющие: - сила лобового сопротивления (она лежит на касательной к траектории и направлена в сторону, обратную по отношению к вектору скорости) и - подъемная сила. В дальнейшем рассматриваем упрощенную схему приложения указанных сил, полагая и считая, что сила направлена по оси снаряда; в этом случае и .

Силы, приложенные к снаряду на траектории

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

где - плотность воздуха; - площадь миделевого сечения снаряда (d - калибр); - скоростной напор; - коэффициент лобового сопротивления; - число Маха; a - скорость звука в данной точке траектории; - число Рейнольдса; - кинематический коэффициент вязкости.

Зависимость определяют опытным путем для снарядов типовой («эталонной») формы. Подобие процессов обтекания снарядов воздушным потоком чаще всего не обеспечивается из-за различия конфигурации головной части, и, для того чтобы иметь возможность использовать имеющиеся опытные данные, вводят коэффициент формы снаряда

,

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

Преобразуя формулу , получаем выражение для «ускорения сопротивления воздуха»

где q- масса снаряда. Далее вводят функцию

,

где y - высота; - плотность воздуха на поверхности Земли в точке выстрела. Кроме того, для получения более удобных для практических расчетов значений соответствующих величин, вводят множитель

где кг/м3 - плотность воздуха для нормальных артиллерийских условий. Тогда (2), с учетом (1), будет иметь вид

В этом выражении фигурирует баллистический коэффициент

..

Обычно принимают, что , т.е.

.

Здесь введен коэффициент относительной массы («поперечная нагрузка») , где d - в дециметрах. Видно, что баллистический коэффициент (а, следовательно, и сила сопротивления воздуха) изменяется обратно пропорционально калибру.

Функция

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

.

Как известно [2], скорость звука

,

где - показатель адиабаты для воздуха, обычно принимаемый равным 1,4; R - универсальная газовая постоянная; - «виртуальная» температура, учитывающая влажность воздуха; T - абсолютная температура; e - давление водяного пара; h - давление влажного воздуха/

Эталонные законы сопротивления воздуха приведены к скорости звука в нормальных условиях м/с, поэтому аргумент преобразуют:

.

Здесь

-

так называемая виртуальная скорость. Таким образом,

Соответственно,

Где

Зависимость обычно задают по [3]:

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

2 Законы сопротивления воздуха

Л. Эйлер при решении задачи о полете снаряда пользовался функцией , установленной Ньютоном и применяемой в основном для дозвуковых скоростей. Одной из первой была степенная функция Маиевского-Забудского [3]

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

Базируясь на работах Маиевского-Забудского и опытах конца XIX века, итальянский баллистик Франческо Сиаччи предложил новую функцию сопротивления воздуха , носящую его имя (1888 г.). Сиаччи также принял за эталонный снаряд старой формы, но сгладил угловые точки на графике . Большой заслугой Сиаччи является эмпирическая аппроксимация закона сопротивления воздуха в предложенной им форме (5) [3]:

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

С развитием артиллерии основным становится современный снаряд дальнобойной формы, имеющий удлиненную головную часть и сравнительно короткую хвостовую часть. Опыты по созданию новой функции проводились после Первой мировой войны в ряде стран, например, в 1921-1923 гг. во Франции (законы Гарнье и Дюпуи).

В нашей стране был создан закон сопротивления воздуха 1930 г. На его основе составлены таблицы внешней баллистики АНИИ, однако выяснилось, что данный закон дает неточные результаты при расчете траектории с большими начальными скоростями; кроме того, коэффициент формы современных снарядов по отношению к функции 1930 г. заметно колеблется при различных скоростях.

3 Закон 1943 г.

Перед Великой Отечественной войной в СССР были начаты работы по установлению новой функции сопротивления воздуха на основе обработки результатов стрельб современными снарядами дальнобойной формы. Эти работы были закончены в 1943 г., новая функция получила название закон Артиллерийской академии им. Ф.Э. Дзержинского, или просто закон 1943 г. При этом. была обнаружена ошибка функции Сиаччи, проявляющаяся при скорости снаряда более 1410 м/с. Закон 1943 г. принят в нашей стране в качестве основного. Применительно к этой функции проводятся все баллистические расчеты, хотя ввиду наличия таблиц находят применение также функции 1930 г. и Сиаччи.

Полная таблица закона 1943 г. содержится в книге [4]; в сокращенном виде она приведена в [5], наряду с законами Сиаччи и 1930 г. В [2] функция 1943 г. задана в пределах ограниченного диапазона (), разбитого на участки:

В [6] имеется следующее описание закона 1943 г.:

Внутри этого диапазона рекомендуется производить пересчет с закона Сиаччи, значения коэффициента согласования (переходного множителя)приведены в таблице 1:

Таблица - Переходный множитель [6]

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

Коэффициенты формы для современных снарядов (ОФ) по отношению к закону 1943 г. изменяются в пределах , а по отношению к функции Сиаччи .

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

4 Аппроксимация закона 1943 г.

Можно предложить такой способ коррекции закона Сиаччи и приведения его к закону 1943 г. [7]. Определив по табличным данным (дискретную) зависимость коэффициента формы от скорости, затем аппроксимировать ее как некую непрерывную функцию скорости и затем производить пересчет следующим образом:

.

Результаты реализации этой идеи в среде пакета MathCAD представлены на рисунке 2, где 1 - ; 2 - табличный закон 1943 г, ; 3 - функция согласования ; 4- аппроксимация функции согласования ; 5 - .

Функция согласования аппроксимируется полиномом 3-го порядка:

,

коэффициенты которого определены с помощью функции MathCADlinfit, относящейся к линейной комбинации аппроксимирующих формул:

; ; ; .

Аппроксимация закона 1943 г.:

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

Таким образом, принимаем следующее эмпирическое описание закона 1943 г.:

В таблице 2 дано сравнение данных, приведенных в работе [5], с получаемыми по предлагаемой аппроксимации: 1 - табличные значения ; 2 - расчет по данной методике; 3 - отклонение, %.

Таблица - Сравнение аппроксимирующих и табличных значений

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

5 Пример расчета

Применение предложенной аппроксимации проиллюстрируем на примере расчета параметров траектории снаряда линкора «Бисмарк», который был проведен автором при математическом моделировании обстрела английского линейного крейсера «Худ» 24 мая 1941 г. Подробное описание «дуэли» двух выдающихся кораблей приведено в [8].
В работе [9] читаем: «…коэффициент формы iследует рассматривать как параметр, позволяющий согласовывать результаты теоретических расчетов с опытными данными. Например, пусть на основании стрельб снарядами определенного типа при фиксированных значениях начальной скорости и угла бросания найдена опытная дальность стрельбы Х. …По величинам Х, и можно определить коэффициент формы снаряда i. Если расчет траектории проводить с использованием коэффициента , удовлетворяющего выражению при тех же значениях и , то получим дальность стрельбы, совпадающую с опытной. Этот способ применяется для определения коэффициента формы при составлении таблиц стрельбы для конкретного орудия».

Соответствующие вычисления проводят, используя известную систему уравнений, описывающих движение снаряда как материальной точки [2]:

где - угол наклона касательной к траектории (вектора скорости) относительно горизонта.

Именно так был определен коэффициент i для снарядов «Бисмарка», На линкоре были установлены восемь 380-мм орудий (по два в каждой из четырех башен) 38cm/52 SKC/34. Известно, что максимальная дальность 35 550 м достигается при массе снаряда 800 кг, дульной скорости 820 м/с и угле возвышения [10]. Методом подбора, пользуясь соответствующей программой численного решения системы (9), было определено и .

В работе [10] приведены параметры траектории при стрельбе с различными углами возвышения; в таблице 3 дается сравнение этих данных с результатами расчета, полученными с помощью закона Сиаччи при (в знаменателе). Расхождение между этими данными составляет единицы и доли процента. Закон Сиаччи был использован, поскольку расчеты, проведенные немцами, могли быть получены только с помощью именно этого закона. Это подтверждается сведениями в статье [11], в которой представлены результаты расчетов внешней баллистики «Бисмарка», выполненные в 1939-1940 гг. во время достройки линкора на судоверфи «Блом и Фосс» [12]. Данные результаты также представлены и в графической форме на Интернет-сайте линейного крейсера «Худ» [13].

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

Считая приемлемость принятой методики определения коэффициентов формы подтвержденной, можно провести расчет параметров траектории снаряда «Бисмарка» при обстреле «Худа», дальность составляла м [87]. В таблице 3 приведены значения скорости падения снаряда , угла падения и времени полета , рассчитанные по закону Сиаччи и предложенной аппроксимации закона 1943 г.

Таблица - Параметры траектории снаряда линкора «Бисмарк»

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

4.6 Как уменьшить сопротивление воздуха

Воздух тормозит летящий снаряд, замедляет его полет.

Можно ли бороться с этим замедлением?

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

А на войне важно иметь возможность забросить снаряд как можно дальше. Поэтому уменьшать его скорость не всегда выгодно.

Поищем, нет ли других, более подходящих способов бороться с замедлением полета снаряда из-за сопротивления воздуха.

Такие способы существуют.

Представьте себе, что вы хотите выбраться из трамвая, битком набитого пассажирами. Попробуйте итти прямо -- грудью вперед; пожалуй, вы не доберётесь до выхода. Но если вы начнете пробираться боком, вам уже не так трудно будет протолкнуться.

Нечто подобное испытывает и снаряд в полете: не безразлично, как он будет пробираться между частицами воздуха.

Был в старину -- во время первой севастопольской обороны -- такой снаряд: светящее ядро к полупудовой медной мортире. Это ядро имело форму цилиндра.

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

Такое ядро летело всего лишь метров на 500.

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

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

Вот почему головную часть снаряда начали заострять, едва лишь научились делать устойчивым в полете продолговатый снаряд,-- еще в середине XIX века.

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

В XX веке резко возросли скорости транспорта всех видов, быстро развилась авиация. Ученые начали внимательно изучать действие сопротивления воздуха на быстро движущиеся предметы разной формы. Оказалось, что не только для самолета, но даже для быстроходного автомобиля или поезда важна такая форма, которая является удобообтекаемой. Если автомобилю придать такую форму, то при большой скорости движения он экономит 10-15 процентов горючего или при том же расходе горючего движется заметно быстрее.

Тем большее значение имеет форма снаряда: ведь снаряд движется во много раз быстрее автомобиля, он встречает огромное сопротивление воздуха. Взгляните на все предыдущие рисунки с раздела. Перед вами четыре снаряда разной формы. На рисунках изображены волны и завихрения воздуха, которые сопровождали бы полет каждого из этих снарядов, если бы скорость их всех была одна и та же, и притом больше, чем скорость звука. Давление на головную часть снаряда тем меньше, чем она острее. Разреженная зона за снарядам также тем меньше, чем больше скошена его донная часть; меньше в этом случае и завихрений позади летящего снаряда.

Очевидно, что наиболее выгодна форма снаряда, изображенная на правом рисунке.

Более подробное изучение этого вопроса показало, что каждой скорости полета соответствует своя наиболее выгодная форма снаряда.

Чем больше скорость снаряда, тем острее должна быть его головная часть.

Допустим, что воздух давит на головную часть снаряда с силой 4 атмосферы, а в разреженной зоне позади снаряда давление составляет всего лишь четверть атмосферы.

Давление на дно снаряда уменьшилось против нормального на три четверти атмосферы: это составляет примерно пятую часть того давления, которое испытывает головная часть снаряда.

А вот другой снаряд: скорость его значительно больше, чем у первого, а потому он испытывает и большее сопротивление воздуха,-- предположим, равное давлению 100 атмосфер. Пусть он летит так быстро, что за ним позади образуется почти полная пустота: частицы воздуха не успевают ее заполнить. Разница с нормальным давлением составляет целую атмосферу.

Но это всего лишь 1 процент -- всего сотая часть -- того давления, которое испытывает головная часть такого снаряда!

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

30 лет тому назад граната 76-миллиметровой пушки могла пролететь около 8,5 километра.

Но стоило только заострить ее головную часть, удлинить и скосить длинную часть, как граната такого же веса стала лететь больше чем на 11 километров; простое изменение формы снаряда увеличило почти на одну треть дальность его полета.

Размещено на Allbest.ru