專題研究心得.docVIP

國立台灣海洋大學 河海工程學系 專題研究(一) 心得 姓名:謝正昌 學號:系級:河工三A 日期:2005/06/30 專題研究心得: 在專題研究(一)課程裡，真的學到了不少東西，每個星期到實驗室總能學到一些東西，總是把老師的研究室比喻為一座寶山，秉持一個信念，入寶山豈能空手而回。在研究室的學習中，我也接觸了Mathematica使用，一開始是自己看書操作，計算專題中的矩陣運算， 真的快速多了，接下來就是Mathematica在工數(四)上的運用，令我覺得無比的神奇，從2-D圖形變成3-D動畫，讓我對波的傳遞行為，有更進一步的認識，雖然在寫運算過程中仍要翻著書對照才能輸入運算符號，不過由自己寫出來地運算過程跑出正確的圖形是一件無比興奮的事。 這學期一開始的重點是國科會大學專題計畫，雖然是利用學長沒成功的計畫，不過還好有嘉俊學長的幫忙，耐心的指導，他給了我一本書，是老師和洪宏基教授合著的邊界元素法，他指定了幾個和我專題有關的章節回去念，一開始很挫折的，所以只好再求助於學長，慢慢的也了解那些數學推導式子，在專題一開始的一維Laplace桿問題，自己也從兩條邊界積分方程找出四個核函數的關係，推導出影響係數矩陣與對偶矩陣，並發現我的答案和之前學長做的不一樣，馬上和嘉俊學長一起討論，發現我做的才是對的，心裡覺得很有成就感，讓我覺得的越來越有興趣；緊接著推廣到Biharmonic方程的樑，這個就更麻煩了，從四條邊界積分方程找出十六個核函數的關係，同樣的推導出影響係數矩陣與對偶矩陣，這也就是對偶積分方程中之直接法。 在我的專題研究裡增加了一種方法，那就是對偶積分方程中之間接法，這個方法是在專題截稿三天前，老師要我增加的新元素，不斷和嘉俊學長討論，參考老師開的碩士班課程—邊界元素法，學長做的作業裡去摸索，由於有嘉俊學長從旁指導也慢慢了解了，利用直接法得到的16個核函數，然後配合不同的核函數，共有單雙層勢能法、單三層勢能法、單四層勢能法、雙三層勢能法、雙四層勢能法、三四層勢能法，找出各解法的矩陣形式即可。 在這一學期的過程中，學到不少東西，尤其深深體認到做事要有謹慎的態度，相信對以後自己在處理事情時會有幫助的，當自己在做專題時，體驗到作研究的辛苦，雖然過程是艱辛的，可是如果擁有強烈的求知欲和極大的熱誠，那就是作研究的最大動力。 接下來是我整理的專題研究內容: 一維Laplace桿問題，邊界積分方程： (1) 由對偶架構的概念將式（1）對s微分可得二條邊界積分方程： 四個核函數的退化核 一維桿對偶邊界積分元素法推導之影響係數矩陣與對偶矩陣 Equation [A] [B] [K]= [B]-1[A] （4）、（5） Rank=1 Rank=2 （8）、（9） Rank=1 NA (12) 直接法 先利用對偶架構的概念可得四條邊界積分方程: (13) (14) (15) (16) 說明十六個核函數表示式 由(13)～(16)式任取2式，並配合核函數退化核表二～五，表為下列矩陣形式，此為對偶積分方程之直接法，藉由此法可得表六，其中空格處即為本研究計畫要完成的。 研究結果 (17) 表六 Equation [A] [B] [K]= [B]-1[A] (13)、(1) Rank=2 Rank=4 (13)、(1) (13)、(1) (14)、(1) (14)、(1) (15)、(1) 間接法 由(13)～(16)式，配合不同的核函數，共有單雙層勢能法、單三層勢能法、單四層勢能法、雙三層勢能法、雙四層勢能法、三四層勢能法。 配合核函數退化核表二～五，表為下列矩陣形式，由此法可得表七，其中空格處即為本研究計畫要完成的。 研究結果 (24) 表七 勢能 [A] [A]-1 單雙層 Rank(A)=4 Rank(A)=4 單三層 Rank(A)=4 單四層 Rank(A)=3 雙三層 Rank(A)=3 雙四層 Rank(A)=3 三四層 Rank(A)=2 由(18)~(23)式，配合核函數退化核表二～五，將其、、、四個未知數和邊界條件表為矩陣形式，即(24)式，其[A]若有反矩陣，則四個未知數即可得解，尤拉樑之撓度、傾角、彎矩、剪力也可分析求得，若[A]無反矩陣，則四個