鋼架及matlab算法-免費(fèi)閱讀
【正文】 3 x1(118,1)。)。3。for jj=1:3:6T(jj:jj+2,jj:jj+2)=[0 1 0。0 0 1]。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)。F=F0+F1。 F0=zeros(180,1)。0 S3 S4]。x(40+30*ii:40+30*ii+2,40+30*ii:40+30*ii+2)=[S1 0 0。S1=EAo/d。 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。0 S3 S5]。*a*T。%K22% enda(150+3*jj2:150+3*jj,150+3*jj2:150+3*jj)=[S1 0 0。0 S2 S3。0 0 2*S4 ]。S3=6*EIc/(h^2)。q=input(39。梁的抗壓EAb=\n39。)。L=input(39。 二． 說明程序本程序使用matlab軟件通過矩陣位移法對(duì)題目進(jìn)行編程設(shè)計(jì)，按照提示輸入柱子的抗彎剛度EIc；柱子的抗壓剛度EAc；梁的抗彎剛度EIb；梁的抗壓剛度 EAb；斜桿的抗彎剛度EIo；斜桿的抗壓剛度EAo；側(cè)向均布荷載集度q，便可得到最右端節(jié)點(diǎn)的位移和彎矩。（5） 坐標(biāo)轉(zhuǎn)換對(duì)豎向坐標(biāo)系和斜桿體系進(jìn)行轉(zhuǎn)置，其坐標(biāo)轉(zhuǎn)換陣為其中。（2） 單元分析所有單元均為平面彎曲式自由單元如圖13所示。（9） 聯(lián)立方程求解通過方程求出位移矩陣。 三． 算法流程開始定義變量EIc,EAc,EIb……輸入變量EIc,EAc,EIb…… 豎向分體系剛度矩陣賦值水平分體系剛度矩陣賦值斜向分體系剛度矩陣賦值力的矩陣賦值求解方程邊界條件處理單元后處理輸出M, 結(jié)束四． 源代碼disp(39。)。梁的抗彎剛度EIb=\n39。EAo=input(39。x=zeros(180,180)。0 S3 S4]。0 S2 S3。%K21% a(30*(ii+1)+3*jj2:30*(ii+1)+3*jj,30*(ii+1)+3*jj2:30*(ii+1)+3*jj)=[2*S1 0 0。1 0 0。S5=2*EIb/L。0 2*S2 0。 e(30*ii+3*(jj+1)2:30*ii+3*(jj+1),30*ii+3*jj2:30*ii+3*jj)=[S1 0 0。0 S3 S4]。for ii=0:1:3x(7+30*ii:7+30*ii+2,40+30*ii:40+30*ii+2)=[S1 0 0。0 S3 S5]。endK3=R39。z(151,1)=*q*h。x1=K\F。 T=zeros(6,6)。for ii=
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1