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

TechExpose

Цыфровые устройства и системы

Технические разделы

Синтез системы автоматического регулирования

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

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

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

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

Данная работа — первый опыт синтеза системы автоматического управления в моей образовательной деятельности. В ней заключены основные этапы проектирования, анализа и корректировки реальной системы автоматического управления. При выполнении синтеза системы автоматического управления необходимо точно знать цель выполняемых действий, и методы её достижения. Целью действий в моем случае является синтез корректирующего устройства, реализующего корректирующую передаточную функцию, которая, при последовательном включении в систему, преобразует имеющуюся АЧХ системы к желаемому виду.

синтез корректирующее устройство температура пар

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

Рис.1.1 (Функциональная схема нескорректированной системы)

Функциональные элементы системы:

ОУ — объект управления (паропровод);

ЗУ — задающее устройство;

СУ — сравнивающее устройство;

qд° — действительное значение температуры пара;

qЗ° — заданное значение температуры;з — напряжение, выдаваемое задающим устройством;т — напряжение на выходе термопары;

DU = Uз-Uт — напряжение, выдаваемое сравнивающим устройством;- напряжение на выходе усилителя;

aд — угол поворота вала двигателя;

aр — угол поворота вала редуктора;

aз — угол поворота заслонки.

Исходные данные ДЛЯ курсового проектирования по ТАУ

ПЕРЕДАТОЧНЫЕ ФУНКЦИИ ЭЛЕМЕНТОВ СИСТЕМЫ:

Wзу(S) = Kзу — передаточная функция задающего устройства;

— передаточная функция электродвигателя,

Тм — механическая постоянная времени;

Тэ — электрическая постоянная времени;

Кдв — коэффициент передачи двигателя;

Wр(S) = Кр — передаточная функция редуктора,

Кр — коэффициент передачи редуктора;

где- передаточное отношение редуктора;= 200;

Wз(S) = Кз — передаточная функция заслонки,

Кз — коэффициент передачи заслонки;

Кз = 5 0С/угл.град;

Wп(S) = Кп´е-t´S — передаточная функция паропровода,

Кп — коэффициент передачи паропровода;

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

— передаточная функция термопары,

Тт — постоянная времени;

Тт = 0,05 с;у(S) = Ку — передаточная функция усилителя;

ТРЕБОВАНИЯ, ПРЕДЪЯВЛЯЕМЫЕ К СИСТЕМЕ:

. Заданное значение температуры пара qЗ° = 80 0С;

. Величина перерегулирования s = 30 %;

. Время регулирования tp = 0,7 с;

. Установившиеся ошибки: по положению — нулевая;

по скорости E = 4 % при Uз = 0,08 В/с.

. Составить структурную схему нескорректированной системы.

.1. Найти общий коэффициент передачи разомкнутой системы.

.2. Определить коэффициент передачи усилителя.

.3. Определить передаточные функции всех элементов системы.

.4. Определить передаточную функцию разомкнутой системы.

.5. Определить передаточные функции замкнутой системы.

. Исследовать устойчивость системы.

.1. По критерию Гурвица.

.2. По критерию Михайлова.

.3. По критерию Найквиста на основе ЛЧХ.

. Определить корректирующее устройство.

.1. Построить желаемые ЛЧХ системы.

.2. Определить передаточную функцию корректирующего устройства.

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

. Техническая реализация корректирующего устройства.

.1. Аппаратная реализация.

.2. Программная реализация.

Еще статьи по технике и технологиям

Система отображения информации
В качестве примера структурной схемы контроллера была выбрана структурная схема автоматизированной системы контроля (АСК) [7]. В таких системах необходимая информация об объекте контроля, формируется с помощью измерительных датчико .

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

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

А ВТОРЕФЕРАТЫ КАНДИДАТСКИХ, ДОКТОРСКИХ ДИССЕРТАЦИЙ

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

Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.

Министерство образования Российской Федерации

Российский химико-технологический университет им. Д. И. Менделеева

Министерство образования Российской Федерации

Российский химико-технологический университет им. Д.И. Менделеева

Утверждено Редакционным советом университета в качестве учебного пособия Москва 2001 УДК 66.01 –52(076) ББК 35 я73 Д44 Рецензенты:

Доктор технических наук, профессор, проректор по учебной работе Московского государственного университета инженерной экологии М.Г. Беренгартен Доктор технических наук, профессор, заведующий кафедрой кибернетики химико-технологических процессов Российского химико-технологического университета им. Д. И. Менделеева Л.С. Гордеев Д44 Динамические звенья. Временные характеристики. Учеб. пособие / А. В. Беспалов, Н. И. Харитонов, С. Е. Золотухин, Л. Н. Финякин, А. С. Садиленко, В. Н. Грунский;

РХТУ им. Д. И. Менделеева. М., 2001. 80 с.

ISBN 5-7237-0328- В пособии рассмотрены свойства типовых динамических звеньев применительно к курсу СУ ХТП, даны их дифференциальные уравнения и метод их решения. Показана плодотворность введения понятия передаточной функции. Даны временные характеристики типовых динамических звеньев. Рассмотрены также соединения типовых динамических звеньев. Представленные в пособии задачи могут быть использованы на практических занятиях по курсу «Системы управления химико технологическими процессами».

Пособие предназначено для студентов химико-технологических cпециальностей.

УДК 66.01 –52(076) ББК 35 я ISBN 5-7237-0328-5 © Российский химико-технологический университет им. Д. И. Менделеева, -3 Содержание 1. Типовые динамические звенья. 1.1. Классификация типовых динамических звеньев. 1.2. Решение дифференциальных уравнений операционным методом. 1.3. Временные характеристики. ЗАДАЧИ. 2. Соединения (комбинации) звеньев. 2.1. Типовые комбинации звеньев. ЗАДАЧИ. Приложение. -4 1. ТИПОВЫЕ ДИНАМИЧЕСКИЕ ЗВЕНЬЯ Элементы автоматической системы управления химико технологическими процессами (АСУ ХТП) можно представить в виде типовых динамических звеньев, а также их соединений (комбинаций).

Динамические звенья называют типовыми, если изменение проходящего через них сигнала можно описать алгебраическим или дифференциальным уравнением не выше второго порядка (как правило, это линейные неоднородные дифференциальные уравнения с постоянными коэффициентами). Ниже приведена классификация основных типовых динамических звеньев. Во всех случаях зависимая переменная (выходной, y, сигнал) обозначена независимая переменная (время) а вынуждающая функция (возмущение, входное воздействие) x.

1.1. Классификация типовых динамических звеньев Статическое звено нулевого порядка ( безынерционное, усилительное, пропорциональное):

y( )= Kx( ), где K – статический коэффициент усиления.

Статическое звено первого порядка (апериодическое, инерционное, устойчивое):

dy + y = Kx, T d где Т – постоянная времени.

-5 Статическое звено второго порядка (инерционное, устойчивое):

d2y dy a2 2 + a1 + a0 y = b0 x, d d b0 a a, T = 2, 2T = 1, K= обозначив a0 a0 a получим d2y dy + 2T + y = Kx, T d d или 1 d 2 y 2 dy + + y = Kx, 0 d 0 d 2 – коэффициент 0 – частота где демпфирования, свободных колебаний (собственная частота) недемпфированной системы.

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

0 1 – колебательное;

Если = 1, то демпфирование называют критическим.

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

dy = x, Ta d где Та – постоянная времени интегрирования.

-6 Реальное интегрирующее звено (нейтральное звено второго порядка):

d2y dy TTa 2 + Ta = x, d d Идеальное дифференцирующее звено dx y =Td, d где Тd – постоянная времени дифференцирования.

Реальное дифференцирующее звено dy dx + y = Td T d d Звено запаздывания y( )=x( зап), где зап – время транспортного (или чистого) запаздывания.

Неустойчивое звено первого порядка dy y = Kx T d 1.2. Решение дифференциальных уравнений операционным методом Решаются указанные дифференциальные уравнения широко распространенным методом преобразования Лапласа, позволяющим ввести понятие передаточной функции и получить частотные характеристики.

-7 Преобразование Лапласа – математический метод, позволяющий сравнительно просто решать линейные дифференциальные уравнения. В результате преобразования дифференциальное уравнение (оригинал) приобретает форму алгебраического уравнения (изображения), в котором в качестве независимой переменной вместо времени используется s.

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

Преобразование Лапласа для временной функции запишем в виде:

F (s) = L[ f ( )] = f ( ) e s d, F(s) – функция s, L – символ где комплексной переменной прямого преобразования Лапласа.

1.2.1. Свойства преобразования Лапласа • Линейность Если функции f1() и f2() являются оригиналами и их изображения есть соответственно L[f1()] и L[f2()] и если величины с1 и с2 не зависят от и s, то справедливо равенство:

L[c1 f1( ) + c2 f 2 ( )] = c1 L[ f1( )] + c2 L[ f 2 ( )].

-8 • Дифференцирование оригинала f() f() Если функция и её производная являются оригиналами, то справедливо равенство:

L[ f ‘ ( )] = s L[ f ( )] f (0), f (0) = lim f ( ).

где + В случае преобразования второй производной:

L[ f ‘ ‘ ( )] = s L[ f ( )] s f (0) f ‘ (0).

