Развитие трибологии Текст научной статьи по специальности «Физика»

УДК 621.891

DOI: 10.30987/article 5c0f8089a49ae4.92243541

В.П. Тихомиров, В.В. Гриб

РАЗВИТИЕ ТРИБОЛОГИИ

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

трении - трибологии.

Ключевые слова: фрактальная размерность, фрактальные поверхности, коэффициент износа, контактная жесткость.

У.Р. тиошкоу, У.У. ОпЬ

TRIBOLOGY DEVELOPMENT

The paper reports the information on the state of the matter of a contact correlation assessment in engineering rough surfaces and their wear taking into account the trend in tribilogy development - modern science on friction.

The existing methods for modeling contact correlation and wear are based on the simplified models of micro-irregularities and do not take into account some significant factors or geometrical peculiarities of the structure of actual friction surfaces. The computer

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

techniques application in contact correlation modeling based on a fractal approach to the formation of a surface layer micro-geometry allows increasing considerably the accuracy of contact parameter assessment and predicting forecasting the resource of a friction (tribo-logical) couple designed and automating a design process and making it more vivid and comprehensible.

Key words: fractal dimensionality, fractal surfaces, wear factor, contact rigidity.

аварий (15 %) и износа материалов (70 %). Решение проблемы повышения долговечности машин во многом связано с грамотным использованием трибологических знаний и развитием представлений о процессах трения и изнашивания [2].

Одним из направлений является применение методологических принципов, основанных на фрактальном описании инженерных поверхностей. В основу положены работы Б. Мандельброта (Mandelbrot) [3], Е. Федера (Feder) [4] и др.

Какие преимущества дают фрактальные представления о контакте инженерных поверхностей? Рассмотрим широко известную модель Гринвуда - Вильямсона (Г-В-модель) [5], которая применима для упругого контакта номинально плоских поверхностей. Модель поверхности представляет собой набор сферических сегментов, имеющих одинаковый радиус закругления верхней части выступов и расположенных на средней плоскости шероховатой поверхности. Модель основана на достаточно четких физических положениях о контактном взаимодействии шероховатых поверхностей при упругом состоянии фрикционных пятен контакта. Принятие постоянного радиуса верхней части выступов упрощает моделирование процесса

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

Устранить недостатки, присущие модели Г-В, удалось А. Маджумдару (Ма-]ишёаг) [6]. Фрактальная модель предполагает, что поверхность обладает самоподобием (часть поверхности отражает весь объект) и скейлингом (часть повторяет свои структурные особенности при разном масштабе измерения). Инженерные поверхности являются мультифрактальными. Под фракталом понимается изломанный объект с дробной размерностью, всегда большей, чем топологическая размерность.

Фрактальная размерность. Фрактальный объект (профиль поверхности) обладает свойством самоподобия: любой участок кривой имеет ту же фрактальную размерность, что и вся кривая. Длина фрактальной кривой (по Б. Мандельброту [1]) определяется по формуле I =

Фрактальная размерность профиля (1<П<2) определяется многими способами [7]. Фрактальная размерность поверхности (по Б. Мандельброту) определяется выражением Пз =0+1. Заметим, что приведенное выражение не имеет пока научного

обоснования.

Фрактальные модели поверхности.

Профиль инженерной поверхности может быть описан уравнением Вейерштрасса -Мандельброта [8]:

7(х) = (О)П-1 ]Г (СОБ2пу"х1 у(2-°)п).

П=П1

Здесь О - фрактальный параметр шероховатости; П - фрактальная размерность профиля; у - масштабный параметр (^>1), f определяет частотный спектр профиля шероховатой поверхности. Масштабный параметр у связан с длиной выборки Ь (расстояние между двумя кроссоверами) соотношением

у71'- = 1/1 .

К параметрам, характеризующим профиль шероховатой поверхности, следует отнести О, П и п. По мнению А. Мад-жумдара, подходящим значением для описания профиля является величина 7=1,5. Нижний предел суммирования в уравнении Вейерштрасса - Мандельброта равен пг — 1п(1/1)/1пу.

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

¿(х)

а)

400

г(х) а

б)

400

40С1

в)

Рис. 1. Профили поверхностей при фрактальной размерности: а - П=1,2; б - 0=1,4; в - П= 1,6

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

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

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

