999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

溫度和濕度環(huán)境下未鍍鋅高強(qiáng)鋼絲的腐蝕速率譜

2018-04-11 09:06:43陳小雨唐茂林
關(guān)鍵詞:測(cè)量環(huán)境

陳小雨,唐茂林

(西南交通大學(xué)土木工程學(xué)院,四川 成都 610031)

溫度和濕度是影響所有金屬材料腐蝕的重要環(huán)境因素[1].在外界環(huán)境的變化下,纜索承載橋梁纜索內(nèi)的水分蒸發(fā)使得溫度和濕度常常高于外界環(huán)境.然而,纜索鋼絲間縫隙小且承受巨大拉力,很難實(shí)現(xiàn)對(duì)其腐蝕實(shí)時(shí)監(jiān)控[2].傳統(tǒng)方法是將主纜楔開,將傳感器埋入主纜,這對(duì)主纜造成傷害[3-4],因此,許多研究采用試驗(yàn)室模擬的方法研究環(huán)境對(duì)鋼絲腐蝕的影響.

文獻(xiàn)[5]中分別用濕毛巾,密封箱等模擬了3種濕度條件下鋼絲的腐蝕,得到3種濕度下鋼絲腐蝕的平均速率.文獻(xiàn)[6-7]中采用恒溫箱將鋼絲浸入酸雨溶液,研究不同溫度和pH值對(duì)鋼絲腐蝕的影響.文獻(xiàn)[8]中采用多功能環(huán)境試驗(yàn)箱,模擬懸索橋服役環(huán)境,測(cè)量得到當(dāng)外界環(huán)境變化時(shí),纜內(nèi)的溫度濕度以及腐蝕速率的變化.文獻(xiàn)[9]中試驗(yàn)研究了溫度、含鹽量、pH值等對(duì)鋼絲的腐蝕影響,得到了三者的回歸關(guān)系.文獻(xiàn)[10]中綜合分析了各腐蝕因素對(duì)鋼絲腐蝕速率的影響,并對(duì)這些因素進(jìn)行了主次關(guān)系的排序.

上述研究已經(jīng)對(duì)鋼絲腐蝕環(huán)境進(jìn)行了討論,但由于測(cè)量方法的限制,并沒(méi)有得到腐蝕速率與溫度、濕度之間的譜,缺乏研究三者之間的回歸方程.

由于普通傳感器測(cè)量原理的限制,需要被測(cè)鋼絲與測(cè)量電極間有溶液作為介質(zhì),才能進(jìn)行電化學(xué)測(cè)量,因此,不同濕度下腐蝕速率測(cè)量一直是難以解決的問(wèn)題.本研究采用極化電阻傳感器,將電極鑲嵌在導(dǎo)電薄膜上,滿足了在不同濕度下的測(cè)量要求.使用恒溫恒濕試驗(yàn)箱控制試驗(yàn)溫度和濕度,采用正交試驗(yàn)方法,試驗(yàn)研究在溫度、濕度共同環(huán)境作用下的未鍍鋅高強(qiáng)鋼絲腐蝕速率譜.

1 腐蝕速率譜測(cè)量與計(jì)算

1.1 腐蝕速率測(cè)量原理

金屬表面的腐蝕實(shí)質(zhì)是金屬離子與空氣或水中的氧分子發(fā)生電化學(xué)反應(yīng),腐蝕過(guò)程中金屬失去電子成為陽(yáng)極,氧氣等得到電子成為陰極,電極之間的離子轉(zhuǎn)移形成了原電池,發(fā)生的化學(xué)反應(yīng)為

其中:z為金屬的電子價(jià).

依據(jù)法拉第定律,電荷的轉(zhuǎn)移與反應(yīng)物質(zhì)的變化量存在等量關(guān)系,可以通過(guò)測(cè)量電極間的極化電流而得到腐蝕速率,具體的電化學(xué)測(cè)量原理見文獻(xiàn)[11].

