Preview

Вестник Концерна ВКО «Алмаз – Антей»

Расширенный поиск

Моделирование формирования порохового заряда, состоящего из зерненых пороховых элементов

Полный текст:

Аннотация

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

Для цитирования:


Хмельников Е.А., Заводова Т.Е., Смагин К.В., Дубинина С.Ф. Моделирование формирования порохового заряда, состоящего из зерненых пороховых элементов. Вестник Концерна ВКО «Алмаз – Антей». 2018;(3):69-81.

For citation:


Khmelnikov E.A., Zavodova T.E., Smagin K.V., Dubinina S.F. Simulation of the formation of powder charge consisting of granulated powder elements. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2018;(3):69-81. (In Russ.)

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

Для решения этой задачи в настоящее время широко используется газодинамический подход. Одна из наиболее развитых моделей горения пороха, движения пороховых газов и пороховых элементов на основе этого подхо­да описана в трудах [1, 2]. В существующей модели распределение пороховых элементов по длине заряда в начальный момент времени считается равномерным, что не вполне соответ­ствует реальным процессам, наблюдаемым при выстреле. Поэтому была разработана модель формирования порохового заряда, состояще­го из зерненых пороховых элементов, которая позволит учесть неравномерное распределение пороховых элементов по длине заряда.

Система уравнений, используемых при газодинамическом подходе, состоит из урав­нения неразрывности, уравнения движения, уравнения энергии, уравнения горения по­рохового заряда, а также дополнительных уравнений для вычисления сил межфазного взаимодействия и других параметров. Таким образом, система уравнений имеет следую­щий вид [1-4]:

где ρ j - плотность газа для части заряда с ин­дексом j;

m - пористость смеси;

S - переменная площадь поперечного се­чения каморы и канала ствола;

t - время;

v, w - скорость движения газа и твердой фазы в канале ствола соответственно; х - координата;

ρ = ρ12 - суммарная плотность смеси продуктов горения;

τ w1, τ w2 - сила взаимодействия, приходяща­яся на единицу объема, вызванная разностью скоростей между фазами для первой и второй частей заряда;

G1, G2 - газоприход с поверхности пороха в единице объема за единицу времени;

Пс - периметр ствола;

qc - тепловой поток, направленный на по­верхность канала ствола;

τc - сила трения газа о поверхность ствола, приходящаяся на единицу площади;

ε - внутренняя энергия единицы массы пороховых газов;

Q1, Q2 - теплотворная способность (потен­циал) пороха;

р - давление;

α - коволюм;

θ = k — 1, k = cp / cv - показатель адиабаты; cp, cv - теплоемкость газа при постоянном давлении и объеме;

aj - счетная концентрация зерненых поро­ховых элементов;

δ - плотность материала пороха;

ψj - относительная доля сгоревшего по­роха для первой и второй частей заряда;

Λ0, S0 - начальный объем и поверхность порохового зерна;

σ(ψ) - отношение текущей поверхности горения к первоначальной;

Uk - линейная скорость горения пороха;

æр - коэффициент формы порохового эле­мента;

e1 - половина начальной толщины горяще­го свода.

Индекс j в уравнениях обозначает номер части заряда для комбинированных зарядов.

Начальные и граничные условия имеют следующий вид [1, 2].

  1. Начальные условия (t = 0):

где р0 - давление форсирования;

ω j - масса j-го полузаряда;

Wj - объем каморы, занятый j-м полуза- рядом;

m0 - начальная пористость заряда;

Wкм - объем каморы орудия;

U1 j - скорость горения пороха при атмо­сферном давлении (единичная скорость горе­ния);

v j - показатель степени в уравнении для скорости горения;

fj - сила пороха.

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

При х = 0 (у дна каморы)

v = w = 0;

при х = хсн (у дна снаряда)

v = vсн

где хсн - координата дна снаряда в данный момент времени;

vсн - скорость снаряда в данный момент времени.

При расчете внутрибаллистических па­раметров в уравнениях газодинамического подхода используется такая величина, как по­ристость - m. Пористость - это объем пустот в единице объема, занятого пороховыми эле­ментами. Таким образом, пороховые элементы будут занимать объем 1 - m. По мере выгорания пороха пористость будет увеличиваться.

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

где Δ - плотность заряжания;

δ - плотность пороха.

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

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

  • при засыпке зерненого пороха в гильзу каждый из пороховых элементов принимает случайное положение и ориентацию в про­странстве;
  • после того как пороховой элемент по­пал на нижележащий слой пороха, он может продолжить движение вдоль боковых поверх­ностей уже сформированного объема пороха (т. е. может происходить «осыпание» заряда).

При моделировании формирования порохового заряда заснарядный объем, за­полняемый порохом, представляется в виде пустотелого цилиндра того же объема, распо­ложенного вертикально. В случае двумерного моделирования рассматривается только одно вертикальное сечение цилиндра, проходящее через его ось симметрии. При этом использу­ется система координат OX0Y0 (рис. 1).

 

Рис. 1. Система координат при моделировании фор­мирования порохового заряда:

OX0, OY0 - оси системы координат; D - диаметр ци­линдра; L - высота цилиндра; I - длина порохового элемента; d - диаметр порохового элемента

 

Цилиндр заполняется пороховыми эле­ментами. Каждый пороховой элемент возникает на верхнем срезе цилиндра в произвольной точ­ке оси OX0 - точка возникновения элемента со­ответствует расположению его центра тяжести в начальный момент времени. Затем он падает вниз до тех пор, пока не достигает дна цилин­дра или других пороховых элементов, располо­женных под ним. После этого на верхнем срезе цилиндра возникает новый пороховой элемент, который имеет произвольную ориентацию в пространстве, т. е. его ось симметрии может составлять любой угол α с осью OX0.

После того как падающий пороховой эле­мент достигает слоя нижележащих элементов, он останавливается или продолжает движение в зависимости от положения центра тяжести элемента относительно точки опоры. При рас­чете движения порохового элемента кроме ос­новной системы координат OX0Y0 для каждого падающего элемента строится еще одна систе­ма координат Oxyz (рис. 2, ось z на рисунке не показана), связанная с пороховым элементом. За начало координат принята точка O - центр тяжести элемента.

 

Рис. 2. Силы, действующие на падающий пороховой элемент:

Ox, Oy - оси системы координат; О - центр тяжести элемента; G - сила тяжести; N - сила реакции опоры; Fтp - сила трения; h - расстояние между осью Oy и ли­нией действия силы реакции опоры; r - расстояние меж­ду осью Ox и линией действия силы трения; α - угол поворота порохового элемента

 

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

qax = G sin α- Fтp;

где q - масса порохового элемента;

ax, ay - проекции ускорения элемента на оси х и у;

α - угол поворота порохового элемента относительно оси OX0;

Iz - момент инерции порохового элемента относительно оси z;

ε - угловое ускорение;

g - ускорение свободного падения;

f - коэффициент трения.

Если проекции ускорений представить в виде вторых производных координаты по вре­мени, можно получить следующую систему дифференциальных уравнений:

где х и у - координаты точки опоры в системе координат Oxyz; t - время.

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

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

  1. Переход к свободному падению эле­мента (рис. 3, а).

 

Рис. 3. Условия прекращения расчета: а - свободное падение элемента; б - остановка элемента; в - поступательное движение элемента; г - вращатель­ное движение элемента; д - изменение положения точки опоры; Ox, Oy - оси системы координат; O - центр тяжести элемента; G - сила тяжести; N1, N2 - силы реакции опор; Fтр1, Fтр2 - силы трения; h - расстояние между осью Oy и линией действия силы реакции опоры N1; ∆r - расстояние между осью Ox и линией дей­ствия силы реакции опоры N2; a - длина падающего элемента вдоль оси Ox; h2 - расстояние между осью Oy и линией действия силы реакции опоры N2; r - расстояние между осью Ox и линиями действия сил трения Fтр1 и Fтp2;  ω, ω′ - угловые скорости поворота элемента; C1 - первая точка опоры; C2 - вторая точка опоры

 

В процессе движения порохового эле­мента он может потерять контакт с другими пороховыми элементами, это возможно при выполнении условия

где a = l или a = d в зависимости от располо­жения в пространстве падающего элемента.

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

  1. Прекращение движения элемента (рис. 3, б).

В этом случае система уравнений изменя­ется с учетом появления второй точки контакта:

Если точки соприкосновения расположе­ны на двух взаимно перпендикулярных сторо­нах порохового элемента, то поступательное движение прекратится. При этом поворот от­носительно оси z может продолжаться в зави­симости от соотношения между величинами Δr, h и α. Можно доказать, что прекращение вращательного движения элемента произойдет в случае выполнения условия:

Если условие выполняется, то расчет дви­жения данного элемента прекращается, генери­руется новый пороховой элемент. Если условие не выполняется, то продолжается расчет пово­рота данного элемента относительно оси z.

  1. Прекращение вращательного движения при продолжении поступательного (рис. 3, в).

