金建海，葉永林，田 超，鄒明松，冷文浩

（1江南大學(xué) 物聯(lián)網(wǎng)學(xué)院，江蘇 無(wú)錫214122；2中國(guó)船舶科學(xué)研究中心，江蘇 無(wú)錫 214082；3中國(guó)船級(jí)社，北京 100007）

三維水彈性集成軟件的設(shè)計(jì)與實(shí)現(xiàn)

金建海1,2，葉永林2，田 超3，鄒明松2，冷文浩2

（1江南大學(xué) 物聯(lián)網(wǎng)學(xué)院，江蘇 無(wú)錫214122；2中國(guó)船舶科學(xué)研究中心，江蘇 無(wú)錫 214082；3中國(guó)船級(jí)社，北京 100007）

為了滿足三維水彈性軟件商用化的應(yīng)用要求，文章通過(guò)對(duì)已有三維水彈性程序的使用經(jīng)驗(yàn)的分析，設(shè)計(jì)并實(shí)現(xiàn)了其前后處理功能模塊，形成統(tǒng)一的商用化集成軟件。并以某18萬(wàn)噸的散貨船為例，利用該商用集成軟件進(jìn)行分析處理，用以進(jìn)一步展示完成的集成軟件的功能。最后給出了軟件進(jìn)一步發(fā)展的設(shè)想。

三維水彈性；前后處理；模態(tài)篩選；主坐標(biāo)響應(yīng)

1 引 言

現(xiàn)代船舶設(shè)計(jì)制造技術(shù)發(fā)展的歷史，也是船舶力學(xué)發(fā)展的歷史。船舶力學(xué)的主要研究?jī)?nèi)容可以概括為兩大學(xué)科：船舶結(jié)構(gòu)力學(xué)和船舶水動(dòng)力學(xué)。到上世紀(jì)70年代末，兩者的交叉與融合促成了一門新的學(xué)科分支—船舶水彈性力學(xué)（Hydroelasticity）的建立。船舶水彈性力學(xué)是從動(dòng)力學(xué)的觀點(diǎn)出發(fā)，將柔性船體與周圍流場(chǎng)作為一個(gè)相互作用的整體系統(tǒng)進(jìn)行描述和分析，并在此基礎(chǔ)上預(yù)報(bào)船體結(jié)構(gòu)的動(dòng)力學(xué)特性，諸如：船體運(yùn)動(dòng)、波浪載荷、結(jié)構(gòu)變形、應(yīng)力應(yīng)變以及疲勞性能等，以更合理地預(yù)報(bào)和評(píng)估船舶及海洋工程結(jié)構(gòu)的運(yùn)動(dòng)特性、可靠性和安全性[1-2]。

三維水彈性理論從1984年建立以來(lái)，發(fā)展了考慮航速等因素的數(shù)值方法[3]。1997年提出了二階非線性水彈性理論，考慮了由船體大型剛體運(yùn)動(dòng)和瞬態(tài)濕表面引起的非線性水動(dòng)力，其后完善了相應(yīng)的數(shù)值分析方法，并應(yīng)用于SWATH船[4]及大型散貨船[5]等。

本文的工作主要介紹在中國(guó)船舶科學(xué)研究中心現(xiàn)有的三維線性頻域水彈性分析程序“THAFTS”基礎(chǔ)上，增加其前后處理功能模塊，形成新的工程實(shí)用化的三維水彈性集成軟件，并利用該集成軟件，對(duì)某18萬(wàn)噸級(jí)的散貨船進(jìn)行主坐標(biāo)響應(yīng)計(jì)算。在文章最后部分，給出了軟件進(jìn)一步發(fā)展的展望。

2 三維線性水彈性分析基本原理

上世紀(jì)70年代中期，隨著大型計(jì)算機(jī)的出現(xiàn)，各種用于分析大型海洋浮動(dòng)結(jié)構(gòu)耐波性問(wèn)題的三維水動(dòng)力學(xué)方法相繼發(fā)展起來(lái)，Wu[6]開(kāi)創(chuàng)性地將三維適航性理論與三維結(jié)構(gòu)動(dòng)力學(xué)理論相結(jié)合，提出了廣義流固界面條件，發(fā)展了一個(gè)適用于分析波浪中任意三維可變形體承受內(nèi)、外激勵(lì)時(shí)動(dòng)響應(yīng)性能的三維水彈性理論。

