Доклады Российской академии наук. Физика, технические науки, 2023, T. 510, № 1, стр. 70-76

СТРУКТУРА И МЕХАНИЧЕСКИЕ СВОЙСТВА КОМПОЗИТОВ С ВЫСОКИМИ ДЕМПФИРУЮЩИМИ СВОЙСТВАМИ

Член-корреспондент РАН Е. В. Ломакин 12*, С. А. Лурье 2**, Л. Н. Рабинский 2

1 Московский государственный университет имени М.В. Ломоносова
Москва, Россия

2 Московский авиационный институт (национальный исследовательский университет)
Москва, Россия

* E-mail: evlomakin@yandex.ru
** E-mail: salurie@mail.ru

Поступила в редакцию 09.08.2022
После доработки 09.08.2022
Принята к публикации 08.12.2022

Полный текст (PDF)

Аннотация

Исследуются демпфирующие свойства композитных материалов, армированных упругими сферическими и волокнистыми включениями с вязкоупругими покрытиями. Показывается, что в композитах с морфологией частиц может быть реализован эффект значительного усиления диссипативных потерь, когда эффективные свойства потерь композита могут превышать диссипативные свойства вязкоупругого покрытия более чем в 20 раз. Предлагаются аналитические оценки оптимальных параметров гибридных композитов. Анализируется влияние возможных несовершенств структуры композита на эффективные диссипативные свойства.

Ключевые слова: дискретные и волокнистые композиты, вязкоупругое покрытие, модуль потерь, оптимальные структуры

ВВЕДЕНИЕ

В аэрокосмической промышленности традиционным является значительный интерес к композитным материалам с высоким демпфированием колебаний и, одновременно, с высокой прочностью/жесткостью. Известно несколько “традиционных” концепций оптимизации структуры материалов для улучшения демпфирующих свойств композитов на микромеханическом или макромеханическом уровне [1, 2]. Весьма перспективной является идея создания композитных материалов, включающих упругие включения с тонкими вязкоупругими покрытиями, которая может быть основой для проектирования нового типа композитов, обладающих одновременно высоким демпфированием и высокими эффективными механическими свойствами [310]. При поиске оптимальных свойств в работе [11] исследовалось влияние толщины слоя вязкоупругого покрытия, окружающего сферические включения, внедренные в эпоксидную матрицу. Отмечено, что эффект усиления диссипативных свойств не является монотонным с точки зрения толщин покрытий. Имеет место ярко выраженный экстремум в значениях модуля потерь композита, который показывает, что оптимальное эффективное значение модуля потерь значительно превышает модули потерь отдельных компонентов композита, включая и модуль потерь вязкоупругого покрытия. Аналогичные эффекты обсуждались в [12] и для слоистых композитов.

В работе исследуются механизмы процессов диссипации при деформировании гибридных композитов, имеющих включения с вязкоупругими покрытиями. Фактически впервые установлены закономерности изменения механизмов диссипации с изменением параметров композитов и предложены методы определения оптимальных характеристик. Показано, что оптимальные значения толщин вязкоупругих покрытий в высокой степени локализованы и, как ни странно, являются фактически одинаковыми для дисперсных и волокнистых композитов, если они имеют одни и те же свойства компонент. Предложены аналитические методы определения эффективных модулей упругости и оптимальных характеристик композитов, при которых реализуются аномально высокие эффективные демпфирующие свойства.

МОДЕЛИРОВАНИЕ ЭФФЕКТИВНЫХ СВОЙСТВ

Для описания диссипативных свойств композитов с упругими включениями, имеющими покрытия из вязкоупругого материала, будем использовать метод комплексных модулей упругости в рамках вязкоупругой аналогии. Тогда эффективный модуль упругости определяется как комплексный модуль E = E ' + i E '', где действительная часть E ' является динамическим модулем упругости, а мнимая часть E '' представляет собой модуль потерь композита. Предполагается, что в композите реализовано поле деформаций, изменяющееся по гармоническому закону $\varepsilon \, = \,{{\varepsilon }_{0}}\,{\text{sin}}\omega t$. Соответствующие напряжения при диссипации энергии меняются по закону $\sigma = {{\sigma }_{0}}\sin (\omega t + \delta )$. Мнимая часть модуля упругости позволяет формально оценить энергию диссипации в единице объема Ù = $1{\text{/}}2E{\kern 1pt} '{\kern 1pt} '\varepsilon _{0}^{2}$. Отношение E ''/E ' называется тангенсом потерь (tg$\delta $), который связан с демпфированием, а $\delta $ – угол сдвига по фазе между напряжением и деформацией при синусоидальной нагрузке $E = E{\kern 1pt} '(1 + i\eta ),$ $E{\kern 1pt} '{\kern 1pt} ' = E{\kern 1pt} '\eta ,$ $\eta = {\text{tg}}\delta $. Подобные оценки обычно проводят для угловой частоты $\omega = 1$ рад/с.