通常測(cè)量腐蝕速率可以采用大氣暴露和試驗(yàn)室模擬兩種方法.大氣暴露是測(cè)量長(zhǎng)時(shí)間內(nèi)大氣條件下金屬腐蝕的平均速率,其中包含了影響腐蝕的眾多因素.本文設(shè)計(jì)的試驗(yàn),是測(cè)量短期內(nèi)不變的試驗(yàn)環(huán)境下樣本鋼絲的瞬時(shí)腐蝕速率,目的在于考察溫度、濕度兩個(gè)因素影響下腐蝕的速率譜,在數(shù)量上與大氣暴露測(cè)量的數(shù)據(jù)存在差異.該試驗(yàn)的目的是將復(fù)雜的腐蝕環(huán)境進(jìn)行分解,獲得每個(gè)因素對(duì)腐蝕速率影響的擬合計(jì)算公式,最后得到綜合、通用的計(jì)算腐蝕速率譜公式,以便對(duì)橋梁的實(shí)際服役環(huán)境進(jìn)行模擬.

1.2 試驗(yàn)溫度、濕度取值

由于橋位所處環(huán)境差異,溫度和濕度的變化范圍與橋位所在環(huán)境有一定關(guān)系.依據(jù)已有的實(shí)際橋溫度和濕度監(jiān)控資料[12],夏季主纜內(nèi)最高溫度可達(dá)到50 ℃以上,取試驗(yàn)的最高溫度為50 ℃;考慮0 ℃為冰點(diǎn),鋼絲結(jié)冰凍融與腐蝕的交疊效應(yīng)復(fù)雜,最低溫度取10 ℃.

主纜除了被浸泡的局部鋼絲外,最大濕度就是99%,則濕度上限為99%.

由主纜除濕系統(tǒng)的濕度設(shè)置可知,在濕度為60%以下的環(huán)境中可以忽略鋼絲的腐蝕[13],為了得到更寬范圍的濕度影響,為安裝除濕系統(tǒng)后的腐蝕速率提供參考,將最低濕度取為40%,恒溫恒濕箱采用的供濕液體為中性純凈水,pH值為6.5~8.5.

溫度取值間隔為10 ℃,一共5個(gè)溫度取值,為了滿足正交試驗(yàn)的整齊可比性,正交表格中濕度也取5個(gè)值,如表1所示.

表1 正交試驗(yàn)各因素取值設(shè)計(jì)表Tab.1 Values of each factor of orthogonal experiment design table

由于濕度大于60%后,濕度對(duì)腐蝕影響將增加,因此在正交試驗(yàn)測(cè)量的基礎(chǔ)上,增加濕度大于60%的濕度取值,為擬合腐蝕速率譜提供更多樣本點(diǎn),在每個(gè)溫度條件下,均測(cè)量8個(gè)不同濕度的腐蝕速率,插入后的濕度取值如表2所示.

按上述表格取值,共有40個(gè)溫度和濕度環(huán)境組合,每個(gè)組合下測(cè)量18根鋼絲的腐蝕速率,每根鋼絲測(cè)量3次,每個(gè)環(huán)境下共有54個(gè)樣本值.所有試驗(yàn)共計(jì)樣本數(shù)據(jù)為2 160個(gè).

表2 濕度因素取值設(shè)計(jì)表Tab.2 Values of humidity factor design table

2 試驗(yàn)環(huán)境及試驗(yàn)測(cè)量方法

本試驗(yàn)采用Micro-sized linear polarization resistance (μLPR)極化電阻腐蝕速率傳感器,測(cè)量未鍍鋅高強(qiáng)鋼絲的腐蝕速率,將傳感器貼于鋼絲表面即可進(jìn)行實(shí)時(shí)的腐蝕速率測(cè)量[14],傳感器及恒溫恒濕環(huán)境如圖1所示.

試驗(yàn)樣本采用強(qiáng)度為1 770 MPa的未鍍鋅高強(qiáng)鋼絲,鋼絲直徑5.1 mm,樣本鋼絲盤條金屬成分及質(zhì)量比分別為:

