摘要:爲(wèi)降低流(liu)體黏度(dù)對渦輪(lun)流量計(ji) 測量精(jing)度的影(yǐng)響，将渦(wō)輪流量(liang)計儀表(biao)系數線(xiàn)性度誤(wù)差最小(xiǎo)值作爲(wèi)目标函(han)數，在運(yùn)用計算(suan)流體力(li)學(CFD)仿真(zhen)的基礎(chǔ)上，先通(tōng)過Plackett－Burman設計(jì)篩選結(jie)構參數(shu)，并根據(ju)幾何結(jié)構對目(mu)标函數(shù)的影響(xiǎng)将其劃(hua)分爲兩(liang)個等級(jí)，即顯著(zhe)影響因(yīn)素和次(ci)顯著影(yǐng)響因素(su);再通過(guò)Box－Behnken設計及(jí)響應面(mian)法對顯(xian)著影響(xiǎng)因素進(jìn)行優化(hua)設計，分(fèn)析結構(gou)參數間(jiān)的交互(hù)作用，得(de)到參數(shu)的設計(ji)點;最後(hou)在響應(yīng)面分析(xi)基礎上(shàng)通過正(zheng)交試驗(yan)對次顯(xiǎn)著影響(xiǎng)因素進(jin)行優化(hua)設計，得(de)到最優(you)參數組(zu)合。對參(can)數組合(he)的渦輪(lún)流量計(jì)進行試(shì)驗研究(jiū)，試驗結(jié)果與CFD計(jì)算值吻(wěn)合，儀表(biǎo)系數線(xian)性度誤(wu)差由1.71%下(xia)降至1.59%，表(biǎo)明優化(huà)後的渦(wo)輪流量(liàng)計測量(liang)精度得(dé)到了顯(xiǎn)著提高(gāo)，基于響(xiang)應面法(fǎ)和正交(jiāo)試驗的(de)優化方(fāng)法可以(yi)用于渦(wo)輪流量(liang)計的結(jié)構設計(jì)。
引言
　　渦(wō)輪流量(liang)計具有(yǒu)精度高(gāo)、重複性(xìng)好、結構(gou)簡單、測(ce)量範圍(wei)廣、體積(ji)小、質量(liang)輕、壓力(lì)損失小(xiǎo)、維修方(fāng)便等優(you)點，但存(cun)在性能(néng)會随被(bèi)測流體(tǐ)黏度增(zēng)大而變(bian)差的問(wen)題。目前(qián)，國内的(de)渦輪流(liú)量計在(zài)出廠時(shi)，其性能(néng)一般都(dōu)是用水(shuǐ)或黏度(dù)比較低(di)的柴油(you)進行鑒(jiàn)定，但很(hěn)多使用(yong)者卻用(yòng)渦輪流(liú)量計來(lái)測量液(yè)壓油、潤(run)滑油等(děng)中黏度(du)甚至高(gāo)黏度液(ye)體的流(liú)量，導緻(zhi)出現很(hěn)大的測(cè)量誤差(chà)。因此，提(tí)高渦輪(lun)流量計(ji)在測量(liàng)黏性介(jie)質時的(de)精度具(ju)有非常(cháng)重要的(de)現實意(yì)義。
　　目前(qián)關于黏(nian)性介質(zhi)對渦輪(lún)流量計(ji)影響的(de)研究主(zhǔ)要集中(zhong)在分析(xī)流量計(jì)内部幾(jǐ)何結構(gòu)和流體(tǐ)介質對(dui)其性能(neng)的影響(xiang)以及儀(yí)表系數(shù)的修正(zheng)方法等(děng)方面，而(er)根據流(liú)體性能(néng)對流量(liang)計進行(hang)結構優(you)化的研(yán)究較少(shao)，在結構(gou)優化時(shi)考慮到(dào)内部幾(ji)何參數(shù)間交互(hù)作用的(de)則更少(shǎo)。由于渦(wo)輪流量(liàng)計幾何(he)參數較(jiao)多，作用(yong)的機理(li)各不相(xiang)同，各個(ge)參數之(zhi)間存在(zai)交互作(zuo)用，因此(cǐ)有必要(yào)研究各(gè)個參數(shu)間的相(xiang)互關系(xì)，确定最(zuì)優參數(shu)組合。以(yi)DN40渦輪流(liu)量計爲(wèi)例，從優(yōu)化幾何(he)結構出(chu)發，探究(jiu)幾何參(can)數對渦(wo)輪流量(liàng)計性能(neng)的影響(xiǎng)，分析顯(xian)著影響(xiǎng)因素之(zhi)間的交(jiao)互作用(yòng)，并在計(ji)算流體(ti)力學(CFD)仿(páng)真的基(ji)礎上通(tōng)過響應(ying)面法和(hé)正交試(shì)驗對結(jié)構進行(hang)優化設(she)計。
1模型(xing)與仿真(zhēn)
1.1模型的(de)建立
　　選(xuan)擇LWGY系列(liè)DN40渦輪流(liu)量計，其(qi)主要參(cān)數爲:葉(yè)輪葉片(pian)數N1=6，葉片(piàn)頂端半(ban)徑Ｒt=9.5mm，葉輪(lún)輪毂半(ban)徑Ｒo=10mm，葉輪(lun)輪毂長(zhǎng)度Lh=8mm，葉輪(lun)導程L=88.5mm，導(dǎo)流體葉(yè)片數N2=4，前(qian)導流體(tǐ)輪毂長(zhang)度H1=54mm，後導(dao)流體輪(lún)毂長度(dù)H2=38mm。按照上(shàng)述幾何(he)參數建(jian)立三維(wei)模型，如(rú)圖1所示(shi)。爲了使(shǐ)流體接(jie)近充分(fen)發展狀(zhuang)态從而(er)形成穩(wen)定的流(liú)速分布(bù)，在渦輪(lún)流量計(ji)前後分(fèn)别加裝(zhuang)10D和5D長直(zhi)管段［10］。

