Thuật toán đã trình bày ở trên, sử dụng công thức tích phân Euler cải biên bao gồm các bước tính toán sau:
Cao học 2012- 2014 67
(1) Tính CĐXL và tính E∠δ0 của nút trong của 3 máy phát và công suất của chúng P0(bằng công suất turbin không đổi PT)
(2) Lập ma trận Y sau khi xảy ra ngắn mạch
(3) Xác định gia số công suất của các máy phát ngay sau khi xảy ra nhiễu
(tính toán như ở bước (6) dưới đây)
Tại bước tích phân thứ k+1 (bắt đầu từ k=0):
(4) Dự đoán gia số tốc độ góc cho bước k+1 từ góc bước hiện tại (k). Với bước 1 (k=0) thì ∆ω0 = 0 và δ0 là góc CĐXL, D là hệ số cản. 1 [ ( ) . ] p k k k k j h P D T ω ω + ω δ ω ∆ = ∆ + ∆ − ∆
(5) Tính góc bước k+1 theo gia số tốc độ góc bước k và giá trị dự đoán ở bước k+1: ( ) 1 2 1 p h k k k k δ + =δ + ∆ω + ∆ω +
(6) Tính gia số công suất mới ∆P(δk+1):
- Giải V từ hệ phương trình: (δk+1):
0
VV VE
Y V+Y E=
Vì YVV hoàn toàn không thay đổi qua các bước lặp, vì vậy một cách giải tốt là tìm ma trận nghịch đảo vv việc này chỉ cần thực hiện một lần sau đó trong mỗi bước lặp chỉ cần tính điện áp: = * (- ve* ).
- Tính dòng máy phát: EE
EV
I=Y V+Y E
- Từ đó tính công suất của máy phát và gia số công suất ∆P = PT - P
(7) Hiệu chỉnh gia số tốc độ góc bước k+1: Tính theo ∆P cũ và ∆P mới Học viên: Phạm Đình Nguyện
Cao học 2012- 2014 68 ( ) 1 ( ) ( 1) 1 2 p k k k k k k j h P P D T ω ω + ω δ δ + ω ω + ∆ = ∆ + ∆ + ∆ − ∆ + ∆
Quay lại bước (4)
(8) Tích phân đến thời điểm cắt NM thì thay đổi ma trận Y theo lưới sau khi cắt NM rồi tiếp tục bước (3)-(6)
Tính toán được thực hiện trong Matlab. Tính cho 2 trường hợp không xét và có xét mômen cản.
Bước thời gian tích phân được chọn là h = 0.01 s Thời gian tính toán quá trình quá độ: 2 s
3.3 Thực hiện tính toán. 3.3. 1 Tính CĐXL của lưới .