Dear Prof Ray..
Case_info:
case_info(mpc)
Checking connectivity ... single fully connected network
Elapsed time is 0.107661 seconds.
================================================================================
Full
System
Number of: ----------
buses 31
loads 25
on 25
off -
fixed 25
dispatchable -
on -
off -
generators 11
on 11
off -
shunt elements 3
branches 38
on 38
off -
ties (off) -
Load
active (MW)
dispatched 13851.9
fixed 13851.9
dispatchable -
nominal 13851.9
on 13851.9
off -
fixed 13851.9
dispatchable -
on -
off -
reactive (MVAr)
dispatched 4300.1
fixed 4300.1
dispatchable -
nominal 4300.1
on 4300.1
off -
fixed 4300.1
dispatchable -
on -
off -
Generation
active (MW)
dispatched 15078.0
max capacity 18086.0
on 18086.0
off -
min capacity 9494.0
on 9494.0
off -
reactive (MVAr)
dispatched 2981.0
max capacity 6630.0
on 6630.0
off -
min capacity -3240.0
on -3240.0
off -
Shunt Injections
active (MW) -
reactive (MVAr) -391.8
Branch Losses
active (MW) -
reactive (MVAr) -
DC line
export (MW)
dispatch -
max capacity -
on -
off -
min capacity -
on -
off -
Reference Buses
num of ref buses 1
ref bus numbers 1
RegardsSabhan K
Pada Senin, 21 Oktober 2019 22.35.13 WIB, Sabhan Kanata
<[email protected]> menulis:
Dear Prof Ray..
How do I check with the case_info command, sir? Please the instructions
myfiledata : casejawabali.m
RegardsSabhan
Pada Senin, 21 Oktober 2019 22.23.34 WIB, Ray Zimmerman <[email protected]>
menulis:
Did you check the output of case_info() for your case?
Ray
> On Oct 21, 2019, at 10:04 AM, Sabhan Kanata <[email protected]> wrote:
>
> Dear Prof Ray,
>
> Thank you for responding
> My branch data, There are no R and X parameters that are 0. Where do you
> think the problem is?
>
> "
OutmailID: 124038130, List: 'matpower-l', MemberID: 28632318
SCRIPT: "branch data > % fbus tbus r x b rateA rateB
rateC ratio angle status > branch = [ > 1 2 0.000062650
0.0007008 0.0000000 9900 0 0 0 0 1; %Suralaya-Cilegon
> 1 24 0.000367768 0.0035333 0.0000000 9900 0 0 0
0 1; %Suralaya-Balaraja > 1 31 0.000006000 0.0005600
0.0000000 9900 0 0 0 0 1; %Suralaya-LBE > 31 24
0.000072160 0.0071762 0.0000000 9900 0 0 0 0 1;
%LBE-Balaraja > 2 5 0.001313332 0.0146926 0.2109700 9900 0
0 0 0 1; %Cilegon-Cibinong > 4 3 0.000151318
0.0016982 0.0000000 9900 0 0 0 0 1; %Gandul-Kembangan
> 18 4 0.000034700 0.0003334 0.0000000 9900 0 0 0
0 1; %Depok-Gandul > 7 5 0.000444188 0.0042675 0.0000000
9900 0 0 0 0 1; %Bekasi-Cibinong > 8 5 0.000621160
0.0059678 0.0000000 9900 0 0 0 0 1;
%Muaratawar-Cibinong > 11 5 0.000411138 0.0045995 0.2652900
9900 0 0 0 0 1; %Saguling-Cibinong > 6 7
0.000197365 0.0018962 0.0000000 9900 0 0 0 0 1;
%Cawang-Bekasi > 8 6 0.000562560 0.0054048 0.0000000 9900
0 0 0 0 1; %Muaratawar-Cawamg > 8 9 0.000282206
0.0027112 0.0000000 9900 0 0 0 0 1; %Muaratawar-Cibatu
> 10 9 0.000273900 0.0026324 0.0000000 9900 0 0 0
0 1; %Cirata-Cibatu > 11 10 0.000147473 0.0014168 0.0000000
9900 0 0 0 0 1; %Saguling-Cirata > 12 11
0.000195780 0.0021902 0.0000000 9900 0 0 0 0 1;
%Bandung Selatan-Saguling > 13 12 0.013981960 0.1343320
0.1924070 9900 0 0 0 0 1; %Mandiracan-Bandung Selatan >
13 26 0.001012900 0.0097320 0.1393940 9900 0 0 0
0 1; %Mandiracan-Ujung Berung > 26 12 0.000385400 0.0037034
0.0000000 9900 0 0 0 0 1; %Ujung Berung-Bandung Selatan >
30 13 0.000200000 0.0022500 0.0000000 9900 0 0 0
0 1; %Pemalang-Mandiracan > 14 30 0.000320000 0.0035300
0.0000000 9900 0 0 0 0 1; %Ungaran-Pemalang > 15 14
0.000067660 0.0075703 0.4366440 9900 0 0 0 0 1;
%Tanjung Djati-Ungaran > 16 14 0.002979200 0.0286229 0.4099740
9900 0 0 0 0 1; %Surabaya Barat-Ungaran > 14 20
0.000903612 0.0086815 0.0000000 9900 0 0 0 0 1;
%Ungaran-Pedan > 16 17 0.000140000 0.0013399 0.0000000 9900
0 0 0 0 1; %Surabaya Barat-Gresik > 23 16 0.000398638
0.0044597 0.0000000 9900 0 0 0 0 1; %Grati-Surabaya
Barat > 18 5 0.000091200 0.0008765 0.0000000 9900 0 0
0 0 1; %Depok-Cibinong > 19 18 0.001403000 0.0156967
0.9015550 9900 0 0 0 0 1; %Tasikmalaya-Depok > 29
19 0.000544600 0.0060933 0.3499780 9900 0 0 0 0
1; %Kesugihan-Tasikmalaya > 27 29 0.000110000 0.0012668
0.0000000 9900 0 0 0 0 1; %Adipala-Kesugihan > 28
27 0.000038100 0.0004387 0.0000000 9900 0 0 0 0
1; %Cilacap-Adipala > 20 29 0.000983900 0.0110073 0.6348830
9900 0 0 0 0 1; %Pedan-Kesugihan > 21 20
0.002053100 0.0229694 0.3298180 9900 0 0 0 0 1;
%Kediri-Pedan > 22 21 0.001026500 0.0114847 0.7239300 9900
0 0 0 0 1; %Paiton-Kediri > 22 23 0.000044100
0.0049343 0.3110350 9900 0 0 0 0 1; %Paiton-Grati >
24 4 0.000297922 0.0028623 0.0000000 9900 0 0 0
0 1; %Balaraja-Gandul > 25 14 0.002347961 0.0225581 0.3244670
9900 0 0 0 0 1; %Ngimbang-Ungaran > 16 25
0.000596665 0.0057324 0.0000000 9900 0 0 0 0 1;
%Surabaya Barat-Ngimbang > ]; > > > Regards > Sabhan K > Pada Senin, 21
Oktober 2019 20.55.34 WIB, Ray Zimmerman <[email protected]> menulis: > > > It
is possible that you have a disconnected network with an island that does not
have a reference bus, or maybe you have a zero impedance branch. I would begin
by looking at the output of case_info() for your case to see if you see
anything unusual. > > Ray > > >> On Oct 20, 2019, at 11:59 PM, Sabhan
Kanata <[email protected] <mailto:[email protected]>> wrote: >>
>> Dear Prof Ray, >> >> I get a warning when running the power flow as below:
>> >> What do you think is the solution? >> >> >> Error using \
>> Matrix is singular to working precision. >> >> Error in newtonpf (line 100)
>> >> dx = -(J \ F); >> >> Error in runpf (line 148) >> [V,
>> success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq, mpopt); >> >>
>> Regards >> Sabhan K > --Apple-Mail=_A33335EE-D265-4CC7-AB1C-A53F457D7D13
>> Content-Transfer-Encoding: quoted-printable Content-Type: text/html;
>> charset="us-ascii" <html><head><meta http-equiv="Content-Type"
>> content="text/html; charset=us-ascii"></head><body style="word-wrap:
>> break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"
>> class="">Did you check the output of <font face="Courier"
>> class="">case_info()</font> for your case?<div class=""><br
>> class=""></div><div class=""> Ray<br class=""><div><br
>> class=""><blockquote type="cite" class=""><div class="">On Oct 21, 2019, at
>> 10:04 AM, Sabhan Kanata <<a href="mailto:[email protected]"
>> class="">[email protected]</a>> wrote:</div><br
>> class="Apple-interchange-newline"><div class=""><div class=""><div
>> class="ydp26f789a1yahoo-style-wrap" style="font-family: Helvetica Neue,
>> Helvetica, Arial, sans-serif; font-size: 13px;"><div class=""></div>
>> <div dir="ltr" data-setdir="false" class="">Dear Prof Ray,</div><div
>> dir="ltr" data-setdir="false" class=""><br class=""></div><div dir="ltr"
>> data-setdir="false" class=""><span class="">Thank you for
>> responding</span><br class=""></div><div dir="ltr" data-setdir="false"
>> class="">My branch data, <span class="">There are no R and X parameters that
>> are 0. <span class="">Where do you think the problem
>> is?</span></span><br class=""></div><div dir="ltr" data-setdir="false"
>> class=""><br class=""></div><div dir="ltr" data-setdir="false" class=""><div
>> class="">"
TCL MERGE ERROR ( 10/21/2019 11:23:31 ): "extra characters after close-quote
branch data<br class="">% fbus
tbus r x
b rateA rateB
rateC ratio angle
status<br class="">branch = [<br class="">
1 2 0.000062650
0.0007008 0.0000000 9900
0 0 0
0 1; %Suralaya-Cilegon<br
class=""> 1 24
0.000367768 0.0035333
0.0000000 9900 0
0 0 0
1; %Suralaya-Balaraja<br class=""> 1
31 0.000006000 0.0005600
0.0000000 9900 0
0 0 0
1; %Suralaya-LBE<br class=""> 31
24 0.000072160 0.0071762
0.0000000 9900 0
0 0 0
1; %LBE-Balaraja<br class=""> 2
5 0.001313332 0.0146926
0.2109700 9900 0
0 0 0
1; %Cilegon-Cibinong<br class=""> 4
3 0.000151318 0.0016982
0.0000000 9900 0
0 0 0
1; %Gandul-Kembangan<br class=""> 18
4 0.000034700 0.0003334
0.0000000 9900 0
0 0 0
1; %Depok-Gandul<br class=""> 7
5 0.000444188 0.0042675
0.0000000 9900 0
0 0 0
1; %Bekasi-Cibinong<br class=""> 8
5 0.000621160 0.0059678
0.0000000 9900 0
0 0 0
1; %Muaratawar-Cibinong<br class=""> 11
5 0.000411138 0.0045995
0.2652900 9900 0
0 0 0
1; %Saguling-Cibinong<br class=""> 6
7 0.000197365 0.0018962
0.0000000 9900 0
0 0 0
1; %Cawang-Bekasi<br class=""> 8
6 0.000562560 0.0054048
0.0000000 9900 0
0 0 0
1; %Muaratawar-Cawamg<br class=""> 8
9 0.000282206 0.0027112
0.0000000 9900 0
0 0 0
1; %Muaratawar-Cibatu<br class=""> 10
9 0.000273900 0.0026324
0.0000000 9900 0
0 0 0
1; %Cirata-Cibatu<br class=""> 11
10 0.000147473 0.0014168
0.0000000 9900 0
0 0 0
1; %Saguling-Cirata<br class=""> 12
11 0.000195780 0.0021902
0.0000000 9900 0
0 0 0
1; %Bandung Selatan-Saguling<br class="">
13 12 0.013981960
0.1343320 0.1924070 9900
0 0 0
0 1; %Mandiracan-Bandung Selatan<br
class=""> 13 26
0.001012900 0.0097320
0.1393940 9900 0
0 0 0
1; %Mandiracan-Ujung Berung<br class=""> 26
12 0.000385400 0.0037034
0.0000000 9900 0
0 0 0
1; %Ujung Berung-Bandung Selatan<br class="">
30 13 0.000200000
0.0022500 0.0000000 9900
0 0 0
0 1; %Pemalang-Mandiracan<br
class=""> 14 30
0.000320000 0.0035300
0.0000000 9900 0
0 0 0
1; %Ungaran-Pemalang<br class=""> 15
14 0.000067660 0.0075703
0.4366440 9900 0
0 0 0
1; %Tanjung Djati-Ungaran<br class=""> 16
14 0.002979200 0.0286229
0.4099740 9900 0
0 0 0
1; %Surabaya Barat-Ungaran<br class=""> 14
20 0.000903612 0.0086815
0.0000000 9900 0
0 0 0
1; %Ungaran-Pedan<br class=""> 16
17 0.000140000 0.0013399
0.0000000 9900 0
0 0 0
1; %Surabaya Barat-Gresik<br class=""> 23
16 0.000398638 0.0044597
0.0000000 9900 0
0 0 0
1; %Grati-Surabaya Barat<br class=""> 18
5 0.000091200 0.0008765
0.0000000 9900 0
0 0 0
1; %Depok-Cibinong<br class=""> 19
18 0.001403000 0.0156967
0.9015550 9900 0
0 0 0
1; %Tasikmalaya-Depok<br class=""> 29
19 0.000544600 0.0060933
0.3499780 9900 0
0 0 0
1; %Kesugihan-Tasikmalaya<br class=""> 27
29 0.000110000 0.0012668
0.0000000 9900 0
0 0 0
1; %Adipala-Kesugihan<br class=""> 28
27 0.000038100 0.0004387
0.0000000 9900 0
0 0 0
1; %Cilacap-Adipala<br class=""> 20
29 0.000983900 0.0110073
0.6348830 9900 0
0 0 0
1; %Pedan-Kesugihan<br class=""> 21
20 0.002053100 0.0229694
0.3298180 9900 0
0 0 0
1; %Kediri-Pedan<br class=""> 22
21 0.001026500 0.0114847
0.7239300 9900 0
0 0 0
1; %Paiton-Kediri<br class=""> 22
23 0.000044100 0.0049343
0.3110350 9900 0
0 0 0
1; %Paiton-Grati<br class=""> 24
4 0.000297922 0.0028623
0.0000000 9900 0
0 0 0
1; %Balaraja-Gandul<br class=""> 25
14 0.002347961 0.0225581
0.3244670 9900 0
0 0 0
1; %Ngimbang-Ungaran<br class=""> 16
25 0.000596665 0.0057324
0.0000000 9900 0
0 0 0
1; %Surabaya Barat-Ngimbang<br class=""> ];<br class=""><br
class=""></div><div class=""><br class=""></div></div><div dir="ltr"
data-setdir="false" class="">Regards</div><div dir="ltr" data-setdir="false"
class="">Sabhan K<br class=""></div>
</div><div id="yahoo_quoted_1773503525" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial,
sans-serif;font-size:13px;color:#26282a;" class="">
<div class="">
Pada Senin, 21 Oktober 2019 20.55.34 WIB, Ray Zimmerman
<<a href="mailto:[email protected]" class="">[email protected]</a>> menulis:
</div>
<div class=""><br class=""></div>
<div class=""><br class=""></div>
<div class=""><div id="yiv6822524658" class=""><div class="">It
is possible that you have a disconnected network with an island that does not
have a reference bus, or maybe you have a zero impedance branch. I would begin
by looking at the output of <font class="yiv6822524658"
face="Courier">case_info()</font> for your case to see if you see anything
unusual.<div class="yiv6822524658"><br class="yiv6822524658"
clear="none"></div><div class="yiv6822524658"> Ray</div><div
class="yiv6822524658yqt3153578238" id="yiv6822524658yqt88331"><div
class="yiv6822524658"><br class="yiv6822524658" clear="none"><div class=""><br
class="yiv6822524658" clear="none"><blockquote class="yiv6822524658"
type="cite"><div class="yiv6822524658">On Oct 20, 2019, at 11:59 PM, Sabhan
Kanata <<a rel="nofollow" shape="rect" class="yiv6822524658"
ymailto="mailto:[email protected]" target="_blank"
href="mailto:[email protected]">[email protected]</a>>
wrote:</div><br class="yiv6822524658Apple-interchange-newline"
clear="none"><div class="yiv6822524658"><div class="yiv6822524658"><div
class="yiv6822524658yahoo-style-wrap" style="font-family:Helvetica Neue,
Helvetica, Arial, sans-serif;font-size:13px;"><div class="yiv6822524658"
dir="ltr"><div class="yiv6822524658"><div class="yiv6822524658" dir="ltr"><div
class="yiv6822524658"><div class="yiv6822524658" dir="ltr">Dear Prof
Ray,</div><div class="yiv6822524658" dir="ltr"><br class="yiv6822524658"
clear="none"></div><div class="yiv6822524658">I get a warning when running the
power flow as below:</div><br class="yiv6822524658" clear="none">What do you
think is the solution?</div><div class="yiv6822524658"><br
class="yiv6822524658" clear="none"></div></div><div class="yiv6822524658"><br
class="yiv6822524658" clear="none"></div><div class="yiv6822524658">Error
using \ </div>Matrix is singular to working precision.<br
class="yiv6822524658" clear="none"><br class="yiv6822524658" clear="none">Error
in newtonpf (line 100)<br class="yiv6822524658" clear="none"> dx =
-(J \ F);<br class="yiv6822524658" clear="none"><br class="yiv6822524658"
clear="none">Error in runpf (line 148)<br class="yiv6822524658"
clear="none">
[V, success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq,
mpopt);</div><div class="yiv6822524658"><br class="yiv6822524658"
clear="none"></div><div class="yiv6822524658" dir="ltr">Regards</div><div
class="yiv6822524658" dir="ltr">Sabhan K<br class="yiv6822524658"
clear="none"></div></div></div></div></div></blockquote></div><br
class="yiv6822524658" clear="none"></div></div></div></div></div>
</div>
</div></div></div></blockquote></div><br class=""></div></body></html>