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

    基于模型試驗(yàn)的高壩泄洪誘發(fā)場(chǎng)地振動(dòng)影響因素研究

    2016-09-18 02:45:36練繼建張文皎李松輝
    振動(dòng)與沖擊 2016年16期
    關(guān)鍵詞:振源消力池孔口

    張 龑, 練繼建, 劉 昉, 張文皎, 李松輝

    (1.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072;2.中國(guó)水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100038)

    ?

    基于模型試驗(yàn)的高壩泄洪誘發(fā)場(chǎng)地振動(dòng)影響因素研究

    張龑1,2, 練繼建1, 劉昉1, 張文皎1, 李松輝2

    (1.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津300072;2.中國(guó)水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100038)

    高壩泄洪引起巨大水流脈動(dòng)荷載產(chǎn)生的振動(dòng),通過(guò)地基沿河谷向上下游傳播,結(jié)果會(huì)導(dǎo)致場(chǎng)地振動(dòng),對(duì)周圍建筑物及居民正常生活造成不利影響。通過(guò)改進(jìn)傳統(tǒng)物理模型,結(jié)合原型觀測(cè)建立關(guān)聯(lián)體系,研究高壩泄洪誘發(fā)場(chǎng)地振動(dòng)的振源特性,以及影響振動(dòng)強(qiáng)度的關(guān)鍵因素。結(jié)果表明:通過(guò)一系列隔振、減振措施建立的完全水彈性模型,能較好的模擬實(shí)際泄洪振動(dòng)情況;通過(guò)線性預(yù)測(cè),能準(zhǔn)確的反映實(shí)際場(chǎng)地振動(dòng)情況。泄流時(shí)孔口處脈動(dòng)荷載是誘發(fā)壩體地基以及場(chǎng)地振動(dòng)的首要振源,導(dǎo)墻和底板荷載次之。泄洪過(guò)程中不同的孔口開(kāi)啟方式,以及不同的上、下游水位是影響場(chǎng)地振動(dòng)強(qiáng)度的主要因素。通過(guò)調(diào)節(jié)泄流方式可以從振動(dòng)源頭上有效的控制場(chǎng)地振動(dòng)強(qiáng)度。地基剛度較小區(qū)域會(huì)對(duì)場(chǎng)地振動(dòng)產(chǎn)生放大效應(yīng)。

    水彈性模型;水流脈動(dòng)荷載;泄流方式;振源;振動(dòng)強(qiáng)度

    高壩大庫(kù)工程是水資源綜合利用和水能資源開(kāi)發(fā)的需要,一直倍受世界各國(guó)重視。由于高壩泄水時(shí)泄洪落差大,水體攜帶能量巨大,泄洪引起的水流脈動(dòng)荷載會(huì)誘發(fā)泄流結(jié)構(gòu)及其他水工建筑物產(chǎn)生振動(dòng),振動(dòng)由大壩上部結(jié)構(gòu)傳遞至大壩基礎(chǔ),進(jìn)而通過(guò)地基傳遞至周邊場(chǎng)地,對(duì)周邊環(huán)境造成影響[1]。一直以來(lái),高壩泄洪引起振動(dòng)問(wèn)題主要局限于泄洪建筑物本身,其引起周邊環(huán)境振動(dòng)的研究鮮有報(bào)道。原因在于高壩大庫(kù)大多距離城市較遠(yuǎn),泄洪引起的振動(dòng)在傳遞過(guò)程中已經(jīng)有很大的衰減;另外,泄洪引起的環(huán)境振動(dòng)振幅和能量相對(duì)較小,其影響形式是一個(gè)長(zhǎng)期累積過(guò)程,一般不會(huì)像地震那樣造成較大而直接的破壞。然而,高壩泄洪引起的水流脈動(dòng)荷載地基向周邊場(chǎng)地傳遞過(guò)程中,遇到特殊的場(chǎng)地條件,會(huì)對(duì)振動(dòng)產(chǎn)生“放大效應(yīng)”,而場(chǎng)地上的房屋或其他建筑物又可能會(huì)對(duì)場(chǎng)地上的振動(dòng)進(jìn)行“二次放大”,經(jīng)“多次放大”的場(chǎng)地振動(dòng)會(huì)造成一定的環(huán)境危害,即對(duì)建筑物的結(jié)構(gòu)安全和人的身體心理產(chǎn)生的不利影響。國(guó)內(nèi)外有關(guān)高壩泄洪引起的周邊環(huán)境振動(dòng),對(duì)建筑物的結(jié)構(gòu)安全和人的身體心理產(chǎn)生影響已有事例。

    俄羅斯的Zhigulevskaya水電站1979年汛期宣泄大洪水曾引起了左岸城市場(chǎng)地和房屋強(qiáng)烈振動(dòng),甚至在距大壩左岸3 km遠(yuǎn)處的居民樓也產(chǎn)生了可感振動(dòng)。長(zhǎng)期的場(chǎng)地振動(dòng)導(dǎo)致大壩附近區(qū)域一些房屋開(kāi)裂,使附近居民出現(xiàn)焦慮,頭暈,惡心等現(xiàn)象[2]。我國(guó)溪洛渡水電站曾在主汛期壩身4個(gè)深孔開(kāi)啟泄洪,壩區(qū)下游右岸混凝土拌合系統(tǒng)的制冰樓(一棟5層鋼結(jié)構(gòu)建筑)發(fā)生明顯振動(dòng)[3]。2012年10月我國(guó)某水電站中孔開(kāi)閘泄洪期間,中孔啟閉機(jī)室聲振、塔帶機(jī)立柱振動(dòng),下游縣城部分門店卷簾門晃動(dòng)、民居及校舍的門窗響動(dòng)、家具顫動(dòng)、吊燈擺動(dòng)等,振動(dòng)現(xiàn)象引起壩區(qū)和城區(qū)民眾的不安。

    因此,對(duì)泄洪誘發(fā)場(chǎng)地振動(dòng)問(wèn)題進(jìn)行研究十分必要。高壩泄洪產(chǎn)生的振動(dòng)由水流引起,與水力學(xué)條件和泄流結(jié)構(gòu)自身的動(dòng)力特性有著密切關(guān)系。由于流固耦合振動(dòng)機(jī)理較為復(fù)雜,通過(guò)完全水彈性模擬,建立流固耦合問(wèn)題的物理模型,進(jìn)行流激振動(dòng)試驗(yàn),是解決高壩泄洪流激振動(dòng)問(wèn)題的有效研究手段[4]。

    本文以中國(guó)某底流消能水電站為研究背景,針對(duì)泄流誘發(fā)場(chǎng)地振動(dòng)現(xiàn)象,提出傳統(tǒng)水彈性模型的改進(jìn)措施,結(jié)合原型觀測(cè)資料,對(duì)模型進(jìn)行修正。同時(shí)建立大壩與場(chǎng)地的關(guān)聯(lián)體系,研究泄流誘發(fā)場(chǎng)地振動(dòng)的振源特性及影響場(chǎng)地振動(dòng)強(qiáng)度主要因素,以期為高壩結(jié)構(gòu)減振研究提供理論基礎(chǔ)。

    1 工程簡(jiǎn)介

    以國(guó)內(nèi)某高壩為研究背景,該壩為混凝土重力壩,壩頂高程384.00 m,最大壩高162.00 m,壩頂長(zhǎng)度909.26 m,正常蓄水位380.00 m,校核洪水位381.86 m,死水位及汛限水位370.00 m。壩身由12個(gè)表孔和10個(gè)中孔組成,采用表中孔間隔布置形式,由中導(dǎo)墻分割成兩個(gè)對(duì)稱的消能區(qū),屬于典型的高壩底流消能的布置形式。

    水電站距下游縣城區(qū)距離最近處僅為0.5 km,最遠(yuǎn)可達(dá)2.5 km。根據(jù)實(shí)際觀測(cè),下游縣城靠近右岸山體一側(cè)場(chǎng)地振感明顯,存在振動(dòng)放大現(xiàn)象。而根據(jù)下游場(chǎng)地的地質(zhì)勘探成果:下游場(chǎng)地表面為砂層覆蓋層,下部基礎(chǔ)為泥巖。下游主城區(qū)靠近右岸山體一側(cè)存在一條古河道,平均寬度200 m左右,古河道上的覆蓋層最大厚度達(dá)80 m。與之相比,主城區(qū)其他區(qū)域的場(chǎng)地覆蓋層厚度明顯較小。不同區(qū)域位置見(jiàn)圖1。

    圖1 水電站壩區(qū)與場(chǎng)地基礎(chǔ)示意圖Fig.1 Schematic diagram of the dam area and ground of the hydropower station

    2 模擬原理

    泄流激振的水彈性試驗(yàn)?zāi)M,是對(duì)“結(jié)構(gòu)-水體-地基-動(dòng)荷載”四位一體的流固耦合振動(dòng)系統(tǒng)的模擬,要求同時(shí)滿足“動(dòng)力荷載”輸入系統(tǒng)相似和結(jié)構(gòu)系統(tǒng)動(dòng)力響應(yīng)相似,即要求同時(shí)滿足水力學(xué)條件和結(jié)構(gòu)動(dòng)力學(xué)條件相似[5]。

    2.1底流消能水力學(xué)條件相似

    水力學(xué)條件相似,即“動(dòng)力荷載”輸入系統(tǒng)相似,實(shí)質(zhì)是按重力相似律設(shè)計(jì)的水力學(xué)模型中脈動(dòng)壓力相似律問(wèn)題。天津大學(xué)通過(guò)系列比尺模型試驗(yàn)和原型觀測(cè)對(duì)比,發(fā)現(xiàn)由整體邊界或水流條件急劇變化引起的水流分離或擴(kuò)散所導(dǎo)致的邊壁壓力脈動(dòng),主要是受低頻大尺度渦漩運(yùn)動(dòng)控制,在雷諾數(shù)足夠大的模型中大尺度渦旋運(yùn)動(dòng)能夠得到較好的模擬,脈動(dòng)壓力可按重力相似定律引伸到原型。該結(jié)論通過(guò)不同比尺水躍模型試驗(yàn)(見(jiàn)圖2~3)結(jié)合二灘水電站原型觀測(cè)資料分析得出(見(jiàn)圖4)。

    圖2 水躍脈動(dòng)壓力強(qiáng)度系數(shù)分布圖Fig.2 Distribution of fluctuating pressure intensity coefficient of hydraulic jump

    圖3 不同比尺歸一化功率譜圖Fig.3 Normalized power spectrums of different scales

    圖4 二灘水電站水墊塘原型與模型數(shù)據(jù)比較(按重力律換算為原型)Fig.4 Comparison between prototype and model data of cushion pool of Ertan hydropower station(Conversion according to the law of gravitation)

    2.2結(jié)構(gòu)動(dòng)力學(xué)條件相似

    結(jié)構(gòu)動(dòng)力學(xué)條件相似指的是結(jié)構(gòu)系統(tǒng)的動(dòng)力響應(yīng)相似,它與結(jié)構(gòu)的頻率、振型、阻尼等因素有關(guān),包括幾何條件相似、物理?xiàng)l件相似、運(yùn)動(dòng)條件相似、邊界條件相似[6]。

    依據(jù)水力學(xué)相似條件和結(jié)構(gòu)動(dòng)力學(xué)相似條件相協(xié)調(diào)的原則水彈性模型模擬時(shí)應(yīng)該滿足:等密度(λγs=1)、等阻尼比(λε=1)、彈性模量小(λE=λL)、等泊松比(λμ=1)等的相似準(zhǔn)則,并合理地選取水域和地基的模擬范圍來(lái)保證結(jié)構(gòu)動(dòng)力響應(yīng)系統(tǒng)的相似[7]。

    3 水彈性模型改進(jìn)

    3.1模型試驗(yàn)精度及影響因素

    場(chǎng)地振動(dòng)及大壩結(jié)構(gòu)振動(dòng)都屬微幅振動(dòng),振動(dòng)信號(hào)采集過(guò)程中易受噪聲和干擾信號(hào)的影響。因此,在開(kāi)展試驗(yàn)前需研究各種干擾因素對(duì)數(shù)據(jù)有效性的影響。首先,對(duì)試驗(yàn)場(chǎng)地背景振動(dòng)情況進(jìn)行24小時(shí)監(jiān)測(cè),各個(gè)時(shí)間段試驗(yàn)場(chǎng)地地面3個(gè)方向的加速度典型的監(jiān)測(cè)結(jié)果見(jiàn)圖5。

    圖5 試驗(yàn)場(chǎng)地地面背景振動(dòng)Fig.5 Background vibration of the test ground

    可以看出,地面的背景振動(dòng)加速度均方根值在凌晨相對(duì)較低,在白天的8∶00-12∶00,14∶00-16∶00及17∶00-19∶00等幾個(gè)時(shí)段振幅較大,主要由于在幾個(gè)交通高峰時(shí)段,場(chǎng)地附近道路交通振動(dòng)引起。原型較小流量(300~3 360 m3/s)時(shí),消力池廊道地面各位置測(cè)點(diǎn)的豎向加速度均方根值范圍為0.01~0.11 gal,因此,水彈模型基礎(chǔ)的微小振動(dòng)信號(hào)容易被試驗(yàn)場(chǎng)地的背景振動(dòng)噪聲淹沒(méi)。

    其次,試驗(yàn)過(guò)程中,水泵運(yùn)行產(chǎn)生振動(dòng)可能對(duì)結(jié)果產(chǎn)生影響,水泵運(yùn)行前后試驗(yàn)場(chǎng)地測(cè)點(diǎn)豎向振動(dòng)情況見(jiàn)圖6,典型測(cè)點(diǎn)頻譜見(jiàn)圖7。

    圖6 水泵運(yùn)行前后試驗(yàn)場(chǎng)地測(cè)點(diǎn)豎向振動(dòng)Fig.6 Test ground vibration before and after the pump operation in vertical direction

    圖7 水泵運(yùn)行前后試驗(yàn)場(chǎng)地測(cè)點(diǎn)豎向振動(dòng)頻譜圖Fig.7 Test ground spectrum before and after the pump operation in vertical direction

    可以看出,水泵開(kāi)啟后產(chǎn)生振動(dòng)通過(guò)地基傳遞至水彈性模型基礎(chǔ),引起模型基礎(chǔ)振動(dòng)增大,8個(gè)基礎(chǔ)測(cè)點(diǎn)的豎向加速度均方根整體提高10%左右。開(kāi)泵前,模型基礎(chǔ)背景振動(dòng)優(yōu)勢(shì)頻帶在0~20 Hz;開(kāi)泵后,基礎(chǔ)測(cè)點(diǎn)譜圖高頻成分比重普遍有所增加,并在24.5 Hz附近出現(xiàn)了一個(gè)峰值,水泵轉(zhuǎn)速為24.5 r/s,因此,增加的振動(dòng)正是由水泵運(yùn)行引起。

    圖8 距模型基礎(chǔ)不同距離人工激勵(lì)加速度衰減Fig.8 Acceleration attenuation from different distance of model based on the artificial excitation

    可以看出,對(duì)于同一種人工激勵(lì)方式,激振力相同,隨著激勵(lì)點(diǎn)距離增大,水彈性模型基礎(chǔ)的振動(dòng)加速度有所減小。當(dāng)距離達(dá)到7 m左右時(shí),隨著距離的進(jìn)一步增大,水彈性模型基礎(chǔ)的加速度衰減較慢。同時(shí),尾水達(dá)到7 m時(shí),水流平順。

    3.2改進(jìn)措施

    通過(guò)對(duì)影響因素的分析,在制作水彈模型時(shí),底部素填土基礎(chǔ)進(jìn)行開(kāi)挖置換加固處理,澆注成C50混凝土基礎(chǔ)(長(zhǎng)7 m×寬7 m×深1.5 m),并在基礎(chǔ)四周預(yù)留一條隔振溝(寬10 cm×深1.5 m),以減小背景振動(dòng)及水泵運(yùn)行產(chǎn)生振動(dòng)對(duì)模型的影響。另外,為了防止水彈性模型上、下游銜接段形成背景振動(dòng)至水彈模型的振動(dòng)傳遞路徑,在上、下游銜接段與水彈模型連接處采用橡膠帶形成軟連接,減小模型周邊振動(dòng)對(duì)水彈模型試驗(yàn)的影響。同時(shí),跌水位置選定在距水彈性模型基礎(chǔ)末端下游7 m處。圖9為模型隔振、減振措施。

    改進(jìn)后的試驗(yàn)?zāi)P透粽駵蟽?nèi)外測(cè)點(diǎn)功率譜圖見(jiàn)圖10。水泵開(kāi)啟前后,隔振系統(tǒng)對(duì)振源的高頻成分皆有一定的隔振效果,隔振溝內(nèi)外地面的低頻背景振動(dòng)情況較為一致。從幅值上看,隔振溝內(nèi)外地面測(cè)點(diǎn)豎向加速度均方根值從0.055 gal減小為0.034 gal。試驗(yàn)進(jìn)行(水泵運(yùn)行)過(guò)程中,隔振溝外地面的高頻成分明顯增大(0~200 Hz),而隔振溝內(nèi)水彈模型地面的高頻成分作用明顯降低。隔振溝內(nèi)外地面測(cè)點(diǎn)豎向加速度均方根值從0.081 gal減小為0.052 gal。因此,該減振隔振系統(tǒng),能有效減小水彈模型外部振動(dòng)源的干擾。

    在此基礎(chǔ)上建立模型比尺1∶80的兩個(gè)消力池、模擬地基(長(zhǎng)500 m×寬400 m×深90 m)的完全水彈性模型,見(jiàn)圖11。振動(dòng)加速度傳感器測(cè)點(diǎn)布置見(jiàn)圖12。

    3.3結(jié)果修正及驗(yàn)證

    某泄洪工況場(chǎng)地豎向振動(dòng)實(shí)測(cè)相對(duì)幅值見(jiàn)圖13。其中,T1、T2和T3為大壩測(cè)點(diǎn),其余為場(chǎng)地測(cè)點(diǎn)。以場(chǎng)地振幅最大測(cè)點(diǎn)T9測(cè)點(diǎn)為例,統(tǒng)計(jì)出13個(gè)典型的原型觀測(cè)工況,對(duì)應(yīng)進(jìn)行模型試驗(yàn)研究。典型工況見(jiàn)表1,模型各結(jié)構(gòu)各測(cè)點(diǎn)與原型T9測(cè)點(diǎn)豎向振動(dòng)振幅對(duì)比見(jiàn)圖14。

    我經(jīng)常在廚房聽(tīng)音頻,與其說(shuō)廚房的活兒干得時(shí)間夠長(zhǎng),足夠聽(tīng)完一期《奇葩說(shuō)》,不如說(shuō)廚房里總有干不完的活兒:洗菜、切菜、炒菜、清理水池、擦洗灶臺(tái)……這些活兒,很難想象不聽(tīng)著點(diǎn)兒啥的話,能堅(jiān)持干完哪怕一件。

    從圖14中可知,水彈性模型各基礎(chǔ)測(cè)點(diǎn)與原型場(chǎng)地T9測(cè)點(diǎn)的豎向振動(dòng)變化趨勢(shì)和規(guī)律較為一致:首先,振動(dòng)強(qiáng)度整體上隨消力池過(guò)流量增大而增大;其次,振動(dòng)強(qiáng)度對(duì)泄洪方式比較敏感,不利的開(kāi)孔方式能明顯增大振動(dòng)強(qiáng)度。水彈性模型測(cè)點(diǎn)振動(dòng)加速度與原型T9測(cè)點(diǎn)振動(dòng)位移的相關(guān)性分析見(jiàn)圖15。

    可以看出,模型消力池基礎(chǔ)內(nèi)、壩踵基礎(chǔ)測(cè)點(diǎn)與原型場(chǎng)地振動(dòng)的相關(guān)系數(shù)分別為0.84和0.86,而水彈模型基礎(chǔ)其余測(cè)點(diǎn)的豎向加速度與原型T9測(cè)點(diǎn)的實(shí)測(cè)豎向位移、豎向加速度的相關(guān)系數(shù)在0.71~0.89內(nèi),相關(guān)性較高。因此,試驗(yàn)結(jié)果可以有效反映實(shí)際泄流狀況。

    圖10 隔振溝內(nèi)外地面測(cè)點(diǎn)豎向振動(dòng)功率譜密度對(duì)比Fig.10 The comparison of power spectrums between inside and outside isolation trench in vertical direction

    圖11 水彈性模型Fig.11 Water elastic model

    圖12 振動(dòng)加速度傳感器測(cè)點(diǎn)布置圖(黑色圓點(diǎn)表示傳感器所在位置)Fig.12 Arrangement of measuring points of vibration acceleration sensor(The black dots indicate the location of sensors)

    圖13 實(shí)測(cè)場(chǎng)地豎向振動(dòng)位移Fig.13 Ground vibration displacements in flow direction

    工況編號(hào)消力池流量/(m3·s-1)泄洪工況13306#、10#中孔局開(kāi)1.3m23601#、5#中孔局開(kāi)1.5m35501#、5#中孔局開(kāi)2.4m46861#、5#中孔局開(kāi)1.4m,2#、4#中孔局開(kāi)1.3m57006#、10#中孔局開(kāi)3.1m69406#、8#、10#中孔局開(kāi)2.4m79601#、5#中孔局開(kāi)4.5m811006#、8#、10#中孔局開(kāi)3.3m929001#、3#、5#中孔局開(kāi)3.0m,6#、8#、10#中孔局開(kāi)3.6m,7#、9#中孔局開(kāi)3.4m1033201#、3#、5#中孔局開(kāi)3.9m,2#、4#中孔局開(kāi)3.4m,6#、8#、10#中孔局開(kāi)2.8m,7#、9#中孔局開(kāi)2.4m1133601#、3#、5#中孔局開(kāi)4.4m,6#、8#、10#中孔局開(kāi)3.6m,7#、9#中孔局開(kāi)3.2m1234301#、3#、5#中孔局開(kāi)3.9m,2#、4#中孔局開(kāi)3.5m,6#、8#、10#中孔局開(kāi)2.8m,7#、9#中孔局開(kāi)2.6m136600左池5中孔局開(kāi)6.0m,右池5中孔局開(kāi)5.8m

    圖14 模型各測(cè)點(diǎn)振動(dòng)與原型T9測(cè)點(diǎn)豎向振動(dòng)對(duì)比Fig.14 Comparison of vertical vibration between the model measuring points and the prototype measuring point T9

    圖15 模型消力池基礎(chǔ)測(cè)點(diǎn)Fig.15 Correlation between measuring points of the base of stilling pool

    研究泄流誘發(fā)場(chǎng)地振動(dòng)相關(guān)問(wèn)題時(shí),建立大壩模型與場(chǎng)地原型觀測(cè)數(shù)據(jù)的聯(lián)系是進(jìn)行各種分析的前提。在模型上進(jìn)行反演2013年和2014年所有原型泄洪工況,根據(jù)試驗(yàn)結(jié)果,對(duì)場(chǎng)地振動(dòng)進(jìn)行線性預(yù)測(cè),T9豎向振動(dòng)加速度預(yù)測(cè)值與實(shí)測(cè)值對(duì)比見(jiàn)圖16。

    圖16 2013年與2014年T9豎向振動(dòng)加速度預(yù)測(cè)值與實(shí)測(cè)值對(duì)比Fig.16 Comparison between the prediction value and the measured value of acceleration in vertical direction in the year of 2013 and 2014

    可以看出,T9測(cè)點(diǎn)的豎向加速度預(yù)測(cè)結(jié)果與原型實(shí)測(cè)結(jié)果吻合較好,不同流態(tài)與泄洪方式對(duì)原型場(chǎng)地振動(dòng)的影響都能通過(guò)水彈性模型得到直接的反映。因此,可通過(guò)以上改進(jìn)措施建立水彈性模型,并根據(jù)各種泄洪工況下水彈性模型的振動(dòng)規(guī)律,來(lái)研究泄洪誘發(fā)下游場(chǎng)地的振動(dòng)問(wèn)題。

    4 泄流誘發(fā)場(chǎng)地振動(dòng)特性

    4.1振源分析

    泄洪誘發(fā)結(jié)構(gòu)及周邊場(chǎng)地振動(dòng)問(wèn)題的研究,首先要正確判斷產(chǎn)生振動(dòng)的機(jī)制,一般包括外部誘發(fā)的振動(dòng)、不確定誘發(fā)的振動(dòng)和運(yùn)動(dòng)誘發(fā)的振動(dòng),其中最重要的是外部誘發(fā)的振動(dòng),它是由水流脈動(dòng)及壓力脈動(dòng)引起的[8]。在泄流結(jié)構(gòu)自身?xiàng)l件一定的情況下,泄流狀態(tài)下不同類型和特征的水流脈動(dòng)荷載對(duì)于結(jié)構(gòu)振動(dòng)的產(chǎn)生、分布和傳遞起到了關(guān)鍵性作用,不同荷載振源對(duì)泄洪誘發(fā)結(jié)構(gòu)及場(chǎng)地振動(dòng)所做的貢獻(xiàn)也不同[9]。該水電站采用表、中孔聯(lián)合的底流泄流消能方式,通過(guò)消力池內(nèi)形成淹沒(méi)水躍產(chǎn)生表面旋滾和水流間強(qiáng)烈的紊動(dòng)剪切來(lái)消除下泄水流的動(dòng)能,同時(shí)在表孔或中孔溢流面和邊墻上產(chǎn)生水流脈動(dòng)荷載,這些荷載都有可能成為場(chǎng)地振動(dòng)振源。在泄流狀態(tài)下可能出現(xiàn)的水流脈動(dòng)荷載有:孔口脈動(dòng)荷載、消力池底板脈動(dòng)荷載、導(dǎo)墻脈動(dòng)荷載、尾坎脈動(dòng)荷載以及跌坎脈動(dòng)荷載[10],荷載位置見(jiàn)圖17。

    圖17 荷載位置圖Fig.17 Loads positions

    試驗(yàn)測(cè)定了不同工況下各個(gè)振源的荷載變化情況。通過(guò)試驗(yàn)發(fā)現(xiàn),不同泄流方式的振源比重不盡相同,結(jié)果見(jiàn)圖18。兩種泄洪方式水流流態(tài)見(jiàn)圖19。

    圖18 不同泄洪工況振源脈壓強(qiáng)度百分比Fig.18 Percentage of energy generated by vibration sources under different discharging conditions

    圖19 不同泄洪工況水流流態(tài)Fig.19 Flow patterns under different discharging conditions

    從柱狀圖可以看出,不同泄流方式作用在溢流壩上的孔口振動(dòng)能量均大于其他位置。表孔泄洪情況下,導(dǎo)墻的振動(dòng)能量?jī)H次于孔口荷載,消力池底板與尾坎振動(dòng)能量相當(dāng)且小于導(dǎo)墻處,跌坎振動(dòng)能量最小。中孔泄洪情況下,消力池底板的振動(dòng)能量增大,僅次于孔口荷載,導(dǎo)墻振動(dòng)能量有所減小,尾坎與跌坎與表孔泄流時(shí)振動(dòng)能量比重相當(dāng)。

    根據(jù)泄水過(guò)程水流流態(tài)可以看出,不同工況泄水過(guò)程中,溢流壩孔口出口段水流條件極為復(fù)雜,紊動(dòng)最為強(qiáng)烈,因此,孔口荷載產(chǎn)生振動(dòng)最為復(fù)雜且能量巨大。中孔泄流時(shí)出流水舌由于出口高程較低,進(jìn)入消力池后形成淹沒(méi)射流,水舌以較大的流速?zèng)_擊到消力池底板,在底板上產(chǎn)生較大的脈動(dòng)壓力,因此其產(chǎn)生脈動(dòng)壓力增大。表孔泄流時(shí)出流水舌出口高程相對(duì)較高,距消力池底板的距離更遠(yuǎn),主流進(jìn)入消力池后主要在表層形成表面旋滾,對(duì)底板的沖擊較小,但導(dǎo)墻的脈動(dòng)壓力增大。兩種泄流方式的下泄水流流經(jīng)尾坎與跌坎時(shí)都已達(dá)到穩(wěn)定狀態(tài),產(chǎn)生能量相當(dāng)。

    4.2場(chǎng)地振動(dòng)影響因素

    4.2.1泄洪方式

    驗(yàn)證模型有效性時(shí)發(fā)現(xiàn),下游場(chǎng)地振動(dòng)強(qiáng)度與泄洪流量及泄洪方式關(guān)系密切,其具體表現(xiàn)為表孔、中孔開(kāi)啟方式,上游水位、下游水位的變化等。因此,通過(guò)控制單一變量原則,研究以上因素對(duì)場(chǎng)地振動(dòng)的影響。試驗(yàn)所得各因素影響下,場(chǎng)地T9測(cè)點(diǎn)振動(dòng)預(yù)測(cè)結(jié)果見(jiàn)圖20~24。

    圖23 不同下游水位T9測(cè)點(diǎn)豎向振動(dòng)加速度預(yù)測(cè)Fig.23 Prediction of vertical vibration acceleration of T9 with different downstream water level

    圖24 不同泄洪方式振動(dòng)情況對(duì)比Fig.24 Vertical vibration acceleration of T9 under different discharge conditions

    從圖20~22可知,不同的孔口開(kāi)啟方式振動(dòng)規(guī)律差異較大。中孔與表孔單獨(dú)運(yùn)行時(shí),振動(dòng)隨開(kāi)度增大而增大,其中,中孔存在不利運(yùn)行區(qū),不利開(kāi)度范圍在5~7 m左右,此時(shí)泄流產(chǎn)生漩滾在孔口處,因此在孔口產(chǎn)生較大脈動(dòng)荷載。小開(kāi)度流量較小,大開(kāi)度時(shí)水流流速與脈動(dòng)能量增大,產(chǎn)生漩滾遠(yuǎn)離孔口,兩者在孔口處產(chǎn)生脈動(dòng)荷載較小。表孔運(yùn)行工況不存在不利運(yùn)行區(qū),但在大開(kāi)度時(shí)振動(dòng)強(qiáng)度增長(zhǎng)幅度有所減小。表中孔聯(lián)合泄洪時(shí),場(chǎng)地振動(dòng)對(duì)中孔的開(kāi)度更為敏感,開(kāi)度為1~4 m時(shí),振幅顯著減小。從圖23可知,下游水位越高,振幅越小。而從圖24可知,上游水位越高,振動(dòng)越強(qiáng),但與流量增長(zhǎng)速度相比,振幅增長(zhǎng)量可以忽略。另外,同流量級(jí)下,表孔泄洪時(shí)場(chǎng)地振動(dòng)強(qiáng)度最大,中孔泄洪時(shí)振動(dòng)次之,且存在明顯的不利運(yùn)行區(qū)間,表中孔聯(lián)合泄洪時(shí)振動(dòng)強(qiáng)度明顯小于同流量其余工況。

    4.2.2場(chǎng)地基礎(chǔ)

    觀測(cè)資料顯示,該水電站下游場(chǎng)地振動(dòng)強(qiáng)度分布有所不同,靠近右岸山體一側(cè)古河道區(qū)域出現(xiàn)振動(dòng)放大現(xiàn)象,左側(cè)非古河道區(qū)域振動(dòng)強(qiáng)度不大,泄洪前后場(chǎng)地古河道與非古河道實(shí)測(cè)工況下典型測(cè)點(diǎn)豎向振動(dòng)功率譜見(jiàn)圖25和圖26。

    圖25 場(chǎng)地古河道區(qū)域T9測(cè)點(diǎn)振動(dòng)功率譜密度譜Fig.25 Vibration power spectrum density of T9 in the area of the ancient river channel

    圖26 場(chǎng)地非古河道區(qū)域T12測(cè)點(diǎn)振動(dòng)功率譜密度譜Fig.26 Vibration power spectrum density of T12 in the area of the non-ancient river channel

    將場(chǎng)地典型頻譜分布與場(chǎng)地地質(zhì)條件結(jié)合不難發(fā)現(xiàn),在大壩右岸下游場(chǎng)地范圍,場(chǎng)地的背景振動(dòng)幅值較小,頻譜特性有明顯的規(guī)律:由于古河道上厚覆蓋層的影響,地基剛度相對(duì)較小,場(chǎng)地較“軟”,自振頻率較低,在低頻振動(dòng)激勵(lì)下,一定程度上激發(fā)了地基的自振,背景振動(dòng)響應(yīng)頻譜反映出了地基的自振特性。在古河道區(qū)域測(cè)點(diǎn)所在位置覆蓋層越厚,背景振動(dòng)響應(yīng)頻譜中的地基自振特性越明顯;其他區(qū)域覆蓋層較淺,地基剛度相對(duì)較大,場(chǎng)地較“硬”,自振頻率較高,表現(xiàn)為低頻振動(dòng)激勵(lì)下的隨機(jī)受迫振動(dòng),其背景振動(dòng)響應(yīng)頻譜反映的是地脈動(dòng)振源的頻率特性。

    5 結(jié) 論

    文章利用水彈性模型試驗(yàn)結(jié)合原型觀測(cè)資料,研究高壩泄流引起周圍場(chǎng)地振動(dòng)振源以及影響振動(dòng)強(qiáng)度因素。研究發(fā)現(xiàn):

    首先,傳統(tǒng)的物理模型受環(huán)境振動(dòng)、試驗(yàn)設(shè)備以及其他因素影響,試驗(yàn)結(jié)果難以反映場(chǎng)地振動(dòng)情況。因此,建立模型前,通過(guò)對(duì)模型底部素填土基礎(chǔ)進(jìn)行開(kāi)挖置換加固處理,并在基礎(chǔ)四周預(yù)留一條隔振溝,并在上、下游銜接段與水彈模型連接處采用橡膠帶形成軟連接可以有效減小水彈模型外部振動(dòng)源的干擾。同時(shí),模型試驗(yàn)反演實(shí)際泄洪工況結(jié)果與實(shí)測(cè)場(chǎng)地振動(dòng)變化規(guī)律一致,相關(guān)性較高。其預(yù)測(cè)結(jié)果真實(shí)反映實(shí)際場(chǎng)地振動(dòng)情況。

    其次,大壩在泄流狀態(tài)下可能出現(xiàn)孔口脈動(dòng)荷載、消力池底板脈動(dòng)荷載、導(dǎo)墻脈動(dòng)荷載、尾坎脈動(dòng)荷載以及跌坎脈動(dòng)荷載,這些荷載成為影響場(chǎng)地振動(dòng)的主要振源。不同泄流方式下,作用在孔口的振動(dòng)能量遠(yuǎn)大于其他位置;中孔泄洪情況下,除了孔口振動(dòng)能量為主要振源外,消力池底板的振動(dòng)能量較大;表孔泄洪情況下,導(dǎo)墻的振動(dòng)能量增大到第二,底板的振動(dòng)能量小于尾坎,減小到第四。

    最后,泄水過(guò)程中,孔口開(kāi)啟方式、上下游水位、消力池內(nèi)流水流流態(tài)等因素對(duì)場(chǎng)地振動(dòng)強(qiáng)度有不同程度的影響,而表孔與中孔聯(lián)合泄洪并避開(kāi)中孔不利運(yùn)行區(qū)是減小場(chǎng)地振動(dòng)的有利運(yùn)行工況;不同的場(chǎng)地土體條件對(duì)振動(dòng)的傳播有較大影響,地基剛度較小區(qū)域易激發(fā)土體自振,從而產(chǎn)生振動(dòng)放大效應(yīng)。

    [1] 崔廣濤,馬斌. 泄流誘發(fā)水工結(jié)構(gòu)振動(dòng)問(wèn)題的綜合集成研討[J]. 水利水電科技進(jìn)展,2010,30(1):1-10.

    CUI Guangtao, MA Bin. Comprehensive compound method for flow-induced vibration of hydraulic structures[J]. Advances in Science and Technology of Water Resources, 2010, 30(1): 1-10.

    [2] SHUMAKOVA E M, KOTLYAKOV A V, SHUMAKOV G V. The effect of vibrations in the Zhigulevskii hydropower structure on soils in the nearby territories of Tolyatti city [J]. Water Resources, 2010,37(3): 306-310.

    [3] 練繼建. 金沙江溪洛渡水電站深孔泄流運(yùn)行相關(guān)建筑物振動(dòng)監(jiān)測(cè)分析報(bào)告[D].天津:天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,2013.

    [4] 崔廣濤,彭新民,苑希民,等. 高拱壩泄洪振動(dòng)水彈性模型[J]. 水力學(xué)報(bào),1996(4):1-9.

    CUI Guangtao, PENG Xinmin, YUAN Ximin, et al. Hydrolastic model study on high arch dam vibration induced by discharge flow[J].Journal of Hydraulic Engineering, 1996(4):1-9.

    [5] 練繼建,王繼敏,辜晉德. 水躍區(qū)水流脈動(dòng)壓力頻譜相似律研究[J].科學(xué)通報(bào),2007,52(15):1832-1839.

    LIAN Jijian, WANG Jimin, GU Jinde. Analysis of the similarity law of fluctuating pressures in the hydraulic jump region[J]. Chinese Science Bulletin,2007,52(15):1832-1839.

    [6] 劉昉.水流脈動(dòng)比尺特性及其相似律研究[D].天津:天津大學(xué),2007.

    [7] 練繼建,楊敏.高壩泄流工程[M].北京:中國(guó)水利水電出版社, 2009.

    [8] 吳一紅,李世琴,謝省宗. 拱壩-庫(kù)水-地基耦合系統(tǒng)壩身泄洪動(dòng)力分析[J]. 水利學(xué)報(bào),1996(11):6-13.

    WU Yihong, LI Shiqing, XIE Shengzong. Dynamic analysis for vibration of dam-water-foundation coupling system during discharge[J]. Journal of Hydraulic Engineering, 1996(11): 6-13.

    [9] 張建偉,崔廣濤,馬斌,等. 基于泄流響應(yīng)的高拱壩振源時(shí)域識(shí)別[J]. 天津大學(xué)學(xué)報(bào),2008,41(9):1124-1129.

    ZHANG Jianwei, CUI Guangtao, MA Bin, et al. Time domain identification of excitation source for high arch dam based on discharge flow vibration response[J]. Journal of Tianjin University, 2008, 41(9):1124-1129.

    [10] 練繼建,崔廣濤,黃錦林.導(dǎo)墻結(jié)構(gòu)的流激振動(dòng)研究[J].水利學(xué)報(bào),1998(11):34-38.

    LIAN Jijian, CUI Guangtao, HUANG Jinlin. Study on flow induced vibration of spillway guide wall[J]. Journal of Hydraulic Engineering, 1998(11):34-38.

    Influencing factors of the ground vibration induced by flood discharge of high dams based on model experiments

    ZHANG Yan1,2, LIAN Jijian1, LIU Fang1, ZHANG Wenjiao1, LI Songhui2

    (1.State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072, China;2.State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin,China Institate of Water Resoures and Hydropower Research, Beijing 100038,China)

    While flood is discharged in a high dam, the vibration generated by huge flow fluctuating load can be transmitted to upstream and downstream along the river through the ground, which results in ground vibration and is harmful to the neighboring residents and buildings. In this paper, the traditional physical model was improved, and then a correlation system was established by combining the model with prototype observations. Further, the characteristics of vibration sources as well as factors influencing vibration intensities were investigated. The results show that the established hydroelastic model with a series of vibration isolation and vibration damping can simulate the flood discharging conditions and reflect the ground vibration accurately through a linear prediction. The vibration energy of orifice is the primary vibration source of ground vibration, and the secondary vibration source is the load of guide wall and stilling pool base slab. The orifice opening modes, upstream and downstream water levels seriously influence the vibration intensity of ground. So it suggests to control the vibration intensity from the sources through regulating the discharge modes. The vibration will be amplified when the foundation stiffness is small.

    hydroelastic model; flow fluctuating load; discharge modes; vibration source; vibration intensity

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51379140;51209158;51379177);中國(guó)水科院科研專項(xiàng)(1513)

    2015-04-28修改稿收到日期:2015-08-18

    張龑 男,博士,1988年11月生

    練繼建 男,教授,博士生導(dǎo)師,1965年8月生

    TV61

    A

    10.13465/j.cnki.jvs.2016.16.006

    猜你喜歡
    振源消力池孔口
    交錯(cuò)起爆下爆炸應(yīng)力波的碰撞機(jī)制與破巖效果
    爆炸與沖擊(2024年6期)2024-10-31 00:00:00
    Effects of Prandtl number in two-dimensional turbulent convection?
    一種筒類零件孔口去毛刺工具
    逆作法孔口邊梁內(nèi)力計(jì)算
    考慮振源相互作用的船舶甲板減振方法
    船海工程(2019年1期)2019-03-04 01:58:04
    U型鋼板樁在干河子壩二級(jí)消力池修復(fù)工程中的應(yīng)用
    某攔河大壩下游消力池結(jié)構(gòu)對(duì)消能防沖的影響
    一種小型化低噪聲微波本振源設(shè)計(jì)
    電子制作(2017年9期)2017-04-17 03:00:36
    沖沙閘滲透穩(wěn)定計(jì)算及消力池防沖設(shè)計(jì)
    底流消力池內(nèi)懸柵消能工數(shù)值模擬研究
    噜噜噜噜噜久久久久久91| 精品无人区乱码1区二区| 桃红色精品国产亚洲av| 亚洲欧美清纯卡通| 国产亚洲av嫩草精品影院| 天天躁日日操中文字幕| 亚洲七黄色美女视频| 亚洲av熟女| 狠狠狠狠99中文字幕| 制服丝袜大香蕉在线| 久久久久久大精品| 亚洲国产欧洲综合997久久,| 久久精品综合一区二区三区| 国产成年人精品一区二区| av在线天堂中文字幕| 一区二区三区四区激情视频 | 搞女人的毛片| 又紧又爽又黄一区二区| 成年女人看的毛片在线观看| 国语自产精品视频在线第100页| 在线观看66精品国产| 变态另类成人亚洲欧美熟女| 一级黄色大片毛片| 日韩中文字幕欧美一区二区| 国产免费一级a男人的天堂| av天堂在线播放| 日韩中字成人| 婷婷六月久久综合丁香| 久久精品国产亚洲av香蕉五月| 亚洲av.av天堂| 久久精品国产清高在天天线| 欧美精品国产亚洲| 男女边吃奶边做爰视频| 美女xxoo啪啪120秒动态图| 亚洲精品国产成人久久av| 乱人视频在线观看| 最近中文字幕高清免费大全6 | 亚洲美女视频黄频| 日韩欧美三级三区| 国产精品三级大全| 国产黄色小视频在线观看| a在线观看视频网站| 国产淫片久久久久久久久| 亚洲欧美激情综合另类| 精品日产1卡2卡| 亚洲第一电影网av| 欧美潮喷喷水| 色5月婷婷丁香| 白带黄色成豆腐渣| 身体一侧抽搐| 免费看a级黄色片| 草草在线视频免费看| 88av欧美| 极品教师在线免费播放| 嫩草影院入口| 国产精品嫩草影院av在线观看 | 窝窝影院91人妻| 日韩欧美在线乱码| 网址你懂的国产日韩在线| ponron亚洲| 男女做爰动态图高潮gif福利片| 亚洲精品粉嫩美女一区| 乱人视频在线观看| 在线观看午夜福利视频| 乱系列少妇在线播放| 欧美中文日本在线观看视频| 精品久久久久久久久久久久久| 欧美潮喷喷水| 婷婷亚洲欧美| 国产老妇女一区| 悠悠久久av| 亚洲精品日韩av片在线观看| 国产精品久久久久久精品电影| 99久久中文字幕三级久久日本| 亚洲精品国产成人久久av| 亚洲av熟女| 人妻少妇偷人精品九色| 男人狂女人下面高潮的视频| 国产综合懂色| 国产精品99久久久久久久久| 又紧又爽又黄一区二区| 亚洲专区中文字幕在线| 97人妻精品一区二区三区麻豆| 精品久久久久久久久亚洲 | 精品久久久久久久久久久久久| 美女被艹到高潮喷水动态| 久久久久国产精品人妻aⅴ院| 我的老师免费观看完整版| 成人三级黄色视频| 亚洲精品色激情综合| 国产亚洲91精品色在线| 校园人妻丝袜中文字幕| 1024手机看黄色片| 婷婷丁香在线五月| 国产av麻豆久久久久久久| 国产精品98久久久久久宅男小说| 欧美+亚洲+日韩+国产| 欧美日韩精品成人综合77777| 国产一区二区在线av高清观看| 久久久精品欧美日韩精品| 日本 欧美在线| 两人在一起打扑克的视频| 直男gayav资源| 中文字幕av在线有码专区| 亚洲精品一卡2卡三卡4卡5卡| 国产 一区精品| 国产亚洲91精品色在线| 日本爱情动作片www.在线观看 | 免费观看的影片在线观看| 欧美丝袜亚洲另类 | 亚洲美女视频黄频| 色综合站精品国产| 欧美xxxx黑人xx丫x性爽| 韩国av在线不卡| 真实男女啪啪啪动态图| 天天一区二区日本电影三级| 日本黄大片高清| 自拍偷自拍亚洲精品老妇| 十八禁网站免费在线| 亚洲av一区综合| 国产美女午夜福利| 人妻久久中文字幕网| 啦啦啦韩国在线观看视频| 亚洲欧美日韩东京热| 亚洲电影在线观看av| 免费av不卡在线播放| 日本与韩国留学比较| 国产精品三级大全| 久久天躁狠狠躁夜夜2o2o| 韩国av在线不卡| 国产探花在线观看一区二区| 在线观看66精品国产| www日本黄色视频网| 最近视频中文字幕2019在线8| 亚洲欧美日韩高清专用| 亚洲av不卡在线观看| 色av中文字幕| 看黄色毛片网站| 白带黄色成豆腐渣| 免费看a级黄色片| 九色成人免费人妻av| 18禁裸乳无遮挡免费网站照片| 一级av片app| 免费看光身美女| 真人做人爱边吃奶动态| 日韩一区二区视频免费看| 欧美激情国产日韩精品一区| 嫩草影视91久久| 国产精品综合久久久久久久免费| 最新中文字幕久久久久| 久久久精品欧美日韩精品| 午夜爱爱视频在线播放| 久久久久精品国产欧美久久久| 亚洲无线观看免费| 国内少妇人妻偷人精品xxx网站| 亚洲精品日韩av片在线观看| 噜噜噜噜噜久久久久久91| 一本一本综合久久| 日本三级黄在线观看| 欧美日韩国产亚洲二区| av在线蜜桃| 香蕉av资源在线| а√天堂www在线а√下载| 色视频www国产| 国产精品不卡视频一区二区| 欧美一区二区精品小视频在线| 国产蜜桃级精品一区二区三区| 欧美日韩乱码在线| 欧美色视频一区免费| 啦啦啦观看免费观看视频高清| av福利片在线观看| 男人舔女人下体高潮全视频| 亚洲狠狠婷婷综合久久图片| 亚洲av免费高清在线观看| 亚洲第一电影网av| 免费看光身美女| 91av网一区二区| 极品教师在线免费播放| 亚洲欧美日韩高清专用| 亚洲第一电影网av| 免费看光身美女| 人人妻,人人澡人人爽秒播| 久久久久久大精品| 1024手机看黄色片| 免费黄网站久久成人精品| 中文资源天堂在线| 悠悠久久av| 高清日韩中文字幕在线| 五月玫瑰六月丁香| 中亚洲国语对白在线视频| 国产单亲对白刺激| 韩国av一区二区三区四区| 12—13女人毛片做爰片一| 国产私拍福利视频在线观看| 99久久无色码亚洲精品果冻| 性欧美人与动物交配| 波多野结衣巨乳人妻| 免费无遮挡裸体视频| 麻豆成人av在线观看| 悠悠久久av| 婷婷精品国产亚洲av在线| 九色成人免费人妻av| 亚洲va在线va天堂va国产| 国产成人a区在线观看| 中文亚洲av片在线观看爽| 男女做爰动态图高潮gif福利片| 欧美黑人巨大hd| 一级黄片播放器| 中文字幕av在线有码专区| 国内少妇人妻偷人精品xxx网站| 亚洲精品国产成人久久av| 人妻少妇偷人精品九色| 啦啦啦韩国在线观看视频| 91久久精品电影网| 我的女老师完整版在线观看| 日韩欧美精品免费久久| 国内久久婷婷六月综合欲色啪| 欧美激情久久久久久爽电影| 亚洲自拍偷在线| 一夜夜www| 老司机福利观看| 国产免费男女视频| 午夜影院日韩av| 偷拍熟女少妇极品色| 国产日本99.免费观看| 色5月婷婷丁香| 亚洲美女黄片视频| 九九热线精品视视频播放| 成人性生交大片免费视频hd| 日日摸夜夜添夜夜添小说| 午夜福利高清视频| 亚洲精品色激情综合| 久久久色成人| 国产精品一区二区免费欧美| 国产一级毛片七仙女欲春2| 久久久成人免费电影| 美女xxoo啪啪120秒动态图| 国产精华一区二区三区| 天堂网av新在线| 欧美三级亚洲精品| 五月玫瑰六月丁香| 午夜免费激情av| 天堂影院成人在线观看| 国产乱人伦免费视频| 亚洲美女视频黄频| 97热精品久久久久久| 亚洲精品乱码久久久v下载方式| 欧美精品啪啪一区二区三区| 国产一区二区三区av在线 | 日韩欧美在线乱码| 日本免费一区二区三区高清不卡| 国产精品三级大全| 久99久视频精品免费| 日韩在线高清观看一区二区三区 | 九色成人免费人妻av| 麻豆成人午夜福利视频| 免费人成在线观看视频色| 日韩欧美 国产精品| 丰满乱子伦码专区| 麻豆精品久久久久久蜜桃| 熟妇人妻久久中文字幕3abv| 色在线成人网| 黄色视频,在线免费观看| 亚洲人与动物交配视频| 欧美最黄视频在线播放免费| 12—13女人毛片做爰片一| or卡值多少钱| 91在线精品国自产拍蜜月| 男人狂女人下面高潮的视频| 人妻制服诱惑在线中文字幕| 国产高清激情床上av| 日本一二三区视频观看| av在线老鸭窝| 精品一区二区三区人妻视频| 五月玫瑰六月丁香| 欧美色视频一区免费| 真人做人爱边吃奶动态| 成人国产一区最新在线观看| 亚洲av美国av| 国产人妻一区二区三区在| 人妻丰满熟妇av一区二区三区| 精品一区二区三区人妻视频| 国产高清视频在线观看网站| 日韩亚洲欧美综合| 亚洲欧美日韩高清在线视频| 午夜福利在线观看免费完整高清在 | 久久精品人妻少妇| 亚洲国产欧洲综合997久久,| 亚洲av免费高清在线观看| 乱码一卡2卡4卡精品| 国产精品亚洲一级av第二区| 国产精品免费一区二区三区在线| 成人午夜高清在线视频| 色av中文字幕| 三级男女做爰猛烈吃奶摸视频| 成人国产麻豆网| 欧美成人一区二区免费高清观看| 久久久久精品国产欧美久久久| 九九在线视频观看精品| www.色视频.com| 欧美色视频一区免费| 中文字幕av在线有码专区| 久久精品国产自在天天线| 性色avwww在线观看| 免费观看精品视频网站| 国产女主播在线喷水免费视频网站 | 如何舔出高潮| 国产精品国产高清国产av| 久久99热这里只有精品18| 久久久久久大精品| 美女黄网站色视频| av在线蜜桃| 尾随美女入室| 亚洲最大成人av| 午夜影院日韩av| h日本视频在线播放| 成人国产一区最新在线观看| 国国产精品蜜臀av免费| 韩国av在线不卡| 国产精品乱码一区二三区的特点| 欧美bdsm另类| 亚洲欧美清纯卡通| 99精品在免费线老司机午夜| 大型黄色视频在线免费观看| 国内毛片毛片毛片毛片毛片| 欧美国产日韩亚洲一区| 99视频精品全部免费 在线| 国产免费av片在线观看野外av| 国产亚洲av嫩草精品影院| 欧美最黄视频在线播放免费| 精品人妻1区二区| 精品免费久久久久久久清纯| 真实男女啪啪啪动态图| av在线蜜桃| 中文字幕熟女人妻在线| 欧美3d第一页| 国产熟女欧美一区二区| 琪琪午夜伦伦电影理论片6080| 最新在线观看一区二区三区| 少妇高潮的动态图| 99热只有精品国产| 精品人妻熟女av久视频| 日韩精品中文字幕看吧| 亚洲精品色激情综合| 欧美丝袜亚洲另类 | 18禁在线播放成人免费| 精品久久久久久久久久久久久| 日韩一本色道免费dvd| 久久久久国内视频| 国产激情偷乱视频一区二区| 神马国产精品三级电影在线观看| av视频在线观看入口| 天堂网av新在线| 亚洲欧美日韩无卡精品| 狂野欧美激情性xxxx在线观看| 一级av片app| 熟女电影av网| 精品久久国产蜜桃| 干丝袜人妻中文字幕| 中文亚洲av片在线观看爽| 国国产精品蜜臀av免费| 成人二区视频| 全区人妻精品视频| 露出奶头的视频| 亚洲一区高清亚洲精品| 看免费成人av毛片| 国产91精品成人一区二区三区| 别揉我奶头~嗯~啊~动态视频| 久久久久免费精品人妻一区二区| 欧美bdsm另类| 午夜福利视频1000在线观看| 国产av麻豆久久久久久久| www日本黄色视频网| 一个人观看的视频www高清免费观看| av在线老鸭窝| 亚洲精品在线观看二区| 亚洲欧美日韩卡通动漫| 真实男女啪啪啪动态图| 亚洲精品亚洲一区二区| 精品免费久久久久久久清纯| 午夜久久久久精精品| 久久精品人妻少妇| 啪啪无遮挡十八禁网站| 亚洲三级黄色毛片| 国产午夜精品久久久久久一区二区三区 | 亚洲av中文字字幕乱码综合| 婷婷色综合大香蕉| 亚洲四区av| 成人高潮视频无遮挡免费网站| 国产精品免费一区二区三区在线| 亚洲欧美日韩高清专用| 国产精品国产高清国产av| 亚洲av美国av| 熟女人妻精品中文字幕| 国产大屁股一区二区在线视频| 日韩欧美国产一区二区入口| 国内毛片毛片毛片毛片毛片| av黄色大香蕉| 男人舔奶头视频| 91久久精品国产一区二区成人| 极品教师在线视频| 亚洲色图av天堂| 久久精品影院6| 国产精品久久久久久久久免| 三级男女做爰猛烈吃奶摸视频| 丰满乱子伦码专区| 直男gayav资源| 国产三级在线视频| 亚洲精华国产精华精| 人人妻人人看人人澡| 99热网站在线观看| 又爽又黄无遮挡网站| 又紧又爽又黄一区二区| 噜噜噜噜噜久久久久久91| 欧美高清性xxxxhd video| 神马国产精品三级电影在线观看| 99riav亚洲国产免费| 内射极品少妇av片p| АⅤ资源中文在线天堂| 直男gayav资源| 久久国产乱子免费精品| 给我免费播放毛片高清在线观看| 免费无遮挡裸体视频| 精品午夜福利在线看| 亚洲久久久久久中文字幕| 国产免费一级a男人的天堂| 日日干狠狠操夜夜爽| 99热精品在线国产| 99久久精品热视频| 最新在线观看一区二区三区| 精品一区二区三区人妻视频| 久久天躁狠狠躁夜夜2o2o| 亚洲国产日韩欧美精品在线观看| 美女黄网站色视频| 特大巨黑吊av在线直播| 亚洲专区中文字幕在线| 免费搜索国产男女视频| 成人一区二区视频在线观看| 国产一区二区在线观看日韩| 三级毛片av免费| 悠悠久久av| 在线观看午夜福利视频| 免费看光身美女| 丰满乱子伦码专区| 99热这里只有是精品50| 少妇被粗大猛烈的视频| 美女 人体艺术 gogo| 日本-黄色视频高清免费观看| 日韩av在线大香蕉| 国产精品人妻久久久久久| 91麻豆精品激情在线观看国产| 国产亚洲欧美98| 又粗又爽又猛毛片免费看| 国产亚洲精品久久久com| 黄色日韩在线| 可以在线观看毛片的网站| 高清日韩中文字幕在线| 免费看av在线观看网站| 中文字幕av成人在线电影| 春色校园在线视频观看| 国内精品一区二区在线观看| 久久久久久国产a免费观看| av在线老鸭窝| 男女做爰动态图高潮gif福利片| 国产亚洲欧美98| 中文字幕av成人在线电影| 国产淫片久久久久久久久| 老熟妇仑乱视频hdxx| 男女边吃奶边做爰视频| 老女人水多毛片| 成人三级黄色视频| 精品一区二区三区视频在线| 小说图片视频综合网站| 亚洲欧美日韩卡通动漫| 在线国产一区二区在线| 天堂√8在线中文| 天堂影院成人在线观看| 午夜a级毛片| 国产精品久久久久久久久免| 日本与韩国留学比较| 亚洲成人久久爱视频| 全区人妻精品视频| 男女视频在线观看网站免费| 欧美人与善性xxx| 国产精品爽爽va在线观看网站| 欧美潮喷喷水| 亚洲18禁久久av| 国产成人a区在线观看| 国产精品免费一区二区三区在线| 别揉我奶头 嗯啊视频| 搡女人真爽免费视频火全软件 | 性插视频无遮挡在线免费观看| 欧美绝顶高潮抽搐喷水| a级毛片免费高清观看在线播放| 99精品在免费线老司机午夜| 久久这里只有精品中国| 欧美极品一区二区三区四区| 最好的美女福利视频网| 国内精品久久久久久久电影| 日本免费一区二区三区高清不卡| 禁无遮挡网站| 欧美三级亚洲精品| 午夜老司机福利剧场| h日本视频在线播放| 日本爱情动作片www.在线观看 | 一个人看的www免费观看视频| 91久久精品国产一区二区成人| 国产探花极品一区二区| 美女xxoo啪啪120秒动态图| 18+在线观看网站| 琪琪午夜伦伦电影理论片6080| 人妻丰满熟妇av一区二区三区| 国产精品国产三级国产av玫瑰| 国内久久婷婷六月综合欲色啪| 色噜噜av男人的天堂激情| 国产v大片淫在线免费观看| 亚洲国产精品sss在线观看| 久久久国产成人精品二区| 村上凉子中文字幕在线| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 国产精品美女特级片免费视频播放器| 亚洲成人久久性| 午夜爱爱视频在线播放| 久久久久久九九精品二区国产| 一本一本综合久久| av在线蜜桃| 精品人妻熟女av久视频| 一区二区三区高清视频在线| 春色校园在线视频观看| 亚洲成人免费电影在线观看| 天天躁日日操中文字幕| 变态另类成人亚洲欧美熟女| 国模一区二区三区四区视频| 在线观看免费视频日本深夜| av在线天堂中文字幕| 91狼人影院| 日韩欧美免费精品| 看免费成人av毛片| 亚洲av一区综合| 特级一级黄色大片| 一夜夜www| 免费无遮挡裸体视频| 嫩草影院新地址| 一级av片app| 国产精品无大码| 国产亚洲欧美98| 日韩欧美 国产精品| 久久久国产成人免费| av在线蜜桃| 18禁黄网站禁片午夜丰满| 97碰自拍视频| 亚洲人成伊人成综合网2020| 亚洲精品一区av在线观看| 国产69精品久久久久777片| 成人三级黄色视频| 成熟少妇高潮喷水视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲图色成人| 国产精品免费一区二区三区在线| 亚洲人成网站高清观看| 国产av不卡久久| 欧美最黄视频在线播放免费| 麻豆国产97在线/欧美| 午夜视频国产福利| 国产精品自产拍在线观看55亚洲| 波多野结衣巨乳人妻| 噜噜噜噜噜久久久久久91| 伦精品一区二区三区| 亚洲国产精品成人综合色| 少妇熟女aⅴ在线视频| 午夜福利在线观看吧| 亚洲一区二区三区色噜噜| 国产淫片久久久久久久久| 男人舔奶头视频| 欧美中文日本在线观看视频| 黄片wwwwww| 午夜久久久久精精品| 欧美精品国产亚洲| 啦啦啦韩国在线观看视频| 亚洲精品国产成人久久av| av在线亚洲专区| 亚洲欧美日韩卡通动漫| 伦理电影大哥的女人| 18禁裸乳无遮挡免费网站照片| 日本黄大片高清| 非洲黑人性xxxx精品又粗又长| 日韩欧美精品v在线| 欧美国产日韩亚洲一区| 国产精品一区二区三区四区久久| a级一级毛片免费在线观看| 极品教师在线免费播放| 久久久久国产精品人妻aⅴ院| 成人无遮挡网站| 一区二区三区高清视频在线| 亚洲欧美激情综合另类| 美女免费视频网站| av黄色大香蕉| 精品久久久久久久久av| 国产伦在线观看视频一区| 免费搜索国产男女视频| 精品一区二区三区视频在线观看免费| 麻豆成人午夜福利视频| 国产三级中文精品| 国产v大片淫在线免费观看| 99久久中文字幕三级久久日本| 日日摸夜夜添夜夜添小说| 免费大片18禁| 精品福利观看| 国内精品宾馆在线| 波多野结衣高清无吗| 91麻豆av在线| 丰满的人妻完整版| 成人一区二区视频在线观看| 观看美女的网站|