李海靜，陳有芳，王麗榮，唐宜健，胡豐梁，章志兵

（1.華中科技大學(xué) 材料科學(xué)與工程學(xué)院，武漢 430074；2.中國(guó)船級(jí)社技術(shù)研發(fā)中心，北京 100007）

疲勞失效是影響船舶壽命的重要因素，也是船舶結(jié)構(gòu)設(shè)計(jì)中不容小覷的部分[1]。目前普遍采用簡(jiǎn)化計(jì)算的方式來(lái)校核船舶結(jié)構(gòu)強(qiáng)度，按照評(píng)估范圍不同，可分為局部結(jié)構(gòu)計(jì)算、艙段結(jié)構(gòu)計(jì)算與全船結(jié)構(gòu)計(jì)算[2—6]。校核時(shí)，面對(duì)數(shù)量龐大、截面形狀尺寸多變、方向各異的筋結(jié)構(gòu)，設(shè)計(jì)人員往往采用抽樣法進(jìn)行筋校核。這種方式受到計(jì)算數(shù)量、樣本分布的制約，經(jīng)驗(yàn)依賴度高，計(jì)算工作量大，效率低下，可靠性差?；诖耍闹刑岢隽艘惶状胺鰪?qiáng)材自動(dòng)理想化系統(tǒng)，實(shí)現(xiàn)對(duì)分布密集、復(fù)雜多樣的船舶扶強(qiáng)材的批量自動(dòng)規(guī)范計(jì)算，降低校核的人工依賴性，提高校核準(zhǔn)確性，從而提高船舶設(shè)計(jì)效率，保障設(shè)計(jì)質(zhì)量。

1 扶強(qiáng)材理想化系統(tǒng)體系結(jié)構(gòu)

文中所提出的船舶扶強(qiáng)材自動(dòng)理想化系統(tǒng)隸屬于船舶工程計(jì)算軟件Compass 3D 集成系統(tǒng)。該系統(tǒng)集CAD 建模模塊[7—9]、SDR 數(shù)據(jù)管理模塊[10]、CAE有限元分析模塊[11—12]、SDP 規(guī)范計(jì)算模塊[13]于一體，設(shè)計(jì)人員可以一站式完成船體建模、分析、校核的全過(guò)程，無(wú)需處理各子系統(tǒng)間的過(guò)渡，操作簡(jiǎn)潔，用戶友好度高。

規(guī)范計(jì)算子系統(tǒng)用于部分三維模型建立后。在規(guī)范計(jì)算子系統(tǒng)中，板結(jié)構(gòu)與筋結(jié)構(gòu)被離散為子板與子筋單元，重組后構(gòu)成基本計(jì)算單元。模型示意圖如圖1 所示。

圖1 模型示意圖Fig.1 Model diagram

2 扶強(qiáng)材理想化系統(tǒng)功能設(shè)計(jì)

系統(tǒng)功能設(shè)計(jì)見(jiàn)圖2。

2.1 單元梁快速自動(dòng)識(shí)別

圖2 系統(tǒng)功能設(shè)計(jì)Fig.2 Design of system function

扶強(qiáng)材，也稱為加強(qiáng)筋，是為了增強(qiáng)船體金屬板結(jié)構(gòu)的剛度和強(qiáng)度所引入的鋼質(zhì)船體結(jié)構(gòu)材料[14]。在船體三維模型中，加強(qiáng)筋依附于板架，由理論線、復(fù)線和筋三維模型組成，如圖3 所示。其中，復(fù)線用于規(guī)范計(jì)算；理論線反應(yīng)了加強(qiáng)筋在CAD 模型中的空間分布；筋三維模型則直觀地呈現(xiàn)出該加強(qiáng)筋的空間形態(tài)，包括截面形狀、截面尺寸、依附狀況、傾斜角度等。

圖3 加強(qiáng)筋模型描述Fig.3 Description of stiffener model

