Герметичность металл-металлических уплотнительных устройств Текст научной статьи по специальности «Физика»

УДК 621.891 DOI: 10.12737/18269

В.П. Тихомиров, М.А. Измеров

ГЕРМЕТИЧНОСТЬ МЕТАЛЛ-МЕТАЛЛИЧЕСКИХ УПЛОТНИТЕЛЬНЫХ УСТРОЙСТВ

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

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

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

V.P. Tikhomirov, M.A. Izmerov

IMPERMEABILITY OF METAL-METAL SEALING ARRANGEMENTS

A contact of machinery real engineering surfaces has a complex geometry and is defined by errors at machining, in particular, by macro-deviations, corrugation and roughness. This paper report the consideration carried out on a fluid consumption through an end ax-isymmetric sealing arrangement formed by a joint of metal surfaces. Sealed environment consumption is defined through a sum of leaks through a slit between surfaces conditioned on corrugation and dripping through a porous layer farmed by surface roughness on contact contour spots. The surface corrugation is modeled by splines according to initial parameters of a corrugation, and as a roughness a surface fractal model is used repeating a real roughness most accurate. On al-

Введение

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

gorithms presented taking into account basic parameters of corrugation and roughness the 3D models of engineering surfaces are created and after their bringing in contact arises a possibility to obtain practically all parameters of contacting defining the impermeability of a sealing arrangement. The application of 3D models allows rendering processes passing in a contact of conjugate surfaces and obtaining some data, for instance, a tortuosity of leaking channels which earlier were chosen according to recommendations which no doubt increases computation accuracy.

Key words: impermeability, passing, leakage, porous environment, fractal object, compacting.

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

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

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

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

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

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

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

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

Рис. 1. Контурные площадки

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

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

поверхностями

Для осесимметричного кольцевого соединения число волн определяется выражением

п,„ = 2тг rm/WSm,

где гт - средний радиус (;;п = (г, г2}/2); WSm - длина волны (обозначение соответствует стандартам ISO 4287: '97 и ISO 3274: '96).

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

где Жр - высота сглаживания волнистого профиля; индексы 1 и 2 относятся к двум сопряженным поверхностям.

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

к = Н0-3.

Расчет утечек выполним на основе уравнения Рейнольдса для щели с гладкими стенками для неподвижного соединения. Граничные условия: Р = Р\ при г = гг, р = р0 при г = т2. Решая уравнение Рейнольдса с учетом граничных условий, найдем в рамках щелевой модели [2] величину утечки Q для кольцевого зазора с гладкими стенками высотой И:

<1

Здесь Лр=р1-р0 - разность давлений уплотняемой жидкости и атмосферы; ц -динамическая вязкость.

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

равных условиях (перепад давления, динамическая вязкость уплотняемой среды, внутренний и наружный радиусы торцевой поверхности, параметры волнистости) можно найти требуемое сближение поверхностей:

