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

    用于 FLAC3D力學(xué)分析的礦床塊段模型優(yōu)化建模技術(shù)研究

    2010-04-27 05:24:26寇向宇賈明濤王李管
    中國鉬業(yè) 2010年3期
    關(guān)鍵詞:塊段空區(qū)分布圖

    寇向宇,賈明濤,王李管,吳 霞

    (1.中南大學(xué)資源與安全工程學(xué)院,湖南 長沙 410083)

    (2.中南大學(xué)數(shù)字礦山研究中心,湖南 長沙 410083)

    0 引言

    FLAC3D是一種基于連續(xù)介質(zhì)快速拉格朗日差分法的三維空間分析軟件,在采礦與巖土工程中廣泛應(yīng)用。但 FLAC3D的前處理能力不強(qiáng),在建立計(jì)算模型時(shí)主要采用鍵入數(shù)據(jù)/命令行文件方式,雖然可以采用內(nèi)置 Fish語言編程實(shí)現(xiàn)復(fù)雜模型的建立,但其獨(dú)特的源代碼表達(dá)方式、加之實(shí)際采礦工程介質(zhì)和工程形體的復(fù)雜性,由 FLAC3D直接構(gòu)建復(fù)雜模型十分困難[1~2]。

    三維礦床建模軟件可以利用地質(zhì)鉆孔、平剖面圖或其他測量數(shù)據(jù),生成復(fù)雜地質(zhì)和工程對象界(表)面的 T IN模型,并在此模型控制下,將其離散為可以逼近形體表面的規(guī)則六面體單元,從而用于對地質(zhì)體屬性(如品位)等的估值。

    侯恩科等[3~4]介紹了將礦床建模軟件與FLAC3D軟件進(jìn)行耦合的應(yīng)用實(shí)例,但由于六面體單元在 2類軟件中的應(yīng)用目的不同,關(guān)于如何在塊段建模過程中對模型進(jìn)行優(yōu)化以使之滿足力學(xué)分析要求,在這些文獻(xiàn)中并沒有詳細(xì)討論。

    本文從 FLAC3D網(wǎng)格劃分的原則出發(fā),根據(jù)礦床建模軟件 SURPAC構(gòu)建塊段模型的方法,對用于力學(xué)分析的塊段模型構(gòu)建的步驟、原則等進(jìn)行研究,使建立的模型得到優(yōu)化,使其能夠直接作為 FLAC3D分析的前處理模型進(jìn)行模擬計(jì)算。

    1 FLAC3D網(wǎng)格劃分原則及 SURPAC塊段模型構(gòu)建方法

    1.1 FLAC3D網(wǎng)格劃分原則

    有限元網(wǎng)格劃分是進(jìn)行數(shù)值模擬分析至關(guān)重要的一步,它直接影響著后續(xù)數(shù)值計(jì)算分析結(jié)果的精確性。網(wǎng)格劃分涉及單元的形狀及其網(wǎng)格的密度、數(shù)量、質(zhì)量,單元階次,網(wǎng)格分界面和分界點(diǎn),位移協(xié)調(diào)性,網(wǎng)格布局,節(jié)點(diǎn)和單元編號等[5]。

    網(wǎng)格數(shù)量的多少將影響計(jì)算結(jié)果的精度和計(jì)算規(guī)模的大小。網(wǎng)格疏密是指在結(jié)構(gòu)不同部位采用大小不同的網(wǎng)格,這是為了適應(yīng)計(jì)算數(shù)據(jù)的分布特點(diǎn)。在計(jì)算數(shù)據(jù)變化梯度較大的部位(如應(yīng)力集中處),為了較好地反映數(shù)據(jù)變化規(guī)律,需要采用比較密集的網(wǎng)格。而在計(jì)算數(shù)據(jù)變化梯度較小的部位,為減小模型規(guī)模,則應(yīng)劃分相對稀疏的網(wǎng)格。選用高階單元可提高計(jì)算精度,因?yàn)楦唠A單元的曲線或曲面邊界能夠更好地逼近結(jié)構(gòu)的曲線和曲面邊界,且高次插值函數(shù)可更高精度地逼近復(fù)雜場函數(shù),所以當(dāng)結(jié)構(gòu)形狀不規(guī)則、應(yīng)力分布或變形很復(fù)雜時(shí)可以選用高階單元。但高階單元的節(jié)點(diǎn)數(shù)較多,在網(wǎng)格數(shù)量相同的情況下由高階單元組成的模型規(guī)模要大得多,因此在使用時(shí)應(yīng)權(quán)衡考慮計(jì)算精度和時(shí)間。網(wǎng)格質(zhì)量好壞將影響計(jì)算精度。結(jié)構(gòu)中的一些特殊界面和特殊點(diǎn)應(yīng)分為網(wǎng)格邊界或節(jié)點(diǎn)以便定義材料特性、物理特性、載荷和位移約束條件。當(dāng)結(jié)構(gòu)形狀對稱時(shí),其網(wǎng)格也應(yīng)劃分對稱網(wǎng)格,以使模型表現(xiàn)出相應(yīng)的對稱特性(如集中質(zhì)矩陣對稱)。不對稱布局會引起一定誤差。節(jié)點(diǎn)和單元的編號影響結(jié)構(gòu)總剛矩陣的帶寬和波前數(shù),因而影響計(jì)算時(shí)間和存儲容量的大小,因此合理的編號有利于提高計(jì)算速度。

    1.2 SURPAC塊段模型構(gòu)建方法

    塊段模型是一種分割礦體、模擬礦體內(nèi)部特征的有效數(shù)據(jù)模型,就是對目標(biāo)地質(zhì)體進(jìn)行三維柵格化處理[6],把礦體劃分為一系列小長方體單元,其方法是,首先將研究的范圍形成最小包絡(luò)長方體,將其定義為原型;然后根據(jù)地質(zhì)勘探網(wǎng)度、采礦方法、地質(zhì)條件以及地質(zhì)統(tǒng)計(jì)學(xué)等方面對塊段要求來確定單元塊尺寸,以此對原型進(jìn)行三維柵格化;然后對地質(zhì)體進(jìn)行塊段劃分時(shí),應(yīng)選用最合適的、且最具代表性的礦體平剖面圖建立地質(zhì)體表面模型(包括礦體、地層和斷面等模型)對實(shí)體模型進(jìn)行邊界約束,在邊界處進(jìn)一步細(xì)分,以逼近地質(zhì)體的空間形態(tài);每個(gè)小長方體單元的內(nèi)部屬性相同,且該屬性代表了礦體該位置的內(nèi)部特征,這樣所有小長方體單元的屬性變化規(guī)律就是整個(gè)礦體的特征變化規(guī)律。因此,在建立好塊段模型之后,就可以進(jìn)行相應(yīng)的統(tǒng)計(jì)分析工作了。

    2 用于力學(xué)分析的 SURPAC塊段模型優(yōu)化建模技術(shù)

    為使 S URPAC建立的塊段模型能夠直接用于FLAC3D力學(xué)分析,保證力學(xué)分析的效率和精度,根據(jù)前文分析,塊段模型建立時(shí)應(yīng)遵從以下步驟和方法。

    2.1 模型分區(qū)

    根據(jù)網(wǎng)格劃分原則,為了適應(yīng)計(jì)算數(shù)據(jù)的分布特點(diǎn),較好地反映數(shù)據(jù)變化規(guī)律,應(yīng)力可能集中的地方應(yīng)密,其他地方可稀疏一些。網(wǎng)格數(shù)量也要考慮精度和計(jì)算規(guī)模而確定。因此,考慮增加網(wǎng)格的經(jīng)濟(jì)性和分析數(shù)據(jù)的類型,按照不同的劃分方法和力學(xué)參數(shù)建立原始塊段模型,并根據(jù)三維模型進(jìn)行約束,最后,利用編制的接口程序,將塊段模型轉(zhuǎn)換為 FLAC3D分析前處理模型,然后對比不同網(wǎng)格劃分方法的模擬結(jié)果,以確定一個(gè)合適的值。這就是為滿足疏密、密度等的要求,在建模時(shí),可根據(jù)實(shí)際情況分成相關(guān)區(qū)段,外圍區(qū)段單元大、往里逐漸減小,只對最里面有形狀變化的區(qū)段細(xì)分。這樣既保證計(jì)算精度,又保證單元數(shù)量不會太多,且虛擬節(jié)點(diǎn)不會很多,找不全的問題,進(jìn)而提高了分析效率,實(shí)現(xiàn)了網(wǎng)格劃分及其優(yōu)化。

    2.2 邊界單元細(xì)分

    邊界單元細(xì)分一是為逼近幾何邊界、二是達(dá)到疏密網(wǎng)格的目的,從而反映應(yīng)力集中區(qū)的應(yīng)力變化。在通過 SURPAC軟件輸出網(wǎng)格時(shí),要進(jìn)行塊段模型的建立,在建立的過程中,塊段模型選取地質(zhì)體線框模型進(jìn)行約束,塊段模型的單元塊在邊界處將自動按照一定的準(zhǔn)則發(fā)生細(xì)分(變塊技術(shù))。同時(shí)在建模時(shí)細(xì)分次數(shù)也可超過 2次以上,導(dǎo)入 FLAC3D中使一個(gè)單元跟相鄰單元的邊長之比小于 1∶1或 1∶2,這樣對復(fù)雜形態(tài)邊界的逼近能力更高,效果更好。另外,在細(xì)分后導(dǎo)入 FLAC3D必須利用Attach FAce或 Gen merge命令找出并處理這些虛節(jié)點(diǎn)。因?yàn)楫?dāng)模型較大時(shí),如細(xì)分次數(shù)過多,會導(dǎo)致接觸面虛節(jié)點(diǎn)數(shù)量多,FLAC3D可能找不全這些點(diǎn),導(dǎo)致計(jì)算結(jié)果出錯(cuò)。

    2.3 結(jié)果驗(yàn)證

    為了驗(yàn)證這種網(wǎng)格劃分實(shí)現(xiàn)前處理模型優(yōu)化方法的正確性,根據(jù)二維彈性力學(xué)的原理[7],設(shè)計(jì)了一個(gè)邊長 48 m×48 m的板中有一個(gè)半徑為 3 m圓孔,在左右兩邊受均布拉力,大小為 1 MPa,然后根據(jù)對模型不同的細(xì)分方法得出疏密不同的單元體網(wǎng)格,進(jìn)行分析,再與理論值作比較,最終得出最優(yōu)的網(wǎng)格劃分方法,進(jìn)而也驗(yàn)證了這種優(yōu)化方法在力學(xué)分析中的重要性和正確性。所以以此方法做出理論應(yīng)力分布圖(11200Zones and 22402 Grid),遠(yuǎn)近疏密一致的應(yīng)力分布圖 (9988zones and 15448 grid points),疏密程度不一致的應(yīng)力分布圖 (9256zones and 18544 grid points),比較應(yīng)力的主要分布,顯示出不同網(wǎng)格劃分對數(shù)值模擬效果的影響,如圖 1~圖 3所示。圖中顯示,理論上的應(yīng)力應(yīng)該集中在中心的兩側(cè),而在不同位置的不同疏密程度模擬出來的效果不同,并不是網(wǎng)格越密效果越好,也不是網(wǎng)格越疏效果越好,只有在應(yīng)力相對應(yīng)該集中的地方密一些,其他地方疏一些,圖 2在兩側(cè)時(shí)應(yīng)力出現(xiàn)了發(fā)散,而圖 3模擬的效果最接近于理論值,效果最好,因此得出了網(wǎng)格劃分的最優(yōu)選擇,圖 4為網(wǎng)格劃分的顯示,進(jìn)而也得到了前處理模型的優(yōu)化。

    圖 1 理論應(yīng)力分布圖

    圖 2 疏密程度一致應(yīng)力分布圖

    圖 3 疏密程度不一致應(yīng)力分布圖

    圖 4 最優(yōu)網(wǎng)格劃分效果圖

    3 工程應(yīng)用

    3.1 工程概況

    云南某銅礦經(jīng)過數(shù)十年的露天及地下大規(guī)模工業(yè)開采,井下采空區(qū)越來越多、越來越大,且許多早期、淺部空區(qū)已經(jīng)發(fā)生垮塌,多個(gè)中段上下相鄰或同一水平相鄰的空區(qū)發(fā)生貫通,甚至部分已經(jīng)貫通地表,對地表各種工業(yè)和社會活動,以及生態(tài)環(huán)境造成一定的威脅,存在較大的隱患。如果沒有有效的采空區(qū)穩(wěn)定性及地表移動規(guī)律的評價(jià)研究,則采空區(qū)垮落不僅加大礦石回采的貧化與損失,也必然對礦井及礦區(qū)相關(guān)工農(nóng)業(yè)的安全生產(chǎn)造成嚴(yán)重的威脅。為保證礦山有序的正常安全生產(chǎn)以及后續(xù)開采的安全,對該礦采空區(qū)穩(wěn)定性及由此引發(fā)的地表移動問題進(jìn)行研究具有非常重要的現(xiàn)實(shí)意義。因此,運(yùn)用FLAC3D數(shù)值方法模擬了影響采場穩(wěn)定性的空區(qū)圍巖應(yīng)力及地表位移變化規(guī)律,并在模型建立時(shí)以不同網(wǎng)格劃分方法進(jìn)行對比,實(shí)現(xiàn)前處理模型的優(yōu)化并進(jìn)行數(shù)值模擬,進(jìn)而有效地提高了采空區(qū)穩(wěn)定性分析和評價(jià)結(jié)果的準(zhǔn)確性和可靠性。

    3.2 力學(xué)模擬及其分析

    在 SURPAC中,在地表三維模型的約束下,建立原始塊段模型。其次,在各中段空區(qū)三維表面模型的約束下,建立相應(yīng)的空區(qū)塊段模型,然后,把所有塊段模型合并,形成一個(gè)統(tǒng)一的模型,該模型內(nèi)單元塊數(shù)目為 638762塊。最后,利用編制的接口程序,將 SURPAC塊段模型轉(zhuǎn)換為 FLAC3D分析前處理模型,該模型節(jié)點(diǎn)數(shù)為 810780個(gè)。

    在前處理模型中,將各采空區(qū)所包含的單元設(shè)置為不同的組 (Group),賦力學(xué)參數(shù)、初始地應(yīng)力、邊界條件,就可以在分析過程中,通過指定某個(gè)組的移出,以達(dá)到模擬開挖,考察在采空區(qū)形成后,周邊圍巖應(yīng)力及位移分布情況的目的,最后得出結(jié)果。本文選取礦區(qū)某中段的空區(qū)作為實(shí)例進(jìn)行分析比較,按照相關(guān)力學(xué)參數(shù)進(jìn)行賦值,做出空區(qū)圍巖應(yīng)力分布圖和地表位移分布圖(圖 5、圖 6)。

    圖 5 中段開采后最小主應(yīng)力空間分布圖

    圖 6 中段開采后地表位移圖

    圖 5為第一步開挖后模型內(nèi)的最小主應(yīng)力分布區(qū)及其大小。由圖 5可見,在采空區(qū)周邊圍巖中出現(xiàn)了拉應(yīng)力,多在0.3MPa以下,少數(shù)單元拉應(yīng)力較大,接近但小于 0.8 MPa。由圖 6可以看出,地表位移量整體上不是很大,在 1.5 cm以內(nèi);最大位移發(fā)生在露天坑坑底中心位置,隨著離坑底中心漸遠(yuǎn),整體位移逐漸減??;然后按采場模擬開挖順序?qū)Ω髦卸芜M(jìn)行開挖,直至結(jié)束,并且分析各中段開挖后模型內(nèi)的最小主應(yīng)力分布區(qū)及其大小,以及采場中各應(yīng)力的變化情況以及開采后的地表位移變化,實(shí)現(xiàn)對該礦區(qū)采空區(qū)穩(wěn)定性的準(zhǔn)確模擬,以制定措施對井下空區(qū)進(jìn)行處理,排除礦區(qū)生產(chǎn)的安全隱患,保障礦區(qū)的正常生活和作業(yè)。

    4 結(jié)論

    為了實(shí)現(xiàn)數(shù)值模擬的效果更加準(zhǔn)確,研究了有限元網(wǎng)格劃分原則和 SURPAC構(gòu)建塊段模型的方法,然后建立了最優(yōu)的 FLAC3D前處理模型,并構(gòu)建了一個(gè)實(shí)例,依據(jù)彈性力學(xué)的原理,判定其有效性和準(zhǔn)確性,同時(shí)也驗(yàn)證此種前處理模型優(yōu)化方法的正確性,從而實(shí)現(xiàn)了 FLAC3D的前處理模型優(yōu)化。以此在云南某銅礦進(jìn)行設(shè)計(jì)實(shí)踐,通過按照礦區(qū)實(shí)際開采順序進(jìn)行模擬,做出采空區(qū)的應(yīng)力分布和地表的位移移動情況,判斷空區(qū)的穩(wěn)定性和地表生產(chǎn)生活的安全性,并制定一系列措施進(jìn)行治理,排除了礦區(qū)生產(chǎn)的安全隱患,保障了礦區(qū)的正常生活和安全作業(yè),在實(shí)際中驗(yàn)證了此模擬的有效性和精確性,同時(shí)也驗(yàn)證了該前處理模型優(yōu)化方法的可行性和必要性。

    [1] 鄧紅衛(wèi),朱和玲,周科平,等.基于 FLAC3D數(shù)值模擬的前后處理優(yōu)化研究 [J].礦業(yè)研究與開發(fā),2008,28 (2):60-62.

    [2] 胡斌,張倬元,黃潤秋,等.FLAC3D前處理程序的開發(fā)及仿真效果檢驗(yàn) [J].巖石力學(xué)與工程學(xué)報(bào),2002, 21(9):1387-1391.

    [3] 劉科偉,李夕兵,宮鳳強(qiáng),等.基于 CALS及 Surpac-FLAC3D耦合技術(shù)的復(fù)雜空區(qū)穩(wěn)定性分析[J].巖石力學(xué)與工程學(xué)報(bào),2008,27(9):1924-1931.

    [4] 侯恩科,吳立新,李建民,等.三維地學(xué)模擬與數(shù)值模擬的耦合方法研究[J].煤炭學(xué)報(bào),2002,27(4):388 -392.

    [5] 杜平安.有限元網(wǎng)格劃分的基本原則[J].機(jī)械設(shè)計(jì)與制造,2000:34-36.

    [6] 馮興隆,王李管,畢林,等.基于三維模擬技術(shù)的礦巖可崩性評價(jià)[J].煤炭學(xué)報(bào),2008,33(9):971-976.

    [7] 徐芝綸.彈性力學(xué)簡明教程 [M].人民教育出版社, 1983:200-203.

    猜你喜歡
    塊段空區(qū)分布圖
    入駐機(jī)構(gòu)樓層分布圖
    分憂(2024年3期)2024-01-01 00:00:00
    資源儲量估算中地質(zhì)塊段法與剖面法對比分析
    貴州十大地質(zhì)公園分布圖
    空區(qū)群結(jié)構(gòu)力學(xué)效應(yīng)模擬分析
    淺析露天礦儲量計(jì)算方法
    中國癌癥分布圖
    左右江水沖石器采集分布圖
    寶藏(2017年6期)2017-07-20 10:01:01
    地質(zhì)塊段法在卡而卻卡銅礦區(qū)資源量估算中應(yīng)用及注意事項(xiàng)
    弓長嶺露天礦含空區(qū)邊坡穩(wěn)定性研究
    金屬礦山(2015年11期)2015-03-20 07:54:01
    地震空區(qū)法在大地震危險(xiǎn)性初判中的應(yīng)用——以青藏高原東南緣為例
    91精品国产国语对白视频| 丁香六月天网| 老司机影院成人| 一级毛片电影观看| 欧美日韩视频精品一区| 成人亚洲精品一区在线观看| 精品人妻一区二区三区麻豆| 女性被躁到高潮视频| 制服人妻中文乱码| 少妇被粗大猛烈的视频| 久久久久久久久久人人人人人人| 国产精品久久久av美女十八| 亚洲欧美中文字幕日韩二区| 久久影院123| 成人国产麻豆网| 最近最新中文字幕大全免费视频 | 成人黄色视频免费在线看| 26uuu在线亚洲综合色| 一级,二级,三级黄色视频| 黄色配什么色好看| 亚洲人成77777在线视频| 久久久久国产网址| 久久亚洲国产成人精品v| 精品亚洲乱码少妇综合久久| 国产日韩欧美在线精品| 久久久a久久爽久久v久久| 久久精品夜色国产| 精品少妇内射三级| 免费av中文字幕在线| 久久精品国产自在天天线| 永久免费av网站大全| 黄色配什么色好看| 国产淫语在线视频| 18禁观看日本| 九色亚洲精品在线播放| 又黄又粗又硬又大视频| 久久 成人 亚洲| 亚洲精品国产一区二区精华液| 日韩中字成人| 久久影院123| 一本—道久久a久久精品蜜桃钙片| 超色免费av| 免费观看在线日韩| 深夜精品福利| 国产精品麻豆人妻色哟哟久久| 国产亚洲一区二区精品| 我要看黄色一级片免费的| 观看美女的网站| 2018国产大陆天天弄谢| 最新中文字幕久久久久| 亚洲av欧美aⅴ国产| 看十八女毛片水多多多| 日韩欧美一区视频在线观看| av线在线观看网站| 亚洲精品第二区| 春色校园在线视频观看| 国产综合精华液| 2018国产大陆天天弄谢| 亚洲美女视频黄频| 亚洲精品aⅴ在线观看| 18禁国产床啪视频网站| 涩涩av久久男人的天堂| 国产老妇伦熟女老妇高清| 免费女性裸体啪啪无遮挡网站| 热99国产精品久久久久久7| 国产午夜精品一二区理论片| 日日撸夜夜添| 国产精品偷伦视频观看了| 国产成人精品福利久久| 久久精品aⅴ一区二区三区四区 | av有码第一页| 夜夜骑夜夜射夜夜干| 成人影院久久| 99久久中文字幕三级久久日本| 国产av码专区亚洲av| 亚洲在久久综合| 热99久久久久精品小说推荐| 国产xxxxx性猛交| 国产在线一区二区三区精| 飞空精品影院首页| 黄网站色视频无遮挡免费观看| 久久久久精品性色| 欧美人与性动交α欧美软件| 精品少妇内射三级| 黄色怎么调成土黄色| 99九九在线精品视频| 免费黄网站久久成人精品| 国产伦理片在线播放av一区| 80岁老熟妇乱子伦牲交| 久久久久精品人妻al黑| 肉色欧美久久久久久久蜜桃| 亚洲一区二区三区欧美精品| 婷婷色综合大香蕉| 777米奇影视久久| 日韩一卡2卡3卡4卡2021年| 极品人妻少妇av视频| 国产97色在线日韩免费| 亚洲av电影在线进入| 国产在线一区二区三区精| 三级国产精品片| 男人爽女人下面视频在线观看| 精品少妇内射三级| 国产日韩一区二区三区精品不卡| 成年女人毛片免费观看观看9 | 老汉色av国产亚洲站长工具| 日韩一区二区视频免费看| 高清在线视频一区二区三区| 亚洲成人一二三区av| 亚洲欧美一区二区三区国产| 一本久久精品| 91精品国产国语对白视频| 青春草亚洲视频在线观看| 欧美精品国产亚洲| 亚洲综合色惰| 中文字幕人妻熟女乱码| 汤姆久久久久久久影院中文字幕| 寂寞人妻少妇视频99o| 亚洲av电影在线进入| 伦精品一区二区三区| 久久久精品免费免费高清| 男女午夜视频在线观看| 久久久久久久大尺度免费视频| 高清欧美精品videossex| 免费观看无遮挡的男女| 精品第一国产精品| 亚洲欧洲日产国产| 一本色道久久久久久精品综合| 国产黄频视频在线观看| 精品人妻偷拍中文字幕| 亚洲欧美精品自产自拍| 亚洲成人av在线免费| 国产在线免费精品| 精品少妇黑人巨大在线播放| 久久久欧美国产精品| 久久婷婷青草| 国产精品麻豆人妻色哟哟久久| 伊人久久大香线蕉亚洲五| 国产成人欧美| 超碰97精品在线观看| 99热全是精品| 水蜜桃什么品种好| 在线 av 中文字幕| 国产成人免费无遮挡视频| 人人妻人人添人人爽欧美一区卜| 亚洲欧洲日产国产| 欧美精品国产亚洲| 中文字幕色久视频| 一本色道久久久久久精品综合| 少妇人妻久久综合中文| 日韩精品有码人妻一区| 婷婷色综合大香蕉| 免费观看a级毛片全部| 国产一区二区三区av在线| 日韩人妻精品一区2区三区| 国产成人精品无人区| 亚洲欧美一区二区三区黑人 | xxxhd国产人妻xxx| 欧美日韩视频高清一区二区三区二| 欧美少妇被猛烈插入视频| 国产精品三级大全| 男女边吃奶边做爰视频| 中国三级夫妇交换| 亚洲国产av影院在线观看| 成人黄色视频免费在线看| 亚洲精品,欧美精品| 欧美亚洲 丝袜 人妻 在线| 日韩av不卡免费在线播放| 精品福利永久在线观看| av网站免费在线观看视频| 最近最新中文字幕大全免费视频 | 高清不卡的av网站| 最近中文字幕高清免费大全6| 亚洲精品久久午夜乱码| 黄色怎么调成土黄色| 精品少妇久久久久久888优播| 日本av手机在线免费观看| 日韩成人av中文字幕在线观看| 青青草视频在线视频观看| 国产成人av激情在线播放| 黄色怎么调成土黄色| 久久97久久精品| 久久 成人 亚洲| 亚洲成人手机| 美国免费a级毛片| av不卡在线播放| 亚洲国产色片| 人妻少妇偷人精品九色| 97精品久久久久久久久久精品| 秋霞在线观看毛片| 国产亚洲一区二区精品| 蜜桃国产av成人99| 午夜精品国产一区二区电影| 亚洲美女搞黄在线观看| 新久久久久国产一级毛片| 国产一区有黄有色的免费视频| 国产日韩欧美在线精品| 久久精品熟女亚洲av麻豆精品| 考比视频在线观看| 国产伦理片在线播放av一区| 高清视频免费观看一区二区| www.av在线官网国产| 90打野战视频偷拍视频| 国产精品av久久久久免费| 成年美女黄网站色视频大全免费| 成人毛片60女人毛片免费| 天天躁夜夜躁狠狠久久av| 国产黄频视频在线观看| 欧美日韩一级在线毛片| 亚洲综合色惰| 又粗又硬又长又爽又黄的视频| 黑人猛操日本美女一级片| 久久久久久久久免费视频了| 久久精品久久精品一区二区三区| 91精品伊人久久大香线蕉| 最近手机中文字幕大全| a 毛片基地| 免费看av在线观看网站| 自线自在国产av| 亚洲欧美色中文字幕在线| 久久久精品免费免费高清| 伦精品一区二区三区| 国产免费现黄频在线看| 欧美xxⅹ黑人| 亚洲av国产av综合av卡| 欧美人与性动交α欧美精品济南到 | 久久久久久久亚洲中文字幕| 日韩视频在线欧美| 熟女av电影| 国产免费现黄频在线看| 国产免费福利视频在线观看| 亚洲精品日韩在线中文字幕| 中文欧美无线码| 精品国产乱码久久久久久小说| 女人精品久久久久毛片| 大片电影免费在线观看免费| 久久热在线av| 丝袜美足系列| 欧美日韩av久久| 超碰97精品在线观看| 国产乱人偷精品视频| 精品一品国产午夜福利视频| 深夜精品福利| 精品国产露脸久久av麻豆| 久久免费观看电影| 五月伊人婷婷丁香| 国产97色在线日韩免费| 多毛熟女@视频| 日日摸夜夜添夜夜爱| 在现免费观看毛片| 久久久久久免费高清国产稀缺| 亚洲伊人色综图| 人妻 亚洲 视频| 黄色一级大片看看| 满18在线观看网站| 两性夫妻黄色片| 成人国产av品久久久| 久久久久久人妻| 日韩电影二区| 亚洲国产日韩一区二区| 亚洲第一区二区三区不卡| 国产黄频视频在线观看| 桃花免费在线播放| 制服人妻中文乱码| 精品久久久精品久久久| 欧美日韩精品网址| 美女视频免费永久观看网站| av在线app专区| 欧美变态另类bdsm刘玥| 成人亚洲精品一区在线观看| 精品人妻偷拍中文字幕| 视频区图区小说| 免费av中文字幕在线| 午夜福利视频精品| 亚洲欧美一区二区三区国产| 国产爽快片一区二区三区| 在线观看免费视频网站a站| 美女高潮到喷水免费观看| 成人毛片a级毛片在线播放| 看非洲黑人一级黄片| 国产激情久久老熟女| 最近中文字幕高清免费大全6| 99久久精品国产国产毛片| 亚洲三区欧美一区| 精品亚洲成国产av| 欧美激情 高清一区二区三区| 久久国产精品大桥未久av| 麻豆av在线久日| 久久精品久久久久久噜噜老黄| 久久久国产一区二区| av在线播放精品| 久久久久精品性色| 亚洲精品国产av蜜桃| h视频一区二区三区| 看十八女毛片水多多多| 一级片免费观看大全| 性少妇av在线| 人人澡人人妻人| 老司机亚洲免费影院| 国产成人欧美| 一级爰片在线观看| 欧美黄色片欧美黄色片| 一区福利在线观看| 99国产精品免费福利视频| 国产精品国产av在线观看| 夫妻性生交免费视频一级片| 丝袜人妻中文字幕| 最近最新中文字幕免费大全7| 黄色怎么调成土黄色| 亚洲一区中文字幕在线| 精品99又大又爽又粗少妇毛片| 80岁老熟妇乱子伦牲交| 老司机影院毛片| 国产精品一国产av| 亚洲内射少妇av| 99久久综合免费| 免费看av在线观看网站| 国产色婷婷99| 色94色欧美一区二区| 国产精品无大码| 国产男人的电影天堂91| 久久国产亚洲av麻豆专区| 久久久久国产精品人妻一区二区| 精品99又大又爽又粗少妇毛片| 日韩电影二区| 久久国产精品大桥未久av| 久久精品亚洲av国产电影网| 亚洲国产欧美在线一区| 久久 成人 亚洲| 国产日韩欧美视频二区| 天天影视国产精品| 女人久久www免费人成看片| 一区二区三区精品91| 亚洲av福利一区| 好男人视频免费观看在线| 丝袜在线中文字幕| 久久久a久久爽久久v久久| 久久国产精品男人的天堂亚洲| 18+在线观看网站| 亚洲av日韩在线播放| 精品一区二区三区四区五区乱码 | 秋霞在线观看毛片| 巨乳人妻的诱惑在线观看| 波野结衣二区三区在线| 在线观看免费日韩欧美大片| √禁漫天堂资源中文www| 美国免费a级毛片| 精品99又大又爽又粗少妇毛片| 亚洲国产精品999| 一级片免费观看大全| 久久久久国产一级毛片高清牌| 亚洲成人手机| 久热久热在线精品观看| 亚洲视频免费观看视频| 日韩精品有码人妻一区| 91久久精品国产一区二区三区| 精品少妇内射三级| 天天影视国产精品| 久久精品亚洲av国产电影网| 又黄又粗又硬又大视频| 国产精品国产av在线观看| 国产成人91sexporn| 九色亚洲精品在线播放| 人人妻人人澡人人爽人人夜夜| 国产精品免费大片| 欧美亚洲日本最大视频资源| 午夜免费鲁丝| 日本-黄色视频高清免费观看| 精品久久久久久电影网| 久久精品久久精品一区二区三区| 人妻系列 视频| 国产伦理片在线播放av一区| 亚洲av中文av极速乱| 国产激情久久老熟女| 欧美日韩av久久| 亚洲人成网站在线观看播放| 天天操日日干夜夜撸| 亚洲国产av新网站| 丰满少妇做爰视频| 欧美激情 高清一区二区三区| 亚洲美女黄色视频免费看| av网站免费在线观看视频| 男女啪啪激烈高潮av片| av免费观看日本| 免费人妻精品一区二区三区视频| 成人免费观看视频高清| 亚洲精品自拍成人| 9色porny在线观看| 久久久久久人人人人人| 国产成人精品在线电影| 色吧在线观看| 曰老女人黄片| 亚洲,欧美精品.| 视频在线观看一区二区三区| 9色porny在线观看| 18禁观看日本| 你懂的网址亚洲精品在线观看| av网站在线播放免费| 国产精品久久久久成人av| 午夜免费观看性视频| 国产精品偷伦视频观看了| 亚洲,一卡二卡三卡| 男女午夜视频在线观看| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| 丝袜在线中文字幕| 亚洲成人av在线免费| 久久久国产一区二区| 三级国产精品片| 女性生殖器流出的白浆| 中文字幕色久视频| 午夜激情av网站| 99国产精品免费福利视频| 日韩制服骚丝袜av| 男人爽女人下面视频在线观看| 熟女av电影| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区久久| 男人添女人高潮全过程视频| 午夜日本视频在线| 一二三四中文在线观看免费高清| 极品少妇高潮喷水抽搐| 国产福利在线免费观看视频| 成人18禁高潮啪啪吃奶动态图| 18禁观看日本| 满18在线观看网站| 黄片播放在线免费| 极品少妇高潮喷水抽搐| 老鸭窝网址在线观看| 日本爱情动作片www.在线观看| videos熟女内射| av网站免费在线观看视频| 男女高潮啪啪啪动态图| 午夜精品国产一区二区电影| 在线看a的网站| 国产综合精华液| 久久久a久久爽久久v久久| 亚洲av中文av极速乱| 美女福利国产在线| 国产免费又黄又爽又色| 亚洲美女黄色视频免费看| 国产免费一区二区三区四区乱码| 秋霞伦理黄片| 久久精品久久久久久噜噜老黄| 国产精品嫩草影院av在线观看| 宅男免费午夜| 久久久精品国产亚洲av高清涩受| 色婷婷久久久亚洲欧美| 人妻一区二区av| 女人精品久久久久毛片| 欧美成人午夜精品| 老司机影院毛片| 夫妻午夜视频| 纯流量卡能插随身wifi吗| 欧美xxⅹ黑人| 校园人妻丝袜中文字幕| 午夜福利乱码中文字幕| 人成视频在线观看免费观看| 久久韩国三级中文字幕| kizo精华| 中文字幕最新亚洲高清| 人妻 亚洲 视频| 国产欧美亚洲国产| 免费在线观看完整版高清| 成人国语在线视频| 中文字幕亚洲精品专区| 丝袜在线中文字幕| 搡老乐熟女国产| 久久久国产欧美日韩av| 精品国产一区二区久久| 亚洲av综合色区一区| 亚洲精品久久久久久婷婷小说| 亚洲欧美一区二区三区国产| 天天躁夜夜躁狠狠久久av| av女优亚洲男人天堂| 免费观看在线日韩| 午夜福利影视在线免费观看| 免费女性裸体啪啪无遮挡网站| 高清欧美精品videossex| 国产成人精品久久二区二区91 | 久久精品久久久久久噜噜老黄| 在线亚洲精品国产二区图片欧美| 免费在线观看黄色视频的| 黑人巨大精品欧美一区二区蜜桃| 日本vs欧美在线观看视频| 亚洲伊人色综图| 又大又黄又爽视频免费| 欧美成人午夜免费资源| 日韩精品有码人妻一区| 亚洲伊人色综图| 免费不卡的大黄色大毛片视频在线观看| 欧美精品高潮呻吟av久久| 菩萨蛮人人尽说江南好唐韦庄| 免费黄网站久久成人精品| 日本vs欧美在线观看视频| av在线老鸭窝| 在线观看三级黄色| 久久久精品免费免费高清| 欧美另类一区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美一区二区三区久久| 一级毛片黄色毛片免费观看视频| 亚洲四区av| 母亲3免费完整高清在线观看 | 国产一区二区三区综合在线观看| xxxhd国产人妻xxx| 伦理电影大哥的女人| 午夜福利视频精品| 国产高清不卡午夜福利| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av高清一级| 欧美国产精品va在线观看不卡| 成人毛片a级毛片在线播放| 人妻系列 视频| 欧美人与性动交α欧美精品济南到 | 亚洲国产精品国产精品| 久久久久视频综合| 一区福利在线观看| √禁漫天堂资源中文www| 尾随美女入室| 欧美日韩视频高清一区二区三区二| 日韩在线高清观看一区二区三区| 一级片'在线观看视频| 午夜免费观看性视频| 97在线视频观看| 欧美+日韩+精品| 丝袜美足系列| 天堂中文最新版在线下载| 曰老女人黄片| 亚洲精品aⅴ在线观看| 水蜜桃什么品种好| 这个男人来自地球电影免费观看 | 激情视频va一区二区三区| 久久精品亚洲av国产电影网| 国产无遮挡羞羞视频在线观看| 日本午夜av视频| 国产爽快片一区二区三区| 99热网站在线观看| 涩涩av久久男人的天堂| 夫妻午夜视频| 亚洲av电影在线进入| 少妇人妻精品综合一区二区| 性少妇av在线| 国产成人一区二区在线| 日本欧美国产在线视频| 欧美人与性动交α欧美精品济南到 | 亚洲人成77777在线视频| 9热在线视频观看99| 亚洲欧洲国产日韩| 免费看不卡的av| 国语对白做爰xxxⅹ性视频网站| 亚洲精品久久久久久婷婷小说| av线在线观看网站| 伦理电影免费视频| 国产有黄有色有爽视频| 久久久精品区二区三区| 一级毛片 在线播放| 亚洲成av片中文字幕在线观看 | 久久鲁丝午夜福利片| 五月天丁香电影| 久久精品久久久久久噜噜老黄| 校园人妻丝袜中文字幕| 成年美女黄网站色视频大全免费| 欧美成人精品欧美一级黄| 欧美精品国产亚洲| 日韩熟女老妇一区二区性免费视频| 午夜激情av网站| 亚洲伊人久久精品综合| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久毛片免费看一区二区三区| 午夜老司机福利剧场| 两性夫妻黄色片| 亚洲三级黄色毛片| 大片免费播放器 马上看| 欧美日韩成人在线一区二区| 国产日韩欧美在线精品| 国产成人免费无遮挡视频| 午夜日本视频在线| 国产亚洲av片在线观看秒播厂| 日韩成人av中文字幕在线观看| 交换朋友夫妻互换小说| 精品亚洲乱码少妇综合久久| 国产综合精华液| 中文天堂在线官网| 十八禁高潮呻吟视频| 综合色丁香网| 美女大奶头黄色视频| 十八禁高潮呻吟视频| 日韩成人av中文字幕在线观看| 欧美精品人与动牲交sv欧美| 亚洲少妇的诱惑av| 日本爱情动作片www.在线观看| a级毛片在线看网站| videos熟女内射| 久久精品久久精品一区二区三区| 日韩一本色道免费dvd| www.自偷自拍.com| 亚洲国产精品一区三区| 亚洲一码二码三码区别大吗| 乱人伦中国视频| 校园人妻丝袜中文字幕| 男女啪啪激烈高潮av片| av在线老鸭窝| 久久韩国三级中文字幕| 国产精品成人在线| h视频一区二区三区| 美国免费a级毛片| 精品国产乱码久久久久久小说| kizo精华| 男人爽女人下面视频在线观看| 日韩一卡2卡3卡4卡2021年| 久久精品国产亚洲av涩爱| 一级a爱视频在线免费观看| 中文字幕色久视频| 成人亚洲精品一区在线观看| 婷婷色av中文字幕|