杆塔塔腿内力计算
clear;
clc;
b=4.070;
G_D=38627;
G_B=4587.5;
P_B=8917.52;
P_D=20106;
hb=44.450;
hd=39.950;
w1=101;
w2=31;
w3=6.2;
w4=7.16;
w5=10.96;
w6=12.24;
w7=13.96;
w8=14.98;
h1=53.50;
h2=42.15;
h3=34.45;
h4=31.45;
h5=27.95;
h6=23.41;
h7=18.85;
h8=14.05;
x1=(2*P_B*hb/2+3*P_D*hd/2+w1*h1/2+w2*h2/2+w3*h3/2+w4*h4/2+w5*h5/2+w6*h6/2+w7*h7/2+w8*h8/2)/b;
x2=(2*G_B/2+3*G_D/2+18254/2)/cos(2*pi*5/360);
u1=(x1+x2)/2
u2=(x1-x2)/2
b2=21.973;
h0=46.294;
x3=(2*P_B*(h0-hb)/2+3*P_D*(h0-hd)/2+w1*(h0-h1)/2+w2*(h0-h2)/2+w3*(h0-h3)/2+w4*(h0-h4)/2+w5*(h0-h5)/2+w6*(h0-h6)/2+w7*(h0-h7)/2+w8*(h0-h8)/2)/b2;
s1=x3/2;
s2=-1*x3/2;
hx=u1*cos(2*pi*85/360)+s1*cos(2*pi*62/360)
N=u1*sin(2*pi*85/360)+s1*sin(2*pi*62/360)
hx2=u2*cos(2*pi*85/360)-s2*cos(2*pi*62/360)
T=u2*sin(2*pi*85/360)-s2*sin(2*pi*62/360)