MATLAB是一款功能強大的科學計算軟件，能夠完成很多任務，具有使用方便，輸入簡潔，上手快等優點。本文提出使用MATLAB進行梁單元模型的創建、計算和后處理。編寫的程序可用于結構優化，由于無需借助商業有限元軟件，計算高效，且無需擔心版權問題。感興趣的同學可以下載代碼，未經允許不得用于商業用途。

1. 有限元模型的建模

       有限元模型可包括五大部分：節點、單元（由節點組成）、材料、屬性（主要是橫截面，如截面尺寸、形狀）、邊界條件（約束和加載力）。定義了單元類型和材料及屬性，即可計算每個單元的剛度矩陣ki，并最終組裝所有單元的剛度矩陣得到整體剛度矩陣K。

2. 有限元模型的計算

       有限元模型的計算其實就是求方程組的過程，對于靜力問題有：F=Kx，F為外力矢量，K為整體剛度矩陣，x為節點位移。有了位移就可以計算其它量，如應變能e=0.5x'*K*x， 應力σ=DBx等

3. 有限元模型的后處理

       本文使用的代碼僅對節點的合位移進行演示，感興趣的同學可以改動代碼，做出自己想要展示的效果圖。