авторефераты диссертаций БЕСПЛАТНАЯ  БИБЛИОТЕКА

АВТОРЕФЕРАТЫ КАНДИДАТСКИХ, ДОКТОРСКИХ ДИССЕРТАЦИЙ

<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ

Моделирование свойств композиционных материалов, дисперсно армированных жесткими короткими волокнами

УЧРЕЖДЕНИЕ РОССИЙСКОЙ АКАДЕМИИ НАУК

ИНСТИТУТ ПРИКЛАДНОЙ МЕХАНИКИ РАН

На правах рукописи

Гордеев Андрей Владимирович

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

МАТЕРИАЛОВ, ДИСПЕРСНО АРМИРОВАННЫХ

ЖЕСТКИМИ КОРОТКИМИ ВОЛОКНАМИ

01.02.04 – Механика деформируемого твердого тела

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата физико-математических наук

Москва 2010 2

Работа выполнена в Учреждении Российской Академии Наук Институте Прикладной Механики РАН

Научный руководитель: доктор технических наук, профессор Лурье Сергей Альбертович

Официальные оппоненты: доктор физико-математических наук, Березин Александр Васильевич, доктор технических наук, профессор Дудченко Александр Александрович

Ведущая организация: «МАТИ» – Российский государственный технологический университет им.

К.Э.Циолковского

Защита состоится “”2010 г. в XX00 часов на заседании диссертационного совета Д212.125.05 в Московском авиационном институте (государственный технический университет), кафедра 902 по адресу: 125871, г.Москва, ГСП, Волоколамское шоссе, дом 4.

С диссертацией можно ознакомиться в научной библиотеке Московского авиационного института.

Автореферат разослан _ _2010 г.

Ученый секретарь Г.В. Федотенков диссертационного совета, кандидат физико-математических наук

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Принципиальные результаты при разработке прикладных моделей новых композиционных конструкционных материалов, методов оценки эффективных свойств, методов оценки разрушения, накопления повреждений композитов, зависимости процессов поврежденности от свойств структуры, а также асимптотических итерационных и иных приближенных методов исследования деформаций элементов конструкций, несущей способности композитных конструкций связаны в первую очередь с именами таких ученых как: Н.С. Азиков, А.В.Бабешко, А.В.Березин, В.А.Бунаков, Г.А.Ванин, В.В.Васильев, И.И.Ворович, Р.В.Гольдштейн, Э.И.Григолюк, А.А. Дудченко, А.Н. Елпатьевский, С.А. Лурье, Ю.М.Новичков, И.Ф.Образцов, А.Н.Полилов, Б.Е.Победря, Н.Н.Рогачева, Салганик Р.Л, Ю.Н.Тарнопольский, K.Б. Устинов, а также B.Budiansky, R. M. Christensen, A.V. Dyskin, L.N.Germanovich, Z.Hashin, M.Kachanov, S. Kanaun, T.Mori, K.

Tanaka, T.Mura, G.M. Odegard и др.

Открытие в 1991 году длинных, цилиндрических углеродных многослойных и однослойных структур получивших название нанотрубок (УНТ) стало отправной точкой в создании нового вида композиционных материалов – нанокомпозитов. Как известно, углеродные нанотрубки обладают механическими характеристиками, значительно более высокими, чем у стали. Одна из наиболее очевидных возможностей использования этих уникальных объектов связана с созданием нанокомпозитов, т. е. полимерных материалов (матриц), содержащих некоторое, весьма небольшое количество УНТ. В этом случае говорят о модификации полимерных матриц. Главная цель такой модификации состоит в целенаправленном изменении свойств матрицы. Ожидается, что малое объемное содержание нановключений изменяя некоторые базовые механические свойства матриц, позволит оставить без изменения их технологические качества, являющиеся существенными при изготовлении неоднородных армированных композитов в целом, что позволит получить улучшенные конструкционные композитные материалы без существенного изменения технологии изготовления. При этом главная трудность состоит в обеспечении хорошей адгезии между поверхностью нанотрубки и молекулами полимерной матрицы. При плохой адгезии нанотрубки внутри матрицы при малых объемных долях могут фактически никак не влиять на свойства матриц, или даже приводить к ухудшению свойств.

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

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

