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

    動(dòng)載擾動(dòng)下深井護(hù)巷煤柱失穩(wěn)破壞機(jī)理

    2023-05-23 04:04:14邢闖闖寧建國
    煤炭科學(xué)技術(shù) 2023年3期
    關(guān)鍵詞:煤柱撓度塑性

    邢闖闖 ,王 俊 ,寧建國 ,沈 圳

    (1.山東科技大學(xué) 礦山災(zāi)害預(yù)防控制省部共建國家重點(diǎn)實(shí)驗(yàn)室培育基地, 山東 青島 266590;2.山東科技大學(xué) 能源與礦業(yè)工程學(xué)院, 山東 青島 266590)

    0 引 言

    近年來,煤柱型沖擊地壓災(zāi)害事故頻率不斷增加,嚴(yán)重威脅礦井的生產(chǎn)安全[1-2]。例如2016 年8月15 日0 時(shí)33 分,山東能源肥礦集團(tuán)梁寶寺能源有限責(zé)任公司發(fā)生了煤柱型沖擊地壓事故,造成2人死亡。2019 年8 月2 日12 時(shí)24 分,開灤(集團(tuán))有限責(zé)任公司唐山礦業(yè)分公司也發(fā)生了煤柱型沖擊地壓事故,造成了12 人傷亡,直接經(jīng)濟(jì)損失高達(dá)614.024 萬元。實(shí)踐表明,煤柱型沖擊地壓一般具有隱蔽性強(qiáng)、無明顯前兆特征等特點(diǎn),造成設(shè)備破壞、巷道變形和人員傷亡等事故,給礦井的安全生產(chǎn)帶來嚴(yán)重威脅。因此,有必要開展煤柱動(dòng)力失穩(wěn)機(jī)理研究,以期對煤柱型沖擊地壓預(yù)警與防治提供理論依據(jù)。

    國內(nèi)學(xué)者針對煤柱穩(wěn)定性問題進(jìn)行了大量的研究。部分學(xué)者從靜力學(xué)角度進(jìn)行了大量的試驗(yàn)。殷志強(qiáng)等[3]通過分離式霍普金森壓桿試驗(yàn),分析了不同靜載下煤樣的動(dòng)態(tài)力學(xué)特性,得出峰值應(yīng)力隨靜載增大而增大,峰值應(yīng)變隨靜載增大而減小的結(jié)論。何永琛[4]通過對煤巖組合體試件進(jìn)行單軸壓縮試驗(yàn),獲得了不同煤巖組合體試件的破壞形態(tài)及峰值應(yīng)變等變化特征。艾迪昊等[5]通過對不同粒度的型煤進(jìn)行單軸壓縮試驗(yàn),得出型煤的破壞可分為壓實(shí),初始破裂,加速破裂和卸壓4 個(gè)階段。李成武等[6]通過室內(nèi)試驗(yàn)和現(xiàn)場試驗(yàn)研究了煤體靜載破壞中低頻磁場時(shí)、頻譜特征,并結(jié)合微震信號提出了低頻磁信號的產(chǎn)生機(jī)制。

    也有學(xué)者從動(dòng)力學(xué)的角度,探討了煤柱穩(wěn)定性問題。徐青云等[7]采用理論分析、數(shù)值模擬和現(xiàn)場實(shí)測等方法,指出窄煤柱在受到劇烈采動(dòng)影響時(shí),主要在其中上部區(qū)域發(fā)生變形破壞。姜學(xué)偉[8]通過分析房采區(qū)煤柱在采動(dòng)影響下圍巖受力、煤柱破壞變形狀態(tài)等得出了煤柱保持穩(wěn)定性所必須的彈性區(qū)寬度計(jì)算公式,并作為判斷煤柱是否破壞失穩(wěn)的依據(jù)。李峰等[9]基于LS-DYNA 有限元分析軟件的HJC 本構(gòu)模型,模擬不同動(dòng)載強(qiáng)度下煤柱的動(dòng)態(tài)損傷過程。結(jié)果表明,應(yīng)變率相同時(shí),煤柱受動(dòng)載荷越大,破壞越嚴(yán)重。王方田等[10]將動(dòng)載荷以應(yīng)力波的形式作用在煤體上,來分析不同動(dòng)載荷烈度下煤柱壩體的動(dòng)力響應(yīng)機(jī)制,發(fā)現(xiàn)巷道圍巖塑性區(qū)發(fā)育面積及深度顯著增加,豎向應(yīng)力積聚強(qiáng)度增大,并與動(dòng)載荷烈度成正比。叢利等[11]通過數(shù)值模擬對相鄰工作面回采中和回采后3 次高強(qiáng)度開采擾動(dòng)對煤柱內(nèi)應(yīng)力變化情況的研究,得出煤柱內(nèi)布置的聯(lián)絡(luò)巷會加劇巷道圍巖的高應(yīng)力集中,在外部動(dòng)載擾動(dòng)下,極易達(dá)到?jīng)_擊載荷并誘發(fā)沖擊的結(jié)論。因煤柱在外部動(dòng)載擾動(dòng)下的破壞失穩(wěn)機(jī)理極其復(fù)雜,當(dāng)前的理論研究和實(shí)踐仍需進(jìn)一步完善。因此針對不同尺寸深井護(hù)巷煤柱在外部動(dòng)載擾動(dòng)下的破壞特征和突變失穩(wěn)機(jī)制開展了詳細(xì)研究。

    首先,采用數(shù)值模擬分析了動(dòng)載沖擊下煤柱失穩(wěn)特征;其次,建立了動(dòng)載作用下高靜載煤柱突變失穩(wěn)力學(xué)模型,推導(dǎo)了煤柱穩(wěn)定性判別準(zhǔn)則,揭示了動(dòng)載作用下深井煤柱突變失穩(wěn)機(jī)理,以期對煤柱穩(wěn)定性控制提供參考。

    1 深井護(hù)巷煤柱數(shù)值模型的建立

    1.1 模型建立與參數(shù)選擇

    利用FLAC3D軟件建立“頂板-煤柱-底板”數(shù)值計(jì)算模型,模型寬高比分別設(shè)置為4、6、8 和10,煤層高度保持為3 m。圖1 為寬高比為10 的數(shù)值計(jì)算模型,模型共計(jì)9 層,7 種巖性,模型尺寸為34 m×34 m×39 m,共劃分62 820 個(gè)單元網(wǎng)格。參考相關(guān)數(shù)值計(jì)算經(jīng)驗(yàn)[12-13],頂?shù)装鍘r層采用各向同性彈性模型,煤層采用應(yīng)變軟化模型。煤層節(jié)理、裂隙發(fā)育,室內(nèi)試驗(yàn)獲得的煤樣物理力學(xué)參數(shù)很難直接用于數(shù)值計(jì)算。因此,采用試錯(cuò)法擬合煤柱強(qiáng)度公式以校核煤柱力學(xué)參數(shù)。相關(guān)學(xué)者提出了多種煤柱強(qiáng)度公式[14],其中Salamon 公式能夠準(zhǔn)確描述寬高比為2~20 時(shí)的煤柱強(qiáng)度,計(jì)算結(jié)果精確度能夠滿足現(xiàn)場要求,可用于計(jì)算現(xiàn)場保護(hù)煤柱的強(qiáng)度[15]。因此,采用Salamon 公式校正數(shù)值模型計(jì)算結(jié)果。

    圖1 數(shù)值計(jì)算模型Fig.1 Numerical calculation model

    式中:Sp為保護(hù)煤柱強(qiáng)度,MPa;W為保護(hù)煤柱寬度,m;H為保護(hù)煤柱高度,m。

    采用試錯(cuò)法時(shí),首先建立“頂板-煤柱-底板”數(shù)值計(jì)算模型;其次,在模型頂部施加速度邊界條件以模擬單軸壓縮,獲得煤柱單軸壓縮強(qiáng)度;再者,通過調(diào)整煤層各項(xiàng)力學(xué)參數(shù),使數(shù)值計(jì)算獲得煤柱單軸壓縮強(qiáng)度與Salamon 公式匹配。圖2 為數(shù)值計(jì)算獲得的煤柱強(qiáng)度與Salamon 公式對比圖,由圖2 可知,數(shù)值計(jì)算獲得的煤柱強(qiáng)度與Salamon 公式基本保持一致,因此,采用表1 所示的巖性參數(shù)模擬煤柱的應(yīng)變軟化行為。

    圖2 數(shù)值模擬和Salamon 公式計(jì)算峰值強(qiáng)度對比Fig.2 Comparison of peak strength between numerical simulation and Salamon formula calculation

    表1 煤柱模擬參數(shù)Table 1 Simulation parameters of protective coal pillar

    1.2 邊界條件及模型方案

    以新河煤礦6302 工作面南部軌道集中巷,6301運(yùn)輸巷和6302 軌道巷之間的保護(hù)煤柱為研究對象。煤層埋深980 m,煤的容重為24 kN/m3,故上覆巖層的自重應(yīng)力為24 MPa。

    數(shù)值計(jì)算分為2 步。①在模型頂板施加恒定載荷,開展靜載計(jì)算;②待靜載計(jì)算平衡后,在模型頂板施加速度邊界條件以模擬動(dòng)載。

    靜載計(jì)算時(shí),模型頂部設(shè)置為應(yīng)力邊界,上覆巖層自重采用均布荷載代替,底部為固定邊界,頂部垂直應(yīng)力取上覆巖層的自重應(yīng)力(24 MPa)。在靜載條件下,模型達(dá)到平衡狀態(tài),再在動(dòng)載條件下開展計(jì)算。動(dòng)力學(xué)計(jì)算時(shí),數(shù)值模型頂部邊界自由,底部設(shè)置為固定邊界。

    為了模擬動(dòng)載擾動(dòng),假設(shè)上覆頂板巖層斷裂產(chǎn)生含有多種子波的壓縮波,而任何復(fù)雜的波形都可以由正弦波、余弦波或三角波合成。為簡化動(dòng)載擾動(dòng),選擇正弦應(yīng)力波作為輸入波形。DUBINSKI[16]等現(xiàn)場監(jiān)測發(fā)現(xiàn)應(yīng)力波PPV 為0.03 m/s。而HADJIGEORGIOU 等[17]通過使用炸藥誘發(fā)巖體振動(dòng)所記錄的應(yīng)力波 PPV 不能大于3.3 m/s,因此,為了模擬不同震動(dòng)強(qiáng)度的應(yīng)力波,振幅設(shè)定為0.6、0.8、1.0、1.2、1.4 m/s。

    2 動(dòng)載擾動(dòng)下煤柱失穩(wěn)破壞變化規(guī)律

    選取寬高比分別為4、6、8、10 四種模型,靜載應(yīng)力下平衡后,在煤柱模型的頂部分別施加振幅為0.6、0.8、1.0、1.2、1.4 m/s 的應(yīng)力波(以應(yīng)力波振幅表示動(dòng)載強(qiáng)度),當(dāng)模型再次達(dá)到平衡后,研究應(yīng)力波輸入保護(hù)煤柱破壞失穩(wěn)規(guī)律。

    2.1 煤柱尺寸對煤柱動(dòng)力失穩(wěn)的影響

    由圖3 可知,動(dòng)載沖擊后煤柱塑性區(qū)由表及里向深部擴(kuò)展最終導(dǎo)致煤柱整體失穩(wěn)。相比而言,隨著煤柱尺寸增加,煤柱發(fā)生整體失穩(wěn)所需的動(dòng)載沖擊強(qiáng)度具有一定的差異。以應(yīng)力波振幅為0.6 m/s的模擬結(jié)果為例,以闡述此種差異。

    圖3 不同應(yīng)力波振幅及不同尺寸保護(hù)煤柱塑性區(qū)分布規(guī)律Fig.3 Distribution law of plastic zone of protective coal pillars with different sizes and different stress wave amplitudes

    煤柱寬高比為6 時(shí),靜載作用后,塑性區(qū)寬度為7~8 m,內(nèi)部存在彈性核心。受到振幅為0.6 m/s 應(yīng)力波沖擊后,煤柱彈性核心完全破壞,塑性區(qū)貫穿整個(gè)煤柱,意味著煤柱整體失穩(wěn)。相比而言,煤柱寬高比增加至8 與10 時(shí),煤柱塑性區(qū)僅向深度擴(kuò)展,內(nèi)部依舊存在穩(wěn)定的彈性核心,且彈性核心所占整個(gè)煤柱的比例由21%增加至47%。此外,對比發(fā)現(xiàn),在相同靜載作用下,不同尺寸的煤柱失穩(wěn)啟動(dòng)所需的動(dòng)載強(qiáng)度亦不同。模擬發(fā)現(xiàn),隨著煤柱寬高比由6增加至10 時(shí),煤柱破壞失穩(wěn)所需的動(dòng)載沖擊強(qiáng)度(以應(yīng)力波振幅表征)由0.6 m/s 增加至1.2 m/s,提高了1 倍。由此說明,增加煤柱尺寸有助于提高煤柱的穩(wěn)定性,減少動(dòng)力失穩(wěn)的幾率。

    2.2 動(dòng)載沖擊強(qiáng)度對煤柱動(dòng)力失穩(wěn)的影響

    以煤柱寬高比為10 為例,分析動(dòng)載沖擊強(qiáng)度為0.6、0.8、1.0、1.2、1.4 m/s 時(shí),煤柱失穩(wěn)破壞特征。如圖3c 所示,當(dāng)應(yīng)力波振幅從0.6 m/s 增加至1.0 m/s,煤柱塑性區(qū)由淺部逐步向著中部擴(kuò)展,范圍由8 m增加至13.5 m,煤柱內(nèi)部彈性核區(qū)逐步減少。當(dāng)動(dòng)載沖擊強(qiáng)度增加至1.2 m/s 時(shí),煤柱整體產(chǎn)生塑性破壞,意味著煤柱整體失穩(wěn)。由此說明,動(dòng)載沖擊強(qiáng)度越大,煤柱發(fā)生失穩(wěn)破壞的概率也相對提高。

    通過對比不同振幅應(yīng)力波0.6、0.8、1.0、1.2、1.4 m/s 作用下,煤柱的失穩(wěn)破壞特征,發(fā)現(xiàn)煤柱的穩(wěn)定性與煤柱尺寸和動(dòng)載沖擊強(qiáng)度有關(guān)。相同條件下,尺寸越大,煤柱發(fā)生失穩(wěn)破壞的幾率越??;動(dòng)載強(qiáng)度越大,煤柱破壞范圍及程度越大,即保護(hù)煤柱失穩(wěn)概率越大。

    3 動(dòng)載作用下煤柱突變失穩(wěn)機(jī)制

    3.1 煤柱承載力學(xué)模型

    一般而言,煤柱載荷主要來源于上覆巖層重量和動(dòng)載擾動(dòng)。上覆巖層重量作為靜載,使煤柱產(chǎn)生塑性區(qū)和彈性核區(qū)[18],而后續(xù)動(dòng)載沖擊使塑性區(qū)向煤柱深處擴(kuò)展,導(dǎo)致煤柱發(fā)生整體性破壞,內(nèi)部彈性區(qū)減小而發(fā)生失穩(wěn),煤柱承載力學(xué)模型如圖4 所示。

    圖4 煤柱承載力學(xué)模型Fig.4 Coal pillar bearing mechanical model

    煤柱埋深為H,上覆巖層容重為γ,煤柱高度為Hm,煤柱寬度為B。假設(shè)外部動(dòng)載q0均勻施加至煤柱,動(dòng)荷載在某一時(shí)刻的運(yùn)動(dòng)狀態(tài)可簡化為靜荷載[19],動(dòng)載簡化系數(shù)可以表示為kd,則煤柱承受的載荷為:

    煤柱塑性區(qū)與彈性區(qū)本構(gòu)關(guān)系有所不同,彈性區(qū)煤體呈線彈性關(guān)系,塑性區(qū)內(nèi)呈非線性關(guān)系,具有應(yīng)變軟化特征,煤柱塑性區(qū)內(nèi)本構(gòu)關(guān)系可以表示為

    式中:E為煤柱的初始彈性模量;ε為應(yīng)變;ε0為極限應(yīng)變。

    假設(shè)煤柱塑性區(qū)寬度為B0,則煤柱塑性區(qū)內(nèi)載荷ps與煤柱壓縮量u的關(guān)系可以表示為

    式中:u0為煤柱最大壓縮量,m。

    彈性區(qū)內(nèi),煤體載荷pt與煤柱壓縮量u的關(guān)系為

    3.2 煤柱失穩(wěn)的尖點(diǎn)突變模型

    整個(gè)煤柱系統(tǒng)的總勢能由煤柱自身的應(yīng)變能和外部載荷所做的功組成。由式(4)可得煤柱塑性區(qū)應(yīng)變能V1(u)為

    由式(5)可得煤柱彈性區(qū)應(yīng)變能V2(u)為

    上覆巖層靜載重力勢能和外部動(dòng)載荷做功的勢能函數(shù)V3(u)為

    根據(jù)上述力學(xué)模型的分析,動(dòng)載作用下煤柱系統(tǒng)的總勢能函數(shù)V(u)可以表示為

    以u為狀態(tài)變量進(jìn)行突變理論分析,對V(u)求導(dǎo)得:

    進(jìn)一步簡化,令:

    將式(12)代入式(11)則有:

    以x為狀態(tài)參量,m、n為控制參量。由尖點(diǎn)突變理論可知,當(dāng)m、n滿足分叉集方程時(shí),頂板-煤柱系統(tǒng)處于臨界平衡狀態(tài)。聯(lián)立上式與其拐點(diǎn)方程3x2+m=0,消去狀態(tài)參量x,得到分叉集方程:

    由突變理論可知,當(dāng)Δ>0 時(shí),系統(tǒng)處于穩(wěn)定狀態(tài),煤柱不會失穩(wěn);當(dāng)Δ=0 時(shí),系統(tǒng)處于臨界狀態(tài),受到微小擾動(dòng)會發(fā)生失穩(wěn);當(dāng)Δ<0 時(shí),是系統(tǒng)發(fā)生失穩(wěn)的充要條件。將式(12)代入式(14)得:

    煤柱失穩(wěn)則方程(13)有解,必須滿足Δ≤0,由此得到煤柱失穩(wěn)的必要條件為:

    3.3 煤柱最大壓縮量計(jì)算

    煤柱失穩(wěn)判別式Δ中還含有一個(gè)與煤柱相關(guān)的未知量u0(煤柱最大壓縮量),下面對煤柱最大壓縮量u0進(jìn)行求解。

    煤柱的最大壓縮量是由頂板的最大下沉量決定,如圖5a 和5b 所示。考慮到煤柱的尺寸較大,擬采用彈性薄板理論研究頂板的最大下沉量[20-21]。圖5c所示為直接頂-煤柱彈性薄板模型,其中,2a為煤柱長度,2b為煤柱寬度,q為直接頂上覆動(dòng)靜載等效均布載荷。直接頂可以視為四邊固支的矩形板,其力學(xué)模型如圖5d 所示。

    圖5 頂板-煤柱協(xié)同變形Fig.5 Theoretical process evolution diagram of elastic thin plate

    該模型滿足以下邊界條件:

    根據(jù)彈性力學(xué)中的伽遼金解法,將板的撓度展開為以下表達(dá)式:

    式中:W為直接頂撓度;Ai為撓度系數(shù);fi(x,y)為撓度函數(shù),要求能滿足薄板撓曲面的形狀和幾何邊界條件,而且是連續(xù)的。

    為便于計(jì)算,取一級近似函數(shù):

    顯然該式能滿足薄板撓曲面的形狀和所有幾何邊界條件。根據(jù)伽遼金法對上式進(jìn)行求解,得到撓度的表達(dá)式:

    對撓度方程求極值,在板的中心處取到撓度的最大值:

    煤柱最大壓縮量等于直接頂最大撓度,如圖5e 所示,即u0=Wmax,至此得到了煤柱失穩(wěn)的最終判別式:

    4 深井煤柱穩(wěn)定性判別

    4.1 深井煤柱穩(wěn)定性判別方法

    通過現(xiàn)場調(diào)研獲得煤柱的物理力學(xué)參數(shù),計(jì)算煤柱塑性區(qū)范圍。

    式中:k為應(yīng)力集中系數(shù);c0為黏聚力;φ0為內(nèi)摩擦角;A為側(cè)壓系數(shù)。

    計(jì)算煤柱最大壓縮量:

    最后,利用式(21)計(jì)算Δ值,進(jìn)行穩(wěn)定性判別。

    4.2 工程背景

    新河煤礦位于山東省濟(jì)寧市唐口勘探區(qū)西南部邊緣,目前開采水平為-980 m 延伸水平,延伸水平位于礦井的東南部,東西寬2 km,南北長2.4 km。新河煤礦正在開采6302 工作面,集中大巷分布在6302工作面以南,大巷之間的煤柱對大巷的穩(wěn)定起到支撐保護(hù)作用。利用前文建立的煤柱穩(wěn)定性判別方法,對新河煤礦6302 工作面南部軌道集中巷,6301 運(yùn)輸巷和6302 軌道巷之間的煤柱進(jìn)行穩(wěn)定性進(jìn)行判別,如圖6 所示。

    圖6 煤柱位置示意Fig.6 Position diagram of coal pillar

    1) 計(jì)算煤柱塑性區(qū)范圍:煤柱高2.9 m,黏聚力1.7 MPa,內(nèi)摩擦角29°,埋深980 m,容重24 kN/m3,側(cè)壓系數(shù)1.12,應(yīng)力集中系數(shù)1.97。將煤柱參數(shù)帶入煤柱塑性區(qū)范圍計(jì)算公式:

    2) 計(jì)算煤柱最大壓縮量:大巷保護(hù)煤柱長度46 m,寬度32 m,等效荷載3.75×107N,泥巖彈性模量6.5 GPa,泊松比0.4,頂板厚度6.63 m。

    3) 利用尖點(diǎn)突變模型計(jì)算Δ值,并進(jìn)行煤柱穩(wěn)定性判別:

    根據(jù)煤柱穩(wěn)定性判別公式可知,Δ=1.58>0,該煤柱能保持穩(wěn)定。

    5 結(jié) 論

    1) 數(shù)值模擬研究表明,煤柱受到外部動(dòng)載擾動(dòng)后,塑性區(qū)向中部擴(kuò)展,彈性核心占比減小,導(dǎo)致煤柱破壞失穩(wěn)。煤柱動(dòng)力失穩(wěn)與煤柱尺寸和外界動(dòng)載沖擊強(qiáng)度有關(guān)。增大煤柱的寬度能夠提高煤柱的穩(wěn)定性,減小煤柱破壞失穩(wěn)的概率。相同條件下,煤柱破壞失穩(wěn)的幾率隨著動(dòng)載強(qiáng)度的增大而增大。

    2)基于彈性薄板理論,建立了直接頂-煤柱彈性薄板模型,推導(dǎo)了直接頂?shù)膿隙确匠?,獲得了煤柱的最大壓縮量表達(dá)式?;诩恻c(diǎn)突變理論,建立了煤柱尖點(diǎn)突變模型,推導(dǎo)了主要受塑性區(qū)寬度和煤柱壓縮量影響的煤柱突變失穩(wěn)判別式,為煤柱突變失穩(wěn)提供了力學(xué)判據(jù)。

    3) 基于煤柱突變失穩(wěn)判別式,對新河煤礦6302工作面南部軌道集中巷,6301 運(yùn)輸巷和6302 軌道巷之間的保護(hù)煤柱進(jìn)行了驗(yàn)證。結(jié)果表明,新河煤礦6302 工作面南部集中大巷保護(hù)煤柱能夠保持穩(wěn)定,與工程現(xiàn)場結(jié)果一致。

    猜你喜歡
    煤柱撓度塑性
    基于應(yīng)變梯度的微尺度金屬塑性行為研究
    硬脆材料的塑性域加工
    Spontaneous multivessel coronary artery spasm diagnosed with intravascular ultrasound imaging:A case report
    采動(dòng)影響下雙巷掘進(jìn)煤柱承載特征研究
    鈹材料塑性域加工可行性研究
    石英玻璃的熱輔助高效塑性域干磨削
    懸高測量在橋梁撓度快速檢測中的應(yīng)用
    建筑物下煤柱巷道穿采分析
    河南科技(2014年11期)2014-02-27 14:16:58
    保護(hù)煤柱寬度的理論值分析
    河南科技(2014年5期)2014-02-27 14:08:25
    收縮徐變在不同鋪裝時(shí)間下對連續(xù)梁橋長期撓度的影響
    韩国av在线不卡| 人妻人人澡人人爽人人| 日本欧美视频一区| 国产高清三级在线| 国国产精品蜜臀av免费| 一级毛片 在线播放| 97超视频在线观看视频| 天堂8中文在线网| 午夜福利视频精品| 亚洲高清免费不卡视频| 国产精品人妻久久久影院| 少妇被粗大猛烈的视频| 免费高清在线观看视频在线观看| 国产午夜精品久久久久久一区二区三区| freevideosex欧美| 美女主播在线视频| 国产精品秋霞免费鲁丝片| 免费播放大片免费观看视频在线观看| 日韩制服骚丝袜av| 亚洲综合色网址| 亚洲情色 制服丝袜| 国产精品久久久久久av不卡| 一本一本综合久久| 韩国av在线不卡| 亚洲精品国产av蜜桃| 丝瓜视频免费看黄片| 欧美精品高潮呻吟av久久| 久久久久久久国产电影| 特大巨黑吊av在线直播| 国产欧美日韩一区二区三区在线 | 亚洲av二区三区四区| 久久人人爽人人片av| 成人18禁高潮啪啪吃奶动态图 | 下体分泌物呈黄色| 亚洲激情五月婷婷啪啪| 人妻夜夜爽99麻豆av| 国产免费一级a男人的天堂| 内地一区二区视频在线| 亚洲欧洲日产国产| 日韩av免费高清视频| 久久ye,这里只有精品| 久久97久久精品| 国产免费一区二区三区四区乱码| 在现免费观看毛片| 亚洲av免费高清在线观看| 人妻少妇偷人精品九色| 国产精品久久久久久精品电影小说| 99久久精品国产国产毛片| videossex国产| 美女内射精品一级片tv| 亚洲人与动物交配视频| 91精品国产九色| 99热6这里只有精品| 久久久久国产网址| 国产黄片视频在线免费观看| 国产熟女午夜一区二区三区 | 一区二区三区四区激情视频| 插阴视频在线观看视频| 日本欧美国产在线视频| 美女cb高潮喷水在线观看| 日韩亚洲欧美综合| 亚洲人成网站在线播| 国产日韩欧美在线精品| 精品一区二区三卡| 亚洲丝袜综合中文字幕| 91精品一卡2卡3卡4卡| 飞空精品影院首页| 亚洲精品亚洲一区二区| 只有这里有精品99| 日本vs欧美在线观看视频| 97在线视频观看| 亚洲国产av新网站| 成人免费观看视频高清| 麻豆成人av视频| 精品人妻在线不人妻| 两个人免费观看高清视频| 欧美人与善性xxx| 免费看不卡的av| 亚洲久久久国产精品| 麻豆成人av视频| 亚洲人成77777在线视频| 亚洲精品亚洲一区二区| 99热国产这里只有精品6| 亚洲av.av天堂| 久久人人爽人人爽人人片va| 成人漫画全彩无遮挡| 狂野欧美激情性xxxx在线观看| 精品久久久精品久久久| 欧美xxxx性猛交bbbb| 亚洲不卡免费看| 欧美+日韩+精品| 亚洲成人av在线免费| 大话2 男鬼变身卡| 2022亚洲国产成人精品| 国产成人免费观看mmmm| 制服人妻中文乱码| 亚洲高清免费不卡视频| 飞空精品影院首页| 国产成人午夜福利电影在线观看| 国产精品一区二区在线观看99| 日韩一本色道免费dvd| 日韩三级伦理在线观看| 日本猛色少妇xxxxx猛交久久| 在线观看国产h片| 好男人视频免费观看在线| 日韩欧美精品免费久久| 国产精品99久久99久久久不卡 | 特大巨黑吊av在线直播| tube8黄色片| 欧美激情国产日韩精品一区| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久久电影| 久久99精品国语久久久| 色婷婷av一区二区三区视频| 久久久久久伊人网av| 亚洲情色 制服丝袜| 狂野欧美白嫩少妇大欣赏| 丝袜喷水一区| 伦精品一区二区三区| 22中文网久久字幕| 少妇人妻久久综合中文| 成人漫画全彩无遮挡| 成年人午夜在线观看视频| 三上悠亚av全集在线观看| 少妇精品久久久久久久| 蜜桃在线观看..| 精品国产露脸久久av麻豆| 国产女主播在线喷水免费视频网站| 久久人人爽人人爽人人片va| 如日韩欧美国产精品一区二区三区 | 免费黄色在线免费观看| 特大巨黑吊av在线直播| 少妇被粗大的猛进出69影院 | 国产日韩欧美在线精品| 99久久人妻综合| 热99久久久久精品小说推荐| 免费观看在线日韩| 女人精品久久久久毛片| 午夜91福利影院| 国产精品国产三级国产专区5o| 国产亚洲av片在线观看秒播厂| 国产在线视频一区二区| 日日摸夜夜添夜夜添av毛片| 女性被躁到高潮视频| 久久久久精品性色| 老熟女久久久| 久久久精品区二区三区| xxxhd国产人妻xxx| 午夜福利视频精品| 久久 成人 亚洲| 视频在线观看一区二区三区| 久久免费观看电影| 观看av在线不卡| 久久久欧美国产精品| 91精品三级在线观看| 日本黄大片高清| 欧美日韩视频精品一区| 女人精品久久久久毛片| 亚州av有码| 中国三级夫妇交换| 亚洲欧洲日产国产| 国国产精品蜜臀av免费| 亚洲欧美清纯卡通| 蜜桃在线观看..| 日韩精品有码人妻一区| 少妇 在线观看| 国产在线视频一区二区| 亚洲欧洲日产国产| 香蕉精品网在线| 国产有黄有色有爽视频| 天美传媒精品一区二区| 日韩电影二区| 一区在线观看完整版| 人妻人人澡人人爽人人| 国产成人一区二区在线| 亚洲国产欧美在线一区| 午夜精品国产一区二区电影| 天天躁夜夜躁狠狠久久av| 午夜激情久久久久久久| 黄色怎么调成土黄色| 国产色婷婷99| 久久久久久久亚洲中文字幕| 一级爰片在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产色婷婷99| 精品酒店卫生间| 超色免费av| 91久久精品国产一区二区成人| 国产色爽女视频免费观看| 蜜臀久久99精品久久宅男| 国精品久久久久久国模美| 日韩视频在线欧美| 日本午夜av视频| 免费观看的影片在线观看| h视频一区二区三区| 美女主播在线视频| 国产成人aa在线观看| 伦理电影大哥的女人| 少妇的逼好多水| 999精品在线视频| 熟女电影av网| 亚洲色图综合在线观看| 午夜免费观看性视频| 婷婷色综合大香蕉| 最近的中文字幕免费完整| 国产白丝娇喘喷水9色精品| 国产精品免费大片| 中文字幕人妻熟人妻熟丝袜美| 午夜日本视频在线| av不卡在线播放| 伦精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 韩国av在线不卡| 欧美日韩综合久久久久久| 狂野欧美激情性xxxx在线观看| 久久精品国产鲁丝片午夜精品| 日韩一本色道免费dvd| 美女福利国产在线| 国产成人一区二区在线| 亚洲欧洲日产国产| 国产欧美日韩一区二区三区在线 | 天美传媒精品一区二区| av卡一久久| 欧美 亚洲 国产 日韩一| 51国产日韩欧美| 蜜桃久久精品国产亚洲av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美日韩卡通动漫| 考比视频在线观看| 亚洲精品乱久久久久久| 国产爽快片一区二区三区| 成人毛片a级毛片在线播放| 亚洲色图综合在线观看| av专区在线播放| 亚洲国产av影院在线观看| 97精品久久久久久久久久精品| 国产免费现黄频在线看| 国产精品女同一区二区软件| 亚洲国产毛片av蜜桃av| 亚洲精品日韩av片在线观看| 日韩一区二区三区影片| 久久久久久久久久久久大奶| 18禁在线无遮挡免费观看视频| 蜜桃在线观看..| 一本色道久久久久久精品综合| 另类精品久久| 亚洲精品自拍成人| 国产精品久久久久成人av| 欧美人与善性xxx| 亚洲内射少妇av| 亚洲av在线观看美女高潮| 国产视频内射| 成人毛片60女人毛片免费| 一二三四中文在线观看免费高清| 亚洲国产精品999| 精品午夜福利在线看| 中文字幕制服av| 天天躁夜夜躁狠狠久久av| 有码 亚洲区| 最后的刺客免费高清国语| 最近的中文字幕免费完整| 97精品久久久久久久久久精品| 日韩视频在线欧美| 久久狼人影院| 看十八女毛片水多多多| 亚洲av成人精品一二三区| 午夜老司机福利剧场| 国产免费又黄又爽又色| 美女大奶头黄色视频| 久久午夜综合久久蜜桃| 久热久热在线精品观看| 亚洲欧美一区二区三区国产| 国产精品女同一区二区软件| 精品亚洲乱码少妇综合久久| 2022亚洲国产成人精品| 久久午夜福利片| 亚洲精品久久成人aⅴ小说 | 多毛熟女@视频| 国产视频内射| 在线精品无人区一区二区三| 亚洲国产av影院在线观看| 成人无遮挡网站| 新久久久久国产一级毛片| 亚洲av二区三区四区| 国产精品人妻久久久影院| 大香蕉97超碰在线| 国产亚洲一区二区精品| 一区二区三区乱码不卡18| 国产爽快片一区二区三区| 色网站视频免费| 亚洲一级一片aⅴ在线观看| 国内精品宾馆在线| 亚洲丝袜综合中文字幕| 国产极品天堂在线| 校园人妻丝袜中文字幕| 九九爱精品视频在线观看| 亚洲欧美精品自产自拍| 日韩免费高清中文字幕av| 伦理电影大哥的女人| 欧美日韩视频高清一区二区三区二| 日韩一区二区视频免费看| 一本—道久久a久久精品蜜桃钙片| 亚洲人成网站在线播| 如日韩欧美国产精品一区二区三区 | 大片免费播放器 马上看| 欧美日本中文国产一区发布| 国产黄片视频在线免费观看| 国产在线视频一区二区| 91在线精品国自产拍蜜月| 日韩成人伦理影院| 久久影院123| 一级毛片aaaaaa免费看小| 亚洲精品色激情综合| av线在线观看网站| 国产一区二区在线观看av| 国产欧美日韩一区二区三区在线 | 丁香六月天网| 亚洲国产av新网站| 十八禁高潮呻吟视频| 成人国产麻豆网| 久热久热在线精品观看| 国产熟女午夜一区二区三区 | 七月丁香在线播放| 中文字幕人妻熟人妻熟丝袜美| 免费播放大片免费观看视频在线观看| 蜜桃久久精品国产亚洲av| 欧美日韩精品成人综合77777| 日本黄色片子视频| 日本与韩国留学比较| 在线免费观看不下载黄p国产| 大话2 男鬼变身卡| 美女国产高潮福利片在线看| 国产免费一级a男人的天堂| 高清欧美精品videossex| 人体艺术视频欧美日本| 中文字幕精品免费在线观看视频 | 亚洲精品乱久久久久久| 亚洲第一av免费看| 一边亲一边摸免费视频| 日韩一区二区视频免费看| 国产国拍精品亚洲av在线观看| 欧美精品高潮呻吟av久久| 十分钟在线观看高清视频www| 国产乱人偷精品视频| 国产欧美日韩综合在线一区二区| 三上悠亚av全集在线观看| 日本色播在线视频| 女的被弄到高潮叫床怎么办| av福利片在线| 欧美成人精品欧美一级黄| 亚洲美女黄色视频免费看| 久久99热6这里只有精品| 免费高清在线观看视频在线观看| 性色av一级| 日产精品乱码卡一卡2卡三| 最近2019中文字幕mv第一页| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩另类电影网站| 久久99蜜桃精品久久| 久久99热6这里只有精品| 免费人妻精品一区二区三区视频| 欧美3d第一页| 少妇丰满av| 成年av动漫网址| 国产国语露脸激情在线看| 亚洲av不卡在线观看| 日韩伦理黄色片| 午夜精品国产一区二区电影| 菩萨蛮人人尽说江南好唐韦庄| 九九爱精品视频在线观看| 搡女人真爽免费视频火全软件| 18禁观看日本| 久久久久久久久大av| 高清在线视频一区二区三区| 十分钟在线观看高清视频www| 日本91视频免费播放| 七月丁香在线播放| 欧美一级a爱片免费观看看| 亚洲综合色网址| 精品视频人人做人人爽| 伦理电影免费视频| 精品少妇黑人巨大在线播放| 国精品久久久久久国模美| 五月开心婷婷网| 国产精品久久久久久av不卡| 日本免费在线观看一区| 美女国产高潮福利片在线看| freevideosex欧美| 91精品国产九色| 国产一区二区在线观看日韩| 久久久久久久久久久丰满| 亚洲av.av天堂| 亚洲国产毛片av蜜桃av| √禁漫天堂资源中文www| av专区在线播放| 又粗又硬又长又爽又黄的视频| 国产精品蜜桃在线观看| 在线精品无人区一区二区三| 99久久精品一区二区三区| 夜夜看夜夜爽夜夜摸| 国产国拍精品亚洲av在线观看| av线在线观看网站| 国产国语露脸激情在线看| 国产69精品久久久久777片| 熟女电影av网| 精品少妇内射三级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品久久久久成人av| 成人免费观看视频高清| 久久人人爽av亚洲精品天堂| 亚洲丝袜综合中文字幕| 色吧在线观看| 久久99热6这里只有精品| 国产亚洲一区二区精品| 看免费成人av毛片| 久久久久人妻精品一区果冻| 久久国产精品男人的天堂亚洲 | 国产欧美另类精品又又久久亚洲欧美| 18禁观看日本| 91国产中文字幕| 亚洲人与动物交配视频| 国产免费福利视频在线观看| 男男h啪啪无遮挡| 国产探花极品一区二区| 欧美另类一区| 精品久久久噜噜| 久久久久久久久大av| 精品99又大又爽又粗少妇毛片| 国产极品天堂在线| 午夜激情av网站| 日本黄大片高清| 最后的刺客免费高清国语| 久久久久国产网址| 9色porny在线观看| av视频免费观看在线观看| 精品久久蜜臀av无| 新久久久久国产一级毛片| 视频区图区小说| 久久精品久久久久久久性| 亚洲av日韩在线播放| 欧美日韩国产mv在线观看视频| 久久久久久久久久成人| 国产片内射在线| 国产成人精品一,二区| 精品久久久久久久久亚洲| 久久精品国产亚洲网站| 中文字幕人妻丝袜制服| 亚洲精品久久午夜乱码| 一区二区日韩欧美中文字幕 | 一个人免费看片子| 日韩av在线免费看完整版不卡| 精品少妇久久久久久888优播| 国产又色又爽无遮挡免| 久久亚洲国产成人精品v| 午夜久久久在线观看| 欧美精品高潮呻吟av久久| 久久久久久久久大av| 自线自在国产av| 国产精品成人在线| 少妇人妻精品综合一区二区| 交换朋友夫妻互换小说| 亚洲综合精品二区| 国产成人精品在线电影| 啦啦啦中文免费视频观看日本| 91精品一卡2卡3卡4卡| 亚洲欧洲日产国产| 亚洲无线观看免费| 国模一区二区三区四区视频| 精品人妻偷拍中文字幕| 欧美 日韩 精品 国产| 七月丁香在线播放| 免费看不卡的av| 日日摸夜夜添夜夜爱| 一级毛片黄色毛片免费观看视频| 亚洲五月色婷婷综合| 看非洲黑人一级黄片| 午夜福利视频在线观看免费| 一级毛片电影观看| 国产精品.久久久| 免费日韩欧美在线观看| 中文欧美无线码| 亚洲少妇的诱惑av| 国产片内射在线| 亚洲中文av在线| 在现免费观看毛片| 男人操女人黄网站| 亚洲av成人精品一区久久| 国产精品一区二区在线不卡| 色94色欧美一区二区| 在现免费观看毛片| 午夜av观看不卡| 亚洲国产精品国产精品| 男女国产视频网站| 国产成人精品在线电影| 国产黄色视频一区二区在线观看| 日本午夜av视频| 少妇精品久久久久久久| 久久综合国产亚洲精品| 亚洲成人av在线免费| 欧美xxxx性猛交bbbb| 国产欧美日韩综合在线一区二区| 青春草视频在线免费观看| 一级二级三级毛片免费看| 亚洲人成77777在线视频| 午夜激情av网站| 精品久久久久久久久亚洲| 国产日韩欧美在线精品| 日韩成人av中文字幕在线观看| 国产熟女欧美一区二区| 亚洲精品成人av观看孕妇| 久久精品国产自在天天线| 精品久久久噜噜| 国产成人av激情在线播放 | 一级二级三级毛片免费看| 成人亚洲欧美一区二区av| 97精品久久久久久久久久精品| 免费播放大片免费观看视频在线观看| 一级毛片aaaaaa免费看小| 久久97久久精品| 成人毛片60女人毛片免费| 亚洲精品中文字幕在线视频| 中国三级夫妇交换| 亚洲性久久影院| 大片免费播放器 马上看| 一级,二级,三级黄色视频| 中文字幕免费在线视频6| 国产男女内射视频| 国产精品.久久久| 黄色视频在线播放观看不卡| 99热国产这里只有精品6| 成年女人在线观看亚洲视频| 天堂8中文在线网| 99国产综合亚洲精品| 少妇人妻精品综合一区二区| www.色视频.com| 国产午夜精品一二区理论片| 亚洲高清免费不卡视频| 国产黄色免费在线视频| 乱码一卡2卡4卡精品| 18禁在线播放成人免费| 亚洲第一区二区三区不卡| 婷婷成人精品国产| xxx大片免费视频| 亚洲精品美女久久av网站| 国产老妇伦熟女老妇高清| 亚洲熟女精品中文字幕| 久久狼人影院| 国产免费视频播放在线视频| 久久久亚洲精品成人影院| 99久久人妻综合| 搡老乐熟女国产| 国产av码专区亚洲av| 欧美丝袜亚洲另类| 免费不卡的大黄色大毛片视频在线观看| 伦精品一区二区三区| 日本欧美国产在线视频| 看非洲黑人一级黄片| 精品人妻在线不人妻| 王馨瑶露胸无遮挡在线观看| 亚洲av国产av综合av卡| videos熟女内射| 国产精品女同一区二区软件| 寂寞人妻少妇视频99o| 亚洲国产欧美日韩在线播放| 日本猛色少妇xxxxx猛交久久| 成人亚洲精品一区在线观看| 少妇 在线观看| 99久久中文字幕三级久久日本| 国产精品99久久久久久久久| a级毛片黄视频| 亚洲美女视频黄频| 尾随美女入室| 在线 av 中文字幕| 搡女人真爽免费视频火全软件| 亚洲精华国产精华液的使用体验| 夜夜骑夜夜射夜夜干| 成人黄色视频免费在线看| 欧美一级a爱片免费观看看| 免费观看无遮挡的男女| 男女啪啪激烈高潮av片| 在线观看人妻少妇| 啦啦啦中文免费视频观看日本| 久久热精品热| 精品亚洲成国产av| 天美传媒精品一区二区| 母亲3免费完整高清在线观看 | 成人无遮挡网站| 亚洲人成网站在线观看播放| 精品卡一卡二卡四卡免费| 日产精品乱码卡一卡2卡三| 国产高清三级在线| 91在线精品国自产拍蜜月| 一级,二级,三级黄色视频| 精品亚洲成国产av| 国产av国产精品国产| 这个男人来自地球电影免费观看 | 18禁在线播放成人免费| 高清黄色对白视频在线免费看| 九九爱精品视频在线观看| 成人漫画全彩无遮挡| 一级毛片黄色毛片免费观看视频| 丝袜在线中文字幕| 一级毛片电影观看| 精品熟女少妇av免费看| 狠狠精品人妻久久久久久综合| 制服人妻中文乱码| 国产免费一级a男人的天堂| 国产在线一区二区三区精| 人人澡人人妻人| 精品人妻偷拍中文字幕| 久久久a久久爽久久v久久| 国产成人免费无遮挡视频| 这个男人来自地球电影免费观看 | 黄色怎么调成土黄色| 久久久亚洲精品成人影院|