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

    考慮粘彈性人工邊界的高聳進水塔結構地震動態(tài)響應分析

    2016-08-08 00:33:14劉云賀鄭曉東張小剛
    西安理工大學學報 2016年2期
    關鍵詞:進水塔塔體粘彈性

    劉云賀, 鄭曉東, 張小剛

    (西安理工大學 水利水電學院,陜西 西安 710048)

    ?

    考慮粘彈性人工邊界的高聳進水塔結構地震動態(tài)響應分析

    劉云賀, 鄭曉東, 張小剛

    (西安理工大學 水利水電學院,陜西 西安 710048)

    邊界條件的選取對結構的計算設計有著顯著的影響。本文通過ANSYS有限元方法建立了三維進水塔模型,接著分別采用無質量固定邊界和粘彈性人工邊界模擬了高聳進水塔結構在地震作用下的位移、應力和接觸的分布規(guī)律,最后對進水塔抗震安全性進行了計算分析。結果表明:與無質量固定邊界相比,在粘彈性人工邊界條件下,塔體頂部峰值位移減小了10%~30%,塔體峰值應力減小了15%~30%,塔體x正向最大張開值和y正向最大張開值均有所減小,塔體穩(wěn)定安全系數時程最小值增大了20%。研究結果表明粘彈性人工邊界在實際工程設計中具有一定的精確性和適用性。

    高聳進水塔; 粘彈性人工邊界; 動態(tài)響應分析; 接觸分析

    進水塔是引水和泄水系統(tǒng)最前端的水工建筑物,是水利樞紐工程宣泄洪水的安全通道。它的整體抗震穩(wěn)定性關系到整個引水、泄水系統(tǒng)甚至大壩的安危[1],因此,進水塔的整體穩(wěn)定性和抗震性能成為結構設計的重點。在地震激勵作用下,進水塔-地基是一個受力整體,兩者之間的運動和變形相互影響,其中邊界條件的選取對進水塔-地基結構的設計影響比較大。

    目前,國內外在對進水塔地震動力響應及抗震安全性的研究方面取得了一定的成果,趙海濤[2]等采用振型分解反應譜法研究了進水塔的地震動力響應,劉亞琴[3]等采用動力時程法對沙牌水電站進水塔進行了震損模擬,李寧[1]等采用時程分析和反應譜法對紫平鋪進水塔進行了動態(tài)響應研究。上述研究成果大多是以無質量固定邊界為基礎的,在研究過程中忽略了地基輻射阻尼效應,這會使地震波在邊界斷面上發(fā)生反射,影響結構-地基相互作用的應力場和位移場,從而最終影響結構-地基模擬結果的準確性和真實性。

    廖振鵬等[4]研究表明,解決邊界斷面上地震動反射問題的有效辦法是采用人工邊界條件。國內很多學者對人工邊界進行了數值研究,并得出以下結論,結構的尺寸越大,跨度越大,地基輻射阻尼的動力反應影響也就越顯著,例如考慮地基輻射阻尼以后,五嘎沖拱壩[5]動應力結果相應降低了25%~52%,大朝山重力壩[6]主應力結果最大降幅可達50%,雙江口土石壩[7]加速度降低了15%~30%。喻虎圻等[8]對河床式廠房的研究表明,廠房的地震動力響應最大降幅可達64.36%,張青等[9]對長河壩進水塔的分析表明,應力結果降低了15%~20%。

    本文通過ANSYS有限元軟件建立了三維有限元模型,分別采用無質量固定邊界和粘彈性人工邊界模擬了高聳進水塔結構在地震狀態(tài)下的位移、應力和接觸的分布規(guī)律,并將兩種邊界下的計算結果進行了比較。研究結果可為進水塔動力分析中邊界條件的選取提供參考。

    1 粘彈性人工邊界模型

    1994年,Deeks和Randolph提出了粘彈性人工邊界。粘彈性人工邊界條件作為一種應力邊界條件,其應力是邊界節(jié)點速度和位移的函數,表達式為:

    (1)

    對于三維有限元模型[7,10],如圖1所示,粘彈性人工邊界節(jié)點l的法向彈性系數和阻尼系數以及切向彈性系數和阻尼系數分別為:

    圖1 三維有限元模型彈簧阻尼單元示意圖Fig.1 Sketch map of spring damping unit of three dimensional finite element models

    因為粘彈性人工邊界具有較好的魯棒性,人工邊界參數αN、αT在一定范圍內取值,阻尼系數均可以得到較好的計算結果。劉晶波等[10]給出了二維和三維問題中αN和αT的取值范圍及推薦值,對于三維問題,αN取值為4/3,αT取值為2/3。

    圖2為進水塔粘彈性人工邊界地震動輸入模型。

    圖2 進水塔粘彈性人工邊界地震動輸入模型Fig.2 Intake tower viscoelastic artificial boundary seismic input model

    本文通過在人工邊界節(jié)點上施加等效荷載的方式來實現地震動的輸入,即將地震波動的位移和速度時程轉換為等效節(jié)點荷載施加于人工邊界上,完成地震動的輸入。地基的波場包括已知的入射自由場和由基礎上附加結構反射形成的散射場。在截斷邊界處,附加結構反射形成的散射場已被人工邊界吸收掉,因此在人工邊界處僅考慮自由場輸入[8]。

    2 工程實例分析

    2.1工程概況

    以某水電站的塔式進水塔為原型進行建模,進水塔塔頂高程為2 721.00 m,進水口底板高程為2 640.0 m,底板厚5.0 m,地基深度取1倍塔體高度,上下游地基長度均取67.5 m,左右側地基長度均取50 m。進水塔塔體-地基模型如圖3所示。

    圖3 塔體-地基有限元模型Fig.3 The tower body-foundation finite element model

    本文以ANSYS有限元軟件為基礎建立了進水塔-地基有限元模型,接著利用非線性動力時程法對該有限元模型進行了研究分析。進水塔塔體和地基均采用solid45實體單元,單元類型為線性六面體、沙漏控制、縮減積分,網格劃分后單元個數為18 080,結點個數為22 176。附加質量采用mass21質量單元。計算坐標系采用笛卡爾坐標系,坐標原點選在進水塔底板上邊前沿中點處,x軸正方向為順水流方向,y軸正方向垂直于水流方向,z軸正方向垂直向上指向塔頂。求解結構自身的動力特性是結構地震動力分析的第一步,通過模態(tài)分析可以得到結構各階的自振周期和自振頻率[11]。通過有限元軟件分析可知,該進水塔前四階的自振頻率分別為:1.762 6 Hz、1.854 9 Hz、7.160 6 Hz、8.465 1 Hz,進水塔前四階振型如圖4所示。

    圖4 進水塔前四階振型圖Fig.4 Intake tower four vibration modes

    參考《水工建筑物荷載設計規(guī)范》和《水工混凝土結構設計規(guī)范》[12-13],選用的混凝土力學參數如表1所示。泄洪洞進水塔塔座(2663.00高程以下)混凝土強度等級為C30;塔筒(2663.00高程以上)混凝土強度等級為C25?;炷恋谋緲嬆P桶础痘炷两Y構設計規(guī)范(GB50010-2010)》中的混凝土本構關系進行計算,巖體的本構模型采用Mohr-Coulomb塑性模型。塔基巖體力學參數如表2所示。

    表1 混凝土力學參數

    表2 塔基巖體力學參數

    由本工程相關資料和《水工建筑物抗震設計規(guī)范(SL203-97)》可知,該進水塔為1級水工建筑物,工程抗震設防類別為甲類,工程場地類別為Ⅰ類,地震設防烈度為8度,基本烈度為7度。

    本文取基準期100年超越概率為2%的地震動參數作為設計地震,相應的水平地震動峰值加速度為0.304m/s2,特征周期Tg=0.2s;設計反應譜最大值βmax=6.71。根據《建筑抗震設計規(guī)范(GB50011-2010)》的規(guī)定可知,三個方向輸入的地震波峰值比為x:y:z=1∶0.85∶0.65。擬合出三條地震人工波,總時長30 s,時間步長為0.01 s,如圖5所示。

    圖5 設計人工地震波Fig.5 Design of artificial seismic waves

    2.2動態(tài)位移分析

    以塔頂中部關鍵點1為代表,其結點單元編號為6788,如圖6所示。同時輸入3個方向的地震動,即順水流x方向地震波、垂直水流y方向地震波、豎向z方向地震波,設計人工地震動見圖5,得到關鍵點1在三個方向的最大位移和出現時間,如表3所示。

    圖6 塔體頂點示意圖Fig.6 The top of tower schematic diagram

    表3 塔體關鍵點1的相對位移極值及出現時間

    由表3可以得出:與無質量固定邊界相比,在粘彈性人工邊界條件下,塔體頂部峰值位移減小了10%~30%。這是由于地基為彈性體所致。地基的地震動力反應會影響塔體的動力反應,無質量地基均勻輸入情況下地基位移反應較小,而在考慮輻射阻尼效應的粘彈性人工邊界條件下,有質量地基產生較大變形和沉降,位移變化較大。因此,地震動能量向遠域地基逸散的輻射阻尼效應對塔體的動力反應有重要影響,在進行結構動力反應分析時,考慮地基輻射阻尼效應是十分必要的。

    2.3動態(tài)應力分析

    在無質量固定邊界和粘彈性人工邊界條件下,塔體應力分布最不利時刻均為12.80 s,其第一主應力σ1和第三主應力σ3分布如圖7和圖8所示,圖中單位為Pa。

    圖7 無質量固定邊界條件下塔體主應力σ1、σ3云圖Fig.7 The tower principal stresses σ1and σ3contour in the condition of massless foundation

    圖8 粘彈性人工邊界條件下塔體主應力σ1、σ3云圖Fig.8 The tower principal stresses σ1and σ3 contour in the condition of viscoelastic artificial boundaries

    由圖7和圖8可知,在這兩種邊界條件下,塔體第一主應力極值點均未出現在塔體與地基接觸面角緣處,而是在進水口應力集中部位,第三主應力極值點均出現在塔體與地基接觸面角緣處。

    表4為塔體應力極值比較。由表4可以看出,在遭遇基準期為100年超越概率為2%的罕遇地震作用時,兩種邊界條件下,第一主應力極大值均超出混凝土動態(tài)軸心抗拉強度標準值;第三主應力極小值均未超出混凝土動態(tài)軸心抗壓強度標準值。與無質量固定邊界條件相比,在粘彈性邊界條件下,塔體應力極值減小了15%~20%。

    表4 塔體應力極值比較

    2.4接觸分析

    塔體與地基采用面-面接觸分析,地基目標面用Targe170單元來模擬,塔體接觸面用Conta173單元來模擬,因為塔體與地基之間沒有滑動,即采用理想粗糙接觸,塔體與地基交界面的各關鍵點如圖9所示。

    在無質量固定邊界條件下,x正向位移最大時刻(12.92s)和y正向位移最大時刻(7.32s)接觸面的接觸狀態(tài)和張開云圖如圖10和圖11所示,圖11中的單位為m。

    在粘彈性人工邊界條件下,x正向位移最大時刻(12.81s)和y正向位移最大時刻(7.33s)接觸面的接觸狀態(tài)和張開云圖如圖12和圖13所示,圖13中的單位為m。

    圖9 塔體和地基交界面各關鍵點Fig.9 The tower and foundation interface each key points

    圖10 無質量固定邊界條件下塔基接觸面在12.92s和7.32s時刻的接觸狀態(tài)Fig.10 Tower-foundation interface contact status at 12.92s and 7.32s in the condition of massless fixed boundary

    圖11 無質量固定邊界條件下塔基接觸面在12.92s和7.32s時刻的張開值Fig.11 Tower-foundation interface open at 12.92s and 7.32s in the condition of massless fixed boundary

    圖12 粘彈性人工邊界條件下塔基接觸面在12.81s和7.33s時刻的接觸狀態(tài)Fig.12 Tower-foundation interface contact status at 12.81s and 7.33s in the condition of viscoelastic artificial boundaries

    由圖10和圖11可以看出,無質量固定邊界條件下,在x正向位移最大時刻(即12.92 s),塔體與地基間薄弱面的上游側部分被拉裂,下游側部分處于粘結狀態(tài),在關鍵點4附近接觸面張開值最大,為0.622 cm;在y正向位移最大時刻(即7.32 s),塔體與地基間薄弱面的左側部分被拉裂,右側部分處于粘結狀態(tài),靠近關鍵點3處接觸面張開值最大,為0.329 cm。

    由圖12和圖13可以看出,粘彈性人工邊界條件下,在x正向位移最大時刻(即12.81 s),塔體與地基間薄弱面的上游側部分被拉裂,下游側部分處于粘結狀態(tài),靠近關鍵點4處接觸面張開值最大,為0.342 cm;在y正向位移最大時刻(即7.33 s),塔體與地基間薄弱面的左側部分被拉裂,右側部分處于粘結狀態(tài),靠近關鍵點3處接觸面張開值最大,為0.136 cm。

    圖13 粘彈性人工邊界下塔基接觸面在12.81s和7.33s時刻的張開值Fig.13 Tower-foundation interface open values at 12.81s and 7.33s in the condition of viscoelastic artificial boundaries

    綜上可知,無質量固定邊界模型考慮了并不存在的反射地震動,與粘彈性人工邊界模型相比,塔體x正向最大張開值和y正向最大張開值均有所增大,這會引起工程造價估算值偏大,而粘彈性人工邊界則更加合理地反應了地震動的作用機理,能更有效的模擬地震動的遠域能量消散。

    3 進水塔整體穩(wěn)定性評價

    3.1進水塔抗滑穩(wěn)定分析

    根據《水電站進水口設計規(guī)范》,進水塔抗滑穩(wěn)定安全系數用抗剪斷強度計算公式計算:

    (2)

    (3)

    (4)

    將式(3)和(4)代入式(2)計算得到進水塔抗滑穩(wěn)定安全系數如表5所示。

    表5 塔體抗滑穩(wěn)定安全系數極值比較

    由表5可以看出,兩種邊界下,塔體抗滑穩(wěn)定安全系數時程最小值全部大于1,說明進水塔在整個地震過程中沒有滑動的可能,整體抗滑穩(wěn)定性滿足規(guī)范要求,且富裕度較大。與無質量固定邊界條件相比,粘彈性人工邊界條件下,塔體抗滑穩(wěn)定安全系數時程最小值增大了5%~20%,這與分析的塔體位移和應力反應規(guī)律相一致。

    3.2進水塔抗傾覆穩(wěn)定分析

    根據《水電站進水口設計規(guī)范》,進水塔抗傾覆穩(wěn)定安全系數為:

    (5)

    (6)

    (7)

    式中,KS表示抗傾覆安全系數,M0為基礎計算面上全部傾覆力矩之和,M表示基礎計算面上全部抗傾覆力矩之和;zj為單元中心z坐標, Mx、My分別為x、y方向上的全部抗傾覆力矩之和。

    當KS>1.0時,塔體沒有傾倒的趨勢,至少在這一時刻是穩(wěn)定的。

    當KS<1.0時,塔體有了傾倒的趨勢,此時塔體是不穩(wěn)定的。

    將式(6)和(7)代入式(5),計算得到進水塔抗傾覆穩(wěn)定安全系數如表6所示。

    表6 塔體抗傾覆穩(wěn)定安全系數極值比較

    由表6可以看出,兩種邊界條件下塔體抗傾覆穩(wěn)定安全系數時程最小值全部大于1,說明進水塔在整個地震過程中沒有傾覆的趨勢,整體抗傾覆穩(wěn)定性滿足規(guī)范要求,且富裕度較大。與無質量固定邊界條件相比,粘彈性人工邊界條件下,塔體抗傾覆穩(wěn)定安全系數時程最小值增大了20%,這與分析的塔體位移和應力反應規(guī)律相一致。

    4 結 論

    本文通過ANSYS有限元軟件建立了三維進水塔模型,并分別采用無質量固定邊界和粘彈性人工邊界,模擬了地震動作用下進水塔結構位移、應力和接觸的分布規(guī)律。通過分析可以得到如下結論:

    1) 無質量固定邊界均勻輸入方法的計算結果是偏保守的;而粘彈性人工邊界則更真實地反應了地震動的作用機理,能更合理的模擬地震動的遠域能量消散,因此,在進行塔體結構動力反應分析時,考慮地基輻射阻尼效應是十分必要的。

    2) 在遭遇基準期為100年超越概率為2%的罕遇地震作用時,進水塔在兩種邊界條件下的抗滑和抗傾覆穩(wěn)定性均滿足規(guī)范要求且富裕度較大。與無質量固定邊界條件相比,粘彈性人工邊界條件下,塔體穩(wěn)定安全系數時程最小值增大了20%,與塔體位移和應力反應規(guī)律相一致。

    [1]李寧, 李騫, 任堂, 等. 汶川地震紫坪鋪進水塔的動態(tài)響應研究[J]. 地下空間與工程學報, 2014, 10(5):1127-1134.

    LI Ning, LI Qian, REN Tang, et al. Dynamic response analysis of Zipingpu intake tower under Wenchuan earthquake action[J]. Chinese Journal of Underground Space and Engineering, 2014, 10(5):1127-1134.

    [2]趙海濤, 羅勇軍,王潘繡,等. 高聳鋼筋混凝土進水塔結構抗震穩(wěn)定性分析與安全評估[J]. 水利與建筑工程學報, 2012, 10(6):106-109.

    ZHAO Haitao, LUO Yongjun, WANG Panxiu, et al. Seismic stability analysis and safety assessment for reinforced concrete intake tower structure [J]. Journal of Water Resources and Architectural Engineering, 2012, 10(6):106-109.

    [3]劉亞琴, 趙蘭浩, 錢文江.沙牌水電站進水塔震損模擬分析[J]. 水利水電技術, 2015, 46(1):30-33.

    LIU Yaqin, ZHAO Lanhao, Qian Wenjiang.Simulation analysis on seismic damage of intake tower of Shapai Hydropower Station[J].Water Resources and Hydropower Engineering, 2015, 46(1):30-33.

    [4]廖振鵬. 工程波動理論導論[M]. 北京:科學出版社, 2002.

    [5]陳平,何蘊龍,陳海霞.地基邊界對五嘎沖拱壩地震動力響應的影響[J].武漢大學學報(工學版),2011,44(5):554-559.

    CHEN Ping, HE Yunlong, CHEN Haixia. Effect of foundation boundary on earthquake dynamic response of Wugachong arch dam[J].Engineering Journal of Wuhan University, 2011,44(5):554-559.

    [6]何建濤,馬懷發(fā),張伯艷,等.黏彈性人工邊界地震動輸入方法及實現[J].水利學報, 2010, 41(8):960-969.

    HE Jiantao, MA Huaifa, ZHANG Boyan, et al. Method and realization of seismic motion input of viscous-spring boundary[J]. Journal of Hydraulic Engineering, 2010, 41(8):960-969.

    [7]李宏恩, 李同春, 田景元, 等.黏-彈性人工邊界在雙江口土石壩動力分析中的應用[J].巖土力學,2008,29(增1):189-192.

    LI Hongen, LI Tongchun, TIAN Jingyuan, et al. Application of viscoelastic artificial boundary to dynamic response analysis of Shuangjiangkou rock-fill dam[J]. Rock and Soil Mechanics, 2008,29(S1):189-192.

    [8]喻虎圻, 何蘊龍, 曹學興, 等. 基于粘彈性邊界的河床式廠房地震動力響應分析[J]. 武漢大學學報(工學版), 2015, 48(1): 27-33.

    YU Huqi, HE Yunlong, CAO Xuexing, et al. Seismic analysis of powerhouse in river channel based on viscoelastic boundary[J]. Engineering Journal of Wuhan University, 2015, 48(1): 27-33.

    [9]張青, 徐元杰, 夏明錟, 等. 長河壩進水口粘彈性人工邊界有限元分析[J]. 中國水運, 2015, 15(3): 180-183.

    ZHANG Qing, XU Yuanjie, XIA Mingtan, et al. The intake of Changhe dam viscoelastic artificial boundary finite element analysis[J]. China Water Transport, 2015, 15(3): 180-183.

    [10]劉晶波, 杜義欣, 閆秋實. 粘彈性人工邊界及地震動輸入在通用有限元軟件中的實現[J]. 防災減災工程學報, 2007, 27(增1):37-42.

    LIU Jingbo, DU Yixin, YAN Qiushi. Viscoelastic artificial boundary and seismic input in the general finite element software realization[J]. Journal of Disaster Prevention and Mitigation Engineering, 2007, 27(S1):37-42.

    [11]陶磊, 張俊發(fā). 遠場大震對深厚覆蓋土層地基上超高層建筑物的影響分析[J]. 西安理工大學學報, 2014, 30(1):79-84.

    TAO Lei, ZHANG Junfa. Influence and analysis of suoer-rise building on soil layer with deep deposit subjected to strong far-field ground motion[J]. Journal of Xi’an University of Technology, 2014, 30(1):79-84.

    [12]中華人民共和國電力工業(yè)部. 水工建筑物荷載設計規(guī)范(DL/5077-1997)[S]. 北京:中國電力出版社, 1997.

    [13]中華人民共和國水利部. 水工混凝土結構設計規(guī)范(SL/191-2008)[S]. 北京:中國水利水電出版社, 2005.

    (責任編輯周蓓)

    Considering the analysis of high intake tower seismic dynamic response at the viscoelastic artificial boundary

    LIU Yunhe, ZHENG Xiaodong, ZHANG Xiaogang

    (School of Water Resources and Hydro-electric Engineering,Xi’an University of Technology,Xi’an 710048, China)

    The selection of boundary conditions has a significant impact on the design of the structure. In this paper, three-dimensional finite element models are established by ANSYS. This paper uses the massless fixed boundary and viscoelastic artificial boundary to simulate high intake tower structure of distribution of displacement, stress and contact in the earthquake state. Finally, the seismic safety of the intake tower is calculated and analyzed. Compared with the massless fixed boundary, in the viscoelastic artificial boundary the tower top peak displacement is reduced by 10% to 30%, the tower peak stress is reduced by 15% to 30%, the towerxpositive direction andypositive direction of the largest open values are reduced, and the minimum value of tower body stability safety coefficient is increased by 20%. Results show that the viscoelastic artificial boundary has certain precision and applicability in practical engineering design.

    high intake tower; the viscoelastic artificial boundary; dynamic response analysis; contact analysis

    1006-4710(2016)02-0134-08

    10.19322/j.cnki.issn.1006-4710.2016.02.002

    2015-10-26

    國家自然科學基金資助項目(51179154/E090801)

    劉云賀,男,教授,博導,博士,研究方向為水工結構抗震。E-mail: liuyunhe1968@163.com

    TV312

    A

    猜你喜歡
    進水塔塔體粘彈性
    外載荷作用下塔器開孔補強計算探討
    二維粘彈性棒和板問題ADI有限差分法
    含啟閉機房進水塔結構靜動力分析
    某水電站進水塔靜力穩(wěn)定性分析
    Preliminary Design of a Submerged Support Structure for Floating Wind Turbines
    復雜施工環(huán)境下的進水塔安裝
    港工技術(2020年5期)2020-10-22 06:39:06
    基于振型分解反應譜法的深孔泄洪洞進水塔動力分析
    水電站設計(2020年4期)2020-07-16 08:23:40
    時變時滯粘彈性板方程的整體吸引子
    應用數學(2020年2期)2020-06-24 06:02:30
    冷卻塔爆破拆除傾倒解體及振動研究
    爆炸與沖擊(2019年2期)2019-02-27 02:25:00
    塔體現場改造技術
    极品教师在线免费播放| 麻豆久久精品国产亚洲av| 欧美中文综合在线视频| 欧美激情高清一区二区三区| 亚洲一区二区三区色噜噜| 欧美另类亚洲清纯唯美| 男女下面插进去视频免费观看| 村上凉子中文字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 老鸭窝网址在线观看| 亚洲视频免费观看视频| 中文亚洲av片在线观看爽| 美女扒开内裤让男人捅视频| 国产高清有码在线观看视频 | 国产高清视频在线播放一区| 一卡2卡三卡四卡精品乱码亚洲| 日本a在线网址| 男人舔女人的私密视频| 久久青草综合色| 精品国产一区二区久久| 女性生殖器流出的白浆| 久久国产精品影院| 精品久久久久久成人av| 精品久久久久久久毛片微露脸| 韩国精品一区二区三区| 黄色丝袜av网址大全| 啪啪无遮挡十八禁网站| 天天躁狠狠躁夜夜躁狠狠躁| 欧美性长视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 香蕉丝袜av| 久久狼人影院| 午夜激情av网站| 亚洲情色 制服丝袜| 午夜老司机福利片| 亚洲黑人精品在线| 他把我摸到了高潮在线观看| 天天躁夜夜躁狠狠躁躁| 黑人欧美特级aaaaaa片| 香蕉国产在线看| 亚洲中文av在线| av天堂在线播放| 成年版毛片免费区| 亚洲一区高清亚洲精品| 日本免费a在线| 久久中文看片网| 九色亚洲精品在线播放| 50天的宝宝边吃奶边哭怎么回事| 如日韩欧美国产精品一区二区三区| 美女高潮喷水抽搐中文字幕| 此物有八面人人有两片| 欧美绝顶高潮抽搐喷水| 国产亚洲欧美精品永久| 欧美成人性av电影在线观看| 欧美黑人欧美精品刺激| 国产私拍福利视频在线观看| 99久久久亚洲精品蜜臀av| 亚洲国产欧美网| 免费不卡黄色视频| 国产成人精品无人区| 亚洲av成人av| 亚洲成人精品中文字幕电影| 国产视频一区二区在线看| 亚洲av片天天在线观看| svipshipincom国产片| 亚洲人成电影免费在线| 日本免费a在线| 一二三四社区在线视频社区8| 亚洲av美国av| 亚洲五月色婷婷综合| 亚洲欧美日韩高清在线视频| 狂野欧美激情性xxxx| 两性夫妻黄色片| 国产av精品麻豆| 中文字幕人妻丝袜一区二区| 波多野结衣高清无吗| 国产亚洲av高清不卡| 亚洲国产欧美网| 熟女少妇亚洲综合色aaa.| 国产精华一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 男女床上黄色一级片免费看| 最好的美女福利视频网| 精品不卡国产一区二区三区| 99国产精品一区二区三区| 91大片在线观看| 久久精品国产99精品国产亚洲性色 | 精品第一国产精品| 欧美色欧美亚洲另类二区 | 又黄又粗又硬又大视频| 高潮久久久久久久久久久不卡| 自拍欧美九色日韩亚洲蝌蚪91| 久久香蕉激情| 夜夜爽天天搞| 叶爱在线成人免费视频播放| 国产一区二区在线av高清观看| 美女午夜性视频免费| ponron亚洲| 999精品在线视频| 黑人巨大精品欧美一区二区mp4| 好看av亚洲va欧美ⅴa在| 国产激情欧美一区二区| √禁漫天堂资源中文www| 久久精品国产99精品国产亚洲性色 | 国产av一区二区精品久久| 此物有八面人人有两片| 老司机靠b影院| 黄网站色视频无遮挡免费观看| 精品福利观看| 91在线观看av| 国产在线精品亚洲第一网站| 操美女的视频在线观看| 久久精品91无色码中文字幕| 欧美激情极品国产一区二区三区| avwww免费| 黄色 视频免费看| 日韩国内少妇激情av| 50天的宝宝边吃奶边哭怎么回事| 亚洲第一欧美日韩一区二区三区| 久9热在线精品视频| 国产精品 国内视频| 91麻豆av在线| 国产欧美日韩一区二区精品| 国产精品久久久久久人妻精品电影| 免费高清视频大片| 午夜影院日韩av| 999精品在线视频| www日本在线高清视频| 久久午夜综合久久蜜桃| 亚洲人成77777在线视频| 亚洲av五月六月丁香网| 一级作爱视频免费观看| 人人妻人人爽人人添夜夜欢视频| 又黄又粗又硬又大视频| 国产国语露脸激情在线看| 欧美黑人精品巨大| 亚洲av成人一区二区三| 一本久久中文字幕| 国产又爽黄色视频| 黄色毛片三级朝国网站| 丝袜人妻中文字幕| 成人三级黄色视频| 亚洲中文日韩欧美视频| 桃红色精品国产亚洲av| 欧美黑人精品巨大| 精品少妇一区二区三区视频日本电影| 一本大道久久a久久精品| 一区二区三区国产精品乱码| 欧美久久黑人一区二区| 高清在线国产一区| 精品欧美一区二区三区在线| 国产在线精品亚洲第一网站| 9色porny在线观看| 亚洲avbb在线观看| 婷婷丁香在线五月| 久久久久久亚洲精品国产蜜桃av| 99在线视频只有这里精品首页| 啪啪无遮挡十八禁网站| 午夜亚洲福利在线播放| 免费看美女性在线毛片视频| 欧美激情极品国产一区二区三区| 一进一出好大好爽视频| av在线播放免费不卡| 黄片小视频在线播放| 国产又色又爽无遮挡免费看| 亚洲第一av免费看| 久久婷婷成人综合色麻豆| 亚洲一区中文字幕在线| 亚洲专区字幕在线| 一本大道久久a久久精品| 亚洲专区字幕在线| 韩国av一区二区三区四区| 亚洲 欧美一区二区三区| 中文字幕久久专区| 免费观看精品视频网站| 久久热在线av| 欧美大码av| 亚洲第一电影网av| √禁漫天堂资源中文www| 精品欧美一区二区三区在线| 女人被狂操c到高潮| av片东京热男人的天堂| 精品无人区乱码1区二区| 欧美日韩福利视频一区二区| 母亲3免费完整高清在线观看| 亚洲九九香蕉| 亚洲性夜色夜夜综合| 欧美日本中文国产一区发布| 99在线人妻在线中文字幕| 欧美老熟妇乱子伦牲交| 亚洲精品久久成人aⅴ小说| 婷婷精品国产亚洲av在线| 侵犯人妻中文字幕一二三四区| 可以在线观看的亚洲视频| 激情视频va一区二区三区| а√天堂www在线а√下载| 狠狠狠狠99中文字幕| 欧美黑人精品巨大| 身体一侧抽搐| 久久久久久免费高清国产稀缺| 亚洲人成网站在线播放欧美日韩| 精品少妇一区二区三区视频日本电影| 亚洲av成人一区二区三| 91九色精品人成在线观看| 国产成人免费无遮挡视频| 黄色a级毛片大全视频| 亚洲五月婷婷丁香| 欧美激情 高清一区二区三区| 手机成人av网站| 操美女的视频在线观看| 亚洲无线在线观看| 亚洲电影在线观看av| 9色porny在线观看| 视频在线观看一区二区三区| 日韩精品免费视频一区二区三区| 丝袜美足系列| 午夜老司机福利片| 麻豆av在线久日| 午夜福利视频1000在线观看 | 好看av亚洲va欧美ⅴa在| 成人国产综合亚洲| 精品国产亚洲在线| 一进一出抽搐gif免费好疼| 午夜影院日韩av| 国产不卡一卡二| 日本精品一区二区三区蜜桃| 国产一区二区在线av高清观看| 香蕉丝袜av| 亚洲精品久久国产高清桃花| 欧美日韩一级在线毛片| 日本撒尿小便嘘嘘汇集6| 色av中文字幕| 久久人妻av系列| 麻豆一二三区av精品| 欧美绝顶高潮抽搐喷水| 黄色成人免费大全| 一区二区三区高清视频在线| 国产av在哪里看| 久久精品国产亚洲av香蕉五月| 国产精品,欧美在线| 久久久国产精品麻豆| 免费人成视频x8x8入口观看| 禁无遮挡网站| 99久久精品国产亚洲精品| 91老司机精品| 欧美激情久久久久久爽电影 | 久久久久久久午夜电影| 国产伦人伦偷精品视频| 一本久久中文字幕| aaaaa片日本免费| 丰满人妻熟妇乱又伦精品不卡| 国产精品野战在线观看| 熟妇人妻久久中文字幕3abv| 激情视频va一区二区三区| 无遮挡黄片免费观看| 午夜免费观看网址| 黄色视频,在线免费观看| 又大又爽又粗| 非洲黑人性xxxx精品又粗又长| 婷婷丁香在线五月| 日韩大尺度精品在线看网址 | 长腿黑丝高跟| 男女之事视频高清在线观看| 免费看十八禁软件| 欧美色欧美亚洲另类二区 | 日韩免费av在线播放| 午夜a级毛片| 亚洲人成77777在线视频| 久久天躁狠狠躁夜夜2o2o| 热99re8久久精品国产| 成熟少妇高潮喷水视频| 一级毛片精品| 9热在线视频观看99| 日韩免费av在线播放| 欧美在线黄色| 亚洲中文av在线| 久久天躁狠狠躁夜夜2o2o| 国产成人一区二区三区免费视频网站| 欧美中文综合在线视频| 国产成人精品久久二区二区免费| 亚洲专区国产一区二区| 神马国产精品三级电影在线观看 | www.www免费av| 亚洲成a人片在线一区二区| 欧美一级a爱片免费观看看 | 国产91精品成人一区二区三区| av在线天堂中文字幕| 这个男人来自地球电影免费观看| 午夜老司机福利片| 亚洲精品国产色婷婷电影| 亚洲自偷自拍图片 自拍| 国产一区二区三区视频了| 搡老岳熟女国产| 久久精品人人爽人人爽视色| 一边摸一边抽搐一进一小说| 久久久久久亚洲精品国产蜜桃av| 免费观看精品视频网站| 亚洲色图av天堂| 亚洲全国av大片| 亚洲视频免费观看视频| 免费观看人在逋| 免费在线观看黄色视频的| 多毛熟女@视频| 91在线观看av| 国产99白浆流出| av视频免费观看在线观看| 国产伦人伦偷精品视频| 可以在线观看的亚洲视频| 久久久久精品国产欧美久久久| 国产精品亚洲美女久久久| 精品国产亚洲在线| 欧美激情 高清一区二区三区| 国产精品久久久久久精品电影 | 久久国产乱子伦精品免费另类| 国产91精品成人一区二区三区| 久久中文字幕人妻熟女| 国产精品久久视频播放| 性少妇av在线| 久久热在线av| 97人妻天天添夜夜摸| 亚洲精品在线观看二区| 精品乱码久久久久久99久播| 男女之事视频高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久精品国产欧美久久久| 香蕉久久夜色| 精品电影一区二区在线| 国产免费av片在线观看野外av| av在线天堂中文字幕| 一级作爱视频免费观看| 天堂影院成人在线观看| 欧美乱码精品一区二区三区| 99久久久亚洲精品蜜臀av| 国产精品久久久人人做人人爽| 日韩欧美免费精品| 99香蕉大伊视频| 国内久久婷婷六月综合欲色啪| 日韩国内少妇激情av| 美女大奶头视频| 淫秽高清视频在线观看| 欧美日韩一级在线毛片| 精品第一国产精品| 一二三四在线观看免费中文在| 无遮挡黄片免费观看| 国产精品av久久久久免费| 人人妻人人澡人人看| 亚洲色图综合在线观看| 国产一区二区三区视频了| 日韩欧美一区二区三区在线观看| 免费人成视频x8x8入口观看| 中文字幕色久视频| 亚洲一区二区三区不卡视频| 国产精品影院久久| 国产精品爽爽va在线观看网站 | 在线观看免费视频网站a站| av在线天堂中文字幕| 美女 人体艺术 gogo| 黄色a级毛片大全视频| 激情视频va一区二区三区| 在线免费观看的www视频| 99riav亚洲国产免费| svipshipincom国产片| 美女国产高潮福利片在线看| 国产又爽黄色视频| 欧美激情高清一区二区三区| 两人在一起打扑克的视频| 欧美最黄视频在线播放免费| 亚洲国产高清在线一区二区三 | 午夜日韩欧美国产| 自线自在国产av| 国产精品精品国产色婷婷| 一区在线观看完整版| 国产成人av激情在线播放| 精品国产超薄肉色丝袜足j| 欧美激情 高清一区二区三区| 免费搜索国产男女视频| 老司机午夜福利在线观看视频| 超碰成人久久| www.自偷自拍.com| 色综合站精品国产| 亚洲五月婷婷丁香| 午夜免费鲁丝| 激情视频va一区二区三区| 久久精品国产清高在天天线| 在线观看午夜福利视频| 国产av一区在线观看免费| 激情在线观看视频在线高清| 久久国产乱子伦精品免费另类| av在线天堂中文字幕| 97超级碰碰碰精品色视频在线观看| 久久 成人 亚洲| 亚洲中文字幕一区二区三区有码在线看 | 国产精品美女特级片免费视频播放器 | 中文字幕人妻丝袜一区二区| or卡值多少钱| 操美女的视频在线观看| 亚洲国产欧美网| 色在线成人网| 高清黄色对白视频在线免费看| 欧美日韩亚洲国产一区二区在线观看| 久久久久国产精品人妻aⅴ院| 97人妻精品一区二区三区麻豆 | 黄色 视频免费看| 嫩草影院精品99| 久久久国产成人精品二区| 少妇熟女aⅴ在线视频| 他把我摸到了高潮在线观看| 午夜福利视频1000在线观看 | 巨乳人妻的诱惑在线观看| 高清黄色对白视频在线免费看| 亚洲精品国产一区二区精华液| 精品福利观看| 久久中文看片网| 给我免费播放毛片高清在线观看| 欧美最黄视频在线播放免费| 一区二区三区高清视频在线| 91成年电影在线观看| 久久久久久亚洲精品国产蜜桃av| 国产熟女xx| 午夜a级毛片| 禁无遮挡网站| 中文字幕av电影在线播放| av天堂在线播放| 两个人视频免费观看高清| 精品久久久久久成人av| 国产激情欧美一区二区| 成人av一区二区三区在线看| 露出奶头的视频| 精品久久蜜臀av无| 老熟妇乱子伦视频在线观看| 男女做爰动态图高潮gif福利片 | 性欧美人与动物交配| 国产精品98久久久久久宅男小说| 亚洲成人国产一区在线观看| 欧美日韩一级在线毛片| 精品国产国语对白av| 日本欧美视频一区| 一边摸一边抽搐一进一小说| 久久久久久久午夜电影| 久久狼人影院| 国产精品爽爽va在线观看网站 | 男男h啪啪无遮挡| 欧美精品亚洲一区二区| av免费在线观看网站| 村上凉子中文字幕在线| 中文字幕人妻丝袜一区二区| 很黄的视频免费| 777久久人妻少妇嫩草av网站| 午夜免费观看网址| 99香蕉大伊视频| 91大片在线观看| 亚洲色图av天堂| 国产精品久久久人人做人人爽| x7x7x7水蜜桃| 欧美绝顶高潮抽搐喷水| 日韩欧美三级三区| 精品一区二区三区四区五区乱码| 国产精品国产高清国产av| 老汉色av国产亚洲站长工具| 性少妇av在线| 精品不卡国产一区二区三区| 日韩 欧美 亚洲 中文字幕| 成年人黄色毛片网站| 中文字幕精品免费在线观看视频| 国产av一区二区精品久久| 亚洲欧美一区二区三区黑人| 国产aⅴ精品一区二区三区波| 国产精品久久久久久亚洲av鲁大| 国产精品 国内视频| 日日摸夜夜添夜夜添小说| 国产不卡一卡二| 欧美在线一区亚洲| 黄片播放在线免费| 午夜视频精品福利| 性色av乱码一区二区三区2| 久久伊人香网站| 老熟妇仑乱视频hdxx| 精品免费久久久久久久清纯| 国产片内射在线| 97碰自拍视频| 亚洲精华国产精华精| 亚洲色图综合在线观看| 日本欧美视频一区| 波多野结衣一区麻豆| 在线十欧美十亚洲十日本专区| 国产免费av片在线观看野外av| 久99久视频精品免费| 亚洲色图av天堂| 亚洲中文av在线| 午夜福利成人在线免费观看| 久久婷婷成人综合色麻豆| 国产高清有码在线观看视频 | 亚洲一区二区三区不卡视频| 1024香蕉在线观看| 女人被狂操c到高潮| 老司机靠b影院| 免费搜索国产男女视频| 精品久久久久久久人妻蜜臀av | videosex国产| 国产高清videossex| 午夜影院日韩av| 成年人黄色毛片网站| 99国产精品一区二区三区| 亚洲av成人av| 国产成人啪精品午夜网站| 男人舔女人的私密视频| av超薄肉色丝袜交足视频| 波多野结衣巨乳人妻| avwww免费| 日本vs欧美在线观看视频| 激情视频va一区二区三区| 制服丝袜大香蕉在线| 免费观看精品视频网站| 一a级毛片在线观看| 日韩一卡2卡3卡4卡2021年| 国产视频一区二区在线看| 亚洲av电影在线进入| 黄色丝袜av网址大全| 最好的美女福利视频网| 一区在线观看完整版| 精品福利观看| 少妇裸体淫交视频免费看高清 | 国产乱人伦免费视频| 黑丝袜美女国产一区| 亚洲精品国产色婷婷电影| 亚洲国产精品999在线| 国产欧美日韩一区二区精品| 欧美精品亚洲一区二区| 香蕉丝袜av| 丝袜美足系列| 在线免费观看的www视频| 国产精品久久电影中文字幕| av中文乱码字幕在线| 久9热在线精品视频| 一级a爱视频在线免费观看| 极品教师在线免费播放| 窝窝影院91人妻| 99精品欧美一区二区三区四区| 精品久久久久久成人av| 久久人人精品亚洲av| 波多野结衣一区麻豆| 欧美成人午夜精品| 亚洲av第一区精品v没综合| 成人永久免费在线观看视频| 午夜免费观看网址| 国产精品久久电影中文字幕| 精品久久久久久,| 欧美成人一区二区免费高清观看 | 97人妻精品一区二区三区麻豆 | 18禁裸乳无遮挡免费网站照片 | 亚洲精品av麻豆狂野| 女生性感内裤真人,穿戴方法视频| 久久青草综合色| 亚洲熟妇熟女久久| 免费一级毛片在线播放高清视频 | 久久 成人 亚洲| 久久久久久大精品| 女人爽到高潮嗷嗷叫在线视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲自偷自拍图片 自拍| 黄频高清免费视频| 欧美国产精品va在线观看不卡| 久久久久久久精品吃奶| 国语自产精品视频在线第100页| 中国美女看黄片| 精品国产超薄肉色丝袜足j| 久久久久国内视频| 亚洲一区二区三区色噜噜| 国产麻豆成人av免费视频| 午夜影院日韩av| 国产高清videossex| 黑丝袜美女国产一区| 中文亚洲av片在线观看爽| 欧美黄色淫秽网站| 亚洲av美国av| 国产精品乱码一区二三区的特点 | 91字幕亚洲| 久久伊人香网站| 午夜福利免费观看在线| 亚洲专区国产一区二区| e午夜精品久久久久久久| 深夜精品福利| 久久久久国产精品人妻aⅴ院| 欧美日韩一级在线毛片| 99re在线观看精品视频| 国产精品电影一区二区三区| av超薄肉色丝袜交足视频| 国产av在哪里看| 91老司机精品| 大型av网站在线播放| 自线自在国产av| 中文字幕最新亚洲高清| 大陆偷拍与自拍| 亚洲中文字幕一区二区三区有码在线看 | 欧美日本视频| 久久精品国产亚洲av高清一级| 777久久人妻少妇嫩草av网站| 中文亚洲av片在线观看爽| 中文字幕另类日韩欧美亚洲嫩草| 夜夜躁狠狠躁天天躁| 午夜福利,免费看| 亚洲情色 制服丝袜| 欧美日本中文国产一区发布| 亚洲精品av麻豆狂野| 国产av精品麻豆| 曰老女人黄片| 日韩大码丰满熟妇| 午夜视频精品福利| 亚洲一区高清亚洲精品| 久久国产精品男人的天堂亚洲| 国产av一区在线观看免费| 欧美日本中文国产一区发布| 欧美乱码精品一区二区三区| 亚洲七黄色美女视频| 久久香蕉精品热|