Рассмотрим композит со слоистым сферическим или цилиндрическим включением и введем неоднородный представительный фрагмент ${{G}_{1}} \cup {{G}_{2}} \cup {{G}_{M}}$, где подобласти ${{G}_{1}} = \left\{ {r < {{r}_{0}}} \right\}$, G2 = = $\{ {{r}_{0}} < r < {{r}_{1}}\} $ и ${{G}_{M}} = \left\{ {{{r}_{1}} < r < {{r}_{2}}} \right\}$ являются включением, промежуточным слоем вязкоупругого полимера и матрицей соответственно. Для определения эффективных свойств таких композитов воспользуемся самосогласованным методом трех фаз Эшелби [13], показавшим высокую точность почти для предельных концентраций.

Считаем, что составной, слоистый фрагмент композита ${{G}_{1}} \cup {{G}_{2}} \cup {{G}_{M}}$ погружен в бесконечную среду с эффективными свойствами GH. Для замыкания задачи по определению эффективных характеристик привлекается интегральное равенство Эшелби, отражающее условие гомогенизации и указывающее на энергетическую эквивалентность составной среды с однородной средой, имеющей эффективные осредненные свойства:

(1)
$U{\kern 1pt} ' = \int\limits_S {(\sigma _{{ij}}^{M}u_{i}^{0} - \sigma _{{ij}}^{0}u_{i}^{M})} \,dS = 0\quad \left( {i,j = 1,2,3} \right),$
где $U{\kern 1pt} '$ – приращение энергии в элементарной ячейке матричного материала, содержащего составное включение; S – поверхность контакта слоя матрицы ${{G}_{M}}$ с эффективной средой ${{G}_{H}}$; $\sigma _{{ij}}^{M}$, $u_{i}^{M}$ – компоненты тензора напряжений и вектора перемещений на контактной поверхности составного включения (включая матрицу) и эффективной среды, находящиеся из решения контактной задачи для составного включения ${{G}_{1}} \cup {{G}_{2}} \cup {{G}_{M}}$; $\sigma _{{ij}}^{0}$, $u_{i}^{0}$ – компоненты тензора напряжений и вектора перемещений на границе контакта с эффективной средой, связанные с условиями задачи на бесконечности в слое ${{G}_{H}}$.

Общее представление для вектора перемещений ${\mathbf{U}} = \left\{ {{{U}_{i}}} \right\}$, как для слоистых сферических, так и для слоистых цилиндрических включений, имеет один и тот же вид, если использовать соотношение Папковича–Нейбера, записываемое через гармонический потенциал ${{{\mathbf{f}}}_{0}} = \{ f_{i}^{0}\} $ [13]:

${\mathbf{U}}(P) = \frac{{{{{\mathbf{f}}}_{0}}(P)}}{\mu } - \frac{{\nabla \left( {{\mathbf{r}}{{{\mathbf{f}}}_{0}}} \right)}}{{[4\mu (1 - \nu )]}},\quad {{\nabla }^{2}}{{{\mathbf{f}}}_{0}}(P) = 0.$

Для построения решения для составного включения следует выполнить граничные условия на сферических (цилиндрических) поверхностях контакта фаз

(2)
$\left[ {{{U}_{i}}} \right] = \left[ {{{p}_{i}}({\mathbf{U}})} \right] = 0,$
где ${{p}_{i}}({\mathbf{U}})$ – векторы напряжений, определенные на поверхности контакта слоев с вектором нормали n.

В общем случае решение в каждой из подобластей ${{G}_{n}}$ определяется с помощью специального разложения искомого гармонического потенциала ${{{\mathbf{f}}}_{0}}(P)$,$P \in {{G}_{n}}$, построенного с привлечением потенциала, описывающего однородное нагружение (однородные деформации) на бесконечности ${\mathbf{f}}_{0}^{{(H)}}(P) = \left\{ {f_{i}^{{(H)}}} \right\}$:

для сферических включений

