蘇懷維 張大有 毛瑞芝

(1.北京理工大學，北京100081； 2.北京航天計量測試技術(shù)研究所，北京 100076)

1 引 言

為準確獲取壓力傳感器的動態(tài)特性，需要對壓力傳感器進行動態(tài)壓力校準。目前國內(nèi)廣泛采用的校準方法主要有階躍壓力校準、脈沖壓力校準和正弦壓力校準。正弦壓力校準方法是最理想的動態(tài)校準方法[1]，是檢驗壓力傳感器穩(wěn)定性、可靠性的最佳手段。在動態(tài)校準中，正弦壓力裝置最關(guān)鍵的問題在于如何產(chǎn)生一個穩(wěn)定、準確，并足以激發(fā)傳感器動態(tài)特性的正弦壓力信號?；诓煌懋a(chǎn)生正弦壓力信號的裝置主要有進口調(diào)制正弦壓力發(fā)生器、出口調(diào)制正弦壓力發(fā)生器以及進出口雙向調(diào)制正弦壓力發(fā)生器等。目前對正弦壓力發(fā)生器的結(jié)構(gòu)優(yōu)化與分析大都基于實際的測試試驗，成本較高，操作復(fù)雜，而且很難通過試驗工具來檢測壓力腔內(nèi)部的流場情況，因此探索一種能夠模擬正弦壓力發(fā)生器工作過程的數(shù)值分析方法具有重要意義。

本文基于所設(shè)計的一種具有濾波特性的新型正弦壓力發(fā)生器，以FLUENT仿真軟件為計算平臺，運用FLUENT動網(wǎng)格技術(shù)和UDF方法對正弦壓力的發(fā)生過程進行動態(tài)數(shù)值仿真計算，分析了正弦壓力腔內(nèi)部流場的流動情況及壓力變化規(guī)律，為正弦壓力發(fā)生器的優(yōu)化設(shè)計提供了一定的理論參考。

2 裝置結(jié)構(gòu)及工作原理

所設(shè)計的具有濾波特性的新型正弦壓力發(fā)生器由正弦壓力發(fā)生腔、壓力室主架、機械濾波器、開孔旋轉(zhuǎn)圓盤等組成。正弦壓力發(fā)生腔結(jié)構(gòu)示意圖如圖1所示，它包括主副結(jié)構(gòu)正弦壓力腔體，其兩側(cè)分別為副正弦壓力腔和主正弦壓力腔，在副正弦壓力腔和主正弦壓力腔之間相連通的通道內(nèi)放置有機械濾波器；副正弦壓力腔頂部和底部分別開了螺紋孔以便于安裝被測傳感器和標準傳感器；主正弦壓力腔頂部和底部分別開通孔形成進氣口和排氣口；主副結(jié)構(gòu)正弦壓力腔體的正弦壓力腔部分固定在壓力室主架中部開有的小孔內(nèi)；壓力室主架的頂部開有氣體入口，底部開有兩個排氣方孔。

圖1 正弦壓力發(fā)生腔結(jié)構(gòu)示意圖Fig.1 Sine pressure generating chamber structure diagram

該新型正弦壓力發(fā)生器屬于出口調(diào)制式，即通過周期性地改變氣體出口面積來實現(xiàn)發(fā)生腔內(nèi)壓力的正弦變化。恒壓穩(wěn)定的氣流通過進氣口進入壓力室，壓力室另一端開有兩個均勻長方孔，并與開有均勻圓孔的旋轉(zhuǎn)圓盤外側(cè)面相接觸。氣體通過壓力室排氣長方孔，經(jīng)旋轉(zhuǎn)盤上的圓孔排出。排氣面積由轉(zhuǎn)盤上圓孔與壓力室長方孔的相對位置決定。當轉(zhuǎn)盤勻速轉(zhuǎn)動時，排氣口開啟和關(guān)閉的面積按正弦規(guī)律變化，進而導(dǎo)致排出氣室的氣體流量也按正弦規(guī)律變化，這就引起了氣室內(nèi)壓力也按相應(yīng)的正弦規(guī)律變化。當改變旋轉(zhuǎn)圓盤的轉(zhuǎn)速時，氣室內(nèi)的正弦壓力信號的頻率也隨之改變。

主正弦壓力腔內(nèi)的正弦壓力信號通過引壓腔進入副正弦壓力腔，引壓腔內(nèi)放置有專用機械濾波器(如圖2所示)。通過機械濾波器的作用，即可抑制高次諧波，避免高頻正弦壓力波形嚴重失真，從而在副正弦壓力腔內(nèi)獲得所需要的穩(wěn)定正弦壓力信號。

圖2 新型正弦壓力發(fā)生器濾波結(jié)構(gòu)Fig.2 New sine pressure generator filter structure

3 模型建立及求解

