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

    氯離子在混凝土中擴(kuò)散效應(yīng)分析的MCA模型

    2023-09-11 03:12:36馬俊軍藺鵬臻
    鐵道學(xué)報 2023年8期
    關(guān)鍵詞:元胞細(xì)觀擴(kuò)散系數(shù)

    馬俊軍,藺鵬臻

    (蘭州交通大學(xué) 土木工程學(xué)院,甘肅 蘭州 730070)

    對于海洋等腐蝕性環(huán)境下的混凝土結(jié)構(gòu)而言,由于氯離子侵蝕引起的鋼筋銹蝕是導(dǎo)致混凝土結(jié)構(gòu)發(fā)生耐久性失效的主要原因[1-2]。研究氯離子在混凝土中的擴(kuò)散機(jī)理和材料組分對混凝土中氯離子侵蝕擴(kuò)散的影響規(guī)律將對海洋等腐蝕性環(huán)境下混凝土結(jié)構(gòu)的耐久性設(shè)計、維修和加固具有重要意義[3]。

    目前,物理試驗由于與實際情況最為接近,已成為國內(nèi)外學(xué)者研究混凝土中氯離子擴(kuò)散性能最常用的方法之一[4]。但由于試驗設(shè)備、場地、條件等因素的限制,使得試驗結(jié)果無法反映真實環(huán)境中結(jié)構(gòu)的退化侵蝕現(xiàn)象[5]。而且試驗研究還存在周期長、成本高、試驗結(jié)果離散性大等問題。針對試驗方法的局限性,文獻(xiàn)[6-8]在Fick第二定律的基礎(chǔ)上,通過考慮不同的影響因素提出一系列用于描述混凝土中氯離子擴(kuò)散過程的預(yù)測模型。雖然這些模型能夠很好地預(yù)測混凝土中氯離子的擴(kuò)散行為,但都將混凝土看成一種宏觀勻質(zhì)材料,且模型采用的理論公式比較簡單,只能考慮單一、確定性因素對氯離子擴(kuò)散效應(yīng)的影響,無法考慮材料組分和隨機(jī)分布對氯離子擴(kuò)散效應(yīng)的影響[9]。為此,文獻(xiàn)[10]從細(xì)觀角度出發(fā),將混凝土看成由砂漿、界面過渡區(qū)(Interfacial Transition Zone,ITZ)、骨料組成的三相復(fù)合材料,并利用有限元法建立用于描述氯離子在混凝土中的細(xì)觀模型,并對混凝土氯離子擴(kuò)散過程的影響因素進(jìn)行分析。文獻(xiàn)[11]利用有限元法研究骨料對再生骨料混凝土中氯離子擴(kuò)散特性的影響。文獻(xiàn)[12-13]將混凝土作為由砂漿和骨料組成的兩相復(fù)合材料,利用有限元法對混凝土中氯離子擴(kuò)散行為進(jìn)行分析。文獻(xiàn)[2,4]提出一種考慮材料非均質(zhì)分布的細(xì)觀化等效方法,有效解決細(xì)觀材料的宏觀描述,提高模型的計算效率。但這些基于細(xì)觀尺度的混凝土氯離子擴(kuò)散分析方法,由于網(wǎng)格尺寸的限制,無法考慮混凝土實際界面過渡區(qū)的尺寸,只適用于小尺度材料層面的研究,無法在整體結(jié)構(gòu)中適用。

    自文獻(xiàn)[14]用元胞自動機(jī)(Cellular Automata,CA)模擬氯離子在混凝土中的擴(kuò)散過程以來,大量學(xué)者[3,15-16]利用元胞自動機(jī)對混凝土氯離子擴(kuò)散過程進(jìn)行研究。與有限元方法相比,該方法對網(wǎng)格尺寸沒有限制,可以根據(jù)需要隨意劃分網(wǎng)格。因此,結(jié)合細(xì)觀等效化方法和CA求解方法的優(yōu)點(diǎn),用來解決細(xì)觀模型在應(yīng)用上的局限性。近來,文獻(xiàn)[3]采用細(xì)觀數(shù)值模擬方法,建立用于描述氯離子在混凝土中擴(kuò)散行為的CA模型,并對元胞尺寸大小、時間步長等求解細(xì)節(jié)進(jìn)行討論和優(yōu)化。

    在上述研究的基礎(chǔ)上,本文利用元胞自動機(jī)和細(xì)觀尺度模擬相結(jié)合的方法對混凝土中氯離子擴(kuò)散行為進(jìn)行研究。利用細(xì)觀尺度模擬方法建立考慮骨料等不均勻分布的混凝土隨機(jī)細(xì)觀模型。根據(jù)氯離子在混凝土中的實際擴(kuò)散特點(diǎn),推導(dǎo)基于細(xì)觀尺度的混凝土氯離子CA演變方程,建立基于細(xì)觀尺度的元胞自動機(jī)(Mesoscale Cellular Automata,MCA)模型。利用試驗結(jié)果對所提出的MCA模型的正確性和描述材料變異性的能力進(jìn)行驗證。利用提出的模型對影響混凝土氯離子擴(kuò)散特性的各個細(xì)觀結(jié)構(gòu)參數(shù)進(jìn)行分析。

    1 混凝土細(xì)觀模型的建立

    在細(xì)觀尺度上,可將混凝土看成由砂漿、界面過渡區(qū)(ITZ)、骨料組成的三相復(fù)合材料[17-18],由于骨料形狀基本不會對混凝土中氯離子擴(kuò)散效應(yīng)產(chǎn)生影響[13,19-20],因此本文骨料采用圓形。骨料粒徑范圍為4.75~31.50 mm,滿足Fuller級配曲線?;诿商乜?Monte-Carlo,MC)方法對生成的骨料顆粒進(jìn)行隨機(jī)投放[18],建立考慮骨料等不均勻分布的混凝土隨機(jī)細(xì)觀模型,見圖1。圖1中,ITZ被模擬成介于骨料和砂漿之間的等厚度薄層,其厚度一般在20~100 μm之間[21],為簡化研究和提高計算效率,本文取ITZ厚度為50 μm。

    圖1 混凝土細(xì)觀模型示意

    2 基于細(xì)觀尺度的混凝土氯離子擴(kuò)散過程的元胞自動機(jī)表達(dá)

    元胞自動機(jī)作為一種數(shù)值計算方法,可以用來求解氯離子在混凝土中的擴(kuò)散行為。在已知元胞領(lǐng)域類型和局部進(jìn)化規(guī)則的情況下,元胞自動機(jī)局部演化方程可描述為[14-15]

    S(i,t+Δt)=f[S(1,t),…,S(i,t),…,S(n,t)]

    ( 1 )

    式中:S(i,t+ Δt)為t+Δt時刻,元胞i的狀態(tài);f為狀態(tài)轉(zhuǎn)移函數(shù);n為元胞總數(shù)。

    元胞C(i,j,k)各表面氯離子擴(kuò)散通量示意見圖2。圖2中,t+ Δt時刻元胞C(i,j,k)中氯離子含量與t時刻元胞C(i,j,k)中氯離子含量和氯離子擴(kuò)散通量之間的關(guān)系可描述為[3]

    圖2 元胞C(i, j, k)各表面氯離子擴(kuò)散通量示意

    ( 2 )

    式中:C(i,j,k,t)、C(i,j,k,t+ Δt)分別為t、t+ Δt時刻單元(i,j)中氯離子含量:k為混凝土中氯離子擴(kuò)散沿平面坐標(biāo)的4個方向;Ak為兩元胞之間沿k方向氯離子擴(kuò)散面的面積,m2;Jout為從元胞(i,j)流出的氯離子擴(kuò)散通量;Jin為從元胞(i,j)流入的氯離子擴(kuò)散通量。

    二維元胞領(lǐng)域類型示意見圖3。

    圖3 二維元胞領(lǐng)域類型示意

    對于平面模型,根據(jù)式( 1 )、圖3,式( 2 )可進(jìn)一步改寫為

    ( 3 )

    根據(jù)擴(kuò)散第一定律,沿某一方向,氯離子擴(kuò)散通量與其濃度成正比,即

    ( 4 )

    式中:Dk(i,j)為元胞(i,j)沿k方向擴(kuò)散到鄰居元胞(i+ik,j+jk)時的氯離子擴(kuò)散系數(shù),m2/s;c(i,j,t)為t時刻元胞(i,j)氯離子濃度值,%;δk為兩相鄰元胞之間的中心距離,m;當(dāng)k=1~4時,ik、jk的取值分別為

    ( 5 )

    結(jié)合式( 4 ),式( 3 )可以進(jìn)一步描述為

    ( 6 )

    當(dāng)采用正方形元胞時,即dk=δk=δ。其中,δ為正方形元胞尺寸。式( 6 )可以簡化為

    [c(i,j,t)-c(i+ik,j+jk,t)}

    ( 7 )

    根據(jù)大數(shù)定理,元胞內(nèi)氯離子含量與氯離子濃度成正比,即

    C(i,j,k,t)=c(i,j,k,t)·Ae

    ( 8 )

    式中:Ae為元胞內(nèi)氯離子傳輸?shù)挠行娣e。

    此時,氯離子在混凝土中的擴(kuò)散過程可以表達(dá)為

    [c(i,j,t)-c(i+ik,j+jk,t)}

    ( 9 )

    在不考慮骨料傳輸?shù)那闆r下,即骨料中氯離子擴(kuò)散系數(shù)DAgg=0,有效面積Ae可以表示為

    Ae=(AMor+AITZ)δ2=(1-AAgg)δ2

    (10)

    式中:AMor、AITZ分別為元胞內(nèi)砂漿、TTZ面積;AAgg為元胞內(nèi)骨料所占的面積。

    至此,式( 9 )可表達(dá)為

    (11)

    利用式(11)可建立混凝土中氯離子擴(kuò)散過程的元胞自動機(jī)模型。

    對于混凝土細(xì)觀模型而言,由于骨料和界面過渡區(qū)在空間上的隨機(jī)分布,使得元胞(i,j)沿任意方向表現(xiàn)出不同的侵蝕特性,主要體現(xiàn)在氯離子擴(kuò)散系數(shù)Dk(i,j) 在各個方向上的取值具有隨機(jī)性。因此,為便于分析這種材料隨機(jī)性,本文提出一種細(xì)觀等效方法來模擬材料隨機(jī)性,以圖1中的混凝土細(xì)觀模型中的元胞(i,j)為例,其等效過程見圖4。主要分為以下幾個步驟:

    圖4 細(xì)觀模型等效過程示意

    Step1根據(jù)元胞內(nèi)砂漿、ITZ的面積分?jǐn)?shù),將元胞內(nèi)砂漿、ITZ的氯離子擴(kuò)散系數(shù)DMor(i,j)、DITZ(i,j)進(jìn)行第一次均勻化等效,得到等效后氯離子擴(kuò)散系數(shù)DMI(i,j)。

    Step2根據(jù)元胞內(nèi)骨料的面積分?jǐn)?shù),將Step1中等效后的氯離子擴(kuò)散系數(shù)DMI(i,j)和骨料中氯離子擴(kuò)散系數(shù)DAgg(i,j)進(jìn)行第二次等效,得到等效后元胞氯離子擴(kuò)散系數(shù)D(i,j)。

    Step3根據(jù)元胞自動機(jī)局部演變規(guī)則,將上述得到的各元胞氯離子擴(kuò)散系數(shù)D(i,j)沿著不同的擴(kuò)散方向進(jìn)行平均化處理,得到平均氯離子擴(kuò)散系數(shù)Dk(i,j)。

    由圖4可知,等效后各元胞的氯離子擴(kuò)散系數(shù)D(i,j) 可以表示為各混凝土組成成分與其氯離子擴(kuò)散系數(shù)的關(guān)系式,即

    (12)

    式中:AITZ(i,j)、AMor(i,j)、AAgg(i,j)分別為元胞(i,j)中界面過渡區(qū)、砂漿、骨料所占的面積分?jǐn)?shù),%;DAgg為骨料中氯離子擴(kuò)散系數(shù),一般取DAgg=0;DMor為砂漿中氯離子擴(kuò)散系數(shù);DITZ為ITZ中氯離子擴(kuò)散系數(shù),其與砂漿中氯離子擴(kuò)散系數(shù)DITZ存在如下關(guān)系[21-22]

    DITZ=αDMor

    (13)

    其中,α為ITZ與砂漿中氯離子擴(kuò)散系數(shù)的比值,根據(jù)文獻(xiàn)[22],α和界面過渡區(qū)厚度tITZ之間的關(guān)系可以表示為

    (14)

    根據(jù)元胞內(nèi)有效擴(kuò)散系數(shù)D(i,j),Step3中各方向等效擴(kuò)散系數(shù)Dk(i,j)表現(xiàn)出不同的等效結(jié)果,分別為:

    (1)當(dāng)D(i,j)=0或D(i+ik,j+jk)=0時,k方向擴(kuò)散系數(shù)Dk(i,j)可以等效為

    Dk(i,j)=0

    (15)

    (2)當(dāng)D(i,j)≠0且D(i+ik,j+jk)≠0時,k方向擴(kuò)散系數(shù)Dk(i,j) 可以等效為

    (16)

    3 氯離子擴(kuò)散過程的時變效應(yīng)

    根據(jù)Thomas等[23]的試驗研究,發(fā)現(xiàn)混凝土中氯離子擴(kuò)散系數(shù)隨時間t的關(guān)系可以用冪函數(shù)來表示,即

    (17)

    式中:DMor(t)為考慮侵蝕時間后的氯離子擴(kuò)散系數(shù),m2/s;D0為表觀氯離子擴(kuò)散系數(shù),m2/s;t0為混凝土參考齡期,一般取t0=28 d;m為時間依賴性參數(shù)。

    4 計算流程

    按照上述提出的混凝土隨機(jī)細(xì)觀模型和元胞自動機(jī)模擬方法,利用Matlab軟件編寫混凝土氯離子擴(kuò)散效應(yīng)分析的細(xì)觀元胞自動機(jī)模型計算程序,具體流程見圖5。圖5中,Ncell為元胞總個數(shù)。

    圖5 混凝土氯離子擴(kuò)散效應(yīng)分析流程示意

    5 元胞尺寸確定方法

    從上述等效過程可知,元胞尺寸太大,模型無法真實反映骨料在截面內(nèi)的實際分布,造成模擬值與實際結(jié)果的偏差。元胞尺寸太小,不利于混凝土細(xì)觀參數(shù)的統(tǒng)計與提取,也降低了計算效率。根據(jù)式(10),在采用圓形骨料的情況下,為使建立的MCA模型能真實反映骨料在截面內(nèi)的實際分布,理論上元胞尺寸應(yīng)滿足如下關(guān)系

    (18)

    式中:δe、δmax分別為模型計算、最大理論尺寸,當(dāng)δe=δmax時,模型計算結(jié)果最準(zhǔn)確、高效,當(dāng)采用尺寸δe無法進(jìn)行網(wǎng)格劃分時,可取δe<δmax,但為保證精度,取值不易太小;dmin為骨料最小粒徑,m。

    為驗證式(18)確定的元胞尺寸是否合理。以圖1中的混凝土細(xì)觀模型為例,考慮混凝土實際尺寸和元胞類型,同時為便于網(wǎng)格劃分,利用數(shù)值試驗,分別比較δ=1、2、4、5 mm時MCA模型的模擬結(jié)果。模擬歷經(jīng)480 d后混凝土截面氯離子濃度的分布情況,選擇距離混凝土表面不同深度氯離子濃度模擬值的平均值作為評價指標(biāo),比較結(jié)果見圖6。同時,為便于比較,圖中還給出了宏觀尺度下Fick第二定律的解析解。

    圖6 不同元胞尺寸下MCA模型的模擬結(jié)果與解析解的比較

    由圖6可知,不同網(wǎng)格尺寸下,截面氯離子濃度模擬值隨深度的變化規(guī)律與解析解隨深度的變化規(guī)律一致。當(dāng)δ=2 mm(最接近δmax=3.4 mm)時,模型模擬值與解析解之間偏差最小。當(dāng)元胞尺寸太小(δ=1 mm),模型將低估氯離子在混凝土中的擴(kuò)散行為。當(dāng)δ>δmax,即δ=4、5 mm,模型將高估氯離子在混凝土中的擴(kuò)散行為。這說明MCA模型元胞尺寸確定方法是合理的。對于常規(guī)混凝土而言,為使MCA模型計算結(jié)果高效、準(zhǔn)確,推薦采用δ=0.707dmin。

    6 數(shù)值模擬與試驗驗證

    6.1 MCA模型的建立

    為驗證MCA模型的準(zhǔn)確性,對文獻(xiàn)[24]中強(qiáng)度等級為C30的混凝土受氯離子侵蝕試驗過程進(jìn)行數(shù)值模擬,將模型模擬結(jié)果與試驗結(jié)果進(jìn)行比較。根據(jù)混凝土細(xì)觀尺度模擬方法,建立截面尺寸為100 mm×100 mm的混凝土細(xì)觀模型,見圖1。圖1中,骨料面積分?jǐn)?shù)(AAgg=40%)和粒徑分布(4.75~31.50 mm)與文獻(xiàn)[24]相同。為便于與試驗結(jié)果進(jìn)行對比,根據(jù)第4節(jié)的研究結(jié)果,在保證計算精度的情況下,δ=2 mm,劃分后的混凝土隨機(jī)細(xì)觀模型見圖7。

    圖7 混凝土細(xì)觀模型網(wǎng)格劃分示意

    在混凝土細(xì)觀模型等效時,D0=3.931×10-12m2/s,m=0.608 3,t=480 d,混凝土界面擴(kuò)散特性參數(shù)利用式(14)確定,等效后各元胞等效氯離子擴(kuò)散系數(shù)見圖8。由圖8可知,等效后截面各位置氯離子擴(kuò)散系數(shù)不再是常數(shù),而是隨骨料等組成成分的變化而變化。通過與圖1對比,發(fā)現(xiàn)顏色越接近紅色代表元胞內(nèi)砂漿的面積分?jǐn)?shù)越高,等效后D0越大;顏色越接近藍(lán)色代表元胞內(nèi)AITZ越高,等效后D0越小。說明提出的氯離子擴(kuò)散系數(shù)等效方法可以模擬混凝土中骨料和ITZ隨機(jī)分布對D0的影響。

    圖8 等效氯離子擴(kuò)散系數(shù)分布

    6.2 氯離子濃度驗證

    為驗證MCA模型的準(zhǔn)確性,將模型的模擬結(jié)果與文獻(xiàn)[24]中試驗結(jié)果進(jìn)行比較。由于基于化學(xué)滴定法得到的氯離子濃度實測值代表的是截面氯離子濃度分布的平均結(jié)果,因此為便于與試驗結(jié)果進(jìn)行比較,選擇距離混凝土表面同一深度氯離子濃度模擬值的平均值作為評價指標(biāo)。模型計算參數(shù)分別為:D0=3.931×10-12m2/s,實際氯離子擴(kuò)散系數(shù)根據(jù)式(12)~式(16)確定,m=0.608 3,Δt=0.1 d,δ=2 mm。不同侵蝕時間后截面平均氯離子濃度模擬值和試驗實測值隨深度h、時間的變化曲線見圖9。

    由圖9可知,模型的模擬結(jié)果與試驗值隨侵蝕深度和時間的變化趨勢一致,均隨侵蝕深度的增加而減小,隨侵蝕時間的增加而增加。表明利用MCA模型預(yù)測的氯離子濃度隨侵蝕深度和時間的變化規(guī)律與試驗結(jié)果一致。

    隨后,將MCA模型預(yù)測的截面氯離子濃度模擬值與文獻(xiàn)[24]中實測值進(jìn)行對比,對比結(jié)果見圖10。由圖10可知,除少數(shù)氯離子濃度模擬值與實測值之間的偏差在35%以內(nèi),所有采用MCA模型模擬的截面氯離子濃度模擬值與實測結(jié)果之間的偏差都在±10%以內(nèi)。說明提出的模型可以用于模擬氯離子在混凝土中的擴(kuò)散過程,且精度較高,這也可以通過圖6中模型模擬值與解析解的偏差得出。通過與文獻(xiàn)[24]對比,發(fā)現(xiàn)造成極少數(shù)氯離子濃度模擬值與試驗值偏差的主要原因是模型采用的氯離子擴(kuò)散系數(shù)是文獻(xiàn)[24]通過試驗數(shù)據(jù)擬合后的結(jié)果,而文獻(xiàn)[24]中的試驗數(shù)據(jù)來自不同的試驗環(huán)境,試驗數(shù)據(jù)的采集也來自不同測試人員,使得各試驗數(shù)據(jù)之間本身存在一定的誤差。

    圖10 氯離子濃度實測值與模擬值的對比

    6.3 空間變異性驗證

    與傳統(tǒng)有限元細(xì)觀模型相比,MCA模型除能準(zhǔn)確描述混凝土截面平均氯離子濃度隨時間和深度的變化外,還應(yīng)能描述材料變異特性,如骨料尺寸大小、空間不均勻分布等對氯離子擴(kuò)散行為的影響。為此,在上述分析的基礎(chǔ)上,利用MCA模型給出不同侵蝕時間后截面氯離子濃度分布示意,見圖11。

    由圖11可知,在同一侵蝕時刻,距離混凝土表面同一深度氯離子濃度值各不相同,在砂漿分布區(qū)域(顏色越接近紅色區(qū)域),元胞代表的氯離子濃度最大,在骨料分布區(qū)域(顏色越接近藍(lán)色區(qū)域),元胞代表的氯離子濃度最小,這與氯離子在混凝土中的實際擴(kuò)散行為相符,即氯離子擴(kuò)散行為僅發(fā)生在砂漿區(qū)域,而不考慮骨料的滲透性。說明提出的氯離子擴(kuò)散系數(shù)等效方法是合理的,也證明模型具有描述材料變異特性的能力。

    在應(yīng)用方面,由于MCA模型在求解過程中不需要求解高階微分方程,因此可模擬任意復(fù)雜結(jié)構(gòu)中氯離子的擴(kuò)散行為。利用MCA模型,可以使工程設(shè)計人員直觀了解混凝土截面內(nèi)最容易引起鋼筋銹蝕的薄弱位置,進(jìn)而對結(jié)構(gòu)的耐久性設(shè)計、壽命預(yù)測、維修和加固等提供一定的理論依據(jù)。

    此外,從計算效率來看,當(dāng)δ=2 mm時,MCA模型所需的元胞數(shù)量僅為50個,遠(yuǎn)遠(yuǎn)小于傳統(tǒng)有限元細(xì)觀模型所需的單元數(shù)[2],相比而言,MCA模型大大提高了模型的計算效率。

    7 參數(shù)分析

    在上述對模型的準(zhǔn)確性進(jìn)行驗證的情況下,利用MCA模型進(jìn)行敏感性分析,以驗證骨料分布形式、骨料含量、界面擴(kuò)散特性等細(xì)觀結(jié)構(gòu)參數(shù)對氯離子擴(kuò)散效應(yīng)的影響。純砂漿截面見圖12(a)。為便于討論,將純砂漿截面中氯離子濃度值作為參考,以此來驗證各細(xì)觀參數(shù)對氯離子宏觀擴(kuò)散特性的影響程度。同時為反映各細(xì)觀結(jié)構(gòu)參數(shù)對氯離子擴(kuò)散的影響規(guī)律,引入截面有效氯離子擴(kuò)散系數(shù)平均值Da,其表達(dá)式為

    圖12 不同骨料分布形式的混凝土細(xì)觀模型和砂漿模型示意

    (19)

    式中:D(i,j)為各元胞的氯離子擴(kuò)散系數(shù),m2/s;m、n分別為截面內(nèi)元胞的行數(shù)、列數(shù)。

    7.1 骨料分布形式

    為研究骨料分布形式對混凝土氯離子擴(kuò)散效應(yīng)的影響,考慮骨料在混凝土中的隨機(jī)分布,建立5種不同骨料分布形式的混凝土細(xì)觀模型,見圖12(b)~圖12(f),截面尺寸為100 mm×100 mm。利用MCA模型對5種不同骨料分布形式下混凝土中氯離子擴(kuò)散過程進(jìn)行數(shù)值模擬,計算參數(shù)分別取:Δt=0.1 d,t=480 d,m=0.608 3,D0=3.931×10-12m2/s,實際氯離子擴(kuò)散系數(shù)根據(jù)式(12)~式(16)確定,AAgg=0.4,tITZ=50 μm,δ= 2 mm。

    截面有效氯離子擴(kuò)散系數(shù)平均值Da隨骨料分布形式的變化關(guān)系見圖13。由圖13可知,不同骨料分布形式下Da基本保持不變。說明骨料分布形式基本不會對混凝土中氯離子宏觀分布規(guī)律產(chǎn)生影響。

    圖13 Da隨骨料分布形式的變化關(guān)系

    同時,為便于同其他細(xì)觀參數(shù)比較,選擇距離混凝土表面同一深度氯離子濃度平均值作為評價指標(biāo)。不同骨料分布形式下截面氯離子濃度模擬值與參考值的比較結(jié)果見圖14。由圖14可知,截面各模擬值與參考值之間最大偏差不超過12%。說明在一定誤差范圍內(nèi),可認(rèn)為骨料分布形式基本不會對截面氯離子濃度分布規(guī)律產(chǎn)生影響,這與圖13得出的結(jié)論一致。

    圖14 不同骨料分布形式下截面內(nèi)氯離子濃度模擬值與參考值的比較結(jié)果

    7.2 骨料含量

    為研究骨料含量對氯離子宏觀擴(kuò)散效應(yīng)的影響,對5種不同骨料含量下混凝土中氯離子侵蝕過程進(jìn)行模擬,除AAgg外,各計算參數(shù)、模型尺寸均與上述模型保持一致。截面有效氯離子擴(kuò)散系數(shù)平均值Da隨骨料含量AAgg的變化曲線見圖15。

    圖15 Da隨AAgg的變化曲線

    由圖15可知,隨著骨料含量的增加,Da呈現(xiàn)出線性減小的變化規(guī)律。結(jié)果表明,隨著AAgg的增加,氯離子在混凝土中的宏觀擴(kuò)散特性逐漸減弱,增加混凝土中骨料的含量可以抑制氯離子在混凝土中的擴(kuò)散行為。同時,為便于同其他細(xì)觀結(jié)構(gòu)參數(shù)進(jìn)行比較,選擇距離混凝土表面同一深度氯離子濃度模擬值的平均值作為評價指標(biāo),將上述氯離子濃度模擬值與其平均值進(jìn)行對比,見圖16。

    圖16 不同AAgg下截面內(nèi)氯離子濃度模擬值與參考值的比較

    由圖16可知,截面各位置氯離子濃度模擬值與參考值之間最大誤差為20%,截面各位置氯離子濃度模擬值與參考值之間的偏差隨氯離子濃度值的增加而減小。

    7.3 界面擴(kuò)散特性

    ITZ作為連接骨料和砂漿的重要組成成分,是混凝土中最為薄弱的部位,也是影響混凝土結(jié)構(gòu)力學(xué)性能和耐久性的主要因素。因此,為研究混凝土界面擴(kuò)散特性對氯離子擴(kuò)散效應(yīng)的影響,對8種不同α(DITZ/DMor)下混凝土中氯離子侵蝕過程進(jìn)行模擬。除α外,各計算參數(shù)、模型尺寸均與上述模型保持一致。截面有效氯離子擴(kuò)散系數(shù)平均值Da隨α的變化曲線見圖17。

    圖17 Da隨 α的變化曲線

    由圖17可知,Da隨α的增加呈現(xiàn)出線性增加的變化規(guī)律。結(jié)果表明,隨著α的增加,氯離子在混凝土中的宏觀擴(kuò)散特性逐漸增強(qiáng)。

    同時,為便于同其他參數(shù)進(jìn)行對比,選擇距離混凝土表面同一深度氯離子濃度平均值作為評價指標(biāo),將上述不同α下截面氯離子濃度模擬值與參考值進(jìn)行比較,結(jié)果見圖18。

    圖18 不同α下截面內(nèi)氯離子濃度模擬值與參考值的對比結(jié)果

    由圖18可知,不同α下截面氯離子濃度模擬值之間的差值隨氯離子濃度的減小而增加。當(dāng)α=15時,兩者之間偏差最大,為28%??梢姦恋娜≈祵β入x子在混凝土中的宏觀擴(kuò)散特性存在很大的影響,在氯離子擴(kuò)散效應(yīng)分析時,合理選擇α值將直接關(guān)乎結(jié)果的正確與否。

    綜上所述,通過對比不同參數(shù)下截面中氯離子濃度模擬值與參考值之間的偏差,發(fā)現(xiàn)界面擴(kuò)散特性對細(xì)觀混凝土中氯離子的宏觀擴(kuò)散特性影響最顯著,兩者之間偏差最大,為28%,體積分?jǐn)?shù)對細(xì)觀混凝土中氯離子的宏觀擴(kuò)散特性影響較為顯著,兩者之間偏差為20%,而骨料分布形式則基本不會對混凝土氯離子擴(kuò)散效應(yīng)產(chǎn)生影響,兩者之間偏差最小,僅為14%,這與既有研究結(jié)果規(guī)律一致[13,19-20]。

    8 結(jié)論

    本文基于細(xì)觀尺度模擬方法,將混凝土看成由砂漿、界面過渡區(qū)和骨料組成的三相復(fù)合材料,建立相應(yīng)的隨機(jī)細(xì)觀模型,提出考慮氯離子時變擴(kuò)散效應(yīng)的MCA模型。利用試驗實測數(shù)據(jù)和理論解對模型進(jìn)行驗證,并在此基礎(chǔ)上利用模型進(jìn)行敏感性分析。最后得出以下結(jié)論:

    (1)根據(jù)混凝土組成特點(diǎn),將混凝土看成由砂漿、ITZ、骨料組成的三相復(fù)合材料,基于MC方法建立考慮骨料等不均勻分布的混凝土隨機(jī)細(xì)觀模型。并在此基礎(chǔ)上,結(jié)合元胞自動機(jī)理論,提出一種均勻化等效方法,建立用于描述氯離子在混凝土中擴(kuò)散行為的MCA模型。

    (2)利用模型對影響模型計算效率和應(yīng)用限制的元胞尺寸大小進(jìn)行討論,為保證計算精度和準(zhǔn)確反映氯離子在各相組成材料中的擴(kuò)散特性,推薦采用δ= 0.707dmin。 與傳統(tǒng)有限元細(xì)觀模型相比,提出的模型有效解決了有限元模型由于網(wǎng)格尺寸限制而無法考慮實際界面過渡區(qū)尺寸和計算效率低的難題。

    (3)通過與實驗結(jié)果和理論解的比較,驗證MCA模型的準(zhǔn)確性和描述材料變異特性的能力。模型能夠反映混凝土細(xì)觀組成材料(如砂漿、骨料、ITZ)中氯離子的實際擴(kuò)散特性。提出的方法為細(xì)觀結(jié)構(gòu)和宏觀模型之間建立了橋梁,這對進(jìn)一步整體截面或結(jié)構(gòu)分析具有重要意義。

    (4)參數(shù)化分析結(jié)果表明,界面擴(kuò)散特性對細(xì)觀混凝土中氯離子的宏觀擴(kuò)散特性影響最顯著,體積分?jǐn)?shù)對細(xì)觀混凝土中氯離子的宏觀擴(kuò)散特性影響較為顯著,而骨料分布形式基本不會對氯離子擴(kuò)散效應(yīng)產(chǎn)生影響。總體規(guī)律是氯離子擴(kuò)散效應(yīng)隨界面擴(kuò)散特性的增加而增加,隨骨料體積分?jǐn)?shù)的增加而減小。

    猜你喜歡
    元胞細(xì)觀擴(kuò)散系數(shù)
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    基于元胞自動機(jī)下的交通事故路段仿真
    智富時代(2018年5期)2018-07-18 17:52:04
    基于元胞數(shù)據(jù)的多維數(shù)據(jù)傳遞機(jī)制
    北京測繪(2016年2期)2016-01-24 02:28:28
    基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴(kuò)散系數(shù)的研究
    上海金屬(2015年5期)2015-11-29 01:13:59
    FCC Ni-Cu 及Ni-Mn 合金互擴(kuò)散系數(shù)測定
    上海金屬(2015年6期)2015-11-29 01:09:09
    非時齊擴(kuò)散模型中擴(kuò)散系數(shù)的局部估計
    基于四叉樹網(wǎng)格加密技術(shù)的混凝土細(xì)觀模型
    基于AIS的航道移動瓶頸元胞自動機(jī)模型
    中國航海(2014年1期)2014-05-09 07:54:25
    PBX炸藥的抗壓強(qiáng)度及抗拉強(qiáng)度細(xì)觀尺度的數(shù)值計算
    開裂混凝土中水分傳輸過程的細(xì)觀模型
    av线在线观看网站| 老师上课跳d突然被开到最大视频| 欧美成人a在线观看| 欧美精品国产亚洲| 黄片wwwwww| 最近中文字幕2019免费版| 久久久久久久亚洲中文字幕| 美女高潮的动态| 三级毛片av免费| 性色avwww在线观看| 免费播放大片免费观看视频在线观看| 九草在线视频观看| 色网站视频免费| 97精品久久久久久久久久精品| 我的女老师完整版在线观看| 免费av观看视频| 久久97久久精品| 国产高清不卡午夜福利| 免费黄网站久久成人精品| 国产乱来视频区| 建设人人有责人人尽责人人享有的 | 最近视频中文字幕2019在线8| 91aial.com中文字幕在线观看| 亚洲av国产av综合av卡| 不卡视频在线观看欧美| 人人妻人人看人人澡| 国产视频内射| 国产综合精华液| 少妇丰满av| 久久久久久久久久久免费av| 久久精品国产亚洲av天美| 亚洲精品久久午夜乱码| 男女下面进入的视频免费午夜| 麻豆成人av视频| 国产乱来视频区| 日韩成人av中文字幕在线观看| av国产免费在线观看| 免费大片黄手机在线观看| 青春草视频在线免费观看| 只有这里有精品99| 亚洲精品久久久久久婷婷小说| 欧美最新免费一区二区三区| 亚洲av成人精品一区久久| www.av在线官网国产| 亚洲国产高清在线一区二区三| 日韩伦理黄色片| 美女xxoo啪啪120秒动态图| 18禁在线无遮挡免费观看视频| 国产视频内射| 日韩在线高清观看一区二区三区| 国产精品久久久久久久电影| 亚洲一级一片aⅴ在线观看| 免费看美女性在线毛片视频| 一级黄片播放器| 久久久a久久爽久久v久久| 岛国毛片在线播放| 亚洲欧洲国产日韩| 久久韩国三级中文字幕| 99久国产av精品| 亚洲人成网站在线观看播放| 日韩av不卡免费在线播放| 亚洲自偷自拍三级| 国产片特级美女逼逼视频| 极品少妇高潮喷水抽搐| 日本-黄色视频高清免费观看| 日韩不卡一区二区三区视频在线| 亚洲aⅴ乱码一区二区在线播放| 色视频www国产| 亚洲在线自拍视频| 你懂的网址亚洲精品在线观看| 99久久精品热视频| 欧美成人一区二区免费高清观看| 亚洲国产成人一精品久久久| 看免费成人av毛片| 搡女人真爽免费视频火全软件| 搡老妇女老女人老熟妇| 国产精品蜜桃在线观看| 久久这里有精品视频免费| 国产午夜福利久久久久久| 久久99热这里只有精品18| 深爱激情五月婷婷| 少妇人妻精品综合一区二区| 国产大屁股一区二区在线视频| 网址你懂的国产日韩在线| 免费观看a级毛片全部| 国产精品爽爽va在线观看网站| 久久久色成人| 国产成人精品一,二区| 97超视频在线观看视频| 三级毛片av免费| 亚洲国产精品国产精品| 久久久久久伊人网av| 夜夜爽夜夜爽视频| 久久久久久久大尺度免费视频| 成人特级av手机在线观看| 丰满人妻一区二区三区视频av| 岛国毛片在线播放| 国语对白做爰xxxⅹ性视频网站| 日韩在线高清观看一区二区三区| 99热6这里只有精品| 男人爽女人下面视频在线观看| 午夜福利在线观看吧| 一本久久精品| 欧美日韩国产mv在线观看视频 | 亚洲成人一二三区av| 在线免费观看的www视频| 大片免费播放器 马上看| 女的被弄到高潮叫床怎么办| .国产精品久久| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 亚洲国产色片| 国产伦理片在线播放av一区| eeuss影院久久| 一级片'在线观看视频| 成人毛片a级毛片在线播放| 午夜福利视频1000在线观看| 午夜亚洲福利在线播放| 一边亲一边摸免费视频| av一本久久久久| 国产成人a∨麻豆精品| 久久韩国三级中文字幕| 听说在线观看完整版免费高清| 中文精品一卡2卡3卡4更新| 免费不卡的大黄色大毛片视频在线观看 | 亚洲不卡免费看| 亚洲av电影不卡..在线观看| 国产精品嫩草影院av在线观看| 校园人妻丝袜中文字幕| 高清av免费在线| av播播在线观看一区| 免费看a级黄色片| 肉色欧美久久久久久久蜜桃 | 男女那种视频在线观看| 国产v大片淫在线免费观看| 国产人妻一区二区三区在| av卡一久久| 中国国产av一级| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久国内精品自在自线图片| 欧美成人一区二区免费高清观看| 国产精品久久视频播放| 亚洲av电影在线观看一区二区三区 | 免费黄色在线免费观看| 97超碰精品成人国产| 亚洲欧美精品专区久久| 国产精品久久视频播放| 波多野结衣巨乳人妻| 欧美区成人在线视频| 欧美潮喷喷水| 七月丁香在线播放| 中文天堂在线官网| 亚洲av二区三区四区| 免费黄频网站在线观看国产| 亚洲国产精品成人久久小说| 国产女主播在线喷水免费视频网站 | 久久久精品免费免费高清| 国产精品一区www在线观看| 麻豆精品久久久久久蜜桃| 水蜜桃什么品种好| 丝袜美腿在线中文| 久久韩国三级中文字幕| 国产高清国产精品国产三级 | 日本免费a在线| 少妇被粗大猛烈的视频| 精品99又大又爽又粗少妇毛片| 人妻系列 视频| 国产人妻一区二区三区在| 亚洲成人一二三区av| 午夜福利视频1000在线观看| 天天躁日日操中文字幕| 美女xxoo啪啪120秒动态图| 成人午夜高清在线视频| 午夜福利视频1000在线观看| 色视频www国产| 一级毛片我不卡| 国产男人的电影天堂91| 国产综合懂色| 乱系列少妇在线播放| 97精品久久久久久久久久精品| 久久97久久精品| 精品久久久久久久久av| 91aial.com中文字幕在线观看| 亚洲真实伦在线观看| 国产精品无大码| 草草在线视频免费看| 国产精品久久久久久精品电影| 最近最新中文字幕免费大全7| 春色校园在线视频观看| 18+在线观看网站| 亚洲成人精品中文字幕电影| 亚洲精品久久午夜乱码| 国产大屁股一区二区在线视频| 国产精品一区二区三区四区免费观看| 中文精品一卡2卡3卡4更新| 亚洲精品第二区| 欧美高清成人免费视频www| 午夜福利在线在线| 国产极品天堂在线| 国产又色又爽无遮挡免| 久久人人爽人人片av| 精品久久久精品久久久| 亚洲国产最新在线播放| 色视频www国产| 秋霞伦理黄片| 亚洲丝袜综合中文字幕| 在线播放无遮挡| 亚洲无线观看免费| 日韩欧美国产在线观看| 精品一区二区三区视频在线| 免费黄色在线免费观看| 国产久久久一区二区三区| 成年av动漫网址| 久久久a久久爽久久v久久| 久久久亚洲精品成人影院| 超碰av人人做人人爽久久| 欧美成人精品欧美一级黄| 美女脱内裤让男人舔精品视频| 国产av不卡久久| 亚洲成人久久爱视频| 久久久久久久久久成人| 97超视频在线观看视频| 成年免费大片在线观看| 身体一侧抽搐| 久久久亚洲精品成人影院| 少妇高潮的动态图| 亚洲自偷自拍三级| 国产色婷婷99| 在线观看美女被高潮喷水网站| 国产成年人精品一区二区| 欧美一区二区亚洲| 一级片'在线观看视频| 3wmmmm亚洲av在线观看| 插逼视频在线观看| 人人妻人人看人人澡| 免费看日本二区| 春色校园在线视频观看| 日本爱情动作片www.在线观看| 精品久久国产蜜桃| 国产精品一及| 精品久久久久久久久亚洲| 亚洲国产av新网站| 淫秽高清视频在线观看| 大香蕉97超碰在线| 国产在视频线在精品| 成人亚洲精品一区在线观看 | 最新中文字幕久久久久| 国产精品久久视频播放| 亚洲在久久综合| xxx大片免费视频| 久久精品久久久久久久性| 成人亚洲欧美一区二区av| 日韩av免费高清视频| 亚洲国产最新在线播放| 人妻制服诱惑在线中文字幕| 欧美另类一区| 国产 一区 欧美 日韩| 免费看a级黄色片| 禁无遮挡网站| 国产成人精品福利久久| 一级毛片电影观看| 欧美成人精品欧美一级黄| 我的老师免费观看完整版| 视频中文字幕在线观看| 国产成人精品一,二区| 亚洲精品成人久久久久久| 水蜜桃什么品种好| 成年女人看的毛片在线观看| 日本免费a在线| 精品一区二区三卡| 久久这里只有精品中国| 久热久热在线精品观看| 婷婷色麻豆天堂久久| 亚洲精品国产成人久久av| 麻豆精品久久久久久蜜桃| 亚洲电影在线观看av| 一个人免费在线观看电影| 极品少妇高潮喷水抽搐| 3wmmmm亚洲av在线观看| 欧美3d第一页| 男女边摸边吃奶| 能在线免费看毛片的网站| 午夜视频国产福利| 日日撸夜夜添| 久久精品国产亚洲av天美| 日韩人妻高清精品专区| 亚洲国产色片| 欧美xxⅹ黑人| 国产精品久久视频播放| 三级经典国产精品| 久久综合国产亚洲精品| 日本一本二区三区精品| 精品久久久久久电影网| 99热这里只有是精品在线观看| 一级毛片黄色毛片免费观看视频| 97超视频在线观看视频| 久久久久网色| 久久精品国产鲁丝片午夜精品| 亚洲精品456在线播放app| 色综合亚洲欧美另类图片| 少妇高潮的动态图| 国产av码专区亚洲av| 国内精品美女久久久久久| 日韩精品青青久久久久久| 亚洲av国产av综合av卡| 亚洲精品日韩av片在线观看| 免费av毛片视频| 草草在线视频免费看| 精品一区在线观看国产| 午夜福利视频1000在线观看| 男人狂女人下面高潮的视频| 九草在线视频观看| 国产精品综合久久久久久久免费| 美女大奶头视频| 看非洲黑人一级黄片| 国产国拍精品亚洲av在线观看| 成人av在线播放网站| 免费少妇av软件| 中文字幕av成人在线电影| 青青草视频在线视频观看| 麻豆乱淫一区二区| 禁无遮挡网站| 尤物成人国产欧美一区二区三区| 久久久精品免费免费高清| 三级男女做爰猛烈吃奶摸视频| 亚洲欧洲国产日韩| 久久午夜福利片| videos熟女内射| 黄色日韩在线| a级毛片免费高清观看在线播放| 国产片特级美女逼逼视频| 日韩欧美国产在线观看| 免费高清在线观看视频在线观看| 国产精品久久久久久av不卡| 在线观看人妻少妇| 2022亚洲国产成人精品| 淫秽高清视频在线观看| 麻豆av噜噜一区二区三区| 老司机影院毛片| 成人特级av手机在线观看| 日韩不卡一区二区三区视频在线| 国国产精品蜜臀av免费| 国产片特级美女逼逼视频| 中文字幕人妻熟人妻熟丝袜美| 精品国产一区二区三区久久久樱花 | 少妇人妻一区二区三区视频| 国产高清国产精品国产三级 | 免费av毛片视频| 亚洲国产色片| 小蜜桃在线观看免费完整版高清| 肉色欧美久久久久久久蜜桃 | 五月伊人婷婷丁香| 99热6这里只有精品| 国产精品久久久久久久电影| 噜噜噜噜噜久久久久久91| 免费黄频网站在线观看国产| 3wmmmm亚洲av在线观看| 成年人午夜在线观看视频 | 在线免费观看的www视频| 久久久久性生活片| 大香蕉久久网| 久久这里只有精品中国| 干丝袜人妻中文字幕| 国产av码专区亚洲av| 久久97久久精品| 亚洲av在线观看美女高潮| 国产午夜精品久久久久久一区二区三区| 激情 狠狠 欧美| 最近最新中文字幕大全电影3| 日韩av不卡免费在线播放| av在线天堂中文字幕| 99久久人妻综合| 看非洲黑人一级黄片| 91久久精品电影网| 黑人高潮一二区| 18禁在线无遮挡免费观看视频| 乱码一卡2卡4卡精品| 国产成人精品久久久久久| 久久久精品欧美日韩精品| 天美传媒精品一区二区| 干丝袜人妻中文字幕| 国产精品一区二区在线观看99 | 日韩制服骚丝袜av| 国产片特级美女逼逼视频| av在线亚洲专区| 久久国产乱子免费精品| 好男人视频免费观看在线| 成人综合一区亚洲| 国产午夜福利久久久久久| 日韩三级伦理在线观看| 精品久久久精品久久久| 亚洲国产日韩欧美精品在线观看| 亚洲精品日韩av片在线观看| 国产在视频线精品| 国产淫语在线视频| 搡女人真爽免费视频火全软件| 亚洲av免费高清在线观看| 日韩欧美一区视频在线观看 | 亚洲人成网站在线播| 亚洲色图av天堂| 免费观看a级毛片全部| 欧美性感艳星| 欧美成人午夜免费资源| 国产大屁股一区二区在线视频| 国产美女午夜福利| 免费人成在线观看视频色| ponron亚洲| 综合色av麻豆| 国产亚洲一区二区精品| kizo精华| 午夜福利在线观看免费完整高清在| 国产黄片视频在线免费观看| 国产精品蜜桃在线观看| 欧美区成人在线视频| 欧美bdsm另类| 少妇猛男粗大的猛烈进出视频 | 美女cb高潮喷水在线观看| 国产色爽女视频免费观看| ponron亚洲| 亚洲高清免费不卡视频| 国产精品麻豆人妻色哟哟久久 | 少妇裸体淫交视频免费看高清| 18禁动态无遮挡网站| 精品国产露脸久久av麻豆 | 一区二区三区四区激情视频| 别揉我奶头 嗯啊视频| 午夜福利高清视频| 国产亚洲午夜精品一区二区久久 | 国产精品国产三级国产av玫瑰| 亚洲美女视频黄频| 精品欧美国产一区二区三| 久久久久久久亚洲中文字幕| 日日啪夜夜撸| 久久久久九九精品影院| 欧美区成人在线视频| 少妇裸体淫交视频免费看高清| 中文字幕免费在线视频6| 亚洲欧美日韩东京热| 亚洲av男天堂| 亚洲最大成人手机在线| 亚洲精品一区蜜桃| 成人毛片60女人毛片免费| 国产成人精品久久久久久| 免费看美女性在线毛片视频| 夜夜爽夜夜爽视频| av播播在线观看一区| 女人十人毛片免费观看3o分钟| 熟妇人妻不卡中文字幕| 国产又色又爽无遮挡免| 日本爱情动作片www.在线观看| 男人爽女人下面视频在线观看| 美女脱内裤让男人舔精品视频| 免费看美女性在线毛片视频| 天堂网av新在线| 国产一级毛片七仙女欲春2| 午夜亚洲福利在线播放| 青春草视频在线免费观看| av网站免费在线观看视频 | 国产爱豆传媒在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲无线观看免费| 久久这里只有精品中国| 韩国av在线不卡| 亚洲国产高清在线一区二区三| 精品一区在线观看国产| 干丝袜人妻中文字幕| 日韩电影二区| 欧美变态另类bdsm刘玥| 80岁老熟妇乱子伦牲交| 日韩一本色道免费dvd| 一本久久精品| 成人美女网站在线观看视频| 97人妻精品一区二区三区麻豆| 丝袜美腿在线中文| 乱码一卡2卡4卡精品| 神马国产精品三级电影在线观看| 一个人免费在线观看电影| 亚洲在线自拍视频| 最新中文字幕久久久久| 插阴视频在线观看视频| 亚洲高清免费不卡视频| 男人爽女人下面视频在线观看| 少妇的逼水好多| 男女边摸边吃奶| 国产视频首页在线观看| 人妻少妇偷人精品九色| 午夜久久久久精精品| 真实男女啪啪啪动态图| 亚洲av国产av综合av卡| 最近中文字幕高清免费大全6| av线在线观看网站| 2022亚洲国产成人精品| 亚洲精品影视一区二区三区av| 国产探花在线观看一区二区| 亚洲成人精品中文字幕电影| 欧美性猛交╳xxx乱大交人| 亚洲久久久久久中文字幕| 永久网站在线| 亚洲综合精品二区| 成人高潮视频无遮挡免费网站| 高清av免费在线| 亚洲在线自拍视频| 亚洲欧美日韩无卡精品| 国产片特级美女逼逼视频| 99re6热这里在线精品视频| av播播在线观看一区| 国产av国产精品国产| 可以在线观看毛片的网站| 亚洲熟女精品中文字幕| 男插女下体视频免费在线播放| 中文字幕av在线有码专区| www.色视频.com| 国产单亲对白刺激| 日本免费在线观看一区| 亚洲精品久久午夜乱码| 国产色爽女视频免费观看| 自拍偷自拍亚洲精品老妇| 亚洲av男天堂| 五月天丁香电影| 伦精品一区二区三区| 国产成人精品婷婷| 久久久久精品性色| 黄色一级大片看看| 国产综合精华液| 国产欧美日韩精品一区二区| 亚洲精品乱码久久久久久按摩| 18禁在线播放成人免费| 亚洲精品一区蜜桃| 国产成人freesex在线| 色综合色国产| 哪个播放器可以免费观看大片| 亚洲av福利一区| 国产乱人偷精品视频| 只有这里有精品99| 国产人妻一区二区三区在| 欧美+日韩+精品| 亚洲最大成人中文| 国产av不卡久久| 午夜精品一区二区三区免费看| 国产国拍精品亚洲av在线观看| 亚洲精品一二三| 天天一区二区日本电影三级| 亚洲自偷自拍三级| 精品人妻熟女av久视频| 国产午夜福利久久久久久| 国产精品美女特级片免费视频播放器| 国产精品.久久久| 国产黄片视频在线免费观看| 波多野结衣巨乳人妻| 亚洲成人av在线免费| 日韩欧美三级三区| 国产成人精品一,二区| 日韩大片免费观看网站| 中文字幕制服av| 久久久久久九九精品二区国产| 国产亚洲av嫩草精品影院| 国产一区二区三区综合在线观看 | 九草在线视频观看| 免费不卡的大黄色大毛片视频在线观看 | 免费观看的影片在线观看| 日韩成人伦理影院| 精品久久久久久久久av| 欧美极品一区二区三区四区| 国产在视频线在精品| 人人妻人人澡人人爽人人夜夜 | 国产爱豆传媒在线观看| 一二三四中文在线观看免费高清| 亚洲久久久久久中文字幕| 少妇丰满av| 性色avwww在线观看| 亚洲欧美一区二区三区黑人 | 又黄又爽又刺激的免费视频.| 2018国产大陆天天弄谢| 大又大粗又爽又黄少妇毛片口| 91久久精品国产一区二区成人| 日韩av在线免费看完整版不卡| 一级黄片播放器| 国产真实伦视频高清在线观看| 午夜爱爱视频在线播放| 国产精品久久久久久精品电影| 99久国产av精品| 精品久久国产蜜桃| 免费看不卡的av| 在线观看美女被高潮喷水网站| 人妻夜夜爽99麻豆av| 97精品久久久久久久久久精品| 亚洲精品亚洲一区二区| 一级av片app| 久久久久精品久久久久真实原创| 免费av观看视频| 日本爱情动作片www.在线观看| 日韩国内少妇激情av| 国产亚洲av嫩草精品影院| 少妇被粗大猛烈的视频| 成人午夜精彩视频在线观看| 天美传媒精品一区二区| 日韩一本色道免费dvd| 老女人水多毛片| 中文精品一卡2卡3卡4更新| 麻豆成人av视频| 欧美一区二区亚洲| 插逼视频在线观看| 国产精品蜜桃在线观看| av一本久久久久| av在线亚洲专区| 禁无遮挡网站| 精品国产三级普通话版| 亚洲在线观看片| 日日摸夜夜添夜夜爱| 日本熟妇午夜| 最近中文字幕高清免费大全6| 国产精品综合久久久久久久免费| 国产精品嫩草影院av在线观看| 亚洲欧洲国产日韩| 国产又色又爽无遮挡免|