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

    基于OpenSEES的RC柱擬靜力實驗數(shù)值分析

    2016-11-03 05:11:51孟麗巖
    黑龍江科技大學學報 2016年1期
    關鍵詞:中柱梁柱個數(shù)

    王 濤, 孫 嚴, 孟麗巖

    (黑龍江科技大學 建筑工程學院, 哈爾濱 150022)

    ?

    基于OpenSEES的RC柱擬靜力實驗數(shù)值分析

    王濤,孫嚴,孟麗巖

    (黑龍江科技大學 建筑工程學院, 哈爾濱 150022)

    采用OpenSEES軟件中基于位移法梁柱單元建模時,需要人為事先定義單元模型參數(shù)。為了驗證參數(shù)取值對數(shù)值模擬結果精度的影響,基于鋼筋混凝土中柱C的擬靜力實驗數(shù)據(jù),進行建模參數(shù)敏感性分析。結果表明:基于位移法梁柱單元建模時,積分點個數(shù)及混凝土纖維劃分數(shù)量對模擬結果的影響并不十分敏感;與纖維截面纖維個數(shù)、積分點個數(shù)相比,單元個數(shù)對模擬精度的影響更大。單元劃分數(shù)量、單元個數(shù)為5、積分點個數(shù)在5和7之間可以提高模擬精度。

    擬靜力實驗; OpenSEES; 位移法; 梁柱單元; 滯回曲線

    0 引 言

    鋼筋混凝土框架結構是應用較廣泛的結構形式。國內外眾多學者對鋼筋混凝土框架結構的抗震性能和地震反應分析方面做了大量深入研究。2002年,陳滔等對國內外鋼筋混凝土框架非彈性地震反應分析模型的研究進展進行了介紹和評述[1]。2003年,羅文斌等對鋼筋混凝土框架彈塑性位移的解構規(guī)則與構件的目標側移角進行了研究[2]。2010年,陸新征等[3]開發(fā)了可以模擬復雜結構倒塌的程序。目前,常用的結構非線性分析軟件種類繁多,如基于宏觀有限元的OpenSEES、CANNY等,以及基于微觀有限元的ABAQUS、ANSYS、MARC等。在軟件應用時,通常需要事先定義結構的數(shù)值模型、單元模型、單元或材料的恢復力模型以及數(shù)值求解方法等。

    研究者們通常利用實驗數(shù)據(jù)來驗證所采用數(shù)值模擬方法的有效性。即在已知實驗數(shù)據(jù)的情況下,通過調整數(shù)值模型或模型參數(shù),可以得到相對較好的模擬結果。然而,數(shù)值模擬更重要的任務是用來預測及評估結構的抗震性能及地震反應,即在未知實驗結果的前提下來預測結構的響應。這就對數(shù)值模擬方法的有效性提出了更大的挑戰(zhàn)。

    為了檢驗現(xiàn)有的結構數(shù)值模擬方法的可靠性,近幾年,國際及國內陸續(xù)舉辦了多種類型的模型結構盲測實驗,并公開了相應的實驗數(shù)據(jù)及比賽結果[4-6]。從比賽的結果可以看出,對于同一個結構實驗,不同的參賽者得到的模擬結果有時差異較大,即使采用相同軟件也會有不同。可見,采用目前的數(shù)值模擬來再現(xiàn)結構震害準確性仍存在一定的問題。

    為探究這些人為定義參數(shù)對模擬結果的影響,筆者采用清華大學鋼筋混凝土框架結構柱構件實驗的中柱C的實驗數(shù)據(jù),分析纖維截面纖維個數(shù)、單元個數(shù)、積分點個數(shù)對模擬結果的影響規(guī)律,給出參數(shù)設置建議。

    1 基于位移法的梁柱單元

    基于位移法的梁柱單元的任意截面在受分布力的作用下,根據(jù)力的分布情況,相應截面進入塑性狀態(tài),能更真實的反應實際構件的受力情況。基于位移法的梁柱單元,高斯積分點個數(shù)一般為5個,采用三階埃爾米特差值函數(shù)構建基于結點位移的形函數(shù),通過虛功原理求剛度矩陣。

    基于位移法的梁柱單元在定義時,首先沿構件柱豎向等間距劃分多個節(jié)點,兩個節(jié)點之間定義為一個單元,單元中設定不同積分點,積分點處設置纖維截面。纖維截面將構件截面劃分成很多小纖維(包括鋼筋纖維和混凝土纖維)對每一根纖維只考慮它的軸向本構關系,纖維可以定義不同的本構關系。然后通過纖維截面確定截面的力-變形關系,用數(shù)值積分方法確定整個單元的力-變形關系,可以同時反映彎曲變形和軸向變形的耦合,截面布置如圖1所示。

    圖1 單元截面布置

    2 鋼筋混凝土柱的低周往復加載

    2.1試件材料與配筋

    2011年,清華大學進行鋼筋混凝土框架的擬靜力倒塌實驗,由于影響結構倒塌的因素很多,所以并不能分析數(shù)值模型的誤差來源[6]。因此,清華大學對底層框架柱和關鍵節(jié)點[7]進行了擬靜力倒塌實驗,得出關鍵柱對結構整體變形起更加重要的作用。實驗所用柱的配筋圖及實驗數(shù)據(jù)公布于抗倒塌協(xié)會官方網(wǎng)站。

    柱采用的混凝土強度等級為C30,混凝土立方體抗壓強度實測值fcu=30.1 MPa??v筋采用HRB335,縱向受力鋼筋沿截面四周均勻對稱布置,箍筋采用HPB235,截面尺寸200 mm×200 mm。柱配筋見圖2。鋼筋力學性能見表1,表中強度和彈性模量仍采用標稱直徑計算[8]。

    圖2 構件柱配筋

    d/mmfy/MPafu/MPaεyE/MPaδ/%64415290.002220394134.2085828550.002028985028.80104817450.002026543323.60

    2.2實驗加載與結果

    框架中柱C豎向通過柱頂端軸壓千斤頂加載,豎向軸力分別為256和276 kN,在實驗過程中保持恒定。水平方向通過水平推拉千斤頂首先按力循環(huán)加載,再按位移控制加載。水平加載點位置在柱頂下方,中柱C的受力如圖3所示。

    圖3 中柱C受力示意

    水平力F循環(huán)加載依次為10、20、30、40 kN,各級荷載循環(huán)一圈。位移控制加載制度,位移s大小依次為10.0、15.0、20.0、25.0、30.0、37.5、45.0、55.0 mm的水平位移控制加載,各級荷載循環(huán)兩圈。

    中柱C實驗滯回曲線如圖4所示。構件柱的滯回曲線形狀比較飽滿,該構件塑性變形能力強,具有很強的抗震性能和耗能能力。滯回曲線形狀呈梭形,說明中柱C受P-△效應的影響。

    圖4 實驗滯回曲線

    3 鋼筋混凝土柱的數(shù)值模擬

    3.1混凝土纖維個數(shù)影響

    矩形纖維截面需要將截面劃分為1~4區(qū)和核心區(qū),如圖5所示。1~4區(qū)統(tǒng)稱為保護層區(qū),在保護層區(qū)、核心區(qū)劃分不同纖維,每一個纖維用來定義混凝土單軸本構。為了分析纖維數(shù)量對數(shù)值模擬結果精度的影響,文中針對中柱C考慮了六種不同的截面纖維劃分方案,如表2所示。

    圖5 纖維截面區(qū)域劃分

    方案核心區(qū)1區(qū)2區(qū)3區(qū)4區(qū)12020204424040408838080801616410010010020205200200200252563003003005050

    在六次的數(shù)值模擬中,柱均采用基于位移法的梁柱單元,單元個數(shù)取5,積分點個數(shù)取6,方案1~方案6的模擬結果與實驗結果對比分別如圖6a至6f所示。圖6a~6f,各圖中虛線為實驗值,實線為模擬值。

    采用地震模擬振動臺波形再現(xiàn)的定量判別方法[9],通過Matlab求解模擬結果與實驗結果的相關系數(shù),評價模擬滯回曲線和實驗滯回曲線的相近程度。相關系數(shù)定義為,

    (1)

    式中:xn——恢復力實驗測量值;

    yn——恢復力數(shù)值模擬值;

    N——實驗加載總步數(shù);

    Kxy——相關系數(shù)。

    Kxy的量綱為1,Kxy越接近1就表示兩組曲線的近似程度越高,可以定量評價在不同參數(shù)設定情況下模擬結果的精度。

    a 方案1               b 方案2                c 方案3

    不同纖維個數(shù)劃分方案對應的相關系數(shù)如圖7所示。由圖7得出,相關系數(shù)隨著方案1~方案4減小,方案4~方案5相關系數(shù)增大,在方案4處取得相關系數(shù)最小值0.974 0,在方案1處取得相關系數(shù)最大值0.974 7。相關系數(shù)最小值與最大值的相對誤差為0.071 9%,因此,保護層區(qū)和核心區(qū)混凝土纖維個數(shù)改變對滯回曲線的模擬結果影響不大。隨著纖維數(shù)量的增加,計算負荷也會大大增加。為了提高計算效率,纖維劃分數(shù)量建議優(yōu)選方案1。

    圖7 不同纖維劃分方案與相關系數(shù)

    3.2積分點個數(shù)影響

    采用基于位移法梁柱單元模擬中柱C時,單元個數(shù)為5,纖維截面劃分采用表2中的方案1。同時,需要確定纖維截面的數(shù)量,即積分點的數(shù)量。下面分別采用多種積分點數(shù)量進行數(shù)值模擬。當積分點個數(shù)為1或者個數(shù)超過11,OpenSEES會停止運算,積分點個數(shù)可以取2~10。積分點個數(shù)為3、5、7、9的模擬結果分別對應圖8a~8d所示。圖8中實線為模擬結果,虛線為實驗結果。由圖8看出,模擬得到的滯回曲線與試驗結果在往復加載過程中吻合度較高,模擬結果能反映構件柱的負剛度現(xiàn)象及P-△效應的影響,積分點個數(shù)的改變對滯回曲線形狀的影響并不明顯。

    OpenSEES中基于位移法的梁柱單元采用的是Gauss-Lobatta數(shù)值積分法,通過積分點應用增量法確定各個截面的曲率,最終建立構件在荷載作用下的剛度矩陣。該積分算法計算精度高,應用比較廣泛。對于中柱C,在應用基于位移的梁柱單元進行模擬時,不同積分點個數(shù)對應的相關系數(shù)如圖9所示。由圖9可見,當積分點個數(shù)取3時相關系數(shù)最小,Kxy=0.973 9;當積分點個數(shù)取7時相關系數(shù)最大,Kxy=0.974 5。相關系數(shù)最小值和最大值的相對誤差為0.061 6%,得出相關系數(shù)隨積分點個數(shù)變化不大。通過圖9知,在積分點個數(shù)由3到7,相關系數(shù)隨著增大;積分點個數(shù)由7到9,相關系數(shù)隨著減小。在積分點個數(shù)取5、6、7時,模擬結果較好。

    a 3個積分點

    b 5個積分點

    c 7個積分點

    d 9個積分點

    Fig. 8Comparison of hysteresis loop with different integration point number

    圖9 積分點個數(shù)與相關系數(shù)曲線

    3.3單元個數(shù)影響

    基于位移法的梁柱單元可以準確地模擬構件在非線性地震反應中的受力狀態(tài)。為了使模擬結果更加準確,根據(jù)中柱C的加載形式,采用細分中柱C節(jié)點。不同的數(shù)值模擬需要劃分為不同節(jié)點,根據(jù)細分的節(jié)點定義不同數(shù)量的基于位移法的梁柱單元,通過對不同梁柱單元個數(shù)建模分析,減少模擬結果與實驗結果誤差,為以后進行構件與框架結構建模分析提供參考。采用基于位移法梁柱單元模擬中柱C時,積分點個數(shù)為5,采用不同單元個數(shù)模擬并與實驗結果對比,單元個數(shù)采用1~6、8時,不同單元個數(shù)對應的相關系數(shù)如圖10所示。不同單元的模擬結果如圖11所示,其中虛線為實驗結果,實線為模擬結果。

    圖10 單元數(shù)量與相關系數(shù)曲線

    由圖10看出,對于中柱C應用基于位移法的梁柱單元,當單元個數(shù)由1到5變化時,相關系數(shù)隨之增大,在單元個數(shù)取5時,相關系數(shù)最大值為0.974 5,單元個數(shù)取1時,相關系數(shù)最小值為0.947 7;當單元個數(shù)由5到8變化時,相關系數(shù)隨之減小,可以得出單元個數(shù)為5時,模擬結果最好。相關系數(shù)最小值和相關系數(shù)最大值的相對誤差為3.14%??梢?,與纖維個數(shù)變化和積分點個數(shù)變化相比,單元個數(shù)變化對模擬結果影響比對結果的影響更大,建議建模時選擇5個單元。

    a 1個單元              b 2個單元              c 4個單元

    d 5個單元               e 6個單元              f 8個單元

    4 結 論

    利用中柱C擬靜力試驗數(shù)據(jù),分析了基于位移法纖維梁柱單元建模過程中人為定義參數(shù)對模擬結果的影響程度。采用不同單元個數(shù)、不同積分點個數(shù)、不同混凝土纖維個數(shù)進行了分析,結論如下:

    (1)基于位移的梁柱單元,單元個數(shù)為5時模擬結果較好。積分點個數(shù)不能為1,也不能超過10,積分點個數(shù)在5和7之間時,模擬結果較好。保護層區(qū)和核心區(qū)混凝土纖維個數(shù)改變對滯回曲線的模擬結果影響不大。

    (2)單元個數(shù)變化、纖維截面纖維個數(shù)變化、積分點個數(shù)變化得到的相關系數(shù)最大值與最小值的相對誤差分別為3.140 0%、0.071 9%、0.061 6%,可以得出,單元個數(shù)的改變對模擬結果的影響敏感性較高。

    [1]陳滔, 黃宗明. 鋼筋混凝土框架非彈性地震反應分析模型研究進展[J]. 世界地震工程, 2002, 18(1): 91-97.

    [2]羅文斌, 錢稼茹. RC框架彈塑性位移的解構規(guī)則與構件的目標側移角[J]. 工程力學, 2003, 20(5): 32-36.

    [3]陸新征, 林旭川, 葉列平, 等. 地震下高層建筑連續(xù)倒塌的數(shù)值模型研究[J].工程力學, 2010, 27(11): 64-70.

    [4]解琳琳. 基于OPENSEES的RC構件抗震性能數(shù)值模擬及驗證[D].合肥:合肥工業(yè)大學, 2012: 6-16.

    [5]中國建筑學會建筑結構防倒塌專業(yè)委員會. 鋼筋混凝上框架結構擬靜力倒塌試驗研究及數(shù)值模擬競賽Ⅲ:模擬結果分析[J].建筑結構, 2012, 42(11): 27-30.

    [6]陸新征, 葉列平, 潘鵬, 等. 鋼筋混凝上框架結構擬靜力倒塌試驗研究及數(shù)值模擬競賽—Ⅱ:關鍵構件試驗[J]. 建筑結構, 2012, 42(11): 23-26.

    [7]許建立, 余瓊, 劉金成. 基于OpenSEES的RC梁柱組合件低周反復試驗數(shù)值模擬[J]. 結構工程師, 2013, 29(2): 82-87.

    [8]陸新征, 葉列平, 潘鵬, 等. 鋼筋混凝上框架結構擬靜力倒塌試驗研究及數(shù)值模擬競賽I:框架試驗[J]. 建筑結構, 2012, 42(11): 19-22, 26.

    [9]劉永昌, 關淑賢, 任貴興, 等. 地震模擬振動臺波形再現(xiàn)定量判別法[J]. 地震工程與工程振動, 1998, 18(1): 110-112.

    (編輯荀海鑫)

    OpenSEES-based on numerical analysis of quasi-static test of reinforced concrete column

    WANGTao,SUNYan,MENGLiyan

    (School of Civil Engineering, Heilongjiang University of Science & Technology, Harbin 150022, China)

    This paper presents an experimental numerical analysis motivated by the need for artificial predefinition of model parameters, such as occurs due to the use of displacement-based beam-column element of the OpenSEES software. This research aimed at verifying the effect of parameters on the numerical simulation accuracy builds on model parameters sensitivity analysis using quasi-static test data of RC middle column C. The study demonstrates that numbers of integration points and concrete fibers are not sensitive to simulation effect when the displacement-based beam-column element is used; and numbers of elements provide a greater effect on simulation accuracy than numbers of integration points and concrete fibers. The improved simulation accuracy may result when both the numbers of integration points and the numbers of elements are five and the numbers integration points range between five and seven.

    quasi-static test; OpenSEES; displacement method; beam-column element; hysteretic curves

    2015-09-15

    黑龍江省自然科學基金項目(QC2013C055)

    王濤(1978-),男,黑龍江省集賢人,副教授,博士,研究方向:結構地震反應分析、結構抗震實驗方法與技術,E-mail: hitwangtao@126.com。

    10.3969/j.issn.2095-7262.2016.01.020

    TU398.2

    2095-7262(2016)01-0089-06

    A

    猜你喜歡
    中柱梁柱個數(shù)
    怎樣數(shù)出小正方體的個數(shù)
    等腰三角形個數(shù)探索
    大截面鋼管勁性柱梁柱節(jié)點施工技術的要點分析
    中華建設(2019年8期)2019-09-25 08:26:38
    怎樣數(shù)出小木塊的個數(shù)
    斜腿剛構拱橋貝雷梁柱式支架的安全性驗算
    工程與建設(2019年2期)2019-09-02 01:34:08
    實驗和理論方法比較TWB和CR420/CFRP復合材料中心柱增強體的碰撞試驗結果
    汽車文摘(2018年7期)2018-11-27 02:07:48
    怎樣數(shù)出小正方體的個數(shù)
    鋼框架腹板雙角鋼連接梁柱子結構抗倒塌性能分析
    淺埋單層三跨地下結構在豎向地震作用下的動力響應
    低屈服點鋼板剪力墻梁柱剛接與梁柱鉸接對比分析
    一进一出好大好爽视频| 99热精品在线国产| 国产一区二区三区视频了| 黄色女人牲交| 99国产精品一区二区蜜桃av| 色视频www国产| 亚洲在线自拍视频| 我的老师免费观看完整版| 亚洲美女视频黄频| 性色avwww在线观看| 精品久久久久久成人av| 老司机在亚洲福利影院| 草草在线视频免费看| 国产精品自产拍在线观看55亚洲| 在线观看免费午夜福利视频| 天堂av国产一区二区熟女人妻| 欧美三级亚洲精品| www日本黄色视频网| 精品久久久久久久毛片微露脸| 久9热在线精品视频| 久久草成人影院| www.熟女人妻精品国产| 久久伊人香网站| xxx96com| 黑人操中国人逼视频| 色老头精品视频在线观看| 欧美三级亚洲精品| 精品人妻1区二区| 午夜影院日韩av| 久久亚洲精品不卡| 天堂动漫精品| 国产精品av久久久久免费| 在线观看美女被高潮喷水网站 | 日韩国内少妇激情av| 国产亚洲av高清不卡| 国产三级中文精品| 九九久久精品国产亚洲av麻豆 | 中出人妻视频一区二区| 最近最新中文字幕大全电影3| 丁香六月欧美| 欧美日韩瑟瑟在线播放| 三级国产精品欧美在线观看 | 精品无人区乱码1区二区| 欧美日韩亚洲国产一区二区在线观看| 久久人妻av系列| 国产亚洲精品av在线| 婷婷亚洲欧美| 12—13女人毛片做爰片一| 午夜福利18| 18禁黄网站禁片免费观看直播| 国内揄拍国产精品人妻在线| 国产成人福利小说| 亚洲av熟女| av女优亚洲男人天堂 | 国产精品影院久久| 男女做爰动态图高潮gif福利片| 免费高清视频大片| 免费观看人在逋| 毛片女人毛片| 成人鲁丝片一二三区免费| 真实男女啪啪啪动态图| 亚洲国产日韩欧美精品在线观看 | 成人国产一区最新在线观看| 久久婷婷人人爽人人干人人爱| 天堂√8在线中文| 一个人免费在线观看电影 | 亚洲 欧美 日韩 在线 免费| 午夜福利在线在线| 色综合亚洲欧美另类图片| 久久精品综合一区二区三区| 夜夜爽天天搞| 狂野欧美白嫩少妇大欣赏| svipshipincom国产片| 少妇的逼水好多| 国产精品1区2区在线观看.| 国产三级在线视频| 欧美黄色片欧美黄色片| 搡老妇女老女人老熟妇| 国产淫片久久久久久久久 | 麻豆国产av国片精品| 怎么达到女性高潮| 午夜影院日韩av| 精品免费久久久久久久清纯| 亚洲avbb在线观看| 亚洲中文字幕日韩| 88av欧美| 国产麻豆成人av免费视频| 99国产精品99久久久久| 日韩成人在线观看一区二区三区| 免费搜索国产男女视频| 国产aⅴ精品一区二区三区波| 99视频精品全部免费 在线 | 无限看片的www在线观看| 日韩成人在线观看一区二区三区| 亚洲精华国产精华精| 国产精品 欧美亚洲| 婷婷亚洲欧美| 久久这里只有精品中国| 精品欧美国产一区二区三| 夜夜躁狠狠躁天天躁| 国产极品精品免费视频能看的| 国产精品久久久久久亚洲av鲁大| 偷拍熟女少妇极品色| 十八禁人妻一区二区| 女人被狂操c到高潮| 搞女人的毛片| 一进一出好大好爽视频| 女警被强在线播放| 亚洲专区字幕在线| 欧美在线黄色| 午夜免费激情av| 夜夜躁狠狠躁天天躁| 91字幕亚洲| 听说在线观看完整版免费高清| 欧美黄色淫秽网站| 在线观看66精品国产| 91麻豆精品激情在线观看国产| 男插女下体视频免费在线播放| xxx96com| 18禁观看日本| 国产成人精品无人区| 高清在线国产一区| 两个人的视频大全免费| 亚洲无线在线观看| 亚洲男人的天堂狠狠| 午夜两性在线视频| 90打野战视频偷拍视频| 国产成人啪精品午夜网站| 精品国产超薄肉色丝袜足j| 天堂√8在线中文| 很黄的视频免费| 天天一区二区日本电影三级| 天堂动漫精品| 亚洲人成伊人成综合网2020| 亚洲自偷自拍图片 自拍| 美女午夜性视频免费| 这个男人来自地球电影免费观看| 欧美高清成人免费视频www| 午夜福利在线在线| 成年女人永久免费观看视频| 亚洲国产欧美一区二区综合| 国产成人av教育| 制服人妻中文乱码| 不卡av一区二区三区| 男人舔奶头视频| 在线免费观看的www视频| 嫩草影院入口| 日韩欧美在线二视频| 国产成+人综合+亚洲专区| 国产乱人视频| 亚洲精品色激情综合| 中国美女看黄片| 国产成人啪精品午夜网站| 国产精品国产高清国产av| 亚洲美女黄片视频| 成人特级av手机在线观看| 麻豆成人午夜福利视频| 久久久久久久久免费视频了| 精品国产亚洲在线| 美女黄网站色视频| 可以在线观看的亚洲视频| 国产亚洲精品一区二区www| 亚洲片人在线观看| 麻豆国产av国片精品| 搡老熟女国产l中国老女人| 国产高清激情床上av| 久久九九热精品免费| 小蜜桃在线观看免费完整版高清| 国产高清激情床上av| 99在线视频只有这里精品首页| 亚洲国产日韩欧美精品在线观看 | 日本免费a在线| 亚洲欧美激情综合另类| 亚洲va日本ⅴa欧美va伊人久久| 欧美日本亚洲视频在线播放| 中国美女看黄片| 国产1区2区3区精品| 日韩欧美精品v在线| 嫩草影院入口| 色综合婷婷激情| 国产不卡一卡二| 1024手机看黄色片| h日本视频在线播放| 日韩三级视频一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人澡欧美一区二区| 久久久色成人| 成人国产综合亚洲| 久久久久久久午夜电影| 亚洲欧洲精品一区二区精品久久久| 国产爱豆传媒在线观看| 日日干狠狠操夜夜爽| 日本一本二区三区精品| 一级a爱片免费观看的视频| 国模一区二区三区四区视频 | 国产精品精品国产色婷婷| 亚洲国产精品久久男人天堂| 亚洲乱码一区二区免费版| 嫩草影院入口| 国产免费男女视频| 久久香蕉国产精品| 国产aⅴ精品一区二区三区波| 国产极品精品免费视频能看的| 99国产综合亚洲精品| 国产人伦9x9x在线观看| 神马国产精品三级电影在线观看| 国产精品久久视频播放| 脱女人内裤的视频| av欧美777| 亚洲成av人片免费观看| 脱女人内裤的视频| 欧美在线黄色| 两人在一起打扑克的视频| 女同久久另类99精品国产91| 俄罗斯特黄特色一大片| 国产aⅴ精品一区二区三区波| 免费看a级黄色片| 国产成人影院久久av| 国产熟女xx| 日韩欧美一区二区三区在线观看| 日韩高清综合在线| 91九色精品人成在线观看| 国产精品免费一区二区三区在线| 欧美丝袜亚洲另类 | 最近在线观看免费完整版| 又黄又粗又硬又大视频| 成年女人永久免费观看视频| 亚洲第一电影网av| 欧美在线黄色| 美女高潮的动态| 亚洲成人免费电影在线观看| 成人午夜高清在线视频| 黑人操中国人逼视频| 美女黄网站色视频| svipshipincom国产片| 国产精品久久久久久亚洲av鲁大| а√天堂www在线а√下载| 亚洲在线自拍视频| 大型黄色视频在线免费观看| 青草久久国产| 黑人操中国人逼视频| 欧美黄色片欧美黄色片| 高清在线国产一区| 色噜噜av男人的天堂激情| 亚洲人与动物交配视频| 国产视频一区二区在线看| 国内毛片毛片毛片毛片毛片| 丁香六月欧美| 精华霜和精华液先用哪个| 国产精华一区二区三区| 国产精品99久久99久久久不卡| 九九热线精品视视频播放| 99久久成人亚洲精品观看| 精品人妻1区二区| 两性夫妻黄色片| 日韩精品中文字幕看吧| 免费看美女性在线毛片视频| 天堂网av新在线| 一a级毛片在线观看| 国产成人aa在线观看| av天堂中文字幕网| 亚洲欧美日韩卡通动漫| 好看av亚洲va欧美ⅴa在| 国产精品一区二区三区四区久久| 国产aⅴ精品一区二区三区波| 精品久久久久久,| 亚洲人成电影免费在线| 欧美日韩黄片免| 天天躁日日操中文字幕| 亚洲国产欧美网| 丰满人妻一区二区三区视频av | 天堂网av新在线| 日本 av在线| 国产亚洲欧美98| av在线蜜桃| 亚洲五月天丁香| 国产午夜精品论理片| 长腿黑丝高跟| 国产精品 欧美亚洲| 亚洲自拍偷在线| 日韩有码中文字幕| 身体一侧抽搐| 国产亚洲精品久久久com| 岛国在线观看网站| 黄频高清免费视频| 91在线观看av| 久久久久久久午夜电影| 人人妻人人看人人澡| 在线免费观看的www视频| 精品99又大又爽又粗少妇毛片 | 桃色一区二区三区在线观看| 成人高潮视频无遮挡免费网站| 久久精品国产亚洲av香蕉五月| 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 一级黄色大片毛片| 久久久久久国产a免费观看| 99热6这里只有精品| 黑人操中国人逼视频| avwww免费| 1024手机看黄色片| 中文资源天堂在线| 日韩 欧美 亚洲 中文字幕| 中亚洲国语对白在线视频| 日本一本二区三区精品| 国产乱人视频| 欧美成人性av电影在线观看| 欧美大码av| 大型黄色视频在线免费观看| 两个人视频免费观看高清| 午夜久久久久精精品| 国产爱豆传媒在线观看| 亚洲七黄色美女视频| 久久国产精品人妻蜜桃| 欧美+亚洲+日韩+国产| 亚洲av熟女| xxx96com| 国产高清视频在线播放一区| 亚洲av第一区精品v没综合| 国产综合懂色| 大型黄色视频在线免费观看| 狠狠狠狠99中文字幕| av黄色大香蕉| 一进一出抽搐gif免费好疼| 51午夜福利影视在线观看| 真实男女啪啪啪动态图| 亚洲av成人不卡在线观看播放网| 在线观看美女被高潮喷水网站 | 亚洲av成人不卡在线观看播放网| 69av精品久久久久久| 日日夜夜操网爽| 丁香欧美五月| 男人的好看免费观看在线视频| 亚洲国产欧美网| 国产精品久久久久久精品电影| 精品欧美国产一区二区三| 一本综合久久免费| 色吧在线观看| 热99re8久久精品国产| 日韩精品中文字幕看吧| 亚洲精品色激情综合| 免费在线观看亚洲国产| 一个人看的www免费观看视频| 亚洲熟妇中文字幕五十中出| 97超级碰碰碰精品色视频在线观看| 亚洲专区国产一区二区| 国产精品爽爽va在线观看网站| 国产成人啪精品午夜网站| 国产乱人视频| 色精品久久人妻99蜜桃| 少妇的逼水好多| 巨乳人妻的诱惑在线观看| 禁无遮挡网站| 亚洲av成人一区二区三| 国产伦精品一区二区三区视频9 | 午夜免费激情av| av福利片在线观看| 久久亚洲真实| 亚洲欧美精品综合一区二区三区| 日本a在线网址| 一级毛片女人18水好多| 亚洲国产精品sss在线观看| www.www免费av| 波多野结衣巨乳人妻| 国内揄拍国产精品人妻在线| 国内揄拍国产精品人妻在线| 十八禁人妻一区二区| 91九色精品人成在线观看| 99精品欧美一区二区三区四区| 国产成人影院久久av| 欧美精品啪啪一区二区三区| 波多野结衣高清作品| 亚洲男人的天堂狠狠| 亚洲熟女毛片儿| 成人特级av手机在线观看| 国产精品久久电影中文字幕| 精品午夜福利视频在线观看一区| 中文资源天堂在线| 色噜噜av男人的天堂激情| 成年人黄色毛片网站| 99久久99久久久精品蜜桃| 香蕉av资源在线| 桃色一区二区三区在线观看| 免费av不卡在线播放| 中文字幕熟女人妻在线| 日韩欧美国产在线观看| 男人舔女人的私密视频| 一边摸一边抽搐一进一小说| 欧美日韩国产亚洲二区| 国产高清视频在线观看网站| 国产v大片淫在线免费观看| 久久久久久九九精品二区国产| 久久精品91无色码中文字幕| 好看av亚洲va欧美ⅴa在| 欧美午夜高清在线| 亚洲精品456在线播放app | 亚洲中文字幕日韩| av女优亚洲男人天堂 | 又爽又黄无遮挡网站| 少妇人妻一区二区三区视频| 久久久久国产一级毛片高清牌| 精华霜和精华液先用哪个| 精品99又大又爽又粗少妇毛片 | 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 级片在线观看| 日韩三级视频一区二区三区| 国产亚洲精品av在线| 制服丝袜大香蕉在线| 国内精品一区二区在线观看| 中国美女看黄片| 在线十欧美十亚洲十日本专区| 国产精品爽爽va在线观看网站| 久久久久免费精品人妻一区二区| 久久久久久久精品吃奶| 国产精品99久久99久久久不卡| 一边摸一边抽搐一进一小说| 悠悠久久av| 午夜激情福利司机影院| 村上凉子中文字幕在线| 中出人妻视频一区二区| 国产欧美日韩一区二区三| 亚洲av成人精品一区久久| 国产黄a三级三级三级人| 色综合婷婷激情| 观看美女的网站| 老熟妇乱子伦视频在线观看| 黄片大片在线免费观看| 大型黄色视频在线免费观看| 国产单亲对白刺激| 在线永久观看黄色视频| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| 日韩欧美在线乱码| 国产成人精品无人区| 欧美日韩综合久久久久久 | 网址你懂的国产日韩在线| 中国美女看黄片| 亚洲一区二区三区色噜噜| bbb黄色大片| 亚洲自拍偷在线| 久久天堂一区二区三区四区| 免费在线观看日本一区| 欧美性猛交╳xxx乱大交人| 免费看光身美女| 18美女黄网站色大片免费观看| 久久香蕉精品热| 99精品久久久久人妻精品| av天堂在线播放| 精品午夜福利视频在线观看一区| 宅男免费午夜| 看黄色毛片网站| 欧美性猛交╳xxx乱大交人| 亚洲自偷自拍图片 自拍| 一个人观看的视频www高清免费观看 | 欧美日韩黄片免| 一进一出抽搐动态| 亚洲美女黄片视频| 黄色丝袜av网址大全| 亚洲av成人av| 欧美乱色亚洲激情| 一边摸一边抽搐一进一小说| 高清在线国产一区| 国产精品久久久人人做人人爽| 欧美一级a爱片免费观看看| 欧美性猛交黑人性爽| 两性午夜刺激爽爽歪歪视频在线观看| 少妇裸体淫交视频免费看高清| 亚洲国产精品合色在线| 十八禁人妻一区二区| 久久久久久九九精品二区国产| 欧美丝袜亚洲另类 | 国产欧美日韩精品一区二区| 性色av乱码一区二区三区2| 亚洲色图av天堂| 两性夫妻黄色片| 婷婷精品国产亚洲av在线| 色综合欧美亚洲国产小说| 91麻豆精品激情在线观看国产| 日韩大尺度精品在线看网址| av在线天堂中文字幕| 少妇熟女aⅴ在线视频| 欧美zozozo另类| 在线观看66精品国产| 国产亚洲av嫩草精品影院| 人人妻,人人澡人人爽秒播| 熟女少妇亚洲综合色aaa.| 日本精品一区二区三区蜜桃| 9191精品国产免费久久| 欧美一区二区国产精品久久精品| 成人国产综合亚洲| 午夜精品在线福利| 2021天堂中文幕一二区在线观| 久久伊人香网站| 国产午夜精品久久久久久| 无限看片的www在线观看| 又大又爽又粗| 国产三级中文精品| 美女扒开内裤让男人捅视频| 日韩三级视频一区二区三区| 国内精品一区二区在线观看| 午夜精品在线福利| 亚洲美女黄片视频| 日韩欧美国产在线观看| 久久久国产精品麻豆| 老熟妇乱子伦视频在线观看| 成人性生交大片免费视频hd| 色噜噜av男人的天堂激情| 国内精品久久久久精免费| 日韩欧美免费精品| 亚洲人成网站在线播放欧美日韩| 无人区码免费观看不卡| 国产成人系列免费观看| 制服人妻中文乱码| 最近视频中文字幕2019在线8| 婷婷精品国产亚洲av| 午夜精品久久久久久毛片777| 久久国产乱子伦精品免费另类| 国产激情久久老熟女| 天堂√8在线中文| 欧美日韩国产亚洲二区| 一本久久中文字幕| 精品电影一区二区在线| 人妻久久中文字幕网| 在线免费观看的www视频| 精品久久久久久久人妻蜜臀av| 一夜夜www| 国内精品美女久久久久久| 久久精品国产99精品国产亚洲性色| 美女黄网站色视频| 国产亚洲精品一区二区www| 久久久久九九精品影院| 最好的美女福利视频网| 国产综合懂色| 日韩欧美免费精品| 亚洲电影在线观看av| 性欧美人与动物交配| 欧美黄色片欧美黄色片| 无限看片的www在线观看| 青草久久国产| 国产成人系列免费观看| 人妻丰满熟妇av一区二区三区| 搡老妇女老女人老熟妇| 老汉色∧v一级毛片| 亚洲欧美精品综合久久99| 男女床上黄色一级片免费看| 国产精品乱码一区二三区的特点| 后天国语完整版免费观看| 亚洲欧美日韩无卡精品| 黑人欧美特级aaaaaa片| 黄色视频,在线免费观看| 国产视频一区二区在线看| 亚洲人成伊人成综合网2020| 久久久成人免费电影| 亚洲午夜理论影院| www.www免费av| 亚洲乱码一区二区免费版| 制服人妻中文乱码| 欧美性猛交╳xxx乱大交人| 亚洲人成网站在线播放欧美日韩| 欧美3d第一页| 免费高清视频大片| 婷婷亚洲欧美| 99热6这里只有精品| 国产91精品成人一区二区三区| 国产精品亚洲美女久久久| 毛片女人毛片| 巨乳人妻的诱惑在线观看| 美女被艹到高潮喷水动态| 欧美成人性av电影在线观看| 国产高清视频在线播放一区| 日本在线视频免费播放| av在线天堂中文字幕| 国产一级毛片七仙女欲春2| 特大巨黑吊av在线直播| 久久香蕉国产精品| 欧美日韩精品网址| 国产视频内射| 国产精品久久久人人做人人爽| 中文字幕高清在线视频| 精品久久久久久久末码| 久久人人精品亚洲av| 免费观看精品视频网站| 夜夜爽天天搞| 欧美日本亚洲视频在线播放| 中文字幕久久专区| a级毛片在线看网站| 村上凉子中文字幕在线| 99热这里只有是精品50| 国产高清有码在线观看视频| 欧美+亚洲+日韩+国产| 国产 一区 欧美 日韩| 日日干狠狠操夜夜爽| www.自偷自拍.com| 中文在线观看免费www的网站| 少妇裸体淫交视频免费看高清| 国产高清视频在线播放一区| 精品久久久久久久久久久久久| 最近视频中文字幕2019在线8| 亚洲成人免费电影在线观看| 精品日产1卡2卡| 一个人免费在线观看的高清视频| 亚洲精品在线美女| www.精华液| 叶爱在线成人免费视频播放| 亚洲中文字幕一区二区三区有码在线看 | 欧美一区二区精品小视频在线| 国产精品永久免费网站| 欧美日韩乱码在线| 亚洲 国产 在线| 久久人妻av系列| 国产成人精品无人区| 国产av在哪里看| 99热精品在线国产| 国模一区二区三区四区视频 |