Hi, Thanks for the reply, I have decreased the load and now my cpf is converged.
Best regards, Cynthia 在 2019年7月11日,19:30,Abhyankar, Shrirang G <shrirang.abhyan...@pnnl.gov<mailto:shrirang.abhyan...@pnnl.gov>> 写道: Cynthia, Please see this MATPOWER FAQ<https://matpower.org/doc/faq/#pfconvergence> on debugging power flow divergence. Also, you forgot to attach your case file. Thanks, Shri From: <bounce-123739424-83373...@list.cornell.edu<mailto:bounce-123739424-83373...@list.cornell.edu>> on behalf of W Ruoyu <cynthiaw1...@hotmail.com<mailto:cynthiaw1...@hotmail.com>> Reply-To: MATPOWER discussion forum <matpowe...@list.cornell.edu<mailto:matpowe...@list.cornell.edu>> Date: Thursday, July 11, 2019 at 10:54 AM To: "MATPOWER-L@cornell.edu<mailto:MATPOWER-L@cornell.edu>" <MATPOWER-L@cornell.edu<mailto:MATPOWER-L@cornell.edu>> Subject: CPF power flow did not converge 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: 83373175 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: 83373175 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: <8f94c2354585e4468b909975c5f41...@namprd08.prod.outlook.com<mailto:8f94c2354585e4468b909975c5f41...@namprd08.prod.outlook.com>> 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=""> <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'); " 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]; " OutmailID: 123739522, List: 'matpower-l', MemberID: 84159315 SCRIPT: "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); </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); </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; " TCL MERGE ERROR ( 07/11/2019 14:30:28 ): "invalid command name "change" 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 ‘ </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 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 should I fix it? </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: <0847624d7891e14eaaca0c3b07854...@namprd08.prod.outlook.com<mailto:0847624d7891e14eaaca0c3b07854...@namprd08.prod.outlook.com>> 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= --
FourAreaNetwork.m
Description: FourAreaNetwork.m