陳正興臺灣大學：土木工程學研究所柯永彥Ko, Yung-YenYung-YenKo2007-11-252018-07-092007-11-252018-07-092006http://ntur.lib.ntu.edu.tw//handle/246246/50405對稱葛勒金邊界元素法(SGBEM)的最大優點為：能建立具對稱形式之系統矩陣，同時由於加權剩餘法之應用，使得奇異積分之處理較容易，且在數值精確度與收斂行為上，較傳統點配置邊界元素法(CBEM)表現為佳，因此在近年來，被廣泛應用在各領域之邊界元素法分析。本研究將針對均質等向材料之三維靜力問題，利用無限域基本解(Kevin’s fundamental solution)，組成對偶邊界積分方程(DBIE)，建立一採四邊形線性元素之SGBEM分析模式。對於SGBEM之雙重邊界積分中，所遭遇之奇異積分問題，包含弱奇異、強奇異與超奇異積分，將利用座標轉換與積分域投影技巧，配合基本解核函數之對稱特性與基本性質處理之。經由範例分析驗證可知，本研究所建立之SGBEM模式，應用在具有良好支承條件之三維靜力問題，能得到相當正確的分析結果。 在數值分析方法中，結合有限元素法(FEM)與邊界元素法(BEM)共同使用，能兼採兩者之優點，有利於分析幾何配置與材料性質複雜的問題。當採用SGBEM與FEM連結使用時，經由適當的處理技巧，能夠得到具對稱形式之整體系統矩陣，而不破壞原本有限元素部分之系統對稱性。在本研究所採用之FE~BE連結模式中，將藉由SGBEM推導式中與自由項係數相關之部分，建立FE與BE間節點力與曳引力之轉換關係。 然而，SGBEM應用於分析力邊界值問題時，不論其屬內域或外域問題，所得位移解均會存在有一剛體運動項，而有解非唯一之現象。本研究將討論此現象之內涵與成因，並彙整一些能去除剛體運動影響的方法，藉由範例分析，討論這些方法運用在SGBEM與FE~BE連結時之適用性。結果顯示，對於外力條件之合力與合力矩為零的一般力邊界值問題，這些方法多能有效的去除剛體運動，得到唯一且正確的解；而在應用自由表面截切技巧所從事之半無限域問題分析上，由於其外力條件之合力或合力矩不為零，且因自由表面之截切造成誤差，因此雖然能夠得到唯一解，但各種方法所得之結果，仍然會存在一些不合理的現象。在所檢驗的方法中，以基於Fredholm理論所建立之修正邊界積分方程，所得之結果相對上較佳。As compared to the conventional collocation boundary element method (CBEM), the primary advantage of the symmetric Galerkin boundary element method (SGBEM) is the ability to produce a symmetric system matrix. Besides, through the application of the weighted residual method in SGBEM, it becomes easier to regularize the singular integrals and shows better numerical behaviors than the conventional CBEM. In this study, a SGBEM procedure to form linear quadrilateral elements for 3-D elastostatic problems is established. The formulation adopts the form of dual boundary integral equations (DBIEs) based on the Kevin’s fundamental solutions. In order to evaluate the singular double integrals in the SGBEM, including the types of weak, strong or hyper-singularity, techniques of coordinate transformation and mappings of integral domains are utilized in company with the basic properties of kernel functions. Through some numerical examples, the validity of this SGBEM procedure in application to well-supported 3-D elastostatic problems is verified. The coupling of FEM and BEM is a profitable result that exploits the advantages of each. When SGBEM is combined with FEM through appropriate techniques, a symmetric global system matrix can be obtained without ruining the symmetric virtue of the FE part. The FE~BE coupling strategy adopted in this study is to use the free-term components in the SGBEM, so that the equilibrium conditions between the nodal forces on the FE part and the nodal tractions on the BE part are reserved. Nevertheless, when SGBEM is applied to interior or exterior Neumann problems, the displacement solution cannot be uniquely determined because of the inevitable rigid body motion term involved. In this study, discussions on this phenomenon and methods that had been used to remove the rigid-body-motion terms are investigated for the problems modeled with the SGBEM and the FE~BE coupling formulations. For general equilibrated Neumann problems, the rigid body motions can be effectively removed by using these approaches. However, for half-space problems in which the free surface are modeled by limited number of elements, the solutions obtained are still not satisfactory because of the errors introduced from the truncation of the free surface. Among the methods investigated, the one using the modified boundary integral equations based on the Fredholm theory is relatively preferable.謝 誌 i 摘 要 ii Abstract iii 目 錄 iv 圖表目錄 vii 第一章 前言 1 1.1 研究動機與目的 1 1.2 研究方向與方法 2 第二章 文獻回顧 4 2.1 邊界積分方程與邊界元素法 4 2.2 對稱化邊界元素法 9 2.3 奇異積分之規則化 13 2.4 有限元素與邊界元素之連結使用 15 2.5 邊界元素法遭遇之秩降問題 17 第三章 線彈性材料三維靜力問題之對稱邊界元素法 24 3.1 對偶邊界積分方程：u-BIE與t-BIE 24 3.2 邊界方程之離散化與加權剩餘法之應用 25 3.3 對稱葛勒金邊界元素法(SGBEM) 28 3.4 域內點位移與應力以及邊界應力之計算 30 第四章 SGBEM中之奇異積分處理 36 4.1 SGBEM中遭遇之各種奇異積分情況 36 4.2 奇異積分規則化之座標轉換技巧概述 37 4.3 各種奇異積分情況之處理 38 4.3.1 完全重合情況(coincident case) 38 4.3.2 邊緣相鄰情況(edge-adjacent case) 41 4.3.3 頂點相鄰情況(vertex-adjacent case) 42 4.4 具良好支承條件之彈性靜力問題分析實例 43 4.4.1 立方體受單軸張力 43 4.4.2 Leon問題 44 第五章 SGBEM與FEM之連結使用 54 5.1 SGBEM與FEM之結合 54 5.2 具良好支承條件之彈性靜力問題分析實例 56 5.2.1 矩形體受單軸張力 56 5.2.2 薄版受彎矩作用 57 第六章 SGBEM之解非唯一問題與處理方法 64 6.1 SGBEM於力邊界值問題中之解非唯一現象 64 6.2 外加點支承 66 6.3 奇異值分解－特徵值分解 67 6.4 奇異系統之規則化 69 6.5 修正邊界積分方程 70 6.5.1 修正強奇異邊界積分方程 71 6.5.2 修正超奇異邊界積分方程 75 6.6 綜合討論 77 第七章 SGBEM於力邊界值問題之分析例 80 7.1 內域問題－立方體受單軸張力 80 7.2 外域問題－無限域中球形空穴受內壓 81 7.3 半無限域問題 83 7.3.1 半無限域問題之分析模式 83 7.3.2 半無限域表面受鉛垂向圓形均佈載重 84 第八章 FEM~SGBEM連結模式於力邊界值問題之分析例 110 8.1 內域問題－單軸拉力問題 110 8.1.1 矩形體受單軸張力 110 8.1.2 立方體受單軸張力 112 8.2 半無限域問題 113 8.2.1 半無限域表面受鉛垂向圓形均佈載重 114 8.2.2 半無限域表面圓形剛性版受鉛垂向載重 115 8.2.3 半無限域中單樁受鉛垂向載重 117 第九章 結論與建議 143 9.1 結論 143 9.2 建議 144 參考文獻 145 附錄一 無限域靜力基本解 153 A1.1 無限域基本解各核函數 153 A1.2 無限域基本解各核函數之對稱特性 154 附錄二 積分領域之等義表示式與投影及分割關係 155 A2.1 三角形積分領域等義表示式之一 155 A2.2 三角形積分領域等義表示式之二 155 A2.3 三角形積分領域之投影關係 155 A2.4 四邊形積分領域之分割 156 附錄三 奇異積分之變數轉換操作 158 A3.1 完全重合情況(coincident case) 158 A3.2 邊緣相鄰情況(edge-adjacent case) 159 A3.3 頂點相鄰情況(vertex-adjacent case) 1611921236 bytesapplication/pdfen-US葛勒金對稱邊界元素法三維靜力分析無限域靜力基本解奇異積分有限元素∼邊界元素連結解非唯一半無限域問題symmetric Galerkin BEM3D elastostaticssingular integralsFE~BE couplingnon-uniquenesshalf-space problemsthesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/50405/1/ntu-95-F90521105-1.pdf