平面三角形單元常應(yīng)變單元matlab程序的編制-資料下載頁
2025-07-07 14:03本頁面
【正文】 SDISP LNODS%ELEDISP(1:6)=0。 %當(dāng)前單元的結(jié)點位移向量l 說明：ELEDIS為單元的結(jié)點位移% for i=1:NELEM for j=1:3 ELEDISP(j*21:j*2)=ASDISP(LNODS(i,j)*21:LNODS(i,j)*2)。 end FORMSMATX(i) %% %調(diào)用子程求應(yīng)力矩陣 i STRESS=SMATX*ELEDISP39。 %求內(nèi)力 endl 說明：通過循環(huán)，依次從結(jié)構(gòu)位移列陣中對號，賦值給第i個單元的單元位移向量ELEDISP。% 程序應(yīng)用舉例【例1】利用三角形三節(jié)點常應(yīng)變單元程序計算圖2所示的結(jié)構(gòu)。設(shè)彈性模量，泊松比為0，厚度為1。 %：6 4 1 6 13 1 25 2 43 2 5 6 3 5 圖2 1 0 11 1 02 1 04 1 15 0 16 0 1%說明：第一行：讀入程序控制信息NPION=fscanf(FP1,39。%d39。,1) %結(jié)點個數(shù)（結(jié)點編碼總數(shù)）NELEM=fscanf(FP1,39。%d39。,1) %單元個數(shù)（單元編碼總數(shù)）NFORCE=fscanf(FP1,39。%d39。,1) ％結(jié)點荷載個數(shù)NVFIX=fscanf(FP1,39。%d39。,1) ％受約束邊界點數(shù)YOUNG=fscanf(FP1,39。%e39。,1) ％彈性模量POISS=fscanf(FP1,39。%f39。,1) ％泊松比 THICK=fscanf(FP1,39。%d39。,1) %厚度第二、三、四、五行：讀入單元連接信息：LNODS=fscanf(FP1,39。%d39。,[3,NELEM])39。 %單元定義數(shù)組，單元結(jié)點號，逆時針輸入第六、七、八、九、十、十一行：讀入結(jié)點坐標(biāo)COORD=fscanf(FP1,39。%f39。,[2，NPION])39。 %結(jié)點坐標(biāo)值，第1列為x坐標(biāo)值，第2列為y坐標(biāo)值第十一行：讀入結(jié)點荷載信息FORCE=fscanf(FP1,39。%f39。,[3,NFORCE])39。 %結(jié)點號，X向結(jié)點荷載數(shù)值，Y向結(jié)點荷載數(shù)值（與坐標(biāo)軸方向一致為正）第十二、十六行：讀入零位移信息FIXED=fscanf(FP1,39。%d39。,[3,inf])39。 %結(jié)點號，X向約束，Y向約束% 程序的改進要點上述三角形三節(jié)點程序反映了有限元的基本思路，可以計算簡單的平面應(yīng)力或平面應(yīng)變問題。在熟練掌握了程序的編制與使用后，可在以下幾方面對程序進行改進，以加深對矩陣位移法及MATLAB語言編程的理解：本程序的彈性模量僅能輸入一個數(shù)值，意味著程序僅能計算由同種材料構(gòu)成的結(jié)構(gòu)?？紤]如何改進使程序可計算由不同材料構(gòu)成的組合結(jié)構(gòu)。本程序僅能計算結(jié)點集中荷載類型，考慮如何編制體積力、表面力、跨中集中力等類型的程序?？紤]如何編制有支座位移的程序。本程序最后的結(jié)果沒有生成輸出文件，可以編制輸出文件。本程序計算的應(yīng)力沒有進行結(jié)果處理，可以編制最后結(jié)果處理的程序?？梢栽诖顺绦蚧A(chǔ)上編制三角形六節(jié)點、四邊形四節(jié)點等程序。綜上所述，本章的三角形三節(jié)點常應(yīng)變程序體現(xiàn)了如何將有限元法的計算方法和過程用MATLAB程序語言表達(dá)出來，重點放在程序架構(gòu)和流程的建立以及算法實現(xiàn)方面，主要依賴手工操作－手工輸入初始數(shù)據(jù)（前處理）、手工繪制計算結(jié)果（后處理），目的是使學(xué)生能夠清晰、明確地掌握有限元法的基本理論和概念，熟練掌握應(yīng)用MATLAB語言編制、修改和調(diào)試簡單程序的能力。13
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1