簡(jiǎn)化的三維線性水彈性力學(xué)理論認(rèn)為船體結(jié)構(gòu)周圍的流體為理想、不可壓縮流體，波浪幅值為小量，此時(shí)可得航行船體的頻域線性水彈性運(yùn)動(dòng)方程為：

3 三維線性水彈性程序簡(jiǎn)介

現(xiàn)有的三維水彈性系列程序都是以吳有生院士的三維線性頻域水彈性力學(xué)程序?yàn)榛A(chǔ)發(fā)展起來(lái)的[1,3,7]。它利用三維水彈性理論，分析任意形狀海洋浮體結(jié)構(gòu)在各種航行條件下，浮體結(jié)構(gòu)在波浪力、系泊力和任意強(qiáng)迫外激勵(lì)力的作用下的剛體運(yùn)動(dòng)和彈性變形等響應(yīng)。目前三維線性頻域水彈性程序包“THAFTS”由三個(gè)模塊程序組成。

其中第一個(gè)程序FLXBD計(jì)算彈性結(jié)構(gòu)濕表面幾何特性以及各階模態(tài)的濕表面邊界條件和廣義流體恢復(fù)力系數(shù)；第二個(gè)程序HYCOF計(jì)算彈性浮體結(jié)構(gòu)周圍流場(chǎng)速度勢(shì)，計(jì)算結(jié)構(gòu)在單位幅值規(guī)則波中剛體運(yùn)動(dòng)和彈性變形模態(tài)的廣義水動(dòng)力系數(shù)和廣義波浪激勵(lì)力。第三個(gè)程序HYELAS求解各階模態(tài)的主坐標(biāo)響應(yīng)。

4 水彈性集成軟件功能模塊設(shè)計(jì)與實(shí)現(xiàn)

為實(shí)現(xiàn)商業(yè)化的水彈性集成軟件，需要增加原有“THAFTS”程序的前處理模塊，以方便地自動(dòng)生成輸入數(shù)據(jù)文件，并自動(dòng)啟動(dòng)程序進(jìn)行計(jì)算；增加后處理模塊，以直觀地顯示計(jì)算結(jié)果。

下文將在理解輸入輸出數(shù)據(jù)文件內(nèi)容的基礎(chǔ)上，結(jié)合水彈性分析經(jīng)驗(yàn)，設(shè)計(jì)并實(shí)現(xiàn)前后處理模塊，并與原有的“THAFTS”程序一起，形成統(tǒng)一的水彈性集成軟件。

4.1 水彈性程序輸入輸出數(shù)據(jù)文件內(nèi)容

水彈性程序“THAFTS”的輸入數(shù)據(jù)文件主要包括flxbd.in、hycof.in以及hyelas.in；輸出文件比較多，但主要是主坐標(biāo)響應(yīng)計(jì)算結(jié)果的輸出。

4.1.1 flxbd.in輸入數(shù)據(jù)文件內(nèi)容

FLXBD程序的數(shù)據(jù)輸入文件，其輸入內(nèi)容分類后可以分成以下6個(gè)部分：

（1）船體干結(jié)構(gòu)有限元模型信息

這類數(shù)據(jù)主要表征原干結(jié)構(gòu)有限元模型的Mesh信息，主要有：

單元總數(shù)、節(jié)點(diǎn)總數(shù)和各個(gè)節(jié)點(diǎn)的ID和坐標(biāo)；

梁?jiǎn)卧倲?shù)、組成梁?jiǎn)卧墓?jié)點(diǎn)ID；

板單元總數(shù)、組成板單元的節(jié)點(diǎn)ID。

干結(jié)構(gòu)坐標(biāo)系與平衡坐標(biāo)系之間的關(guān)系（原點(diǎn)偏移、坐標(biāo)軸夾角余弦）。其中平衡坐標(biāo)系示意圖如下：

