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

    基于最小耗能原理的巖石損傷本構(gòu)模型研究

    2018-09-07 09:56:00孫夢(mèng)成徐衛(wèi)亞王蘇生王如賓王偉
    關(guān)鍵詞:角礫熔巖巖樣

    孫夢(mèng)成,徐衛(wèi)亞,王蘇生,王如賓,王偉

    ?

    基于最小耗能原理的巖石損傷本構(gòu)模型研究

    孫夢(mèng)成1, 2,徐衛(wèi)亞1, 2,王蘇生1, 2,王如賓1, 2,王偉1, 2

    (1. 河海大學(xué) 巖土力學(xué)與堤壩工程教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京,210098;2. 河海大學(xué) 巖土工程科學(xué)研究所,江蘇 南京,210098)

    在連續(xù)損傷理論框架下引入最小耗能原理,將巖石統(tǒng)一能量屈服準(zhǔn)則作為耗能約束條件,基于各向同性基本假定構(gòu)建一種新型損傷本構(gòu)模型。假定以巖石初始屈服點(diǎn)作為損傷閾值,模型考慮損傷閾值對(duì)損傷演化規(guī)律的影響;推導(dǎo)巖石三軸壓縮強(qiáng)度和變形參量與模型參數(shù)的理論關(guān)系。為驗(yàn)證模型的合理性,對(duì)角礫熔巖進(jìn)行不同圍壓下的常規(guī)三軸壓縮試驗(yàn),并分析破壞過(guò)程中巖石變形各階段的能量演化特征。研究結(jié)果表明:該模型能較好地反映巖石材料在復(fù)雜應(yīng)力條件下的非線性力學(xué)行為,且考慮損傷閾值影響的模型對(duì)巖石后屈服段的擬合精度更高,對(duì)復(fù)雜應(yīng)力條件下的巖石工程安全分析具有參考價(jià)值。

    巖石力學(xué);最小耗能原理;巖石統(tǒng)一能量屈服準(zhǔn)則;損傷閾值;損傷演化方程

    漫長(zhǎng)的地質(zhì)演變過(guò)程及復(fù)雜的三向應(yīng)力環(huán)境使巖石內(nèi)部產(chǎn)生大量非連續(xù)且形態(tài)不規(guī)則的隱裂隙、微裂紋、微孔洞及節(jié)理等微缺陷,微缺陷在不同應(yīng)力條件下的力學(xué)響應(yīng)決定其顯著的宏觀非線性力學(xué)特性[1]。大量巖石高邊坡及大型地下洞室開(kāi)挖工程實(shí)踐表明,巖石的力學(xué)特性對(duì)巖土結(jié)構(gòu)的安全設(shè)計(jì)和穩(wěn)定性分析尤為重要[2?3],因此,構(gòu)建巖石在復(fù)雜應(yīng)力狀態(tài)下的非線性力學(xué)本構(gòu)模型歷來(lái)是巖土界的研究熱點(diǎn)之一。為了準(zhǔn)確反映巖石的非線性本構(gòu)關(guān)系,基于有效應(yīng)力概念和LEMAITRE等效應(yīng)變假說(shuō)[4],引入損傷變量構(gòu)建損傷本構(gòu)模型已發(fā)展成研究巖石本構(gòu)關(guān)系的一種重要方法。韋立德等[5?7]假定巖石內(nèi)部微缺陷隨機(jī)分散,巖石微元強(qiáng)度服從概率分布,結(jié)合連續(xù)損傷理論和統(tǒng)計(jì)規(guī)律,構(gòu)建了可描述巖石應(yīng)變軟化特性的彈塑性統(tǒng)計(jì)損傷本構(gòu)模型;曹瑞瑯等[8?10]基于體積變化和峰后殘余強(qiáng)度對(duì)損傷變量的影響,修正模型參數(shù)并考慮損傷閾值,使巖石統(tǒng)計(jì)損傷本構(gòu)研究得到進(jìn)一步發(fā)展。但巖石損傷同時(shí)受材料內(nèi)的微缺陷和巖石微元所受應(yīng)力狀態(tài)的影響,上述模型存在一定的局限性。巖石變形破壞的實(shí)質(zhì)是其內(nèi)部能量耗散和釋放的過(guò)程,能量耗散是損傷逐步發(fā)展的內(nèi)動(dòng)力[11],為更深層次地分析巖石破壞過(guò)程的能量耗散與損傷演化的關(guān)系,眾多研究者從能量角度對(duì)巖石非線性力學(xué)特性展開(kāi)研究并取得了許多成果。劉天為等[12?13]開(kāi)展了巖石室內(nèi)三軸壓縮試驗(yàn),通過(guò)計(jì)算加載各階段的能量轉(zhuǎn)換量,研究了巖石加載變形過(guò)程的能量耗散規(guī)律;謝和平等[14?17]在不可逆熱力學(xué)框架下,基于能量耗散原理推導(dǎo)了損傷演化方程并建立了巖石損傷本構(gòu);依據(jù)巖石損傷演化的能量耗散過(guò)程符合最小耗能原理的客觀規(guī)律,高瑋 等[18?19]運(yùn)用最小耗能原理研究了花崗巖在單軸壓縮下的損傷演變規(guī)律,為巖石損傷研究提供了一條新路徑??梢?jiàn),研究復(fù)雜應(yīng)力條件下巖石的非線性本構(gòu)關(guān)系時(shí),很少將最小耗能原理運(yùn)用于巖石損傷演化方程的建立,已有基于能量原理的損傷本構(gòu)中較少考慮損傷閾值的影響。鑒于此,本文作者基于最小耗能原理和連續(xù)損傷理論,引入巖石統(tǒng)一能量屈服準(zhǔn)則為耗能約束條件,建立基于最小耗能原理的新型損傷本構(gòu)模型。此外,對(duì)角礫熔巖開(kāi)展不同圍壓下的常規(guī)三軸壓縮試驗(yàn),結(jié)合試驗(yàn)結(jié)果分析巖石變形各階段的能量轉(zhuǎn)換關(guān)系,在此基礎(chǔ)上驗(yàn)證所建損傷本構(gòu)模型的合理性。

    1 巖石損傷本構(gòu)模型

    1.1 損傷演化方程推導(dǎo)

    巖石內(nèi)部從萌生初始裂隙到貫通為宏觀裂紋導(dǎo)致整體破壞,其實(shí)質(zhì)是損傷變量逐漸增大,直至臨界值的過(guò)程。巖石損傷演化伴隨著能量的不可逆耗散,在封閉無(wú)熱交換的系統(tǒng)內(nèi),巖石單元耗散能d滿足:

    式中:為巖石單元總應(yīng)變能,由試驗(yàn)機(jī)壓縮巖樣作功1和巖樣沿徑向變形壓縮液壓油作功2組成;e為巖石單元可釋放彈性應(yīng)變能;(=1,2,3)為巖石單元的名義主應(yīng)力和主應(yīng)變。

    在常規(guī)三軸應(yīng)力狀態(tài)下(2=3),假定加載中巖石損傷為各向同性且泊松比不受損傷影響,即1=2=3=,1=2=3=(其中,為損傷巖石的卸載彈性模量,為材料損傷前的名義泊松比)。將泊松比=?3/1代入式(1)~(3)得

    在巖石三軸壓縮試驗(yàn)中,通常先施加圍壓到預(yù)設(shè)值,再施加軸向應(yīng)力至巖樣破壞?,F(xiàn)階段計(jì)算單元應(yīng)變能時(shí)一般不考慮施加軸向應(yīng)力前巖石在圍壓作用下產(chǎn)生的初始應(yīng)變10和初始應(yīng)變能0,為精確反映巖樣在各階段的能量轉(zhuǎn)換特征,結(jié)合式(4)對(duì)巖石應(yīng)變能進(jìn)行修正,得到巖石單元實(shí)際應(yīng)變能1:

    式中:10為初始?jí)好芏蔚某跏紤?yīng)變;0為巖石損傷前的初始彈性模量。在主應(yīng)力空間中,巖石單元可釋放彈性應(yīng)變能[17]為

    式中:為巖石的卸載彈性模量。為簡(jiǎn)化數(shù)值分析,取初始彈性模量0代入計(jì)算。

    結(jié)合式(6)和式(7),整理得常規(guī)三軸應(yīng)力狀態(tài)下的巖石單元耗散能為

    在巖石損傷破壞各階段,將系統(tǒng)的熱力學(xué)變化看作是一種平衡被打破和重建的過(guò)程。當(dāng)巖石內(nèi)產(chǎn)生裂隙時(shí),當(dāng)前平衡狀態(tài)被打破;當(dāng)能量耗散趨于穩(wěn)定時(shí),系統(tǒng)處于新的平衡狀態(tài)。因此,巖石的損傷演化是符合最小耗能原理的能量耗散過(guò)程,可以通過(guò)研究能量耗散率分析損傷演化規(guī)律。

    最小耗能原理[20]指出:所有耗能過(guò)程都在對(duì)應(yīng)的約束條件下沿最小耗能路徑進(jìn)行。“沿最小耗能路徑進(jìn)行”指瞬時(shí)耗能率對(duì)應(yīng)于耗能過(guò)程的任何時(shí)刻均為最小值,物理變量需滿足的控制方程即耗能約束條件。

    在彈性損傷模型中,假定損傷產(chǎn)生的不可逆應(yīng)變?yōu)閹r石破壞過(guò)程中的唯一耗能機(jī)制,材料耗能率定義為

    將式(11)代入式(9),耗能率可表示為

    式中:和為材料參數(shù),表征巖石初始屈服時(shí)剪切應(yīng)變能與體積變形能的數(shù)值關(guān)系;為內(nèi)摩擦角;和分別為剪切模量和體積模量。

    根據(jù)最小耗能原理,在材料耗能損傷過(guò)程的任意時(shí)刻,耗能率都在對(duì)應(yīng)耗能約束條件下取駐值,引入待定Lagrange乘子*,可表達(dá)為

    結(jié)合式(13)~式(15),推導(dǎo)得損傷演化方程為

    巖石損傷存在閾值問(wèn)題[9],巖石損傷當(dāng)應(yīng)變達(dá)到閾值后開(kāi)始發(fā)展,在損傷閾值前處于線彈性變形階段。因而,需要考慮損傷演化規(guī)律的閾值影響,引入損傷閾值應(yīng)變0后的損傷演化方程表示為

    對(duì)于常規(guī)三軸壓縮試驗(yàn),試驗(yàn)實(shí)測(cè)軸向偏應(yīng)力1t為實(shí)際軸向表觀應(yīng)力1與圍壓3的差值,表示為

    為真實(shí)反映應(yīng)力?應(yīng)變關(guān)系,需要考慮巖石初始?jí)好苓^(guò)程,根據(jù)式(5)和實(shí)測(cè)應(yīng)變1t將軸向表觀應(yīng)變1修正為

    將式(18)~(20)代入式(10),得基于最小耗能原理的巖石損傷本構(gòu)為

    1.2 損傷閾值確定和模型參數(shù)辨識(shí)

    損傷閾值通常需要借助CT實(shí)時(shí)細(xì)觀試驗(yàn)或聲發(fā)射技術(shù)來(lái)進(jìn)行定量分析。EBERHARDT等[22]研究花崗巖的脆性破壞特性時(shí),定義損傷閾值為體積應(yīng)變?軸向應(yīng)變曲線駐點(diǎn)處的應(yīng)力水平,并分析了花崗巖損傷階段的力學(xué)參數(shù)演化規(guī)律。鑒于巖石在體積擴(kuò)容起始點(diǎn)之后,巖石內(nèi)部產(chǎn)生塑性變形且力學(xué)性質(zhì)開(kāi)始弱化,為簡(jiǎn)化分析,將初始屈服點(diǎn)與體積擴(kuò)容應(yīng)力點(diǎn)建立聯(lián)系,并確定體積擴(kuò)容點(diǎn)為損傷閾值。根據(jù)角礫熔巖在不同圍壓下的三軸壓縮試驗(yàn)結(jié)果,確定體積應(yīng)變?軸向應(yīng)變曲線駐點(diǎn)處的軸向應(yīng)變?yōu)閾p傷閾值應(yīng)變0。

    考慮角礫熔巖在低圍壓下的應(yīng)變軟化特征,采用巖石應(yīng)力?應(yīng)變曲線的峰值點(diǎn)強(qiáng)度辨識(shí)模型參數(shù),使參數(shù)具有明確的物理含義。假定圍壓下巖石的應(yīng)力?應(yīng)變曲線峰值點(diǎn)處應(yīng)力、軸向應(yīng)變和徑向應(yīng)變?yōu)閟,1s和3s;由應(yīng)力?應(yīng)變曲線的極值性,在峰值點(diǎn)處,有

    2 角礫熔巖三軸壓縮試驗(yàn)及能量演化特征

    2.1 角礫熔巖三軸壓縮試驗(yàn)

    角礫熔巖三軸壓縮試驗(yàn)在河海大學(xué)巖石流變力學(xué)實(shí)驗(yàn)室?guī)r石三軸伺服測(cè)試系統(tǒng)上完成,對(duì)角礫熔巖巖樣開(kāi)展圍壓為6,8,10 MPa的常規(guī)三軸壓縮試驗(yàn)。按照國(guó)際巖石力學(xué)學(xué)會(huì)(ISRM)標(biāo)準(zhǔn)[23],將角礫熔巖原狀樣打磨成直徑×高度為50 mm×100 mm的標(biāo)準(zhǔn)圓柱體,所選巖樣的均勻性和完整性均滿足試驗(yàn)規(guī)范,對(duì)應(yīng)于圍壓6,8和10 MPa下的巖樣編號(hào)分別為B-J-01,B-J-02和B-J-03。角礫熔巖試樣如圖1所示。

    試驗(yàn)以軸向位移控制方法實(shí)施,先按靜水壓力條件施加圍壓至預(yù)設(shè)值,待圍壓施加完畢并穩(wěn)定后按0.02 mm/min速率施加軸向偏壓至巖樣破壞。在整個(gè)加載過(guò)程中,伺服測(cè)試系統(tǒng)自動(dòng)記錄施加的軸向荷載,并由安裝在壓力室底座的2個(gè)高精度LVDT和固定在試樣上的1個(gè)應(yīng)變測(cè)量環(huán)分別記錄軸向和徑向變形量。

    圍壓/MPa:(a) 6;(b) 8;(c) 10

    圖2所示為破壞后的巖石試樣。從圖2可看出:巖樣在三軸壓縮應(yīng)力路徑下均為單截面剪切破壞,且在低圍壓狀態(tài)下由于圍壓對(duì)巖樣徑向約束作用較弱,破壞時(shí)巖樣兩端易于產(chǎn)生張拉性裂紋;在高圍壓下,因靜水壓力約束試樣徑向膨脹變形,抑制了巖樣內(nèi)部張性裂隙的發(fā)育和擴(kuò)展,其破壞形式主要為剪切滑移破壞。

    圍壓/MPa:(a) 6;(b) 8;(c) 10

    按照上述試驗(yàn)方案,對(duì)角礫熔巖開(kāi)展常規(guī)三軸壓縮試驗(yàn)得到應(yīng)力?應(yīng)變關(guān)系曲線和體積應(yīng)變?軸向應(yīng)變曲線,分別如圖3和圖4所示。從圖3可見(jiàn):隨著圍壓增加,抗壓強(qiáng)度分別為147.04,170.10和185.95 MPa,表明初始圍壓的增加能夠顯著提高巖石的強(qiáng)度;在初始圍壓下,巖樣內(nèi)部微裂隙逐漸閉合且顆粒骨架更密實(shí);此外,加載時(shí)圍壓作用限制巖石內(nèi)原生裂隙面間的摩擦滑移,使內(nèi)部裂紋貫通更困難,增強(qiáng)了巖石抗壓能力。根據(jù)室內(nèi)三軸壓縮試驗(yàn)結(jié)果,整理得出角礫熔巖的力學(xué)強(qiáng)度參數(shù)如表1所示。

    圍壓/MPa:(a) 6;(b) 8;(c) 10

    據(jù)圖4所示應(yīng)變關(guān)系曲線可確定體積應(yīng)變?軸向應(yīng)變曲線駐點(diǎn)為初始屈服點(diǎn),計(jì)算得到不同圍壓下的角礫熔巖在初始屈服處的剪切應(yīng)變能與體積應(yīng)變能近似符合線性關(guān)系,表明巖石統(tǒng)一能量屈服準(zhǔn)則能較好地描述角礫熔巖的屈服特性,結(jié)合試驗(yàn)結(jié)果擬合式(14),確定屈服準(zhǔn)則(1,3)的材料參數(shù)=0.032,=7.960。

    圍壓/MPa:1—6;2—8;3—10。

    表1 角礫熔巖常規(guī)三軸壓縮試驗(yàn)力學(xué)參數(shù)

    注:黏聚力為15.68 MPa;內(nèi)摩擦角為53.39°。

    2.2 三軸壓縮角礫熔巖的能量特征

    基于能量數(shù)值計(jì)算原理,采用表1所示的力學(xué)參數(shù),計(jì)算得到角礫熔巖壓縮破壞過(guò)程中能量演化與偏應(yīng)力的關(guān)系如圖5所示。根據(jù)能量演化規(guī)律和應(yīng)力?應(yīng)變關(guān)系,巖石破壞過(guò)程可分為4個(gè)階段。

    1) 初始非線性段()。巖石單元可釋放彈性能U增大速度緩慢,表明加載初期外荷壓縮巖樣做功主要用于材料內(nèi)部孔隙閉合和顆粒骨架彈性變形。

    2) 線彈性段()。隨著偏應(yīng)力逐漸增大,巖樣內(nèi)部原生裂隙閉合完全,應(yīng)力?應(yīng)變關(guān)系趨于線性且?guī)r石內(nèi)部顆粒骨架產(chǎn)生線彈性變形并將外荷做功儲(chǔ)存為可釋放彈性勢(shì)能,巖石吸收總能量和可釋放彈性能快速增加,耗散能仍不增大。

    圍壓/MPa:(a) 6;(b) 8;(c) 10

    3) 破裂發(fā)展段()。當(dāng)應(yīng)力水平超過(guò)初始屈服點(diǎn)后,巖石內(nèi)部微裂紋快速擴(kuò)展并產(chǎn)生塑性變形,損傷開(kāi)始發(fā)展導(dǎo)致巖石強(qiáng)度逐漸弱化,巖石呈現(xiàn)出非線性力學(xué)特征??舍尫艔椥阅茉鏊僮兙彶⒃诜逯祽?yīng)力處達(dá)到極值而耗散能快速增加,巖石單元總應(yīng)變能以塑性能、熱能等其他形式被耗散。

    4) 峰后應(yīng)變軟化段()。隨著巖石損傷逐漸加強(qiáng),內(nèi)部裂紋不斷擴(kuò)展并貫通為宏觀破壞面,能量達(dá)到巖石破壞所需的耗散能時(shí)快速釋放導(dǎo)致巖樣整體破壞。

    表2所示為不同圍壓下的角礫熔巖三軸壓縮過(guò)程的能量特征統(tǒng)計(jì)值,對(duì)應(yīng)于峰值強(qiáng)度的巖石單元可釋放彈性能和吸收總應(yīng)變能分別定義為極限彈性能和破壞總應(yīng)變能。分析表2可知:峰值處耗散能和極限彈性能表現(xiàn)出明顯的圍壓效應(yīng),且初始屈服應(yīng)力也隨著圍壓的增大而增高,B-J-01巖樣在初始圍壓6 MPa下破壞時(shí)峰值處耗散能為0.092 MJ·cm?3,極限彈性能為0.393 MJ·cm?3;B-J-02巖樣施加初始圍壓8 MPa,破壞時(shí)峰值處耗散能為0.237 MJ·cm?3,極限彈性能為0.434 MJ·cm?3;B-J-03巖樣在初始圍壓10 MPa下,巖樣破壞時(shí)峰值處耗散能為0.303 MJ·cm?3,極限彈性能為0.521 MJ·cm?3;對(duì)應(yīng)于6,8,10 MPa初始圍壓下,巖樣在屈服階段損傷發(fā)展耗散的能量占破壞總應(yīng)變能的比例分別為30.75%,48.86%,56.39%。

    結(jié)合圖5分析角礫熔巖的能量演化規(guī)律,能量主要耗散于巖石的后屈服變形段,巖石內(nèi)部損傷隨著耗散能的快速增加而開(kāi)始發(fā)展,巖石剛度逐漸弱化;當(dāng)耗散的能量達(dá)到巖石破壞所需的耗散能時(shí),其內(nèi)部可釋放彈性能導(dǎo)致巖樣整體破壞,且初始圍壓愈高,巖石破壞時(shí)耗散的能量愈多,在損傷發(fā)展階段耗散的能量占總應(yīng)變能的比例愈大;在高圍壓條件下,巖石破壞時(shí)內(nèi)部積累的可釋放彈性能越多,發(fā)生破壞的程度越劇烈。因此,能量耗散是巖石內(nèi)部裂隙發(fā)展的內(nèi)動(dòng)力,能夠在分析能量演化規(guī)律的基礎(chǔ)上,通過(guò)研究能量耗散率以分析損傷演化規(guī)律。

    表2 角礫熔巖三軸壓縮的能量特征值

    3 模型驗(yàn)證及損傷演化規(guī)律分析

    表3 模型參數(shù)

    將表3中的模型參數(shù)代入式(18),得出不同圍壓下的巖石損傷演化曲線如圖7所示。從圖7可見(jiàn):在加載初期,隨著應(yīng)力水平增加,巖石內(nèi)部的微裂隙和微空洞逐漸閉合且顆粒骨架產(chǎn)生彈性變形,但產(chǎn)生的軸向應(yīng)變1較小且低于損傷閾值應(yīng)變0,巖石應(yīng)力?應(yīng)變關(guān)系基本呈線性關(guān)系,損傷變量=0,即巖樣在線彈性階段不產(chǎn)生損傷;隨著應(yīng)力水平進(jìn)一步提高,當(dāng)巖樣所受應(yīng)力超過(guò)一定水平時(shí),軸向應(yīng)變1高于損傷閾值應(yīng)變后損傷開(kāi)始發(fā)展,且損傷變量與軸向應(yīng)變呈非線性關(guān)系;巖石內(nèi)部的微裂隙快速擴(kuò)展,局部裂隙密集并貫通產(chǎn)生宏觀裂紋,裂紋間經(jīng)階梯狀連接成主破壞面。根據(jù)角礫熔巖的能量演化規(guī)律,巖石內(nèi)部的裂隙發(fā)展需要耗散能量,且在屈服階段耗散能快速增大,當(dāng)應(yīng)力水平趨近峰值強(qiáng)度時(shí)損傷變量急劇增大,致使巖石剛度顯著弱化,表明損傷演化規(guī)律與能量耗散規(guī)律和試驗(yàn)曲線的非線性特征相符。對(duì)比不同圍壓下的損傷演化曲線發(fā)現(xiàn):隨著初始圍壓增加,對(duì)應(yīng)于不同圍壓下巖石的損傷閾值逐漸增大,表明圍壓作用改善了巖石內(nèi)部受力狀態(tài),從而抑制了微裂隙 發(fā)展。

    圍壓/MPa:(a) 6;(b) 8;(c) 10

    圍壓/MPa:1—6;2—8;3—10

    4 結(jié)論

    1) 基于連續(xù)損傷理論和各向同性假定,運(yùn)用最小耗能原理導(dǎo)出巖石損傷演化方程,建立了復(fù)雜應(yīng)力條件下的巖石損傷本構(gòu)模型。對(duì)比分析角礫熔巖常規(guī)三軸壓縮試驗(yàn)結(jié)果與理論結(jié)果,驗(yàn)證了損傷模型的合理性,表明模型能較好地反映巖石在不同圍壓下的應(yīng)力應(yīng)變關(guān)系。

    2) 巖石損傷演化的實(shí)質(zhì)是能量耗散。隨著圍壓增加,巖石在損傷發(fā)展階段耗散的能量占總應(yīng)變能的比例越大,且其破壞時(shí)需耗散的能量越多;破壞時(shí),巖石內(nèi)積累的可釋放極限彈性能具有明顯的圍壓效應(yīng),初始圍壓越高,破壞釋放的彈性能越多,破壞程度越劇烈。

    3) 結(jié)合巖石損傷特性考慮損傷演化規(guī)律的閾值問(wèn)題,將巖石的初始屈服點(diǎn)與體積擴(kuò)容應(yīng)力點(diǎn)建立聯(lián)系,以巖石體積擴(kuò)容點(diǎn)作為損傷閾值,本模型能合理地反映巖石變形前期的線彈性特性和后屈服段的非線性力學(xué)行為,更加接近實(shí)際情況。

    [1] 徐志英. 巖石力學(xué)[M]. 3版. 北京: 水利水電出版社, 2007: 69?78.XU Zhiyin. Rock mechanics[M]. 3rd ed. Beijing: Water Resources and Electric Power Press, 2007: 69?78.

    [2] 李建林, 王樂(lè)華. 節(jié)理巖體卸荷非線性力學(xué)特性研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2007, 26(10): 1968?1975.LI Jianlin, WANG Lehua. Study on unloading nonlinear mechanical characteristics of jointed rock mass[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(10): 1968?1975.

    [3] 徐衛(wèi)亞, 宋曉晨, 周維垣. 水電站進(jìn)水口巖石高邊坡及壩坡與洞室相互作用的三維數(shù)值分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 2004, 23(16): 2712?2717.XU Weiya, SONG Xiaochen, ZHOU Weiyuan. Numerical analysis on deformation and stress of high slope of intake and interaction of dam-slope and tunnel at huge hydropower project[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(16): 2712?2717.

    [4] LEMAITRE J. How to use damage mechanics[J]. Nuclear Engineering & Design, 1984, 80(2): 233?245.

    [5] 韋立德, 徐衛(wèi)亞, 楊春和, 等. 具有統(tǒng)計(jì)損傷的巖石彈塑性本構(gòu)模型研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2004, 23(12): 1971?1975. WEI Lide, XU Weiya, YANG Chunhe, et al. Study on elastoplastic constitutive model of rock with statistical damage[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(12): 1971?1975.

    [6] 劉齊建, 楊林德, 曹文貴. 巖石統(tǒng)計(jì)損傷本構(gòu)模型及其參數(shù)反演[J]. 巖石力學(xué)與工程學(xué)報(bào), 2005, 24(4): 616?621.LIU Qijian, YANG Linde, CAO Wengui. Statistical damage constitutive model for rock and back analysis of its parameters[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(4): 616?621.

    [7] 楊圣奇, 徐衛(wèi)亞, 蘇承東. 考慮尺寸效應(yīng)的巖石損傷統(tǒng)計(jì)本構(gòu)模型研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2005, 24(24): 4484?4490. YANG Shengqi, XU Weiya, SU Chengdong. Study on statistical damage constitutive model of rock considering scale effect[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(24): 4484?4490.

    [8] 曹瑞瑯, 賀少輝, 韋京, 等. 基于殘余強(qiáng)度修正的巖石損傷軟化統(tǒng)計(jì)本構(gòu)模型研究[J]. 巖土力學(xué), 2013, 34(6): 1652?1667. CAO Ruilang, HE Shaohui, WEI Jing, et al. Study of modified statistical damage softening constitutive model for rock considering residual strength[J]. Rock and Soil Mechanics, 2013, 34(6): 1652?1667.

    [9] 曹文貴, 趙衡, 張玲, 等. 考慮損傷閥值影響的巖石損傷統(tǒng)計(jì)軟化本構(gòu)模型及其參數(shù)確定方法[J]. 巖石力學(xué)與工程學(xué)報(bào), 2008, 27(6): 1148?1154. CAO Wengui, ZHAO Heng, ZHANG Ling, et al. Damage statistical softening constitutive model for rock considering effect of damage threshold and its parameters determination method[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(6): 1148?1154.

    [10] 王偉, 田振元, 朱其志, 等. 考慮孔隙水壓力的巖石統(tǒng)計(jì)損傷本構(gòu)模型研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2015, 34(S2): 3676?3682.WANG Wei, TIAN Zhenyuan, ZHU Qizhi, et al. Study of statistical damage constitutive model for rock considering pore water pressure[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(S2): 3676?3682.

    [11] 謝和平, 鞠楊, 黎立云, 等. 巖體變形破壞過(guò)程的能量機(jī)制[J]. 巖石力學(xué)與工程學(xué)報(bào), 2008, 27(9): 1729?1740. XIE Heping, JU Yang, LI Liyun, et al. Energy mechanism of deformation and failure of rock masses[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(9): 1729?1740.

    [12] 劉天為, 何江達(dá), 徐文杰. 大理巖三軸壓縮破壞的能量特征分析[J]. 巖土工程學(xué)報(bào), 2013, 35(2): 395?400. LIU Tianwei, HE Jiangda, XU Wenjie. Energy properties of failure of marble samples under triaxial compression[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(2): 395?400.

    [13] 田勇, 俞然剛. 不同圍壓下灰?guī)r三軸壓縮過(guò)程能量分析[J]. 巖土力學(xué), 2014, 35(1): 118?122. TIAN Yong, YU Rangang. Energy analysis of limestone during triaxial compression under different confining pressures[J]. Rock and Soil Mechanics, 2014, 35(1): 118?122.

    [14] 謝和平, 鞠楊, 黎立云. 基于能量耗散與釋放原理的巖石強(qiáng)度與整體破壞準(zhǔn)則[J]. 巖石力學(xué)與工程學(xué)報(bào), 2005, 24(17): 3003?3010. XIE Heping, JU Yang, LI Liyun. Criteria for strength and structural failure of rocks based on energy dissipation and energy release principles[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(17): 3003?3010.

    [15] LIU Xuesheng, NING J G, TAN Yunliang, et al. Damage constitutive model based on energy dissipation for intact rock subjected to cyclic loading[J]. International Journal of Rock Mechanics & Mining Sciences, 2016, 85: 27?32.

    [16] 郭建強(qiáng), 劉新榮, 王景環(huán), 等. 基于能量原理鹽巖的損傷本構(gòu)模型[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 44(12): 5045?5050. GUO Jianqiang, LIU Xinrong, WANG Jinghuan, et al. Damage constitutive model of rock salt based on energy principles[J]. Journal of Central South University (Science and Technology), 2013, 44(12): 5045?5050.

    [17] 劉曉明, 熊力, 劉建華, 等. 基于能量耗散原理的紅砂巖崩解機(jī)制研究[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 42(10): 3143?3149. LIU Xiaoming, XIONG Li, LIU Jianhua, et al. Slacking mechanism of red sandstone based on energy dissipation principle[J]. Journal of Central South University(Science and Technology), 2011, 42(10): 3143?3149.

    [18] 高瑋, 汪磊, 楊大勇. 巖石損傷演化的能量方法研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2011, 30(S2): 4087?4092.GAO Wei, WANG Lei, YANG Dayong. Rock damage evolution based on energy principle[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(S2): 4087?4092.

    [19] GAO Wei, WANG Xu, DAI Shuang, et al. Numerical study on stability of rock slope based on energy method[J]. Advances in Materials Science and Engineering, 2016, 2016: 1?10.

    [20] 周筑寶. 最小耗能原理及其應(yīng)用[M]. 北京: 科學(xué)出版社, 2001: 7?10. ZHOU Zhubao. The least energy consumption principle and its applications[M]. Beijing: Science Press, 2001: 7?10.

    [21] 周輝, 李震, 朱國(guó)金, 等. 基于巖石統(tǒng)一能量屈服準(zhǔn)則的硬巖損傷模型[J]. 巖土力學(xué), 2016, 37(3): 609?615. ZHOU Hui, LI Zhen, ZHU Guojin, et al. A damage model for hard rock based on unified energy yield criterion of rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 37(3): 609?615.

    [22] EBERHARDT E, STEAD D, STIMPSON B, et al. Identifying crack initiation and propagation thresholds in brittle roc[J]. Canadian Geotechnical Journal, 1998, 35(2): 222?233.

    [23] ULUSAY R. The ISRM suggested methods for rock characterization, testing and monitoring: 2007—2014[J]. Springer International Publishing, 2014, 15(1): 47?48.

    (編輯 陳燦華)

    Study on damage constitutive model of rock based on principle of minimum dissipative energy

    SUN Mengcheng1, 2, XU Weiya1, 2, WANG Susheng1, 2, WANG Rubin1, 2, WANG Wei1, 2

    (1. Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering, Hohai University, Nanjing 210098, China; 2. Institute of Geotechnical Engineering, Hohai University, Nanjing 210098, China)

    The principle of minimum dissipative energy was introduced based on the continuum damage theory. Taking unified energy yield criterion of rock as the constraints of energy dissipation, a new damage constitutive model was proposed based on the assumption of isotropy. The new damage constitutive model took into account the influence of damage threshold which was assumed to be the starting point of initialyield phase. Furthermore, the theoretical correlation between model parameter and deformation moduli of rock under triaxial compression was deduced. A series of conventional triaxial tests on breccia lava samples were carried out under different confining pressures, and evolution characteristics of energy at different stages of rock deformation and failure were analyzed to verify the rationality of model. The results show that the model can not only reflect the nonlinear mechanical characteristics of rock materials under complicated stress condition, but also has high fitting precision on post-yield stage of rock, which can provide reference for safety analysis of rock engineering under complicated stress condition.

    rock mechanics; the principle of minimum dissipative energy; unified energy yield criterion of rock; damage threshold; damage evolution equation

    TU452

    A

    1672?7207(2018)08?2067?09

    10.11817/j.issn.1672?7207.2018.08.029

    2017?09?08;

    2017?11?12

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51479049,51679069);江蘇省研究生科研與實(shí)踐創(chuàng)新計(jì)劃項(xiàng)目(KYCX17_0471, 2017B700X14)(Projects(51479049, 51679069) supported by the National Natural Science Foundation of China; Project(KYCX17_0471, 2017B700X14) supported by the Postgraduate Research & Practice Innovation Program of Jiangsu Province)

    徐衛(wèi)亞,博士,教授,從事巖土力學(xué)與工程研究;E-mail:wyxu@hhu.edu.cn

    猜你喜歡
    角礫熔巖巖樣
    熔巖管道
    裂隙幾何特征對(duì)巖體力學(xué)特性的模擬試驗(yàn)研究
    煤礦安全(2023年7期)2023-08-04 03:39:06
    鄂爾多斯盆地馬家溝組準(zhǔn)層狀巖溶角礫巖特征及其成因*
    預(yù)制裂紋砂巖樣力學(xué)特性與能量演化規(guī)律研究
    內(nèi)蒙古林東普通球粒隕石的角礫特征
    天津薊縣晚元古代冰磧巖的發(fā)現(xiàn)
    頻繁動(dòng)力擾動(dòng)對(duì)圍壓卸載中高儲(chǔ)能巖體的動(dòng)力學(xué)影響*
    爆炸與沖擊(2019年6期)2019-07-10 01:37:50
    泡泡熔巖燈
    離開(kāi)熔巖島
    熔巖燈
    美女内射精品一级片tv| 日本爱情动作片www.在线观看| 亚洲av男天堂| 人人妻人人澡欧美一区二区| 99视频精品全部免费 在线| 一级毛片aaaaaa免费看小| 国产乱来视频区| 亚洲国产精品成人久久小说| 国产亚洲av嫩草精品影院| 最近中文字幕高清免费大全6| av在线天堂中文字幕| 久久精品国产自在天天线| 91久久精品电影网| 国产有黄有色有爽视频| 免费av不卡在线播放| 色网站视频免费| 亚洲成色77777| av福利片在线观看| 熟妇人妻久久中文字幕3abv| 免费av不卡在线播放| 你懂的网址亚洲精品在线观看| 国产高潮美女av| 街头女战士在线观看网站| 人人妻人人澡欧美一区二区| 亚洲电影在线观看av| 天堂av国产一区二区熟女人妻| 女的被弄到高潮叫床怎么办| 女人被狂操c到高潮| 成人鲁丝片一二三区免费| 亚洲av电影不卡..在线观看| 九九爱精品视频在线观看| 极品教师在线视频| 精品久久久久久成人av| 成人亚洲精品av一区二区| 欧美极品一区二区三区四区| 成人亚洲精品av一区二区| 内地一区二区视频在线| 久久久久久久久久久丰满| 国产人妻一区二区三区在| 人妻少妇偷人精品九色| 亚洲精品第二区| 国产精品蜜桃在线观看| 极品教师在线视频| 最近最新中文字幕大全电影3| 乱码一卡2卡4卡精品| 欧美极品一区二区三区四区| 国产美女午夜福利| 久久韩国三级中文字幕| 91精品伊人久久大香线蕉| 久久午夜福利片| 好男人在线观看高清免费视频| 晚上一个人看的免费电影| 特级一级黄色大片| 国产伦精品一区二区三区视频9| 女人十人毛片免费观看3o分钟| 又爽又黄无遮挡网站| 国产午夜精品论理片| 亚洲精品,欧美精品| 青春草亚洲视频在线观看| 午夜激情欧美在线| 亚洲无线观看免费| 蜜臀久久99精品久久宅男| 久久久精品欧美日韩精品| 卡戴珊不雅视频在线播放| 欧美成人一区二区免费高清观看| av卡一久久| 国产精品久久视频播放| 一二三四中文在线观看免费高清| 欧美xxxx黑人xx丫x性爽| 精品久久国产蜜桃| 国产不卡一卡二| 午夜免费男女啪啪视频观看| 熟妇人妻久久中文字幕3abv| 国产午夜精品一二区理论片| 中文字幕免费在线视频6| 免费看av在线观看网站| 麻豆av噜噜一区二区三区| 青春草国产在线视频| 亚洲婷婷狠狠爱综合网| 麻豆国产97在线/欧美| 噜噜噜噜噜久久久久久91| 亚洲精品成人久久久久久| 日韩一区二区视频免费看| 色视频www国产| 2021少妇久久久久久久久久久| 大香蕉97超碰在线| 老司机影院成人| 久99久视频精品免费| 成年免费大片在线观看| 女的被弄到高潮叫床怎么办| 午夜爱爱视频在线播放| 如何舔出高潮| 午夜免费男女啪啪视频观看| 精品国产露脸久久av麻豆 | 亚洲国产精品国产精品| 国产高清三级在线| 国产色爽女视频免费观看| 亚洲人成网站高清观看| 直男gayav资源| 美女被艹到高潮喷水动态| 久久久久久久久久久免费av| 又爽又黄a免费视频| 特大巨黑吊av在线直播| 国产在线一区二区三区精| 国产激情偷乱视频一区二区| h日本视频在线播放| 久久99精品国语久久久| 五月玫瑰六月丁香| 少妇猛男粗大的猛烈进出视频 | 在线免费观看不下载黄p国产| 亚洲精品中文字幕在线视频 | 免费人成在线观看视频色| 久久久久久久久久人人人人人人| 亚洲熟女精品中文字幕| or卡值多少钱| 久久久久精品性色| 乱码一卡2卡4卡精品| 国产成人freesex在线| 精品酒店卫生间| 亚洲人成网站在线播| 国产成人aa在线观看| 在线观看人妻少妇| 人体艺术视频欧美日本| 精品久久久久久久久亚洲| 国产一区二区亚洲精品在线观看| 亚洲高清免费不卡视频| av线在线观看网站| 夜夜爽夜夜爽视频| 日本免费在线观看一区| 男女下面进入的视频免费午夜| 亚洲精品久久午夜乱码| 熟女人妻精品中文字幕| 99热全是精品| 蜜桃久久精品国产亚洲av| 免费看光身美女| 国产免费福利视频在线观看| 日韩制服骚丝袜av| 国产高清国产精品国产三级 | 特级一级黄色大片| 色5月婷婷丁香| 日韩av在线免费看完整版不卡| 精品久久久久久电影网| 校园人妻丝袜中文字幕| 一本一本综合久久| 成人鲁丝片一二三区免费| 天堂网av新在线| 亚洲av二区三区四区| 亚洲真实伦在线观看| 国产成年人精品一区二区| 亚洲国产av新网站| 国产精品久久久久久精品电影| 亚洲aⅴ乱码一区二区在线播放| 国产精品嫩草影院av在线观看| 国产精品国产三级国产专区5o| 丝瓜视频免费看黄片| 日韩在线高清观看一区二区三区| 在线天堂最新版资源| 自拍偷自拍亚洲精品老妇| 亚洲一级一片aⅴ在线观看| 亚洲av一区综合| 久久午夜福利片| 亚洲成人一二三区av| 免费观看精品视频网站| 白带黄色成豆腐渣| 少妇被粗大猛烈的视频| 成人漫画全彩无遮挡| 国产美女午夜福利| 中文字幕久久专区| 国产成人freesex在线| 国产高潮美女av| 国产精品av视频在线免费观看| 有码 亚洲区| 97热精品久久久久久| 亚洲精品第二区| 国产精品麻豆人妻色哟哟久久 | av在线观看视频网站免费| 免费看a级黄色片| 久久这里只有精品中国| 国产国拍精品亚洲av在线观看| 亚洲av一区综合| 男人和女人高潮做爰伦理| 极品教师在线视频| 可以在线观看毛片的网站| 国产精品.久久久| 亚洲欧美一区二区三区黑人 | 欧美区成人在线视频| 久久国产乱子免费精品| 欧美精品国产亚洲| 午夜老司机福利剧场| 建设人人有责人人尽责人人享有的 | 18禁在线播放成人免费| 最后的刺客免费高清国语| 在线 av 中文字幕| 亚洲不卡免费看| 尾随美女入室| 三级男女做爰猛烈吃奶摸视频| ponron亚洲| 免费播放大片免费观看视频在线观看| 男女边吃奶边做爰视频| 亚洲va在线va天堂va国产| 超碰97精品在线观看| 日本三级黄在线观看| 亚洲欧美成人精品一区二区| 肉色欧美久久久久久久蜜桃 | 纵有疾风起免费观看全集完整版 | 亚洲精品国产成人久久av| 午夜福利在线观看吧| 蜜桃久久精品国产亚洲av| 久久精品国产亚洲av涩爱| 日韩欧美 国产精品| 国产精品一及| 校园人妻丝袜中文字幕| 水蜜桃什么品种好| 乱码一卡2卡4卡精品| 久久久久国产网址| 日产精品乱码卡一卡2卡三| 国产亚洲午夜精品一区二区久久 | 成人av在线播放网站| 亚洲一级一片aⅴ在线观看| 69人妻影院| 国产乱人偷精品视频| 伦理电影大哥的女人| 99热网站在线观看| 在线免费十八禁| 午夜精品一区二区三区免费看| 全区人妻精品视频| a级毛片免费高清观看在线播放| 日本一二三区视频观看| av女优亚洲男人天堂| 大又大粗又爽又黄少妇毛片口| 亚洲精品日韩在线中文字幕| 亚洲熟妇中文字幕五十中出| 日本一二三区视频观看| 国产探花在线观看一区二区| 亚洲最大成人中文| 日韩 亚洲 欧美在线| 91精品一卡2卡3卡4卡| 男女边摸边吃奶| 一级片'在线观看视频| 亚洲综合精品二区| videos熟女内射| 国产伦在线观看视频一区| 男的添女的下面高潮视频| 精品午夜福利在线看| 能在线免费看毛片的网站| 高清毛片免费看| 18+在线观看网站| 日韩一本色道免费dvd| 国产国拍精品亚洲av在线观看| 亚洲精品日本国产第一区| 国产成人a区在线观看| 国产熟女欧美一区二区| 亚洲怡红院男人天堂| 精品久久久久久久末码| 男人狂女人下面高潮的视频| 亚洲性久久影院| 人人妻人人澡欧美一区二区| 少妇被粗大猛烈的视频| 国产一区二区亚洲精品在线观看| 亚洲高清免费不卡视频| 午夜激情福利司机影院| 人人妻人人看人人澡| 国内少妇人妻偷人精品xxx网站| 亚洲av成人av| 五月天丁香电影| 国产成人精品一,二区| 国产精品一区二区在线观看99 | 亚洲国产精品专区欧美| 一二三四中文在线观看免费高清| 亚洲久久久久久中文字幕| 大香蕉久久网| 少妇丰满av| av卡一久久| 91精品国产九色| 一个人看视频在线观看www免费| 中文字幕av成人在线电影| 久久久国产一区二区| 成人性生交大片免费视频hd| 亚洲三级黄色毛片| 国内少妇人妻偷人精品xxx网站| 美女黄网站色视频| 大片免费播放器 马上看| 女的被弄到高潮叫床怎么办| 国产伦在线观看视频一区| 最近的中文字幕免费完整| 成人一区二区视频在线观看| 欧美日韩综合久久久久久| 国产黄频视频在线观看| 国产成人福利小说| 婷婷色综合www| 又黄又爽又刺激的免费视频.| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影| 免费看美女性在线毛片视频| 非洲黑人性xxxx精品又粗又长| 99热这里只有是精品50| 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩av片在线观看| 美女xxoo啪啪120秒动态图| 亚洲精品第二区| 国产精品人妻久久久久久| 国产午夜精品久久久久久一区二区三区| 激情 狠狠 欧美| 少妇人妻精品综合一区二区| 亚洲在线观看片| 高清午夜精品一区二区三区| 精品一区二区免费观看| 精品久久久久久久末码| 日本三级黄在线观看| 伦理电影大哥的女人| 又大又黄又爽视频免费| 午夜福利视频精品| 亚洲精华国产精华液的使用体验| 精品久久久久久电影网| 2021少妇久久久久久久久久久| 精品国内亚洲2022精品成人| 中文欧美无线码| 一级毛片我不卡| 亚洲,欧美,日韩| 精华霜和精华液先用哪个| 五月天丁香电影| 99久久九九国产精品国产免费| a级毛片免费高清观看在线播放| 内地一区二区视频在线| 男女啪啪激烈高潮av片| 在线观看免费高清a一片| 久久久久久久久中文| 免费播放大片免费观看视频在线观看| 欧美xxxx性猛交bbbb| 不卡视频在线观看欧美| 日日撸夜夜添| 男女视频在线观看网站免费| 欧美成人午夜免费资源| 汤姆久久久久久久影院中文字幕 | 一区二区三区免费毛片| 国产淫语在线视频| 91午夜精品亚洲一区二区三区| 日韩av免费高清视频| 乱系列少妇在线播放| 人人妻人人看人人澡| 免费看a级黄色片| 18禁在线播放成人免费| 美女高潮的动态| 国产精品人妻久久久久久| 日产精品乱码卡一卡2卡三| 嫩草影院精品99| 国产精品一二三区在线看| 国产精品不卡视频一区二区| 亚洲av中文av极速乱| 搡老乐熟女国产| 建设人人有责人人尽责人人享有的 | 午夜精品国产一区二区电影 | 麻豆精品久久久久久蜜桃| 人人妻人人看人人澡| 2022亚洲国产成人精品| 激情五月婷婷亚洲| 一级黄片播放器| 亚洲精华国产精华液的使用体验| 国产一区二区三区av在线| 美女cb高潮喷水在线观看| 国产午夜精品论理片| 免费不卡的大黄色大毛片视频在线观看 | 日韩在线高清观看一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本免费a在线| 天堂av国产一区二区熟女人妻| 少妇被粗大猛烈的视频| 少妇人妻精品综合一区二区| 肉色欧美久久久久久久蜜桃 | 国产av国产精品国产| 精品久久久久久成人av| 夫妻性生交免费视频一级片| 爱豆传媒免费全集在线观看| 久久人人爽人人片av| 黄色配什么色好看| 成人漫画全彩无遮挡| 国产综合懂色| 黄片wwwwww| 建设人人有责人人尽责人人享有的 | 黄色配什么色好看| a级毛片免费高清观看在线播放| 欧美成人午夜免费资源| 你懂的网址亚洲精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 好男人在线观看高清免费视频| av在线亚洲专区| 黄片wwwwww| 男人和女人高潮做爰伦理| 国产欧美日韩精品一区二区| 国产精品综合久久久久久久免费| 搡老妇女老女人老熟妇| 国产一区有黄有色的免费视频 | 久久精品久久久久久噜噜老黄| 亚洲国产最新在线播放| 成人美女网站在线观看视频| 精品人妻偷拍中文字幕| 亚洲精品久久午夜乱码| 国产乱人偷精品视频| 你懂的网址亚洲精品在线观看| 亚洲自偷自拍三级| 简卡轻食公司| 国产亚洲精品av在线| 婷婷六月久久综合丁香| 国产亚洲av嫩草精品影院| 小蜜桃在线观看免费完整版高清| 少妇人妻一区二区三区视频| 少妇猛男粗大的猛烈进出视频 | 午夜免费男女啪啪视频观看| 亚洲精品日韩av片在线观看| 亚洲婷婷狠狠爱综合网| 大又大粗又爽又黄少妇毛片口| 久久99蜜桃精品久久| 一本久久精品| 久久这里有精品视频免费| 国产成人91sexporn| 超碰97精品在线观看| 亚洲av中文av极速乱| 国产精品嫩草影院av在线观看| 色网站视频免费| 一二三四中文在线观看免费高清| 看黄色毛片网站| 国产精品一区二区三区四区久久| 成人二区视频| 精品熟女少妇av免费看| 国产成人福利小说| 国产男女超爽视频在线观看| 婷婷色综合大香蕉| 午夜福利视频1000在线观看| 欧美变态另类bdsm刘玥| 久久精品国产自在天天线| 亚洲欧美清纯卡通| 国产午夜精品论理片| 少妇被粗大猛烈的视频| 最近手机中文字幕大全| 国产精品人妻久久久久久| 精华霜和精华液先用哪个| 最近手机中文字幕大全| 久久久精品免费免费高清| 国产精品无大码| 一本久久精品| 免费av观看视频| 久久国产乱子免费精品| 插阴视频在线观看视频| a级一级毛片免费在线观看| 只有这里有精品99| 午夜福利视频精品| 啦啦啦韩国在线观看视频| 亚洲自拍偷在线| 春色校园在线视频观看| 国产成人一区二区在线| 美女国产视频在线观看| 狂野欧美白嫩少妇大欣赏| 天天一区二区日本电影三级| 搡女人真爽免费视频火全软件| 简卡轻食公司| 国产白丝娇喘喷水9色精品| 高清av免费在线| 99热全是精品| 插逼视频在线观看| 国产精品爽爽va在线观看网站| 高清在线视频一区二区三区| 在线观看av片永久免费下载| 丝袜美腿在线中文| 啦啦啦韩国在线观看视频| eeuss影院久久| 男人爽女人下面视频在线观看| 日韩国内少妇激情av| 美女黄网站色视频| 国产69精品久久久久777片| 在线天堂最新版资源| 日韩强制内射视频| 丝瓜视频免费看黄片| 亚洲精品久久久久久婷婷小说| 久久久色成人| 精品熟女少妇av免费看| 国产精品久久视频播放| 亚洲国产日韩欧美精品在线观看| 精品人妻视频免费看| 久久久色成人| 黄色配什么色好看| 国产黄a三级三级三级人| 欧美日韩亚洲高清精品| 久久久午夜欧美精品| 男插女下体视频免费在线播放| 人妻一区二区av| 国产亚洲91精品色在线| 国产中年淑女户外野战色| 插逼视频在线观看| 欧美性感艳星| 99热这里只有是精品50| 黄片wwwwww| 少妇人妻精品综合一区二区| 色综合色国产| 狂野欧美激情性xxxx在线观看| 国产 一区 欧美 日韩| 只有这里有精品99| 久久久色成人| 亚洲经典国产精华液单| 中国国产av一级| 99热这里只有是精品50| 3wmmmm亚洲av在线观看| 国产麻豆成人av免费视频| 国产一级毛片七仙女欲春2| 青青草视频在线视频观看| 成人综合一区亚洲| 日韩制服骚丝袜av| 成人二区视频| 亚洲精品乱码久久久久久按摩| 午夜免费激情av| 免费大片黄手机在线观看| 综合色av麻豆| 高清在线视频一区二区三区| 日本色播在线视频| 成人鲁丝片一二三区免费| 日本免费在线观看一区| 亚洲国产最新在线播放| 不卡视频在线观看欧美| 成人特级av手机在线观看| 天堂影院成人在线观看| 日本午夜av视频| 亚洲精品成人av观看孕妇| 夜夜看夜夜爽夜夜摸| 国产免费视频播放在线视频 | 久久精品久久久久久噜噜老黄| 亚洲精品中文字幕在线视频 | 免费少妇av软件| 日本wwww免费看| 欧美成人精品欧美一级黄| 神马国产精品三级电影在线观看| 亚洲欧美清纯卡通| 亚洲经典国产精华液单| 最近中文字幕高清免费大全6| 成人一区二区视频在线观看| 色吧在线观看| 丝瓜视频免费看黄片| 五月伊人婷婷丁香| 亚洲欧美中文字幕日韩二区| 国产精品一区二区在线观看99 | 久久久色成人| 肉色欧美久久久久久久蜜桃 | 一个人免费在线观看电影| 成年女人在线观看亚洲视频 | 精华霜和精华液先用哪个| 婷婷六月久久综合丁香| 成人二区视频| 亚洲真实伦在线观看| 三级毛片av免费| 人人妻人人澡人人爽人人夜夜 | 国内精品宾馆在线| 精品不卡国产一区二区三区| 少妇猛男粗大的猛烈进出视频 | 亚洲在久久综合| 激情 狠狠 欧美| 欧美变态另类bdsm刘玥| 久久久久久国产a免费观看| 夜夜爽夜夜爽视频| 爱豆传媒免费全集在线观看| 乱码一卡2卡4卡精品| 亚洲国产av新网站| 熟妇人妻不卡中文字幕| eeuss影院久久| 晚上一个人看的免费电影| 亚洲欧美日韩无卡精品| 色视频www国产| 国产不卡一卡二| 尤物成人国产欧美一区二区三区| 日韩伦理黄色片| 干丝袜人妻中文字幕| 精品久久久久久久久亚洲| 老师上课跳d突然被开到最大视频| 国产一区二区亚洲精品在线观看| 中国美白少妇内射xxxbb| 久久99热这里只频精品6学生| 人妻制服诱惑在线中文字幕| 两个人视频免费观看高清| 大香蕉久久网| 91精品伊人久久大香线蕉| 波野结衣二区三区在线| 午夜福利在线观看吧| 三级男女做爰猛烈吃奶摸视频| 性色avwww在线观看| 日本免费在线观看一区| 91久久精品国产一区二区成人| 欧美人与善性xxx| 91在线精品国自产拍蜜月| 搞女人的毛片| 国产男人的电影天堂91| 一级av片app| 美女主播在线视频| 成人av在线播放网站| 看十八女毛片水多多多| 女人十人毛片免费观看3o分钟| 天天一区二区日本电影三级| 两个人视频免费观看高清| 成人高潮视频无遮挡免费网站| 一级毛片电影观看| 亚洲最大成人中文| 日韩欧美精品免费久久| 高清午夜精品一区二区三区| 日本一二三区视频观看| 人妻夜夜爽99麻豆av| 亚洲欧美中文字幕日韩二区| 日韩av不卡免费在线播放| 久99久视频精品免费| 91久久精品电影网| 成人av在线播放网站| 亚洲av国产av综合av卡| 美女高潮的动态| 日韩不卡一区二区三区视频在线| 国产乱人偷精品视频| 亚州av有码| 一级毛片aaaaaa免费看小| 丰满乱子伦码专区| 尾随美女入室| 亚洲aⅴ乱码一区二区在线播放|