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

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

УДАЛОВ Сергей Владимирович

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

Специальность 05.02.04 - Трение и износ в машинах

АВТОРЕФЕРАТ

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

□ □347283 1

Тверь 2009

003472831

Работа выполнена в Тверском государственном техническом университете.

Научный руководитель - доктор технических наук, профессор БОЛОТОВ А.Н.

Официальные оппоненты:

- доктор технических наук, профессор ЛУКЬЯНЧИКОВ А.Н.

- доктор технических наук, профессор ПОГОНЫШЕВ В.А.

Ведущая организация - ЗАО НО «Тверской институт вагоностроения»

Защита диссертации состоится (МО'аЛ 2009 г., в часов на

заседании диссертационного совета Д262.262.02 Тверского государственного технического университета по адресу: 170026, г. Тверь, наб. Афанасия Никитина, 22.

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

/

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

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

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

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

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

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

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

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

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

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

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

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

, Провести сопоставление математической модели с экспериментом.

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

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

Научная новизна.

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

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

критерий (2К, т.е. при относительном давлении (2К<~ 0,2 расчеты по

традиционной методике и предложенной модели сходятся, а при ()к> 0,2

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

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

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

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

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

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

Апробация работы. Основные положения работы были доложены и обсуждены на:

V Международном семинаре «Физико-математическое моделирование систем», Воронеж, 2008;

научно-техническом семинаре «Контактные задачи в трибологии», Тверь, кафедра «Физики», ТГСХА, 2008;

научно-технических семинарах «Механика и физика фрикционного контакта», кафедра «Физики», Тверской государственный технический универсигет, 2006 - 2008 г.;

научно-технических конференциях профессорско-преподавательского состава и заседаниях кафедры «Сопротивления материалов», Тверской государственный технический университет, 2008 г.; объединенном семинаре «Контакт шероховатых, волнистых поверхностей» кафедры «Физики» и «Сопротивления материалов», Тверской государственный технический университет, 2008; научно-техническом семинаре «Методы расчета посадок с натягом», Тверской институт вагоностроения, 2008.

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

Структура работы. Диссертациогшая работа состоит из введения, пяти глав, выводов но работе, списка основных обозначений и списка литературы и приложений. Полный объем диссертации вместе с иллюстрациями составляет 135 страниц, включая 32 рисунка и 6 таблиц. Список литературы включает 118 наименований.

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

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

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

Основное внимание уделено работам в области изучения закономерностей формирования фрикциогаюго контакта и теории контактного взаимодействия шероховатых поверхностей выполненными исследователями научных школ Института проблем механики РАН, Института машиноведения РАН, Гомельского института металлополимерных систем, Брянского государствешюго технического университета, Тверского государственного технического университета, Волгоградского государственного технического университета и других организациях.

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

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

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

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

S[q(x)] = C[q(x)r , (1)

где q{x) - функция распределения давления, С и т- коэффициенты, зависящие от распределения материала в поверхностном слое и его механических свойств.

Примените функции в виде (1) ограничено случаем, когда

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

^ 1

Яа

глубина сглаживания шероховатости, /

(2)

V - параметры опорной

где

кривой, НВ - твердость материала мшфовыступов, ца - номинальное

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

а=Кр.

I

у1тНВ

V «

1+

Ча

\ Я*

(3)

Результаты расчета по формулам (2) и (3) сопоставлены на рис. 1.

Рис. 1. Графики а ' | при

КЛНВ

зависимости различных

значениях параметра V с учетом и без учета эффекта взаимного влияния

микронеровностей: Где (!) - V - 2,

(2) - V - 3 , (4) - V = 2 (с учетом эффекта взаимного влияния), (3) -V — 3 (с учетом эффекта взаимного влияния).

относительное давление — НВ

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

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

что при Ял- <02. характерном для процесса трения и изнашивания, НВ '

приведенные кривые практически совпадают, а при >02 наблюдаются

НВ '

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

(4)

где Л~ коэффициент, зависящий от физико-механических свойств материалов и микрорельефа контактирующих поверхностей, для упрощенных расчетов можно принять Я = 0,5.

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

(3):

На рис. 2 представлены результаты расчетов характеризующих зависимость сближения сферы радиусом Я =10 мм с шероховатым полупространством (/?Р=10 мкм, = 0.5, V ~ 2) и его составляющих от действующей на нее силы. Здесь же приведены аналогичные зависимости, полученные без учета эффект взаимного влияния микронеровностей, а также зависимость сближения сферы с гладким упругим полупространством, рассчитанные по формулам Герца, и зависимость сближения сферы с шероховатым слоем на недеформируемом основании. В последнем случае эффект взаимного влияния микронеровностей не учитывался и сближение рассчитывалось по формуле (6)

1

( (у + 1)ЛГ

/ 2

1+ (5)

V 1 НВ) j

а=Яр.

ч2 яЖр1тНВ

(6)

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

значительным погрешностям.

Рис. 2. Графики зависимости сближения сферы с шероховатым упругим полупространством (общее, полупространства, шероховатого слоя) от действующей на нее силы с учетом и без учета эффекта взаимного влияния микронеровностей ( 1,2,3 - сближения (общее, полупространства, шероховатого слоя) без учета эффекта взаимного влияния, 4,5,6 - сближения (общее, полупространства, шероховатого слоя) с учетом эффекта взаимного влияния, 7 - сближение по Герцу, 8 - сближение по формуле

(5))-

Рис. 3. Распределения контактного Рис. 4. Графики зависимости давления, рассчитанные с учетом (кривая относительной площади контакта от 1) и без учета (кривая 2) эффекта относительного давления с учетом (кривая взаимного влияния. 1) и без учета (кривая 2) эффега

взаимного влияния.

На рис. 4 представлены зависимость относительной площади контакта от относительного давления. Относительная фактическая площадь контакта была рассчитана по формуле:

л=— ] ч(х)2жхсЬс (7)

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

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

индентора, представлена в виде

гр{х,у)=Ьх У (7), гДе расстояние 2 Ь а

от оси Ох и у по оси Оу от центра, а, Ь и с полуоси эллиптического индешора.

Полагая, что микровыступы шероховатого слоя деформируются пластически и функция Я[д(х, у)]> учитывая (3) примет вид:

S[q(x,y)]=Rp

\ tmHB

\ «

1+

f , Ч(х,у)

у tjm, ^

(8)

В предположении упругих деформаций

Рис. 5 Схема контакта упругого шероховатого эллипсоидального сегмента с упругим полупространством

волн перемещение

полупространства в направлении оси Oz может быть найдено по контактным

давлениям из соотношения:

тЯ

где

У

Для нахождения распределения давления по площадке контакта и внедрения эллипсоида используем следующее интегральное уравнение:

лАА,' ^ IV Ь'х'+а'у2 (д(х,у) 4-\с!(р \ч(х',у )с!г1 =И\--—2--К I-

71 й о

26 а

\ ¡т1Ш

1 +

ф,у)

1-НВ

(9)

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

V р

где V №= 2/„

Ъ.

-1

(10)

Полагая, что параметры и '-¡с незначительно зависят от направления профилирования, то выражение для вычисления среднего радиуса кривизны волнистой поверхности имеет вид:

(11)

Я = 0,22*1

Зр ПР^к поп.

ТГ.

где ПГ !ЮП - шаги волн в продольном и поперечном направлениях

относительно следов обработки.

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

(12)

о

где # - нагрузка на единичную волну, пт№ - число волн поверхности на

высоте средней линии, (рш - функция распределения волн по высоте.

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

индентора с полупространством представим нагрузку, воспринимаемую стыком в виде:

2лК№р -" ~'ЦгрГг1'

где ^ - еближехше волнистых поверхностей под действием

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

формуле: "г г,-^ Ас= I жр,, апс

о

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

На рис. 6 представлены зависимости сближения волнистых поверхностей от параметров

микрогеометрии и механических свойств материалов контактирующих поверхностей.

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

Контурная площадь контакта волнистых поверхностей была рассчитана по формуле:

(14)

где дл - единичная расчетная

площадь контакта.

Полагая, что глубине уу

число волн на по площади контакта

выражается

И>г=Ип

уг\

'1¥р зависимость

формулой и учитывая,

16 УУр^икм

Рис. 6 Зависимости сближения поверхностей от параметров микрогеометрии

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

^ _ Лрг 2 = ^ АЛ получим следующее выражение для контурной

<=1 м

площади контакта:

где Л'™», у1г

2

АСг - *Р: Ч, ("„"»М* * А)"

П5)

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

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

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

У=АаИр+А^Ур-}У0

12 14 мкм

0.2

0.15

0.1

0.05

10

12

14 16 №р, мкм

18

20

Е,1011Па

Рис. 7 Влияние параметров микрогеометрии и механических свойств материалов

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

_A.-AV.-bV.

(16)

где дуа - уменьшение объема материала вследствие упругой деформации шероховатого слоя; д у - уменьшение объема материала вследствие упругой деформации волн.

В качестве характеристики объема зазора удобно использовать среднюю высоту зазора , _ V .

На рис. 8 представлены зависимости, характеризующие влияние

геометрических характеристик на величину высоты зазора, рассчитанную но формуле (16). Данные полученные в результате расчета показывают большое влияние параметров шероховатости и волнистости на величину объема зазора.

В четвертой главе приведены результаты экспериментальных

исследований влияния шероховатости на формирование контурной площади контакта и сближения поверхностей постоянной кривизны. Исследования производились при получении контактных деформаций и площади касания гладкого твердого сферического индентера с плоскими поверхностями различной степени шероховатости. Индентер был изготовлен из стали У7А с последующей закалкой до НЫС 60. Радиус индентсра составлял 100 мм, что соответствует типичным значениям ириведегаюго радиуса кривизны волн.

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

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

Рис. 8 Влияние параметров на расчет высоты зазора с помощью предлагаемой модели.

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

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

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

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

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

ла Ы

К- коэффициент, зависящий от способа получения соединений с натягом; при сборке под прессом % - 2У»; при сборке тепловым способом (нагревом втулки или охлаждением вала) # =1; £, и Ег - модули упругости вала и втулки; ц и - коэффициенты пуансона вала и втулки; ^ и ^ ■ диаметры вала и втулки до соединения натягом; с1 - номинальный диаметр посадочной поверхности, с - параметр, рассчитанный согласно решения задачи Ламме.

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

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

Был проведен расчет натяга для запрессовки вала на ротор для прядильно-мотального автомата и запрессовки колеса на ось для колесной пары РУ1Ш-950 (Рис 9,10).

Рис 9. Виды роторов, Рис 10. Колесная пара РУIШ-950.

вал запрессованный на ротор.

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

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

Для анализа силы трения представим коэффициент трения в виде суммы адгезионной и деформационной составляющих:

/ = /,+/л

Зависимость усилия выпресовки от вида материала

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

коэффициент трения и его адгезионная и деформационная составляющие.

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

Экспериментальные данные

приведенные на рисунке также подтверждают адекватность

предлагаемой модели.

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

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

экспериментом и классической теорией контакта реальных поверхностей.

ОД МО2-^2)

О

/Дз+ъэ-

дЯ *д3У2

У

(20)

НВ1 *[р' - А1

Параметр А характеризующий топографию детали можно определить следующим образом:

Л ЬЩ)\*НВУ' (21)

Г У

г дс/у

где ^ - контурное давление на контакте, г - радиус микронеровности, 5|д(;с)] -

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

ОСНОВНЫЕ НАУЧНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ.

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

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

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

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

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

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

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

1. Демкин Н.Б., Удалов C.B. Влияние волнистости цилиндрической формы на контурную площадь стыков деталей машин //Вестник Тверского Государственного Технического университета. 2005. № 8. С. 47-51

2. Болотов А.Н.,Удалов C.B. Влияние скольжения структурированной магнитной жидкости на трение в подшипниках и уплотнениях //Механика и физика процессов на поверхности и в контакте твердых тел и деталей машин: Межвуз. сб. науч. тр. Тверь: ТГТУ, 2006. С. 52-54

3. Алексеев В.М., Удалов C.B. Оценка эффекта взаимного влияния поверхностных неровностей на контакте сферы с шероховатым упругим полупространством //Вестник Тверского Государственного Технического университета. 2006. № 9, С. 6-11

4. Удалов C.B. Контакт эллипсоида с шероховатым упругим полупространством //Механика и физика процессов на поверхности и в контакте твердых тел и деталей машин: Межвуз. сб. науч. тр. Тверь: ТГТУ, 2007. С. 25-31

5. Демкин Н.Б., Удалов C.B. Моделирование контактного взаимодействия поверхностей с учетом эффекта взаимного влияния поверхностных неровностей, имеющих волнистость цилиндрической формы //Трение и смазка в машинах и механизмах. 2007. № 9. С. 10-14

6. Демкин Н.Б., Удалов C.B., Алексеев В.А., Измайлов В.В., Болотов А.Н. Контакт шероховатых волнистых поверхностей с учетом взаимного влияния неровностей //Трение и износ. 2008. № 3. С. 231-238

Выражаю глубокую признательность д.т.н., профессору Дсмкину Н.Б. и к.т.н., доценту Алексееву В.М. за консультации по теме диссертации и критические замечания.

Подписано в печать 13.05.09 Физ.печ.л. 1,25 Тираж 100 экз. Заказ №50 Типография Тверского государственного технического университета 170026, г. Тверь, наб. Афанасия Никитина, 22

Введение

1. Глава 1. Контактное взаимодействие твердых тел.

1.1 Характеристики поверхностей, определяющих контактное взаимодействие твердых тел.

1.2 Роль контактное взаимодействия в процессах трения и 5 изнашивания.

1.3 Контакт тел сферической формы.

1.4 Расчет контактных деформаций.

1.5 Последние достижения в области контактного взаимодействия.

1.5.1 Приближенное решение двумерной контактной задачи теории 13 упругости о вдавливании параболического штампа в упругое слоистое полупространство.

1.5.2 Исследование влияния зазора на статическое трение в полимер- 14 полимерных микроподшипниках.

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

1.5.4 Напряженное состояние упругих тел с включением при 21 контактном взаимодействии.

1.5.5 Механика фрикционного металлополимерного контакта. 22 2.

Глава 2. Контакт единичной волны с шероховатой поверхностью.

2.1 Контакт единичного сферического индентора с 26 полупространством.

2.2 Взаимное влияние микронеровностей при деформации.

2.3 Использование эллипсоидов для моделирования микронеровностей.

2.4 Расчет контурной площади контакта.

2.5 Влияние параметров качества контактирующих поверхностей 44 на их деформацию.

3 Глава 3. Контактное взаимодействие волнистых поверхностей.

3.1 Геометрические характеристики волнистых поверхностей, определяющих их контактное взаимодействие.

3.2 Сближение поверхностей под действием нормального 47 нагружения.

3.2.1 Контактные деформации волнистых поверхностей.

3.2.2 Контурная площадь сжатых поверхностей. 53 3.3 Объем зазора между контактирующими поверхностями.

Глава 4. Экспериментальная проверка полученных результатов.

4.1 Исследование влияния шероховатости на свойства контакта 62 криволинейных поверхностей.

4.1.1 Определение площади контакта сферического индентора с 64 шероховатой поверхностей.

4.1.2 Определение сближения при контактировании сферы с 66 шероховатой поверхностей.

4.2 Определение контактных деформаций плоских стыков.

4.3 Экспериментальное определение объема зазора между 85 контактирующими поверхностями.

4.3.1 Методика определение объема зазора и ее обоснование.

4.3.2 Определение объема зазора между контактирующими 87 поверхностями.

4.4 Экспериментальное определение усилия выпресовки для 101 определения силы трения покоя.

5. Глава 5. Применение предложенной модели контактного 105 взаимодействия. Соединения с натягом. Трение в резьбовых соединениях и в соединениях с гарантированным натягом.

5.1 Характеристика соединений с натягом.

5.2 Основные методы расчета соединений с натягом.

5.3 Расчет сближения и фактической площади в соединениях с 106 натягом.

5.4 Пример расчета прочности соединений собранных с натягом.

5.4.1 Методика расчета.

5.4.2 Пример расчета на прочность соединения для запрессовки 112 колеса на ось для колесной пары РУ1Ш-950.

5.5 Трение в резьбовых соединениях и в соединениях с 117 гарантированным натягом.

5.5.1 Определение сил трения покоя при пластических деформациях 117 в зонах касания.

5.5.2 Нагрузки в резьбовых соединениях.

5.6 Влияние волнистости и эффекта взаимного влияния на 120 износостойкость.

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

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

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

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

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

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

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

Выводы по работе.

На основании проведённых исследований можно сделать следующие выводы:

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

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

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

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

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

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

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

- V Международном семинаре «Физико-математическое моделирование систем», Воронеж, 2008;

- научно-техническом семинаре «Контактные задачи в трибологии», Тверь, кафедра «Физики», ТГСХА, 2008;

- научно-технических семинарах «Механика и физика фрикционного контакта», кафедра «Физики», Тверской государственный технический университет, 2006 - 2008 г.;

- научно-технических конференциях профессорско-преподавательского состава и заседаниях кафедры «Сопротивления материалов», Тверской государственный технический университет, 2008 г.;

- объединенном семинаре «Контакт шероховатых, волнистых поверхностей» кафедры «Физики» и «Сопротивления материалов», Тверской государственный технический университет, 2008;

- научно-техническом семинаре «Методы расчета посадок с натягом», Тверской институт вагоностроения, 2008.

Разработанные методы и основные результаты опубликованы в шести статьях:

1. Демкин Н.Б., Удалов С.В. Влияние волнистости цилиндрической формы на контурную площадь стыков деталей машин //Вестник Тверского осударственного Технического университета. 2005. № 8. С. 47-51

Л°ТОВ А-Н"Удалов Св- Влияние скольжения структурированной

Н°И ЖИДК0СТИ на тРение в подшипниках и уплотнениях //Механика и

Г:ГРГ6СС0В На "ОВерХНОСТИ и в контакте твердых тел и деталей машин: Межвуз. сб. науч. тр. Тверь: ТГТУ, 2006. С. 52-54 пп*^ГСееВ В'М" С'В- УдаЛОВ СВ- °ценка эффекта взаимного влияния Г"™ Нер°В™СТеЙ на ^нтакте сферы с шероховатым упругим

ZseZ7aH7nT Тверского Государственного Технического университета. 2006. № 9, С. 6-11 пиДШТОВ °'В' /лКоНХакт эллипсоида с шероховатым упругим полупространством //Механика и физика процессов на поверхности и в

2007асе25-31™ Т6Л И Д6ТаЛеЙ МаШИН: МбЖВу3- Сб- наУч- тр- Тверь: ТГТУ,

5. Демкин Н.Б., Удалов С.В. Моделирование контактного взаимодействия поверхностей с учетом эффекта взаимного влияния поверхностных неровностей, имеющих волнистость цилиндрической формы //Трение и смазка в машинах и механизмах. 2007. № 9. С. 10-14

6. Демкин Н.Б., Удалов С.В., Алексеев В.А., Измайлов В.В., Болотов А.Н. кхштакх шероховатых волнистых поверхностей с учетом взаимного влияния неровностей //Трение и износ. 2008. № 3. С. 231-238

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

Заключение.

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

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

1. Авдеев М.В, Воловик E.J1. Ульман И.С. Технология ремонта машин и оборудования. М.: Агропромиздат, 1986. - 247.

2. Аксёнов А.Ф. Трение и изнашивание в углеводородных жидкостях. М.: Машиностроение, 1977. 147.

3. Алюшин Ю.А, Еленев С.А. Применение энергетического метода для расчёта и анализа процессов пластического формоизменения металлов II Исследование процессов пластического деформирования металлов. М.: Наука, 1965. С. 106-133.

4. Атопов В.И., Сердобинцев Ю.П., Славин O.K. Моделирование контактных напряжений. Москва: Машиностроение (1988).

5. Ахматов А.С. Молекулярная физика граничного трения. М.: Физматгиз, 1963. 742 с.

6. Барт В.Е., Васильев Ю.Н. и др. Триботехнические свойства антифрикционных самосмазывающихся пластмасс.— М., ВНИУ ГСССД.— 1982. -62 с.

7. Березин И.С., Жидков Н.П. Методы вычислений, тт.1,2,М: Наука, 1966.-670,750.

8. Близнец М.М., Богданович П.Н., Стухляк П.Д. Особенности изнашивания эпоксидных полимеров.- Трение и износ. 1988, т.9, №3.-с.549-553.

9. Борисов М.В. Павлов И.Ф. Постников В.И. Ускоренные испытания машин на износостойкость как основа повышения качества. М.: Издательство стандартов, 1976. — 351.

10. Болотов А.Н., Горлов И.В. Восстановление изношенных поверхностей методом пластического деформирования // Механика и физика фрикционного контакта. Межвузовский сборник научных трудов. Выпуск 9. Тверь. ТГТУ, 2002. С. 39-43.

11. Болотов А.Н., Горлов И.В. Сравнительные триботехничекие испытания материалов // Механика и физика фрикционного контакта. Межвузовский сборник научных трудов. Выпуск 10. Тверь. ТГТУ, 2003. С 40-44.

12. Браун Э.Д., Евдокимов Ю.А., Чичинадзе А.В. Моделирование трения и изнашивания в машинах. М.: Машиностроение, 1982. 190 с.

13. Буше Н.А. Трение, износ и усталость в машинах. Транспортная техника. М.: Транспорт, 1987. 223.

14. Воловик E.JI. Справочник по восстановлению деталей.— М.: Колос, 1981.351 с.

15. Восстановление деталей футерованием композиционными полимерными материалами на основе эпоксидных смол. Минск, БелНИИНТИ. 1988.-50 с.

16. Горлов И.В. Применение восстановленных деталей и новые технологии восстановления // Актуальные проблемы развития машиностроительного комплекса Тверской области. Материалы научно-практической конференции. Тверь: ТГТУ, 2001. С.79-81.

17. Джонсон К. Механика контактного взаимодействия. Перевод с английского. М.: Мир, 1989. 500.

18. Евдокимов Ю.А., Колесников В.И., Тетерин А.И. Планирование и анализ экспериментов при решении задач трения и износа. М.: Наука, 1980. 228 с.

19. Зубчанинов В.Г. Математическая теория пластичности. Тверь: ТГТУ, 2002. 300с.

20. Избирательный перенос при трении и его экономическая эффективность/ Под ред. Д.Н. Гаркунова и А.А. Полякова. М.: ВИНИТИ, 1972. 252 с.

21. Избирательный перенос при трении/ Под ред. Д.Н. Гаркунова и Ю.С. Симакова. М.: Наука, 1975. 88 с.

22. Лавринович М.Ф., Шустерняк М.Ф. Повышение износостойкости деталей автомобилей. Минск: Беларусь, 1985. 142.

23. Канарчук Ф.Д., Чигринец А.Д., Голяк О.Л., Шоцкий П.М. Восстановление автомобильных деталей. Технология и оборудование. М.: Транспорт, 1995.-303.

24. Карасик И.И. Прирабатываемость материалов для подшипников скольжения. М.: Наука, 1978. 136 с.

25. К вопросу нормирования протяжённости фактического контакта и шага микроканавок поверхностей с частично регулярным микрорельефом (ЧРМР) / B.C. Комбалов и др. // Трение и износ. 1992.Т 13. №1. С.110-115.

26. Композиционные материалы. Т.З. Пер. с англ./ Под ред. Л. Браутмана, Р. Крока.— М.: Машиностроение, 1978.

27. Корн Г, Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука. Главная редакция физико-математической литературы, 1984. 831.

28. Крагельский И.В. Трение и износ. Изд.2-е. М.: Машиностроение, 1968. 480 с.

29. Lockett F.J. Indentation of a rigid-plastic material by a conical indenter.— J. Mech. and Phys. of Solids, 1963, p. 345.

30. Масино M.A., Алексеев B.H., Мотовилин Г.В. Автомобильные материалы. М.: Транспорт, 1979. 169.

31. Маслов Н.Н. Эффективность и качество ремонта автомобилей.— М.: Транспорт, 1981. 304 с.

32. Машины и стенды для испытаний деталей / Под ред. Д.Н. Решетова. -М.: Машиностроение, 1979. -343.

33. Mellor Р.В., Johnson W. Plasticity for Mechanical Engineers, London, 1962.

34. Методы профилактики и ремонта промышленного оборудования: Учеб. Для сред. проф. Образования / Ю.Н. Воронкин; Н.В. Поздняков. М.: Академия, 2002. - 240 е.: ил.

35. Механика пластических сред / Д.Д. Ивлев; Н.Б. Бартошевич-Жегаль. Т. 1: Теория идеальной пластичности; Ред. Н.Б. Бартошевич-Жегель. -М.: Физмат, 2001.-445 с. : ил.

36. Муляр Ю.Н. Антифрикционные смазывающие материалы в машиностроении // Машиностроительное производство. Серия Прогрессивные технологические процессы в машиностроении. Выпуск 3. М.: ВНИИТЭМР, 1991. 64.

37. Нелинейное деформирование твёрдых тел / С.Н. Коробейников; Рос. АН, Сиб. отд-ние, Ин-т гидродинамики им. М.А. Лаврентьева. -Новосибирск: Изд-во СО РАН, 2000. 261 с.

38. О моделировании трения и износа/ Под ред. А.В. Чичинадзе. М.: НИИМАШ, 1970,318 с.

39. Операции соединения и механической обработки композиционных материалов// Dastie S.I.— Handbook of Composites-SI.— 1982. P. 602-632./ВЦП. №1. 58270. Пер. 86/17806.

40. Пат. РФ, МКП В 23 Р 6/00, 6/02. Способ восстановления локально изношенной поверхности детали /А.Н. Болотов, И.В. Горлов № 2228247 С1; Заявл. 20.12.2002; Опубл. 10.05.2004; Бюл. № 13. - 8с.

41. Полимеры в узлах трения машин и приборов / Под ред. А.В. Чичинадзе. М.: Машиностроение, 1980. 208.

42. Полимеры в узлах трения машин и приборов / Справочник. А.В. Чичинадзе, А.Л.Левин и др.- М.: Машиностроение, 1988. - 328 с.

43. Польцер Г. Майснер Ф. Основы трения и изнашивания. Перевод с немецкого. М.: Машиностроение, 1984. -263.

44. Применение композиционных материалов в машиностроении. Тезисы докладов научно-технической конференции.- Минск, Бел НИИНТИ, 1988.-196 с.

45. Прохоров B.C. Трибологические методы испытания масел и присадок. М.: Машиностроение, 1983. 183.

46. Расчётно-экспериментальный метод решения контактных задач / А.Г. Кузьменко и др. // Трение и износ. 1992.Т 13. №2. С. 257-264.

47. Ремонт автомобилей: Теорет. Курс: Учеб. пособие / А.В. Коробейник. Ростов н/Д.: Феникс, 2003. - 283 е.: ил.

48. Ремонт автомобилей и двигателей: Учеб. для сред. проф. Образования / В.И. Карагодин; Н.Н. Митрохин. 2-е изд., стереотип. -М.: Академия, 2002. - 496 е.: ил.

49. Семенов А.П., Савицкий М.Э. Металлофторопластовые подшипники. М.: Машиностроение, 1976. 190 с.

50. Сиренко Г.А., Свидерский В.П. Износостойкость композиционных материалов на основе термостойких полимеров в условиях предельных нагрузок и ограниченной смазки.— Трение и износ.— 1988, т.9, № 5.- с. 841-861.

51. Справочник по композиционным материалам. В 2-х томах/ Под ред. Дж. Любина.— М.: Машиностроение, 1988. Т.1.— 448 с.

52. Справочник по композиционным материалам. В 2-х томах/ Под ред. Дж. Любина.— М.: Машиностроение, 1988. Т.2.— 582 с.

53. Справочник по триботехнике. Смазочные материалы, техника смазки, опоры скольжения и качения тт, 1,2 / Под ред. М. Хебды, А.В. Чичинадзе. М.: Машиностроение, Варшава. ВКЛ. 1990. 820.

54. Shield R.T. On plastic flow of metals under conditions of axial symmetry.— Proc. Roy. Soc., 1955, p. 267.

55. Струк B.A. Композиционные материалы для узлов трения. Опыт разработки и применения. Минск: БелНИИНТИ, 1986.-46 с.

56. Суслов А.Г. Технологическое обеспечение параметров состояния поверхностного слоя деталей. М., 1987.

57. Рыжов Э.В., Суслов А.Г., Фёдоров В.П. Технологическое обеспечение эксплуатационных свойств деталей машин. М., 1979.

58. Теоретические основы процессов поверхностного пластического деформирования / Под ред. В.И. Беляева. Минск. Наука и техника, 1988.-360.

59. Томлёнов А.Д. Теория пластического деформирования металлов. М.: Металлургия, 1972. 408.

60. Томленов А.Д. Механика процессов обработки металлов давлением. Машгиз, 1963.

61. Трение изнашивание и смазка. Справочник тт,1,2 / Под ред. И.В. Крагельского, В.В. Алисина. М.: Машиностроение, 1978. 812.

62. Упрочнение машиностроительных материалов. Справочник / С.Н. Полевой. В.Д. Евдокимов. М.: Машиностроение. 1994. 471.

63. Федорченко Н.М. Влияние механической обработки на работоспособность композиционных антифрикционных материалов. Тез. докл. 3 Всесоюзного семинара по конструкционным спеченным материалам и изделиям.— Киев: ОНТИ ИПМ АН УССР, 1975.

64. Фляйшер Г. К вопросу о количественном определении трения и износа.— В кн.: Теоретические и прикладные задачи трения, износа и смазки машин. М.: Наука, 1982, с. 285-296.

65. Хрущев М.М, Бабичев М.А. Абразивное изнашивание. М.: Наука, 1970. 252 с.

66. Шилько С.В. Расчетно-экспериментальный метод определения контактных давлений. Динамика и прочность транспортных машин / Сб. научных трудов. Брянск: БИТМ (1994), 116-121.

67. Шилько С.В. Расчетно-экспериментальные методы исследования контактных напряжений // Труды Белорусского конгресса по теор. и прикл. механике. Минск (1995), 248-249.

68. Экспериментальные и расчётно-экспериментальные методы определения контактных напряжений /С.В. Шилько. // Трение и износ. 1996.Т 17. №3. С.402-407.

69. Экспериментальная механика. В 2-х кн. Пер. с англ./ Под ред. А.Кобаяси. Москва. Мир (1990).