LGJ-150/35型导线应力弧垂计算程序（c语言）
main() {int l;
float g1=0.0365117,f,g2,A,B,t1=15,t2,c=85.025,x,x1,x0,dx; printf(\"请输入第二气象状态下的比载g2,温度t2:\\n\"); scanf(\"%f%f\for(l=50;l<=800;l++) {if(l%50==0) {A=c-(80000*g1*g1*l*l)/(24*c*c)-80000*0.0000178*(t2-10); B=80000*g2*g2*l*l/24; x0=1.0; do{ dx=(x0*x0*x0-A*x0*x0-B)/(3*x0*x0-2*A*x0); x1=x0; x0=x0-dx; }while(x0!=x1); //printf(\"档距为%d\ printf(\"档距l=%dm时,求的应力σ=%fMpa/m,\ f=(g2*l*l)/(8*x0); printf(\"求的弧垂f=%fm;\\n\ } }
例如最高气温气象条件下的计算结果如下：比载=0.0365117，温度=40℃：
请输入第二气象状态下的比载g2,温度t2: 最高气温： 0.0365117 40
档距l=50m时,求的应力σ=46.010757Mpa/m,求的弧垂f=0.247983m; 档距l=100m时,求的应力σ=52.360615Mpa/m,求的弧垂f=0.8710m; 档距l=150m时,求的应力σ=58.093681Mpa/m,求的弧垂f=1.7678m; 档距l=200m时,求的应力σ=62.791565Mpa/m,求的弧垂f=2.907373m; 档距l=250m时,求的应力σ=66.562653Mpa/m,求的弧垂f=4.285401m; 档距l=300m时,求的应力σ=69.579750Mpa/m,求的弧垂f=5.903393m; 档距l=350m时,求的应力σ=71.999672Mpa/m,求的弧垂f=7.765111m; 档距l=400m时,求的应力σ=73.950539Mpa/m,求的弧垂f=9.874627m; 档距l=450m时,求的应力σ=75.533440Mpa/m,求的弧垂f=12.235672m; 档距l=500m时,求的应力σ=76.827019Mpa/m,求的弧垂f=14.851424m; 档距l=550m时,求的应力σ=77.2105Mpa/m,求的弧垂f=17.724501m; 档距l=600m时,求的应力σ=78.775711Mpa/m,求的弧垂f=20.857019m; 档距l=650m时,求的应力σ=79.514259Mpa/m,求的弧垂f=24.250672m; 档距l=700m时,求的应力σ=80.136055Mpa/m,求的弧垂f=27.906810m; 档距l=750m时,求的应力σ=80.663231Mpa/m,求的弧垂f=31.826508m; 档距l=800m时,求的应力σ=81.113167Mpa/m,求的弧垂f=36.010628m; Press any key to continue