楊永斌臺灣大學：土木工程學研究所簡秀伊Chien, Hsiu-YiHsiu-YiChien2007-11-252018-07-092007-11-252018-07-092006http://ntur.lib.ntu.edu.tw//handle/246246/50606本文旨在探討以迭代法改善有限元素法求解之精度，以曲梁為例，分成線性與幾何非線性分析兩個部分加以研究。 一般以有限元素法分析平面曲梁線性問題時，將曲梁切割成許多元素來表示。然而這些元素因為受到諧和條件及形狀函數選擇的限制，無法完全滿足曲梁的各項條件。因此在分析時，常會造成很大的誤差。本文將引用由前人所提出的一套迭代流程，來改善分析結果之精準度。由前人之研究可知，在迭代流程中，預測子(predictor)只影響迭代的次數或收斂的速度，而校正子(corrector)卻會影響計算結果的精準度。為驗證此一道理，本文將採用品質並不完美的曲梁元素之幾何勁度矩陣作為預測子，並以準確的校正子，即準確之力與位移關係式，透過迭代流程以消除分析結果之誤差。 另一方面，針對幾何非線性分析部分，前人已經運用剛體運動法則推導出直梁元素之幾何勁度矩陣，本文將進一步延伸應用剛體法則於平面曲梁元素之幾何勁度矩陣的推導，並以外顯式表示之。對於剛體而言，其幾何勁度矩陣可視為模擬元素在受力狀態下進行剛體位移的行為，與材料係數及斷面特性無關聯，故當元素兩端具有相同剛體模態時，不論它是直梁或曲梁元素，都具有相同的幾何勁度矩陣。根據此一原理，我們可以將直梁元素之幾何勁度矩陣，經座標系統轉換後，而得到曲梁元素之幾何勁度矩陣。此一矩陣的可靠性將在實例中加以驗證。In this study, the concept of iterations will be adopted to improve the accuary of finite element solutions. We shall take the planar curved beam as an example. Both linear and geometrically nonlinear analyses will be discussed. In analyzing the linear behavior of planer curved beams by the finite element method, the beam is usually discretized into a number of elements. However, owing to the constraints of compatibility conditions and selection of shape functions, some errors can still occur for the linear problems. In this study, theiterative procedure proposed by previous researchers will be employed to improve the accuracy of the solution obtained. According to previous researches, in an iterative procedure, the predictor affects only the number of iterations or the speed of convergence, while the corrector determines solely the accuracy of the iterative solution. To verify such a concept, stiffness matrices with some defects will be used in this study as the predictor, while accurate force-displacement relations will be used as the corrector. It will be demonstrated thatthrough the iterative procedure, the accuracy of the solution can be significantly improved.. For the nonlinear part, it is realized that the rigid body rule was successfully applied to derivation of the geometric stiffness matrix for the planer straight beam element. Such a procedure will be followed in this study to derive the geometric stiffness matrix for the planer curved beam problem, which will be presented in explicit form. As for a rigid body, the geometric stiffness matrix can be regard as the ability of the element with initial forces in undergoing the rigid motion. For this reason, the geometric stiffness matrix is only related to the shape of the element, but not the properties of the materials or cross sections. Thus, the geometric stiffness matrices for the straight and curved beams will be the same, as long as they have identical nodal degrees of freedom. It follows that the geometric stiffness matrix for a curved beam can derived from that for the corresponding straight beam merely by transformation of the coordinates from the Cartesian system to the curvilinear system. The capability of such an element will be verified in the numerical studies.致謝 / 一 摘要 / 二 目錄 / 六 表目錄 / 一O 圖目錄 / 一一 1 緒論/1 1.1 研究動機與目的 / 1 1.2 論文架構 / 5 2 文獻探討/7 2.1 前言 / 7 2.2 以迭代法修正梁元素之缺陷 / 8 2.3 迭代流程 / 14 2.4 直梁之幾何勁度矩陣 / 16 3 迭代法應用於曲梁線性分析/23 3.1　前言 / 23 3.2　迭代法 / 24 3.3 預測子之推導 / 25 3.3.1 形狀函數 / 26 3.3.2 預測子 / 28 3.4 校正子之推導 / 30 3.4.1 形狀函數 / 31 3.4.2 校正子 / 34 3.5 迭代流程應用於平面曲梁之線性分析 / 35 3.6 數值實例分析 / 38 3.6.1 單對稱斷面之1/4圓曲梁 / 38 3.6.2 雙對稱斷面之1/4圓曲梁 / 40 3.6.3 二端點均為鉸接的圓拱 / 41 3.6.4 二端點均為固定端的圓拱 / 42 3.7 小結 / 43 4 平面曲梁幾何勁度矩陣之簡易推導法/63 4.1　增量平衡方程式 / 63 4.2 非線性元素之剛體測試 / 65 4.2.1 初始受力元素之剛體運動法則 / 65 4.2.2 剛體運動測試 / 68 4.3 自然變形和剛體位移 / 72 4.4 以剛體法則推導幾何勁度矩陣 / 73 4.4.1 桁架元素 / 74 4.4.2 平面直梁元素 / 76 4.4.3 平面曲梁元素 / 78 4.5 平面曲梁之彈性勁度矩陣 / 82 4.6 非線性曲梁元素之剛體運動測試 / 83 4.6.1 剛體移動 之測試 / 84 4.6.2 剛體移動 之測試 / 85 4.6.3 剛體轉動 之測試 / 86 4.7 線性化挫屈 / 88 4.7.1 線性化挫屈理論 / 88 4.7.2 實例分析 / 91 4.8 小結 / 93 5 有限元素之增量非線性分析/105 5.1　簡介 / 105 5.1.1 模擬方法 / 104 5.2 增量非線性分析程序 / 106 5.2.1 預測階段－預測矩陣 / 107 5.2.2 校正階段－元素內力計算 / 107 5.3 非線性增量－迭代方法 / 110 5.3.1 廣義位移控制法 / 111 5.3.2 非線性增量-迭代分析流程 / 116 5.4 實例分析 / 120 5.4.1 受中央垂直載重作用之深拱 / 121 5.4.2 受中央垂直載重作用之圓拱結構 / 121 5.4.3 受中央荷重作用之淺圓拱 / 122 5.4.4 受中央荷重作用之圓拱 / 123 5.5 小結 / 124 6 結論與展望/133 6.1　結論 / 133 6.2 未來展望 / 134 參考文獻/1371136651 bytesapplication/pdfen-US迭代法幾何非線性幾何勁度矩陣Iterative MethodCurved BeamGeometric Stiffness Matrixthesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/50606/1/ntu-95-R93521218-1.pdf