MATLAB在土木工程中的應(yīng)用

申請(qǐng)免費(fèi)試用、咨詢電話：400-8352-114

MATLAB是由美國MathWorks公司于 1984年發(fā)布的、面向科學(xué)計(jì)算、數(shù)據(jù)可視化以及交互程序設(shè)計(jì)的高級(jí)計(jì)算機(jī)語言。MATLAB軟件開發(fā)的初衷是方便地進(jìn)行矩陣運(yùn)算，如今的MATLAB已經(jīng)把功能延伸到了科學(xué)研究和工程應(yīng)用的諸多領(lǐng)域。在國外，MATLAB已經(jīng)成為數(shù)值分析、數(shù)理統(tǒng)計(jì)、系統(tǒng)識(shí)別、信號(hào)處理、動(dòng)態(tài)仿真等領(lǐng)域的基本工具。同傳統(tǒng)的土木工程計(jì)算機(jī)語言FORTRAN和C相比，MATLAB 更具有以下幾大優(yōu)勢(shì)：
a)功能強(qiáng)大。除了數(shù)值計(jì)算上的絕對(duì)優(yōu)勢(shì)外，MATLAB5．2以上版本中還開發(fā)了自己的符號(hào)計(jì)算功能。用戶只要擁有一個(gè)MAT ，就可以方便的處理諸如矩陣運(yùn)算，線性與非線性方程求解，特征值問題，統(tǒng)計(jì)及優(yōu)化問題了。
b)語言簡單。MATLAB允許用戶以數(shù)學(xué)形式的語言編寫程序，比FORTRAN和C更接近計(jì)算公式的書寫思維方式。它的操作和功能函數(shù)指令就是常用計(jì)算機(jī)和數(shù)學(xué)上一些簡單的英文單詞。
c)擴(kuò)充能力強(qiáng)、可開發(fā)性強(qiáng)。MATLAB本身的函數(shù)就是以源代碼形式出現(xiàn)，用戶可以仿照其寫法，創(chuàng)建自己的程序“庫”。另外，MATLAB可以方便的與FORTRAN 或C語言接口，充分利用現(xiàn)有資源。
d)編程易，效率高。MATLAB程序文件是純文本文件，用任何字處理文件都可以對(duì)它進(jìn)行編寫和修改。因此程序易調(diào)試，人機(jī)交互性強(qiáng)。另外，MATLAB自己也擁有視窗環(huán)境下的調(diào)試系統(tǒng)，程序調(diào)試簡單、方便。隨著科研水平的不斷提高，科研領(lǐng)域的不斷擴(kuò)展，多學(xué)科并行、交叉發(fā)展已成必然，認(rèn)識(shí)到MATLAB的強(qiáng)大功能，使得MATLAB在土木工程領(lǐng)域中得到充分應(yīng)用，達(dá)到利用MATLAB來快速實(shí)現(xiàn)科研構(gòu)想和提高工作效率的效果。
1 MATLAB在結(jié)構(gòu)分析中的應(yīng)用
土木工程中的結(jié)構(gòu)分析主要是指結(jié)構(gòu)在靜力和動(dòng)力荷載作用時(shí)結(jié)構(gòu)內(nèi)力和位移的求解問題。由于結(jié)構(gòu)的復(fù)雜多變和對(duì)求解精度的嚴(yán)格要求，采用有限元法(finite element method)程序是常用的分析手段。有限元程序中非線性方程組的求解和解的收斂問題，是困擾許多科研人員的基本問題之一，采用FORTRAN或C不但需要很多語句，而且需要研究者有較高的計(jì)算機(jī)理論和實(shí)際操作水平，程序的編制與調(diào)試又要占用近一半的研究時(shí)間。MATLAB的出現(xiàn)，可以在用簡單的語句完成基本算法程序后，方程組的求解利用MATLAB自身的命令實(shí)現(xiàn)，整個(gè)程序不但小巧、便于調(diào)試與操作，而且解的收斂問題易于解決。例如，經(jīng)典的非線性方程組Newton解法MATLAB求解程序僅有l(wèi)4行，而用C語言編寫的Newton法子函數(shù)程序接近30行，而且在每個(gè)迭代步當(dāng)Newton方程建立后，還要再調(diào)用其他子函數(shù)程序求解線性方程組，較MATLAB程序要繁瑣許多。再有，對(duì)圖1所示的為整體坐標(biāo)系中二維一次四結(jié)點(diǎn)單元，單元可以是任意形狀的四邊形。這種單元比較容易滿足結(jié)構(gòu)實(shí)際分析的需要，該單元在整體坐標(biāo)系下(x0y坐標(biāo)系)的形函數(shù)表達(dá)式非常復(fù)雜，用FORTRAN或C語言將其表達(dá)出來，是非常困難的。而用MATLAB語言就可以很容易的將其寫出。如二維一次四結(jié)點(diǎn)單元的MATLAB程序命令如下：