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

    基于吸脫附理論的RDX驅(qū)酸洗滌動力學(xué)模型研究

    2023-11-16 01:28:14陳厚和
    火炸藥學(xué)報 2023年10期
    關(guān)鍵詞:磁力云圖超聲波

    謝 強,陳厚和,李 慧

    (1.安徽理工大學(xué) 化學(xué)工程學(xué)院,安徽 淮南 232001;2.南京理工大學(xué) 化學(xué)與化工學(xué)院,江蘇 南京 210094)

    引 言

    黑索今(RDX)是一種應(yīng)用廣泛的單質(zhì)高能炸藥[1-3],其生產(chǎn)能力與本質(zhì)安全化程度直接影響著其應(yīng)用。目前常用的RDX制造方法為直接硝解法與醋酸酐法。在RDX的直接硝解法生產(chǎn)制造過程中,由于硝化液的氧化過程與RDX的結(jié)晶過程同時進(jìn)行,極易在結(jié)晶過程中于RDX晶體周圍包裹著含酸的氣泡[4-5]。此外,有學(xué)者[5-6]對RDX等一類硝胺炸藥的晶體形貌進(jìn)行表征發(fā)現(xiàn),其晶體存在著孔洞、凹坑麻點、溝槽、裂紋等缺陷,而且在這些缺陷中均充滿著硝酸的母液[7-8]。由此綜合導(dǎo)致了直接硝解法制造的RDX晶體含酸量較高。酸性粒子催化使得RDX貯存的危險性有所增加[9]。因此,研究RDX的驅(qū)酸洗滌動力學(xué)可為RDX驅(qū)酸洗滌的工藝參數(shù)優(yōu)化與選擇提供理論依據(jù)。

    目前針對RDX驅(qū)酸洗滌過程的研究主要集中在機理、數(shù)值模擬[10]、動力學(xué)等方面。如張幺玄等[11]對超聲波驅(qū)酸洗滌RDX的機理進(jìn)行了探討,認(rèn)為超聲波清洗器內(nèi)的空化氣泡在RDX界面上發(fā)生了非對稱性潰滅,從而產(chǎn)生局部高速射流與渦流作用,強化了RDX驅(qū)酸清洗效果。2014年羅志龍等[12]基于RDX的超聲驅(qū)酸洗滌數(shù)據(jù),假設(shè)驅(qū)酸過程符合一級動力學(xué)規(guī)律,采用最小二乘法進(jìn)行了線性擬合,獲得了相關(guān)動力學(xué)參數(shù)。在此基礎(chǔ)上,2016年李婷婷等[13]對比分析了磁力常規(guī)清洗與超聲波清洗的動力學(xué)參數(shù),常規(guī)清洗的表觀活化能大于超聲波清洗的表觀活化能,表明超聲洗滌更容易驅(qū)酸。然而,對其驅(qū)酸洗滌的動力學(xué)研究僅停留于假設(shè)的一級反應(yīng),由擬合效果來看,動力學(xué)模型描述的準(zhǔn)確性不夠高。

    基于此,本研究以吸附、擴散等理論為出發(fā)點,基于相應(yīng)理論的動力學(xué)模型進(jìn)行RDX驅(qū)酸洗滌的動力學(xué)探討,以獲得較為精確的動力學(xué)模型來較好地描述RDX驅(qū)酸洗滌過程,以期為RDX制造工藝參數(shù)優(yōu)化提供較為準(zhǔn)確的理論依據(jù),為相關(guān)的工藝改進(jìn)與連續(xù)化設(shè)計提供基礎(chǔ)。

    1 驅(qū)酸洗滌動力學(xué)模型的建立

    基于物理吸附與化學(xué)吸附的定義可知,由于RDX晶體表面存在著孔洞、凹坑麻點、溝槽、裂紋等缺陷且其晶體存在著非球形的菱與角,在含有硝酸的水溶液中RDX表面存在著不飽和力場。因此RDX表面自發(fā)地吸附硝酸與水。此外,因RDX與水、硝酸不存在化學(xué)反應(yīng),無任何形式的化學(xué)鍵,所以其吸附過程應(yīng)為物理吸附,吸附力為RDX表面分子與硝酸分子、水分子間的分子間作用力。此外,陳剛[14]在對RDX表面吸附溶劑分子的研究過程中發(fā)現(xiàn),溶劑分子與RDX晶面之間的吸附能主要是由靜電作用能和范德華作用能組成,且溶劑分子對RDX晶面的吸附是一個放熱過程。由于硝酸水溶液可作為RDX的溶劑,此研究結(jié)果同樣說明了RDX與硝酸水溶液之間的吸附不存在任何的化學(xué)變化。

    由此可以推測,RDX表面酸脫除機制主要由兩個步驟組成:(1)RDX表面吸附的酸在吸收一定能量(能量大于吸附能)后從固體表面脫離;(2)脫離的酸由RDX晶體附近借助一定濃度梯度擴散遷移至溶液中。因此,在驅(qū)酸洗滌過程中,速率控制步驟有可能是其中一種或是兩種的綜合。

    目前,國內(nèi)外對固液吸附動力學(xué)模型與機理的研究有很多。在探索與實踐的歷程中,逐漸形成了6種典型且常用的動力學(xué)模型[15-17],包括基于溶液濃度的擬一階與擬二階動力學(xué)模型、Elovich方程模型、Freundlich方程模型、Avrami方程模型、H-S分?jǐn)?shù)階動力學(xué)模型。結(jié)合脫附與吸附之間的可逆關(guān)系,此6種常用的脫附動力學(xué)模型如表1所示。

    表1 常用的脫附動力學(xué)模型Table 1 Common desorption kinetic models

    目前對質(zhì)量擴散的動力學(xué)機理及其函數(shù)已有很多,如表2所示[18]。其擴散動力學(xué)模型為:

    表2 常用的質(zhì)量擴散動力學(xué)機理函數(shù)Table 2 Commonmass diffusion kinetic mechanism functions

    (1)

    (2)

    式中:α為轉(zhuǎn)化率,即RDX表面酸的脫除率,見下式:

    α=(qmax-qi)/(qmax-qe)

    (3)

    式中:qe為RDX表面吸附與脫附達(dá)到動態(tài)平衡時的酸度。

    此外,不論是以脫附模型為主、還是以擴散模型為主,假設(shè)RDX的驅(qū)酸洗滌速率常數(shù)與溫度的關(guān)系服從阿倫尼烏斯方程,即:

    k=Aexp(-E/RT)

    (4)

    式中:E為RDX表面脫酸的活化能。

    在上述動力學(xué)模型進(jìn)行非線性擬合的同時,需建立評選準(zhǔn)則以選擇較優(yōu)的RDX表面脫附動力學(xué)模型。因此采用全局優(yōu)化的通用算法(UGO)將實驗數(shù)據(jù)與脫附動力學(xué)模型進(jìn)行擬合。相關(guān)系數(shù)是選擇的首要衡量標(biāo)準(zhǔn)。其次的衡量標(biāo)準(zhǔn)為均方根誤差(RMSE)、確定系數(shù)(DC)和殘差平方和(RSS)。相關(guān)系數(shù)與DC越大,RMSE與RSS越小,則說明模型擬合程度越高。這些參數(shù)可由式(5)~(7)獲得:

    (5)

    (6)

    (7)

    式中:qexp,i為第i個實測轉(zhuǎn)化率值;qpre,i為第i個預(yù)測的轉(zhuǎn)化率值;N為觀測的樣本數(shù)。

    2 實 驗

    2.1 實驗材料

    RDX粗品,平均粒徑1.3mm,酸度約1.072%,自制;氫氧化鈉,分析純,國藥集團(tuán)化學(xué)試劑有限公司;酚酞,分析純,上海三愛思試劑有限公司。

    2.2 實驗裝置

    經(jīng)過實驗室直接硝解法合成RDX的步驟與抽濾后,得到RDX粗品。然后將裝有一定量的RDX粗品與去離子水的燒杯放置于集熱式磁力攪拌加熱器與超聲波清洗機內(nèi),在不同的洗滌溫度下分別進(jìn)行不同清洗時間的實驗。隨后將驅(qū)酸洗滌實驗后的RDX經(jīng)抽濾后得到含酸濾餅。最后,將濾餅均勻分散于去離子水中,滴加酚酞指示劑,并使用裝有氫氧化鈉標(biāo)準(zhǔn)溶液的微量滴定管進(jìn)行酸度測定。濾餅酸度計算公式為:

    (8)

    式中:q為濾餅酸度(%);c為氫氧化鈉標(biāo)準(zhǔn)溶液的濃度(mol/L);V1為空白實驗所需氫氧化鈉標(biāo)準(zhǔn)溶液的體積(L);V2為滴定濾餅所需氫氧化鈉標(biāo)準(zhǔn)溶液的體積(L);m為濾餅質(zhì)量(g)。

    實驗流程與裝置如圖1所示。

    圖1 RDX驅(qū)酸洗滌動力學(xué)實驗流程與裝置Fig.1 The kinetic test process and device of washing RDX

    2.3 實驗過程

    為了研究RDX驅(qū)酸洗滌動力學(xué),在20、25、30、35、40℃的洗滌溫度下分別進(jìn)行不同清洗時間的實驗(超聲清洗時間分別為10、20、40、60、80min,加熱攪拌清洗時間分別為10、40、70、100、130min)。為了與李婷婷等[13]研究的結(jié)果進(jìn)行對比分析,其他實驗條件均為一致:兩種清洗方式的液固比均為2∶1,超聲清洗頻率為25kHz,超聲功率為85%,攪拌速率為150r/min。然后基于不同溫度下不同清洗時間的酸度值及6個脫附模型和9個擴散模型,采用共享參數(shù)模式下非線性擬合方法進(jìn)行驅(qū)酸洗滌動力學(xué)參數(shù)的求解。其中,針對脫附模型中可轉(zhuǎn)換為積分形式的模型1、2、5、6來說,采用共享參數(shù)模式下非線性擬合方法進(jìn)行求解;針對不可轉(zhuǎn)換積分形式的模型3來說,可采用共享參數(shù)模式下常微分方程參數(shù)擬合方法進(jìn)行求解;針對不可轉(zhuǎn)換積分形式的模型4來說,微分方程求解器無法進(jìn)行計算。此外,針對9個擴散模型來說,基于不同機理函數(shù)的積分形式將實測數(shù)據(jù)進(jìn)行轉(zhuǎn)化,隨后采用共享參數(shù)模式下非線性擬合方法進(jìn)行求解。

    共享參數(shù)模式下非線性擬合方法,即在具有共同的擬合參數(shù)且含有隨不同數(shù)據(jù)文件而變化的常數(shù)情況下采用的一種非線性擬合方法,可通過1stOpt軟件實現(xiàn)。以脫附模型6與磁力攪拌加熱清洗方式獲得的5種溫度條件下的數(shù)據(jù)為例進(jìn)行說明,其共有4個擬合參數(shù)(脫附活化能E、指前因子A、脫附時間指數(shù)n、脫附模型常數(shù)m),洗滌溫度(20、25、30、35、40℃)即為隨不同數(shù)據(jù)而變化的常數(shù)。

    3 結(jié)果與討論

    3.1 脫附模型的優(yōu)選

    不同溫度下兩種清洗方式的脫附模型擬合狀態(tài)參數(shù)間的對比如圖2所示?;赗與DC參數(shù)均為最大且RMSE與RSS均為最小的選擇原則,黑色柱形則表示最優(yōu)的模型。擬合較優(yōu)模型的動力學(xué)參數(shù)見表3。

    圖2 不同溫度下兩種清洗方式的脫附模型擬合狀態(tài)參數(shù)對比Fig.2 Comparison of fitting parameters of desorption models for two washing methods at different temperatures

    表3 兩種清洗方式的較優(yōu)脫附模型擬合參數(shù)Table 3 The better fitting parameters of desorption models for two washing methods

    由圖2可知,在20、25、30、35、40℃的清洗溫度下,磁力攪拌清洗的最優(yōu)脫附模型分別為模型6、6、5、6、6,而超聲波清洗的最優(yōu)脫附模型分別為模型1、1、1、6、6。此外,還可以看出,對磁力攪拌清洗來說,模型5、6的擬合狀態(tài)參數(shù)大都相近;而對超聲波清洗來說,模型1、5、6的擬合狀態(tài)參數(shù)大都相近,由表3可看出,模型5、6的脫附時間指數(shù)n均接近于1且模型6的脫附常數(shù)m也接近于1。綜上所述,模型6(即H-S分?jǐn)?shù)階脫附動力學(xué)模型)能夠更好地、統(tǒng)一地描述兩種清洗方式的驅(qū)酸洗滌過程。對磁力攪拌來說,擬合決定系數(shù)R與DC均在0.999以上、RMSE均小于0.012、RSS均小于0.0008;對超聲波清洗來說,擬合的R均在0.998以上、DC均在0.993以上、RMSE均小于0.031、RSS均小于0.006。

    3.2 擴散模型的優(yōu)選

    不同溫度下兩種清洗方式的擴散模型擬合狀態(tài)參數(shù)(R、DC、RMSE、RSS)分布云圖如圖3與圖4所示。其中模型9的機理函數(shù)積分形式將擬合的因變量進(jìn)一步擴大,基于RMSE與RSS的定義與計算公式,因變量的擴大直接導(dǎo)致RMSE與RSS的值擴大,使得RMSE與RSS的衡量準(zhǔn)則失效,故在圖3與圖4中,舍去模型9的RMSE與RSS值。此外,由于在擴散模型擬合中,因變量的值與機理函數(shù)的積分形式關(guān)系密切,因此擴散模型的優(yōu)選性僅能使用作為相對值的R與DC來衡量。

    圖3 不同溫度下磁力攪拌驅(qū)酸洗滌擴散模型的擬合狀態(tài)參數(shù)云圖Fig.3 Cloud map of fitting parameters of diffusion models for the magnetic stirring cleaning at different temperatures

    圖4 不同溫度下超聲驅(qū)酸洗滌擴散模型的擬合狀態(tài)參數(shù)云圖Fig.4 Cloud map of fitting parameters of diffusion models for the ultrasonic cleaning at different temperatures

    如圖3所示,針對不同洗滌溫度下的磁力攪拌清洗來說,相關(guān)系數(shù)R的云圖中僅在模型6位置出現(xiàn)“紅色縱向貫穿”的現(xiàn)象,而在DC的云圖中出現(xiàn)“紅色縱向貫穿”的位置位于模型4、6、7,在RMSE、RSS的云圖中出現(xiàn)“藍(lán)色縱向貫穿”的位置位于模型6、7、8。因此用來描述磁力攪拌驅(qū)酸過程的最優(yōu)模型為模型6。

    由圖4可知,在R云圖中,雖未出現(xiàn)“紅色縱向貫穿”,但在模型6的位置時縱向紅色區(qū)域較多,在DC的云圖中出現(xiàn)“紅色縱向貫穿”的位置位于模型4、6、7,在RMSE的云圖中藍(lán)色區(qū)域的現(xiàn)象與圖4相應(yīng)的云圖一致,在RSS云圖中,出現(xiàn)“藍(lán)色縱向貫穿”的位置位于模型4、6、7、8。因此,最優(yōu)的超聲波驅(qū)酸洗滌擴散模型為模型6。

    綜上所述,擴散模型6,即Jander方程V類,能夠更好且統(tǒng)一地描述兩種清洗方式驅(qū)酸洗滌過程,其動力學(xué)相關(guān)參數(shù)為:當(dāng)洗滌方式為磁力攪拌時,活化能為13.49kJ/mol,指前因子為0.973min-1;當(dāng)洗滌方式為超聲波清洗時,活化能為11.02kJ/mol,指前因子為0.569min-1。

    3.3 最優(yōu)驅(qū)酸洗滌模型的選擇

    由上述對脫附模型與擴散模型的分析可得,兩種較優(yōu)的動力學(xué)方程如下式所示:

    (9)

    (10)

    將上述兩式變化為:

    (11)

    qi=qmax{1-[Aexp(-E/RT)t]1/2}3

    (12)

    基于兩種清洗方式的實驗數(shù)據(jù),將以脫附模型6與擴散模型6得到的預(yù)測曲線進(jìn)行對比分析,其結(jié)果如圖5和圖6所示。

    圖5 磁力攪拌清洗數(shù)據(jù)與脫附模型及擴散模型對比Fig.5 Comparison of data from experiments, the desorption model and the diffusion model for the magnetic stirring cleaning

    圖6 超聲波清洗數(shù)據(jù)與脫附模型及擴散模型對比Fig.6 Comparison of data from experiments, the desorption model and the diffusion model for the ultrasonic cleaning

    由圖5和圖6可知,脫附模型6的預(yù)測曲線均更接近于實驗數(shù)據(jù),且擬合的相關(guān)系數(shù)均在0.998以上,高于李婷婷等[13]的擬合結(jié)果。結(jié)果表明,用來描述的RDX驅(qū)酸洗滌過程的最優(yōu)動力學(xué)模型為脫附模型6,即H-S分?jǐn)?shù)階脫附動力學(xué)模型。此外,還可以看出,隨著溫度的升高,擴散模型越接近于脫附模型。這是因為洗滌溫度增加,由于脫附活化能高于擴散活化能,RDX表面脫附速率增長的幅度較大,從而導(dǎo)致其驅(qū)酸洗滌過程的速率控制步驟由脫附為主的綜合過程逐步轉(zhuǎn)化為酸的擴散過程。

    綜上所述,RDX的驅(qū)酸洗滌動力學(xué)模型如式(13)所示:

    (13)

    式中:當(dāng)洗滌方式為磁力攪拌時,E=20.05kJ/mol、A=32.070min-1、n=1.4682、m=1.2984;當(dāng)洗滌方式為超聲波清洗時,E=21.42kJ/mol、A=329.516min-1、n=0.9856、m=1.0545。

    此外,在式(13)中,脫附時間指數(shù)n體現(xiàn)了RDX驅(qū)酸洗滌過程中可能存在的復(fù)雜脫附機理綜合作用。

    將超聲波清洗的動力學(xué)參數(shù)與磁力攪拌的情況相比(見表3),其表觀活化能的值僅增加了1.37kJ/mol,變化較小,而指前因子的值約提高了9倍。雖然超聲波清洗的脫附活化能稍高,但其指前因子的增加非常明顯,兩種清洗方式的酸脫除速率常數(shù)對比如圖7所示。

    圖7 兩種清洗方式的酸脫除速率常數(shù)對比Fig.7 Comparison of rate constants of removal acid for two washing methods

    由圖7可知,當(dāng)洗滌溫度為20~50℃時,與磁力攪拌清洗相比,超聲波清洗時的速率常數(shù)增加了4.86~5.17倍。結(jié)果表明,雖然超聲清洗的脫附能壘稍高,但即將脫附的酸活性位點數(shù)目較大,總體來說,反應(yīng)速率常數(shù)較大。此結(jié)果主要因為,超聲波的引入使其空化效應(yīng)在含RDX粗品的水溶液中進(jìn)行作用,從而引起溶液的湍動效應(yīng)與微擾效應(yīng),即空化效應(yīng)產(chǎn)生的聲沖擊波引起溶液宏觀湍流并于RDX晶體缺陷邊界發(fā)生氣泡破裂引起局部的微射流。而磁力攪拌僅是通過攪拌葉輪轉(zhuǎn)動引起溶液的湍流效應(yīng)。兩者相比,超聲波清洗較大幅度地提高了RDX表面吸附的酸脫效率。

    4 結(jié) 論

    (1)H-S分?jǐn)?shù)階脫附模型能夠較為準(zhǔn)確地描述RDX驅(qū)酸洗滌過程,其擬合相關(guān)系數(shù)均在0.998以上、擬合決定系數(shù)均在0.993以上。

    (2)在本研究的實驗條件下,當(dāng)洗滌方式為磁力攪拌時,脫附活化能為20.05kJ/mol、指前因子為32.070min-1、脫附時間指數(shù)為1.4682、脫附模型常數(shù)為1.2984;當(dāng)洗滌方式為超聲波清洗時,脫附活化能為21.42kJ/mol、指前因子為329.516min-1、脫附時間指數(shù)為0.9856、脫附模型常數(shù)為1.0545。

    (3)當(dāng)洗滌溫度為20~50℃時,與磁力攪拌清洗相比,超聲波清洗時RDX表面吸附的酸脫除速率約增大了4.86~5.17倍,且隨著洗滌溫度的升高,驅(qū)酸洗滌過程的速率控制步驟逐漸由脫附為主的綜合過程逐步轉(zhuǎn)化為酸的擴散過程。

    猜你喜歡
    磁力云圖超聲波
    磁力文件夾
    磁力珠
    制作磁力小車
    磁力不怕水
    成都云圖控股股份有限公司
    基于Niosll高精度超聲波流量計的研究
    電子制作(2018年17期)2018-09-28 01:56:38
    黃強先生作品《雨后松云圖》
    名家名作(2017年3期)2017-09-15 11:13:37
    基于TV-L1分解的紅外云圖超分辨率算法
    蝙蝠的超聲波
    超聲波流量計的研究
    電子制作(2016年21期)2016-05-17 03:52:47
    男人的好看免费观看在线视频 | 亚洲五月色婷婷综合| 一区二区三区国产精品乱码| 成人亚洲精品av一区二区| 久久人人爽av亚洲精品天堂| 久久人妻av系列| 免费一级毛片在线播放高清视频 | 国产高清激情床上av| 国产亚洲欧美精品永久| 午夜久久久在线观看| 欧美大码av| 中文字幕高清在线视频| 日韩精品中文字幕看吧| 国产精品乱码一区二三区的特点 | 91麻豆精品激情在线观看国产| 国产精品久久久av美女十八| 两个人免费观看高清视频| av视频免费观看在线观看| 国产高清激情床上av| 久久精品国产清高在天天线| 99国产精品免费福利视频| 色综合站精品国产| 日本精品一区二区三区蜜桃| 精品乱码久久久久久99久播| 精品人妻在线不人妻| 久久久水蜜桃国产精品网| 黄频高清免费视频| 女警被强在线播放| 真人做人爱边吃奶动态| 人人妻人人澡欧美一区二区 | 黄色a级毛片大全视频| 校园春色视频在线观看| 中文字幕精品免费在线观看视频| 午夜福利欧美成人| 最近最新免费中文字幕在线| 级片在线观看| 国产一区二区三区视频了| 亚洲,欧美精品.| 1024视频免费在线观看| 国产成人av激情在线播放| 动漫黄色视频在线观看| 欧美黑人欧美精品刺激| 女警被强在线播放| 一本大道久久a久久精品| 国产三级在线视频| 99国产精品一区二区蜜桃av| 午夜福利成人在线免费观看| 变态另类成人亚洲欧美熟女 | 国产成人av激情在线播放| 国产午夜福利久久久久久| 国产精品免费一区二区三区在线| 18禁裸乳无遮挡免费网站照片 | 成人国产一区最新在线观看| 啦啦啦 在线观看视频| 欧美色视频一区免费| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 久久天堂一区二区三区四区| 香蕉国产在线看| 欧美久久黑人一区二区| 成人av一区二区三区在线看| or卡值多少钱| 校园春色视频在线观看| 777久久人妻少妇嫩草av网站| 男人舔女人下体高潮全视频| 日韩大码丰满熟妇| 亚洲欧美日韩高清在线视频| 日韩精品青青久久久久久| 成人亚洲精品一区在线观看| 中文字幕色久视频| 又紧又爽又黄一区二区| 亚洲精品美女久久久久99蜜臀| 欧美色视频一区免费| 99精品在免费线老司机午夜| 午夜免费成人在线视频| 男女下面插进去视频免费观看| 成人三级做爰电影| 九色亚洲精品在线播放| www.999成人在线观看| 久9热在线精品视频| 国产在线精品亚洲第一网站| netflix在线观看网站| 人妻丰满熟妇av一区二区三区| 欧美黑人欧美精品刺激| 18禁国产床啪视频网站| 中文字幕人妻丝袜一区二区| 视频区欧美日本亚洲| av视频在线观看入口| 波多野结衣一区麻豆| 亚洲欧美精品综合久久99| 久久人妻av系列| 国产成人免费无遮挡视频| av免费在线观看网站| 午夜视频精品福利| 精品人妻1区二区| 亚洲,欧美精品.| av在线天堂中文字幕| 中亚洲国语对白在线视频| 国产成人av教育| 国产av一区在线观看免费| 免费观看精品视频网站| 国产精品免费视频内射| 亚洲少妇的诱惑av| 国产精品 国内视频| 精品欧美一区二区三区在线| 国产成人精品久久二区二区免费| 精品久久久久久久人妻蜜臀av | 午夜久久久在线观看| 午夜日韩欧美国产| 国产一区二区三区视频了| 亚洲一区二区三区不卡视频| 亚洲中文字幕一区二区三区有码在线看 | 成人国语在线视频| 日韩欧美国产在线观看| 亚洲 欧美一区二区三区| 午夜激情av网站| 大型黄色视频在线免费观看| 久热这里只有精品99| 国产精品九九99| av网站免费在线观看视频| av福利片在线| 成人国语在线视频| 精品第一国产精品| 在线观看舔阴道视频| 一区二区三区激情视频| 变态另类丝袜制服| 欧美中文日本在线观看视频| 香蕉久久夜色| 一区二区日韩欧美中文字幕| 十八禁人妻一区二区| 免费观看精品视频网站| 亚洲第一av免费看| 99在线视频只有这里精品首页| 精品卡一卡二卡四卡免费| 免费少妇av软件| 一边摸一边抽搐一进一小说| 久久久久久久久中文| 亚洲第一av免费看| 久久性视频一级片| 亚洲av片天天在线观看| 免费av毛片视频| 欧美乱妇无乱码| 国产99久久九九免费精品| 国产精品久久久久久亚洲av鲁大| 免费看十八禁软件| 国产成人精品久久二区二区免费| 男女做爰动态图高潮gif福利片 | 精品国内亚洲2022精品成人| 精品人妻1区二区| 免费人成视频x8x8入口观看| 国产1区2区3区精品| 成人三级黄色视频| 一区在线观看完整版| 亚洲午夜理论影院| 亚洲视频免费观看视频| av福利片在线| 国产精品乱码一区二三区的特点 | 丝袜美足系列| 丰满的人妻完整版| 国产成人av教育| 日韩精品中文字幕看吧| 亚洲精品粉嫩美女一区| 日韩欧美在线二视频| 身体一侧抽搐| 欧美在线一区亚洲| 精品国内亚洲2022精品成人| 一级a爱片免费观看的视频| 久久久久九九精品影院| 欧美 亚洲 国产 日韩一| 欧美老熟妇乱子伦牲交| 国产成人欧美在线观看| 欧美亚洲日本最大视频资源| 日日干狠狠操夜夜爽| 9色porny在线观看| 中文字幕人成人乱码亚洲影| 男人操女人黄网站| 老司机在亚洲福利影院| 午夜福利免费观看在线| 亚洲精品国产一区二区精华液| 激情视频va一区二区三区| 亚洲狠狠婷婷综合久久图片| 久久伊人香网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产欧美一区二区综合| 黑人操中国人逼视频| 亚洲精品久久成人aⅴ小说| 久久久国产成人免费| 不卡av一区二区三区| 亚洲激情在线av| 很黄的视频免费| 一级a爱视频在线免费观看| 国产三级黄色录像| 午夜成年电影在线免费观看| 男女午夜视频在线观看| 99在线视频只有这里精品首页| 国产免费男女视频| 久久香蕉国产精品| 久久影院123| 日本五十路高清| 1024视频免费在线观看| 99国产精品一区二区三区| 美女高潮喷水抽搐中文字幕| 欧美另类亚洲清纯唯美| 婷婷精品国产亚洲av在线| 一级片免费观看大全| 一级毛片精品| av欧美777| 国产精品美女特级片免费视频播放器 | 国产一区二区在线av高清观看| 操美女的视频在线观看| 精品国产国语对白av| 老司机在亚洲福利影院| 免费在线观看日本一区| 亚洲国产精品合色在线| 亚洲成人精品中文字幕电影| 黄色视频,在线免费观看| 两性夫妻黄色片| 久久热在线av| 国产乱人伦免费视频| 国产又色又爽无遮挡免费看| 日本免费a在线| 久久久久国产一级毛片高清牌| 搡老熟女国产l中国老女人| 国产亚洲精品第一综合不卡| 国产一区二区激情短视频| av福利片在线| 久久国产亚洲av麻豆专区| 不卡一级毛片| 19禁男女啪啪无遮挡网站| 午夜精品久久久久久毛片777| 两个人看的免费小视频| 亚洲一区二区三区色噜噜| 国产精品一区二区三区四区久久 | 两个人视频免费观看高清| 国产av一区在线观看免费| 午夜福利影视在线免费观看| 久久热在线av| 老司机在亚洲福利影院| 免费少妇av软件| 亚洲av熟女| 亚洲色图 男人天堂 中文字幕| 欧美性长视频在线观看| 色综合站精品国产| 一级毛片女人18水好多| 可以免费在线观看a视频的电影网站| 精品国产亚洲在线| 亚洲精品久久国产高清桃花| 又黄又爽又免费观看的视频| 国产av精品麻豆| 欧美一级a爱片免费观看看 | 亚洲成a人片在线一区二区| 9191精品国产免费久久| 久久久久亚洲av毛片大全| 免费观看人在逋| 久久精品国产99精品国产亚洲性色 | 好看av亚洲va欧美ⅴa在| 搡老岳熟女国产| 亚洲五月婷婷丁香| 色播在线永久视频| 亚洲国产看品久久| 男女下面插进去视频免费观看| 精品欧美一区二区三区在线| 曰老女人黄片| 久热这里只有精品99| 欧美乱妇无乱码| 最近最新中文字幕大全电影3 | 丝袜美足系列| 亚洲精品一区av在线观看| 男女之事视频高清在线观看| 老司机午夜福利在线观看视频| 黑人操中国人逼视频| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 法律面前人人平等表现在哪些方面| 一级黄色大片毛片| 亚洲无线在线观看| 亚洲avbb在线观看| 亚洲熟女毛片儿| 成人亚洲精品一区在线观看| 久久精品国产99精品国产亚洲性色 | 激情视频va一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 波多野结衣av一区二区av| www.999成人在线观看| 久99久视频精品免费| 香蕉久久夜色| 亚洲国产毛片av蜜桃av| 亚洲av电影在线进入| 精品一区二区三区四区五区乱码| 岛国在线观看网站| 无遮挡黄片免费观看| 中文字幕人妻熟女乱码| 国语自产精品视频在线第100页| 国产成年人精品一区二区| 亚洲人成电影免费在线| 咕卡用的链子| 亚洲国产毛片av蜜桃av| 亚洲av熟女| 国产片内射在线| 精品国产亚洲在线| 成人免费观看视频高清| 久久久国产成人精品二区| 女性被躁到高潮视频| 亚洲成人国产一区在线观看| 可以在线观看的亚洲视频| 精品国产国语对白av| 欧美在线黄色| 欧美成人一区二区免费高清观看 | 国产亚洲av高清不卡| 精品人妻在线不人妻| 后天国语完整版免费观看| 丁香六月欧美| www.自偷自拍.com| 中文字幕人妻丝袜一区二区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲三区欧美一区| 亚洲精品粉嫩美女一区| 亚洲免费av在线视频| 777久久人妻少妇嫩草av网站| 国产私拍福利视频在线观看| 久久久久九九精品影院| 日本在线视频免费播放| 757午夜福利合集在线观看| 欧美激情 高清一区二区三区| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久久久毛片| 露出奶头的视频| www国产在线视频色| 黑人巨大精品欧美一区二区mp4| 国产欧美日韩精品亚洲av| 女人精品久久久久毛片| 精品不卡国产一区二区三区| 国产黄a三级三级三级人| 亚洲欧美激情综合另类| 欧美日韩精品网址| 亚洲熟妇中文字幕五十中出| 很黄的视频免费| 成人av一区二区三区在线看| 一进一出抽搐动态| 欧美+亚洲+日韩+国产| 国产亚洲av嫩草精品影院| 欧美乱色亚洲激情| avwww免费| 老汉色∧v一级毛片| 亚洲av成人av| 欧美黄色淫秽网站| 激情在线观看视频在线高清| 免费高清视频大片| 人妻久久中文字幕网| 亚洲国产欧美日韩在线播放| 国产亚洲精品久久久久5区| 9色porny在线观看| 巨乳人妻的诱惑在线观看| 国产欧美日韩一区二区精品| 国产免费男女视频| 19禁男女啪啪无遮挡网站| 国产精品久久电影中文字幕| 精品电影一区二区在线| av超薄肉色丝袜交足视频| 午夜日韩欧美国产| 亚洲激情在线av| 首页视频小说图片口味搜索| 精品国产超薄肉色丝袜足j| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩一区二区三区精品不卡| 成人av一区二区三区在线看| 欧美色视频一区免费| 国产99白浆流出| 一边摸一边抽搐一进一小说| 午夜两性在线视频| tocl精华| 天天躁夜夜躁狠狠躁躁| av片东京热男人的天堂| 在线视频色国产色| 欧美中文综合在线视频| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| av视频免费观看在线观看| 黄色丝袜av网址大全| 午夜免费观看网址| 免费在线观看日本一区| 亚洲人成77777在线视频| 亚洲国产日韩欧美精品在线观看 | 禁无遮挡网站| 午夜福利高清视频| 老司机午夜十八禁免费视频| 亚洲午夜理论影院| 精品国产乱码久久久久久男人| 久久精品亚洲精品国产色婷小说| 国产成人精品无人区| 老司机在亚洲福利影院| 桃红色精品国产亚洲av| 久久久久久人人人人人| 欧美色欧美亚洲另类二区 | 久久伊人香网站| 丰满的人妻完整版| 人妻丰满熟妇av一区二区三区| 亚洲熟妇熟女久久| 韩国av一区二区三区四区| 制服丝袜大香蕉在线| 欧美性长视频在线观看| 午夜福利在线观看吧| 黄网站色视频无遮挡免费观看| 欧美成狂野欧美在线观看| 999精品在线视频| 操出白浆在线播放| 亚洲美女黄片视频| 免费久久久久久久精品成人欧美视频| 91麻豆精品激情在线观看国产| 久久久久久久久免费视频了| 久久久久久国产a免费观看| 久久久久久久久中文| 丝袜美腿诱惑在线| 亚洲午夜精品一区,二区,三区| 久久久久久久精品吃奶| 亚洲欧美精品综合一区二区三区| 国产亚洲av嫩草精品影院| 亚洲三区欧美一区| 亚洲一卡2卡3卡4卡5卡精品中文| 此物有八面人人有两片| 在线av久久热| 亚洲片人在线观看| 国产在线精品亚洲第一网站| 亚洲无线在线观看| 黑人操中国人逼视频| 久久人妻熟女aⅴ| 亚洲最大成人中文| 亚洲精品中文字幕一二三四区| 香蕉久久夜色| 亚洲五月天丁香| 黄网站色视频无遮挡免费观看| 欧美成狂野欧美在线观看| ponron亚洲| 非洲黑人性xxxx精品又粗又长| 午夜影院日韩av| 人妻久久中文字幕网| 他把我摸到了高潮在线观看| 亚洲色图综合在线观看| 免费观看人在逋| 老司机在亚洲福利影院| 精品欧美国产一区二区三| 91字幕亚洲| www国产在线视频色| 国产av一区二区精品久久| 精品久久久久久久久久免费视频| 一个人免费在线观看的高清视频| 女性被躁到高潮视频| 又大又爽又粗| 一进一出抽搐gif免费好疼| 亚洲成a人片在线一区二区| 中文字幕人成人乱码亚洲影| 9色porny在线观看| 亚洲狠狠婷婷综合久久图片| 国产精品久久久久久人妻精品电影| 亚洲一区中文字幕在线| 国产色视频综合| 韩国精品一区二区三区| 成人亚洲精品av一区二区| 久久中文看片网| 看片在线看免费视频| 日本免费a在线| 在线观看66精品国产| 后天国语完整版免费观看| 婷婷精品国产亚洲av在线| 精品午夜福利视频在线观看一区| 欧美黄色片欧美黄色片| 咕卡用的链子| 中国美女看黄片| 欧美一区二区精品小视频在线| 欧美丝袜亚洲另类 | 人妻丰满熟妇av一区二区三区| 露出奶头的视频| 又黄又爽又免费观看的视频| 香蕉国产在线看| 亚洲av成人不卡在线观看播放网| 成人亚洲精品av一区二区| 日韩欧美三级三区| 精品久久久久久久久久免费视频| 大陆偷拍与自拍| 国语自产精品视频在线第100页| 亚洲色图 男人天堂 中文字幕| 婷婷精品国产亚洲av在线| 国产激情久久老熟女| 国产精品一区二区三区四区久久 | 又紧又爽又黄一区二区| 国产麻豆成人av免费视频| 亚洲第一av免费看| 久久久久久久久免费视频了| 久久中文字幕一级| 午夜两性在线视频| 性少妇av在线| 精品人妻在线不人妻| 国产99久久九九免费精品| av在线天堂中文字幕| 午夜激情av网站| 久久性视频一级片| 午夜福利18| 亚洲国产中文字幕在线视频| 岛国视频午夜一区免费看| 黑人操中国人逼视频| 色精品久久人妻99蜜桃| 久久久久久人人人人人| av片东京热男人的天堂| 啦啦啦免费观看视频1| 黄片播放在线免费| 国内毛片毛片毛片毛片毛片| 亚洲成人免费电影在线观看| 性少妇av在线| 国产私拍福利视频在线观看| 91在线观看av| 十八禁人妻一区二区| 国产高清视频在线播放一区| 可以免费在线观看a视频的电影网站| 亚洲第一青青草原| 成人永久免费在线观看视频| 又大又爽又粗| 乱人伦中国视频| 亚洲国产中文字幕在线视频| 久久香蕉国产精品| 欧美日韩瑟瑟在线播放| 欧美色视频一区免费| 丁香欧美五月| 黄色a级毛片大全视频| 午夜亚洲福利在线播放| 国产一区二区三区视频了| 国产区一区二久久| 国产不卡一卡二| 亚洲国产精品久久男人天堂| 久久婷婷人人爽人人干人人爱 | 国产亚洲精品第一综合不卡| 久久亚洲真实| 欧美国产精品va在线观看不卡| 少妇裸体淫交视频免费看高清 | 国产一区二区三区综合在线观看| 又大又爽又粗| 亚洲欧美日韩高清在线视频| 变态另类丝袜制服| 久久 成人 亚洲| 久久香蕉国产精品| 亚洲中文日韩欧美视频| 国产成年人精品一区二区| svipshipincom国产片| 男女下面进入的视频免费午夜 | 国产欧美日韩精品亚洲av| 夜夜躁狠狠躁天天躁| 宅男免费午夜| 亚洲情色 制服丝袜| 黄色女人牲交| 人妻久久中文字幕网| 亚洲精品久久成人aⅴ小说| 亚洲avbb在线观看| 亚洲欧美日韩无卡精品| 首页视频小说图片口味搜索| 无遮挡黄片免费观看| 操美女的视频在线观看| 91字幕亚洲| x7x7x7水蜜桃| 午夜免费成人在线视频| 极品人妻少妇av视频| 99国产极品粉嫩在线观看| 女人被狂操c到高潮| 美女高潮喷水抽搐中文字幕| 又黄又粗又硬又大视频| 一边摸一边抽搐一进一小说| 午夜福利18| 午夜久久久久精精品| 制服丝袜大香蕉在线| 一边摸一边做爽爽视频免费| 国产单亲对白刺激| 88av欧美| 成人国产综合亚洲| 国产精品美女特级片免费视频播放器 | 中文字幕人妻熟女乱码| 久久午夜综合久久蜜桃| 欧美激情高清一区二区三区| 色哟哟哟哟哟哟| 欧美另类亚洲清纯唯美| 亚洲人成77777在线视频| 国产精品亚洲一级av第二区| 久久久久久亚洲精品国产蜜桃av| 日韩 欧美 亚洲 中文字幕| 一级毛片高清免费大全| 一级片免费观看大全| 最好的美女福利视频网| 中文字幕最新亚洲高清| 18禁美女被吸乳视频| 亚洲 欧美 日韩 在线 免费| 人妻久久中文字幕网| 女性生殖器流出的白浆| 又紧又爽又黄一区二区| 极品人妻少妇av视频| 男女床上黄色一级片免费看| 久久精品亚洲精品国产色婷小说| 欧美人与性动交α欧美精品济南到| 搡老岳熟女国产| 成人精品一区二区免费| 日韩欧美免费精品| 麻豆久久精品国产亚洲av| 久久久国产成人精品二区| 久久香蕉国产精品| 一级黄色大片毛片| 老司机午夜福利在线观看视频| 色综合站精品国产| 一二三四在线观看免费中文在| 性欧美人与动物交配| 99国产精品免费福利视频| 国产精品98久久久久久宅男小说| 免费不卡黄色视频| 亚洲av成人av| 国产精品一区二区精品视频观看| 夜夜夜夜夜久久久久| 午夜福利,免费看| 色老头精品视频在线观看| 日本vs欧美在线观看视频| 黄片小视频在线播放|