亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        往復動載下明挖隧道結(jié)構(gòu)可靠性評價

        2023-11-06 01:44:18張祖迪王玉鎖盧雅欣肖宗揚何鎖宋
        四川建筑 2023年5期
        關鍵詞:抗力劣化剪應力

        張祖迪, 王玉鎖, 盧雅欣, 肖宗揚, 何鎖宋, 張 瑞

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

        0 引 言

        隨著我國經(jīng)濟的發(fā)展,快速交通網(wǎng)絡正處于高速發(fā)展時期,由市政公路、軌道交通及飛機等交通方式組成的復雜城市交通線路在交織時,難免會產(chǎn)生重疊、近接等情況,而無可避免會出現(xiàn)多動載耦合情況,這無疑會給隧道結(jié)構(gòu)安全運營帶來極大挑戰(zhàn),為此,評價多動載往復作用下隧道結(jié)構(gòu)的可靠性,以滿足運營安全需求,是隧道結(jié)構(gòu)設計亟待解決的難題。

        針對運營期間動載作用下隧道結(jié)構(gòu)可靠性問題,國內(nèi)外學者已開展了一定的研究。謝楠等[1]采用隨時間衰減的抗力模型模擬二次襯砌抗力的劣化,采用連續(xù)介質(zhì)模型+響應面法求解荷載效應,用不等寬的脈沖隨機過程近似描述抗力,提出100年設計使用期內(nèi)可靠度的計算公式。楊建國等[2]基于熵權(quán)的模糊物元分析方法,把隧道結(jié)構(gòu)的狀態(tài)級別作為物元的事物,用襯砌裂縫、混凝土強度、襯砌背后空洞等評價指標及其相應的模糊量值構(gòu)造復合模糊物元,通過關聯(lián)度計算,得到服役公路隧道可靠性評價結(jié)果。李小勇[3]基于隧道安全監(jiān)測系統(tǒng)選取5個指標作為城市隧道結(jié)構(gòu)可靠性評價依據(jù),建立了基于單指標可靠性等級評價的多指標綜合評價模型。歐爾峰[4]考慮圍巖荷載和列車振動荷載共同作用,分析研究高速鐵路黃土隧道列車振動響應的極值分布,采用JC法計算分析高速鐵路黃土隧道結(jié)構(gòu)的可靠度指標,來評價其安全性。

        綜上所述,對于運營期間隧道結(jié)構(gòu)可靠性研究,目前已取得了一定的成果,但主要集中在列車荷載作用下的單體隧道結(jié)構(gòu)可靠性研究,對于多種類往復動載作用下多跨復雜斷面明挖隧道結(jié)構(gòu)的長期可靠性評價方法,還缺乏相關研究,而這無疑也給相關地區(qū)地鐵隧道長期運營安全性評估及線路的維保帶來了諸多不便。在此背景下,本文以成都天府國際機場為例,結(jié)合隧道結(jié)構(gòu)特點,提出評價多種類往復動載下結(jié)構(gòu)可靠性評價方法,以滿足運營安全需求。

        1 往復動載下基底累積沉降預測模型

        土體在長期循環(huán)荷載作用下的變形特性是研究土體在多種運行荷載作用下的累積沉降以及動力穩(wěn)定性的基礎。為了研究土體在多種運行荷載作用下的累積沉降規(guī)律,首先必須研究土體在循環(huán)荷載作用下的長期變形特性,這樣才能給實際工程設計提供比較可靠地設計依據(jù)。國內(nèi)外許多學者根據(jù)土體循環(huán)動力特性試驗提出了基于第一次循環(huán)變形的計算模型,最常用的是Monismith[5]在總結(jié)前人研究成果基礎上提出的簡單指數(shù)模型見式(1)。

        εp=A×Nb

        (1)

        式中:εp為累積塑性應變,N為循環(huán)加載次數(shù),A和b為地基特性參數(shù),對于復雜多層土體時,參數(shù)取值離散性大。

        Li 和 Selig[6]改進了對參數(shù)A的確定見式(2)。

        (2)

        于是Li和Selig改進后的累積塑性應變公式為式(3)。

        (3)

        Chai和Miura依據(jù)Samnag試驗結(jié)果,考慮土體地應力影響,假設累積塑性應變隨初始地應力增加而增大,在Li和Selig模型基礎上進行了發(fā)展,提出新的修正公式[7],見式(4)。

        (4)

        式中:qs為靜偏應力,參數(shù)a、b和m與土體的類型和塑性指數(shù)有關,按建議值n取常數(shù)1[7],qf為土體靜止破壞應力,可通過室內(nèi)試驗或原位試驗獲得,也可根據(jù)有效固結(jié)應力理論[9],利用強度指標c和φ計算,見式(5)、式(6)。

        (5)

        K0=1-sinφ

        (6)

        式中:K0為靜止土壓力系數(shù),σ1為第一主應力。

        根據(jù)以上理論,利用分層總和法計算循環(huán)荷載累積塑性變形的具體步驟如下:

        1.1 應力偏量

        根據(jù)不同荷載折減系數(shù)下的數(shù)值模擬結(jié)果,利用Midas后處理提取單元的正應力、剪應力,按式(7)計算應力偏量。

        (7)

        (8)

        1.2 靜偏應力

        利用有限元分析的方法,得到土體在僅存在自重荷載時各單元的應力值,將其代入式(9),計算響應的靜偏應力。

        (9)

        1.3 動偏應力

        利用動力有限元分析方法,得到土體在動荷載作用下各單元對應應力值,然后將其代入式(10),計算響應動偏應力。

        (10)

        1.4 累積沉降

        將式qf、qd代入計算式(4)中計算得到各層土的塑性應變εp。首先按一定的厚度對隧道襯砌下方土體進行分層,假設各層的塑性應變不變;然后計算出各層關鍵單元處的累積塑性應變εp;最后用分層總和法計算土體的長期累積沉降,見式(11)[6]。

        (11)

        式中:εpi為第i層土體的累積塑性應變;Hi為第i層土的厚度,共n層土。

        2 可靠性評價方法

        2.1 非概率可靠性指標計算方法

        多種類往復動載作用下的隧道結(jié)構(gòu),所涉及的動載速度、基底材料的力學參數(shù)及其在動載作用下隨時間劣化無明顯統(tǒng)計特征規(guī)律,在此條件下,傳統(tǒng)的概率可靠性計算方法難以實現(xiàn)結(jié)構(gòu)可靠性分析,因此引入非概率可靠性計算方法,基于非概率可靠性計算原理,將動載速度及土體力學參數(shù)看作均勻分布,定義不確定參數(shù)區(qū)間,研究結(jié)構(gòu)可靠性。

        區(qū)間分析法可以說是凸集合的一種特殊情況[8]。假定在實數(shù)域R內(nèi),對于任意給定的2個實數(shù)x1,xu存在以下關系:x1,xu∈R,且x1

        xI=[x1,xu]={x∈R|xl≤x≤xu}

        (12)

        式中:x為區(qū)間變量;xI為區(qū)間變量x所在區(qū)間;x1,xu為區(qū)間變量x的上下界。

        當功能函數(shù)中自變量是有界變量,且功能函數(shù)Z是連續(xù)函數(shù)時,那么Z也是一區(qū)間變量,由此可以求得zmin、zmax見式(13)、式(14)。

        zmin=R-Fmin

        (13)

        zmax=R-Fmax

        (14)

        設變量Z均值和離差分別為Rc和RΓ,可求得其對應的均值Rc和離差RΓ見式(15)、式(16)。

        (15)

        (16)

        代入式(17),可求得結(jié)構(gòu)非概率可靠指標μ。

        (17)

        按照結(jié)構(gòu)可靠性理論,當Z=0時結(jié)構(gòu)失穩(wěn),則:當μ>1時,則對于任意自變量均有Z>0,結(jié)構(gòu)安全;當μ<1時,則對于任意自變量均有Z<0,結(jié)構(gòu)失效;當-1≤μ≤1時,Z>0和Z<0的情況均存在,此時結(jié)構(gòu)狀態(tài)無法確定。

        2.2 功能函數(shù)的建立

        由第三強度理論可知,剪應力是綜合考慮第一、第三主應力的結(jié)果,且最大剪應力是引起材料屈服的原因,因此對于明挖復雜隧道結(jié)構(gòu)的可靠性分析以結(jié)構(gòu)所受剪應力為判斷依據(jù)。由建立的結(jié)構(gòu)最大剪應力τmax與其他基本變量的顯示公式,可得到其功能函數(shù)Z,并求得Z的取值范圍Y=[YL=zmin,YR=zmax],則由該功能函數(shù)的可靠性分析,得到3種情況,如圖1所示。

        根據(jù)第2節(jié)可知,在長期荷載作用下基底的累積塑性應變滿足εp=A×Nb的關系,但考慮到荷載的長期作用可能會導致基底材料的劣化,故需要對上式進行修正。考慮到材料的劣化是由多種因素公共作用的結(jié)果,本文參考強度折減法的原理,模擬材料的劣化,其中c、φ值的劣化系數(shù)為F1,彈性模量E的劣化系數(shù)為F2,即式(18)~式(20)。

        (18)

        (19)

        (20)

        (21)

        (22)

        由于基底的沉降會導致結(jié)構(gòu)變形,從而引起結(jié)構(gòu)最大剪應力變化,本節(jié)擬采用基底沉降作為結(jié)構(gòu)在長期荷載作用下的豎向變形值,故結(jié)構(gòu)最大剪應力τmax與F1,F2也成一定的數(shù)學關系,即往復動載下隧道結(jié)構(gòu)可靠度功能函數(shù)為式(23)。

        τmax=f(F1,F2,t)

        (23)

        2.3 計算流程

        首先確定劣化系數(shù),建立多種類往復動載作用下的數(shù)值模型,計算不同劣化系數(shù)下土體的累積沉降,將土體沉降看作隧道結(jié)構(gòu)在動載作用下的變形,對隧道結(jié)構(gòu)施加強制位移,模擬由于變形而導致的隧道結(jié)構(gòu)受力,利用基于凸集理論的穩(wěn)健可靠性原理,引入非概率可靠指標,提出一套考慮長期往復動載作用的隧道結(jié)構(gòu)可靠性評價的方法(圖2)。

        圖2 計算流程

        3 工程應用

        3.1 工程概況

        成都天府國際機場北垂滑隧道群包括3條飛機滑行道,9條下穿機場滑行道明挖隧道,形成了以下穿機場北側(cè)滑行道的地鐵13號、18號線及機場輕型自動導向軌道交通(PRT)的五跨及中跨兩層的“品”字形明挖隧道結(jié)構(gòu)(圖3)。

        圖3 復雜斷面明挖隧道

        3.2 數(shù)值模擬

        3.2.1 動載激勵的模擬

        3.2.1.1 地鐵列車動載激勵的模擬

        根據(jù)調(diào)研相關文獻的做法[10],參考英國德比鐵道技術中心根據(jù)試驗及理論研究成果,用一個激振力函數(shù)來模擬列車荷載,其表達式為式(24)。

        F(t)=p0+p1sinω1t+p2sinω2t+p3sinω3t

        (24)

        根據(jù)國內(nèi)地鐵列車行駛速度,行車速度v取100 km/h,列車簧下質(zhì)量M0取75 kN·s2/m,列車靜載p0取80 kN,可得地鐵列車運行荷載時程曲線如圖4所示。

        圖4 地鐵列車運行荷載時程曲線

        3.2.1.2 PRT列車動載激勵的模擬

        對于PRT列車,由于相關荷載激勵的模擬較少,其荷載作用方式同汽車荷載類似,可將車輛荷載分解成一系列的正弦波荷載[11],其表達式為式(25)。

        P(t)=p0+psinωt

        (25)

        式中:P(t)為PRT動載(kN);p0為PRT靜載(kN);p為振動荷載幅值,p=M0aω2;M0為汽車模型簧下質(zhì)量(kN·s2/m);a為幾何不平順矢高,反映了路面的不平整度狀況,a取2 mm;ω為振動圓頻率,ω=2πv/l;v為車速(km/h);l為路面幾何曲線波長(m)。

        PRT列車靜載取p04.73 kN,簧下質(zhì)量M0取60 kN·s2/m,車速v取40 km/h,路面幾何曲線波長l取3.7 m。得到PRT列車運行荷載時程曲線,見圖5。

        圖5 PRT列車運行荷載時程曲線

        3.2.1.3 飛機動載激勵的模擬

        以A380載客飛機為主要研究對象,飛機產(chǎn)生的最不利動荷載采用相關文獻[9]計算公式進行計算,飛機前輪動荷載F1=84.3 kN,飛機后輪動荷載F2=272.6 kN。

        F11=F1+10%F1·sin(ωt)=84.3+8.43sin(22t)kN

        (26)

        F22=F2+10%F2·sin(ωt)=272.6+87.26sin(22t)kN

        (27)

        可得到前后輪的動載時程曲線如圖6所示。

        圖6 飛機動荷載時程

        3.2.2 模型建立

        本節(jié)基于成都天府國際機場下穿機場滑行道多洞并行明挖隧道群建立數(shù)值模型,將3.2.1節(jié)確定的動載以集中力的形式,通過控制荷載的到達時間和作用時間,施加在鋼軌及機場跑道相應位置處,通過定義接觸界面的法向剛度模量kn、切向剛度模量kt來模擬隧道襯砌與軌道板之間的接觸界面,并在模型四周添加粘彈性邊界,建立數(shù)值模型,如圖7所示。

        圖7 數(shù)值模型

        下穿機場滑行道“品”字形明挖隧道結(jié)構(gòu)復雜、橫向多跨、中部上下重疊,在地鐵列車、PRT以及飛機多動載耦合作用下,結(jié)構(gòu)及基礎受力機理更為復雜,因此本次研究主要以此斷面為對象進行分析,如圖8所示。

        圖8 “品”字形隧道結(jié)構(gòu)模型細部

        數(shù)值模擬所涉及的相關材料參數(shù)如表1所示。

        表1 材料參數(shù)

        3.3 動載組合對“品”字形隧道結(jié)構(gòu)力學響應影響

        在實際運營過程中,結(jié)構(gòu)所受往往是多種動載組合作用,通過對PRT、地鐵、飛機動載的兩兩或者三者組合,荷載組合工況如表2所示,來研究不同動載組合下“品”字形隧道結(jié)構(gòu)最大剪應力響應。

        表2 交通荷載組合

        通過數(shù)值模擬,得到隧道結(jié)構(gòu)在不同荷載組合下的最大剪應力值,如表3所示。

        表3 不同組合下結(jié)構(gòu)最大剪應力 單位:MPa

        對比可知,組合3下結(jié)構(gòu)所受剪應力最大,為1.17 MPa,故組合3為最不利荷載組合。

        3.4 累積沉降計算

        由于實際土體材料的劣化過程是一個多因素耦合作用的結(jié)果,本次研究考慮采用區(qū)間分析的理論,劣化系數(shù)F1、F2的取值區(qū)間參考相關文獻[13-14],本次數(shù)值模擬中的取值F=[1,1.3],通過控制劣化系數(shù)上下限,來保證土體的實際劣化處于假定區(qū)間(圖9)。

        圖9 土體材料劣化示意

        通過對劣化系數(shù)F1、F2進行不同組合(表4),計算考慮地基材料劣化條件下隧道結(jié)構(gòu)可靠度。

        以組合5為例進行計算,基于第3節(jié)所述往復動載作用下明挖回填土體累積沉降預測模型,研究隧道在PRT+地鐵列車+飛機(圖10)長期累積荷載作用下的豎向變形值。

        圖10 數(shù)值模型

        表5 “品”字形隧道特征點土體累積塑性變形 單位:mm

        (1)隧道襯砌下方土體累積沉降隨深度的增加而逐漸減小,在深度為10 m時,累積沉降值趨近于0 mm。

        (2)隨著時間增長,各監(jiān)測點的沉降速率呈現(xiàn)先是增加較快,之后逐漸趨于平穩(wěn)規(guī)律。

        (3)結(jié)構(gòu)兩側(cè)土體的累積沉降值明顯大于結(jié)構(gòu)中部土體的累積沉降值。

        3.5 功能函數(shù)的建立

        基于3.4節(jié)研究在不同折減系數(shù)以及荷載作用時間下土體的豎向變形,利用有限元軟件,對隧道襯砌施加強制位移,如圖12所示,強制位移量取各測點下方土體的沉降量。由此可以得到了不同折減系數(shù)下結(jié)構(gòu)最大剪應力隨時間變化的曲線,如圖13所示。

        圖12 施加強制位移

        圖13 變形導致的最大剪應力示意

        用同樣的方法確定不同劣化系數(shù)組合下隧道襯砌各測點隨時間變化的最大剪應力,如圖14所示。

        圖14 不同劣化系數(shù)下最大剪應力

        在不同劣化系數(shù)下,結(jié)構(gòu)所受的最大剪應力不相等,且劣化程度越高,結(jié)構(gòu)所受最大剪應力值越大。隨著時間增長,結(jié)構(gòu)的剪應力增大速率呈現(xiàn)先是增加較快,之后逐漸趨于平穩(wěn)規(guī)律。

        通過線性回歸擬合圖14中的最大剪應力變化曲線,得到各曲線對應的數(shù)學表達式,如表6所示。

        表6 不同劣化系數(shù)組合下的功能函數(shù)

        表7 可靠指標μ隨時間變化

        A=-1.217F1-4.750F2+5.142

        (28)

        B=1.417F1+5.767F2-6.185

        (29)

        1.417F1+5.76F2-6.185

        (30)

        式中:τmax為最大剪應力(MPa);F1為c、φ值的劣化系數(shù);F2為彈性模量E值的劣化系數(shù);t為荷載作用時間(a)。

        運營期隧道結(jié)構(gòu)的可靠度分析需考慮抗力衰減引起的時變特性。目前,結(jié)構(gòu)抗力的時變模型常表示為初始抗力與衰減函數(shù)的乘積[17]。結(jié)構(gòu)抗力的隨機過程可表示為式(31)。

        R=R0φ(t)

        (31)

        式中:R0為抗力初始值,φ(t)為確定性衰減函數(shù)。

        在針對隧道襯砌結(jié)構(gòu)的抗力進行時變衰減分析時,衰減函數(shù)的選擇是模型分析的關鍵。參考相關文獻[16],衰減函數(shù)為式(32)。

        φ(t)=1-k1t+k2t2

        (32)

        式中:k1、k2為衰減系數(shù),取值與抗力退化速率有關,隧道襯砌的抗力退化速率屬于中等退化情況,取參數(shù)k1=0.0005,k2=0[17]。

        綜上,參考本文第2節(jié)非概率可靠性計算方法,建立功能函數(shù)式(33)。

        Z=R-S=[τ]·φ(t)-τmax=1.35(1-0.0005t)+

        1.417F1-5.767F2+6.185

        (33)

        當劣化系數(shù)區(qū)間取F=[1,1.3]時,代入式(17)得結(jié)構(gòu)非概率可靠指標μ隨時間變化情況如表7所示。

        擬合可靠指標μ與時間t的曲線如圖17所示。

        圖17 可靠指標μ隨時間變化

        從圖17可知,考慮隧道底部地基材料參數(shù)劣化,結(jié)構(gòu)可靠指μ隨荷載作用時間的增加呈現(xiàn)逐漸降低的趨勢,且可靠指標μ降低速率呈現(xiàn)先是降低較快,之后逐漸趨于平穩(wěn)規(guī)律,在前30年,可靠指標急劇降低,之后降低速率逐漸減緩;同時可以看出在累積荷載作用約42年后,結(jié)構(gòu)可靠指標小于1,故在考慮土體材料劣化的條件下,當劣化系數(shù)區(qū)間取F=[1,1.3]時,在累積荷載作用約42年后結(jié)構(gòu)可能發(fā)生失效風險,此時還未達到設計使用年限,因此有必要對結(jié)構(gòu)進行優(yōu)化,以保證結(jié)構(gòu)在運營期間的安全。

        3.6 優(yōu)化措施

        為了保證結(jié)構(gòu)在長期累積荷載作用下的運營安全,有必要對結(jié)構(gòu)采取一定的優(yōu)化措施,減小結(jié)構(gòu)在動載作用下的受力變形,具體優(yōu)化方案(圖18):

        圖18 數(shù)值模型

        方案1:在PRT通道下方增設底板。

        方案2:機械可進入范圍內(nèi)采用土石回填,其余部分采用由泡沫混凝土對PRT空腔實施填充。

        方案3:在方案2基礎上增設底板。

        方案4:在結(jié)構(gòu)下方布設樁基,樁基長15.6 m。

        根據(jù)3.5節(jié)研究可知,原方案下結(jié)構(gòu)在運營42年后會發(fā)生失效風險,因此,本節(jié)對提出的4種優(yōu)化方案下結(jié)構(gòu)在長期累積荷載作用下的可靠性進行研究,按照計算流程計算得出劣化系數(shù)F=[1,1.3]時可靠指標μ隨時間變化情況如表8所示。擬合可靠指標μ與時間t的曲線,如圖19所示。

        表8 F=[1,1.3]時可靠指標μ隨時間變化

        圖19 可靠指標μ隨時間變化

        根據(jù)圖19可以看出,與原方案相比,方案1~方案3對延長結(jié)構(gòu)在往復動載下的安全運營年限影響不大,而方案4條件下(圖18(e)),在累積荷載作用約100年后,結(jié)構(gòu)可靠指標仍大于1,故在考慮土體材料劣化的條件下,當劣化系數(shù)區(qū)間取F=[1,1.3]時,在結(jié)構(gòu)設計使用年限100年內(nèi),結(jié)構(gòu)都不存在失效風險,因此方案4可以有效延長結(jié)構(gòu)的安全運營年限。

        4 結(jié)論

        基于非概率可靠性計算原理,建立考慮回填材料劣化的隧道結(jié)構(gòu)最大剪應力預測模型,提出了一套考慮長期往復動載作用的隧道結(jié)構(gòu)可靠性評價新方法,并應用于成都天府國際機場復雜斷面明挖隧道結(jié)構(gòu)可靠性評價中,得到主要結(jié)論:

        (1)針對最不利荷載組合,考慮隧道底部地基材料劣化,對往復動載下結(jié)構(gòu)底部地基沉降隨時間變化的規(guī)律進行了分析,發(fā)現(xiàn)結(jié)構(gòu)兩側(cè)累積沉降明顯大于結(jié)構(gòu)中部,同時,累積沉降隨時間逐漸增大,10年內(nèi)增長速率較快后趨于平穩(wěn)。

        (2)利用強制位移法,對地基差異沉降引起的結(jié)構(gòu)受力進行了分析,發(fā)現(xiàn)結(jié)構(gòu)所受最大剪應力增長規(guī)律與累積沉降變化趨勢一致。

        (3)通過研究結(jié)構(gòu)可靠指標與荷載時間作用時間的關系,發(fā)現(xiàn)在考慮土體材料劣化以及襯砌材料結(jié)構(gòu)抗力衰減的條件下,結(jié)構(gòu)可靠指標隨時間的增加而逐漸降低,呈現(xiàn)先降低較快,后隨著時間的增加逐漸趨于平穩(wěn)。

        (4)通過在“品”字形隧道結(jié)構(gòu)下方布設樁基,可以有效延長結(jié)構(gòu)的安全運營年限,在結(jié)構(gòu)設計使用年限100年內(nèi),結(jié)構(gòu)都不存在失效風險,故在類似工程中可以參考該優(yōu)化方案對結(jié)構(gòu)進行施工,以保證結(jié)構(gòu)的運營安全。

        猜你喜歡
        抗力劣化剪應力
        橋梁低樁承臺基底豎向抗力作用效應的計算研究
        交通科技(2021年4期)2021-09-03 09:47:22
        變截面波形鋼腹板組合箱梁的剪應力計算分析
        基于S形試件五軸數(shù)控機床動態(tài)性能劣化評價
        跨海橋梁混凝土劣化分析及預養(yǎng)護
        上海公路(2017年1期)2017-07-21 13:38:33
        巖塊的彈性模量及巖體單位彈性抗力系數(shù)的確定方法
        八鋼2500m3高爐適應焦炭質(zhì)量劣化的生產(chǎn)實踐
        新疆鋼鐵(2016年3期)2016-02-28 19:18:52
        30MnSi熱變形抗力的數(shù)學模型
        上海金屬(2014年2期)2014-12-18 06:52:51
        裂紋齒輪嚙合剛度的劣化特性
        瀝青路面最大剪應力分析
        河南科技(2014年13期)2014-02-27 14:11:25
        復合式路面層間最大剪應力影響因素研究
        亚洲综合新区一区二区| 天天天综合网| 久久天天躁狠狠躁夜夜中文字幕| 亚洲福利网站在线一区不卡| 我和丰满妇女激情视频| 性色av无码中文av有码vr| 91国视频| 蜜桃av噜噜一区二区三区香| 99久久精品国产91| 99久久免费只有精品国产| 成人无码区免费a片www| 无码专区亚洲avl| 高清国产国产精品三级国产av| 精品视频一区二区三区在线观看 | 小sao货水好多真紧h视频| 国产高清在线精品一区αpp| 在线亚洲妇色中文色综合| 亚洲综合欧美色五月俺也去| 天堂中文在线资源| 国产三级黄色在线观看| 永久免费看黄网站性色| 国产一区二区三区在线观看黄| 日本三级片在线观看| 国产亚洲精品久久久久久| 国产目拍亚洲精品一区二区| 国产色视频在线观看了| 美女不带套日出白浆免费视频| 中文字幕亚洲无线码| 91久久精品一二三区蜜桃| 日本久久久免费观看视频| 国产裸体xxxx视频在线播放| 欧美视频九九一区二区| 亚洲av区一区二区三区| 在线播放五十路熟妇| 粗一硬一长一进一爽一a级| 国产一级r片内射免费视频| 国产av剧情刺激对白| 国产亚洲精品aaaa片小说| 亚洲一区二区在线视频播放| 少妇被粗大猛进进出出男女片 | 久久精品国产亚洲av影院毛片|