Математическая модель деформирования стали при мартенситных переходах

Аннотация


Задача описания деформации при фазовых переходах является составной частью проблемы построения двухуровневой модели неупругого деформирования стали при термомеханическом воздействии, основанной на подходах физической теории пластичности. Градиент трансформационной деформации представительного мезообъема построен на основе кристаллографической теории мартенситного перехода и описывает деформирование с инвариантной плоскостью. Он выражается через вектор нормали к инвариантной плоскости и вектор сдвига, которые задают трансформационную систему (по аналогии с системой кристаллографического сдвига для дислокационного скольжения). Согласно ориентационным соотношениям Курдюмова - Закса, существуют 24 трансформационные системы, выражения для компонент которых приводятся в работе. Рассматриваются особенности применения меры неупругого деформирования представительного объема мезоуровня (кристаллита) при мартенситном превращении, в частности определяются условия активации этих систем при разных типах нагружения.

Полный текст

Современные технологии металлургии широко применяют термомеханическую обработку для получения требуемых механических свойств металлических заготовок и конструкций. Процесс обработки включает сложные взаимосвязанные физические явления, такие как теплоперенос, диффузия компонентов сплава, неупругое деформирование, полиморфные превращения в стали. Описание этих процессов вовлекает в рассмотрение различные масштабные уровни, для которых характерны те или иные физические механизмы. Настоящая работа является составной частью задачи построения двухуровневой модели неупругого деформирования стали под действием приложенных механических нагрузок и изменения температуры. Модель описывает поведение стали при деформировании на уровне представительного объема макроуровня и на мезоуровне (уровне отдельного зерна, кристаллита); математическая постановка для мезоуровня включает в себя конститутивные соотношения физической теории пластичности (закон Гука в скоростной релаксационной постановке), критерии полиморфного превращения, эволюционные уравнения для доли сосуществующих фаз и кинематические соотношения для меры деформации [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 трансформационных систем. Исследуется последовательность активации трансформационных систем под действием приложенного напряжения: одно-, двух-, трехосного растяжения и простого сдвига. Установлено, что в первую очередь активизируются системы, благоприятно ориентированные по отношению к направлению приложенных нагрузок, и еще несколько систем из других групп.

Об авторах

Н. Д. Няшина

Пермский национальный исследовательский политехнический университет

Список литературы

  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).

Статистика

Просмотры

Аннотация - 32

PDF (Russian) - 15

Ссылки

  • Ссылки не определены.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах