Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Основы пневмоавтоматики

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
20.51 Mб
Скачать

Ро и /?2, то давление в междроссельной камере р \

можно найти,

используя расходные характеристики 1 -го

и 2 -го дросселей

(рис. 48, а, б), полученные для р0 = const

и р2 =

const в зави­

симости от р\, путем их наложения друг на друга и определения

абсциссы

точки

пересечения (рис. 48, в).

Другими

словами,

способ сводится

к нахождению

точки пересечения расходных

характеристик дросселей, соответствующей

равенству

расходов

в статике

через

1-й и 2-й дроссели. Применяя графическое по­

строение,

можно

также получить

статическую характеристику

пневматического усилителя сопло — заслонка. Такое построение проведено на рис. 48, г.

2. ДИНАМИКА ПНЕВМАТИЧЕСКИХ КАМЕР

Линеаризованное уравнение динамики проточной пневмати­ ческой камеры с турбулентными дросселями [24]. При выводе

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

сил

инерции, обусловленных

измене­

 

 

 

нием

скорости

течения во

времени

 

 

 

в каждой данной точке, и мгновенное

 

 

 

значение

расхода

 

принимают

равным

 

 

 

тому расходу, который был бы в ста­

 

 

 

тике при тех же давлениях до и после

 

 

 

дросселя. Другими

словами, процессы

 

 

 

течения воздуха через дроссели рас­

 

 

 

сматривают

как

квазистатические.

 

 

 

Второе допущение связано с пред­

 

 

 

положением о том, что изменение со­

 

 

 

стояния воздуха в проточных камерах

 

 

 

происходит

по изотермическому зако­

Рис. 49.

Схема

проточной

ну. На самом деле

процесс изменения

пневматической

камеры, у

состояния воздуха в камерах является

которой

могут быть измене­

средним

между

 

изотермическим и

ны все

входные

величины,

адиабатическим процессами.

Расход­

характеризующие состояние

ные

характеристики

турбулентных

воздуха в камере

дросселей

 

существенно

нелинейны.

 

 

 

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

Выведем линеаризованное дифференциальное уравнение проточной пневматической камеры, представленной на рис. 49.

91

Допустим, что в динамическом процессе могут изменяться все величины, характеризующие состояние воздуха в междроссель­ ной камере, а именно: питающее давление р0, давление на вы­ ходе из камеры p2t эффективные проходные сечения дросселей /1 и / 2 и объем камеры V за счет перемещения поршня (или уп­ ругого элемента) х. За время dt масса воздуха в камере изме­

няется на величину

 

йМ = йМх— йМ2,

 

(35)

где dM\ — масса

воздуха,

поступившего в междроссельную ка­

меру за время dt через 1 -й дроссель;

камеры через 2 -й дроссель

dM2— масса воздуха, ушедшего

из

за время dt.

 

уравнение

(35)

по времени, будем

Продифференцировав

иметь

dMx

dM2

г

г

dSM .

 

dM

(36)

dt

dt

dt

 

 

dt

 

 

 

здесь G\ и G2— соответственно

массовые

расходы воздуха в

камеру и из камеры.

Таким образом, расход воздуха, идущий на увеличение дав­ ления в камере, равен разности этих расходов. Масса воздуха в камере в любой момент времени М = Vp\.

Дифференциал этого выражения

