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

    堆石料的三維應(yīng)力分?jǐn)?shù)階本構(gòu)模型

    2022-05-14 21:40:31吉華孫逸飛
    土木建筑與環(huán)境工程 2022年4期
    關(guān)鍵詞:狀態(tài)參數(shù)

    吉華 孫逸飛

    摘 要:為合理反映粗粒土的狀態(tài)依賴非關(guān)聯(lián)應(yīng)力應(yīng)變特性,提出應(yīng)力分?jǐn)?shù)階塑性力學(xué)模型。已有模型基于三軸試驗(yàn)結(jié)果,無法對堆石料真三軸條件下的應(yīng)力應(yīng)變特性進(jìn)行預(yù)測,為解決這一問題,基于特征應(yīng)力法對已有分?jǐn)?shù)階塑性力學(xué)模型進(jìn)行完善。進(jìn)一步選取不同初始狀態(tài)條件下堆石料的真三軸壓縮試驗(yàn)數(shù)據(jù)對模型進(jìn)行驗(yàn)證,結(jié)果表明,三維化后的分?jǐn)?shù)階巖土塑性力學(xué)模型可以合理地模擬不同初始狀態(tài)的堆石料在真三軸壓縮條件下的應(yīng)力應(yīng)變行為。與傳統(tǒng)塑性力學(xué)模型相比,提出的分?jǐn)?shù)階塑性模型在描述堆石料非關(guān)聯(lián)流動時(shí)不需要額外引入塑性勢函數(shù),僅需對已有屈服面求解分?jǐn)?shù)階導(dǎo)數(shù)。此外,模型在特征應(yīng)力空間中推導(dǎo)完成再映射到原應(yīng)力空間,可描述土體的三維強(qiáng)度特性,無需額外采用三維強(qiáng)度準(zhǔn)則。

    關(guān)鍵詞:堆石料;三維應(yīng)力;塑性力學(xué);分?jǐn)?shù)階微積分;狀態(tài)參數(shù);本構(gòu)模型

    中圖分類號:TU452 文獻(xiàn)標(biāo)志碼:A 文章編號:2096-6717(2022)04-0027-08

    Three-dimensional stress-fractional constitutive model for rockfill

    JI Hua1, 2, SUN Yifei3

    (1. Huadong Engineering Corporation Limited, PowerChina, Hangzhou 310014, P. R. China; 2. College of Civil and Tranportation Engineering, Hohai University, Nanjing 210098, P.R. China; 3. Faculty of Civil and Environmental Engineering, Ruhr University Bochum, Bochum 44801, Germany)

    Abstract: The stress-fractional plasticity model has been proposed to model the state-dependent stress-strain behavior of granular soils. However, the previous model was established based on triaxial test, which cannot be easily used to predict the three-dimensional (3D) stress-strain characteristics of rockfill. To solve this problem, an attempt is made to modify the fractional plasticity model for rockfill by using a novel concept of characteristic stress to make the model capable of predicting the 3D stress-strain behavior of rockfill. A series of true triaxial test results of rockfills under different initial states are simulated, from which a good agreement between the model predictions and the corresponding testing results is found. Comparing with the classical plasticity model, no additional plastic potential is required; in addition, the model can describe 3D strength characteristics by mapping from the characteristic stress space to the normal stress space, where no specific 3D strength criterion is needed.

    Keywords:rockfill; three-dimensional stress; plasticity; fractional calculus; state parameter; constitutive model

    堆石料常常被用于水利工程,諸如堆石壩、堤壩等的建造[1]。因此,在此類工程的設(shè)計(jì)、建造和運(yùn)營全過程中,堆石料的強(qiáng)度和變形特性顯得尤為重要。由于河流水荷載的復(fù)雜性,堆石料在其全壽命周期內(nèi)經(jīng)常會遭受不同應(yīng)力路徑的外部荷載作用。為了分析不同應(yīng)力路徑下堆石料的強(qiáng)度和變形行為,學(xué)者們[2-6]開展了一系列試驗(yàn)和理論研究,并建立了相應(yīng)的本構(gòu)模型。Xiao等[5]通過對不同圍壓和初始密度條件下堆石料的強(qiáng)度和剪脹行為進(jìn)行研究,提出了狀態(tài)依賴的堆石料邊界面塑性力學(xué)本構(gòu)模型。陳生水等[7]通過研究顆粒破碎對堆石料臨界狀態(tài)特性等的影響,建立了一個(gè)可以考慮堆石料顆粒破碎特性的廣義塑性力學(xué)模型。此外,在實(shí)際工程中,堆石料等粗粒土常遭受三維荷載作用,中主應(yīng)力對其強(qiáng)度和變形行為影響較大[8]。為此,Xiao等[9]提出了基于統(tǒng)一破壞準(zhǔn)則的堆石料三維邊界面塑性本構(gòu)模型;姚仰平等[10-11]、Tian等[12]提出了基于轉(zhuǎn)換應(yīng)力法[13]的軟土和粗粒土三維統(tǒng)一硬化模型。為了考慮材料的非關(guān)聯(lián)塑性流動特性,以上這些模型在建立過程中均額外引入了不同于屈服面函數(shù)的塑性勢面函數(shù)。這些塑性勢面函數(shù)往往基于經(jīng)驗(yàn)假定。那么,是否可以在不引入塑性勢的條件下建立堆石料的三維本構(gòu)模型呢?為此,Sun等[14]通過分析巖土材料的三軸試驗(yàn)特性,提出了土體的應(yīng)力分?jǐn)?shù)階塑性力學(xué)方法[14-15]。Lu等[16]基于特征應(yīng)力法,提出了軟土的三維分?jǐn)?shù)階塑性本構(gòu)模型。Qu等[17]建立了適用于巖石的分?jǐn)?shù)階塑性力學(xué)模型。李海潮等[18]進(jìn)一步引入溫度效應(yīng),拓展建立了軟巖的分?jǐn)?shù)階熱塑性模型。這些模型均處于“應(yīng)力分?jǐn)?shù)階塑性力學(xué)”框架,還有一些諸如“時(shí)間分?jǐn)?shù)階塑性力學(xué)”框架下的模型[19],仍然采用屈服面和塑性勢面不一致的假定,這里不再贅述。不同于傳統(tǒng)塑性力學(xué),應(yīng)力分?jǐn)?shù)階塑性力學(xué)基于非局域性的分?jǐn)?shù)階微分算子,在某一應(yīng)力點(diǎn)求導(dǎo)所得塑性流動方向不僅與該點(diǎn)的應(yīng)力狀態(tài)有關(guān),還與到達(dá)該點(diǎn)的加載歷史或相對臨界狀態(tài)的距離有關(guān)。在實(shí)際建模過程中,僅需對屈服面進(jìn)行一階和分?jǐn)?shù)階求導(dǎo),即可建立堆石料的非關(guān)聯(lián)彈塑性本構(gòu)模型。然而,目前提出的分?jǐn)?shù)階塑性力學(xué)模型[14, 20]僅對粗粒土的三軸應(yīng)力應(yīng)變特性進(jìn)行了分析,對土體真三軸應(yīng)力應(yīng)變特性尚未進(jìn)行預(yù)測和模擬。58794B6A-604A-455F-A538-39F6135C5C96

    為解決這一問題并同時(shí)完善現(xiàn)有分?jǐn)?shù)階巖土塑性力學(xué)理論體系,有必要建立一個(gè)三維應(yīng)力分?jǐn)?shù)階塑性力學(xué)本構(gòu)模型。筆者擬通過特征應(yīng)力法[21]對現(xiàn)有分?jǐn)?shù)階塑性力學(xué)模型進(jìn)行修正,使其可以合理地預(yù)測土體的三維應(yīng)力應(yīng)變關(guān)系。

    1 基本定義

    1.1 應(yīng)力應(yīng)變的定義

    根據(jù)臨界狀態(tài)土力學(xué)[22],約定堆石料是均勻各向同性材料;受壓為正,受拉為負(fù);模型推導(dǎo)所用到的應(yīng)力均為有效應(yīng)力。有效增量應(yīng)力張量(σ·′)定義為

    式中:i、j = 1、2、3;σ·′ij為應(yīng)力分量;ei、ej為基底矢量。相應(yīng)的增量變張量(ε·)定義為

    式中:ε·ij為應(yīng)變分量。此外,約定在主應(yīng)變空間中ε·1、ε·2和ε·3分別代表第一、第二和第三主應(yīng)變增量;相應(yīng)的主應(yīng)力空間中σ·′1、σ·′2和σ·′3分別是第一、第二和第三有效主應(yīng)力增量。

    除此之外,還有平均有效主應(yīng)力(p′)、偏應(yīng)力(q)、中主應(yīng)力系數(shù)(b)這幾個(gè)需要用到的定義。

    式中:偏應(yīng)力張量s=σ′-p′I,I為單位張量。

    1.2 分?jǐn)?shù)階導(dǎo)數(shù)

    Caputo的分?jǐn)?shù)階左導(dǎo)數(shù)和右導(dǎo)數(shù)分別定義為[23]

    式中:D =α/xα;α∈(n-1,n],是求導(dǎo)階數(shù),其中n是正整數(shù),取值范圍為1~2;σ′是加載應(yīng)力;Γ(x)=∫∞0e-ττx-1dτ。

    1.3 特征應(yīng)力法

    根據(jù)文獻(xiàn)[21],用特征應(yīng)力來代替原有強(qiáng)度準(zhǔn)則中的常規(guī)應(yīng)力,使得修正后的強(qiáng)度準(zhǔn)則可以反映土體的三維強(qiáng)度特性。特征應(yīng)力(σ′ij)定義為

    σ′ij=paσ′ijpaβ (7)

    式中:pa = 100 kPa是大氣壓強(qiáng)。β∈(0,1]是特征應(yīng)力參數(shù),反映了中主應(yīng)力比對土體強(qiáng)度的影響。根據(jù)Lu等[21]的研究,β的大小不受土體密度和圍壓等因素影響,對于給定材料,它是一個(gè)定參數(shù)。據(jù)此,修正后的特征有效平均主應(yīng)力(p′)和特征偏應(yīng)力(q)分別為

    2 三維分?jǐn)?shù)階應(yīng)力應(yīng)變關(guān)系

    增量主應(yīng)變可以分解為彈性應(yīng)變和塑性應(yīng)變兩部分,即

    式中:ε·e為彈性應(yīng)變張量;

    ε·p為塑性應(yīng)變張量,根據(jù)非線性K-G彈性理論,可得

    式中:G是剪切模量;K是體積模量,分別定義為

    式中:G0是彈性常數(shù);ν是泊松比;e是孔隙比。此外,塑性應(yīng)變張量(ε·p)可以計(jì)算為

    式中:Kp為硬化模量,其定義將在后面的本構(gòu)模型部分給出。n和m分別為加載方向和塑性流動方向,定義為[24]

    在知道屈服面的情況下,分?jǐn)?shù)階塑性力學(xué)模型無需再額外假定一個(gè)塑性勢面函數(shù),而是通過分?jǐn)?shù)階微分直接得到一個(gè)非關(guān)聯(lián)的流動法則,從而使模型得到簡化。由文獻(xiàn)[15]可知,當(dāng)α越大時(shí),所得到的非關(guān)聯(lián)性越明顯。具體參見圖1,當(dāng)α由1.05增加到1.3時(shí),同一應(yīng)力水平下加載方向與流動方向的夾角變大,即非關(guān)聯(lián)程度增加。

    3 三維分?jǐn)?shù)階塑性力學(xué)模型

    3.1 屈服面

    由傳統(tǒng)塑性力學(xué)可知,加載方向是屈服面的法向。選取特征應(yīng)力條件下的修正劍橋模型屈服面函數(shù)來描述土體的受荷屈服。

    式中:p′0是特征應(yīng)力空間屈服面與p′軸的交點(diǎn)。M是特征應(yīng)力空間的臨界狀態(tài)參數(shù),可定義為

    式中:φc是常規(guī)三軸壓縮時(shí)的臨界狀態(tài)摩擦角。如圖2所示,特征應(yīng)力法修正后的屈服面在常規(guī)主應(yīng)力空間的屈服面不再是一個(gè)圓錐,而是一個(gè)在π平面的投影,是曲邊三角形的曲面,如圖3所示。因此,可以反映土體在真三軸應(yīng)力條件下的屈服特性。

    3.2 加載張量與流動張量

    將式(19)分別代入式(17)、式(18),進(jìn)一步求得df和dg為將式(21)和式(22)代回至式(15)和式(16),即可得到塑性加載張量和塑性流動張量的具體表達(dá)式

    比較式(23)和式(24)可知,當(dāng)α=1時(shí),模型退化為傳統(tǒng)關(guān)聯(lián)流動法則。由圖4可知,對屈服面進(jìn)行分?jǐn)?shù)階求導(dǎo)后,所得子午平面上塑性流動方向不再垂直于屈服面,而是與加載方向存在一個(gè)夾角。

    此外,為了可以反映受荷變形的狀態(tài)依賴特性,Sun等[14-15]建議采用狀態(tài)依賴的分?jǐn)?shù)階數(shù)

    式中:Δ>0,是材料常數(shù)。式(25)是經(jīng)驗(yàn)公式,來源于基于試驗(yàn)現(xiàn)象的考慮:土體所處狀態(tài)會影響其塑性流動特性,而分?jǐn)?shù)階數(shù)決定了塑性流動方向。因此,運(yùn)用式(25)直觀地考慮了加載過程中土體所處狀態(tài)對其塑性流動方向的影響。ψ是狀態(tài)參數(shù),可以定義為

    式中:ec是臨界狀態(tài)孔隙比,定義為[25]

    式中:eΓ、λ和ξ是e-p′空間的臨界狀態(tài)參數(shù)。

    3.3 硬化模量

    硬化模量描述了土體的硬化特性。Sun等[14]采用Li等[25]提出的常規(guī)三軸應(yīng)力條件下的硬化模量。為了使建立的模型可以反映真三維應(yīng)力條件下的硬化特性,基于特征應(yīng)力法修正了Li等[25]的硬化模量,見式(28)。

    式中:特征峰值應(yīng)力比(Mp)定義為

    式中:h1和h2是材料常數(shù),描述了硬化參數(shù)對孔隙比的依賴;k是材料常數(shù),描述了特征峰值應(yīng)力比對土體物質(zhì)狀態(tài)的依賴。由式(28)可知,硬化模量的大小受多重因素影響(e、e0、p′、ψ、G),且在剪切尚未開始時(shí)(η=0),Kp=+∞,在臨界狀態(tài)和峰值應(yīng)力點(diǎn)時(shí),Kp=+∞。因此,所提三維硬化模量合理。為進(jìn)一步驗(yàn)證所提出的模型,選擇某堆石料的真三維試驗(yàn)數(shù)據(jù)對模型進(jìn)行驗(yàn)證。

    4 模型參數(shù)識別與模型驗(yàn)證

    4.1 模型參數(shù)識別

    模型含有11個(gè)參數(shù),分別為:φc、λ、eΓ、ξ、Δ、β、k、h1、h2、G0以及ν。這些模型參數(shù)均可以通過三軸試驗(yàn)獲得。

    參數(shù)φc、λ、eΓ和ξ為臨界狀態(tài)參數(shù)。φc可以通過擬合p′-q平面臨界狀態(tài)線的斜率得到;λ、eΓ、ξ可以通過擬合e-p′空間的臨界狀態(tài)數(shù)據(jù)得到[5]。圖5所示為臨界狀態(tài)參數(shù)對模型預(yù)測結(jié)果的影響,可見,隨著φc和eΓ數(shù)值的增加,模型預(yù)測得到的同一應(yīng)變水平下的應(yīng)力比逐漸增加;隨著λ的增加,模型預(yù)測的同一應(yīng)變水平下的應(yīng)力比逐漸減小。58794B6A-604A-455F-A538-39F6135C5C96

    參數(shù)Δ決定了土體的塑性流動特性,可以由試樣狀態(tài)轉(zhuǎn)換時(shí)的剪脹關(guān)系得到,即d=0。將式(25)代入式(22),得

    式中:ψd和ηd分別是ψ和η在狀態(tài)轉(zhuǎn)換點(diǎn)的數(shù)值。參數(shù)k可以通過試樣在峰值應(yīng)力狀態(tài)點(diǎn)的應(yīng)力應(yīng)變關(guān)系得到,即Kp=0時(shí)

    式中:ψp和ηp分別是ψ和η在峰值應(yīng)力點(diǎn)的數(shù)值。

    參數(shù)Δ和k對模型預(yù)測的影響見圖6。分析發(fā)現(xiàn):隨著參數(shù)Δ的增加,模型預(yù)測得到的同一應(yīng)變水平下的應(yīng)力比逐漸增加;隨著參數(shù)k的增加,模型預(yù)測得到的同一應(yīng)變水平下的應(yīng)力比逐漸降低。

    特征應(yīng)力參數(shù)β可由土體相關(guān)三軸壓縮(φc)和拉伸(φe)時(shí)的臨界狀態(tài)摩擦角關(guān)系式獲得[21]。β≤0.1時(shí),模型在π平面的屈服軌跡非常接近于Matsuoka-Nakai準(zhǔn)則[26]。因此,β近似取0.1。

    參數(shù)h1和h2可以通過擬合土體受荷時(shí)ε1-q關(guān)系曲線獲得。具體參見文獻(xiàn)[15, 25],這里不再贅述。圖7所示為參數(shù)h1和h2對模型預(yù)測結(jié)果的影響,分析可知:當(dāng)參數(shù)h1增加時(shí),同一應(yīng)變水平下模型預(yù)測的應(yīng)力比增加;當(dāng)參數(shù)h2增加時(shí),同一應(yīng)變水平下模型預(yù)測的應(yīng)力比減小。

    兩個(gè)彈性參數(shù)(G0、ν)可以通過擬合堆石料初始加載階段的應(yīng)力應(yīng)變曲線獲得,具體參見文獻(xiàn)[27]。由于受荷后堆石料彈性變形相對塑性變形較小,因此,不再探討彈性參數(shù)對模型預(yù)測的影響。

    4.2 模型驗(yàn)證

    選取不同初始狀態(tài)的堆石料[28]在真三軸壓縮條件下的試驗(yàn)數(shù)據(jù)進(jìn)行擬合和對比,以驗(yàn)證所提出三維分?jǐn)?shù)階巖土塑性力學(xué)模型的合理性。

    所選堆石料為兩河口心墻土石壩筑壩料,試驗(yàn)初始孔隙比為0.68。試驗(yàn)時(shí)固定第三主應(yīng)力和中主應(yīng)力比不變,第二和第一主應(yīng)力按照設(shè)定中主應(yīng)力比值的大小比例增加。其他有關(guān)材料和試樣準(zhǔn)備資料見文獻(xiàn)[28],不再贅述。根據(jù)前述方法,確定具體模型參數(shù)為:G0 = 90、ν=0.25、φc=46°、λ=0.11、eΓ=0.404、ξ=0.1、Δ=0.2、β = 0.1、k = 0.1、h1 = 1.2、h2 = 0.3。

    圖8~圖10為所提出分?jǐn)?shù)階塑性力學(xué)模型對真三軸試驗(yàn)數(shù)據(jù)的模擬。對比分析可以發(fā)現(xiàn):模型可以合理地模擬堆石料在不同圍壓和中主應(yīng)力比條件下的真三維應(yīng)力應(yīng)變特性。具體來說,隨著第一主應(yīng)變的增加,試樣的應(yīng)力比不斷增加;第二主應(yīng)變只有在中主應(yīng)力比等于0時(shí),即常規(guī)三軸試驗(yàn)條件下,數(shù)值為負(fù),其他中主應(yīng)力比條件下均為正;第三主應(yīng)變在任何中主應(yīng)力比條件下均為負(fù)值。這些應(yīng)力應(yīng)變特性均可被較好地模擬。

    此外,需要指出的是,中主應(yīng)力系數(shù)較小時(shí),圖8~圖10中試樣在同一偏應(yīng)力水平所產(chǎn)生的第一主應(yīng)變大于第二、第三主應(yīng)變,但隨著中主應(yīng)力系數(shù)的增加,第一主應(yīng)變小于第二、第三主應(yīng)變;這是因?yàn)殡S著中主應(yīng)力系數(shù)從0變化到1,原先的第二主應(yīng)力大小慢慢接近于第一主應(yīng)力,對側(cè)向應(yīng)變的主導(dǎo)作用逐漸變大。從而在第一主應(yīng)力方向產(chǎn)生的應(yīng)變會逐漸小于同等偏應(yīng)力水平下其他兩個(gè)主應(yīng)力方向的應(yīng)變。

    5 結(jié)論

    基于特征應(yīng)力法,修正了已提出的分?jǐn)?shù)階應(yīng)力塑性力學(xué)模型,使其可以對堆石料在一般三維應(yīng)力條件下的應(yīng)力應(yīng)變進(jìn)行預(yù)測。并選擇堆石料的真三軸壓縮試驗(yàn)數(shù)據(jù),對模型進(jìn)了驗(yàn)證。主要結(jié)論如下:

    1)基于特征應(yīng)力修正的屈服面,其在主應(yīng)力空間不再是一個(gè)圓錐,而是一個(gè)可以反映土體三維屈服特性的曲面,該面在π平面的投影為曲邊三角形。

    2)無需額外假定塑性勢面,只需要對已有屈服面進(jìn)行分?jǐn)?shù)階微分求導(dǎo),便可以得到非關(guān)聯(lián)的流動法則。進(jìn)一步得出剪脹方程受分?jǐn)?shù)階求導(dǎo)階數(shù)、特征應(yīng)力參數(shù)以及中主應(yīng)力系數(shù)的影響。

    3)基于特征應(yīng)力法修正了已有常規(guī)三軸應(yīng)力條件下的硬化模量,使其可以反映三維應(yīng)力條件下堆石料的硬化特性。

    4)對比試驗(yàn)數(shù)據(jù),驗(yàn)證了所提模型的合理性;模型可以對不同初始狀態(tài)試樣在不同中主應(yīng)力比條件下的真三軸應(yīng)力應(yīng)變特性進(jìn)行較好的模擬。

    5)建立模型僅考慮了完全飽和堆石料的變形特性,對于干燥或半干燥堆石料遇水飽和后存在的濕化或者軟化性能還有待進(jìn)一步研究。此外,模型也未考慮堆石料在填筑或者運(yùn)維過程中的顆粒破碎,對于考慮顆粒破碎影響的堆石料三維分?jǐn)?shù)階塑性模型也有待進(jìn)一步研究。

    參考文獻(xiàn):

    [1] 劉漢龍, 孫逸飛, 楊貴, 等. 粗粒料顆粒破碎特性研究述評[J]. 河海大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 40(4): 361-369.

    LIU H L, SUN Y F, YANG G, et al. A review of particle breakage characteristics of coarse aggregates [J]. Journal of Hohai University (Natural Sciences), 2012, 40(4): 361-369. (in Chinese)

    [2] 劉漢龍, 秦紅玉, 高玉峰, 等. 堆石粗粒料顆粒破碎試驗(yàn)研究[J]. 巖土力學(xué), 2005, 26(4): 562-566.

    LIU H L, QIN H Y, GAO Y F, et al. Experimental study on particle breakage of rockfill and coarse aggregates [J]. Rock and Soil Mechanics, 2005, 26(4): 562-566. (in Chinese)

    [3] 賈宇峰, 王丙申, 遲世春. 堆石料剪切過程中的顆粒破碎研究[J]. 巖土工程學(xué)報(bào), 2015, 37(9): 1692-1697.

    JIA Y F, WANG B S, CHI S C. Particle breakage of rockfill during triaxial tests [J]. Chinese Journal of Geotechnical Engineering, 2015, 37(9): 1692-1697. (in Chinese)58794B6A-604A-455F-A538-39F6135C5C96

    [4] XIAO Y, LIU H L, CHEN Y M, et al. Strength and deformation of rockfill material based on large-scale triaxial compression tests. I: Influences of density and pressure [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2014, 140(12): 04014070.

    [5] XIAO Y, LIU H L, CHEN Y M, et al. Bounding surface model for rockfill materials dependent on density and pressure under triaxial stress conditions [J]. Journal of Engineering Mechanics, 2014, 140(4): 04014002.

    [6] DAFALIAS Y F, TAIEBAT M. SANISAND-Z: Zero elastic range sand plasticity model [J]. Géotechnique, 2016, 66(12): 999-1013.

    [7] 陳生水, 傅中志, 韓華強(qiáng), 等. 一個(gè)考慮顆粒破碎的堆石料彈塑性本構(gòu)模型[J]. 巖土工程學(xué)報(bào), 2011, 33(10): 1489-1495.

    CHEN S S, FU Z Z, HAN H Q, et al. An elastoplastic model for rockfill materials considering particle breakage [J]. Chinese Journal of Geotechnical Engineering, 2011, 33(10): 1489-1495. (in Chinese)

    [8] XIAO Y, LIU H L, CHEN Y M, et al. Influence of intermediate principal stress on the strength and dilatancy behavior of rockfill material [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2014, 140(11): 04014064.

    [9] XIAO Y, LIU H L, ZHU J G, et al. Modeling and behaviours of rockfill materials in three-dimensional stress space [J]. Science China Technological Sciences, 2012, 55(10): 2877-2892.

    [10] YAO Y P, YAMAMOTO H, WANG N D. Constitutive model considering sand crushing [J]. Soils and Foundations, 2008, 48(4): 603-608.

    [11] 姚仰平, 劉林, 羅汀. 砂土的UH模型[J]. 巖土工程學(xué)報(bào), 2016, 38(12): 2147-2153.

    YAO Y P, LIU L, LUO T. UH model for sands [J]. Chinese Journal of Geotechnical Engineering, 2016, 38(12): 2147-2153. (in Chinese)

    [12] TIAN Y, YAO Y P. Modelling the non-coaxiality of soils from the view of cross-anisotropy [J]. Computers and Geotechnics, 2017, 86: 219-229.

    [13] YAO Y P, WANG N D. Transformed stress method for generalizing soil constitutive models [J]. Journal of Engineering Mechanics, 2014, 140(3): 614-629.

    [14] SUN Y F, XIAO Y. Fractional order plasticity model for granular soils subjected to monotonic triaxial compression [J]. International Journal of Solids and Structures, 2017, 118/119: 224-234.

    [15] SUN Y F, SHEN Y. Constitutive model of granular soils using fractional-order plastic-flow rule [J]. International Journal of Geomechanics, 2017, 17(8): 04017025.

    [16] LU D C, LIANG J Y, DU X L, et al. Fractional elastoplastic constitutive model for soils based on a novel 3D fractional plastic flow rule [J]. Computers and Geotechnics, 2019, 105: 277-290.58794B6A-604A-455F-A538-39F6135C5C96

    [17] QU P F, ZHU Q Z, SUN Y F. Elastoplastic modelling of mechanical behavior of rocks with fractional-order plastic flow [J]. International Journal of Mechanical Sciences, 2019, 163: 105102.

    [18] 李海潮, 馬博, 張升, 等. 基于分?jǐn)?shù)階熱彈塑性理論的軟巖力學(xué)特性描述[J]. 巖石力學(xué)與工程學(xué)報(bào), 2020, 39(7): 1311-1320.

    LI H C, MA B, ZHANG S, et al. Mechanical behaviors of soft rocks based on the fractional thermal elastic-plastic theory [J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(7): 1311-1320. (in Chinese)

    [19] SUN Y, XIAO Y. Fractional order model for granular soils under drained cyclic loading [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2017, 41(4): 555-577.

    [20] SUN Y, GAO Y, SHEN Y. Mathematical aspect of the state-dependent stress-dilatancy of granular soil under triaxial loading [J]. Géotechnique, 2019, 69(2): 158-165.

    [21] LU D C, MA C, DU X L, et al. Development of a new nonlinear unified strength theory for geomaterials based on the characteristic stress concept [J]. International Journal of Geomechanics, 2017, 17(2): 04016058.

    [22] SCHOFIELD A, WROTH P. Critical state soil mechanics [M]. New York, USA: McGraw-Hill London, 1968.

    [23] PODLUBNY I. Fractional differential equations: An introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications [M]. San Diego, California: Academic press, 1998.

    [24] SUN Y F, NIMBALKAR S. Stress-fractional soil model with reduced elastic region [J]. Soils and Foundations, 2019, 59(6): 2007-2023.

    [25] LI X S, DAFALIAS Y F. Dilatancy for cohesionless soils [J]. Géotechnique, 2000, 50(4): 449-460.

    [26] MATSUOKA H, NAKAI T R. Stress-deformation and strength characteristics of soil under three different principal stresses [J]. Proceedings of the Japan Society of Civil Engineers, 1974, 232: 59-70.

    [27] SUN Y F, WICHTMANN T, SUMELKA W, et al. Karlsruhe fine sand under monotonic and cyclic loads: Modelling and validation [J]. Soil Dynamics and Earthquake Engineering, 2020, 133: 106119.

    [28] XIAO Y, LIU H L, ZHU J G, et al. A 3D bounding surface model for rockfill materials [J]. Science China Technological Sciences, 2011, 54(11): 2904-2915.

    (編輯 王秀玲)

    收稿日期:2020-07-23

    基金項(xiàng)目:中央高?;究蒲袠I(yè)務(wù)費(fèi)(2017B05214);中國博士后科學(xué)基金(2017M621607)

    作者簡介:吉華(1987- ),博士(后),主要從事巖土材料力學(xué)特性與本構(gòu)理論研究,E-mail:ji_h@ecidi.com。

    孫逸飛(通信作者),博士,E-mail:yifei.sun@rub.de。

    Received:2020-07-23

    Foundation items:Fundamental Research Funds for the Central Universities (No. 2017B05214); China Postdoctoral Science Foundation (No. 2017M621607)

    Author brief:JI Hua (1987- ), postdoctoral fellow, main research interests: mechanical properties and constitutive theory of geomaterials, E-mail: ji_h@ecidi.com.

    SUN Yifei (corresponding author), PhD, E-mail: yifei.sun@rub.de.58794B6A-604A-455F-A538-39F6135C5C96

    猜你喜歡
    狀態(tài)參數(shù)
    基于Transformer的飛機(jī)狀態(tài)預(yù)測
    空壓機(jī)運(yùn)行狀態(tài)參數(shù)智能監(jiān)測方法設(shè)計(jì)
    基于CKF的大型拖拉機(jī)狀態(tài)參數(shù)估計(jì)研究
    基于HHT算法的分布式電力系統(tǒng)狀態(tài)自動監(jiān)測方法
    中國農(nóng)機(jī)院科研成果:飼草料作物收獲技術(shù)與裝備研發(fā)+農(nóng)田提質(zhì)工程技術(shù)與裝備研發(fā)+機(jī)器作業(yè)狀態(tài)參數(shù)測試方法研究項(xiàng)目
    鋰電池異常狀態(tài)監(jiān)測的FEKF方法研究
    適用于計(jì)算機(jī)程序計(jì)算的濕空氣計(jì)算方法
    基于Visual Basic 6.0制冷劑狀態(tài)參數(shù)模擬的比較與分析
    居業(yè)(2019年3期)2019-06-07 10:54:51
    基于Android操作系統(tǒng)的智能手機(jī)下濕空氣參數(shù)計(jì)算
    基于大數(shù)據(jù)風(fēng)電場狀態(tài)參數(shù)相關(guān)性分析
    国产精品久久久久久精品古装| videosex国产| 考比视频在线观看| 国产精品国产av在线观看| 你懂的网址亚洲精品在线观看| 一级,二级,三级黄色视频| www.色视频.com| av国产精品久久久久影院| 少妇人妻精品综合一区二区| 久久久精品免费免费高清| 波野结衣二区三区在线| 免费观看a级毛片全部| av天堂久久9| 曰老女人黄片| 中文字幕久久专区| 欧美xxⅹ黑人| 国产成人精品婷婷| 午夜激情av网站| 亚洲第一av免费看| 久久久久久久久久久久大奶| 国产片内射在线| 亚洲国产色片| 国产精品国产av在线观看| 交换朋友夫妻互换小说| 有码 亚洲区| 热99国产精品久久久久久7| 国产亚洲精品第一综合不卡 | 久久精品人人爽人人爽视色| 国产一区二区在线观看日韩| 美女福利国产在线| av在线app专区| 亚洲av二区三区四区| 欧美日韩精品成人综合77777| 黄片播放在线免费| 欧美亚洲日本最大视频资源| 十八禁网站网址无遮挡| 草草在线视频免费看| 久久国产亚洲av麻豆专区| 午夜激情福利司机影院| 日韩免费高清中文字幕av| 国产精品久久久久久av不卡| 大码成人一级视频| 欧美人与善性xxx| 精品亚洲成国产av| 18禁观看日本| 亚洲美女搞黄在线观看| 精品亚洲成国产av| a级片在线免费高清观看视频| 国产精品99久久久久久久久| 日本免费在线观看一区| 制服丝袜香蕉在线| 亚洲一区二区三区欧美精品| 如日韩欧美国产精品一区二区三区 | 亚洲国产精品成人久久小说| a级毛片黄视频| 曰老女人黄片| 国产亚洲一区二区精品| 人妻少妇偷人精品九色| 卡戴珊不雅视频在线播放| 亚洲内射少妇av| 免费人妻精品一区二区三区视频| 亚洲av国产av综合av卡| 国产熟女午夜一区二区三区 | 亚洲少妇的诱惑av| 久久婷婷青草| 亚洲人成网站在线观看播放| 18禁裸乳无遮挡动漫免费视频| 国产成人精品一,二区| 国产色婷婷99| 成人免费观看视频高清| 国产精品国产三级专区第一集| 大码成人一级视频| 国产成人a∨麻豆精品| 中文乱码字字幕精品一区二区三区| 欧美成人精品欧美一级黄| 99久久人妻综合| 亚洲精华国产精华液的使用体验| 久久久久久久久久人人人人人人| 成人亚洲精品一区在线观看| 亚洲精品乱码久久久久久按摩| 美女国产视频在线观看| 丰满迷人的少妇在线观看| 成人影院久久| 嫩草影院入口| 熟女人妻精品中文字幕| 人妻少妇偷人精品九色| 日韩一区二区三区影片| 国产成人91sexporn| 久久久久久久久久久丰满| 丝瓜视频免费看黄片| 新久久久久国产一级毛片| 又黄又爽又刺激的免费视频.| av在线老鸭窝| 日韩电影二区| 免费人成在线观看视频色| av视频免费观看在线观看| 老司机影院毛片| 大香蕉久久成人网| 18禁裸乳无遮挡动漫免费视频| 国产高清有码在线观看视频| 国产片特级美女逼逼视频| 九九爱精品视频在线观看| 另类亚洲欧美激情| 免费av不卡在线播放| 91国产中文字幕| 欧美精品国产亚洲| 少妇猛男粗大的猛烈进出视频| 水蜜桃什么品种好| 一区二区日韩欧美中文字幕 | 老司机影院毛片| 97超视频在线观看视频| 久久精品国产亚洲网站| 久久精品人人爽人人爽视色| 久久久久视频综合| 久久精品国产亚洲网站| 精品酒店卫生间| 最新中文字幕久久久久| 亚洲av免费高清在线观看| 午夜免费观看性视频| av播播在线观看一区| 久久国产亚洲av麻豆专区| 涩涩av久久男人的天堂| 亚洲一区二区三区欧美精品| 免费观看性生交大片5| 自线自在国产av| 99久久精品一区二区三区| 精品久久久噜噜| 在线观看国产h片| 免费观看性生交大片5| 国产精品久久久久成人av| 晚上一个人看的免费电影| 国产成人精品一,二区| 人妻一区二区av| 精品酒店卫生间| 精品国产乱码久久久久久小说| 国产男女内射视频| 成人国产av品久久久| 亚洲av二区三区四区| 久久99蜜桃精品久久| 精品午夜福利在线看| videossex国产| 欧美亚洲日本最大视频资源| 国产午夜精品一二区理论片| 交换朋友夫妻互换小说| 成人国产av品久久久| 在线精品无人区一区二区三| 女人精品久久久久毛片| 国产视频内射| 黄色一级大片看看| 亚洲精品色激情综合| 午夜激情av网站| 人人妻人人爽人人添夜夜欢视频| 97超视频在线观看视频| 国产一级毛片在线| 国产成人精品无人区| 久久久精品94久久精品| av线在线观看网站| 人成视频在线观看免费观看| 国产免费视频播放在线视频| 另类亚洲欧美激情| 一区二区av电影网| 午夜激情av网站| 国产在线视频一区二区| 男女边吃奶边做爰视频| 亚洲精品一区蜜桃| 只有这里有精品99| 99久国产av精品国产电影| 熟女av电影| 国产一区二区在线观看日韩| 欧美xxxx性猛交bbbb| 久久久久国产网址| 国产精品三级大全| 欧美xxxx性猛交bbbb| 丰满饥渴人妻一区二区三| 爱豆传媒免费全集在线观看| 男女边吃奶边做爰视频| 久久午夜综合久久蜜桃| 美女国产高潮福利片在线看| 十八禁高潮呻吟视频| 国产精品人妻久久久久久| 伦理电影大哥的女人| videosex国产| 大陆偷拍与自拍| av有码第一页| 久久国产精品大桥未久av| 汤姆久久久久久久影院中文字幕| 婷婷色综合大香蕉| 中国三级夫妇交换| 伦理电影免费视频| 亚洲人与动物交配视频| 成年女人在线观看亚洲视频| 欧美日韩一区二区视频在线观看视频在线| 色婷婷av一区二区三区视频| 欧美 亚洲 国产 日韩一| 亚洲人成网站在线观看播放| 国产极品天堂在线| 欧美另类一区| 日韩 亚洲 欧美在线| 精品视频人人做人人爽| 在线看a的网站| 大码成人一级视频| 免费日韩欧美在线观看| 国产亚洲精品第一综合不卡 | 国产在线视频一区二区| 久久精品久久久久久噜噜老黄| 美女福利国产在线| 伊人久久国产一区二区| 美女脱内裤让男人舔精品视频| 丰满乱子伦码专区| 国产高清三级在线| 久久久久精品久久久久真实原创| 国产视频内射| 国产精品 国内视频| 午夜福利视频精品| 少妇的逼水好多| 精品国产一区二区久久| 亚洲精品中文字幕在线视频| 国产有黄有色有爽视频| 91在线精品国自产拍蜜月| 美女中出高潮动态图| av福利片在线| 麻豆乱淫一区二区| 看十八女毛片水多多多| 91aial.com中文字幕在线观看| av免费在线看不卡| 一区二区三区免费毛片| 久久午夜福利片| 97在线视频观看| 99热全是精品| 免费看光身美女| 3wmmmm亚洲av在线观看| 18禁在线无遮挡免费观看视频| 久久久久久久亚洲中文字幕| 狠狠婷婷综合久久久久久88av| 亚洲内射少妇av| 97精品久久久久久久久久精品| 欧美老熟妇乱子伦牲交| 五月伊人婷婷丁香| 91久久精品国产一区二区成人| 高清毛片免费看| 熟女电影av网| 欧美少妇被猛烈插入视频| 韩国高清视频一区二区三区| 五月开心婷婷网| 色吧在线观看| 狂野欧美激情性bbbbbb| 精品人妻熟女毛片av久久网站| 精品久久久久久电影网| 精品久久久精品久久久| 成人午夜精彩视频在线观看| 日本与韩国留学比较| 最后的刺客免费高清国语| 国产精品麻豆人妻色哟哟久久| 日韩人妻高清精品专区| 精品久久蜜臀av无| 欧美精品亚洲一区二区| 亚洲欧美成人精品一区二区| 在线看a的网站| 国产男人的电影天堂91| 午夜福利影视在线免费观看| 国产淫语在线视频| 亚洲一区二区三区欧美精品| 91在线精品国自产拍蜜月| 99热6这里只有精品| 天天影视国产精品| 99视频精品全部免费 在线| 亚洲伊人久久精品综合| 久久 成人 亚洲| 伊人久久国产一区二区| 免费观看在线日韩| 一边亲一边摸免费视频| 亚洲五月色婷婷综合| 欧美激情国产日韩精品一区| 少妇丰满av| 老司机影院成人| 久久久欧美国产精品| 国产免费一级a男人的天堂| 日韩制服骚丝袜av| 久久精品久久精品一区二区三区| 18禁裸乳无遮挡动漫免费视频| 精品久久久精品久久久| 男女无遮挡免费网站观看| 亚洲激情五月婷婷啪啪| av国产精品久久久久影院| 国产日韩欧美视频二区| 国产又色又爽无遮挡免| 国产精品久久久久久久电影| 日韩av免费高清视频| 我要看黄色一级片免费的| 精品亚洲成a人片在线观看| 在线精品无人区一区二区三| 亚洲精品自拍成人| 国产高清国产精品国产三级| 欧美 亚洲 国产 日韩一| 久久久久久久久久久丰满| 久久热精品热| 简卡轻食公司| 久久99一区二区三区| 18禁在线无遮挡免费观看视频| 亚洲久久久国产精品| 欧美少妇被猛烈插入视频| 天天影视国产精品| xxx大片免费视频| 亚洲av在线观看美女高潮| 国产精品一国产av| 丰满饥渴人妻一区二区三| 亚洲av欧美aⅴ国产| 亚洲精品国产av成人精品| 五月开心婷婷网| 色视频在线一区二区三区| 国产午夜精品一二区理论片| 男女高潮啪啪啪动态图| 少妇的逼水好多| 大片电影免费在线观看免费| 欧美精品高潮呻吟av久久| 日韩不卡一区二区三区视频在线| 纯流量卡能插随身wifi吗| 国产精品无大码| 97精品久久久久久久久久精品| 满18在线观看网站| 能在线免费看毛片的网站| 日韩熟女老妇一区二区性免费视频| 日本爱情动作片www.在线观看| 日日摸夜夜添夜夜添av毛片| 国产精品一区二区在线不卡| 日日摸夜夜添夜夜添av毛片| 国产精品女同一区二区软件| 亚洲,欧美,日韩| 中国国产av一级| 80岁老熟妇乱子伦牲交| 伦理电影免费视频| 欧美老熟妇乱子伦牲交| 久久综合国产亚洲精品| 欧美人与善性xxx| 日韩中文字幕视频在线看片| 亚洲欧美日韩另类电影网站| 在线精品无人区一区二区三| 蜜臀久久99精品久久宅男| 久久午夜综合久久蜜桃| 纯流量卡能插随身wifi吗| 欧美最新免费一区二区三区| 少妇人妻 视频| 五月天丁香电影| 两个人免费观看高清视频| 欧美人与性动交α欧美精品济南到 | 久久狼人影院| 国产精品免费大片| 午夜免费观看性视频| 亚洲内射少妇av| 少妇精品久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品麻豆人妻色哟哟久久| 99久久中文字幕三级久久日本| 高清毛片免费看| 视频区图区小说| 中文字幕亚洲精品专区| 国产欧美日韩综合在线一区二区| 女的被弄到高潮叫床怎么办| 搡女人真爽免费视频火全软件| 18禁观看日本| www.色视频.com| 久久精品国产鲁丝片午夜精品| 日韩不卡一区二区三区视频在线| 亚洲少妇的诱惑av| 亚洲精品乱久久久久久| 欧美97在线视频| 精品久久蜜臀av无| 亚洲精品久久成人aⅴ小说 | 午夜影院在线不卡| 亚洲一级一片aⅴ在线观看| 国产乱人偷精品视频| 草草在线视频免费看| av线在线观看网站| 97超视频在线观看视频| 精品视频人人做人人爽| 亚洲欧洲精品一区二区精品久久久 | 日韩一本色道免费dvd| 下体分泌物呈黄色| 2022亚洲国产成人精品| 国产一级毛片在线| 精品国产一区二区久久| videos熟女内射| 成人国产av品久久久| 亚洲国产精品专区欧美| 日本黄色片子视频| 日韩人妻高清精品专区| 一级黄片播放器| av.在线天堂| 黄色怎么调成土黄色| 91精品一卡2卡3卡4卡| 国产精品女同一区二区软件| 黄色一级大片看看| 性色avwww在线观看| 狠狠婷婷综合久久久久久88av| 亚洲四区av| 三级国产精品欧美在线观看| 久久精品夜色国产| 国产一区二区三区av在线| 国产精品久久久久久久久免| 一区二区av电影网| 国产高清三级在线| 夫妻午夜视频| 伦精品一区二区三区| 日韩成人av中文字幕在线观看| 人妻制服诱惑在线中文字幕| 午夜91福利影院| 欧美三级亚洲精品| 51国产日韩欧美| 最新的欧美精品一区二区| 日本黄色片子视频| 视频区图区小说| av黄色大香蕉| 最近最新中文字幕免费大全7| 国产日韩欧美在线精品| 伦理电影免费视频| 亚洲国产日韩一区二区| 能在线免费看毛片的网站| 亚洲欧美精品自产自拍| 少妇 在线观看| 成人手机av| 欧美日韩亚洲高清精品| 菩萨蛮人人尽说江南好唐韦庄| 午夜福利视频在线观看免费| 一级a做视频免费观看| 午夜免费男女啪啪视频观看| 久久精品久久精品一区二区三区| 亚洲成人av在线免费| 2021少妇久久久久久久久久久| 交换朋友夫妻互换小说| 另类精品久久| 精品一区二区免费观看| 中文字幕久久专区| 丰满迷人的少妇在线观看| 看免费成人av毛片| 婷婷色综合www| av网站免费在线观看视频| 少妇 在线观看| 欧美日韩视频精品一区| 亚洲国产精品国产精品| 亚洲综合色网址| 国产精品人妻久久久久久| 久久99热这里只频精品6学生| 国产精品一区二区在线观看99| 大香蕉97超碰在线| 亚洲av国产av综合av卡| 国产欧美日韩一区二区三区在线 | 两个人免费观看高清视频| 欧美性感艳星| 久久国产精品大桥未久av| 在现免费观看毛片| 91国产中文字幕| 蜜桃久久精品国产亚洲av| 亚洲一级一片aⅴ在线观看| h视频一区二区三区| 久久人人爽av亚洲精品天堂| 夫妻性生交免费视频一级片| 精品一区二区三区视频在线| 中文欧美无线码| 亚洲图色成人| 国产成人免费观看mmmm| 人妻 亚洲 视频| 激情五月婷婷亚洲| 日本黄色日本黄色录像| 亚洲欧美清纯卡通| 视频在线观看一区二区三区| 特大巨黑吊av在线直播| 五月玫瑰六月丁香| 色5月婷婷丁香| 午夜福利视频精品| 国产免费一级a男人的天堂| 麻豆精品久久久久久蜜桃| 国产国语露脸激情在线看| 在线观看一区二区三区激情| 中文字幕最新亚洲高清| 国产精品一区二区三区四区免费观看| 欧美日韩国产mv在线观看视频| 永久网站在线| 免费观看av网站的网址| 中文字幕精品免费在线观看视频 | 久久久久人妻精品一区果冻| 五月开心婷婷网| 99久国产av精品国产电影| 汤姆久久久久久久影院中文字幕| 欧美成人午夜免费资源| 午夜久久久在线观看| 伊人久久精品亚洲午夜| 我要看黄色一级片免费的| 这个男人来自地球电影免费观看 | 欧美精品国产亚洲| 久久久久久久久久成人| 国产在视频线精品| 亚洲精品视频女| 飞空精品影院首页| 五月玫瑰六月丁香| 99国产精品免费福利视频| 亚洲国产最新在线播放| 日韩一本色道免费dvd| 亚洲av不卡在线观看| 亚州av有码| 久久亚洲国产成人精品v| 久久国产亚洲av麻豆专区| 高清黄色对白视频在线免费看| 久久99热6这里只有精品| 精品一区二区免费观看| 热99久久久久精品小说推荐| 精品卡一卡二卡四卡免费| 边亲边吃奶的免费视频| 国产精品蜜桃在线观看| 亚洲国产欧美在线一区| 久久精品国产亚洲av涩爱| 美女cb高潮喷水在线观看| 欧美亚洲 丝袜 人妻 在线| 自线自在国产av| 一级a做视频免费观看| 人人妻人人爽人人添夜夜欢视频| 啦啦啦中文免费视频观看日本| 岛国毛片在线播放| 成人无遮挡网站| 丝袜喷水一区| 日本-黄色视频高清免费观看| 亚洲综合精品二区| 国产成人精品久久久久久| 亚洲欧美清纯卡通| 久久人妻熟女aⅴ| 啦啦啦中文免费视频观看日本| 成人免费观看视频高清| 一本大道久久a久久精品| 久久久国产精品麻豆| 大片免费播放器 马上看| 国产有黄有色有爽视频| 中文乱码字字幕精品一区二区三区| 精品熟女少妇av免费看| av在线播放精品| 亚洲中文av在线| 成人无遮挡网站| 亚洲欧美精品自产自拍| 精品人妻熟女av久视频| 国产精品一区二区在线观看99| 人妻人人澡人人爽人人| 亚洲国产精品一区二区三区在线| 一本一本综合久久| 人成视频在线观看免费观看| 少妇 在线观看| 国产黄频视频在线观看| 飞空精品影院首页| 九九久久精品国产亚洲av麻豆| 国产日韩欧美亚洲二区| 国产一区亚洲一区在线观看| 精品视频人人做人人爽| 久久ye,这里只有精品| 99国产精品免费福利视频| 少妇猛男粗大的猛烈进出视频| 美女脱内裤让男人舔精品视频| 热99国产精品久久久久久7| av网站免费在线观看视频| 国产乱来视频区| 一区在线观看完整版| 黄色视频在线播放观看不卡| 午夜日本视频在线| 91精品国产九色| 久久久久久久久大av| 99久久精品一区二区三区| 在线观看美女被高潮喷水网站| 岛国毛片在线播放| 午夜老司机福利剧场| 亚洲精品aⅴ在线观看| 少妇 在线观看| 街头女战士在线观看网站| 18禁动态无遮挡网站| 亚洲av.av天堂| 又黄又爽又刺激的免费视频.| 精品久久久久久久久av| 久久久欧美国产精品| 久久青草综合色| 最近中文字幕2019免费版| 99久国产av精品国产电影| 亚洲精品乱码久久久久久按摩| 另类精品久久| 日韩制服骚丝袜av| 久久国产精品男人的天堂亚洲 | 国产成人一区二区在线| 午夜激情久久久久久久| 乱码一卡2卡4卡精品| 亚州av有码| 午夜视频国产福利| 青青草视频在线视频观看| 精品人妻熟女毛片av久久网站| av卡一久久| 在线看a的网站| 日韩一区二区三区影片| 最黄视频免费看| 久久久a久久爽久久v久久| 国产精品欧美亚洲77777| 欧美激情国产日韩精品一区| 精品卡一卡二卡四卡免费| 男女免费视频国产| 天堂俺去俺来也www色官网| 日本色播在线视频| 久久久久久久久久久丰满| 午夜激情久久久久久久| 久久毛片免费看一区二区三区| 日本爱情动作片www.在线观看| 老司机影院毛片| 久久久久精品性色| 亚洲美女黄色视频免费看| 久久久久久久久久久丰满| 亚洲内射少妇av| 在线播放无遮挡| 99久久综合免费| 我的老师免费观看完整版| 国产日韩欧美视频二区| 视频区图区小说| 成人综合一区亚洲| 久久综合国产亚洲精品| 人人妻人人澡人人爽人人夜夜|