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

книги из ГПНТБ / Дмитриев, В. Н. Основы пневмоавтоматики

.pdf
Скачиваний:
25
Добавлен:
22.10.2023
Размер:
17.55 Mб
Скачать

дросселе. Давление в междроссельной камере р\ может быть определено также и с помощью графика, представленного на рис. 46. Аналогичным способом можно получить уравнения для случая, когда первым по потоку располагается турбулентный дроссель (рис. 45, б ) :

—РI а2ри

где

d't У Р2

а, = 0,174

 

.и д /

V RT

и при а2 = 1, V Р о — р \ = р і.

 

равен корню квадратному

Избыток давления над р 2 здесь

из перепада давлений на турбулентном дросселе.

Рис. 48. Графическое построение статической характеристи­ ки пневматического усилителя сопло — заслонка:

а — расходная характеристика

1-го

дросселя;

б

— расходная х а ­

рактеристика 2-го дросселя;

в

— совмещение

расходных характе­

ристик 1-го и 2-го дросселей;

с

графический

способ построения

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

усилителя

 

 

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

90

Po и p2 , то давление в междроссельной камере р\ можно найти,

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

и

2 -го

дросселей

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

р0 = const и

р2 = const в зави­

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

абсциссы точки

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

Другими

словами,

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

к нахождению

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

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

равенству

расходов

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

1-й и 2-й дроссели. Применяя

графическое по­

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

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

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

сил

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

измене­

 

 

нием

скорости

течения во

времени

 

 

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

 

 

значение

расхода

 

принимают

равным

 

 

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

 

 

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

 

 

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

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

 

 

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

 

 

сматривают

как

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

 

 

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

 

 

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

 

 

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

 

 

происходит

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

Рис. 49. Схема

проточной

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

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

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

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

камеры, у

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

средним

между

 

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

ны все входные

величины,

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

Расход­

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

ные

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

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

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

дросселей

 

существенно

нелинейны.

 

 

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

цессе.

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

91

Допустим, что в динамическом процессе могут изменяться все величины, характеризующие состояние воздуха в междроссель­ ной камере, а именно: питающее давление р0, давление па вы­ ходе из камеры ро, эффективные проходные сечения дросселей fi и І2 и объем камеры У за счет перемещения поршня (или уп­

ругого элемента) х. За время dt масса воздуха в камере изме­ няется на величину

 

 

dM = <Ш,—гі/ѴЬ,

 

 

 

(35)

где dM1 — масса

воздуха, поступившего в

междросселыіую ка­

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

 

 

-й дроссель

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

из камеры через 2

за время dt.

 

 

 

по

времени,

будем

Продифференцировав уравнение (35)

иметь

__

 

 

 

 

 

 

d.\t

dM

— G

d&M .

 

(36)

dt

~

dt

dt

 

 

 

 

 

 

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

массовые

расходы

воздуха в

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

 

 

 

 

 

 

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

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

Масса

воздуха

в камере в любой момент времени М = Урь

 

 

 

 

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

 

 

 

 

 

dM = У f/p, + P| dV.

Так как У = öl7 + У,. и щ = f»!t. + брі. то