В случае преобразования производной порядка n:

L[ f (n) ( )] = s n L[ f ( )] s n 1 f (0). s f (n 1) (0) f (n 1) (0), d n f ( ) ( n) (0) = lim f где.

+0 d n Таким образом, дифференцированию оригиналов отвечает умножение s.

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

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

• Интегрирование оригинала f() является оригиналом и L[f()] – его изображение, Если функция то справедливо равенство:

L[ f ( )] L [ f (t ) dt ] =, s -9 т.е. интегрированию в области оригиналов соответствует деление изображения на s.

• Теорема запаздывания f() является оригиналом и L[f()] – его изображение, Если функция L[f(-c)], c0, то изображение смещённого оригинала где определяется равенством L [ f ( c)] = L [ f ( )] ecs 1.2.2. Передаточная функция W(s) Передаточной функцией звена называется отношение L[y()] изображения выходного сигнала к изображению входного сигнала L[x()] при нулевых начальных условиях:

L[ x ( )] Пример 1: Получить дифференциальное уравнение термоэлектрического преобразователя (термопары), связывающее ЭДС термопары с измеряемой температурой:

Решение: Термоэлектрический преобразователь (первичный измерительный преобразователь, часто называемый термопарой) широко применяется в системах диагностики и управления. Входной величиной t термоэлектрического преобразователя является температура измеряемой среды, а выходной величиной – термоэлектродвижущая сила Е.

Сделаем следующие допущения:

а) рабочий спай термопары обменивается теплом (за счёт конвекции) только со средой, температура которой измеряется;

количество тепла, которое передается к рабочему спаю термопары или от него за счёт -10 теплопроводности (по проводам термопары) и за счёт лучистого теплообмена пренебрежимо мало;

б) температура внутри рабочего спая термопары везде одинакова;

в) температура свободных концов термопары равна 0°C;

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

В равновесном состоянии, когда температура рабочего спая термопары Е равна температуре измеряемой среды, термоэлектродвижущая сила термопары определённа и неизменна. Предположим, что в какой-то момент t. Вследствие разности времени температура среды ступенчато возросла на температур измеряемой среды и рабочего спая термопары последний нагревается до температуры среды. Нагрев рабочего спая термопары происходит не мгновенно, а в течение некоторого времени, зависящего от тепловой ёмкости рабочего спая и сопротивления передаче тепла.

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

dQ = A (t )] d, (1.1) – коэффициент A – поверхность где теплоотдачи конвекцией;

t – температура – температура теплопередачи рабочего спая;

рабочего спая термопары;

dQ, полученного термопарой при её нагреве на d Количество тепла градусов, равно:

dQ = mc d, (1.2) где m – масса рабочего спая термопары;

с – удельная теплоемкость спая.

Приравнивая выражения (1.1) и (1.2), получим:

-11 mc d + = t (1.3) A d Е Зависимость термоэлектродвижущей силы термоэлектрического преобразователя от температуры рабочего спая считаем линейной и аппроксимируем статическим звеном нулевого порядка (усилительным), т.е.:

E=K, (1.4) где K – коэффициент пропорциональности.

Е Из уравнения (1.4) выражаем через и подставляем в уравнение (1.3):

mc dE +E=Kt A d или, обозначая mc/A через T, получаем окончательное выражение:

dE +E=Kt T (1.5) d Полученное дифференциальное уравнение показывает, что по динамическим свойствам термоэлектрический преобразователь соответствует статическому звену первого порядка.

Пример 2: Найти передаточную функцию термоэлектрического преобразователя (статического звена первого порядка), уравнение динамики которого получено в предыдущем примере.

Решение: Преобразуем дифференциальное уравнение (1.5) по Лапласу, воспользовавшись свойством линейности:

dE T L + L ( E ) = K L (t ).

d Затем, применяя теорему о дифференцировании оригинала, получим:

-12 T [ s L ( E ) E ( 0 )] + L ( E ) = K L ( t ).

Учитывая, что по определению передаточной функции начальные условия являются нулевыми, т.е. E(0) = 0, получим алгебраическое уравнение:

L( E ) (Ts + 1) = K L (t ), откуда следует отношение Лапласовых изображений выходной и входной величин:

L(t ) Ts + 1.3. Временные характеристики О динамических свойствах звена (или системы) судят по реакции на типовые входные воздействия (единичное ступенчатое, единичное импульсное, с постоянной скоростью – рамповое и т. д.).

Единичное ступенчатое воздействие (функция Хевисайда).

x=1(). L[1()] = 1/s Аналитическое выражение Его изображение находим в таблице преобразований Лапласа Единичное импульсное воздействие (дельта-функция Дирака).

x = (). L[()] = Аналитическое выражение Его изображение находим в таблице преобразований Лапласа.

Единичное рамповое воздействие (воздействие с постоянной x = ·1().

скоростью). Аналитическое выражение Его изображение L[·1()] = 1/s2 находим в таблице преобразований Лапласа.

Изменения во времени выходной величины звена (системы) называют временными характеристиками звена (системы) или кривыми отклика.

-13 Переходная характеристика – аналитическое выражение отклика звена на единичное ступенчатое входное воздействие при нулевых начальных условиях:

s Произведение передаточной функции на изображение единичного ступенчатого воздействия, стоящее в квадратных скобках, соответствует L-1 – символ изображению переходной характеристики, а знак обратного преобразования Лапласа. Графическое изображение переходной функции – кривая разгона.

Импульсная переходная характеристика (весовая функция) – аналитическое выражение отклика звена на единичное импульсное входное воздействие при нулевых начальных условиях:

Графическое изображение импульсной переходной функции – кривая веса.

Рамповая переходная характеристика – аналитическое выражение отклика звена на единичное рамповое входное воздействие при нулевых начальных условиях:

y ( ) = L [ W ( s ) 2 ].

s Переходная, импульсная переходная и рамповая переходная характеристики являются частными случаями временных характеристик.

-14 ЗАДАЧИ 1.1. Получите дифференциальное уравнение стеклянного ртутного термометра, если известно, что объём находящейся в нём ртути 10-7 м3, площадь поверхности теплопередачи 10-4 м2, удельная теплоёмкость ртути 140 Дж/(кгК), плотность ртути 13600 кг/м3, а коэффициент теплопередачи от воздуха к ртути в условиях работы термометра равен 16 Вт/(м2К) (теплоёмкостью стекла можно пренебречь).

1.2. Дано дифференциальное уравнение дифманометра:

d2y dy 10 2 + 5 + 3 y = 5x, d d где х – измеряемая разность давления, а у – показания дифманометра.

Получите передаточную функцию дифманометра. Какому типовому динамическому звену соответствует дифманометр по своим динамическим свойствам?

1.3. Постоянная времени манометрического термометра при измерении температуры газового потока равна 6 мин.

Напишите дифференциальное уравнение манометрического термометра. Пользуясь преобразованием Лапласа, найдите динамическую погрешность манометрического термометра через 6 мин после ступенчатого изменения температуры газового потока от 20°C до 100°С.

Полагая, что по своим динамическим свойствам манометрический термометр соответствует статическому звену первого порядка, запишем его дифференциальное уравнение в виде:

-15 d + = t, d где – показания термометра, а t – измеряемая температура.

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

= t = 20°C Чтобы перейти к нулевым начальным условиям, примем за начало отсчёта температуру 20°С, введем новые переменные в виде отклонений параметров от их начальных значений:

t = t = и перепишем дифференциальное уравнение манометрического термометра, используя введённые переменные:

d ( 20 + ) + ( 20 + ) = ( 20 + t ).

d После простейших преобразований получим приведённое к нулевым начальным условиям уравнение манометрического термометра в отклонениях параметров от начальных значений:

d Передаточная функция манометрического термометра:

L ( t ) 6 s + Для ступенчатого изменения температуры на t:

-16 1 t L[ ( )] = W ( s ) L[ t 1( )] =, 6 s +1 s откуда, выполнив обратное преобразование Лапласа, найдем изменение показаний термометра во времени:

1 t ( ) = L = t (1 e 6 ).

6 s +1 s Через 6 мин после ступенчатого изменения температуры от 20°С до 100°С изменение показаний термометра будет следующим:

( ) = (100 20) (1 e Показание манометрического термометра через 6 мин:

Динамическая погрешность («недоход») манометрического термометра:

t = 70,6 100 = 29,4°C.

Определите динамическую погрешность манометрического термометра = 2Т, 3Т, 4Т, 5Т.

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

При сжигании 1900 м3/ч природного газа в печи устанавливается 1.4.

температура 1300°С. Найти статической коэффициент усиления печи, если при увеличении расхода газа до 2400 м3/ч температура в печи повышается до 1400°С.

1.5. Дана передаточная функция динамического звена:

W (s) = 5s + -17 Назовите звено и получите дифференциальное уравнение. Найдите переходную функцию этого звена и нарисуйте кривую разгона. Найдите импульсную переходную функцию этого звена и нарисуйте кривую веса.

Найдите рамповую переходную функцию этого звена, т.е. отклик звена на х = ·1() входное возмущающее воздействие и нарисуйте соответствующую кривую отклика.

1.6. Функция веса статического звена первого порядка изображена на рис. 1.1. Найдите параметры передаточной функции.

Рис. 1. 1.7. На вход статического звена первого порядка с постоянной времени Т = 15 с и коэффициентом усиления К = 2 подали единичное ступенчатое воздействие при нулевых начальных условиях. Чему равна выходная величина звена через 15 с, 30 с, 45 с, 60 с ?

-18 1.8. Покажите, что кривая разгона статического звена первого порядка обладает следующими свойствами:

а) разница между новым установившимся значением выходной величины y и её начальным значением yo равна К;