（2）船體干結(jié)構(gòu)各階模態(tài)的模態(tài)信息

該類數(shù)據(jù)主要有彈性模態(tài)數(shù)、各階模態(tài)的固有頻率、廣義質(zhì)量、廣義阻尼、振型放大系數(shù)、左右對(duì)稱性，以及各階彈性模態(tài)下的節(jié)點(diǎn)振型等。

（3）船體濕表面網(wǎng)格信息

這類數(shù)據(jù)主要表征原船體濕表面網(wǎng)格的Mesh信息，包括濕面元節(jié)點(diǎn)數(shù)、單元數(shù)，以及節(jié)點(diǎn)坐標(biāo)，組成各個(gè)濕面元的節(jié)點(diǎn)ID等。

（4）干濕單元組對(duì)應(yīng)關(guān)系

FLXBD程序計(jì)算彈性結(jié)構(gòu)濕表面幾何特性以及各階模態(tài)的濕表面邊界條件，需要在輸入文件中指定干單元的分組情況、濕單元的分組情況，以及每組干單元對(duì)應(yīng)哪幾組濕單元。

（5）船體特征參數(shù)

船體特征參數(shù)包括以下數(shù)據(jù)內(nèi)容項(xiàng)：對(duì)稱特性參數(shù)、船長(zhǎng)、艏柱至重心的水平距離、排水體積、重心距水線面的距離、橫搖慣性矩、縱搖慣性矩、搖艏慣性矩以及橫搖/艏搖慣性矩。

（6）計(jì)算控制參數(shù)

控制計(jì)算過(guò)程的相關(guān)參數(shù)，包括是否計(jì)入剛體位移參數(shù)等，這部分的控制參數(shù)可以全部采用默認(rèn)值固定。

4.1.2 hycof.in輸入數(shù)據(jù)文件內(nèi)容

HYCOF程序的數(shù)據(jù)輸入文件，可分為海況參數(shù)和求解控制參數(shù)。

其中海況參數(shù)主要有波幅、波角、速度、波浪頻率等。

求解控制參數(shù)主要有求解方法選擇控制參數(shù)、壓力輸出選擇控制參數(shù)等。

4.1.3 hyelas.in輸入數(shù)據(jù)文件內(nèi)容

HYELAS程序的數(shù)據(jù)輸入文件，內(nèi)容主要包括船體干結(jié)構(gòu)各階模態(tài)的模態(tài)信息、船體特征參數(shù)和求解控制參數(shù)。其中求解控制參數(shù)需要考慮的主要有：

要求輸出廣義附加質(zhì)量、廣義附加阻尼和廣義波浪激勵(lì)力的種類數(shù)；

要求輸出的廣義附加質(zhì)量、廣義附加阻尼或者廣義波浪激勵(lì)力對(duì)應(yīng)的角標(biāo)號(hào)。而其他內(nèi)容在前面兩個(gè)數(shù)據(jù)輸入文件中都已出現(xiàn)過(guò)。

4.1.4 ResponseTC2.dat輸出數(shù)據(jù)文件內(nèi)容

該輸出文件存儲(chǔ)的是計(jì)算頻率點(diǎn)的響應(yīng)結(jié)果，包括各階模態(tài)，各個(gè)頻率下的主坐標(biāo)響應(yīng)的幅值和相位。

4.2 水彈性集成軟件功能模塊設(shè)計(jì)

通過(guò)對(duì)THAFTS三個(gè)模塊程序輸入數(shù)據(jù)文件內(nèi)容的分析，以及工程中的使用經(jīng)驗(yàn)，三維水彈性集成軟件可設(shè)計(jì)成如圖2所示的功能模塊圖：

（1） 三維建模

使用目前成熟的商用軟件，或者參數(shù)化建模等方式，建立船體干結(jié)構(gòu)三維有限元模型，以及用于水動(dòng)力分析的船體濕表面網(wǎng)格模型[8]。并通過(guò)開(kāi)發(fā)接口軟件的方式，將干結(jié)構(gòu)模型和濕面元模型都導(dǎo)入到集成軟件中。

