CC BY
13
УДК 621.771.23
DOI: 10.14529/met170410
ИССЛЕДОВАНИЕ МАКРО- И МИКРОМЕХАНИКИ ДЕФОРМИРОВАНИЯ СТАЛИ В ХОДЕ ПРОЦЕССА ВОЛОЧЕНИЯ НА ОСНОВЕ МУЛЬТИМАСШТАБНОГО МОДЕЛИРОВАНИЯ
1 12 Д.В. Константинов', А.Г. Корчунов', О.П. Ширяев2,
А.С. Кузнецова1, Д.Г. Емалеева1
1 Магнитогорский государственный технический университет им. Г.И. Носова, г. Магнитогорск, Россия,
2 ОАО «Магнитогорский метизно-калибровочный завод «ММК-МЕТИЗ», г. Магнитогорск, Россия
Описан метод мультимасштабного моделирования, позволяющий учесть микроструктурное строение обрабатываемой стали и динамические структурно-фазовые превращения. На примере необходимости учета микроструктурного строения TRIP сталей сформулированы актуальные проблемы традиционных методов проектирования технологических режимов с использованием компьютерного конечно-элементного моделирования. Апробация предложенного метода моделирования проводилась на примере внедрения стали TRIP 700 в традиционный процесс волочения и оценки возможности управления распределением свойств посредством параметров процесса. Метод также позволил изучить деформационное взаимодействие элементов микроструктуры TRIP стали между собой, тем самым объяснив наличие более высоких (в 2-3 раза превышающих средние значения) радиальных деформаций в пластичных фазах вблизи большего скопления зерен более прочных фаз бейнита и мартенсита.
Ключевые слова: волочение, сталь, трип-эффект, мультимасштабное моделирование, репрезентативный объем, напряженно-деформированное состояние, микроструктура.
Введение
Современная машиностроительная индустрия с каждым годом предъявляет все более высокие требования к механическим и эксплуатационным свойствам деталей машин и агрегатов. Элементы механизмов и машин должны обладать не только достаточным ресурсом надежности, но и быть в состоянии адаптироваться под возможные неблагоприятные условия эксплуатации.
Основными технологическими процессами производства заготовок для производства деталей и элементов агрегатов являются методы обработки металлов давлением (ОМД). Получение новых свойств и технологических характеристик заготовок становится возможным за счет либо расширения номенклатуры используемых при ОМД материалов, либо разработки новых режимов или способов ОМД [1, 2]. В обоих случаях возникает необходимость в современных инженерных методах проектирования, обеспечивающих предельно высокую достоверность прогноза результатов исследуемого технологического процесса.
Наименее затратным методом проектирования технологических процессов с позиции
необходимых ресурсов является компьютерное моделирование с использованием коммерческих программных продуктов, предоставляющих инженеру широкий спектр возможностей по решению задач анализа и оптимизации. В целом это позволяет значительно снизить затраты на тестирование результатов исследуемых технологических режимов и технологий.
Большинство компьютерных моделей в современных исследованиях процессов ОМД не учитывают особенности структурного строения обрабатываемого металла. Как правило, проектирование режимов обработки осуществляется с позиции изотропного материала.
Однако известно, что главным «индикатором» технологических преобразований в обрабатываемом металле является его микроструктурная композиция. Параметры микроструктуры металла в совокупности наиболее полно отражают результаты деформационных, термических и комбинированных воздействий [3-5].
В современных металлургических технологиях стали активно использоваться материалы со сложной, метастабильной микро-
структурой: например, стали классов TRIP (Transformation-Induced Plasticity) и TWIP (Twinning-Induced Plasticity) [6-11]. Динамические структурно-фазовые превращения, которые происходят под действием деформации в этих сталях, обеспечивают большую свободу инженерам при выборе дизайна, оптимизации веса и общей технологии производства. Поэтому моделирование обработки давлением данных классов сталей без учета поведения их микроструктуры под нагрузкой в принципе невозможно.
По этой причине развитие методов компьютерного моделирования с позиции повышения их прогнозирующей точности должно осуществляться в направлении репрезентации микроструктурного строения обрабатываемого металла и напряженно-деформированного состояния (НДС) на его микроуровне.
В то же время переход на микроуровень в компьютерных моделях приводит к неизбежному уменьшению размера расчетной сетки и, как следствие, росту количества элементов, требуемых расчетных ресурсов и времени. Поэтому имеющиеся на данный момент методы моделирования в своей реализации не всегда являются технически доступными и оперативными.
Следовательно, совершенствование подходов компьютерного моделирования с учетом особенностей поведения микроструктуры металла при исследовании процессов ОМД является современным актуальным направлением в инженерном анализе. При этом используемые методы должны быть оптимизированы с позиции расчетного времени и необходимых калькуляционных ресурсов.
Исходные данные
Одним из базовых процессов ОМД для производства длинномерных металлических изделий, используемых в качестве передельной заготовки, является волочение калиброванной стали круглого сечения. В ходе этого процесса горячекатаная круглая заготовка относительно большого диаметра (от 6 до 40-50 мм) подвергается однократному волочению через монолитную волоку с обжатиями, как правило, 5-15 %.
Данный процесс широко используется для получения промежуточных заготовок для стержневых деталей (валы, оси) в машиностроении, приборостроении, сельхозмашиностроении. С технологической точки зрения
процесс волочения предоставляет широкие возможности по управлению характеристиками обрабатываемой калиброванной стали. Геометрические параметры волоки (угол рабочего канала волоки, длина калибрующей зоны) в совокупности со степенью обжатия и варьируемой скоростью волочения позволяют управлять глубиной проработки сечения металла, скоростью его деформации, качеством поверхности, точностью размеров.
Поэтому расширение номенклатуры используемых в данном процессе материалов за счет более технологичных сталей классов TRIP и TWIP потенциально позволит получать традиционные виды изделий с расширенным комплексом технологических и эксплуатационных характеристик.
В частности, одними из наиболее широко используемых сталей в современном машиностроении являются стали TRIP 600-700. Главной отличительной особенностью данных сталей является так называемый трип-эффект, который заключается в способности остаточного аустенита в их микроструктуре аккумулировать энергию и превращаться в более твердый мартенсит под действием высокоскоростной деформации.
В ходе эксплуатации стержневой детали из TRIP стали она будет обладать свойством поверхностного упрочнения в условиях вибраций и ударных нагрузок, так как остаточный аустенит в поверхностных слоях будет трансформироваться в более прочный мартенсит. При использовании элементов механизмов из таких сталей в условиях абразивного износа будет реализовываться свойство сверхвязкости, которое позволяет за счет формирования зоны упрочнения уменьшить и замедлить отрыв частиц металла с поверхности стержня.
После процесса волочения калиброванной TRIP стали результатом является заготовка для стержневых деталей с повышенными точностью размеров, качеством поверхности, износостойкостью и запасом прочности. То есть получается конечный продукт с адаптирующимися под негативные воздействия рабочей среды механическими свойствами или передельная заготовка с высокой технологической гибкостью.
Особенности процесса волочения калиброванной TRIP стали потенциально определяют возможность получения продукта с широким спектром механических свойств:
1. При малых скоростях деформации результатом волочения будет холоднотянутая калиброванная TRIP сталь с высокими пластическими характеристиками и повышенным ресурсом потенциального упрочнения за счет ^трансформированного остаточного аусте-нита в микроструктуре. Такая заготовка представляет интерес для производителей стержневых конструкций, защитных рам или стержневых крепежных изделий. Большой ресурс нетрансформированного остаточного аустенита позволяет в данных видах производства управлять локализацией механических свойств в необходимых элементах готовой продукции (например, локально упрочнять головку или резьбу болта).
2. При высоких скоростях деформации при волочении становится возможным получение стержней с высокими классами прочности, превышающими традиционные углеродистые стали.
3. При варьировании угла волоки и степени деформации в результате волочения можно получить калиброванную сталь с упроченным поверхностным слоем заданной толщины.
Следовательно, для проектирования режимов деформационной обработки и прогнозирования конечных свойств заготовки необходима разработка методов моделирования, которые обеспечивают учет динамических структурно-фазовых превращений в металле, но в то же время позволяют осуществлять оперативный анализ свойствообразующих
параметров НДС при относительно высоких объемах деформируемого металла.
Моделирование
Макромоделирование
Для оценки технологического потенциала процесса было выбрано однократное волочение калиброванной стали с исходного диаметра 38 мм. Данный диаметр наиболее востребован во всех областях машиностроения.
Для изучения влияния обжатия на свойства металла финальный диаметр варьировался от 37 до 35 мм. Угол раскрытия волоки варьировался в пределах от 18 до 12°. Скорость волочения составляла 100 мм/с.
В качестве исходного материала была выбрана сталь TRIP 700, химический состав которой представлен в таблице.
Механические свойства стали TRIP 700 представлены на рис. 1.
Микромоделирование (Representative
Volume Element - RVE)
Для субмоделирования была использована фотография микроструктуры стали TRIP 700 (рис. 2). Исходная микроструктура исследуемой стали имела следующий состав: остаточный аустенит (~ 22 %), феррит (~ 66 %) и бейнит (~ 12 %).
Бинаризированное изображение микроструктуры было покрыто КЭ сеткой с треугольными элементами типа CAX3 (A 3-node linear axisymmetric triangle) из стандартной библиотеки Abaqus.
Химический состав стали марки TRIP 700, % (масс.)
C Mn Si Cr Mo Ni Cu Al P S
0,177 1,54 0,32 0,02 - 0,03 0,02 1,35 0,075 0,007
Рис. 1. Механические свойства стали TRIP 700
Рис. 2. Микроструктура стали TRIP 700
Для симуляции трип-эффекта была выбрана следующая экспериментально установленная математическая зависимость [12, 13]:
Vy = Vy-,n,t,al - 26,708s - 0,003S +
+ 53,516s2 - 0,018ss -1,06-10-5ё2, (1) где vy-initial - исходное количество остаточного аустенита в микроструктуре TRIP стали; ё - скорость деформации; s - степень деформации.
С учетом зависимости (1) была разработана надстройка для пакета Abaqus, реализующая следующий алгоритм:
1. Расчет объемной доли трансформированного аустенита после каждого шага расчета.
2. Выбор рассчитанного процента из группы элементов остаточного аустенита, которые имеют максимальные значения деформации и скорости деформации.
3. Замена свойств выбранных элементов на свойства мартенсита.
Результаты
Сравнение результатов моделирования
с симуляцией TRIP эффекта и без него
На рис. 3 представлено распределение напряжений Мизеса и эквивалентной деформации в центральных слоях металла в деформационной зоне.
Трансформация относительно пластичного остаточного аустенита в более прочный мартенсит приводит к появлению локализованных высоких значений эквивалентных напряжений. При этом, поскольку только что образовавшийся мартенсит является наиболее прочной структурной составляющей во всей микроструктуре TRIP стали, его скопления перестают деформироваться. Поэтому количество локализаций околонулевых значений эквивалентных деформаций в модели с симу-
ляцией TRIP-эффекта значительно больше.
Особенно ярко эта закономерность прослеживается на поверхности деформируемой калиброванной стали (рис. 4), где скорость и степень деформации выше, чем в центральных слоях, а значит выше и количество преобразованного остаточного аустенита.
Анализ деформированного состояния и микроструктурных преобразований в ходе процесса волочения Значения радиальных деформаций сжатия и глубина их проникновения в металл в деформационной зоне главным образом определяют распределение механических свойств в конечном изделии. Особенно это важно для исследуемых TRIP сталей, в которых из-за протекания TRIP-эффекта динамика изменения механических свойств стали несоизмеримо выше, чем в сталях без метастабильных микроструктурных составляющих.
Анализ радиальных деформаций (рис. 5) осуществлялся в центральном и поверхностном слое металла в середине деформационной
зоны волоки.
Локализации околонулевых деформаций обусловлены высокой прочностью бейнита и образовавшихся зародышей мартенсита. Более пластичные фазы (феррит и остаточный аустенит) испытывают сжимающие радиальные деформации в диапазонах значений от -0,03 до -0,15 (в центре) и от -0,07 до -0,22 (на поверхности). При этом отмечено, что локализации наиболее высоких значений сжимающих радиальных деформаций пластичных фаз находятся в областях высокого скопления более прочных микроструктурных составляющих, то есть часть микродеформаций в металле обусловлена деформационным взаимодействием элементов микроструктуры между собой.
Рис. 3. Распределение эквивалентных напряжений и деформаций в центральном слое металла в модели с симуляцией TRIP эффекта и без него
Рис. 4. Распределение эквивалентных напряжений и деформаций в поверхностном слое металла в модели с симуляцией TRIP эффекта и без него
Рис. 5. Распределение радиальных деформаций на поверхности и в центре деформационной зоны
Рис. 6. Влияние угла волоки на количество мартенсита в микроструктуре
калиброванной стали
Рис. 7. Влияние обжатия на количество мартенсита в микроструктуре калиброванной стали (исходный диаметр 38 мм)
Данный факт должен приниматься во внимание при проектировании технологий производства калиброванной стали, включающих пластическую деформацию и последующую термическую обработку, как правило, в виде операций рекристаллизационного
отжига. В этом случае необходимо избегать так называемых «критических» степеней деформации при назначении режимов волочения калиброванной стали. После так называемой критической степени деформации не происходит процесс рекристаллизации по ме-
ханизму образования новых зерен и их роста. Нагрев в таком случае вызывает быстрый рост одних исходных нерекристаллизованных зерен за счет поглощения соседних. Поэтому первичная рекристаллизация осуществляется в ограниченных объемах, а рост зерен при вторичной рекристаллизации практически отсутствует. При критических степенях деформации возникает заметное проявление неодинаково деформированных соседних зерен и отсутствуют условия для появления центров первичной рекристаллизации в отдельных объемах металла. Все это приводит к укрупнению одних слабодеформированных зерен за счет других, а впоследствии - к ухудшению механических свойств стали. Также неоднородность свойств отрицательно сказывается на стабильности процесса дальнейшей переработки и качестве готовой продукции.
Было исследовано влияние геометрических параметров инструмента на протекание мартенситного превращения в центральных и поверхностных слоях калиброванной стали (рис. 6, 7).
Отмечено, что при увеличении угла волочильного инструмента количество образовавшегося мартенсита также увеличивается вследствие роста скорости деформации. Однако при больших углах волоки неравномерность деформации растет, что приводит к увеличению разницы процентной доли мартенсита в микроструктуре между центральными и периферийными слоями. Уменьшение угла волоки позволяет снизить эту разницу, но при почти на треть меньшем количестве образовавшегося мартенсита.
Степень деформации оказывает относительно меньшее влияние на образование мартенсита и распределение его процентного соотношения по сечению калиброванной стали.
Заключение
Результаты расчета НДС в мультимас-штабных моделях продемонстрировали не только перспективность данного инженерного инструмента, но и неотъемлемую необходимость учета микроструктуры и ее поведения под нагрузкой при проектировании даже традиционных процессов с использованием современных материалов.
Предложенный метод создания мульти-масштабных моделей в случае исследования процесса волочения калиброванной TRIP стали позволил:
1. Доказать важность учета TRIP эффекта с позиции НДС, так как формирующиеся в ходе его реализации зародыши высокопрочного мартенсита являются локализаторами высоких значений напряжений Мизеса.
2. Обнаружить высокие сжимающие деформации в диапазонах значений от -0,03 до -0,15 (в центре) и от -0,07 до -0,22 (на поверхности) в фазах феррита и остаточного аустенита.
3. Изучить деформационное взаимодействие элементов микроструктуры между собой, тем самым объяснив наличие более высоких радиальных деформаций в пластичных фазах вблизи большего скопления зерен более прочных фаз бейнита и мартенсита.
В совокупности это позволило не просто повысить прогнозирующую способность моделей процесса и скорость их расчета, но и расширить номенклатуру материалов традиционного процесса волочения калиброванной стали путем использования современных, технологически гибких TRIP сталей.
Работа проведена при финансовой поддержке Минобрнауки России в рамках реализации комплексного проекта по созданию высокотехнологичного производства с участием высшего образовательного учреждения (Договоры № 02.G25.31.0178 от 01.12.2015 г.; № МК204895 от 27.07.2015 г.).
^HTepaTypa/References
1. Polyakova M.A., Gulin A.V., Constantinov D.V. Investigation of Microstructure and Mechanical Properties of Carbon Steel Wire after Continuous Method of Deformational Nanostructuring. Applied Mechanics and Materials, 2013, vol. 436, pp. 114-120.
2. Korchunov A.G., Polyakova M.A Gulin A.V. Konstantinov D.V. Techno-logical Inherited Connections in Continuous Method of Deformational Nanostructuring. Applied Mechanics and Materials, 2014, vol. 555, pp. 401-405. DOI: 10.4028/www.scientific.net/AMM.555.401
3. Sarma G.B., Radhakrishnan B., Zacharia T. Finite Element Simulations of Cold Deformation at the Mesoscale, Computational Materials Science, 1998, vol. 12, pp. 105-123.
4. Son S.B., Roh H., Kang S.H. Relationship between Microstructure Homogeneity and Bonding Stability of Ultrafine Gold Wire. Journal of Materials Science, 2010, vol.45, pp. 236-244.
5. Watanabe I., Setoyama D., Nagasako N., Iwata N., Nakanishi K. Multiscale Prediction of Mechanical Behavior of Ferrite-Pearlite Steel with Numerical Material Testing. International Journal for Numerical Methods in Engineering, 2012, vol. 89, pp. 829-845.
6. Fischer F.D., Reisner G., Werner E., Tanaka K., Cailletaud G., Antretter T. A New View on Transformation Induced Plasticity (TRIP). International Journal of Plasticity, 2000, vol. 16 (7-8), pp. 723-748.
7. Bleck W., Schael I. Determination of Crash-Relevant Material Parameters by Dynamic Tensile Tests. Steel Research, 2000, vol. 71, pp. 173-178.
8. Choi I.D., Bruce D.M., Kim S.J. Deformation Behavior of Low Carbon TRIP Sheet Steels at High Strain Rates. ISIJInternational, 2002, vol. 42, pp. 1483-1489.
9. Kim YG, Kim TW, Han JK, Chang RW. Development of Newaustenitic Fe-Mn-Al-C Steels for Automotive Applications. Key Engineering Materials, 1993, vol. 84, pp. 461-72.
10. Grässel O, Krüger L, Frommeyer G, Meyer LW. High Strength Fe-Mn-(Al,Si) TRIP/TWIP Steels Development - Properties - Application. International Journal of Plasticity, 2000, vol. 16, pp. 391-409.
11. Frommeyer G, Brüx U, Neumann P. Supra-Ductile and High-Strength Manganese-TRIP/TWIP Steels for High Energy Absorption Purposes. ISIJ International, 2003, vol. 43, pp. 43-46.
12. Wiewiorowska S., Determination of Content of Retained Austenite in Steels with TRIP Effect Deformed at Different Strain Rates. Steel Research International, 2010, vol. 81, pp. 262-265.
13. Wiewiorowska S. The Influence of Strain Rate and Strain Intensity on Retained Austenite Content in Structure of Steel with TRIP Effect. Solid State Phenomena, 2010, vol. 165, pp. 216-222.
Константинов Дмитрий Вячеславович, канд. техн. наук, ведущий специалист по вопросам международных научных исследовательских проектов, Магнитогорский государственный технический университет им. Г.И. Носова, г. Магнитогорск; const_dimon@mail.ru.
Корчунов Алексей Георгиевич, д-р техн. наук, профессор, проректор по международной деятельности, Магнитогорский государственный технический университет им. Г.И. Носова, г. Магнитогорск; international@magtu.ru.
Ширяев Олег Петрович, директор, ОАО «Магнитогорский метизно-калибровочный завод «ММК-МЕТИЗ», г. Магнитогорск.
Кузнецова Алла Сергеевна, научный сотрудник научно-исследовательского сектора, Магнитогорский государственный технический университет им. Г.И. Носова, г. Магнитогорск; allakuznetsova.mgtu@mail.ru.
Емалеева Динара Гумаровна, канд. техн. наук, доцент кафедры технологий обработки материалов, Магнитогорский государственный технический университет им. Г.И. Носова, г. Магнитогорск; emaleevadg@mail.ru.
Поступила в редакцию 19 мая 2017 г.
DOI: 10.14529/met170410
RESEARCH OF MACRO- AND MICROMECHANICS OF STEEL DEFORMATION DURING DRAWING BY MULTISCALE COMPUTER SIMULATION
D.V. Konstantinov1, const_dimon@mail.ru, A.G. Korchunov1, international@magtu.ru, O.P. Shiryaev2,
A.S. Kuznetsova1, allakuznetsova.mgtu@mail.ru, D.G. Emaleeva1, emaleevadg@mail.ru
1 Nosov Magnitogorsk State Technical University, Magnitogorsk, Russian Federation,
2 OJSC Magnitogorsk Hardware and Sizing Plant "MMK-METIZ", Magnitogorsk, Russian Federation
The article compares the results of modelling allowing for the TRIP effect and overlooking it. The comparison shows that martensite, which forms in the microstructure, causes a significant rise of equivalent stresses across the entire representative volume, which is extremely critical in the drawing process engineering. The method applied also gave a better understanding of how the microstructural elements interact in TRIP steel under strain, which helped explain more intense (2 or 3 times as high compared with the average values) radial strains in the plastic phases close to bigger grain clusters of stronger bainite and martensite phases.
Keywords: drawing, steel, TRIP effect, multiscale simulation, representative volume, stressstrain state, microstructure.
Received 19 May 2017
ОБРАЗЕЦ ЦИТИРОВАНИЯ
Исследование макро- и микромеханики деформирования стали в ходе процесса волочения на основе мультимасштабного моделирования / Д.В. Константинов, А.Г. Корчунов, О.П. Ширяев и др. // Вестник ЮУрГУ. Серия «Металлургия». - 2017. - Т. 17, № 4. - С. 92-100. DOI: 10.14529/теП 70410
FOR CITATION
Konstantinov D.V., Korchunov A.G., Shiryaev O.P., Kuznetsova A.S., Emaleeva D.G. Research of Macro- and Micromechanics of Steel Deformation during Drawing by Multiscale Computer Simulation. Bulletin of the South Ural State University. Ser. Metallurgy, 2017, vol. 17, no. 4, pp. 92-100. (in Russ.) DOI: 10.14529/met170410
