有限元鋼架結(jié)構(gòu)分析-手算matlabansys模擬-閱讀頁
2025-07-14 02:31本頁面
【正文】 的本質(zhì)是將每個單元的矩陣按照矩陣所對應(yīng)節(jié)點的信息進行分塊，并累加到總體矩陣中，由于矩陣實現(xiàn)較為麻煩，所以我想利用元胞數(shù)組來存儲節(jié)點信息，并成功實現(xiàn)了所需功能，又利用cat拼接指令，將元胞數(shù)組展開成一個大型矩陣完成拼接。 我認為在學習過程中，提高編程能力是很重要的一個方面，有著良好的編程能力，可以讓很多工程問題得以用計算機解決，更容易獲得結(jié)果。 我認為編程思想指的是，要對所做的事情有一個過程性與結(jié)構(gòu)性的認識，并根據(jù)使用對象進行相關(guān)的調(diào)整。 以上是我對這次大作業(yè)的心得體會，希望老師輔導(dǎo)校正第七章 附錄一、matlab程序clcclearformat pactformat shortGjd=input(39。)。請輸入單元數(shù)：39。E=input(39。)。請輸入慣性矩I：39。L=input(39。)。請輸入單元截面積：39。FAI=input(39。)。請輸入單元與節(jié)點對應(yīng)關(guān)系：39。%輸入力與扭矩約束[值 作用節(jié)點 作用類型]（轉(zhuǎn)矩為3 x方向為1 y方向為2）lys=input(39。)。結(jié)構(gòu)約束矩陣：39。%原始數(shù)據(jù)% L=1。% P=1000。% dy=2。LL=[L 2*L]。% dy_jd=[1 1 2。% FAI=[pi/2 0]。M=P*L/10。12*P*L/125 1 3。P 2 2。P 3 2。% wys=[1 1。1 3。3 2。%對力約束與位移約束式子分別進行編號處理wys(:,3)=(wys(:,1)1)*3+wys(:,2)。%對力約束與位移約束式子進行排序lys=sortrows(lys,4)。%單元剛度矩陣syms fai e a i l realk=[e*a/l 0 0 e*a/l 0 0。 0 6*e*i/l^2 4*e*i/l 0 6*e*i/l^2 2*e*i/l。 0 12*e*i/l^3 6*e*i/l^2 0 12*e*i/l^3 6*e*i/l^2。t=[ cos(fai), sin(fai), 0。 0, 0, 1]。%總體坐標系下的單元剛度矩陣K=T39。%帶入每個單元的數(shù)，生成單元剛度矩陣kk，其每一頁對應(yīng)相應(yīng)頁數(shù)的單元的剛度矩陣for j=1:dy。 i=I。 a=A。 kk(:,:,j)=eval(K)。 for jj=1:jd。 endendling2=ling1。 kk1=kk(1:3,1:3,j)。 kk3=kk(4:6,1:3,j)。 ling2{dy_jd(j,2),dy_jd(j,2)}=kk1+ling2{dy_jd(j,2),dy_jd(j,2)}。 ling2{dy_jd(j,3),dy_jd(j,2)}=kk3+ling2{dy_jd(j,3),dy_jd(j,2)}。end%將元胞數(shù)組進行拼接，形成總體剛度矩陣for j=1:jd。endKK=ling3(:,:,1)。 KK=[KK。end%消去有已知位移的行與列b=KK。b(wys(:,3),:)=[]。%提取對應(yīng)外力lyss=lys。 for jj=1:size(lys,1)。 lyss(jj,:)=0。 break end endendlyss(all(lyss==0,2),:)=[]。weiyijie(:,1)=weiyijie。weiyijie(:,3)=lyss(:,3)。%計算不計作用在約束方向上時的支反力lysjiee=[作用值 作用節(jié)點 作用類型（轉(zhuǎn)角為3 x方向為1 y方向為2） 序列]lysjie(:,1)=KK(wys(:,3),lyss(:,4))*weiyijie(:,1)。%將作用在約束方向上時的支反力加在上面的求解結(jié)果上for j=1:size(lysjie,1) for jj=1:size(lys,1)。 lysjie(j,1)=lysjie(j,1)lys(jj,1)
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1