Mesoscopic models for definition of the large-scale elastic properties of the soft magnetic elastomers

Full Text

Введение. Мягкие магнитные эластомеры (ММЭ) выделяются среди других композитных материалов своей способностью к сильному изменению как формы, так и упругих свойств в ответ на приложение магнитного поля. Первопричиной возникновения всех этих особенностей являются микрочастицы магнитомягкого ферромагнетика, которыми наполнена высокоэластичная матрица ММЭ. Силовое взаимодействие частиц, вызванное их намагничиванием, заставляет ММЭ проявлять интересные эффекты - магнитоиндуцированные деформацию и упрочнение, эффект магнитной памяти формы, которые имеют большой практический потенциал. На сугубо качественном уровне картина проста. Приложенное внешнее поле намагничивает частицы и между их индуцированными магнитными моментами возникают пондермоторные силы, которые стремятся организовать частицы в структуру, отвечающую минимуму магнитостатической энергии. Однако в ММЭ этой тенденции противодействуют упругие силы (высокоэластичность), возрастающие по мере удаления частиц от их исходных положений внутри сетки. Чем выше степень сшивки полимерной матрицы, тем выше ее модуль упругости, а значит, тем сильнее она может ограничивать смещения частиц. Однако, как показывают теоретические исследования [1-3], даже небольшие изменения в пространственном распределении частиц внутри эластомера могут приводить к качественно различному поведению и свойствам моделируемого образца композита. Таким образом, критически важным становится понимание магнитоупругого взаимодействия частиц в мягком магнитном эластомере на мезоструктурном уровне как основы эффектов, проявляемых ММЭ. 1. Мезоскопический элемент. Ввиду важности межчастичных взаимодействий для моделирования и предсказания свойств ММЭ необходимо располагать как можно более точным представлением о магнитных и упругих силах, возникающих в композите на мезоскопическом уровне. Для изучения этих вопросов в наших работах [4-6] использовалась модельная система, представляющая собой пару сферических намагничивающихся частиц, помещенных в эластомерную матрицу (см. рис. 1). Такой мезоскопический образец является наименьшим представительным элементом ММЭ, позволяющим изучать его магнитомеханику. 123 Б и л л е р А. М., С т о л б о в О. В. В настоящей работе мы модифицировали нашу модель следующим образом. Размер эластомерного цилиндра, внутрь которого заключены две нелинейно намагничивающихся частицы, подобран так, что частицы занимают 30% геометрического объема образца модельного ММЭ. Именно такова объемная концентрация частиц наполнителя в композитах, представляющих наибольший интерес для приложений. Частицы размещены вдоль оси цилиндра; они имеют сферическую форму и одинаковый радиус Рис. 1. Схематическое изобра- a. Внешнее магнитное поле с напряженножение мезоскопического элемента стью H 0 приложено вдоль оси образца. Под [Figure 1. The mesoscopic element влиянием поля частицы намагничиваются и schematic representation] вступают в магнитостатическое (пондеромоторное) взаимодействие. Магнитная энергия пары Umag в такой конфигурации зависит от расстояния между их центрами l = |l| и величины приложенного поля H0 = |H 0 |. Поскольку энергия пары достигает минимума при плотном контакте намагничивающихся частиц, действующие между ними силы - это притяжение, стремящееся передвинуть частицы навстречу друг другу. В то же время смещение частиц вызывает деформацию эластомерной матрицы, в которую они встроены на условиях полной адгезии. Накапливаемая в матрице упругая энергия Uel зависит от начального l0 и текущего l расстояний между центрами частиц; она минимальна в отсутствие деформаций матрицы, то есть при l = l0 . Таким образом, под действием внешнего магнитного поля в представительном элементе эластомера возникают два типа сил, конкурирующих между собой: магнитные, стремящиеся сблизить частицы, и упругие, противодействующие любой деформации матрицы. Общая энергия U рассматриваемого цилиндра с двумя частицами, включающая магнитную Umag и упругую Uel составляющие, имеет вид U (l, l0 , H0 ) = Umag (l, H0 ) + Uel (l, l0 ). В целом алгоритм получения общей энергии мезоскопического элемента в настоящей работе аналогичен тому, что описан в работе [6]. Приведем его коротко, сосредоточив внимание на выполненной модификации модели и на характерных свойствах рассматриваемой системы. Для магнитной части Umag общей энергии представительного элемента использовались данные численного решения задачи магнитостатики для двух сферических частиц, которое было получено в работе [5]. В этом решении учитывались как неоднородность, так и нелинейность намагничивания частиц. Связь между намагничиванием M и напряженностью поля внутри сферы H определялась эмпирическим законом Фрелиха-Кеннелли [7]: M (H) = 124 χ 0 Ms H = Ms + χ 0 H χ0 H при χ0 H < Ms , Ms (1 - Ms /χ0 H) при χ0 H > Ms . Мезоскопические модели для определения упругих характеристик. . . Здесь χ0 - начальная магнитная восприимчивость частиц, а Ms - намагниченность насыщения. Для частиц карбонильного железа, чаще всего используемых в качестве наполнителя для ММЭ, характерные значения этих величин составляют 104 и 1500 кА/м соответственно. В дальнейшем изложении мы будем использовать безразмерные единицы: напряженность магнитного поля h0 = H0 /Ms , расстояние между центрами частиц q = l/a и магнитную ˜mag (q, h0 ) = Umag (l, H0 )/(µ0 Ms2 a3 ). энергию U Упругая составляющая Uel энергии представительного элемента, накапливаемая в эластомерной матрице при перемещении частиц, была получена аналогично тому, как это было сделано в работе [6], но расчеты выполнялись для цилиндрического образца эластомера с указанными выше пропорциями. Задача нелинейной теории упругости формулировалась в осесимметричной постановке и решалась методом конечных элементов. Поскольку эластомер может подвергаться значительным деформациям, его определяющее соотношение было выбрано в форме закона Муни-Ривлина для несжимаемой среды [8]: ˜ (C) = c1 [(I1 (C) - 3) + c˜2 (I2 (C) - 3)] , W (C) = c1 W (1) где I1 (C) и I2 (C) - инварианты меры деформации Коши-Грина C, c1 и c2 - числовые константы модели. Здесь также введен обезразмеренный потенциал ˜ = W/c1 и параметр c˜2 = c2 /c1 . W В расчетах расстояние l между частицами изменялось дискретно от начального l0 до l = 2a, соответствующего плотному контакту частиц. На каждом шаге методом конечных элементов определялось поле перемещений, затем вычислялась объемная плотность энергии. В результате был сформирован дискретный массив данных, представляющий зависимость упругой энергии Uel от начального l0 и текущего l расстояний между центрами частиц. Для построения и верификации интерполяционной формулы, описывающей зависимость Uel от l, использовалась эвристическая схема, родственная той, что была предложена в работе [6]. Наша конструкция представляет собой систему цилиндрических стержней, которая по общей длине и по расстоянию между частицами идентична пропорциям расчетного образца. Нумерация стержней в этой схеме начинается с 2; номер 1 зарезервирован для обозначения в дальнейшем относительного изменения расстояния между частицами. Рис. 2. Эвристическая схема для интерполяции упругой энергии [Figure 2. The heuristic scheme for the elastic energy interpolation] 125 Б и л л е р А. М., С т о л б о в О. В. Стержни подчиняются тому же соотношению Муни-Ривлина (1), что и эластомерная матрица в численном расчете. Коэффициентами интерполяции (2) служат радиусы стержней rk и длина l0 внутреннего стержня 2. Формула для Uel при каждом конкретном q0 имеет вид 4 (k) ˜ (Ck )r2 l . (2 - δk4 )W k 0 Uel (l) = πc1 (2) k=2 Интерполяция строилась с учетом того, что при малых деформациях образца, вызванных перемещением частиц внутри него или внешней нагрузкой, разложение энергии в ряд совпадает с таким же разложением энергии стержневой модели. Этот факт позволяет найти часть интерполяционных коэффициентов по коэффициентам разложения, а оставшиеся подобрать из соображений наилучшего соответствия между энергией стержневой модели и расчетного образца в области больших деформаций. На следующем шаге решения коэффициенты формулы (2) интерполируются посредством различных степенных функций. В результате мы получили непрерывную функцию 4 ˜ (Ck )r2 (l0 )l(k) (l0 ), (2 - δk4 )W k 0 Uel (l, l0 ) = πc1 (3) k=2 2 r2 = 1.003 - 0.2255/(1 + e(0.8989l0 -6.311) ), r3 = 1.103 - 0.01318l02 - 3.301 · 10-6 l08 , r4 = 0.1637 + 0.02213(l0 - 3.045)2 , (2) l0 = 1.464 - 4.232l0 + 2.314l02 - 0.2781l03 , (3) (4) l0 = (l0 - l0 - 2a)/2, (4) l0 = 6.4a, описывающую упругую энергию, накапливаемую в эластомерном цилиндре с погрешностью не более 5%. Суммирование интерполяционных формул для магнитной и упругой составляющих энергии дает аналитическую функцию ˜mag (q, h0 ) + U ˜el (q, q0 ), U (q, q0 , h0 )/(c1 a3 ) = β U являющуюся хорошим приближением точной магнитоупругой энергии мезоскопического элемента; в ней относительная величина магнитного и упругого вкладов определяется параметром β = µ0 Ms2 /c1 . Анализ функции (3) показывает, что в некотором диапазоне полей система переходит в режим бистабильности: в ней одновременно существуют два положения устойчивого равновесия. Следствием бистабильности является возникновение гистерезиса межчастичного расстояния при изменении величины приложенного магнитного поля (см. рис. 3). В системах, где частицы намагничиваются с насыщением, возможность появления гистерезиса при заданном начальном расстоянии между частицами q0 = l0 /a контролируется параметром β, который имеет смысл эффективной податливости композита. 126 Мезоскопические модели для определения упругих характеристик. . . Отметим, что гистерезисное поведение возникает вне зависимости от того, находятся ли частицы в неограниченной эластомерной среде [5] или в образце конечного размера [6]. 2. Мезоскопический элемент под действием внешних сил. В случае отсутствия сторонних механических сил интерполяционная формула (3) для энергии Uel цилиндрического образца позволяет найти зависимость изменения его общего («внешнего») размера u4 от изменения параметра u1 - расстояния между центрами частиц по оси цилиндра. В присутствии сторонней механической нагрузки (растяжение или сжатие, приложенное к торцам цилиндра) относительные изменения расстояния между частицами λ1 и общей длины образца λ4 становятся независимыми переменными. Соответственно, энергетическая функция превращается в U (λ1 , λ4 , q0 , h0 ), что дает возможность получить упругие характеристики модельного образца эластомера. Введем систему координат так, чтобы ось Oz была направлена вдоль вектора l, соединяющего центры частиц в образце, упругий отклик которого мы описываем системой стержней (см. рис. 2). Под действием внешнего магнитного поля, направленного по Oz, образец деформируется так, что относительные изменения межчастичного расстояния (воображаемый стержень 1) и стержня 4 принимают значения ε∗1 = λ∗1 - 1 и ε∗4 = λ∗4 - 1 соответственно. В отсутствие внешней нагрузки состояние системы, задаваемое параметрами (ε∗1 , ε∗4 ), является равновесным. В зависимости от величины приложенного поля равновесному состоянию системы могут соответствовать как небольшое сближение частиц, так и их фактический коллапс: q 2. Рассмотрим, к примеру, мезоскопический элемент с q0 = 3.5, магнитомеханический гистерезис которого в числе других представлен на рис. 3. Такой гистерезис является «латентным» [5] в том смысле, что, хотя парный кластер является одним из двух возможных устойчивых положений системы, перейти в это состояние только под действием магнитного поля она не может. Однако коллапс частиц можно реализовать при добавлении внешней механической нагрузки, и это, конечно, отразится Рис. 3. Магнитомеханический гистерезис в системе с β = 110 и начальным межцентровым расстоянием q0 , равным 2.5 (линия 1), 3 (линия 2), 3.5 (линия 3), 4 (линия 4) [Figure 3. Magnetomechanical hysteresis in the system with β = 110 and initial interparticle distance q0 equals 2.5 (line 1), 3 (line 2), 3.5 (line 3), and 4 (line 4)] 127 Б и л л е р А. М., С т о л б о в О. В. на общей длине образца. Приложим к торцам образца силы величиной f , сжимающие или растягивающие его. Изменяя величину приложенной силы и вычисляя равновесную относительную деформацию мезоскопического элемента ε4 , то есть решая задачу оптимизации U (ε1 , ε4 ) - f u4 → min, получаем диаграмму «напряжение-деформация», которая представлена на рис. 4, a. Рис. 4. Диаграммы «напряжение-деформация» мезоскопических элементов с относительной податливостью β = 110 и начальными расстояниями между центрами частиц: a) q0 = 3.5 в поле с напряженностью h0 = 0.32; b) q0 = 3.25 в полях c напряженностью h0 = 0.16 (линия 1), 0.24 (линия 2), 0.32 (линия 3) [Figure 4. Stress-strain diagrams for mesoscopic elements with relative compliance β = 110 and initial interparticle distance: a) q0 = 3.5 under applied field intensity h0 = 0.32 and b) q0 = 3.25 under applied field intensity h0 = 0.16 (line 1), 0.24 (line 2), and 0.32 (line 3)] Как видно из рис. 4, a, если мезоскопический элемент намагничен, то он уже немного деформирован даже в отсутствие внешней механической силы (см. точку 1). Это сжатие соответствует известному магнитострикционному эффекту, присущему ММЭ. Пусть к образцу приложена сжимающая нагрузка, стремящаяся уменьшить его длину. С ростом силы f сжатие образца будет сначала идти плавно, но по достижении силой некоторого конечного значения произойдет коллапс частиц, и деформация изменится скачком (см. точку 2). На диаграмме рис. 4, a это соответствует перескоку представляющей точки с нижней ветки кривой на верхнюю и дальнейшему движению вдоль последней. Примечательно, что при снятии нагрузки система не возвращается в начальное положение (см. точку 3), а продолжает оставаться в сжатом состоянии («магнитная скрепка»). Разрушить эту структуру можно только с помощью усилия противоположного направления - растягивающего (см. точку 4). В результате графическое представление цикла «нагрузка-разгрузка» отображается на рис. 4, a в виде петли, то есть на определенном участке функция обладает двузначностью. Приведенные рассуждения показывают, что если в мезоскопическом элементе существует магнитомеханический гистерезис, то на кривой «напряжение-деформация» петля будет появляться в магнитном поле любой напряженности. Этот вывод иллюстрирует серия кривых на рис. 4, b, представляющая деформационные циклы для мезоскопического элемента с q0 = 3.25 и β = 110. Как видно, положение петли зависит от величины приложенного поля. Так, для элемента, находящегося в поле с напряженностью h0 = 0.16, 128 Мезоскопические модели для определения упругих характеристик. . . она целиком расположена в области сжимающей нагрузки, а для того же элемента в поле h0 = 0.32 петля имеет место только при растягивающих усилиях. Из этого следует, что, регулируя величину приложенного магнитного поля, можно моделировать свойства мезоскопического элемента: устанавливать величину сил, при которых будет происходить коллапс или расцепление частиц, то есть увеличиваться или уменьшаться способность элемента сопротивляться деформации. Необходимо отметить, что эффективная податливость мезоскопического элемента β = 110 (ед. СИ) соответствует достаточно жесткой полимерной матрице. Для реальных ММЭ с модулем Юнга порядка 10-30 кПа (см. например [9]) величина параметра β достигает нескольких тысяч. В этом случае в рассматриваемых системах с начальным межчастичным расстоянием q0 от 2 до 4 магнитомеханический гистерезис существует всегда, и коллапс частиц происходит в полях меньших h0 = 0.1 (см. рис. 5). Рис. 5. Магнитомеханический гистерезис в системе с β = 2500 и начальным межцентровым расстоянием q0 , равным 2.5 (линия 1), 3 (линия 2), 3.5 (линия 3), 4 (линия 4) [Figure 5. Magnetomechanical hysteresis in the system with β = 2500 and initial interparticle distance q0 equals 2.5 (line 1), 3 (line 2), 3.5 (line 3), and 4 (line 4)] Разрушить «магнитную скрепку», образовавшуюся в элементе с высокой податливостью, непросто. Для этого требуется растянуть образец в несколько раз. Для примера в таблице приведены характерные значения параметра f /(c1 S), необходимые для образования и разрушения «магнитной скрепки» в мезоскопическом элементе с относительной податливостью β = 2500, а также значения деформации элемента до (εs ) и после (εe ) по обе стороны от скачкообразного перехода. Данные в таблице соответствуют точкам, отмеченным на рис. 4, a: 2 - коллапс частиц и 4 - разрушение парного кластера. Характерные точки диаграммы «напряжение-деформация» для мезоскопического образца с β = 2500 и расстоянием между центрами частиц q0 = 3.5 [Typical points of the stress-strain diagram for the mesoscopic element with β = 2500 and interparticle distance q0 = 3.5] h0 f /(c1 S) = 0 ε 0.02 0.04 0.06 -0.021 -0.080 -0.080 Point 2, see Fig. 4, a f /(c1 S) εs εe -1.58 3.40 22.95 -0.109 0.416 5.373 -0.140 0.285 3.981 Point 4, see Fig. 4, a f /(c1 S) εs εe 22.95 60.90 95.40 5.494 14.727 21.946 3.981 11.288 17.933 129 Б и л л е р А. М., С т о л б о в О. В. 3. Упругие свойства мезоскопического элемента. Для определения механических характеристик мезоскопического элемента разложим его энергию в ряд в окрестности положения равновесия (ε∗1 , ε∗4 ): U (ε∗1 , ε∗4 ) = 1 ∂2U 2 ∂ε21 ε4 = ε∗4 , ε1 = ε∗1 ε21 + 1 ∂2U 2 ∂ε24 ε4 = ε∗4 , ε1 = ε∗1 ε24 + ∂2U ∂ε1 ∂ε4 ε4 = ε∗4 , ε1 ε4 . ε1 = ε∗1 По понятным причинам производные первого порядка здесь обращаются в нуль. Приложенная малая нагрузка вызывает дополнительную малую деформацию образца. Возникающие при этом силы должны уравновешивать как магнитные силы взаимодействия частиц, так и внешнюю механическую нагрузку. Следовательно, энергия образца в сумме с работой этих сил должна стремиться к минимуму: U (ε1 , ε4 ) - f u4 → min. (4) (4) Запишем работу внешних сил в виде f u4 = f ε4 l0 , используя соотношение (4) между деформацией и перемещением ε4 = u4 /l0 . Тогда задача оптимизации (4) запишется в виде 1 1 (4) (4) U - f ε4 l0 = A11 ε21 + A22 ε24 + A12 ε1 ε4 - f ε4 l0 → min, 2 2 где Aij - вторые частные производные энергии. В этих обозначениях экстремум рассматриваемой функции описывается решениями системы уравнений ∂U = A11 ε1 + A12 ε4 = 0, ∂ε1 ∂U (4) = A22 ε4 + A12 ε1 = f l0 . ∂ε4 (5) Введем обозначение f /ε4 = K, которое определяет способность рассматриваемой модели сопротивляться деформации. Из системы (5) видно, что f A2 1 = A22 - 12 (4) . ε4 A11 l 0 ˜ как отноОпределим относительную жесткость мезоскопического образца K шение найденной K системы в магнитном поле к начальной жесткости K0 в его отсутствие. Равновесное положение частиц (ε∗1 , ε∗4 ) в элементе изменяется в зависимости от приложенного магнитного поля гистерезисным образом. Аналогично изменяется и жесткость рассматриваемой системы. На рис. 3 было представлено изменение расстояния между частицами в мезоскопических элементах с различными q0 под действием внешнего магнитного поля. Соответствующее ˜ приведено на рис. 6. изменение относительной жесткости элементов K С ростом внешнего поля частицы плавно сближаются, а модуль системы уменьшается. Затем при определенной величине поля в системе появляется K= 130 Мезоскопические модели для определения упругих характеристик. . . вторая жесткость, что соответствует возникновению в системе второго положения равновесия. С дальнейшим увеличением поля происходит коллапс частиц, который влечет за собой резкий рост жесткости мезоскопического элемента. Чем больше начальное расстояние между частицами q0 , тем сильнее эффект упрочнения двухчастичного образца в магнитном поле. 4. Процедура гомогенизации. Двухчастичная модель позволяет разобраться в основных закономерностях перестройки внутренней структуры ММЭ и следующим за ней изменением свойств композита. Однако для предсказания свойств макроскопических образцов необходимо перенести на них результаты мезоскопического моделирования. Для этого используются процедуры гомогенизации, некоторые варианты которых описаны в работах [10, 11]. Поскольку наша модель учитывает лишь одну конфигурацию частиц в эластомере во внешнем поле, хотя и вносящую серьезный вклад, мы ограничимся здесь простой процедурой осреднения. Представим некоторый макроскопический образец ММЭ как систему наших мезоскопических элементов, каждый из которых содержит пару частиц с различным начальным межцентровым расстоянием q0 . Предположим, что к образцу приложена растягивающая сила в направлении вдоль длин цилиндров. Будем считать, что каждый мезоскопический элемент под действием этой силы будет растягиваться независимо. Таким образом, мы не учитываем никакого взаимодействия на границах между отдельными элементами. Для вычисления модуля упругости макроскопического образца или напряжений, возникающих при его растяжении/сжатии, требуется только усреднить соответствующие характеристики по всем мезоскопическим элементам. Подобные модели материалов называют статистическими, основанными на гипотезе Фойгта. Описываемый подход можно сравнить с построением феноменологической модели материала, где элементами схемы выступают наши мезоскопические образцы, соединенные параллельно. Поскольку число частиц в макроскопических образцах ММЭ значительно, то и число элементов в схеме необ- Рис. 6. Изменение относительной жесткости мезоскопических элементов с β = 2500 в зависимости от приложенного магнитного поля при начальном расстоянии между частицами q0 , равным 2.5 (линия 1), 3 (линия 2), 3.5 (линия 3), 4 (линия 4) [Figure 6. Relative stiffness of mesoscopic elements with β = 2500 depend on the apllied magnetic field for the initial interparticle distance q0 equals 2.5 (line 1), 3 (line 2), 3.5 (line 3), and 4 (line 4)] 131 Б и л л е р А. М., С т о л б о в О. В. ходимо устремить к бесконечности и усреднение проводить по непрерывной функции распределения расстояния между центрами частиц. Таким образом, жесткость макроскопического образца с заданным распределением расстояний между частицами вычисляется по формуле q0max K(H0 ) = q0min q0max K(q0 , h0 )f (q0 )dq0 q0min f (q0 )dq0 , где f (q0 ) - функция плотности распределения q0 . Для ММЭ логично предположить нормальное распределение расстояний между частицами наполнителя. Максимальное межчастичное расстояние, которое допускает наш мезоскопический образец, составляет 4.25a. Зададим распределение со средним значением m = 3a и среднеквадратичным отклонением σ = 0.2. В результате получим кривую зависимости жесткости ММЭ от величины приложенного поля h0 , которая представляет собой простое усреднение множества подобных кривых для мезоскопических элементов. Жесткость макроскопического образца ММЭ сначала слегка уменьшается, а затем плавно растет с увеличением поля и в полях порядка ∼ 0.6 в безразмерных единицах достигает насыщения. В обратном направлении при снятии внешнего поля наблюдается гистерезис жесткости. Такое поведение легко объяснить, опираясь на графики изменения модуля для мезоскопических элементов с различным q0 , представленные на рис. 6. При включении поля частицы сближаются и во всех элементах без исключения жесткость уменьшается. По достижении напряженностью внешнего поля значения 0.02 на графике рис. 7 возникают два различных значения жесткости. Это соответствует области гистерезиса в мезоскопических элементах с межчастичным расстоянием около 2.5. С ростом поля первыми в плотный контакт приходят пары с наименьшим q0 . Увеличение жесткости в таких парах не настолько велико, чтобы уравновесить его уменьшение в элементах с большим q0 . В рассматриваемом ММЭ наибольшее число пар расположено на расстоянии q0 = 3 между центрами частиц, поэтому в поле h0 = 0.045, когда эти пары коллапсируют, начинается значительный рост жесткости образца. В поле h0 = 0.06 даже частицы, находящиеся на большом расстоянии друг от друга, схлопываются, рост жесткости ММЭ останавливается и эффект упрочнения достигает насыщения. Приведенное выше описание позволяет без труда разобраться, как вид распределения q0 влияет на получаемую в результате жесткость макроскопи˜ от величины ческого образца ММЭ. На рис. 8 представлены зависимости K приложенного поля для других распределений. Смещение пика распределения в область больших межчастичных расстояний приводит к увеличению максимальной жесткости ММЭ. При этом уменьшение разброса q0 увеличивает разницу между наибольшим и наименьшим значением жесткости образца. Заключение. В настоящей работе показано, как гистерезисное поведение частиц в мезоскопическом образце влияет на его механические характеристики. На качественном уровне эффект магнитоиндуцированного упрочнения ММЭ может быть объяснен образованием плотных кластеров из частиц внутри композита при его намагничивании, что проиллюстрировано в работе с помощью процедуры осреднения. Несмотря на то, что рассматрива132 Мезоскопические модели для определения упругих характеристик. . . Рис. 7. Относительная жесткость образца ММЭ, имеющего нормальное распределение межчастичных расстояний с m = 3 и σ = 0.2 в зависимости от внешнего магнитного поля h0 [Figure 7. Relative stiffness of the SME sample, that has normal distribution with m = 3 and σ = 0.2, depend on the extrenal magnetic field h0 ] Рис. 8. Зависимость относительной жесткости образцов ММЭ от внешнего магнитного поля h0 . Образцы имеют нормальное распределение межчастичных расстояний: a) с различным m и σ = 0.2; вид распределений представлен на части b рисунка; c) с m = 3 и различным σ; вид распределений представлен на части d рисунка [Figure 8. Dependence of relative stiffness of the SME samples on the external magnetic field h0 . Samples have normal distribution of interparticle distances: a) with different m and σ = 0.2; distributions presented on the part b of the plot; c) with m = 3 and different σ; distributions presented on the part d of the plot] 133 Б и л л е р А. М., С т о л б о в О. В. емая модель учитывает только те пары частиц, которые расположены вдоль приложенного магнитного поля, представленные зависимости относительной жесткости ММЭ от внешнего магнитного поля качественно совпадают с кривыми изменения упругих модулей структурированных образцов композита, полученных экспериментально и опубликованных в работах [9, 12]. Конкурирующие интересы. Мы не имеем конкурирующих интересов. Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.

About the authors

Anastasiya M Biller

Institute of Continuous Media Mechanics, Ural Branch of RAS

Email: kam@icmm.ru
Juniour Researcher; Physics and Mechanics of Soft Matter Laboratory 1, Akad. Korolyova st., Perm, 614013, Russian Federation

Oleg V Stolbov

Institute of Continuous Media Mechanics, Ural Branch of RAS

Email: sov@icmm.ru
Cand. Phys. & Math. Sci.; Researcher; Physics and Mechanics of Soft Matter Laboratory 1, Akad. Korolyova st., Perm, 614013, Russian Federation

References

Supplementary files