Modeling the stress-strain state of a stamping tool with a wear-resistant coating on sheet metal stamping separation operations using the LS-DYNA software package

Abstract


This study examined the modeling process in the Ansys LS-DYNA software package of the stress-strain state of the working zone of deformation in sheet metal stamping separation operations using a tool with a wear-resistant ion-plasma coating. The object of the study is a working tool (punches and dies) of die tooling for sheet metal stamping processes with a wear-resistant coating based on titanium nitride deposited by ion-plasma spraying and cathodic-ion cleaning methods. The aim of the study was to determine the effect of a wear-resistant coating applied to the working edge of the stamping tool on the stress-strain state in the deformation zone. As part of the study, physical and mathematical models of cutting processes using software packages were proposed, and the stress state of the tool materials in the deformation working zone was evaluated and analyzed. In the modeling process, the finite element simulation of dynamic interactions of solids was used in the Ansys LS-DYNA software package. The following results were obtained: mathematical models of the processes of cutting down the sheet material were constructed (models of the punch, matrix and blanks), a transition was made from 3D-models to simplified coaxial symmetric models that allow one to evaluate the stress-strain state in the working area by separate sections, models of the processes of cutting with different types of tools (with wear-resistant coating and without coating), visual interpretations of the distribution of stress fields in the tool material are obtained, stress graphs are constructed eny and potential energy in the cell surface layer at the critical points of contact on the characteristic process steps of sheet material cutting.

Full Text

