【文章內(nèi)容簡介】 NO=INT(TM/T)+1LET M0=*H0*Q0/(N1*P0)DIM Q(0 TO 550,0 TO 10) ,H(0 TO 550,0 TO 10) ,X(0 TO 10),Z(0 TO 6),F(55),Y(55)DIM R1(0 TO 550),V1(0 TO 550),B(0 TO 550),TI(0 TO 550),TT(0 TO 550,0 TO 10)DIM U(0 TO 150),V2(0 TO 550)! 輸 入 原 始 數(shù) 據(jù) H/(R^2+V^2) , M/(R^2+V^2)FOR M=1 TO J READ F(M)NEXT MFOR M=1 TO J READ Y(M)NEXT MFOR K=0 TO BO READ U(K)NEXT K! 穩(wěn) 態(tài) 工 況 下 的 流 量,揚(yáng) 程 計 算FOR I=0 TO N LET H(0,0)=H0 LET H(0,I)=H(0,0)I*(H(0,0)H1)/N LET Q(0,I)=Q0*Z1 LET Z(I)=H(0,I) LET X(I)=Q(0,I)NEXT I! 賦 迭 代 初 值LET TI(0)=0LET R1(0)=1LET V1(0)=1LET V2(0)=1LET A=0LET D1=0LET D2=0LET B(0)=1! 暫 態(tài) 工 況 下 的 流 量,揚(yáng) 程 計 算FOR K=1 TO NO LET S=0 ! 中 間 斷 面 的 流 量,揚(yáng) 程 計 算 FOR I=1 TO N1 LET C3=G*A0/C LET C4=(F0*T)/(2*D0*A0) LET C5=C3*(H0h3) LET C6=(15*T*M0*G)/(PI*W2*N1) LET C1=X(I1)+C3*Z(I1)C4*X(I1)*ABS(X(I1)) LET C2=X(I+1)C3*Z(I+1)C4*X(I+1)*ABS(X(I+1)) LET Q(K,I)=*(C1+C2) LET H(K,I)=(C1Q(K,I))/C3 LET C5=C3*(H0h3) LET C6=(15*T*M0*G)/(PI*W2*N1) NEXT I ! 出 口 斷 面 的 流 量,揚(yáng) 程 計 算 LET H(K,N)=H1 LET C7=X(N1)+C3*Z(N1)C4*X(N1)*ABS(X(N1)) LET Q(K,N)=C7C3*H(K,N) ! 水 泵 斷 面 的 流 量,揚(yáng) 程 計 算 CALL SB ! 為 下 一 計 算 時 段 儲 存 數(shù) 據(jù) FOR I=0 TO N LET X(I)=Q(K,I) LET Z(I)=H(K,I) NEXT INEXT K! 繪Q－T，H－T，M－T 曲線CALL HTPRINTPRINT DT=。T。NO=。NO! 繪制表格CALL DGFOR K=0 TO NO STEP 1 FOR I=0 TO N ! 把對應(yīng)數(shù)據(jù)填入表格中 CALL SC NEXT INEXT K! 水 泵 斷 面 的 流 量,揚(yáng) 程 計 算 子 程 序!SUBSUB SB IF ZZ=Z1 THEN LET TI(K)=TI(K1)+T LET C8=X(1)C3*Z(1)C4*X(1)*ABS(X(1)) IF TI(K)(T3T) THEN LET Q(K,0)=0 LET H(K,0)=C8/C3 IF R1(K1)0 THEN LET R1(K)=(1SQR(1+4**C6*(R1(K1)+C6*B(K1))))/(2*()*C6) LET B(K)=*R1(K)^2 ELSE LET R1(K)=(1SQR(14**C6*(R1(K1)+C6*B(K1))))/(2**C6) LET B(K)=*R1(K)^2 END IF ELSE LET R1(K)=R1(K1)+A LET V1(K)=V1(K1)+D1 DO IF R1(K)=0 THEN LET E=180 LET E1=PI END IF IF V1(K)=0 THEN LET E=90 LET E1=PI/2 END IF LET E1=ATN(R1(K)/V1(K)) LET E=E1*180/PI IF E0 THEN IF R1(K)0 OR V1(K)0 THEN LET E=E LET E1=E1 END IF IF R1(K)0 OR V1(K)0 THEN LET E=180+E LET E1=PI+E1 END IF ELSE LET E=180+E LET E1=PI+E1 END IF LET M=INT(E/E0+1) LET A1=F(M)*MF(M+1)*(M1) LET A2=(F(M+1)F(M))/(E0*) LET A3=Y(M)*MY(M+1)*(M1) LET A4=(Y(M+1)Y(M))/(E0*) LET VT=V1(K)*V1(K)+R1(K)*R1(K) CALL SUB2 LET F8=C8+C3*H3C3*CV*Q0*Q0*V1(K)*ABS(V1(K)) LET F1=C5*A1*VT+C5*A2*VT*E1Z1*Q0*V1(K)+F8 LET F9=R1(K1)+C6*B(K1) LET F2=R1(K)C6*A3*VTC6*A4*VT*E1F9 LET F3=C5*(2*A1*R1(K)+A2*V1(K)+2*A2*R1(K)*E1) LET F4=C5*(2*A1*V1(K)A2*R1(K)+2*A2*V1(K)*E1)Z1*Q02*C3*cv*Q0^2*ABS(V1(K)) LET F5=1C6*(2*A3*R1(K)+A4*V1(K)+2*A4*R1(K)*E1) LET F6=C6*(2*A3*V1(K)+A4*R1(K)2*A4*V1(K)*E1) LET F7=F3*F6F4*F5 LET A=(F2*F4F1*F6)/F7 LET D1=(F1*F5F2*F3)/F7 LET R1(K)=R1(K)+A LET V1(K)=V1(K)+D1 LET S=S+1 IF S40 THEN EXIT DO LOOP UNTIL ABS(A) AND ABS(D1) LET A=R1(K)R1(K1) LET D1=V1(K)V1(K1) LET Q(K,0)=Q0*V1(K)*Z1 LET H(K,0)=(Q(K,0)C8)/C3 LET B(K)=(A3+A4*E1)*VT END IF END IF IF ZZZ1 AND ZZ0 THEN LET TI(K)=TI(K1)+T LET C10=X(1)C3*Z(1)C4*X(1)*ABS(X(1)) IF TI(K)(T3T) THEN LET Q(K,0)=Q0*(Z1ZZ) LET H(K,0)=(Q(K,0)C10)/C3 IF R1(K1)0 THEN LET R1(K)=(1SQR(1+4**C6*(R1(K1)+C6*B(K1))))/(2*()*C6) LET B(K)=*R1(K)*R1(K) ELSE LET R1(K)=(1SQR(14**C6*(R1(K1)+C6*B(K1))))/(2**C6) LET B(K)=*R1(K)*R1(K) END IF ELSE LET R1(K)=R1(K1)+A LET V1(K)=V1(K1)+D1 LET V2(K)=V2(K1)+D2 DO IF R1(K)=0 THEN LET E=180 LET E1=PI END IF IF V1(k)=0 THEN LET E=90 LET E1=PI/2 END IF LET E1=ATN(R1(K)/V1(K)) LET E=E1*180/PI IF E0 THEN IF R1(K)0 OR V1(K)0 THEN LET E=E LET E1=E1 END IF IF R1(K)0 OR V1(K)0 THEN LET E=180+E LET E1=PI+E1 END IF ELSE LET E=180+E LET E1=PI+E1 END IF LET M=INT(E/E0+1) LET A1=F(M)*MF(M+1)*(M1) LET A2=(F(M+1)F(M))/(E0*) LET A3=Y(M)*MY(M+1)*(M1) LET A4=(Y(M+1)Y(M))/(E0*) LET VT=V1(K)^2+R1(K)^2 CALL SUB2 LET R2=1 IF V2(K)=0 THEN LET EE=90 LET E2=PI/2 END IF LET E2=ATN(1/V2(K)) LET EE=E2*180/PI IF EE0 THEN LET EE=EE LET E2=E2 ELSE LET EE=180+EE LET E2=PI+E2 END IF