Научная статья на тему 'Математическое моделирование движения снаряда в периоде последействия пороховых газов при выстреле из артиллерийского орудия'

Математическое моделирование движения снаряда в периоде последействия пороховых газов при выстреле из артиллерийского орудия Текст научной статьи по специальности «Механика и машиностроение»

CC BY
632
87
i Надоели баннеры? Вы всегда можете отключить рекламу.
i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Математическое моделирование движения снаряда в периоде последействия пороховых газов при выстреле из артиллерийского орудия»

Практическое выполнение расчетов по формуле Очевидно, что предлагаемый порядок расчета ни

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

ЛИТЕРАТУРА

1. Федеральный закон «О пожарной безопасности» от 21 декабря 1994 г. № 69-ФЗ). М., 2006 (с изменениями и дополнениями).

2. Постановление Правительства РФ от 7 апреля 2009 г. N 304 «Об утверждении правил оценки соответствия объектов защиты (продукции) установленным требованиям пожарной безопасности путем независимой оценки пожарного риска (в ред. Постановления Правительства РФ от 02.10.2009 N 777).

3. Приказ МЧС РФ от 30 июня 2009 г. N 382 «Об утверждении методики определения расчетных величин пожарного риска в зданиях, сооружениях и строениях различных классов функциональной пожарной опасности (в ред. Приказа МЧС РФ от 12.12.2011 N 749).

4. Приказ МЧС РФ от 16 марта 2007 г. № 141 «Об утверждении Инструкции о порядке согласования отступлений от требований пожарной безопасности, а также не установленных нормативными документами дополнительных требований пожарной безопасности» (с изменениями от 7 февраля 2008 г.).

5. Абдурагимов И. М. Еще раз о принципиальной невозможности выполнения расчетов пожарных рисков детерминированными методами.// Пожаровзрывобезопасность. -2013.- Т. 22, № 6. - С. 13-23.

6. Федеральный закон Российской Федерации от 22 июля 2008 г. N 123-ФЗ "Технический регламент о требованиях пожарной безопасности" (с изменениями на 13 июля 2015 года).

7. Светушенко С.Г. Аудит пожарной безопасности. Специальные технические условия и расчет пожарного риска. Три сомнительных кита пожарной безопасности. "Алгоритм Безопасности" № 5, 2011 год, с. 72-76.

8. Правила пожарной безопасности в Российской Федерации. ППБ 01-03.

9. ГОСТ 12.1.004-91 ССБТ. Пожарная безопасность. Общие требования (с Изменением N 1).

10. Пожары и пожарная безопасность в 2009 г. : стат. Сб. / Под общей редакцией Н.П. Копылова. -М.: ФГУ ВНИИПО, 2010. -137 с.

УДК 623.412 Букаси Амин

Филиал ФГКВОУ ВО «Военная академия материально-технического обеспечения имени генерала армии А.В. Хрулева» в Пензе, Алжир, Алжир

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДВИЖЕНИЯ СНАРЯДА В ПЕРИОДЕ ПОСЛЕДЕЙСТВИЯ ПОРОХОВЫХ ГАЗОВ ПРИ ВЫСТРЕЛЕ ИЗ АРТИЛЛЕРИЙСКОГО ОРУДИЯ

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

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

— снаряд рассматривается в виде тела конечной геометрии и массы, состоящего из активной части и ведущего устройства;

— снаряд рассматривается как твёрдое тело со смещенным (в общем случае) положением центра масс относительно геометрического центра;

— начальное (на дульном срезе ствола) положение снаряда в периоде последействия пороховых газов определяется параметрами поведения (динамики) системы «пушка-снаряд» при выстреле в момент потери им механической связи со стволом;

— движение снаряда рассматривается относительно невращающейся Земли.

Для анализа боевой эффективности танкового - снаряд рассматривается в виде тела конечной

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

ную задачу с учётом максимального числа факто- - снаряд рассматривается как твёрдое тело со

ров. смещенным (в общем случае) положением центра

