Nghiên cứu phương pháp điều khiển cho thiết bị bù lõm điện áp nguồn trong công nghiệp

Nhan đề : Nghiên cứu phương pháp điều khiển cho thiết bị bù lõm điện áp nguồn trong công nghiệp Tác giả : Nguyễn Vũ Anh Tuấn Người hướng dẫn: Phạm Quang Đăng Từ khoá : Thiết bị bù lõm; Điện áp Năm xuất bản : 2020 Nhà xuất bản : Trường đại học Bách Khoa Hà Nội Tóm tắt : Tổng quan về khắc phục sự cố lồilõm điện áp ngắn hạn cho phụ tải; thuật toán điều chế vector không gian; thuật toán điều khiển; mô phỏng Matlab.

TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI LUẬN VĂN THẠC SĨ Nghiên cứu phương pháp điều khiển cho thiết bị bù lõm điện áp nguồn công nghiệp NGUYỄN VŨ ANH TUẤN tuannguyen694@gmail.com Ngành Kỹ thuật điều khiển tự động hóa Giảng viên hướng dẫn: TS Phạm Quang Đăng Viện: Điện HÀ NỘI, 11/2020 Chữ ký GVHD CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc BẢN XÁC NHẬN CHỈNH SỬA LUẬN VĂN THẠC SĨ Họ tên tác giả luận văn: Nguyễn Vũ Anh Tuấn Đề tài luận văn: Nghiên cứu phương pháp điều khiển cho thiết bị bù lõm điện áp nguồn công nghiệp Chuyên ngành: Kỹ thuật điều khiển tự động hóa Mã số SV: CB170308 Tác giả, Người hướng dẫn khoa học Hội đồng chấm luận văn xác nhận tác giả sửa chữa, bổ sung luận văn theo biên họp Hội đồng ngày 30/10/2020 với nội dung sau: Bổ sung kịch mô thấy khả bù lồi điện áp thiết bị AVC Sửa lại số hình vẽ bị mờ (hình 3.7, 4.5, 4.6) Bổ sung thêm trích dẫn tài liệu tham khảo Ngày 26 tháng 11 năm 2020 Giáo viên hướng dẫn Tác giả luận văn CHỦ TỊCH HỘI ĐỒNG BỘ GIÁO DỤC & ĐÀO TẠO TRƯỜNG ĐH BÁCH KHOA HÀ NỘI CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự - Hạnh phúc LUẬN VĂN THẠC SĨ Họ tên học viên: Nguyễn Vũ Anh Tuấn Khóa: 2017B Viện: Điện Tên đề tài: Ngành: Kỹ thuật điều khiển tự động hóa Nghiên cứu phương pháp điều khiển cho thiết bị bù lõm điện áp nguồn công nghiệp Nội dung đề tài: Nghiên cứu, tìm hiểu giải pháp khắc phục cố lồi, lõm điện áp ngắn hạn cho phụ tải, từ đề xuất phương pháp điều khiển phù hợp Giáo viên hướng dẫn: TS Phạm Quang Đăng Thời gian giao đề tài: 18/07/2018 Thời gian hoàn thành: 29/09/2020 Ngày 26 tháng 11 năm 2020 GIÁO VIÊN HƯỚNG DẪN LỜI CẢM ƠN Em xin bày tỏ lòng biết ơn sâu sắc tới TS Phạm Quang Đăng thầy cô Viện Kỹ thuật điều khiển tự động hóa, Trường Đại học Bách Khoa Hà Nội nhiệt tình đào tạo, hướng dẫn, tạo điều kiện thuận lợi cho em nghiên cứu khoa học, giúp em hồn thành luận văn cách tốt TĨM TẮT NỘI DUNG LUẬN VĂN Trong thực tiễn công nghiệp, xảy biến động điện áp ngắn hạn lưới, điều khiển biến tần tác động làm dừng động cấp nguồn biến tần để bảo vệ điều làm gián đoạn sản xuất Để đảm bảo sản xuất liên tục, phụ tải nhạy cảm cần phải trang bị thiết bị bảo vệ riêng để đảm bảo hoạt động liên tục có cố biến động điện áp ngắn hạn Các giải pháp chống biến động điện áp bao gồm sử dụng lọc tích cực, cấp nguồn liên tục, chống sụt áp bù điện áp tích cực AVC Trong giải pháp sử dụng AVC nghiên cứu với ưu điểm cho phép truyền lượng theo hai chiều, chất lượng làm việc tốt hơn, cho phép đáp ứng nhanh giữ cho điện áp tải gần khơng có biến động Để khắc phục trường hợp cố hệ thống điện ba pha biến thiên đối xứng không đối xứng, luận văn nghiên cứu giải pháp sử dụng hệ thống nghịch lưu ba pha bốn nhánh van Ở phương pháp điều chế vector ba chiều 3D-SVM sử dụng nhờ ưu điểm sử dụng liên kết DC cao, phối hợp hài hịa tổn thất chuyển mạch so với phương pháp PWM hình sin Hệ thống nghịch lưu ba pha bốn nhánh van mơ hình hóa tốn học khơng gian ba chiều, vector điện áp dao động không gian ba chiều Do đó, với việc sử dụng mơ hình hóa khơng gian ba chiều số lượng tính tốn lớn địi hỏi trợ giúp phần mềm tính tốn mơ Matlab Ngồi ra, cơng cụ Matlab Simulink Plecs sử dụng để mô hoạt động hệ thống nghịch lưu ba pha bốn nhánh van Kết mô cho thấy nghịch lưu đáp ứng tốt có cố điện áp bị dâng/ sụt ba pha không cân Việc bù biến động điện áp vấn đề cấp bách vận hành lưới điện kết luận văn vừa có ý nghĩa đóng góp cho lý thuyết điều khiển bù điện áp ngắn hạn cho phụ tải, vừa có ý nghĩa thực tiễn Luận văn bao gồm nội dung sau: Chương 1: Tổng quan khắc phục cố lồi/ lõm điện áp ngắn hạn cho phụ tải Chương 2: Thuật tốn điều chế vector khơng gian Chương 3: Thuật tốn điều khiển Chương 4: Mơ Matlab Chương 5: Kết luận Mặc dù em cố gắng thực tốt có thể, luận văn khơng thể tránh khỏi sai sót, mong nhận ý kiến đóng góp thầy Em xin trân thành cảm ơn! Hà Nội, ngày 26 tháng 11 năm 2020 HỌC VIÊN THỰC HIỆN MỤC LỤC CHƯƠNG TỔNG QUAN VỀ KHẮC PHỤC SỰ CỐ LỒI/ LÕM ĐIỆN ÁP NGẮN HẠN CHO PHỤ TẢI 1.1 Hiện tượng biến động điện áp ngắn hạn 1.2 Các giải pháp chống biến động điện áp ngắn hạn 1.3 1.2.1 Giải pháp dùng lọc tích cực (Active filter) 1.2.2 Giải pháp sử dụng cấp nguồn liên tục (UPS) 1.2.3 Giải pháp sử dụng chống sụt áp (Sag fighter) 1.2.4 Giải pháp sử dụng bù điện áp tích cực (AVC) Nguyên tắc hoạt động AVC 1.3.1 Nguyên tắc hoạt động 1.3.2 Cấu trúc chung bù điện áp tích cực CHƯƠNG THUẬT TỐN ĐIỀU CHẾ VECTOR KHƠNG GIAN 2.1 Máy biến áp nối tiếp biến đổi phía tải 2.2 Bộ biến đổi phía tải dùng cầu ba pha bốn nhánh 2.3 Điều chế vector, sơ đồ ba pha bốn nhánh van (3-D SVM) CHƯƠNG THUẬT TOÁN ĐIỀU KHIỂN 22 3.1 Phương pháp phân tích điện áp thành phần đối xứng 22 3.2 Mơ hình hóa nghịch lưu ba pha bốn nhánh 23 3.2.1 Mơ hình trung bình tín hiệu lớn 25 3.2.2 Mơ hình tín hiệu nhỏ 28 3.3 Phương pháp điều khiển 29 3.4 Cấu trúc điều khiển 30 3.5 Thiết kế điều khiển 32 CHƯƠNG MÔ PHỎNG MATLAB 34 4.1 Tính tốn mạch lực 34 4.2 Tính tốn điều khiển kênh d q 34 4.3 Tính tốn điều khiển kênh 36 4.4 Mô 37 4.5 4.4.1 Khối mạch lực 37 4.4.2 Khối mô 37 Kịch mô 38 4.6 Kết mô 39 4.6.1 Thời gian nạp tụ 40 4.6.2 Đáp ứng với cố sụt áp pha 41 4.6.3 Đáp ứng với cố sụt áp hai pha không cân 42 4.6.4 Đáp ứng với cố sụt áp ba pha không cân 43 4.6.5 Đáp ứng với cố dâng điện áp ba pha không cân 44 CHƯƠNG KẾT LUẬN 45 5.1 Kết luận 45 5.2 Hướng phát triển luận văn tương lai 45 TÀI LIỆU THAM KHẢO 46 PHỤ LỤC 47 PHỤ LỤC A1 Code mô Matlab a) Code tạo xung điều khiển van function XungRa = fcn(Valpha,Vbeta,Vgama,VA,VB,VC,Xung) % Dien_ap_DC Vg=520; % Ma tran tinh cac ty so TD11 = [1 1/2 -sqrt(3)/2 -1 sqrt(3) 0]; TD12 = [3/2 -sqrt(3)/2 -1/2 sqrt(3)/2 1/2 sqrt(3)/2 -1]; TD13 = [3/2 -sqrt(3)/2 0 sqrt(3) -1/2 -sqrt(3)/2 1]; TD14 = [3/2 -sqrt(3)/2 0 sqrt(3) -1 -1]; TD21 = [1 1/2 sqrt(3)/2 -1 -3/2 sqrt(3)/2 0]; TD22 = [3/2 sqrt(3)/2 -1/2 sqrt(3)/2 -1 -1]; TD23 = [3/2 sqrt(3)/2 -3/2 sqrt(3)/2 -1/2 -sqrt(3)/2 1]; TD24 = [3/2 sqrt(3)/2 -3/2 sqrt(3)/2 1/2 -sqrt(3)/2 -1]; TD31 = [-1/2 sqrt(3)/2 1/2 sqrt(3)/2 -1 -3/2 -sqrt(3)/2 0]; TD32 = [0 sqrt(3) -1/2 -sqrt(3)/2 -1 -1]; TD33 = [0 sqrt(3) 47 -3/2 -sqrt(3)/2 1]; TD34 = [0 sqrt(3) -3/2 -sqrt(3)/2 -1/2 sqrt(3)/2 1]; TD41 = [-1/2 sqrt(3)/2 -1 -1 -sqrt(3) 0]; TD42 = [-3/2 sqrt(3)/2 -1/2 -sqrt(3)/2 1/2 -sqrt(3)/2 -1]; TD43 = [-3/2 sqrt(3)/2 0 -sqrt(3) 1]; TD44 = [-3/2 sqrt(3)/2 0 -sqrt(3) 1/2 sqrt(3)/2 -1]; TD51 = [-1/2 -sqrt(3)/2 -1 -1 3/2 -sqrt(3)/2 0]; TD52 = [-3/2 -sqrt(3)/2 1 1/2 -sqrt(3)/2 -1]; TD53 = [-3/2 -sqrt(3)/2 3/2 -sqrt(3)/2 -1/2 sqrt(3)/2 1]; TD54 = [-3/2 -sqrt(3)/2 3/2 -sqrt(3)/2 1/2 sqrt(3)/2 -1]; TD61 = [-1/2 -sqrt(3)/2 1/2 -sqrt(3)/2 -1 3/2 sqrt(3)/2 0]; TD62 = [0 -sqrt(3) 1 1/2 sqrt(3)/2 -1]; TD63 = [0 -sqrt(3) 3/2 sqrt(3)/2 -1/2 sqrt(3)/2 1]; TD64 = [0 -sqrt(3) 3/2 sqrt(3)/2 48 -1 -1]; % -Xac dinh vi tri vector tham chieu -d1 = 0; d2 = 0; d3 = 0; t1 = 0; t2 = 0; t3 = 0; t4 = 0; t = 0; if ((Valpha >= && Vbeta >= 0) || (Valpha < && Vbeta < 0)) if (Vbeta > 0) if (abs(Valpha) >= 1/sqrt(3)*abs(Vbeta)) % -Lang Tru -if (VA >= && VB < && VC < 0) %Tu dien pnnn pnnp ppnp A11 = (1/Vg)*TD11*[Valpha;Vbeta;Vgama]; d1 = A11(1,1); d2 = A11(2,1); d3 = A11(3,1); t1 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t2 = d3 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d2 + d3 + (1-d1-d2-d3)/2; t = 1; end if (VA >= && VB >= && VC < 0) %Tu dien pnnn ppnn ppnp A12 = (1/Vg)*TD12*[Valpha;Vbeta;Vgama]; d1 = A12(1,1); d2 = A12(2,1); d3 = A12(3,1); t1 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t2 = d2 + d3 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d3 + (1-d1-d2-d3)/2; t = 1.1; end if (VA >= && VB >= && VC >= 0) %Tu dien pnnn ppnn pppn A13 = (1/Vg)*TD13*[Valpha;Vbeta;Vgama]; d1 = A13(1,1); d2 = A13(2,1); d3 = A13(3,1); 49 t1 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t2 = d2 + d3 + (1-d1-d2-d3)/2; t3 = d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 1.2; end if (VA < && VB < && VC < 0) %Tu dien pnnp ppnp nnnp A14 = (1/Vg)*TD14*[Valpha;Vbeta;Vgama]; d1 = A14(1,1); d2 = A14(2,1); d3 = A14(3,1); t1 = d1 + d2 + (1-d1-d2-d3)/2; t2 = d2 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 1.3; end else % -Lang Tru -if (VA >= && VB >= && VC < 0) %Tu dien ppnn ppnp npnn A21 = (1/Vg)*TD21*[Valpha;Vbeta;Vgama]; d1 = A21(1,1); d2 = A21(2,1); d3 = A21(3,1); t1 = d1 + d2 + (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d2 + (1-d1-d2-d3)/2; t = 2; end if (VA < && VB >= && VC < 0) %Tu dien ppnp npnn npnp A22 = (1/Vg)*TD22*[Valpha;Vbeta;Vgama]; d1 = A22(1,1); d2 = A22(2,1); d3 = A22(3,1); t1 = d1 + (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; 50 t3 = (1-d1-d2-d3)/2; t4 = d1 + d3 + (1-d1-d2-d3)/2; t = 2.1; end if (VA >= && VB >= && VC >= 0) %Tu dien ppnn npnn pppn A23 = (1/Vg)*TD23*[Valpha;Vbeta;Vgama]; d1 = A23(1,1); d2 = A23(2,1); d3 = A23(3,1); t1 = d1 + d3 + (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 2.2; end if (VA < && VB < && VC < 0) %Tu dien ppnp npnp nnnp A24 = (1/Vg)*TD24*[Valpha;Vbeta;Vgama]; d1 = A24(1,1); d2 = A24(2,1); d3 = A24(3,1); t1 = d1 + (1-d1-d2-d3)/2; t2 = d1 + d2 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 2.3; end end else if (abs(Valpha) >= 1/sqrt(3)*abs(Vbeta)) % -Lang Tru -if (VA < && VB >= && VC >= 0) %Tu dien nppn nppp nnpn A41 = (1/Vg)*TD41*[Valpha;Vbeta;Vgama]; d1 = A41(1,1); d2 = A41(2,1); d3 = A41(3,1); t1 = (1-d1-d2-d3)/2; t2 = d1 + d2 + (1-d1-d2-d3)/2; 51 t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = d2 + (1-d1-d2-d3)/2; t = 4; end if (VA < && VB < && VC >= 0) %Tu dien nppp nnpn nnpp A42 = (1/Vg)*TD42*[Valpha;Vbeta;Vgama]; d1 = A42(1,1); d2 = A42(2,1); d3 = A42(3,1); t1 = (1-d1-d2-d3)/2; t2 = d1 + (1-d1-d2-d3)/2; t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = d1 + d3 + (1-d1-d2-d3)/2; t = 4.1; end if (VA >= && VB >= && VC >= 0) %Tu dien nppn nnpn pppn A43 = (1/Vg)*TD43*[Valpha;Vbeta;Vgama]; d1 = A43(1,1); d2 = A43(2,1); d3 = A43(3,1); t1 = d3 + (1-d1-d2-d3)/2; t2 = d1 + d3 + (1-d1-d2-d3)/2; t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 4.2; end if (VA < && VB < && VC < 0) %Tu dien nppp nnpp nnnp A44 = (1/Vg)*TD44*[Valpha;Vbeta;Vgama]; d1 = A44(1,1); d2 = A44(2,1); d3 = A44(3,1); t1 = (1-d1-d2-d3)/2; t2 = d1 + (1-d1-d2-d3)/2; t3 = d1 + d2 + (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 4.3; end 52 else % -Lang Tru -if (VA < && VB < && VC >= 0) %Tu dien nnpn nnpp pnpp A51 = (1/Vg)*TD51*[Valpha;Vbeta;Vgama]; d1 = A51(1,1); d2 = A51(2,1); d3 = A51(3,1); t1 = d3 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = d2 + d3 + (1-d1-d2-d3)/2; t = 5; end if (VA >= && VB < && VC >= 0) %Tu dien nnpn pnpn pnpp A52 = (1/Vg)*TD52*[Valpha;Vbeta;Vgama]; d1 = A52(1,1); d2 = A52(2,1); d3 = A52(3,1); t1 = d2 + d3 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = d3 + (1-d1-d2-d3)/2; t = 5.1; end if (VA >= && VB >= && VC >= 0) %Tu dien nnpn pnpn pppn A53 = (1/Vg)*TD53*[Valpha;Vbeta;Vgama]; d1 = A53(1,1); d2 = A53(2,1); d3 = A53(3,1); t1 = d2 + d3 + (1-d1-d2-d3)/2; t2 = d3 + (1-d1-d2-d3)/2; t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 5.2; end if (VA < && VB < && VC < 0) %Tu dien nnpp pnpp nnnp 53 A54 = (1/Vg)*TD54*[Valpha;Vbeta;Vgama]; d1 = A54(1,1); d2 = A54(2,1); d3 = A54(3,1); t1 = d2 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + d2 + (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 5.3; end end end else if (Vbeta > 0) if (abs(Valpha) >= 1/sqrt(3)*abs(Vbeta)) % -Lang Tru -if (VA < && VB >= && VC < 0) %Tu dien npnn npnp nppp A31 = (1/Vg)*TD31*[Valpha;Vbeta;Vgama]; d1 = A31(1,1); d2 = A31(2,1); d3 = A31(3,1); t1 = (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = d3 + (1-d1-d2-d3)/2; t4 = d2 + d3 + (1-d1-d2-d3)/2; t = 3; end if (VA < && VB >= && VC >= 0) %Tu dien npnn nppn nppp A32 = (1/Vg)*TD32*[Valpha;Vbeta;Vgama]; d1 = A32(1,1); d2 = A32(2,1); d3 = A32(3,1); t1 = (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = d2 + d3 + (1-d1-d2-d3)/2; t4 = d3 + (1-d1-d2-d3)/2; t = 3.1; end 54 if (VA >= && VB >= && VC >= 0) %Tu dien npnn nppn pppn A33 = (1/Vg)*TD33*[Valpha;Vbeta;Vgama]; d1 = A33(1,1); d2 = A33(2,1); d3 = A33(3,1); t1 = d3 + (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = d2 + d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 3.2; end if (VA < && VB < && VC < 0) %Tu dien npnp nppp nnnp A34 = (1/Vg)*TD34*[Valpha;Vbeta;Vgama]; d1 = A34(1,1); d2 = A34(2,1); d3 = A34(3,1); t1 = (1-d1-d2-d3)/2; t2 = d2 + d3 + (1-d1-d2-d3)/2; t3 = d2 + (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 3.3; end else % -Lang Tru -if (VA >= && VB >= && VC < 0) %Tu dien ppnn ppnp npnn A21 = (1/Vg)*TD21*[Valpha;Vbeta;Vgama]; d1 = A21(1,1); d2 = A21(2,1); d3 = A21(3,1); t1 = d1 + d2 + (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d2 + (1-d1-d2-d3)/2; t = 2; end if (VA < && VB >= && VC < 0) %Tu dien ppnp npnn npnp 55 A22 = (1/Vg)*TD22*[Valpha;Vbeta;Vgama]; d1 = A22(1,1); d2 = A22(2,1); d3 = A22(3,1); t1 = d1 + (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d1 + d3 + (1-d1-d2-d3)/2; t = 2.1; end if (VA >= && VB >= && VC >= 0) %Tu dien ppnn npnn pppn A23 = (1/Vg)*TD23*[Valpha;Vbeta;Vgama]; d1 = A23(1,1); d2 = A23(2,1); d3 = A23(3,1); t1 = d1 + d3 + (1-d1-d2-d3)/2; t2 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t3 = d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 2.2; end if (VA < && VB < && VC < 0) %Tu dien ppnp npnp nnnp A24 = (1/Vg)*TD24*[Valpha;Vbeta;Vgama]; d1 = A24(1,1); d2 = A24(2,1); d3 = A24(3,1); t1 = d1 + (1-d1-d2-d3)/2; t2 = d1 + d2 + (1-d1-d2-d3)/2; t3 = (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 2.3; end end else if (abs(Valpha) >= 1/sqrt(3)*abs(Vbeta)) % -Lang Tru -if (VA >= && VB < && VC >= 0) %Tu dien pnpn pnpp pnnn 56 A61 = (1/Vg)*TD61*[Valpha;Vbeta;Vgama]; d1 = A61(1,1); d2 = A61(2,1); d3 = A61(3,1); t1 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + d2 + (1-d1-d2-d3)/2; t4 = d2 + (1-d1-d2-d3)/2; t = 6; end if (VA >= && VB < && VC < 0) %Tu dien pnpp pnnn pnnp A62 = (1/Vg)*TD62*[Valpha;Vbeta;Vgama]; d1 = A62(1,1); d2 = A62(2,1); d3 = A62(3,1); t1 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + (1-d1-d2-d3)/2; t4 = d1 + d3 + (1-d1-d2-d3)/2; t = 6.1; end if (VA >= && VB >= && VC >= 0) %Tu dien pnpn pnnn pppn A63 = (1/Vg)*TD63*[Valpha;Vbeta;Vgama]; d1 = A63(1,1); d2 = A63(2,1); d3 = A63(3,1); t1 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t2 = d3 + (1-d1-d2-d3)/2; t3 = d1 + d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 6.2; end if (VA < && VB < && VC < 0) %Tu dien pnpp pnnp nnnp A64 = (1/Vg)*TD64*[Valpha;Vbeta;Vgama]; d1 = A64(1,1); d2 = A64(2,1); d3 = A64(3,1); 57 t1 = d1 + d2 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 6.3; end else % -Lang Tru -if (VA < && VB < && VC >= 0) %Tu dien nnpn nnpp pnpp A51 = (1/Vg)*TD51*[Valpha;Vbeta;Vgama]; d1 = A51(1,1); d2 = A51(2,1); d3 = A51(3,1); t1 = d3 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = d2 + d3 + (1-d1-d2-d3)/2; t = 5; end if (VA >= && VB < && VC >= 0) %Tu dien nnpn pnpn pnpp A52 = (1/Vg)*TD52*[Valpha;Vbeta;Vgama]; d1 = A52(1,1); d2 = A52(2,1); d3 = A52(3,1); t1 = d2 + d3 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = d3 + (1-d1-d2-d3)/2; t = 5.1; end if (VA >= && VB >= && VC >= 0) %Tu dien nnpn pnpn pppn A53 = (1/Vg)*TD53*[Valpha;Vbeta;Vgama]; d1 = A53(1,1); d2 = A53(2,1); d3 = A53(3,1); t1 = d2 + d3 + (1-d1-d2-d3)/2; t2 = d3 + (1-d1-d2-d3)/2; 58 t3 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t4 = (1-d1-d2-d3)/2; t = 5.2; end if (VA < && VB < && VC < 0) %Tu dien nnpp pnpp nnnp A54 = (1/Vg)*TD54*[Valpha;Vbeta;Vgama]; d1 = A54(1,1); d2 = A54(2,1); d3 = A54(3,1); t1 = d2 + (1-d1-d2-d3)/2; t2 = (1-d1-d2-d3)/2; t3 = d1 + d2 + (1-d1-d2-d3)/2; t4 = d1 + d2 + d3 + (1-d1-d2-d3)/2; t = 5.3; end end end end % -So sanh tao xung thoi gian -u = 0; v = 0; w = 0; f = 0; if (Xung < (t1)/2 || Xung >= 1-(t1)/2) u = 1; else u = 0; end if (Xung < (t2)/2 || Xung >= 1-(t2)/2) v = 1; else v = 0; end if (Xung < (t3)/2 || Xung >= 1-(t3)/2) w = 1; else w = 0; end if (Xung < (t4)/2 || Xung >= 1-(t4)/2) f = 1; else f = 0; end XungRa = [u;v;w;f]; y2 = t; 59 b) Code tính tốn điều khiển %%% Thong_so_mach -%%% C = 1e-6; L = 110e-3; L1= 110e-3; R = 100; fc = 10000; s = tf('s'); Vg = 520; Vref=311; w = 2*pi*50; %%% Ham_truyen_doi_tuong_dq -%%% Gsdq = (Vg*4*(R+L1*s))/(4*C*L*L1*s^3+4*C*L*R*s^2+(L+4*L1)*s+4 *R); bode(Gsdq); [mag1,phase1]=bode(Gsdq,2*pi*fc); %%% Ham_truyen_PI_dq %%% f1=fc/20; PM = 55; theta = PM -(phase1 + 180); fz = fc*sqrt((1sin(theta*pi/180))/(1+sin(theta*pi/180))); fp = fc*sqrt((1+sin(theta*pi/180))/(1sin(theta*pi/180))); wz = 2*pi*fz; wp = 2*pi*fp; numc = [1/wz 1]; denc = [1/wp 1]; Gc1 = tf(numc,denc)*tf([1 2*pi*f1],[1 0]); [mag2,phase2]=bode(Gc1,2*pi*fc); kc=1/(mag1*mag2); Gc=kc*Gc1; %bo-dk-dq Gdq=Gc*Gsdq; bode(Gdq); %%% Ham_truyen_doi_tuong_0 -%%% 60 Gs0 = (Vg*4*(R+L1*s))/(4*4*C*L*L1*s^3+4*4*C*L*R*s^2+(4*L+4*L 1)*s+4*R); bode(Gs0); [mag0,phase0]=bode(Gs0,2*pi*fc); %%% Ham_truyen_PI_0 %%% f2=fc/20; PM_0 = 55; theta = PM_0 - (phase0 + 180); fz_0 = fc*sqrt((1sin(theta*pi/180))/(1+sin(theta*pi/180))); fp_0 = fc*sqrt((1+sin(theta*pi/180))/(1sin(theta*pi/180))); wz_0 = 2*pi*fz_0; wp_0 = 2*pi*fp_0; numc = [1/wz_0 1]; denc = [1/wp_0 1]; Gc0 = tf(numc,denc)*tf([1 2*pi*f2],[1 0]); [mag00,phase00]=bode(Gc0,2*pi*fc); kc0=1/(mag0*mag00); Gc2=kc0*Gc0; %bo-dk-0 G0=Gc2*Gs0; bode(G0); 61 ... Viện: Điện Tên đề tài: Ngành: Kỹ thuật điều khiển tự động hóa Nghiên cứu phương pháp điều khiển cho thiết bị bù lõm điện áp nguồn công nghiệp Nội dung đề tài: Nghiên cứu, tìm hiểu giải pháp khắc... văn: Nguyễn Vũ Anh Tuấn Đề tài luận văn: Nghiên cứu phương pháp điều khiển cho thiết bị bù lõm điện áp nguồn công nghiệp Chuyên ngành: Kỹ thuật điều khiển tự động hóa Mã số SV: CB170308 Tác giả,... sụt áp Cụ thể cho phép bù điện áp lưới 10% mức điện áp định mức (giữ 60% điện áp) điện áp pha sụt tới 40% điện áp sụt ba pha (giữ 70% điện áp) Hoạt động online AVC cho phép ? ?áp ứng nhanh giữ cho