б) в точке = 0 скорость изменения выходной величины у равна К/Т;

в) касательная к кривой разгона при = 0 пересекает асимптоту в точке = Т;

г) величина подкасательной (проекция касательной на асимптоту) равна Т и не зависит от точки проведения касательной.

д) за время = Т выходная величина звена достигает 63% от установившегося значения.

В инженерных расчетах считается, что выходная величина процесса достигает своего установившегося состояния через = 4,6Т.

Какова величина фактического недохода при = 4,6Т?

1.9. Резервуар, постоянная времени которого равна 1 мин, а жидкость вытекает через отверстие с гидравлическим сопротивлением 1/90 м/(л/мин), работает в статическом режиме при расходе жидкости 10 л/мин. В момент = времени подача жидкости в резервуар увеличена до 100 л/мин, поддерживалась такой в течение 0,1 мин, после чего снова уменьшена до 10 л/мин. За счёт этого в резервуар поступило дополнительно 9 литров жидкости. Представьте графически изменение уровня жидкости в баке и сравните с импульсной переходной характеристикой бака, которая была бы получена при мгновенном добавлении в бак 9 литров воды.

-19 1.10. Стеклянный термометр с постоянной времени 2 мин, показывает 36°С через 2 мин после ступенчатого изменения температуры газового потока от 20°С.

Сколько времени потребуется для достижения того же показания, если вместо стеклянного термометра использовать манометрический термометр с постоянной времени 6 мин? Сравните динамические погрешности термометров через 2 мин после изменения температуры газового потока.

1.11. Для измерения температуры используется стеклянный ртутный термометр, у которого резервуар для ртути имеет очень тонкие стенки. В одном из экспериментов термометр нагрели до постоянной температуры в муфельной печи, после чего быстро перенесли в воздух при комнатной температуре и отмечали показания термометра во времени. Полученные результаты приведены в табл. 1. Таблица 1. Время, с 2 4 6 8 10 15 20 25 Показания 101,9 87,0 74,9 64,9 56,8 42,3 33,5 28,2 20, термометра, °С В другом эксперименте использовали тот же самый термометр для наблюдения за нагревом воздуха. Через некоторое время после начала измерения термометр показывал линейное увеличение температуры со скоростью 0,5°С/с. Определить, какова на самом деле температура воздуха, если термометр показывает 90°С.

1.12. Магнитный газоанализатор для определения содержания кислорода в газовых смесях по динамическим свойствам соответствует статическому звену первого порядка с постоянной времени, равной 10 с, и коэффициентом усиления, равным 1.

-20 Напишите дифференциальное уравнение магнитного газоанализатора.

Получите, пользуясь преобразованием Лапласа, передаточную функцию магнитного газоанализатора. Определите показания магнитного газоанализатора через 10, 20 и 40 с после ступенчатого изменения концентрации кислорода от 15 до 21%.

1.13. Первичный измерительный преобразователь магнитного газоанализатора (см. задачу 1.12) соединен с реактором трубкой для отбора проб анализируемой газовой смеси, из-за чего появляется транспортное запаздывание, равное 5 с.

Какими будут показания магнитного газоанализатора через 10, 20 и 40 с после такого же ступенчатого изменения концентрации кислорода, как в задаче 1.12? Рассчитайте длину трубки для отбора проб, если известны её диаметр (0.01 м) и расход газового потока (60 л/ч).

1.14. Термокондуктометрический газоанализатор ИВА–1В11 предназначен для измерения концентрации водорода в контуре охлаждения турбогенераторов. Благодаря некоторым оригинальным техническим решениям газоанализатор обладает уникальным быстродействием.

Предполагая, что по своим динамическим свойствам газоанализатор ИВА–1В11 соответствует статическому звену первого порядка, определите его постоянную времени, если время установления показаний газоанализатора 0,9 = 15 с (время установления показаний 0,9 определяют как время, в течение которого изменение показаний газоанализатора достигнет 90% полного изменения показаний, вызванного ступенчатым изменением концентрации водорода).

1.15. Дана передаточная функция трубопровода:

-21 W ( s ) = e 5 s Назовите звено. Получите переходную функцию трубопровода и нарисуйте кривую разгона для него.

1.16. По динамическим свойствам газоанализатор соответствует статическому звену первого порядка с постоянной времени T = 15 с, и коэффициентом усиления K = 1. Для отбора проб анализируемой газовой смеси первичный измерительный преобразователь газоанализатора соединен с реактором трубкой длиной 10 м и диаметром 0,006 м. Расход газа, отбираемого из реактора, равен 100 л/ч.

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

1.17. Установка состоит из подогревателя и потенциометра КСП3 с регулятором, поддерживающим температуру в подогревателе, равную 120°С.

В подогреватель поместили два теплоприемника: терморезистор термометра сопротивления и термобаллон манометрического термометра. После того, как температура теплоприемников достигла 120°С и некоторое время оставалась постоянной, теплоприемники быстро вынули из подогревателя и поместили в среду с комнатной температурой, тем самым имитируя скачкообразное понижение измеряемой температуры. Полученные в результате такого эксперимента кривые разгона термометров показаны на Т рис. 1.2. По кривым разгона определите значение постоянной времени каждого термометра четырьмя способами и сравните полученные результаты (статической погрешностью термометров пренебречь).

-22 Рис. 1. Способ 1. Провести касательную к кривой разгона в точке, 1 (часто соответствующей некоторому произвольному моменту времени выбирают 1 = 0);

пересечение касательной с асимптотой даст второй момент времени 2. Постоянная времени Т = 2 — 1.

Способ 2. Исходя из соотношения ( ) ( ) = e, ( +T ) ( ) () – показание, () – конечное где термометра в момент времени установившееся показание термометра.

Способ 3. По котангенсу угла наклона прямой линии:

( ) ( ) =, ln ( ) ( 0) T где (0) – показание термометра в момент времени = 0.

Способ 4. По соотношению -23 (T ) (0) = 0,63 [ () (0)], т.е. за время, равное постоянной времени, изменение температуры достигает 63% полного изменения.

1.18. Температуру в печи измеряют термоэлектрическим термометром, в комплект которого входят:

— первичный хромель-копелевый термоэлектрический преобразователь, допускаемое отклонение термо-ЭДС которого от градуировочного значения составляет ± 0,2 мВ в измеряемом интервале температур;

— милливольтметр класса точности 1 с диапазоном измерения от 0°С до 600°С. Внутреннее сопротивление милливольтметра равно 100 Ом и он отградуирован для внешнего сопротивления (термоэлектрический преобразователь + соединительные провода) Rвн = 15 Ом. Действительная температура в печи 400°С.

а) Определите предельное значение абсолютной основной погрешности измерительного комплекта при измерении данной температуры.

б) Определите абсолютную дополнительную погрешность, вызванную отличием внешнего сопротивления милливольтметра от градуировочного значения Rвн, если известно, что сопротивление термоэлектрического преобразователя 0,5 Ом и сопротивление соединительной линии 4,5 Ом в) Какова чувствительность хромель-копелевого термоэлектрического преобразователя вблизи измеряемой температуры, если его номинальную статическую характеристику можно описать полиномом вида 2 E (t ) = c + c t + c t +. + c t, 0 1 2 E – ЭДС t – температура где термоэлектрического преобразователя, мкВ;

с1, с2, …, с9 – рабочего спая, °C (температура «холодного спая» 0°C);

коэффициенты полинома, значения которых приведены в табл. 1.2?

-24 Таблица 1. с0 с 0 1.5425922111E- с1 с 5.8695857799E1 -2.4850089136E- с2 с 4.3110945462E-2 2.3389721459E- с3 с 5.7220358202E-5 -1.1946296815E- с4 с -5.4020668085E-7 2.5561127497E- г) Какую ЭДС развивает хромель-копелевый термоэлектрический преобразователь при температуре в печи 400°С, если температура его свободных концов равна 20°С?

д) Чему равна динамическая погрешность термоэлектрического термометра через 3 мин после того, как хромель-копелевый термоэлектрический преобразователь, находившийся в производственном помещении при температуре 20°С, быстро поместили в печь, где температура 400 С? По своим динамическим свойствам хромель-копелевый термоэлектрический преобразователь в условиях измерений соответствует Т = 3 мин, статическому звену первого порядка с постоянной времени а милливольтметр можно считать звеном нулевого порядка. Чему равнялась бы в тех же условиях динамическая погрешность манометрического термометра с постоянной времени Т = 6 мин?

е) Какова динамическая погрешность термоэлектрического термометра через 3 мин после пуска печи, если известно, что до пуска температура в печи была 20°С, а после пуска температура в печи поднимается с постоянной скоростью, равной 30°С/мин?

