Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

INFLUENCE OF THE HIGH-STRENGTH COVERING, PUT WITH HVOF SPRAYING, ON FRACTURE OF THE STEEL SUBSTRATE AT IMPACT

Radchenko P.A. 1, 2 Balokhonov R.R. 1 Radchenko A.V. 2
1 Institute of Strength Physics and Materials Science of SB RAS
2 Tomsk State University of Architecture and Building
1321 KB
In work features of deformation and material fracture with a multilayered covering are investigated at high-speed shock interaction. The boundary dynamic task in three-dimensional statement is solved numerically by a method of finite elements. The multilayered covering is set in calculations obviously. When modeling mechanical reaction of a steel substrate, the projectile and a covering intermediate layer from NiAl ratios of Prandtl-Reuss and the condition equation in the form by Mi-Gruneisen are used. The coat layer from WC-Co shows elastic and fragile behavior. For the description of fracture of materials of a multilayered barrier Hoffman´s criterion which considers various durability on tension and compression is used. It is shown that existence of a multilayered covering promotes increase of dynamic strength of a steel basis. It is established that change of amplitude of waves of compression and the unloading, connected with existence of a two-layer covering, nonlinearly depends on a ratio of thickness of coverings.
fracture
deformation
composites with layer
numerical modeling
mechanic of media with structure

Введение

Высокоскоростное (сверхзвуковое) газопламенное напыление (HVOF) является сравнительно недавним дополнением к семейству процессов газотермического напыления. Данная методика по праву считается наиболее современной из технологий напыления. В странах Европы и Северной Америки HVOF практически вытеснило гальванику и методы вакуумного напыления во многих отраслях. Из-за малой скорости частиц при газоплазменном напылении поверхность частиц успевает окислиться, что приводит к низкой плотности формирования покрытий. Увеличение скорости частиц при меньшей их температуре позволило снизить уровень окисленности частиц и повысить плотность порошкового покрытия. В порошковых распылителях HVOF первого и второго поколений использовалось цилиндрическое сопло, тогда как в третьем поколении используются расширяющиеся профильные сопла Лаваля. Такой подход позволяет скорости потока газа превышать 2000м/с, а скорости частиц порошка 800м/с. При HVOF методике достигается высокая адгезия, а пористость покрытия не превышает 1% [8].

В последнее время получено много результатов по упругим и прочностным свойствам многослойных покрытий, полученных методом сверхзвукового газопламенного напыления [2-4]. Но основная часть работ по исследованию покрытий, полученных методом HVOF, посвящена исследованию их напряженно-деформируемого состояния при статических нагрузках. Анализ же деформирования и разрушения подобных сред при динамических нагрузках сложен и требует, как правило, трехмерного подхода и адекватных моделей поведения не только каждого отдельного слоя покрытия, но и подложки.

В данной работе исследуется поведение многослойного покрытия со стальной подложкой при динамическом нагружении конструкции ударником при различных скоростях взаимодействия. Задача решалась численно, методом конечных элементов с использованием явной конечно-разностной схемы Г. Джонсона [5].

Основные уравнения математической модели

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

,

,

.

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

Поведение стального изотропного цилиндра, являющегося ударником, при высокоскоростном ударе описывается упругопластической средой, в которой связь между компонентами тензора скоростей деформации и компонентами девиатора напряжений определяется соотношениями Прандтля–Рейса [6]:

,

где , — модуль сдвига, — параметр, который определяется с помощью условия Мизеса: , при этом — динамический предел текучести (при упругой деформации , а при пластической ). Шаровая часть тензора напряжений (давление) рассчитывается по уравнению Ми–Грюнайзена как функция удельной внутренней энергии и плотности :

,

где — константы материала, и — начальный и текущий удельные объёмы материала соответственно. Данной моделью также описывается поведения стальной подложки и покрытия NiAl.

Поведение материала покрытия из WC-Co не выходит за рамки упруго-хрупкой модели. Компоненты тензора напряжений в материале преграды до момента разрушения удовлетворяют соотношениям обобщенного закона Гука, записанного в терминах скоростей деформаций:

В качестве критерия разрушения материала подложки используется критерий Хоффмана. Данный критерий позволяет использовать различные прочностные характеристики материала на растяжение и сжатие, и, выраженный через скалярные функции компонент тензора напряжений, имеет вид:

,

где – постоянные материала.

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

Как уже отмечалось ранее, численное моделирование проводится методом конечных элементов в трехмерной постановке. Расчет контактных границ осуществляется по алгоритму, предложенному в работе [6], основанному на применении законов сохранения количества движения и момента количества движения к узлам контактных поверхностей.

Для проверки адекватности, как моделей, так и конечно-разностной схемы был сделан ряд сравнений численных и экспериментальных данных, который показал их удовлетворительное согласование [7]. Также исследована сходимость численного решения, которая послужила основой для дискретизации областей конечными элементами. При численном моделировании максимальное количество конечных элементов (тетраэдров), на которые далее разбивается расчетная область, составляет 2,35·106, при этом количество узлов равняется 4,3·105 (рис.1). Шаг сетки отдельных материалов конструкции был выбран таким образом, чтобы масса всех узлов расчетной области была приблизительно равной. Данное условие значительно повышает точность численного решения. В зависимости от начальной скорости и наличия многослойного покрытия на расчет соответствующего начальным данным варианта задачи требуется приблизительно 1,3·104 временных шагов.

