Studi aliran daya adalah salah satu topik yang dibahas dalam mata kuliah Analisis Sistem Tenaga Listrik (Power System Analysis) yang merupakan mata kuliah wajib jurusan teknik tenaga listrik di hampir semua universitas di seluruh dunia. Teori2 yang ada di dalamnya akan cukup membuat pusing, namun salah satu cara untuk membantu memahaminya adalah dengan memecahkan suatu persoalan.

Gambar di atas adalah one line diagram sistem tenaga listrik sederhana 3 bus, dengan generator di bus 1 dan 2. Tegangan di bus 1 adalah 1 + j0 per unit. Besar tegangan di bus 2 tetap, sebesar 1.03 pu dengan daya nyata generator sebesar 400 MW. Beban sebesar 500 MW dan 350 MVAR ada di bus 3. Admitansi saluran pada gambar dalam besaran per unit dan berbasis 100 MVA. Resistansi saluran dan suseptansi line charging diabaikan.
- Cari solusi aliran daya dengan menggunakan metode Newton-Raphson, tentukan besar phasor V2 dan V3. Lakukan 2 iterasi.
- Ulangi no 1 dengan algoritma Decoupled.
- Ulangi no 1 dengan algoritma Fast Decoupled.
Jawaban:
Klik untuk memperbesar gambar2 ini.
1. Untuk Newton-Raphson, langkah2nya sbb:
- Langkah pertama adalah merubah semua nilai ke dalam besaran per unit (pu).
- Kemudian menyusun Y matriks.
- Identifikasi nilai2 yang diketahui dan yang tidak diketahui.
- Perhatikan jenis bus, apakah slack bus/reference bus, load / PQ bus atau voltage controlled / PV bus.
- Jika tegangan di suatu bus tidak diketahui, asumsikan tegangannya 1 + j0 pu. Asumsi ini sering disebut sebagai flat start. Hal ini disebabkan, biasanya besar tegangan suatu bus tidak akan jauh dari 1 pu.
- (♣) Substitusi nilai2 yang diketahui untuk mendapatkan persamaan2 P dan Q yang diketahui.
- Hitung turunan parsial dalam Jacobian matrix yang akan kita susun.
- Hitung daya residu ΔP dan ΔQ.
- Setelah semuanya dihitung, kita siap mencari penyelesaian dari iterasi pertama, yaitu dengan menyelesaikan persamaan linier di atas. Disini terlihat persamaan tsb. mengandung matriks Jacobian 3 x 3, sehingga jika matriks tsb berpindah ruas, maka tentu saja kita harus mencari invers dari matriks tsb.
- Dengan mendapatkan solusi dari persamaan ini, maka didapatkan lah besar dan sudut fasa baru tegangan pada bus-bus untuk iterasi yang pertama.
- Langkah berikutnya adalah mengulang langkah2 di atas, mulai dari (♣) sampai sebanyak iterasi yang diinginkan.
2. Untuk metode Decoupled Load Flow, langkah2nya sama, kecuali J2 dan J3 dalam matriks Jacobian tidak dihitung, atau dianggap 0.
3. Untuk Fast Decoupled Load Flow,
- Tidak memerlukan inversi matriks Jacobian. Sebagai gantinya, dari matriks admitansi Y, dibentuk matriks B’ dan B”.
- Persamaannya sendiri,
- Langkah2 selanjutnya sama dengan metode NR/Decoupled.
Untuk mencari besar daya di slack bus, nilai2 yang diperoleh dari masing2 algoritma, disubstitusi kembali ke persamaan daya. Sedangkan rugi2 sistem (losses), dapat dengan mudah dicari dengan prinsip hukum kekekalan energi, yaitu daya yang masuk ke suatu saluran sama dengan daya yang keluar ditambah rugi2.
Jika anda jeli, maka hasil perhitungan dengan kalkulator tangan di atas, ada yang bisa dikritisi, seperti losses yang negatif. Hal ini akibat adanya pembulatan nilai pada setiap langkah perhitungan, misal hasil yang sebenarnya 0.9333…dst tapi dibulatkan jadi 0.933. Dengan software2 komersial, komputer dengan mudah akan memberikan solusi yang lebih akurat.
Jika anda sudah cukup terlatih dengan pemecahan soal2 semacam ini, maka anda dengan mudah akan memahami cara kerja software2 komersial yang biasa digunakan dalam suatu sistem tenaga listrik.















