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

Яковлев Максим Григорьевич

ПОВЫШЕНИЕ ПРОИЗВОДИТЕЛЬНОСТИ ТОКАРНОЙ ОБРАБОТКИ

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

МОДЕЛИРОВАНИЯ ДИНАМИКИ ПРОЦЕССА РЕЗАНИЯ

Специальность:

05.03.01 Технологии и оборудование механической

и физико-технической обработки.

АВТОРЕФЕРАТ

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

Москва, 2009

Работа выполнена на кафедре «Высокоэффективные технологии обработки» Государственного образовательного учреждения Московский государственный технологический университет «СТАНКИН».

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

Официальные оппоненты: доктор технических наук, доцент Воронов Сергей Александрович кандидат технических наук, доцент Максимов Анатолий Дмитриевич Ведущее предприятие: ОАО «Национальный институт авиационных технологий» (НИАТ), г. Москва.

Защита состоится « » 2009 года в часов на заседании диссертационного совета Д212.142.01 в ГОУ ВПО МГТУ «СТАНКИН» по адресу: 127994, ГСП-4, Москва, Вадковский пер., 3а.

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

С диссертацией можно ознакомиться в библиотеке ГОУ МГТУ «СТАНКИН».

Автореферат разослан « » 2009 года.

Ученый секретарь диссертационного совета, кандидат технических наук Волосова М.А.

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

Актуальность вопроса. Эффективность обработки деталей из жаропрочных материалов резанием существенно зависит от динамического поведения технологической системы “станок-приспособление-инструментдеталь” (ТС). Эта зависимость подтверждается многочисленными исследованиями и практическим опытом технологов предприятий. Однако отсутствие количественных зависимостей между динамическими характеристиками технологической системы и изнашиванием инструмента не позволяет назначать наиболее эффективные режимы резания для конкретных технологических процессов. Указанная проблема машиностроительного производства в большей степени сказывается на операциях механической обработки с применением маложесткого инструмента: глубокого сверления, растачивания, концевого фрезерования и др. Вместе с тем особенно актуальна эта проблема в технологии авиадвигателестроения, где широко применяются тонкостенные детали из жаропрочных труднообрабатываемых материалов. В процессе токарной обработки таких деталей возникают повышенные силы резания и автоколебания технологической системы, приводящие к интенсивному изнашиванию, сколам и поломкам инструмента и, как следствие, к снижению производительности и качества обработки.

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

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

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

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

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

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

- разработка методики определения стойкости инструмента на базе термо – силовой модели расчета стойкости для процессов токарной обработки с колебаниями.

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

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

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

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

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

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

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

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

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

Апробация работы. Основные результаты работы докладывались и обсуждались на следующих научно-технических конференциях: “Актуальные проблемы Российской космонавтики” (ХХХI чтения по космонавтике, Москва, МГТУ им. Н.Э. Баумана, 2007), “Новые материалы и технологии НМТ-2008” (Москва, МАТИ, 2008), “Будущее машиностроения России” (Москва, МГТУ им. Н.Э. Баумана, 2008).

Публикации. По теме диссертации опубликовано 7 печатных работ.

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

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

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

В современном авиадвигателестроении широко применяются детали из специальных жаропрочных сплавов, относящихся к группе труднообрабатываемых материалов. Трудоемкость их обработки связана не только с высокими температурами и напряжениями при резании, но и с большими вибрациями инструмента относительно обрабатываемых поверхностей детали. Многочисленные исследования различных авторов свидетельствуют о том, что вибрации, возникающие при резании, оказывают значительное влияние на изнашивание инструмента и качество обработанной поверхности. Исследованиями ученых: Подураева В.Н., Маркова А.И., Ярославцева В.М., Туктанова А.Г. и др., установлено положительное влияние низкочастотных (до 200Гц) и высокочастотных (от 20кГц) вынужденных колебаний на интенсификацию и повышение качества обработки резанием труднообрабатываемых материалов. С другой стороны в работах Жаркова И.Г., Кондратова А.С., Соколовского А.П., Решетова Д.Н. и др. показано, что автоколебания от 0,3 до 10кГц понижают производительность и качество обработки резанием. Проблеме, связанной со снижением вибраций при обработке резанием, посвящены научные работы многих авторов: Каширина А.И., Кудинова В.А., Мурашкина Л.С., Козочкина М.П., Altintatas Y., Budak E., Armarego E.J.A. и др.

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

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

