盛慶武

(上海佳豪船舶工程設(shè)計股份有限公司 上海 201612)

工程船船體濕表面網(wǎng)格NAPA軟件自動生成技術(shù)*

盛慶武

(上海佳豪船舶工程設(shè)計股份有限公司 上海 201612)

NAPA;鐘型隸屬度函數(shù);縱傾;自動劃分

利用NAPA軟件中二次開發(fā)工具宏語言,通過引入模糊數(shù)學(xué)中的鐘型隸屬度函數(shù)和二次函數(shù),開發(fā)出對平浮和帶縱傾狀態(tài)的船體濕表面網(wǎng)格自動劃分模塊。該模塊不但可以快速和準(zhǔn)確地劃分船體網(wǎng)格,還可以生成直接被HydroStar軟件調(diào)用的輸入文件,并且通過控制輸出文件的格式可以生成符合其他水動力軟件調(diào)用的船體濕表面網(wǎng)格格式。

0 引 言

在船舶耐波性計算中,三維勢流理論目前已成為工程船舶解決耐波性問題的重要理論。通常,以三維勢流理論開發(fā)的耐波性計算軟件需要劃分多個不同吃水和縱傾下的船體濕網(wǎng)格,因此,船體濕網(wǎng)格劃分需要快速而準(zhǔn)確地完成,特別是在工程船舶設(shè)計中,其設(shè)計周期一般比較短,如果建立的新方法雖然具有很強的理論性,但實用性差,依然不能在工程船設(shè)計中占據(jù)一席之地。所以,本文擬依據(jù)工程船舶設(shè)計的特點,建立一套船體濕網(wǎng)格快速和準(zhǔn)確劃分的實用方法。

NAPA軟件[1]具有比較強的二次開發(fā)功能——宏語言,應(yīng)用NAPA中的宏,可以有針對地編寫適合的宏代碼,達到預(yù)期的功能。

在船舶設(shè)計中,開發(fā)船舶的型線是最早開展的工作之一,在初步確定船體主尺度后,通過改造適合的一條或幾條母型船,應(yīng)用NAPA軟件可以方便、可靠和快速地建立新船的型線模型,進而才能開展后續(xù)的工作。

在NAPA軟件中型線設(shè)計完成后,即可以獲得相應(yīng)的型線圖或肋骨型線,如果需要劃分網(wǎng)格,一般情況下,需要在NAPA輸出的型線圖或肋骨型線圖(大多為DXF文件格式)上處理。這樣處理可以通過其他的軟件或者AUTOCAD軟件自帶的二次開發(fā)語言完成。如果能夠直接在NAPA軟件中完成網(wǎng)格自動劃分,這自然是非常好的方法,可以省卻各種軟件格式之間的轉(zhuǎn)化。

鑒于NAPA軟件具有比較強的二次開發(fā)功能,本文借助NAPA軟件中二次開發(fā)語言——宏語言,在已建成的NAPA型線模型基礎(chǔ)上,開發(fā)出可以直接生成耐波性計算軟件HydoStar要求的網(wǎng)格形式,輸出的網(wǎng)格形式文件可以直接被HydroStar軟件[2]調(diào)用。這種網(wǎng)格劃分方法可以非常簡便、準(zhǔn)確和快速地完成船體濕表面網(wǎng)格的自動生成工作,并且網(wǎng)格的精細和疏密可以任意控制:既可以達到船體耐波性計算需要的網(wǎng)格數(shù)量:幾百塊至幾千塊,還可以達到計算流體力學(xué)諸如FLUENT等的幾萬塊甚至十幾萬塊的網(wǎng)格數(shù)量等。因此可認(rèn)為這種方法不失為比較好的船體濕網(wǎng)格劃分方法。

船舶在實際航運的工況很少有平浮狀態(tài),這在船舶的典型裝載手冊中可以看出,該文件需要由船級社審查。在建立好各種裝載工況后,通過船舶不同裝載工況下的耐波性能計算,可以更加準(zhǔn)確地預(yù)報該船的耐波性的優(yōu)劣,更加準(zhǔn)確地預(yù)報該船錨泊定位能力,這樣可以從船舶耐波性能的角度比較便捷地指導(dǎo)船舶設(shè)計。應(yīng)用本文的方法,可以非常方便地對已計算的各種裝載工況進行網(wǎng)格劃分,這樣各種裝載工況相當(dāng)于船體吃水和縱傾的函數(shù)。如N個工況,就有N個吃水和縱傾,NAPA中很容易讀出這些量的數(shù)據(jù),通過一個循環(huán)語句,在NAPA宏語言程序中運行一次,就可獲得所有裝載工況下的網(wǎng)格劃分,真正做到了快速和準(zhǔn)確的效果。

