Упругое последействие в деталях при высокоскоростной обработке резанием Текст научной статьи по специальности «Физика»

ОБРАБОТКА МАТЕРИАЛОВ РЕЗАНИЕМ

УДК 539.219.2:621.91.01

Упругое последействие в деталях

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

Д. В. Васильков, Т. Б. Кочина

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

Ключевые слова: упругое последействие, остаточные напряжения, резание, высокоскоростная обработка.

Деформирующая способность

остаточных напряжений

как мера упругого последействия

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

1) остаточные деформации или коробление маложестких деталей;

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

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

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

глубине залегания остаточных напряжений Нм и учитывающей значения остаточных напряжений на глубине 0 < Н < Им , целесообразно использовать деформирующую способность остаточных напряжений [2]

Нм

1 (Нм) = / sостdh,

(1)

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

Деформирующая способность остаточных напряжений 1(Нм) может быть физически интерпретирована как момент пары сил, приведенной к единице площади в рассматриваемой точке Ь на поверхности детали ЛБСБ (рис. 1). Данный подход предусматривает рассмотрение действия деформирующей способности в виде распределенной моментной нагрузки. Тогда момент пары сил от действия деформирующей способности остаточных напряжений, приведенной к поверхности изделия хОу в системе координат Охуг, на участке а X Ь определяется как совокупность моментов тх и ту в направлениях осей Ох и Оу (рис. 2).

Моменты тх и ту определяются в соответствии с выражениями:

тх = Тх (Нм)аЬ;

ту = 1у (Нм)аЬ.

(2) (3)

ОБРАБОТКА МАТЕРИАЛОВ РЕЗАНИЕМ

МЕТШЕЕ РАБОТКА

с

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

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

Дифференцирование функции 1(Н) позволяет построить эпюру остаточных напряжений по глубине:

йТ (К)

ан

(4)

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

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

в материалах, склонных к таким изменениям [3, 4].

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

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

В качестве обобщенной количественной меры степени дефектности н-слоя может рассматриваться значение деформирующей способности К-слоя Т(ам). Для идеального образца без остаточных напряжений I (ам) = 0 при отсутствии дефектного слоя или равновесии эпюры напряжений. При наличии дефектного слоя, как правило, I(ам) Ф 0, причем напряженное состояние от действия остаточных напряжений неоднородное и неравновесное, поэтому в зависимости от преобладающего знака остаточных напряжений возможно выполнение неравенств I(ам) < 0 или I(ам) > 0. Деформирующая способность может быть принята в качестве параметрической характеристики НДС. С учетом выражения (1) она представляет собой функционал вида

з = I (Км).

(5)

Отметим свойства, которыми обладает принятый критериальный параметр 3 [5]:

1) величина 3 функционально определяет упругое последействие детали;

№ 2 (80)/2014

3"

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

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

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

3 = 3 (-Х^ Х2, Хз, Х4), (6)

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

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

Представление о действии деформирующей способности остаточных напряжений в виде приведенного к поверхности момента пары сил дает возможность рассмотреть его влияние на формирование остаточных деформаций тонкостенных деталей. Для этого рассмотрим задачу об остаточных деформациях тонкостенной широкополой балки постоянного поперечного сечения. Балка имеет ширину Ь, толщину а и длину l. Представим ее положение в системе координат Оxyz так, как это показано на рис. 3. Балка имеет две плоские горизонтальные поверхности — верхнюю (в) и нижнюю (н). Данные поверхности механически обработаны по различным технологиям. Необходимо определить значение и направление перемещения точки Ь конечного поперечного сечения балки.

По аналогии с (2)-(3) определим моментные нагрузки, действующие на указанных поверхностях. Для этого на каждой из них выделим элементарные участки размерами dx X dy.

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

Деформирующие способности остаточных напряжений, параметризованные в соответствии с (5) со стороны верхней 3^ и нижней 3^ поверхностей будем рассматривать как приведенные моментные нагрузки. Тогда производные крутящих моментов dmxв и dmxн на элементарных участках dx X dy для верхней и нижней поверхностей будут:

dmxв = -3xBdxdy; (7)

dmxн = 3^ dxdy. (8)

Знак минус в выражении (7) определяется принятым правилом знаков.

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

I Ь

^в = -¡¡3xвdxdy; (9)

00

I Ь

mxя = (10)

00

x iz

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

Если предположить, что деформирующие способности остаточных напряжений на верхней 3хв и нижней 3хн поверхностях постоянны, то выражения (9), (10) могут быть существенно упрощены:

тхв 3хвЬ1;

3хнЬ1.

т.

(11) (12)