（2）干結(jié)構(gòu)模態(tài)分析

使用Nastran等成熟的商用CAE求解器，或者其他自研方法，依靠其有限元模型，計(jì)算船體干結(jié)構(gòu)的振動(dòng)模態(tài)。并讀取模態(tài)計(jì)算結(jié)果文件，包括各階模態(tài)的固有頻率、節(jié)點(diǎn)振型等。

（3） 模態(tài)篩選

模態(tài)分析之后，要選擇進(jìn)行水彈性計(jì)算的模態(tài)。主要是篩選總模態(tài)，以及熱點(diǎn)區(qū)域局部模態(tài)。前者通過(guò)強(qiáng)結(jié)構(gòu)節(jié)點(diǎn)變形提取，后者通過(guò)局部區(qū)域節(jié)點(diǎn)位移提取。軟件應(yīng)支持自動(dòng)/手動(dòng)兩種方式。

自動(dòng)方式主要包括“關(guān)鍵節(jié)點(diǎn)位移法”和“局部模態(tài)提取法”。

關(guān)鍵節(jié)點(diǎn)位移法是在強(qiáng)框架上選擇n個(gè)關(guān)鍵節(jié)點(diǎn)，并指定變形值，如果某階模態(tài)下這些節(jié)點(diǎn)的變形值大于指定的變形值，則認(rèn)為該階模態(tài)符合要求，自動(dòng)將其挑選出來(lái)。

局部模態(tài)提取法就是剔除其中的局部模態(tài)，保留剩下的整體模態(tài)，以用于水彈性后續(xù)的計(jì)算。其中局部模態(tài)是指存在部分節(jié)點(diǎn)的振型特別突出的情況的模態(tài)，否則視為整體模態(tài)。其算法描述如下：

設(shè)節(jié)點(diǎn)坐標(biāo)為（X k（0,0），Y k（0,0），Z k（0,0））第 k（0,0）個(gè)節(jié)點(diǎn)，該節(jié)點(diǎn)周圍第一圈節(jié)點(diǎn)為 k（1,j），有 M 個(gè)；第二圈節(jié)點(diǎn)為 k（2,j），有 N 個(gè)，第三圈節(jié)點(diǎn)為 k（3,j），有 T個(gè)。節(jié)點(diǎn)模態(tài)歸一化振型為A k（0,0），篩選規(guī)則可采用下式

式中，λ為靈敏度系數(shù)，可設(shè)置，默認(rèn)設(shè)為5。滿足上述不等式，則表明該模態(tài)為局部模態(tài)，可剔除。

手動(dòng)方式設(shè)計(jì)成讓用戶可以在軟件界面查看任一階模態(tài)的節(jié)點(diǎn)振型云圖，由用戶自己根據(jù)變形云圖來(lái)判斷是否要選擇該模態(tài)。手動(dòng)方式可以單獨(dú)使用，也可以在自動(dòng)方式的基礎(chǔ)上繼續(xù)。

（4）設(shè)定模態(tài)參數(shù)

模態(tài)分析之后，各相關(guān)模態(tài)的固有頻率、廣義質(zhì)量都已經(jīng)具備，但還需要設(shè)定界面讓用戶可以直接輸入各模態(tài)的廣義阻尼、振型放大系數(shù)以及對(duì)稱性。

（5）設(shè)定船體特征參數(shù)

該模塊提供界面讓用戶可以方便地輸入船體特征參數(shù)。

（6）設(shè)定海況參數(shù)

該模塊提供界面讓用戶可以方便地輸入海況參數(shù)。

（7）定義坐標(biāo)系

在三維建模過(guò)程中，干結(jié)構(gòu)模型參考坐標(biāo)系與水彈性計(jì)算的平衡坐標(biāo)系未必一致，因此需要在干結(jié)構(gòu)模型中定義平衡坐標(biāo)系，用以計(jì)算兩個(gè)坐標(biāo)系之間的原點(diǎn)偏移和坐標(biāo)軸夾角余弦。

同時(shí)濕網(wǎng)格模型采用的坐標(biāo)系也未必就是水彈性計(jì)算所需要的濕面元坐標(biāo)系，因此需要在已有的模型中定義濕面元坐標(biāo)系，在生成flxbd.in文件時(shí)，需要將節(jié)點(diǎn)坐標(biāo)轉(zhuǎn)換為濕面元坐標(biāo)系下的坐標(biāo)。

（8）干濕網(wǎng)格單元對(duì)應(yīng)

支持用戶自動(dòng)或者手動(dòng)設(shè)計(jì)干單元組和濕單元組，并設(shè)定干濕單元組的對(duì)應(yīng)關(guān)系。

（9） 求解計(jì)算

該模塊用以設(shè)定求解控制參數(shù)，在用戶確認(rèn)之后，依照之前完成的內(nèi)容，生成輸入數(shù)據(jù)文件，并依次調(diào)用三個(gè)模塊程序，完成水彈性計(jì)算程序的計(jì)算。

（10） 查看結(jié)果

顯示計(jì)算結(jié)果的后處理階段，目前可先將完成主坐標(biāo)響應(yīng)的曲線顯示，并能結(jié)合模型節(jié)點(diǎn)振型，以云圖/動(dòng)畫(huà)的形式顯示模型的運(yùn)動(dòng)響應(yīng)。

其中主坐標(biāo)響應(yīng)可以從結(jié)果文件中直接讀取，結(jié)構(gòu)運(yùn)動(dòng)響應(yīng)可用主振型位移的迭加表示：

4.3 水彈性集成軟件的實(shí)現(xiàn)

按照上面的設(shè)計(jì)思想，我們?cè)谥袊?guó)船舶科學(xué)研究中心自研的Smart集成定制平臺(tái)基礎(chǔ)上，進(jìn)行集成開(kāi)發(fā)，其實(shí)現(xiàn)的軟件界面如圖3所示。

5 應(yīng)用實(shí)例分析

下面以某18萬(wàn)噸級(jí)的散貨船為例，使用水彈性集成軟件對(duì)其進(jìn)行分析。其主尺度如下表所示。

?

該船的干結(jié)構(gòu)有限元模型已經(jīng)采用通用大型有限元軟件Msc.Patran建立（存儲(chǔ)于180K.bdf文件中），其用于水動(dòng)力分析的濕表面網(wǎng)格也已經(jīng)建立（存儲(chǔ)于180K_wet.bdf）。

下文介紹其使用過(guò)程如下：

（1）導(dǎo)入模型：選擇“導(dǎo)入干結(jié)構(gòu)有限元模型文件”選項(xiàng)，干結(jié)構(gòu)有限元模型文件選擇文件“180K.bdf”，選擇濕表面網(wǎng)格文件“180K_wet.bdf”。成功導(dǎo)入后，可看出該有限元模型文件包括11 427個(gè)節(jié)點(diǎn)，15 652個(gè)板單元。

（2）模態(tài)分析：根據(jù)之前選擇的BDF，通過(guò)后臺(tái)調(diào)用Nastran進(jìn)行模態(tài)分析，一般需要位移歸一化，這里共輸出100階模態(tài)。

（3）模態(tài)篩選：根據(jù)4.2節(jié)中所述的關(guān)鍵節(jié)點(diǎn)位移法，先自動(dòng)進(jìn)行總體模態(tài)的初步篩選（共選出20階模態(tài)），然后在此基礎(chǔ)上，通過(guò)人工查看各階模態(tài)的變形圖，并根據(jù)對(duì)稱、反對(duì)稱變形，分別挑選出前四階總體模態(tài)，共8階，見(jiàn)圖4和圖5。

（4）設(shè)定模態(tài)參數(shù)：在挑選出的各階模態(tài)基礎(chǔ)上，有關(guān)模態(tài)信息的主要數(shù)據(jù)都已經(jīng)具備（從模態(tài)分析結(jié)果文件中讀取獲得），但還要再手工設(shè)定各階模態(tài)的廣義阻尼和振型放大系數(shù)參數(shù)。其中廣義阻尼按經(jīng)驗(yàn)取值0.05，由于已經(jīng)實(shí)現(xiàn)位移歸一化，振型放大系數(shù)取默認(rèn)值1.0。

