Полный текст

Введение Актуальность темы обусловлена использованием высоконаполненных полимерных композитов в ответственных конструкциях аэрокосмической техники и других отраслях [1-3]; действием сложных гармонических нагрузок на подобные конструкции вследствие всевозможных ударных воздействий, транспортировки, пульсаций внутрикамерного давления и т.д. [4, 5]; необходимостью развития методов экспериментального исследования и определения деформационных свойств материалов, а также методов расчета конструкций, работающих в экстремальных условиях. В целом работа направлена на разработку методик проведения динамического опыта, определение вязкоупругих параметров высоконаполненных полимерных композитов при стационарных двухчастотных воздействиях и процедуры идентификации используемой многофакторной математической модели для оценки напряженно-деформированного состояния вязкоупругих конструкций аэрокосмической техники. 1. Математическая модель описания поведения вязкоупругих материалов при стационарных двухчастотных нагрузках Одним из наиболее универсальных для описания линейных вязкоупругих свойств является способ, основанный на применении линейных интегральных уравнений Больцмана-Вольтерра [6, 7]. В случае одноосного растяжения-сжатия можно записать , (1) где τ - время, предшествующее моменту времени наблюдения t; Е(t) - функция релаксации (ядро релаксации). Бесконечный нижний предел интегрирования подразумевает учет полной истории деформирования. На практике бывают случаи, когда вязкоупругие конструкции работают в условиях стационарных гармонических колебаний. Описанные выше соотношения между напряжениями и деформациями применимы и при таких условиях эксплуатации. Однако при описании поведения вязкоупругого материала удобнее использовать комплексные параметры [5-8]. Пусть далее деформации от времени изменяется по гармоническому закону вида , (2) где εа - амплитуда деформации; ω - угловая частота. Подставим (2) в (1), в результате получим . (3) Путем замены переменной t - τ = η (3) можно записать в виде . (4) Воспользуемся формулой Эйлера [9] (5) и перепишем (4) в форме . (6) Далее будем рассматривать такие ядра релаксации, для которых интегралы из уравнения (6) сходятся (например, экспоненциальное ядро при β > 0). Таким образом, ограничения, накладываемые на ядра, фактически сводятся к тому, что все полученные уравнения будут справедливы для асимптотически затухающей модели вязкоупругого деформирования при σ = const. Примем , , , (7) , (8) , (9) , , где - комплексный параметр; Е' - действительная часть комплексного параметра; Е" - мнимая часть комплексного параметра, φE; E* - вязкоупругие параметры. Таким образом, (6) можно записать в виде . (10) Подставив (8) и (9) в (7) и, воспользовавшись формулой Эйлера, получим . (11) Затем перепишем (10) с учетом (11): . (12) Интерпретация (12) важна с физической точки зрения. При стационарных гармонических колебаниях деформация запаздывает от напряжения на величину, определяемую углом потерь φE [7]. Наполненным полимерам свойственно нелинейное поведение: отличие петли гистерезиса от эллиптической формы [6, 10], размягчение материала (эффект Маллинза-Патрикеева) [11, 12], влияние амплитуды деформации (эффект Пэйна) [12, 13] и предварительной статической деформации на определяемые вязкоупругие характеристики [14-16]. Общее описание нелинейной математической модели вязкоупругого поведения материала было получено Вольтерра с использованием представления, разработанного Фреше. Уравнение Вольтерра-Фреше [8, 17, 18] является общей формой записи физически нелинейных операторов вязкоупругой среды. Для одного измерения его можно записать в следующем виде: (13) Вследствие того что зависимости вязкоупругих параметров от предварительной статической деформации для данного типа материалов с достаточной точностью описываются полиномами второй степени [19], в дальнейшем будем ограничиваться только первыми тремя членами ряда (13). В случае действия стационарных двухчастотных колебаний деформации от времени изменяются по гармоническому закону вида , (14) где индекс 1 соответствует низкочастотной гармонике; индекс 2 - высокочастотной гармонике. Проделаем действия, аналогичные тем, которые были сделаны ранее (2)-(5). В результате всех преобразований получим (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) . (42) Анализ (15)-(42) показывает, что и , и , и , и , и , и , и , и , и имеют однотипные зависимости соответственно. Примем К = ω1 / ω2. При К = 1 двухчастотная нагрузка преобразуется в одночастотную нагрузку: ω1 = ω2 = ω, . При К = 0 двухчастотная нагрузка преобразуется в одночастотную нагрузку с предварительной статической деформацией: εст = εа1, ω1 = 0, ω2 = ω. В случае если то и и и и и и и и и имеют однотипные зависимости соответственно, которые сдвинуты относительно друг друга в полулогарифмических координатах на величину и на величину относительно и и и и и и и и и соответственно. В качестве функции релаксации можно использовать различные зависимости [8, 20]. Однако определение функции релаксации по результатам гармонических испытаний достаточно сложный и трудоемкий процесс, зависящий от условий эксперимента. С другой стороны, вязкоупругие параметры могут быть определены в условиях сложных воздействий, близких к эксплуатационным, например, путем регистрации зависимостей напряжений и деформаций в образце, испытывающем гармоническое деформирование вида , (43) где ν - частота (ν = ω / 2π); φg - угол сдвига начальных фаз. Зависимости вязкоупругих параметров от частоты и температуры могут быть описаны полиномами. Полиномиальные модели достаточно просты в практическом применении. Описание объекта с помощью таких моделей легко уточнить, повысив степень полинома. Для описания поведения вязкоупругого материала при различных значениях температуры может быть использована температурно-временная аналогия (суперпозиция, эквивалентность) [21-24] (в случае стационарных гармонических воздействий ее можно назвать температурно-частотной аналогией), называемая также гипотезой о термореологически простом поведении. Суть данного принципа заключается в том, что вязкоупругие характеристики (функции релаксации и ползучести, комплексные параметры) в полулогарифмических координатах (логарифм времени или частоты) при различных постоянных температурах часто имеют одинаковую форму кривых, но сдвинуты относительно друг друга по логарифмической шкале абсцисс. Если, приняв одну из температур за температуру приведения, сдвигом по логарифмической шкале можно сместить все кривые в одну без существенных погрешностей, тогда возможно использование данного принципа. Величины сдвига описываются функцией температурно-временного смещения. Значения вязкоупругих параметров могут быть определены с помощью преобразований Фурье аналогично работам [25, 26]. В условиях действия двухчастотных нагрузок (43) зависимость (15) может быть записана в следующем виде: (44) После раскрытия скобок получим (45) Таким образом, поведение материала при двухчастотных нагрузках описывается двенадцатью гармониками. В результате, были разработаны полиномиальные зависимости вязкоупругих параметров от температуры и частоты: , (46) , (47) , (48) , (49) , (50) , (51) , (52) , (53) , (54) , (55) , (56) , (57) , (58) , (59) , (60) , (61) , (62) , (63) , (64) , (65) , (66) , (67) где n1, …, n13 - степени полиномов, b1/0, …, b1/n1, …, b13/n13 - коэффициенты полиномов. Количество коэффициентов, их значения и степени полиномов определяются экспериментально. После выявления полиномиальных зависимостей (46)-(67) необходимо разработать оптимальный план экспериментальных исследований с использованием аппарата теории планирования экспериментов и на основе предложенного плана идентифицировать коэффициенты (46)-(67). После данной процедуры возможно оценить влияние различных вязкоупругих параметров на точность описания поведения материала и в случае их незначительности пренебречь ими. В более ранних работах [19, 26, 27] была предложена математическая модель описания поведения вязкоупругого материала при двухчастотных нагрузках: (68) для которой было принято допущение, что динамические нагрузки, действующие на фоне квазистатических, не приводят к большим деформациям, поэтому существенных искажений петли гистерезиса от эллиптической формы не происходит. Исходя из этого при одно- и двухчастотном нагружениях дополнительных частот (2ν1, 3ν1, 2ν2, 3ν2), помимо основных ν1 и ν2, не возникает. При этом остается существенной зависимость вязкоупругих параметров от амплитуд деформации εа1 и εа2. В работах [19, 28] данные зависимости были описаны линейными функциями. Отсутствие дополнительных частот не позволяет учесть различное сопротивление при растяжении и сжатии. Для учета разносопротивляемости материала был введен дополнительный параметр E*10, который смещает напряжения в сторону большего сопротивления (в область сжатия) [26]. Кроме того, при стационарных гармонических нагрузках для данного типа материалов наблюдаются зависимости вязкоупругих параметров от предварительной статической деформации [19, 29]. В связи с этим в случае двухчастотных воздействий появляются дополнительные гармоники (параметры E*23, E*33, E*34). Значениями параметра E*34 в (68) пренебрегаем. При εа1 > εа2 E*34 пренебрежимо мал, а при εа1 < εа2 отсутствие параметра E*34 компенсировалось изменениями E*10. Допущения о линейной зависимости вязкоупругих параметров от амплитуд деформации и об отсутствии искажений петли гистерезиса, а также пренебрежение параметром E*34 могут быть не характерны для других материалов. Так, например, если в результате исследования будут наблюдаться существенные искажения формы петли гистерезиса, модель (68) не сможет это учесть. В этом смысле математическая модель (44), представленная в данной статье, более предпочтительна к использованию. Резюмируя, скажем, что для исследуемого типа вязкоупругих материалов в той или иной степени наблюдается нелинейное поведение, а также различное сопротивление растяжению и сжатию (разносопротивляемость) [30]. Поведение таких материалов при двухчастотных нагрузках может быть описано восемнадцатью вязкоупругими параметрами E*1/1, E*2/1, E*3/1, E*2/3, E*3/3, E*3/4, E*1/2, E*2/2, E*3/2, φE1/1, φE2/1, φE3/1, φE2/3, φE3/3, φE3/4, φE1/2, φE2/2, φE3/2. В случае если нелинейные эффекты и разносопротивляемость не существенны, значения параметров E*2/1, E*3/1, E*2/3, E*3/3, E*3/4, E*2/2, E*3/2, φE2/1, φE3/1, φE2/3, φE3/3, φE3/4, φE2/2, φE3/2 будут стремиться к нулю. 2. Экспериментальные исследования Экспериментальные исследования были проведены в ЦКП «Центр экспериментальной механики ПНИПУ» на электродинамической испытательной системе Instron ElectroPuls E10000 с использованием температурной камеры Instron 3119. В испытаниях использовался высоконаполненный эластомер марки «ПДИ» (имитатор твердого топлива, НИИ полимерных материалов, г. Пермь) [27]. Вязкоупругие параметры материала определялись для установившегося режима деформирования (после выхода на заданный режим испытательной системы). Для реализации изотермического нагружения образца температурная камера работала в течение всего времени проведения экспериментов. Для уменьшения влияния разогрева образца на определяемые параметры время проведения эксперимента было ограничено по количеству циклов (в зависимости от вида воздействия) [27]. При определении вязкоупругих параметров материала учитывалась жесткость элементов экспериментальной установки вследствие того, что перемещения, фиксируемые используемым датчиком, складываются из перемещений образца и перемещений элементов этой системы [27]. Для подобных материалов характерна разносопротивляемость, поэтому перемещения элементов системы при сжатии образца больше, чем при его растяжении. С понижением температуры опыта и увеличением частоты нагружения перемещения элементов экспериментальной системы и, следовательно, искажение задаваемого сигнала становятся существеннее. Поэтому, при необходимости вносились корректировки в закон деформирования образца. Нужно также отметить, что геометрия образцов существенно влияет на определяемые вязкоупругие параметры, так как при нагружении в образце возникает сложное неоднородное напряженно-деформированное состояние. Для учета влияния геометрического фактора на определяемые параметры были решены динамические осесимметричные задачи в программном обеспечении ANSYS, которые соответствуют экспериментальной ситуации. В результате были определены поправочные коэффициенты [27]. Вязкоупругие параметры материала определялись с использованием преобразования Фурье аналогично работам [25-27]. Как уже отмечалось, для описания зависимости вязкоупругих параметров от частоты нагружения и температуры опыта были использованы полиномы (46)-(67). Для определения этих полиномиальных зависимостей были проведены экспериментальные исследования (таблица), результаты которых представлены на рис. 1-11. План экспериментальных исследований The experimental design Амплитуды деформации Частоты Температура εа1, % εа1, % ν1 ν2 Т1 Т2 Т3 Т4 Т5 3 3 0,1 1 50 30 3 -20 -50 0,01 0,1 Рис. 1. Зависимости параметров E*1/1, E*1/2 от частоты при температуре 3 °С Fig. 1. The dependences of parameters E*1/1, E*1/2 on frequency at the temperature of (T = 3 °С) Рис. 2. Зависимости параметров φE1/1, φE1/2 от частоты при температуре 3 °С Fig. 2. The dependences of parameters φE1/1, φE1/2 on frequency at the temperature of (T = 3 °С) Рис. 3. Зависимости параметров E*2/1, E*2/2 от частоты при температуре 3 °С Fig. 3. The dependences of parameters E*2/1, E*2/2 on frequency at the temperature of (T = 3 °С) Рис. 4. Зависимости параметра E*2/3 от частоты ν2 при температуре 3 °С Fig. 4. The dependences of parameter E*2/3 on frequency at the temperature of (T = 3 °С) Рис. 5. Зависимости параметров φE2/1, φE2/2 от частоты при температуре 3 °С Fig. 5. The dependences of parameters φE2/1, φE2/2 on frequency at the temperature of (T = 3 °С) Рис. 6. Зависимости параметра φE2/3 от частоты ν2 при температуре 3 °С Fig. 6. The dependences of parameter φE2/3 on frequency at the temperature of (T = 3 °С) Рис. 7. Зависимости параметров E*3/1, E*3/2 от частоты при температуре 3 °С Fig. 7. The dependences of parameters E*3/1, E*3/2 on frequency at the temperature of (T = 3 °С) Рис. 8. Зависимости параметров E*3/3, E*3/4 от частоты ν2 при температуре 3 оС Fig. 8. The dependences of parameters E*3/3, E*3/4 on frequency at the temperature of (T = 3 оС) Рис. 9. Зависимости параметров φE3/1, φE3/2 от частоты при температуре 3 оС Fig. 9. The dependences of parameters φE3/1, φE3/2 on frequency at the temperature of (T = 3 оС) Рис. 10. Зависимости параметров φE3/3, φE3/4 от частоты ν2 при температуре 3 оС Fig. 10. The dependences of parameters φE3/3, φE3/4 on frequency at the temperature of (T = 3 оС) Рис. 11. Зависимость lgaT от температуры Fig. 11. The dependence of lgaT on the temperature Заключение Рассмотрена математическая модель, описывающая линейное поведение вязкоупругих материалов в форме интегральных операторов, а также математическая модель, описывающая нелинейное поведение вязкоупругих материалов в виде интегрального ряда Вольтерра-Фреше. На основе нелинейного представления Вольтерра-Фреше предложена модель для описания нелинейного поведения материала в условиях действия стационарных одночастотных и двухчастотных нагрузок с использованием комплексных параметров. Проведен анализ математической модели для описания зависимостей вязкоупругих параметров от частот нагружения и температуры, предложены полиномиальные зависимости с использованием температурно-временной аналогии (суперпозиции, эквивалентности). Проведены экспериментальные исследования, после обработки которых были выявлены зависимости вязкоупругих параметров от частот нагружения и температуры. Полученные результаты позволяют разработать оптимальный план экспериментальных исследований (с использованием аппарата теории планирования экспериментов), с помощью которого с достаточной точностью будут определены коэффициенты предложенной модели. Далее возможна оценка влияния различных вязкоупругих параметров на точность описания поведения материала в условиях действия гармонических нагрузок.

Об авторах

А С Янкин

Пермский национальный исследовательский политехнический университет

Р В Бульбович

Пермский национальный исследовательский политехнический университет

С В Словиков

Пермский национальный исследовательский политехнический университет

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

Статистика

Просмотры

Аннотация - 168

PDF (Russian) - 64