Курс_АСУ_2013
.pdf3.5. Уравнения динамики САРЧ.
Задачей этого раздела является выбор уравнений динамики звеньев
САРЧ [3, 4].
Для выполнения расчета необходимо представить уравнение динамики заданного объекта регулирования. Уравнение динамики двигателя можно принять в виде
Tа (d φ) / (d t) + βд φ = μ – λ (t) , |
(1) |
где Та – постоянная времени объекта регулирования, характеризующая его инерционные свойства (по ГОСТ 10511 этот параметр называется временем разгона); βд - коэффициент самовыравнивания; φ, μ, λ –
относительные координаты частоты вращения, регуляторного и нагрузочного воздействий на объект, соответственно. В дальнейшем расчете примем во внимание только регуляторное воздействие, исключив влияние нагрузки λ (t).
Выбрав уравнение в таком упрощенном виде, будем полагать, что влияние агрегатов наддува, впускного ресивера и инерционных свойств топливной аппаратуры при сбросе нагрузки пренебрежимо мало. Элемен-
ты регулятора представляются в виде типовых динамических звеньев,
описываемых типовыми уравнениями с входными и выходными
координатами, выраженными в о.е.. Например:
а) чувствительного элемента (датчика частоты вращения)
Tr2 (d2η) / (dt2) + Тк (dη) / (dt) + r η = kφφ, |
(2) |
где Tr - постоянная времени, характеризующая влияние массы измерителя;
Tк - постоянная времени, характеризующая силы вязкого трения; φ -
координата частоты вращения; η - координата муфты измерителя; r
11
Таблица 2
Исходные данные к определению момента инерции КШМ
№ |
Тип дизеля |
Размеры деталей КШМ дизеля, см |
|
|
|
Масса, кг, |
|||
варианта |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Радиус |
Диаметр |
Длина |
Ширина |
Толщина |
Приведенный |
Масса |
Масса |
|
|
кривошипа, |
шатунной |
шатунной |
щеки, b |
щеки, q |
радиус |
поршня |
шатуна |
|
|
r |
шейки, d2 |
шейки, l2 |
|
|
вращения, R |
Gп |
Gш |
1 |
6ЧН32/35 |
17.5 |
23 |
21 |
34 |
8.7 |
32 |
134 |
110 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
2 |
4ЧН26/26 |
13 |
19 |
17 |
28 |
7.9 |
24 |
40 |
56,5 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
3 |
6ЧРН36/45 |
22.5 |
30 |
28 |
36.5 |
8.2 |
33.5 |
143 |
243 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
4 |
12ЧН15/18 |
9 |
12.5 |
13.5 |
22 |
4 |
19 |
2,38 |
05,5 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
5 |
6ЧН |
16.5 |
24 |
22 |
31 |
8.1 |
30.5 |
60,8 |
102,5 |
|
31,8/33 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
6 |
8ЧН19/21 |
10.5 |
13.5 |
11 |
25 |
4.5 |
21.8 |
21,2 |
48 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
7 |
8ЧН16,5/18 |
11 |
11.5 |
10.5 |
14 |
3.9 |
15 |
17,4 |
16,7 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
8 |
6ЧН30/38 |
19 |
22 |
27 |
20 |
7.2 |
31 |
76,4 |
70,6 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
9 |
8Ч23/30 |
15 |
16 |
14 |
36 |
8.2 |
27.5 |
33,4 |
24,5 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
0 |
8ЧН 25/34 |
17 |
15 |
12.5 |
28 |
6.7 |
29 |
36,2 |
56,5 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12 |
13
Таблица 3
Исходные данные для уравнений динамики регулятора (по предпоследней цифре шифра)
№ |
Давление |
|
|
|
Параметры регулятора |
|
|
||||
|
наддува, |
|
|
|
|
|
|
|
|
|
|
|
рк, мПа |
|
|
|
|
|
|
|
|
|
|
|
|
Тr, |
Тk, |
Тs, |
Тi, |
, |
k , |
ks, |
kжос, |
ki, |
|
|
|
с |
с |
с |
с |
r |
|
|
|
о.е. |
|
|
|
о.е |
о.е |
о.е. |
о.е. |
||||||
1 |
0.17 |
|
0,08 |
0,20 |
0,20 |
0,02 |
1,5 |
8 |
5 |
7 |
|
2 |
0.18 |
0 |
0,120 |
0,18 |
0,60 |
0,03 |
1,4 |
9 |
3 |
6 |
|
|
|
|
|
|
|
|
|
|
|
||
3 |
0.11 |
0,04 |
0,25 |
0,50 |
0,04 |
1,3 |
10 |
4 |
5 |
||
|
|||||||||||
4 |
0.14 |
|
0,05 |
0,30 |
0,40 |
0,05 |
1,2 |
11 |
5 |
4 |
|
5 |
0.14 |
0,04 |
0 |
0,30 |
0 |
0 |
1,0 |
12 |
0 |
3 |
|
6 |
0.16 |
0 |
0,06 |
0,15 |
0,30 |
0,05 |
0,8 |
11 |
4 |
2 |
|
7 |
0.13 |
0,06 |
0 |
0,22 |
0 |
0,04 |
1,0 |
10 |
5 |
3 |
|
8 |
0.16 |
0 |
0,07 |
0,12 |
0,40 |
0,03 |
1,2 |
9 |
6 |
4 |
|
9 |
0.13 |
0 |
0,05 |
0,22 |
0,50 |
0,02 |
1,4 |
8 |
7 |
5 |
|
0 |
0.15 |
0,09 |
0,15 |
0,60 |
0,03 |
1,5 |
7 |
8 |
6 |
||
|
|||||||||||
х |
– нулевое значение постоянных времени понижает порядок |
||||||||||
|
|||||||||||
дифференциального уравнения. |
|
|
|
|
|
|
– степень неравномерности ИЧВ; kφ - коэффициент неравномерности измерителя, примем равным 1.
б) сервомотора, состоящего из управляющего элемента – золотника
μ =kμ η, |
(3) |
где μ - координата золотника; kμ – коэффициент передачи
и исполнительного механизма (ИМ) |
|
Tѕ (d )/( dt) =ksμ, |
(4) |
где - координата исполнительного механизма; Ts - постоянная времени ИМ.
в) жесткой обратной связи
жос = −kжос |
(5) |
где жос - координата на выходе ЖОС; kжос – коэффициент усиления ЖОС;
в) гибкой обратной связи |
|
Ti (di ) /( dt)= −ki |
(6) |
13
14
где i – координата на выходе ГОС; Ti - постоянная времени ГОС
(изодрома); ki – коэффициент усиления ГОС; |
|
г) уравнение передаточного звена |
|
h = кмп , |
(7) |
где кмп – коэффициент передачи, принимаемый для всех вариантов равным
0,66; h – координата рейки ТНВД.
Исходные данные с вариантами постоянных величин в уравнениях динамики элементов регулятора приведены в приложении 6.
3.6. Построение структурной схемы регулятора и САРЧ Для получения структурной схемы регулятора необходимо
преобразовать принципиальную схему регулятора (Приложение 3) в
структурную схему (рис.1), обозначив входные и выходные параметры звеньев схемы, используя соответствующие знаки сложения и вычитания входных величин и полученные статические характеристики элементов.
Первоначально следует определить в принципиальной схеме регулятора основные функциональные элементы: датчики частоты вращения,
усилительные элементы, передаточные механизмы, упруго-
присоединенные катаракты и др., - и определить последовательность передачи сигнала от измерителя частоты вращения к выходному элементу регулятора. В этой последовательности и следует расположить элементы структурной схемы в виде прямоугольников, обозначив стрелками направление прохождения сигнала. Таким образом, будет получена разомкнутая, а затем и замкнутая, схемы регулятора. Внутри этой схемы могут быть дополнительные контуры при наличии прямых и обратных связей.
3.7. Определение передаточных функций (ПФ).
Для получения ПФ необходимо выполнить преобразование дифференциальных уравнений, используя преобразование Лапласа.
14
15
При этом полагают, что d/dt=p. Тогда, к примеру, ПФ ОР получим следующим образом:
Рис.1. Структурная схема САРЧ ГД
от Tа (d φ) / (d t) + βд φ = μ перейдем к Tа ·p·φ(t)+ βд φ(t)= μ(t) или (Tа ·p· + βд)· φ(t)= μ(t). Тогда ПФ ОР примет вид
Wор(p)= |
( ) |
= |
1 |
|
. |
(8) |
( ) |
∙ |
+1 |
||||
|
|
|
|
|
|
|
Далее, также для примера, покажем преобразование дифференциального уравнения ИЧВ: от уравнения Tr2 (d2η) / (dt2) + Тк (dη) / (dt) + η = kφφ
перейдем у уравнению в операторной форме Tr2 ·p2·η(t)+ Тк ·p · η(t)+ η(t)= kφφ(t) или
(Tr2 ·p2+ Тк ·p + 1) η(t)= kφφ(t). Тогда ПФ ИЧВ примет вид
Wичв(p)= |
( ) |
= |
|
|
|
|
. |
(9) |
( ) |
2 |
∙2+ |
|
∙+1 |
||||
|
|
|
|
|
|
|
|
15
16
Подобным же образом получим ПФ остальных звеньев.
Затем, используя правила получения передаточных функций цепей при различных видах соединений звеньев структурной схемы (рис.1),
последовательно определим ПФ регулятора частоты вращения, а затем,
зная ПФ ОР, и ПФ САРЧ.
Определение передаточных функций регулятора и САРЧ вначале следует выполнить в общем виде. И только после их получения подставить в ПФ их ранее полученные выражения (8), (9) и другие, а также заданные значения параметров регуляторов (табл.3).
3.8. Определение параметров пропульсивного комплекса В данном разделе следует произвести вычисления, необходимые для
определения параметров пропульсивного комплекса. Это позволит определить постоянную времени Ta и коэффициент самовыравнивания β двигателя – величины, необходимые для использования уравнения динамики двигателя.
3.8.1. Динамические режимы работы двигателя как объекта регулирования определяются дифференциальным уравнением (1)
Tа (d φ) / (d t) + βд φ = μ – λ (t)
|
|
J |
|
|
Мс |
|
Мд |
|
||
где |
Та |
|
- время разгона двигателя, с; |
F |
|
|
|
|
|
- фактор |
|
|
F |
|
|
|
|
|
|
|
устойчивости двигателя.
Физический смысл коэффициентов дифференциального уравнения состоит в следующем:
Время разгона двигателя (Tа) - характеризует его инерционные свойства и зависит от фактора устойчивости. Она определяет время условного переходного процесса, в котором изменение угловой частоты вращения (ω) происходило бы с постоянной скоростью, равной ее величине в начальный период реального переходного процесса.
16
17
Вуравнениях переменные представлены в относительных единицах
ввиде безразмерных параметров.
φ= ∆ωi /ωr, η =∆hi /hr, f= ∆λi/λr,,
где ωr, hr, λr - номинальные значения координат/номинальные значения соответствующих параметров. В безразмерном виде после деления на
коэффициент k при координате η уравнение динамики двигателя без
наддува, работающего на винт фиксированного шага (ВФШ) примет вид
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
kэ f |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
Та dt |
|
|
|
|
m |
f |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
М |
с |
|
|
|
Мд |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
где |
β= |
|
|
|
|
0 |
|
|
|
|
|
|
0 |
|
|
|
(tg 1 - tg 2 ) 0/ |
Мд0 |
коэффициент |
||||
|
|
|
|
|
М |
|
|
|
|
|
|
|
|
|
|
|
- |
||||||
|
|
|
|
|
|
|
д0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
самовыравнивания двигателя; |
|
|
k э |
= k |
f |
/ k |
- коэффициент |
усиления по |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
|
|
|
|
|
|
||
каналу |
нагрузки; |
T T / k |
|
= |
J 0 |
|
- |
|
постоянная |
времени двигателя, |
|||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
a |
d |
|
|
М д0 |
|
|
|
|
|
|
|
|
определяющая длительность переходного процесса при ступенчатом единичном воздействии по каналу топливоподачи.
Расчет значения коэффициента β уравнения динамики двигателя сводится к определению частных производных (рис.2), определяемых либо графически построением касательных к построенным на рис. кривым в точке их пересечения и определением тангенсов углов α1 и α2 их наклона либо аналитически [3, 5].
3.8.2. Определение коэффициентов и постоянных времени уравнения динамики ГД. Для определения искомых значений необходимо построить статические характеристики двигателя и нагрузки.
Для определения постоянной времени двигателя Та необходимо:
- в относительной системе координат построить зависимости
Рc=f(ω), Ре=f(ω) (рис.1),
17
18
Рис.2. К определению координат точки О установившегося режима работы ГД
- определить на их пересечении точку рабочего режима 0 и ее
координаты (Ро, ωо);
|
- построить в точке O установившегося режима касательные к |
|||||||
характеристикам Рс(ω) и Рд(ω), |
|
|||||||
|
- в равносторонней системе координат определить величины |
|||||||
углов α1 и α2 |
и их тангенсы; |
|
|
|||||
|
- определить величину коэффициента самовыравнивания [5] с |
|||||||
учетом полученного фактора устойчивости F |
|
|||||||
|
Рe0 |
F |
|
Рe0 |
(tg tg ) |
(1) |
||
|
2 |
|||||||
|
Р |
|
1 |
2 |
|
|||
|
|
|
|
|
||||
|
0 |
|
0 |
|
|
|
18
19
где Ре - номинальное значение мощности ГД, в относительных единицах; ω
-значение угловой частоты вращения вала ГД, в относительных единицах;
-вычислить значение постоянной времени двигателя Tд, c,
Та=J/FP,
где J - значение приведенного момента инерции вращающихся масс,
которое определяется в относительных единицах в последующих
расчетах.
Для построения зависимости мощности ГД от частоты вращения
Pе=f(ω) используется эмпирическая формула, действительная в диапазоне
0<ω<1.1ωr [2] при условии предположения Pe=k·h
|
i |
|
P P 0,5 |
||
|
||
e ei |
r |
|
|
|
|
|
i |
2 |
|
i |
3 |
|
3 |
|
1,5 |
|
|
|
|
|||
|
r |
|
|
r |
|
|
или Рi=Ре·(ωi/ ωr) , |
|
|
|
|
|
|
|
|||
|
|
|
|
|
где Pе - мощность ГД при заданном положении рейки ТНВД, кВт или о.е.; ωi - текущее (задаваемое) значение угловой частоты вращения вала ГД, c-1; ωr - номинальное значение угловой частоты вращения вала ГД, c-1.
Результаты расчетов для Pе=f(ω) при h=1.0·P0 сводятся в таблицу 4.
Теперь в соответствии с исходными данными (см. табл.4)
необходимо определить зависимость мощности сопротивления – нагрузки
– от угловой частоты вращения (винтовая характеристика). Она определяется выражением для частичной винтовой характеристики
Р0,8=0,8∙Ре· (ni/nr)3;
где P0,8 – значение координаты частичной характеристики двигателя при
ni=0,8·nr.
При построении зависимостей использованы характеристики СЭУ,
приведенные в [17].
3.8.3 Определение параметров установившегося режима. По исходным данным, приведенным в таблице 1, могут быть построены статические характеристики ГД и нагрузки в виде винтовой характеристики (рис. 4). Точка пересечения характеристик О характеризует установившийся режим ГД по частоте вращения (и только
19
20
по частоте вращения). Определив координаты точки О (Ро, ωо), можем произвести вычисление одной из постоянных величин уравнения динами ГД - коэффициента самовыравнивания β - по формуле (1). Для вычисления же второй постоянной величины - Ta – знания координат точки О недостаточно. Необходимо определить момент инерции вращающихся масс пропульсивного комплекса.
Таблица 4
Расчет координат статической характеристики (скоростной) судового дизеля
Наименование |
параметра, |
|
|
Численное значение |
||||||||||
единица измерения |
|
|
|
|
|
|
|
|
|
|
|
|||
ИСХОДНЫЕ ДАННЫЕ |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Марка дизеля |
|
|
|
|
|
|
|
|
|
|
|
|||
Номинальная мощность, кВт |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Номинальная частота вращения |
|
|
|
|
|
|
|
|
|
|
||||
коленчатого вала nr, об/мин |
|
|
|
|
|
|
|
|
|
|
||||
|
|
РАСЧЕТНЫЕ ДАННЫЕ |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Относительная |
величина |
1.0 |
0.9 |
0.8 |
0.7 |
0.6 |
0.5 |
0.4 |
0.3 |
0.2 |
0.1 |
|||
номинальной |
частоты |
|
|
|
|
|
|
|
|
|
|
|||
вращения коленчатого вала, ωi |
|
|
|
|
|
|
|
|
|
|
||||
Долевая |
частота |
вращения |
|
|
|
|
|
|
|
|
|
|
||
коленчатого вала ni, об/мин |
|
|
|
|
|
|
|
|
|
|
||||
|
|
КООРДИНАТЫ ХАРАКТЕРИСТИКИ ДВИГАТЕЛЯ |
|
|
||||||||||
Внешней |
номинальной |
|
|
|
|
|
|
|
|
|
|
|||
мощности Ре |
|
|
|
|
|
|
|
|
|
|
|
|||
Относительная |
величина |
|
|
|
|
|
|
|
|
|
|
|||
номинальной мощности, Рei |
|
|
|
|
|
|
|
|
|
|
||||
|
|
КООРДИНАТЫ ВИНТОВОЙ ХАРАКТЕРИСТИКИ |
|
|
||||||||||
Частичной |
|
винтовой |
|
|
|
|
|
|
|
|
|
|
||
характеристики Рс |
|
|
|
|
|
|
|
|
|
|
|
|||
Относительная |
величина |
|
|
|
|
|
|
|
|
|
|
|||
частичной |
|
винтовой |
|
|
|
|
|
|
|
|
|
|
||
характеристики, Рсi |
|
|
|
|
|
|
|
|
|
|
|
3.8.4.Определение момента инерции пропульсивного комплекса.
20