$\begin{gathered} {{{\mathbf{f}}}_{0}}(P) = {{A}_{n}}{\mathbf{f}}_{0}^{{(H)}} + {{B}_{n}}{{r}^{{ - 3}}}{\kern 1pt} {\mathbf{f}}_{0}^{{(H)}} + {{C}_{n}}\nabla div({{r}^{{ - 3}}}{\kern 1pt} {\mathbf{f}}_{0}^{{(H)}}) + \\ \, + {{D}_{n}}{{r}^{7}}\nabla div({{r}^{{ - 3}}}{\mathbf{f}}_{0}^{{(H)}}),\quad r = \sqrt {{{x}^{2}} + {{y}^{2}} + {{z}^{2}}} ; \\ \end{gathered} $

для цилиндрических включений

$\begin{gathered} {{{\mathbf{f}}}_{0}}(P) = {{A}_{n}}{\mathbf{f}}_{0}^{{(H)}} + {{B}_{n}}{{r}^{{ - 2}}}{\kern 1pt} {\mathbf{f}}_{0}^{{(H)}} + {{C}_{n}}\nabla div({{r}^{{ - 2}}}{\kern 1pt} {\mathbf{f}}_{0}^{{(H)}}) + \\ \, + {{D}_{n}}{{r}^{6}}\nabla div({{r}^{{ - 2}}}{\mathbf{f}}_{0}^{{(H)}}),\quad r = \sqrt {{{x}^{2}} + {{y}^{2}}} . \\ \end{gathered} $

Здесь ${{A}_{n}},$ ${{B}_{n}},$ ${{C}_{n}},$ ${{D}_{n}}$ – неизвестные постоянные, определяемые из условий контакта (2).

На сферической (и цилиндрической) поверхности контакта векторный потенциал может быть представлен в виде разложения на нормальную ${\kern 1pt} f_{n}^{{(H)}} = ({\mathbf{nf}}_{0}^{{(H)}})$ и касательную ${\kern 1pt} {\mathbf{f}}_{s}^{{(H)}} = {\mathbf{f}}_{0}^{{(H)}} - {\mathbf{n}}({\mathbf{nf}}_{0}^{{(H)}})$ компоненты. Также могут быть представлены вектор перемещений и вектор напряжений в контактных граничных условиях (2). В результате в общем случае в каждой точке поверхности контакта имеем четыре граничных условия.

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

Для определения эффективного модуля объемного сжатия композит со сферическими включениями подвергается действию напряжений, соответствующих гармоническому потенциалу ${\kern 1pt} {\mathbf{f}}_{0}^{{(H)}}$ = = {x, y, z}. Для волокнистого композита, армированного в направлении оси 0x, в случае плоской объемной деформации следует принять ${\kern 1pt} {\mathbf{f}}_{0}^{{(H)}}$ = {0, y, z}. В результате задача упрощается и сводится к определению шести неизвестных постоянных ${{A}_{i}},{{B}_{i}},$ $i = 1{\kern 1pt} - {\kern 1pt} 3$ из двух условий на границе контакта $\sigma _{{rr}}^{n}({{r}_{i}})\, = \,\sigma _{{rr}}^{{n + 1}}({{r}_{i}}),$ $U_{r}^{n}({{r}_{i}}) = U_{r}^{{n + 1}}({{r}_{i}}),$ $(U_{r}^{{}} = {{U}_{n}},{{p}_{n}} = \sigma _{{rr}}^{{}})$. Следуя самосогласованному методу Эшелби, к числу неизвестных добавляется эффективный модуль объемной деформации четвертой фазы ${{K}_{H}}$ и учитывается интегральное равенство Эшелби (1), которое сводится к требованию равенства нулю постоянной ${{B}_{4}} = 0$. В результате модуль упругости объемной деформации всегда находится явно. Отметим, что использование интегрального равенства Эшелби всегда приводит к необходимости выполнения условия ${{B}_{4}} = 0$ [13].

Для девиатора деформаций одному векторному уравнению соответствует два алгебраических уравнения, соответственно для нормальной и касательной составляющих векторного поля, и четыре граничных условия в каждой точке поверхности контакта для 12 постоянных ${{A}_{n}},$ ${{B}_{n}},$ ${{C}_{n}},$ ${{D}_{n}}$ и одной дополнительной постоянной – искомого эффективного модуля сдвига среды. Например, задача определения эффективного модуля ${{\mu }_{H}}$ в поперечной плоскости является полноценной векторной задачей и решается с использованием состояния чистого сдвига ${\mathbf{f}}_{0}^{{(H)}} = \left\{ {x, - y,0} \right\}$. В этом случае контактные условия приводят к системе 12 алгебраических уравнений относительно двенадцати постоянных, которые зависят от неизвестного модуля ${{\mu }_{H}}$. В данной задаче зависимость коэффициента ${{B}_{4}}$ от ${{\mu }_{H}}$ нелинейная, поэтому требуется решать трансцендентное уравнение итерационным способом. В результате для выполнения условия ${{B}_{4}} = 0$ устраивается итерационный процесс, в котором изменяется эффективный мо-дуль ${{\mu }_{H}}$.

