999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

葉頂間隙大小對螺旋離心泵內部壓力脈動的影響

2017-01-17 15:14:44王炳祺史長彪沈登武
農業工程學報 2017年1期

張 華,陳 斌,王炳祺,史長彪,沈登武

(1. 合肥工業大學機械工程學院,合肥 230009;2. 江蘇大學流體機械工程技術研究中心,鎮江 212013;3. 藍深集團股份有限公司,南京 211500)

葉頂間隙大小對螺旋離心泵內部壓力脈動的影響

張 華1,2,陳 斌1※,王炳祺1,史長彪3,沈登武3

(1. 合肥工業大學機械工程學院,合肥 230009;2. 江蘇大學流體機械工程技術研究中心,鎮江 212013;3. 藍深集團股份有限公司,南京 211500)

葉頂間隙大小是影響半開式螺旋離心泵內外特性的重要因素之一,為充分探究葉頂間隙大小對其內外特性的影響程度及作用機理,該文以一臺比轉速237的單葉片螺旋離心泵為研究對象,通過設計葉頂間隙專用調節機構,分別調整葉頂間隙大小至0.3、0.5和0.8 mm,與葉輪外徑的比值分別為0.13%、0.22%和0.35%。針對3種葉頂間隙大小情況,同時進行外特性試驗和壓力脈動試驗研究。壓力脈動的監測點包括葉輪進口邊、葉輪中部、蝸殼第三斷面和蝸殼隔舌附近。結果表明葉頂間隙大小與葉輪平均直徑的比值應介于0.13%~0.22%之間。同時,結合壓力脈動試驗獲得的實測壓力波形圖、壓力系數波形圖及頻域圖,分析了單葉片螺旋離心泵內部壓力脈動的規律,為半開式螺旋離心泵的減振降噪提供參考和借鑒。

離心泵;葉片;葉輪;螺旋葉輪;葉頂間隙;外特性;壓力脈動

0 引 言

單葉片螺旋離心葉輪具有較強的通過能力和抗纏繞能力,主要應用于無損運輸和污水輸送等領域[1-3]。該類型葉輪具有葉片包角大、扭曲程度高[4]等特點,多采用半開式設計[5],葉輪輪緣與底蓋板之間存在一定的間隙,該間隙不僅直接影響機組運行效率,對泵的內部流場物理結構及壓力脈動、振動等特性也有較大的影響[6-7]。施衛東等[8]以軸流泵為研究對象,討論了4種不同葉頂間隙大小對軸流泵外特性、葉輪出口軸面速度及環量的影響,發現葉頂間隙超過0.5 mm時,泵性能下降幅度明顯增大,并成功捕獲了葉頂泄漏渦。因此,研究葉頂間隙大小對單葉片螺旋離心泵性能的影響不僅具有實際工程意義,也具有較高的學術價值。單葉片螺旋離心葉輪由于非對稱葉片結構,在不同工況下,徑向力變化較大[9-10],國內外學者逐漸對其展開深入的研究[11-14]。Han等[15]通過數值計算的方法,研究了單葉片螺旋離心泵的軸向力和徑向力,并給出了極值情況下的葉輪與蝸殼相對位置。Quan等[16]對單葉片螺旋離心泵的能量傳遞機理進行了分析。Scherbatenko等[17]研究了螺旋離心葉輪出口的絕對速度分布。Pei等[18-20]采用流固耦合與試驗相結合的方法,研究了單葉片葉輪進口的振動及變形情況,并引入波動量標準差概念分析了內部壓力脈動強度。Nishi等[21-23]研究了不同單葉片出口安放角對葉輪徑向力的影響,并通過設計專用試驗測量裝置,獲得了單葉片葉輪所受徑向力的情況。

