破坏科学研讨会2004Keeting.PPT

第八屆破壞科學研討會 海洋大學力學聲響振動實驗室 對偶邊界元素法在疲勞破壞之工程應用 時間：2004/03/26.27 地點：墾丁福華 報告者：陳正宗.沈文成 有限元素法網格 傳統邊界元素法網格 傳統邊界元素法網格 以邊界元素法發表的文章 對偶邊界元素法網格 大綱 對偶邊界元素積分方程式 對偶邊界元素法 簡單測試例 成功應用的工程實例 結論 功能互換原理 今有二個彈性體，其變位場、面力場 與徹體力場分別為(ui ,ti ,bi)與(ui* ,ti* ,bi*) 由 Betti 功能互換原理，可知 域內與域外點對偶邊界積分式 考慮 ui* , ti* , bi* 為基本解系統 對偶邊界積分方程式 將場點 x 推到邊界 對偶邊界積分方程式 對偶邊界元素法 簡單測試例 成功應用的工程實例 結論 對偶邊界元素法 將邊界離散化 矩陣式 對偶式 二維裂縫問題模型 對偶邊界積分方程式 對偶邊界元素法 簡單測試例 成功應用的工程實例 結論 簡單測試例 特性資料分述如下： 楊氏係數=10×106 psi 卜松比 n =0.3 負荷重 P =1000 kgf/mm2 合理地束制剛體運動。 MSC/NASTRAN 有限元素法分析 二維有限元素綱目分佈圖(CRAC2D) 三維有限元素綱目分佈圖(CRAC3D) 邊界元素法求解 邊界元素綱目圖 內點分佈圖 應力強度因子結果比較 有限元素法與邊界元素法簡介 對偶邊界積分方程式 對偶邊界元素法 簡單測試例 成功應用的工程實例 結論 含應力消除塊固體火箭發動機 固體火箭發動機 裂縫發生處與成長預估 V 型環 邊界元素法網格 有限元素法網格 V 型環 V 型環結構變形圖 裂縫成長路徑 應力強度因子 結論 對偶邊界積分方程式與對偶邊界元素法，已成功地應用含應力消除塊固體火箭發動機、飛彈組件 V 型環。 對偶邊界元素法與傳統的邊界元素法不同的是，無需多領域切割即可進行破壞分析。 在裂縫成長預估，更無需網格重切割，在前處理建模時相當方便，對於實際工程應用將有很大的發揮空間。 * * 有限元素法網格 裂縫成長後之有限元素法網格 傳統邊界元素法網格 裂縫成長後之傳統邊界元素法網格 傳統邊界元素法網格 裂縫成長後之傳統邊界元素法網格 對偶邊界元素法網格 裂縫成長後之對偶邊界元素法網格 其中，D 表示定義域，B 表示邊界 其中， U(s,x) 與 T(s,x) 分別代表位移與面力基本解 Di、De 分別表示域內與域外 s 為源點，x 為場點 B Di De 其中， CPV、RPV、HPV 分別表示柯西主值、 黎曼主值與阿達馬主值 x 其中， N 表示元素的點數，Bm 表示第 m 個元素 UT UT 式：位移式 LM 式：面力式 裂縫尖端 LM 裂縫上緣 裂縫下緣 裂縫上緣 裂縫下緣 UT 式 LM 式 4mm 12mm 1mm P 2066 （0.482%） 2021 （2.7%） 2138 （2.9%） KI* （誤差值） 36 37 137 節點數 18 18 18 元素數 對偶邊界元素法 一維元素模擬 有限元素法 二維元素模擬 有限元素法 三維元素模擬 *KI（正確解）=2077.5 實線：數值 虛線：實驗 6.561 8.748 KII 34.874 36.119 KI 邊界元素法 有限元素法