Vì điều kiện thời gian, nên luận văn mới dừng lại ở mức độ thành công về mặt
điều khiển theo phương pháp LQ cho tốc độ Turbine thuỷ điện tốt hơn khi dùng
phương pháp điều khiển truyền thống (bộ điều khiển PID), trong đó hệ được xét ở lân
cận điểm làm việc xác lập. Việc kết hợp kỹ thuật chuyển mạch hay Gain-Scheduling để
có thể điều khiển đối tượng trên toàn dải tải, nghiên cứu xây dựng chương trình điều
khiển để thử nghiệm trên thiết bị thuỷ lực là vấn đề mà tác giả dự định tiếp tục phát
triển nghiên cứu trong thời gian tới.
101 trang |
Chia sẻ: lylyngoc | Lượt xem: 2323 | Lượt tải: 2
Bạn đang xem trước 20 trang tài liệu Luận văn Điều khiển turbine thuỷ điện, để xem tài liệu hoàn chỉnh bạn click vào nút DOWNLOAD ở trên
Dầu bổ xung cho buồng hơi được lấy từ thùng xả. Từ đây máy bơm dầu sẽ bơm
dầu qua van một chiều đi vào nồi hơi. Khi áp lực dầu quá tăng so với giá trị cho
phép thì van an toàn mở không cho phép dầu vào nồi hơi nữa. Thiết bị dầu áp lực
thường có hai máy bơm, một chiếc làm ở chế độ dự phòng.
2.5.6. Tính toán thông số chính của điều tốc Turbine.
2.5.6.1.Các thông số đặc trưng.
- Công suất Turbine:
N = 9,81
QH
N = 0,9.1489,13.30 = 41100 (KW) (2-13)
- Số vòng quay đặc trưng:
ns =
4 HH
Nn
(2-14)
ns =
4 3030
41100.4,136
= 95,371
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
54
Theo Fapurt với ns = 95,371 tra biểu đồ quan hệ ns với kiểu Turbine và cột áp,
ta được kiểu Turbine thiết kế là loại Turbine KapLan với các thông số sau:
n’Itư
= 58 (v/p); n’Itt = 62 (v/p)
Q’Imin = 0,1 (m
3/s); Q’Imax = 0,18 (m
3
/s)
Từ các thông số trên ta tính được đường kính bánh xe công tác của Turbine theo
công thức:
D
2
1
=
..'.81,9 2
3
ttI
tt
HQ
N
(2-15)
D
2
1
=
9316,0.30.14,0.81,9
41100
2
3
= 2,37
D1 = 1,54 (m)
Chọn D1 = 1,5 (m)
Từ loại Turbine đã chọn có được các thông số hình học của Turbine thực như
sau:
b0/D0 = 0,1 b0 = 0,1. D0 =0,1.1,5 = 0,15 (m) (2-16)
D0/D0 = 1,24 D0 = 1,24.D1 = 1,24.1,5 = 1,86 (m) (2-17)
Chọn D0 = 1,8 (m)
2.5.6.2. Tính toán cánh hướng.
Chọn cánh hướng:
Vì dòng chất lỏng đi qua cánh hướng không trao đổi năng lượng do vậy nó tuân
theo quy luật dòng chế
Cu.r = const (2 - 18)
Chọn cánh hướng theo tài liệu
Đƣờng kính bánh công tác Số cánh hƣớng Z0
nhỏ hơn 2,5 m 16
2,5m – 6,5m 24
Lớn hơn 6,5 m 32
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
55
Với D1 = 1,5 (m) chọn z0 =16 cánh
Bước cánh hướng:
T =
16
5,114,3
0
0
z
D
= 0,294 (m) (2 - 19)
Để đảm bảo đóng kín hoàn toàn cánh hướng khi dừng Turbine chọn chiều dìa
cánh hướng là L = 600mm
Chọn loại frofin đối xứng có các kích thước như sau:
Các kích thƣớc Mẫu Thực
L0 172,6 750
R 11,5 50,0
28 121,7
L1 83,6 363,3
L2 89 64,3
C1 14,8 64,3
C2 14,8 64,3
C3 14,8 64,3
C4 14,8 64,3
C5 14,8 64,3
D1 1,5 6,5
D2 1,5 6,5
E 16,2 70,4
F 17,9 77,8
G 18,6 80,8
H 14,7 63,9
17,9 77,8
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
56
2.5.6.3. Tính lực đóng mở cánh hướng.
Khi cánh hướng đóng hoàn toàn, ta có:
- Lực tác dụng lên cánh hướng:
F’0 =
0
00
Z
bD
H’ (2-20)
Trong đó:
H’: Cột áp khi đóng cánh hướng có tính đến hệ số tăng áp.
H’ = (1+
)Hmax
= 1,3 - 1,5 (2-21)
- Mômen tác dụng lên cánh hướng:
M’0 =
00
2
0
0 )( nb
Z
D
H’ (2-22)
n0 – độ lệch của cánh hướng, n0 = 0,05 (m)
M’0 =
13658.5,2.1000.05,0.15,0)
16
8,1.14,3
( 2
(kGm)
- Lực lớn nhất tác dụng lên cánh hướng (tính đến ma sát thuỷ lực và ma sát
tại các ổ đỡ):
F0 = (1,2-1,3)F’0 = (9220-9988) (kG) (2-23)
Chọn F0 = 9600(kG)
- Mômen lớn nhất tác dụng lên cánh hướng (tính đến ma sát thuỷ lực và ma
sát tại các ổ đỡ):
M0 = (1,2 -1,3)M’0 = (163 - 177) (kGm)
Chọn M0 = 170 (kGm)
- Từ M0 ta tính được lực cần thiết của động cơ Servo tác dụng lên cánh
hướng khi đóng hoàn toàn:
Fc = Z0. K1. K2.K3.
pl
M 0
(2-24)
- Ta xác định được các thông số theo tài liệu “Turbine thuỷ lực – các thiết bị
thuỷ lực và công trình trạm thuỷ điện”.
D = 3500 (mm)
D0 = 3800 (mm)
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
57
Dc = 3300 (mm)
lc = 200 (mm)
lp = 200 (mm)
Z0 = 16 cánh
= 70
o
= 80
o
; K1 = 1,2
K2 =
yp
cc
Rl
lR
.
.
= 0,897
K3 =
sin
cos
= 0,35
FC = 16.1,2.0,897.0,35.
2,0
170
= 5124 (kG)
Chọn FC = 5000 (kG)
2.5.6.4. Tính toán thông số động cơ secvô và cụm bơm.
a. Tính động cơ tiếp lực:
Ta có: dH =
max
1
0
1
1
H
D
b
D
ZH
(2-25)
Trong đó: dH là đường kính của xilanh.
ZH : số xilanh của động cơ.
: hệ số phụ thuộc vào số cánh hướng.
Số cánh hƣớng Z0 16 24 32
Hệ số
0,034 0,03 0,028
dH = 0,5.0,034.1,5.
58.1,0
= 0,0614 (m)
Chọn dH = 60 (mm) theo tiêu chuẩn
SHmax = (1,4 1,6)a0max (2-26)
SHmax: độ dài lớn nhất của ĐCTL.
a0max: độ mở của cánh hướng dòng.
Ta có công thức:
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
58
a0max =
00
max
vZb
Q
(2-27)
Với v =
tb
bx
tt
F
Q
(2-28)
v =
25,0.
806,10
=13,77 (m/s)
Suy ra: a0max =
16.77,13.15,0
139,11
= 0,34 (m) = 340 (mm)
SHmax = (0,476 0,544) (m)
b. Chọn thiết bị áp lực dầu.
Thể tích ĐCTL: VH =
4
d
2
H
.ZH.
H
.SHmax (2-29)
Với
H
: hệ số (0,7
1)
VH =
4
14,3
.0,3
2
.2.1.0,5 = 0,07 (m
3
)
+ Lưu lượng dầu có áp đi vào van trượt chính của máy điều tốc
Qvt =
10
70
s
H
T
V
= 7 (l/s) (2-30)
+ Đường kính dẫn ống từ van trượt đến ĐCT
dvt =
d
vt
V
Q
4
(2-31)
Vd = 5 (m/s): vận tốc dầu có áp trong ống dẫn dầu.
dvt =
5.14,3
10.7.4 3 = 0,042 (mm)
Chọn dvt = 0,05 (mm)
Kích thước của thiết bị dầu áp lực phụ thuộc vào thể tích nồi hơi, dung tích
của nó phải đủ để đóng BPHN trong điều kiện bất lợi khi áp suất nồi hơi thấp
hơn áp suất định mức từ 35%
40% đồng thời đủ để bổ xung lượng dầu rò rỉ.
Thiết bị dầu có áp tính theo dung tích cần thiết của nồi hơi V, khi đó dung tích
được tính như sau:
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
59
V = (18
20 )VH + (4 5)Vbx (cho turbine cánh quay)
V = (18
20)VH + (9 10)Vbxk (cho turbine tâm trục có van xả không tải)
V = (18
20)VH (cho turbine tâm trục)
Trong đó:
VH : tổng dung tích ĐCTL của BPHN.
Vbx : tổng dung tích các ĐCTL của BXCT.
Vbxk: tổng dung tích các ĐCTL van xả không tải.
V = (1,26
1,4) (m
3
)
Chọn các thông số kích thước cơ bản của TBDAL như sau:
Thiết bị áp lực: TBAL 1,6
Nồi hơi: thể tích nồi hơi là 1,6 m3
kích thước nồi hơi H = 2460 (mm)
D = 1040 (mm)
h1 = 385 (mm)
trọng lượng: 1,7 tấn
Thùng dầu: thế tích: 2,5 (m3)
Kích thước: H = 990 (mm)
L = 1600 (mm)
B = 1800 (mm)
A = 1250 (mm)
Trọng lượng: 2,8 tấn
Máy bơm dầu:
Lưu lượng: 3,5 (l/s)
Công suất: 19 (KW)
Số vòng quay 2930 (v/p)
Đường kính ống 50 (mm)
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
60
CHƯƠNG III: TỔNG HỢP HỆ THỐNG ĐIỀU TỐC TURBINE
3.1. Đặt vấn đề.
Hệ thống điều tốc Turbine là một hệ thống điều khiển đa biến có mối liên hệ
phức tạp giữa điện cơ và điện từ. Khi nghiên cứu bài toán với những điều kiện ban
đầu, thông qua khai thác cơ sở và đo đạc thực tế ta thấy hệ mang tính phi tuyến
mạnh. Trước hết nghiên cứu quá trình điện từ tác động lên quá trình điện cơ ta thấy:
Với mỗi vùng tải (tải từ phía điện) sai số tốc độ khác nhau là rất lớn. Sai số tốc độ
tại vùng tải nhỏ là nhỏ nhất (với điều tốc dải rộng), vùng tải nằm trong mức 85% -
95% tải định mức là tương đối nhỏ, vùng tải định mức thì sai số lớn hơn và vùng tải
lớn nhất Nmax thì sai số lớn nhất. Như vậy, các quan hệ điện từ lên điện cơ không
hoàn toán tuyến tính. Khi nghiên cứu quá trình điện cơ lên điện từ cũng vậy với các
giá trị tốc độ khác nhau (trong lân cận cho phép sai số tốc độ). Khi tốc độ nằm dưới
tốc độ định mức , sự sai khác điện áp là rất lớn. Trong khi đó cũng với sự sai khác
này, nếu tốc độ nằm ở phía dưới vùng tốc độ trên vùng tốc độ định mức thì sai khác
là nhỏ hơn nhiều. Như vậy hệ mang tính chất phi tuyến mạnh.
Về sự tương tác dây chuyền giữa các quá trình ta cũng thấy rõ. Giả sử nghiên
cứu quá trình nhận tải của máy phát, ngay tại thời điểm ban đầu điện áp ngoài việc
bị sụt mạnh do điện kháng đồng bộ, do phản ứng gây nên. còn do lúc đó tốc độ sơ
cấp bị sụt đột ngột vì hệ thống bơm nhiên liệu chưa kịp phản ứng. Điện áp tại thời
điểm đó có độ sụt áp rất lớn dẫn đến hệ thống tự động điều chỉnh điện áp phải hoạt
động tạo tín hiệu cường kích để tăng nhanh điện áp hao hụt. Cứ như vậy sau một
thời gian quá độ đủ để hệ thống tự ổn định.
Sự tác động của hệ thống điện từ lên điện cơ đó là do mômen máy phát,
mômen này chính là mômen phụ tải tác động lên hệ cơ. Mặt khác hệ tác động lên
điện từ chính là sự tác động của tốc độ chính, tốc độc cũng sẽ làm thay đổi điện áp
của máy phát.
Như vậy sự tương tác giữa quá trình điện cơ và điện từ là một quá trình phức
tạp, sự ràng buộc mang tính chất dây chuyền.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
61
Một mặt do mối quan hệ giữa các thông số mang tính phi tuyến và mặt khác
các thông số luôn thay đổi, chính vì vậy sẽ rất khó khăn khi xác định được bộ điều
khiển theo phương pháp truyền thống.
Nội dung của chương này sẽ mặc định các thông số ổn định, các quan hệ phi
tuyến sẽ được tuyến tính hoá chính vì vậy sẽ dễ dàng cho việc tổng hợp và thiết lập
bộ điều chỉnh hệ thống.
3.2. Mô hình toán học.
3.2.1. Khâu turbine.
Ta có : Md – Mc = Jd /dt (3 -1)
J là mô men quán tính của Rôto tổ máy thuỷ lực.
tốc độ góc của Rôto tổ máy.
Md là mô men chuyển động của turbine.
Mc là Rôto chuyển động của rôto tổ máy.
t là thời gian.
Md = Nt / =
QH
/
(3- 2)
Md(p) – Mc(p) = J.p. (p)
H.
.Q(p)/
- Mc = Jp (p)
Kn.Q(p)- Mc = J.p. (p)
Mặt khác theo công thức (2-3) ta có:
Q =
22
2
01
01
2
2
cot2
cot.
gF
U
br
gU
U
gHlt
(3-3)
Q = f(
0
)= KQ.
0
= dQ/d
0
(
0
=
dm
)
Turbine là khâu quán tính vì vậy ta có thể mô hình toán học khâu Turbine.
Mc
KQ
KnQ
Tnp + 1 Jp
1M
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
62
Hình 3-1. Mô hình toán học khối Turbine.
Ta có: Kn=
H
/
đm =
27,14
30.81,9
= 20,62 (3-4)
J = GD
2
/4g =
(kGm) 1050.103
4,0.27,2.10.4
7,0.17.10.40.182
..4
..182
2
3
0
2
0
1 ng
TN S
(3-5)
Q =
CgB
A
cot.
(3-6)
A =
520 14,27.1,52
27,14
30.10.9,0
.
.. 2
2 r
Hg
B =
02
1
b
=
0,106
5,1.14,3.2
1
C =
07,075cot.
25,6
5,1
cot 02
2
2 gg
F
r
520
)106,075cot.106,0(77sin
106,0.51
)cot..(sin
.
200222
gCgB
BA
d
Q
Vậy KQ = 520.
3.2.2. Khâu khuếch đại.
Hình 3-2. Hình vẽ động cơ thuỷ lực.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
63
Động cơ thuỷ lực có công thức:
x
lRA
PP
dt
d es
.
....2
(3-7)
p.
(p) = Ks. p(x) (p) =
)(. px
p
K s
Tính ra Ks = 5
Để làm dịch chuyển kim van trượt một đoạn là x cần phải có hệ thống điện khí
nén khi tác động điện áp U(p) sẽ làm dịch chuyển đoạn x. Ta có:
x(p) =
)(
1.
pU
pT
K
x
x
(3-8)
Như vậy hàm truyền đạt động cơ thuỷ lực
Ks
pTxp + 1
(p)U(p) Kx X(p)
Hình 3-3. Hàm truyền khối khuếch đại.
Trong đó Kx = 0,2/10 = 0,02
Tx = 0,1
3.2.3. Các khâu đo.
3.2.3.1. Khâu đo vị trí.
Sử dụng cảm biến đo vị trí là công tắc hành trình. Độ lớn mở cánh hướng sẽ tỉ lệ
với điện áp đầu ra của khâu đo. Ta có hàm truyền khâu đo:
Fu(p) =
1. pT
K
(3-9)
K
= 10/0,3 = 33; T
= 0,1
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
64
3.2.3.2. Khâu đo tốc độ.
Là cảm biến quang bao gồm một đầu phát và một đầu thu, khi phát thì ánh sáng
sẽ phản xạ từ trục Turbine đến đầu thu. Tốc độ càng cao thì tần số lượng xung điện
áp phát ra càng nhiều sau đó các xung này vào đầu vào của bộ biến đổi D/A chuyển
về tín hiệu điện áp, điện áp càng cao ứng với tốc độ càng lớn. Ta có hàm truyền:
1
)(
pT
K
pF
(3 -10)
K
= 10/50 = 0,2
T
= 0,001
3.3. Tổng hợp hệ thống.
3.3.1. Tổng hợp mạch vòng vị trí.
Ks
pTxp + 1
U® Kx x
K
Tp + 1
R
U
Hình 3-4. Sơ đồ cấu trúc mạch vòng vị trí.
Với các tham số T
= 0,001; Ks = 5; Kx = 0,02; Tx = 0,005; K
= 33;
Ta có: S01(p) =
)1)(1()1)(1(
.. 1
pTpTp
K
TpTp
KKK
xx
sx
(3- 11)
Do Tx và T
là những hằng số thời gian bé vì vậy áp dụng theo tiêu chuẩn
modul tối ưu thì bộ điều chỉnh R
là một khấu khuếch đại R
= k. Ta cũng có:
S01(p) =
)1(
1
pTp
K
, trong đó T
= Tx + T
Hàm truyền đạt hệ kín:
F1k =
1.
.
1
.
1
.
.
.1
11
1
2
1
01
01
p
kkkk
TkkppT
kk
SR
SR
(3-12)
Mặt khác áp dụng theo tiêu chuẩn modul đối xứng ta có:
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
65
F1k=
12..2
1
22 TppT
(3-13)
Suy ra:
1.
1
kk
= 2.T; T
= T
2
1
2
.
T
kk
T
k =
1.2
1
KT
Tính toán:
T
= Tx + T
= 0,005+ 0,001 = 0,006
K1 = Ks.Kx.K
= 0,002.5.37 = 0,37
K =
37,0.006,0.2
1
= 252,5
3.3.2. Tổng hợp mạch vòng điều chỉnh tốc độ.
Ta có hàm truyền đạt kín mạch vòng vị trí
F1k =
12..2
1
22 TppT 12
1
Tp
=
1012,0
1
p
(3-14)
Như vậy ta có mạch vòng phản hồi tốc độ
Rf
Tp + 1
K
KnUf®
Tnp + 1 Jp
1
2Tp + 1
1
KQ
Q M
Hình 3-5. Mạch vòng phản hồi tốc độ.
Ta có: S02 =
)1)(12)(1()1)(12)(1(.
..
pTTppTp
K
pTTppTpJ
KKK
nn
nQ
(3-15)
Trong đó: K =
310.1050
62,20.520.2,0..
J
KKK nQ
= 0,02 (3-16)
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
66
Áp dụng tiêu chuẩn tối ưu đối xứng ta có bộ điều chỉnh tốc độ là bộ điều chỉnh
PID có công thức:
Rf = Kp + Kd.p +
p
K i
Hàm chuẩn tối ưu đối xứng có dạng:
FDX(p) =
3322 8841
41
pTpTpT
pT
(3-17)
Ta có:
Fh(p) = Rf.S02 = Rf.
)1)(12)(1( pTTppTp
K
n
(3-18)
Fk(p) =
h
h
F
F
1 )1)(12)(1(.
.
pTTppTpKR
KR
nf
f
(3-19)
Hàm truyền đạt kín của bộ điều tốc áp dụng theo tiêu chuẩn tối ưu đối
xứng:
Fk(p)= FĐX
)1)(12)(1(.
.
pTTppTpKR
KR
nf
f
=
3322 8841
41
pTpTp
pT
)41)(1)(12)(1()1(.8.. 2 TTTppTpTTKR nf
)41)(12)(1(.8.. 2 TTppTpTKR f
(Kp + Kd.p +
p
K i
)
)41)(12)(1(.8. 2 TTppTpTK
Thực hiện phép biến đổi ta được phương trình hệ số của phương trình đặc tính:
TTTKK
TKKTTT
TKKTT
i
p
d
6..8
..8.6.8
..8..8
2
22
22
(3-20)
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
67
5,437
10.10.2.8
1010.6
..8
6
875
)10.(10.2.8
10.10.6)10.(8
.8
.6.8
5,0
10.2
10
33
33
2
233
3323
2
2
3
3
TK
TT
K
TK
TTT
K
K
T
K
i
p
d
(3-21)
Vậy bộ điều chỉnh tốc độ tìm được có hàm truyền là:
Rf = 875+
p
5,437
+ 0,5.p (3-22)
3.3.3. Mô phỏng hệ thống điều chỉnh Turbine.
Mô phỏng hệ thống điều chỉnh turbine bằng phần mềm Mathlab Simulink của:
“ The math Works”.
3.3.3.1. Sơ đồ cấu trúc bộ điều tốc trong Simulink.
Mạch vòng điều khiển vị trí
Hàm truyền tương đương
Hình 3-6a. Sơ đồ mô phỏng bộ điều khiển vị trí.
Hình 3 -6b. Sơ đồ cấu trúc hệ thống điều tốc
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
68
3.3.3.2. Kết quả mô phỏng.
Hình 3 – 7. Kết quả mô phỏng với thông số PID như tính toán
Hình 3 – 8: Khi tăng hệ số Ki
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
69
Hình 3 – 9: Khi tăng hệ số Kp.
Hình 3-10: Đồ thị khởi động Turbine và cho ăn tải sau 40 giây.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
70
Kết luận chương
Từ kết quả mô phỏng đã phản ánh đúng quá trình làm việc của bộ điều tốc.
Quá trình khởi động bộ điều tốc turbine sau một thời gian hệ sẽ đạt tới trạng thái
xác lập, sau khi hệ ổn định tại tần số định mức f = 50Hz sẽ đóng tải vào khi đó tốc
độ sẽ bị giảm xuống tới một giá trị cho phép (< 2%) sau đó dưới tác động của bộ
điều chỉnh tần số sẽ về trạng thái định mức nhưng đặc tính tốc độ không đạt như
ban đầu. Từ kết quả mô phỏng cũng cho phép ta có thể hiệu chỉnh lại thông số của
bộ điều khiển PID nhằm đạt được đáp ứng như mong muốn, một phương pháp điều
khiển tốt hơn đó là phương pháp điều chỉnh tuyến tính LQ (Linear Quadratic).
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
71
CHƢƠNG IV: NÂNG CAO CHẤT LƢỢNG ĐIỀU CHỈNH TURBINE
4.1. Giới thiệu chung:
Lý thuyết điều khiển tối ưu tuyến tính kinh điển đã cải thiện được đặc tính động
học của các hệ thống, tuy nhiên ở một mô hình chính xác, các mô hình dưới dạng hàm
truyền trở nên phức tạp nếu số bậc đa thức quá lớn. Khi đơn giản hóa mô hình để phù
hợp với bài toán tổng hợp bộ điều chỉnh thì đặc tính động học của bộ điều chỉnh cũng
sẽ giảm. Gần đây, lý thuyết điều khiển tuyến tính hiện đại đã được ứng dụng rộng rãi
cho nhiều bài toán điều khiển, với sự trợ giúp mạnh mẽ của máy tính số, thời gian tính
toán đã được rút ngắn rất nhiều.
Phương pháp điều chỉnh tuyến tính LQ (Linear Quadratic) cho đối tượng có mô
hình trạng thái tuyến tính sẽ được tác giả ứng dụng để thiết kế bộ điều chỉnh tốc độ
Turbine thủy lực NMTĐ để nâng cao chất lượng điều chỉnh Turbine NMTĐ, cụ thể
cho các đối tượng là mô hình SS_LN_0 và SS_LN_1. Đây là các đối tượng phi tuyến
đã được tuyến tính hóa, xét ở lân cận điểm làm việc xác lập, có cấu trúc SISO và mang
đặc điểm của một hệ pha không cực tiểu.
Nội dung chương này sẽ trình bày tiếp cận phương pháp và tính toán, thiết kế các
khâu điều chỉnh và mạch vòng điều chỉnh LQ, bao gồm: điều chỉnh phản hồi âm trạng
thái LQR, điều chỉnh phản hồi đầu ra với tác động của nhiễu Gaussian LQG. Kỹ thuật
LTR áp dụng cho LQG cũng sẽ được tác giả đề cập để xét đến khả năng bộ điều chỉnh
LQG tiệm cận được chất lượng động học của LQR.
4.2. Cơ sở lý thuyết về điều chỉnh LQ.
4.2.1. Bộ điều chỉnh LQR.
Tính chất động học của bộ điều chỉnh LQR phải đảm bảo đưa hệ trở lại điểm cân
bằng (hoặc điểm làm việc) cũ nếu hệ bị tác động của nhiễu tức thời đánh bật hệ ra khỏi
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
72
điểm đó đến một điểm trạng thái khác, với yêu cầu chỉ tiêu chất lượng toàn phương
(Phiếm hàm mục tiêu):
t
t
TT
dttFuttExt ux
0
min)]()()()([
2
1
u(t))Q(x(t),
(4.1)
phải là nhỏ nhất khi hệ thống quay trở về trạng thái cũ.
E
R
nxn
xác định bán dương và F
R
rxr
xác định dương, tương ứng là ma trận
trọng số trạng thái E và điều khiển F, được quyết định bởi người thiết kế để mô tả mối
quan hệ giữa các trạng thái của đối tượng và chi phí của quá trình điều khiển.
Tuy nhiên cũng lưu ý rằng ta đang xét ở hệ lân cận điểm làm việc, tương ứng với
việc dịch chuyển gốc tọa độ của hệ đến điểm làm việc và xét ổn định với tín hiệu bé tác
động vào hệ này.
DuCx y
BuAx
dt
dx (4.2)
Ta đã có:
A
R
nxn
, B
R
nxm
, B11 R
nxm
, C
R
rxn
,D
Θ
Trong đó:
- Mô hình SS_LN_0: n = 5
- Mô hình SS_LN_1: n = 8
- Số biến đầu vào (biến điều khiển): m = 1.
- Số đầu ra: r = 1
Hệ số khuyếch đại của bộ điều khiển phản hồi âm trạng thái được xác định:
.K.BFR T-1LQR
(4.3)
Với K là nghiệm xác định bán dương của phương trình đại số Riccati (CARE).
E K.A -.K A -.K .BK.B.F TT-1
(4.4)
Bộ điều chỉnh LQR đảm bảo mạch vòng kín ổn định với chất lượng được quyết
định bởi các ma trận trọng số E, F đã chọn. Có một số phương pháp được nhiều tác giả
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
73
đề nghị cho việc lựa chọn giá trị khởi đầu của các ma trận E, F. Tuy nhiên, E và F vẫn
phải được hiệu chỉnh một cách phù hợp để có được đáp ứng điều chỉnh tốt nhất bằng
cách thử nghiệm đáp ứng.
Hình 4.1 mô tả sơ đồ nguyên lý điều chỉnh LQR.
Từ sơ đồ nguyên lý bộ điều chỉnh LQR, ta có hàm truyền mạch vòng hở giữa u và
u~
của sơ đồ điều chỉnh LQR được xác định:
.BA)-.(sIK(S)H -1LQRLQR
(4.5)
Phương pháp điều khiển LQR đạt được chất động học rất tốt, tuy nhiên nó đòi hỏi
tất cả các biến trạng thái của quá trình đều phải đo được, điều này vừa không thực tế,
vừa không kinh tế trong mọi trường hợp.
4.2.2. Bộ điều khiển LQG.
4.2.2.1. Bài toán tuyến tính có nhiễu.
Mục tiêu của bài toán là làm sao có thể ước lượng được các trạng thái của mô
hình MNTĐ ở (4.2), loại được tác động của nhiễu vào hệ, đồng thời phiếm hàm mục
tiêu (4.1) đạt giá trị nhỏ nhất. Mô hình trạng thái hệ tuyến tính với sự tác động của
nhiễu được mô tả:
y
x
nCx y
nBuAx
dt
dx
(4.6)
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
74
Trong đó nx(t), ny(t) lần lượt là nhiễu tác động vào quá trình (hệ thống) và nhiễu
đo ở đầu ra. Giả thuyết chúng là nhiễu ồn trắng, có kỳ vọng bằng không. ny(t) không
tương quan chéo với x(t) và nx(t) không tương quan chéo với x() với < t. Các nhiễu
nx(t), ny(t) có ma trận hàm tương quan là Nx R
nxn
xác định bán dương, Ny R
rxr
.
4.2.2.2. Bộ quan sát trạng thái (lọc) Kalman.
Bộ điều chỉnh LQR đòi hỏi tất cả các biến trạng thái đều phải được phản hồi, điều
này không thực tế và vì thế cần phải có các bộ quan sát. Bộ quan sát trạng thái Kalman
có nhiệm vụ ước lượng trạng thái của đối tượng tuyến tính (4.2). Phương trình trạng
thái của bộ quan sát Kalman:
xˆC =yˆ
)ˆ(ˆ
ˆ
yyLBuxA
dt
xd (4.7)
Trong đó
xˆ
là một ước lượng của trạng thái x và
yˆ
là ước lượng đầu ra y mà bộ
quan sát thực hiện được. (4.7) có thể viết lại:
y
u
LBxLCAxCyLBuxA
dt
xd
)(ˆ)()ˆ(ˆ
ˆ (4.8)
Hệ số L được tính:
1.. y
T NCPL
(4.9)
Với P là nghiệm xác định bán dương của CARE.
x
T-1
y
T N = A.P - P.A - .C.P.NP.C
(4.10)
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
75
Hình 4.2 cho thấy cấu trúc của bộ quan sát trạng thái Kalman.
Hình 4.2 Cấu trúc của bộ quan sát trạng thái Kalman
Ma trận Nx , Ny mô tả phổ nhiễu quá trình và nhiễu đo đầu ra, chúng chỉ là các
tham số của quá trình và được lựa chọn bởi người thiết kế. Nx phải là ma trận xác định
bán dương và Ny xác định dương.
4.2.2.3. Bộ điều chỉnh phản hồi đầu ra LQG.
Bộ điều chỉnh LQG là sự kết hợp giữa bộ điều chỉnh LQR và bộ quan sát trạng
thái Kalman để điều chỉnh và kháng nhiễu cho hệ tuyến tính (4.6), với đầu ra trạng thái
của bộ quan sát là đầu vào của bộ điều chỉnh LQR. Cấu trúc của bộ điều chỉnh LQR
được trình bày ở hình (4.3).
Hình 4.3: Sơ đồ nguyên lý bộ điều chỉnh LQG
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
76
Như vậy, về bản chất thì bộ điều chỉnh LQR và LQG đều có chung khâu điều
chỉnh RLQR. Sự khác biệt giữa chúng là ở bộ điều chỉnh LQG các biến trạng thái không
phải đo trực tiếp như LQR, mà chúng được ước lượng nhờ bộ quan sát Kalman. Chính
vì điều đó mà bộ điều chỉnh LQG có biên ổn định kém hơn bộ điều chỉnh LQR, nhưng
lại có khả năng kháng nhiễu tốt hơn LQR.
Hàm truyền mạch vòng hở giữa u và û của sơ đồ điều chỉnh LQR trên hình 4.3
được xác định:
.BA)-.LC.(sIL.C)B.KA-.(sIK(s)H -1-1LQRLQRLQG
(4.11)
4.2.2.4. Loop Transfer Recovery.
Phương pháp điều chỉnh phản hồi trạng thái LQR khó thực hiện được do nó yêu
cầu các trạng thái của hệ thống phải đo được. Việc sử dụng bộ quan sát Kalman để ước
lượng trạng thái của hệ thống từ tín hiệu ra là một giải pháp rất tốt để giải quyết vấn đề
trên, nhưng nó cũng làm cho đặc tính động học của bộ điều chỉnh LQR giảm chất
lượng.
Phương pháp Loop Transfer Recovery (LTR) dựa trên cấu hình của bộ điều chỉnh
LQG, nhằm mang lại cho LQG đặc tính điều chỉnh tiệm cận đến đặc tính điều chỉnh
của bộ điều chỉnh LQR, dù có sự tồn tại của nhiễu, trong đó các ma trận nhiễu Nx, Ny
có vai trò quyết định đến phương pháp LTR.
Viết lại hàm truyền mạch vòng hở của LQR và LQG từ (4.5) và (4.11):
.BA)-.(sIK(s)H -1LQRLQR
(4.12)
.BA)-(sI .LC.L.C)B.KA-.(sIK(s)H -1-1LQRLQRLQR
(4.13)
Với việc hiệu chỉnh hệ số P để Ny = p.No (No là ma trận phổ nhiễu đo) đạt giá trị
hợp lý, sao cho:
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
77
(s)H(s)H LQRLQG
0
Lim
(4.14)
Ta sẽ có được LTR, tức là bộ điều chỉnh LQG/LTR sẽ có đặc tính điều chỉnh của
LQR với cùng các ma trận trọng số E,F.
Để đạt được LTR, có thể kết hợp chọn
BBNx ,.
, σ càng lớn thì hệ càng dễ đạt
đến LTR, nhưng cũng có thể có những tác động không mong muốn do hệ số khuyếch
đại của hệ. Tuy nhiên với hệ pha không cực tiểu, khi các điểm không phía bên phải trục
ảo không quá xa trục ảo, thì kỹ thuật LTR ít có hiệu quả.
4.3. Phân tích tính điều khiển đƣợc và quan sát đƣợc.
4.3.1. Phân tích tính điều khiển được.
Để thuận tiện trong việc trình bày, các hệ số của ma trận A được tính, có trị số
khác không sẽ được thể hiện bằng ký hiệu chỉ số của nó trong ma trận, các hệ số có trị
số bằng không sẽ được ghi rõ bằng chính trị số “0” của nó.
Ma trận hệ thống của mô hình S_LN_1:
888786858281
78
666361
55
454441
3332
2321
17
SS_LN_1
aaaa00aa
a0000000
00a00a0a
000a0000
000aa00a
00000aa0
00000a0a
00a00000
A
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
78
Ma trận hệ thống của mô hình S_LN_0:
55
45444
3332
2321
151211
SS_LN_0
a0000
aa00a
00aa0
00a0 a
a00aa
A
Với quan hệ của các biến và trị số của các phần tử ma trận, ta dễ dàng kiểm tra và
khẳng định được, ma trận hệ số
SS_LN_1A
và ma trận
SS_LN_0A
không suy biến với mọi
giá trị của biến quá trình, tức là
5 A SS_LN_0 rank
và
8 A SS_LN_1rank
.
Tiêu chuẩn để xét tính điều khiển được cho hệ tuyến tính tham số hằng, thì điều
kiện cần và đủ để hệ tuyến tính (4.2) điều khiển được là:
nBAsIRank ,
C, trong đó n là hạng của ma trận A.
Ta có:
*
|, SS_LN_1BAsIRank
8888786858281
78
666361
555
454441
3332
2321
17
ba-saaa00aa
0as000000
000a-s00a0a
b000a-s0000
0000aa-s00 a
000000a-sa0
000000asa
00a00000 s
Rank
*
|, SS_LN_0BAsIRank
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
79
555
45444
3332
2321
151211
ba-s0000
0aa-s00a
000a-sa0
000as a
0a00aa-s
Rank
Dễ thấy rằng:
8 |, SS_LN_1 BAsIRank
5 |, SS_LN_0 BAsIRank
Như vậy, các hệ đang xét đều điều khiển được trong toàn miền phức.
4.3.2. Phân tích tính quan sát được.
Tiêu chuẩn để xét tính quan sát được cho hệ tuyến tính tham số hằng, thì điều kiện
cần và đủ để hệ tuyến tính (4.2) điều khiển được là:
Csn
C
BAsI
Rank
,
Ta có:
SS_LN_1|
,
C
BAsI
Rank
0000c000
a-saaa00aa
as000000
00a-s00a0a
000a-s0000
000aa-s00 a
00000a-sa0
00000asa
0a00000 s
4
888786858281
78
666361
55
454441
3332
2321
17
Rank
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
80
SS_LN_0|
,
C
BAsI
Rank
0c000
a-s0000
aa-s00a
00a-sa0
00asa
a00aa-s
4
55
45444
3332
2321
151211
Rank
Cũng dễ thấy rằng:
Cs
C
BAsI
Rank
8|
,
SS_LN_1
Cs
C
BAsI
Rank
5|
,
SS_LN_0
Như vậy các hệ đang xét đều quan sát được trong toàn miền phức.
4.4. Thiết kế bộ điều chỉnh LQ.
4.4.1. Thiết kế bộ điều chỉnh LQR.
Phần này sẽ xây dựng cấu trúc điều chỉnh LQR cho hai mô hình SS_LN_0 (các
biến trạng thái có thể đo được) và SS_LN_1 (một số biến trạng thái không đo được).
Với giả thiết các mô hình đang làm việc ổn định với tải
puPload 5,0
. Sau đó mô
phỏng với sự thay đổi tải của hệ để đánh giá chất lượng động học của bộ điều chỉnh.
4.4.1.1. Xây dựng cấu trúc bộ điều chỉnh LQR.
Nhiệm vụ của khâu điều chỉnh LQR là làm cho hệ ổn định với nhiễu tức thời tác
động vào hệ, nó mang lại những yếu tố quyết định đến chất lượng của hệ là tính ổn
định, bền vững và thời gian đáp ứng. Tuy nhiên, với tác động của tín hiệu điều chỉnh
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
81
với dạng hàm đơn vị w(t)=1(t), thì sai lệch tĩnh vẫn có thể tồn tại ngay cả khi hệ kín đã
ổn định. Nguyên nhân dẫn đến sai số tĩnh là do hàm truyền hở không có điểm cực,
ngoài ra trong hệ thống thực thì có thể do không có được bộ tham số chính xác tuyệt
đối khi tổng hợp bộ điều khiển.
Bài toán thiết kế LQR đã đưa ra mục tiêu là tìm khâu LQR để hệ kín ổn định với
ma trận phản hồi âm của hệ kín
LQRKBA .
có các giá trị riêng nằm bên trái trục ảo, nên
rõ ràng hệ kín với bộ điều chỉnh phản hồi âm LQR sẽ tồn tại sai lệch tĩnh. Vì vậy, để
khử sai lệch tĩnh ta cần phải bổ xung thêm thành phần tích phân vào tín hiệu sai lệch
e(t) giữa giá trị đặt và giá trị trạng thái (đầu ra) cần ổn định. Đây cũng là khâu quyết
định đến chất lượng của bộ điều chỉnh, đó là sai lệch tĩnh và quá điều chỉnh.
Các hình (4.4), (4.5) trình bày hai cấu trúc bộ điều chỉnh LQR có cùng bản chất,
chỉ khác nhau về việc lấy tín hiệu trạng thái hoặc đầu ra để so sánh với giá trị đặt.
Hình 4.4: Cấu trúc LQR với việc khử sai lệch tĩnh bằng tín hiệu trạng thái
Hình 4.5: Cấu trúc LQR với việc khử sai lệch tĩnh bằng tín hiệu đầu ra
Sơ đồ cấu trúc bộ điều chỉnh LQR với mạch vòng khử sai lệch tĩnh như hình (4.7).
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
82
Các hình (4.6), (4.7) trình bày hai cấu trúc mạch vòng điều chỉnh LQR, tương
ứng với hai cấu trúc bộ điều chỉnh đã xét ở hình (4.4) và (4.5).
Hình 4.6: Mạch vòng điều chỉnh LQR phản hồi trạng thái
Hình 4.7: Mạch vòng điều chỉnh LQR phản hồi trạng thái và đầu ra
4.4.1.2. Tính chọn tham số của bộ điều khiển.
Từ giả thiết hệ đang làm việc xác lập với mức tải
puPload 5,0
, đồng thời chọn
ma trận trọng số E, F như sau:
a. Mô hình SS_LN_1:
0,05F
4)j(i 1j)(i,E
4) j j,(i10j)(i,E
j) (i 0j)(i,E
8x8
3-8x8
8x8
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
83
Ta tính được hệ số của khâu LQR:
0.0081] 0.6566- 0.0002 4.8013 4.4612 0.0092 1.9395 [-3.6578KLQR
Các điểm cực của hệ kín:
0,0092]- 0,0274j;-0,0036- 0,0274j;0,0036- 1,2115;- 0,4757j;-2,4291-
0,4757j;2,4291- 13,2526j;-12,7312- 13,2526j;[-12,7312LQR_pol
b. Mô hình SS_LN_0:
0,03F
4)j(i 1j)(i,E
4)j j,(i 10j)(i,E
j) (i 0j)(i,E
5x5
3-5x5
5x5
Ta tính được hệ số của khâu LQR:
0.4477] 5.7595 0.0125 1.0935 [0.5287KLQR
Các điểm cực của hệ kín:
0.0273j]-0.0033- 0.0273j;0.0033- 1.1428;- 1.7986;- [-5.0158;LQR_pol
4.4.1.3. Kết quả mô phỏng.
Các hình (4.8) và (4.9) dưới đây trình bày kết quả mô phỏng đáp ứng điều chỉnh tốc
độ wr của Turbine tương ứng với các mô hình SS_LN_1 và SS_LN_0, khi có sự tăng
mức tải 0,1pu, giảm 0,05pu và tiếp tục giảm 0,02pu tương ứng ở các thời điểm t=50s,
150s và 250s.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
84
a: Đáp ứng tốc độ
r
b: Đáp ứng vị trí van
G
và công suất cơ
mecP
Hình 4.8: Đáp ứng điều chỉnh LQR của mô hình SS_LN_1 theo tải
a: Đáp ứng tốc độ
r
b: Đáp ứng vị trí van
G
và công suất cơ
mecP
Hình 4.9: Đáp ứng điều chỉnh LQR của mô hình SS_LN_0 theo tải
Nhận xét:
Trong cả hai mô hình SS_LN_0 và SS_LN_1, đáp ứng tốc độ có sự sụt tốc chỉ
khoảng 0,02pu (2%) khi tải tăng đột biến với mức tăng lên đến 0,1pu (10% định mức,
20% so với mức tải xác lập). Hệ ổn định bền vững, có sai lệch tĩnh bằng và thời gian
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
85
xác lập khoảng 10s với mạch vòng điều chỉnh cho đối tượng SS_LN_1 hay khoảng 15s
với đối tượng SS_LN_0
Tuy nhiên, với đáp ứng Pmec cho thấy sự bất hợp lý với Pload về điều kiện cân
bằng. Thực ra, đây là hệ quả của việc xét hệ tuyến tính hóa ở một điểm làm việc cách
xa điểm làm việc xác lập cũ (là cơ sở để tuyến tính hóa) sự vô lý này sẽ không tồn tại
khi ta xét ở lân cận gần với điểm làm việc xác lập.
4.4.2. Thiết lập bộ điều chỉnh LQG.
Phần này sẽ xây dựng cấu trúc điều chỉnh LQG cho hai mô hình SS_LN_0 và
SS_LN_1, với cấu trúc mô hình, các ma trận trọng số E, F, các giả thiết về điểm xác
lập ban đầu và các mức tải gia tăng trong quá trình mô phỏng như đã đề cập ở phần
thiết kế bộ điều chỉnh LQR (mục 4.4.1).
4.4.2.1. Xây dựng cấu trúc bộ điều chỉnh khiển LQG.
Nhiệm vụ của bộ quan sát Kalman trong điều chỉnh LQG làm nhằm mục đích ước
lượng các trạng thái của hệ thống, và không làm ảnh hưởng đến các điểm cực cũ của hệ
thống (điểm cực của bộ điều chỉnh LQR). Vì vậy, các kết tính toán các khâu LQR và
phân tích về sai lệch tĩnh đã được đề cập ở 4.4.1 cũng được áp dụng khi xây dựng cấu
trúc mạch vòng điều chỉnh LQG.
Ngoài ra, trong điều khiển LQG ta còn xét đến tác động của nhiễu hệ thống và
nhiễu đo vào hệ. Cấu trúc của bộ điều chỉnh LQG được điều chỉnh ở hình (4.10).
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
86
Hình 4.10: Cấu trúc bộ điều chỉnh LQG
Hình 4.11: Mạch vòng điều chỉnh LQG phản hồi đầu ra
4.4.2.2. Tính chọn tham số bộ điều khiển.
a. Mô hình SS_LN_1:
Chọn Nx =I
8x8
, Ny = 10
-1
.
Chọn nhiễu quá trình và nhiễu đo có biên độ xấp xỉ ±0,005pu (0,5%).
b. Mô hình SS_LN_0:
Chọn Nx =I
5x5
.10
-1
, Ny = 10
-1
.
Chọn nhiễu quá trình và nhiễu đo có biên độ xấp xỉ ± 0,005pu (0,5%).
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
87
4.4.2.3. Kết quả mô phỏng.
Kết quả mô phỏng đáp ứng điều chỉnh được trình bày ở hình (4.12) và (4.13).
a: Đáp ứng tốc độ
r
b: Đáp ứng vị trí van
G
và công suất cơ
mecP
Hình 4.12: Đáp ứng điều chỉnh LQG của mô hình SS_LN_1 theo tải và nhiễu
a: Đáp ứng tốc độ
r
b: Đáp ứng vị trí van
G
và công suất cơ
mecP
Hình 4.13: Đáp ứng điều chỉnh LQG của mô hình SS_LN_0 theo tải và nhiễu
4.4.2.4. Loop Transfer Recovery.
Trên cơ sở các bộ điều chỉnh LQR, LQG đã thiết kế, bây giờ ta xét đến đặc điểm
LTR đối với bộ điều chỉnh LQG. LTR thường được ứng dụng và có kết quả rất tốt đối
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
88
với hệ pha cực tiểu, với hệ pha không cực tiểu thì tác dụng của LTR không rõ ràng
hoặc chỉ có đáp ứng tốt nếu các điểm không phía bên phải trục ảo nằm xa trục ảo.
Áp dụng kỹ thuật LTR cho bộ điều chỉnh LQG đối với mô hình SS_LN_1, ta
chọn các ma trận hàm tương quan của nhiễu như sau:
Nx =σ.B.B
’
= B.B
’
(σ =1); Ny =ρ.N0 = 10
-5
(ρ=10-5; N0=1)
Gọi bộ điều chỉnh LQG tương ứng với các tham số này là LQG/LTR.
Gọi bộ điều chỉnh LQG đã xét ở 4.4.2.2.a là LQG(1).
Ta cũng sẽ xét thêm một bộ điều chỉnh LQG với các ma trận nhiễu khác nhằm có
cái nhìn đầy đủ hơn về tác động của các ma trận này đến chất lượng điều chỉnh. Chọn
bộ điều chỉnh LQG(2): Nx = 10
-1
.I
8x8
, Ny=10
-3
.
Xét đáp ứng điều chỉnh tốc độ của LQR, LQG(1), LQG(2), và LQG/LTR khi có
sự tăng/giảm mức tải 0,1pu ở điểm làm việc 0,5pu, nhiễu hệ thống, nhiễu đo có biên độ
dao động khoảng ±0,005pu.
Kết quả mô phỏng được trình bày ở hình (4.14). Hình (4.15) mô tả chi tiết các
đáp ứng ở hình (4.14).
Đáp ứng tốc độ cho thấy bộ điều chỉnh LQG/LTR đã có hình dáng đáp ứng gần
giống với LQR, tuy nhiên nhiễu tác động vẫn chưa được triệt tiêu. Hai bộ điều chỉnh
LQG(1) và LQG(2) có hình dáng đáp ứng giống nhau, nhưng với sự giảm nhỏ Nx và Ny
nên tác động của nhiễu lên LQG(2) đã giảm, đường đặc tính của nó mịn hơn rất nhiều
so với LQG(1) và LQG/LTR.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
89
Hình 4.14: Đáp ứng điều chỉnh tốc độ của LQR, LQG(1), LQG(2) và LQG/LTR
Hình 4.15: Chi tiết hình 4.14
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
90
* Nhận xét:
Với sự tác động của nhiễu hệ thống và nhiễu đo, trong cả hai mô hình SS_LN_0
và SS_LN_1, đáp ứng tốc độ của các bộ điều chỉnh LQG nói chung vẫn bị sụt tốc chỉ
khoảng 0,02pu (2%) khi tải tăng đột biến với mức tăng lên đến 0,1pu. Hệ ổn định bền
vững và thời gian xác lập khoảng 40s (LQG(1), LQG(2)), thời gian quá điều chỉnh 1%
chỉ tồn tại trong khoảng 10s. LQG(2) cho đáp ứng mịn hơn nhưng vẫn không cải thiện
được thời gian xác lập.
Đáp ứng dạng răng cưa vẫn tồn tại với biên độ khoảng ±0,0005pu (0,05%) ở các
bộ LQG(1), LQG/LTR, nhưng bộ điều chỉnh LQG/LTR có chất lượng động học rất tốt,
với thời gian đáp ứng và quá điều chỉnh đã đạt được chất lượng của bộ điều chỉnh
LQR.
Từ kết quả mô phỏng, có thể nói các bộ điều chỉnh LQG đã loại được đến 90%
tác động của nhiễu hệ thống và nhiễu đo vào hệ.
4.5. Điều khiển LQ cho mô hình phi tuyến.
Nhằm đánh giá kết quả xây dựng các mô hình trạng thái phi tuyến SS_LN và chất
lượng động học của các bộ điều khiển tuyến tính đã được xây dựng, mục này sẽ trình
bày các kết quả mô phỏng mạch vòng điều chỉnh LQR, LQG, xét ở lân cận điểm làm
việc xác lập cho các đối tượng phi tuyến SS_NL_1 và LL_NL_0
Quá trình mô phỏng ứng với tải xác lập 0,5pu. Các bộ điều chỉnh LQR, LQG
được xây dựng cho mô hình tuyến tính hóa tại điểm làm việc 0,5pu.
4.5.1. Điều khiển LQR cho mô hình phi tuyến.
Ở thời điểm t=25s, tăng tải cho hệ thêm 0,05pu và giảm 0,01pu ở thời điểm
t=175s (tương ứng +5% và -1% so với tải định mức, hay +10% và -2% so với tải xác
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
91
lập). Đáp ứng tốc độ
r
và
mecP
của hai mô hình SS_NL_1 và SS_NL_0 được so
sánh với kết quả của các mô hình tuyến tính hóa SS_LN_1 và SS_LN_0 (tương ứng
cùng chung bộ điều khiển), trình bày ở các hình (4.16) và (4.17).
a: Đáp ứng tốc độ
r
b: Đáp ứng công suất cơ
mecP
Hình 4.16: Đáp ứng điều chỉnh LQR của mô hình SS_LN_1 theo tải
a: Đáp ứng tốc độ
r
b: Đáp ứng công suất cơ
mecP
Hình 4.17: Đáp ứng điều chỉnh LQR của mô hình SS_LN_0 theo tải
Kết quả mô phỏng cho thấy:
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
92
- Mô hình SS_NL_1: Ứng với thay đổi tải 0,05pu, có mức chênh lệch sụt tốc so
với mô hình SS_LN_1 khoảng 0,0035pu, quá điều chỉnh khoảng 0,0015pu, thời gian
xác lập khoảng 18s. Với mức thay đổi tải 0,01pu, sự chênh lệch độ sụt tốc chỉ còn
khoảng 0,0008pu và quá điều chỉnh không đáng kể (0,00025pu).
- Mô hình SS_NL_0: Đáp ứng tốc độ và công suất cơ rất chính xác so với đáp
ứng của mô hình SS_LN_0. Chênh lệch về đáp ứng tốc độ gần như bằng 0, chênh lệch
về công suất cơ chỉ khoảng 0,0003pu ứng với mức thay đổi tải 0,05pu. Các chỉ tiêu
khác như thời gian xác lập, sai lệch tĩnh và sai lệch bám đều giống nhau.
4.5.2. Điều khiển LQG cho mô hình phi tuyến.
Nhiễu hệ thống và nhiễu đo tác động vào hệ được chọn có biên độ khoảng
±0,005pu (0,5%). Ở thời điểm t = 25s, giảm tải cho hệ 0,01pu và tăng 0,04pu ở thời
điểm t=100s. Đáp ứng tốc độ
r
và
mecP
của hai mô hình SS_NL_1 và SS_NL_0
được so sánh với kết quả của các mô hình tuyến tính hóa SS_LN_1 và SS_LN_0
(tương ứng cùng chung bộ điều khiển), trình bày ở các hình (4.18) và (4.19).
a: Đáp ứng tốc độ
r
b: Đáp ứng công suất cơ
mecP
Hình 4.18: Đáp ứng điều chỉnh LQG của mô hình SS_LN_1 theo tải và nhiễu
Kết quả mô phỏng cho thấy:
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
93
- Mô hình SS_NL_1: Ứng với thay đổi 0,01pu, độ sụt tốc chỉ khoảng 0,0035pu
và lệch so với mô hình SS_LN_1 khoảng 0,0005pu, thời gian xác lập khoảng 8s và
không có sự khác biệt lớn giữa đáp ứng của hai mô hình. Khi có sự thay đổi tải 0,04pu,
chênh lệch về độ sụt tốc khoảng 0,0022pu và quá điều chỉnh cũng không lớn (0,0055pu
so với 0,0028 với mô hình SS_LN_1). Thời gian xác lập không khác biệt, khoảng 20s.
a: Đáp ứng tốc độ
r
b: Đáp ứng công suất cơ
mecP
Hình 4.19: Đáp ứng điều chỉnh LQG của mô hình SS_LN_0 theo tải và nhiễu
- Mô hình SS_NL_0: Đáp ứng tốc độ và công suất cơ rất chính xác so với đáp
ứng của mô hình SS_LN_0. Chênh lệch về đáp ứng tốc độ gần như bằng 0, chênh lệch
về công suất cơ chỉ khoảng 0,0003pu ứng với mức thay đổi tải 0,04pu. Các chỉ tiêu
khác như thời gian xác lập, sai lệch tĩnh và sai lệch bám đều giống nhau.
4.6. Kết luận chƣơng.
Qua việc xây dựng cấu trúc các khâu điều chỉnh LQ, mạch vòng điều chỉnh, phân
tích kết quả mô phỏng đáp ứng điều chỉnh tốc độ cho hai mô hình SS_LN_0 và
SS_LN_1 đã cho thấy:
Các cấu trúc điều chỉnh LQ mà tác giả đã xây dựng thể hiện chất lượng động học
rất tốt, đáp ứng tốc độ với mạch điều chỉnh LQG có dao động rất nhỏ, khoảng
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
94
±0,005pu khi có nhiễu hệ thống và nhiễu đo (nhiễu Gaussian) với biên độ xấp xỉ
±0,005pu tác động vào hệ, hay nói cách khác, các bộ điều chỉnh LQG đã loại đến 90%
tác động của nhiễu hệ thống và nhiễu đo vào hệ khi xác lập. Chất lượng của các bộ
điều chỉnh phụ thuộc nhiều bởi các ma trận trọng số E, F ở khâu LQR và ma trận hàm
tương quan Nx, Ny của nhiễu khi thiết kế khâu quan sát Kalman. Chúng được lựa chọn
bằng cách thử đáp ứng và phân tích điểm cực của hệ kín trong quá trình thiết kế, xuất
phát các ma trận khởi điểm thường được chọn là ma trận đơn vị. Cũng thấy rằng chất
lượng của bộ điều chỉnh LQG không thể tốt hơn bộ điều chỉnh LQR, sự giảm nhỏ các
ma trận hàm tương quan của nhiễu có thể làm cho đáp ứng ít bị tác động của nhiễu,
nhưng vẫn không thể vượt qua được LQR về các chỉ tiêu như thời gian đáp ứng, quá
điều chỉnh và sai lệch tĩnh.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
95
KẾT LUẬN
Sau một thời gian nghiên cứu và làm việc nghiêm túc, được sự giúp đỡ nhiệt tình
của PGS -TS. Nguyễn Văn Liễn và các thầy giáo trong khoa, đến nay luận văn đã
hoàn thành đúng thời gian.
Tác giả đã giải quyết hoàn chỉnh các nội dung trong luận văn:
Đề cập đến những vấn đề tổng quan về cấu trúc, đặc điểm, chức năng làm việc
của các phần tử thuỷ lực, máy điện trong NMTĐ. Phân tích nguyên lý điều chỉnh tốc
độ Turbine và tác động của các yếu tố động học chất lỏng trên các phần tử thuỷ lực,
máy điện vào chúng.
Đi sâu vào nghiên cứu, phân tích quá trình động học phi tuyến của từng phần tử
thuỷ lực, máy điện, từ đó mô tả mô hình toán và xay dựng cấu trúc phi tuyến theo quan
hệ động học của chúng, trong đó hầu hết các hiệu ứng động lực học của nước xảy ra
trên các phần tử thuỷ lực đều được xem xét, chỉ bỏ qua hiệu ứng đàn hồi cột nước
trong đường ống dẫn do áp lực dòng chảy thấp.
Từ kết quả mô phỏng đã phản ánh đúng quá trình làm việc của bộ điều tốc. Quá
trình khởi động bộ điều tốc turbine sau một thời gian hệ sẽ đạt tới trạng thái xác lập,
sau khi hệ ổn định tại tần số định mức f = 50Hz sẽ đóng tải vào khi đó tốc độ sẽ bị
giảm xuống tới một giá trị cho phép (< 2%) sau đó dưới tác động của bộ điều chỉnh tần
số sẽ về trạng thái định mức. Với mô hình điều khiển tốc độ Turbine NMTĐ ở chương
3 dùng bộ điều khiển PID sau khi khởi động rất tốt nhưng sau khi cho ăn tải (sau 40
giây) tốc độ Turbine đã ổn định nhưng không đạt được tốc độ như ban đầu.
Xây dựng các cấu trúc mạch vòng điều chỉnh và thiết kế các khâu điều chỉnh LQ
tương ứng cho các mô hình SS_LN_0 và SS_LN_1, mô phỏng đáp ứng điều chỉnh để
so sánh giữa các bộ điều chỉnh LQR, LQG và LQG/LTR, đồng thời so sánh với đáp
ứng điều chỉnh theo phương pháp Lyapunov và các phương pháp dựa trên cấu trúc PID
để có đánh giá khách quan hơn về chất lượng động học của các bộ điều chỉnh đã xây
dựng.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
96
Về phân tích, đánh giá kết quả mô phỏng thực hiện trong quá trình xây dựng cấu
trúc phi tuyến, mô hình trạng thái và các mạch vòng điều chỉnh:
Sự đồng dạng đáp ứng giữa mô hình SS_LN và SS_NL cho thấy khi xét ở một
lân cận bé của điểm làm việc xác lập, các mô hình tuyến tính SS_LN thể hiện đầy đủ
và khá chính xác bản chất động học của mô hình phi tuyến SS_NL. Mô hình SS_NL_0
và SS_LN_0 sẽ rất hiệu quả để nghiên cứu ổn NMTĐ khi đường ống áp lực được chế
tạo bằng những vật liệu cứng (xem như cứng tuyệt đối và không có đàn hồi), khi đó có
thể sử dụng cấu trúc điều chỉnh LQR cho hệ, bởi 5 biến trạng thái của mô hình
SS_NL_0 đều đo được.
Các cấu trúc điều chỉnh LQ cho mô hình tuyến tính hoá cũng đã được tác giả xây
dựng thành công, đáp ứng tốc độ của tuabin khi có nhiễu tải tác động với độ lớn đến
0,1pu vẫn đảm bảo sụt tốc chỉ xấp xỉ 0,02pu với thời gian xác lập rất nhanh, bền vững
và sai lệch tĩnh rất bé. Các bộ điều chỉnh LQG và LQG/LTR ngoài việc đáp ứng nhanh
với sự thay đổi của tải, khả năng loại bỏ đến 90% tác động của nhiễu hệ thống và nhiễu
đo cũng là một kết quả tốt. Đặc biệt, kỹ thuật LTR đã mang lại cho bộ điều chỉnh LQG
chất lượng động học rất gần với bộ điều chỉnh LQR, tuy nhiên tác động của nhiễu lên
đầu ra vẫn tồn tại do đặc điểm của hệ pha không cực tiểu, khó tìm được LQG/LTR có
thể khử hoàn toàn tác động của nhiễu.
Cuối cùng, các bộ điều chỉnh được xây dựng từ việc tuyến tính hoá mô hình phi
tuyến được sử dụng để điều chỉnh cho đối tượng phi tuyến, sự sai lệch bé về độ sụt tốc
(lượng quá điều chỉnh <0,8% trong khi đó thông số này với bộ PID truyền thống là
5%) và sự đồng dạng về đáp ứng, cũng như các chỉ tiêu chất lượng khác như sai lệch
tĩnh, sai lệch bám và thời gian quá độ nhỏ, thời gian xác lập hầu như không thay đổi
đã khẳng định chất lượng động học của các bộ điều chỉnh đã được xây dựng.
Với kết quả đạt được như vậy, cho phép khẳng định rằng: luận văn đã tiếp cận và
xây dựng đúng phương pháp đối với mô hình toán, mô hình trạng thái của đối tượng.
Việc phân tích vai trò các khâu điều chỉnh, nguyên nhân tồn tại sai lệch tĩnh và đưa ra
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
97
giải pháp khử sai lệch tĩnh là hoàn toàn hợp lý, từ đó đã xây dựng đúng cấu trúc các
khâu điều chỉnh cũng như mạch vòng điều chỉnh, điều này đã được khẳng định một lần
nữa từ kết quả điều chỉnh mô hình trạng thái phi tuyến ở mục 4.5. Kết quả này hoàn
toàn có thể đáp ứng được yêu cầu của các hệ thống thực, khi NMTĐ là một hệ phi
tuyến có quá trình động học phức tạp và các hằng số thời gian lớn.
Vì điều kiện thời gian, nên luận văn mới dừng lại ở mức độ thành công về mặt
điều khiển theo phương pháp LQ cho tốc độ Turbine thuỷ điện tốt hơn khi dùng
phương pháp điều khiển truyền thống (bộ điều khiển PID), trong đó hệ được xét ở lân
cận điểm làm việc xác lập. Việc kết hợp kỹ thuật chuyển mạch hay Gain-Scheduling để
có thể điều khiển đối tượng trên toàn dải tải, nghiên cứu xây dựng chương trình điều
khiển để thử nghiệm trên thiết bị thuỷ lực là vấn đề mà tác giả dự định tiếp tục phát
triển nghiên cứu trong thời gian tới.
Xin chân thành cảm ơn sự giúp đỡ quý báu của Thầy hướng dẫn và các thầy trong
khoa đã tạo điều kiện giúp đỡ tôi hoàn thành luận văn. Kính mong các Thầy và các bạn
đồng nghiệp đóng góp ý kiến để kết quả nghiên cứu được tốt hơn.
Luận văn tốt nghiệp Cao học
Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên
98
Tµi liÖu tham kh¶o
[ 1] NguyÔn Do·n Ph•íc ( 2007), Lý thuyÕt ®Òu khiÓn tuyÕn tÝnh, NXB khoa häc
vµ kÜ thuËt Hµ Néi.
[ 2] NguyÔn Do·n Ph•íc ( 2007), Lý thuyÕt ®iÒu khiÓn phi tuyÕn, NXB khoa häc
vµ kÜ thuËt Hµ Néi.
[ 3] NguyÔn Do·n Ph•íc ( 2007), Lý thuyÕt ®Òu khiÓn n©ng cao, NXB khoa häc
vµ kÜ thuËt Hµ Néi.
[ 4] QuyÕt ®Þnh vÒ Quy ®Þnh ®Êu nèi vµo hÖ thèng ®iÖn quèc gia sè 37/2006/Q§-
BCN ngµy 16/10/2006 cña Bé C«ng nghiÖp.
[ 5] Bé m«n Tù ®éng ho¸, §¹i häc B¸ch khoa Hµ Néi (2006), Lý thuyÕt ®iÒu
khiÓn tù ®éng.
[ 6] Bé m«n Tù ®éng ho¸, §¹i häc B¸ch khoa Hµ Néi (2006), M« h×nh ho¸ vµ m«
pháng hÖ thèng ®iÒu khiÓn tù ®éng.
[ 7] §¹i häc B¸ch khoa §µ N½ng (2006), Turbine thuû lùc – C¸c thiÕt bÞ thuû lùc
vµ C«ng tr×nh tr¹m thuû ®iÖn.
[ 8] Bïi Quèc Kh¸nh, NguyÔn V¨n LiÔn, NguyÔn ThÞ HiÒn(2006), TruyÒn ®éng
®iÖn, NXB Khoa häc vµ Kü thuËt.
[ 9] Bïi Quèc Kh¸nh, NguyÔn V¨n LiÔn, Ph¹m Quèc H¶i, D•¬ng V¨n
Nghi(2006), §iÒu chØnh tù ®éng truyÒn ®éng ®iÖn, NXB Khoa häc vµ Kü thuËt.
TiÕng Anh:
[ 10] Control System Toolbox 8.3, LQR Control Design.
emos/shipping/control/dcdemo.html#12
[ 11] Gl;enn O.Brown, “The history of the Darcy-Weisbach Equation for Pipe
Flow Resistance”, Proceedings of the 150th Anniversary Conference of ASCE
Washington, D.C, November, 3-6, 2002, A. Fredrich, and J. Rogers eds.,
American Society of Civil Engineer, Reston, VA
[ 12] Harrry L. trentelman - Anton A. Stoorvogel - Malo hautus (2001), Control
theory for linear sytems, Springer London.
[ 13] H Goyal, M hanmandlu and D>P>Kothari (2006), “ A New Optimal Flow
Control Approach for Automatic Control of Small Hydro Power Plants”, JIE(I),
Vol 87.
[ 14] “ Hydraulic turbine and turbine control moldels for sytsem dynamic
studies”. IEEE Transon Power Systems Vol.7. No.1 pp. 167179, Februaary
1992.
Các file đính kèm theo tài liệu này:
- Đề tài- ĐIỀU KHIỂN TURBINE THUỶ ĐIỆN.pdf