Контактная жесткость плоского стыка Текст научной статьи по специальности «Физика»

УДК 621.891

КОНТАКТНАЯ ЖЕСТКОСТЬ ПЛОСКОГО СТЫКА © 2016 В.П. Тихомиров, М.А. Измеров, Д.Я. Антипин Брянский государственный технический университет

Статья поступила в редакцию 28.03.2016

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

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

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

¿а

Рс = Ра

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

^ = 2,4г(|)2.

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

Е*

■ +

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

G = Rqo-1

4 (D - 2 )1пу

2(0-1)

2(D-2) _ 2(D-2)> max min )

Параметр G по данным Д. Павелеску и А. Тудора [4] изменяется в пределах от G=9,9-10-16 до 1,2-10-2 мкм. Радиус закругления верхней части неровности в соответствии с моделью А. Маджумдара [2, 5] равен

1 а0'2

г = -

d2z

dx2 х=0

л 2GD~r

Е± Е2 где [I— коэффициент Пуассона.

Для оценки деформации шероховатого слоя используем фрактальные представления о поверхности твердого тела. Профиль поверхности опишем функцией Вейерштрасса-Ман-дельброта [2, 3, 5]

^ ^ rD-1 V cos2nYnx

= G° 1 Z y(2-0)n ■

n=n-i

где z(x) - высота профиля; D - фрактальная размерность профиля (1<D<2); у - скейлинговый параметр для определения спектральной плотности (у>1); G -высотный скейлинговый параметр, не зависящий от частоты; yni =1/L, где L — длина выборки профиля.

Тихомиров Виктор Петрович, доктор технических наук, профессор, заведующий кафедрой «Детали машин». E-mail: [email protected]

Измеров Михаил Александрович, кандидат технических наук, доцент кафедры «Детали машин». E-mail: m. izmerov@yandex. ru

Антипин Дмитрий Яковлевич, кандидат технических наук, доцент, заведующий кафедрой «Подвижной состав железных дорог». E-mail: [email protected]

При известных данных и можно определить фрактальные параметры рассматриваемой поверхности G и Б.

Модель Маджумдара в отличие от модели Гринвуда-Вильямсона (Г-В) принимает критерий перехода в следующем виде: при деформации, превосходящей критическую имеет месс-то

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

Рис. 1. Модель выступа фрактального профиля

Высота отдельно взятого выступа определяется выражением

г(х = 0) = к = С°-Ч2-°. Пусть площадь основания выступа на уровне средней линии равна а ь = 2 ш к, а площадь среза на уровне 8 от

вершины выступа а = 2 ш 8. Так как а ь =п (^ , то

а ь = у. Тогда

1

/4ah\2

1=Ы-

Соотношение между деформацией 8 имеет вид

высотой выступа и

2-D

4ah\ 2

_ . /4ап\

Анализ полученного выражения при G=const и й=сапй, показывает, что отношение ^ должно быть

о

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

Единичный контакт. Для единичного упругого контакта запишем основные соот-ношения

£ = гЫ, 5 < 8„. Подставив значение радиуса в формулу Герца и, учитывая соотношение 5 = ^ , получим

4л/7Г

F=-

-EG0-1^3-0^2, 0 < а < а„

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

, . а — ае Fep = Fe + [Fp ~ Fe) -

Для пластического

Fv = аН, а > ар.

^р ^есг

состояния

У-есг — & — О-р-

контакта

Рис. 2. Площади упругого контакта: деформируемая неровность и деформируемое полупространство

Множественный контакт. Для множественного контакта используем размерное распределение площадей пятен контакта относительно площади наибольшего пятна ai [5]. Число пятен контакта, имеющих площадь, которая превосходит а, равно

N(aL > а) =

Тогда размерное распределение площадей пятен контакта (по модулю) имеет вид

dN _D aDJ2

П(а) "й"2й(в+2)/2' Нагрузка, воспринимаемая всеми пятнами контакта, определяется выражением

Г аесг

F = Fe(a')n(a')da' + + Г* V а') п(a')dа' + Г Fp(a')п( а') dа'. (1)

аесг ар

Критерии перехода от одного вида деформации к другому (от упругого к упругопластическому и к пластическому) выражаются следующими

зависимостями [6] :

аесг = 2Аш2(Н/Е)2-, ар = 5,35лг2(Н/Е)2. Для количественной оценки в качестве примера используем экспериментальные данные сочетания, для которого эффективный радиус за-кругления верхней части неровности = 5 0 мкм-, Raef=1,76 мкм.

Примем в соответствии с [6] Rq Jf мкм.

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

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

G = Rqo-

4(D - 2)lny

Г,,2(0-2) _ 2(D-. max min

aD/2

2(0 — 1)

л2С)-1'

Для рассматриваемого сочетания получены £>=1,214 и G=8,758•10-9 мкм=8,758 10-12мм. Приведем в табл. 1 экспериментальные и расчетные данные, характеризующие зависимость сближения от давления при первом нагружении для рассматриваемого сочетания.

Таблица 1. Зависимость сближения от давления

Сближение, мкм Давление (номинальное), МПа

3 6 9 12 15

первое нагружение (эксперимент) 1,8 4,4 6,7 7,5 8,1

сближение за счет волнистости (расчет) 0,6 0,8 1,0 1,2 1,4

деформация шероховатого слоя (эксперимент + расчет) 1,2 3,6 5,7 6,3 6,7

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

1. Метод, основанный на анализе экспериментальных данных (табл. 1). Зависимость деформации от давления в плоском стыке представлена на рис. 3. Прямая линия представляет собой уравнение регрессии

р = 1.9365.

Контактная жесткость шероховатого слоя плоского стыка равна

Лр МПа Кш=-^= 1,936-.

Л5 мкм

2. Аналитический метод. В основу метода положена зависимость (1), в которой переменной величиной принята максимальная площадь пятна . Эта максимальная площадь определяет деформацию наибольшего выступа. Связь деформации и площади приведена для рассматриваемого сочетания в табл. 2. Номинальная площадь контакта Аа принята равной 0,6 мм2. На рис. 4 представлены зависимости давления от деформации шероховатого слоя плоского стыка.

Рис. 3. Зависимость контактного давления от деформации

Уравнение имеет вид

регрессии фрактальной модели р = 1,373 8 + 1.254.

Контактная жесткость шероховатого слоя равна

МПа

Кмп — 1,373 -.

мкм

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

Рассмотрим другое сочетание стальных поверхностей, для которого эффективные параметры шероховатости были следующими: Яа г = 1,23 мкм; г/ = 37 мкм; Гк е/г = 2242 мкм. В табл. 3 приведены расчетно-экспериментальные данные, отражающие связь давления с деформацией стыка и выделенной деформацией шероховатого слоя мкм.

Таблица 2. Основные параметры контактного взаимодействия

Деформация шероховатого слоя 8В, мкм Площадь наибольшего пятна а, мм2 Нагрузка Г, Н Контактное давление р=¥/Аа, МПа

1,5 2,356- 10-4 2,6 4,3

2,0 3,142 • 10-4 2,9 4,8

3,0 4,712- 10-4 3,4 5,7

4,0 6,283- 10-4 3,9 6,5

5,0 7,854- 10-4 4,5 7,5

Рис.

У -

о 1 : з 4 з Контактам деформация 5. ыкм

4. Зависимость «давление - деформация»: - эксперимент; У1 - фрактальная модель; у(х) - уравнение регрессии модели

Уравнение регрессии фрактальной модели имеет вид

р = 1,373 8 + 1.254. Контактная жесткость шероховатого слоя равна

МПа

Кмд — 1,373 -.

мкм

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

Приняв Яц =1,54 мкм, найдем для рассматриваемого сочетания Б=1,195 и 0=1,367 10-13 мм. Эти фрактальные параметры были применены для получения количественной оценки контактной жесткости рассматриваемого плоского стыка.

Таблица 3. Зависимость деформации от давления

Параметр Давление в стыке р, МПа

3 6 9 12 15

деформация стыка 8 Бит 3,6 5,8 8,1 9,8 11,2

деформация шероховатого слоя 8а 3,0 5,0 7,1 8,1 9,8

На рис. 5 (а) приве-дена зависимость «давление - деформация», получен-ная на основании экспериментальных данных и здесь же построена кривая уравнения регрессии. Уравнение регрессии имеет вид

р = 0.7988к.3.

Нормальная контактная жесткость шероховатого слоя для данного сочетания по экспериментальным данным определяется выражением

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

р: МПа

б)

К NR —

faL

J *

J о

fenn(a') da'.

пластическое состояние пятен касания. Произведя интегрирования с учетом выражений для и размерного распределения пятен контакта, авторы работы [6, 8] получили следующее уравнение

К MR —

(3 — DS)(DS — 1)

Рис. 5. Механическая характеристика стыка: а) зависимость «давление - деформация», б) зависимость нормальной контактной жесткости от нагрузки (давления): сплошной линией показан эксперимент, пунктирной - модель

