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

    雙軸受壓狀態(tài)下的高延性纖維增強(qiáng)水泥基復(fù)合材料本構(gòu)模型

    2020-12-18 03:08:02陳善富陳靜芬楊鳳祥劉志明
    工程力學(xué) 2020年12期
    關(guān)鍵詞:主壓雙軸單軸

    陳善富,陳靜芬,楊鳳祥,劉志明

    (暨南大學(xué)力學(xué)與建筑工程學(xué)院“重大工程災(zāi)害與控制”教育部重點(diǎn)實(shí)驗(yàn)室,廣州510632)

    高延性纖維增強(qiáng)水泥基復(fù)合材料(Engineered Cementitious Composites,簡稱ECC)是一種高性能水泥基復(fù)合材料。ECC纖維摻量低且在受拉和受壓時具有高延性的特點(diǎn),同時因其良好的準(zhǔn)應(yīng)變硬化特性和多縫開裂特性,被廣泛應(yīng)用于橋面板、橋面連接板、建筑防震抗震構(gòu)件、混凝土保護(hù)層等實(shí)際工程中[1]。

    1992年,密歇根大學(xué)的Li等[2]根據(jù)細(xì)觀力學(xué)和斷裂力學(xué)的理論,在水泥基體中添加短纖維研制出具有高延性的水泥基復(fù)合材料ECC。在我國,眾多研究ECC的學(xué)者將其命名為高延性纖維增強(qiáng)水泥基復(fù)合材料[3?7]。

    自ECC面世以來,國內(nèi)外學(xué)者對其力學(xué)性能進(jìn)行了大量的試驗(yàn)研究[8?17],主要集中于單軸拉伸、單軸壓縮、彎曲和剪切等方面。試驗(yàn)研究[8?17]表明,ECC因纖維橋聯(lián)作用而具有良好的拉伸,壓縮和彎曲韌性,同時具有良好的裂縫控制能力和耐久性。

    有限元分析是對工程結(jié)構(gòu)進(jìn)行設(shè)計和研究的重要手段,要對一種材料的工程結(jié)構(gòu)進(jìn)行有限元分析需要該材料的本構(gòu)模型。但是,目前的大型通用有限元程序(如ANSYS,ABAQUS等)的材料庫中,尚未擁有ECC的本構(gòu)模型。在對ECC 結(jié)構(gòu)構(gòu)件進(jìn)行有限元分析時,一些學(xué)者仍然采用有限元程序材料庫中現(xiàn)有的混凝土本構(gòu)模型[5,18?19]。因此,開發(fā)一種能夠較準(zhǔn)確地體現(xiàn)ECC力學(xué)行為的本構(gòu)模型對推動ECC在工程結(jié)構(gòu)中的廣泛應(yīng)用具有重要意義。

    為了對ECC結(jié)構(gòu)構(gòu)件進(jìn)行有限元分析,許多學(xué)者提出了能夠用于預(yù)測ECC 力學(xué)行為的本構(gòu)模型。Han 等[20]提出了一種基于總應(yīng)變的正交各向異性二維旋轉(zhuǎn)裂縫模型模擬ECC在循環(huán)荷載下的力學(xué)響應(yīng),該模型分別定義了受拉和受壓時的滯回行為,且被吸收進(jìn)OpenSees地震工程模擬軟件的材料庫中,便于對ECC構(gòu)件在承受地震作用時的力學(xué)行為進(jìn)行模擬。但是,二維模型存在受形狀、加載方式和邊界條件限制的缺點(diǎn)。為了克服上述缺點(diǎn),Hung等[21]在Han 模型[20]的基礎(chǔ)上發(fā)展了一種ECC的三維本構(gòu)模型,并編寫了本構(gòu)模型的用戶自定義材料子程序UMAT嵌入LS-DYNA中描述ECC的力學(xué)行為,同時建立ECC沖切板與兩跨連續(xù)梁的三維有限元模型進(jìn)行數(shù)值計算,數(shù)值計算結(jié)果與試驗(yàn)結(jié)果較吻合。張曉悅等[22]基于Rots[23]旋轉(zhuǎn)裂縫理論,提出了一種ECC的三維本構(gòu)模型,該模型通過斷裂能來判斷裂縫平面的旋轉(zhuǎn),通過模擬ECC單軸受拉和四點(diǎn)彎曲試驗(yàn),得到了與試驗(yàn)較吻合的結(jié)果。

    上述三種模型都是基于彌散裂縫理論,假定ECC 由于在加載方向開裂后的材料性質(zhì)與未加載方向的材料性質(zhì)不同而轉(zhuǎn)化為正交各向異性的材料,較好地模擬出ECC在多軸受力狀態(tài)下的力學(xué)響應(yīng)。但由于目前對ECC的壓縮試驗(yàn)研究還是主要集中于單軸壓縮試驗(yàn),上述模型判定材料受壓破壞的強(qiáng)度參數(shù)也僅來源于從ECC單軸壓縮試驗(yàn)得到的單軸應(yīng)力-應(yīng)變曲線中:Han 等[20]根據(jù)Kesner等[24]的ECC單軸壓縮試驗(yàn)結(jié)果采用雙線性的應(yīng)力-應(yīng)變關(guān)系描述ECC的單軸受壓行為;Hung 等[21]根據(jù)ECC單軸壓縮試驗(yàn)結(jié)果[25]采用三線性的應(yīng)力-應(yīng)變關(guān)系描述ECC的單軸受壓行為;張曉悅等[22]用一條非線性的Hognestad 拋物曲線[26]來描述ECC的單軸壓縮應(yīng)力-應(yīng)變關(guān)系的上升段,其強(qiáng)度參數(shù)取自Yoo等[27]的ECC單軸壓縮試驗(yàn)。

    目前,對于ECC受壓行為的試驗(yàn)研究主要集中于單軸壓縮試驗(yàn)。除了上述三種模型所引用的ECC 單軸壓縮試驗(yàn)外,一些學(xué)者對不同配合比的ECC 進(jìn)行了單軸壓縮試驗(yàn),得到了不同配合比的ECC 單軸壓縮應(yīng)力應(yīng)變曲線,并采用不同的曲線方程進(jìn)行擬合:徐世烺等[28]采用二次拋物線形式描述其上升段,采用雙折線形式描述其下降段;姜海軍等[29]采用過鎮(zhèn)海[30]提出的混凝土本構(gòu)方程描述其上升段,采用Sargin[31]提出的混凝土本構(gòu)方程描述其下降段;李艷等[11]用兩個非線性的方程分別擬合了其上升段和下降段;Zhou 等[15]用線性和非線性組合的分段形式描述其上升段,用雙折線形式描述其下降段。

    上述三種模型中判定ECC 受壓破壞的強(qiáng)度參數(shù)皆取自ECC單軸壓縮試驗(yàn),認(rèn)為ECC在某一方向達(dá)到單軸抗壓強(qiáng)度則該方向發(fā)生受壓破壞。然而,現(xiàn)有的ECC雙軸試驗(yàn)結(jié)果都表明[32?34],ECC在承受雙軸壓應(yīng)力時,一個方向上的抗壓強(qiáng)度是受其垂直方向上壓應(yīng)力的大小影響的。這是由于ECC 在雙軸受壓時,一個方向上的壓應(yīng)力會對其垂直方向上的壓應(yīng)力所引起的側(cè)向變形產(chǎn)生一定程度的約束,限制了內(nèi)部微裂縫的發(fā)展,從而提高了其垂直方向的抗壓強(qiáng)度,導(dǎo)致其垂直方向的應(yīng)力應(yīng)變關(guān)系與單軸受壓時不同。另外,文獻(xiàn)[33]研究了ECC 在單軸受壓和雙軸受壓時的應(yīng)力-應(yīng)變曲線的非線性特性,并擬合了ECC在不同雙軸壓應(yīng)力比下的應(yīng)力-應(yīng)變?nèi)€方程。因此,若在本構(gòu)模型中不考慮ECC在雙軸受壓時的非線性應(yīng)力-應(yīng)變關(guān)系和抗壓強(qiáng)度變化,就不能準(zhǔn)確地預(yù)測ECC在承受不同雙軸壓應(yīng)力時的力學(xué)行為。

    要在本構(gòu)模型中同時考慮ECC在雙軸受壓狀態(tài)下因應(yīng)力比的不同而產(chǎn)生的非線性應(yīng)力-應(yīng)變關(guān)系和抗壓強(qiáng)度的變化,顯然各向同性本構(gòu)模型無法展現(xiàn)上述行為。Darwin 和Pecknold[35?36]為了展現(xiàn)混凝土在雙軸受壓時的非線性應(yīng)力-應(yīng)變關(guān)系和抗壓強(qiáng)度變化,提出了一種混凝土的二維非線性正交各向異性本構(gòu)模型。該模型假定混凝土在相互垂直的兩個加載方向上為正交各向異性材料,并采用Kupfer 雙軸強(qiáng)度包絡(luò)線[37]來預(yù)測不同雙軸應(yīng)力比下混凝土兩個正交各向異性方向上的抗壓強(qiáng)度的變化,通過引入等效單軸應(yīng)變的概念將混凝土單軸的應(yīng)力-應(yīng)變關(guān)系擴(kuò)展到正交各向異性兩個方向上,能夠同時考慮混凝土在雙軸應(yīng)力狀態(tài)下的強(qiáng)度變化和非線性應(yīng)力-應(yīng)變關(guān)系,從而能準(zhǔn)確地描述混凝土在雙軸應(yīng)力狀態(tài)下的力學(xué)行為。文獻(xiàn)[23]的試驗(yàn)研究表明,在雙軸受壓狀態(tài)下ECC 兩個加載方向上的非線性應(yīng)力-應(yīng)變關(guān)系和抗壓強(qiáng)度是隨著雙軸壓應(yīng)力比的變化而變化的,且變化規(guī)律與混凝土相似。所以,為了在本構(gòu)模型中展現(xiàn)ECC上述力學(xué)行為,可以建立ECC在兩個相互垂直的加載方向上的正交各向異性本構(gòu)關(guān)系,然后通過引入雙軸強(qiáng)度準(zhǔn)則的方式考慮ECC兩個加載方向上的強(qiáng)度變化。

    因此,本文的主要目的是開發(fā)一種能夠同時體現(xiàn)ECC 在雙軸受壓狀態(tài)下的非線性應(yīng)力-應(yīng)變關(guān)系和抗壓強(qiáng)度變化的本構(gòu)模型,嵌入有限元分析程序ABAQUSv6.14中,用于預(yù)測ECC在雙軸受壓狀態(tài)下的非線性應(yīng)力-應(yīng)變關(guān)系和抗壓強(qiáng)度變化。本文建立了本構(gòu)模型的控制方程,推導(dǎo)了本構(gòu)模型的顯式數(shù)值算法,采用Fortran 語言編寫該算法的用戶自定義材料子程序UMAT,并嵌入大型通用有限元分析軟件ABAQUSv6.14中,通過對ECC試件的雙軸受壓試驗(yàn)進(jìn)行有限元分析驗(yàn)證了所開發(fā)的本構(gòu)模型的有效性。

    1 雙軸受壓狀態(tài)下的ECC本構(gòu)模型

    本文提出的雙軸受壓狀態(tài)下的ECC本構(gòu)模型能夠描述ECC在雙軸受壓時的非線性力學(xué)行為和抗壓強(qiáng)度變化的特點(diǎn)。模型包含:1)ECC在兩個相互垂直的加載方向的二維增量型正交各向異性本構(gòu)關(guān)系,適合于在有限元程序中采用增量-迭代法求解非線性問題[38];2)一條非線性的受壓應(yīng)力-應(yīng)變曲線用于描述ECC雙軸受壓時兩個加載方向上展現(xiàn)的非線性力學(xué)行為;3)一條雙軸強(qiáng)度包絡(luò)線用于描述ECC 雙軸受壓時兩個加載方向上的抗壓強(qiáng)度變化。

    1.1 ECC二維增量型本構(gòu)關(guān)系

    針對ECC雙軸受壓時展現(xiàn)的非線性特性,在每一增量步內(nèi)建立ECC 在兩個垂直加載方向上的二維增量型正交各向異性本構(gòu)關(guān)系為:

    為了將單軸應(yīng)力-應(yīng)變曲線擴(kuò)展為雙軸應(yīng)力狀態(tài)下的應(yīng)力-應(yīng)變曲線,這里采用文獻(xiàn)[35]的方法,引入等效單軸應(yīng)變的概念。等效單軸應(yīng)變是由實(shí)際應(yīng)變消除泊松比效應(yīng)得到的假定值,沒有實(shí)際的物理意義。但有必要指出的是,本研究通過引入等效單軸應(yīng)變有利于采用應(yīng)力-等效單軸應(yīng)變的形式來分別描述ECC 在1、2方向上的受壓應(yīng)力-應(yīng)變關(guān)系,并通過該關(guān)系式求解任意應(yīng)力狀態(tài)下1、2方向的切線模量。

    式中,k為增量步總數(shù)。

    式(9)表明在模型的數(shù)值算法中,1、2方向上等效單軸總應(yīng)變 εui是加載過程中每一增量步等效單軸應(yīng)變增量dεui的總和。

    1.2 ECC受壓應(yīng)力-應(yīng)變曲線

    圖1 ECC 壓縮應(yīng)力-應(yīng)變曲線示意圖Fig.1 Illustration of compressive stress-strain curve of ECC

    本文將文獻(xiàn)[33]擬合的ECC受壓應(yīng)力-應(yīng)變多項式曲線方程擴(kuò)展到引入等效單軸應(yīng)變后各方向上的應(yīng)力-等效單軸應(yīng)變關(guān)系,如圖2所示,表達(dá)式為:

    通過將式(11)的 σi對 εui求導(dǎo)可求得式(2)中1、2方向的切線模量Ei(i=1,2):

    圖2 ECC應(yīng)力-等效單軸應(yīng)變曲線示意圖Fig.2 Illustration of stress-equivalent uniaxial strain curve of ECC

    1.3 ECC雙軸強(qiáng)度包絡(luò)線

    圖3 ECC 雙軸強(qiáng)度包絡(luò)線示意圖Fig.3 Illustration of biaxial strength envelope of ECC

    為了展示在本文提出的雙軸受壓狀態(tài)下的ECC本構(gòu)模型(下文標(biāo)記為“模型1”)中考慮雙軸受壓狀態(tài)下ECC抗壓強(qiáng)度變化的必要性,本文還建立了不考慮雙軸受壓狀態(tài)下ECC抗壓強(qiáng)度變化的本構(gòu)模型(下文標(biāo)記為“模型2”),即模型2沒有引入第1.3節(jié)中所述的理論,然后將模型1和模型2應(yīng)用于不同配合比的ECC試件在不同應(yīng)力比下的雙軸受壓有限元分析中,并對采用模型1和模型2進(jìn)行數(shù)值計算得到的結(jié)果進(jìn)行對比分析。

    2 數(shù)值實(shí)現(xiàn)方法

    本文分別推導(dǎo)了模型1和模型2的數(shù)值算法,并采用Fortran 語言編寫了兩個模型的用戶自定義材料子程序UMAT以更新應(yīng)力、1和2方向的彈性模量、剪切模量以及相關(guān)狀態(tài)變量,將子程序嵌入ABAQUSv6.14中對ECC 雙軸受壓試驗(yàn)進(jìn)行有限元分析。

    模型1的數(shù)值算法流程如下所示(上標(biāo)t表示當(dāng)前增量步,t?1表示上一增量步):

    a)初始條件:

    3 模型驗(yàn)證

    本文通過對文獻(xiàn)[33]中兩種配合比的ECC試件在不同應(yīng)力比下的雙軸加載試驗(yàn)進(jìn)行有限元分析驗(yàn)證模型1的有效性。

    兩種配合比的ECC 試件均由普通硅酸鹽水泥、粉煤灰、精制石英砂、水、高效減水劑和PVA 纖維制作而成。ECC-I 試件的配合比為水泥∶粉煤灰∶石英砂∶水∶高效減水劑=1∶2∶0.6∶0.66∶0.024,纖維體積率為2.0%;ECC-II試件的配合比為水泥∶粉煤灰∶石英砂∶水∶高效減水劑=1∶3∶0.8∶0.92∶0.024,纖維體積率為2.0%。

    通過單軸和雙軸受壓試驗(yàn)獲得的材料屬性見表1。本文假定ECC-I和ECC-II試件的等效泊松比為ν=0.2;模型參數(shù)l、m和n取自文獻(xiàn)[35],分別為?1.6、2.25和0.35。試件尺寸為100 mm×100 mm×100 mm,分別進(jìn)行了應(yīng)力比為0、0.25、0.5、0.75 和1時的雙軸受壓試驗(yàn)。

    本文利用有限元程序ABAQUSv6.14進(jìn)行數(shù)值分析,建立了圖4所示的有限元模型分別進(jìn)行ECC-I和ECC-II試件在上述五種應(yīng)力比下的雙軸受壓分析。模型尺寸為100 mm×100 mm。在模型的左側(cè)和下側(cè)分別施加1、2方向的約束,在模型的右側(cè)和上側(cè)分別施加均布應(yīng)力 σ1和 σ2,并令σ2/σ1=0、0.25、0.5、0.75和1以模擬ECC試件在五種應(yīng)力比下的雙軸受壓試驗(yàn)。為了保證在各應(yīng)力比下均能達(dá)到抗壓強(qiáng)度,本文在ECC-I試件和ECC-II試件的1方向施加的初始荷載分別為75 MPa 和50 MPa,根據(jù)應(yīng)力比可得2 方向上施加的初始荷載。如圖4所示,為了考察不同網(wǎng)格劃分方案對模型計算結(jié)果的影響,以ECC-I試件在單軸受壓(β=0)和等壓雙軸受壓(β=1)時為例,分別采用兩種網(wǎng)格劃分方案進(jìn)行計算:圖4(a)1×1網(wǎng)格和圖4(b)5×5網(wǎng)格。有限單元類型采用線性減縮積分平面應(yīng)力單元(CPS4R)。

    表1 材料屬性[33]Table 1 Material properties

    圖4 ECC試件雙軸受壓有限元模型Fig.4 Finite element model of ECCspecimens under biaxial compression

    如圖5所示為ECC-I試件在應(yīng)力比β=0和β=1時分別采用1×1和5×5網(wǎng)格劃分方案的主壓應(yīng)力方向的應(yīng)力-應(yīng)變曲線計算結(jié)果,可以看出在應(yīng)力比β=0和β=1下,兩種網(wǎng)格劃分方案計算得到的主壓應(yīng)力方向的應(yīng)力-應(yīng)變曲線重疊在一起,并且在β=0時兩種網(wǎng)格劃分方案計算得到的極限抗壓強(qiáng)度和極限抗壓應(yīng)變均一致,分別為?57.657 MPa和?0.298%;在β=1時也一致,分別為?66.527 MPa和?0.327%。由此可知,采用圖4所示的不同網(wǎng)格劃分方案并不影響數(shù)值計算結(jié)果。這是由于圖4中的模型在兩個方向上施加均布應(yīng)力(σ1,σ2)的情況下,當(dāng)采用1×1網(wǎng)格劃分方案,即只有1個單元時,顯然該單元的外力為(σ1,σ2);當(dāng)采用5×5網(wǎng)格劃分方案時,根據(jù)靜力平衡條件,每個單元的外力同樣為(σ1,σ2),因此可以得到相同的計算結(jié)果。為了節(jié)省數(shù)值計算時間,后文算例統(tǒng)一采用1×1的網(wǎng)格劃分方案。

    圖5 應(yīng)力比β=0和β=1的ECC-I 試件采用1×1和5×5網(wǎng)格劃分方案的主壓應(yīng)力方向的預(yù)測應(yīng)力-應(yīng)變曲線Fig.5 Predicted stress-strain curvesof ECC-I specimens with stress ratios β=0 and β=1 in the major compressive stress direction with 1×1 and 5×5 meshing scheme

    如圖6所示為ECC-I試件在應(yīng)力比β=0和β=1時分別采用三種最大時間增量步長計算得到的主壓應(yīng)力方向的預(yù)測應(yīng)力應(yīng)變曲線,表明所采用的三種最大時間增量步長均獲得高度一致的計算結(jié)果。表2所示為ECC-I試件在應(yīng)力比β=0和β=1時分別采用三種最大時間增量步長時的極限抗壓強(qiáng)度和極限抗壓應(yīng)變計算結(jié)果,可知隨著增量步長的減小,計算得到的極限抗壓強(qiáng)度和極限抗壓應(yīng)變越接近于精確解。在應(yīng)力比β=0和β=1時,采用2×10?4s 進(jìn)行計算時得到的極限抗壓強(qiáng)度和極限抗壓應(yīng)變與精確解的誤差分別小于0.1%和1%。為了保證其余算例與精確解的誤差也能在可接受的范圍內(nèi),本文將最大時間增量步長統(tǒng)一設(shè)為1×10?4s。

    圖6 應(yīng)力比β=0和β=1的ECC-I 試件在不同最大時間增量步長下的主壓應(yīng)力方向的預(yù)測應(yīng)力-應(yīng)變曲線Fig.6 Predicted stress-strain curvesof ECC-I specimens with stress ratios β=0 and β=1 in the major compressive stress direction under variousmaximum time increments

    表2 應(yīng)力比β=0和β=1的ECC-I 試件在不同最大時間增量步長下的主壓應(yīng)力方向預(yù)測極限抗壓強(qiáng)度和極限抗壓應(yīng)變Table 2 Predicted ultimate compressive strengthsand ultimate compressive strains of ECC-I specimens with stress ratios β=0 and β=1 in the major compressive stress direction under variousmaximum time increments

    圖7和圖8分別為ECC-I和ECC-II 試件在雙軸受壓時不同應(yīng)力比β(β=σ2/σ1, σ1<σ2<0)下主壓應(yīng)力方向,即1方向的σ1?ε1曲線的試驗(yàn)和數(shù)值計算結(jié)果。

    由圖7和圖8可知,在單軸受壓(β=0)時,采用模型1和模型2分別進(jìn)行有限元分析得到的兩種ECC試件的σ1?ε1曲線與試驗(yàn)曲線吻合良好,且模型1和模型2的計算曲線重合。這是由于當(dāng)試件處于單軸受壓狀態(tài),模型1用雙軸強(qiáng)度準(zhǔn)則可確定σcp1=fc′p和εcp1=ε′cp,而模型2中令σcp1=fc′p和εcp1=ε′cp,因此每一增量步中模型1和模型2的E1相等,從而計算得到了重合的σ1?ε1曲線。

    對于雙軸受壓(β≠0)的情況,采用模型1計算得到的σ1?ε1曲線與試驗(yàn)結(jié)果吻合良好,而采用模型2時則吻合程度較差。試驗(yàn)曲線和采用模型1進(jìn)行有限元分析得到的σ1?ε1曲線都表現(xiàn)為在加載開始階段,曲線的斜率變化幅度較小,近似為直線,隨后越接近破壞點(diǎn)曲線的斜率下降越明顯,到達(dá)破壞時曲線斜率為零。表明模型1 能夠較準(zhǔn)確地預(yù)測各雙軸壓應(yīng)力比下ECC試件在主壓應(yīng)力方向加載過程中的剛度變化特征。而模型2在各雙軸壓應(yīng)力比下的σ1?ε1曲線的初始斜率均相同,說明模型2無法準(zhǔn)確預(yù)測ECC試件在不同雙軸壓應(yīng)力狀態(tài)下的剛度變化趨勢。

    圖7 各雙軸壓應(yīng)力比下ECC-I 試件主壓應(yīng)力方向的應(yīng)力-應(yīng)變曲線Fig.7 Stress-strain curvesof ECC-I specimensin the major compressive stress direction under different biaxial stress ratios

    表3為ECC-I和ECC-II試件在各應(yīng)力比下主壓應(yīng)力方向的試驗(yàn)和預(yù)測極限抗壓強(qiáng)度和極限抗壓應(yīng)變值。由表3可知模型1預(yù)測得到的ECC-I試件的極限抗壓強(qiáng)度值與試驗(yàn)值相比,在β=1時最小,僅為0.03%,當(dāng)β=0.25時的誤差亦小于5%,模型1預(yù)測得到的ECC-II試件的極限抗壓強(qiáng)度值與試驗(yàn)值相比,在β=1時誤差最小為?0.00%,在β=0.5時最大,僅為?6.17%。表明模型1能夠準(zhǔn)確地預(yù)測雙軸受壓時ECC主壓應(yīng)力方向的抗壓強(qiáng)度。

    模型1預(yù)測得到的ECC-I和ECC-II試件在各應(yīng)力比下的極限抗壓應(yīng)變也同樣與試驗(yàn)結(jié)果較吻合,但預(yù)測到得到的ECC-I試件在β=0.25和0.5時和ECC-II試件在β=0.25時的誤差分別達(dá)到了16.07%,15.34%和?14.89%,這是由于試驗(yàn)數(shù)據(jù)離散性較大而式(17)假定抗壓強(qiáng)度σcpi與等效單軸抗壓應(yīng)變εcpi為線性關(guān)系,在上述三種情況下對應(yīng)的抗壓強(qiáng)度與等效單軸應(yīng)變的坐標(biāo)點(diǎn)剛好位于線性關(guān)系直線的兩側(cè),故僅此三組有限元預(yù)測值與試驗(yàn)結(jié)果相差較大。

    圖8 各雙軸壓應(yīng)力比下ECC-II 試件主壓應(yīng)力方向的應(yīng)力-應(yīng)變曲線Fig.8 Stress-strain curves of ECC-II specimens in the major compressive stress direction under different biaxial stress ratios

    由表3可知模型2預(yù)測得到的ECC-I和ECCII試件在各應(yīng)力比下主壓應(yīng)力方向極限抗壓強(qiáng)度和極限抗壓應(yīng)變除了在單軸受壓(β=0)時較準(zhǔn)確,在雙軸受壓時皆與試驗(yàn)值相比相差巨大。表明模型2無法預(yù)測ECC抗壓強(qiáng)度在雙軸壓應(yīng)力狀態(tài)下的變化。

    從數(shù)值分析算例可知,模型1同時考慮了ECC在雙軸受壓時的非線性應(yīng)力-應(yīng)變關(guān)系和應(yīng)力比的變化對抗壓強(qiáng)度的影響,能更準(zhǔn)確地描述不同配合比的ECC在不同應(yīng)力比雙軸受壓狀態(tài)下的力學(xué)行為。

    表3 各雙軸壓應(yīng)力比下ECC-I 和ECC-II 試件主壓應(yīng)力方向的試驗(yàn)和預(yù)測極限抗壓強(qiáng)度和極限抗壓應(yīng)變值Table 3 Experimental and predicted ultimatecompressive strengthsand ultimatecompressive strains of ECC-I and ECC-II specimensin the major compressivestress direction under various biaxial stress ratios

    4 結(jié)論

    (1)建立了一個能夠同時考慮ECC在雙軸受壓時的非線性力學(xué)行為和抗壓強(qiáng)度變化的二維本構(gòu)模型,開發(fā)了本構(gòu)模型的顯式積分?jǐn)?shù)值算法,采用Fortran 語言將模型算法編寫成用戶自定義材料子程序UMAT,并嵌入有限元分析軟件ABAQUS v6.14中,用于預(yù)測ECC在雙軸受壓時的非線性力學(xué)行為及抗壓強(qiáng)度變化。

    (2)由數(shù)值計算結(jié)果可知,本文模型數(shù)值計算得到的兩種不同配合比的ECC 試件在五組雙軸壓應(yīng)力比下的應(yīng)力-應(yīng)變曲線與試驗(yàn)曲線吻合良好,并且能夠較準(zhǔn)確地預(yù)測各應(yīng)力比下兩種不同配合比的ECC試件在主壓應(yīng)力方向上的抗壓強(qiáng)度和抗壓應(yīng)變,表明本文開發(fā)的ECC 本構(gòu)模型能夠較準(zhǔn)確地預(yù)測不同配合比ECC試件在雙軸受壓狀態(tài)下的力學(xué)行為。

    (3)作為對比,建立了不考慮ECC在雙軸受壓狀態(tài)下抗壓強(qiáng)度變化的本構(gòu)模型,并對上述兩種不同配合比的ECC試件在五組應(yīng)力比下進(jìn)行數(shù)值分析,結(jié)果表明為了準(zhǔn)確描述ECC雙軸受壓狀態(tài)下的力學(xué)行為,有必要考慮其雙軸受壓狀態(tài)下的抗壓強(qiáng)度變化。本文提出的本構(gòu)模型為準(zhǔn)確預(yù)測ECC 在雙軸受壓狀態(tài)下的非線性應(yīng)力-應(yīng)變關(guān)系和破壞強(qiáng)度提供了一種有效的分析方法。

    猜你喜歡
    主壓雙軸單軸
    喜馬拉雅晚期構(gòu)造應(yīng)力場及其與油氣分布的關(guān)系
    ——以準(zhǔn)噶爾盆地南緣為例
    斷塊油氣田(2021年4期)2021-07-28 07:56:42
    單軸壓縮條件下巖石峰后第Ⅱ種類型應(yīng)力——應(yīng)變曲線的新解釋
    高密度棉花打包機(jī)節(jié)能設(shè)計的分析
    CFRP-鋼復(fù)合板的單軸拉伸力學(xué)性能
    單軸應(yīng)變Si NMOS電流模型研究
    電子測試(2017年12期)2017-12-18 06:35:42
    簡易雙軸立銑頭裝置設(shè)計與應(yīng)用
    基于SolidWorks對雙軸攪拌機(jī)的靜力學(xué)分析
    自動打包機(jī)液壓系統(tǒng)的設(shè)計
    斜單軸跟蹤式光伏組件的安裝傾角優(yōu)化設(shè)計
    雙軸太陽能跟蹤與市電互補(bǔ)的路燈控制系統(tǒng)
    又爽又黄无遮挡网站| 最近中文字幕2019免费版| 真实男女啪啪啪动态图| 午夜免费观看性视频| 亚洲av成人精品一区久久| 久久精品熟女亚洲av麻豆精品| 99热这里只有是精品50| 亚洲国产精品999| 欧美+日韩+精品| 国产成人午夜福利电影在线观看| 五月开心婷婷网| 久久久久国产网址| 黄色一级大片看看| 国产av码专区亚洲av| av在线蜜桃| 亚洲精品乱码久久久久久按摩| 只有这里有精品99| 亚洲av男天堂| 国产精品爽爽va在线观看网站| 亚洲,一卡二卡三卡| 18禁裸乳无遮挡免费网站照片| 男人爽女人下面视频在线观看| 久久久久性生活片| av线在线观看网站| 赤兔流量卡办理| 国产成人精品福利久久| 免费观看在线日韩| 99re6热这里在线精品视频| 蜜桃久久精品国产亚洲av| 精品久久国产蜜桃| 国产av国产精品国产| 国产免费一区二区三区四区乱码| 麻豆乱淫一区二区| 亚洲av在线观看美女高潮| 亚洲丝袜综合中文字幕| 岛国毛片在线播放| 久久久久久九九精品二区国产| 精华霜和精华液先用哪个| 亚洲成人中文字幕在线播放| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产色婷婷电影| 午夜激情久久久久久久| 国产欧美日韩一区二区三区在线 | 性插视频无遮挡在线免费观看| 免费观看a级毛片全部| 国产有黄有色有爽视频| 欧美97在线视频| 国产69精品久久久久777片| 最近手机中文字幕大全| 99热网站在线观看| 精品人妻偷拍中文字幕| av线在线观看网站| 国产老妇伦熟女老妇高清| av播播在线观看一区| 国产男女超爽视频在线观看| 精品视频人人做人人爽| 欧美 日韩 精品 国产| 久久久久九九精品影院| 我的老师免费观看完整版| 小蜜桃在线观看免费完整版高清| 亚洲国产成人一精品久久久| 18禁裸乳无遮挡免费网站照片| 国产精品偷伦视频观看了| 亚洲成人av在线免费| 在线精品无人区一区二区三 | 别揉我奶头 嗯啊视频| 夜夜爽夜夜爽视频| 亚洲真实伦在线观看| 男女国产视频网站| 又爽又黄无遮挡网站| 伊人久久精品亚洲午夜| 26uuu在线亚洲综合色| 视频中文字幕在线观看| 亚洲精品乱码久久久v下载方式| 欧美人与善性xxx| 国产精品无大码| 亚洲av电影在线观看一区二区三区 | 亚洲va在线va天堂va国产| 中文字幕免费在线视频6| 少妇高潮的动态图| 97超视频在线观看视频| 久久人人爽人人片av| 免费av观看视频| 亚洲国产日韩一区二区| 99re6热这里在线精品视频| 成人毛片a级毛片在线播放| 一级片'在线观看视频| 麻豆成人午夜福利视频| 久久久欧美国产精品| 日日啪夜夜撸| 最近中文字幕高清免费大全6| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 免费不卡的大黄色大毛片视频在线观看| 高清欧美精品videossex| 亚洲三级黄色毛片| 国产精品成人在线| 性插视频无遮挡在线免费观看| 国产日韩欧美亚洲二区| av在线蜜桃| 成人特级av手机在线观看| 一级黄片播放器| 国产毛片在线视频| 九九久久精品国产亚洲av麻豆| av国产久精品久网站免费入址| 亚洲激情五月婷婷啪啪| 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区三区| 秋霞在线观看毛片| 五月开心婷婷网| 99热网站在线观看| 日韩亚洲欧美综合| 国产成人一区二区在线| 最新中文字幕久久久久| 美女cb高潮喷水在线观看| 午夜精品一区二区三区免费看| 午夜福利网站1000一区二区三区| 欧美一区二区亚洲| 亚洲不卡免费看| 精品国产一区二区三区久久久樱花 | 街头女战士在线观看网站| 校园人妻丝袜中文字幕| 国产高清三级在线| 汤姆久久久久久久影院中文字幕| 久久久久久久精品精品| 丝袜喷水一区| 国产午夜精品久久久久久一区二区三区| 欧美zozozo另类| 国产在线一区二区三区精| 成年女人看的毛片在线观看| 嫩草影院精品99| 国产精品伦人一区二区| 成年女人看的毛片在线观看| 欧美日韩视频高清一区二区三区二| 又爽又黄无遮挡网站| 久久久久久久久久久免费av| 97在线视频观看| 韩国高清视频一区二区三区| 在线a可以看的网站| xxx大片免费视频| 欧美xxxx性猛交bbbb| 国产伦理片在线播放av一区| 亚洲av.av天堂| 久久久久精品性色| 国产成人一区二区在线| 男女下面进入的视频免费午夜| 一个人观看的视频www高清免费观看| av在线app专区| 1000部很黄的大片| 女人十人毛片免费观看3o分钟| 国产综合懂色| 深夜a级毛片| 欧美zozozo另类| 亚洲人成网站高清观看| 97人妻精品一区二区三区麻豆| 天天躁日日操中文字幕| 亚洲三级黄色毛片| 亚洲天堂av无毛| 欧美成人一区二区免费高清观看| 免费av观看视频| 免费播放大片免费观看视频在线观看| 午夜激情福利司机影院| 色网站视频免费| 一级片'在线观看视频| 久久99热6这里只有精品| 久久女婷五月综合色啪小说 | 欧美成人午夜免费资源| 欧美精品一区二区大全| 色播亚洲综合网| 99久久人妻综合| 在线观看一区二区三区激情| 五月玫瑰六月丁香| 日韩制服骚丝袜av| 欧美成人a在线观看| 亚洲在线观看片| 亚洲精品亚洲一区二区| 免费观看性生交大片5| 国产乱人偷精品视频| 亚洲欧美成人精品一区二区| av在线蜜桃| 久久99热这里只有精品18| 久久久色成人| 免费av不卡在线播放| 欧美国产精品一级二级三级 | 国产免费视频播放在线视频| 亚洲精品一二三| 亚洲av成人精品一二三区| 国产精品嫩草影院av在线观看| 联通29元200g的流量卡| 视频区图区小说| 国产老妇伦熟女老妇高清| 内射极品少妇av片p| 国产成人免费观看mmmm| 女人被狂操c到高潮| 最近中文字幕2019免费版| 国产精品爽爽va在线观看网站| 丝瓜视频免费看黄片| 成人毛片60女人毛片免费| 久久综合国产亚洲精品| 国产日韩欧美亚洲二区| 日本午夜av视频| 亚洲精品日本国产第一区| 麻豆成人av视频| 亚洲一级一片aⅴ在线观看| 久久久精品免费免费高清| 国产毛片a区久久久久| 国产欧美另类精品又又久久亚洲欧美| 成年av动漫网址| 国内少妇人妻偷人精品xxx网站| 免费观看性生交大片5| 男人添女人高潮全过程视频| 亚洲欧美成人综合另类久久久| 亚洲国产精品成人久久小说| 免费av不卡在线播放| 国产视频首页在线观看| 一级爰片在线观看| 日韩精品有码人妻一区| 性色av一级| 国产av码专区亚洲av| 亚洲av电影在线观看一区二区三区 | 汤姆久久久久久久影院中文字幕| 好男人在线观看高清免费视频| 十八禁网站网址无遮挡 | 一区二区av电影网| 免费看av在线观看网站| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 在线精品无人区一区二区三 | 国产人妻一区二区三区在| 成人高潮视频无遮挡免费网站| 亚洲内射少妇av| 久久久久久久国产电影| 国产精品国产三级专区第一集| 一区二区三区免费毛片| 国产探花极品一区二区| 91久久精品国产一区二区三区| 新久久久久国产一级毛片| 一区二区三区精品91| 国内精品宾馆在线| 欧美日韩综合久久久久久| 免费大片18禁| 精品久久久久久久人妻蜜臀av| 少妇人妻久久综合中文| 国产91av在线免费观看| 亚洲熟女精品中文字幕| 国产精品一二三区在线看| 制服丝袜香蕉在线| 亚洲va在线va天堂va国产| 特级一级黄色大片| 男女国产视频网站| 成人综合一区亚洲| 久久久久久伊人网av| a级毛色黄片| 亚洲电影在线观看av| 午夜免费鲁丝| 午夜免费观看性视频| 亚洲国产精品专区欧美| 亚洲国产成人一精品久久久| 特大巨黑吊av在线直播| 成人亚洲欧美一区二区av| 国产精品爽爽va在线观看网站| 久久久a久久爽久久v久久| 国产一区二区三区av在线| 国产精品一区www在线观看| 亚洲av二区三区四区| 深爱激情五月婷婷| 人妻 亚洲 视频| 少妇人妻久久综合中文| 欧美成人午夜免费资源| 看非洲黑人一级黄片| 欧美xxxx性猛交bbbb| 全区人妻精品视频| 亚洲欧美日韩卡通动漫| 国产精品不卡视频一区二区| 国产精品三级大全| 天堂网av新在线| 日本wwww免费看| 超碰97精品在线观看| 国产免费一级a男人的天堂| 卡戴珊不雅视频在线播放| 国产伦理片在线播放av一区| 亚洲精品国产成人久久av| 国产亚洲av片在线观看秒播厂| 国产真实伦视频高清在线观看| 亚洲欧美中文字幕日韩二区| 天堂中文最新版在线下载 | 欧美老熟妇乱子伦牲交| 亚洲欧美成人精品一区二区| 日本wwww免费看| 亚洲国产色片| 全区人妻精品视频| 亚洲综合精品二区| 亚洲成色77777| 偷拍熟女少妇极品色| 国产成人a区在线观看| 欧美日韩视频精品一区| 国产人妻一区二区三区在| 啦啦啦啦在线视频资源| 中国美白少妇内射xxxbb| 国产成人免费观看mmmm| 又爽又黄a免费视频| 欧美bdsm另类| 少妇 在线观看| 国产av码专区亚洲av| 在线观看免费高清a一片| 特级一级黄色大片| 美女xxoo啪啪120秒动态图| 国产日韩欧美亚洲二区| 亚洲av.av天堂| 欧美高清性xxxxhd video| 国产黄频视频在线观看| 欧美丝袜亚洲另类| 成人毛片60女人毛片免费| 18禁裸乳无遮挡免费网站照片| 看十八女毛片水多多多| 日日摸夜夜添夜夜爱| 中国美白少妇内射xxxbb| 国产亚洲91精品色在线| 少妇 在线观看| 看免费成人av毛片| 亚洲国产av新网站| 好男人视频免费观看在线| 在线观看人妻少妇| 波野结衣二区三区在线| 菩萨蛮人人尽说江南好唐韦庄| 午夜福利视频1000在线观看| 亚洲人成网站高清观看| av网站免费在线观看视频| 亚洲av成人精品一二三区| 亚洲av电影在线观看一区二区三区 | 秋霞伦理黄片| 日日啪夜夜撸| 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品古装| 人体艺术视频欧美日本| 国产精品爽爽va在线观看网站| 国产一级毛片在线| 国产精品无大码| 99热这里只有是精品在线观看| 欧美精品人与动牲交sv欧美| 狂野欧美白嫩少妇大欣赏| 中文天堂在线官网| 亚洲精品日本国产第一区| 亚洲第一区二区三区不卡| 一本久久精品| 国产男女超爽视频在线观看| 国语对白做爰xxxⅹ性视频网站| 免费av不卡在线播放| 久久韩国三级中文字幕| 网址你懂的国产日韩在线| 啦啦啦中文免费视频观看日本| 特级一级黄色大片| 国精品久久久久久国模美| 人人妻人人爽人人添夜夜欢视频 | 只有这里有精品99| 又爽又黄无遮挡网站| 国产午夜精品久久久久久一区二区三区| 韩国av在线不卡| 国产黄色免费在线视频| 啦啦啦啦在线视频资源| 欧美日韩精品成人综合77777| 91精品国产九色| 久久国产乱子免费精品| 国产v大片淫在线免费观看| 国产黄频视频在线观看| 香蕉精品网在线| av国产免费在线观看| 欧美成人精品欧美一级黄| 中文欧美无线码| 热re99久久精品国产66热6| 日产精品乱码卡一卡2卡三| 久久久久性生活片| 成人黄色视频免费在线看| 91精品国产九色| 五月伊人婷婷丁香| 欧美成人午夜免费资源| 免费电影在线观看免费观看| 国产色婷婷99| 男女边吃奶边做爰视频| 日韩成人伦理影院| 视频中文字幕在线观看| 伦理电影大哥的女人| 免费观看无遮挡的男女| 成人亚洲欧美一区二区av| 精品久久久精品久久久| 我的老师免费观看完整版| 亚洲欧美中文字幕日韩二区| 亚洲精品日韩在线中文字幕| 婷婷色av中文字幕| 国产一区二区在线观看日韩| 97在线人人人人妻| 各种免费的搞黄视频| 精品久久久噜噜| 国产男女内射视频| 另类亚洲欧美激情| 精品少妇久久久久久888优播| 久久鲁丝午夜福利片| 少妇猛男粗大的猛烈进出视频 | 国产 精品1| 亚洲一区二区三区欧美精品 | 亚洲综合色惰| 中文字幕久久专区| 搡老乐熟女国产| 高清视频免费观看一区二区| av国产精品久久久久影院| 色综合色国产| 最后的刺客免费高清国语| 国产av码专区亚洲av| 99久国产av精品国产电影| 永久网站在线| 99热6这里只有精品| av播播在线观看一区| 欧美 日韩 精品 国产| 视频区图区小说| 新久久久久国产一级毛片| av又黄又爽大尺度在线免费看| 免费人成在线观看视频色| 日本猛色少妇xxxxx猛交久久| a级一级毛片免费在线观看| 丰满乱子伦码专区| 日本熟妇午夜| 一区二区三区乱码不卡18| 国产精品久久久久久久电影| 亚洲成人久久爱视频| 欧美亚洲 丝袜 人妻 在线| 又爽又黄无遮挡网站| 最新中文字幕久久久久| tube8黄色片| 国产亚洲91精品色在线| 日韩成人伦理影院| 少妇猛男粗大的猛烈进出视频 | 99热全是精品| 精品国产乱码久久久久久小说| 亚洲美女搞黄在线观看| 亚洲av福利一区| 日韩一区二区视频免费看| 日韩强制内射视频| 国产毛片a区久久久久| 一区二区三区免费毛片| 国产女主播在线喷水免费视频网站| 午夜激情福利司机影院| 3wmmmm亚洲av在线观看| tube8黄色片| 欧美国产精品一级二级三级 | 久久久欧美国产精品| 亚洲最大成人av| 99re6热这里在线精品视频| 日日摸夜夜添夜夜爱| 插逼视频在线观看| 在线免费十八禁| 午夜日本视频在线| 在线亚洲精品国产二区图片欧美 | 国产又色又爽无遮挡免| 亚洲精品久久午夜乱码| 国产高清不卡午夜福利| 毛片女人毛片| 欧美亚洲 丝袜 人妻 在线| 成人亚洲精品av一区二区| 亚洲国产精品专区欧美| 亚洲精品自拍成人| 少妇 在线观看| 免费观看a级毛片全部| 精品人妻一区二区三区麻豆| 国产一区二区亚洲精品在线观看| 人人妻人人爽人人添夜夜欢视频 | 久久久久久伊人网av| 国产淫语在线视频| 一级二级三级毛片免费看| 白带黄色成豆腐渣| 国内精品宾馆在线| 国产精品福利在线免费观看| 精品亚洲乱码少妇综合久久| 国产一区二区亚洲精品在线观看| 国产欧美亚洲国产| 亚洲av免费高清在线观看| 午夜福利在线在线| 别揉我奶头 嗯啊视频| 国产成人精品婷婷| 在线精品无人区一区二区三 | 又黄又爽又刺激的免费视频.| 精品一区在线观看国产| 女人久久www免费人成看片| 99久国产av精品国产电影| 汤姆久久久久久久影院中文字幕| 久久精品国产亚洲av天美| 五月玫瑰六月丁香| 高清视频免费观看一区二区| 一本久久精品| 蜜臀久久99精品久久宅男| 在线亚洲精品国产二区图片欧美 | 在线精品无人区一区二区三 | 国产日韩欧美在线精品| 超碰av人人做人人爽久久| 久久女婷五月综合色啪小说 | av.在线天堂| 51国产日韩欧美| 亚洲成人中文字幕在线播放| 久久久久久久亚洲中文字幕| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看| 亚洲成人久久爱视频| 精品久久久久久久久亚洲| 久久久久久久亚洲中文字幕| 在线观看美女被高潮喷水网站| 大陆偷拍与自拍| 2021少妇久久久久久久久久久| 免费看a级黄色片| 亚洲不卡免费看| 男女边吃奶边做爰视频| 国产白丝娇喘喷水9色精品| 国产精品久久久久久精品电影小说 | 九九久久精品国产亚洲av麻豆| 精品视频人人做人人爽| av在线播放精品| 免费观看在线日韩| 欧美3d第一页| 亚洲自偷自拍三级| 国产成人一区二区在线| 美女cb高潮喷水在线观看| 亚洲精品日韩在线中文字幕| 九九爱精品视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 日韩成人伦理影院| 男插女下体视频免费在线播放| av在线观看视频网站免费| 99久久人妻综合| 一级片'在线观看视频| 国产成人freesex在线| 男插女下体视频免费在线播放| 免费高清在线观看视频在线观看| 中文在线观看免费www的网站| 搡老乐熟女国产| 国产老妇女一区| 免费黄色在线免费观看| 国产成人91sexporn| 成人欧美大片| av卡一久久| 午夜福利高清视频| 丝瓜视频免费看黄片| a级毛片免费高清观看在线播放| 大香蕉久久网| 久久久久精品性色| 亚洲欧美成人综合另类久久久| 成人午夜精彩视频在线观看| 亚洲婷婷狠狠爱综合网| 亚洲一区二区三区欧美精品 | 午夜激情久久久久久久| 免费看光身美女| 性色avwww在线观看| 国产黄频视频在线观看| 国产成人a区在线观看| 听说在线观看完整版免费高清| 全区人妻精品视频| 国产精品伦人一区二区| 九色成人免费人妻av| 久久久久久久久久久免费av| 秋霞在线观看毛片| av在线观看视频网站免费| 性色av一级| 欧美bdsm另类| 乱码一卡2卡4卡精品| 在线观看国产h片| 高清午夜精品一区二区三区| 成年免费大片在线观看| 日本wwww免费看| 另类亚洲欧美激情| 亚洲综合精品二区| 欧美日韩视频精品一区| 免费人成在线观看视频色| 国产成人a∨麻豆精品| 日日啪夜夜撸| 亚洲av成人精品一区久久| 日韩成人av中文字幕在线观看| 成人毛片60女人毛片免费| 国产精品蜜桃在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲av二区三区四区| 三级国产精品片| av又黄又爽大尺度在线免费看| www.av在线官网国产| 亚洲国产精品成人久久小说| 国产伦精品一区二区三区视频9| 国产午夜精品一二区理论片| 亚洲欧美清纯卡通| eeuss影院久久| 看黄色毛片网站| 禁无遮挡网站| 久久久久精品久久久久真实原创| 成人午夜精彩视频在线观看| 午夜老司机福利剧场| 自拍偷自拍亚洲精品老妇| tube8黄色片| 丝袜美腿在线中文| 男女边摸边吃奶| tube8黄色片| 91久久精品国产一区二区成人| 中文字幕av成人在线电影| 国产亚洲av片在线观看秒播厂| 久久99热这里只频精品6学生| 国产永久视频网站| 国产大屁股一区二区在线视频| av国产精品久久久久影院| 男人舔奶头视频| 老女人水多毛片| 免费观看无遮挡的男女| 日产精品乱码卡一卡2卡三| 日本一本二区三区精品| 91久久精品国产一区二区成人| 免费高清在线观看视频在线观看| 日韩视频在线欧美| 卡戴珊不雅视频在线播放| 日本与韩国留学比较| 国产欧美日韩精品一区二区| 国语对白做爰xxxⅹ性视频网站| 国产一区二区在线观看日韩|