Для того, чтобы ответить на поставленные вопросы, следует вспомнить следующие термины и определения, используемые в метрологии:

-25 Абсолютной погрешностью измерительного прибора называют A разность между показанием прибора и истинным значением измеряемой величины Aи:

= A A и Когда истинное значение измеряемой величины неизвестно, вместо него Ад, полученное используют действительное значение измеряемой величины с помощью образцового прибора:

= A A д Абсолютная погрешность имеет наименование и знак.

Относительной погрешностью называют отношение абсолютной погрешности к истинному значению измеряемой величины, или, если оно неизвестно, к действительному значению:

Aи Aд Приведенной погрешностью называют отношение абсолютной погрешности к нормирующему значению Анорм:

прив. = 100% Aнорм Аmin Для приборов, нижний предел диапазона измерения которых выше нуля, нормирующее значение принимается равным верхнему пределу Аmax, измерения если же нижний предел диапазона измерения равен нулю или ниже нуля, то нормирующее значение принимается равным диапазону измерения Аmax — Аmin.

-26 Класс точности – обобщённая характеристика средства измерений, определяемая пределами допускаемых основной и дополнительных погрешностей, а также другими свойствами средства измерений, влияющих на точность. Основная погрешность средства измерений – погрешность, получаемая при измерении в нормальных условиях, установленных стандартами на прибор;

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

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

max К. Т. = 100%, Amax Аmin где max – допускаемая основная абсолютная погрешность.

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

Динамические погрешности связаны с изменением измеряемой величины и обусловлены запаздыванием изменений показаний измерительного прибора.

1.19. Концентрация продукта реакции на выходе из реактора с мешалкой (c, моль/м3) зависит от расхода подаваемого в реактор реагента (F, кг/мин) в соответствии с уравнением:

-27 dc ( ) + c ( ) = 4 F ( 2) d Определить, как будет меняться концентрация продукта после того, как произошло ступенчатое изменение расхода реагента от 2 кг/мин до 2,4 кг/мин, если перед этим реактор находился в статическом режиме.

Решить задачу двумя способами:

— применить преобразование Лапласа к исходному уравнению, учитывая ненулевые начальные условия;

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

d ( z ) d ( y ) + z = 2 y + 5, d d y – входное где воздействие для регулятора (соответствует отклонению z – изменение текущего значения регулируемого параметра от заданного);

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

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

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

Решение: В пространство над мембраной 3 поступает воздух под Р, давлением создавая усилие на мембрану. Под действием этого усилия -28 мембрана прогибается, перемещая вниз шток 1 с затвором 5. Через клапан проходит жидкость в количестве, зависящем от степени открытия клапана.

h Перемещению штока клапана с затвором препятствует сила упругости пружины 2, трение штока в сальнике 4 и затвора 5 в протекающей жидкости.

Рис. 1. Рассмотрим в качестве входной величины пневматического мембранного исполнительного механизма с регулирующим органом Р, а выходной – перемещение h штока 1 с затвором изменение давления 5 (отсчет ведем в малых приращениях от равновесного состояния).

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

F +F +F = F, и тр пр F – входное где усилие (равно произведению давления воздуха на площадь Fи – сила мембраны);

инерции (равна произведению массы подвижной Fтр – сила системы на ускорение);

вязкого трения (пропорциональна Fпр – сила скорости перемещения подвижной системы);

противодействия пружины (пропорциональна её сжатию).

-29 Подставив соответствующие выражения сил в уравнение равновесия, получим дифференциальное уравнение движения штока с затвором в следующем виде:

d 2 ( h ) d ( h ) +r + c h = A P, m d d m – масса с– где движущихся частей с учетом половины массы пружины;

r – коэффициент A – площадь коэффициент жёсткости пружины;

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

m d 2 ( h ) r d ( h ) A + + h = P c d c d 2 c m r A = 2 T, = T 2, = K, Обозначив получим уравнение c c c динамики пневматического мембранного исполнительного механизма с регулирующим органом с учетом силы вязкого трения и инерции его подвижной системы:

d 2 ( h ) d ( h ) + 2 T + h = K P T d d Для нахождения передаточной функции звена, воспользуемся свойствами преобразования Лапласа (см. приложение, табл. 1, свойства 2–4):

T s L( h) + 2 T s L( h) + L ( h) = K L ( P), откуда следует отношение изображений по Лапласу выходной и входной величин -30 L( h ) K W ( s) = =2.

L ( P ) T s + 2 T s + Пример 4: Получить уравнение динамики U-образного жидкостного h дифманометра (рис. 1.4), связывающее показания дифманометра c измеряемой разностью давлений P = P1 – P2.

Решение: Запишем уравнение равновесия сил, для упрощения пренебрегая плотностью паров над уровнем жидкости и предположив, что движение жидкости в трубке поршневое, а падение давления, вызванное трением, пропорционально Рис. 1.4 скорости потока:

F +F + F = F, и тр т 1 d 2h F = LA А– где – сила инерции ( – плотность жидкости, 2 d и площадь поперечного сечения трубки);

32 Lµ 1 dh = P A = 2 A – сила F трения (Pтр – потери 2 d d тр тр µ– давления на трение, определяемые по уравнению Хагена-Пуазейля, вязкость жидкости, d – внутренний диаметр трубки дифманометра);

F = hAg – сила тяжести, (g – ускорение свободного падения);

т -31 F = A ( P P2 ) – сила, приложенная к жидкости со стороны измеряемого давления.

Подставляя эти выражения в уравнение равновесия сил, получим LA d 2 h 16 µLA dh + ghA = A P.

+ d 2 d 2 d gA, Разделив все члены уравнения на получим уравнение динамики U-образного жидкостного дифманометра L d 2 h 16µL dh + +h= P, gd d g 2 g d 2 которое обычно записывают в следующем стандартном виде:

1 d 2 h 2 dh + + h = K P, 0 d 0 d 2 0 – собственная где частота недемпфированной системы, рад/с;

— коэффициент 0 демпфирования. Физический смысл и становится понятным после анализа решения уравнения при единичном ступенчатом изменении входного давления. При коэффициенте демпфирования выходной сигнал будет совершать затухающие колебания (рис. 1.5). Такую систему называют слабо демпфированной. При коэффициенте = демпфирования (недемпфированная система) выходной сигнал представляет собой незатухающие колебания с частотой 0 и амплитудой K.

= 1, Если коэффициент демпфирования то имеет место критическое демпфирование и жидкость в дифманометре приходит к состоянию равновесия без перерегулирования. При коэффициенте демпфирования изменение уровня жидкости в дифманометре происходит по S-образной кривой (без колебательного процесса) тем медленнее, чем больше -32 коэффициент демпфирования. Такие системы называют сильно демпфированными.

Рис. 1. 1.21. Получите переходную функцию дифманометра (см. задачу 1.2), рассчитайте коэффициент демпфирования и нарисуйте кривую разгона в общем виде.

1.22. Зависимость между концентрацией продукта реакции на выходе из реактора (с, моль/м3) и расходом реагента, поступающего в реактор (F, м3/ч), описывается уравнением:

d 2c dc 4 2 +6 + c = 5F d d Получите передаточную функцию реактора. Какому типовому динамическому звену соответствует реактор по своим динамическим -33 свойствам? Рассчитайте коэффициент демпфирования и нарисуйте кривую разгона реактора в общем виде.

1.23. Дано дифференциальное уравнение динамического звена:

d3y d2y dy dx 5 3 + 10 2 + 3 = d d d d Получите передаточную функцию. Рассчитайте коэффициент демпфирования и нарисуйте кривую разгона в общем виде.

1.24. Экспериментально получена кривая разгона объекта регулирования (рис. 1.6). Найдите передаточную функцию объекта и запишите его дифференциальное уравнение.

Рис. 1. Порядок решения может быть следующим (рис. 1.7):

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

• определяем новое установившееся значение выходной величины -34 y = lim y ( ), что позволяет найти статический коэффициент усиления K;

• проводим касательную в точке перегиба i;

• определяем соотношение если оно больше 0,74, то относим объект к статическим звеньям второго порядка (в противном случае считаем, что порядок объекта выше второго и пользуемся подсказками к задаче 1.34);

TCD/TBD • по соотношению и по графику (рис. 1.8) находим коэффициент демпфирования ;

-35 • определяем постоянную времени T, исходя из соотношения:

• для полученного значения коэффициента демпфирования по TAB/TBD, графику (рис. 1.8) находим соотношение используя его, рассчитываем TAB и определяем время транспортного запаздывания:

Если время запаздывания оказывается очень малым или отрицательным, можно попробовать скорректировать положение касательной для того, чтобы точно получить зап = 0.

Рис. 1. -36 1.25. Объект состоит из ёмкости и насоса с постоянной производительностью, установленного на выходе из ёмкости (рис. 1.9). Связь Fвх между изменением подачи жидкости в бак и изменением уровня жидкости в ёмкости L описывается передаточной функцией:

5s Рис. 1. Какому типовому динамическому звену соответствует объект?

Получите переходную функцию объекта и нарисуйте кривую разгона.

Найдите импульсную переходную функцию объекта и нарисуйте кривую веса. Найдите отклик объекта на входное возмущающее воздействие Fвх = 2·1() и нарисуйте соответствующую кривую отклика.

