УДК 539.5
DOI: 10.18698/2308-6033-2022-1-2142
Моделирование эффекта мартенситной неупругости в плоских композиционных пружинах, выполненных из сплава с эффектом памяти формы
© А.А. Бутрина, С.М. Ганыш, С.С. Гаврюшин МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
Разработана математическая модель плоской композиционной пружины, работающей в области малых перемещений, обладающей эффектом памяти формы. Для описания поведения материала применен феноменологический подход, основанный на диаграмме фазовых переходов. В качестве модели материала выбрана феноменологическая модель Бринсон. Предложены соотношения для определения приведенных механических характеристик монослоя, эквивалентного композиционному слою в пружине. Эффект памяти учтен с помощью дополнительного внутреннего силового фактора — момента памяти формы при изгибе, являющегося результатом ориентации мартенсита по сечению. Для получения упругой характеристики плоской пружины использован конечно-элементный подход, в котором момент памяти формы выступает в качестве дополнительной узловой нагрузки. Представлен алгоритм построения зависимости между изгибающим моментом в сечении и моментом памяти формы для изотермического нагружения в зоне стабильности мартенсита. Получены зависимости момента памяти формы при изгибе от изгибающего момента для различных конфигураций поперечного сечения плоской пружины.
Ключевые слова: сплав с эффектом памяти формы, плоская пружина, композиционная пружина
Введение. В настоящее время наблюдается повышенный интерес к разработке интеллектуальных материалов и созданию на их основе конструкций с адаптивными возможностями. Он особенно активно проявляется при совершенствовании применяющихся в авиационно-космической отрасли, робототехнике и медицине сенсорных и актю-аторных компонентов датчиков, исполнительных механизмов и микроконтроллеров, имеющих упругие элементы, выполненные из сплавов с эффектом памяти формы [1, 2].
Рис. 1. Примеры поперечного сечения плоской пружины
Среди геометрических форм таких изделий широко распространены плоские пружины, которые могут быть выполнены в виде как монослоя, так и многослойной (композитной) конструкции [3-5] (рис. 1).
Цель работы — разработка алгоритма расчета плоских композиционных пружин, работающих в условиях мартенситной неупругости.
Математическая модель сплава с эффектом памяти формы.
Такие сплавы обладают рядом особенностей, в том числе способностью подвергаться большим обратимым деформациям и способностью восстанавливать первоначальную форму после термоциклиро-вания. Эффект мартенситной неупругости проявляется в результате накопления структурной деформации, сохраняющейся в процессе разгрузки, которая, однако, может быть снята после окончания термоцикла [6] (рис. 2).
Рис. 2. Зависимость напряжения с от структурной деформации е:
ст5, сту — постоянные параметры диаграммы фазовых переходов для сплава с эффектом памяти формы
Для практических расчетов наиболее удобны феноменологические модели макроуровня, основанные на диаграмме фазовых переходов сплава [7]. В процессе исследования была применена модифицированная модель Бринсон [8, 9], в которой для описания эффекта мартенситной неупругости в качестве внутренних переменных используются объемная доля хаотического (неориентированного мартенсита) £м и объемная доля ориентированного мартенсита £5, связанные соотношением баланса долей:
^ + £м =1. (1)
Соотношение, определяющее объемную долю ориентированного мартенсита, может быть записано в дискретной форме, что позволяет учитывать разгрузку и повторное нагружение:
ПРи а к+1 > а к
0 , если а к+1 < а^;
.(к+1) _ =5
ПРи ак+1 < ак
РМ5 (ак+1), если РМ5 (ак+1) > ^) и а5 < ак+1 < а/;
йк), если РМ5 (ак+1) < £5) и а5 < ак+1 < а/; (2) 1 , если а к+1 ^ а г.
V-
^(к+1) _ ),
где +1), #) — объемные доли ориентированного мартенсита на текущем и предыдущем шаге соответственно; ак+1, ак — напряжения на текущем и предыдущем шаге соответственно; РМ5 (а) — функция, аппроксимирующая фазовый переход:
РМ5 ( а)_ 2С08
а-а
п-
Г
а5 - а
Г
1
(3)
Для деформаций предполагается справедливым аддитивное разложение на упругую и неупругую составляющие:
а
8 _ 8упр + гЬ^5 _ Е + 8^5,
(4)
где 8Ь — механическая характеристика материала, представляющая
собой максимальную деформацию, возникающую в результате полного перехода мартенсита из хаотического состояния в ориентированное; Е — модуль упругости.
Момент памяти формы. Для учета эффекта ориентации мартенсита вводится дополнительный внутренний силовой фактор: момент памяти формы при изгибе [9]. Он приводит к тем же изменениям кривизны осевой линии плоской пружины, что и дополнительные деформации, вызванные ориентацией мартенсита:
М5М _\аЕ8^^(У)У*А,
М + М 5МА
д^ _ изг изг
Е1и
(5)
(6)
где М^^ — момент памяти формы в сечении; Дда — изменение кривизны осевой линии; Мизг — изгибающий момент в сечении;
Е1изг — изгибная жесткость поперечного сечения, определяемая соотношением
£/шг = | Е(у)у2ёА.
(7)
Соотношения (5)-(7) получены для симметричного поперечного сечения в предположении справедливости гипотезы плоских сечений.
Приведенные механические характеристики слоя. Если в пределах слоя присутствуют материалы с различными механическими характеристиками, такой слой заменяется слоем из мономатериала. Для получения приведенных свойств монослоя используется двух-компонентная схема слоя (рис. 3) [10].
Е\, аар гь
Е, а8, с/, еь
Рис. 3. Переход от двухкомпонентного слоя к монослою
Приведенные модуль упругости ЕЕ и максимальная неупругая деформация 8Ь определяются из условий равенства нормальных усилий N и удлинений АЬ слоя, состоящего из двух материалов (двух подслоев в случае шпоночной модели) и эквивалентного монослоя:
N = оЪИ = о1ЪИ1 + о2ЪЬ2; АЬ = гЬ = 81Ь = в2 Ь,
(8)
где Ъ1, Ъ2, Ъ, Ъ\, Ь, Ь — геометрические размеры слоя и подслоев (см. рис. 3); а1, 81 — напряжение и деформация в подслое из сплава
с памятью формы; 02, 82 — напряжение и деформация в подслое из традиционного материала; о, 8 — напряжение и деформация в эквивалентном монослое.
Физические соотношения для подслоев и монослоя имеют следующий вид:
°1 = Е18+Е18^;
°2 = Е28;
о = Е 8 + Е 8 ь ^,
(9) (10)
где Е1, 8Ь — механические характеристики подслоя с памятью формы; Е2 — модуль упругости в подслое из традиционного материала.
Подставляя соотношения (9), (10) в соотношения (8) и проводя ряд преобразований, получаем формулы для определения приведенных модуля упругости Е и максимальной неупругой деформации в£ монослоя:
Е _ Е1/1 +Е2Л; (11)
8 _ Е1^1 8 (12) ^ Е1/1 + Е2/2 ^ У }
где /], /2 — объемные доли сплава с памятью формы и традиционного сплава соответственно в пределах рассматриваемого слоя, определяемые соотношениями
г _ Ъ1. 11 _ И ;
И (13)
г _ И2
И
Приведенные напряжения начала а 5 и окончания а у фазового
перехода определяем из условия равенства деформаций в монослое и подслое из сплава с памятью формы. В момент начала фазового перехода (объемная доля ориентированного мартенсита £5 = 0) деформация 8 5 в монослое и подслое из сплава с памятью формы находится из физических соотношений (9), (10):
85 _Е _Е. (14)
Е1 Е
Аналогично в момент окончания фазового перехода (объемная доля ориентированного мартенсита £5 =1) деформация, соответствующая окончанию фазового перехода,
а у а у
8у _~ЕТ + 8Ь _+ 81. (15)
Е Е
Выполнив преобразования, получаем выражения для приведенных напряжений начала и окончания фазового перехода:
ЕЕ
а ^ =—1^;
1
(16)
Численная реализация. Для численного определения момента памяти формы сечение пружины дискретизируется (рис. 4), т. е. заменяется набором слоев. Предполагается, что напряженно-деформированное и фазовое состояния в пределах слоя постоянны.
Для построения одной точки зависимости между изгибающим моментом и моментом памяти формы действия выполняются в следующей последовательности:
1) задаем некоторое значение деформации втах на внешнем слое;
2) в предположении о линейном распределении деформаций по сечению определяем деформацию в каждом из слоев:
3) для каждого из слоев, если он состоит из двух материалов, вычисляем приведенные механические характеристики по соотношениям (11)-(15);
4) определяем составляющую момента памяти формы в каждом из слоев. Если слой не содержит сплав с памятью формы, то составляющая момента памяти формы М^мл равна нулю. Если слой выполнен из сплава с памятью формы или он композитный, то составляющая момента памяти формы определяется соотношением
у
Рис. 4. Дискретизация поперечного сечения плоской пружины
(17)
(18)
5) рассчитываем момент памяти формы по формуле
2п
(19)
6) находим изгибную жесткость сечения и изгибающий момент в сечении, соответствующий вычисленному моменту памяти формы и заданной деформации на внешнем слое:
2п
Е/изг =2 Е^Ъ, Ак; г=1
Мизг = Е/изгА^ -ММ = Е/изг ММ.
к
(20)
(21)
Упругая характеристика плоской пружины. Для построения упругой характеристики плоской пружины, работающей в области малых перемещений, применяется конечно-элементный подход [11, 12]. Каноническое уравнение метода конечных элементов для элемента из сплава с памятью формы имеет вид
Е/„
\(е)
Е3
12 6 Е -12 6Е
6 Е 4 Е2 -6 Е 2 Е2
-12 -6 Е 12 -6Е
6 Е 2 Е2 -6 Е 4 Е2
{ } = {(е)}
0
-М,
БМЛ
М,
0
8МЛ
(22)
Для расчета плоской пружины, работающей в области малых перемещений, можно также применять дифференциальное уравнение изогнутой оси стержня [11]:
Е/и
ё 2У ( 7 ) ёг 2
= М „
( г) + ММ ( г ).
(23)
Результаты численного моделирования. На основе предложенных соотношений и алгоритмов в среде МЛТЬЛВ был создан программный комплекс для расчета плоских пружин, работающих в условиях мартенситной неупругости. Расчет проводился для следующих параметров материала [13]:
Е1 = 20 000 МПа; Е2 = 1500 МПа; вЕ = 0,05; а £ = 25 МПа; а г = 78 МПа.
(24)
Зависимости момента памяти формы от изгибающего момента для различных типов поперечных сечений представлены на рис. 5.
Полученные зависимости для момента памяти формы позволяют моделировать эффект мартенситной неупругости в плоских пружинах. На рис. 6 представлено семейство упругих характеристик для консольной плоской пружины длиной 80 мм с поперечным сечением № 1 (рис. 5, а), нагруженной сосредоточенной нагрузкой.
500
400 — =
300
200
100
50 100 150
Изгибающий момент, Н-мм
| 1000 д
» 800
600
400
200
50 100 150 200 250 300 Изгибающий момент, Н-мм
Рис. 5. Зависимость момента памяти формы от изгибающего момента:
а — сечение № 1 (Ь = 10 мм, к = 2 мм); б — сечение № 2 (Ь = 10 мм, к = 2 мм, к = 0,1 мм); в — сечение № 3 (Ь = 10 мм, к = 2 мм)
. -Ь\ = 0,1 мм 0,2 мм
-Ai = 0,3 мм
О 0,5 1,0 1,5 2,0 2,5 3,0
Сила, Н
Рис. 6. Семейство упругих характеристик для плоской пружины
Справедливость полученных результатов была подтверждена сравнением с теоретическими и экспериментальными результатами, представленными в [3-5, 11-13].
Заключение. Разработана математическая модель плоской композиционной пружины, обладающей эффектом памяти формы, работающей при малых перемещениях. Получены соотношения для приведенных характеристик монослоя: модуля упругости, напряжений начала и окончания фазового перехода, максимальной неупругой деформации, позволяющие учитывать эффект ориентации мартенсита с помощью виртуального дополнительного внутреннего силового фактора — момента памяти формы при изгибе. Предложен алгоритм построения зависимости момента памяти формы от изгибающего момента для композиционного сечения. Написана пользовательская программа на языке MATLAB, с помощью которой получены результаты расчета плоских пружин с различными поперечными сечениями. Представленные в работе методика, алгоритм и соотношения могут быть рекомендованы для проектирования термодатчиков и предохранителей с исполнительным элементом в форме плоских пружин, обладающих эффектом памяти формы.
ЛИТЕРАТУРА
[1] Jani J.M., Leary M., Subic A., Gibson M.A. A review of shape memory alloy research, applications and opportunities. Materials & Design, 2015, vol. 56, pp. 1078-1113.
[2] Narahari P.S., Wanhill R.J.H. Aerospace materials and material technologies. Springer Science+BusinessMedia, 2017, vol. 1, pp. 1-14.
[3] Zhang Y., Zhao Y.-P. A study of composite beam with shape memory alloy arbitrarily embedded under thermal and mechanical loadings. Materials & Design,, 2007, vol. 28 (4), pp. 1096-1115.
[4] Rasid Z.A., Zahari R. The deflection of the shape memory alloy composite beams using finite element method. Applied Mechanics and Materials, 2014, vol. 695, pp. 135-138.
[5] Saraswat Y., Yadav S., Parihar H.S. The effect of shape memory alloy in composite beam. International Journal of Innovative Technology and Exploring Engineering, 2019, vol. 8, pp. 130-134.
[6] Тихомирова К.А. Изотермическое деформирование сплава с памятью формы в разных температурных интервалах. Случай одноосного растяжения. Механика композиционных материалов и конструкций, 2017, т. 23, № 2, c. 263-282.
[7] Cisse C., Zaki W., Zineb T.B. A review of modeling techniques for advanced effects in shape memory alloy behavior. Smart Materials and Structures, 2016, vol. 25 (10), pp. 1-36.
[8] Poorasadion S., Arghavani J., Naghdabadi R., Sohrabpour S. An improvement on the Brinson model for shape memory alloys with application to two-dimensional beam element. Journal of Intelligent Material Systems and Structures, 2013, vol. 25, pp. 1905-1920.
[9] Гаврюшин С.С., Ганыш С.М. Расчет элементов коммутационных и исполнительных устройств, выполненных из сплава с памятью формы. Проблемы машиностроения и надежности машин, 2019, № 7, c. 6-14.
[10] Кузьмин М.А., Лебедев Д.Л., Попов Б.Г. Прочность, жесткость, устойчивость элементов конструкций. Теория и практикум. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2012, 344 с.
[11] Гаврюшин С.С., Ганыш С.М. Построение упругих характеристик для плоских и винтовых цилиндрических пружин, выполненных из сплава с эффектом памяти формы. Сб. науч. трудов: Математическое моделирование и экспериментальная механика деформируемого твердого тела. Тверь, 2018, т. 2, с. 21-25.
[12] Гаврюшин С.С., Ганыш С.М. Применение метода конечных элементов для расчета плоских пружин, выполненных из сплавов с памятью формы. Сб. науч. трудов: Современные вопросы устойчивости, пластичности и ползучести в механике деформируемого твердого тела, Тверь, 2020, с. 94-105.
[13] Sayyaadi H., Zakerzadeh M.R., Salehi H. A comparative analysis of some one-dimensional shape memory alloy constitutive models based on experimental tests. Scientia Iranica, 2012, vol. 19 (2), pp. 249-257.
Статья поступила в редакцию 03.01.2022
Ссылку на эту статью просим оформлять следующим образом:
Бутрина А. А., Ганыш С.М., Гаврюшин С. С. Моделирование эффекта мартен-ситной неупругости в плоских композиционных пружинах, выполненных из сплава с эффектом памяти формы. Инженерный журнал: наука и инновации, 2022, вып. 1. http://dx.doi.org/10.18698/2308-6033-2022-1-2142
Бутрина Анастасия Андреевна — студентка кафедры «Прикладная механика» МГТУ им. Н.Э. Баумана. Область научных интересов: композиционные материалы, метод конечных элементов, робототехника. e-mail: butrina_99@mail.ru
Ганыш Святослав Мирославович — канд. техн. наук, доцент кафедры «Прикладная механика» МГТУ им. Н.Э. Баумана. Область научных интересов: метод конечных элементов, материалы с памятью формы. e-mail: s.ganysh@bmstu.ru
Гаврюшин Сергей Сергеевич — д-р техн. наук, заведующий кафедрой «Компьютерные системы автоматизации производства» МГТУ им. Н. Э. Баумана. e-mail: gss@bmstu.ru
Modeling the martensitic inelasticity effect in flat composite springs made of a shape memory alloy
© A.A. Butrina, S.M. Ganysh, S.S. Gavryushin Bauman Moscow State Technical University, Moscow, 105005, Russia
The purpose of the study was to develop a mathematical model of a flat composite spring operating in the region of small displacements and having a shape memory effect. To describe the behavior of the material, we applied a phenomenological approach based on the phase transition diagram; as the material model, we chose the Brinson phenomenological model. As a result, we proposed relationships for determining the reduced mechanical characteristics of a monolayer equivalent to a composite layer in the spring. The memory effect was taken into account by an additional internal force factor, i. e. the shape memory moment during bending, which resulted from the cross-section orientation of the martensite. To obtain the elastic characteristic of the spring, we used a finite element approach, in which the shape memory moment acts as an additional nodal load. Eventually, we introduced an algorithm for the dependence between the bending moment in the section and the shape memory moment for isothermal loading in the martensite stability zone and obtained the dependences of the shape memory moment during bending on the bending moment for various configurations of the cross-section of the spring.
Keywords: shape memory alloy, flat spring, composite spring
REFERENCES
[1] Jani J.M., Leary M., Subic A., Gibson M.A. A review of shape memory alloy research, applications and opportunities. Materials & Design, 2015, vol. 56, pp. 1078-1113.
[2] Narahari P.S., Wanhill R.J.H. Aerospace materials and material technologies.
Springer Science+BusinessMedia, 2017, vol. 1, pp. 1-14.
[3] Zhang Y., Zhao Y.-P. A study of composite beam with shape memory alloy arbitrarily embedded under thermal and mechanical loadings. Materials & Design, 2007, vol. 28 (4), pp. 1096-1115.
[4] Rasid Z.A., Zahari R. The deflection of the shape memory alloy composite beams using finite element method. Applied Mechanics and Materials, 2014, vol. 695, pp. 135-138.
[5] Saraswat Y., Yadav S., Parihar H.S. The effect of shape memory alloy in composite beam. International Journal of Innovative Technology and Exploring Engineering, 2019, vol. 8, pp. 130-134.
[6] Tikhomirova K.A. Mekhanika kompozitsionnykh materialov i konstruktsii — Journal on Composite Mechanics and Design, 2017, vol. 23, no. 2, pp. 263-282.
[7] Cisse C., Zaki W., Zineb T.B. A review of modeling techniques for advanced effects in shape memory alloy behavior. Smart Materials and Structures, 2016, vol. 25 (10), pp. 1-36.
[8] Poorasadion S., Arghavani J., Naghdabadi R., Sohrabpour S. An improvement on the Brinson model for shape memory alloys with application to two-dimensional beam element. Journal of Intelligent Material Systems and Structures, 2013, vol. 25, pp. 1905-1920.
[9] Gavryushin S.S., Ganysh S.M. Problemy mashinostroeniya i nadezhnosti mashin — Journal of Machinery Manufacture and Reliability, 2019, no. 7, pp. 6-14.
[10] Kuzmin M.A., Lebedev D.L., Popov B.G. Prochnost, zhestkost, ustoychivost elementov konstruktsii. Teoriya i praktikum [Strength, rigidity, stability of structural elements. Theory and practice]. Moscow, BMSTU Publ., 2012, 344 p.
[11] Gavryushin S.S., Ganysh S.M. Postroenie uprugikh kharakteristik dlya ploskikh i vintovykh tsilindricheskikh pruzhin, vypolnennykh iz splava s effektom pa-myati formy [Building elastic characteristics for flat and helical coil springs made of shape memory alloy]. In: Sb. nauch. tr. Matematicheskoe modeliro-vanie i eksperimentalnaya mekhanika deformiruemogo tverdogo tela [Collection of scientific papers Mathematical modeling and experimental mechanics of a de-formable solid body]. Tver, 2018, vol. 2, pp. 21-25.
[12] Gavryushin S.S., Ganysh S.M. Primenenie metoda konechnykh elementov dlya rascheta ploskikh pruzhin, vypolnennykh iz splavov s pamyatyu formy [Application of the finite element method for the calculation of flat springs made of shape memory alloys]. In: Sb. nauch. tr. Sovremennye voprosy ustoychivosti, plastichnosti i polzuchesti v mekhanike deformiruemogo tverdogo tela [Collection of scientific papers. Modern issues of stability, plasticity and creep in the mechanics of a deformable solid body]. Tver, 2020, pp. 94-105.
[13] Sayyaadi H., Zakerzadeh M.R., Salehi H. A comparative analysis of some one-dimensional shape memory alloy constitutive models based on experimental tests. Scientia Iranica, 2012, vol. 19 (2), pp. 249-257.
Butrina A.A., student, Department of Applied Mechanics, Bauman Moscow State Technical University. Research interests: composite materials, finite element method, robotics. e-mail: butrina_99@mail.ru
Ganysh S.M., Cand. Sc. (Eng.), Assoc. Professor, Department of Applied Mechanics, Bauman Moscow State Technical University. Research interests: finite element method, shape memory materials. e-mail: s.ganysh@bmstu.ru
Gavryushin S.S., Dr. Sc. (Eng.), Professor, Head of Department of Computer Systems of Manufacturing Automation, Bauman Moscow State Technical University. e-mail: gss@bmstu.ru