Если точки соприкосновения расположе­ны на одной стороне порохового элемента и линия действия силы тяжести проходит между точками опоры, то поступательное движение элемента вдоль оси х будет продолжаться, а вра­щение элемента прекратится. Для доказатель­ства этого утверждения достаточно составить два уравнения моментов относительно первой и второй точек опоры; решая их совместно, можно получить значение ε = 0. Для этого слу­чая справедлива следующая система уравнений:

В случае равенства нулю углового уско­рения угловая скорость поворота вокруг оси z остается постоянной величиной: ω = const. Так как к моменту появления второй точки опоры угловая скорость больше нуля, то по­ворот элемента будет продолжаться (рис. 3, г) и произойдет потеря контакта с первой точкой опоры. На пороховой элемент действует сила тяжести G, вызывающая замедление вращения элемента, поэтому угловая скорость ω посте­пенно уменьшится до нуля, а затем начнется поворот элемента в противоположном направ­лении с угловой скоростью ω', и он вернется в исходное положение. В некоторых случаях, при больших значениях угла α, угловая ско­рость вращения может остаться положитель­ной в момент достижения значения α= 90°, в этом случае происходит переворот элемента.

Для проверки переворота элемента ис­пользуются следующие уравнения:

где ω0 - значение угловой скорости ω в мо­мент отрыва от второй точки опоры;

εс - среднее значение углового ускорения при повороте элемента на угол απ0;

tп - промежуток времени, за который про­изойдет поворот на угол απ - α0;

αп = 90° - значение угла α в момент пере­ворота элемента;

α0 - значение угла α в момент отрыва по­рохового элемента от второй точки опоры.

Решая полученные уравнения, можно найти значение угловой скорости ω в момент переворота элемента т. е. когда α = απ = 90°. При этом если угловая скорость ω окажется больше нуля ω>0, то произойдет переворот порохового элемента; если ω ≤ 0, то перево­рот невозможен, т. е. пороховой элемент вер­нется в исходное положение с двумя точками опоры.

  1. Изменение положения точки опоры (см. рис. 3, д).

Если точки соприкосновения располо­жены на одной стороне порохового элемента и линия действия силы тяжести не проходит между точками опоры, то поступательное и вращательное движения элемента продолжат­ся, но зона контакта с нижним слоем элементов переместится из точки C1 в точку C2. После по­явлении второй точки контакта под действием силы тяжести угол α будет продолжать увели­чиваться, в результате элементы в зоне точки C1 перестанут соприкасаться. В этом случае единственной точкой контакта между порохо­выми элементами будет точка C2, а для расчета движения элемента можно использовать систе­му уравнений (1), приняв в качестве начальных координат точки контакта координаты C2.

Задача решалась в двумерной постанов­ке с помощью численного моделирования на ЭВМ. При этом рассматриваемое сечение раз­бивалось на отдельные элементы с размерами Δ х и Δу. Таким образом, был сформирован двумерный массив элементов (рис. 4).

 

 

В начальный момент времени всем эле­ментам массива присваивалось значение 0 - оно соответствует пустому, незаполненному элементу объема. Затем в некотором элементе массива появлялся пороховой элемент, его рас­положение (угол α ) определялось случайным образом. Для того чтобы получить координату центра тяжести элемента и его расположение, использовалась стандартная функция для ге­нерации случайных чисел.

После того как координата центра тяже­сти и ориентация порохового элемента опреде­лены, начинается расчет процесса его падения. При этом всем элементам массива, в которых в данный момент он располагается, присваива­ется значение, равное 2. Падение продолжает­ся до тех пор, пока под пороховым элементом остается свободное пространство.

В модели этот расчет реализован в виде проверки значений массива, располагающихся под участком, занятым пороховым элементом. Если хотя бы один из элементов массива под по­роховым элементом имеет значение, равное 1, то перемещение порохового элемента прекра­щается и процесс расчета падения элемента заканчивается. Значение 1 соответствует эле­ментам массива, в которых располагается не­подвижный пороховой элемент.

При моделировании движения порохово­го элемента уравнения основной системы ре­шаются численно с помощью метода Эйлера. В данном случае используется регулярная вре­менная сетка. Она образуется линиями tn = n% (n = 0, 1, 2, ...), где τ - шаг по времени. Значения рассчитываемых параметров (vx, ω, х, α) определяются в целых точках по времени т . Таким образом, можно вычислить линейную скорость vx и угловую скорость по­ворота ω в любой момент времени, а также значения координаты х и угла α:

При расчете по уравнениям для каждого момента времени определяются координаты центра тяжести порохового элемента и других точек в основной системе координат ОХ0 Y0:

где Хi, Yi - координаты точки с номером i в системе координат ОХ0 Y0.

Схема расчета представлена на рис. 4.