1.26. Дано дифференциальное уравнение тахогенератора:

d U вых = 5, d Uвых – напряжение, – угол где снимаемое с клемм тахогенератора;

поворота вала тахогенератора.

Какому типовому динамическому звену соответствует тахогенератор?

Получите передаточную функцию тахогенератора. Определите переходную функцию тахогенератора и нарисуйте кривую разгона. Найдите отклик -37 = 3·1() тахогенератора на входное воздействие и нарисуйте соответствующую кривую отклика.

1.27. Дана передаточная функция динамического звена:

3s W (s) = 3s + Назовите звено и напишите его дифференциальное уравнение. Найдите переходную функцию звена и нарисуйте соответствующую кривую разгона.

Найдите импульсную переходную функцию звена и постройте кривую веса.

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

1.28. Один из модулей электрического регулятора описывается уравнением:

dU 2 dU + U2 = 3 d d Какому типовому динамическому звену соответствует указанный модуль регулятора? Получите его передаточную функцию. Постройте график U изменения во времени выходного напряжения после ступенчатого U1 = 2·1().

изменения напряжения на входе модуля Найдите изменение U2, выходного напряжения если напряжение на входе модуля меняется с постоянной скоростью U1 = 3·1().

1.29. Дана передаточная функция динамического звена:

W (s) = 2s -38 Напишите дифференциальное уравнение звена. Найдите переходную функцию звена и нарисуйте соответствующую кривую разгона. Найдите весовую функцию звена и изобразите её графически.

1.30. Дана передаточная функция динамического звена:

3s W (s) = 2s Получите дифференциальное уравнение звена, найдите его переходную функцию и постройте кривую разгона. По виду полученной кривой разгона определите, что это за звено.

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

б) h( ) = 5(1 + e10 ) 1( ) ;

в) h( ) = 0,1 1( ) г) h( ) = 2 [1 (1 + 0,5 ) e0,5 ]1( ) д) Найдите параметры этих звеньев и запишите выражения передаточных функций.

1.32. Определите, какие динамические звенья имеют импульсную переходную функцию (функцию веса) следующего вида:

а) -39 w( ) = 8 e0,5 1( ) ;

б) w( ) = 50(e 5 e 10 ) 1( ) ;

в) w( ) = 0,11( ) г) w( ) = 2 ( ) e0,5 1( ) д) Найдите параметры этих звеньев и запишите выражения передаточных функций.

1.33. Определите, какие динамические звенья имеют рамповую переходную характеристику следующего вида:

б) y( ) = 50 (e5 + 5 1) 1( ) ;

в) y( ) = 0,11( ) г) y( ) = 2 2 1( ) д) Найдите параметры этих звеньев и запишите выражения передаточных функций.

1.34. Экспериментально получены кривые разгона четырёх объектов (рис. 1.10). Обработайте их, как рекомендовано в задаче 1.24 (рис. 1.7) и, используя данные табл. 3, определите передаточные функции объектов в виде:

e зап s K W (s) = (Ts + 1)n -40 Рис. 1. Подсказки:

• определяем новое установившееся значение выходной величины y = lim y ( ), что позволяет найти статический коэффициент усиления K;

• проводим касательную в точке перегиба i;

yi/y, • по графику определяем округляем его до ближайшего меньшего табличного значения и находим порядок объекта n;

TBD • по графику определяем и, используя табличное значение TBD/T для найденного порядка n, вычисляем постоянную времени Т;

ТВD • по найденному из графика значению и табличному значению ТАВ/ТВD вычисляем ТАВ;

• рассчитываем время запаздывания зап:

-41 зап = TB TAB Если время запаздывания оказывается очень малым или отрицательным, можно попробовать скорректировать положение касательной для того, чтобы точно получить зап = Таблица n yi/y ТАВ/ТВD ТВD/T ТAB/T ТAC/T 1 0 0 1 0 2 0,264 0,104 2,718 0,282 3 0,323 0,218 3,695 0,805 4 0,353 0,319 4,463 1,425 5 0,371 0,410 5,119 2,100 1.35. В изотермическом реакторе идеального вытеснения протекает без изменения объёма химическая реакция первого порядка. Определить изменение концентрации во времени на выходе из реактора:

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

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

1.36. На выходе химического реактора установлен металлический трубопровод с тонкими стенками длиной 10 м и внутренним диаметром 5 см.

Разбавленный водный раствор протекает через него со скоростью 10 см/с.

Для измерения температуры жидкости в конце трубопровода установлен термометр, по своим динамическим свойствам соответствующий статическому звену первого порядка с постоянной времени 3 с. Температура окружающей среды равна 20°С, а коэффициент теплопередачи между -42 жидкостью, протекающей по трубопроводу, и окружающей средой равен 50 Вт м-2 К-1. Можно предположить, что движение жидкости в трубопроводе поршневое, и пренебречь продольной теплопроводностью.

• Записать уравнение теплопередачи в тонком слое трубопровода.

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

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

• Определить изменение показаний термометра и представить графически это изменение в двух указанных ниже случаях:

а) температура жидкости на выходе из реактора ступенчато изменяется от 70°С до 80°С;

б) температура жидкости на выходе из реактора меняется синусоидально между 70°С и 80°С с периодом 18 с.

-43 2. СОЕДИНЕНИЯ (КОМБИНАЦИИ) ЗВЕНЬЕВ Структурной схемой в теории автоматического управления называют графическое изображение математической модели автоматической системы управления в виде соединений звеньев. Звено на структурной схеме условно обозначают в виде прямоугольника с указанием входных и выходных величин, а также передаточной функции внутри него.

Сумматоры и элементы сравнения изображают в виде круга, разделенного на секторы (рис. 2.1). В элементе сравнения сектор, на который подается «вычитаемое», затемняют.

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

2.1. ТИПОВЫЕ КОМБИНАЦИИ ЗВЕНЬЕВ 2.1.1. Последовательное соединение При последовательном соединении звеньев выходная величина каждого предшествующего звена является входным воздействием последующего звена. При преобразовании структурных схем цепочку из последовательно -44 соединенных звеньев можно заменить одним звеном с передаточной функцией, равной произведению передаточных функций отдельных звеньев:

L( y ) W ( s) = = W ( s ) W ( s ). W ( s ), L( x) 1 2 n L(x) – изображение L(y) – где по Лапласу входного воздействия;

W1(s), W2(s),…Wn(s) – изображение по Лапласу выходного сигнала;

передаточные функции отдельных звеньев.

Структурная схема последовательного соединения звеньев дана на рис. 2.2:

Рис. 2. 2.1.2. Параллельное соединение При параллельном соединении на вход всех звеньев подается один и тот же сигнал, а выходные величины складываются. Систему из нескольких параллельно соединенных звеньев можно заменить одним звеном с передаточной функцией, равной сумме передаточных функций входящих в нее звеньев:

L( y) W (s) = = W (s) + W (s) +. + W (s) 1 2 n L( x) Структурная схема параллельного соединения звеньев дана на рис. 2.3:

-45 Рис. 2. 2.1.3. Соединение с обратной связью Принято считать, что звено охвачено обратной связью, если его выходной сигнал через какое-либо другое звено подается на вход. Если сигнал обратной связи вычитается из входного воздействия, то обратную связь называют отрицательной. Если сигнал обратной связи складывается с входным воздействием, то обратную связь называют положительной.

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

L( y) W1(s) W ( s) = = L( x) 1 + W1(s) W2 (s) Структурная схема соединения звеньев с отрицательной обратной связью дана на рис. 2.4:

-46 Рис. 2. Передаточная функция замкнутой цепи с положительной обратной связью равна дроби, числитель которой – передаточная функция прямой цепи, знаменатель – единица минус произведение передаточных функций прямой и обратной цепи:

L( y) W1(s) W ( s) = = L( x) 1 W1(s) W2 (s) Структурная схема соединения звеньев с положительной обратной связью дана на рис. 2.5:

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

-47 ЗАДАЧИ 2.1. Структурная схема устойчивого объекта первого порядка имеет вид (рис. 2.6):

Рис. 2. Назовите звенья. Получите передаточную функцию объекта. Напишите дифференциальное уравнение объекта. Постройте кривую разгона объекта.

Т К 2.2. Как изменятся постоянная времени и коэффициент усиления статического звена первого порядка, если его охватить отрицательной обратной связью с передаточной функцией цепи обратной связи Woc(s) = Td s.

2.3. Дана комбинация звеньев (рис. 2.7). Получите передаточную функцию комбинации. Какому звену эквивалентна комбинация? Нарисуйте кривую разгона полученного звена.

Рис. 2. -48 2.4. Дана структурная схема объекта регулирования (рис. 2.8):

Рис. 2. Определите его передаточную функцию. Найдите переходную характеристику объекта и постройте кривую разгона, обратив особое внимание на её поведение при малых значениях.

2.5. Дана структурная схема объекта регулирования (рис. 2.9):

Рис. 2. Определите его передаточную функцию. Найдите переходную характеристику объекта и постройте кривую разгона, обратив особое внимание на её поведение при малых значениях.

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

-49 2.7. Найдите передаточные функции двух таких статических звеньев первого порядка, при параллельном соединении которых передаточная функция полученной системы была бы тождественна передаточной функции статического звена второго порядка:

