Моделирование температурного поля и напряженно-деформированного состояния алмазосодержащих инструментальных композитов на полимерной матрице тема автореферата и диссертации по химии, 02.00.06 ВАК РФ

4857586

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

лш

ЯХУТЛОВА Марианна Разиуановна

МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРНОГО ПОЛЯ И НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ АЛМАЗОСОДЕРЖАЩИХ ИНСТРУМЕНТАЛЬНЫХ КОМПОЗИТОВ НА ПОЛИМЕРНОЙ МАТРИЦЕ

Специальность: 02.00.06 - Высокомолекулярные соединения

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

2 О ОПТ 2077

НАЛЬЧИК 2011

Работа выполнена в Кабардино-Балкарском государственном университете имени Х.М. Бербекова.

Научный руководитель: доктор химических наук,

профессор Лигидов Мухамед Хусенович

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

профессор Никитин Лев Николаевич

доктор технических наук, профессор Языев Батыр Меретович

Ведущая организация: Российский химико-технологический

университет имени Д.И. Менделеева

Защита состоится « 3 » ноября 2011 г. в 13.00 часов на заседании дис сертационного совета Д 212.076.09 в Кабардино-Балкарском государствен ном университете им. Х.М. Бербекова по адресу: 360004, г. Нальчик, ул. Чернышевского, 173, КБГУ, ауд. 203.

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

Автореферат разослан » сентября 2011 г.

Учёный секретарь диссертационного совета

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

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

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

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

Поставленная цель определила необходимость решения следующих задач:

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

- разработка математической модели и программно-алгоритмического обеспечения для конечноэлементного моделирования напряженно-деформированного состояния системы алмаз - переходный слой - матрица при силовых и тепловых возмущениях;

- численное моделирование температурного поля системы алмаз - переходный слой - полимерная матрица;

- численное моделирование напряженно-деформированного состояния системы алмаз - переходный слой - полимерная матрица при силовых и тепловых возмущениях;

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

Научная новизна полученных результатов заключается в:

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

- установлении влияния теплопроводности матрицы, переходного слоя и его толщины на температурное поле в системе алмазное зерно - переходный слой - полимерная матрица;

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

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

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

Практическая значимость работы заключается в:

- программно-алгоритмическом обеспечении расчётов температурного поля и напряженно-деформированного состояния в системе алмазное зерно - переходный слой - полимерная матрица при силовых и тепловых возмущениях;

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

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

- Апробация работы. Материалы диссертации доложены:

- на Международных научно-технических конференциях: «Наука, техника и технология XXI века» (Нальчик, 2007 г.); «Компьютерное моделирование 2008» (Санкт-Петербург, 2008 г.); «Новые полимерные композиционные материалы» (Нальчик, 2010,2011 годы);

- на Всероссийских научно-технических конференциях: «Наука, техника и технология нового века» (Нальчик, 2003, 2005 годы); «Новые полимерные композиционные материалы» (Нальчик, 2005 г.).

Публикации. Основные результаты исследований изложены в 10 научных работах, из которых одна - в рецензируемом журнале.

Структура и объем работы. Диссертация состоит из введения и четырёх глав, выводов и списка использованной литературы. Работа изложена на 135 страницах машинописного текста, содержит 54 рисунка, 18 таблиц и список литературы, включающий 90 наименований.

Автор выражает глубокую благодарность доктору физико-математических наук, профессору, заведующему кафедрой вычислительной математики КБГУ Шханукову-Лафишеву Мухамеду Хабаловичу за помощь и научную консультацию при выполнении диссертационной работы.

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

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

Наиболее широкое применение для изготовления алмазных инструментов нашли новолачный фенолформальдегидный олигомер СФ-012, фе-нольное связующее СФП-012А, резольная безаммиачная анилинфенолфор-мальдегидная смола СФ-342А. В связи с появлением новых областей применения алмазного инструмента и повышением требований к его функциональным характеристикам постоянно расширяется круг используемых полимерных связующих. Основные направления развития полимерных связующих связаны с обеспечением большей теплостойкости, механической прочности, хороших антифрикционных свойств, а также с повышением адгезии к поверхности алмаза и наполнителей.

Отмечается лучшее смачивание и более высокая работа адгезии полимера к алмазу, чем к другим материалам, что связывают с наличием на поверхности алмаза гидроксильных, карбонильных, карбоксильных и эфирных групп, с которыми реакционноспособные группы полимеров могут образовывать физические и химические связи. Зависимость краевого угла смачивания от температуры носит экстремальный характер. Температуры, при которых наблюдается наибольшее смачивание, для распространённых полимерных связующих находятся в диапазоне порядка 50-150 °С. Ухудшение смачиваемости при более высоких температурах объясняют интенсивной адсорбцией низкомолекулярных фракций и продуктов термоокисления молекулярных цепей, формированием при этом монослоя с низким поверхностным натяжением.

Анализ литературных данных по формированию структуры и свойств алмазосодержащих инструментальных композитов на полимерной матрице пока-

5

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

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

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

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

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

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

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

Рис. 1. Расчётная схема

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

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

