Курсовая работа: Расчет системы водоснабжения

for i=1 to N

s(i)=0.001735*(1+a)*length(i)/diam(i)^5.3

next i

[A]

dH1=s(3)*flux(3)*Abs(flux(3))-s(7)*flux(7)*Abs(flux(7))-s(4)*flux(4)*Abs(flux(4))

dH2=s(4)*flux(4)*Abs(flux(4))-s(5)*flux(5)*Abs(flux(5))-s(8)*flux(8)*Abs(flux(8))

dH3=s(2)*flux(2)*Abs(flux(2))+s(8)*flux(8)*Abs(flux(8))-s(6)*flux(6)*Abs(flux(6))-s(9)*flux(9)*Abs(flux(9))

if Abs(dH1)<dHmin and Abs(dH2)<dHmin and Abs(dH3)<dHmin goto [B]

dq1=0.5*dH1/(s(4)*Abs(flux(4))+s(7)*Abs(flux(7))+s(3)*Abs(flux(3)))

dq2=0.5*dH2/(s(4)*Abs(flux(4))+s(5)*Abs(flux(5))+s(8)*Abs(flux(8)))

dq3=0.5*dH3/(s(8)*Abs(flux(8))+s(2)*Abs(flux(2))-s(9)*Abs(flux(9))+s(6)*Abs(flux(6)))

flux(1)=flux(1)

flux(2)=flux(2)-dq3

flux(3)=flux(3)-dq1

flux(4)=flux(4)+dq1-dq2

flux(5)=flux(5)+dq2

flux(6)=flux(6)+dq3

flux(7)=flux(7)+dq1

flux(8)=flux(8)+dq2-dq3

flux(9)=flux(9)-dq3

goto [A]

[B]

for i=1 to N

press(i)=s(i)*flux(i)*Abs(flux(i))

next i

for i=1 to N

PRINT "i="; i;

PRINT ";lenght="; length(i);

PRINT ";diam="; diam(i);