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

Размещено на http://www.allbest.ru/

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

В.П.Тихомиров, О.А.Горленко, М.А.Измеров, А.Н.Прокофьев

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

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

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

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

торцовый уплотнение герметичность зазор

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

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

Рис. 1. Визуализация волнистой поверхности

где r_m - средний радиус (; W_Sm - длина волны (обозначение соответствует стандартам ISO4287: '97 и ISO3274: '96).

В качестве модели волнистой поверхности примем набор радиально расположенных волн, имеющих в верхней части цилиндрическую форму. При этом учитывается случайный разброс амплитуд волн. При оценке уровня герметичности длина волны W_Sm и ее радиус r_w приняты постоянными. Полагаем, что все вершины волн находятся в упругом контакте с сопряженной поверхностью. Уровень деформации связан с амплитудой волн, имеющей логарифмически нормальное распределение [3, с. 9]. Начальный зазор ненагруженного соединения принимается равным расстоянию между средними линиями профиля волнистости:

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

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

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

Здесь Дp=p₁?p₀ - разность давлений уплотняемой жидкости и атмосферы; з - динамическая вязкость.

Рассмотрим алгоритм решения задачи герметичности торцового уплотнения. Поверхность торцового уплотнения представим в виде цилиндрических волн (рис. 2).

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

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

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

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

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

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

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

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

Здесь L - длина линии контакта (L= r₂ - r₁); где м - коэффициент Пуассона, Е - модуль упругости; r_w - радиус волны; F - нагрузка, приходящаяся на волну, которая деформирована до сближения д.

При предварительно рассчитанной величине сближения д_max деформация i-й волны, согласно рис. 3, оказывается равной

2. Найдем реакцию i-й волны F_i, соответствующую деформации д_i. Сравним сумму реакций, приходящихся на волн, УF_i с внешней нагрузкой F. Если УF_i >F, то следует уменьшить сближение: д_max = д_max - Д, где Д = д_max /2. Если УF_i <F, то следует соответственно увеличить сближение: д_max = д_max + Д. Если происходит смена неравенств с большего на меньшее или наоборот, то методом половинного деления уменьшаем Д: Д = Д/2. Расчет следует закончить, если выполняется условие

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

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

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

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

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

5. Если отношение , то используем как точечную оценку и завершаем процедуру прогона для данной нагрузки. В итоге получаем соотношение F₁ ~ . Погрешность оценим по формуле где г - относительная погрешность (0<г<1), доверительный интервал - 100 (1-б) %. В противном случае число прогонов следует увеличить.

6. Изменив нагрузку на волны, найдем в соответствии с предложенной процедурой соотношения F₂ ~ F₃ ~ На рис. 5 представлена блок-схема приведенного алгоритма. По предложенной методике в среде C++ Builder была разработана программа, позволяющая смоделировать контакт волнистых поверхностей и рассчитать величину сближения и контурную площадь контакта волны с гладкой поверхностью при заданной величине нагрузки или определить величину нагрузки при заданном сближении (рис. 4).

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

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

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

X^LN=1,60?10^-3exp(9,78?X^R)+7,06.

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

С помощью полученного уравнения регрессии можно найти набор вершин волн (h_wi~W_max/2) и провести имитационное моделирование в соответствии с ранее предложенной процедурой.

Результаты моделирования. Результаты расчёта по предложенному алгоритму представлены в виде графиков на рис. 6 и 7. Исходные данные: W_p = 10 мкм, Дp = 0,9 МПа, r₂ - r₁ = 5 мм, динамическая вязкость з = 0,01 Па•с (вода).

Рис. 5. Блок-схема основных модулей программы

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

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

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

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

· число волн в соответствующих направлениях по осям;

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

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

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

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

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

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

Рис. 10. Модель торцового уплотнения

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

Список литературы

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

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

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

Размещено на Allbest.ru