Конструкция может находиться под действием возмущений от процесса резания - статической сосредоточенной нагрузки Г, приложенной под углом (3 к вертикальной оси алмазного зерна, и удельного теплового потока 2, приложенного к свободному контуру зерна, непосредственно контактирую-

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

Расчет поля температур в конструкции вёлся на основе решения нелинейного уравнения нестационарной теплопроводности в плоской изотропной среде

дх

ЦТ)

дТ(х,у,1) дх

ду

Х(Т)

дТ(х,уЛ) ду

-р-с(Т)

дТ(х,у,1) дt

= 0,

где Т(х,у,1) - температура; к(т) - коэффициент теплопроводности; р -плотность; с(т) - удельная теплоемкость; Г - время.

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

1) Если на границе задан поток тепла <3 (<3 считается положительным, если тепло теряется телом):

дтЩ

К \ дх х ду

I

+ 0 = 0.

2) Если на границе происходит конвективный теплообмен, характеризуемый коэффициентом теплоотдачи А:

дТ{х,у,1) 5Т(х,у,()

1-е

дх ду у

+ к[т{х,у,1)-т<а}= 0,

где Т - температура окружающей среды; 1Х - ¡то.:1у = соза - направляющие косинусы.

3) Если граница представляет собой теплоизолированный участок:

дТ{х,у,1)1 | дТ{х,у,1)1 =() дх х ду у

4) Если на некоторой части границы температура известна:

т(х,у,г)=г/5),

где - координаты точек границы.

Начальные условия при решении уравнений нестационарной теплопроводности имеют вид

т(х,у,0) = 0.

Компоненты деформаций и напряжений находятся из соотношений двухмерной задачи теории упругости

ди дх'

ду

ди дv

ау

ду дх

оу =

1-V^

дх ду ду дх

ди Эу

+ ■

где - линейные и угловая деформации; и у - перемещения вдоль

осей х и у соответственно; ох ,иу ,хху - нормальные и касательное напряжения; е,V — модуль упругости и коэффициент Пуассона; О - коэффициент линейного расширения.

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

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

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

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

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

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

9

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

Препроцессор обеспечивает: формирование и визуализацию геометрической модели конструкции и модели внешних воздействий; формирование методических параметров задачи; ввод и диагностику файлов исходных данных; подготовку и выпуск расчётной (табличной и графической) информации; преобразование исходных данных из внешнего представления во внутреннее.

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

При формировании базы данных материалов задаются следующие характеристики: Е[Н/м2] - модуль упругости; V - коэффициент Пуассона; р[кг/мъ] - плотность; [ст] [Н/м2] - допускаемое напряжение; па - число узлов на графике {т+а(т)} (рис. 2), где а [град'1] - коэффициент линейного расширения; Тх,ахХг,аг,—,Т„,(х„ - пары координат узлов на графике {тЧа^)}. Аналогично задаются и зависимости

Рис. 2. График для задания зависимости {т 4- а(т)}

Для формирования файла расчетной схемы конструкции необходимо задать следующие характеристики (см. рис. 1): А к ~ толщина конструкции; а, Ь -полуоси алмазного зерна; ■ иГ- число эллипсов (слоев), на которое разбивается 10

зерно; п2 — число секторов, на которое разбивается четверть зерна; 5 - толщина переходного слоя; Ьк - глубина погружения зерна в матрицу; пЪ - число эллипсов (слоев), на которое разбивается матрица; б - коэффициент разрежения сетки матрицы (1 < е < 2); при е = 1 толщина слоев будет одинакова, при 8 = 2 толщина каждого следующего слоя будет увеличиваться в два раза.

Задаются параметры, определяющие нагрузку на систему: сила /■"[//], угол р , номер узла, к которому приложена сила кр ; удельный тепловой поток 0[Вт/мг] и число граней кд свободного контура зерна, к которым он приложен симметрично (на рис. 1 кд = 3). Также задаются коэффициенты

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

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

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

Главная погрешность решения задач определяется методом конечных элементов, точность которого напрямую зависит от степени дискретизации рассматриваемой конструкции, т.е. от размеров конечных элементов. Поэтому на первом этапе испытаний разработанного программного комплекса исследована сходимость метода конечных элементов по температуре, перемещениям и напряжениям в зависимости от числа конечных элементов. Рассматривались следующие расчётные схемы, соответствующие различным характеристикам разбиения конструкции на конечные элементы: п1 = п2 = пЗ = 6 (488 элементов); п1 = п2 = пЗ = 8 (840 элементов); п1 = п2 = пЗ = 10 (1332 элемента); п1 =п2 = пЗ = 12 (1936 элементов).

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

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

Проверена сходимость указанных выше итерационных процессов. Результаты многочисленных расчётов показали, что число итераций, необходимых для достижения относительной погрешности е = 10'6, заложенной в алгоритмы, не превышает 10-12.

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

1. Использование свойств симметрии рассчитываемой конструкции:

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

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

Отклонения от симметрии в случае а) и от антисимметрии в случае; б) является следствием ошибки во внутреннем построении самого алгоритма.

