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

книги из ГПНТБ / Петров И.К. Технологические измерения и приборы в пищевой промышленности учебник

.pdf
Скачиваний:
30
Добавлен:
25.10.2023
Размер:
18.47 Mб
Скачать

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

Кс=

ХВЫХл

(18)

 

х в х 0

и в общем случае различна в разных точках. Дифференциальная чувствительность

&XBXq "-^вхо

и также в общем случае различна в разных точках.

Для линейных элементов

АГс = /Сд = const.

(20)

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

 

K = K1Kr..Kn_v

(21)

Следовательно, чувствительность измерительного

прибора

с незамкнутой

цепью равна

произведению чувствительности

всех звеньев,

за исключением

отсчетного устройства.

 

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

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

Существует два подхода к решению задачи определения ди­ намических характеристик, в соответствии с которыми уравне­

ния

динамики

составляются

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

а) структу­

ры

 

измерительных устройств

(объектов)

и происходящих

в них

физических и химических явлений; б) реакции

измерительных

устройств на входные возмущения.

 

 

 

 

 

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

тельных

устройств

дают

графики

переходных

процессов

(кривые разгона),

возникаю­

 

 

 

 

щих

при

единичном

скачко­

 

 

 

 

образном

изменении

 

(возму­

 

 

 

 

щении)

входной

величины

 

 

 

 

(лгвх=1) и нулевых

начальных

 

 

 

 

координатах (рис. 6,а) .

Функ­

 

 

 

 

ция,

определяющая

изменение

 

 

 

 

выходной

величины

во

време­

 

 

 

 

ни

при единичном возмущении

 

 

 

 

на

входе,

называется

п е р е -

 

 

 

 

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

XSt,

 

 

 

 

 

 

Хвх9

X*,

 

 

 

 

 

 

Рис. 5. Статические характеристики

Рис. 6. Виды переходных функций при

релейных звеньев.

 

 

 

единичном возмущении.

 

 

х о д н о й ф у н к ц и е й .

Графики переходных

функций

могут

иметь различный вид (рис. 6,6, і и 2), который определяется составом входящих в устройство звеньев.

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

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

гулирования, которые описываются с помощью

передаточных

функций и частотных

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

Поэтому

целесообразно

использовать эти понятия и при рассмотрении

измерительных

устройств.

 

W(p) называется отноше-

П е р е д а т о ч н о й

ф у н к ц и е й

ниє операторного изображения выходной величины к оператор­ ному изображению входной величины (при нулевых начальных условиях), полученным согласно прямому преобразованию Лап­ ласа

 

 

со

 

 

 

 

f ( P ) > I

f{t)e-P(dt,

 

(22)

 

 

о

 

 

 

где f(p) — функция

комплексного

переменного

р,

называемая изображением

оригинала

или оператором

Лапласа;

 

 

f(t) — функция

вещественного

переменного

t

(времени), называемая ори­

гиналом;

 

 

 

 

 

р = а + / ю — комплексная переменная.

 

 

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

уравнение измерительного

устрой­

ства получено в виде:

 

 

 

4ых + • • • + ai *в'ых + • • • + « ! *В ыя +

% *вых = К *вх

+

+ b1'xBX+---

+ b i ^ + . . . +

b m 4 f l

(23)

где а и Ь — постоянные коэффициенты;

 

 

п=т,

 

 

 

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

прп+

...

+ a i P

r + ...

+ a l p

+ а0)

х в ы х (р) =

 

= (&mPm +

. . . +

Ь 1 Р 1 +

. . . +

b l P +

b0)xBX(p).

(24)

На основании этого уравнения в соответствии с определени­ ем передаточная функция устройства

 

Хвых(Р)

= 0

 

 

=

 

 

 

