---------- Forwarded message ---------
From: akhil <[email protected]>
Date: Thu, 13 Jun 2019 at 19:52
Subject: Maximum load limit for radial distribution cases.
To: <[email protected]>
I am running the radial distribution power flow which is based on backward
forward sweep algorithm. for the network data which is provided with
matpower the radial load flow converges successfully. But to get the
maximum load limit of the system, when I am running the load flow with
increased load at each bus in a step of 1 % the radial load flow does not
converge after some point.
The change in Load is as per Pnew=Lambda*P0 and Qnew=Lambda*Q0 with 1 %
change in load at all buses simultaneously.
I am Trying to get the minimum bus voltage of the system after each change
in load and running corresponding radial flow successively until the
minimum voltage of the system (bus with minimum voltage) getting reduced
smoothly.
The code which I run is attached herewith for your reference.
Thank you
Yours
Akhilesh K Barnwal
Ph.D. Research Scholar
IIT BHU Varanasi.
mail: [email protected]
--
@khil
clc
define_constants
mpc0 = loadcase(case33bw);
mpc0.bus(:,VMIN)=0.2
status = find(mpc0.branch(:,BR_STATUS)==0);
mpc0.branch(status,:)=[];
Lambda=0;i=0;
Vnew=1; Vold=5;
while Vnew<=Vold % i<=15
i=i+1
mpopt=mpoption('model', 'ac', 'pf.alg', 'isum', 'pf.radial.max_it', 1500);
% increase system load at each bus P or/and Q in small steps of 1% (1+Lambda*P0)
mpc.bus(:,PD)=Lambda*mpc0.bus(:,PD);
mpc.bus(:,QD)=Lambda*mpc0.bus(:,QD);
results=runpf(mpc,mpopt);
Vmin=min(results.bus(:,VM));
Vold=Vnew;
Vnew=Vmin;
Lambda=Lambda+0.1;
end
Lambda
Vold
Vnew