В современном машиностроении большое значение имеет стойкость инструментальной оснастки, ее износостойкость и ресурс рабочей поверхности. Одним из факторов, влияющих на износ рабочей поверхности штампового инструмента, является напряженно-деформированное состояние инструментального материала в зоне деформации. Рабочие поверхности штампового инструмента (пуансонов и матриц) работают в специфических условиях: ударные нагрузки, высокие значения коэффициентов трения, что в совокупности обусловливает существенное снижение эксплуатационного ресурса пуансонов и матриц, а также образование неоднородных схем напряженно-деформированного состояния. В Ульяновском государственном техническом университете на кафедрах «Материаловедение и обработка металлов давлением» и «Инновационные технологии в машиностроении» совместно с АО «Ульяновский НИАТ» и АО «Ульяновский патронный завод» проводятся исследования технологии повышения стойкости штампового инструмента с использованием износостойких покрытий, наносимых методами ионно-плазменного напыления [1-10]. Одним из инструментов для изучения особенностей НДС является построение математических моделей с использованием программных средств, позволяющих производить анализ и прогнозирование развития НДС в исследуемой зоне (зоне контакта рабочих поверхностей инструмента и материала заготовки). Сотрудниками кафедр «Материаловедение и обработка металлов давлением» и «Инновационные технологии в машиностроении» совместно с предприятием АО «Ульяновский НИАТ» с целью определения влияния износостойкого покрытия на основе нитрида титана на напряженно-деформированное состояние рабочей зоны инструмента были проведены исследования по построению математической модели процесса вырубки. Для построения математической модели был выбран метод конечных элементов, реализованный в программном пакете Ansys LS-DYNA. LS-DYNA - это многоцелевая программа, предназначенная для анализа нелинейного динамического отклика трехмерных неупругих структур. Программа использует следующие методы: дискретные методы, бессеточные методы [11-17]. Моделирование процесса вырубки осуществлялось поэтапно в следующей последовательности: - построение исходных моделей процессов вырубки-пробивки (оценка НДС); - построение моделей инструмента с износостойким покрытием (оценка НДС); - сравнение НДС моделей без нанесенного износостойкого покрытия и с нанесенным на инструмент износостойким покрытием на основе нитрида титана. На первом этапе была построена конечно-элементная модель процесса вырубки со следующими параметрами: диаметр пуансона 10 мм, диаметр отверстия в матрице 10,05 мм, толщина вырубаемой заготовки 1 мм, межинструментальный зазор 0,05 мм, материал инструмента - инструментальная сталь (сталь марки Х12М), материал заготовки - алюминиевый сплав Д16 (рис. 1). В программе LS-DYNA были заданы следующие параметры процесса вырубки: тип материала - упругопластический, размер конечно-элементной сетки - (0,01×0,01) мм. При задании параметров материала инструмента использовался тип материала инструмента без учета разрушения; пластический тип материала вырубаемой заготовки; применялись справочные и экспериментальные данные, полученные на предыдущем этапе НИР, - кривые течения материала, результаты испытаний на растяжение, модуль упругости, коэффициент Пуассона и др. Рис. 1. Конечно-элементная 3D-модель процесса вырубки: 1 - пуансон (упрощенная модель); 2 - заготовка; 3 - матрица Представленная модель позволяет оценить НДС рабочей зоны, однако вследствие наличия большого количества расчетных элементов и наложенных ограничений требует длительного времени для расчетов (около двух суток) и является достаточно сложной для последующей модификации и корректировки. В то же время из-за большого размера ячеек конечно-элементной сетки данная модель не позволяет с достаточной степенью точности определить НДС в отдельных участках рабочей зоны - размер ячеек сетки превышает толщину покрытия на два порядка, поэтому для уменьшения времени расчета при одновременном уменьшении размеров ячеек конечно-элементной сетки был осуществлен переход к построению двухмерной модели в сечении и наложены соответствующие ограничения (рис. 2). Согласно данным справочной литературы [16, 17], данный переход позволяет повысить точность получаемых значений НДС в области рабочей зоны без повышения времени расчета. Рис. 2. 2D-модель процесса вырубки: 1 - заготовка; 2 - пуансон; 3 - матрица Оптимизированная модель позволяет уменьшить размер ячеек конечно-элементной сетки до значений 0,030 мм, а модель, выполненная в виде сечения, позволяет без внесения большого числа изменений преобразовывать данную модель под конкретные операции и габаритные размеры инструмента, определить НДС в отдельном сечении и распределение полей напряжений и деформаций на протяжении всего процесса штамповки с достаточно высокой точностью (в элементах размером 30×30 мкм) (рис. 3), что позволяет предугадывать образование критических зон и износ инструмента. На рис. 3, а, иллюстрирующем напряженно-деформированное состояние в рабочей зоне деформирования, градацией желтого и красного цветов показаны поля напряжений (по Мизесу) в начале процесса вырубки. Необходимо отметить, что данная модель при необходимости преобразуется в 3D-модель путем поворота вокруг центральной оси на заданный угол и позволяет рассчитывать поля напряжений по трем осям координат, однако в таком случае время расчета увеличивается значительно. а б Рис. 3. Модифицированная 2D-модель процесса вырубки: а - визуальное отображение полей напряжений в процессе вырубки; б - модель, модифицированная с учетом износостойкого покрытия На следующем этапе разработанная модель была скорректирована с учетом наличия износостойкого покрытия на основе нитрида титана толщиной около 30 мкм (см. рис. 3, б). Для этого часть конечно-элементной сетки заданной толщины была отделена от основного материала модели и к ней были приложены следующие параметры и ограничения: тип материала - хрупкий, вид контакта с инструментальной основой - абсолютный, модуль упругости - 300 ГПа. Параметры материала брались из справочной литературы, а также по результатам испытаний режущего инструмента с покрытием, проводимым на кафедре «Инновационные технологии в машиностроении» УлГТУ [1-10, 18, 19]. После добавления в модель износостойкого покрытия для увеличения точности расчета на рабочей кромке был добавлен искусственный радиус скругления R = 0,010 мкм, что обеспечило [19-25] приближение модели к реальному процессу и предотвратило искажение показаний модели в результате возможного появления концентратора напряжений. На рис. 4 представлено напряженно-деформированное состояние рабочей зоны в процессе вырубки (32-й шаг расчета из 110) - инструмент без покрытия демонстрирует более высокие и неравномерные значения поля напряжений. а б Рис. 4. НДС рабочей зоны инструмента в процессе вырубки: а - без износостойкого покрытия; б - с износостойким покрытием На рис. 5 представлены графики изменения напряженно-деформированного состояния в угловом элементе пуансона (в зоне очага деформации) на всем протяжении процесса вырубки (для инструмента без покрытия значение пиковых напряжений в элементе габаритными размерами 30×30 мкм составляет 2,2 ГПа, для упрочненного инструмента - не превышает 1,7 ГПа). Как видно из графических зависимостей и визуальной интерпретации поля напряжений, представленных на рис. 4 и 5, при использовании износостойкого покрытия наблюдается снижение пиков напряжений и сглаживание НДС в рабочей зоне (площадь поля напряжений уменьшается 1,5-2 раза, пиковые значения напряжений снижаются на 20 %). Этот эффект объясняется снижением контактного трения и изменением характера взаимодействия рабочей поверхности инструмента и листовой заготовки под влиянием износостойкого покрытия, в том числе изменением физико-механических характеристик поверхностного слоя, которые происходят в процессе нанесения износостойкого покрытия. Рис. 5. Графики изменения напряжения в угловом элементе на поверхности рабочей кромки пуансона: 1 - инструмент с покрытием нитрид титана; 2 - инструмент без покрытия На основании анализа поведения математической модели и полученных графических зависимостей было установлено, что с использованием износостойкого покрытия НДС в зоне деформации в процессе вырубки изменяется благоприятным образом, полученные модели напряженно-деформированного состояния в рабочей зоне деформации могут позволить как варьировать схемой НДС, тем самым обеспечивая благоприятные условия нагружения, так и способствовать реализации более «мягких» схем нагружения, тем самым снижая уровень растягивающих напряжений в поверхностном слое инструмента, увеличивая адгезионную способность, а следовательно, и износостойкость покрытия. Сделан вывод, что данные модели могут позволить прогнозировать возможные области износа рабочей поверхности инструмента, а следовательно, более точно определять необходимость поверхностного упрочнения конкретных локальных областей, снижая затраты на производственный процесс.