在經(jīng)過(guò)預(yù)處理的目標(biāo)船體模型中，原始的整段加強(qiáng)筋經(jīng)過(guò)一定規(guī)則的切割，被離散為多段子筋，每個(gè)子筋都繼承了其父筋的相關(guān)屬性。將共用端點(diǎn)且該端點(diǎn)未終止于結(jié)構(gòu)相交處的子筋集合記錄為一個(gè)筋計(jì)算單元，稱為單元梁，如圖4 所示。

圖4 單元梁模型描述Fig.4 Description of unit beam model

單元梁識(shí)別的總體流程如圖5 所示。①步驟1為初始化子筋相交關(guān)系，提取分析范圍內(nèi)的全部子筋，按照依附板架對(duì)子筋進(jìn)行分組，以分組為單位，對(duì)子筋、子筋端點(diǎn)進(jìn)行編號(hào)，并記錄子筋之間的相交關(guān)系；② 步驟2 為識(shí)別可構(gòu)成單元梁的子筋集合，以初始化分組為分析單位，首先任意獲取一個(gè)子筋（Curve1），任取一側(cè)端點(diǎn)（Point11），根據(jù)初始化相交關(guān)系，搜索到與其相連的子筋（Curve2），并判斷相連子筋的非共用端點(diǎn)（Point22）是否與板架結(jié)構(gòu)相交，如果相交，則該端終止，轉(zhuǎn)到初始子筋的另一側(cè)端點(diǎn)（Point12）繼續(xù)進(jìn)行上述操作，如果相連子筋的非共用端點(diǎn)（Point22）不與板架結(jié)構(gòu)相交，則以該點(diǎn)作為基點(diǎn)，繼續(xù)搜索相連子筋并進(jìn)行板架相交性判斷，直至非共用端點(diǎn)與板架相交或者已經(jīng)到達(dá)板邊界，如圖6 所示，在這個(gè)過(guò)程中，如果通過(guò)一個(gè)端點(diǎn)可以搜索到多個(gè)相連子筋，則根據(jù)相連子筋與用于搜索的上層子筋的夾角進(jìn)行篩選，選擇夾角最接近平角的子筋作為唯一的目標(biāo)子筋；③步驟3 為搜索到的子筋集合記錄相同的標(biāo)識(shí)，構(gòu)成一個(gè)單元梁，采用自定義數(shù)據(jù)對(duì)象進(jìn)行存儲(chǔ)，當(dāng)分析范圍內(nèi)所有的子筋都加入且只加入到唯一的單元梁后，即完成了單元梁的快速自動(dòng)識(shí)別，通過(guò)單元梁的快速自動(dòng)識(shí)別，將離散后的加強(qiáng)筋重新組合簡(jiǎn)化為若干單元梁，從而實(shí)現(xiàn)相關(guān)規(guī)范參數(shù)的快速準(zhǔn)確計(jì)算。

圖5 單元梁自動(dòng)識(shí)別流程Fig.5 Automatic recognition flowchart of unit beam

圖6 識(shí)別相連曲線示意圖Fig.6 Identifying connected curves

2.2 船體梁結(jié)構(gòu)強(qiáng)度校核參數(shù)分析

船體梁結(jié)構(gòu)強(qiáng)度校核分析是以單元梁為單位對(duì)船體局部強(qiáng)度模型進(jìn)行相關(guān)校核。對(duì)于加強(qiáng)筋而言，需要校核的參數(shù)包含壓力載荷計(jì)算點(diǎn)、彎曲應(yīng)力載荷計(jì)算點(diǎn)、單元梁間距、帶板與水平面夾角、單元梁兩側(cè)的艙室環(huán)境等共13 項(xiàng)參數(shù)。下面以單元梁間距為例，進(jìn)行簡(jiǎn)要介紹。

