CPF power flow did not converge

Thu, 11 Jul 2019 10:54:18 -0700 

Hi all,

I am trying to analyse the influence of topology changes (like disconnect 
branches) on system voltage stability using the following code:

define_constants;



mpc =loadcase('FourAreaNetwork');   "
OutmailID: 123739424, List: 'matpower-l', MemberID: 82861091
SCRIPT: "load basecase
 
 mac.branch(:,11)=[0,0,1,0,1,1,1,1,1,1];"
TCL MERGE ERROR ( 07/11/2019 13:54:03 ): "invalid command name 
"0,0,1,0,1,1,1,1,1,1"change branch status

mpopt = mpoption('out.all',0,'verbose',2,'out.bus',1);
mpopt=mpoption(mpopt,'cpf.stop_at','NOSE','cpf.step',0.05);
mpct = mpc;          "
OutmailID: 123739424, List: 'matpower-l', MemberID: 82861091
SCRIPT: "set up target case with basecase
 mpct.gen(:,[PG QG]) = mpc.gen(:,[PG QG])*2
 mpct.bus(:,[PD QD]) = mpc.bus(:,[PD QD])*2
 results = runcpf(mpc, mpct,mpopt);
 
 After running this, I got the warning message ‘ CPF TERMINATION: Base case 
power flow did not converge.’
 Then I decreased the size of cpf step, but still got the same warning.
 
 Could you please suggest me how should I fix it? And I have attached the 
FourAreaNetwork case file.
 
 Thanks in advance!
 
 Best regards,
 Cynthia
 
 
  --_000_CY4PR08MB2808DEF11DE25F2332D96BACDFF30CY4PR08MB2808namp_ Content-Type: 
text/html; charset="utf-8" Content-ID: 
<[email protected]> 
Content-Transfer-Encoding: base64  <html>
 <head>
 <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
 </head>
 <body>
 <div class="" style="word-wrap:break-word; line-break:after-white-space">Hi 
all,
 <div class=""><br class="">
 <div class=""><font face="Times" size="4" class="">I am trying to analyse the 
influence of topology changes (like disconnect branches) on system voltage 
stability using the following c</font><font face="Times" size="4" 
class="">ode:</font></div>
 </div>
 <div class=""><font face="Times" size="4" class=""><br class="">
 </font></div>
 <div class="">
 <div class="" style="margin:0px; font-stretch:normal; 
line-height:normal"><font face="Times" size="4" 
class="">define_constants;</font></div>
 <p class="" style="margin:0px; font-stretch:normal; line-height:normal; 
min-height:12px">
 <font face="Times" size="4" class="">&nbsp;<br 
class="x_webkit-block-placeholder">
 </font></p>
 <div class="" style="margin:0px; font-stretch:normal; 
line-height:normal"><font face="Times" size="4" class="">mpc 
=loadcase('FourAreaNetwork'); &nbsp;"
TCL MERGE ERROR ( 07/11/2019 13:54:03 ): "wrong # args: should be "set varName 
?newValue?" load basecase</font></div>
</div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
face="Times" size="4" class=""><br class="">
</font></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
face="Times" size="4" class=""><span class="" 
style="">mac.branch(:,11)=[0,0,1,0,1,1,1,1,1,1]; %%change branch 
status</span></font></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
face="Times" size="4" class=""><span class="" style=""><br class="">
</span></font></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal">
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style=""><font face="Times" size="4" class="">mpopt = 
mpoption('out.all',0,'verbose',2,'out.bus',1);&nbsp;</font></span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style=""><font size="4" face="Times" 
class="">mpopt=mpoption(mpopt,'cpf.stop_at','NOSE','cpf.step',0.05);&nbsp;</font></span></div>
</div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
face="Times" class="">
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style=""><font size="4" class="">mpct = mpc; &nbsp; &nbsp; &nbsp; 
&nbsp; &nbsp;%% set up target case with basecase</font></span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style=""><font size="4" class="">mpct.gen(:,[PG QG]) = mpc.gen(:,[PG 
QG])*2</font></span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style=""><font size="4" class="">mpct.bus(:,[PD QD]) = mpc.bus(:,[PD 
QD])*2</font></span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style=""><font size="4" class="">results = runcpf(mpc, 
mpct,mpopt);</font></span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style=""><font size="4" class=""><br class="">
</font></span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
size="4" class=""><span class="" style="">After running this, I got the warning 
message ‘&nbsp;</span></font><font size="4" class=""><span class="" 
style="">CPF TERMINATION: Base case
 power flow did not converge.’</span></font></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
size="4" class=""><span class="" style="">Then I&nbsp;decreased the size of cpf 
step, but still got the same warning.</span></font></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
size="4" class=""><span class="" style=""><br class="">
</span></font></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
size="4" class=""><span class="" style="">Could you please suggest me 
how&nbsp;should I fix it?&nbsp;</span></font><span class="" 
style="font-size:large">And I have attached the FourAreaNetwork
 case file.</span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style="font-size:large"><br class="">