dM = (6 У + УсѴ/(бр, + Pu) + (бр, + pu)</(бУ + Ус).

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

Пренебрегая приращениями объема 5У и плотности брі по сравнению с Ус и pf. при исходном статическом режиме, получим

dM = ycrföp, + р1 сгібУ,

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

 

 

брі = бр 1 , 6

У = —бxF и р1с = Ріс

RT

 

RT

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

гібД* = ■

■йбщ

P cF döx\

RT

RT

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

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

92

любого другого упругого

элемента

(сильфона,

мембраны n

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

по времени,

получим

 

 

 

 

 

rfö/V! _

Vc

döpt

picF

döx

(37)

~dt

RT

dt

RT

dt

 

С другой стороны, изменение расхода, идущего на изменение давления в камере, является функцией р0, ри р2, fi и h, т. е.

 

 

 

G,—G2

döM

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

dt

 

 

G{Pa, Pi, P ,

 

fь f )-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

 

 

 

döM

dG

&Ро +

dG

 

p,

 

ÖCj

 

e

Ö(J

 

er

 

âG

(38)

 

 

 

 

6

 

 

 

 

/ ? 2

H--- -

 

6

/, +

 

6 / 2

 

dt

дро

dp1

 

 

 

dpi

 

б

 

 

 

âf2

 

 

 

 

 

 

 

 

 

Oh

 

 

 

 

 

 

 

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

сокра­

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

= р0с, р і = Pu-,

г

г

 

 

 

ÜG

 

и т. д.

 

 

 

 

 

 

 

 

1 1 с при производном

----

 

 

 

 

 

 

 

 

 

 

 

Учитывая,

 

 

 

dpО

 

 

 

 

 

 

(37)

и

(38)

равны,

 

что левые части уравнений

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dSp!

plcF döx

 

 

dG

c

 

 

dG

<-

 

 

 

 

 

RT

dt

RT

dt

 

 

- — lVü + ——

6 p, +

 

 

 

 

 

 

dpo

 

 

 

dpK

 

 

 

 

 

 

 

 

dG

<

 

 

â G f - r

dG

 

sc

 

 

(39)

 

 

 

+ T — bP2 + ——

6

/, + ——

 

6

/, .

 

 

 

 

 

 

ap2

 

 

 

dp

 

 

d/ 2

 

 

 

 

 

 

 

 

Поделим

все члены уравнения

(39) на

 

 

dG

P l c

II V M H O -

 

 

döx dpI

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

 

член с

 

dt

 

на А'с,

член

с

ÖP2 на P2

c, член с 6

/і на /,с п

 

член с 6 /

 

 

 

 

 

 

 

2 на /2с

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

- Д й ;

 

вр'

- Д й ;

 

 

 

 

 

 

 

 

 

Рос

 

 

 

 

 

Ріс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- Д й ;

 

б/|

- д / -1

>

 

 

 

 

 

 

 

 

 

Ріс

 

 

 

 

 

J.

-V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/1C

 

 

 

 

 

 

 

 

 

 

 

 

f

 

-А /,

 

 

бд*

 

 

 

 

 

 

 

 

 

 

 

 

 

;

---- = Д*,

 

 

 

 

 

 

 

 

 

 

 

/ 2

С

 

 

 

 

Д'с

д

 

 

 

 

 

 

 

получим дифференциальное уравнение проточной камеры

в

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ту

--- Ь

Л р , =

тх ~~~~ +

 

fco^Pü + k2&P2+

kjAfl +

(40)

93

где постоянные времени и коэффициенты усиления, ственно, равны:

Тѵ

T1 X

=

FXc

d G

d G

 

 

R T -

 

 

R T

dpI

dG

 

dpi

 

 

 

 

dh

hc

 

 

dG

РIC

 

 

dpi

 

 

 

dG

 

 

ko —

,dpt

 

 

dG

 

 

 

 

 

 

dpi

 

 

 

dG

Л c

 

kr, =

dh

 

dG

 

 

dp\

P lc

 

 

 

 

 

dG

 

 

kn

dpo

 

 

dG

 

 

 

 

 

 

dp 1

 

 

соответ-

(41)

(42)

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

dV= —Fdx—xdF

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

d m

- G ,

G2—- Vc_

döp,

 

 

 

 

 

dt

 

RT

dt

 

 

... VT

döp:1

Plc

/

 

R T

d i

RT

' F z

P IC

/ r

dx

 

RT

{ F c ~ d T

+ Xc

+ x c

d F

\

d \ x

d x

) c

d t

 

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

 

F с + Л'с

d F

Тх =

d x / с

 

/?Г- 0 G

 

 

dp I

Причем

в послеДнем

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

ляется как тангенс угла наклона касательной к характеристике F = } {х) в точке, соответствующей исходному статическому режиму.

94

Для примера определим значение коэффициента /?/, для со­

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

для

сочетания режимов истечения НН. Напишем выражение

для

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

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

/ 2

/ о

G = G,—G2 = / , j /

— Pi (Po—Pi) /2| ' -^гР2 (Рі—Pz)

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

 

OG

, /

2

,

,

 

 

і й = V

 

 

 

 

i)G

^|с \ /

RT

^ 0с

 

^!е \ /

^іс

дР\

2

КРіДРос— Pie)

 

2 ( ' ( р , с — р2с)

Воспользуемся выражением (41) для коэффициента kj„ подставив в него найденные частные производные,У*

 

| /

~jjf Р.с о с - Р . с)______________________

 

/ 1 с (Рос

2Ріс) V Ріс

Р2 С

he

І^РгсРіс (Рос

Pi с)

^ lc

 

У РісіРос

РісНріс

Ргс)

 

 

 

Переходя к отношениям

давлений

/"2 = p2J p іс и

г = ргс/Рос

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

 

 

 

 

и

2

(г2 —г) / і —Г,

 

 

 

 

= —f—

Z Z Z I ---------------

■'

 

 

-J- V r2r(r2— г)— (r2-2r) y

' l - г ,

 

 

 

/]

Для определения /г0 при сочетании режимов истечения Н—Н запишем соответствующую разность расходов

G = G i - G 2 = fip„ | / - L - - f aPl ^ - ~ г

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

На основании выражения (42) для коэффициента /г0 можно написать

___ / і С

1

___ / І С

Г 2

/гс

г,

/2С

г

95

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

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

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

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

