• <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)格劃分
    亚洲,欧美,日韩| 国产av一区在线观看免费| 如何舔出高潮| 亚洲不卡免费看| 国产成人a区在线观看| 在线国产一区二区在线| 色综合色国产| 国内精品美女久久久久久| 欧美日韩一区二区视频在线观看视频在线 | a级一级毛片免费在线观看| 99久久精品一区二区三区| 乱码一卡2卡4卡精品| 亚洲在线自拍视频| 欧美高清性xxxxhd video| 成人av在线播放网站| 超碰av人人做人人爽久久| 精品乱码久久久久久99久播| 日日撸夜夜添| 91麻豆精品激情在线观看国产| 国产成人91sexporn| 欧美日韩在线观看h| 日本免费一区二区三区高清不卡| 国产熟女欧美一区二区| 午夜福利高清视频| 欧美日韩在线观看h| 人人妻人人澡欧美一区二区| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美三级三区| 精品久久国产蜜桃| 两个人视频免费观看高清| 日本欧美国产在线视频| 亚洲人成网站在线播| 观看免费一级毛片| 精品久久久久久久久久久久久| 啦啦啦韩国在线观看视频| 国产亚洲精品av在线| 国产精品乱码一区二三区的特点| 亚洲经典国产精华液单| 免费看光身美女| 少妇高潮的动态图| 人妻制服诱惑在线中文字幕| 波野结衣二区三区在线| 亚洲第一区二区三区不卡| 日韩一本色道免费dvd| 免费看a级黄色片| 午夜a级毛片| 亚洲熟妇中文字幕五十中出| 女人十人毛片免费观看3o分钟| 夜夜夜夜夜久久久久| 看片在线看免费视频| 午夜福利在线观看免费完整高清在 | 丰满的人妻完整版| 嫩草影院精品99| 男插女下体视频免费在线播放| 九九热线精品视视频播放| 波多野结衣高清作品| 亚洲婷婷狠狠爱综合网| 男人狂女人下面高潮的视频| 国产伦精品一区二区三区四那| 久久久久国产网址| 寂寞人妻少妇视频99o| 一个人看的www免费观看视频| 国产精品乱码一区二三区的特点| 少妇高潮的动态图| 国产高清不卡午夜福利| av.在线天堂| 一级a爱片免费观看的视频| 欧美日韩精品成人综合77777| 日韩欧美三级三区| 亚洲国产欧美人成| 在线播放无遮挡| 久久精品国产鲁丝片午夜精品| 国产aⅴ精品一区二区三区波| 免费观看精品视频网站| 最近最新中文字幕大全电影3| 亚洲图色成人| 欧美日韩在线观看h| 欧美区成人在线视频| 少妇熟女aⅴ在线视频| 搡老熟女国产l中国老女人| 婷婷亚洲欧美| 波野结衣二区三区在线| 国产亚洲精品久久久久久毛片| 成人亚洲欧美一区二区av| 成人av一区二区三区在线看| 久久欧美精品欧美久久欧美| 男女下面进入的视频免费午夜| 免费高清视频大片| 亚洲无线在线观看| 嫩草影院入口| 能在线免费观看的黄片| 女的被弄到高潮叫床怎么办| 国产精品一区二区免费欧美| 免费av观看视频| 免费无遮挡裸体视频| 成人av一区二区三区在线看| 深夜a级毛片| 男女做爰动态图高潮gif福利片| 国产成年人精品一区二区| 亚洲av免费高清在线观看| 黄色欧美视频在线观看| 村上凉子中文字幕在线| 一区二区三区免费毛片| 成人欧美大片| 精华霜和精华液先用哪个| 欧美激情在线99| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 一级毛片电影观看 | 99久久精品热视频| 免费搜索国产男女视频| 免费在线观看成人毛片| 黄片wwwwww| 男女下面进入的视频免费午夜| 天堂√8在线中文| 午夜免费激情av| 直男gayav资源| 小蜜桃在线观看免费完整版高清| 欧美日韩综合久久久久久| 久久精品国产清高在天天线| 精品久久久久久久末码| 日本黄色片子视频| 九九爱精品视频在线观看| 亚洲人成网站高清观看| 亚洲第一区二区三区不卡| 国产精品av视频在线免费观看| 麻豆国产97在线/欧美| 久久6这里有精品| 最后的刺客免费高清国语| 亚洲七黄色美女视频| 村上凉子中文字幕在线| 国产男靠女视频免费网站| 99久久久亚洲精品蜜臀av| 亚洲自拍偷在线| 深夜精品福利| 校园人妻丝袜中文字幕| 亚洲va在线va天堂va国产| 美女高潮的动态| 国产高清三级在线| av专区在线播放| 三级男女做爰猛烈吃奶摸视频| av天堂在线播放| 亚洲自偷自拍三级| 欧美激情在线99| 久久精品夜色国产| 国产综合懂色| 亚洲,欧美,日韩| 久久婷婷人人爽人人干人人爱| 毛片女人毛片| 搞女人的毛片| 天堂√8在线中文| 国产精品一二三区在线看| 日韩国内少妇激情av| 香蕉av资源在线| 乱系列少妇在线播放| 成人永久免费在线观看视频| 高清日韩中文字幕在线| 亚洲一区高清亚洲精品| 欧美日韩一区二区视频在线观看视频在线 | 亚洲成av人片在线播放无| 国产真实伦视频高清在线观看| 91久久精品国产一区二区成人| 观看免费一级毛片| 91久久精品电影网| 成人午夜高清在线视频| 日日摸夜夜添夜夜添小说| 国产午夜精品久久久久久一区二区三区 | 中文字幕免费在线视频6| 99视频精品全部免费 在线| 国产一区二区在线av高清观看| 非洲黑人性xxxx精品又粗又长| 色综合色国产| 国产av不卡久久| 久久久国产成人精品二区| 波野结衣二区三区在线| 久久综合国产亚洲精品| 中文字幕久久专区| 五月伊人婷婷丁香| 热99re8久久精品国产| 国内精品一区二区在线观看| 成年女人永久免费观看视频| 色尼玛亚洲综合影院| 午夜福利高清视频| 久久天躁狠狠躁夜夜2o2o| av专区在线播放| 婷婷色综合大香蕉| 美女 人体艺术 gogo| 国产精品国产高清国产av| 深夜精品福利| 高清毛片免费看| 一区二区三区高清视频在线| 国产精品久久久久久av不卡| 97超碰精品成人国产| 久久精品国产亚洲av香蕉五月| 人妻少妇偷人精品九色| 成人亚洲欧美一区二区av| 91久久精品国产一区二区三区| 97超级碰碰碰精品色视频在线观看| 真人做人爱边吃奶动态| 伦精品一区二区三区| 男人狂女人下面高潮的视频| 亚洲第一电影网av| 国产精品一区二区三区四区免费观看 | 亚洲在线自拍视频| 插逼视频在线观看| 国产av在哪里看| 色综合站精品国产| 日韩高清综合在线| 亚洲国产精品久久男人天堂| 国产精品乱码一区二三区的特点| 无遮挡黄片免费观看| 美女黄网站色视频| 免费人成视频x8x8入口观看| 国产69精品久久久久777片| 国产69精品久久久久777片| 国产激情偷乱视频一区二区| 国产精品野战在线观看| 欧美不卡视频在线免费观看| 亚洲欧美成人综合另类久久久 | 欧美性感艳星| 老师上课跳d突然被开到最大视频| av天堂在线播放| 久久久久九九精品影院| 久久久久免费精品人妻一区二区| 麻豆国产97在线/欧美| 亚洲av不卡在线观看| 国产精品一区二区性色av| 一夜夜www| 99热这里只有是精品50| 深夜a级毛片| 99在线视频只有这里精品首页| 日韩精品中文字幕看吧| 俄罗斯特黄特色一大片| 亚洲丝袜综合中文字幕| 黑人高潮一二区| 欧美成人免费av一区二区三区| 国产爱豆传媒在线观看| 国产精品久久视频播放| 欧美激情久久久久久爽电影| 免费观看人在逋| 中文字幕久久专区| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 22中文网久久字幕| 欧美成人一区二区免费高清观看| 日韩,欧美,国产一区二区三区 | 国产成人91sexporn| 精品无人区乱码1区二区| 国产亚洲精品久久久久久毛片| 一级毛片久久久久久久久女| 小蜜桃在线观看免费完整版高清| 国产成年人精品一区二区| 成人亚洲精品av一区二区| 欧美+亚洲+日韩+国产| 一本久久中文字幕| 欧美激情国产日韩精品一区| 日韩强制内射视频| 天堂av国产一区二区熟女人妻| 国产精品av视频在线免费观看| 在线看三级毛片| av.在线天堂| 成人午夜高清在线视频| 特级一级黄色大片| 亚洲性夜色夜夜综合| 又黄又爽又刺激的免费视频.| 亚洲精品日韩av片在线观看| 久久99热6这里只有精品| 国产免费一级a男人的天堂| 嫩草影院入口| 免费无遮挡裸体视频| 黄色视频,在线免费观看| 国产色婷婷99| 国产精品国产三级国产av玫瑰| 国产精品久久久久久久久免| 日本与韩国留学比较| 亚洲欧美中文字幕日韩二区| av在线蜜桃| 六月丁香七月| 大香蕉久久网| 一夜夜www| 午夜免费男女啪啪视频观看 | 久久亚洲国产成人精品v| 亚洲av成人av| 午夜精品一区二区三区免费看| 日韩精品有码人妻一区| 中出人妻视频一区二区| 国产 一区精品| 午夜激情福利司机影院| 日韩av不卡免费在线播放| 毛片女人毛片| 成年女人看的毛片在线观看| 亚洲精品影视一区二区三区av| av在线亚洲专区| 国产黄a三级三级三级人| 老师上课跳d突然被开到最大视频| 3wmmmm亚洲av在线观看| 欧美中文日本在线观看视频| 村上凉子中文字幕在线| 欧美丝袜亚洲另类| 日日啪夜夜撸| 热99re8久久精品国产| 寂寞人妻少妇视频99o| 97热精品久久久久久| 久久久久国产精品人妻aⅴ院| 乱人视频在线观看| 中国美白少妇内射xxxbb| 国产国拍精品亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 欧美日韩精品成人综合77777| 久久99热6这里只有精品| 国产乱人视频| 亚洲人与动物交配视频| 男女啪啪激烈高潮av片| 美女 人体艺术 gogo| 悠悠久久av| 亚洲欧美清纯卡通| 精品熟女少妇av免费看| 日本欧美国产在线视频| 色哟哟哟哟哟哟| 日韩欧美精品免费久久| av天堂中文字幕网| 成人无遮挡网站| 白带黄色成豆腐渣| 免费人成视频x8x8入口观看| 亚洲经典国产精华液单| 一区二区三区免费毛片| 欧美日韩综合久久久久久| 毛片一级片免费看久久久久| 搞女人的毛片| 搡老妇女老女人老熟妇| 国产成人一区二区在线| 偷拍熟女少妇极品色| 日韩在线高清观看一区二区三区| 黄色视频,在线免费观看| 中文字幕av在线有码专区| 国产欧美日韩精品亚洲av| 久久久a久久爽久久v久久| 色综合亚洲欧美另类图片| 免费黄网站久久成人精品| 日本色播在线视频| 晚上一个人看的免费电影| 国产真实伦视频高清在线观看| 麻豆精品久久久久久蜜桃| 97碰自拍视频| 麻豆乱淫一区二区| 亚洲精品一区av在线观看| 久久人人爽人人片av| 91午夜精品亚洲一区二区三区| 成人av在线播放网站| 国产免费一级a男人的天堂| 国产亚洲av嫩草精品影院| 日本成人三级电影网站| 禁无遮挡网站| 12—13女人毛片做爰片一| 老司机影院成人| 麻豆国产av国片精品| 国产精品一区二区三区四区免费观看 | 啦啦啦观看免费观看视频高清| 国产成年人精品一区二区| 久久久a久久爽久久v久久| 国产精品日韩av在线免费观看| 联通29元200g的流量卡| 97在线视频观看| 国产精品电影一区二区三区| 亚洲无线在线观看| 亚洲性久久影院| 我要看日韩黄色一级片| 国产成人一区二区在线| 两个人视频免费观看高清| 看免费成人av毛片| 美女内射精品一级片tv| 欧美日韩一区二区视频在线观看视频在线 | 激情 狠狠 欧美| 人人妻人人看人人澡| 淫秽高清视频在线观看| 成人永久免费在线观看视频| 亚洲av熟女| 亚洲av五月六月丁香网| 插阴视频在线观看视频| 精品一区二区三区人妻视频| 国产麻豆成人av免费视频| 日韩人妻高清精品专区| 色视频www国产| 国产午夜精品久久久久久一区二区三区 | 欧美日韩国产亚洲二区| 两个人视频免费观看高清| 又黄又爽又免费观看的视频| 国产精品久久久久久av不卡| 全区人妻精品视频| 99在线视频只有这里精品首页| 国产免费一级a男人的天堂| 亚洲欧美中文字幕日韩二区| 日产精品乱码卡一卡2卡三| 最新在线观看一区二区三区| 精品久久久久久久久亚洲| 欧美极品一区二区三区四区| 久久久久久久久大av| 国产aⅴ精品一区二区三区波| 亚洲av美国av| 精品久久久久久久人妻蜜臀av| 91久久精品国产一区二区三区| 亚洲国产日韩欧美精品在线观看| 日韩欧美免费精品| 国产毛片a区久久久久| 日韩av在线大香蕉| 精品不卡国产一区二区三区| 国产片特级美女逼逼视频| 久久热精品热| 亚洲成a人片在线一区二区| 国产精品无大码| 精品日产1卡2卡| 看非洲黑人一级黄片| 亚洲国产精品国产精品| 精品午夜福利在线看| 99精品在免费线老司机午夜| 国产精品精品国产色婷婷| 国产精品不卡视频一区二区| av在线播放精品| .国产精品久久| 插阴视频在线观看视频| 22中文网久久字幕| 国产精品一二三区在线看| 精品人妻视频免费看| 成年av动漫网址| av免费在线看不卡| 免费一级毛片在线播放高清视频| 最近视频中文字幕2019在线8| 色5月婷婷丁香| 一级毛片久久久久久久久女| 网址你懂的国产日韩在线| 亚洲欧美清纯卡通| 亚洲激情五月婷婷啪啪| 亚洲成a人片在线一区二区| 美女大奶头视频| 日本在线视频免费播放| www.色视频.com| 精品日产1卡2卡| 此物有八面人人有两片| 久久这里只有精品中国| 舔av片在线| 午夜精品国产一区二区电影 | 久久国内精品自在自线图片| 婷婷亚洲欧美| 嫩草影院精品99| 啦啦啦观看免费观看视频高清| 久久精品国产亚洲av香蕉五月| 欧美色视频一区免费| 成人鲁丝片一二三区免费| 中文亚洲av片在线观看爽| 插阴视频在线观看视频| 最新在线观看一区二区三区| 身体一侧抽搐| 69av精品久久久久久| 人人妻人人看人人澡| 99在线人妻在线中文字幕| 精品99又大又爽又粗少妇毛片| 国产三级中文精品| 99久久精品国产国产毛片| 亚洲国产精品合色在线| 淫秽高清视频在线观看| 亚洲在线自拍视频| 99久久中文字幕三级久久日本| 欧美成人精品欧美一级黄| 人人妻,人人澡人人爽秒播| 麻豆成人午夜福利视频| av中文乱码字幕在线| 午夜福利成人在线免费观看| 国产成人91sexporn| 免费人成在线观看视频色| 成年女人永久免费观看视频| 国产精品亚洲美女久久久| 99久久无色码亚洲精品果冻| 成人特级黄色片久久久久久久| 国产黄色视频一区二区在线观看 | 国产亚洲精品av在线| 一区二区三区四区激情视频 | 免费电影在线观看免费观看| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| 精品不卡国产一区二区三区| 国产精品一区www在线观看| 久久欧美精品欧美久久欧美| 亚洲精品乱码久久久v下载方式| 婷婷六月久久综合丁香| 特大巨黑吊av在线直播| 国产一级毛片七仙女欲春2| 色在线成人网| 国产中年淑女户外野战色| 免费在线观看影片大全网站| 美女cb高潮喷水在线观看| 日本黄大片高清| 美女cb高潮喷水在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲精品国产av成人精品 | 99九九线精品视频在线观看视频| 亚洲av中文字字幕乱码综合| 床上黄色一级片| 国产精品野战在线观看| 亚洲熟妇熟女久久| 男女那种视频在线观看| 日本精品一区二区三区蜜桃| 99热这里只有是精品50| 午夜影院日韩av| 国产精品嫩草影院av在线观看| 99视频精品全部免费 在线| 插阴视频在线观看视频| 插逼视频在线观看| 久久韩国三级中文字幕| 欧美绝顶高潮抽搐喷水| 国产精品精品国产色婷婷| 尾随美女入室| 国产精品免费一区二区三区在线| 国产单亲对白刺激| 日产精品乱码卡一卡2卡三| 午夜激情欧美在线| 精品人妻偷拍中文字幕| 午夜免费激情av| av福利片在线观看| 天天躁日日操中文字幕| 欧美+日韩+精品| 激情 狠狠 欧美| 欧美潮喷喷水| 99在线人妻在线中文字幕| 极品教师在线视频| 欧美日韩在线观看h| 欧美色欧美亚洲另类二区| 国产精品久久久久久亚洲av鲁大| 老司机影院成人| 亚洲一区二区三区色噜噜| 亚洲五月天丁香| 欧美不卡视频在线免费观看| 亚洲五月天丁香| 黄色欧美视频在线观看| 久久精品影院6| 亚洲av中文字字幕乱码综合| 成年女人永久免费观看视频| 久久久欧美国产精品| 亚洲精品国产成人久久av| 99热精品在线国产| 哪里可以看免费的av片| 精品久久久久久久久久久久久| 啦啦啦韩国在线观看视频| 国产不卡一卡二| 女的被弄到高潮叫床怎么办| 97人妻精品一区二区三区麻豆| 免费看av在线观看网站| 此物有八面人人有两片| 好男人在线观看高清免费视频| 国产一级毛片七仙女欲春2| 欧美xxxx黑人xx丫x性爽| 人人妻人人澡人人爽人人夜夜 | 最新中文字幕久久久久| 特大巨黑吊av在线直播| 九九久久精品国产亚洲av麻豆| 蜜桃久久精品国产亚洲av| 日本 av在线| 国产av一区在线观看免费| 国产精品一及| 一本精品99久久精品77| 国产私拍福利视频在线观看| 男女那种视频在线观看| 国产精品一区二区三区四区久久| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| 波多野结衣巨乳人妻| 九九热线精品视视频播放| 中国美白少妇内射xxxbb| 国产精品爽爽va在线观看网站| 国产精品不卡视频一区二区| 乱人视频在线观看| 一个人看视频在线观看www免费| 1024手机看黄色片| www日本黄色视频网| 免费观看精品视频网站| 色播亚洲综合网| 22中文网久久字幕| 高清毛片免费看| 成人特级av手机在线观看| av天堂中文字幕网| 韩国av在线不卡| 99热只有精品国产| 久久精品国产亚洲av涩爱 | 好男人在线观看高清免费视频| 日韩av不卡免费在线播放| 最近2019中文字幕mv第一页| 嫩草影视91久久| 欧美3d第一页| 美女黄网站色视频| 亚洲,欧美,日韩| 亚洲欧美中文字幕日韩二区| 在线观看一区二区三区| 亚洲成av人片在线播放无| 午夜老司机福利剧场| av在线观看视频网站免费| 久久久a久久爽久久v久久| av视频在线观看入口| ponron亚洲| 中文在线观看免费www的网站| 在线观看美女被高潮喷水网站| 成人毛片a级毛片在线播放| 亚洲国产精品久久男人天堂| 久久久精品欧美日韩精品| 亚洲一区二区三区色噜噜| 成人亚洲精品av一区二区| a级毛片免费高清观看在线播放| h日本视频在线播放| 亚洲一区高清亚洲精品| 国产高清视频在线观看网站| 亚洲av二区三区四区| aaaaa片日本免费| 国产一区亚洲一区在线观看| 久久精品国产亚洲av涩爱 | 亚洲av中文字字幕乱码综合| 精品久久久久久久末码| 亚洲精品日韩av片在线观看|