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

    爆炸地震動(dòng)下儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)試驗(yàn)研究

    2022-11-21 03:38:56張浩天趙雪川宋春明吳紅曉鄭際鏡岳松林程怡豪
    振動(dòng)與沖擊 2022年21期
    關(guān)鍵詞:動(dòng)水儲(chǔ)液儲(chǔ)水

    張浩天, 趙雪川, 宋春明, 吳紅曉, 鄭際鏡, 岳松林, 程怡豪

    (1. 陸軍工程大學(xué) 爆炸沖擊防災(zāi)減災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 南京 210007; 2. 96911部隊(duì), 北京 100010)

    世界格局與戰(zhàn)爭(zhēng)形態(tài)不斷演變,各國(guó)爭(zhēng)相運(yùn)用高新技術(shù)研發(fā)和改進(jìn)武器裝備。在持續(xù)緊張的安全形勢(shì)之下,國(guó)防工程建設(shè)與防護(hù)技術(shù)研究面臨著嚴(yán)峻的挑戰(zhàn)。作為工程防護(hù)關(guān)注的特種結(jié)構(gòu)之一,儲(chǔ)液結(jié)構(gòu)(liquid storage structures ,LSS)在給排水系統(tǒng)和石油化工等方面有著重要的應(yīng)用,關(guān)乎水源、油料等戰(zhàn)略資源的儲(chǔ)備安全[1-3]。在面臨突發(fā)武器打擊和偶然爆炸襲擊,特別是核武器爆炸時(shí),強(qiáng)烈的沖擊地震動(dòng)將導(dǎo)致儲(chǔ)液結(jié)構(gòu)產(chǎn)生動(dòng)力響應(yīng),嚴(yán)重時(shí)將造成結(jié)構(gòu)的損傷破壞,乃至引發(fā)巨大的次生災(zāi)害,對(duì)人民群眾的生命和財(cái)產(chǎn)安全造成極大的威脅[4-6]。因此,有必要開展爆炸地震動(dòng)下儲(chǔ)液結(jié)構(gòu)的動(dòng)力響應(yīng)研究,對(duì)于補(bǔ)充和完善國(guó)防工程理論體系,提高工程防護(hù)能力具有重要意義。

    當(dāng)前,國(guó)內(nèi)外學(xué)者圍繞儲(chǔ)液結(jié)構(gòu)開展了大量研究工作。其中,儲(chǔ)液結(jié)構(gòu)動(dòng)力特性研究是分析結(jié)構(gòu)動(dòng)力響應(yīng)的基礎(chǔ),魏發(fā)遠(yuǎn)等[7-9]將儲(chǔ)液容器簡(jiǎn)化為梁式結(jié)構(gòu),分別基于懸臂梁彎曲、剪切、彎剪變形理論建立了儲(chǔ)液結(jié)構(gòu)的振動(dòng)方程,并指出工程應(yīng)用時(shí)可根據(jù)結(jié)構(gòu)高寬比選用不同變形理論進(jìn)行近似計(jì)算。Kim等[10]對(duì)儲(chǔ)液結(jié)構(gòu)單一壁板進(jìn)行分析,建立了固支-自由-對(duì)邊簡(jiǎn)支板和固支-自由-對(duì)邊固支板兩種計(jì)算模型,發(fā)現(xiàn)隨著壁板長(zhǎng)高比的增大,結(jié)構(gòu)基頻逐漸趨近于懸臂梁模型計(jì)算結(jié)果。Hashemi等[11]提出了儲(chǔ)液結(jié)構(gòu)四面壁板的振動(dòng)分析方法,在水平和豎直方向分別利用連續(xù)梁和懸臂梁變形理論,組合得到了儲(chǔ)液結(jié)構(gòu)折板振型。楊鳴等[12]基于聲固耦合算法分析儲(chǔ)液容器固有特性,發(fā)現(xiàn)液體的附加質(zhì)量效應(yīng)使得結(jié)構(gòu)的固有頻率下降,增大了結(jié)構(gòu)在低頻荷載下發(fā)生強(qiáng)烈振動(dòng)的可能性。杜永峰等通過(guò)數(shù)值算例得到類似的結(jié)論,指出液動(dòng)壓力使結(jié)構(gòu)自振頻率明顯下降,同時(shí)發(fā)現(xiàn)在儲(chǔ)液條件下結(jié)構(gòu)振型變化不大,在分析其振動(dòng)響應(yīng)時(shí)可使用無(wú)液結(jié)構(gòu)振型近似替代。程選生等[13]針對(duì)Winkler彈性地基和彈性構(gòu)造底板的情況,推導(dǎo)了儲(chǔ)液結(jié)構(gòu)液-固耦合振動(dòng)方程,討論了無(wú)量綱參數(shù)對(duì)自振頻率的影響。

    關(guān)于儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng),現(xiàn)有研究多以自然地震動(dòng)作為輸入載荷,進(jìn)行理論推導(dǎo)、數(shù)值分析與振動(dòng)臺(tái)模型試驗(yàn)。地震激勵(lì)下需要考慮液體與結(jié)構(gòu)的相互作用問(wèn)題,以達(dá)到科學(xué)安全的設(shè)計(jì)要求。Housner[14]提出了剛性儲(chǔ)液結(jié)構(gòu)計(jì)算簡(jiǎn)化模型,將液動(dòng)壓力分為兩部分:一是隨結(jié)構(gòu)作同步運(yùn)動(dòng)的液體產(chǎn)生的脈沖壓力;二是結(jié)構(gòu)內(nèi)液體晃動(dòng)產(chǎn)生的對(duì)流壓力,并將液體等效為與結(jié)構(gòu)相連的彈簧-質(zhì)量系統(tǒng),在儲(chǔ)液結(jié)構(gòu)設(shè)計(jì)中被廣泛應(yīng)用。Veletsos[15]將儲(chǔ)液結(jié)構(gòu)考慮為單自由度體系,并假定結(jié)構(gòu)按照給定形式發(fā)生撓曲變形,提出了彈性儲(chǔ)液結(jié)構(gòu)簡(jiǎn)化計(jì)算方法,發(fā)現(xiàn)考慮結(jié)構(gòu)柔性時(shí)液動(dòng)壓力有所增加。Chen等[16]基于廣義單自由度理論并考慮懸臂梁變形邊界條件進(jìn)行儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)案例分析,發(fā)現(xiàn)單自由度方法用于儲(chǔ)液結(jié)構(gòu)設(shè)計(jì)具有較高的精度,同時(shí)應(yīng)考慮高階模態(tài)的影響。Kim等、Hashemi等提出了考慮流固耦合作用和壁板彈性的柔性儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)計(jì)算模型,在計(jì)算脈沖壓力時(shí)疊加了壁板變形引起的脈沖壓力。特別地,針對(duì)帶有彈性壁板的鋼筋混凝土儲(chǔ)液結(jié)構(gòu),程選生等[17]推導(dǎo)了液動(dòng)壓力計(jì)算公式。程選生等[18]采用ADINA軟件建立儲(chǔ)液結(jié)構(gòu)有限元模型,考慮了液體表面重力波的影響,探討了地震烈度和液位高度對(duì)液晃波高和結(jié)構(gòu)壁板變形的影響。張如林等[19]建立流固耦合運(yùn)動(dòng)方程并利用ANSYS軟件建模,分析了地震波頻譜特性對(duì)儲(chǔ)液結(jié)構(gòu)地震響應(yīng)的影響,發(fā)現(xiàn)當(dāng)?shù)卣鸩ㄗ吭街芷诮咏鞴恬盥?lián)振動(dòng)周期時(shí),結(jié)構(gòu)響應(yīng)顯著增強(qiáng)。Radnic等[20]利用振動(dòng)臺(tái)進(jìn)行地震激勵(lì)下儲(chǔ)液結(jié)構(gòu)的動(dòng)力響應(yīng)試驗(yàn)研究,討論了激勵(lì)特性、壁板剛度與儲(chǔ)液水位等參數(shù)對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)的影響。程選生等[21-22]采用人工黏彈性邊界模擬地基效應(yīng)進(jìn)行有限元分析,發(fā)現(xiàn)考慮土-結(jié)構(gòu)相互作用后,內(nèi)部水體液動(dòng)壓力和儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)均有較大程度的減小。

    綜上所述,儲(chǔ)液結(jié)構(gòu)在生命線工程中具有重要應(yīng)用,其安全性備受關(guān)注,地震激勵(lì)下儲(chǔ)液結(jié)構(gòu)的動(dòng)力響應(yīng)研究已經(jīng)取得了一些成果。然而,以往研究側(cè)重于自然地震動(dòng)下儲(chǔ)液結(jié)構(gòu)響應(yīng),在方法上則偏重于理論計(jì)算和模擬仿真。相比之下,爆炸地震動(dòng)加速度信號(hào)具有峰值大、持時(shí)短等特點(diǎn),該震動(dòng)環(huán)境下儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)研究較少,特別是在試驗(yàn)方面,由于操作難度大、危險(xiǎn)系數(shù)高、成本高等,成功案例鮮有,缺少能夠指導(dǎo)工程實(shí)踐的規(guī)律性發(fā)現(xiàn),也不足以有效驗(yàn)證和補(bǔ)充相關(guān)理論與數(shù)值研究結(jié)果。因此,本文開展爆炸地震動(dòng)下儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)試驗(yàn)研究,對(duì)不同強(qiáng)度地震動(dòng)下,處于無(wú)水、淺水、深水液位狀態(tài)的儲(chǔ)液結(jié)構(gòu)進(jìn)行振動(dòng)臺(tái)試驗(yàn),同時(shí)結(jié)合理論方法與相關(guān)設(shè)計(jì)規(guī)范,分析結(jié)構(gòu)振動(dòng)加速度、結(jié)構(gòu)變形、動(dòng)水壓力等動(dòng)力響應(yīng)的變化規(guī)律,以期為國(guó)防工程中儲(chǔ)液結(jié)構(gòu)的研究與設(shè)計(jì)提供有益參考。

    1 試驗(yàn)概況

    1.1 試驗(yàn)設(shè)備

    結(jié)構(gòu)試驗(yàn)于爆炸沖擊震動(dòng)模擬平臺(tái)進(jìn)行,如圖1所示。該平臺(tái)由陸軍工程大學(xué)爆炸沖擊防災(zāi)減災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室研制,用以模擬核武器爆炸沖擊震動(dòng)環(huán)境[23]。平臺(tái)試驗(yàn)時(shí),用掛鉤吊住擺錘,卷?yè)P(yáng)機(jī)通電提升擺錘至預(yù)定高度后釋放,自由下落的擺錘撞擊振動(dòng)臺(tái),提供單次近半正弦加速度脈沖,通過(guò)調(diào)整擺錘下落高度與緩沖墊塊厚度,可改變輸入加速度波的峰值與脈寬。其中,下落高度定義為擺錘自釋放至水平撞擊振動(dòng)臺(tái),其錘頭重心下降高度,示意圖如圖2所示。需要說(shuō)明的是,爆炸震動(dòng)與爆炸參數(shù)、介質(zhì)條件、結(jié)構(gòu)特性等因素有關(guān),該平臺(tái)以沖擊荷載對(duì)爆炸震動(dòng)環(huán)境進(jìn)行模擬,究其原理:一是反應(yīng)譜等效[24-25],將結(jié)構(gòu)系統(tǒng)簡(jiǎn)化為若干單自由度系統(tǒng),爆炸震動(dòng)對(duì)結(jié)構(gòu)的作用效果與對(duì)各單自由度系統(tǒng)作用效果相同時(shí),可對(duì)爆炸震動(dòng)加速度信號(hào)進(jìn)行反應(yīng)譜分析,并以半正弦加速度信號(hào)反應(yīng)譜包絡(luò),則半正弦加速度脈沖峰值和脈沖持時(shí)即為相應(yīng)爆炸震動(dòng)加速度等效值;二是波形等效,爆炸地震動(dòng)的波形難以準(zhǔn)確預(yù)測(cè),通常采用某一較符合實(shí)際情況的假定波形進(jìn)行參數(shù)分析。設(shè)計(jì)中可近似將自由場(chǎng)地運(yùn)動(dòng)參數(shù)取為結(jié)構(gòu)運(yùn)動(dòng)參數(shù),而半正弦類脈沖及其主導(dǎo)的加速度信號(hào)是爆炸試驗(yàn)典型波形之一[26-27]。

    圖1 爆炸沖擊震動(dòng)模擬平臺(tái)Fig.1 Explosion shock and vibration simulation platform

    圖2 擺錘自由下落Fig.2 Free falling of the pendulum bob

    1.2 模型結(jié)構(gòu)及傳感器布置

    試驗(yàn)?zāi)P蜑榈孛媸骄匦螣o(wú)頂蓋儲(chǔ)液結(jié)構(gòu),其尺寸特征參考地下工程典型儲(chǔ)液結(jié)構(gòu),如圖3所示。模型材料為201不銹鋼,長(zhǎng)1 200 mm,寬900 mm,高750 mm,壁厚10 mm,存儲(chǔ)液體為常溫常壓下的自然常用水。儲(chǔ)液結(jié)構(gòu)底板與震動(dòng)平臺(tái)對(duì)應(yīng)布置φ30 mm孔洞,通過(guò)螺栓連接,結(jié)構(gòu)長(zhǎng)邊平行于加載方向。動(dòng)力試驗(yàn)過(guò)程中,布置一系列傳感器來(lái)記錄各種參數(shù),如圖4所示。其中,加速度傳感器(記為A)與動(dòng)水壓力傳感器(記為P)沿加載方向布置,分別測(cè)量加載后結(jié)構(gòu)的振動(dòng)響應(yīng)、結(jié)構(gòu)受到的動(dòng)水壓力,應(yīng)變片(記為S)沿結(jié)構(gòu)短邊方向布置,反映結(jié)構(gòu)的變形情況,以上傳感器均為單向式。

    圖3 矩形儲(chǔ)液結(jié)構(gòu)Fig.3 Rectangular LSS

    圖4 傳感器布置 (mm)Fig.4 Sensors layout (mm)

    1.3 試驗(yàn)工況

    分別進(jìn)行爆炸地震動(dòng)下無(wú)水、300 mm儲(chǔ)水、500 mm儲(chǔ)水結(jié)構(gòu)動(dòng)力響應(yīng)試驗(yàn),墊塊厚度50 mm,通過(guò)調(diào)整擺錘下落高度H1~H5,提供5種幅值的輸入加速度。其中,H1~H5分別為1 140 mm,1 340 mm,1 540 mm,1 740 mm, 1 940 mm。500 mm儲(chǔ)水、擺錘高度H3時(shí)得到的典型輸入加速度時(shí)程曲線,如圖5所示。該工況下獲得幅值15.8g、持續(xù)時(shí)間15.5 ms的單次加速度脈沖。可以看出,結(jié)構(gòu)底部A1處與振動(dòng)臺(tái)臺(tái)面A0處采集的加速度時(shí)程曲線一致性較好,說(shuō)明螺栓連接緊實(shí)可靠,保證了加載過(guò)程中振動(dòng)臺(tái)與結(jié)構(gòu)之間無(wú)相對(duì)滑移,試驗(yàn)的有效性得到驗(yàn)證。

    圖5 典型輸入加速度時(shí)程曲線Fig.5 Typical time history curve of input acceleration

    對(duì)于每種工況進(jìn)行至少兩次平行試驗(yàn),在試驗(yàn)結(jié)果差值小于兩者平均值5%的前提下取平均值作為可信數(shù)據(jù)進(jìn)行后續(xù)分析,否則補(bǔ)充第三次及以上試驗(yàn)直至滿足可信條件。不同工況下得到的地震動(dòng)參數(shù)如表1所示,地面峰值加速度(peak ground acceleration,PGA)范圍為10.0~26.4g,持續(xù)時(shí)間(duration time of ground acceleration, TGA)平均值為14.8 ms。經(jīng)查閱文獻(xiàn)資料[28-29],了解到加速度型地震動(dòng)強(qiáng)度參數(shù)可以較好地表征地震動(dòng)對(duì)結(jié)構(gòu)的潛在破壞作用,同時(shí)考慮到GB 50011—2010《建筑抗震設(shè)計(jì)規(guī)范》[30]將峰值加速度作為地震動(dòng)強(qiáng)度指標(biāo),本文亦選取峰值加速度這一參數(shù)描述地震動(dòng)強(qiáng)度特性,以便于分析爆炸地震動(dòng)對(duì)儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)的影響。此外,由于篇幅有限,而不同地震動(dòng)輸入條件下動(dòng)力響應(yīng)時(shí)程曲線變化規(guī)律基本一致,僅展示500 mm儲(chǔ)水、15.8g地震動(dòng)強(qiáng)度(擺錘高度H3)時(shí)得到的典型時(shí)程數(shù)據(jù),并分析儲(chǔ)液狀態(tài)、地震動(dòng)強(qiáng)度對(duì)響應(yīng)峰值等參數(shù)的影響。

    表1 不同工況下地震動(dòng)參數(shù)Tab.1 Seismic load parameters under different working conditions

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

    2.1 結(jié)構(gòu)動(dòng)力特性分析

    結(jié)構(gòu)的動(dòng)力特性包括固有頻率、阻尼比、振型等,通常認(rèn)為固有頻率對(duì)擋墻等短周期結(jié)構(gòu)的地震響應(yīng)有很大影響[31]。本文振動(dòng)臺(tái)模型試驗(yàn)中,通過(guò)傳遞函數(shù)分析法確定結(jié)構(gòu)的動(dòng)力特性參數(shù)[32],傳遞函數(shù)Ha(ω,zj)計(jì)算為

    (1)

    式中:GXX(ω)為振動(dòng)臺(tái)臺(tái)面A0處加速度時(shí)程信號(hào)的自功率密度函數(shù);GXY(ω,zj)為結(jié)構(gòu)測(cè)點(diǎn)Aj(j=2,3,4)處加速度時(shí)程與振動(dòng)臺(tái)臺(tái)面A0處加速度時(shí)程信號(hào)的互功率密度函數(shù)。描繪復(fù)函數(shù)Ha(ω,zj)的幅值-頻率曲線,其首個(gè)峰值對(duì)應(yīng)的頻率fn即為結(jié)構(gòu)固有頻率。

    不同儲(chǔ)液狀態(tài)下得到的加速度傳遞函數(shù)幅值-頻率曲線,如圖6所示。由圖6可知,結(jié)構(gòu)壁板不同高度測(cè)點(diǎn)A2,A3,A4得到的結(jié)構(gòu)固有頻率基本一致,取3處測(cè)點(diǎn)計(jì)算結(jié)果的平均值為結(jié)構(gòu)固有頻率。相比于無(wú)水結(jié)構(gòu)(fn=54.15 Hz),300 mm儲(chǔ)水結(jié)構(gòu)(fn=48.84 Hz)和500 mm儲(chǔ)水結(jié)構(gòu)(fn=38.56 Hz)自振頻率分別下降9.82%,28.80%,說(shuō)明儲(chǔ)液深度的增加使儲(chǔ)液結(jié)構(gòu)固有頻率下降,與以往杜永峰等和楊鳴等研究的結(jié)論一致。

    圖6 結(jié)構(gòu)加速度傳遞函數(shù)幅值-頻率特性曲線Fig.6 Amplitude-frequency characteristic curve of structural acceleration transfer function

    2.2 結(jié)構(gòu)振動(dòng)加速度響應(yīng)分析

    500 mm儲(chǔ)水、15.8g加載條件下,測(cè)點(diǎn)A2,A3,A4處結(jié)構(gòu)振動(dòng)加速度響應(yīng)典型時(shí)程曲線,如圖7所示。由圖7可知,受沖擊荷載后結(jié)構(gòu)振動(dòng)加速度迅速達(dá)到峰值,在低阻尼作用下,往復(fù)振動(dòng)幅值逐漸減小。加載后120 ms內(nèi)輸入加速度與不同測(cè)點(diǎn)處結(jié)構(gòu)加速度隨時(shí)間變化的規(guī)律,如圖8所示。由圖8可知,由振動(dòng)臺(tái)臺(tái)面A0測(cè)點(diǎn)至頂部A4測(cè)點(diǎn),加速度曲線峰值依次出現(xiàn),說(shuō)明地震動(dòng)影響由底部逐漸傳至頂部。為了描述地震動(dòng)從基礎(chǔ)傳遞到上部結(jié)構(gòu)的變化情況,采用傳遞系數(shù)β表征結(jié)構(gòu)加速度沿模型高度的放大或者衰減效應(yīng),β的定義為測(cè)點(diǎn)加速度峰值與臺(tái)面輸入加速度峰值的比值,即

    (a)

    (b)

    (c)圖7 結(jié)構(gòu)加速度響應(yīng)典型時(shí)程曲線Fig.7 Typical time history curve of structural acceleration response

    圖8 結(jié)構(gòu)加速度響應(yīng)規(guī)律對(duì)比Fig.8 The comparison of the law of structural acceleration response

    β=ai/a0

    (2)

    式中:ai為壁板上A2~A4測(cè)點(diǎn)加速度峰值;a0為臺(tái)面A0輸入加速度峰值。

    結(jié)構(gòu)加速度峰值與加速度傳遞系數(shù)沿壁板高度變化規(guī)律,如圖9、圖10所示。由圖9、圖10可知,由壁板底部至頂部,加速度峰值呈現(xiàn)非線性增大趨勢(shì),傳遞系數(shù)均大于1,說(shuō)明地震動(dòng)在由基礎(chǔ)向上部結(jié)構(gòu)傳遞的過(guò)程中存在遞增的放大效應(yīng),在結(jié)構(gòu)頂部放大效應(yīng)最為明顯。無(wú)水條件下,結(jié)構(gòu)加速度傳遞系數(shù)沿結(jié)構(gòu)高度呈加速增大趨勢(shì);儲(chǔ)水條件下,液面以下由于液體的附加質(zhì)量效應(yīng),結(jié)構(gòu)加速度放大效應(yīng)更加顯著。因此,以儲(chǔ)水液面為劃分界面,加速度傳遞系數(shù)沿壁板底部至頂部呈現(xiàn)先加速增大、后減速增大的趨勢(shì)。

    (a)

    (b)

    (c)圖9 加速度峰值沿壁板高度變化Fig.9 Peak acceleration varies along the height of the wall

    (a)

    (b)

    (c)圖10 結(jié)構(gòu)加速度傳遞系數(shù)沿壁板高度變化Fig.10 Acceleration transfer coefficient varies along the height of the wall

    不同工況下,結(jié)構(gòu)加速度傳遞系數(shù)隨輸入加速度變化曲線,如圖11所示??偟膩?lái)說(shuō),不同儲(chǔ)水條件下,隨著輸入加速度的增大,結(jié)構(gòu)振動(dòng)加速度放大效應(yīng)均得到一定提高。無(wú)水條件下,曲線變化較平緩,說(shuō)明隨著輸入加速度的增加,結(jié)構(gòu)加速度傳遞系數(shù)增幅不大;儲(chǔ)水條件下,曲線斜率增加,說(shuō)明結(jié)構(gòu)加速度傳遞系數(shù)受輸入加速度的影響增大,地震動(dòng)強(qiáng)化效應(yīng)更加明顯。為了直觀體現(xiàn)儲(chǔ)液狀態(tài)對(duì)加速度傳遞系數(shù)的影響,定義16.0g等級(jí)輸入加速度(無(wú)水15.9g、300 mm儲(chǔ)水16.3g、500 mm儲(chǔ)水15.8g)和19.2g等級(jí)輸入加速度(無(wú)水19.3g、300 mm儲(chǔ)水19.2g、500 mm儲(chǔ)水19.0g),該兩種等級(jí)輸入加速度下加速度傳遞系數(shù)隨儲(chǔ)水狀態(tài)變化規(guī)律,如圖12所示。由圖12可知,同一強(qiáng)度等級(jí)地震動(dòng)加載條件下,隨著儲(chǔ)液深度的增加,不同測(cè)點(diǎn)的結(jié)構(gòu)振動(dòng)加速度放大效應(yīng)均得到增強(qiáng)。

    圖11 結(jié)構(gòu)加速度傳遞系數(shù)隨輸入加速度變化Fig.11 Acceleration transfer coefficient varies with the input acceleration

    圖12 結(jié)構(gòu)加速度傳遞系數(shù)隨儲(chǔ)液狀態(tài)變化Fig.12 Acceleration transfer coefficient varies with the liquid storage conditions

    2.3 結(jié)構(gòu)應(yīng)變響應(yīng)分析

    500 mm儲(chǔ)水、15.8g加載條件下,測(cè)點(diǎn)S1,S2,S3處應(yīng)變響應(yīng)時(shí)程曲線,如圖13所示,可以發(fā)現(xiàn),受沖擊荷載后結(jié)構(gòu)迅速達(dá)到最大變形,經(jīng)阻尼消耗能量,往復(fù)振動(dòng)變形幅值逐漸減小。加載后120 ms內(nèi)不同測(cè)點(diǎn)處應(yīng)力隨時(shí)間變化的規(guī)律,如圖14所示,可以看出,不同測(cè)點(diǎn)處應(yīng)變的變化趨勢(shì)基本一致,由于應(yīng)變片布置于結(jié)構(gòu)內(nèi)側(cè),故加載后應(yīng)變首先達(dá)到負(fù)向極值,即結(jié)構(gòu)振動(dòng)變形首先呈現(xiàn)內(nèi)側(cè)受壓、外側(cè)受拉的狀態(tài)。

    (a)

    (b)

    (c)圖13 應(yīng)變響應(yīng)典型時(shí)程曲線Fig.13 Typical time history curve of strain response

    圖14 不同測(cè)點(diǎn)處應(yīng)變規(guī)律對(duì)比Fig.14 The comparison of the law of strain response

    不同儲(chǔ)液深度、不同輸入加速度下結(jié)構(gòu)應(yīng)變峰值沿壁板高度變化情況,如圖15所示??梢苑从吵觯虞d過(guò)程中結(jié)構(gòu)動(dòng)態(tài)應(yīng)力壁板底部至頂部呈現(xiàn)減速增大的趨勢(shì),在結(jié)構(gòu)頂部達(dá)到最大值。相比于無(wú)水結(jié)構(gòu),儲(chǔ)水結(jié)構(gòu)應(yīng)變分布發(fā)生變化,液面以下應(yīng)變?cè)鏊偬岣?,液面以上?yīng)變?cè)鏊俳档?,?dāng)儲(chǔ)液深度達(dá)到500 mm時(shí),結(jié)構(gòu)頂部S3測(cè)點(diǎn)處應(yīng)變與500 mm高度S2測(cè)點(diǎn)處應(yīng)變相差不大。

    (a)

    (b)

    (c)圖15 結(jié)構(gòu)應(yīng)變峰值沿壁板高度變化Fig.15 Peak strain varies along the height of the wall

    不同儲(chǔ)水條件下,結(jié)構(gòu)應(yīng)變隨輸入加速度變化曲線,如圖16所示。由圖16可知:隨著輸入加速度的提高,結(jié)構(gòu)應(yīng)變響應(yīng)逐漸增大,其中結(jié)構(gòu)頂部增幅更加明顯;與無(wú)水條件相比,儲(chǔ)水條件下曲線斜率增大,說(shuō)明應(yīng)變對(duì)于地震動(dòng)強(qiáng)度變化的敏感性提高。16.0g等級(jí)和19.2g等級(jí)輸入加速度下,結(jié)構(gòu)應(yīng)變隨儲(chǔ)液狀態(tài)變化,如圖17所示,說(shuō)明在同一等級(jí)地震動(dòng)加載條件下,儲(chǔ)液深度的增加使結(jié)構(gòu)變形值增大。以19.2g輸入加速度為例,無(wú)水、300 mm儲(chǔ)水、500 mm儲(chǔ)水條件下結(jié)構(gòu)頂部最大應(yīng)變值分別提高30.4%,45.9%。

    圖16 結(jié)構(gòu)應(yīng)變峰值隨輸入加速度變化Fig.16 Peak strain varies with the input acceleration

    圖17 結(jié)構(gòu)應(yīng)變峰值隨儲(chǔ)液狀態(tài)變化Fig.17 Peak strain varies with the liquid storage conditions

    2.4 結(jié)構(gòu)動(dòng)水壓力響應(yīng)分析

    由2.2節(jié)、2.3節(jié)分析可知,儲(chǔ)液狀態(tài)對(duì)于儲(chǔ)液結(jié)構(gòu)的動(dòng)力響應(yīng)影響很大,動(dòng)水壓力是關(guān)鍵因素。500 mm儲(chǔ)水、15.8g加載條件下,測(cè)點(diǎn)P1,P2,P3處動(dòng)水壓力響應(yīng)時(shí)程曲線,如圖18所示。由圖18可知,受沖擊荷載后作用于結(jié)構(gòu)壁板的動(dòng)水壓力迅速達(dá)到峰值,隨著結(jié)構(gòu)振動(dòng)衰減,動(dòng)水壓力幅值逐漸減小。此外,在800~1 000 ms段,動(dòng)水壓力傳感器監(jiān)測(cè)到的響應(yīng)值為零,這是由于爆炸沖擊震動(dòng)模擬平臺(tái)缺少限位裝置,振動(dòng)臺(tái)晃動(dòng)致結(jié)構(gòu)內(nèi)儲(chǔ)液體流向?qū)?cè)壁板。加載后120 ms內(nèi)不同測(cè)點(diǎn)處動(dòng)水壓力隨時(shí)間變化的規(guī)律,如圖19所示。由圖19可知,由結(jié)構(gòu)底部P1測(cè)點(diǎn)至上部P3測(cè)點(diǎn),動(dòng)水壓力首個(gè)峰值依次出現(xiàn),說(shuō)明動(dòng)水壓力響應(yīng)自結(jié)構(gòu)底部向頂部傳遞的過(guò)程中產(chǎn)生了相位差。

    (a)

    (b)

    (c)圖18 動(dòng)水壓力響應(yīng)典型時(shí)程曲線Fig.18 Typical time history curve of hydrodynamic pressure

    圖19 不同測(cè)點(diǎn)處動(dòng)水壓力規(guī)律對(duì)比Fig.19 The comparison of the law of hydrodynamic pressure

    不同儲(chǔ)液深度、不同輸入加速度下結(jié)構(gòu)動(dòng)水壓力峰值沿壁板高度變化情況,如圖20所示。由圖20可知,在500 mm儲(chǔ)水條件下,動(dòng)水壓力隨水深增加表現(xiàn)出減速增大的趨勢(shì),結(jié)構(gòu)底部動(dòng)水壓力最大,而2/5水位處與結(jié)構(gòu)底部動(dòng)水壓力相差不大,結(jié)構(gòu)臨水部分中下段承受了水平相當(dāng)且幅值較高的動(dòng)水壓力。

    (a)

    (b)圖20 結(jié)構(gòu)動(dòng)水壓力沿壁板高度變化Fig.20 Hydrodynamic pressure varies along the height of the wall

    不同工況下,結(jié)構(gòu)動(dòng)水壓力隨輸入加速度變化曲線,如圖21所示。由圖21可知:隨著輸入加速度的增加,結(jié)構(gòu)承受的動(dòng)水壓力逐漸增大;結(jié)構(gòu)儲(chǔ)水由300 mm增至500 mm,在不同加速度加載的條件下,動(dòng)水壓力均有明顯的增幅。

    圖21 結(jié)構(gòu)動(dòng)水壓力隨輸入加速度變化Fig.21 Hydrodynamic pressure varies with the input acceleration

    3 理論分析與討論

    3.1 儲(chǔ)液結(jié)構(gòu)壁板振型探討

    由引言部分可知,以往學(xué)者對(duì)于儲(chǔ)液結(jié)構(gòu)自振特性開展了相關(guān)研究,關(guān)于壁板振型尚未得到一致結(jié)論。本文基于試驗(yàn)數(shù)據(jù),對(duì)于儲(chǔ)液結(jié)構(gòu)壁板振型進(jìn)行初步的探討。儲(chǔ)液結(jié)構(gòu)計(jì)算模型示意圖如圖22所示,圖22中:地震動(dòng)沿x方向;Lx為矩形儲(chǔ)液結(jié)構(gòu)沿地震作用方向的邊長(zhǎng);Ly為結(jié)構(gòu)垂直地震作用方向的邊長(zhǎng);Hs為結(jié)構(gòu)高度;Hw為儲(chǔ)液深度。

    圖22 儲(chǔ)液結(jié)構(gòu)計(jì)算模型示意圖Fig.22 Schematic diagram of LSS calculation model

    由板殼振動(dòng)理論[33]可知,矩形板固有振動(dòng)分析可采用雙向梁函數(shù)組合級(jí)數(shù)逼近方法,一般性矩形薄板的撓度振型W(y,z)可表示為

    (3)

    式中:Ym(y)與Zn(z)分別為與壁板水平y(tǒng)方向,豎直z方向兩端邊界條件對(duì)應(yīng)的第m與n階梁振型函數(shù);Amn為待定系數(shù),用于調(diào)整梁函數(shù)的組合以逼近矩形板振型真實(shí)解。矩形薄板的振動(dòng)方程w(y,z,t)可表示為

    w(y,z,t)=W(y,z)T(t)

    (4)

    式中,T(t)為時(shí)間相關(guān)的動(dòng)力函數(shù)。因此,得到板內(nèi)橫向應(yīng)變分量εy為

    (5)

    將式(3)、式(4)代入可得

    (6)

    式中,t為壁板厚度。考慮到試驗(yàn)中單向應(yīng)變片沿壁板水平y(tǒng)方向布置,且均粘貼處于壁板中軸位置(y=0),故由式(6)可知S1,S2,S3測(cè)點(diǎn)應(yīng)變值對(duì)比結(jié)果為

    εy(S3):εy(S2):εy(S1)=Zn(S3):Zn(S2):Zn(S1)

    (7)

    該結(jié)果反映了結(jié)構(gòu)壁板豎向振型。因此,將圖15應(yīng)變測(cè)量結(jié)果進(jìn)行歸一化處理,結(jié)果如圖23所示。

    圖23 結(jié)構(gòu)應(yīng)變峰值沿壁板高度變化歸一化結(jié)果Fig.23 Normalized results of peak strain along the height of wall

    圖23中除歸一化的試驗(yàn)數(shù)據(jù)外,標(biāo)注了彎曲懸臂梁一階振型曲線與剪切懸臂梁一階振型曲線。可以看出,不同工況下試驗(yàn)值與剪切懸臂梁振型曲線更接近,說(shuō)明儲(chǔ)液結(jié)構(gòu)壁板振動(dòng)豎向振型更加符合剪切型變形理論。由無(wú)水至300 mm儲(chǔ)水、500 mm儲(chǔ)水,液體附加質(zhì)量效應(yīng)出現(xiàn)并逐漸增強(qiáng),壁板振型在剪切懸臂變形的基礎(chǔ)上發(fā)生局部變化,表現(xiàn)為液面以下壁板臨水部分變形得到增強(qiáng),上部無(wú)水部分變形相對(duì)弱化。

    3.2 動(dòng)水壓力理論計(jì)算與規(guī)范對(duì)比

    考慮到動(dòng)水壓力對(duì)于儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)的重要影響,抗震設(shè)計(jì)時(shí)需對(duì)爆炸地震動(dòng)下的動(dòng)水壓力進(jìn)行驗(yàn)算。首先進(jìn)行理論方法求解,計(jì)算模型參考圖22。基于程選生等推導(dǎo)的動(dòng)水壓力計(jì)算方法,在考慮儲(chǔ)液結(jié)構(gòu)彈性變形的條件下,壁板受到的脈沖動(dòng)水壓力pi由兩部分組成:剛體運(yùn)動(dòng)產(chǎn)生的脈沖壓力pg與考慮壁板側(cè)向變形產(chǎn)生的脈沖壓力pf,可按式(8)~式(11)計(jì)算

    (8)

    (9)

    (10)

    (11)

    (12)

    式中,βmax為單質(zhì)點(diǎn)彈性結(jié)構(gòu)的最大水平加速度與地震動(dòng)最大加速度的統(tǒng)計(jì)平均值的比值[34],即

    (13)

    對(duì)于鋼筋混凝土儲(chǔ)液結(jié)構(gòu)取值2.25。對(duì)于本文不銹鋼儲(chǔ)液結(jié)構(gòu),以2.2節(jié)加速度傳遞系數(shù)統(tǒng)計(jì)平均值進(jìn)行修正,取為2.85。

    此外,基于Hashemi等提出的流固耦合理論方法,壁板對(duì)流壓力pc可按式(14)~式(19)計(jì)算

    (14)

    (15)

    (16)

    λj=(2j+1)π/2

    (17)

    αj=λj/Lx

    (18)

    (19)

    以500 mm儲(chǔ)水、擺錘高度H3工況為例,通過(guò)Maple軟件進(jìn)行理論計(jì)算前,以Asym2Sig模型[見式(20)]對(duì)地震動(dòng)加速度信號(hào)進(jìn)行單峰值擬合,結(jié)果如圖24所示,擬合曲線與試驗(yàn)曲線吻合較好。

    圖24 典型輸入加速度信號(hào)與Asym2Sig模型擬合結(jié)果Fig.24 The typical input acceleration signal and the fitting results using Asym2Sig model

    (20)

    將擬合函數(shù)代入式(8)~式(19),計(jì)算得到該工況下動(dòng)水壓力理論值,計(jì)算結(jié)果如圖25~圖26所示。圖25中呈現(xiàn)了地震動(dòng)加速度峰值時(shí)刻脈沖壓力pi、剛體運(yùn)動(dòng)脈沖壓力pg、壁板變形脈沖壓力pf、以及對(duì)流壓力pc沿臨水壁板分布情況??梢钥闯觯阂好嬉韵码S著水深增加,pi,pf先增大后減?。籶g逐漸增大,在結(jié)構(gòu)底部達(dá)到最大值;pc逐漸減小,最大值出現(xiàn)在液面處。相比之下,液體晃動(dòng)產(chǎn)生的對(duì)流壓力不足脈沖壓力的1%,對(duì)動(dòng)水壓力影響可忽略。脈沖壓力則是動(dòng)水壓力主要組成部分,對(duì)儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)起主要作用。

    注:脈沖壓力pi,pf,pg參照下軸;對(duì)流壓力pc參照上軸。圖25 理論計(jì)算動(dòng)水壓力沿壁板分布Fig.25 Theoretically calculated hydrodynamic pressure distribution along the wall

    (a) 300 mm儲(chǔ)水

    (b) 500 mm儲(chǔ)水圖26 動(dòng)水壓力理論計(jì)算與試驗(yàn)結(jié)果對(duì)比Fig.26 Comparisons between theoretical calculations and experimental results of hydrodynamic pressure

    對(duì)比不同測(cè)點(diǎn)處動(dòng)水壓力最大值理論計(jì)算與試驗(yàn)結(jié)果,如圖26所示??梢园l(fā)現(xiàn)兩者存在一定偏差,分析其原因,一是式(8)~式(11)計(jì)算pg與pf時(shí)對(duì)結(jié)構(gòu)振型進(jìn)行簡(jiǎn)化,弱化了結(jié)構(gòu)響應(yīng)的非線性;二是式(12)計(jì)算pi時(shí)作了簡(jiǎn)化處理,以便于工程設(shè)計(jì)。但總的來(lái)說(shuō),理論計(jì)算值在結(jié)構(gòu)底部P1,上部P3處吻合較好,不同測(cè)點(diǎn)處計(jì)算偏于保守,有利于結(jié)構(gòu)安全設(shè)計(jì)。加之工程設(shè)計(jì)往往關(guān)注計(jì)算求解的便利性,上述理論計(jì)算方法仍具有較好的適用性,下面結(jié)合現(xiàn)有設(shè)計(jì)規(guī)范進(jìn)行討論。

    現(xiàn)行規(guī)范中提供了盛水構(gòu)筑物動(dòng)水壓力的抗震驗(yàn)算方法。其中,GB 50032—2003《室外給水排水和燃?xì)鉄崃こ炭拐鹪O(shè)計(jì)規(guī)范》[35]在進(jìn)行矩形水池抗震驗(yàn)算時(shí),將壁板受到的動(dòng)水壓力簡(jiǎn)化為沿平面、沿高度等值分布,并給出地面式矩形水池在水平地震作用下動(dòng)水壓力標(biāo)準(zhǔn)值Fwr,c的計(jì)算方法

    Fwr,c=KH·γw·Hw·fwr

    (21)

    式中:KH為水平地震加速度與重力加速度的比值;γw為池水重力密度;fwr為地面式矩形水池動(dòng)水壓力系數(shù),可按表2采用。

    表2 地面式水池動(dòng)水壓力系數(shù)fwrTab.2 Hydrodynamic pressure coefficient of above-ground pool, fwr

    參照規(guī)范計(jì)算式(21)對(duì)動(dòng)水壓力理論計(jì)算式(12)進(jìn)行改進(jìn),得到

    Fwr,c2=KH·γw·Hw·fwr2

    (22)

    (23)

    將式(23)進(jìn)一步化簡(jiǎn),得到

    (24)

    式(24)表達(dá)了fwr2與Lx/Hw的函數(shù)關(guān)系,對(duì)應(yīng)函數(shù)曲線如圖27所示??梢钥闯?,對(duì)于本文地面式不銹鋼矩形儲(chǔ)液結(jié)構(gòu),GB 50032—2003《室外給水排水和燃?xì)鉄崃こ炭拐鹪O(shè)計(jì)規(guī)范》系數(shù)fwr偏差較大且不利,需要進(jìn)行修正。參考規(guī)范中常見做法進(jìn)行離散化處理,得到修正后的動(dòng)水壓力系數(shù)fwr2取值表,見表3,可供相關(guān)工程按照線性差值法使用。

    圖27 動(dòng)水壓力系數(shù)修正結(jié)果Fig.27 Modified result of hydrodynamic pressure coefficient

    表3 修正后的地面式不銹鋼儲(chǔ)液結(jié)構(gòu)動(dòng)水壓力系數(shù)fwr2Tab.3 Modified hydrodynamic pressure coefficient for ground-mounted stainless steel rectangular LSS, fwr2

    此外,在進(jìn)行儲(chǔ)液結(jié)構(gòu)抗震設(shè)計(jì)時(shí),需要考慮動(dòng)水壓力所產(chǎn)生的傾覆力矩。在推導(dǎo)動(dòng)水壓力計(jì)算式的基礎(chǔ)上,作用于儲(chǔ)液結(jié)構(gòu)壁板上的動(dòng)水壓力所導(dǎo)致的傾覆力矩可計(jì)算為

    (25)

    聯(lián)立式(12)并參照規(guī)范計(jì)算式(21)進(jìn)行改進(jìn),得到

    (26)

    (27)

    式(27)表達(dá)了傾覆力矩系數(shù)mover與Lx/Hw的函數(shù)關(guān)系,對(duì)應(yīng)函數(shù)曲線如圖28所示,經(jīng)離散化處理得到的動(dòng)水壓力傾覆力矩系數(shù)mover取值表,見表4,可供相關(guān)工程按照線性差值法使用。

    圖28 動(dòng)水壓力傾覆力矩系數(shù)計(jì)算結(jié)果Fig.28 Calculation result of overturning moment coefficient

    表4 地面式不銹鋼儲(chǔ)液結(jié)構(gòu)傾覆力矩系數(shù)moverTab.4 Overturning moment coefficient for ground-mounted stainless steel rectangular LSS, mover

    4 結(jié) 論

    本文利用爆炸沖擊震動(dòng)模擬平臺(tái)對(duì)固有頻率為54.15 Hz的地面式矩形無(wú)頂蓋儲(chǔ)液結(jié)構(gòu)開展了動(dòng)力響應(yīng)試驗(yàn)研究,得到了峰值加速度范圍為10.0~26.4g、平均持時(shí)為14.8 ms的地震動(dòng)荷載,分析了儲(chǔ)液狀態(tài)、地震動(dòng)強(qiáng)度對(duì)儲(chǔ)液結(jié)構(gòu)動(dòng)力響應(yīng)的影響,結(jié)合試驗(yàn)與理論方法揭示了儲(chǔ)液結(jié)構(gòu)振型特點(diǎn)和動(dòng)水壓力變化規(guī)律,主要結(jié)論有:

    (1) 通過(guò)分析傳遞函數(shù)得到不同儲(chǔ)液狀態(tài)下儲(chǔ)液結(jié)構(gòu)的固有頻率,發(fā)現(xiàn)儲(chǔ)液狀態(tài)影響液-固體系固有頻率,由無(wú)水到300 mm儲(chǔ)水、500 mm儲(chǔ)水狀態(tài),液-固體系固有頻率分別下降9.82%,28.80%。

    (2) 儲(chǔ)液結(jié)構(gòu)振動(dòng)加速度沿壁板高度呈放大效應(yīng),地震動(dòng)強(qiáng)度的提高使放大效應(yīng)更加明顯。相比于無(wú)水條件,儲(chǔ)水條件下加速度傳遞系數(shù)對(duì)地震動(dòng)強(qiáng)度的敏感性提高。同一等級(jí)地震動(dòng)強(qiáng)度下,隨著儲(chǔ)液深度增加,結(jié)構(gòu)加速度放大效應(yīng)增強(qiáng)。

    (3) 儲(chǔ)液結(jié)構(gòu)應(yīng)變沿壁板高度呈現(xiàn)減速增大的趨勢(shì),壁板振動(dòng)豎向振型符合懸臂梁剪切變形理論。儲(chǔ)水條件下,壁板臨水部分變形由于液體附加質(zhì)量效應(yīng)得到加強(qiáng),無(wú)水部分變形相對(duì)弱化。結(jié)構(gòu)變形隨地震動(dòng)增強(qiáng)而增加,結(jié)構(gòu)頂部應(yīng)變值的地震動(dòng)強(qiáng)度敏感性較強(qiáng)。同一等級(jí)地震動(dòng)強(qiáng)度下,儲(chǔ)水深度的增加使結(jié)構(gòu)變形幅度提高。

    (4) 儲(chǔ)液結(jié)構(gòu)動(dòng)水壓力隨水深增加呈現(xiàn)先增大后減小的趨勢(shì),相較于對(duì)流壓力,脈沖壓力是動(dòng)水壓力主要組成部分。分別以流固耦合理論方法和現(xiàn)行設(shè)計(jì)規(guī)范計(jì)算得到了水平地震作用下動(dòng)水壓力標(biāo)準(zhǔn)值,發(fā)現(xiàn)現(xiàn)有規(guī)范對(duì)地面式矩形不銹鋼儲(chǔ)液結(jié)構(gòu)設(shè)計(jì)不利,結(jié)合理論方法修正了動(dòng)水壓力系數(shù)。推導(dǎo)了動(dòng)水壓力所致傾覆力矩的計(jì)算公式,并得到了傾覆力矩系數(shù)。

    本文采用模型試驗(yàn)方法進(jìn)行參數(shù)研究,并結(jié)合理論分析得到了一些結(jié)論。未來(lái)還需要對(duì)針對(duì)不同結(jié)構(gòu)尺寸、形狀等進(jìn)一步研究,以擴(kuò)充認(rèn)知體系,得出更具普遍指導(dǎo)意義的結(jié)論。

    猜你喜歡
    動(dòng)水儲(chǔ)液儲(chǔ)水
    蝶閥動(dòng)水力矩計(jì)算方法辨析
    一種橡膠制品的硫化裝置
    只給1升水強(qiáng)迫節(jié)約
    基于ADINA的矩形儲(chǔ)液池流固耦合動(dòng)力響應(yīng)分析
    山西建筑(2016年25期)2016-09-28 09:59:00
    9省市12批次儲(chǔ)水式電熱水器抽查不合格
    直液式儲(chǔ)水器原理
    糯扎渡水電站筒閥動(dòng)水關(guān)閉試驗(yàn)與分析
    錫林郭勒盟牧區(qū)儲(chǔ)水窖工程建設(shè)探究
    機(jī)載儲(chǔ)液罐研制方法研究
    動(dòng)水壓力對(duì)深水橋梁地震響應(yīng)的影響
    小蜜桃在线观看免费完整版高清| 又爽又黄无遮挡网站| 又紧又爽又黄一区二区| 老司机深夜福利视频在线观看| 免费大片18禁| 成年版毛片免费区| 免费不卡的大黄色大毛片视频在线观看 | 天天躁日日操中文字幕| 国产高清视频在线播放一区| av女优亚洲男人天堂| 一区福利在线观看| 免费黄网站久久成人精品| 国产精品人妻久久久影院| 全区人妻精品视频| 国产精品伦人一区二区| 久久人人爽人人爽人人片va| 亚洲无线观看免费| 热99re8久久精品国产| 日韩欧美三级三区| 久久久久国产精品人妻aⅴ院| 国产精品女同一区二区软件 | 伦精品一区二区三区| 日韩av在线大香蕉| 五月伊人婷婷丁香| 国产探花极品一区二区| 少妇人妻精品综合一区二区 | 天堂动漫精品| 精品久久久久久久久久免费视频| 国产午夜精品久久久久久一区二区三区 | 极品教师在线免费播放| 高清毛片免费观看视频网站| 免费av观看视频| 春色校园在线视频观看| 99国产极品粉嫩在线观看| 亚洲,欧美,日韩| 人妻少妇偷人精品九色| 搡女人真爽免费视频火全软件 | av黄色大香蕉| 香蕉av资源在线| 亚洲性久久影院| 欧美激情在线99| 日本免费一区二区三区高清不卡| 亚洲av不卡在线观看| 亚洲精品粉嫩美女一区| 欧美人与善性xxx| 国产人妻一区二区三区在| 国产色婷婷99| 亚洲中文字幕日韩| 国产麻豆成人av免费视频| 一个人免费在线观看电影| 观看免费一级毛片| 我要搜黄色片| 可以在线观看毛片的网站| 欧美黑人欧美精品刺激| 天天一区二区日本电影三级| 国产欧美日韩一区二区精品| 精品人妻1区二区| 欧美国产日韩亚洲一区| 国产精品爽爽va在线观看网站| 午夜视频国产福利| 一级毛片久久久久久久久女| 成年女人毛片免费观看观看9| 国产真实伦视频高清在线观看 | 69av精品久久久久久| 波多野结衣巨乳人妻| 淫秽高清视频在线观看| 欧美+日韩+精品| 亚洲av第一区精品v没综合| 日本成人三级电影网站| 国内精品宾馆在线| 成人毛片a级毛片在线播放| 久久久久精品国产欧美久久久| 久久精品影院6| 国产精品99久久久久久久久| 国产精华一区二区三区| 国产精品人妻久久久影院| 精品日产1卡2卡| 欧美性猛交╳xxx乱大交人| 国产三级在线视频| 日韩一区二区视频免费看| 日韩 亚洲 欧美在线| av国产免费在线观看| 小蜜桃在线观看免费完整版高清| 夜夜爽天天搞| 日韩欧美三级三区| 久久精品夜夜夜夜夜久久蜜豆| 人妻夜夜爽99麻豆av| 中文亚洲av片在线观看爽| 老熟妇乱子伦视频在线观看| 欧美3d第一页| 两个人视频免费观看高清| 国内揄拍国产精品人妻在线| 国产视频一区二区在线看| 色尼玛亚洲综合影院| 国产精品女同一区二区软件 | 看十八女毛片水多多多| 如何舔出高潮| 在线免费十八禁| 久久精品影院6| 国产三级在线视频| 女的被弄到高潮叫床怎么办 | 婷婷色综合大香蕉| 99国产精品一区二区蜜桃av| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久国内视频| 神马国产精品三级电影在线观看| 内地一区二区视频在线| 舔av片在线| 女同久久另类99精品国产91| 老司机福利观看| 麻豆国产97在线/欧美| 精品不卡国产一区二区三区| 99久久成人亚洲精品观看| 日本 欧美在线| 亚洲精品粉嫩美女一区| 男女啪啪激烈高潮av片| 久久国产精品人妻蜜桃| 亚洲aⅴ乱码一区二区在线播放| 午夜福利欧美成人| 国产熟女欧美一区二区| 久久99热6这里只有精品| 免费电影在线观看免费观看| 婷婷精品国产亚洲av| 久久久久精品国产欧美久久久| 精品久久久久久久久av| 久久久久久九九精品二区国产| 黄色丝袜av网址大全| 久久6这里有精品| 丰满人妻一区二区三区视频av| 欧美三级亚洲精品| 欧美最黄视频在线播放免费| 久久久色成人| 欧美日韩精品成人综合77777| 精品久久久久久久人妻蜜臀av| 久99久视频精品免费| 国产欧美日韩精品亚洲av| 熟女电影av网| 日韩国内少妇激情av| 久久人妻av系列| 美女高潮的动态| 岛国在线免费视频观看| 3wmmmm亚洲av在线观看| 色综合婷婷激情| 中文字幕免费在线视频6| 别揉我奶头~嗯~啊~动态视频| 日韩中字成人| 三级国产精品欧美在线观看| 免费看光身美女| 九九爱精品视频在线观看| 久久精品国产鲁丝片午夜精品 | 老熟妇乱子伦视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 不卡视频在线观看欧美| 国产精品1区2区在线观看.| 在线天堂最新版资源| 精品一区二区三区视频在线观看免费| 欧美人与善性xxx| 夜夜看夜夜爽夜夜摸| 人人妻人人看人人澡| 99热这里只有是精品50| 国产69精品久久久久777片| 国产日本99.免费观看| 一级av片app| 久久久精品大字幕| 成人一区二区视频在线观看| 精品人妻一区二区三区麻豆 | 变态另类丝袜制服| 欧美bdsm另类| 亚洲美女视频黄频| 制服丝袜大香蕉在线| 成年人黄色毛片网站| 麻豆国产av国片精品| 男人狂女人下面高潮的视频| 中文字幕高清在线视频| 在线观看66精品国产| 国产精品一区二区性色av| 成人三级黄色视频| 91麻豆精品激情在线观看国产| 男女做爰动态图高潮gif福利片| 村上凉子中文字幕在线| 婷婷精品国产亚洲av在线| 婷婷丁香在线五月| 中文字幕免费在线视频6| 亚洲欧美日韩高清在线视频| 在线观看舔阴道视频| 看黄色毛片网站| 色5月婷婷丁香| 成人特级黄色片久久久久久久| 黄色日韩在线| 国产一级毛片七仙女欲春2| 精品人妻1区二区| 男女下面进入的视频免费午夜| 五月玫瑰六月丁香| 琪琪午夜伦伦电影理论片6080| 欧美日韩乱码在线| www.www免费av| 国内精品久久久久久久电影| 蜜桃亚洲精品一区二区三区| 观看免费一级毛片| 国产在线男女| 国语自产精品视频在线第100页| 伊人久久精品亚洲午夜| 99久久九九国产精品国产免费| 99久久九九国产精品国产免费| 91狼人影院| 在线天堂最新版资源| 精品人妻偷拍中文字幕| 久久这里只有精品中国| 五月玫瑰六月丁香| 亚洲七黄色美女视频| 啪啪无遮挡十八禁网站| 99在线视频只有这里精品首页| 国产亚洲精品久久久com| 18禁黄网站禁片午夜丰满| 婷婷丁香在线五月| 久久久久久久亚洲中文字幕| 成人欧美大片| 国语自产精品视频在线第100页| 中文亚洲av片在线观看爽| 日本黄色视频三级网站网址| 国产精品亚洲一级av第二区| 国产极品精品免费视频能看的| 日本黄大片高清| av中文乱码字幕在线| 日本 av在线| 国产一区二区三区在线臀色熟女| 69av精品久久久久久| 午夜激情福利司机影院| 日本在线视频免费播放| 人妻丰满熟妇av一区二区三区| 少妇的逼水好多| 一区二区三区激情视频| 3wmmmm亚洲av在线观看| 亚洲av第一区精品v没综合| 一进一出抽搐动态| 两人在一起打扑克的视频| 国产精品永久免费网站| 日韩欧美在线乱码| 成年版毛片免费区| 亚洲精品日韩av片在线观看| 高清日韩中文字幕在线| 我要看日韩黄色一级片| 岛国在线免费视频观看| 最近在线观看免费完整版| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久久久| 色哟哟哟哟哟哟| 午夜免费成人在线视频| 美女黄网站色视频| 国产高清激情床上av| 国产精品精品国产色婷婷| 日韩亚洲欧美综合| 亚洲人成网站在线播| 91狼人影院| 国产极品精品免费视频能看的| 日韩精品青青久久久久久| 日日撸夜夜添| 免费大片18禁| 国产成人aa在线观看| 久久久久久久亚洲中文字幕| 99视频精品全部免费 在线| 在线观看一区二区三区| 超碰av人人做人人爽久久| 在线看三级毛片| 亚洲三级黄色毛片| 亚洲,欧美,日韩| 日本与韩国留学比较| 亚洲av美国av| 欧美在线一区亚洲| 欧美色欧美亚洲另类二区| 亚洲欧美日韩高清在线视频| av福利片在线观看| 精品国产三级普通话版| 日本爱情动作片www.在线观看 | 男女边吃奶边做爰视频| 久久久精品欧美日韩精品| 亚洲成人久久爱视频| 亚洲三级黄色毛片| 日本成人三级电影网站| 五月伊人婷婷丁香| 狂野欧美激情性xxxx在线观看| 春色校园在线视频观看| 国产精品永久免费网站| 欧美日韩瑟瑟在线播放| 悠悠久久av| 长腿黑丝高跟| 免费人成视频x8x8入口观看| 欧美丝袜亚洲另类 | 夜夜看夜夜爽夜夜摸| 3wmmmm亚洲av在线观看| 国产精品三级大全| 给我免费播放毛片高清在线观看| 天堂网av新在线| 久久久久性生活片| 可以在线观看毛片的网站| 一区二区三区高清视频在线| 国产亚洲91精品色在线| 亚洲av免费高清在线观看| 国产色爽女视频免费观看| 成人特级黄色片久久久久久久| 亚洲美女黄片视频| 成年女人永久免费观看视频| 亚洲精品久久国产高清桃花| 如何舔出高潮| 国产白丝娇喘喷水9色精品| 一级a爱片免费观看的视频| 国产精品爽爽va在线观看网站| 日日干狠狠操夜夜爽| 欧美xxxx黑人xx丫x性爽| 别揉我奶头 嗯啊视频| 婷婷亚洲欧美| 欧美bdsm另类| 亚洲性夜色夜夜综合| 夜夜夜夜夜久久久久| 亚洲国产日韩欧美精品在线观看| 亚洲中文日韩欧美视频| 亚洲中文字幕一区二区三区有码在线看| 国内揄拍国产精品人妻在线| 精品免费久久久久久久清纯| 亚洲人成伊人成综合网2020| 免费观看精品视频网站| 亚洲电影在线观看av| 黄色女人牲交| 最近在线观看免费完整版| 联通29元200g的流量卡| 夜夜看夜夜爽夜夜摸| 日本精品一区二区三区蜜桃| 中亚洲国语对白在线视频| 亚洲av电影不卡..在线观看| 国产精品99久久久久久久久| 午夜影院日韩av| а√天堂www在线а√下载| 亚洲欧美精品综合久久99| 亚洲色图av天堂| 久久婷婷人人爽人人干人人爱| 日韩欧美免费精品| 99热6这里只有精品| 日日摸夜夜添夜夜添av毛片 | 亚洲精华国产精华精| 麻豆精品久久久久久蜜桃| 麻豆一二三区av精品| 12—13女人毛片做爰片一| 亚洲内射少妇av| 可以在线观看的亚洲视频| 亚洲性久久影院| 婷婷色综合大香蕉| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美精品综合久久99| 在线播放国产精品三级| 精品99又大又爽又粗少妇毛片 | 舔av片在线| 性欧美人与动物交配| 日韩一区二区视频免费看| 成熟少妇高潮喷水视频| 久久久久精品国产欧美久久久| 国产乱人伦免费视频| 精品午夜福利在线看| 狂野欧美白嫩少妇大欣赏| 少妇熟女aⅴ在线视频| 深夜精品福利| 97超级碰碰碰精品色视频在线观看| 亚洲国产日韩欧美精品在线观看| 国产av在哪里看| 少妇裸体淫交视频免费看高清| 国产亚洲欧美98| 又黄又爽又免费观看的视频| 国产91精品成人一区二区三区| 国产欧美日韩一区二区精品| 搞女人的毛片| 美女高潮喷水抽搐中文字幕| 久9热在线精品视频| 91在线观看av| 久久这里只有精品中国| 麻豆成人午夜福利视频| 亚洲无线观看免费| 丰满的人妻完整版| bbb黄色大片| 内射极品少妇av片p| 久久中文看片网| 麻豆av噜噜一区二区三区| 精品一区二区三区人妻视频| 丝袜美腿在线中文| 国内精品久久久久精免费| 国产男靠女视频免费网站| 男女视频在线观看网站免费| 黄色欧美视频在线观看| 国产三级在线视频| 99riav亚洲国产免费| 国产精品野战在线观看| 婷婷精品国产亚洲av在线| av视频在线观看入口| 乱码一卡2卡4卡精品| 99国产极品粉嫩在线观看| 男女视频在线观看网站免费| 成人一区二区视频在线观看| 日韩欧美一区二区三区在线观看| 欧美日韩乱码在线| 最新在线观看一区二区三区| 日本成人三级电影网站| 中文字幕熟女人妻在线| 直男gayav资源| 小说图片视频综合网站| 久久久久久大精品| 成年女人毛片免费观看观看9| 国产精品1区2区在线观看.| 国产黄a三级三级三级人| 欧美丝袜亚洲另类 | 村上凉子中文字幕在线| 偷拍熟女少妇极品色| 色av中文字幕| 全区人妻精品视频| 久久久久久大精品| 欧美极品一区二区三区四区| 两个人的视频大全免费| 精品免费久久久久久久清纯| 久久久精品大字幕| 色在线成人网| 久久午夜亚洲精品久久| 亚洲欧美日韩东京热| 欧美不卡视频在线免费观看| 12—13女人毛片做爰片一| 中文字幕熟女人妻在线| 国产精品伦人一区二区| 亚洲第一区二区三区不卡| 淫秽高清视频在线观看| 国产在视频线在精品| 日日啪夜夜撸| a在线观看视频网站| 亚洲中文日韩欧美视频| 黄片wwwwww| 简卡轻食公司| 欧美日韩瑟瑟在线播放| 国产一区二区在线av高清观看| 美女 人体艺术 gogo| 一个人看视频在线观看www免费| 91麻豆精品激情在线观看国产| 男插女下体视频免费在线播放| 久久久精品大字幕| av国产免费在线观看| 麻豆成人午夜福利视频| 国产精品不卡视频一区二区| 亚洲综合色惰| 美女高潮的动态| 人妻夜夜爽99麻豆av| 九九热线精品视视频播放| 国产v大片淫在线免费观看| 成人国产一区最新在线观看| 在线观看美女被高潮喷水网站| 亚洲欧美激情综合另类| 伦理电影大哥的女人| 国产一区二区激情短视频| 久久久久性生活片| 日韩欧美 国产精品| 狠狠狠狠99中文字幕| 国产色婷婷99| 欧美+亚洲+日韩+国产| 99久久无色码亚洲精品果冻| 在线观看免费视频日本深夜| 91午夜精品亚洲一区二区三区 | 亚洲av电影不卡..在线观看| 成人av一区二区三区在线看| 亚洲七黄色美女视频| 麻豆一二三区av精品| 亚洲av美国av| 国产久久久一区二区三区| 国产三级在线视频| 人妻少妇偷人精品九色| 黄色配什么色好看| 国产一区二区亚洲精品在线观看| 91精品国产九色| 97碰自拍视频| 亚洲四区av| 成年版毛片免费区| 露出奶头的视频| 中文在线观看免费www的网站| 1000部很黄的大片| 中国美女看黄片| 伊人久久精品亚洲午夜| 国产精品日韩av在线免费观看| 成人鲁丝片一二三区免费| 亚洲va日本ⅴa欧美va伊人久久| 久久久成人免费电影| 又爽又黄无遮挡网站| 国产爱豆传媒在线观看| 性色avwww在线观看| 欧美精品啪啪一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美在线二视频| 国产精品无大码| 国产成人aa在线观看| 不卡一级毛片| 国产成人影院久久av| 伦理电影大哥的女人| 美女xxoo啪啪120秒动态图| 九色国产91popny在线| 成人国产麻豆网| 最新在线观看一区二区三区| 成人高潮视频无遮挡免费网站| 日本 av在线| 成人鲁丝片一二三区免费| 黄色一级大片看看| 亚洲18禁久久av| 中文在线观看免费www的网站| 国产精品无大码| 国产视频内射| 国内少妇人妻偷人精品xxx网站| 色播亚洲综合网| 久久精品国产亚洲av涩爱 | 日日摸夜夜添夜夜添小说| 久久精品夜夜夜夜夜久久蜜豆| 精品无人区乱码1区二区| 成人av一区二区三区在线看| 国产精华一区二区三区| 国产视频一区二区在线看| 在线观看午夜福利视频| 亚洲av二区三区四区| 丝袜美腿在线中文| 少妇熟女aⅴ在线视频| 极品教师在线免费播放| 赤兔流量卡办理| 哪里可以看免费的av片| 别揉我奶头 嗯啊视频| www.色视频.com| 内地一区二区视频在线| 一级a爱片免费观看的视频| 日本撒尿小便嘘嘘汇集6| 亚洲精品亚洲一区二区| 国产中年淑女户外野战色| 成人av一区二区三区在线看| 亚洲四区av| 熟女电影av网| 欧美成人a在线观看| 日韩精品中文字幕看吧| 观看美女的网站| 尾随美女入室| 亚洲性夜色夜夜综合| 十八禁国产超污无遮挡网站| 色吧在线观看| 精品99又大又爽又粗少妇毛片 | 国产人妻一区二区三区在| 中文亚洲av片在线观看爽| 最近最新免费中文字幕在线| 伊人久久精品亚洲午夜| 我要搜黄色片| 久久精品国产亚洲av涩爱 | 赤兔流量卡办理| 亚洲国产欧美人成| 成人欧美大片| 久久久久久国产a免费观看| 少妇的逼水好多| 少妇被粗大猛烈的视频| 人妻久久中文字幕网| 精品久久久久久久久亚洲 | 在线天堂最新版资源| 国产精品久久视频播放| 亚洲av五月六月丁香网| 99国产精品一区二区蜜桃av| 亚洲国产精品合色在线| 日本爱情动作片www.在线观看 | 精品人妻熟女av久视频| 一个人看的www免费观看视频| 免费观看人在逋| 内射极品少妇av片p| 偷拍熟女少妇极品色| 亚洲,欧美,日韩| 男女之事视频高清在线观看| 亚洲性久久影院| 制服丝袜大香蕉在线| 天天一区二区日本电影三级| 欧美性感艳星| 国产伦一二天堂av在线观看| 我要搜黄色片| 成人高潮视频无遮挡免费网站| 深爱激情五月婷婷| 窝窝影院91人妻| 狂野欧美激情性xxxx在线观看| 国产色婷婷99| 变态另类丝袜制服| 亚洲欧美日韩高清在线视频| 十八禁网站免费在线| 成人一区二区视频在线观看| 国产精品三级大全| 成人欧美大片| 久久久精品欧美日韩精品| 成人永久免费在线观看视频| 狠狠狠狠99中文字幕| 老熟妇乱子伦视频在线观看| 欧美xxxx黑人xx丫x性爽| 波野结衣二区三区在线| 午夜老司机福利剧场| 日本成人三级电影网站| 99久久精品国产国产毛片| 91久久精品电影网| 午夜激情欧美在线| 日日啪夜夜撸| 免费看日本二区| 97人妻精品一区二区三区麻豆| 在线免费十八禁| 国产av在哪里看| 亚洲欧美精品综合久久99| 久久精品夜夜夜夜夜久久蜜豆| 亚洲美女搞黄在线观看 | 日本与韩国留学比较| 亚洲欧美精品综合久久99| www.色视频.com| 亚洲在线自拍视频| 欧美日韩中文字幕国产精品一区二区三区| 99精品在免费线老司机午夜| 丰满乱子伦码专区| 亚洲成av人片在线播放无| 免费观看人在逋| 婷婷精品国产亚洲av| 亚洲av免费高清在线观看| 亚洲第一区二区三区不卡| 国产人妻一区二区三区在|