Задача анализа боевой эффективности ТВ может масс относительно геометрического центра; быть решена только в случае наличия полной и - начальное (на дульном срезе ствола) поло-

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

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

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

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

математическая модель движения снаряда на стемы координат). Движение снаряда будем рас-внешнебаллистическом участке траектории; сматривать относительно не вращающейся Земли.

математическая модель действия снаряда по При выводе уравнений движения снаряда исполь-

цели. зованы правые системы координат (рисунок 1):

В основу построения системы уравнений мате- - инерциалвная система OgXgYgZg с единичными

матической модели движения снаряда в периоде последействия пороховых газов при выстреле из тан-

зекторами , образующими базис [ eg ] , и

ковой пушки положены следующие допущения: началом в некоторой фиксированной точке про-

странства, где ось ОдХд направлена по геометрической оси канала ствола в направлении стрельбы, ось ОgYg направлена перпендикулярно геометрической оси канала ствола вверх, а ось ОgZg составляет с осями ОдХд и ОдУд правую систему координат;

- вязанную со снарядом ОХУЕ с триадой единичны:': векторов 1,],К образующие: базис { е }, и началом в некоторой фиксированной точке, снаряда (геометрическом центре), в общем случае не совпадающем с её центром масс; ось ОХ направлена к вершине снаряда параллельно его геометрической оси симметрии.

Мгновенное положение связанной со снарядом система координат ОXYZ по отношению к инерциаль-ной ОgXgYgZg определяется пятью независимыми координатами: тремя координатами Х^^ положения начала системы координат ОXYZ, образующими радиус-вектор

и двумя углами Эйлера-Крылова Ц/,в , (углами, соответственно, тангажа и рысканья). Единичные векторы систем координат ОgXgYgZg и ОXYZ связаны линейными соотношениями преобразования в виде

[е„ ] = [а] [е] , (2)

где

[] -

матрица направляющих косинусов

а11 а12 а13

[а] = а21 а22 а23

а31 а32 а33 _

где (1,3 = 1,2,3) являются направляющими ко-

синусов углов Эйлера-Крылова:

аи = соъ^соъв ап=— ъту ап = зт^созв

а21 = со505т^ а22 = соъу а23 = вт^тб (3)

а31 = —

а = 0

= соъв

г =Х,

-У-,

■'в

■ г.

г

(1)

Углы ф , Д будем считать фиксированными и

определяемыми положением снаряда в момент потери им механической связи со стволом.

Рисунок 1 - Взаимное расположение систем координат, принятых в математической модели

В соответствии с формулами (2) и (3) единичные векторы систем координат ОдХдУг^д и 0±Х±У±И± связаны линейными соотношениями преобразования в общем виде

[е„] = [а] [а,][е]

( 4 )

обозначим через

р = А1 + ¥]+гк

(5)

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

Тогда радиус-вектор этой же частицы в инер-циальной системе координат можно представить в виде

К = г + р . (6)

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

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

,в . обозначим скорость изменения эти:: углов

через, соответственно, ЦУ, 0 (здесь и в дальнейшем точка над параметром обозначает производную этого параметра по времени). Поэтому, если принять, что - вектор) мгновенной скорости

вращения связанной со снарядом системы координат относительно инерциальной системы, то

О = Ц/+6 (7 )

Тогда, проектируя данное векторное равенство на оси системы координат, связанной со снарядом, получим

о = ад7+о^/ + п_к (в)

= -у/ыпв',

= О',

О2 = \yCOS0.

(9)

Дифференцируя равенство (8) по времени, находим угловое ускорение связанной со снарядом системы координат в виде

О = Од- /' + О + О.- к ,

(10)

а

Ох = - \//всо&в - \//&тв',

пу =в;

Ог = у/ сое в -цг в эт в.

(11)

Радиус-вектор произвольной элементарной частицы снаряда в момент времени Ь в инерциалвной

системе координат определяется вектором Я , дифференцируя который по времени, находим ско-роств и ускорение элементарной частицы снаряда:

