Проблемы машиностроения и надежности машин, 2019, № 7, стр. 6-14
РАСЧЕТ ЭЛЕМЕНТОВ КОММУТАЦИОННЫХ И ИСПОЛНИТЕЛЬНЫХ УСТРОЙСТВ, ВЫПОЛНЕННЫХ ИЗ СПЛАВОВ С ПАМЯТЬЮ ФОРМЫ
С. С. Гаврюшин 1, 2, *, С. М. Ганыш 1, **
1 Московский государственный технический университет имени Н.Э. Баумана
г. Москва, Россия
2 Институт машиноведения им. А.А. Благонравова
г. Москва, Россия
* E-mail: gss@bmstu.ru
** E-mail: s.ganysh-rk5@yandex.ru
Поступила в редакцию 22.07.2019
Принята к публикации 26.08.2019
Аннотация
В статье разработаны математические модели для стержневых элементов, плоских и винтовых цилиндрических пружин, выполненных из сплава с эффектом памяти формы. Математическая модель материала и алгоритм расчета построены на основе дифференциальных уравнений гибких стержней и диаграмм фазовых переходов. Для учета эффекта памяти формы в систему дифференциальных уравнений аддитивно вводятся фиктивные силовые факторы, вызывающие эквивалентный эффект. Приведены результаты численного расчета прикладных задач, подтверждающие полученные практические результаты и проясняющие влияние различных параметров упругих элементов на вид упругой характеристики.
Материалы с памятью формы относятся к классу интеллектуальных материалов, т.е. материалов которые могут контролируемым образом изменять свои свойства в ответ на изменения окружающей среды. Подобная особенность связана с наличием таких свойств как сверхупругость и эффект памяти формы. В настоящее время коммутационные и исполнительные устройства из сплава с эффектом памяти формы широко применяются в различных областях науки и техники [1, 3]. Среди конструктивных исполнений упругих элементов из сплава с эффектом памяти формы широкое распространение получили элементы, выполненные в форме нитей, стержней, плоских и винтовых цилиндрических пружин.
Процесс деформирования упругих элементов, выполненных из сплава с эффектом памяти формы, сложен и зависит от параметров материала, текущих параметров нагружения, а также от истории нагружения, при этом на процесс деформирования влияют не только история изменения температуры и усилий, но и первоначальное фазовое состояние материала. Анализируя комплекс проблем, связанный с разработкой и применением упругих стержневых элементов из сплавов с эффектом памяти формы следует признать, что задача создания эффективных практических методов их расчета до конца не решена. Несмотря на целый ряд разработанных математических моделей для решения задач расчета и синтеза элементов из сплавов с памятью формы применяются эмпирические подходы [4, 5]. При этом, если вопросы изотермического деформирования упругих элементов достаточно проработаны, то деформирование в условиях изменяющейся температурной нагрузки исследовано недостаточно. Вопросы проектирования упругих элементов при больших перемещениях, а также краевые задачи рассматриваются лишь в малом числе узкоспециализированных работ. Разработка методики расчета и проектирования упругих элементов из сплава с эффектом памяти формы является актуальной задачей.
Для описания особенностей сплава с эффектом памяти формы предложена модель материала на основе диаграммы фазовых переходов [3, 4], в которой внутренними переменными являются объемные доли аустенита ${{\xi }_{A}}$, ориентированного мартенсита ${{\xi }_{S}}$ и неориентированного мартенсита ${{\xi }_{M}}$, связанные между собой соотношением
Для определения фазовых долей строятся итерационные соотношения, описывающие заданный фазовый переход в соответствии с экспериментально полученной диаграммой фазовых переходов. Итерационное соотношение для эффекта мартенситной неупругости имеет вид
(1)
$\left\{ \begin{gathered} {\text{If}}\;\;{{\sigma }_{{k + 1}}} > {{\sigma }_{k}} \hfill \\ {\text{ }}\xi _{S}^{{(k + 1)}} = \left\{ \begin{gathered} 0{\text{ , if }}{{\sigma }_{{k + 1}}} \leqslant \sigma _{{start}}^{{MS}} \hfill \\ \left\{ \begin{gathered} {{F}_{{MS}}}({{\sigma }_{{k + 1}}}){\text{, if }}{{F}_{{MS}}}({{\sigma }_{{k + 1}}}) > \xi _{S}^{{(k)}}\;\;{\text{and}}\;\;\sigma _{{start}}^{{MS}} < {{\sigma }_{{k + 1}}} < \sigma _{{finish}}^{{MS}} \hfill \\ \xi _{S}^{{(k)}}{\text{, if }}{{F}_{{MS}}}({{\sigma }_{{k + 1}}}) \leqslant \xi _{S}^{{(k)}}\;\;{\text{and}}\;\;\sigma _{{start}}^{{MS}} < {{\sigma }_{{k + 1}}} < \sigma _{{finish}}^{{MS}} \hfill \\ \end{gathered} \right. \hfill \\ 1{\text{ , if }}{{\sigma }_{{k + 1}}} \geqslant \sigma _{{finish}}^{{MS}} \hfill \\ \end{gathered} \right. \hfill \\ {\text{If}}\;\;{{\sigma }_{{k + 1}}} \leqslant {{\sigma }_{k}} \hfill \\ {\text{ }}\xi _{S}^{{(k + 1)}} = \xi _{S}^{{(k)}}, \hfill \\ \end{gathered} \right.$Для описания закона изменения долей (функция ${{F}_{{MS}}}(\sigma )$) могут быть использованы как аналитические функции (например, тригонометрические, экспоненциальные или линейные), так и дискретный набор экспериментальных данных. В случае использования, тригонометрической аппроксимации, функция ${{F}_{{MS}}}(\sigma )$ имеет вид
Итерационное соотношение для обратного мартенситного превращения принимается в виде
Функция ${{F}_{{SA}}}(T,\sigma )$ для тригонометрической аппроксимации имеет вид
В статье рассматривается кручение стержней с круглым поперечным сечением, что позволяет записать физические соотношения для упрощенного плоского напряженного состояния в виде [5]
(2)
$\left\{ \begin{gathered} \varepsilon = \sigma {{E}^{{ - 1}}} + {{\varepsilon }_{L}}{{\xi }_{S}}\sigma {{({{\sigma }^{2}} + 3{{\tau }^{2}})}^{{ - 0.5}}} \hfill \\ \gamma = \tau {{G}^{{ - 1}}} + 3{{\varepsilon }_{L}}{{\xi }_{S}}\tau {{({{\sigma }^{2}} + 3{{\tau }^{2}})}^{{ - 0.5}}} \hfill \\ \end{gathered} \right.,$Для учета эффекта памяти формы при плоском изгибе вводится дополнительный момент памяти формы изгиба. Учитывая гипотезу плоских сечений соотношения для момента памяти формы изгиба имеют вид
Для учета эффекта памяти формы при кручении стержня круглого поперечного сечения вводится дополнительный момент памяти формы кручения. Учитывая гипотезу о неискривляемости радиуса соотношения для момента памяти формы кручения имеют вид
Дополнительные внутренние силовые факторы позволяют учитывать эффект памяти формы в известных соотношениях механики стержней. Для расчета нитей можно использовать конечно-элементное соотношение
Для построения упругой характеристики плоской пружины при малых перемещениях справедливо соотношение
Для плоской пружины при больших перемещениях справедлива система нелинейных дифференциальных уравнений [6]
При решении геометрически линейных задач плоского изгиба для построения упругой характеристики пружины из сплава с эффектом памяти формы справедливы соотношения
Для построения упругой характеристики винтовой плоской пружины при малых перемещениях справедливо конечно-элементное соотношение
При больших перемещениях для расчета винтовой цилиндрической пружины справедливы соотношения
Для учета эффекта изменения механических свойств материала, вызванного изменением долей мартенситной и аустенитной фаз, справедливо соотношение
Алгоритм построения зависимостей для дополнительных внутренних силовых факторов рассмотрим на примере задачи плоского изгиба с растяжением. Для построения зависимостей момента памяти формы и усилия памяти формы от изгибающего момента и нормального усилия, сечение стержня разбивается на некоторое количество слоев (рис. 1), в пределах каждого из которых напряжения, деформации и объемные фазовые доли считаются постоянными.
Для прямоугольного поперечного сечения параметры слоя определяются соотношениями
(3)
${{\varepsilon }_{i}} = \varepsilon ({{y}_{i}}) = \Delta \unicode{230} {{y}_{i}} + {{\varepsilon }_{0}} = ({{\varepsilon }_{1}} - {{\varepsilon }_{2}}){{y}_{i}}{{h}^{{ - 1}}} + 0.5({{\varepsilon }_{1}} + {{\varepsilon }_{2}}).$Соотношения для дополнительных внутренних силовых факторов после дискретизации принимают вид
(4)
$\begin{gathered} {{N}^{{SMA}}} = E{{\varepsilon }_{L}}\sum\limits_{i = 1}^{2n} {{{b}_{i}}{{\xi }_{S}}({{y}_{i}})\Delta y} \\ M_{b}^{{SMA}} = E{{\varepsilon }_{L}}\sum\limits_{i = 1}^{2n} {{{b}_{i}}{{\xi }_{S}}({{y}_{i}}){\text{sign}}({{\varepsilon }_{i}}){{y}_{i}} \cdot \Delta h} \\ N = 0.5({{\varepsilon }_{1}} + {{\varepsilon }_{2}})EA - {{N}^{{SMA}}} \\ {{M}_{b}} = {{h}^{{ - 1}}}({{\varepsilon }_{1}} - {{\varepsilon }_{2}})E{{I}_{b}} - M_{b}^{{SMA}}. \\ \end{gathered} $Для получения одной точки трехмерных зависимостей для дополнительных внутренних силовых факторов необходимо выполнить последовательность действий: задаемся значениями деформаций ${{\varepsilon }_{1}}$ и ${{\varepsilon }_{2}}$ на верхнем и нижнем слое сечения соответственно; по соотношению (3) определяем деформации в каждом слое; с помощью соотношений (1), (2) определяем величины фазовых долей ориентированного мартенсита в каждом слое ${{\xi }_{S}}({{y}_{i}})$; с помощью соотношений (4) определяем внутренние силовые факторы и дополнительные внутренние силовые факторы, соответствующие заданным деформациям ${{\varepsilon }_{1}}$, ${{\varepsilon }_{2}}$.
На рис. 2 представлены зависимости для дополнительных внутренних силовых факторов, полученные при: ${{E}_{M}} = 30\,000$ МПа, ${{\varepsilon }_{L}} = 0.016$, ${{\sigma }_{S}} = 1$ МПа, ${{\sigma }_{f}} = 100$ МПа, $B = 2$ мм, $H = 2$ мм.
Представленные в статье соотношения и алгоритмы легли в основу авторской программы, для расчета и проектирования упругих элементов из сплава с эффектом памяти формы.
На рис. 3 представлена упругая характеристика консольно-защемленной по краям плоской пружины нагруженной посередине сосредоточенной силой. Рассмотрено изотермическое нагружение в зоне стабильности мартенсита. Использованы геометрические и физико-механические параметры: $L = 30$ мм, $B = 2$ мм, $H = 2$ мм, ${{\varepsilon }_{L}} = 0.016$, ${{\sigma }_{S}} = 1$ МПа, ${{\sigma }_{f}} = 100$ МПа, ${{E}_{M}} = 30\,000$ МПа, ${{E}_{A}} = 73200$ МПа, ${{C}_{A}} = 10$ = 10 МПа/С°, ${{A}_{S}} = 0$°С, ${{A}_{f}} = 20$°С, ${{M}_{S}} = - 10$°С, ${{M}_{f}} = - 20$°С.
На рис. 4 представлены результаты расчета при нагреве плоской пружины при постоянной нагрузке F = 30 Н.
На рис. 5 представлены семейства упругих характеристик винтовой цилиндрической пружины, нагруженной растягивающей силой, при изотермическом нагружении в зоне стабильности мартенсита. Параметры пружины: ${{H}_{0}} = 12$ мм, $i = 10$, параметры материала такие же, как для плоской пружины.
На рис. 6 представлены результаты расчета при нагреве винтовой цилиндрической пружины при постоянной нагрузке F = 0.16 Н.
Получаемые с помощью разработанной программы результаты хорошо согласуются с теоретическими и экспериментальными результатами, представленными в литературе [7–9].
Представленные в статье математическая модель и алгоритм, реализованные в виде прикладной программы, показали свою эффективность и могут быть рекомендованы для расчета и проектирования коммутационных и исполнительных устройств с элементами в форме стержней, плоских и винтовых цилиндрических пружин, изготовленных из материала с эффектом памяти формы.
Конфликт интересов. Авторы заявляют, что у них нет конфликта интересов.
Список литературы
Jani J.M., Leary M., Subic A., Gibson M.A. A review of shape memory alloy research, applications and opportunities. Materials & Design. 2014. V. 56. P. 1078.
Narahari P. S., Wanhill R.J.H. Aerospace Materials and Material Technologies. Volume 1: Aerospace Materials. Springer Science + Business Media, Singapore, 2017. P. 594.
Khandelwal A., Buravalla V. Models for Shape Memory Alloy Behavior: An overview of modeling approaches. International Journal of Structural Changes in Solids . Mechanics and Applications. 2009. V. 1. P. 111.
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// J. of Intelligent Material Systems and Structures. 2013. V. 25. P. 1905.
Azadi B., Rajapakse R.K.N.D., Maijer D.M. Multi-dimensional constitutive modeling of SMA during unstable pseudoelastic behavior // International Journal of Solids and Structures. 2007. 44 (20). P. 6473.
Гаврюшин С.С., Барышникова О.О., Борискин О.Ф. Численный анализ элементов конструкций машин и приборов. Москва: Изд-во МГТУ им. Н.Э. Баумана, 2014. С. 480.
Eshghinejad A., Elahinia M. Exact solution for bending of shapememory alloy beams // Mech. Adv. Mat Struct. 2015. V. 22. (10). P. 829.
Heidari B., Kadkhodaei M., Barati M., Karimzadeh F. Fabrication and modeling of shape memory alloy springs. Smart Materials and Structures. 2016. V. 25 (12). P. 1.
Cheng S.S., Kim Y., Desai J.P. Modeling and characterization of shape memory alloy springs with water cooling strategy in a neurosurgical robot // Journal of Intelligent Material Systems and Structures. 2017. V. 28 (16). P. 2167.
Шишкин С.В., Махутов Н.А. Расчет и проектирование силовых конструкций на сплавах с эффектом памяти формы. М. Ижевск: Регулярная и хаотическая динамика, 2007. С. 412.
Жечков Д.П., Гаврюшин С.С. Расчет стержневых элементов, выполненных из сплавов с памятью формы, при больших перемещениях. Ученые записки Казанского университета. Серия: Физико-математические науки. 2010. Т. 152. № 4. С. 135.
Гаврюшин С.С., Ганыш С.М. Численное моделирование винтовой цилиндрической пружины из сплава с эффектом памяти формы // Сборник: Материалы XIX международного симпозиума “Динамические и технологические проблемы механики конструкций и сплошных сред” имени А.Г. Горшкова Московский авиационный институт. 2013. С. 62.
Гаврюшин С.С., Ганыш С.М. Простейшая математическая модель пространственного стержня, выполненного из сплава с эффектом памяти формы//Сборник: Материалы ХХ международного симпозиума “динамические и технологические проблемы механики конструкций и сплошных сред” имени А.Г. Горшкова Московский авиационный институт. 2014. С. 48.
Ганыш С.М., Гаврюшин С.С. Простейшая математическая модель пространственного стержня, выполненного из сплава с эффектом памяти формы // Инженерный вестник. 2014. № 10. С. 3.
Ганыш С.М., Гаврюшин С.С. Расчет и проектирование упругих элементов из сплавов с эффектом памяти формы // Сборник: Механика и математическое моделирование в технике. II Всероссийская научно-техническая конференция, посвященная юбилеям основателей кафедры “Прикладная механика” МГТУ им. Н.Э. Баумана. Сборник трудов. 2017. С. 33.
Мовчан А.А., Мовчан И.А., Сильченко Л.Г. Микромеханическая модель нелинейного деформирования сплавов с памятью формы при фазовых и структурных превращениях // Известия Российской академии наук. Механика твердого тела. 2010. № 3. С. 118.
Мишустин И.В., Мовчан А.А. Аналог теории пластического течения для описания деформации мартенситной неупругости в сплавах с памятью формы // Известия Российской академии наук. Механика твердого тела. 2015. № 2. С. 78.
Дополнительные материалы отсутствуют.
Инструменты
Проблемы машиностроения и надежности машин