完成模態(tài)參數(shù)設(shè)定后的表格如下表所示：

Mode No. 固有頻率(Hz) 廣義質(zhì)量 廣義阻尼 振型放大系數(shù) 左右對(duì)稱性7 8 9 1 0 11 12 13 14 0.648 43 0.817 23 0.865 72 1.320 4 1.693 1.858 4 1.951 8 2.469 3 2.351 702E+07 1.875 475E+07 2.585 143E+07 1.845 968E+07 1.821 342E+07 2.602 281E+07 4.938 913E+06 4.288 284E+06 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0對(duì)稱反對(duì)稱反對(duì)稱對(duì)稱反對(duì)稱反對(duì)稱對(duì)稱對(duì)稱

（5）設(shè)定船體特征參數(shù)：結(jié)合船體示意圖，手動(dòng)輸入該散貨船的Lpp、艏柱距重心的水平距離、排水體積等參數(shù)。

（6）設(shè)定海況參數(shù)：考慮一個(gè)航速，速度為0；一個(gè)入射波波角，180°（迎浪）；波幅為1.0m，入射波頻率考慮從0.05rad/s開(kāi)始，并以0.05遞增到20，共計(jì)400個(gè)頻率。

（7）定義坐標(biāo)系：這里建立的干結(jié)構(gòu)模型坐標(biāo)系與平衡坐標(biāo)系是平行的，只有原點(diǎn)不一致，因此定義平衡坐標(biāo)系時(shí)只要指定干坐標(biāo)系中重心的坐標(biāo)即可。而濕表面網(wǎng)格模型已經(jīng)在濕面元坐標(biāo)系下建立，因此不需要再額外定義。

（8）干濕網(wǎng)格對(duì)應(yīng)：選擇干結(jié)構(gòu)模型外表面的網(wǎng)格單元為一組（干單元組1），選擇所有的濕表面網(wǎng)格單元為一組（濕單元組1），并確定干單元組1對(duì)應(yīng)濕單元組1。

（9）求解計(jì)算：求解計(jì)算時(shí)需要先指定求解控制參數(shù)，其中輸出數(shù)據(jù)類型，用戶可以指定自己關(guān)心的廣義附加質(zhì)量、廣義附加阻尼的矩陣對(duì)角元素，其他均采用軟件推薦的默認(rèn)值即可。確認(rèn)后，點(diǎn)擊求解計(jì)算按鈕，即可進(jìn)行自動(dòng)的分析求解，軟件將依次調(diào)用原THAFTS程序包中的FLXBD、HYCOF和HYELAS模塊程序。軟件運(yùn)行結(jié)果文件與原有的手工過(guò)程完成的結(jié)果一致。

（10）結(jié)果查看：可進(jìn)行主坐標(biāo)響應(yīng)的曲線顯示、船體運(yùn)動(dòng)響應(yīng)云圖和動(dòng)畫(huà)顯示等，如圖6a和圖6b分別給出的是曲線和云圖顯示的界面。

6 結(jié)論與進(jìn)一步展望

本文在中國(guó)船舶科學(xué)研究中心現(xiàn)有的三維線性頻域水彈性分析程序“THAFTS”基礎(chǔ)上，根據(jù)輸入輸出文件的內(nèi)容和工程使用經(jīng)驗(yàn)，增加其前后處理功能，形成新的工程實(shí)用化的三維水彈性集成軟件，并利用該集成軟件，對(duì)某180k噸散貨船進(jìn)行分析計(jì)算。其使用過(guò)程符合工程軟件的要求，大大提高了三維線性水彈性分析軟件的易用性。

除了本文介紹的三維線性水彈性程序，杜雙興在此基礎(chǔ)上建立的更為嚴(yán)格的帶航速的三維線性水彈性力學(xué)頻域數(shù)值分析方法[9]，田超建立的大浪中航行船舶的三維非線性水彈性分析方法[1]，以及“十一五”期間中國(guó)船舶科學(xué)研究中心發(fā)展的聲介質(zhì)中船舶結(jié)構(gòu)三維水彈性計(jì)算方法[10]，都將統(tǒng)一集成到本文介紹的集成軟件中。