Разработанная модель может описывать перемещения режущей кромки инструмента относительно детали в двух взаимно перпендикулярных направлениях, совпадающих с действием составляющих силы резания (рис.1).

Рис.1. Схема операции обработки диска компрессора высокого давления а) и расчетная динамическая схема при резании б).

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

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

где Kp,x,z – эмпирический коэффициент резания, b – глубина резания.

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

Так как в рассматриваемом случае инструмент колеблется и в направлении Z, то T не равно времени оборота То, а отличается на величину t, которая появляется за счет вибрационной составляющей перемещения режущий кромки b в направлении оси Z за время оборота. Величина b будет равняться разности перемещений на этом обороте Z(t) и на предыдущем Z(t-Tо), а t определяется, как:

где R-радиус заготовки, no-частота вращения заготовки.

Тогда время запаздывания T в плоскости Х с учетом вибрационной составляющей в плоскости Z будет:

а составляющие силы резания для направлений X и Z запишутся как:

В процессе резания всегда присутствует сила трения, которая возникает составляющие силы трения в направлениях X и Z зависят от относительной скорости скольжения между инструментом и заготовкой и описываются нелинейной функцией, например полиномом третьей степени:

где -эмпирический коэффициент, (мин/м)2.

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

где Mx, Mz- обобщенные массы колебательной системы, x, z - скорость перемещений колебательной системы, cx, cz – обобщенные коэффициенты трения в системе (коэффициент демпфирования), kx, kz – обобщенные коэффициенты жесткости систем.

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

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

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

где ki = hi, hi – износ по задней грани, [hз]– допустимый износ по задней грани, G, r, a - эмпирические коэффициенты.

Так же во второй главе предложено в термо-силовой модели за оценку критерия износа инструмента взять допустимую силу резания [ Pд ], а не износ инструмента по задней грани. Так как силы резания можно измерять в процессе резания, а для измерения износа по задней грани необходимо прерывать процесс обработки. Тогда модель (7) будет иметь такой вид:

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

В предложенной методике для исследования ТС была разработана линейная модель (9), учитывающая регенеративный механизм возбуждении колебаний при резании по следу поверхности обработанной на предыдущем обороте заготовки и позволяющая строить диаграммы виброустойчивых режимов резания. Динамические параметры технологической системы Mx,z, cx,z, kx,z и коэффициенты модели силы резания Kp идентифицировались эмпирическим путем. По полученным значениям Mx,z, cx,z, kx,z, Kp рассчитывались частотные характеристики колебаний технологической системы в процессе резания. В нелинейной модели (6) учитывалась зависимость силы трения от скорости резания с помощью эмпирического фрикционного коэффициента. Данная методика позволила рассчитать параметры колебаний режущего инструмента относительно обрабатываемой поверхности в зависимости от режимов резания, жесткостных характеристик технологической системы и свойств обрабатываемого материала.

Экспериментальные исследования проводились при продольном точении жаропрочного никелевого сплава ХН73МБТЮ (ЭИ698). В качестве инструмента применялись токарные резцы, оснащенные твердым сплавом ВК10хом. Изменяя жесткость державки, варьировали динамическими параметрами данной технологической системы.

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

Методы определения коэффициентов Kx, Kz модели сил резания разработаны на базе экспериментов по прямоугольному свободному резанию исследуемого обрабатываемого материала при точении детали с постоянной специальные образцы трубчатой формы. В процессе точения измерялись значения составляющих силы резания в направлении скорости резания Pz, и в направлении подачи Px при различных значениях подачи s0. Из полученных результатов были найдены коэффициенты модели сил резания:

Для идентификации параметров Mx,z, cx,z, kx,z динамической модели тестируемой технологической системы с помощью виброударного молоточка и измерении возникающих вибраций. Полученные в результате измерений данные анализировались с помощью специальных программ обработки сигналов, которые используют Фурье-преобразование и описаны в работе.

динамические характеристики нашей системы:

ТС№1: Mx=0.074 кг, kx=7.74*107Н/м, cx=3%, Mz=0.57 кг, kz=5.9*108Н/м, cz=4%.

ТС№2: Mx=0.063 кг, kx=1.3*107Н/м, cx=9%, Mz=0.095 кг, kz=1.7*107Н/м, cz=9%.

ТС№3: Mx=0.063 кг, kx=1*107Н/м, cx=14%, Mz=0.096 кг, kz=1.67*107Н/м, cz=5%.

ТС№4: Mx=0.1 кг, kx=8.4*106Н/м, cx=7.9%, Mz=0.087 кг, kz=6.7*106Н/м, cz=7%.

ТС№5: Mx=0.134 кг, kx=5.5*106Н/м, cx=7%, Mz=0.092 кг, kz=4.1*106Н/м, cz=7%.

Используя полученные коэффициенты: Mx,z, kx,z, cx,z, ТС №2, были сделаны расчеты перемещения режущей кромки инструмента на различных режимах резания. Пример полученной зависимости приведен на рис.2.

Рис.2. График изменения рассчитанных перемещений режущей кромки инструмента во времени при точении сплава ХН73МБТЮ на режиме резания b =1 мм, V = 20 м/мин, s0 = 0,1 мм/об, для ТС №2.

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

Для идентификации параметров нелинейной модели (6) установим связь между скоростью резания и амплитудой перемещений кромки инструмента через эмпирический фрикционный коэффициент. Результаты испытания с различной скоростью резания 10-30 м/мин и замером амплитуды вибраций кромки инструмента изображены на рис.3.

Амплитуда, мкм С использованием полученных экспериментальных данных были подобраны значения коэффициента. Зависимость (V) аппроксимировали экспоненциальной функцией = 10325 e0,2525V.

колебания технологической системы в двух взаимно перпендикулярных направлениях при заданных жесткостных характеристиках технологической системы: Mx,z, kx,z, cx,z на различных режимах резания при токарной обработке.

Рассчитанные амплитудно-частотные характеристики процесса резания при изменении жесткостных характеристик технологической системы (№1,2,3,4,5) на режиме резания b =1 мм, V = 28 м/мин, s0 = 0,12 мм/об, представлены на рис.4.

инструмента в процессе токарной обработки с учетом вибраций.

Для этого сначала были определены коэффициенты термо-силовой модели стойкости инструмента: G = 3.1 1014, r=4.1, a=0.0332.

жесткостных характеристик технологической системы (№1,2,3,4,5), вычисленные ранее, рассчитаем динамические значений ширины, толщины срезаемого слоя и относительную скорость движения между инструментом и заготовкой. Ширина резания и толщина меняются незначительно для всех пяти случаев, поэтому их изменениями можно пренебречь, а изменения скорости для пяти ТС будут следующими:

ТС№1 28±3 м/мин, ТС№2 28±5 м/мин, ТС№3 28±6 м/мин, ТС№4 28± м/мин, ТС№5 28±10 м/мин.

Подставляя эти значения в формулу (7), получим следующие значения стойкости инструмента для пяти ТС: ТС№1 21 мин., ТС№2 20,5 мин., ТС№ 19 мин., ТС№4 16 мин., ТС№5 11 м/мин.