Для композита со сферическими включениями проблема определения эффективных свойств заканчивается определением эффективного объемного модуля упругости и эффективного модуля сдвига. После вычисления ${{K}_{H}}$ и ${{\mu }_{H}}$ модуль Юнга и коэффициент Пуассона определяются равенствами: ${{E}_{H}}\, = \,9{{K}_{H}}{{\mu }_{H}}{\text{/}}(3{{K}_{H}}\, + \,{{\mu }_{H}}),$ ${{\nu }_{H}} = (3{{K}_{H}}$ – ‒ 2μH)/2(3KH + μH). Для композитных материалов с цилиндрическими включениями эффективная среда является трансверсально-изотропной. Предлагаемый подход позволяет вычислить “плоский” объемный модуль упругости KH = ${{({{C}_{{2222}}} + {{C}_{{3333}}})} \mathord{\left/ {\vphantom {{({{C}_{{2222}}} + {{C}_{{3333}}})} 2}} \right. \kern-0em} 2}$, определяемый функцией ${\mathbf{f}}_{0}^{{(H)}} = \left\{ {0,y,z} \right\}$, и эффективный модуль сдвига в поперечной плоскости, ортогональной волокнам, $\mu _{H}^{{(1)}} = {{C}_{{1212}}}$, когда композит подвергается деформации чистого сдвига ${\mathbf{f}}_{0}^{{(H)}} = \left\{ {0, - y,z} \right\}$. Дополнительно вычисляется продольный модуль растяжения/сжатия ${{k}_{H}} = {{C}_{{1111}}}$ (по правилу смеси) и поперечный модуль Ламе $\lambda _{H}^{{(2)}} = {{C}_{{1122}}} = {{C}_{{1133}}}$, когда композит подвергается растяжению вдоль оси многослойного цилиндра ${\mathbf{f}}_{0}^{{(H)}} = \left\{ {x,0,0} \right\}$. Аналогично решается проблема вычисления модуля сдвига в продольной плоскости $\mu _{H}^{{(2)}} = {{\mu }_{{13}}} = {{C}_{{1212}}} = {{C}_{{1313}}}$ по условию чистого сдвига вдоль образующей многослойного цилиндра ${\mathbf{f}}_{0}^{{(H)}} = \left\{ {x, - {\kern 1pt} {y \mathord{\left/ {\vphantom {y 2}} \right. \kern-0em} 2},{{ - z} \mathord{\left/ {\vphantom {{ - z} 2}} \right. \kern-0em} 2}} \right\}$. Изложенный метод определения эффективных свойств композита дает аналитические оценки эффективных модулей упругости и поэтому является привлекательным для метода комплексных модулей упругости. Модули вязкоупругих компонент композита задаются в комплексной форме. В результате, эффективные модули упругости также вычисляются аналитически как комплексные модули, определяющие эффективные свойства демпфирования.

Определим условия, при которых реализуются механизмы демпфирования, связанные со сдвиговыми и объемными деформациями [11]. Предварительный анализ приведем на примере модуля потерь поперечного сдвига $\mu _{{23}}^{{''}}$ волокнистого композита с осью изотропии $x = {{x}_{1}}$, направленной вдоль волокон (плоскость ${{x}_{2}}{{x}_{3}}$ является плоскостью изотропии трансверсально изотропного эффективного материала).

На рис. 1 представлены зависимости для модуля потерь поперечного сдвига $\mu _{{23}}^{{''}}$, построенные для различных значений модуля сдвига ${{\mu }_{2}}$ вязкоупругого слоя. Представленные результаты позволяют сделать ряд выводов.

Рис. 1.

