船體水動力網格動態劃分技術及應用

朱雄軍，李紫麟

(1.武漢職業技術學院 計算機學院，武漢 430074； 2.中國船舶科學研究中心，江蘇 無錫 214082)

船體水動力網格動態劃分技術及應用

朱雄軍1，李紫麟2

(1.武漢職業技術學院 計算機學院，武漢 430074； 2.中國船舶科學研究中心，江蘇 無錫 214082)

摘要：考慮到船舶與波浪耦合作用下，船體濕表面時刻變化，引起非線性動力效應，根據網格單元與波浪的不同接觸類型，分別采用拆分、合并和局部調整得到新的船體濕表面網格，選取波浪環境中線性規則波2種情況，對簡體船舶模型在波浪環境下橫搖和縱搖進行算法驗證，并對線性波中集裝箱船濕表面網格動態劃分過程進行數值模擬，驗證算法的正確性與實用性，為波浪中船體非線性動力特性的時域計算提供技術支持。

關鍵詞：船舶；時域分析；動態劃分；網格重構；數值模擬

近年來，網格技術廣泛應用于各種工程和科學的數值計算。網格生成[1]，網格劃分[2]，網格優化[3]，動網格技術[4]等諸方面關鍵技術得到發展和完善。模型的網格生成與劃分是有限元前處理的重要部分，在數值計算中通常采用初始劃分好的網格。在某些研究領域中，網格單元和網格形狀需要隨時間進行調整，進行網格重新劃分(remeshing)，完成從舊網格到新網格的信息轉換[5]。在網格重新劃分的研究中，張玉新等[6]應用有限元分析軟件(MSC1SuperForm)中的網格再劃分功能，對鋁型材的擠壓變形過程的三維有限元模擬；丁永祥和夏巨[7]基于等參轉換概念及應用小區理論對剛塑性有限元網格重新劃分算法進行了研究；吳淑芳等[8]首次考慮了四邊形網格被直線切割后各種幾何特征，并給出相應的混合調整算法。利用勢流理論進行船舶與海洋結構物水動力分析時，常采用格林函數法，計算時需對船體平均濕表面進行網格劃分[9-10]。但考慮物體濕表面積及吃水變化引起的非線性動力效應時，需要在每一時間步重新計算濕表面，對變化的水線面附近網格進行動態的重新劃分。為此，針對船舶水動力計算中網格動態劃分和重構的關鍵技術進行研究。波浪環境分別選取線性規則波2種情況。船體結構由四邊形單元與三角形單元構成?？紤]波浪與船體單元間的接觸判斷確定需要調整的網格。根據網格單元與波浪的不同接觸類型，分別采用拆分、合并和局部調整獲取新的船體濕表面網格。選取簡化的船體模型在波浪環境下橫搖和縱搖進行算法驗證。最后應用網格動態劃分技術，對時域線性波浪下集裝箱船濕表面網格動態劃分過程進行數值模擬。

1網格動態劃分的理論依據

在運算之前，需要讀取船體單元和水線面的參數數據。如圖1所示，波浪環境選取直線表示靜水狀態，可通過首尾吃水計算出水線位置；曲線表示線性規則波浪，通過波高，頻率，相位描述波浪的位置信息。

圖1　二維船體在波浪中的簡化幾何模型

對船體外表面采用四邊形網格和少部分三角形網格進行劃分，因此船體表面信息由節點和單元的集合構成。在算法中采用有限元網格編號方法，即先尋找一個邊角節點，使該節點的編號為1，新編號節點沿著已經編號的節點的前沿向外擴展，直到所有的節點編完為止。有限元網格生成以后，單元的信息即可獲得。這里用x(i)和y(i)分別表示節點位置屬性，elementinf(i,1∶4)表示構成i號四邊形單元網格的節點編號，按照逆時針順序記錄。此外，構成船體表面的少量三角形單元，可由四邊形退化而成，將該單元的第4個節點與第3個節點重合即可。

船體四邊形網格單元與水線面的接觸判斷是決定網格是否需要重構的直接因素。四邊形與水線面的位置關系一般包括3種情況，如圖2示：l1與四邊形ABCD不相交；l2與四邊形ABCD的兩條邊相交；l3過四邊形ABCD頂點并與其相交?？梢园l現，當四邊形中的任意一組對角頂點在直線的兩側的時候，直線穿過四邊形，將其切割成兩個部分。否則，直線未將四邊形剖分。因此，水線面與四邊形單元的接觸判斷可以轉化成單元的4個節點與水線面的位置關系進行判斷。