dM = V dp{+ PJdV.

Так как V = 6V + Vc и pi = p!c + 6pi, TO

dM = (61/ + Vc)d(89l + pic) + (6pj + pic)d(8V + Vc) .

Здесь и далее значком б будут обозначаться абсолютные приращения величин, отсчитываемых от их значений на исход­ ном статическом режиме, а индексом «с» — значения, относя­ щиеся к исходному статическому режиму.

Пренебрегая приращениями объема 8V и плотности 6pi по сравнению с Vc и рс при исходном статическом режиме, получим

dM = Vcd8p{ + plcd8V,

учитывая, что

 

 

 

 

брх=

, 8V = — 8xF и р1с =

RT

^

RT

 

приводим правую часть равенства к виду

 

d8M = - Vc - d8p{

PisIL d§x\

 

 

RT

^

RT

 

здесь F — площадь поршня или площадь поршневого действия 1

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

92

любого

другого

упругого

элемента

(сильфона,

мембраны

и

т. д.). Продифференцировав последнее равенство

по времени,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dbM

_

Vc

dbpx

 

p}cF

d&x

 

 

 

(37)

 

 

 

dt

 

RT

dt

 

RT

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С другой стороны, изменение расхода, идущего на изменение

давления в камере, является функцией р0, Ри Рч>f\ и /г, т. е.

 

 

 

 

Gi — G2 = —37— = G(Po, Ри Р%,

fu h)-

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

Следовательно, в линейном приближении правую часть

уравнения (37) можно представить в следующем виде:

 

 

 

d m

dG

Ьро

да

бР1 +

4

^ б р 2 + 4

^ б г 1 + 4

^ б / 2.

(38)

dt

дро

 

др

 

 

dpi

 

 

dh

 

dfi

 

 

 

В этом уравнении для частных производных введена

сокра­

щенная запись, при которой опущены индексы ро =

рос, Рi =

Pic,

/1 = flc при производной

dG

и т. д.

 

 

 

 

 

 

 

 

-----

 

 

 

 

 

 

 

 

Учитывая, что левые

дро

уравнений

(37)

и

(38)

равны,

части

можно записать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ddp1

plcF dSx

 

dG

о

,

dG

*

 

 

 

 

 

RT

dt

RT

 

dt

 

 

' &PoH— :— ^Pi +

 

 

 

 

 

 

dp0

 

dp 1

 

 

 

 

 

 

 

, dG

я

 

dG

6/,

dG

б/2-

 

 

 

(39)

 

 

 

H— :—

&Рч '

df{

dh

 

 

 

 

 

 

dpi

 

 

 

 

 

 

 

 

 

Поделим

все

члены уравнения (39) на

 

dG

p ic

и умно-

dpi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

жим и разделим член с бр<> на рос,

член с

dbx

на xc,

член

с

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

бр2 на р2с, член с 6fi на / 1с и член с б/ 2 на f2c. Переходя затем

к

относительным приращениям и обозначая

 

 

 

 

 

 

 

 

 

 

- ^ - А р о ,

 

- ^ - = Др,;

 

 

 

 

 

 

 

 

Рос

 

 

 

Р1C

 

 

 

 

 

 

 

 

 

 

 

&Р2 = Ар2,

 

б / ,

= АЛ ;

 

 

 

 

 

 

 

 

 

Р2С

 

 

 

he

 

 

 

 

 

 

 

 

 

 

 

 

б/,

: Af^

- Ах,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

хс

 

 

 

 

 

 

 

 

получим

дифференциальное

уравнение

проточной

камеры

 

в

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dbpj +

Api — Тх dkx

ko&Po + k2&P2 +

 

i 4- &f2Af2,

(40)

 

dt

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

93

где постоянные времени

и коэффициенты

усиления,

ственно, равны:

 

 

 

 

 

Tv =

 

T

=

F x c

 

 

 

d G

1 X

 

RT

d G

RT-

 

 

 

dpi

 

 

 

dpi

 

 

d G

 

 

 

kr =

 

d h

he

 

 

 

d G

Pic

 

 

 

 

 

 

 

 

dp i

 

 

 

 

 

d G

 

 

 

*2 =

dp2

r2;

 

 

 

 

d G

 

 

 

 

 

dpi

 

 

 

 

 

d G

 

 

 

k f =

 

dh

he .>

 

 

 

dG

Pic

 

 

dp i d G

dpo

kn =

d G dp\

соответ­

(41)

(42)

Площадь F может быть связана с х функциональной зависи­ мостью F = f(x). Тогда

dV = — F dx— xdF

и на основании уравнения (37) можно записать

d m

— Си

dt

U 1

11

и

 

RT

и 2

— ■

 

RT

 

«о

 

1

dt

d8pi

Pic

( F

dx

+

X C

dF

 

dt

RT

V е~dt~

 

dhx

dt

Р

1C

/' F c + x c

dF '

 

 

 

 

)

X c

 

dt

 

RT

'

dx ,)

C

 

 

Следовательно, постоянная времени

 

Fс ~Ь

dF

 

dx / с

 

Т =

 

 

1 X

d G

 

 

RT-

 

 

 

dpi

 

Причем (-4£-)

в последнем

выражении для Тх опреде­

ляется как тангенс угла наклона касательной к характеристике F = f(x) в точке, соответствующей исходному статическому

режиму.

94

Для примера определим значение коэффициента kfx для со­ четаний режимов истечения Д —Д и значение коэффициента k0

для сочетания режимов истечения Я—Я. Напишем выражение для разности расходов воздуха в камеру и из камеры для соче­ тания режимов истечения Д Д:

G = G,— G2 = fi j^/ - j^ p d p o — p t) — f2^/~ £ j:P 2 (P i— p2)

и возьмем частные производные от этой разности по эффектив­ ной площади f 1 и по р\ и, подставив в эти выражения значения

параметров на исходном статическом режиме, получим

 

 

dG

, /

2

.

 

ч

 

 

 

<9/,

V

RT

Plc{p0c

Plc)’

 

dG _

^1с 1

- W

^ ~

2P ^

f2c j ^ R T Ргс

 

dp\

2

^ Р 1с(Рос—Pic)

2 П рю—Рас)

 

