相對于現(xiàn)有的流量檢測技術(shù)[1-5], 如:電磁流量計、渦街流量計、渦輪流量計、轉(zhuǎn)子流量計、差壓流量計、質(zhì)量流量計等, 超聲波流量計具有與被檢測介質(zhì)的非接觸性、通用性廣、適應(yīng)性強以及安裝與維護(hù)簡便等特點, 因此超聲波流量計被廣泛地應(yīng)用于現(xiàn)代工業(yè)的各個領(lǐng)域[6-9], 尤其是石油、天然氣、水利、化工等行業(yè).其發(fā)展經(jīng)歷了漫長的過程, 1931年德國科學(xué)家Ruttgen[1-2]首次提出了利用時差法測量流量的概念, 并進(jìn)行了可行性的理論預(yù)測, 但在當(dāng)時的儀器精度下難以實現(xiàn).以此為基礎(chǔ), 此后出現(xiàn)了大量的理論與實驗驗證, 但直到1955年美國科學(xué)家使用“鳴環(huán)法”[1-2], 即:利用聲音在液體中傳播過程中的多次反射所產(chǎn)生的時間差計算液體的流量, 才將超聲波流量計從實驗室應(yīng)用到工業(yè)領(lǐng)域, 然而由于此檢測電子系統(tǒng)的過于復(fù)雜, 使得該流量計并未被廣泛推廣應(yīng)用.20世紀(jì)70年代隨著電子技術(shù)的快速發(fā)展, 超聲波流量計的穩(wěn)定性、可靠性以及測量精度得到了很大提高, 同時也促進(jìn)了流量測量原理的迅速發(fā)展, 如:時差法、多普勒法、波束偏移法、相關(guān)法、空間濾波法及噪聲法等[4,7].目前常用的算法主要是時差法和多普勒法, 測量度也相對較高[8-10].時差法測量原理適用于純凈的或含有顆粒直徑小于1mm的流體材料檢測, 而多普勒法檢測原理適合于含有較大顆粒的流體檢測.目前, 超聲波流量計的研究主要集中超聲流量計的設(shè)計以及算法的修正上[10-12], 而對管道內(nèi)的流場分布特性上尚有不足、針對流場測試中的不足之處, 本文利用有限元方法, 構(gòu)建圓形管道中流場的理論分析模型, 并數(shù)值分析管道中的空氣流場在風(fēng)機驅(qū)動下管道內(nèi)流速的分布特性;實驗上利用時差法測量流量原理, 實驗測定相同管道中空氣在風(fēng)機驅(qū)動下流場的平均流速, 并與理論分析結(jié)果相印證.

使用2對中心頻率均為40kHz的超聲波傳感器作為探測器, 用于檢測口徑為18cm的有機玻璃導(dǎo)管中的流體速度.2對超聲波傳感器的安裝方式如圖1所示, 將4個超聲波傳感器兩兩配對構(gòu)成兩組測試探測器, 其中每組中各含1個接收和發(fā)射傳感器, 將2組傳感器分別安裝在被測導(dǎo)管的兩側(cè), 且每組中的收/發(fā)超聲探測器和另外一組中的發(fā)/收超聲探測器分別一一相對, 并保持聲程一致.這樣當(dāng)管道中有速度為v0的流體流動時, 經(jīng)過相應(yīng)聲程所需要的時間為

其中, D為管徑, θ為超聲發(fā)射方向與流體流速方向的夾角, k為由于儀器誤差所加的修正值, ±分別代表超聲傳播的聲速矢量方向與流體速度矢量方向相同或相反.

測試獲得的兩路超聲信號間的時間差為

由于v0遠(yuǎn)小于c, 化簡 (2) 式可得

實驗中采用直徑為18cm的有機玻璃管作為測量管道, 在管道兩側(cè)成64°夾角放置2組超聲波探頭, 實驗中采用函數(shù)信號發(fā)生器驅(qū)動超聲波傳感器, 分別將采集到的信號輸入前置放大電路, 將經(jīng)前置放大的信號輸入數(shù)字示波器顯示并采集記錄到計算機, 由計算機進(jìn)行***后的信號處理, 以獲取管道中流體的流速.

在如圖1所示的實驗裝置的一端使用小型風(fēng)機向測試管內(nèi)送風(fēng), 使用DF601熱式智能風(fēng)速變送器標(biāo)定風(fēng)機的輸出風(fēng)速為2.860m/s.使用示波器監(jiān)控2組超聲探頭的輸出信號, 其結(jié)果如圖2所示, 其中虛線所示為超聲傳播的聲速矢量方向沿管軸的分量與流體速度矢量方向相同的超聲信號, 實線為逆向信號.實線所示信號滯后于虛線所示信號, 其滯后的時間量反映了超聲波信號被管道中流場所調(diào)制而導(dǎo)致的超聲波信號在相同聲程條件下傳輸時間的變化特性.