Зависимости эффективного модуля потерь $\mu _{{23}}^{{''}}$ волокнистого композита от относительной толщины $\Delta {\text{/}}R$ вязкоупругого покрытия волокон радиуса $R$ со свойствами волокна ${{K}_{1}} = 40$ МПа, ${{\mu }_{1}} = b = 30$ МПа, матрицы ${{K}_{m}} = 4$ МПа, ${{\mu }_{m}} = 3$ МПа и различными значениями модуля сдвига вязкоупругого покрытия ${{\mu }_{2}}$ ($( \cdot \cdot \cdot \cdot \cdot )$ ${{\mu }_{2}} = (1 + i)$, $( - - - )$ ${{\mu }_{2}} = 0.1(1 + i)$, $( - \cdot - \cdot - )$ ${{\mu }_{2}} = 0.01(1 + i)$, ${{\mu }_{2}} = 0.001(1 + i)$, и ${{\mu }_{2}} = 0.0005(1 + i)$).

1. Вид зависимостей модуля потерь модуля поперечного сдвига $\mu _{{23}}^{{''}}$ существенно зависит от отношения $a{\text{/}}b$, где ${{\mu }_{1}} = b,$ ${{\mu }_{2}} = a(1 + i\eta )$ – модули упругости жесткой фазы (волокна) и вязкоупругого слоя. Для малых значений a/b первый максимум на рассматриваемых зависимостях появляется при достаточно большой толщине вязкоупругого слоя, и даже может отсутствовать, если $0.5 < a{\text{/}}b < 1$, указывая на достаточно очевидный путь увеличения эффективного модуля потерь за счет увеличения объемной доли вязкоупругого слоя. С ростом объемного содержания “слабого” слоя, очевидно, имеет место значительное снижение эффективных механических характеристик композита.

С уменьшением a/b (${{\mu }_{1}} = b$ является постоянной, изменяется параметр a) максимум модуля потерь сдвигается в область малых значений толщин вязкоупругого слоя $\Delta {\text{/}}R$, что представляется весьма привлекательным при проектировании композитов с высокой степенью диссипации и высокими механическими характеристиками. Анализ зависимости модуля потерь от параметра a/b и толщины покрытия показывает, что в исследуемой композитной структуре на разных толщинах покрытия [11] проявляются два различных механизма диссипации. Оказывается, что на малых толщинах покрытия, т.е. при малых концентрациях вязкоупругого “слабого” слоя, могут реализоваться высокие диссипативные свойства, связанные со сдвигом, когда эффективный модуль потерь почти на два порядка превышает модуль потерь вязкоупругого слоя, являющегося первопричиной диссипативных процессов. Действительно, достаточно сравнить модули потерь вязкоупругого слоя $\mu _{2}^{{''}} = 0.0005$ со значением эффективного модуля потерь, близкого к 0.05. Проявляется синергетический эффект. Удивительно, что этот пик и соответствующие оптимальные значения толщин покрытия не были обнаружены прежде в исследованиях механиков. Оценка оптимальных значений толщин покрытия является главной целью исследования.

2. Из представленных зависимостей видно, что значение максимума второго пика для малых a/b фактически не зависит от собственной вязкости $\eta $, ${{\mu }_{2}} = a(1 + i\eta )$, вязкоупругого слоя. Отметим, что характеристики сдвига реального вязкоупругого полимера при температуре выше температуры стеклования Tg и при температуре стеклования Tg равны, соответственно, значениям μ2 = = $0.001(1 + 0.1i)$ и ${{\mu }_{2}} = 0.01(1 + i)$.

3. Первый пик, для которого отношение a/b является относительно большим, даже в случае среды, близкой к вязкой жидкости, соответствует механизму объемного деформирования. Второй пик соответствует механизму сдвига. Для него характерна малая величина оптимальных толщин покрытий, что подтверждается путем прямого расчета энергии диссипации в единице объема Ù = $1{\text{/}}E{\kern 1pt} '{\kern 1pt} 'e_{0}^{2}$. Расчет показывает, что в тонких покрытиях имеет место колоссальное увеличение сдвиговых деформаций, которым соответствует высокая плотность энергии потерь.

ПРИБЛИЖЕННАЯ ОЦЕНКА ПАРАМЕТРОВ СТРУКТУРЫ КОМПОЗИТОВ С ВЫСОКОЙ ДИССИПАТИВНОСТЬЮ

