CC BY
35
ДИНАМИКА, ПРОЧНОСТЬ, НАДЕЖНОСТЬ ¡
УДК 539.3:539.377
К. Ф. Комков
ОСОБЕННОСТИ УПРУГИХ СВОЙСТВ ВЫСОКОНАПОЛНЕННЫХ ПОЛИМЕРНЫХ МАТЕРИАЛОВ
Приведены результаты экспериментальных исследований низкомодульных высоконаполненных полимеров. Выполнен анализ, позволяющий представить механизм деформации сред, механические свойства которых имеют некоторые особенности: аномально высокие значения коэффициента поперечной деформации, разномодуль-ность и очень малое значение отношения модуля сдвига к модулю объемного сжатия. Для описания деформации высоконаполненных полимерных материалов использованы тензорно-нелинейные уравнения. Результаты расчетов по предложенной модели сопоставлены с результатами испытаний образцов.
Известно [1], что высоконаполненные полимерные материалы (ВНП) проявляют зависимость свойств от времени, поэтому под характеристиками упругости следует понимать значения мгновенных модулей или податливостей. Для проведения расчета напряжено-деформированного состояния изделий из ВНП, например при нагрузках, возрастающих с большой скоростью, наличие упругих характеристик является одним из главных условий. Однако при обработке результатов испытаний возникают трудности в установлении связи между модулями упругости из-за отсутствия надежных сведений о значении коэффициента поперечной деформации.
Поэтому их приходится определять по результатам лабораторных испытаний для конкретных видов напряженного состояния. Например, для определения модуля объемного сжатия разработана и изготовлена специальная установка, несколько видоизмененная схема которой приведена в работе [2]. Нагружение осуществлялось в камере, заполняемой жидкостью или сжатым воздухом. Изменение объемной деформации фиксировалось визуально. Результат испытания — график зависимости объемной деформации от давления, меняющегося в пределах 0... 5 МПа. Искомый модуль K определяли по начальному участку графика 1, изображенного на рис. 1.
Модуль упругости при растяжении и сжатии определяли по графикам напряжение-деформация (as), полученным при испытаниях образцов. На рис. 2, а приведена схема приспособления, которая включает в себя образец 1, датчик продольных деформаций 3, датчик
AF/F-10
Е(р)1Е( 0)
1
2
1,2
1,0
силы 2 и систему шарниров для исключения изгиба образца. Образец с тензометрическими датчиками продольной и поперечной деформации показан на рис. 2, б.
Схема регистрации, включающая в себя датчики с усилителем и двухкоординатными потенциометрами, позволила записать диаграммы а-е и графики зависимости поперечной деформации от продольной (е2-е1). Подобная же схема использовалась для испытаний на сжатие. По этой схеме проводились испытания на растяжение в камере при давлении 0... 5 МПа, что позволило определить зависимость модуля упругости Е от давления (см. рис. 1, кривая 1).
Графики е2-е1 дали возможность установить, что при малых деформациях растяжения (е1 < 3... 4 %) коэффициент поперечной деформации ^12 может принимать значения большие 0,5. На
1
Рис. 1. Влияние внешнего давления на свойства ВНП:
1 — график изменения модуля упругости при растяжении от гидростатического сжатия; 2 — график зависимости объемной деформации от давления
Рис. 2. Схема приспособления для испытания образцов ВНП на растяжение:
1 — образец, 2 — датчик силы, 3 — датчик для измерения продольной деформации; 4 — усилитель, 5 — двухкоординатный потенциометр
рис. 3 приведен график (кривая 1), усредненный по результатам испытаний пяти образцов, вырезанных из заготовки в различных направлениях. Для повышения достоверности результатов тарировка датчиков дополнялась испытаниями тестовых материалов.
График 3 относится к образцам из плотной, наполненной сажей, резины. Для нее (по литературным данным) ^12 = 0,48 ... 0,49. График 2 относится к образцам из пластичного материала, не подвергавшегося полимеризации. Для них, как и ожидалось, ^12 = 0,5. В результате испытаний образцов из других ВНП, отличающихся составом связующего, получены кривые е2-еь расположенные над тестовыми кривыми 2 и 3.
Поскольку отмеченные особенности
Рис. 3. Зависимости поперечной деформации от продольной:
1 — при испытании образцов ВНП, 2 — образцов не полиме-ризованного ВНП, 3 — образцов резины
упругих свойств ВНП проявляются при деформациях, существенно меньших, чем при разрушении, то на рис.4 приведены начальные участки диаграмм а-е (кривые 1 и 2) и графиков и-е (кривые 3 и 4) для растяжения и сжатия.
Полученные результаты стимулировали поиск сведений о подобных исследованиях таких или близких по механическим свойствам, материалов. Выяснилось, что аномально высокие значения ^12 и описаны в работе [3] для условий растяжения образцов из полиэтилена и фторопласта. Там же высказано предположение, что этот эффект связан с деформационной анизотропией и изменением объема.
Для изотропного материала
условие V > 0,5 равнозначно утверждению, что деформация растяжения сопровождается уменьшением объема образца. Чтобы убедиться в этом, провели испытания с измерением объемной деформации посредством специально изготовленного дилатометра. Эксперименты на образцах из ВНП не обнаружили отрицательной объемной деформации. Уровень жидкости в измерительной трубке оставался неизменным до
(т/Я-Ю2
5.0 3,75 2.5 1,25
3
(Л —- \\
■■■ 4 2\
5 N
V
-0.5 0,48 -0,46
-6 -4 -2 0 2 4 6 £-ю2
Рис. 4. Начальные участки диаграмм а—е (кривые 1 и 2) и графиков V—е (кривые 3 и 4):
справа — при растяжении, слева — при сжатии; сплошные кривые — эксперимент, штрих-пунктирные — расчет (кривые 5 и 6)
е = 0,6... 0,8 %, после е > 2 % с возрастающей скоростью уровень жидкости поднимался (подобные результаты можно найти в обзорной работе [1]).
Полученные сведения дали основание для поиска механизма деформации, способствующего эффекту повышенной поперечной деформации. Основной вклад в описанную гипотезу дала работа [4], в которой проведено исследование, раскрывающее причины зависимости модуля упругости от внешнего давления. В этой работе при растяжении образцов ненаполненных полимеров исследовалось состояние макроструктуры с помощью микроскопа. Обнаружено, что при повышенной температуре, когда деформация достигает & 1 %, в поле зрения микроскопа появляются полосы сдвига под углом примерно 45° к оси образца, сопутствующие существенному повышению податливости. Такое поведение материала названо "сдвиговым деформированием". Эти полосы, незаметные для невооруженного глаза, есть следствие структурных изменений. Они аналогичны линиям скольжения, возникающим в образцах пластичных материалов при напряжениях, превышающих предел текучести [5].
Предполагается, что нечто подобное происходит при сравнительно малой деформации низкомодульных ВНП. При высокой степени наполнения полимера минеральными частицами (рис.5,а) малую жесткость можно получить при слабом межмолекулярном взаимодействии в клеевом слое и высокой подвижности их цепочек на поверхности наполнителя [6].
Если за поле наблюдения принять лунку из-под крупной частицы, то при большом увеличении наблюдается примерно такая же структура, что и на рис. 5, а. Пограничные цепочки (связи) могут скользить (позиция 1, рис.5,б) или отслаиваться (позиции 2 и 3) при сохране-
Рис. 5. Поверхность разрушения ВНП (х30):
а — элементы структуры, 1 — зерно наполнителя, 2 — лунка из-под зерна; б — взаимодействие цепочек молекул полимера (связей) с поверхностью наполнителя, 1 — микродеформация со скольжением, 2 и 3 — с отслоением
а
б
3
нии общего числа точек контакта за счет энергии теплового движения молекул.
Исходной причиной особенностей упругого поведения ВНП является существенная неоднородность поля деформаций из-за наличия в нем твердых частиц наполнителя. Макродеформация материала включает в себя существенные микродеформации и перемещения связей на поверхности наполнителя и молекул клеевого слоя, имеющего большую склонность к изменению формы, чем объема.
Если к тому же сопротивление связей на отрыв больше сопротивления скольжению, то отмеченное способствует приобретению регулярности в структуре связей и появлению сдвигов между поверхностью наполнителя и клеевого слоя. Совокупность таких сдвигов при их слиянии образуют "полосы", не обязательно прямолинейные, в направлении касательных напряжений. На макроуровне этот процесс можно представить в форме "деформирование сдвигом".
Несмотря на то, что видимых признаков наличия полос сдвига на поверхности образца не обнаруживается из-за сложности макроструктуры, его поведение в процессе активной деформации можно сравнить с поведением плоского образца, испытывающего развитую пластическую деформацию [7, 8], что схематично представлено на рис. 6, б.
При разгрузке заметных остаточных деформаций не наблюдается, и остаточные напряжения при таком состоянии связующего тоже следует считать пренебрежимо малыми. Структура материала при сравнительно малой деформации изменяется и он приобретает признаки анизотропного тела, упругие свойства которого зависят от направления и знака главных напряжений.
Дальнейший рост макродеформации ведет к массовому отрыву связей, увеличению объемной деформации (дилатансии), изменению возникшей регулярности, с которой связаны отмеченные эффекты. Проявляется уже другой, более заметный механизм, ведущий не только к повышению податливости, но и к разрушению материала [9]. Более
Рис. 6. Диаграмма интенсивность напряжений-интенсивность деформаций (а) и схема образца при дефомировании сдвигом (б):
1 — растяжение, 2 — сжатие
глубокое изучение описанного механизма деформации ВНП требует целенаправленных исследований.
Для оценки состояния среды с учетом отмеченных эффектов используются тензорно-нелинейные уравнения, которые для главных деформаций приводятся к виду [10]
ег = а,а, {г,з = 1, 2, 3), (1)
где ог, £г — главные напряжения и деформации; а, — переменные коэффициенты деформаций [11], которые связаны с используемыми в технических приложениях характеристиках следующим образом:
агг = аг, = = 1, 2, 3, « = 3! ^ «,3 (2)
Здесь Ег — модули упругости в направлении главных напряжений; V, — коэффициенты поперечных деформаций. Их можно определить по соотношениям
Е = 9/(3Фт + Фк + Ф^СУ; (3)
V, = (3/2Фт - Фк - ФаСгз)Е/9, 1,3 = 1, 2, 3, (4)
где Фт, Ф^, Фк — обобщенные характеристики среды. Первая из них представляет собой среднее арифметическое значение податливостей в направлении главных касательных напряжений, вторая — пропорциональна среднему квадратичному отклонению этих податливостей, а третья характеризует податливость при объемном растяжении (сжатии). Последняя из них при малых деформациях ВНП значительно меньше Фт и Ф^. При оценочных расчетах можно принять Фк = 1/К (где К — модуль объемной деформации), Фк = 0:
Сгг = Ci(1 + ffi); (5)
= Ca + i, j, a = 1, 2, 3, i = j = a; ^ i, j, a; (6)
C1 = 2 cos 9, C2 = V^sin 9 — cos 9, C3 = — (v^3 sin 9 + cos 9); (7)
9 = 3 ™ , 0 < 9 < n/3; (8)
Sj = Gij — , S0 = (3/2SijSij)1/2, = Gii/3.
Здесь Gj, Sij — тензор и девиатор напряжений; S0, g0 — интенсивность напряжений и среднее напряжение соответственно; 9 — угол вида напряженного состояния; ж — параметр, характеризующий интенсивность снижения плотности связей и роста положительной объемной деформации (дилатансии).
Определяя интенсивность деформаций e0 = (2/3eij-e*j)1/2, уравнения (1) позволяют установить зависимость между инвариантами e0
и So
eo = So/3« + 4/3ФтФ^п + 4/9Ф2)1/2. (9)
В этих соотношениях ej = ej — eo5j; eo = £ü/3; п = cos3#; ej, ej — тензор и девиатор деформаций соответственно, eo — средняя деформация.
Аналогично можно найти параметры Лоде [10]
А, = — (2Ti/T2 + 1) = 2(C2 — Cs)/(Ci — Сз) — 1; (10)
Ае = — (271/72 + 1) = 2(А, + 1)(3Фт — 2Ф^С1)/(3Фт — 2Ф^С2) — 1, (11)
а также их разницу, характеризующую отклонение от подобия девиа-торов напряжений и деформаций
ДА = А, — А£ = 2ФДА, + 1)(С1 — С2)/(3Фт — 2Ф^С2), (12)
где тг = (aj — aa)/2, y% = е — £«, (i, j, a = 1, 2, 3; i = j = a) — главные касательные напряжения и главные деформации сдвига соответственно.
Принимая допущение, что Фт и Ф^ не зависят от вида напряженного состояния, из соотношений (9)-(12) можно определить податливости формоизменения и разницу параметров Лоде [12] при: растяжении (С1 = 2, С2 = —1, С3 = —1, п = 1)
Фр = 3eo/So = Фт + 2/3Ф^, ДА = 0; (13)
сжатии (С1 = 1, С2 = 1, С3 = —2, п = —1)
Фс = Фт — 2/3Ф^, ДА = 0; (14)
чистом сдвиге (С1 = v^3, С2 = 0, С3 = — v^3, П = 0)
Фт = (Фт + 4/9Ф^)1/2, ДА = 2Ф^/(У3Фт). (15)
Чтобы получить представление о Фр и Фс по экспериментальным данным, необходимо располагать диаграммами So—eo для растяжения и сжатия. При этом следует учесть, что диаграммы a-e (см. рис. 4) не являются независимыми от ao. Испытания образцов в камере показывают, что модуль упругости на начальном участке возрастает практически пропорционально давлению (см. рис. 1, кривая 1), т.е. E(p) = E(0)(1 + a1p), где E(0), E(p) — модули упругости до и после приложения внешнего давления p, a1 — коэффициент, определяемый по результатам испытаний.
Учитывая принятое допущение и тот факт, что при простом растяжении (сжатии) гидростатическое напряжение переменно, для коэффициента a1 принимается линейная зависимость a1 = a,ao. Это позволяет уточнить напряжения:
So = a(0) = a(ao)/[1 — a, (signaos/E2]. (16)
АЛ-
1,5 1,0 0.5
0.75
3\ """у
5 "7,5
0 0,48 ^2 -
0 1,25 2.5 3,75 (50/£>102
Рис. 7. Зависимости обобщенных податливостей Фт и Ф^ (кривые 1 и 2) и разницы параметров Лоде (кривая 3) от интенсивности напряжений
Интенсивность деформации вычисляется с учетом значений V на кривых ^-е, а именно е0 = 2(1 + V)е/3.
Для фиксированных значений е0 по графикам 1 и 2 (см. рис. 6, а) находят 50 для растяжения и сжатия и определяют соответствующие податливости Фр и Фс, которые позволяют найти
Фт = (Фр + Фс)/2, Ф^ = 3(Фр - Фс)/4. (17)
Поскольку последние приняты независимыми от 0, то, исходя из соотношений (15), можно вычислить ДА = /(50), Фт и 50т = = л/2Б0рБ0с/(50р + 50С)1/2. На рис. 7 приведены графики Фт, Ф^ и ДА в зависимости от 50/Е. Соответствующие кривые 1, 2 и 3 при 50 ^ 0 стремятся к исходному ненулевому значению и тем самым указывают на наличие разномодульности материала.
Описанная модель предполагает, что начальная анизотропия, связанная, например, с технологией изготовления материала или изделия из него, вполне возможна, но при широком диапазоне размеров округлых зерен наполнителя (0... 300 мк) она мала по сравнению с анизотропией, вернее с ее проявлением, при деформации. Последняя имеет то преимущество, что коэффициенты деформации могут быть определены по функциям Фт и Ф^, найденным согласно (17) по результатам испытаний на простое растяжение и сжатие. Для этого достаточно воспользоваться соотношениями (3) и (4).
Как уже отмечалось ранее, для мягких ВНП при малых деформациях допустимо принять Фк = 0. Тогда можно оценить начальные значения модулей упругости по соотношениям (3), а именно
Е = 3^/(1 + С/), (18)
и коэффициентов поперечной деформации по соотношениям (4):
^2 = ^з = (0,5 - Сз/)/(1 + С/) = (0,5 + /)/(1 + 2/) = 0,5 (19) — при растяжении;
V31 = V32 = (0,5 - С2/)/(1 + Сз/) = (0,5 - /)(1 - 2,1) = 0,5 (20)
— при сжатии;
^ = (0,5 + УЭ/)/(1 + ^3/) > 0,5; vзl = 0,5/(1 - ^3/) > 0,5 (21)
— при чистом сдвиге, где принято / = Ф^/3Фт, ^ ~ 1/Фт и ж = 0.
Из соотношений (19)-(21) следует, что на характеристики Ег и Vij влияют тригонометрические функции Сг, отражающие вид напряженного состояния. Результаты вычислений по соотношениям (19) и (20) расходятся с результатами испытаний, представленными на рис. 4 начальными точками графиков 3 и 4, тогда как расхождения с результатами вычислений по выражению (21) более приемлемы.
Если принять во внимание тот факт, что при деформировании сдвигом существенную роль играют главные касательные напряжения, то допустимо состояние чистого сдвига дополнять для каждого исходного с помощью величин Сг. Не нарушая соосности тензоров аг, и ег, и вводя в рассмотрение функцию
д = 1/(1+ в ехр(-а£)), (22)
где а и в — коэффициенты, а £ = ао/50, можно внести изменения в соотношения (19) и (20) с помощью следующего уравнения:
С* = Сг(1 - д) + Стгд, (23)
где Стг — состояние чистого сдвига.
Выполненное наложение следует рассматривать как дополнение к алгоритму, учитывающее наличие рассмотренного механизма деформации и позволяющее согласовать теорию с опытом. Результаты расчета начальных значений vгj с учетом уравнения (23) в виде табл. 1 дают представление о их различии. Среднее значение vс = 0,5 при всех напряженных состояниях. Коэффициенты а и в подбирали, исходя из условия лучшего согласования с опытными данными для начальных значений V при растяжении и сжатии.
Таблица 1
Коэффициенты поперечной деформации
Вид напряженного состояния V12 V21 V13 V31 V23 V32
Растяжение 0,517 0,554 0,481 0,536 0,445 0,463
Чистый сдвиг 0,520 0,542 0,479 0,522 0,457 0,478
Сжатие 0,532 0,537 0,467 0,504 0,463 0,495
Выше отмечалось, что снижение плотности внутренних связей учитывается переменной величиной ж. Как показано в работе [10], она может быть принята в следующем виде:
ж = Жи^и , (24)
где Би = $0/$ь; Бь и т — постоянные величины. Интенсивность этого процесса существенно зависит от вида напряженного состояния и гидростатического напряжения. Поэтому жи принята в виде функции, аналогичной (22):
Жи = Жо/[1 + в1 ехр(-а1 £)], (25)
где ж0, в1, а1 — постоянные величины. Их значения определяются из условия лучшей аппроксимации экспериментальных графиков ^е. На рис. 4 эти графики представлены кривыми 5 и 6, построенными по данным из табл. 2.
Таблица 2
£ • 102 1,5 2,5 3,5 4,0 4,5 5,0 5,5 6,0
Растяжение (v12) 0,511 0,505 0,496 0,490 0,484 0,478 0,472 0464
Сжатие (v31) 0,503 0,504 0,505 0,506 0,507 0,508 0,510 0,512
Необходимо принять во внимание, что учет рассмотренных эффектов по разработанной модели может внести поправки в расчет напряженно-деформированного состояния изделий из ВНП только при их малой деформации. Традиционные подходы, описывающие эту среду тензорно-линейными уравнениями, вполне оправданы.
Трудности с определением модулей упругости, как показано в работе [1], вполне преодолимы, если использовать не экспериментальные данные с прямым измерением поперечной деформации, а интегральные результаты, которые следуют из испытаний с определением объемной деформации. Однако при этом останутся скрытыми отмеченные эффекты и связанные с ними особенности упругих свойств. Опытные и расчетные данные о начальных значениях Е и V, показывают, что разномодульность ВНП реальна и ее следует учитывать при анализе упругих свойств.
В заключение следует отметить, что рассматриваемые особенности упругого поведения ВНП могут быть отражены в феноменологической модели, включающей в себя тензорно-нелинейные уравнения связи деформаций с напряжениями и дополнения (23) и (24). Необходимость подобных дополнений при описании поведения сложных сред подтверждают многие исследователи. Например, в работе [13] для учета внутренних сдвиговых процессов главные компоненты обобщенного
тензора напряжений включают в себя интенсивность напряжений. Это говорит еще и о том, что работы по данной проблеме нельзя назвать законченными.
СПИСОК ЛИТЕРАТУРЫ
1. S chwarzl F. R. On mechanical properties of unfilled and filled elastomers // Mechanics and chemistry of solid propellants. Proceeding of the fourth symposium of naval structural mechanics. - 1965. - P. 503-538.
2. Kiekbusch M. Elektrisches Messen von Volumenanderungen beim Triaxialversuch // Mitt. Bl. d. BAW. - 1977. - № 41. - P. 111-119.
3. Уржумцев Ю. С., П у т а н с А. В., К а л н р о з е З. В. Термоползучесть полиэтилена при циклических температурных воздействиях // Механика полимеров. - 1968. - № 3. - C. 421-427.
4. Matsushige K., Radcliffe S. V., B a e r E. The mechanical behavior of poly (methylmethacrylate) under Pressure // J. Polym. Sci.: Polymer Physics Edition.
- 1976. -V. 14. - P. 703-721.
5. Феодосьев В. И. Сопротивление материалов. - М.: Наука, 1972. - 544 с.
6. M u 11 i n s L. Softening of rubber by deformation // Rubber Chem. Technol. -1966. -41. - P. 339-362.
7. Onat E.,Prager W. The necking of tension specimen in plane plastic flow // J. Appl. Phys. - 1954. - V. 25. - № 4. - P. 491-493.
8. Хромов А. И. Разрушение жесткопластичных тел, константы разрушения // Изв. АН. МТТ. - 2005. - № 3. - C. 137-152.
9. Wiegand James H. The failure mechanism of solid propellant grains // Mechanics and chemistry of solid propellants. - 1965. - P. 539-574.
10. Комков К. Ф. Об использовании тензорно-нелинейных уравнений для анализа поведения пластических сред // Вестник МГТУ им. Н.Э. Баумана. Сер. "Машиностроение". - 2007. - № 1. - C. 46-56.
11. Лехницкий С. Г. Теория упругости анизотропного тела. - М.: Наука, 1977.
- 416 с.
12. Комков К. Ф.К определению параметров Лоде при обработке результатов испытаний // Изв. РАН. МТТ. - 2005. - № 2. - C. 126-135. 13. А р т е м о в А. Н., З и б р о в Л. А., М а т ч е н к о Н. М. Определяющие соотношения для нелинейных разно сопротивляющихся материалов // Проблемы прочности. - 1989. - № 4. - C. 39-42.
Статья поступила в редакцию 24.11.2006
Кирилл Федотович Комков родился в 1938г., окончил МВТУ им. Н.Э. Баумана в 1966 г. Канд. техн. наук, доцент военного технического университета. Автор более 30 научных работ в области механики материалов.
K.F. Komkov (b. 1938) graduated from the Bauman Moscow Higher Technical School in 1966. Ph. D. (Eng.), assoc. professor of the Military and Technical University. Author of more than 30 publications in the field of mechanics of materials.