J *вх (0

e~ptdt

 

 

 

 

о

 

 

 

 

= ЬтРт+

••• +ЬіР1+

. . . + b l P

+

b0 =

Q(p)

 

а п р п +

. . . +аср1+

. . . + a l P

+

a0

Р (р) '

(

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

W(0) = K = — .

(26)

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

*вх (О = *вх0 sin (Grf -J- ф в х ) .

(27)

Уравнение (27) может быть выражено в комплексной форме

*вх (0 = *вх0 <?Л и '+ф вх> ,

(28)

где со— круговая частота

колебаний;

 

Фвх — фазовый сдвиг.

,

 

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

Выходные колебания для линейных систем также будут гар­ моническими синусоидальными:

*вых (/) =

х В Ы Х о sin (oaf + ф в ы х ) .

(29)

Следует иметь в виду,

что указанное соотношение зависит

только от частоты колебаний.

 

Различают несколько видов частотных характеристик:

1) амплитудно-частотная (амплитудная),

представляющая

собой отношение амплитуды колебаний на выходе к амплитуде

колебаний на входе

при прохождении

гармонического

сигнала

через

измерительное

устройство

при

разных частотах

(от О

до оо) :

 

 

 

 

 

 

*вых

 

 

 

 

А (со) =

(ш);

(30)

 

 

 

*вх

 

 

2)

фазо-частотная

(фазовая), представляющая собой

зависи­

мость сдвига фаз колебаний между выходной и входной гармо­

никами от круговой частоты

(0<.<о<;оо):

 

 

Ф (ю) =

Фвых («>) — Ф и (ю);

 

(31)

3) амплитудно-фазовая,

объединяющая

обе предыдущие

ха­

рактеристики — амплитудную и фазовую. Аналитически она

оп­

ределяется отношением преобразованной по Фурье выходной ко­

ординаты дгв ы х

(/) к входной

xBX(t):

 

 

 

J

*вых ( O e - ^ ' + W

at

 

 

W (/со) =

л

 

= /?(©) + / / (ш) = А (<о)вМю >,

(32)

 

 

j

*вх (0 e - ' ^ + ' W dt

 

 

 

0

 

 

 

где

R(<£>)

и /(со) — вещественная

и мнимая части амплитудно-фазовой

ха­

 

 

 

рактеристики;

 

 

 

 

А(к>) — амплитудно-частотная характеристика;

 

 

 

ср(со) фазо-частотная

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

 

Частотные характеристики связаны между собой, кроме ука­ занных, следующими соотношениями:

А (со) = VlR (*>)Р + [/ И Р ;

(33)

Ф (со) = arctg і? (и)

(34)

Амплитудно-фазовая характеристика может быть также по­ лучена путем подстановки вместо р переменного /со в выраже­ ние передаточной функции:

U

Р(/со)

а „ ( / с о ) " + . . .

+ a f ( / ( B ) ' + . . .

+ a i / c o + a 0

Так как W(ja>) является дробно-рациональной функцией, ее модуль равен отношению модулей числителя и знаменателя. На основании сказанного

Л(со)

= IQ (/со) |

_ У"(&0&2to2H-fe4CQ4— • • .)2+(6iCQ-&8to3+&5co5-

• • - ) 2 .(36)

 

1 Р ( / с о )

1

V(ao—а2со24сй*

— . . .)2+(f l l co—as cus +aB co5—.. . ) а

 

 

 

Ф (со) =

arg Q (/со) — arg Р (/со) =

arctg

бхсо — 63 со3

+

• • •

 

 

 

 

V •62co2

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

• arctg

a

3 coj +

• • •

 

 

 

 

 

(37)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a 0 — a 2

M 2 +

 

 

 

 

 

 

 

 

Амплитудно-фазовая характеристика представляет собой го­

дограф, описываемый концом вектора

W(ja> ) при изменении ча­

стоты от 0 до оо. Он строится на плоскости комплексного

перемен­

ного

(рис. 7). Каждой частоте

соответствует

определенная

точ­

 

 

 

 

 

 

ка

 

на

 

годографе,

 

j |

 

 

 

 

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

век­

 

 

 

 

 

 

тор,

проведенный

из

 

 

 

 

 

 

начала

координат

в

 

 

 

 

 

 

какую-либо

точку,

 

 

 

 

 

 

характеризует

соот­

 

 

 

 

 

 

ношение

амплитуд

и.

 

 

 

 

 

 

сдвиг фаз при часто­

 

 

 

 

 

 

те,

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

 

 

 

 

 

 

этой

точке.

 

 

 

 

 

 

 

 

 

 

Таким

 

образом,

 

 

 

 

 

 

динамические харак­

 

 

 

 

 

 

теристики

 

могут

 

 

 

 

 

 

быть

получены

не­

 

 

 

 

 

 

сколькими

способа­

Рис. 7. Амплитудно-фазовая характеристика.

ми и представлены

в

разной

форме.

При­

 

 

 

 

 

 

веденные

соотноше­

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

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

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

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

 

(01 — 92 ) a

Fdt =

cmdQ2,

(38)

где «—коэффициент теплопередачи

от среды к гильзе, Вт/(м2 -К);

F — площадь поверхности гильзы, м2 ;

 

 

с—удельная

теплоемкость материала

гильзы,

Дж/(кг - К);

т — масса гильзы, кг.

 

 

 

Из уравнения (38) следует, что

 

 

 

 

с/л

d9,

(39)

 

« і — "2

= ' aF

dt

 

 

Приняв за начальные условия

 

Защитная

 

 

 

 

9V

е 2 п = о,

 

 

гильза

получаем

( 0 i - O i o ) - ( 0 2 - 02 „)

~aF Ш

Примем

ЛЄі: :Bl— 9; А9» = 9 , - 9 ,

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

0 i o = \ -

 

 

Рис. 8. Схема термопары в од­

В результате получим

 

 

нослойной защитной гильзе.

 

 

 

 

 

 

 

 

А9 а

 

А 9 ^ 6 Г

А9 2

cm

S

(42)

9

9,,

aF

dt

 

J4

~*о

 

Введем следующие обозначения:

A9j.

Хвх = "7;— —входная величина;

К

Ав,2

*вых= а— — выходная величина;

К

cm

Г = ^ 7 — коэффициент, называемый постоянной времени, с.

При подстановке этих значений уравнение

(42) записывается

в виде

 

 

1 В Х — '

j a Т Л В Ы Х -

(43)

 

dt

 

25

kn Защитные
гильзы

Или, принимая

p — d/dt,

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

стику термопары, приведенную к операторной форме:

 

 

*вх =

( 7 > +

1)*вых-

(44)

Динамическая

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

может быть выражена

также

в виде передаточной функции:

 

 

Тр + 1

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

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

быть вторым, третьим и т. д.

 

 

 

 

Перейдем

к рассмотрению

ди­

Рис. 9. Схема термопары в двух­

намической

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

тер­

мопары, заключенной

в двухслой­

слойном чехле.

 

 

 

 

 

 

 

ный

чехол (рис. 9). Запишем

 

 

 

 

dt

вых^.

 

 

(46)

 

 

 

 

 

 

С другой стороны,

 

 

 

 

 

 

 

 

 

„ ^*вых .

 

 

(47)

 

 

 

• П —

+ * в

 

 

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

значением хв

,

необходимо

установить

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

этим значением и значением хвх,

характеризующим

температуру

среды. Путем несложных преобразований получаем

выражение

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

термопары в

двухслойном

чехле:

 

 

 

 

 

 

 

 

ЛВХ

' 1 * 2

dt2

( Г 1 +

Г 2 ) % 1 .

 

(48)

 

 

 

 

at

 

 

 

 

 

 

 

 

 

•"•вых^;

 

(49)

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

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

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

тельного устройства, переходные

функции

их имеют различ­

ный вид.

 

 

Для обработки переходных функций с целью получения соот­

ветствующих дифференциальных

уравнений

применяется ряд

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

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

W (р) = е~Рх.

(51)

где т — время чистого запаздывания, с.

При введении в правую часть дифференциального уравнения измерительного устройства значения времени чистого запазды­ вания ( т ^ О ) уравнение принимает вид:

(п)

+

. . . + а .

+

. . . +

«"

+ ап х.

=

а х,

1

' • '

1

"г"вых

і

••' і

""вых

1

"О" вых

 

"л ~вых

 

 

 

 

= Ъп х„

(52)

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

W(P)

ь т

Р т + ... +Ьср1 + ... + b l P + b0 -рт

(53),

 

ап

Р" + • • • + at Р1 + • • • + % р + а0

 

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

 

 

 

 

 

 

W(P)=

Тр+

*

,

,

 

 

 

 

 

 

 

(54)

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W(p)

= (TlP+l)(T2p

К

 

 

 

 

 

 

 

 

 

(55)

 

 

 

 

 

 

 

+

l)'

 

 

 

 

 

 

или при наличии чистого

запаздывания

 

 

 

 

 

 

 

 

 

 

 

 

W(p)-

 

 

К

 

 

 

 

 

 

 

 

 

 

(56)

 

 

 

 

(TiP+

1) (Т2 р +

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

Ті

и Г 2 -

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

называемые

постоянными

времени,

с;

 

 

 

 

К-

передаточный коэффициент

объекта.

было

сказано выше,

 

 

 

 

 

 

 

 

 

 

 

Как

 

 

(

г

'

 

 

 

 

имеется ряд методов

опреде­

 

 

 

 

Т

 

 

 

ления

динамических

харак­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

теристик измерительных

уст­

йВЗ

 

 

 

 

 

 

 

 

ройств

и их обработки

с це­

 

7 ^

 

I

 

 

 

лью

получения

конкретных

 

 

 

 

 

 

 

 

 

числовых значений

величин,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\

 

 

 

входящих

в

рассматривае­

 

 

 

 

 

 

 

 

 

 

мые

уравнения.

Одним из

 

 

 

 

 

 

 

 

 

 

наиболее простых

и

распро­

 

г

а

 

 

 

 

 

t

 

страненных

является

графи­

 

 

 

 

 

 

 

ческий

метод

нахождения

 

 

 

 

 

 

 

 

 

 

Рис.

10.

Переходная

функция

(кривая

 

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

дифференци­

 

ального

уравнения.

 

 

разгона)

измерительного

устройства

при

 

 

 

скачкообразном возмущении.

 

 

 

 

На

рис.

10

приведена

 

 

 

 

 

 

 

 

 

 

простейшая

переходная фун­

кция апериодического устройства

(звена)

с чистым запаздывани­

ем

и показан графический метод

определения

времени

 

чистого

запаздывания и постоянной времени, который ясен из

произве­

денных построений. В р е м я ч и с т о г о з а п а з д ы в а н

и я рав-

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

единичного возмущения

на входе устройства до появления соот­

ветствующего сигнала на выходе. П о с т о я н н а я

в р е м е н и

оп­

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

проекции на ось времени

t отрезка,

от­

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

на уровне

*вых (°°) .

вертикалью, проведенной из точки с (или другой,

проходящей

через произвольно выбранную точку

экспоненты)

и касательной, восстановленной к этой точке. Приближенно ве­ личина постоянной времени апериодического звена может быть определена по переходной функции как время, за которое вы­ ходная величина достигнет значения х в ы х (/) = 0,63 К.

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

 

 

 

 

Т = —

tg

ю (СО)

 

 

 

(57)

 

 

 

 

ё

Ч

М

; .

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

Связь

с

другими

частотными

 

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

имеет

вид:

• = _ 1

/

- L

_ _ i

1 l -

K

l

-

[/(со)]2 - _

1

, /

1

-1

со

V

Я (со)

со

 

 

/(со)

со

К

[Л(с1 (со)]2

 

 

 

 

 

 

 

 

 

 

 

 

(58)

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

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

ва

подаются

гармонические

или

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

колебания

хвх

(t),

а на

выходе регистрируются

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

колебания

*вых

{t)-

Для

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

вычисляются

амплитуды xB X (t)

и хвых

{t) и временные сдвиги

между основной гармоникой входных и выходных колебаний. Эксперимент проводится для нескольких значений периодов (частот) входных колебаний.

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

§5. ПОГРЕШНОСТИ ИЗМЕРЕНИИ

ИИЗМЕРИТЕЛЬНЫХ ПРИБОРОВ. ОЦЕНКА РЕЗУЛЬТАТОВ ИЗМЕРЕНИЙ

Для оценки

достоверности измерения необходимо знать его

п о г р е ш н о с т и

— отклонения результата измерения от истин­

ного значения

измеряемой величины. Различают абсолютную

и относительную

погрешности измерений.

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