Для качественного описания диссипативных эффектов изучим диссипативные свойства слоистой двухфазной системы, состоящей из упругого и вязкоупругого слоя. В этом случае трансверсальные модули определяются формулами Рейса. Для модуля сдвига имеем ${{({{\mu }_{{{\text{eff}}}}})}^{{ - 1}}}\, = \,f{{({{\mu }_{2}})}^{{ - 1}}}$ + (1 – ‒ f  )(μ1)–1, где $f$ – объемное содержание (относительная толщина) вязкоупругого слоя, ${{\mu }_{1}} = b$ – модуль сдвига упругого слоя и ${{\mu }_{2}} = a(1 + i\eta )$ – модуль сдвига вязкоупругого слоя. Запишем явное выражение для модуля потерь. С использованием формулы Рейса найдем μ'' = $\frac{{{{b}^{2}}s\eta }}{{{{{(1\, + \,bs\, - \,f)}}^{2}}\, + \,{{{(\eta \, - \,f\eta )}}^{2}}}}$, $f{\text{/}}a = s$. Используя в последнем выражении разложение по малому параметру f, после несложных преобразований получим

(3)
$\mu {\kern 1pt} '{\kern 1pt} ' = B(f)\frac{{{{b}^{2}}s\eta }}{{{{{(1 + bs)}}^{2}} + {{\eta }^{2}}}},$
где                      $B(f) = 1 + \frac{{2(1 + bs + {{\eta }^{2}})f}}{{{{{(1 + bs)}}^{2}} + {{\eta }^{2}}}}$  +            + $\left( {\frac{{4{{{(1 + bs + {{\eta }^{2}})}}^{2}}}}{{{{{({{{(1 + bs)}}^{2}} + {{\eta }^{2}})}}^{2}}}} - \frac{{1 + {{\eta }^{2}}}}{{{{{(1 + bs)}}^{2}} + {{\eta }^{2}}}}} \right){{f}^{2}} + ...$

Функция $B(f),$ $f < 1$ в (3) является монотонной положительной функцией для всех физически обоснованных параметров $a,\;b,\;\eta $. В таком случае экстремальное значение для модуля потерь $\mu ''$ реализуется при ${{f}_{{\max }}}$, определяющем экстремум функции $\varphi (f) = \frac{{{{b}^{2}}s\eta }}{{{{{(1 + bs)}}^{2}} + {{\eta }^{2}}}},$ $s = f{\text{/}}a$. Точное значение ${{f}_{{\max }}}$, при котором имеет место экстремум функции $\varphi (f)$, равно

(4)
${{f}_{{\max }}} = {{(1 + {{\eta }^{2}})}^{{1/2}}}a{\text{/}}b.$

Оно совпадает с асимптотическим равенством, полученным в работе [11] при $a{\text{/}}b \to 0$. Следовательно, доказано, что диапазон использования асимптотического выражения (4) расширяется и формулу (4) можно применять для оценок оптимальных значений толщины вязкоупругого слоя во всем диапазоне параметров $a,\;b,\;\eta $. Непосредственно из (4) и условия ${{f}_{{\max }}} < 1$ следует, что максимум для модуля потерь может реализоваться только при условии $a{\text{/}}b < {{(1 + {{\eta }^{2}})}^{{ - 1/2}}}$. В противном случае рост диссипативных свойств наблюдается лишь с увеличением толщины вязкоупругого слоя.

На рис. 2а, б представлены зависимости модуля потерь, построенные с привлечением формул (4), которые позволяют судить об особенностях и механизмах демпфирования. Данные зависимости демонстрируют смещение максимумов модуля потерь в область меньших концентраций ${{f}_{{\max }}}$ при уменьшении отношения a/b (рис. 2 a, б) с фактическим сохранением значения максимумов, которые пропорциональны только величине параметра $b$ (рис. 2б, кривые 1, 2).

Рис. 2.

Зависимости модуля потерь $\mu {\kern 1pt} '{\kern 1pt} '$ для диапазона $0 < {{f}_{{\max }}} < 0.2$ (а) и диапазона $0 < {{f}_{{\max }}} < 0.05$ (б), полученные для следующих параметров: $( \cdot \cdot \cdot \cdot \cdot )$ $a{\text{/}}b = 0.001$, $( - - - )$ $a{\text{/}}b = 0.01$ и $ - a{\text{/}}b = 0.1$, $\eta = 0.1$; $b = 40$ ГПа (1) и $b = 20$ ГПа (2).