Тогда полный крутящий момент тх с учетом (11, 12) может быть определен следующим образом

^^^х ттхв I ^^^х'н.

= Ы(3хн -3хв). (13)

Вертикальное смещение точки Ь конечного поперечного сечения балки (рис. 3) будет определяться зависимостью

УЬ =

тх12 2Шх '

(14)

где Е - модуль упругости второго рода; Jx — момент инерции поперечного сечения балки относительно оси Ох.

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

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

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

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

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

1. Функция /1 характеризует напряжения сжатия, уменьшающиеся от некоторого значения А1 = ап1 на поверхности монотонно асимптотически до точки выбега в ноль Е1 в глубину поверхностного слоя Нм1. Поскольку функция /1 при выбеге в ноль асимптотически приближается к оси ОН, то таким приближением можно считать значение Е1 = АН = = 2... 5 МПа. С достаточной для практических целей точностью можно ограничиться значением АН = 5 МПа. При дальнейшем изложении под точкой выбега будем понимать именно данное значение.

2. Функция /2 характеризует напряжения растяжения, уменьшающиеся от некоторого значения А2 = ап2 на поверхности монотонно асимптотически до точки выбега в ноль Е2 в глубину поверхностного слоя Нм2.

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

а2/3

ан2

= 0.

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

в

'Н

Рис. 4. Характерные эпюры остаточных напряжений в поверхностном слое детали при воздействии отдельных физических полей

МЕШПООБМБОТК|»

dh2

= 0.

Функции /1 - /4 могут быть аппроксимированы уравнением экспоненциальной регрессии вида

где а^, Ь^, е^ - коэффициенты регрессии.

В выражении (15) коэффициенты а¿, Ь¿, ег принимают значения:

а1 = ап

ний. При этом возникают сжимающие остаточные напряжения в направлении действия силы (функция /1, рис. 4, напряжения от действия физических полей, обозначенных соответственно Spzz, Spxx) и растягивающие — в перпендикулярном направлении (функция /2, рис. 4, напряжения от действия физических полей, / = а1 ехр -(ЬН) 1 ; г = 1, ..., 4, (15) обозначенных Spzx, Spxz). Интенсивность данных воздействий можно обобщенно свести к соотношению составляющих силы резания Рг : Ру : Рг = 1,0 : 0,6 : 0,4 [7].

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

{2

, г = 1, ...,4;

при Ь = 1,2; при Ь = 3, 4;

ч =

1 1п

^ уДи

1 1п

Л2. Ш1 чДа

при I = 1, 2;

при I = 3, 4.

(16)

(17)

(18)

Выражение (18) получено из (15) подстановкой значений (16), (17).

Действие физических полей будем рассматривать независимо на основе установленной закономерности аддитивности последействия в объемных частях и поверхностных слоях пар трибосопряжения при резании [6].

Нормальная составляющая силы резания Ру формирует в зоне резания объемное действие (рис. 5), уменьшая его объем и вызывая в поверхностном слое обработанной детали сжимающие остаточные напряжения а^ аг, имеющие одностороннюю эпюру, которая асимптотически уменьшается в глубину от некоторого значения на поверхности (функция /1, рис. 4, напряжения от действия физических полей, обозначенные соответственно аруХ, аpуZ). Интенсивность данных воздействий по отношению к прямому воздействию по нормали к обработанной поверхности детали в формообразующей точке А обратно пропорциональна коэффициенту усадки стружки

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

ф-с

зованием зависимости [3]

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

тфс (Н) = 1

Е

Ц

1 -

Уе ум

(Н)

(19)

где Е — модуль упругости 2-го рода материала; ц — коэффициент Пуассона; уе(Н) — функция распределения удельного объема материала по глубине поверхностного слоя, определяемая удельными объемами структурных составляющих: цементит, мартенсит, перлит, тростит, феррит, аустенит; ум — удельный объем основного материала.

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

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

2) поле фазово-структурных изменений ослабляется в глубину поверхностного слоя.

3) скорость деформации в зоне резания при высокоскоростной обработке приближается к шлифованию.

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

5) поскольку температура в зоне резания превышает 600 °С, эту точку можно принять в качестве критериальной при сопоставлении остаточных напряжений с более высокими температурами на основе поправочного коэффициента

к =' ат

ф-с

ав )раб ЧавУб00

(20)

6) механической обработке подвергается, как правило, исходный нетермообработан-ный материал;

7) физико-механические характеристики материала зависят в основном от температуры

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

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

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

