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

Bystrov_Grigorye_Pershin_-_Sintez_raspredelennykh_regulyatorov

.pdf
Скачиваний:
14
Добавлен:
21.03.2016
Размер:
3.83 Mб
Скачать

m

 

 

 

N

2

~2

i

(1)

i

4

 

(x, y)=

 

 

Ei

(ψη

+ϕγ

)

 

Cη,γ ,ξ (s) Bη,γ ,ξ

i=1

 

η,γ =1

 

 

 

 

 

ξ =1

 

 

 

N

 

 

4

 

 

 

 

 

N

4

= Kη,γ

Cη,γ ,ξ (s) Bη,γ ,ξ (x, y)E0

∑ ∑ Cη,γ ,ξ (s) Bη,γ ,ξ (x, y).

η,γ =1

ξ =1

 

 

 

 

η,γ =1ξ =1

Преобразуя, придем к следующему результату:

 

 

N

4

 

 

 

 

 

 

m

2

~2

)

i

 

i

 

 

 

∑ ∑ Cη

γ ξ (s) Bη

 

 

 

 

 

(1)

 

 

γ ξ (x, y)

(ψη

+ϕγ

 

 

 

 

 

,

,

,

,

 

 

 

 

 

 

 

 

 

 

 

η,γ =1ξ =1

 

(η,γ =

 

).

i=1

 

 

 

 

 

 

 

где

 

 

= Kη,γ

E0 ,

 

 

 

 

 

 

 

 

 

 

 

η,γ

1, N

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

Kη,γ = 0 ,

Приравнивая выражение в фигурных скобках к нулю, получим систему уравнений

 

 

 

 

2

 

~2

)

i

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(η,γ =1, N ).

 

 

 

 

 

 

 

 

 

 

(1)

Kη,γ = 0,

 

 

 

 

 

(2.39)

Ei

(ψη

+ϕγ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Представим (2.39) в матричной форме:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ E K = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.40)

 

 

 

 

 

где

 

+ϕ1 );

L (ψ1

+ϕ1 )

 

 

(1)

 

 

; K

(ψ1

+ϕ1

)

 

(1)

 

;

 

(ψ1

i

 

i

m

m

 

 

2

~2

 

 

 

 

 

 

 

 

2

 

~2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

~2

 

 

 

 

2

~2

 

 

 

 

 

 

 

 

 

2

 

~2

 

i

 

 

i

 

 

 

 

 

 

 

 

 

 

 

2

~2

 

m

 

m

 

Ψ= (ψ1

+ϕ2 );

K (ψ1

+ϕ2 )

 

 

(1)

 

 

; K (ψ1

+ϕ2 )

 

(1)

 

;

;

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

2

~2

);

K (ψ

2

 

 

~2

)

i

(1)

i

; K

 

(ψ

2

~2

)

m

(1)

m

 

