鋼架及matlab算法-預(yù)覽頁
【正文】 碼以桿件為單元，將結(jié)構(gòu)拆分，建立整體坐標(biāo)系，并對節(jié)點(diǎn)位移按圖12所示編碼。（5） 坐標(biāo)轉(zhuǎn)換對豎向坐標(biāo)系和斜桿體系進(jìn)行轉(zhuǎn)置，其坐標(biāo)轉(zhuǎn)換陣為其中。（8） 邊界條件處理因?yàn)橹着c大地剛接，所以支座處不發(fā)生位移，利用置零置一法對剛度矩陣進(jìn)行處理，使位移滿足條件。 二． 說明程序本程序使用matlab軟件通過矩陣位移法對題目進(jìn)行編程設(shè)計(jì)，按照提示輸入柱子的抗彎剛度EIc；柱子的抗壓剛度EAc；梁的抗彎剛度EIb；梁的抗壓剛度 EAb；斜桿的抗彎剛度EIo；斜桿的抗壓剛度EAo；側(cè)向均布荷載集度q，便可得到最右端節(jié)點(diǎn)的位移和彎矩。該程序借用matlab軟件中大量的矩陣賦值和矩陣運(yùn)算語句，通過對桿件單元和節(jié)點(diǎn)的數(shù)據(jù)化處理，實(shí)現(xiàn)計(jì)算機(jī)對結(jié)構(gòu)的計(jì)算，改程序經(jīng)過簡單修改后即可按照輸入指令，計(jì)算n列n跨任意位置加斜桿的鋼架體系，具有一定的普遍性。L=input(39。層高h(yuǎn) =\n39。)。EIb=input(39。梁的抗壓EAb=\n39。)。q=input(39。e=zeros(180,180)。S3=6*EIc/(h^2)。0 S2 S3。0 0 2*S4 ]。 a(30+3*jj2:30+3*jj,3*jj2:3*jj)=[S1 0 0。0 S2 S3。0 S3 S5]。%K22% enda(150+3*jj2:150+3*jj,150+3*jj2:150+3*jj)=[S1 0 0。for ii=1:3:178T(ii:ii+2,ii:ii+2)=[0 1 0。*a*T。S4=4*EIb/L。0 S3 S5]。 e(30*ii+4:30*ii+6,30*ii+4:30*ii+6)=[2*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 S2 S3。S1=EAo/d。S5=2*EIo/d。x(40+30*ii:40+30*ii+2,40+30*ii:40+30*ii+2)=[S1 0 0。0 S2 S3。0 S3 S4]。0 0 1]。 F0=zeros(180,1)。z(1,1)=*q*h。F=F0+F1。F(1:30,1)=0。for ii=1:5M1(ii,1)=[0 6*EIb/(L^2) 2*EIb/L 0 6*EIb/(L^2) 4*EIb/L ]*x1(25+30*ii:30+30*ii,1)。x2(4:6,1)=x1(58+30*ii:58+30*ii+2,1)。0 0 1]。x3=zeros(6,1)。for jj=1:3:6T(jj:jj+2,jj:jj+2)=[0 1 0。M3(ii+1,1)=[0 6*EIc/(h^2) 4*EIc/h 0 6*EIc/(h^2) 2*EIc/h]*x3。3。Q(1:5,3)=M2。)。從左至右依次為，節(jié)點(diǎn)層數(shù),節(jié)點(diǎn)側(cè)向位移39。3 x1(118
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1