Теория вероятностей в артиллерии. Часть 1. Разрушение одиночной ДЗОТ

Введение

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

Но в погоне за точностью, дальнобойностью и пробивным действием не следует забывать и об экономности: от трёх первых факторов не будет никакого толка, если снаряды требуют экзотичных и редких материалов для своего изготовления, а орудия можно будет сделать только с привлечением самых что ни на есть высококвалифицированных рабочих в количестве пары штук в месяц. Здесь в игру вступают инженеры-технологи, которые по мере возможностей стараются объединить идеи учёных и конструкторов с суровой реальностью производственных возможностей. В итоге применительно к артиллерийской системе и боеприпасам для неё получается некоторый компромисс между изначальным проектом, и теми орудиями и снарядами, которые поставляются в войска. Причём далеко не всегда получается так, чтобы создать изобилие боеприпасов и не оглядываться на их расход на учениях и в бою. Здесь скорее имеет место обратная ситуация, когда этот расход нужно тщательно планировать, дабы сберечь не такой уж бесконечный запас снарядов и не беспредельный ресурс орудий. Не случайно именно это весьма настойчиво рекомендуется в наставлениях, руководствах службы и таблицах стрельбы. Поскольку на движение снаряда влияет целый сонм случайных изменений многих параметров, детерминистические методы расчёта дают только усреднённый результат, а на практике реальные значения величин имеют некоторый разброс около полученных из теории значений. Эти отклонения применительно к координатам точки падения снаряда называются рассеянием. Из-за него даже в случае максимально точной наводки орудия поражение цели вовсе не является гарантированным результатом. Как следствие, возникает необходимость оценки количества снарядов, требуемых для выполнения учебной или боевой задачи. Тут на помощь артиллеристу вновь приходит математика, точнее её раздел, называемый теорией вероятности — учение о случайных событиях и закономерностях в их мире. В данной статье автор как раз и предлагает рассмотреть пример её использования на практическом примере, который является частью профессиональной подготовки артиллеристов.

Постановка задачи

Рассмотрим ту же мысленную учебно-боевую задачу, которая обсуждалась в предыдущей статье «Что скрывается за внешней простотой» — 152-мм гаубица обр. 1943 г. (Д-1) обстреливает две дерево-земляных огневых точки (ДЗОТ) условного противника. Первая из них (цель №1) удалён на 2500 метров, вторая (цель №2) — на 5000 метров. Будем считать, что наводка орудия с учётом угла места, деривации, метеорологических условий, состояния артиллерийской системы и боеприпасов уже точно выполнена — то есть центр эллипса рассеивания координат точек разрывов совмещён с целью. Для ведения огня назначена осколочно-фугасная граната ОФ-530 с установкой взрывателя на замедленное (фугасное) действие, метательный заряд можно использовать любой. Необходимо оценить расход гранат, требуемых для разрушения обеих ДЗОТ.

Очевидно, что для решения поставленной задачи нужна и некоторая информация о цели и разрушительных свойствах используемых боеприпасов. Пусть обитаемое пространство вражеской ДЗОТ в плане представляет собой прямоугольник со сторонами 3×4 м. Также известно, что при взрыве гранаты ОФ-530 в грунте средней плотности образуется воронка диаметром 3,5 м и глубиной 1,2 м. Поэтому для разрушения ДЗОТ в общем случае прямое попадание в него не требуется. Будем считать, что если точка разрыва удалена не далее чем на 1 метр по каждой из своих координат в горизонтальной плоскости от любой точки обитаемого пространства ДЗОТ, последняя гарантированно разрушается совместным действием механического воздействия, ударной волны и осколков. Иными словами, граната должна попасть тоже в прямоугольную область размером 5×6 м, центр которой совпадает с центром обитаемого пространства ДЗОТ.

Определения и допущения

Строгие математические определения нужных нам для решения поставленной задачи понятий и величин можно найти в любом учебнике по теории вероятностей. На субъективный взгляд автора, в этом вопросе ничто не заменит замечательную книгу Елены Сергеевны Вентцель «Теория вероятностей», выдержавшую целый ряд переизданий, да ещё и написанную по материалам лекций по теории стрельбы в одном из военных высших учебных заведений. Здесь же поневоле для краткости изложения придётся воспользоваться некоторыми упрощениями и приближениями, которые, тем не менее, дадут достаточно точный результат.

Математические величины

Введём следующие обозначения:

