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

    基于MATLAB和FLAC3D的開采沉陷仿真教學實驗系統(tǒng)研發(fā)

    2021-04-07 10:11:22滕超群江克貴
    金屬礦山 2021年3期
    關(guān)鍵詞:巖層機理可視化

    王 磊 滕超群 江克貴

    (安徽理工大學空間信息與測繪工程學院,安徽淮南232001)

    《礦山開采沉陷學》是煤炭行業(yè)特色高校測繪工程專業(yè)的必修課程,也是數(shù)字礦山與沉陷控制學科的專業(yè)基礎(chǔ)課程,課程包括變形監(jiān)測、開采沉陷規(guī)律、開采沉陷災(zāi)害預(yù)測和防控以及礦山環(huán)境修復(fù)和治理等內(nèi)容,涉及地質(zhì)、采礦、測量、環(huán)境、力學、建筑、計算機等多個學科,屬于典型的交叉學科課程。由于課程涵蓋的知識理論體系龐雜,并且存在部分概念抽象、機理與規(guī)律復(fù)雜、理論與方法深奧等授課障礙,這就要求學生除了掌握多學科交叉理論之外,還應(yīng)具備較高的空間想象和邏輯思維能力。上述特點使得大部分學生感到課程理論枯燥、生澀難懂、失去興趣,導(dǎo)致學生對基本概念和理論認識不清,機理與規(guī)律的科學內(nèi)涵理解膚淺,進而難以達到培養(yǎng)學生解決復(fù)雜工程問題能力的目標,創(chuàng)新性應(yīng)用更是無從談起。在此背景下,探索適合課堂教學規(guī)律和課程特點的新教學方法,是當下《礦山開采沉陷學》教學亟需解決的問題。仿真教學是一種具有綜合作用的教學手段,使學生置身于仿真環(huán)境中,可以充分調(diào)動感覺、運動和思維[1-2]。有關(guān)研究表明[3],在仿真教學模式下,學生可以記憶約70%的內(nèi)容,而傳統(tǒng)的“教師講、學生聽”的教學模式,學生只能記憶約30%的內(nèi)容。仿真教學方法將課程的重點、難點內(nèi)容基于模擬、虛擬等相關(guān)技術(shù)制作成仿真教學內(nèi)容,通過教師課堂的引導(dǎo)仿真訓練和課后學生的自主仿真訓練,將會極大地提升學生對課程的學習效果。

    近些年來,MATLAB和FLAC3D軟件在科學研究和工程應(yīng)用領(lǐng)域得到了廣泛應(yīng)用,其強大的計算、數(shù)據(jù)分析、非線性模擬求解、圖形可視化等功能克服了物理實驗觀測困難和重復(fù)性差等缺陷。文獻[4-8]將FLAC3D軟件運用到開采沉陷領(lǐng)域,分別對高速公路下伏采空區(qū)殘余沉降、矸石充填開采沉陷、復(fù)雜開挖過程力學仿真、開采沉陷預(yù)計、似膏體充填開采沉陷開展了數(shù)值模擬研究;文獻[9-13]基于MATLAB軟件系統(tǒng)地開展了開采沉陷預(yù)計系統(tǒng)研發(fā)、開采沉陷預(yù)計參數(shù)求解、移動變形可視化等研究,取得了良好的工程應(yīng)用結(jié)果。綜合上述研究,為提高《礦山開采沉陷學》課程的教學質(zhì)量,文獻[2,14]基于MATLAB開展了開采沉陷預(yù)計控制仿真教學研究;文獻[15]利用FLAC3D進行了開采沉陷數(shù)值模擬教學研究,教學實踐表明,將計算機仿真教學納入礦山開采沉陷教學過程,可以極大地提高學生理解能力,培養(yǎng)學生學習興趣。然而,復(fù)雜三維開采沉陷模型構(gòu)建一直是FLAC3D軟件前處理中的難題[16],且上述教學系統(tǒng)未涉及到巖層移動機理、特殊地質(zhì)采礦和預(yù)計參數(shù)求解等關(guān)鍵技術(shù)的研發(fā)。

    為此,本研究針對《礦山開采沉陷學》中巖層移動與地表沉陷機理及規(guī)律的教學重點和難點,耦合MATLAB平臺和FLAC3D數(shù)值模擬技術(shù),研發(fā)礦山開采沉陷仿真教學實驗系統(tǒng)。通過該系統(tǒng)的實訓,可加深學生對開采沉陷機理和規(guī)律的理解和認識,有助于高效解決課程教學中的瓶頸問題。

    1 系統(tǒng)總體框架設(shè)計

    巖層移動與地表沉陷機理及規(guī)律是《礦山開采沉陷學》的核心基礎(chǔ)理論,其相關(guān)內(nèi)容復(fù)雜且抽象,是教學的重點和難點,這也是導(dǎo)致常規(guī)方法教學效果不理想的主要原因。為此,本研究設(shè)計并開發(fā)礦山開采沉陷仿真教學實驗系統(tǒng),為廣大師生提供交互式圖形用戶界面,幫助學生更深刻地理解開采沉陷學中有關(guān)內(nèi)容。該系統(tǒng)主要分為兩部分,一是揭示有關(guān)巖層移動機理及規(guī)律仿真,二是地表沉陷規(guī)律仿真。系統(tǒng)總體框架如圖1所示。

    針對巖層移動基本理論和相關(guān)教學內(nèi)容,本研究通過MATLAB和FLAC3D數(shù)值模擬技術(shù)耦合實現(xiàn)對巖層移動機理及規(guī)律過程的仿真,揭示常規(guī)開采、特殊開采(充填開采、部分開采等)巖層移動機理以及覆巖應(yīng)力場、裂隙場、位移場的時空演化規(guī)律。對于地表沉陷基本理論的相關(guān)教學內(nèi)容,本研究基于MATLAB數(shù)值計算平臺開發(fā)地表沉陷規(guī)律仿真實驗系統(tǒng),從開采沉陷規(guī)律、開采沉陷預(yù)計原理、開采沉陷預(yù)計參數(shù)反演3個方面開發(fā)用于仿真教學研究的可視化模塊[7]。共開發(fā)了二維動態(tài)預(yù)計仿真、三維動態(tài)預(yù)計仿真、地質(zhì)采礦因素影響規(guī)律仿真、預(yù)計原理仿真、任意工作面剖分及預(yù)計仿真、預(yù)計參數(shù)智能反演仿真6個可視化仿真子模塊,各個模塊具有一定的可擴展性,與地表沉陷規(guī)律相關(guān)的其他課程內(nèi)容可在本系統(tǒng)的基礎(chǔ)上擴展。系統(tǒng)主界面如圖2所示。

    2 仿真教學模塊研發(fā)及應(yīng)用

    2.1 基于MATLAB和FLAC3D耦合的巖層移動規(guī)律仿真

    開采沉陷的實質(zhì)為復(fù)雜的巖層時空力學變形過程,數(shù)值模擬的本質(zhì)是對研究對象進行力學建模、求解和可視化分析的過程,因此用FLAC3D數(shù)值模擬方法仿真開采沉陷過程從機理上是契合的[6]。FLAC3D軟件具有豐富的本構(gòu)模型庫和結(jié)構(gòu)單元庫、多樣化的邊界條件、內(nèi)置復(fù)雜模型生成器和fish程序編譯器以及強大的可視化分析等系列模塊及功能,可實現(xiàn)對常規(guī)開采、特殊開采(充填開采、部分開采等)巖層移動機理及規(guī)律的精細模擬和仿真。然而當面對復(fù)雜地層和地質(zhì)構(gòu)造時,復(fù)雜三維地質(zhì)體模型建立一直是FLAC3D軟件前處理中的難題,而且該款軟件圖形后處理功能匱乏,無法有效滿足實際工程應(yīng)用要求。

    考慮到MATLAB具備數(shù)值計算與可視化表達優(yōu)勢和FLAC3D擅長數(shù)值模擬的特點,本研究提出了基于MATLAB和FLAC3D耦合的巖層移動規(guī)律仿真設(shè)計方法。主要思路為借助MATLAB軟件強大的矩陣功能,通過編程計算并形成FLAC3D命令語言,實現(xiàn)復(fù)雜三維地質(zhì)體建模(前處理),F(xiàn)LAC3D模型運行結(jié)果分析時,提取FLAC3D模型巖層移動應(yīng)力、位移等成果數(shù)據(jù),利用MATLAB強大的數(shù)據(jù)處理和可視化功能(后處理)進行規(guī)律和機理分析?;贛ATLAB和FLAC3D耦合的巖層移動規(guī)律仿真流程如圖3所示。

    在準備相關(guān)仿真教學內(nèi)容時,考慮到耦合求解過程耗時長、計算環(huán)境配備要求高,可課前準備典型地質(zhì)采礦條件下開采沉陷數(shù)值模擬試驗結(jié)果,然后通過錄屏、視頻制作、PPT等方式提前準備相關(guān)仿真訓練文件。通過MATLAB和FLAC3D耦合數(shù)值模擬仿真訓練不僅可以加深學生對巖層理論的理解,還可以提升利用MATLAB和FLAC3D耦合解決復(fù)雜工程問題的創(chuàng)新能力。利用MATLAB和FLAC3D耦合數(shù)值模擬仿真某工作面開采過程中巖層與地表位移場及覆巖應(yīng)力場演化特征的部分成果如圖4所示。

    2.2 開采沉陷預(yù)計原理仿真

    概率積分法是我國礦山開采沉陷中應(yīng)用較為成熟和廣泛的預(yù)計方法之一。該方法基于隨機介質(zhì)理論,用概率積分表示地下單元開采引起地表移動和變形的預(yù)計公式,結(jié)合疊加原理可計算出整個開采引起的地表移動和變形[17]。地表移動盆地主斷面的半無限開采和有限開采預(yù)計是概率積分法預(yù)計的基礎(chǔ)。半無限開采指在平面問題中,開切眼一側(cè)的煤層未被開采,而其他方向的煤層已全部采出,即沿垂直于工作面推進方向的開采尺寸足夠大,達到充分采動;有限開采是根據(jù)半無限疊加原理,認為等效于兩個半無限開采之差。由于預(yù)計公式復(fù)雜抽象,學生很難理解計算原理,故本系統(tǒng)設(shè)計了矩形工作面下的開采沉陷預(yù)計原理仿真可視化模塊,通過對兩種預(yù)計類型下的基本原理及其地表任意點的移動和變形進行演示,有助于加深學生對開采沉陷預(yù)計原理的理解。

    本模塊的界面如圖5所示,主要由預(yù)計類型和任意點預(yù)計兩個面板組成。在預(yù)計類型面板中,有半無限開采和有限開采兩個按鈕,分別對兩種情況下地表主斷面的下沉形態(tài)進行演示。由圖5可看出,有限開采等效于兩個半無限開采疊加。在任意點預(yù)計面板中,設(shè)計5個單選按鈕,分別是地表下沉、傾斜、曲率、水平變形和水平移動5種移動變形,同時模擬煤層開采單元B(s,t)處引起地表任意點A(x,y)的移動變形規(guī)律,借此化繁為簡,有助于加深學生對預(yù)計公式的理解。

    近年來,矩形工作面開采沉陷移動與變形預(yù)計方法已經(jīng)相當成熟,但由于實際采礦條件的限制,開采工作面會出現(xiàn)各種不規(guī)則形狀,而采用常規(guī)矩形工作面開采沉陷預(yù)計方法對不規(guī)則工作面進行預(yù)計效果不理想。故本研究設(shè)計了任意形狀工作面剖分仿真實驗?zāi)K,針對復(fù)雜不規(guī)則形狀工作面,應(yīng)用三角剖分和二重積分預(yù)計相結(jié)合的方法著重解決其預(yù)計困難的問題,具體原理可參考文獻[18-20]。仿真模塊界面及運算結(jié)果如圖6所示。

    由圖6可知:該可視化模塊有工作面剖分和開采沉陷預(yù)計兩個面板。工作面剖分面板中列舉了三角剖分和矩形剖分兩種方法。系統(tǒng)以三角剖分算法為基礎(chǔ),將不規(guī)則工作面剖分成多個三角形,并將剖分后的工作面顯示在右側(cè)圖形顯示框中,再對剖分后的工作面進行開采沉陷預(yù)計。

    2.3 開采沉陷預(yù)計參數(shù)智能反演仿真

    概率積分法關(guān)于地表移動盆地主要有8個預(yù)計參數(shù),分別為下沉系數(shù)q,主要影響角正切tanβ,開采影響傳播角θ,水平移動系數(shù)b,矩形工作面上、下、左、右拐點偏移距S1、S2、S3、S4。下沉系數(shù)指充分采動條件下,地表最大下沉值W0與煤層法線采厚m在鉛垂方向投影長度的比值,其預(yù)測的準確程度直接影響到地表移動變形及采動影響區(qū)建筑物破壞程度的精確預(yù)計。主要影響角正切為開采深度與主要影響半徑之比。開采影響傳播角是傾向主斷面特有的預(yù)計參數(shù),是指在移動盆地傾向主斷面上,按拐點偏移距求得的計算開采邊界和地表下沉曲線拐點的連線與水平線在下山方向的夾角。水平移動系數(shù)指充分采動條件下,單一煤層開采最大水平移動與最大下沉之比。

    預(yù)計參數(shù)由實測數(shù)據(jù)反演得來,概率積分函數(shù)模型具有高度非線性、參數(shù)眾多并具有相關(guān)性,基于何種方法精確反演預(yù)計參數(shù)則成為重點和難點。本系統(tǒng)將目前主流的智能算法[21-25]應(yīng)用于概率積分函數(shù)中,構(gòu)建基于智能算法的概率積分預(yù)計參數(shù)反演模型。系統(tǒng)界面如圖7所示。

    在智能算法參數(shù)反演仿真實驗可視化模塊中,本研究設(shè)計了智能算法、預(yù)計參數(shù)輸出兩個面板。在智能算法面板中,共選取了遺傳算法、蟻群算法、蛙跳算法、狼群算法、煙花算法5種智能算法,以實測下沉值和水平移動值分別與它們的預(yù)計值之差的平方和作為適應(yīng)度函數(shù)[9],在界面右側(cè)添加兩個圖形顯示區(qū)域,以便學生更直觀地理解算法運行時每一次反演預(yù)計值與實測值的擬合情況。此次實驗是以淮南市顧橋礦某工作面為背景展開,為了定量分析智能算法反演的精確程度,本模塊可將反演的最后結(jié)果輸出在界面上,對于實際工程應(yīng)用有一定的參考價值。

    2.4 開采沉陷規(guī)律仿真

    為了仿真不同地質(zhì)采礦條件下地表沉陷的時空演化規(guī)律,基于MATLAB開發(fā)了多工作面、不同頂板管理方法等條件下的地表沉陷任意點動態(tài)預(yù)計模塊,可滿足地表任意點、主斷面、全盆地的時空演化規(guī)律仿真演示需求,也可實現(xiàn)各因素對地表沉陷影響規(guī)律的形象仿真。地表走向主斷面動態(tài)預(yù)計和全盆地動態(tài)演化模塊界面如圖8和圖9所示。

    在開采沉陷動態(tài)預(yù)計仿真實驗和三維動態(tài)預(yù)計仿真實驗?zāi)K中,本研究利用MATLAB的二維、三維繪圖及動畫功能,實現(xiàn)了工作面在不同采動狀態(tài)下,地表下沉、傾斜、曲率、水平移動和水平變形在主斷面和全盆地的動態(tài)演化,使學生比較形象、生動地掌握礦山開采在經(jīng)歷啟動、非充分采動、充分采動和超充分采動過程中地表移動盆地主斷面、全盆地移動與變形的時空演化基本規(guī)律和基本概念。

    開采沉陷動態(tài)預(yù)計仿真實驗?zāi)K界面由顯示方向面板、二維顯示面板以及圖形顯示框組成。學生可以自行選擇地表移動盆地主斷面的走向和傾向。在二維顯示面板中,模擬工作面沿走向方向推進,每次推進100 m,直至達到工作面走向長度800 m,共推進8次,結(jié)合MATLAB動態(tài)繪圖功能在圖形顯示框中繪制出沿走向(傾向)方向的動態(tài)預(yù)計曲線,同時,圖形顯示框的下方繪制出模擬工作面推進時煤層的采動程度,以此讓學生清晰地了解煤層開采過程中地表移動盆地的變形。

    因地表移動盆地三維變形涉及到方位角度數(shù),故三維動態(tài)預(yù)計仿真實驗?zāi)K添加了顯示方向面板,在確定方位角度數(shù)后,以上述模塊中的模擬工作面為背景,在圖形顯示框演示地表移動盆地三維動態(tài)形變場(圖9),可在此基礎(chǔ)上添加等值線和工作面推進程度等信息輔助學生理解。

    地質(zhì)和采礦因素對開采沉陷分布規(guī)律有一定的影響,為加深學生對不同地質(zhì)采礦條件下開采沉陷規(guī)律的理解,開發(fā)出地質(zhì)采礦因素影響規(guī)律仿真模塊,研究在不同的地質(zhì)和采礦條件下,地表主斷面沉陷規(guī)律及三維下沉場、水平移動場分布規(guī)律。模塊設(shè)計界面如圖10所示。

    模塊仿真設(shè)計時,在圖10所示的模塊界面左側(cè)設(shè)計覆巖巖性,松散層厚度,工作面幾何特征參數(shù)(采深度、采厚、走向長度、傾向斜長),頂板管理方法等開采沉陷影響因素可視化控件,在其右側(cè)設(shè)計不同地質(zhì)采礦條件下地表移動盆地主斷面和全盆地下沉的定量可視化結(jié)果。仿真訓練時,當改變界面左側(cè)的不同地質(zhì)采礦條件時,右側(cè)可實時、定量地顯示主斷面和全盆地的下沉空間分布,可有效實現(xiàn)各因素對地表沉陷影響規(guī)律的定量仿真。

    3 結(jié) 語

    《礦山開采沉陷學》是一門典型的交叉學科課程,其涵蓋的知識理論體系龐雜,存在部分概念抽象、機理與規(guī)律復(fù)雜、理論與方法深奧等授課障礙,導(dǎo)致常規(guī)教學方法效果不理想。仿真教學作為新式的教學手段,有助于克服常規(guī)方法的弊端,針對該課程的核心基礎(chǔ)理論、教學的重點和難點,即巖層移動與地表沉陷機理及規(guī)律,開發(fā)了礦山開采沉陷仿真教學實驗系統(tǒng)。該系統(tǒng)通過耦合MATLAB和FLAC3D研發(fā)了巖層移動機理及規(guī)律仿真模塊,融合概率積分法原理開發(fā)了基于MATLAB GUI地表沉陷規(guī)律仿真模塊,實現(xiàn)了對礦山開采從巖層到地表的移動變形機理和規(guī)律的仿真,通過該系統(tǒng)的實訓,有助于加深學生對開采沉陷機理和規(guī)律的理解和認識。該系統(tǒng)的開發(fā)與應(yīng)用,有助于解決《礦山開采沉陷學》課程教學的瓶頸問題,對于類似特點的交叉學科教學實踐也有一定的借鑒價值。

    猜你喜歡
    巖層機理可視化
    基于CiteSpace的足三里穴研究可視化分析
    高應(yīng)力巖層巷道鉆孔爆破卸壓技術(shù)
    基于Power BI的油田注水運行動態(tài)分析與可視化展示
    云南化工(2021年8期)2021-12-21 06:37:54
    隔熱纖維材料的隔熱機理及其應(yīng)用
    地球故事之復(fù)理石巖層
    基于CGAL和OpenGL的海底地形三維可視化
    煤層氣吸附-解吸機理再認識
    中國煤層氣(2019年2期)2019-08-27 00:59:30
    “融評”:黨媒評論的可視化創(chuàng)新
    傳媒評論(2019年4期)2019-07-13 05:49:14
    霧霾機理之問
    DNTF-CMDB推進劑的燃燒機理
    火炸藥學報(2014年1期)2014-03-20 13:17:27
    两性夫妻黄色片| 国产xxxxx性猛交| 精品国产亚洲在线| 一二三四在线观看免费中文在| 美国免费a级毛片| 欧美大码av| 桃花免费在线播放| 老司机午夜十八禁免费视频| 久久香蕉激情| 午夜老司机福利片| 狂野欧美激情性xxxx| 欧美人与性动交α欧美精品济南到| 18在线观看网站| 亚洲av电影在线进入| 美女高潮到喷水免费观看| 国产精品麻豆人妻色哟哟久久| 亚洲熟女毛片儿| 欧美久久黑人一区二区| 国产区一区二久久| av免费在线观看网站| 在线观看免费高清a一片| 中文字幕人妻熟女乱码| 国产成人精品无人区| 亚洲第一青青草原| 女人高潮潮喷娇喘18禁视频| 999精品在线视频| 一个人免费看片子| 亚洲 国产 在线| 麻豆乱淫一区二区| 久久精品亚洲av国产电影网| 久久久久精品国产欧美久久久| 女人精品久久久久毛片| 久久这里只有精品19| aaaaa片日本免费| 亚洲成国产人片在线观看| 亚洲免费av在线视频| 777米奇影视久久| 在线播放国产精品三级| 精品一区二区三区av网在线观看 | 岛国毛片在线播放| 亚洲色图av天堂| 久久久精品区二区三区| 国产免费av片在线观看野外av| 免费女性裸体啪啪无遮挡网站| 夜夜夜夜夜久久久久| 亚洲av片天天在线观看| 99久久人妻综合| av有码第一页| 国产97色在线日韩免费| 成年人午夜在线观看视频| 啦啦啦免费观看视频1| 狠狠婷婷综合久久久久久88av| 亚洲国产av新网站| 99久久99久久久精品蜜桃| 免费在线观看影片大全网站| 日韩三级视频一区二区三区| 午夜精品久久久久久毛片777| 人人妻人人添人人爽欧美一区卜| 成人av一区二区三区在线看| 人人妻人人澡人人爽人人夜夜| 亚洲一卡2卡3卡4卡5卡精品中文| 日日夜夜操网爽| 老司机靠b影院| 美女高潮喷水抽搐中文字幕| 久9热在线精品视频| 亚洲情色 制服丝袜| 国产精品亚洲一级av第二区| 亚洲精品美女久久久久99蜜臀| 欧美精品一区二区免费开放| 亚洲精品av麻豆狂野| 国产亚洲一区二区精品| 日本精品一区二区三区蜜桃| 精品国产乱码久久久久久小说| 久久99热这里只频精品6学生| 精品一区二区三区视频在线观看免费 | 纵有疾风起免费观看全集完整版| 国产成人啪精品午夜网站| 曰老女人黄片| 一本一本久久a久久精品综合妖精| 免费人妻精品一区二区三区视频| 成人18禁高潮啪啪吃奶动态图| 中文字幕精品免费在线观看视频| 12—13女人毛片做爰片一| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕一级| 国产精品自产拍在线观看55亚洲 | 老熟女久久久| 精品视频人人做人人爽| 国产一卡二卡三卡精品| 亚洲,欧美精品.| 一级毛片女人18水好多| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久欧美国产精品| 亚洲精品美女久久久久99蜜臀| 91成年电影在线观看| 中文欧美无线码| 人妻久久中文字幕网| 成人18禁高潮啪啪吃奶动态图| 精品久久久精品久久久| 欧美亚洲日本最大视频资源| 国产高清激情床上av| 最近最新中文字幕大全免费视频| 国产成人免费观看mmmm| 精品国产国语对白av| 欧美 日韩 精品 国产| 日韩欧美一区视频在线观看| 19禁男女啪啪无遮挡网站| 最近最新中文字幕大全免费视频| 亚洲情色 制服丝袜| 纯流量卡能插随身wifi吗| 国产精品.久久久| 精品国产一区二区久久| 欧美日韩精品网址| 亚洲精品美女久久久久99蜜臀| 一进一出抽搐动态| 欧美人与性动交α欧美软件| 国产不卡av网站在线观看| 在线观看免费日韩欧美大片| 精品少妇内射三级| 精品人妻在线不人妻| 青草久久国产| 黄色毛片三级朝国网站| 国产淫语在线视频| 欧美亚洲 丝袜 人妻 在线| 色视频在线一区二区三区| 999久久久国产精品视频| 俄罗斯特黄特色一大片| 男女高潮啪啪啪动态图| 12—13女人毛片做爰片一| 日韩视频一区二区在线观看| 美女扒开内裤让男人捅视频| 悠悠久久av| 岛国毛片在线播放| 亚洲精品国产一区二区精华液| 久久国产精品影院| 国产午夜精品久久久久久| 精品久久久久久电影网| 国产精品久久久久成人av| 国产高清激情床上av| 久久九九热精品免费| 精品高清国产在线一区| 在线永久观看黄色视频| 精品国产一区二区久久| 色婷婷久久久亚洲欧美| 一二三四社区在线视频社区8| 久久久精品国产亚洲av高清涩受| 少妇 在线观看| 99在线人妻在线中文字幕 | 一区二区三区国产精品乱码| 国产xxxxx性猛交| 超色免费av| 国精品久久久久久国模美| 在线观看免费视频日本深夜| 岛国毛片在线播放| √禁漫天堂资源中文www| 热re99久久精品国产66热6| 性高湖久久久久久久久免费观看| 国产亚洲午夜精品一区二区久久| 欧美+亚洲+日韩+国产| 天堂中文最新版在线下载| 国产高清国产精品国产三级| 三上悠亚av全集在线观看| 中文字幕人妻丝袜制服| 一级毛片电影观看| av线在线观看网站| 大香蕉久久成人网| 中文字幕精品免费在线观看视频| 久久久久精品人妻al黑| 国产亚洲欧美精品永久| 中文字幕制服av| 91av网站免费观看| 欧美日韩一级在线毛片| 国产日韩一区二区三区精品不卡| 不卡av一区二区三区| 欧美亚洲日本最大视频资源| 日韩欧美国产一区二区入口| 国产精品久久久久久精品古装| 亚洲精品国产区一区二| 无遮挡黄片免费观看| 十八禁网站免费在线| 如日韩欧美国产精品一区二区三区| 高清在线国产一区| 欧美国产精品va在线观看不卡| 亚洲国产欧美日韩在线播放| 老汉色∧v一级毛片| 美女午夜性视频免费| 丁香欧美五月| 99在线人妻在线中文字幕 | 男女床上黄色一级片免费看| 天堂动漫精品| 日韩大码丰满熟妇| 男女下面插进去视频免费观看| 亚洲国产av新网站| 丁香六月欧美| 亚洲专区字幕在线| 男女无遮挡免费网站观看| 大陆偷拍与自拍| 99热国产这里只有精品6| 亚洲精品av麻豆狂野| 黄色成人免费大全| 国产精品久久久人人做人人爽| 这个男人来自地球电影免费观看| 欧美成人免费av一区二区三区 | 一个人免费在线观看的高清视频| 午夜福利在线观看吧| 99久久精品国产亚洲精品| 黑人猛操日本美女一级片| 一本久久精品| 精品亚洲成国产av| 久久人妻福利社区极品人妻图片| 免费不卡黄色视频| 午夜91福利影院| 无遮挡黄片免费观看| 最近最新中文字幕大全免费视频| 亚洲avbb在线观看| 亚洲色图av天堂| 国产在线一区二区三区精| 18禁国产床啪视频网站| 免费女性裸体啪啪无遮挡网站| netflix在线观看网站| 国产亚洲精品第一综合不卡| 一夜夜www| 国产精品免费大片| 在线观看66精品国产| 亚洲国产毛片av蜜桃av| 亚洲精品中文字幕一二三四区 | 1024香蕉在线观看| 91麻豆av在线| 国产欧美日韩一区二区三区在线| 欧美激情高清一区二区三区| 亚洲情色 制服丝袜| 国产成人系列免费观看| 黄色片一级片一级黄色片| 不卡av一区二区三区| 成年动漫av网址| 中亚洲国语对白在线视频| 精品久久久久久电影网| 久久国产精品人妻蜜桃| 黑人操中国人逼视频| 一区二区三区激情视频| 波多野结衣av一区二区av| 久久人妻福利社区极品人妻图片| av又黄又爽大尺度在线免费看| 日日摸夜夜添夜夜添小说| 欧美成人免费av一区二区三区 | 男女床上黄色一级片免费看| 欧美乱妇无乱码| 香蕉久久夜色| 亚洲国产看品久久| 欧美亚洲 丝袜 人妻 在线| 午夜日韩欧美国产| av在线播放免费不卡| 免费观看人在逋| 熟女少妇亚洲综合色aaa.| 亚洲精品在线美女| 性高湖久久久久久久久免费观看| 久久久国产一区二区| 久久婷婷成人综合色麻豆| 亚洲av电影在线进入| 亚洲欧洲精品一区二区精品久久久| 精品国产超薄肉色丝袜足j| 中文字幕人妻丝袜制服| 久久影院123| 成人特级黄色片久久久久久久 | 老司机福利观看| 精品国产超薄肉色丝袜足j| 久久中文字幕人妻熟女| www.精华液| 国产亚洲一区二区精品| 日韩熟女老妇一区二区性免费视频| 亚洲人成伊人成综合网2020| 亚洲色图av天堂| videosex国产| 精品午夜福利视频在线观看一区 | 欧美日韩中文字幕国产精品一区二区三区 | 欧美乱码精品一区二区三区| 18在线观看网站| 一区二区av电影网| 欧美精品av麻豆av| 亚洲精品成人av观看孕妇| 亚洲欧美色中文字幕在线| 丁香欧美五月| 天堂动漫精品| 成人手机av| 国产成人精品久久二区二区免费| 搡老岳熟女国产| 肉色欧美久久久久久久蜜桃| 99riav亚洲国产免费| 国产aⅴ精品一区二区三区波| 999精品在线视频| 国产免费福利视频在线观看| 国产aⅴ精品一区二区三区波| 三上悠亚av全集在线观看| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 亚洲色图av天堂| 久久天堂一区二区三区四区| 国产av又大| 91九色精品人成在线观看| 99re在线观看精品视频| 大香蕉久久成人网| 成人精品一区二区免费| 国产精品麻豆人妻色哟哟久久| cao死你这个sao货| 欧美久久黑人一区二区| 久久久久久人人人人人| 叶爱在线成人免费视频播放| 免费av中文字幕在线| 精品国产超薄肉色丝袜足j| 国产男女超爽视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 大型av网站在线播放| 中文字幕高清在线视频| 国产亚洲欧美在线一区二区| 国产99久久九九免费精品| 一进一出好大好爽视频| 人妻久久中文字幕网| 免费看a级黄色片| 欧美日韩福利视频一区二区| 午夜精品国产一区二区电影| 欧美人与性动交α欧美精品济南到| 侵犯人妻中文字幕一二三四区| 免费女性裸体啪啪无遮挡网站| 美女福利国产在线| 亚洲久久久国产精品| 高清在线国产一区| 在线观看一区二区三区激情| 国产97色在线日韩免费| 午夜免费成人在线视频| 热99re8久久精品国产| 成人av一区二区三区在线看| 日韩有码中文字幕| 亚洲国产欧美日韩在线播放| 波多野结衣一区麻豆| 国产亚洲av高清不卡| 在线看a的网站| 欧美人与性动交α欧美软件| 久久国产精品大桥未久av| 亚洲精品av麻豆狂野| 久久久久精品人妻al黑| 国产精品一区二区在线不卡| 狠狠精品人妻久久久久久综合| 国产精品九九99| 久久九九热精品免费| 91麻豆av在线| 电影成人av| 亚洲avbb在线观看| 国产老妇伦熟女老妇高清| 久热这里只有精品99| 交换朋友夫妻互换小说| 国产亚洲欧美精品永久| 在线观看www视频免费| 精品国产乱码久久久久久男人| 黄色怎么调成土黄色| 国产淫语在线视频| 久久久国产成人免费| 久久久久网色| 亚洲九九香蕉| 国产精品欧美亚洲77777| 国产精品影院久久| 久久九九热精品免费| 欧美性长视频在线观看| 人人妻人人爽人人添夜夜欢视频| 人成视频在线观看免费观看| 久久久国产一区二区| 欧美人与性动交α欧美软件| av在线播放免费不卡| 如日韩欧美国产精品一区二区三区| 国产在视频线精品| 久久av网站| 国产97色在线日韩免费| 日韩中文字幕欧美一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 国产有黄有色有爽视频| 午夜成年电影在线免费观看| 老熟妇仑乱视频hdxx| 亚洲欧美一区二区三区黑人| 两个人看的免费小视频| videos熟女内射| 一区二区三区精品91| 王馨瑶露胸无遮挡在线观看| 热99久久久久精品小说推荐| 精品久久久久久电影网| 国产视频一区二区在线看| 国产精品自产拍在线观看55亚洲 | 国产免费视频播放在线视频| 另类亚洲欧美激情| 日本vs欧美在线观看视频| 精品一区二区三区视频在线观看免费 | 午夜久久久在线观看| 亚洲第一av免费看| www.精华液| 他把我摸到了高潮在线观看 | 国产三级黄色录像| 色综合婷婷激情| 人人妻人人澡人人爽人人夜夜| 美女福利国产在线| 俄罗斯特黄特色一大片| 操出白浆在线播放| 多毛熟女@视频| 一区二区日韩欧美中文字幕| 黑人猛操日本美女一级片| 一区福利在线观看| 久久久久精品国产欧美久久久| avwww免费| 法律面前人人平等表现在哪些方面| 一区福利在线观看| 搡老乐熟女国产| 国产欧美日韩一区二区三| 搡老乐熟女国产| 天堂动漫精品| 国产精品电影一区二区三区 | 三级毛片av免费| 蜜桃国产av成人99| 久久午夜亚洲精品久久| 欧美亚洲日本最大视频资源| 一进一出抽搐动态| 欧美精品高潮呻吟av久久| 欧美成人免费av一区二区三区 | 乱人伦中国视频| 黄色怎么调成土黄色| 五月天丁香电影| 国产不卡av网站在线观看| 欧美精品一区二区大全| 成人三级做爰电影| 国精品久久久久久国模美| 午夜精品国产一区二区电影| 日韩欧美免费精品| 精品国内亚洲2022精品成人 | 久久久精品国产亚洲av高清涩受| 丰满少妇做爰视频| 中文字幕人妻丝袜一区二区| 国产精品亚洲av一区麻豆| 一二三四在线观看免费中文在| 精品一品国产午夜福利视频| 国产欧美日韩综合在线一区二区| 午夜福利乱码中文字幕| 黄片大片在线免费观看| 日日爽夜夜爽网站| 丰满少妇做爰视频| 午夜福利免费观看在线| 国产精品一区二区免费欧美| 日韩欧美一区二区三区在线观看 | 免费观看av网站的网址| 国产99久久九九免费精品| 亚洲男人天堂网一区| 亚洲国产欧美在线一区| 丝瓜视频免费看黄片| 制服诱惑二区| 一级黄色大片毛片| 一级,二级,三级黄色视频| 99国产精品一区二区三区| 视频区欧美日本亚洲| 国产精品欧美亚洲77777| 精品一区二区三卡| 女人精品久久久久毛片| 伊人久久大香线蕉亚洲五| 青草久久国产| av免费在线观看网站| 久久中文字幕人妻熟女| 在线观看免费高清a一片| 十八禁网站免费在线| 黄色丝袜av网址大全| 欧美黑人精品巨大| 欧美成人免费av一区二区三区 | 久久精品亚洲精品国产色婷小说| 久久久国产精品麻豆| 亚洲精品av麻豆狂野| 成年人免费黄色播放视频| 久久影院123| 精品国产一区二区三区四区第35| av片东京热男人的天堂| 极品少妇高潮喷水抽搐| 老司机午夜十八禁免费视频| 少妇粗大呻吟视频| 十八禁人妻一区二区| 国产精品国产高清国产av | 欧美日韩亚洲综合一区二区三区_| 精品国内亚洲2022精品成人 | 国产精品久久久av美女十八| 国产成人免费观看mmmm| 欧美亚洲日本最大视频资源| 电影成人av| 美女国产高潮福利片在线看| 天天添夜夜摸| 亚洲精品久久午夜乱码| 国产伦理片在线播放av一区| 亚洲成国产人片在线观看| 久久人人97超碰香蕉20202| 久久九九热精品免费| 黄色成人免费大全| 国产一区二区三区在线臀色熟女 | 亚洲一码二码三码区别大吗| 成人精品一区二区免费| 国产欧美日韩一区二区三| 国产激情久久老熟女| 天天影视国产精品| 99热国产这里只有精品6| 国产精品98久久久久久宅男小说| 亚洲av第一区精品v没综合| 天天躁日日躁夜夜躁夜夜| 国产欧美亚洲国产| 免费久久久久久久精品成人欧美视频| 久久久国产欧美日韩av| bbb黄色大片| 久久久久久久大尺度免费视频| 亚洲五月婷婷丁香| 午夜两性在线视频| 久久久久国内视频| 一级黄色大片毛片| 99riav亚洲国产免费| 亚洲欧美精品综合一区二区三区| 黄色视频不卡| 久久青草综合色| 日本黄色日本黄色录像| 国产亚洲欧美精品永久| 天天影视国产精品| 可以免费在线观看a视频的电影网站| 18禁观看日本| 国产欧美日韩一区二区精品| 亚洲中文字幕日韩| 亚洲人成电影观看| 国产精品1区2区在线观看. | 国产日韩欧美亚洲二区| 欧美日韩福利视频一区二区| 黄网站色视频无遮挡免费观看| 电影成人av| 黑人巨大精品欧美一区二区mp4| 欧美精品一区二区大全| 免费日韩欧美在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品秋霞免费鲁丝片| 丝袜在线中文字幕| 波多野结衣一区麻豆| 欧美亚洲日本最大视频资源| 精品国产乱码久久久久久小说| 婷婷成人精品国产| 日本av免费视频播放| 超碰97精品在线观看| 麻豆国产av国片精品| 超碰97精品在线观看| 久久久久国内视频| 日韩中文字幕欧美一区二区| 一级片免费观看大全| 天天躁日日躁夜夜躁夜夜| 国产av精品麻豆| 美女福利国产在线| 久久精品人人爽人人爽视色| 中文字幕av电影在线播放| 日韩三级视频一区二区三区| 精品少妇久久久久久888优播| 少妇被粗大的猛进出69影院| 熟女少妇亚洲综合色aaa.| 国产日韩欧美在线精品| 欧美一级毛片孕妇| 男女午夜视频在线观看| 女警被强在线播放| 久久人人97超碰香蕉20202| 国产熟女午夜一区二区三区| 亚洲欧美日韩另类电影网站| 亚洲九九香蕉| 欧美 亚洲 国产 日韩一| 岛国毛片在线播放| 国产精品一区二区在线不卡| 一进一出好大好爽视频| 欧美成人免费av一区二区三区 | 久久午夜综合久久蜜桃| 日韩视频在线欧美| 欧美成人免费av一区二区三区 | 一二三四在线观看免费中文在| 视频在线观看一区二区三区| 免费日韩欧美在线观看| 国产福利在线免费观看视频| 成年人午夜在线观看视频| 久久人妻av系列| 国产淫语在线视频| 91麻豆av在线| 日韩成人在线观看一区二区三区| 亚洲国产av新网站| 丝瓜视频免费看黄片| 一区二区av电影网| 涩涩av久久男人的天堂| 最新在线观看一区二区三区| 人人妻,人人澡人人爽秒播| 精品一区二区三区四区五区乱码| 久久久久网色| 亚洲av日韩在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av美国av| 国产精品欧美亚洲77777| 久久久国产欧美日韩av| 精品久久久久久电影网| 男人舔女人的私密视频| tocl精华| 老熟妇仑乱视频hdxx| 一区在线观看完整版| 亚洲精品在线美女| 国产在线视频一区二区| 亚洲天堂av无毛| 成人特级黄色片久久久久久久 | a在线观看视频网站| 在线观看66精品国产| 亚洲一区二区三区欧美精品| 又紧又爽又黄一区二区| 亚洲人成电影观看| 国产精品99久久99久久久不卡| 最近最新中文字幕大全免费视频| av电影中文网址| 成在线人永久免费视频| 亚洲自偷自拍图片 自拍| 大陆偷拍与自拍| 亚洲av成人一区二区三| 90打野战视频偷拍视频| 91精品国产国语对白视频| 老汉色av国产亚洲站长工具|