(ψ

N

+ϕ

N

N1

+ϕ

N

 

 

 

 

N

+ϕ

N

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E =

E

2

;

 

 

 

 

 

K =

 

 

K

2,1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

K

N,N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определитель матрицы Ψ сводится к определителю Вандермонда, который отличен от нуля. Следовательно, Матрица Ψ неособенная, а значит решение уравнения (2.40) существует, что и требовалось доказать.

Таким образом, используя уравнение (2.40), можно определить значения коэффициентов (Ei (i =1,m)) распределенного блока, преобразующего воздействие (2.36) в (2.37).

70

Рис. 2.14. График K.

Следствие из доказанного утверждения

Положим, что необходимо реализовать статический коэффициент усиления, график которого приведен на Рис. 2.14. При этом входное воздействие полагается заданным в виде (2.36). Определим спектр

значений координаты G при изменении η,γ =1, N

G1 (η =1, γ =1); G2 (η =1, γ = 2); KGm (η = N, γ = N).

Для вычисленных значений Gi (1,m) по графику определим значения

K(Gi ) (i =1,m) . Выберем E0 = K(G1) .

Подставляя вектор значений K в уравнение (2.40), определим вектор

коэффициентов

 

 

 

При этом матрица Ψ с использованием

 

 

 

 

 

 

Ei (i =1,m) .

обобщенной координаты G может быть представлена в виде

 

G (1);

K Gi

(1)i ;

K Gm

(1)m;

 

 

1

(1);

1

(1)i ;

1

 

 

 

G

2

K Gi

K Gm

(1)m;

 

 

 

 

 

2

 

 

 

2

 

.

 

 

 

 

 

 

K

 

 

 

K

 

 

 

 

 

 

 

 

i

 

i

;

m

(1)

m

 

Gm (1);

K Gm (1)

 

K Gm

;

Распределенный блок с вычисленными параметрами (Еi ) реализует коэффициент усиления K p (G) . В узлах Gi (i =1,m) значения K p (G) и K(G) (см. Рис. 2.14) совпадают.

2.2.2 Распределенный фильтр

Пусть требуется определить структуру и параметры распределенного блока, статический коэффициент усиления которого приведен на Рис. 2.15.

71

Рассматриваемый распределенный блок выделяет из входного воздействия совокупность пространственных мод, для которых выполняется условие G2 G G3 , т.е. по существу является

распределенным фильтром.

Рис. 2.15. График коэффициента усиления.

Структуру распределенного фильтра будем формировать из рассмотренных в п. 3.1 распределенных звеньев. Для реализации статического коэффициента в виде кривой 1 (см. Рис. 2.16) будем использовать пространственно-усилительное звено:

~

 

n

 

1

 

1

2

 

 

W

= E

 

1

 

 

 

 

.

 

 

 

n

1

1

 

n

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

Значения параметров E1 и n1 выбираются из следующих условий:

 

 

 

 

 

 

 

n1 1

 

G2

 

 

 

 

 

20 lg E =

 

,

20

lg E

+

 

=

 

,

1

n

n

 

2

1

 

 

 

1

 

 

 

 

 

где 1 и

 

 

 

 

 

 

1

1

 

 

 

 

 

2 - заданные числа (см. Рис. 2.16).

 

 

Определив из графиков (см. Рис. 2.16)

1= - 10 дБ;

 

2 =20 дБ;

lgG = 3,2, получим: E1 = 0,316 ; n1 = 51,8 .

Статические коэффициенты усиления, графики которых заданы в виде кривых 2 и 3 (см. Рис. 2.16), будем реализовывать с помощью устройства,

~

рассмотренного в п. 3.1 (см. Рис. 2.17), где передаточные функции W2 и

~

W3 имеют вид:

72

~

 

 

 

n

2

1

 

1

2

 

 

~

 

 

n

 

1

 

1

2

 

 

W

= E

 

 

 

 

 

 

 

 

 

;

W

= E

 

 

3

 

 

 

 

.

2

 

n

 

 

n

 

 

 

 

 

n

2

 

 

2

 

 

2

 

 

3

3

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

 

 

 

 

Рис. 2.16. Синтез фильтра.

Выберем произвольные точки на кривых 2 и 3 (см. Рис. 2.16)

A{lgG4 , 4}, B {lg G 5 , 5

}, (G4 > G2; G5 > G3) .

Используя статический коэффициент передачи усилительного звена, включенного в обратную соотношения между параметрами точек A и B:

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K2

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

4 = 20 lg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

n2

1

 

 

G4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

+ K2

E2

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2

 

 

 

 

n2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K3

 

 

 

 

 

 

 

 

5

= 20 lg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

~

 

 

 

 

 

n

1

 

G

 

 

 

 

 

 

1 + K

 

E

 

3

 

 

 

+

5

 

 

 

 

 

 

3

3 n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

3

 

 

К полученным

уравнениям

 

допишем

 

условие,

20 lg K(G) = 0 или

пространственносвязь, запишем

(2.41)

что при G 0

73

 

 

 

~

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

K2

 

 

 

 

 

=1,

 

 

K3

 

 

=1.

 

~

 

 

 

 

n

2

1

~

 

 

 

n

 

1

1 + K

2

E

2

 

 

 

 

 

1 + K

3

E

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2

 

 

3

 

n3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.42)

Рассматривая совместно (2.41) и (2.42), можно определить значения

~

параметров Ki , Ei , ni (i=2,3).

Полагая в соотношениях (2.42) n2 >>1 и

n3 >>1, преобразуем их к

виду

 

~

 

 

 

~

 

 

 

 

K

 

K

 

 

 

~ 2

=1,

~ 3

=1.

 

 

1 + K2 E2

 

1 + K3 E3

 

 

 

 

 

~

~

= 2 ,

