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

    功率超聲作用下鋼液中空化泡尺寸的演變特性

    2023-11-17 08:48:34程偉浩
    鋼鐵釩鈦 2023年5期
    關(guān)鍵詞:鋼液聲壓空化

    王 順,陳 敏,郭 慶,程偉浩

    (東北大學(xué)冶金學(xué)院,遼寧 沈陽 110819)

    0 引言

    功率超聲因其產(chǎn)生的空化效應(yīng)可提高化學(xué)反應(yīng)速率而廣泛應(yīng)用于冶金領(lǐng)域,Rayleigh 建立的空化泡動(dòng)力學(xué)模型為超聲空化的研究提供了理論依據(jù)[1-3]。目前,空化泡在功率超聲作用下的研究主要集中于水溶液模擬體系,以直觀的方式觀察空化泡運(yùn)動(dòng)過程中尺寸的變化,并不能定量描述空化泡尺寸的演變行為[4-5]。此外,目前關(guān)于空化泡的研究主要集中在常溫下水溶液中,或功率超聲作用下對(duì)鋼液中夾雜物去除的影響,缺少功率超聲對(duì)鋼液中空化泡尺寸演變的影響的系統(tǒng)研究[5-6]。

    由于空化泡不是規(guī)則的球形,且受力情況與能量的交換過程較為復(fù)雜,目前存在的方程難以準(zhǔn)確描述空化泡尺寸的變化規(guī)律[7-8]。基于這一問題,筆者利用Matlab 軟件運(yùn)用四階龍格庫(kù)塔算法[9],在RP 方程的基礎(chǔ)上引入能量粘滯損耗及輻射阻尼項(xiàng),對(duì)功率超聲作用下鋼液中單個(gè)空化泡尺寸的演變進(jìn)行研究,分析聲壓幅值、頻率、空化泡初始平衡半徑和氣體多變指數(shù)等因素對(duì)空化泡尺寸的影響,探究穩(wěn)態(tài)空化及瞬態(tài)空化條件下的空化效應(yīng),旨在得出功率超聲參數(shù)對(duì)鋼液中空化泡尺寸演變行為的影響規(guī)律,為探究功率超聲產(chǎn)生的空化效應(yīng)對(duì)鋼液中流動(dòng)及傳質(zhì)行為提供理論依據(jù)。

    1 模型建立

    1.1 基本假設(shè)

    功率超聲在鋼液中傳播產(chǎn)生的空化泡發(fā)生膨脹與收縮,為了簡(jiǎn)化計(jì)算做以下假設(shè)[10-11]:

    1)鋼液為不可壓縮流體,忽略溫度與密度的變化;

    2)空化泡為球形,且僅做徑向運(yùn)動(dòng),運(yùn)動(dòng)只與徑向軸有關(guān);

    3)空化泡的半徑遠(yuǎn)小于距鋼包壁的距離;

    4)不考慮重力和浮力對(duì)空化泡的影響;

    5)不考慮空化泡外液體邊界層的化學(xué)反應(yīng)。

    1.2 控制方程

    Rayleigh-Plesset 方程:

    功率超聲作用下鋼液與蒸氣處于相對(duì)平衡狀態(tài),則:

    由于空化泡演變過程中存在能量粘滯損耗及輻射阻尼,則:

    功率超聲作用在空化泡壁上的壓力:

    功率超聲周期T:

    式中,μ為粘滯系數(shù),Pa·s;R為空化泡瞬時(shí)半徑,μm;R0為空化泡初始平衡半徑,μm;νl為液相運(yùn)動(dòng)粘度,取值0.006 7 m2/s;σ為表面張力系數(shù),取值1.4 N/m;P0為靜壓力,取值101 325 Pa;PB(t)為空化泡內(nèi)壓力,取值6 Pa;P∞(t)為無窮遠(yuǎn)處壓力,Pa;PA為超聲波聲壓幅值,Pa;Pv為飽和蒸汽壓,Pa;ρl為液相密度,取值7 000 kg/m3;t為功率超聲作用時(shí)間,s;T為功率超聲作用周期,s;k 為氣體多變指數(shù);ω為功率超聲角頻率,rad/s;f為頻率,kHz。

    1.3 計(jì)算求解

    R-P 方程為非線性二階常微分方程,常規(guī)方法無法求解此方程的數(shù)值解,采用四階龍格庫(kù)塔法進(jìn)行迭代求解[4]。

    假設(shè)時(shí)間步長(zhǎng)為h,四階龍格庫(kù)塔法計(jì)算所涉及的方程如下。

    其中:

    1.4 研究方案

    通過改變聲壓幅值、頻率、空化泡初始平衡半徑及氣體多變指數(shù)考察功率超聲參數(shù)對(duì)鋼液中空化泡半徑的影響,以空化泡半徑的大小為指標(biāo)評(píng)價(jià)空化效應(yīng)的強(qiáng)弱。通過方程(7)與方程(8)可以得到空化泡半徑,以R/R0值表述空化泡尺寸的變化規(guī)律。

    所涉及的參數(shù)及初始條件取值為:聲壓幅值PA取值范圍是1~100P0;超聲頻率f的取值范圍是20~80 kHz;空化泡初始平衡半徑R0的取值范圍為1~50 μm;氣體多變指數(shù)k的取值范圍為1~1.65。其中,氣體多變指數(shù)定義為:在實(shí)際熱力內(nèi)平衡過程中,如果P乘以v的n次方為一個(gè)定值,則稱這個(gè)過程為多變過程,以k描述氣體多變指數(shù)。

    2 結(jié)果與討論

    2.1 聲壓幅值對(duì)鋼液中空化泡半徑演變過程的影響

    在頻率為24 kHz,空化泡初始平衡半徑為10 μm,氣體多變指數(shù)為1.35 條件下,探究聲壓幅值對(duì)空化泡半徑演變的影響。利用Matlab 軟件計(jì)算得到的空化泡半徑與初始平衡半徑的比值(R/R0)隨聲壓幅值變化的曲線如圖1 和圖2 所示。

    圖1 1~2P0 范圍內(nèi)聲壓幅值對(duì)鋼液中空化泡半徑的影響Fig.1 Influence of sound pressure amplitude on cavitation bubble radius in the liquid steel in the range of 1~2P0

    圖2 3~100P0 范圍內(nèi)聲壓幅值對(duì)鋼液中空化泡半徑的影響Fig.2 Influence of sound pressure amplitude on cavitation bubble radius in the liquid steel in the range of 3~100P0

    由圖1(a)可見,聲壓幅值在1~1.5P0范圍內(nèi)變化時(shí),鋼液中空化泡半徑變化趨勢(shì)為正弦曲線,隨著聲壓幅值的增加,空化泡的半徑變化幅度增強(qiáng)。當(dāng)聲壓幅值為2P0時(shí),在nT/2 周期內(nèi),空化泡半徑的演變僅表現(xiàn)為振動(dòng)幅度明顯增強(qiáng),但在(n+1)T/2 周期內(nèi),空化泡半徑演變的過程中出現(xiàn)振蕩現(xiàn)象,但空化泡未發(fā)生崩潰,為穩(wěn)態(tài)空化。聲壓幅值在1~2P0范圍內(nèi),空化泡半徑做周期性振蕩,由于聲壓幅值的增加,空化效應(yīng)對(duì)空化泡的作用增強(qiáng),表現(xiàn)為負(fù)壓區(qū)功率超聲對(duì)空化泡的拉伸作用增強(qiáng),空化泡的半徑增大;在正壓區(qū)功率超聲對(duì)空化泡的壓縮作用增強(qiáng),空化泡的半徑減小。空化泡的峰值半徑隨聲壓幅值的變化如圖1(b)所示,聲壓幅值由1P0增加至1.5P0時(shí),空化泡的峰值半徑由11.18 μm 增加至12.26 μm。當(dāng)聲壓幅值由1.5P0增加至2P0時(shí),空化泡的峰值半徑增加至17.03 μm,且斜率增加,表明空化泡的峰值半徑急劇增加,而空化泡半徑的劇烈變化易使空化現(xiàn)象由穩(wěn)態(tài)空化向瞬態(tài)空化轉(zhuǎn)變。因此有利于穩(wěn)態(tài)空化的聲壓幅值為1.5P0。另外,穩(wěn)態(tài)空化時(shí)產(chǎn)生的空化泡振蕩周期較長(zhǎng),持續(xù)上浮,可在鋼液鈣處理過程中有效促進(jìn)夾雜物上浮[10]。

    由圖2(a)可見,聲壓幅值在3~100P0范圍變化時(shí),隨著聲壓幅值的增加,空化泡半徑與初始平衡半徑的比值增加,空化效應(yīng)增強(qiáng)。空化泡振動(dòng)周期縮短,空化泡的半徑變化趨勢(shì)為非正弦曲線,空化泡經(jīng)歷若干次膨脹與收縮過程后發(fā)生崩潰,為瞬態(tài)空化。另外,隨著聲壓幅值的增大,空化泡經(jīng)歷的膨脹與收縮至崩潰次數(shù)逐漸減小,當(dāng)聲壓幅值為10P0時(shí),空化泡經(jīng)歷3 次膨脹與收縮后發(fā)生崩潰,而當(dāng)聲壓幅值大于30P0時(shí),空化泡在經(jīng)歷2 次膨脹與收縮過程后發(fā)生崩潰??栈荼罎⑶芭蛎浥c收縮過程的減少亦表明空化效應(yīng)的增強(qiáng)。由圖2(b)可見,聲壓幅值在3~30P0范圍內(nèi),空化泡的峰值半徑先增加后減小,在聲壓幅值由3P0增加至10P0時(shí),空化泡的峰值半徑由33.37 μm 增加至57.18 μm,當(dāng)聲壓幅值增加至30P0時(shí),空化泡的峰值半徑減小為34.86 μm。隨后,當(dāng)聲壓幅值由30P0增加至100P0時(shí),空化泡的峰值半徑增加至65.78 μm??栈莘逯蛋霃降脑黾佑欣诋a(chǎn)生瞬態(tài)空化,但過高的聲壓幅值脫離工業(yè)實(shí)際操作過程。且聲壓幅值為10P0與80P0條件下所對(duì)應(yīng)的空化泡峰值半徑分別為57.18 μm 與58.02 μm,這表明10P0與80P0條件下所產(chǎn)生的空化效應(yīng)強(qiáng)度基本一致,考慮工業(yè)實(shí)際應(yīng)用,瞬態(tài)空化適宜的聲壓幅值為10P0。另外,由于瞬態(tài)空化條件下空化泡產(chǎn)生至崩潰所持續(xù)的時(shí)間較短,且較高的聲壓幅值可加速空化泡的運(yùn)動(dòng),可有效提高鋼液內(nèi)的傳質(zhì)速率。

    2.2 穩(wěn)態(tài)空化(1.5P0)條件下各參數(shù)對(duì)鋼液中空化泡半徑的影響

    2.2.1 初始平衡半徑對(duì)鋼液中空化泡半徑的影響

    由圖3(a)可見,空化泡初始平衡半徑在5~50 μm的范圍內(nèi)變化時(shí),隨著空化泡初始平衡半徑增大,空化泡膨脹與收縮的程度增強(qiáng)。這是由于初始平衡半徑增大,可提高空化泡與鋼液的接觸面積,促進(jìn)空化泡對(duì)周圍能量的吸收,導(dǎo)致空化泡的峰值半徑增大。當(dāng)空化泡初始平衡半徑為50 μm 時(shí),空化泡半徑變化周期延長(zhǎng),這是由于空化泡初始平衡半徑較大時(shí),需要更充足的時(shí)間吸收能量,導(dǎo)致空化泡膨脹與收縮周期增加。由圖3(b)可見,隨著空化泡半徑的增加,空化泡的峰值半徑增加,當(dāng)空化泡初始平衡半徑由5 μm 增加至50 μm 時(shí),空化泡的峰值半徑由5.35 μm 增加至141.75 μm,與空化泡初始平衡半徑相比,空化泡的峰值半徑分別增加了0.35 μm 和91.75 μm。因此,為提高穩(wěn)態(tài)空化條件下的空化效果,空化泡初始平衡半徑應(yīng)控制在50 μm。

    圖3 5~50 μm 范圍內(nèi)初始平衡半徑對(duì)鋼液中空化泡半徑的影響Fig.3 Influence of initial equilibrium radius on cavitation bubble radius in the molten steel in the range of 5~50 μm

    2.2.2 頻率對(duì)鋼液中空化泡半徑的影響

    由圖4(a)可見,頻率在20~32 kHz 范圍內(nèi)變化時(shí),空化泡半徑變化趨勢(shì)基本一致,此頻率范圍內(nèi)對(duì)空化泡半徑的影響較小。當(dāng)頻率為50 kHz 和80 kHz 時(shí),空化泡半徑的振蕩幅度提高,表明較高的頻率可提高空化泡的半徑,且在一個(gè)周期內(nèi)空化泡半徑出現(xiàn)多次振蕩。圖4(b)可見,隨著功率超聲頻率的增加,空化泡的峰值半徑增大。當(dāng)功率超聲的頻率由20 kHz 增加至80 kHz 時(shí),空化泡的峰值半徑由11.21 μm 增加至12.29 μm,與功率超聲頻率增加60 kHz 的幅度相比,空化泡的峰值半徑僅增加了1.08 μm。說明功率超聲的頻率對(duì)空化泡半徑演變的影響較小。

    圖4 20~80 kHz 范圍內(nèi)頻率對(duì)空化泡半徑的影響Fig.4 Influence of frequency on cavitation bubble radius in the range of 20~80 kHz

    2.2.3 氣體多變指數(shù)對(duì)鋼液中空化泡半徑的影響

    由圖5(a)可見,隨著氣體多變指數(shù)k的增加,空化泡膨脹與收縮的周期基本保持一致,空化泡半徑的振蕩幅度逐漸降低,且同一周期內(nèi)的振蕩次數(shù)逐漸減少,說明氣體多變指數(shù)的增加對(duì)空化泡半徑產(chǎn)生負(fù)面影響。由圖5(b)可見,隨著k值的增加,空化泡的峰值半徑逐漸減小。當(dāng)氣體多變指數(shù)由1增加至1.65 時(shí),空化泡的峰值半徑由12.68 μm 減小至11.08 μm,表明增加氣體多變指數(shù)將削弱功率超聲的空化強(qiáng)度。主要原因在于單原子氣體空化泡的峰值半徑大于多原子氣體,且振蕩周期基本一致。上述結(jié)論與多變指數(shù)較小的單原子氣體的空化效應(yīng)優(yōu)于多原子氣體的結(jié)論相吻合[6]。

    圖5 氣體多變指數(shù)對(duì)空化泡半徑的影響Fig.5 Influence of gas multiplicity index on cavitation bubble radius

    2.3 瞬態(tài)空化(10P0)條件下各參數(shù)對(duì)鋼液中空化泡半徑的影響

    2.3.1 初始平衡半徑對(duì)鋼液中空化泡半徑的影響

    由圖6(a)可見,與穩(wěn)態(tài)空化條件下空化泡半徑的變化曲線相比,空化泡半徑變化幅度明顯提高。當(dāng)空化泡初始平衡半徑為5 μm 時(shí),空化泡僅經(jīng)歷2 次膨脹與收縮的過程后發(fā)生崩潰。而空化泡初始平衡半徑為10、15 和20 μm 時(shí),空化泡膨脹與收縮次數(shù)增加,分別為3、5 和8 次膨脹與收縮過程后發(fā)生崩潰。當(dāng)空化泡初始平衡半徑為50 μm 時(shí),空化泡經(jīng)歷13 次膨脹與收縮過程后發(fā)生崩潰,且空化泡膨脹幅度明顯減小,表明隨著空化泡初始平衡半徑的增加,空化效應(yīng)減弱。另外,空化泡初始平衡半徑為5 μm 時(shí),空化泡發(fā)生崩潰時(shí)峰值半徑達(dá)到423.01 μm,空化泡最大半徑膨脹率為846%,而崩潰時(shí)膨脹率越大,釋放的能量越多,空化效應(yīng)增強(qiáng)。由圖6(b)可見,空化泡的峰值半徑隨著初始平衡半徑的增加先增大后減小。當(dāng)空化泡初始平衡半徑由5 μm增加至20 μm 時(shí),空化泡的峰值半徑由423.01 μm增加至896.12 μm,此時(shí)的膨脹率為448%;當(dāng)初始平衡半徑增加至50 μm 時(shí),空化泡的峰值半徑減小到544.16 μm,此時(shí)空化泡膨脹率為108.8%。半徑由5 μm 增大至20 μm 過程中峰值半徑不斷增加,但發(fā)生崩潰時(shí)的膨脹率由846%降低為448%,瞬態(tài)空化效應(yīng)減弱。因此,空化泡初始平衡半徑為5 μm 時(shí)有利于發(fā)生瞬態(tài)空化。

    圖6 5~50 μm 范圍內(nèi)初始平衡半徑對(duì)鋼液中空化泡半徑的影響Fig.6 Influence of initial equilibrium radius on cavitation bubble radius in the molten steel in the range of 5~50 μm

    2.3.2 頻率對(duì)鋼液中空化泡半徑的影響

    由圖7(a)可見,頻率在20~28 kHz 范圍內(nèi)變化時(shí),空化泡經(jīng)歷2 次膨脹與收縮的過程后發(fā)生崩潰,且發(fā)生崩潰時(shí)空化泡半徑的膨脹率逐漸減小。頻率為32、50 和80 kHz 時(shí),空化泡分別經(jīng)歷4、5和8 次膨脹與收縮過程后崩潰。瞬態(tài)空化過程中,發(fā)生崩潰所經(jīng)歷的膨脹與收縮次數(shù)越多,其空化強(qiáng)度越弱。由圖7(b)可見,在低頻率20~28 kHz 范圍內(nèi),隨著頻率的增加,空化泡的峰值半徑由488.05 μm 減小至354.43 μm;當(dāng)頻率為32 kHz 時(shí),空化泡的峰值半徑為593.78 μm,為各頻率條件下半徑最大值,但其膨脹與收縮次數(shù)較多,空化泡崩潰時(shí)間較長(zhǎng),空化效應(yīng)減弱。當(dāng)頻率為80 kHz 時(shí),空化泡的峰值半徑減小至172.65 μm。表明高頻率的條件下空化泡膨脹與收縮時(shí)間延長(zhǎng),空化核未能及時(shí)發(fā)生崩潰,空化效應(yīng)減弱。因此,為提高超聲的瞬態(tài)空化強(qiáng)度,功率超聲的頻率應(yīng)選擇為20 kHz。

    圖7 20~80 kHz 范圍內(nèi)頻率對(duì)空化泡半徑的影響Fig.7 Influence of frequency on cavitation bubble radius in the range of 20~80 kHz

    2.3.3 氣體多變指數(shù)對(duì)鋼液中空化泡半徑的影響

    圖8(a)可見,瞬態(tài)空化條件下,空化泡僅經(jīng)歷1 次膨脹與收縮過程即發(fā)生崩潰,且隨著氣體多變指數(shù)k的增加,空化泡崩潰所經(jīng)歷的周期基本一致,表明瞬態(tài)空化條件下氣體多變指數(shù)對(duì)空化泡膨脹與收縮過程影響較小。主要原因在于瞬態(tài)空化頻率高、能量大,空化泡崩潰過程中單原子及多原子氣體的作用性能減弱,對(duì)空化效應(yīng)強(qiáng)弱的影響并不明顯。由圖8(b)可見,隨著氣體多變指數(shù)的增加,空化泡的峰值半徑逐漸減小。當(dāng)氣體多變指數(shù)由1 增加至1.65 時(shí),空化泡的峰值半徑由18.22 μm 減小至17.87 μm,僅減小了0.35 μm,亦表明氣體多變指數(shù)對(duì)瞬態(tài)空化強(qiáng)弱的影響較小。但考慮能量的吸收與釋放過程,多變指數(shù)較低的單原子氣體有利于產(chǎn)生瞬態(tài)空化。

    圖8 氣體多變指數(shù)對(duì)空化泡半徑的影響Fig.8 Influence of gas multiplicity index on cavitation bubble radius

    3 結(jié)論與展望

    1)功率超聲作用下引起鋼液中空化泡的膨脹與收縮,不同聲壓幅值可產(chǎn)生不同的空化效應(yīng)。聲壓幅值在1~2P0的范圍內(nèi)產(chǎn)生穩(wěn)態(tài)空化,可促進(jìn)夾雜物上浮過程;聲壓幅值在3~100P0的范圍內(nèi)產(chǎn)生瞬態(tài)空化,可加速鋼液內(nèi)的傳質(zhì)過程。

    2)穩(wěn)態(tài)空化條件下,隨著頻率和初始平衡半徑的增加,空化泡的峰值半徑增大。初始平衡半徑由5 μm 增加至50 μm 時(shí),空化泡的峰值半徑由5.35 μm增加至141.75 μm;頻率由20 kHz 增加至80 kHz 時(shí),空化泡的峰值半徑由11.45 μm 增加至12.29 μm。

    3)瞬態(tài)空化條件下,頻率為20 kHz 時(shí),空化泡經(jīng)歷2 次膨脹與收縮過程后崩潰,對(duì)應(yīng)空化泡的峰值半徑為488.05 μm。初始平衡半徑對(duì)空化泡尺寸的影響顯著,由5 μm 增加至20 μm 時(shí),空化泡的峰值半徑由423.01 μm 增加至896.12 μm,當(dāng)初始平衡半徑增加至50 μm 時(shí),空化泡的峰值半徑減小至544.16 μm。

    4)氣體多變指數(shù)對(duì)穩(wěn)態(tài)和瞬態(tài)空化條件下空化泡尺寸的影響較小,但氣體多變指數(shù)較小的單原子氣體產(chǎn)生的空化效應(yīng)優(yōu)于多原子氣體。當(dāng)氣體多變指數(shù)由1 增加至1.65 時(shí),空化泡的峰值半徑僅由18.22 μm 減小至17.87 μm。

    5)本研究成果可為探究功率超聲產(chǎn)生的空化效應(yīng)對(duì)鋼液中流動(dòng)及傳質(zhì)行為提供理論依據(jù),但空化泡的體積、密度及空化泡之間的作用對(duì)空化泡膨脹與收縮產(chǎn)生影響,后續(xù)需進(jìn)一步修正空化泡的尺寸變化方程,并結(jié)合實(shí)際設(shè)計(jì)試驗(yàn),對(duì)模擬結(jié)果進(jìn)行驗(yàn)證。

    猜你喜歡
    鋼液聲壓空化
    基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
    基于格子Boltzmann方法的雙空化泡遠(yuǎn)壁區(qū)潰滅規(guī)律研究
    真空熔煉過程中定向?qū)Я餮b置內(nèi)鋼液運(yùn)動(dòng)現(xiàn)象的研究
    日鋼板坯連鑄SPA-H鋼液位波動(dòng)原因分析及控制
    山東冶金(2019年6期)2020-01-06 07:45:56
    LF深脫硫過程中的鋼液增碳增氮的控制研究
    車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
    汽車工程(2018年12期)2019-01-29 06:46:36
    三維扭曲水翼空化現(xiàn)象CFD模擬
    對(duì)煉鋼生產(chǎn)中轉(zhuǎn)爐煉鋼脫氧工藝的探討
    科技尚品(2016年1期)2016-05-30 06:01:59
    不同運(yùn)動(dòng)形式下水物相互作用空化數(shù)值模擬
    基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
    黑人操中国人逼视频| 欧美xxxx黑人xx丫x性爽| 丝袜人妻中文字幕| 97超级碰碰碰精品色视频在线观看| 法律面前人人平等表现在哪些方面| 中文字幕精品亚洲无线码一区| 国产日本99.免费观看| 亚洲欧美精品综合久久99| 国产精品电影一区二区三区| 亚洲最大成人中文| 99久久久亚洲精品蜜臀av| av国产免费在线观看| 99精品欧美一区二区三区四区| 色综合站精品国产| 又大又爽又粗| 色综合欧美亚洲国产小说| 精品国产乱子伦一区二区三区| 国产精品99久久久久久久久| 18禁观看日本| 少妇人妻一区二区三区视频| 51午夜福利影视在线观看| 欧美色视频一区免费| 狂野欧美白嫩少妇大欣赏| 午夜福利视频1000在线观看| 老汉色∧v一级毛片| 国产高清视频在线播放一区| 变态另类成人亚洲欧美熟女| 久久久久久久久中文| 真人做人爱边吃奶动态| 在线国产一区二区在线| 国产亚洲精品综合一区在线观看| 欧美日本视频| 91老司机精品| 中文字幕最新亚洲高清| 国产日本99.免费观看| 日本黄色片子视频| 国产69精品久久久久777片 | av女优亚洲男人天堂 | 国模一区二区三区四区视频 | 国产单亲对白刺激| 国产欧美日韩一区二区精品| 午夜成年电影在线免费观看| 国产成人精品久久二区二区免费| 精品不卡国产一区二区三区| 一a级毛片在线观看| 亚洲狠狠婷婷综合久久图片| 99在线视频只有这里精品首页| 色av中文字幕| 搡老熟女国产l中国老女人| 国产亚洲精品久久久久久毛片| 日韩欧美精品v在线| 男女床上黄色一级片免费看| 这个男人来自地球电影免费观看| 免费看日本二区| 可以在线观看毛片的网站| 久久香蕉精品热| 黑人操中国人逼视频| 天堂√8在线中文| 亚洲成人精品中文字幕电影| 亚洲真实伦在线观看| 国产精品综合久久久久久久免费| 精品不卡国产一区二区三区| 丝袜人妻中文字幕| 精品欧美国产一区二区三| 亚洲av电影不卡..在线观看| 日韩国内少妇激情av| 岛国在线观看网站| 美女午夜性视频免费| 国产一区二区激情短视频| 男人的好看免费观看在线视频| 婷婷丁香在线五月| 国产成+人综合+亚洲专区| 国产私拍福利视频在线观看| 一边摸一边抽搐一进一小说| 在线播放国产精品三级| 国产精品99久久99久久久不卡| 国内少妇人妻偷人精品xxx网站 | 狂野欧美激情性xxxx| 波多野结衣巨乳人妻| 小说图片视频综合网站| 日韩 欧美 亚洲 中文字幕| 19禁男女啪啪无遮挡网站| 这个男人来自地球电影免费观看| 99久久成人亚洲精品观看| 亚洲专区字幕在线| 五月伊人婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 99热6这里只有精品| 亚洲av电影不卡..在线观看| 亚洲国产色片| 女同久久另类99精品国产91| 高潮久久久久久久久久久不卡| 在线免费观看不下载黄p国产 | 欧美成狂野欧美在线观看| 久久久国产成人免费| 男女视频在线观看网站免费| 午夜亚洲福利在线播放| 日本免费a在线| 亚洲狠狠婷婷综合久久图片| 国产黄色小视频在线观看| 久久精品91无色码中文字幕| 又爽又黄无遮挡网站| 一个人看视频在线观看www免费 | 美女高潮喷水抽搐中文字幕| 久久性视频一级片| 日本一本二区三区精品| 欧美日韩黄片免| 天堂√8在线中文| 国产亚洲av嫩草精品影院| 亚洲国产精品久久男人天堂| 曰老女人黄片| 欧美在线黄色| 一边摸一边抽搐一进一小说| 午夜久久久久精精品| 老熟妇乱子伦视频在线观看| 国产成人系列免费观看| а√天堂www在线а√下载| 免费看日本二区| 亚洲无线观看免费| 在线视频色国产色| or卡值多少钱| 欧美色欧美亚洲另类二区| 天堂av国产一区二区熟女人妻| 国产极品精品免费视频能看的| 久久久成人免费电影| 亚洲av熟女| 又爽又黄无遮挡网站| 精品久久久久久久人妻蜜臀av| 男女午夜视频在线观看| 国产成人精品久久二区二区91| 国产成人欧美在线观看| 欧美日韩瑟瑟在线播放| 美女午夜性视频免费| 91在线精品国自产拍蜜月 | 男女之事视频高清在线观看| aaaaa片日本免费| 亚洲国产中文字幕在线视频| 亚洲最大成人中文| 国产亚洲av高清不卡| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久久免费视频了| 真实男女啪啪啪动态图| 级片在线观看| 国产高清三级在线| 午夜久久久久精精品| xxx96com| 亚洲专区字幕在线| 成人高潮视频无遮挡免费网站| 免费在线观看视频国产中文字幕亚洲| 97碰自拍视频| 综合色av麻豆| 欧美丝袜亚洲另类 | 全区人妻精品视频| АⅤ资源中文在线天堂| 91麻豆av在线| 一区二区三区激情视频| 床上黄色一级片| 久久午夜亚洲精品久久| 亚洲av中文字字幕乱码综合| 亚洲片人在线观看| 一个人观看的视频www高清免费观看 | 国产麻豆成人av免费视频| 日日夜夜操网爽| 老熟妇仑乱视频hdxx| 国产高清视频在线观看网站| 免费观看的影片在线观看| 久久午夜综合久久蜜桃| 午夜视频精品福利| 国产一区二区三区在线臀色熟女| 亚洲av电影不卡..在线观看| 熟女电影av网| 国产精品日韩av在线免费观看| 国产麻豆成人av免费视频| 人人妻人人看人人澡| 亚洲国产日韩欧美精品在线观看 | www.999成人在线观看| 51午夜福利影视在线观看| 两性夫妻黄色片| 久久久久久久久免费视频了| 亚洲天堂国产精品一区在线| 久久这里只有精品19| 一卡2卡三卡四卡精品乱码亚洲| 亚洲 欧美一区二区三区| 成人高潮视频无遮挡免费网站| 日韩欧美 国产精品| 他把我摸到了高潮在线观看| 国产精品精品国产色婷婷| 久久精品91无色码中文字幕| 老司机午夜十八禁免费视频| 精品一区二区三区av网在线观看| 国产1区2区3区精品| 99久久精品一区二区三区| 99久久精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 黄色丝袜av网址大全| 香蕉久久夜色| 身体一侧抽搐| av欧美777| 90打野战视频偷拍视频| 日本 欧美在线| 美女高潮的动态| 国产精品电影一区二区三区| 亚洲自偷自拍图片 自拍| 日本成人三级电影网站| 97碰自拍视频| 他把我摸到了高潮在线观看| 国产欧美日韩精品亚洲av| 九色成人免费人妻av| 亚洲成人久久爱视频| 久久精品国产综合久久久| 一级a爱片免费观看的视频| 免费观看人在逋| 99久久精品国产亚洲精品| 无限看片的www在线观看| 久久久国产欧美日韩av| 欧美xxxx黑人xx丫x性爽| 中文字幕av在线有码专区| 国产野战对白在线观看| 日本a在线网址| 亚洲精品乱码久久久v下载方式 | 亚洲无线观看免费| 男女下面进入的视频免费午夜| 色播亚洲综合网| 午夜久久久久精精品| 一区二区三区高清视频在线| 亚洲国产欧美人成| 亚洲人与动物交配视频| 欧美日韩瑟瑟在线播放| 中文字幕人妻丝袜一区二区| 国产视频内射| 网址你懂的国产日韩在线| 亚洲片人在线观看| 国产成人aa在线观看| 在线观看免费视频日本深夜| 色吧在线观看| 欧美中文综合在线视频| 国产精品综合久久久久久久免费| 99riav亚洲国产免费| 啦啦啦观看免费观看视频高清| 亚洲国产精品999在线| 久久久久九九精品影院| 黑人操中国人逼视频| 色噜噜av男人的天堂激情| 一边摸一边抽搐一进一小说| 亚洲av成人一区二区三| 久久精品影院6| 精品久久久久久成人av| 欧美av亚洲av综合av国产av| 无人区码免费观看不卡| 亚洲成人久久性| 欧美在线一区亚洲| 亚洲无线观看免费| 老汉色av国产亚洲站长工具| 亚洲国产中文字幕在线视频| 色在线成人网| 无人区码免费观看不卡| 午夜视频精品福利| 中国美女看黄片| 老司机深夜福利视频在线观看| 十八禁网站免费在线| 亚洲精品乱码久久久v下载方式 | 色综合亚洲欧美另类图片| 99国产精品一区二区蜜桃av| 欧美丝袜亚洲另类 | 亚洲五月天丁香| 法律面前人人平等表现在哪些方面| 欧美日本视频| 亚洲 国产 在线| 哪里可以看免费的av片| 狠狠狠狠99中文字幕| 熟女人妻精品中文字幕| 国内精品一区二区在线观看| 香蕉丝袜av| 国产又黄又爽又无遮挡在线| 男女午夜视频在线观看| 国产亚洲精品久久久com| 91在线观看av| 色吧在线观看| 国产成人欧美在线观看| 亚洲熟女毛片儿| 嫩草影院精品99| 色综合欧美亚洲国产小说| 国产一区二区三区视频了| 国产精品自产拍在线观看55亚洲| 欧美大码av| 国产一区在线观看成人免费| 久久久色成人| 手机成人av网站| АⅤ资源中文在线天堂| 成人三级做爰电影| 精品福利观看| 亚洲国产精品999在线| 麻豆国产97在线/欧美| 免费搜索国产男女视频| 国内精品久久久久久久电影| 91老司机精品| 岛国在线观看网站| 嫩草影视91久久| 亚洲av美国av| 18美女黄网站色大片免费观看| 精品日产1卡2卡| 超碰成人久久| 成年女人永久免费观看视频| 国产高清videossex| 久久精品人妻少妇| 亚洲黑人精品在线| 香蕉久久夜色| 88av欧美| 老司机福利观看| 亚洲五月天丁香| 国产精品亚洲美女久久久| 日韩三级视频一区二区三区| 国产又黄又爽又无遮挡在线| 国产一区在线观看成人免费| 啪啪无遮挡十八禁网站| 亚洲黑人精品在线| 欧美极品一区二区三区四区| 日韩人妻高清精品专区| 18禁裸乳无遮挡免费网站照片| 99久久99久久久精品蜜桃| 国产91精品成人一区二区三区| 很黄的视频免费| 欧美日韩精品网址| 亚洲精品中文字幕一二三四区| 国产在线精品亚洲第一网站| 97人妻精品一区二区三区麻豆| a级毛片a级免费在线| 欧美一级毛片孕妇| 日本黄大片高清| 97超视频在线观看视频| 成人鲁丝片一二三区免费| 亚洲欧美日韩东京热| 国产精品,欧美在线| 色综合欧美亚洲国产小说| 亚洲成人久久性| 母亲3免费完整高清在线观看| 黑人操中国人逼视频| 日韩欧美三级三区| 国产成人欧美在线观看| 欧美日韩黄片免| 香蕉丝袜av| 亚洲 欧美一区二区三区| 国产一区在线观看成人免费| 免费在线观看亚洲国产| 亚洲精品在线美女| 色播亚洲综合网| 丁香六月欧美| 免费电影在线观看免费观看| 国产v大片淫在线免费观看| 小说图片视频综合网站| 国产精品久久久久久精品电影| 一边摸一边抽搐一进一小说| 亚洲第一电影网av| 亚洲最大成人中文| 色噜噜av男人的天堂激情| 欧美黑人巨大hd| 久久热在线av| 老司机在亚洲福利影院| 精华霜和精华液先用哪个| 九色成人免费人妻av| 日韩有码中文字幕| 免费一级毛片在线播放高清视频| 又粗又爽又猛毛片免费看| 19禁男女啪啪无遮挡网站| 日韩高清综合在线| 色尼玛亚洲综合影院| 又黄又爽又免费观看的视频| 国产不卡一卡二| 一个人观看的视频www高清免费观看 | 国产午夜精品论理片| 天天躁狠狠躁夜夜躁狠狠躁| 淫秽高清视频在线观看| 一进一出抽搐动态| 一二三四在线观看免费中文在| 十八禁网站免费在线| 亚洲国产精品久久男人天堂| 啦啦啦免费观看视频1| 久久这里只有精品19| 成人永久免费在线观看视频| 国产男靠女视频免费网站| 国产黄色小视频在线观看| 日韩欧美在线二视频| 精品国产亚洲在线| 国产av不卡久久| 成人鲁丝片一二三区免费| 日日夜夜操网爽| 两人在一起打扑克的视频| 色播亚洲综合网| 欧美国产日韩亚洲一区| 久久精品国产综合久久久| 色在线成人网| 母亲3免费完整高清在线观看| 麻豆国产av国片精品| 香蕉久久夜色| 欧美日韩综合久久久久久 | 老熟妇仑乱视频hdxx| 成年女人永久免费观看视频| 香蕉国产在线看| 亚洲人成网站高清观看| 亚洲精品在线美女| 岛国视频午夜一区免费看| a级毛片a级免费在线| 在线观看午夜福利视频| 夜夜夜夜夜久久久久| 国产免费男女视频| 动漫黄色视频在线观看| 欧美黑人欧美精品刺激| 国产精品爽爽va在线观看网站| 亚洲性夜色夜夜综合| 中亚洲国语对白在线视频| 精品久久蜜臀av无| 欧美最黄视频在线播放免费| 国产一区二区三区视频了| 九色成人免费人妻av| 欧美中文日本在线观看视频| 黑人巨大精品欧美一区二区mp4| 97超视频在线观看视频| 一本一本综合久久| 窝窝影院91人妻| www日本黄色视频网| av片东京热男人的天堂| 成人国产一区最新在线观看| 成年人黄色毛片网站| 一区二区三区激情视频| 久9热在线精品视频| 午夜福利高清视频| 亚洲av熟女| 99久久成人亚洲精品观看| 九九久久精品国产亚洲av麻豆 | 成人av一区二区三区在线看| 亚洲成a人片在线一区二区| 巨乳人妻的诱惑在线观看| 99热这里只有精品一区 | 中出人妻视频一区二区| 久久精品aⅴ一区二区三区四区| 国产综合懂色| 国产成人福利小说| 欧美乱码精品一区二区三区| 日本免费a在线| www.熟女人妻精品国产| 五月玫瑰六月丁香| 99国产精品99久久久久| 日本免费a在线| 婷婷精品国产亚洲av| 嫁个100分男人电影在线观看| 国产高清视频在线播放一区| 美女高潮的动态| avwww免费| 真实男女啪啪啪动态图| 亚洲国产欧美人成| 国产免费男女视频| 国产伦人伦偷精品视频| 国产视频一区二区在线看| 精品欧美国产一区二区三| 不卡一级毛片| 99riav亚洲国产免费| 网址你懂的国产日韩在线| а√天堂www在线а√下载| 波多野结衣巨乳人妻| 变态另类丝袜制服| 国产精品99久久久久久久久| 嫩草影视91久久| a在线观看视频网站| 成人国产综合亚洲| 精品国产乱码久久久久久男人| 国产在线精品亚洲第一网站| 久久精品aⅴ一区二区三区四区| 久久精品国产99精品国产亚洲性色| 午夜福利在线在线| 在线观看日韩欧美| 不卡一级毛片| 午夜精品在线福利| 亚洲精品在线美女| 欧美乱妇无乱码| 日韩中文字幕欧美一区二区| 这个男人来自地球电影免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产高清有码在线观看视频| 欧美乱色亚洲激情| 天天添夜夜摸| 欧美一级a爱片免费观看看| 国产精品免费一区二区三区在线| 日本三级黄在线观看| 精品国产三级普通话版| 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦免费观看视频1| 狂野欧美激情性xxxx| 免费看a级黄色片| 在线观看日韩欧美| 成人av一区二区三区在线看| 成年人黄色毛片网站| 亚洲专区国产一区二区| 亚洲国产欧美一区二区综合| 亚洲午夜精品一区,二区,三区| 午夜免费激情av| 波多野结衣高清无吗| 蜜桃久久精品国产亚洲av| 亚洲国产欧洲综合997久久,| 亚洲一区二区三区不卡视频| 日韩高清综合在线| 婷婷精品国产亚洲av在线| 亚洲avbb在线观看| 国产午夜精品久久久久久| 日本熟妇午夜| 后天国语完整版免费观看| 日韩欧美精品v在线| 久久午夜综合久久蜜桃| av在线天堂中文字幕| 在线a可以看的网站| 在线观看日韩欧美| 最新中文字幕久久久久 | 人人妻人人澡欧美一区二区| 色尼玛亚洲综合影院| 欧美一区二区国产精品久久精品| 亚洲av五月六月丁香网| 后天国语完整版免费观看| 日韩成人在线观看一区二区三区| 久久国产精品人妻蜜桃| 国产精品av久久久久免费| 国产精品1区2区在线观看.| 精品一区二区三区视频在线 | 性色av乱码一区二区三区2| 亚洲电影在线观看av| 欧美丝袜亚洲另类 | 人人妻人人看人人澡| 国产精品久久久久久人妻精品电影| 不卡一级毛片| 婷婷丁香在线五月| 色综合亚洲欧美另类图片| 男女那种视频在线观看| 亚洲国产欧洲综合997久久,| 午夜成年电影在线免费观看| 在线观看舔阴道视频| 日本a在线网址| 国产精品1区2区在线观看.| 观看美女的网站| 老司机深夜福利视频在线观看| 99国产精品一区二区蜜桃av| 麻豆成人av在线观看| 男人的好看免费观看在线视频| 日韩中文字幕欧美一区二区| 一边摸一边抽搐一进一小说| 亚洲精品在线美女| 欧美成人一区二区免费高清观看 | 69av精品久久久久久| 午夜两性在线视频| 国产淫片久久久久久久久 | 国产97色在线日韩免费| 中文字幕人成人乱码亚洲影| 99热这里只有精品一区 | 黄片大片在线免费观看| 亚洲国产高清在线一区二区三| 久久久久国内视频| 日韩中文字幕欧美一区二区| 99在线视频只有这里精品首页| 男人舔女人的私密视频| 国产亚洲欧美98| 巨乳人妻的诱惑在线观看| 又大又爽又粗| 夜夜看夜夜爽夜夜摸| 国产69精品久久久久777片 | 搡老熟女国产l中国老女人| 亚洲五月天丁香| 婷婷精品国产亚洲av在线| 观看免费一级毛片| 亚洲av电影在线进入| 麻豆成人av在线观看| 日本黄色片子视频| 亚洲av成人av| 国产精品一区二区精品视频观看| 亚洲人成伊人成综合网2020| 91av网一区二区| 丰满人妻熟妇乱又伦精品不卡| 黑人操中国人逼视频| 国产一区二区激情短视频| 欧美日本视频| 久久亚洲真实| 亚洲精品在线美女| 无人区码免费观看不卡| 亚洲成人中文字幕在线播放| 老司机午夜十八禁免费视频| 久久久国产欧美日韩av| 亚洲av电影不卡..在线观看| 久久精品aⅴ一区二区三区四区| 最新美女视频免费是黄的| 亚洲精品国产精品久久久不卡| 激情在线观看视频在线高清| 很黄的视频免费| 99热只有精品国产| 日韩精品中文字幕看吧| 一进一出抽搐动态| 丰满人妻一区二区三区视频av | netflix在线观看网站| 午夜福利视频1000在线观看| 五月伊人婷婷丁香| 在线国产一区二区在线| 毛片女人毛片| 午夜福利欧美成人| 波多野结衣巨乳人妻| 日本一二三区视频观看| 国产又色又爽无遮挡免费看| 在线国产一区二区在线| 久久午夜亚洲精品久久| 婷婷六月久久综合丁香| 欧美日本视频| 国产淫片久久久久久久久 | 久久精品91蜜桃| 美女扒开内裤让男人捅视频| 国产精品自产拍在线观看55亚洲| 三级国产精品欧美在线观看 | av黄色大香蕉| 国内毛片毛片毛片毛片毛片| 一进一出抽搐gif免费好疼| 日韩欧美国产在线观看| 欧美日韩乱码在线|