wC=0.830%;

wMn=0.750%;

wSi=0.007%;

wP=0.008%;

wCr=0.140%;

wCu=0.030%;

其余為Fe.

盤條金相顯微組織為索氏體(92%)、片狀珠光體及晶界上少量條狀滲碳體.

試驗(yàn)所用樣本鋼絲均來(lái)自同一盤條,將樣本切割成長(zhǎng)為30 cm,6個(gè)樣本鋼絲為一組,每組測(cè)量10 min,同一批樣本反復(fù)測(cè)量3次.樣本用清水洗凈,丙酮除油污,再用清水沖洗后風(fēng)干.將傳感器貼于樣本表面后,放入恒溫恒濕箱進(jìn)行測(cè)量.

圖1 μLPR極化電阻腐蝕傳感器及恒溫恒濕環(huán)境Fig.1 Miniature linear polarization resistance corrosion sensors (μLPR); temperature and humidity chamber

3 試驗(yàn)數(shù)據(jù)處理

3.1 樣本數(shù)據(jù)取值方法

測(cè)量得到樣本數(shù)據(jù)量巨大,為了得到合理的試驗(yàn)數(shù)據(jù),按圖2的流程進(jìn)行數(shù)據(jù)處理,圖中:xi為任意樣本試驗(yàn)值.

統(tǒng)計(jì)測(cè)量樣本后,將每一個(gè)環(huán)境組合下的數(shù)據(jù)進(jìn)行分組統(tǒng)計(jì).計(jì)算組間最大最小值的數(shù)據(jù)全距,依據(jù)全距將數(shù)據(jù)分為5組,略去數(shù)據(jù)個(gè)數(shù)少于5的組,求出其余數(shù)據(jù)的平均值,該值作為該溫度和濕度條件下的試驗(yàn)測(cè)量均值.

3.2 試驗(yàn)數(shù)據(jù)結(jié)果

按上述試驗(yàn)測(cè)量結(jié)果,試驗(yàn)數(shù)據(jù)取試驗(yàn)的均值,如表3所示.

由表3可知,在濕度低于60%時(shí)腐蝕速率不為0,腐蝕速率隨濕度增高而呈現(xiàn)增大的趨勢(shì),并且較明顯,但隨溫度的增高并沒(méi)有呈現(xiàn)出一致的增大趨勢(shì).在40、50 ℃的高溫環(huán)境中,濕度超過(guò)70%,腐蝕速率明顯加大.

圖2 試驗(yàn)數(shù)據(jù)處理流程Fig.2 Experimental data processing flow chart

表3 不同溫度濕度組合條件下鋼絲腐蝕速率值Tab.3 Corrosion rate values of steel wires for various temperature and humidity conditions ×10-4g/(dm2·h)

4 試驗(yàn)結(jié)果分析

4.1 溫度、濕度耦合關(guān)系判斷

考察溫度、濕度兩個(gè)因素在對(duì)未鍍鋅高強(qiáng)鋼絲腐蝕的作用中是否存在耦合關(guān)系.采用因子設(shè)計(jì)法對(duì)數(shù)據(jù)進(jìn)行處理[15].

A為溫度因素:

極小值=10 ℃;

極大值=50 ℃.

B為濕度因素:

極小值=40%~50%;

極大值=90%.

溫度、濕度與腐蝕速率的因子如表4所示,表中,k為影響均值.

表4 溫度、濕度與腐蝕速率的2k因子Tab.4 Temperature,humidity and corrosion rate of 2k factor table

依據(jù)2k因子的設(shè)計(jì)方法[15],計(jì)算如下:

溫度對(duì)腐蝕速率的影響為

濕度對(duì)腐蝕速率的影響為

濕度和計(jì)算溫度耦合作用的影響為

由計(jì)算結(jié)果可知,B>A>AB,即濕度對(duì)腐蝕速率的影響要大于溫度對(duì)腐蝕速率的影響,且兩者交互作用對(duì)速率的影響小于兩者獨(dú)立的影響,但由于交互效應(yīng)影響與A因素單獨(dú)影響的結(jié)果相差不大,仍需要考慮其交互作用.

