許銳敏，潘 駿，陳文華，王 孟，許成彬

（浙江理工大學(xué)浙江省機(jī)電產(chǎn)品可靠性技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，杭州310018）

筒子染色機(jī)染缸碟形封頭的可靠性優(yōu)化設(shè)計(jì)

許銳敏，潘 駿，陳文華，王 孟，許成彬

（浙江理工大學(xué)浙江省機(jī)電產(chǎn)品可靠性技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，杭州310018）

染缸的碟形封頭對筒子染色機(jī)使用的安全可靠性有著極其重要的影響。以筒子染色機(jī)染缸碟形封頭為研究對象，采用ANSYS Workbench軟件對其進(jìn)行參數(shù)化建模和有限元分析；通過靈敏度分析得到了影響其等效應(yīng)力和總變形的關(guān)鍵參數(shù)，利用實(shí)驗(yàn)數(shù)據(jù)法和多目標(biāo)驅(qū)動(dòng)優(yōu)化方法對封頭的結(jié)構(gòu)尺寸進(jìn)行優(yōu)化設(shè)計(jì)；建立考慮應(yīng)力和變形相關(guān)性的可靠性計(jì)算模型，運(yùn)用6西格瑪設(shè)計(jì)準(zhǔn)則和拉丁超立方體抽樣方法，得到了最大等效應(yīng)力和最大總變形的概率分布，得出了優(yōu)化前后封頭的可靠度，并且優(yōu)化后可靠度有明顯提高。

染缸；碟形封頭；ANSYS Workbench；靈敏度分析；可靠性優(yōu)化設(shè)計(jì)

0 引 言

筒子染色機(jī)是紗線染色機(jī)中的一種，通過循環(huán)泵控制染液在染缸內(nèi)循環(huán)流動(dòng)，對紗錠進(jìn)行染色。筒子染色機(jī)的主要結(jié)構(gòu)包括染缸、循環(huán)泵、熱交換裝置和控制系統(tǒng)等，而其中染缸是由缸蓋和缸身組成的封閉壓力容器，由于惡劣的工作環(huán)境，其使用的安全可靠性引起了廣泛的關(guān)注。碟形封頭又稱帶折邊的球形封頭，是壓力容器凸形封頭的一種，由于制造簡單而被廣泛應(yīng)用于中、低壓壓力容器。隨著染整行業(yè)的不斷發(fā)展，碟形封頭的需求量越來越大，對碟形封頭使用的安全可靠性要求也日益提高，人們也逐步開始重視碟形封頭在參數(shù)不確定情況下的可靠性研究。

傳統(tǒng)的壓力容器封頭設(shè)計(jì)大多采用常規(guī)的安全系數(shù)設(shè)計(jì)法［1］，由于常規(guī)設(shè)計(jì)未考慮影響因素的不確定性情況，因此許多學(xué)者對壓力容器封頭可靠性設(shè)計(jì)進(jìn)行了研究。胡宜生［2］假設(shè)各設(shè)計(jì)參數(shù)都服從正態(tài)分布，利用靜力學(xué)建立應(yīng)力的計(jì)算公式，并得到其概率分布，從而利用應(yīng)力-強(qiáng)度干涉理論計(jì)算其可靠度。農(nóng)琪［3］以應(yīng)力-強(qiáng)度干涉理論為基礎(chǔ)對壓力容器橢圓形封頭進(jìn)行了可靠性設(shè)計(jì)。周鵬飛［4］利用ANSYS軟件的可靠性分析功能和優(yōu)化設(shè)計(jì)功能對異形封頭的可靠性分析和優(yōu)化進(jìn)行了研究。彭翠玲［5］運(yùn)用ANSYS的概率設(shè)計(jì)功能，以壓力容器的壁厚、壓力載荷及彈性模量為隨機(jī)輸入變量，利用蒙特卡羅抽樣方法，獲得了壓力容器封頭有限元分析模型的應(yīng)力概率分布特征。上述方法僅考慮了應(yīng)力和強(qiáng)度的隨機(jī)性，未考慮變形的影響及其與應(yīng)力之間的相關(guān)性。

