船舶結(jié)構(gòu)有限元模型局部修正方法研究

， ,，，

(1.中國(guó)船舶科學(xué)研究中心, 江蘇 無(wú)錫 214082;2.中國(guó)船級(jí)社, 北京 100007)

船體結(jié)構(gòu)強(qiáng)度的有限元校核多采用平行中體區(qū)域的三艙段模型。在進(jìn)行艙段的有限元計(jì)算后，需要對(duì)某些關(guān)鍵部件或區(qū)域使用更精細(xì)的網(wǎng)格進(jìn)行再次分析，而設(shè)計(jì)單位送審的模型，一般不考慮細(xì)部特征的建模，如開孔、趾端等。驗(yàn)船師在進(jìn)行規(guī)范校核時(shí)，通常的做法是根據(jù)初次計(jì)算的應(yīng)力結(jié)果，篩選出需要進(jìn)行開孔和趾端細(xì)化建模的區(qū)域，依據(jù)圖紙，在模型中手工創(chuàng)建開孔和趾端并細(xì)化，再進(jìn)行細(xì)化后的規(guī)范校核。這種方法適合針對(duì)典型區(qū)域的規(guī)范驗(yàn)證，但隨著HCSR規(guī)范的即將生效[1]，對(duì)細(xì)化網(wǎng)格的質(zhì)量要求更加嚴(yán)格，規(guī)范中要求校核的開孔區(qū)域和趾端數(shù)量也有所增加。對(duì)有限元模型進(jìn)行手工局部修正的方式已經(jīng)不適合新規(guī)范體系校核的需求。

此前依據(jù)CSR規(guī)范，討論過基于結(jié)構(gòu)有限元模型進(jìn)行開孔的方法，實(shí)現(xiàn)了借助硬點(diǎn)、硬線結(jié)合撒種的方法創(chuàng)建“雙圓弧形”開口[2]，某種程度上解決了工程實(shí)際中對(duì)典型區(qū)域抽樣校核的需要。但是隨著新規(guī)范的生效，開孔類型并不局限于雙圓弧形一種，需要校核的開孔區(qū)域范圍增多，并且對(duì)于趾端的細(xì)化和局部修正也有了明確的要求，因此原系統(tǒng)已不能適應(yīng)新的規(guī)范體系。所以，開發(fā)改良開孔細(xì)化以及趾端細(xì)化修正工具，能夠準(zhǔn)確地創(chuàng)建開孔和趾端，提高審圖驗(yàn)船人員的分析效率顯得尤為重要。

1 程序框架

針對(duì)開孔的創(chuàng)建方法、趾端細(xì)化以及局部建模展開討論，所提及的開孔創(chuàng)建基于局部細(xì)化后的模型[3]。根據(jù)設(shè)置的開孔和趾端的位置、尺寸參數(shù)等信息，在已有的船體結(jié)構(gòu)有限元模型上完成開孔的創(chuàng)建和趾端的局部細(xì)化和修正;將研究的成果基于Patran平臺(tái)，利用PCL語(yǔ)言編程進(jìn)行驗(yàn)證，并作為獨(dú)立的工具功能模塊，納入到CCS的HCSR規(guī)范校核軟件體系中。程序框架見圖1。

圖1 程序框架

2 開孔

開孔分為4種常用類型，分別為雙圓弧形、導(dǎo)角形、圓形和方形，多出現(xiàn)在雙層底實(shí)肋板和雙殼垂直桁，見圖2。

圖2 開孔類型和位置

參數(shù)設(shè)置界面及參數(shù)說(shuō)明以及雙圓弧形示例見圖3、表1。界面中給出規(guī)范中的開孔示意圖，并標(biāo)注相應(yīng)參數(shù)。根據(jù)開孔類型不同，界面示意圖及參數(shù)項(xiàng)做相應(yīng)的切換。

圖3 開孔功能界面

