架空输电线路设计第8章习题答案c语言精确求解各档应力和弧垂 3页

某220kV线路通过我国典型气象区Ⅶ区，导线采用LGJ—300/35，线路的一个耐张段如图所示。假定架线竣工时（初伸长未放出）悬垂串铅垂，气象条件变化引起悬垂串偏斜后各档应力趋于一致，试按连续档编写程序精确求解各档的应力和弧垂。已知悬垂串串长λ=1582mm，串重GJ=438.7N。（导线安全系数取k=2.5，E=78400MPa，α=18.8×10－61/℃）。只要学过输电线路设计的人都应该看得懂前面那些数据，那些作为已知条件直接定义了，学线路的肯定会算。连续档架空线求弧垂，应力(C语言程序代码)#defineE78400#defineri31.514/1000#defineRi31.128/1000#defineQ063.939#definecos10.9974587#definecos20.9960726#definecos30.9937313#definecos40.9950372#definet1.582#defineGj438.7#defineA333.31#include#includevoidmain(){doublel[4];/*绝缘子的偏移量*/doubleQ[4];/*应力*/doubleLP[4];/*档距的变化量*/floatL1=350.0,L2=450.0,L3=400.0,L4=300.0;/*档距*/doublef[4];/*弧垂*/floatH1=25.0,H2=40.0,H3=45.0,H4=30.0;Q[0]=48.0;do{LP[0]=(pow(L1,2)*pow(cos1,2)/24*(pow(ri/Q0,2)-pow(Ri/Q[0],2))+(Q[0]-Q0)/E/cos1+9.4/10000)*L1/pow(cos1,2)/(1+pow((L1*Ri),2)/(8*Q[0]*Q[0]));l[0]=LP[0]; Q[1]=(l[0]*(Gj/2/A+Ri*L1/2/cos1+Q[0]*H1/L1+Ri*L2/2/cos2)+Q[0]*sqrt(t*t-l[0]*l[0]))/(sqrt(t*t-l[0]*l[0])+l[0]*H2/L2);LP[1]=(pow(L2,2)*pow(cos2,2)/24*(pow(ri/Q0,2)-pow(Ri/Q[1],2))+(Q[1]-Q0)/E/cos2+9.4/10000)*L2/pow(cos2,2)/(1+pow((L2*Ri),2)/(8*Q[1]*Q[1]));l[1]=LP[1]+LP[0];Q[2]=(l[1]*(Gj/2/A+Ri*L2/2/cos2+Q[1]*H2/L2+Ri*L3/2/cos3)+Q[1]*sqrt(t*t-l[1]*l[1]))/(sqrt(t*t-l[1]*l[1])+l[1]*H3/L3);LP[2]=(pow(L3,2)*pow(cos3,2)/24*(pow(ri/Q0,2)-pow(Ri/Q[2],2))+(Q[2]-Q0)/E/cos3+9.4/10000)*L3/pow(cos3,2)/(1+pow((L3*Ri),2)/(8*Q[2]*Q[2]));l[2]=LP[2]+LP[1]+LP[0];Q[3]=(l[2]*(Gj/2/A+Ri*L3/2/cos3+Q[2]*H3/L3+Ri*L4/2/cos4)+Q[2]*sqrt(t*t-l[2]*l[2]))/(sqrt(t*t-l[2]*l[2])+l[2]*H4/L4);LP[3]=(pow(L4,2)*pow(cos4,2)/24*(pow(ri/Q0,2)-pow(Ri/Q[3],2))+(Q[3]-Q0)/E/cos4+9.4/10000)*L4/pow(cos4,2)/(1+pow((L4*Ri),2)/(8*Q[3]*Q[3]));l[3]=LP[3]+LP[2]+LP[1]+LP[0];Q[0]=Q[0]+0.01;}while(fabs(l[3])>0.001);f[0]=Ri*pow(L1,2)/8/Q[0]/cos1;f[1]=Ri*pow(L2,2)/8/Q[1]/cos2;f[2]=Ri*pow(L3,2)/8/Q[2]/cos3;f[3]=Ri*pow(L4,2)/8/Q[3]/cos4;printf("第一档绝缘子的偏移量L1=%fn",l[0]);printf("第二档绝缘子的偏移量L2=%fn",l[1]);printf("第三档绝缘子的偏移量L3=%fn",l[2]);printf("第四档绝缘子的偏移量L4=%fn",l[3]);printf("*******************************nn");printf("第一档应力q1=%fn",Q[0]);printf("第二档应力q2=%fn",Q[1]);printf("第三档应力q3=%fn",Q[2]);printf("第四档应力q4=%fn",Q[3]);printf("*******************************nn");printf("第一档档距的变化量LP1=%fn",LP[0]);printf("第二档档距的变化量LP2=%fn",LP[1]);printf("第三档档距的变化量LP3=%fn",LP[2]);printf("第四档档距的变化量LP4=%fn",LP[3]);printf("*******************************nn");printf("第一档弧垂f1=%fn",f[0]);printf("第二档弧垂f2=%fn",f[1]);printf("第三档弧垂f3=%fn",f[2]);printf("第四档弧垂f4=%fn",f[3]);printf("*******************************nn");}