Please use this identifier to cite or link to this item:
http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/28617
Full metadata record
???org.dspace.app.webui.jsptag.ItemTag.dcfield??? | Value | Language |
---|---|---|
dc.contributor.advisor | 陳俊杉(Chuin-Shan Chen) | |
dc.contributor.author | Wei-Lin Lo | en |
dc.contributor.author | 羅威麟 | zh_TW |
dc.date.accessioned | 2021-06-13T00:14:25Z | - |
dc.date.available | 2007-07-30 | |
dc.date.copyright | 2007-07-30 | |
dc.date.issued | 2007 | |
dc.date.submitted | 2007-07-26 | |
dc.identifier.citation | Akanda, M.A.S., Ahmed, S.R. and Uddin, M.W. (2000), A finite-difference scheme for mixed boundary value problems of arbitrary-shaped elastic bodies, Advances in Engineering Software, Vol. 31(3),173-184.
Mitchell A.R. and Griffiths D.F., (1980), The Finite Difference Method in Partial Differential Equations. John Wiley & Sons. Liggett, J.A. and Liu, P. L.-F. (1983), The Boundary Integral Element Method for Porous Media Flow. Allen and Unwin, London. Thompson, J.F., Warsi, Z.U.Z. and Mastin, C.W. (1985), Numerical Grid Generation: Foundations and Applications. North Holland, New York. Timoshenko, S.P. and Goodier, J.N. (1986), Theory of Elasticity. 3rd ed. McGraw-Hill Book Company, New York. Trim, D.W. (1995), Complex Analysis and Its Applications. PWS Publishing Company. Tsay, T.K., Ebersole B.A. and Liu, P. L.-F. (1989), “Numerical Modeling og Wave Propagation Using Parabolic Approximation with a Boundary-Fitted Co-ordinate System,” International Journal for Numerical Methods Engineering, Vol. 27, 37-55. Tsay, T.K. and Hsu, F.S. (1997), “Numerical Grid Generation of an Irregular Region,” International Journal for Numerical Methods in Engineering, Vol. 40, 343-356. Tsay, T.K., Yeh, G.T., Wilson, G.V. and Toran, L.E. (1990), GRIDMAKER: ”a Grid Generator for Two- and Three-dimensional Finite Element Subsurface Flow Models.” Oak Ridge National Laboratory, Tennessee. Tsay, T.K., Wang, J. and Huang, Y.T. (2006), “Numerical generation and grid controls of boundary-fitted conformal grids in multiply connected regions,” International Journal for Numerical Methods in Engineering, Vol. 67, 1045-1062. Wang J. (2004), On the Application of BEM to Boundary-fitted Conformal Grid Generation and Associated Singularity Problems. PhD Thesis, National Taiwan University. Wang, J., Tsay, T.K. (2005), “Analytical evalution and application of singularities in boundary element method,” Engineering Analysis with Boundary Elements, Vol. 29, 241-256 王存欵,「土石流之數值模擬」,1996,國立台灣大學土木工程學研究所碩士論文。 吳萬崧,「數值網格在流場計算的應用」,1999, 國立台灣大學土木工程學研究所碩士論文。 林齊堯,「暴潮數值模式之硏究」,2005,國立台灣大學土木工程學研究所碩士論文。 徐福盛,「區域數值轉換之研究」,國立台灣大學土木工程學研究所碩士論文,1995。 黃永德,「複連通區域多區塊之符合邊界正交網格建立」,1998,國立台灣大學土木工程學研究所碩士論文。 劉宜華,「應用區域轉換法反算不規則區域之邊界不恆定溫度分佈」,2000,國立台灣大學土木工程學研究所碩士論文。 鄭正德,「邊界符合座標在暴潮模式之應用」,2002,國立台灣大學土木工程學研究所碩士論文。 蔡季陸,「有限差分法結合邊界符合正交網格於三維地下水流模擬之應用」,2001,國立台灣大學土木工程學研究所碩士論文。 賴俊達,「自動調整微小網格在污染物傳輸計算之研究」,1993,國立台灣大學土木工程學研究所碩士論文。 | |
dc.identifier.uri | http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/28617 | - |
dc.description.abstract | 本研究旨在應用Tsay等人(Tsay and Hsu, 1997; Wang and Tsay, 2005; Tsay, Wang and Huang, 2006; 黃永德,1998)發展出來的「數值全域轉換法」(Numerical Global Domain Transformation Method, NGDTM),在任意二維幾何區域求解彈性力學的問題(Airy stress function),並發展出一套通用的數值求解策略。
「數值全域轉換法」為一個將整體解題空間由物理發生的平面轉換至另一個平面求解的方法。利用複變映射原理(Complex mapping)以及邊界元素積分法(Boundary integral element method)求解Laplace方程式,將一個任意二維幾何形狀轉換至一矩形區域,使得有限差分法(Finite difference method)得以在此矩形區域方便地運用。 Airy stress function為求解彈性力學的控制方程式之ㄧ(Timoshenko and Goodier, 1986),由於定義的關係,traction邊界條件只能提供二階的偏導數值。過去在矩形空間應用FDM求解Airy stress function時,必須透過沿著幾何邊界對traction邊界條件積分獲得一階及零階偏導數的數值邊界條件(Timoshenko and Goodier, 1986)。然而當幾何邊界複雜時,這個積分的動作會變得相當地困難,並且程式不容易撰寫成一般化,因此未被廣泛地應用在數值方法。 本研究透過Tsay等人發展出來的區域數值轉換理論,將任意的幾何邊界條件及traction邊界條件轉換至矩形區域,然後在該矩形區域上將traction邊界條件透過數值積分轉換成運用有限差分法需要之數值邊界條件。控制方程式亦透過此區域數值轉換理論由不規則區域轉換至矩形區域。運用轉換過的控制方程式和邊界條件在矩形區域進行有限差分法之計算。進行控制方程式和邊界條件的轉換時,將透過「柯西-里曼方程式」(Cauchy-Riemann conditions)化簡方程式,並且轉換後的Jacobian, 、 不為常數。在矩形區域進行離散時,本研究採取在區域內(不含邊界)設置以中央差分離散之控制方程式,邊界點上的自由度每個點設置兩條邊界條件,四個角落分別設置三條邊界條件的策略,使得未知數的數目等於方程式的數目,獲得唯一且收歛的解。以一兩端承受彎矩的栱型樑為例,轉換至矩形區域後,在網格點為 ,網格大小為 的情況下,取中間的剖面的數值結果與解析解相較,勢能(Potential)的誤差在 之內,應力 的誤差在 之內,應力 的誤差在 之內,應力 最大誤差為0.004(解析解為0)。 本研究推導出控制方程式、邊界條件及應力式正確的轉換關係,並提出一套通用的數值求解策略,奠定Tsay等人發展出來的「數值全域轉換法」在彈性力學領域應用的基礎。同時經由本研究建立之求解模式,有限差分法可以很方便地運用在二維彈性力學的求解上,提昇了運用有限差分法求解二維彈性力學時適用的廣度。 | zh_TW |
dc.description.abstract | The objective of this research is to develop a general numerical scheme to solve elasticity problems in an arbitrary two dimensional domain by “Numerical global domain transformation method (NGDTM),” a domain transformation theorem developed by Tsay et al. (Tsay and Hsu, 1997; Wang and Tsay, 2005; Tsay, Wang and Huang, 2006; 黃永德,1998).
NGDTM transforms the problem solving domain from an arbitrary two dimensional domain, the physical domain, to a rectangular domain by applying complex mapping theorem and solving Laplace equations by boundary integral element method (BIEM). When the domain is transformed to a rectangle domain, it becomes very convenient to solve the problems with finite difference method (FDM). Airy stress function is one of the governing equations in elasticity (Timoshenko and Goodier, 1986). Based on the Airy stress function, the traction boundary conditions can only provide the second order derivative of the primary variable. It is thus necessary to integrate the traction boundary condition along the geometry to obtain the numerical boundary condition when one tries to solve the Airy stress function in a rectangular domain (Timoshenko and Goodier, 1986). Such integration becomes relatively difficult when the geometry is complex. Partly due to such shortcoming, FDM is not popular in solving elasticity problems. In this work, an arbitrary two dimensional domain and traction boundary conditions were transformed to a rectangular domain through the domain transformation theorem. Traction boundary conditions in the rectangular domain were integrated to obtain the numerical boundary condition. The governing equation was also transformed form the arbitrary domain to the rectangular one by the domain transformation theorem. After applying the transformed governing equation and numerical boundary conditions, the finite difference method calculation was performed in a rectangular domain. During the transformation derivative, Cauchy-Riemann conditions were applied in which the Jacobians and , were not constants. The governing equation and numerical boundary conditions were both expressed in the first-order central difference. The governing equation was discreted to a 13 point formulation and applied on the internal meshes excluding the boundary. Each node on the boundary has two boundary conditions. Each corner has three boundary conditions. This scheme makes the number of unknowns equal to the number of equations and result in an unique and convergent solution. An arc with a couple of moment applied was used to verify the proposed scheme. The rectangle domain was discreted in nodes and the mesh size was. . In comparison with the analytical solution on the middle of the arc, the error of the potential was less than , the error of the normal stress was less than , the error of the normal stress was less than , and the maximum error of the shear stress was 0.004 (the analytical solution is zero). A NGDTM for the Airy stress function has been developed in this work. A transformed governing equation, boundary conditions, stress formulation and a general numerical scheme have been derived. Through the machinery of NGDTM developed herein, the FDM can be used conveniently to solve two-dimensional elasticity problems. | en |
dc.description.provenance | Made available in DSpace on 2021-06-13T00:14:25Z (GMT). No. of bitstreams: 1 ntu-96-R92521604-1.pdf: 4748849 bytes, checksum: 2bb2345ffe4bd83c2974dbb55d17d7e0 (MD5) Previous issue date: 2007 | en |
dc.description.tableofcontents | 致謝 ……………………………………………………………………………… i
摘要 ……………………………………………………………………………… v Abstract …………………………………………………………………………… vii 目錄 ……………………………………………………………………………… ix 圖目錄 …………………………………………………………………………… xi 表目錄 …………………………………………………………………………… xiii 第一章 緒論 …………………………………………………………………… 1 1-1 研究背景與動機 ……………………………………………………… 1 1-2 研究目的 ……………………………………………………………… 2 1-3 文獻回顧與相關研究 ………………………………………………… 3 第二章 二維彈性力學分析 …………………………………………………… 9 2-1控制方程式、邊界條件與應力表示式 ……………………………… 9 2-2 控制方程式、邊界條件與應力表示式之轉換 ……………………… 11 第三章 利用全域轉換法的解析求解 ………………………………………… 23 3-1 推導目的 ……………………………………………………………… 23 3-2 解析求解範例 ………………………………………………………… 23 第四章 數值方法 ……………………………………………………………… 31 4-1 求解流程 ……………………………………………………………… 31 4-2 不規則區域和矩形區域的物理邊界條件轉換關係 ………………… 32 4-3 數值積分求一階及零階偏導數邊界條件 …………………………… 36 4-4 控制方程式與邊界條件的離散 ……………………………………… 39 4-5 自由度的分布 ………………………………………………………… 42 4-6 不規則區域的應力式在矩形區域的表示式 ………………………… 43 4-7 矩形區域計算策略 …………………………………………………… 47 第五章 數值方法之驗證 ……………………………………………………… 57 5-1 解析解 ………………………………………………………………… 57 5-2 本文數值解的前置作業與流程 ……………………………………… 59 5-3 數值解之驗證 ………………………………………………………… 61 第六章 結論與建議 …………………………………………………………… 67 6-1 結論 …………………………………………………………………… 67 6-4 建議 …………………………………………………………………… 68 參考文獻 ………………………………………………………………………… 69 附錄A 偏導數的相關推導 ……………………………………………………… 73 附錄B 本文極座標下偏導數的轉換推導 ……………………………………… 79 附錄C 本文環狀區域之位移 …………………………………………………… 83 附錄D 控制方程式的差分式 …………………………………………………… 87 作者簡曆……………………………………………………93 | |
dc.language.iso | zh-TW | |
dc.title | 數值全域轉換法於二維彈性力學之應用 | zh_TW |
dc.title | Applications of Numerical Global Domain Transformation Method to 2D Elasticity Problems | en |
dc.type | Thesis | |
dc.date.schoolyear | 95-2 | |
dc.description.degree | 碩士 | |
dc.contributor.coadvisor | 蔡丁貴(Ting-Kuei Tsay) | |
dc.contributor.oralexamcommittee | 葉超雄(Chau-Shioung Yeh),吳光鐘(Kuang-Chong Wu) | |
dc.subject.keyword | 數值全域轉換法,區域轉換,邊界符合保角網格,保角網格,彈性力學,有限差分法,數值計算, | zh_TW |
dc.subject.keyword | Numerical global domain transformation method,NGDTM,Boundary-fitted conformal grid,Conformal grid,Airy stress,Biharmonic,Finite difference method,FDM,Numerical solution, | en |
dc.relation.page | 71 | |
dc.rights.note | 有償授權 | |
dc.date.accepted | 2007-07-27 | |
dc.contributor.author-college | 工學院 | zh_TW |
dc.contributor.author-dept | 土木工程學研究所 | zh_TW |
Appears in Collections: | 土木工程學系 |
Files in This Item:
File | Size | Format | |
---|---|---|---|
ntu-96-1.pdf Restricted Access | 4.64 MB | Adobe PDF |
Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.