林浩宇,蔣仲安,楊 斌,陳建武
(1. 北京科技大學(xué) 土木與資源工程學(xué)院,北京100083;2. 中國安全生產(chǎn)科學(xué)研究院,北京100029)
密集拋光打磨作業(yè)場所一直都是粉塵發(fā)生的重災(zāi)區(qū),由于工藝過程的特殊性,粉塵產(chǎn)生量大,每臺砂磨機平均產(chǎn)塵量為0.2 t/a,粉塵控制效果欠佳,現(xiàn)場粉塵濃度遠遠超過國家安全標(biāo)準(zhǔn)[1]。為了控制和治理打磨拋光粉塵,目前采取的主要手段有改進打磨工藝、使用除塵打磨設(shè)備,以及使用除塵排風(fēng)固定打磨臺等[2-4]。由于打磨產(chǎn)塵的必然性,治理粉塵從改進打磨工藝角度出發(fā)收效甚微;除塵打磨設(shè)備要兼顧除塵和打磨2種功用,其體積和風(fēng)量必受限制,除塵效率低;在固定打磨臺上進行拋光打磨作業(yè),可通過固定風(fēng)道將絕大部分粉塵排出,粉塵濃度降低顯著,但是對于配風(fēng)參數(shù)的設(shè)置卻缺少相關(guān)依據(jù)[5]。基于上述現(xiàn)狀,本文以數(shù)值模擬為主要研究手段,針對日益普及的組合式通風(fēng)打磨臺進行最佳除塵配風(fēng)參數(shù)的研究,以期為相關(guān)的工業(yè)實踐提供參考和借鑒。
粉塵運動和分布規(guī)律的研究,本質(zhì)上都屬于氣固兩相流的研究范疇,其數(shù)值模擬主要有2種方法:一種是將氣、固2種介質(zhì)看作1種混合流體;另一種是將2者分別用不同的方法研究,2者彼此獨立又相互作用[6-8]。本文將車間內(nèi)的風(fēng)流和打磨中產(chǎn)生的粉塵分別看作背景流體和顆粒相,分別運用歐拉法和拉格朗日法研究背景流體和顆粒相的運動。
空氣的流動選擇三維不可壓縮穩(wěn)態(tài)的Navier-Stokes方程作為控制方程,采用標(biāo)準(zhǔn)k-ε雙方程模型建立封閉控制方程組,只考慮動量的傳輸,忽略熱量傳遞[9-12],其控制方程如下。
1)N-S方程
(1)
2)標(biāo)準(zhǔn)k-ε雙方程
①k方程:
(2)
②ε方程:
(3)
對空氣中的粉塵顆粒進行量級分析,除重力、浮力、氣動阻力外,其他形式的作用力數(shù)量級都很小,可忽略不計,粉塵顆粒在空氣中的作用力平衡方程為:
(4)
式中:CD為阻力系數(shù);μ為空氣相的速度,m·s;μp為粉塵顆粒的運動速度,m·s;ρ為空氣密度,kg·m-3;ρp為粉塵顆粒的密度,kg·m-3;dp為粉塵顆粒的直徑,m。
為科學(xué)提出組合式通風(fēng)打磨臺最佳除塵配風(fēng)參數(shù),利用AutoCAD按照1∶1的比例建立三維幾何模型,模型導(dǎo)入ANSYS workbench劃分非結(jié)構(gòu)網(wǎng)格,如圖1所示。
根據(jù)實際應(yīng)用對模型進行合理簡化,打磨臺臺面為1.3 m × 0.6 m(長×寬)的矩形,距地面高0.8 m,并布有直徑φ=15 mm,圓心距d=25 mm的排風(fēng)圓孔;打磨臺后壁側(cè)為1.3 m × 0.9 m(長×高)的矩形,并在下方的2/3后壁面均勻布有100 mm×10 mm的條形排風(fēng)孔;打磨臺兩側(cè)為上底0.25 m,下底0.5 m,高0.9 m的直角梯形擋板;打磨臺臺面以下為粉塵收集裝置,簡化為0.5 m×0.3 m×0.25 m(長×寬×高)的規(guī)則長方體。打磨作業(yè)點假設(shè)位于臺面中央,粉塵源簡化為0.1 m×0.1 m的面源。計算域簡化為4 m×3 m×2.5 m(長×寬×高)的立方體,打磨臺位于計算域中央位置。
圖1 組合式通風(fēng)打磨臺三維模型及網(wǎng)格劃分Fig.1 Three-dimensional model and mesh of combined ventilation grinding platform
根據(jù)打磨作業(yè)實際情況及相關(guān)實測數(shù)據(jù),結(jié)合Fluent數(shù)值模擬方法和所建立的數(shù)學(xué)模型[13-17],確定固定打磨作業(yè)過程相關(guān)參數(shù)。求解參數(shù)及邊界條件設(shè)置,如表1所示;離散相模型及粉塵源參數(shù)設(shè)置,如表2所示。
本文是針對固定打磨的特定作業(yè),作業(yè)者位置固定,動作幅度變化小,操作人員擾動氣流對粉塵分布的影響可忽略。根據(jù)現(xiàn)場調(diào)研,打磨工件尺寸較小,打磨作業(yè)多針對工件棱邊進行,作業(yè)人員動作幅度變化小,并且單次打磨周期時間短、速度快,過程中產(chǎn)生的熱量少,所以忽略打磨位置動態(tài)變化及打磨產(chǎn)熱的影響。其中,作業(yè)過程中使用的手持打磨工具普遍采用砂輪直徑為125 mm的角磨機,標(biāo)準(zhǔn)轉(zhuǎn)速為12 000 rad·min-1,故砂輪邊緣線速度為12.5 m·s-1,并以打砂輪邊緣線速度作為金屬打磨粉塵產(chǎn)生的初速度。以固定打磨車間內(nèi)粉塵為樣本,過200目篩后,采用馬爾文激光粒度儀,干法測定粉塵粒度和分散度,如圖2所示。
表1 求解參數(shù)及邊界條件參數(shù)設(shè)置Table 1 Parameter setting of boundary conditions and solve
表2 離散相模型及粉塵源參數(shù)設(shè)置Table 2 Parameter setting of discrete phase and injection
圖2 固定打磨粉塵粒度分布Fig.2 Particle size distribution of fixed grinding dust
以角磨機砂輪傾斜45°打磨工況為例進行數(shù)值模擬,打磨臺除塵總風(fēng)量為2 000 m3·h-1,操作間平均風(fēng)速0.3 m·s-1,以臺面圓孔排風(fēng)量與后壁面條縫口排風(fēng)量之比為K值。
風(fēng)量分配比例分別為K=1,2,3,4,5時,粉塵源中心豎直向上粉塵質(zhì)量濃度沿程變化如圖3所示,圖4為粉塵最大濃度隨K值的變化。
圖3 距臺面不同高度粉塵濃度的變化Fig.3 The change of dust concentration different height from the work surface
圖4 最大粉塵濃度隨K值的變化Fig.4 The change of maximum dust concentration with the K value
由圖3和圖4可知:
1)在豎直方向,粉塵濃度隨高度增加先增大后減小,變化趨勢與K值無關(guān)。距臺面0.15~0.2 m高度時,粉塵濃度達到最大值;距臺面0.35 m時,粉塵濃度降低至粉塵最大濃度的10%,且粉塵濃度增加速率大于減小速率。
2)隨K值增大,粉塵最大濃度先增大后減小,在K=3時,粉塵最大濃度達到最大。粉塵最大濃度的上升速率(K<3)明顯大于下降速率(K>3),所以實際應(yīng)用中為了避開K=3的配風(fēng)比,減小K值效果顯著。
呼吸帶高度通常被認(rèn)為是距地面0.5~1.5 m左右,即人口鼻處相對于腳底高度,由于打磨操作特殊性,作業(yè)人員需要長時間適當(dāng)彎腰操作,所以呼吸帶最大高度需要適當(dāng)降低,以距臺面高度0.6 m(距地面1.4 m)作為人的呼吸帶高度。不同K值時,粉塵源正上方呼吸帶高度粉塵質(zhì)量濃度沿程變化如圖5所示。
圖5 不同K值時粉塵濃度分布Fig.5 Dust concentration change regularities under different the K values
由圖5可知:在呼吸帶高度,隨著K值不斷增大,粉塵濃度整體會呈現(xiàn)先減小后增加的變化規(guī)律,在K=1.67時,粉塵濃度達到了最小值,即此時的風(fēng)量分配為臺面圓孔排風(fēng)量1 250 m3·h-1,壁面條縫口排風(fēng)量750 m3·h-1。
組合式通風(fēng)打磨臺配風(fēng)比一定(K=1.67),總排風(fēng)量(Q)不同時,粉塵源中心豎直向上粉塵質(zhì)量濃度沿程變化如圖6所示,最大粉塵濃度隨總排風(fēng)量的變化規(guī)律如圖7所示。
圖6 距臺面不同高度粉塵濃度的變化Fig.6 The change of dust concentration different height from the work surface
圖7 最大粉塵濃度隨排風(fēng)量的變化Fig.7 The change of maximum dust concentration with exhaust air rate
由圖6和圖7可知:粉塵最大濃度出現(xiàn)在距臺面高度h=0.15 m附近,與總排風(fēng)量大小無關(guān),隨著總排風(fēng)量減小,粉塵主要聚集高度范圍會以距臺面0.15 m為中心逐漸擴大。在配風(fēng)比一定時(K=1.67),打磨臺的總排風(fēng)量越大,粉塵最大濃度越小,呈線性變化趨勢。
組合式通風(fēng)打磨臺不同總排風(fēng)量時,粉塵源正上方呼吸帶高度粉塵質(zhì)量濃度沿程變化如圖8所示。
圖8 不同總排風(fēng)量時呼吸帶粉塵濃度分布Fig.8 Dust concentration change regularities of breathing zone under different total exhaust air
由圖8可知:總排風(fēng)量從2 000 m3·h-1開始逐漸減小的過程中,呼吸帶粉塵濃度呈現(xiàn)出先減小后增大的變化趨勢,當(dāng)總排風(fēng)量Q=1 850 m3·h-1時,呼吸帶粉塵濃度達到最低。
模擬結(jié)果中,粉塵源中心豎直向上粉塵質(zhì)量濃度沿程變化顯示,粉塵最大濃度出現(xiàn)的位置與打磨臺通風(fēng)參數(shù)無關(guān),且低于作業(yè)人員的呼吸帶位置。對于控制粉塵最大濃度和呼吸帶粉塵濃度,考慮的控制對象不同,為達到最佳的粉塵控制效果,則應(yīng)采取不同的通風(fēng)參數(shù)。
無論是粉塵濃度最大的位置,還是作業(yè)人員呼吸帶位置,控制點的排風(fēng)風(fēng)速是打磨臺下排風(fēng)和側(cè)排風(fēng)在控制點形成的排風(fēng)風(fēng)流的矢量和。配風(fēng)比K值發(fā)生變化時(總排風(fēng)量不變),控制點的風(fēng)速大小發(fā)生變化,且風(fēng)流方向與垂直線的夾角也會發(fā)生變化,就造成了不同配風(fēng)比時的控制點排風(fēng)流的除塵效果不同。本文通過對不同K值條件下的通風(fēng)除塵工況進行數(shù)值模擬,找到最佳的K值,使得粉塵濃度滿足國家規(guī)定的安全限值,且濃度值最小。
打磨臺下排風(fēng)會加劇重力作用對打磨金屬粉塵的沉降,在加速粉塵沉降方面的作用明顯大于側(cè)排風(fēng)產(chǎn)生的降塵效果,側(cè)排風(fēng)主要是起到改變并控制粉塵沉降路徑的作用,使得粉塵沿著遠離作業(yè)人員的方向,即向打磨臺內(nèi)側(cè)沉降,將粉塵污染范圍控制在一定區(qū)域,從而大幅度降低了作業(yè)者呼吸帶位置(口鼻處)的粉塵濃度,能夠更好地保護作業(yè)人員的職業(yè)健康。
為了驗證FLUENT模擬結(jié)果的準(zhǔn)確性,選擇實際應(yīng)用中未進行參數(shù)優(yōu)化時(總排風(fēng)量2 000 m3/h,K=2),進行粉塵濃度的實際測量。根據(jù)GBZ/T 192.1—2007《工作場所空氣中粉塵測定》及相關(guān)文獻中粉塵采樣點布置方法,在固定打磨場所布置粉塵質(zhì)量濃度測點,取粉塵源上方呼吸帶高度水平方向和粉塵源豎直向上2條測點線,按照一定間距各布置12個測點。采用濾膜采樣器對粉塵質(zhì)量濃度進行測量,每個測點均進行不少于3次的數(shù)據(jù)測定,并取平均值,模擬結(jié)果與實測數(shù)據(jù)對比如圖9所示。
由圖9可知:數(shù)值模擬與現(xiàn)場實測數(shù)據(jù)結(jié)果在粉塵濃度曲線走向、波峰數(shù)目及其位置基本一致,數(shù)值誤差在10%以內(nèi),采用數(shù)值模擬的方法研究組合式通風(fēng)打磨臺最佳通風(fēng)參數(shù)是可行的,能夠反映一般規(guī)律。無論是水平還是豎直方向,現(xiàn)場測試的粉塵濃度都會略高于數(shù)值模擬的結(jié)果,這主要是由于現(xiàn)場有操作者站在臺面前進行打磨作業(yè),對空間風(fēng)流斜向打磨臺里側(cè)流動造成干擾,削弱了降塵效果導(dǎo)致的。
圖9 數(shù)值模擬與現(xiàn)場測塵結(jié)果對比Fig.9 Comparison of numerical simulation and field respirable dust measurement
1)打磨作業(yè)過程中,為控制空間內(nèi)粉塵最大濃度,應(yīng)采用配風(fēng)比K<3和較大的總排風(fēng)量。
2)以作業(yè)者呼吸帶高度粉塵質(zhì)量濃度為標(biāo)準(zhǔn),采用配風(fēng)比K=1.67,總排風(fēng)量Q=1 850 m3·h-1時,通風(fēng)除塵效果最佳。
3)本文的研究建立在該種尺寸的打磨臺的基礎(chǔ)上,揭示了1種從粉塵最大濃度和呼吸帶粉塵濃度與風(fēng)量分配和總排風(fēng)量的關(guān)系出發(fā),利用數(shù)值模擬的方法和手段,對組合式通風(fēng)打磨臺進行通風(fēng)除塵系統(tǒng)參數(shù)優(yōu)化的方法思路,對所有組合式通風(fēng)打磨臺的參數(shù)設(shè)計提供參考和借鑒。