本文以一臺比轉速ns=237的單葉片螺旋離心泵為研究對象,通過設計葉頂間隙專用可調裝置,設定3種葉頂間隙大小,進行外特性試驗,分析泵外特性性能與葉頂間隙大小之間的映射關系。同時,對葉輪進口邊、葉輪中部、蝸殼第三斷面及蝸殼隔舌附近進行壓力脈動監測[24-25],通過時域圖和頻域圖的對比,分析葉頂間隙大小對單葉片螺旋離心泵內部脈動特性的影響。為單葉片螺旋離心泵的減振降噪提供參考和借鑒。

1 研究對象

1.1 水力設計

單葉片螺旋離心泵的設計流量(以下標“d”代表設計工況)Qd=190 m3/h,設計揚程Hd=9 m,轉速n=1 470 r/min,比轉速ns=237。葉輪采用保角變換法設計,葉輪與蝸殼的水力設計結果見表1。

1.2 結構設計

試驗樣機的結構設計如圖1所示,采用潛水排污泵機電一體化結構設計[26-27]。葉輪為錐度軸傳動,并通過螺釘緊固。

表1 葉輪及蝸殼設計結果Table 1 Impeller and volute design results

圖1 單葉片螺旋離心式潛水排污泵的結構總裝圖Fig.1 Structural assembly diagram of screw-centrifugal sewage pump with single blade

底蓋與蝸殼在徑向方向采用O型圈密封,并通過調節螺栓連接,葉頂間隙大小通過調節螺栓調整(此間隙可調結構已應用于實際工程,主要在水泵出廠安裝檢驗和維修保養時調節,在水泵正常運行時不會產生軸向移動,因此,不會對O型圈密封形成磨損。),其具體調節步驟為:

1)葉輪及蝸殼安裝到位后,將雙頭螺柱擰入蝸殼;

2)將調節螺栓擰入底蓋,使得調節螺栓的螺紋頭部伸出一定長度;

3)將底蓋裝到蝸殼上,安裝時注意使雙頭螺柱穿過調節螺栓的內孔;

4)將緊固螺母擰入雙頭螺柱的下端;

5)利用塞尺測量葉片頂部至底蓋的最小距離,調節螺栓控制底蓋至蝸殼的距離,從而控制葉頂間隙大小,當距離調節到需要大小時,擰緊緊固螺母。

2 試驗方法

在符合國家標準的2級精度開式試驗臺上進行外特性試驗。首先將試驗泵的葉頂間隙調整到較大值,水泵通電,點動試車,觀察空開電流及葉輪旋轉時聲音的變化,判斷葉輪與底蓋之間是否發生刮擦,并利用塞尺檢查葉頂間隙大小。最終,確定最小葉間隙為0.3 mm。將泵安裝上管路系統,從流量調節閥關死至敞開重復試驗3次,每次間隔采集13個數據點以上,在確定試驗的精度及可重復性后,確定試驗泵的最優工況點。

安裝高頻潛水壓力傳感器,對0.3、0.5、0.8 mm 3種葉頂間隙情況下(葉頂間隙大小與葉輪平均直徑的比值分別為0.13%、0.22%、和0.35%)的壓力脈動進行研究。傳感器從泵外部通過螺紋孔連接至底蓋上,并控制傳感器端面不超過底蓋內部錐面,保證對流場不造成破壞。每次僅測量8個相同或相近流量下的數據,即0、0.2、0.4、0.6、0.8、1.0、1.2和1.4Qd工況。壓力脈動數據采集儀為華中科技大學開發的HSJ-2010型號水力機械綜合測試儀[28]。壓力傳感器采樣頻率設置為8192,各工況采樣時間為30~50 s[29-30]。利用鍵相技術,完成對整周期數據樣本的分離。4只壓力傳感器的布置方式如圖2所示,傳感器Pre1安裝在葉輪進口邊附近,傳感器Pre2靠近葉輪中部,傳感器Pre3處于蝸殼第三斷面的中間截面上,傳感器Pre4為蝸殼隔舌附近。

圖2 四只壓力傳感器的布置位置Fig.2 Location of four pressure sensors

3 結果與分析

3.1 外特性結果

圖3為3種葉頂間隙大小情況下的外特性性能對比。

圖3 不同葉頂間隙外特性性能對比Fig.3 External characteristic comparison of different tip clearances

試驗中測量得到的最優工況(以下標“opt”代表最優工況)為:流量Qopt=196.9 m3/h,揚程Hopt=9.05 m,功率Popt=6.84 kW,效率ηopt=70.92%。從圖3中可以看出,葉頂間隙大小的改變對流量-揚程、流量-功率和流量-效率曲線的走勢規律影響不大,但隨著葉頂間隙的增加,揚程、功率及效率均出現一定程度的下降,通過數據統計發現,在間隙由0.3 mm變為0.5 mm時,1.0Qd工況下,揚程及效率分別下降0.74 m和3.6%;當間隙從0.5 mm變為0.8 mm時,1.0Qd工況下,揚程及效率分別下降1.16 m和5.3%。其中,0.3 mm為試驗所用泵在旋轉時不發生刮擦的最小間隙;當間隙超過0.5 mm時,泵性能下降幅度過大。因此,試驗研究所用螺旋離心泵的葉頂間隙大小應不超過0.5 mm,即葉頂間隙的許用范圍為0.3~0.5 mm,與葉輪平均直徑的比值范圍為0.13%~0.22%。

3.2 壓力脈動結果

定義壓力脈動系數P*為

式中Pi為任意時刻的實測壓力,kPa;Pave為所取樣本長度上的平均壓力值,kPa。

定義時間系數t*為

式中ti為任意時刻,s;T為葉輪轉過一周的時間長度,s。

定義頻率系數NF為

式中f為傅立葉變化求得的任意頻率值,Hz;F為葉輪的轉頻(因為葉片為1枚,因此,轉頻與葉頻相同),Hz。

對每個壓力傳感器獲得的數據中,截取50個整周期的數據進行傅立葉變化,獲得壓力脈動頻譜圖;并從該數據中截取2個整周期進行壓力脈動實測波形圖和壓力脈動系數波形圖的繪制。由于泵內部各點處的壓力大小不同,無法反應出各點的壓力脈動強度。因此,本文選擇壓力系數波形圖表述,壓力波形系數越平緩,代表壓力脈動誘導水力振動和噪聲的程度越小。

圖4為壓力傳感器Pre1處的測試及分析結果。從圖4a可以看出3種葉頂間隙大小情況下,實測壓力的波形走勢基本一致。在葉頂間隙大小為0.3和0.5 mm時,實測壓力波形圖上壓力的極大值與極小值分別相互接近,但葉頂間隙大小為0.8 mm時,壓力的極大值變小,而壓力的極小值變大,其主要原因為同一圓周角度上,葉片頂部的半徑最大,做功能力最強,隨著葉片的旋轉,當壓力面逼近測壓點時,所測壓力達到極大值。葉頂間隙增大時,會加劇葉片壓力面高壓流體向吸力面的泄漏程度,從而導致壓力極大值降低,壓力極小值增大。因此,葉頂間隙增大雖然導致外特性性能的下降,但對葉頂間隙內的汽蝕能夠起一定的抑制作用(受研究條件的限制,該文未對葉頂間隙內的汽蝕現象進行深入研究,尚不清楚葉頂間隙內的汽蝕程度,無法得出汽蝕程度與葉頂間隙大小之間的關系,但該方向將會是以后研究的重點方向之一)。從圖4a還可以看出,在葉頂間隙大小為0.3和0.5 mm時,時間系數t*值介于0~0.6范圍內,實測壓力值均小于0,結合圖4b和公式(1)可知,這2種葉頂間隙下,葉輪旋轉一周內的平均壓力值大于0。隨著葉頂間隙增加至0.8 mm,由于葉頂附近壓力泄漏嚴重,使得壓力的平均值增加,從而導致壓力系數波動分布最為平緩。從壓力脈動頻域圖中發現,葉頂間隙大小的變化對前幾階主頻的頻率值基本無影響,前四階主頻分別為1F、2F、3F和4F,從各階主頻對應的幅值分析,葉頂間隙大小的增加使得主頻的幅值逐漸降低,但對二階、三階、四階主頻的幅值影響不大。

圖4 傳感器Pre1所測壓力脈動結果Fig.4 Pressure fluctuation test results of sensor Pre1

圖5為壓力傳感器Pre2處的測試及分析結果。在實測壓力波形圖及壓力系數波形圖中,間隙為0.3與0.5 mm時,二者的曲線基本重合,當間隙增加至0.8 mm時,正斜率曲線段內壓力實測值出現明顯下降,實測壓力的平均值也隨之下降,但間隙的變化并未對實測壓力中的極值產生明顯影響,因此,在壓力系數波形圖中,0.8 mm間隙情況下,壓力系數的極大值最大,極小值最小。根據鍵相標記點和橫坐標時間軸位置,圖6給出了時間系數t*1=0.418和t*2=0.603時葉輪與傳感器Pre2(標記點“×”)之間的相對位置[31],t*1=0.418時對應了壓力極大值,可知當測壓點極度逼近葉片壓力面時,壓力達到最大值;t*2=0.603時對應壓力極小值,此時測壓點處于間隙內部。t3=0.750時,間隙剛好轉過測壓點,隨著葉頂間隙大小的增加,葉頂泄漏流動增強,泄漏流與主流發生卷吸形成泄漏渦,泄漏渦繞過測壓點,造成壓力脈動在經歷從極小值逐漸上升的過程中,出現壓力突然小幅下降后回升的現象,如圖5a中橢圓虛線區域所示。

從圖5c頻域圖中還可看出,除主頻幅值隨著間隙增加略有下降外,其他階次主頻的幅值基本不變。對比壓力傳感器Pre1位置,Pre2位置在對應葉輪倍頻下的幅值更大,且5~8整數倍倍頻被明顯捕捉,其主要原因為Pre2位置為葉輪中部區域,葉片的直徑大,做功能力強,且易受周邊葉片的影響。

圖5 傳感器Pre2所測壓力脈動結果Fig.5 Pressure fluctuation test results of sensor Pre2

圖6 、下葉輪與測壓點Pre2的相對位置Fig.6 Relative position between impeller and pressure sensor Pre2 atand

圖7、圖8所示分別給出了蝸殼第三斷面(傳感器Pre3位置)和蝸殼隔舌附近(傳感器Pre4位置)的壓力脈動測試及分析結果,兩者的實測壓力波形圖走勢相反,但2個監測點位置存在以下共性:1)0.3與0.5 mm情況下的實測壓力曲線基本重合,當間隙達到0.8 mm時,實測壓力曲線整體向下平移;2)3種間隙大小情況下的壓力系數波形基本重合,若以壓力系數來判斷壓力脈動程度,則間隙大小對蝸殼第三端面與隔舌基本沒有影響;3)1F為主頻,2F為次頻,當間隙從0.3 mm變化到0.5 mm時,主頻幅值稍有下降,當間隙從0.5 mm變化到0.8 mm時,主頻幅值基本穩定不變;次頻幅值基本不隨間隙大小而變化。

圖7 傳感器Pre3所測壓力脈動結果Fig.7 Pressure fluctuation test results of sensor Pre3

圖8 傳感器Pre4所測壓力脈動結果Fig.8 Pressure fluctuation test results of sensor Pre4

4 結 論

1)間隙由0.3 mm變為0.5 mm時,設計工況下,揚程及效率分別下降0.74 m和3.6%;當間隙從0.5 mm變為0.8 mm時,設計工況下,揚程及效率分別下降1.16 m和5.3%。因此,結合壓力脈動特性與外特性結果,針對本文試驗用泵,認為可允許的間隙范圍為0.3~0.5 mm,與葉輪平均直徑的比值范圍為0.13%~0.22%,為以后的研究提供了參考和借鑒。

2)葉輪進口邊附近受負壓作用較大,增大葉頂間隙,壓力泄漏嚴重,增加了周期內底蓋表面監測點處的平均壓力,導致0.8 mm情況下,壓力系數波形圖最為平緩。