Расхождение растет с увеличением давления. Этот факт подтверждается и исследованием Х. Янга. В работе [7] приведены данные о контактной жесткости плоского стыка стальных поверхностей, имеющих следующие физико-механические и фрактальные параметры: Н=1960 МПа; G=5•10-11мм; модуль упругости Е1 =Е2 = 2,07105 МПа; коэффициент Пуассона ^1=^2=0,29. Янг отмечает, что до 1,2 МПа экспериментальные данные совпадают с расчетом. При увеличении давления в стыке наблюдается расхождение теории с экспериментом, причем это расхождение также растет с увеличением давления.

Определенный интерес для случая преимущественно упругого контакта представляют работы [6, 8]. Для отдельно взятой фрактальной неровности контактная жесткость определяется следующей упрощенной зависимостью

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

При оценке контактной жесткости стыка (множественный контакт) в работе [6, 8] использовалось размерное распределение пятен контакта, предложенное А. Маджумдаром [5]. В этом случае контактная жесткость стыка, выражаемая в Н/мм, определяется выражением

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

Зл/27Г " (2 -DSY Дадим оценку контактной жесткости шероховатого слоя по этой формуле. Пусть, как в предыдущем случае, D=1,195; aL=n- 0.0 3 7 ■ 0.00 1 = 1,162 ■ 1 0 " 4 мм2. Тогда контактная жесткость шероховатого слоя оказывается равной

