杜文毅，高紅利

(廣東石油化工學(xué)院 機(jī)電工程學(xué)院，廣東 茂名 525000)

?

Workbench子模型技術(shù)在球罐應(yīng)力分析中的應(yīng)用*

杜文毅，高紅利

(廣東石油化工學(xué)院 機(jī)電工程學(xué)院，廣東 茂名 525000)

以3 000 m3液化烴球罐為例，對(duì)Workbench子模型技術(shù)在球罐應(yīng)力分析中的應(yīng)用進(jìn)行了討論。利用綁定接觸(基于MPC算法)形成過渡區(qū)域的方法在Workbench中實(shí)現(xiàn)了殼到體的子模型技術(shù)。根據(jù)粗糙模型的總體位移云圖確定了子模型的切割邊界，在子模型切割邊界上均勻地選取了8個(gè)點(diǎn)，通過對(duì)比粗糙模型網(wǎng)格加密前后這8個(gè)點(diǎn)的總體位移解確認(rèn)了切割邊界上位移解的合理性。對(duì)比結(jié)果可知，基于Workbench的子模型技術(shù)在占用較少計(jì)算資源的情況下能獲得精度較高的計(jì)算結(jié)果。

球罐；子模型；有限元；Workbench

0　引言

隨著石油化工工業(yè)的高速發(fā)展，球罐由于自身突出的優(yōu)點(diǎn)(和其他形式的壓力容器相比球罐占地面積更小、所需鋼材更少)得到了進(jìn)一步的應(yīng)用[1]。為了節(jié)約材料、降低成本采用基于有限元法的球罐分析設(shè)計(jì)已經(jīng)成為趨勢(shì)。對(duì)球罐的高應(yīng)力區(qū)域(如球殼和支腿的連接處、球殼和U型托板的連接處、U型托板和球殼的連接處等)的有限元處理是應(yīng)力分析中的重要環(huán)節(jié)。工程上為了提高高應(yīng)力區(qū)域應(yīng)力分析的準(zhǔn)確性，常用的方法是增加該區(qū)域及附近的網(wǎng)格密度[2]。由于球罐結(jié)構(gòu)的特點(diǎn)和載荷工況的復(fù)雜性，這種對(duì)高應(yīng)力區(qū)域的處理方法往往需要占用較多的計(jì)算資源和時(shí)間。而子模型技術(shù)可以在保證關(guān)鍵區(qū)域計(jì)算精度的同時(shí)有效的減少計(jì)算資源的占用。

文獻(xiàn)[3]～[5]利用ANSYS軟件中的子模型技術(shù)對(duì)球罐進(jìn)行了應(yīng)力分析。在此類分析中由于粗糙模型使用殼單元建模，子模型使用實(shí)體單元建模，在采用殼到體特殊子模型技術(shù)時(shí)就需要在ANSYS中進(jìn)行復(fù)雜的插值計(jì)算和文件操作。與ANSYS經(jīng)典版(Mechanical APDL)相比，ANSYS Workbench作為ANSYS公司集成了各種與仿真相關(guān)的API而提出的協(xié)同仿真環(huán)境，其界面友好、操作相對(duì)簡(jiǎn)單、更加適用于工程應(yīng)用。文章主要對(duì)利用Workbench子模型技術(shù)進(jìn)行球罐應(yīng)力分析進(jìn)行了討論。

1　子模型技術(shù)及其分析步驟

子模型技術(shù)又稱為切割邊界位移法，是一種從原模型(粗糙模型)上截取局部結(jié)構(gòu)(子模型)，并對(duì)該局部結(jié)構(gòu)加以細(xì)化操作(如網(wǎng)格加密或更改單元類型等)后進(jìn)行二次分析從而提高局部計(jì)算精度的一種方法。子模型技術(shù)在理論上基于圣維南原理。圣維南原理認(rèn)為彈性體表面的某部分的外力，可用作用在該部分上的等效力系來代替，這種代替只會(huì)使外力作用區(qū)域附近的應(yīng)力有顯著改變，但對(duì)于較遠(yuǎn)處的影響則可以忽略[5]。

在ANSYS Workbench軟件中子模型分析一般包含了下面幾個(gè)步驟。

1)建立粗糙模型。由于子模型的結(jié)果是根據(jù)粗糙模型在切割邊界上的位移解插值得到的，因此須保證粗糙模型具有足夠的網(wǎng)格密度以獲得合理的位移解。

