• <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)性分析
    99久久综合免费| 久久久久久久亚洲中文字幕| 美女高潮到喷水免费观看| 中文字幕人妻丝袜制服| 成人国语在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美在线黄色| 99久久综合免费| 高清欧美精品videossex| 亚洲av免费高清在线观看| www日本在线高清视频| 下体分泌物呈黄色| 成年人午夜在线观看视频| 超碰成人久久| 亚洲精品日本国产第一区| 97精品久久久久久久久久精品| 亚洲成人一二三区av| 免费黄网站久久成人精品| av在线app专区| 久久国产精品男人的天堂亚洲| 人妻 亚洲 视频| 人妻一区二区av| 久久鲁丝午夜福利片| 亚洲经典国产精华液单| 亚洲国产精品一区二区三区在线| 大香蕉久久成人网| 国产一区二区 视频在线| 国产黄色视频一区二区在线观看| 赤兔流量卡办理| 人妻系列 视频| 青春草视频在线免费观看| 一级毛片黄色毛片免费观看视频| 免费在线观看黄色视频的| 日韩成人av中文字幕在线观看| 一二三四在线观看免费中文在| 男人添女人高潮全过程视频| a级毛片在线看网站| 亚洲av男天堂| 最近手机中文字幕大全| 青春草视频在线免费观看| 成年美女黄网站色视频大全免费| 免费久久久久久久精品成人欧美视频| 久久久久网色| 国产不卡av网站在线观看| 久久久亚洲精品成人影院| 国产亚洲欧美精品永久| 国产片特级美女逼逼视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国产网址| 亚洲伊人久久精品综合| 久久国产亚洲av麻豆专区| 久久久久久人人人人人| 国产男女超爽视频在线观看| 在线观看三级黄色| 精品少妇一区二区三区视频日本电影 | 永久免费av网站大全| 国产精品亚洲av一区麻豆 | 亚洲伊人色综图| 五月开心婷婷网| 五月伊人婷婷丁香| 一本久久精品| 精品一区二区三卡| 国产成人精品久久二区二区91 | 黑人欧美特级aaaaaa片| 欧美日韩亚洲高清精品| 熟女电影av网| 亚洲精品国产一区二区精华液| 中文字幕精品免费在线观看视频| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 久久久久人妻精品一区果冻| 久久久久久久国产电影| 一级毛片黄色毛片免费观看视频| 免费久久久久久久精品成人欧美视频| 国产亚洲一区二区精品| 久久久久久久亚洲中文字幕| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黄色配什么色好看| 一级毛片黄色毛片免费观看视频| 久久99一区二区三区| a级毛片在线看网站| 最新中文字幕久久久久| 在线观看三级黄色| 成人18禁高潮啪啪吃奶动态图| 久久99热这里只频精品6学生| 亚洲精品国产av成人精品| 精品少妇一区二区三区视频日本电影 | 自拍欧美九色日韩亚洲蝌蚪91| 色视频在线一区二区三区| 天天躁日日躁夜夜躁夜夜| 国产精品免费视频内射| 国产在线一区二区三区精| 人妻 亚洲 视频| 老司机影院毛片| 国产精品 国内视频| 亚洲激情五月婷婷啪啪| 国产成人免费观看mmmm| 国精品久久久久久国模美| 欧美日韩精品网址| 欧美人与善性xxx| 69精品国产乱码久久久| 性色avwww在线观看| 久久久a久久爽久久v久久| 天堂中文最新版在线下载| 中文精品一卡2卡3卡4更新| 亚洲欧美一区二区三区国产| 日韩在线高清观看一区二区三区| 久久久久精品久久久久真实原创| 99久国产av精品国产电影| 亚洲国产看品久久| 欧美精品一区二区免费开放| 一级毛片我不卡| 午夜福利视频在线观看免费| 婷婷色麻豆天堂久久| 青青草视频在线视频观看| 久久久久久人人人人人| 最近最新中文字幕大全免费视频 | 国产欧美日韩综合在线一区二区| 国产国语露脸激情在线看| 国产淫语在线视频| 午夜影院在线不卡| 两个人免费观看高清视频| 国产不卡av网站在线观看| 人人妻人人澡人人爽人人夜夜| 中文天堂在线官网| 国产麻豆69| 综合色丁香网| 曰老女人黄片| 九九爱精品视频在线观看| 亚洲视频免费观看视频| 久久人妻熟女aⅴ| 亚洲人成77777在线视频| 青草久久国产| 欧美日韩亚洲国产一区二区在线观看 | √禁漫天堂资源中文www| 国产成人aa在线观看| 亚洲男人天堂网一区| 人人妻人人澡人人看| 一本久久精品| a级片在线免费高清观看视频| 国产 一区精品| 久久精品国产综合久久久| 亚洲欧美色中文字幕在线| 日本av手机在线免费观看| 精品99又大又爽又粗少妇毛片| 99精国产麻豆久久婷婷| 久久热在线av| 丝袜喷水一区| 日本av免费视频播放| 嫩草影院入口| 热99久久久久精品小说推荐| 精品久久久久久电影网| 日韩一区二区三区影片| 99国产综合亚洲精品| 午夜久久久在线观看| 亚洲成人一二三区av| 波野结衣二区三区在线| 久久ye,这里只有精品| 女的被弄到高潮叫床怎么办| 午夜福利网站1000一区二区三区| 两个人看的免费小视频| 18禁动态无遮挡网站| 久久久久久久久久人人人人人人| 免费久久久久久久精品成人欧美视频| 自线自在国产av| 亚洲第一区二区三区不卡| 欧美国产精品一级二级三级| 日日啪夜夜爽| 成人手机av| 欧美日韩国产mv在线观看视频| 一级爰片在线观看| 久久久久精品久久久久真实原创| 免费黄色在线免费观看| 久久精品人人爽人人爽视色| av一本久久久久| 三上悠亚av全集在线观看| 免费观看av网站的网址| 成人18禁高潮啪啪吃奶动态图| 九九爱精品视频在线观看| 永久免费av网站大全| 天美传媒精品一区二区| 男女无遮挡免费网站观看| 在线观看免费日韩欧美大片| 超碰成人久久| 看十八女毛片水多多多| 成人漫画全彩无遮挡| 亚洲伊人色综图| 精品人妻偷拍中文字幕| 天天躁日日躁夜夜躁夜夜| 久久毛片免费看一区二区三区| freevideosex欧美| 国产精品欧美亚洲77777| 性色avwww在线观看| 久久久久久人人人人人| 成人免费观看视频高清| 少妇人妻久久综合中文| 婷婷色综合www| 国产深夜福利视频在线观看| 久久国产精品大桥未久av| 一区在线观看完整版| 自拍欧美九色日韩亚洲蝌蚪91| xxx大片免费视频| 熟女少妇亚洲综合色aaa.| xxx大片免费视频| 精品第一国产精品| 欧美日韩国产mv在线观看视频| 亚洲精品成人av观看孕妇| 国产乱来视频区| 国产亚洲精品第一综合不卡| 亚洲精华国产精华液的使用体验| 亚洲久久久国产精品| 老司机影院成人| 亚洲精品第二区| 国产亚洲一区二区精品| 亚洲伊人色综图| 国产免费又黄又爽又色| 国产又爽黄色视频| 伦精品一区二区三区| 日韩,欧美,国产一区二区三区| 男的添女的下面高潮视频| 一边亲一边摸免费视频| 99久国产av精品国产电影| 亚洲第一av免费看| 天堂中文最新版在线下载| 日韩不卡一区二区三区视频在线| 久久99热这里只频精品6学生| 一级a爱视频在线免费观看| 日韩视频在线欧美| 亚洲国产精品一区二区三区在线| 爱豆传媒免费全集在线观看| 精品久久久精品久久久| 狠狠精品人妻久久久久久综合| 性高湖久久久久久久久免费观看| 国产精品秋霞免费鲁丝片| 亚洲国产欧美网| 91成人精品电影| 欧美xxⅹ黑人| 欧美精品国产亚洲| 免费黄色在线免费观看| 中文字幕人妻丝袜一区二区 | 中文字幕人妻熟女乱码| 人妻系列 视频| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一区二区在线不卡| 日韩电影二区| 99热网站在线观看| 18禁观看日本| 亚洲熟女精品中文字幕| 国产乱来视频区| 制服丝袜香蕉在线| av在线播放精品| 日韩欧美一区视频在线观看| 日韩大片免费观看网站| 九九爱精品视频在线观看| 亚洲一区二区三区欧美精品| 免费在线观看视频国产中文字幕亚洲 | 免费观看a级毛片全部| av女优亚洲男人天堂| 国产免费现黄频在线看| 色视频在线一区二区三区| 美女中出高潮动态图| 免费在线观看黄色视频的| 日韩 亚洲 欧美在线| 看十八女毛片水多多多| 观看美女的网站| 久久精品熟女亚洲av麻豆精品| 精品国产一区二区久久| 18禁动态无遮挡网站| 亚洲精品自拍成人| 亚洲国产精品成人久久小说| 最近中文字幕2019免费版| 午夜免费鲁丝| 亚洲男人天堂网一区| av免费在线看不卡| 精品99又大又爽又粗少妇毛片| 国产精品 国内视频| av电影中文网址| 在线天堂中文资源库| 亚洲精品国产色婷婷电影| 久久久a久久爽久久v久久| 看非洲黑人一级黄片| 99热全是精品| 1024视频免费在线观看| 亚洲一区中文字幕在线| 高清视频免费观看一区二区| 丰满乱子伦码专区| 秋霞伦理黄片| 高清av免费在线| 日韩一区二区三区影片| 日韩av不卡免费在线播放| 精品久久久精品久久久| 岛国毛片在线播放| 亚洲精品久久午夜乱码| 亚洲在久久综合| 亚洲内射少妇av| 午夜精品国产一区二区电影| 只有这里有精品99| 韩国av在线不卡| 永久免费av网站大全| 人妻人人澡人人爽人人| 国产成人精品在线电影| 久久韩国三级中文字幕| 女人久久www免费人成看片| 又粗又硬又长又爽又黄的视频| 久久久久久久大尺度免费视频| 久久久亚洲精品成人影院| 国产成人a∨麻豆精品| 高清视频免费观看一区二区| 久久综合国产亚洲精品| 久久精品aⅴ一区二区三区四区 | 成年av动漫网址| 国产精品国产三级专区第一集| 五月天丁香电影| 叶爱在线成人免费视频播放| 综合色丁香网| 欧美精品高潮呻吟av久久| 男女边摸边吃奶| 午夜免费鲁丝| 精品亚洲成a人片在线观看| 亚洲一码二码三码区别大吗| 国产人伦9x9x在线观看 | 国产精品一国产av| 美女福利国产在线| 精品人妻在线不人妻| 在线观看免费日韩欧美大片| 夜夜骑夜夜射夜夜干| 男人添女人高潮全过程视频| 美女主播在线视频| 午夜激情av网站| 亚洲成av片中文字幕在线观看 | 欧美 亚洲 国产 日韩一| 久久这里有精品视频免费| 国产熟女欧美一区二区| 免费观看av网站的网址| 另类精品久久| 色播在线永久视频| 一本大道久久a久久精品| 飞空精品影院首页| 久久av网站| 国产免费福利视频在线观看| 制服诱惑二区| 久久久精品94久久精品| 2021少妇久久久久久久久久久| 99久国产av精品国产电影| 狂野欧美激情性bbbbbb| 亚洲图色成人| 久久亚洲国产成人精品v| 在线看a的网站| 一级毛片 在线播放| 卡戴珊不雅视频在线播放| 久久久精品区二区三区| 久久久精品免费免费高清| 亚洲av日韩在线播放| 欧美精品国产亚洲| 99九九在线精品视频| www.自偷自拍.com| 欧美黄色片欧美黄色片| 1024视频免费在线观看| 18禁裸乳无遮挡动漫免费视频| 色视频在线一区二区三区| 欧美日韩精品网址| 精品视频人人做人人爽| 国产成人精品久久久久久| 午夜激情久久久久久久| 欧美 亚洲 国产 日韩一| 97在线人人人人妻| 夫妻性生交免费视频一级片| 一级毛片 在线播放| 国产精品熟女久久久久浪| 精品国产一区二区三区四区第35| 免费女性裸体啪啪无遮挡网站| 日日撸夜夜添| 99九九在线精品视频| av在线老鸭窝| 高清黄色对白视频在线免费看| 18禁裸乳无遮挡动漫免费视频| 美国免费a级毛片| 精品亚洲成a人片在线观看| 波野结衣二区三区在线| 极品人妻少妇av视频| 在线观看免费高清a一片| 午夜福利视频精品| 久久精品国产a三级三级三级| 国产男女超爽视频在线观看| 性少妇av在线| 春色校园在线视频观看| 欧美在线黄色| 女人久久www免费人成看片| 国产一区二区在线观看av| 亚洲综合精品二区| 国产在线一区二区三区精| 免费大片黄手机在线观看| 亚洲美女黄色视频免费看| 日本av手机在线免费观看| 国产成人欧美| 成人免费观看视频高清| 日本vs欧美在线观看视频| 亚洲情色 制服丝袜| 欧美亚洲日本最大视频资源| 国产一区二区三区综合在线观看| 成人国产麻豆网| 咕卡用的链子| 国产色婷婷99| av网站免费在线观看视频| 伦理电影免费视频| 久久人人爽人人片av| 人人妻人人爽人人添夜夜欢视频| 日本色播在线视频| 色视频在线一区二区三区| 久久精品国产a三级三级三级| 交换朋友夫妻互换小说| √禁漫天堂资源中文www| 久久精品夜色国产| 国产午夜精品一二区理论片| 2022亚洲国产成人精品| 亚洲国产精品一区三区| 香蕉精品网在线| 成人亚洲精品一区在线观看| 美女国产视频在线观看| 国产在视频线精品| 国产xxxxx性猛交| 久久久亚洲精品成人影院| 男女国产视频网站| 飞空精品影院首页| 亚洲人成77777在线视频| 免费不卡的大黄色大毛片视频在线观看| 国产精品熟女久久久久浪| 久久久欧美国产精品| 国产精品无大码| 好男人视频免费观看在线| 国产成人精品一,二区| 这个男人来自地球电影免费观看 | 亚洲情色 制服丝袜| 精品人妻一区二区三区麻豆| 满18在线观看网站| 国产综合精华液| 精品亚洲成a人片在线观看| 日韩三级伦理在线观看| 久久午夜综合久久蜜桃| 亚洲成人一二三区av| 国产熟女欧美一区二区| 久久久精品94久久精品| 国产一级毛片在线| 最近中文字幕2019免费版| 高清av免费在线| 在线观看www视频免费| 免费大片黄手机在线观看| 中文字幕色久视频| 啦啦啦啦在线视频资源| 国产精品av久久久久免费| 精品一区在线观看国产| 九九爱精品视频在线观看| 永久免费av网站大全| 久久精品久久精品一区二区三区| 国产又爽黄色视频| 亚洲天堂av无毛| 超碰成人久久| 免费看不卡的av| 一区二区av电影网| 一级a爱视频在线免费观看| 在线天堂最新版资源| 精品一区在线观看国产| 成人漫画全彩无遮挡| 亚洲国产精品一区三区| 精品久久久久久电影网| 亚洲国产精品一区二区三区在线| 欧美精品高潮呻吟av久久| 日日爽夜夜爽网站| 国产精品一国产av| 尾随美女入室| 一级毛片我不卡| 免费av中文字幕在线| 日韩一区二区三区影片| 深夜精品福利| 色吧在线观看| 亚洲av综合色区一区| 男女下面插进去视频免费观看| 国产精品女同一区二区软件| 男女午夜视频在线观看| 国产极品天堂在线| 9191精品国产免费久久| 亚洲男人天堂网一区| 爱豆传媒免费全集在线观看| 另类精品久久| 永久免费av网站大全| 中文字幕av电影在线播放| 久久精品熟女亚洲av麻豆精品| 午夜av观看不卡| 97人妻天天添夜夜摸| 精品国产一区二区三区久久久樱花| 黑人欧美特级aaaaaa片| 精品亚洲成a人片在线观看| 大香蕉久久成人网| 国产老妇伦熟女老妇高清| 久久久久国产网址| 精品久久蜜臀av无| 热re99久久精品国产66热6| 欧美xxⅹ黑人| 人人澡人人妻人| 成年动漫av网址| 日本wwww免费看| 国产精品免费视频内射| 老司机影院毛片| 国产综合精华液| 亚洲国产毛片av蜜桃av| 日韩不卡一区二区三区视频在线| 亚洲精品成人av观看孕妇| 国产男女超爽视频在线观看| 中国国产av一级| 天天影视国产精品| 人人妻人人澡人人看| 最近的中文字幕免费完整| 最近中文字幕2019免费版| 免费观看性生交大片5| 日韩一区二区视频免费看| 日本vs欧美在线观看视频| 久热久热在线精品观看| 一级片免费观看大全| 高清欧美精品videossex| 黄色一级大片看看| 多毛熟女@视频| 波多野结衣av一区二区av| 日韩中字成人| 性色av一级| 日韩av不卡免费在线播放| 久久精品国产综合久久久| 捣出白浆h1v1| av天堂久久9| 日韩欧美一区视频在线观看| 亚洲欧美成人综合另类久久久| 精品国产超薄肉色丝袜足j| 满18在线观看网站| 好男人视频免费观看在线| 国产在线视频一区二区| 亚洲精品国产色婷婷电影| 国产日韩欧美亚洲二区| 国产熟女午夜一区二区三区| 老司机影院成人| 午夜福利在线免费观看网站| 日本爱情动作片www.在线观看| 欧美变态另类bdsm刘玥| 国产成人精品婷婷| 精品人妻一区二区三区麻豆| 精品国产乱码久久久久久小说| 少妇人妻精品综合一区二区| 亚洲第一区二区三区不卡| 亚洲成人手机| 亚洲第一青青草原| 男女边吃奶边做爰视频| 精品99又大又爽又粗少妇毛片| 精品国产露脸久久av麻豆| 国产精品女同一区二区软件| 免费黄频网站在线观看国产| 欧美精品一区二区大全| 丝袜美足系列| 久久鲁丝午夜福利片| 在线观看国产h片| 欧美变态另类bdsm刘玥| 国产男女内射视频| 亚洲国产毛片av蜜桃av| 国产成人91sexporn| 久久久久久久亚洲中文字幕| 精品午夜福利在线看| 欧美97在线视频| 久久久国产一区二区| 成年动漫av网址| 男女国产视频网站| 亚洲国产精品999| 黄片播放在线免费| 精品少妇黑人巨大在线播放| 久久青草综合色| 伦理电影大哥的女人| 午夜福利视频精品| 99精国产麻豆久久婷婷| 国产精品熟女久久久久浪| 亚洲av欧美aⅴ国产| 日韩一区二区三区影片| 91aial.com中文字幕在线观看| 亚洲精品在线美女| 美女主播在线视频| 自线自在国产av| av福利片在线| 黑人巨大精品欧美一区二区蜜桃| 十分钟在线观看高清视频www| 国产成人精品一,二区| 免费av中文字幕在线| 欧美+日韩+精品| 午夜福利视频精品| 国产精品av久久久久免费| 老司机亚洲免费影院| 这个男人来自地球电影免费观看 | 乱人伦中国视频| 99久久中文字幕三级久久日本| 欧美黄色片欧美黄色片| 男人爽女人下面视频在线观看| 中文字幕人妻熟女乱码| 80岁老熟妇乱子伦牲交| 最近手机中文字幕大全| 波多野结衣一区麻豆| 久久久精品免费免费高清| 一区二区三区乱码不卡18| 五月伊人婷婷丁香| 一区二区日韩欧美中文字幕| av电影中文网址| 亚洲熟女精品中文字幕| 国产av一区二区精品久久| 又黄又粗又硬又大视频| 久久久久精品久久久久真实原创| 国产午夜精品一二区理论片| 久久国内精品自在自线图片| 日韩中字成人| 亚洲精品一二三| a级毛片黄视频| 国产精品秋霞免费鲁丝片| 熟女av电影|