吳光鐘臺灣大學：應用力學研究所張宗堯Chang, Tsung-YaoTsung-YaoChang2007-11-292018-06-282007-11-292018-06-282006http://ntur.lib.ntu.edu.tw//handle/246246/62355本文所要推導的是，無邊際（ Unbounded ）均質( Homogeneous )異向彈性體( Anisotropic elastic materials ) 於二維（ Two - dimension ）狀況下承受線性時間週期性荷重（ Line time - harmonic load ），或於三維（ Three - dimension ）狀況下承受點時間週期性荷重（ Point time - harmonic load ）之格林函數（ Green’s function ）。此格林函數需要處理積分表示和邊界積分的問題，可以運用邊界元素法( Boundary element method ) 推導此格林函數。由於邊界元素法在處理等向彈性體分析波傳( Wave propagation) 方面問題時，扮演著相當重要的角色，所以希望藉由邊界元素法在處理異向彈性體於彈性動力方面問題時，能與處理等向彈性體一樣的容易。本文將以推導所得到的格林函數，於等向彈性體材料與參考文獻上的解析解比對，以確認推導所得到的格林函數其正確性。於二維異向彈性體材料，推估格林函數其數據結果；於三維異向彈性體材料，與參考文獻比較數據結果是否吻合。另外運用傅立葉轉換的方式，解決二維及三維異向彈性體受時間週期性荷重狀況下所會遭遇到的問題，可以證明除了運用Radon轉換來推導格林函數，當然也可以使用傅立葉轉換來處理相同狀況的問題。本文期望運用座標轉換方式，先將材料性質轉換到新座標系，在得到轉換成新座標系後的結果，再次運用座標轉換的方式，轉換回原始座標系，即可得到原始座標系下，所欲求得的結果，亦即經由座標轉換，來處理觀察點方向不同時的狀況。Abstract Generally speaking , no matter which kind of material want to ask the displacement or the strain , it is the way that tries to get the Green’s function under this state . This paper is to determine the Green’s function in an unbounded homogeneous anisotropic media generated by the application of a two-dimensional time-harmonic point load , or the application of a three-dimensional time-harmonic line load . This paper based on the Fourier transform is presented to determine the Green’s function in a anisotropic media due to the application of a time-harmonic load .This Green's function needs to deal with the problem of integral representations and boundary integral equations by the boundary element method . Because the boundary element method has been playing a very important role in analyzing wave propagation problems in isotropic media . It hopes to deal the elastodynamic problems for anisotropic media as easy as isotropic media . In the isotropic media , it needs to estimate its data result of the Green's function with the data result of the analytical solution . This paper expects to use the coordinate transformed to deal the state with different viewpoint directions.目錄 摘要 I Abstract III 目錄 IV 圖目錄 VI 表目錄 VIII 符號表 IX 第一章、導論 - 1 - 1-1 研究動機 - 1 - 1-2 文獻回顧 - 2 - 1-3 研究目的與內容大綱 - 4 - 第二章、二維格林函數理論推導 - 7 - 2-1 異向彈性體 - 7 - 2-2 等向彈性體 - 16 - 2-3 無因次化 - 18 - 第三章、三維格林函數理論推導 - 23 - 3-1 異向彈性體 - 23 - 3-2 等向彈性體 - 40 - 3-3 無因次化 - 49 - 第四章、格林函數數值分析 - 54 - 4-1 高斯N點積分公式 - 54 - 4-2 二維等向彈性體之格林函數數值分析 - 55 - 4-3 二維異向彈性體之格林函數數值分析 - 57 - 4-4 三維等向彈性體之格林函數數值分析 - 59 - 4-5 三維異向彈性體之格林函數數值分析 - 61 - 第五章、結論與未來展望 - 65 - 5-1 結論 - 65 - 5-2 未來展望 - 68 - 參考文獻 - 70 - 圖目錄 圖2.1線性時間週期性荷重與座標系關係圖 - 72 - 圖2.2積分輪廓圖 - 72 - 圖3.1點時間週期性荷重與座標系關係圖 - 73 - 圖3.2 k座標系與原始座標系關係圖 - 73 - 圖3.3點時間週期性荷重與新座標系關係圖 - 74 - 圖3.4兩座標系、x與點時間週期性荷重的關係圖 - 74 - 圖4.1二維等向彈性體格林函數於X-X實數部分圖 - 75 - 圖4.2二維等向彈性體格林函數於X-X虛數部分圖 - 75 - 圖4.3二維等向彈性體格林函數於Y-Y實數部分圖 - 76 - 圖4.4二維等向彈性體格林函數於Y-Y虛數部分圖 - 76 - 圖4.5二維等向彈性體格林函數於Z-Z實數部分圖 - 77 - 圖4.6二維等向彈性體格林函數於Z-Z虛數部分圖 - 77 - 圖4.7二維異向彈性體格林函數於X-X實數部分圖 - 78 - 圖4.8二維異向彈性體格林函數於X-X虛數部分圖 - 78 - 圖4.9三維等向彈性體格林函數於Z-Z， 圖 - 79 - 圖4.10三維等向彈性體格林函數於Z-Z， 圖 - 79 - 圖4.11三維等向彈性體格林函數於Z-Z， 圖 - 80 - 圖4.12三維等向彈性體格林函數於Z-Z， 圖 - 80 - 圖4.13三維異向彈性體（冰）於Z-Z， 圖 - 81 - 圖4.14三維異向彈性體（冰）於Z-Z， 圖 - 81 - 圖4.15三維異向彈性體（冰）於Z-Z， 圖 - 82 - 圖4.16三維異向彈性體（冰）於Z-Z， 圖 - 82 - 圖4.17三維異向彈性體（石墨環氧）於Z-Z， 圖 - 83 - 圖4.18三維異向彈性體（石墨環氧）於Z-Z， 圖 - 83 - 圖4.19三維異向彈性體（石墨環氧）於Z-Z， 圖 - 84 - 圖4.20三維異向彈性體（石墨環氧）於Z-Z， 圖 - 84 - 表目錄 表4.1二維等向彈性體與解析解於X-X數據結果比較 - 85 - 表4.2二維等向彈性體與解析解於Y-Y數據結果比較 - 86 - 表4.3二維等向彈性體與解析解於Z-Z數據結果比較 - 87 - 表4.4碳環氧材料彈性材料常數 - 88 - 表4.5三維等向彈性體格林函數與參考文獻[1]比較， - 88 - 表4.6三維等向彈性體格林函數與參考文獻[1]比較，N＝10 - 89 - 表4.7三維等向彈性體格林函數與參考文獻[1]比較，N＝20 - 90 - 表4.8三維等向彈性體格林函數與參考文獻[1]比較，N＝30 - 91 - 表4.9冰彈性材料常數 - 91 - 表4.10石墨環氧複合材料彈性材料常數 - 92 -1097165 bytesapplication/pdfen-US時間週期性荷重格林函數傅立葉轉換邊界元素法座標轉換高斯 N點積分法the time-harmonic loadthe Green’s functionthe Fourier transformthe boundary element methodthe coordinate transformedthe Gauss-Legendre N-point quadrature formulathesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/62355/1/ntu-95-R93543051-1.pdf