About the authors

O. I Morozov

Ulyanovsk State Technical University

V. P Tabakov

Ulyanovsk State Technical University

V. N Kokorin

Ulyanovsk State Technical University

M. V Ilyushkin

PC ”Ulyanovsk NIAT”

Yu. A Titov

Ulyanovsk State Technical University

D. I Sagitov

Ulyanovsk Institute of Civil Aviation named after the Chief Marshal of Aviation B.P. Bugaev

References

  1. Табаков В.П. Формирование износостойких ионно-плазменных покрытий режущего инструмента. - М.: Машиностроение, 2008. - 311 с.
  2. Хван А.Д., Хван Д.В., Осинцев A.JI. Повышение стойкости инструментальной стали Х12М // Вестник Воронеж. гос. техн. ун-та. - 2012. - Т. 8, № 5. - С. 131-134.
  3. Табаков В.П., Чихранов А.В. Определение механических характеристик износостойких ионно-плазменных покрытий на основе нитрида титана // Известия Самар. науч. центра РАН. - 2010. - Т. 12, № 4. - С. 292-297.
  4. Згалат-Лозинский О.Б., Рагуля А.В., Дуб Н.В. Свойства нанокристаллических материалов на основе нитрида титана // Наноструктурное материаловедение. - 2010. - № 1. - С. 30-38.
  5. Моделирование процесса интенсивной деформации в роликах с использованием программного пакета Ansys LS-DYNA / О.И. Морозов, М.В. Илюшкин, В.П. Табаков, В.Н. Кокорин, В.В. Марковцева, Ю.Н. Гаврилина // Упрочняющие технологии и покрытия. - 2016. - № 12 (144). - С. 43-47.
  6. Повышение стойкости рабочих поверхностей деталей штампов и пресс-форм из теплостойких сталей / О.И. Морозов, В.П. Табаков, В.Н. Кокорин, Ю.А. Титов // Известия Тул. гос. ун-та. Технические науки. - 2017. - Вып. 11, ч. 1. - С. 64-68.
  7. Повышение стойкости рабочих поверхностей деталей штампов при использовании комплексной модификации / О.И. Морозов, В.Н. Кокорин, В.П. Табаков, Д.И. Сагитов, М.В. Илюшкин, Н.А. Ширманов // Известия Тул. гос. ун-та. Технические науки. - 2019. - № 3. - С. 317-326.
  8. Физическая модель структурирования системы «подложка-покрытие» в процессах комплексного модифицирования поверхностного слоя рабочих частей штампов и пресс-форм / О.И. Морозов, В.Н. Кокорин, В.П. Табаков, Д.И. Сагитов, М.В. Илюшкин, Н.А. Ширманов // Наукоемкие технологии в машиностроении. - 2018.- № 11 (89). - С. 35-41.
  9. Исследование процесса комплексной модификации поверхностного слоя рабочих деталей штампов и пресс-форм, работающих в условиях повышенного износа / О.И. Морозов, В.Н. Кокорин, В.П. Табаков, Д.И. Сагитов, М.В. Илюшкин, Н.А. Ширманов // Живучесть и конструкционное материаловедение (ЖивКоМ-2018): материалы 4-й Междунар. науч.-техн. конф., посвящ. 80-летию ИМАШ РАН, г. Москва, 4-6 декабря 2018 г. / Ин-т комп. исслед. - М.; Ижевск, 2018. - С. 181-183.
  10. Физическая модель механической активации поверхностного слоя металлов и сплавов в процессе одноосного сжатия / О.И. Морозов, В.Н. Кокорин, В.П. Табаков, Д.И. Сагитов, М.В. Илюшкин, Н.А. Ширманов // Живучесть и конструкционное материаловедение (ЖивКоМ-2018): материалы 4-й Междунар. науч.-техн. конф., посвящ. 80-летию ИМАШ РАН, г. Москва, 4-6 декабря 2018 г. / Ин-т комп. исслед. - М.; Ижевск, 2018. - С. 183-184.
  11. Чернявский А.О. Метод конечных элементов. Основы практического применения: монография. - М.: Машиностроение, 2003. - 24 с.
  12. Каплун А.Б., Морозов Е.М., Олферьева М.А. Ansys в руках инженера: Практическое руководство. - М.: Едиториал УРСС, 2003. - 272 с.
  13. Чигарев А.В., Кравчук, А.С., Смалюк А.Ф. Ansys для инженеров: справ. пособие. - М.: Машиностроение, 2004. - 512 с.
  14. Курненков А.В., Шурыгин А.Ю. Численное моделирование процессов резания: учеб. пособие / Нижегород. гос. техн. ун-т им. Р.Е. Алексеева. - Н. Новгород, 2017. - Ч. 1. - 184 с.
  15. LS-DYNA analysis for structural mechanics. Predictive engineering / LSTC Livermore Software Technology Corp. - Livermore, 2014. - 115 p.
  16. Криворучко Д.В., Залога В.О., Корбач В.Г. Основы 3D-моделирования процессов механической обработки методом конечных элементов: учеб. пособие. - Сумы: Изд-во СумДУ, 2009. - 208 с.
  17. Илюшкин М.В. Моделирование процессов обработки металлов давлением (осадка цилиндрической заготовки): учеб.-метод. пособие. - Ульяновск: Изд-во УлГУ, 2013. - 112 с.
  18. Повышение стойкости штампового инструмента с износостойким покрытием на формоизменяющих операциях / В.П. Табаков, В.Н. Кокорин, Е.Л. Корняков, О.И. Морозов, А.С. Алешин, Д.И. Сагитов // Известия Тул. гос. ун-та. Технические науки. - 2020. - № 6. - С. 352-358.
  19. Физическая модель структурирования в процессах комплексного модифицирования поверхностного слоя рабочих частей штампов и пресс-форм / О.И. Морозов, Н.В. Мишов, Н.П. Шиллер, А.М. Бузыцкий // МИКМУС-2018: материалы XXX Междунар. инновац. конф. молодых ученых и студентов, г. Москва, 20-23 ноября 2018 г. - М., 2019. - С. 39-41.
  20. Tobias Erhart Review of solid element formulation in LS-DYNA // LS-DYNA Forum 2011, Stuttgart, Germany, 12 October 2011. - Stuttgart, Germany, 2011. - P. 36.
  21. Haufe A., Schweizerhot K., Du Bois P. Review of shell element formulations // Developer Forum, Filderstand, Germany, 24 September 2013. - Filderstand, Germany, 2013. - P. 35.
  22. Maresca G., Milella P.P., Pino G. A critical review of triaxiality based failure criteria. - ANPA - Via V. Brancati, Roma, 1997.
  23. Bradley N. Maker, Xinhai Zhu Input parameters for metal forming simulation using LS-DYNA / LSTC. - 2000. - 12 p.
  24. Муйземнек А.Ю. Описание поведения материалов в системах автоматизированного инженерного анализа: учеб. пособие / ИИЦ ПГУ. - Пенза, 2005. - 152 с.
  25. Смирнов-Аляев Г.А. Сопротивление материалов пластическим деформациям. Краткие основы. - М.; Л.: Машгиз, 1949. - 248 с.

Statistics

Views

Abstract - 100

PDF (Russian) - 55

Refbacks

  • There are currently no refbacks.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies