• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    帶式輸送機(jī)轉(zhuǎn)載溜槽的抑塵設(shè)計(jì)

    2020-11-25 01:55:14宋冠霆宋偉剛
    礦山機(jī)械 2020年11期
    關(guān)鍵詞:方法模型設(shè)計(jì)

    宋冠霆,宋偉剛

    1中煤科工集團(tuán)沈陽(yáng)設(shè)計(jì)研究院有限公司 遼寧沈陽(yáng) 110015

    2東北大學(xué)機(jī)械工程與自動(dòng)化學(xué)院 遼寧沈陽(yáng) 110819

    帶式輸送機(jī)轉(zhuǎn)載溜槽系統(tǒng)是散料輸送過(guò)程中粉塵的主要來(lái)源,由于環(huán)保和人員健康的要求,有關(guān)誘導(dǎo)空氣與粉塵生成機(jī)理和抑塵設(shè)計(jì)的研究成為熱點(diǎn)[1-3]。隨著對(duì)轉(zhuǎn)載溜槽抑塵機(jī)理的深入研究,抑塵方法逐漸從除塵器的動(dòng)力抑塵向無(wú)動(dòng)力抑塵方式轉(zhuǎn)變;對(duì)逃逸的粉塵處理從灑水到噴霧、干霧方式的轉(zhuǎn)變。英國(guó)紐卡斯?fàn)柎髮W(xué)建立了轉(zhuǎn)載溜槽粉塵測(cè)試實(shí)驗(yàn)臺(tái),研究抑塵的計(jì)算方法和合理抑塵結(jié)構(gòu)參數(shù)分析[3-5];澳大利亞臥龍崗大學(xué)(University of Wollongong)的 Arnold、Wypych 和國(guó)內(nèi)的夏博實(shí)等人對(duì)粉塵下落的誘導(dǎo)空氣量進(jìn)行理論研究,提出了誘導(dǎo)空氣的計(jì)算模型[1,6,7];劉澤勤等人[8-9]研究了自由下落物料誘導(dǎo)空氣的計(jì)算方法并進(jìn)行試驗(yàn);李小川等人也對(duì)誘導(dǎo)空氣量及影響因素進(jìn)行了研究[10-11]。

    本文是帶式輸送機(jī)轉(zhuǎn)載溜槽設(shè)計(jì)計(jì)算方法研究系列之一,將對(duì)轉(zhuǎn)載溜槽誘導(dǎo)空氣的計(jì)算模型與方法進(jìn)行評(píng)述,進(jìn)而對(duì)抑塵設(shè)計(jì)方法進(jìn)行分析,對(duì)無(wú)動(dòng)力抑塵采用的引風(fēng)管、沉降室的粉塵沉降過(guò)程和空氣刮刀結(jié)構(gòu)進(jìn)行仿真分析,最后給出無(wú)動(dòng)力抑塵導(dǎo)料槽結(jié)構(gòu)參數(shù)的選擇方法。

    1 轉(zhuǎn)載過(guò)程的空氣流及誘導(dǎo)空氣量

    轉(zhuǎn)載點(diǎn)處的溜槽外形及其他部分幾何形狀的良好設(shè)計(jì)對(duì)減少粉塵的產(chǎn)生有越來(lái)越重要的作用,要很好地理解溜槽的幾何形狀與粉塵逃逸的關(guān)系,必須理解其機(jī)理,包括溜槽內(nèi)部的粉塵生成。

    粉塵的生成隨空氣速度增大而增加,隨物料尺寸、內(nèi)聚力的增大而減小。在每個(gè)轉(zhuǎn)載點(diǎn),物料流的方向發(fā)生變化,導(dǎo)致不同顆粒尺寸(其質(zhì)量也不同)的物料軌跡分離。無(wú)動(dòng)力抑塵轉(zhuǎn)載溜槽系統(tǒng)如圖1 所示,其設(shè)計(jì)的關(guān)鍵是將轉(zhuǎn)載溜槽全部密封,對(duì)粉塵進(jìn)行控制。一種觀(guān)點(diǎn)認(rèn)為,當(dāng)在轉(zhuǎn)載溜槽系統(tǒng)上設(shè)置風(fēng)機(jī)引流除塵時(shí),粉塵量將增加 4 倍。

    圖1 無(wú)動(dòng)力抑塵轉(zhuǎn)載溜槽系統(tǒng)Fig.1 Unpowered dust suppression transfer chute system

    誘導(dǎo)空氣是由于運(yùn)動(dòng)物料曳引作用而隨物料流動(dòng)的空氣,帶入設(shè)備或防塵密閉罩中的空氣量為誘導(dǎo)空氣量。在大量的顆粒自由下落時(shí),周?chē)目諝怆S著物料的加速運(yùn)動(dòng)而擴(kuò)散,在氣固兩相間曳力、壓力梯度力等的相互作用和影響下,會(huì)被卷吸到料流的空隙中去(見(jiàn)圖2)[3],這也是粉體或細(xì)顆粒物料在轉(zhuǎn)載過(guò)程充氣后密度明顯降低的原因。由于誘導(dǎo)空氣流和料流存在的速度差,使吸附在料塊表面的粉塵在摩擦力作用下剝離飛揚(yáng),并向四周飄逸,形成污染。圍繞著料流中心會(huì)形成一個(gè)粉塵邊界層,邊界層的半徑隨著落料高度的增加而增大,誘導(dǎo)空氣量的大小和很多因素有關(guān),例如顆粒的密度和粒徑分布、落料高度、顆粒質(zhì)量流量、物料流進(jìn)口的大小等。一些試驗(yàn)是在實(shí)驗(yàn)室進(jìn)行的,采用的測(cè)試設(shè)備為小型設(shè)備,且試驗(yàn)用物料為極細(xì)的氧化鋁粉、沙石等,因而很難代表實(shí)際物料轉(zhuǎn)載的誘導(dǎo)空氣。下面闡述幾種基于典型的誘導(dǎo)空氣模型的誘導(dǎo)空氣量的估算方法。

    圖2 誘導(dǎo)氣流產(chǎn)塵原理示意Fig.2 Sketch of dust generation principle of induced air flow

    2 誘導(dǎo)空氣計(jì)算模型與誘導(dǎo)空氣量估算方法

    2.1 Hemeon 模型方法

    劉澤勤[8]給出 Hemeon 模型的推導(dǎo)過(guò)程,假設(shè)在靜止的空氣中,誘導(dǎo)氣流具有一定的末速度所需要的功等于顆粒自由下落過(guò)程中空氣阻力所做的功。將誘導(dǎo)氣流看成單顆粒物料在靜止氣流中自由下落產(chǎn)生的氣流的總和,誘導(dǎo)空氣量

    式中:Qind為誘導(dǎo)空氣量,m3/s;g為重力加速度,m/s2;Qm為料流質(zhì)量流量,kg/s;h為物料下落高度,m;A1為料流截面積,m2;dp為顆粒名義直徑,m;ρp為顆粒密度,kg/m3。

    此計(jì)算方法僅適用于單顆粒。假設(shè)物料開(kāi)始下落時(shí)顆粒速度為零,且無(wú)法計(jì)入實(shí)際料流內(nèi)部顆粒與顆粒、顆粒與空氣之間的相互作用,此計(jì)算結(jié)果比實(shí)際值偏大,不適合轉(zhuǎn)載溜槽的誘導(dǎo)空氣量預(yù)測(cè)。

    2.2 Tooker 模型方法

    Tooker[12]指出 Hemeon 功率方程中的總壓力,系統(tǒng)的總壓力是速度壓力和系統(tǒng)阻力壓力之和,沒(méi)有直接給出誘導(dǎo)空氣的計(jì)算式,但給出了輸送物料為煤和巖石的帶式輸送機(jī)轉(zhuǎn)載溜槽中產(chǎn)生的誘導(dǎo)空氣量。按照 Tooker 的計(jì)算原理,誘導(dǎo)空氣

    式中:Effe為誘導(dǎo)效率,考慮物料間的碰撞與摩擦,一般為0.2~0.3,對(duì)于大下落高度,取0.5~0.6;A2為空氣進(jìn)口截面積,m2;Ce為空氣進(jìn)入系數(shù),一般取0.6。

    Tooker 建議,輸送煤和巖石時(shí)顆粒的名義粒度為1.0 和12.5 mm,該方法考慮了空氣流入的阻力,當(dāng)入口較大時(shí),1 m2數(shù)量級(jí)上誘導(dǎo)空氣量過(guò)大;當(dāng)入口面積減小2 個(gè)數(shù)量級(jí)時(shí),誘導(dǎo)空氣量急劇減少,這也說(shuō)明了減小入口面積可以有效地減少粉塵。

    2.3 CEMA 計(jì)算方法

    CEMA[13]認(rèn)為:粉塵逃逸(排放)的主要原因是溜槽內(nèi)部的正壓,使得空氣中漂浮的粒子逃逸到外面。因?yàn)榱锊鄄⒉煌耆芊?,物料的進(jìn)入使壓力升高是直接原因。在轉(zhuǎn)載點(diǎn)有3 個(gè)運(yùn)動(dòng)空氣來(lái)源,分別為置換空氣Qdis、誘導(dǎo)空氣Qind和生成空氣Qgen。

    (1)置換空氣 當(dāng)一定體積的物料進(jìn)入溜槽時(shí)同樣體積的空氣被推出,流出空氣的量等于進(jìn)入溜槽的物料的體積,即

    例如,當(dāng)轉(zhuǎn)載溜槽的輸送量(輸送煤)為3 500 t/h,堆積密度為900 kg/m3,其置換空氣量將達(dá)到1.08 m3/s=3 888.88 m3/h。置換空氣只有在系統(tǒng)啟動(dòng)前溜槽中沒(méi)有物料時(shí)才有作用,而在系統(tǒng)穩(wěn)定運(yùn)行過(guò)程中溜槽中的物料量基本保持不變,因而在實(shí)際系統(tǒng)設(shè)計(jì)時(shí)不需要考慮置換空氣問(wèn)題。

    (2)誘導(dǎo)空氣 誘導(dǎo)空氣是物料離開(kāi)頭部滾筒后沿物料軌跡運(yùn)動(dòng)而帶動(dòng)的空氣。輸送物料時(shí)含有一定量的夾帶空氣,在下落過(guò)程中顆粒分離,且它們每一個(gè)都拖動(dòng)空氣,被拖動(dòng)空氣的總量除了與物料下落的速度有關(guān)外,還與顆粒間距和顆粒尺寸有關(guān)。當(dāng)物料轉(zhuǎn)載且大部分的壓縮空氣釋放時(shí),引起轉(zhuǎn)載區(qū)空氣壓力的增加,向各方向推動(dòng)粉狀物料。

    誘導(dǎo)空氣總量的近似計(jì)算方法為

    式中:k為轉(zhuǎn)換系數(shù),k=0.078。

    (3)生成空氣 是向輸送機(jī)轉(zhuǎn)載區(qū)給料的設(shè)備運(yùn)行而產(chǎn)生的空氣,比如粉碎機(jī)、錘式粉碎機(jī)、木材削片機(jī)等,應(yīng)根據(jù)具體設(shè)備確定這部分空氣的總量,或者實(shí)際測(cè)得空氣流量。在帶式輸送機(jī)之間轉(zhuǎn)載時(shí)沒(méi)有生成空氣。

    CEMA 的氣流計(jì)算方法來(lái)源于馬丁工程公司的安全維護(hù)指南,從量綱關(guān)系看,該方法是在 Hemeon計(jì)算式基礎(chǔ)上的修正。但將料流截面積替換為進(jìn)口截面積時(shí),此計(jì)算方法在入口不采取密封措施的情況下過(guò)高地估計(jì)了誘導(dǎo)空氣量;而在入口采取密封措施后,誘導(dǎo)空氣量隨截面積減小而減少,若能在進(jìn)料區(qū)完全密封,可將誘導(dǎo)空氣降至 0。

    2.4 Cooper 和 Arnold 模型方法

    Cooper 和 Arnold[1]對(duì)不同粒度物料下落過(guò)程中的行為提出 2 種分析模型。

    2.4.1 模型 1

    第1 種模型假設(shè)從漏斗中下落的大顆粒物料全部以重力加速度的方式運(yùn)動(dòng),與下落物料的動(dòng)量相比,忽略了流動(dòng)的空氣動(dòng)量。誘導(dǎo)空氣為物料截面積乘以下落過(guò)程的速度差,誘導(dǎo)空氣量

    式中:up0為物料自由落體初速度,m/s。

    該方法估計(jì)的誘導(dǎo)空氣量雖然相對(duì)較小,但沒(méi)有計(jì)入系統(tǒng)的截面積和入口阻力的影響,因此所估計(jì)的誘導(dǎo)空氣量仍相對(duì)較大。

    2.4.2 模型 2(羽流模型)

    第2 種模型認(rèn)為:從漏斗流出的微粒羽流運(yùn)動(dòng)與自由下落過(guò)程中粉塵排放情況最為接近,且這種流動(dòng)被認(rèn)為是以高斯速度分布的混相羽流運(yùn)動(dòng)。根據(jù)羽流運(yùn)動(dòng)的方程式,推導(dǎo)出了小顆粒粉末從高度h下落時(shí)的誘導(dǎo)空氣量

    式中:C為與誘導(dǎo)作用α有關(guān)的系數(shù),α=ve/vplume,對(duì)于氧化鋁粉,C=0.1;ve為羽流邊沿誘導(dǎo)氣流速度,m/s;vplume為誘導(dǎo)氣流豎直方向的速度,m/s。

    值得注意的是,Wypych[6]所提到的誘導(dǎo)空氣量與下落高度的 5/3 次冪、與物料的質(zhì)量流量成正比,此說(shuō)法來(lái)源于Cooper 和Arnold[1],而非 Wypych首次提出。

    Cooper 和 Arnold 模型是研究粉體物料下落過(guò)程中產(chǎn)生誘導(dǎo)空氣的生成機(jī)理,不能直接應(yīng)用于轉(zhuǎn)載溜槽的誘導(dǎo)空氣預(yù)測(cè)。

    2.5 《火力發(fā)電廠(chǎng)及變電所供暖通風(fēng)空調(diào)設(shè)計(jì)手冊(cè)》方法

    《火力發(fā)電廠(chǎng)及變電所供暖通風(fēng)空調(diào)設(shè)計(jì)手冊(cè)》[14]給出了輸送機(jī)轉(zhuǎn)載過(guò)程中誘導(dǎo)空氣量的計(jì)算方法,誘導(dǎo)空氣量

    式中:ξ為雙密封導(dǎo)料槽所采用的系數(shù),取值范圍為 1.0~1.5,單密封時(shí)取 3;A為落料管橫截面積,m2;φ為與料流末速度有關(guān)的系數(shù);ρm為物料的堆積密度,t/m3;vend為料流末速度,對(duì)于直接下落情況,假設(shè)下落前物料速度為 0,則下落后的速度vend=(m/s)。

    該方法還給出了料流末速度的計(jì)算方法和針對(duì)不同結(jié)構(gòu)參數(shù)的誘導(dǎo)空氣數(shù)據(jù)表,《DL/T 51872—2004火力發(fā)電廠(chǎng)運(yùn)煤設(shè)計(jì)技術(shù)規(guī)程——2 煤塵防治》[15]引用并擴(kuò)展了此誘導(dǎo)空氣數(shù)據(jù)表,并規(guī)定:“當(dāng)運(yùn)煤系統(tǒng)落煤管上不加任何鎖氣裝置時(shí),可參照附錄 A 數(shù)據(jù)選用;當(dāng)運(yùn)煤系統(tǒng)各部位密封較好,落煤管上加裝鎖氣器、導(dǎo)料槽內(nèi)安裝二級(jí)橡皮擋簾時(shí),其除塵抽風(fēng)量可參照附錄 A 數(shù)據(jù)的 1/3 選用”。同時(shí)在條文說(shuō)明中介紹:“運(yùn)煤系統(tǒng)各部位密封較好,在落煤管末端加裝緩沖鎖氣器,導(dǎo)料槽內(nèi)加裝雙層滿(mǎn)擋式中間擋簾,除塵抽風(fēng)量是本規(guī)程附錄 A 所查風(fēng)量的 22%~26%。考慮輸送帶寬度、落差及帶速等因素的影響,盡量使理論計(jì)算與實(shí)際相吻合,使設(shè)計(jì)者能在不同條件下直接準(zhǔn)確地確定風(fēng)量,因此,在第 2 條中推薦采用附錄A 中數(shù)據(jù)的 1/3?!?/p>

    可見(jiàn),DL/T 51872—2004[15]所規(guī)定的抽風(fēng)量是為了保守起見(jiàn)而推薦的較大值,這也是轉(zhuǎn)載溜槽抑塵設(shè)計(jì)的一個(gè)誤區(qū)。

    2.6 誘導(dǎo)空氣的算例分析

    夏博實(shí)等人[7]對(duì)現(xiàn)場(chǎng)溜槽的出口空氣流速進(jìn)行了測(cè)試,物料質(zhì)量流Qm=3 500 t/h=972.22 kg/s,導(dǎo)料槽出口面積為 0.221 6 m2,無(wú)護(hù)罩情況下落高度為7.36 m,氣體流速為 9.8 m/s,空氣流量為 2.17 m3/s=7 818.048 m3/h;有護(hù)罩下落高度為 35.1 m,氣體流速為 1.8 m/s,空氣流量為 0.399 m3/s=1 436 m3/h,帶速為 3.8 m/s,物料顆粒密度為 1 500 kg/m3,堆積密度為 900 kg/m3,下落高差為 5 m。

    采用前述計(jì)算方法的計(jì)算結(jié)果如表1 所列,計(jì)算中所用其他參數(shù)為:Ce=0.6,A2=0.5 m2,Effe=0.5,dp=0.01 m,A=4A1,vend==9.9 m/s,A1==0.170 6 m2。

    表1 不同計(jì)算方法的誘導(dǎo)空氣量Tab.1 Induced air volume obtained with various calculation methods m3/s

    如前所述,Hemeon 模型以及 Cooper 和 Arnold模型不適用于轉(zhuǎn)載溜槽誘導(dǎo)空氣的估算。從表1 計(jì)算結(jié)果看:Tooker 模型和 CEMA 計(jì)算方法的計(jì)算結(jié)果比“設(shè)計(jì)手冊(cè)”的結(jié)果大 1 個(gè)數(shù)量級(jí)。但在 Tooker模型中,當(dāng)A2=0.000 5 m2時(shí),Qind=0.09 m3/s=325.27 m/h。減小進(jìn)料口的面積(通過(guò)設(shè)置擋簾)可有效地減少誘導(dǎo)空氣,方法是卸料輸送機(jī)頭部滾筒處的良好密封。對(duì)于設(shè)計(jì)手冊(cè)方法,當(dāng)下落高度為 35 m時(shí),空氣流量為 1.37 m3/s,也遠(yuǎn)大于實(shí)測(cè)值??梢?jiàn),現(xiàn)存誘導(dǎo)空氣的估計(jì)方法均遠(yuǎn)大于實(shí)際值,這在工程上偏于保守,也是工程設(shè)計(jì)的誤區(qū):即為了達(dá)到轉(zhuǎn)載溜槽的除塵效果,通過(guò)增加除塵器風(fēng)機(jī)的抽風(fēng)量,盡量使導(dǎo)料槽內(nèi)形成負(fù)壓,從而在除塵器內(nèi)除塵,但結(jié)果卻發(fā)現(xiàn),隨著抽風(fēng)量增大,導(dǎo)料槽的出口粉塵超標(biāo)問(wèn)題卻更加嚴(yán)重,隨著空氣的流動(dòng),粉塵也流動(dòng)。而設(shè)計(jì)的目標(biāo)是減緩空氣流動(dòng),使其所含的能量不足以承受粉塵的負(fù)荷,因此必須中斷空氣流動(dòng)。降低誘導(dǎo)空氣量的其他措施包括:在溜槽形狀上應(yīng)盡量減小物料流的截面積,以使輸送物料更加聚集;在溜槽的幾何形狀上也需要防止落下的加速度過(guò)大,或者最好完全限制自由落體的量;另一個(gè)重要且實(shí)用的措施是開(kāi)放區(qū)A2的減小。因此采用無(wú)動(dòng)力抑塵方法優(yōu)于風(fēng)機(jī)的動(dòng)力除塵方式,且工程實(shí)踐已經(jīng)證實(shí),良好的無(wú)動(dòng)力抑塵方式可使粉塵量降低 98%。

    3 轉(zhuǎn)載溜槽抑塵的計(jì)算機(jī)仿真方法

    粉塵的產(chǎn)生不可避免,但可以抑制。轉(zhuǎn)載溜槽的粉塵抑制問(wèn)題可以通過(guò)間接或直接的仿真方法去評(píng)價(jià)粉塵的抑制效果。當(dāng)前有關(guān)轉(zhuǎn)載溜槽的設(shè)計(jì)計(jì)算方法包括連續(xù)模型、DEM 和 CFD 以及他們之間的耦合方法。

    3.1 間接仿真方法

    間接仿真方法是根據(jù)工程經(jīng)驗(yàn)給出具有良好導(dǎo)流的轉(zhuǎn)載溜槽的結(jié)構(gòu)設(shè)計(jì),通過(guò)連續(xù)模型和 DEM 仿真,驗(yàn)證轉(zhuǎn)載過(guò)程中的料流集中度高、料流速度變化小,來(lái)說(shuō)明系統(tǒng)的設(shè)計(jì)具有較低程度的抑塵效果,這種方法可用于不同設(shè)計(jì)方案的對(duì)比。

    (1)離散元方法(DEM)是溜槽設(shè)計(jì)的輔助工具,以解決堵塞、溢出和磨損的問(wèn)題,也可以間接說(shuō)明所設(shè)計(jì)的系統(tǒng)具有較低的粉塵,但是由于沒(méi)有氣流,因此無(wú)法對(duì)轉(zhuǎn)載溜槽的空氣流動(dòng)及粉塵進(jìn)行直接預(yù)測(cè)。

    (2)單相流 CFD 仿真方法 是針對(duì)轉(zhuǎn)載溜槽系統(tǒng)建立幾何模型,單獨(dú)對(duì)空氣流動(dòng)仿真,其基本假設(shè)為:較大的出口空氣流速會(huì)導(dǎo)致更多的粉塵產(chǎn)生,而較大的空氣質(zhì)量流意味著會(huì)產(chǎn)生更多的粉塵。例如通過(guò) CFD 仿真氣體流場(chǎng)的壓力分布,當(dāng)導(dǎo)料槽內(nèi)的氣壓低于外部靜壓,說(shuō)明粉塵不會(huì)從導(dǎo)料槽溢出。

    (3)連續(xù)-CFD 方法 是采用 Roberts 理論計(jì)算料流,再在其基礎(chǔ)上進(jìn)行單相 CFD 仿真。單相 CFD仿真和連續(xù)-CFD 方法都無(wú)法準(zhǔn)確地對(duì)粉塵控制性能進(jìn)行預(yù)測(cè),只能作為定性分析。

    (4)間接 DEM-CFD 耦合法 Katterfeld 等人[16]提出間接 DEM-CFD 耦合法,通過(guò) 1 次 DEM 分析后導(dǎo)出散狀材料的幾何輪廓和孔隙率分布,2 次 CFD 分析分別對(duì)多孔介質(zhì)進(jìn)行 CFD 求解,沿散狀物料流輸出空氣質(zhì)量流量,再對(duì)周邊區(qū)域進(jìn)行 CFD 求解,最后得出結(jié)果。間接耦合的 DEM-CFD 可以部分預(yù)測(cè)溜槽性能,并具有減少計(jì)算時(shí)間的優(yōu)點(diǎn),因?yàn)樵摲椒ǖ膽?yīng)用僅限于簡(jiǎn)單的幾何形狀,其準(zhǔn)確性在很大程度上取決于對(duì)邊界條件的假設(shè)。

    3.2 直接仿真方法

    直接方法是在仿真過(guò)程中含有粉塵,可以得到粉塵的分布與粉塵的溢出量。

    (1)多相 CFD(TEM)方法 Donohue[16]采用 3相 CFD 模型研究粉塵抑制問(wèn)題,氣體相、固體相用1.0 mm 的顆粒表示糧食,用 0.1 mm 的顆粒來(lái)表示粉塵,通過(guò)仿真可以看到,粉塵在糧食顆粒的下面,驗(yàn)證了所改進(jìn)的曲線(xiàn)溜槽具有抑塵效果。

    Chen(陳曉玲)[3]對(duì) 6 種轉(zhuǎn)載溜槽結(jié)構(gòu)進(jìn)行試驗(yàn)研究,并通過(guò) 2 相流 CFD 仿真,采用歐拉-歐拉方法,驗(yàn)證了 CFD 兩相流算法的有效性。其固相顆粒尺寸為1.0 mm,認(rèn)為導(dǎo)料槽上設(shè)引風(fēng)管對(duì)粉塵控制效果不明顯,增加導(dǎo)料槽處沉降室的高度對(duì)抑塵意義不大。沉降室的高度模型試驗(yàn)是在 90 mm 高度(按比例 5 對(duì)應(yīng)實(shí)際的沉降室高度 450 mm)和增加 1 倍(實(shí)際沉降室高度 900 mm)的情況下進(jìn)行的,如圖3 所示。這說(shuō)明馬丁公司沉降室的高度設(shè)計(jì)過(guò)高(馬丁公司的無(wú)動(dòng)力除塵的沉降室高度約為 900 mm);延長(zhǎng)沉降室長(zhǎng)度可以有效抑塵,但當(dāng)?shù)竭_(dá)一定長(zhǎng)度后再延長(zhǎng)效果不明顯。

    圖3 沉降室高度變化的仿真模型[17]Fig.3 Simulation model of height variation of dust settling chamber

    (2)直接 DEM-CFD 耦合法 Goniva 等人[5]對(duì)Chen 的同一試驗(yàn)結(jié)果采用開(kāi)源 CFD 軟件和 LIGGGHTS的 DEM 軟件直接耦合,對(duì) Chen 的試驗(yàn)?zāi)P瓦M(jìn)行仿真,給出計(jì)算結(jié)果。該方法由于采用高效的 CFD-DEM 耦合,在20 個(gè)2.1 GHz 的CPU 上進(jìn)行計(jì)算,模擬時(shí)間為1 s,計(jì)算時(shí)間約為3.5 h,顆粒數(shù)約為150 000 個(gè)。

    3.3 不同仿真方法比較

    不同仿真方法的比較如表2 所列[17]。直接耦合的 CFD-DEM 法相比其他方法更有競(jìng)爭(zhēng)優(yōu)勢(shì),因?yàn)樵摲椒梢苑治鲱w粒物料在轉(zhuǎn)載溜槽的流動(dòng)性以及對(duì)周?chē)諝獾挠绊?,但是,更?fù)雜的數(shù)值模型和更長(zhǎng)的編程時(shí)間限制了該方法的應(yīng)用,其準(zhǔn)確性介于計(jì)算時(shí)間和使用難度之間。兩相流體模型 CFD(TFM)方法是最實(shí)用的方法,可幫助設(shè)計(jì)溜槽。采用沒(méi)有粉塵的DEM 和 CFD 耦合的方法,可以對(duì)物料通過(guò)轉(zhuǎn)載溜槽時(shí)的誘導(dǎo)空氣量進(jìn)行預(yù)測(cè)。

    表2 不同仿真方法比較Tab.2 Comparison of various simulation methods

    4 無(wú)動(dòng)力抑塵的幾種結(jié)構(gòu)設(shè)計(jì)分析

    4.1 引風(fēng)管設(shè)計(jì)

    在一些無(wú)動(dòng)力抑塵轉(zhuǎn)載溜槽系統(tǒng)中設(shè)置引風(fēng)管路,其目的是將空氣導(dǎo)流到落料管中,以減少導(dǎo)料槽出口處的空氣流量。而實(shí)際上設(shè)置引風(fēng)管后的效果可以在不考慮粉塵的條件下進(jìn)行仿真,這里采用不包含粉塵的 DEM 和單相 CFD 耦合方法進(jìn)行仿真。所建立的轉(zhuǎn)載溜槽系統(tǒng)模型如圖4 所示,仿真 20 s 后的氣流速度場(chǎng)分布如圖5 所示。分別對(duì)直徑為 300、400、500、600 mm 的引風(fēng)管進(jìn)行仿真,經(jīng)過(guò)引風(fēng)管的空氣流量?jī)H占從導(dǎo)料槽出口流出空氣的 5% 以下,這說(shuō)明引風(fēng)管結(jié)構(gòu)基本上不起作用,這與文獻(xiàn)[3]所得結(jié)果一致。在文獻(xiàn)[3]中,引風(fēng)管連接到落料管的位置比圖4 的位置還要高,這也說(shuō)明引風(fēng)管位置幾乎不起作用。

    圖4 引風(fēng)管抑塵轉(zhuǎn)載系統(tǒng)模型Fig.4 Model of dust suppression transfer system with inducing pipe

    圖5 20 s 時(shí)氣流速度場(chǎng)分布Fig.5 Distribution of air flow velocity at 20th s

    4.2 導(dǎo)料槽中的粉塵沉降

    通常顆粒沉降速度[18]

    式中:ρL為空氣密度,kg/m3;η為空氣動(dòng)力黏性系數(shù),常溫下一般為 1.8×10-5Pa/s。

    可見(jiàn)沉降室高度過(guò)高可能造成粉塵顆粒需要更長(zhǎng)的沉降時(shí)間,從而要求更長(zhǎng)的沉降室。導(dǎo)料槽中的粉塵沉降規(guī)律研究對(duì)于導(dǎo)料槽的長(zhǎng)度以及截面設(shè)計(jì)很有意義。在假設(shè)導(dǎo)料槽入口處的空氣流速前提下,可以采用 CFD 的多相流(DPM)模型進(jìn)行仿真,仿真中的顆粒分布如表3 所列。

    表3 仿真顆粒的分散度Tab.3 Dispersion degree of simulated particles

    對(duì)粒徑為 10~90 μm 的粉塵顆粒的沉降情況進(jìn)行對(duì)比仿真。在 FLUENT 的輸入?yún)?shù)中將粉塵顆粒的初始速度設(shè)為 0,入口速度設(shè)為 2 m/s,在 DPM 模型中粉塵粒徑的設(shè)定選擇 uniform 模型。仿真結(jié)束后觀(guān)察粉塵顆粒在風(fēng)流帶動(dòng)下的運(yùn)動(dòng)及沉降情況。

    圖6 中從上往下依次是粒徑為 10、20、30、40、50、60、70、80、90 μm 的粉塵顆粒在氣流帶動(dòng)下的沉降情況。從圖6 可以看出,當(dāng)粉塵顆粒在10 μm 以下時(shí),在氣流作用下不易沉降,會(huì)漂浮在導(dǎo)料槽中;當(dāng)粉塵顆粒在20 μm 以上時(shí),逐漸沉降在導(dǎo)料槽中,粉塵粒子的漂移距離隨粒子直徑增大而減小;當(dāng)粒徑大于 90 μm 時(shí),粉塵粒子基本都會(huì)在達(dá)到導(dǎo)料槽出口前降到受料輸送帶上。因此,導(dǎo)料槽防塵和降塵的重點(diǎn)應(yīng)該放在小顆粒粉塵上。盡管以上仿真可能存在參數(shù)選取問(wèn)題,但從定性上可以得出結(jié)論:當(dāng)導(dǎo)料槽設(shè)計(jì)得足夠長(zhǎng),僅有20 μm 以下的粉塵顆粒會(huì)隨空氣流出,對(duì)于這樣的粉塵需要采用干霧抑塵方法進(jìn)行降塵。

    圖6 10~90 μm 粉塵顆粒沉降情況Fig.6 Settlement status of dust particles of 10~90 μm

    4.3 空氣刮刀抑塵的 CFD 仿真

    空氣刮刀抑塵(Airscrape)是由德國(guó)工程師根據(jù)現(xiàn)場(chǎng)維護(hù)經(jīng)驗(yàn)發(fā)明的,其基本原理是在導(dǎo)料槽外側(cè)設(shè)置傾斜布置(與輸送帶運(yùn)行方向存在銳角)的空氣刮刀,且導(dǎo)料槽裙板和空氣刮刀均不與輸送帶接觸,如圖7所示。在輸送帶運(yùn)行過(guò)程中引起的氣流,通過(guò)傾斜空氣刮刀將導(dǎo)料槽外的空氣和粉塵吸進(jìn)導(dǎo)料槽內(nèi),圖7給出了空氣流動(dòng)方向。可以使用單相流 CFD 仿真,通過(guò)空氣壓力分布得到氣流流向的定性分析結(jié)果,其仿真結(jié)果如圖8 所示。整個(gè)導(dǎo)料槽模型內(nèi)部的壓力均小于 0,內(nèi)部完全是負(fù)壓。導(dǎo)料槽負(fù)壓狀態(tài)能夠保證導(dǎo)料槽內(nèi)的粉塵在內(nèi)外壓力差的作用下不會(huì)向外噴發(fā)。同時(shí),刮刀和輸送帶的0.5 mm 間距處的壓力比模型內(nèi)部壓力值更小,因此在更大的壓力差作用下,模型外部的空氣更容易通過(guò)0.5 mm 的間距涌入到模型中,因而能夠使轉(zhuǎn)載過(guò)程中的粉塵在空氣吸力的作用下保留在輸送帶上隨物料轉(zhuǎn)載出去。

    圖7 空氣刮刀布置簡(jiǎn)圖Fig.7 Layout sketch of air scraper

    圖8 空氣刮刀的空氣壓力分布Fig.8 Distribution of air pressure of air scraper

    4.4 生成粉塵量估算方法

    文獻(xiàn)[17]給出了散狀物料的儲(chǔ)存、轉(zhuǎn)載和運(yùn)輸過(guò)程中的粉塵排放量估算方法,粉塵排放量

    式中:vab為物料卸料到?jīng)_擊點(diǎn)的速度,m/s;qnorm為標(biāo)準(zhǔn)排放因子;a為物料粉塵化的加權(quán)因子,由表4確定。

    表4 物料粉塵化的加權(quán)因子Tab.4 Weighting factor of material dusting

    該排放量是物料通過(guò)轉(zhuǎn)載溜槽系統(tǒng)的最大排放量,而且在輸送系統(tǒng)中大部分的粉塵可能在物料下部隨其一起運(yùn)動(dòng),不能用于評(píng)價(jià)系統(tǒng)。Katterfeld[4]給出算例:沖擊點(diǎn)的物料速度為 6.3 m/s,輸送量為 774 t/h,低粉塵生成的加權(quán)因子為a=,標(biāo)準(zhǔn)排放因子qnorm=97.7 g·m3/t2,轉(zhuǎn)載過(guò)程生成粉塵qdust=48 g/t,每小時(shí)粉塵量qdust=37 kg/h。

    5 轉(zhuǎn)載溜槽的抑塵設(shè)計(jì)

    根據(jù)上述分析,可以得出轉(zhuǎn)載溜槽抑塵設(shè)計(jì)要點(diǎn):

    (1)通過(guò)在入口處設(shè)置密封擋簾來(lái)盡量縮小轉(zhuǎn)載溜槽的物料進(jìn)口,同時(shí)需要考慮入口其他部位(回程段)的密封;

    (2)通過(guò)設(shè)置 U 形、圓形截面的導(dǎo)流板和曲線(xiàn)溜槽來(lái)控制料流的導(dǎo)向,使料流集中,并減少物料與空氣的接觸面;

    (3)粉塵擴(kuò)散的關(guān)鍵是物料流經(jīng)落料管下落后沖擊造成的,為避免沖擊,物料應(yīng)落到下部曲線(xiàn)溜槽后,在其導(dǎo)向作用下使物料的速度接近受料帶式輸送機(jī)的速度;

    (4)適當(dāng)加大導(dǎo)料槽的長(zhǎng)度和高度構(gòu)成沉降室,以減緩空氣的流動(dòng)速度,即將空氣總壓力中的靜壓增加、動(dòng)壓減小,以利于降低粉塵的流動(dòng)速度,導(dǎo)料槽密封空間高度一般為:物料外輪廓最大高度的約 1.5倍,將裙板系統(tǒng)延伸到物料與輸送帶接觸區(qū)前部至少10 m 處;

    (5)在沉降室內(nèi)加設(shè)多道擋簾以起到阻尼的作用,進(jìn)一步減小空氣的動(dòng)壓;

    (6)在受料點(diǎn)附近設(shè)置緩沖床并使導(dǎo)料槽密封,避免粉塵的逃逸;

    (7)導(dǎo)料槽的后部也需適當(dāng)加長(zhǎng)并在后部密封。

    圖9 所示為導(dǎo)料槽的沉降室和擋簾的布置。

    圖9 導(dǎo)料槽的沉降室和擋簾布置Fig.9 Settling chamber and curtain layout of guide chute

    6 結(jié)論

    對(duì)轉(zhuǎn)載溜槽抑塵設(shè)計(jì)的關(guān)鍵問(wèn)題誘導(dǎo)空氣量的幾種計(jì)算模型進(jìn)行了分析對(duì)比,分析了轉(zhuǎn)載溜槽抑塵設(shè)計(jì)的計(jì)算機(jī)仿真方法,對(duì)引風(fēng)管設(shè)計(jì)采用 DEM-CFD方法進(jìn)行了仿真,對(duì)導(dǎo)料槽沉降室中的粉塵沉降采用CFD 兩相流仿真,對(duì)空氣刮刀設(shè)計(jì)采用 CFD 仿真,給出了轉(zhuǎn)載溜槽的設(shè)計(jì)原則,所得結(jié)論如下:

    (1)基于模型假設(shè)的轉(zhuǎn)載溜槽的誘導(dǎo)空氣估計(jì)方法所得結(jié)果偏大,且 CEMA 空氣量的計(jì)算中所包含的置換空氣在穩(wěn)定運(yùn)行狀態(tài)下是不存在的;

    (2)通過(guò) DEM-CFD 耦合仿真,證明采用引風(fēng)管設(shè)計(jì)對(duì)減小導(dǎo)料槽的空氣流量意義不大,導(dǎo)料槽沉降室應(yīng)具有足夠長(zhǎng)度,但過(guò)高的沉降室高度是不必要的,驗(yàn)證了采用空氣刮刀設(shè)計(jì)可以將空氣流導(dǎo)入料槽;

    (3)所總結(jié)的轉(zhuǎn)載溜槽的抑塵設(shè)計(jì)要點(diǎn)可以有效降低轉(zhuǎn)載溜槽的粉塵逃逸,有經(jīng)驗(yàn)表明可以降低轉(zhuǎn)載溜槽 98% 的粉塵。

    猜你喜歡
    方法模型設(shè)計(jì)
    一半模型
    重要模型『一線(xiàn)三等角』
    重尾非線(xiàn)性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    瞞天過(guò)?!律O(shè)計(jì)萌到家
    設(shè)計(jì)秀
    海峽姐妹(2017年7期)2017-07-31 19:08:17
    有種設(shè)計(jì)叫而專(zhuān)
    Coco薇(2017年5期)2017-06-05 08:53:16
    可能是方法不對(duì)
    3D打印中的模型分割與打包
    用對(duì)方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    国产人伦9x9x在线观看| 国产成人免费观看mmmm| 国产三级黄色录像| 中文字幕制服av| 啦啦啦视频在线资源免费观看| 一区在线观看完整版| 日韩三级视频一区二区三区| 亚洲精品国产av蜜桃| 久久 成人 亚洲| 中文字幕人妻熟女乱码| 亚洲精品国产精品久久久不卡| 亚洲全国av大片| 欧美 日韩 精品 国产| 男男h啪啪无遮挡| 99久久精品国产亚洲精品| 久久人妻熟女aⅴ| 国产主播在线观看一区二区| 一边摸一边做爽爽视频免费| 我的亚洲天堂| 黄色片一级片一级黄色片| 国产精品自产拍在线观看55亚洲 | 欧美黑人欧美精品刺激| 国产一级毛片在线| 精品高清国产在线一区| 日韩大片免费观看网站| 欧美精品人与动牲交sv欧美| 美女午夜性视频免费| 九色亚洲精品在线播放| 精品一区在线观看国产| 国产一区二区三区av在线| 成人亚洲精品一区在线观看| av在线老鸭窝| 99热全是精品| 亚洲精品av麻豆狂野| 黄色片一级片一级黄色片| 男女免费视频国产| 一区二区日韩欧美中文字幕| 一级黄色大片毛片| 精品国产乱子伦一区二区三区 | 九色亚洲精品在线播放| 高潮久久久久久久久久久不卡| 老司机午夜十八禁免费视频| 黄频高清免费视频| 国产免费av片在线观看野外av| 欧美 日韩 精品 国产| 中文字幕精品免费在线观看视频| 2018国产大陆天天弄谢| 久久久久国内视频| 女人久久www免费人成看片| 亚洲欧美色中文字幕在线| 黑人巨大精品欧美一区二区mp4| 啦啦啦中文免费视频观看日本| 成人国产av品久久久| 99国产极品粉嫩在线观看| 女人爽到高潮嗷嗷叫在线视频| 女人精品久久久久毛片| 精品福利永久在线观看| 国产精品久久久人人做人人爽| a 毛片基地| 国产日韩一区二区三区精品不卡| 99香蕉大伊视频| 91成人精品电影| 在线十欧美十亚洲十日本专区| 欧美另类亚洲清纯唯美| 亚洲av国产av综合av卡| 精品人妻熟女毛片av久久网站| 久久人人97超碰香蕉20202| 咕卡用的链子| 午夜久久久在线观看| 大香蕉久久成人网| 色94色欧美一区二区| 搡老岳熟女国产| 97精品久久久久久久久久精品| av在线播放精品| 最新在线观看一区二区三区| av网站免费在线观看视频| 亚洲成av片中文字幕在线观看| 成人av一区二区三区在线看 | 啦啦啦免费观看视频1| 激情视频va一区二区三区| 久久 成人 亚洲| 亚洲精品国产一区二区精华液| 国产一区二区三区av在线| 免费观看人在逋| 久久精品aⅴ一区二区三区四区| 纯流量卡能插随身wifi吗| 日韩免费高清中文字幕av| 男女高潮啪啪啪动态图| 别揉我奶头~嗯~啊~动态视频 | 精品第一国产精品| 黄色片一级片一级黄色片| av超薄肉色丝袜交足视频| 成人免费观看视频高清| av欧美777| 亚洲精品国产av蜜桃| 国产精品久久久久久精品电影小说| 久久99热这里只频精品6学生| 成年女人毛片免费观看观看9 | 免费av中文字幕在线| bbb黄色大片| 亚洲精品自拍成人| 国内毛片毛片毛片毛片毛片| 欧美黑人欧美精品刺激| 精品一品国产午夜福利视频| 欧美日韩亚洲高清精品| 97人妻天天添夜夜摸| 99久久99久久久精品蜜桃| 91大片在线观看| av又黄又爽大尺度在线免费看| 国产熟女午夜一区二区三区| 纵有疾风起免费观看全集完整版| 成人亚洲精品一区在线观看| 国产伦理片在线播放av一区| 午夜免费成人在线视频| 麻豆乱淫一区二区| 热99国产精品久久久久久7| 啦啦啦中文免费视频观看日本| 正在播放国产对白刺激| 日韩欧美一区二区三区在线观看 | 午夜福利,免费看| 精品国产一区二区三区四区第35| 亚洲欧美色中文字幕在线| 日韩精品免费视频一区二区三区| 女人高潮潮喷娇喘18禁视频| 亚洲精品在线美女| 狂野欧美激情性bbbbbb| 我要看黄色一级片免费的| 免费一级毛片在线播放高清视频 | 日韩大片免费观看网站| 操出白浆在线播放| 久久久久国内视频| 国产伦人伦偷精品视频| 国产精品九九99| 亚洲国产精品999| 日本欧美视频一区| av线在线观看网站| 亚洲精品中文字幕在线视频| 一级黄色大片毛片| 老熟妇乱子伦视频在线观看 | 在线亚洲精品国产二区图片欧美| 亚洲五月色婷婷综合| 99九九在线精品视频| 午夜免费鲁丝| 国产高清国产精品国产三级| 午夜激情久久久久久久| 少妇人妻久久综合中文| 一区在线观看完整版| 久久久久国内视频| 电影成人av| 高清黄色对白视频在线免费看| 黑丝袜美女国产一区| 色94色欧美一区二区| 岛国毛片在线播放| 免费在线观看完整版高清| 热99久久久久精品小说推荐| 国产主播在线观看一区二区| 国产真人三级小视频在线观看| 午夜福利视频在线观看免费| 亚洲av成人不卡在线观看播放网 | 中国美女看黄片| 久久99一区二区三区| 免费观看人在逋| 亚洲九九香蕉| 亚洲av电影在线进入| 国产成人精品在线电影| 欧美日韩中文字幕国产精品一区二区三区 | 99国产精品99久久久久| 免费一级毛片在线播放高清视频 | 高清视频免费观看一区二区| 久久精品aⅴ一区二区三区四区| 精品一区二区三卡| 男女无遮挡免费网站观看| 日韩 欧美 亚洲 中文字幕| 国产成人系列免费观看| av在线老鸭窝| 欧美精品高潮呻吟av久久| 国产成人系列免费观看| 天天影视国产精品| 美女午夜性视频免费| 一二三四社区在线视频社区8| 国产一卡二卡三卡精品| 精品乱码久久久久久99久播| 亚洲视频免费观看视频| 久久天堂一区二区三区四区| 国产一卡二卡三卡精品| a级片在线免费高清观看视频| 亚洲视频免费观看视频| a级片在线免费高清观看视频| 亚洲第一青青草原| 欧美黄色片欧美黄色片| 黄片小视频在线播放| 十八禁网站免费在线| 国产精品九九99| 国产精品一区二区在线观看99| 亚洲av日韩精品久久久久久密| 美女高潮到喷水免费观看| 久久精品人人爽人人爽视色| 精品久久久久久久毛片微露脸 | 久久精品亚洲av国产电影网| 如日韩欧美国产精品一区二区三区| 久久九九热精品免费| 如日韩欧美国产精品一区二区三区| 麻豆国产av国片精品| 18禁黄网站禁片午夜丰满| 亚洲欧美成人综合另类久久久| 极品人妻少妇av视频| 18禁裸乳无遮挡动漫免费视频| 91麻豆av在线| 男女国产视频网站| 青春草亚洲视频在线观看| 日本av手机在线免费观看| 久久国产亚洲av麻豆专区| 在线观看免费高清a一片| 女人爽到高潮嗷嗷叫在线视频| 在线十欧美十亚洲十日本专区| 97人妻天天添夜夜摸| 国产亚洲一区二区精品| 久久久久久久大尺度免费视频| 亚洲国产成人一精品久久久| 人人妻,人人澡人人爽秒播| 波多野结衣av一区二区av| 一区福利在线观看| 他把我摸到了高潮在线观看 | 少妇被粗大的猛进出69影院| 欧美黑人欧美精品刺激| 国产三级黄色录像| 中文字幕另类日韩欧美亚洲嫩草| 亚洲一码二码三码区别大吗| 亚洲精品日韩在线中文字幕| 亚洲色图 男人天堂 中文字幕| 黄片大片在线免费观看| 一区二区三区精品91| 后天国语完整版免费观看| 国产欧美日韩综合在线一区二区| 国产成人av教育| 亚洲成国产人片在线观看| 熟女少妇亚洲综合色aaa.| 中文字幕最新亚洲高清| 啦啦啦在线免费观看视频4| 美国免费a级毛片| 9热在线视频观看99| 国产无遮挡羞羞视频在线观看| 老熟妇仑乱视频hdxx| 女人高潮潮喷娇喘18禁视频| 国产淫语在线视频| 亚洲久久久国产精品| 91成年电影在线观看| 亚洲欧美一区二区三区黑人| 亚洲av国产av综合av卡| 日韩有码中文字幕| 18禁裸乳无遮挡动漫免费视频| 高清在线国产一区| 美女福利国产在线| 99国产精品一区二区三区| 一边摸一边做爽爽视频免费| 亚洲精品美女久久久久99蜜臀| 美女大奶头黄色视频| 正在播放国产对白刺激| 黄色视频在线播放观看不卡| 久久久欧美国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久精品精品| 激情视频va一区二区三区| 人妻一区二区av| 久久久久国内视频| 精品欧美一区二区三区在线| 国产老妇伦熟女老妇高清| 欧美精品高潮呻吟av久久| 国产精品九九99| 日韩一卡2卡3卡4卡2021年| 色婷婷久久久亚洲欧美| 久久精品国产亚洲av香蕉五月 | 一级a爱视频在线免费观看| 久久精品国产亚洲av香蕉五月 | 新久久久久国产一级毛片| 大型av网站在线播放| √禁漫天堂资源中文www| 丰满饥渴人妻一区二区三| 久久久水蜜桃国产精品网| 夜夜夜夜夜久久久久| 亚洲va日本ⅴa欧美va伊人久久 | 国产精品偷伦视频观看了| 中文字幕高清在线视频| 欧美在线黄色| 国产色视频综合| 亚洲欧美一区二区三区久久| 他把我摸到了高潮在线观看 | 久久国产精品大桥未久av| 狠狠婷婷综合久久久久久88av| 高清黄色对白视频在线免费看| 伦理电影免费视频| 建设人人有责人人尽责人人享有的| 在线天堂中文资源库| 国产精品亚洲av一区麻豆| 一边摸一边抽搐一进一出视频| 欧美午夜高清在线| 男人操女人黄网站| 人妻人人澡人人爽人人| h视频一区二区三区| 十八禁人妻一区二区| 91九色精品人成在线观看| 一二三四在线观看免费中文在| 五月开心婷婷网| 亚洲精品国产一区二区精华液| 手机成人av网站| 蜜桃在线观看..| 久久久国产成人免费| 9191精品国产免费久久| 久久狼人影院| 久久久久久久久久久久大奶| 黄色怎么调成土黄色| 黑人巨大精品欧美一区二区蜜桃| 91成年电影在线观看| 欧美中文综合在线视频| 婷婷成人精品国产| 美女高潮到喷水免费观看| 国产男人的电影天堂91| 国产亚洲av高清不卡| 国产欧美日韩精品亚洲av| 欧美日韩一级在线毛片| 人妻人人澡人人爽人人| 久久精品熟女亚洲av麻豆精品| 国产成人影院久久av| 啦啦啦视频在线资源免费观看| 在线观看免费视频网站a站| 后天国语完整版免费观看| 咕卡用的链子| 人妻 亚洲 视频| tocl精华| 午夜福利影视在线免费观看| 午夜免费观看性视频| 国产成人一区二区三区免费视频网站| 老司机深夜福利视频在线观看 | e午夜精品久久久久久久| 黑人巨大精品欧美一区二区mp4| 亚洲中文av在线| 一级片'在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美精品自产自拍| 精品国产乱子伦一区二区三区 | 18禁裸乳无遮挡动漫免费视频| 精品少妇一区二区三区视频日本电影| 国产精品二区激情视频| 少妇人妻久久综合中文| 亚洲一区中文字幕在线| 99国产精品免费福利视频| 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| 欧美日韩中文字幕国产精品一区二区三区 | 中文字幕高清在线视频| 一区二区日韩欧美中文字幕| 欧美人与性动交α欧美精品济南到| 免费日韩欧美在线观看| 免费在线观看完整版高清| 亚洲精品乱久久久久久| 777久久人妻少妇嫩草av网站| 在线天堂中文资源库| 国产区一区二久久| 9热在线视频观看99| 19禁男女啪啪无遮挡网站| 各种免费的搞黄视频| 国产亚洲av高清不卡| 91成年电影在线观看| 伦理电影免费视频| 日韩欧美国产一区二区入口| 99国产精品一区二区蜜桃av | 国产福利在线免费观看视频| 国产精品一区二区免费欧美 | 国产一区二区在线观看av| 99热全是精品| 俄罗斯特黄特色一大片| 丰满迷人的少妇在线观看| 亚洲国产毛片av蜜桃av| 国产黄频视频在线观看| 在线观看舔阴道视频| 美女高潮喷水抽搐中文字幕| 国产97色在线日韩免费| 少妇 在线观看| 国产有黄有色有爽视频| 久久国产精品人妻蜜桃| 少妇精品久久久久久久| 中亚洲国语对白在线视频| 久久国产亚洲av麻豆专区| 久久精品aⅴ一区二区三区四区| 国产精品秋霞免费鲁丝片| 久久久精品区二区三区| 满18在线观看网站| 久久 成人 亚洲| 亚洲视频免费观看视频| 国产精品秋霞免费鲁丝片| 日韩一区二区三区影片| 欧美日韩亚洲综合一区二区三区_| 亚洲av电影在线观看一区二区三区| 午夜免费成人在线视频| 中亚洲国语对白在线视频| 亚洲欧美一区二区三区久久| 免费在线观看完整版高清| 成人国语在线视频| 超碰97精品在线观看| 香蕉丝袜av| 青青草视频在线视频观看| 一个人免费看片子| 后天国语完整版免费观看| 欧美性长视频在线观看| www.999成人在线观看| 777米奇影视久久| 中文字幕制服av| a级毛片在线看网站| 在线观看舔阴道视频| 久久天堂一区二区三区四区| 桃花免费在线播放| 黄网站色视频无遮挡免费观看| 欧美激情高清一区二区三区| 精品国内亚洲2022精品成人 | 免费在线观看黄色视频的| 成在线人永久免费视频| 老司机影院毛片| 丰满人妻熟妇乱又伦精品不卡| 999精品在线视频| 天堂中文最新版在线下载| 免费在线观看日本一区| 国产xxxxx性猛交| bbb黄色大片| 美国免费a级毛片| 国产在线视频一区二区| 国产男人的电影天堂91| videos熟女内射| 成人手机av| 亚洲va日本ⅴa欧美va伊人久久 | 久久久久国内视频| 欧美人与性动交α欧美软件| 91精品三级在线观看| 免费久久久久久久精品成人欧美视频| 12—13女人毛片做爰片一| 日韩人妻精品一区2区三区| 国产成人精品久久二区二区免费| 亚洲精品一二三| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看 | 国产亚洲一区二区精品| 欧美xxⅹ黑人| 亚洲av片天天在线观看| 精品人妻熟女毛片av久久网站| 大片电影免费在线观看免费| 黑人巨大精品欧美一区二区蜜桃| 成人手机av| 黄色视频在线播放观看不卡| 黑丝袜美女国产一区| 精品国产乱子伦一区二区三区 | 精品欧美一区二区三区在线| 18禁观看日本| 中国美女看黄片| www日本在线高清视频| 黑丝袜美女国产一区| 人人澡人人妻人| 亚洲精品久久成人aⅴ小说| 欧美日韩一级在线毛片| 动漫黄色视频在线观看| a级毛片在线看网站| 亚洲专区字幕在线| 欧美日本中文国产一区发布| 国产精品 欧美亚洲| 亚洲国产毛片av蜜桃av| 操美女的视频在线观看| 三上悠亚av全集在线观看| 亚洲av国产av综合av卡| 精品高清国产在线一区| 国产精品二区激情视频| 国产精品欧美亚洲77777| 丁香六月天网| tocl精华| 日韩 亚洲 欧美在线| 日本a在线网址| 搡老熟女国产l中国老女人| 90打野战视频偷拍视频| 两个人看的免费小视频| 精品人妻在线不人妻| 可以免费在线观看a视频的电影网站| 不卡一级毛片| 男女边摸边吃奶| 午夜福利影视在线免费观看| 久久久久精品人妻al黑| 亚洲avbb在线观看| tube8黄色片| 亚洲国产av新网站| 亚洲国产成人一精品久久久| 欧美成狂野欧美在线观看| 少妇 在线观看| 精品福利永久在线观看| 国产亚洲精品一区二区www | 亚洲中文日韩欧美视频| 国产伦理片在线播放av一区| 欧美另类亚洲清纯唯美| 欧美中文综合在线视频| 久久精品国产a三级三级三级| 在线观看免费午夜福利视频| 精品少妇一区二区三区视频日本电影| 国产成人精品无人区| 丰满少妇做爰视频| 精品久久久久久久毛片微露脸 | 妹子高潮喷水视频| 国产黄色免费在线视频| 曰老女人黄片| 亚洲精品中文字幕一二三四区 | 老司机亚洲免费影院| 久久中文字幕一级| 欧美日韩国产mv在线观看视频| 午夜福利在线免费观看网站| 日日爽夜夜爽网站| 国产欧美日韩精品亚洲av| 天堂8中文在线网| 午夜福利,免费看| 淫妇啪啪啪对白视频 | 久久香蕉激情| 纯流量卡能插随身wifi吗| 欧美人与性动交α欧美精品济南到| 老汉色∧v一级毛片| 一区二区三区精品91| 一二三四在线观看免费中文在| 久久精品成人免费网站| 黑人巨大精品欧美一区二区mp4| 嫁个100分男人电影在线观看| 国产亚洲一区二区精品| 久久性视频一级片| 欧美成人午夜精品| 欧美在线一区亚洲| 久久国产精品影院| 大香蕉久久成人网| 老汉色∧v一级毛片| 免费高清在线观看视频在线观看| 97精品久久久久久久久久精品| 亚洲精品久久久久久婷婷小说| 一区二区三区激情视频| 成人手机av| 性色av一级| 麻豆国产av国片精品| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩综合在线一区二区| 嫩草影视91久久| av片东京热男人的天堂| 国产亚洲精品第一综合不卡| 欧美黄色淫秽网站| 精品亚洲成a人片在线观看| 国产97色在线日韩免费| 亚洲国产毛片av蜜桃av| 一级片'在线观看视频| h视频一区二区三区| 男女国产视频网站| 18禁国产床啪视频网站| 制服人妻中文乱码| 黄片小视频在线播放| 91老司机精品| 欧美日本中文国产一区发布| 91精品国产国语对白视频| 黄色视频,在线免费观看| 亚洲自偷自拍图片 自拍| 亚洲国产精品一区二区三区在线| 黄色视频在线播放观看不卡| 国产成人影院久久av| 每晚都被弄得嗷嗷叫到高潮| 99香蕉大伊视频| 亚洲少妇的诱惑av| 精品久久蜜臀av无| 国产精品久久久人人做人人爽| 两个人看的免费小视频| 日韩三级视频一区二区三区| 中文字幕高清在线视频| 成在线人永久免费视频| 99国产精品一区二区三区| 国产精品 国内视频| 日韩中文字幕欧美一区二区| 99精品久久久久人妻精品| 日韩欧美国产一区二区入口| 在线av久久热| 青春草视频在线免费观看| 国产高清videossex| 日本黄色日本黄色录像| 色播在线永久视频| 久久精品国产亚洲av高清一级| 国产精品一区二区精品视频观看| 亚洲一区二区三区欧美精品| 亚洲精品国产区一区二| 热99re8久久精品国产| 在线观看一区二区三区激情| 国产精品 欧美亚洲| 国产精品国产av在线观看| 高潮久久久久久久久久久不卡| 久久女婷五月综合色啪小说| 一本综合久久免费| 国产精品 国内视频| 麻豆乱淫一区二区| 亚洲人成电影免费在线| 午夜福利视频在线观看免费| 亚洲精品国产一区二区精华液| 18禁国产床啪视频网站| 真人做人爱边吃奶动态| 欧美激情高清一区二区三区| 9热在线视频观看99| 两性夫妻黄色片| 自线自在国产av| 欧美一级毛片孕妇| 国产不卡av网站在线观看| 久久青草综合色| 久久中文字幕一级| 97精品久久久久久久久久精品| tocl精华| videos熟女内射| 高清在线国产一区| 国产一区二区激情短视频 | av有码第一页| 精品久久久精品久久久| 亚洲性夜色夜夜综合| 老司机深夜福利视频在线观看 | 国产精品99久久99久久久不卡| a在线观看视频网站|