Я = г + О.Х р ; Д = ? + Йх />+йх(Йх р).

(12)

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

Уравнения снаряда в векторной форме имеют следующий вид:

- уравнение поступателвного движения снаряда

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

п

(13)

1=1

где Рв - главный вектор внешних сил, действующих "г„ - вектор ускорения центра масс

на снаряд; г¥в снаряда, равный

(14)

т - масса снаряда.

- уравнение вращателвного движения снаряда

©хЙ + Йх©Й=7Йгг , (15)

где

©

тензор инерции снаряда;

М

в

главный

момент внешних сил, действующих на снаряд.

Принимая во внимание выражения (12), уравнения (13), (15) принимают соответственно

т(г + Г2хг+Г2х р + С1х (С1 х р) = Рв;

(16)

©•Й + Пх © С1=М

в-

На выходе снаряда из области пороховых газов образовавшейся в периоде последействия на него действует распределённая по поверхности тела сила сопротивления воздуха, т.е. к каждому элементу снаряда с1<Ув приложена элементарная сила

ёв=-Рв"с1еТв' (17)

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

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

воздуха р

и

вектор суммарного момента Мп

от-

носительно начала связанной со снарядом системы координат, математические выражения которых имеют вид:

Ра=-\р&с1<тв-, (18)

Ма

{а = - J (£ х и)рвс1(тв

где С

площадь поверхности снаряда; рв -

(19)

дав

ление воздуха, определяемое по законам аэродинамики.

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

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

1<а = -1<х + 1<у + 1<\2

I J к

сила лобового сопротивления;

(20)

аэроди-

аэродинамическая подъёмная сила;

намическая боковая сила.

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

Рх = С^;

РУ = С^Фт; (21)

Р = С2^т,

коэффициент силы лобового сопротивле-коэффициент аэродинамической подъем-

где Сха -

ния; Суа

ной силы;

С

коэффициент аэродинамической бо-

ковой силы;

ч = -

2

скоростной напор;

Я-

площадв Миделевого сечения снаряда,

МаС =масх7 +Ма:к ; (22)

Стандартная форма записи аэродинамических моментов определяется выражением.

Мас = тхдБтЬ;

МУс = тудБтЬ;

(23)

где

М'с = т^тЦ

коэффициенты соответствующих мо-

ментов.

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

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

В общем случае сила, с которой пороховые газы действуют на поверхность снаряда С , равна сумме элементарных сил dFp , действующих на элементарные площадки с1 СУ составляющие эту поверхности. Если принятв за п единичный вектор внешней нормали к поверхности с1 СУ , то сила р = —рГпс1з, где рг - давление пороховых газов на площадке (1СУ .

Суммируя силы р , получим выражения главного

относительно

вектора Р и главного момента Мр

начала связанной со снарядом системы координат силы давления пороховых газов в виде.

Рр =-^ргМ(Т ; (24)

X

Мр =-^(1-п)ргс!(Т , (25)

х

где I = XI -(У + Ур+ЯцХ)] ~(2 + 2р+аъ1Х)к - плечо

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

Момент в векторной форме определяются разложением на соответствующие оси связанной системы координат снаряда:

7 I пу 7 , п?Т.

О' =Ох1 +Оу) +а2к ; (26)

Ма =Мах1 +Мау7 +Ма*к . (27)

тарная сила.

где

где

т , т , т

Проекция силы тяжести и возникающий момент действующей на снаряд в системе ОgXgYgZg определяются из выражений:

(28)

(29)

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

Обобщая полученные выше выражения внешних сил, действующих на снаряд, для определения

главного вектора Рш всех внешних сил, входящих в правые части уравнений движения снаряда (16), получим

(30)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

относительно друг снаряда, но и другие параметры, определяющие их поведение с момента нарушения механического контакта со стволом, без знания которых невозможно корректное определение его дальнейшего движения, а также оценка в дальнейшем характеристик рассеивания снаряда. В частности, значение дульной скорости снаряда рассчитывается в соответствии с теоремой о количестве движения обозначив через Уд = {Ух,У¥,Уг} вектор абсолютной скорости центра масс снаряда, проекции вектора скорости БПС на оси связанной с ним системы координат определяются выражениями:

где Ух, УТ, У2

Р№=Р + Р+0;

=Мас +мрс +МСс .

(31)

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

В процессе решения системы дифференциальных уравнений движения снаряда в периоде последействия пороховых газов определяются не только координаты пространственного перемещения и скорость элементов БПС активной части и секторов ВУ

Ух=Уд-1 ; УГ=УД-] -,У2=Уд-к ; (32)

что позволяет рассчитывать результирующую скорость (V) движения снаряда в виде зависимости

V =у/у \ + V \ + V \ . (33)

Определение значений (начальных) углов вылета и бросания снаряда основано на зависимостях.

С учётом проведенных исследований, индивидуальный угол вылета снаряда отдельно для вертикальной и горизонтальной плоскостей в периоде последействия пороховых газов вычисляется суммированием значений угла вылета ( Ц/в , 0 ) и по-

правки угла ( А^ , А© ), последействия.

возникающей в периоде

А.

Комочков

А.

ЛИТЕРАТУРА

Баллистика ракетного и ствольного оружия. Волгоград 2010.

под ред. Л.Н.

1. Королева А. 470 с.

2. Баллистика ствольных систем / В. В. Бурлов, В. В. Грабин, А. Ю. Козлов [и др. Лысенко, А.М. Липанова. - М.: Машиностроение, 2006. - 460 с,

3. Ассовский И.Г., Кудрявцев О.А., Расходов В.С. Современные модели и методы расчета внутрибал-листических процессов в ракетно— ствольных системах. - Тула: ТулГУ, 2004. - 96 с.

УДК 623.412 Букаси Амин

Филиал ФГКВОУ ВО «Военная академия материально-технического обеспечения имени генерала армии А.В. Хрулева» в Пензе, Алжир, Алжир

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ФУНКЦИОНИРОВАНИЯ СИСТЕМЫ «ТАНКОВАЯ ПУШКА - СНАРЯД - ЦЕЛЬ»

Поражение цели бронебойно-подкалиберным снарядом представляет собой выходную характеристику сложного динамического процесса, основными этапами которого являются: - период взаимодействия снаряда со стволом в процессе выстрела; период движения снаряда в спутном потоке пороховых газов (период последействия); период движения снаряда на внешнетраекторном участке; процесс поражения цели.

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

В целях алгоритмизации процесса поражения цели бронебойно-подкалиберным снарядом определена структура основных блоков модели функционирования системы, «танковая — пушка - снаряд — цель».

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

Поражение цели БПС представляет собой выходную характеристику сложного динамического процесса, основными этапами которого являются:

- период взаимодействия снаряда со стволом в процессе выстрела;

- период движения снаряда в спутном потоке ПГ (период последействия);

- период движения снаряда на внешнетраектор-ном участке;

- процесс поражения цели.

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

В целях алгоритмизации процесса поражения цели БПС необходимо определить структуру основ-

ных блоков модели функционирования системы «танковая - пушка - снаряд - цель». При формировании структуры определены главные этапы, представленные на рисунке 1.

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

Уравнение поступательного прямолинейного движения БПС относительно ствола:

т 5

сн п пг соп

(1)

здесь - сила давления ПГ на БПС, Рсоп =жКен рсн

^ - сила сопротивления поступательному движению БП^ ^ = 2жЯвнкобпМ^х ; тс - масса БПС; Бп - поступательное движение БПС;Л - внутренний радиус канала ствола; р - давление пороховых газов на БПС; кобп - ширина обтюрирующего пояска;

^ - коэффициент трения между стволом и обтюрирующим пояском; и к- контактные напряжения между стволом и обтюрирующим пояском.

i Надоели баннеры? Вы всегда можете отключить рекламу.