应用对称葛勒金边界元素法求解静力纽曼问题之技巧UseofSymmetric.PDFVIP

九十六年電子計算機於土木水利工程應用研討會 應用對稱葛勒金邊界元素法求解靜力紐曼問題之技巧 柯永彥 1 陳正興 2 摘要 應用對稱葛勒金邊界元素法 (SGBEM)分析紐曼問題(Neumann problems)時，所得位 移解將會存在有一剛體運動項， 而有解非唯一之現象。本文將討論此現象之內涵與成 因，並彙整一些能去除剛體運動影響的方法，以實例分析探討其效能與適用性。在大地 工程領域常見之半無限域問題，多屬紐曼問題，本研究藉由所彙整可獲唯一解之方法， 利用以無限域基本解建立之 SGBEM模式，採自由表面截切技巧來分析半無限域問題， 但因截去遠域自由 表面造成之誤差，使得結果有若干不合理之現象。其中，以基於 Fredholm理論建立之修正邊界積分方程，所得結果相對較佳。 關鍵字 : SGBEM ，三維靜力分析，紐曼問題，解非唯一，半無限域問題 Use of Symmetric Galerkin Boundary Element Method for Elastostatic Neumann Problems Yung-Yen Ko 1 and Cheng-Hsing Chen 2 ABSTRACT When Symmetric Galerkin boundary element method (SGBEM) is applied to Neumann problems, the displacement solution cannot be uniquely determined because of the inevitable rigid body motion term involved. In this paper, discussions on this phenomenon and methods for the removal of the rigid-body-motion terms are investigated by case study. In geotechnical engineering, half-space problems are dealt with quite often, and most of them are Neumann problems. In this study, half-space problems are analyzed by SGBEM by modeling free surface with limited number of elements. Nevertheless, the solutions obtained are still not satisfactory because of the errors due to the truncation of the free surface. Among the methods investigated, the modified boundary integral equations based on Fredholm theory is relatively preferable. Keywords: SGBEM, 3D elastostatics, Neumann problem, non-uniqueness, half-space problems 一、前言 素法(CBEM)具較佳精確性與收斂性[2] 。 柯永彥與陳正興 (2005, 2006)提出一 對稱葛勒金邊界元素法(SGBEM)採 套可分析三維彈力問題之 SGBEM 模式 對偶邊界積分方程(dual BIEs)[1] ，配合加 [3] ，其中整合[4][5]所採方法，利用座標轉 權剩餘法之葛勒金法 (Galerkin method) ，得 換與積分域投影技巧，配合基本解特性來 到雙重邊界積分