圖1文件連接關(guān)系

2)建立文件連接關(guān)系。在Workbench軟件中靜力分析系統(tǒng)(static structure)一共包含了工程數(shù)據(jù)(Engineering Data)、幾何模型(Geometry)、有限元模型(Model)、設(shè)置(Setup)、解(Solution)、結(jié)果(Results)這六個(gè)單元。其中設(shè)置單元包含了求解所需的所有數(shù)據(jù)；有限元模型單元?jiǎng)t包含了網(wǎng)格、接觸、坐標(biāo)系和名稱選擇數(shù)據(jù)。建立子模型與粗糙模型的文件連接關(guān)系只需要把原靜力分析系統(tǒng)(粗糙模型靜力分析系統(tǒng))復(fù)制一份，然后連接原系統(tǒng)的解單元與復(fù)制系統(tǒng)的設(shè)置單元即可，具體如圖1所示。

3)建立子模型。首先清空(Reset)子模型系統(tǒng)的有限元模型單元，編輯其幾何模型(Geometry)單元即進(jìn)入DesignModeler(DM)。在DM中建立子模型切割邊界，并刪除邊界外的模型。然后編輯其有限元模型單元即進(jìn)入Mechanical，在Mechanical中對(duì)切割邊界加載位移約束(該位移約束文件由軟件根據(jù)切割邊界及子模型的網(wǎng)格自動(dòng)生成)。最后對(duì)子模型進(jìn)行網(wǎng)格劃分即可。

2　Workbench子模型技術(shù)在球罐應(yīng)力分析中的應(yīng)用

2.1球罐結(jié)構(gòu)參數(shù)及工況

以3 000 m3球罐為例進(jìn)行研究。其結(jié)構(gòu)參數(shù)，如表1所示。

表1　球罐結(jié)構(gòu)參數(shù)表

應(yīng)力分析涉及的工況為球罐的操作工況，即操作狀態(tài)重力載荷+計(jì)算壓力[7]。

2.2粗糙模型的建立與分析

由于載荷的對(duì)稱性使用1/10球罐模型進(jìn)行建模。粗糙模型采用殼單元shell181劃分網(wǎng)格，節(jié)點(diǎn)總數(shù)為4 395，單元總數(shù)為4 450。位移云圖及網(wǎng)格劃分如圖2所示。

為了使切割邊界上的位移解更加合理，應(yīng)盡量使切割邊界經(jīng)過位移響應(yīng)變化均勻的區(qū)域。故根據(jù)粗糙模型的位移云圖，可將圖3所示曲線在球殼上的投影設(shè)為子模型切割邊界。其中曲線圓弧部分的半徑為支腿內(nèi)徑的4倍(1 300 mm)，直線部分的長度為U型托板圓心到赤道面的距離(3 042 mm)。

圖2　粗糙模型的位移云圖和網(wǎng)格劃分　　　圖3　用于投影的空間曲線

為了確認(rèn)切割邊界位移解的正確性，可在切割邊界上均勻地取8個(gè)點(diǎn)，并將粗糙模型的網(wǎng)格密度增加一倍，最后對(duì)比網(wǎng)格加密前后這8個(gè)點(diǎn)的總體位移情況?？梢园l(fā)現(xiàn)粗糙模型網(wǎng)格加密后切割邊界上這8個(gè)點(diǎn)的位移偏差最大不超過3%，故認(rèn)為在此網(wǎng)格密度下切割邊界上的位移解是合理的[4]。選取點(diǎn)的位移情況對(duì)比，如圖4所示。

2.3子模型的建立與分析

Workbench并不支持由殼到體的子模型技術(shù)，為使子模型能在切割邊界上與粗糙模型連接，需要用殼單元shell181創(chuàng)建過渡區(qū)域。然后通過綁定接觸(基于MPC算法)連接過渡區(qū)域和遠(yuǎn)端區(qū)域(遠(yuǎn)端區(qū)域指子模型中遠(yuǎn)離切割邊界的區(qū)域，遠(yuǎn)端區(qū)域采用實(shí)體單元solid186建模)。遠(yuǎn)端區(qū)域與過渡區(qū)域的邊界生成方法與子模型切割邊界的類似，遠(yuǎn)端區(qū)域邊界圓弧部分的半徑變?yōu)橹葍?nèi)徑的3倍，直線部分不變。過渡區(qū)域和切割邊界上的位移約束，如圖5所示。