K W (s) = (T s +1) (T s +1) 1 Т, коэффициент демпфирования 2.8. Как изменятся постоянная времени К, а также время и форма переходного и статический коэффициент усиления процесса статического звена второго порядка, если его охватить отрицательной обратной связью, передаточная функция которой Woc(s) = Кос ?

2.9. Объект представляет собой последовательное соединение статического Т = 5 с, коэффициент усиления звена первого порядка (постоянная времени К = 2) и звена запаздывания (зап = 2 с). Получите передаточную функцию объекта и постройте кривую разгона.

2.10. Объект управления можно представить следующей структурной схемой (рис. 2.10).

Рис. 2. -50 Назовите звенья. Получите передаточную функцию комбинации звеньев. Напишите дифференциальное уравнение объекта. Постройте кривую разгона объекта в общем виде.

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

W (s) = Ta s Постройте кривые разгона и веса для одного идеального интегрирующего звена и для системы.

2.12. Объектом, обладающим запаздыванием, является ленточный транспортёр 2 для транспортировки сыпучего материала из бункера 1 в приёмный бункер 3 (рис. 2.11). Входной величиной ленточного транспортёра l = 30 м и скоростью движения ленты v =6 м/мин является расход длиной сыпучего материала из бункера 1 (F1, м3/мин), а выходной величиной – подача сыпучего материала в приёмный бункер 3 (F2, м3/мин). Зависимость между уровнем в приёмном бункере (L, м) и подачей в него сыпучего W(s) = 1/Tas, материала описывается передаточной функцией где Тa = 5 мин.

Как будет меняться уровень сыпучего материала в приемном бункере 3, если расход сыпучего материала при его подаче из бункера 1 на ленточный транспортер 2 скачком изменился от 2 м3/мин до 4 м3/мин?

-51 Рис. 2. 2.13. Дана структурная схема объекта управления (рис. 2.12):

Рис. 2. Назовите звенья. Получите передаточную функцию объекта. Напишите его дифференциальное уравнение. Постройте кривую разгона объекта.

2.14. Нейтральный объект представлен в виде комбинации звеньев (рис. 2.13). Назовите звенья. Получите передаточную функцию объекта.

Напишите его дифференциальное уравнение. Постройте кривую разгона объекта.

Рис. 2. -52 2.15. Объект состоит из двух цилиндрических аппаратов, соединенных трубопроводом, на котором установлен вентиль (рис. 2.14).

Рис. 2. R Сопротивление потоку жидкости из первого аппарата во второй значительно. Если жидкость из объекта отводится самотеком через R2, гидравлическое сопротивление то объект устойчив, а его переходная F характеристика, полученная при скачкообразном изменении притока в первый аппарат от 1 м3/мин до 2 м3/мин, показана на рис. 2.15.

Определите передаточную функцию объекта, связывающую изменение L2 c изменением притока жидкости в уровня жидкости во втором аппарате первый аппарат F (см. задачи 1.24 и 1.34).

Определите передаточную функцию каждого аппарата, считая входной F, величиной первого аппарата изменение притока выходной величиной первого аппарата и входной величиной второго – F1, и выходной величиной второго аппарата – изменение уровня во втором аппарате L2.

Как будет меняться уровень жидкости в первом аппарате, если приток F в первый аппарат мгновенно изменился от 1 м3/мин до 2 м3/мин?

Используя переходную характеристику объекта, выберите закон управления.

-53 Рис. 2. 2.16. В объекте, свойства которого описаны в задаче 2.15, аппараты соединили так, как показано на рис. 2.16.

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

Рис. 2. -54 Определите, как будут изменяться во времени уровни жидкости в F в первый аппарат мгновенно изменился от каждом аппарате, если приток 1.0 до 2.0 м3/мин.

2.17. Для объекта, состоящего из двух аппаратов и насоса с постоянной производительностью (рис. 2.17), получены кривые отклика на мгновенное F в первый аппарат на 1 м3/мин, приведенные на изменение притока L2 – рис. 2.18 (L1 – изменение уровня жидкости в первом аппарате, изменение уровня жидкости во втором аппарате).

Рис. 2. Напишите дифференциальное уравнение объекта, пользуясь кривыми разгона (рис. 2.18), определите коэффициенты дифференциального уравнения и представьте структурную схему объекта.

Как будет изменяться уровень жидкости в первом и во втором аппаратах, если приток жидкости в первый аппарат мгновенно изменили на 2 м3/мин, поддерживали его постоянным в течение 10 мин, а затем вернули к прежнему значению?

-55 Рис. 2. 2.18. В объекте, свойства которого описаны в задаче 2.17, аппараты соединили так, как показано на рис. 2.19.

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

-56 Определите, как будут меняться во времени уровни жидкости в каждом F аппарате при мгновенном изменении притока в первый аппарат на 1 м3/мин.

2.19. Реальный химический реактор можно представить в виде соединения V реактора идеального вытеснения объёмом с реактором идеального V2, помещённым в линию рецикла так, как это показано смешения объёмом на рис. 2.20:

Рис. 2. Получите передаточную функцию такого соединения реакторов, у() связывающую изменение концентрации реагента на выходе с x() изменением его концентрации на входе при условии, что расходы реакционной смеси F и F1 (см. рис. 2.20) не изменяются, объём реакционной массы в реакторе идеального смешения постоянен, а химическая реакция не протекает.

2.20. Определите передаточную функцию объекта регулирования, структурная схема которого дана на рис. 2.21.

-57 Рис. 2. 2.21. В каскаде из трех реакторов идеального смешения (рис. 2.22) проводят жидкофазную реакцию. Объем реакционной смеси не изменяется.

Отклонение концентрации реагента на входе в реактор от номинального сi-1 F значения или отклонение расхода реакционной смеси от номинального значения вызывает изменение концентрации на выходе из сi реактора в соответствии со структурной схемой, приведенной на рис. 2.22.

с Определите изменение концентрации на выходе из третьего реактора:

а) если произошло единичное ступенчатое изменение расхода F = 1(), а концентрация реагента на входе в первый реактор не изменилась (cо = 0);

б) если произошло единичное ступенчатое изменение концентрации с0 = 1(), а расход реакционной смеси не изменился (F = 0).

-58 Рис. 2. 2.22. Получите передаточную функцию системы, структурная схема которой изображена на рис. 2.23. Какому динамическому звену она эквивалентна?

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

Рис. 2. -59 2.23. Регулятор (звено 1 на рис. 2.24) охвачен обратной связью (звено 2).

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

Рис. 2. 2.24. В замкнутой системе управления (рис. 2.25) из-за ошибки при подключении управляющего устройства (звено 2) к реактору идеального смешения (звено 1) вместо отрицательной обратной связи появилась положительная обратная связь.

Рис. 2. Изменяя параметр настройки управляющего устройства К (0,5;

2,0) в цепи обратной связи, получите передаточные функции системы.

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

-60 2.25. Идеальное дифференцирующее звено (его постоянная времени Td = 2) дифференцирования охвачено стопроцентной отрицательной обратной связью.

Нарисуйте блок-схему и получите передаточную функцию комбинации. Какому динамическому звену она соответствует? Как изменился отклик регулятора на рамповое входное воздействие х = 5 ·1() после добавления обратной связи?

2.26. В цепь обратной отрицательной связи пневматического двухвходового усилителя 1 с большим коэффициентом усиления (К 1) включена пневматическая RC-цепочка, состоящая из пневматического сопротивления и пневматической ёмкости 2 (рис. 2.26).

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

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

Рис. 2. 2.27. Дана комбинация звеньев (рис. 2.27). Назовите звенья, стоящие в прямой и обратной цепи. Найдите точную и приближённую передаточную -61 функцию комбинации звеньев. Какому звену (точно и приблизительно) соответствует эта комбинация?

Рис. 2. 2.28. Структурная схема промышленного И-регулятора дана на рис. 2.28.

Передаточные функции звеньев, представленных на схеме:

W ( s ) = K = 10 ( K 1) W ( s) = Ts + Найдите передаточную функцию регулятора. Постройте графики переходных характеристик (кривые разгона) для значений постоянной времени интегрирования Ta = 2 мин и Ta = 10 мин.

Рис. 2. -62 2.29. На рис. 2.29 дана структурная схема промышленного регулятора.

Получите его передаточную функцию. Какой закон регулирования реализуется? Постройте кривую разгона регулятора.

Рис. 2. 2.30. Учитывая, что коэффициент усиления звена 1 величина большая, найдите передаточную функцию регулятора, блок-схема которого представлена на рис. 2.30. Определите закон регулирования и постройте кривую разгона.

Рис. 2. 2.31. Даны структурные схемы двух промышленных регуляторов (рис. 2.31 и рис. 2.32).

В первой структурной схеме (рис. 2.31):

W ( s) = 5s 1 Во второй структурной схеме (рис. 2.32):

5 s +1 2,5 s 1 2 Получите передаточные функции промышленных регуляторов. Какие законы управления реализуют регуляторы? Постройте кривые разгона промышленных регуляторов.