1.2網格劃(huà)分
　　将三(sān)維模型(xíng)導入網(wǎng)格劃分(fèn)軟件ICEM中(zhong)，考慮圖(tú)1渦輪流(liu)量計三(san)維模型(xing)Fig．13Dmodeloftheturbineflowmeter到流量(liang)計葉輪(lún)部分和(he)導流件(jiàn)部分結(jie)構複雜(za)，而且是(shì)仿真計(ji)算的關(guan)鍵部件(jiàn)，因此在(zai)ICEM中均采(cai)用非結(jié)構化的(de)四面體(ti)網格對(duì)葉輪流(liu)域和導(dǎo)流件流(liú)域進行(hang)劃分;而(er)前、後直(zhi)管段流(liu)域結構(gou)相對簡(jian)單但尺(chǐ)寸較大(dà)，采用四(si)面體網(wang)格劃分(fen)會使網(wang)格數量(liàng)大大增(zeng)加，爲了(le)減少仿(pang)真時間(jian)，采用結(jié)構化的(de)六面體(tǐ)網格對(dui)該流域(yu)進行劃(hua)分，劃分(fen)後的網(wǎng)格數爲(wei)1474621個，其Quali-ty最(zuì)小值爲(wèi)0.36。通過增(zeng)加整體(tǐ)網格數(shu)進行網(wǎng)格無關(guān)性檢驗(yan)，網格尺(chǐ)度符合(he)計算要(yào)求。
1.3邊界(jiè)條件定(dìng)義
邊界(jie)條件如(rú)下:
(1)仿真(zhēn)介質采(cai)用實際(ji)狀況下(xia)的原油(yóu)，其運動(dòng)黏度爲(wèi)2.64×10-5m2/s，密度爲(wèi)887kg/m3，流量範(fan)圍2～20m3/h;
(2)入口(kǒu)采用速(su)度入口(kou)，選取2m3/h、4m3/h、8m3/h、14m3/h、20m3/h這(zhè)5個體積(jī)流量下(xia)的入口(kǒu)速度;出(chu)口采用(yòng)壓力出(chu)口，設置(zhì)爲1個标(biao)準大氣(qì)壓;
(3)管壁(bi)，上、下導(dǎo)流體和(hé)葉輪表(biǎo)面均采(cǎi)用無滑(hua)移壁面(mian)邊界條(tiáo)件;
(4)渦輪(lun)流量計(jì)葉輪部(bù)分流域(yu)設置爲(wei)旋轉流(liú)域，前後(hòu)導流件(jiàn)部分設(shè)置爲固(gù)定流域(yù)，旋轉流(liú)域與固(gu)定流域(yù)之間采(cǎi)用交界(jie)面進行(háng)連接。
1.4湍(tuan)流模型(xing)的選擇(ze)
　　由于渦(wō)輪流量(liang)計葉輪(lun)在流體(ti)中處于(yú)高速旋(xuan)轉狀态(tai)，其表面(miàn)曲率變(bian)化非常(chang)大，而雷(léi)諾應力(li)模型(ＲSM)考(kǎo)慮到了(le)流體旋(xuán)轉或流(liu)線彎曲(qu)所帶來(lái)的應力(lì)張量的(de)急劇變(biàn)化，可以(yi)更好地(di)模拟渦(wō)輪流量(liang)計在複(fú)雜流場(chǎng)狀況下(xià)的運行(hang)規律，因(yin)此選用(yòng)ＲSM湍流模(mó)型［11］。
1.5仿真(zhen)儀表系(xi)數和線(xian)性度誤(wu)差的計(jì)算
　　儀表(biao)系數爲(wèi)渦輪感(gǎn)應放大(dà)器産生(shēng)的脈沖(chong)數與流(liu)過傳感(gan)器流體(tǐ)體積的(de)比值［12］。在(zài)計算仿(páng)真儀表(biǎo)系數之(zhi)前需要(yào)計算流(liu)量計葉(yè)輪在該(gāi)流量下(xia)的穩定(dìng)轉速θ。通(tōng)過監測(cè)發現，當(dang)葉輪驅(qu)動力矩(jǔ)與阻力(li)矩的差(chà)值小于(yú)10-8時，可認(rèn)爲葉輪(lun)所受力(li)矩達到(dao)平衡，則(ze)此時的(de)葉輪轉(zhuan)速即爲(wèi)穩定轉(zhuan)速。葉輪(lún)穩定轉(zhuan)速确定(dìng)後，根據(jù)葉片個(ge)數、入口(kou)流速與(yǔ)管道截(jie)面積可(ke)以得到(dao)此時的(de)渦輪流(liu)量計仿(pang)真儀表(biao)系數，其(qi)計算公(gōng)式爲

