Процесс создания математической модели объекта

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

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

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


Задание на курсовую работу.

1. По экспериментальным данным найти математическую модель объекта (системы) методом площадей в виде звена 2-го или 3-го порядка. Оценить точность аппроксимации.

2. Найти и построить частотные характеристики объекта (АЧХ, ФЧХ, АФЧХ) и провести подробный анализ этих характеристик.

3. По заданному закону регулирования найти математическую модель ЗСАУ.

4. Определить устойчивость ЗСАУ по одному из критериев. Если система неустойчива, то добиться ее устойчивости.

5. Найти переходную функцию ЗСАУ и построить ее. Найти по кривой основные ПКР.

6. Привести структурную схему САУ в соответствии с требованиями ГОСТ.

7. Дать выводы по работе.

t00.250.30.350,50,7511.251,51.7522.252,5
s00.050.070.090.110.130.160.190.210.250.290.350.4
t2.733,24,254.54.7555,255.55.756
s0.450.50.550.60.70.750.80.850.90.951
Регулятор ПИД: Кп = 1;  Ки = 0,5;  Тд = 15сек; К=2

Объект регулирования АПЗ-2 


1. Обработка исходных данных Методом площадей.

Данный метод был разработан М.П. Симою. Метод служит для определения передаточной функции объекта по  экспериментальной кривой разгона.

В основе метода  лежит предположение, что исследуемый объект может быть описан линейным дифференциальным уравнением с постоянными коэффициентами:

 (2)

W(p) = S bipi / Sajpj (-),  b0 = 0,  a0 = 1

Задача состоит в том,  чтобы определить неизвестные коэффициенты

а1¸аn и 1¸bm.

Коэффициенты ai будут определяться по следующим формулам:

а1 = F1 + 1

а2 = F2 + 2 + F11

а3 = F3 + 3 + 2F1 + 1F2

…………….

аi = Fi + bi + SbjFi-j

В системе уравнений, приведенной выше i = m + n. Составляющие элементы системы  определяются из следующих формул:

F1 = Dt{S(1-) – 0.5}

F2 = F12DQ{S(1 - )*(1 - Q) – 0.5}

F3 = F13DQ{S(1 - )*(1 – 2*Q + Q2/2) – 0.5}

и т.д.

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

По исходной кривой значения Yi для каждого значения времени заносим в таблицу Exсel и находим значения, необходимые для вычисления значений Fi.

Исходя из полученных данных, имеем:

F1 = 3,2875,  F2 = 5,31953,  F3 = 7,30796. F4 = -7,61321

По полученным значениям видно, что разница между F3 и F4 существенная, при этом F4 является числом отрицательным, что дает нам основание говорить о том, что значение коэффициента а4 = 0.

Исходя из приведенных выше формул нахождения аi , получаем коэффициенты 1, a1, a2, а3:

1 = 1,042; a1 = 4,32927;  a2 = 8,74435, а3 = 12,8497.

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

W(p) = (1,042 + 1)/(12,8497р3 + 8,74435р2 + 4,32927р + 1).

Построим данную передаточную функцию в пакете VisSim, получим характеристику и найдем все ошибки (среднеквадратическое отклонение, абсолютную и относительную (приведенную) ошибки). График полученной характеристики приведен в приложении.

tУэУр

hi = Уэ – Ур

Dhi2

10,250,050,060,010,0001
20,50,110,130,020,0004
310,160,180,020,0004
41,50,210,240,030,0009
520,290,30,010,0001
630,50,5180,0180,000324
74,250,60,730,130,0169
84,50,70,80,10,01
950,80,850,050,0025
105,50,90,940,040,0016
115,750,950,9530,0030,000009
Актуально: