CC BY
86
УДК 004.942
ИССЛЕДОВАНИЕ УДАРНОГО ВЗАИМОДЕЙСТВИЯ ИНДЕНТОРА С ТКАНЕВЫМИ БРОНЕПЛАСТИНАМИ, РАСПОЛОЖЕННЫМИ НА ПЛАСТИЛИНОВОМ ОСНОВАНИИ
H.Ю. Долганина
Рассматривается суперкомпьютерное моделирование динамического взаимодействия индентора и тканевых бронепластин, состоящих из 1 - 5 слоев ткани размером 10х10 см с возможностью разрушения нитей, расположенных на пластилиновом основании. Определены механические характеристики регистрирующей среды (технического пластилина). Проведено сравнение расчетных и экспериментальных данных.
Ключевые слова: суперкомпьютерное моделирование, конечно-элементная .модель, удар, регистрирующая среда, тканевая бронепластина.
Введение
Проверка тканевого бронежилета на прочность в стандартах разных стран - это обстрел бронепластин стандартными пулями контрольного оружия на основании, имитирующем тело человека (специальный пластилин) [1, 2]. Оценка последствий удара пули в бронежилет, надетый на человека, осуществляется путем сопоставления размеров вмятины, оставленной пластиной бронежилета на пластилине, с уровнем травмирования тела человека.
При ударе пули в тканевый бронежилет возникают сложные физические явления, связанные с большими перемещениями, скольжениями, фрикционными контактами, повышением температуры, пластическими деформациями, учесть которые в простых экспериментально-аналитических подходах ранее было невозможно [3]. Для более адекватного описания этого процесса в настоящее время используют численные методы, в частности, МКЭ-пакет LS-DYNA [4-7].
В данной работе проведено сопоставление расчетных (LS-DYNA) и экспериментальных данных по взаимодействию ударника с тканевыми бронепластинами, расположенными на пластилиновом основании.
Статья организована следующим образом. В разделе 1 приведена постановка задачи. В разделе 2 описываются методы исследования, и приводится описание задачи. В разделе 3 обсуждаются результаты проведенных исследований. В заключении суммируются основные результаты, полученные в данной работе.
I. Постановка задачи
Проведено суперкомпьютерное моделирование натурных экспериментов: удар ин-дентором в тканевую бронепластину, расположенную на пластилиновом основании. Определены механические свойства регистрирующей среды (технического пластилина). Проведено сравнение расчетных и экспериментальных данных.
Были рассмотрены бронепластины, состоящие из 1 - 5 слоев ткани полотняного переплетения размером 10х10 см. В работе использована ткань из арамидных нитей Русар. В расчетной модели нити имеют относительную свободу перемещения с возможностью
вытягивания с учетом сухого трения. Рассматривали нити, которые имеют прямоугольное поперечное сечение и были представлены одним оболочечным элементом по ширине с одной точкой интегрирования по толщине и выполнены из ортотропного материала с малыми поперечно - сдвиговыми свойствами. Нити в расчетной модели могли разрушаться. Тканевая бронепластина располагалась на неподвижном пластилиновом основании. В расчете индентор имел сферическую форму диаметром 4,5 мм и массой 0,5 г, с начальной скоростью 100 м/с и был выполнен из абсолютно жесткого материала.
2. Методы исследования
Чтобы снизить размерность задачи при расчете, геометрия нитей была предельно упрощена. Криволинейная ось нити была заменена ломаной, с прямолинейными горизонтальными участками и наклонными (34° к горизонтальному участку нити). Повторяющийся элемент ткани полотняного переплетения показан на рис. 1, где отмечены номера ключевых точек, координаты которых были введены в программу ANSYS, после чего по ключевым точкам были заданы соответствующие поверхности. Далее набор поверхностей был размножен до получения необходимых размеров модели.
Для решения задачи динамического взаимодействия тканевой бронепластины, расположенной на регистрирующей среде и жесткого индентора с помощью пакета программ LS-DYNA геометрия и сетка конечных элементов (рис. 2) (к-файл) были созданы в пакете программ ANSYS с последующей коррекцией.
1 А 79
Рис. 1. Повторяющийся элемент Рис. 2. Сетка конечных элементов
Реальная арамидная нить Русар номинальной линейной плотности 58,8 текс состоит из нескольких сотен элементарных волокон диаметром 12 - 15 мкм, которые имеют крутку 100 м-1 (зарубежные нити, как правило, некрученые). Для материала нити наиболее близкой является модель ортотропного упругого тела, в которой высокая жесткость имеет место в продольном направлении нити и на несколько порядков меньшая в поперечном и сдвиговом направлениях. Нити моделировали оболочечными элементами с одной точкой интегрирования по толщине Белычко-Вонга-Чанга (тип 10) [8].
Принято, что разрушение нити происходит при достижении в ней величины первого главного напряжения сг1 = 3 ГПа. Эта величина соответствует средней прочности исследованных арамидных нитей.
В пакете программ Ь8-БУКЛ представлено более сотни моделей материала [8], наиболее интересными для нас были малопараметрические модели, их анализ показал, что для нитей наиболее подходит материал *MAT_ENHANCED_COMPOSITE_DA-
MAGE, позволяющий задать ортотропные свойства нитей и учесть разрушение при достижении первого главного напряжения равным 3 ГПа (табл. 1).
Таблица 1
Механические характеристики нитей
Параметр Обозначение Величина
Толщина нити, мкм Т 100
Модуль упругости, МПа Ех 1,4-105
Ey ■ О 3
Плотность, кг/м3 Р 1440
Коэффициент Пуассона ^xy 0,3
Модуль сдвига, МПа Gxy ■ О 3
Gyz ■ О 3
Gzx ■ О 3
В конечно-элементной модели ткани нити имеют относительную свободу перемещения и возможность вытягивания с учетом сухого трения. Контакт объектов моделировали командой *CONTACT_AUTOMATГС_SURFACE_TO_SURFACE с коэффициентом сухого трения 0,4, характерным для арамидных нитей [9].
Стальной индентор при взаимодействии с тканевой бронепластиной при исследуемых скоростях до 110 м/с не разрушается и не имеет пластических деформаций, поэтому для сокращения времени расчетов для него был выбран материал *MAT_ЩGID (жесткое тело) с плотностью р = 7800 кг/м3 и модулем упругости Е = 2,1-Ю11 Па.
Для того чтобы построить адекватную модель регистрирующей среды (технического пластилина), необходимо экспериментально определить ее механические свойства. Модуль упругости регистрирующей среды (технического пластилина) был определен из испытаний на квазистатическое сжатие цилиндрических образцов (диаметр 20 мм, высота 40 мм) на машине Щ8ТБ,ОК 5882 со скоростями нагружения 5, 50, 250 и 500 мм/мин (рис. 3). Обработка начальных участков кривых деформирования дала значение Е = 0,7 ± 0,1 МПа. Плотность пластилина р « 1800 кг/м3.
Для определения зависимости предела текучести от скорости деформирования был поставлен эксперимент на низкоскоростной удар (рис. 4). Цилиндрический образец из пластилина, соединенный соосно со стальным грузом-цилиндром, бросали с разных высот Н = 0,5 - 2 м на жесткое неподвижное основание. После удара была измерена остаточная длина образца. Скорость соударения V) была определена из равенства потенциальной и кинетической энергий (тепловые эффекты были незначительными) тдН = шУ02/2, где д - ускорение свободного падения, м/с2; т - масса груза и пластилина, кг.
Движение цилиндра (груза) считали равнозамедленным. С учетом введенных допущений среднюю скорость деформации образца из пластилина можно определить по зависимости: В ~ ^ ^ ■ 2д7 = ^ 1 гДе к) ~ начальная длина пластилинового цилин-
дра, м; / - время, с.
0,10 0,16 0,14 0,12 0,1 0,08
0,06
0,04 0,02
0 Ю 20 30 40 50 60
Деформация, %
-------500 мм/мин;
Скорость захвата _____________50мм/1ин|
Рис. 3. Диаграммы деформирования цилиндрических образцов из пластилина с различными скоростями нагружения при сжатии
Так как в пластилиновом образце при ударе возникают большие пластические деформации, то для определения предела текучести упругим участком диаграммы деформирования пренебрегали (материал - жесткопластический). Предел текучести можно найти из равенства потенциальной энергии работе внутренних сил: тдН = ОтДБ, где От - предел текучести; Б - площадь поперечного сечения цилиндра из пластилина.
070
:г1
Р 048
А
Рис. 4. Схема испытания на сжатие цилиндрических образцов из пластилина:
1 - жесткое неподвижное основание; 2 - цилиндр из пластилина; 3 - стальной цилиндр
По полученным в эксперименте значениям предела текучести От была построена степенная функция зависимости От (МПа) от скорости деформирования £ (1/с) с показателем меньше единицы ar(s) = 0.94(é)015. Температура пластилина при испытаниях была (21 ± 1) °С.
Материал технического пластилина в модели считали упругопластическим с зависимостью предела текучести от скорости деформирования. Из списка материалов, заложенных в библиотеку пакета LS-DYNA, для технического пластилина был выбран *MAT_STRAIN_RATE_DEPENDENT_PLASTICITY, который позволяет учесть зависимость предела текучести от скорости деформирования в табличном виде [8].
3. Результаты исследований и их анализ
Расчеты динамического взаимодействия индентора с тканевыми бронепластинами (полотняное переплетение), состоящих из одного, двух, трех, четырех и пяти слоев ткани размером 10х10 см с возможностью разрушения нитей, расположенных на пластилиновом основании размером 10х10х5 см были проведены на суперкомпьютере «СКИФ -Аврора ЮУрГУ» [10], который используется для вычисления большого количества индустриальных задач [11].
В натурном эксперименте выстрелы производились из пневматического пистолета ИЖ-53 в центр бронепластины под прямым углом. Скорость индентора составляла (100 ± 5) м/с, ее измеряли хронографом S044 с погрешностью 1 м/с. Индентор имел сферическую форму диаметром 4,5 мм и массой 0,5 г, материал - сталь. По каждой бронепластине было сделано не менее 3-х выстрелов. После каждого выстрела были проведены измерения глубины отпечатков, оставленных бронепластиной на основании из пластилина (табл. 2). Различие расчетных и экспериментальных глубин вмятин не превышает 4 %.
Таблица 2
Сравнение расчетных и экспериментальных данных. Г лубина вмятин в регистрирующей
среде после выстрела
Количество слоев в бронепластине Глубина отпечатка, мм Погрешность, %
Среднее экспериментальное значение Расчетное значение
1 10,0 9,6 4
2 6,0 5,8 3
3 5,0 4,9 2
4 4,4 4,4 0
5 5,2 5,2 0
Сравнение расчетных и экспериментальных картин деформирования бронепластины из двух слоев ткани показано на рис. 5.
Как видно из рис. 5 характер деформирования тканевой бронепластины в расчетном и экспериментальном исследованиях одинаков. На рис. 5 видно, что в расчете и
эксперименте происходит вытягивание нитей, которые оказались расположенными под индентором.
Расчет
Эксперимент
Рис. 5. Сравнение расчетных и экспериментальных картин деформирования тканевой
бронепластины (два слоя ткани)
Сравнение расчетных и экспериментальных картин деформирования основания (формы и глубины вмятин) при деформировании бронепластины из двух слоев ткани показано на рис. б.
Форма вмятины, расчет
Форма вмятины, эксперимент
Глубина вмятины, расчет
Глубина вмятины, эксперимент
Рис. б. Сравнение расчетных и экспериментальных картин деформирования основания (формы и глубины вмятин) при деформировании бронепластины из двух слоев ткани
Заключение
Рассмотрено локальное динамическое деформирование многослойной тканевой бронепластины, расположенной на пластилиновом основании (МКЭ, пакет LS-DYNA) и проведено сопоставление с экспериментом. Показано, что форма и размеры вмятин, полученные расчетным методом, отлично согласуются с экспериментальными данными.
Показана возможность использования модели упругопластического тела с пределом текучести, зависящим от скорости нагружения для замены технического пластилина (регистрирующей среды, имитирующей реакцию тела человека на локальный удар).
Работа выполнена при поддержке грантов РФФИ № 12-08-31184 мол_а и № 10-07-96007-р_урал_а.
Литература
1. ГОСТ Р 50744 - 95. Бронеодежда. Классификация и общие технические требования. Прин. Постановлением Госстандарта России от 27.02.95 № 82. Введ. с изм. № 1, 2. Утвержден в сент. 1998 г., мае 2002. - М.: ИПК изд-во стандартов, 2003. -8 с.
2. National Institute of Technology Standard, NIJ Standard 0101.06 Ballistic Resistance of Body Armor. - July, 2008.
3. Бова, В.Г. Концепция построения мягкого бронежилета для максимальной реализации свойств арамидных нитей / В.Г. Бова // Рабочая одежда и средства индивидуальной защиты. - 2001. - № 2(10). - С. 11 - 15.
4. Nilakantan, G. Using LS-DYNA® to Computationally Assess the Уо-Ущо Impact Response of Flexible Fabrics Through Probabilistic Methods / G. Nilakantan, M. Keefe, E.D. Wetzel, T.A. Bogetti, R. Adkinson, J.W. Gillespie // 11th International LS-DYNA® Users Conference: Proceedings of the International Users Conference (06-08 June 2010, Detroit). - 2010. - P. 5-1 - 5-12.
5. Долганина, Н.Ю. Оценка баллистического предела и прогиба многослойных тканевых пластин при ударе индентором / Н.Ю. Долганина // Вестник ЮжноУральского государственного университета. Серия "Машиностроение". - 2010. -Вып. 15. - № 10(186). - С. 17-23.
6. Долганина, Н.Ю. Моделирование ударных процессов в тканевых бронежилетах и теле человека на вычислительном кластере "СКИФ Урал" / Н.Ю. Долганина, С.Б. Сапожников, А.А. Маричева // Вычислительные методы и программирование: Новые вычислительные технологии. - 2010. - Т. 11. - C. 117-126.
7. Долганина, Н.Ю. Проектирование новых конструкций тканевых бронепанелей с использованием суперкомпьютерных вычислений / Н.Ю. Долганина, С.Б. Сапожников // Вестник Южно-Уральского государственного университета. Серия "Математическое моделирование и программирование". - 2011. - Вып. 10. -№ 37(254). - С. 71-81.
8. Hallquist, J.O. LS-DYNA Keyword User’s Manual v.970 / J.O. Hallquist - Livermore Software Technology Corporation, 2003. - 1564 p.
Исследование ударного взаимодействия индентора с тканевыми...
9. Martinez, M.A. Friction and wear behaviour of Kevlar fabrics / M.A. Martinez, C. Navarro, R. Cortes, J. Rodriguez, V. Sanchez-Galvez // Journal of materials science. - 1993. - Vol. 28. - P. 1305-1311.
10. Суперкомпьютер «СКИФ-Аврора ЮУрГУ»: URL: http://supercomputer.susu.ac.ru/ computers/skif_avrora/ (дата обращения: 01.09.2012).
11. Исследование производительности суперкомпьютеров семейства "СКИФ Аврора” на индустриальных задачах/ А.А. Московский, М.П. Перминов, Л.Б. Соколинский и др. // Вестник Южно-Уральского государственного университета. Серия "Математическое моделирование и программирование". - Вып. 6. - 2010. - № 35(211). - С. 66-78.
Наталья Юрьевна Долганина, кандидат технических наук, доцент кафедры системного программирования, начальник отдела поддержки и обучения пользователей Лаборатории суперкомпьютерного моделирования, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), dnu001@mail.ru.
Investigation of Impact Interaction of the Indenter With Fabric Armor Plate Which Located on the Clay Basis
N. Yu. Dolganina, South Ural State University (Chelyabinsk, Russian Federation)
Supercomputer simulation of the dynamic interaction between the indenter and textile armor plate is considered. Armor plate made of 1 - 5 layers of fabric 10x10 cm and located on the clay bottom. The fibers in textiles have the ability to break down. Mechanical characteristics of the backing material (technical clay) are defined. Comparison of the calculated and experimental results conducted.
Keywords: supercomputer simulation, FEA model, impact, backing material, textile armor plate.
References
1. GOST R 50744 - 95. Broneodezhda. Klassifikatsiya i obshchiye tekhnicheskiye trebo-vaniya [Armor Clothes. Classification and General Technical Requirements]. Prin. Postanovleniyem Gosstandarta Rossii ot 27.02.95 № 82. Vved. s izm. № 1, 2. Utver-zhden v sent. 1998 g., maye 2002. Moscow: IPK izd-vo standartov, 2003. 8 p.
2. National Institute of Technology Standard, NIJ Standard 0101.06 Ballistic Resistance of Body Armor. July, 2008.
3. Bova V.G. Kontseptsiya postroyeniya myagkogo bronezhileta dlya maksimalnoy reali-zatsii svoystv aramidnykh nitey [The Concept of Building a Soft Body Armor to Maximize the Properties of Aramid Fibers] Rabochaya odezhda i sredstva individualnoy zashchity [Clothing and personal protective equipment]. 2001. No. 2(10). P. 11-15.
4. Nilakantan G., Keefe M., Wetzel E.D., Bogetti T.A., Adkinson R., Gillespie J.W. Using LS-DYNA® to Computationally Assess the Vo-V1oo Impact Response of Flexible Fabrics Through Probabilistic Methods // 11th International LS-DYNA® Users Conference: Proceedings of the International Users Conference (06-08 June 2010, Detroit). 2010. P. 5-1 - 5-12.
5. Dolganina N.Yu. Otsenka ballisticheskogo predela i progiba mnogosloynykh tkanevykh plastin pri udare indentorom [Evaluation of Ballistic Limit and Multilayer Fabric Plate Deflection Under Indenter Impact] Vestnik Yuzhno-Uralskogo gosudarstvennogo uni-versiteta. Seriya "Mashinostroyeniye" [Bulletin of South Ural State University. Series: Mechanical Engineering]. 2010. No. 10. P. 17-23.
6. Dolganina N.Yu., Sapozhnikov S.B., Maricheva A.A. Modelirovaniye udarnykh protsessov v tkanevykh bronezhiletakh i tele cheloveka na vychislitelnom klastere "SKIF Ural” [Modeling of Impact Processes in a Texile Armor and in a Human Body on a SKIF Ural Computational Cluster] Vychislitelnyye metody i programmirovaniye: Novyye vychislitelnyye tekhnologii [Numerical Methods and Programming: New Computing Technologies] 2010. Vol. 11. P. 117-126.
7. Dolganina N.Yu., Sapozhnikov S.B. Proyektirovaniye novykh konstruktsiy tkanevykh bronepaneley s ispolzovaniyem superkompyuternykh vychisleniy [Design of new Constructions of Textile Armor Panel Using Supercomputing] Vestnik Yuzhno-Uralskogo gosudarstvennogo universiteta. Seriya "Matematicheskoye modelirovaniye i programmi-rovaniye" [Bulletin of South Ural State University. Series: Mathematical Modeling, Programming & Computer Software]. 2011. No. 37(254). Vol. 10. P. 71-81.
8. Hallquist J.O. LS-DYNA Keyword User’s Manual v.970 / Livermore Software Technology Corporation, 2003. 1564 p.
9. Martinez M.A., Navarro C., Cortes R., Rodriguez J., Sanchez-Galvez V. Friction and Wear Behaviour of Kevlar Fabrics // Journal of materials science. 1993. Vol. 28. P. 1305 - 1311.
10. Supercomputer «SKIF-Avrora YuUrGU»: URL: http://supercomputer.susu.ac.ru/
computers/skif_avrora/ (accessed: 01.09.2012).
11. Moskovskiy A.A., Perminov M.P., Sokolinskiy L.B., Cherepennikov V.V., Shamakina A.V. Issledovaniye proizvoditelnosti superkompyuterov semeystva "SKIF Avrora" na industrialnykh zadachakh [Performance Analysis of the Family of Supercomputers "SKIF Aurora" on Industrial Problems] Vestnik Yuzhno-Uralskogo gosudarstvennogo universiteta. Seriya "Matematicheskoye modelirovaniye i programmirovaniye" [Bulletin of South Ural State University. Series: Mathematical Modeling, Programming & Computer Software]. 2010. No. 35(211). P. 66-78.
Поступила в редакцию 5 ноября 2012 г.
