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

МУ к лаб_работе 3

.pdf
Скачиваний:
7
Добавлен:
22.05.2015
Размер:
103.58 Кб
Скачать

Методические указания к лабораторной работе 1.5.11 "Экспериментальное определение динамической характеристики прибора" 2010-2011 учебный год

Нормативные ссылки В методических указаниях к лабораторной работе использованы

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

ГОСТ Р 8.585-2001 ГСИ. Термопары. Номинальные статические характеристики преобразования.

ГОСТ 2.301-68 ЕСКД. Форматы.

ГОСТ 21.101-97 СПДС. Основные требования к проектной и рабочей документации.

ГОСТ 6616-74 Преобразователи термоэлектрические ГСП. Общие технические условия.

РМГ 29-99 ГСИ. Метрология. Термины и определения.

1 Цели работы Изучение методики экспериментально-аналитического определения

динамической характеристики термоэлектрического термометра. Изучение методики обработки экспериментальных данных.

2 Общие сведения При исследовании и разработке систем автоматического регулиро-

вания (САР) решаются задачи анализа систем или их синтеза. В первом случае имеется конкретная система и требуется определить ее свойства. Во втором случае сначала задаются свойства, которыми должна обладать система, т.е. требования к ней, а затем создается САР, удовлетворяющая этим требованиям [1].

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

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

В лабораторной работе объектом исследований является термоэлектрический термометр, состоящий из датчика – термоэлектрического преобразователя и подключенного к нему термоэлектродным кабелем вторичного прибора – электронного автоматического потенциометра.

2

Термоэлектрический преобразователь ТХК-0515 по ГОСТ 6616 имеет номинальную статическую характеристику типа L по ГОСТ Р 8.585.

Электронный автоматический потенциометр типа КСП4 реализует функции показаний и регистрации измеряемой температуры. Прибор имеет следующие метрологические характеристики: диапазон измерений по температуре от минус 50 до плюс 200 ОС, класс точности 0,5; номинальную статическую характеристику (НСХ) типа L.

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

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

T

2

d 2 N (t)

+ T

d N (t)

+ N (t) = K θ (t),

(2.1)

 

 

2

dt

2

1

dt

 

 

 

 

где T22 - постоянная времени, с2; T1 - постоянная времени, с;

К - коэффициент передачи (чувствительность), ОС/ОС;

N(t) - временная функция изменения показаний вторичного прибора (выходной величины), ОС;

θ (t) - временная функция изменения температуры (входной величины), ОС.

3

Передаточная функция, соответствующая уравнению (2.1), имеет

вид

W ( p) =

N ( p) =

 

 

 

K

,

(2.2)

2

 

2

 

 

 

θ ( p)

 

+ T1 p + 1

 

 

 

T2

p

 

 

 

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

Коэффициент передачи (чувствительность) К и постоянные времени T22 и T1 - это параметры математической модели, подлежащие определению при выполнении лабораторной работы.

Если быстро перенести датчик - термоэлектрический преобразователь из нагретого термошкафа в сосуд с холодной водой, то его входная величина - температура θ изменится ступенчато от начального значения

θн до конечного значения θк, а выходная величина термометра - показа-

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

чений: Δθi = θi - θн;

Ni = Ni – N н.

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

Δθ0 = θк - θн;

Nуст = Nк – N н.

Для термоэлектрического термометра изменение показаний на выходе после окончания переходного процесса (новое установившееся значение Nуст) будет равно по размеру ступенчатому изменению темпера-

туры на входе Δθ0, поэтому коэффициент передачи (чувствительность) термометра К=1 ОС/ОС.

Одним из методов определения значений параметров математической модели динамики термометра по его кривой разгона является "метод площадей", изложенный в монографии [3]. В соответствии с ним ось абсцисс над кривой разгона необходимо разбить на n равных по длине отрезков с интервалом времени t - так, чтобы на этих интервалах кривая разгона мало отличалась от прямой линии. В графу 1 таблицы 2.1 необхо-