此外，在前處理方面，將逐步實(shí)現(xiàn)常規(guī)船舶（潛艇）的三維參數(shù)化有限元建模，以進(jìn)一步提高效率；在后處理方面，將結(jié)合應(yīng)用領(lǐng)域的不同，如船體運(yùn)動(dòng)、波浪載荷、結(jié)構(gòu)變形、應(yīng)力應(yīng)變以及疲勞性能分析等，集成這些領(lǐng)域的專業(yè)處理方法，形成更加直接的預(yù)報(bào)結(jié)果，并可結(jié)合虛擬現(xiàn)實(shí)技術(shù)，提供情景化的船舶在波浪中運(yùn)動(dòng)的性能分析。

[1]田 超.航行船舶的非線性水彈性理論與應(yīng)用研究[D].上海:上海交通大學(xué),2007.

[2]田 超,吳有生.船舶水彈性力學(xué)理論的研究進(jìn)展[J].中國(guó)造船,2008,49（4）:1-11.

[3]杜雙興.完善的三維航行船體線性水彈性力學(xué)頻域分析方法[D].無(wú)錫:中國(guó)船舶科學(xué)研究中心,1996.

[4]Tian Chao,Wu Yousheng.The second-order hydroelastic analysis of a SWATH ship moving in large-amplitude waves[J].Journal of Hydrodynamics,Ser.B,2006,18(6):631-639.

[5]Tian Chao,Wu Yousheng,Chen Yingqiu.Numerical predictions on the hydroelastic responses of a large bulker in waves[C]//Proceedings of the 5th International Conference on Hydroelasticity in Marine Technology.University of Southampton,UK,8th-10th Sept.,2009:205-213.

[6]Wu Y S.Hydroelasticity of floating bodies[D].Brunel University,UK,1984.

[7]王志軍.箱式超大型浮體結(jié)構(gòu)的水彈性響應(yīng)研究[D].上海:上海交通大學(xué),2001.

[8]戴愚志,余建星.一種船體及周圍自由面的網(wǎng)格自動(dòng)生成方法[J].船舶工程,2006,28(5):118-123.

[9]Du S X,Wu Y S,Price W G.Forward speed effect on the structure responses of a ship traveling in waves[C]//Proceedings of 2nd International Conference on Hydroelasticity in Marine Technology.Fukuoka,Janpan,1998:401-410.

[10]金建海,謝基榕,鄒明松等.潛艇定量聲學(xué)設(shè)計(jì)計(jì)算體系設(shè)計(jì)說(shuō)明書(shū)[R].無(wú)錫:中國(guó)船舶科學(xué)研究中心報(bào)告,2010.

Design and implementation for 3D hydroelasticity integrated software

JIN Jian-hai1,2,YE Yong-lin2,TIAN Chao3,ZOU Ming-song2,LENG Wen-hao2
(1 School of IoT Engineering,JiangNan University,Wuxi 214122,China;2 China Ship Scientific Research Center,Wuxi 214082,China;3 China Classification Society,Beijing 100007,China)

In order to satisify the requiments of the commercial application of 3D hydroelasticity program,the pre-post processing modules for 3D hydroelasticity were designed and realized,and the integrated software was developed based on analyzing the experience of using the existing program.Then the example of a 180kton bulk cargo was given to illustrate how to use the Integrated Software and the basic functions of the software.Finally,the further function of the software was recommended.

3D hydroelasticity;pre-post processing;mode selection;principal coordinate response

U662.9

A

1007-7294（2011）05-0521-09

2011-02-22

“十一五”國(guó)家科技支撐計(jì)劃制造業(yè)信息化工程重大項(xiàng)目（2006BAF01A43）

金建海（1978-），男，江南大學(xué)博士研究生，中國(guó)船舶科學(xué)研究中心高級(jí)工程師；

葉永林（1974-），男，中國(guó)船舶科學(xué)研究中心高級(jí)工程師。