ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH
TRƯỜ NG
NG ĐẠI HỌC BÁCH KHOA
KHOA ĐIỆN – ĐIỆN TỬ
BỘ MÔN THIẾT BỊ ĐIỆN
BK
TP.HCM
LUẬN VĂN TỐT NGHIỆP:
ĐIỀU KHIỂN MÁY PHÁT ĐIỆN CẢM Ứ NG
NG CẤP NGUỒN TỪ HAI
HAI PHÍA
GVHD: PGS.TS DƯƠNG HOÀI NGHĨA
SVTH: Nguyễn Lê Huy Bằng
MSSV: 40400160
Tp. Hồ Chí Minh, tháng 01/2010
Nhận Xét C ủa Giáo Viên H ướ
ướ ng
ng Dẫ n
Nhận Xét C ủa Giáo Viên Phản Biện
LỜ
I C Ả
M Ơ
N
Em xin trân tr ọng gử i lờ i cám ơ n chân thành và t ố
ố t đẹ p nhấ t đế nn::
Thầ y Dươ ng
ng Hoài Nghĩ a,
a, ngườ i đ ã t ận tâm, nhiệt tình hướ ng
ng d ẫ
ẫ n và cung cấ p
các tài liệu cần thiế t em trong suố t thờ i gian làm luận văn.
Các thầ y, cô trong khoa
kho a Điện – Điện t ử
ử đ ã t ận tình chỉ b
bảo và truyề n thụ nhữ ng
ng
kiế n thứ c vô giá cho em trong suố t thờ i gian học đại học.
Cuố i cùng, em r ấ
mẹ và anh chị đ ã chăm lo, cám ơ n bạn bè đ ã ủng hộ
ất biế t ơ n bố m
em hoàn thành luận văn này.
Tp. H ồ Chí Minh, tháng 01 năm 2010.
M Ụ
Ụ C L Ụ C
LỜ
I C
Ả M Ơ N
..........................
.............
...........................
...........................
...........................
............................
...........................
...........................
......................i
........i
M Ụ
C ...........................
........................................
...........................
............................
...........................
...........................
...........................
........................
........... ii
ỤC
LỤ C
M ục lục hình
...........................
..............
...........................
...........................
...........................
............................
...........................
...........................
....................iv
......iv
Ký hiệu dùng trong luận văn
...........................
..............
...........................
...........................
...........................
............................
......................vi
........vi
Tóm t ắ
ắt luận văn ............
..........................
...........................
...........................
............................
...........................
...........................
..........................
............ viii
Chươ ng
ng 1: HỆ THỐNG CHUYỂN ĐỔI NĂNG LƯỢNG GIÓ. .........
......................
.........................
............1
1.1
Năng luợ ng
ng gió............................
.........................................
............................
...........................
...........................
...........................
................
...1
1.1.1
Giớ i thiệu về phát điện bằng sức gió, những thuận lợ i và khó khăn. .........1
1.1.2
Tiềm năng về năng lượ ng
ng gió ở Vi
Việt Nam.............
..........................
............................
.......................
.........2
1.2
Hệ thống chuyển đổi năng lượ ng
ng gió........................................
....................................................
.........................
............2
1.2.1
Các thành phần của một hệ thống chuyển đổi năng lượ ng
ng gió. ..................2
1.2.2 Các loại hệ thống chuyển đổi năng lượ ng
ng gió. ............
.........................
...........................
...................
..... 3
1.2.2.1 Hệ thống turbine gió tốc độ cố định. ..............
...........................
...........................
............................
..............4
1.2.2.2 Hệ thống turbine gió tốc độ thay đổi biển đổi toàn bộ công suất. ...........5
1.2.2.3 Hệ thống turbine gió tốc độ thay đổi biến đổi một phần công suất.........5
1.3
Tốc độ gió và công suất của turbine. ..........................
........................................
...........................
.........................
............6
1.3.1
Mật độ gió. ...........................
.........................................
. mạ
mạch từ
từ, v.v…
Các thông số
số củ
của bộ
bộ điều khiể
khiển cũ
cũng như
như giá trị
trị đặ
đặtt đượ c giữ
giữ không đổ
đổii như
như trong
tr
trườ
ườ ng
ng hợ
hợ p mô phỏ
phỏng ở chế
chế độ
độ danh
danh đị
định
nh cho tấ
tất cả
cả các trườ
trườ ng
ng hợ
hợ p khả
khảo sát dướ
dướ i đây.
•
Kết quả 1: Khi tă
tăng điện trở
trở stator
stator và rotor lên 20%
ễ n Lê Huy B ằ nngg
ễ
Trang 67
Chươ ng
ng 5: S ơ
ơ đồ mô phỏng và k ế
ế t quả.
GVHD: PGS.TS Dươ ng
ng Hoài Nghĩ a
ng của dòng rotor
Đáp ứ ng
Các mặt tr ượ
ượ t S1, S2
Công suấ t tác d ụng và phản kháng đầu cự c stator
ễ n Lê Huy B ằ nngg
ễ
Trang 68
Chươ ng
ng 5: S ơ
ơ đồ mô phỏng và k ế
ế t quả.
GVHD: PGS.TS Dươ ng
ng Hoài Nghĩ a
Dòng đ iện stator và rotor pha a
Hình 5.15: Đ áp ứ
ng hệ thố
ng khi tă ng đ iệ n trở stator
stator và rotor lên 20%
•
Kết quả 2: Khi giả
giảm điện trở
trở stator
stator và rotor xuố
xuống 20%
ng của dòng rotor
Đáp ứ ng
Các mặt tr ượ
ượ t S1, S2
ễ n Lê Huy B ằ nngg
ễ
Trang 69
Chươ ng
ng 5: S ơ
ơ đồ mô phỏng và k ế
ế t quả.
GVHD: PGS.TS Dươ ng
ng Hoài Nghĩ a
Công suấ t tác d ụng và phản kháng đầu cự c stator
Dòng đ iện stator và rotor pha a
Hình 5.16: Đ áp ứ
ng hệ thố
ng khi giả m đ iệ n trở stator
stator và rotor xuố
ng 20%
•
Kết quả 3: Khi tă
tăng điện cả
cảm tả
tản stator và hổ
hổ c
cảảm giữ
giữa stator và rotor lên 20%
Đáp ứ ng
ng của dòng rotor
ễ n Lê Huy B ằ nngg
ễ
Trang 70
Chươ ng
ng 5: S ơ
ơ đồ mô phỏng và k ế
ế t quả.
GVHD: PGS.TS Dươ ng
ng Hoài Nghĩ a
Các mặt tr ượ
ượ t S1, S2
Công suấ t tác d ụng và phản kháng đầu cự c stator
Dòng đ iện stator và rotor pha a
Hình 5.17: Đ áp ứ
ng hệ thố
ng khi tă ng đ iệ n cả m và hỗ c
cả m
Nhận xét kết quả: Qua 3 trườ
trườ ng
ng hợ
hợ p tă
tăng, giả
giảm điện áp stator và rotor; tă
t ăng điện cả
cả m
tản cuộ
cuộđnáng
stator
và vớ
hớ
ổi cả
cmô
ảm phỏ
giữỏangcuộ
giữ
cuở ộ chế
n stator
và rotor.
KNhư
ết ưquả
qungả mô
ph
ng
th
ấyểnsaicósố
số
không
kể so
kể
vhổ
ph
ch
ế độ
độ danh
danh
đị
định.
nh.Kế
Nh
vậyphỏ
vậ
hệỏ s
hệ
số
ố đcho
iều thấ
khi
khiể
tính bề
bền vữ
vững cao vớ
vớ i sự
sự thay đổ
đổii các thông số
số điện trở
trở và
và điện cả
cảm stator và rotor.
ễ n Lê Huy B ằ nngg
ễ
Trang 71
Chươ ng
ng 6: K ế
ế t luận.
GVHD: PGS.TS Dươ ng
ng Hoài Nghĩ a
Chươ ng
ng 6: KẾT LUẬN.
6.1 Các vấn đề đượ c thự c hiện trong luận văn.
•
•
•
•
•
•
Tìm hiể
hiểu về
về l ĩ ĩ nh
nh vự
vực phát điện bằ
bằng sứ
sức gió, nhữ
những thuậ
thuận lợ
lợ i và khó khă
khăn trong
vi
việệc phát triể
triển nă
năng lượ
lượ ng
ng gió, tiề
tiềm nă
năng về
về n
năăng lượ
lượ ng
ng gió củ
của Việ
Việt Nam.
đổii nă
Tìm hiể
hiểu các hệ
hệ th
thốống chuyể
chuyển đổ
năng lượ
lượ ng
ng gió: sự
sự phân bố
bố gió, hiệ
hiệu suấ
suất củ
của
tuabin gió, vậ
vận hành hệ
hệ th
thốống nă
năng lượ
lượ ng
ng gió tố
tốc độ
độ c
cốố đị
định
nh và thay đổ
đổi.i.
Tìm hiể
hiểu các loạ
loại máy phát điện sử
sử dụ
dụng trong hệ
hệ th
thốống chuyể
chuyển đổ
đổii nă
năng lượ
lượ ng
ng
gió.
Tìm hiể
hiểu về
về mô hình toán họ
học củ
của máy phát điện cả
cảm ứng cấ
cấp nguồ
nguồn từ
từ hai
phía (DFIG), mô phỏ
phỏng mô hình DFIG trên MATLAB/SIMULINK.
Tìm hiể
hiểu về
về ph
phươ
ươ ng
ng pháp điều khiể
khiển trượ
trượ tt,, ứng dụ
dụng vào điều khiể
khiển độ
độcc llậập
công suấ
suất tác dụ
dụng và công suấ
suất phả
phản kháng phía stator củ
của máy phát DFIG.
Xây dự
dựng mô hình hệ
hệ th
thốống điều khiể
khiển trên MATLAB/SIMULINK, kế
kết quả
quả mô
phỏ
phỏng cho thấ
thấy có thể
thể điều khiể
khiển độ
độcc lậ
lập công suấ
suất tác dụ
dụng và công suấ
suất phả
phản
kháng phía stator thông qua việ
việc điều khiể
khiển hai thành phầ
phần dòng rotor. Từ
Từ đó ta
có thể
thể điều khiể
khiển độ
độcc lậ
lập công suấ
suất tác dụ
dụng để
để t tốối ưu công suấ
suất nhậ
nhận đượ c ttừừ
gió, đồ
đồng
ng thờ
thờ i có thể
thể điều khiể
khiển thu hay phát công suấ
su ất phả
phản kháng góp phầ
phần
điều chỉ
chỉnh hệ
hệ th
thốống điện, đây là ưu điểm lớ
lớ n nhấ
nhất củ
của DFIG so vớ
vớ i các máy
phát điện gió khác.
• Bộ điều khiể
khiển trượ
trượ t cho kế
kết quả
quả bám theo giá trị
trị đặ
đặtt tố
tốt,t,đđáp ứng nhanh, không
có vọ
vọt llốố, ảnh hưở
hưở ng
ng qua lạ
lại giữ
giữa hai thành phầ
phần dòng rotor không đáng kể
kể nên
có thể
thể điều khiể
khiển độ
độcc lậ
lập đượ c.
c. Bộ
Bộ điều khiể
khiển có tính bề
bền vữ
vững cao khi có sự
sự
thay đổ
đổii củ
của thông số
số máy phát điện trở
trở , điện cả
cảm stator và rotor.
6.2 Đề nghị và các hướ ng
ng phát triển của luận văn.
•
Thực hiệ
Thự
hiện hệ
hệ th
thốống điều khiể
khiển trên hệ
hệ th
thựực, có thể
thể dùng độ
động
ng cơ
cơ DC
DC tạ
tạo
momen đầu
đầu vào.
• Ti
Tiếếp tụ
tục xây dự
dựng mô hình mô phỏ
phỏng hệ
hệ th
thốống điều khiể
khiển từ
từ lúc khở
khở i độ
động,
ng, hòa
đồng
đồ
ng bộ
bộ, đến
đến khi xác lậ
lập. Chi tiế
tiết hóa mô hình mô phỏ
phỏng: mô phỏ
phỏng thêm các
bộ bi
biếến đổ
đổii công suấ
suất phía rotor và áp dụ
dụng hệ
hệ th
thốống điều khiể
khiển thông qua điều
khiểển các bộ
khi
bộ bi
biếến đổ
đổii công suấ
suất này.
• Sử dụ
dụng các hệ
hệ th
thốống điều khiể
khiển thông minh như
như m
mạạng neural, fuzzy logic…
vào phươ
phươ ng
ng pháp điều khiể
khiển trượ
trượ t nhằ
nhằm nâng cao hơ
hơ n nữ
nữa chấ
chất lượ
lượ ng
ng điều
khiểển.
khi
ễ n Lê Huy B ằ nngg
ễ
Trang 72
Tài Liệu Tham Khảo
[1] Phan Qu
Quốốc Dũ
D ũng, Tô Hữ
Hữu Phúc,“Truyề n Động Điện”. NXB ĐHQG Tp.HCM –
2008.
[2]
J.F. Manwell, J.G. McGowan , A.L. Rogers , “Wind energy explained. Theory,
Design and Application”
Applicatio n”. University of Massachusetts, Amherst, USA. JOHN WLEY
& SONS, LTD.
[3] Mukund R. Patel Ph.D., P.E, “Wind and Solar Power Systems”.
U.S.Merchant Marine Academy Kings Point, New York. CRC Press LLC-1999.
[4] Ronald Hirschorn, ”Sliding mode control design for nonlinear stytems”,
Chapter 2. June 3, 2003.
[5]
Hugh Piggott,“Windpower Workshop”.British Wind Energy Association.
[6]
S. K Salman and Babak Badrzadeh, ” New
New Approach for modelling Doubly-Fed
Induction Generator (DFIG) for grid-connection studies”.School of Engineering, The
Robert Gordon University, Schoolhill, Aberdeen, AB10 1FR, Scotland, U.K.
Truyề n Động Điện – Truyề n động động cơ không
không
[7] Nguyễ
Nguyễn Vă
Văn Nhờ
Nhờ , “ C ơ
ơ ssở Truy
đồng bộ”. Đạ
Đạii họ
học Bách Khoa Tp. Hồ
Hồ Chí Minh.
không đồng
[8] Dươ
Dươ ng
ng Hoài Ngh ĩ
Ngh ĩ a,
a, Trươ
Trươ ng
ng Vă
Văn Hiề
Hiền,”
n,” Đ
Điề u khiể n t ốố c độ động cơ không
bộ dùng mô hình nội”.Tạp chí khoa họ
học và công nghệ
nghệ,pp 42-46, số
số 38+39,2002.
[9]
A.Petersson, “ Analisys, modeling and control of doub
doubly
ly - fed induction
generators for wind turbines”. Ph.D thesis, Chalmers University of Technology,
Sweden, 2005.
[10]
Simulati on and
Balasubramaniam
Balasubramani am Babypriya ,Rajapalan Anita,” Modelling, Simulation
Analysis of douby fed
f ed induction generator ffor
or wind turbines” Journal of Electrical
Engineering, vol. 60, no. 2, 2009, 79–85.
[11]
Ahmed G. Abo-Khalil, Dong-Choon Lee, Seung-Pyo Ryu, “Synchronization of
DFIG Output Voltage to
t o Utility Grid in Wind Power System”. Dept. of Electrical. Eng.
Yeungnam Univ., 214-1, Daedong, Gyeongsan, Gyeongbuk, Korea.
[12] Mai Tuan Dat, Mai Ba Loc, Pham Dinh Truc, Nguyen Huu Phuc, “ Control of
Doubly-fed induction generator using DSpace R&D Controller Board - an application
of rapid control coordinated with matlab/simulink”. Division of Electrical Machines &
Equipment Engineering, HCMUT, Vietnam
[13]
Guofeng Yuan, Jianyun Chai, Yongdong Li, “Vector Control and
Synchronization of Doubly Fed Induction Wind Generator system”. Department of
Electrical Engineering, Tsinghua University, Beijing, P.R.China.
[14]
Ahmed G. Abo-Khalil,Dong-Cho
Abo-Khalil,Dong-Choon
on Lee, Se-Hyun Lee, “Grid Connection of
Doubly-Fed Induction Generators
Generator s in Wind Energy Conve
Conversion
rsion System”. Dept. of
Electrical. Eng. Yeungnam Univ., 214-1, Daedong, Gyeongsan, Gyeongbuk, Korea.
Ph ụ Lục
File cài đặ
đặtt tham số
số và vẽ
vẽ đồ
đồ th
thịị trong matlab/simulink
%Thong so may dien khong dong bo rotor day quan(DFIG) 20HP
Pn = 4000;
Vn = 400;
p= 2;
fn = 50;
Vrn = 950;
(V)
Irn = 2.7;
n = 1440;
Sb = Pn;
wb = 2*pi*fn;
we = wb;
wbm = wb/p;
Tb = Pn/wbm;
Zb = Vn*Vn/Sb;
Vm = Vn*sqrt(2/3);
%
%
%
%
Dien ap dinh muc (line to line) (V)
so cap cuc
tan so dinh muc (Hz)
Dien ap dinh muc rotor (line-to-line)
%
% toc do dinh muc (rpm)
% Momen co ban
% tro khang co ban
% bien do ap pha
Vb
= Vm;=
(3/2)*p;
Tfactor
Rs = 1.07;
Rr = 1.32;
Lls = 0.0066;
Llr = 0.0098;
Lm = 0.1601;
Ls = Lls+Lm;
Lr = Llr+Lm;
Rs1 = Rs;
Rr1 = Rr;
Lls1 = Lls;
Llr1 = Llr;
Lm1 = Lm;
Ls1 = Lls1+Lm1;
Lr1 = Llr1+Lm1;
LM = 1/(1/Lm + 1/Lls + 1/Llr);
J = 0.32;
H = J*wbm*wbm/(2
J*wbm*wbm/(2*Sb);
*Sb);
Domega = 0;
f = 0;
% thoi gian mo phong
tstop = 4.25;
% Momen vao
Tm_time = [0 1.25 1.25 tstop];
Tm_value = [0 0 -0.7 -0.7]*Tb;
% Thong so cua bo dieu khien truot
Ts = Ls/Rs;
Tr = Lr/Rr;
%
%
%
%
%
dien tro
dien tro
diem cam
dien cam
ho cam
stator
rotor
stator
rotor
% momen quan tinh
gamma = 1/Tr;
sigma = 1-(Lm*Lm)/(
1-(Lm*Lm)/(Ls*Lr);
Ls*Lr);
a1 = Lm1-(Ls1*Lr1)/Lm1;
a2 = Lr1/Lm1;
a3 = Ls1/Lm1;
a4 = Rs1/Lm1;
a5 = a3*Rr1/a1+a2*Rs1/a1;
k1 = 1800;
k2 = 800;
Ti = 0.01;
% Gia tri dat cho idr & iqr
idr_time = [0 0.99 0.99 1.25 1.25 tstop];
idr_value = [-15 -15 0 0 6 6];
iqr_time = [1 1.75 1.75 2.5 2.5 3.5 3.5 tstop];
iqr_value = [-6.5 -6.5 -4 -4 -9 -9 -6.5 -6.5];
thi');
);
disp('Go return de ve do thi'
disp('Go
keyboard
clf;
figure
plot(y(:,1),y(:,2),'-k'
plot(y(:,1),y(:,2),
'-k')
)
grid on
ylabel('i*dr
ylabel('i*dr
axis([1 4 -8
'Thoi
xlabel('Thoi
xlabel(
title('Thanh
title(
'Thanh
[A]'
[A]')
)
8]);
[s]')
)
gian [s]'
phan dong dat idr*'
idr*')
)
figure
'-k')
)
plot(y(:,1),y(:,3),'-k'
plot(y(:,1),y(:,3),
grid on
ylabel('idr
ylabel(
'idr [A]'
[A]')
)
axis([1 4 -8 8]);
'Thoi gian [s]')
xlabel('Thoi
xlabel(
[s]')
'Thanh phan dong idr'
idr')
)
title('Thanh
title(
figure
'-k')
)
plot(y(:,1),y(:,4),'-k'
plot(y(:,1),y(:,4),
grid on
[A]')
)
ylabel('iqr* [A]'
ylabel('iqr*
axis([1 4 -12 -2]);
xlabel('Thoi
xlabel(
'Thoi gian [s]'
[s]')
)
'Thanh phan dong dat iqr*'
iqr*')
)
title('Thanh
title(
figure
plot(y(:,1),y(:,5),'-k'
plot(y(:,1),y(:,5),
'-k')
)
grid on
ylabel('iqr
ylabel(
'iqr [A]'
[A]')
)
axis([1 4 -12 -2]);
'Thoi gian [s]')
xlabel('Thoi
xlabel(
[s]')
title('Thanh
title(
'Thanh phan dong iqr'
iqr')
)
figure
'-k')
)
plot(y(:,1),y(:,6),'-k'
plot(y(:,1),y(:,6),
grid on
ylabel('Te [N.m]'
ylabel('Te
[N.m]')
)
axis([1 4 -30 30]);
'Thoi gian [s]'
[s]')
)
xlabel('Thoi
xlabel(
[N.m]')
)
title('Moment dien tu, Te [N.m]'
title('Moment
figure
plot(y(:,1),y(:,7),'-k'
plot(y(:,1),y(:,7),
'-k')
)
grid on
'Ps [W]'
[W]')
)
ylabel('Ps
ylabel(
axis([1 4 -4e3 4e3]);
'Thoi gian [s]')
xlabel('Thoi
xlabel(
[s]')
title('Cong
title(
'Cong suat tac dung dau cuc stator, Ps [W]'
[W]')
)
figure
'-k')
)
plot(y(:,1),y(:,8),'-k'
plot(y(:,1),y(:,8),
grid on
ylabel('Qs
ylabel(
'Qs [VAr]'
[VAr]')
)
axis([1 4 -3e3 3e3]);
'Thoi gian [s]'
[s]')
)
xlabel('Thoi
xlabel(
'Cong suat phan khang dau cuc stator, Qs [VAr)]'
[VAr)]')
)
title('Cong
title(
figure
plot(y(:,1),y(:,9),'-k'
plot(y(:,1),y(:,9),
'-k')
)
grid on
ylabel('S1'
ylabel(
'S1')
)
axis([1 4 -8 8]);
'Thoi gian [s]'
[s]')
)
xlabel('Thoi
xlabel(
'Mat truot S1'
S1')
)
title('Mat
title(
figure
plot(y(:,1),y(:,10), '-k'
'-k')
)
grid on
ylabel('S2'
ylabel(
'S2')
)
axis([1 4 -12 12]);
'Thoi gian [s]'
[s]')
)
xlabel('Thoi
xlabel(
'Mat truot S2'
S2')
)
title('Mat
title(
figure
plot(y(:,1),y(:,11), '-k'
'-k')
)
grid on
'ias [A]'
[A]')
)
ylabel('ias
ylabel(
axis([1.2 1.4 -15 15]);
xlabel('Thoi
xlabel(
[s]')
'Thoi gian [s]')
'Dong dien stator pha - a, ias'
ias')
)
title('Dong
title(
figure
'-k')
)
plot(y(:,1),y(:,12), '-k'
grid on
'iar [A]'
[A]')
)
ylabel('iar
ylabel(
axis([1 4 -20 20]);
'Thoi gian [s]')
xlabel('Thoi
xlabel(
[s]')
title('Dong
title(
'Dong dien rotor pha - a, iar'
iar')
)
figure
'-k')
)
plot(y(:,1),y(:,13), '-k'
grid on
ylabel('vas
ylabel(
'vas [V]'
[V]')
)
axis([1 4 -400 400]);
'Thoi gian [s]')
xlabel('Thoi
xlabel(
[s]')
'Dien ap stator pha - a, vas'
vas')
)
title('Dien
title(
figure
'-k')
)
plot(y(:,1),y(:,14), '-k'
grid on
[V]')
)
ylabel('var [V]'
ylabel('var
axis([1.5 1.7 -200 200]);
'Thoi gian [s]')
xlabel('Thoi
xlabel(
[s]')
title('Dien
title(
'Dien ap rotor pha - a, var'
var')
)
figure
plot(y1(:,1),y1(:,2),'-k'
plot(y1(:,1),y1(:,2),
'-k')
)
grid on
ylabel('vdr
ylabel(
'vdr [V]'
[V]')
)
axis([1.5 1.7 -200 200]);
xlabel('Thoi
xlabel(
'Thoi gian [s]')
[s]')
vdr')
)
title('Thanh phan dien ap
title('Thanh
rotor vdr'
figure
plot(y1(:,1),y1(:,3),'-k'
plot(y1(:,1),y1(:,3),
'-k')
)
grid on
'vqr [V]'
[V]')
)
ylabel('vqr
ylabel(
axis([1.5 1.7 -200 200]);
xlabel('Thoi
xlabel(
[s]')
'Thoi gian [s]')
title('Thanh
title(
'Thanh phan dien ap rotor vqr'
vqr')
)