蘇懷維 張大有 毛瑞芝

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

1 引 言

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

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

2 裝置結構及工作原理

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

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

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

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

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

3 模型建立及求解

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

3.1 模型建立

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

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

3.2 網格劃分

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

圖4 活動閥網格Fig.4 Active valve mesh

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

3.3 求解設置

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

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

b)氣體為理想氣體；

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

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

3.3.1 流動屬性設置

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

3.3.2 材料屬性和邊界條件設置

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

3.3.3 動網格參數(shù)設置

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

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

圖6 動網格更新過程Fig.6 Dynamic mesh update process

3.3.4 求解器設置與求解

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

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

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

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

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

4.1 排氣口質量流量

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

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

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

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

圖8 恒壓150kPa,分析頻率分別為10Hz，100Hz，1000Hz，2000Hz時排氣口質量流量變化曲線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 結束語

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