/_ -и1/*

VI/ . — Г"|и

С=1

4тгЛр

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

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

кл

1

V™

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

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

1,а

4(1± 4 12)/ V

4 2,33529

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

5 - Р (Я 4 я

Здесь Ь - длина линии контакта (Ь= г2 ~ г\)\ Д( — {1 — где /и - коэф-

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

2. Полагая, что вершины волн имеют

логарифмически нормальное распределение [4], найдем реакцию 1-й волны Е, соответствующую деформации 37. Сравним сумму реакций, приходящихся на п^ волн, с внешней нагрузкой Е. Если ¿Е >Е, то следует уменьшить сближение: 3тах = 3тах - Л, где Л = 3тах /2. Если ¿¥1 <Е, то следует увеличить сближение: 3тах = 3тах + Л. Если происходит смена неравенства с большего на меньшее или наоборот, то методом половинного деления уменьшаем Л до Л/2. Расчет следует закончить, если выполняется условие

р - "и

-=1 :< [4 р 1 J

Здесь [е] - заданная точность (0,01).

3. Выполним N прогонов моделирования волн (предварительно примем N=20) и определим в каждом случае сближение 3,7=1,...N.

4. По результатам N прогонов вычисляем среднее арифметическое отклонение

и половину доверительного интерва-по формулам

1 *

Здесь дисперсия отклика (сближения)

равна

ОГ

л-

Табличные значения параметра I (критерия Стьюдента) приведены в справочных пособиях по математической статистике. Так, при N=20 и а=0,10 имеем = 2,09.

5. Если отношение

С ОД «))/№))< Г.

то используем <Г(Л/} как точечную оценку и завершаем процедуру прогона для данной нагрузки. В итоге получаем соотношение ~ 5(й0 Погрешность оценим по формуле у'=у/(И-у}, где у - относительная погрешность (0<у<1), доверительный интервал - 100 (1-а) %. В противном случае число прогонов следует увеличить.

6. Изменив нагрузку на волны, найдем в соответствии с предложенной процедурой соотношения Б? ~ Бз ~

По предложенной методике в среде

C++ Builder была разработана программа,

позволяющая смоделировать контакт вол-э

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

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

Моделирование выступов с разной высотой расположения проводилось по логарифмически нормальному закону распределения. Пусть X - случайная величина, распределенная равномерно на отрезке [0,1]. Тогда можно найти зависимость случайной величины (максимальной высоты волны), распределенной по логарифмически нормальному закону, от Согласно [4], после сглаживания данных уравнение регрессии имеет вид

Хм=1,60^10"3ехр(9,78Х)+7,06.

Здесь Xм =Жтах имеет размерность, выражаемую в мкм.

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

С помощью полученного уравнения регрессии можно найти набор вершин волн (Иw■i~Wmaxl2) и провести имитационное моделирование в соответствии с ранее предложенной процедурой.

Расчёты по предложенному алгоритму представлены в виде графиков на рис. 4 и 5. Исходные данные: Wp = 10 мкм, Лp = 0,9 МПа, т2 - т1 = 5 мм, динамическая вязкость п = 0,01 Па-с (вода). Для оценки влияния макрогеометрии на герметичность уплотнения каждый график построен в трех вариантах при разных числах волн.

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

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

• общая длина участка;

• вид волнистости - в ортогональной системе координат или в полярной;

• число волн в соответствующих направлениях по осям или шаг волны;

• амплитуда волны в соответствующих направлениях.

-3 волны ■ 5 волн * 7 волн

Рис. 4. Зависимость утечки Q от нагрузки ^

3 волны И 5 волн волн

Рис. 5. Зависимость утечки Q от радиуса кривизны волны тк

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

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

Рис. 6. Построение профиля волны сплайнами и пример поверхности

Моделируя волны одновременно в двух направлениях (по осям X и Y в ортогональной системе координат или по окружности и радиусу в полярной), можно получить волнистые поверхности с практически любыми геометрическими харак-

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

Фильтрация уплотняемой среды через пористый слой

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

В качестве модели шероховатой поверхности мы используем фрактальную модель [5], позволяющую более точно описать всю сложность структуры поверхности металла после её технологической обработки. Тогда контакт фрактальных поверхностей уплотнительного устройства, являющийся дискретным и состоящим из множества замкнутых и соединённых в цепочки пор, можно представить в виде фрактальной пористой среды. Длина цепочки соединённых между собой пор в плоскости разъёма, по которой проходит утечка, может быть определена (по Ман-дельброту [6]) по следующей зависимости: 1 = Х81~1>г, (1)

где Ь - длина кривой; X - параметр, определяемый экспериментально; 3 - масштаб измерения (раствор циркуля при измерении длины кривой); БТ - фрактальная размерность кривой (цепочки пор).

Фрактальный объект (кривая) обладает свойством самоподобия: любой участок кривой имеет ту же фрактальную размерность, что и вся кривая. С математической точки зрения свойство самоподобия определяется формулой

Здесь к - масштаб увеличения.

Для измерения длины фрактальной кривой в к раз длиннее участка, на котором проводились измерения (стандартный раз-

мер участка - 0,8 х 0,8 мм), необходимо измерять (в частности, методом yard stick) длину раствором циркуля в k раз больше предыдущего. Процедура определения длины фрактальной кривой представлена на рис. 7 (участок AB).

Рис. 7. Определение фрактальной размерности

Длина определяется выражением Ь=Ы(3) 3. С уменьшением размера 3 длина нелинейно возрастает. Прологарифмировав уравнение (1), запишем

1п 1*= ¿тс Л. (1 — £?т) 1п

Зависимость 1пЬ = f(l■гl 5} может быть представлена в виде графика (рис. 8).

График имеет три характерных

участка, разграниченных кроссоверами. В правом участке раствор циркуля 3 боль-

Рис. 8. Определение фрактальной размерности

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

Наклон прямой под углом а на участке 31...32 позволяет найти фракталь-

£?г = 1 — Ьдв = 1 — ■

ную размерность:

1 < От < 2.

Прологарифмировав уравнение (1) при Ь = Ь1 и 3 = 31 и преобразовав полу-

1п ¿± = 1п А +

ченное выражение, запишем

1п А = 1п ¿1 —

,

,

откуда

¿1

Растворы циркуля от 32 до 31 определяют соответствующие длины волн 3 = 1 / а, где а - частота. На рис. 9 представлена спектральная плотность типичного профиля шероховатой поверхности.

Рис. 9. Спектральная плотность

Спектральная плотность профиля фрактальной поверхности определяется зависимостью

5(» =

Здесь в - угловой коэффициент (наклон) прямой зависимости ^ Б(а) от ^ а; с - коэффициент пропорциональности.

Фрактальная размерность оценивается выражением

2

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

е у*

Здесь Ас - контурная площадь исследуемого участка.

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

Магистральный канал протекания длиной Ьт

^ Ас

(контурная площадь)

(инжекция жидкости)

Рис. 10. Протекание через пористую среду

Ь

0

Тогда коэффициент извилистости определяется как

Lr fS

if™ —

= л

С другой стороны, для площадки Ь00 х Ь00 и диаметра канала (2г) имеем [7; 8]

v ° - Z1 - т1^'1*

НгЕ _ г _ 1

О

При DT = 1 Кизв = 1. Учитывая самоподобие (_KasE = Ka.s°\ запишем

L of- i 2 я

1 -BrM

с

гу-°г = X

gX-Dr

д¥гг £>Г-1

откуда средний диаметр канала

Таким образом, средний гидравлический радиус равен

2Л1?Г/С2[1-0Г]>

Расход жидкости через прямолинейный цилиндрический канал, по Пуа-зейлю, определяется уравнением

7ГГ

где ёр - перепад давления. Дифференциал

найдём из уравнения (2), приняв за

переменную величину Ь0. Тогда --£?г

¿г у " Д ^4] ^ О ] £д ] f

L-DJ 2

'С

dLT = к Ас* д1'^ (DT - dLa.

Выражение для утечки через пори-

стую среду имеет вид

-f »>1

q*Cr

о

Здесь гтп и гтах - минимальный и максимальный размеры пор; Ы(г) - число пор, радиус которых находится в диапазоне от г до г+ёг.

X Л г-С*/41)