диагностическом стенде, разработанном на базе токарного станка фирмы Jesco Machinery 1650ENC с ЧПУ системы Fagor 800TGI, укомплектованного трехкомпонентным динамометром «Kistler 9257ВА», акселерометром «Kistler 8614A1000M1» с усилителем и фильтром сигнала «Kistler 5127В» и трех компонентным акселерометром “ГлобалТест АР2043-50” с высокочастотным фильтром сигнала. В качестве УМСД в системе используется модуль USB3000 производства НПГ "R-Technology". Универсальный модуль сбора преобразований в секунду.

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

В главе проведен сравнительный экспериментальный анализ значений амплитудо-частотных характеристик, рассчитанных по линейной (9) и нелинейной (10) моделям.

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

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

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

Результаты экспериментальных исследований подтвердили с достаточной точностью данные, рассчитанные по модели (6) (рис.7).

Рис.6. Пример численного решения системы нелинейных уравнений описывающих вибрационные перемещения кромки инструмента.

А,мкм.

экспериментальные Аэксп результаты влияния а) глубины резания и б) подачи на амплитуду колебаний.

В этой главе было произведено сравнение экспериментальных и расчетных данных стойкости инструмента из твердого сплава ВК10хом при обработке ХН73МБТЮ на режиме резания b =1 мм, V = 28 м/мин, s0 = 0, мм/об при изменении жесткости технологической системы, которые представлены на рис. 8.

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

Из рис. 8 можно сделать вывод о хорошем (15%) совпадении адекватности предложенных методик и формул.

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

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

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

На базе этих методик определены нормативы на эффективные режимы токарной обработки деталей из группы жаропрочных сплавов (ХН73МБТЮ, ХН62ВМЮТ, ХН60ВТ, ХН77ТЮ) учитывающие динамику процесса резания.

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

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

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

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

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

6. На основании полученных научно-технических результатов был разработан и внедрен на ФГУП ММПП «Салют» автоматизированный стенд, измеряющий силы и вибрации в процессе резания для ускоренного выбора технологических условий токарной обработки деталей.

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

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Горелов В.А., Яковлев М.Г. Выбор оптимальных режимов резания акустико-эмиссионным методом. // Актуальные проблемы Российской космонавтики. Материалы ХХХ1(31) академических чтений по космонавтике. Москва: МГТУ им. Н.Э. Баумана, 2007. с.473-474.

2. Горелов В.А., Черкасова Н.Ю., Яковлев М.Г. Применение автоматизированных систем диагностики на ФГУП “ММПП Салют” для ускоренного выбора оптимальных условий обработки резанием деталей газотурбинных двигателей. // Крылья родины. Национальный авиационный журнал. 2007. №11, с.27-29.

3. Яковлев М.Г., Черкасова Н.Ю. Автоматизированная система диагностики для ускоренного выбора оптимальных условий обработки резанием деталей из жаропрочных сплавов. // Новые материалы и технологии НМТ-2008.

Материалы Всероссийской научно-технической конференции т.2. Москва, 2008. с.58-59.

4. Черкасова Н.Ю., Яковлев М.Г. Моделирование вибраций режущей кромки инструмента при обработке никелевых сплавов в условиях автоколебаний. // Сборник трудов всероссийской конференции молодых ученых и специалистов “Будущее машиностроения России”. Москва: МГТУ им. Н.Э. Баумана, 2008. с. 56-57.

5. Григорьев С.Н., Яковлев М.Г. Повышение эффективности обработки жаропрочных материалов на токарных станках с ЧПУ. // Инструмент технология оборудование. Информационно-аналитический журнал. Москва.

2009. №6, с. 42-43.

6. Григорьев С.Н., Яковлев М.Г. Исследование динамической системы токарного станка при обработке деталей авиационных двигателей. // Справочник. Инженерный журнал. Москва. 2009. №8, с. 62-64.

7. Яковлев М.Г. Исследование динамики процесса резания при обработке жаропрочных материалов. // Наука и образование: электронное научнотехническое издание. Москва. 2009. №8.