Используя полученные соотношения, выберем значения K2

= K3

E2 = E3 = 0,5 .

 

 

Рис. 2.17. Структурная схема фильтра

 

 

 

 

Определив из графиков (см. Рис. 2.16)

4 = −10дб; 5 = −12дб;

G

 

 

~

 

,

4

=104,8 ; G =105,8 и подставив вместе с вычисленными значениями K

i

 

5

 

 

Ei

(i=2, 3) в уравнения (2.41), получим n2 = 9983,

n3 = 79246 .

 

 

Моделирование работы распределенного фильтра было осуществлено на ЭВМ. По результатам моделирования были определены значения коэффициентов усиления фильтра в зависимости от Gi (i=1,2 …) (см. Рис.

2.18).

74

Рис. 2.18. Результаты моделирования.

2.2.3 Распределенный регулятор прямого действия

Область применения распределенных регуляторов прямого действия может быть достаточно широкой. Эти регуляторы управляют «живым» сечением канала текущей жидкости и используются в системах подвода (отвода) тепловой энергии.

Положим, что требуется отвести энергию (q) от поверхности оболочки, показанной на Рис. 2.19, при этом интенсивность тепловой энергии меняется случайным образом и в различных точках поверхности оболочки она различна. Требуется разработать систему управления, обеспечивающую отвод тепловой энергии от оболочки.

Покроем оболочку с внутренней стороны специальными «секциями». В зависимости от интенсивности теплового потока (q) теплоноситель каждой секции нагревается (охлаждается), а следовательно регулятор (биметаллическая пластинка), помещенный в поток теплоносителя, увеличивает (уменьшает) живое сечение каналов, осуществляя регулирование расхода теплоносителя через заданную секцию. Блокпоглотитель энергии перерабатывает тепловую энергию в другие виды энергии (см. Рис. 2.19). Блок-насос поддерживает заданное значение

перепада давления P теплоносителя.

Как известно, биметаллическая пластинка (Б.П.) изготовлена из материалов А – активного (обладающим большим коэффициентом теплового расширения) и В – пассивного (обладающим малым коэффициентом расширения). Пусть Б.П. помещена в поток теплонесущей жидкости. При этом омывается потоком одна плоскость Б.П., а другая

75

находится в «тени» потока (см. Рис. 2.19). У обтекаемой потоком поверхности Б.П. формируется пограничный слой, толщина которого зависит от физических свойств теплонесущей жидкости. Изменение температурного поля на границе пограничного слоя приводит к изменению температуры внутри слоя, а следовательно изменяется температурное поле Б.П.

Пусть один конец Б.П. закреплен, а второй свободен. Изменение температурного поля Б.П. приводит к перемещению свободного конца пластинки на величину Ls (x,τ) (см. Рис. 2.19).

В качестве входного воздействия в регуляторе (Б.П.) выберем температуру теплонесущей жидкости на границе пограничного слоя.

Tср(xc , yc , zc ,τ) =T1(xc , yc , zc ,τ) ,

где T1(xc , yc , zc ,τ) - температурное поле пограничного слоя; xc , yc , zc - параметры пограничного слоя, толщина которого равна L.

Рис. 2.19. Распределенный регулятор прямого действия

76

В качестве функции выхода будет перемещение Ls (x,τ) . Структурная схема рассматриваемого регулятора приведена на Рис. 2.20.

Рис. 2.20. Структурная схема регулятора прямого действия.

Рассмотрим другие конструкции регуляторов прямого действия, управляющих «живым» сечением каналов текущих жидкостей с целью регулирования скорости подвода (отвода) тепловой энергии.

1. Конструкции на основе эластичных резервуаров (см. Рис. 2.21, а). Изменение температуры теплоносителя на величину T приводит к изменению объема эластичного резервуара на величину V, что в свою очередь приводит к изменению площади живого сечения на величину S. 2. Конструкции на основе биметаллических пластин (см. Рис. 2.21, б). Канал - 2, по которому движется теплоноситель-3, изготовлен из материала с низким коэффициентом линейного расширения. Снаружи канала нанесены пластины-1 активного материала, имеющего большой коэффициент линейного расширения. При изменении температуры теплоносителя происходит изменение профиля поперечного сечения

канала (изменение «живого» сечения канала).

Рис. 2.21. Конструкции регуляторов прямого действия.