Рис. 2. Рис. 2. 2.32. ПИ-закон регулирования может быть реализован устройством, блок схема которого приведена на рис. 2.33. Определите передаточную функцию W2(s) звена, стоящего в цепи обратной связи.

-64 Рис. 2. 2.33. В регулирующем блоке пневматической системы «СТАРТ», реализующем ПИ-закон регулирования, имеется цепочка пневмоэлементов (рис. 2.34), которую можно представить структурной схемой, изображённой на рис. 2.35.

Чему в цепочке пневмоэлементов соответствуют обозначения на структурной схеме?

Получите передаточную функцию цепочки пневмоэлементов. Какому звену она соответствует? Напишите дифференциальное уравнение пневматической цепочки.

Рис. 2. -65 Рис. 2. 2.34. В пневматических регулирующих блоках системы «СТАРТ» заданное и текущее значения регулируемой величины вводятся в виде давления Рзад Рх, сжатого воздуха и а управляющее воздействие, формируемое регулирующим блоком (регулятором), передается в виде давления сжатого Рвых.

воздуха Структурная схема одного из регулирующих блоков системы «СТАРТ» приведена на рис. 2.36.

Рис. 2. Найдите передаточные функции регулирующего блока по каналу Wx(s) = L(Pвых)/L(Px) и по текущего значения регулируемой величины Wзад(s) = L(Pвых)/L(Pзад). Какие законы каналу заданного значения регулирования реализуются по каждому из каналов?

-66 Определите, как будет меняться давление на выходе регулирующего блока Рвых:

Pзад = 1() а) при единичном ступенчатом изменении задания и постоянном значении регулируемого параметра Px = const;

б) при единичном ступенчатом изменении регулируемого параметра Pх = 1() и постоянном значении задания Рзад = const.

2.35. Пневматический приборный регулятор может быть представлен структурной схемой, приведённой на рис. 2.37:

Рис. 2. где y – отклонение регулируемого параметра от заданного значения;

Pвых – управляющее воздействие, формируемое регулятором и передаваемое в виде давления сжатого воздуха на исполнительный механизм системы регулирования.

Получите передаточную функцию регулятора. Какой закон регулирования он реализует? Постройте кривую разгона регулятора.

-67 2.36. Структурная схема системы регулирования, состоящей из реактора идеального смешения 1 и пропорционального управляющего устройства (регулятора) 2, представлена на рис. 2.38. Постройте кривые разгона реактора без регулятора и системы регулирования. Как влияет величина параметра Кр настройки регулятора на динамические и статические свойства системы Кр выберите регулирования (для сравнения параметр настройки регулятора равным 1/Ко, 2/Ко, 10/Ко)?

Рис. 2. 2.37. В системе регулирования (рис. 2.38) П-регулятор заменили ПД-регулятором. Получите передаточные функции и постройте кривые разгона системы регулирования с П-регулятором и с ПД-регулятором. Как изменились динамические свойства системы регулирования при замене П-регулятора на ПД-регулятор?

2.38. В реакторе идеального смешения проводят жидкофазную реакцию.

c1, Отклонение концентрации продукта от заданного значения вызванное cо от номинального значения, отклонением концентрации реагента в сырье устраняется регулятором, изменяющим расход сырья F. Уравнение объекта регулирования (реактора) выглядит так:

dc + c = 0,5 c + 1,5 F.

d 1 -68 Определите остаточное отклонение в системе регулирования после единичного ступенчатого изменения концентрации реагента в сырье cо = 1(), если регулятор реализует:

а) П-закон регулирования (Кр = 2);

б) ПД-закон регулирования (Кр = 2, Тd = 8).

Для каждого варианта можно рекомендовать следующий порядок решения:

• Нарисуйте структурную схему системы регулирования;

• Определите её передаточную функцию;

c1(). c1() • Постройте кривую регулирования Значение является остаточным отклонением;

Проверьте, получится ли тот же результат, если воспользоваться свойством 10 из табл. 1 (см. приложение). Как влияет Д-составляющая на процесс регулирования?

2.39. Дана структурная схема промышленного регулятора (рис. 2.39).

Получите его передаточную функцию. Какой закон регулирования реализует регулятор с данной структурной схемой? Постройте кривую разгона регулятора.

Рис. 2. -69 2.40. При исследовании системы регулирования температуры реактора, структурная схема которой приведена на рис. 2.40, провели следующий опыт.

В момент, когда система регулирования находилась в статическом режиме, tзад мгновенно изменили заданное значение температуры на 10°C и с помощью регистрирующего прибора записали изменение температуры реактора t во времени (рис. 2.41).

Рис. 2. Рис. 2. Wo(s), Определите передаточную функцию объекта регулирования если известна передаточная функция регулятора Wp(s) = 1.

-70 2.41. Для определения свойств системы регулирования провели эксперимент, в ходе которого, убедившись, что система регулирования находится в статическом режиме, быстро изменили заданное значение регулируемого параметра на 20% и с помощью регистрирующего прибора y() записали изменение регулируемого параметра во времени (рис. 2.42).

Для сокращения времени эксперимента не дожидались завершения переходного процесса, предполагая определить новое установившееся y значение регулируемого параметра по его текущим значениям в переходном процессе.

Рис. 2. Определите передаточную функцию системы регулирования по каналу задающего воздействия и время переходного процесса.

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

-71 Подсказка 2: Новое установившееся значение регулируемого параметра y можно определить с помощью формулы:

y y y y = lim y ( ) = A C B y + y 2 y A C B 2.42. На рис. 2.43 приведена структурная схема компенсационной системы автоматического регулирования, где:

L(x) – изображение по Лапласу возмущающего воздействия;

• L(z) – изображение по Лапласу управляющего воздействия;

• L(y) – изображение по Лапласу регулируемого параметра;

• Wox(s) – передаточная функция объекта регулирования по • каналу возмущающего воздействия;

Woz(s) – передаточная • функция объекта регулирования по каналу регулирующего воздействия;

Wk(s) – передаточная функция компенсатора.

• Рис. 2. Определить передаточную функцию «идеального» компенсатора, т.е.

такого компенсатора, который обеспечивает равенство нулю реакции объекта y() = 0 на любое возмущающее воздействие x().

-72 2.43. На рис. 2.44 и 2.45 даны структурные схемы комбинированных автоматических систем регулирования, где:

L(yзад) – изображение по Лапласу задающего воздействия;

• Wр(s) – передаточная функция регулятора.

• Остальные обозначения такие же, как в предыдущей задаче (рис. 2.43).

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

Рис. 2. -73 Рис. 2. 2.44. На рис. 2.46 представлена структурная схема регулирования объекта (на схеме отмечен пунктирной рамкой) с взаимосвязанными параметрами:

(z1 и z2), каждый из входных параметров воздействует на оба выходных параметра (y1 и y2).

При расчёте систем регулирования в этом случае необходимо учитывать внутренние связи и другие контуры регулирования. На рис. 2. приведена часть структурной схемы, эквивалентная «объекту R1.

регулирования» по отношению к регулятору Можно привести похожую структурную схему, эквивалентную «объекту регулирования» для регулятора R2.

Определить передаточные функции «эквивалентных объектов» для регуляторов R1 и R2:

L ( y1 ) L ( y2 ) W (s) = W (s) = ;

L( z1 ) L ( z2 ) 1 -74 Рис. 2. Рис. 2. -75 ПРИЛОЖЕНИЕ Таблица Свойства преобразования Лапласа № Оригинал Изображение s d f ( ) e F ( s ) = L[ f ( )] = f ( ) c1 f1 ( ) + c2 f 2 ( ) c1 L[ f1 ( )] + c2 L[ f 2 ( )] s L[ f ( )] f (0) f ‘ ( ) ( n1) n n f (n) ( ) 4 L[ f ( )] s f ( 0 ). f s (0) L[ f ( )] f (t ) dt s f ( c) L[ f ( )] e cs (c 0 ) F (s c) f ( ) e cs 1 s f (c ) (c 0) F c c f f = f (t ) f ( t ) dt F1 ( s) F2 ( s) 9 12 1 lim f ( ) lim s L[ f ( )] 10 s lim f ( ) lim s L[ f ( )] 11 0 s -76 Таблица Преобразование Лапласа некоторых функций № Оригинал Изображение ( ) 1 1( ) s s n!

n s n + ec s+c 1 e T s(Ts + 1) ec + c 1 c2 ( s + c) s s cos( ) s2 + sin( ) s2 + ea eb 10 ( s + a) ( s + b) ba n1e c 11 ( s + c) n (n 1)!

(b c) e a + (a c) e b + (a b) e c (s + a) ( s + b) ( s + c) (a b) (b c) (c a ) -77 Таблица Правила преобразования структурных схем Исходная схем Эквивалентная схема -78 Продолжение таблицы Исходная схем Эквивалентная схема -79 Литература 1.Жукова Г. С., Митрохин С. И., Дарсалия В. Ш. Дифференциальные уравнения / РХТУ им. Д. И. Менделеева. М., 1999. 366 с.

2.Коваль Ж. А., Харитонов Н. И., Шмульян И. К. Сборник упражнений и задач по курсу «Автоматика и автоматизация производства» / МХТИ им. Д. И. Менделеева. М., 1982. 64 с.

3.Полоцкий Л. М., Лапшенков Г. И. Автоматизация химических производств.

