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

    CATIA-FLAC3D耦合建模方法及其應(yīng)用

    2014-03-05 07:25:06喬世范謝濟仁郭麒麟許文龍
    關(guān)鍵詞:網(wǎng)格程序建模

    喬世范,謝濟仁,郭麒麟,許文龍

    (1.中南大學(xué) 土木工程學(xué)院,長沙410075;2.長江委巖石工程總公司(武漢),武漢430010)

    FLAC3D(Fast Lagrangian Analysis of Continua)是以拉格朗日法為基礎(chǔ)編寫的有限差分程序,三維快速拉格朗日法是基于三維顯式的一種有限差分法的數(shù)值分析方法,它可以模擬巖土或其他材料的三維力學(xué)行為。自20世紀(jì)90年代以來,已經(jīng)在土建、交通、采礦、地質(zhì)、水利等工業(yè)部門廣泛應(yīng)用,逐漸成為巖土工程界的重要工具之一。

    然而FLAC3D軟件在模型前處理上存在很大缺陷,主要表現(xiàn)為:建模過程不直觀、檢查不方便、工作量大耗費時間長。為了解決這些問題,一些學(xué)者做了一些工作:王樹仁等[1]基于 MIDAS/GTS軟件,采用 MATLAB語言編寫 MIDAS/GTS-FLAC3D接口程序。廖秋林等[2]基于 ANSYS軟件,采用ANSYS語言編寫了FLAC3D-ANSYS接口程序,借助有限元軟件ANSYS相對便捷的前處理功能,實現(xiàn)了層狀地質(zhì)體FLAC3D模型的建立。胡斌等[3]采用FORTRAN語言編寫了FLAC3D的前處理程序,對于簡單地質(zhì)體實現(xiàn)了快速、便捷的建模;Aringoli等利用Surfer和GID處理,建立地質(zhì)體模型,通過Fish編 程 將 建 好 的 網(wǎng) 格 讀 入 到 FLAC3D中[4-7]。ANSYS、MIDAS/GTS等有限元計算軟件,針對規(guī)則的地質(zhì)體能夠快速建立地質(zhì)體的幾何模型和有限元模型,然而對于復(fù)雜的地質(zhì)體,其建模過程復(fù)雜、工作量大、花費時間長[8-11]。法國 Dassault公司開發(fā)的CATIA軟件是集CAD/CAM/CAE于一體的優(yōu)秀三維設(shè)計軟件,在機械、電子、航空、航天和汽車等行業(yè)獲得很廣泛應(yīng)用。它從V5R19版本開始,增加了地形地質(zhì)建模功能,用于土木工程的規(guī)劃與設(shè)計以及基礎(chǔ)工程等。一些單位在相關(guān)課題中已應(yīng)用CATIA進行三維地質(zhì)建模,開發(fā)了相應(yīng)的地質(zhì)插件,實現(xiàn)了地質(zhì)點的轉(zhuǎn)換、鉆孔數(shù)據(jù)的導(dǎo)入、鉆孔信息的分層、平面圖的導(dǎo)入以及剖面圖的導(dǎo)入等操作,從而簡化了建模過程,具有比較好的實用性。然而,CATIA雖然有比較強的建模能力,但其有限元分析模塊比較薄弱,所以需要開發(fā)與巖土工程專業(yè)軟件FLAC3D的接口程序。

    針對FLAC3D前處理建模存在的技術(shù)不足[12],筆者借助CATIA進行復(fù)雜地質(zhì)體及工程結(jié)構(gòu)的幾何建模和網(wǎng)格劃分,通過編寫CATIA-FLAC3D數(shù)據(jù)轉(zhuǎn)換接口程序,將建好的模型導(dǎo)入到FLAC3D中,從而降低FLAC3D前處理建模的難度,實現(xiàn)復(fù)雜工程的FLAC3D三維模型快速、準(zhǔn)確構(gòu)建,彌補了CATIA在有限元計算方面的不足。

    1 復(fù)雜地質(zhì)體的建立和網(wǎng)格劃分

    1.1 CATIA軟件簡介

    CATIA軟件包括多個功能模塊,用于三維地質(zhì)建模的模塊主要包括有數(shù)字曲面編輯(DSE)、創(chuàng)成式外形設(shè)計(GSD)、快速曲面重建(QSR)和雕刻(Shape Sculptor)模塊等。采用CATIA建立地質(zhì)模型,建模方便,模型準(zhǔn)確。

    1.2 CATIA中建立地質(zhì)模型

    1.2.1 點云和地表Mesh面的形成 根據(jù)測繪部門所提供的工程地形圖,通過轉(zhuǎn)化軟件,生成地形圖等高線所有點的.asc文件。通過CATIA中DSE模塊導(dǎo)入,形成點云圖,見圖1。

    將點云進行過濾、刪除明顯錯誤的數(shù)據(jù)點,通過mesh creation命令形成3D三角小網(wǎng)格面片組,采用DSE平臺中的網(wǎng)格面修改功能對三角網(wǎng)格面片組進行補點、刪點、修改網(wǎng)格形狀等操作,生成合理網(wǎng)格面,見圖2。

    圖1 導(dǎo)入的點云

    1.2.2 形成地表曲面 進入QSR模塊,點擊Power Fit命令,將網(wǎng)格面強制實體化,調(diào)整界面參數(shù),生成接近真實情況的實體面。如圖3。

    1.2.3 生成地質(zhì)體模型 在地表面范圍內(nèi)建立凸臺,進入零件設(shè)計模塊,點擊分割命令(插入/基于曲面的特征/分割),其中,分割元素為上步生成的地表面,生成的地質(zhì)實體模型如圖4所示。

    圖4 地質(zhì)體模型

    1.2.4 生成網(wǎng)格體模型 進入分析與模擬模塊的高級網(wǎng)格劃分工具(Advanced Meshing Tools)模塊,點擊 Meshing Methods工具欄下的網(wǎng)格劃分工具(OCTREE Tetrahedron Mesh)命令,設(shè)定單元體大小,最后生成地質(zhì)體有限元網(wǎng)格模型[13-14],如圖5所示。

    圖5 網(wǎng)格體模型

    2 CATIA-FLAC3D接口程序

    CATIA有限元模型可采用采用四面體、六面體進行劃分,其單元形狀與FLAC3D所采用的四面體形網(wǎng)格和矩形網(wǎng)格形狀相似,但其每一單元節(jié)點編制的規(guī)則和節(jié)點坐標(biāo)(單元數(shù)據(jù))有差異。通過深入分析兩個軟件數(shù)據(jù)文件的內(nèi)在聯(lián)系,在VS2010的平臺下,使用Visual Basic語言編寫了CATIA-FLAC3D接口程序,實現(xiàn)了在CATIA中建模和網(wǎng)格劃分,在FLAC3D中計算,使兩個軟件的優(yōu)勢得到發(fā)揮。

    2.1 FLAC3D與CATIA單元數(shù)據(jù)關(guān)系

    CATIA主要提供四面體(Ⅰ)、四面體(Ⅱ)和六面體3種實體單元,可分別對應(yīng)于FLAC3D中的Tetrahedron和Brick單元[15],其單元節(jié)點編號對應(yīng)關(guān)系見表1(以四面體(Ⅰ)為例)。

    表1 CATIA與FLAC3D單元數(shù)據(jù)關(guān)系對照

    CATIA導(dǎo)出的節(jié)點、單元數(shù)據(jù)格式與FLAC3D可識別的節(jié)點、單元數(shù)據(jù)格式對照分別見表2(以四面體(I)為例)和表3(以四面體(I)為例)所列內(nèi)容。

    表2 CATIA與FLAC3D節(jié)點數(shù)據(jù)格式對照表

    表3 CATIA與FLAC3D單元數(shù)據(jù)格式對照表

    2.2 CATIA與FLAC3D數(shù)據(jù)轉(zhuǎn)換思路

    基于CATIA與FLAC3D單元數(shù)據(jù)關(guān)系,應(yīng)用Visual Basic編寫了CATIA-FLAC3D接口程序。將CATIA導(dǎo)出的節(jié)點、單元數(shù)據(jù)通過編程存儲到幾個動態(tài)數(shù)組中,再利用CATIA-FLAC3D數(shù)據(jù)轉(zhuǎn)換接口程序轉(zhuǎn)換成FLAC3D可識別的數(shù)據(jù)文件,并存入一個.txt文件中,通過調(diào)用這個文件,在FLAC3D中生成有限元模型;依次施加邊界條件和初始條件、對材料參數(shù)賦值等,進行計算[13-15]。

    2.3 CATIA與FLAC3D數(shù)據(jù)轉(zhuǎn)換實現(xiàn)過程

    CATIA與FLAC3D數(shù)據(jù)轉(zhuǎn)換接口程序?qū)崿F(xiàn)過程如下:

    2.3.1 節(jié)點數(shù)據(jù)轉(zhuǎn)換 將CATIA的節(jié)點數(shù)據(jù)轉(zhuǎn)換為FLAC3D數(shù)據(jù)文件,需刪掉非數(shù)據(jù)列(即“GRID*”列)以及多余的節(jié)點號列,當(dāng)遇到關(guān)鍵字“GR”,開始存儲節(jié)點數(shù)據(jù),將節(jié)點號和節(jié)點坐標(biāo)分別存儲到動態(tài)數(shù)組node number()和axis()數(shù)組中。部分關(guān)鍵代碼如下:

    2.3.2 單元數(shù)據(jù)轉(zhuǎn)換 將CATIA的單元數(shù)據(jù)轉(zhuǎn)換為FLAC3D數(shù)據(jù)文件。主要流程如下:

    1)讀取單元編號,將單元編號存儲到動態(tài)數(shù)組element number()。

    2)讀取單元的所有節(jié)點,將這些節(jié)點存入動態(tài)數(shù)組element contact()。

    3)將網(wǎng)格組存入動態(tài)數(shù)組element group()。

    部分關(guān)鍵代碼如下:

    2.3.3 單元分組 FLAC3D可識別的單元分組信息格式為GROUP組名。由表3可見,CATIA導(dǎo)出的單元數(shù)據(jù)中包括單元所屬的網(wǎng)格組信息。每個單元都會與網(wǎng)格組類型對應(yīng),此信息存儲于數(shù)組element group()內(nèi),由此,可以通過網(wǎng)格組信息對單元進行分組,這樣就可實現(xiàn)單元分組信息的識別。筆者利用Visual Basic編寫的轉(zhuǎn)換程序,對每個單元進行分組,知道所有單元分組完成,程序終止。部分關(guān)鍵代碼如下:

    2.3.4 轉(zhuǎn)換接口的流程圖及界面 執(zhí)行完以上操作后,將轉(zhuǎn)換完畢的節(jié)點、單元數(shù)據(jù)及單元分組信息輸出到一個.txt文 件中,在 FLAC3D調(diào)用[16-18]。流程見圖6。CATIA-FLAC3D的接口程序界面如圖7所示。

    圖6 算法流程圖

    圖7 CATIA-FLAC3D接口程序界面

    3 應(yīng)用實例

    3.1 工程背景

    以某已建成的偏壓隧道為依托,驗證該方法的可行性和有效性。該地區(qū)的特點是地表形態(tài)比較復(fù)雜,巖層主要以石灰?guī)r為主,其力學(xué)參數(shù)見表4。

    表4 巖石參數(shù)

    使用ANASYS,MIDAS等軟件建模,工程量大,建模繁瑣,而直接使用FLAC3D建模,建模周期長,費時費力[16,19]。采用該方法,能夠方便快速的建立幾何模型和有限元模型,利用CATIA-FLAC3D接口程序很方便的將模型導(dǎo)入FLAC3D中進行分析。

    3.2 有限元計算模型及計算結(jié)果

    根據(jù)以上流程,在CATIA中生成的有限元模型如圖8。導(dǎo)入到FLAC3D中的有限元模型如圖9。該模型上表面為自由面,不需要施加約束,其他5個面都施加法向約束,施加重力場,材料參數(shù)如表4所示,本構(gòu)模型采用摩爾庫倫,見圖10。

    圖8 CATIA有限元計算模型

    圖9 FLAC3D有限元計算模型

    圖10 Mohr-Coulomb屈服面

    式中:c為凝聚力;φ為內(nèi)摩擦角;s和θf分別為破壞面上的剪切應(yīng)力和法向應(yīng)力。

    如果用主應(yīng)力表示,式(1)可以改寫為:

    式中:θ1和θ3分別為第一主應(yīng)力和第三主應(yīng)力;c表示凝聚力;φ為內(nèi)摩擦角。計算結(jié)果如圖11[12,15]所示。

    圖11 Z方向應(yīng)力云圖

    工程實例驗證表明,CATIA-FLAC3D耦合建模方法是有效、可行的。

    4 結(jié) 論

    1)基于CATIA平臺進行三維建模,生成有限元模型,通過 Visual Basic語言編寫了 CATIAFLAC3D的接口程序,實現(xiàn)了在CATIA中建模,在FLAC3D中進行有限元計算的功能。

    2)CATIA作為一款優(yōu)秀的設(shè)計建模軟件,但它的有限元分析計算模塊并不強大;FLAC3D是巖土方向的專業(yè)有限元計算軟件,它在前期處理上存在著一定的缺陷;通過CATIA-FLAC3D接口程序,能夠利用CATIA方便、快速、準(zhǔn)確地建立地質(zhì)模型,并利用FLAC3D進行有限元計算,減少了建模所需要的時間和精力,提高了效率,縮短了計算周期。

    [1]王樹仁,張海清.MIDAS/GTS-FLAC3D耦合建模新方法及其應(yīng)用[J].土木建筑與環(huán)境工程,2010,32(1):12-17.Wang S R,Zhang H Q.A coupling modeling method with MIDAS/GTS-FLAC3Dand its application [J].Journal of Civil, Architectural & Environmental Engineering,2010,32(1):12-17.

    [2]廖秋林,曾錢幫.基于ANSYS平臺復(fù)雜地質(zhì)體FLAC3D模型的自動生成[J].巖石力學(xué)與工程學(xué)報,2005,24(6):1010-1013.Liao Q L,Zeng Q B.Automatic model Generation of complex geologic body with FLAC3Dbased on ANSYS platform [J].Chinese Journal of Rock Mechanics and Engineering,2005,24(6):1010-1013.

    [3]胡斌,張倬元,黃潤秋,等.FLAC3D前處理程序的開發(fā)及仿真效果檢驗[J].巖石力學(xué)與工程學(xué)報,2002,21(9):1387-1391.Hu B,Zhang Z Y,Huang R Q,et al.Development of pre-processing package for FLAC3Dand verification of its simulating effects [J].Chinese Journal of Rock Mechanics and Engineering,2002,21(9):1387-1391.

    [4]羅周全,吳亞斌,劉曉明,等.基于SURPAC的復(fù)雜地質(zhì)體FLAC3D模型生成技術(shù)[J].巖土力學(xué),2008,29(5):1334-1338.Luo Z Q,Wu Y B,Liu X M,et al.FLAC3Dmodeling for complex geologic body based on SURPAC [J].Rock and Soil Mechanics,2008,29(5):1334-1338.

    [5]Jenck O,Dias D,Kastner R.Three-dimensional numerical modeling of a piled Embankment [J].International Journal of Geomechanics,2009,9(3):102-112.

    [6]Chugh A K,Stark T D,Dejong K A.Reanalysis of a municipal landfill slope failure near Cincinnati,Ohio,USA [J].Canadian Geotechnical Journal,2007,44(1):33-53.

    [7]Aringoli D, Calista M, Gentili B, et al.Geomorphological features and 3Dmodelling of Montelparo mass movement (Central Italy) [J].Engineering Geology,2008,99(1/2):70-84.

    [8]侯恩科,吳立新,李建民,等.三維地學(xué)模擬與數(shù)值模擬的耦合方法研究[J].煤炭學(xué)報,2002,27(4):388-392.Hou E K,Wu L X,Li J M,et al.Study on the coupling of 3Dgeoscience modeling with numerical simulation[J].Journal of China Coal Society,2002,27(4):388-392.

    [9]王明華,白云.層狀巖體三維可視化構(gòu)模與數(shù)值模擬的集成研究[J].巖土力學(xué),2005,26(7):1123-1126.Wang M H,Bai Y.Study on integration of threedimensional modeling and numerical simulation for stratified rock mass [J].Rock and Soil Mechanics,2005,26(7):1123-1126.

    [10]李明超,鐘登華,秦朝霞,等.基于三維地質(zhì)模型的工程巖體結(jié)構(gòu)精細數(shù)值建模[J].巖石力學(xué)與工程學(xué)報,2007,26(9):1893-1898.Li M C,Zhong D H,Qin Z X,et al.Refined modeling for numerical simulation of engineering rock mass structures based on 3Dgeological model[J].Chinese Journal of Rock Mechanics and Engineering,2007,26(9):1893-1898.

    [11]Papakonstantinou S,Schubert W.Para-metric study of displacements on sequential and full face tunnel excavation [C]//Proceedings of the 33rd ITA-AITES World Tunnel Congress-Underground Space-The 4th Dimension of Metropolises,2007:639-643.

    [12]Itasca Consulting Group,Inc.FLAC3D,F(xiàn)ast Lagrangian Analysis of Continua in 3Dimensions,version 2.0,User’s Manual[R].USA Itasca Consulting Group,Inc.,1997.

    [13]劉曉明,羅周全,楊彪,等.復(fù)雜礦區(qū)三維地質(zhì)可視化及數(shù)值模型構(gòu)建[J].巖土力學(xué),2010,31(12):4006-4015.Liu X M,Luo Z Q,Yang B,et al.Numerical modeling and geological body visualization for complex mine[J].Rock and Soil Mechanics,2010,31(12):4006-4015.

    [14]劉秀軍.基于GOCAD的復(fù)雜地質(zhì)體FLAC3D模型生成技術(shù)[J].中國地質(zhì)災(zāi)害與防治學(xué)報,2011,22(4):41-45.Liu X J.FLAC3Dmodeling for complex geologic body based on GOCAD [J].The Chinese Journal of Geological Hazard and Control,2011,22(4):41-45.

    [15]伍永平,高永剛,解盤石.基于AutoCAD的FLAC3D地下工程快速建模方法研究[J].煤炭工程,2011(12):61-64.Wu Y P,Gao Y G,Xie P S.Study on rapid modeling method of FLAC3Dunderground project based on AutoCAD [J].Coal Engineering,2011(12):61-64.

    [16]王陽平,崔穎輝.基于 Midas/Gts的FLAC3D的建模方法[J].北方工業(yè)大學(xué)學(xué)報,2010,22(3):78-81.Wang Y P,Cui Y H.Method of modeling the FLAC3Dbased on midas/Gts [J].Journal of North China University of Technology,2010,22(3):78-81.

    [17]Lin H,Liu T Y,Li J T.A simple generation technique of complex geotechnical computational model[J].The Scientific World Journal,2013,37(8):1-8.

    [18]Zhang X D,Zhang B,Jiang L Z,et al.Research on modeling technology for three-dimension analysis of irregular slope stability in surface mines and its application[J].Advanced Materials Research,2013,796:838-841.

    [19]盧志剛,劉興權(quán),唐義宏.基于GIS和FLAC3D礦山地表沉陷可視化[J].現(xiàn)代礦業(yè),2011(11):44-46.

    猜你喜歡
    網(wǎng)格程序建模
    用全等三角形破解網(wǎng)格題
    聯(lián)想等效,拓展建?!浴皫щ娦∏蛟诘刃鲋凶鰣A周運動”為例
    反射的橢圓隨機偏微分方程的網(wǎng)格逼近
    試論我國未決羈押程序的立法完善
    基于PSS/E的風(fēng)電場建模與動態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    不對稱半橋變換器的建模與仿真
    “程序猿”的生活什么樣
    重疊網(wǎng)格裝配中的一種改進ADT搜索方法
    英國與歐盟正式啟動“離婚”程序程序
    基于曲面展開的自由曲面網(wǎng)格劃分
    久热这里只有精品99| 美女扒开内裤让男人捅视频| 亚洲 欧美一区二区三区| 美女国产高潮福利片在线看| 午夜影院在线不卡| 亚洲一码二码三码区别大吗| 亚洲精品乱久久久久久| 日韩欧美精品免费久久| 91成人精品电影| 成人免费观看视频高清| 精品一区二区免费观看| 精品少妇内射三级| 精品亚洲乱码少妇综合久久| 精品国产一区二区三区久久久樱花| 制服诱惑二区| 国精品久久久久久国模美| 秋霞在线观看毛片| 亚洲精华国产精华液的使用体验| 久久久亚洲精品成人影院| 777米奇影视久久| 嫩草影院入口| 狠狠精品人妻久久久久久综合| 亚洲欧美成人精品一区二区| 青草久久国产| 免费观看性生交大片5| 国产精品一二三区在线看| 午夜福利一区二区在线看| 亚洲国产精品一区二区三区在线| 夫妻午夜视频| 久久人人爽人人片av| 日韩一区二区视频免费看| 欧美xxⅹ黑人| www.精华液| 日本欧美国产在线视频| 午夜福利乱码中文字幕| 亚洲av中文av极速乱| 美女扒开内裤让男人捅视频| 欧美成人午夜精品| 国产一区二区三区综合在线观看| 亚洲欧美成人综合另类久久久| 我的亚洲天堂| 亚洲av综合色区一区| 日韩视频在线欧美| 国产日韩欧美在线精品| 久久人人爽人人片av| 精品亚洲成a人片在线观看| www.av在线官网国产| 亚洲美女黄色视频免费看| 黑人欧美特级aaaaaa片| 久久人人爽人人片av| 久久天堂一区二区三区四区| 亚洲精品视频女| √禁漫天堂资源中文www| 美女国产高潮福利片在线看| a级片在线免费高清观看视频| 日韩一本色道免费dvd| √禁漫天堂资源中文www| 97人妻天天添夜夜摸| 日韩精品免费视频一区二区三区| 亚洲第一区二区三区不卡| 最近最新中文字幕免费大全7| 99九九在线精品视频| av线在线观看网站| 欧美日韩亚洲综合一区二区三区_| 亚洲av成人精品一二三区| 国产1区2区3区精品| 亚洲精品国产av蜜桃| 国产一区二区三区综合在线观看| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| 成年动漫av网址| 国产成人精品久久二区二区91 | 一区二区三区四区激情视频| 妹子高潮喷水视频| 国产一区二区激情短视频 | 母亲3免费完整高清在线观看| 2018国产大陆天天弄谢| 一级片'在线观看视频| 亚洲国产av新网站| 天天躁夜夜躁狠狠久久av| 这个男人来自地球电影免费观看 | 午夜激情av网站| 国产片特级美女逼逼视频| 亚洲欧美中文字幕日韩二区| 麻豆乱淫一区二区| 黄色 视频免费看| videos熟女内射| 午夜影院在线不卡| 尾随美女入室| 成年美女黄网站色视频大全免费| 成人手机av| 中文字幕最新亚洲高清| 欧美 日韩 精品 国产| av有码第一页| 国产亚洲欧美精品永久| 热re99久久精品国产66热6| 在线精品无人区一区二区三| 亚洲国产av新网站| 亚洲国产精品成人久久小说| 亚洲,欧美,日韩| 哪个播放器可以免费观看大片| 久久99热这里只频精品6学生| 飞空精品影院首页| 精品国产乱码久久久久久男人| 少妇人妻久久综合中文| 国产精品熟女久久久久浪| 久久精品aⅴ一区二区三区四区| 亚洲精品aⅴ在线观看| 中文字幕色久视频| 久久久久国产精品人妻一区二区| 熟女av电影| 日本午夜av视频| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 精品国产露脸久久av麻豆| 久久久久视频综合| 最近的中文字幕免费完整| 宅男免费午夜| 视频区图区小说| 中文字幕高清在线视频| 国产日韩欧美亚洲二区| 日韩av不卡免费在线播放| 国产高清国产精品国产三级| 午夜精品国产一区二区电影| 久久精品国产亚洲av高清一级| 熟女少妇亚洲综合色aaa.| 免费黄色在线免费观看| 成人免费观看视频高清| 男女边吃奶边做爰视频| 九九爱精品视频在线观看| 黑丝袜美女国产一区| xxx大片免费视频| 国产精品一区二区精品视频观看| 国产在线视频一区二区| 精品午夜福利在线看| 老司机在亚洲福利影院| 国产在线视频一区二区| 久久久久国产一级毛片高清牌| 国产 精品1| a级片在线免费高清观看视频| 中文字幕人妻熟女乱码| 免费日韩欧美在线观看| 亚洲精华国产精华液的使用体验| 少妇人妻久久综合中文| 国产精品久久久久久精品电影小说| 黑人巨大精品欧美一区二区蜜桃| 丝袜喷水一区| 久久毛片免费看一区二区三区| 51午夜福利影视在线观看| 一级,二级,三级黄色视频| 精品一区二区免费观看| 亚洲情色 制服丝袜| 久久 成人 亚洲| 五月天丁香电影| 欧美亚洲日本最大视频资源| 黑人巨大精品欧美一区二区蜜桃| 人人妻,人人澡人人爽秒播 | 伊人久久大香线蕉亚洲五| 久久天堂一区二区三区四区| 亚洲av欧美aⅴ国产| 久久这里只有精品19| 国产在视频线精品| 国产成人一区二区在线| 亚洲美女视频黄频| 最近的中文字幕免费完整| 久久鲁丝午夜福利片| 成人国产av品久久久| 欧美日韩一区二区视频在线观看视频在线| av在线观看视频网站免费| 99热国产这里只有精品6| 亚洲av男天堂| 大片免费播放器 马上看| 日韩视频在线欧美| 又大又黄又爽视频免费| 欧美精品亚洲一区二区| 捣出白浆h1v1| 多毛熟女@视频| 欧美成人精品欧美一级黄| 叶爱在线成人免费视频播放| 亚洲情色 制服丝袜| 高清在线视频一区二区三区| 1024视频免费在线观看| 天天躁日日躁夜夜躁夜夜| 秋霞伦理黄片| 亚洲美女视频黄频| 丰满饥渴人妻一区二区三| 人人妻人人澡人人爽人人夜夜| 色综合欧美亚洲国产小说| 精品国产超薄肉色丝袜足j| 久久精品人人爽人人爽视色| 成人18禁高潮啪啪吃奶动态图| 成年av动漫网址| 纵有疾风起免费观看全集完整版| 99久久综合免费| 精品一品国产午夜福利视频| 久久这里只有精品19| 国产黄色视频一区二区在线观看| 在线天堂最新版资源| 中文字幕最新亚洲高清| 国产免费又黄又爽又色| bbb黄色大片| 男人爽女人下面视频在线观看| 欧美黄色片欧美黄色片| 国产一区亚洲一区在线观看| 中文字幕制服av| 国产成人精品无人区| 观看美女的网站| 国产亚洲精品第一综合不卡| 男人操女人黄网站| 欧美日韩精品网址| 成年人免费黄色播放视频| 在线观看三级黄色| 日日啪夜夜爽| av免费观看日本| 五月天丁香电影| 日韩,欧美,国产一区二区三区| 色综合欧美亚洲国产小说| 欧美亚洲日本最大视频资源| 午夜日本视频在线| 免费在线观看视频国产中文字幕亚洲 | 热99久久久久精品小说推荐| 各种免费的搞黄视频| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看| 亚洲国产精品一区三区| 亚洲国产看品久久| 亚洲伊人久久精品综合| 精品国产超薄肉色丝袜足j| 天堂中文最新版在线下载| 我的亚洲天堂| 久久天堂一区二区三区四区| 国产亚洲欧美精品永久| 天天躁夜夜躁狠狠躁躁| 两性夫妻黄色片| 国产免费视频播放在线视频| 欧美亚洲日本最大视频资源| 国产精品二区激情视频| 成人手机av| 一本—道久久a久久精品蜜桃钙片| 男女国产视频网站| 丰满少妇做爰视频| 三上悠亚av全集在线观看| 另类亚洲欧美激情| 国产免费现黄频在线看| 免费在线观看黄色视频的| 精品第一国产精品| 美女国产高潮福利片在线看| 日韩精品免费视频一区二区三区| 午夜福利网站1000一区二区三区| 成年女人毛片免费观看观看9 | 欧美乱码精品一区二区三区| 国产在线免费精品| 人人澡人人妻人| 日韩视频在线欧美| 一本—道久久a久久精品蜜桃钙片| 下体分泌物呈黄色| 欧美人与善性xxx| 国产有黄有色有爽视频| 欧美 日韩 精品 国产| 午夜福利,免费看| 亚洲欧美一区二区三区黑人| 亚洲国产欧美在线一区| 亚洲五月色婷婷综合| 欧美激情 高清一区二区三区| 国产 一区精品| 精品少妇内射三级| 亚洲精品国产av蜜桃| 交换朋友夫妻互换小说| 欧美日韩亚洲国产一区二区在线观看 | 无限看片的www在线观看| 亚洲少妇的诱惑av| 看十八女毛片水多多多| 成人黄色视频免费在线看| 亚洲免费av在线视频| 亚洲 欧美一区二区三区| 亚洲图色成人| 美女中出高潮动态图| 伊人亚洲综合成人网| 午夜久久久在线观看| 中文乱码字字幕精品一区二区三区| 国产日韩欧美在线精品| 天天添夜夜摸| av网站在线播放免费| 亚洲精品日本国产第一区| av视频免费观看在线观看| 国产精品久久久久久人妻精品电影 | 国产xxxxx性猛交| 亚洲欧美色中文字幕在线| 天天躁夜夜躁狠狠久久av| 国产亚洲午夜精品一区二区久久| 免费黄色在线免费观看| 午夜免费男女啪啪视频观看| av福利片在线| 日韩欧美精品免费久久| 亚洲国产毛片av蜜桃av| 亚洲国产精品一区二区三区在线| 国产精品无大码| svipshipincom国产片| 成年人午夜在线观看视频| 国产淫语在线视频| 国产成人精品久久二区二区91 | 久久国产亚洲av麻豆专区| www.自偷自拍.com| 99精国产麻豆久久婷婷| 各种免费的搞黄视频| 超碰成人久久| 国产日韩欧美在线精品| 黄色视频不卡| 美女主播在线视频| 午夜日韩欧美国产| 十八禁高潮呻吟视频| 中文字幕人妻丝袜制服| 一级毛片 在线播放| 亚洲欧美精品综合一区二区三区| 亚洲国产精品999| 男女下面插进去视频免费观看| 高清av免费在线| 激情视频va一区二区三区| 另类亚洲欧美激情| 最黄视频免费看| 卡戴珊不雅视频在线播放| 黑人巨大精品欧美一区二区蜜桃| 叶爱在线成人免费视频播放| 国产成人91sexporn| 国语对白做爰xxxⅹ性视频网站| 国产又爽黄色视频| 天天添夜夜摸| 视频区图区小说| 午夜免费鲁丝| 久久人人97超碰香蕉20202| 18禁国产床啪视频网站| 成人手机av| 卡戴珊不雅视频在线播放| 亚洲国产精品国产精品| 精品国产一区二区三区久久久樱花| 久久久国产一区二区| 在线免费观看不下载黄p国产| 伦理电影大哥的女人| 黄片播放在线免费| 尾随美女入室| 丰满乱子伦码专区| 亚洲成人国产一区在线观看 | 亚洲一区中文字幕在线| 国产成人精品无人区| 超碰成人久久| 国产精品av久久久久免费| 欧美少妇被猛烈插入视频| av在线老鸭窝| 又粗又硬又长又爽又黄的视频| 一个人免费看片子| bbb黄色大片| 亚洲精品国产av蜜桃| 99精品久久久久人妻精品| 美女视频免费永久观看网站| 日本午夜av视频| 老汉色∧v一级毛片| 黄片无遮挡物在线观看| 国产一级毛片在线| 精品国产乱码久久久久久小说| 女人被躁到高潮嗷嗷叫费观| 视频区图区小说| 18禁观看日本| 国产一区二区在线观看av| 国产深夜福利视频在线观看| 七月丁香在线播放| 一个人免费看片子| 久久久国产一区二区| 日韩 欧美 亚洲 中文字幕| 精品久久蜜臀av无| 日韩精品免费视频一区二区三区| 夫妻午夜视频| 黄片小视频在线播放| 国产人伦9x9x在线观看| 国产伦理片在线播放av一区| 色视频在线一区二区三区| 久久精品国产亚洲av涩爱| 亚洲精品久久午夜乱码| 精品酒店卫生间| 午夜福利视频精品| 99国产综合亚洲精品| 午夜91福利影院| 叶爱在线成人免费视频播放| 免费日韩欧美在线观看| 日韩一本色道免费dvd| 免费黄色在线免费观看| 91aial.com中文字幕在线观看| 妹子高潮喷水视频| 久久久久精品久久久久真实原创| 黑人巨大精品欧美一区二区蜜桃| 成人国产av品久久久| 操美女的视频在线观看| 亚洲国产毛片av蜜桃av| 美女扒开内裤让男人捅视频| 久久毛片免费看一区二区三区| 成人影院久久| 捣出白浆h1v1| 精品卡一卡二卡四卡免费| 精品国产一区二区三区久久久樱花| 天天躁日日躁夜夜躁夜夜| 在线天堂最新版资源| 欧美人与性动交α欧美软件| 高清黄色对白视频在线免费看| 丝袜脚勾引网站| 日韩一本色道免费dvd| 久久人人爽人人片av| 亚洲精品,欧美精品| 999精品在线视频| 在线观看免费午夜福利视频| 新久久久久国产一级毛片| 少妇人妻久久综合中文| 天天躁狠狠躁夜夜躁狠狠躁| 美女扒开内裤让男人捅视频| 国产av一区二区精品久久| 多毛熟女@视频| 国产精品久久久人人做人人爽| 青春草亚洲视频在线观看| 亚洲,欧美,日韩| 啦啦啦中文免费视频观看日本| 日韩熟女老妇一区二区性免费视频| 亚洲国产欧美网| av有码第一页| 国产精品久久久久成人av| 成人毛片60女人毛片免费| 欧美日韩av久久| 汤姆久久久久久久影院中文字幕| 欧美 日韩 精品 国产| 精品国产一区二区久久| 亚洲欧美色中文字幕在线| 亚洲精品国产一区二区精华液| 少妇被粗大的猛进出69影院| xxxhd国产人妻xxx| 久久女婷五月综合色啪小说| 国产精品嫩草影院av在线观看| 国产精品 国内视频| 亚洲三区欧美一区| 精品国产乱码久久久久久小说| 亚洲伊人色综图| av又黄又爽大尺度在线免费看| 午夜福利乱码中文字幕| 欧美日韩视频高清一区二区三区二| 国产成人精品福利久久| 人妻人人澡人人爽人人| a级毛片在线看网站| 在线观看www视频免费| 欧美日韩视频精品一区| 国产乱人偷精品视频| 爱豆传媒免费全集在线观看| 亚洲国产欧美一区二区综合| 亚洲欧洲国产日韩| 777久久人妻少妇嫩草av网站| 国产男女超爽视频在线观看| av一本久久久久| 三上悠亚av全集在线观看| 久久天躁狠狠躁夜夜2o2o | 好男人视频免费观看在线| 亚洲久久久国产精品| 国产亚洲午夜精品一区二区久久| 国产精品麻豆人妻色哟哟久久| 国产麻豆69| 欧美 日韩 精品 国产| 亚洲精品,欧美精品| 亚洲欧美精品综合一区二区三区| 男女免费视频国产| 天天躁日日躁夜夜躁夜夜| 亚洲av日韩精品久久久久久密 | 中文字幕人妻丝袜制服| 亚洲精品av麻豆狂野| 久久久精品免费免费高清| 中国三级夫妇交换| 日日啪夜夜爽| 自拍欧美九色日韩亚洲蝌蚪91| 欧美变态另类bdsm刘玥| av天堂久久9| 中文字幕人妻丝袜一区二区 | 亚洲av日韩精品久久久久久密 | 天美传媒精品一区二区| 欧美最新免费一区二区三区| 热99国产精品久久久久久7| 永久免费av网站大全| 97人妻天天添夜夜摸| 免费高清在线观看视频在线观看| 如何舔出高潮| 国产一区二区在线观看av| 天天躁狠狠躁夜夜躁狠狠躁| 多毛熟女@视频| 尾随美女入室| 国精品久久久久久国模美| 午夜福利在线免费观看网站| 亚洲国产最新在线播放| 九草在线视频观看| 国产成人精品无人区| 黄片无遮挡物在线观看| 不卡av一区二区三区| 啦啦啦 在线观看视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av福利一区| 三上悠亚av全集在线观看| 中文欧美无线码| 亚洲欧美精品综合一区二区三区| 制服丝袜香蕉在线| netflix在线观看网站| 精品国产乱码久久久久久小说| 精品少妇一区二区三区视频日本电影 | 人人妻人人爽人人添夜夜欢视频| 国产亚洲午夜精品一区二区久久| 日韩精品免费视频一区二区三区| 丰满迷人的少妇在线观看| 久久鲁丝午夜福利片| av国产久精品久网站免费入址| 亚洲欧美一区二区三区黑人| 国产亚洲av高清不卡| 亚洲熟女精品中文字幕| 亚洲伊人色综图| 啦啦啦在线观看免费高清www| 国产av一区二区精品久久| 综合色丁香网| 国产av精品麻豆| 国产片特级美女逼逼视频| 免费高清在线观看视频在线观看| 亚洲av成人不卡在线观看播放网 | 叶爱在线成人免费视频播放| 黄色视频在线播放观看不卡| 青草久久国产| 老熟女久久久| 国产男女内射视频| 欧美变态另类bdsm刘玥| av国产精品久久久久影院| 久久国产精品大桥未久av| 在现免费观看毛片| 在线观看免费日韩欧美大片| 99久久人妻综合| 国产男人的电影天堂91| 欧美老熟妇乱子伦牲交| 精品少妇内射三级| 人成视频在线观看免费观看| 蜜桃国产av成人99| 色婷婷av一区二区三区视频| 久久精品国产a三级三级三级| 久久久久精品久久久久真实原创| 少妇人妻久久综合中文| 欧美成人午夜精品| av网站免费在线观看视频| 久久久久网色| 亚洲一码二码三码区别大吗| 久久久久精品性色| 亚洲精品久久久久久婷婷小说| 久久精品久久久久久久性| av免费观看日本| 国产精品蜜桃在线观看| 如日韩欧美国产精品一区二区三区| 亚洲国产欧美日韩在线播放| 在线观看免费午夜福利视频| 男女午夜视频在线观看| 国产成人欧美在线观看 | 丰满少妇做爰视频| 国产成人精品福利久久| 免费观看性生交大片5| 色94色欧美一区二区| 午夜福利视频精品| 五月开心婷婷网| 日韩成人av中文字幕在线观看| 亚洲精品美女久久久久99蜜臀 | 校园人妻丝袜中文字幕| 欧美老熟妇乱子伦牲交| 国产高清国产精品国产三级| 黄片小视频在线播放| 国产成人av激情在线播放| 亚洲男人天堂网一区| 天堂8中文在线网| 亚洲精品美女久久久久99蜜臀 | www.精华液| 在线观看免费视频网站a站| 免费av中文字幕在线| 熟女av电影| 国精品久久久久久国模美| 午夜福利在线免费观看网站| 亚洲四区av| 人成视频在线观看免费观看| 韩国av在线不卡| 成人手机av| 嫩草影视91久久| 色吧在线观看| 男人添女人高潮全过程视频| 在线观看国产h片| 亚洲久久久国产精品| 精品卡一卡二卡四卡免费| 香蕉丝袜av| 国产精品免费视频内射| 十八禁人妻一区二区| 午夜免费男女啪啪视频观看| 亚洲av电影在线观看一区二区三区| 欧美黄色片欧美黄色片| 王馨瑶露胸无遮挡在线观看| 丝袜美腿诱惑在线| 亚洲欧美一区二区三区国产| 日韩制服骚丝袜av| 亚洲成av片中文字幕在线观看| 欧美国产精品va在线观看不卡| 最黄视频免费看| 亚洲精华国产精华液的使用体验| 日韩免费高清中文字幕av| 别揉我奶头~嗯~啊~动态视频 | 九九爱精品视频在线观看| 国产亚洲午夜精品一区二区久久| 最近中文字幕2019免费版| 亚洲精品国产区一区二| 一边亲一边摸免费视频| 国产精品一区二区在线观看99| 亚洲五月色婷婷综合| 少妇的丰满在线观看| 国产探花极品一区二区| 国产精品一区二区精品视频观看| 日本vs欧美在线观看视频| 在线精品无人区一区二区三| 曰老女人黄片|