Кинетика зарождения локальных микродефектов при квазихрупком разрушении полимеров и композитов на их основе Текст научной статьи по специальности «Физика»

УДК 539.3

Кинетика зарождения локальных микродефектов при квазихрупком разрушении полимеров и композитов на их основе

© А. А. Валишин1,2, Т.С. Степанова1

1 Московский государственный университет тонких химических технологий им. М.В. Ломоносова, Москва, 119571, Россия 2 МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

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

Ключевые слова: полимеры, композиционные материалы на основе полимеров, трещина, зона вынужденной эластичности.

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

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

Локальные микроповреждения интенсивно образуются на медленной, флуктуационной, «спокойной» стадии развития трещины. Фрак-

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

Образование и накопление локальных микроповреждений происходит при всех температурах с разной интенсивностью, но при температурах выше температуры хрупкости перед фронтом трещины параллельно развиваются два процесса: накопление локальных повреждений и вынужденная эластическая ползучесть. Эти процессы взаимосвязаны, но до некоторой степени и независимы. Они протекают с разными скоростями, так как связаны с различными типами межчастичных взаимодействий и разными формами теплового движения. Характер теплового движения в полимерах подробно рассмотрен в работе [3]. Многочисленные экспериментальные исследования, подытоженные в [1], свидетельствуют о том, что элементарные акты разрушения в линейных полимерах связаны с термофлуктуационным разрывом химических связей, прежде всего, связей главной валентности полимерных макромолекул, не означающих, что межмолекулярное взаимодействие не вносит вклад в процесс разрушения. Оба типа взаимодействий формируют сложный потенциальный рельеф, на котором происходят элементарные акты разрушения. Кроме того, межмолекулярное взаимодействие играет определяющую роль в процессах вынужденного эластического деформирования.

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

Таким образом, выше температуры хрупкости в нагруженном полимерном образце перед фронтом трещины развиваются два процесса — накопление локальных повреждений и вынужденная эластическая ползучесть, приводящая к образованию перед трещиной разрушения зоны вынужденной эластичности. В работе [1] определены форма и размеры эластической зоны. Внешние напряжения (механического, термического или какого-либо другого происхождения) на микроскопическом уровне распределяются по химическим и межмолекулярным связям. За образование локальных повреждений ответственны химические связи, за вынужденную эластическую деформацию — межмолекулярные силы.

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

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

Наблюдаемое макроскопическое однородное, континуальное напряжение является результатом усреднения по временным промежуткам, много большим характерного времени теплового движения (например, периода колебаний около положения равновесия), и усреднения по пространственным объемам, много большим, чем характерный размер микронеоднородностей.

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

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

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

Не все элементы микроструктуры материала в одинаковой степени воспринимают внешнюю механическую нагрузку. В микроструктуре, как в строительной конструкции, есть некоторые «активные» элементы — «несущие балки, держащие напряжение», и есть «пассивные» элементы, которые непосредственно внешнюю нагрузку не воспринимают и играют роль «внутренних стен и перегородок». Так называемые несущие элементы микроструктуры — это, прежде всего, химические связи главной валентности полимерных макромолекул и группы таких связей. Разрушение таких элементов и составляет суть элементарных актов. Несущие элементы объединены в связный непрерывный «несущий каркас», как бы армирующий напряженный материал. Именно этот «несущий каркас» воспринимает и держит внешнее напряжение. Представления о каркасной связанности микроструктуры полимеров развивал Г. А. Патрикеев [4]. С математической точки зрения несущий каркас — это своеобразный случайный материальный граф, т. е. разветвленная сеть, охватывающая весь объем тела. Несущие элементы микроструктуры и несущий каркас в целом находятся в постоянном непрерывном тепловом движении. Микронеоднородность по различным физическим и механическим параметрам проявляется на всех уровнях, вплоть до химических связей. Экспериментально установлено [2], что химические связи в макромолекулах неодинаковы по энергии диссоциации. Так, в ПММА (полиметил-метакрилат) встречаются слабые химические связи С—С с энергией

кДж кДж ио = 134-и прочные — с энергией ио = 345- [2, 5]. Очевидно,

моль моль

что тепловыми флуктуациями в первую очередь рвутся слабые связи. Кроме того, и внешнее напряжение распределяется по химическим связям неравномерно. Методом ИК-спектроскопии установлено, что

вблизи трещины имеются перенапряженные химические связи, растягивающая нагрузка на которых в десятки раз превышает среднее макроскопическое напряжение [2].

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