димо занести значения моментов времени ti, при этом t0 = 0, t1 =

t, t2 =

2 t и т.д. В графу 2 нужно записать отклонения показаний прибора

Ni от

начального значения в конце каждого интервала.

 

4

Затем значения отклонений DNi в конце каждого интервала следу-

ет разделить на установившееся значение DNуст и полученные значения

 

σ (i ×

t) =

Ni

 

 

 

 

(2.3)

 

N уст

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

занести в графу 3 таблицы 2.1.

 

 

 

 

Таблица 2.1 - Промежуточные результаты расчетов

 

 

ti, с

DN , ОС

 

s(i×Dt)

 

1-s(i×Dt)

i×DQ

1- i×DQ

 

[1-s(i×Dt)](1- i×DQ)

 

i

 

 

 

 

 

 

 

 

 

 

1

2

 

 

3

 

4

5

6

 

7

0

0

 

 

0

 

1

0

1

 

1

t1

DN1

 

s(Dt)

 

1-s(Dt)

Dt/T1

1-DQ

 

[1-s(Dt)](1-DQ)

t2

DN2

 

s(2×Dt)

 

1-s(2×Dt)

2Dt/T

1-2×DQ

 

[1-s(2×Dt)](1-2×DQ)

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tn

DN

 

s(n×Dt)

 

1-s(n×Dt)

nDt/T

1- n×DQ

 

[1-s(n×Dt)](1- n×DQ)

 

уст

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

В графу 4 таблицы 2.1 требуется записать значения разностей [1-s(i×Dt)]. После этого может быть рассчитано значение постоянной времени Т1 по формуле

 

n

-σ (i

 

 

T = Dt ∑[1

× Dt)] - 0,5 .

(2.4)

1

i=0

 

 

 

Затем

необходимо

ввести новую

безразмерную переменную

Q = t/T1, рассчитать значения произведений i×DQ и разностей (1- i×DQ), в

которых DQ=Dt/T1, и записать их в графу 6 таблицы 2.1. Далее нужно рассчитать значения произведений [1-s(i×Dt)](1- i×DQ) и записать результаты в графу 7 таблицы 2.1. После этого может быть рассчитано значение постоянной времени T22 по формуле

T22 = T1

n

 

 

× Dt × [1 - σ (i × Dt)](1 - i × DQ) - 0,5 .

(2.5)

 

i=0

 

 

Проверка адекватности математической модели динамики термометра проводится сравнением ординат кривой разгона и переходной характеристики, являющейся решением уравнения (2.1), при значении входной величины Dq(t) = Dqо = q к - q н.

Вид переходной характеристики зависит от вида корней характеристического уравнения

 

5

T22λ2 + T1λ +1 = 0.

(2.6)

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

 

 

 

 

 

λ

2

 

λ1t -

DN

расч

(t) = К × Dθ

0

1 +

 

е

 

 

 

 

 

λ1 - λ2

 

 

 

 

 

 

 

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

λ

λ

 

t

 

1

е

2

.

(2.7)

 

λ1 - λ2

 

 

 

 

 

 

 

λ1 = −α1;

λ2 = −α2 .

то переходная функция может быть записана как

 

 

(t) = К × Dθ

 

 

 

α 2 exp(-α1t)

 

α1 exp(-α 2t)

 

DN

расч

0

1

-

+

.

 

 

 

 

 

 

α

2

- α

 

α

2

- α

1

 

 

 

 

 

 

 

 

1

 

 

 

 

Для проверки адекватности математической модели динамики термометра необходимо рассчитать значения расхождений "ε" между переходной характеристикой и кривой разгона по формуле

ε =

N расч(ti ) − Ni

100

.

(2.8)

 

 

N уст

 

 

 

При выполнении условия εмакс < ±5 % математическая модель динамики термометра считается адекватной.

3 Описание лабораторной установки Лабораторная установка содержит: термоэлектрический термометр,