本文以筒子染色機(jī)的碟形封頭為研究對象，利用有限元分析軟件ANSYS Workbench對碟形封頭進(jìn)行參數(shù)化建模和有限元分析，利用多目標(biāo)驅(qū)動(dòng)優(yōu)化方法對封頭的結(jié)構(gòu)尺寸進(jìn)行優(yōu)化設(shè)計(jì)；運(yùn)用6西格瑪設(shè)計(jì)準(zhǔn)則，來得到最大等效應(yīng)力和總變形的概率分布；建立考慮應(yīng)力和變形相關(guān)性的可靠性計(jì)算模型，驗(yàn)證優(yōu)化方案的可行性。

1 筒子染色機(jī)染缸碟形封頭結(jié)構(gòu)簡介

碟形封頭主要由球冠、過渡圓弧環(huán)和短圓筒三部分組成。GF241XL-21型筒子染色機(jī)染缸的筒體內(nèi)徑Di=210 mm，據(jù)GB 150-2011《壓力容器》和JB/T 4746-2002《鋼制容器用封頭》；其碟形封頭的球冠半徑Ri=0.9 Di=189 mm，過渡半徑r= 0.17 Di=35.7 mm，圓整后取r=36 mm；碟形封頭計(jì)算厚度δ1=0.64 mm，考慮到鋼板負(fù)偏差及腐蝕裕量，參考鋼材標(biāo)準(zhǔn)規(guī)格取碟形封頭名義厚度［6］δ1n=3 mm；染缸外徑和內(nèi)徑的比值K=Di+2δ/Di= 210+6/210=1.028＜1.2，屬于薄壁容器［7］，薄壁筒體的計(jì)算厚度δ2=0.54 mm；同理取筒體的名義厚度δ2n=3 mm。當(dāng)筒體長度大于時(shí)，可以不考慮筒體過渡區(qū)存在的邊緣應(yīng)力的影響，取筒體長度為90 mm。根據(jù)ASMEⅧ-1［8］可知過渡區(qū)均需要1∶3的斜面過渡，筒子染色機(jī)的碟形封頭的結(jié)構(gòu)如圖1所示。

圖1 碟形封頭結(jié)構(gòu)

2 碟形封頭的有限元分析

2.1 碟形封頭參數(shù)化建模

碟形封頭屬于軸對稱結(jié)構(gòu)，采用2D建模和面體分析，可以節(jié)約計(jì)算時(shí)間。在對模型進(jìn)行網(wǎng)格劃分時(shí)，采用自動(dòng)劃分網(wǎng)格的方式，取單元尺寸為0.5 mm時(shí)可以得到劃分細(xì)密的網(wǎng)格，此時(shí)的單元個(gè)數(shù)為3 570個(gè)。碟形封頭的材料設(shè)置為奧氏體不銹鋼，密度ρ=7.93×103kg/m3，彈性模量E=200 GPa，泊松比μ=0.3。

施加合理邊界條件是保證仿真結(jié)果正確性的關(guān)鍵。在進(jìn)行仿真分析時(shí)，需結(jié)合碟形封頭的實(shí)際受力情況設(shè)置邊界條件。如圖2所示，（a）在筒體上端面約束Y方向的位移，（b）在球冠區(qū)對稱面上約束X方向的位移，（c）在封頭內(nèi)壁上施加均布內(nèi)壓載荷，p=0.5 MPa。從圖3仿真分析結(jié)果可知施加的載荷和約束與碟形封頭的實(shí)際工作情況是相符的。

圖2 碟形封頭邊界條件

2.2 有限元仿真分析

通過有限元仿真分析得到的碟形封頭應(yīng)力應(yīng)變云圖如圖3所示。圖3（a）顯示等效應(yīng)力最大值出現(xiàn)在球冠區(qū)與過渡區(qū)連接處，最大值為33.08 MPa；圖3（b）顯示總變形最大值出現(xiàn)在球冠區(qū)，最大值為0.04 mm。材料的許用值為［σ］tφ=115.2×0.85= 97.92 MPa，一般總變形不得超過3%，即許可變形量為［ξ］=0.09 mm，應(yīng)力和總變形符合壓力容器設(shè)計(jì)規(guī)范要求。

