导线应力弧垂计算的BASIC程序

导线应力弧垂计算的BASIC程序

一、前言

架空线路设计和施工都需要进行导线力学计算．笔者编制了导线应力、弧垂计算的BASIC程序，用户只需按屏幕显示的表格键入导线参数、气象条件，计算机即能完成计算全过程，并将计算结果打印制表。各种计算项目采用菜单选择，用户使用非常方便。本文就该程序的设计方法及特点作一简单介绍，以供参考．二、架空导线应力、孤垂的计算机算法

1．导线比载计算

导线的综合比载是垂直比载（自重、冰重）、水平比载（风压）的矢量和．对各种气象情况的综合比载可用下式表示：

式中：q——导线的单位重量（千克／千米）

S——导线的计算截面（毫米2）

d——导线的计算外径（毫米）

b——导线覆冰厚度（毫米）

v——设计风速（米／秒）

C——风荷载体形系数，当线径d＜17毫米时，C=1.2，当线径≥17毫米时，C=1.1；覆冰时不论线径大小C=1.2

α——风速不均匀系数，根据不同风速取值。（程序框图略）

2．临界档距计算及有效临界档距判别

根据工程需要，导统应力孤垂的计算项目有时多达十种，即最大风速、覆冰情况、安装情况、事故断线、最低气温、最高气温、外过电压（有风、无风）、内过电压、平均气温。这十种情况对应十种气象条件．但导线选用应力的控制条件只可能是其中的4种情况，即最低气温、最大风速、覆冰情况和平均气温．这4种控

制条件的两两组合有6个临界档距。一般地n种控制条件有=n(n－1)／2个

临界档距，其中有效临界档距有0～（n—1）个。两个控制条件的临界档距为

式中：E——导线弹性模数（千克／毫米２）

a——导线温度线膨胀系数（l／℃）

δi、δj——两种控制条件的限定应力（最大使用应力或年平均运行应力上限）（千克／毫米２）

ti、tj——两种控制条件的气温（℃）

gi、gj——两种控制条件的比载（千克／米·毫米２）。

由式（2－1）可知，若将n个控制条件的g／δ值由小到大排列，再比较各δ

＋ aEt，并满足下式：

不满足式（2－2）的控制条件不起作用舍去。当两种控制条件的 g／δ相同时，舍去δ＋ aEt较大者；若两者的δ＋ aEt相同，舍去g／δ较小者，则所有满足式（2－2）的控制条件均有实数解的临界档距，把满足（2－2）式的控制条件由小到大编为序号1、2、3、…c（c≤n），并相应建立C－l个临界档距数栏。

由式（2－l）按序号依次求出控制条件i与其后的i＋1、i＋2、i＋3、…，C

控制条件间的临界档距li－(i+1)、li－(i+2)、…li－c，并填入i栏。首先从第1栏开始，找出其中的最小临界档距l1－i，即是序号为1、i控制条件的有效临界档距。其中1条件控制的档距区间为〔0,li－j〕,舍去1～i序号间的控制条件，这些条件不起控制作用．再从第i栏中找出最小者li－j即为 i、j控制条件的有效临界档距，i条件控制的档距区间为[l1－i,li－j序号间的控制条件舍去。如此直至求出有效临界档距lK－C，则C条件控制的档距区间为[lK－C,∞]。

上程序设计框图：

当只需要计算几个档距的应力、孤垂时，则可不必求解临界档距，直接用下式进行判断：

将n个控制条件分别代人（2－3）式，求出最大者Di，i控制条件即是该档距的有效控制条件。

3．导线应力、弧垂计算

导线状态方程为

式中：δm、δn——分别为己知控制条件的限定应力和待求条件的应力（千

克／毫米2）；

gm、gn——分别为已知控制条件的比载和待求气象条件下的比我

千克／米·毫米2）；

tm、tn——分别为已知控制条件和待来条件的气温（℃）；

a——导线温度线膨胀系数（ l／℃）；

E——导线弹性模数（千克／毫米2）；

l ——档距（米）．

（3－1）式可改写为一般形式的方程式：

用迭代法解导线状态方程（3－2），求出应力δ，即可解得导线弧垂了：

解导线状态方程除用以上方法外，还可用代数法求解精确根，其方法如下：