Воспользуемся

выражением

(41)

для

коэффициента

 

подставив в него найденные частные производные,

 

 

 

 

|/

Pic (рос

Pic)

 

 

Ьи =

 

/1 с (Рос

2pic)

У Pic

Р2с

/2с V^P2cPic(Poc Pic)

Pic

 

 

|/

2RT

 

 

К Plc(Poc Pic) (p1C

P2c)

 

Переходя

к отношениям

давлений

r2 = P2clP\c и г = р2с/Рос

на исходном статическом режиме, получим

 

 

^_ _________ 2 (г2—г) У 1 —г2__________

/ г 2г(г2 —Л)—(г2 2 г) У 1 —г2

/1

 

Для определения &0 при сочетании режимов истечения Я—Я

запишем соответствующую разность расходов

О - в , - О , - 1 ,Р » у / - ^ r - U e , у - ± г

и найдем частные производные

dG

_ 1 __.

 

1

дро

2RT

<*Pi

2/?Г

На основании выражения (42) для коэффициента k0 можно

написать

/1C ^

flC Г2

Более подробно линеаризованные уравнения динамики про­ точных пневматических камер разобраны в работе [24]. В этой же работе приведены графики, которые дают возможность оп­ ределить коэффициенты дифференциальных уравнений проточ­ ных камер без проведения предварительных расчетов.

Нелинейные дифференциальные уравнения проточной пнев­ матической камеры с турбулентными дросселями. Если принять

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

Выведем дифференциальные уравнения изменения давления во времени в проточной пневматической камере (рис. 33). Как и ранее, будем считать, что Т0 = Тх = Т2 = Т и что процесс

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

Масса воздуха в междроссельной камере

М = 1/р.

(43)

Дифференцируя уравнение (43) по времени и используя известное уравнение состояния газа, получим

<*М

= V

dpt

/4 4 ч

dt

RT

dt

к }

В уравнении (44) величина —— представляет собой мае* dt

совый расход воздуха в междроссельную камеру (или из нее)

в динамическом режиме, т. е.

= Gx— G2.

 

 

 

dt

 

Переписав уравнение (44) с учетом сказанного, получим

V

dp\

= 0 , — G2.

(45)

RT

dt

 

 

Массовые расходы G\ и G2 можно рассчитать по приближен­

ным формулам (4) и (5)

для

адиабатического

процесса. Как и

в статике, здесь также возможны четыре различных сочетания

режимов истечения через первый

и второй дроссели:

Д —Д,

НД, Д — Н и НН. Подставляя

в уравнение (45)

соответ­

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

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

96

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

Для сочетания Д —Д

t

= -----larctg

2г1— 1

 

2k

6 2 К г,(1 — г,)

2r + p f-r ,

-arctg

2 V (r + p?)(r,-r)

 

P it ^ - D

+ r,

PlP2 arctg

 

2p,p2 //- ,(! -r .)

-M ^ ln

Рз 1^1-Г ,+

Ур4Г, X

2cp

Рз K l- Г , -

14г,

y K p 3 - r+ y v i - r

+ С,

(46)

 

где С — постоянная интегрирования, которую

находят из

начальных условий;

 

 

Для сочетания режимов истечения Н—Д

/ =

л

i /

1

1--*-1п

 

 

(47)

 

 

V

г*

