Re: Effect of loads in IEEE 39 bus system
Hi Shri, Sorry for not stating my problem clearly. I need to implement ZIP load model in matpower to evaluate it's effect, instead of constant power load in runcpf. P = P0[ a1 (v/v0)2 + a2(v/v0) + a3 ] Vo is the rated voltage of the device, and Po and Qo is the power consumed at rated voltage. This is the equation I need to implement. I tried understanding the runcpf code, But I couldn't make out where to make those changes. Is it possible to tell where to make those changes or the name of the matlab file used for this? Thanking you and With Regards, Arun Nair On Mon, Feb 23, 2015 at 7:46 AM, Abhyankar, Shrirang G. abhy...@mcs.anl.gov wrote: Arun, Please help us to understand the meaning of stuck by providing relevant details of (a) equations you are trying to implement (b) what have you attempted till now c) code snippets (d) errors/warnings. Shri From: arun s nair arunsnair...@gmail.com Reply-To: MATPOWER discussion forum matpowe...@list.cornell.edu Date: Sat, 21 Feb 2015 18:10:11 -0700 To: MATPOWER discussion forum matpowe...@list.cornell.edu Subject: Effect of loads in IEEE 39 bus system Hi All, I am trying to implement the different load models((like ZIP load model, composite load model, exponential load model etc ) in matpower to analyse it's effect on the load buses in an IEEE 39 bus system (using cpf ). But now I am stuck as to how to implement it in matpower. Can anyone please give me some guidelines as to how I should proceed so to get it implemented. Your help will be greatly appreciated, its for my graduate research. Thanking you and With Regards, Arun Nair
Re: runcpf problem
Thank you, it now works. I have some questions: * How to run cpf for other load models (composite, constant current and constant impedance)? * How to calculate and plot stability index? Thank you in advance Amal On Monday, February 23, 2015 3:56 PM, Abhyankar, Shrirang G. abhy...@mcs.anl.gov wrote: Please check if your base case and target case have the same number of buses. Shri From: Amal Elramly amal_...@yahoo.com Reply-To: MATPOWER discussion forum matpowe...@list.cornell.edu Date: Mon, 23 Feb 2015 12:58:12 + To: MATPOWER Discussion Forum matpowe...@list.cornell.edu Subject: runcpf problem Dear, runcpf function does not work for 33 radial distribution system Error Message below: MATPOWER Version 5.0, 17-Dec-2014 -- AC Continuation Power FlowError using - Matrix dimensions must agree.Error in runcpf (line 213)Sxfr = Sbust - Sbusb; could you help me Best Regards, Amal Hassan
Re: Convergence Issue for Distribution Network Power Flow Analysis
rem On Tue, Feb 24, 2015 at 2:03 AM, #ZHANG TIAN# tzhang...@e.ntu.edu.sg wrote: Hello! Yes, by changing the LV impedances I can get the correct results. This is worth noting that the data in the “Model Explorer” in PowerWorld Simulator is in per unit based on the whole system, rather than its real value. Many thanks. Tian *From:* bounce-118850794-72551...@list.cornell.edu [mailto: bounce-118850794-72551...@list.cornell.edu] *On Behalf Of *Carlos E Murillo-Sanchez *Sent:* Monday, February 23, 2015 22:56 *To:* MATPOWER discussion forum *Subject:* Re: Convergence Issue for Distribution Network Power Flow Analysis Hi; You have used the TAP column of the branch table to indicate two transformers with a 55:1 ratio, and then you used the LV impedances, I think, still in PU with a base voltage of 22KV; they are pretty large. The TAP column should be used to specify small variations about 1.0 for true tap-changing transformers. All branch impedances should be specified using a single pu system. I believe that the following is what you are actually trying to run: define_constants; mpc = loadcase('caseSG'); mpc.branch(3:end, BR_R:BR_X) = mpc.branch(3:end, BR_R:BR_X) / 55^2; mpc.branch(3:4, TAP) = 0; runpf(mpc); #ZHANG TIAN# wrote: Dear All, I am using Matpower for a distribution network power flow analysis. The data I used are taken from a real network, and the network is working fine using PowerWorld simulator. Basically, it is a radial network with 14 buses. Bus 1 is the connection point to the main grid, which I modelled it as a slack bus (type 3). However, when converting the model into a Matpower case, I can only get a failure notification “Did NOT converge (0.01 seconds)”. Does anyone have any idea about what might be the problem? Or what are the options for me to debug? I attached my case data in this email for your reference. Thank you all very much. Best regards, Tian
Re: Convergence Issue for Distribution Network Power Flow Analysis
Please make sure that the loads i.e. PD and QD in mpc.bus are in p.u. Run the program again. It'll works. good luck. Andervazh On Tue, Feb 24, 2015 at 1:09 PM, Carlos Gonzalez Almeida cgonzalezalme...@gmail.com wrote: rem On Tue, Feb 24, 2015 at 2:03 AM, #ZHANG TIAN# tzhang...@e.ntu.edu.sg wrote: Hello! Yes, by changing the LV impedances I can get the correct results. This is worth noting that the data in the “Model Explorer” in PowerWorld Simulator is in per unit based on the whole system, rather than its real value. Many thanks. Tian *From:* bounce-118850794-72551...@list.cornell.edu [mailto: bounce-118850794-72551...@list.cornell.edu] *On Behalf Of *Carlos E Murillo-Sanchez *Sent:* Monday, February 23, 2015 22:56 *To:* MATPOWER discussion forum *Subject:* Re: Convergence Issue for Distribution Network Power Flow Analysis Hi; You have used the TAP column of the branch table to indicate two transformers with a 55:1 ratio, and then you used the LV impedances, I think, still in PU with a base voltage of 22KV; they are pretty large. The TAP column should be used to specify small variations about 1.0 for true tap-changing transformers. All branch impedances should be specified using a single pu system. I believe that the following is what you are actually trying to run: define_constants; mpc = loadcase('caseSG'); mpc.branch(3:end, BR_R:BR_X) = mpc.branch(3:end, BR_R:BR_X) / 55^2; mpc.branch(3:4, TAP) = 0; runpf(mpc); #ZHANG TIAN# wrote: Dear All, I am using Matpower for a distribution network power flow analysis. The data I used are taken from a real network, and the network is working fine using PowerWorld simulator. Basically, it is a radial network with 14 buses. Bus 1 is the connection point to the main grid, which I modelled it as a slack bus (type 3). However, when converting the model into a Matpower case, I can only get a failure notification “Did NOT converge (0.01 seconds)”. Does anyone have any idea about what might be the problem? Or what are the options for me to debug? I attached my case data in this email for your reference. Thank you all very much. Best regards, Tian
Re: Convergence Issue for Distribution Network Power Flow Analysis
No. The loads, PD and QD, are NOT in p.u., they are in MW and MVAr. Check the documentation for caseformat or Table B-1 in the User’s Manual http://www.pserc.cornell.edu/matpower/manual.pdf. Ray On Feb 24, 2015, at 7:32 AM, mohammad reza Andervazh mr.anderv...@gmail.com wrote: Please make sure that the loads i.e. PD and QD in mpc.bus are in p.u. Run the program again. It'll works. good luck. Andervazh On Tue, Feb 24, 2015 at 1:09 PM, Carlos Gonzalez Almeida cgonzalezalme...@gmail.com mailto:cgonzalezalme...@gmail.com wrote: rem On Tue, Feb 24, 2015 at 2:03 AM, #ZHANG TIAN# tzhang...@e.ntu.edu.sg mailto:tzhang...@e.ntu.edu.sg wrote: Hello! Yes, by changing the LV impedances I can get the correct results. This is worth noting that the data in the “Model Explorer” in PowerWorld Simulator is in per unit based on the whole system, rather than its real value. Many thanks. Tian From: bounce-118850794-72551...@list.cornell.edu mailto:bounce-118850794-72551...@list.cornell.edu [mailto:bounce-118850794-72551...@list.cornell.edu mailto:bounce-118850794-72551...@list.cornell.edu] On Behalf Of Carlos E Murillo-Sanchez Sent: Monday, February 23, 2015 22:56 To: MATPOWER discussion forum Subject: Re: Convergence Issue for Distribution Network Power Flow Analysis Hi; You have used the TAP column of the branch table to indicate two transformers with a 55:1 ratio, and then you used the LV impedances, I think, still in PU with a base voltage of 22KV; they are pretty large. The TAP column should be used to specify small variations about 1.0 for true tap-changing transformers. All branch impedances should be specified using a single pu system. I believe that the following is what you are actually trying to run: define_constants; mpc = loadcase('caseSG'); mpc.branch(3:end, BR_R:BR_X) = mpc.branch(3:end, BR_R:BR_X) / 55^2; mpc.branch(3:4, TAP) = 0; runpf(mpc); #ZHANG TIAN# wrote: Dear All, I am using Matpower for a distribution network power flow analysis. The data I used are taken from a real network, and the network is working fine using PowerWorld simulator. Basically, it is a radial network with 14 buses. Bus 1 is the connection point to the main grid, which I modelled it as a slack bus (type 3). However, when converting the model into a Matpower case, I can only get a failure notification “Did NOT converge (0.01 seconds)”. Does anyone have any idea about what might be the problem? Or what are the options for me to debug? I attached my case data in this email for your reference. Thank you all very much. Best regards, Tian
Errors On Radial Distribution Network OPF
To Whom It May Concern, I built the network of RBTS and it is shown below. Now, the OPF of the network is running correctly. However, when I disconnected the line 70 which links Supply Point 1 (SP1) and Supply Point 2 (SP2), the SP2 and SP3 cannot get the power supply. The SP2 and SP3 suppose to get supply due to the line 71 which links SP2 and SP3 at that condition. Because the power can flow to SP3 via line 68 and 69 from SP1, then the SP2 can get power via line 71 at that condition. Do you know what happened to the network and how can I fix it ? I list the case lines at the end of content body. Best Regards, Hui [cid:a8f7b1a4-2893-4e1c-946d-68212f2fbd64] function mpc = caseRBTS %CASEcase33kv-11kvPower flow data of Distribution System for.RBTS BUS 4 %% MATPOWER Case Format : Version 4.1 mpc.version = '2.0'; %%- Power Flow Data -%% %% system MVA base mpc.baseMVA = 100; %% bus data % bus_i type Pd Qd Gs Bs Area Vm Va baseKV zone Vmax Vmin mpc.bus = [ 1 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 2 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 3 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 4 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 5 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 6 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 7 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 8 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 9 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94; 10 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 11 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94; 12 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94; 13 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 14 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94; 15 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 16 1 2.445 0.8036 0 0 1 1 360 11 1 1.06 0.94; 17 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 18 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94; 19 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 20 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 21 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 22 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 23 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 24 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 25 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94; 26 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 27 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94; 28 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 29 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94; 30 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94; 31 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 32 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 33 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 34 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 35 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 36 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 37 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 38 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94; 39 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 40 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94; 41 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 42 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94; 43 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94; 44 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 45 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94; 46 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 47 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94; 48 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 49 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94; 50 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 51 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94; 52 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 53 1 1.63 0.5358 0 0 1 1 360 11 1 1.06 0.94; 54 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 55 1 2.445 0.8036 0 0 1 1 360 11 1 1.06 0.94; 56 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 57 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 58 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 59 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 60 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 61 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 62 1 0.8869 0.2915 0 0 1 1 360 11 1 1.06 0.94; 63 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 64 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94; 65 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 66 1 0.8137 0.2675 0 0 1 1 360 11 1 1.06 0.94; 67 1 0.6714 0.2207 0 0 1 1 360 11 1 1.06 0.94; 68 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 69 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 70 1 0 0 0 0 1 1 360 11 1 1.06 0.94; 71 3 0 0 0 0 1 1 360 33 1 1.03 0.97; 72 1 0 0 0 0 1 1 360 33 1 1.06 0.94; 73 1 0 0 0 0 1 1 360 33 1 1.06 0.94; ]; %% generator data % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf mpc.gen = [ 2 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0 0 00 0 0 0 0 0; 4 0.8869 0.2915 0.2915 -0.2915 1.00 100 1 0.8869 0 0 0 0