/И = Е(>.

 

(43)

Дифференцируя уравнение

(43)

по времени и

используя

известное уравнение состояния газа, получим

 

ііМ =

V

d p ,

/4 4 ч

dt

RT

dt

V

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

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

междросселыіую

камеру (или из нее)

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

сМІ

п

п '

----

= Сц — и2.

Переписав уравнение

(44)

dt

 

 

с учетом сказанного, получим

V'

 

—G, —Go

(45)

RT dt

 

 

 

Массовые расходы G\ и G2

можно рассчитать по приближен­

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

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

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

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

и второй

дроссели:

Д Д,

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

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

соответ­

ствующие выражения (4) п

(5), можно получить четыре диф­

ференциальных уравнения,

относящихся к

различным сочета­

ниям режимов истечения. Полученные таким образом выраже­ ния при условии скачкообразного изменения любой из входных величии (ро, 11 , f2 и Р2 ) являются уравнениями первого порядка

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

96

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

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

, А

\

,

 

2г, — і

 

ßißs

 

 

ßf (2г, — І)+ г,

 

t = ---- Carets: —

------

 

arctg

Уг,(I —г,)

 

2/г

 

&2

Уг,(1 —г[)

 

 

 

 

2ß,ß2

 

—arctg

 

2r + ß?-rt

+

Mi_ ln

 

ß3 Ѵ\-Гі+ VWx

х

 

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

2cp

 

 

Рз ѵт=7[- V W i

 

 

 

 

 

V Рз— г + > Г1— г

+ с ,

 

 

 

(46)

 

 

 

 

Vßl—r— \/ гх— г

 

 

 

 

где С — постоянная

 

 

 

 

 

 

интегрирования,

которую

находят

из

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

 

 

 

 

 

 

 

 

 

 

 

cp =

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

ИГ'

 

 

[

 

 

 

 

 

 

 

 

 

 

 

= j /

— ^ + - ^ г + (р;

Р* —і /

2 +

Ш + ф;

 

ß3 =

 

— + ф; ß4 —

 

— Ч— —

 

 

 

V

2

 

 

 

 

 

 

 

2k"-

 

 

 

 

 

2

2£2

 

 

 

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

 

 

 

 

t = A

- г / Щ - .

А

1п

л / - — 1— —

+ С,

(47)

 

У

Г2

 

К

Г,

 

 

J

 

где г2 = Р2 ІР\-

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

t = A-

 

2k arctg j /

 

 

 

 

 

 

 

 

+ с.

4k"-+ 1

 

 

 

 

 

 

 

 

 

 

 

 

(48)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

(49)

 

 

 

 

t = А ln I k г1

1+ С.

 

 

 

 

Вслучае, когда fі или /2 обращаются в нуль, выражения

(46)— (49) вырождаются соответственно в уравнения опо­ рожнения или заполнения глухой камеры.

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

t = —А ------

1+С

7 З а к а з 993

97

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

t = —А In-----1- С.

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

/ = Лі arctg I /

'-+С -,

где

V

т±г 1

 

 

1

V

/ 2

1

І\ | /

RT

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

t —АJгI + С.

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

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

ной камеры. Кроме

того, при

расчете

переходного

процесса,

протекающего в междроссельной

камере,

необходимо

знать

также величины коэффициентов

расхода ц, и р2

дросселирую­

щих органов. Проточная пневматическая камера

в устройствах

пневмоавтоматики

выполняется

чаще всего

в виде

усилителя

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

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

расхо­

да

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

дросселя

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

входящего

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

задается

графически.

Кри­

вые,

выражающие

указанную

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

 

приведены

на рис. 15.

 

 

 

 

 

 

 

 

 

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

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

ный дроссель типа жиклер с Щ= 0,427 мм, переменный дроссель сопло—за­

слонка с с12 = 2,02 мм и объем междроссельной

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

лить временную характеристику рі = /(/) при

изменении расстояния между

соплом и заслонкой скачком от Лс = 0.037 мм до /іуі.т = 0,083 мм. Абсолют­

ное давление питания р0 = 0.25 МПа, абсолютное давление за

вторым дрос­

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

 

Прежде чем приступить к построению переходного процесса, следует опре­

делить давление в междроссельной камере /щ при исходном

статическом

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

98

Найдем значение ріь. Учитывая, что на'исходном

статическом режиме

/і;с = 0,037 мм, а давление в междросселыюп камере

лежит в

пределах

0.25 МПа йз ріс

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

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

считать,

что рі =

—const = 0,8.

 

 

 

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

 

 

:rid?

 

( / 1/ / 2), =

^1 — ^— :

= 1,22.

 

 

 

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

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

определяем уточненное зна­

По графику,

представленному на рис. 15, б,

чение рі 2 ' = 0,48.

 

 

 

 

 

 

Тогда (fi/j2)u

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

приближении

для

Найдем установившееся

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

первом

//уст = 0,083 мм по графику

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

определяем

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

находим г2 = 0,9 и р \ уст

=

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

находим

уточненное

значение .Ug1= 0,755 и рассчитываем (/і//2)ц = 0,29.

По графику па рис. 34 определяем г2 = 0,877 и, следовательно, pt ,-ст =

= 0.114 МПа.

 

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

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

Вычислим также /2 для нового установившегося режима,

р2 = 0,755:

 

/2 = р2ги/2йс = 0,004 см2*.

 

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

те сочетания

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

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

(г2 = 0,5,

Г| = 0.8)

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

замечаем,

что вре­

менная

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

сочетание

режимов

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

намике.

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

ответствующее моменту перехода от сочетания режимов истечения Д — Д кН — Д. Данные расчета сведены в табл. Г.

На рис. 50 представлена экспериментальная кривая переходного процесса. Там же нанесены расчетные точки (по табл. 1).

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

7*

99

 

Соседние файлы в папке книги из ГПНТБ