【正文】 0││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00││ .00││ .00││┃ │ .00││ .00││ .00│。 NEXT I PRINT2 PRINT2: ┠───┼。 PRINT2, USING.│: TI(K)。 FOR K1=1 TO 9 PRINT2: ━━━┿。 PRINT2,USING │:0,0,0,1,1,2,2,3,3,4。 FOR K1=1 TO 9 PRINT2:━━━┿。 NEXT K1 PRINT2: M/S │。 FOR K1=1 TO 4 PRINT2: M/S │。 揚 程│ PRINT2:┃。 揚 程│。 NEXT K1 PRINT2: ━━━┓ PRINT2: ┃時 間│。 NEXT K PLOT LINES:TT(NO,I)+5,B(NO)*2 PRINT PAUSE 1END SUBDATA 309,4,1026,1470,DATA ,15,5,55,0DATA 8,86,18,! PUMP NS=77 DATADATA , ,.246,DATA ,DATA ,DATA , ,DATA ,DATA ,DATA , 0 ,DATA , , , ,DATA , , ,.586DATA ,0707 ,DATA ,data ,! 閥門流量系數(shù)DATA , ,DATA ,SUB DG PRINT2: ┏━━━┯。5,10 PLOT LINES:5,4。T$(I)amp。25,15 FOR K=1 TO NO IF K=1 THEN LET TT(K,I)=0 ELSE LET TT(K,I)=TT(K1,I)+T END IF PLOT LINES:TT(K,I)+5,H(K,I)/。HY$ PLOT TEXT ,AT 15,55:A$ PLOT TEXT ,AT 3,55:H(M) PLOT TEXT ,AT 23,18:T(S) PLOT TEXT ,AT 3,18:(0,0) PLOT LINES:5,20。 NEXT K PLOT LINES:TT(NO,I)+5,Q(NO,I)/20 PAUSE 1 CLEAR NEXT I SET MODE EGAHIRES SET WINDOW 0,26,25,60 FOR I=0 TO N LET HY$= ) 的 HT 圖 LET A$=X$amp。5,10 PLOT LINES:5,0。T$(I)amp。NO=。(21)蝶閥開度5度等分點數(shù):18(22)蝶閥全開設計流量時的水頭損失:(23)水泵比轉數(shù):60(由于水泵全特性曲線資料的短缺采用相近的數(shù)據(jù),比轉數(shù)為77)二). ! 用水錘特征線法求解停泵水錘程序 CLEARINPUT PROMPT s NAME? :s$OPEN2:NAME s$,CREATE NEWOLDSET 2: MARGIN 120ERASE 2INPUT PROMPT 并聯(lián)水泵臺數(shù)Z1=:Z1INPUT PROMPT 事故停泵臺數(shù)ZZ=:ZZINPUT PROMPT 蝶閥二階段關閉時間T2=:T2! 程序中所用字母的有關說明! L0___管路總長度(m) D0___管道直徑(m)! A0___管道斷面積(m^2) N ___管路分段數(shù)! C ___水擊波速(m/s) F0___管道摩阻系數(shù)! H0___水泵額定揚程(m) Q0___水泵額定流量(m^3/s)! N0___水泵額定轉速(rpm) P0___水泵額定功率(Kw)! W ___水泵機組慣性矩(N^) TM___最大計算時段(s)! H1___水泵凈揚程(m) Hs___水泵吸出高度(m)! E0___水泵全特性曲線分點間隔度數(shù)(度) J ___水泵全特性曲線分點數(shù)! G ___水的容重(N/m^3) Z1___并聯(lián)水泵臺數(shù)! T1___蝶閥一階段關閉時間(s) T2___蝶閥二階段關閉時間(s) ! T3___蝶閥三階段關閉時間(s) BB___蝶閥二階段關閉完畢時的關閉角(度)! B0___蝶閥開度5度等分點數(shù) Hfr__蝶閥全開(TO=1)設計流量時的水頭損失(m)! Q ___二維數(shù)組中的流量(m^3/s) H ___二維數(shù)組中的揚程(m)! X ___一維數(shù)組中的流量(m^3/s) Z ___一維數(shù)組中的揚程(m)! F ___水泵全特性曲線表中的H/(R^2+V^2) Y ___水泵全特性曲線表中的M/(R^2+V^2)值! R1___停泵水泵的轉速比 V1___停泵水泵的流量比! R2___未停水泵的轉速比 V2___未停水泵的流量比! B1___停泵水泵的轉矩 TF___蝶閥流量系數(shù)! CV___蝶閥水頭損失 ! 定義變量和數(shù)組READ L0,D0,A0,N,C,F0,H0,Q0,N1,P0READ W2,TM,H1,H3,E0,J,G,T1,T3,BB,BO,HFRLET T=L0/(N*C)LET 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) 工 況 下 的 流 量,揚 程 計 算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) 工 況 下 的 流 量,揚 程 計 算FOR K=1 TO NO LET S=0 ! 中 間 斷 面 的 流 量,揚 程 計 算 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 ! 出 口 斷 面 的 流 量,揚 程 計 算 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) ! 水 泵 斷 面 的 流 量,揚 程 計 算 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=。一） 原始數(shù)據(jù)的準備：（1） 水錘波速a平均 鋼管 a=1425/(1+KD/Eδ) 鋼筋混凝土管 a=1425/[1+KD/Eδ(1+)] K—水的彈性模量， E—管材彈性模量，E鋼管=210 Gpa，E鋼混=21 Gpa D—管徑（mm） δ—管壁厚(mm) μ—泊松比 初設δ壓支=5mm， δ并鋼=8mm， δ并鋼混=55mm 在壓力水管的結構計算時進一步驗算壁厚是否能滿足要求 a平均=(l1a1+a2l2)/(l1+l2)=1026米/秒（2） 管路分段數(shù)：4（3） 主管直徑： (4) 管的長度: 309米 (5) 管道摩阻系數(shù): n鋼管= n鋼混管= n=(173+136)/309