专题研究书面报告-NTOU-海洋大学力学声响振动实验室.docVIP

PAGE 專題研究書面報告 摘要 在本研究計畫中，探討利用「對偶積分元素法」來解決尤拉樑之變形分析的問題。 首先對一維Laplace桿問題為範例來探討邊界積分方程。藉由桿的分析，從已知控制方程式，利用分段法推得其基本解及退化核，再藉由功能互換定理及分部積分可得到對偶積分方程，由對偶架構的概念對s（source point）微分，其中，我們可得核函數的關係及其退化核，將之代入積分方程，整理成矩陣形式，而成為我們可知的影響係數矩陣與對偶矩陣。 再將Laplace桿問題推廣至Biharmonic方程的樑問題。在研究計畫中將推導求得樑之勁度矩陣，進而建構出場解表示式。 此計畫以邊界積分方程來學習一維之Laplace方程的桿及Biharmonic方程的樑問題，深入淺出、內容完整，將來可成為許多初學者研習邊界積分方程或邊界元素法相關領域之入門教材，此研究主題深具教育與學術價值，且可將材料力學與結構學中所得的結果與邊界積分方程作一比較，及其在數學與物理方面之意義加以洞悉，對於學生有溫故知新之效，亦可成為搭建工程數學與結構學的橋樑。 研究動機與研究問題 在修習過材料力學及結構學後，發現在求解尤拉樑之傾角、撓度、彎矩及剪力時，可由很多方法求得，如：傾角撓度法、共軛樑法與直接勁度法…等。在此我們將提供一套由數學推導結合物理觀念的方法來求解尤拉樑之的變形分析。 此研究在於利用邊界積分方程來求解尤拉樑之變形分析，我們可建構出尤拉樑之各種場解表示式。結合基本解及退化核之數學技巧，我們將做解析推導，由推導的過程中我們亦可求得樑之勁度矩陣；配合材料力學與結構學所學的知識，將尤拉樑的勁度矩陣與場解表示式做一完整的分析討論，並與教科書做一對照。 本計畫將針對樑之自然邊界條件(natural boundary condition) 及基本邊界條件(essential boundary condition) 及混合型邊界條件來探討固定端、自由端與簡支端的邊界。藉著計畫中理論建立與公式推導之過程，對於邊界積分方程之相關背景知識，可做一初步的入門與了解，為往後研究生涯做準備，並培養獨立思考的能力，學習研究的方法，建立日後研究能力的良好基礎。 研究方法及步驟 首先，以一維Laplace問題為範例來學習邊界積分方程，說明如下圖一： 圖一(a)位移自由度圖一(b)軸力自由度 圖一(a)位移自由度 圖一(b)軸力自由度 考慮一桿之問題，其控制方程式為： , 為桿之軸向變形量,為欲解問題之領域，並引入一輔助系統： 其中，為Dirac Delta 由分段法可得基本解 及退化核 藉由功能互換定理及分部積分技巧，可得： 　　 　　 可得邊界積分方程： （1） 由對偶架構的概念將式（1）對s微分可得二條邊界積分方程： （2） （3） 其中，，我們可得四個核函數的關係，如下所示： 四個核函數之函數圖形分別如下所示 四個核函數的退化核，如下所示 由（1）式（） 當，可知 整理可得 (4) 當，可知 整理可得 (5) 由(4)、(5)式可得如下的矩陣形式 (6) 由(6)式 (7) 由（3）式（） 當，可知 整理可得 (8) 當，可知 整理可得 (9) 由(8)、(9)式可得 (10) 將(6)、(10)式表為下列形式： 可得表一 表一：一維桿對偶邊界積分元素法推導之影響係數矩陣與對偶矩陣 Equation [A] [B] [K]= [B]-1[A] （4）、（5） Rank=1 Rank=2 （8）、（9） Rank=1 NA 將Laplace方程的桿問題推廣至Biharmonic方程的樑問題，圖一的自由度變成圖二 圖二(b)樑之廣義力自由度圖二(a)樑之位移自由度 圖二(b)樑之廣義力自由度 圖二(a)樑之位移自由度 考慮一尤拉樑問題，其控制方程式為： ，其中，為樑長 為樑之側向位移,為欲解問題之領域,並引入一輔助系統： 藉由功能互換定理及分部積分之技巧： 可得邊界積分方程： (11) 由對偶架構的概念可得四條邊界積分方程： (12) (13) (14) (15) 為簡化問題，此處將EI設為1，則 其中，分別代表撓度、傾角、彎矩、剪力。 因此，我們可得16個核函數的關係，如下所示： 16個核函數的退化核，如下所示，可得表二～五 表二：說明十六個核函數表示式 表三：說明十六個核函數表示式 表四：說明十六個核函數表示式 表五：說明十六個核函數表示式 由(12)～(15)式任取2式，並配合核函數退化核表二～五，表為下列矩陣形式，可得表六，其中空格處即為本研究計畫要完成的 研究結果 表六 Equation [A] [B] [K]= [B]-1[A] (12)、(13) Rank=2 Rank=4 (12)、(14) Rank=2 Rank=4 (12)、(1