Математическое моделирование пульсаций давления и шума ручного механизированного пневмоинструмента

УДК 621.001.63

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПУЛЬСАЦИЙ ДАВЛЕНИЯ И ШУМА РУЧНОГО МЕХАНИЗИРОВАННОГО ПНЕВМОИНСТРУМЕНТА

© 2006 А. Н. Крючков Самарский государственный аэрокосмический университет

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

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

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

- зависимости изменения объемов этих камер от времени;

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

Схема ротационного пневмопривода приведена на рис. 1. В статоре А, имеющем ограниченные углами у^ и ук отверстия С и Б для входа сжатого воздуха из магистрали с давлением Рм и для выхода его в атмосферу

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

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

1. Наполнение;

2. Расширение;

3. Выталкивание воздуха из рабочей камеры.

Для описания процесса сжатия (расширения) получено выражение из первого зако -на термодинамики с учетом постоянства энтальпии и уравнения Клапейрона:

кЯТмётм=У1ёр1 + кр1ёУ1 ,

где к - показатель адиабаты; Я - универсальная газовая постоянная; 01 - температура в

магистрали; У1 - объем камеры; р1 - давление в камере.

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

Рис. 1. Схема ротационного пневмопривода: I,..., IV - номера пластин ротора; 1, 2, 3, 4 - номера камер; А - статор; В - ротор; С - входное отверстие; Б - выхлопное отверстие; е - эксцентриситет; г - радиус ротора; Я - радиус статора; ра - атмосферное

давление; у^ 1 - угол начала впуска; ук1 - угол конца впуска; у^ 2 - угол начала выхлопа; у^ 2 - угол конца

