欧美日本高清一区_成人福利网站在线观看11_欧美日本在线视频中文字字幕_国产在线精品一区免费香蕉

高架軌道交通引起環(huán)境振動(dòng)的數(shù)值模擬

   2006-04-29 中國(guó)路橋網(wǎng) 佚名 8410
高架軌道交通引起環(huán)境振動(dòng)的實(shí)測(cè)與數(shù)值模擬摘要:從現(xiàn)場(chǎng)實(shí)測(cè)和數(shù)值模擬2個(gè)方面分析研究上海高架軌道交通明珠線引起的環(huán)境振動(dòng)實(shí)況、發(fā)生機(jī)理及振動(dòng)傳播規(guī)律.根據(jù)實(shí)測(cè)資料提出了預(yù)測(cè)明珠線引起環(huán)境振動(dòng)振級(jí)的統(tǒng)計(jì)回歸公式.建立了車-橋子結(jié)構(gòu)和框架-地基基礎(chǔ)子結(jié)構(gòu)2個(gè)體系的分析模型及計(jì)算方法.計(jì)算了典型區(qū)段在列車通過(guò)時(shí)引起的沿線環(huán)境振動(dòng).數(shù)值模擬分析結(jié)果表明,用提出的分析模型和計(jì)算方法算得的環(huán)境振動(dòng)振級(jí)與統(tǒng)計(jì)回歸公式符合較好關(guān)鍵詞:城市高架軌道交通;振動(dòng)實(shí)測(cè);車-橋體系振動(dòng)分析;土-結(jié)構(gòu)動(dòng)相互作用 我國(guó)第一條城市高架軌道交通線———明珠線已在上海建成并投入運(yùn)行.城市軌道交通引起的環(huán)境振動(dòng)是一種新的城市公害.城市的環(huán)境振動(dòng)除了會(huì)影響人體健康,干擾人們的正常生活及降低工作效率外,還會(huì)影響精密設(shè)備和儀器的正常使用,嚴(yán)重者甚至?xí)?duì)老舊建筑物造成損傷. 城市軌道交通的建設(shè)在我國(guó)起步較晚,故關(guān)于城市軌道交通引起環(huán)境振動(dòng)的實(shí)測(cè)和理論分析研究較少.20世紀(jì)80年代中期,我國(guó)曾對(duì)4個(gè)城市的鐵路交通引起的環(huán)境振動(dòng)實(shí)況進(jìn)行過(guò)實(shí)測(cè)[1],并于1988年制定了《城市區(qū)域環(huán)境振動(dòng)標(biāo)準(zhǔn)》[2].90年代末,在我國(guó)城市軌道交通進(jìn)入快速發(fā)展階段后,北方交通大學(xué)等利用二維有限元模型對(duì)高架軌道交通引起的環(huán)境振動(dòng)進(jìn)行過(guò)理論分析[3,4],但無(wú)實(shí)測(cè)資料印證. 筆者從現(xiàn)場(chǎng)實(shí)測(cè)和數(shù)值模擬分析2個(gè)方面研究上海高架軌道交通———明珠線引起的環(huán)境振動(dòng)的實(shí)況、發(fā)生機(jī)理及振動(dòng)傳播規(guī)律.1高架軌道交通引起環(huán)境振動(dòng)的現(xiàn)場(chǎng)實(shí)測(cè) 在明珠線沿線選擇若干振動(dòng)觀測(cè)場(chǎng)地及與線路垂直的觀測(cè)線,在列車通過(guò)時(shí)對(duì)各觀測(cè)線上諸觀測(cè)點(diǎn)的地面振動(dòng)進(jìn)行了實(shí)測(cè).列車通過(guò)的車速約為60km·h-1.觀測(cè)儀器為愛(ài)華電子研究所生產(chǎn)的AWA6256A型壓電晶體式環(huán)境振動(dòng)分析儀,使用前與AltusK2型強(qiáng)震儀進(jìn)行了同時(shí)觀測(cè)校驗(yàn),驗(yàn)證了儀器的可靠性.1.1沿線地面振動(dòng)的實(shí)測(cè)結(jié)果 振動(dòng)加速度測(cè)量的典型結(jié)果參見(jiàn)圖1.圖1a表示距離明珠線中心0,6,10,16,20m處地面的豎向加速度時(shí)程.圖1b表示圖1a加速度時(shí)程的傅立葉譜.圖1c表示圖1a加速度時(shí)程的振級(jí)(V,dB)的頻率分布,加速度振級(jí)的定義根據(jù)我國(guó)《城市區(qū)域環(huán)境振動(dòng)標(biāo)準(zhǔn)》[2]為式中:αω為頻率計(jì)權(quán)加速度有效值,m·s-2;α0為基準(zhǔn)加速度值,10-6m·s-2. 由圖1可見(jiàn),觀測(cè)線上各測(cè)點(diǎn)的振動(dòng)加速度、加速度時(shí)程的傅里葉譜和1/3倍頻程中心頻率處的振級(jí)隨著與中心線距離的增加而衰減,且當(dāng)頻率愈高時(shí)衰減愈快.圖1c還給出了列車通過(guò)前業(yè)已存在的環(huán)境振動(dòng),即高架軌道交通以外的其他振源引起的振動(dòng)(本文稱本底振動(dòng))對(duì)總的振動(dòng)振級(jí)的影響.由圖可見(jiàn),當(dāng)頻率高于20Hz或距離R小于10m時(shí),本底振動(dòng)的影響很小.這是由于當(dāng)不存在特定振源時(shí),本底振動(dòng)的主要頻率范圍低于20Hz且強(qiáng)度不大.當(dāng)有必要考慮本底振動(dòng)影響評(píng)價(jià)僅由高架軌道交通引起的環(huán)境振動(dòng)時(shí),在距離線路中心10~20m范圍內(nèi)可根據(jù)實(shí)測(cè)結(jié)果將測(cè)得的總振級(jí)減少0.5~3dB[5].圖2表示對(duì)應(yīng)于不同距離測(cè)點(diǎn)1/3倍頻程中心頻率處的振級(jí)的分布情況.實(shí)測(cè)表明:頻率小于16Hz的振動(dòng)隨距離的衰減很小,其倍距離衰減量約為4dB.20~25Hz的振動(dòng)隨距離的衰減較大,其倍距離衰減量約為8dB.而64Hz和80Hz的振動(dòng)隨距離衰減更大,其倍距離衰減量約為8~16dB.1.2振級(jí)V的統(tǒng)計(jì)回歸經(jīng)驗(yàn)公式 根據(jù)明珠線沿線約10個(gè)觀測(cè)區(qū)段的振動(dòng)實(shí)測(cè)結(jié)果,振級(jí)V與距明珠線中心距離R的對(duì)數(shù)值logR的關(guān)系接近于線性關(guān)系.參照鮑尼茲[6]關(guān)于彈性半空間土體在表面擾力作用下地表振動(dòng)衰減的一般公式及日本國(guó)鐵技術(shù)研究所用于預(yù)測(cè)鐵路高架區(qū)間振級(jí)的經(jīng)驗(yàn)公式[7],對(duì)實(shí)測(cè)數(shù)據(jù)利用回歸方法進(jìn)行擬合,得到明珠線高架軌道交通沿線地面振動(dòng)振級(jí)經(jīng)驗(yàn)公式為式中:R為觀測(cè)點(diǎn)與線路中心的距離,m.式(2)與實(shí)測(cè)結(jié)果的對(duì)比參見(jiàn)圖3.由圖可見(jiàn),二者符合良好.2 高架軌道交通引起環(huán)境振動(dòng)的數(shù)值模擬分析 高架軌道交通引起環(huán)境振動(dòng)的數(shù)值分析模型是一復(fù)雜的系統(tǒng),本文將總體分析模型分解為如圖4所示的2個(gè)子結(jié)構(gòu)分別求解.2.1總體分析模型的分解 如圖4所示,將總體分析模型分解為車-橋體系的子結(jié)構(gòu)1和框架-樁基-土體系的子結(jié)構(gòu)2,由子結(jié)構(gòu)1可算得列車過(guò)橋時(shí)作用在框架頂部的動(dòng)支座反力,將此動(dòng)反力作用在子結(jié)構(gòu)2上并考慮框架結(jié)構(gòu)-地基基礎(chǔ)的動(dòng)相互作用,即可算得高架線路沿線任意位置地面的振動(dòng).2.2 子結(jié)構(gòu)的分析模型 在圖5所示的車-橋體系子結(jié)構(gòu)1的振動(dòng)分析中,考慮了列車軌道不平順及梁支座彈性的影響.如圖所示,每一節(jié)車廂假設(shè)為剛體并有2個(gè)輪對(duì)、支承在轉(zhuǎn)向架上由二系彈簧合并成的一系彈簧上.設(shè)車廂具有豎向沉浮運(yùn)動(dòng)Z和回轉(zhuǎn)點(diǎn)頭運(yùn)動(dòng)θ這2個(gè)自由度.設(shè)轉(zhuǎn)向架的間距為S,第i個(gè)車輪至第1個(gè)車輪的距離為bi,列車速度為v,考慮軌道不平順的條件建立車體運(yùn)動(dòng)與轉(zhuǎn)向架運(yùn)動(dòng)的關(guān)系.本文采用頻域法(FDM)擬合軌道不平順功率譜密度函數(shù)生成軌道不平順的樣本[5]. 考慮橋梁支座下設(shè)置的減振橡膠墊及框架、樁基和地基土的彈性變形,可建立彈性支座上梁的動(dòng)平衡方程.利用彈性支承梁的基本振型代替剛性支承梁的基本振型,可求解車-橋結(jié)構(gòu)在列車通過(guò)時(shí)的響應(yīng),進(jìn)而可算得支座動(dòng)反力[5]. 圖6所示的框架-樁基-土體系子結(jié)構(gòu)2在支座動(dòng)反力作用下的地面響應(yīng)可用有限元法計(jì)算.此時(shí)必須在土體有限元模型側(cè)面和底面邊界設(shè)置粘性邊界單元,以消除有限元模型邊界對(duì)土體內(nèi)傳播波的反射.利用粘性邊界單元計(jì)算了下述土-結(jié)構(gòu)動(dòng)相互作用經(jīng)典問(wèn)題:豎向集中簡(jiǎn)諧荷載作用下地面位移響應(yīng),地表和埋置基礎(chǔ)的豎向阻抗函數(shù),樁基的豎向阻抗函數(shù)等,并將其與相應(yīng)的解析解或其他數(shù)值解進(jìn)行對(duì)比,確認(rèn)了粘性邊界單元的有效性[5].3數(shù)值模擬計(jì)算結(jié)果及分析 組合上述2個(gè)子結(jié)構(gòu),即可計(jì)算明珠線高架軌道交通引起的沿線觀察點(diǎn)處的地面振動(dòng).3.1數(shù)值計(jì)算的模型及分析條件 圖6表示列車通過(guò)3個(gè)支承框架時(shí)的框架-樁基-地基分析模型,為進(jìn)行對(duì)比還分析計(jì)算了列車通過(guò)單個(gè)支承框架時(shí)的情況.圖中給出了有限元的分割情況,計(jì)算中利用了問(wèn)題的對(duì)稱性.標(biāo)號(hào)C30的框架混凝土柱截面為1.2m×1.4m、高6m.鋼筋混凝土灌注樁直徑為800mm、樁長(zhǎng)37m、樁的間距為2m,承臺(tái)尺寸為6m×8m×2m、下設(shè)12根樁,承臺(tái)及樁的混凝土標(biāo)號(hào)為C25.場(chǎng)地土特性根據(jù)上海典型層狀地基剖面[8]取用.在有限元模型中,承臺(tái)及土采用三維實(shí)體單元,框架柱及樁采用梁?jiǎn)卧?動(dòng)力荷載P1(t),P2(t),P3(t)是根據(jù)6車廂的列車通過(guò)5跨橋梁的條件算得的3個(gè)相鄰支座的動(dòng)反力,它們互有相位差.對(duì)單框架的計(jì)算模型,僅取單個(gè)支座的動(dòng)反力.明珠線典型區(qū)段的鋼筋混凝土預(yù)應(yīng)力雙孔箱型簡(jiǎn)支梁參數(shù)如下:跨度30m、C50混凝土、單位長(zhǎng)度的質(zhì)量為17980kg·m-1、抗彎剛度為100GPa、阻尼比為0.05.列車參數(shù)如下:各輪對(duì)質(zhì)量為2808kg,空車質(zhì)量24t、轉(zhuǎn)動(dòng)慣量1825.2t·m2、載重12.3t(考慮定員的一半).二系彈簧合并后的彈簧剛度為684.909kN·m-1,阻尼系數(shù)314kN·s·m-1,車廂的長(zhǎng)度為22.8m、轉(zhuǎn)向架的間距為15.7m、輪距為2.5m,車速取60km·h-1.考慮到高架軌道的橋梁結(jié)構(gòu)與鐵路橋梁結(jié)構(gòu)相似,計(jì)算中采用鐵路預(yù)應(yīng)力混凝土橋梁上的軌面不平順功率譜密度函數(shù)[9]作為高架軌道交通軌道不平順功率譜密度函數(shù).3.2數(shù)值計(jì)算與實(shí)測(cè)結(jié)果的對(duì)比 圖7a和b分別表示根據(jù)單框架和三框架模型算得的與線路中心線垂直線上不同距離各觀察點(diǎn)處地面振動(dòng)的振級(jí),并與振級(jí)經(jīng)驗(yàn)公式(2)進(jìn)行比較.圖形顯示,根據(jù)三框架模型算得的結(jié)果比單框架模型更接近于公式(2).這是由于三框架模型考慮了相鄰框架振動(dòng)引起的波動(dòng)對(duì)地面振動(dòng)的影響.當(dāng)用單框架模型評(píng)價(jià)離開(kāi)線路中心線10m以外觀察點(diǎn)的振級(jí)時(shí),計(jì)算結(jié)果與經(jīng)驗(yàn)公式的誤差加大.圖8a和b分別表示單框架和三框架模型算得的1/3倍頻程中心頻率處的振級(jí)隨距離的衰減.圖形顯示,三框架模型的計(jì)算結(jié)果普遍大于單框架模型計(jì)算結(jié)果,且算得的沿線地面振動(dòng)的振級(jí)與實(shí)測(cè)結(jié)果接近.因此本文建議的數(shù)值分析模型可應(yīng)用于預(yù)測(cè)高架軌道交通引起的環(huán)境振動(dòng)振級(jí).3.3 車速和載客量對(duì)環(huán)境振動(dòng)的影響利用以上三框架計(jì)算模型,分析計(jì)算了列車運(yùn)行速度對(duì)沿線環(huán)境振動(dòng)的影響.計(jì)算了額定車速60km·h-1和假定車速80,100km·h-1時(shí)沿線振動(dòng)振級(jí).圖9表示車速為60,80,100km·h-1時(shí)與線路中心線不同距離各觀測(cè)點(diǎn)處的振級(jí)隨距離的衰減情況.由圖可見(jiàn),隨著車速的提高,振級(jí)不斷加大,車速每增加20km·h-1,振級(jí)約增加4dB.利用同樣的計(jì)算模型,分析計(jì)算了列車載客量對(duì)沿線環(huán)境振動(dòng)的影響.計(jì)算了定員半數(shù)載客量12.3t和滿員載客量24.6t兩種情況下的沿線振動(dòng)振級(jí).圖10給出了半載和滿載時(shí)與線路中心線不同距離各觀測(cè)點(diǎn)處的振級(jí)隨距離的衰減情況.由圖可見(jiàn),列車載客量對(duì)振級(jí)的影響很小.4結(jié)論 上海市高架軌道交通明珠線引起環(huán)境振動(dòng)的現(xiàn)場(chǎng)實(shí)測(cè)和數(shù)值分析的研究結(jié)論如下:(1)高架軌道交通沿線環(huán)境振動(dòng)的振級(jí)隨著與線路中心線距離的增大而衰減,振動(dòng)高頻成分的衰減比低頻更快.(2)根據(jù)實(shí)測(cè)結(jié)果的統(tǒng)計(jì)回歸分析,得出了上海高架軌道交通沿線環(huán)境振動(dòng)振級(jí)的經(jīng)驗(yàn)公式,經(jīng)驗(yàn)公式與實(shí)測(cè)數(shù)據(jù)吻合良好.(3)提出了高架軌道交通引起環(huán)境振動(dòng)的數(shù)值模擬模型和計(jì)算方法.總體系可分解成車-橋二維子結(jié)構(gòu)模型和框架-樁基礎(chǔ)-地基三維子結(jié)構(gòu)模型.(4)將上述2個(gè)子結(jié)構(gòu)結(jié)合,對(duì)明珠線典型區(qū)段沿線的環(huán)境振動(dòng)進(jìn)行了實(shí)例計(jì)算.結(jié)果表明,三框架模型算得的沿線地面振動(dòng)振級(jí)衰減規(guī)律比單框架模型更接近于實(shí)測(cè)經(jīng)驗(yàn)公式.該模型可用于預(yù)測(cè)高架軌道交通引起環(huán)境振動(dòng)的振級(jí).(5)利用三框架計(jì)算模型,分析計(jì)算了列車運(yùn)行速度和載客量對(duì)沿線環(huán)境振動(dòng)的影響.計(jì)算表明,車速每增加20km·h-1,振級(jí)約增加4dB.列車載客量對(duì)振級(jí)的影響很小.參考文獻(xiàn):[1] 馬 筠.我國(guó)鐵路環(huán)境振動(dòng)現(xiàn)狀及傳播規(guī)律[J].中國(guó)環(huán)境科學(xué),1987,17(6):70-74.[2] GB10070—88,城市區(qū)域環(huán)境振動(dòng)標(biāo)準(zhǔn)[S].[3] 張 彌,夏 禾.輕軌列車高架軌道系統(tǒng)的動(dòng)力響應(yīng)分析[J].北方交通大學(xué)學(xué)報(bào),1999,18(1):1-8.[4] 于大明,魏 芳.輕軌列車作用下的高架橋梁系統(tǒng)振動(dòng)及其對(duì)周圍環(huán)境的影響[R].北京:北方交通大學(xué)土木建筑系,1996.[5] 張 昕.軌道交通引起環(huán)境振動(dòng)的實(shí)測(cè)和理論研究[D].上海:同濟(jì)大學(xué)結(jié)構(gòu)工程與防災(zāi)研究所,2002.[6] BornitzG. berdieAusbreitungdervonGroszkolbenmashinenErzeugtenBodenschw ingungenindieTiefe[M].Berlin:JSpring,1931.[7] 高木興一,鹽田正純.環(huán)境アセスメント[J].騷音制御,1996,20(4):242-253.[8] 同濟(jì)大學(xué)結(jié)構(gòu)理論研究所.上海市地震小區(qū)劃研究報(bào)告[R].上海:同濟(jì)大學(xué)結(jié)構(gòu)理論研究所,1992.[9] 羅 林.我國(guó)特大橋軌道不平順功率譜的研究[R].北京:中國(guó)鐵道科學(xué)研究院,1998.

 
舉報(bào)收藏 0打賞 0評(píng)論 0
更多>相關(guān)評(píng)論
暫時(shí)沒(méi)有評(píng)論,來(lái)說(shuō)點(diǎn)什么吧
更多>同類論文
推薦圖文
推薦論文
點(diǎn)擊排行
網(wǎng)站首頁(yè)  |  隱私政策  |  版權(quán)隱私  |  使用協(xié)議  |  聯(lián)系方式  |  關(guān)于我們  |  排名推廣  |  廣告服務(wù)  |  積分換禮  |  網(wǎng)站留言  |  幫助中心  |  網(wǎng)站地圖  |  違規(guī)舉報(bào)
津ICP備20006083號(hào)-112010502100290
 