圖4　切割邊界上的點(diǎn)位移　　　　圖5　切割邊界位移約束和過渡區(qū)域

圖6等效應(yīng)力云圖和網(wǎng)格劃分

子模型包含了174 603個(gè)節(jié)點(diǎn)、58 872個(gè)單元。子模型遠(yuǎn)端區(qū)域的最大等效應(yīng)力為302.47 MPa，位置位于支腿和球殼連接處的頂部(粗糙模型得到的最大應(yīng)力為347.28 MPa，位置位于U型托板與球殼連接處的直線上)。子模型遠(yuǎn)端區(qū)域的等效應(yīng)力云圖和網(wǎng)格劃分情況，如圖6所示。

3　結(jié)語

(1)通過建模過程的對(duì)比可以發(fā)現(xiàn)，與ANSYS經(jīng)典界面相比，在Workbench中利用子模型技術(shù)進(jìn)行應(yīng)力分析要更加簡(jiǎn)便。

(2)對(duì)比最大應(yīng)力點(diǎn)的位置可以發(fā)現(xiàn)，在子模型中通過增加單元密度、改變單元類型可顯著的提高局部分析的精度。而如果為了在局部獲得同樣精度的解，直接地對(duì)整體模型進(jìn)行這種高密度的網(wǎng)格劃分將占用大量的計(jì)算資源。

[1] 徐英，楊一凡.球罐和大型儲(chǔ)罐[M].北京：化學(xué)工業(yè)出版社，2005.

[2] 萬興，張群，向玲.以Workbench進(jìn)行球罐設(shè)計(jì)的建模分析方案[J].化工裝備技術(shù)，2015，36(2)：35-40.

[3] 國宏斌，張景.ANSYS子模型技術(shù)在LPG球罐設(shè)計(jì)檢驗(yàn)中的應(yīng)用[J].化工機(jī)械，2009，36(6)：582-584.

[4] 高炳軍，高燕紅，李金紅.子模型法在球罐應(yīng)力分析中的應(yīng)用[J].壓力容器，2009，26(5)：27-31.

[5] 高艷紅，董俊華，高炳軍.子模型技術(shù)在大型球罐支撐區(qū)應(yīng)力分析及結(jié)構(gòu)優(yōu)化中的應(yīng)用[J].機(jī)械強(qiáng)度，2010，32(5)：735-739.

[6] 小魚兒和加加.JB4732-1995鋼制壓力容器-分析設(shè)計(jì)標(biāo)準(zhǔn)[EB/OL].(2012-08-28)[2016-05-11].http://wenku.baidu.com/link?url=SLJwyJGUZjcjGEYXHyQJshOrXKmCcb-dvoFxO5KLstsKM9r0pvNkFtere76PefZAT1ayLs4amBGeepR2iBXzjvnONORsZhUv0BrHA74gH3q.

[7] 中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局,中國國家標(biāo)準(zhǔn)化管理委員會(huì).GB12337—2014鋼制球形儲(chǔ)罐[S].北京：中國標(biāo)準(zhǔn)出版社,2015.

(責(zé)任編輯：黃容)

Stress Analysis of Spherical Tank based on Workbench Submodeling

DU Wenyi， GAO Hongli

(College of Machinery and Electronic Engineering, Guangdong University of Petrochemical Technology, Maoming 525000, China)

Taking 3 000 m3liquefied hydrocarbon spherical tank for example, the application of Workbench submodeling technology in spherical tank stress analysis is discussed. Bonded contact (based on the MPC algorithm) is used to form transition area to fulfill the submodeling technology from shell to solid in Workbench. Based on the total displacement contours of coarse model, the cutting boundary of the submodeling is determined, and 8 points on the cutting boundary is selected uniformly. The total displacements of the 8 points before and after mesh refining are compared so as to confirm the rationality of total displacement solution on cutting boundary. The results show that submodeling technique based on the Workbench can obtain better calculation results with less computational resources.

Spherical tank; Submodeling; Finite element; Workbench

2016-04-07；

2016-05-11

杜文毅(1987—)，男，江西南昌人，碩士，助教，主要從事結(jié)構(gòu)分析方面的科研工作。

TQ50.2

A

2095-2562(2016)04-0044-03