為通過數(shù)值模擬的方法真實地反映出正弦壓力腔內(nèi)部壓力變化過程，采用FLUENT動網(wǎng)格技術(shù)和UDF方法對所建立的正弦壓力腔三維簡化模型進行求解計算。

3.1 模型建立

利用UG三維建模軟件根據(jù)新型正弦壓力發(fā)生腔的氣道尺寸建立三維流場計算區(qū)域。為減少運算量，節(jié)約運算時間，同時便于后期動網(wǎng)格處理，主正弦壓力腔的排氣口簡化為一個通孔，建立的氣道三維幾何模型如圖3所示。

圖3 正弦壓力腔氣道三維幾何模型Fig.3 Three-dimensional geometric model of sine pressure chamber airway

3.2 網(wǎng)格劃分

利用ANSYS Workbench環(huán)境中的網(wǎng)格劃分組件Mesh生成有限元網(wǎng)格。采用混合網(wǎng)格技術(shù)對三維幾何模型進行網(wǎng)格劃分，即對正弦壓力腔體部分采用非結(jié)構(gòu)化網(wǎng)格劃分，活動閥采用結(jié)構(gòu)化網(wǎng)格進行劃分(圖4)。為提升計算精度，運用局部控制選項Sizing控制排氣口管道和活動閥的網(wǎng)格尺寸為0.1mm。整個流場計算區(qū)域網(wǎng)格劃分結(jié)果如圖5所示。

圖4 活動閥網(wǎng)格Fig.4 Active valve mesh

圖5 氣道流場區(qū)域網(wǎng)格劃分圖Fig.5 Airway flow field area meshing map

3.3 求解設(shè)置

旋轉(zhuǎn)圓盤勻速旋轉(zhuǎn)時，主正弦壓力腔排氣口面積以正弦規(guī)律周期性變化，為便于對整個流場區(qū)域進行動網(wǎng)格分析以及運用UDF定義動網(wǎng)格的運動，將旋轉(zhuǎn)圓盤的勻速旋轉(zhuǎn)運動狀態(tài)簡化為一方形活動閥的往復(fù)勻速直線運動狀態(tài)。同時，F(xiàn)LUENT求解時做如下假設(shè)：

a)不考慮摩擦力和氣體泄露；

b)氣體為理想氣體；

c)滿足氣體流量的連續(xù)性條件；

d)忽略溫度因素的影響。

3.3.1 流動屬性設(shè)置

啟動FLUENT三維單精度求解器，F(xiàn)LUENT求解算法包括基于壓力(Pressure-Based)和基于密度(Density-Based)兩種，本文采用壓力求解器。由于存在移動的壁面，因此將時間項設(shè)置為Transient，即非定常流。

3.3.2 材料屬性和邊界條件設(shè)置

由于正弦壓力腔內(nèi)的流動屬于可壓縮流動，因此設(shè)置材料參數(shù)時應(yīng)選擇ideal-gas(理想氣體)。將發(fā)生腔進氣口和排氣口邊界條件分別設(shè)置為壓力進口和壓力出口，入口壓力根據(jù)所分析的壓力點設(shè)置為不同值，出口壓力默認為0kPa，保留其它默認設(shè)置。

3.3.3 動網(wǎng)格參數(shù)設(shè)置

正弦壓力發(fā)生器工作時，壓力腔排氣口周期性開啟與閉合，即排氣口固體壁面邊界是隨時間變化的，因此需采用動網(wǎng)格技術(shù)來實現(xiàn)排氣口狀態(tài)的周期性改變。由于方形活動閥采用結(jié)構(gòu)化網(wǎng)格進行劃分，因此在動網(wǎng)格計算過程中，網(wǎng)格的動態(tài)變化過程應(yīng)采用Smoothing(彈簧光滑模型)和Layering(動態(tài)分層模型)進行計算。在Dynamic Mesh Zones(動態(tài)區(qū)域)定義動網(wǎng)格的相關(guān)區(qū)域。定義活動閥的下壁面為Rigid Body(剛體運動)，四周壁面及上壁面為Deforming(變形)，動網(wǎng)格區(qū)域運動過程通過編寫用戶自定義函數(shù)UDF來控制[2]。

在預(yù)覽網(wǎng)格運動時，根據(jù)分析頻率設(shè)置相應(yīng)的時間步長和時間步數(shù)，即可動態(tài)顯示活動閥周期性開啟關(guān)閉過程中的網(wǎng)格更新狀態(tài)，從而動態(tài)模擬正弦壓力信號的產(chǎn)生過程，活動閥動網(wǎng)格區(qū)域周期性更新過程如圖6所示。

圖6 動網(wǎng)格更新過程Fig.6 Dynamic mesh update process

3.3.4 求解器設(shè)置與求解

基于壓力求解器的求解算法主要包括分離式算法和耦合式算法。對于分離式算法，F(xiàn)LUENT提供了SIMPLE、SIMPLEC、PISO和Fractional Step四種壓力速度耦合算法，本文采用SIMPLE算法。各流場迭代的亞松弛因子保持默認設(shè)置，因為這些默認設(shè)置都是根據(jù)各種算法的特點優(yōu)化得出的。

