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

    考慮多因素作用下的RC梁抗力時空退化模型

    2017-11-02 10:46:40陽逸鳴彭建新張建仁
    關(guān)鍵詞:抗力氯離子變異性

    陽逸鳴, 彭建新, 張建仁

    (橋梁工程安全控制省部共建教育部重點實驗室(長沙理工大學(xué)), 長沙 410114)

    考慮多因素作用下的RC梁抗力時空退化模型

    陽逸鳴, 彭建新, 張建仁

    (橋梁工程安全控制省部共建教育部重點實驗室(長沙理工大學(xué)), 長沙 410114)

    為研究鋼筋混凝土結(jié)構(gòu)銹蝕后的力學(xué)性能和抗力退化規(guī)律,基于RC梁加速銹蝕試驗結(jié)果,首先綜合考慮均勻銹蝕和點蝕同時發(fā)生的影響,對鋼筋截面面積的時變模型進(jìn)行分析,其次探討裂縫寬度與氯離子擴(kuò)散系數(shù)及RC梁抗力的關(guān)系,隨后考慮受拉主筋銹蝕不均衡的影響,對結(jié)構(gòu)抗力計算模型進(jìn)行討論,最后基于材料性能和結(jié)構(gòu)尺寸的空間變異性,建立RC梁抗力時空退化模型,并對橋梁服役100 a后的抗力進(jìn)行預(yù)測. 研究發(fā)現(xiàn):當(dāng)考慮各參數(shù)的空間變異性時,不考慮裂縫影響得到的抗力值比考慮裂縫影響得到的抗力值高了9.91%;考慮各參數(shù)空間變異性時的抗力值比不考慮時下降了14.68%;適當(dāng)?shù)倪x取波動系數(shù)和單元尺寸對考慮結(jié)構(gòu)各參數(shù)的空間變異性十分重要;潮汐區(qū)和浪濺區(qū)環(huán)境下結(jié)構(gòu)的剩余抗力分別為初始抗力的43.23%和36.45%,較海岸線大氣區(qū)環(huán)境下的抗力值分別下降了6.73%和21.35%,因此在結(jié)構(gòu)服役期間,應(yīng)盡量做好防腐措施,減少氯離子侵蝕對結(jié)構(gòu)耐久性的破壞.

    銹蝕梁;抗力退化;鋼筋截面面積;裂縫;銹蝕不均衡;空間變異性

    氯離子侵蝕引起的鋼筋混凝土結(jié)構(gòu)抗力及耐久性下降的問題逐漸成為研究熱點[1-3]. 在結(jié)構(gòu)抗力研究中,銹蝕鋼筋截面面積計算模型主要采用均勻銹蝕模型和點蝕模型[4-5],并未綜合考慮均勻銹蝕和坑蝕兩種狀態(tài)同時發(fā)生時對鋼筋截面面積變化的影響. 同時,關(guān)于裂縫影響的研究大多集中于建立裂縫寬度或鋼筋銹蝕率的聯(lián)系,如文獻(xiàn)[6]對鋼筋開裂后的銹蝕產(chǎn)物體積進(jìn)行計算,基于修正后的裂縫寬度計算模型,得到混凝土開裂后鋼筋的銹蝕率;文獻(xiàn)[7]基于45片RC構(gòu)件腐蝕試驗結(jié)果,建立了銹脹開裂寬度的預(yù)測模型,并對各參數(shù)進(jìn)行敏感性分析. 很少有學(xué)者在理論上建立銹脹裂縫與結(jié)構(gòu)抗力的聯(lián)系. 以往鋼筋混凝土結(jié)構(gòu)可靠度的研究大多假定材料性能及結(jié)構(gòu)參數(shù)在空間上分布均勻,只考慮材料參數(shù)的時變性和不確定性,不考慮各參數(shù)的空間變異性對結(jié)構(gòu)性能退化的影響. 但實際結(jié)構(gòu)中,不同位置的材料性能和結(jié)構(gòu)尺寸可能存在差異,最后造成結(jié)構(gòu)的安全評價偏向保守. 隨著可靠度研究的深入,考慮參數(shù)的空間變異性對結(jié)構(gòu)可靠度的研究越來越迫切. 目前,有極少數(shù)學(xué)者對結(jié)構(gòu)的空間變異性展開研究. 文獻(xiàn)[8]基于混凝土抗壓強(qiáng)度、表面氯離子濃度質(zhì)量和保護(hù)層厚度三者的空間變異性,建立一維時空模型,并對鋼筋混凝土結(jié)構(gòu)的承載能力和銹脹開裂損傷進(jìn)行系統(tǒng)的分析,同時針對性的提出相關(guān)維護(hù)方法;文獻(xiàn)[9]綜合考慮多種參數(shù)的空間變異性,研究了氯離子侵蝕環(huán)境下結(jié)構(gòu)的開裂比例;文獻(xiàn)[10]基于氯離子擴(kuò)散試驗數(shù)據(jù)得到氯離子擴(kuò)散系數(shù)計算模型,通過引入?yún)?shù)空間變異性,建立RC結(jié)構(gòu)初始銹蝕模型. 然而以上學(xué)者并未考慮同一RC梁中不同受拉主筋銹蝕不均衡以及裂縫對RC梁抗力的影響.

    因此本文以RC梁加速銹蝕試驗為基礎(chǔ),首先綜合考慮均勻銹蝕和點蝕的影響,得到改進(jìn)的鋼筋截面面積時變模型;隨后考慮裂縫的影響,建立裂縫寬度與氯離子擴(kuò)散系數(shù)和RC梁結(jié)構(gòu)抗力的關(guān)系;其次,通過引入失效模式特征系數(shù)k和銹蝕不平衡系數(shù)ψ,得到包含nm根銹蝕不均勻主筋的RC梁的抗力計算公式;最后基于隨機(jī)場理論,建立考慮材料性能和結(jié)構(gòu)尺寸空間變異性時RC梁的抗力時空退化模型.

    1 RC梁抗力時空退化模型

    1.1 鋼筋截面面積時變改進(jìn)模型

    在實際銹蝕中,由于材料和環(huán)境的不確定性,鋼筋很少單獨發(fā)生均勻銹蝕或者坑蝕,往往兩種銹蝕狀態(tài)同時發(fā)生,然而這兩種銹蝕狀態(tài)均會導(dǎo)致鋼筋截面面積的損失. 本文綜合考慮以上兩種銹蝕狀態(tài)的影響,得到截面面積計算公式為

    A(t)=M1Apit(t)+M2Auni(t).

    (1)

    式中:M1、M2分別為坑蝕和均勻銹蝕兩種銹蝕狀態(tài)發(fā)生的權(quán)重,M2=1-M1;Apit(t)和Auni(t)分別為鋼筋在坑蝕和均勻銹蝕下的剩余截面面積.

    運用文獻(xiàn)[11]中的面積損失計算公式得到鋼筋發(fā)生坑蝕時的截面面積損失,從而參數(shù)M1可由坑蝕造成的實際截面面積損失和理論截面面積損失的比值表示為

    (2)

    式中:Auni可近似于均勻銹蝕在鋼筋截面邊緣的非光滑區(qū)域引起的面積損失,Apit為除上述區(qū)域外由蝕坑影響引起的面積損失,bpit為蝕坑寬度,D0為鋼筋初始直徑.

    本文基于實驗所測數(shù)據(jù)及式(2)計算參數(shù)M1,部分結(jié)果如表1所示. 對參數(shù)M1進(jìn)行統(tǒng)計分析發(fā)現(xiàn):1)當(dāng)最大截面損失率C<9.5%時,M1服從均值為0.32、標(biāo)準(zhǔn)差為0.076的正態(tài)分布;2)當(dāng)C≥9.5%時,M1服從均值為0.65、標(biāo)準(zhǔn)差為0.078的正態(tài)分布. 對上述正態(tài)分布進(jìn)行K-S檢驗和S-W檢驗:當(dāng)C<9.5%時,K-S檢驗和S-W檢驗的結(jié)果分別為0.192和0.717;當(dāng)C≥9.5%時,K-S檢驗和S-W檢驗的結(jié)果分別為0.137和0.311. 以上結(jié)果表明:當(dāng)顯著性水平α=0.05時,參數(shù)M1服從正態(tài)分布.

    表1 參數(shù)M1(部分?jǐn)?shù)據(jù))

    1.2 銹蝕鋼筋屈服強(qiáng)度退化模型

    目前,相關(guān)研究表明鋼筋屈服強(qiáng)度fy的退化與截面面積損失可近似為線性相關(guān)[12-14],即

    (3)

    式中:ΔAst(t)為銹蝕引起的鋼筋截面面積損失,由式(1)推導(dǎo)可得;A0、fy0分別為鋼筋初始的截面面積和屈服強(qiáng)度;αy為經(jīng)驗系數(shù),由于實驗測試環(huán)境不同,國內(nèi)外學(xué)者對αy的取值存在不同的意見,文獻(xiàn)[15]認(rèn)為αy的值為0.017,而一些學(xué)者將αy取值為0.005,本文對實測的銹蝕鋼筋屈服強(qiáng)度進(jìn)行統(tǒng)計分析,得到αy的平均值為0.003 5.

    1.3 考慮裂縫的影響

    在服役過程中,混凝土結(jié)構(gòu)的開裂會導(dǎo)致混凝土強(qiáng)度降低,同時也會加快氯離子對結(jié)構(gòu)的侵蝕,導(dǎo)致結(jié)構(gòu)承載力下降. 而結(jié)構(gòu)承載力的下降會進(jìn)一步導(dǎo)致混凝土開裂,最終導(dǎo)致結(jié)構(gòu)承載力大幅降低.

    裂縫的存在為氯離子和其他物質(zhì)到達(dá)鋼筋提供了最優(yōu)路線,加速了銹蝕. 盡管氯離子擴(kuò)散特性受到裂縫寬度、溫度、濕度和表面氯離子質(zhì)量濃度變化等影響,但在實際結(jié)構(gòu)中,由于裂縫寬度較易獲得,且其與氯離子擴(kuò)散系數(shù)存在明顯的聯(lián)系,故本文只考慮裂縫寬度w對氯離子侵蝕的影響. 文獻(xiàn)[16]對服役8 a和11 a的鋼筋混凝土結(jié)構(gòu)早期銹脹裂縫研究,得到了不同裂縫寬度對應(yīng)的平均氯離子擴(kuò)散系數(shù)D(w)和表面氯離子質(zhì)量濃度,通過引入裂縫影響函數(shù)f(w)對開裂混凝土中氯離子的侵蝕行為描述為

    (4)

    式中D為非開裂混凝土中氯離子擴(kuò)散系數(shù),w≥0.1 mm.

    當(dāng)不考慮裂縫影響時,結(jié)構(gòu)服役過程中的混凝土強(qiáng)度fc′可認(rèn)為與初始混凝土強(qiáng)度fc相同;當(dāng)考慮裂縫影響時,混凝土強(qiáng)度fc′具有時變性. 本文基于銹蝕速度和裂縫寬度,采用文獻(xiàn)[17]建立的模型計算混凝土強(qiáng)度的降低值,從而得到任一時刻的混凝土強(qiáng)度fc′為

    (5)

    式中:Δfc為混凝土強(qiáng)度降低值,fc為初始混凝土強(qiáng)度,K是與鋼筋粗糙度和直徑相關(guān)的系數(shù),文獻(xiàn)[18]認(rèn)為當(dāng)結(jié)構(gòu)中的鋼筋為帶肋鋼筋時,K=0.1,ε1、εc0分別為開裂混凝土在直角方向上的壓應(yīng)力峰值和平均拉伸應(yīng)力,ε1的計算公式為

    (6)

    式中bf、b0分別為銹脹裂縫增加的寬度和初始狀態(tài)的截面寬度,bf與b0之差[19]為

    bf-b0=nbarsw=

    (7)

    式中:nbars為受壓鋼筋的數(shù)量,w為裂縫寬度,可由上式計算,νc為混凝土的泊松比,α為切向剛度折減系數(shù),具體計算方法見文獻(xiàn)[19],ft為混凝土的抗拉強(qiáng)度,Eef為混凝土的有效彈性模量,d1、d2分別為厚壁圓柱的內(nèi)、外半徑,且d1=(D0+2d0)/2,d2=S/2,D0為鋼筋初始直徑,d0為混凝土空隙的環(huán)形層厚度,S為鋼筋間距,ds(t)為腐蝕產(chǎn)物形成的厚度,計算公式[20]為

    (8)

    式中:αrust為與生銹類型相關(guān)的系數(shù);ρrust、ρst分別為鐵銹和鋼筋的密度;Wrust為單位長度的鋼筋銹蝕量,計算式[20]為

    (9)

    1.4 空間變異性

    為得到更為精準(zhǔn)的結(jié)構(gòu)抗力退化模型,本文將混凝土抗壓強(qiáng)度、保護(hù)層厚度、表面氯離子質(zhì)量濃度和臨界氯離子質(zhì)量濃度視為空間隨機(jī)變量,將鋼筋混凝土結(jié)構(gòu)離散成N個立體單元,運用一維隨機(jī)場理論對整個RC梁的抗力退化進(jìn)行研究. 隨機(jī)場的分析主要有以下4種方法:1)中點法;2)形函數(shù)法;3)空間均值法;4)Karhunen-Loeve分解法. 由于實際結(jié)構(gòu)中的材料性能、荷載和抗力等參數(shù)并非全部服從高斯分布,而中點法對于高斯分布和非高斯分布均適用,同時以任意單元質(zhì)心處的參數(shù)值代表該參數(shù)在該單元的分布值可大大簡化計算,故本文采用中點法進(jìn)行隨機(jī)場分析. 關(guān)于中點法的介紹見文獻(xiàn)[21].

    1.4.1 研究假定條件

    本文采用離散梁單元的方法對RC梁各參數(shù)的空間變異性進(jìn)行研究,所示算例的單元劃分如圖1所示. 基于劃分得到的N個梁單元作如下假定:1)對N個梁單元進(jìn)行研究時,采用的退化模型和材料參數(shù)的分布形態(tài)具有一致性; 2)各個梁單元已經(jīng)充分小,不再考慮各參數(shù)在單元內(nèi)部不同位置處的空間變異性.

    1.4.2 自相關(guān)函數(shù)

    隨機(jī)場分析的關(guān)鍵在于選取合適的自相關(guān)函數(shù)對各單元參數(shù)的空間變異性進(jìn)行描述. 兩相關(guān)單元的距離越小,自相關(guān)函數(shù)越大,反之越小. 目前,不同的隨機(jī)領(lǐng)域存在不同的自相關(guān)函數(shù),如三角函數(shù)、高斯函數(shù)、矩形函數(shù)等,而高斯(平方指數(shù))相關(guān)函數(shù)在工程實際研究中得到廣泛運用,因此,本文采用一維高斯函數(shù)作為自相關(guān)函數(shù),即

    (10)

    在分析材料性能和結(jié)構(gòu)尺寸的空間變異性時,一旦結(jié)構(gòu)單元劃分完畢,相關(guān)系數(shù)矩陣就隨之被確定. 假設(shè)結(jié)構(gòu)被劃分為N單元,則相關(guān)系數(shù)矩陣為

    (11)

    式中ρi j為單元i和單元j之間的相關(guān)系數(shù),有

    (12)

    1.5 銹蝕RC梁抗力時空退化模型

    考慮銹蝕RC梁各單元之間參數(shù)的空間變異性時,各單元中鋼筋的銹蝕速度、強(qiáng)度和材料性能退化值均不相同,從而每個單元中的極限承載能力也不盡相同. 不考慮空間變異性時,單筋RC梁的極限彎曲能力Mu計算式[14]為

    (13)

    式中:Q為模型不確定系數(shù),Ebond,i為第i個單元中鋼筋和混凝土之間的粘結(jié)行為,具體見文獻(xiàn)[14],b、ho分別為梁的寬度和有效高度,fc′為混凝土抗壓強(qiáng)度.

    然而在實際結(jié)構(gòu)中,RC梁底部一般包含多根(nm)受拉主筋,且由試驗分析可知,同一RC梁中各受拉主筋的銹蝕程度存在較大差異,則nm根鋼筋的抗力也隨之不同. 假設(shè)nm根鋼筋中銹蝕程度最小的鋼筋截面損失率為ρ0,將其余鋼筋的銹蝕程度用ρ0表示,則梁中nm根受拉主筋的銹蝕程度從小到大排列為

    {ρ0,ρ1,ρ2,…,ρm-1}={ρ0,B1ρ0,B2ρ0, …,Bm-1ρ0}.

    (14)

    同時定義極限損失量Mlimit和nm根受拉主筋的最大銹蝕率max{ρi}與ρ0的比值分別為失效模式特征系數(shù)k和銹蝕不均衡系數(shù)ψ,具體表示為

    (15)

    (16)

    則RC梁中nm根鋼筋的失效模式及結(jié)構(gòu)抗力可分為以下幾種情況.

    1)當(dāng)ρi≥Mlimit,即k≤1時,全部的受拉主筋均發(fā)生脆性失效,此時RC梁抗力為[14]

    r1(t)≤r2(t)≤rnm(t).

    (17)

    式中:r1(t)≤r2(t)≤……rnm(t),ri(t)=fy, i(t)Ast, i(t),fy, i(t)、Ast, i(t)分別為編號為i的鋼筋在時刻t時的剩余強(qiáng)度和截面面積.

    2)當(dāng)ρi

    (18)

    3)當(dāng)ρ0

    (19)

    式中L為集合{k≤Bi

    考慮結(jié)構(gòu)參數(shù)的空間變異性時,RC梁抗力的退化分析包含大量隨機(jī)參數(shù),如結(jié)構(gòu)材料、保護(hù)層厚度和銹蝕鋼筋的力學(xué)性能等,這些隨機(jī)變量均具有時間和空間變化的分布特征和退化過程,具體計算過程十分復(fù)雜且難以推導(dǎo)其閉合解. 因此,本文使用Monte-Carlo模擬的方法分析結(jié)構(gòu)抗力的退化. 只要隨機(jī)場單元確定,就能基于各單元之間的自相關(guān)函數(shù),得到各參數(shù)的空間隨機(jī)變量,繼而預(yù)測RC梁的抗力退化. 基于Monte-Carlo模擬計算RC梁抗力時空退化的步驟:1)基于隨機(jī)場理論,將結(jié)構(gòu)離散為多個空間單元,明確各單元內(nèi)所需考慮的隨機(jī)變量;2)選定合適的自相關(guān)函數(shù)模型,定義相關(guān)系數(shù)矩陣ρ; 3)將定義的系數(shù)矩陣ρ進(jìn)行Choleski分解得到一上三角矩陣S;4)計算考慮裂縫影響時任意時刻的混凝土強(qiáng)度值;5)基于上三角矩陣S及各單元初始隨機(jī)變量值得到所需考慮的各材料性能和結(jié)構(gòu)尺寸的空間隨機(jī)分布矩陣;6)綜合考慮所有隨機(jī)變量,對抗力退化模型進(jìn)行多重循環(huán)計算,得到各單元力學(xué)性能及抗力退化,得到RC梁抗力時空退化曲線.

    2 算例分析

    2.1 RC結(jié)構(gòu)

    選取一長度為12.5 m的RC簡支梁,其截面b×h=250 mm×500 mm,結(jié)構(gòu)底部為兩根Φ22鋼筋. 選擇單元尺寸時,若選取的單元個數(shù)越多,則計算結(jié)果越準(zhǔn)確,但此時運算量極大,導(dǎo)致程序計算速度大大降低. 綜合考慮以上因素,本文取單元尺寸Δ=0.5 m對RC梁進(jìn)行劃分,得到25個立體單元,此時RC梁中的主筋也隨之被離散,具體如圖1所示.

    圖1 梁單元劃分

    其中,用于結(jié)構(gòu)抗力計算的表面氯離子質(zhì)量濃度(%)、臨界氯離子質(zhì)量濃度(%)和混凝土保護(hù)層厚度(mm)均服從正態(tài)分布,均值和方差分別為0.114和0.1, 0.042和0.108, 及29.8和0.07;模型不確定系數(shù)Q和腐蝕電流密度(μA/cm2)均服從對數(shù)正態(tài)分布,均值和方差分別為0.98和0.5,及0.53和0.2;混凝土初始強(qiáng)度(MPa)服從均值為20.72和方差為0.177的正態(tài)分布;鋼筋初始屈服強(qiáng)度值為366 MPa. 同時,對計算裂縫影響的相關(guān)變量進(jìn)行如下取值:νc、ft、Eef、ρst、S分別取值為0.2 MPa、2.01 MPa、 3.0×104MPa、7 800 kg/m3、190 mm,d0、ρrust、αrust的取值分別為12.5 μm、3 600 kg/m3、0.57[19].

    本文假定各參數(shù)的變化均為平穩(wěn)隨機(jī)過程,考慮混凝土抗壓強(qiáng)度、保護(hù)層厚度、表面氯離子質(zhì)量濃度和臨界氯離子質(zhì)量濃度的空間變異性對結(jié)構(gòu)抗力的影響,并對結(jié)構(gòu)正常服役100 a后的抗力進(jìn)行預(yù)測. 在隨機(jī)場分析中,選取波動系數(shù)θ=2.0 m來計算各單元之間的相關(guān)函數(shù)值. 盡管文獻(xiàn)[8]研究發(fā)現(xiàn):當(dāng)波動系數(shù)由1 m變化到3 m時,發(fā)生銹蝕破壞的概率并沒有明顯的變化,但波動系數(shù)對結(jié)構(gòu)抗力的影響仍需進(jìn)一步展開研究.

    2.2 結(jié)果分析

    1)銹蝕模型的影響.當(dāng)結(jié)構(gòu)位于海洋線大氣環(huán)境下且考慮裂縫對抗力的影響時,基于改進(jìn)后的面積計算模型、點蝕模型和均勻銹蝕模型3種面積計算模型,研究了不同銹蝕模型對RC梁抗力退化規(guī)律的影響. 如圖2所示,不論采用哪種面積計算模型,RC梁的抗力均隨時間t的增大而減小. 當(dāng)不考慮參數(shù)的空間變異性時,結(jié)構(gòu)正常服役100 a后,采用本文改進(jìn)的面積計算模型得到的抗力值較采用均勻銹蝕模型時低了5.60%,比采用點蝕模型時高3.12%;同時采用改進(jìn)模型得到的抗力變化規(guī)律在服役前段更接近均勻銹蝕模型,而服役后段更靠近點蝕模型,這表明以往采用的點蝕模型低估了結(jié)構(gòu)的抗力,改進(jìn)的模型更好地反應(yīng)RC梁銹蝕狀態(tài). 當(dāng)考慮相關(guān)參數(shù)空間變異性時,以上3種面積計算模型對應(yīng)的RC梁的抗力值比不考慮各參數(shù)空間變異性時分別下降了14.68%、17.84%、14.13%,且初銹時間較不考慮各參數(shù)空間變異性時提高了9.28%. 因此不考慮各參數(shù)的空間變異性將低估RC梁的抗力退化.

    圖2 不同面積模型下抗力退化(Δ=0.5 m、θ=2.0 m)

    Fig.2 Resistance degradation under different area models(Δ=0.5 m,θ=2.0 m)

    2)裂縫的影響.圖3、4分別為裂縫與抗力退化關(guān)系圖和不同裂縫寬度w對應(yīng)的抗力退化圖. 由圖3可知,結(jié)構(gòu)正常服役100 a后,在考慮各參數(shù)空間變異性的前提下,不考慮裂縫影響得到的抗力值比考慮裂縫影響的抗力值高9.91%;而不考慮各參數(shù)的空間變異性時,不考慮裂縫影響得到的抗力值比考慮裂縫影響的抗力值高5.75%. 本文以沿海地區(qū)混凝土結(jié)構(gòu)在正常使用條件下的最大允許裂縫寬度0.1 mm為基礎(chǔ),研究不同裂縫寬度對應(yīng)的抗力退化:當(dāng)考慮空間變異性且裂縫寬度為0.1、0.15、0.2、0.25 mm時,結(jié)構(gòu)剩余抗力較不考慮裂縫影響時分別下降了2.41%、3.61%、4.60%、5.56%,且隨著裂縫寬度的增大,增加相同裂縫寬度值對結(jié)構(gòu)抗力的影響逐漸減小,鋼筋初始銹蝕時間的提前幅度也隨之減小. 表明裂縫對結(jié)構(gòu)正常服役行為影響較大,不考慮裂縫的影響會較大程度高估結(jié)構(gòu)的剩余抗力.

    3)隨機(jī)場參數(shù)的影響.由圖5可知,結(jié)構(gòu)正常服役100 a后,考慮結(jié)構(gòu)各參數(shù)的空間變異性且波動系數(shù)θ取1、2、3 m時,結(jié)構(gòu)抗力相比初始抗力分別下降51.57%、53.65%、55.47%. 與波動系數(shù)取2 m時相比,波動系數(shù)為1 m時結(jié)構(gòu)的初銹時間延緩了4.47%,而波動系數(shù)為3 m時結(jié)構(gòu)的初銹時間卻提前了4.91%. 這說明適當(dāng)?shù)脑黾硬▌酉禂?shù)能夠使結(jié)構(gòu)的抗力進(jìn)一步退化,但影響效果不顯著.

    當(dāng)單元尺寸Δ=12.5 m時,單元尺寸與結(jié)構(gòu)的實際尺寸相同,此時可視為不考慮各參數(shù)的空間變異性對抗力的影響. 由圖6可知,當(dāng)考慮各參數(shù)的空間變異性時,RC 梁的抗力隨著單元尺寸劃分的不同有較大的差異,其中在單元尺寸Δ=0.5、1.25 m兩種情形下,RC梁的抗力比不考慮空間變異性的影響時分別下降了14.68%、9.62%. 因此,當(dāng)考慮結(jié)構(gòu)各參數(shù)的空間變異性時,適當(dāng)?shù)倪x取單元尺寸對RC梁的安全評估十分重要.

    圖3 裂縫與抗力退化圖

    圖4 不同裂縫寬度w對應(yīng)的抗力

    圖5 波動系數(shù)θ與抗力退化

    4)環(huán)境參數(shù)的影響.在不同環(huán)境中,由于氯離子質(zhì)量濃度和腐蝕電流密度的不同,RC結(jié)構(gòu)的抗力退化規(guī)律也存在較大的差異. 圖7為RC結(jié)構(gòu)在海岸線大氣區(qū)、潮汐區(qū)和浪濺區(qū)3種環(huán)境下的抗力退化圖. 如圖7所示,不同環(huán)境對RC結(jié)構(gòu)抗力退化的影響顯著,在海岸線大氣區(qū)環(huán)境下,結(jié)構(gòu)正常服役100 a后,其剩余抗力值為初始結(jié)構(gòu)抗力值的46.35%,而在潮汐區(qū)和浪濺區(qū)環(huán)境下結(jié)構(gòu)剩余抗力分別為初始抗力值的43.23%和36.45%,比海岸線大氣區(qū)環(huán)境下的對應(yīng)值分別下降了6.73%和21.35%. 結(jié)果表明以上3種環(huán)境中,浪濺區(qū)環(huán)境對結(jié)構(gòu)抗力退化的影響最大,因此,在結(jié)構(gòu)服役期間,要盡量做好防腐措施,減少氯離子侵蝕對結(jié)構(gòu)耐久性的破壞.

    圖6 單元尺寸Δ與抗力退化

    圖7 不同環(huán)境下的抗力退化

    3 結(jié) 論

    1)以往計算鋼筋截面面積的點蝕模型低估了結(jié)構(gòu)的抗力,本文改進(jìn)后的面積時變模型能更好的反應(yīng)RC梁實際銹蝕狀態(tài).

    2)不考慮各參數(shù)的空間變異性和裂縫影響將較大程度上高估RC結(jié)構(gòu)的抗力:當(dāng)考慮各參數(shù)的空間變異性時,不考慮裂縫影響的抗力值比考慮裂縫影響時高9.91%,而不考慮各參數(shù)的空間變異性時,不考慮裂縫影響的抗力值比考慮裂縫影響時高5.75%,并且是否考慮各參數(shù)的空間變異性對預(yù)測RC梁的抗力值影響較大,兩種情況得到的抗力值最大相差高達(dá)14.68%.

    3)適當(dāng)?shù)卦黾硬▌酉禂?shù)能夠使結(jié)構(gòu)的抗力進(jìn)一步退化,但效果不顯著;選取適當(dāng)?shù)膯卧叽鐚C梁的安全評估十分重要;浪濺區(qū)環(huán)境對結(jié)構(gòu)抗力退化的影響最大,要盡量做好防腐措施,減少氯離子侵蝕對結(jié)構(gòu)耐久性的破壞.

    [1] 馬亞麗,張愛林. 基于規(guī)定可靠指標(biāo)的混凝土結(jié)構(gòu)氯離子侵蝕耐久壽命預(yù)測[J].土木工程學(xué)報,2006,39(2):36-41. DOI: 10.3321/j.issn:1000-131X.2006.02.008.

    MA Yali, ZHANG Ailin. Durablity life prediction of concrete structure based on the regulated reliability index under chloride environment[J]. China Civil Engineering Journal, 2006, 39(2):36-41. DOI:10.3321/j.issn:1000-131X.2006.02.008.

    [2] 劉文軍,王軍強(qiáng).氯離子對鋼筋混凝土結(jié)構(gòu)的侵蝕分析[J].混凝土,2007,27(4):20-22. DOI:10.3969/j.issn.1002-3550.2007.04.008.

    LIU Wenjun, WANG Junqiang. Analysis of chloride ion diffusion in reinforced concrete structures [J]. Concrete, 2007, 27(4):20-22. DOI:10.3969/j.issn.1002-3550.2007.04.008.

    [3] 孫叢濤.基于氯離子侵蝕的混凝土耐久性與壽命預(yù)測研究[D].西安:西安建筑科技大學(xué),2011.DOI:10.7666/d.d243557.

    SUN Congtao. Study on concrete durability and service life prodiction based on chloride corrosion [D]. Xi’an: Xi’an University of Architecture and Technology, 2011. DOI:10.7666/d.d243557.

    [4] 張俊芝,周巧萍,伍亞玲,等.基于混凝土氯離子擴(kuò)散系數(shù)時變性的鋼筋坑蝕面積隨機(jī)時變模型[J]. 自然災(zāi)害學(xué)報,2014, 23(2):254-263. DOI:10.13577/j.jnd.2014.0232.

    ZHANG Junzhi, ZHOU Qiaoping, WU Yaling, et al. A random time-varying model of corrosion pit area of steel bars based on the time variability of chloride ion diffusion coefficient in concrete[J]. Journal of Natural Disasters, 2014, 23(2):254-263. DOI:10.13577/j.jnd.2014.0232.

    [5] 王磊,馬亞飛,張建仁. 模糊及隨機(jī)信息下橋梁構(gòu)件抗力演化特征分析[J].土木工程學(xué)報, 2013, 46(1):76-83. DOI: 10.15951/j.tmgcxb.2013.01.005.

    WANG Lei, MA Yafei, ZHANG Jianren. Analysis for evolutionary characteristics of resistance degradation of bridge members under fuzzy and random information[J]. China Civil Engineering Journal,2013,46(1):76-83. DOI: 10.15951/j.tmgcxb.2013.01.005.

    [6] 吳鋒,張章,龔景海.基于銹脹裂縫的銹蝕梁鋼筋銹蝕率計算[J]. 建筑結(jié)構(gòu)學(xué)報,2013,34(10):144-150. DOI: 10.14006/j.jzjgxb.2013.10.017.

    WU Feng, ZHANG Zhang, GONG Jinghai. Calculation of steel corrosion rate based on corrosive crack of beams [J]. Journal of Building Structures, 2013,34(10):144-150. DOI: 10.14006/j.jzjgxb.2013.10.017.

    [7] 彭建新,胡守旺,張建仁,等.鋼筋混凝土結(jié)構(gòu)銹脹開裂寬度的試驗研究及預(yù)測模型[J]. 實驗力學(xué),2014, 29(1):33-41. DOI: 10.7520/1001-4888-13-003.

    PENG Jianxin, HU Shouwang, ZHANG Jianren, et al. Experimental study and prediction model of corrosion-induced crack width in RC structure[J].Journal of Experimental Mechanics, 2014,29(1):33-41. DOI: 10.7520/1001-4888-13-003.

    [8] STEWART M G,MULLARD J A.Spatial time-dependent reliability analysis of corrosion damage and the timing of first repair for RC structures[J].Engineering Structures, 2007, 29(7):1457-1464. DOI: 10.1016/j.engstruct. 2006.09.004.

    [9] VU K A, STEWART M G. Predicting the likelihood and extent of reinforced concrete corrosion-induced cracking [J]. Journal of Structural Engineering, 2005, 131(11): 1681-1689. DOI: 10.1061/(ASCE)0733-9445.

    [10]彭建新,胡守旺,張建仁.考慮溫室效應(yīng)的氯鹽環(huán)境下RC橋梁銹脹開裂性能預(yù)測[J].工程力學(xué),2013,30(8):103-110. DOI: 10.6052/j.issn.1000-4750.2012.04.0261.

    PENG Jianxin, HU Shouwang, ZHANG Jianren. Corrosion-induced crack performance prediction of rc bridge under chloride attack considering effect of global warming [J].Engineering Mechanics, 2013, 30(8): 103 -110. DOI: 10.6052/j.issn.1000-4750.2012.04.0261.

    [11]彭建新,陽逸鳴,張建仁. 銹蝕鋼筋截面面積分析及強(qiáng)度概率模型[J].中外公路,2015,35(5):124-129. DOI: 10.14048/j.issn.1671-2579.2015.05.028.

    PENG Jianxin, YANG Yiming, ZHANG Jianren. Probabilistic cross-sectional area and strength model of corroded steel bar [J]. Journal of China and Foreign Highway, 2015,35(5):124-129. DOI: 10.14048/j.issn.1671-2579. 2015.05.028.

    [12]STEWART M G. Mechanical behaviour of pitting corrosion of flexu-ral and shear reinforcement and its effect on structural reliability of corroding RC beams [J]. Structural Safety, 2009, 31(1): 19-30. DOI: 10.1016/j.strusafe.2007. 12.001.

    [13]MAYafei, ZHANG Jianren, WANG Lei, et al. Probabilistic prediction with Bayesian updating for strength degradation of RC bridge beams[J]. Structural Safety, 2013, 44(2334):102-109. DOI: 10.1016/j.strusafe.2013.07.006.

    [14]彭建新,張建仁,張克波,等. 銹蝕RC橋梁彎曲抗力時變概率模型與試驗研究[J]. 工程力學(xué),2012, 29(6):125-132. DOI: 10.6052/j.issn.1000-4750.2010.08.0614.

    PENG Jianxin, ZHANG Jianren, ZHANG Kebo, et al. Probabilistic resistance model and experimental investigation for corroded RC bridges[J]. Engineering Mechanics, 2012, 29(6): 125-132. DOI: 10.6052/j.issn.1000-4750.2010. 08.0614.

    [15]DU Y G, CLARK L A, CHAN A H C. Residual capacity of corroded reinforcing bars [J]. Magazine of Concrete Research, 2005, 57(3): 135-147. DOI: 10.1680/macr.2005.57.3.135.

    [16]KWON S J, NA U J, PARK S S, et al. Service life prediction of concrete wharves with early-aged crack: probabilistic approach for chloride diffusion[J]. Structure Safety,2009, 31(1):75-83. DOI: 10.1016/j.strusafe.2008.03.004.

    [17]VECCHIO F, COLLINS M P. The modified compression field theory for reinforced concrete elements subjected to shear [J]. Proc ACI Structures,1986, 83(2):219-231.

    [18]CAPE M. Residual service-life assessment of existing RC structures [D]. Goteborg:Chalmers University of Technology, 1999.

    [19]LI C Q, ZHENG J J, LAWANWISUT W, et al. Concrete delamination caused by steel reinforcement corrosion [J]. ASCE Mater Civil Eng, 2007,19(7):591-600. DOI: 10.1061/(ASCE)0899-1561(2007)19:7(591).

    [20]LIU Y, WEYERS R E. Modelling the time-to-corrosion cracking in chloride contaminated reinforced concrete structures [J]. ACI Mater J, 1998, 95(6):675-81.

    [21]VANMARCKE E H. Random field: analysis and synthesis [M]. London: The MIT Press, 1983.

    ResistancedegradationmodelofRCbeamundertheinfluenceofmanyfactors

    YANG Yiming, PENG Jianxin, ZHANG Jianren

    (Key Laboratory of Bridge Engineering Safety Control by Department of Education (Changsha University of Science and Technology), Hunan Province and Ministry of Education, Changsha 410114, China)

    To study the decline of RC structural resistance and durability caused by chloride ion erosion, based on RC beams accelerated corrosion test, firstly, the model of steel bar section area is analyzed considering the two kinds of uniform corrosion and pitting corrosion at the same time. Secondly, the relationship between the crack width and chloride diffusion coefficient and the resistance of RC beam is established. Then the structural resistance calculation model is discussed under the influence of corroded main bars with different corrosion degree in the same RC beam. Finally, considering spatial variability of material properties and structure size, the time-spatial-varying degradation model of RC beams resistance is built and the resistance of the bridge after 100 years of service is predicted. The results show that when the spatial variability of the parameters is considered, the resistance value wihtout considering the influence of the crack is 9.91% higher than the corresponding value considering the influence of crack. Considering the spatial variability of the parameters, the resistance of RC beam is decreased by 14.68% when compared with that without consideration. The appropriate choice of the fluctuation coefficient and the unit size to consider the spatial variability of the parameters of RC structure is very important for the safety evaluation of RC beam. The residual resistance of the structure in tidal area and splash zone is 43.23% and 36.45% of the initial value, compared with the atmospheric environment of the coastline, decreased by 6.73% and 21.35%, respectively. Therefore, in the service of the structure, corrosion protection measures should be done well to reduce the damage of the structure durability of chloride ion erosion.

    corroded beams; resistance degradation; cross-sectional area; crack; corrosion imbalance; spatial variability

    10.11918/j.issn.0367-6234.201611077

    TU375

    A

    0367-6234(2017)09-0058-07

    2016-11-17

    國家重點基礎(chǔ)研究發(fā)展計劃(2015CB057701);國家自然科學(xué)基金面上項目(51178060,51378081);湖南省教育廳優(yōu)秀青年項目(15B015);長沙理工大學(xué)交通基礎(chǔ)設(shè)施安全風(fēng)險管理行業(yè)重點實驗室開放基金(16BCX08)

    陽逸鳴(1991—),男,博士研究生;張建仁(1958—),男,教授,博士生導(dǎo)師

    張建仁,jianrenz@hotmail.com

    (編輯魏希柱)

    猜你喜歡
    抗力氯離子變異性
    橋梁低樁承臺基底豎向抗力作用效應(yīng)的計算研究
    交通科技(2021年4期)2021-09-03 09:47:22
    咳嗽變異性哮喘的預(yù)防和治療
    海洋水下區(qū)纖維混凝土中氯離子的擴(kuò)散性能
    低滲透黏土中氯離子彌散作用離心模擬相似性
    巖塊的彈性模量及巖體單位彈性抗力系數(shù)的確定方法
    水化硅酸鈣對氯離子的吸附
    30MnSi熱變形抗力的數(shù)學(xué)模型
    上海金屬(2014年2期)2014-12-18 06:52:51
    咳嗽變異性哮喘的中醫(yī)治療近況
    清肺止咳湯治療咳嗽變異性哮喘40例
    冬病夏治止咳貼貼敷治療小兒咳嗽變異性哮喘40例
    成人国产一区最新在线观看| 麻豆乱淫一区二区| 别揉我奶头~嗯~啊~动态视频| 三级毛片av免费| 久久久久精品国产欧美久久久| 国产亚洲欧美精品永久| 50天的宝宝边吃奶边哭怎么回事| 精品国产亚洲在线| 欧美日本中文国产一区发布| 亚洲人成电影观看| 大码成人一级视频| 国产老妇伦熟女老妇高清| 亚洲 国产 在线| 91精品三级在线观看| 18在线观看网站| 亚洲免费av在线视频| 亚洲中文日韩欧美视频| 亚洲伊人色综图| 国产午夜精品久久久久久| 两性夫妻黄色片| 欧美日韩亚洲国产一区二区在线观看 | 好男人电影高清在线观看| 午夜福利,免费看| 色精品久久人妻99蜜桃| 成人18禁高潮啪啪吃奶动态图| 老司机在亚洲福利影院| 在线观看免费视频网站a站| 精品久久久久久久毛片微露脸| 亚洲熟女精品中文字幕| 国产av国产精品国产| av不卡在线播放| 老熟女久久久| 一二三四社区在线视频社区8| 99riav亚洲国产免费| 高清在线国产一区| 激情视频va一区二区三区| 91麻豆精品激情在线观看国产 | 日韩视频一区二区在线观看| 精品国产一区二区久久| 香蕉久久夜色| 日本欧美视频一区| 1024视频免费在线观看| 国产片内射在线| 欧美日韩精品网址| 1024香蕉在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产亚洲av麻豆专区| 丝袜在线中文字幕| 成在线人永久免费视频| 91麻豆av在线| 国产精品98久久久久久宅男小说| 在线观看免费视频网站a站| 亚洲黑人精品在线| 高清在线国产一区| 麻豆成人av在线观看| 91精品三级在线观看| a级毛片在线看网站| 国产精品免费大片| 亚洲精品一卡2卡三卡4卡5卡| 久久精品熟女亚洲av麻豆精品| 18禁美女被吸乳视频| 18禁裸乳无遮挡动漫免费视频| 女人爽到高潮嗷嗷叫在线视频| 午夜免费成人在线视频| 久久精品人人爽人人爽视色| 久久久久久亚洲精品国产蜜桃av| 免费在线观看视频国产中文字幕亚洲| 精品人妻在线不人妻| 母亲3免费完整高清在线观看| 亚洲国产中文字幕在线视频| 国产日韩欧美视频二区| 91精品三级在线观看| 男人操女人黄网站| 亚洲av成人一区二区三| 少妇精品久久久久久久| 久久天躁狠狠躁夜夜2o2o| av一本久久久久| 亚洲欧美精品综合一区二区三区| 亚洲久久久国产精品| 香蕉国产在线看| 我的亚洲天堂| 丝瓜视频免费看黄片| 精品久久久久久久毛片微露脸| e午夜精品久久久久久久| 一区二区三区乱码不卡18| 日韩成人在线观看一区二区三区| 电影成人av| 久久ye,这里只有精品| 女同久久另类99精品国产91| 国产一卡二卡三卡精品| 国产成人av教育| 午夜91福利影院| 国产日韩欧美亚洲二区| 久久久久久久久久久久大奶| 淫妇啪啪啪对白视频| 亚洲国产欧美日韩在线播放| 99精国产麻豆久久婷婷| 亚洲精品国产一区二区精华液| 精品福利观看| 精品人妻熟女毛片av久久网站| 别揉我奶头~嗯~啊~动态视频| 一级毛片精品| 十八禁高潮呻吟视频| 久久人妻福利社区极品人妻图片| 欧美中文综合在线视频| 久久久久久久大尺度免费视频| 国产精品一区二区在线观看99| 女人精品久久久久毛片| 在线观看免费视频日本深夜| 男女床上黄色一级片免费看| 操美女的视频在线观看| 777久久人妻少妇嫩草av网站| 国产不卡av网站在线观看| 国产精品成人在线| 一级毛片精品| 一区福利在线观看| 波多野结衣一区麻豆| 国产成人免费无遮挡视频| 亚洲伊人色综图| 亚洲七黄色美女视频| 法律面前人人平等表现在哪些方面| 侵犯人妻中文字幕一二三四区| 精品乱码久久久久久99久播| 国产一区二区在线观看av| 国产在视频线精品| 久久天堂一区二区三区四区| 一区二区av电影网| 9热在线视频观看99| 1024香蕉在线观看| 女人精品久久久久毛片| 中文亚洲av片在线观看爽 | 男人操女人黄网站| 热99国产精品久久久久久7| 亚洲中文字幕日韩| 一本一本久久a久久精品综合妖精| 久9热在线精品视频| 9191精品国产免费久久| 两个人免费观看高清视频| 制服诱惑二区| 热99国产精品久久久久久7| 久久久欧美国产精品| av国产精品久久久久影院| 日本黄色视频三级网站网址 | 无人区码免费观看不卡 | 午夜福利在线观看吧| 人人澡人人妻人| 久久久久久久国产电影| 国产精品久久久久久人妻精品电影 | 三上悠亚av全集在线观看| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全电影3 | 国产又色又爽无遮挡免费看| 成人手机av| 亚洲中文日韩欧美视频| 少妇被粗大的猛进出69影院| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久成人av| av视频免费观看在线观看| 黄片大片在线免费观看| 久久国产亚洲av麻豆专区| 国产精品久久久av美女十八| 欧美日韩一级在线毛片| 欧美日韩av久久| 十八禁网站网址无遮挡| 法律面前人人平等表现在哪些方面| 亚洲精品中文字幕在线视频| 露出奶头的视频| 老司机福利观看| 丰满饥渴人妻一区二区三| 久久久国产精品麻豆| 国产单亲对白刺激| 另类精品久久| 国产福利在线免费观看视频| 人人澡人人妻人| 国产日韩欧美亚洲二区| 亚洲天堂av无毛| 香蕉久久夜色| 日本黄色日本黄色录像| 国产区一区二久久| 亚洲三区欧美一区| 天堂中文最新版在线下载| 亚洲欧洲精品一区二区精品久久久| 欧美日韩亚洲高清精品| 亚洲七黄色美女视频| 国产在线精品亚洲第一网站| 国产又爽黄色视频| 国产色视频综合| 精品国产国语对白av| 欧美变态另类bdsm刘玥| 欧美日韩一级在线毛片| 麻豆乱淫一区二区| 天天添夜夜摸| 黄色a级毛片大全视频| 精品亚洲成国产av| 波多野结衣一区麻豆| 国产成+人综合+亚洲专区| 久久亚洲精品不卡| 日韩有码中文字幕| 亚洲精品成人av观看孕妇| 国产三级黄色录像| 午夜久久久在线观看| 久久久久久久国产电影| 一区二区三区精品91| 久久精品国产综合久久久| 国产精品香港三级国产av潘金莲| 午夜老司机福利片| 建设人人有责人人尽责人人享有的| 久久国产精品大桥未久av| 免费av中文字幕在线| 亚洲全国av大片| 老鸭窝网址在线观看| 中文字幕av电影在线播放| 久久精品熟女亚洲av麻豆精品| 最黄视频免费看| 757午夜福利合集在线观看| 国产在线视频一区二区| 免费av中文字幕在线| 亚洲av成人不卡在线观看播放网| 两性夫妻黄色片| 在线十欧美十亚洲十日本专区| 欧美黑人欧美精品刺激| 一本—道久久a久久精品蜜桃钙片| 麻豆乱淫一区二区| 亚洲黑人精品在线| 中文字幕精品免费在线观看视频| 女性生殖器流出的白浆| 人人妻人人添人人爽欧美一区卜| 国产av一区二区精品久久| 国产欧美日韩一区二区精品| 久热爱精品视频在线9| 成人亚洲精品一区在线观看| 精品少妇久久久久久888优播| 啦啦啦免费观看视频1| 美女主播在线视频| 日韩中文字幕欧美一区二区| 纵有疾风起免费观看全集完整版| 亚洲欧美一区二区三区久久| 日韩欧美三级三区| 亚洲精品美女久久久久99蜜臀| www.自偷自拍.com| 亚洲欧美日韩高清在线视频 | 黄片小视频在线播放| 欧美在线黄色| 国产日韩欧美亚洲二区| avwww免费| 日韩精品免费视频一区二区三区| 欧美+亚洲+日韩+国产| 亚洲中文日韩欧美视频| 一级片'在线观看视频| aaaaa片日本免费| 久久久精品国产亚洲av高清涩受| 久久香蕉激情| 一个人免费在线观看的高清视频| 久久中文字幕人妻熟女| 交换朋友夫妻互换小说| 最新的欧美精品一区二区| 国产黄频视频在线观看| 国产极品粉嫩免费观看在线| 亚洲成a人片在线一区二区| 十分钟在线观看高清视频www| 亚洲色图综合在线观看| 两个人免费观看高清视频| 国产国语露脸激情在线看| 欧美日韩av久久| 美女高潮到喷水免费观看| 亚洲精品国产一区二区精华液| 亚洲三区欧美一区| 国产成人精品无人区| 日日爽夜夜爽网站| 久久中文看片网| 亚洲成a人片在线一区二区| 动漫黄色视频在线观看| 国产精品久久久久久人妻精品电影 | 亚洲va日本ⅴa欧美va伊人久久| 捣出白浆h1v1| 欧美久久黑人一区二区| h视频一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲成人免费电影在线观看| 国产成人精品久久二区二区91| 99精品在免费线老司机午夜| 人人妻人人爽人人添夜夜欢视频| 法律面前人人平等表现在哪些方面| 99久久人妻综合| 日韩大码丰满熟妇| 色婷婷久久久亚洲欧美| 国产一区二区三区在线臀色熟女 | 亚洲第一欧美日韩一区二区三区 | 1024视频免费在线观看| 一区二区三区乱码不卡18| 深夜精品福利| 国产片内射在线| 丁香六月天网| 热99re8久久精品国产| 亚洲熟女毛片儿| 性色av乱码一区二区三区2| 看免费av毛片| 精品乱码久久久久久99久播| 夫妻午夜视频| 天天影视国产精品| 99热网站在线观看| 亚洲专区中文字幕在线| 香蕉国产在线看| 国产精品国产高清国产av | 黄色片一级片一级黄色片| 一级黄色大片毛片| 久久久久久久精品吃奶| 又黄又粗又硬又大视频| 黄色视频,在线免费观看| 精品少妇久久久久久888优播| avwww免费| 亚洲免费av在线视频| 免费在线观看影片大全网站| 免费人妻精品一区二区三区视频| 欧美精品亚洲一区二区| 久久99一区二区三区| 丝袜美足系列| 国产在线免费精品| 黄色怎么调成土黄色| 久久天堂一区二区三区四区| 熟女少妇亚洲综合色aaa.| 国产精品影院久久| 操出白浆在线播放| 婷婷成人精品国产| 日本av免费视频播放| 水蜜桃什么品种好| 国产成人精品久久二区二区91| 欧美日韩国产mv在线观看视频| 精品人妻熟女毛片av久久网站| 国产人伦9x9x在线观看| 精品人妻在线不人妻| 一本—道久久a久久精品蜜桃钙片| 少妇粗大呻吟视频| 制服诱惑二区| 精品国产亚洲在线| 亚洲人成电影观看| 日韩中文字幕欧美一区二区| 午夜91福利影院| 久久av网站| 国产精品欧美亚洲77777| 免费黄频网站在线观看国产| aaaaa片日本免费| 18禁观看日本| 最黄视频免费看| 建设人人有责人人尽责人人享有的| 99国产精品一区二区蜜桃av | 91成人精品电影| 免费在线观看黄色视频的| 757午夜福利合集在线观看| 一个人免费看片子| 日本精品一区二区三区蜜桃| 亚洲男人天堂网一区| 肉色欧美久久久久久久蜜桃| 麻豆av在线久日| 亚洲国产av新网站| 久久久久国产一级毛片高清牌| 老司机福利观看| 在线观看免费高清a一片| 欧美激情久久久久久爽电影 | 国产男女内射视频| 国产高清videossex| 两个人看的免费小视频| 啦啦啦中文免费视频观看日本| 啦啦啦在线免费观看视频4| 亚洲欧美日韩另类电影网站| 亚洲av成人不卡在线观看播放网| 他把我摸到了高潮在线观看 | 亚洲成国产人片在线观看| 高潮久久久久久久久久久不卡| 九色亚洲精品在线播放| 一级黄色大片毛片| 亚洲成人国产一区在线观看| 纵有疾风起免费观看全集完整版| av超薄肉色丝袜交足视频| 18在线观看网站| 国产欧美日韩综合在线一区二区| 女警被强在线播放| 一级片'在线观看视频| 精品高清国产在线一区| 欧美日韩av久久| 老熟女久久久| 人妻 亚洲 视频| 国产精品九九99| 麻豆av在线久日| 亚洲第一欧美日韩一区二区三区 | 精品一区二区三区视频在线观看免费 | 高清视频免费观看一区二区| 亚洲五月婷婷丁香| 黄片播放在线免费| 超碰97精品在线观看| 美女午夜性视频免费| 黄色视频不卡| 视频在线观看一区二区三区| 国产麻豆69| 最近最新中文字幕大全电影3 | 在线十欧美十亚洲十日本专区| 久久国产精品影院| 欧美 日韩 精品 国产| 午夜福利欧美成人| 69av精品久久久久久 | 天堂俺去俺来也www色官网| 久久狼人影院| 欧美激情极品国产一区二区三区| 免费不卡黄色视频| 精品少妇内射三级| 每晚都被弄得嗷嗷叫到高潮| 80岁老熟妇乱子伦牲交| 99国产精品免费福利视频| 国产福利在线免费观看视频| 久久香蕉激情| 高潮久久久久久久久久久不卡| 在线 av 中文字幕| 国产男靠女视频免费网站| 啦啦啦 在线观看视频| 久久久久国内视频| 午夜久久久在线观看| 午夜91福利影院| 国产精品久久久av美女十八| av超薄肉色丝袜交足视频| 亚洲情色 制服丝袜| 在线亚洲精品国产二区图片欧美| 俄罗斯特黄特色一大片| 大码成人一级视频| 精品少妇一区二区三区视频日本电影| 欧美在线黄色| 国产伦人伦偷精品视频| 亚洲伊人久久精品综合| 国产成人av激情在线播放| 在线十欧美十亚洲十日本专区| 久久人人97超碰香蕉20202| 宅男免费午夜| 久久亚洲真实| 一本—道久久a久久精品蜜桃钙片| 窝窝影院91人妻| 精品少妇一区二区三区视频日本电影| 天天躁狠狠躁夜夜躁狠狠躁| 一个人免费在线观看的高清视频| 久久久久久亚洲精品国产蜜桃av| 热99国产精品久久久久久7| tube8黄色片| 后天国语完整版免费观看| 色尼玛亚洲综合影院| 亚洲精品一卡2卡三卡4卡5卡| 19禁男女啪啪无遮挡网站| 亚洲精品乱久久久久久| 真人做人爱边吃奶动态| 建设人人有责人人尽责人人享有的| 涩涩av久久男人的天堂| 久久久久精品人妻al黑| 欧美日韩成人在线一区二区| 亚洲精品久久午夜乱码| 精品高清国产在线一区| 午夜免费成人在线视频| 黄色视频,在线免费观看| 日本av手机在线免费观看| 侵犯人妻中文字幕一二三四区| 久久亚洲精品不卡| 亚洲一码二码三码区别大吗| 精品久久蜜臀av无| 成年人午夜在线观看视频| 日韩 欧美 亚洲 中文字幕| 国产亚洲精品久久久久5区| 法律面前人人平等表现在哪些方面| 妹子高潮喷水视频| 国产午夜精品久久久久久| 亚洲免费av在线视频| 亚洲欧美精品综合一区二区三区| 黄频高清免费视频| 精品久久久精品久久久| 日韩视频一区二区在线观看| 国产免费视频播放在线视频| 99国产精品免费福利视频| 人人妻人人添人人爽欧美一区卜| 国产精品久久久av美女十八| 国产成人一区二区三区免费视频网站| 亚洲国产欧美日韩在线播放| 一区二区三区精品91| 久热爱精品视频在线9| 岛国毛片在线播放| 欧美国产精品一级二级三级| 国产日韩欧美亚洲二区| 久久午夜亚洲精品久久| 成人永久免费在线观看视频 | 国产精品久久电影中文字幕 | 真人做人爱边吃奶动态| 午夜福利乱码中文字幕| 99香蕉大伊视频| 在线 av 中文字幕| 亚洲精品美女久久av网站| 欧美成人免费av一区二区三区 | 亚洲国产欧美网| 日韩欧美一区二区三区在线观看 | 欧美精品一区二区免费开放| 侵犯人妻中文字幕一二三四区| 免费看十八禁软件| 国产视频一区二区在线看| 免费看十八禁软件| 亚洲成人免费电影在线观看| 久久中文字幕一级| 青青草视频在线视频观看| 亚洲欧美色中文字幕在线| 久久午夜亚洲精品久久| 女性被躁到高潮视频| 中文字幕人妻丝袜一区二区| 中文字幕色久视频| 亚洲三区欧美一区| 狠狠婷婷综合久久久久久88av| 国产精品一区二区免费欧美| 乱人伦中国视频| 免费看十八禁软件| 亚洲五月色婷婷综合| 99riav亚洲国产免费| 激情视频va一区二区三区| 亚洲五月婷婷丁香| 999精品在线视频| 69精品国产乱码久久久| 一二三四在线观看免费中文在| 久久99一区二区三区| 色老头精品视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 又黄又粗又硬又大视频| videosex国产| 50天的宝宝边吃奶边哭怎么回事| 两个人免费观看高清视频| 亚洲成国产人片在线观看| 十八禁高潮呻吟视频| 久久九九热精品免费| 一本久久精品| 国产亚洲精品第一综合不卡| a级毛片在线看网站| 搡老熟女国产l中国老女人| 精品乱码久久久久久99久播| 久久久久国产一级毛片高清牌| 97人妻天天添夜夜摸| 欧美激情极品国产一区二区三区| 久久久久精品国产欧美久久久| 侵犯人妻中文字幕一二三四区| 男女下面插进去视频免费观看| 国产精品1区2区在线观看. | 国产精品自产拍在线观看55亚洲 | 久久婷婷成人综合色麻豆| 国产精品美女特级片免费视频播放器 | 满18在线观看网站| 免费一级毛片在线播放高清视频 | 色94色欧美一区二区| 淫妇啪啪啪对白视频| 精品人妻在线不人妻| 婷婷丁香在线五月| 免费看a级黄色片| 国产成人av教育| 久久久久久久大尺度免费视频| 在线观看免费午夜福利视频| 亚洲精品在线观看二区| 一级片'在线观看视频| 日韩欧美一区二区三区在线观看 | 亚洲人成伊人成综合网2020| www.自偷自拍.com| 亚洲伊人久久精品综合| 丝袜喷水一区| 少妇的丰满在线观看| 国产亚洲精品第一综合不卡| 女人精品久久久久毛片| 午夜免费成人在线视频| 日本vs欧美在线观看视频| 狠狠精品人妻久久久久久综合| 超色免费av| 丰满饥渴人妻一区二区三| a在线观看视频网站| 亚洲色图av天堂| 男女之事视频高清在线观看| netflix在线观看网站| 女人久久www免费人成看片| 国产av国产精品国产| 在线天堂中文资源库| 色视频在线一区二区三区| 午夜福利,免费看| 国产aⅴ精品一区二区三区波| 91麻豆精品激情在线观看国产 | 欧美日本中文国产一区发布| 怎么达到女性高潮| 人人妻人人澡人人爽人人夜夜| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美国产精品一级二级三级| 久久久国产成人免费| av线在线观看网站| av网站在线播放免费| 91麻豆精品激情在线观看国产 | 亚洲色图 男人天堂 中文字幕| 91字幕亚洲| 久久精品国产99精品国产亚洲性色 | 热re99久久国产66热| 视频区欧美日本亚洲| 中文字幕精品免费在线观看视频| 交换朋友夫妻互换小说| 亚洲三区欧美一区| 99热国产这里只有精品6| 日本五十路高清| 亚洲国产欧美在线一区| 国产av一区二区精品久久| 国产精品一区二区免费欧美| 午夜福利免费观看在线| www日本在线高清视频| 亚洲人成电影观看| 好男人电影高清在线观看| 欧美一级毛片孕妇| 国产精品国产av在线观看| 日日夜夜操网爽| 精品免费久久久久久久清纯 | 国产欧美日韩一区二区三区在线| 一本色道久久久久久精品综合| 考比视频在线观看| 日日爽夜夜爽网站| 男女免费视频国产| 久久久久久免费高清国产稀缺| 三级毛片av免费|