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

Курс_АСУ_2013

.pdf
Скачиваний:
93
Добавлен:
15.02.2015
Размер:
1.27 Mб
Скачать

3.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= ∆λir,,

где ω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