• <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ù)
    鋼框架腹板雙角鋼連接梁柱子結構抗倒塌性能分析
    淺埋單層三跨地下結構在豎向地震作用下的動力響應
    低屈服點鋼板剪力墻梁柱剛接與梁柱鉸接對比分析
    啦啦啦啦在线视频资源| 国产深夜福利视频在线观看| 中文字幕最新亚洲高清| 亚洲精品日韩在线中文字幕| 国产精品国产av在线观看| 亚洲人成网站在线播| 成人无遮挡网站| 一边亲一边摸免费视频| 久久99精品国语久久久| 性高湖久久久久久久久免费观看| 丁香六月天网| 精品人妻熟女av久视频| 日韩av免费高清视频| 日本av免费视频播放| 国产一区二区三区综合在线观看 | 免费观看无遮挡的男女| 最新中文字幕久久久久| 久久久久久久久久成人| 国产精品一二三区在线看| 丰满乱子伦码专区| av播播在线观看一区| 亚洲在久久综合| 久久久久久人妻| 一区二区三区乱码不卡18| 大片电影免费在线观看免费| 国产女主播在线喷水免费视频网站| 国产不卡av网站在线观看| 婷婷成人精品国产| 一区二区三区免费毛片| 一区二区av电影网| 国产精品.久久久| 99久久人妻综合| 狂野欧美白嫩少妇大欣赏| 国产成人av激情在线播放 | 国产av国产精品国产| 少妇的逼好多水| 哪个播放器可以免费观看大片| 九九在线视频观看精品| 久久久久久久久久成人| 精品人妻一区二区三区麻豆| av不卡在线播放| 综合色丁香网| 中国三级夫妇交换| 亚洲精品视频女| 久久久国产精品麻豆| 蜜桃在线观看..| 日日啪夜夜爽| 国产乱来视频区| 精品久久久精品久久久| 午夜激情久久久久久久| 国产 一区精品| 高清午夜精品一区二区三区| 午夜免费观看性视频| 如何舔出高潮| 热99国产精品久久久久久7| 国产日韩一区二区三区精品不卡 | 一级二级三级毛片免费看| 亚洲伊人久久精品综合| 最新的欧美精品一区二区| 久久精品人人爽人人爽视色| 久热这里只有精品99| 丁香六月天网| 精品人妻在线不人妻| 妹子高潮喷水视频| 人妻人人澡人人爽人人| 国产高清国产精品国产三级| 亚洲美女视频黄频| 日韩av在线免费看完整版不卡| 一个人看视频在线观看www免费| 狠狠精品人妻久久久久久综合| 日本欧美国产在线视频| 国产精品久久久久久av不卡| 亚洲欧洲日产国产| 精品久久久久久电影网| 2022亚洲国产成人精品| 国产精品久久久久久久电影| 精品人妻熟女毛片av久久网站| 日韩视频在线欧美| 99久久人妻综合| 涩涩av久久男人的天堂| 999精品在线视频| 精品人妻熟女av久视频| 天堂8中文在线网| 国产男女内射视频| 亚洲国产av新网站| 高清欧美精品videossex| 黄片播放在线免费| 2018国产大陆天天弄谢| 欧美bdsm另类| 交换朋友夫妻互换小说| 色网站视频免费| www.av在线官网国产| 人人妻人人爽人人添夜夜欢视频| 又黄又爽又刺激的免费视频.| 韩国高清视频一区二区三区| 女性被躁到高潮视频| 欧美日韩国产mv在线观看视频| 一级毛片我不卡| 18禁在线播放成人免费| 日本猛色少妇xxxxx猛交久久| 亚洲国产av影院在线观看| 国产精品国产三级专区第一集| 日韩av在线免费看完整版不卡| 肉色欧美久久久久久久蜜桃| 国产一区二区三区综合在线观看 | 精品久久久噜噜| 亚洲欧美色中文字幕在线| 中文字幕制服av| 国产视频内射| 久久精品夜色国产| 久久久久网色| 能在线免费看毛片的网站| 亚洲婷婷狠狠爱综合网| 国产精品人妻久久久影院| 菩萨蛮人人尽说江南好唐韦庄| 制服诱惑二区| 欧美精品一区二区大全| 国产男女超爽视频在线观看| 久久久国产精品麻豆| 亚洲久久久国产精品| 日韩一区二区视频免费看| 日本欧美视频一区| 欧美精品高潮呻吟av久久| 国产精品久久久久久久久免| 久久99热6这里只有精品| 少妇被粗大猛烈的视频| 性高湖久久久久久久久免费观看| 亚洲欧美成人精品一区二区| 日韩电影二区| 免费少妇av软件| 精品亚洲成国产av| 建设人人有责人人尽责人人享有的| 欧美日韩国产mv在线观看视频| 人人妻人人添人人爽欧美一区卜| 日韩强制内射视频| 五月玫瑰六月丁香| 七月丁香在线播放| 99热网站在线观看| 日韩制服骚丝袜av| 七月丁香在线播放| 久久久久视频综合| 18禁裸乳无遮挡动漫免费视频| 亚洲怡红院男人天堂| 亚洲国产av新网站| 女人精品久久久久毛片| 我的女老师完整版在线观看| 亚洲精品色激情综合| 亚洲欧洲日产国产| 日本猛色少妇xxxxx猛交久久| 一二三四中文在线观看免费高清| 日日摸夜夜添夜夜爱| 能在线免费看毛片的网站| 久久韩国三级中文字幕| 久久午夜福利片| 亚洲婷婷狠狠爱综合网| 亚洲婷婷狠狠爱综合网| 能在线免费看毛片的网站| 欧美xxxx性猛交bbbb| 欧美日韩av久久| 亚洲精品成人av观看孕妇| 蜜桃国产av成人99| 丝瓜视频免费看黄片| 久久久久久久大尺度免费视频| 久久久久久久久久久丰满| 日韩视频在线欧美| 大又大粗又爽又黄少妇毛片口| 国产高清三级在线| 日韩一区二区视频免费看| 欧美人与善性xxx| 女性生殖器流出的白浆| 国产成人精品福利久久| 免费av不卡在线播放| 免费观看在线日韩| 日韩av不卡免费在线播放| 欧美精品人与动牲交sv欧美| 母亲3免费完整高清在线观看 | 亚洲精品国产av成人精品| 国产精品熟女久久久久浪| 亚洲欧美日韩另类电影网站| 韩国av在线不卡| 天堂俺去俺来也www色官网| 91精品国产九色| 婷婷成人精品国产| 自线自在国产av| 伊人久久国产一区二区| 日本欧美视频一区| 天天影视国产精品| 成人午夜精彩视频在线观看| 麻豆成人av视频| 人人妻人人澡人人爽人人夜夜| 免费高清在线观看日韩| 亚洲国产日韩一区二区| 国产免费又黄又爽又色| 91成人精品电影| 丰满饥渴人妻一区二区三| 国产日韩欧美视频二区| 熟妇人妻不卡中文字幕| 青春草视频在线免费观看| 日韩在线高清观看一区二区三区| 国产精品一区二区三区四区免费观看| 亚洲av二区三区四区| 亚洲国产精品一区二区三区在线| 成人毛片a级毛片在线播放| 亚洲av中文av极速乱| 亚洲av成人精品一二三区| 国产伦精品一区二区三区视频9| 色婷婷av一区二区三区视频| 2018国产大陆天天弄谢| 亚洲av电影在线观看一区二区三区| 亚洲精品久久午夜乱码| av电影中文网址| 亚洲av成人精品一区久久| 国产高清不卡午夜福利| 成人手机av| 97在线视频观看| 国产av国产精品国产| 狂野欧美白嫩少妇大欣赏| 一个人免费看片子| videosex国产| 国产精品偷伦视频观看了| 各种免费的搞黄视频| 欧美性感艳星| 曰老女人黄片| 亚洲图色成人| 午夜福利,免费看| 人妻夜夜爽99麻豆av| 一区二区av电影网| 日韩成人av中文字幕在线观看| 国产亚洲精品久久久com| 亚洲人成网站在线播| 国产精品人妻久久久久久| 中文天堂在线官网| 国产精品秋霞免费鲁丝片| 精品久久蜜臀av无| 国产免费视频播放在线视频| 久久久久网色| 极品少妇高潮喷水抽搐| 亚洲精品一二三| 国产成人精品无人区| 久久久久久伊人网av| 啦啦啦中文免费视频观看日本| 少妇 在线观看| 中文字幕久久专区| 国产有黄有色有爽视频| 久久99蜜桃精品久久| 日韩av免费高清视频| 国语对白做爰xxxⅹ性视频网站| 国产无遮挡羞羞视频在线观看| 免费黄频网站在线观看国产| 女人久久www免费人成看片| 街头女战士在线观看网站| 黑人高潮一二区| 成人18禁高潮啪啪吃奶动态图 | 午夜激情福利司机影院| 亚洲人成网站在线播| 两个人免费观看高清视频| 美女脱内裤让男人舔精品视频| 亚洲av.av天堂| 国产精品熟女久久久久浪| 亚洲少妇的诱惑av| 伊人久久精品亚洲午夜| 91午夜精品亚洲一区二区三区| 草草在线视频免费看| 色视频在线一区二区三区| 久久精品夜色国产| 精品亚洲成国产av| 成人午夜精彩视频在线观看| 亚洲经典国产精华液单| 国产精品国产三级专区第一集| 成人国产av品久久久| 久久精品国产鲁丝片午夜精品| 日本91视频免费播放| 日韩熟女老妇一区二区性免费视频| 91aial.com中文字幕在线观看| 午夜久久久在线观看| 亚洲成人手机| 欧美成人精品欧美一级黄| 国产精品国产av在线观看| 欧美精品人与动牲交sv欧美| 色吧在线观看| 午夜福利影视在线免费观看| 99热这里只有精品一区| 大香蕉久久网| 日本午夜av视频| 国产免费视频播放在线视频| 亚洲不卡免费看| 午夜福利影视在线免费观看| 欧美成人精品欧美一级黄| 在线精品无人区一区二区三| 午夜av观看不卡| 一级二级三级毛片免费看| 国产无遮挡羞羞视频在线观看| 亚洲欧美成人精品一区二区| 人成视频在线观看免费观看| 国产男人的电影天堂91| 国产精品久久久久成人av| 乱码一卡2卡4卡精品| 久久久久久久久久久免费av| 观看av在线不卡| 91精品伊人久久大香线蕉| av黄色大香蕉| 99国产综合亚洲精品| 日韩不卡一区二区三区视频在线| 一区二区三区免费毛片| 国产一区有黄有色的免费视频| 成年美女黄网站色视频大全免费 | 国产精品国产三级国产专区5o| 一本大道久久a久久精品| 51国产日韩欧美| 插阴视频在线观看视频| 美女视频免费永久观看网站| 国产精品久久久久久久久免| av视频免费观看在线观看| 99精国产麻豆久久婷婷| 汤姆久久久久久久影院中文字幕| 亚洲综合色网址| 日韩 亚洲 欧美在线| 精品一品国产午夜福利视频| 天堂中文最新版在线下载| 国产精品无大码| 精品午夜福利在线看| 韩国高清视频一区二区三区| 亚洲精品美女久久av网站| 亚洲av男天堂| 久久人人爽人人片av| 91精品三级在线观看| 极品人妻少妇av视频| 91精品伊人久久大香线蕉| 亚洲三级黄色毛片| 免费播放大片免费观看视频在线观看| 亚洲国产精品成人久久小说| 亚洲国产日韩一区二区| 有码 亚洲区| 水蜜桃什么品种好| 91精品三级在线观看| 中国三级夫妇交换| 日韩亚洲欧美综合| 日韩在线高清观看一区二区三区| 如日韩欧美国产精品一区二区三区 | 永久网站在线| 欧美性感艳星| tube8黄色片| 汤姆久久久久久久影院中文字幕| 日本爱情动作片www.在线观看| av黄色大香蕉| 制服丝袜香蕉在线| 日本wwww免费看| 久久午夜综合久久蜜桃| 水蜜桃什么品种好| 国产深夜福利视频在线观看| 日韩在线高清观看一区二区三区| 亚洲av日韩在线播放| 成人国产麻豆网| 久久久久久久精品精品| 狂野欧美白嫩少妇大欣赏| 成人综合一区亚洲| 久久鲁丝午夜福利片| 久久久亚洲精品成人影院| 精品午夜福利在线看| 人妻人人澡人人爽人人| 精品视频人人做人人爽| 亚洲精品日韩av片在线观看| 国产精品熟女久久久久浪| 99re6热这里在线精品视频| 精品人妻偷拍中文字幕| 久久97久久精品| 九九爱精品视频在线观看| 99国产综合亚洲精品| 午夜激情久久久久久久| 国产精品久久久久久精品古装| 99国产综合亚洲精品| 成人国语在线视频| 狂野欧美白嫩少妇大欣赏| 在线精品无人区一区二区三| 日本vs欧美在线观看视频| 国产精品99久久99久久久不卡 | 肉色欧美久久久久久久蜜桃| 天美传媒精品一区二区| 青春草视频在线免费观看| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| 十八禁网站网址无遮挡| 在线观看三级黄色| 亚洲成人一二三区av| 国产成人精品久久久久久| 成人漫画全彩无遮挡| 内地一区二区视频在线| 亚洲精品自拍成人| 午夜91福利影院| 80岁老熟妇乱子伦牲交| 成人手机av| 美女内射精品一级片tv| 精品一区二区三卡| 制服诱惑二区| 久久午夜福利片| 大又大粗又爽又黄少妇毛片口| 国产精品国产三级专区第一集| 久久午夜福利片| 午夜免费男女啪啪视频观看| 亚洲国产精品专区欧美| 狂野欧美激情性bbbbbb| 日韩一区二区视频免费看| 亚洲精品,欧美精品| 国模一区二区三区四区视频| 又黄又爽又刺激的免费视频.| 日韩一本色道免费dvd| 国产高清有码在线观看视频| 久久久国产欧美日韩av| 亚洲精品视频女| 王馨瑶露胸无遮挡在线观看| 天天影视国产精品| 91国产中文字幕| 亚洲欧美清纯卡通| 少妇被粗大的猛进出69影院 | 国产一区有黄有色的免费视频| 夜夜看夜夜爽夜夜摸| 一本色道久久久久久精品综合| 久久午夜福利片| 精品一区二区三卡| 在线天堂最新版资源| 色视频在线一区二区三区| 国产视频内射| 亚洲性久久影院| 少妇被粗大猛烈的视频| 国精品久久久久久国模美| 精品亚洲乱码少妇综合久久| 午夜日本视频在线| 插阴视频在线观看视频| 国产不卡av网站在线观看| 日韩大片免费观看网站| 亚洲精品中文字幕在线视频| 成人毛片60女人毛片免费| 日韩成人av中文字幕在线观看| 在线观看人妻少妇| 91国产中文字幕| 最后的刺客免费高清国语| 国产成人精品一,二区| 91aial.com中文字幕在线观看| 国内精品宾馆在线| 国产毛片在线视频| 亚洲精品乱久久久久久| 欧美xxxx性猛交bbbb| 久久女婷五月综合色啪小说| 一级a做视频免费观看| 欧美一级a爱片免费观看看| 草草在线视频免费看| 亚洲av不卡在线观看| 国产精品偷伦视频观看了| 插阴视频在线观看视频| 人人妻人人爽人人添夜夜欢视频| 亚洲精品乱久久久久久| 欧美日韩亚洲高清精品| 日韩人妻高清精品专区| 男男h啪啪无遮挡| 亚洲欧洲精品一区二区精品久久久 | 欧美三级亚洲精品| 精品久久蜜臀av无| 亚洲精品成人av观看孕妇| 日日摸夜夜添夜夜添av毛片| 黄片无遮挡物在线观看| 91精品国产国语对白视频| 老女人水多毛片| 91午夜精品亚洲一区二区三区| 欧美国产精品一级二级三级| 久久久精品区二区三区| 看十八女毛片水多多多| 久久久久久伊人网av| 男人爽女人下面视频在线观看| 99re6热这里在线精品视频| 欧美人与性动交α欧美精品济南到 | 国产探花极品一区二区| 高清av免费在线| 欧美日韩亚洲高清精品| 熟女电影av网| 狂野欧美激情性xxxx在线观看| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 日日摸夜夜添夜夜爱| 免费人成在线观看视频色| 国产精品无大码| 插阴视频在线观看视频| 国产又色又爽无遮挡免| 久久久久久久久大av| 韩国高清视频一区二区三区| 国产白丝娇喘喷水9色精品| 久久人人爽人人爽人人片va| 欧美亚洲日本最大视频资源| 夜夜看夜夜爽夜夜摸| 国产综合精华液| 夜夜爽夜夜爽视频| 欧美日韩成人在线一区二区| 久久午夜福利片| 一级毛片aaaaaa免费看小| 欧美精品一区二区免费开放| 91精品伊人久久大香线蕉| 中文字幕久久专区| 欧美日韩视频精品一区| 人妻少妇偷人精品九色| 国产精品.久久久| 91久久精品国产一区二区三区| 丁香六月天网| 只有这里有精品99| 一本久久精品| 免费看光身美女| 久久久久久久久大av| 九色亚洲精品在线播放| 大香蕉久久网| 激情五月婷婷亚洲| 精品亚洲乱码少妇综合久久| 18在线观看网站| 嫩草影院入口| 午夜视频国产福利| 久久久久久久久大av| 免费看av在线观看网站| 国产免费一级a男人的天堂| 一区二区三区四区激情视频| 日韩欧美一区视频在线观看| 亚洲av成人精品一二三区| 久久精品夜色国产| 国产成人精品久久久久久| 国产成人av激情在线播放 | 中文乱码字字幕精品一区二区三区| 99久久精品一区二区三区| 亚洲欧洲国产日韩| 久久久久国产精品人妻一区二区| 国产成人精品在线电影| 久久影院123| 亚洲国产毛片av蜜桃av| 亚洲av成人精品一二三区| 蜜桃久久精品国产亚洲av| 免费高清在线观看视频在线观看| av女优亚洲男人天堂| 永久网站在线| 亚洲高清免费不卡视频| 黄色一级大片看看| 青春草亚洲视频在线观看| 亚洲丝袜综合中文字幕| 国产男人的电影天堂91| 丝袜喷水一区| 亚洲av.av天堂| 中文欧美无线码| 一级黄片播放器| 22中文网久久字幕| 老司机影院成人| 欧美日韩一区二区视频在线观看视频在线| 精品卡一卡二卡四卡免费| 久久鲁丝午夜福利片| 久久午夜综合久久蜜桃| 一级爰片在线观看| 人人妻人人澡人人看| 亚洲少妇的诱惑av| 一本色道久久久久久精品综合| 国产一区有黄有色的免费视频| www.色视频.com| 性色av一级| 亚洲一区二区三区欧美精品| 亚洲四区av| 夜夜看夜夜爽夜夜摸| 精品少妇内射三级| 免费大片黄手机在线观看| 久久综合国产亚洲精品| 欧美激情 高清一区二区三区| 精品一区在线观看国产| 两个人免费观看高清视频| 国产精品免费大片| 成年女人在线观看亚洲视频| 国产亚洲欧美精品永久| 天天影视国产精品| 大码成人一级视频| 黄色欧美视频在线观看| 在线观看www视频免费| av.在线天堂| 国产亚洲精品久久久com| 国产女主播在线喷水免费视频网站| 久久久久网色| 久久热精品热| 国产片特级美女逼逼视频| 免费观看av网站的网址| 春色校园在线视频观看| 国产精品久久久久久久电影| 精品久久久精品久久久| 在线天堂最新版资源| 一区二区三区免费毛片| 久久女婷五月综合色啪小说| 91在线精品国自产拍蜜月| 国产探花极品一区二区| 人人妻人人添人人爽欧美一区卜| 久久av网站| 成人漫画全彩无遮挡| 色5月婷婷丁香| 国产成人av激情在线播放 | 精品国产一区二区三区久久久樱花| 黄色毛片三级朝国网站| av有码第一页| 精品少妇内射三级| 夫妻性生交免费视频一级片| 美女内射精品一级片tv| 国产男女超爽视频在线观看| 亚洲婷婷狠狠爱综合网| 老司机影院毛片| 性色avwww在线观看| 久久久久久久久久人人人人人人| 亚洲高清免费不卡视频| 国产高清不卡午夜福利| 性色avwww在线观看| 亚洲第一区二区三区不卡| 国内精品宾馆在线| 成人毛片60女人毛片免费| 曰老女人黄片| 亚洲精品乱码久久久久久按摩| 在线免费观看不下载黄p国产| 欧美日韩成人在线一区二区| 亚洲av综合色区一区| 制服诱惑二区| 熟妇人妻不卡中文字幕|