
книги из ГПНТБ / Кутьин Л.И. Автоматизация судовых дизельных и газотурбинных установок учебник
.pdfВ соответствии с двумя входными координатами в объекте можно рассматривать два канала воздействий или проводимости.
Координаты, изменяющиеся в цепочке передачи воздействий по каналу проводимости, могут быть сосредоточенными или распреде ленными. Поэтому все объекты разделяются на две группы: группу с распределенными параметрами и группу с сосредоточенными параметрами. Если в объекте параметр распределен (как, например, температура вдоль длины трубок в теплообменнике), и это распре деление учитывается, то движение такого объекта описывается диф ференциальным уравнением в частных производных.
Объекты судовых дизельных и газотурбинных установок обычно рассматриваются как объекты с сосредоточенными параметрами. Динамика линейной математической модели таких объектов описы вается обыкновенными дифференциальными уравнениями, порядок которых зависит от числа учитываемых емкостей (аккумуляторов).
Перейдем к математическому описанию динамики двигателя как объекта регулирования скорости вращения вала, полагая сначала, что запаздывание отсутствует.
Ранее, при оценке статических свойств, были получены аналити ческие зависимости для линейной модели объекта, связывающие приращения выходной и входных координат на установившихся режимах, когда АМА— ДМС0Пр = 0, а изменения координат не зависят от времени.
Такие уравнения были получены в предыдущем параграфе. Они могут описывать также движение линейной модели безынерцион ного объекта относительно некоторого равновесного состояния, принятого за исходное, либо движение инерционного объекта при квазистатическом изменении координат. При обычном изменении координат во времени динамика двигателя как одноемкостного объекта регулирования запишется уравнением вида
AMA- A M COnp^ J ^ , |
(1.19) |
где J — приведенный к оси вращения вала момент инерции всех вращающихся и присоединенных масс.
Двигатель внутреннего сгорания может рассматриваться как одноемкостный объект в случае, если он не имеет наддува или оборудован механической системой наддува, в которой нагнетатель жестко связан с коленчатым валом двигателя.
Так как |
= d |
то исходное уравнение (1.19) может быть |
переписано в приращениях относительно некоторого равновесного состояния, принятого за исходное:
у Д ^ І = ДМд- Д Л 1„пр. |
(1-20) |
38
Чтобы уравнение раскрывало зависимость между приращениями координат А/іт и АС, представим приращения моментов в соответ ствии с (1.10):
АМД— АМсопр = Кн А/іх -(- /Спод А©— Кс КС ^отв А<*>;
получим
J |
I (Котв ~ Кп0д) Дсо - КнАА, - Кс АС. |
(1.21) |
Перейдя к относительным величинам и координатам, разделим все слагаемые уравнения (1-21) на номинальную величину нагрузки и в соответствии с (1.15) получим
|
Таф "I" Z<p = КѴов — kcK |
( 1 |
. 2 2 |
) |
|
|
|||
где Та — |
У(0НОМ |
|
|
|
- -----постоянное время аккумулятора, или емкости. |
|
Для одноемкостного объекта оно является постоянным временем
объекта; |
z = k0TB— &под— коэффициент |
самовыравнивания, уже |
||
известный |
по статическим свойствам; kh и |
kc — известные по ста |
||
тическим |
свойствам угловые коэффициенты статических |
характе- |
||
|
|
„ |
„ |
До) |
ристик регуляторного и нагрузочного воздействии; ср = |
——---- |
относительное приращение (изменение) регулируемого параметра: Дсо = со — со0; со0 — значение регулируемого параметра на исход ном установившемся режиме, принятом за начало отсчета координат (за начало движения). Исходный установившийся режим может быть охарактеризован нагрузкой М°, при которой вводится возмущение.
Коэффициенты z, kh и kc определяются в точке исходного равно весного состояния, относительно которого рассматривается движе ние объекта.
Обратим внимание на величину Та:
J0)НОМ
уц н о м
(1.23)
Как видно, это время представляет величину, пропорциональную моменту инерции вращающихся масс J. Оно измеряет ту же емкость (инерционность) объекта, но в единицах времени. Если момент инерции — постоянная величина, то и время Тя постоянно и харак теризует инерционные качества объекта. В двигателях внутреннего сгорания, работающих на гребной винт, момент инерции J для каждого скоростного режима имеет свое значение вследствие влия ния различных по величине присоединенных масс воды. Однако часто его считают постоянной величиной, не зависящей от скорост ного режима судна.
Уравнение динамики линейной модели одноемкостного объекта можно представить и в иной форме записи — через коэффициенты усиления 6ц и kx. Для этого достаточно все члены уравнения (1.22)
39
разделить на коэффициент самовыравнивания г. Тогда уравнение примет вид
|
Т’оФ + |
Ф = КѴчл — |
(1.24) |
т |
— постоянная |
времени разгона аккумулятора, или |
|
здесь То = — |
|||
емкости (для одноемкостного объекта она является |
постоянной |
времени разгона объекта); kß и k%— коэффициенты усиления объекта по каждому каналу проводимости, определяемые статическими свойствами объекта в точке исходного установившегося режима.
Временные параметры Та и Т „ характеризуют инерционные свойства одноемкостного объекта. Однако время Та определяется только инерционными свойствами и для объекта является обязатель ной величиной. Оно постоянно, если постоянна величина емкости J. Постоянная времени разгона Т 0 имеет физический смысл лишь для объектов с положительным самовыравниванием. Она зависит не только от инерционных свойств, но и от величины коэффициента самовыравнивания г.
В частном случае, когда на исходном режиме z = 1 или когда с целью упрощения приближенно считают г = 1, постоянная вре мени разгона Т 0 равна постоянной времени объекта Та.
В общем случае для объекта с положительным самовыравнива нием существует соотношение Та= zT0 = const.
Поскольку коэффициент самовыравнивания у одного и того же объекта может быть различным на разных установившихся режимах, то и постоянная времени разгона Т 0 будет зависеть от начального режима, относительно которого рассматривается движение объекта. Для выбранного начального режима она постоянна.
Для раскрытия физической сущности понятия постоянной вре мени разгона объекта Т 0 обратимся к решению уравнения движения (1.24). Решение этого уравнения описывает кривую разгона (кривую переходной проводимости) и при скачкообразном возмущении, как известно, представляет экспоненту, уравнение которой записывается так:
Ф = (бцРоб — М-) (1 — е~ т° 0 , |
(1.25) |
|
или раздельно для канала |
регуляторного воздействия (при X = 0): |
|
ф = |
бцРоб (1 — е~ т° 0 |
(1.26) |
и для канала нагрузочного воздействия (при роб = |
0): |
|
Ф = - М < (і |
(1-27) |
Как видно, постоянная времени разгона есть время экспоненты Т 0. На рис. 13 представлены кривые разгона объекта, обладающего на исходном режиме некоторым положительным самовыравнива
40
нием г — z v Они построены для канала регуляторного воздействия и отличаются одна от другой величиной скачкообразного возмуще ния роб.
Для каждой кривой разгона существует свое конечное значение
выходной координаты сркон, которое наступает при t - - - > 0 0 |
и зависит |
от величины скачка. В соответствии с (1.26) имеем |
|
Фкон = lim ф = /фроб, |
(І.28а) |
t-> СО |
|
и для возмущения со стороны нагрузки (отвода) из (1.24) получим
Фкон = |
(1.286) |
Рис. 13. Кривые разгона двигателя как одноемкостного объекта с самовыравниванием.
Кривые разгона (или кривые переходной проводимости) могут быть построены и экспериментальным путем. По кривым разгона весьма просто определяется постоянная времени разгона объекта Т 0 и раскрывается его физическая сущность. Постоянная времени разгона объекта представляет время, в течение которого выходная координата объекта достигает своего конечного значения фкон при скачкообразном возмущении, изменяясь с постоянной скоростью,
равной скорости в начальный момент ф<=о- Скорость в начальный момент измеряется тангенсом угла наклона касательной, прове денной к экспоненте в начале координат. Эти касательные пересе каются с соответствующими линиями конечных значений выходной координаты в точках F х, F2, Fh, . . ., расположенных на общей вер тикали, которая удалена от начала координат на величину постоян ной времени разгона Т 0, как показано на рис. 13. Действительное же время разгона объекта равно бесконечности.
Так как постоянная времени разгона не зависит от величины возмущения, то она определяется по любой из кривых разгона.
41
При этом скачкообразное возмущение в равной мере может быть и со стороны нагрузки Я, а не только со стороны подвода роб.
С точки зрения управления и регулирования объекта определя ющим динамическим показателем служит скорость изменения регу
лируемого параметра <р, для которой из (1.26) получим следующее выражение:
|
Ф = |
£|#об - 4 - е т° 1- |
(1.29) |
|||
|
kh |
1О |
|
|
|
|
Учитывая, что |
а T 0z = |
Та, |
будем |
иметь |
||
kß =-j-> |
||||||
|
Ф = |
1 |
- у - |
1 |
(1.30) |
|
|
*№>б^г-е |
а |
• |
Как следует из (1.30), скорость изменения регулируемого пара метра для данного объекта обратно пропорциональна постоянной времени Та, прямо пропорциональна величине скачкообразного возмущения роб и угловому коэффициенту kh статической харак теристики регуляторного воздействия. С течением времени она уменьшается и стремится к нулю. Темп уменьшения скорости зави сит от коэффициента самовыравнивания; он тем быстрее, чем больше
самовыравнивание. В начальный момент скорость ф<=0 не зависит от самовыравнивания и равна
ф<=0 = £/іРоб ~ у |
• |
(1.31) |
1 |
а |
|
Если на исходном режиме объект не обладает самовыравниванием (г = 0), то его динамика описывается лишь уравнением вида (1.22), так как постоянная времени разгона Т 0 для такого объекта теряет физический смысл, превращаясь в бесконечность.
Уравнение динамики нейтрального объекта запишется так:
Ta(f = khiio6 — kcl. |
(1.32) |
Решение этого уравнения, представляющее линию разгона объекта при скачкообразном возмущении, есть уравнение прямой, а не экспоненты:
ф = khno6-^r-i. |
(1.33) |
1 а
Скорость изменения регулируемого параметра в этом случае постоянна во времени, так как определяется выражением
Ф = 6/гИоб-ж-- |
(1-34) |
* а |
|
Интересно отметить, что эта скорость равна скорости в началь ный момент для объекта, обладающего самовыравниванием, что видно из сопоставления выражений (1.31) и (1.34). Поэтому при
42
веденные на рис. 14 линии разгона нейтрального объекта будут совпадать с касательными, построенными на рис. 13 для объекта
сположительным самовыравниванием, если одинаковы возмущение
иинерционность объектов.
По линиям разгона нейтрального объекта можно определить его постоянную времени Та. Для этого необходимо на оси регулируе
мого параметра отложить вели- |
|
<р |
|
|
|
||||||
чину, |
равную &/гЦоб‘> тогда на |
|
|
/ |
|
||||||
соответствующих |
прямых |
раз- |
|
|
|
|
|||||
гона будут получены точки Fly |
|
|
|
/ ^ |
|
||||||
|
|
Fj А/ |
|
||||||||
F2. . ., |
абсциссы которых опре- |
|
|
|
7 |
f ' |
|||||
деляют время Та. Действи- |
|
|
к |
||||||||
тельно, |
из (1.33) |
следует, |
что |
|
|
£---.У. |
|
||||
при <р = |
kh[io6 имеем t = |
Та. |
~ |
__ |
/ / |
|
|
||||
В том случае, |
когда |
на ис- |
|
|
|||||||
ходном режиме объект неустой- |
а. |
о |
|
|
|
||||||
чив, его движение может быть |
|
а. |
|
|
|
||||||
|
А? ІО |
|
|
|
|||||||
описано, |
как и для |
нейтраль- |
|
£ |
|
|
|
||||
ного объекта, лишь уравнением |
|
|
|
|
|
||||||
вида (1.22): |
|
|
|
|
|
|
|
|
|
||
Таф + |
2ф = kh\lo6— ксК, |
|
(1.35) |
|
0 . |
\ _ |
|
t,c |
|||
так как |
при г < |
0 |
постоянная |
|
|
|
|
|
|||
времени |
разгона |
Т0 Приобре- |
Рис. 14. Кривые разгона |
двигателя как |
|||||||
тает отрицательный |
знак |
и те- |
|
нейтрального одноемкостного объекта. |
ряет физический смысл.
Решение этого уравнения, включающее абсолютное значение
коэффициента |
самовыравнивания, запишется |
так: |
||
|
|
|
|
(1.36) |
Конечное |
значение |
параметра |
|
|
|
|
Фкон = lim tp = |
оо. |
(1.37) |
|
|
І ->СО |
|
|
Скорость |
изменения |
параметра |
|
|
|
|
1 |
Т ~ 1 |
С1'38) |
|
|
<Р= £лРоб~т~е |
а . |
|
|
|
*а |
|
|
и в начальный момент она составляет |
|
|
||
|
|
Фі=о = ^ Р о б - ж - • |
(1-39) |
|
|
|
|
1 а |
|
С течением времени скорость увеличивается и стремится к беско нечности. Темп возрастания скорости зависит от значения коэффи циента самовыравнивания г. темп тем интенсивнее, чем больше абсолютная величина г.
43
Как можно заметить, начальная скорость изменения регулируе
мого параметра <рі=0 не зависит от свойства самовыравнивания объ екта и во всех случаях определяется одинаковыми выражениями: (1.31), (1.34) и (1.39).
Таким образом, для одного и того же объекта при одинаковом скачкообразном возмущении можно получить принципиально различ ные кривые разгона в зависимости от выбранного исходного устано вившегося режима, относительно которого рассматривается движе ние. Это обусловлено свойством самовыравнивания на исходном ре жиме, т. е. величиной коэффициента г. Однако начальная скорость
изменения ф<=0 будет одинакова.
В практике теоретического анализа динамических свойств одно емкостных объектов их уравнения динамики представляют в виде, соответствующем уравнениям типовых элементарных динамических
звеньев, и записывают вместо (1.22) так: |
|
|
* |
7Ѵр + гф = роб — Я*, |
(1.40) |
. * |
пропорциональ |
|
где Цоб = |
kh\iоб и А = kck являются величинами, |
ными истинному изменению координат регуляторного и нагрузоч
ного воздействия, а вместо (1.24): |
|
|
|
||
|
|
ГоФ-f ф = |
цГб - |
Г , |
(1.41) |
где Цоб |
&цЦоб и Я |
— Я?,Я. |
такой |
формы записи |
уравнения |
Непосредственное |
применение |
динамики объекта совместно с уравнением движения регулятора при анализе свойств системы регулирования в целом может привести к заметным неточностям, если не учитывать отличия истинного из менения входных координат роб и Я от пропорциональных им вели
чин роб, X* или роб и Я • Простая замена их означает, что коэффи циенты kh и kc или и kx приняты равными единице. Однако в пер вом приближении часто допускают такую замену и записывают урав нение динамики объекта так:
Таф-f гф = |
роб — Я, |
(1.42) |
или |
|
|
7\>ф + ф = |
Роб— А,- |
(1.43) |
При описании динамических свойств дизеля до сих пор мы не учитывали явления запаздывания в передаче воздействий от входной координаты регулирующего органа (ТНВД) р к входной коорди нате в объект роб. Мы полагали, что изменение координаты роб во времени точно соответствует изменению р, а именно: роб (t) = р (t), и поэтому в уравнении динамики под роб понималось воздействие, приложенное к регулирующему органу.
Если учесть присущее дизелю запаздывание в передаче воздей ствия р к роб, то связь между ними будет записана так:
Еоб (0 = Е (* — Тзап). |
(1.44) |
где т3ап — время запаздывания, с.
44
Координата роб может принимать следующие значения: при
t <С Т-зап Роб |
П рИ |
I |
Тзап р од |
р,. |
В общем |
случае |
любое динамическое звено с запаздыванием |
можно разбить на два: обыкновенное линейное звено и элемент чи
стого |
запаздывания. |
Уравнение |
ф |
|
||||
динамики в этом случае записы- |
‘ 1 |
Z=Z) |
||||||
вается так (для канала регулятор- ^ |
[ |
|||||||
ного |
воздействия): |
|
|
|
F |
|
||
|
Т’аф + |
= |
kh\lo6, |
J |
|
|
||
|
(^об = |
ПК— Тзап)- |
I |
|
|
|||
Кривая |
разгона |
объекта при |
|
|
||||
наличии запаздывания будет иметь |
|
|
||||||
вид, представленный на рис. 15. |
|
|
||||||
Перейдем к описанию |
динами |
0 ТзапІ Т0 |
|
|||||
ческих свойств линейной модели |
и |
|||||||
двигателя, |
когда |
он |
оборудован |
Рис. 15. Кривая разгона двигателя |
||||
системой |
наддува. |
Из множества |
как одноемкостного объекта с запазды |
|||||
возможных схем наддува |
двигате |
|
ванием. |
|||||
лей, |
которые |
изучаются |
в специ |
|
|
альном курсе ДВС, рассмотрим наддув со свободно вращающимся нерегулируемым турбонагнетателем в соответствии со схемой, изобра женной на рис. 16. Из схемы можно заметить, что в составе объекта регулирования имеются два основных аккумулятора энергии: один —
О х л а д и т е л ь В о з д у х
Рис. 16. Принципиальная схема двигателя с газотурбинным наддувом.
вращающиеся и присоединенные к ним массы двигателя и нагрузоч ного агрегата; второй — вращающиеся массы турбокомпрессора. Кроме отмеченных основных имеются и другие аккумуляторы: ре сивер воздуха с трубопроводами, где аккумулируется масса воздуха и изменяется при этом его давление; выпускной коллектор, где акку мулируется газ; объем газа в турбине, являющейся аккумулятором
45
потенциальной энергии; стенки камеры турбины и выпускного кол лектора, в металле которых аккумулируется тепловая энергия. Влия нием этих аккумуляторов будем пренебрегать, пренебрежем также запаздыванием. Параметры воздушного ресивера будем считать со средоточенными и, так как его емкостью пренебрегаем, то давление наддува перед двигателем ps можно считать равным давлению за компрессором рк в любой момент времени.
Выведем уравнения динамики линейной модели двигателя с газо турбинным наддувом как объекта регулирования угловой скорости вала.
Д вигатель — винт. Исходное уравнение движения в прираще ниях для двигателя, работающего на винт, запишется аналогично
уравнению для двигателя без наддува (1.22): J d — AMд—
ЛМСОпр.
Однако в рассматриваемом случае момент, развиваемый двига телем, будет функцией времени t и трех величин:
Мд = Мд(/гТ) со, ps). |
(1.46) |
Функциональная зависимость для момента сопротивления не ме няется (1.4): Л4С0ПР = Л4сопр (С, со).
Выражая приращения моментов ДМд и ДЛ4С0пр в линейном при ближении через приращения определяющих их координат, получаем
или, переходя к относительным величинам,
Гаф + zfp = /гЛроб — k c X + |
р , |
(1.47) |
где Та, z, kh, kc определяются так же, как для двигателя без над дува; kpÄ— относительный угловой коэффициент статической харак теристики Мл = Л4д (ps) при постоянных со и С:
kPH |
(1.48) |
Фр— относительное изменение во времени приращения давления над дува:
Арз Р*-Р*
Турбокомпрессор. Запишем исходное уравнение динамики турбо компрессора в приращениях:
4 ^ Р ± = АМТ- А М К> |
(1.49) |
46
где Дсок — приращение угловой скорости вала турбокомпрессора:
Дсок = (ок— Gv, / к — момент инерции вращающихся масс турбо компрессора; АМт — приращение момента, развиваемого газовой турбиной; АЛ4К— приращение момента сопротивления компрессора.
Момент, развиваемый газовой турбиной, является функцией ко
личества газа Gr, которое определяется подачей топлива |
и воздуха |
в двигатель, а также угловой скоростью ротора турбины |
сок. Коли |
чество подведенного в двигатель топлива можно представить завися
щим от координаты регулирующего органа |
hT и угловой скорости |
|
вала двигателя, поэтому |
|
|
Мт = Мт(hT, со, |
сок). |
(1.50) |
Момент сопротивления на валу компрессора |
|
|
Мк = Мк (сок, |
ps). |
(1.51) |
Представляя приращения моментов АМТи АМк в линейном при ближении через приращения координат, определяющих их, уравне
ние (1.49) перепишем так: |
|
|
j |
d (Дсок) |
|
JK |
dt |
|
(.V<?<ок /о Асок— |
(1.52) |
После перехода к относительным величинам и деления всех сла гаемых на Л4?ом получим
|
|
7 > К |
~Т ^ к Ф к |
*АкНч>6 |
|
^ ю к Ф |
^ р к Ф р - |
|
(1.53) |
||||
|
|
|
|
|
|
|
|
|
|
|
I |
, ч Н С |
|
Здесь: Тк — постоянная |
времени |
компрессора 1Т] |
JA |
||||||||||
M"ot |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
фк — относительное |
изменение угловой скорости вала турбокомпрес- |
||||||||||||
|
|
<а |
|
|
|
|
|
|
|
|
|
||
сора |
срк |
— со“ |
; zK— коэффициент |
самовыравнивания |
турбо |
||||||||
|
|
||||||||||||
компрессора • |
z = |
(JM*\ |
< Ы ( дМг ) |
; khK— относитель |
|||||||||
|
|
к |
V |
<Эй)к /0 м н о м |
\ |
д с о к |
/ о |
Ml |
|
|
|||
ный |
угловой |
коэффициент статической |
характеристики |
регулятор |
|||||||||
ного |
воздействия |
Л4Т= Мт(hT) при постоянных значениях |
со и сок: |
||||||||||
= |
( дМ \ |
hH0M |
*«* — относительный |
угловой коэффициент |
|||||||||
( “0*7 )0 |
- ф я г ; |
||||||||||||
|
|
м; |
|
|
Мт— Мт(со) |
|
|
|
|||||
статической характеристики |
при постоянных |
значе- |
|||||||||||
|
|
|
............................чном |
|
|
|
|
|
|
|
|||
ниях |
hT и сок: kaK ■ |
/ |
дМт\ |
— - ; |
|
knK— относительный |
угловой |
||||||
\ |
дел )с |
|
|||||||||||
|
|
|
|
/о Л1"ом |
|
Р |
|
|
|
Мк (ps) |
|||
коэффициент статической характеристики компрессора Мк = |
|||||||||||||
|
|
|
|
|
|
|
( |
-4 W |
,.Н О М |
|
|
||
|
|
|
|
|
|
|
дМк \ |
|
|
|
■*57)0 -^ s r*
47