漁船水上運(yùn)動(dòng)特性研究【畢業(yè)作品】

漁船水上運(yùn)動(dòng)特性研究專(zhuān)業(yè)：船舶與海洋工程目錄 5 5 51.3sesam軟件的簡(jiǎn)介說(shuō)明 61.4課題研究背景、內(nèi)容 7 7 71.5本課題的主要研究?jī)?nèi)容 9 2.2基本方程及邊界條件 2.3線(xiàn)形波浪理論 2.4斯托克斯波理論 2.4.1斯托克斯三階波 2.4.2斯托克斯五階波 2.5.1隨機(jī)波浪波面高度的分布特性 2.5.2隨機(jī)波浪的譜特性 2.6船舶在波浪上的運(yùn)動(dòng)理論基礎(chǔ) 2.7波浪載荷的研究綜述 2.8波浪的長(zhǎng)期預(yù)報(bào) 三船模參數(shù)計(jì)算 3.1.1主尺度及主要參數(shù) 3.1.2裝載工況 3.1.3波浪散布圖的選取 3.2波浪誘導(dǎo)載荷參數(shù) 四全船有限元模型的建立 4.1船體外表面模型（PanelModel） 4.2質(zhì)量模型 4.3水動(dòng)力分析計(jì)算 4．4本章小結(jié) 五數(shù)據(jù)分析結(jié)果 5.1六個(gè)自由度剛體運(yùn)動(dòng)情況 5.1.1在有外龍骨下的六個(gè)自由度的運(yùn)動(dòng)情況 5.1.2在無(wú)外龍骨下的六個(gè)自由度運(yùn)動(dòng)的情況 5.1.3有無(wú)外龍骨的比較 2 參考文獻(xiàn) 3本科畢業(yè)論文摘要用挪威船級(jí)社SESAM軟件HYDRO程序結(jié)合工作海域的波浪散布圖，使用POSTRESP程序輸出不同裝載工況下的不同波浪入射角時(shí)的有無(wú)外龍骨對(duì)整船操縱與搖擺的性能的影響。[關(guān)鍵詞]8米寬魷釣漁船；波浪誘導(dǎo)載荷；外龍骨4本科畢業(yè)論文摘要Waterperformancestudyoffishingboatsf5船舶運(yùn)動(dòng)與波浪載荷的準(zhǔn)確預(yù)報(bào)是進(jìn)行合理的船舶結(jié)構(gòu)設(shè)計(jì)的基礎(chǔ)。在給定了船舶型值的結(jié)構(gòu)強(qiáng)度而又經(jīng)濟(jì)實(shí)用的船體。所以船舶結(jié)構(gòu)生產(chǎn)設(shè)計(jì)水平的進(jìn)步與波浪載荷理論的不斷發(fā)展與完善是密不可分的。船舶在海上經(jīng)歷六個(gè)自由度的剛體運(yùn)動(dòng):縱蕩、橫蕩、垂蕩、縱搖、橫搖和首搖。在這些運(yùn)動(dòng)過(guò)程中，船舶波浪載荷將以垂向、橫向的剪力(彎矩)以及扭矩的形式作用于船體。在船舶的結(jié)構(gòu)分析中，最初的步驟是確定船體的靜水和波浪誘導(dǎo)載荷。船體將在其運(yùn)營(yíng)期中承受這兩種載荷的極端情況，這將導(dǎo)致船體可能的疲勞及破損，這正是我們所關(guān)心的問(wèn)題。船舶的設(shè)計(jì)載荷被定義為船舶運(yùn)營(yíng)期間一定概率水平下的載荷。為了得到這些設(shè)計(jì)載荷，波浪環(huán)境被離散為各種固定的海況。船舶的統(tǒng)計(jì)特性，如船速、浪向角以及波浪干擾力由短期值整合而成。以極值表示的結(jié)果表現(xiàn)為響應(yīng)峰值的分布。對(duì)于疲勞問(wèn)題，此結(jié)果將表現(xiàn)為基于船舶遭受的載荷周期數(shù)目，以及相關(guān)的應(yīng)力幅值的累積損傷。對(duì)于頻域方法，載荷的預(yù)報(bào)涉及響應(yīng)幅度算子的確定。給定波浪頻譜，根據(jù)響應(yīng)幅度算子可以確定船體的運(yùn)動(dòng)響船體在海浪中六個(gè)自由度的運(yùn)動(dòng)均是相對(duì)于其在隨船平動(dòng)坐標(biāo)系中的靜平衡位置而言般水動(dòng)力及波浪載荷的計(jì)算均先以中縱剖面一側(cè)的船體為研究對(duì)象，然后利用對(duì)稱(chēng)關(guān)系得到橫搖和首搖)兩組互不禍合的運(yùn)動(dòng)來(lái)分別討論。6擁有極其豐富的水運(yùn)資源，這為發(fā)展我國(guó)船舶工業(yè)提供了有利條件。合特定航行區(qū)域的波浪散布圖，采用三維勢(shì)流理論對(duì)有外龍骨的漁船進(jìn)行船體運(yùn)動(dòng)響應(yīng)計(jì)合特定航行區(qū)域波浪散布圖，采用三維勢(shì)流理論在頻域內(nèi)對(duì)漁船進(jìn)行隨機(jī)波浪載荷長(zhǎng)期預(yù)報(bào)，并與相應(yīng)規(guī)范計(jì)算值進(jìn)行比較，為合理選取波浪誘導(dǎo)載荷提供依據(jù)。度計(jì)算、疲勞強(qiáng)度計(jì)算方面的研究工作。結(jié)構(gòu)強(qiáng)度計(jì)算方面，隨著計(jì)算機(jī)軟件、硬件能力的結(jié)構(gòu)強(qiáng)度計(jì)算方法的不斷進(jìn)展對(duì)外載荷計(jì)算提出了更高的要求。當(dāng)采用梁理論評(píng)估船體疲勞和極限強(qiáng)度設(shè)計(jì)中使用，而在船體結(jié)構(gòu)直接設(shè)計(jì)計(jì)算中，整船有限元分析使水動(dòng)壓力成船體表面波動(dòng)壓力和波浪誘導(dǎo)的剖面剪力彎矩是在強(qiáng)度計(jì)算中要考慮的最重要的外載荷?；诖?，尋求合理而實(shí)用的船舶運(yùn)動(dòng)、波浪載荷和水動(dòng)壓力預(yù)報(bào)方法對(duì)船舶工程的發(fā)展具有重要意義。本文擬在這方面做一些工作。SESAM是挪威船級(jí)社（DNV）針對(duì)海洋平臺(tái)、船舶等海上結(jié)構(gòu)物開(kāi)發(fā)的一個(gè)大型有限元靜力與動(dòng)力分析程序及顯示計(jì)算結(jié)果的后處理器等四個(gè)部分。其中GENIE模塊用于建模，Hydrod模塊主要進(jìn)行水動(dòng)力計(jì)算,它是基于三維頻率的船舶運(yùn)動(dòng)和波浪載荷預(yù)報(bào)模塊，解決荷可以轉(zhuǎn)換為有限元進(jìn)行結(jié)構(gòu)分析。準(zhǔn)備數(shù)據(jù)的時(shí)間；計(jì)算船舶波浪載荷程序可計(jì)算并可直接傳遞到結(jié)構(gòu)有限元模型上；結(jié)構(gòu)分析程序可解決各種結(jié)構(gòu)在個(gè)種類(lèi)型載荷下的應(yīng)力示計(jì)算結(jié)果，而且提供了用語(yǔ)作最終評(píng)估的程序。該系統(tǒng)在不同程序模塊之間的數(shù)據(jù)自動(dòng)傳7遞是通過(guò)界面文件實(shí)現(xiàn)，使數(shù)據(jù)傳遞非常靈活而有效。的選取應(yīng)根據(jù)不同結(jié)構(gòu)的受力特性分別選取不同的有限元類(lèi)型。由于整船結(jié)構(gòu)有限元模型需工數(shù)據(jù)準(zhǔn)確及計(jì)算機(jī)計(jì)時(shí)。小型遠(yuǎn)洋漁船的穩(wěn)性安全，是漁船設(shè)計(jì)中首先必須考慮的一個(gè)重要問(wèn)題。設(shè)計(jì)人員通常在漁船魚(yú)艙底部設(shè)置固體或液體壓載來(lái)降低漁船空船的重心高度，以達(dá)到漁船穩(wěn)性安全的目了漁船的艙容，而且壓載效果也不一定最佳。通過(guò)在小型遠(yuǎn)洋漁船上設(shè)置箱形外龍骨，并在又可減少漁船在大風(fēng)浪中的橫搖擺幅，改善漁民的生產(chǎn)作業(yè)條件，還能省出本該由壓載物占用的艙容，以裝載更多的燃油或淡水，有效的提高漁船的經(jīng)濟(jì)性能。采用線(xiàn)性隨機(jī)波浪理論：線(xiàn)性隨機(jī)波浪理論假定波浪是有多個(gè)波系疊加而相互不引起干簡(jiǎn)諧波疊加而成的不規(guī)則波系.頻域三維匯源理論：船舶在規(guī)則波中運(yùn)動(dòng)的流體動(dòng)力問(wèn)題的關(guān)鍵是求解流場(chǎng)中的速度勢(shì)，包括入射波的貢獻(xiàn)、物體存在和運(yùn)動(dòng)對(duì)流場(chǎng)擾動(dòng)的貢獻(xiàn)以及相互之間的耦合作用。(1）多普勒效應(yīng)（2）流場(chǎng)速度勢(shì)分解（3）船舶運(yùn)動(dòng)方程。隨著計(jì)算機(jī)軟硬件技術(shù)的發(fā)展，全船有限元分析技術(shù)使船體總強(qiáng)度分析有了革命性的突的實(shí)際變形與應(yīng)力。采用全船有限元分析法和波浪載荷的直接計(jì)算可以進(jìn)一步提高計(jì)算結(jié)果8在漁船上外龍骨的橫剖面形式主要選用矩形或方形的鑄造件或鋼板焊接件（內(nèi)設(shè)生鐵塊）板焊接件（內(nèi)澆混凝土）三種?？v剖面形式主要選用等高式，有些拖網(wǎng)漁船因需要設(shè)置為長(zhǎng)度L，寬度B和高度H，在漁船上外龍骨的長(zhǎng)度L設(shè)計(jì)應(yīng)盡可能地貫通全船，高度H一漁船外龍骨的形式以選用鑄造式橫剖面方、矩形或焊接式橫剖面方、矩形為好。鑄造式外龍骨在一艘長(zhǎng)約50m的漁船上，至少會(huì)有一個(gè)端接頭，端接頭設(shè)計(jì)結(jié)構(gòu)形式必須合理，焊接工藝必須規(guī)范。焊接式外龍骨在每道肋位處其內(nèi)必須設(shè)有橫向加強(qiáng)板右、底板須三面連續(xù)焊。外龍骨內(nèi)若設(shè)置生鐵塊，生鐵塊的形狀必須規(guī)則，排列縫隙宜小；若設(shè)置鐵砂，選用的鐵砂粒直徑宜小為好，以增加鐵砂比重；若設(shè)置混凝土，為方便維修，混凝土中的水泥比例不宜過(guò)大。在設(shè)置上述壓載物時(shí)，壓載物的上表面與船底K行板間須留有3~5mm的空隙，以保證船底K行板與外龍骨側(cè)板間的焊接質(zhì)量。本文主要運(yùn)用了挪威船級(jí)社(DNV)的SESAM(SuperElemerntStructureAnalysisModules)軟件系統(tǒng)對(duì)8m寬魷釣漁船整船建立有限元模型，以外龍骨為變量，對(duì)其進(jìn)行了波浪統(tǒng)計(jì)預(yù)報(bào)，確定設(shè)計(jì)波，計(jì)算船體應(yīng)力分布。主要內(nèi)容有：(1)整理8m寬魷釣漁船的圖紙和資料，分析其結(jié)構(gòu)特點(diǎn)，利用SESAM型導(dǎo)入SESAM軟件系統(tǒng)中。(2)利用挪威船級(jí)社(DNV)的SESAM軟件系統(tǒng)的波浪載荷分析與響應(yīng)計(jì)算子模塊WADAM，計(jì)算全船在一系列規(guī)則波上的流體動(dòng)力載荷，包括動(dòng)壓力在濕表面上的分布及其六個(gè)自由度的運(yùn)動(dòng)特性，形成載荷傳遞函數(shù)。軟件系統(tǒng)的進(jìn)行計(jì)算結(jié)果統(tǒng)計(jì)處理的交互式圖形后處理器POSTRESP,對(duì)計(jì)算結(jié)果進(jìn)行統(tǒng)計(jì)處理和預(yù)報(bào)。通過(guò)長(zhǎng)期預(yù)報(bào)，確定設(shè)計(jì)波。構(gòu)分析的子模塊SESTRA對(duì)全船結(jié)構(gòu)進(jìn)行計(jì)算求解。(5)利用挪威船級(jí)社(DNV)的SESAM軟件系統(tǒng)的數(shù)據(jù)轉(zhuǎn)換處理器PREPOST和有限元分析圖形處理器，對(duì)其結(jié)果以圖形方式顯示出來(lái)，分析全船或局部有限元模型及其運(yùn)動(dòng)分布狀況。9在實(shí)際海況中，船舶運(yùn)動(dòng)受到諸多因素的影響，嚴(yán)重的搖蕩運(yùn)動(dòng)不僅會(huì)影響船舶的運(yùn)計(jì)耐波性能良好的船舶提供可靠的依據(jù)。在船舶耐波性的適居性、實(shí)際使用性、生命力的三特征值、加速度特征值、運(yùn)動(dòng)速度和加速度的極值、甲板淹濕、砰擊和砰擊負(fù)荷、波浪誘導(dǎo)船舶安全性和運(yùn)營(yíng)經(jīng)濟(jì)性具有十分重要的意義。和船舶運(yùn)動(dòng)研究息息相關(guān)的一個(gè)研究領(lǐng)域是船舶在波浪中的載荷研究，船舶在波浪上的運(yùn)動(dòng)理論是研究波浪載荷的理論基礎(chǔ)之一。進(jìn)行船舶結(jié)構(gòu)分析時(shí)，首先要確定作用在船體上結(jié)構(gòu)研究中非常重要的一個(gè)問(wèn)題。對(duì)船舶設(shè)計(jì)者來(lái)說(shuō)，提供一種合理的、接近實(shí)際情況的預(yù)報(bào)工具是非常必要的。在船舶運(yùn)動(dòng)和載荷預(yù)報(bào)領(lǐng)域，雖然基于頻域和概率的線(xiàn)性切片理論已經(jīng)取得了巨大成功，但是仍有首尾砰擊等問(wèn)題便是明顯的例子。由于忽略了非線(xiàn)性載荷，在某些種類(lèi)船型的載荷計(jì)算中，尼，預(yù)測(cè)運(yùn)動(dòng)的幅值也將不準(zhǔn)確。本文根據(jù)船舶在波浪上的運(yùn)動(dòng)理論，運(yùn)用SESAM計(jì)算了船體在波浪中的線(xiàn)性運(yùn)動(dòng)響應(yīng)，（1）計(jì)算8米寬魷釣漁船在不同波浪入射角狀態(tài)下的六個(gè)自由度上運(yùn)動(dòng)的頻率響應(yīng)曲線(xiàn)；（2）計(jì)算8米寬魷釣漁船在規(guī)則波中的垂向波浪剪力和垂向波浪彎矩頻率響應(yīng)曲線(xiàn)，然后選渤海海浪譜，計(jì)算8米寬魷釣漁船在不同海況下的運(yùn)動(dòng)響應(yīng)譜。（3）最后分析計(jì)算得出8米寬魷釣漁船在有無(wú)外龍骨的情況下對(duì)漁船的操縱與搖擺的性能理論基礎(chǔ)。波浪理論也已得到廣泛的研究，主要有線(xiàn)性理論和非線(xiàn)性理論，線(xiàn)性波浪理論(Airy波)是假定波浪振幅足夠小，這樣就可以基本忽略非線(xiàn)性項(xiàng)而得到速度勢(shì)的近似解。駐波理論、流函數(shù)波理論等。以下主要介紹作為基礎(chǔ)的線(xiàn)性波浪理論和本論文重點(diǎn)應(yīng)用的Stokes五階波理論和隨機(jī)波浪譜。如圖2-01表示一在海底平坦、光滑、靜水深度為d的海域向前傳播的波浪。設(shè)波浪兩個(gè)相鄰波峰經(jīng)過(guò)一特定點(diǎn)的時(shí)間間隔T稱(chēng)為波浪的周期。設(shè)c為波浪的傳播速度，ZX波浪周期T-L/Cd▽2Φ2+2應(yīng)滿(mǎn)足Laplace1、海域底部的運(yùn)動(dòng)邊界條件在海底上的流體質(zhì)點(diǎn)不能超過(guò)固體邊界，只能沿著邊界的切線(xiàn)方向運(yùn)動(dòng)，即在z=-d處垂直于固體邊界的法向速度為零，即2、自由表面的運(yùn)動(dòng)邊界條件由表面的水質(zhì)點(diǎn)的垂直于該表面的速度，等于自由表面在該方向的運(yùn)動(dòng)速度。--3、自由表面的動(dòng)力邊界條件數(shù)。因此把貝努利方程式應(yīng)用到自由表面處，得到自由表面的動(dòng)力邊界條件為：(δΦ)2]---(δΦ)2]4、波浪的周期性條件:實(shí)際情況有一定的出入，但這些假定仍然在各種波浪理論中采用。所以對(duì)它們可能產(chǎn)生的從上述可以看出要想精確求解出波浪的速度勢(shì)是非常困難的:一是自由表面條件的非線(xiàn)條件僅在自由表面z=η上滿(mǎn)足，而η本身又是未知量，就使得拉普拉斯方程的求解區(qū)域也是可變的。這就是至今還不能建立起一種波浪理論可普遍適用于各種水深、波高和波長(zhǎng)的海鑒于對(duì)上述各種非線(xiàn)性關(guān)系式求解的困難，必須進(jìn)一步作出一些假定來(lái)簡(jiǎn)化這些關(guān)系式，從而建立起幾種確定性的波浪理論以適用于各自特定的海況條件。簡(jiǎn)單的波動(dòng)。滿(mǎn)足線(xiàn)性波浪理論的波動(dòng)面是水面呈簡(jiǎn)諧形式的起伏運(yùn)動(dòng)。水質(zhì)點(diǎn)的運(yùn)動(dòng)是以平衡位置為圓心的圓周運(yùn)動(dòng)，即以圓頻率負(fù)作簡(jiǎn)諧振動(dòng)。假定波幅或波高相對(duì)于波長(zhǎng)是無(wú)限小，因此可以忽略波動(dòng)自由表面引起的非線(xiàn)性影響，即邊界條件中的乘積項(xiàng)和平方項(xiàng)都可以忽略。并且由于假定波高H足夠小，條件式(2-3)和式(2-4)可近似的在波浪平均位置z=0處滿(mǎn)足。此時(shí)自由表面的邊界條件線(xiàn)性化為2然后通過(guò)分離變量法得出有限水深的線(xiàn)性波速度勢(shì)為─+g2從而得到2負(fù)2這就是線(xiàn)性波的色散關(guān)系。方向的運(yùn)動(dòng)速度如下v=─= 水質(zhì)點(diǎn)在水平及垂直方向的運(yùn)動(dòng)加速度為xz2H-T22HT2位，在海洋工程設(shè)計(jì)中常常采用。這就是所謂的斯托克斯波浪理論。Stokes根據(jù)勢(shì)波理論在推演中考慮了波陡的影響，證明波面不再為簡(jiǎn)單的余弦形式，而是呈波峰較窄而波谷較寬的接近于擺線(xiàn)的形狀，這是和實(shí)際余波的波面頗為相近的。移的近似于圓或橢圓的軌道上運(yùn)動(dòng)。波浪運(yùn)動(dòng)中伴隨有“質(zhì)量的遷移”的這一情況，也是符合于波浪運(yùn)動(dòng)的實(shí)際現(xiàn)象的。所謂的Stokes波是用有限個(gè)簡(jiǎn)單的頻率成比例的余弦波來(lái)逼近具有單一周期的規(guī)則的有限振幅波。為解決自由表面邊界條件的非線(xiàn)性問(wèn)題，假定速度勢(shì)可按某一小參量攝動(dòng)展開(kāi): Φ=εΦ+ε2Φ+因?yàn)椴ǚ邢?，與Φ相關(guān)的其他變量如η等，亦可作相同的小參數(shù)攝動(dòng)展開(kāi):nη都滿(mǎn)足Laplace方程(2-1)及邊界條件;n222n+n困難問(wèn)題。自由表面總是在靜水面附近。將Φ在自由表面z=η處用Taylor級(jí)數(shù)展開(kāi)為+.n+2將上式代入自由表面邊界條件(2-3)及(2-4)，可得)))]2)]2理合并，得2|1)212122(2kd))]于ε的偏微分方程組:一階: 二階:──η1-11=0──η1-11=022由表面邊界條件(2-26)、(2-27)式可以看出，在Stokes高階波勢(shì)的邊界條件中包含根據(jù)以上推導(dǎo)，可得到Stokes三階波的有關(guān)計(jì)算結(jié)果13二、色散關(guān)系2三、波面方程四、水質(zhì)點(diǎn)的運(yùn)動(dòng)速度和加速度|4(3|4(3水質(zhì)點(diǎn)水平速度v c3水質(zhì)點(diǎn)垂直速度v c3水質(zhì)點(diǎn)的水平加速度TT12x=水質(zhì)點(diǎn)的垂直加速度.12.1-5-52233(π)2(d)（L)33.（L)3f6(kd)]f3|Stokes五階波的推導(dǎo)和三階波完全一樣，近些年來(lái)五階波得到越來(lái)越多的關(guān)注，其非線(xiàn)性更加顯著。海洋工程的各類(lèi)規(guī)范經(jīng)常要求用五階波進(jìn)行設(shè)計(jì)。表2.1是Stokes五階波表2-1有關(guān)Stokes五階波計(jì)算公式nn=Σn5nnn22pxxxzz3A5A3A5AA1+λ4B+λ5B5B表達(dá)式系數(shù)中的λ是一個(gè)比值，且λ=ka，a是須對(duì)每一個(gè)波確定的達(dá)式如下:4)642)7)64216426424223)6429)66)86424292642424)6421ijiji定。在求得A,B后，λ及kd可由下式確定:ijij1215)H2gT2應(yīng)用上述Stokes五階波理論的計(jì)算公式，可以用計(jì)算機(jī)編寫(xiě)相應(yīng)的程序，只要輸入水深、波高、周期等基本參數(shù)數(shù)據(jù)，便可得到五階波波長(zhǎng)、波的傳播速度、和波面升高等結(jié)果，并可計(jì)算出波浪水質(zhì)點(diǎn)的速度、加速度和壓力值。計(jì)理論研究海浪現(xiàn)象的理論稱(chēng)為隨機(jī)波浪理論。海底管道立管部分在海洋中所遭遇的海浪一性的變化，因此，海洋中的波浪作為一個(gè)隨機(jī)過(guò)程從長(zhǎng)期而言并不具備平穩(wěn)性。但對(duì)較短的一段時(shí)間來(lái)說(shuō)，可以認(rèn)為波浪是一個(gè)平穩(wěn)的隨機(jī)過(guò)程。此外，觀(guān)測(cè)證實(shí)，波面升高η(t)又歷經(jīng)性保證某一隨機(jī)過(guò)程的一個(gè)具體樣本能代替總體，通過(guò)一樣本能推求出隨機(jī)過(guò)程總體的統(tǒng)計(jì)特性。每一短期海況由表征波浪特性的參數(shù)以及該海況出現(xiàn)的頻率來(lái)描述。常用的波浪參數(shù)為有義波高Hs和平均跨零周期T。Hs定義為所有波浪中波高最大的三分之一波浪的平Z13T定義為波面升高在相鄰兩次以正斜率跨越零均值線(xiàn)之間Z的平均時(shí)間間隔。它是波浪過(guò)程的跨零率的倒數(shù)，故有:1Zf0mm0m2余弦波分量疊加而成的，可以假定其振幅、頻率、相位與方向都是隨機(jī)量，但為處理簡(jiǎn)便，當(dāng)前常僅假定相位是隨機(jī)量。海面上某點(diǎn)的升高為η(t)=Σ偽iiiiiiiiii易證明η(t)的均值為0，即E[η(t)]=0。ii222式中σ2(2(t))表示波面高度的方差。η(2-30)(2-31)構(gòu)由它的各組成波所提供的能量來(lái)體現(xiàn)。海浪譜從數(shù)學(xué)意義上講就是一個(gè)函數(shù)。所謂譜分析對(duì)振幅為a的規(guī)則波，單位波面(單位波長(zhǎng)*單位波寬)的波浪內(nèi)所具有的波能為n12n或12na2可見(jiàn)波能量與波幅平方成正比。因此把d負(fù)可見(jiàn)波能量與波幅平方成正比。因此把d負(fù)范圍內(nèi)的各子波的2迭加起來(lái)，并除以η22η負(fù)顯然函數(shù)S(負(fù))比例于頻率位于間隔η顯然函數(shù)S(負(fù))比例于頻率位于間隔η內(nèi)的各組成波提供的能量。如取η于波頻的分布，故又稱(chēng)頻譜。η海浪波面高度的方差σ2(=η2(t))也比例于單位波面內(nèi)各組成波的總能量，即η=Σn2η1-21-2a=率。這時(shí)，波幅服從Rayleigh分布。由此可計(jì)算波浪過(guò)程的標(biāo)準(zhǔn)差σ、跨零率f、峰值率η0n，以及不規(guī)則系數(shù)a和帶寬系數(shù)ε分別為:0ση0f0112m0n01mmε=2,目前求海浪譜的主要方法有:(1)利用定點(diǎn)觀(guān)測(cè)到的波面記錄η(t)，計(jì)算波面高度的自相關(guān)函數(shù)，然后經(jīng)Fourier變換求得頻譜.(2)由觀(guān)測(cè)資料得出波高與周期的聯(lián)合分布函數(shù)，經(jīng)過(guò)理論推導(dǎo)，得出能量相對(duì)于頻率下面介紹P-M譜:譜有如下的表達(dá)式:24]|式中無(wú)因次常數(shù)a=0.0081，上式所示的波譜僅含一個(gè)參數(shù)即海面上19.5米高處的風(fēng)aU42-2U2將（2-41）代將負(fù)代入（2-42）可得：0-4P-M譜為經(jīng)驗(yàn)譜，由于所依據(jù)的資料比較充分，分析方法比較合理，使用也比較方便，因此在海洋工程和船舶工程中得到了廣泛的應(yīng)用。用力的大小;反過(guò)來(lái)，受力又影響了船體的運(yùn)動(dòng)。因此，二者是禍合的，原則上講，船體受動(dòng)規(guī)律和流體運(yùn)動(dòng)規(guī)律都是可以假定的，因而可以將流體運(yùn)動(dòng)和船體運(yùn)動(dòng)分開(kāi)處理。本文在研究船舶在波浪上的運(yùn)動(dòng)過(guò)程中引入如下假定：1)流體假定:假定流體是沒(méi)有粘性的理想流體，同時(shí)認(rèn)為流體是無(wú)旋有勢(shì)、均勻不可壓2)剛體假定:假定船體是具有六個(gè)自由度的剛體，在外力作用下不變形；3)微幅波假定:假定入射波的波幅A遠(yuǎn)小于波長(zhǎng)兄，由此可以認(rèn)為船體在波浪上的搖蕩運(yùn)動(dòng)也是相對(duì)于平衡位置的小振幅運(yùn)動(dòng)，船體運(yùn)動(dòng)引起的濕表面積變化量遠(yuǎn)小于其靜浮在靜4)等速直線(xiàn)運(yùn)動(dòng)假定:假定波浪引起的船體搖蕩運(yùn)動(dòng)參數(shù)隨時(shí)間的變化要比船體操縱運(yùn)動(dòng)參數(shù)隨時(shí)間的變化快得多。如圖2-1所示，本文定義三個(gè)右手坐標(biāo)系:于描述波浪；2)參考坐標(biāo)系oxyz。原點(diǎn)o位于船體重心，ox軸與船舶航行方向一致，oz軸豎直向3)隨船平動(dòng)坐標(biāo)系o'x'y'z'。原點(diǎn)o'位于未擾動(dòng)的靜水面上，o'x'軸與船舶航行方向一致，o'z'軸豎直向上。在參考坐標(biāo)系中，動(dòng)坐標(biāo)的原點(diǎn)的位置行(,,)稱(chēng)為船體相對(duì)參考為縱蕩，為橫蕩，為垂蕩。oxyz分別繞x,y,z;軸旋轉(zhuǎn)a,β,Y角度和oxy是ZZ0ZUtx'oxxo'OU0圖2-2顯示了船舶在波浪上運(yùn)動(dòng)六種運(yùn)動(dòng)模態(tài)：縱蕩、橫蕩、垂蕩、縱搖、橫搖、艏搖。ZZoxpitchrollyox對(duì)于航行于海上的船舶而言,波浪載荷（即波浪誘導(dǎo)載荷）是所有船舶載荷中最重要的洋波浪的隨機(jī)性、船體結(jié)構(gòu)的復(fù)雜性及航行工況的多樣性，計(jì)算波浪載荷并非易事。一種比較強(qiáng)度。克雷洛夫在1896年發(fā)表的《船舶在波浪上的縱搖新理論及由此運(yùn)動(dòng)而產(chǎn)生的應(yīng)力》及1955年Korvin-Kroukovsky提出了處理船舶搖蕩問(wèn)題的切片理論，后經(jīng)Korvin-Kroukovsky并將最初的迎浪縱向計(jì)算推廣到斜浪5個(gè)自由度運(yùn)動(dòng)的計(jì)算。由此產(chǎn)生了很多類(lèi)似的計(jì)算的計(jì)算方法、水動(dòng)力系數(shù)耦合項(xiàng)及尾端修正項(xiàng)的處理、繞射力的確定等，但所得的數(shù)值結(jié)果及慣性力。在相對(duì)較低的遭遇頻率或高傅氏數(shù)的情況下，線(xiàn)性切片理論的計(jì)算結(jié)果與實(shí)際情致的流體靜力和動(dòng)力的非線(xiàn)性以及大幅度運(yùn)動(dòng)可能發(fā)生的砰擊現(xiàn)象。這些非線(xiàn)性計(jì)算方法大致分為頻率和時(shí)域方法。最具有代表性的頻域非線(xiàn)性方法是基于攝動(dòng)原理建立的二階理論。此理論中二階項(xiàng)源于激勵(lì)波浪的非線(xiàn)性、船舶的非直舷以及非線(xiàn)性的流體動(dòng)力。采用時(shí)域的流體力的計(jì)算，考慮了瞬時(shí)水下剖面的真實(shí)形狀以及水動(dòng)力系數(shù)隨吃水的變化；Meyerhoff和Schlachter通過(guò)模擬不規(guī)則海浪中船舶的航行狀態(tài)求解運(yùn)動(dòng)方程，考慮了波浪沖擊力對(duì)船舶運(yùn)動(dòng)的反饋；Borresen和Tellsgard提出了一種預(yù)報(bào)幅度規(guī)則波中迎浪蕩和縱搖耦合運(yùn)動(dòng)以及垂向載荷的非線(xiàn)性響應(yīng)方法。在時(shí)域中步進(jìn)求解斜浪中船舶運(yùn)動(dòng)微分方程時(shí)存在數(shù)值發(fā)散問(wèn)題。由于船舶在水平面內(nèi)的運(yùn)舵及船體上流體作用力的計(jì)算非常復(fù)雜。目前國(guó)際上通用的方法是建立可計(jì)及操舵效應(yīng)的船此外，針對(duì)經(jīng)典切片理論預(yù)報(bào)高速船的運(yùn)動(dòng)響應(yīng)效果欠佳，充分考慮高速船的高速細(xì)長(zhǎng)理論（二維半理論）應(yīng)運(yùn)而生。在高速細(xì)長(zhǎng)理論中，流場(chǎng)速度勢(shì)的控制方程和物面條件仍是論僅適合高航速下的船舶耐波性的預(yù)報(bào)。時(shí)域問(wèn)題的完善解法是通過(guò)建立初、邊值問(wèn)題來(lái)解決。Finkelstein系統(tǒng)的推導(dǎo)了各種憶項(xiàng)兩部分分別求解，從而將船體的幾何形狀與船舶運(yùn)動(dòng)相分離；Ogilvie總結(jié)了上述理論并推廣到有航速情形。目前應(yīng)用時(shí)域格林函數(shù)求解非線(xiàn)性問(wèn)題已經(jīng)取得了很大的進(jìn)展。林函數(shù)法采用在船體表面布置源的方法確定速度勢(shì)，點(diǎn)源速度勢(shì)－格林函數(shù)滿(mǎn)足除船體表面以外的所有邊界條件，而流場(chǎng)速度勢(shì)可由格林函數(shù)沿船濕表面積分得到。此方法處理無(wú)航速船舶及其他海洋結(jié)構(gòu)物的運(yùn)動(dòng)和載荷問(wèn)題非常成功，但對(duì)于有航速問(wèn)題，由于有航速頻域格速假定，利用無(wú)航速格林函數(shù)并加上關(guān)于航速的修正項(xiàng)來(lái)處理此問(wèn)題。Rankine源法在近些力問(wèn)題。Bishop和Price首先運(yùn)用水彈性力學(xué)方法系統(tǒng)的研究撓性船體在波浪中的動(dòng)力響應(yīng)問(wèn)題。將船體模型簡(jiǎn)化成一根兩端自由的彈性梁，與經(jīng)典的切片理論結(jié)合，采用結(jié)構(gòu)模態(tài)分析法求解船體在波浪中的動(dòng)力響應(yīng)。由于在模態(tài)中引入剛體模態(tài)，因而這種方法可以同時(shí)得到剛體位移和彈性變形。其后許多學(xué)者將水彈性方法與擴(kuò)展的切片理論相結(jié)合，用于非線(xiàn)形成了一種適用于在波浪中或水下運(yùn)動(dòng)的任意形狀彈性體的三維線(xiàn)性水彈性理論。果，對(duì)于弄清載荷的特征及其統(tǒng)計(jì)分布規(guī)律，具有重要的意義。但受種種客觀(guān)及主觀(guān)因素的限制，實(shí)船海上試驗(yàn)結(jié)果不可避免的存在著局限性。為此作為實(shí)船海上試驗(yàn)的補(bǔ)充和完善，知的現(xiàn)象，反映問(wèn)題的力學(xué)機(jī)理和實(shí)質(zhì)，對(duì)理論計(jì)算方法進(jìn)行改進(jìn)和完善。分析的方法得到不規(guī)則波中船舶運(yùn)動(dòng)或波浪載荷的響應(yīng)譜2的響應(yīng)譜。因?yàn)閷?shí)際頻率為遭遇頻率。，所以輸入波浪譜S()應(yīng)變換為Swo根據(jù)統(tǒng)計(jì)時(shí)間的長(zhǎng)短，對(duì)系統(tǒng)響應(yīng)的極值預(yù)報(bào)可以分為短期預(yù)報(bào)和長(zhǎng)期預(yù)報(bào)。短期預(yù)報(bào)是指統(tǒng)計(jì)時(shí)間在半小時(shí)到數(shù)小時(shí)之間。在此期間，假定船舶的裝載、航速、應(yīng)服從Rayleigh分布。該分布只有方差σ2一個(gè)參數(shù)，可由響應(yīng)譜S(負(fù))得到，即y2進(jìn)而可得到船舶運(yùn)動(dòng)與波浪載荷短期預(yù)報(bào)的各種統(tǒng)計(jì)值，包括均值和有義值等。其中，均值為0過(guò)程的組合來(lái)處理。長(zhǎng)期預(yù)報(bào)的概率密度函數(shù)可由很多短期的概率密度函數(shù)以其出現(xiàn)概率為在與確定船體總載荷設(shè)計(jì)值有關(guān)的問(wèn)題中，波浪海況的選取或許是最為重要的。確定該傳遞函數(shù)為基礎(chǔ)導(dǎo)出相應(yīng)的長(zhǎng)期分布。目前使用的方法是在所考慮海域內(nèi)選取一個(gè)實(shí)際波高和波浪周期組合的散射圖譜。這個(gè)散射圖譜也可用來(lái)進(jìn)行極限強(qiáng)度和疲勞強(qiáng)度分析。度分析，而世界范圍的組合波譜已經(jīng)被用于統(tǒng)計(jì)資料為依據(jù)，并結(jié)合對(duì)于波高一些重要修改，以確保最大波峰不在破碎波的極限之外，而且25年一遇的最大的預(yù)計(jì)波高波浪周期與波高同樣重要。早期的波浪氣候資料大都不包含對(duì)目測(cè)波浪周期中誤差的適當(dāng)修使預(yù)報(bào)的長(zhǎng)期波浪載荷誤差達(dá)到7％。而與用波浪散射圖譜進(jìn)行預(yù)報(bào)的長(zhǎng)期波浪載荷相比，表3-18米魷釣漁船的主尺度LOAB△Dd46.81m無(wú)外龍骨六種裝載工況穩(wěn)性參數(shù)見(jiàn)表3-2不同海況和每一海況出現(xiàn)概率組成。每一海況常用有義波高和跨零周期表征。經(jīng)過(guò)幾十年的此采用了渤海海域的全年波浪散布圖。P使得按它計(jì)算出來(lái)的船體應(yīng)力能代表實(shí)際船體航行過(guò)程中的最嚴(yán)重的狀態(tài)。正確的途徑是通過(guò)船體在波浪上受力的長(zhǎng)期統(tǒng)計(jì)預(yù)報(bào)，得到在108概率下預(yù)報(bào)的船體彎矩、剪力等代表性根據(jù)實(shí)踐經(jīng)驗(yàn)以及研究目的的需要，本文選取以下條件進(jìn)行計(jì)算：波浪入射角是波浪傳播方向與船舶航向間的夾角,并由波浪傳播方向按順時(shí)針進(jìn)行度量。當(dāng)波浪入射角在左右舷0°~15°之間時(shí)稱(chēng)為順浪.波浪入射角在左右舷165°~180°之間時(shí)稱(chēng)為頂浪。順浪和頂浪統(tǒng)稱(chēng)為縱向?qū)?縱向?qū)酥饕a(chǎn)生縱向運(yùn)動(dòng),它包括縱搖、縱蕩和垂蕩,其中主要是縱搖和垂蕩。波浪入射角在左右舷75°~105°之間時(shí)稱(chēng)為橫浪。橫浪主要產(chǎn)生橫向運(yùn)動(dòng),它包括橫波浪入射角在左右舷105°~165°之間時(shí)稱(chēng)為首斜浪。尾斜浪和首斜浪既產(chǎn)生縱向運(yùn)動(dòng),也產(chǎn)生橫向運(yùn)動(dòng),如圖3-2所示。首尾浪順G斜斜浪浪圖3-2波浪入射角示意圖波頻ω和波長(zhǎng)λ之間的關(guān)系按λ=2πg(shù)/ω24.1船體外表面模型（PanelMode去，再由點(diǎn)及線(xiàn)，由線(xiàn)建面，初步完成全船的框架。最后應(yīng)當(dāng)對(duì)每一個(gè)面建立相應(yīng)的有限元波浪載荷分析子模塊HydroD中再將先前的Genie中的模型轉(zhuǎn)化成FEM格式并導(dǎo)入到水動(dòng)力程序中進(jìn)行加載分析和計(jì)算。分。對(duì)于水線(xiàn)面以上的部分認(rèn)為不受波浪水動(dòng)壓力。本文采用三維勢(shì)流理論計(jì)算8米寬魷釣漁船迎浪時(shí)運(yùn)動(dòng)響應(yīng)。根據(jù)型線(xiàn)圖建立船體外表正方向指向左舷，z正方向向上。panel模型的建立：建立結(jié)構(gòu)有限元模型要應(yīng)用結(jié)構(gòu)力學(xué)的知識(shí)對(duì)次要構(gòu)件進(jìn)行簡(jiǎn)化處理,而主要構(gòu)件要予以正確表達(dá)。根據(jù)8米寬魷釣漁船母型船的資料，利用型線(xiàn)圖將所需的數(shù)據(jù)讀取，通過(guò)SESAM（1）先從8米寬魷釣漁船的型線(xiàn)圖中讀出個(gè)站位的三維坐標(biāo)值，并將其通過(guò)幾何（Geometry）中的點(diǎn)（Point）功能輸入坐標(biāo)值，再用曲線(xiàn)（curve）功能將他們光順合理的將每個(gè)站位的曲線(xiàn)連接起來(lái)。（2）利用Elements功能建立圖形的有限元網(wǎng)格的劃分。首先，在每個(gè)站位的曲線(xiàn)上對(duì)布置網(wǎng)格劃分節(jié)點(diǎn)，再通過(guò)Mesh功能將各個(gè)站位的網(wǎng)格劃分出來(lái)。其次，利用verify、Equivalence功能來(lái)對(duì)模型進(jìn)行檢查，消除重復(fù)節(jié)點(diǎn)。在處理網(wǎng)格劃分時(shí)往往會(huì)產(chǎn)生重復(fù)節(jié)點(diǎn)，這樣的節(jié)點(diǎn)需要檢查并處理，否則析運(yùn)算。并且，網(wǎng)格劃分好之后，應(yīng)該對(duì)網(wǎng)格的質(zhì)量加以評(píng)估，如果覺(jué)得不甚滿(mǎn)意，則可以（4）定義有限元模型單元的物理特性（Properties）Properties的部分功能就是根據(jù)實(shí)際結(jié)構(gòu)和分析程序創(chuàng)建物理特性，并將其應(yīng)用于對(duì)一項(xiàng)中選擇materials，然后右擊選擇newmaterials一項(xiàng)，要定義其中的各項(xiàng)屬性首先需要命名，然后再是定義Yield，Density，Young，poission，等等其他的相關(guān)數(shù)據(jù)。之后再定義殼單元屬性，厚度取12mm。具體的方法是在propertie擊然后Newthicknesses建立板厚選項(xiàng)，先命名，然后再點(diǎn)Alloweidt，選擇全船，取板圖4-18米寬有龍骨的魷釣漁船panel模型（5）在定義好了8米寬有龍骨的魷釣漁船panel模型之后，要進(jìn)行8米寬無(wú)龍骨的魷釣漁圖4-28米寬無(wú)龍骨的魷釣漁船panel模型量,因此直接計(jì)算得到的結(jié)構(gòu)重量總是小于實(shí)際的空船重量,重心位置也與空船重心不相符。構(gòu)模型直接計(jì)算得到的結(jié)構(gòu)自重與空船重量相符。模型可以與實(shí)船的真實(shí)重量分布基本一致。眾所周知,外載荷、結(jié)構(gòu)響應(yīng)和強(qiáng)度標(biāo)準(zhǔn)是評(píng)定艦船與海洋工程結(jié)構(gòu)物安全的3個(gè)相互關(guān)聯(lián)的重要因素。其中,合理的確定載荷是正確評(píng)定結(jié)構(gòu)安全性的基礎(chǔ)和關(guān)鍵。艦船與海洋程結(jié)構(gòu)所受的波浪載荷,進(jìn)行了長(zhǎng)期的研究,并形成了經(jīng)驗(yàn)公式、數(shù)值預(yù)報(bào)與模型試驗(yàn)3種技術(shù)手段水動(dòng)力分析可獲得海洋結(jié)構(gòu)物在波浪中的運(yùn)動(dòng)、特征剖面上的載荷和濕表面的水動(dòng)壓特征剖面上的載荷響應(yīng)用于結(jié)構(gòu)物的總體強(qiáng)度分析。不同的海洋結(jié)構(gòu)物,所關(guān)注的特征剖面上的載荷響應(yīng)也有所不同。如,關(guān)注“船”形浮體的橫向剖面的剪力、彎矩、扭矩等。濕表面的水動(dòng)壓力響應(yīng)是艙段與全船有限元分析及其疲勞分析中必不可少的載荷參數(shù)。水動(dòng)力分析的目的是研究浮體在海浪中的動(dòng)態(tài)響應(yīng)。海浪是各態(tài)歷經(jīng)的平穩(wěn)隨機(jī)過(guò)程,由無(wú)限多個(gè)頻率不等、方向不同、振幅變化而且相位雜亂的簡(jiǎn)諧微幅波疊加而成。浮體在隨機(jī)不規(guī)則波中的響應(yīng)譜可由單位波幅波中的響應(yīng)(傳遞函數(shù))和海浪譜來(lái)確定。由此可知,浮體在規(guī)則波中響應(yīng)的分析是它在隨機(jī)不規(guī)則海浪中響應(yīng)分析的基礎(chǔ)。水動(dòng)力分析中,通常假定浮體所處的海洋是均勻、不可壓縮、無(wú)粘和無(wú)旋的理想流場(chǎng)。因而,流場(chǎng)的速度勢(shì)滿(mǎn)足拉普拉斯方程,并滿(mǎn)足流場(chǎng)自由表面和底部條件、浮體物面條件以及初始條件與輻射條件。數(shù)學(xué)求解浮體水動(dòng)力的困難在于完全滿(mǎn)足上述條件非線(xiàn)性定解問(wèn)題。為此,人們引入了一定的假定,并作了簡(jiǎn)化,形成了水動(dòng)力分析理論。運(yùn)用SESAM程序系統(tǒng)中HydroD進(jìn)行水動(dòng)力分析。擇，輸入數(shù)據(jù)。如圖4-3、圖4-4、圖4-5所示。查靜水彎矩計(jì)算書(shū)讀取船舯水線(xiàn)處的吃水，再通過(guò)船舯水線(xiàn)處（Z軸方向）的吃水、縱傾和尾傾來(lái)定義裝載工況（Loadcondition）。并將所讀取的數(shù)值輸入。如圖4-8所示。通過(guò)SESAM建立panel模型，定義六種有外龍骨下的不同的工況和擇正確的波浪方向及頻率、水動(dòng)力模型、裝載工況和海洋狀況。輸入8米寬魷釣漁船的垂線(xiàn)間長(zhǎng)，通過(guò)MaximumPrint輸出信息，如圖4-10所示。界面上的窗口就會(huì)顯示出做好的模對(duì)話(huà)框，選定需要計(jì)算的工況，點(diǎn)擊“ExecuteWadam”，運(yùn)行得出波浪載荷預(yù)報(bào)結(jié)果。如復(fù)雜船體結(jié)構(gòu)的模型化技術(shù)是進(jìn)行全船有限元分析的基礎(chǔ)。本章以8米寬魷釣漁船為研最終的計(jì)算結(jié)果能較準(zhǔn)確地反映出實(shí)船的重力和慣性力的真實(shí)分布，為水動(dòng)力分析做了奠定，使得能夠較為準(zhǔn)確的反映出8米寬魷釣漁船在有無(wú)外龍骨下的六中自由度的長(zhǎng)期預(yù)報(bào)。圖5-1到圖5-6為船體在不同波浪入射角作用下對(duì)不同裝載狀態(tài)的六個(gè)自由度剛體運(yùn)和空載結(jié)冰抵達(dá)基地港（工況56）4.00E+03)圖5-1在不同波浪入射角情況下的垂蕩(HEAVE)4.00E+03)單(圖5-2在不同波浪入射角情況下的橫蕩(SWAY)4.00E+00)單(L圖5-3在不同波浪入射角情況下的橫搖(ROLL)4.00E+03)圖5-4在不同波浪入射角情況下的縱蕩(SURGE)4.00E-01)圖5-5在不同波浪入射角情況下的艏搖(YAW)4.50E-014.50E-014.00E-01)圖5-6在不同波浪入射角情況下的縱搖(PITCH)4.00E+03)圖5-7在不同波浪入射角情況下的垂蕩(HEAVE)4.00E+03)單(圖5-8在不同波浪入射角情況下的橫蕩(SWAY)4.00E+00)單(L圖5-9在不同波浪入射角情況下的橫搖(ROLL)4.00E+03)圖5-10在不同波浪入射角情況下的縱蕩(SURGE)4.00E-01)圖5-11在不同波浪入射角情況下的艏搖(YAW)4.50E-014.50E-014.00E-01)圖5-12在不同波浪入射角情況下的縱搖(PITCH)在LC1工況下的六個(gè)自由度比較有無(wú)外龍骨對(duì)船舶操縱搖擺性能的影響。4.00E+03有外龍骨)圖5-13在垂蕩（heave）運(yùn)動(dòng)下的有無(wú)外龍骨比較4.50E-014.50E-014.00E-01有外龍骨)單(p圖5-14在縱搖（pitch）運(yùn)動(dòng)下的有無(wú)外龍骨比較(4.00E+00有外龍骨)單圖5-15在橫搖（roll）運(yùn)動(dòng)下的有無(wú)外龍骨比較4.00E+03有外龍骨) 圖5-16在縱蕩（surge）運(yùn)動(dòng)下的有無(wú)外龍骨比較4.00E+03(有外龍骨)單4.00E-01有外龍骨y)圖5-18在艏搖（yaw）運(yùn)動(dòng)下的有無(wú)外龍骨比較為90°處達(dá)到最大值，之后再逐漸的減少。而縱蕩(SURGE)、縱搖(PITCH)這兩個(gè)自由度下的不同裝載工況下船體運(yùn)動(dòng)會(huì)隨著波浪入射角的增加在波浪入射角為90°處達(dá)到最小值，(2)無(wú)論在有無(wú)外龍骨的情況下，在垂蕩(HEAVE)、橫蕩(SWAY)、縱蕩(SURGE)、縱搖(PITCH)四種自由度下，裝載工況的不同所產(chǎn)生的運(yùn)動(dòng)響應(yīng)變化不大。而對(duì)于橫搖(ROLL)、艏搖（YAW）裝載工況變化了，所產(chǎn)生的運(yùn)動(dòng)響應(yīng)變化非常大。（3）當(dāng)在固定裝載工況（LC1）下，對(duì)于垂蕩(HEAVE)、橫蕩(SWAY)、縱蕩(SURGE)、縱搖(PITCH)四種自由度，無(wú)論有無(wú)外龍骨，其運(yùn)動(dòng)響應(yīng)變化均不大。（4）當(dāng)在固定裝載工況（LC1）下，對(duì)于橫搖(ROLL)、艏搖（YAW沒(méi)有外龍骨下所產(chǎn)生的運(yùn)動(dòng)響應(yīng)要遠(yuǎn)遠(yuǎn)大于有外龍骨所產(chǎn)生的運(yùn)動(dòng)響應(yīng)。在SESAM程序中用POSTRESP子模塊進(jìn)行波浪載荷長(zhǎng)期預(yù)報(bào)。在建立panel模型時(shí)，對(duì)結(jié)構(gòu)時(shí)進(jìn)行適當(dāng)?shù)暮?jiǎn)化。本次研究的主要是船體在波浪運(yùn)動(dòng)中產(chǎn)生的波浪剪力和波浪彎矩對(duì)船體結(jié)構(gòu)強(qiáng)度的影響，因此正確計(jì)算和分析船體的波浪誘導(dǎo)載荷預(yù)報(bào)是關(guān)鍵。作為船體結(jié)構(gòu)強(qiáng)度的衡準(zhǔn)，需要找到船舶服務(wù)周期內(nèi)(20—25年)所受到的一次最大波載荷的長(zhǎng)期預(yù)報(bào)來(lái)尋求，研究表明，這相當(dāng)于108超越概率水平的波浪載荷極值。壓載狀態(tài)下船體波浪誘導(dǎo)垂直波浪彎矩、垂直波浪剪力在這截面上的傳遞函數(shù)，得到了相應(yīng)本研究以8米魷釣漁船為研究對(duì)象，應(yīng)用挪威船級(jí)社(DNV)的SESAM軟件系統(tǒng)，進(jìn)行了全船有限元水動(dòng)力分析。在完成了建模以及波浪載荷的長(zhǎng)期預(yù)報(bào)后，應(yīng)用設(shè)計(jì)波法確定了波來(lái)對(duì)比不同入射角情況下的波浪載荷預(yù)報(bào)。在有龍骨情況下的六個(gè)自由度不同工況運(yùn)動(dòng)特性表1：不同裝載工況下橫搖的運(yùn)動(dòng)響應(yīng)ROLL00.00E+000.00E+002.71E-038.95E-040.00E+000.00E+003.04E+002.46E+008.60E-015.91E-015.65E-016.05E+002.29E+004.90E+008.93E+003.44E+007.25E+002.58E+004.60E+009.36E+003.30E+002.27E+002.16E+005.54E+003.78E+002.59E+002.47E+005.93E+003.96E+002.72E+002.58E+005.58E+003.81E+002.63E+002.50E+004.65E+009.46E+003.35E+002.32E+002.21E+009.04E+003.48E+007.34E+002.62E+006.12E+002.31E+004.95E+003.07E+002.48E+008.73E-016.01E-015.74E-010.00E+000.00E+002.85E-030.00E+000.00E+000.00E+00表2：不同裝載工況下縱搖的運(yùn)動(dòng)響應(yīng)PITCHPITCH3.84E-014.07E-013.67E-013.61E-012.98E-018.42E-0304.08E-014.04E-013.90E-013.63E-013.11E-012.00E-012.22E-024.10E-014.05E-013.89E-013.58E-013.01E-018.81E-034.10E-014.05E-013.89E-013.58E-013.02E-018.64E-034.03E-013.99E-013.86E-013.60E-013.09E-012.01E-012.59E-022.02E-013.82E-024.12E-013.43E-013.04E-013.79E-013.92E-012.98E-012.79E-012.94E-013.08E-013.07E-013.07E-013.58E-013.34E-013.52E-013.63E-013.61E-013.61E-013.87E-013.62E-013.82E-013.93E-013.91E-013.92E-014.01E-013.76E-013.96E-014.08E-014.06E-014.07E-014.06E-013.80E-014.01E-014.13E-014.11E-014.11E-01表3：不同裝載工況下艏搖的運(yùn)動(dòng)響應(yīng)YAWYAW00.00E+000.00E+000.00E+000.00E+000.00E+000.00E+009.63E-027.34E-026.70E-026.58E-023.79E-013.25E-015.58E-012.57E-014.72E-017.17E-013.14E-015.97E-012.38E-018.34E-013.52E-016.87E-012.55E-018.91E-013.70E-017.31E-012.61E-018.71E-013.67E-017.18E-012.71E-012.03E-017.73E-013.36E-016.43E-012.60E-012.10E-012.02E-016.13E-012.76E-015.16E-012.18E-014.18E-013.56E-012.11E-017.96E-027.10E-026.93E-020.00E+000.00E+000.00E+000.00E+000.00E+000.00E+00表4：不同裝載工況下橫蕩的運(yùn)動(dòng)響應(yīng)5.49E-013.20E+034.68E+0303.03E+034.42E+033.21E+034.70E+033.22E+034.70E+034.48E-013.07E+034.49E+033.06E+034.48E+035.59E+035.70E+035.69E+035.96E+035.96E+035.96E+036.42E+036.56E+036.51E+036.87E+036.83E+036.84E+036.74E+036.87E+036.80E+037.21E+037.15E+037.16E+036.45E+036.54E+036.47E+036.89E+036.84E+036.85E+035.64E+035.66E+035.62E+036.01E+035.98E+035.97E+034.46E+034.44E+034.42E+034.72E+034.72E+034.71E+033.06E+033.03E+033.02E+033.22E+033.23E+033.22E+033.33E-012.63E-01表5：不同裝載工況下縱蕩的運(yùn)動(dòng)響應(yīng)SURGESURGE06.24E+036.15E+036.22E+036.38E+036.44E+036.43E+036.07E+035.99E+036.05E+036.21E+036.27E+036.26E+035.56E+035.49E+035.55E+035.69E+035.74E+035.73E+034.68E+034.62E+034.67E+034.79E+034.83E+034.82E+033.43E+033.39E+033.42E+033.50E+033.53E+033.53E+037.38E+018.11E+012.43E+013.41E+033.36E+033.40E+033.50E+033.54E+033.53E+034.66E+034.60E+034.65E+034.79E+034.83E+034.83E+035.54E+035.46E+035.52E+035.69E+035.74E+035.73E+036.05E+035.96E+036.03E+036.21E+036.27E+036.26E+036.22E+036.12E+036.20E+036.38E+036.44E+036.43E+03表6：不同裝載工況下垂蕩的運(yùn)動(dòng)響應(yīng)HEAVEHEAVE6.52E+0306.46E+036.45E+036.38E+036.59E+036.59E+036.53E+036.45E+036.51E+036.59E+036.65E+036.65E+036.72E+036.65E+036.70E+036.76E+036.83E+036.83E+037.01E+036.97E+037.00E+037.04E+037.10E+037.10E+037.38E+037.36E+037.38E+037.38E+037.43E+037.44E+037.75E+037.76E+037.76E+037.71E+037.76E+037.76E+037.97E+038.01E+037.97E+037.88E+037.91E+037.92E+037.83E+037.87E+037.84E+037.73E+037.76E+037.77E+037.46E+037.47E+037.46E+037.39E+037.44E+037.44E+037.07E+037.05E+037.07E+037.05E+037.10E+037.10E+036.76E+036.72E+036.76E+036.77E+036.83E+036.83E+036.57E+036.52E+036.56E+036.59E+036.65E+036.65E+036.50E+036.45E+036.49E+036.53E+036.59E+036.59E+03在無(wú)龍骨情況下的六個(gè)自由度不同工況運(yùn)動(dòng)特性表7：不同裝載工況下橫搖的運(yùn)動(dòng)響應(yīng)ROLLROLL0.00E+002.14E+004.26E+006.33E+008.22E+009.63E+009.71E+008.34E+006.44E+004.33E+002.17E+0000.00E+009.73E-012.00E+003.04E+003.96E+004.61E+004.86E+004.67E+004.05E+003.11E+002.04E+009.94E-010.00E+008.58E-012.53E+003.26E+003.79E+004.00E+003.83E+003.31E+002.57E+008.70E-010.00E+006.07E-012.29E+002.65E+002.79E+002.67E+002.33E+006.16E-010.00E+006.58E-012.04E+002.66E+003.09E+003.26E+003.13E+002.72E+002.09E+006.71E-010.00E+007.21E-012.25E+002.94E+003.42E+003.61E+003.47E+003.01E+002.31E+007.36E-010.00E+000.00E+000.00E+000.00E+000.00E+000.00E+000.00E+000.00E+000.00E+000.00E+000.00E+000.00E+00表8：不同裝載工況下縱搖的運(yùn)動(dòng)響應(yīng)PITCH04.16E-014.06E-014.13E-014.10E-014.07E-014.08E-014.11E-014.01E-014.08E-014.05E-014.02E-014.03E-013.97E-013.88E-013.95E-013.89E-013.86E-013.87E-013.68E-013.62E-013.66E-013.58E-013.55E-013.55E-013.13E-013.10E-013.13E-013.02E-012.98E-012.99E-012.01E-012.48E-028.89E-039.80E-033.04E-012.95E-013.01E-013.07E-013.06E-013.07E-013.65E-013.55E-013.62E-013.61E-013.59E-013.59E-013.95E-013.84E-013.91E-013.91E-013.89E-013.89E-014.10E-013.98E-014.06E-014.06E-014.04E-014.04E-014.15E-014.03E-014.11E-014.11E-014.08E-014.09E-01表9：不同裝載工況下艏搖的運(yùn)動(dòng)響應(yīng)YAW00.00E+000.00E+000.00E+000.00E+000.00E+000.00E+006.64E-027.39E-026.96E-026.82E-027.83E-022.78E-014.12E-012.25E-015.33E-012.40E-012.15E-012.24E-012.83E-016.25E-012.59E-012.25E-012.44E-013.21E-016.70E-012.65E-012.28E-012.54E-013.41E-012.01E-016.53E-012.70E-012.43E-012.65E-013.48E-012.04E-015.75E-012.56E-012.40E-012.52E-013.22E-014.52E-012.13E-012.04E-012.