2. Реакции в опорах: сумма проекций на оси х и у реакций в опорах должна быть равна сумме проекций на оси х и у внешних нагрузок.

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

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

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

Таблица

Параметры модели

Размеры, мм Нагрузка

а b А h Ö F, Н ß, град в, Вт/м2 К=К. Вт /{м1 ■ град)

0,200 0,250 0,200 0,100 0,010 10 30 3,0х106 5x10"

На первом этапе исследований решением уравнения стационарной теплопроводности проведено моделирование влияния свойств матрицы и переходного слоя на температурное поле системы.

На рис. 3 представлено распределение температуры (здесь и далее речь идёт об избыточной температуре - выше 20 °С) в сечении по- вертикальной оси зерна (см. рис. 1) при различных коэффициентах теплопроводности материала матрицы А3. В этих численных экспериментах свойства покрытия принимались идентичными свойствам матрицы. Эксперименты показали, что увеличение теплопроводности матрицы приводит к существенному уменьшению температур в системе. При одинаковой тепловой нагрузке базовому варианту с Къ =0,18 Я/и/( м-град) (использованию бакелитового связующего) соответствует средняя температура в зерне 300 °С, "/.3 = 0,5 - 220°, л3 =1,0 - 170°, л3 =2,0 — 130°. При этом при всех значениях Л3 алмазное зерно прогревается практически равномерно. Небольшой градиент температур в зерне отмечен лишь в сечении вдоль вертикальной оси зерна. Эксперименты показали, что температура в переходном слое при заданном значении Лз меняется несущественно - имеет место её незначительное уменьшение при приближении к охлаждаемой поверхности.

Моделирование влияния теплопроводности переходного слоя на температурное поле показало, что применение наиболее распространенных покрытий - натрийборосиликатного стекла (Х2=1,0 Вт/{м-град)) и никеля (Х2 ~ 65,0), а также титана (Х2 = 15,1), не приводит к существенному уменьшению температуры матрицы (см. рис. 4). Также численные эксперименты показали, что толщина покрытия практически не влияет на температурное поле. Эти результаты моделирования опровергают распространённое мнение о том, что указанные покрытия служат тепловым барьером, резко уменьшающим тепловые потоки в полимерное связующее.

Рис. 3. Изменение температуры в сечении по вертикальной оси зерна при различных коэффициентах теплопроводности матрицы Х3 (Вт/{м-град))

Рис. 4. Изменение температуры в сечении по горизонтальной оси зерна при различных коэффициентах теплопроводности материала покрытия к2 (Вт/(м-град))

На втором этапе исследований проведено моделирование напряжённо-деформированного состояния системы алмаз - переходный слой - полимерная матрица при силовых и температурных возмущениях.

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

25 20 15 10

а5

-5

10 ■15 -20

2' 2 0 4 4'

зона

Рис. 5. Распределение напряжений в сечении по горизонтальной оси зерна при силовом воздействии

Напряжения в контактной области со стороны матрицы увеличиваются по абсолютной величине по мере приближения к областям А и С (см. рис. 1), являющимся зонами концентрации напряжений.

При заданных исходных данных проведены также расчёты напряжённо-деформированного состояния системы при силовых возмущениях для трёхмерной модели с использованием программного комплекса COSMOS

Works. Численные эксперименты показали, что расчётные напряжения, полученные на двух- и трехмерных моделях, различаются не более чем на 5-10 %.

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

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

Как и в случае силового нагружения, максимальная интенсивность напряжений в системе имеет место в зонах концентрации напряжений А и С (см. рис. 1).

Исследовано влияние теплопроводности матрицы Хз на температурные напряжения в системе. Отмечена общая тенденция к повышению интенсивности температурных напряжений в конструкции с уменьшением Х3. При этом зависимости носят нелинейный характер - с возрастанием градиента а, при уменьшении Х3, что особенно выражено при его низких значениях, характерных для полимерных матриц. Анализ результатов показывает, что уровень температурных напряжений в переходном слое существенно зависит от материала покрытия. Наименьший уровень интенсивности напряжений отмечен для покрытия из стекла. Далее следуют в порядке возрастания напряжений бакелит, свинец, медь, титан, никель (см. рис. 7).

Рис. 7. Зависимость интенсивности температурных напряжений в зоне В (см. рис. 1) от теплопроводности матрицы при различных материалах покрытия на зерно

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

а, = 1,05 -10 град' , т. е. случай равенства коэффициентов теплового расширения всех элементов системы. Значениям коэффициентов теплового расширения бакелитовой матрицы а3 =70-10~6 град'1 соответствуют оптимальные значения а2 = (1-П.5,) Ю-6 град'1.

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

3

Рис. 8. Зависимость интенсивности напряжений в зоне В (см. рис. 1) от коэффициентов теплового расширения переходного слоя а2 и матрицы а3

