單元剛度矩陣(等參元)matlab編程-資料下載頁
2025-07-26 11:57本頁面
【正文】 l 算例2：如圖10所示曲四邊形等參單元，已知8個結(jié)點整體坐標(biāo)系內(nèi)的坐標(biāo)為1(5,4)。 2(6,12)。 3(10,8)。 4(8,12)。 5(0,6)。 6(5,1)。 7(0,4)。 8(4,2)。彈性模量為200GPa，、。試求單元剛度矩陣。圖10 算例2曲四邊形等參單元解：根據(jù)如圖10所示曲四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc。clear。%*********輸入結(jié)點坐標(biāo)數(shù)組********xy8=[5,4。6,12。10,8。8,12。0,6。5,1。0,4。4,2]。 %****輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)****mat=[2e11,]。 K8 = ele_mat_quad8(xy8,mat)運行程序，得到單元剛度矩陣K8（GPa）如下： l 算例3：如圖11所示曲四邊形等參單元，已知8個結(jié)點整體坐標(biāo)系內(nèi)的坐標(biāo)為1(6,4)。 2(6,4)。 3(6,4)。 4(6,4)。 5(0,3)。 6(3,0)。 7(0,3)。 8(3,0)。彈性模量為200GPa，、。試求單元剛度矩陣。圖11 算例3曲四邊形等參單元解：根據(jù)如圖11所示曲四邊形等參單元及其幾何和材料參數(shù)，編制主程序如下：clc。clear。%*********輸入結(jié)點坐標(biāo)數(shù)組********xy8=[6,4。6,4。6,4。6,4。0,3。3,0。0,3。3,0]。 %****輸入材料參數(shù)矩陣(彈性模量，泊松比，壁厚)****mat=[2e11,]。 K8 = ele_mat_quad8(xy8,mat)運行程序，得到單元剛度矩陣K8（GPa）如下： 總結(jié)與備注這次有限元編程實驗給了我很多收獲，不僅鍛煉了Matlab編程的能力，還加強了Word排版的能力，最重要的是加深了對有限單元法的理解，加深了自己的力學(xué)功底。本次實驗的所有算例，不僅都通過我的程序得出結(jié)果，還經(jīng)過宋戎妝同學(xué)的程序驗證，故可以保證程序與算例的正確性。在這里衷心地感謝周博教授的指導(dǎo)，同時感謝幾位幫助過我的同學(xué)
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1