P₁ — вероятность разрушения ДЗОТ одним выстрелом. Её можно определить как отношение числа попавших в определённую выше окрестность ДЗОТ гранат к общему числу израсходованных боеприпасов N, если условия каждого выстрела и полёта снаряда совершенно идентичны, а N стремится к бесконечности. В реальности, конечно, количество выстрелов ограничено и это отношение является так называемой относительной частотой требуемого нам события (попадания гранаты в заданную область и разрушения ДЗОТ как следствие). Но уже при N в диапазоне 10—20 отличие вероятности от относительной частоты уже не столь велико; а при N порядка 100 им уже можно пренебречь для любых практических расчётов.

P_i — вероятность того, что ДЗОТ окажется разрушенным после i = 2, 3, 4, … выстрелов; при условии, что все выпущенные ранее гранаты — 1-я, … , (i-1)-я — не попали в заданную окрестность цели. Это является следствием того факта, что для разрушения ДЗОТ достаточно только одного близкого попадания и после него продолжать стрельбу просто незачем. Соответственно математическое выражение для P_i для этого условия согласно правилам сложения и умножения вероятностей независимых друг от друга событий будет следующим:

P_i = P₁ + (1 - P₁) × P₁ + (1 - P₁)² × P₁ + … + (1 - P₁)^{i - 1} × P₁

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

F(x, y) — двумерная плотность вероятности разрушения ДЗОТ одним выстрелом. Это действительная функция от координат x и y в горизонтальной плоскости, обладающая тем свойством, что при достаточно малых размерах прямоугольной области Δx × Δy в окрестности точки с координатами (x₀, y₀) вероятность P₁ может быть вычислена по формуле:

P₁ = F(x₀, y₀) × Δx × Δy

Критерием малости Δx и Δy можно считать небольшое изменение функции F (например менее 1%) при любом изменении её аргументов в диапазоне

(x₀ - ½Δx … x₀ + ½Δx, y₀ - ½Δy … y₀ + ½Δy)

Более строгое определение плотности вероятности уже требует знания высшей математики, поскольку P₁ будет ни чем иным, как двумерным интегралом от F(x, y) по требуемой нам окрестности. Однако для наших целей вполне пригодится и сильно упрощённый его вариант.

Алгоритм решения и используемые приближения

В артиллерийском деле принято считать функцию F(x, y) подчиняющимся так называемому нормальному или гауссовому (в честь исследовавшего подобные явления знаменитого немецкого математика Карла Гаусса) закону распределения случайной величины — когда отклонение от её среднего значения вызывается совместным действием многих и независимых друг от друга факторов. Математически она записывается следующим образом:

F(x, y) = exp(-x² / 2В_Д) × exp(-y² / 2В_б) / (2 × π × В_Д × В_б) ,

где начало координат совмещено с целью; ось X считается направленной вдоль проекции линии цели на горизонтальную плоскость уровня цели по высоте, а ось Y перпендикулярна оси X и лежит в той же плоскости. Величины В_Д и В_б называются срединными отклонениями по дальности и фронту соответственно. Под exp(z) понимается показательная функция с возводимым в степень z основанием натурального логарифма e ≈ 2,7182…

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

Для практического применения задание F(x, y) в аналитическом виде подходит плохо. Только ЭВМ способна быстро рассчитать нам вероятность по известным Δx и Δy, В_Д и В_б (последние приведены в таблицах стрельбы); даже с помощью инженерного калькулятора это займёт достаточно много времени с высокой вероятностью ошибки при выполнении операций под действием стресса. Однако нормальное распределение обладает рядом замечательных свойств, строго математический вывод которых превосходно изложен в учебнике Е. С. Вентцель. Нам будет полезно следующее из них: при большом числе выстрелов 50% попаданий придутся в эллипс с полуосью В_Д вдоль оси X и полуосью В_б вдоль оси Y. Если увеличить полуоси этого эллипса в два раза, то внутри него очутится 82% всех попаданий, а если в три — то 96%. Оставшимися 4% на практике обычно пренебрегают, причём львиная их доля из них уложится в ещё раз расширенный эллипс с полуосями в 4В_Д и 4В_б.

Без особого ущерба для точности в каждом из интервалов с шагом в одно срединное отклонение по обеим координатам функцию F(x, y) можно считать постоянной, равной некому среднему значению, которое даёт заданную вероятность при перемножении его на длину интервала. Таким образом, её плавные колоколообразные сечения в координатных плоскостях (х меняется, y = 0 и наоборот) подменяются чем-то вроде ступенчатого пьедестала. Это существеннейшим образом упрощает счёт, позволяя обойтись только тривиальными арифметическими действиями.