/ ВШ1

¿'Г.

> Г

и:

Это уравнение получено дифференцированием известного выражения

\Df max > '

Г > '

где Df - фрактальная размерность пор.

Приняв градиент (перепад) давления и подставив приведенные зависимости в предыдущее выражение, найдём

Dtг г ■ г -1 / max

(T*-Df v*~Df) 1

\ тиг miii / 4 _ '

'max 'max mi« J

Здесь п - динамическая вязкость (в расчётах принимают гт;п / гтах = 10" ).

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

— т)

т 4(2 — '

где т - пористость; - площадь пор.

Для осесимметричного уплотнения , откуда (в соответствии

с рис. 11)

Аг

ЬАс =

Ъ - г,

Тогда, обозначив высоту сглажива ния волнистости Жр, найдём

Окончательно запишем:

а) б)

Рис. 11. Пористая среда: а - схема; б - 3Б-модель межконтактного зазора

Фрактальную размерность пористости Б/ можно определить с учётом числа пор фрактальной пористой среды:

N (

= (:

г ч-РГ

Здесь Б/ - фрактальная размерность пористой среды (1 < Б/ < 2). Дифференцируя представленное выражение, запишем

= щ ит

гг.

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

N.

> г

_ (гтах

тсыА Г

В1!Н

Разделив предпоследнее выражение на последнее, получим -¡т в,

К*

= О* г г

Г В11Я

-(»+»/) ¿Г =

гг.

В работе [8] фрактальная размерность пористой среды определяется формулой

1п-т + —--1.

ог = ои

Ь5^

Здесь Бе - евклидова размерность (например, для трёхмерного пространства Бе = 3); т - пористость (0<т<1).

В расчётах принимают тт7п / Ттах = 10" 2 [8]. Тогда фрактальная размерность пористой среды окажется равной

1П7Я

о, = з +-

- 1 2+ 0,215 III т.

Фрактальную размерность извилистости Бт можно определить, представив капиллярную модель протекания через фрактальную пористую среду в виде цепочки соединённых между собой пор. В этом случае капилляр представляет собой фрактальный объект. Фрактальная размерность капилляра характеризует его извилистость, которая, в свою очередь, связана с диаметром капилляра (рис. 12).

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

11-Яг

Рис. 12. Каналы протекания

Прологарифмировав это выражение, получим

1 nLt = DT ln£0 + (l — £>г}1п(2г}, откуда

Тогда

1 L 1 (2 )

Здесь г находится в диапазоне от гм„ до гтах. Значение Ь можно определить, например, с помощью компьютерного моделирования в трёхмерной пористой среде со стороной Ь0 в направлении протекания (рис. 13).

Тогда

Канал прогекания Рис. 13. 3Б-модель пористой среды и канал протекания

- Iяг-1,

Ц.

) lllf

1-ör.

Г ) о С г) С } С г) о С г.