3)相同整數倍轉頻下,葉輪中部比葉片進口邊的幅值大。在0.8 mm間隙情況下,壓力脈動在經歷從極小值逐漸上升的過程中,出現壓力突然小幅下降后回升的現象。

4)蝸殼第三斷面和蝸殼隔舌附近的實測壓力波形走勢相反,葉頂間隙大小對兩個位置下的壓力系數波形圖影響不大。

[1] 李仁年,陳昊,韓偉,等. 大顆粒下螺旋離心泵內運動特性數值模擬與分析[J]. 流體機械,2016,44(7):18-21. Li Rennian,Chen Hao,Han Wei,et al. Numerical simulation and analysis of particle motion trajectory in the screw centrifugal pump with large particles[J]. Fluid Machinery,2016,44(7):18-21.(in Chinese with English abstract)

[2] 劉宜,白茂寧,程效銳,等. 螺旋離心泵的雙向流固耦合分析[J]. 西華大學學報:自然科學版,2015,34(6):85-91. Liu Yi,Bai Maoning,Cheng Xiaorui,et al. Bidirectional fluid structure interaction analysis for the screw centrifugal pump[J]. Journal of Xihua University:Natrual Science Edition,2015,34(6):85-91.(in Chinese with English abstract)

[3] 朱榮生,林鵬,龍云,等. 螺旋軸流泵的固液兩相流動數值模擬[J]. 排灌機械工程學報,2014,32(1):6-11. Zhu Rongsheng,Lin Peng,Long Yun,et al. Numerical simulation of solid-liquid two-phase flow in screw axial-flow pump[J]. Journal of Drainage and Irrigation Machinery Engineering,2014,32(1):6-11.(in Chinese with English abstract)

[4] 林鵬. 螺旋軸流泵的固液兩相數值模擬及試驗研究[D].鎮江:江蘇大學,2013. Lin Peng. Numerical Simulation and Experimental Research of Solid-Liquid Two-Phase in A Screw Axial Flow Pump[D]. Zhenjiang:Jiangsu University,2013.

[5] 張華. 單葉片螺旋離心式潛水排污泵的優化設計及試驗研究[D]. 鎮江:江蘇大學,2014. Zhang Hua. Optimal Design and Experimental Research on Screw-Centrifugal Submersible Pump with Single Blade [D]. Zhenjiang:Jiangsu University,2014.

[6] 梁開洪,張克危,許麗. 軸流泵葉頂間隙流動的計算流體動力學分析[J]. 華中科技大學學報:自然科學版,2004,32(9):36-38. Liang Kaihong,Zhang Kewei,Xu Li. Analysis of the flow through the blade tip clearances of axial pumps by CFD[J]. Journal of Huazhong University of Science and Technology:Nature Science Edition,2004,32(9):36-38.(in Chinese with English abstract)

[7] 戴辰辰,郭鵬程,羅興锜,等. 軸流泵端壁間隙流動特性的數值分析[J]. 流體機械,2009,37(6):32-35. Dai Chenchen,Guo Pengcheng,Luo Xingqi,et al. Numerical analysis of tip clearance flow characteristic in axial flow pump[J]. Fluid Mechinery,2009,37(6):32-35.(in Chinese with English abstract)

[8] 施衛東,張華,陳斌,等. 不同葉頂間隙下的軸流泵內部流場數值計算[J]. 排灌機械工程學報,2010,28(5):374-377. Shi Weidong,Zhang Hua,Chen Bin,et al. Numerical simulation of internal flow field in axial-flow pump with different blade tip clearance sizes[J]. Journal of Drainage and Irrigation Machinery Engineering,2010,28(5):374-377.(in Chinese with English abstract)

[9] 鄧育軒,李仁年,韓偉,等. 螺旋離心泵內回流渦空化特性[J]. 農業工程學報,2015,31(1):86-90. Deng Yuxuan,Li Rennian,Han Wei,et al. Characteristics of backflow vortex cavitation in screw centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2015,31(1):86-90.(in Chinese with English abstract)