В работе [1] показано, что в зоне вынужденной эластичности наблюдаемое макроскопическое напряжение выше, чем в окружающем объеме, однородно и равно пределу вынужденной эластичности оь. Однородность макроскопического напряженного состояния — это, как было сказано выше, результат сглаживающего пространственно-временного усреднения. В пространственных и временных микромасштабах, сравнимых с неоднородностями структуры, напряжение представляет собой сложное мозаичное случайное поле, но в среднем оно выше, чем в окружающем объеме. В разных местах микронеоднородной эластической зоны происходят элементарные термофлуктуационные акты разрыва и рекомбинации межчастичных связей. Каждый элементарный акт вызывает в малой окрестности дополнительное возмущение полей напряжений и деформаций. В результате, элементарные разрывы накапливаются в слабых местах объема эластической зоны, пока не образуется кластер, т. е. микрополость (дырка), которую уже можно зафиксировать экспериментально, например, методом малоугловой рентгеновской дифракции [2]. Вообще говоря, каждый отдельный акт макродеструкции предваряется множеством более мелких явлений, масштаб которых растет по мере развития процесса. Концентрация элементарных микроразрывов в отдельных местах эластической зоны и возникновение здесь микрополостей (дырок) дополнительно подтверждается экспериментальными данными по акустической эмиссии в полимерах под нагрузкой (и не только в полимерах) [6, 7].

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

на малую окрестность дырки и определяется ее размерами. Вдали от дырки поля напряжений и деформаций остаются невозмущенными. Поэтому, если вдали от первой дырки образовалась вторая, то эти два события независимы. Упругие поля этих дырок не перекрываются и не взаимодействуют. Однако, если две дырки возникли рядом, их упругие поля перекрываются на достаточно высоком уровне, то их совокупное упругое поле уже имеет больший радиус действия. Кроме того, в пространстве между дырками напряжения значительно выше, чем вдали от них, и «перемычка», разделяющая две близко расположенные дырки может «лопнуть» и тогда образуется большая дырка. Она создает еще большее возмущение с большим дальнодействием, и по ее периметру возможно возникновение новых дырок. Скопление дырок создает упругое поле с еще большим радиусом действия и т. д. В подтверждение сказанного уместно описать результаты, полученные в работе [8] на напряженных плоских оптически активных средах при различной концентрации дефектов. В круговой плоской области диаметром Ь создавались искусственно разрезы размером I « Ь, имитирующие трещины. Пока разрезов было мало, каждый из них создавал индивидуальное упругое поле с дальнодействием порядка I, характерное для изолированной трещины. При увеличении числа разрезов, сначала начинается взаимодействие полей соседних разрезов внутри области, а при концентрации разрезов, превышающих определенное значение, вся область создает внешнее поле, характерное для разреза размера Ь, т. е. скопление дефектов, каждое из которых имеет размер I « Ь, при достаточной их концентрации начинает действовать на расстоянии порядка размера всего скопления Ь.

Таким образом, до некоторого момента дырки возникают независимо друг от друга в разных местах эластической зоны. Когда дырок становится достаточно много, их упругие поля, перекрываясь, взаимно усиливают друг друга, увеличивается дальнодействие совокупного поля, что инициирует появление новых дырок, и вновь увеличивает радиус дальнодействия и т. д. Скопление (кластер) дырок стимулирует возникновение новых дырок, и процесс идет по нарастающей.

Началом быстрого распада (коллапса) эластической зоны является достижение некоторой критической концентрации дырок, при которой происходит так называемый перколяционный пробой, т. е. потеря устойчивости зоны, и трещина, которая «подпирает» клювом зону, быстро ее «проскакивает», впереди формируется новая эластическая зона, и все начинается сначала. Таким образом, рост трещины происходит микроскачками по мере достижения критической концентрации дырок в эластической зоне.

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

Механизм распада слабых узлов несущего каркаса. Все элементы микроструктуры материала находятся в состоянии постоянного теплового движения, интенсивность которого характеризуется температурой. Тепловые флуктуации — это спонтанные, случайные в пространстве и во времени выбросы энергии этого движения (подробно об этом в работе [3]). Результатом некоторых таких флуктуаций может быть разрыв межчастичных связей (а также их рекомбинация). В отсутствие внешнего поля между актами разрыва и рекомбинации устанавливается термодинамическое (т. е. статистическое) равновесие. Во внешнем поле напряжений (механического или термического происхождения), которое воспринимают элементы несущего каркаса, преимущество получают элементарные акты разрыва несущих элементов каркаса. Они происходят с разной интенсивностью в различных участках каркаса, но прежде всего в его слабых узлах. Элементарные акты разрыва несущих элементов в этих узлах постепенно накапливаются и, в конце концов, узел распадается, и на его месте образуется микрополость (дырка). Таким образом, возникновение дырок в эластической зоне связано с термофлуктуационным распадом узлов несущего каркаса. Дырка на месте такого узла — это конечный результат многих актов разрыва и рекомбинации несущих элементов этого узла.

При каждом элементарном акте разрыва или рекомбинации в напряженном материале изменяется деформация в месте происшествия и совершается силами упругих напряжений работа. Пусть в некоторой точке М произошел элементарный акт разрыва. В этом месте изменилось деформационное состояние (изменился тензор деформации), а именно возникла добавочная деформация типа растяжения. Эта деформация локализована в малой окрестности точки М. По мере удаления от этой точки возмущение поля деформаций затухает, так что существует радиус зоны возмущения и ее объем. Обозначим его как V 0 (М ). Именно в этом объеме и сосредоточена добавочная элементарная работа упругих сил, связанная с возмущением деформационного поля. Обозначим ее как

5А +(М)=-а1к (М)П +к(М) V +(М),

