Использование принципов наименьшего действия при определении закономерностей разрушения (стр. 3 )

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

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

УДК 599.4.015

А., д-р техн. наук, профессор, чл.-корр. РААСН

Московский автомобильно-дорожный институт

, д-р техн. наук, профессор

Государственный технический университет

компьютерноЕ МОДЕЛИРОВАНИЕ И ИССЛЕДОВАНИЕ СВОЙСТВ

КОМПОЗИЦИОННЫХ МАТЕРИАЛОВ

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

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

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

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

Однако до последнего времени моделирование свойств КМ наталкивалось на трудности связанные с отсутствием математического аппарата. В большинстве случаев изменение свойств КМ матричного типа в зависимости от свойств и концентрации наполнителя пытаются описывать теорией «эффективной среды», суть которой заключается в том, что каждый отдельный элемент наполнителя КМ считается помещенным в однородную «эффективную среду» со свойствами, совпадающими с их истинными величинами для КМ в целом. В этом случае в результате расчетов получают «эффективные» значения различных характеристик КМ, исходя из свойств отдельных компонентов. Используя теорию «эффективной среды» многие исследователи обнаружили, что все наиболее известные методы и расчётные формулы данной теории (Оделевского, Лихтенекера, Бруггемана, Дульнева и т. п.) дают серьёзные отличия от экспериментальных данных начиная с некоторых значений объёмной концентрации наполнителя. На кривой изменения свойств КМ от объёмной концентрации наполнителя присутствует характерный излом при некотором её значении, зависящий от свойств и гранулометрического состава последнего (чем и объясняется отличия от эксперимента). Значение этой объёмной концентрации получило название «критической концентрации», а области на зависимости свойств от объёмной концентрации заполнителя получили название «докритических» и «послекритических». Теория «эффективной среды» не способна объяснить природу пороговых явлений.

Пороговые явления в изменении свойств КМ описываются теорией «просачивания» (percolation) [3], изучающей свойства связных компонент случайных графов. Эта теория, используя методику кластерных структур, описывает поведение неоднородных материалов в зависимости от концентрации компонентов в них. В силу гипотез универсальности и подобия характеристические параметры теории «просачивания» считаются независимыми от конкретной природы фазового перехода, т. е. различных материалов входящих в состав КМ.

Основным результатом теории «просачивания» для матричных материалов (состоящих из заполнителя и связующей матрицы) является достаточно простой степенной характер концентрационного поведения свойства композитного материала С.

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

Вероятно, единственным точным результатом в теории «просачивания» является результат, полученный для двухфазной системы Эфросом и Шкловским [4], которые аналитически показали, что в докритической и послекритической областях изменение свойств идет с одинаковым темпом.

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

Используя формулу Ландауэра-Бруггемана теории «эффективной среды» и накладывая на неё степенные характеристики теории «просачивания», были получены расчётные выражения для определения эффективных значений свойств КМ в докритической и послекритической областях объёмных концентраций наполнителя:

(1)

где Сэф – эффективное значение свойства КМ;

Сн – свойство наполнителя;

Vкон – текущая объемная концентрация наполнителя;

Vккон – критическая объемная концентрация наполнителя.

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

Используя решение выражения (1) относительно Vкон можно определить объёмную концентрацию наполнителя необходимую для получения КМ с заданными свойствами, т. е. осуществить расчёт его состава.

Единственной неизвестной величиной в (1) является значение Vккон. Критическая концентрация является наиболее важным и трудноопределимым критическим индексом теории «просачивания». Для ее определения обычно используют физические модели на основе решеток того или иного типа, но значения Vккон, полученные таким образом, справедливы только для решеточных задач и не справедливы для континуальных, для которых данные значения практически не получали.