Численное моделирование процесса управления фазовым переходом при кручении полого цилиндра из сплава Гейслера
- Авторы: Роговой АА1, Столбова ОС1
- Учреждения:
- Институт механики сплошных сред УрО РАН
- Выпуск: № 3 (2019)
- Страницы: 75-87
- Раздел: Статьи
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/526
- DOI: https://doi.org/10.15593/perm.mech/2019.3.08
- Цитировать
Аннотация
Рассматриваются ферромагнитные сплавы с эффектом памяти формы (сплавы Гейслера), испытывающие фазовый переход из высокотемпературной кубической фазы (аустенит) в низкотемпературную тетрагональную фазу (мартенсит) в ферромагнитном состоянии. В этих сплавах в процессе прямого температурного фазового перехода из аустенитного состояния в мартенситное возможны генерации значительных макроскопических деформаций за счет приложения механических напряжений. Поскольку критические температуры процесса зависят от магнитного поля и полей напряжений, в таких сплавах возможно управление процессом фазового перехода аустенит-мартенсит с помощью магнитного поля. Представлена модель, позволяющая в рамках конечных деформаций описывать процесс управления прямым (аустенит-мартенсит) и обратным (мартенсит-аустенит) фазовыми переходами с помощью магнитного поля в ферромагнитных поликристаллических материалах с памятью формы при действии внешних силовых, тепловых и магнитных полей. Поскольку магнитное поле влияет на деформацию материала, которая, в свою очередь, изменяет магнитное поле, поставлена связанная краевая задача и рассмотрен пример об азимутальном кручении длинного полого цилиндра (плоская деформация) из сплава Гейслера. Задача реализована численно методом конечных элементов с использованием процедуры пошагового нагружения. Накопление фазовых деформаций происходило в процессе прямого фазового перехода в цилиндре, внешняя поверхность которого была предварительно в аустенитном состоянии закручена в азимутальном направлении относительно фиксированной внутренней. При этом величина накопленных фазовых деформаций зависела от того, какие граничные условия (силовые или кинематические), производящие закручивание, заданы на внешней поверхности цилиндра. Полное снятие накопленных деформаций и раскручивание внешней поверхности цилиндра обратно происходило как при нагревании образца в интервале температур обратного фазового перехода в отсутствие магнитного поля, так и при постоянной температуре, превышающей этот интервал, при снятии магнитного поля, приложенного предварительно в мартенситном состоянии.
Полный текст
Введение Ферромагнитные сплавы Гейслера, в которых наблюдается эффект памяти формы, относятся к классу функциональных (интеллектуальных) материалов [1-14]. Конструкции из таких материалов могут существенно изменять свою конфигурацию под действием внешних тепловых, магнитных или электрических полей. Среди сплавов Гейслера с эффектом памяти формы наиболее изученной является система Ni-Mn-Ga, которая демонстрирует структурный (мартенситный) фазовый переход из высокотемпературной кубической фазы (аустенит) в низкотемпературную тетрагональную фазу (мартенсит) в ферромагнитном состоянии. Выделяют два различных физических механизма, ответственных за инициирование больших деформаций (до 10 %) в ферромагнитных сплавах с памятью формы под действием внешних сил и магнитных полей. Один из них связан с восстановлением структурных доменов в мартенситной фазе и наблюдается в монокристаллических материалах или сильно текстурированных поликристаллических образцах [1-9]. При отсутствии приложенной механической нагрузки прямое фазовое превращение из кубической высокотемпературной фазы (аустенит) в тетрагональную низкотемпературную фазу (мартенсит) приводит к образованию сдвойникованных вариантов мартенситного состояния (структурные домены). При этом изменение симметрии кристаллической решетки вызывает деформации, которые при пространственном усреднении по репрезентативному объему материала, состоящему из определенного числа структурных доменов, приводят только к несущественной макроскопической объемной деформации. Структурный домен, в свою очередь, делится на магнитные домены, в которых векторы намагниченности имеют противоположные направления: в каждом магнитном домене этот вектор направлен вдоль или против оси легкого намагничивания данного домена, при этом домены согласованы таким образом, что минимизируют магнитостатическую энергию. Приложение внешнего магнитного поля к ферромагнитным сплавам вызывает движение стенок магнитных доменов, вращение вектора намагниченности и, при условии высокой магнитной анизотропии, переориентацию мартенсита. Первые два процесса также происходят в обычных ферромагнитных материалах, в то время как последний процесс присущ только сплавам с памятью формы: приложение внешнего магнитного поля к материалу в мартенситном состоянии (так же, как приложение силового поля) приводит к вращению структурных доменов (движению границ двойников) таким образом, что оси легкого намагничивания пытаются выстроиться по индуцированному в теле магнитному полю. Совместная реориентация определенного числа мартенситных вариантов приводит к раздвойникованию мартенсита и сопровождается макроскопической деформацией, которая в некоторых ферромагнитных сплавах может быть чрезвычайно сильной. В случае снятия магнитного поля деформация, создаваемая вращением мартенситных вариантов, не восстанавливается, но частично или полностью исчезает, когда материал переводится в аустенитное состояние в интервале температур обратного фазового перехода. Другой механизм генерации значительных макроскопических деформаций в материалах связан с прямым температурным фазовым переходом из аустенитного состояния в мартенситное в присутствии механических напряжений. Этот механизм имеет место как в монокристаллах, так и в поликристаллах [1, 2, 10-14]. В соответствии с ним формирование типа мартенситного варианта полностью зависит от направления поля напряжений, так что возникающие деформации ориентированы в соответствии с этим напряжением. Эти деформации называются фазовыми и дополняют обычные деформации, возникающие в материалах из-за приложенных сил. Усреднение фазовых деформаций по репрезентативному объему материала приводит к макроскопической деформации, которая может быть значительно выше, чем обычная деформация в мартенситном состоянии, и, в отличие от обычной деформации, не исчезает после силовой разгрузки материала. С помощью магнитного поля можно эффективно управлять процессом фазового перехода, о чем свидетельствуют экспериментальные исследования [10-12]. Если прямые или обратные фазовые переходы в ферромагнитном материале осуществляются под действием магнитного поля, критические температуры процесса изменяются в соответствии с обобщенным законом Клаузиуса-Клапейрона и зависят от магнитного поля и полей напряжений [13, 14]. Такое влияние магнитного поля и поля напряжений позволяет варьировать температуры прямых и обратных фазовых превращений и тем самым управлять процессом фазового перехода аустенит-мартенсит. В серии работ [15-17] задача магнитоупругости исследуется в рамках теории конечных деформаций. Внутреннее магнитное поле, индуцированное внешним магнитным полем в образце конечных размеров, генерирует как массовые (пондеромоторные), так и поверхностные силы в дополнение к обычным силовым или температурным внешним воздействиям. В свою очередь индуцированное внутреннее магнитное поле и деформация образца оказывают сильное влияние на внешнее магнитное поле, вызывая его возмущение. В результате внутреннее поле изменяется, что приводит к изменению массовых и поверхностных сил, генерируемых этим полем, и к изменению формы образца. Кроме того, вектор магнитной индукции и вектор напряжённости магнитного поля имеют на поверхности тела конечной геометрии разрывы в их тангенциальной и нормальной составляющих соответственно. Это порождает поверхностный эффект, который проявляется в определенной области и не учитывается в большинстве работ, посвященных рассматриваемой проблеме. Таким образом, необходимо решить связанную магнитомеханическую задачу для намагничиваемого тела конечных размеров, которое находится в пространстве, где действует магнитное поле. Для построения физико-механической модели поведения ферромагнитного сплава с памятью формы (ФСПФ) в процессах прямого и обратного фазовых превращений, которая будет учитывать влияние теплового, магнитного и силового внешних полей на разных стадиях этих процессов, будем использовать три группы соотношений: 1) соотношения механики деформируемого твердого тела; 2) соотношения, описывающие прямой и обратный температурно-индуцированные управляемые фазовые переходы аустенит-мартенсит; 3) соотношения, описывающие изменение магнитного поля в кусочно-однородной среде (уравнения Максвелла). Эти соотношения с краткими пояснениями и необходимыми ссылками на первоисточники приведены ниже в разд. 1-3 для удобства читателя. В разд. 4 разработанная физико-механическая модель поведения ферромагнитного сплава с памятью формы используется для описании магнитомеханических процессов, происходящих в длинном полом цилиндре из ферромагнитного материала с закрепленной внутренней поверхностью при прямом и обратном фазовых переходах. В аустенитном состоянии внешняя поверхность этого цилиндра сдвигается (закручивается) относительно внутренней в азимутальном направлении (азимутальное кручение) заданным угловым перемещением или усилием. Затем цилиндр охлаждается, в результате чего осуществляется прямой фазовый переход, и сдвиговые усилия, приложенные к внешней поверхности (заданные или соответствующие заданному угловому перемещению этой поверхности), снимаются. Цилиндр находится в мартенситном состоянии, в деформированной за счет возникших в нем фазовых деформаций, но устойчивой форме. Включается внешнее магнитное поле, и цилиндр нагревается до температуры, превышающей температуру окончания обратного фазового перехода, соответствующую случаю отсутствия механических и магнитных внешних полей. Но обратный фазовый переход при этом не происходит, так как приложенное магнитное поле увеличивает температуру начала этого перехода. При уменьшении магнитного поля до нуля осуществляется обратный фазовый переход, и закрученная в азимутальном направлении внешняя поверхность цилиндра раскручивается обратно в свое первоначальное положение. 1. Основные соотношения Кинематические соотношения. В рамках конечных деформаций введем согласно [18, 19] следующие понятия и обозначения, используемые в рассматриваемом процессе: и - начальная (недеформированная) и текущая конфигурации, - градиент места, переводящий конфигурацию и , С - мера деформации Коши-Грина, ; Е - тензор деформации Коши-Грина, ; g - метрический (единичный) тензор. Согласно [20] введем в рассмотрение промежуточную конфигурацию близкую к текущей, и F* - градиент места, переводящий конфигурацию в . Градиент места, мера и тензор деформации Коши-Грина относительно промежуточной конфигурации представляются в виде [20] (1) где - градиент относительно промежуточного состояния вектора перемещений из конфигурации в близкую текущую ; е - тензор малых деформаций относительно промежуточной конфигурации, ; - малый положительный параметр, характеризующий близость конфигураций и . Значения всех величин, помеченных здесь и ниже в статье «звездочкой», относятся к промежуточной конфигурации. В результате предельного перехода при стремлении промежуточной конфигурации к текущей () из соотношений (1) следует, что (2) где , D - тензор деформации скорости, , совпадающий в данном случае с тензором скорости деформаций . Определяющее уравнение. Из эквивалентных форм представления определяющих соотношений для простого материала, удовлетворяющих принципу объективности [18], используем форму (3) в которой - функция отклика материала на упругую деформацию - совпадает со вторым (симметричным) тензором напряжений Пиола-Кирхгоффа: , - упругий потенциал. В соотношении (3) T - тензор истинных напряжений, - третий инвариант F, определяющий относительное изменение объёма; - абсолютная температура. Из соотношения (3) с учетом (1) вытекает следующее представление второго тензора напряжений Пиола-Кирхгоффа, записанного относительно промежуточной конфигурации с точностью до линейных по слагаемых: (4) где и - тензоры четвертого ранга, определяющие отклик материала на приращение упругих деформаций и соответственно. Эти тензоры представляются следующим образом: где - позиционное скалярное умножение тензора второго ранга слева на третий базисный вектор тензора четвертого ранга [21, 22]. В результате предельного перехода в соотношении (4) получаем, что где тензоры и имеют тот же вид, что и , но с величинами, определенными не в промежуточной, а в текущей конфигурации. Уравнения магнитного поля. Как уже отмечалось, внутреннее магнитное поле, наведенное в образце конечных размеров внешним магнитным полем, возмущает последнее. В результате внутреннее поле изменяется тоже. В отсутствие электрического тока изменения внутреннего и внешнего полей описывается зависимостью [26] где Н и - напряженности возмущенного и приложенного магнитных полей; - оператор Гамильтона в текущей конфигурации; - скалярная функция, характеризующая возмущение. Вектор индукции магнитного поля В определяется в теле текущая конфигурация) соотношением а в пространстве вне тела текущая конфигурация) - соотношением Здесь - магнитная постоянная - вектор намагниченности тела, определяемый соотношением в котором - магнитная восприимчивость. Вектор индукции магнитного поля В удовлетворяет условию в результате чего имеем два уравнения и позволяющих определить функцию вне тела и в теле. Причем при и на поверхности S, разделяющей тело и его внешность, и где N - внешняя единичная нормаль к поверхности объема в текущей конфигурации, а Из последнего соотношения следует, что вектор напряженности магнитного поля Н испытывает скачок в нормальной составляющей на поверхности раздела S. Внешнее магнитное поле и магнитное поле, наведенное в образце конечных размеров, генерируют объемные (пондеромоторные) K и поверхностные Q силы в дополнение к внешним силам другой физической природы: где 2. Описание фазового перехода Для описания кинематики фазового перехода за основу принята модель, предложенная для полимеров с памятью формы в работах [23-25] и адаптированная к сплавам с памятью формы в работе [26]. Поскольку в ФСПФ имеет место фазовый переход аустенит-мартенсит, определим кинематику процесса градиентами места для аустенитной фазы и для мартенситной. Тогда соотношения (2) запишутся в следующем виде: Предполагая аддитивность скоростей изменения упругих (аустенитных и мартенситных), температурных и фазовых деформаций, представим кинематическое уравнение (аналогично [26]) в следующем виде: (5) Здесь и - параметры, определяющие долю мартенситной и аустенитной фаз в объёме материала (); и - скорости изменения фазовых и температурных деформаций. В работах [27-33] предложены различные соотношения, описывающие возникновение и развитие фазовых деформаций. Ниже будет использоваться модель, описанная в [30-33]. Зависимость от может быть аппроксимирована различными соотношениями. Согласно работе А.А. Мовчана [30], Здесь - температуры начала и завершения прямого и обратного мартенситных превращений в свободном от напряжений и от магнитного поля материале; - эти же температуры в нагруженном материале и/или в приложенном магнитном поле, которые определяются, согласно соотношению Клаузиуса-Клапейрона [26], следующим образом: Здесь и - критические температуры прямого и обратного фазовых переходов в нагруженном и ненагруженном материале; Н - вектор напряженности магнитного поля; , и - выделившееся (поглощенное) тепло, изменение фазовых деформаций и намагниченности в процессе фазового перехода. Предложенные в [31] соотношения, описывающие изменения фазовых деформаций, обобщаются в работе [26] на конечные деформации следующим образом: Здесь - параметры материала; - девиатор тензора ; и - значения параметра мартенситной фазы и фазовой деформации в начальной точке процесса обратного превращения. Скорость изменения температурных деформаций зададим, обобщая закон линейного температурного расширения на конечные деформации: где - коэффициент линейного температурного расширения. В работах [26, 34] рассмотрены требования, которые накладывает термодинамика на процесс фазового перехода в ФСПФ. В результате термодинамического анализа получено следующее равенство: , в котором , , и - упругие потенциалы аустенитной и мартенситной фаз. 3. Постановка краевой задачи Для ФСПФ необходимо решать связанную задачу: магнитное поле вызывает деформирование тела, которое, в свою очередь, искажает деформирующее это тело магнитное поле. Для решения такой краевой задачи используем вариационную постановку [26, 34]. Вариационное уравнение Лагранжа строится с применением стандартной процедуры Галеркина к уравнениям равновесия в текущей конфигурации и силовым граничным условиям действующим на части текущей поверхности тела с учетом кинематических граничных условий заданных на оставшейся части текущей поверхности тела Вариационное уравнение для магнитного поля строится совершенно аналогичным образом, с применением стандартной процедуры Галеркина к уравнениям для магнитного поля и выписанным в конце первого раздела. Поскольку до решения задачи поверхность и объем тела в текущем состоянии неизвестны, то, используя соотношения, связывающие объемы текущего и начального состояний, а также соотношения, связывающие поверхности, ограничивающие эти объемы, вариационную постановку для краевой задачи можно привести к начальной конфигурации: (6) Здесь n - внешняя единичная нормаль к поверхности тела с объемом в начальной конфигурации; - объем окружающей тело внешней среды в начальной конфигурации; u - вектор перемещений из начальной конфигурации в текущую; - оператор Гамильтона в начальной конфигурации; - символ вариации. Для описания упругого поведения аустенитной и мартенситной фаз используется упрощенный закон Синьорини [18]: где и - параметры материала, имеющие смысл параметра Ламе и модуля сдвига линейной теории упругости; - первый главный инвариант соответствующего тензора. Закон намагничивания зададим в виде , предполагая, что для магнитной восприимчивости справедливо соотношение , при этом и - магнитные восприимчивости мартенситной и аустенитной фаз, которые задаются согласно формуле Фрёлиха-Кенелли [35] следующим образом: Здесь - намагниченность насыщения; - постоянная материала. 4. Решение задачи о кручении полого цилиндра Приведённые выше соотношения использовались для решения задачи о кручении длинного полого цилиндра (плоская деформация) из сплава Гейслера. В начальный момент времени образец из ФСПФ в виде длинного полого цилиндра находится в аустенитном состоянии. Внутренняя поверхность цилиндра закреплена, а на внешней поверхности неизменного радиуса задаются перемещения или усилия, приводящие к её сдвигу относительно внутренней в азимутальном направлении (азимутальное кручение). Деформированный таким образом цилиндр охлаждается (при фиксированном перемещении/усилии на внешней поверхности) в интервале температур прямого фазового перехода, что приводит к накоплению фазовых деформаций. В мартенситном состоянии внешняя поверхность цилиндра разгружается (снимаются приложенные к ней заданные усилия или усилия, соответствующие заданному угловому перемещению этой поверхности), при этом конфигурацию образца определяют накопленные в нем фазовые деформации. Последующее снятие фазовых деформаций в образце возможно двумя способами: 1. Производится нагрев в интервале температур обратного фазового перехода, и образец возвращается к своей начальной ненагруженной конфигурации. 2. В мартенситном состоянии осуществляются включение магнитного поля и нагрев образца до температуры, при которой материал перешел бы в аустенитное состояние, будучи ненамагниченным. Однако приложенное магнитное поле сдвигает в соответствии с законом Клаузиуса-Клапейрона температуры обратного фазового перехода, и образец остается в мартенситном состоянии. Обратный фазовый переход и восстановление начальной конфигурации будет происходить в условиях постоянной температуры при снятии магнитного поля. Отметим, что когда внешнее магнитное поле постоянно и направлено вдоль оси бесконечно длинного цилиндра, образец намагничивается однородно, при этом магнитное поле внутри цилиндра равно внешнему: . Действительно, в единичном базисе цилиндрической системы координат функция зависит только от ρ, и уравнения, представленные в конце первого раздела и описывающие поведение магнитного поля как в цилиндре, так и вне его, принимают вид Отсюда следует, что где и - константы; - полость цилиндра, - цилиндр; - внешность цилиндра. Используя условия стыковки этих трех решений на внешней и внутренней поверхностях цилиндра, условия нормального разрыва на этих поверхностях и условие затухания до нуля при приходим к тому, что т.е. и в теле, и в пространстве. Отсюда следует, что в теле а в окружающем пространстве Поэтому скачок вектора магнитной индукции происходящий на внутренней и внешней поверхностях цилиндра, направлен вдоль этих поверхностей. Решение, приведенное выше, является точным аналитическим решением уравнений магнитостатики, а значит, и второго вариационного уравнения в (6). Поэтому в данном случае это уравнение уже не нужно решать, в отличие от задачи, рассмотренной в работе [26]. Поставленная в начале этого раздела задача решалась численно методом конечных элементов с помощью пакета FEnicCS (http://fenicsproject.org). Использовалась процедура пошагового нагружения для всех этапов решения, для чего осуществлена линеаризация вариационных уравнений. На каждом шаге задавались приращения перемещения/усилия на внешней поверхности, температуры или внешнего магнитного поля. Считалось, что распределение температуры однородно по образцу, т.е. он выдерживался необходимое для этого время при заданном на шаге приращении температуры окружающей среды. Внутренний радиус цилиндра r, внешний радиус R. Задача решалась в цилиндрической системе координат , при этом ось z совпадала с осью цилиндра. Ненулевые компоненты вектора перемещений зависели только от координаты ρ: и Для численного расчета использовались следующие параметры реального материала, приведенные или определенные в [26]: модули упругости , ; коэффициент Пуассона для аустенитной и мартенситной фазы , характерные температуры фазового перехода при отсутствии напряжений и магнитного поля: , , , ; параметры материала, описывающие эволюцию фазовых деформаций, , , ; коэффициент линейного теплового расширения ; магнитные константы , , , На рисунках показаны результаты численного расчета для цилиндра с внутренним радиусом и внешним радиусом Охлаждение после кручения в аустенитном состоянии происходило либо при фиксированном перемещении на внешней поверхности , либо при соответствующем ему усилии (в аустенитном состоянии) При этом Отметим, что результаты, полученные при решении на первом этапе (только упругие деформации в аустенитном состоянии) в рамках конечных деформаций, практически совпадают с аналитическим решением, которое можно получить при решении задачи в рамках малых деформаций. Для рассматриваемой задачи об азимутальном кручении полого цилиндра при плоской деформации единственной ненулевой компонентой вектора перемещений в цилиндрической системе координат будет при этом граничные условия запишутся в виде Используя соотношения, представленные в работах [36, 37], получаем Из уравнения равновесия следует, что Решением полученного дифференциального уравнения будет функция в которой и находятся из граничных условий. В результате получаем, что , где Тогда и На рис. 1 представлено сечение цилиндра и полученные распределения перемещений по радиусу (умноженные на 20) для задачи с заданными перемещениями на внешней поверхности (показаны красными линиями) и с заданными усилиями (показаны синими линиями). Упругое решение в аустенитном состоянии соответствует пути O-A перемещения точки на внешней поверхности цилиндра. Рис. 1. История перемещения точек сечения цилиндра Fig. 1. Displacement history of the points of the cylinder section Путь A-B получается после прямого фазового перехода при фиксированном перемещении внешней поверхности и последующем полном снятии усилия, соответствующего этому фиксированному перемещению в мартенситном состоянии (силовая разгрузка). Путь B-O соответствует обратному фазовому переходу в вышерассмотренном процессе. Путь A-C получается после прямого фазового перехода при фиксированном усилии на внешней поверхности. Последующая полная силовая разгрузка этой поверхности в мартенситном состоянии соответствует пути C-D, а D-O - путь при обратном фазовом переходе в таком процессе. Зависимость от в виде обычного графика представлена на рис. 2. Здесь красной сплошной линией показано распределение углового перемещения по радиусу при кручении в аустенитном состоянии, заданном перемещением или соответствующем ему усилием (вначале возникают только упругие деформации), красной пунктирной - распределение углового перемещения после прямого фазового перехода при фиксированном перемещении внешней поверхности , и последующей полной силовой разгрузке этой поверхности в мартенситном состоянии. Синей сплошной линией представлено распределение углового перемещения после прямого фазового перехода при фиксированном усилии на внешней поверхности синей пунктирной - распределение углового перемещения после прямого фазового перехода при фиксированном усилии на внешней поверхности и последующей полной силовой разгрузке этой поверхности в мартенситном состоянии. В процессе прямого фазового перехода при фиксированном усилии фазовые деформации растут. В силу того что упругие деформации в этом процессе должны оставаться неизменными, полные деформации тоже растут, все время превышают фазовые на величину упругих деформаций, и в конце такого процесса прямого фазового перехода перемещение соответствует точке С на рис. 2. При силовой разгрузке снимаются упругие деформации, полные деформации уменьшаются до фиксированных фазовых, образец раскручивается, и перемещение соответствует точке D на рис. 2. Рис. 2. Распределение перемещений по радиусу Fig. 2. Distribution of displacements along the radius Процесс прямого мартенситного превращения при фиксированном перемещении внешней поверхности цилиндра эквивалентен процессу сдвига этой поверхности усилием, но не постоянным, как в предыдущем случае, а изменяющимся по вполне определенному закону. Так как полные деформации фиксированы, а фазовые растут, упругие деформации уменьшаются, и это усилие затухает (релаксирует) до нуля при не полностью завершенном процессе прямого мартенситного превращения, Приращение фазовых деформаций при этом сохраняет свой знак (продолжается ориентированное превращение), хотя и становится меньше, чем в предыдущий момент времени. Последующий процесс охлаждения осуществляется уже при отрицательном сдвиговом усилии, так как накапливающиеся фазовые деформации превышают полные фиксированные деформации. Это отрицательное сдвиговое усилие будет только нарастать, пока его влияние на приращение фазовых деформаций не сравнится с влиянием накопленных фазовых деформаций (см. соотношение для в конце разд. 2). Поэтому при силовой разгрузке, когда полные деформации становятся равными фазовым, перемещение испытывает скачок из точки А в точку В и образец докручивается, реализуя эти накопленные фазовые деформации. В результате этого процесса угловое перемещение внешней поверхности цилиндра в разгруженном мартенситном состоянии (точки В и D) возрастает по сравнению с заданным в 10 раз для перехода с фиксированным угловым перемещением и в 13 раз для перехода с фиксированным касательным усилием за счет возникновения фазовых деформаций, которые сохраняются и в разгруженном мартенситном состоянии. На рис. 3 представлена зависимость углового перемещения на внешней поверхности от температуры (прямой фазовый переход при охлаждении при фиксированном усилии). Начальное и конечное значения перемещения отмечены точками A и C, которые соответствуют тем же точкам на рис. 1 и 2. При уменьшении температуры происходит дополнительное к упругому закручивание цилиндра за счет возникающих фазовых деформации. Рис. 3. Зависимость перемещения на внешней поверхности от температуры (охлаждение при фиксированном усилии) Fig. 3. Dependence of outer surface displacement on temperature (cooling under fixed force) На рис. 4 представлены зависимости углового перемещения на внешней поверхности цилиндра от приложенного магнитного поля. Начальные значения , помеченные точками B и D, соответствуют приложенному магнитному полю в (при этом магнитная индукция равна - размерность величины в тесла) и нагреванию цилиндра до температуры выше , соответствующей концу обратного фазового перехода мартенсит-аустенит при отсутствии силового и магнитного полей. Значения перемещений в точках B и D почти не отличаются от перемещений разгруженного состояния, помеченных на рис. 1 и 2 теми же буквами, так как температурные деформации и деформации, возникающие от приложенного магнитного поля, малы по сравнению с уже существующими в теле. При уменьшении магнитного поля до нуля образец полностью раскручивается (путь показан стрелками; красная линия соответствует задаче, в которой прямой фазовый переход аустенит-мартенсит осуществлялся при фиксированном перемещении, синяя - задаче, в которой прямой фазовый переход аустенит-мартенсит осуществлялся при фиксированном усилии). Этот рисунок демонстрирует возможность осуществления обратного фазового перехода при постоянной температуре, превышающей , при снятии магнитного поля, приложенного в мартенситном состоянии. На рис. 5 красной линией показано распределение сдвиговой компоненты тензора истинных напряжений по радиусу при кручении в аустенитном состоянии (упругое решение), а синей - распределение напряжений после прямого фазового перехода (в мартенситном состоянии) при фиксированном перемещении на внешней поверхности. Рис. 4. Раскручивание образца при снятии магнитного поля Fig. 4. Untwisting of the sample when removing the magnetic field Видно, что после прямого фазового перехода напряжение становится отрицательным по всему радиусу цилиндра, так как возникающие фазовые деформации при фиксированном перемещении на внешней поверхности приводят к отрицательным упругим деформациям (об этом говорилось при обсуждении рис. 2). Точки соответствуют сдвиговым усилиям на внешней и внутренней поверхностях. Рис. 5. Распределение сдвигового усилия по радиусу Fig. 5. Distribution of shear force along the radius На рис. 6 показаны зависимости сдвигового усилия на внешней (синей линией) и внутренней (красной линией) поверхностях от температуры (прямой фазовый переход при охлаждении при фиксированном на внешней поверхности перемещении и последующей разгрузке). Напряжения, отмеченные точками , соответствуют тем же точкам на рис. 5. При уменьшении температуры сдвиговые усилия уменьшаются и становятся отрицательными. При последующем снятии нагрузки на внешней поверхности (, соответствующее точке b, стремится к нулю) при фиксированной внутренней поверхности, на внутренней поверхности, соответствующее точке d, тоже стремится к нулю. На рис. 7 представлено распределение сдвиговой компоненты тензора фазовых деформаций по радиусу. Красная линия - результат решения, полученный при фазовом переходе при фиксированном перемещении на внешней поверхности, синяя - при фиксированном усилии. Рис. 6. Зависимость сдвигового усилия от температуры Fig. 6. Dependence of shear force on temperature Из рисунка видно, что охлаждение при фиксированном усилии дает большие фазовые деформации, чем охлаждение при фиксированном перемещении, поскольку во втором случае, как уже говорилось при обсуждении рис. 2, происходит уменьшение напряжений в процессе охлаждения по всему образцу (см. рис. 5), что замедляет развитие фазовых деформаций. Рис. 7. Распределение фазовых деформаций по радиусу Fig. 7. Distribution of phase strains along the radius Рис. 8. Распределение осевой составляющей Bz вектора магнитной индукции Fig. 8. Distribution of the axial component of the magnetic inductance vector На рис. 8 показан разрыв в осевой (касательной к поверхности цилиндра) составляющей вектора магнитной индукции В для рассматриваемой задачи, о чем говорилось в начале этого раздела. Здесь . Все остальные составляющие вектора В равны нулю. Разрыва в нормальной к поверхности цилиндра составляющей вектора напряженности магнитного поля Н в данной задаче не происходит, что уже отмечалось ранее. Заключение В рамках модели, описывающей поведение сплавов с памятью формы при конечных деформациях, решена краевая задача о накоплении и снятии фазовых деформаций в длинном полом цилиндре из ферромагнитного сплава с памятью формы (сплаве Гейслера). Показано, что величина накопленных фазовых деформаций в процессе прямого перехода зависит от граничных условий. Полное снятие накопленных деформаций может происходить как при нагревании образца в интервале температур обратного фазового перехода, так и при постоянной температуре, превышающей температуру конца обратного фазового перехода мартенсит-аустенит при отсутствии силового и магнитного полей, в случае снятия магнитного поля, предварительно приложенного в мартенситном состоянии. Последнее продемонстрировано на примере раскручивания внешней поверхности цилиндра из ферромагнитного материала, предварительно закрученной в аустенитном состоянии и значительно увеличившей угол закрутки в результате фазового перехода в мартенситное состояние. Отметим, что с помощью описанной модели в работе [38] представлено решение задачи о натягивания предварительно сдвинутой в осевом направлении внешней поверхности цилиндра из ферромагнитного сплава с памятью формы как при нагревании образца в интервале температур обратного фазового перехода, так и при постоянной температуре при снятии магнитного поля, предварительно приложенного в мартенситном состоянии.Об авторах
А А Роговой
Институт механики сплошных сред УрО РАН
О С Столбова
Институт механики сплошных сред УрО РАН
Список литературы
- Ферромагнетики с памятью формы / А.Н. Васильев [и др.] // Успехи физических наук. - 2003. - Т. 173, № 6. - С. 577-608. doi: 10.3367/UFNr.0173.200306a.0577
- Магнитные сплавы с памятью формы: фазовые переходы и функциональные свойства / В.Д. Бучельников [и др.] // Успехи физических наук. - 2006. - Т. 176, № 8. - С. 900-906. doi: 10.3367/UFNr.0176.200608j.0900
- Roubicek T., Stefanell U. Magnetic shape-memory alloys: thermomechanical modelling and analysis // Continuum Mechanics and Thermodynamics. - 2014. - Vol. 26. - No. 6. - P. 783-810. doi: 10.1007/s00161-014-0339-8
- Lagoudas D.C. Shape Memory Alloys. Modeling and Engineering Applications. - New York: Springer Science+Business Media, 2008. - 455 p. doi: 10.1007/978-0-387-47685-8
- Haldar K., Lagoudas D.C., Karaman I. Magnetic field-induced martensitic phase transformation in magnetic shape memory alloys: Modeling and experiments // Journal of the Mechanics and Physics of Solids. - 2014. - Vol. 69. - P. 33-66. doi: 10.1016/j.jmps.2014.04.011
- Kiefer B., Lagoudas D.C. Modeling the coupled strain and magnetization response of magnetic shape memory alloys under magnetomechanical loading // Journal of Intelligent Material Systems and Structures. - 2009. - Vol. 20. - P. 143-170. doi: 10.1177/1045389X07086688
- A three-dimensional model of magneto-mechanical behaviors of martensite reorientation in ferromagnetic shape memory alloys / X. Chen, Z. Moumni, Y. He, W. Zhang // Journal of the Mechanics and Physics of Solids. - 2014. - Vol. 64. - P. 249-286. doi: 10.1016/j.jmps.2013.11.005
- Thermal effects on high-frequency magnetic-field-induced martensite reorientation in ferromagnetic shape memory alloys: An experimental and theoretical investigation / S. Zhang, X. Chen, Z. Moumni, Y. He // International Journal of Plasticity. - 2018. - Vol. 108. - P. 1-20. doi: 10.1016/j.ijplas.2018.04.008
- Peng Q., He Y.J., Moumni Z. A phase-field model on the hysteretic magneto-mechanical behaviors of ferromagnetic shape memory alloy // Acta Materialia. - 2015. - Vol. 88. - P. 13-24. doi: 10.1016/j.actamat.2015.01.044
- Shape memory effect due to magnetic field induced thermoelastic martensitic transformation in polycrystalline Ni-Mn-Fe-Ga alloy / A.A. Cherechukin [et al.] // Physics Letters A. - 2001. - Vol. 291. - No. 2-3. - P. 175-183. doi: 10.1016/S0375-9601(01)00688-0
- Training of the Ni-Mn-Fe-Ga ferromagnetic shape-memory alloys due cycling in high magnetic field / A.A. Cherechukin [et al.] // Journal of Magnetism and Magnetic Materials. - 2003. - Vol. 258-259. - P. 523-525. doi: 10.1016/S0304-8853(02)010
- Наблюдение одностороннего эффекта памяти формы, обусловленного магнитоиндуцированным мартенситным фазовым переходом в сплаве Ni-Mn-Fe-Ga [Электронный ресурс] / В.Г. Шавров [и др.] // Журнал радиоэлектроники. - 2001. - № 4. - URL: http://jre.cplire.ru/jre/may01/1/text_r.html (дата обращения: 25.03.2019).
- Малыгин Г.А. Теория эффектов магнитной памяти формы и псевдоупругой деформации в сплавах Ni-Mn-Ga // Физика твердого тела. - 2009. - Т. 51, № 8. - С. 1599-1603.
- Малыгин Г.А. Моделирование эффекта магнитной памяти формы в пленках сплава Ni2MnGa // Журнал технической физики. - 2007. - Т. 77, № 2. - С. 136-139.
- Bustamante R., Dorfman A., Ogden R.W. Universal relations in isotropic nonlinear magnetoelasticity // The Quarterly Journal of Mechanics and Applied Mathematics. - 2006. - Vol. 59. - No. 3. - P. 435-450. doi: 10.1093/qjmam/hbl010
- Bustamante R., Dorfman A., Ogden R.W. A nonlinear magnetoelastic tube under extension and inflation in an axial magnetic field: numerical solution // Journal of Engineering Mathematics. - 2007. - Vol. 59. - No. 1. - P. 139-153. doi: 10.1007/s10665-006-9088-4
- Bustamante R., Dorfman A., Ogden R.W. Numerical solution of finite geometry boundary-value problems in nonlinear magnetoelasticity // International Journal of Solids and Structures. - 2011. - Vol. 48. - No. 6. - P. 874-883. doi: 10.1016/j.ijsolstr.2010.11.021
- Лурье А.И. Нелинейная теория упругости. - М.: Наука, 1980. - 512 с.
- Трусделл К. Первоначальный курс рациональной механики сплошных сред. - М.: Мир, 1975. - 585 с.
- Rogovoy A.A. Formalized approach to construction of the state equations for complex media under finite deformations // Continuum Mechanics and Thermodynamics. - 2012. -Vol. 24. - No. 2. - P. 81-114. doi: 10.1007/s00161-011-0220-y
- Роговой А.А. Дифференцирование скалярных и тензорных функций тензорного аргумента // Вестн. Перм. техн. ун-та. Динамика и прочность машин. - 2001. - № 2. - С. 83-90.
- Rogovoy A.A. Differentiation of scalar and tensor functions of tensor argument // IOSR Journal of Mathematics. - 2019. - Vol. 15. - No. 2. - P. 1-20. doi: 10.9790/5728-1502010120
- A thermodynamically-consistent 3D constitutive model for shape memory polymers / M. Baghani, R. Naghdabadi, J. Arghavani, S. Sohrabpour // International Journal of Plasticity. - 2012. - Vol. 35. - P. 13-30. doi: 10.1016/j.ijplas.2012.01.007
- Baghani M., Naghdabadi R., Arghavani J. A large deformation framework for shape memory polymers: Constitutive modeling and finite element implementation // Journal of Intelligent Material Systems and Structures. - 2013. - Vol. 24. - No. 1. - P. 21-32. doi: 10.1177/1045389X12455728
- Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modeling / Y. Liu [et al.] // International Journal of Plasticity. - 2006. - Vol. 22. - P. 279-313. doi: 10.1016/j.ijplas.2005.03.004
- Rogovoy A., Stolbova O. Modeling the magnetic field control of phase transition in ferromagnetic shape memory alloys // International Journal of Plasticity. - 2016. - Vol. 85. - P. 130-155. doi: 10.1016/j.ijplas.2016.07.006
- Theoretical and numerical modeling of shape memory alloys accounting for multiple phase transformations and martensite reorientation / F. Auricchio [et al.] // International Journal of Plasticity. - 2014. - Vol. 59. - P. 30-54. doi: 10.1016/j.ijplas.2014.03.008
- Cisse C., Zaki W., Zineb T.B. A review of constitutive models and modeling techniques for shape memory alloys // International Journal of Plasticity. - 2016. - Vol. 76. - P. 244-284. doi: 10.1016/j.ijplas.2015.08.006
- Analytical model for the torsional response of superelastic shape memory alloy circular sections subjected to a loading-unloading cycle / N.V. Viet [et al.] // International Journal of Solids and Structures. - 2019. - Vol. 156-157. - P. 49-60. doi: 10.1016/j.ijsolstr.2018.08.001
- Мовчан А.А., Чжо Ту Я. Решение начально-краевых задач о прямом и обратном превращении в рамках нелинейной теории деформирования сплавов с памятью формы // Механика композиционных материалов и конструкций. - 2007. - Т. 13, № 4. - С. 452-468.
- Мовчан А.А., Шелымагин П.В., Казарина С.А. Определяющие уравнения для двухэтапных термоупругих фазовых превращений // ПМТФ. - 2001. - Т. 42, № 5. - С. 152-160.
- Мовчан А.А., Сильченко Л.Г. Аналитическое решение связной задачи об устойчивости пластины из сплава с памятью формы при обратном мартенситном превращении // Изв. РАН. Механика твердого тела. - 2004. - № 5. - С. 164-178.
- Мовчан А.А., Чжо Ту Я. Решение связанной термоэлектромеханической задачи для стержня из сплава с памятью формы в рамках теории нелинейного деформирования этих материалов // Механика композиционных материалов и конструкций. - 2008. -Т. 14, № 3. - С. 443-460.
- Роговой А.А., Столбова О.С. Конечные деформации в сплавах и полимерах с памятью формы // Учён. зап. КнАГТУ. - 2018. - Т. 1, № 3 (35). - С. 6-17.
- Бозорт Р. Ферромагнетизм. - М.: Изд-во иностр. лит., 1956. - 784 с.
- Лурье А.И. Теория упругости. - М.: Наука, 1970. - 940 с.
- Kachanov M., Shafiro B., Tsukrov I. Handbook of Elasticity Solutions. - Dordrecht, Springer Science+Business Media, 2003. - 324 p. doi: 10.1007/978-94-017-0169-3
- Rogovoy A.A., Stolbova O.S. Numerical simulation of the phase transition control in a cylindrical sample made of ferromagnetic shape memory alloy // Computation. - 2019. - Vol. 7. - No. 3. - P. 38.1-38.15. doi: 10.3390/computation7030038