Mathematical model of the steel deformation at martensitic transformation

Abstract


The task of the deformation description at phase transformation is a part of a problem involving creation of two-level model of steel inelastic deformation at thermomechanical treatment based on approaches of the physical plasticity theory. The gradient of transformational deformation of representative mesovolume is constructed on the basis of the crystallographic theory of martensitic transformation and describes deformation with the invariant plane. It is expressed through a vector of a normal to the invariant plane and a shear vector which set transformational system (by analogy with system of crystallographic shear for dislocation sliding). According to the Kurdyumov - Sachs’s orientation relationships, there are 24 transformational systems, expressions for which component are given in the article. Features associated with application of a measure of inelastic deformation of mesolevel (crystallite) representative volume at martensitic transformation are considered. In particular conditions of activation of these systems at different types of loading are defined.

Full Text

Современные технологии металлургии широко применяют термомеханическую обработку для получения требуемых механических свойств металлических заготовок и конструкций. Процесс обработки включает сложные взаимосвязанные физические явления, такие как теплоперенос, диффузия компонентов сплава, неупругое деформирование, полиморфные превращения в стали. Описание этих процессов вовлекает в рассмотрение различные масштабные уровни, для которых характерны те или иные физические механизмы. Настоящая работа является составной частью задачи построения двухуровневой модели неупругого деформирования стали под действием приложенных механических нагрузок и изменения температуры. Модель описывает поведение стали при деформировании на уровне представительного объема макроуровня и на мезоуровне (уровне отдельного зерна, кристаллита); математическая постановка для мезоуровня включает в себя конститутивные соотношения физической теории пластичности (закон Гука в скоростной релаксационной постановке), критерии полиморфного превращения, эволюционные уравнения для доли сосуществующих фаз и кинематические соотношения для меры деформации [1-4]. В работе [5] приводится подробный вывод и обоснование применения в многоуровневых моделях градиента деформации, связанной с мартенситными переходами, для представительного объема мезоуровня. Выкладки базируются на использовании кристаллографической теории мартенситных превращений в сталях [6-9], в предположении, что трансформационная деформация является деформацией с инвариантной плоскостью. Градиент трансформационной деформации для представительного мезообъема родительской фазы аустенита, в котором появляется пластинка мартенсита, для одной трансформационной системы имеет вид [7-9] (1) где m и s - вектор нормали к инвариантной плоскости габитуса и вектор сдвига соответственно (задают трансформационную систему); векторы в общем случае не перпендикулярны. Эта мера описывает деформацию с инвариантной плоскостью, поскольку все плоскости с нормалью m не искажаются и не поворачиваются, так как смещения в каждой из них одинаковы и пропорциональны расстоянию от габитусной плоскости. Соотношение (1) для градиента деформации записано для одной трансформационной системе. Однако согласно ориентационному соотношению Курдюмова - Закса, в зерне аустенита существуют 24 эквивалентные трансформационные системы. Тогда градиент трансформационной деформации можно записать в виде (2) где - доля мартенситной фазы, образовавшейся по b-й трансформационной системе. Все 24 трансформационные системы могут быть получены из одной, выведенной в рамках кристаллографического описания мартенситного перехода, путем всевозможных перестановок компонент и изменения знаков. В табл. 1 приведены значения для компонент векторов всех трансформационных систем. Следует отметить, что они группируются по 8 систем в множества, связанные с ориентацией тензора Бейна. Таблица 1 Компоненты векторов 24 трансформационных систем № п/п m s 1 -0,837 -0,533 0,122 -0,286 0,426 0,112 2 -0,837 0,533 0,122 -0,286 -0,426 0,112 3 -0,837 0,122 -0,533 -0,286 0,112 0,426 4 -0,837 0,122 0,533 -0,286 0,112 -0,426 5 -0,837 -0,122 0,533 -0,286 -0,112 -0,426 6 -0,837 -0,122 -0,533 -0,286 -0,112 0,426 7 -0,837 0,533 -0,122 -0,286 -0,426 -0,112 8 -0,837 -0,533 -0,122 -0,286 0,426 -0,112 9 -0,533 -0,837 0,122 0,426 -0,286 0,112 10 0,533 -0,837 0,122 -0,426 -0,286 0,112 11 0,122 -0,837 -0,533 0,112 -0,286 0,426 12 0,122 -0,837 0,533 0,112 -0,286 -0,426 13 -0,122 -0,837 0,533 -0,112 -0,286 -0,426 14 -0,122 -0,837 -0,533 -0,112 -0,286 0,426 15 0,533 -0,837 -0,122 -0,426 -0,286 -0,112 16 -0,533 -0,837 -0,122 0,426 -0,286 -0,112 Окончание табл. 1 № п/п m s 17 0,122 -0,533 -0,837 0,112 0,426 -0,286 18 0,122 0,533 -0,837 0,112 -0,426 -0,286 19 0,533 0,122 -0,837 0,426 0,112 -0,286 20 -0,533 0,122 -0,837 -0,426 0,112 -0,286 21 0,533 -0,122 -0,837 -0,426 -0,112 -0,286 22 -0,533 -0,122 -0,837 0,426 -0,112 -0,286 23 -0,122 0,533 -0,837 -0,112 -0,426 -0,286 24 -0,122 -0,533 -0,837 -0,112 0,426 -0,286 Тензор деформации Бейна описывает преобразование элементарной ячейки ГЦК аустенита в ячейку ОЦТ мартенсита в предположении минимальности перемещений атомов. Как видно из рис. 1, в двух смежных ячейках ГЦК-решетки можно формально выделить ячейку ОЦТ (рис. 1, а), параметры которой отличаются от наблюдаемых в мартенсите (рис. 1, б); таким образом, для получения реальной ячейки ОЦТ мартенсита формально выделенную необходимо сжать вдоль одной оси и растянуть вдоль другой (рис. 1, в) [10]. Поскольку ГЦК-решетка симметрична относительно осей [100], [010], [001] и подобную ОЦТ-ячейку можно выделить в любом из этих направлений, деформация Бейна, соответственно, может быть выражена одним из трех тензоров В1, В2 или В3, компоненты которых в базисе кристаллографической системы координат решетки аустенита имеют вид (3) Здесь а0 и а - периоды исходной решетки аустенита и мартенсита соответственно. a б в Рис. 1. Геометрическая интерпретация деформации Бейна: а - формально выделенная ячейка ОЦТ-решетки (жирная линия) в ГЦК-решетке аустенита; б - реальная ячейка ОЦТ-решетки мартенсита; в - деформация Бейна, превращающая формально выделенную в аустените ячейку ОЦТ-решетки (жирная линия) в реальную ячейку ОЦТ-решетки мартенсита [10] На рис. 2 изображены стереографические проекции указанных трансформационных систем, полученные для каждого из трех тензоров Бейна. Для сравнения приведены проекции нескольких характерных для кубических решеток плоскостей (типа [100], [101], [110] и т.д.). Заметно, что они группируются определенным образом: вблизи плоскостей [101]-[-101] для тензора деформации Бейна В1, [011]-[0-11] - для В2, [-101]-[0-11] - для В3. а б в Рис. 2. Стереографические проекции инвариантных плоскостей для 24 трансформационных систем (см. табл. 1): а - плоскости, связанные с деформацией Бейна В1; б - В2; в - В3 (черным цветом отмечены некоторые кристаллографические плоскости, характерные для ОЦТ кристалла) (авторские результаты) Для анализа условий активации трансформационных систем под действием приложенных нагрузок воспользуемся простейшим термодинамическим критерием. Как известно из экспериментов [6, 11], значение температуры, при которой начинается мартенситный переход, зависит от приложенного напряжения: чем выше уровень приложенных нагрузок, тем выше температура, при которой начинается переход. Такую зависимость можно объяснить тем, что при высокой температуре нет достаточной движущей силы, чтобы преодолеть потенциальный барьер, необходимый для фазового перехода, т.е. разность между свободной энергией мартенситной фазы и аустенитной фазы недостаточна для активизации фазового перехода. Таким образом, для начала мартенситного перехода при температурах выше равновесной системе необходимо передать дополнительную энергию, например, совершив работу над системой [8]: Здесь U - работа приложенного напряжения на деформациях, вызванных мартенситным переходом. При достижении U критического значения при заданной температуре происходит фазовое превращение. Проводились исследования поведения указанных трансформационных систем под действием простого нагружения (одно-, двух-, трехосное растяжение, простой сдвиг) представительного объема мезоуровня. Результаты расчетов приведены на рис. 3. ?i N а N ?i б ?i N в Рис. 3. Число активных трансформационных систем в зависимости от величины нагрузки при разных схемах нагружения (авторские результаты) N Как показали расчеты, при одно- и двухосном растяжении даже при небольшом уровне нагрузки сразу активизируются 16 трансформационных систем - это в первую очередь 8 благоприятно ориентированных систем, полученных с использованием тензора Бейна, «сонаправленного» с направлением нагрузки. В табл. 2 отмечены номера трансформационных систем, которые активизируются в первую очередь (подчеркиванием отмечены номера систем, полученные из одного тензора деформации Бейна). Таблица 2 Номера активных трансформационных систем при различных видах нагружения Вид нагружения Число активных систем Номера активных систем (см. табл. 1) Одноосное растяжение s11 ? 0 16 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15,16, 19, 20, 21, 22 Одноосное растяжение s22 ? 0 16 1, 2, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 23, 24 Простой сдвиг s12 = s21 ? 0 8 1, 2, 7, 8, 9, 10, 15, 16 16 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,13, 14, 15, 16 Простой сдвиг s13 = s31 ? 0 8 3, 4, 5, 6, 19, 20, 21, 22 16 1, 2, 3, 4, 5, 6, 7, 8, 17, 18, 19, 20, 21, 22, 23, 24 Двухосное растяжение s11 = s22 ? 0 16 3, 4, 5, 6, 11, 12, 13, 14, 17, 18, 19, 20, 21, 22, 23, 24 Для простого сдвига вначале активизируются 8 систем, а затем практически сразу еще 8. Эти 16 систем получены для двух тензоров Бейна, «сонаправленых» с приложенными нагрузками. Для активизации всех 24 систем приходится прикладывать нагрузку, в разы большую, чем для первично активированных систем. Для трехосного растяжения активируются сразу 24 системы, но это происходит при уровнях нагрузки, значительно меньших, чем для одно-, двухосного растяжения и простого сдвига, но значительно больших, чем необходимо для первичной активации 16 трансформационных систем в рассмотренных примерах. Выводы В работе приведены выражения для компонент векторов 24 трансформационных систем. Исследуется последовательность активации трансформационных систем под действием приложенного напряжения: одно-, двух-, трехосного растяжения и простого сдвига. Установлено, что в первую очередь активизируются системы, благоприятно ориентированные по отношению к направлению приложенных нагрузок, и еще несколько систем из других групп.

About the authors

N. D. Niashina

Perm National Research Polytechnic University

References

  1. Исупова И.Л., Трусов П.В. Математическое моделирование фазовых превращений в сталях при термомеханической нагрузке // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2013. - № 3. - С. 126-156.
  2. Трусов П.В., Ашихмин В.Н., Швейкин А.И. Двухуровневая модель упругопластического деформирования поликристаллических материалов // Механика композиционных материалов и конструкций. - 2009. - Т. 15, № 3. - С. 327-344.
  3. Трусов П.В., Швейкин А.И. Многоуровневые физические модели моно- и поликристаллов. Статистические модели // Физическая мезомеханика. - 2011. - Т. 14, № 4. - С. 17-28.
  4. Трусов П.В., Швейкин А.И. Многоуровневые физические модели моно- и поликристаллов. Прямые модели // Физическая мезомеханика. - 2011. - Т. 14, № 5. - С. 5-30.
  5. Няшина Н.Д., Трусов П.В. Моделирование мартенситных превращений в сталях: кинематика мезоуровня // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 4. - С. 118-151. doi: 10.15593/perm.mech/2014.4.05.
  6. Курдюмов Г.В., Утевский Л.М., Энтин Р.И. Превращения в железе и стали. - М.: Наука,1977. - 240 c.
  7. Bhadeshia H.K.D.H., Kundu S., Abreu H. Mathematics of crystallographic texturein martensitic and related transformations // Microstructure and Textures in steel and other materials. - 2009. - Vol. XIX. - P. 19-31.
  8. Perdahcioğlu E.S. Constitutive Modeling of Metastable Austenitic Stainless Steel. Ph.D. Thesis. - Enschede, The Netherlands, 2008. - 145 р.
  9. Wechsler M.S., Liebermann D.S., Read Т.A. On the theory of the formation of martensite // J. Metals. - 1953. - Vol. 5, sect. 2, no. 11. - Р. 1503-1515.
  10. Физика твердого тела / Г.Н. Елманов, А.Г. Залужный, В.И. Скрытный, Е.А. Смирнов, В.Н. Яльцев // Физическое материаловедение: учебник для вузов: в 6 т.; под ред. Б.А. Калина. - М.: Изд-во Моск. физ.-инжен. ин-та, 2007. - Т. 1. - 236 с.
  11. Центральный металлический портал РФ. Мартенсит - образование и превращения. Ч. 1. [Электронный ресурс]. - URL: http://metallicheckiy-portal.ru/articles/chermet/fazovie_sostoyania/martensit _obrazovanie_i_prevrashenia/1 (дата обращения: 7.09.2014).

Statistics

Views

Abstract - 15

PDF (Russian) - 8

Refbacks

  • There are currently no refbacks.

This website uses cookies

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

About Cookies