运用格林函数求解力学问题-海洋大学.PDF

運用格林函數求解力學問題 海洋大學河海工程學系 報 告 者：柯佳男 授課教授：曹登皓教授 專題討論四 (結構組) 大綱 引言 修改勢能法 零場積分方程 結論 大綱 引言 修改勢能法 零場積分方程 結論 文獻的回顧 MMP法稱為修改勢能法 MMP法第一次被採用在西方的有關計算力學的 文獻上[Melnikov(1977)], 那次結果的觀念較後 被應用[Melnikov(1982)] MMP被證明更有效在許多的實例中相較於傳統 的勢能法MCP[Brebbia(1978),Kupradze(1965), Lovitt(1950), Smirnow(1964), Tikhonov and Samarskii(1963)] 動機 當計算格林函數對於邊界值問題，陳述 Laplace 和 biharmonic 方程在二維複雜 結構的區域上，使用勢能的表現是值得 研究的 勢能的修改，在某些邊界線段不需要去 解邊界積分方程 大綱 引言 修改勢能法 零場積分方程 結論 積分方程 ∫k (x, s)φ(s)dB(s) ϕ(x) 核函數 未知函數 1.基本解 2.格林函數 3.譜積分之形式 修改勢能法 1. 在多連通區域尋找一個格林函數 G( z,ζ) 2. 將格林函數分解成奇異部份和規則部分 G( z,ζ) G ( z,ζ) +g (z ,ζ) 0 修改勢能法 * 3. g (z ,ζ ) 必須滿足邊界條件 * g (z ,ζ ) 0, z ∈Λ g (z , * ) G (z , * ), z ζ − 0 ζ ∈Γ 4. 將它帶入積分的方程就得 到新的勢能表現 * g (z ,ζ ) ∫Γ G (z ,t )q (t )d Γ, z ∈Ω 0 0 t 0 * ∫Γ G ( z, t)q(t)d Γ =−G( z,ζ ), z ∈Ω 0 0 t 0 比較 MMP法 MCP法 核函數不是控制微分 核函數是控制微分方 方程的基本解 程的基本解 場點和源點在不同項 在數值計算處理上較 為直觀 力學應用上較多的發 展 數值算例(熱傳導問題) 在二維的平面上，將物體分割為n條狀且含一個洞L 熱傳導問題(續) 溫度u在條狀上滿足下列關析式 2 ∂ ∂u ∂u (λ (x ) k ) +λ (x ) k f (x , y ) (k 1,2, ,n) ∂x k ∂x k ∂y 2 k Au 0; A u 0;