Максимум модуля потерь для $a{\text{/}}b = 0.1$ соответствует характеристикам поперечного деформирования по отношению к плоскости слоев (например, при ${{E}_{1}} = b = 40$ ГПа, ${{E}_{2}} = 4(1 + 0.1i)$ ГПа) и реализуется для сравнительно больших значений толщины вязкоупругого слоя ${{f}_{{\max }}} \approx 0.1$. В этом случае динамический модуль падает в два раза по сравнению с модулем упругого слоя. С другой стороны, для сдвиговых деформаций характерны малые оптимальные значения толщин вязкоупругого слоя, при которых реализуются аномально высокие диссипативные характеристики, – эффективный модуль потерь превышает модуль потерь вязкого слоя фактически в 100 раз при сохранении высоких механических свойств. Возможность появления двух пиков для дисперсных композитов связано с тем, что для них могут реализовываться оба механизма диссипации, связанные и с объемным, и со сдвиговым характером деформирования. Удивительно, что наиболее эффективный сдвиговой механизм диссипации был упущен исследователями, вероятно, потому, что он проявлялся на очень малых концентрациях – толщинах вязкоупругого покрытия.

Наконец, отметим эффект, имеющий отношение к проектированию композитных структур с высокими демпфирующими свойствами. Сравнение результатов для модулей потерь, полученных для композитов со сферическими и цилиндрическими включениями на основе самосогласованного метода Эшелби, с результатами моделирования слоистых структур показывает, что аналитическое равенство (4), полученное в результате исследования слоистых структур, может быть использовано и для прогноза толщин вязкоупругих покрытий дисперсных и волокнистых композитов. При этом соотношение (4) является фактически точным для случая, когда включение и матрица являются упругими и выполнены из одного материала. В общем случае, когда свойства жестких упругих включений отличны от свойств матрицы, формула (4) также дает хорошие оценки оптимальных толщин покрытий, если в них вместо $b$ использовать $\bar {b} = b{{{v}}_{f}}$, где $b$ – модуль упругости упругих включений, ${{{v}}_{f}}$ – объемная концентрация волокон в композите. Оценка эффективного модуля потерь должна проводиться по характеристикам, найденным в рамках самосогласованного метода трех фаз Эшелби.

ВЛИЯНИЕ ПЕРЕМЕННОСТИ ТОЛЩИНЫ ВЯЗКОУПРУГОГО СЛОЯ НА ДЕМПФИРУЮЩИЕ СВОЙСТВА

Оценим влияние возможной переменности толщины покрытия на диссипативные свойства композитов. Считаем, что для этого достаточно рассмотреть слоистую структуру, ибо модель слоистого композита с высокой точностью позволяет моделировать и диссипативные свойства гибридных дисперсных и волокнистых композитов. Полагаем, что относительная толщина (концентрация) вязкоупругого слоя записывается в виде $f = \bar {f} + \varepsilon \sin \omega x$, где $\bar {f}$ – среднее значение концентрации вязкоупругого слоя, а параметры $\omega $ и $\varepsilon $ определяют периодическое отклонение концентрации от средней величины. Относительная толщина упругих слоев может быть записана в виде ${{h}_{1}} = 1 - f = 1 - \bar {f} + \varepsilon \sin \omega x$. Тогда выражение для эффективного модуля сдвига $\varepsilon $ представляется выражением

$\begin{gathered} \mu = {{(A + B\varepsilon \sin \omega x)}^{{ - 1}}}, \\ A = \left[ {\frac{{\bar {f}}}{{{{\mu }_{2}}}} + \frac{{1 - \bar {f}}}{{{{\mu }_{1}}}}} \right] = \mu _{0}^{{ - 1}},\quad B = \left[ {\frac{1}{{{{\mu }_{1}}}} - \frac{1}{{{{\mu }_{2}}}}} \right]. \\ \end{gathered} $

Усредняя последнее выражение по периоду неоднородности, получим

(5)
$\begin{gathered} {{\mu }_{{av}}} = {{A}^{{ - 1}}}{{(1 - {{B}^{2}}{{\varepsilon }^{2}}{\text{/}}{{A}^{2}})}^{{ - 1/2}}}, \\ {{\mu }_{{av}}} = {{(2\pi {\text{/}}\omega )}^{{ - 1}}}\int\limits_0^{2\pi /\omega } {{{{(A + B\varepsilon \sin \omega x)}}^{{ - 1}}}dx} . \\ \end{gathered} $