[10] 鄧育軒,李仁年,韓偉,等. 螺旋離心泵內非定常漩渦空化的計算與分析[J].華中科技大學學報:自然科學版,2014,42(9):118-121. Deng Yuxuan,Li Rennian,Han Wei,et al. Unsteady vortex cavitation in a screw centrifugal pump[J]. Journal of Huazhong University of Science and Technology:Natrual Science Edition,2014,42(9):118-121.(in Chinese with English abstract)

[11] 權輝,李仁年,蘇清苗,等. 含沙水下粒徑對螺旋離心泵磨蝕效應的數值分析[J]. 西華大學學報:自然科學版,2014,33(3):91-94. Quan Hui,Li Rennian,Su Qingmiao,et al. Numerical analysis of erosion effect of sand-laden water on screw centrifugal pump[J]. Journal of Xihua University:Natrual Science Edition,2014,33(3):91-94.(in Chinese with English abstract)

[12] 張華,陳斌,施衛東,等. 單葉片螺旋離心泵內部流場數值計算及油膜試驗研究[J]. 排灌機械工程學報,2016,34(5):381-385. Zhang Hua,Chen Bin,Shi Weidong,et al. Numerical simulation and oil film observation on flow field in impeller of screw centrifugal pump with single-blade[J]. Journal of Drainage and Irrigation Machinery Engineering,2016,34(5):381-385.(in Chinese with English abstract)

[13] 權輝,李仁年,蘇清苗,等. 基于PIV測試的螺旋離心泵內部流動特性研究[J]. 農業機械學報,2015,46(4):28-32. Quan Hui,Li Rennian,Su Qingmiao,et al. Internal flow characteristic of screw centrifugal pump based on PIV[J]. Transactions of the Chinese Society for Agricultural Machinery,2015,46(4):28-32.(in Chinese with English abstract)

[14] Cheng X R,Li R N,Gao Y,et al. Numerical research on the effects of impeller pump-out vanes on axial force in a solid-liquid screw centrifugal pump[C]// 6th International Conference on Pumps and Fans with Compressors and Wind Turbines. IOP Conference Series:Materials Science &Engineering,2013,52:062008.

[15] Han Wei,Ma Wei,Li Rennian,et al. The numerical analysis of radial thrust and axial thrust in the screw centrifugal pump[J]. Procedia Engineering,2012(31):176-181.

[16] Quan H,Li R N,Su Q M,et al. Research on energy conversion mechanism of a screw centrifugal pump under the water[C]// 6th International Conference on Pumps and Fans with Compressors and Wind Turbines. IOP Conference Series:Materials Science &Engineering,2013,52:032017. [17] Scherbatenko I V,Surikov V I. Distribution of the absolute flow velocity behind the screw impeller of a centrifugal pump[J]. Chemical and Petroleum Engineering,2005,41(3):199-206.

[18] Pei J,Yuan S Q,Yuan J P. Numerical analysis of periodic flow unsteadiness in a single-blade centrifugal pump[J]. Science China-Technological Sciences,2013,56(1):212-221.

[19] Pei Ji,Yuan Shouqi,Yuan Jianping. Fluid-structure effects on periodically transient flow of a single-blade sewage centrifugal pump[J]. Journal of Mechanical Science and Technology,2013,27(7):2015-2023.

[20] 裴吉,袁壽其,袁建平,等. 單葉片離心泵壓力脈動強度多工況對比研究[J]. 華中科技大學學報:自然科學版,2013,41(12):29-33. Pei Ji,Yuan Shouqi,Yuan Jianping,et al. Comparative study of pressure fluctuation intensity for a single-blade pump under multiple operating conditions[J]. Journal of Huazhong University of Science and Technology:Natural Science Edition,2013,41(12):29-33.(in Chinese with English abstract)

[21] Nishi Y,Fukutomi J,Fujiwara R. Radial thrust of single-blade centrifugal pump[J]. International Journal of

10.11975/j.issn.1002-6819.2017.01.011

TH313

A

