摘要:本文对PC-1500机上高压送电线路电线力学计较法式进行了修改,并对把PC-1500机上的法式移植到微机上的方式和技巧进行了一些说明。
关头词:力学计较法式电线状态方程式修改移植TurboBasic语言
PC-1500机由于它功能比力齐全,开发容易,使用简洁,体积小,投资少,易学,生效快,在80年月比力流行,但随着计较机技术的成长和PC-1500机的自然消耗,微机将慢慢取代PC-1500机。近年来各部门在PC-1500机上开发了许多秀的法式,取得了可喜的功效和效益,若是把这些秀法式充实哄骗起来,把它们移植到微机上加以使用,就能够加速计较机的运用步伐,取得事倍功半的效果。 电线力学计较是送电线路设计中复杂的一项工作,其内容包括:1、计较各控制条件之间的临界档距:2、判别有用临界档距;3、计较各类气象条件和分歧档距下的电线应力、弧垂;4、计较各类代表档距下在分歧气温时的电线应力、弧垂(即安装曲线)。这些工作完成以后,方能列出线路杆塔的荷重条件,继而进行杆塔结构计较。 本法式分法式1、法式二,法式一解决上述工作1、2、3;法式二解决工作4。 PC-1500机上电线力学法式的主要计较公式的解法及法式的毛病。 电线状态方程式:就是从已知气象条件下的应力,求另外一条件下应力的计较式,悬挂于两固定点间的导线,当气象条件发生变化时,导线的应力亦将随之变化。设已知气象条件下的温度为Tm,比载为Gm,此时应力为Sm,L为变(档距),令待求气象条件下的温度是T,比载是Gn,应力是S,若何求得S呢?可用下面的导线的状态方程式对S求解 (1) 令上式中: (2) 把(2)代进(1)得: (3) 式中S为待求应力,注脚m暗示已知条件下的数值,注脚n暗示待求条件下的数值。 人工计较时,先求出与每L相对应的A、B值,然后用试算法求解上述三次方程。本方程系采用直接求解法,故将方程作以下变换: 对于(4) 令(5) 原方程化为(6) 又令(7) 得:(8) 因而(9) 据韦达定理,为二次方程:的两个根(10) 故:(11) 当时 (12) 当时,令: (13) 时, (14) 时, (15) 凭据以上公式算出各类控制条件下的应力S,然后又凭据弧垂公式F=GL2/8S,得出各类控制条件下的弧垂F。这样就可计较各类代表档距下在分歧温度时的电线应力、弧垂(即安装曲线) 而在PC-1500机输电线路力学计较法式编写中把判别式编写错了, 即:时, (16) 时, (17) 因判别式遍写错,故法式得出的应力(S)、弧垂(F)有误,计较出的安装曲线误差更年夜,不能用于设计,本法式不能用。我对这一法式进行了修改并把它成功地移植到微机中。在移植中,因PC-1500机上的Basic语言比力接近微机上的一般Basic语言,而TurboBasic语言与一般Basic语言相比,是向下兼容的,即TurboBasic语言包括一般Basic语言,并所扩大,下面就移植进程中需要注重的一些问题作一些介绍。两种机型的Basic语言几近是不异的,只是在语句中,微机的语句界说符(如IF、THEN、FOR、TO、STEP等)与前后部门的内容一般都用空格来分格开,以便与其它成份进行区分,而PC-1500机上没有这类要求,这是首先要注重的问题。 其次,在PC-1500机中,函数后所跟的自变可以不用圆括号括起来,而在微机上函数后所跟的自变必需用圆括号括起来,这也是必需引发注重的地方。 此外,常、变、语句、函数中,也有一些差异需求注重,下面划分进行说明。 常:PC-1500机上有PI(丌)常;微机上没有,只能用“3.1415927”这个具体的数来暗示。变:PC-1500机的变用一个英文字母、或二个英文字母的组合,或一个英文字母加上一个数字的组合组成的字符串来暗示;而微机则是作英文字母开首的字母串或字母数字串来暗示变名,这里需要说明的是,与语言界说符不异的字符串不能做变名。赋值语句:LET…,PC-1500机在多个赋值语句之间可以用逗号“,”分隔,而微机上只能以冒号“:”分隔。在微机上,当赋值号“=”双方的数据类型分歧时,自动转换成赋值号左侧变的数据类型。键盘输进语句:INPUT…,PC-1500机的INPUT语句中可以包括多个“说明”,用来说明每一个期待输进的变,微机上则只能有一个说明。对于PC-1500的输进语句,INPUT“说明”;变,“说明”;变,…的形式,在微机上可以改写为多个INPUT“说明”;变语句来实现,也能够写为一个INPUT“说明”;变,变,…语句的形式。输出语句:PRINT/LPRINT;输出语句分屏幕显示输出语句PRINT,打印机打印输出语句LPRINT。PC-1500机显示输出语句PRINT会中断法式运行,一般需要与PAUSE语句、或WAIT语句配合使用,而微机的PRINT语句不会中断法式运行,则可以把PC-1500机中的PAUSE语句或WAIT语句直接删失落。在用USING界说的花式输出语句中,PC-1500机可以把USING作为一个界说花式的语句零丁使用,而微机上则不能零丁作为语句使用,只能配合PRINT(LPRINT)使用,即:PRINTUSING“花式”;变,变,…。在输出语句中使用了TAB位置函数时,PC-1500机TAB后所跟的位置参数可以不用括号,而微机上TAB后所跟的位置参数必需用圆括号括起来,即要把TAB15改写为TAB(15)。清屏语句:CLS;PC-1500机屏幕只有一行,显示输出成效和信息时,需屡次频仍地使用清屏语句,而微机屏幕共25行,一般情况下,只需在法式开首使用,PC-1500机法式中过剩的清屏语句楞直接删除失落。条件语句:IF…THEN…;PC-1500机上IF…LET形式的语句,若是要在微机上使用,需改成IF…THEN或IF…THENLET语句形式。轮回语句:FOR…NEXT;两种机械语句形式基底细同,只是在执行FOR语句时,微机是先判别后轮回,而PC-1500机是先轮回,后判别。转移语句:GOTO…;PC-1500机GOTO语句后可所以表达式,而微机上只能跟行号。转子语句GOSUB以后在PC-1500机上可所以表达式或肯定的行号。而微机上只能跟行号。 PC-1500机随机配CE-150绘图式打印机,具有绘图、打印两种功能,其绘图语句和部门打印语句与该打印有关,一般不能直接用于其它计较机上。由于微机用打印机绘图时,通常为采用拷贝屏幕的方式实现,是以,对于绘图语句,在移植法式时,由于实现方式分歧,好用注释语句注释起来,或直接删除的法子处置。等法式初步伐通后,再进一步斟酌绘图功能的实现。 法式移植技巧 把PC-1500机法式移植到微机上使用,一般需要经过以下几个步骤:首先将PC-1500机与微机联机,经由过程RC-232串行口输进、或手工键盘输进,把PC-1500机法式以ASCⅡ码形式寄存在微机上,然后逐行检查法式,对需要注释或删除的语句(如显示绘图语句、绘图打印语句)进步前辈行处置,对两种机械上有差异的语句、函数进行修改(如语句界说符加空格、函数加括号)。在这些工作做完后,就能够起头调试法式。在调试进程中,多输出一些中心成效,以便检查犯错缘由,直到算出准确的成效。 微机法式在初步伐试经由过程后,可凭据微机的特点,扩充法式功能,完善打印输出,增加和设计用户界面,直到编写出使人满意的法式为止。 总之,只要弄清晰了上面所叙述的两种计较机之间的语言差异,并在法式移植进程中加以注重,PC-1500机法式很容易就能够移植到微机上使用。本高压送电线路电线力学计较法式经修改、移植到微机中,现也经由过程算例,并可以用于设计了。 附修改、移植后的法式以下: 法式一: 1:CLS 2:LPRINTTAB(10);"高压送电线路力学计较法式计较成效" 3:INPUT"电线型号B$=";B$ 4:INPUT"电线标称截面C$=";C$ 6:LPRINTTAB(4);B$;"-"C$ 8:LETK=50 9:LETD=1 10:DIML(6),GM(4),TM(4),SM(4) 15:INPUT" 电线弹性模E=";E 16:LPRINTTAB(4);" 电线弹性模=";E 17:INPUT" 电线温度线胀系数W="W 18:LPRINTTAB(4);" 电线温度线胀系数="W 20:DATA3.55E-3,-30,11.2,4.96E-3,-5,11.2,5.03E-3,-5,11.2,3.55E-3,5,7 25:INPUT"肇端计较档距L(0)="L(0) 26:INPUT"年夜计较档距L(4)="L(4) 35:LPRINTTAB(4);" 肇端计较档距=";L(0);" 年夜计较档距=";L(4) 40:FORI=1TO4 50:READGM(I),TM(I),SM(I) 60:NEXTI 70:FORI=1TO3 80:X=1:Y=I 1:GOSUB900 90:IFZ<0GOTO110 100:L(I)=INT(.5 SQR(Z)) 105:NEXTI:GOTO120 110:L(I)=-1:GOTO105 120:FORI=4TO5 130:X=2:Y=I-1:GOSUB900 140:IFZ<0GOTO160 150:L(I)=INT(.5 SQR(Z)) 155:NEXTI:GOTO170 160:L(I)=-1:GOTO155 170:X=3:Y=4:GOSUB900 180:IFZ<0GOTO200 190:L(6)=INT(.5 SQR(Z)):GOTO210 200:L(6)=-1 210:FORI=1TO6 220:LPRINTTAB(4);"L(";I;")=";L(I);"(m)":NEXTI 225:LPRINT"L(AB)=";L(1)"(m)";"L(AC)=";L(2)"(m)";"L(AD)=";L(3);"(m)"; 227:LPRINT"L(BC)=";L(4)"(m)";"L(BD)=";L(5)"(m)";"L(CD)=";L(6)"(m)" 230:IFL(1)<=0ORL(2)<=0ORL(3)<=0GOTO340 235:G1=GM(1):T1=TM(1):S1=SM(1) 240:IFL(1)<L(2)ANDL(1)<L(3)THENL1=L(1):GOTO270 250:IFL(2)<L(1)ANDL(2)<L(3)THENL1=L(2):GOTO300 260:L1=L(3):L2=0:L3=0:G2=GM(4):T2=TM(4):S2=SM(4):GOTO360 270:IFL(4)<=0ORL(5)<=0GOTO330 275:G2=GM(2):T2=TM(2):S2=SM(2) 280:IFL(4)<L(5)THENL2=L(4):GOTO310 290:L2=L(5):L3=0:GOTO360 300:IFL(6)>L1THENL2=L(6):G2=GM(3):T2=TM(3):S2=SM(3) 305:L3=0:G3=GM(4):T3=TM(4):S3=SM(4):GOTO360 310:IFL(6)>L2THENL3=L(6):G3=GM(3):T3=TM(3):S3=SM(3) 315:G4=GM(4):T4=TM(4):S4=SM(4):GOTO360 320:L3=0:G3=GM(3):T3=TM(3):S3=SM(3):GOTO360 330:IFL(6)>0THENL(2)=L(6):G2=GM(3):T2=TM(3):S2=SM(3):L3=0 335:G3=GM(4):T3=TM(4):S3=SM(4):GOTO360 340:IFL(4)<=0ORL(5)<=0GOTO354 342:IFL(4)<L(5)THENL1=L(4):G1=GM(2):T1=TM(2):S1=SM(2):GOTO350 344:L1=L(5):G1=GM(2):T1=TM(2):S1=SM(2):L2=0:L3=0:GOTO360 350:IFL(6)>L1THENL(2)=L(6):G2=GM(3):T2=TM(3):S2=SM(3):L3=0 351:G3=GM(4):T3=TM(4):S3=SM(4):GOTO360 352:L2=0:L3=0:G2=GM(3):T2=TM(3):S2=SM(3):GOTO360 354:IFL(6)<=0GOTO598 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):GOTO360 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:DIMA$(10),G(10),T(10),C(L4/K 3),H(L4/K 3) 505:RESTORE550 510:FORI=1TO10 520:READA$(I),G(I),T(I) 530:LPRINTA$(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:DATA3.59E-3,15,"安装",3.59E-3,-15,"操作过电压",3.74E-3,-15,"事故",3.55E-3 580:DATA-10 590:GOSUB600 595:NEXTI 598:END 600:GM=G1:TM=T1:SM=S1 610:FORL=L(0)TOK*INT(L1/K)STEPK 612:GOSUB700 615:NEXTL 620:L=L1:GOSUB700 625:GM=G2:TM=T2:SM=S2 630:FORL=K K*INT(L1/K)TO600STEPK 631:GOSUB700 635:NEXTL 671:GOSUB700 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:IFQ*Q 4*P^(3)<0GOTO850 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:LPRINTUSING"######";L; 801:LPRINTUSING"##########.##";S; 802:LPRINTUSING"###########.##";F 820:RETURN 850:O=ATN((SQR((-Q^2-4*P^(3))))/(-Q)) 860:IFQ>0THENS=(COS(O/3)-SQR(3)*SIN(O/3))*SQR(-P)-A/3 870:IFQ<0THENS=2*SQR(-P)*COS(O/3)-A/3 880:GOTO780 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:LPRINTTAB(8)"安装曲线表" 810:K=50 815:X=15 816:E=8490 817:W=18.86E-6 818:PRINT" W=";W;" E=";E 820:DIMX(11),A(11) 830:restore921 840:FORI=1TO11 850:READX(I),A(I) 851:L=X(I):SM=A(I) 921:DATA100,9.01,150,9.21,200,8.88,250,8.48,300,8.17,350,7.93,400,7.76 922:DATA450,7.63,500,7.53,550,7.46,600,7.40 930:GOSUB955 935:nextI 950:END 955:G=.00355:TM=-30:GM=.00355 980:LPRINTtab(4);" L=";L;"(m)" 985:Lprint"气温T(0C)应力S(kg/mm2)弧垂F(m)" 990:forT1=-5TO30STEP5 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:IFQ*Q 4*P^(3)<0GOTO1850 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:LPRINTUSING"######";T2; 1793LPRINTUSING"##########.##";S;F 1794NEXTT1 1795IFL=0GOTO950 1820:RETURN 1850:O=ATN((SQR((-Q^2-4*P^(3))))/(-Q)) 1860:IFQ>0THENS=(COS(O/3)-SQR(3)*SIN(O/3))*SQR(-P)-A/3 1870:IFQ<0THENS=2*SQR(-P)*COS(O/3)-A/3 1880:GOTO1780 参考文献:1、《水利水电工程PC—1500法式集》 2、《Basic语言》谭浩强田淑清编著科学普及出书社 3、《计较机技术与计较机运用》1995.1