Форум режимщиков

Set tvetv=Rastr.Tables("vetv")

Set tnode=Rastr.Tables("node")

Set ip=tvetv.Cols("ip")

Set iq=tvetv.Cols("iq")

Set bsh=tnode.Cols("bsh")

Set gsh=tnode.Cols("gsh")

Set b=tvetv.Cols("b")

Set g=tvetv.Cols("g")

Set tip=tvetv.Cols("tip")



'Обнуляем шунты

GrCor "node","bsh","","0"

GrCor "node","gsh","","0"



tvetv.SetSel("tip=1") 'трансформаторные ветви

k=tvetv.FindNextSel(-1)

While k <>-1

i=ip.Z(k)

j=iq.Z(k)

bt=abs(b.Z(k))

gt=abs(g.Z(k))



rastr.printp "трансформаторная ветвь " & i & "-" & j

'ищем индексы узлов начала и конца ветви

tnode.SetSel("ny=" & i) 

uzel1=tnode.FindNextSel(-1)

tnode.SetSel("ny=" & j)

uzel2=tnode.FindNextSel(-1)



'задаем половинные шунты

bsh.Z(uzel1) = bsh.Z(uzel1)+-bt/2

bsh.Z(uzel2) = bsh.Z(uzel2)+-bt/2

gsh.Z(uzel1) = gsh.Z(uzel1)+gt/2

gsh.Z(uzel2) = gsh.Z(uzel2)+gt/2



'обнуляем проводимости ветви

b.Z(k)=0

g.Z(k)=0

k=tvetv.FindNextSel(k)

Wend

Rastr.rgm("p")