欧美日本高清一区_成人福利网站在线观看11_欧美日本在线视频中文字字幕_国产在线精品一区免费香蕉
国产精品三级美女白浆呻吟| 色偷偷噜噜噜亚洲男人| 亚洲国产精品www| 欧美激情亚洲国产| 国产精品久久久久久五月尺| 欧美成年人视频网站| 欧美精品制服第一页| 中文字幕精品一区日韩| 亚洲免费视频一区| 日韩亚洲一区在线播放| 日韩欧美亚洲日产国产| 欧美精品在线观看| 亚洲欧美日韩国产成人综合一二三区| 亚洲人一区二区| 日本午夜精品电影| 欧美在线视频网| 国产综合在线看| 国产精品一区二区免费| 久久综合久久久久| 久久青草福利网站| 国产精品天天狠天天看| 久久成人精品电影| 亚洲一区二区三区乱码| 日本精品一区二区三区不卡无字幕 | 久久99热精品| 亚洲一区在线免费| 人体内射精一区二区三区| 品久久久久久久久久96高清| 精品午夜一区二区三区| 8050国产精品久久久久久| 国产精品偷伦免费视频观看的| 中文字幕欧美日韩一区二区三区| 日本精品视频网站| 国产欧美日韩伦理| 国产av天堂无码一区二区三区| 国产精品久久久久久影视| 亚洲乱码日产精品bd在线观看| 欧美精品一区二区三区四区五区 | 欧美激情一级精品国产| 欧美一区二区.| 精品无码久久久久久久动漫| 久久人91精品久久久久久不卡| 国产精品高清网站| 日本最新高清不卡中文字幕 | 日韩精品无码一区二区三区| 国产视频一区二区三区在线播放| 69**夜色精品国产69乱| 欧美成年人视频| 日本不卡久久| 国产欧美日韩小视频| 俺也去精品视频在线观看| 亚洲精品日韩av| 国产又黄又猛视频| 久久久久久久久久久久久国产精品 | 日本亚洲导航| 国产精品一国产精品最新章节| 色噜噜狠狠狠综合曰曰曰88av| 美女av一区二区| 欧美日韩第二页| 久久久水蜜桃| 一级日韩一区在线观看| 国精产品一区一区三区视频| www.国产一区| 日本精品一区二区三区高清 久久 日本精品一区二区三区视频 | 久久综合福利| 一区二区精品视频| 国产一区二区不卡视频| 色偷偷91综合久久噜噜| 熟女少妇在线视频播放| 91精品国产91久久久久麻豆 主演| 久久中文字幕一区| 免费看欧美黑人毛片| 国产成人鲁鲁免费视频a| 欧美一区二区三区四区在线观看地址| 国产精品香蕉国产| 欧美激情视频在线免费观看 欧美视频免费一 | 色妞一区二区三区| 性色av香蕉一区二区| 国产伦精品一区二区三区照片| 久久亚洲成人精品| 国产一区二区视频免费在线观看| 久久中文字幕一区| 国产天堂视频在线观看| 免费97视频在线精品国自产拍| 欧美高清性xxxxhd| 国产精品久久久久久久久久| 麻豆久久久9性大片| 久久亚洲精品成人| 国产日韩亚洲欧美在线| 欧美精品在线极品| 成人做爽爽免费视频| 亚洲综合av一区| 久久久女人电视剧免费播放下载| 亚洲一区二区三区视频| 91久久精品久久国产性色也91| 亚洲一区在线直播| 久久久伊人欧美| 日本a级片在线观看| 久久久久久国产免费| 欧美精品久久| 久久av.com| 成人福利网站在线观看11| 中国人体摄影一区二区三区| 91精品视频在线| 视频一区二区视频| 日韩有码在线视频| 黄色动漫在线免费看| 欧美人与性动交| 97精品视频在线| 日本久久久久亚洲中字幕| 久久激情五月丁香伊人| 国产日韩欧美在线视频观看| 亚洲视频精品一区| 久久99国产精品99久久| 精品免费视频123区| 九九热这里只有精品6| 91国产在线播放| 欧美精品一区三区在线观看| 欧美日韩国产成人| 国产xxxxx视频| 国产一区免费在线| 亚洲二区三区四区| 国产精品色午夜在线观看| 国产深夜男女无套内射| 性亚洲最疯狂xxxx高清| 国产成人精品在线播放| 蜜桃av噜噜一区二区三| 亚洲不卡中文字幕无码| 久久99国产精品一区| 国产欧美日韩亚洲| 日韩wuma| 久久久久国产视频| 日韩一区在线视频| 国产精品一区专区欧美日韩| 日本精品久久久久久久久久| 国产精品久久久久久久av电影| 成人免费毛片在线观看| 欧美少妇一区| 亚洲在线免费视频| 国产精品视频精品| 国产精品a久久久久久| 免费观看美女裸体网站| 熟女少妇精品一区二区| 操91在线视频| www日韩欧美| 8050国产精品久久久久久| 国产原创精品| 欧美亚洲国产成人精品| 天堂资源在线亚洲资源| 欧美激情精品久久久久久变态| 日韩在线视频中文字幕| 97久久精品午夜一区二区| 国产专区欧美专区| 欧美性久久久久| 日本久久中文字幕| 亚洲va久久久噜噜噜久久狠狠 | 欧美日韩三区四区| 日产中文字幕在线精品一区| 亚洲最大福利视频| 欧美猛交ⅹxxx乱大交视频| 国产精品私拍pans大尺度在线| 国产精品69精品一区二区三区| 国产日韩av在线播放| 精品www久久久久奶水| 日韩久久不卡| 日本高清一区| 日本一区二区三区免费看| 亚洲国产精品久久久久婷蜜芽| 久久国产精品网站| 久久精品国产亚洲7777| 久久久久久这里只有精品| 久久久免费视频网站| 91精品国产91久久| 91精品久久久久久久| 亚洲 国产 日韩 综合一区| 亚洲熟妇无码一区二区三区导航| 在线视频一区观看| 一区二区三区电影| 亚洲国产精品久久久久婷婷老年 | 亚洲成人午夜在线| 亚洲在线免费看| 亚洲精品高清视频| 中文字幕在线亚洲精品| 亚洲五码在线观看视频| 亚洲人成网站在线播放2019| 一本二本三本亚洲码| 亚洲中文字幕无码专区| 一区二区冒白浆视频| 亚洲熟妇av日韩熟妇在线| 亚洲精品日韩激情在线电影| 亚州成人av在线| 日本中文字幕久久看| 日本免费黄视频| 青青影院一区二区三区四区| 欧美专区一二三| 国模精品一区二区三区| 国产日韩二区| 成人国产精品日本在线| 久久免费看av| 久久视频中文字幕|