本文介紹了2 400 t起重船帶平浮狀態(tài)和縱傾(艉傾和艏傾)狀態(tài)下的網(wǎng)格劃分的原理和方法,使用NAPA軟件中的宏語言開發(fā)網(wǎng)格自動劃分模塊后,可以方便、快速和準(zhǔn)確地研究船體耐波性能和錨泊定位系統(tǒng)設(shè)計。

1 NAPA在處理船體線型的方法

NAPA軟件在處理船體型線上是先建立若干條基本控制線:如平邊線、平底線、首尾輪廓線、平行中體首尾肋位線、尾封板線等,然后通過其他控制線:如肋骨線、水線、縱剖線或者斜剖線等輔助控制線組成船體的骨架,然后通過使用曲線插值成曲面的方法建立整個船體曲面,也就是由點即線,由線到面的方法。

在NAPA軟件中,使用三次樣條對各控制點進行擬合:三次樣條插值曲線的方法有兩種:M1是曲線的曲率不連續(xù),即二次微分;M2是更連續(xù)的曲率。圖1顯示了兩種樣條插值的方法。可以看出,M1是按各段進行分段擬合曲線,而M2則所有的點均用以擬合曲線。

圖1 NAPA軟件中樣條擬合的兩種方法

由于NAPA軟件處理船體曲面的方法是通過控制線形成光順的曲面?？梢奛APA軟件處理船體型線具有三維性,這樣獲得的三維船體曲面可以很方便地剖出各種剖面線:橫剖線、縱剖線、水線面和斜剖線。由于NAPA軟件在光順后的船體曲面上可以剖出任何位置的剖面線,因此這樣的處理方法可以獲得較精確的船體網(wǎng)格線。圖2為NAPA的船體編輯器(HULL EDITOR)里的一個模型。圖3為該船體幾何窗口(GEOMETRY WINDOW)下顯示的船體三維圖形。

圖2 NAPA船體編輯器里的一個模

圖3 船體幾何窗口顯示的船體三維圖形

由圖3可以看出,對于已經(jīng)在NAPA模型中生成的三維船體模型,可以對其任意切割相交成橫剖線、縱剖線和斜剖線,本文正是利用NAPA軟件的這個特點,依托NAPA軟件中二次開發(fā)宏語言,對船體劃分網(wǎng)格,原則上可以任意劃分,不僅可以劃分平浮吃水,而且還可以考慮船體縱傾的影響,劃分成帶縱傾的船體網(wǎng)格[4]。為了能更充分利用HydroStar軟件,經(jīng)過合理的開發(fā),可以使生成的網(wǎng)格文件直接被HydroStar軟件調(diào)用,不僅可以大大節(jié)省劃分網(wǎng)格的時間,而且還可以保證網(wǎng)格劃分的精度。

2 平浮吃水船體網(wǎng)格劃分

由上節(jié)可知,對于NAPA軟件,如果在縱向給定坐標(biāo)可以獲得該處的橫剖線(如圖4所示),圖中x1-x12分別為x=1m至x=12m,且間隔為1m處的橫剖線。

由于服務(wù)對象身份及年齡的多樣性，比如學(xué)生、社會在職人士等，其閑暇時間是不用的。因此，在課程的時間安排上要做到靈活，能夠最大程度上滿足不同群體的需求。

圖4 船體幾何窗口下船體橫剖線

船體網(wǎng)格劃分依據(jù)NAPA軟件船體編輯器里的船體建立的思想:將船體也分成船首、船舯和船尾三部分分段,分別對每部分分段進行劃分。

在NAPA編寫宏時,設(shè)定船首分段的橫剖面數(shù)為NF、船舯分段的橫剖面數(shù)NM和船尾分段的橫剖面數(shù)NA。在NAPA宏語言中,有一個PLENGTH語句,可以獲得曲線上給定長度的點,進而獲得該點的三個軸的坐標(biāo)。眾所周知,波浪能量隨水深的變化成指數(shù)級衰減,所以近水面的網(wǎng)格要求細些,遠離水面網(wǎng)格可劃分粗些,因此,對于任一已剖出的橫剖線,采取如下的方法來選取其上的點作為劃分網(wǎng)格的點:

設(shè)該橫剖線長度為L,這可以由NAPA軟件中PLENGTH語句直接獲得,N為劃分該橫剖線的網(wǎng)格點,ΔLi為橫剖線上第i段的劃分的長度,它可表示為:

對于橫剖線第i-1段的劃分的長度,i=2表示起點處橫剖線劃分的長度,i=N+1表示終點處橫剖線劃分的長度?？梢钥闯鰴M剖線上劃分的長度以(2·i-1)/N/N數(shù)值遞增,這樣做也可以認(rèn)為沿水深方向逐漸增大,在水線面處的網(wǎng)格尺寸最小。圖5為船體某橫剖線劃分圖。

圖5 船體某橫剖線劃分圖

眾所周知,通常船舶的首尾處線型變化較大,而在船的中間部位線型不變或變化很緩,因此,縱向網(wǎng)格劃分時通常首尾部網(wǎng)格密些,船的中間部位保持相同較疏的網(wǎng)格。如何在NAPA中實現(xiàn)呢?

本文應(yīng)用模糊邏輯控制理論中模糊語言變量的隸屬度函數(shù)[3]—鐘型隸屬度函數(shù)實現(xiàn)船體網(wǎng)格數(shù)沿縱向的變化。

一般的鐘型隸屬度函數(shù)表達式為:

其中:參數(shù)x指定變量的論域范圍,a、b和c分別為指定鐘型函數(shù)的形狀,y為鐘型隸屬度函數(shù)。通過選取不同的a,b,c,可以適應(yīng)不同的船體沿船長方向的網(wǎng)格劃分。

以下以2 400 t起重船船體濕表面網(wǎng)格劃分為例說明應(yīng)用NAPA軟件對船體網(wǎng)格劃分的方法。

設(shè)船長為96m,這樣艏部分段長度為10 m,即艏部距平行中體前段的距離為10 m;艉部分段長度為12m,即艉部距平行中體后段的距離為12m;舯部分段為74m。擬在縱向使用42條橫剖線:艉部分段在0～12m間用6條橫剖線劃分成5段,艏部在86 m～96m間用5條橫剖線劃分成4段,中間12 m～86 m用31條橫剖線劃分成30段。這樣可選擇a=44、b=8、c=48,則鐘型隸屬度函數(shù)曲線如圖6所示。

因為本船艉部縱向范圍是從0～12m,從圖6可知鐘型隸屬度函數(shù)y的值不為1,所以應(yīng)按照此函數(shù)來劃分0～12m范圍內(nèi)的間隔

這里n=5,由此公式則可將0～12m范圍內(nèi)分隔成表1所示。如果想加密縱向網(wǎng)格數(shù),則取n=11,如表2所示。

圖6 鐘型隸屬度函數(shù)曲線

表1 艉部0～12 m范圍內(nèi)的橫剖線位置(分隔為5段)

表2 艉部0～12m范圍內(nèi)的橫剖線位置(分隔為11段)

從表1和表2中可以看到,通過選擇不同的n,可對艉部劃分為不等間距;同樣,可以將12 m～86m區(qū)間及86 m～96 m區(qū)間劃分成不同間隔的縱坐標(biāo)值。

對于艉封板的處理,也可以通過對艉部的不等距處理而獲得相應(yīng)的網(wǎng)格。圖7為該船艉封板的橫剖線劃分圖。圖中粗實線為艉封板的邊界線,按照一定的縮小比例形成幾個橫剖線。最后的橫剖線收縮成一點。

經(jīng)過NAPA軟件宏語言的處理,可以很輕松地劃分出任意吃水和任意網(wǎng)格數(shù)的船體網(wǎng)格,只要給定橫剖線的網(wǎng)格點數(shù)以及總的橫剖線數(shù)及艏艉橫剖線數(shù),并對鐘型隸屬度函數(shù)中的a、b、c選取適合的值,即可劃分網(wǎng)格成功。圖8和圖9為本船在4.357 m吃水下,兩種不同網(wǎng)格數(shù)的網(wǎng)格?？梢钥闯?對于選取不同的網(wǎng)格點數(shù)以及總的橫剖線數(shù)及艏艉橫剖線數(shù),可以劃分出任意精度的船體網(wǎng)格。

圖7 艉封板的橫剖線劃分圖

3 縱傾狀態(tài)下船體網(wǎng)格劃分

圖8 1 008塊船體網(wǎng)格劃分

圖9 1 728塊船體網(wǎng)格劃分

通常情況下,在船舶耐波性性能計算中,一般只考慮平浮吃水下的網(wǎng)格劃分,很少考慮到縱傾的影響。通常舶舶典型裝載工況的浮態(tài)很少有平浮狀態(tài),對于工程船舶,其縱傾可能會很大,因此有必要考慮縱傾狀態(tài)與平浮狀態(tài)下對船體六個自由度的耐波性能影響程度,這方面的工作將在下一節(jié)展開。因此,本節(jié)以2 400 t起重船縱傾狀態(tài)下船體網(wǎng)格劃分為例闡述縱傾狀態(tài)下船體網(wǎng)格劃分方法。

