楊大偉 李程 孫淼

摘要：以優(yōu)化艦船穩(wěn)性計算為出發(fā)點，提出了一種基于STL模型的艦船排水體積與浮心位置計算方法。與以往二維計算方法相比大大減少了計算的誤差，提高了穩(wěn)性計算的準(zhǔn)確性與可靠性。并使計算結(jié)果與NAPA軟件的靜水力計算結(jié)果對比，驗證了算法的準(zhǔn)確性與可行性。

關(guān)鍵詞：艦船穩(wěn)性;STL模型;排水體積;浮心

引言

艦船排水體積與浮心位置的計算具有十分重要的意義，是艦船穩(wěn)性、艦船浮態(tài)等計算的基礎(chǔ)。在傳統(tǒng)的穩(wěn)性計算過程中，設(shè)計人員利用二維的線框模型通過插值及數(shù)值積分的方法進行計算，數(shù)值積分點的有限性及插值法都會造成計算精度不高，整個計算過程也十分復(fù)雜。

隨著現(xiàn)代數(shù)字造船技術(shù)的不斷發(fā)展，目前在船舶設(shè)計軟件當(dāng)中，大多可以直接輸出船舶建模的STL（Stereo Lithography）三角網(wǎng)格數(shù)據(jù)。這種數(shù)據(jù)能夠精準(zhǔn)的描述船體模型信息，相比于數(shù)字型值表的表達方式，其在穩(wěn)性計算等方面將有無可比擬的優(yōu)勢。

1 STL模型的建立

STL（STereo Lithography）文件是一種文件格式，用來表現(xiàn)三維物體的幾何信息，在3D打印、快速成形等領(lǐng)域中的應(yīng)用非常廣泛。STL格式表示的模型文件是用三角形面片的信息構(gòu)成物體的表面，可以認為是將模型的曲面進行了三角化的處理。

STL模型三維數(shù)據(jù)格式簡單、跨平臺性能優(yōu)良、在計算幾何特性時簡單和與平面、曲面求交時計算速度快。船舶的外殼及艙室STL模型可由以下兩個途徑獲得：（1）對于存在設(shè)計數(shù)據(jù)的船舶，可由船舶設(shè)計軟件導(dǎo)出，比如可用NAPA導(dǎo)出IGES格式，然后再轉(zhuǎn)成STL格式;（2）對于得不到船舶設(shè)計數(shù)據(jù)的船舶，就需要根據(jù)船舶的型值表、靜水力數(shù)值表、總布置圖等數(shù)據(jù)利用建模工具對其進行三維重建。

2 排水體積與浮心位置計算

2.1 排水體積計算

在傳統(tǒng)的排水體積計算方法中，通常采用二維積分法進行計算，此方法計算過程較為復(fù)雜。而利用STL三角網(wǎng)格模型進行排水體積計算時，采用的是投影法，通過將模型曲面的三角面片向船體水線面進行投影，三角形的頂點在水線面得到相應(yīng)的投影點，再依次連接曲面中三角面片頂點和水線面上的投影點構(gòu)造五面體，最后計算出五面體的體積進行疊加，即可得到整個船體模型的排水體積。

假設(shè)三角面片上的某一個頂點為點A，與其相對應(yīng)的投影點為點B，其中水線面的法向量為n，則有如下關(guān)系式：

式中，dAB為A點到水線面的距離，dn為平面法向量的模長。

因此只需找到A、B兩點之間的關(guān)系，便可得出B點的坐標(biāo)，但需要注意在計算過程中，首先需要判斷法向量n是否為零向量。

由于構(gòu)造出的五面體并不一定是標(biāo)準(zhǔn)正五面體，計算過程會較為不便，因此在計算時，為提高計算精度，可以將構(gòu)造的五面體進行劃分，分為三個三棱錐分別進行體積計算。計算三棱錐體積的公式為：

式中，h為三棱錐高，S為底面面積。

如圖2-1所示，模型中的一個三角面片頂點123與其相對應(yīng)的投影點1’2’3’依次連接構(gòu)成五面體123-1’2’3’，將五面體進行劃分得到三棱錐2-1'2'3'、三棱錐1-323’和三棱錐3’-211’。

利用五面體體積疊加的方法在計算整個模型體積時會出現(xiàn)某一部分體積重復(fù)計算的情況，因此，需要對五面體的體積進行區(qū)分，根據(jù)模型中原三角面片的法向量與其對應(yīng)的投影面法向量夾角，計算結(jié)果分正、負、零三種情況：假設(shè)投影面法向量為正方向，則若兩法向量夾角為鈍角，該五面體是需要計算保留的部分，此時體積取正值;若兩法向量夾角為銳角，該五面體在計算過程中是需要去掉的部分，此時五面體體積取負值;若兩法向量夾角為直角，此時五面體體積取零值，示意圖如圖2-2所示。整個模型的體積為所有體積相加之和。

