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

    基于響應(yīng)面法的斜拉橋動力模型修正方法

    2022-10-27 00:41:20張立業(yè)王兵見呂竟銘
    公路交通科技 2022年8期
    關(guān)鍵詞:特征參數(shù)斜拉橋修正

    張立業(yè),汪 波,王兵見,陳 可,呂竟銘

    (1.交通運(yùn)輸部公路科學(xué)研究院,北京 100088;2.安徽省公路管理服務(wù)中心,安徽 合肥 230022)

    0 引言

    建立能夠精確反映實(shí)際橋梁結(jié)構(gòu)特性的基準(zhǔn)有限元模型是進(jìn)行結(jié)構(gòu)分析、損傷識別和狀態(tài)評估的關(guān)鍵[1-2]。由于實(shí)際橋梁結(jié)構(gòu)的物理參數(shù)存在隨機(jī)性和離散性,有限元模型采用的物理參數(shù)很難與實(shí)際橋梁結(jié)構(gòu)一致,且建模過程中往往需要對邊界條件和連續(xù)條件進(jìn)行簡化,致使橋梁結(jié)構(gòu)有限元模型計(jì)算的特征參數(shù)與健康監(jiān)測的特征參數(shù)存在一定的差異,進(jìn)而影響橋梁結(jié)構(gòu)分析、損傷識別和狀態(tài)評估的準(zhǔn)確率。獲得能夠精確反映實(shí)際橋梁結(jié)構(gòu)的有限元基準(zhǔn)模型需要借助健康監(jiān)測數(shù)據(jù)和模型修正(Model Updating)技術(shù)。模型修正是通過優(yōu)化結(jié)構(gòu)有限元模型的幾何物理等參數(shù),使有限元計(jì)算的特征參數(shù)接近于結(jié)構(gòu)健康監(jiān)測的特征參數(shù),其核心是最優(yōu)化問題[3-4]。

    常用的橋梁結(jié)構(gòu)模型修正方法包括靜力模型修正法[5-6]和動力模型修正法[7-8]。通過靜力法測得的特征參數(shù)具有準(zhǔn)確率高、抗干擾性強(qiáng)、適應(yīng)性好等優(yōu)點(diǎn),但需要對橋梁結(jié)構(gòu)進(jìn)行荷載試驗(yàn),獲得靜力響應(yīng)數(shù)據(jù),其加載方式單一,實(shí)測數(shù)據(jù)量有限,通常需要中斷交通,因此,采用靜力法進(jìn)行橋梁結(jié)構(gòu)模型修正是比較困難的。動力模型修正法以橋梁結(jié)構(gòu)的模態(tài)參數(shù)作為特征參數(shù),因模態(tài)參數(shù)反映了橋梁結(jié)構(gòu)自身的特性,與加載方式無關(guān),可通過環(huán)境激勵的方法測得,無須中斷交通[9],且大量研究表明:基于模態(tài)參數(shù)的損傷識別和狀態(tài)評估方法對實(shí)際橋梁結(jié)構(gòu)比較有效[10]。因此,動力模型修正方法更適用于實(shí)際橋梁結(jié)構(gòu),在橋梁結(jié)構(gòu)健康監(jiān)測中具有很好的應(yīng)用前景[11-12]。

    有限元模型修正的核心是最優(yōu)化問題,構(gòu)建合理的包含結(jié)構(gòu)主要響應(yīng)信息的目標(biāo)函數(shù)是有限元模型修正的關(guān)鍵問題[13-14],目前已有不少學(xué)者基于實(shí)際橋梁結(jié)構(gòu)或者室內(nèi)試驗(yàn)?zāi)P偷撵o動力測試數(shù)據(jù)[15-16],提出了比較有效的模型修正方法[8,17]。任偉新等[18]提出了一種基于響應(yīng)面的模型修正方法;宗周紅等[19-20]提出了一種基于響應(yīng)面模型修正的橋梁結(jié)構(gòu)識別方法,并應(yīng)用于橋梁健康監(jiān)測的有限元模型確認(rèn)中;馬印平等[21]提出了基于響應(yīng)面法的鋼管混凝土組合桁梁橋多尺度有限元模型修正方法;周林仁等[22]提出了基于徑向基函數(shù)響應(yīng)面方法的斜拉橋有限元模型修正方法;王曉光等[23]將穩(wěn)健估計(jì)法引入響應(yīng)面優(yōu)化求解過程,提高了基于響應(yīng)面模型修正的可靠性;Sanayei等[24]提出了一種基于多響應(yīng)面數(shù)據(jù)的有限元模型修正和損傷識別方法。本研究在斜拉橋健康監(jiān)測和參數(shù)化有限元分析技術(shù)的基礎(chǔ)上,基于響應(yīng)面法建立了斜拉橋動力模型修正方法,解決了斜拉橋有限元動力模型難于修正的問題,研究結(jié)果表明:修正后的有限元模型能夠精確模擬實(shí)際斜拉橋結(jié)構(gòu),對斜拉橋的模型修正具有很好的適應(yīng)性。

    1 模型修正的基本原理

    為了使有限元模型能夠精確地反映實(shí)際橋梁結(jié)構(gòu)的力學(xué)特性,除了建模過程中盡可能準(zhǔn)確地選用符合實(shí)際的結(jié)構(gòu)有限元模型參數(shù),合理地處理邊界條件和連續(xù)條件外,更重要的是要借助試驗(yàn)測試或橋梁監(jiān)測結(jié)果修正有限元模型,這是獲得高精度橋梁結(jié)構(gòu)有限元模型的關(guān)鍵。

    結(jié)構(gòu)健康監(jiān)測數(shù)據(jù)包括結(jié)構(gòu)位移、應(yīng)變、加速度等監(jiān)測指標(biāo),特征參數(shù)指監(jiān)測指標(biāo)的極值、均值等時域特征,以及頻率、振型等頻域特征。由橋梁結(jié)構(gòu)健康監(jiān)測數(shù)據(jù)獲得的1組特征參數(shù)為yi,m(i=1,2,3,…n),由橋梁結(jié)構(gòu)有限元模型計(jì)算得到的1組與yi,m相對應(yīng)的特征參數(shù)為yi,c(i=1,2,3,…n)。由于橋梁結(jié)構(gòu)監(jiān)測指標(biāo)是實(shí)際橋梁結(jié)構(gòu)的真實(shí)響應(yīng),其特征參數(shù)反映了結(jié)構(gòu)的真實(shí)狀態(tài),因此,可以認(rèn)為橋梁結(jié)構(gòu)健康監(jiān)測數(shù)據(jù)獲得的特征參數(shù)yi,m是準(zhǔn)確且可靠的。模型修正的主要任務(wù)是優(yōu)化有限元模型參數(shù)xj(j=1,2,3,…,k),使yi,m和yi,c之間的誤差最小。

    橋梁結(jié)構(gòu)有限元模型計(jì)算得到的yi,c可以看成是有限元模型參數(shù)xj的函數(shù)yi,c(xj),可以將yi,c(xj)和yi,m之間的誤差函數(shù)作為有限元模型修正優(yōu)化問題的目標(biāo)函數(shù)。構(gòu)造目標(biāo)函數(shù)的途徑有最小二乘法、Bayesian概率方法等,其中廣泛應(yīng)用的是最小二乘法。

    (1)

    式中,γi為各特征參數(shù)的權(quán)重,反映了各特征參數(shù)對有限元模型的影響程度,通??梢酝ㄟ^重要度或靈敏度方法來確定。

    2 動力模型修正方法

    振動頻率、振型等模態(tài)參數(shù)是橋梁結(jié)構(gòu)的固有特性,反映了橋梁結(jié)構(gòu)幾何尺寸、質(zhì)量和剛度等主要幾何物理量的分布特征。這些模態(tài)參數(shù)與加載工況無關(guān),避免了監(jiān)測過程中荷載工況差異所帶來的誤差。動力模型修正方法的一般過程包括:(1)建立待修正橋梁結(jié)構(gòu)的初始有限元模型(FEM);(2)通過重要度或靈敏度等分析方法確定各幾何物理參數(shù)權(quán)重并選取待修正參數(shù)xj;(3)根據(jù)待修正參數(shù)的重要度或者靈敏度確定待修正參數(shù)的變化倍率;(4)進(jìn)行FEM計(jì)算,得到與實(shí)測結(jié)果對應(yīng)的特征參數(shù);(5)采用最小二乘法等回歸算法分析FEM的計(jì)算結(jié)果,得出響應(yīng)面方程;(6)建立形如式(1)的目標(biāo)函數(shù),確定各特征參數(shù)的權(quán)重系數(shù);(7)采用最優(yōu)化算法得出目標(biāo)函數(shù)最優(yōu)解,確定待修正參數(shù)的修正值。有限元模型修正的具體流程如圖1所示。

    圖1 有限元模型修正流程Fig.1 Flowchart of finite element modification

    由于橋梁結(jié)構(gòu)動力響應(yīng)與待修正參數(shù)通常存在著難以表達(dá)的隱性函數(shù)關(guān)系,響應(yīng)面法可以通過建立響應(yīng)面方程近似模擬這種隱性函數(shù)關(guān)系,其中,多元二次響應(yīng)面方程是常用的響應(yīng)面方程形式:

    Y=AX,

    (2)

    其中

    Y=[Y1Y2…Yi]T,

    (3)

    X=[1x1x2…xjx1x2x1x3…x1xj…

    (4)

    (5)

    式中,A為響應(yīng)面特征點(diǎn)參數(shù);Y為橋梁結(jié)構(gòu)特征參數(shù)向量;X為待修改正參數(shù)向量;Yi為橋梁結(jié)構(gòu)特征參數(shù);xj為待修正參數(shù);aik為待定系數(shù);i,j,k分別為橋梁結(jié)構(gòu)特征參數(shù)、待修正參數(shù)和多元二次響應(yīng)面方程待定系數(shù)的個數(shù)。

    上述響應(yīng)面方程反映了橋梁結(jié)構(gòu)特征參數(shù)與待修正參數(shù)的函數(shù)關(guān)系,待定系數(shù)需要通過最優(yōu)化方法得到,從而得出響應(yīng)面方程。結(jié)合實(shí)際橋梁結(jié)構(gòu)監(jiān)測的特征參數(shù),便能夠建立起模型修正的目標(biāo)函數(shù)。

    3 基于頻率監(jiān)測的斜拉橋模型修正

    頻率是斜拉橋監(jiān)測的重要指標(biāo),測試成本低且精度高,因此,基于頻率監(jiān)測指標(biāo)來修正斜拉橋有限元模型的可行性較高。這里選用東海大橋主航道斜拉橋作為研究對象,其主跨為420 m,采用雙向六車道加緊急停車帶的高速公路標(biāo)準(zhǔn)。橋?qū)?1.5 m,設(shè)計(jì)時速80 km/h。

    3.1 修正參數(shù)的選取及其變化倍率

    斜拉橋的主梁、橋塔和斜拉索的參數(shù)都會影響到斜拉橋的振動頻率,因此選用主梁鋼材密度、主梁鋼材彈性模量、主梁混凝土密度、主梁混凝土彈性模量、主塔彈性模量和斜拉索彈性模量6個待修正參數(shù)。相對而言,斜拉橋頻率對主梁鋼材的參數(shù)更為敏感,故主梁鋼材密度和彈性模量的變化倍率為1±0.2,其余參數(shù)的變化倍率為1±0.3,具體如表1所示。

    表1 待修正參數(shù)及其變化倍率Tab.1 Parameters to be corrected and their change rates

    3.2 試工況及有限元分析

    采用參數(shù)化有限元分析技術(shù),建立全橋初始有限元模型,選用Beam4單元模擬主梁混凝土,通過模態(tài)分析,得出各階振動頻率計(jì)算結(jié)果,由于低階振動頻率能夠較好地反映橋梁結(jié)構(gòu)狀態(tài),其中豎向振動頻率影響最大,橫向和扭轉(zhuǎn)振動頻率影響次之,且現(xiàn)實(shí)工程中很少發(fā)生。因此,選用前5階振動頻率進(jìn)行修正,包括前3階豎向振動頻率,一階橫向和一階扭轉(zhuǎn)振動頻率。恰當(dāng)?shù)剡x擇頻率修正階數(shù),可使模型修正精度既滿足工程需要,又節(jié)約計(jì)算資源。實(shí)測和修正前頻率如表2所示。

    表2 實(shí)測頻率與修正前頻率Tab.2 Measured frequency and frequency before modification

    有限元分析的振動頻率和振型如圖2所示。

    圖2 模態(tài)分析振型圖Fig.2 Mode shapes for model analysis

    為了建立響應(yīng)面方程,采用中心復(fù)合設(shè)計(jì)法來設(shè)計(jì)響應(yīng)面的試驗(yàn)工況,由于待修正參數(shù)為6個,所以因子數(shù)為6,立方點(diǎn)數(shù)為64,軸向點(diǎn)數(shù)為12,α值為2.828。各代表性試驗(yàn)工況的有限元分析結(jié)果如表3所示,其中工況1~64為立方點(diǎn),工況65~76為軸向點(diǎn),工況77~82為中心點(diǎn)。

    表3 代表性試驗(yàn)工況的有限元分析結(jié)果Tab.3 Finite element analysis result of representative test cases

    由表3可以看出,各階頻率隨著修正參數(shù)的倍率變化而變化,變化倍率越大、變化參數(shù)越多,對各階頻率影響越大,計(jì)算結(jié)果用于構(gòu)建響應(yīng)面方程,包括立方點(diǎn)(工況1~64)、軸向點(diǎn)(工況65~76)和中心點(diǎn)(工況77~82)。

    3.3 響應(yīng)面方程及目標(biāo)函數(shù)

    (6)

    利用響應(yīng)面方程的計(jì)算值與橋梁結(jié)構(gòu)的實(shí)測值,建立的目標(biāo)函數(shù)如下:

    (7)

    式中γi為各特征參數(shù)的權(quán)重值。

    3.4 動力響應(yīng)權(quán)重系數(shù)的確定方法

    由式(7)可以看出,權(quán)重系數(shù)γi是影響目標(biāo)函數(shù)優(yōu)化結(jié)果的重要因素。本研究采用固定變量法分析各權(quán)重系數(shù)對目標(biāo)函數(shù)的影響規(guī)律,用以選擇最佳的動力響應(yīng)權(quán)重系數(shù),進(jìn)而得出模型修正的最優(yōu)化結(jié)果。

    敏感性理論分析可知,式(7)中各動力響應(yīng)權(quán)重系數(shù)對目標(biāo)函數(shù)值的敏感性是不盡相同的,對于不同的工程應(yīng)用情況,由于動力特征不同,其最優(yōu)的權(quán)重系數(shù)也會有所不同,一個比較可行的方法是采用數(shù)值模擬的方法,具體工程具體分析。

    擬定權(quán)重系數(shù)區(qū)間為0~1,采用固定變量法逐個分析γ1,γ2,γ3,γ4,γ5與振動頻率變化率的相關(guān)系數(shù),分析結(jié)果如表4所示。

    表4 權(quán)重系數(shù)與振動頻率變化率的相關(guān)系數(shù)Tab.4 Correlation coefficient between weight coefficient and vibration frequency change rate

    由表4可知,權(quán)重系數(shù)與各階振動頻率的相關(guān)性存在較大差異。對于該橋梁,權(quán)重系數(shù)γ1,γ2,γ3,γ4和γ5分別與豎向二階、橫向一階、橫向一階、豎向二階和豎向三階的相關(guān)性較高,其相關(guān)系數(shù)分別為:-0.756 3,-0.358 3,-0.563 6,-0.236 7和0.190 2。選用相關(guān)系數(shù)的最大值作為權(quán)重系數(shù)可以得到均衡的全局最優(yōu)化結(jié)果,對于有特殊要求的情況,可以選擇關(guān)心階次頻率的相關(guān)系數(shù)最大值,以此保證在關(guān)心階次頻率上得出最優(yōu)化的結(jié)果。

    3.5 模型修正結(jié)果對比

    由前面模態(tài)分析可知,前5階模態(tài)中包含了豎向、橫向和扭轉(zhuǎn)振動頻率。根據(jù)前面分析結(jié)果,這里選用相關(guān)系數(shù)的最大值[-0.756 3 -0.358 3 -0.563 6 -0.236 7 0.190 2]作為權(quán)重系數(shù),得到均衡的全局最優(yōu)化結(jié)果。

    采用最優(yōu)化方法求解目標(biāo)函數(shù),得到的優(yōu)化結(jié)果為x1=0.923,x2=0.922,x3=0.773,x4=0.927,x5=1.384,x6=0.949。修正前與修正后的頻率對比結(jié)果如表5所示??芍?jīng)過模型修正,有限元模型計(jì)算的振動頻率值更接近于實(shí)際橋梁結(jié)構(gòu)監(jiān)測的振動頻率值。

    由圖3可知,通過模型修正,有限元模型計(jì)算的振動頻率值與實(shí)際監(jiān)測的振動頻率值的誤差絕對值由2.03%~7.95%變化到0.05%~1.27%,其中一階振動頻率由-4.59%變化到-0.05%,幾乎與實(shí)際監(jiān)測的振動頻率值一致。

    表5 模型修正前與修正后的頻率對比Tab.5 Comparison of frequencies before and after model modification

    圖3 模型修正結(jié)果Fig.3 Model modification result

    4 結(jié)論

    本研究建立了基于響應(yīng)面法的斜拉橋動力模型修正方法,實(shí)現(xiàn)了斜拉橋有限元模型修正和健康監(jiān)測數(shù)據(jù)快速分析。解決了響應(yīng)面方程建立、目標(biāo)函數(shù)求解和最優(yōu)化算法等斜拉橋動力模型修正的關(guān)鍵問題,得出了斜拉橋動力模型修正的一般流程。實(shí)際橋梁結(jié)構(gòu)算例表明:修正后有限元模型計(jì)算的振動頻率值與實(shí)際橋梁結(jié)構(gòu)監(jiān)測的振動頻率值的誤差由2.03%~7.95%變化到0.05%~1.27%,其中一階振動頻率由-4.59%變化到-0.05%,幾乎與實(shí)際監(jiān)測的振動頻率值一致,修正后的有限元模型能夠精確的模擬實(shí)際斜拉橋結(jié)構(gòu),表明該方法對斜拉橋結(jié)構(gòu)的有限元模型修正具有很好的有效性。

    猜你喜歡
    特征參數(shù)斜拉橋修正
    斜拉橋風(fēng)致振動控制及其未來發(fā)展趨勢
    Some new thoughts of definitions of terms of sedimentary facies: Based on Miall's paper(1985)
    修正這一天
    快樂語文(2021年35期)2022-01-18 06:05:30
    故障診斷中信號特征參數(shù)擇取方法
    基于特征參數(shù)化的木工CAD/CAM系統(tǒng)
    矮塔斜拉橋彈塑性地震響應(yīng)分析
    (112+216+112)m部分斜拉橋設(shè)計(jì)
    合同解釋、合同補(bǔ)充與合同修正
    法律方法(2019年4期)2019-11-16 01:07:28
    基于PSO-VMD的齒輪特征參數(shù)提取方法研究
    軟件修正
    国产精华一区二区三区| 久久国产精品影院| 中文字幕人成人乱码亚洲影| 久久久久久大精品| 精品人妻熟女av久视频| 国产av麻豆久久久久久久| 动漫黄色视频在线观看| 在线观看av片永久免费下载| 18禁黄网站禁片免费观看直播| 色综合亚洲欧美另类图片| 97超视频在线观看视频| 国产精品嫩草影院av在线观看 | 国产亚洲av嫩草精品影院| 老司机午夜福利在线观看视频| 免费一级毛片在线播放高清视频| 永久网站在线| 丰满人妻一区二区三区视频av| 毛片一级片免费看久久久久 | 国产综合懂色| 久久香蕉精品热| 午夜a级毛片| 日本免费a在线| 美女 人体艺术 gogo| 美女被艹到高潮喷水动态| 国产极品精品免费视频能看的| 日韩欧美一区二区三区在线观看| 免费人成视频x8x8入口观看| 高清在线国产一区| 精品国产亚洲在线| 国产国拍精品亚洲av在线观看| 一本久久中文字幕| 又黄又爽又刺激的免费视频.| 在线观看美女被高潮喷水网站 | 亚洲成人免费电影在线观看| 国产精品一及| 国产真实乱freesex| 国产精品一区二区免费欧美| 每晚都被弄得嗷嗷叫到高潮| 久久久久久大精品| 亚洲精品456在线播放app | 日韩成人在线观看一区二区三区| 美女xxoo啪啪120秒动态图 | 亚洲18禁久久av| 欧美黑人巨大hd| 成人av在线播放网站| 国产精品亚洲av一区麻豆| 欧美成人免费av一区二区三区| 国产大屁股一区二区在线视频| 免费在线观看日本一区| 免费看光身美女| 精品国产亚洲在线| 一级a爱片免费观看的视频| 免费高清视频大片| 日韩有码中文字幕| 国产午夜精品论理片| 日本五十路高清| 色综合站精品国产| 欧美丝袜亚洲另类 | 久久久久久九九精品二区国产| 久久国产乱子免费精品| 香蕉av资源在线| 日本a在线网址| 最近在线观看免费完整版| 成年免费大片在线观看| 男女那种视频在线观看| 欧美+亚洲+日韩+国产| 亚洲av日韩精品久久久久久密| 18禁黄网站禁片午夜丰满| 丁香六月欧美| 少妇人妻一区二区三区视频| 亚洲aⅴ乱码一区二区在线播放| 国产精品亚洲av一区麻豆| 韩国av一区二区三区四区| 精品不卡国产一区二区三区| 神马国产精品三级电影在线观看| 中文字幕人成人乱码亚洲影| 美女高潮喷水抽搐中文字幕| 久久人人爽人人爽人人片va | 老司机福利观看| 九色国产91popny在线| 色5月婷婷丁香| 欧美日韩综合久久久久久 | 在线播放无遮挡| av天堂中文字幕网| 欧美日本亚洲视频在线播放| 国产主播在线观看一区二区| 色综合站精品国产| 亚洲欧美精品综合久久99| 中文字幕精品亚洲无线码一区| 精品久久久久久成人av| 最近最新中文字幕大全电影3| av专区在线播放| 99国产精品一区二区三区| 在线天堂最新版资源| 亚洲欧美日韩无卡精品| 免费观看人在逋| 国产淫片久久久久久久久 | 天堂av国产一区二区熟女人妻| 色综合婷婷激情| 欧美精品啪啪一区二区三区| 免费看a级黄色片| 国产一区二区激情短视频| 亚洲av熟女| 久99久视频精品免费| 成年女人毛片免费观看观看9| 中文资源天堂在线| av福利片在线观看| 国产高潮美女av| 中文字幕人成人乱码亚洲影| 亚洲熟妇熟女久久| 性插视频无遮挡在线免费观看| 最近在线观看免费完整版| 99久久精品一区二区三区| 日本免费一区二区三区高清不卡| 欧美午夜高清在线| 国产精品精品国产色婷婷| 亚洲五月婷婷丁香| aaaaa片日本免费| 亚洲国产精品久久男人天堂| 人人妻人人看人人澡| 亚洲 欧美 日韩 在线 免费| www日本黄色视频网| 欧美一区二区国产精品久久精品| 午夜老司机福利剧场| 久久精品国产99精品国产亚洲性色| 中文字幕久久专区| 国产精品1区2区在线观看.| 最好的美女福利视频网| 少妇的逼好多水| www.www免费av| 欧美日韩亚洲国产一区二区在线观看| 人人妻人人看人人澡| 亚洲av二区三区四区| 最好的美女福利视频网| 床上黄色一级片| 久久性视频一级片| 久久久久精品国产欧美久久久| 久久久久国产精品人妻aⅴ院| 麻豆成人av在线观看| 日本与韩国留学比较| 九色成人免费人妻av| 99riav亚洲国产免费| 校园春色视频在线观看| 亚洲专区中文字幕在线| 性插视频无遮挡在线免费观看| 亚洲在线自拍视频| 精品久久久久久久久亚洲 | 国产欧美日韩精品亚洲av| 国内精品久久久久久久电影| 三级男女做爰猛烈吃奶摸视频| 我的女老师完整版在线观看| 97热精品久久久久久| 精品欧美国产一区二区三| 日韩大尺度精品在线看网址| 亚洲精品一区av在线观看| 两个人视频免费观看高清| 一级av片app| 好看av亚洲va欧美ⅴa在| 国产免费一级a男人的天堂| 高清日韩中文字幕在线| 亚洲欧美日韩卡通动漫| 久久精品国产99精品国产亚洲性色| 亚洲av成人精品一区久久| 波野结衣二区三区在线| 亚洲国产欧美人成| 色视频www国产| 赤兔流量卡办理| 一级黄色大片毛片| eeuss影院久久| 亚洲性夜色夜夜综合| 国产欧美日韩一区二区精品| 一级av片app| 中文资源天堂在线| 午夜影院日韩av| 十八禁网站免费在线| 久久午夜福利片| 老司机福利观看| 国产精品亚洲美女久久久| 每晚都被弄得嗷嗷叫到高潮| 亚洲av美国av| 俺也久久电影网| 国产av一区在线观看免费| 成人国产综合亚洲| 亚洲片人在线观看| 一级毛片久久久久久久久女| 青草久久国产| 亚洲男人的天堂狠狠| 特级一级黄色大片| 人人妻,人人澡人人爽秒播| 一二三四社区在线视频社区8| 免费人成视频x8x8入口观看| 欧美极品一区二区三区四区| 日本免费一区二区三区高清不卡| 亚洲精品久久国产高清桃花| 国产成+人综合+亚洲专区| 国产欧美日韩一区二区三| 国产中年淑女户外野战色| 99久久九九国产精品国产免费| 亚洲av美国av| 国产成人啪精品午夜网站| 欧美一区二区精品小视频在线| 搡老妇女老女人老熟妇| 可以在线观看毛片的网站| 久久国产乱子伦精品免费另类| 日本免费a在线| 欧美黄色淫秽网站| 久久精品国产亚洲av涩爱 | 他把我摸到了高潮在线观看| 亚洲一区高清亚洲精品| 91久久精品电影网| 中国美女看黄片| 91在线观看av| 天堂√8在线中文| 亚洲精品456在线播放app | 有码 亚洲区| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品一区av在线观看| 亚洲国产精品成人综合色| 久久伊人香网站| 亚洲经典国产精华液单 | 亚洲欧美日韩高清在线视频| 欧美区成人在线视频| 在线免费观看的www视频| 色精品久久人妻99蜜桃| 国产aⅴ精品一区二区三区波| 亚洲精品在线美女| 国内毛片毛片毛片毛片毛片| 99精品在免费线老司机午夜| 久久久成人免费电影| 很黄的视频免费| 全区人妻精品视频| 久久亚洲真实| 亚洲,欧美,日韩| 亚洲av美国av| 丰满人妻熟妇乱又伦精品不卡| 欧美又色又爽又黄视频| 别揉我奶头~嗯~啊~动态视频| 亚洲色图av天堂| 女生性感内裤真人,穿戴方法视频| 欧美高清性xxxxhd video| 日韩欧美 国产精品| 天美传媒精品一区二区| 9191精品国产免费久久| 天堂√8在线中文| 亚洲精品乱码久久久v下载方式| 午夜福利在线在线| 国产美女午夜福利| 999久久久精品免费观看国产| 国产黄色小视频在线观看| 国产精品电影一区二区三区| 我的女老师完整版在线观看| 日本三级黄在线观看| 日本五十路高清| 一本一本综合久久| 全区人妻精品视频| 男女之事视频高清在线观看| 韩国av一区二区三区四区| 亚洲片人在线观看| 日韩欧美精品免费久久 | 亚洲精品一区av在线观看| 国产黄色小视频在线观看| 99久国产av精品| 欧美zozozo另类| 国产日本99.免费观看| 好看av亚洲va欧美ⅴa在| 18禁黄网站禁片免费观看直播| 嫩草影院入口| 日本在线视频免费播放| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区三区在线臀色熟女| 国产精品永久免费网站| 国产一级毛片七仙女欲春2| 亚洲国产欧美人成| 一进一出抽搐动态| 老司机深夜福利视频在线观看| 长腿黑丝高跟| 欧美日韩亚洲国产一区二区在线观看| 看十八女毛片水多多多| 少妇被粗大猛烈的视频| 午夜福利视频1000在线观看| 国产激情偷乱视频一区二区| 午夜影院日韩av| 天美传媒精品一区二区| 国产精品野战在线观看| 婷婷精品国产亚洲av在线| 久久久久国内视频| 色5月婷婷丁香| 成年女人永久免费观看视频| 脱女人内裤的视频| 久久精品影院6| 国产av一区在线观看免费| 美女免费视频网站| 中文字幕精品亚洲无线码一区| 日本黄大片高清| 中亚洲国语对白在线视频| 久久久久久久久久成人| 黄色视频,在线免费观看| 欧美潮喷喷水| 久久精品国产清高在天天线| 色尼玛亚洲综合影院| 亚洲成人中文字幕在线播放| 久久精品国产99精品国产亚洲性色| 国产精品一区二区性色av| 免费观看人在逋| 成人三级黄色视频| 国产精品免费一区二区三区在线| 99国产综合亚洲精品| 久久中文看片网| 欧美激情国产日韩精品一区| av国产免费在线观看| 97热精品久久久久久| 啦啦啦韩国在线观看视频| 丁香欧美五月| 毛片女人毛片| 亚洲中文字幕一区二区三区有码在线看| 最近最新免费中文字幕在线| 变态另类成人亚洲欧美熟女| 久久天躁狠狠躁夜夜2o2o| 日日摸夜夜添夜夜添av毛片 | 不卡一级毛片| 中文字幕熟女人妻在线| 99热这里只有精品一区| 在线观看av片永久免费下载| 一级作爱视频免费观看| 国产爱豆传媒在线观看| 欧美xxxx性猛交bbbb| 亚洲中文字幕一区二区三区有码在线看| 久久国产精品人妻蜜桃| 国产不卡一卡二| 毛片一级片免费看久久久久 | 免费大片18禁| 日本在线视频免费播放| 老熟妇仑乱视频hdxx| 亚洲七黄色美女视频| 老女人水多毛片| 亚洲成人久久性| 欧美成人一区二区免费高清观看| 免费在线观看日本一区| 精品一区二区免费观看| 国产成人av教育| 在线观看免费视频日本深夜| 精品久久久久久久久av| 午夜福利成人在线免费观看| 中文亚洲av片在线观看爽| 国产精品女同一区二区软件 | 午夜视频国产福利| 3wmmmm亚洲av在线观看| 亚洲人成电影免费在线| 亚洲专区国产一区二区| 少妇的逼好多水| 亚洲不卡免费看| 成人国产一区最新在线观看| 亚洲五月婷婷丁香| 又黄又爽又免费观看的视频| 九色成人免费人妻av| 亚洲avbb在线观看| 国产精品综合久久久久久久免费| 国产精品久久久久久人妻精品电影| 国产在线精品亚洲第一网站| 搡老岳熟女国产| 成人无遮挡网站| 婷婷亚洲欧美| 五月玫瑰六月丁香| 在线播放国产精品三级| 亚洲精品日韩av片在线观看| a在线观看视频网站| 成人永久免费在线观看视频| 99久久精品热视频| 老司机午夜福利在线观看视频| 亚洲第一区二区三区不卡| 听说在线观看完整版免费高清| 老司机福利观看| 99久久久亚洲精品蜜臀av| 天美传媒精品一区二区| 中文字幕久久专区| 免费av毛片视频| av黄色大香蕉| 88av欧美| 日日夜夜操网爽| 欧美成人性av电影在线观看| 国产精品电影一区二区三区| 女生性感内裤真人,穿戴方法视频| 色尼玛亚洲综合影院| 久久久久亚洲av毛片大全| 老鸭窝网址在线观看| 国产精品一区二区免费欧美| 免费在线观看影片大全网站| www.色视频.com| 简卡轻食公司| 婷婷精品国产亚洲av| 久久久久久久午夜电影| 国产探花在线观看一区二区| 国产精品爽爽va在线观看网站| 赤兔流量卡办理| 舔av片在线| 91狼人影院| 国内毛片毛片毛片毛片毛片| 又黄又爽又免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 欧美日韩综合久久久久久 | 一本久久中文字幕| 老熟妇仑乱视频hdxx| 午夜激情福利司机影院| 色在线成人网| 精品日产1卡2卡| 热99在线观看视频| 精品日产1卡2卡| av在线观看视频网站免费| 嫁个100分男人电影在线观看| 在线观看舔阴道视频| 亚洲激情在线av| 国产v大片淫在线免费观看| av天堂在线播放| 欧美色视频一区免费| 三级男女做爰猛烈吃奶摸视频| a级毛片a级免费在线| 亚洲男人的天堂狠狠| 久久欧美精品欧美久久欧美| 国内精品久久久久久久电影| 精品久久久久久久久亚洲 | av天堂中文字幕网| 国产亚洲精品久久久com| 一进一出抽搐gif免费好疼| 美女 人体艺术 gogo| 日日摸夜夜添夜夜添小说| 成年女人永久免费观看视频| 9191精品国产免费久久| 久久6这里有精品| 一级黄色大片毛片| 国产真实乱freesex| 人人妻,人人澡人人爽秒播| 男人狂女人下面高潮的视频| 美女免费视频网站| 欧美3d第一页| 国产老妇女一区| 国产视频一区二区在线看| 久久精品综合一区二区三区| 亚洲人成伊人成综合网2020| 给我免费播放毛片高清在线观看| 天天一区二区日本电影三级| 美女黄网站色视频| 日本三级黄在线观看| 久久久久亚洲av毛片大全| 狂野欧美白嫩少妇大欣赏| 999久久久精品免费观看国产| 亚洲五月婷婷丁香| 麻豆成人av在线观看| 久久久久久久精品吃奶| 成人高潮视频无遮挡免费网站| 九九热线精品视视频播放| 99热精品在线国产| 亚洲,欧美精品.| 黄色丝袜av网址大全| 成人三级黄色视频| 狠狠狠狠99中文字幕| 国产乱人伦免费视频| 禁无遮挡网站| 欧美又色又爽又黄视频| 亚洲男人的天堂狠狠| 美女高潮喷水抽搐中文字幕| 国产亚洲精品综合一区在线观看| 69人妻影院| 九九在线视频观看精品| 可以在线观看的亚洲视频| 久久国产乱子伦精品免费另类| 日本在线视频免费播放| 国产黄a三级三级三级人| 亚洲一区二区三区不卡视频| 中国美女看黄片| 一区二区三区高清视频在线| 成年女人毛片免费观看观看9| 丰满的人妻完整版| 午夜视频国产福利| 亚洲aⅴ乱码一区二区在线播放| 亚洲中文字幕日韩| 免费在线观看亚洲国产| 婷婷六月久久综合丁香| 99久久精品一区二区三区| 性色av乱码一区二区三区2| 欧美日韩综合久久久久久 | 久久天躁狠狠躁夜夜2o2o| 12—13女人毛片做爰片一| 久久久国产成人精品二区| 久久99热这里只有精品18| 可以在线观看毛片的网站| 国产黄色小视频在线观看| 九九在线视频观看精品| 免费在线观看亚洲国产| 国产免费av片在线观看野外av| 欧美潮喷喷水| 亚洲av熟女| 日本黄大片高清| 最新在线观看一区二区三区| 91麻豆精品激情在线观看国产| 国产激情偷乱视频一区二区| 欧美zozozo另类| 身体一侧抽搐| 又爽又黄无遮挡网站| 色吧在线观看| 免费av观看视频| 久久亚洲真实| 中文字幕av成人在线电影| 亚洲av不卡在线观看| 国产精品一区二区三区四区久久| 成人永久免费在线观看视频| 午夜福利视频1000在线观看| 亚洲一区二区三区色噜噜| 在线免费观看的www视频| 老司机午夜十八禁免费视频| 网址你懂的国产日韩在线| 亚洲精华国产精华精| 在线观看66精品国产| av视频在线观看入口| 黄色日韩在线| a在线观看视频网站| 精品一区二区三区人妻视频| 51国产日韩欧美| 大型黄色视频在线免费观看| 国产成人av教育| 欧美又色又爽又黄视频| 久久久久性生活片| 神马国产精品三级电影在线观看| 午夜激情福利司机影院| 久久国产乱子免费精品| 午夜精品一区二区三区免费看| 欧美性感艳星| 久久精品91蜜桃| 久久久国产成人精品二区| 亚州av有码| 好看av亚洲va欧美ⅴa在| 久久久久九九精品影院| 亚洲精品一区av在线观看| 午夜久久久久精精品| 最近视频中文字幕2019在线8| 婷婷精品国产亚洲av| 午夜激情欧美在线| 欧美午夜高清在线| 99久久99久久久精品蜜桃| 成人无遮挡网站| av中文乱码字幕在线| 国产又黄又爽又无遮挡在线| 少妇裸体淫交视频免费看高清| 男女下面进入的视频免费午夜| 亚洲一区高清亚洲精品| 久久精品久久久久久噜噜老黄 | 亚洲欧美激情综合另类| 999久久久精品免费观看国产| 国产成人aa在线观看| 热99在线观看视频| 国产亚洲精品综合一区在线观看| 国产成人影院久久av| 午夜福利18| 夜夜躁狠狠躁天天躁| 亚洲av熟女| 一进一出抽搐gif免费好疼| 中文字幕免费在线视频6| 久久久久久久午夜电影| 日本成人三级电影网站| 色视频www国产| 亚洲国产欧洲综合997久久,| 国产精品一区二区三区四区免费观看 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 麻豆国产av国片精品| 久久精品国产自在天天线| 欧美潮喷喷水| 国产精品久久视频播放| 村上凉子中文字幕在线| 啦啦啦观看免费观看视频高清| 欧美激情在线99| 嫩草影视91久久| 一本精品99久久精品77| 少妇裸体淫交视频免费看高清| 麻豆久久精品国产亚洲av| 老司机福利观看| 丰满乱子伦码专区| 九九在线视频观看精品| 亚洲第一欧美日韩一区二区三区| 麻豆国产97在线/欧美| 三级男女做爰猛烈吃奶摸视频| 夜夜夜夜夜久久久久| 亚洲熟妇中文字幕五十中出| 看十八女毛片水多多多| 97超级碰碰碰精品色视频在线观看| 深夜精品福利| 老熟妇乱子伦视频在线观看| 美女 人体艺术 gogo| 美女黄网站色视频| 欧美激情久久久久久爽电影| 国产中年淑女户外野战色| 久久久久久久久久黄片| 少妇高潮的动态图| 国产日本99.免费观看| 波多野结衣高清作品| 白带黄色成豆腐渣| 欧美色欧美亚洲另类二区| 波多野结衣高清作品| 婷婷精品国产亚洲av| 亚洲一区高清亚洲精品| 色哟哟·www| 日本成人三级电影网站| 中出人妻视频一区二区| 国产精品日韩av在线免费观看| 一级黄色大片毛片| 成年免费大片在线观看| 国产欧美日韩精品亚洲av| av天堂中文字幕网| 女人被狂操c到高潮| 亚洲av免费在线观看| 看十八女毛片水多多多| 国产69精品久久久久777片| 麻豆国产97在线/欧美| 在线天堂最新版资源| 欧美另类亚洲清纯唯美| 俄罗斯特黄特色一大片| 每晚都被弄得嗷嗷叫到高潮|