吳利瑞 孟濤 韓磊 胡強

1 同濟大學(xué)機械與能源工程學(xué)院

2 同濟大學(xué)建筑設(shè)計研究院（集團）有限公司

實驗動物房空調(diào)系統(tǒng)的任務(wù)除了要控制室內(nèi)的溫度、濕度和風(fēng)速外，還要控制室內(nèi)污染物濃度。對室內(nèi)污染物控制的主要手段之一是維持各個動物房間合理的壓力梯度，保持氣流由高潔凈區(qū)到低潔凈區(qū)的流向，避免室外或低潔凈級別區(qū)的空氣回流造成對工作區(qū)的污染。如何判斷實驗動物房各區(qū)域間空氣的流動值得研究。本文通過一種有效而簡便的辦法研究單走廊SPF實驗動物房各個區(qū)域間的壓力波動情況，從而為研究實驗動物房各區(qū)域間空氣的流動提供幫助。

1 動物房壓力模型的建立

SPF實驗動物房按照潔凈走廊的布局方式分為單走廊實驗動物房和雙走廊實驗動物房兩種。實驗動物房一般包括動物飼養(yǎng)室、動物檢疫室、緩沖室和潔凈走廊。本文建立了一個簡化的單走廊實驗動物房潔凈區(qū)域模型，如圖1所示。包括三個動物飼養(yǎng)室，一條潔凈走廊，一個動物檢疫室，兩個緩沖室。每個房間都有一個通往潔凈走廊的自動門，緩沖室1和緩沖室2分別各有一個通往室外的自動門，各個門窗的編號如圖所示。

圖1 單走廊動物房壓力模型

空氣流經(jīng)兩房間之間的縫隙時流量和壓差之間的關(guān)系可以由下式來表示[1]：式中：Q為通過縫隙的氣體體積流量，m3/h；P為縫隙兩側(cè)的壓差，Pa；C 為流量系數(shù)，m3/(h·Pa)；n 為流量指數(shù)，范圍是0.5～1，當(dāng)氣體流動為充分發(fā)展的紊流時取極限值0.5，當(dāng)流動為層流時取1。工程中較大尺寸的縫隙往往起著主導(dǎo)作用，因此為方便使用，一般取流量指數(shù)n=0.5，可靠性達到90%以上，因此本文為了方便討論，取 n=0.5[2]。

除了圖中標(biāo)明位置的縫隙之外，每個房間實際上還有其它縫隙，如有些動物飼養(yǎng)室會有外窗，但對于同一類型的縫隙，縫隙兩側(cè)的壓力都相同，可以歸結(jié)為一處。由于潔凈室的壓力變化相對大氣壓而言很小，因此由于壓力變化引起的空氣密度的改變忽略不計，即認為空氣是不可壓縮的。本文的研究目的是分析各個實驗動物房間在影響因素變化的條件下壓力的波動情況，因此可以忽略各潔凈房間內(nèi)部壓力不均勻性，即認為每個潔凈區(qū)只有一個靜壓值。

由式（1），對于每個房間而言其進入的空氣量等于流出的空氣量，可以得到如下的方程組：

2 壓力模型的參數(shù)設(shè)定與求解

2.1 壓力模型參數(shù)的設(shè)定

對于上述方程組，已知流量系數(shù)C1～C8和壓差風(fēng)量Q1～Q7可以求解出各個房間的靜壓值；同樣的道理，已知流量系數(shù)C1～C8和各個動物房間要求的靜壓值就可以求解出各個動物房的壓差風(fēng)量。本文的目的是得出在動物房流量系數(shù)變化和靜壓風(fēng)量變化的情況下，各個動物房間壓力的變化情況。所以需要給定模型中的流量系數(shù) C1～C8和壓差風(fēng)量 Q1～Q7。

對于SPF實驗動物房，要求不同等級的動物房間壓差保持在10Pa。對于單走廊實驗動物房而言，靜壓值可以設(shè)置如下[3]：動物飼養(yǎng)室保持30Pa靜壓值，潔凈走廊保持20Pa靜壓值，與潔凈走廊相連接的緩沖室保持10Pa靜壓值，檢疫室保持5Pa的靜壓值。假定所有的門窗流量滲透系數(shù)（如表1），就可以解出所有的壓差風(fēng)量（如表2）。

表1 門窗流量系數(shù)（m3/(h·Pa1/2)）

表2 各個動物房的靜壓風(fēng)量（m3/h）

有了所有的流量系數(shù)和壓差風(fēng)量，只要改變其中任意一個條件，就可以通過求解方程組，得出各個房間靜壓值與各個房間壓差風(fēng)量和流量系數(shù)間的變化關(guān)系。

