架空输电线路设计第9章习题3答案C语言精确求解断线张力.docx

某110kV架空线路，导线为LGJ?95/20，某耐张段共有10档，无高差，档距基本相等为lr=250 m，在档距l7内跨越Ⅰ级通讯线，靠输电线路档距中央侧的通讯线高7m，位于距6号杆40 m处，如图所示。直线杆塔悬点高14m，挠度系数B=0.0003m/N，悬垂串长=1.322m，重393N。在气温+15℃、无风时的导线应力为σ0=66MPa，试核验邻档断线后的交叉垂直距离。题3 耐张段示意其它数据，都直接定义了，学过架空输电线路设计的人肯定会算，我就不再列出。C语言求断线张力代码 #define L 250.0 #define t 1.322 #define B 0.0003 #define Gj 393.0 #define A 113.96 #define r 35.1873/1000 #define Pi r*A*L //导线作用在悬挂点的垂直荷载 #define E 76000.0 #define T0 7521.36 //断线前的水平张力 #define p 4.0099 //架空线的荷载 #includestdio.h #includemath.hvoid main (){ double T[5]; //第i档架空线的水平张力 double l[4]; //第i档距的变化量 double k[4]; //绝缘子的偏移量 double x,x1,eps=1e-4,f,f1; //求解方程所需的变量，误差为eps double o,fx;//断线应力，弧垂,T[0]=4939.78;do{ l[0]=(pow(p,2)*pow(L,2)/24*(1/T[0]/T[0]-1/T0/T0)+(T0-T[0])/E/A)*L/(1+p*p*L*L/8/T[0]/T[0]);k[0]=l[0]; x=1.0; //x=1.0是初值do {x1=x; f=k[0]-t*(T[0]-x1)/sqrt(pow(Pi+Gj/2,2)+pow(T[0]-x1,2))-B*(T[0]-x1); //f为f(x)函数 f1=t*(pow((Pi+Gj/2),2)+pow((T[0]-x1),2)-pow((T[0]-x1),2))/pow(pow((Pi+Gj/2),2)+pow((T[0]-x1),2),1.5)+B; x=x1-f/f1; f=k[0]-t*(T[0]-x)/sqrt(pow(Pi+Gj/2,2)+pow(T[0]-x,2))-B*(T[0]-x);}while(fabs(f)=eps fabs(x-x1)=eps);T[1]=x; l[1]=(pow(p,2)*pow(L,2)/24*(1/T[1]/T[1]-1/T0/T0)+(T0-T[1])/E/A)*L/(1+p*p*L*L/8/T[1]/T[1]);k[1]=l[0]+l[1]; x=1.0; do {x1=x; f=k[1]-t*(T[1]-x1)/sqrt(pow(Pi+Gj/2,2)+pow(T[1]-x1,2))-B*(T[1]-x1); f1=t*(pow((Pi+Gj/2),2)+pow((T[1]-x1),2)-pow((T[1]-x1),2))/pow(pow((Pi+Gj/2),2)+pow((T[1]-x1),2),1.5)+B; x=x1-f/f1; f=k[1]-t*(T[1]-x)/sqrt(pow(Pi+Gj/2,2)+pow(T[1]-x,2))-B*(T[1]-x);}while(fabs(f)=eps fabs(x-x1)=eps); T[2]=x; l[2]=(pow(p,2)*pow(L,2)/24*(1/T[2]/T[2]-1/T0/T0)+(T0-T[2])/E/A)*L/(1+p*p*L*L/8/T[2]/T[2]);k[2]=l[0]+l[1]+l[2]; x=1.0; do { x1=x; f=k[2]-t*(T[2]-x1)/sqrt(pow(Pi+Gj/2,2)+pow(T[2]-x1,2))-B*(T[2]-x1);