在船舶結(jié)構(gòu)中，存在許多大跨度加筋板結(jié)構(gòu)，比如甲板、橫向艙壁等，板架上加強(qiáng)筋的分布、尺寸對(duì)板架強(qiáng)度的影響很大[15]。單元梁間距是單元梁與其鄰近兩側(cè)單元梁的垂直距離。由于船舶三維模型在轉(zhuǎn)換過(guò)程中經(jīng)過(guò)了多次切割，剖面模型中存在許多碎邊、懸臂邊、不閉合區(qū)域，如圖7 所示。如果采用直接搜索兩側(cè)單元梁求解的方式，計(jì)算結(jié)果和計(jì)算速度都會(huì)受到影響。在校核誤差允許的條件下，最終采用等效簡(jiǎn)化的方式來(lái)計(jì)算單元梁間距。

圖7 間距示意圖Fig.7 Pitch diagram

據(jù)單元梁端部形式，將單元梁分為屈曲梁和非屈曲梁兩類。其中，由非屈曲梁、板相交縫圍成的區(qū)域稱作屈服板格；由屈服梁、屈曲梁和板相交縫圍成的區(qū)域稱作屈曲板格。采用矩形等效的方法，將非矩形板格等效為矩形板格，如圖8 所示。分別在單元梁B1B2兩端（B1，B2）搜索板格，并獲取到非矩形板格1 和板格2 的面積，用板格面積除以單元梁跨距即可計(jì)算出對(duì)應(yīng)的等效間距。

圖8 板格矩形等效示意圖Fig.8 Rectangle equivalent diagram of panel

采用矩形等效法簡(jiǎn)化后的理想化板格如圖9 所示。對(duì)每個(gè)結(jié)構(gòu)板而言，只需要進(jìn)行一次板格分組，即可計(jì)算出整個(gè)板架上所有單元梁之間的間距。

圖9 理想化板格示意圖Fig.9 Diagram of idealized gird

3 應(yīng)用實(shí)例

基于NX11.0 平臺(tái)，使用C++語(yǔ)言開發(fā)上述船舶扶強(qiáng)材自動(dòng)理想化系統(tǒng)，并在以油船和散貨船為代表的船型三維模型上進(jìn)行測(cè)試。測(cè)試結(jié)果表明，在不同縱向位置處，扶強(qiáng)材自動(dòng)理想化系統(tǒng)均可快速、準(zhǔn)確地輸出剖面中所有筋的校核參數(shù)分析結(jié)果。

圖10 是180000DWT 散貨船F(xiàn)R176 位置處剖面模型的扶強(qiáng)材自動(dòng)理想化系統(tǒng)計(jì)算結(jié)果示意圖。其中，模型的結(jié)構(gòu)建模公差為2.0 mm。在Intel(R) Core(TM) i5-6400、CPU 2.70 GHz、8 G 內(nèi)存的測(cè)試機(jī)上測(cè)試發(fā)現(xiàn)，剖面模型中共220 個(gè)單元梁、每個(gè)單元梁13 項(xiàng)參數(shù)的計(jì)算過(guò)程共耗時(shí)約15 s。

圖10 180000DWT 散貨船F(xiàn)R161 剖面自動(dòng)理想化計(jì)算結(jié)果Fig.10 Automatic idealization calculation result of FR161 section of 180000DWT bulk carrier

4 結(jié)論

主要研究了基于三維幾何模型的船舶扶強(qiáng)材自動(dòng)理想化系統(tǒng)。本系統(tǒng)可對(duì)船舶三維模型中分布密集、截面多變、依附情況復(fù)雜、人工難以計(jì)算的扶強(qiáng)材規(guī)范校核參數(shù)進(jìn)行快速自動(dòng)的理想化計(jì)算。結(jié)果表明，本系統(tǒng)可以快速準(zhǔn)確地計(jì)算出船舶剖面模型中所有扶強(qiáng)材的壓力載荷計(jì)算點(diǎn)、彎曲應(yīng)力載荷計(jì)算點(diǎn)、間距、兩側(cè)的艙室環(huán)境等共13 項(xiàng)規(guī)范校核參數(shù)，解決了人工校核的多項(xiàng)不便，有效提高了船舶加強(qiáng)筋結(jié)構(gòu)校核的精度和效率。