УДК 623.4.082.6; 533.95
ЭЛЕКТРОДИНАМИЧЕСКАЯ ЗАЩИТА ОТ КУМУЛЯТИВНЫХ СРЕДСТВ ПОРАЖЕНИЯ: ФИЗИЧЕСКИЕ АСПЕКТЫ ФУНКЦИОНИРОВАНИЯ
С.В. Федоров
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]
Рассмотрен нетрадиционный способ защиты различных объектов от кумулятивных средств поражения посредством реализации мощного электрического разряда через кумулятивную струю для уменьшения глубины ее проникания в преграду. В качестве возможных физических механизмов снижения пробивной способности кумулятивных струй при токовом воздействии проанализированы развитие магнитогидродинамической неустойчивости перетяжечного типа и диспергирование материала струи с его радиальным рассеиванием. На основании обобщения результатов численного моделирования данных эффектов разработана методика расчета пробивного действия кумулятивных зарядов в условиях пропускания по кумулятивной струе мощного импульса электрического тока. Проведено сопоставление результатов расчетов с имеющимися экспериментальными данными.
Ключевые слова: взрыв, кумулятивный заряд, кумулятивная струя, пробивное действие, электрический ток, электродинамическая защита.
ELECTRODYNAMIC PROTECTION AGAINST SHAPED CHARGE WEAPONS. PHYSICAL ASPECTS OF FUNCTIONING
S.V. Fedorov
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]
A nonconventional method is considered for protecting different objects against the shaped-charge weapons by means of passing a powerful electric discharge through a shaped-charge jet to decrease a depth of its penetration into the target. Development of the necking-type magnetohydrodynamic instability and dispersion of jet material with its radial scattering are analyzed as possible physical mechanisms of a reduction in the penetrative ability of shaped-charge jets under exposure to the electric current. Based on generalization of results of numerical simulation of these effects, a computational procedure is developed for determining the penetration action of shaped charges under condition ofpassing a powerful pulse of electric current along the shaped-charge jet. The calculation results are compared with experimental data.
Keywords: explosion, shaped charge, shaped-charge jet, penetration action, electric current, electrodynamic protection.
Эффективным средством поражения сильно защищенных бронированных объектов являются кумулятивные боеприпасы [1], формирующие при взрыве высокоскоростные металлические кумулятивные струи (КС), имеющие высокую пробивную способность. В настоящее время в качестве основного средства обеспечения противокумулятивной стойкости объектов военной техники выступает динамическая защита [2], принцип действия которой основан на взрывном метании
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 79
металлических пластин, сбивающих КС. Метание пластин элемента динамической защиты осуществляется заключенным между ними слоем взрывчатого вещества, инициируемого при пробитии КС. Несмотря на свою достаточно высокую эффективность, широкое применение динамическая защита получила главным образом на тяжело бронированных объектах военной техники (танках). Ее использование на легко бронированной технике ограничивается проблемами, связанными с собственным достаточно мощным взрывным воздействием элемента динамической защиты на защищаемый объект и на его экипаж.
В последнее время в связи с активизацией деятельности различного рода террористических группировок, широко использующих кумулятивные средства ближнего боя (реактивные противотанковые гранаты) вследствие их относительной доступности и простоты применения, приобрела актуальность противокумулятивная защита и ряда гражданских объектов. К их числу относятся, прежде всего, хранилища и транспортные контейнеры для перевозки ядовитых, горючих, взрывоопасных и радиоактивных материалов. Разрушение или повреждение подобных объектов в результате террористической атаки способно привести к серьезным последствиям — многочисленным человеческим жертвам, химическому и радиоактивному заражению территории. Использование в данном случае для противодействия кумулятивным средствам поражения динамической защиты практически исключается из-за опасности поражения при ее функционировании находящихся вблизи защищаемого объекта других гражданских объектов и людей разлетающимися с высокой скоростью фрагментами пластин. Кроме того, размещение на защищаемом объекте большого количества взрывчатого вещества (в элементах динамической защиты) требует дополнительных специальных мер по обеспечению его пожаровзрывобезопасности.
В качестве возможной альтернативы динамической защите при противодействии кумулятивным средствам поражения может рассматриваться так называемая электродинамическая защита (ЭДЗ) [3], принцип действия которой основан на разрушении КС кратковременным мощным импульсом электрического тока. В простейшем варианте устройство ЭДЗ представляет собой две металлических пластины (электрода), разделенные слоем диэлектрика и соединенные с источником электрической энергии (рис. 1). Такой элемент ЭДЗ устанавливается перед защищаемой преградой и после пробития струей слоя диэлектрика и замыкания электрической цепи осуществляет токовое воздействие на участки КС при их движении в межэлектродном промежутке. В отличие от динамической защиты функционирование ЭДЗ не представляет опасности для объектов, расположенных вблизи защищаемого, что существенно расширяет область ее применения.
80 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
Рис. 1. Схема противокумулятивной защиты с реализацией токового воздействия на кумулятивную струю
Первые исследования по ЭДЗ были выполнены более 30 лет назад в Институте гидродинамики им. М.А. Лаврентьева СО РАН [4, 5], где и было предложено использовать электродинамическое воздействие (ЭДВ) на КС для снижения ее пробивной способности. За прошедшие годы в рамках исследований ЭДЗ накоплен довольно обширный материал в основном экспериментального характера. В нашей стране эксперименты по разрушению КС импульсом тока проводились в Институте гидродинамики им. М.А. Лаврентьева [4-7], РФЯЦ-ВНИИЭФ [7-9] и ОАО “НИИ Стали” [7, 10]. В США подобные работы велись фирмой Maxwell Laboratories (Калифорния) и институтом Southwest Research Institute (Техас), в Англии — центром Defence Research Agency [11, 12]. В экспериментах исследовалась в основном двухэлектродная схема ЭДЗ (см. рис. 1), в качестве источника электрической энергии в большинстве случаев использовалась конденсаторная батарея, в [7] была также исследована возможность применения для токового воздействия на КС магнитокумулятивного генератора. Экспериментальные результаты получены для кумулятивных зарядов (КЗ) диаметром от 30 до 100 мм и представляют собой в основном данные по снижению пробивного действия зарядов в сочетании с фиксируемыми в экспериментах кривыми изменения разрядного тока в межэлектродном промежутке, а также рентгенограммы КС, испытавших мощное токовое воздействие [4, 6].
Для решения вопросов, связанных с практическим применением ЭДЗ, необходимо иметь соответствующие расчетные методики, позволяющие прогнозировать пробивное действие КЗ при пропускании через КС мощного импульса электрического тока. Очевидно, создание подобных методик невозможно без выявления физических причин уменьшения проникающей способности КС при токовом воздействии и разработки соответствующих физико-математических моделей поведения струи при ЭДВ.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 81
Рис. 2. Рентгенограммы разрушения кумулятивных струй мощным импульсом электрического тока:
а — без токового воздействия; б, в — при токовом воздействии различной интенсивности
Анализ данных рентгенографических исследований КС после мощного токового воздействия (рис. 2) [13] позволяет в качестве одной из основных причин снижения их проникающей способности рассматривать развитие МГД-неустойчивости перетяжечного типа [14]. Вследствие МГД-неустойчивости ускоряется процесс распада удлиняющейся КС на отдельные движущиеся друг за другом элементы, что приводит к уменьшению эффективной длины струи и, соответственно, ее пробивного действия. Как свидетельствуют рентгенограммы, элементы КС, на которые она распадается после мощного ЭДВ, имеют существенно меньшее удлинение, чем в случае естественного распада струи вследствие развития пластической неустойчивости [15]. При этом радиальные размеры элементов распавшейся в результате токового воздействия КС могут увеличиваться.
Расчетно-теоретическое исследование МГД-неустойчивости металлических КС при пропускании по ним мощного импульса электрического тока проводилось на основе модели, в которой элементы струи рассматривались как участки удлиняющегося осесимметричного несжимаемого жесткопластического стержня переменного радиуса в предположении сохранения плоской формы его поперечных сечений в течение всего процесса деформирования. Использование гипотезы плоских сечений позволяет выделить на участке стержня произвольное поперечное сечение толщиной dz (рис. 3) и записать для него уравнение движения в осевом направлении в виде
2 dv z dFz
nR p~dt = aZ"
„ „ dR
+ 2nRpm —,
dz
(1)
где R — радиус выделенного сечения; p — плотность материала стержня; vz — осевая скорость движения поперечного сечения; Fz — полная осевая сила, действующая в сечении, pm — приложенное к поверхности стержня магнитное давление, учитывающее действие на него электромагнитных сил, обусловленных протеканием по стержню электрического тока.
82 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
Действующая в сечении осевая сила Fz определяется распределением по нему осевой компоненты tz тензора напряжений:
R
Fz = 2п az rdr.
(2)
о
При этом значения компонент девиатора тензора напряжений Sj,
необходимые для расчета механических напряжений (в том числе, и осевого напряжения tz), находились с использованием соотношений Сен-Венана-Леви-Мизеса [16] для несжимаемой жесткопластической среды:
2 ту .
Sij = eij,
(3)
где ту — предел текучести материала стержня; e»j — компоненты тензора скоростей деформации; ё» — интенсивность скоростей деформации.
Отметим, что интегрирование уравнения неразрывности для произвольного поперечного сечения стержня дает линейный закон распределения по его радиусу r радиальной скорости vr частиц стержня
Vr = —£z r/2,
(4)
где ez = dvz/dz — текущая осевая скорость деформации в данном сечении. С учетом данного обстоятельства, а также следующего из принятой гипотезы плоских сечений условия равномерности распределения осевой скорости vz по радиусу сечения dvz/dr = 0 получаем радиальную er, окружную ёв и сдвиговую ezr компоненты тензора скоростей деформации в виде
dvr ez . vr ez . 1 dvr 1 dez
r dr 2 ’ в r 2 ’ zr 2 dz 4 dz ’
что, в свою очередь, для определения интенсивности скоростей деформации
ё» =
ёв )2 + (ee
ez )2 + (er
дает соотношение
r2 (de
ё \ ez +12\ dz
2
ez)2 + 6e2r
В результате нормальные компоненты девиатора тензора напряжений sz, sr и se с использованием (3) представляются выражениями
Sz =
ту ez
,ё2 + r! (де±
‘ + 12 \ dz
Sr = Se = —
ту ez
'«+s (f2
2
(5)
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 83
Действующее в материале стержня среднее давление а (дающее в сумме со значением sz осевое напряжение az) определялось с привлечением уравнения движения частиц стержня в радиальном направлении, записываемого в пренебрежении влиянием касательных напряжений Tzr в виде
dvr dsr да
^ dt дг дг
Проинтегрировав данное уравнение по радиусу стержня от некоторого его текущего значения г до радиуса наружной поверхности R, получим радиальное распределение среднего напряжения
R
dv
а (г) = Sr (R) - Sr (г) + a(R) - р J -t dr, (6)
r
в котором среднее напряжение на боковой поверхности стержня a(R) можно определить из динамического граничного условия (записываемого, по-прежнему, в пренебрежении влиянием касательных напряжений Tzr) равенства полного напряжения на поверхности стержня магнитному давлению pm, взятому с противоположным знаком:
(sr(R) + a(R)) cos2 a + (sz(R) + a(R)) sin2 a = -pm,
где a — угол наклона касательной к боковой поверхности стержня в данном сечении (рис. 3).
Приняв во внимание приведенное граничное условие, а также то, что в соответствии с (4)
dvr
dt
dvr dvr г
dt + Vr дг 4
2
z
2 d^
2 dt
1
можно переписать (6) в виде
а (г) = -sr (г) (sz (R) - sr (R)) sin2 a-8
dez
2-z
dt
R2 - г2
pm.
Рис. 3. Расчетная схема деформирования участка кумулятивной струи при пропускании по нему электрического тока
84 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
Определяя далее с использованием (5) радиальное распределение осевых напряжений oz (r) = sz (r) + a(r) в стержне и интегрируя его для вычисления осевой силы Fz (2) по площади поперечного сечения, преобразуем уравнение движения стержня (1) к виду
R2 dvz dS р д
^2 dt dz 32 dz где
ё2 -2 § >R4
R2 dpm 2 dz
(7)
12aY £^
(d£z/dz)2
o\ 1/2
R (d£z/dzf\ _ . 12 £2 )
aYR2 sin2 a ( R2 (d£z/dz)2
2 1 +l2 £2
1/2
характеризует прочностную составляющую осевой силы (так как зависит от предела текучести материала стержня).
Дополнив уравнение (7) вытекающим из (4) законом изменения радиуса произвольного сечения стержня dR/dt = -£zR/2 и учитывая, что sin2 a = (dR/dz)2 / (1 + (dR/dz)2), приходим к описанию динамики деформирования участка КС с возмущениями боковой поверхности фактически в рамках одномерной (квазидвумерной) задачи, где все искомые функции vz, £z и R зависят от времени t и лишь одной осевой координаты z (см. рис. 3).
В представляемой модели принимался также во внимание эффект термического разупрочнения материала стержня (снижения его предела текучести с ростом температуры). При этом предполагалось, что стержень нагревается за счет выделения джоулевой теплоты и пластической диссипации энергии. Учитывающее данные тепловые факторы уравнение энергии, записываемое для произвольного сечения стержня, имеет вид
nR2pc
dT
dt
R R
2n
nj 2rdr + 2n
oY £,trdr,
0 0
где c — теплоемкость материала стержня; T — усредненная по сечению температура; n — удельное сопротивление материала; j — плотность тока в сечении.
Вычисляя интегралы в правой части уравнения энергии в предположении равномерности распределения тока по сечениям стержня j = J/nR2, где J — сила тока в стержне, получаем дифференциальное
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 85
соотношение, описывающее рост его температуры:
dT 1 I n J2 2oY \ё2\
dt pc | n2R4 + (0iz/Oz)2 R2
1 +
R2
12i|
Принималось, что значение предела текучести материала стержня с ростом его температуры снижается по линейному закону [17] вплоть до полного исчезновения прочности при достижении температуры плавления материала oY — оY0 (Ts — T) / (Ts — T0), где oY0 — предел текучести материала стержня при нормальной температуре T0; Ts — температура плавления. Учитывалась также зависимость удельного сопротивления материала стержня от температуры в виде n = П0 (1 + в (T — T0)) [18], где n0 — удельное сопротивление при нормальной температуре T0; в — температурный коэффициент сопротивления.
Задаваемое на поверхности стержня для учета силового действия электрического тока магнитное давление (см. рис. 3) зависит от силы тока и размера поперечного сечения стержня [18]:
P0J 2
Pm —
8n2R2 ’
(8)
где ^0 — 4п 10-7 (Гн/м) — магнитная постоянная. В соответствии с данным соотношением при наличии на стержне неоднородностей размеров его поперечных сечений магнитное давление будет возрастать в областях перетяжек (участков с меньшим радиусом). Данное обстоятельство должно приводить к прогрессирующему утончению стержня в областях появления перетяжек, в чем, собственно говоря, и состоит суть эффекта МГД-неустойчивости.
Процесс растяжения КС является неустойчивым и в отсутствие токового воздействия. Развитие на КС перетяжек с последующим разрывом струи на отдельные элементы связано в естественных условиях с проявлением пластической неустойчивости [15]. На основе выведенного уравнения движения удлиняющегося стержня (7) можно получить оценку параметров электродинамического воздействия, при которых МГД-неустойчивость способна существенно повлиять на процесс развития естественной пластической неустойчивости. В предположении малости возмущений на поверхности стержня и равномерности распределения осевой скорости деформирования £z по его длине уравнение (7) может быть переписано в виде
R2 dvz _ udR R2 dpm
Р 2 dt °Y dz 2 dz ’
где первое слагаемое в правой части определяет темп развития пластической неустойчивости, а второе — МГД-неустойчивости. Принимая
86 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
во внимание (8), из условия равенства абсолютных значений этих слагаемых находим, что влияние токовой неустойчивости должно заметно проявляться при силе тока, сопоставимой с ее критическим значением JKp = 2nRy/2aY/^о. Нетрудно установить, что при критическом токе обеспечивается равенство магнитного давления пределу текучести материала стержня.
Особенности проявления МГД-неустойчивости и ее влияния на процесс растяжения КС исследовались в рамках сформулированной модели на основе подхода, предложенного в работе [19] по отношению к естественным образом деформирующейся КС и заключающегося в определении эволюции изначально задаваемых малых поверхностных возмущений на струе. Задача решалась численно с использованием конечно-разностного метода, представленного в [20]. В качестве материала высокоградиентного стержня рассматривалась медь.
Путем численных расчетов было установлено, что ЭДВ может не только ускорить развитие естественной пластической неустойчивости КС, приводя к ее ускоренному распаду на отдельные элементы, но и при достаточной интенсивности способно привести к явлению “дискообразования”, при котором элементы струи становятся плоскими, сжимаясь в осевом направлении (рис. 4). Следствием дискообразования является существенное сокращение эффективной длины КС. Следует отметить, что в процессе развития МГД-неустойчивости происходит резкое возрастание плотности тока в областях перетяжек на
2 I
0,4
0,2
-0,2
R, мм t = 0
£„ МКС
10 мкс ' -
- -
' 20 мкс'
А ) \
' 30 мкс' '
0 5 10 0 5 10 0 5 10 15 0 5 10 z, мм
0,4
0,2
-0,2
- ^ /W\ Л. - 'I А /| \ A h
_ - ■J W V - / У \1 - / У У
Лл/ \j
FkF
0 5 10 0 5 10 0 5 10 15 0 5 10
Рис. 4. Расчетное изменение формы участка кумулятивной струи и распределения осевой скорости деформирования по его длине:
а — при удлинении в естественных условиях; б — при пропускании импульса тока силой 400 кА в течение 5 мкс
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. №3 87
струе. Следствием этого в реальной ситуации может быть электрический взрыв материала данных областей (за счет его перехода в газообразное состояние вследствие мощного джоулева тепловыделения). В рассматриваемой модели возможность проявления данного эффекта не учитывалась, так как предполагалось, что его влияние на пробивную способность испытавшей токовое воздействие КС несущественно. В областях перетяжек при достаточно мощном ЭДВ практически не остается материала, поэтому его фазовое состояние с точки зрения вклада в пробитие КС не играет большой роли.
Для того чтобы иметь возможность использовать результаты численного моделирования развития МГД-неустойчивости КС при прогнозировании их пробивного действия, были введены в рассмотрение количественные показатели состояния струи, подвергающейся токовому воздействию. В число этих показателей входят коэффициент формы kf, показывающий отклонение формы участка КС от цилиндрической, коэффициент скоростного инвертирования kv, характеризующий перераспределение осевой скорости по длине участка струи, и коэффициент термического разупрочнения kT, учитывающий снижение прочностных свойств материала струи вследствие его нагрева [21]. Введенные коэффициенты определяются соотношениями
kf = R; kv =
€2
R2
‘"m
kY —
ау„
(9)
€zu &Y 0
где Rm, €zm и aYm — средние по массе участка КС радиус, градиент осевой скорости и предел текучести соответственно; Ru, €zu — радиус и градиент осевой скорости при условии равномерного растяжения участка струи с сохранением цилиндрической формы; aY0 — начальное значение предела текучести.
Входящие в (9) средние значения Rm, €zm и aYm рассчитывались по формулам:
i i i
R3dz
€z R2dz
aY R2dz
Rm —
j €zm
R2dz
"j &Ym
R2dz
R2dz
0 0 0 где l — текущая длина участка КС. Значения Ru и €zu в зависимости от коэффициента удлинения участка струи n — l/l0 (l0 — начальная длина участка) определяются по соотношениям Ru — R0 Д/n и €zu — €z0/n, где R0 и €z0 — начальные значения радиуса и осевой скорости деформирования.
88 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
При естественном деформировании (рис. 4, а) введенные коэффициенты (9), уменьшаясь от начального единичного значения, к моменту распада КС на безградиентные элементы оказываются следующими: kf = 0,6... 0,7; kv = 0; kT ~ 1 (коэффициент kT становится чуть меньше единицы за счет разупрочнения, обусловленного нагревом вследствие пластической диссипации энергии). При достаточно интенсивном ЭДВ (рис. 4, б) коэффициент формы резко уменьшается, коэффициент скоростного инвертирования становится отрицательным (соответствует осевому сжатию материала в областях формирующихся дисков), а коэффициент термического разупрочнения приближается к нулевому значению.
Обобщение результатов численных расчетов показало, что для описания развития МГД-неустойчивости КС при токовом воздействии с приемлемой точностью могут быть использованы три безразмерных определяющих комплекса:
и =
P£ZQR0 Vy 0
s =
Vo J2
8n2 R^ay 0 ’
q =
По (1+ в(Tro - To)) J2 n2izoRtpc (Ta - Tro) ’
(10)
где Tro — начальная температура участка КС. Параметр и, равный отношению действующих в струе инерционных сил к пределу текучести, определяет процесс развития естественной пластической неустойчивости КС [15]. Параметр s характеризует силовое действие тока (отношение магнитного давления pm (8) к пределу текучести), а параметр q — термическое действие тока (отношение количества джоу-левой теплоты, выделяющейся в стержне радиусом Ro за характерное время деформирования 1/ezo, к количеству теплоты, необходимому для нагрева стержня до температуры плавления). При построении аппроксимационных зависимостей использовался также безразмерный параметр
= q = 8 aY ono (1 + в (Tro - To)) qi s po£zoR‘2pc (Ts - Tro) ’
удобный тем, что он зависит только от характеристик участка КС (не включает интенсивности токового воздействия).
В проведенных расчетах параметр и изменялся в диапазоне от 0,05 до 0,5, что соответствует окончанию инерционной и переходу к шеечной стадии деформирования КС [15] (в реальной ситуации токовое воздействие на КС может быть осуществлено как раз на данном этапе ее эволюции). Показатель динамического действия тока варьировался от 0 до 5.
Одной из важнейших характеристик процесса растяжения КС является коэффициент предельного удлинения ее элементов, определяемый отношением длины участка КС в момент его распада на отдельные элементы к начальной длине участка. При деформировании в естественных условиях значение коэффициента предельного удлине-
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 89
ния nlim зависит от параметра и. С учетом результатов [15] данная зависимость может быть принята в виде nlim = 5, 8 и1/3. Токовое воздействие на КС приводит к уменьшению предельного удлинения ее элементов. Достаточно близкое соответствие с результатами численных расчетов по значению коэффициента предельного удлинения nfim в условиях ЭДВ дает соотношение
n
e
\im
1 . n\im 1
+ 1 + 2, 5 s0,5,
(11)
включающее в себя из безразмерных комплексов, характеризующих интенсивность ЭДВ (10), лишь параметр s (влияние параметра термического действия тока q на величину neim не очень существенно).
В расчетах по определению nim предполагалось, что ток по участку КС протекает вплоть до его разделения на отдельные элементы (достижения радиусом перетяжки нулевого значения). На этот момент времени фиксировались также значения коэффициентов формы kfb и скоростного инвертирования kvb. Характер их изменения в зависимости от безразмерных определяющих комплексов и, s и q1 проиллюстрирован на рис. 5. Значение параметра q 1 = 0 соответствует идеальной проводимости материала КС (п0 = 0) и, соответственно, отсутствию джоулева тепловыделения и эффекта термического разупрочнения. При q 1 = 0 с увеличением параметра s сначала происходит возрастание коэффициента формы kfb на момент распада участка КС (рис. 5, а). Это связано с тем, что при небольших значениях s (s < 1) токовое воздействие приводит лишь к уменьшению предельного удли-
Рис. 5. Зависимость коэффициентов формы (а) и скоростного инвертирования (б) на момент разрыва участка кумулятивной струи от определяющих параметров процесса
90 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
нения элемента, не вызывая эффекта дискообразования. В результате элемент КС в момент ее разрыва получается менее вытянутым по сравнению с естественным деформированием, и его форма меньше отклоняется от цилиндрической. При дальнейшем увеличении параметра s начинает проявляться эффект дискообразования (рис. 4, б) и значение kfb уменьшается.
Уменьшение значения kfb при токовом воздействии происходит также с увеличением параметра и (см. рис. 5, а). Это связано с тем, что с ростом значения и увеличивается коэффициент предельного удлинения nfim (11) элементов КС, и, соответственно, до момента разрыва они испытывают токовое воздействие более длительное время и в большей степени наследуют последствия этого воздействия. При фиксированных значениях параметров и и s степень проявления МГД-неустойчивости КС (со снижением коэффициента kfb) повышается при увеличении параметра q1, что обусловлено усилением влияния фактора термического разупрочнения (со снижением прочности облегчается формоизменение участка КС под действием электромагнитных сил). При достаточно больших значениях qi изменяется и вид кривой зависимости коэффициента kfb от параметра s (см. рис. 5, а). Наличие максимума в области малых значений s с ростом qi становится менее выраженным, и в дальнейшем зависимость kfb(s) переходит в монотонно убывающую, либо даже с наличием минимума на начальном участке.
Что касается коэффициента скоростного инвертирования kvb, то на момент распада участков КС в результате ЭДВ он лежит в отрицательной области, и его абсолютное значение с увеличением параметра s монотонно возрастает (см. рис. 5, б). К усилению степени осевого сжатия элементов КС (увеличению |kvb|) ведет также и увеличение безразмерных комплексов и и q1. Причины здесь те же, что приводят и к снижению коэффициента формы kfb с ростом данных параметров.
Для аппроксимации представленных на рис. 5 данных по коэффициентам kfb и kvb могут быть использованы зависимости
kfb = 0,65 +
0,37a1s
1 + ais
0,48a2s 1 + a2 s
; kvb = -3, 4 (u0’5
+ ° 2qi)
0,5 s0’7,
где a1 = 2, 55 + 19u2; a2 = и (3, 2 + 4, 2q1 + 9, 4u2q^). Отметим, что входящие в приведенные соотношения комплексы и, s и q1 (10) должны определяться с использованием параметров участка КС (iz0, R0, aY0, Tr0) на момент начала токового воздействия.
Обработка результатов численного моделирования позволила также получить расчетные соотношения для определения текущих значений коэффициентов формы kf и скоростного инвертирования kv для участка КС в процессе токового воздействия (от момента его начала
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 91
до момента распада участка струи на отдельные элементы):
(12)
где n — текущий коэффициент удлинения участка КС от момента начала ЭДВ на него (1 < n < nfim), определяемый в зависимости от времени t как n = 1 + £z0t.
Изменение коэффициента термического разупрочнения kT (9) участка КС в процессе ЭДВ, как показал анализ результатов численных расчетов, определяется лишь значением безразмерного показателя термического действия тока q (10), причем для описания их взаимосвязи получается зависимость линейного вида
Поскольку предельным значением коэффициента kT является нулевое, то в случае получения с использованием приведенного соотношения отрицательного результата полагается kT = 0.
Построенные зависимости для текущих значений показателей развития МГД-неустойчивости КС позволяют определить состояние участка струи на момент окончания ЭДВ (выхода участка КС из межэлектродного промежутка). Однако, как показывают расчеты, при типичных режимах ЭДВ существенное изменение формы элементов струи с приближением ее к дискообразной происходит не в межэлектродном промежутке, а в силу инерционности материала КС уже после ее выхода из области воздействия, что согласуется и с данными рентгенографических исследований (рис. 2, в). Поэтому для описания эволюции формы элементов КС после прекращения токового воздействия была дополнительно разработана простая модель, удовлетворительно согласующаяся с результатами численного моделирования.
Динамика формоизменения вышедших из межэлектродного промежутка элементов КС может быть приближенно определена на основе задачи об инерционном сжатии цилиндра из несжимаемого жесткопластического материала. Формулировка задачи является следующей. Элемент КС цилиндрической формы длиной le и радиусом Re в начальный момент времени сжимается в осевом направлении со скоростью деформирования £ze. Предполагается, что в процессе осевого сжатия элемента струи скорость осевого деформирования £z остается равномерно распределенной по его длине, и элемент сохраняет цилиндрическую форму с текущим радиусом R и длиной l. Эволюция формы элемента в рамках данной модели может быть описана на основе закона сохранения энергии:
kT(n) = 1 — 1, 5q(n — 1).
(13)
E:
ikin ' ze
+ E
kin
're
E
kin
z
+ Erkin + Epl,
(14)
92 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
где Ekin и E^g” — начальные значения кинетической энергии, осевого и радиального движения материала; Е^*” и Egin — текущие значения данных видов кинетической энергии; Epl — энергия пластического деформирования материала.
Определяя составляющие энергетического баланса с использованием результатов [22], энергетическое соотношение (14) можно переписать в виде:
m^e^ mcR^£2,e mcl2il mcR2s\ mcaYe ^ f l \
24 + 16 = 24 + 16 Р~ П\Те) )
где шс — масса элемента КС. Учитывая выражение для определения текущего значения скорости осевой деформации £z = (dl/dt)/l, из уравнения энергетического баланса получаем дифференциальное уравнение, описывающее динамику инерционного осевого сжатия элемента КС, испытавшего ЭДВ:
d = _l /s2g (l| + (3/2)R2) + 24 (aYe/p) In (l/gj dt ^ l2 + (3/2)R2 '
Полученное уравнение позволяет определить изменение эффективной длины участка КС на этапе его движения после выхода из межэлектродного промежутка до момента взаимодействия с преградой. Входящие в него эффективные характеристики участка струи на момент окончания токового воздействия le, Re, £ze, aYe с использованием (9) определяются соотношениями:
le kfe(1 + £z0te )l0j Re Ru/\/ kfe j (15)
£ze kve ^Ye kTe^Y 0)
где te — продолжительность токового воздействия, а k/e, kve и kTe — значения коэффициентов формы, скоростного инвертирования и термического разупрочнения участка КС в момент его выхода из межэлектродного промежутка, рассчитываемые на основании выражений (12), (13), в которых полагается n = 1 + £z0te.
Таким образом, обобщение результатов численного моделирования МГД-неустойчивости КС при токовом воздействии дает возможность прогнозирования их пробивного действия. После определения на основании полученных соотношений эффективной длины элементов КС, испытавших ЭДВ, глубина их проникания в преграду может быть рассчитана в соответствии с гидродинамической теорией [23].
С использованием представленных результатов была разработана методика определения пробивного действия КЗ в условиях токового воздействия на КС [24]. Геометрические и кинематические параметры КС при этом рассчитывались с помощью методики [25]. В расчетах дополнительно принимался во внимание фактор объемного разрушения материала КС при ее выходе из межэлектродного проме-
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 93
жутка [13, 26, 27]. Условия для объемного разрушения КС возникают при ее выходе из межэлектродного промежутка вследствие резкого исчезновения сжимающего действия электромагнитных сил, что может повлечь за собой радиальное диспергирование термически разу-прочненного материала струи. Реализация данного эффекта ведет к уменьшению средней плотности материала КС, что, в свою очередь, снижает ее пробивное действие. При рассмотрении этого механизма разрушения КС была использована модель равномерно удлиняющегося цилиндрического сжимаемого упругопластического стержня, по которому протекает изменяющийся по определенному закону осевой электрический ток (возможностью развития МГД-неустойчивости в данной модели пренебрегалось) [26]. По результатам численных расчетов было установлено, что скорость радиального рассеивания материала струи можно оценивать на основе простого энергетического соотношения в предположении, что потенциальная энергия объемного сжатия струи магнитным давлением pm (8), после прекращения его действия, расходуется на совершение работы по разрушению материала, “гашение” кинетической энергии радиального сходящегося к оси симметрии движения частиц струи и, напротив, на придание им кинетической энергии радиального разлета. При этом работа разрушения материала струи принималась равной потенциальной энергии его всестороннего растяжения со средним напряжением, определяемым значением предела текучести (снижающимся по мере джоулева нагрева материала). Для определения скорости радиального рассеяния vR поверхности КС после ее выхода из межэлектродного промежутка использовалось соотношение
где K — модуль объемного сжатия материала струи [28, 29].
Из (16) следует, что радиальное рассеяние КС реализуется при интенсивности токового воздействия (магнитном давлении pm), обеспечивающей неотрицательность подкоренного выражения. Пробивное действие объемно разрушенной части КС рассчитывалось на основе гидродинамической теории проникания [23] в сочетании с концепцией критической скорости проникания [30]. Предполагалось, что после выхода из области ЭДВ средняя плотность материала струи непрерывно уменьшается в зависимости от рассчитываемой на основании (16) радиальной скорости поверхности КС, а с уменьшением плотности растет нижний порог скорости струи, необходимой для проникания в преграду.
Факторы объемного разрушения и развития МГД-неустойчивости начинают заметно проявляться при примерно одинаковых параметрах токового импульса. Поэтому в расчетной методике определения про-
94 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
(16)
бивного действия КЗ при ЭДВ учитывалась возможность их совместного проявления. Расчет скорости радиального рассеяния объемно разрушающихся элементов КС проводился с учетом изменения их формы и перераспределения скорости, обусловленных развитием МГД-неустойчивости. Входящие в(16) характеристики элементов КС Re, ize и aye на момент окончания ЭДВ определялись с использованием (15).
Типичные результаты расчетов действия КЗ по стальной преграде в отсутствие токового воздействия на КС и при пропускании по ней мощного импульса электрического тока, полученные с использованием разработанной методики, проиллюстрированы на рис. 6. Приведенные результаты были получены для КЗ диаметром 70 мм. Перед стальной преградой размещался элемент ЭДЗ с межэлектродным расстоянием 5е =30 мм. Данные рис. 6, а соответствуют случаю, когда элемент ЭДЗ не связан с источником электрической энергии (нет токового воздействия на КС). Действие КЗ по преграде, представленное
а б
Рис. 6. Расчетные результаты действия кумулятивного заряда по стальной преграде: а — в отсутствие токового воздействия на кумулятивную струю; б — при пропускании по струе мощного импульса электрического тока
ISSN 0236-3941. Вестник МГТУ им. Н.Э.Баумана. Сер. “Машиностроение” 2014. №3 95
на рис. 6, б, зафиксировано для случая запитки элемента ЭДЗ от емкостного накопителя с энергией 300 кДж. Как следует из сравнения приведенных на рис. 6, а, б данных, токовое воздействие рассмотренной интенсивности позволяет снизить глубину проникания КС в преграду примерно в 3 раза. Характерной чертой мощного ЭДВ на КС является то, что, как свидетельствуют расчеты, остаточное пробитие в преграде в этом случае обеспечивается головным участком струи (см. рис. 6, б). Данное обстоятельство связано с тем, что после замыкания струей межэлектродного промежутка протекающий по ней электрический ток нарастает не мгновенно, а в течение некоторого промежутка времени, зависящего от параметров разрядного контура. В результате головной участок КС испытывает слабое токовое воздействие, не приводящее к его разрушению. Как следует из рис. 6, б, эффект дискообразования не затрагивает головных элементов КС, успевающих проскочить межэлектродный промежуток в начале токового разряда. В то же время интенсивность токового воздействия на последующие участки КС оказывается достаточной для их полного разрушения, так что дополнительного вклада в пробитие эти участки не дают.
На рис. 7 проиллюстрировано сравнение расчетных данных снижения глубины пробития КЗ в стальной преграде при различных режимах токового воздействия на КС с результатами экспериментов, выполненных в Институте гидродинамики им. М.А. Лаврентьева СО РАН [14, 21] (рис. 7, а, б) и РФЯЦ-ВНИИЭФ [8] (рис. 7, в, г). Данные, приведенные на рис. 7, а, б, в, получены для зарядов с диаметром
Рис. 7. Сравнение экспериментальных (светлые диаграммы) и расчетных (затемненные диаграммы) данных снижения пробивного действия кумулятивных зарядов при различных режимах токового воздействия на кумулятивную струю (цифры над диаграммами пробития соответствуют номерам кривых изменения разрядного тока, индексом “0” обозначены диаграммы пробития в отсутствие токового воздействия)
96 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
кумулятивной облицовки 45 мм, на рис. 7, г — для зарядов с диаметром облицовки 100 мм. Ширина межэлектродного промежутка 8е (см. рис. 1) составляла 21 мм для экспериментов на рис. 7, а, б, 30 мм — на рис. 7, в и 45 мм — на рис. 7, г. Как видно из рис. 7, между результатами расчетов и экспериментальными данными наблюдается достаточно близкое соответствие в широком диапазоне изменения параметров токового импульса, пропускаемого по КС.
ЛИТЕРАТУРА
1. Средства поражения и боеприпасы / под ред. В.В. Селиванова. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 984 с.
2. Защита танков / под ред. В.А. Григоряна. М.: Изд-во МГТУ им. Н.Э. Баумана, 2007. 327 с.
3. Бабкин А.В., Ладов С.В., Федоров С.В. Электрическая защита перспективной боевой машины XXI века // Оборонная техника. 2000. № 1-2. С. 19-25.
4. Швецов Г.А., Матросов А.Д. Экспериментальное исследование токовой неустойчивости кумулятивных струй // Труды VII Междунар. конф. по генерации мегагауссных магнитных полей и родственным экспериментам. Саров: ВНИИЭФ, 1997. С. 979-986.
5. Швецов Г.А., Матросов А.Д. Токовая неустойчивость кумулятивных струй // Тез. докл. IV Междунар. конф. “Лаврентьевские чтения по математике, механике и физике”. Новосибирск: ИГ СО РАН, 1995. С. 132.
6. Матросов А.Д., Швецов Г.А. Экспериментальное исследование токовой неустойчивости кумулятивных струй // Прикладная механика и техническая физика. 1996. Т. 37. № 4. С. 9-14.
7. Швецов Г.А., Башкатов Ю.Л., Матросов А.Д. и др. Использование МКГ в экспериментах по токовому разрушению струй кумулятивных зарядов // Труды VII Междунар.конф. по генерации мегагауссных магнитных полей и родственным экспериментам. Саров: ВНИИЭФ, 1997. С. 987-991.
8. Павловский А.И., Пляшкевич Л.Н., Шувалов А.М., Бродский А.Я.Эксперимен-тальные исследования разрушения кумулятивной струи импульсом тока // ЖТФ. 1994. Т. 64. Вып. 2. С. 76-82.
9. Павловский А.И., Пляшкевич Л.Н., Шувалов А.М., Бродский А.Я.Исследование некоторых особенностей процесса разрушения кумулятивной струи в сильноточном режиме // ЖТФ. 1994. Т. 4. Вып. 5. С. 43-48.
10. Бабкин А.В., Кружков В.А., Ладов С.В., Маринин В.М., Федоров С.В. Поведение металлических кумулятивных струй под действием импульса тока // Труды VII Междунар. конф. по генерации мегагауссных магнитных полей и родственным экспериментам. Саров: ВНИИЭФ, 1997. С. 992-997.
11. Pollock C.E. Electromagnetic effect of the natural hydrodynamic instability of stretching, high velocity, metallic jets // Megagauss Magnetic Field Generation and Pulsed Power Applications / Ed. M. Cowan and R.B. Spielman. N.Y.: Nova Sci. Publ., 1994. P. 309-316.
12. Ogorkiewicz R.M. Future tank armors revealed // Janes Int. Defense Review. 1997. No. 5. P. 50-51.
13. Appelgren P, Westerling L., Skoglund M.et al. Radial jet dispersion due to current interaction in an electric armour application // Proc. of the 24th Int. Symp. on Ballistics. New Orleans, USA, 2008. P. 935-943.
14. Shvetsov G.A., Matrosov A.D., Babkin A.V., Ladov S.V., Fedorov S.V Behavior of metallic shaped-charge jets with passage of a pulsed electric current through them // Journal of Applied Mechanics and Technical Physics. 2000. Vol. 41. No. 3. P. 394400.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 97
15. Бабкин А.В., Ладов С.В., Маринин В.М., Федоров С.В. Закономерности растяжения и пластического разрушения металлических кумулятивных струй // Прикладная механика и техническая физика. 1999. Т 40. № 4. С. 25-35.
16. Аркулис Г.Э., Дорогобид В.Г. Теория пластичности. М.: Металлургия, 1987. 352 с.
17. Johnson G.R., Cook W.N. A constitutive model and data for metals subjected to large strains, high rates and high temperatures // Proc. of the 7th Int. Symp. on Ballistics. Hague, Netherlands, 1983. P. 541-547.
18. Кнопфель Г Сверхсильные импульсные магнитные поля. М.: Мир, 1972. 392 с.
19. Chou P.C., Carleone J. The stability of shaped-charge jets // Journal of Applied Physics. 1977. Vol. 48. No. 10. P. 4187-4194.
20. Федоров С.В., Бабкин А.В., Ладов С.В. Развитие магнитогидродинамической неустойчивости на подвергающейся электродинамическому воздействию кумулятивной струе // Оборонная техника. 1998. № 1-2. С. 49-56.
21. Федоров С.В., Бабкин А.В., Ладов С.В., Швецов Г.А., Матросов А.Д. О возможностях управления кумулятивным эффектом взрыва с помощью электромагнитных воздействий // Физика горения и взрыва. 2000. Т. 36. № 6. С. 126-145.
22. Бабкин А.В., Ладов С.В., Маринин В.М., Федоров С.В. Особенности инерционного растяжения кумулятивных струй в свободном полете // Прикладная механика и техническая физика. 1997. Т 38. № 2. С. 3-9.
23. Лаврентьев М.А. Кумулятивный заряд и принципы его работы // Успехи математических наук. 1957. Т. 12. Вып. 4. С. 41-56.
24. Федоров С.В., Бабкин А.В., Ладов С.В., Швецов Г.А., Матросов А.Д.Прогнози-рование пробивной способности металлических кумулятивных струй при пропускании по ним мощного импульса электрического тока // ЖТФ. 2003. Т. 73. Вып. 7. С. 28-36.
25. Маринин В.М., Бабкин А.В., Колпаков В.И.Методика расчета параметров функционирования кумулятивного заряда // Оборонная техника. 1995. № 4. С. 34-39.
26. Бабкин А.В., Кружков В.А., Луговой Э.В., Федоров С.В. Математическое моделирование растяжения кумулятивной струи при пропускании через нее электрического тока // Оборонная техника. 1993. № 9. С. 36-39.
27. Бабкин А.В., Колычев М.Е., Ладов С.В., Федоров С.В. О возможном механизме разрушения кумулятивной струи импульсом тока // Оборонная техника. 1995. № 4. С. 47-54.
28. Федоров С.В., Омельченко Е.Ю. Термические и механические эффекты при пропускании по металлическим кумулятивным струям мощного импульса электрического тока // Сб. докл. VI научн. конф. Волжского регионального центра РАРАН “Современные методы проектирования и отработки ракетноартиллерийского вооружения”. В 2-х томах. Саров: РФЯЦ-ВНИИЭФ, 2010. Т. 2 С. 182-192.
29. Федоров С.В. Об эффекте рассеивания металлических кумулятивных струй при пропускании по ним мощного импульса электрического тока // ЖТФ. 2012. Т. 82. Вып. 10. С. 18-30.
30. Физика взрыва. В 2 Т. / под ред. Л.П. Орленко. М.: Физматлит, 2004. Т. 2. 656 с.
REFERENCES
[1] Selivanov V.V. Sredstva porazhenija i boepripasy [Destruction means and explosives]. Moscow, MGTU im. N.E. Baumana Publ., 2008. 984 p.
[2] Grigorjan V.A. Zashhita tankov [Protection for tanks]. Moscow, MGTU im. N.E. Baumana Publ., 2007. 327 p.
[3] Babkin A.V., Ladov S.V., Fedorov S.V. Electrical protection of advanced technology combat vehicle for the XXI century. Oboronnaja tehnika [Defence Technology ], 2000, no. 1-2, pp. 19-25 (in Russ.).
98 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3
[4] Shvecov G.A., Matrosov A.D. Tr. VII Mezhdunar. Konf. “Generacija megagaussnyh magnitnyh polej i rodstvennye jeksperimenty” [Proc. VII Int. Conf. "Generation of megagauss magnetic fields and the related experiments"], Sarov, VNIIJeF Publ., 1997, pp. 979-986 (in Russ.).
[5] Shvecov G.A., Matrosov A.D. The current instability of shaped-charge jet. Tezisy Dokl. IV Mezhd. Konf. “Lavrent’evskie chtenija po matematike, mehanike i fizike” [Summ. Rep. IV Int. Conf. “Lavrentyev Readings on Mathematics, Mechanics and Physics”], Novosibirsk, IG SO RAN Publ., 1995, p. 132 (in Russ.).
[6] Matrosov A.D., Shvecov G.A. Experimental research of the current instability of shaped-charge jets. Prikl. Mekh. Tekh. Fiz. [J. Appl. Mech. Techn. Phys.], 1996, vol. 37, no. 4, pp. 9-14 (in Russ.).
[7] Shvecov G.A., Bashkatov Ju.L., Matrosov A.D. Using of magnetic cumulative generation (ICG) in experiments on the current destruction of the shaped-charges jets. Tr. VII Mezhdunar. Konf. “Generacija megagaussnyh magnitnyh polej i rodstvennye jeksperimenty” [Proc. VII Int. Conf. “Generation of megagauss magnetic fields and the related experiments”], Sarov, VNIIJeF Publ., 1997, pp. 987-991 (in Russ.).
[8] Pavlovskij A.I., Pljashkevich L.N., Shuvalov A.M., Brodskij A.Ja. Experimental researches for the destruction of the shaped-charge jet by pulsed electric current. Zh. Tekh. Fiz. [Techn. Phys. The Russ. J. Appl. Phys.], 1994, vol. 64, iss. 2, pp. 76-82 (in Russ.).
[9] Pavlovskij A.I., Pljashkevich L.N., Shuvalov A.M., Brodskij A.Ja. Some peculiarities investigation of the destruction process of the shaped-charge jet in a heavy-current mode. Zh. Tekh. Fiz. [Techn. Phys. The Russ. J. Appl. Phys.], 1994, vol. 64, iss. 5, pp. 43-48 (in Russ.).
[10] Babkin A.V., Kruzhkov V.A., Ladov S.V., Marinin V.M., Fedorov S.V. Behavior of metallic shaped-charge jet under the action of the current pulse. Tr. VII Mezhdunar. Konf. “Generacija megagaussnyh magnitnyh polej i rodstvennye jeksperimenty” [Proc. VII Int. Conf. “Generation of megagauss magnetic fields and the related experiments”], Sarov, VNIIJeF Publ., 1997, pp. 992-997 (in Russ.).
[11] Pollock C.E. Electromagnetic effect of the natural hydrodynamic instability of stretching, high velocity, metallic jets. In book “Megagauss Magnetic Field Generation and Pulsed Power Applications”. Ed. Cowan M., Spielman R.B. N.Y., Nova Sci. Publ., 1994, pp. 309-316.
[12] Ogorkiewicz R.M. Future tank armors revealed. Janes Int. Defense Review, 1997, no. 5, pp. 50-51.
[13] Appelgren P., Westerling L., Skoglund M. Radial jet dispersion due to current interaction in an electric armour application. Proc. of the 24th Int. Symp. on Ballistics, New Orleans, USA, 2008, pp. 935-943.
[14] Shvetsov G.A., Matrosov A.D., Babkin A.V., Ladov S.V., Fedorov S.V. Behavior of metallic shaped-charge jets with passage of a pulsed electric current through them. Journal of Applied Mechanics and Technical Physics, 2000, vol. 41, no. 3, pp. 394400.
[15] Babkin A.V., Ladov S.V., Marinin V.M., Fedorov S.V. Regularities of stretching and ductile failure of metallic shaped-charge jets. Prikl. Mekh. Tekh. Fiz. [J. Appl. Mech. Techn. Phys.], 1999, vol. 40, no. 4, pp. 25-35 (in Russ.).
[16] Arkulis G.Je., Dorogobid V.G. Teorija plastichnosti [Plasticity theory]. Moscow, Metallurgija Publ., 1987. 352 p.
[17] Johnson G.R., Cook W.N. A constitutive model and data for metals subjected to large strains, high rates and high temperature. Proc. of the 7th Int. Symp. on Ballistics, Hague, Netherlands, 1983, pp. 541-547.
[18] Knopfel G. Pulsed High Magnetic Fields. Amsterdam, North-Holland Publ. Company, 1970, 373 p. (Russ. Ed.: Knopfel’ G. Sverhsil’nye impul’snye magnitnye polja. Moscow, Mir Publ., 1964. 392 p.)
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 3 99
[19] Chou P.C., Carleone J. The stability of shaped-charge jets. J. Appl. Phys., 1977, vol. 48, no. 10, pp. 4187-4194.
[20] Fedorov S.V., Babkin A.V, Ladov S.V. Formation of magnetohydrodynamic instability of the shaped-charge jet exposing the electrodynamic impact. Oboronnaja tehnika [Defence Technology], 1998, no. 1-2, pp. 49-56 (in Russ.).
[21] Fedorov S.V., Babkin A.V, Ladov S.V, Shvecov G.A., Matrosov A.D. On the control possibility of a shaped-charge effect of explosion using electromagnetic impact. Fiz. Goreniya Vzryva [Combust., Explos., Shock Waves], 2000, vol. 36, no. 6, pp. 126145 (in Russ.).
[22] Babkin A.V., Ladov S.V., Marinin V.M., Fedorov S.V. Inertial stretching features of shaped-charge jets in free flight. Prikl. Mekh. Tekh. Fiz. [J. Appl. Mech. Techn. Phys.], 1997, vol. 38, no. 2, pp. 3-9 (in Russ.).
[23] Lavrent’ev M.A. Shaped-charge and principles of its operation. Usp. Fiz. Nauk [Sov. Phys.-Usp.], 1957, vol. 12, iss. 4, pp. 41-56 (in Russ.).
[24] Fedorov S.V., Babkin A.V, Ladov S.V Shvecov G.A., Matrosov A.D. Penetration prediction of metallic shaped-charge jets by passing through them a powerful pulse of electric current. Zh. Tekh. Fiz. [Techn. Phys. The Russ. J. Appl. Phys.], 2003, vol. 73, iss. 7, pp. 28-36 (in Russ.).
[25] Marinin VM., Babkin A.V., Kolpakov VI. Analysis technique of functioning parameters of the shaped-charge. Oboronnaja tehnika [Defence Technology], 1995, no. 4, pp. 34-39 (in Russ.).
[26] Babkin A.V., Kruzhkov V.A., Lugovoj Je.V, Fedorov S.V Mathematical simulation of stretching of the shaped-charge jet when electric current is passing through it. Oboronnaja tehnika [Defence Technology ], 1993, no. 9, pp. 36-39 (in Russ.).
[27] Babkin A.V., Kolychev M.E., Ladov S.V., Fedorov S.V. On a possible mechanism of destruction of the shaped-charge jet by current pulse. Oboronnaja tehnika [Defence Technology], 1995, no. 4, pp. 47-54 (in Russ.).
[28] Fedorov S.V., Omel’chenko E.Ju. Thermal and mechanical effects when powerful pulse of electric current is passing through metal shaped-charge jets. Sb. Dokl. VI Nauchn. konf. Volzhskogo regional’nogo centra RARAN “Sovremennye metody proektirovanija iotrabotki raketno-artillerijskogo vooruzhenija”. V2-h tomah [Proc. 4th Sc. Conf. of Volga Regional Center of Russian Academy of Rocket and Artillery Sciences (RARAS) “Advanced methods of design and development of missile and artillery weapons”. In two volumes], Sarov, FGUP “RFJaC-VNIIJeF” Publ., 2010, vol. 2, pp. 182-192 (in Russ.).
[29] Fedorov S.V. On the dispersion effect of metallic shaped-charge jets when a powerful pulse of electric current is passing through them. Zh. Tekh. Fiz. [Technical Physics. The Russ. J. Appl. Phys.], 2012, vol. 82, iss. 10, pp. 18-30 (in Russ.).
[30] Orlenko L.P. Fizika vzryva. V 2-h tomah [Physics of explosion]. In 2 volumes. 3rd Ed. Moscow, Fizmatlit Publ., 2004. 656 p. (vol. 2).
Статья поступила в редакцию 13.11.2013
Сергей Владимирович Федоров — старший преподаватель кафедры “Высокоточные летательные аппараты” МГТУ им. Н.Э. Баумана. Автор более 250 научных работ в области физики взрыва и высокоскоростного удара.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
S.V. Fedorov — senior teacher of “High Precision Airborne Devices” department of the Bauman Moscow State Technical University. Author of more than 250 publications in the field of physics of explosion and high velocity impact.
Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
100 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №3