鋼架及matlab算法(編輯修改稿)
2024-09-19 19:23 本頁面
【文章內(nèi)容簡介】 :30*ii+6,30*ii+4:30*ii+6)=[2*S1 0 0。0 2*S2 0。0 0 2*S4 ]。 e(30*ii+4:30*ii+6,30*ii+1:30*ii+3)=[S1 0 0。0 S2 S3。0 S3 S5]。 for jj=2:1:9 e(30*ii+3*jj2:30*ii+3*jj,30*ii+3*(jj+1)2:30*ii+3*(jj+1))=[S1 0 0。0 S2 S3。0 S3 S5]。 e(30*ii+3*(jj+1)2:30*ii+3*(jj+1),30*ii+3*jj2:30*ii+3*jj)=[S1 0 0。0 S2 S3。0 S3 S5]。 e(30*ii+3*(jj+1)2:30*ii+3*(jj+1),30*ii+3*(jj+1)2:30*ii+3*(jj+1))=[2*S1 0 0。0 2*S2 0。0 0 2*S4 ]。 end e(30+30*ii2:30+30*ii,30+30*ii2:30+30*ii)=[S1 0 0。0 S2 S3。0 S3 S4]。endK2=e。 d=(h^2+L^2)^(1/2)。S1=EAo/d。S2=12*EIo/(d^3)。S3=6*EIo/(d^2)。S4=4*EIo/d。S5=2*EIo/d。for ii=0:1:3x(7+30*ii:7+30*ii+2,40+30*ii:40+30*ii+2)=[S1 0 0。0 S2 S3。0 S3 S5]。x(40+30*ii:40+30*ii+2,40+30*ii:40+30*ii+2)=[S1 0 0。0 S2 S3。0 S3 S4]。x(40+30*ii:40+30*ii+2,7+30*ii:7+30*ii+2)=[S1 0 0。0 S2 S3。0 S3 S5]。x(7+30*ii:7+30*ii+2,7+30*ii:7+30*ii+2)=[S1 0 0。0 S2 S3。0 S3 S4]。endR=zeros(180,180)。for ii=1:3:178R(ii:ii+2,ii:ii+2)=[L/d h/d 0。h/d L/d 0。0 0 1]。endK3=R39。*x*R。K=K1+K2+K3。 F0=zeros(180,1)。z=zeros(180,1)。z(3,1)=q*(h^2)/12。z(153,1)=q*(h^2)/12。z(1,1)=*q*h。z(151,1)=*q*h。for ii=1:2z(1+30*ii,1)=q*h。endF1=z。F=F0+F1。 K(1:30,:)=0。K(:,1:30)=0。K(1:30,1:30)=eye(30)。F(1:30,1)=0。x1=K\F。disp(x1)。 M1=zeros(5,1)。for ii=1:5M1(ii,1)=[0 6*EIb/(L^2) 2*EIb/L 0 6
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1