資源描述:
《浮體水動力邊界元分析中線性方程組迭代求解方法研究》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。
1、哈爾濱工程大學(xué)碩士學(xué)位論文浮體水動力邊界元分析中線性方程組迭代求解方法研究姓名:蘇焱申請學(xué)位級別:碩士專業(yè):流體力學(xué)指導(dǎo)教師:陳曉波201106浮體水動力學(xué)邊界元分析中線性方程組迭代求解方法研究摘要隨著海洋運(yùn)輸業(yè)的不斷發(fā)展,船舶與海洋工程結(jié)構(gòu)物尺寸越來越大,結(jié)構(gòu)越來越復(fù)雜,而且出現(xiàn)了很多雙體以及多體問題?;谶吔缭椒ǎ糜嬎銠C(jī)對結(jié)構(gòu)物進(jìn)行輻射繞射計算過程中線性方程組的求解時間隨著物體濕表面面元數(shù)目增加而不斷增加。對于物面面元數(shù)數(shù)以萬計的大型離岸結(jié)構(gòu)物,傳統(tǒng)的直接方法需要花費(fèi)數(shù)天時間求解該線性方程組,這對于以快速立足的
2、邊界元方法是不能接受的。為了提高線性方程組求解效率,本文將對其結(jié)構(gòu)特點(diǎn)進(jìn)行研究并采用合適的迭代方法求解。線性方程組系數(shù)矩陣的結(jié)構(gòu)特性對于迭代求解方法的選擇至關(guān)重要,本文將從條件數(shù)的角度對浮體水動力學(xué)邊界元分析中產(chǎn)生的線性方程組系數(shù)矩陣進(jìn)行研究,分析其在不同計算模型下性態(tài)。間接方法主要是指迭代方法,通過對各種迭代方法以及本領(lǐng)域線性方程組結(jié)構(gòu)特點(diǎn)的研究,在理論上確定適合的迭代求解方法。編寫迭代法求解線性方程組求解程序并與直接法求解結(jié)果對比驗證;為了加速迭代求解方法的收斂,我們采用不完全LU分解預(yù)處理技巧降低求解過程中的迭代次
3、數(shù),并通過對不同形狀離岸結(jié)構(gòu)進(jìn)行計算確定預(yù)處理技巧的有效性。對于不同結(jié)構(gòu)的線性方程組,迭代方法并不一定快于直接法,通過計算并對比兩種方法的浮點(diǎn)運(yùn)算量獲得包含兩種求解程序的優(yōu)化算法。首先本文簡要介紹浮體在波浪中運(yùn)動的勢流理論,然后引出論文所要研究的輻射繞射線性方程組并給出其各項含義。編寫矩陣條件數(shù)計算程序分析三維浮體無航速頻域水動力分析中線性方程組系數(shù)矩陣性態(tài),結(jié)果表明該矩陣條件數(shù)隨著入射波頻率升高和物體濕表面網(wǎng)格劃分?jǐn)?shù)目增加而增大;通過在物體內(nèi)部自由面劃分網(wǎng)格方法消除非規(guī)則頻率,可以明顯緩解非規(guī)則頻率處條件數(shù)突變,但是其
4、它頻率處條件數(shù)將會普遍增大。對于給定線性方程組,編寫基于廣義最小殘量迭代求解方法程序,結(jié)果表明求解效率明顯高于基于LU分解的直接求解方法;在上述迭代程序基礎(chǔ)上,采用不完全LU分解作為預(yù)處理可以顯著緩解高頻階段迭代次數(shù)升高問題,而且對于含有較大內(nèi)部自由面的物體,需要加大預(yù)處理中不完全LU分解的程度才可以使迭代次數(shù)降低到合適的水平;通過對比基于LU分解的直接法和廣義最小殘量的迭代法浮點(diǎn)運(yùn)算量,我們獲得一個優(yōu)化判斷公式,對于不同的線性方程組,該公式可以在兩種方法之間選取運(yùn)算量最小的方法進(jìn)行求解。關(guān)鍵詞:勢流理論;頻域:GMRE
5、S方法;線性方程組;條件數(shù)浮體水動力學(xué)邊界元分析中線性方程組迭代求解方法研究AbstractWiththedevelopmentofmarinetransportation,thedimensionsofoffshorestructuresbecomelargerandthestructuresaremoreandmorecomplicated.Andalsotherehavebeenalotofcatamaranandmulti-bodyproblem.Duetotheincreasingnumberofpanels,
6、itwillcostUSmoreandmoreCPUtimetosolvethelinearsystemsarisingfromthediffraction-radiationanalysisoffloatingbodieswithboundaryelementmethods.Asforlargeoffshorestructures,itwillconsumeUSseveraldaystosolvethelinearsystemsusingtraditionaldirectmethod,whichisnotaccepta
7、bleforengineeringapplication.Inordertoimprovetheeffectiveness,wewillanalysisthefeaturesoflinearsystemsandselectthesuitablesolver.Thestructuresofcoefficientmatricesinlinearsystemsplayaveryimportantroleinchoosingtheappropriatemethod.Thispaperwilldosomeresearchesont
8、hebehaviorofthematricesarisingfromtheboundarymethodofdifferentmodelsfromtheviewpointofconditionnumber.Throughtheresearchesondifferentiterativemethods,謝也thechar