鋼架及matlab算法(完整版)
【正文】 2*EIc/h。0 S2 S3。%K12% a(30*(ii+1)+3*jj2:30*(ii+1)+3*jj,30*ii+3*jj2:30*ii+3*jj)=[S1 0 0。0 S3 S4]。S2=12*EIb/(L^3)。0 S2 S3。 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 0 2*S4 ]。S3=6*EIo/(d^2)。0 S3 S4]。for ii=1:3:178R(ii:ii+2,ii:ii+2)=[L/d h/d 0。z(3,1)=q*(h^2)/12。K(:,1:30)=0。x2=zeros(6,1)。M2(ii+1,1)=[0 6*EIc/(h^2) 2*EIc/h 0 6*EIc/(h^2) 4*EIc/h]*x2。0 0 1]。5]。鋼架側(cè)移量（最右側(cè)）39。5 x1(178,1)]。y=[1 x1(58,1)。disp(39。Q(1:5,1)=[1。x3(4:6,1)=x1(88+30*ii:88+30*ii+2,1)。for jj=1:3:6T(jj:jj+2,jj:jj+2)=[0 1 0。disp(x1)。for ii=1:2z(1+30*ii,1)=q*h。*x*R。x(7+30*ii:7+30*ii+2,7+30*ii:7+30*ii+2)=[S1 0 0。0 S2 S3。endK2=e。0 S2 S3。0 0 2*S4 ]。for ii=1:1:5 e(30*ii+1:30*ii+3,30*ii+4:30*ii+6)=[S1 0 0。0 0 1]。0 2*S2 0。0 S3 S5]。 a(30+3*jj2:30+3*jj,30+3*jj2:30+3*jj)=[2*S1 0 0。 S1=EAc/h。斜桿的抗壓剛度EAo =\n39。)。EIc=input(39。五層，九跨，第三跨有斜桿鋼架求解：39。則58，88，118，148，178位移編號所對應(yīng)的位移值為所求點位移。桿件的單元剛度矩陣其中 （3） 建立局部坐標(biāo)系分別建立如圖14所示的豎向分體系局部坐標(biāo)系，水平分體系局部坐標(biāo)系和斜桿分體系坐標(biāo)系。（4） 建立分體系剛度分別建立三個分體系的105105的剛度矩陣，引入循環(huán)變量，依次對相應(yīng)位移剛度賦值。（10） 通過位移求解彎矩將位移帶回到單元方程中，求出對應(yīng)的力，豎直分體系的位移需進(jìn)行一個轉(zhuǎn)置再代入計算最后的。)。柱的抗彎剛度EIc=\n39。EAb=input(39。)。S2=12*EIc/(h^3)。0 2*S2 0。 for ii=1:1:4 a(30*ii+3*jj2:30*ii+3*jj,30*(ii+1)+3*jj2:30*(ii+1)+3*jj)=[S1 0 0。0 0 2*S4 ]。endK
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1