После определения координат точек 3, 4, 5 и 6 строится контур движущегося порохового элемента в момент времени tn .

Если более 50 % площади какого-либо элемента массива находится внутри постро­енного контура 4-6-3-5, то считается, что

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

После определения положения падаю­щего порохового элемента проводится про­верка выполнения условий прекращения рас­чета. Если ни одно из условий не выполняется, расчет продолжается и определяется положе­ние порохового элемента в следующий мо­мент времени tn+1. Если какое-либо условие выполняется, порядок расчета изменяется в соответствии с правилами, описанными выше. После полной остановки порохового элемента элементам массива, в которых в данный мо­мент он находится, присваивается значение 1 и процесс расчета перемещений порохового элемента заканчивается, генерируется новый пороховой элемент.

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

При расчете пористости в каждом из гори­зонтальных сечений определяется число пустых элементов массива (т. е. имеющих значение 0) и число заполненных элементов (имеющих значе­ние 1). Затем вычисляется значение пористости для данного горизонтального сечения:

где mK - значение пористости для горизон­тального сечения, имеющего координату К;

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

На основе составленной модели был про­веден ряд расчетов по формированию поро­хового заряда для пушки АК-230. На рис. 5 показано заполнение цилиндра пороховыми элементами. Рисунок соответствует одному из расчетов и изображает часть заснарядного про­странства (длиной 48 мм) в начале (рис. 5, а) и в конце (рис. 5, б) его заполнения.

 

Рис. 5. Начало (а) и окончание (б) заполнения цилиндра:

D = 40 мм, L = 48 мм, d = 2,5 мм, l = 2,9 мм, ∆х = ∆у = 0,1 мм

 

В табл. 1 приведены результаты 20 расче­тов, выполненных при начальных данных, со­ответствующих конструктивным параметрам 30-милиметровой пушки АК-230.

 

Таблица 1

Результаты расчетов

№ расчета

mснар

mдк

расчета

mснар

mдк

1

0,4649

0,9540

0,3324

11

0,4623

0,5618

0,5466

2

0,5419

0,7042

0,3485

12

0,5330

0,8011

0,3738

3

0,4934

0,6947

0,3406

13

0,4955

0,8136

0,4399

4

0,5194

0,8019

0,5116

14

0,5000

0,7120

0,3904

5

0,5422

0,8878

0,3732

15

0,5042

0,8687

0,4839

6

0,5060

0,8088

0,4876

16

0,4879

0,8363

0,3283

7

0,4811

0,7179

0,3192

17

0,5060

0,7029

0,3457

8

0,4836

0,5658

0,4794

18

0,4889

0,7079

0,4597

9

0,5232

0,8496

0,3256

19

0,4815

0,7501

0,3477

10

0,4977

0,6772

0,3933

20

0,5433

0,8872

0,3501

Примечание:  — средняя по объему пористость, mснар — пористость на верхнем срезе рассматриваемого цилиндра (у дна снаряда), mдк — пористость у дна каморы.

Анализ полученных данных показывает, что значения пористости для различных горизон­тальных сечений цилиндра могут значительно отличаться друг от друга. Средняя пористость заряда также изменяется в некоторых преде­лах - она зависит от уровня заполнения гильзы пороховыми элементами (вдоль оси OY0).

На рис. 6 представлены графики распре­деления пористости по длине каморы. Первый график соответствует расчету 3 (см. табл. 1) с резким изменением пористости в районе дна каморы, второй - расчету 8 с наиболее равно­мерным из всех расчетов распределением по­ристости по длине заряда.

 

Рис. 6. Распределение пористости по длине каморы: а - расчет 3; б - расчет 8; Y0 - координата (значение «0» соответствует дну каморы), m(Y0) - среднее значение пористости в данном сечении

 

При моделировании формирования за­ряда наблюдалась значительная разница в по­ристости для соседних сечений. Например, в расчете № 3 для сечений с координатами от 26 мм до 32 мм значение пористости изменя­ется от 0,2217 до 0,7322. Такое неравномерное распределение пористости может приводить к значительным перепадам давления пороховых газов для соседних сечений каморы при расче­те внутрибаллистических параметров.

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

Графики изменения давлений пороховых газов в районе дна боеприпаса - рс (х), в райо­не дна каморы - рк (х) , а также графики измене­ния скорости боеприпаса V показаны на рис. 7. Кроме того, приведены кривые изменения давления пороховых газов в случае постоян­ной пористости по длине заряда. Все кривые на графиках представляют собой зависимости давлений пороховых газов от координаты дна боеприпаса х. Нулевое значение по оси х со­ответствует дну каморы орудия.

 

 