圖3 仿真分析結(jié)果

3 碟形封頭的可靠性優(yōu)化設(shè)計(jì)

3.1 碟形封頭相關(guān)參數(shù)的靈敏度分析

由于材料屬性、幾何尺寸、載荷等不確定性參數(shù)，對碟形封頭的應(yīng)力和變形有不同程度的影響，因此通過靈敏度分析來找出敏感性較大的參數(shù)。參考機(jī)械設(shè)計(jì)手冊，并假設(shè)該碟形封頭的相關(guān)設(shè)計(jì)參數(shù)均服從正態(tài)分布，可確定設(shè)計(jì)參數(shù)的平均值和標(biāo)準(zhǔn)差如表1所示。

表1 輸入?yún)?shù)取值

靈敏度的目標(biāo)函數(shù)是設(shè)計(jì)變量的導(dǎo)數(shù)，即隨機(jī)變量分布函數(shù)的變化引起失效概率變化的比率，可以用來評(píng)定輸入?yún)?shù)對輸出變量的影響程度。通過ANSYS Workbench仿真分析所得的靈敏度柱狀圖［10］如圖4所示。從圖4可以直觀地看出輸入?yún)?shù)對輸出參數(shù)的敏感程度：其中輸入?yún)?shù)球冠區(qū)過渡半徑（P7）、球冠半徑（P8）和封頭厚度（P11），內(nèi)壓載荷（P13）對輸出參數(shù)最大等效應(yīng)力和最大總變形的影響較大。

圖4 碟形封頭靈敏度分析結(jié)果

3.2 碟形封頭多目標(biāo)驅(qū)動(dòng)優(yōu)化設(shè)計(jì)結(jié)果

在ANSYSWorkbench中，根據(jù)上述結(jié)構(gòu)參數(shù)和靈敏度分析結(jié)果，可將靈敏度較高的關(guān)鍵結(jié)構(gòu)參數(shù)優(yōu)化的重要性設(shè)置為高，靈敏度較低的結(jié)構(gòu)參數(shù)設(shè)置為低。運(yùn)用實(shí)驗(yàn)數(shù)據(jù)法（DOE）進(jìn)行響應(yīng)面擬合，即在設(shè)置好的各結(jié)構(gòu)參數(shù)的取值范圍內(nèi)，利用數(shù)理統(tǒng)計(jì)方法和正交原理，從大量的設(shè)計(jì)點(diǎn)中挑選出具有代表性和典型性的設(shè)計(jì)點(diǎn)，并通過響應(yīng)面法擬合計(jì)算出45個(gè)更新的設(shè)計(jì)點(diǎn)。在多目標(biāo)驅(qū)動(dòng)優(yōu)化（GDO）時(shí)，取初始樣本數(shù)為1 000，以等效應(yīng)力最大值最小化和總變形最大值最小化為優(yōu)化目標(biāo)，程序基于此目標(biāo)從上述更新的45個(gè)設(shè)計(jì)點(diǎn)中自動(dòng)選取最合適的3個(gè)候選點(diǎn)，如表2所示。

表2 碟形封頭相關(guān)參數(shù)優(yōu)化候選點(diǎn)

對比表2中的三個(gè)候選點(diǎn)，綜合考慮等效應(yīng)力最大值和總變形最大值以及結(jié)構(gòu)尺寸參數(shù)取值，選擇B點(diǎn)作為新的設(shè)計(jì)點(diǎn)。與未優(yōu)化前的等效應(yīng)力最大值和總變形最大值對比，等效應(yīng)力最大值減小了32.3%，總變形最大值減小了45.0%。

3.3 碟形封頭可靠性計(jì)算模型

碟形封頭的可靠度可以定義為應(yīng)力S小于強(qiáng)度δ且同時(shí)總變形D小于允許變形量ξ的概率。假設(shè)應(yīng)力、強(qiáng)度、總變形量和允許變形量均服從正態(tài)分布，即

