6)"(m)"
230:IF L(1)<=0 OR L(2)<=0 OR L(3)<=0 GOTO 340
235:G1=GM(1):T1=TM(1):S1=SM(1)
240:IF L(1)<L(2) AND L(1)<L(3) THEN L1=L(1) :GOTO 270
250:IF L(2)<L(1) AND L(2)<L(3) THEN L1=L(2) :GOTO 300
260:L1=L(3):L2=0:L3=0:G2=GM(4):T2=TM(4):S2=SM(4) :GOTO 360
270:IF L(4)<=0 OR L(5)<=0 GOTO 330
275:G2=GM(2):T2=TM(2):S2=SM(2)
280:IF L(4)<L(5) THEN L2=L(4): GOTO 310
290:L2=L(5):L3=0:GOTO 360
300:IF L(6)>L1 THEN L2=L(6):G2=GM(3):T2=TM(3):S2=SM(3)
305:L3=0:G3=GM(4):T3=TM(4):S3=SM(4): GOTO 360
310:IF L(6)>L2 THEN L3=L(6):G3=GM(3):T3=TM(3):S3=SM(3)
315:G4=GM(4):T4=TM(4):S4=SM(4) :GOTO 360
320:L3=0:G3=GM(3):T3=TM(3):S3=SM(3): GOTO 360
330:IF L(6)>0 THEN L(2)=L(6):G2=GM(3):T2=TM(3):S2=SM(3):L3=0
335:G3=GM(4):T3=TM(4):S3=SM(4): GOTO 360
340:IF L(4)<=0 OR L(5)<=0 GOTO 354
342:IF L(4)<L(5) THEN L1=L(4):G1=GM(2):T1=TM(2):S1=SM(2) :GOTO 350
344:L1=L(5):G1=GM(2):T1=TM(2):S1=SM(2):L2=0:L3=0: GOTO 360
350:IF L(6)>L1 THEN L(2)=L(6):G2=GM(3):T2=TM(3):S2=SM(3):L3=0
351:G3=GM(4):T3=TM(4):S3=SM(4): GOTO 360
352:L2=0:L3=0:G2=GM(3):T2=TM(3):S2=SM(3): GOTO 360
354:IF L(6)<=0 GOTO 598
356:L1=L(6):G1=GM(3):T1=TM(3):S1=SM(3):L2=0:L3=0
358:G2=GM(4):T2=TM(4):S2=SM(4): GOTO 360
360:LPRINT "L1=";L1;"(m)";
361:LPRINT "L2=";L2;"(m)";
365:LPRINT "L3=";L3;"(m)"
400:print "S1=";S1;"S2="S2;"S3="S4;"S4="S4
410:PRINT "G1=";G1;"G2="G2;"G3="G3;"G4="G4
420:PRINT "T1=";T1;"T2="T2;"T3="T3;"T4="T4
500:DIM A$(10),G(10),T(10),C(L4/K+3),H(L4/K+3)
505:RESTORE 550
510:FOR I=1 TO 10
520:READ A$(I), G(I),T(I)
530:LPRINT A$(I)
540:LPRINT "档距L(m) 应力S(kg/mm2) 弧垂F(m)"
550:DATA "最高温度",3.55E-3,40,"最低温度",3.55E-3,-30,"最大风速",5.03E-3,-5,"最大复冰",4.96E-3,-5
560:DATA "平均温度",3.55E-3,5,"大气过电压<1>:V=0",3.55E-3,15,"大气过电压<2>"
570:DATA 3.59E-3,15,"安装",3.59E-3,-15,"操作过电压",3.74E-3,-15,"事故",3.55E-3
580:DATA -10
590:GOSUB 600
595:NEXT I
598:END
600:GM=G1:TM=T1:SM=S1
610:FOR L=L(0) TO K*INT(L1/K) STEP K
612:GOSUB 700
615:NEXT L
620:L=L1:GOSUB 700
625:GM=G2:TM=T2:SM=S2
630:FOR L=K+K*INT(L1/K) TO 600 STEP K
631:GOSUB 700
635:NEXT L
671:GOSUB 700
675:RETURN
700:G=G(I):T=T(I)
710:A=(L^2*GM^2*E/(24*SM^2))-SM+E*W*(T-TM)
720:B=L^2*G^2*E/24
730:P=-A^2/9:Q=A^3*2/27-B
740:IF Q*Q+4*P^(3)<0 GOTO 850
750:U=((-Q+SQR (Q^2+4*P^(3)))/2)^(1/3)
760:V=((-Q-SQR (Q^2+4*P^(3)))/2)^(1/3)
770:S=U+V-A/3
780:F=L*L*G/8/S
790:S=INT ((S+0.005)*1000)/1000:F=INT ((F+.005)*1000)/1000
800:LPRINT USING "######"; L ;
801:LPRINT USING "##########.##";S;
802:LPRINT USING "###########.##";F
820:RETURN
850:O=ATN ((SQR ((-Q^2-4*P^(3))))/(-Q))
860:IF Q>0 THEN S=(COS (O/3)-SQR (3)*SIN (O/3))*SQR(-P)-A/3
870:IF Q<0 THEN S=2*SQR (-P)*COS(O/3)-A/3
880:GOTO 780
900:Z1=(24/E)*(SM(X)-SM(Y))+24*W*(TM(X)-TM(Y))
902:Z2=(GM(X)/SM(X))^2-(GM(Y)/SM(Y))^2
905:Z=Z1/Z2
910:RETURN
程序二:
800:LPRINT TAB(8) "安装曲线表"
810:K=50
815:X=15
816:E=8490
817:W=18.86E-6
818:PRINT "W=";W;"E=";E
820:DIM X(11), A(11)
830:restore 921
840:FOR I=1 TO 11
850:READ X(I),A(I)
851:L=X(I):SM=A(I)
921:DATA 100,9.01,150,9.21,200,8.88,250,8.48,300,8.17,350,7.93,400,7.76
922:DATA 450,7.63,500,7.53,550,7.46,600,7.40
930:GOSUB 955
935:next I
950:END
955:G=.00355:TM=-30:GM=.00355
980:LPRINT tab(4); "L=";L;"(m)"
985:Lprint "气温T(0C) 应力S(kg/mm2) 弧垂F(m)"
990:for T1=-5 TO 30 STEP 5
995:T=T1-X
1000:T2=T1-5
1710:A=(L^2*GM^2*E/(24*SM^2))-SM+E*W*(T-TM)
1720:B=L^2*G^2*E/24
1730:P=-A^2/9:Q=A^3*2/27-B
1740:IF Q*Q+4*P^(3)<0 GOTO 1850
1750:U=((-Q+SQR (Q^2+4*P^(3)))/2)^(1/3)
1760:V=((-Q-SQR (Q^2+4*P^(3)))/2)^(1/3)
1770:S=U+V-A/3
1780:F=L*L*G/(8*S)
1790:S=INT ((S+0.005)*1000)/1000:F=INT ((F+.005)*1000)/1000
1792:LPRINT USING "######";T2;
1793 LPRINT USING "##########.##";S;F
1794 NEXT T1
1795 IF L=0 GOTO 950
1820:RETURN
1850:O=ATN ((SQR ((-Q^2-4*P^(3))))/(-Q))
1860:IF Q>0 THEN S=(COS (O/3)-SQR (3)*SIN (O/3))*SQR(-P)-A/3
1870:IF Q<0 THEN S=2*SQR (-P)*COS(O/3)-A/3
1880:GOTO 1780
参考文献:1、《水利水电工程PC—1500程序集》
2、《Basic语言》谭浩强 田淑清编著 科学普及出版社
3、《计算机技术与计算机应用》 1995.1
上一页 [1] [2]
|