Отсюда найдём фрактальную размерность извилистости />,■ :

ГЦ

1гА

_

LT l 2 г DT = 1 - In — /In —.

III

Ljgf

Средний размер цилиндрических пор (капилляров) можно определить как

V

¿ч>=2гч> = 4*10э- (им), где V - объём пор

(см3); А - площадь их поверхности (м2) в одном грамме.

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

О. мм3/с

0,2 0,15

0,05 О

М)С Л Л

г9

1 г

трансформаторное масло (динамическая вязкость п = 0,0316 Па-с - штриховая линия) и нефть (п = 0,2 Па-с - сплошная линия), при перепаде давления Ар = 106 Па для осесимметричного уплотнения со следующими параметрами: Ь0 = 0,8 мм; Ас = 10 мм2; 5 = 2 мкм; Ь/ = 1,5; Ьт = 1,5.

Полученные зависимости хорошо согласуются с данными литературных источников [10].

О. мм3/с 0,3

0,25 0,2 0,15 0,1 0,05 О

V

l'max- МКМ

1,1

1,6

щ

1,1

1,6

DT

а)

б)

в)

Рис. 14. Зависимость утечки Q: а - от гтах; б - от Ь; в - от Ьт

Вывод

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

можным смоделировать протекание и шероховатых поверхностей с оценкой

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

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

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

1. Greenwood, J.A. Contact of nominally flat surfaces / J.A. Greenwood, J.B.P. Williamson // Proc. Royal Soc. London. Ser. A. - 1966.-V. 293.-P. 300-319.

2. Богомолов, Д.Ю. Математическое моделирование течения жидкости в щелевых каналах с учетом реальной микротопографии поверхности их стенок / Д.Ю. Богомолов, В.В. Поро-шин, В.Ю. Радыгин, А.А. Сыромятникова, А.А. Шейпак. - М.: МГИУ, 2010. - 162 с.

3. Norden, B.N. On the compression of a cylinder in contact with a plane surfaces / B.N. Norden.-Washington: Institute for Basic Standards National Bureau of Standards, 1973. - 47 p.

4. Дианов, А.А. Технологическое обеспечение качества деталей с износостойкими покрытиями за счет управления параметрами точности основы и покрытия: автореф. дис.... канд. техн. наук / А.А. Дианов. - Барнаул: АлтГТУ, 2010.-16 с.

5. Маджумдар, А. Фрактальная модель упруго-

1. Greenwood, J.A. Contact of nominally flat surfaces / J.A. Greenwood, J.B.P. Williamson // Proc. Royal Soc. London. Ser. A. - 1966.-V. 293.-P. 300-319.

2. Bogomolov, D.Yu. Simulation of Fluid Flow in Slit Channels Taking into Account Surface Real Micro-topography of Their Walls / D.Yu. Bogomolov, V.V. Poroshin, V.Yu. Ryndin, A.A. Syromyatnikova, A.A. Sheipak. - М.: MSIM, 2010. - pp. 162.

3. Norden, B.N. On the compression of a cylinder in contact with a plane surfaces / B.N. Norden.-Washington: Institute for Basic Standards National Bureau of Standards, 1973. - 47 p.

4. Dianov, A.A. Technological Support of Parts Quality with Antiwear Coating at the Expense of Accuracy Control of Base and Coating: Abstract of Can.Eng. Degree Thesis / A.A. Dianov. - Barnaul: AltSTU, 2010. - pp. 16.

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

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

7. Kou, J. Fractal analysis of effective thermal conductivity for three-phase (unsaturated) porous media / J. Kou, Y. Liu, F. Wu, J. Fan, H. Lu, Y. Xu / Journal of applied physics. - 2009. - V.106.

8. Yu, B. M. Fractals / B. M. Yu, J. H. Le // Fractals. - 2001. - V.9. - P.365.

9. Bhushan, B. Contact mechanics of rough surfaces in tribology: Multiple asperity contact / B. Bhu-shan / Tribology Letters.-1998.-V.4.-P. 1-35.

10. Тихомиров, В.П. Имитационное моделирование герметичности плоских стыков / В.П. Тихомиров, Л.В. Вольпер // Машиноведение. -1986. - №2. - С. 91 - 94.

5. Majumdar, A. Fractal model of narrow-plate contact of rough surfaces / A. Majumdar, B. Bhushan / Modern Mechanical Engineering. Series. B.-1991.-№6.- pp.11-23.

6. Mandelbrot, B. Nature fractal geometry: [translated from Engl.] / B. Mandelbrot. - M.: Institute of Computer Researches, 2002. - pp. 656.

7. Kou, J. Fractal analysis of effective thermal conductivity for three-phase (unsaturated) porous media / J. Kou, Y. Liu, F. Wu, J. Fan, H. Lu, Y. Xu / Journal of applied physics. - 2009. - V.106.

8. Yu, B. M. Fractals / B. M. Yu, J. H. Le // Fractals. - 2001. - V.9. - P.365.

9. Bhushan, B. Contact mechanics of rough surfaces in tribology: Multiple asperity contact / B. Bhushan / Tribology Letters.-1998.-V.4.-P. 1-35.

10. Tikhomirov, V.P. Simulation of Hermetically Flat Joints / V.P. Tikhomirov, L.V. Volper // Engineering Science. - 1986. - №2. - pp. 91 - 94.

Статья поступила в редколлегию 3.06.2015.

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

Погонышев В.А.

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

Тихомиров Виктор Петрович, д.т.н., профессор, зав. кафедрой «Детали машин» Брянского государственного технического университета, e-mail: dm-bgtu@yandex.ru.

Tikhomirov Victor Petrovich, D.Eng., Prof., Head of the Dep. "Machinery", Braynsk State Technical University e-mail: dm-bgtu@yandex.ru.

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

Izmerov Mikhail Alexandrovich, Can.Eng., Assistant Prof. of the Dep. "Machinery", Bryansk State Technical University, e-mail: m.izmerov@yandex.ru.