S～N（μS，），δ～N（μδ），D～N（μD），ξ～N（μξ），

令x=S-δ，y=D-ξ，則x和y分別服從均值為μx=μS-μδ和μy=μD-μξ，標(biāo)準(zhǔn)差為σx=和的正態(tài)分布，則隨機(jī)變量x和y的概率密度函數(shù)可表示為：

當(dāng)x＜0且y＜0時(shí)，碟形封頭正常工作，其可靠度R可表示為：

其中，g（x，y）表示變量x和y的聯(lián)合分布密度函數(shù)，由于變量x和y均服從正態(tài)分布，則其聯(lián)合密度函數(shù)為：

其中，ρ稱為相關(guān)系數(shù)，其定義為：

若變量x和y不相關(guān)，則ρ=0；若變量x和y呈線形關(guān)系，則｜ρ｜=1。

利用ANSYSWorkbench中的6西格瑪設(shè)計(jì)準(zhǔn)則，假設(shè)輸入的各結(jié)構(gòu)參數(shù)均為服從正態(tài)分布的的變量，并且優(yōu)化前后的結(jié)構(gòu)參數(shù)的方差均取表1中的方差值；根據(jù)碟形封頭實(shí)際承載情況，假設(shè)施加的內(nèi)載荷也是一個(gè)均值為0.5 MPa，方差為0.1 MPa的正態(tài)分布。運(yùn)用拉丁超立方體抽樣［11］方法隨機(jī)抽取800個(gè)樣本，得到優(yōu)化前后碟形封頭最大等效應(yīng)力和最大總變形的概率密度，結(jié)果如圖5所示。

圖5 碟形封頭6西格瑪分析結(jié)果

從圖5可見，最大等效應(yīng)力和最大總變形近似服從正態(tài)分布。優(yōu)化前最大等效應(yīng)力的均值估計(jì)值=39.3 MPa，標(biāo)準(zhǔn)差估計(jì)值=17.8 MPa，優(yōu)化前最大總變形量估計(jì)值=0.046 mm，標(biāo)準(zhǔn)差估計(jì)值=0.024 mm；優(yōu)化后最大等效應(yīng)力的均值估計(jì)值=27.4 MPa，標(biāo)準(zhǔn)差估計(jì)值=12.2 MPa，優(yōu)化后總變形量估計(jì)值=0.026 mm，標(biāo)準(zhǔn)差估計(jì)值=0.013 mm。取強(qiáng)度和許可總應(yīng)變的變異系數(shù)為0.2，即可認(rèn)為強(qiáng)度δ是服從均值為μδ =97.9 MPa，標(biāo)準(zhǔn)差為σδ=19.6 MPa的正態(tài)分布，許可變形量ξ是服從均值為μξ=0.09 mm，標(biāo)準(zhǔn)差為σξ=0.018 mm的正態(tài)分布。相關(guān)系數(shù)ρ可由模擬數(shù)據(jù)利用式（5）估計(jì)得出，其估計(jì)值?ρ=0.97。

那么優(yōu)化前隨機(jī)變量（x0，y0）～N（-58.3，-0.044，26.472，0.0302，0.97），優(yōu)化后隨機(jī)變量（x，y）～N（-70.5，-0.064，23.092，0.0222，0.97），通過式（3）可靠度計(jì)算公式可以得出優(yōu)化前碟形封頭的可靠度R0=0.922 8，優(yōu)化后可靠度R= 0.998 0?？梢?，優(yōu)化后可靠度有明顯的提高。

4 結(jié) 論