K*NR = 311,4 Н/мм = 0,31 Н/мкм.

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

Выводы:

1. Из общей деформации стыка выделена деформация шероховатого слоя.

2. Приведена методика определения удельной контактной жесткости плоского стыка как фрактального объекта.

3. Определены фрактальные параметры приведенной поверхности за счет замены контакта двух шероховатых поверхностей на контакт гладкой поверхности с шероховатой, имеющей эффективные параметры. Оценка фрактальных параметров D и G произведена путем решения системы уравнений, содержащих эффективные параметры шероховатости, измеренные в соответствии с ГОСТ 2789-73.

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

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

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

1. Тихомиров, В.П. Механика контактного взаимо-действия плоских волнистых поверхностей / В.П.Тихомиров, ОА.Горленко, М.А.Измеров, А.Н.Прокофьев // «Вестник» БГТУ. 2013. № 4 (40). С. 87-94.

2. Majumdar, A. Fractal characterization and simulation of rough surfaces / A. Majumdar, C.L. Tien // Wear. 1990. V.136. P. 313-327.

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

4. 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.

5. Маджумдар, А. Фрактальная модель упругопла-стического контакта шероховатых поверхностей / А. Маджумдар, Б. Бхушан // Современное машиностроение. Сер. Б. 1991. №6. С.11-23.

6. Крагельский, И.В. Основы расчетов на трение и износ / И.В. Крагельский, М.Н. Добычин, В.С. Комбалов. - М.: Машиностроение, 1977. 526 с.

7. Yang, Н. Modeling and Analysis of Normal Contact Stiffness of Machined Joint Surfaces // International Journal of Control and Automation. 2014. Vol. 7, No.6. Р. 21-32.

8. Buczkowski, R. A fractal stiffness model for elasto-plastic contact analysis in press joint/ R. Buczkowski, M. Kleiber // CMM-2011 Computer Methods in Mechanics. 9-12 May 2011, Warsaw, Poland.

CONTACT RIGIDITY OF THE FLAT JOINT

© 2016 V.P. Tikhomirov, M.A. Izmerov, D.Ya. Antipin

Bryansk State Technical University

In work the method of determination the specific contact rigidity of a flat joint of rough surfaces as fractal object with assessment of fractal parameters of the given surface due to replacement of contact of two rough surfaces by contact of smooth surface with rough, having effective parameters is presented. Influence of fractal parameters is shown and consistent patterns of change of specific contact rigidity from nominal pressure are determined, the corresponding examples are given.

Key words: contact rigidity, contact pressure, roughness, fractal dimension, deformation, microledges, rapprochement

Viktor Tikhomirov, Doctor of Technical Sciences, Professor, Head of the Department"Machine Details". E-mail: [email protected]

Mikhail Izmerov, Candidate of Technical Sciences, Associate Professor at the Department "Machine Details". E-mail: [email protected]

Dmitriy Antipin, Candidate of Technical Sciences, Associate Professor, Head of the Department"Rolling Stock of the Railroads". E-mail: [email protected]