Full Text

Введение Вязкоупругие полимерные композиты широко используются в современной авиационной и ракетной технике. Данные материалы должны выдерживать не только статические, но и высокие динамические нагрузки, обладать хорошими демпфирующими свойствами. Статья посвящена изучению влияния сложных динамических (бигармонических) нагрузок на поведение вязкоупругих наполненных полимеров. Для описания поведения вязкоупругого материала используются интегральные и дифференциальные операторы, а также комплексные модули [1-3]. В условиях действия стационарных гармонических колебаний для описания деформационных свойств вязкоупругих материалов удобно использовать комплексные модули [1-5]. Поведение материала в таком случае определяется динамическим модулем E* и углом потерь φE. При действии бигармонических нагрузок количество компонент увеличивается до четырех: и (индексы 1 и 2 соответствуют первой и второй гармоникам). Механические характеристики полимеров обычно существенно зависят от температуры проведения опыта и частоты нагружения [1]. Кроме того, для наполненных полимеров в условиях действия одночастотных нагрузок также характерны некоторые нелинейные эффекты, например отличие петли гистерезиса от эллиптической формы, размягчение материала (эффект Маллинза-Патрикеева), влияние амплитуды деформации εа (эффект Пэйна) и предварительной статической деформации εst на динамические модули [6-13]. При двухчастотных нагрузках также наблюдаются существенные зависимости динамических модулей от амплитуд деформации первой и второй гармоник нагружения [11, 14, 15]. Для изучения влияния частот бигармонического нагружения на поведение вязкоупругого материала было предложено ввести относительный частотный фактор ν1/ν2 [7, 14]. Тогда отношение ν1/ν2 имеет два предельных значения: ν1/ν2 = 0 (частота первой гармоники ν1 равна 0, амплитуда деформации первой гармоники преобразуется в предварительную статическую деформацию εst) и ν1/ν2 = 1 (частоты первой и второй гармоник равны, ν1 = ν2). И при ν1/ν2 = 0, и при ν1/ν2 = 1 двухчастотная нагрузка преобразуется в одночастотную. При подобном описании вязкоупругих свойств материала возникает ряд вопросов. Как осуществить переход от одночастотного нагружения к двухчастотному нагружению при граничных значениях ν1/ν2? Как описать поведение материала при близком значении частот (например, ν1 = 1 Гц, ν2 = 0,9 Гц)? Как определить динамические свойства при нецелом отношении частот гармоник? 1. Описание поведения материала при двухчастотных нагрузках Соотношение напряжений и деформаций при стационарных гармонических колебаниях можно записать в следующем виде: (1) ; где Е¢ - действительная часть комплексного модуля; Е²- мнимая часть комплексного модуля; εа - амплитуда деформации; σа - амплитуда напряжения, t - время; ν - частота. Обобщая соотношение (1) для многочастотного нагружения, запишем (2) где m - количество гармоник воздействия; φg - угол сдвига начальных фаз. Формулы (1), (2) подходят для описания линейного вязкоупругого поведения материала. При нелинейном поведении (зависимость динамического модуля от предварительной статической деформации, амплитуды деформации) использование уравнения (2) приводит к отличиям реального поведения материала от расчетного. Динамические нагрузки, действующие на фоне квазистатических, часто не приводят к большим деформациям, поэтому существенных искажений формы петли гистерезиса от эллиптической не происходит [9]. Исходя из этого возможен прием гармонической линеаризации [16] динамического нагружения в окрестности некоторого квазистатического нагружения. Экспериментальные исследования показывают, что при стационарных гармонических нагрузках в окрестности каких-либо предварительных статических деформаций для данного типа материалов наблюдаются различные значения динамического модуля E* [7, 8, 11]. Зависимость динамического модуля E* от предварительной статической деформации εst могут быть описаны полиномом второй степени: (3) Например, для режима нагружения с частотой 1 Гц и амплитудой деформации 1 % коэффициенты А1 = 0,155; А2 = -0,45; А3 = 11,45. Предположим, что предварительная статическая деформация εst изменяется во времени периодически. Тогда будут наблюдаться изменения динамического модуля E*(t) с периодом 1/ν. В случае одночастотного (моногармонического) нагружения предварительная статическая деформация εst играет роль первой (низкочастотной) гармоники при бигармоническом нагружении: (4) Тогда при бигармонических воздействиях значение в зависимости от деформации будет изменяться подобным образом, как в функции (3). Подставив (4) в (3), получим (5) Примем Между напряжениями и деформациями существует запаздывание на величину φE, поэтому запишем уравнение (5) в виде (6) Изменениями угла потерь в зависимости от предварительной статической деформации εst при одночастотных воздействиях пренебрегаем [11], следовательно, можно записать Для высоконаполненных полимерных сред, как правило, характерно различное сопротивление при растяжении и сжатии [17]. Исходя из этого в условиях симметричного цикла деформирования (рис. 1) возможны различные значения динамического модуля при растяжении и сжатии (σ1 ≠ - σ2, ε1 = -ε2). Рис. 1. Один цикл одночастотного нагружения Амплитуды деформации и напряжения, а также среднее значение напряжения определяются с помощью формул Далее определим две составляющие динамического модуля: Соответственно, динамический модуль при моногармонических нагрузках запишем в виде при бигармонических нагрузках Таким образом поведение материала при двухчастотных нагружениях можно описать семью составляющими: , , , , , , ; при одночастотных нагружениях - тремя составляющими: E*(0), E*(1), φE. 2. Методические вопросы проведения динамического эксперимента Некоторые методические аспекты были затронуты в работе [18]. Экспериментальные исследования проводились на электродинамической испытательной системе Instron ElectroPuls E10000 в Центре экспериментальной механики Пермского национального исследовательского политехнического университета [19]. Для определения динамических деформационных свойств использовался низкомодульный вязкоупругий материал марки ПДИ, из которого были изготовлены образцы круглого сечения диаметром 36,5±0,1 мм и высотой 39,5±1,5 мм (рис. 2). Всего было изготовлено 8 образцов. Данный материал обладает ярко выраженными вязкоупругими свойствами. а б Рис. 2. Чертеж образца (а) и металлической оснастки (б), выполненных в виде «грибка» Образцы подвергались термостатированию в температурной камере Instron 3119 в течение 3 часов с целью установления равномерной температуры по всему объему образца. На изготовленных образцах проводились одноосные динамические испытания по закону деформирования вида (7) После проведения испытаний для расчета брались данные установившегося режима деформирования, так как экспериментальной системе необходимо время для выхода сигнала на заданные значения деформации tзад. Кроме того, на начальном этапе деформирования наблюдается размягчение материала (эффект Маллинза-Патрикеева). Обычно эффект ограничен несколькими первыми циклами нагружения, после которых им можно пренебречь. Для примера на рис. 3 и 4 приведены зависимости напряжений от времени для одночастотного режима нагружения с частотой 1 Гц, амплитудой деформации 1 % и температурой 30 °С и для двухчастотного режима нагружения с частотами 1/7 и 1 Гц, амплитудами деформации 3 и 1 % и температурой 30 °С. Было отмечено, что для одночастотной нагрузки эффектом размягчения можно пренебречь после 5-го цикла (максимальные отличия между определяемыми вязкоупругими характеристиками для 5-го и 6-го циклов не превысили 1 %), для двухчастотной нагрузки - после 3-го цикла низкочастотной гармоники (максимальные отличия между определяемыми вязкоупругими характеристиками для 3-го и 4-го циклов не превысили 2 %). Рис. 3. Зависимость напряжения от времени для одночастотного нагружения Рис. 4. Зависимость напряжения от времени для двухчастотного нагружения Нужно отметить, что за каждый цикл нагружения материал выделяет тепло, которое идет на нагрев образца. Повышение температуры образца, в свою очередь, способствует изменению значений динамического модуля и угла потерь. Исходя из этого температурная камера должна работать в процессе нагружения. В результате тепло уносится с поверхности образца. Помимо этого, также было ограничено время проведения эксперимента (по количеству циклов). Материал почти полностью восстанавливает свои механические свойства в течение 24 часов. Далее возможно использовать образец снова, без существенных отличий в результатах. На рис. 5 представлены диаграммы нагружения образца до и после 24 часов выдержки для режима нагружения с параметрами ν1 = 1/11 Гц; ν2 = 1 Гц; Максимальное отличие данных составило 2,8 %. Рис. 5. Зависимость напряжения от времени для двухчастотного нагружения По умолчанию для определения перемещений система ElectroPuls использует встроенный датчик перемещений траверсы. Это значит, что фиксируемое данным датчиком перемещение складывается из реального перемещения на образце uис и перемещения системы uс. Исходя из этого перед проведением экспериментальных исследований был произведен замер перемещений элементов нагружающей системы (захватных приспособлений). Для этого в захватных приспособлениях испытательной системы ElectroPuls E10000 была закреплена стальная пластина толщиной 12 мм с минимальным зазором между захватными приспособлениями. Затем были проведены испытания на растяжение и сжатие. Поскольку расстояние между захватами мало по сравнению с их длиной, удлинением пластины можно пренебречь. Соответственно, перемещения, зафиксированные по встроенному датчику системы, и будут перемещениями элементов нагружающей системы. Зависимость перемещения элементов системы от приложенной нагрузки можно описать функцией вида Для нахождения «истинных» перемещений образца нужно из перемещений, фиксируемых испытательной системой uэк, вычесть перемещения испытательной системы uс при соответствующих нагрузках. 3. Методика определения механических свойств вязкоупругих материалов, представление их в виде комплексных модулей После проведения двухчастотных испытаний результаты эксперимента (напряжения, деформации и т.д.) сохраняются в виде электронных таблиц (табл. 1). Таблица 1 Результаты динамического эксперимента Время ti, с Деформации (εэк)i, % Напряжения σi, МПа t1 (εэк)1 σ1 t2 (εэк)2 σ2 … … … tn (εэк)n σn Для учета жесткости нагружающей системы необходимо определить перемещения, фиксируемые испытательной системой ((εэк)i - деформации из табл. 1); нагрузку перемещения элементов испытательной системы «истинные» перемещения «истинные» деформации Затем построим табл. 2 аналогично табл. 1. Таблица 2 Результаты динамического эксперимента Время ti, с Деформации εi, % Напряжения σi, МПа t1 ε1 σ1 t2 ε2 σ2 … … … tn εn σn Между соседними точками значений напряжений и деформаций (например, ε2 и ε3), взятых из табл. 2, проведем отрезки. Совокупность этих отрезков можно описать функциями вида Далее данные эксперимента разбиваются на временные интервалы через каждый период колебаний первой низкочастотной гармоники 1/ν1, и выделяется первый низкочастотный временной интервал В результате такого разбиения получим n интервалов. Каждый низкочастотный временной интервал разложим в ряд Фурье [20] с помощью формул (8) Далее с помощью формул получим первую низкочастотную гармонику по деформациям (9) После аналогичных преобразований для напряжений выделим первую низкочастотную гармонику по напряжениям (10) Аппаратная реализация закона нагружения по деформации в общем случае предполагает наличие угла сдвига фазы вследствие запаздывания экспериментального оборудования. Предложенная методика обработки позволяет учитывать запаздывание экспериментального оборудования. Значение при симметричном цикле деформирования будет стремиться к нулю. Таким образом, для первой низкочастотной гармоники можно определить значения составляющих динамического модуля и угла потерь: С помощью аналогичных преобразований можно определить значения составляющих динамического модуля и угла потерь при одночастотных нагрузках. Подставим теперь в выражения (9) и (10) вместо t значения времени ti из табл. 2 и вычтем выражения (9) и (10) из табличных значений напряжений и деформаций соответственно В результате предложенных преобразований можно получить значения напряжений и деформаций только от высокочастотной гармоники. Данные значения разбиваются на высокочастотные временные интервалы периодом 1/ν2 следующим способом: - 1-й высокочастотный временной интервал, - 2-й, …, - последний и т.д. до тех пор, пока j + 1 - количество высокочастотных временных интервалов). Затем с помощью аналогичных преобразований (8), описанных выше, можно определить зависимости для каждого высокочастотного временного интервала: Для каждого высокочастотного временного интервала определяются значения составляющих динамического модуля и угла потерь для высокочастотной гармоники (табл. 3). Таблица 3 Значения динамического модуля и угла потерь для второй гармоники Номер временного интервала Время с Высокочастотный временной интервал Угол потерь град. Динамический модуль МПа 1 2 3 … … … … … j + 1 Так как значения динамического модуля и угла потерь определяются для конкретного интервала времени, то они присваиваются среднему значению времени этого интервала. Таким образом, в результате обработки получается j + 1 значение угла потерь и динамического модуля высокочастотной гармоники. Угол потерь следовательно, Динамический модуль второй гармоники характеризуется уравнением (6), поэтому для определения составляющих Е2*(0), Е2*(1), Е2*(2) воспользуемся преобразованиями Фурье: С помощью формул определяются значения составляющих динамического модуля. Затем выделяется второй низкочастотный интервал, третий и т.д. Всего в итоге получаем n значений составляющих динамического модуля и угла потерь для первой низкочастотной и второй высокочастотной гармоник. Далее определяются их средние значения. Геометрия используемых образцов существенно влияет на определяемые вязкоупругие характеристики вследствие возникновения сложного неоднородного вида напряженно-деформированного состояния при нагружении (при растяжении или сжатии торцы образца, приклеенные к металлической оснастке, не могут свободно перемещаться в поперечном направлении). Поэтому для определения «истинных» вязкоупругих характеристик композита были решены соответствующие динамические осесимметричные задачи на растяжение-сжатие в программном обеспечении ANSYS. В расчетах принимались следующие характеристики материала грибка (сталь): плотность ρс = 7800 кг/м3; коэффициент Пуассона μс = 0,3; модуль упругости Ес = 200 ГПа. Для полимерного композита была принята плотность ρк = 1720 кг/м3. Постоянство коэффициента Пуассона наполненных полимеров часто не подтверждается экспериментально [6, 10, 21, 22]. Исходя из этого, помимо варьирования «истинных» значений динамического модуля при определении напряженно-деформированного состояния также варьировались «истинные» значения коэффициента Пуассона (табл. 4). На основе анализа напряженно-деформированного состояния определялись значения динамического модуля, соответствующие экспериментальной ситуации его (динамического модуля) определения где F - площадь поперечного сечения образца; N - нагрузка, возникающая при деформировании образца. При сопоставлении «истинного» и экспериментального значений динамического модуля определялся коэффициент корректировки (см. табл. 4). Таблица 4 Значения коэффициентов корректировки динамического модуля № п/п Динамический модуль, МПа Коэффициент Пуассона Коэффициент корректировки KE 1 1 1,21 0,495 0,83 2 10 11,91 0,84 3 100 112,15 0,89 4 1 1,02 0,4 0,98 5 10 10,12 0,99 6 100 95,31 1,05 7 1 0,93 0,3 1,08 8 10 9,23 1,08 9 100 86,97 1,15 Данные из табл. 4 показывают, что корректировочный коэффициент KE зависит от значений динамического модуля и коэффициента Пуассона. Данная зависимость с высокой точностью описана функцией вида С использованием зависимости определялись «истинные» значения составляющих динамического модуля. 4. Особенности обработки экспериментальных данных при граничных значениях отношения частот бигармонического нагружения При отношении частот ν1/ν2 = 0 частота первой гармоники ν1 равна 0, амплитуда деформации первой гармоники преобразуется в предварительную статическую деформацию εst. С помощью уравнений (3)-(6) можно определить значения составляющих динамического модуля и угла потерь в бигармонической постановке задачи. При отношении частот ν1/ν2 = 1 частоты первой и второй гармоник равны, ν1 = ν2 (одночастотное нагружение), при этом принимаем, что Одночастотную нагрузку можно представить в виде двухчастотной следующим образом: Амплитуда деформации εа(n) зависит от угла сдвига начальных фаз Рассмотрим два предельных значения: а) примем тогда С увеличением амплитуды εа значения динамического модуля E* уменьшаются: следовательно, тогда можно принять Далее принимаем для первой гармоники Тогда для второй гармоники б) примем тогда После аналогичных преобразований получаем Эксперимент показывает, что при примерно равен при Таким образом, не оказывает существенного влияния на динамические модули и Данный подход позволяет определить значения динамического модуля и угла потерь в бигармонической постановке задачи при отношении ν1/ν2 = 0 и ν1/ν2 = 1. 5. Экспериментальные исследования поведения вязкоупругого материала при различном соотношении частот бигармонического нагржения Проведены экспериментальные исследования поведения вязкоупругого полимерного композита марки ПДИ при различном соотношении частот (табл. 5) бигармонического нагружения (7) при температуре 30 °С. С помощью описанной выше методики определялись значения составляющих динамических модулей и углов потерь (табл. 6, рис. 6, 7). Номера режимов нагружения из табл. 6 соответствуют номерам режимов нагружения из табл. 5. Таблица 5 Параметры нагружения образцов Номер режима нагружения Частота гармоник, Гц Амплитуда деформации гармоник, % ν1/ν2 ν2 1 0 1 3 1 2 0,01 3 0,02 4 0,05 5 1/11 6 0,1 7 1/8 8 1/7 9 1/6 10 0,2 11 1/4 12 1/3 13 0,5 14 1 Таблица 6 Значения динамических деформационных свойств материала Номер режима нагружения МПа МПа МПа град. МПа МПа град 1 12,15 -1,35 0,70 17,87 - - - 2 11,82 -1,35 0,73 18,28 4,44 -0,43 17,72 3 11,68 -1,29 0,79 18,50 5,39 -0,46 18,38 4 11,38 -1,32 0,80 18,69 5,91 -0,49 18,31 5 11,09 -1,35 0,88 19,14 6,53 -0,50 18,41 6 11,57 -1,40 0,87 18,9 6,93 -0,53 18,6 7 11,10 -1,40 0,87 19,24 6,95 -0,49 18,60 8 11,10 -1,31 0,83 19,72 7,11 -0,53 18,82 9 10,88 -1,32 0,90 19,46 7,22 -0,51 18,51 10 11,38 -1,49 0,93 19,83 7,62 -0,55 18,90 11 10,84 -1,38 0,81 19,97 7,72 -0,53 18,61 12 10,89 -1,33 0,62 19,72 8,21 -0,61 18,81 13 10,33 -1,20 0,34 19,23 9,08 -0,68 18,74 14 9,20 -1,16 0 17,9 10,69 -0,74 17,9 Также были проведены бигармонические испытания, аналогичные испытаниям, результаты которых отражены в табл. 5, при ν2 равной 0,1 Гц и 10 Гц. Результаты исследования представлены на рис. 7. а б в Рис. 6. Зависимости компонент (а); (б) и (в) от логарифма частоты lg(ν1) а б в г Рис. 7. Зависимости компонент (а); (б); (в) и (г) от логарифма отношения частот lg(ν1/ν2) при различных значениях частоты ν2 (* экспериментальные данные при lg(ν1/ν2) = = -3 соответствуют одночастотному режиму нагружения с варьированием предварительной статической деформации (табл. 5, режим нагружения 1)) Анализ экспериментальных данных показывает, что отношение частот ν1/ν2 практически не оказывает влияния на компоненты , при ν1/ν2 ≤ 0,2 и при ν1/ν2 ≤ 0,02. Максимальное отличие значений угла потерь φE2 при варьировании lg(ν1/ν2) составили 10 % при ν2 = 1 Гц, 7 % при ν2 = 0,1 Гц, 7 % при ν2 = 10 Гц. Можно принять Статическая деформация практически не оказывает влияния на угол потерь, следовательно, Частота ν2 также существенно влияет на составляющие , , , . Данные зависимости можно описать полиномами второй степени: (11) (12) (13) (14) Также было установлено, что частота ν2 не оказывает существенного влияния на компоненты низкочастотной гармоники а частота ν1 существенно влияет на (15) (16) (17) Зависимости, приведенные выше, характерны только для исследованного диапазона частот. При расширении диапазона частот (например, с использованием принципа температурно-временной аналогии) данные зависимости могут стать более сложными. Следует отметить, что статические испытания (см. табл. 5, режим нагружения № 1) сложны для экспериментатора. Для того чтобы определить зависимость динамического модуля от статической деформации, необходимо провести как минимум три испытания. Материал почти полностью восстанавливает свои свойства через 24 часа поле проведения испытаний, поэтому нельзя использовать один образец несколько раз в течение дня. Кроме того, на величины динамического модуля также влияет скорость выхода на заданный уровень статической деформации. С другой стороны, достаточно провести всего один двухчастотный эксперимент, чтобы определить зависимость динамического модуля от статической деформации. Низкочастотная гармоника в данном случае - аналог статической деформации, а высокочастотная гармоника - аналог моногармонической нагрузки. Скорость выхода на заданный уровень статической деформации будет определяться амплитудой и частотой низкочастотной гармоники. Таким образом, бигармонические испытания являются более предпочтительными для определения зависимостей динамического модуля от предварительной статической деформации. Представляет интерес также проведение испытаний при нецелом отношении частот ν1/ν2 с последующим сравнением результатов эксперимента (σэк) с результатами, рассчитанными по предложенным зависимостям (11)-(17) (σм). Для этого были проведены два эксперимента: при ν1/ν2 = 1/4,25 (рис. 8) и при ν1/ν2 = 0,9 (рис. 9) (когда значения частот гармоник близки друг к другу, наблюдается периодическое увеличение и уменьшение амплитуды - биение), остальные параметры нагружения соответствуют параметрам из табл. 5 при температуре 30 °С. Рис. 8. Зависимость напряжения от времени проведения испытания при отношении частот ν1/ν2 = 1/4,25 При ν1/ν2 = 1/4,25 максимальные отличия между σэк и σм составляют 5,9 %. При ν1/ν2 = 0,9 максимальные отличия между σэк и σм составляют 6,8 %. Помимо этого, были проведены экспериментальные исследования при различных значениях угла сдвига начальных фаз φg2. Анализ результатов показал, что изменение φg2 практически не оказывает влияния на определяемые характеристики (составляющие динамических модулей и углы потерь). Рис. 9. Зависимость напряжения от времени проведения испытания при отношении частот ν1/ν2 = 0,9 Заключение В ходе работы затронуты некоторые методические вопросы проведения экспериментальных исследований при стационарных гармонических нагрузках, описана методика определения динамических деформационных свойств (в виде комплексных модулей) вязкоупругих материалов при двухчастотных нагрузках, проведены экспериментальные исследования поведения низкомодульного вязкоупругого полимерного композита при бигармонических (двухчастотных) нагружениях и моногармонических (одночастотных) нагружениях. В рамках одночастотного нагружения реализовывались режимы с изменением предварительной статической деформации и угла сдвига начальных фаз. В рамках бигармонического нагружения реализовывались режимы с изменением частот первой (низкочастотной) и второй (высокочастотной) гармоник. По результатам испытаний в соответствии с предложенной моделью описания вязкоупругого поведения материала определялись составляющие динамических модулей и углы потерь. Предложена методика представления результатов эксперимента при граничных значениях отношения частот первой и второй гармоник ν1/ν2 = 0 и ν1/ν2 = 1, а также при любом нецелом отношении частот. Данный подход позволяет определить значения динамического модуля и угла потерь в бигармонической постановке задачи. Анализ результатов проведенных экспериментальных исследований показал, что изменение отношения частот ν1/ν2 практически не оказывает влияния на компоненты при ν1/ν2 ≤ 0,2; при ν1/ν2 ≤ 0,02 и на угол потерь Изменение частоты ν2 существенно влияет на , но не оказывает существенного влияния на компоненты низкочастотной гармоники Частота ν1 существенно влияет на Были предложены зависимости, описывавшие поведение материала при различных значениях частот первой и второй гармоник. Показано, что бигармонические испытания являются более предпочтительными для определения зависимостей динамического модуля от предварительной статической деформации, чем одночастотные испытания при варьировании предварительной статической деформации. Угол сдвига начальных фаз практически не оказывает влияния на поведение вязкоупругого материала. В дальнейшем необходимо проведение испытаний в широком диапазоне температур, разработка многофакторной математической модели, описывающей зависимости вязкоупругих характеристик (составляющих динамических модулей и углов потерь) от различных параметров (амплитуды деформации и частоты) и условий (температура) нагружения, определение коэффициентов этой модели, дисперсии повторяемости результатов эксперимента, а также проверка адекватности модели.

About the authors

A S Yankin

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 166

PDF (Russian) - 55