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

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

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

апериодических звеньев в таких приборах, как интеграторы и

дифференциаторы.

Поэтому

они содержат

или постоянные,

или регулируемые

дроссели,

проводимости

которых остаются

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

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

Рис. 51. К расчету глухой камеры:

а

схема камеры; б

переходная

функция

при опорожнении;

в — переходная

 

 

функция при

наполнении камеры

 

дросселями за исходную модель примем камеру,

представлен­

ную

на рис. 44, а,

и будем считать,

что расход

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

дроссель линейно зависит от перепала давления. Поэтому здесь

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

отклонений

от

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

 

 

 

 

Выведем уравнение проточной пневматической

камеры,

со­

держащей

линейные

дроссели.

Для

этого в выражение

(45)

вместо

G1

и G2 подставим Gi = ai(p0 — Pi),

G2 =

a2(Pi— Рг).

После соответствующих преобразований получим

 

 

 

 

 

Т п ' at ' + Pi =

k 0p0 + &2Р2>

 

 

 

 

rp

 

V

— постоянная

времени

апериодического

где ^ п = ~^ц—^—“

звена;

k2 = ---- ——

и k0 = — —---------коэффициенты

усиления,

 

 

сх 1 + 0-2

СЦ+ а2

 

 

 

 

 

меньшие единицы (k0 = 1 — &i).

 

 

0,

то

уравнение

Если проводимость второго дросселя аг =

проточной

камеры вырождается

в уравнение

глухой камеры

(рис. 51, а ):

 

 

 

 

 

 

 

 

 

 

Tn^ -

+ Pi=Po,

 

 

 

(50)

 

 

 

at

 

 

 

 

 

 

101

где

v

RTax

