有限元鋼架結(jié)構(gòu)分析-手算matlabansys模擬-資料下載頁
2025-06-29 02:31本頁面
【正文】 l^2。 0 6*e*i/l^2 2*e*i/l 0 6*e*i/l^2 4*e*i/l]。t=[ cos(fai), sin(fai), 0。 sin(fai), cos(fai), 0。 0, 0, 1]。%坐標(biāo)變換矩陣T=blkdiag(t,t)。%總體坐標(biāo)系下的單元剛度矩陣K=T39。*k*T。%帶入每個單元的數(shù)，生成單元剛度矩陣kk，其每一頁對應(yīng)相應(yīng)頁數(shù)的單元的剛度矩陣for j=1:dy。 e=E。 i=I。 l=LL(j)。 a=A。 fai=FAI(j)。 kk(:,:,j)=eval(K)。end%生成總體剛度矩陣KK%采用元胞數(shù)組的方式對各項進行保存%生成空元胞數(shù)組，元胞數(shù)組的行列大小與節(jié)點數(shù)相同for j=1:jd。 for jj=1:jd。 ling1{j,jj}=zeros(3)。 endendling2=ling1。%將對單元剛度矩陣部分分成4分加入元胞數(shù)組中for j=1:dy。 kk1=kk(1:3,1:3,j)。 kk2=kk(1:3,4:6,j)。 kk3=kk(4:6,1:3,j)。 kk4=kk(4:6,4:6,j)。 ling2{dy_jd(j,2),dy_jd(j,2)}=kk1+ling2{dy_jd(j,2),dy_jd(j,2)}。 ling2{dy_jd(j,2),dy_jd(j,3)}=kk2+ling2{dy_jd(j,2),dy_jd(j,3)}。 ling2{dy_jd(j,3),dy_jd(j,2)}=kk3+ling2{dy_jd(j,3),dy_jd(j,2)}。 ling2{dy_jd(j,3),dy_jd(j,3)}=kk4+ling2{dy_jd(j,3),dy_jd(j,3)}。end%將元胞數(shù)組進行拼接，形成總體剛度矩陣for j=1:jd。 ling3(:,:,j)=cat(2,ling2{j,:})。endKK=ling3(:,:,1)。for j=2:jd。 KK=[KK。ling3(:,:,j)]。end%消去有已知位移的行與列b=KK。b(:,wys(:,3))=[]。b(wys(:,3),:)=[]。kjiejuzhen=inv(b)。%提取對應(yīng)外力lyss=lys。for j=1:size(wys,1)。 for jj=1:size(lys,1)。 if lyss(jj,4)==wys(j,3)。 lyss(jj,:)=0。 end if jj==size(lyss,1)。 break end endendlyss(all(lyss==0,2),:)=[]。%求解weiyijie=[作用值 作用節(jié)點 作用類型（轉(zhuǎn)角為3 x方向為1 y方向為2） 序列]weiyijie=kjiejuzhen*lyss(:,1)。weiyijie(:,1)=weiyijie。weiyijie(:,2)=lyss(:,2)。weiyijie(:,3)=lyss(:,3)。weiyijie(:,4)=lyss(:,4)。%計算不計作用在約束方向上時的支反力lysjiee=[作用值 作用節(jié)點 作用類型（轉(zhuǎn)角為3 x方向為1 y方向為2） 序列]lysjie(:,1)=KK(wys(:,3),lyss(:,4))*weiyijie(:,1)。lysjie(:,2:4)=wys(:,1:3)。%將作用在約束方向上時的支反力加在上面的求解結(jié)果上for j=1:size(lysjie,1) for jj=1:size(lys,1)。 if lysjie(j,4)==lys(jj,4)。 lysjie(j,1)=lysjie(j,1)lys(jj,1)。 end endend%答案weiyijie lysjie
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1