結構分析理論早在電子計算機問世之前即已廣泛採用矩陣之形式求解.docVIP

結構工程 林聰悟 結構分析理論早在電子計算機問世之前即已廣泛採用矩陣之形式求解線性聯立方程式。例如以虛功原理、最小位能原理、虛補功原理等方法求解靜不定結構；求解靜定結構則最常採用力法以各種簡單之自由體直接求算單一力素或以二至三階之聯立方程式求算二至三個力素。 靜不定結構之贅餘力素之數目等於結構之靜不定度數，求解贅餘力素必須求解與靜不定度數同階之線性聯立方程式。在電腦問世之前，要求得10階之線性聯立方程式之解至少得花費一星期的時間。且隨著未知數的數目增加，其求解時間與未知數之數目成三次方的比例遞增。因此早期結構分析主要用於分析桁架或跨度數目不多的連續梁，或僅限於間數或層數不多之剛架結構。 結構分析在1929年哈氏(Hardy Cross)提出彎矩分配法之後進入一個斬新而實用的階段。此時層數或間數較多的剛架結構得以有效的加以分析。其實彎矩分配法屬於一種聯立方程式之反覆試算法。由於其收歛速度快，運算量亦相當少，可說是當時結構分析之不二法寶。 電子計算機於1950年代問世後，彎矩分配法仍然可以利用電腦的快速運算能力而發展出有關的分析程式。由於結構分析需要求解數量龐大的聯立方程式，而電腦的程式又很容易處理聯立方程式之矩陣運算，因此從事結構分析的學者專家們很快就意識到電腦的超強功能，並積極開發程式加以充分利用。因此電腦在工程上的應用發展最快的可能為結構分析。 結構分析由於所處理的對象相當複雜，未知數目又多，為了尋找有效的方法，因此研究發展出的各種解題方法可說是五花八門。不過隨著電腦問世，解題方法上則傾向於二大矩陣結構分析法：勁度矩陣法與柔度矩陣法。其中尤以勁度矩陣法之矩陣較易獲得較小之帶寬而較廣泛用於電腦分析程式中。 由於結構分析在電腦應用的快速發展而導致結構分析學者發展出有限元素分析法。因此而將傳統結構分析所處理的梁（柱）元素，一舉擴展應用於板元素、立體元素、甚至應用於各種多維元素。 有限元素法的提出確實對應用數學造成很大的衝擊。以應用數學的傳統觀點，對於像板或三維結構體等所謂聯體的分析，亦即彈性力學或聯體力學的數學解法，通常須求解偏微分方程式。而其所能求得的數學解析解則通常僅限於簡單的分析域及邊界條件。解題技巧多數引用級數解做近似分析。其所能求解的問題亦同樣受限於簡單的問題及邊界條件。在數值分析上之傳統做法為將偏微分方程式改以差分方程式近似。因此求解上較受限於規則的網點。 有限元素法自提出後，雖然在結構分析上已經廣泛而成功地用於各種彈性力學與聯體力學的分析領域，但對於應用數學家而言，仍多存有懷疑而排斥的態度。這其間，有限元素法亦以較為嚴謹的數學理論探討（雖然當初的提出也許謹憑直覺），証明有限元素法實為數學上泛函分析的一種推廣應用。傳統的基底函數常定義於整個分析域，因此會受限於簡單的分析域；而有限元素法則承傳結構分析法採用以一個桿件為局部範圍的基底函數，因此使得基底函數可以用比較簡單的形狀（基底）函數，如一次、二次或高次的多項式，而使得數值化的處理較為簡單。同時各元素可以較為隨意的分割而不若差分法須用較規則化的網格。且由於形狀函數的局部化，亦使得所獲致的系統矩陣之帶寬較小。由於有限元素法的諸多優點均為傳統差分法所無法比擬，且又有泛函分析的嚴謹數學理論可以佐證，因此在經過數十年的延遲之後，有限元素法即廣受應用數學界之學者所接受而推廣應用於各種領域上並已成為解偏微分方程的數值方法之主要利器。 電腦計算速度即使遠非人腦所能比擬，但以電腦分析大型結構則永遠有記憶容量不足與計算時間太久之問題。早期記憶體在百萬位元(Maga Bytes)以下當然存在記憶容量不足的困境。今日即使記憶容量已達20億位元(2Giga Bytes)，仍然有記憶容量不足的問題。而隨著記憶容量增大，所能分析的問題也隨之加大。電腦運算速度的增加亦往往趕不上需求運算量的增加，因此計算時間也可能長到不能接受的程度。 以目前電腦的記憶容量(2GB)與運算速度(200MFLOPS)而言，分析一般的平面有限元素通常不會有太大的問題。以最不利的正方形規則網點計算：設每邊網點數為，每一網點之未知變數之數目為，則矩陣階數為，半帶寬，所須記憶容量以倍精度實數計算為。設，，得。足可應付一般的分析需求。其矩陣分解所需運算量（以一次乘算與和算為運算單位）為，運算時間約為560秒。 但如欲分析三維的有限元素問題，仍以最不利的立方形規則網點計算：設每邊網點數為，每一網點之未知變數之數目為，則矩陣階數為，半帶寬，所須記憶容量為。設，，僅僅得。實不足以應付實際應用的分析需求。且即使以該等大小之矩陣計算，其矩陣分解所需運算量為，運算時間約為1860秒。值得注意的是每邊網點數如加倍，記憶容量增加為32倍，運算量增加為128倍。 如要突破上述限制，必須以結構分段方式處理或以機外儲存方式處理。以不失運算效率的情形