где о,к (М) — локальные напряжения в точке М, действующие на несущий элемент, участвующий в элементарном акте; п+к (М) — добавочная деформация в этой точке, обусловленная элементарным актом. Аналогично, если в точке М произошел элементарный акт рекомбинации, то здесь тоже возникнет возмущение деформационного поля — добавочная деформация п-к (М), но уже типа стягивания (сжатия). Она так же сосредоточена в некотором объеме V - (М). Силами упругих напряжений при этом совершается дополнительная работа

5А" (М)=-о 1к (М) п -к(М К (М).

Элементарные акты разрыва и рекомбинации несущих элементов происходят в узлах несущего каркаса. Дырка на месте узла возникает в результате множества элементарных актов, происходящих в узле. Узлы имеют различный объем, который определяется количеством несущих элементов в нем, как следствие, и дырки имеют различный объем и различное время зарождения. Обозначим объем узла в точке М, а значит, и объем возникающей на его месте дырки, как 5У(М). Объем узла (дырки) значительно больше объема элементарного возмущения деформационного поля, т. е. 5 У (М) » v± (М).

При каждом элементарном акте изменяется свободная энергия эластической зоны и даже всего образца. Но это изменение локализовано — оно сосредоточено в месте происшествия акта. Свободная энергия Г является полевой величиной с плотностью /(Т(М), ик(М)), в каждой точке М. Эта плотность зависит от локальной температуры и деформации в точке М. Свободная энергия какого-либо объема получается интегрированием ее плотности по этому объему. Плотность свободной энергии можно записать в виде

/(Т(М), ил(М)) = 8(М) - Т(М) Б(М),

где входящие сюда величины — локальные значения в точке М плотности внутренней энергии, температуры и плотности энтропии. Температура является интенсивным параметром, а все остальные величины — экстенсивные. Поэтому внутренняя энергия и энтропия какого-либо объема получаются интегрированием их плотностей по этому объему. Малое изменение плотности свободной энергии при изменении ее аргументов будет

А/(М)=-5(М)Д Т(М) + о±к(М, Т)Аи,к (М).

Изменение свободной энергии узла в точке М в результате одного элементарного акта сосредоточено в зоне возмущения деформационного поля объемом V0 и

Л^ (М) = -5 (М)ЛТ± (МК (М)-с% (М, /X (МК (М). (1)

Знак плюс в верхнем индексе относится к элементарному акту разрыва, а минус — к элементарному акту рекомбинации. Первое слагаемое определяет тепловой эффект элементарного акта, а второе—работу упругих сил при элементарном акте. В этих формулах (М, /) — это локальные микронапряжения в точке М на отдельных несущих элементах, участвующих в элементарном акте. Напряженное состояние в данной точке микронеоднородной среды, каковой является эластическая зона, можно представить состоящим из слагаемых:

а0к (М, /)=о'к (М) оа'к (М, /), (2)

где а0к — макроскопическое наблюдаемое напряжение, измеряемое в обычном макроскопическом эксперименте; ъ\к (М) — напряжение более низкого уровня — локальное напряжение на неоднородностях структуры, усредненное по времени. Оно образует случайное поле и описывается моментными функциями различного порядка (средними значениями, дисперсией, корреляционной функцией). Наконец, с]к (М, /) — напряжения наинизшего уровня, а именно, переменное, флуктуирующее во времени напряжение на отдельных несущих элементах структуры. Наблюдаемое макроскопически напряжение о0к рассчитывается обычными методами теории упругости. В работе [1] показано, что в эластической зоне впереди трещины нормального отрыва или поперечного сдвига устанавливается однородное напряжение, равное пределу вынужденной эластичности, т. е. а0к =ав. Элементарные акты разрыва и рекомбинации определяются последней составляющей тензора напряжений о]к (М, /). Поэтому, выполнив свертку по тензорным индексам в формуле (1) можем записать

ЛF ± (М ) = - 5 (М ) V 0 (М ) ЛТ ± (М ) - 2 к (М, / )овц± (М )v0 (М ), (3)

где безразмерный коэффициент к (М, /) характеризует распределение номинального напряжения ов по несущим элементам узла. Он зависит от точки М и от времени /, так как напряжение на несущих элементах узла меняется во времени по мере распада узла.

Структурные элементы несущего каркаса, воспринимающие внешнюю нагрузку, распределены по объему эластической зоны У0 с плот-

ностью р0(М). Это означает, что количество несущих элементов в узле объемом 5У(М) будет

5п0(М) = Ор0(М )5 У(М). (4)

Полное число слабых несущих элементов во всей эластической зоне обозначим ЛГ0(У0), оно зависит от объема эластической зоны и получается суммированием по всем узлам несущего каркаса, или интегрированием по объему эластической зоны:

*0(У0) = |5п0(М) = 0|р0(М)5У(М).

У0 У0

Отсюда нормировочный множитель О получим

О=ЗД)

|р0(М )5У (М)

