Анализ и баллистическое проектирование системы с присоединенной камерой подгона Текст научной статьи по специальности «Физика»

УДК 629.7.018.3:623.5

АНАЛИЗ И БАЛЛИСТИЧЕСКОЕ ПРОЕКТИРОВАНИЕ СИСТЕМЫ С ПРИСОЕДИНЕННОЙ КАМЕРОЙ ПОДГОНА

©2009 А. И. Сафронов1, В. В. Потапенко2

1Тольяттинский госуниверситет 2 Тольяттинский военный технический институт

На основе модели двухфазной двухскоростной среды проводится численное моделирование нетрадиционной схемы с присоединённой камерой подгона. Возможно повышение начальных скоростей метаемых элементов с ^=4,55 кГ/дм3 примерно на 7% при неизменных максимальных давлениях на дно закрытого торца цилиндрического канала по сравнению с классической схемой метания.

Баллистическое проектирование, баллистическая система, присоединённая камера подгона

Схема с присоединённой камерой подгона (рис.1) позволяет повышать начальные скорости метания элементов, что необходимо для практики.

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

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

/— поршень /— метаемый /___________ / элемент

I п

—> !—> ,

Рис. 1. Схема метания с присоединённой камерой подгона:

I - область основного заряда; II - область присоединенной камеры подгона

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

ния топлива. В начальный момент времени в движение приходит вся сборка, состоящая из присоединённой камеры подгона (поршнь + топливо во второй области) и метаемого элемента. Через некоторый промежуток времени ^ (время задержки зажигания) воспламеняется заряд в присоединённой камере подгона (будем называть его присоединённым), в результате давление во второй области повышается; и в момент, когда сила сопротивления становится больше силы, ускоряющей сборку, последняя разделяется: метаемый элемент отделяется и далее движется вперед самостоятельно, а поршень притормаживается. При достижении метаемым элементом конца ствола расчет заканчивается. Основным положительным свойством этого процесса является то, что за счет работы присоединённого заряда происходит перераспределение энергии: часть энергии идет на дополнительное ускорение, вследствие которого достигается выигрыш в скорости по сравнению с системами использующими классическую схему метания.

Математическое моделирование схемы с присоединённой камерой подгона основано на подходе механики гетерогенных сред [1,2], модифицированного для внутренней баллистики ствольных систем в [3], и проводится при следующих основных допущениях:

- движение камеры подгона (поршень + присоединённый заряд) и метаемого элемента начинается при достижении давления форсирования;

- начальный период для основного и присоединённого зарядов не учитывается;

- горение частиц пороха происходит по геометрическому закону;

- вязкость и теплопроводность существенны только в процессах взаимодействия фаз;

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

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

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

- присоединённый заряд воспламеняется мгновенно в момент времени, соответствующий времени задержки зажигания;

- после воспламенения присоединённого заряда перетоки между областями I и II отсутствуют.

Система уравнений, описывающая

процессы в области I, записывается в инер-циальной системе координат (0, х), в области II - в неинерциальной системе координат (0', х'), связанной с поршнем (рис. 1). Далее приводится система уравнений, которая при N=0 описывает процессы в области I, а при N=1 и замене переменных и параметров на переменные и параметры со штрихом «'» (и=и' и т.д.) описывает процессы в области присоединенного заряда (области II) после воспламенения всего присоединенного заряда.