符號(hào)說(shuō)明o開孔所在平面的局部坐標(biāo)O開孔中心,Node或坐標(biāo)[x y z]r兩端半圓弧的半徑h0開孔的高度值l0開孔的長(zhǎng)度值,l0要求大于h0h開孔區(qū)域的高度值l開孔區(qū)域的長(zhǎng)度值

本文研究的開孔和趾端建模方法，均基于細(xì)化后網(wǎng)格創(chuàng)建，細(xì)化方法可利用獨(dú)立的通用網(wǎng)格細(xì)化方法[4]，將網(wǎng)格細(xì)化至100 mm×100 mm以內(nèi)，開孔區(qū)域不能有構(gòu)件與其相交，并且開孔區(qū)域單元必須在同一個(gè)平面內(nèi)。見圖4、5。

圖4 規(guī)范中的開孔圖示

圖5 網(wǎng)格細(xì)化

2.1 雙圓弧形

雙圓弧形開孔即是兩端為兩個(gè)等半徑的圓弧，形狀似足球場(chǎng)。創(chuàng)建步驟為確定開孔范圍,創(chuàng)建輔助建模特征及創(chuàng)建開孔處細(xì)化單元。

開孔范圍確定方法如下。

1)雙圓弧的一個(gè)圓弧圓心作為參考坐標(biāo)系o的原點(diǎn)。

2)圖6中開孔區(qū)域①為局部坐標(biāo)系o中，坐標(biāo)[x1+,y1+,0][x1-,y1-,0]覆蓋的矩形區(qū)域；

x1+=l0-r,x2+=l-r

(1)

y1+=h0/2,y2+=h/2

(2)

x1-= -r,x2-= -r-2Δ

(3)

y1-=-h0/2,y2-=-h0/2-2Δ

(4)

圖6 雙圓弧形開孔

3)開孔邊緣區(qū)域②為局部坐標(biāo)系o中，坐標(biāo)[x1++Δ,y1++Δ,0][x1--Δ,y1--Δ,0]覆蓋的矩形區(qū)域扣除區(qū)域①；規(guī)范中要求開孔區(qū)域的網(wǎng)格需要細(xì)化到50×50，開孔邊緣兩圈單元均勻并盡量不使用三角形單元，因此，取Δ=100。

4)開孔過渡區(qū)域③為局部坐標(biāo)系o中，坐標(biāo)[x2+,y2+,0][x2-,y2-,0]覆蓋的矩形區(qū)域扣除區(qū)域①和②。

5)需挖空的區(qū)域?yàn)棰?②+③，若單元的任何一個(gè)節(jié)點(diǎn)坐標(biāo)滿足式(5)、式(6)，即判定該單元落入需要挖空的區(qū)域。

NodePx={x|x∈(x2+,x2-)}

(5)

NodePy={y|y∈(y2+,y2-)}

(6)

2.1.1 輔助建模特征創(chuàng)建方法

1)根據(jù)節(jié)點(diǎn)的坐標(biāo)和單元拓?fù)潢P(guān)系，依據(jù)式(5)～(6)判定所在區(qū)域，得到要挖空的開孔附近單元。

2)刪除開孔區(qū)域單元。

3)以輪廓線的一個(gè)圓弧圓心為原點(diǎn)，x指向開孔長(zhǎng)度方向，y指向開孔高度方向，創(chuàng)建局部坐標(biāo)系。

4)創(chuàng)建開孔輪廓線b，為兩端半徑r的圓弧，兩端圓弧的端點(diǎn)由直線相連；

5)創(chuàng)建開孔邊緣區(qū)域輪廓線a，形狀同b，與開孔輪廓線之間的距離為100 mm。

2.1.2 創(chuàng)建開孔處單元方法

1)獲取挖空區(qū)域邊界的節(jié)點(diǎn)，可通過Free Edge方式獲取。

2)將x正方向節(jié)點(diǎn)按y值大小進(jìn)行排序。

3)根據(jù)圓弧長(zhǎng)度，單元尺寸要求為50×50，計(jì)算出需劃分的單元個(gè)數(shù)，并在輔助輪廓線上撒種。