以筒子染色機(jī)染缸的碟形封頭為研究對象，利用ANSYS Workbench軟件建立了其參數(shù)化模型，并進(jìn)行了有限元分析；通過靈敏度分析，得到了影響碟形封頭最大等效應(yīng)力和最大總變形的關(guān)鍵參數(shù)；利用實(shí)驗(yàn)數(shù)據(jù)法（DOE）和多目標(biāo)驅(qū)動(dòng)優(yōu)化方法（GDO）從45個(gè)設(shè)計(jì)點(diǎn)中確定了最佳設(shè)計(jì)點(diǎn)；建立了考慮應(yīng)力和變形相關(guān)性的可靠性計(jì)算模型，運(yùn)用6西格瑪設(shè)計(jì)準(zhǔn)則和拉丁超立方體抽樣方法，得到了優(yōu)化前后最大等效應(yīng)力和總變形的概率分布，得出了優(yōu)化前后封頭的可靠度，表明優(yōu)化后可靠度有明顯提高，驗(yàn)證了本文提出方法的可行性。

［1］胡宜生.壓力容器封頭的可靠性設(shè)計(jì)［J］.安慶師范學(xué)院學(xué)報(bào)，2006（12）：25-26.

［2］王心明.工程壓力容器設(shè)計(jì)與計(jì)算［M］.北京：國防工業(yè)出版社，2011：24-25.

［3］農(nóng) 琪.壓力容器橢圓形封頭的可靠性設(shè)計(jì)［J］.裝備制造技術(shù)，2010（9）：53-54.

［4］周鵬飛，董金善，楊益清，等.基于ANSYS的異性封頭可靠性分析與優(yōu)化設(shè)計(jì)［J］.機(jī)械設(shè)計(jì)與制造，2012（1）：28-30.

［5］彭翠玲，艾華宇，劉青松，等.基于ANSYS的壓力容器可靠性分析［J］.核動(dòng)力工程，2009，30（1）：109-111.

［6］金秋紅.關(guān)于壓力容器設(shè)計(jì)時(shí)厚度的討論［J］.價(jià)值工程，2010（7）：116-117.

［7］李福寶，李 勤.壓力容器及過程設(shè)備設(shè)計(jì)［M］.北京：化學(xué)工業(yè)出版社，2009：36-38.

［8］Code P V.Rules for construction of pressure vessels［S］.

［9］Anghel C，Laz A I.Risk assessment for vessels affected by corrosion［J］.Chemical Engineering，2005，49（2）：103-118.

［10］王金龍.高壓力容器筒體與封頭連接區(qū)可靠性分析［J］.設(shè)計(jì)與研究，2007（8）：4-6.

［11］劉紀(jì)濤.基于拉丁超立方抽樣及響應(yīng)面的結(jié)構(gòu)模糊分析［J］.機(jī)械強(qiáng)度，2011，33（1）：73-76.

Reliability Optimization Design for Dished Head of Dye Vat of Cone Dyeing Machine

XU Rui-min，PAN Jun，CHEN Wen-hua，WANG Meng，XUCheng-bin
（Zhejiang Key Laboratory of Reliability Technology Research for Mechanical and Electrical Products，Zhejiang Sci-Tech University，310018，China）

Dished head of dye vat has extremely important influence on the safe reliability of the use of cone dyeing machine.With dished head of dye vat of cone dyeing machine as the research object，this paper conducts parametric modeling and finite element analysis on it with ANSYS Workbench；obtains key parameters influencing its equivalent stress and total deformation through sensitivity analysis and conducts optimization design for the structure size of the head with experimental data method and multi-objective driven optimization method；establishes a reliability calculation model considering stress and deformation correlation and obtains the probability distribution of maximum equivalent stress and maximum total deformation and the reliability of head before and after optimization with Six Sigma design criteria and Latin Hypercube Sampling method.The result shows that the reliability is greatly improved after optimization.

dye vat；dished head；ANSYS Workbench；sensitivity analysis；reliability optimization design

TH136

A

（責(zé)任編輯：張祖堯）

1673-3851（2014）01-0006-05

2013-04-26

國家自然科學(xué)基金（51075370，51275480），浙江省重點(diǎn)科技創(chuàng)新團(tuán)隊(duì)（2010R50005），浙江理工大學(xué)科研啟動(dòng)基金（1103824-Y）資助項(xiàng)目。

許銳敏（1989-），女，湖北咸寧人，碩士研究生，主要從事機(jī)電產(chǎn)品可靠性研究。