состоящий из смонтированного на стенде электронного потенциометра типа КСП4 и датчика типа ТХК-0515, подключенного к потенциометру гибким термоэлектродным кабелем. Скорость перемещения диаграммы потенциометра установлена равной 2400 мм/ч. Рядом со стендом на полу размещен термошкаф, на столе находится сосуд с холодной водой.

4 Техника безопасности при выполнении лабораторной работы Пакетные выключатели, с помощью которых подается напряжение

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

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

6

Запрещается выдвигать шасси прибора КСП4 при поданном напряжении питания.

На лабораторном столе не должно быть лишних посторонних предметов (пакетов, сумок, одежды, пищевых продуктов).

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

5 Порядок выполнения работы

5.1Подготовьте лабораторную установку к проведению эксперимента, выполните для этого следующие операции:

- начертите заготовку таблицы 2.1 для записей результатов наблюдений и расчетов, укажите в графе 1 значения времени с интервалом 10 с; - поместите датчик ТХК-0515 в вертикальное отверстие термошка-

фа;

- после включения преподавателем питания потенциометра КСП4

проследите за возрастающими показаниями прибора до достижения ими значений 80…90 ОС.

5.2Проведите эксперимент по снятию кривой разгона термометра, выполните для этого следующие операции:

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

- напишите карандашом на диаграмме рядом с прямой линией соответствующее значение показания прибора Nн (что необходимо из-за не-

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

-быстро перенесите датчик ТХК-0515 из нагретого термошкафа в сосуд с холодной водой, проследите за регистрацией потенциометром уменьшающихся значений температуры до тех пор, пока на диаграмме не будет вычерчиваться прямая линия;

-напишите карандашом на диаграмме рядом с прямой линией соответствующее значение показания прибора Nк ;

-после выключения преподавателем электропитания прибора КСП4 осторожно оторвите часть диаграммы с записью кривой разгона.

5.3 Проведите обработку экспериментальных данных, выполните для этого следующие операции:

-выберите интервал разбиения t=10 секунд на оси абсцисс под

кривой разгона – при их общем количестве n=6…15;

определите, сколько

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

t (при скорости про-

тяжки 2400 мм/ч диаграмма перемещается на 6,7 мм за 10 с) – разметьте карандашом интервалы разбиения на диаграмме;

7

- определите значения DNi на кривой разгона, соответствующие моментам времени t=i×Dt (i=0,…, n); при этом, очевидно, что при tо=0 показание Ni будет равно начальному значению Nн; значения ti запишите в графу 1 таблиц 2.1 и 2.2, значения DNi запишите в графу 2 таблиц 2.1 и 2.2;

Таблица 2.2 - Исходные данные и результаты проверки адекватности

ti, c

DNi, ОС

DNрасч(ti), ОС

ε, %

1

2

3

4

t0

 

 

 

t1

 

 

 

t2

 

 

 

 

 

 

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

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

дели динамики термометра; при этом значения DNрасч(ti), рассчитанные

формуле (2.7), запишите в графу 3 таблицы 2.2; значения расхождения "ε" между переходной характеристикой и кривой разгона, рассчитанные по формуле (2.8), запишите в графу 4 таблицы 2.2;

-проверьте выполнение условия εмакс < ±5 %; сделайте заключение об адекватности математической модели динамики термометра;

-постройте график переходной характеристики (2.7) рядом (совмещенно) с кривой разгона.

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

T1.1

d N (t)

+ N (t) = K θ (t),

(2.9)

 

 

dt

 

в котором постоянная времени Т1.1 найдена упрощенным способом: как время, за которое отклонение DNi достигнет 0,632 от своего установившегося значения. Для этого случая рассчитайте значения переходной характеристики DNрасч.1(ti), по формуле

8

 

 

 

 

 

-

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

DN

 

(t) = К × Dθ

 

 

- е

 

 

расч.1

0

1

1.1

 

(2.10)

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишите значения этой переходной характеристики в отдельную таблицу 3, аналогичную таблице 2.2, рассчитайте значения отклонений переходной характеристики от кривой разгона, сделайте заключение об адекватности математической модели (2.9). Постройте график переходной характеристики (2.10) рядом (совмещенно) с кривой разгона.

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

T1.2

d N (t)

+ N (t) = K θ (t),

(2.11)

 

 

dt

 

в котором постоянная времени Т1.2 найдена упрощенным способом: как одна третья часть времени, за которое отклонение Ni достигнет 0,95 от своего установившегося значения. Для этого случая рассчитайте значения переходной характеристики Nрасч.2(ti), по формуле

 

 

 

 

 

-

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

DN

 

(t) = К × Dθ

 

 

- е

 

 

расч.2

0

1

1.2

 

(2.12)

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишите значения этой переходной характеристики в отдельную таблицу 4, аналогичную таблице 2.2, рассчитайте значения отклонений переходной характеристики (2.12) от кривой разгона, сделайте заключение об адекватности математической модели (2.11).

6 Содержание отчета Отчет о лабораторной работе должен содержать:

-название работы;

-три раздела со следующими заголовками: 1 Цели работы

2 Технические данные используемых средств измерений

3 Порядок выполнения работы, полученные результаты и выводы по работе.

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

9

зать метрологические характеристики (тип, диапазон измерений, диапазон показаний или диапазон изменения выходного сигнала, НСХ, класс точности или пределы погрешности). В разделе 3 должна быть дана информация

овыполнении каждого пункта лабораторной работы, а также:

-таблицы 2.1, 2.2, 2.3, 2.4;

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

-график ступенчатого изменения температуры Δθ(t); под ним график переходной характеристики (2.7), построенный рядом (совмещено) с кривой разгона;

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

нения (2.1);

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

(2.2);

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

-еще один график ступенчатого изменения температуры Δθ(t); под ним график переходной характеристики (2.10), построенный рядом (совмещено) с кривой разгона;

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

нения (2.9);

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

(2.10);

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

нения (2.11);

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

(2.12).

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

Отчет должен быть оформлен на листах бумаги формата А4 по ГОСТ 2.301, снабженных рамкой и основной надписью по ГОСТ 21.101.

В приложении А на рисунке А.1 показана форма основной надписи для первого листа отчета, на рисунке А.2 – форма основной надписи для последующих листов.

10

7 Контрольные вопросы

1)Какие функции выполняет прибор типа КСП4?

