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

    基于內(nèi)聚力模型的配筋鋼纖維混凝土結(jié)構(gòu)承載性能研究

    2025-08-15 00:00:00孟棟棟黃逸群林本清
    重慶大學學報 2025年8期
    關(guān)鍵詞:本構(gòu)界面纖維

    中圖分類號:TU318 文獻標志碼:A 文章編號:1000-582X(2025)08-054-13

    doi:10.11835/j.issn.1000-582X.2024.275

    Study on the bearing capacity performance of steel fiber reinforced concrete structures based on the cohesive zone model

    MENG Dongdong',HUANG Yiqun2,LIN Benqing?

    (1.China Communications Construction Rail Transit Subsidiary,Beijing10l3oo,P.R.China;2.SchoolofCivil

    Engineering,F(xiàn)ujian UniversityofTechnology,F(xiàn)uzhou 350118,P.R.China; 3.Collge ofCivil Engineering, Fuzhou University, Fuzhou 350108,P.R. China)

    Abstract: This study proposes a numerical modeling method based on the cohesive zone model to investigate the mechanical response and fracture mechanism of steel fiber reinforced concrete (SFRC) structures.In the proposed model,cohesiveelementsareused tostimulate potential fracturesurfacesand rebar-concreteinterfaces.A constitutive model for SFRC fracture surfaces is developed by considering mixed-mode damage evolution, interfacial friction,and the fiber bridging effect.Aditionally,a modified bond-slipconstitutive model for the rebar-concrete interfaceis proposed,accounting fornormal separation.Tovalidate theproposed model,asriesof four-point bending experiments on SFRC specimens are conducted.The simulation results closely align with experimental observation,confirming the model’sability to accurately capture both mechanical response and fracture behavior.Parametric analysis reveals that inadequate fiber content or improper friction coefficients significantly reduce structural bearing capacity and ductility.

    Keywords:steel fiber reinforced concrete;cohesive zone model;constitutive model; fracture mechanics; fibei bridging effect; bond-slip relation

    鋼纖維混凝土作為普通混凝土的改進材料,指在常規(guī)混凝土內(nèi)部加入一定摻量的鋼纖維,使混凝土在抗拉強度、裂縫擴展控制以及韌性等方面得到了很好的提高,受到了廣泛的關(guān)注,并且已在工程中有所應(yīng)用。對于鋼纖維混凝土的本構(gòu)模型,目前研究主要可分為兩大類,一類是基于大量鋼纖維混凝土材料抗拉/抗壓試驗結(jié)果而得到相應(yīng)的單軸經(jīng)驗計算公式[或宏觀力學本構(gòu)[5;另一類則是從細觀出發(fā),專注研究纖維與混凝土間的相互作用,并結(jié)合試驗結(jié)果建立纖維在混凝土內(nèi)部拔出的黏結(jié)滑移計算方法[2]。上述本構(gòu)模型雖然能在某些具體問題中得到較好的預測效果,但由于鋼纖維混凝土結(jié)構(gòu)實際在斷裂破壞過程中可能同時會存在受壓/受拉破壞、受剪破壞以及復合破壞,且纖維的橋接作用還會增加這些破壞的復雜性,現(xiàn)有的本構(gòu)模型需要進一步改進才能表征這些破壞模式,相關(guān)研究仍有待開展。

    內(nèi)聚力模型,是一種基于內(nèi)聚力單元的有限元計算模型。內(nèi)聚力單元作為在實體單元中的內(nèi)插過渡單元,只需要定義法向及切向的應(yīng)力-應(yīng)變(或應(yīng)力-位移)關(guān)系,可較好地模擬界面斷裂損傷行為,已被廣泛的應(yīng)用于宏、細觀混凝土結(jié)構(gòu)的計算中,如Lopez等[1-4將雙線性損傷本構(gòu)引入內(nèi)聚力單元中,考慮砂漿基體、骨料、ITZ以及混凝土內(nèi)部的潛在斷裂面,建立了基于界面單元的混凝土細觀結(jié)構(gòu)計算模型,并較好地模擬了混凝土的細觀開裂行為。Wang等[15]、Huang等[6-1以及De Maio等[8]在這一模型的基礎(chǔ)上,做了大量參數(shù)分析,研究了骨料形狀、骨料含量、細觀缺陷、以及網(wǎng)格尺寸等因素對混凝土宏觀力學特性的影響。此外,內(nèi)聚力單元還可被用于表征鋼筋-混凝土界面,Chen等[19]、Syroka等[20]、De Maio等[2將規(guī)范中的經(jīng)典界面黏結(jié)-滑移關(guān)系引入內(nèi)聚力單元,并在二維模型中模擬了鋼筋混凝土梁的彎曲斷裂問題,且發(fā)現(xiàn)鋼筋-混凝土界面對模擬結(jié)果影響較大,不考慮黏結(jié)-滑移性能會導致計算結(jié)果被夸大。已有的研究表明,界面元模型能夠清晰地表現(xiàn)出混凝土內(nèi)部裂縫的形態(tài)及擴展規(guī)律,相較于經(jīng)典有限元方法,該方法克服了在界面附近區(qū)域網(wǎng)格容易畸形和模型網(wǎng)格數(shù)量大的問題,節(jié)省了計算資源,具有良好的應(yīng)用性。

    綜上,筆者基于內(nèi)聚力模型,提出了一種配筋鋼纖維混凝土結(jié)構(gòu)的建模方法,并且建立了相應(yīng)的鋼纖維混凝土潛在斷裂面本構(gòu)模型,該模型綜合考慮了損傷關(guān)系、縫間摩擦以及纖維的橋接效應(yīng)。此外,在補充考慮界面的法相分離作用基礎(chǔ)上,還建立了鋼筋-混凝土界面本構(gòu)模型。為驗證所建立模型正確性,設(shè)計并開展了1組配筋鋼纖維混凝土梁四點彎曲試驗。在驗證模型準確性后,針對所建立本構(gòu)模型的特點,還研究了不同纖維摻量,不同混凝土縫間摩擦系數(shù)對配筋鋼纖維混凝土梁承載能力的影響。

    1配筋鋼纖維混凝土結(jié)構(gòu)數(shù)值模型的建立

    為準確的模擬配筋鋼纖維混凝土結(jié)構(gòu),有兩個方面需要重點考慮,一是鋼纖維混凝土斷裂損傷的表征,包括混凝土的損傷開裂以及纖維的橋接作用(bridging effect);一是鋼筋與混凝土間的黏結(jié)滑移作用?;谝陨蟽煞矫?,結(jié)合內(nèi)聚力模型,建立了一種二維的配筋鋼纖維混凝土結(jié)構(gòu)的數(shù)值模型。

    1.1鋼纖維混凝土的建模及表征

    鋼纖維混凝土模型主要是基于內(nèi)聚力單元(cohesive element,COH2D4,ABAQUS)建立的,零厚度的內(nèi)聚力單元被插入于所有的實體單元(CPS3)交界面中,用以表征纖維混凝土的潛在斷裂面,如圖1所示。除此之外,所有的實體單元在模型中被假設(shè)為滿足線彈性關(guān)系且不可損傷。建模時,在單元選取上主要采用三角形單元而非四邊形單元,是因為采用三角形單元可以讓模型在開裂時產(chǎn)生比四邊形更為隨機平滑的開裂路徑,從而更好地模擬實際情況。

    圖1鋼纖維混凝土數(shù)值模型的建立

    1.2鋼筋與鋼纖維混凝土界面的建模及表征

    零厚度的內(nèi)聚力單元(COH2D4)在模型中同時還被用于表征鋼筋-混凝土界面,如圖2所示。鋼筋及鋼筋-混凝土界面的建模方式為:1)在混凝土模型中預先劃分出埋置鋼筋的區(qū)域;2)采用1.1節(jié)的方法,拆分基底混凝土,并全局插入零厚度內(nèi)聚力元,用以表征混凝土潛在斷裂面;3)基于重建后的混凝土模型以及預先劃分的鋼筋埋置區(qū)域,讀取并復制該區(qū)域的單元及結(jié)點信息,生成新的實體單元(CPS3),覆蓋于原混凝土模型之上(鋼筋與混凝土是重疊的,但二者間的單元互不干涉);4)通過零厚度的內(nèi)聚力單元(COH2D4),連接混凝土單元結(jié)點與覆蓋于其上的后建立的鋼筋單元結(jié)點,如圖2(a)所示。整個配筋鋼纖維混凝土建構(gòu)的建模流程可參照圖2(b)。綜上,模型的建立主要包含3部分,即鋼纖維混凝土建模(含潛在斷裂面)、鋼筋建模,以及鋼筋-混凝土界面的建模,其中,鋼纖維混凝土基體和鋼筋采用實體單元(CPS3),而鋼纖維混凝土潛在斷裂面和鋼筋-混凝土界面則在用零厚度內(nèi)聚力單元(COH2D4)表征。上述建模步驟難以直接在商用有限元軟件的窗口中直接實現(xiàn),因此,通過二次開發(fā),對有限元模型的前處理input文件進行修改而完成建模的。

    Fig.1 Establishmentofanumericalmodel forsteelfiberreinforcedconcrete

    圖2配筋鋼纖維混凝土建模方法

    Fig.2Modelingmethod forreinforced steel fiberconcrete

    上述建模方式建立的模型,可用于模擬分析配筋鋼纖維混凝土結(jié)構(gòu)平面斷裂破壞的全過程。對于在鋼纖維混凝土基礎(chǔ)上添加的鋼筋和鋼筋-混凝土界面,厚度均被設(shè)定為 n×πr(n 為鋼筋根數(shù), r 為鋼筋半徑),可以保證鋼筋上部和下部的界面的單位長度(延軸向)面積為 n×2πr ,從而與實際情況的界面面積一致。此外,為保證鋼筋的截面面積與實際情況一致,鋼筋的高度在模型中被設(shè)定為 r ,模型中的鋼筋截面積同樣為 n× πr2 。這種針對平面問題的鋼筋和鋼筋-混凝土界面幾何尺寸選取方式在已有研究中[19-21也得到了檢驗。需要注意的是,所建立的模型,基體的混凝土是連續(xù)的,并沒有被鋼筋所隔斷,因此,不同于傳統(tǒng)的連續(xù)單元有限元模型,本模型允許擴展中的裂縫在混凝土內(nèi)部自由貫穿鋼筋所在區(qū)域。

    2各關(guān)鍵組分的本構(gòu)模型

    2.1 鋼纖維混凝土的本構(gòu)模型

    鋼纖維混凝土在斷裂破壞時有3個主要方面需要著重考慮,即基體混凝土的斷裂損傷、混凝土界面間的摩擦作用,以及裂縫擴展后縫間纖維的橋接效應(yīng)。此外,由于鋼纖維混凝土的斷裂破壞主要通過模型中潛在斷裂面(內(nèi)聚力單元COH2D4)的力學行為來表示,因此,從上述3個方面出發(fā),基于內(nèi)聚力單元建立相應(yīng)的鋼纖維混凝土潛在斷裂面本構(gòu)模型。

    2.1.1 復合損傷關(guān)系

    對于混凝土基體,采用基于PL能量準則的雙線性復合損傷本構(gòu)模型[18.2-23]。該本構(gòu)模型基于如圖3(a)~(b)所示的單軸應(yīng)力-位移損傷關(guān)系,表達式為

    式中: σ 與 τ 為法向和切向應(yīng)力; δn 和 δs 為法向和切向位移; kn 和 ks 為法向和切向剛度,在描述混凝土界面時[13,16,18]假定 kn=ks ,根據(jù)已有研究[1-8],該值取值范圍通常在 105~106GPa/m δn0 和 δs0 為初始損傷位移; δnf 和 δsf 為失效位移; D 為損傷因子,可由圖3中幾何關(guān)系計算得到。

    圖3單軸條件下混凝土損傷關(guān)系

    Fig.3Damage relationshipofconcreteunderuniaxialconditions

    需要注意的是,在本模型中,當法向受壓時,潛在斷裂面被假定位滿足線彈性關(guān)系,即界面受壓不會產(chǎn)生損傷。

    在復合損傷模式下,采用平方準則判斷損傷的發(fā)生。

    式中: Ft 為抗拉強度; τ?0 為抗剪強度; 為MacaulayBracket算子,當算子內(nèi)數(shù)值小于0,對應(yīng)結(jié)果為0。在單調(diào)加載條件下,將式(2)代入式(3),則復合損傷模式下相對法向和切向初始損傷位移 (δn0r,δs0r) 為

    模型中潛在斷裂面的失效采用基于能量的PL準則進行判斷,判別式可寫為

    式中: Gr 和 GII 為法向和切向的斷裂能,可通過圖3中的幾何關(guān)系計算得到; GIr 和 GIIr 為加載過程中的能量釋放率,計算方式與 GrI,Gr 類似,僅需將相關(guān)變量替換為相對值[18,22-23]。

    基于PL準則,將能量及能量釋放率的表達式代入其中,則復合損傷模式下相對失效位移 (δnfr,δf0r) 為

    綜上,假定加載過程總位移 δ 、總初始損傷位移 δ?0 以及總失效位移 δf

    則復合損傷模式下?lián)p傷因子 D 可表示為

    2.1.2縫間摩擦作用

    當潛在斷裂面演化為裂縫并處于受壓狀態(tài)時,縫間會受到摩擦作用,在壓力較大的情況下,摩擦會極大地影響混凝土承載性能。因此,需要對摩擦作用進行細化考慮。模型主要考慮縫間不滑移以及縫間有滑移2種情況,將滑動摩擦產(chǎn)生的滑移量作用實時更新的計算變量,計算縫間的摩擦應(yīng)力[6,18]。

    1)不滑移情況。在縫間無滑移的狀態(tài)下,縫間的摩擦應(yīng)力必然小于最大摩擦應(yīng)力 Tfmax Tfmax 的表達式可根據(jù)摩擦應(yīng)力寫為

    式中:f為在斷裂面間的摩擦系數(shù)。

    摩擦應(yīng)力表達式可寫為

    Tf=ks(δssslide),(ksssslide|?Tfmax),

    式中: δsslide 為已產(chǎn)生的縫間相對滑移,該值在滑移情況下需要更新計算。

    2)有滑移情況。當縫間發(fā)生了滑移,則此時摩擦應(yīng)力值即為最大摩擦應(yīng)力,其表達式為

    需要注意的是,在有滑移情況下,當摩擦應(yīng)力計算完成后,還需更新已產(chǎn)生的滑移位移,更新表達式為

    式中: δsslide* 為更新后的滑移位移。

    2.1.3 橋接效應(yīng)

    當混凝土裂縫產(chǎn)生后,縫間的纖維會產(chǎn)生一定的“橋接效應(yīng)”,裂縫間此時仍會存在一定的牽引作用,如圖4所示。為方便計算,纖維的橋接效應(yīng)被等效為縫間的法向橋接應(yīng)力以及切向橋接應(yīng)力。

    圖4縫間纖維橋接效應(yīng)的等效

    綜合已有文獻[69-1024],將纖維混凝土受拉的應(yīng)力-位移關(guān)系扣去普通混凝土的應(yīng)力位移關(guān)系,斷裂面上端勾型纖維的宏觀橋接效應(yīng)可由圖5中的應(yīng)力-位移關(guān)系表示。在該應(yīng)力-位移關(guān)系中,橋接效應(yīng)被分成了3部分,分別為上升段、持平段以及下降段,這一關(guān)系與常見的鋼筋-混凝土黏結(jié)滑移關(guān)系[2)類似。橋接效應(yīng)的總應(yīng)力-位移關(guān)系可近似表示為

    式中: T? 為纖維橋接效應(yīng)總應(yīng)力,方向假設(shè)與斷裂面分離方向一致; δ 為潛在斷裂面總位移,采用式(6)計算;Tbmax 為最大橋接效應(yīng)應(yīng)力; δ1,δ2,δ3 為計算橋接應(yīng)力的關(guān)鍵參數(shù),為式(12)各計算段的邊界位移值; β 為函數(shù)下降段的形狀因子,是經(jīng)驗參數(shù),取 4δ2/(δ23) 。

    Fig.4Equivalence of fiberbridgingeffectbetweencracks

    圖5纖維在混凝土內(nèi)的宏觀橋接效應(yīng)應(yīng)力-位移關(guān)系

    Fig.5Stress-displacement relationship of the macroscopic bridging effect offibers in concrete

    為便于計算,假設(shè)纖維的橋接效應(yīng)只與混凝土縫間間距有關(guān),且橋接效應(yīng)的應(yīng)力與縫間相對位移的方向一致(通常纖維橋接效應(yīng)的承載力與裂縫分離的角度有關(guān),且受力方向并不完全與縫間相對位移的方向一致)。則法向橋接應(yīng)力 T?,n 和切向的橋接應(yīng)力 T?,s 可最終表示為

    2.1.4總應(yīng)力計算表達式

    綜合上述損傷、摩擦以及橋接效應(yīng)3個部分,則鋼纖維混凝土潛在斷裂面本構(gòu)模型最終可表示為

    需要注意的是,式(15)中 D?Tf 的物理含義為摩擦作用僅發(fā)生在已經(jīng)損傷的區(qū)域,而未損傷的區(qū)域裂縫尚未產(chǎn)生,自然沒有縫間摩擦應(yīng)力。

    2.2鋼筋-混凝土界面黏結(jié)滑移本構(gòu)模型

    對于鋼筋-混凝土界面的黏結(jié)滑移本構(gòu)模型,盡管采用的是纖維混凝土,纖維的作用會一定程度上提高鋼筋-混凝土界面的黏結(jié)性能,但從工程安全的角度,仍采用經(jīng)典的分段式本構(gòu)模型[20],如圖6所示,其表達式為

    式中: τmax 為黏結(jié)強度; τres 為殘余黏結(jié)強度; δ?s1,δ?s2 為與計算有關(guān)的位移參數(shù); a 為曲線上升段的形狀因子參數(shù)。

    圖6鋼筋-混凝土界面黏結(jié)滑移的應(yīng)力-位移關(guān)系

    Fig.6Stress-displacementrelationshipofbond-slipatthe steel-concrete interface

    對于帶肋鋼筋,根據(jù)CEP-FIB的推薦值[25],上述參數(shù)的取值為 τmax=fck0.5 =fck0.5,τres=0.15?τmax,δsl=0.6mm,δs2=2.5mm (204 a=0.4 。其中, fck 為基體混凝土的抗壓強度,單位為 MPa 。對于光圓鋼筋,則相應(yīng)的參數(shù)取值為 τmaxres=

    此外,由于經(jīng)典的黏結(jié)滑移本構(gòu)模型只給出了純切向的黏結(jié)-滑移關(guān)系,而法向分離對切向黏結(jié)性能的影響沒有被考慮,這可能會導致在計算過程中鋼筋對混凝土的約束作用被高估,從而使得計算誤差被放大。因此,在本構(gòu)模型中考慮法向分離的影響,令鋼筋-混凝土界面在法向上滿足式(1)中的變化規(guī)律,允許鋼筋-混凝土界面有法向的分離,如圖7(a)所示。此外,法向分離對界面切向黏結(jié)作用的削弱作用則通過黏結(jié)性能損失因子 D?2 表示,如圖7(b)所示,其表達式為

    式中: hrib 為鋼筋的肋高。該表達式的物理含義為,當鋼筋-混凝土界面的法向分離位移值超過鋼筋肋高時,則界面無任何黏結(jié)力。

    綜上,則考慮了法相分離后鋼筋-混凝土界面的黏結(jié)滑移應(yīng)力-位移關(guān)系式最終可以表示為

    圖7鋼筋-混凝土界面法向分離對黏結(jié)性能的影響

    Fig.7Effect of normal separation at the steel-concrete interface on bond performance

    圖8依據(jù)所建立本構(gòu)模型所編制材料子程序的程序框圖

    2.3本構(gòu)模型的編程實現(xiàn)

    需要注意的是,由于所建立的鋼纖維混凝土潛在斷裂面以及鋼筋-混凝土界面的本構(gòu)模型在ABAQUS中無法被直接提供,因此,上述本構(gòu)模型通過基于FORTRAN編程的用戶材料子程序VUMAT來實現(xiàn)。所編制的子程序?qū)?yīng)的程序流程圖如圖8所示。

    3模型的驗證及討論

    3.1配筋鋼纖維混凝土梁抗彎性能試驗

    為驗證所建立的配筋鋼纖維混凝土數(shù)值模型,基于端勾型鋼纖維,設(shè)計了一組配筋鋼纖維混凝土梁四點彎曲斷裂試驗,試件的設(shè)計及加載方式如圖9所示,基體為鋼纖維混凝土的試件命名為FPB-SFRC。此外,為體現(xiàn)鋼纖維的作用,在實驗中設(shè)計了素混凝土(C35強度)基體的鋼筋混凝土梁,命名為FPB-C35。

    在配筋方面,如圖9(a)(b)所示,梁上部配有2根直徑為 8mm 的帶肋鋼筋,下部配3根直徑為 14mm 的帶肋鋼筋,箍筋采用間距為 100mm 直徑為 8mm 的帶肋鋼筋。為更加明顯的體現(xiàn)鋼纖維對梁承載性能的影響,梁在純彎段中部不設(shè)箍筋,且純彎段上部同樣不設(shè)鋼筋。鋼筋標號統(tǒng)一選用HRB400。

    在材料方面,所選用的端勾型鋼纖維直徑為 0.75mm ,長度為 35mm ,彈性模量及抗拉強度分別為210GPa 和 1 100MPa 。纖維在混凝土中的摻量以纖維與混凝土的質(zhì)量比為度量,取質(zhì)量比 1% 。此外,經(jīng)測試,所澆筑的C35混凝土彈性模量、抗壓強度以及抗拉強度分別為30GPa,45MPa和 3MPa 。鋼筋的彈性模量以及屈服強度分別為 210GPa 和 400MPa 。

    Fig.9Testspecimendesignand loadingmethod

    3.2模型建立及驗證

    采用章節(jié)1的方法建立了配筋鋼纖維混凝土梁的數(shù)值模型,如圖10所示。模型采用基于ABAQUS/EXPLICITE求解器的準靜態(tài)方式進行求解。基于試驗結(jié)果、已有的參考文獻[13,6,8以及反復試算,模型所采用的主要計算參數(shù)如表1所示。除表1中主要參數(shù)外,用于計算鋼纖維橋接效應(yīng)的參數(shù)取值為: Tbmax= 1.16MPa δ1=0.2mm δ2=0.7mm δ?3=10mm ;用于計算鋼筋-混凝土界面黏結(jié)滑移關(guān)系的參數(shù)取值為: fck= 45MPa 。除了所提出的本構(gòu)模型外,鋼纖維混凝土的實體單元采用線彈性本構(gòu)模型,其屬性為: E=30 GPa,ν=0.2 。鋼筋實體單元采用理想彈塑性模型,其屬性為: E=210GPa,ν=0.3 ,屈服強度 fy=400MPa

    圖9試件設(shè)計及加載方式

    圖10配筋鋼纖維混凝土梁數(shù)值模型

    Fig.10Numerical model ofareinforced steel fiber concrete beam

    表1數(shù)值模型計算所用的主要參數(shù)

    Table1 Mainparametersused inthenumericalmodelcalculation

    經(jīng)過計算,模型計算結(jié)果以及試驗所得到的荷載位移曲線如圖11所示。由于鋼纖維的橋接效應(yīng),鋼纖維混凝土梁的承載能力以及延性均優(yōu)于普通混凝土梁,特別是當混凝土在頂部碎裂時(見圖12(b)),普通混凝土梁會瞬間喪失承載能力,而鋼纖維混凝土梁在頂部由于有纖維的限制,其在頂部壓碎時仍能保持原有的承載能力。此外,從圖11中可以看出,試驗所得的力學響應(yīng)在數(shù)值模型中也得到了較好的還原。

    圖11配筋鋼纖維混凝土梁四點彎曲模擬與試驗得到的荷載-位移曲線 Fig.11Load-displacement curve obtained from four-point bending simulation and testing ofa reinforced steel fiberconcrete beam

    圖12為數(shù)值模型得到的配筋鋼纖維混凝土梁裂縫分布與試驗結(jié)果的對比。由于在斷裂形態(tài)上普通混凝土與纖維混凝土所得到的結(jié)果類似,因此,選擇鋼纖維混凝土梁FPB-SFRC作為典型試件。對于數(shù)值模型,裂縫的表示方式為刪除損傷因子為1的內(nèi)聚力單元。從圖12中可以看出,模擬得到的梁最終斷裂形態(tài)與試驗結(jié)果類似,均表現(xiàn)為試件在底部純彎段發(fā)展出多條主裂縫并向上擴展,且頂部在加載到一定程度后會出現(xiàn)壓碎的現(xiàn)象。綜上,可以得出,所建立的數(shù)值模型可較好地模擬配筋鋼纖維混凝土結(jié)構(gòu)的承載性能以及相應(yīng)的斷裂形態(tài)。

    圖12配筋鋼纖維混凝土梁的裂縫分布

    Fig.12 Crackdistributioninareinforcedsteelfiberconcretebeam

    3.3 參數(shù)分析

    為進一步研究材料參數(shù)對配筋鋼纖維混凝土梁承載性能的影響,結(jié)合文中本構(gòu)模型的特點(纖維混凝土本構(gòu)模型中重點考慮了摩擦以及纖維的橋接效應(yīng)),以下將主要研究纖維摻量和混凝土摩擦系數(shù)取值對承載性能的影響。

    3.3.1 纖維摻量的影響

    根據(jù)已有研究[626],纖維摻量一般與橋接效應(yīng)中的最大橋接應(yīng)力 Tbmax 成正比。因此,可通過改變最大橋接應(yīng)力 Tbmax 來研究纖維摻量對結(jié)構(gòu)承載能力的影響,取3.2節(jié)中 0.0.5,1.1.5 倍的 Tbmax 進行計算,用以表示纖維質(zhì)量摻量為 0%0.5%1%1.5% 的情況。不同纖維摻量的計算結(jié)果如圖13所示,纖維摻量對結(jié)構(gòu)承載性能的影響十分顯著。首先,纖維摻量會直接影響梁結(jié)構(gòu)的最大承載力,且二者近似成正比關(guān)系。其次,纖維摻量還會影響梁結(jié)構(gòu)的延性和塑性,從圖中纖維質(zhì)量摻量為 0.5% 的結(jié)果可以看出,當混凝土頂部第一次發(fā)生碎裂,由于纖維摻量過少,頂部縫間的摩擦力和橋接效應(yīng)難以維持結(jié)構(gòu)承載力,梁結(jié)構(gòu)荷載會出現(xiàn)一定下降,且低于不摻纖維的試件,而當加載繼續(xù)進行到一定程度后,由于纖維摻量本身較少且橋接效應(yīng)所提供的應(yīng)力大幅降低,則試件承載力會再一次大幅下降;而對于纖維摻量合理的試件 (1%,1.5% 質(zhì)量比),當混凝王頂部碎裂時(見圖12(a)),頂部纖維提供的橋接效應(yīng)應(yīng)力能夠繼續(xù)維持結(jié)構(gòu)的結(jié)構(gòu)承載力,因此具有良好的延性。

    Fig.13Effect of different fibercontents on the mechanical properties of reinforced steel fiber concrete beams

    3.3.2混凝土摩擦系數(shù)

    為了解鋼纖維混凝土潛在斷裂面內(nèi)部摩擦作用對力學承載性能的影響,對不同摩擦系數(shù)(f為0、0.15、0.25、0.35、0.45)下配筋鋼纖維混凝土梁的四點彎曲試驗進行了計算,結(jié)果如圖14所示。從圖14可知,混凝土摩擦系數(shù)的取值會極大的影響梁結(jié)構(gòu)的斷裂形態(tài)和相應(yīng)的力學響應(yīng)。當摩擦系數(shù)取值過小,則梁加載至塑性階段時,隨著加載的進行,梁的撓度持續(xù)增大,頂部縫間盡管縫間存在纖維橋接效應(yīng),但由于縫間摩擦作用過低,二者的綜合作用仍不足以限制頂部已壓碎的混凝土塊產(chǎn)生滑移,而碎塊一旦滑移,頂部的局部抗剪及承壓能力就會隨之下降,無法與下部的抗拉能力相匹配,因此,梁會出現(xiàn)明顯的承載力下降。反之,當摩擦系數(shù)取值足夠大時,頂部由縫間摩擦作用以及橋接效應(yīng)所提供的局部抗剪及承壓能力能持續(xù)保持較大值,則此時結(jié)構(gòu)的承載力主要由鋼筋骨架決定,因此,圖中摩擦系數(shù) f=0.35 與 f=0.45 計算所得到的曲線具有高承載力與高延性,且二者幾乎重合。上述結(jié)果表明,在涉及計算或模擬中,纖維混凝土縫間的摩擦作用在纖維混凝土結(jié)構(gòu)中不可被忽視或簡化,特別是在分析結(jié)構(gòu)延性或峰后力學響應(yīng)的情況下。

    圖13不同纖維摻量對配筋鋼纖維混凝土梁力學性能的影響

    圖14不同摩擦系數(shù)取值對配筋鋼纖維混凝土梁力學性能的影響 Fig.14Effect of different friction coefficient values on the mechanical properties of reinforced steel fiberconcrete beams

    需要注意的是,在圖13和圖14中,荷載位移曲線均出現(xiàn)了承載能力下降后能在一定程度再次上升,這是由于在混凝土頂部壓碎后,頂部抗壓能力迅速下降導致了梁的承載能力迅速下降;但隨著加載點位移的進一步施加,混凝土頂部被壓碎的部分與原結(jié)構(gòu)形成了新的平衡關(guān)系,從而能夠繼續(xù)傳遞一部分的壓力和摩擦力,因此梁在荷載下降后能夠恢復少量的承載能力。

    4結(jié)論

    首先,基于內(nèi)聚力單元建立了一種配筋鋼纖維混凝土梁的有限元模型,該模型中內(nèi)聚力單元被用于表征鋼纖維混凝土的潛在斷裂面以及鋼筋-混凝土界面。在該模型的基礎(chǔ)上,建立了相應(yīng)的鋼纖維混凝土潛在斷裂面本構(gòu)模型以及改進的鋼筋-混凝土界面黏結(jié)滑移本構(gòu)模型。所建立的鋼纖維混凝土潛在斷裂面本構(gòu)綜合考慮了損傷作用、縫間摩擦作用以及纖維的橋接效應(yīng);而在改進的鋼筋-混凝土界面本構(gòu)模型中,補充考慮了界面的法向分離以及其對切向黏結(jié)滑移性能的影響。

    為驗證所建立數(shù)值模型的準確性,設(shè)計并開展了一組配筋鋼纖維混凝土梁四點彎曲試驗,鋼纖維選用端勾型纖維。經(jīng)過對比模擬與試驗結(jié)果,所建立的模型可較好地模擬配筋鋼纖維混凝土的全過程力學響應(yīng)及斷裂形態(tài),包括裂縫擴展、裂縫分布以及頂部混凝土的壓碎現(xiàn)象等。

    基于所建立材料本構(gòu)模型的特點,著重分析了纖維摻量以及混凝土摩擦系數(shù)取值對計算結(jié)果的影響,分析不同纖維摻量以及不同摩擦系數(shù)計算得到的結(jié)果,得到主要結(jié)論如下:

    1)纖維混凝土中纖維摻量會直接影響結(jié)構(gòu)的承載能力及延性,表現(xiàn)為在合理范圍內(nèi)纖維摻量越高,則結(jié)構(gòu)承載性能和延性就越好;此外,當纖維摻量不足時,可能會導致試件承載能力在塑性階段出現(xiàn)大幅下降,甚至低于不摻纖維的試件。

    2)混凝土縫間摩擦系數(shù)取值對纖維混凝土結(jié)構(gòu)的承載性能影響極大,主要體現(xiàn)在塑性段,當摩擦系數(shù)取值較小、或不考慮縫間摩擦作用時,與纖維摻量較少的情況類似,結(jié)構(gòu)的承載能力在塑性段會出現(xiàn)大幅下降,同時延性也會隨之減弱,因此,摩擦作用在分析這類結(jié)構(gòu)的延性時不可被忽視或被簡單考慮。

    參考文獻

    [1]范向前,胡少偉,朱海堂,等.齡期和試件尺寸對鋼纖維混凝土抗壓強度的影響[J].混凝土與水泥制品,2011(4):42-45. FanXQ,HuSW,ZhuHT,etalEfectsofageandspecimensizeonthecompresive strengthofsteelfiberreforced concrete[J]. China Concrete and Cement Products,2011(4): 42-45. (in Chinese)

    [2]關(guān)虓,牛獲濤,吳博,等.纖維混凝土受壓細觀統(tǒng)計損傷本構(gòu)模型研究[J].西安建筑科技大學學報(自然科學版),2015, 47(5):678-683. Guan X,NiuDTWuB,etal.Researchonesotochasticdamagecontitutivemodeloffberoncretedercomsion[J]. JournalofXi'anUniversityofArchitectureamp;Technology(NaturalScienceEdition),2015,47(5):678-683.(inChinese)

    [3]池寅,黃樂,余敏.基于ABAQUS的鋼-聚丙烯混雜纖維混凝土損傷塑性本構(gòu)模型取值方法研究[J].工程力學,2017, 34(12):131-142. ChiY,HuangL/y),YuM.Calibrationmetodof damageplasticitymodelforstel-polypropylenehybridfiberreiforced concrete based onAbaqus[J].EngineeringMechanics,2017,34(12):131-142.(in Chinese)

    [4]何雅杰.靜動載作用下鋼纖維混凝土力學性能試驗研究[D].淮南:安徽理工大學,2021. HeYJ.Experimental study on mechanical properties of steel fiber reinforced concrete under static anddynamic loads[D]. Huainan: Anhui University of Science amp; Technology,2021.(in Chinese) concrete[J].Engineering Fracture Mechanics,2016,159:35-62.

    [6]畢繼紅,霍琳穎,喬浩玥,等.單向受拉狀態(tài)下的鋼纖維混凝土本構(gòu)模型[J].工程力學,2020,37(6):155-164. BiJH,HuoLY,QiaoHY,etal.Aconstitutive modelofstelfiberreinforedconcreteunderuniaxialtension[J].Engieering Mechanics,2020,37(6):155-164.(in Chinese)

    [7]SujivorakulC,WaasA M,Naaman AE.Puloutresponseofasmothfiberwithan end anchorage[J].JournalofEngineering Mechanics,2000,126(9): 986-993.

    [8]RobinsP,AustinS,Jones PPul-outbehaviourofhookedstelfibes[J].Materialsand Structures,2002,35(7):434-442.

    [9]LaranjeiraF,AguadoA,MolinsC.Predicting thepulout responseofinclinedstraightsteelfbers[J].Materialsand Structures, 2010,43(6): 875-895.

    [10]LaranjeiraF,MolinsC,AguadoA.Predictingthepulloutresponseof inclinedhookedsteel fibers[J].Cementand Concrete Research,2010,40(10): 1471-1487.

    [11]SoetensT,VanGyseA,MathsS,etal.Asemi-analytical modeltopredictthepul-outbhaviourofinclinedhooked-end steel fibres[J]. Constructionand BuildingMaterials,2013,43:253-265.

    [12]Feng J,Sun WW, Wang X M,etal. Mechanical analysesof hooked fiber pullout performanceinultra-high-performance concrete[J].Construction and BuildingMaterials,2014,69: 403-410.

    [13]LópezCM,CarolI,AguadoA.Meso-structuralstudyofconcrete fractureusinginterfaceelements.I: numericalmodeland tensile behavior[J]. Materials and Structures.2008,41(3): 583-599.

    [14]LópezC,CarolI,AguadoA.Meso-structuralstudyofoncetefractreuinginterfaceelements.I:compresioiaial and Brazilian test[J].Materials and Structures,2008,41(3): 601-620.

    [15]Wang XF,Yang ZJ,YatesJR,et al.Monte Carlo simulations of mesoscale fracture modelingofconcrete withrandom aggregates and pores[J]. Construction and Building Materials,2015,75:35-45.

    [16]Huang YQ,HuSW.Acohesive modelforconcrete mesostructureconsidering friction efect betweencracks[J].Computers and Concrete,2019,24(1): 51-61.

    [17]Huang YQ,Zhang W,Liu X.Asessentofdiagonal macrocrack-induced debonding mechanisms inFRP-strengthened RC beams[J]. Journal of Composites for Construction,2022,26(5): 04022056.

    [18]DeMaio U,F(xiàn)abbrocinoF,GrecoF,etal.AstudyofconcretecoverseparationfailureinFRP-platedRCbeams viaanterelement fracture approach[J]. Composite Structures,2019,212: 625-636.

    [19]Chen GM,TengJG,ChenJF.Finite-elementmodelingof ntermediate crackdebondinginFRP-platedRCbeams[J].Joural of Composites for Construction,2011,15(3): 339-353.

    [20]SyrokaE,BobinskiJ,TejchmanJ.FEanalysisofreinforcedconcretecorbelswithenhancedcontinuummodels[J].Finite Elements in Analysis and Design,2011,47(9): 1066-1078.

    [21]DeMaioU,GrecoF,LeonetiL,etal.AninvestigationaboutdebondingmechanismsinFRP-strengthenedRCstructural elementsbyusingacohesive/volumetricmodeling technique[J].TheoreticalandApliedFractureMechanics,2022,117: 103199.

    [22]Zhang W,TangZZ,YangY,etal.AsessmentofFRP-conreteinterfacialdebondingwithcoupledmixed-modechsivezone model[J]. Journal of Composites for Construction,2021,25(2): 04021002.

    [23]張偉,唐站站,楊艷,等.復合受力模式下CFRP混凝土界面剝離分析方法[J].建筑結(jié)構(gòu)學報,2022,43(6):257-264,283. Zhang W,TangZZ,YangY,etal.Aalysis method fordebondingofCFRP-toconcrete interfaceundermixed-mode loading[J]. Journal of Building Structures,2022,43(6): 257-264,283.(in Chinese)

    [24]Lee Y,Kang ST,Kim JK.Pulloutbehaviorof inclinedsteefiber inanultra-highstrengthcementitiousmatrix[J]. Construction and Building Materials,2010,24(10): 2030-2041.

    [25]ComitéEuro-International du Béton.CEB-FIPmodelcode 1990:Design code[M].Thomas Telford Publishing,1993.

    [26]SoetensT,MathysS.Diferentmethodstomodelthepost-crackingbhaviourofhooked-endstelfibrereinforcedcocete[J]. Construction and Building Materials,2014,73: 458-471.

    (編輯 鄭潔)

    猜你喜歡
    本構(gòu)界面纖維
    挑選抽紙須注意幾點
    酶促打漿對竹漿性能的影響研究
    中國造紙(2025年7期)2025-08-18 00:00:00
    彈塑性力學內(nèi)容體系梳理和課程思政建設(shè)探究
    基于細觀仿真建模的CFRP縱-扭超聲銑削材料去除機理研究
    C/SiC復合材料激光超聲復合加工螺紋缺陷控制及拉伸性能分析
    國企黨委前置研究的“四個界面”
    當代陜西(2020年13期)2020-08-24 08:22:02
    基于FANUC PICTURE的虛擬軸坐標顯示界面開發(fā)方法研究
    人機交互界面發(fā)展趨勢研究
    手機界面中圖形符號的發(fā)展趨向
    新聞傳播(2015年11期)2015-07-18 11:15:04
    日韩强制内射视频| 欧美xxxx黑人xx丫x性爽| 亚洲精品第二区| 看十八女毛片水多多多| 亚洲成人中文字幕在线播放| 国产精品女同一区二区软件| 欧美xxxx性猛交bbbb| 一级爰片在线观看| 伊人久久精品亚洲午夜| 国内少妇人妻偷人精品xxx网站| 亚洲av.av天堂| 一个人观看的视频www高清免费观看| 超碰97精品在线观看| 午夜视频国产福利| 日日啪夜夜撸| av在线天堂中文字幕| 狂野欧美激情性bbbbbb| 国产欧美日韩精品一区二区| 1000部很黄的大片| 有码 亚洲区| 国产亚洲精品久久久com| 欧美高清性xxxxhd video| 日本熟妇午夜| 搡女人真爽免费视频火全软件| 嫩草影院精品99| 下体分泌物呈黄色| 欧美97在线视频| 欧美区成人在线视频| 成人国产麻豆网| 亚洲aⅴ乱码一区二区在线播放| 成人鲁丝片一二三区免费| 国产午夜精品一二区理论片| 一区二区三区精品91| 亚洲精品国产成人久久av| 18禁裸乳无遮挡免费网站照片| 丝瓜视频免费看黄片| 国产欧美亚洲国产| 久久久久精品久久久久真实原创| 91久久精品国产一区二区成人| 亚洲精品中文字幕在线视频 | 国产精品精品国产色婷婷| 国产精品99久久99久久久不卡 | 国产精品女同一区二区软件| a级毛色黄片| 久久久精品免费免费高清| 人妻夜夜爽99麻豆av| 免费不卡的大黄色大毛片视频在线观看| 性插视频无遮挡在线免费观看| 亚洲av不卡在线观看| 欧美激情国产日韩精品一区| 看非洲黑人一级黄片| 日本爱情动作片www.在线观看| 国产精品国产三级国产专区5o| 午夜视频国产福利| 午夜福利在线在线| 最近中文字幕高清免费大全6| 国产69精品久久久久777片| 亚洲精品一二三| 亚洲精品视频女| 夫妻午夜视频| 99热全是精品| 欧美激情久久久久久爽电影| 26uuu在线亚洲综合色| 国产av国产精品国产| 91精品国产九色| 亚洲在久久综合| 国产午夜精品久久久久久一区二区三区| 激情 狠狠 欧美| 97精品久久久久久久久久精品| 日韩一区二区视频免费看| 日韩欧美 国产精品| 欧美精品国产亚洲| 国产精品一区www在线观看| 免费人成在线观看视频色| 午夜福利高清视频| 国产精品女同一区二区软件| 丝瓜视频免费看黄片| 国产男女超爽视频在线观看| 日韩三级伦理在线观看| 色婷婷久久久亚洲欧美| 久热这里只有精品99| 国产爱豆传媒在线观看| av一本久久久久| 国产中年淑女户外野战色| 国产在线男女| 亚洲欧美精品专区久久| 国产免费一区二区三区四区乱码| 美女视频免费永久观看网站| 成年版毛片免费区| 国产一区亚洲一区在线观看| 一级黄片播放器| 国产毛片在线视频| 99久久精品热视频| 国内精品宾馆在线| 亚洲无线观看免费| 亚洲av成人精品一二三区| 午夜老司机福利剧场| 国产大屁股一区二区在线视频| 一级毛片黄色毛片免费观看视频| 各种免费的搞黄视频| 日韩亚洲欧美综合| 国产精品国产av在线观看| 国产精品国产三级专区第一集| 亚洲av一区综合| 成人鲁丝片一二三区免费| 岛国毛片在线播放| 久久久久久久久大av| 春色校园在线视频观看| 人人妻人人看人人澡| 亚洲真实伦在线观看| 亚洲av一区综合| 99热这里只有精品一区| 国产精品一区二区性色av| 男女边吃奶边做爰视频| 欧美 日韩 精品 国产| 伦理电影大哥的女人| 亚洲国产av新网站| 熟女人妻精品中文字幕| 国产精品人妻久久久久久| 国产精品无大码| 色播亚洲综合网| 观看免费一级毛片| 欧美精品一区二区大全| 欧美激情在线99| 久久久久久久久久久丰满| 夜夜看夜夜爽夜夜摸| 亚洲内射少妇av| 亚洲精品第二区| 日韩,欧美,国产一区二区三区| 亚洲欧美成人综合另类久久久| 春色校园在线视频观看| 麻豆成人av视频| 久久精品久久久久久噜噜老黄| 日韩免费高清中文字幕av| 在线免费观看不下载黄p国产| kizo精华| 美女国产视频在线观看| 丝袜脚勾引网站| 亚洲精品一二三| 26uuu在线亚洲综合色| 日韩免费高清中文字幕av| 美女被艹到高潮喷水动态| 神马国产精品三级电影在线观看| 久久精品久久精品一区二区三区| 久久久久久久久久久丰满| 麻豆久久精品国产亚洲av| 国产国拍精品亚洲av在线观看| 国产美女午夜福利| 亚洲色图av天堂| 日韩大片免费观看网站| 亚洲欧美一区二区三区国产| 欧美+日韩+精品| 少妇丰满av| 最近手机中文字幕大全| 欧美高清成人免费视频www| 精华霜和精华液先用哪个| 国模一区二区三区四区视频| a级一级毛片免费在线观看| 亚洲精品日韩av片在线观看| 日本欧美国产在线视频| 亚洲欧美日韩东京热| 免费电影在线观看免费观看| 久久人人爽人人片av| 女人久久www免费人成看片| 亚洲精品乱久久久久久| 成人二区视频| 成人毛片a级毛片在线播放| 观看美女的网站| 精品人妻一区二区三区麻豆| 欧美日韩在线观看h| 日韩欧美精品免费久久| 国产探花在线观看一区二区| 婷婷色综合大香蕉| 久久久久精品性色| 夫妻性生交免费视频一级片| 丰满乱子伦码专区| 最近最新中文字幕大全电影3| 亚洲激情五月婷婷啪啪| 18禁裸乳无遮挡免费网站照片| 国产在线男女| 我的女老师完整版在线观看| 日韩一区二区三区影片| 日本免费在线观看一区| 91久久精品国产一区二区三区| 如何舔出高潮| 高清午夜精品一区二区三区| 成人特级av手机在线观看| 欧美 日韩 精品 国产| 精品熟女少妇av免费看| 久久99热这里只频精品6学生| videos熟女内射| 神马国产精品三级电影在线观看| 久久午夜福利片| 欧美潮喷喷水| 男女下面进入的视频免费午夜| 老司机影院成人| 久久久久久久久久成人| 五月玫瑰六月丁香| 亚洲精品第二区| 波野结衣二区三区在线| 在线观看国产h片| 欧美成人午夜免费资源| 日韩人妻高清精品专区| 午夜老司机福利剧场| 亚洲精品一二三| 99热网站在线观看| 亚洲精品中文字幕在线视频 | 亚洲在线观看片| av播播在线观看一区| 午夜老司机福利剧场| 熟女av电影| 免费av毛片视频| 狂野欧美白嫩少妇大欣赏| 国产成人精品福利久久| 在线观看国产h片| 涩涩av久久男人的天堂| 国产 一区精品| 日韩亚洲欧美综合| 久久久久久久久久人人人人人人| 国产综合精华液| 老师上课跳d突然被开到最大视频| 男人舔奶头视频| 午夜老司机福利剧场| 91狼人影院| 亚洲精品色激情综合| 干丝袜人妻中文字幕| 国产av码专区亚洲av| 国产高清国产精品国产三级 | av.在线天堂| 午夜亚洲福利在线播放| 久热这里只有精品99| 九九久久精品国产亚洲av麻豆| 一级黄片播放器| 99视频精品全部免费 在线| 国产精品爽爽va在线观看网站| 国产日韩欧美亚洲二区| 在线观看一区二区三区| 国产亚洲一区二区精品| 亚洲激情五月婷婷啪啪| 精品一区二区三区视频在线| 你懂的网址亚洲精品在线观看| 综合色丁香网| 国产淫语在线视频| 中国三级夫妇交换| 日本三级黄在线观看| 91午夜精品亚洲一区二区三区| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品999| 亚州av有码| 欧美极品一区二区三区四区| 亚洲欧美日韩无卡精品| 久久久久国产精品人妻一区二区| 在线天堂最新版资源| 青青草视频在线视频观看| 夜夜爽夜夜爽视频| 少妇人妻久久综合中文| 女人被狂操c到高潮| 三级国产精品欧美在线观看| 亚洲va在线va天堂va国产| 十八禁网站网址无遮挡 | 亚洲国产精品专区欧美| 日韩一区二区三区影片| 男女边摸边吃奶| 亚洲综合精品二区| 人人妻人人澡人人爽人人夜夜| 舔av片在线| 国产高清三级在线| 国产伦精品一区二区三区视频9| 成年免费大片在线观看| 国产爱豆传媒在线观看| 国产片特级美女逼逼视频| 成人亚洲精品一区在线观看 | 天美传媒精品一区二区| 久久久欧美国产精品| 我的女老师完整版在线观看| 亚洲美女搞黄在线观看| 街头女战士在线观看网站| 春色校园在线视频观看| 大陆偷拍与自拍| 免费观看性生交大片5| 国产伦精品一区二区三区四那| 日韩中字成人| 建设人人有责人人尽责人人享有的 | 欧美丝袜亚洲另类| 色网站视频免费| av在线天堂中文字幕| 欧美97在线视频| 日本wwww免费看| 超碰97精品在线观看| 日韩人妻高清精品专区| 国产亚洲5aaaaa淫片| 久久久久九九精品影院| 亚洲成色77777| 嫩草影院新地址| 新久久久久国产一级毛片| 国产精品偷伦视频观看了| 午夜爱爱视频在线播放| 三级国产精品欧美在线观看| 禁无遮挡网站| 九九在线视频观看精品| 校园人妻丝袜中文字幕| 国产精品国产av在线观看| xxx大片免费视频| 直男gayav资源| 2022亚洲国产成人精品| 下体分泌物呈黄色| 身体一侧抽搐| 99久国产av精品国产电影| 国产黄片美女视频| 69人妻影院| 在线播放无遮挡| 国产亚洲午夜精品一区二区久久 | 国产高清国产精品国产三级 | 亚洲av欧美aⅴ国产| 亚洲va在线va天堂va国产| 网址你懂的国产日韩在线| 美女xxoo啪啪120秒动态图| 成人欧美大片| 国产男女内射视频| 国产成人a区在线观看| 亚洲综合色惰| 极品少妇高潮喷水抽搐| 国产精品av视频在线免费观看| 久热久热在线精品观看| 亚洲精品日韩av片在线观看| 极品教师在线视频| 精华霜和精华液先用哪个| 青春草国产在线视频| 亚洲精品亚洲一区二区| 少妇 在线观看| 国内少妇人妻偷人精品xxx网站| 午夜福利视频1000在线观看| 全区人妻精品视频| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区三区四区免费观看| 91在线精品国自产拍蜜月| 亚洲国产精品999| 亚洲最大成人av| 美女被艹到高潮喷水动态| 国产成人精品婷婷| 男人添女人高潮全过程视频| 麻豆国产97在线/欧美| 天堂网av新在线| 亚州av有码| 精品人妻一区二区三区麻豆| 99热网站在线观看| 久久久a久久爽久久v久久| 99久久精品热视频| 乱码一卡2卡4卡精品| 精品国产三级普通话版| 欧美成人午夜免费资源| 美女国产视频在线观看| 国产大屁股一区二区在线视频| 精品一区二区免费观看| 亚洲图色成人| 成人国产av品久久久| 3wmmmm亚洲av在线观看| 亚洲国产成人一精品久久久| 欧美bdsm另类| 欧美高清性xxxxhd video| 国产一区二区在线观看日韩| 蜜桃亚洲精品一区二区三区| 深夜a级毛片| 麻豆精品久久久久久蜜桃| 亚洲精品乱码久久久久久按摩| 黄色视频在线播放观看不卡| 成年女人看的毛片在线观看| 美女cb高潮喷水在线观看| 一级毛片黄色毛片免费观看视频| 大香蕉97超碰在线| 麻豆久久精品国产亚洲av| 亚洲最大成人手机在线| 中文天堂在线官网| 亚洲欧美中文字幕日韩二区| 91午夜精品亚洲一区二区三区| 最近2019中文字幕mv第一页| 汤姆久久久久久久影院中文字幕| 超碰97精品在线观看| 中文欧美无线码| 免费看日本二区| 欧美日韩视频精品一区| 中文字幕亚洲精品专区| 久久久成人免费电影| 99久久人妻综合| 国产精品秋霞免费鲁丝片| 国产久久久一区二区三区| 日韩制服骚丝袜av| 国产大屁股一区二区在线视频| 精品一区二区免费观看| 一级爰片在线观看| 精品午夜福利在线看| 欧美一级a爱片免费观看看| 一级片'在线观看视频| 午夜精品一区二区三区免费看| 亚洲自偷自拍三级| 国产免费又黄又爽又色| 综合色av麻豆| 在线精品无人区一区二区三 | 性色av一级| 国产精品福利在线免费观看| 欧美极品一区二区三区四区| 精品久久久久久电影网| 在线观看一区二区三区激情| 一级片'在线观看视频| 日韩欧美一区视频在线观看 | 国产精品久久久久久精品电影| 亚洲国产精品国产精品| 国产成人91sexporn| av线在线观看网站| xxx大片免费视频| 精品一区二区三卡| 搡老乐熟女国产| 91狼人影院| av国产精品久久久久影院| 九草在线视频观看| 国模一区二区三区四区视频| 麻豆久久精品国产亚洲av| av女优亚洲男人天堂| 亚洲精品乱码久久久久久按摩| 一级黄片播放器| 精品人妻视频免费看| 三级国产精品片| 免费观看a级毛片全部| 久热久热在线精品观看| 国产伦精品一区二区三区视频9| 国产大屁股一区二区在线视频| 亚洲精品乱久久久久久| 秋霞伦理黄片| 天美传媒精品一区二区| 亚洲av二区三区四区| 丰满少妇做爰视频| 欧美精品人与动牲交sv欧美| 久久99热6这里只有精品| 久久热精品热| 六月丁香七月| 99re6热这里在线精品视频| 又粗又硬又长又爽又黄的视频| 又爽又黄a免费视频| 亚洲美女视频黄频| 青春草国产在线视频| 亚洲av国产av综合av卡| 国产亚洲精品久久久com| 日韩国内少妇激情av| 日韩强制内射视频| 免费人成在线观看视频色| 一本色道久久久久久精品综合| 免费看av在线观看网站| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| av在线亚洲专区| 麻豆国产97在线/欧美| 又爽又黄a免费视频| 深夜a级毛片| 丝袜脚勾引网站| 亚洲成人av在线免费| 国产高清有码在线观看视频| 国产成人一区二区在线| 国产成人福利小说| 久久精品国产亚洲av天美| 国产精品国产三级专区第一集| 久久午夜福利片| .国产精品久久| 插阴视频在线观看视频| 亚洲精品乱码久久久久久按摩| 啦啦啦啦在线视频资源| 国产伦精品一区二区三区视频9| 久久久久久久亚洲中文字幕| 国国产精品蜜臀av免费| 十八禁网站网址无遮挡 | 我的女老师完整版在线观看| 亚洲精品一区蜜桃| 亚洲精品乱久久久久久| 国产精品国产av在线观看| 午夜亚洲福利在线播放| 高清午夜精品一区二区三区| 狂野欧美激情性bbbbbb| 精品国产乱码久久久久久小说| videos熟女内射| 久久久久久国产a免费观看| 亚洲,欧美,日韩| 能在线免费看毛片的网站| 高清午夜精品一区二区三区| 韩国av在线不卡| 欧美日韩精品成人综合77777| 免费黄频网站在线观看国产| 欧美国产精品一级二级三级 | 赤兔流量卡办理| 少妇裸体淫交视频免费看高清| 少妇高潮的动态图| 国产片特级美女逼逼视频| 亚洲真实伦在线观看| 别揉我奶头 嗯啊视频| 欧美精品国产亚洲| 国产成人一区二区在线| 18禁裸乳无遮挡免费网站照片| 午夜福利高清视频| 午夜精品国产一区二区电影 | 99久久人妻综合| 午夜激情久久久久久久| 不卡视频在线观看欧美| 久久久久久久精品精品| 国产伦在线观看视频一区| 国产精品爽爽va在线观看网站| 可以在线观看毛片的网站| 国产成人免费观看mmmm| 亚洲不卡免费看| 又大又黄又爽视频免费| 热re99久久精品国产66热6| 成人毛片a级毛片在线播放| 在线 av 中文字幕| 国产午夜福利久久久久久| 好男人视频免费观看在线| 十八禁网站网址无遮挡 | 亚洲一区二区三区欧美精品 | 亚洲av福利一区| av在线亚洲专区| 性色avwww在线观看| 午夜福利高清视频| 免费黄网站久久成人精品| 在线观看三级黄色| 久久久久久久久久成人| 大香蕉97超碰在线| 黄片无遮挡物在线观看| 一级爰片在线观看| 成人漫画全彩无遮挡| 老司机影院毛片| 国产中年淑女户外野战色| 哪个播放器可以免费观看大片| 在线观看一区二区三区激情| 欧美一区二区亚洲| 三级男女做爰猛烈吃奶摸视频| 少妇人妻精品综合一区二区| 久久久久精品久久久久真实原创| 日韩欧美一区视频在线观看 | 免费大片18禁| 国产精品久久久久久久电影| 久久人人爽人人片av| 国产精品久久久久久精品电影| 午夜视频国产福利| av卡一久久| 午夜福利网站1000一区二区三区| 欧美成人精品欧美一级黄| 国产成人精品婷婷| 99re6热这里在线精品视频| 国内精品宾馆在线| 街头女战士在线观看网站| 中文在线观看免费www的网站| 搡老乐熟女国产| 伦理电影大哥的女人| 国产美女午夜福利| 久久精品国产自在天天线| 国产探花在线观看一区二区| 亚洲国产精品成人综合色| 国产av码专区亚洲av| 欧美丝袜亚洲另类| 亚洲综合色惰| 久久亚洲国产成人精品v| 春色校园在线视频观看| 亚洲精品,欧美精品| tube8黄色片| 在线观看美女被高潮喷水网站| 寂寞人妻少妇视频99o| 免费看av在线观看网站| 七月丁香在线播放| 国产成人一区二区在线| 一本久久精品| 一区二区三区精品91| 最近最新中文字幕大全电影3| 99九九线精品视频在线观看视频| 免费观看av网站的网址| 在线免费观看不下载黄p国产| 亚洲国产av新网站| 久久久欧美国产精品| 国模一区二区三区四区视频| 街头女战士在线观看网站| 亚洲人成网站在线播| 免费电影在线观看免费观看| 久久精品熟女亚洲av麻豆精品| 观看美女的网站| 大话2 男鬼变身卡| 免费在线观看成人毛片| 精品人妻视频免费看| 久久精品国产a三级三级三级| 听说在线观看完整版免费高清| 最近中文字幕2019免费版| 超碰97精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 免费大片黄手机在线观看| 国产精品99久久99久久久不卡 | 欧美老熟妇乱子伦牲交| 亚洲欧美清纯卡通| 亚洲精品日韩av片在线观看| 国产伦精品一区二区三区四那| 日韩欧美精品免费久久| av在线观看视频网站免费| 观看免费一级毛片| 日本爱情动作片www.在线观看| 国产乱人视频| 亚洲天堂国产精品一区在线| 日本爱情动作片www.在线观看| 欧美精品一区二区大全| 高清在线视频一区二区三区| 久久人人爽av亚洲精品天堂 | 麻豆国产97在线/欧美| .国产精品久久| 久久精品久久久久久久性| 一边亲一边摸免费视频| 精品酒店卫生间| 青春草国产在线视频| xxx大片免费视频| 亚洲成色77777| 免费看av在线观看网站| 亚洲内射少妇av| 国产熟女欧美一区二区| 女人被狂操c到高潮|