圖2　四邊形與直線位置關系

通過接觸判斷，依次記錄與水線面有接觸的四邊形單元，將其單元編號記錄在mark(i)(i=1～nmark)，nmark為與水線面接觸的單元總數。接觸單元的交點位置可通過對水線面與網格的接觸邊的方程組求解得出。圖3，通過mark(i)即可獲得對應目標單元的存儲信息，包括接觸的單元信息、單元節點信息、單元鄰居列表、單元接觸類型與單元接觸點的位置。

圖3　與水線面接觸的四邊形單元的存儲方式

動態網格劃分的實質，是對與水線面接觸的船體表面單元進行處理，通常采用拆分、合并或區域重組的方法加以實現?？刹捎酶F舉法，將線段EF與四邊形網格ABCD可能的全部剖分形式列舉出來, 建立接觸類型數據庫。然后針對不同的重構方式，將計算得到的某節點、邊或四邊形單元附近的區域的單元信息，與給定的“類型”的定義進行匹配，以確定該局部區域網格是否需要重構及相應的重構類型。最后調用所對應的操作，對目標網格進行重構。

2動態網格算法實現

可視化處理有多種方式，對二維數據場的可視化處理方法包括X-Y(曲線)圖法、極坐標圖法、二維網格法、二維矢量圖法、二維等值線法及柱狀圖法等，對三維數據場的可視化處理方法包括三維網格法、三維矢量圖法、三維等值線法及離散點圖法等。為了能夠清晰的表達三維仿真數據的結構，可視化處理還包括透明度處理、燈光處理、陰影處理和邊界處理等方法。此外，數據消隱、坐標軸設置、可視化結果局部控制、動畫仿真等一系列高級功能能讓仿真數據可視化的結果更形象直觀、也更易理解。

Tecplot提供的可視化處理模塊作為系統的可視化處理引擎。Tecplot其提供的宏編程接口允許用戶對其進行二次開發。

分流程圖見圖4，可分為4個階段：①程序初始化，在運算開始前獲得水線面的函數表達式及網格信息；②計算出與水線面接觸的四邊形單元的信息，包括接觸點，接觸類型兩部分；③對需要調整和優化的網格進行重構，并對網格單元信息優化，去掉重復節點和多余單元；④計算時間增量，更新水線面位置信息，進入下一個時間步。

圖4　動態網格劃分流程

3算例驗證

3.1水線面切割下的船體網格單元重構

在算例中，選用如圖5所示的長方體模型進行計算方法驗證。使用ANSYS建模并對其進行網格剖分，生成252個四邊形單元及250個節點。將水線面簡化為水平面，垂向高度為模型高的一半。圖中為初始時刻網格的劃分情況。以水線面為界將模型分為2個部分，包括水線面以上的單元用3層網格表示，以及中間1層和最底下2層組成的為水下單元，其中中間1層為水線面附近經過動態劃分后重新構成的單元。

圖5　簡化的模型單元網格初始狀態

考慮靜水情況，保持水線面不變，令模型分別進行強迫橫搖和縱搖運動，重新劃分后的網格效果如圖6a)和b)所示。進一步考慮到規則波情況，假定物體固定不動，頂浪作用下網格的變化情況如圖6c)和d)所示。此處余弦波可采用差分法，將波函數曲線離散成若干與單元網格相交的線段。這樣波浪與模型仍然可采用直線與物體的接觸判斷準則。

圖6　模型單元網格重構結果

3.2非線性波作用下的集裝箱船網格動態劃分

試驗選取集裝箱船為研究對象，模擬集裝箱船垂蕩與縱搖運動下，船體網格的動態劃分過程。選取規則波，頂浪狀態。船體的運動響應可以通過水動力分析計算，分別測定垂蕩和縱搖運動的固有周期、幅值和相位。具體參數見表1。

表1　相關參數