21пу

■й

Ш-Б

Для поверхности эта функция имеет

вид

-вц

Здесь а - частота, мкм-1, эквивалентное значение частоты определяется частотами во взаимно перпендикулярных направлениях:

Чч = ^ +

Параметры Си С определяются из

- (5 - 2

2В3)1дывц-,

< 3.

уравнении

ь

1д5{ы) = 1д—-

2 ту

1 < О < 2; 2 < Оь Для профиля и самой поверхности фрактальная размерность определяется своими угловыми коэффициентами К и К' (наклоном прямой, построенной в координатах ^ - и ^Р - Тогда фрактальная размерность для изотропной поверхности, параметры шероховатости которой можно определить по одной профи-лограмме, равна

5-К О =-.

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

3 -К* 1 2

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

ДГ-1

м

п = 0

™= 1

/2пт\ /2тгт\ хсоб (- -Ь узт -

\ М ) ■ V М ).

На рис. 2 представлена модель поверхности, построенная при следующих данных: д=2,7; К=1; Ы=Ы=3; 0n,ш~Rav[O,л] - случайное число, равномерно распределенное на отрезке от 0 до п; =2,17.

Рис. 2. Модель поверхности 16

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

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

\ /Нк

а)

б)

в)

г)

Рис. 3. Построение выступа микронеровности поверхности: а, б, в - этапы построения, г - результат

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

Функция распределения вероятности определяет величину смещения средней точки и степень изрезанности поверхности. В результате указанных процедур получается модель выступа шероховатой поверхности (рис. 3г).

На рис. 4 приведены поверхности (реальные и модели) элементов пары трения после разных видов обработки.

Поверхность после токарной обработки

Моделирование по Вейерштрассу

Рис. 4. Реальные поверхности и их модели

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

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

Таблица

на трибологические показатели

Узел трения и процесс Амплитудные параметры Sa, Форма распределения Шаговые параметры За* Гибридные параметры 3 Дц, 88С

Подшипники • • • •

Уплотнения • • • •

Контактная жесткость • • • •

Изнашивание • • • •

Контактная выносливость • • • о

Электроконтакты • • • •

Примечание. • - сильное влияние, • - среднее, о - слабое.

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

Пример 1. Для имеющихся инженерных поверхностей требуется оценить нормальную контактную жесткость.

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

Я =

1Р аР?2

Здесь а - площадь рассматриваемого сечения выступа; О - фрактальный параметр, определяемый по формуле

1

[Д-1

С -

Ы

гп-ь

\_20 — 4

где Яд - оценка среднего квадратического отклонения ординат профиля; ш^, -наивысшая и низшая частоты профиля как случайного процесса.

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

гл-4'

1/2

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

ас =

где Зс — критическая деформация [11; 12], /тг кН\2

Здесь

к = 0,454+0,41^; 1/Е = £" .. - модуль упругости и коэффициент Пуассона; Н - твердость деформируемого тела. Выразим критерий перехода от упругого состояния контакта к упругопластическому в виде _ /7ткН\2 а°с/г

~ \ 2Ё~) п2.-о' С учетом критической деформации и соотношения для радиуса запишем:

/пкН\2

п

2Е ) [2Е )

/тткН\2

аг

После преобразования получим:

71

ас =

Основные соотношения между нагрузкой и деформацией. Связь между на-

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

1 3

232, а. = пЕ8, 8<8С.

0 3 413

= -ЕИ Ш^2, = ^Л*3'2,

С ^ к- в и

д

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

4 1 I 2 пткн\1

{ а уУ2 а ас

Площадь пятна контакта равна номинальной.

Жесткость сферической неровности:

гЗР 7Р

Т, 0 - - Д П 1/2 *1/2 **1/2

^-А.Я-Л.*

или

5*< 1,

к0

2 Еа

а<п,.

АаЛ/пас

Здесь Аа - площадь контакта. Множественный контакт. Число пятен контакта относительно максимального пятна определяется законом Корчака

[4]:

на >«) = №.

у а 1

Размерное распределение площадей пятен (по модулю) определяется выражением

, dN(A > а) В §

п{ а) =-= —а, а 2 /.

<1.а 2

Найдем силу, приложенную к сопряжению:

ас