2.2 浮心位置計算

在傳統(tǒng)的浮心位置計算中，通常采用的是計算體積矩的方法。在基于STL模型的方法計算浮心位置時，需要計算出所有五面體的形心位置，在計算時，同樣需要按上節(jié)的方法把五面體分成三個三棱錐，計算出每個三棱錐的形心位置，然后求出每個三棱錐的體積矩，最后將三個體積矩相加并除以該五面體體積，即得到五面體的形心位置。同理可得出整個模型形心位置。

下面以實例進行簡要介紹。在該例中，使用二維坐標(biāo)系進行闡述，三維坐標(biāo)系同理。如圖2-3所示，計算矩形ABCD的形心位置，已知AB長度為2，AD長度為1，投影線距離X軸為2。

依據(jù)常識我們可判斷出形心位置坐標(biāo)為（1，0.5）。若使用STL網(wǎng)格模型方法進行計算時，首先將矩形的AB邊向投影線投影，，聯(lián)系地得到相應(yīng)的投影點EF，則矩形ABEF的面積為4，形心坐標(biāo)為（1，1）。然后將CD邊向投影線投影，得到相應(yīng)的投影點MN，CDMN為需要減掉的部分且面積為2，取值-2，形心坐標(biāo)為（1，1.5）。整個矩形的形心位置計算公式如下：

可以看出利用STL模型投影計算形心位置時，為避免出現(xiàn)較大誤差，投影體積需考慮正負值。

3 排水體積和浮心位置算法驗證

本文基于Geomagic Studio軟件平臺，以一艘散貨船為例，運用上述算法完成了三角網(wǎng)格模型計算程序的編制。表3.1、表3.5為基于STL模型的計算方法與NAPA軟件計算數(shù)據(jù)結(jié)果。

計算結(jié)果表明，基于STL網(wǎng)格模型的計算結(jié)果與NAPA計算結(jié)果相比，型排水體積與排水量最大相對誤差為0.31%，最小為0.08%;浮心縱向坐標(biāo)與浮心橫向坐標(biāo)計算結(jié)果最大相對誤差為0.23%，最小為0.12%。差別非常小，驗證了算法的準(zhǔn)確性。

4結(jié)束語

本文提出了一種基于STL網(wǎng)格模型的艦船排水體積與型心位置計算方法，其中艦船STL網(wǎng)格模型來源于艦船設(shè)計數(shù)據(jù)，可保證計算精度。文中提出的算法在保證計算精度的基礎(chǔ)上力求保證程序的實時性，計算結(jié)果證明了新算法的可行性與實用性，具有一定的工程實用價值。

參考文獻

[1]張明霞. 基于NURBS曲面的船舶破艙穩(wěn)性計算方法研究[D]. 大連理工大學(xué)，2002.

[2]孫殿柱，孫永偉，田中朝. 三角網(wǎng)格曲面模型快速求交算法[J]. 北京工業(yè)大學(xué)學(xué)報，2012，08： 1121-1125.

[3]盛振邦，劉應(yīng)中.船舶原理（上冊）[M].上海：上海交通大學(xué)出版社，2003：20-24.

[4]姜玲.船舶穩(wěn)性三維計算方法研究及軟件開發(fā)[D].大連：大連海事大學(xué)，2015.

[5]張小青.基于三角網(wǎng)格模型的文物幾何信息計算[D]北京：北京建筑工程學(xué)院，2012：17-37.

[6]曹晶，陳明，孫永剛?；贑ATIA的三維船舶靜水力計算研究[J].中國艦船研究，2011，6（2）：25-28.

[7]王智洲，孫霄峰，尹勇，劉春雷. 散貨船配載儀中完整穩(wěn)性計算方法研究[J]. 大連理工大學(xué)學(xué)報，2017，05： 476-481.

[8]孫殿柱，孫永偉，田中朝. 三角網(wǎng)格曲面模型快速求交算法[J]. 北京工業(yè)大學(xué)學(xué)報，2012，08： 1121-1125.

[9]張晞，孟立新.基于三角網(wǎng)格模型的三維CAD模型的建立[J]. 制造技術(shù)與機床，2008，07.

[10]IMO SLF53/INF.10. Information collected by the Correspondence Group on Intact Stability[R]. 2011.

[11]賈敬蓓，宗智，師華強，馬駿. 不對稱型雙體船初穩(wěn)性研究[J]. 中國艦船研究，2008，01.

作者簡介：楊大偉，出生年月：1979.5.3，性別：男，單位：海軍航空大學(xué)，技術(shù)職稱：副教授，主要從事工作或研究方向：航海技術(shù)，船舶結(jié)構(gòu)力學(xué)