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

    基于ANSYS平臺的橋梁時域顫振分析

    2011-01-24 01:37:20卿前志
    土木工程與管理學(xué)報 2011年4期
    關(guān)鍵詞:脈沖響應(yīng)表達(dá)式時域

    卿前志

    (上海市政工程設(shè)計研究總院(集團)有限公司, 上海 200092)

    隨著橋梁結(jié)構(gòu)跨徑的日益增加,結(jié)構(gòu)的顫振穩(wěn)定性能成為設(shè)計過程中重要的控制因素之一。對結(jié)構(gòu)顫振穩(wěn)定性能研究的主要目的是確定結(jié)構(gòu)的顫振臨界狀態(tài),揭示顫振發(fā)生的機理,從而保證結(jié)構(gòu)的抗風(fēng)穩(wěn)定性。目前評價橋梁的顫振性能主要有三種方法,即經(jīng)典理論法、風(fēng)洞試驗測試法和試驗加理論法[1, 2]。

    至今為止,結(jié)構(gòu)顫振分析求解方法包含頻域和時域兩大類方法。頻域范圍內(nèi)對顫振問題進行求解的基本思路為:首先將氣動力用Scanlan自激力表達(dá)式[3]進行表達(dá),然后通過求解結(jié)構(gòu)運動特征方程得到結(jié)構(gòu)顫振臨界風(fēng)速和顫振頻率。頻域分析屬于線性分析方法,與之相比,時域分析能方便地考慮各類非線性因素的影響,并反映結(jié)構(gòu)在顫振后的振幅演變規(guī)律,有利于基于過程性能的橋梁氣動性能研究,因而近年來發(fā)展迅速。在時域范圍對橋梁顫振穩(wěn)定性進行求解首先需要得到時域化的氣動自激力表達(dá)式,再在動力有限元中進行求解[2]。對于明顯存在流動分離的鈍體橋梁斷面而言,以顫振導(dǎo)數(shù)表示的自激力是一種時頻混合表達(dá)式,通常不能直接用于時域顫振分析,需要轉(zhuǎn)換成等效的時域表達(dá)。目前,顫振時域分析氣動力的表達(dá)方法有兩種。其一是通過脈沖響應(yīng)函數(shù)結(jié)合有理函數(shù)進行表達(dá)[4, 5];其二是采用階躍函數(shù)結(jié)合橋梁斷面的氣動導(dǎo)數(shù)進行自激力模擬[6]。本文采用第一種形式,首先得到橋梁斷面的脈沖響應(yīng)函數(shù)所表達(dá)的氣動自激力,再對其進行動力有限元求解。

    由于ANSYS軟件具有良好的二次開發(fā)功能,可以應(yīng)用于橋梁頻域顫振分析中[7],但相關(guān)文獻(xiàn)對ANSYS應(yīng)用于橋梁時域顫振分析的報道較少。文獻(xiàn)[8]中雖然提出了一種在ANSYS中實現(xiàn)顫振時程分析的方法,但是,其必須同時對頻率和風(fēng)速進行搜索,因此實質(zhì)上仍然是一種頻域方法。本文提供了在ANSYS中實現(xiàn)顫振時程分析的另一方法,在瞬態(tài)動力學(xué)分析過程中嵌入氣動自激力的遞推算法來考慮其運動狀態(tài)歷史記憶特性,得到結(jié)構(gòu)的顫振臨界風(fēng)速。本文的方法也適合其他風(fēng)振問題,如顫抖振時域分析中的氣動自激力模擬。

    1 氣動自激力脈沖響應(yīng)函數(shù)表達(dá)

    結(jié)構(gòu)體系的運動方程為:

    (1)

    式中[M]、[C]、[K]分別為結(jié)構(gòu)質(zhì)量矩陣、阻尼矩陣和剛度矩陣,{Fa}為結(jié)構(gòu)上受到的荷載列向量,當(dāng)用于顫振分析時,其可以用氣動自激力的形式表達(dá)。忽略斷面?zhèn)认蛘駝拥挠绊?,斷面單位長度氣動自激力可采用如下形式:

    Lse=

    (2)

    Mse=

    (3)

    當(dāng)采用脈沖響應(yīng)函數(shù)對斷面自激力進行表達(dá)時,同樣,忽略斷面?zhèn)认蛘駝拥挠绊?,將豎向運動和扭轉(zhuǎn)運動對斷面自激力的貢獻(xiàn)分開,則斷面單位長度氣動自激力可表示為:

    Lae=Lsea(t)+Lseh(t)

    (4)

    Mae=Msea(t)+Mseh(t)

    (5)

    式中Ifx(f=M,L;x=α,h)為相應(yīng)的脈沖響應(yīng)函數(shù),其直觀意義為單位脈沖位移引起的氣動自激力。

    根據(jù)脈沖響應(yīng)函數(shù)表達(dá)的氣動自激力與Scanlan氣動自激力表達(dá)式兩者頻譜特性一致,可得脈沖響應(yīng)函數(shù)與顫振導(dǎo)數(shù)之間的關(guān)系式:

    (6)

    (7)

    (8)

    (9)

    (10)

    式中,ALh1、ALh2、ALh3、ALhi、dLhi(dLhi≥0;i=4,…,m)都是待擬合的參數(shù)。第一項和第二項分別表示由位移項引起的氣動力和速度項引起的氣動力,第三項表示由加速度項引起的氣動力,該項通常較小而被忽略,第四項用于描述滯后于速度項的氣動力非定常部分,m的大小決定了這種近似的精度和附加方程的數(shù)量,所有這些參數(shù)可以通過試驗提取的顫振導(dǎo)數(shù)進行非線性擬合而確定[9]。

    當(dāng)?shù)玫接欣砗瘮?shù)的相關(guān)系數(shù)后,對式(10)做傅里葉變換,可得到脈沖響應(yīng)函數(shù)的表達(dá)式:

    (11)

    代入式(4)中,有:

    Lseh(t)=

    (12)

    對上式中的積分項進行一次分步積分,可以得到Lseh(t)的最終表達(dá)式如下:

    (13)

    有了式(13)的表達(dá)式后,在進行計算時,必須將卷積積分表達(dá)式轉(zhuǎn)化為遞推關(guān)系表達(dá)式。從式(13)易知,只需要解決Lseh2(t)的計算問題,該式可通過以下方法完成遞推:

    (14)

    Zi(t+Δt)

    (15)

    同理可以得出Lseα(s)、Mseα(s)、Mseh(s)相應(yīng)的表達(dá)式。

    得到斷面的氣動自激力表達(dá)式后,即可方便地在ANSYS中實現(xiàn)顫振時程分析。在ANSYS中節(jié)點速度和加速度的獲得,參照ANSYS幫助文件,對于采用Newmark法的時程分析,節(jié)點加速度和速度可表示為如下差分格式:

    (16)

    (17)

    式中α=1/4(1+γ)2,δ=1/2+γ,γ為振幅衰減系數(shù),即數(shù)值阻尼因子,ANSYS中默認(rèn)為0.005,{un}及{un+1}為結(jié)構(gòu)前一時刻運動狀態(tài)量和結(jié)構(gòu)當(dāng)前時刻運動狀態(tài)量,Δt為計算時間步長,根據(jù)式(16)、(17)即可通過前一時刻的速度、加速度及當(dāng)前時刻的位移求得當(dāng)前時刻的速度和加速度大小,逐步提高風(fēng)速,便可得到結(jié)構(gòu)的顫振臨界風(fēng)速。圖1給出了在ANSYS有限元分析軟件中實現(xiàn)顫振時程分析的計算流程。

    圖1 ANSYS顫振時程分析流程

    2 數(shù)值算例

    為驗證本文計算模型和求解方法的正確性,采用文獻(xiàn)[7]中具有理想平板斷面的簡支板梁進行ANSYS顫振時程分析驗證。模型參數(shù):長度L=300 m,寬40 m,約束兩端扭轉(zhuǎn)自由度。平板斷面豎向和橫向剛度分別為EIz=2.1×106MPa·m4,EIy=1.8×107MPa·m4,扭轉(zhuǎn)剛度GIx= 4.1×105MPa·m4。每延米長度質(zhì)量m=20000 kg/m,質(zhì)量慣矩Im=4.5×106kg·m2/m,空氣的密度ρ=1.248 kg/m3。建模時質(zhì)量矩陣采用集中質(zhì)量矩陣形式。橋面主梁采用BEAM4單元,質(zhì)量慣矩采用MASS21單元模擬,整個模型具有30個橋面梁單元和29個扭轉(zhuǎn)質(zhì)量單元。為便于和文獻(xiàn)[7]對比,結(jié)構(gòu)阻尼設(shè)為零。

    對理論平板氣動導(dǎo)數(shù)進行有理函數(shù)參數(shù)擬合,表1列出了擬合結(jié)果。圖2為擬合曲線和理論曲線的對照。從圖2可以看出,由擬合有理函數(shù)參數(shù)反算顫振導(dǎo)數(shù)和理論顫振導(dǎo)數(shù)之間差別較小,擬合結(jié)果具有較高的精度。

    表1 理想平板有理函數(shù)擬合參數(shù)

    圖2 理想平板顫振導(dǎo)數(shù)擬合值與理論值的比較

    圖3 跨中位移時程(U=135.5 m/s)

    對有理函數(shù)參數(shù)擬合后,給結(jié)構(gòu)施加初始激勵進行動力響應(yīng)分析,逐步提高風(fēng)速,便可得到結(jié)構(gòu)的顫振臨界風(fēng)速。圖3、圖4、圖5分別為風(fēng)速135.5 m/s、136.6 m/s及138 m/s簡支梁跨中點位移響應(yīng)曲線。由圖3可知風(fēng)速為135.5 m/s時跨中點豎向及扭轉(zhuǎn)位移是衰減的。由圖4可知風(fēng)速為136.6 m/s時跨中點豎向及扭轉(zhuǎn)位移具有等幅特性。從圖5可知當(dāng)風(fēng)速為138 m/s時,結(jié)構(gòu)振動位移出現(xiàn)明顯的發(fā)散現(xiàn)象。得到結(jié)構(gòu)的顫振臨界風(fēng)速136.6 m/s,對顫振臨界狀態(tài)位移響應(yīng)時程做頻譜分析得到結(jié)構(gòu)的顫振頻率為f=0.3906 Hz,與文獻(xiàn)[7]中提供的顫振臨界風(fēng)速精確解U=136.3 m/s及顫振頻率f=0.3914 Hz非常接近。

    圖4 跨中位移時程(U=136.6 m/s)

    圖5 跨中位移時程(U=138 m/s)

    3 結(jié) 論

    橋梁時域顫振分析在非線性特性模擬、后顫振形態(tài)以及結(jié)構(gòu)風(fēng)振全過程再現(xiàn)等方面具有頻域方法不可代替的優(yōu)勢,因此時域顫振分析在將來的大跨橋梁抗風(fēng)研究中將會得到越來越廣泛的應(yīng)用。本文的算例表明,基于擬合的氣動自激力,在ANSYS中能方便地實現(xiàn)顫振自激力的施加,計算得到的顫振臨界風(fēng)速和相關(guān)文獻(xiàn)中的報道結(jié)果吻合良好。本文的研究表明,利用ANSYS平臺的二次開發(fā)功能進行橋梁氣動穩(wěn)定性能的時域分析是可行的,這一方法可為橋梁顫振分析或顫抖振理論分析提高效率,降低自主開發(fā)大型軟件過程中的算法錯誤或疏漏所帶來的風(fēng)險。最后,需要指出的是,采用本文時域方法對顫振臨界風(fēng)速求解的關(guān)鍵在于脈沖響應(yīng)函數(shù)參數(shù)的擬合,文中采用最小二乘法。若斷面顫振導(dǎo)數(shù)隨折算風(fēng)速變化非常劇烈,盡管能得到折算風(fēng)速范圍內(nèi)全局最優(yōu)參數(shù),仍可能引起自激力的模擬失真,可嘗試其他參數(shù)擬合方法,繼而利用ANSYS數(shù)值計算平臺進行時域顫振分析。

    [1] 陳政清. 橋梁風(fēng)工程[M]. 北京:人民交通出版社,2005.

    [2] 項海帆. 現(xiàn)代橋梁抗風(fēng)理論與實踐[M]. 北京:人民交通出版社, 2005.

    [3] Scanlan R H, Tomko J J. Airfoil and bridge deck flutter derivatives[J]. Journal of Engineering Mechanics, ASCE, 1971,97(6): 1717-1737.

    [4] Li Q C, Lin Y K. New stochastic theory for bridge stability in turbulent flow[J]. Journal of Engineering Mechanics, ASCE, 1995, 121(1): 102-116.

    [5] Lin Y K, Yang J N. Multimode bridge response to wind excitations [J]. Journal of Engineering Mechanics, ASCE,1983, 109(2): 586-603.

    [6] Scanlan R H. Motion-related body force functions in two-dimensional low-speed flow[J]. Journal of Fluids and Structures, 2000, 14: 49-63.

    [7] Hua X G, Chen Z Q, Ni Y Q, et al. Flutter analysis of long-span bridges using ANSYS[J]. Wind and Structures, 2007, 10(1): 61-82.

    [8] 華旭剛,陳政清,祝志文.一種在ANSYS中實現(xiàn)顫振時程分析的方法[J].中國公路學(xué)報,2002, 15(4): 32-34.

    [9] 張志田. 大跨度橋梁非線性抖振及其對抗風(fēng)穩(wěn)定性影響的研究[D]. 上海:同濟大學(xué),2004.

    猜你喜歡
    脈沖響應(yīng)表達(dá)式時域
    基于重復(fù)脈沖響應(yīng)的發(fā)電機轉(zhuǎn)子繞組匝間短路檢測技術(shù)的研究與應(yīng)用
    一個混合核Hilbert型積分不等式及其算子范數(shù)表達(dá)式
    表達(dá)式轉(zhuǎn)換及求值探析
    淺析C語言運算符及表達(dá)式的教學(xué)誤區(qū)
    基于時域信號的三電平逆變器復(fù)合故障診斷
    基于極大似然準(zhǔn)則與滾動時域估計的自適應(yīng)UKF算法
    基于時域逆濾波的寬帶脈沖聲生成技術(shù)
    脈沖響應(yīng)函數(shù)下的我國貨幣需求變動與決定
    基于有限元素法的室內(nèi)脈沖響應(yīng)的仿真
    電大理工(2015年3期)2015-12-03 11:34:12
    基于時域波形特征的輸電線雷擊識別
    電測與儀表(2015年2期)2015-04-09 11:28:50
    精品少妇一区二区三区视频日本电影| 天天添夜夜摸| 久久精品国产综合久久久| 亚洲aⅴ乱码一区二区在线播放 | 在线看三级毛片| 亚洲av五月六月丁香网| 久久性视频一级片| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美在线二视频| 欧美激情高清一区二区三区| 成人特级黄色片久久久久久久| 三级毛片av免费| av免费在线观看网站| av福利片在线| 男人舔女人下体高潮全视频| 欧美乱色亚洲激情| 欧美日韩福利视频一区二区| 男女下面进入的视频免费午夜 | 久久精品国产亚洲av高清一级| 不卡av一区二区三区| 成人手机av| 欧美激情 高清一区二区三区| 亚洲欧美精品综合一区二区三区| 夜夜夜夜夜久久久久| 国产三级黄色录像| 国产亚洲精品av在线| 亚洲精品久久国产高清桃花| 在线视频色国产色| 成人永久免费在线观看视频| 1024手机看黄色片| www日本黄色视频网| 黄色毛片三级朝国网站| 欧美人与性动交α欧美精品济南到| 日韩欧美三级三区| 欧美日韩黄片免| 国产av一区在线观看免费| 亚洲aⅴ乱码一区二区在线播放 | 国产成年人精品一区二区| 国产视频一区二区在线看| 亚洲免费av在线视频| 久久天堂一区二区三区四区| 欧美性长视频在线观看| 亚洲avbb在线观看| 国产精品久久视频播放| www日本黄色视频网| 熟女电影av网| 日韩大尺度精品在线看网址| 真人一进一出gif抽搐免费| 国产精品,欧美在线| 国产99白浆流出| 最近最新中文字幕大全电影3 | 国产成人欧美| 中文在线观看免费www的网站 | 熟妇人妻久久中文字幕3abv| 一本一本综合久久| 亚洲精品久久成人aⅴ小说| 久久久久国产精品人妻aⅴ院| 非洲黑人性xxxx精品又粗又长| 亚洲无线在线观看| 国产在线观看jvid| 少妇熟女aⅴ在线视频| 特大巨黑吊av在线直播 | 成人手机av| 男人舔女人下体高潮全视频| 精品国产美女av久久久久小说| 成人永久免费在线观看视频| av片东京热男人的天堂| 国产一区二区三区在线臀色熟女| 岛国在线观看网站| 亚洲一区高清亚洲精品| 国产一区二区激情短视频| 亚洲专区中文字幕在线| 美国免费a级毛片| 色综合欧美亚洲国产小说| 久久精品国产亚洲av高清一级| 精品第一国产精品| 久久婷婷人人爽人人干人人爱| 午夜精品在线福利| 欧美最黄视频在线播放免费| 久久精品国产亚洲av高清一级| 12—13女人毛片做爰片一| 日韩大码丰满熟妇| 国产精品久久久久久精品电影 | 波多野结衣巨乳人妻| 久久精品91无色码中文字幕| 中文字幕人成人乱码亚洲影| 午夜久久久在线观看| 搡老妇女老女人老熟妇| 伊人久久大香线蕉亚洲五| 高潮久久久久久久久久久不卡| 久9热在线精品视频| 国产私拍福利视频在线观看| 好男人电影高清在线观看| 日日夜夜操网爽| 国产精品二区激情视频| 国产精品国产高清国产av| 法律面前人人平等表现在哪些方面| 久久久久久久精品吃奶| 此物有八面人人有两片| 国产精品亚洲美女久久久| 久久精品国产99精品国产亚洲性色| 免费女性裸体啪啪无遮挡网站| 亚洲国产毛片av蜜桃av| 日本熟妇午夜| 欧美另类亚洲清纯唯美| 久久精品成人免费网站| 成人亚洲精品一区在线观看| 他把我摸到了高潮在线观看| 日本成人三级电影网站| 精品国产乱子伦一区二区三区| 国产成人av教育| 露出奶头的视频| 日韩欧美一区视频在线观看| 国产亚洲av高清不卡| 国产成人精品久久二区二区免费| 人妻久久中文字幕网| 亚洲七黄色美女视频| 久久精品国产亚洲av香蕉五月| 欧美成人午夜精品| 伦理电影免费视频| bbb黄色大片| 欧美在线一区亚洲| 久久久久久大精品| 亚洲激情在线av| 欧美又色又爽又黄视频| 亚洲人成网站在线播放欧美日韩| 午夜日韩欧美国产| 亚洲第一电影网av| 很黄的视频免费| 这个男人来自地球电影免费观看| 国产精品自产拍在线观看55亚洲| 亚洲激情在线av| 国产又爽黄色视频| 神马国产精品三级电影在线观看 | 欧美性猛交黑人性爽| 国产三级在线视频| 日本黄色视频三级网站网址| 美女 人体艺术 gogo| 18禁裸乳无遮挡免费网站照片 | 色哟哟哟哟哟哟| 无遮挡黄片免费观看| 自线自在国产av| 中文字幕高清在线视频| 男女那种视频在线观看| 日韩大尺度精品在线看网址| 啦啦啦韩国在线观看视频| 男女那种视频在线观看| 久久人妻av系列| 久久人妻av系列| 亚洲欧美日韩高清在线视频| xxx96com| 精品久久久久久久久久免费视频| 国产aⅴ精品一区二区三区波| 日韩欧美国产在线观看| 三级毛片av免费| 色精品久久人妻99蜜桃| 国产黄a三级三级三级人| 国产成人精品久久二区二区免费| 欧美国产日韩亚洲一区| 成人亚洲精品av一区二区| 国产97色在线日韩免费| 人人妻人人澡人人看| 亚洲一码二码三码区别大吗| 日本精品一区二区三区蜜桃| 身体一侧抽搐| 久久国产精品影院| 757午夜福利合集在线观看| 麻豆一二三区av精品| 亚洲精品国产一区二区精华液| 亚洲欧美激情综合另类| 国产精品美女特级片免费视频播放器 | 欧美乱色亚洲激情| 久久人妻av系列| 国产精品久久电影中文字幕| 欧美日韩福利视频一区二区| 一边摸一边抽搐一进一小说| 大型av网站在线播放| 国产私拍福利视频在线观看| 一a级毛片在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 婷婷精品国产亚洲av在线| 少妇粗大呻吟视频| 女同久久另类99精品国产91| 激情在线观看视频在线高清| 91九色精品人成在线观看| 九色国产91popny在线| 亚洲人成网站高清观看| 黄片大片在线免费观看| 天堂动漫精品| 老熟妇乱子伦视频在线观看| 久久 成人 亚洲| 人人妻人人澡人人看| 俄罗斯特黄特色一大片| 2021天堂中文幕一二区在线观 | 久久久久久久久中文| 国产三级黄色录像| 男人舔奶头视频| 超碰成人久久| 国产aⅴ精品一区二区三区波| 久久天堂一区二区三区四区| 国产又黄又爽又无遮挡在线| 国产成人一区二区三区免费视频网站| 免费观看人在逋| 日本 av在线| 亚洲精品在线美女| 色av中文字幕| 精品久久久久久久末码| 天堂动漫精品| 狠狠狠狠99中文字幕| 草草在线视频免费看| 亚洲精品中文字幕在线视频| 欧美精品啪啪一区二区三区| 男男h啪啪无遮挡| 成人国语在线视频| 欧美日韩一级在线毛片| 最近在线观看免费完整版| 久久精品国产亚洲av香蕉五月| 国产又爽黄色视频| 久久久久亚洲av毛片大全| 变态另类成人亚洲欧美熟女| 亚洲午夜精品一区,二区,三区| 中文字幕最新亚洲高清| 国产成人精品久久二区二区免费| 淫秽高清视频在线观看| 久久天堂一区二区三区四区| or卡值多少钱| 亚洲精品久久国产高清桃花| bbb黄色大片| 精品午夜福利视频在线观看一区| 黄色a级毛片大全视频| 中文资源天堂在线| 亚洲成人免费电影在线观看| 18禁观看日本| 人人妻人人澡欧美一区二区| 国产精品,欧美在线| 国产精品av久久久久免费| 最近最新免费中文字幕在线| 最近在线观看免费完整版| 国产v大片淫在线免费观看| 久久久久久久久免费视频了| 中文字幕精品免费在线观看视频| 熟女电影av网| 男女午夜视频在线观看| 亚洲,欧美精品.| 国内毛片毛片毛片毛片毛片| 日韩高清综合在线| xxxwww97欧美| xxx96com| 每晚都被弄得嗷嗷叫到高潮| 午夜福利一区二区在线看| 搡老妇女老女人老熟妇| 欧洲精品卡2卡3卡4卡5卡区| 中文资源天堂在线| 黄色丝袜av网址大全| 国产精品1区2区在线观看.| 岛国视频午夜一区免费看| 真人一进一出gif抽搐免费| 国产精品永久免费网站| 天堂影院成人在线观看| 日韩免费av在线播放| 两个人视频免费观看高清| 国产aⅴ精品一区二区三区波| 波多野结衣高清无吗| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩亚洲国产一区二区在线观看| 成人特级黄色片久久久久久久| videosex国产| 国产av一区在线观看免费| ponron亚洲| 欧美性猛交╳xxx乱大交人| 国产激情久久老熟女| 精品国产超薄肉色丝袜足j| 91av网站免费观看| 国产日本99.免费观看| netflix在线观看网站| 日韩 欧美 亚洲 中文字幕| 淫秽高清视频在线观看| 欧美乱妇无乱码| 久久精品亚洲精品国产色婷小说| 国产精品久久久av美女十八| 国产精品久久视频播放| 亚洲中文日韩欧美视频| 无遮挡黄片免费观看| 欧美成人一区二区免费高清观看 | 韩国精品一区二区三区| 视频区欧美日本亚洲| 可以免费在线观看a视频的电影网站| 欧美激情久久久久久爽电影| 99热只有精品国产| 男人的好看免费观看在线视频 | 国产精品98久久久久久宅男小说| www.自偷自拍.com| АⅤ资源中文在线天堂| 亚洲一码二码三码区别大吗| 欧美在线一区亚洲| 一区二区三区精品91| 99精品久久久久人妻精品| 无人区码免费观看不卡| 国内久久婷婷六月综合欲色啪| 亚洲av成人一区二区三| 身体一侧抽搐| 亚洲第一青青草原| 久久精品国产亚洲av高清一级| 欧美人与性动交α欧美精品济南到| 99热6这里只有精品| 国产成人av激情在线播放| 777久久人妻少妇嫩草av网站| 99国产极品粉嫩在线观看| 国产精品二区激情视频| 亚洲国产欧美一区二区综合| 亚洲第一欧美日韩一区二区三区| 免费一级毛片在线播放高清视频| 日本黄色视频三级网站网址| 黄色女人牲交| 亚洲av日韩精品久久久久久密| 香蕉丝袜av| 中文字幕人妻丝袜一区二区| 亚洲av片天天在线观看| 国产亚洲精品久久久久5区| 日本成人三级电影网站| 天堂动漫精品| 一个人免费在线观看的高清视频| 久9热在线精品视频| 日本一本二区三区精品| 欧美精品亚洲一区二区| 久久久久亚洲av毛片大全| 波多野结衣高清作品| 啪啪无遮挡十八禁网站| 无人区码免费观看不卡| 两个人看的免费小视频| 欧美日韩精品网址| 午夜免费成人在线视频| 成人免费观看视频高清| 成人手机av| 丰满人妻熟妇乱又伦精品不卡| 天天一区二区日本电影三级| 国产视频一区二区在线看| 国产野战对白在线观看| 天堂动漫精品| 18美女黄网站色大片免费观看| 成人国语在线视频| 观看免费一级毛片| 欧美日韩瑟瑟在线播放| 男人舔女人的私密视频| 丰满的人妻完整版| 亚洲七黄色美女视频| 久9热在线精品视频| 午夜久久久在线观看| 久久精品亚洲精品国产色婷小说| 每晚都被弄得嗷嗷叫到高潮| avwww免费| 日韩视频一区二区在线观看| 1024手机看黄色片| 伦理电影免费视频| 一本大道久久a久久精品| 日本五十路高清| 欧美激情极品国产一区二区三区| 真人一进一出gif抽搐免费| xxx96com| 久久精品91蜜桃| 欧美午夜高清在线| 国产精品久久视频播放| 精品久久久久久久久久久久久 | 动漫黄色视频在线观看| 亚洲成av片中文字幕在线观看| 久久久国产成人免费| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻1区二区| 亚洲狠狠婷婷综合久久图片| 人人妻人人澡人人看| 怎么达到女性高潮| 国产一级毛片七仙女欲春2 | 国产精华一区二区三区| 最新美女视频免费是黄的| 精品少妇一区二区三区视频日本电影| 丰满的人妻完整版| 又大又爽又粗| 九色国产91popny在线| 黄网站色视频无遮挡免费观看| 88av欧美| 女警被强在线播放| 国产熟女xx| 午夜影院日韩av| 久久人妻福利社区极品人妻图片| 久久久国产成人精品二区| 国产亚洲精品一区二区www| 久久久国产欧美日韩av| 大型av网站在线播放| 成人三级做爰电影| 婷婷精品国产亚洲av| 免费av毛片视频| 久久性视频一级片| 国产激情偷乱视频一区二区| 精品电影一区二区在线| or卡值多少钱| 狂野欧美激情性xxxx| 在线观看66精品国产| 他把我摸到了高潮在线观看| 久久午夜亚洲精品久久| 国产又黄又爽又无遮挡在线| 色哟哟哟哟哟哟| 亚洲精品美女久久久久99蜜臀| 非洲黑人性xxxx精品又粗又长| 好看av亚洲va欧美ⅴa在| 美女大奶头视频| 日本a在线网址| 精品欧美一区二区三区在线| 国产熟女xx| 在线视频色国产色| 99riav亚洲国产免费| 丁香六月欧美| 无遮挡黄片免费观看| 国产精品 国内视频| 精品久久久久久久久久免费视频| 99久久久亚洲精品蜜臀av| 国产精品自产拍在线观看55亚洲| 成人三级做爰电影| 国产黄片美女视频| 十八禁网站免费在线| 日韩有码中文字幕| 天天一区二区日本电影三级| 久久精品国产清高在天天线| 后天国语完整版免费观看| 色播在线永久视频| 美女大奶头视频| 亚洲精华国产精华精| 日日夜夜操网爽| 久久久国产欧美日韩av| 午夜福利在线观看吧| 精品久久久久久久人妻蜜臀av| 欧美av亚洲av综合av国产av| 色精品久久人妻99蜜桃| 啦啦啦韩国在线观看视频| 国产亚洲精品第一综合不卡| 视频在线观看一区二区三区| 十八禁人妻一区二区| 欧美丝袜亚洲另类 | 亚洲久久久国产精品| 在线看三级毛片| 女同久久另类99精品国产91| 女生性感内裤真人,穿戴方法视频| 91麻豆av在线| 亚洲中文字幕日韩| 日韩视频一区二区在线观看| 国产亚洲欧美精品永久| 精品国产亚洲在线| 国内毛片毛片毛片毛片毛片| 亚洲天堂国产精品一区在线| 亚洲av成人不卡在线观看播放网| 一夜夜www| 久久久久久久久免费视频了| 亚洲真实伦在线观看| 人成视频在线观看免费观看| 日日爽夜夜爽网站| av在线播放免费不卡| 成年人黄色毛片网站| 国产成人精品久久二区二区免费| 欧美在线黄色| 欧美成狂野欧美在线观看| 国产亚洲精品一区二区www| 他把我摸到了高潮在线观看| 亚洲人成伊人成综合网2020| 精品不卡国产一区二区三区| 国产精品九九99| 亚洲精品国产一区二区精华液| 欧美亚洲日本最大视频资源| 免费看美女性在线毛片视频| 一级毛片女人18水好多| 久久久国产欧美日韩av| 久久性视频一级片| 欧美性长视频在线观看| 久久精品影院6| 国产激情久久老熟女| 日韩免费av在线播放| 久久亚洲真实| 男男h啪啪无遮挡| 久久精品国产亚洲av香蕉五月| 久久久精品国产亚洲av高清涩受| 777久久人妻少妇嫩草av网站| 亚洲熟妇熟女久久| 精品久久久久久成人av| 大型黄色视频在线免费观看| 亚洲av熟女| 日日夜夜操网爽| 国产精品久久久人人做人人爽| 老司机午夜十八禁免费视频| 国产精品亚洲一级av第二区| 亚洲九九香蕉| 国产又黄又爽又无遮挡在线| 亚洲精品在线观看二区| 大香蕉久久成人网| 日韩国内少妇激情av| 亚洲精品美女久久av网站| 99热这里只有精品一区 | 欧美激情 高清一区二区三区| 在线观看日韩欧美| 欧美日韩瑟瑟在线播放| 欧美激情极品国产一区二区三区| 日本在线视频免费播放| 日本五十路高清| 久久性视频一级片| 老司机深夜福利视频在线观看| www.999成人在线观看| 精品欧美一区二区三区在线| 每晚都被弄得嗷嗷叫到高潮| 日韩av在线大香蕉| 久久国产亚洲av麻豆专区| 黄片小视频在线播放| 欧美成人免费av一区二区三区| 国产成人欧美| 最近最新免费中文字幕在线| 国产伦在线观看视频一区| 19禁男女啪啪无遮挡网站| 欧美国产日韩亚洲一区| 男女午夜视频在线观看| 制服诱惑二区| 欧美不卡视频在线免费观看 | 国产男靠女视频免费网站| 一本久久中文字幕| 欧美日韩瑟瑟在线播放| or卡值多少钱| 久久青草综合色| 两个人免费观看高清视频| 十分钟在线观看高清视频www| 午夜视频精品福利| 国产精品自产拍在线观看55亚洲| 中文字幕久久专区| 久久中文看片网| 亚洲av电影不卡..在线观看| 日本五十路高清| 国产精品美女特级片免费视频播放器 | 久热爱精品视频在线9| 日韩视频一区二区在线观看| 精品国内亚洲2022精品成人| 欧美午夜高清在线| 精品久久久久久成人av| 波多野结衣高清无吗| 国产熟女午夜一区二区三区| 变态另类成人亚洲欧美熟女| 最近最新中文字幕大全免费视频| 99久久99久久久精品蜜桃| 欧美激情极品国产一区二区三区| 九色国产91popny在线| 一a级毛片在线观看| xxxwww97欧美| 好男人电影高清在线观看| 久9热在线精品视频| 亚洲精品在线美女| 又大又爽又粗| 成人国产综合亚洲| 久久中文字幕一级| 99国产精品99久久久久| 一二三四社区在线视频社区8| 国产人伦9x9x在线观看| 亚洲国产看品久久| 最新在线观看一区二区三区| 大型黄色视频在线免费观看| 99精品久久久久人妻精品| 久久久久久免费高清国产稀缺| 91国产中文字幕| 亚洲 欧美 日韩 在线 免费| 18美女黄网站色大片免费观看| 手机成人av网站| 满18在线观看网站| 欧美精品啪啪一区二区三区| 免费女性裸体啪啪无遮挡网站| 热re99久久国产66热| 天堂√8在线中文| 国内精品久久久久精免费| 国产成年人精品一区二区| 午夜久久久在线观看| 久久国产乱子伦精品免费另类| 久久久精品国产亚洲av高清涩受| 中文亚洲av片在线观看爽| 久久亚洲真实| 国产精品一区二区精品视频观看| 日韩av在线大香蕉| 老司机深夜福利视频在线观看| 久久中文字幕一级| 99久久99久久久精品蜜桃| 免费高清视频大片| 亚洲avbb在线观看| 一a级毛片在线观看| 亚洲av电影不卡..在线观看| 久久久国产欧美日韩av| 首页视频小说图片口味搜索| 欧美又色又爽又黄视频| 一级作爱视频免费观看| 99久久久亚洲精品蜜臀av| 久9热在线精品视频| 悠悠久久av| 国产极品粉嫩免费观看在线| 男人舔女人的私密视频| 国产极品粉嫩免费观看在线| 亚洲全国av大片| 国产极品粉嫩免费观看在线| 俺也久久电影网| 久久久久久人人人人人| 精品无人区乱码1区二区| 欧美在线黄色| 日本一本二区三区精品| 欧美激情极品国产一区二区三区| av电影中文网址| 18美女黄网站色大片免费观看| 亚洲中文字幕一区二区三区有码在线看 | 正在播放国产对白刺激| 亚洲熟妇熟女久久| 久久精品亚洲精品国产色婷小说| av超薄肉色丝袜交足视频| 老司机靠b影院| 久久久国产成人精品二区| 亚洲片人在线观看| 一级毛片女人18水好多| 国产亚洲精品久久久久5区| 亚洲专区国产一区二区|