有限元鋼架結(jié)構(gòu)分析-手算matlabansys模擬(參考版)
2025-07-02 02:31本頁面
【正文】 end endend%答案weiyijie lysjie。 if lysjie(j,4)==lys(jj,4)。lysjie(:,2:4)=wys(:,1:3)。weiyijie(:,4)=lyss(:,4)。weiyijie(:,2)=lyss(:,2)。%求解weiyijie=[作用值 作用節(jié)點 作用類型（轉(zhuǎn)角為3 x方向為1 y方向為2） 序列]weiyijie=kjiejuzhen*lyss(:,1)。 end if jj==size(lyss,1)。 if lyss(jj,4)==wys(j,3)。for j=1:size(wys,1)。kjiejuzhen=inv(b)。b(:,wys(:,3))=[]。ling3(:,:,j)]。for j=2:jd。 ling3(:,:,j)=cat(2,ling2{j,:})。 ling2{dy_jd(j,3),dy_jd(j,3)}=kk4+ling2{dy_jd(j,3),dy_jd(j,3)}。 ling2{dy_jd(j,2),dy_jd(j,3)}=kk2+ling2{dy_jd(j,2),dy_jd(j,3)}。 kk4=kk(4:6,4:6,j)。 kk2=kk(1:3,4:6,j)。%將對單元剛度矩陣部分分成4分加入元胞數(shù)組中for j=1:dy。 ling1{j,jj}=zeros(3)。end%生成總體剛度矩陣KK%采用元胞數(shù)組的方式對各項進(jìn)行保存%生成空元胞數(shù)組，元胞數(shù)組的行列大小與節(jié)點數(shù)相同for j=1:jd。 fai=FAI(j)。 l=LL(j)。 e=E。*k*T。%坐標(biāo)變換矩陣T=blkdiag(t,t)。 sin(fai), cos(fai), 0。 0 6*e*i/l^2 2*e*i/l 0 6*e*i/l^2 4*e*i/l]。 e*a/l 0 0 e*a/l 0 0。 0 12*e*i/l^3 6*e*i/l^2 0 12*e*i/l^3 6*e*i/l^2。wys=sortrows(wys,3)。lys(:,4)=(lys(:,2)1)*3+lys(:,3)。3 3]。3 1。1 2。P*L/3 3 3]。67/750*P*L 2 3。81/125*P 2 1。% lys=[44/125*P 1 1。% q=P/L。2 2 3]。I=20*A。jd=3。% A=。% E=3*10^10。)。%輸入結(jié)構(gòu)約束[作用節(jié)點 作用類型]（轉(zhuǎn)角為3 x方向為1 y方向為2）wys=input(39。力與轉(zhuǎn)矩約束矩陣：39。)。%輸入對應(yīng)關(guān)系時，小節(jié)點放前面[單元 節(jié)點1 節(jié)點2]dy_jd=input(39。請輸入單元相對旋轉(zhuǎn)角度：39。)。A=input(39。請輸入單元長度L：39。)。I=input(39。請輸入楊氏模量E：39。)。dy=input(39。請輸入節(jié)點數(shù)：39。過程性指的是在編程要有一個完整的流程圖，將復(fù)雜問題轉(zhuǎn)化成若干簡單的小問題，針對小問題進(jìn)行求解；結(jié)構(gòu)性指的是在編程中要利用條件if，循環(huán)for等結(jié)構(gòu)簡化程序，同時也要對每部分程序的輸入、輸出以及執(zhí)行作用有一個明確的認(rèn)知；最后根據(jù)使用對象，對輸入、輸出，進(jìn)行調(diào)整，也可以用GUI編制界面方便使用。而想要提高編程能力，首先要有一定的編程思想與數(shù)學(xué)建模能力。再獲得總體剛度矩陣后，只要找到對應(yīng)的位移與外力的參數(shù)行，再進(jìn)行乘除計算就相對簡單了。再通過模板獲得每個單元的矩陣。深化我對有限元計算流程的理解，提升了我的matlab的編程水平。八、編程：見附錄一第四章 有限元求解一、預(yù)處理選擇單元類型：ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete… →Add… →beam：2D elastic 3→OK (返回到Element Types 窗口) →Close 選擇單元類型定義材料參數(shù)：ANSYS Main Menu: Preprocessor →Ma
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1