2)Что называется динамической характеристикой средства измере-

ний?

3) В каких формах могут быть представлены динамические характеристики элементов автоматических систем?

5)Как могут быть записаны в аналитической форме теоретические динамические характеристики термоэлектрического термометра?

6)Что называется чувствительностью (коэффициентом передачи) средства измерений с линейной статической характеристикой?

7)В чем состоит суть экспериментального активного метода получения математической модели динамики элемента по кривой разгона?

8)Что называется кривой разгона?

9)Что понимается под аппроксимацией кривой разгона?

10)Какая математическая модель динамики считается адекватной?

11)Как реализуется процедура проверки адекватности математической модели динамики элемента?

12)Как аналитически решается дифференциальное уравнение первого (второго) порядка?

13)Как изменится передаточная функция термометра, если при снятии кривой разгона датчик переместить из сосуда с холодной водой в нагретый термошкаф?

14)Как может быть записана упрощенная передаточная функция термометра, представленного в виде инерционного звена первого порядка?

15)Как может быть построен график изменения показаний термометра во времени при импульсном изменении температуры?

Список литературы

1Фарзане Н.Г., Илясов Л.В., Азим-Заде А.Ю. Технологические измерения и приборы. М. : Высшая школа, 1989. 456 с.

2Соколов В.А. Автоматизация технологических процессов пищевой промышленности. М. : Агропромиздат, 1991. 445 с.

3Стефани Е.П. Основы расчета настройки регуляторов теплоэнергетических процессов. М. : Энергия, 1972. 376 с.