В (5) множитель ${{A}^{{ - 1}}}$ определяет эффективный модуль ${{A}^{{ - 1}}} = \mu _{0}^{{}}$, записанный в терминах постоянных – осредненных геометрических характеристик структуры. Для ${{A}^{{ - 1}}} = \mu _{0}^{{}}$ оптимальное значение слоистой структуры дается соотношением (4), ${{f}_{{\max }}} = {{(1 + {{\eta }^{2}})}^{{1/2}}}a{\text{/}}b$. Функция ${{(1 - {{B}^{2}}{{\varepsilon }^{2}}{\text{/}}{{A}^{2}})}^{{ - 1/2}}}$ содержит параметр, учитывающий возможную переменность параметров вязкоупругих слоев. Численный анализ поведения этой функции показал, что действительная часть этой функции $G(\bar {f},\varepsilon ) = \operatorname{Re} [{{(1 - {{B}^{2}}{{\varepsilon }^{2}}{\text{/}}{{A}^{2}})}^{{ - 1/2}}})$ во всем диапазоне объемных содержаний является монотонной положительной функцией, близкой к единице (см. рис. 3a), даже для сравнительно больших значений параметра $\alpha = \varepsilon {\text{/}}\bar {f}$, характеризующего возможные относительные отклонения толщин вязкоупругого слоя.

Рис. 3.

Зависимость функции $G(\bar {f},\varepsilon )$ (а) и модуля потерь $\mu {\kern 1pt} '{\kern 1pt} '$ (б) от средней постоянной концентрации $\bar {f}$ и различных значений относительных отклонений толщин вязкоупругого слоя $\alpha = \varepsilon {\text{/}}\bar {f}$ при $a{\text{/}}b = 0.00025,$ $d = 80$ ГПа, $\eta = 1$.

Функция $F(\bar {f},\varepsilon ) = \operatorname{Im} [{{(1 - {{B}^{2}}{{\varepsilon }^{2}}{\text{/}}{{A}^{2}})}^{{ - 1/2}}})$ ведет себя при $\varepsilon \to 0$ как $O({{\varepsilon }^{2}})$ во всем диапазоне объемных содержаний, а ее максимальное значение при $a{\text{/}}b < 0.01$ и $f < 0.1$ не превышает величины 0.001. Приведенные оценки показывают (рис. 3б), что в реальном диапазоне свойств материалов переменность параметров вязкоупругого слоя не оказывает существенного влияния на диссипативные характеристики композита.

ЗАКЛЮЧЕНИЕ

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

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

  1. Finegan I., Gibson R. Recent research on enhancement of damping in polymer composites // Compos. Struct. 1999. V. 44. P. 89–98.

  2. Chandra R., Singh S., Gupta K. Damping studies in fiber-reinforced composites – a review // Compos. Struct. 1999. V. 46. P. 41–51.

  3. Remillit C. Damping mechanism of polymers filled with elastic particles // J. Mech. Mater. 2007. V. 39. P. 525–537.

  4. Meaud J., Sain T., Hulbert G.M., Waas A.M. Analysis and optimal design of layered composites with high stiffness and high damping // Int. J. Solids Structures. 2013. V. 50. P. 1342–1353.

  5. Lakes R.S. High damping composite material. Effect of structural hierarchy // J. Compos. Mater. 2003. V. 36(3). P. 287–297.

  6. Berthelot J.M., Assarar M., Sefrani Y., Mahi A.E. Damping of composite materials and structures // Compos. Struct. 2008. V. 85. P. 189–204.

  7. Yu L., Ma Y., Zhou C., Xu H. Damping efficiency of the coating structure // Int. J. Solids Structures. 2005. V. 42. P. 3045–3058.

  8. Chandra R., Singh S.P., Gupta K. A study of damping in fiber-reinforced composites // J. Sound Vib. 2003. V. 262. P. 475–496.

  9. Fisher F.T., Brinson L.C. Viscoelastic interphases in polymer-matrix composites: theoretical model and finite element analysis // Compos. Sci. Technol. 2001. V. 61. P. 731–748.

  10. Wei P.J., Huang Z.P. Dynamic effective properties of the particle reinforced composites with viscoelastic interphase // Int. J. Solids Structures. 2004. V. 41. P. 6993–7007.

  11. Gusev A.A., Lurie S.A. Loss amplification effect in multiphase materials with viscoelastic interfaces // Macromolecules. 2009. V. 42 (14). P. 5372–5377.

  12. Sain T., Meaud J., Hulbert G., Arruda E.M., Waas A.M. Simultaneously high stiffness and damping in a class of wavy layered composites // Compos. Struct. 2013. V. 101. P. 104–110.

  13. Christensen R.M., Lo K.H. Solutions for effective shear properties in three phase sphere and cylinder models // J. Mech. Phys. Solids. 1979. V. 27 (4). P. 315–330.

  14. Папкович П.Ф. Теория упругости. М., Л.: Оборонгиз, 1939. 640 с.

Дополнительные материалы отсутствуют.

Инструменты

Доклады Российской академии наук. Физика, технические науки