Р = I ^(а'Ма.'^а'.

о

Фактическая площадь контакта равна

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

случае, когда максимальная площадь пятна контакта аь=ас, имеем:

с) \аГ 2 1 сЪ-0 с

п-з з-д

2

Для случая, когда а = ас, запишем: О

Г = Р-,

0 с 3 - О

При а = ас и .0=1,5 полученное выражение соответствует значению критической силы. При другом значении фрактальной размерности критическая сила отличается на величину 0/(3-0). Это связано

а,

О'

а < От = ас.

с тем, что модель фрактальной поверхно-

^ (О+

Л™ = С1 а'-а?а^—1 ¿а' = — а,, а, < аг.

ге 1 £ 2—Л

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

Площадь фактического преимущественно упругого контакта оценивается соотношением

При аь = ас получим АТв = — ас.

Фактическое давление в зависимости от сближения равно

О

д-з з-д

2 „ 2

Рг

с3 -

о

Д 2 —д |

2 2 а:а.

а;

2 - О с ь

Предельное значение давления при

упругом контакте равно

П1

\3 - О/а.

V,

Контактная жесткость фрактальных поверхностей. Контактная жесткость фрактальных поверхностей выражается зависимостью

Г 4Е О

Кы = = а,

У ^„л/т^ -

Я£

О

При а£ К,

ас получим 4Е £1

а.

1/2

" Аа^(0- 1} Для расчета контактной жесткости в качестве исходных данных примем сред-

нее квадратическое отклонение ординат профиля 7?<7=0,0015 мм и фрактальную размерность О = 1,5. Найдем параметр фрактальной шероховатости:

0,094 ™т 2

Критическая площадь контакта, соответствующая предельному упругому состоянию:

3 + 1,5

п15-1 С2

а.

= 7,12 ■ 10~4мм

/п- 0,577-3000^1.5-1

I 2 ■ 105 /

Здесь Н=3000 МПа - твердость деформируемой поверхности; £=105 МПа -приведенный модуль упругости.

Номинальная площадь равна

л 4т/я \ 1,5 — 1/ 4

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

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

Полагаем, что оценка объема изношенного материала производится на основе закона Арчарда [13]:

Аа = (1 - 4 мм2.

Контактная жесткость фрактальной поверхности равна

= — ■ :: ; МП,я/:ч:ч=3,44 МПа/мкм.

40 =^-¿(0,

П

где н> — объемный износ; — коэффициент износа; ^ —нормальная нагрузка; Н -твердость изнашиваемого (менее твердого) тела; I — путь трения.

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

Однако в уравнение Арчарда коэффициент трения не входит в явном виде. Введем в

уравнение силу трения (по Амонтону)

=

где / — коэффициент трения.

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

w(t)

w H

L(t).

Здесь kw* будем называть приведенным коэффициентом износа.

Для многих случаев интенсивность изнашивания Ih связана с коэффициентом трения степенной зависимостью : 1К ос Г, п > 1.

тп i к

f = — -+Р+ к -.

1 HБ и AR

Здесь т0 — напряжение среза фрикционной связи при нулевом контактном давлении; H В — твердость по Бринеллю; р — фактор упрочнения фрикционной связи; к - коэффициент; h — глубина внедрения неровности в полупространство; R -радиус закругления верхней части неровности.

Объем изношенного материала.

Модель Арчарда основывается на следующем допущении: объемный износ равен произведению номинальной площади на сближение поверхностей пары трения. Так, для схемы «палец - диск» номинальная площадь равна ттг2, где г - радиус пальца; h - линейный износ. Таким образом, изношенный объем равен w(t) =nr2h(t).

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

ругом и/или пластическом состояниях. Полагая форму верхней части неровности в виде сферического сегмента, найдем объем износа по формуле 1

= -пк(Ър2 + к2). 6

Здесь к — высота сегмента (линейный износ); р — радиус основания сегмента.

Ввиду того что квадрат радиуса основания сегмента существенно больше, чем высота, выражение в скобках вышеприведенного уравнения упростим, оставив первое слагаемое. Заменив р2. = а, /тг, где и, — площадь пятна кон-

такта, перепишем уравнение в виде

1

1^(0 = -к^.

Высота к износа в соответствии с представлением профиля как фрактальной кривой определяется выражением

/ачС2-В)/2 к; — в0-1 --

5

л*

Тогда объем изношенного материала на отдельном выступе оценивается соотношением

где С — параметр фрактальной шероховатости; П - фрактальная размерность (1<П<2).

Уравнение суммарного износа материала имеет вид

w(t) =

2тг

(2-D)/2

4-Е0/2

Здесь пг — число пятен контакта, определяемое по фундаментальному закону Корчака относительно максимальной площади пятна а^ т.е.

п.

=

а

Оту

d/2

П - фрактальная размерность (1<П<2); О -фрактальная шероховатость.

Приведенный коэффициента износа имеет вид

' "г

к» 2K(2~W

: = 1

H

VjAaL

Перепишем полученное выражение, введя площадь максимального пятна кон-

такта и учитывая, что а, * а, а//.

(4-Д)/2

LI О.,

t* = _^_

пг :=1

н

VjAaL

Сумму вычислим с помощью имитационного моделирования случайных чисел, распределенных по закону F ■; = л'1 Для данной площади aL найдем целое число пятен контакта mod nr. Тогда случайное число в соответствии с

приведенным законом найдем из соотно-1/ъ

шения ха* = xR , где хд - случайное число, равномерно распределенное на отрезке [0,1]. В этом случае имеем: ^ С ^ D)

yaf-Dif2 = xiJrDif2b+ - + x / . 1 ai

2b

mo а П7-

Зависимость износа от пути трения (кривая изнашивания) выражается соотношением

РаГ1

h = к

w HB

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

Развитие молекулярно-

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Карюкин, С. Подход к обеспечению живучести стволов артиллерийских орудий / С. Карюкин, О. Митрохин // Военная мысль. - 2012. - № 1. -С. 72-78.

2. Суслов, А.Г. Инженерия поверхности деталей: монография / под ред. А.Г. Суслова. - М.: Машиностроение, 2008. - 320 с.

3. Мандельброт, Б. Фрактальная геометрия природы / Б. Мандельброт. - М.: Ин-т компьютер. ис-след., 2002. - 656 с.

4. Федер, Е. Фракталы / Е. Федер. - М.: Мир, 1991. - 248 с.

5. Greenwood, J.A. Contact of nominally flat surfaces / J.A. Greenwood, J.B.P. Williamson // Proc. Roy. Soc. London. Ser. A. - 1966. - V. 295. - P. 300-319.

6. Majumdar, A. Role of fractal geometry in roughness characterization and contact mechanics of surfaces / A. Majumdar, B. Bhushan // Trans. ASME: Journal of Tribology. - 1990. - № 112. - P. 205-216.

7. Klinkenberg, B. A Review of Methods Used to Determine the Fractal Dimension of Linear Features/B. Klinkenberg // Mathematical Geology. -1994. - Vol. 26. - № 1. - P. 23-46.

8. Wang, S. A fractal theory of the interfacial tem-

perature distribution in the slow sliding regime. Part I. Elastic contact and heat transfer analysis / S. Wang, K. Komvopoulos // Trans. ASME: Journal of Tribology. - 1994. - № 116. - P. 812-823.

9. Pavelescu, D. On the roughness fractal character, the tribological parameters and the error factors / D. Pavelescu, A. Tudor // Proceedings of the Romanian Academy. Ser. A. - 2004. - Vol. 5. - № 2.

10. Потапов, А. А. Теория рассеяния волн фрактальной анизотропной поверхностью / А. А. Потапов, А.В. Лактюнькин // Нелинейный мир. -2008. - Т. 6. - № 1. - С. 3-36.

11. Kogut, L. Elastic-plastic contact analysis of a sphere and a rigid flat / L. Kogut, I. Etsion // Journal of Applied Mechanics. - 2002. - Vol. 69. - № 5.

- P. 657-662.

12. Chang, W.R. An Elastic-Plastic Model for the Contact of Rough Surfaces / W.R. Chang, I. Etsion, D.B. Bogy // ASME J. Tribol. - 1987. - № 109. - P. 257-263.

13. Archard, J.F. Wear theory and mechanisms / J.F. Archard; M.B. Peterson, W.O. Winer (eds.) // Wear Control Handbook. ASME. - 1980.

14. Суслов, А.Г. К вопросу трения и изнашивания деталей машин / А.Г. Суслов // Трение и износ. -1990. - № 5. - С. 801-807.

15. Шалыгин, М.Г. Изнашивание субшероховатости поверхностей трения в водородсодержащей среде: монография / М.Г. Шалыгин. - М.: Инновационное машиностроение, 2018. - 92 с.

1. Karyukin, S. Approach to support barrel life / S. Karyukin, O. Mitrokhin // Military Thought. 2012. - No.1. - pp. 72-78.

2. Suslov, A.G. Parts Surface Engineering: monograph / under the editorship of A.G. Suslov. - M.: Mechanical Engineering, 2008. - pp. 320. .

3. Mandelbot, B. Fractal Geometry of Nature / B. Mandelbrot. - M.: Institute of Computer Investigations, 2002. - pp. 656.

4. Feder, E. Fractals / E. Feder. - M.: Mir, 1991. - pp. 248.

5. Greenwood, J.A. Contact of nominally flat surfaces / J.A. Greenwood, J.B.P. Williamson // Proc. Roy. Soc. London. Ser. A. - 1966. - V. 295. - P. 300-319.

6. Majumdar, A. Role of fractal geometry in roughness characterization and contact mechanics of surfaces / A. Majumdar, B. Bhushan // Trans. ASME: Journal of Tribology. - 1990. - № 112. - P. 205-216.

7. Klinkenberg, B. A Review of Methods Used to Determine the Fractal Dimension of Linear Features/B. Klinkenberg // Mathematical Geology. -1994. - Vol. 26. - № 1. - P. 23-46.

8. Wang, S. A fractal theory of the interfacial temperature distribution in the slow sliding regime. Part I. Elastic contact and heat transfer analysis / S. Wang, K. Komvopoulos // Trans. ASME: Journal of Tribology. - 1994. - № 116. - P. 812-823.

9. Pavelescu, D. On the roughness fractal character, the tribological parameters and the error factors / D. Pavelescu, A. Tudor // Proceedings of the Romanian Academy. Ser. A. - 2004. - Vol. 5. - № 2.

10. Potapov, A.A. Theory of wave scattering by fractal anisotropic surface / A.A. Potapov, A.V. Lakyunkin // Non-Linear World. - 2008. - Vol.6. - No.1. - pp. 3-36.

11. Kogut, L. Elastic-plastic contact analysis of a sphere and a rigid flat / L. Kogut, I. Etsion // Journal of Applied Mechanics. - 2002. - Vol. 69. - № 5. - P. 657-662.

12. Chang, W.R. An Elastic-Plastic Model for the Contact of Rough Surfaces / W.R. Chang, I. Etsion, D.B. Bogy // ASME J. Tribol. - 1987. - № 109. - P. 257-263.

13. Archard, J.F. Wear theory and mechanisms / J.F. Archard; M.B. Peterson, W.O. Winer (eds.) // Wear Control Handbook. ASME. - 1980.

14. Suslov, A.G. To the problem of friction and wear of machinery / A.G. Suslov // Friction and Wear. -1990. - No.5. - pp. 801-807.

15. Shalygin, M.G. Sub-roughness wear of friction surfaces in hydrogen environment: monograph / M.G. Shalygin. - M.: Innovation Mechanical Engineering, 2018. - pp. 92.

Статья поступила в редакцию 03.10.18. Рецензент: д.т.н., профессор Брянского государственного

технического университета Горленко А. О. Статья принята к публикации 12.11.18.

Сведения об авторах

Тихомиров Виктор Петрович, д.т.н., профессор, зав. кафедрой «Детали машин» Брянского государственного технического университета, тел. 8 (4832) 58-82-12, e-mail: dm-b [email protected].

Tikhomirov Victor Petrovich, Dr. Sc. Tech., Prof., Head of the Dep. "Machinery", Bryansk State Technical University, e-mail: [email protected]. Grib Vladimir Vasilievich, Dr. Sc. Tech., Prof. of the

Гриб Владимир Васильевич, д.т.н., профессор кафедры «Производство и ремонт автомобилей и дорожных машин» Московского автомобильно-дорожного института, тел.: 8 (499) 155-08-95.

Dep. "Production and Repair of Motor Cars and Road Machines", Moscow Auto-Road Institute, phone: 8 (499) 155-08-95.