鋼架及matlab算法-全文預(yù)覽
【正文】 0鋼架側(cè)移量（最右側(cè)）從左至右依次為，節(jié)點(diǎn)層數(shù),節(jié)點(diǎn)側(cè)向位移 。4 x1(148,1)。)。disp(Q)disp(39。Q(1:5,4)=M3。4。endQ=zeros(5,4)。1 0 0。for ii=0:1:3x3(1:3,1)=x1(58+30*ii:58+30*ii+2,1)。endx2=T*x2。 T=zeros(6,6)。endM2=zeros(5,1)。x1=K\F。 K(1:30,:)=0。z(151,1)=*q*h。z=zeros(180,1)。endK3=R39。endR=zeros(180,180)。0 S3 S5]。0 S2 S3。for ii=0:1:3x(7+30*ii:7+30*ii+2,40+30*ii:40+30*ii+2)=[S1 0 0。S2=12*EIo/(d^3)。0 S3 S4]。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 S5]。0 2*S2 0。 e(30*ii+1:30*ii+3,30*ii+1:30*ii+3)=[S1 0 0。S5=2*EIb/L。 S1=EAb/L。1 0 0。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。0 S3 S5]。0 S2 S3。 a(3*jj2:3*jj,30+3*jj2:30+3*jj)=[S1 0 0。0 S3 S4]。S4=4*EIc/h。x=zeros(180,180)。水平均布荷載集度q =\n39。EAo=input(39。)。梁的抗彎剛度EIb=\n39。EAc=input(39。)?？缍萀=\n39。 三． 算法流程開始定義變量EIc,EAc,EIb……輸入變量EIc,EAc,EIb…… 豎向分體系剛度矩陣賦值水平分體系剛度矩陣賦值斜向分體系剛度矩陣賦值力的矩陣賦值求解方程邊界條件處理單元后處理輸出M, 結(jié)束四． 源代碼disp(39。在建立剛度矩陣，引入ii和jj作為循環(huán)變量，根據(jù)剛度矩陣特點(diǎn)對(duì)節(jié)點(diǎn)逐列逐跨進(jìn)行矩陣塊的循環(huán)賦值。（9） 聯(lián)立方程求解通過方程求出位移矩陣。（6） 建立整體剛度方程將三個(gè)剛度矩陣相加，得到整體剛度方程。（2） 單元分析所有單元均為平面彎曲式自由單元如圖13所示。（1） 位移編
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1