输电杆塔导线地线比载计算,应力计算

比载计算
clear;
clc;
%绝缘子串的片数选择
alph=(1.82+2.27)/2;%爬电比距
Un=500;%标称电压
h=40;%爬电距离cm
ns=0;%设计变化绝缘子数
n_j=ceil(alph*Un/h)+ns;%绝缘子数
hj=15.5;%绝缘子的公称高度cm
ljy=n_j*hj;%悬垂绝缘子总长度cm
%导线的比载计算
q=1349;%单重kg/km
g=9.80665;%重力加速度m/s2
A=425.24;%总截面
r1=q*g/A*1.0e-3;%自重比载mpa/m
b=5;%覆冰厚度mm
d=26.82;%架空线外径mm
V=pi*b*(d+b);%覆冰体积mm
r2=27.728*b*(d+b)/A*1.0e-3;%冰重比载
r3=r1+r2;%垂直总比载
v0_10=10;%安装风速
v0_35=35;%最大风速
v0_15=15;%雷电过电压
v0_18=18;%操作过电压
hf=42;%500kv导线平均高度20m
beita=1.000;%修正系数
z=0.15;%粗糙度指数
v_10=beita*(hf/10)^z*v0_10;
v_35=beita*(hf/10)^z*v0_35;
v_15=beita*(hf/10)^z*v0_15;
v_18=beita*(hf/10)^z*v0_18;
baitac_xy20=1.00;%风载调整系数
baitac_dy20xy27=1.10;
baitac_dy27xy315=1.20;
baitac_dy315=1.30;
alph_xy20_qd=1.00;%风速不均匀系数(强度)
alph_dy20xy27_qd=0.85;
alph_dy27xy315_qd=0.75;
alph_dy315_qd=0.70;
alph_xy20_fp=1.00;%风速不均匀系数(风偏)
alph_dy20xy27_fp=0.75;
alph_dy27xy315_fp=0.61;
alph_dy315_fp=0.61;
usc=1.1;%架空线体形系数
theat=pi/2;%风向与线路的夹角
r4_10_qd=baitac_xy20*alph_xy20_qd*usc*d*0.625*v_10^2/A*sin(theat)^2*1.0e-3;%r4(0,10)强度
r4_10_fp=baitac_xy20*alph_xy20_fp*usc*d*0.625*v_10^2/A*sin(theat)^2*1.0e-3;%r4(0,10)风偏
r4_35_qd=baitac_dy315*alph_dy315_qd*usc*d*0.625*v_35^2/A*sin(theat)^2*1.0e-3;%r4(0,35)强度
r4_35_fp=