г. Москва и Московская область, Россия
ГРНТИ 67.09 Строительные материалы и изделия
ББК 383 Строительные материалы и изделия
Композиционные материалы, состоящие из нескольких фаз, широко используются в современном строительстве. В статье рассмотрены механические характеристики упруго-ползучих слоистых композиционных материалов. Каждая из составляющих фаз обладает свойствами упругости, вязкости или ползучести. Многочисленные эксперименты показали, что свойства структурно неоднородных материалов могут значительно отличаться от свойств отдельных компонент, входящих в их состав. Кроме того, быстро осциллирующие коэффициенты дифференциальных уравнений, описывающих такие композиционные материалы, значительно усложняют решение краевых задач даже с помощью современных компьютерных методов. Поэтому для решения таких задач используется метод усреднения, после применения которого исходный неоднородный материал, содержащий изотропные слои, превращается в однородную трансверсально-изотропную среду. Ядра ползучести слоев задаются суммой конечно-го числа убывающих экспоненциальных функций с коэффициентами одного знака. В этом случае для иллюстрации можно использовать последовательное и параллельное соединения простейших, “наивных” механических моделей вязкоупругого поведения Максвелла и Фойхта. Экспериментально обосновано использование ядер ползучести этого типа. В данной работе показано, что эффективная (усредненная) модель для композиционного материала с указанными свойствами представляет однородную (гомогенизированную) среду, которая описывается ядром ползучести, также представленным конечной суммой экспоненциальных функций. Для слоистого композита, состоящего из слоев двух изотропных упруго-ползучих материалов с ядрами ползучести, соответствующих “наивной” вязкоупругости, усредненный (эффективный) модуль растяжения или сжатия в направлении, перпендикулярном слоям, также будет соответствовать “наивной” вязкоупругости. Предложен алгоритм быстрого и точного расчета усредненных ядер ползучести гомогенизированного материала.
композиционный материал, усреднение, гомогенизация, слоистый упруго-ползучий материал, спектральное ядро ползучести.
Введение. Задачи для неоднородных (гетерогенных) материалов возникают во многих областях строительства. Многие строительные материалы являются слоистыми, что является причиной их специфического поведения при деформировании. С помощью метода асимптотического усреднения такой материал преобразуется в однородный, который подчиняется усредненным уравнениям. Этот метод усреднения был разработан в [1, 2], одним из первых его применений были задачи теории упругости [3, 4]. В качестве примера, задача усреднения неоднородной упругой среды с использованием теоретических результатов была решена в [5]. Кроме того, напряженно-деформированное состояние композиционных материалов изменяется в течение длительных интервалов времени. Это следует учитывать при исследовании прочности таких материалов. Проблема изучения реологических свойств актуальна для строительной практики, поскольку ползучесть действует на перераспределение напряжений и в некоторых случаях может привести к недопустимому увеличению деформации. В [6–9] изложены основы теории наследственной механики Больцмана-Вольтерра. Эта теория описывает процессы, в которых состояние механической системы зависит от всего периода совершаемых над ней действий. Приведены примеры решения задач линейной наследственной теории упругости и нелинейной теории ползучести. В работе [10] построена усредненная модель, описывающая совместное движение слоев, составленных из упругого и вязкоупругого материалов. Получены усредненные уравнения состояния слоистого композита, состоящего из слоев попарно чередующихся ползучих материалов [11]. В работе [12] даны эффективные упругие модули слоистой упруго-ползучей среды.
В статье предлагается метод моделирования напряженно-деформированного состояния слоистых материалов на основе сочетания методов теории усреднения и теории ползучести. В качестве основных реологических соотношений приняты соотношения теории линейной наследственности Больцмана-Вольтерра, согласно которой деформация в момент времени t имеет вид:
, (1)
где E – модуль Юнга, – напряжение,
– ядро ползучести, являющееся убывающей функцией [8].
Параметры реологической модели определяются на основе лабораторных испытаний. Ядра ползучести для каждого слоя являются конечными суммами экспоненциальных функций времени. Эксперименты доказывают, что эти функции хорошо описывают ползучесть материалов в течение длительного периода времени. Такие ядра называются спектральными. Кроме того, в ряде работ используются и другие типы ядер, например, ядра Абеля.
Очевидно, что модели могут быть различными для различных материалов, используемых в инженерной практике, таких, как металлы и сплавы, полимеры, волоконные композиты (с полимерной или металлической матрицей), бетон и дерево. Тем не менее можно использовать те же самые принципы для установления уравнений состояния для этих материалов.
Модели вязкоупругих материалов могут быть построены путем объединения упругих и вязких элементов. Каждый упругий элемент может быть представлен пружинкой с параметром E, который является его модулем Юнга. Каждый вязкий элемент может быть представлен демпфером с параметром , который представляет собой вязкость материала. Можно построить модель линейной вязкоупругости, рассматривая комбинации упругой пружины и демпфера. Такая модель вязкоупругого тела называется “наивной” [6].
Например, модель Максвелла может быть представлена чисто вязким демпфером и чисто упругой пружинкой, соединенных последовательно, как показано на рисунке 1. Другая двухэлементная модель, модель Фойгта, состоит из пружинки и демпфера, соединенных параллельно, как показано на рисунке 2.
Рис. 1. Модель Максвелла Рис.2. Модель Фойгта
Модель Максвелла представляет собой вязкоупругую жидкость, модель Фойгта представляет вязкоупругое твердое вещество. Таким образом, модели Фойгта и Максвелла являются простейшими вязкоупругими телами. Модель Максвелла представляет экспоненциальную релаксацию напряжений и линейную ползучесть деформации; ее также называют релаксационным элементом.
Модель Фойгта демонстрирует экспоненциальную ползучесть деформации. Модель применяется при моделировании ползучести в материалах. Чтобы описать реологическое поведение сложного материала, можно комбинировать разными способами модели этих простейших идеальных тел. Существуют также обобщения механических моделей, позволяющих описать физическую нелинейность. В работе [14] модель Максвелла, состоящая из нелинейной пружины, соединенной последовательно с нелинейным демпфером, подчиняющимся степенному закону с постоянными параметрами материала, предложена для представления зависящих от времени свойств различных вязкоупругих материалов. В работе [15] описано численное моделирование напряженно-деформированного состояния в грунте, подвергнутом волновым процессам. Модели Фойгта и Максвелла используются для описания деформаций грунтового основания и контактирующей с ним бетонной конструкции. Авторы создали свой собственный компьютерный код на основе метода конечных элементов (FEM).
Хорошо известно, что при моделировании упруго-ползучих или вязкоупругих материалов с помощью комбинаций пружин и демпферов передаточная функция (множитель, связывающий преобразования Лапласа напряжений и деформаций) представляет собой сумму простейших дробей с коэффициентами одного знака [6]. Целью настоящей статьи является исследование проблемы сохранения свойства передаточной функции для усредненной модели слоистого композита. Это исследование ответит на следующий вопрос: для каких деформаций слоистого композита мы можем использовать аппроксимирующую модель, состоящую из пружин и демпферов, упомянутых выше.
Постановка задачи и решение. Рассмотрим область , состоящую из попарно чередующихся слоев упруго-ползучих изотропных материалов. Пусть
- это отношение длины двух соседних слоев к характерному размеру образца, а толщины каждого отдельного слоя соответственно равны:
и
,
. Выберем декартову систему координат
. Все слои рассматриваемой области
параллельны координатной плоскости
. Компоненты тензора напряжений
определяются не только деформацией в данный момент, но и всей предшествующей историей деформации тела. Поэтому уравнения состояния, связывающие компоненты тензоров деформации и напряжений для каждого слоя, следующие [6, 7]:
(2)
где
,
- номер слоя),
– компоненты вектора перемещений,
– компоненты тензора напряжений,
– дельта-функция Дирака,
– интегральные операторы Вольтерра, то есть
, (3)
и переменная t обозначает время. (Используется соглашение Эйнштейна для повторяющихся индексов.)
Ядра релаксации зависят от разности
. Это следует из условия инвариантности величины
по отношению к началу отсчета времени t. В данной статье ядра релаксации представлены в виде сумм экспоненциальных функций, так как такие ядра, как показывают эксперименты, хорошо описывают длительные деформационные процессы.
На горизонтальных поверхностях слоев принимаются условия идеального контакта: компоненты перемещений и компоненты напряжений, параллельные оси , непрерывны, т.е.
.
Уравнения равновесия элемента сплошной среды имеют вид [9]:
, (4)
где , и
- компоненты вектора внешних сил.
Для изотропных материалов компоненты тензора упругости и ядер релаксации в (2) имеют следующий вид [6]:
– постоянные Ламе,
и
– регулярные части сдвиговой и объемной релаксации соответственно,
– символ Кронекера. Предположим, что для каждого слоя амплитуда тензора объемной релаксации пропорциональна амплитуде сдвиговой с коэффициентом пропорциональности
, то есть:
,
const,
. Далее
обозначим
.
В задаче все модули упругости и коэффициенты ядер релаксации являются периодическими функциями координаты 
(
- относительный период ячейки) и являются кусочно-постоянными функциями этой переменной, т. е. имеют следующий вид [12]:
.
Выберем ядро ползучести экспоненциального типа для каждого слоя (s – номер слоя, ):
,
- постоянные,
,
[6].
Применим к уравнениям (4) с учетом (2) преобразование Лапласа во временной области
, (5)
В результате получим систему уравнений теории упругости с комплексным параметром р, к которой применяем метод усреднения, описанный в [1, 2]. После этого имеем однородную анизотропную среду с усредненными (эффективными) модулями [12]. Изучим некоторые качественные свойства полученных эффективных модулей. Для этого необходимы следующие два утверждения.
Утверждение 1. Пусть
, (6)
где, как известно из теории вязкоупругости [3],
Нули и полюса
функции
изображены на Рис.3. Тогда обратная функция имеет подобный вид
, где
простые нули функции
на отрезке
,
– простой нуль на отрезке
, как изображено на рис.4, и
- положительные постоянные.
Рис. 3. Нули и полюса
функции
Рис. 4. Нули и полюса
функции
Утверждение 2. Пусть две функции вида (6). Рассмотрим функцию
, где h постоянная,
. Эту функцию можно представить в следующем виде:
, где
.
– нули функции
на отрезках, у которых концы – это точки множеств
,
– нули функции
между ее полюсами, и один из нулей находится между полюсом, ближайшем к началу координат и началом координат
. То есть мы берем отрезки, у которых нет внутренних точек из набора
. Внутри каждого из этих отрезков имеется один корень функции
, как показано на
рис. 5.
Рис. 5. Нули функций
Кроме того, , где
- такая же постоянная для функции
, как
для
.
Теперь мы используем приведенные выше утверждения для анализа качественных свойств усредненного тензора упруго-ползучей среды. После применения преобразования Лапласа система уравнений равновесия для усредненной среды принимает вид
где усредненные модули определены в [2]. В частности, если материалы двух упруго-ползучих слоев изотропны, то модуль, определяющий эффективные упруго-ползучие свойства композиционного материала для перемещений, перпендикулярных плоскости слоев, имеет вид [13]:
,
где – функции вида (6), принимающие постоянные значения в каждом слое, переменная z заменяется на p,
- среднее значение функции
, которая принимает постоянные значения
и
на соседних слоях соответственно.
Для сравнения один из модулей сдвига имеет вид [13]: , но если перемножить две функции вида (6) и разложить результат на простые дроби, то коэффициенты этих дробей могут иметь разные знаки. Это противоречит “наивной” теории вязкоупругости.
Прямое использование утверждений 1 и 2 приводит к следующим выводам.
Выводы. Для слоистого композита, состоящего из слоев двух изотропных упруго-ползучих материалов с ядрами ползучести, соответствующих “наивной” вязкоупругости, усредненный (эффективный) модуль растяжения или сжатия в направлении, перпендикулярном слоям, также будет соответствовать “наивной” вязкоупругости. Для таких эффективных модулей новые показатели экспонент можно легко найти, например, методом половинного деления, а коэффициенты экспоненциальных функций из решения систем линейных уравнений. Следует отметить, что для других усредненных модулей утверждение о неизменяемости свойства “наивности” неверно. Предлагаемый метод допускает любые напряжения внутри композитного материала. Слоистая структура с параллельными слоями используется только при построении усредненной (эффективной) модели. После того, как коэффициенты уже найдены в явном аналитическом виде, напряжения внутри рассматриваемого композита не обязаны удовлетворять никаким дополнительным требованиям, в частности перпендикулярности слоям рассматриваемого материала. Предложенная модель может быть использована в любых условиях многоосного напряженно-деформированного состояния.
1. Олейник О.А., Иосифьян Г.А., Шамаев А.С. Математические задачи теории сильно неоднородных сред. М. Изд-во МГУ, 1990. 311 с.
2. Бардзокас Д.И., Зобнин А.И. Математическое моделирование физических процессов в композиционных материалах периодической структуры. М. Едиторал УРСС, 2003. 376 с.
3. Победря Б.Е. Механика композиционных материалов. М. Изд-во МГУ, 1984. 336 с.
4. Christensen R.M. Mechanics of Composite Materials. New York. J. Wiley-Int.-science, 1979. 348 p.
5. Bobyleva T.N. Approximate Method of Calculating Stresses in Layered Array // Procedia Engineering. 2016. V. 153. P. 103-106. DOI:https://doi.org/10.1016/j.proeng.2016.08.087
6. Ильюшин А.А., Победря Б.Е. Основы математической теории термовязкоупругости. М. Наука, 1970. 281 с.
7. Christensen R.M. Theory of viscoelasticity. New York. Dover Publ, 2010. 384 p.
8. Работнов Ю.Н. Элементы наследственной механики твердых тел. М. Наука, 1977. 384 с.
9. Работнов Ю.Н. Механика деформируемого твердого тела. М. Наука, 1988. 712 с.
10. Шамаев А.С., Шумилова В.В. Асимптотическое поведение спектра одномерных колебаний в среде из слоев упругого материала и вязкоупругого материала Кельвина-Фойгта // Современные проблемы механики. Сборник статей. Тр. МИАН. М. МАИК Наука / Интерпериодика. 2016. Т. 295. С. 218-228. DOI: https://doi.org/10.1134/S0371968516040130
11. Шамаев А.С., Шумилова В.В. Усреднение уравнений состояния для гетерогенной среды, состоящей из слоев двух ползучих материалов // Современные проблемы механики. Сборник статей. Тр. МИАН. М. МАИК Наука/Интерпериодика. 2016. Т. 295. 229-240. DOI: https://doi.org/10.1134/S0371968516040142
12. Бобылева Т.Н., Шамаев А.С. Эффективный алгоритм расчета усредненных характеристик слоистых грунтов, состоящих из упруго-ползучих материалов // Основания, фундаменты и механика грунтов. 2017. № 4. С. 7-12.
13. Bobyleva T.N. Method of Calcu-lation of Stresses in the Layered Elastic-Creeping Arrays. MATEC Web of Conf. 2016. V. 86. P. 01024. DOI: https:// doi.org/10.1051/matecconf/20168601024
14. Monsia M.D. A Simplified Nonlinear Generalized Maxwell Model for Predicting the Time Dependent Behavior of Viscoelastic Materials // World Journal of Mechanics. 2011. № 1. P. 158-167. DOI:https://doi.org/10.4236/wjm.2011.13021
15. Шешенин С.В., Закалюкина И.М., Коваль С.В. Численная реализация моделей Фойгта и Максвелла для моделирования волн в грунте // Вестник МГСУ. 2014. № 11. С.82-89. DOI:https://doi.org/10.22227/1997-0935.2014.11.82-89.