　　式(shi)中,K爲渦(wo)輪流量(liang)計仿真(zhen)儀表系(xi)數,L-1;N爲葉(yè)輪葉片(piàn)個數;?爲(wèi)葉輪穩(wěn)定轉速(sù),rad/s;V爲入口(kǒu)流速，m/s;A爲(wèi)前直管(guan)段入口(kou)截面積(ji),m2
　　儀表系(xì)數線性(xìng)度誤差(cha)可以反(fan)映渦輪(lún)流量計(jì)的測量(liàng)精度，儀(yí)表系數(shù)線性度(dù)誤差越(yue)小，則流(liu)量計的(de)測量精(jīng)度越高(gao),反之則(zé)測量精(jīng)度越低(di)。
　　通過式(shi)(1)計算出(chū)2m3/h、4m3/h、8m3/h、14m3/h、20m3/h這5個點(dian)的仿真(zhen)儀表系(xì)數後，便(bian)可以得(de)到渦輪(lún)流量計(ji)儀表系(xi)數線性(xìng)度誤差(chà)8,其計算(suan)公式爲(wei)

　　式中Kmin,i爲(wei)流量計(ji)在5個流(liu)量點處(chu)得到的(de)儀表系(xì)數最大(da)值;Kmin,i，爲流(liu)量計在(zai)5個流量(liàng)點處得(dé)到的儀(yí)表系數(shu)最小值(zhi)。
2Plackett-Burman設計
　　根(gen)據.Plackett-Burman(PB)試驗(yàn)設計，選(xuan)取8個試(shi)驗因素(su)(葉輪頂(dǐng)端半徑(jìng)、葉輪葉(yè)片數、葉(ye)輪輪毂(gū)半徑、葉(ye)輪輪毂(gu)長度、葉(yè)輪導程(cheng)、前導流(liú)件長度(dù)、後導流(liú)件長度(dù)、導流體(tǐ)葉片數(shu))和3個空(kong)白因素(sù),每個因(yin)素設高(gāo)、低兩個(ge)水平，以(yi)儀表系(xì)數線性(xìng)度誤差(cha)爲響應(ying)值,共計(jì)12個試驗(yan)，試驗設(she)計因素(su)及水平(píng)見表1。