</span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style="font-size:large">Thanks in advance!</span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style="font-size:large"><br class="">
</span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style="font-size:large">Best regards,</span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><span 
class="" style="font-size:large">Cynthia</span></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"></div>
</font></div>
</div>
<div class="" style="word-wrap:break-word; line-break:after-white-space">
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
face="Times" class="">
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"></div>
<div class="" style="margin:0px; font-stretch:normal; line-height:normal"><font 
size="4" class=""><span class="" style=""><br class="">
</span></font></div>
<div class="" style=""><br class="">
</div>
</font></div>
</div>
</body>
</html>

--_000_CY4PR08MB2808DEF11DE25F2332D96BACDFF30CY4PR08MB2808namp_--

--_004_CY4PR08MB2808DEF11DE25F2332D96BACDFF30CY4PR08MB2808namp_
Content-Type: application/octet-stream; name="FourAreaNetwork.m"
Content-Description: FourAreaNetwork.m
Content-Disposition: attachment; filename=FourAreaNetwork.m; size=2432;
        creation-date="Thu, 11 Jul 2019 17:53:49 GMT";
        modification-date="Thu, 11 Jul 2019 17:53:49 GMT"
Content-ID: <[email protected]>
Content-Transfer-Encoding: base64

ZnVuY3Rpb24gbXBjID0gRm91ckFyZWFOZXR3b3JrCQkJCQkJCQkJCQkJCQkJCQkJCQkJDQptcGMu
dmVyc2lvbiA9ICcyJzsJCQkJCQkJCQkJCQkJCQkJCQkJCQkNCm1wYy5iYXNlTVZBID0gMTAwOwkJ
CQkJCQkJCQkJCQkJCQkJCQkJCQ0KDQolJSBCdXMgRGF0YQkJCQkJCQkJCQkJCQkJCQkJCQkJCQ0K
JQlidXNfaQl0eXBlCVBkCVFkCUdzCUJzCWFyZWEJVm0JVmEJYmFzZUtWCXpvbmUJVm1heAlWbWlu
CQkJCQkJCQkNCm1wYy5idXMgPSBbCQkJCQkJCQkJCQkJCQkJCQkJCQkJDQogICAgMSAgIAkyICAg
CTAJMAkwCTAJMQkxCTAJMjIJMQkxLjA1CTAuOTU7DQogICAgMiAgIAkxICAgCTI4MDAJMAkwCTAJ
MQkxCTAJNDAwCTEJMS4wNQkwLjk1Ow0KICAgIDMgICAJMSAgIAkxNzAwCTAJMAkwCTEJMQkwCTQw
MAkxCTEuMDUJMC45NTsNCiAgICA0ICAgCTIgICAJMAkwCTAJMAkxCTEJMAkyMgkxCTEuMDUJMC45
NTsNCiAgICA1ICAgCTIgICAJMAkwCTAJMAkxCTEJMAkxMy44CTEJMS4wNQkwLjk1Ow0KICAgIDYg
ICAJMSAgIAkxMzAwCTAJMAkwCTEJMQkwCTQwMAkxCTEuMDUJMC45NTsNCiAgICA3ICAgCTEgICAJ
MTcwMAkwCTAJMAkxCTEJMAk0MDAJMQkxLjA1CTAuOTU7DQogICAgOCAgIAkzICAgCTAJMAkwCTAJ
MQkxCTAJMTMuOAkxCTEuMDUJMC45NTsJCQkJCQkJCQ0KXTsJCQkJCQkJCQkJCQkJCQkJCQkJCQkN
CiUlIEdlbmVyYXRvciBEYXRhCQkJCQkJCQkJCQkJCQkJCQkJCQkJDQolCWJ1cwlQZwlRZwlRbWF4
CVFtaW4JVmcJbUJhc2UJc3RhdHVzCVBtYXgJUG1pbglQYzEJUGMyCVFjMW1pbglRYzFtYXgJUWMy
bWluCVFjMm1heAlyYW1wX2FnYwlyYW1wXzEwCXJhbXBfMzAJcmFtcF9xCWFwZg0KbXBjLmdlbiA9
IFsJCQkJCQkJCQkJCQkJCQkJCQkJCQkNCiAgICAxCTE3NDAJMAkxMjYzCS0xMjYzCTEJMTAwCTEJ
MjYwOC4yCTAJMAkwCTAJMAkwCTAJMAkwCTAJMAkwOw0KICAgIDQJMTgwMAkwCTEyNjMJLTEyNjMJ
MQkxMDAJMQkyNjA4LjIJMAkwCTAJMAkwCTAJMAkwCTAJMAkwCTA7DQogICAgNQkxOTYwCTAJMTE0
Mi4wMwktMTE0Mi4wMwkxCTEwMAkxCTIzNTgJMAkwCTAJMAkwCTAJMAkwCTAJMAkwCTA7DQogICAg
OAkyMDE5CTAJMTI1Ni4yMwktMTI1Ni4yMwkxCTEwMAkxCTI1OTMuOAkwCTAJMAkwCTAJMAkwCTAJ
MAkwCTAJMDsNCl07CQkJCQkJCQkJCQkJCQkJCQkJCQkJDQolJSBUcmFuc21pc3Npb24gTGluZXMg
RGF0YQkJCQkJCQkJCQkJCQkJCQkJCQkJCQ0KJQlmYnVzCXRidXMJcgl4CWIJcmF0ZUEJcmF0ZUIJ
cmF0ZUMJcmF0aW8JYW5nbGUJc3RhdHVzCWFuZ21pbglhbmdtYXgJCQkJCQkJCQ0KbXBjLmJyYW5j
aCA9IFsJCQkJCQkJCQkJCQkJCQkJCQkJCQkNCiANCiAgICAyCTMJMC4wMDE4NzUJMC4wMjY1NjI1
CTAJNDAwCTQwMAk0MDAJMAkwCTEJLTM2MAkzNjAuMDA7DQogICAgMgk2CTAuMDA3NQkwLjEwNjI1
CTAJNDAwCTQwMAk0MDAJMAkwCTEJLTM2MAkzNjAuMDA7DQogICAgMgk3CTAuMDA3NQkwLjEwNjI1
CTAJNDAwCTQwMAk0MDAJMAkwCTEJLTM2MAkzNjAuMDA7DQogICAgDQogICAgMwk2CTAuMDAxODc1
CTAuMDI2NTYyNQkwCTQwMAk0MDAJNDAwCTAJMAkxCS0zNjAJMzYwLjAwOw0KICAgIDMJNwkwLjAw
MTg3NQkwLjAyNjU2MjUJMAk0MDAJNDAwCTQwMAkwCTAJMQktMzYwCTM2MC4wMDsNCiAgIA0KICAg
IDYJNwkwLjAwNzUJMC4xMDYyNQkwCTQwMAk0MDAJNDAwCTAJMAkxCS0zNjAJMzYwLjAwOw0KICAg
DQogICAgDQogICAgMQkyCTAJMC4wMDEJMAkzMDAwCTMwMDAJMzAwMAkwCTAJMQktMzYwCTM2MC4w
MDsNCiAgICAzCTQJMAkwLjAwMQkwCTMwMDAJMzAwMAkzMDAwCTAJMAkxCS0zNjAJMzYwLjAwOw0K
ICAgIDUJNgkwCTAuMDAxCTAJMzAwMAkzMDAwCTMwMDAJMAkwCTEJLTM2MAkzNjAuMDA7DQogICAg
Nwk4CTAJMC4wMDEJMAkzMDAwCTMwMDAJMzAwMAkwCTAJMQktMzYwCTM2MC4wMDsNCl07CQkJCQkJ
CQkJCQkJCQkJCQkJCQkJDQptcGMuYXJlYXMgPVsJCQkJCQkJCQkJCQkJCQkJCQkJCQkNCjEgMTsJ
CQkJCQkJCQkJCQkJCQkJCQkJCQkNCl07CQkJCQkJCQkJCQkJCQkJCQkJCQkJDQolJSBHZW5lcmF0
b3IgQ29zdCBEYXRhIChNYWRlIHVwKQkJCQkJCQkJCQkJCQkJCQkJCQkJCQ0KJQkxCXN0YXJ0dXAJ
c2h1dGRvd24Jbgl4MQl5MQkuLi4JeG4JeW4JCQkJCQkJCQkJCQkNCiUJMglzdGFydHVwCXNodXRk
b3duCW4JYyhuLTEpCS4uLgljMAkJCQkJCQkJCQkJCQkJDQptcGMuZ2VuY29zdCA9IFsJCQkJCQkJ
CQkJCQkJCQkJCQkJCQkNCjIJMTUwMAkwCTMJMAkwLjEJMC4wNTsJCQkJCQkJCQkJCQkJCQkNCjIJ
MTUwMAkwCTMJMAkwLjEJMC4wMTsJCQkJCQkJCQkJCQkJCQkNCjIJMTUwMAkwCTMJMAkwLjEJMC41
OwkJCQkJCQkJCQkJCQkJCQ0KMgkxNTAwCTAJMwkwCTAuMQkwLjAyOwkJCQkJCQkJCQkJCQkJCQkJ
CQkJCQkJCQkJCQkNCl07CQkJCQkJCQkJCQkJCQkJCQkJCQkJDQo=

--

Reply via email to