由互相關(guān)理論可知, 互相關(guān)函數(shù)反映了2個相互間存在時延信號的相似程度, 而該相似度不僅可以用于判斷信號的存在與否, 還可以推斷信號到達(dá)的時間和距離, 即可以依據(jù)互相關(guān)函數(shù)***大值時的延遲量, 計算出兩接收信號之間的時間差.依據(jù)此原理, 計算逆流和順流時所獲取超聲波信號的互相關(guān)特性, 其結(jié)果如圖3所示.依據(jù)兩通路信號的互相關(guān)系數(shù)得到逆流方向所獲取的聲信號比順流方向所獲取的聲信號滯后約26個采樣點, 由采樣周期可得逆流信號比順流信號共滯后5.2μs, 將該結(jié)果代入式 (3) , 可得管道中流體的流速v0=0.814m/s.與實際測量得風(fēng)機出口的流速相比, 此結(jié)果遠(yuǎn)小于風(fēng)機出口的風(fēng)速.

以圖1所示的測量管道為模型, 利用有限元分析方法構(gòu)建在標(biāo)準(zhǔn)大氣壓下測量管道中的流體模型, 在管道中心處加載口徑為50 mm、流速為2.860 m/s的流場, 模擬實驗中所使用的風(fēng)機.在此條件下數(shù)值分析管道內(nèi)流場, 其模擬結(jié)果如圖4所示.取管道中與實際測試位置相同處的管道橫切面, 該橫切面的流速分布如圖5所示, 管道中流速的等速分布圖呈中心對稱, 且流場中的速度分布呈中心高、邊沿低的分布特性, 其中心***大流速約為2.781m/s, 然后逐漸降低, 在邊緣處約為0.080m/s.

根據(jù)實驗公式[13], 管內(nèi)流速分布可描述為

將圖5中的模擬結(jié)果代入式 (4) 中, 使用***小二乘法擬合管道中的數(shù)值流場速度分布特性, 則可得到模擬流場中的平均流速為0.859m/s.此結(jié)果與實驗中所得到的流場速度0.814 m/s接近, 偏差為5.2%.由于在實驗分析中, 采用的聲速按照340m/s的標(biāo)準(zhǔn)值進(jìn)行數(shù)值分析, 未對實際測試環(huán)境下的聲速進(jìn)行標(biāo)定[14];再加上忽略了有機玻璃管道對流動空氣的黏滯性影響.綜合上述誤差因素, 數(shù)值模擬證明, 實驗中所獲取的管道中流場的流速為管道中整個流場速度的平均值, 且與管道中的流場均值速度基本相符.

從實驗原理以及測試的過程中的信號處理上看, 如在以下幾方面進(jìn)行適當(dāng)改進(jìn)將提高實驗結(jié)果的準(zhǔn)確性.首先, 如在實驗進(jìn)行前使用標(biāo)準(zhǔn)流速進(jìn)行預(yù)先的實驗標(biāo)定將有助于實驗結(jié)果的準(zhǔn)確性;其次, 采用高質(zhì)量高頻超聲傳感器也將有效地提高測量結(jié)果, 但由于高頻超聲在流體中的衰減較快, 因而在選擇傳感器時應(yīng)兼顧測量量程的需求, 在測程和準(zhǔn)確性之間尋求***佳的匹配;此外, 適當(dāng)?shù)靥岣卟蓸铀俾室部商岣邷y量結(jié)果的準(zhǔn)確性, 但其對高頻采樣的硬件電路要求較高, 易導(dǎo)致成本的非線性增長.

采用基于有限元的數(shù)值分析方法分析了管道中空氣流場在外界風(fēng)機作用下的空氣流速分布特性, 并使用時差法建立了針對圓柱形管道的超聲波流場測試計量實驗系統(tǒng), 利用該測試系統(tǒng)可獲取管道內(nèi)流場的平均流速, 且該流速與基于有限元的數(shù)值模擬結(jié)果基本一致.

致謝:在本課題的完成過程中, 南京大學(xué)聲學(xué)所為本項目提供了相關(guān)的實驗儀器、場地, 南京大學(xué)聲學(xué)所張仲寧老師為實驗裝置設(shè)計提出了建設(shè)性的改進(jìn)措施.