М.: Химия, 1982. 296 с.

4.Эрриот П. Регулирование производственных процессов. М.: Энергия, 1967.489 с.

5.Сборник задач по теории автоматического регулирования и управления / Под ред. В. А. Бесекерского. М.: Наука, 1978. 512 с.

6.Справочник по теории автоматического управления / Под ред. А. А. Красовского. М.: Наука, 1987. 712 с.

БЕСПАЛОВ Александр Валентинович ХАРИТОНОВ Николай Иванович ЗОЛОТУХИН Сергей Егорович ФИНЯКИН Леонид Николаевич САДИЛЕНКО Алевтина Степановна ГРУНСКИЙ Владимир Николаевич ДИНАМИЧЕСКИЕ ЗВЕНЬЯ ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ Редактор Н. А. Заходякина Лицензия ПР N 020714 от 02.02.98.

Подписано в печать 2001. Формат 60х84 1/16. Бумага типографская.

Отпечатано на ризографе.

Усл. печ. л.. Уч.-изд. л.. Тираж 500 экз. Заказ.

Российский химико-технологический университет им. Д.И.Менделеева Издательский центр Адрес университета и издательского центра: 125047 Москва, Миусская пл.,

2013 www.netess.ru — «Бесплатная библиотека авторефератов кандидатских и докторских диссертаций»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.

О термопарах: что это такое, принцип действия, подключение, применение

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

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

Устройство и принцип действия

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

Рис. 1. Схема строения термопары

Красным цветом выделено зону горячего спая, синим – холодный спай.

Электроды состоят из разных металлов (металл А и металл В), которые на схеме окрашены в разные цвета. С целью защиты термоэлектродов от агрессивной горячей среды их помещают в герметичную капсулу, заполненную инертным газом или жидкостью. Иногда на электроды надевают керамические бусы, как показано на рис. 2).

Рис. 2. Термопара с керамическими бусами

Принцип действия основан на термоэлектрическом эффекте. При замыкании цепи, например милливольтметром (см. рис. 3) в точках спаек возникает термо-ЭДС. Но если контакты электродов находятся при одинаковой температуре, то эти ЭДС компенсируют друг друга и ток не возникает. Однако, стоит нагреть место горячей спайки горелкой, то согласно эффекту Зеебека возникнет разница потенциалов, поддерживающая существование электрического тока в цепи.

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

В тех случаях, когда необходимо добиться высокой точности показателей, холодные спайки измерительных преобразователей помещают даже в специальные камеры, в которых температурная среда поддерживается на одном уровне специальными электронными устройствами, использующими данные термометра сопротивления (схема показана на рис. 4). При таком подходе можно добиться точности измерений с погрешностью до ± 0,01 °С. Правда, такая высокая точность нужна лишь в немногих технологических процессах. В ряде случаев требования не такие жёсткие и погрешность может быть на порядок ниже.

Рис. 4. Решение вопроса точности показаний термопар

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

Типы термопар и их характеристики

Различные сплавы, используемые для изготовления термопар, обладают разными коэффициентами термо-ЭДС. В зависимости от того, из каких металлов изготовлены термоэлектроды, различают следующие основные типы термопар:

  • ТПП13 – платинородий-платиновые (тип R);
  • ТПП10 – платинородий-платиновые (тип S);
  • ТПР – платинородий-платинродиевые (тип B);
  • ТЖК – железо-константановые (тип J);
  • ТМКн – медь-константановые (тип T);
  • ТНН – нихросил-нисиловые (тип N);
  • ТХА – хромель-алюмелевые (тип K);
  • ТХКн – хромель-константановые (тип E);
  • ТХК – хромель-копелевые (тип L);
  • ТМК – медь-копелевые (тип M);
  • ТСС – сильх-силиновые (тип I);
  • ТВР – вольфрамрениевые (типы A-1 – A-3).

Технические требования к термопарам задаются параметрами определёнными ГОСТ 6616-94, а их НСХ (номинальные статические характеристики преобразования), оптимальные диапазоны измерений, установленные классы допуска регулируются стандартами МЭК 62460, и определены ГОСТ Р 8.585-2001. Заметим, также, что НСХ в вольфрам-рениевых термопарах отсутствовали в таблицах МЭК до 2008 г. На сегодняшний день указанными стандартами не определены характеристики термопары хромель-копель, но их параметры по прежнему регулируются ГОСТ Р 8.585-2001. Поэтому импортные термопары типа L не являются полным аналогом отечественного изделия ТХК.

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

Типы спаев

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

Рис. 5. Типы спаев

Буквами обозначено:

  • И – один спай, изолированный от корпуса;
  • Н – один соединённый с корпусом спай;
  • ИИ – два изолированных друг от друга и от корпуса спая;
  • 2И – сдвоенный спай, изолированный от корпуса;
  • ИН – два спая, один из которых заземлён;
  • НН – два неизолированных спая, соединённых с корпусом.

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

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

Многоточечные термопары

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

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

Таблица сравнения термопар

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

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

Для сравнения термопар составлены таблицы, в которых отображены основные параметры измерительных преобразователей. В качестве примера приводим один из вариантов таблицы для сравнения распространённых термопар.

БЕСПЛАТНАЯ БИБЛИОТЕКА
Тип термопарыKJNRSBTE
Материал положительного электродаCr—NiFeNi—Cr—SiPt—Rh (13 % Rh)Pt—Rh (10 % Rh)Pt—Rh (30 % Rh)CuCr—Ni
Материал отрицательного электродаNi—AlCu—NiNi—Si—MgPtPtPt—Rh (6 % RhCu—NiCu—Ni
Температурный коэффициент40…4155.268
Рабочий температурный диапазон, ºC0 до +11000 до +7000 до +11000 до +16000 до 1600+200 до +1700−185 до +3000 до +800
Значения предельных температур, ºС−180; +1300−180; +800−270; +1300– 50; +1600−50; +17500; +1820−250; +400−40; +900
Класс точности 1, в соответствующем диапазоне температур, (°C)±1,5 от −40 °C до 375 °C±1,5 от −40 °C до 375 °C±1,5 от −40 °C до 375 °C±1,0 от 0 °C до 1100 °C±1,0 от 0 °C до 1100 °C±0,5 от −40 °C до 125 °C±1,5 от −40 °C до 375 °C
±0,004×T от 375 °C до 1000 °C±0,004×T от 375 °C до 750 °C±0,004×T от 375 °C до 1000 °C±[1 + 0,003×(T − 1100)] от 1100 °C до 1600 °C±[1 + 0,003×(T − 1100)] от 1100 °C до 1600 °±0,004×T от 125 °C до 350 °C±0,004×T от 375 °C до 800 °C
Класс точности 2 в соответствующем диапазоне температур, (°C)±2,5 от −40 °C до 333 °C±2,5 от −40 °C до 333 °C±2,5 от −40 °C до 333 °C±1,5 от 0 °C до 600 °C±1,5 от 0 °C до 600 °C±0,0025×T от 600 °C до 1700 °C±1,0 от −40 °C до 133 °C±2,5 от −40 °C до 333 °C
±0,0075×T от 333 °C до 1200 °C±0, T от 333 °C до 750 °C±0,0075×T от 333 °C до 1200 °C±0,0025×T от 600 °C до 1600 °C±0,0025×T от 600 °C до 1600 °C±0,0075×T от 133 °C до 350 °C±0,0075×T от 333 °C до 900 °C
Цветовая маркировка выводов по МЭКЗелёный — белыйЧёрный — белыйСиреневый — белыйОранжевый — белыйОранжевый — белыйОтсутствуетКоричневый — белыйФиолетовый — белый

Способы подключения

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

Некоторые измерительные приборы содержат схемы корректировки показаний на основе встроенного термистора. К таким приборам просто подключаются провода, соблюдая их полярность (см. рис. 6).

Рис. 6. Компенсационные провода

Часто используют схему подключения «на разрыв». Измерительный прибор, подключают через проводник того же типа что и клеммы (чаще всего медь). Таким образом, в местах соединения отсутствует холодный спай. Он образуется лишь в одном месте: в точке присоединения провода к электроду термопары. На рисунке 7 показана схема такого подключения.

Рис. 7. Схема подключения на разрыв

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

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

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

Применение

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

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

Преимущества

  • высокая точность измерений;
  • достаточно широкий температурный диапазон;
  • высокая надёжность;
  • простота в обслуживании;
  • дешевизна.

Недостатки

Недостатками изделий являются факторы:

  • влияние свободных спаев на показатели приборов;
  • ограничение пределов рабочего диапазона нелинейной зависимостью ТЭДС от степени нагревания, порождающей сложности в разработке вторичных преобразователей сигналов;
  • при длительной эксплуатации в условиях перепадов температур ухудшаются градуировочные характеристики;
  • необходимость в индивидуальной градуировке для получения высокой точности измерений, в пределах погрешности в 0,01 ºC.

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


источники:

http://netess.ru/3metodichki/1252337-1-ministerstvo-obrazovaniya-rossiyskoy-federacii-rossiyskiy-himiko-tehnologicheskiy-universitet-d-mendeleeva-dinamicheskie-zv.php

http://www.asutpp.ru/termopary.html