1002-6819(2017)-01-0084-06

張 華,陳 斌,王炳祺,史長彪,沈登武. 葉頂間隙大小對螺旋離心泵內部壓力脈動的影響[J]. 農業工程學報,2017,33(1):84-89.

10.11975/j.issn.1002-6819.2017.01.011 http://www.tcsae.org

Zhang Hua,Chen Bin,Wang Bingqi,Shi Changbiao,Shen Dengwu. Influence of tip clearance on internal pressure fluctuation of screw centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2017,33(1):84-89.(in Chinese with English abstract)doi:10.11975/j.issn.1002-6819.2017.01.011 http://www.tcsae.org

2016-06-29

2016-11-11

安徽省合肥工業大學“秋實計劃”A類計劃項目(JZ2015 QSJH0246);中央高校基本科研業務費專項資金資助(JZ2015HGBZ 0129);江蘇大學國家水泵及系統工程技術研究中心開放基金資助(NR CP-201603)。

張 華,男(漢族),江蘇南京人,助理研究員,博士,主要從事水力機械內部流動研究。合肥 合肥工業大學機械工程學院,230009。Email:zhh091088@126.com

※通信作者:陳 斌,男(漢族),江蘇靖江人,研究員,博士,博士生導師,主要從事流體機械及工程研究。合肥 合肥工業大學機械工程學院,230009。Email:chenbin21cn@126.com

主站蜘蛛池模板: 国产一二视频| 国产美女精品在线| 毛片一区二区在线看| 国产精品亚洲精品爽爽| 精品夜恋影院亚洲欧洲| 波多野结衣视频一区二区| 国产美女在线免费观看| 在线视频一区二区三区不卡| 婷婷色狠狠干| 天天综合网站| 99精品久久精品| 日韩在线播放中文字幕| 性色生活片在线观看| 自拍偷拍欧美| 99这里只有精品在线| 久久精品国产在热久久2019 | 色老头综合网| 亚洲一区二区日韩欧美gif| 婷婷在线网站| 久久久久亚洲精品成人网| www.youjizz.com久久| 伊人成人在线| 亚洲色无码专线精品观看| 在线一级毛片| 国产激情国语对白普通话| 99精品视频播放| 亚洲成在线观看| 成人午夜福利视频| 国产高清色视频免费看的网址| 伊人天堂网| 国产免费网址| 国产一级精品毛片基地| 国产午夜看片| 久久精品人妻中文系列| 日韩精品成人网页视频在线| 日本三级黄在线观看| 67194亚洲无码| 精品亚洲麻豆1区2区3区| 亚洲视频a| 伊人国产无码高清视频| 国产丝袜无码精品| 精品国产香蕉在线播出| 国产精品无码AV中文| 日韩在线欧美在线| 国产网站黄| 欧美在线一二区| 在线欧美日韩| 国产精品专区第一页在线观看| 国产永久免费视频m3u8| 国产精品久久久久久久久久久久| 国产成人AV男人的天堂| 国产嫩草在线观看| 婷婷色一二三区波多野衣| 一级高清毛片免费a级高清毛片| 久99久热只有精品国产15| 欧美国产综合色视频| 亚洲欧美在线综合图区| 91精品人妻互换| 中文字幕欧美日韩| 欧美色伊人| 成人福利在线免费观看| 超碰精品无码一区二区| 亚洲精品无码久久久久苍井空| 欧美啪啪网| 国产无码精品在线| 精品国产一区91在线| 亚洲成aⅴ人片在线影院八| 中文字幕第1页在线播| 久久精品午夜视频| 久久99国产乱子伦精品免| 国产精品亚洲一区二区三区z | 免费人欧美成又黄又爽的视频| 国产在线视频二区| 亚洲男人天堂2018| 久久免费观看视频| 无码一区二区波多野结衣播放搜索| 久久性妇女精品免费| 白浆免费视频国产精品视频 | 2020最新国产精品视频| 九九热视频在线免费观看| 99在线视频免费观看| 波多野结衣无码AV在线|