Герметичность металл-металлических уплотнительных устройств

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

УДК 621.891

DOI: 10.12737/18269

Герметичность металл-металлических уплотнительных устройств

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

Аннотация

металлический волнистость уплотнительный стык

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

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

Annotation

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 axisymmetric 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 algorithms 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. Контурные площадки

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

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

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

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

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

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

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

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

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



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

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

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

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

Здесь L - длина линии контакта (L= r₂ - r₁);

,

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



2. Полагая, что вершины волн имеют логарифмически нормальное распределение [4], найдем реакцию 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₃ ~

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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



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

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

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

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

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

, (1)

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

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

.

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

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

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

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



Зависимость может быть представлена в виде графика (рис. 8).

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

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

Наклон прямой под углом б на участке д₁…д₂ позволяет найти фрактальную размерность:

Прологарифмировав уравнение (1) при L = L₁ и д = д₁ и преобразовав полученное выражение, запишем

,

,

откуда

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

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

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

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

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

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

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

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

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

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

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

При D_T = 1 К_изв = 1. Учитывая самоподобие ⁰), запишем

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

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

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

где dp - перепад давления. Дифференциал найдём из уравнения (2), приняв за переменную величину L₀. Тогда

Выражение для утечки через пористую среду имеет вид

Здесь r_min и r_max - минимальный и максимальный размеры пор; N(r) - число пор, радиус которых находится в диапазоне от r до r+dr.

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

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

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

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

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

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

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

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

.

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



а) б)

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

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

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

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

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

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

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

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

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

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



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

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

откуда

.

Тогда

Здесь r находится в диапазоне от r_min до r_max. Значение L_t можно определить, например, с помощью компьютерного моделирования в трёхмерной пористой среде со стороной L₀ в направлении протекания (рис. 13).



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

Тогда

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

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

где V - объём пор (см³); А - площадь их поверхности (м²) в одном грамме.

Учитывая изложенное, можно смоделировать пористую среду, через которую будет просачиваться уплотняемая жидкость. Так, на рис. 14 представлены графики утечки таких уплотняемых сред, как трансформаторное масло (динамическая вязкость з = 0,0316 Па·с - штриховая линия) и нефть (з = 0,2 Па·с - сплошная линия), при перепаде давления Дp = 10⁶ Па для осесимметричного уплотнения со следующими параметрами: L₀ = 0,8 мм; А_с = 10 мм²; д = 2 мкм; D_f = 1,5; D_T = 1,5.

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

а) б) в)

Рис. 14. Зависимость утечки Q: а - от r_max; б - от D_f; в - от D_T

Вывод

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

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

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. Маджумдар, А. Фрактальная модель упругопластического контакта шероховатых поверхностей / А. Маджумдар, Б. Бхушан / Современное машиностроение. Сер. Б.?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. Bhushan / Tribology Letters.?1998.?V.4.?P. 1-35.

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

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.

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. - М.: 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.

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

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

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

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