выхлопа; Ь - угол между пластинами; а. - текущий угол поворота; р - радиус-вектор; ( - круговая

частота вращения; р - давление в магистрали

кЯГ, ёт, + dQдл - dQ = Vdp + kpdV .

Заменим в этом уравнении массу сжатого воздуха dmм, поступающего в полость V1 в течение времени dt, соответствующим значением Омрасхода

dm, = О1 dt ,

тогда

кО, ЯТ . dt + dQTp - dQ = Vdp + kpdV .

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

О, = т 11 pi

2 к +1

2 к 1 ( P1 1к ( Pl 1 к

к - 1 ЯТ, 1 Pl 0 1 Pl 0

где т1 - коэффициент расхода отверстий; /1 -площадь отверстия; Тм - температура воздуха в магистрали; Я - газовая постоянная (Я = 287 Дж/кг-К для воздуха при Тм=290 К), и надкритический, при котором имеет место максимальный критический расход воздуха [2]:

О*= 0,0404 pм .

,

Расход воздуха, вытекающего из камер в полость за статором, будем рассчитывать по формулам (5) и (6), но с учетом возможного «обратного течения» воздуха. Это связано с тем, что с полостью за статором могут соединиться сразу три камеры. И возможна ситуация, когда воздух из области повышенного давления полости начнет заполнять камеру. Для того, чтобы учесть возможность " обратного течения", введем безразмерный коэффициент у. На рис. 2 представлен график зави-

симости безразмерной величины у от перепада давления:

У =

2 ё+1 Г „XI Г „ХГ

р± V Ро 0

Р±

V Ро 0

р1

,апее—^ > а* Ро

2 ё++1 р а* ё - а* ё , апее— < а * Ро

Тогда уравнение расхода, используя выражение (7), можно записать в виде

О 01 = т/Ро

2 ё 1

ё -1 ЯТ

У

Работа сил трения определяется выражением

^0ТР = ЛТ (к - 1)Мхйа,

А = г1Ък р 1 /ддЯ

Я - к с 2

циент трения пластины, г - число пластин, I - длина пластины.

С учетом утечек сжатого воздуха из камер повышенного давления в камеры с более низким давлением, а также теплообмена с окружающей средой, для любой /-той камеры можно записать следующее уравнение для определения в ней давления [2]:

^ ^ко!ЯТм + Л Т (к - 1)М! ^ ± ± ко:ЯТг - кО2ЯТг - кР 1 1,

йУ

где М=АС - момент сил трения пластин о статор, А - постоянный коэффициент, со -

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

где О1 - расход, поступающий в камеру;

О2 - расход, вытекающий из камеры; Оу -расход, связанный с утечками между соседними камерами.

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

йР = (кО м ЯТм - ко 2ЯТг )у.

где Ъ и к - толщина и высота пластины, рП - Расход воздуха через каналы выхлопной

плотность материала пластины,/Тр - коэффи- системы определим из уравнения Бернулли

с

Рис. 2. Зависимость безразмерной величины у от перепада давления

dG dt

Dp - X

G2

• P

v F2 •г у

F L

где Dp - перепад давления на длинном канале; X - коэффициент местных потерь; F -

площадь каналов; L - длина каналов; р -плотность рабочего тела.

С целью определения крутящего момента, развиваемого пневмодвигателем, рассмотрим уравнение движения ротора, которое в общем случае имеет вид [3, 4]:

, d w w2 dJ м- I, A

J-+--= M - M1 sign w - M 2,

dt 2 da

где J - момент инерции вращающихся масс, приведенный к ротору; a и w - текущий угол поворота и его угловая скорость; М- движущий момент; М} - момент силы трения; М2 -момент сил полезного сопротивления.

Движущий момент М определяется силой, действующей на i -ю пластину (рис. 1), при вращении ротора по часовой стрелке:

P = (Pi - Рг+1 )(Р, - г)1,

мг = l(pt - рг+1 )(р2 - Г2)/2.

Jdw / dt =

' "z-1

= il

Z

i=1

- M1 signw - M2.

pt RC -e2 sin2 ai - ecosai,

где Яс - радиус статора; а1 - угол поворота .-й пластины.

Зависимости объемов вычисляем по формуле

V = Fl

где F = 0,5 J (р2 - r2 )da .

Радиус-вектор р . может быть найден на основании теоремы косинусов (рис. 1):

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

M „ =

- M1.

где pi ,p+1 - давление сжатого воздуха в 1-й и (+1)-й камерах; pi - текущий радиус-вектор, определяющий положение .-й пластины; г, I- радиус ротора и его длина. Движущий момент М от приложенных в середине пластин сил Р:

I

Диаметр выхлопных отверстий гораздо меньше длины излучаемой звуковой волны, поэтому данные отверстия можно представить как точечный источник, который создает сферические звуковые волны. При этом не учитываются нелинейные процессы в виде слабых ударных волн, имеющих в случае выхлопа РМИ малую интенсивность. С учетом этого допущения звуковое давление в момент времени t на расстоянии г от выхлопного отверстия площадью £ составит

Для упрощения задачи примем постоянными момент инерции и момент сил полезного сопротивления в уравнении (14). Тогда это уравнение примет вид

P (r,t ) =

P о S

д v

i r А

t--

c

'0 0

4p r

д t

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

Из уравнения следует, что звуковое давление зависит от производной по времени скорости исходящего потока, т. е. чем рав-

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

Р(г,1) =

Р о

а

4яр дг

О

Г г >

г--

с

V

-о 0

Ь = Ю ;

где р - эффективное значение реального зву-

кового давления; ро = 2 • 1о-5 Па - пороговое

значение давления для воздушного тракта.

При определении эффективного значения пользуются выражением

Р = — Т

1 Т 2

-1 р 2 (г)йг .

где р - плотность воздуха, выбрасываемого из выхлопных отверстий; q - объемный расход.

Для анализа акустических характеристик машин и оборудования используются логарифмические уровни звукового давления Ь (децибел), поэтому в модели предусмотрен переход от полученных по формуле (22) пульсаций давления к уровням шума в децибелах

[5]:

Схематично основные компоненты математической модели представлены на рис. 3.

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

- закон открытия выхлопных отверстий статора кусочно-линейный;

- отсутствуют утечки по торцам пластин;

- отсутствует теплообмен с окружающей средой;

- распределенностью параметров можно пренебречь.

о

Рис. 3. Структура компонента математической модели пневмоинструмента

Для определения динамических характеристик пневматического инструмента запишем систему уравнений:

={кОЯм +ЛТ (к-1М1 йУ ±

к&гЯТг -ОЯТ-кЦг^У)• г = 1..5

= (кОмаЯТк + ОЯТ2 + кОтЯТкъ + кОттЯТы +

йг

+ОЯТ, -кОЯТ)1 •

У,'

%=(ЩЯТ5 - ЩЯТ)-1 •

йг

= { Т(Р1 -Р++1)(Р2 -г1 )+(Рг -Р)Р -гд2) ]/2 }-1 1,

/=1

р8 д . ,

Р(х,у,г,г) = J--q(t - г/с).

* 4яг дГ 7

В результате расчета получены временные зависимости давления в камерах пнев-модвигателя и его выхлопных полостях, а также крутящего момента. Удобно анализировать процессы изменения давления в камерах пневмодвигателя совместно с графиком изменения объема камеры (рис. 4). Во всех

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

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

Движущий момент изменяется следующим образом. Когда одна из пластин начина-

V, см.

СХ г град

Рис. 4. Зависимость объема камеры и давления в камере от угла поворота

Рис. 5. Зависимости пульсации давления в полости за статором

0.01 ' ' 'Шо5 'иОгП ' ' *0.СИ 15 ' ' 0.012-

Рис. 6. Зависимость движущего момента от времени

ет открывать входное отверстие, сжатый воздух из магистрали поступает не только в камеру, в которой давление должно быть близким магистральному, но и в соседнюю, создавая "отрицательный" момент, действующий на соседнюю пластину в сторону, обратную движению мотора, поэтому в это время момент резко падает (рис. 6). Затем по мере увеличения объема при дальнейшем вращении ротора этот момент начинает расти. Падения момента повторяются периодически через угол ¡3 при прохождении каждой пластиной входного отверстия. Участки стабильного момента приходятся на периоды, когда

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

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

давления в камерах выхлопной системы РМИ. На базе разработанной модели можно учитывать влияние различных конструктивных параметров (объемов камер выхлопной системы, площади и длины соединительных каналов, угловых размеров впускных и выпускных отверстий) на момент, развиваемый гайковертом, и его акустические характеристики. На базе модели разработана программа расчета характеристик РМИ, используемая конструкторами инструментального производства АО «АвтоВАЗ». Данная разработка выполнена в рамках Программы сотрудничества СГАУ и АО «АвтоВАЗ».

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

1. Зеленецкий С. Б., Рябков Е. Д., Ми-керов А. Г. Ротационные пневматические двигатели. - Л.: Машиностроение (Ленинградское отделение), 1976.

2. Герц Е. В. Динамика пневматических систем машин. - М.: Машиностроение, 1985.

3. Герц Е.В. Пневматические приводы. Теория и расчет. - М.: Машиностроение, 1969.

4. Зиневич В. Д., Гешлин Л. А. Поршневые и шестеренные пневмодвигатели горно-шахтного оборудования. - М.: Недра, 1982.

5. Справочник по технической акустике / Пер. с нем. Под редакцией М. Хекла и Х. А. Мюллера. - Л.: Судостроение, 1980.

MATHEMATICAL MODELLING OF PRESSURE PULSATIONS AND NOISE OF A MANUAL MECHANIZED PNEUMATIC TOOL

© 2006 A. N. Krutchkov

Samara State Aerospace University

A mathematical model for analyzing dynamic characteristics of a manual mechanized pneumatic tool is developed. The model makes it possible to assess the influence of various design parameters on the acoustic characteristics and the torque created by the tool in operation at the stage of designing.