Разработке модели, позволяющей предсказать конечные свойства нанокомпозитов, и посвящена данная работа.

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

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



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

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

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

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

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

- Предложена процедура оценки свойств «изотропного»

нанокомпозита, армированного нановолокнами с изотропным распределением по углам с использованием градиентной модели.

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

Разработанный в диссертации метод и алгоритмы могут быть рекомендованы для проектных и научно-исследовательских организаций.

Реализация результатов работы. Результаты, полученные в диссертации, используются в Учреждении Российской Академии Наук Институте Прикладной механики РАН, ОАО НИАТ, ВИАМ.

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

Апробация работы. Результаты диссертации были представлены на всероссийской конференции, приуроченной 20-летию ИПРИМ РАН:

«Механика и наномеханика структурно-сложных и гетерогенных сред.

Успехи, проблемы, перспективы», 2010г, Москва;

XVIII Всероссийской школы-конференции молодых ученых и студентов «Математическое моделирование в естественных науках», 2009, Пермь.

Основные результаты диссертации были доложены на семинарах:

- Семинар по механике гетерогенных сплошных сред ИПРИМ РАН;

- Научный семинар им. А. Г. Горшкова «Проблемы механики деформируемого твердого тела и динамики машин» в МАИ (под руководством д.ф.-м.н., проф. Д. В. Тарлаковского, д.т.н., проф.

Ф. Н. Шклярчука, д.т.н., проф. В. В. Фирсанова).

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

На защиту выносятся:

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

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

- Алгоритм расчёта эффективного модуля Юнга, в рамках модели, учитывающей сдвиговые адгезионные взаимодействия.

- Алгоритм расчёта эффективного модуля Юнга, в рамках модели, учитывающей адгезию сдвига и поверхностного натяжения на границе контакта матрица/включение.

- Алгоритм расчёта эффективного модуля Юнга, в рамках градиентной модели композита, учитывающей когезионные взаимодействия.

- Алгоритм расчёта эффективного модуля Юнга для изотропного нанокомпозита с разноориентированными волокнами, построенный в рамках градиентной модели.

- Алгоритм расчёта эффективного модуля Юнга, в рамках градиентной модели, учитывающей когезионные взаимодействия, адгезию сдвига и поверхностного натяжения.

Объём и структура работы. Диссертационная работа состоит из введения, шести глав, заключения и списка литературы. Она содержит страницы, из них 10 занимает список использованных источников. Список используемой литературы включает 81 наименование (из них 46 на иностранном языке).

СОДЕРЖАНИЕ РАБОТЫ

Во введении обосновывается актуальность научных исследований, изложенных в диссертации, а также сформулированы: цель исследования, научная новизна, практическая и теоретическая ценность работы.

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

Рассматриваются и анализируются методы моделирования нанообъектов и нанокомпозитов. Приведена краткая характеристика метода молекулярного моделирования, область применения и возможности модели.

Рассматриваются и анализируются методы, основанные на механике сплошной среды.

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

Во второй главе рассмотрена задача по расчету модуля Юнга композиционного материала, армированного короткими продольно ориентированными волокнами, в рамках классической теории упругости.

В главе строятся, а затем сравниваются общая постановка задачи в осесимметричной (рис.1) и плоской постановке (рис.2), где в соответствии с методом В.З.Власова, распределение перемещений вдоль оси волокна задается в виде линейной функции x : R( x, y ) xR1 ( y ).





Y y x RM rd Z X l l l 2 Рис. 1. Цилиндрический Рис. 2. Фрагмент композита, с фрагмент композита, с длиной длиной волокна l. Плоская волокна l. Осесимметричная постановка задачи.

постановка задачи.

Записывается вариационное уравнение, определяющее расчетную модель для контактной задачи волокно-матрица в осесимметричной постановке:

rD GD l 2 G l3 rD (rR1 ) E D rR1 rP}R1dr 2 D rR1R L 2l { 12 (1) RМ GM l 2 G l3 RМ ) E M rR1 rP}R1dr 2 M rR1R 2l (rR { rD 12 rD Индексом M обозначены параметры, относящиеся к матрице, а индексом D - параметры, относящиеся к волокну (дефекту). В вариационном уравнении l - длина включения, EM, D, GM, D, - модуль Юнга и модуль сдвига, P - кусочно постоянная нагрузка, rD - радиус волокна, RM – внешний радиус слоя матрицы (рис.1).

В качестве упрощенной постановки приводится плоская постановка рассматриваемой проблемы. В соответствии с методом В.З.Власова вариационное уравнение приводится к виду:

hD GD l 2 G l3 hD R1 E D R1 P}R1dy D R1R L l { 12 (2) hD hM G l2 G l3 hD hM { M R1 E M R1 P}R1dy M R1R l hD 12 hD Здесь hM,D - толщина матрицы и толщина включения (рис.2).

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

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

В результате, для классической постановки получено следующее соотношение для расчёта эффективного модуля эквивалентного континуума:

1 1 f * f * (3) E EM ED * Здесь f - относительная объемная доля для эквивалентного континуума. Для осесимметричной постановки f * будет:

rd С J 12ED ir C Y 12ED ir rdr f f 1D 0 D G D l 2 0D 2D GD l 21 1 RМ P E E D M (4) RМ С J 12EM ir C Y 12EM ir rdr 0M 0M 1M 2M 2 GM l GM l 2 1 1 rd RМ P E E D M r f d Здесь - относительная объемная доля включений, RM 12 E D 12 E D 12 E M 12 E M J 0D ir, Y0 D ir, J 0 M ir, Y0 M - функции Бесселя G l 2 ir Gl Gl Gl 2 2 D D M M первого и второго рода мнимого аргумента, С1D, C 2 D, С1M, C 2 M - постоянные интегрирования в решениях для волокна и матрицы соответственно.

Для плоской постановки получаем:

hf f* f (5) (hM hD ) hD f (hD hM ) Здесь: - относительная объемная доля включений в плоской постановке;

величина h f, в дальнейшем называемая как «обменная»

объемная доля, вычисляется по формуле:

( E D E M ) D h D M hM hf (6) [ E M M hM E D D hD ] 1 12 E M, D - показатели затухания классического краевого эффекта;

bM, D l GM,D th(bM, D hM, D ) - относительная объемная доля однородного M,D (b M, D h M, D ) состояния.

Чтобы пояснить физический смысл величины h f выражение для эффективного модуля упругости переписывается в виде:

hM hD E hM h f hD h f (7) EM ED Физический смысл h f состоит в том, что объемная доля слабой фазы уменьшается - матрица «отдает» включению часть своей объемной доли.

Далее приводятся трактовки различных форм полученного решения.

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

1 1 f f * (8) E EM ED где:

* ED ED ( E EM ) h f (9) 1 D EM hD Ниже соотношение для эффективного модуля представлено в форме, соответствующей гипотезе эквивалентной матрицы. Здесь также как и ранее формула расчёта эффективного модуля также принимает вид осреднения по Рейсу:

1 (1 f ) f (10) * E EM ED где * E M EM ( E EM ) h f (11) 1 D ED hM Аналогично предыдущим случаям получается представление для эффективного модуля в виде осреднения по Рейсу с точки зрения гипотезы трех фаз:

hM hD hM hMf hD hDf (hDf hMf ), (12) E EM ED Ef где:

1 hDf E 1 D[ ] E D D hD EM M hM (13) hMf 1 1 EM [ ] E D D hD EM M hM (E D EM ) Ef (14) Здесь величина связана с суммарной толщиной межфазного слоя hf (E D E M ) следующим соотношением. Полученный (hDf hMf ) h f ( h Df h Mf ) (ED EM ) результат соответствует классическому осреднению Рейсса трехфазного композита с модулем третьей фазы E f.

Далее приводится доказательство того, что в предельных случаях при больших длинах волокон формула осреднения принимает вид классического осреднения по Фойхту, а при малых длинах классического осреднения по Рейсу.

Результат расчета величины эффективного модуля нанокомпозита в осесимметричной постановке представлен в виде графиков на рис. 3. На рисунке показаны зависимости значений модуля Юнга от длины волокна, полученные по классическим формулам Фойхта и Рейсса, показано решение, полученное на модели, рассмотренной во второй главе, а также приведены данные из работы Г.М.Одегарда. Исходя из решений, показанных на рис.3, делается вывод о том, что результат, полученный на предложенной модели, не объясняет эффект усиления, описанного Одегардом и др.

Е - Эффективный модуль Юнга композита, ГПа Данные Одегард 8 и др.

Осреднеие по Фойгту Решение в Ev=5.2507ГПа осесимметричной постановке Осреднеие по Рейссу Er=0.9091ГПа 0 50 100 150 200 250 300 350 400 l - длина волокна, нм Рис. 3. Зависимость эффективного модуля Юнга композиционного материала от длины волокна. Осесимметричная постановка, а также данные численного эксперимента.

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

Осреднение по Фойгту Данные Ev=11.015ГПа Одегард и др.

Е - Эффективный модуль Юнга композита, ГПа Решение в осесимметричной постановке Осреднение по Рейссу Er=0.9091ГПа 0 50 100 150 200 250 300 350 400 l - длина волокна, нм Рис. 4. Зависимость эффективного модуля Юнга композиционного материала от длины волокна с примененимем эквивалентного модуля включения Ed*=1012,4ГПа Фактически тем самым предложена процедура приведения полученных теоретических результатов к результатам численного эксперимента, найденных группой Одегарда.

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

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

На рис. 5 показаны решения в осесимметричной постановке при фиксированной объемной доле включений и разных диаметрах включений.

Осреднение по Фойгту Данные Ev=11.015 ГПа Одегард и др.

Е - Эффективный модуль Юнга композита, ГПа Решение в осесимметричной постановке при фиксированной объемной доле и различных радиусах включений rd от 0,3 до 3 нм Осреднение по Рейссу Er=0.9091 ГПа 0 50 100 150 200 250 300 350 400 l - длина волокна, нм Рис. 5. Зависимость эффективного модуля Юнга композиционного материала от длины волокна при фиксированной объемной доле включений 1%, эквивалентном модуле включения Ed*=1012.4 ГПа и различных радиусах включений. Осесимметричная постановка.

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

Аналогичный эффект проявляется и для модели в плоской постановке, что позволяет подобрать толщину hD таким образом, чтобы решение совпадало с полученным для этой модели в осесимметричной постановке (рис.6).

Осреднение по Фойгту Ev=11.015 ГПа Е - Эффективный модуль Юнга композита, ГПа Решение в осесимметричной постановке Данные Одегард Решение в плоской и др.

постановке Осреднение по Рейссу Er=0.9091 ГПа 0 50 100 150 200 250 300 350 400 l - длина волокна, нм Рис. 6. Зависимость эффективного модуля Юнга композиционного материала от длины волокна при фиксированной объемной доле и эквивалентном модуле Ed*=1012.4ГПа. Плоская постановка совпадает с осесимметричной при толщине hD 0.133 нм В главе были получены следующие результаты:

1. При анализе результатов обнаружен эффект изменения эффективного модуля композита в рамках вилки Фойгта-Рейсса с изменением диаметра включения при фиксированной объемной доле и длине волокон. Точнее: эффективный модуль является монотонно убывающей функцией от hD / l. Появление этого эффекта объясняется тем, что предложенная в главе модель, учитывает краевые эффекты, возникающие поперек оси волокна и играющие роль «межфазного слоя». При уменьшении диаметра включений растет относительная длина краевого эффекта, а значит и суммарная толщина «межфазного слоя».

2. На основании этого эффекта предложены технологические рекомендации: повышение/понижение эффективного модуля композита путем использования волокон меньшего/большего диаметра при фиксированной длине и объемной доле.

3. В качестве параметра «межфазного погранслоя» введена «обменная доля» композита h f h f ( ED, EM, GD, GM, hD, hM, l ), которая определяет его толщину. Необходимо отметить, что в предложенной модели межфазного погранслоя отсутствуют параметры, связанные с изменением морфологии материалов. Однако, очевидно, что изменение морфологии материала в областях, прилегающих к зоне контакта матрица/волокно, является существенным фактором, влияющим на повышение эффективного модуля. Следовательно, дальнейшие разработки моделей должны проводиться в этом направлении.

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

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

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

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

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

Модель композита с адгезионным взаимодействием волокна и матрицы рассматривается в плоской постановке. Решение также строится методом В.З.Власова, вариационное уравнение при этом имеет вид:

hD hD hM GDl 3 GM l rD ]rD dy L rM ]rM dy [ PD l E D lrD [ PM l E M lrM 12 0 hD l3 l3 (15) G D rD (0)rD (0) [G D rD (hD ) g g rD (hD ) g g rM (hD )]rD (hD ) 12 l3 l [G M rM (hD ) g g rM (hD ) g g rD (hD )]rM (hD ) G M rM (hD hM )rM (hD hM ) 12 Отличие этого вариационного уравнения от рассмотренных ранее состоит в том, что в уравнение добавляются новые члены с дополнительными параметрами, соответствующие модели типа винклера на сдвиг для адгезионного слоя. Величина G g определяет модуль сдвига фиктивного слоя, величина g g G g / h g определяет удельный адгезионный сдвиговой модуль адгезионного слоя, hg - толщина адгезионного слоя.

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

hM hD E hM h f hD h f (16) EM ED где ( E D E M ) D h D M hM hf 12 E M M hM E D D hD (17) [ E D D h D E M M hM ] ggl Из полученных уравнений можно сделать следующие выводы:

1. При g g величина h f стремится снизу к соответствующей величине h f, полученной при идеальном контакте (см. (6)), когда волокно не проскальзывает относительно матрицы.

2. При g g 0 величина h f стремится к нулю и, соответственно, эффективный модуль стремится к величине, определяемой осреднением по Рейссу.

3. При эффективный модуль будет удовлетворять 0 gg неравенствам: для любой длины волокна.

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

Кривые на рис.8 демонстрируют зависимости осредненного модуля композита от длины включения при идеальном контакте, при контакте с учётом адгезии и экспериментальные данные группы Одегарда.

Осреднение Данные по Фойгту Одегард Ev=11.015 ГПа и др.

g=0. Решение ГПа/нм без учета адгезии g=0. ГПа/нм Е - Эффективный модуль Юнга композита, ГПа g=0. ГПа/нм g=0. ГПа/нм Решение с учетом адгезии g=0. Осреднение ГПа/нм по Рейссу Er=0.9091 ГПа 0 50 100 150 200 250 300 350 400 l - длина волокна, нм Рис. 8. Зависимость эффективного модуля Юнга композиционного материала от длины волокна при фиксированной объемной доле 1%, эквивалентном модуле включения Ed*=1012.4ГПа. Адгезионный параметр g изменяется от 0,0005 ГПа/нм до 0,3125 ГПа/нм Как видно из приведенных графиков, перебор параметра g g позволяет описать различные условия взаимодействия матрицы и волокна. На рисунке показано, что полученные кривые при различных значениях сдвигового адгезионного параметра g g лежат между кривой, определенной полным геометрическим контактом, рассмотренной в главе 2 и осреднением по Рейссу.

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

1. В качестве параметра межфазного слоя предложена «обменная»

объемная доля композита в виде h f h f ( E D, E M, G D, G M, hD, hM, l, g g ), которая определяет его толщину. Необходимо отметить, что в предложенной модели межфазного слоя присутствует параметр g g, связанный с взаимодействием пары матрица/волокно на границе контакта. Таким образом, выражение (16) для вычисления эффективного модуля совпадает с выражением (7), полученным во второй главе. Выражения (6) и (17) для вычисления величины h f в этих моделях отличаются.

2. Адгезионный параметр g g может играть роль параметра качества технологии при изготовлении композитов. Нужно отметить, что существуют технологические приемы (например, предварительная активация поверхности волокон, функционализация), которые привели бы к более высоким значениям эффективного модуля при той же объемной доле волокон. И наоборот, «загрязние» поверхности волокон, может свести эффект усиления к нулю, какая бы большая объемная доля волокон не использовалась. Здесь учет адгезии влияет на эффективные свойства, как поврежденность и снижает эффективный модуль. Следовательно, при исследовании перспективных пар матрица-волокно следует иметь в виду как принципиальную возможность дополнительного увеличения эффективного модуля композита за счет целенаправленного увеличения адгезии волокна к матрице, так и возможность его снижения.

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

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

Решение строится методом В.З. Власова. Вариационное уравнение получено в соответствии с принципом Лагранжа из требования стационарности лагранжиана:

hD hD hM GDl 3 GM l L rD ]rD dy rM ]rM dy [PD l ED lrD [PM l EM lrM 12 0 hD l GD rD (0)rD (0) l {GM rM (hD ) / 2 GD rD (hD ) / 2 12eg / l 2 [rM (hD ) rD (hD )]} [rM (hD ) rD (hD )] (18) l {GM rM (hD ) / 2 GD rD (hD ) / 2 g g [rM (hD ) rD (hD )]} [rM (hD ) rD (hD )] l GM rM (hD hM )rM (hD hM ) Здесь g g Gg / hg - удельный адгезионный модуль сдвига, eg E g hg / 4 адгезионнная жесткость адгезионной пленки на растяжение.

Эффективный модуль композита определяется из условия равенства потенциальных энергий для композита и эквивалентной однородной среды.

Формула для расчета имеет вид:

hM hD E hM h f hD h f (19) EM ED Здесь величина h f («обменная» объемная доля) имеет вид:

h f (C Db D hD C Mb M hM ) [ l2 E E E E l ( M M hM e g )( D D hD g g ) ( M M hM g g )( D D hD e g ) 2 2 12 2 12 P P ( M D) l2 l2 E EM ED E D D hD g g ) g g ( D D hD e g )} D hD {e g ( P P 2 12 12 ( M D) (20) EM E D l2 E E E E l ( M M hM e g )( D D hD g g ) ( M M hM g g )( D D hD e g ) 2 2 12 2 12 P P ( M D) l lE EM ED EM { g g ( M M hM e g ) e g M hM g g )} M hM ] ( P P 12 2 2 ( M D) EM E D На рис. 9 приведены данные, полученные Одегардом и др. для двух разных случаев адгезионных взаимодействий (функционализированное и нефункционализированное волокно), а также решения модели, сформулированной в данной главе. При значении параметра g g и изменении e g, все кривые эффективных модулей лежат между кривой, не учитывающей адгезию, кривой определяемой параметрами и gg e g e g _ max. В модели, показанной на рис. 9, параметр e g изменялся в диапазоне от 2 до 24 ГПа нм. На графике видно, что силы поверхностного натяжения приводят к усилению композита в диапазоне длинных волокон (рис.9), если пролагать возможность увеличения приведенной жесткости адгезионной пленки на сдвиг и на растяжение за счет соответсвующих технологических приемов (функционализации). При этом, модель не может объяснить эффект аномального усиления, возникающий на коротких волокнах.

gg=inf.

eg= Данные Одегард Е - Эффективный модуль и др.

Юнга композита, ГПа адгезия Данные Одегард и др.

адгезия gg=inf.

eg= 6 gg=inf.

eg= Осреднение по Фойгту gg=inf.

Ev=5.2507 ГПа eg= Решение с учетом адгезии gg=inf.

eg= Осреднение по Рейссу Решение Er=0.9091 ГПа без учета 0 50 100 150 200 250 300 350 400 адгезии l - длина волокна, нм Рис. 9. Зависимость эффективного модуля Юнга от длины волокна с учетом адгезии. Учитываются только поверхностное натяжение g g =inf, e g произвольно Полученная модель и ее решение позволяют сделать следующие выводы:

1. В качестве параметра межфазного слоя предложена величина h f h f ( E D, E M, G D, G M, hD, hM, l, g g, e g ), названная «обменной долей»

композита, которая определяет его приведенную толщину. В предложенной модели межфазного слоя присутствуют параметры g g, e g, связанные с взаимодействием пары матрица/волокно на границе контакта. Таким образом, выражения для вычисления эффективного модуля совпадают с полученными во второй (7) и четвертой главе (16), отличаются лишь выражения для вычисления величины h f.

2. Адгезионные параметры g g и e g могут играть роль параметров качества технологии изготовления композитов.

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

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

В пятой главе также описана достаточно простая детерминированная модель композита с изотропным распределением волокон но углам ориентации.

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

Для градиентной модели вариационное уравнение выглядит следующим образом:

hD hM l2 G2 l r r }rdy L l { p Er G C (21) l2 G2 l lG r )Gr (G r r )r} l{( y hD 12 C C Здесь r ( y) - искомая неизвестная функция. В вариационном уравнении имеется дополнительное слагаемое с дополнительным градиентным модулем С.

Градиентный (когезионный) параметр С, имеет смысл градиентного модуля упругости, отличающегося на квадрат длины от модуля упругости.

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

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

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

Таким образом, формула для вычисления эффективного модуля по-прежнему имеет вид:

hM hD E hM h f hD h f (22) EM ED Где h f :

hD C ch D y cos D y C D 2 sh D y sin D y dy hf D 1 P( ) EM E D hD hM C ch M ( y hM hD ) cos M ( y hM hD )dy M 1 1 (23) P( ) hD E M ED hD hM C M 2 sh M ( y hM hD ) sin M ( y hM hD ) dy 1 P( ) hD E M ED Здесь C D1, C D 2, C M 1, C M 2 - постоянные интегрирования.

На рис.10 показаны решения при изменении когезионного модуля матрицы.

Осреднение по Фойгту Ev=11.015 ГПа Данные Одегард и др.

Решение с учетом когезии Cm от 0.00009 ГПа/нм до 0.06561 ГПа/нм Решение Е - Эффективный модуль без учета когезии Юнга композита, ГПа Осреднение по Рейссу Er=0.9091 ГПа 0 50 100 150 200 250 300 350 400 l - длина волокна, нм Рис. 10. Зависимость эффективного модуля Юнга композиционного материала от длины волокна при фиксированной объемной доле 1%, эквивалентном модуле Ed*=1012.4ГПа. Когезионный модуль С M менялся в пределах от 0,00009 ГПа/нм2 до 0,06561ГПа/нм Решение задачи моделирования композиционного материала в рамках градиентной модели принципиально не отличается от задачи с растяжением вдоль оси волокна. Значение осредненного значения модуля Юнга показывает хорошее согласование с данными, приведенными Одегардом и др., позволяя моделировать эффект усиления для композита, армированного короткими тонкими волокнами (микро- и нано-волокнами).

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

Основные выводы по главе:

1. В качестве расчетного параметра межфазного слоя используется величина h f h f ( E D, E M, G D, G M, hD, hM, l, С D, С M ), которая определяет его приведенную толщину. В предложенной модели межфазного слоя присутствуют параметры СD,СM, определяющие когезионные взаимодействия как в волокне, так и в матрице. Выражение для вычисления эффективного модуля совпадает с полученными в предыдущих моделях (7), (16), (22), отличается лишь выражение для вычисления величины h f.

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

Наибольшее влияние на эффективный модуль оказывает изменение морфологии матрицы в области, близкой к зоне контакта с волокном.

Материал матрицы с измененной морфологией образует дополнительный слой, характерная ширина которого определяется по формуле l M G M / C M, и для рассматриваемого композита превышает радиус включения в ~6 раз. Рассмотренный эффект не связан с геометрией структуры композита и носит характер масштабного эффекта.

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

Решение строится методом В.З. Власова. Распределение перемещений выбирается в виде R( x, y ) xr ( y ). Вариационное уравнение в данной постановке имеет вид:

hD GD l 2 G l 2 IV rD D L l [PD ED rD rD ]rD dy C D hD hM G l2 G l 2 IV [PM EM rM M rM M rM ]rM dy l CM hD l 2 GD l2 G l2 y hD rD )GD rD (GD rD D rD )rD } l{( (24) y 12 CD C D l 2 GM l2 G l2 y hD hM )GM rM (GM rM M rM )rM } l{( rM y hD 12 CM CM ggl 3 ggl [rM (hD ) rD (hD )]rD (hD ) [rM (hD ) rD (hD )]rM (hD ) 12 e g l[(rM (hD ) rD (hD ))]rD (hD ) e g l[(rM (hD ) rD (hD ))]rM (hD ) Эффективный модуль композита определяется из условия равенства потенциальных энергий для композита и эквивалентной однородной среды.

Формула для расчета эффективного модуля композита по-прежнему имеет вид:

hM hD E hM h f hD h f (25) EM ED Здесь «обменная» объемная доля h f выражена через постоянные интегрирования C D1, C D 2, C M 1, C M 2 следующим образом:

hD C ch D y cos D y C D 2 sh D y sin D y dy hf D 1 P( ) EM E D hD hM C ch M ( y hM hD ) cos M ( y hM hD )dy M 1 1 (26) P( ) hD E M ED hD hM C sh M ( y hM hD ) sin M ( y hM hD ) dy M 1 P( ) hD E M ED На рис. 11 приведены решения для двух композиционных материалов с разными адгезионными параметрами e g 19 ГПа*нм и e g 14,7 ГПа*нм (функционализированные и нефункционализированные нанотрубки).

Остальные параметры композитов зафиксированы: С M 0.0088 ГПа/нм2, g g =inf, E D 435,97 ГПа, E M 0.9 ГПа, f=1%. Также на графике представлены решения для этих композитов, полученные Одегардом и др.

Данные Одегард Осреднение и др.

по Ф ойгту адгезия1, Ev=10.964 адгезия ГПа Осреднение по Ф ойгту Ev=9. ГПа Е - Эффективный модуль Юнга композита, ГПа Решение с 8 учетом когезии и адгезии См=0, ГПа/нм eg= 6 ГПа*нм g=inf См=0, ГПа/нм eg=14. ГПа*нм g=inf Осреднение по Ф ойгту Ev=5. ГПа Решение без учета когезии Осреднение и по Рейссу адгезии Er=0. ГПа 0 50 100 150 200 250 300 350 400 l - длина волокна, нм Рис. 11. Зависимость эффективного модуля Юнга от длины волокна с учетом когезионных взаимодействий при С M от 0,0088 ГПа/нм2 и различных значениях e g Модель показывает хорошее соответствие данным численного эксперимента и позволяет объяснить эффект аномального усиления в нанокомпозитах. Основные выводы по главе:

1. В качестве параметра межфазного слоя используется эффективная толщина межфазного слоя h f h f ( E D, E M, G D, G M, hD, hM, l, С D, С M, g g, e g ).

Данная модель обобщает модели, построенные в предыдущих главах.

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

3. Адгезионные параметры g g и e g могут играть роль параметров качества технологии изготовления композитов.

ОСНОВНЫЕ ВЫВОДЫ ПО РАБОТЕ Наиболее существенные научные результаты и выводы, полученные в диссертации:

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

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

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

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

3. Установлены численные оценки для эффекта изменения эффективного модуля композита в рамках вилки Фойгта-Рейсса с изменением диаметра включения при фиксированной объемной доле и длине волокон. Показано, что эффективный модуль является монотонно убывающей функцией от hD / l. Установлено, что диапазон значений hD / l, при котором эффективный модуль существенно меняется, ограничен. При уменьшении диаметра включений растет относительная длина краевого эффекта в матрице, а значит и толщина межфазного слоя. Указанная зависимость может быть использована при формулировке рекомендаций по получению композитов повышенной жесткости, если удастся снять частично технологические ограничения по толщине используемых волокон.

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

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

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

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

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

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

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

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

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

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

Основное содержание диссертации отражено в публикациях:

1. Гордеев А.В. Моделирование свойств композиционного материала, армированного короткими волокнами // Механика композиционных материалов и конструкций, т. 16, №1, 2010 г.

с.106-116.

2. Гордеев А.В. Моделирование эффектов усиления дисперсных композитов, армированных ориентированнми микро- и нанотрубками // «Механика и наномеханика структурно-сложных и гетерогенных сред. Успехи, проблемы и перспективы». Материалы Всероссийской конференции, приуроченной к 20-летию ИПРИМ РАН. Тезисы докладов. Москва, 30 ноября – 2 декабря 2009 г. М.:Альянстрансатом, 2009, с.95-96.

3. Белов П.А., Гордеев А.В. Моделирование свойств композиционного материала, армированного короткими волокнами.

Учет адгезионных взаимодействий // Композиты и наноструктуры №1, 2010, с.40-46.

4. Соляев Ю., Гордеев А. Моделирование эффектов усиления дисперсных композитов, армированных короткими микро- и нановолокнами с использованием прикладной теории межфазного слоя // «Математическое моделирование в естественных науках».

Тезисы докладов XVIII Всероссийской школы-конференции молодых ученых и студентов. Пермь, издательство Пермского государственного технического университета, 2009, с.94-95.



 

Похожие работы:





 
2013 www.netess.ru - «Бесплатная библиотека авторефератов кандидатских и докторских диссертаций»

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