CC BY
26
Хасанов А.Р. Динамическое расширение цилиндрической полости в сжимаемой упругопластической среде. Анализ сопротивления среды динамическому прониканию остроконечного ударника // Вестник Пермского национального исследовательского политехнического университета. Механика. 2020. № 1. С. 126-137. БО!: 10.15593/регш.шесЬ/2020.1.10
Khasanov A.R. A dynamic expansion of a cylindrical cavity in a compressible elastic-plastic medium. The analysis of medium resistance to dynamic penetration of a sharp-nosed impactor. PNRPU Mechanics Bulletin, 2020, no. 1, pp. 126-137. DOI: 10.15593/perm.mech/2020.1.10
ВЕСТНИК ПНИПУ. МЕХАНИКА № 1,2020 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
БО1: 10.15593/регш.шесЬ/2020.1.10 УДК 539.3
ДИНАМИЧЕСКОЕ РАСШИРЕНИЕ ЦИЛИНДРИЧЕСКОЙ ПОЛОСТИ В СЖИМАЕМОЙ УПРУГОПЛАСТИЧЕСКОЙ СРЕДЕ. АНАЛИЗ СОПРОТИВЛЕНИЯ СРЕДЫ ДИНАМИЧЕСКОМУ ПРОНИКАНИЮ ОСТРОКОНЕЧНОГО УДАРНИКА
А.Р. Хасанов
Пермский государственный национальный исследовательский университет, Пермь, Россия
АННОТАЦИЯ
Приведено решение задачи о динамическом расширении цилиндрической полости в сжимаемой упругопластической среде с учетом конечных деформаций, нелинейной сжимаемости, зависимости предела текучести от давления. Главная цель настоящего исследования - разработка новой инженерной модели проникания остроконечных ударников в диапазоне средних скоростей удара на основе полученных результатов анализа задачи о расширении цилиндрической полости в полупространстве («цилиндрическая» аналогия).
На базе аналитического подхода получена модель, определяющая сопротивление среды динамическому расширению полости. Ключевые параметры модели зависят от механических свойств среды, для этих зависимостей предложены аппроксимирующие соотношения на основе обработки механических характеристик ряда материалов (некоторые сплавы и грунты).
Для вывода модели динамического проникания использовалось основное допущение А.Я. Сагомоняна о радиальном расширении отверстия: частицы материала преграды движутся в радиальном направлении от поверхности проникающего в преграду ударника. Такое допущение будет оправданным для класса ударников в виде тонких заостренных тел вращения.
На базе введенных допущений получена модель сопротивления среды динамическому прониканию тонкого заостренного тела вращения. Новая модель, помимо «стандартных» прочностной и инерционной составляющих, содержит «присоединенную массу», которая меняется в ходе процесса проникания.
Проведена верификация новой модели проникания с использованием ряда экспериментальных работ по прониканию ударников различной формы в алюминиевые сплавы. Изучено влияние «присоединенной массы» и инерционных сил сопротивления среды прониканию. Получены условия применимости новой модели: модель проникания предлагается использовать при оценке сопротивления сжимаемой среды прониканию тонкого заостренного тела вращения при скоростях удара 200-800 м/с.
©ПНИПУ
О СТАТЬЕ
Получена: 17 ноября 2019 г. Принята: 30 января 2020 г. Опубликована: 30 марта 2020 г.
Ключевые слова:
расширение цилиндрической полости, упругопластическая среда, нелинейная сжимаемость, прикладная теория проникания, металлы, грунты, модель динамического проникания, модель радиального расширения отверстия, диапазон средних скоростей удара, критерий Кулона, тонкий заостренный ударник.
© Хасанов Артур Раисович - ассистент, e-mail: artur_raisovich@rambler.ru, : 0000-0003-0776-6909.
Artur R. Khasanov-Assistant, e-mail: artur_raisovich@rambler.ru, : 0000-0003-0776-6909.
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International
License (CC BY-NC 4.0)
© 0®
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
A DYNAMIC EXPANSION OF A CYLINDRICAL CAVITY IN A COMPRESSIBLE ELASTIC-PLASTIC MEDIUM. THE ANALYSIS OF MEDIUM RESISTANCE TO DYNAMIC PENETRATION OF A SHARP-NOSED IMPACTOR
A.R. Khasanov
Perm State National Research University, Perm, Russian Federation
ABSTRACT
This paper presents the solution of the problem related to the dynamic cylindrical cavity expansion in a compressible elastic-plastic medium. Finite strains, nonlinear compressibility and dependence of the yield stress versus pressure are taken into account in the problem formulation. The study targets at developing a new engineering model on the penetration of a sharp-nosed impactor in the range of middle impact velocities based on the problem analysis results of the cylindrical cavity expansion in a half-space (cylindrical cavity expansion approximation).
Based on the analytical approach a model is obtained that determines the resistance of the medium to dynamic cavity expansion. The main parameters of the model depend on the mechanical properties of the medium. For these dependences we proposed approximating relations based on manipulation of the mechanical properties of a number of materials (some alloys and soils).
To derive the dynamic penetration model the A.Ya. Sagomonyan assumption of the radial expansion of the hole is used. It is assumed that particles of the medium material move in a radial direction from the surface of the impactor penetrating into the shield. Such assumption can be applied for the class of impactors in the form of slender sharp-nosed bodies of revolution.
Based on the assumptions we obtained a model of the medium resistance to the dynamic penetration of a slender sharp-nosed body of revolution. The new model, in addition to the "standard" strength and inertial components, contains the "attached mass", which changes during the penetration process.
The experimental validation of the new penetration model using a series of experimental studies on the penetration of various forms of impactors into aluminum alloys is considered. The influence of the "attached mass" and inertial forces of the medium resistance to the penetration is estimated. The conditions of applicability of the new model are obtained: the penetration model is applicable for estimation of the resistance of a compressible medium to penetration of a thin sharp-nosed body of revolution at impact velocities of 200-800 m/s.
©PNRPU
ARTICLE INFO
Received: 17 November 2019 Accepted: 30 January 2020 Published: 30 March 2020
Keywords:
cylindrical cavity expansion, elastic-plastic medium, nonlinear compressibility, applied theory of penetration, metals, soils, dynamic penetration model, hole radial expansion model, middle impact velocities range, Mohr-Coulomb criterion, slender sharp-nosed impactor.
Введение
Настоящая работа является логическим продолжением исследования [1], в котором была поставлена и решена задача о статическом расширении цилиндрической полости в сжимаемой упругопластической среде. Также подчеркнем, что аналогичная задача о расширении сферической полости была рассмотрена в работах [2, 3].
Заметим, что вывод моделей, описывающих расширение полости в среде, - один из эффективных способов оценки сопротивления среды прониканию. Развитию этого подхода посвящено большое число работ, отметим классические работы: исследования советских [4-6] и зарубежных ученых [7-10]; обширные монографии [11, 12]; исследования последних лет, среди которых выделяются работы нижегородских [13-21] и израильских специалистов [22, 23], также отметим обзорные работы [24, 25] и монографию [26], которые содержат важные выводы о связи моделей локального взаимодействия с моделями, описывающими расширение полости. Более подробно обзорная часть представлена в работе [1].
В настоящей статье рассмотрена задача о динамическом расширении цилиндрической полости в сжимаемой упругопластической среде с учетом конечных де-
формаций, нелинейной сжимаемости, зависимости предела текучести от давления. Получена модель сопротивления среды динамическому прониканию тонкого заостренного тела вращения, проведена верификация модели на базе ряда экспериментальных результатов, взятых из справочника [27].
1. Постановка задачи и построение аналитического решения
Приведенная в работе [1] постановка задачи о статическом расширении цилиндрической полости в полупространстве содержит основные соотношения, которые также используются в настоящей работе, отличие возникает только для уравнения движения, в которое включена инерционная составляющая. Поэтому здесь ограничимся перечислением основных формул и уравнений.
Рассматривается пространственно-одномерная задача (с осевой симметрией) о динамическом расширении цилиндрической полости в полупространстве под действием нормальной распределенной нагрузки ст0 на границе. Вводятся следующие обозначения: г, х - ла-гранжевая и эйлеровая координата; Яо, Я - начальный
и конечный радиус полости (рис. 1). Расширение цилиндрического отверстия в материале осуществляется от начального Я0 до конечного радиуса Я, в области Я < х, < х(г) расположена зона пластических деформаций, за границей х > х(г„) - упругая зона. Кроме того, справедливы соотношения:
„ „ д u д x
= -дцт, x(r,t) = r + u(r, t) ,
d v dt
(1)
где u(r, t) - перемещение частицы с координатой r.
Рис. 1. Схема расширения полости: 1 - полость; 2 - область пластических деформаций; 3 - область упругих деформаций
Fig. 1. The scheme of the cavity expansion: 1 is a cavity; 2 is a plastic regime; 3 is anelastic regime
Закон сохранения массы выделенного материального объема pdV = р0d¥0 принимается в форме [28]
Р0 Р
X Ш
r )дг
(2)
где р0, р - начальная и текущая плотность; dV0, dV -начальный и текущий элементарный объем частицы.
Связь между главными компонентами лагранжевого (отмечены символом «*») и эйлерового тензоров напряжений [28] записывается в виде:
* x * дх °r =ar~x , = а) дХ .
Уравнение движения примем в форме [5, 28]
d v
да
дх
Ро r— = х^+( стг-а) дх.
dt
дr
(3)
(4)
В упругой области девиаторные компоненты тензора напряжений в главных осях определяются из соотношений закона Гука в дифференциальной форме:
^ = + — 5Р1 + + = 5Р (5)
5/ [ 5/ 3р ) '5/ 5/ 5/ р '
где О - модуль сдвига; , е, - главные компоненты девиатора тензора напряжений и деформаций.
Условие пластичности принимается в форме Кулона [4]:
х = c+atg9,
(6)
где т, с - касательное и нормальное напряжение на площадке; с - сцепление; ф - угол внутреннего трения. В компонентах тензора напряжений аг и а0 для
пространственно-одномерного случая с осевой симметрией условие (6) принимает вид [4]
ar-0e = ±то +Kar+°e) ,
(7)
где т0 = 2c cos ф, ц = sin ф, знак перед т0 - «минус», если а < а, в противном случае «плюс».
Уравнение состояния запишем в единой форме (для металлов и пористых сред):
р = 3K 'е- 3 K" £
S-S„
(8)
где КК" - коэффициенты (объемные модули) при линейном и нелинейном слагаемом.
Начальный (касательный) объемный модуль материала имеет вид
к=1 др
3 д£
£ = О
= к" + к" .
(9)
Выражение (8), записанное в общем виде, позволяет переходить к частным случаям.
Вариант А. Поскольку металлы являются пластически несжимаемыми, то К" ^ 0 и из (8) следует известное уравнение состояния (рис. 2, кривая 1):
р = 3K £ .
(10)
Рис. 2. Зависимость уравнения состояния: 1 - для металла; 2 - для пористого материала
Fig. 2. The dependence for the equation of state in the metal (1) and porous material (2)
Вариант Б. Для описания сжимаемости пористой среды (K' ^ 0) необходимо учесть экспериментальный факт существования некоторой степени уплотнения р0 / р , начиная с которой сопротивление уплотнению
существенно возрастает. Этот эффект связан с выборкой начальной пористости материала. В этом случае уравнение состояния принимаем в форме [2]:
р = -3K£ р
(11)
£
£
£ — £
где б - предельная деформация пористой среды (рис. 2, кривая 2).
Существует предельный объем V , связанный с предельной деформацией зависимостью
(V - V)/ V = Збр . Переходя в этой формуле от объемов
к плотностям, получим связь между предельной деформацией грунта б и предельной плотностью р в виде
соотношения р = р0 / (1 + Зб ).
Используя второе соотношение (5), можно получить связь между текущей плотностью, предельной деформацией и давлением для пористых сред:
Р = -3к б р
1п р0 - 1п р 1п ро - 1п р- Зб р
и металлов:
р = К (ро/ р-1).
(12)
(13)
Предполагаем отсутствие деформации в направлении оси цилиндрической полости б2 = 0, поэтому в упругой области справедлива связь ст2=у(стг+сте). Также полагаем, что это соотношение сохраняется и в пластической области для у = 1/2 [29]. Используя условие пластичности (7) в условиях цилиндрической симметрии, можно выразить все компоненты тензора напряжений в пластической области через давление:
X X
ст,. = (1 + ц)р, ст0= (1 -ц)р + х0 . (14)
Полагая, что т0 и ц - постоянные величины, преобразуем уравнение (4) к виду
х |Р + р % = г ^. (15)
дг 1 + ц дг 1 + цдг 1 + ц 3/
» ц
Обозначим ц =-, умножим уравнение (15) на
1 + ц
х2ц 1 и проинтегрируем по лагранжевой координате от Ro до г с учетом граничных условий р(Я) = —р0 и х(Я) = Я. Получим выражение для распределения давления в пластической области при Я < г < г (или
Я < х < х,):
Х0 Г _ъХ0 VЯ , ро —2ц' г 2ц*—1 д2и
р = 2ГГ +2цхяJ +ТТТ {гх , (16)
где х = х(г) - внешняя граница пластической области. Необходимо отметить, что в выражение (16) входит
интеграл (обозначим его I (г )= | гх2ц —1 ёг), кото-
рый отражает влияние сил инерции.
При вычислении интеграла будем опираться на предположение о постоянстве плотности (это предположение применяется только при оценке сил инерции) р ~ р0, откуда следует:
х2 = г2 + Я2 - Я2.
Дифференцируем два раза последнее выражение, получим:
2х-х = 2 ЯЯ,
2х-х+2х = 2Я-Я+2Я2.
Вьфазим х:
Я-Я+Я2-\ 1-^-
(17)
С учетом (17) преобразуем интеграл I (г) и перейдем к интегрированию по х:
I (х) = }х
Я
2ц*-1
¿X .
(18)
Вычислим интеграл (18):
2ц
Предполагаем, что во внешней (упругой) области при г > г (или х > х,) справедливо обычное упругое статическое решение при малых деформациях. Общее решение упругой задачи, следуя [29]:
Б
Б
Б
= 2Е г
(20)
где D - неизвестная постоянная.
Следуя методике [2, 3] получения решения задачи о расширении сферической полости, требуем на границе упругой и пластической областей [29] непрерывности напряжений стг, ст0 и радиального перемещения u, откуда находим значение неизвестной D и значение
ро( х) = т^
2 ц
-1 +
+ я2ц' I (х). 1 + ц ** *
Подставляя последнее выражение в (16) и используя связь (18), получим:
р <«=2ц(1
1
1 + ц
1 + у 2ц'
^ ')-фг (!)" (^ --)
(21)
где ^ =
р„Я2 Ж - Я х _ р -> У = К= — Л = — , Р = ~-
^ о Л Л** Л*л< Х/->
X = X
ст = -
0
2
2
г
г
>
Ккаапоу А.Я. / РЖРи ЫесЬапс Би1Шп 1 (2020) 126-137
Аналогичным образом определяется удельное сопротивление ст0 = —стг (Я) материала расширению полости:
2ц
—1 + (1 + ц)
+ Р0 Я 2ц I (х), (22)
где ст0 - нормальное напряжение на границе полости (внешнее воздействие).
Связь между параметрами х. и г имеет вид
х =г 11 + ^ 4Е
Чтобы получить соответствующие (21), (22) выражения для металлов, осуществим предельный переход при ц ^ 0. Воспользовавшись правилом Лопиталя, получим:
р (|) = 1П £ + 0
о0=^-т01пЯ + р0Я2
(1 + у) 1П1К2 [1 — -2
(23)
']п1 + 1(Л2-1)"
Я 21 ;
+ р0ЯЯ1п —. (24) Я
Из соотношения (12), выразив р / р и подставив его в уравнение неразрывности, получим окончательное нелинейное уравнение для определения радиуса цилиндрической полости Я в пористой среде с заданными характеристиками под действием внутреннего давления ст0:
(
|ехр
3ер ■ Р (х) р ( х ) + 3 К е
л
xdx =1 (г2 — Я2). (25)
Проводя аналогичные выкладки для соотношения (13), получим нелинейное уравнение для определения радиуса цилиндрической полости в металле:
х* Кх('1х 1 I 2 -А
=1 (г*2 — Я2).
я р (х) + 3К 2 4 '
(26)
Выражение для а0 (22) и (24) позволяет оценить сопротивление среды динамическому расширению полости:
(27)
где первое слагаемое стс отражает собственное сопротивление, второе и третье слагаемые характеризуют инерционное сопротивление. Безразмерные параметры определяются соотношениями: - для металлов:
_ <5п 1
То 2
= 1 — 1ПЯ, к = к +1 (Я2 —1), к = (28) 2 1 2 2У ' 2 Я
- для пористых сред (грунтов):
= 2Г7С1 -Р + Ц)**' ],
2ц
я-2ц+2 — 1 , _я-2ц — 1
к — к ж , к — "
2ц — 2
2ц
Оценка параметра Я, входящего в соотношения (28) и (29), получена в [1], где предложена следующая аппроксимирующая зависимость:
Я = с
4К'
(30)
где К = К / т0 . В работе [1] также был проведен комплекс расчетов для ряда сплавов и геологических сред на основе большого числа данных и получены следующие параметры аппроксимирующей зависимости (30): С = 1,300 (магниевые сплавы); С = 1,435 (алюминиевые сплавы); С = 1,359 (титановые сплавы); С = 1,283 (стали); С = 1,319 (пески); С = 1,313 (супеси); С = 1,460 (суглинки); С = 1,888 (глины). Отметим, что относительная ошибка аппроксимации для всех зависимостей не превосходит 1 %.
2. Модель динамического проникания
Построим инженерную модель динамического проникания жесткого заостренного тела в пластичную преграду при условии вязкого образования кратера. В качестве основного допущения используем известную модель А.Я. Сагомоняна [5] - модель радиального расширения отверстия, справедливую для заостренных ударников.
Предполагаем, что частицы материала преграды движутся в радиальном направлении в процессе проникания и образования кратера.
Будем предполагать, что радиус кратера Я (у) ограничен поверхностью ударника, за исключением лицевых слоев преграды, где наблюдается инерционное расширение отверстия. Пренебрегая последним, при средних скоростях проникания (при ударе по нормали) примем радиус кратера равным границе поверхности ударника / (С):
Я(у) = / (С) ,
(31)
где /- уравнение образующей ударника; С - координата, отсчитываемая от кончика ударника вдоль его оси (рис. 3).
Уравнение образующей кратера зависит от текущего положения ударника при проникании, поэтому в зависимости (31) удобно перейти к координате у, отсчитываемой от лицевой поверхности, и к текущей глубине проникания ударника Ь: у = Ь — С. Тогда (31) перейдет к виду
Я( у) = / (Ь—у).
(32)
с
Т
0
я
ас =
Рис. 3. Схема проникания ударника в преграду Fig. 3. The scheme of the impactor penetration into the shield
Отметим, что для цилиндрического ударника с конической головной частью высотой h и углом полураствора конуса при вершине а уравнение образующей имеет вид
f ( L - У ) =
f( L - y )• tga, (L - y)< h, I h • tga, (L - y) > h.
(33)
Найдем производные первого и второго порядка от выражения (32):
к=Г(1-у)-1=Г(1-у)-у, (34)
к=Г{ь-у\Ё+г{ь-ууь= = /"(Ь - у). V2 + 1г{Ь - у)-(35) Отметим, что в общем случае (см. рис. 3)
/'(Ь - у) = . (36)
Для ударника с конической головной частью с учетом (33) получим:
- 1 d(v2) 2 dL
(37)
Формулы (37) справедливы при условии (Ь - у) < к ,
что соответствует скорости и ускорению расширения радиуса кратера для конической части ударника. Цилиндрическая часть ударника (согласно допущению) не приводит к расширению кратера, и производная в этом случае будет равна нулю.
Используем основное предположение модели А.Я. Сагомоняна о радиальном расширении частиц материала преграды при проникании жесткого заостренного ударника, а именно: полагаем, что сопротивление материала ри образованию кратера (см. рис. 3) совпадает с удельным сопротивлением расширению цилиндрической полости стд, определяемым формулой (27):
В этом случае нормальное на границе удельное сопротивление прониканию для цилиндрического ударника с конической головной частью согласно (37): - для конической части
2 2 1 т „ d (v2
Pn = CTC + ^lPoV tg a + "k2P0tgaR—ГГ 2 dL
(39)
- для цилиндрической части р„ =стс.
Сопротивление прониканию ударника, обусловленное удельным сопротивлением материала р ,
и трением -
L
Fc = 2л[ f (С) pn (Qtga(Qd С:
FT = 2n|f (СК (QdС .
(40)
(41)
Тогда уравнение движения ударника массы m примет вид
1 d (v2) — m——1
2 dL
= -Fc - Ft = -2rcjf (С)Pn (C)[tga(C) + К ]dС . (42)
В выражении (42) принято, что хл = , где -коэффициент трения.
Уравнение движения (42) с учетом соотношений (34), (35), (27) можно привести к виду (43), по форме совпадающему с известным уравнением, основанным на эмпирической двухчленной зависимости сопротивления прониканию [30]:
1
d (v2)
2 <m + ^ L
= -B(L) - C(L)v2.
(43)
Принципиальное отличие уравнения (43) от эмпирической зависимости состоит в наличии «присоединенной массы» A(L), которая меняется в ходе процесса проникания. В результате этого получаем уравнение движения материальной точки с переменной массой.
В уравнении (43)
L
A(L) = f (C)[tga(C) + кТ ]k2PoRf'(C)dС , (44)
0
L
B(L) = f (C)[tga(C) + kT]CTC(C)dС , (45)
0
L
C( L) = 2л| f (Q[tga(Q + kT ]x
0
x[ kP0Cf' (С))2 + k2P0 Rf "(С)] ¿С
Pn =СТД.
(38)
(46)
Следуя [23, 26], можно вывести формулы в явном виде для «интегральных характеристик» проникания, а именно - баллистического предела (БП) ^ :
0
0
0
(47)
и остаточной скорости (ОС) v :
Vi =1 v2 -
/2 2 \ (Vo - Vbi ) exP -i
4 ' П
2 • C (L)
o m
+ A(L)
dL
(48)
где ^ - начальная скорость (скорость удара); к - высота головной (конической или оживальной) части ударника; Ь - толщина плиты.
Под баллистическим пределом уы понимается такая скорость удара, при которой ударник пробивает преграду конечной толщины Ь и покидает ее с нулевой остаточной скоростью V = 0.
Приведем также более простую модель проникания конического ударника без учета влияния присоединенной массы и трения на сопротивление прониканию:
d(V2) _ 4^2 а
dL
m
хЦС^с (С)d С + ^2а} С • к (С) • Р0 (С)d С^. (49)
В случае проникания в однородную пластину в соотношении (49) можно вычислить интегралы и получить следующую зависимость:
d(v2) _ 4rctg2a dL m
L2
(ctc + ^1Potg2 av 2) —
(50)
где выражение в круглых скобках соответствует удельному сопротивлению проникания (27) при нулевом третьем слагаемом.
3. Верификация модели проникания на основе натурных экспериментов
Для верификации модели из справочника [27] были выбраны пять экспериментальных работ по прониканию ударников с конической и оживальной головной частью в алюминиевую плиту [8-10, 31, 32]. Выбор экспериментов для верификации был продиктован, во-первых, соответствием условий проникания допущениям модели (тонкий заостренный ударник, средние скорости проникания и т.д.), а во-вторых, доступностью информации по механическим характеристикам испытываемых материалов, в частности, по динамическому пределу текучести. Отметим, что в рамках данной статьи верификация модели на геологических средах не рассматривается, т.е. формулы (29) не использовались, все расчеты проведены с применением формул (28).
Еще одно замечание касается учета трения. Для оценки коэффициента трения использовался результат экспериментальной работы [33]. Основной вывод данной работы: при ударном проникании влияние трения на расход кинетической энергии ударника составляет 3 %. Поскольку расход начальной кинетической энергии ударника равен работе силы сопротивления преграды прониканию ударника на конечную глубину по мере его движения в полубесконечной плите, то на базе экспериментальных данных можно оценить коэффициент трения.
Рис. 4. Зависимость скорости ударника от глубины проникания, расчет на основе экспериментальных данных [32]: 1 - без переменной массы; 2 - с переменной массой; 3 - глубина проникания (эксперимент)
Fig. 4. The dependence of the impactor velocity on the penetration depth, the calculation is based on the experimental data [32]: 1 - without variable mass; 2 - with variable mass; 3 - penetration depth (experiment line)
V
bl
Результаты вычислений по данным экспериментальных работ [8-10, 31, 32] составили ^ = 0,01...0,02.
В упомянутых экспериментальных работах в качестве материала преграды использовались следующие алюминиевые сплавы: сплав 6061-Т651 в [8], сплав 5083-Н131 в [9], сплав 7075-Т651 в [10], сплав 6061-Т6511 в [31] и сплав 65-8 в [32]. Для этих сплавов в расчетах применялись следующие значения динамического предела текучести: 500 МПа для сплавов 6061-Т651 и 6061-Т6511 (согласно работам [34, 35]), 700 МПа для сплава 7075-Т651 (согласно исследованию [36]) и 430 МПа для сплава 5083-Н131 (согласно [37]). В перечисленных работах приведены эксперименты и численные расчеты по оценке динамического предела текучести указанных сплавов.
К сожалению, информация по динамическому пределу текучести для сплава марки 65-8 отсутствует, однако работа [32] представляет интерес, поскольку в ней варьировались углы при вершине конической части ударника. Поэтому динамический предел текучести подбирался для самого острого ударника (а = 10°). Таким способом оценивался диапазон углов полураствора конуса, при котором справедлива используемая модель проникания. Также был оценен вклад «присоединенной массы», образуемой в процессе движения ударника в преграде. Исходя из графиков на рис. 4, можно заключить, что влияние «присоединенной массы» увеличивается с ростом угла полураствора конуса, однако остается незначительным - разница между двумя кривыми на графиках составляет несколько процентов. Учитывая основное допущение о радиальном расширении отверстия, справедливое для заостренных ударников, рассматривать большие углы при вершине конуса не имеет смысла, поэтому вкладом «присоединенной массы» можно пренебречь, полагая А(Ь) = 0. Этот результат подтверждается и по экспериментальным данным других работ. Отметим, что влияние «присоединенной массы» может оказаться значительным для затупленных форм ударника, в таком случае в процессе проникания возможно образование пробки, однако эти варианты проникания выходят за рамки настоящей статьи.
В табл. 1 содержатся результаты сравнения теоретических и экспериментальных значений глубины проникания (ГП). Ошибка отклонения теоретических результатов от экспериментальных вычислялась по формуле
L - L
teor_ex
L
•100%,
(51)
где - ГП, полученная в эксперименте; - ГП, рассчитанная по формулам (28), (30), (43)-(46), коэффициент трения в расчетах принят ^ = 0,01. Согласно полученным результатам относительная ошибка не превосходит 10 % для ударников с углом полураствора не
выше 20°. Отметим, что при а = 25° ошибка примерно
равна 10 %, а при а = 30° уже превосходит 20 %. Предлагается следующее условие применимости модели для конических ударников: угол полураствора не выше 25°.
Таблица 1
Сравнение результатов расчета глубины проникания конического ударника в алюминиевый сплав 65-S с экспериментальными данными [32]
Table 1
Comparing the calculated and measured penetration depths of a conical impactor into the 65-S aluminum shield [32]
Угол а, град. Скорость удара, м/с ГП (экспериментальная), мм ГП (теоретическая), мм ю, %
10 650 105,2 103,6 -1,5
762 >137,2 132,6 -3,3
15 644 94,7 92,6 -2,2
761 127,0 119,8 -5,7
20 661 92,7 87,8 -5,3
760 118,6 107,6 -9,3
25 654 85,3 78,2 -8,3
768 111,5 97,2 -12,8
30 670 84,6 72,3 -14,5
793 113,8 88,9 -21,9
40 657 76,2 54,8 -28,1
771 99,6 64,5 -35,3
50 646 71,4 39,9 -44,2
752 90,9 45,2 -50,1
60 632 67,1 26,9 -60,0
774 92,2 30,7 -66,7
В табл. 2 приведены результаты испытаний, расчетов и их сравнение по глубине проникания (ГП) в полубесконечную преграду оживальных ударников с отношением радиуса дуги окружности оживальной формы к максимальному диаметру (калибру), равным трем. Коэффициент трения в расчетах кТ принят равным 0,02. Как и ранее, относительная ошибка вычисляется по формуле (51).
На рис. 5, 6 сравниваются остаточная скорость (ОС) и баллистический предел (БП) для цилиндрического ударника диаметром 8,31 мм с конической головной
частью высотой 14,83 мм (а ® 15,7°) на базе работы [9]. Здесь маркеры - экспериментальные результаты, а кривые - теоретические БП и ОС, вычисленные по формулам (47) и (48) соответственно.
Относительная ошибка вычисления значений ОС и БП, представленных на рис. 5, 6, определялась по формуле, аналогичной (51). В целом можно отметить, что, за исключением случаев сравнения по остаточной скорости для скоростей удара, близких к баллистическому пределу, ошибка отклонения лежит в пределах 10 %.
Отметим, что сопоставление результатов проводилось по значениям, которые удовлетворяли допущениям модели, т.е. высокие скорости проникания опускались. Сравнение с результатами других экспериментальных работ подтвердило непригодность модели в случае про-
ю =
никания конических ударников с углами полураствора более 30°. Также модель непригодна при низких ударных скоростях (менее 200 м/с) и при проникании в тонкие преграды.
Таблица 2
Сравнение результатов расчета глубины проникания оживального ударника в алюминиевые сплавы с экспериментальными данными: 6061-T651 [8], 7075-T651 [10], 6061-T6511 [31]
Table 2
Comparing the calculated and measured penetration depths of an ogive impactor into an aluminum shield: 6061-T651 [8], 7075-T651 [10], 6061-T6511 [31]
Источник Скорость удара, м/с ГП (экспериментальная), мм ГП (теоретическая), мм ю, %
[8] 440 44 43 -3,2
730 94 91 -3,2
790 105 102 -2,8
[10] 372 26 28 5,7
695 70 70 0,4
[31] 569 58 56 -3,9
570 55 56 1,6
679 72 73 1,0
821 102 96 -5,6
794 103 93 -9,4
Рис. 5. Зависимость баллистического предела от толщины преграды, расчет на основе экспериментальных данных [9]
Fig. 5. The dependence of the ballistic limit velocity on the shield thickness, the calculation is based on the experimental data [9]
Сравнительный анализ полученной модели проникания с данными натурных экспериментов позволяет
заключить, что новую модель можно применять для тонких конических ударников (и их аналогов в форме оживала) с углом полураствора не более 25°, проникающих в металлические плиты средней толщины и более, при скоростях проникания 200-800 м/с.
уге5, м/с 800
600
400
200
/ 2 __
0
200 300 400
Vi), м/с
500
600
700 800
Рис. 6. Зависимость остаточной скорости ударника от скорости удара, расчет на основе экспериментальных данных [9]: 1 - преграда толщиной 12,7 мм; 2 - преграда толщиной 50,8 мм; 3 - преграда толщиной 76,2 мм
Fig. 6. The dependence of the residual velocity of the impactor on the impactor velocity, the calculation is based on the experimental data [9]: 1 - the shield thickness is 12,7 mm; 2 - the shield thickness is 50,8 mm; 3 - the shield thickness is 76,2 mm
Заключение
Представлена новая постановка задачи о динамическом расширении цилиндрической полости в упруго-пластической сжимаемой среде с учетом конечных деформаций. Приведено разрешающее выражение в виде интегрального уравнения, определяющего сопротивление динамическому расширению полости и размер пластической зоны для различных материалов. Предложены аппроксимирующие соотношения для определения основных параметров, характеризующих сопротивление среды динамическому расширению полости. С учетом полученного решения разработана новая модель проникания при основном допущении радиального расширения полости, предложенном А.Я. Сагомоняном. Новая модель применима при оценке сопротивления сжимаемой среды прониканию тонкого заостренного тела при скоростях удара 200-800 м/с.
Библиографический список
1. Аптуков В.Н., Хасанов А.Р. Расширение цилиндрической полости в сжимаемой упругопластической среде // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 1. - С. 5-23. Б01: 10.15593/регш.шесЬ/2017.1.01
2. Аптуков В.Н. Расширение сферической полости в упругопластической среде при конечных деформациях. Сообщение 1. Влияние механических характеристик, свободной поверхности, слойности // Проблемы прочности. - 1991. -№ 12. - С. 1262-1268.
3. Аптуков В.Н. Расширение сферической полости в упругопластической среде при конечных деформациях. Сообщение 2. Влияние инерционных характеристик. Температурные эффекты // Проблемы прочности. - 1991. - № 12. -С. 1269-1274.
4. Рахматулин Х.А., Сагомонян А.Я., Алексеев H.A. Вопросы динамики грунтов. - М.: Изд-во Моск. ун-та, 1964. - 237 с.
5. Сагомонян А.Я. Проникание. - М.: Изд-во МГУ, 1974.
6. Сагомонян А.Я. Динамика пробивания преград. - М.: Изд-во Моск. ун-та, 1988.
7. Bishop R.F., Hill R., Mott N.F. The theory of indentation and hardness tests // Proceedings of the Physics Society. - 1945. -Vol. 57, № 3. - Р. 147-159.
8. Forrestal M.J., Okajima K., Luk V.K. Penetration of 6061-T651 aluminum target with rigid long rods // Journal of Applied Mechanics. - 1988. - Vol. 55(4). - Р. 755-760.
9. Forrestal M.J., Luk V.K., Brar N.S. Perforation of aluminum armor plates with conical-nose projectiles // Mechanics of Materials. - 1990. - Vol. 10(1-2). - Р. 97-105.
10. Penetration of 7075-T651 aluminum targets with ogival-nose rods / M.J. Forrestal, V.K. Luk, Z. Rosenberg, N.S. Brar // Int. J. of Solids and Structures. - 1992. - Vol. 29(14-15). -Р. 1729-1736.
11. Satapathy S. Application of cavity expansion analysis to penetration problems / The Univ. of Texas at Austin, Inst. Adv. Technol., 1997. Report No. IAT.R-0136.
12. Yu H.S. Cavity Expansion Methods in Geomechanics // Kluwer Acad. Publ., Dordrecht. - 2000.
13. Котов В.Л. Исследование применимости автомодельного решения задачи о расширении сферической полости в сжимаемой среде для определения давления на поверхности контакта «ударник - грунт» // Проблемы прочности и пластичности: межвуз. сб. - Н. Новгород: Изд-во ННГУ, 2008. -Вып. 70. - С. 123-131.
14. Анализ приближенных решений задачи о расширении сферической полости в грунтовой среде / В.Л. Котов [и др.] // Проблемы прочности и пластичности: межвуз. сб. - Н. Новгород, 2011. - Вып. 73. - С. 58-63.
15. Решение задачи о расширении сферической полости в грунтовой среде в предположении несжимаемости за фронтом ударной волны / Е.Ю. Линник [и др.] // Проблемы прочности и пластичности: межвуз. сб. - Н. Новгород, 2012. - Вып. 74. -С. 49-57.
16. Котов В.Л., Линник Е.Ю., Тарасова А.А. Определение параметров квадратичной модели локального взаимодействия при внедрении сферического ударника в мягкий грунт // Проблемы прочности и пластичности: межвуз. сб. - Н. Новгород, 2013. - Вып. 75(1). - С. 47-55.
17. Котов В.Л., Константинов А.Ю. Численное моделирование плоскопараллельного движения конических ударников в грунтовой среде на основе модели локального взаимодействия // Вычислительная механика сплошных сред. - 2014. -Т. 7, № 3. - С. 225-233. DOI: 10.7242/1999-6691/2014.7.3.22
18. Применение модели локального взаимодействия для определения силы сопротивления внедрению ударников в песчаный грунт / В.Л. Котов [и др.] // Прикладная механика и техническая физика. - 2013. - Т. 54, № 4.- С. 114-125.
19. Линник Е.Ю., Котов В.Л., Константинов А.Ю. Моделирование процессов динамического внедрения пространственных тел в сжимаемую упругопластическую среду // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 4. - С. 92108. DOI: 10.15593/perm.mech/2017.4.07
20. Сравнительный анализ методов моделирования проникания и плоскопараллельного движения конических ударников в грунтовой среде / В.Г. Баженов [и др.] // Прикладная механика и техническая физика. - 2015. - Т. 56, № 3. - С. 4454. DOI: 10.15 372/PMTF20150306
21. Баженов В.Г., Котов В.Л., Линник Е.Ю. Методика численного расчета оптимальных форм тел вращения при движении в грунтовой среде // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. -№ 2. - С. 5-20. DOI: 10.15593/perm.mech/2015.2.01
22. Ben-Dor G., Dubinsky A., Elperin T. Analytical Solution for Penetration by Rigid Conical Impactors Using Cavity Expansion Models // Mech. Res. Commun. - 2000. - Vol. 27(2). -Р. 185-189. DOI: 10.1016/S0093-6413(00)00080-X
23. Ben-Dor G., Dubinsky A., Elperin T. Analysis of Ballistic Properties of Layered Targets Using Cavity Expansion Model // Int J. Fract. - 1998. - Vol. 90(4). - Р. 63-67.
24. Ben-Dor G., Dubinsky A., Elperin T. Analytical engineering models of high speed normal impact by hard projectiles on metal shields // Cent. Eur. J. Eng. - 2013. - Vol. 3(3). - Р. 349373. DOI: 10.2478/s13531-013-0108-7
25. Ben-Dor G., Dubinsky A., Elperin T. Engineering models of high speed penetration into geological shields // Cent. Eur. J. Eng. -2014. - Vol. 4(1). - Р. 1-19. DOI: 10.2478/s13531-013-0135-4
26. Ben-Dor G., Dubinsky A., Elperin T. High-Speed Penetration Dynamics // Engineering Models and Methods. - 2013.
27. Ben-Dor G., Dubinsky A., Elperin T. World Scientific Handbook of Experimental Results on High Speed Penetration into Metals, Concrete and Soils // World Scientific. - 2016.
28. Нигматулин Р.И. Динамика многофазных сред. Ч. 1. -М.: Наука, 1987. - 464 с.
29. Качанов Л.М. Основы теории пластичности. - М.: Наука, 1969. - 420 с.
30. Витман Ф.Ф., Златин Н.А., Иоффе Б.С. Сопротивление деформированию металлов при скоростях 10-6-102 м/с // Журнал технической физики. - 1949. - Т. XIX, вып. 3. - С. 123-128.
31. Penetration of 6061-T6511 aluminum targets by ogive-nose steel projectiles with striking velocities between 0.5 and 3.0 km/s / A.J. Piekutowski, M.J. Forrestal, K.L. Poormon, T.L. Warren // Int. J. of Impact Engineering. - 1999. - Vol. 23 (1, part 2). -Р. 723-744. DOI: 10.1016/S0734-743X(99)00117-7
32. Brooks P.N. On the prediction of crater profiles produced in ductile targets by the impact of rigid penetrators at ballistic velocities. Rep. DREV R-686/73. Defence Research Establishment Valcartier, Quebec, Canada, 1973.
33. Krafft J.M. Surface Friction in ballistic penetration // J. Appl. Phys. - 1955. - Vol. 26, no. 10. - P. 1248-1253.
34. Yadav S., Chichili D.R., Ramesh K.T. The mechanical response of a 6061-T6 A1/AkO3 metal matrix composite at high rates of deformation // Acta Metallurgica et Materialia. - 1995. -Vol. 43(12). - P. 4453-4464. DOI: 10.1016/0956-7151(95)00123-D
35. Banerjee B., Bhawalkar A.S. An extended mechanical threshold stress plasticity model: modeling 6061-T6 aluminum alloy // J. of Mechanics of Materials and Structures. - 2008. - Vol. 3(3). - P. 391-424. DOI: 10.2140/jomms.2008.3.391
36. High strain rate tensile behavior of aluminum alloy 7075 T651 and IS 2062 mild steel / R. Pothnis Jayaram, Y. Perla, H. Arya, N.K. Naik // ASME Journal of Engineering Materials and Technology. -2011. - Vol. 133(2): 021026. - P. 1-9. DOI: 10.1115/1.4003113
37. Rosenberg Z., Dekel E. Revisiting the perforation of ductile plates by sharp-nosed rigid projectiles // Int. J. of Solids and Structures. - 2010. - Vol. 47(22-23). - P. 3022-3033. DOI: 10.1016/j.ijsolstr.2010.07.003
References
1. Aptukov V.N., Khasanov A.R. Rasshirenie tsilindricheskoi polosti v szhimaemoi uprugoplasticheskoi srede [Expansion of a cylindrical cavity in a compressible elastic-plastic medium]. PNRPU Mechanics Bulletin, 2017, no. 1, pp. 5-23.
2. Aptukov V.N. Rasshirenie sfericheskoi polosti v uprugoplasticheskoi srede pri konechnykh deformatsiiakh. Soobshchenie
1. Vliianie mekhanicheskikh kharakteristik, svobodnoi poverkh-nosti, sloinosti [Expansion of a spherical cavity in a compressible elasto-plastic medium. I. The influence of mechanical characteristics, free surface, and lamination]. Problemy prochnosti, 1991, no. 12, pp. 1262-1268.
3. Aptukov V.N. Rasshirenie sfericheskoi polosti v uprugo-plasticheskoi srede pri konechnykh deformatsiiakh. Soobshchenie
2. Vliianie inertsionnykh kharakteristik. Temperaturnye effekty [Expansion of a spherical cavity in a compressible elasto-plastic medium. II. Effect of inertial forces. Temperature effects]. Problemy prochnosti, 1991, no. 12, pp. 1269-1274.
4. Rakhmatulin K.A., Sagomonyan A.Ya., Alekseev N.A. Vo-prosy dinamiki gruntov [Soils Dynamics]. Moscow, Moscow State University Press, 1964, 237 p.
5. Sagomonyan A.Ya. Pronikanie [Penetration]. Moscow, Moscow State University Press, 1974.
6. Sagomonyan A.Ya. Dinamika probivaniia pregrad [Dynamics of Barriers Perforation]. Moscow, Moscow State University Press, 1988.
7. Bishop R.F., Hill R., Mott N.F. The theory of indentation and hardness tests. Proceedings of the Physics Society, 1945, vol. 57, no. 3, pp. 147-159.
8. Forrestal M.J., Okajima K., Luk V.K. Penetration of 6061-T651 aluminum target with rigid long rods. Journal of Applied Mechanics, 1988, vol. 55(4), pp. 755-760.
9. Forrestal M.J., Luk V.K., Brar N.S. Perforation of aluminum armor plates with conical-nose projectiles. Mechanics of Materials, 1990, vol. 10(1-2), pp. 97-105.
10. Forrestal M.J., Luk V.K., Rosenberg Z., Brar N.S. Penetration of 7075-T651 aluminum targets with ogival-nose rods. Int. J. of Solids and Structures, 1992, vol. 29(14-15), pp. 1729-1736.
11. Satapathy S. Application of cavity expansion analysis to penetration problems. The Univ. of Texas at Austin, Inst. Adv. Technol., 1997, Report No. IAT.R-0136.
12. Yu H.S. Cavity Expansion Methods in Geomechanics. Kluwer Acad. Publ., Dordrecht, 2000.
13. Kotov V.L. Issledovanie primenimosti avtomodel'nogo resheniia zadachi o rasshirenii sfericheskoi polosti v szhimaemoi srede dlia opredeleniia davleniia na poverkhnosti kontakta «udar-nik-grunt» [Investigation of applicability of self-similarity analysis to problem of expansion of a spherical cavity in a compressible medium for determining pressure on striker-soil contact surface]. Problemy prochnosti i plastichnosti, 2008, vol. 70, pp. 123-131.
14. Kotov V.L., Linnik E.Yu., Makarova A.A., Tarasova A.A. Analiz priblizhennykh reshenii zadachi o rasshirenii sfericheskoi polosti v gruntovoi srede [The analysis of approximate solutions of a spherical cavity expansion in a soil medium problem] Problemy prochnosti iplastichnosti, 2011, vol. 73, pp. 58-63.
15. Linnik E.Yu., Kotov V.L., Tarasova A.A., Gonik E.G. Reshenie zadachi o rasshirenii sfericheskoi polosti v gruntovoi srede v predpolozhenii neszhimaemosti za frontom udarnoi volny [The solution of the problem of the expansion of a spherical cavity in a soil medium assuming incompressibility beyond the shock font]. Problemy prochnosti i plastichnosti, 2012, vol. 74, pp. 49-57.
16. Kotov V.L., Linnik E.Yu., Tarasova A.A. Opredelenie parametrov kvadratichnoi modeli lokal'nogo vzaimodeistviia pri vnedrenii sfericheskogo udarnika v miagkii grunt [Determining the parameters of a quadratic model of local interaction for a spherical striker penetrating into soft soil]. Problemy prochnosti i plastichnosti, 2013, vol. 75(1), pp. 47-55.
17. Kotov V.L., Konstantinov A.Yu. Chislennoe modelirova-nie ploskoparallel'nogo dvizheniia konicheskikh udarnikov v grun-tovoi srede na osnove modeli lokal'nogo vzaimodeistviia [Numerical modeling of the planeparallel motion of conical projectiles in the soil environment based on a local interaction model]. Computational Continuum Mechanics. 2014, vol. 7, no 3, pp. 225-233.
18. Kotov V.L. Primenenie modeli lokal'nogo vzaimodeistviia dlia opredeleniia sily soprotivleniia vnedreniiu udarnikov v peschanyi grunt [Using a local-interaction model to determine the resistance to penetration of projectiles into sandy soil]. J. of Applied Mechanics and Technical Physics, 2013, vol. 54, no 4, pp. 114-125.
19. Linnik E.Yu., Kotov V.L., Konstantinov A.Yu. Modeling the processes of dynamic penetration of dimensional bodies in a compressible elastoplastic medium. PNRPU Mechanics Bulletin, 2017, no. 4, pp. 92-108.
20. Bazhenov V.G. Sravnitel'nyi analiz metodov modeliro-vaniia pronikaniia i ploskoparallel'nogo dvizheniia konicheskikh udarnikov v gruntovoi srede [Comparative Analysis of Methods for Modeling the Penetration and Plane-Parallel Motion of Conical Projectiles in Soil]. Journal of Applied Mechanics and Technical Physics, 2015, vol. 56, no. 3, pp. 44-54.
21. Bazhenov V.G., Kotov V.L., Linnik E.Yu. Metodika chis-lennogo rascheta optimal'nykh form tel vrashcheniia pri dvizhenii v gruntovoi srede [Method of numerical calculation of optimal forms of bodies of revolution at movement in soil medium]. PNRPU Mechanics Bulletin, 2015, no. 2, pp. 5-20.
22. Ben-Dor G., Dubinsky A., Elperin T. Analytical solution for penetration by rigid conical impactors using cavity expansion models. Mech. Res. Commun., 2000, vol. 27(2), pp. 185-189.
23. Ben-Dor G., Dubinsky A., Elperin T. Analysis of ballistic properties of layered targets using cavity expansion model. Int J. Fract., 1998, vol. 90(4), pp. 63-67.
24. Ben-Dor G., Dubinsky A., Elperin T. Analytical engineering models of high speed normal impact by hard projectiles on metal shields. Cent. Eur. J. Eng., 2013, vol. 3(3), pp. 349-373.
25. Ben-Dor G., Dubinsky A., Elperin T. Engineering models of high speed penetration into geological shields. Cent. Eur. J. Eng., 2013, vol. 4(1), pp. 1-19.
26. Ben-Dor G., Dubinsky A., Elperin T. High-Speed Penetration Dynamics. Engineering Models and Methods, 2013.
27. Ben-Dor G., Dubinsky A., Elperin T. World Scientific Handbook of Experimental Results on High Speed Penetration into Metals, Concrete and Soils. World Scientific, 2016.
28. Nigmatulin R.I. Dinamika mnogofaznykh sred. Chast' 1 [Dynamics of Multiphase Media. Part 1]. Moscow, Nauka, 1992, 464 p.
29. Kachanov L.M. Osnovy teorii plastichnosti [Fundamentals of the theory of plasticity]. Moscow, Nauka, 1969, 420 p.
30. Vitman F.F., Zlatin N.A., Ioffe B.S. Soprotivlenie defor-mirovaniju metallov pri skorostjah udara 10-6-102 m/s [The resistance of metals to deformation at impact velocities of 10-6-102 m/s]. Zhurnal tekhnicheskoi fiziki, 1949, vol. XIX., no. 3, pp. 123-128.
31. Piekutowski A.J., Forrestal M.J., Poormon K.L., Warren T.L. Penetration of 6061-T6511 aluminum targets by ogive-nose steel projectiles with striking velocities between 0.5 and 3.0 km/s. Int. J. of Impact Engineering, 1999, vol. 23(1, part 2), 723-744.
32. Brooks P.N. On the prediction of crater profiles produced in ductile targets by the impact of rigid penetrators at ballistic velocities. Rep. DREV R-686/73. Defence Research Establishment Valcartier, Quebec, Canada, 1973.
33. Krafft J.M. Surface Friction in ballistic penetration. J. Appl.Phys, 1955, vol.26, no. 10, pp. 1248-1253.
34. Yadav S., Chichili D.R., Ramesh K.T. The mechanical response of a 6061-T6 A1/AkO3 metal matrix composite at high
rates of deformation. Acta Metallurgica et Materialia, 1995, vol. 43(12), pp. 4453-4464.
35. Banerjee B., Bhawalkar A.S. An extended mechanical threshold stress plasticity model: modeling 6061-T6 aluminum alloy. J. of Mechanics of Materials and Structures, 2008, vol. 3(3), pp. 391-424.
36. Pothnis Jayaram R., Perla Y., Arya H., Naik N.K. High strain rate tensile behavior of aluminum alloy 7075 T651 and IS 2062 mild steel. ASME Journal of Engineering Materials and Technology, 2011, vol. 133(2): 021026, pp. 1-9.
37. Rosenberg Z., Dekel E. Revisiting the perforation of ductile plates by sharp-nosed rigid projectiles. Int. J. of Solids and Structures, 2010, vol. 47(22-23), pp. 3022-3033.