縱傾狀態(tài)下船體網(wǎng)格劃分同平浮吃水船體網(wǎng)格劃相似,差別是縱傾狀態(tài)下船體網(wǎng)格劃分需要對平浮吃水下船體做坐標(biāo)變換,變換成縱傾吃水的狀態(tài),因為船舶縱傾狀態(tài)的艏艉吃水不一樣,所以由需要在船體平浮狀態(tài)下船體旋轉(zhuǎn)一個縱傾角?？v傾坐標(biāo)變換圖如圖10所示:

圖10 縱傾坐標(biāo)變換圖

圖10中X-O-Z為平浮吃水時的坐標(biāo)系,X1-O1-Z1為縱傾Φ角度時的坐標(biāo)系,兩坐標(biāo)系變換關(guān)系為:

其中,x1,z1為縱傾Φ角度時的坐標(biāo),x,z為平浮吃水時的坐標(biāo)。

圖11～圖13分別為平浮吃水、艏傾角為Φ時和艉傾角為Φ時的網(wǎng)格劃分圖(三種情況下的排水量是相同的,所劃分的網(wǎng)格數(shù)相同)。

圖11 平浮吃水時的船體網(wǎng)格劃分圖

圖12 艏傾角為時的船體網(wǎng)格劃分圖

圖13 艉傾角為時的船體網(wǎng)格劃分圖

從圖11～圖13中可以看出,使用NAPA軟件二次開發(fā)語言——宏語言開發(fā)出的網(wǎng)格劃分模塊,可以迅速準(zhǔn)確地劃分出平浮吃水和縱傾狀態(tài)下艏、舯和艉部以及沿橫剖線方向上不同精度的船體網(wǎng)格。如某鋪管起重船的典型裝載工況一共有32個,其吃水在5m～8m之間,縱傾狀態(tài)包括了艏傾和艉傾,借助于NAPA軟件中NAPA BASIC下的循環(huán)語句For和條件判斷語句IF等函數(shù),以及LD SERVICE下的函數(shù)LD.QNT等,只需要在NAPA軟件中運行一次(不需要任何人工干預(yù)),在目前普通配置的微機中,不超過兩分鐘,則32個典型裝載工況的船體網(wǎng)格劃分文件已經(jīng)完成,其生成的文件可直接被HydroStar軟件調(diào)用。

4 結(jié) 語

通過本節(jié)運用NAPA軟件對船體網(wǎng)格快速劃分的方法的建立,可以看出,用這種方法劃分船體網(wǎng)格時非常方便、快速和準(zhǔn)確,劃分網(wǎng)格的精度可控。原則上可以劃分任意吃水下和任意縱傾下和任意精度的船體網(wǎng)格,本文只是對工程船駁船艏艉作了劃分,對于有球首和帶螺旋槳的尾部來說,這種劃分也是適用的,只是在球首和尾部線型分離處做分段處理(限于篇幅,這里就不再累述了),最后仍然可用本文的方法方便、快速和準(zhǔn)確地劃分出船體網(wǎng)格。本方法的建立,為工程船耐波性設(shè)計提供了必要的手段和便捷的工具。

[1] NAPA MANUALS[S].2009.2.

[2] BUREAU VERITAS.HYDROSTAR FOR EXPERTS USER MANUAL[R].March 2009.

[3] 李國勇.智能控制及其MATLAB實現(xiàn)[M].電子工業(yè)出版社,2005.

[4] 劉應(yīng)中,繆國平.船舶在波浪上的運動理論[M].上海交通大學(xué)出版社,1987.

[5] 張志軍.NAPA概述及NAPA船體模型的建立[J].船舶設(shè)計通訊,2003.

[6] 王彩蓮.淺談船舶應(yīng)用軟件NAPA[J].船舶,2003.

Automatic Grid Generation by NAPA for the Hull Wetted Surface of a Work boat

Sheng Qingwu

NAPA;bell function;longitudinal trim;automatic plotting

By using the macro-language,which is the secondary development tool in NAPA,the bell function and quadratic function in fuzzy ma the matics are introduced into the automatic grid distribution module for wetted surface of ships on even keel or in trim.This module can not only plot the grid rapidly and exactly,butal so generate the input file that can be called directly by Hydro Star,as well as control the output format so as to generate variou shull wetted surface grid formats that can be called by other hydrodynamic software.

U 674.3

A

1001-9855(2010)06-0053-05

2010-08-23

盛慶武(1971-),男,漢族,博士,高級工程師,主要從事船舶設(shè)計及流體性能計算研究工作。