When using AC modeling, the power flow is complex and it is not that same at 
the two ends of a branch, so you may want to use the average MVA as the flow.

define_constants;
r = runpf(…);
Af = sqrt(r.branch(:, PF)^2 + r.branch(:, QF)^2);
At = sqrt(r.branch(:, PT)^2 + r.branch(:, QT)^2);
flow = (Af + At)/2;
loading  = flow ./ r.branch(:, RATE_A);

The total active power losses can be computed by simply adding the injections 
at each end of the branch.

losses = r.branch(:, PF) + r.branch(:, PT);

-- 
Ray Zimmerman
Senior Research Associate
419A Warren Hall, Cornell University, Ithaca, NY 14853
phone: (607) 255-9645




On Oct 1, 2012, at 12:43 PM, Evangelos Galinas <vagelisgali...@msn.com> wrote:

> Hallo Basem,
> 
> I guess, except if it's already somewhere else calculated, you can just grab 
> the abs(mpc.branch(:,14)) (power injected at one end)and divide by the 
> mpc.branch(:,6) (rating A). I
> 
> Fot the losses just calculate the current and multiply it by the losses  and 
> save it at some vector
> 
> cheers
> Date: Mon, 1 Oct 2012 18:33:33 +0200
> From: basem.id...@iwes.fraunhofer.de
> To: matpowe...@cornell.edu
> Subject: Line loading
> 
> Hi all,
>  
> I am a beginner user of Matpower, and I would like to know how I can 
> calculate the loding of a branch (e.g. in percentage). I also need to know 
> how I can save the losses of the branches in a matrix.
>  
> Thanks in advance for any answer,
> Basem Idlbi
>  
>  
>  
> Basem Idlbi
> Division Systems Engineering and Distribution Grids
> Fraunhofer Institute for Wind Energy and Energy System Technology, IWES
> Königstor 59 | 34119 Kassel / Germany
> Phone +49 561 7294-166 | Fax +49 561 7294-200
>  

Reply via email to