4.2 試驗(yàn)數(shù)據(jù)極差

極差分析通過(guò)對(duì)比相同因素不同取值間的差異,得到影響腐蝕速率最大的取值組合.根據(jù)不同因素相同取值間的差異,可以確定兩個(gè)因素的主次關(guān)系.表5為正交試驗(yàn)結(jié)果極差分析,表中:k1~k8為A、B因素在各個(gè)取值條件下的測(cè)量數(shù)據(jù)均值;R為最大最小值之差的均值[15].

表5 正交試驗(yàn)結(jié)果極差分析Tab.5 Range analysis table of orthogonal test results

由表5可知,B在不同取值水平對(duì)應(yīng)的腐蝕速率之間差異最大,RB=2.570;A在不同取值水平的腐蝕速率之間差異最小,RA=1.845.依據(jù)R可以確定兩個(gè)影響因素中,濕度為主要因素,溫度為次要因素.

由k可以看出,當(dāng)溫度因素取40 ℃以上時(shí),k=2.009~2.701;濕度因素取80%以上時(shí),k=1.869~2.570,對(duì)腐蝕速率的影響平均值均比較大,這兩種條件組合是鋼絲腐蝕最惡劣的環(huán)境.

4.3 試驗(yàn)數(shù)據(jù)方差

通過(guò)極差分析,可以考察因素對(duì)試驗(yàn)結(jié)果影響的重要程度,但是不能估計(jì)試驗(yàn)過(guò)程中以及試驗(yàn)結(jié)果測(cè)量中存在的誤差大小,不能區(qū)分某因素各水平所對(duì)應(yīng)的試驗(yàn)結(jié)果間的差異究竟是由因素水平不同所引起的,還是由于試驗(yàn)誤差引起的,因此對(duì)數(shù)據(jù)進(jìn)行方差分析和假設(shè)檢驗(yàn),計(jì)算試驗(yàn)過(guò)程中的誤差.

方差分析的基本方法是將總的變差平方和分解為因素的變差平方與誤差平方兩個(gè)部分,并計(jì)算出因素及誤差的平均平方和,用F檢驗(yàn)法對(duì)因素進(jìn)行顯著性檢驗(yàn).由于方差分析要求因素之間整齊可比,因此,選擇5個(gè)濕度因素,即

40%~50%, 50%~60%, 70%, 80%, 90%.

溫度值也只選取上述5個(gè)濕度對(duì)應(yīng)溫度工況下的測(cè)量值進(jìn)行方差分析,如表6所示,表中:K1~K5為A、B因素各取值條件下的試驗(yàn)值之和;Q為K1~K5的平方和均值;P為K的平方取均值;S=Q-P為平方和之差;W為A、B因素各取值條件下的試驗(yàn)值平方之和.

計(jì)算原理和方法見文獻(xiàn)[15].

方差分析顯著性如表7所示.表中:e為誤差平方和.

表6 正交試驗(yàn)結(jié)果方差分析計(jì)算Tab.6 Variance analysis table of orthogonal test results

表7 方差分析顯著性計(jì)算Tab.7 Significant calculation table of variance analysis

由表7可知,在可靠性為95%的概率下,A對(duì)腐蝕速率的影響不顯著,B對(duì)腐蝕速率的影響顯著,這種顯著可靠性有95%的概率.

4.4 溫度、濕度與腐蝕速率計(jì)算關(guān)系擬合

使用Matlab對(duì)測(cè)量得到的數(shù)據(jù)進(jìn)行曲面作圖,如圖3所示.

由圖3可知,在濕度大于60%后,曲面變陡,與低濕度情況下的斜率出現(xiàn)較大差異,為了得到更好的擬合公式,將溫度濕度-腐蝕速率譜進(jìn)行分段擬合.使用數(shù)據(jù)統(tǒng)計(jì)軟件,采用最小二乘法回歸分析對(duì)插值后的數(shù)據(jù)進(jìn)行回歸擬合,得到分段計(jì)算式為