37 tanggapan so far ↓
riski // 25 Agustus 2008 pada 11:45 |
wew walaupun saya anak komputer,
kalo ngeliat ada anak power kayak gini jadi salut saya….
salam kenal mas, dari anak komputer elektro unila.
bagus // 25 Agustus 2008 pada 12:17 |
wow
Analisa sistem Tenaga…( 6_6)
keren
ini topik yang aku suka,
dua semester yang lalu aku ikut kelas analisa sistem tenaga
sepertinya dari buku Pak Saadat ,
Keren sekali ini
Muhammad Imaduddin // 22 September 2008 pada 05:13 |
Thanks to riski dan bagus buat komennya.
Benar yang dibilang bagus, saya mengambil referensi utama-nya dari Power System Analysis-nya Hadi Saadat.
Buku ini cukup bagus. Untuk perhitungan analisis gangguan tidak seimbang, buku Power System Analysis dari Grainger-Stevenson juga sangat bagus.
Heri Budi Utomo // 18 November 2008 pada 04:34 |
mengapa alariran daya persamaanya S=P-JQ bukan S=P+JQ
terima kasih
Heri Budi Utomo // 18 November 2008 pada 04:36 |
bagaimana mendapatkan sofware powerworld?
Muhammad Imaduddin // 20 November 2008 pada 05:05 |
Pak Heri, untuk mendapatkan software PowerWorld, bapak dapat mendownloadnya gratis di situs yang telah saya sebutkan diatas. Di situs itu bapak akan diminta mengisi form yang minta nama, alamat email. Setelah bapak isi, bapak akan mendapatkan email balasan yang berisi link untuk mendowload software ini.
Mengenai mengapa S =(sigma) P – (sigma) jQ akan saya tulis dalam satu artikel tersendiri, namun kira2 jawabannya adalah bahwa persamaan ini merupakan turunan dari persamaan S=P+jQ.
tanjung // 23 Desember 2008 pada 15:23 |
wow… hebat bner mas analisisnya…. walaupun saya hanya d1 PLN untk calon oprtor GI, ingin bisa untk menganalisis STL kyak gini…
adecib3r // 27 Januari 2009 pada 06:32 |
Salam buat anak power (tegangan tinggi) saya juga mahasiswa Teknik Elektro Spesialisasi Sistem Tenaga.. saya juga sudah membuat analisa secara komputasi menggunakan python metode gauss-seidel silakan di lihat di http://admin.inherent-unib.net
Salam kenal,
jhe // 8 Februari 2009 pada 15:36 |
bagus.. metodenya .. thank infonya besok saya pendadaran…lumayan ngiget materi pas semester 5, STL ..
lam kenal, mampir ya di :
fuzairin.blogspot.com
jhe // 8 Februari 2009 pada 15:37 |
thank, dari saya
yuni andri // 10 Februari 2009 pada 16:21 |
wah…hebat ya bisa load flow..
saya masih belajar ni…maksih buat artikelnya ya mas….
yuni andri // 10 Februari 2009 pada 16:28 |
wah…hebat bisa load flow…..
makasih buat artikelnya mas….
membantu banget….
agus rohadi // 19 Februari 2009 pada 11:47 |
saya sangat tertarik pembahasan ini. tks
Muhammad Imaduddin // 20 Februari 2009 pada 08:36 |
Terima kasih atas tanggapan teman2 semua, saya senang saya bisa menulis sesuatu yang bermanfaat untuk kita semua.
Kemarin kunjungan ke situs ini mencatat rekor baru yang sebelumnya belum pernah terjadi.
Kemarin
Cari Tampilkan
analisis sistem tenaga listrik 46 More stats
metode aliran daya newton raphson 5 More stats
Kemarin
Judul Tampilkan
Contoh Penyelesaian Aliran Daya Listrik 50 More stats
PowerWorld: Cara Mudah untuk Menganalisi 7
Statistik Blog Summary Tables
Total Ditampilkan: 6.376
Hari Tersibuk: 93 — Thursday, February 19, 2009
Untuk pertama kalinya dalam sehari situs ini ditengok 93 kali. Tentu saja, hal ini positif, mengingat blog ini sebenarnya terlalu spesifik dan teknis (menurut opini saya). Di tengah begitu beragamnya situs2 lain yang jauh lebih menarik, ternyata minat belajar teman2 web surfer dari Indonesia(saya berasumsi yang mengunjungi blog saya orang Indonesia, dilihat dari statistik dari wordpress) sangat lah bagus.
Sekali lagi, saya mengucapkan terima kasih
wen // 26 Februari 2009 pada 08:18 |
Salam kenal mas.
Nak Power dari TE USK.
Jadi ingat AST 1 dan 2 ni.
Salut buat analisanya.
Taufik // 20 Maret 2009 pada 01:03 |
Assalamu alaikum wr.wb.
terima kasih banyak atas contoh yang diberikan. Sangat bermanfaat buat kami. Tapi kami masih bingung dengan soal asli dari model ine line diagram karena nilainya sudah Y (admitansi) bukan Z (impedansi). Kami ingin uji coba dengan software lain, misalnya Power world. Bagaimana data lengkapnya ?
terima kasih.
Hormat kami,
Taufik
Muhammad Imaduddin // 23 Maret 2009 pada 23:06 |
Pak Taufik,
Di PowerWorld, bapak cukup ganti nilai admitansi dengan impedansi, z = 1/y.
Misal y = -j50, maka z = 1 / (-j50) = j 0.02, jadi nilai 0.02 ini dimasukkan di Series Reactance (X) pada Transmission Line Option.
Namun harus diingat, jika admitans mengandung nilai real dan imaginer, maka impedans-nya harus lah invers bilangan kompleks tsb.
Jangan lupa merubah SetPoint Voltage = 1.03 pu pada Generator yang mengontrol tegangan di bus no 2.
Wassalam,
Selamat mencoba.
ita // 24 Maret 2009 pada 01:28 |
Saya juga sangat terbantu dengan karya pak Imaduddin.
Mohon bantuannya data lengkap yang akan diinput ke ETAP. saya ingin bandingkan hasilnya antara manual, matlab, power world, dan ETAP.
terima kasih
Hormat saya,
ITA
dante // 15 April 2009 pada 08:14 |
perkenalkan,,
nama dante, sedang bikin TA tentang loadflow menggunakan matlab dengan standart IEEE 18 bus.
ada yang tau????????????????
pake metode Gauss-seidel….
makasih…..:)
Muhammad Imaduddin // 15 April 2009 pada 13:40 |
Untuk Ita:
saya belum pernah coba di ETAP, saya lihat di situsnya, untuk sebuah software trial ETAP kita harus download 400MB
Untuk dante:
mungkin masud Dante IEEE 118 bus ya. Saya copy kan source code dari salah satu toolbox-nya Matlab, Matpower, data IEEE 118 bus, untuk metode Gauss Seidel. Menurut saya sudah banyak yang nulis di buku2 untuk metode ini, dante tinggal ganti data busnya sesuai yang diinginkan
function [baseMVA, bus, gen, branch, areas, gencost] = case118
%CASE118 Power flow data for IEEE 118 bus test case.
% Please see ‘help caseformat’ for details on the case file format.
% This data was converted from IEEE Common Data Format
% (ieee118cdf.txt) on 20-Sep-2004 by cdf2matp, rev. 1.11
% See end of file for warnings generated during conversion.
%
% Converted from IEEE CDF file from:
% http://www.ee.washington.edu/research/pstca/
% With baseKV data take from the PSAP format file from the same site,
% added manually on 10-Mar-2006.
%
% 08/25/93 UW ARCHIVE 100.0 1961 W IEEE 118 Bus Test Case
% MATPOWER
% $Id: case118.m,v 1.7 2007/09/17 16:07:48 ray Exp $
%%—– Power Flow Data —–%%
%% system MVA base
baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
bus = [
1 2 51 27 0 0 1 0.955 10.67 138 1 1.06 0.94;
2 1 20 9 0 0 1 0.971 11.22 138 1 1.06 0.94;
3 1 39 10 0 0 1 0.968 11.56 138 1 1.06 0.94;
4 2 39 12 0 0 1 0.998 15.28 138 1 1.06 0.94;
5 1 0 0 0 -40 1 1.002 15.73 138 1 1.06 0.94;
6 2 52 22 0 0 1 0.99 13 138 1 1.06 0.94;
7 1 19 2 0 0 1 0.989 12.56 138 1 1.06 0.94;
8 2 28 0 0 0 1 1.015 20.77 345 1 1.06 0.94;
9 1 0 0 0 0 1 1.043 28.02 345 1 1.06 0.94;
10 2 0 0 0 0 1 1.05 35.61 345 1 1.06 0.94;
11 1 70 23 0 0 1 0.985 12.72 138 1 1.06 0.94;
12 2 47 10 0 0 1 0.99 12.2 138 1 1.06 0.94;
13 1 34 16 0 0 1 0.968 11.35 138 1 1.06 0.94;
14 1 14 1 0 0 1 0.984 11.5 138 1 1.06 0.94;
15 2 90 30 0 0 1 0.97 11.23 138 1 1.06 0.94;
16 1 25 10 0 0 1 0.984 11.91 138 1 1.06 0.94;
17 1 11 3 0 0 1 0.995 13.74 138 1 1.06 0.94;
18 2 60 34 0 0 1 0.973 11.53 138 1 1.06 0.94;
19 2 45 25 0 0 1 0.963 11.05 138 1 1.06 0.94;
20 1 18 3 0 0 1 0.958 11.93 138 1 1.06 0.94;
21 1 14 8 0 0 1 0.959 13.52 138 1 1.06 0.94;
22 1 10 5 0 0 1 0.97 16.08 138 1 1.06 0.94;
23 1 7 3 0 0 1 1 21 138 1 1.06 0.94;
24 2 13 0 0 0 1 0.992 20.89 138 1 1.06 0.94;
25 2 0 0 0 0 1 1.05 27.93 138 1 1.06 0.94;
26 2 0 0 0 0 1 1.015 29.71 345 1 1.06 0.94;
27 2 71 13 0 0 1 0.968 15.35 138 1 1.06 0.94;
28 1 17 7 0 0 1 0.962 13.62 138 1 1.06 0.94;
29 1 24 4 0 0 1 0.963 12.63 138 1 1.06 0.94;
30 1 0 0 0 0 1 0.968 18.79 345 1 1.06 0.94;
31 2 43 27 0 0 1 0.967 12.75 138 1 1.06 0.94;
32 2 59 23 0 0 1 0.964 14.8 138 1 1.06 0.94;
33 1 23 9 0 0 1 0.972 10.63 138 1 1.06 0.94;
34 2 59 26 0 14 1 0.986 11.3 138 1 1.06 0.94;
35 1 33 9 0 0 1 0.981 10.87 138 1 1.06 0.94;
36 2 31 17 0 0 1 0.98 10.87 138 1 1.06 0.94;
37 1 0 0 0 -25 1 0.992 11.77 138 1 1.06 0.94;
38 1 0 0 0 0 1 0.962 16.91 345 1 1.06 0.94;
39 1 27 11 0 0 1 0.97 8.41 138 1 1.06 0.94;
40 2 66 23 0 0 1 0.97 7.35 138 1 1.06 0.94;
41 1 37 10 0 0 1 0.967 6.92 138 1 1.06 0.94;
42 2 96 23 0 0 1 0.985 8.53 138 1 1.06 0.94;
43 1 18 7 0 0 1 0.978 11.28 138 1 1.06 0.94;
44 1 16 8 0 10 1 0.985 13.82 138 1 1.06 0.94;
45 1 53 22 0 10 1 0.987 15.67 138 1 1.06 0.94;
46 2 28 10 0 10 1 1.005 18.49 138 1 1.06 0.94;
47 1 34 0 0 0 1 1.017 20.73 138 1 1.06 0.94;
48 1 20 11 0 15 1 1.021 19.93 138 1 1.06 0.94;
49 2 87 30 0 0 1 1.025 20.94 138 1 1.06 0.94;
50 1 17 4 0 0 1 1.001 18.9 138 1 1.06 0.94;
51 1 17 8 0 0 1 0.967 16.28 138 1 1.06 0.94;
52 1 18 5 0 0 1 0.957 15.32 138 1 1.06 0.94;
53 1 23 11 0 0 1 0.946 14.35 138 1 1.06 0.94;
54 2 113 32 0 0 1 0.955 15.26 138 1 1.06 0.94;
55 2 63 22 0 0 1 0.952 14.97 138 1 1.06 0.94;
56 2 84 18 0 0 1 0.954 15.16 138 1 1.06 0.94;
57 1 12 3 0 0 1 0.971 16.36 138 1 1.06 0.94;
58 1 12 3 0 0 1 0.959 15.51 138 1 1.06 0.94;
59 2 277 113 0 0 1 0.985 19.37 138 1 1.06 0.94;
60 1 78 3 0 0 1 0.993 23.15 138 1 1.06 0.94;
61 2 0 0 0 0 1 0.995 24.04 138 1 1.06 0.94;
62 2 77 14 0 0 1 0.998 23.43 138 1 1.06 0.94;
63 1 0 0 0 0 1 0.969 22.75 345 1 1.06 0.94;
64 1 0 0 0 0 1 0.984 24.52 345 1 1.06 0.94;
65 2 0 0 0 0 1 1.005 27.65 345 1 1.06 0.94;
66 2 39 18 0 0 1 1.05 27.48 138 1 1.06 0.94;
67 1 28 7 0 0 1 1.02 24.84 138 1 1.06 0.94;
68 1 0 0 0 0 1 1.003 27.55 345 1 1.06 0.94;
69 3 0 0 0 0 1 1.035 30 138 1 1.06 0.94;
70 2 66 20 0 0 1 0.984 22.58 138 1 1.06 0.94;
71 1 0 0 0 0 1 0.987 22.15 138 1 1.06 0.94;
72 2 12 0 0 0 1 0.98 20.98 138 1 1.06 0.94;
73 2 6 0 0 0 1 0.991 21.94 138 1 1.06 0.94;
74 2 68 27 0 12 1 0.958 21.64 138 1 1.06 0.94;
75 1 47 11 0 0 1 0.967 22.91 138 1 1.06 0.94;
76 2 68 36 0 0 1 0.943 21.77 138 1 1.06 0.94;
77 2 61 28 0 0 1 1.006 26.72 138 1 1.06 0.94;
78 1 71 26 0 0 1 1.003 26.42 138 1 1.06 0.94;
79 1 39 32 0 20 1 1.009 26.72 138 1 1.06 0.94;
80 2 130 26 0 0 1 1.04 28.96 138 1 1.06 0.94;
81 1 0 0 0 0 1 0.997 28.1 345 1 1.06 0.94;
82 1 54 27 0 20 1 0.989 27.24 138 1 1.06 0.94;
83 1 20 10 0 10 1 0.985 28.42 138 1 1.06 0.94;
84 1 11 7 0 0 1 0.98 30.95 138 1 1.06 0.94;
85 2 24 15 0 0 1 0.985 32.51 138 1 1.06 0.94;
86 1 21 10 0 0 1 0.987 31.14 138 1 1.06 0.94;
87 2 0 0 0 0 1 1.015 31.4 161 1 1.06 0.94;
88 1 48 10 0 0 1 0.987 35.64 138 1 1.06 0.94;
89 2 0 0 0 0 1 1.005 39.69 138 1 1.06 0.94;
90 2 163 42 0 0 1 0.985 33.29 138 1 1.06 0.94;
91 2 10 0 0 0 1 0.98 33.31 138 1 1.06 0.94;
92 2 65 10 0 0 1 0.993 33.8 138 1 1.06 0.94;
93 1 12 7 0 0 1 0.987 30.79 138 1 1.06 0.94;
94 1 30 16 0 0 1 0.991 28.64 138 1 1.06 0.94;
95 1 42 31 0 0 1 0.981 27.67 138 1 1.06 0.94;
96 1 38 15 0 0 1 0.993 27.51 138 1 1.06 0.94;
97 1 15 9 0 0 1 1.011 27.88 138 1 1.06 0.94;
98 1 34 8 0 0 1 1.024 27.4 138 1 1.06 0.94;
99 2 42 0 0 0 1 1.01 27.04 138 1 1.06 0.94;
100 2 37 18 0 0 1 1.017 28.03 138 1 1.06 0.94;
101 1 22 15 0 0 1 0.993 29.61 138 1 1.06 0.94;
102 1 5 3 0 0 1 0.991 32.3 138 1 1.06 0.94;
103 2 23 16 0 0 1 1.001 24.44 138 1 1.06 0.94;
104 2 38 25 0 0 1 0.971 21.69 138 1 1.06 0.94;
105 2 31 26 0 20 1 0.965 20.57 138 1 1.06 0.94;
106 1 43 16 0 0 1 0.962 20.32 138 1 1.06 0.94;
107 2 50 12 0 6 1 0.952 17.53 138 1 1.06 0.94;
108 1 2 1 0 0 1 0.967 19.38 138 1 1.06 0.94;
109 1 8 3 0 0 1 0.967 18.93 138 1 1.06 0.94;
110 2 39 30 0 6 1 0.973 18.09 138 1 1.06 0.94;
111 2 0 0 0 0 1 0.98 19.74 138 1 1.06 0.94;
112 2 68 13 0 0 1 0.975 14.99 138 1 1.06 0.94;
113 2 6 0 0 0 1 0.993 13.74 138 1 1.06 0.94;
114 1 8 3 0 0 1 0.96 14.46 138 1 1.06 0.94;
115 1 22 7 0 0 1 0.96 14.46 138 1 1.06 0.94;
116 2 184 0 0 0 1 1.005 27.12 138 1 1.06 0.94;
117 1 20 8 0 0 1 0.974 10.67 138 1 1.06 0.94;
118 1 33 15 0 0 1 0.949 21.92 138 1 1.06 0.94;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
gen = [
1 0 0 15 -5 0.955 100 1 100 0;
4 0 0 300 -300 0.998 100 1 100 0;
6 0 0 50 -13 0.99 100 1 100 0;
8 0 0 300 -300 1.015 100 1 100 0;
10 450 0 200 -147 1.05 100 1 550 0;
12 85 0 120 -35 0.99 100 1 185 0;
15 0 0 30 -10 0.97 100 1 100 0;
18 0 0 50 -16 0.973 100 1 100 0;
19 0 0 24 -8 0.962 100 1 100 0;
24 0 0 300 -300 0.992 100 1 100 0;
25 220 0 140 -47 1.05 100 1 320 0;
26 314 0 1000 -1000 1.015 100 1 414 0;
27 0 0 300 -300 0.968 100 1 100 0;
31 7 0 300 -300 0.967 100 1 107 0;
32 0 0 42 -14 0.963 100 1 100 0;
34 0 0 24 -8 0.984 100 1 100 0;
36 0 0 24 -8 0.98 100 1 100 0;
40 0 0 300 -300 0.97 100 1 100 0;
42 0 0 300 -300 0.985 100 1 100 0;
46 19 0 100 -100 1.005 100 1 119 0;
49 204 0 210 -85 1.025 100 1 304 0;
54 48 0 300 -300 0.955 100 1 148 0;
55 0 0 23 -8 0.952 100 1 100 0;
56 0 0 15 -8 0.954 100 1 100 0;
59 155 0 180 -60 0.985 100 1 255 0;
61 160 0 300 -100 0.995 100 1 260 0;
62 0 0 20 -20 0.998 100 1 100 0;
65 391 0 200 -67 1.005 100 1 491 0;
66 392 0 200 -67 1.05 100 1 492 0;
69 516.4 0 300 -300 1.035 100 1 805.2 0;
70 0 0 32 -10 0.984 100 1 100 0;
72 0 0 100 -100 0.98 100 1 100 0;
73 0 0 100 -100 0.991 100 1 100 0;
74 0 0 9 -6 0.958 100 1 100 0;
76 0 0 23 -8 0.943 100 1 100 0;
77 0 0 70 -20 1.006 100 1 100 0;
80 477 0 280 -165 1.04 100 1 577 0;
85 0 0 23 -8 0.985 100 1 100 0;
87 4 0 1000 -100 1.015 100 1 104 0;
89 607 0 300 -210 1.005 100 1 707 0;
90 0 0 300 -300 0.985 100 1 100 0;
91 0 0 100 -100 0.98 100 1 100 0;
92 0 0 9 -3 0.99 100 1 100 0;
99 0 0 100 -100 1.01 100 1 100 0;
100 252 0 155 -50 1.017 100 1 352 0;
103 40 0 40 -15 1.01 100 1 140 0;
104 0 0 23 -8 0.971 100 1 100 0;
105 0 0 23 -8 0.965 100 1 100 0;
107 0 0 200 -200 0.952 100 1 100 0;
110 0 0 23 -8 0.973 100 1 100 0;
111 36 0 1000 -100 0.98 100 1 136 0;
112 0 0 1000 -100 0.975 100 1 100 0;
113 0 0 200 -100 0.993 100 1 100 0;
116 0 0 1000 -1000 1.005 100 1 100 0;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status
branch = [
1 2 0.0303 0.0999 0.0254 9900 0 0 0 0 1;
1 3 0.0129 0.0424 0.01082 9900 0 0 0 0 1;
4 5 0.00176 0.00798 0.0021 9900 0 0 0 0 1;
3 5 0.0241 0.108 0.0284 9900 0 0 0 0 1;
5 6 0.0119 0.054 0.01426 9900 0 0 0 0 1;
6 7 0.00459 0.0208 0.0055 9900 0 0 0 0 1;
8 9 0.00244 0.0305 1.162 9900 0 0 0 0 1;
8 5 0 0.0267 0 9900 0 0 0.985 0 1;
9 10 0.00258 0.0322 1.23 9900 0 0 0 0 1;
4 11 0.0209 0.0688 0.01748 9900 0 0 0 0 1;
5 11 0.0203 0.0682 0.01738 9900 0 0 0 0 1;
11 12 0.00595 0.0196 0.00502 9900 0 0 0 0 1;
2 12 0.0187 0.0616 0.01572 9900 0 0 0 0 1;
3 12 0.0484 0.16 0.0406 9900 0 0 0 0 1;
7 12 0.00862 0.034 0.00874 9900 0 0 0 0 1;
11 13 0.02225 0.0731 0.01876 9900 0 0 0 0 1;
12 14 0.0215 0.0707 0.01816 9900 0 0 0 0 1;
13 15 0.0744 0.2444 0.06268 9900 0 0 0 0 1;
14 15 0.0595 0.195 0.0502 9900 0 0 0 0 1;
12 16 0.0212 0.0834 0.0214 9900 0 0 0 0 1;
15 17 0.0132 0.0437 0.0444 9900 0 0 0 0 1;
16 17 0.0454 0.1801 0.0466 9900 0 0 0 0 1;
17 18 0.0123 0.0505 0.01298 9900 0 0 0 0 1;
18 19 0.01119 0.0493 0.01142 9900 0 0 0 0 1;
19 20 0.0252 0.117 0.0298 9900 0 0 0 0 1;
15 19 0.012 0.0394 0.0101 9900 0 0 0 0 1;
20 21 0.0183 0.0849 0.0216 9900 0 0 0 0 1;
21 22 0.0209 0.097 0.0246 9900 0 0 0 0 1;
22 23 0.0342 0.159 0.0404 9900 0 0 0 0 1;
23 24 0.0135 0.0492 0.0498 9900 0 0 0 0 1;
23 25 0.0156 0.08 0.0864 9900 0 0 0 0 1;
26 25 0 0.0382 0 9900 0 0 0.96 0 1;
25 27 0.0318 0.163 0.1764 9900 0 0 0 0 1;
27 28 0.01913 0.0855 0.0216 9900 0 0 0 0 1;
28 29 0.0237 0.0943 0.0238 9900 0 0 0 0 1;
30 17 0 0.0388 0 9900 0 0 0.96 0 1;
8 30 0.00431 0.0504 0.514 9900 0 0 0 0 1;
26 30 0.00799 0.086 0.908 9900 0 0 0 0 1;
17 31 0.0474 0.1563 0.0399 9900 0 0 0 0 1;
29 31 0.0108 0.0331 0.0083 9900 0 0 0 0 1;
23 32 0.0317 0.1153 0.1173 9900 0 0 0 0 1;
31 32 0.0298 0.0985 0.0251 9900 0 0 0 0 1;
27 32 0.0229 0.0755 0.01926 9900 0 0 0 0 1;
15 33 0.038 0.1244 0.03194 9900 0 0 0 0 1;
19 34 0.0752 0.247 0.0632 9900 0 0 0 0 1;
35 36 0.00224 0.0102 0.00268 9900 0 0 0 0 1;
35 37 0.011 0.0497 0.01318 9900 0 0 0 0 1;
33 37 0.0415 0.142 0.0366 9900 0 0 0 0 1;
34 36 0.00871 0.0268 0.00568 9900 0 0 0 0 1;
34 37 0.00256 0.0094 0.00984 9900 0 0 0 0 1;
38 37 0 0.0375 0 9900 0 0 0.935 0 1;
37 39 0.0321 0.106 0.027 9900 0 0 0 0 1;
37 40 0.0593 0.168 0.042 9900 0 0 0 0 1;
30 38 0.00464 0.054 0.422 9900 0 0 0 0 1;
39 40 0.0184 0.0605 0.01552 9900 0 0 0 0 1;
40 41 0.0145 0.0487 0.01222 9900 0 0 0 0 1;
40 42 0.0555 0.183 0.0466 9900 0 0 0 0 1;
41 42 0.041 0.135 0.0344 9900 0 0 0 0 1;
43 44 0.0608 0.2454 0.06068 9900 0 0 0 0 1;
34 43 0.0413 0.1681 0.04226 9900 0 0 0 0 1;
44 45 0.0224 0.0901 0.0224 9900 0 0 0 0 1;
45 46 0.04 0.1356 0.0332 9900 0 0 0 0 1;
46 47 0.038 0.127 0.0316 9900 0 0 0 0 1;
46 48 0.0601 0.189 0.0472 9900 0 0 0 0 1;
47 49 0.0191 0.0625 0.01604 9900 0 0 0 0 1;
42 49 0.0715 0.323 0.086 9900 0 0 0 0 1;
42 49 0.0715 0.323 0.086 9900 0 0 0 0 1;
45 49 0.0684 0.186 0.0444 9900 0 0 0 0 1;
48 49 0.0179 0.0505 0.01258 9900 0 0 0 0 1;
49 50 0.0267 0.0752 0.01874 9900 0 0 0 0 1;
49 51 0.0486 0.137 0.0342 9900 0 0 0 0 1;
51 52 0.0203 0.0588 0.01396 9900 0 0 0 0 1;
52 53 0.0405 0.1635 0.04058 9900 0 0 0 0 1;
53 54 0.0263 0.122 0.031 9900 0 0 0 0 1;
49 54 0.073 0.289 0.0738 9900 0 0 0 0 1;
49 54 0.0869 0.291 0.073 9900 0 0 0 0 1;
54 55 0.0169 0.0707 0.0202 9900 0 0 0 0 1;
54 56 0.00275 0.00955 0.00732 9900 0 0 0 0 1;
55 56 0.00488 0.0151 0.00374 9900 0 0 0 0 1;
56 57 0.0343 0.0966 0.0242 9900 0 0 0 0 1;
50 57 0.0474 0.134 0.0332 9900 0 0 0 0 1;
56 58 0.0343 0.0966 0.0242 9900 0 0 0 0 1;
51 58 0.0255 0.0719 0.01788 9900 0 0 0 0 1;
54 59 0.0503 0.2293 0.0598 9900 0 0 0 0 1;
56 59 0.0825 0.251 0.0569 9900 0 0 0 0 1;
56 59 0.0803 0.239 0.0536 9900 0 0 0 0 1;
55 59 0.04739 0.2158 0.05646 9900 0 0 0 0 1;
59 60 0.0317 0.145 0.0376 9900 0 0 0 0 1;
59 61 0.0328 0.15 0.0388 9900 0 0 0 0 1;
60 61 0.00264 0.0135 0.01456 9900 0 0 0 0 1;
60 62 0.0123 0.0561 0.01468 9900 0 0 0 0 1;
61 62 0.00824 0.0376 0.0098 9900 0 0 0 0 1;
63 59 0 0.0386 0 9900 0 0 0.96 0 1;
63 64 0.00172 0.02 0.216 9900 0 0 0 0 1;
64 61 0 0.0268 0 9900 0 0 0.985 0 1;
38 65 0.00901 0.0986 1.046 9900 0 0 0 0 1;
64 65 0.00269 0.0302 0.38 9900 0 0 0 0 1;
49 66 0.018 0.0919 0.0248 9900 0 0 0 0 1;
49 66 0.018 0.0919 0.0248 9900 0 0 0 0 1;
62 66 0.0482 0.218 0.0578 9900 0 0 0 0 1;
62 67 0.0258 0.117 0.031 9900 0 0 0 0 1;
65 66 0 0.037 0 9900 0 0 0.935 0 1;
66 67 0.0224 0.1015 0.02682 9900 0 0 0 0 1;
65 68 0.00138 0.016 0.638 9900 0 0 0 0 1;
47 69 0.0844 0.2778 0.07092 9900 0 0 0 0 1;
49 69 0.0985 0.324 0.0828 9900 0 0 0 0 1;
68 69 0 0.037 0 9900 0 0 0.935 0 1;
69 70 0.03 0.127 0.122 9900 0 0 0 0 1;
24 70 0.00221 0.4115 0.10198 9900 0 0 0 0 1;
70 71 0.00882 0.0355 0.00878 9900 0 0 0 0 1;
24 72 0.0488 0.196 0.0488 9900 0 0 0 0 1;
71 72 0.0446 0.18 0.04444 9900 0 0 0 0 1;
71 73 0.00866 0.0454 0.01178 9900 0 0 0 0 1;
70 74 0.0401 0.1323 0.03368 9900 0 0 0 0 1;
70 75 0.0428 0.141 0.036 9900 0 0 0 0 1;
69 75 0.0405 0.122 0.124 9900 0 0 0 0 1;
74 75 0.0123 0.0406 0.01034 9900 0 0 0 0 1;
76 77 0.0444 0.148 0.0368 9900 0 0 0 0 1;
69 77 0.0309 0.101 0.1038 9900 0 0 0 0 1;
75 77 0.0601 0.1999 0.04978 9900 0 0 0 0 1;
77 78 0.00376 0.0124 0.01264 9900 0 0 0 0 1;
78 79 0.00546 0.0244 0.00648 9900 0 0 0 0 1;
77 80 0.017 0.0485 0.0472 9900 0 0 0 0 1;
77 80 0.0294 0.105 0.0228 9900 0 0 0 0 1;
79 80 0.0156 0.0704 0.0187 9900 0 0 0 0 1;
68 81 0.00175 0.0202 0.808 9900 0 0 0 0 1;
81 80 0 0.037 0 9900 0 0 0.935 0 1;
77 82 0.0298 0.0853 0.08174 9900 0 0 0 0 1;
82 83 0.0112 0.03665 0.03796 9900 0 0 0 0 1;
83 84 0.0625 0.132 0.0258 9900 0 0 0 0 1;
83 85 0.043 0.148 0.0348 9900 0 0 0 0 1;
84 85 0.0302 0.0641 0.01234 9900 0 0 0 0 1;
85 86 0.035 0.123 0.0276 9900 0 0 0 0 1;
86 87 0.02828 0.2074 0.0445 9900 0 0 0 0 1;
85 88 0.02 0.102 0.0276 9900 0 0 0 0 1;
85 89 0.0239 0.173 0.047 9900 0 0 0 0 1;
88 89 0.0139 0.0712 0.01934 9900 0 0 0 0 1;
89 90 0.0518 0.188 0.0528 9900 0 0 0 0 1;
89 90 0.0238 0.0997 0.106 9900 0 0 0 0 1;
90 91 0.0254 0.0836 0.0214 9900 0 0 0 0 1;
89 92 0.0099 0.0505 0.0548 9900 0 0 0 0 1;
89 92 0.0393 0.1581 0.0414 9900 0 0 0 0 1;
91 92 0.0387 0.1272 0.03268 9900 0 0 0 0 1;
92 93 0.0258 0.0848 0.0218 9900 0 0 0 0 1;
92 94 0.0481 0.158 0.0406 9900 0 0 0 0 1;
93 94 0.0223 0.0732 0.01876 9900 0 0 0 0 1;
94 95 0.0132 0.0434 0.0111 9900 0 0 0 0 1;
80 96 0.0356 0.182 0.0494 9900 0 0 0 0 1;
82 96 0.0162 0.053 0.0544 9900 0 0 0 0 1;
94 96 0.0269 0.0869 0.023 9900 0 0 0 0 1;
80 97 0.0183 0.0934 0.0254 9900 0 0 0 0 1;
80 98 0.0238 0.108 0.0286 9900 0 0 0 0 1;
80 99 0.0454 0.206 0.0546 9900 0 0 0 0 1;
92 100 0.0648 0.295 0.0472 9900 0 0 0 0 1;
94 100 0.0178 0.058 0.0604 9900 0 0 0 0 1;
95 96 0.0171 0.0547 0.01474 9900 0 0 0 0 1;
96 97 0.0173 0.0885 0.024 9900 0 0 0 0 1;
98 100 0.0397 0.179 0.0476 9900 0 0 0 0 1;
99 100 0.018 0.0813 0.0216 9900 0 0 0 0 1;
100 101 0.0277 0.1262 0.0328 9900 0 0 0 0 1;
92 102 0.0123 0.0559 0.01464 9900 0 0 0 0 1;
101 102 0.0246 0.112 0.0294 9900 0 0 0 0 1;
100 103 0.016 0.0525 0.0536 9900 0 0 0 0 1;
100 104 0.0451 0.204 0.0541 9900 0 0 0 0 1;
103 104 0.0466 0.1584 0.0407 9900 0 0 0 0 1;
103 105 0.0535 0.1625 0.0408 9900 0 0 0 0 1;
100 106 0.0605 0.229 0.062 9900 0 0 0 0 1;
104 105 0.00994 0.0378 0.00986 9900 0 0 0 0 1;
105 106 0.014 0.0547 0.01434 9900 0 0 0 0 1;
105 107 0.053 0.183 0.0472 9900 0 0 0 0 1;
105 108 0.0261 0.0703 0.01844 9900 0 0 0 0 1;
106 107 0.053 0.183 0.0472 9900 0 0 0 0 1;
108 109 0.0105 0.0288 0.0076 9900 0 0 0 0 1;
103 110 0.03906 0.1813 0.0461 9900 0 0 0 0 1;
109 110 0.0278 0.0762 0.0202 9900 0 0 0 0 1;
110 111 0.022 0.0755 0.02 9900 0 0 0 0 1;
110 112 0.0247 0.064 0.062 9900 0 0 0 0 1;
17 113 0.00913 0.0301 0.00768 9900 0 0 0 0 1;
32 113 0.0615 0.203 0.0518 9900 0 0 0 0 1;
32 114 0.0135 0.0612 0.01628 9900 0 0 0 0 1;
27 115 0.0164 0.0741 0.01972 9900 0 0 0 0 1;
114 115 0.0023 0.0104 0.00276 9900 0 0 0 0 1;
68 116 0.00034 0.00405 0.164 9900 0 0 0 0 1;
12 117 0.0329 0.14 0.0358 9900 0 0 0 0 1;
75 118 0.0145 0.0481 0.01198 9900 0 0 0 0 1;
76 118 0.0164 0.0544 0.01356 9900 0 0 0 0 1;
];
%%—– OPF Data —–%%
%% area data
areas = [
1 1;
];
%% generator cost data
% 1 startup shutdown n x1 y1 … xn yn
% 2 startup shutdown n c(n-1) … c0
gencost = [
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.0222222 20 0;
2 0 0 3 0.117647 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.0454545 20 0;
2 0 0 3 0.0318471 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 1.42857 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.526316 20 0;
2 0 0 3 0.0490196 20 0;
2 0 0 3 0.208333 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.0645161 20 0;
2 0 0 3 0.0625 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.0255754 20 0;
2 0 0 3 0.0255102 20 0;
2 0 0 3 0.0193648 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.0209644 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 2.5 20 0;
2 0 0 3 0.0164745 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.0396825 20 0;
2 0 0 3 0.25 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.277778 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
];
return;
% Warnings from cdf2matp conversion:
%
% ***** area data conversion not yet implemented (creating dummy area data)
% ***** negative Pg at bus 4 treated as Pd
% ***** negative Pg at bus 8 treated as Pd
% ***** negative Pg at bus 24 treated as Pd
% ***** negative Pg at bus 27 treated as Pd
% ***** negative Pg at bus 40 treated as Pd
% ***** negative Pg at bus 42 treated as Pd
% ***** negative Pg at bus 72 treated as Pd
% ***** negative Pg at bus 73 treated as Pd
% ***** negative Pg at bus 90 treated as Pd
% ***** negative Pg at bus 91 treated as Pd
% ***** negative Pg at bus 99 treated as Pd
% ***** negative Pg at bus 107 treated as Pd
% ***** negative Pg at bus 112 treated as Pd
% ***** negative Pg at bus 113 treated as Pd
% ***** negative Pg at bus 116 treated as Pd
% ***** Insufficient generation, setting Pmax at slack bus (bus 69) to 805.2
% ***** MVA limit of branch 1 – 2 not given, set to 9900
% ***** MVA limit of branch 1 – 3 not given, set to 9900
% ***** MVA limit of branch 4 – 5 not given, set to 9900
% ***** MVA limit of branch 3 – 5 not given, set to 9900
% ***** MVA limit of branch 5 – 6 not given, set to 9900
% ***** MVA limit of branch 6 – 7 not given, set to 9900
% ***** MVA limit of branch 8 – 9 not given, set to 9900
% ***** MVA limit of branch 8 – 5 not given, set to 9900
% ***** MVA limit of branch 9 – 10 not given, set to 9900
% ***** MVA limit of branch 4 – 11 not given, set to 9900
% ***** MVA limit of branch 5 – 11 not given, set to 9900
% ***** MVA limit of branch 11 – 12 not given, set to 9900
% ***** MVA limit of branch 2 – 12 not given, set to 9900
% ***** MVA limit of branch 3 – 12 not given, set to 9900
% ***** MVA limit of branch 7 – 12 not given, set to 9900
% ***** MVA limit of branch 11 – 13 not given, set to 9900
% ***** MVA limit of branch 12 – 14 not given, set to 9900
% ***** MVA limit of branch 13 – 15 not given, set to 9900
% ***** MVA limit of branch 14 – 15 not given, set to 9900
% ***** MVA limit of branch 12 – 16 not given, set to 9900
% ***** MVA limit of branch 15 – 17 not given, set to 9900
% ***** MVA limit of branch 16 – 17 not given, set to 9900
% ***** MVA limit of branch 17 – 18 not given, set to 9900
% ***** MVA limit of branch 18 – 19 not given, set to 9900
% ***** MVA limit of branch 19 – 20 not given, set to 9900
% ***** MVA limit of branch 15 – 19 not given, set to 9900
% ***** MVA limit of branch 20 – 21 not given, set to 9900
% ***** MVA limit of branch 21 – 22 not given, set to 9900
% ***** MVA limit of branch 22 – 23 not given, set to 9900
% ***** MVA limit of branch 23 – 24 not given, set to 9900
% ***** MVA limit of branch 23 – 25 not given, set to 9900
% ***** MVA limit of branch 26 – 25 not given, set to 9900
% ***** MVA limit of branch 25 – 27 not given, set to 9900
% ***** MVA limit of branch 27 – 28 not given, set to 9900
% ***** MVA limit of branch 28 – 29 not given, set to 9900
% ***** MVA limit of branch 30 – 17 not given, set to 9900
% ***** MVA limit of branch 8 – 30 not given, set to 9900
% ***** MVA limit of branch 26 – 30 not given, set to 9900
% ***** MVA limit of branch 17 – 31 not given, set to 9900
% ***** MVA limit of branch 29 – 31 not given, set to 9900
% ***** MVA limit of branch 23 – 32 not given, set to 9900
% ***** MVA limit of branch 31 – 32 not given, set to 9900
% ***** MVA limit of branch 27 – 32 not given, set to 9900
% ***** MVA limit of branch 15 – 33 not given, set to 9900
% ***** MVA limit of branch 19 – 34 not given, set to 9900
% ***** MVA limit of branch 35 – 36 not given, set to 9900
% ***** MVA limit of branch 35 – 37 not given, set to 9900
% ***** MVA limit of branch 33 – 37 not given, set to 9900
% ***** MVA limit of branch 34 – 36 not given, set to 9900
% ***** MVA limit of branch 34 – 37 not given, set to 9900
% ***** MVA limit of branch 38 – 37 not given, set to 9900
% ***** MVA limit of branch 37 – 39 not given, set to 9900
% ***** MVA limit of branch 37 – 40 not given, set to 9900
% ***** MVA limit of branch 30 – 38 not given, set to 9900
% ***** MVA limit of branch 39 – 40 not given, set to 9900
% ***** MVA limit of branch 40 – 41 not given, set to 9900
% ***** MVA limit of branch 40 – 42 not given, set to 9900
% ***** MVA limit of branch 41 – 42 not given, set to 9900
% ***** MVA limit of branch 43 – 44 not given, set to 9900
% ***** MVA limit of branch 34 – 43 not given, set to 9900
% ***** MVA limit of branch 44 – 45 not given, set to 9900
% ***** MVA limit of branch 45 – 46 not given, set to 9900
% ***** MVA limit of branch 46 – 47 not given, set to 9900
% ***** MVA limit of branch 46 – 48 not given, set to 9900
% ***** MVA limit of branch 47 – 49 not given, set to 9900
% ***** MVA limit of branch 42 – 49 not given, set to 9900
% ***** MVA limit of branch 42 – 49 not given, set to 9900
% ***** MVA limit of branch 45 – 49 not given, set to 9900
% ***** MVA limit of branch 48 – 49 not given, set to 9900
% ***** MVA limit of branch 49 – 50 not given, set to 9900
% ***** MVA limit of branch 49 – 51 not given, set to 9900
% ***** MVA limit of branch 51 – 52 not given, set to 9900
% ***** MVA limit of branch 52 – 53 not given, set to 9900
% ***** MVA limit of branch 53 – 54 not given, set to 9900
% ***** MVA limit of branch 49 – 54 not given, set to 9900
% ***** MVA limit of branch 49 – 54 not given, set to 9900
% ***** MVA limit of branch 54 – 55 not given, set to 9900
% ***** MVA limit of branch 54 – 56 not given, set to 9900
% ***** MVA limit of branch 55 – 56 not given, set to 9900
% ***** MVA limit of branch 56 – 57 not given, set to 9900
% ***** MVA limit of branch 50 – 57 not given, set to 9900
% ***** MVA limit of branch 56 – 58 not given, set to 9900
% ***** MVA limit of branch 51 – 58 not given, set to 9900
% ***** MVA limit of branch 54 – 59 not given, set to 9900
% ***** MVA limit of branch 56 – 59 not given, set to 9900
% ***** MVA limit of branch 56 – 59 not given, set to 9900
% ***** MVA limit of branch 55 – 59 not given, set to 9900
% ***** MVA limit of branch 59 – 60 not given, set to 9900
% ***** MVA limit of branch 59 – 61 not given, set to 9900
% ***** MVA limit of branch 60 – 61 not given, set to 9900
% ***** MVA limit of branch 60 – 62 not given, set to 9900
% ***** MVA limit of branch 61 – 62 not given, set to 9900
% ***** MVA limit of branch 63 – 59 not given, set to 9900
% ***** MVA limit of branch 63 – 64 not given, set to 9900
% ***** MVA limit of branch 64 – 61 not given, set to 9900
% ***** MVA limit of branch 38 – 65 not given, set to 9900
% ***** MVA limit of branch 64 – 65 not given, set to 9900
% ***** MVA limit of branch 49 – 66 not given, set to 9900
% ***** MVA limit of branch 49 – 66 not given, set to 9900
% ***** MVA limit of branch 62 – 66 not given, set to 9900
% ***** MVA limit of branch 62 – 67 not given, set to 9900
% ***** MVA limit of branch 65 – 66 not given, set to 9900
% ***** MVA limit of branch 66 – 67 not given, set to 9900
% ***** MVA limit of branch 65 – 68 not given, set to 9900
% ***** MVA limit of branch 47 – 69 not given, set to 9900
% ***** MVA limit of branch 49 – 69 not given, set to 9900
% ***** MVA limit of branch 68 – 69 not given, set to 9900
% ***** MVA limit of branch 69 – 70 not given, set to 9900
% ***** MVA limit of branch 24 – 70 not given, set to 9900
% ***** MVA limit of branch 70 – 71 not given, set to 9900
% ***** MVA limit of branch 24 – 72 not given, set to 9900
% ***** MVA limit of branch 71 – 72 not given, set to 9900
% ***** MVA limit of branch 71 – 73 not given, set to 9900
% ***** MVA limit of branch 70 – 74 not given, set to 9900
% ***** MVA limit of branch 70 – 75 not given, set to 9900
% ***** MVA limit of branch 69 – 75 not given, set to 9900
% ***** MVA limit of branch 74 – 75 not given, set to 9900
% ***** MVA limit of branch 76 – 77 not given, set to 9900
% ***** MVA limit of branch 69 – 77 not given, set to 9900
% ***** MVA limit of branch 75 – 77 not given, set to 9900
% ***** MVA limit of branch 77 – 78 not given, set to 9900
% ***** MVA limit of branch 78 – 79 not given, set to 9900
% ***** MVA limit of branch 77 – 80 not given, set to 9900
% ***** MVA limit of branch 77 – 80 not given, set to 9900
% ***** MVA limit of branch 79 – 80 not given, set to 9900
% ***** MVA limit of branch 68 – 81 not given, set to 9900
% ***** MVA limit of branch 81 – 80 not given, set to 9900
% ***** MVA limit of branch 77 – 82 not given, set to 9900
% ***** MVA limit of branch 82 – 83 not given, set to 9900
% ***** MVA limit of branch 83 – 84 not given, set to 9900
% ***** MVA limit of branch 83 – 85 not given, set to 9900
% ***** MVA limit of branch 84 – 85 not given, set to 9900
% ***** MVA limit of branch 85 – 86 not given, set to 9900
% ***** MVA limit of branch 86 – 87 not given, set to 9900
% ***** MVA limit of branch 85 – 88 not given, set to 9900
% ***** MVA limit of branch 85 – 89 not given, set to 9900
% ***** MVA limit of branch 88 – 89 not given, set to 9900
% ***** MVA limit of branch 89 – 90 not given, set to 9900
% ***** MVA limit of branch 89 – 90 not given, set to 9900
% ***** MVA limit of branch 90 – 91 not given, set to 9900
% ***** MVA limit of branch 89 – 92 not given, set to 9900
% ***** MVA limit of branch 89 – 92 not given, set to 9900
% ***** MVA limit of branch 91 – 92 not given, set to 9900
% ***** MVA limit of branch 92 – 93 not given, set to 9900
% ***** MVA limit of branch 92 – 94 not given, set to 9900
% ***** MVA limit of branch 93 – 94 not given, set to 9900
% ***** MVA limit of branch 94 – 95 not given, set to 9900
% ***** MVA limit of branch 80 – 96 not given, set to 9900
% ***** MVA limit of branch 82 – 96 not given, set to 9900
% ***** MVA limit of branch 94 – 96 not given, set to 9900
% ***** MVA limit of branch 80 – 97 not given, set to 9900
% ***** MVA limit of branch 80 – 98 not given, set to 9900
% ***** MVA limit of branch 80 – 99 not given, set to 9900
% ***** MVA limit of branch 92 – 100 not given, set to 9900
% ***** MVA limit of branch 94 – 100 not given, set to 9900
% ***** MVA limit of branch 95 – 96 not given, set to 9900
% ***** MVA limit of branch 96 – 97 not given, set to 9900
% ***** MVA limit of branch 98 – 100 not given, set to 9900
% ***** MVA limit of branch 99 – 100 not given, set to 9900
% ***** MVA limit of branch 100 – 101 not given, set to 9900
% ***** MVA limit of branch 92 – 102 not given, set to 9900
% ***** MVA limit of branch 101 – 102 not given, set to 9900
% ***** MVA limit of branch 100 – 103 not given, set to 9900
% ***** MVA limit of branch 100 – 104 not given, set to 9900
% ***** MVA limit of branch 103 – 104 not given, set to 9900
% ***** MVA limit of branch 103 – 105 not given, set to 9900
% ***** MVA limit of branch 100 – 106 not given, set to 9900
% ***** MVA limit of branch 104 – 105 not given, set to 9900
% ***** MVA limit of branch 105 – 106 not given, set to 9900
% ***** MVA limit of branch 105 – 107 not given, set to 9900
% ***** MVA limit of branch 105 – 108 not given, set to 9900
% ***** MVA limit of branch 106 – 107 not given, set to 9900
% ***** MVA limit of branch 108 – 109 not given, set to 9900
% ***** MVA limit of branch 103 – 110 not given, set to 9900
% ***** MVA limit of branch 109 – 110 not given, set to 9900
% ***** MVA limit of branch 110 – 111 not given, set to 9900
% ***** MVA limit of branch 110 – 112 not given, set to 9900
% ***** MVA limit of branch 17 – 113 not given, set to 9900
% ***** MVA limit of branch 32 – 113 not given, set to 9900
% ***** MVA limit of branch 32 – 114 not given, set to 9900
% ***** MVA limit of branch 27 – 115 not given, set to 9900
% ***** MVA limit of branch 114 – 115 not given, set to 9900
% ***** MVA limit of branch 68 – 116 not given, set to 9900
% ***** MVA limit of branch 12 – 117 not given, set to 9900
% ***** MVA limit of branch 75 – 118 not given, set to 9900
% ***** MVA limit of branch 76 – 118 not given, set to 9900
dante // 15 April 2009 pada 18:06 |
makasih mas……….
tapi bisa ngasih cara bikin prorgamnya pake matlab mas…………………
dante // 20 April 2009 pada 13:40 |
oiya mas,,,,,
saya kan pemula, pemula pake matlab juga….
jadi ya lom begitu paham gtu……
mohon di ajari ya….
makasih………….
sosokmanusia // 22 Mei 2009 pada 14:33 |
Assalam..
Perkenalkan nama saya Radian…
saya sedang tugas akhir optimisasi daya reaktif..
mau nanya nih mas..
klo kita mau nyari batas kestabilan tegangan berdasarkan singularitas dari Jacobian gimana caranya??
bukannya klo dari load flow yd kt dapet udah hasilnya, gmn cara ngeluarin matrik Jacobiannya?
ada contohnya jika kita ingin menggunakan Matlab?? trims..
Muhammad Imaduddin // 23 Mei 2009 pada 13:12 |
Wa alaikum salaam Radian,
Untuk mencari batas kestabilan tegangan, biasanya kita menggunakan analisis statik :
1. V-Q Sensitivity Analysis
2. Q-V Modal Analysis
Dari load flow saja tidak cukup, karena kita tahu, pada saddle point / batas kestabilan tegangan matriks Jacobian menjadi singular.
Paling gampang, pakai saja toolbox Matlab yang namanya PST Power System Toolbox,
J. Chow, “Power System Toolbox Version 2.0,” Ontario, Canada: Cherry Tree Scientific Software, 2002
Sebelum menggunakannya, lebih baik dipahami dulu konsepnya, bacaan yang bagus di,
A general method for small signal stability analysis
Makarov, Yu.V.; Zhao Yang Dong; Hill, D.J.
Power Systems, IEEE Transactions on
Volume 13, Issue 3, Aug 1998 Page(s):979 – 985
Digital Object Identifier 10.1109/59.709086
Bifurcation theory and its application to nonlinear dynamical phenomena in an electrical power system Ajjarapu, V. Lee, B. Dept. of Electr. Eng. & Comput. Eng., Iowa State Univ., Ames, IA;
Saya beri clue, jika Radian sudah baca yang Modal Analysis, dan melakukan sensitivity analysis dengan PSTV2, coba perhatikan eigen values matriks Jacobian sistem yang diuji, kemana eigen values tsb bergerak menjelang voltage collapse.
sosokmanusia // 23 Mei 2009 pada 22:51 |
terima kasih y mas,…
saya pelajari dulu..
wassalam..
sosokmanusia // 24 Mei 2009 pada 00:02 |
Assalam mas..
boleh saya minta emailnya??
Kalo email saya radiannugraha@gmail.com
radian // 3 Juni 2009 pada 21:36 |
assalam mas..
say mau nanya, kalo kita mw melakukan aliran daya tapi matrik jacobiannya didapat dari persamaan daya yang dalam bentuk rectangular gmn caranya?? soalnya yang umum dipakai kan persamaan daya dalam bentuk polar (P cos x – Q sin x)..
terima kasih sebelumnya…
Muhammad Imaduddin // 8 Juni 2009 pada 15:16 |
Wassalam..
Bisanya sih bisa, namun Radian akan melakukan pekerjaan repetitif/berulang-ulang yang tidak perlu, misal jika kita ingin mengalikan (3 + j4) dengan (3 – j4) maka kita harus mengalikan 3 x 3, ditambah 3 x -j4, ditambah j4 x 3, ditambah j4 x -j4.
Panjang kan…
radian // 3 Juni 2009 pada 23:07 |
ohy mas..
ada ga contoh listing yang newton rhapsonnya pake matlab??
trim’s
Muhammad Imaduddin // 8 Juni 2009 pada 15:19 |
Contoh listing NR dapat dilihat di bagian ke-2 tulisan ini.
nur wahyu h. // 28 Juni 2009 pada 16:29 |
Assalamu’alaikum…
Pak Imaduddin, sy mau tanya:
1. Studi aliran daya itu kan mempelajari aliran daya dan tegangan sistem dalam kodisi operasi tunak. Kondisi operasi tunak itu definisinya apa sih? Sy masih blm paham.
2. Sy sudah mencoba menyimpan gambar perhitungan aliran daya metode newton raphson &fast decouple anda. Tapi ketika dicetak tulisannya pecah-pecah. Bisa minta tolong kirimkan perhitungan manual aliran dayanya ke e-mail saya: wah_u2@yahoo.co.id?
Terimakasih… Smga ilmu anda bermanfaat…!
Muhammad Imaduddin // 29 Juni 2009 pada 14:22 |
1. Kutipan dari Wiki:
steady state:
Any state that does not change over time, or is in a dynamic equilibrium
Keadaan tunak (bahasa Inggris: steady state) adalah kondisi sewaktu sifat-sifat suatu sistem tak berubah dengan berjalannya waktu atau dengan kata lain, konstan.
Artinya, Sistem Tenaga Listrik (STL) yang kita pelajari aliran dayanya, kita anggap tidak berubah secara drastis dengan berubahnya waktu. Perubahan konfigurasi sistem bisa terjadi, misalnya akibat:
-hilangnya beban secara tiba-tiba
-hilangnya suplai dari pembangkit secara tiba-tiba
-gangguan pada transmisi secara tiba-tiba
Yang perlu diingat, Power System Analysis (PSA) tidak cuma berisi Steady State Analysis (SSA). SSA cuma bagian kecil dari PSA.
Yang sering jadi pertanyaan, jika SSA cuma bagian kecil dari PSA, kenapa kita “cuma” mempelajari SSA di mata kuliah PSA?
Jawabannya:
1. SSA memberi gambaran umum (general knowledge) akan sebuah STL. Bahkan dengan bantuan software sekalipun, SSA lah yang pertama kali dilihat ketika merencanakan operasi STL.
2. Belajar SSA aja puyenk, gimana mo belajar Dynamic Analysis hehe… (maaf hanya guyon saja ko’). Intinya, kalau anda mau menggali lebih dalam lagi, ada banyak hal di studi ini yang sampai sekarang masih berkembang.
2. Tips untuk mencetak gambar2 di posting ini:
- Klik kanan gambar yang akan dicetak, pilih Open Link in a New Tab/Window
- Save gambar yang dibuka pada Tab/Window yang baru itu
- Cetak gambar dari file yang sudah anda simpan tadi
Jika di-save-as langsung dari posting ini, hasilnya memang jelek, hal ini saya sengaja, mempertimbangkan disain blog ini agar tidak terlalu dipenuhi oleh gambar yang terlalu besar.
Wira // 10 Juli 2009 pada 11:00 |
mohon bantuannya dong…ttg analisa aliran daya tiga fasa pada sistem distribusi dengan menggunakan metode fast decoupled. saya perlu perhitungan manualnya untuk membuat program dengan matlab. karena waktu kuliah semester lima dulu yang dipelajari hanya 1 fasa saja (diasumsikan semua fasa seimbang). terima kasih banget sebelumnya….
wuri // 13 Juli 2009 pada 23:40 |
salam kenal mas
saya kemarin ujian disuruh ngerjakan sama persis, tapi bus nya ada enam paralel. Kalo jadi matriks 8X8. 100 menit, manual.
Gila ga tuh
Skybers // 3 Agustus 2009 pada 17:50 |
mas…ada postingan teori tentang program yang ada di buku power system analysis g’mas???
makasih sebelumnya mas…
nur wahyu h. // 4 September 2009 pada 16:40 |
Assalamu’alaikum….
Sy mhsiswa t.elektro. Sy sedang mencari definisi tentang Heuristic dan tentang metode New Heuristic …
Mohon penjelasannya secara detail…
Penting untuk revisian skripsi sy rekonfigurasi jaringn distribsi .
trim’s atas bantuannya…
Muhammad Imaduddin // 6 September 2009 pada 10:41 |
Wa’alaikum salam,
Topik yang dicari pak Wahyu untuk skripsi sangat lah menarik. Saya memang berencana ingin membahasnya di blog saya ini, karena ada banyak hal yang masih akan berkembang di ranah ini.
Inti dari teknik algoritma ini adalah melakukan optimisasi dengan cara-cara yang baru yang menghasilkan solusi yang jauh lebih baik. Salah satu karakteristiknya adalah mendapatkan global optima dan terhindar dari jebakan solusi berupa local minima.
New Heuristic yang dimaksud mungkin Metaheuristic, yang menurut Blum dan Roli belum ada definisi yang baku. Khusus untuk topik rekonfigurasi jaringan distribusi, pak Wahyu bisa fokus saja di teknik Evolutionary Computation (termasuk di dalamnya ada Genetic Algorithm) yang inspirasinya diambil dari evolusi biologis.
Dalam jaringan distribusi existing seringkali kita dihadapkan pada kondisi klasik, bahwa kondisi sekarang dan akan datang, sudah tidak sesuai dengan perencanaan/master plan distribusi di masa lalu. Akibatnya bisa saja diperlukan pertukaran saluran, bahkan saluran yang ada dihapus diganti saluran lain. Itu dulu clue awalnya, pak Wahyu bisa search di Google mengenai dua hal ini, Evolutionary Computation dan Genetic Algorithm. Semoga bermanfaat.