Для принятой схематизации деформирующей способности технологических остаточных напряжений (см. рис. 1 и 2) можно показать схему действия сил резания и деформирующей способности, приведенных к формообразующей точке А заготовки (рис. 5). На рис. 5, видно, что на упругое последействие оказывают определяющее влияние деформирующие способности остаточных напряжений 3х и 3г . При этом действием 3у можно пренебречь, так как влияние на упругое последействие оказывают изменения в поверхностном слое в глубину изделия (по нормали к формообразующей точке А на рис. 5). С учетом выражения (1) деформирующие способности остаточных напряжений 3х и 3г определяются зависимостями:

Кг

3х

аг(Н)йН;

Зг ох (Н)йН.

(21)

(22)

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

На основе выражения (15) с учетом (16-18) можно реализовать расчетно-эксперименталь-

н

МЕШПООБМБОТК|»

ный подход к определению технологических остаточных напряжении:

¿=3 )=5

щ х О) = Е ° хь(к) + Е ° х)(ь) =

¿=1 )=4

¿=3 ) =5 2

= Е °ИХ1 еХР [- (^¿^ + Е Стпх) ехР - (Ъх^Г

I=1 )=4

(23)

¿=3 ) =5

Щ2 (й) = Е С2г (Ь) + Е °г) (к) =

¿=1 )=4

¿=3 )=5 - 2

= Е CTП2¿ ехР [- (Ъz¿h)] + Е СТпг) ехР - ())

¿=1 )=4

(24)

Элементами выражении (23, 24) являются остаточные напряжения, представленные в таблице.

Из выражении (23, 24) могут быть выделены компоненты, в результате интегрирования которых по глубине получаем компоненты деформирующей способности остаточных напряжений:

3 ■ =

О

I г

1п

пг1 ъг1 при I = 1,2, 3; от 1

ДО

О Н2 ■

^пгу ^мг1

(25)

1п

До

при I = 4, 5;

3 ■ = 21

а

^мя

ПХ1 "МХ1 1п|

А а а к2 ■

при I = 1, 2, 3;

(26)

1п

п ап

Аа

при ! = 4, 5.

Полученные в (25, 26) компоненты деформирующей способности З^ и Зг£ являются функционалами и позволяют интегрально оценивать упругое последеиствие остаточных напряжений. Однако для управления техно-логиеи обработки целесообразно знать характер изменения компонентов деформирующей способности остаточных напряжений Зх и Зг

2

2

Компоненты остаточных напряжений

Компонент остаточных напряжений Знак Элемент физического поля Остаточные напряжения на поверхности ап Глубина залегания остаточных напряжений hм

°х1 - Орхх 0,4°П2 0,4^2

°х2 + аРгх hмz

°х3 - аРух 1 ^ ЩП2 0,6hм2

°х4 +

°х5 ± °ф-с °Пф-с hМф-с

Щ2

- аРгг °пг = г|0т ^2 =Аз ^

+ аРх2 0,4цоп2 0,4^2

°г3 - ®Руг 1 ^ ЩП2 0,6^2

°г4 +

°г5 ± °ф-с °Пф-с hМф-с

(21, 22) с учетом выражений (23, 26) в соответствии с выражениямии:

I=5

Зх = ^ Зх1; ¿=1

1=5

Зг = ^ Зг1 •

I=1

(27)

(28)

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

Пример реализации модельных решений

В качестве примера реализации полученных модельных решений были приняты следующие условия: обработкаточением торцевойповерх-ности диск а диаметром 300 мм и толщиной 2 мм из сплава ХН73МБТЮ. Для обработки использовался сборный инструмент с твердо-