Он имеет размерность обратного объема, т. е. [О] = да 3, а плотность распределения р0(М) — безразмерная функция. Плотность распределения несущих элементов в объеме эластической зоны представляет собой непрерывную функцию трехмерной точки М(х, у, г). Каждый узел несущего каркаса, на месте которого может возникнуть дырка, характеризуется двумя параметрами: мощностью р0(М) и объемом 5У(М), который тем больше, чем больше несущих элементов в узле. Случайный характер распределения слабых несущих элементов в пространстве обуславливает стохастический характер их разрушения и распада узлов несущего каркаса. Конкретизируем теперь коэффициент к(М, ¿) в формуле (3). Величина к(М, ¿)ов — это локальное напряжение, приходящееся на отдельный несущий элемент узла, участвующий в элементарном акте, т. е. то, что составляет последнее слагаемое в формуле (2). Интенсивность элементарных актов регулируется растягивающими составляющими напряжений на несущих элементах (например, силами, растягивающими химические связи и валентные углы). Поэтому напряжение к(М, ¿)ов — это всегда растягивающее напряжение. Макроскопическое напряженное состояние эластической зоны до возникновения дырок однородно, и напряжение всюду равно ов. В каждом узле идет процесс распада, число неповрежденных несущих элементов уменьшается со временем, а нагрузка на остающихся элементах возрастает. В какой-то момент количество неповрежденных элементов достигает минимального критического значения, а нагрузка на них достигает максимального критического уровня, тогда узел мгновенно

«лопается» и на его месте возникает дырка. И, как следствие, каждый узел через определенное время распадается, т. е. имеет некоторое «время жизни» /р(М), которое одновременно является временем зарождения дырки. Таким образом, узлы характеризуются тремя параметрами: мощностью р0(М), объемом ЪУ(М) и временем распада /р(М).

Начальное число несущих элементов в узле определяется формулой (4), поэтому начальное напряжение, приходящееся на один несущий элемент узла будет

Л(М)=^((М)ЪУ(М)). (5)

Обозначим количество неповрежденных несущих элементов узла М в произвольный момент / как

Ъп(М, /) = Ор(М, /)ЪУ(М) ^ Ъп0(М),

где р(М, /) — текущая плотность распределения неповрежденных несущих элементов в объеме эластической зоны. Текущая нагрузка на неповрежденных несущих элементах узла будет

/ (М, /)(р(М, /)ЪУ (М) )-1 > У0(М). (6)

Из этих формул вытекает соотношение между начальной р0(М) и текущей р(М, /) плотностями

р( М, /)=Р0(М), / = 0, р(М, /) <Р0(М), / > 0.

Из формулы (6) получим, что коэффициент к(М, /), характеризующий нагрузку на один неповрежденный несущий элемент узла, будет

|Р0(М )ЪУ (М)

к (М, /)=(Ор (М, /) ЪУ (М) )-1 = —У-,

^ ; #0(У0)р(М ,Т)ЪУ (М)

т. е. равен обратному количеству неповрежденных несущих элементов узла в момент /. Из формулы (6) следует, что текущее напряжение на один неповрежденный несущий элемент узла

/ (М, /)=к (М, / )ов = У0(М) .

р(М, /)

Вернувшись теперь к формуле (3) получим, что изменение свободной энергии при одном элементарном акте можно записать в виде