4 (р + 4 (ршр) = М;

— 4

9 ..4 2 ч -Р

— (рвфн) +-(ряри + рр) = М®-тТр + р-----------

а а а

N0* аиП;

йЬ

9 ( 0 Е) + 4 /о ( Е + чч 4(1 - р)Ою

— (рОрЕ) + — (Ори( рЕ + р)) = - р------------------

а ах ах

(1)

(2)

Ттрю + М

ю

- NрБщ

йи

(3)

П .

а а

- ( р 2(1 - р)0) +—(р 2(1 - Ф)0ю) = -м ; (4)

а ах

а а

— ( р2(1 - р) 0ю) + а ( р2(1 - р) 0ю2) +

а ах

(1 - р)о 4р = тТР - Мю - Nр2 (1 - р)0—П; 9х йЬ

а а а р

— + ю— = -^~

— -х еь

р

= ЯТ

(5)

(6)

(7)

Е = е + -

Ф = 1 - ПЛ о(1 - ф(ъ))

у(ъ) = К^(1 + 1^) М = Бп802 р 2 а(7)а1р

0(7) = 1 + 211 ъ 1

(8)

жй

''ТР

=- сх р(и - ю) |и - ю1 °п—т- *(1 -у( 2 ))2/3;

Сх =

Яе =

24

— + 0.48; Яе

0.1,

р I и - ю

0 < Яе < 3 • 105;

Яе > 3 • 105;

|Рл/°

02

Здесь Ь - время; х - координата; Р - давление; Р - плотность; Р2 - плотность вещества топлива; Т - температура; Т0 - температура продуктов горения топлива;

и - скорость газа; ю - скорость частиц; ф - пористость; Е,е - полная и внутренняя энергия единицы объёма газа; Б - площадь поперечного сечения канала; Ъ - относительная толщина сгоревшего свода; М - скорость массоприхода от горения топлива; тТР -сила взаимодействия между фазами; N -признак системы координат; ёиП/& - ускорение поршня; Q - тепловой эффект горения топлива; Я - универсальная газовая постоянная; а - коволюм; а1 - коэффициент в законе

2

2

2

скорости горения; ев - толщина горящего свода зерна топлива; п - концентрация; Л0 -начальный объём частицы топлива; ^(ъ) -относительный сгоревший объём частицы топлива; к1, 11 - коэффициенты формы частиц топлива; Б02 - начальная площадь частиц топлива; о(ъ) - относительная горящая поверхность частицы топлива; Сх - коэффициент сопротивления; ё0Р - диаметр шара эквивалентного по объёму частице топлива; р -число Пи; Яе - число Рейнольдса; т - вязкость газа.

Начальные условия в области I: Т(х,0)=То; Р(х,0)=Рф; и(х,0)=ш(х,0)=0; ф(х,0)=фн;

ёит

Ун =

1/А- 1/р2

7 г;

+ а------

2н -

2ун

+ ан)

ф

2

(9)

К1

н

ёит

т

СБ

dt

-

к

dt

- 0

t < tt

(11)

dt

- СР1 - Р2^

(12)

dt dt Здесь Р1 - давление на поршень слева; Р2 - давление на поршень справа; Р3 - давление на дно метаемого элемента; тСБ - масса сборки; тп - масса поршня и осевших на него частиц; - масса метаемого элемента;

tR - момент разделения сборки.

Момент разделения сборки tR определяется из условия разделения:

Рз > —(Р1- рг).

4 п

(13)

Здесь А - плотность заряжания; { - сила топлива; РФ - давление форсирования.

Индексы:

н - означает начальные значения переменных и параметров; о -относится к газу.

Граничные условия:

и(0Д) = ш(0Д) = 0; и(хпД) = ип;

(10)

и '(0,1) = о '(0,1) = 0; и' (х' бД) = и' б.

Здесь хп - координата положения левой границы поршня; хБ - координата положения дна метаемого элемента; ип - скорость поршня; иБ - скорость метаемого элемента.

Индексы:

п- означают, что параметры относятся к поршню; Б - к метаемому элементу.

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

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

Начальные условия во второй области для системы уравнений (1) - (7) с замыкающими соотношениями (8) определяются в момент достижения времени задержки зажигания присоединённого заряда 1;З для параметров этого заряда и давления форсирования РФ2 .

Система уравнений (1) - (12) решается численно методом С.К. Годунова [4], где потоки массы, импульса и энергии в уравнениях (1) - (3) определяются из решения задачи распада произвольного разрыва параметров газа на скачке площади сечения [5]. В уравнениях (4) - (5) потоки массы и импульса находятся из решения задачи распада произвольного разрыва для среды, не имеющей «собственного» давления. Уравнение (6) решается методом, аппроксимирующим конвективный член разностями против потока. Уравнения (11) - (12) решаются методом Эйлера, где давления Р1; Р2, Р3 определяются из решения задачи распада произвольного разрыва параметров газа на неподвижной непроницаемой поверхности [4].

Расчётные области I и II покрываются разностной сеткой с шагами И и И' соответственно.

Для проведения расчётов была выбрана система с диаметром канала 0,1 м с относительной длиной канала 7,0; Сч=4,55 кГ/дм3 для метаемого элемента,

Я

где Сч = -^33, й - диаметр канала баллистиче-й

ской системы; Сч=1,28 кГ/дм3 для поршня. Параметры составов основного и присоединённого зарядов были выбраны по данным

[3].

В таблице приведены максимальные безразмерные давления на дно закрытого торца цилиндрического канала Р!шах; в области II - Рпшах; на дно метаемого элемента

- РЭшах; конечные импульсы топлив в области I - 1к и в области II - I ' к; время задержки зажигания - 1:З для нескольких расчётных вариантов схемы с присоединённой камерой подгона, при различных значениях масс основного заряда ш1 и масс присоединённого заряда ш2. В строках 7,8 таблицы

При дальнейшем увеличении ш2 дульная скорость начинает падать, так как при обеспечении ограничения максимального давления в области II не выше 0,5937 появляется недогорание присоединённого заряда, а если уменьшать I \, обеспечивая полное сгорание присоединённого заряда, нарушается ограничение на максимальное давление в области II. Увеличение массы основного заряда позволяет увеличить дульную скорость до 1,718 при заданных ограничениях на максимальное давление (таблица, строка 6).

Результаты расчётов по двухфазной двухскоростной модели схемы выстрела с присоединённым зарядом хорошо согласуются с расчётными данными по модели газопороховой смеси при фиксированных пара-

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

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

При увеличении массы присоединённого заряда (за счёт увеличения плотности заряжания в области II при постоянном её начальном объёме) дульная скорость увеличивается (таблица, строки 1, 2, 3, 4) при близких значениях максимальных давлений в области II и на дно метаемого элемента.

метрах процесса метания.

Таким образом, реализация схемы с присоединённой камерой подгона на модельной установке с диаметром цилиндрического канала 0,1 м позволит повысить начальные скорости метаемых элементов с Сч=4,55 кГ/дм3 примерно на 7% при неизменных максимальных давлениях на дно закрытого торца цилиндрического канала по сравнению с классической схемой метания элементов.

Библиографический список

1. Рахматуллин, Х.А. Основы газодинамики взаимопроникающих движений сжимаемых сред/ Х.А. Рахматуллин //ПММ. -1956. - Т.20. - Вып. 1. - С.184-195.

Таблица2-Расчётные значения баллистических параметров схемы с присоединённой камерой подгона

№ п/п ю1 кг Ю2 кг 1з мс Р*мдх б/р Римлх б/р РЭшах б/р МПа-с I' к МПа-с Уд б/р

1 6,5 0,8 6,552 0,5934 0,5904 0,5813 1,02113 0,19426 1,657

2 6,5 0,9 6,553 0,5946 0,5910 0,5788 1,02546 0,23000 1,677

3 6,5 1,0 6,554 0,5936 0,5896 0,5851 1,03129 0,26888 1,690

4 6,5 1,1 6,578 0,5933 0,6009 0,5904 1,03645 0,30927 1,698

5 6,5 1,2 6,579 0,5929 0,5831 0,5754 1,04167 0,36019 1,691

6 7,0 1,1 6,504 0,5944 0,5924 0,5836 1,14897 0,32842 1,718

7 6,5 — — 0,5916 — 0,3076 0,95030 — 1,610

8 7,0 — — 0,5939 — 0,3004 0,99044 — 1,631

Примечание: б/р - означает, что параметры приведены в безразмерном виде.

2. Нигматулин, Р.И. Динамика многофазных сред. Ч.1. / Р.И Нигматулин - М.: Наука, 1987. - 464 с.

3. Русяк, И.Г. Внутрикамерные гетерогенные процессы в ствольных системах/ И.Г. Русяк, В.М. Ушаков. Екатеринбург: УрО РАН, 2001.-259 с.

4. Численное решение многомерных задач газовой динамики/ С.К. Годунов [и др.].,-М.:Наука.1976.- 400с.

5. Якушев, И.К. Распад произвольного разрыва в канале со скачком площади сечения/ И.К. Якушев // Известия СО АН СССР. -Сер. техн. наук.- 1967.-Вып. 2.- №8.-С.109-120.

References

1. Rahmatulin H.A. Bases of gazodimamiki

of interperetrative motions of the squeezed envi-ronments//Applied mathematics and mechanics-1956.-v.20-p. 184-195.

2. Nigmatulin R.I. Dynamics of many-phased environments. Part I. - M.: Science, 1987, p. 464

3. Rusyak I. G., Ushakov V. M. Inside-

chamber gcterogennye processes in the barrel systems. Ekaterenburg: Uralski department

RAS, 2001, p. 259

4. Godunov S.K., Zabrodin A.V., Ivanov M.L, Kraiko A.I., Prokopov G.I. Numeral decision of problem of gas dynamics.-M..:Science.-1976, p.400

5. Yakushev I.K. Disintegration of arbitrary gap in canal with bound of area section// News sD AS USSR, issue 2, № 8, series of technical sciencs.-1967, p. 109-120

ANALYSIS AND BALLISTICS PROJECT THE GUN TUBE SYSTEM WITH A JOINED

DRIVE ON CHAMBER

©2009 A. I. Safronov1, V. V. Potapenko2

1Togliatty State University,

2Togliatty military enginiring Institute

On the basis of the mathematical model of internal gun's ballistics with two phases and two velocities of phases the numerical simulation of non traditional scheme of throwing was carried out. The possibility of increasing of initial velocities of missiles with Cq=4,55 kilogram - force per cubic decimeter at about 7% comparison with classical scheme of throwing on conditions that maximum pressures are invariable.

Ballistics Project, Gun Tube System, Joined Drive on Chamber

Информация об авторах

Сафронов Александр Иванович, доктор физико-математических наук, старший научный сотрудник Тольяттинского государственного университета. Тел. 637401. E-mail: safr.a@mail.ru. Область научных интересов: механика гетерогенных систем, численное моделирование процессов в импульсных системах.

Потапенко Вячеслав Викторович, заместитель начальника факультета, начальник учебной части Тольяттинского военного технического института. Тел. +79608431311, Email: pww@mail.ru. Область научных интересов: механика гетерогенных систем, численное моделирование процессов в импульсных системах.

Safronov Alexander Ivanovich, Doctor, Head of Math. Physics Departament of Togliatty State University. Phone: 637401. E-mail: safr.a@mail.ru. Area of research: mechanics heterogeneous Systems, numerical modeling Processes in Impulses Systems.

Potapenko Vyacheslav Viktorovich, Deputy of battalion commander chief of study department of Togliatty military engineering institute. Phone: +79608431311. E-mail: pww@mail.ru. Area of research: mechanics heterogeneous Systems, numerical modeling Processes in Impulses Systems.