在時域模擬過程中，選取0、5、 10 s共3個不同時刻，觀察在對應的狀態下船體進行網格動態劃分過程。圖7展示了四邊形網格被規則波余弦曲面切割后，經過重構生成的新四邊形網格劃分圖。通過使用拆分、合并和局部調整的方法，根據水線面的位置，對水線面附近的網格單元進行重構。中間曲線向上部分為水面以上的單元，曲線部分是水下固定不變的單元，曲線邊緣鋸齒狀部分為水線面附近重新構成的單元。在進行非線性水動力分析時，即可選用重新生成的鋸齒狀單元和曲線下的固定不變單元組成的船體濕表面網格進行計算。

圖7　動態劃分后重構模型單元網格效果

4結論

給出船體表面二維四邊形網格被直線切割的各種幾何拓撲特征，并提供相應的處理方法，對船體表面網格被靜水線面以及線性規則波面切割的區域進行網格重構，建立接觸類型數據庫。建立船體水動力網格的動態劃分方法，通過靜水情況模型的強迫橫搖和縱搖運動，以及頂浪狀態下模型的強迫橫搖和縱搖運動，清晰地展現水線面切割下的船體網格單元重構的應用和效果。仿真模擬集裝箱船垂蕩與縱搖運動下，船體網格的動態劃分過程，驗證了算法的正確性與實用性。本文工作可為波浪中船體非線性動力特性的時域計算提供技術支持。

參考文獻

[1] 關振群,宋超,顧元憲,等.有限元網格生成方法研究的新進展[J].計算機輔助設計與圖形學學報,2003,15(1):1-14.

[2] 古成中,吳新躍.有限元網格劃分及發展趨勢[J].計算機科學與探索,2008,2(3):248-259.

[3] 李迎華,吳寶山,張華.CFD動態網格技術在水下航行體非定常操縱運動預報中的應用研究[J].2010,14(10):1100-1108.

[4] JO E B. Quadrilateral mesh generation in polygonal regions[J]. Computer aided design,1995,27:194-199.

[5] 楊曉松,顧元憲,李云鵬,等.有限元網格體繪制中的剖切算法[J].中國圖象圖形學報,2002:7(A)(1):55-62.

[6] 張玉新,鐘毅,王華昆,等.塑性成形中的六面體網格劃分與重劃分技術[J].鍛壓技術,2004(6):33-35.

[7] 丁永祥,夏巨湛.剛塑性有限元模擬過程中的網格重新劃分算法[J].華中理工大學學報,1995,23(I):83-87.

[8] 吳淑芳,李占國,管力銳,等.四邊形網格被直線切割的混合調整法[J].長春光學精密機械學院學報,1998,21(2):20-24.

[9] FALTINSEN O M. Sea loads on ships and offshore structures[M]. Cambridge: The Cambridge University Press,1990.

[10] 王俊榮,謝彬.半潛式平臺水動力性能及運動響應研究綜述[J].中國造船,2009,50(增刊):255-260.

[11] 李鵬飛,徐敏義,王飛飛,精通CFD工程仿真與案例實戰[M],北京：人民郵電出版社,2011.

[12] 藍霄峰,陳娓,向波,等.Tecplot在大鏡山水庫洪水演進動態可視化的應用[J].人民珠江,2010(1):60-61

Ship Hull Dynamic Remeshing Technique and Its Application

ZHU Xiong-jun1, LI Zi-lin2

(1 Dept. of Computer Science, Wuhan Polytechnic, Wuhan 430074, China;2 China Ship Scientific Research Center, Wuxi Jiangsu 214082, China)

Abstract:When a ship motions in waves, the wet surface is changing continually, causing the nonlinear dynamic effect,. According to the different contact types of grid cell and the waves, respectively adopt split, merge and partial adjustment to get the new mesh of the wet surface. Two kinds of linear regular wave are selected for the wave environment, the dynamic remeshing algorithm is tested and verified by a simplified ship model with motions of rolling and pitching. And the dynamic meshing process for the wet surface of a container ship in the linear wave is simulated to verify the practicability of the algorithm. It can provide the technical support for the time domain calculation of the nonlinear dynamic characteristics of the ship hull.

Key words:ship; time domain analysis; dynamic remeshing; grid reconstruction; numerical simulation

DOI:10.3963/j.issn.1671-7953.2016.03.010

收稿日期：2015-10-26

第一作者簡介：朱雄軍(1966—)，男，碩士，副教授 E-mail:zhuxiongjun@126.com

中圖分類號：U662.9

文獻標志碼：A

文章編號：1671-7953(2016)03-0043-04

修回日期：2015-12-11

研究方向:計算機軟件與理論