АF± (М) = - £(М)у± (М)АТ± (М) - 2к(М, Г)авл± (МК (М) = =-5 (М К (М )АТ( М)-2/(М, г )п± (М, г К (М)=

= -£(МК (М)АГ± (М) - 2/0(М)п± (МК (М). (7)

р(М, г)

Каждый элементарный акт разрыва или рекомбинации несущих элементов структуры (например, химических связей) сопровождается изменением свободной энергии образца на величину АF ± (М), локализованным в месте происшествия акта. Частоты элементарных актов можно записать в виде

ч Г АF1 (М)

ш±(М,г)=у0С (М)ехр - ; /

V кТ(М)

\

Здесь, как и раньше, знак плюс относится к элементарному акту разрыва, а минус — к элементарному акту рекомбинации. Подставив выражение (7) получим

ш± (М, г) =

^С±(М)ехр{-—^7 ((М)у±(М)АТ±(М)-2/(М,г)п±(М))).

V кТ (М) )

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

С ± (М)=ехр

{ и± ^

кТ (М)

или с поправкой на температурную зависимость энергии активации С + (М)=ехр V-кТМ) ((- Я+кТ (М))).

С учетом сказанного формулу (8) для частот элементарных актов представим в окончательном виде:

ш+(М, г)=у 0 ехр

(

1

ш (М, г)=у0ехр

кТ (М)

1

и+ - д+кТ (М) -- 5 (М X (М ) к А Т+ (М ) -

V

-2 /о(М )

Ро(М ) р(м , г)

п(М х (М)

кТ (М)

и0- - д - кТ (М)+

+ 2/о(М )

Ро(М ) р(М, г)

п(М X (М)

(9)

Кинетическое уравнение распада слабых узлов несущего каркаса. Получим кинетическое уравнение, описывающее процесс распада слабого узла. Обозначим количество несущих элементов узла в точке М в начальный момент г = 0 как

Ъп (М,0) = Ъп0(М) = Ор0(М) ЪУ(М).

Количество неповрежденных элементов этого же узла в произвольный момент г обозначим

Ъп(М, г) = Ор(М, г)ЪУ(М).

Изменение количества неповрежденных несущих элементов узла за малое время Дг происходит по двум причинам: во-первых, вследствие их разрыва за это время и, во-вторых, из-за рекомбинации поврежденных элементов за это же время. Составляя баланс поврежденных элементов за время Дг, разделив обе части на Дг и, перейдя к пределу при Дг ^ 0, получим уравнение

^(М,г)=-ш+(М,г)р(М,г)+ш-(М,г)[(М)-р(М,г)] ОЪУ(М). (10) д г

Это и есть кинетическое уравнение, описывающее процесс распада слабого узла каркаса в точке М. Частоты элементарных актов разрыва и рекомбинации несущих элементов определены формулами (9). Искомой функцией в уравнении (10) является р (М, г) — текущая плотность распределения несущих элементов в объеме эластической зоны. Начальное условие для этого уравнения будет

р (М,0) = р0( М ).

(11)

Исследование кинетического уравнения начнем с того, что выясним, существуют ли стабильные, нераспадающиеся узлы. Стабильность узла означала бы, что текущая плотность распределения несущих элементов р(М, г) в точке М не изменяется со временем, сохраняя

др

начальное значение р0(М). В этом случае производная — = 0 в точке М для любого г ^ 0. Из уравнения (10) тогда получим

ш+(М г)р0(М) = 0

и значит р0(М) = 0, т. е. нет распада там, где нечему распадаться. Отсюда следует, что стабильных узлов в поле механических напряжений не существует — все узлы несущего каркаса под действием напряжений и тепловых флуктуаций распадаются со временем, однако с различной скоростью.

Для дальнейшего исследования кинетического уравнения удобно ввести сокращенные обозначения в формулах (9):

а + (М)=и0+ - я+кТ(М) - £(МX (М)кАТ + (М), а - (М ) = и-- я - кТ (М ), Ь + (М ) = 2Л(М К (М ), Ь - (М ) = 2п(М )у0- (М ).

В этих обозначениях частоты элементарных актов записываются в следующем виде:

= V0 ехР

=У0ехр

ш+(М, г) ^ 0 ехр

ю+(М, г) кТ (М)

1

кТ (М) ш- (М, г) = v0 ехр

а+ (М ) - /0(М ) Ь+ (М )

ш (М, г) кТ(М)

1

кТ(М)

а- (М ) + /0(М ) Ь (М )

р0(( )

р(М, г)

рр(()

р(М, г).

(12)

С течением времени, по мере распада узла, частота элементарных актов разрыва ш+ увеличивается, а частота актов рекомбинации ш- уменьшается. Частота разрыва ш+ принимает наибольшее значение, когда показатель экспоненты в формуле (12) обращается в нуль. Это

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

^ММ^М® (13)

Р0(М ) а + (М )

Решив кинетическое уравнение (10) с начальным условием (11) найдем функцию р(М, г). Тогда из уравнения (13) можно найти время распада узла гр(М), т. е. время зарождения дырки в точке М. Таким образом, искать решение кинетического уравнения следует на отрезке г е [0, гр], а искомая функция р(М, г) должна быть монотонно убывающей от начального значения р0(М). Из формулы (13) следует, что время распада узла определяется его мощностью р0(М) — более мощные узлы «живут» дольше.

Рассмотрим теперь условие, при котором узел распадается сразу атермически, т. е. дырка зарождается практически мгновенно. Это произойдет в случае, когда уже в начальный момент частота актов разрыва ш + максимальна. Из (12) получим

г а+ (М) , ч

'к=ь+М). <14)

Формула (14) определяет критическое значение напряжения на несущих элементах узла в точке М. Если начальная нагрузка на несущих элементах /0(М) ^ /к, то дырка на месте этого узла образуется сразу же. Критическая нагрузка' не зависит от мощности узла, т. е. она одинакова для всех узлов. Кроме того, величины а+(М) и Ь+(М) в формулах (12), скорее всего, слабо зависят от точки М. Поэтому критическая нагрузка /к является константой материала, она определяется его микроструктурой и межчастичными силами взаимодействия.

Если начальная нагрузка на несущих элементах узла /0(М) < /к, то процесс распада узла и образование дырки растягивается во времени. Нагрузка на неповрежденных элементах узла возрастает пока не достигает критической /к, после чего оставшиеся несущие элементы практически мгновенно атермически рвутся, узел «лопается» и образуется дырка. Сопоставляя формулы (14) и (5), получим, что мгновенно «лопаются» узлы, мощность которых

Р0(М )ЪУ (М) —.

О а

Это наиболее слабые среди слабых узлов несущего каркаса. Решить кинетическое уравнение (10) можно, несколько его преобразовав, прежде всего, заменив искомую функцию р(М, г) на

1 (М, г) =

рр(М ) р(М, г).

В формулах (12) для частот элементарных актов введем дополнительные обозначения:

В ± (М ) ^0ехр Ь ± (М )

Г а ± (М ) ^

кТ (М)

а± (М)=

кТ (М)'

В новых обозначениях кинетическое уравнение и начальное условие к нему принимают вид

=В+1 ехр (((М ) а+1 )-В - (1 -1)2 5^(М )ехр (-/(М )а~г ),

1 (М ,0)=1. (15)

Искомая функция этого уравнения 1(М, г) монотонно возрастает на отрезке г е [0,гр]. Слагаемые в правой части уравнения (15) не равноценны, а именно второй член, описывающий рекомбинацию поврежденных несущих элементов, всегда меньше первого. Сопоставимость этих слагаемых зависит от значения /0(М) — начального напряжения на несущих элементах узла. Раньше было показано, что /0(М) определяется мощностью узла — чем более мощный узел, тем меньше /0(М). Отсюда следует, что рекомбинационный член в кинетическом уравнении вносит ощутимый вклад для высокомощных узлов в начале процесса распада, когда 1 близко к начальному значению 1 = 1. В дальнейшем с ростом 1 даже для высокомощных узлов вклад рекомбинации быстро убывает. Поэтому вполне оправданно искать решение кинетического уравнения пока без учета рекомбинационного члена, а затем определить, какую поправку вносит учет рекомбинации. Решая в этом приближении уравнение (15), получим

А+г = Е1 (-$1) - Е1 (-Р), где безразмерный параметр

Р(М) = а+/0(М).

Для дальнейшего упрощения используем асимптотику интегральной экспоненты. Тогда

е-р е-рг

Аи=е--—. (16)

в № У 7

Отсюда можно найти время распада узла t Получим

а+

1(М, tp) = /-1(М)—. (17)

Подставляя (17) в (16), после некоторых преобразований и упрощений, получим

кТ С а+ - /0(М )Ь+Л\ tp (М)=——— exp -; . (18)

/о(М

кТ

В (18) определяющую роль играет величина /0(М) — начальная нагрузка, приходящаяся на отдельные несущие элементы узла. Она находится с помощью формулы (5), согласно которой /0(М) обратно пропорциональна мощности узла р0(М). Это значит, что время распада узла (т. е. время зарождения дырки) полностью определяется мощностью узла — чем больше мощность, тем больше «время жизни» узла.

Параметр /0(М) входит в формулу (18) дважды: в предэкспонен-циальный множитель и в показатель экспоненты. Зависимость пред-экспоненты от/0(М) значительно более слабая, чем зависимость экспоненты. Поэтому обозначим

А = кТ

2 /о(м )л1>;

Здесь у0 — частота тепловых колебаний атомов, образующих несущие химические связи около положения равновесия, у0 ~ 1012-1013 с-1. Множители в этой формуле несопоставимы по значению: если у(-1 ~ 10-12-10-13 с, то первый множитель имеет порядок единицы. Время распада узла теперь окончательно можно записать так:

tp (Р0(М) ) = А ехр

С а+- /0 (р0( М ) )Ь + ^ кТ

(19)

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

ется в начале распада высокомощных узлов. Из формулы (19) видно, что с уменьшением /0(р0(М)), т. е. с увеличением мощности узла его «время жизни» возрастает. Экстраполяция в этой формуле на значение /0 = 0 дает конечное «время жизни» узла, в то время как оно, казалось бы, должно неограниченно возрастать. Но это не так. Значения мощности узлов р0 ограничены сверху, не может быть у узла сколь угодно высокой мощности. Значит, начальная нагрузка на несущие элементы узла /0 ограничена снизу, и «время жизни» любого узла конечно. Поэтому упомянутая экстраполяция не имеет смысла. Отступление от формулы (19) проявляется при малых близких к минимальному значениях /0 (высокомощные узлы) в сторону завышения «времени жизни» гр, но при любых обстоятельствах оно остается конечным. Численное решение кинетического уравнения (15) подтвердило эти выводы.

Кинетика накопления дырок в эластической зоне. Флуктуации теплового движения в поле механических напряжений приводят к постепенной деструкции несущего каркаса с образованием в его слабых узлах микрополостей (дырок). В ненагруженном теле флуктуаци-онно поврежденные несущие элементы каркаса (химические связи) с большой вероятностью тут же возвращаются в исходное состояние. Внешнее механическое напряжение вызывает, во-первых, снижение активационного барьера в направлении силы и, во-вторых, резкое уменьшение вероятности рекомбинации поврежденных элементов. Образно говоря, напряжение играет роль «вентиля», придающего хаотическому процессу повреждения и рекомбинации направленный характер.

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

Для того чтобы описать процесс накопления дырок, нужно ответить, по крайней мере, на три вопроса. Сколько имеется дырок в объеме эластической зоны в произвольный момент времени? Какая часть объема эластической зоны повреждена дырками в произвольный момент? Когда можно считать, что эластическая зона «переполнилась» дырками, т. е. каков критерий потери устойчивости зоны?

Каждый узел несущего каркаса имеет конечное время распада, после чего на его месте образуется дырка. Время зарождения дырки определяется формулой (19) и зависит, прежде всего, от мощности узла, т. е. от количества несущих элементов в нем. Несущие элементы в на-

чальный момент времени распределены в объеме эластической зоны с плотностью р0(М). Время зарождения дырки tр(р0(M)) является монотонно возрастающей функцией мощности узла. Плотность распределения несущих элементов р0(М), будучи непрерывной функцией трехмерной точки М в ограниченном объеме эластической зоны, достигает там наибольшего и наименьшего значения, которые определяют узлы с наибольшим и наименьшим временем распада. Отсюда следует, что множество (спектр) мощностей узлов дискретно, конечно и ограничено сверху и снизу, также и множество самих узлов дискретно и конечно. В произвольный момент времени t дырки будут в тех узлах, время распада которых не превышает t, т. е.

tр(Po(M)) < t или, имея в виду формулу (19),

Г а+- /о(Ро(М )) кТ

А ехр

< t.

У

Решение этого неравенства относительно плотности р0(М) дает

5«о(М) <ав-Ь—[=% ),

а+ - кТ 1п — А

где 8п0(М) = Ор0(М )8У(М) — количество несущих элементов в узле М. Это неравенство определяет те узлы, в которых к моменту времени t образовались дырки. Кроме того, оно определяет в объеме эластической зоны множество точек М, занятых дырками, т. е. поврежденную часть зоны. Поврежденный объем определяется как сумма объемов дырок, которую можно записать также в виде интеграла по поврежденному объему:

V ^) = | 8У (М).

8П0(М )< я)

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

N ^) = | ^Р0(М )8У (М).

8п0( М )< я)

Интегралы в этих формулах — это объемные, тройные интегралы.

Как было сказано, множество узлов несущего каркаса эластической зоны дискретно и конечно, хотя и велико. Пусть полное число узлов равно р(У0), оно зависит от объема эластической зоны, увеличиваясь с ростом последнего. Каждый узел характеризуется своей мощностью. Множество значений мощности — спектр мощности, тоже дискретен и конечен, поэтому его можно упорядочить по возрастанию мощности. Получим конечную числовую последовательность мощностей узлов: р1, р2, ..., р где р1 — наименьшее значение мощности, а рр — наибольшее. По признаку равномощности множество всех узлов разбивается на непересекающиеся классы. К каждому отдельному классу принадлежат узлы, имеющие одинаковую мощность; класс характеризуется кратностью ^ 1, равной числу узлов мощностью рк. Число таких классов £ не превосходит полного числа узлов, т. е. £ =5 р(У0). Таким образом, каждый узел помимо мощности имеет еще и свою кратность. Сумма всех кратностей равна количеству узлов:

£

=рю

к=1

Каждый узел несущего каркаса характеризуется, таким образом, четырьмя параметрами: мощностью рк, объемом 8Ук, временем распада грк и кратностью дк.

Построим функцию распределения узлов по их мощности, определив ее следующим образом:

Рк <р И

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

и гистограмму относительных кратностей —-, которая тоже описыва-

р

ет распределение узлов по их мощности.

Время распада узла гр является монотонно возрастающей функцией его мощности (см. формулу (19)), точнее функцией количества несущих элементов узла 5п0 = Ор 05У( р 0). Поэтому на основе возрастающей последовательности мощностей получим возрастающую последовательность времен распада гр1, гр2, ..., гр. Можно построить и функцию распределения узлов по временам их распада. Определим ее соотношением

фр(0 = I

'р к^

р ю'

Суммирование здесь идет по всем узлам, время распада которых меньше При каждом значении аргумента I функция Фр(0 определяет относительную долю узлов с временем распада, меньшим Количество дырок в произвольный момент времени I равно числу узлов, время распада которых меньше Следовательно, р Фр(^) в каждый момент равна числу дырок в этот момент.

Функции распределения Fр(р) = ^р(Ъп), где Ъп = ОрЪК(р), и Фр(0 связаны друг с другом, поскольку их аргументы связаны взаимно-обратными функциями:

С

7=g (Ъп)=А ехр

а+- V

Ъп

кТ

Ъп=g)=ав

а+ - кТ 1п

(20)

Эта связь функций распределения Fр(р) и Фр(0 выражается формулами

^р(р) = Фр( g(p)),

Фр(0 = ^(0).

Функции распределения Fр(р) = ^р(Ъп) и Фр(0, как видно из их определения, представляют собой кусочно-постоянные разрывные функции и являются полными аналогами эмпирических функций распределения математической статистики. С большой точностью их можно аппроксимировать непрерывными и даже гладкими функциями ^(р), Ф(0. Обе эти функции также назовем функциями распределения. Первая ^(р) описывает распределение слабых узлов по их мощности. Вторая Ф(^) описывает распределение дырок во времени и определяет относительное количество дырок в момент времени Связь между этими функциями распределения такая же, как и между их дискретными аналогами Fр(р) и Фр(0, т. е.

F(р) = Ф^(р)),

Ф(0 = ^40).

(21)

+

ь

I

Взаимно-обратные функции g( р ) и ^(г) обе монотонно возрастают, поэтому их производные положительны.

Итак, число дырок в произвольный момент г равно р(¥0)Ф(г), где р(У0) — число слабых узлов несущего каркаса эластической зоны, т. е. количество «вакансий» для дырок; оно растет с увеличением объема эластической зоны. Средняя концентрация дырок, т. е. число дырок в единице объема эластической зоны будет

р (V,) С (г)=Ф(г).

Она растет с течением времени. Когда концентрация дырок достигает некоторого критического значения Ск, эластическая зона уже неспособна «держать» напряжение ов и происходит ее коллапс, т. е. потеря устойчивости. Это наступает в момент гк, который определяется из условия

р(^)

V

Ф(гк)=Ск. (22)

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

Процесс накопления дырок в эластической зоне описывается функцией распределения времен распада слабых узлов Ф(г), которая, в свою очередь, определяется через функцию распределения узлов по их мощности ^(р) в соответствии с формулой (21). Функцию распределения ^(р) или ее плотность /(р) нужно выбрать из разумных физических предположений так, чтобы полученные дальнейшие результаты согласовывались с экспериментальными данными. Во-первых, плотность распределения /(р) отлична от нуля только на положительной полуоси р ^ 0; во-вторых, при р = 0 она должна обращаться в нуль, так как узлы нулевой мощности, т. е. несодержащие несущих элементов, отсутствуют; в-третьих, при больших значениях она должна достаточно быстро убывать, так как слишком высокомощных узлов в эластической зоне должно быть мало; в-четвертых, наибольшее количество узлов имеет некоторую «типичную» мощность, а остальные группируются около этого одного или нескольких типичных значений мощности. Соответственно кривая распределения /(р) является одно- или полимодальной. Распределение, удовлетворяющее всем этим условиям — это известное в теории надежности распределение Вейбулла и различные его частные случаи. Простейшим из них является распределение Рэлея, которое и выберем. Удобно будет плотность распределения Рэлея применительно к распределению узлов по мощности и записать в виде

I (Р) =

0, Р<0,

ехр

V "т У

, Р^0,

где р — мода распределения. Интегральная функция Р(р) будет

Р (Р) =

0, р<0,

(

1- ехр

\гт У

р^ 0.

Для определения текущего числа дырок, используя формулы (21) и (20), получим

Ф(Г) = 1- ехр

2 рт

а+ - кТ 1п —

V А)

Эта функция монотонно возрастает с течением времени от нуля до единицы и определяет относительное число дырок, т. е. долю дырок относительно общего количества слабых узлов. Для получения абсолютного числа дырок нужно умножить на число слабых узлов р(У0). Подставив функцию Ф(^) в уравнение (22) и решив его, получим «время жизни» зоны вынужденной эластичности:

(

гк = А ехр

1 о 1 +

а+ --^ Ъ+

кТ рт

V V

г

1п

1-

У

\

р(У0)

-а

2 ^ 2

У)

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

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

Получено кинетическое уравнение распада слабых узлов несущего молекулярного каркаса.

Исследована кинетика накопления дырок в зоне вынужденной эластичности перед трещиной разрушения. Найдено «время жизни» зоны вынужденной эластичности.

ЛИТЕРАТУРА

[1] Валишин А.А., Степанова Т.С. Особенности квазихрупкого разрушения полимеров и композитов на их основе. Инженерный журнал: наука и инновации, 2012, вып. 2. URL: http://engjournal.ru/catalog/mathmodel/material/52.html

[2] Регель В.Р., Слуцкер А.И., Томашевский Э.Е. Кинетическая природа прочности твердых тел. Москва, Наука, 1974, 560 с.

[3] Валишин А.А., Карташов Э.М., Тишаева С.Д. Статистические характеристики теплового движения в полимерах. Ученые записки МИТХТ, 2002, № 5, с. 11-18.

[4] Патрикеев Г.А. Механика полимеров, 1971, № 2, с. 221-231.

[5] Бартенев Г.М. Прочность и механизм разрушения полимеров. Москва, Химия, 1984, 280 с.

[6] Журков С.Н., Куксенко В.С., Петров В.А. О прогнозировании разрушения горных пород. Известия АН СССР. Физика Земли, 1977, № 6, с. 1-18.

[7] Половников Г.В., Трофимов В.В. Механика композиционных материалов, 1981, № 3, с. 542-546.

[8] Осокина Д.И., Мячкин В.И. и др. Физические процессы в очаге землетрясений. Москва, Наука, 1980, с. 68-78.

Статья поступила в редакцию 27.06.2013

Ссылку на эту статью просим оформлять следующим образом:

Валишин А. А., Степанова Т.С. Кинетика зарождения локальных микродефектов при квазихрупком разрушении полимеров и композитов на их основе. Инженерный журнал: наука и инновации, 2013, вып. 9. URL: http://engjournal. ru/catalog/mathmodel/material/1119.html

Валишин Анатолий Анатольевич — д-р физ.-мат. наук, профессор кафедры «Математическое моделирование» и «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана, профессор кафедры «Высшая и прикладная математика» Московского государственного университета тонких химических технологий им. М.В. Ломоносова. e-mail: enf@mail.ru

Степанова Татьяна Сергеевна — аспирант Московского государственного университета тонких химических технологий им. М.В. Ломоносова.