4)根據(jù)規(guī)范中開孔單元特征要求，依次完成x正方向、x負(fù)方向、y正方向、y負(fù)方向的開孔單元和過渡單元的創(chuàng)建，見圖7、8。

圖7 雙圓弧形x方向創(chuàng)建過程

圖8 雙圓弧形y方向創(chuàng)建過程

5)刪除創(chuàng)建的幾何點(diǎn)、線和坐標(biāo)系等輔助元素。

2.2 導(dǎo)角形

導(dǎo)角形開孔即是矩形區(qū)域?yàn)?個(gè)頂點(diǎn)作等半徑圓弧導(dǎo)角，創(chuàng)建步驟同雙圓弧形。局部坐標(biāo)系取在一端兩個(gè)導(dǎo)角圓心連線的中點(diǎn)，開孔區(qū)域的判定類似式(1)～(4)，輔助建模特征創(chuàng)建和網(wǎng)格創(chuàng)建的方法與雙圓弧形開孔類似。見圖9～11。

圖9 導(dǎo)角形開孔示意圖

圖10 導(dǎo)角形創(chuàng)建過程

圖11 導(dǎo)角形開孔效果

2.3 圓形

圓形開孔即是開孔輪廓為圓，創(chuàng)建步驟同雙圓弧形。局部坐標(biāo)系取在圓心，開孔區(qū)域即是圓外接矩形覆蓋的區(qū)域，輔助建模特征創(chuàng)建和網(wǎng)格創(chuàng)建的方法均與雙圓弧形開孔類似，其中過渡區(qū)域的三角區(qū)域依據(jù)規(guī)范中示意的網(wǎng)格特征，計(jì)算節(jié)點(diǎn)位置，并連接成單元。

圖形開孔示意見圖12～14。

2.4 方形

方形開孔即是開孔輪廓為矩形。直接刪除開孔覆蓋的矩形區(qū)域內(nèi)單元形成開孔的方法，若孔的邊界恰巧位于單元上，而不是單元的邊上，則會(huì)造成開孔尺寸有誤差，為避免該問題，矩形開孔的創(chuàng)建分為確定開孔范圍、創(chuàng)建輔助建模特征、創(chuàng)建開孔邊緣區(qū)域兩圈細(xì)化單元，創(chuàng)建過渡區(qū)細(xì)化單元幾個(gè)部分。

圖12 圓形開孔示意

圖13 圓形開孔X方向創(chuàng)建過程

圖14 圓形開孔Y方向創(chuàng)建過程

1)局部坐標(biāo)系取在矩形輪廓的型心，開孔區(qū)域即是[l/2+100，h/2+100，0]和[-l/2-100，-h/2-100，0]覆蓋的矩形區(qū)域。根據(jù)開孔的尺寸創(chuàng)建矩形輪廓線，以100 mm的間距創(chuàng)建4條直線作為輔助元素,見圖15。

圖15 方形開孔示意

2)將開孔區(qū)域的長(zhǎng)和高按近似50 mm等分，用等分?jǐn)?shù)在輪廓線上撒種，連接開孔邊緣區(qū)域內(nèi)兩圈單元,見圖16。

圖16 方形開孔邊緣區(qū)域

3)將過渡區(qū)域分為8個(gè)小區(qū)域進(jìn)行處理，

區(qū)域1的處理方法為：取與過渡直線x負(fù)方向最近的單元，在過渡邊上撒上種子，完成區(qū)域1的單元?jiǎng)?chuàng)建。

區(qū)域2的處理方法為：將過渡直線與y正向直線撒上種子，完成區(qū)域2處的單元?jiǎng)?chuàng)建。

區(qū)域3的處理方法為：取與過渡直線x正方向的最近的單元，將邊撒上種子，完成區(qū)域3的單元?jiǎng)?chuàng)建。見圖17～19。

其他區(qū)域按位置與區(qū)域1、2、3處理方法類似。

圖17 方形開孔過渡區(qū)域劃分