x1∈[10,50],

x2∈[40,60],

y=25.549-0.298x1-0.549x2+

3.17×10-3x1x2,

x1∈[10,50],

x2∈[61,90],

式中:

x1為溫度,℃;

x2為相對(duì)濕度,%;

y為未鍍鋅高強(qiáng)鋼絲腐蝕速率,g/(dm2·h).

為了擬合曲面計(jì)算公式,使用樣條插值對(duì)結(jié)果進(jìn)行擬合,得到的擬合曲面如圖4所示.

由回歸方程可以看出,在濕度大于60%以后,二次項(xiàng)的系數(shù)變大,表明在高濕度環(huán)境下,溫度與濕度的關(guān)系存在部分耦合和非線性效應(yīng).

需要注意,由于腐蝕速率與溫度、濕度在物理意義上沒(méi)有關(guān)系,上述計(jì)算公式僅是在擬合的溫度、濕度范圍的統(tǒng)計(jì)回歸公式,沒(méi)有具體的物理意義.

圖3 測(cè)量結(jié)果擬合曲面Fig.3 Fitting surface of experimental data

圖4 插值結(jié)果擬合曲面Fig.4 Fitting surface of interpolation results

擬合結(jié)果的可信度指標(biāo)如表8所示.由表8可以看出:

(1) 復(fù)相關(guān)系數(shù)用來(lái)衡量自變量x1,x2與y之間的相關(guān)程度的大小,本試驗(yàn)中,復(fù)相關(guān)系數(shù)分別為0.952、0.937,表明R1、R2之間高度正相關(guān);

(2) 復(fù)測(cè)定系數(shù)用來(lái)說(shuō)明自變量解釋因變量y變差的程度,以測(cè)定因變量y的擬合效果,復(fù)測(cè)定系數(shù)分別為0.906、0.878,表明可用自變量解釋因變量變差的90.6%、87.8%;

(3) 標(biāo)準(zhǔn)誤差用來(lái)衡量擬合程度的大小,也用于計(jì)算與回歸相關(guān)的其它統(tǒng)計(jì)量,本試驗(yàn)中標(biāo)準(zhǔn)差分別為0.039、0.031,此值越小,說(shuō)明擬合程度越好;

(4)Q用來(lái)描述擬合回歸公式與測(cè)量數(shù)據(jù)之間的偏離程度,試驗(yàn)中分別為0.016、0.009,值越小,說(shuō)明擬合程度越高.

表8 擬合回歸統(tǒng)計(jì)結(jié)果Tab.8 Table of fitting regression statistics

5 結(jié) 論

本文基于正交方法,采用恒溫恒濕箱控制試驗(yàn)環(huán)境溫度和濕度,用極化電阻腐蝕傳感器測(cè)量了未鍍鋅高強(qiáng)鋼絲在不同溫濕組合環(huán)境下的腐蝕速率,以腐蝕速率為1×10-4g/(dm2·h)、長(zhǎng)度為1 m、直徑為5 mm的未鍍鋅鋼絲為例,按表面積上發(fā)生均勻腐蝕計(jì)算,每年的腐蝕深度約為0.01 mm,100 a 后腐蝕深度為1 mm,達(dá)到直徑的20%,若假設(shè)鋼絲的腐蝕僅受到溫度、濕度的影響,則認(rèn)為腐蝕速率小于1×10-4g/(dm2·h)滿足橋梁的設(shè)計(jì)壽命,由此可以將鋼絲的腐蝕環(huán)境區(qū)間劃分為3個(gè)區(qū)域:

(1) 弱腐蝕區(qū),環(huán)境濕度70%、溫度10 ℃以下區(qū)域,該區(qū)域的腐蝕速率很小,數(shù)量級(jí)1×10-6g/(dm2·h)左右,可以認(rèn)為在這一條件下,鋼絲腐蝕速率很小,可以不進(jìn)行腐蝕控制.