[/ г2

 

 

где г2 =

р2/р1.

 

 

 

 

 

 

Для сочетания режимов истечения Д —Н

 

 

 

t = A

2

 

 

 

j------In | 1/г, — 2

Л / 1 —г,

,+ С.

 

 

 

 

4£2+ 1

 

 

-г.

 

 

 

 

 

 

 

 

 

 

 

(48)

Для сочетания режимов истечения Н— Н

 

 

 

 

 

 

 

t = A \ n \ k —rJ + C.

 

 

(49)

В случае, когда fi

или / 2 обращаются

в

нуль, выражения

(46) — (49)

вырождаются соответственно

в

уравнения

опо­

рожнения или заполнения глухой камеры.

 

 

 

При f I =. О для докритического истечения

 

 

 

7 Заказ 993

,97

и для надкритического истечения

t =

— А

In ----- \- С .

 

 

Г 2

При /2 = 0 для докритического истечения

t = А , arctg

+ С,

где

 

 

L

i =

- l / A

 

fi

У RT'

и для надкритического истечения

 

t = А

+ С.

Отметим, что уравнения переходных процессов в междрос­ сельных камерах зависят не от абсолютных значений давлений р0, ри р2, а от их отношений гь г2, г. Для расчета переходного

процесса в междроссельной камере необходимо заранее знать, какие сочетания режимов истечения будут иметь место во время его протекания (возможны три различных сочетания режимов истечения через первый и второй дроссели), к какому сочетанию режимов истечения относится исходный статический режим и каковы его параметры. Эта задача может быть решена с по­ мощью графика, представленного на рис. 34 и построенного с использованием уравнений статики пневматической и проточ­ ной камеры. Кроме того, при расчете переходного процесса, протекающего в междроссельной камере, необходимо знать также величины коэффициентов расхода |Xi и р,2 дросселирую­

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

сопло — заслонка

(рис. 35). Зависимость коэффициента

расхо­

да

от параметров

дросселя

сопло — заслонка,

входящего

в состав такого усилителя, обычно

задается

графически.

Кри­

вые,

выражающие

указанную

зависимость,

приведены

на рис. 15.

 

 

 

 

 

 

 

 

 

 

Пример 6. Для усилителя сопло — заслонка

(рис. 35), имеющего постоян­

ный дроссель типа жиклер с d\ = 0,427 мм,

переменный дроссель сопло — за­

слонка с d2 = 2 ,0 2 мм и объем

междроссельной

камеры

V = 462 см3, опреде­

лить

временную характеристику

pi =f ( t )

при

изменении расстояния

между

соплом и заслонкой скачком от hc = 0,037 мм до Нуст= 0,083 мм.

Абсолют­

ное давление питания

р0 = 0,25 МПа,

абсолютное давление за

вторым

дрос­

селем р2 = 0,1 МПа.

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

98

Найдем

значение р\с. Учитывая,

что

на исходном

статическом

режиме

hc = 0,037 мм,

а давление

в

междроссельной

камере

лежит

в

 

пределах

0,25 МПа ^

pi с

^

0,1 МПа, берем в первом приближении \i\ = 0,4 (рис. 15, б).

Постоянным дросселем служит жиклер, поэтому можно

считать,

что щ =

= const = 0 ,8 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем отношение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(f 1/ ^2 )i =

M-i

nd\

 

 

 

 

 

 

 

 

 

 

 

 

1А2 ЗХ6?2 ЛС= 1 ,2 2 .

 

 

 

 

 

При помощи графика на рис. 34, учитывая, что г = 0,4, находим г2 = 0,47.

Следовательно, р1с = 0,212 МПа.

 

рис. 15, ~б,

определяем

уточненное

зна­

По графику, представленному на

чение \ i \ l =

0,48.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда (/1//2 ) 11

~ 1 и, следовательно, р\с = 0,2 МПа.

приближении

для

Найдем

установившееся

значение р\ уст. В

первом

Луст = 0,083

мм по графику

на

рис. 15, б находим р, 2 = 0,8

и

определяем

(f1//2)I = 0,218. Пользуясь графиком на рис. 34, находим

г2 = 0,9

и

Р iycT

=

= 0,108 МПа. Снова с помощью графика

на рис. 15, б

находим

уточненное

значение р^ 1 = 0,755 и рассчитываем

(/1//2)и = 0,29.

 

 

 

 

 

 

 

По графику на рис. 34

определяем г2 = 0,877 и,

следовательно,

 

pi уСт =

= 0,114 МПа.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В новом установившемся режиме отношение / 1//2 = 0,29.

учитывая,

что

Вычислим также f2 для

 

нового

установившегося

режима,

р2 = 0,755:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ 2 = \i2nd2hz = 0,004 см2*.

 

 

 

 

 

 

 

Теперь есть все данные для расчета.

Остается определить

те сочетания

режимов истечения

через первый

и второй дроссели,

которые

будут

сущест­

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

Намечая на кривой для г = 0,4 точку исходного статического режима

2 = 0,5,

г! =0,8)

и нового установившегося значения (г2 = 0,877),

замечаем,

что вре­

менная

характеристика будет проходить сначала через

сочетание

режимов

Д Д, а затем Н Д и граничное значение будет г2 = 0,8, что соответствует гх = 0,5. Это следует из равенства г = гхг2, которое действительно также в ди­ намике.

Для первого участка Д) расчет ведут по формуле (46). Постоянную интегрирования определяют из условия, что Г\ = 0,8, t = 0. На втором участке

Д)

расчет следует вести по формуле (47),

причем постоянную интегри­

рования

определяют из условия, что гх = 0,8 при

t = tlmЗдесь t\ —время, со­

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

режимов истечения Д Д

к Н Д. Данные расчета сведены в табл. 1.

 

На рис. 50 представлена экспериментальная кривая переходного процесса.

Там же нанесены расчетные точки (по табл. 1).

 

* Так как ц2 зависит от давления перед соплом, в данном случае наибо­ лее правильным было бы разбить весь диапазон изменения р\ на участки и взять средние значения коэффициентов расхода для каждого участка. Если же ради простоты взять одно значение р»2>то его следует выбирать соответст­ вующим окончанию переходного процесса, так как только в этом случае рас­ четное установившееся значение будет наиболее близко к экспериментально­ му. Физически это означает, что вместо медленно уменьшающегося в период переходного процесса эффективного проходного сечения f2 из-за уменьшения М2»обусловленного падением давления рх (см. график на рис. 15, б), взято его наименьшее возможное значение. Вследствие этого расчетная кривая переход­ ного процесса должна пройти несколько выше экспериментальной.

7*

99

 

 

 

Таблица

1

 

 

Расчетные данные для построения

 

 

 

 

переходного процесса

 

 

 

 

 

дг- Д

 

 

Н - Д

 

 

 

 

 

 

к.

 

 

 

 

 

 

t, с

Г\

ю

t, с

Г2

~\

С О

 

 

 

 

 

 

® а

 

 

°1

с

 

 

 

 

II с

 

 

1

^

 

 

 

 

 

 

 

£

 

 

 

0

0,800

0,200

и = 6 ,7 2

0,800

0,125

 

 

1,13

0,700

0,175

7,25

0,815

0,123

 

 

2,76

0,600

0,150

7,84

0,830

0,120

Рис. 50.

Экспериментальная

3,91

0,550

0,137

8,94

0,850

0,118

кривая

переходного процес­

fj = 6,72

0,500

0,125

12,25

0,877

0,114

са в междроссельной камере

Динамика пневматических камер, содержащих ламинарные

дроссели. При выводе дифференциального уравнения (40) про­

точной камеры не было сделано никаких предположений отно­

сительно зависимости расхода

воздуха через

дроссели камеры

от их

геометрических

и

термодинамических

параметров.

Поэтому уравнение (40)

и выражения для

его коэффициентов

могут

быть использованы также и для камеры,

содержащей

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

G\ = a[(pl — p\), G2= a'2(p2i — pl),

где для капилляра

 

 

 

, _

ndA

 

 

 

 

 

 

 

 

256|хдIRT 9

 

 

для щелевого цилиндрического дросселя

 

 

 

 

а

'

яОб3

 

 

 

 

 

 

=

-----------.

 

 

 

 

 

 

 

24\хлЩТ

 

 

Уравнение

(40) действительно

и для случая

включения

различных

(смешанных) дросселей,

например,

турбулентного

и ламинарного.

При этом

необходимо окончательный вывод

выражений коэффициентов

дифференциального

уравнения де­

лать с

учетом соответствующей

формулы

для

разности

расходов G.

В устройствах пневмоавтоматики проточные камеры с лами­ нарными дросселями чаще всего используют в качестве точных

100