Введение Известные фрактальные модели [3 − 5] основаны на допущении о том, что для геометрического описания поверхности достаточно знать её фрактальную размерность D и параметр фрактальной шероховатости G. При этом фрактальная размерность является фактором, от которого зависит соотношение между нагрузкой и площадью среза выступа. Для оценки параметров множественного контакта используют размерное распределение площадей среза. Модели [8 − 10] учитывают, что соотношение площади среза и несущей способности пятна контакта зависит от того, в каком состоянии находится деформируемая неровность. Размерное распределение площадей пятен контакта (среза при рассмотрении контакта деформируемой фрактальной поверхности с гладкой твердой поверхностью) связывают с той же фрактальной размерностью, что и для поверхности.В данной работе рассматривается случай, когда фрактальная размерность D как поверхности, так и фрактальная размерность, входящая в размерное распределение площадей среза DXY, имеют разные значения. Приведена методика определения фрактальной размерности DXY. Определение фрактальной размерности DXY Полагаем, что число неровностей, площадь которых больше а, определяется как: NA>a=c aLaD2,#1  где aL  − наибольшая площадь среза выступа; D  − фрактальная размерность (обозначим ее через DXY  в случае, когда она отличается от фрактальной размерности профиля).Распределение площадей среза неровностей фрактальной шероховатой поверхности представлено на рис. 1.  Распределение площадей среза шероховатой поверхности плоскостью, параллельной средней плоскости, зависит от структурных особенностей рассматриваемой поверхности. Процедура определения закона распределение площадей среза поверхности потребовала разработать соответствующее программное обеспечение [7] (рис. 2). Анализ площадей среза неровностей на некотором расстоянии от линии выступов позволяет найти количество площадок, их размеры и количество площадей, превышающих заданное значение Ni A>a'i.  Разработанное программное обеспечение дает [9] оценку размеров площадей среза (в условных единицах) (рис. 2).Найденные значения площадей среза неровностей ранжируем в списке Excel по убыванию {1,080205; 0,478377; 0,478377; 0,451371; … }T. Затем определяем число площадей среза, превышающих заданное значение Ni A>ai'.  Этому значению соответствует отношение aL'ai',    гдеaL'- максимальная площадь среза из рассматриваемого списка(aL'=1,080205 у.е.) .В качестве примера полученные площади среза подвергнем сортировке по убыванию (выделена выборка из 57 данных): N A>a'= 3; 5; 10; 15; … T;#2  aL'a'= 2,3;2,7;5,0;9,3;…T.#3  Обозначив через y число площадей среза, а через x – отношение aL'a', выразим их связь в виде графика (рис. 3). Зависимость числа неровностей N A>a=f aLa  характеризуется степенной зависимостью  N A>a'=1,9453aL'a'0,8705.#4 Тогда DXY=2∙0,8705=1,741.  Следует сказать, что фрактальная размерность профиля (поверхности) для рассматриваемой поверхности равна D = 1,2868 (DS = D + 1 = 2,2868).Профиль фрактальной поверхности описывается функцией Вейерштрасса-Мандельброта: где D −   фрактальная размерность профиля;G − параметр фрактальной шероховатости.В соответствии с рис. 4, и, учитывая уравнение профиля выступа, запишем: δ=GD-1a'2-D2, #6  где a'=πr'2. Радиус закругления верхней части выступа равен [9]:R=a'D22πGD-1.#7  Критическая площадь В качестве критерия перехода от пластического состояния выступа к упругому, согласно [1], примем:δc=πH2E2R. #8 С другой стороны: δc=ac'2πR. Учитывая значение радиуса, найдем:  ac'=G2π232 HE2D-1.#9  Зависимость нормальной нагрузки от площади среза в условиях упругой деформации в соответствии с теорией Герца, имеет вид [9]: Fea'=432πEGD-1a'3-D2,#10  где E=1-μ12E1+1-μ22E2-1,  E1, E2, μ1,μ2  − модули упругости и коэффициенты Пуассона сопряженных поверхностей.Для пластической деформации выступа соотношение между нагрузкой и площадью среза запишем так:  Fpa'=Ha',#11  где H  − твердость материала выступа.В случае упругой деформации связь между площадями среза и пятна контакта определяется соотношением a'=2a.  Критерий перехода от пластической деформации выступа к упругой представим соотношением ac'=2ac. В соответствии с [1], запишем: ac=2EH2D-1G2.#12  Множественный контакт Контактное взаимодействие определяется следующим соотношением: Найдем зависимости фактической площади контакта от нормальной нагрузки и сравним результаты для случаев, когда D=DXY и D≠DXY  при следующих исходных данных: D = 1,459; c = 1; твердость Н = 3000 МПа.Ra = 0,8 мкм; Rq=Raπ/2=1,003 мкм;  DXY =1,700 ;Зависимости фактической площади контакта от нормальной нагрузки для разных случаев представлены на рис. 5. В рассматриваемом диапазоне изменения нагрузки принятие одинаковыми фрактальную размерность поверхности и размерность (показатель степени), определяющую размерное распределение площадей среза неровностей фрактальной поверхности D=DXY  описывает зависимость с меньшим угловым коэффициентом по сравнению с зависимостью, когда D≠DXY .  Зависимость относительной факти-ческой площади контакта от нормального давления (отношения нагрузки к номинальной площади контакта) представлена на рис. 6. Номинальная площадь контакта определялась [1] как: ArAa=0,5 erfcdRq2,#19  где erfc…  − дополнительная функция ошибки; d  − зазор между средней линией профиля фрактальной поверхности и положением гладкого жесткого тела (рис. 7).Зависимости сближения от нагрузки приведены на рис. 8. Таким образом можно сделать вывод, что фрактальные модели имеют разное решение контактной задачи для случая одинакового значения фрактальной размерности профиля и распределения пятен контакта, и разного значения. Обсуждение Функциональные свойства поверхности определяются ее качеством – физико-механическими свойствами материала и топографическими особенностями. Для оценки эксплуатационных показателей используют те или иные модели исследуемой поверхности и зависимости, характеризующие контактное взаимодействие как единичной неровности, так и множественного контакта. Фрактальные модели отражают слабонагруженный контакт, когда существенное влияние на параметры контактного взаимодействия оказывают субмикронеровности. Моделирование взаимодействия фрактальных поверхностей дает возможность получить количественные результаты с требуемой точностью. Допущения, принимаемые в модели, должны быть так или иначе обоснованы. Уточненная фрактальная модель, предложенная в статье, отличается от известных фрактальных моделей тем, что размерность DXY (показатель степени в выражении для числа пятен среза на определенном уровне шероховатой поверхности N A>a'=c aL'a'DXY/2) отличается от фрактальной размерности поверхности. Для обоснования этого допущения было разработано программное обеспечение и представлена методика оценки DXY. Выводы Рассмотрена модель контактного взаимодействия слабонагруженного стыка, учитывающая условие, что фрактальная размерность поверхности отличается от показателя степени DXY, определяющей размерное распределение площадей среза неровностей.Предложена методика оценки DXY с использованием разработанного в данной работе программного обеспечения.Уточненная модель потребовала введение размерности (показателя степени) распределения числа площадей среза и нормировочного коэффициента с, аппроксими-рующего зависимость N A>a=f aLa .Приведено сравнение результатов, полученных с помощью предлагаемой и известной фрактальной модели.