сплавной пластиной марки 004235 (8ап^1к Согошап^ и пластиной из минералокерами-ки М1018 (Тип§а1оу). Пластина ромбическая с задним углом ССМТ 11 Т3 08^Е, державка 8Б.ЛСК 2020К11. Основные углы: угол при вершине е = 55°; главный угол в плане ф = 92°; вспомогательный угол в плане ф1 = 33°, задний угол а = 8°, передний угол у = 8°. Режимы резания: скорость резания м = 25 у 260 м/мин, глубина ре занбя О = 1 мм,подача = 0,19 мм/об.

Выполнен расчет компонентов остаточных напряжение соответственно (еадиальных <тх и ркруакных гл,,, а текнкг дефодмируютцей снособ-ности остаточных напряжений Зх, Зг. Упругое последействие рассматриваемой детали определяет коробление полотна диска, которое не должно превышать 0,03 мм. Если принять, что диск расположен в плоскости хОу (рис. 5), то определяющими с точки зрения формирования упругого последействия являются радиальные остаточные напряжения сх и деформирующая способность остаточных напряжений Зг. Окружные остаточные напряжения! ог оказывают второстепенное действие, ио они обязательно должны быть сжимающи-

50

100

150

200 250 /г, мкм

300

200

®пхх св 100

®пгх ^п ух к § о" 0

— — Ч -100 -200

— -300

Рис. 6. Компоненты эпюры остаточных напряжений сх

—♦- ®пгг

-Ш- ®пхг

-*-

-*- — Ч

—

—

50

100

150

200 250 Л, мкм

1*110. 7. Компоненты эпюры остаточных напряжений а,

15 10

§ 5 й 5

К 0

к

м -5 -10 -15

к—| Ни

- Зг5

50

100

150

200 250 Л, мкм

Рис. 8. Компоненты эпюры остаточных напряжений З,

15 10

й 5 К 0

к

п -5 -10 -15

е—*—* 6—*—» (-x-к

с ж ж ж x

¡=Х=!

50

100

150 200 250

Л, мкм

Рис. 9. Компоненты эпюры остаточных напряжений З х

О

О

О

О

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

Экспериментальные исследования проводились на токарном обрабатывающем центре фирмы OKUMA (Япония) модели GENOS L300E-M с ЧПУ OSP-P200L. В качестве измерительной аппаратуры был использован Измерительно-вычислительный комплекс «СИТОН» [8], который позволяет неразрушающим методом измерять остаточные напряжения в поверхиостном слое деталей из токопроводящих металлов и сплавов. В комплексе реализован резистивный электроконтактный метод неразрушающего контроля остаточных напряжений, который основан на корреляции между интемральны-ми электрическими и механическимихарак-теристиками металлов и сплавов — удельным электрическим сопротивлением в слоа глубиной h и деформирующей способностьюоста-точных напряжений.

Сопоставление расчетных (см. рис. 6) и экспериментальных данных по определению остаточных напряжений sx при скорости резания v = 25 м/мин представлено на рис. 10, а при скорости резания v = 260 м/мин — на рис. 11. Здесь средние значения по экспериментальным данным получены на основе выборки из восьми реализаций.

Сопоставление расчетных и экспериментальных данных показало, что максимальное отклонение между ними не превышает 25 %. Причем напряжения на поверхности, точка перехода нуля и подповерхностный максимум да-

h, мкм

Рис. 10. Расчетные и экспериментальные данные по определению остаточных напряжений ах прискорости резания V = 25 м/мин:

--расчетная кривая; ♦ — расчетныеточки; ■ — средние

значения по экспериментальным данным

к, мкм

Рис. 11. Расчетные и экспериментальные данные по определению остаточных напряжений ах при скорости резания V = 260 м/мин:

--расчетнаякривая; ♦ — расчетныеточки; ■ — средние

значения по экспериментальным данным

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

Литература

1. Кочина Т. Б. Теплофизические характеристики процесса резания при высокоскоростной обработке сплавов на никелевой основе инструментами из минерало-керамики // Науч.-техн. вед. СПбГПУ: сер. «Наука и образование». 2009. № 4-2 (89). С. 63-68.

2. Комплексное исследование состояния поверхностного слоя ответственных деталей машинных агрегатов /Д. В. Васильков [и др.] // Тяжелое машиностроение. 1998. № 3. С. 31-34.

3. Овсеенко А. Н., Клауч Д. Н., Кущева М. Е. Технологическое обеспечение качества изделий машиностроения / Моск. гос. технол. ун-т «Станкин». М.: Янус-К, 2006. 118 с.

4. Максаров В. В., Ольт Ю., Халимоненко А. Д. Управление качеством процесса точения инструментом из режущей керамики. СПб.: Изд-во ГОУ ВПО «СЗТУ», 2010. С. 160.

5. Васильков Д. В., Иванов С. Ю., Хитрик В. Э. Физические основы метода АФЧХ-тестирования технологических остаточных напряжений / / Диагностика, эксплуатация, ремонт, восстановление, модернизация оборудования. Современные технологии: сб. докл. 3-й Междунар. науч.-практ. конф. СПб.: Изд-во ПИМаш, 2005. С. 80-90.

6. Закономерность аддитивности сорбционного последействия в объемных частях и поверхностных слоях пар трения: науч. открытие (диплом № 466) / Д. В. Васильков, Г. С. Ивасышин. М: РАЕН. МААНОиИ, 2014.

7. Васильков Д. В., Кочина Т. Б. Силы резания при высокоскоростной обработке сплавов на никелевой основе инструментами из минералокерамики // Металлообработка. 2011. № 4 (64). С. 10-13.

8. Васильков Д. В. Вейц В. Л., Шевченко В. С. Динамика технологической системы механической обработки. СПб.: Инструмент, 1997. 230 с.