　　對(duì)表1試驗(yan)中各因(yīn)素進行(háng)顯著性(xing)分析，分(fen)析結果(guo)如表2所(suo)示，模型(xing)顯著差(chà)異水平(ping)p=0.004,說明.回(huí)歸方程(chéng)關系顯(xiǎn)著;決定(ding)系數R2=0.9927,說(shuō)明回歸(gui)有效,試(shi)驗結果(guo)可靠。由(you)表2還可(kě)以看出(chu)8個因素(sù)均對流(liu)量計線(xiàn)性度誤(wù)差影響(xiang)顯著,其(qi)中葉輪(lún)葉片數(shu)N1、葉輪頂(ding)端半徑(jing)R1、葉片輪(lun)毂半徑(jing)R.。、葉輪輪(lun)毂長度(dù)Lh這4個爲(wei)顯著影(yǐng)響參數(shu),在後文(wén)中運用(yong)響應面(miàn)法進行(háng)優化;而(er)葉輪導(dao)程L、前導(dǎo)流件長(zhǎng)度H1、後導(dǎo)流件長(zhang)度H2、導流(liu)體葉片(pian)數N2這4個(ge)爲次顯(xian)著影響(xiang)參數,在(zài)後文中(zhōng)運用正(zheng)交試驗(yan)進行優(you)化。

3結構(gou)參數優(you)化.
3.1顯著(zhe)影響參(cān)數的響(xiang)應面法(fa)優化
　　采(cǎi)用Box-Behnken中心(xīn)組合設(shè)計方法(fa)，以葉輪(lún)頂端半(bàn)徑Rt、葉片(pian)數N1、葉片(pian)輪毂半(bàn)徑R.。、葉輪(lun)輪毂長(zhǎng)度Lh這4個(ge)顯著影(yǐng)響因素(sù)爲自變(biàn)量,儀表(biǎo)系數線(xiàn)性度誤(wù)差爲響(xiǎng)應值,其(qi)餘結構(gou)參數保(bǎo)持不變(bian),設計四(sì)因素三(san)水平29個(gè)試驗點(dian)的響應(ying)面優化(huà)試驗。因(yin)素與水(shui)平見表(biao)3,試驗設(she)計見表(biǎo)4。

對模型(xíng)進行方(fang)差分析(xi)得到的(de)響應面(miàn)分析結(jie)果如

　　表(biǎo)5所示,p<0.0001<0.05,表(biǎo)明該模(mo)型是顯(xiǎn)著的，具(jù)有統計(ji)學意義(yi)。由表5數(shù)據可得(de),自變量(liàng)R、N、R。、Lh均顯著(zhe)(p<0.05),按照對(dui)響應值(zhi)的影響(xiǎng)程度排(pai)序爲葉(yè)輪頂端(duan)半徑R1>葉(yè)輪葉片(piàn)數N1>葉輪(lun)輪毂長(zhang)度Lh>葉輪(lun)輪毂半(ban)徑R。失拟(nǐ)項P=0.056>0.05,此值(zhi)不顯著(zhe),說明在(zai)試驗範(fàn)圍内預(yu)測值和(he)實測值(zhi)的拟合(hé)度較高(gao)，能夠選(xuǎn)擇該回(hui)歸方程(chéng)對試驗(yan)結果進(jin)行相關(guan)分析，線(xiàn)性度誤(wù)差R的回(hui)歸方程(cheng)爲
R=17.22+3.37Rt+2.82N1-0.91R。-0.93Lh-0.95RtN1+0.75R1R。+0.2R1Lh+1.38N1R。-0.25N1Lh-1.70R。Lh+4.22R21+5.16N21+2.26R2。+2.97Lh2