初始化流場，設(shè)置殘差監(jiān)視器。本文主要關(guān)心的是在活動閥運動過程中，正弦壓力發(fā)生腔內(nèi)壓力和排氣口質(zhì)量流量隨時間的變化，因此還需設(shè)置相應(yīng)的壓力和流量監(jiān)視器。采用壓力進口邊界條件對流場進行初始化，設(shè)置靜壓計算初始流場。流場初始化后，根據(jù)各分析正弦信號頻率點設(shè)置相應(yīng)的時間步長和迭代步數(shù)，單擊Calculate開始迭代計算。

4 數(shù)值計算結(jié)果分析

FLUENT動網(wǎng)格技術(shù)可以對正弦壓力發(fā)生腔排氣口面積周期性變化時其內(nèi)部流場的流體動力學特性進行分析，從而捕捉流場特性的瞬時變化，為正弦壓力發(fā)生器設(shè)計及其改進提供參考。由此根據(jù)分析需要，設(shè)置五個面平均壓力分析面A、B、C、D、F和質(zhì)量流量分析面E，圖7是各分析面所在位置。在FLUENT完成計算后，將生成的數(shù)據(jù)文件運用MATLAB即可繪出排氣口質(zhì)量流量變化曲線以及各分析面壓力變化曲線。

圖7 分析面定義Fig.7 Definition of Analysis surface

4.1 排氣口質(zhì)量流量

由正弦壓力發(fā)生器工作原理可知，排出氣室的氣體質(zhì)量的改變直接影響著正弦壓力腔內(nèi)壓力的變化[3]。利用動網(wǎng)格技術(shù)和UDF控制排氣口E面面積以正弦規(guī)律變化，即可通過監(jiān)測該分析面處質(zhì)量流量隨時間的變化關(guān)系來驗證該處氣體質(zhì)量流量是否按正弦規(guī)律變化，進而得出其內(nèi)部流場的壓力變化關(guān)系。

圖8是進氣口恒壓150kPa(絕壓)時分析頻率為10Hz、100Hz、1000Hz、2000Hz的排出氣體質(zhì)量流量變化曲線。從圖中可以看出，忽略迭代過程中殘差變化的影響，排出氣體的質(zhì)量流量近似是按給定頻率的正弦規(guī)律變化的。仿真過程中，當網(wǎng)格更新至排氣口面積為0時，流出氣體的質(zhì)量流量為0；當網(wǎng)格更新至排氣口面積最大時，流出氣體質(zhì)量流量達到最大值,這與實際正弦壓力發(fā)生器的工作機制較為吻合。

4.2 主副正弦壓力腔壓力波形

圖9是進氣口恒壓150kPa(絕壓)，排氣口開閉頻率為1500Hz時分析面A、B、C、D、F處的壓力波形。從圖中可以看出，由于排氣口處氣體回流現(xiàn)象和氣體在正弦壓力腔內(nèi)的反射擾動，使得主正弦壓力腔內(nèi)分析面A和B處的正弦壓力波形含有高次諧波成分而出現(xiàn)波形失真，傳感器感測面D和F處能夠產(chǎn)生正弦壓力信號。數(shù)值結(jié)果表明，主正弦壓力腔內(nèi)的失真信號經(jīng)過機械濾波器后，高次諧波被抑制在副正弦壓力腔外，因此副正弦壓力發(fā)生腔內(nèi)能夠產(chǎn)生穩(wěn)定的正弦壓力信號。

圖8 恒壓150kPa,分析頻率分別為10Hz，100Hz，1000Hz，2000Hz時排氣口質(zhì)量流量變化曲線Fig.8 Constant pressure 150kPa, the analysis of the mass flow rate of the exhaust port when the analysis frequency is 10Hz, 100Hz, 1000Hz, 2000Hz

圖9 排氣口開閉頻率1500Hz，進氣口恒壓150kPa(絕壓)時各分析面處波形 Fig.9 Waveform of the opening and closing frequency of the exhaust port is 1500 Hz, and the waveform of each analysis surface is at the constant pressure of the inlet port of 150 kPa (absolute pressure)

5 結(jié)束語

本文通過FLUENT動網(wǎng)格技術(shù)和UDF方法對正弦壓力發(fā)生器進行了動態(tài)模擬，得到了其主副正弦壓力腔內(nèi)流場的壓力以及壓力腔出口面流出氣體質(zhì)量流量隨時間的變化規(guī)律。數(shù)值仿真結(jié)果表明，本方法能夠較為真實地模擬出正弦壓力發(fā)生器的工作過程，驗證了該方法對正弦壓力發(fā)生器研究分析的可行性，為正弦壓力發(fā)生器的性能預(yù)測和結(jié)構(gòu)優(yōu)化提供了一定的參考依據(jù)。