Курсовая работа: Расчет системы водоснабжения
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);
PRINT ";flux="; using ("#.###", flux(i));
PRINT ";press="; using ("##.###", press(i))
next i
PRINT "dH1="; using ("##.#####", dH1) ;";dH2="; using ("##.#####",dH2); ";dH3="; using("##.#####", dH3)