　　從圖(tú)2(a)~(f)可形象(xiang)地看出(chū)影響儀(yí)表系數(shu)線性度(du)誤差的(de)幾何參(cān)數間交(jiāo)互作用(yong)。比較圖(tu)2各個分(fèn)圖可知(zhi),葉輪頂(dǐng)端半徑(jìng)R1對儀表(biǎo)系數線(xiàn)性度誤(wu)差的影(ying)響最爲(wèi)顯著，表(biǎo)現爲曲(qǔ)線最陡(dou)，其餘因(yīn)素影響(xiǎng)大小順(shùn)序爲葉(ye)輪葉片(piàn)數N1>葉輪(lun)輪毂長(zhǎng)度Lh>葉輪(lún)輪毂半(ban)徑R。這也(yě)與表5的(de)方差分(fen)析結果(guo)相吻合(he)。
　　使用DesignExpert軟(ruǎn)件在表(biǎo)3變量的(de)高低水(shuǐ)平範圍(wei)内尋優(yōu),以葉輪(lún)葉片數(shu)是整數(shù)爲前提(ti)，取其中(zhong)一個最(zuì)優組合(hé)進行CFD仿(páng)真計算(suan),并與顯(xian)著因素(su)的響應(ying)面回歸(guī)方程預(yù)測值進(jìn)行比較(jiao),比較結(jie)果如表(biao)6所示。可(ke)以看出(chu)，對于優(yōu)化後的(de)流量計(jì)模型,其(qi)儀表系(xì)數線性(xing)度誤差(cha)拟合公(gōng)式的預(yu)測值與(yǔ)CFD計算值(zhí)非常接(jie)近,誤差(chà)僅爲0.6%,說(shuō)明響應(ying)面法可(ke)以很好(hao)地用于(yú)渦輪流(liú)量計結(jie)構優化(huà)。
3.2次顯著(zhe)影響參(can)數的正(zhèng)交試驗(yan)設計
　　在(zài)對顯著(zhe)參數進(jìn)行響應(ying)面優化(hua)後,選擇(ze)葉輪導(dao)程L、前導(dao)流件長(zhang)度H、導流(liu)體葉片(pian)數N2、後導(dao)流件長(zhǎng)度H2這4個(ge)次顯著(zhe)影響因(yīn)素爲自(zì)變量,以(yi)流量計(ji)線性，度(dù)誤差爲(wei)響應值(zhi)進行正(zhèng)交試驗(yan)設計,根(gēn)據因素(sù)和水平(ping)數選擇(ze)正交表(biao)L9(34),一共9組(zu)仿真計(ji)算模型(xíng),因素與(yu)水平見(jian)表7。
　　正交(jiao)試驗結(jie)果與均(jun1)值如表(biao)8所示,因(yin)素L對應(yīng)的均值(zhi)2最小,表(biǎo)明L取第(di)二水平(píng)上的值(zhi)時線性(xìng)度誤差(chà)最小，同(tong)理可以(yǐ)得到H1、N2和(he)H2的取值(zhí)分别爲(wei):
　　第三水(shui)平、第二(er)水平和(hé)第二水(shuǐ)平,因此(cǐ)理論上(shang)的最優(yōu)水平組(zǔ)合爲L2(H)3(N2)2(H2)22。

　　圖(tu)3和圖4分(fen)别爲原(yuan)始流量(liang)計與優(you)化後流(liú)量計的(de)截面速(su)度分布(bu)雲圖。對(dui)比圖3和(he)圖4可以(yǐ)看出，優(yōu)化後的(de)流量計(ji)在後導(dǎo)流件.上(shàng)下部分(fèn)的流場(chang)速度分(fèn)布較原(yuán)始流量(liàng)計更加(jiā)均勻,說(shuo)明優化(hua)後流量(liang)計的葉(yè)輪結構(gòu)在流場(chǎng)中的旋(xuan)轉穩定(dìng)性更好(hǎo),從而使(shǐ)得測量(liang)精度得(de)到提高(gāo)。