圖18 方形開孔過渡區(qū)創(chuàng)建過程

圖19 方形開孔效果

4)刪除創(chuàng)建的幾何點(diǎn)、線和坐標(biāo)系。

3 趾端建模

趾端分為曲線型(二次曲線，默認(rèn)取曲率為0.5)和圓弧形兩種常用類型，多出現(xiàn)在強(qiáng)框架和艙壁交界處。參數(shù)設(shè)置界面及參數(shù)說(shuō)明以及曲線形示例見圖20，界面中給出規(guī)范中的示意圖，并標(biāo)注相應(yīng)參數(shù)。根據(jù)類型不同，界面示意圖及參數(shù)項(xiàng)作相應(yīng)的切換。見表2。

圖20 趾端建模界面示意

符號(hào)說(shuō)明o趾端所在平面的局部坐標(biāo)b2趾端高度Ra趾端長(zhǎng)度s趾端延長(zhǎng)上部水平長(zhǎng)度Rb趾端延長(zhǎng)下部水平長(zhǎng)度h趾端延長(zhǎng)右側(cè)高度c趾端延長(zhǎng)部分相對(duì)坐標(biāo)原點(diǎn)的位置Kp曲線的曲率α收縮角度

趾端細(xì)化原理流程如下。

1)獲取梁?jiǎn)卧?，針?duì)需要收縮的梁?jiǎn)卧M(jìn)行板元化。

①獲取所有自由邊上的梁?jiǎn)卧?根據(jù)位置信息過濾掉非趾端上的梁?jiǎn)卧?/p>

②獲取梁?jiǎn)卧系墓?jié)點(diǎn)，根據(jù)節(jié)點(diǎn)順序確定板元化后單元的大小。

③獲取節(jié)點(diǎn)最小處對(duì)應(yīng)的梁?jiǎn)卧?，獲取該單元的截面屬性。

④獲取單元截面屬性，梁?jiǎn)卧獎(jiǎng)?chuàng)建成板單元。

⑤反向板元化一個(gè)梁?jiǎn)卧?/p>

2)粗網(wǎng)格細(xì)化。

3)創(chuàng)建趾端。

①獲取細(xì)化范圍中x最大的兩個(gè)節(jié)點(diǎn)。

②根據(jù)獲取到的節(jié)點(diǎn)與參數(shù)創(chuàng)建趾端網(wǎng)格。

4)移動(dòng)節(jié)點(diǎn)，網(wǎng)格修正。趾端模型見圖21、圖22。

5)刪除輔助元素。

圖21 趾端模型示意

圖22 趾端建模效果

4 結(jié)論

本文開發(fā)的功能模塊適用于HCSR規(guī)定的一道或者兩道縱艙壁雙殼油船、單舷側(cè)或雙舷側(cè)散貨船結(jié)構(gòu)強(qiáng)度校驗(yàn)中的模型局部修正，已經(jīng)作為一個(gè)功能模塊納入CCS的HCSR規(guī)范校核軟件體系中。實(shí)測(cè)證明，該程序能夠快速準(zhǔn)確地創(chuàng)建開孔和趾端，在一定程度上減輕審圖驗(yàn)船人員的繁瑣重復(fù)勞動(dòng)，提高審圖效率。

[1] IACS.雙殼油船和散貨船協(xié)調(diào)共同結(jié)構(gòu)規(guī)范(HCSR)[J].船舶標(biāo)準(zhǔn)化工程師,2012,45(5):16-16.

[2] 單威俊,伊金秀,李 鋒,等.艙口角隅處網(wǎng)格劃分研究[J].船海工程,2008,37(1):10-14.

[3] 中國(guó)船級(jí)社.鋼質(zhì)海船入級(jí)規(guī)范[S],北京：人民交通出版社,2009.

[4] 單威俊.船舶結(jié)構(gòu)分析中網(wǎng)格劃分技術(shù)及應(yīng)用[D].無(wú)錫：江南大學(xué),2008.