摘要：本文对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) 凭据以上公式算出各类控制条件下的应力Ｓ，然后又凭据弧垂公式Ｆ＝ＧＬ２／８Ｓ，得出各类控制条件下的弧垂Ｆ。这样就可计较各类代表档距下在分歧温度时的电线应力、弧垂（即安装曲线） 而在ＰＣ－１５００机输电线路力学计较法式编写中把判别式编写错了， 即：时， (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