，

(華南理工大學(xué) 交通學(xué)院，廣州 510641)

現(xiàn)階段船體結(jié)構(gòu)強(qiáng)度的直接計算一般運(yùn)用通用有限元軟件。其中，有限元模型的建模主要是由手工建模完成，這種建模方法的效率及重復(fù)利用率較低，直接應(yīng)用通用有限元軟件建模需要較高的有限元應(yīng)用技巧且較費(fèi)時，拉長了分析周期。這種建模方式不能滿足現(xiàn)實(shí)船舶設(shè)計建造的要求，也不具備處理突發(fā)事件的能力。

因此，提出參數(shù)化建模的方法，所謂參數(shù)化建模，也即自動化建模，是指由用戶在界面上輸入必要的參數(shù)，程序根據(jù)這些參數(shù)進(jìn)行分析和計算，最終完成模型創(chuàng)建作為輸出。整個建模的主要過程由程序完成，期間只需少量的用戶操作。用參數(shù)化建模的方法來實(shí)現(xiàn)船體有限元模型的自動建模，能顯著地縮短船體結(jié)構(gòu)有限元分析的周期，提高效率。

在船體結(jié)構(gòu)參數(shù)化建模的探索中，國內(nèi)外的許多船舶工作者都為之作出了不懈的努力[1-4]。

在綜合前人研究的基礎(chǔ)上，筆者提出了一個參數(shù)化建模研究的思路。

1 研究思路

用PCL語言編寫程序代碼，并導(dǎo)入MSC.Patran中，同時充分利用MSC.Patran的功能，最終實(shí)現(xiàn)三維船體艙段的有限元自動建模。其中，難點(diǎn)與關(guān)鍵技術(shù)在于有限元模型自動生成方法的設(shè)計，即由程序自動生成高質(zhì)量的有限元網(wǎng)格的方法。

船體艙段結(jié)構(gòu)是由桿、梁、板、殼等構(gòu)件組成的彈性體，是一個大型空間薄壁板梁組合結(jié)構(gòu)。根據(jù)船體艙段結(jié)構(gòu)的特點(diǎn)，可以預(yù)先對其有限元模型的網(wǎng)格進(jìn)行規(guī)劃。

首先，從單元類型上講，船體艙段有限元模型所用的單元主要是一維單元和二維單元。其中，一維單元主要是2節(jié)點(diǎn)的梁單元和2節(jié)點(diǎn)的桿單元，二維單元主要是3節(jié)點(diǎn)的三角形殼單元和4節(jié)點(diǎn)的四邊形殼單元。

其次，從船體艙段的結(jié)構(gòu)特征看，可以把船體艙段上的單元?dú)w類為縱向單元和橫向單元，或是把整個網(wǎng)格分為縱向網(wǎng)格和橫向網(wǎng)格兩部分。其中，把沿著船長方向的網(wǎng)格稱為縱向網(wǎng)格，其上的單元即縱向單元；把全部橫向構(gòu)件上的網(wǎng)格稱為橫向網(wǎng)格，其上的單元即橫向單元。

先通過輸入的信息直接生成縱向有限元網(wǎng)格，再以縱向網(wǎng)格為基礎(chǔ)，在相應(yīng)的肋位上利用已有的節(jié)點(diǎn)作為邊界，進(jìn)而生成橫向的有限元網(wǎng)格。

1.1 縱向網(wǎng)格的自動生成方法

遵循一定的規(guī)則輸入若干個典型橫剖面上的插值點(diǎn)的插值信息(如點(diǎn)的幾何信息，選用的插值類型，等等)，程序根據(jù)輸入信息用幾何線構(gòu)造一個空間的幾何骨架，然后在此幾何骨架的基礎(chǔ)上生成縱向的有限元網(wǎng)格。

用分段的三次非均勻B樣條曲線對船體外殼或甲板在每個典型橫剖面上的輸入點(diǎn)進(jìn)行插值，并根據(jù)實(shí)際建模的需要對所構(gòu)造的三次非均勻B樣條曲線進(jìn)行插值計算以獲取建模所需的其它點(diǎn)[5]。在每個剖面上對全部的插值點(diǎn)進(jìn)行連線，然后對相鄰剖面上的插值點(diǎn)進(jìn)行連線，這2組線構(gòu)成了一個空間的幾何骨架，最后，在此幾何骨架上生成縱向的有限元網(wǎng)格[6，7]。

對于縱向構(gòu)件的其它部分，如內(nèi)殼、縱艙壁等，也用相同的方法完成，所不同的是這些部分在橫剖面上的投影表現(xiàn)為直線段，因此只需進(jìn)行線性插值計算而無需用樣條曲線進(jìn)行插值。

縱向網(wǎng)格的整個創(chuàng)建過程如圖1所示。

圖1 縱向網(wǎng)格創(chuàng)建的流程

縱向網(wǎng)格生成的同時，在程序中根據(jù)輸入的信息對各部分構(gòu)件進(jìn)行分組，以便對模型進(jìn)行檢查以及進(jìn)行各部分構(gòu)件屬性賦值的操作。從總體上看，縱向網(wǎng)格生成的整個過程是一個由點(diǎn)連線，進(jìn)而在由線所構(gòu)成的幾何骨架上生成縱向網(wǎng)格的過程，在這一過程中，由幾何元素所承載的信息最終傳遞到有限元網(wǎng)格的單元之中?？v向網(wǎng)格的生成也為接下的橫向網(wǎng)格的生成奠定了基礎(chǔ)。

選取一個雙底雙殼的三維船體艙段結(jié)構(gòu)進(jìn)行分析，其縱向網(wǎng)格的生成如圖2～4所示。

1.2 橫向網(wǎng)格的自動生成方法