Примем давление ро за постоянный уровень отсчета, а дав­ ление, отсчитанное от этого уровня, отметим чертой — р {*.

Тогда

Разделяя переменные и интегрируя, получим

 

Гп Ша = - / + С ,

(51)

где С — постоянная интегрирования.

Для определения постоянной интегрирования обозначим

давление при t = О через рь. Тогда из выражения (51) следует,

что

С = Тп1прь.

И окончательно для рассматриваемого процесса получим

t

(52)

Р\=Рье ?п

Таким образом, падение давления в камере происходит по экспоненциальному закону (рис. 51,6). Время, в течение кото­

рого давление в камере упадет от рь до ра:

/,в = - Г п1п-& -

(53)

Рь

Формулы (52) и (53) были выведены для случая опорожне­ ния камеры, когда р\ > ро. Если же ро > ри то будет происхо­

дить наполнение камеры. Дифференциальное уравнение при этом будет иметь вид

Т„- ^ j - + Pi=Po-

(54)

at

 

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

Пусть входное давление в момент t = 0 претерпевает скач­

кообразное изменение от 0 до р0. Решение неоднородного диф­

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

* Принимается, что в момент t 0 давление изменяется скачком от

рь ДО 0.

102

часть уравнения — величина

постоянная,

то частное

решение

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

 

 

 

 

 

 

и = В,

 

 

 

 

где В — константа.

 

 

 

 

 

Подставляя В в уравнение (54) вместо ри находим

В = р0.

Корень характеристического

уравнения

г = — 1 /Тп и, следова­

тельно, решение

уравнения

(54)

 

 

 

 

можно записать в виде

 

 

 

 

 

 

t

 

 

р Л Я Л V

[

' Pz

Pi — A- e

гп +ро,

Pi

 

_____1/

 

 

 

 

где А = р0 — произвольная по­

 

| Л Л А Г л

 

 

стоянная, определяемая из усло­

 

1

 

 

 

X

 

 

вия, что при t =

0 р\ = 0.

 

 

 

 

 

 

 

 

 

Таким образом

Рис.

52.

Схема пневматического

 

 

 

 

конденсатора

 

 

Pi = А) (1 — е

7")-

(55)

 

 

 

 

График функции (55) представлен на рис. 51, в. Постоянную времени определяют по графику переходной функции как проек­

цию касательной на линию установившегося значения р[ = /?о.

Диаметр капилляра rf, при котором реализуется заданное значение постоянной времени глухой камеры (рис. 51, а) Ти для заданных объема V и длины /, определяют по формуле

4 Г 128,идУ7

d =

У ■лроГп

полученной подстановкой в выражение для Тп значения

Л^4Ро а ^ ------——.

128,Ид№Г

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

лений,

а также

сообщающиеся

сосуды

типа

жидкостного

U-образного дифференциального

манометра с

одинаковым и

неизменными диаметрами.

 

перепада давлений

При

отсутствии

упора под действием

р\ Р2 перемещение сильфона

 

 

 

 

 

х = — (р1— р2) = — Ар,

 

 

 

 

с,

с,

 

 

где F — эффективная площадь

сильфона;

с{— жесткость

сильфона.

 

 

 

юз

 

 

 

 

 

Площадь поршневого действия сильфона примем равной эффективной площади К Тогда масса воздуха, находящегося

в сильфоне:

 

М = Vop + xFр,

 

где

у 0— начальный объем; р — плотность.

 

 

Дифференцируя последнее равенство по времени при

допущении, что давление р\

постоянное12, получим выражение

для расхода воздуха в сильфон:

 

 

 

dM

п п

dx

 

 

------=

G = Fр ----- .

 

 

dt

 

dt

 

 

Исключая х, получим

 

 

 

 

Q _ ^2Р

d&p _

Q dkp

(56)

 

Ci

dt

dt

 

 

где

С = —— — пневматическая емкость.

 

 

Cl

 

есть аналог тока /, а пере­

Учитывая, что массовый расход и

пад давлений Ар — аналог напряжения {У, и принимая во вни­

мание зависимость, имеющую место в электротехнике:

.

r dU

i =

С ---- ,

 

dt

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

3. СТАТИЧЕСКИЕ ХАРАКТЕРИСТИКИ КОММУНИКАЦИОННЫХ КАНАЛОВ СИСТЕМ ПНЕВМОАВТОМАТИКИ

В системах пневмоавтоматики различают два вида

кана­

лов— короткие коммуникационные каналы,

необходимые

для

соединения струйных или иных элементов

в блоках пневмати­

ческих устройств, и длинные трубопроводы

или каналы

связи

отдельных блоков, агрегатов между собой.

 

глав­

При установившихся режимах работы пневмоканалов

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

1 Произведение плсщади поршневого действия на перемещение донышка сильфона дает приращение объема.

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

104

цию канала, сочетание сечений последовательно включенных коммуникационных каналов и т. д. [25].

Рассмотрим местные потери и потери по длине в коммуника­ ционных пневматических каналах.

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

Потери по длине каналов определяют по известной формуле

 

 

 

 

 

(5?)

где :тр — коэффициент

сопротивления

 

трения;

р — плотность

воздуха; v — средняя скорость течения

потока;

/ — длина кана­

ла; X— коэффициент трения, зависящий

от режима

течения,

формы сечения и относительной шероховатости;

dv — гидравли­

ческий диаметр, который находят как

отношение учетверенной

площади поперечного

сечения канала

к периметру,

т. е. dT =

= 4//П.

 

 

 

dTj4. Для лами­

Напомним, что гидравлический радиус RT =

нарного режима течения (Re ^ 2300)

в

канале круглого сече­

ния коэффициент трения X = 64/Re. При неразвитом турбулент­

ном режиме течения этот коэффициент

определяют

по

формуле Блазиуса: X = 0,316 Re- 1 4 . В общем

случае для труб

круглого сечения для всех режимов течения X

определяют

по

графику Никурадзе [9]. Если форма сечения канала отлична от

круглой, то X зависит не только от

Re, но и от отношения h/by

т. е. X =

f(h/b, Re), где h u b — соответственно высота и ширина

канала.

 

 

Вычислим число Рейнольдса для наиболее часто встречаю­

щихся

каналов — прямоугольного

и круглого сечений. Если

известна средняя скорость течения потока v, то нетрудно найти

числа Re, соответствующие разным формам коммуникационных каналов. Считая, что число Re = vdT/л\ можно найти:

для канала прямоугольного сечения

Re = 2vbh =

2Q

v (6 + /z)

v(b + h)

где v = Q/bh\ Q — объемный расход,

 

Ю5

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

 

4Q

V

d 9

тде

v= 4Q/nd2.

А.ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПНЕВМАТИЧЕСКИХ КОММУТАЦИОННЫХ КАНАЛОВ

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

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

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

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

При

обработке экспериментальных данных было

условно

принято,

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

собой

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

соединение

звена чистого запаздывания и

апериодического

звена, т. е. что выходное давление р\

связано

с давлением на входе в канал ро следующими уравнениями:

 

 

dt

+ р \ = kp'i

(58)

106

где т — время чистого запаздывания; р'0 — промежуточное

давление (рис. 53).

На основании опытных данных для величины т принято сле­ дующее выражение:

/

т = Л — ,

а

где I — длина канала; а — скорость звука в воздухе; ц — коэф­ фициент пропорциональности, равный 1 ,2 , согласно опытным

данным для канала сечением 1,5 X 1,5 мм.

Опыты показали также, что для каналов указанного выше сечения при наличии на их концах дроссельных сопротивлений,

Рис. 53. Упрощенная Ро

сг

Pi

и

Pi

структурная схема для

 

 

 

 

расчета коммуникацией-

 

 

 

Та

ного канала

 

 

 

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

лов струйных элементов, постоянная

времени апериодического

звена Тп = 15,5 • 10- 1 0 1с.

получено

для

каналов

длиной

/ =

Это соотношение было

= 100 -т- 900 мм.

времени

Т\

передачи

сигнала

по

При оценке суммарного

коммуникационному каналу было принято, что

Ti = Tn + т.

5. ПНЕВМАТИЧЕСКИЕ ДЛИННЫЕ ЛИНИИ

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

 

+ RQ + L —

= 0;

 

 

д х

dt

 

(59)

 

W- + C-*P- = о,

 

 

 

 

 

д х

dt

 

 

где R = —Л^ д ;

L = - у ;

С =

/ — площадь поперечно­

го сечения канала;

рд — динамический

коэффициент

вязкости

воздуха; р — плотность воздуха; а — скорость звука.

изменение

Уравнения (59)

получены

в предположении, что

давления в линии мало по сравнению с давлением при исходном

107

статическом режиме и, следовательно, газ можно считать не­ сжимаемым. Кроме того, было принято, что потери на трение определяются так же, как при ламинарном установившемся течении газа или жидкости по трубопроводу.

Систему (59) можно привести к одному уравнению 2 -го

порядка в частных производных относительно объемного рас­ хода Q или давления р :

, d 2Q

dQ

1

d 2Q

= 0

L*

+ R

C

dx\

d t2

dt

(60)

 

 

 

 

д 2р

■LC- d 2p

RC

dp = 0

д х 2

d t2

 

dt

 

Методам интегрирования уравнения (59) и (60) при различ­ ных значениях R, L и С посвящено много работ в области элек­

тротехники, акустики, гидродинамики, теории упругости и т. д. Приведем результаты решения уравнений (59), полученные с применением методов операционного исчислений [24].

Методы операционного исчисления удобны в тех случаях,

когда переход от изображения

к оригиналу

можно осуществить

с использованием табличных

изображений.

Применительно к

каналам пневматических приборов результаты решения уравне­ ний (59) и (60) методами операционного исчисления приведем для нескольких частных случаев.

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

Вобщем виде уравнение, связывающее объемный расход Q

иизменение давления в пневмокамере, имеет вид

Q = CKJP~. dt

Для глухой камеры с постоянным объемом

Ск

_ V

У

 

Ро

 

а для глухой камеры с переменным объемом

СК

Vo

F 2

у

 

Ро

С

где V — постоянный объем камеры; У0 и pQ— объем камеры и давление на исходном статическом режиме; F — площадь пор­ шневого действия;- с — жесткость сильфона или пружины, ко­

торой нагружена мембрана, поршень или другой элемент, изменяющий объем.

108

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

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

что

Ск =

0. Давление на входе изменяется скачком,

в

момент

t =

0 от

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

за

р = 0,

до значения рх=о, которой сохраняется затем постоянным. Опу­

ская

промежуточные

выкладки,,

приведем

 

решение уравнений

(59),

полученное методами операционного

исчисления [24].

Давление в трубе на расстоянии х в некоторый

момент

време­

ни t

 

 

 

 

 

3it

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р(ху () = рх=о

 

 

 

 

е

1

 

 

 

X

 

 

 

а

+

 

Р /

— а — Р/

 

 

 

 

 

 

 

 

X (2/ + 1)я COS

2 / + 1

/— х

Я

,

/ =

0,

1 , 2 , ..

(61)

где

. т

2

2

 

 

 

 

 

 

 

 

 

 

 

R .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(62)

Зная значения R f L и С, по

формулам

 

(62)

можно

рассчи­

тать а и Pj при у = 0, 1 , 2 ,... и для постоянного х вычислить чле­ ны ряда для различных t и, наоборот, при постоянном t вычислить члены ряда для различных х в выражении (61) для давления р(ху t). В первом случае найдем изменение давления

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

2 . Рассмотрим случай, когда к трубопроводу присоединена

пневмокамера, а отношение длины трубопровода к его диаметру невелико. При этом условии можно пренебречь влиянием вяз­ кого трения и принять R = 0. Давление на входе в трубопровод изменяется скачком при t = 0 на величину рх=0 и затем остает­

ся постоянным. Решение уравнений

(59), полученное с помощью

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

следующий вид [24]:

 

 

 

Л х

 

 

 

 

р(х, t) = px=0 1— 4

sin Р/ —

 

Р / _

t

(63)

------------ ;-------- COS

 

 

2ру +

s:n 2,3/

/

у LC

 

 

Значение

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

Р/,

определяемых

уравнением

ctg*P/ = —&У—,

находят по специальным

графикам

и

табли­

цам [24].

 

 

 

 

 

 

 

109

Т а б л и ц а 2

Данные расчета коэффициента

Pi

Р*

Рз

р4

Р.

Р2

Рз

р4

 

0

1,5708

4,7124

7,8540

10,9956

50

0,8057

3,3846

6,4133

9,5127

 

5

1,4451

4,3488

7,2865

10,2639

55

0,7506

3,3478

6,3923

9,4983

 

10

1,3390

4,0879

6,9665

9,9432

60

0,6939

3,3141

6,3735

9,4856

 

15

1,2481

3,9044

6,7860

9,7890|

65

0,6341

3,2827

6,3564

9,4740

:

20

1,1686

3,7712

6,6737

9,7008

70

0,5690

3,2530

6,3405

9,4633

25

1,0977

3,6704

6,5975

9,6436

75

0,4956

3,2245

6,3255

9,4532

 

30

1,0330

3,5910

6,5420

9,6033

80

0,4079

3,1968

6,3111

9,4435

 

35

0,9728

3,5264

6,4995

9,5729

85

0,2915

3,1692

6,2970

914341

I

40

0,9157

3,4722

6,4655

9,5490

90

0,000

3,1416

6,2832

9,4248

45

0,8603

3,4256

6,4373

9,5294

 

 

 

 

 

 

1

В табл. 2 приведены значения Pj для четырех корней урав­ нения (63) в зависимости от

0 = arctg - L ,

где

6 = Cl

с*

При больших значениях / коэффициент PJ можно определять по формуле

б

Р/= (/ ■ -!)* • +(/'—1