• <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
    国产亚洲精品一区二区www| 亚洲乱码一区二区免费版| 成人国产综合亚洲| 国产爱豆传媒在线观看| 一区二区三区高清视频在线| 午夜a级毛片| 亚洲精品色激情综合| 亚洲无线在线观看| 免费看美女性在线毛片视频| 99国产精品一区二区蜜桃av| 亚洲精品影视一区二区三区av| 最近在线观看免费完整版| x7x7x7水蜜桃| 色综合欧美亚洲国产小说| 国产av一区在线观看免费| 国产高清三级在线| 成年女人永久免费观看视频| 亚洲欧美一区二区三区黑人| 午夜福利视频1000在线观看| 欧美成人免费av一区二区三区| 国产成人啪精品午夜网站| 亚洲人成网站高清观看| 国产成+人综合+亚洲专区| 欧美日韩瑟瑟在线播放| 亚洲精品在线观看二区| 亚洲欧美日韩无卡精品| 波多野结衣高清作品| 成人特级av手机在线观看| 51午夜福利影视在线观看| 午夜福利在线观看吧| 一本久久中文字幕| 亚洲国产色片| 国产精品久久电影中文字幕| 国产av麻豆久久久久久久| 欧美乱色亚洲激情| 亚洲成av人片在线播放无| 国语自产精品视频在线第100页| 成年女人看的毛片在线观看| 久久精品人妻少妇| 18禁国产床啪视频网站| 狠狠狠狠99中文字幕| 日本一二三区视频观看| 99久久精品国产亚洲精品| 亚洲男人的天堂狠狠| 亚洲av成人精品一区久久| 18+在线观看网站| 欧美日韩精品网址| 97碰自拍视频| 国产精品,欧美在线| 国产黄a三级三级三级人| 日本免费a在线| 午夜老司机福利剧场| 熟女人妻精品中文字幕| 99久国产av精品| 怎么达到女性高潮| 女同久久另类99精品国产91| 久久欧美精品欧美久久欧美| 九色成人免费人妻av| 一个人免费在线观看电影| 国内精品久久久久久久电影| 蜜桃亚洲精品一区二区三区| 女人被狂操c到高潮| 亚洲国产精品久久男人天堂| 日本黄大片高清| 88av欧美| 亚洲国产色片| 欧美乱色亚洲激情| 色综合亚洲欧美另类图片| 国产精品 国内视频| 蜜桃久久精品国产亚洲av| 国产精品99久久99久久久不卡| 少妇丰满av| 69av精品久久久久久| 欧美日韩亚洲国产一区二区在线观看| 又紧又爽又黄一区二区| 波多野结衣高清作品| 麻豆成人午夜福利视频| a级一级毛片免费在线观看| ponron亚洲| 国产黄a三级三级三级人| 亚洲五月天丁香| bbb黄色大片| 久久精品国产99精品国产亚洲性色| 狂野欧美白嫩少妇大欣赏| 欧美极品一区二区三区四区| 两个人看的免费小视频| 窝窝影院91人妻| 欧美日韩瑟瑟在线播放| 美女大奶头视频| 男女床上黄色一级片免费看| 成人无遮挡网站| 国产中年淑女户外野战色| 亚洲第一电影网av| 手机成人av网站| 亚洲精品亚洲一区二区| 床上黄色一级片| 久久久久国内视频| 内地一区二区视频在线| 无限看片的www在线观看| 91久久精品电影网| 美女高潮喷水抽搐中文字幕| 69av精品久久久久久| 在线看三级毛片| 亚洲av成人精品一区久久| 搡老熟女国产l中国老女人| 黄色日韩在线| 少妇丰满av| 99国产综合亚洲精品| 给我免费播放毛片高清在线观看| 尤物成人国产欧美一区二区三区| 一级毛片高清免费大全| 日本黄大片高清| 99久久精品热视频| 久久6这里有精品| 在线免费观看不下载黄p国产 | 国产伦在线观看视频一区| 久久香蕉精品热| 日韩成人在线观看一区二区三区| 国产精品日韩av在线免费观看| 精品久久久久久久末码| 欧美中文综合在线视频| 国产精品国产高清国产av| 日韩欧美国产一区二区入口| 哪里可以看免费的av片| 啪啪无遮挡十八禁网站| 亚洲av不卡在线观看| 可以在线观看毛片的网站| 波野结衣二区三区在线 | 九九热线精品视视频播放| www.www免费av| 午夜福利成人在线免费观看| 国产欧美日韩一区二区三| 18禁在线播放成人免费| 老司机福利观看| 搞女人的毛片| 日韩欧美精品v在线| 有码 亚洲区| 免费观看的影片在线观看| 国产精品久久电影中文字幕| 女同久久另类99精品国产91| 1000部很黄的大片| 黄色女人牲交| 一边摸一边抽搐一进一小说| 国产69精品久久久久777片| 精品国内亚洲2022精品成人| 国产亚洲精品久久久久久毛片| 五月伊人婷婷丁香| 老汉色av国产亚洲站长工具| 国产美女午夜福利| 啦啦啦观看免费观看视频高清| 亚洲专区国产一区二区| 欧美激情久久久久久爽电影| 麻豆国产av国片精品| 高清在线国产一区| 欧美日韩综合久久久久久 | 91久久精品电影网| 丰满人妻一区二区三区视频av | 无限看片的www在线观看| 国产三级中文精品| 亚洲熟妇熟女久久| 一进一出抽搐动态| 欧美日韩国产亚洲二区| 亚洲精品亚洲一区二区| 久久久国产精品麻豆| 禁无遮挡网站| 中文字幕人妻熟人妻熟丝袜美 | 婷婷精品国产亚洲av在线| 久久亚洲精品不卡| 国产国拍精品亚洲av在线观看 | 最近最新中文字幕大全电影3| 香蕉久久夜色| 久久国产乱子伦精品免费另类| 国产精品自产拍在线观看55亚洲| e午夜精品久久久久久久| 真人做人爱边吃奶动态| 男人舔奶头视频| 欧美一区二区国产精品久久精品| 国产成人av教育| 国产激情偷乱视频一区二区| 又紧又爽又黄一区二区| 91麻豆精品激情在线观看国产| 一a级毛片在线观看| 制服丝袜大香蕉在线| 很黄的视频免费| av专区在线播放| 国内精品美女久久久久久| 美女高潮喷水抽搐中文字幕| 一个人免费在线观看电影| x7x7x7水蜜桃| 欧美日本视频| 久久欧美精品欧美久久欧美| 男女之事视频高清在线观看| 精品99又大又爽又粗少妇毛片 | or卡值多少钱| 黄色片一级片一级黄色片| 亚洲av中文字字幕乱码综合| 成人国产一区最新在线观看| 国产精品女同一区二区软件 | 一区二区三区高清视频在线| 国产一区二区亚洲精品在线观看| bbb黄色大片| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 国产成人系列免费观看| 国产高清有码在线观看视频| 狂野欧美激情性xxxx| 国产亚洲精品一区二区www| 亚洲精品在线观看二区| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清专用| 啦啦啦观看免费观看视频高清| 亚洲精品在线美女| 19禁男女啪啪无遮挡网站| 亚洲av一区综合| 夜夜夜夜夜久久久久| 无限看片的www在线观看| 亚洲国产精品999在线| 日本免费a在线| 99久久综合精品五月天人人| 99久久精品一区二区三区| 午夜福利高清视频| 黄色女人牲交| 国产伦人伦偷精品视频| 亚洲人成网站在线播| 手机成人av网站| 久久国产乱子伦精品免费另类| 欧美在线一区亚洲| 中文字幕人妻熟人妻熟丝袜美 | 欧美一区二区亚洲| aaaaa片日本免费| 一本精品99久久精品77| 精品一区二区三区视频在线观看免费| 亚洲国产欧洲综合997久久,| 可以在线观看的亚洲视频| 久久久精品大字幕| 日韩人妻高清精品专区| 久久久精品欧美日韩精品| 很黄的视频免费| 亚洲国产日韩欧美精品在线观看 | aaaaa片日本免费| 看黄色毛片网站| 欧美3d第一页| 久久久久国产精品人妻aⅴ院| 久久久久久人人人人人| 乱人视频在线观看| 一级黄片播放器| 欧美高清成人免费视频www| 国产精品,欧美在线| 国产欧美日韩一区二区精品| 亚洲国产欧美网| 人人妻人人澡欧美一区二区| 18禁在线播放成人免费| 法律面前人人平等表现在哪些方面| 欧美区成人在线视频| 成人国产一区最新在线观看| 99国产精品一区二区蜜桃av| 国产亚洲精品一区二区www| www.999成人在线观看| 国产高潮美女av| 天堂影院成人在线观看| 人妻夜夜爽99麻豆av| 国产精品电影一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久成人av| 一卡2卡三卡四卡精品乱码亚洲| 国产麻豆成人av免费视频| 真人一进一出gif抽搐免费| 亚洲欧美激情综合另类| 老司机福利观看| 99热6这里只有精品| 操出白浆在线播放| 蜜桃亚洲精品一区二区三区| 桃色一区二区三区在线观看| 99热这里只有精品一区| 美女 人体艺术 gogo| 日本在线视频免费播放| 国产久久久一区二区三区| 尤物成人国产欧美一区二区三区| 高清在线国产一区| 国产又黄又爽又无遮挡在线| 国产精品电影一区二区三区| 天天一区二区日本电影三级| 亚洲国产高清在线一区二区三| 欧美一区二区精品小视频在线| 久久欧美精品欧美久久欧美| 五月伊人婷婷丁香| 日本 av在线| 香蕉av资源在线| 欧美午夜高清在线| 亚洲av日韩精品久久久久久密| 淫秽高清视频在线观看| 超碰av人人做人人爽久久 | 免费av毛片视频| 他把我摸到了高潮在线观看| 99精品欧美一区二区三区四区| 一二三四社区在线视频社区8| 少妇的逼好多水| 69人妻影院| 欧美色视频一区免费| 国产成人a区在线观看| 看免费av毛片| 熟妇人妻久久中文字幕3abv| 波野结衣二区三区在线 | 久久精品人妻少妇| 热99在线观看视频| 久久久国产精品麻豆| 俄罗斯特黄特色一大片| 亚洲精华国产精华精| 非洲黑人性xxxx精品又粗又长| 亚洲avbb在线观看| 国产av麻豆久久久久久久| 无限看片的www在线观看| 少妇高潮的动态图| 三级国产精品欧美在线观看| 亚洲成人精品中文字幕电影| 午夜福利欧美成人| 精品国产美女av久久久久小说| 亚洲国产欧美人成| 天天躁日日操中文字幕| 久久久久性生活片| 免费在线观看日本一区| 亚洲欧美日韩东京热| 精品国产亚洲在线| 欧美一区二区精品小视频在线| 18禁国产床啪视频网站| www.www免费av| 蜜桃久久精品国产亚洲av| 国产精品98久久久久久宅男小说| 亚洲美女黄片视频| 男女做爰动态图高潮gif福利片| 黄片大片在线免费观看| 亚洲在线自拍视频| 亚洲电影在线观看av| 中文字幕人妻丝袜一区二区| 国产欧美日韩一区二区精品| 国产成人a区在线观看| 国产精品99久久99久久久不卡| 婷婷六月久久综合丁香| 在线观看午夜福利视频| 国产aⅴ精品一区二区三区波| 亚洲 欧美 日韩 在线 免费| 国产私拍福利视频在线观看| 欧美成人一区二区免费高清观看| 丰满人妻熟妇乱又伦精品不卡| 一级黄色大片毛片| 国产高清视频在线观看网站| 桃红色精品国产亚洲av| 亚洲av成人精品一区久久| 国产高潮美女av| 欧美午夜高清在线| 99在线视频只有这里精品首页| 可以在线观看毛片的网站| 成人国产一区最新在线观看| 国产一区二区激情短视频| 操出白浆在线播放| 麻豆一二三区av精品| 亚洲乱码一区二区免费版| 天美传媒精品一区二区| 法律面前人人平等表现在哪些方面| 在线a可以看的网站| 国产精品一及| 久久久国产成人免费| 9191精品国产免费久久| 亚洲自拍偷在线| av视频在线观看入口| 欧美大码av| 伊人久久大香线蕉亚洲五| av在线天堂中文字幕| 亚洲男人的天堂狠狠| 欧美日韩精品网址| 国产精品一及| bbb黄色大片| 久久精品亚洲精品国产色婷小说| 一a级毛片在线观看| 久久久久国内视频| 成人一区二区视频在线观看| 男女午夜视频在线观看| 中国美女看黄片| 99久久精品一区二区三区| 少妇丰满av| 国产亚洲精品久久久com| 18+在线观看网站| 在线视频色国产色| 母亲3免费完整高清在线观看| 一级毛片女人18水好多| 日本与韩国留学比较| 免费人成视频x8x8入口观看| 观看免费一级毛片| 此物有八面人人有两片| 18禁美女被吸乳视频| 最近最新中文字幕大全电影3| 精品电影一区二区在线| 女人十人毛片免费观看3o分钟| 免费看日本二区| 在线a可以看的网站| 综合色av麻豆| 91麻豆av在线| 不卡一级毛片| 极品教师在线免费播放| 国产久久久一区二区三区| 欧美日韩黄片免| 国产伦人伦偷精品视频| 男女那种视频在线观看| 日本五十路高清| 最近最新免费中文字幕在线| 国内精品久久久久精免费| 亚洲欧美日韩无卡精品| 男插女下体视频免费在线播放| 国产色婷婷99| 精品久久久久久,| 桃红色精品国产亚洲av| 欧美日本视频| 精品久久久久久,| 日本一二三区视频观看| 91九色精品人成在线观看| 男人舔女人下体高潮全视频| 最近在线观看免费完整版| 欧美日本视频| 丰满人妻一区二区三区视频av | 欧美国产日韩亚洲一区| 亚洲成人久久性| 日日干狠狠操夜夜爽| 啦啦啦观看免费观看视频高清| 国产三级在线视频| 色精品久久人妻99蜜桃| 人人妻人人看人人澡| 久久久成人免费电影| 精品一区二区三区视频在线 | 女生性感内裤真人,穿戴方法视频| 91麻豆av在线| 每晚都被弄得嗷嗷叫到高潮| 51午夜福利影视在线观看| 中文字幕高清在线视频| 精品国产美女av久久久久小说| 一级黄片播放器| 女同久久另类99精品国产91| 欧美性感艳星| 九色国产91popny在线| 婷婷精品国产亚洲av在线| 久久久久国内视频| 韩国av一区二区三区四区| 99热只有精品国产| 国产单亲对白刺激| 国产免费一级a男人的天堂| 黄色视频,在线免费观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成人亚洲精品av一区二区| 每晚都被弄得嗷嗷叫到高潮| 亚洲黑人精品在线| 亚洲午夜理论影院| 日韩欧美三级三区| 91麻豆av在线| 又粗又爽又猛毛片免费看| 久久天躁狠狠躁夜夜2o2o| 欧美日本视频| 国产av不卡久久| 日韩欧美国产在线观看| 亚洲美女黄片视频| 1024手机看黄色片| 网址你懂的国产日韩在线| 成人特级av手机在线观看| 亚洲男人的天堂狠狠| 波多野结衣高清无吗| tocl精华| 国产乱人伦免费视频| 俄罗斯特黄特色一大片| 天堂网av新在线| 最新中文字幕久久久久| 亚洲成人久久爱视频| 午夜精品一区二区三区免费看| 看片在线看免费视频| 黄片小视频在线播放| 欧美成人免费av一区二区三区| 成人国产一区最新在线观看| 男人和女人高潮做爰伦理| av福利片在线观看| 两个人视频免费观看高清| 欧美乱码精品一区二区三区| 国产黄a三级三级三级人| av黄色大香蕉| 日韩 欧美 亚洲 中文字幕| 非洲黑人性xxxx精品又粗又长| 国产免费一级a男人的天堂| tocl精华| 12—13女人毛片做爰片一| 午夜精品久久久久久毛片777| 夜夜看夜夜爽夜夜摸| 国产成人影院久久av| 全区人妻精品视频| 一个人观看的视频www高清免费观看| 搡老熟女国产l中国老女人| 国产一区二区激情短视频| 国产欧美日韩精品一区二区| 美女高潮喷水抽搐中文字幕| 国产aⅴ精品一区二区三区波| 国产主播在线观看一区二区| 国产单亲对白刺激| 村上凉子中文字幕在线| 无人区码免费观看不卡| 亚洲成人久久性| 男插女下体视频免费在线播放| 又黄又爽又免费观看的视频| 18禁美女被吸乳视频| 久久久久免费精品人妻一区二区| 国产av不卡久久| 脱女人内裤的视频| 久久精品国产清高在天天线| 熟女少妇亚洲综合色aaa.| av国产免费在线观看| 免费在线观看成人毛片| 一区二区三区高清视频在线| 真人一进一出gif抽搐免费| 欧美区成人在线视频| 国产精品久久久久久人妻精品电影| www日本在线高清视频| 亚洲性夜色夜夜综合| 91在线观看av| 久久久精品欧美日韩精品| 亚洲无线观看免费| 欧美成人性av电影在线观看| 国产黄片美女视频| 在线观看午夜福利视频| 999久久久精品免费观看国产| 毛片女人毛片| 国产极品精品免费视频能看的| 禁无遮挡网站| 国产一区二区三区视频了| 亚洲五月婷婷丁香| 亚洲国产高清在线一区二区三| 午夜老司机福利剧场| av福利片在线观看| 中亚洲国语对白在线视频| 婷婷亚洲欧美| 亚洲欧美日韩卡通动漫| 亚洲最大成人中文| 观看美女的网站| 一本综合久久免费| 黄色视频,在线免费观看| 亚洲国产欧美网| 露出奶头的视频| 亚洲午夜理论影院| 乱人视频在线观看| 在线观看舔阴道视频| 18禁在线播放成人免费| 亚洲一区二区三区不卡视频| 免费看美女性在线毛片视频| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 免费人成在线观看视频色| 免费高清视频大片| 成人欧美大片| 久久精品综合一区二区三区| 丁香六月欧美| 国产亚洲欧美在线一区二区| 真人做人爱边吃奶动态| 成年女人永久免费观看视频| 一区二区三区国产精品乱码| 久久香蕉国产精品| 美女被艹到高潮喷水动态| 国产三级黄色录像| 日韩欧美在线乱码| 国产精品精品国产色婷婷| 欧美乱色亚洲激情| 91在线精品国自产拍蜜月 | 国产黄a三级三级三级人| 99精品在免费线老司机午夜| 免费高清视频大片| 夜夜躁狠狠躁天天躁| 亚洲国产欧洲综合997久久,| 怎么达到女性高潮| 亚洲熟妇熟女久久| 国产欧美日韩精品一区二区| 国产精品一区二区三区四区久久| 午夜精品一区二区三区免费看| 有码 亚洲区| 波多野结衣巨乳人妻| 最近视频中文字幕2019在线8| 国产伦在线观看视频一区| 日日干狠狠操夜夜爽| h日本视频在线播放| 日韩高清综合在线| 不卡一级毛片| 在线观看免费午夜福利视频| 97碰自拍视频| 国产亚洲欧美98| 一级a爱片免费观看的视频| 精品久久久久久成人av| 国产精品久久久久久精品电影| 国产在线精品亚洲第一网站| 哪里可以看免费的av片| 看黄色毛片网站| 日本免费a在线| 欧美精品啪啪一区二区三区| 男女视频在线观看网站免费| 亚洲人成网站在线播放欧美日韩| 免费观看的影片在线观看| 久久久久久人人人人人| 在线看三级毛片| www日本黄色视频网| 亚洲男人的天堂狠狠| 丰满的人妻完整版| 久久人妻av系列| 长腿黑丝高跟| 午夜视频国产福利| 国产一区二区三区在线臀色熟女| 又黄又粗又硬又大视频| 91在线精品国自产拍蜜月 | 叶爱在线成人免费视频播放| 欧美精品啪啪一区二区三区| 亚洲五月天丁香| 精品国产亚洲在线| 国产精品美女特级片免费视频播放器| 嫩草影视91久久| 欧美丝袜亚洲另类 | x7x7x7水蜜桃|