對于三維船體艙段上的橫向網(wǎng)格，由于船體橫向構(gòu)件所在的區(qū)域通常是一些細(xì)帶狀的二維區(qū)域，且其邊界上的節(jié)點(diǎn)的位置和個數(shù)一般都已預(yù)先確定，因此，盡管在二維區(qū)域上生成有限元網(wǎng)格的方法有多種，但要在這些細(xì)帶狀的二維區(qū)域中自動生成高質(zhì)量的二維網(wǎng)格仍然困難。

為了解決此問題，在現(xiàn)有的網(wǎng)格生成理論基礎(chǔ)上提出了一種新的二維有限元網(wǎng)格的鋪設(shè)算法，用于在細(xì)帶狀二維區(qū)域自動生成包含三角形單元和四邊形單元的有限元網(wǎng)格，其中，四邊形單元占主要部分，而少量的三角形單元則是為了加強(qiáng)算法對邊界的適應(yīng)性。

基本原理主要基于Blacker和文獻(xiàn)[8]和[9]。同時，筆者提出一些新的算法規(guī)則，使算法所生成的單元更能有效適應(yīng)細(xì)帶狀的二維區(qū)域。

圖5 橫向網(wǎng)格生成算法的流程圖

主要步驟包括單元和節(jié)點(diǎn)的生成、相交、縫合、把網(wǎng)格節(jié)點(diǎn)移至內(nèi)部固定節(jié)點(diǎn)、網(wǎng)格局部調(diào)整、均勻化共6個部分。算法的主要過程見圖5。

以大型綜合有限元軟件MSC.Patran為應(yīng)用平臺，用MSC.Patran提供給用戶的二次開發(fā)語言PCL實(shí)現(xiàn)了該算法，所得到的有限元網(wǎng)格的質(zhì)量優(yōu)于直接用MSC.Patran自身的Paver命令劃分所得到的非結(jié)構(gòu)化網(wǎng)格[7]。見圖6，這是對于一個艙段肋位上的強(qiáng)肋骨，用本算法對所選定的的邊界節(jié)點(diǎn)以及內(nèi)部固定節(jié)點(diǎn)所進(jìn)行的網(wǎng)格自動生成的結(jié)果。其中，外部邊界上的節(jié)點(diǎn)由縱向網(wǎng)格在指定位置上的節(jié)點(diǎn)所確定，內(nèi)部的固定節(jié)點(diǎn)是根據(jù)實(shí)際需要而添加的。從圖中可以看出，程序可以自動生成網(wǎng)格的主要部分。當(dāng)然，用算法所生成的網(wǎng)格還需要后續(xù)的一些人為的修改與調(diào)整工作。

a)

b)

c)

2 程序的數(shù)據(jù)結(jié)構(gòu)圖

程序中主要的數(shù)據(jù)結(jié)構(gòu)見圖7。

圖7 主要的數(shù)據(jù)結(jié)構(gòu)

3 自動建模模塊的界面

根據(jù)本文所提出的參數(shù)化建模思路，三維船體艙段有限元模型自動生成所依賴的交互界面見圖8～9。

圖8 自動生成縱向網(wǎng)格的菜單選項(xiàng)以及命令框

圖9 鋪設(shè)方法的菜單選項(xiàng)以及命令框

4 結(jié)論與展望

對于一般的油船或散貨船，用這種參數(shù)化建模方法可以由程序生成整個三維船體艙段有限元模型中接近70%的組成部分，甚至更多，能夠減少接近50%的建模時間。當(dāng)然，本課題目前只處于研究階段，所完成的只是基礎(chǔ)部分的工作，對于程序生成的網(wǎng)格還需要一定的人為調(diào)整。整個參數(shù)化建模系統(tǒng)的完成將會進(jìn)一步提高建模的自動化程度。

[1] Forde B W R, Foschi R O, Stiemer S F.Objectoriented finite element analysis[J].Comp & Struct, 1990, 34(3): 355-374.

[2] 羅金炎, 陳慶強(qiáng).船舶面向?qū)ο笥邢拊膽?yīng)用研究[J].計算機(jī)輔助工程, 2004(1): 18-22.

[3] 張偉星, 徐 立, 趙成璧，等.基于MSC.Acumen和PCL開發(fā)船舶快速建模系統(tǒng)[J].CAD/CAM與制造業(yè)信息化，2001: 81-82.

[4] In-II Kim.A development of data structure and mesh generation algorithm for whole ship analysis modeling system[J].Advances in Engineering Software，2006,37: 85-96.

[5] 施法中.計算機(jī)輔助幾何設(shè)計與非均勻有理B樣條(CAGD & NURBS)[M].北京：北京航空航天大學(xué)出版社, 1994.

[6] 林 慰, 趙成璧.基于MSC.Patran的三維船體有限元混合網(wǎng)格自動剖分方法的研究[J].船海工程, 2005，10(特刊): 51-54.

[7] 林 慰.船體艙段有限元模型自動建模方法的研究[D].廣州：華南理工大學(xué), 2006.

[8] Blacker, Ted D Michael B.Stephenson.paving: A new approach to automated quadrilateral mesh generation[J].International Journal for Numerical Methods in Engineering,1991, 32: 811-847.

[9] White, David R Paul Kinney.Redesign of the paving algorithm: Robustness enhancements through element by element meshing[J], Proceedings.6th International Meshing Roundtable.Sandia National Laboratories,1997: 323-335.

[10] 中國船級社.散貨船結(jié)構(gòu)強(qiáng)度直接計算指南[M].北京：人民交通出版社， 2003.

[11] 中國船級社.油船結(jié)構(gòu)強(qiáng)度直接計算指南[M].北京：人民交通出版社，2003.