Рис. 9. Распределение интенсивности напряжений в переходном слое при действии нагрузки: силовой (Р)\ тепловой 0 = 3 ■ 10б (вт / м2); силовой и тепловой (7^0; тепловой 0 = 2Ю6 ; тепловой = 106

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

Как видно из рис. 10, время переходного теплового процесса практически не зависит от теплопроводности матрицы. Она определяет установившуюся температуру в системе. Для указанной зоны В' (см. рис. 1), максимальная установившаяся температура, равная примерно 270 °С имеет место при использовании бакелитового связующего (Х3= 0,18 Вт Ам • градI).

I, с

Рис. 10. Нестационарный тепловой процесс в зоне В'(см. рис. 1) при различных коэффициентах теплопроводности матрицы Х3 (Вт/[м -град))

Удельная теплоёмкость матрицы с3 не влияет на уровень установившейся температуры в системе, но влияет на время переходного теплового процесса. При заданных в численных экспериментах значениях удельной теплоёмкости матрицы с3 от 100 до 2400 Дж/(кг-град) стационарная температура устанавливается за

время от 4 до 16 секунд. При этом для бакелитового связующего (с3-1600 Дж/(кг ■ град)) время переходного процесса составляет 11 секунд.

Коэффициент теплоотдачи Иа = Ит во внешнюю среду существенно влияет на переходный тепловой процесс в исследуемой системе (См. рис. 11). С его увеличением снижается уровень установившейся температуры. Время переходного процесса увеличивается с уменьшением коэффициента теплоотдачи, оно составляло в экспериментах от 5 до 104 секунд. Расчёты показали,

что при базовом значении й0 =йт= 5 ■ 103 Вт/{м2 -град) установившаяся температура 280 °С достигается за 14 секунд.

Рис. 11. Нестационарный тепловой процесс в зоне В' (см. рис. 1) при различных коэффициентах теплоотдачи

в окружающую среду К =кт (Вт/[и2 -град))

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

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

выводы

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

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

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

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

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

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

7. Исследования показали, что нестационарный тепловой процесс в системе алмаз — переходный слой - полимерная матрица определяется, в основном, теплопроводностью и теплоёмкостью матрицы, а также теплоотдачей во внешнюю среду. Определены зависимости установившейся температуры и времени переходного теплового процесса от этих параметров.

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

1. Карданова (Яхутлова) М.Р. Задачи стационарной и нестационарной теплопроводности для трехслойной конструкции // Матер. Всерос. научн.-техн. конф.: Наука, техника и технология нового века (НТТ-2003). - Нальчик,

2003.-С. 76-81.

2. Карданова М.Р. Решение двумерной динамической задачи теории упругости методом конечных элементов // Информационные технологии в промышленности и экономике: сборник научных трудов Института конст-рукторско-технологической информатики РАН. Вып. 4 / под ред. Ю.М. Со-ломенцева. - М„ 2004. - С. 235-237.

3. Яхутлов М.М., Батыров У.Д., Карданова М.Р. Моделирование напряжённо-деформированного состояния абразивного инструмента при динамическом нагружении // Информационные технологии в промышленности и экономике: сборник научных трудов Института конструкторско-техноло-гической информатики РАН. Вып. 4 / под ред. Ю.М. Соломенцева. - М.,

2004.-С. 238-241.

4. Карданова М.Р. Математическая модель абразивного инструмента на органической матрице // Матер. II Всерос. научн.-практ. конф.: Новые полимерные композиционные материалы. - Нальчик, 2005. - С. 307-311.

5. Карданова М.Р. Алгоритм решения нелинейной нестационарной задачи термоупругости для трехслойного композита // Матер. II Всерос. научн.-техн. конф.: Наука, техника и технология XXI века (НТТ-2005). - Нальчик,

2005. - С. 54-58.

6. Карданова М.Р. О физико-механических свойствах элементов алмазного инструмента на полимерной связке // Матер. III межд. научн.-техн. конф.: Наука, техника и технология XXI века (НТТ-2007). Т. 1. - Нальчик, 2007.-С. 157-159.

7. Яхутлов М.М., Батыров У.Д., Карданова М.Р. Программно-алгоритмическое обеспечение моделирования напряжённо-деформированного состояния в системе алмаз-матрица // Труды межд. научн.-техн. конф.: Компьютерное моделирование 2008. - СПб., 2008. - С. 272-275.

8. Яхутлов М.М., Карамурзов Б.С., Батыров У.Д., Ошхунов М.М., Карданова М.Р. Напряжённо-деформированное состояние системы зерно-матрица алмазного инструмента при силовых возмущениях // Сверхтвердые материалы. - Киев, 2009. -№ 6. - С. 81-88.

9. Карданова М.Р., Лигцдов MX, Шхануков-Лафишев М.Х., Яхутлов М.М. Тепловой режим в алмазосодержащем инструментальном композите на полимерном связующем // Матер. VI Межд. научн.-практ. конф.: Новые полимерные композиционные материалы. - Нальчик, 2010. - С. 178-185.

10. Карданова М.Р., Лигвдов М.Х., Шхануков-Лафишев МЛ., Яхутлов М.М. Моделирование температурного поля в композиционном алмазосодержащем материале на органической матрице // Пласт, массы. - 2010. - № 5. - С. 34-3 8. 22

В печать 28.09.2011. Тираж 100 экз. Заказ № 6379. Полиграфический участок ИПЦ КБГУ 360004, г. Нальчик, ул. Чернышевского, 173.

ВВЕДЕНИЕ.

1. АНАЛИЗ СОСТОЯНИЯ ВОПРОСА И ПОСТАНОВКА

ЦЕЛИ И ЗАДАЧ РАБОТЫ.

1.1. Алмазосодержащие инструментальные композиты на полимерной матрице.

1.2. Структура и свойства полимерных связующих.

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

1.4. Формирование структуры и свойств композитов при прессовании и термообработке.

1.5. Моделирование температурного поля и напряжённо-деформированного состояния алмазосодержащих композитов.

1.6. Цель и задачи работы.

2. РАЗРАБОТКА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ

И АЛГОРИТМОВ РАСЧЁТА ТЕМПЕРАТУРНОГО ПОЛЯ И НАПРЯЖЁННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ.

2.1. Уравнение нестационарной теплопроводности.

Постановка краевых задач.

2.2. Конечно-элементная формулировка задачи.

2.3. Постановка задачи термоупругости и её конечно-элементная формулировка.

2.4. Алгоритмы решения задач стационарной и нестационарной термоупруго сти.

2.5. Выводы.

3. РАЗРАБОТКА ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ МОДЕЛИРОВАНИЯ И ОБОСНОВАНИЕ

ДОСТОВЕРНОСТИ АЛГОРИТМОВ.

3.1. Общие принципы построения программного комплекса.

3.2. Принципы построения препроцессора и подготовка исходных данных для расчёта.

3.3. Обоснование достоверности алгоритмов.

3.4. Выводы.

4. МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРНОГО ПОЛЯ И НАПРЯЖЁННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ СИСТЕМЫ АЛМАЗ - ПЕРЕХОДНЫЙ СЛОЙ - ПОЛИМЕРНАЯ МАТРИЦА.

4.1. Расчётная схема и свойства элементов.

4.2. Моделирование температурного поля в системе алмаз - переходный слой - полимерная матрица.

4.3. Моделирование напряжённо-деформированного состояния системы алмаз - переходный слой - полимерная матрица.

4.4. Моделирование нестационарного теплового процесса в системе.

4.5. Выводы.

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

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

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

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

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

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

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

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

Для достижения поставленной цели в диссертационной работе решены следующие задачи:

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

2. Разработка математической модели и программно-алгоритмического обеспечения для конечноэлементного моделирования напряжённо-деформированного состояния системы алмаз - переходный слой - матрица при силовых и тепловых возмущениях.

3. Численное моделирование температурного поля системы алмаз - переходный слой - полимерная,матрица.

4. Численное моделирование напряжённо-деформированного состояния системы алмаз - переходный слой - полимерная матрица при силовых и тепловых возмущениях.

5. Численное моделирование нестационарного теплового процесса в системе алмаз - переходный слой - полимерная матрица.

Автор защищает:

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

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

- результаты численного моделирования стационарного температурного поля системы алмаз - переходный слой - полимерная матрица; результаты численного моделирования напряжённо-деформированного состояния системы алмаз - переходный слой - полимерная матрица при действии силовых и температурных возмущений;

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

Диссертационная работа выполнена на кафедре "Органическая химия и высокомолекулярные соединения» Кабардино-Балкарского государственного университета.

ОСНОВНЫЕ ВЫВОДЫ ПО РАБОТЕ

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

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

3. Численное моделирование1 стационарного-температурного поля системы показало, что оно существенно зависит от теплопроводности матрицы и практически^ не зависит от толщины и теплопроводности переходного слоя: Отмечено; что из-за низкой» теплопроводности бакелитовой матрицы алмазное зерно в процессе работы прогревается практически равномерно.

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

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

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

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

1. Абразивная и алмазная обработка материала. Справочник. / Под ред. А.Н. Резникова. -М.Машиностроение, 1977. -392 с.

2. Александров В.А., Жуковский А.Н., В.П., Мечник В.А. Температурное поле, термоупругое состояние и износ алмазного круга при резании с охлаждением. 4.11 // Трение и износ. —1991. —№3.

3. Александров В.А., Жуковский А.Н., Карагодов В.П., Мечник В.А. Расчёт нестационарного температурного поля алмазного круга при резании с охлаждением // ИФЖ. 1989. -Т.56. -№4. -С. 690-691.

4. Александров В.А., Мечник В.А., Верхоярный A.B. Изучение нестационарного температурного поля алмазного круга при резании с охлаждением. // Сверхтвёрдые материалы. -1989. -№1. -С. 40-45.

5. Алешин В.Г., Смехнов A.A., Богатырёва Г.П., Крук Б.Б. Химия поверхности алмаза. -Киев:Наук. думка, 1990. -200 с.

6. Аммерал JL Принципы программирования в машинной графике: Пер. с англ. -М.:Сол Систем, 1992. -224 с.

7. Байкалов А.К. Введение в теорию шлифования материалов. -Киев: Наукова думка, 1978. -207 с.

8. Баничук Н.В., Кобелев В.В., Рикардс Р.Б. Оптимизация элементов конструкций из композиционных материалов. -М.Машиностроение, 1998. —224 с.

9. Березин И.С., Жидков Н.П. Методы вычислений. В 2-х т. -М.:Физматгиз, 1962. Т.1. -464 с.

10. Берлин A.A., Басин В.Е. Основы адгезии полимеров. -М.:Химия, 1974. -392 с.

11. Богатырёва Г.П. Исследование гидрофильности и гидрофобности поверхности синтетических алмазов // Сверхтвёрдые материалы, 1980, №2. С. 23-27.

12. Богатырёва Г.П., Невструев Г.Ф., Ильницкая Г.Д., Коновалов В.А., Ткач В.Н. Возможность повышения прочности удержания алмазов в связке // Сверхтвёрдые материалы, 2001, №2. С. 21-25.

13. Брандо» Д., Каплан У. Микроструктура материалов. Методы исследования и контроля. -М.¡Техносфера, 2006. -384 с.

14. Васильев Л.А., Белых З.П. Алмазы, их свойства и применение. -М.:Недра, -1983. -102 с.

15. Воронин Г.А., Шило А.Е. Термические напряжения в материалах на основе абразива и связующего. // Сверхтвёрдые материалы. — 1982. —№5. -С. 19-22.

16. Годовский Ю.К. В кн. Энциклопедия полимеров, т.З. М.¡Советская энциклопедия, 1977. -С. 599.

17. Джордж А., Лю Дж. Численное решение больших разреженных систем уравнений: Пер. с англ. -М.:Мир, 1984. -334 с.

18. Дувакина Н.И., Ткачёва Н.И. Выбор наполнителей для придания специальных свойств полимерным материалам // Пластические массы, 1989, №11.-С. 46-48.

19. Журавлёв В.В. Влияние металлизации на прочность алмаза и величину внутренних напряжений системы алмаз-металл // В сб. Повышение эффективности применения алмазных инструментов. Труды ВНИИалмаза. М., -1986. -С. 50-56.

20. Захаренко И.П., Ахундов Э.А. Об устойчивости зёрен в связке шлифовального круга. // Синтетические алмазы. -1978. -Вып.6. -С. 24-28.

21. Захидов С.И. Исследование прочности удержания зерна в связке при температурно-силовых воздействиях. Автореф. дисс. . канд. техн. наук. -М, 1973.

22. Зенкевич О. Метод конечных элементов в технике: Пер. с англ. -М.: Мир, 1975. -544 с.

23. Зубчанинов В.Г. Основы теории упругости и пластичности. —М.:Высш. шк., 1990.-368 с.

24. Инструменты из сверхтвёрдых материалов / Под ред. Н.В. Новикова. -М.Машиностроение, 2005. -555 с.

25. Канторович Л.В., Акилов Г.Р. Функциональный анализ в нормированных пространствах. -М.:Физматгиз, 1959. -684 с.

26. Киреев В.В .Высокомолекулярные соединения. -М.:Высш. шк., 1992. -512 с.

27. I-Сноп А., Шейб В., Фенольные смолы и материалы на их основе. Пер. с англ. -М.:Мир, 1983. -320 с.

28. Кныш C.B., Склепчук В.А. Выбор формы зерна при моделировании процессов шлифования. // Резание и инструмент. -1988. -Вып. 39. -С. 95-98.

29. Композиционные материалы. Справочник / Под ред. Д.И. Карпинос. -Киев:Наук. думка, 1985. -592 с.

30. Композиционные материалы: В 8-ми Т. Пер. с англ. / Под ред. JI. Бра-утмана и Р. Крока. Т.4 Композиционные материалы с металлической матрицей / Под ред. К.Крейдера. -М.Машиностроение, 1978. -503 с.

31. Корн Г., Корн Т. Справочник по математике: Пер. с англ. -М.:Наука, 1977.-836 с.

32. Курдюков В.И. Научные основы проектирования абразивного инструмента. -Курган:Изд-во Курганского гос. ун-та, 2005. -159 с.

33. Кущ В.И., Шило А.Е., Чалый В.Т. Математическая модель теплопере-носа в системе зерно-покрытие-полимерная связка // Сверхтвёрдые материалы. -2003. —№5. -С. 49-59.

34. Лавриненко В.И., Кулаковский В.Н., Ломашевская Н.В. и др. Напряженное состояние в зоне взаимодействия зерна со связкой круга // Сверхтвёрдые материалы. -1995. -№4 -С. 46 49.

35. Лавриненко В.И., Шепелев A.A., Петасюк Г.А. Модели формы зерен СТМ // Сверхтвёрдые материалы. -1994. -№5-6. -С. 18-21.

36. Липатов-Ю.С. Физическая химия наполненных полимеров. -М.:Химия, 1977.-304 с.

37. Лошак М.Г., Полотняк С.Б., Александров Л.И. Численное моделирование напряжённо-деформированного состояния вольфрамовых твёрдых сплавов после спекания // Сверхтвёрдые материалы. -2005. -№4. -С. 30-40.

38. Лыков A.B. Теория теплопроводности. —М.:Высш. шк., 1967. -600 с.

39. Малышев В.И., Янюшкин Ю.М. Контактные температуры при алмазной правке шлифовальных кругов // Сверхтвёрдые материалы. -1986. —№5. -С. 48-54.

40. Микитаев А.К., Козлов Г.В. Фрактальная механика полимерных материалов. -Нальчик:Каб.-Балк. ун-т, 2008. -312 с.

41. Михеев М.А., Михеева И.М. Основы теплопередачи. -М.:Энергия, 1977.-342 с.

42. Михлин С.Г. Вариационные методы в математической физике. -М.:Наука, 1970.-512 с.

43. Мяченков В.И., Ольшанская Г.Н., Чеканин A.B. Автоматизация конструирования и прочностных расчётов тонкостенных осесимметричных конструкций: Общее описание. -М.:МГТУ "Станкин", 1994. -64 с.

44. Мяченков В.И., Ольшанская Т.Н., Чеканин A.B. Автоматизация конструирования и прочностных расчётов тонкостенных осесимметричных конструкций. KIPR-IBM-PC/AT 2.0: Формирование расчётных схем. -М.:МГТУ "Станкин", 1994. -64 с.

45. Мяченков В.И., Ольшанская Г.Н., Чеканин A.B. Автоматизация конструирования и прочностных расчётов тонкостенных осесимметричных конструкций. KIPR-IBM-PC/AT 2.0: Технология работы с системой. -М.:МГТУ "Станкин", 1994.-96 с.

46. Найдич Ф.В., Колесниченко Г.А., Лавриненко И.А., Моцак Я.Ф. Пайка и металлизация сверхтвёрдых инструментальных материалов. —Киев:Наук. думка, 1977. —187 с.

47. Наполнители для полимерных композиционных материалов. Справочное пособие. Под ред. Г.С. Каца. Пер. с англ. Под ред. П.Г. Бабаевского. —М.:Химия, 1981.

48. Никулин Н.И. Исследование сил, возникающих при микрорезании хрупких материалов // Синтетические алмазы. -1978. №5. -С. 52-57.

49. Нильсен Л. Механические свойства полимеров и полимерных композиций. М.:Химия, 1978. -312 с.

50. Носов Н.В. Повышение эффективности и качества абразивных инструментов путём направленного регулирования их функциональных показателей. Автореф. дисс. докт. техн. наук. -Самара, 1997. -46 с.

51. Олейников А.Б. Работоспособность алмазных эластичных кругов при шлифовании газотермических покрытий. Автореф. дисс. . канд. техн. наук. -Киев, 1990. .

52. Олейников А.Б., Сенченков И.К., Рубцова И.Г. Влияние напряжённо-деформированного состояния контакта зерно-связка на работоспособность кругов с режущим слоем из АЛШЛ // Сверхтвёрдые материалы. -1987. —№5. -С. 45-49.

53. Определение контактной температуры при правке абразивных кругов алмазным инструментом. Х.Г. Тхагапсоев, М.Х. Шхануков, Б.С. Хапачев, М.Х. Абрегов // Сверхтвёрдые материалы. -1983. -№4. -С. 44-48.

54. Основы теплофизики и реофизики полимерных материалов / Привалко В.П., Новиков. В.В., Яновский Ю.Г. ; Отв. ред. Романкевич О.В. -Киев:Наук. думка, 1991.-232 с.

55. Петаскж О.У. О поведении алмазного зерна в связке под действием нагрузки. // В сб. Исследование и применение сверхтвёрдых и тугоплавких материалов.-Киев:ИСМ АН УССР, 1981.-С. 106-110.

56. Поляков В.П., Ножкина A.B., Чириков Н.В. Алмазы и сверхтвёрдые материалы. -М. Металлургия, 1990. -327 с.

57. Применение метода конечных элементов к расчёту конструкций / P.A. Хечумов, X. Кепплер, В.И. Прокопьев; Под общ. ред. P.A. Хечумова. -М. ¡Издательство АСВ, 1994. -353 с.

58. Расчёты машиностроительных конструкций методом конечных элементов: Справочник / В.И.Мяченков, В.П.Мальцев, В.П.Майборода и др.; Под общ. ред. В.И. Мяченкова. -М. Машиностроение, 1989. —520 с.

59. Рындин Н.И. Краткий курс теории упругости и пластичности / Под ред. B.C. Постоева. -Л.:Изд-во Ленингр. ун-та, 1974. -136 с.

60. Сагарда A.A., Химач О.В. Контактная температура и силовые зависимости при резании алмазным зерном. // Синтетические алмазы. -1972. —№2. -С.5-9.

61. Сегерлинд Л. Применение метода конечных элементов. -М.:Мир, 1979. -392 с.

62. Сендецкий Дж. Механика композиционных материалов. -М.:Мир,1970.-308 с.

63. Сердюк В.М., Коновалов В.А., Чалый В.Т. Исследование прочности удержания алмазных зёрен в органической связке // Синтетические алмазы,1971, вып. 4. -С. 33-35.

64. Синтетические сверхтвёрдые материалы: В 3-х т. Т.2, Композиционные инструментальные сверхтвёрдые материалы. / Под ред. Н.В. Новикова. -Киев:Наук.думка, 1986. -264 с.

65. Соложенко В.JI., Дуб С.Н., Новиков Н.В. К вопросу о твёрдости кубического карбонитрида бора // Сверхтвёрдые материалы. -2001. -№4. -С. 7378.

66. Справочник по алмазной обработке металлорежущего инструмента / Бакуль В.Н., Захаренко И.П., Кунькин Я.А., Мильштейн М.З. Под общей редакцией Бакуль В.Н. Киев:Техника, 1971. —208 с.

67. Тимошенко С.П., Гудьер Дж. Теория упругости. -М.:Наука, 1979: -560 с.

68. Тихонов A.B., Самарский A.A. Уравнения математической физики. -М.:Наука, 1977. -736 с.

69. Химач О.В., Ковыженко Г.И., Колмогоров П:В. Контактная температура при микрорезании твёрдого сплава ВК8. // Сверхтвёрдые материалы. -1981. -№2. -С. 59-61.

70. Чалый1 В.Т. Рациональные режимы прессования алмазоносного слоя инструмента на органической связке. В кн.: Полимерабразивные технологические материалы и инструменты в металлообработке. —Киев:Наук. думка, 1981.-С. 57-62.

71. Чалый В.Т., Гороховский Г.А., Малоголовец В.Г. и др. Раздельная термическая обработка алмазного инструмента на органической связке. В кн.: Синтетические сверхтвёрдые материалы и твёрдые сплавы. -Киев:ИСМ АН УССР, 1973. -С. 47-59.

72. Чалый В.Т., Карпович Н.С., Кислый П.С. и др. Алмазный хонинговаль-ный инструмент для финишной обработки гильз. // Сверхтвёрдые материалы, 1980.-№6. -С. 26-31.

73. Чистяков Е.М., Коробко В.Р., Мазур К.И. Влияние металлизации на напряжённо-деформированное состояние алмазоносного слоя инструмента // Сверхтвёрдые материалы. -1989. -№4. -С. 30-34.

74. Чистяков Е.М., Шепелев A.A., Дуда Т.М., Черных В.П. -Инструмент из металлизированных сверхтвёрдых материалов. Киев: Наук, думка, 1982. -204 с.

75. Шило А.Е., Кущ В.И., Дудка В.А. и др. Анализ теплового режима рабочей поверхности однослойного инструмента из СТМ методом конечных элементов. // Сверхтвёрдые материалы. -1989. -№5. -С. 38-41.

76. Яхутлов. М.М. Исследование напряжённо-деформированного состояния системы зерно связка алмазных инструментов // Станки и инструмент. -2001.-№11.-С 23-26.

77. Яхутлов М.М. Исследование теплового режима в системе зерно-матрица алмазного инструмента // Вестник машиностроения. -2001. -№8. -С 48-52.

78. Яхутлов М.М. Повышение работоспособности алмазных инструментов направленным изменением физических характеристик их режущей части. Автореф. дисс. докт. техн. наук. -Москва, 2001. —38 с.

79. Дутка В.А., Колодніцький В.М., Заболотний С.Д., Свешиніков І.А., Лукаш В.А. Моделювання рівня температури в породоруйнівних елементах бурових коронок // Сверхтвёрдые материалы. -2004. -№2. -С. 66-73.

80. Дутка И.А., Колодніцький В.М., Мельничук О.В., Заболотный С.Д. Математична модель теплофізичних процесів при взаимодіі породоруйнівних елементів бурових коронок з массивом гірськоі породи // Сверхтвёрдые материалы. -2005. -№1. -С. 67-77.

81. Klimenko S.A., Mukovoz Yu.A., Polonsky L.G. Ch. 1. Cutting Tools of Su-perhard Materials // Advanced Ceramics Tool for Machining Application-2 / Ed. By I. M. Low and X. S. Li. Switzerland: Trans Tech Publications, 1996. P. 1-66.

82. New surface engineering techniques can make metals and plastics more resistant to heat and corrosion. Chem. Eng. 1994. April. 35.

83. Structure & Properties of Shock-Wave Sintered Diamond Composites / N. Novikov, V. Trefilov, A. Maystrenko, V. Kovtun // Ind. Diamond Rev. 1993. 53. №5. P. 278-281.

84. Tanaka T. Thermal erosion and wetting of diamond coated with pure metals // Bull. Jap. Soc. Prec. Eng. -1980. -14, N2. -P. 107-108.

85. Yamamoto Y., Horike M., Hoshina N., Kabayashi A. A study on the temperature variation of workpieces during cylindrical plung grinding process. -Annals of the CIRP, 1977. -Vol. 26/1. -P. 151-154.

86. Yang Dong Y., Seired Ali A. Model for predicting residual stresses in metal cutting // Proc. Jap. Int. Tribol. Conf. Nagoya, Oct.29-Nov. I. 1990,-Tokyo, 1990.-Vol. I.-P. 439-444.