На основании полученных данных уста­новлено, что характер распределения пори­стости по длине заряда значительно влияет на процесс горения порохового заряда и, соот­ветственно, на давление пороховых газов, на­блюдаемых в процессе выстрела. Для зарядов с одинаковой средней пористостью, но разным характером распределения пористости, макси­мальные давления пороховых газов могут от­личаться весьма значительно. Так, например, в расчетах 7 и 8 при близких значениях средней пористости разница между максимальными давлениями пороховых газов составляла око­ло 70 МПа.

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

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

 

Таблица 2

Результаты расчета горения порохового заряда

№ расчета

Pmax м,
МПа

Pсм,
МПа

Pкм,
МПа

V, м/с

№ расчета

Pmax м,
МПа

Pсм,
МПа

Pкм,
МПа

V, м/с

1

351,4

331,4

344,9

937,3

11

358,1

344,4

353,8

941,9

2

354,9

349,8

347,5

941,9

12

381,9

367,8

374,1

958,8

3

366,6

346,2

352,9

949,9

13

398,9

381,3

397,8

962,9

4

343,6

340,9

340,8

936,2

14

360,1

347,9

355,4

948,1

5

343,7

336,2

325,5

932,1

15

386,0

369,2

386,0

960,5

6

375,2

366,1

374,4

952,9

16

393,3

357,7

384,2

957,2

7

406,8

364,8

395,4

966,5

17

347,2

342,7

347,2

936,8

8

336,9

331,3

316,1

926,1

18

376,1

357,4

354,1

949,2

9

323,9

323,9

318,2

916,2

19

345,6

334,0

332,2

929,5

10

357,7

343,1

353,4

941,7

20

333,2

328,2

328,6

924,9

Примечание: Pmax м — максимальное давление пороховых газов за все время расчета для всего заснарядного пространства; Pсм — мак­симальное давление пороховых газов за все время расчета у дна боеприпаса; Ркм — максимальное давление пороховых газов за все время расчета у дна каморы; V — скорость боеприпаса при вылете.

Для расчета при условии постоянной по­ристости значения максимальных давлений по­роховых газов и скорости боеприпаса при вы­лете составили: Pmax м = 377 8, МПа, Pсм = 369 9, МПа, Pкм = 353 3, МПа, V = 946 5602, м/с.

Экспериментальные значения макси­мальных давлений пороховых газов при вы­стреле для этого же орудия находятся в пре­делах от 307,034 до 362,012 МПа (данные получены при испытаниях артиллерийской системы на полигоне ФКП «НТИИМ») - один из экспериментальных графиков показан на рис. 7.

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

Список литературы

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

2. Русяк И.Г., Липанов А.М., Ушаков В.М. Физические основы и газовая динамика горения порохов в артиллерийских системах. М.; Ижевск: Институт компьютерных исследований, 2016. 456 с.

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

4. Русяк И.Г. Математические модели проектирования боеприпасов на основе уравнений механики гетерогенных реагирующих сред // Вопросы оборонной техники. 2011. № 14. С. 3-11.


Об авторах

Е. А. Хмельников
Нижнетагильский технологический институт (филиал) федерального государственного автономного образовательного учреждения высшего образования «Уральский федеральный университет имени первого Президента России Б.Н. Ельцина»
Россия


Т. Е. Заводова
Нижнетагильский технологический институт (филиал) федерального государственного автономного образовательного учреждения высшего образования «Уральский федеральный университет имени первого Президента России Б.Н. Ельцина»
Россия


К. В. Смагин
Нижнетагильский технологический институт (филиал) федерального государственного автономного образовательного учреждения высшего образования «Уральский федеральный университет имени первого Президента России Б.Н. Ельцина»
Россия


С. Ф. Дубинина
Нижнетагильский технологический институт (филиал) федерального государственного автономного образовательного учреждения высшего образования «Уральский федеральный университет имени первого Президента России Б.Н. Ельцина»
Россия


Для цитирования:


Хмельников Е.А., Заводова Т.Е., Смагин К.В., Дубинина С.Ф. Моделирование формирования порохового заряда, состоящего из зерненых пороховых элементов. Вестник Концерна ВКО «Алмаз – Антей». 2018;(3):69-81.

For citation:


Khmelnikov E.A., Zavodova T.E., Smagin K.V., Dubinina S.F. Simulation of the formation of powder charge consisting of granulated powder elements. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2018;(3):69-81. (In Russ.)

Просмотров: 38


Creative Commons License
Контент доступен под лицензией Creative Commons Attribution 4.0 License.


ISSN 2542-0542 (Print)