В рассматриваемой конструкции активный материал может быть нанесен как на наружную поверхность канала, так и на внутреннюю поверхность.

Для построения динамических характеристик рассматриваемого регулятора составим его математическую модель.

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

77

Опишем тепловые процессы, протекающие в пограничном слое и в

Б.П.

При описании тепловых процессов будем полагать:

-боковые грани Б.П., за исключением поверхностей 1 и 2 (см. Рис. 2.19), не влияют на процессы, протекающие внутри Б.П.. Отметим, что рассматриваемые боковые грани имеют малую площадь боковой поверхности по сравнению с площадью поверхностей 1 и 2;

-тепловые процессы вблизи плоскости Б.П., находящейся в «тени» потока теплоносителя (см. Рис. 2.19, поверхность 2), не влияют на тепловые процессы Б.П.;

-Б.П. имеет прямоугольную форму.

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

T

 

 

2T

+

2T

+

2T

 

, (i=1,2,3),

i = a

 

 

i

 

i

 

i

T

i

 

x

2

 

y

2

 

z

2

 

 

 

 

 

 

 

 

 

 

 

(2.43)

0 < x < Lx ; 0 < y < Ly ; Zi+1 < z < Zi (z0 = 0) ,

где Ti (x, y, z,τ) - температурное поле i-ой среды (i =1 – пограничный слой, i=2 – материал А, i=3 – материал В); ai коэффициенты температуропроводности i-ой среды; Lx , Ly , Zi - заданные значения

(координаты границ раздела сред).

Граничные условия для системы уравнений (2.43), отражающие неразрывность тепловых полей и тепловых потоков на границах раздела сред, записывается в виде следующих соотношений

T1(x, y,0,τ) =Tср(x, y, zc ,τ), (zc = 0) ;

 

Ti (x, y, zi ,τ) =Ti+1(x, y, zi ,τ), (i =1,2).

(2.44)

 

 

 

λi

Ti (x, y, zi ,τ)

 

= λi+1

Ti+1 (x, y, zi ,τ)

, (i=1,2);

(2.45)

 

 

z

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 < x < Lx ; 0 < y < Ly ;

 

 

 

T2

(x, y, z,τ)

 

= 0

 

 

, z1 < z < z2; 0 < y < Ly ;

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

x=0, x=Lx

 

 

 

 

 

 

 

 

 

 

 

78

 

 

 

 

T3(x, y, z3,τ)

= 0 , 0 < x < Lx ; 0 < y < Ly ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

T3

(x, y, z,τ)

= 0

 

 

 

 

, z2 < z < z3; 0 < y < Ly ;

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

x=0, x=Lx

T3

(x, y, z3

,τ)

 

 

= 0

 

 

, 0 < x < Lx ; z2 < z < z3,

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

y=0, y=Ly

 

 

 

 

 

 

 

 

 

 

 

где λi - коэффициент теплопроводности i-ой среды (i =1,3) .

Для частотного анализа блока 1 (см. Рис. 2.20) положим, что входное воздействие Tcp (x, y, zcp ,τ) может быть представлено в виде

 

 

 

 

Tср(x, y, z,τ) =

exp( jϖτ) cos(ψη x) cos(ϕγ

y) ,

(2.48)

 

η,γ =1

 

 

где ψη =πη

, ϕγ

=πγ , (η,γ =

 

) , ω - круговая частота.

 

1,

 

Lx

 

Ly

 

 

(Вид разложения в ряд Фурье согласуется с граничными условиями).

Реакцию температурного поля на каждую составляющую ряда Фурье (2.48) будем искать в виде

Ti,η,γ (x, y, z,τ) = Hi,η,γ (z, iω) exp( jωτ) cos(ψη x) cos(ϕγ y) , (2.49)

(i =1,3; η,γ =1, ),

где Hi,η,γ (z, iω) - функция, подлежащая определению.

Подставляя (2.49) в (2.43) и преобразуя, получим

 

 

 

2 Hi,η,γ (z, jω)

 

jω H i,η,γ (z, jω) = ai

(ψη2

ϕγ2 ) Hi,η,γ (z, jω) +

 

 

, (2.50)

z 2

 

 

 

 

 

 

 

 

 

 

 

(i =1,2;η,γ =1, ) .

Преобразуя (2.50), придем к следующему результату:

79

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]