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

    渡槽結(jié)構(gòu)橫向動(dòng)力響應(yīng)分析

    2015-01-20 03:29:08杜應(yīng)吉
    關(guān)鍵詞:槽體渡槽極值

    馮 超,杜應(yīng)吉

    (西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院,陜西 楊凌 712100)

    渡槽結(jié)構(gòu)橫向動(dòng)力響應(yīng)分析

    馮 超,杜應(yīng)吉

    (西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院,陜西 楊凌 712100)

    結(jié)合某渡槽安全鑒定工作,建立單跨渡槽結(jié)構(gòu)在設(shè)計(jì)水位和無水兩種工況下的ANSYS三維有限元模型。選用Housner流-固耦合模型模擬了槽內(nèi)水體與槽體側(cè)壁之間的相互作用,并設(shè)定槽墩高度為8.3,10.3,12.3和14.3m,分別進(jìn)行模態(tài)分析和動(dòng)力響應(yīng)分析,以觀察在不同槽墩高度下渡槽結(jié)構(gòu)的動(dòng)力響應(yīng)(應(yīng)力、位移和速度)變化。分析結(jié)果表明:渡槽結(jié)構(gòu)在設(shè)計(jì)水位工況下的自振頻率小于結(jié)構(gòu)在無水工況下的自振頻率;隨著槽墩高度的增加,結(jié)構(gòu)在設(shè)計(jì)水位工況和無水工況下的自振頻率均呈減小的趨勢(shì),而槽墩頂部、槽體跨中及槽體頂部關(guān)鍵點(diǎn)處的動(dòng)力響應(yīng)值有總體增大的趨勢(shì);但渡槽結(jié)構(gòu)不同位置的響應(yīng)值不同,在地震作用下,高墩渡槽的動(dòng)力響應(yīng)值總體大于矮墩渡槽的動(dòng)力響應(yīng)值。

    渡槽;有限元;動(dòng)力響應(yīng);槽墩高度;流-固耦合;地震作用

    馮超,杜應(yīng)吉.渡槽結(jié)構(gòu)橫向動(dòng)力響應(yīng)分析[J].水利水運(yùn)工程學(xué)報(bào),2015(1):68-73.(FENG Chao,DU Ying?Ji.Analysis of transverse dynamic response for aqueduct structure[J].Hydro?Science and Engineering,2015(1):68-73.)

    渡槽是輸送渠道水流跨越河渠、山?jīng)_、谷口等的輸水建筑物,是渠系建筑物中應(yīng)用最廣的交叉建筑物之一[1]。渡槽在正常工作狀態(tài)時(shí),其槽身結(jié)構(gòu)承受與結(jié)構(gòu)自身重量相當(dāng)?shù)乃w,而如此大質(zhì)量的水體在地震作用下引起晃動(dòng)勢(shì)必會(huì)對(duì)槽體自身的動(dòng)力響應(yīng)產(chǎn)生影響,因此如何反映水體的存在以及水體和槽體側(cè)壁之間的相互作用力是研究渡槽結(jié)構(gòu)動(dòng)力響應(yīng)的關(guān)鍵[2]。以往文獻(xiàn)中,對(duì)渡槽結(jié)構(gòu)動(dòng)力響應(yīng)的研究,多關(guān)于槽體過水截面的高寬比、流固耦合模型的選取對(duì)結(jié)構(gòu)動(dòng)力性能的影響,而針對(duì)槽墩高度的變化對(duì)渡槽結(jié)構(gòu)動(dòng)力響應(yīng)的影響則研究較少。但是,實(shí)際工程中,渡槽結(jié)構(gòu)的槽墩高度往往不同。因此,采用ANSYS軟件,并結(jié)合某渡槽安全鑒定工程實(shí)例,依據(jù)Housner流-固耦合模型考慮水體的動(dòng)力作用并建立不同槽墩高度的渡槽結(jié)構(gòu)單跨有限元模型,分別對(duì)結(jié)構(gòu)進(jìn)行模態(tài)分析和輸入地震動(dòng)荷載,探求不同槽墩高度下結(jié)構(gòu)關(guān)鍵點(diǎn)處的動(dòng)力響應(yīng)變化規(guī)律,具有一定的實(shí)際意義。

    1 Housner流-固耦合模型

    水體與槽體之間的流-固耦合作用是影響渡槽結(jié)構(gòu)體系動(dòng)力特性的重要因素。在強(qiáng)震作用下,流體會(huì)產(chǎn)生強(qiáng)烈的晃蕩,同時(shí)會(huì)對(duì)槽體產(chǎn)生很大的橫向作用力,該橫向力的作用不容忽視[3],故需將槽內(nèi)水體作用納入槽體在強(qiáng)震作用下受力分析的考慮中。在流-固耦合模型中選用Housner流-固耦合模型[4],該模型依據(jù)勢(shì)流理論建立的彈簧-質(zhì)量模型,便于工程應(yīng)用同時(shí)避開了求解Laplace方程和無窮級(jí)數(shù)的困難,同Westergaard附加質(zhì)量模型相比,Housner流-固耦合模型更精確[2]。此外,該模型能夠考慮流體的晃動(dòng)作用,已被多位學(xué)者運(yùn)用并取得了一定的成果[5-6]。本文所研究渡槽的槽墩為非柔性槽墩,可以滿足Housner流-固耦合模型的適用條件[7],故本文選用Housner流-固耦合模型來模擬水體同槽體之間的流-固耦合作用。

    依據(jù)Housner流-固耦合模型,流體對(duì)槽體的作用力由脈動(dòng)壓力和對(duì)流壓力兩部分組成,流體在橫向地震激勵(lì)的作用下,其脈動(dòng)壓力用固定于側(cè)壁的等效質(zhì)量塊M0PR來代替,M0PR按下式計(jì)算:,其中:?jiǎn)挝缓穸人w的重量M1=ρhl,h為槽內(nèi)水深;l為過水截面寬度;ρ為流體密度。質(zhì)量塊距離渡槽底板的距離

    流體對(duì)槽體的對(duì)流壓力以水體的奇數(shù)階振動(dòng)來表示,將這些振動(dòng)簡(jiǎn)諧力等效為一系列和側(cè)壁連接的彈簧-質(zhì)量系統(tǒng)并按下式計(jì)算相應(yīng)奇數(shù)階振動(dòng)的水體質(zhì)量、彈簧剛度、距離底板的距離(n=1,3,5…):

    其中:n為流體對(duì)流壓力的奇數(shù)階階數(shù)。

    2 工程概況及有限元模型的建立

    本文以某引水渡槽工程為例,該渡槽全長(zhǎng)84 m,共7跨,單跨長(zhǎng)12 m,槽體各跨之間采用橡膠止水帶連接。渡槽混凝土材料彈性模量為2.2×104N/mm2,混凝土的泊松比為0.167,密度為2.36×103kg/m3,阻尼比為0.05。渡槽過水截面為矩形截面,過水截面面積為11.78 m2(3.1 m×3.8 m),槽體壁厚0.4 m,槽內(nèi)設(shè)計(jì)水位3.3 m,槽墩高度為10.3 m,槽體和槽墩之間由橡膠支座連接。由于渡槽結(jié)構(gòu)各跨的長(zhǎng)度基本相等,各槽墩高度相差不大且結(jié)構(gòu)所處的地基條件基本一致,其振動(dòng)周期基本相等,地震動(dòng)力反應(yīng)也基本一致,因此可選取單跨渡槽進(jìn)行有限元建模研究[8]。

    運(yùn)用ANSYS軟件建立渡槽結(jié)構(gòu)單跨有限元模型,用附加質(zhì)量塊mass21單元來模擬鄰跨的槽體和水體質(zhì)量的影響,針對(duì)設(shè)計(jì)水位和無水工況分別將水體和槽體總質(zhì)量的一半施加在該跨模型的槽墩上[5];對(duì)渡槽結(jié)構(gòu)施加二維、垂直結(jié)構(gòu)軸線方向的地震動(dòng)進(jìn)行分析,故不考慮槽體各跨之間的三維碰撞作用[9],僅對(duì)槽身各跨連接處施加槽身軸線方向的約束。其中渡槽槽身結(jié)構(gòu)選用solid45[10]單元進(jìn)行離散,Housner流-固耦合模型的質(zhì)量塊則用mass21單元進(jìn)行模擬,用beam4單元模擬質(zhì)量塊和槽體側(cè)壁的剛性連接以實(shí)現(xiàn)脈動(dòng)壓力作用,用彈簧單元combin14單元模擬質(zhì)量塊和槽體側(cè)壁連接以實(shí)現(xiàn)對(duì)流壓力作用。質(zhì)量塊mass21單元通過combin14單元和beam4單元同槽體側(cè)壁連接,以實(shí)現(xiàn)Housner流固-耦合模型,其中彈簧-質(zhì)量模型如圖1所示。有限元軟件進(jìn)行模擬分析的過程中,一般選用彈簧單元模擬支座[11],因此槽墩和槽體連接的支座采用彈簧單元combin14[12]模擬,彈簧單元combin14剛度取值采用橋梁支座的取值方法[13]。因?qū)Χ刹劢Y(jié)構(gòu)進(jìn)行橫向地震響應(yīng)分析,故鎖定該單元豎直方向和槽身軸線方向的自由度。單跨渡槽結(jié)構(gòu)有限元模型如圖2所示。

    圖1 彈簧振子-質(zhì)量塊模型Fig.1 A spring?massmodel

    圖2 渡槽結(jié)構(gòu)有限元模型Fig.2 A finite elementmodel for aqueduct structure

    3 計(jì)算結(jié)果及分析

    對(duì)建立的渡槽結(jié)構(gòu)三維有限元模型進(jìn)行模態(tài)分析,比較槽體在設(shè)計(jì)水位工況和無水工況下的自振特性。由于該結(jié)構(gòu)不建在河流中,所以不承受洪水和漂浮物的沖擊作用。在合理的范圍內(nèi)改變槽墩的高度同時(shí)滿足結(jié)構(gòu)的抗滑穩(wěn)定、抗傾覆穩(wěn)定驗(yàn)算要求;由于槽墩截面面積的變化會(huì)改變結(jié)構(gòu)的剛度,從而對(duì)結(jié)構(gòu)的動(dòng)力特性產(chǎn)生影響[14],因此在原模型的基礎(chǔ)上,僅通過改變槽墩高度而不改變槽墩的截面面積相關(guān)參數(shù),同時(shí)結(jié)構(gòu)滿足穩(wěn)定驗(yàn)算的要求,研究比較槽墩高度對(duì)槽體自振特性的影響;對(duì)槽體輸入地震動(dòng)荷載,同時(shí)在實(shí)際工程中,結(jié)合渡槽安全鑒定工程,分析在槽墩高度變化的情況下渡槽結(jié)構(gòu)相應(yīng)的槽墩底部中心、槽墩頂部、槽體跨中及槽體頂端關(guān)鍵點(diǎn)處的動(dòng)力響應(yīng)值的變化情況。

    3.1 自振特性的分析

    渡槽結(jié)構(gòu)的模態(tài)是作為結(jié)構(gòu)動(dòng)力分析的基礎(chǔ)和判斷共振的依據(jù)[14]。文章采用完全法對(duì)渡槽結(jié)構(gòu)進(jìn)行模態(tài)分析,計(jì)算出在不同的槽墩高度(8.3,10.3,12.3和14.3m)下渡槽結(jié)構(gòu)的前15階振動(dòng)頻率(見表1)。由表1可知,在槽墩高度相同的條件下,設(shè)計(jì)水位工況下結(jié)構(gòu)的自振頻率要大于無水工況下結(jié)構(gòu)的自振頻率,說明由于水體的存在增加了結(jié)構(gòu)的總質(zhì)量,減小了結(jié)構(gòu)的自振頻率,符合物理概念;水體的存在對(duì)結(jié)構(gòu)不同階模態(tài)的頻率改變幅度不同;結(jié)構(gòu)在設(shè)計(jì)水位工況和無水工況下的自振頻率均表現(xiàn)出隨著槽墩高度的增加而自振頻率減小的特征,表明增加槽墩高度能夠延長(zhǎng)結(jié)構(gòu)自振周期。通過比較不同槽墩高度下,渡槽結(jié)構(gòu)在設(shè)計(jì)水位工況和無水工況下的前10階振型,槽體的前3階振型(橫向,橫向,縱向)基本不隨槽墩高度的增加而改變;在無水工況下,隨著槽墩高度的增加,槽墩彎曲的振型出現(xiàn)(5,6階振型),表明槽墩高度的增加,使得槽墩更具柔性,故在模態(tài)分析中,出現(xiàn)槽墩彎曲的振型;在設(shè)計(jì)水位工況下,隨著槽墩高度的增加,槽墩彎曲的振型和槽體縱向變形的振型在前10階基本不出現(xiàn),槽體的模態(tài)以橫向振型和豎向振型為主,表明水體的存在和槽墩高度的變化,對(duì)槽體的振型會(huì)產(chǎn)生一定影響。

    表1 不同槽墩高度下的自振頻率Tab.1 Natural frequencies of vibration at different heights of piers Hz

    3.2 地震動(dòng)輸入及結(jié)構(gòu)動(dòng)力響應(yīng)

    工程區(qū)場(chǎng)址為陜西某地區(qū),地震設(shè)防烈度為8度,加速度峰值為0.2g,所在場(chǎng)地相當(dāng)于二類場(chǎng)地,故對(duì)設(shè)計(jì)水位工況和無水工況下的不同槽墩高度的渡槽結(jié)構(gòu)輸入同槽身垂直X向,持時(shí)15 s,加速度峰值為0.2g的ElCentro地震波,其方向與渡槽軸線垂直走向,以便進(jìn)一步了解槽墩高度對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)的影響,并在ANSYS計(jì)算結(jié)果中觀察渡槽結(jié)構(gòu)關(guān)鍵點(diǎn)的動(dòng)力響應(yīng)時(shí)程曲線,地震波如圖3所示。

    圖3 南北方向ElCentro地震動(dòng)Fig.3 ElCentro north and south

    為表示出渡槽結(jié)構(gòu)在設(shè)計(jì)水位工況下其水體對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)的影響以及驗(yàn)證數(shù)值模擬的正確性,作出槽體在無水工況和設(shè)計(jì)水位工況下槽體側(cè)壁中心節(jié)點(diǎn)的位移響應(yīng)時(shí)程曲線圖,如圖4所示(槽墩高度10.3 m)。從圖中可知,在地震動(dòng)的作用下,設(shè)計(jì)水位工況下節(jié)點(diǎn)位移時(shí)程響應(yīng)大于相應(yīng)的無水工況下的節(jié)點(diǎn)位移時(shí)程響應(yīng),表明水體晃蕩對(duì)槽體產(chǎn)生了動(dòng)力作用且這種作用是不可忽視的,模擬結(jié)果同實(shí)際相符合;同時(shí)設(shè)計(jì)水位工況和無水工況下的節(jié)點(diǎn)位移時(shí)程曲線表現(xiàn)出同時(shí)性,也即運(yùn)動(dòng)規(guī)律的相似性。采用擬靜力法,加速度取峰值0.2g,對(duì)渡槽結(jié)構(gòu)在設(shè)計(jì)水位工況和無水工況進(jìn)行抗震計(jì)算。針對(duì)設(shè)計(jì)水位工況,將水體作為附加質(zhì)量考慮,計(jì)算槽墩底部剪力并除以槽墩底部截面面積,將計(jì)算結(jié)果同數(shù)值模擬單元面力極值相對(duì)比,擬靜力法結(jié)果約大于數(shù)值模擬單元面力極值20%~25%。因?yàn)閿M靜力法的計(jì)算結(jié)果偏大于實(shí)際值,所以數(shù)值模擬結(jié)果在合理范圍內(nèi)。

    圖4 側(cè)壁節(jié)點(diǎn)位移時(shí)程Fig.4 Node displacement of aqueductwall along time

    表2~3分別列出渡槽結(jié)構(gòu)在設(shè)計(jì)水位工況和無水工況下,其關(guān)鍵點(diǎn)在整個(gè)動(dòng)力時(shí)程響應(yīng)下的極值,表中的正負(fù)號(hào)只表示方向。表2表明:在不同槽墩高度下,設(shè)計(jì)水位工況下的槽墩底部中心節(jié)點(diǎn)剪力極值均大于無水工況下的墩底剪力極值;槽墩高度的改變對(duì)于兩種工況產(chǎn)生了不同影響;對(duì)于設(shè)計(jì)水位工況下,底部剪力極值隨著槽墩高度的增加而出現(xiàn)相對(duì)減小的趨勢(shì),而在無水工況下則出現(xiàn)隨著槽墩高度增加而剪力極值增加的現(xiàn)象,但增加幅度較小,究其原因,主要在于水體的晃蕩作用對(duì)結(jié)構(gòu)的動(dòng)力特性產(chǎn)生了一定影響。

    表2 墩底中心節(jié)點(diǎn)剪力、槽體底板跨中節(jié)點(diǎn)位移和剪力Tab.2 Shear force of central node of pier bottom and displacement and shear force ofmid?span bottom node

    表2還給出了槽體底板跨中節(jié)點(diǎn)的位移和剪力極值??芍灰坪图袅O值都隨著槽墩高度的增加而增加,如設(shè)計(jì)水位工況下位移極值由墩高8.3 m的0.310 mm變?yōu)槎崭?4.3 m的0.988 mm;在同一槽墩高度下的設(shè)計(jì)水位工況的位移和剪力極值大于無水工況下的位移和剪力極值,如墩高10.3 m時(shí):無水工況下剪力極值為16 566 N,設(shè)計(jì)水位工況下的剪力極值為31 684 N。因此,對(duì)底板可采取相應(yīng)的加固措施來減小水體晃動(dòng)的影響;通過計(jì)算比較可知隨著槽墩高度的增加,節(jié)點(diǎn)位移的增幅減弱。

    表3列出了槽墩頂部節(jié)點(diǎn)和槽體側(cè)壁中部頂端節(jié)點(diǎn)的速度極值。從表中可知,節(jié)點(diǎn)速度極值在兩種不同工況下都隨槽墩高度的增加而變大,且槽墩高度的增加,槽墩頂部節(jié)點(diǎn)速度值的增幅要遠(yuǎn)遠(yuǎn)大于槽體側(cè)壁中部頂端節(jié)點(diǎn)速度值的增幅。對(duì)于不同位置的節(jié)點(diǎn),如槽墩頂部和槽體跨中關(guān)鍵點(diǎn)處的節(jié)點(diǎn)速度極度值相差約5~6倍,這表明槽墩頂部節(jié)點(diǎn)和槽體頂部節(jié)點(diǎn)的速度時(shí)程響應(yīng)不同,槽體頂部節(jié)點(diǎn)的速度時(shí)程響應(yīng)要大于槽墩頂部的速度時(shí)程響應(yīng);也表明對(duì)于渡槽結(jié)構(gòu),上部節(jié)點(diǎn)的速度響應(yīng)極值往往大于下部節(jié)點(diǎn)的速度響應(yīng)極值。

    表3 墩頂和側(cè)壁中部頂端節(jié)點(diǎn)速度Tab.3 Velocity of pier top node and side wall node (cm/s)

    4 結(jié) 語

    (1)從渡槽結(jié)構(gòu)在兩種工況下的模態(tài)分析可知,設(shè)計(jì)水位工況下的槽體自振頻率要小于無水工況下的槽體自振頻率,表明水體的存在相當(dāng)于增加了結(jié)構(gòu)體系質(zhì)量從而延長(zhǎng)結(jié)構(gòu)的自振周期,采用Housner流-固耦合模型能起到合理模擬水體的效果并且物理概念明確;改變槽墩高度后對(duì)結(jié)構(gòu)進(jìn)行模態(tài)分析,可知在無水和有水工況下結(jié)構(gòu)都表現(xiàn)出隨著槽墩高度的增加而自振頻率降低,槽墩高度影響結(jié)構(gòu)的自振周期。在無水工況下,槽墩高度增加,槽墩彎曲的振型開始出現(xiàn)(5,6階振型);在設(shè)計(jì)水位工況下,槽墩高度增加,槽體振型以橫向振型和豎向振型為主。

    (2)比較槽體在無水工況和設(shè)計(jì)水位工況下的結(jié)構(gòu)動(dòng)力響應(yīng)可知,有水工況下的動(dòng)力響應(yīng)大于無水工況下的動(dòng)力響應(yīng),說明由水體晃蕩產(chǎn)生的作用于槽體側(cè)壁的動(dòng)力荷載會(huì)加大渡槽結(jié)構(gòu)的變形量;兩種工況下的時(shí)程響應(yīng)曲線的趨勢(shì)在時(shí)間上表現(xiàn)出較好的一致性。

    (3)比較不同槽墩高度的動(dòng)力響應(yīng)時(shí),墩底中心節(jié)點(diǎn)的剪力極值的變化隨工況不同而出現(xiàn)不同的變化趨勢(shì),無水工況,槽墩高度增加同時(shí)墩底剪力值增加;設(shè)計(jì)水位工況下,槽墩高度增加而墩底剪力極值減小。但是,針對(duì)水體的存在對(duì)槽墩底部剪力極值的影響仍有待于進(jìn)一步的研究,其他位置關(guān)鍵點(diǎn)處的動(dòng)力響應(yīng)極值均隨著槽墩高度的增加而增大,但不同關(guān)鍵點(diǎn)的動(dòng)力響應(yīng)極值的增幅不同。槽體的底板跨中節(jié)點(diǎn)的剪力、位移和槽體側(cè)壁中部頂端節(jié)點(diǎn)的位移增幅較大。渡槽結(jié)構(gòu)的速度極值出現(xiàn)在槽體側(cè)壁頂端節(jié)點(diǎn),應(yīng)對(duì)槽體頂部采用設(shè)置拉桿等方式加固。因此,在實(shí)際工程中,可以通過渡槽結(jié)構(gòu)的槽墩高度,推求結(jié)構(gòu)在承受地震動(dòng)荷載時(shí)關(guān)鍵點(diǎn)處的動(dòng)力響應(yīng)值及其變化規(guī)律。

    [1]季日臣,夏修身,陳堯隆.水體晃蕩作用對(duì)渡槽橫向抗震影響的研究[J].水力發(fā)電學(xué)報(bào),2007,26(6):30?34.(JIRi?chen,XIA Xiu?shen,CHEN Yao?long.Research on influence ofwater shake acting on aqueduct transverse seismic response[J]. Journal of Hydroelectric Engineering,2007,26(6):30?34.(in Chinese))

    [2]張多新,王清云,白新理.大型渡槽結(jié)構(gòu)動(dòng)力學(xué)研究進(jìn)展[J].自然災(zāi)害學(xué)報(bào),2011,20(4):22?30.(ZHANG Duo?xin,WANG Qing?yun,BAIXin?li.Research progress in structural dynamics of large scale aqueduct[J].Journal of Natural Disasters,2011,20(4):22?30.(in Chinese))

    [3]王海萍,倪修全,吳軍中.考慮流固耦合的大型渡槽動(dòng)力分析[J].水電能源科學(xué),2013,31(10):90?93.(WANG Hai?ping,NI Xiu?quan,WU Jun?zhong.Dynamic analysis of large?scale aqueduct considering fluid solid coupling[J].Water Resources and Power,2013,31(10):90?93.(in Chinese))

    [4]居榮初,曾心傳.彈性結(jié)構(gòu)與液體的耦聯(lián)振動(dòng)理論[M].北京:地震出版社,1983:115?122.(JU Rong?chu,ZENG Xin?chuan.Elastic structure with fluid vibration theory[M].Beijing:Earthquake Press,1983:115?122.(in Chinese))

    [5]李彥軍,王賀,吳迪.大型渡槽結(jié)構(gòu)抗震分析研究[J].人民黃河,2010,32(1):100?102.(LIYan?jun,WANG He,WUDi.Large?scale aqueduct structure seismic analysis[J].Yellow River,2010,32(1):100?102.(in Chinese))

    [6]季日臣,蘇小鳳,嚴(yán)娟.水體質(zhì)量對(duì)大型梁式渡槽橫向抗震性能影響研究[J].地震工程工程學(xué)報(bào),2013,35(3):569?574.(JIRi?Chen,SU Xiao?feng,YAN Juan.Transverse seismic performance of beam aqueduct considering watermass influence[J].China Earthquake Engineering Journal,2013,35(3):569?574.(in Chinese))

    [7]劉云賀,胡寶柱,閏建文,等.Housner模型在渡槽抗震計(jì)算中的適用性[J].水利學(xué)報(bào),2002(9):94?99.(LIU Yun?he,HU Bao?zhu,YAN Jian?Wen,etal.Applicability of Housermodel to aseismic characteristics calculation of aqueduct[J].Journal of Hydraulic Research,2002(9):94?99.(in Chinese))

    [8]張伯艷,鄧迎,李德玉.渡槽抗震計(jì)算中幾個(gè)關(guān)鍵問題的簡(jiǎn)化處理方法[J].南水北調(diào)與水利科技,2005,3(2):46?48.(ZHANG Bo?yan,DENG Ying,LIDe?yu.Simplification treatment of some key problems in aqueduct aseismic calculation[J]. South?to?North Water Transfers and Water Science&Technology,2005,3(2):46?48.(in Chinese))

    [9]周浩.大型渡槽結(jié)構(gòu)三維地震碰撞反應(yīng)分析[D].鄭州:鄭州大學(xué),2013.(ZHOU Hao.Three?dimensional seismic pounding response analysis of large?scale aqueduct structures[D].Zhengzhou:Zhengzhou University,2013.(in Chinese))

    [10]王新敏.ANSYS工程結(jié)構(gòu)數(shù)值分析[M].北京:人民交通出版社,2007:378?384.(WANG Xin?min.ANSYS engineering structure numerical simulation[M].Beijing:China Communication Press,2007:378?384.(in Chinese))

    [11]謝旭.橋梁結(jié)構(gòu)地震響應(yīng)分析與抗震設(shè)計(jì)[M].北京:人民交通出版社,2006:20?23.(XIE Xu.Bridge structure seismic response analysis and seismic design[M].Beijing:China Communication Press,2006:20?23.(in Chinese))

    [12]包陳,王呼佳.ANSYS工程分析進(jìn)階實(shí)例[M].北京:中國水利水電出版社,2009:249?253.(BAO Chen,WANG Hu?jia. Examples of ANSYS engineering analysis[M].Beijing:ChinaWaterPower Press,2009:249?253.(in Chinese))

    [13]羅如登,葉梅新,莫朝慶.橋梁支座水平靜力約束方向抗震中的彈簧剛度取值方法對(duì)比研究[J].鐵道科學(xué)與工程學(xué)報(bào),2008,5(2):23?28.(LUO Ru?deng,YEMei?xin,MO Chao?qin.The comparison study on valuingmethod of the stiffness on the direction of horizontal static constraint of support in seismic finite element analysis on bridges[J].Journal of Railway Science and Engineering,2008,5(2):23?28.(in Chinese))

    [14]王博,李杰.大型渡槽結(jié)構(gòu)模態(tài)分析[J].地震工程與工程振動(dòng),2000,20(3):60?66.(WANG Bo,LI Jie.Large scale aqueductmodal analysis[J].Earthquake Engineering and Engineering Vibration,2000,20(3):60?66.(in Chinese))

    Analysis of transverse dynam ic response for aqueduct structure

    FENG Chao,DU Ying?ji
    (College ofWater Resources and Architectural Engineering,Northwest Agricultural and Forestry University,Yangling 712100,China)

    For safety appraisal of an aqueduct,an ANSYS three?dimensional finite elementmodel of a single?span aqueduct structure has been developed with design water level and withoutwater.The fluid?solid interaction Housner model is selected to simulate the interaction between water and the aqueduct wall and different pier heights of the aqueduct are set:8.3 m,10.3 m,12.3 m,14.3 m respectively.The dynamic response of the aqueduct structure under different pier heights(stress,displacement,velocity)based on modal analysis and dynamic analysis is observed.Analysis results show that the natural vibration frequencies of the aqueduct structure under the design water level are less than those without water.With increase of the pier height of the aqueduct,the frequencies of the structure show the trend of decrease;the dynamic response values of the key points of the top of the pier,themid?span and the top of the mid?span have the tendency of increase;however,the dynamic response values of the different position of the aqueduct structure are different,therefore,the dynamic response values of the high?pier aqueduct is generally greater than that of the short?pier aqueduct under the practical conditions of earthquake.

    aqueduct; finite element method; dynamic response; piers height; fluid?solid interaction;earthquake action

    TV672+3;TU311.3

    A

    1009-640X(2015)01-0068-06

    10.16198/j.cnki.1009-640X.2015.01.010

    2014-07-15

    國家“十一五”科技支撐計(jì)劃重點(diǎn)資助項(xiàng)目(2006BAD-11B03)

    馮 超(1989-),男,黑龍江齊齊哈爾人,碩士研究生,主要從事水工材料耐久性和水工建筑物安全與修復(fù)技術(shù)研究。E?mail:avalon.f@163.com通信作者:杜應(yīng)吉(E?mail:dyj@nwsuaf.edu.cn)

    猜你喜歡
    槽體渡槽極值
    隧道式食品油炸機(jī)槽體熱變形研究
    淺談水利工程渡槽基礎(chǔ)承臺(tái)施工
    極值點(diǎn)帶你去“漂移”
    一種用于TPE塑料條的冷卻槽
    極值點(diǎn)偏移攔路,三法可取
    傍山道路泥石流渡槽設(shè)計(jì)
    退鉻槽有限元分析及結(jié)構(gòu)優(yōu)化
    退鉻槽有限元分析及結(jié)構(gòu)優(yōu)化
    一類“極值點(diǎn)偏移”問題的解法與反思
    袁曉峰攝影作品欣賞
    河北水利(2017年6期)2017-09-05 11:04:13
    捣出白浆h1v1| 日韩 亚洲 欧美在线| 日日爽夜夜爽网站| 亚洲色图 男人天堂 中文字幕| 欧美xxⅹ黑人| av有码第一页| 国产精品不卡视频一区二区| 国产在线视频一区二区| 免费观看av网站的网址| 成人毛片a级毛片在线播放| 欧美变态另类bdsm刘玥| 亚洲国产精品999| 国产免费福利视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品一区二区大全| 汤姆久久久久久久影院中文字幕| 亚洲天堂av无毛| 亚洲av男天堂| 久久影院123| 久久国内精品自在自线图片| 美女视频免费永久观看网站| 伊人久久国产一区二区| 日韩欧美一区视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费观看性视频| 日本黄色日本黄色录像| 日韩av不卡免费在线播放| 亚洲伊人久久精品综合| 黑人猛操日本美女一级片| 国产日韩一区二区三区精品不卡| 美女福利国产在线| 女性被躁到高潮视频| 丝瓜视频免费看黄片| 久久精品aⅴ一区二区三区四区 | 亚洲色图综合在线观看| 国产成人免费无遮挡视频| 777久久人妻少妇嫩草av网站| 亚洲美女搞黄在线观看| 深夜精品福利| 你懂的网址亚洲精品在线观看| 在线观看免费日韩欧美大片| 久久久精品国产亚洲av高清涩受| 一级毛片电影观看| 国产在视频线精品| 国产一级毛片在线| av视频免费观看在线观看| 亚洲av电影在线进入| 国产黄色免费在线视频| 一本—道久久a久久精品蜜桃钙片| 啦啦啦中文免费视频观看日本| 丝袜脚勾引网站| 欧美国产精品va在线观看不卡| 日本欧美国产在线视频| 欧美 日韩 精品 国产| 99久久人妻综合| 国产日韩欧美在线精品| 高清视频免费观看一区二区| 寂寞人妻少妇视频99o| 一区二区三区精品91| 日日撸夜夜添| 亚洲国产精品一区三区| 少妇猛男粗大的猛烈进出视频| 中文字幕精品免费在线观看视频| 三级国产精品片| 伊人久久大香线蕉亚洲五| 久久狼人影院| 黄色配什么色好看| 在线看a的网站| 婷婷色麻豆天堂久久| 18禁裸乳无遮挡动漫免费视频| 青青草视频在线视频观看| 黄色视频在线播放观看不卡| 亚洲国产精品国产精品| 久久人人97超碰香蕉20202| 一区二区日韩欧美中文字幕| 久久久久国产网址| 久久精品国产亚洲av高清一级| 五月伊人婷婷丁香| 日韩中文字幕视频在线看片| 久久青草综合色| 2021少妇久久久久久久久久久| 香蕉精品网在线| 成人二区视频| 午夜精品国产一区二区电影| 日韩 亚洲 欧美在线| 婷婷色av中文字幕| 99久久人妻综合| 少妇人妻 视频| 大香蕉久久成人网| 男女下面插进去视频免费观看| 波野结衣二区三区在线| 精品福利永久在线观看| 一本—道久久a久久精品蜜桃钙片| 国产精品熟女久久久久浪| 国产xxxxx性猛交| 亚洲av日韩在线播放| 久久久久国产一级毛片高清牌| 中文乱码字字幕精品一区二区三区| 久久国产精品男人的天堂亚洲| 亚洲精品aⅴ在线观看| 精品午夜福利在线看| 午夜激情av网站| 18+在线观看网站| 欧美日韩成人在线一区二区| 高清不卡的av网站| 国产日韩欧美在线精品| 街头女战士在线观看网站| 最近中文字幕高清免费大全6| 精品国产超薄肉色丝袜足j| 久久精品aⅴ一区二区三区四区 | 超碰成人久久| 极品少妇高潮喷水抽搐| 国产成人精品无人区| 欧美国产精品一级二级三级| 999精品在线视频| 九草在线视频观看| 欧美精品国产亚洲| 色婷婷久久久亚洲欧美| 国产精品国产av在线观看| 日产精品乱码卡一卡2卡三| 国产日韩一区二区三区精品不卡| 国产乱人偷精品视频| 欧美日韩综合久久久久久| 亚洲成国产人片在线观看| 亚洲精品日本国产第一区| 午夜免费男女啪啪视频观看| 观看av在线不卡| 丝袜美足系列| 午夜福利乱码中文字幕| 久久精品国产自在天天线| 国产免费视频播放在线视频| 免费黄频网站在线观看国产| 1024视频免费在线观看| 看免费av毛片| 日本欧美视频一区| 91aial.com中文字幕在线观看| 人人澡人人妻人| 男人舔女人的私密视频| 欧美精品人与动牲交sv欧美| 在线观看国产h片| 女人被躁到高潮嗷嗷叫费观| 亚洲av日韩在线播放| av天堂久久9| 亚洲国产欧美日韩在线播放| 国产免费福利视频在线观看| 一区二区三区四区激情视频| 成年女人在线观看亚洲视频| 国产国语露脸激情在线看| 国产在视频线精品| 国产精品二区激情视频| 日韩三级伦理在线观看| 黄色 视频免费看| 女性生殖器流出的白浆| 高清av免费在线| 一二三四在线观看免费中文在| 国产精品二区激情视频| 777米奇影视久久| 久久精品夜色国产| 人体艺术视频欧美日本| 夫妻午夜视频| 十八禁高潮呻吟视频| 精品福利永久在线观看| 亚洲精品av麻豆狂野| 在线观看人妻少妇| 成人黄色视频免费在线看| 久久精品国产综合久久久| 久久狼人影院| 一区二区av电影网| 亚洲五月色婷婷综合| 亚洲在久久综合| 精品久久久久久电影网| 一二三四在线观看免费中文在| 日韩,欧美,国产一区二区三区| 亚洲综合色惰| 最近2019中文字幕mv第一页| 夫妻性生交免费视频一级片| 亚洲内射少妇av| 免费黄色在线免费观看| 丝袜在线中文字幕| 纯流量卡能插随身wifi吗| a级毛片黄视频| 国产爽快片一区二区三区| 高清视频免费观看一区二区| 黄频高清免费视频| 日韩av在线免费看完整版不卡| 99香蕉大伊视频| 久久午夜综合久久蜜桃| 亚洲欧美日韩另类电影网站| 老汉色∧v一级毛片| 乱人伦中国视频| 日本91视频免费播放| 久热这里只有精品99| 国产有黄有色有爽视频| 国产av码专区亚洲av| 免费观看av网站的网址| 男女边摸边吃奶| 国产精品人妻久久久影院| 一边亲一边摸免费视频| 男女国产视频网站| 中文乱码字字幕精品一区二区三区| 中文字幕精品免费在线观看视频| 97在线人人人人妻| 成人午夜精彩视频在线观看| 亚洲久久久国产精品| 最近中文字幕2019免费版| 久久久精品免费免费高清| 制服人妻中文乱码| 亚洲国产欧美在线一区| 熟女少妇亚洲综合色aaa.| 人人妻人人添人人爽欧美一区卜| 搡老乐熟女国产| 久久久精品94久久精品| 最近最新中文字幕免费大全7| 日本猛色少妇xxxxx猛交久久| 99久久综合免费| 水蜜桃什么品种好| 五月天丁香电影| 美女视频免费永久观看网站| 亚洲av在线观看美女高潮| 岛国毛片在线播放| 午夜福利在线观看免费完整高清在| 亚洲经典国产精华液单| 国产精品无大码| 青青草视频在线视频观看| 在线观看美女被高潮喷水网站| 国产精品欧美亚洲77777| 久久99精品国语久久久| 成年美女黄网站色视频大全免费| 亚洲精品美女久久av网站| 免费播放大片免费观看视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 下体分泌物呈黄色| 最近的中文字幕免费完整| 亚洲图色成人| 国产亚洲av片在线观看秒播厂| 久久久精品区二区三区| 亚洲精品中文字幕在线视频| 波多野结衣av一区二区av| 日韩欧美一区视频在线观看| 欧美av亚洲av综合av国产av | 精品午夜福利在线看| 91aial.com中文字幕在线观看| 夫妻性生交免费视频一级片| 国产男女超爽视频在线观看| 99热国产这里只有精品6| 一级毛片我不卡| 丁香六月天网| 精品福利永久在线观看| 日本爱情动作片www.在线观看| av在线app专区| 国产av码专区亚洲av| 两个人看的免费小视频| 侵犯人妻中文字幕一二三四区| 人妻一区二区av| 国产成人91sexporn| 一二三四在线观看免费中文在| 大香蕉久久网| 18禁裸乳无遮挡动漫免费视频| 欧美日韩精品网址| 熟女av电影| 性少妇av在线| 曰老女人黄片| 1024香蕉在线观看| 亚洲色图 男人天堂 中文字幕| 99久久精品国产国产毛片| 亚洲少妇的诱惑av| 不卡av一区二区三区| 欧美日韩国产mv在线观看视频| 啦啦啦中文免费视频观看日本| 9色porny在线观看| 午夜福利视频精品| 精品卡一卡二卡四卡免费| 亚洲四区av| 尾随美女入室| 99re6热这里在线精品视频| 美女脱内裤让男人舔精品视频| 午夜福利,免费看| 女性被躁到高潮视频| 又粗又硬又长又爽又黄的视频| 高清黄色对白视频在线免费看| 国产成人精品在线电影| 国产无遮挡羞羞视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产伦理片在线播放av一区| 免费在线观看黄色视频的| 亚洲国产色片| 亚洲精华国产精华液的使用体验| 在线观看免费日韩欧美大片| 午夜日韩欧美国产| 国产又爽黄色视频| 街头女战士在线观看网站| 一区二区三区乱码不卡18| 午夜福利一区二区在线看| 日韩熟女老妇一区二区性免费视频| 美女国产视频在线观看| 黄色毛片三级朝国网站| 精品久久久精品久久久| 欧美成人午夜精品| 国产深夜福利视频在线观看| 国产午夜精品一二区理论片| 少妇的逼水好多| av网站免费在线观看视频| 男女边吃奶边做爰视频| 国产一区二区 视频在线| av国产精品久久久久影院| 精品卡一卡二卡四卡免费| 国产片特级美女逼逼视频| 国语对白做爰xxxⅹ性视频网站| 成人影院久久| 色哟哟·www| 成人影院久久| 国产精品女同一区二区软件| 亚洲四区av| 久久久欧美国产精品| 少妇被粗大猛烈的视频| 好男人视频免费观看在线| av视频免费观看在线观看| 成人国产麻豆网| 亚洲欧美中文字幕日韩二区| 尾随美女入室| 亚洲av男天堂| 最近的中文字幕免费完整| 九草在线视频观看| 97在线视频观看| 18禁动态无遮挡网站| 国产又色又爽无遮挡免| 五月伊人婷婷丁香| a 毛片基地| 亚洲成色77777| 亚洲欧洲精品一区二区精品久久久 | 亚洲av男天堂| 男人爽女人下面视频在线观看| 十八禁网站网址无遮挡| 国产精品欧美亚洲77777| 亚洲欧美日韩另类电影网站| 最近中文字幕2019免费版| 99久久中文字幕三级久久日本| 黑人欧美特级aaaaaa片| 99九九在线精品视频| 国产精品免费视频内射| 少妇精品久久久久久久| 国产毛片在线视频| 最近中文字幕2019免费版| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 欧美变态另类bdsm刘玥| 国产一区二区在线观看av| 一级a爱视频在线免费观看| 久久久国产精品麻豆| 欧美国产精品一级二级三级| 久久人人97超碰香蕉20202| 91午夜精品亚洲一区二区三区| 嫩草影院入口| 宅男免费午夜| 少妇的丰满在线观看| 香蕉精品网在线| 人成视频在线观看免费观看| 在线观看免费高清a一片| 精品一区二区三卡| 午夜av观看不卡| 在线观看免费高清a一片| 另类亚洲欧美激情| 高清黄色对白视频在线免费看| 久久人人爽av亚洲精品天堂| 精品99又大又爽又粗少妇毛片| 久久女婷五月综合色啪小说| 久久久精品国产亚洲av高清涩受| 精品少妇内射三级| 另类精品久久| 国产极品天堂在线| 黑人猛操日本美女一级片| 2021少妇久久久久久久久久久| 如何舔出高潮| 高清视频免费观看一区二区| 亚洲一码二码三码区别大吗| 91aial.com中文字幕在线观看| √禁漫天堂资源中文www| 啦啦啦视频在线资源免费观看| 婷婷色av中文字幕| 久久久精品免费免费高清| 中文字幕人妻丝袜制服| 亚洲国产av新网站| 啦啦啦中文免费视频观看日本| 天堂俺去俺来也www色官网| 久久99一区二区三区| 天天躁日日躁夜夜躁夜夜| 18禁国产床啪视频网站| 亚洲精品国产一区二区精华液| 91成人精品电影| 欧美日韩视频精品一区| a 毛片基地| 亚洲中文av在线| 各种免费的搞黄视频| 一区二区三区激情视频| 亚洲国产精品一区二区三区在线| 午夜免费鲁丝| 亚洲精品一区蜜桃| 国产黄频视频在线观看| 国产免费视频播放在线视频| 人妻 亚洲 视频| 欧美日韩成人在线一区二区| 我要看黄色一级片免费的| 成人午夜精彩视频在线观看| 欧美人与性动交α欧美精品济南到 | 亚洲精品日韩在线中文字幕| 亚洲精品国产色婷婷电影| 下体分泌物呈黄色| 天天影视国产精品| 男女下面插进去视频免费观看| 亚洲欧洲日产国产| 精品久久久精品久久久| 日本欧美视频一区| 国产麻豆69| 免费在线观看视频国产中文字幕亚洲 | 日韩制服丝袜自拍偷拍| 我要看黄色一级片免费的| 欧美成人精品欧美一级黄| 欧美变态另类bdsm刘玥| 1024视频免费在线观看| 日韩av免费高清视频| www日本在线高清视频| 久久久久久久亚洲中文字幕| 最近最新中文字幕大全免费视频 | 999久久久国产精品视频| 黄片小视频在线播放| 久久久久视频综合| 超碰97精品在线观看| 国产精品国产三级专区第一集| 五月天丁香电影| 成人亚洲精品一区在线观看| 国产精品一区二区在线不卡| 久久婷婷青草| 日韩一卡2卡3卡4卡2021年| 国产成人aa在线观看| 国产精品一区二区在线观看99| 欧美日韩视频高清一区二区三区二| 成年动漫av网址| 国产麻豆69| 久久人人爽人人片av| 一级毛片黄色毛片免费观看视频| 亚洲国产精品一区三区| 成人午夜精彩视频在线观看| 人成视频在线观看免费观看| 我要看黄色一级片免费的| 国产成人精品在线电影| 亚洲色图综合在线观看| 国产野战对白在线观看| 人成视频在线观看免费观看| 综合色丁香网| 亚洲精品国产av蜜桃| 午夜福利在线免费观看网站| 七月丁香在线播放| 水蜜桃什么品种好| 欧美+日韩+精品| 亚洲欧美成人综合另类久久久| 王馨瑶露胸无遮挡在线观看| 亚洲第一av免费看| 人人妻人人澡人人看| 男人舔女人的私密视频| 国产一区二区在线观看av| 新久久久久国产一级毛片| 国产精品一国产av| 国产精品一二三区在线看| 高清黄色对白视频在线免费看| 男男h啪啪无遮挡| 久久亚洲国产成人精品v| 国产成人免费观看mmmm| 国产在线一区二区三区精| 视频区图区小说| 精品国产超薄肉色丝袜足j| 天天躁夜夜躁狠狠躁躁| 一级片免费观看大全| 亚洲四区av| 叶爱在线成人免费视频播放| 欧美精品国产亚洲| 中文字幕制服av| 日韩中字成人| 99久久精品国产国产毛片| 日韩,欧美,国产一区二区三区| 午夜久久久在线观看| 三上悠亚av全集在线观看| 国产熟女欧美一区二区| 国产亚洲午夜精品一区二区久久| 如何舔出高潮| 男女下面插进去视频免费观看| 欧美日韩视频高清一区二区三区二| 乱人伦中国视频| 午夜日本视频在线| 少妇猛男粗大的猛烈进出视频| 亚洲中文av在线| 毛片一级片免费看久久久久| 90打野战视频偷拍视频| 男男h啪啪无遮挡| 夫妻性生交免费视频一级片| av国产久精品久网站免费入址| 亚洲av电影在线进入| 午夜91福利影院| 你懂的网址亚洲精品在线观看| 人妻人人澡人人爽人人| 国产熟女欧美一区二区| 欧美人与性动交α欧美软件| 亚洲久久久国产精品| 一区二区三区四区激情视频| 日韩欧美一区视频在线观看| 亚洲欧洲国产日韩| 免费黄色在线免费观看| 国产熟女欧美一区二区| 午夜免费鲁丝| 在线观看免费视频网站a站| 久久这里只有精品19| 亚洲精品成人av观看孕妇| 最近中文字幕2019免费版| 黄片播放在线免费| 久久国产精品男人的天堂亚洲| 久久精品久久久久久噜噜老黄| 少妇人妻精品综合一区二区| 国产高清不卡午夜福利| 国产综合精华液| 老汉色av国产亚洲站长工具| av网站在线播放免费| 午夜免费观看性视频| 9色porny在线观看| 性高湖久久久久久久久免费观看| 亚洲精品国产av成人精品| 国产成人精品无人区| 国产精品99久久99久久久不卡 | 电影成人av| 高清欧美精品videossex| 久久精品国产亚洲av高清一级| 老司机亚洲免费影院| 熟女电影av网| 国产成人免费无遮挡视频| 成人免费观看视频高清| 国产精品久久久久久av不卡| 十八禁网站网址无遮挡| 我的亚洲天堂| 少妇被粗大的猛进出69影院| 亚洲欧洲日产国产| 精品亚洲成国产av| 国产极品粉嫩免费观看在线| 老汉色∧v一级毛片| 日产精品乱码卡一卡2卡三| 曰老女人黄片| 日本免费在线观看一区| 精品少妇内射三级| 国产欧美亚洲国产| 欧美日韩视频高清一区二区三区二| 久久久久精品人妻al黑| 99久久综合免费| 亚洲欧美精品自产自拍| 欧美日韩精品网址| 国产精品成人在线| 男人舔女人的私密视频| 一本—道久久a久久精品蜜桃钙片| 美女中出高潮动态图| 国产视频首页在线观看| 老汉色av国产亚洲站长工具| 亚洲国产欧美网| 满18在线观看网站| 亚洲精品成人av观看孕妇| 国产精品 国内视频| 黄色 视频免费看| 观看av在线不卡| 女的被弄到高潮叫床怎么办| 欧美日韩视频精品一区| 在线观看www视频免费| 亚洲人成电影观看| 亚洲五月色婷婷综合| 1024视频免费在线观看| 麻豆乱淫一区二区| 中文天堂在线官网| 亚洲国产色片| 欧美精品一区二区大全| 国产老妇伦熟女老妇高清| 午夜福利在线观看免费完整高清在| 美国免费a级毛片| 制服丝袜香蕉在线| 免费日韩欧美在线观看| 亚洲国产精品一区二区三区在线| 热99国产精品久久久久久7| 国产熟女欧美一区二区| 老司机影院成人| 亚洲精品自拍成人| 成年av动漫网址| 自拍欧美九色日韩亚洲蝌蚪91| kizo精华| 色播在线永久视频| 男女免费视频国产| 久久久精品区二区三区| 亚洲国产av新网站| 91成人精品电影| 大片电影免费在线观看免费| 亚洲男人天堂网一区| 精品卡一卡二卡四卡免费| 亚洲一码二码三码区别大吗| 久久人人爽人人片av| 在线观看免费日韩欧美大片| 香蕉精品网在线| 99热网站在线观看| 青春草视频在线免费观看| 久久久精品94久久精品| 亚洲,欧美,日韩| 午夜免费鲁丝| 人人妻人人澡人人爽人人夜夜| 日韩欧美精品免费久久| 寂寞人妻少妇视频99o| 国产老妇伦熟女老妇高清| a级毛片在线看网站| 一区二区三区激情视频| 最黄视频免费看| 十八禁高潮呻吟视频| 精品午夜福利在线看| 成人午夜精彩视频在线观看| 在线观看美女被高潮喷水网站| 久久精品人人爽人人爽视色|