4試驗驗(yan)證
　　将優(yōu)化前後(hòu)的渦輪(lún)流量計(ji)在流量(liàng)技術檢(jian)測試驗(yàn)台。上進(jin)行試驗(yàn)，試驗裝(zhuang)置如圖(tú)5所示。試(shi)驗介質(zhi)采用由(you)機油和(he)柴油按(àn)照一定(dìng)比例混(hun)合的密(mi)度爲887kg/m3、運(yùn)動黏度(du)爲2.64×10-5m2/s的混(hùn)合液,采(cǎi)用靜态(tài)容積法(fǎ)原理,利(lì)用泵爲(wei)流體提(ti)供動力(lì),流體經(jīng)過流量(liang)控制閥(fá)和被測(ce)渦輪流(liu)量計後(hòu)直接流(liú)回容積(jī)池中。分(fèn)别選取(qǔ)2m3/h、4m3/h、8m3/h、14m3/h、16m3/h、20m3/h這6個體(ti)積流量(liang)點,通過(guo)計算機(ji)控制台(tai)采集每(mei)個流量(liàng)點下試(shi)驗流量(liang)計産生(shēng)的脈沖(chòng)個數N,從(cong)而得到(dào)渦輪流(liú)量計在(zài)6個流量(liàng)點下的(de)儀表系(xi)數K。試驗(yan)中每個(ge)流量點(diǎn)分别進(jìn)行3次重(zhòng)複性試(shi)驗，試驗(yàn)誤差均(jun1)小于0.025%。

　　将(jiang)原始流(liu)量計儀(yí)表系數(shù)與優化(huà)後的流(liu)量計儀(yí).表系數(shu)進行對(dui)比，結果(guo)如圖6所(suo)示。流量(liang)計儀表(biǎo)系數CFD仿(pang)真值與(yǔ)試驗值(zhi)吻合,證(zhèng)明了CFD仿(pang)真的準(zhun)确性。對(dui)比原始(shǐ)渦輪流(liú)量計和(he)優化後(hou)渦輪流(liu)量計儀(yi)表系數(shù)試驗值(zhí)曲線得(de)出,優化(hua)後的流(liu)量計在(zai)不同流(liú)量處的(de)儀表系(xì)數變化(hua)情況較(jiào)原來更(geng)爲平穩(wěn)。經計算(suan)，儀表系(xì)數線性(xing)度誤差(chà)由原來(lái)的1.71%下降(jiàng)到了1.59%，顯(xiǎn)著提高(gāo)了渦輪(lún)流量計(ji)的測量(liang)精度。

5結(jie)論
(1)Plackett-Burman設計(ji)分析表(biǎo)明，對渦(wō)輪流量(liang)計測量(liang)精度影(ying)響顯著(zhe)的參數(shu)爲葉輪(lún)頂端半(ban)徑、葉輪(lún)葉片數(shù)、葉輪輪(lun)毂半徑(jing)和葉輪(lún)輪毂長(zhang)度,影響(xiǎng)次顯著(zhe)的參數(shu)有葉輪(lún)導程、前(qian)導流件(jian)長度、後(hou)導流件(jian)長度和(he)導流體(ti)葉片數(shu)。
(2)運用Box-Behnken設(she)計方法(fa)對篩選(xuǎn)出來的(de)顯著影(yǐng)響參數(shu)進行試(shi)驗設計(ji)，建立了(le)渦輪流(liú)量計線(xiàn)性度誤(wù)差的多(duo)元回歸(guī)模型,并(bing)檢驗了(le)預測模(mo)型的拟(nǐ)合度。結(jie)果表明(míng),回歸模(mó)型對實(shí)際情況(kuang)拟合較(jiào)好，能夠(gòu)運用響(xiang)應面法(fa)對渦輪(lún)流量計(jì)結構參(cān)數進行(háng)優化。
(3)在(zai)響應面(mian)法優化(huà)的基礎(chǔ)上,對篩(shai)選出來(lái)的次顯(xian)著影響(xiang)參數進(jin)行正交(jiao)試驗設(shè)計,得到(dao)了最優(you)結構組(zu)合。試驗(yan)驗證結(jié)果表明(míng)優化後(hou)的渦輪(lún)流量計(ji)測量精(jing)度得到(dao)了顯著(zhe)提高。

以(yǐ)上内容(róng)源于網(wang)絡，如有(yǒu)侵權聯(lián)系即删(shān)除!