(2) 低腐蝕區(qū),環(huán)境濕度75%以下、溫度20~50 ℃區(qū)域,該區(qū)域腐蝕速率變大,數(shù)量級(jí)小于1×10-4g/(dm2·h),在這一環(huán)境下,鋼絲腐蝕速率較低,這一區(qū)域可以為懸索橋主纜抽濕送風(fēng)系統(tǒng)的控制目標(biāo)區(qū)域.

(3) 強(qiáng)腐蝕區(qū),環(huán)境濕度75%以上、溫度10~50 ℃ 區(qū)域,該區(qū)域腐蝕速率大于1×10-4g/(dm2·h),在這一環(huán)境下,鋼絲腐蝕較快,需要預(yù)防主纜處于這一腐蝕環(huán)境區(qū)域,加強(qiáng)主纜的防腐蝕防護(hù)措施.

參考文獻(xiàn):

[1]孫躍,胡津.金屬腐蝕與控制[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,2003: 1-9 .

[2]BETTI R,YANEV B.Corrosion monitoring research of new york city bridges[R].Washington D.C.: Transportation Research Board,2013.

[3]National Cooperative Highway Research Program.NCHRP-Report[R].New York: [s.n.],2004.

[4]陳小雨,沈銳利,張培炎.懸索橋主纜的內(nèi)部檢測(cè)[J].世界橋梁,2013,41(1):86-90.

CHEN Xiaoyu,SHEN Ruili,ZHANG Peiyan.Internal inspection of suspension bridge main cable[J].World Bridges,2013,41(1):86-90.

[5]SUZUMURA K,NAKAMURA S I.Environmental factors affecting corrosion of galvanized steel wires[J].Journal of Materials in Civil Engineering,2004,16:1-7.

[6]黎學(xué)明,陳大華.模擬酸雨溶液中溫度對(duì)橋梁索纜鍍鋅鋼絲腐蝕行為影響[J].腐蝕科學(xué)與防護(hù)技術(shù),2010,22(1):14-17.

LI Xueming,CHEN Dahua.Effect of temperature on corrosion behavior of galvanized steel bridge wires in simulated acid rain[J].Corrosion Science and Protection Technology,2010,22(1):14-17.

[7]楊文靜,黎學(xué)明,周建庭.大跨度橋梁纜索模擬酸雨加速腐蝕行為研究[J].腐蝕科學(xué)與防護(hù)技術(shù),2011,23(1):65-68.

YANG Wenjing,LI Xueming,ZHOU Jianting.Research of accelerated corrosion behavior of long-span bridge cable simulated acid rain[J].Corrosion Science and Protection Technology,2011,23(1):65-68.

[8]BETTI R, MARCONI G, KHAZEM D, et al.Experimental analysis of a nondestructive corrosion monitoring system for main cables of suspension bridges[J].Journal of Bridge Engineering,2013,18(7):653-662.

[9]繆長(zhǎng)青,尉廷華.大跨橋梁纜索鋼絲腐蝕速率的研究[J].西南交通大學(xué)學(xué)報(bào),2014,49(3):513-518.

MIAO Changing,WEI Tinghua.Corrosion rate test of cable wires of large span bridge[J].Journal of Southwest Jiaotong University,2014,49(3):513-518.

[10]CHEN Xiaoyu,TANG Maolin.Influence analysis of the corrosion factors effect on the corrosion rate of steel wires of main cable of suspension bridge[C]∥IABSE Conference.Guangzhou:[s.n.],2016:155-163.

[11]施彥彥.典型金屬材料大氣腐蝕的模擬電化學(xué)研究[D].杭州:浙江大學(xué),2008.

[12]俞明德,沈銳利,唐茂林,等.西堠門大橋主纜橫斷面溫度場(chǎng)研究[J].建筑科學(xué)與工程學(xué)報(bào),2010,27(3):53-58.

YU Mingde,SHEN Ruili,TANG Maolin,et al.Research on temperature field of main cable section of xihoumen bridge[J].Journal of Architecture and Civil Engineering,2010,27(3):53-58.

[13]楊寧.大跨度懸索橋的主纜除濕防護(hù)系統(tǒng)[D].上海:同濟(jì)大學(xué),2006.

[14]DOUGLAS W B,RICHARD J.Linear polarization resistance sensor using the structure as a working electrode[C]∥European Conference of The Prognostics and Health Management Society.Santa Clara:[s.n.],2014:1-7.

[15]何為,薛衛(wèi)東.優(yōu)化試驗(yàn)設(shè)計(jì)方法與數(shù)據(jù)分析[M].北京:化學(xué)工業(yè)出版社,2012:29-284.

猜你喜歡
測(cè)量環(huán)境
長(zhǎng)期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
孕期遠(yuǎn)離容易致畸的環(huán)境
不能改變環(huán)境,那就改變心境
把握四個(gè)“三” 測(cè)量變簡(jiǎn)單
滑動(dòng)摩擦力的測(cè)量和計(jì)算
環(huán)境
孕期遠(yuǎn)離容易致畸的環(huán)境
滑動(dòng)摩擦力的測(cè)量與計(jì)算
測(cè)量的樂(lè)趣
主站蜘蛛池模板: 欧美日韩一区二区三区在线视频| 极品av一区二区| 午夜视频日本| 国产成人精品一区二区| 亚洲成人在线网| 538精品在线观看| 老司机精品99在线播放| 亚洲成人一区二区三区| 青青热久麻豆精品视频在线观看| 在线观看免费国产| 人人妻人人澡人人爽欧美一区| 亚洲第一成年网| 手机在线免费不卡一区二| 成人综合久久综合| 亚洲精品天堂自在久久77| 97成人在线视频| 日本不卡在线播放| 国产精品亚洲五月天高清| 婷婷久久综合九色综合88| 色综合成人| 国产成a人片在线播放| 激情综合婷婷丁香五月尤物| 奇米影视狠狠精品7777| 久草视频福利在线观看| 日韩第一页在线| 99re这里只有国产中文精品国产精品 | 午夜日韩久久影院| 亚洲成综合人影院在院播放| 国产成人夜色91| 欧美一区二区自偷自拍视频| 色网站在线免费观看| 久久9966精品国产免费| 亚洲毛片一级带毛片基地| 国产在线98福利播放视频免费| 精品国产毛片| 久久中文字幕av不卡一区二区| 91精品日韩人妻无码久久| 亚洲精品爱草草视频在线| 亚洲国产成人自拍| 免费人成视网站在线不卡| 亚洲中文字幕97久久精品少妇| 2020精品极品国产色在线观看| 波多野结衣一区二区三区88| 手机在线免费不卡一区二| 日韩精品专区免费无码aⅴ| 国产97视频在线| 天天色天天综合网| 国产成人盗摄精品| 久久久久久久久久国产精品| 亚洲日本在线免费观看| 久久久久免费看成人影片| 久久精品无码一区二区日韩免费| 国产无码精品在线播放| 国产精品v欧美| 美女免费精品高清毛片在线视| 国产亚洲精品91| 伊人成色综合网| 亚洲天堂视频在线播放| 日韩天堂网| 国产男人的天堂| 中文精品久久久久国产网址| 无码精品国产dvd在线观看9久| 91网红精品在线观看| 99人妻碰碰碰久久久久禁片| 欧美色图久久| 亚洲欧美另类色图| 亚洲最猛黑人xxxx黑人猛交| 亚洲无码日韩一区| 精品视频第一页| 国产成人AV综合久久| 欧美午夜视频| 国产在线97| 国产美女免费| 亚洲视频影院| 国产在线欧美| 亚洲中文字幕国产av| 国产精品香蕉| 香蕉久久国产超碰青草| 国产麻豆福利av在线播放| 欧美精品一区在线看| 在线观看91精品国产剧情免费| 亚洲熟妇AV日韩熟妇在线|