2.2 壓力模型的求解

針對上述的非線性方程組，本文利用MATLAB軟件，用牛頓下山法編寫程序計算。牛頓下山法的算法如下[4～5]：

1）給定初值 x0，控制精度 ε，令 k=0；

2）如果函數(shù)值滿足 |f(xk)|＜ε，返回近似解 xk，迭代終止；

3）計算 dk=-f(xk)/f'(xk)，δ=1；

4）如果 |f(xk+δdk)|＜|f(xk)|，令 xk+1=xk+δdk，否則 δ折半，重新計算第三步；

5）令迭代步k增加1，轉(zhuǎn)到第一步。

3 求解結(jié)果分析

3.1 流量系數(shù)變化對潔凈區(qū)域的影響

流量系數(shù)C2在潔凈走廊和動物房1之間，從圖2中可以看出，隨著C2的變大，只有動物房1的靜壓值不斷變小，其他各個區(qū)域的靜壓值不會改變。從方程組也可以看出C2改變只是引起了靜壓值P2的變化，動物房1和潔凈走廊的靜壓風(fēng)量都沒有改變。與此類似，流量系數(shù)C3～C5的變化也只會改變各自房間所對應(yīng)房間的靜壓值，不會對整個區(qū)域的壓力造成影響。

圖2 流量系數(shù)C2變化對各個房間靜壓的影響

流量系數(shù)C6在潔凈走廊和緩沖室2之間，從圖3中可以看出，隨著C6的變大，整個潔凈區(qū)域的靜壓值都發(fā)生相應(yīng)的改變。這是由于，C6的變化引起了潔凈走廊靜壓值的變化，所有與潔凈走廊相連的區(qū)域靜壓值都受到影響。而且從所有曲線變化趨勢可以看出，隨著流量系數(shù)的增大到一定程度，靜壓值變化趨于平緩；各個區(qū)域所對應(yīng)的流量系數(shù)越小，如檢疫室和動物房3，靜壓值變化越快（曲線斜率越大），即區(qū)域壓力抵抗外界干擾能力越弱。

圖3 流量系數(shù)C6變化對各個房間靜壓的影響

流量系數(shù)C8在緩沖室2和室外之間，從圖4中可以看出，與C6變化情況類似，隨著C8的變大，整個潔凈區(qū)域的靜壓值都發(fā)生相應(yīng)的改變。

圖4 流量系數(shù)C8變化對各個房間靜壓得影響

3.2 滲透風(fēng)量變化對潔凈區(qū)域的影響

從圖5、圖6可以看出，滲透風(fēng)量的變化都會引起整個區(qū)域靜壓值的變化，各個區(qū)域靜壓值變化的快慢程度與該區(qū)域所對應(yīng)的滲透流量正相關(guān)。但總體上可以看出潔凈走廊滲透凈風(fēng)量變化對整個區(qū)域壓力都有較大的影響，而其它某個動物房的滲透凈風(fēng)量變化對其它潔凈區(qū)域的靜壓值影響不大。

圖5 走廊壓差造成的滲透風(fēng)量變化對各個房間靜壓的影響

圖6 動物房1壓差造成的滲透風(fēng)量變化對各個房間靜壓的影響

4 結(jié)論

通過對壓力模型的求解和分析結(jié)果可知，對于單走廊實驗動物房：

1）任意房間滲透風(fēng)量的變化都會對整個潔凈區(qū)域的壓力分布產(chǎn)生影響，隨著某個房間滲透風(fēng)量的加大，各個房間靜壓值都會有不同程度的變大。

2）與潔凈走廊相連的房間所對應(yīng)的流量系數(shù)改變只對該房間的靜壓值有影響，不會改變其它房間的靜壓值，對應(yīng)房間靜壓值隨著流量系數(shù)的增大而減小。

3）其它區(qū)域的流量系數(shù)改變會引起整個區(qū)域靜壓值的變化。

[1]Brain Wiseman.Room Pressure for Critical Environments[J].ASHRAE Journal,2003,45(2):34-39

[2]杜世元.潔凈室壓力控制研究[D].上海:同濟大學(xué),2007

[3]全國實驗動物標(biāo)準(zhǔn)化技術(shù)委員會.實驗動物環(huán)境及設(shè)施(GB14925-2010)[S].2010

[4]吳鵬.MATLAB高效編程技巧與應(yīng)用[M].北京:北京航空航天大學(xué)出版社,2010

[5]同濟大學(xué)計算數(shù)學(xué)教研室.現(xiàn)代數(shù)值數(shù)學(xué)和計算[M].上海:同濟大學(xué)出版社,2003