Постановка задачи

В работе рассматривается задача влияния двухслойного покрытия на поведение стальной подложки при динамическом нагружении. Стальная подложка из стали марки 316L (0,03% C, 1% Si, 2% Mn, 0,045% P, 0,03% S, 16–18% Cr, 2–3% Mo и 10–14% Ni) толщиной 1,5мм имеет двухслойное покрытие, нанесенное методом HVOF (сверхзвуковое газоплазменное напыление). Первый слой из сплава NiAl (95% Ni и 5% Al), второй слой содержит 87% WC и 13% Co [9]. Данная конструкция нагружается стальным компактным цилиндрическим ударником из стали марки Ст3. Начальная скорость нормального удара составляет 200 и 400м/с. Толщина первого слоя варьируется от 105 до 355нм, второго – соответственно от 350 до 100нм таким образом, чтобы суммарная толщина всей конструкции оставалась постоянной. Для удобства введем коэффициент , где – толщина соответствующего слоя. Механические характеристики материалов представлены в табл. 1.

 

Рис. 1 Схема взаимодействия ударника с мишенью. Расчетная сетка.
Начальный момент времени.

Табл. 1.

 

Сталь 3

316L

NiAl

WC-Co

Плотность, кг/м3

7850

8031

5900

13900

Модуль Юнга, ГПа

204

200,38

169

398

Коэффициент Пуассона

0,3

0,29

0,32

0,25

Предел текучести, МПа

1010

262

1453

Предел прочности, МПа

4500

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

.

Здесь – единичный вектор нормали к поверхности в рассматриваемой точке; и – единичные векторы, касательные к поверхности в этой точке; — вектор силы на площадке с нормалью ; – вектор скорости. Нижние индексы у векторов и обозначают проекции на соответствующие векторы базиса; знак «+» указывает на соответствие значений параметров материала верхней границе контактной поверхности, знак «–» – нижней границе.

Обсуждение результатов

Одним из важнейших показателей является амплитуда волн, возникающих в материале преграды после взаимодействия с ударником. На рис. 2 представлены зависимости значений напряжений sz в центральной точке продольной оси стальной подложки для различных значений коэффициента K. Видно, что для случая начальной скорости 200м/с амплитуда как волны сжатия, так и волны разгрузки минимальна для случаев K=0,28 и К=0,49 (максимальная разность составляет 81%), в то время как для случая начальной скорости 400м/с амплитуда для К=0,28 оказывается максимальной, постепенно уменьшаясь с возрастанием коэффициента К и достигая разности в 23,5%.

 

Рис. 2. Зависимость амплитуды напряжений sz от коэффициента К.
Начальная скорость 200м/с (слева) и 400м/с (справа).

На рис. 3 представлено распределение напряжений sz в момент времени мкс для случая начальной скорости 200м/с. Видно, что картина распределения в данной случае меняется и является менее равномерной, и концентраторы сжимающих напряжений сосредоточены вблизи центральной оси преграды, а также в области, располагающей в месте контакта преграды и боковой поверхности ударника. Остаточная скорость центра масс ударника для случая м/с находится в пределах 39,2м/с () и 43,2м/с ().

 

Рис. 3. Распределение напряжений sz по преграде в плоскости для случаев (слева) и (справа) при мкс. м/с.

На рис. 4 представлено распределение напряжений sz в момент времени мкс для случая начальной скорости 400м/с. Видно, что картина распределения напряжений в материале подложки иная и концентрация напряжений равномерно сосредоточена в области пятна контакта.

 

Рис. 4. Распределение напряжений sz по преграде в плоскости для случаев (слева) и (справа) при мкс. м/с.

Заключение

В работе проведен сравнительный анализ проникающей способности ударника в стальную монолитную преграду и преграду с многослойным NiAl - WC-Co покрытием. Исследования показали, что наличие высокопрочного слоя, полученного методом HVOF, повышает сопротивление разрушению конструкции при ударном нагружении. Повышение прочности не пропорционально соотношению объемов покрытия и преграды, а достигается за счет изменения характера развития волновых процессов и связанного с этим снижения уровня растягивающих напряжений в зоне вероятного откола. В результате происходит более интенсивное торможение и деформирование ударника на начальной стадии процесса, что приводит к уменьшению его проникающей способности. Данный эффект наиболее ярко выражен для низких и средних скоростей ударника.

Исследование выполнено в рамках ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (Государственное соглашение № 14.A18.21.2108 по обобщенной теме «Современные методы исследования микроструктуры и механических свойств перспективных материалов») на оборудовании Центра коллективного пользования «Диагностика структуры и свойств наноматериалов» НИУ «БелГУ».

Рецензенты:

Кайбышев Рустам Оскарович, д.ф.м.н., профессор кафедры «Материаловедения и нанотехнологий» ФГАО ВПО «Белгородский государственный национальный исследовательский университет», руководитель лаборатории механических свойств наноструктурных и жаропрочных материалов ФГАО ВПО «Белгородский государственный национальный исследовательский университет», г.Белгород.

Иванов Олег Николаевич, д.ф.-м.н., директор Центра коллективного пользования «Диагностика структуры и свойств наноматериалов» ФГАО ВПО «Белгородский государственный национальный исследовательский университет», г.Белгород.