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

    摩托車車架模態(tài)特性的數(shù)值模擬與試驗

    2014-06-27 05:46:36鄒喜紅石曉輝昌詩力
    關(guān)鍵詞:模態(tài)有限元振動

    熊 鋒,鄒喜紅,石曉輝,昌詩力

    (重慶理工大學(xué)汽車零部件制造及檢測技術(shù)教育部重點實驗室,重慶 400054)

    摩托車車架模態(tài)特性的數(shù)值模擬與試驗

    熊 鋒,鄒喜紅,石曉輝,昌詩力

    (重慶理工大學(xué)汽車零部件制造及檢測技術(shù)教育部重點實驗室,重慶 400054)

    摩托車車架的模態(tài)特性是車架其他動力學(xué)特性的基礎(chǔ)。為滿足某新款摩托車車架設(shè)計開發(fā)的需要,采用數(shù)值模擬和試驗驗證相結(jié)合的方法研究了車架的模態(tài)特性。利用Hyperworks有限元分析軟件進行了車架模態(tài)特性的數(shù)值模擬,提取和分析了車架的前10階固有頻率和振型,通過LMS.TEST.LAB數(shù)據(jù)采集分析軟件對車架進行了試驗?zāi)B(tài)分析,并采用模態(tài)擬合法和模態(tài)置信法(MAC)進行了模態(tài)試驗結(jié)果驗證,為車架的減振降噪和結(jié)構(gòu)優(yōu)化設(shè)計提供了理論依據(jù)。

    摩托車車架;模態(tài)特性;數(shù)值模擬;模態(tài)試驗;減振降噪

    摩托車車架是整車最關(guān)鍵的承載部件。由于路面激勵和發(fā)動機激勵兩大激勵源的存在,摩托車車架長期工作在劇烈的振動環(huán)境中,車架本身的振動特性對整車的操縱穩(wěn)定性、乘員舒適性起著決定性作用[1-2]。

    針對某摩托車企業(yè)某款新車型車架,本文結(jié)合解析模態(tài)分析(FEM)和試驗?zāi)B(tài)分析(EMA)兩種方法[3]來研究關(guān)鍵承載部件摩托車車架的模態(tài)特性,為車架新產(chǎn)品的結(jié)構(gòu)設(shè)計和性能評估提供了一個強有力的工具,也為圍繞其展開的后續(xù)的車架靜力學(xué)和動力學(xué)分析提供了可靠的模型和數(shù)據(jù),并為車架結(jié)構(gòu)優(yōu)化設(shè)計和減振降噪,提高穩(wěn)定性和舒適性提供了一種行之有效的方法。

    1 摩托車車架模態(tài)特性的數(shù)值模擬

    1.1 車架有限元建模

    車架的有限元建模是將車架離散為多自由度系統(tǒng),因此車架有限元模型既要如實反映車架實際結(jié)構(gòu)的力學(xué)特性,又要進行適當(dāng)?shù)暮喕蕴岣哂嬎阈?,主要考慮以下方面:

    1)去除對車架結(jié)構(gòu)剛度和力學(xué)特性影響不大的支耳和線束安裝夾。除此之外應(yīng)盡量保留原有結(jié)構(gòu),使車架的簡化不至對車架的模態(tài)特性有較大的影響。

    2)車架的結(jié)構(gòu)主要是圓管及鋼板,包括轉(zhuǎn)向立管、主梁、左右后管等,其厚度分布在1~5 mm,與圓管及鋼板的長度或其他截面尺寸相比較小,故采取抽中面和劃分SHELL板殼單元方法來模擬實際結(jié)構(gòu)。

    3)車架各組成部件之間通過焊接聯(lián)系在一起,劃分共節(jié)點板殼單元來模擬焊接關(guān)系。

    選取單元尺寸3~4 mm來進行車架網(wǎng)格劃分,整個有限元模型單元總數(shù)為96 160,節(jié)點數(shù)為95 836,其中QUADS單元數(shù)為92 620,TRIAS單元數(shù)為3 540,三角形單元所占比例為3.68%,符合三角形單元比例要求。車架材料為優(yōu)質(zhì)碳素鋼,彈性模量為2.06e11 N/m2,泊松比為0.3,密度為7 800 kg/m3。圖 1為最終構(gòu)建的車架有限元模型。

    圖1 車架有限元模型

    1.2 車架自由模態(tài)的數(shù)值仿真

    摩托車車架是一個多自由度彈性結(jié)構(gòu),假設(shè)車架已按某種方法離散為一個n自由度系統(tǒng),其動力學(xué)方程如下[4-5]:

    彈性體的自由振動總可以分解為一系列簡諧振動的疊加,設(shè)解為

    其中,{φ}是自由響應(yīng)的幅值列陣。將上式代入式(2)得:

    當(dāng){φ}非零時,這是一個廣義特征值問題,ω為特征值,{φ}為特征矢量。該方程有非零解的充要條件是其系數(shù)矩陣行列式為零,即

    式(5)是關(guān)于ω2的n次代數(shù)方程。設(shè)無重根,解此方程得ω的n個互異正根ωi(i=1,2,...n),通常將其按升序排列,則有:

    式(6)中,ωi分別對應(yīng)車架的i階模態(tài)頻率。此時由于車架是小阻尼結(jié)構(gòu),近似為無阻尼振動系統(tǒng),故模態(tài)頻率即固有頻率。將ωi(i=1,2,…,n)代入式(4)得到關(guān)于 {φi}中元素的具有n-1個獨立方程的代數(shù)方程組,共解得n個線性無關(guān)非零矢量{φi}的比例解,經(jīng)過歸一化即為車架主振型,即固有振型。將n個特征矢量 {φi}組成一個n階方陣,即為車架模態(tài)矩陣,見式(7)。

    用Hyperworks自帶的RADIOSS求解器求取方程式(4)時,有跟蹤法(tracking method)、變換法(transformation method)以及蘭索士(Lanczos method)法等多種方法[6]。考慮到蘭索士法是兼顧跟蹤法和變換法的矢量正交化方法,本文采取Lanczos方法進行車架自由模態(tài)的數(shù)值模擬。對于中小排量摩托車,感興趣的頻率范圍為0~300 Hz。限于篇幅,僅提取此頻率范圍內(nèi)車架的前5階固有頻率和前4階固有振型,并對各階次振型進行特征描述。表1為車架前5階固有頻率和振型描述,圖2為振型云圖。

    表1 車架固有頻率和振型

    圖2 車架1~4階固有振型云圖

    2 摩托車車架試驗?zāi)B(tài)分析

    2.1 測試原理

    試驗?zāi)B(tài)分析是基于系統(tǒng)響應(yīng)和激振力的動

    式(8)中:Xe(ω)為各測點響應(yīng)信號的傅里葉變換;Ff(ω)為激勵點激勵信號的傅里葉變換,通過LMS分析軟件,即可識別車架的模態(tài)參數(shù)和獲取一組模態(tài)振型。

    2.2 測試系統(tǒng)

    測試系統(tǒng)包括ICP型加速度傳感器、力錘、ICP型力傳感器、LMS.SCADAS.Ⅲ數(shù)采前端以及LMS.TEST.LAB模態(tài)分析軟件。鑒于LMS試驗設(shè)備采集通道的限制,對車架布置2個三向加速度傳感器和1個單向加速度傳感器。傳感器在車架位置布置上相互錯開,以便準(zhǔn)確獲取車架力錘激勵下不同位置不同方向的動態(tài)響應(yīng)。另外,為方便試驗結(jié)果與有限元結(jié)果相互驗證,試驗邊界條件應(yīng)與有限元模型自由邊界條件相一致,用彈性較好的橡皮繩將車架懸掛起來,使其處于接近自由支撐狀態(tài)。由于橡皮繩具有較低的支撐剛度和阻尼,車架的剛體共振頻率雖從0 Hz偏移到稍高一些的頻率上,但仍遠(yuǎn)低于車架1階彈性共振頻率,因此車架的彈性模態(tài)不會受到很大影響。另外,懸掛點應(yīng)選擇處于車架縱向?qū)ΨQ平面上并避免選擇柔性端,并使其處于或接近于盡可能多的模態(tài)節(jié)點上,這樣有利于準(zhǔn)確獲取車架的模態(tài)參數(shù)。車架模態(tài)試驗系統(tǒng)實物如圖3所示。態(tài)測試,即通過振動測試,經(jīng)信號處理和參數(shù)識別確定系統(tǒng)的模態(tài)參數(shù),建立以模態(tài)參數(shù)表示的運動方程。從模態(tài)分析理論可知,這些參數(shù)可以通過傳遞函數(shù)或頻響函數(shù)曲線進行分析求得。在時域內(nèi),系統(tǒng)的極點、固有頻率和阻尼比由系統(tǒng)的特征方程決定,即由傳遞函數(shù)的分母決定;在頻域內(nèi),根據(jù)輸入-輸出信號在頻域中的關(guān)系,可以采用模態(tài)參數(shù)識別法求得動態(tài)參數(shù)。

    由于車架是中小型結(jié)構(gòu),故采取單輸入多輸出方式(SIMO)對車架進行試驗?zāi)B(tài)分析。具體試驗時采取固定傳感器,移動力錘的方式,即對車架一次只激勵一個點的一個方向,而在其他任何坐標(biāo)上均沒有激勵作用。由頻響函數(shù)物理意義可知,通過激勵各測點f、測量一點e的響應(yīng),可得頻響函數(shù)矩陣中的一行元素[7-8]

    圖3 車架模態(tài)試驗系統(tǒng)

    2.3 測點布置和試驗參數(shù)設(shè)置

    單點激勵時,選擇測點應(yīng)以能有效激發(fā)起各階模態(tài)為原則,因此應(yīng)避免將測點選在結(jié)構(gòu)模態(tài)的節(jié)點(節(jié)線)處,也應(yīng)適當(dāng)遠(yuǎn)離結(jié)構(gòu)低階模態(tài)的反節(jié)點(節(jié)線)處。測點宜選擇車架自由端附近,且不宜選擇對稱面上。另外,車架結(jié)構(gòu)各部分剛度變化較大,測點宜選擇剛度較大的部位。除此之外,可以選取幾個激振點進行激勵試驗,測量若干個頻響函數(shù),觀察若哪個激振點激勵得到的頻響函數(shù)不丟失重要模態(tài),則此點為理想測點。除在所感興趣的區(qū)域進行多布測點之外,使測點在車架試件上呈某種程度的均勻分布,可減少漏掉模態(tài)的機會,也能夠得到較好的振型動畫,方便解析模態(tài)和試驗?zāi)B(tài)結(jié)果對比。參考解析模態(tài)分析結(jié)果,一共布置42個測點,如圖4所示。

    圖4 車架測點布置

    實際模態(tài)試驗時,由于車架空間結(jié)構(gòu)復(fù)雜多變,傳感器粘貼試件表面時,測量方向與部分模型節(jié)點坐標(biāo)系存在角度偏差,因此有必要調(diào)整部分節(jié)點歐拉角,使傳感器測量方向與模型節(jié)點坐標(biāo)系一致,如圖5所示。這樣能真實地反映車架各測點的振動情況,便于獲取準(zhǔn)確真實的模態(tài)參數(shù)和振型。

    圖5 車架歐拉角調(diào)整

    選取采用頻率為1 024 Hz,采樣時間為1 s,為消除誤差因素,每組測點采取3次取平均值。為抑制信號中的噪聲干擾,對力錘信號和響應(yīng)信號進行加指數(shù)窗處理。試驗時,應(yīng)保證力錘的輸入信號是單一時域脈沖信號,并且考慮低頻段(0~300 Hz)脈沖輸入信號和測點響應(yīng)信號之間的相干函數(shù)頻率幅值,保證其盡量大于0.8。由于車架是小阻尼試驗系統(tǒng),結(jié)合時域法和頻域法各自的優(yōu)缺點,采用時域法中的多參考點最小二乘復(fù)指數(shù)法(LSCE)來識別極點(頻率和阻尼),再利用最小二乘頻域法(LSFD)估計振型[9-10]。

    2.4 模態(tài)驗證

    模態(tài)模型的參數(shù)包括固有頻率、阻尼比、模態(tài)振型和模態(tài)參與因子,分別采用模態(tài)擬合法和模態(tài)置信判據(jù)(MAC)進行模態(tài)參數(shù)識別驗證。

    模態(tài)擬合法即對測量的FRF和擬合的FRF評估其相關(guān)性和誤差量。相關(guān)性是對擬合和測量的FRF值復(fù)共軛乘積的歸一化,見式(9)。

    其中:Si為擬合的FRF在譜線i處的復(fù)值;Mi為測量的FRF在譜線i處的復(fù)值;*表示共軛。

    對各個FRF,列出相關(guān)性低于某規(guī)定百分比的FRF和誤差高于規(guī)定百分比的FRF,對于評價擬合的FRF的質(zhì)量能提供有效的信息。限于篇幅,本文只列舉測點4(+Z)激勵,測點8 3個方向振動響應(yīng)的模態(tài)擬合相關(guān)性和誤差數(shù)據(jù),如表2??梢钥闯?,測點3個方向振動的頻響函數(shù),擬合值和實際測量值相關(guān)性高,誤差量較小,估計模態(tài)疊加得到的頻響函數(shù)可以很好地擬合出實測頻響函數(shù)。圖6為僅列舉的測點4(+Z)和測點8(-X)之間實測頻響函數(shù)和擬合的頻響函數(shù)曲線對比圖。其中實心線為實際測量FRF,雙點劃線為擬合FRF。從圖6可以看出,兩曲線在0~300 Hz頻帶范圍內(nèi)基本吻合,幅值和相位基本一致,只是在某些局部頻段出現(xiàn)少量偏差,但在合理的誤差范圍內(nèi)。另外,在每個極點固有頻率處兩曲線均吻合較好,很好地驗證了模態(tài)參數(shù)識別的準(zhǔn)確性和有效性。

    表2 模態(tài)擬合相關(guān)性和誤差

    圖6 頻響函數(shù)擬合圖

    響應(yīng)的模態(tài)置信判據(jù)(MAC)可以表示模態(tài)的可信程度,其公式為

    其中:{Rj},{Rl}分別表示第j階和第l階計算振型在n個測點所對應(yīng)的自由度值。如果復(fù)向量Rj與Rl之間線性相關(guān),則MAC值接近于1;相反,如果二者線性獨立,則MAC值將接近于零。MAC的概念可施加于兩個任意的復(fù)向量,即用于比較兩個有任意標(biāo)尺的模態(tài)振型向量,相似的模態(tài)振型具有高值的MAC。對于兩相對應(yīng)的模態(tài)而言,MAC值接近1,不同的模態(tài)MAC值接近于0。限于篇幅,僅提取車架前5階模態(tài)進行MAC置信驗證,結(jié)果見表3。

    表3 模態(tài)置信準(zhǔn)則(MAC)

    由于車架具有較好的線性特性,不同的模態(tài)應(yīng)該是線性無關(guān)的,因此車架的各階模態(tài)之間應(yīng)滿足:用百分?jǐn)?shù)表示時,同一模態(tài)的MAC值應(yīng)為100,不同模態(tài)向量的MAC值一般比較小。表3中,前5階模態(tài)的MAC值很好地滿足了這一特點,不過2階和3階、3階和4階之間MAC值稍大,這是由于該3階固有頻率值比較接近造成的。以上兩種驗證方法很好地驗證了模態(tài)試驗估計得到的模態(tài)參數(shù)是準(zhǔn)確有效的。

    3 數(shù)值模擬與模態(tài)試驗結(jié)果對比

    車架有限元數(shù)值模擬是實模態(tài),不考慮阻尼作用,車架試驗?zāi)B(tài)結(jié)果包含了車架阻尼作用,然而車架是小阻尼結(jié)構(gòu),因此車架的實模態(tài)和復(fù)模態(tài)頻率值應(yīng)比較接近,表4中前5階頻率結(jié)果對比很好地驗證了這一點。仔細(xì)觀察對比結(jié)果發(fā)現(xiàn),兩種方法在固有頻率值上吻合較好,頻率差最大為10.9 Hz,低于8%,在合理誤差范圍。振型方面,前10階振型均吻合較好,車架振動變形趨勢和幅度較為一致,數(shù)值模擬的結(jié)果和模態(tài)試驗的結(jié)果吻合較好,二者得到了很好的相互驗證,共同呈現(xiàn)了車架的振動特性。

    表4 數(shù)值模擬與模態(tài)試驗結(jié)果對比

    圖7僅列舉了前2階振型,用于數(shù)值模擬仿真結(jié)果和試驗?zāi)B(tài)分析結(jié)果對比。圖中虛線為試驗?zāi)B(tài)變形前車架,實線為試驗?zāi)B(tài)變形后車架,網(wǎng)格為解析模態(tài)變形前車架,云圖為解析模態(tài)變形后車架。

    圖7 數(shù)值模擬與模態(tài)試驗振型對比

    4 結(jié)論

    1)數(shù)值模擬結(jié)果和試驗結(jié)果吻合較好,誤差值較小,相互得到了很好的驗證。所建立的有限元模型正確有效,可以以此模型為基礎(chǔ)進行車架后續(xù)的靜、動態(tài)特性分析。

    2)從車架固有頻率上看,車架各階固有頻率均遠(yuǎn)離路面載荷激勵頻率(0~50 Hz),可以有效地規(guī)避共振失效,也驗證了車架結(jié)構(gòu)總體設(shè)計的合理性。

    3)從模態(tài)振型看,車架的振動形式以彎曲振動和扭轉(zhuǎn)振動為主,且車架尾部較車架前部振動幅度稍大,尤其是車架主、副駕駛員底部支管以及腳踏板支架處,結(jié)構(gòu)強度設(shè)計時應(yīng)予以重視。另外主駕駛員乘坐處支板局部振動幅度較大,可以通過改善焊接質(zhì)量、加厚支板或者添加加強筋的方式來改善局部的振動特性從而提升騎乘舒適性。

    [1] 石曉輝,王銳利,鄒喜紅,等.摩托車車架三通道道路模擬試驗裝置設(shè)計[J].機械設(shè)計與制造,2013(3): 194-196.

    [2] 王勇.摩托車車架強度有限元分析與驗證[J].機械設(shè)計與研究,2012,28(5):114-116.

    [3] Kumar V,Ranjan R,Sewak R,et al.Experimental Modal Analysis of Stator Overhangs of a Large Turbogenerator[J].Energy and Power Engineering,2011(3).

    [4] 毛睿甜,石曉輝.LX200型摩托車車架的模態(tài)特性分析[J].機械設(shè)計與制造,2007(7):40-42.

    [5] 柳小桐.摩托車車架模態(tài)特性分析[J].機械工程與自動化,2010(4):74-76.

    [6] 李偉,涂奎,李先文.摩托車車架模態(tài)分析與驗證[J].重慶理工大學(xué)學(xué)報:自然科學(xué)版,2011(3):11-14.

    [7] 曹樹謙,張文德,蕭龍翔.振動結(jié)構(gòu)模態(tài)分析[M].天津:天津大學(xué)出版社,2001.

    [8] 張力,林建龍,項輝宇.模態(tài)分析與試驗[M].北京:清華大學(xué)出版社,2011.

    [9] 于江,楊振冬,汪先國.摩托車車架試驗?zāi)B(tài)與解析模態(tài)分析[J].重慶工學(xué)院學(xué)報:自然科學(xué)版,2009(9): 25-32.

    [10]徐中明,張志飛,汪先國,等.摩托車車架模態(tài)特性分析[J].小型內(nèi)燃機與摩托車,2006(3):29-31.

    (責(zé)任編輯 何杰玲)

    Numerical Simulation and Test Research on Modal Characteristics of Motorcycle Frame

    XIONG Feng,ZOU Xi-h(huán)ong,SHI Xiao-h(huán)ui,CHANG Shi-li
    (Key Laboratory of Manufacture and Test Techniques for Automobile Parts,Ministry of Education,Chongqing University of Technology,Chongqing 400054,China)

    The modal characteristics of motorcycle frame are the basics of other dynamic characteristics.In order to meet the need of new frame design and development,method on the combination of numerical simulation and experimental verification was adopted to research on modal characteristics of the new type of motorcycle frame.Hyperworks finite element analysis software helped complete the numerical simulation of frame modal characteristics and the first ten order natural frequencies and vibration modes were extracted;at the same time experimental modal analysis of motorcycle frame was done through the aid of LMS.TEST.LAB data acquisition and analysis software and the modal test results were validated through modal synthesis and modal assurance method(MAC),thus providing theoretical basis for frame vibration noise reduction and structural optimization.

    motorcycle frame;modal characteristics;numerical simulation;experimental modal;vibration and noise reduction

    U467.3

    A

    1674-8425(2014)07-0015-06

    10.3969/j.issn.1674-8425(z).2014.07.004

    2014-04-18

    重慶市科技攻關(guān)計劃項目(cstc2011ggB0030,cstc2011ggB60010);重慶市教委資助項目(KJZH10211);重慶理工大學(xué)研究生創(chuàng)新基金資助項目(YCX2013201)

    熊鋒(1989—),男,碩士研究生,主要從事汽車測試技術(shù)研究。

    熊鋒,鄒喜紅,石曉輝,等.摩托車車架模態(tài)特性的數(shù)值模擬與試驗[J].重慶理工大學(xué)學(xué)報:自然科學(xué)版,2014(7):15-20.

    format:XIONG Feng,ZOU Xi-h(huán)ong,SHI Xiao-h(huán)ui,et al.Numerical Simulation and Test Research on Modal Characteristics of Motorcycle Frame[J].Journal of Chongqing University of Technology:Natural Science,2014 (7):15-20.

    猜你喜歡
    模態(tài)有限元振動
    振動的思考
    振動與頻率
    中立型Emden-Fowler微分方程的振動性
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    磨削淬硬殘余應(yīng)力的有限元分析
    UF6振動激發(fā)態(tài)分子的振動-振動馳豫
    計算物理(2014年2期)2014-03-11 17:01:44
    由單個模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
    計算物理(2014年2期)2014-03-11 17:01:39
    基于SolidWorks的吸嘴支撐臂有限元分析
    箱形孔軋制的有限元模擬
    上海金屬(2013年4期)2013-12-20 07:57:18
    国产免费视频播放在线视频 | 国产老妇伦熟女老妇高清| 三级男女做爰猛烈吃奶摸视频| 在线观看美女被高潮喷水网站| 亚洲成色77777| 国产乱人偷精品视频| 亚洲图色成人| 中文字幕制服av| 中文乱码字字幕精品一区二区三区 | 99热精品在线国产| 最近视频中文字幕2019在线8| 中文亚洲av片在线观看爽| 欧美xxxx性猛交bbbb| 成年女人看的毛片在线观看| 久久精品国产鲁丝片午夜精品| 国产精华一区二区三区| 成年女人永久免费观看视频| 成人三级黄色视频| 成人亚洲欧美一区二区av| 成人综合一区亚洲| 男人狂女人下面高潮的视频| 丝袜美腿在线中文| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品久久久久久久久久久久久| 亚洲欧洲日产国产| 男人舔女人下体高潮全视频| av免费观看日本| 日韩欧美精品v在线| 欧美日本亚洲视频在线播放| eeuss影院久久| 男人舔奶头视频| 特大巨黑吊av在线直播| 免费大片18禁| 欧美一区二区精品小视频在线| 欧美精品一区二区大全| 成人高潮视频无遮挡免费网站| 久99久视频精品免费| 最近手机中文字幕大全| 男插女下体视频免费在线播放| 好男人视频免费观看在线| 精品久久久久久久末码| 最近手机中文字幕大全| 乱系列少妇在线播放| 亚洲成人久久爱视频| 精品人妻熟女av久视频| 桃色一区二区三区在线观看| 国产精品一二三区在线看| 特大巨黑吊av在线直播| 日本免费一区二区三区高清不卡| 亚洲av二区三区四区| 欧美97在线视频| 亚洲精品乱久久久久久| 天堂√8在线中文| 日本免费在线观看一区| 亚洲天堂国产精品一区在线| 男女那种视频在线观看| 成人无遮挡网站| 26uuu在线亚洲综合色| av.在线天堂| 亚洲,欧美,日韩| 十八禁国产超污无遮挡网站| 男女那种视频在线观看| 欧美97在线视频| 国产成人午夜福利电影在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久国产网址| 色综合色国产| 欧美日韩综合久久久久久| 亚洲国产日韩欧美精品在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品永久免费网站| 好男人视频免费观看在线| 国产不卡一卡二| 真实男女啪啪啪动态图| 99热这里只有是精品在线观看| 久久午夜福利片| 国产精品爽爽va在线观看网站| 久久久成人免费电影| 亚洲乱码一区二区免费版| 国产午夜福利久久久久久| 18禁裸乳无遮挡免费网站照片| 亚洲欧洲日产国产| www.色视频.com| 国产av在哪里看| 插逼视频在线观看| 欧美97在线视频| 国产真实乱freesex| 日韩成人av中文字幕在线观看| 丝袜美腿在线中文| 久久久色成人| 非洲黑人性xxxx精品又粗又长| 非洲黑人性xxxx精品又粗又长| 超碰97精品在线观看| 69人妻影院| 亚洲欧洲日产国产| 国产av不卡久久| 日韩大片免费观看网站 | 99在线人妻在线中文字幕| 国内揄拍国产精品人妻在线| 亚洲av.av天堂| 国产一区二区亚洲精品在线观看| 一夜夜www| 一级毛片我不卡| 亚洲av免费高清在线观看| 久久精品国产亚洲网站| 精品99又大又爽又粗少妇毛片| 日本一本二区三区精品| 天堂av国产一区二区熟女人妻| 亚洲成av人片在线播放无| 欧美不卡视频在线免费观看| 欧美丝袜亚洲另类| 性插视频无遮挡在线免费观看| 深爱激情五月婷婷| 搡女人真爽免费视频火全软件| 国产乱人视频| 午夜福利网站1000一区二区三区| 人人妻人人看人人澡| 国产亚洲av片在线观看秒播厂 | 婷婷六月久久综合丁香| 波野结衣二区三区在线| 午夜福利高清视频| 国产极品天堂在线| 免费观看a级毛片全部| 国产私拍福利视频在线观看| 国产国拍精品亚洲av在线观看| 听说在线观看完整版免费高清| 91久久精品国产一区二区成人| 男女视频在线观看网站免费| 亚洲欧美清纯卡通| 久久99热这里只有精品18| 人体艺术视频欧美日本| 99热这里只有是精品50| 国产乱来视频区| 午夜激情欧美在线| 亚洲国产色片| 国产亚洲精品久久久com| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩卡通动漫| 一级二级三级毛片免费看| 国产一区二区亚洲精品在线观看| 男女边吃奶边做爰视频| 久久久久久久久大av| 久久久久久伊人网av| 中文欧美无线码| 精品国产三级普通话版| 黄片无遮挡物在线观看| 人人妻人人澡欧美一区二区| 久久精品国产亚洲av涩爱| 级片在线观看| 天堂影院成人在线观看| 99久久精品一区二区三区| 国语自产精品视频在线第100页| 男的添女的下面高潮视频| 26uuu在线亚洲综合色| 精品午夜福利在线看| 一个人免费在线观看电影| 毛片女人毛片| 国产高潮美女av| 日本-黄色视频高清免费观看| 国产精品1区2区在线观看.| 亚洲精品乱码久久久久久按摩| 黄片无遮挡物在线观看| 亚洲自拍偷在线| 久久99蜜桃精品久久| 免费播放大片免费观看视频在线观看 | 网址你懂的国产日韩在线| 日韩一本色道免费dvd| 国产一区二区在线观看日韩| 精品午夜福利在线看| av在线天堂中文字幕| 最后的刺客免费高清国语| 国产亚洲一区二区精品| 亚洲在线观看片| 秋霞在线观看毛片| 能在线免费观看的黄片| 黄片wwwwww| 国产高潮美女av| 亚洲欧美精品自产自拍| 日本黄大片高清| 大香蕉久久网| 欧美日韩在线观看h| 久久久久久九九精品二区国产| 一级黄色大片毛片| 99视频精品全部免费 在线| 美女脱内裤让男人舔精品视频| 日本五十路高清| 国产高清视频在线观看网站| 两个人视频免费观看高清| 午夜爱爱视频在线播放| 又爽又黄a免费视频| 国产免费男女视频| 国产精品综合久久久久久久免费| 热99在线观看视频| 最近视频中文字幕2019在线8| 午夜老司机福利剧场| 一级毛片我不卡| 国产精品久久电影中文字幕| 久久久a久久爽久久v久久| 噜噜噜噜噜久久久久久91| 男女下面进入的视频免费午夜| 精品熟女少妇av免费看| 永久网站在线| 国产爱豆传媒在线观看| 91午夜精品亚洲一区二区三区| 日产精品乱码卡一卡2卡三| 国产 一区精品| 国产伦精品一区二区三区视频9| 国产成人免费观看mmmm| 黑人高潮一二区| 欧美性猛交黑人性爽| 免费观看人在逋| 18禁在线无遮挡免费观看视频| 亚洲va在线va天堂va国产| 久久久久网色| 久久精品久久久久久久性| 亚洲精品乱码久久久v下载方式| 久久精品熟女亚洲av麻豆精品 | 久久精品久久精品一区二区三区| 午夜爱爱视频在线播放| 日韩一区二区三区影片| 国产真实乱freesex| 欧美激情在线99| 乱人视频在线观看| 丰满少妇做爰视频| 亚洲精品乱码久久久v下载方式| 日日摸夜夜添夜夜添av毛片| 国产av码专区亚洲av| 少妇丰满av| 18禁在线播放成人免费| 麻豆乱淫一区二区| 看黄色毛片网站| 午夜福利成人在线免费观看| 国产精品蜜桃在线观看| 久久精品国产鲁丝片午夜精品| 午夜久久久久精精品| 久久精品国产亚洲av涩爱| 一区二区三区高清视频在线| av播播在线观看一区| 久久婷婷人人爽人人干人人爱| 日本av手机在线免费观看| 久久精品国产亚洲网站| 美女被艹到高潮喷水动态| 国产精品久久久久久久电影| 日韩精品有码人妻一区| 成人午夜高清在线视频| 亚洲欧美日韩高清专用| 久久久久久久国产电影| 内射极品少妇av片p| 国产三级中文精品| 亚洲欧洲日产国产| 欧美xxxx性猛交bbbb| 99热精品在线国产| 两个人的视频大全免费| 亚洲成色77777| 日本色播在线视频| 少妇丰满av| 国产亚洲午夜精品一区二区久久 | 久久精品国产99精品国产亚洲性色| 欧美成人精品欧美一级黄| 高清毛片免费看| 久久人妻av系列| 亚洲人成网站在线观看播放| 精品熟女少妇av免费看| 男人舔奶头视频| www.色视频.com| 精品久久久久久久末码| 午夜精品国产一区二区电影 | 汤姆久久久久久久影院中文字幕 | 久久精品人妻少妇| 欧美性猛交黑人性爽| 亚洲五月天丁香| 久久99热6这里只有精品| 高清在线视频一区二区三区 | 国语自产精品视频在线第100页| 国内精品美女久久久久久| 国产亚洲5aaaaa淫片| 成人鲁丝片一二三区免费| 欧美xxxx黑人xx丫x性爽| 在线免费观看不下载黄p国产| 欧美成人午夜免费资源| 赤兔流量卡办理| 国产老妇女一区| 午夜福利在线在线| 国产精品,欧美在线| 我要搜黄色片| 国产淫片久久久久久久久| 高清毛片免费看| 久久亚洲国产成人精品v| 精品熟女少妇av免费看| 久久精品国产自在天天线| 乱人视频在线观看| 欧美色视频一区免费| 亚洲天堂国产精品一区在线| 亚洲在久久综合| 男人舔奶头视频| 精品久久国产蜜桃| av卡一久久| 人妻少妇偷人精品九色| 亚洲精品乱码久久久v下载方式| 18+在线观看网站| 久久热精品热| 狂野欧美激情性xxxx在线观看| 日韩 亚洲 欧美在线| 欧美变态另类bdsm刘玥| 日产精品乱码卡一卡2卡三| 深爱激情五月婷婷| 久久草成人影院| 成年版毛片免费区| 久久这里只有精品中国| 好男人视频免费观看在线| 亚洲国产成人一精品久久久| 欧美97在线视频| 夫妻性生交免费视频一级片| 久久99热这里只有精品18| 亚洲综合精品二区| 欧美97在线视频| 欧美精品一区二区大全| 精品久久国产蜜桃| 色播亚洲综合网| 尤物成人国产欧美一区二区三区| 大香蕉久久网| 亚洲欧美成人综合另类久久久 | 可以在线观看毛片的网站| 汤姆久久久久久久影院中文字幕 | 亚洲精品,欧美精品| 中国国产av一级| 国产v大片淫在线免费观看| 国内精品一区二区在线观看| 在线观看一区二区三区| 免费搜索国产男女视频| 久久草成人影院| 欧美一级a爱片免费观看看| .国产精品久久| 午夜日本视频在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 69av精品久久久久久| 午夜精品在线福利| 99热精品在线国产| www.av在线官网国产| 国产伦精品一区二区三区视频9| 亚洲aⅴ乱码一区二区在线播放| 精品无人区乱码1区二区| 啦啦啦观看免费观看视频高清| 日韩av在线免费看完整版不卡| 最近视频中文字幕2019在线8| 免费看日本二区| 亚洲成色77777| 又粗又硬又长又爽又黄的视频| 免费看av在线观看网站| 岛国毛片在线播放| 精品久久久久久电影网 | 亚洲,欧美,日韩| 国产伦精品一区二区三区视频9| 日本av手机在线免费观看| 亚洲va在线va天堂va国产| 一级爰片在线观看| 一区二区三区高清视频在线| 亚洲av福利一区| 午夜爱爱视频在线播放| 国产又黄又爽又无遮挡在线| 欧美+日韩+精品| 亚洲精品亚洲一区二区| 国产av不卡久久| 最近手机中文字幕大全| 乱码一卡2卡4卡精品| 精品国产露脸久久av麻豆 | 日韩成人伦理影院| 2021天堂中文幕一二区在线观| 中国国产av一级| 干丝袜人妻中文字幕| 久久久久免费精品人妻一区二区| 亚洲精品,欧美精品| 狠狠狠狠99中文字幕| 国产高清三级在线| 久久热精品热| 欧美激情在线99| 又黄又爽又刺激的免费视频.| 国产乱人视频| 一区二区三区四区激情视频| 午夜久久久久精精品| 久久欧美精品欧美久久欧美| 一区二区三区四区激情视频| 国产成人精品一,二区| 亚洲第一区二区三区不卡| 日本wwww免费看| 国产淫片久久久久久久久| 久久久久久久久久成人| 久久亚洲精品不卡| 国产人妻一区二区三区在| 搞女人的毛片| 97超视频在线观看视频| 午夜老司机福利剧场| videos熟女内射| 午夜福利高清视频| 日韩av在线大香蕉| 免费观看精品视频网站| 日本与韩国留学比较| 男插女下体视频免费在线播放| 国产成人精品一,二区| 麻豆久久精品国产亚洲av| 深爱激情五月婷婷| 亚洲精品456在线播放app| 白带黄色成豆腐渣| 久久人人爽人人爽人人片va| 超碰97精品在线观看| 久久精品91蜜桃| 国产黄a三级三级三级人| 亚洲欧美精品综合久久99| 日韩制服骚丝袜av| 青春草亚洲视频在线观看| 久久国内精品自在自线图片| 国产一级毛片在线| 久久精品国产鲁丝片午夜精品| 国产欧美日韩精品一区二区| 成人高潮视频无遮挡免费网站| 亚洲中文字幕一区二区三区有码在线看| 看非洲黑人一级黄片| 久久久精品欧美日韩精品| 日本一二三区视频观看| a级毛色黄片| 一级爰片在线观看| 免费av不卡在线播放| 91av网一区二区| 国产精品爽爽va在线观看网站| 一本久久精品| 国产在视频线在精品| 欧美激情国产日韩精品一区| 天天躁夜夜躁狠狠久久av| 日韩三级伦理在线观看| 午夜福利在线观看吧| 久热久热在线精品观看| 久久久久免费精品人妻一区二区| 欧美又色又爽又黄视频| 亚洲国产精品成人综合色| 日韩一区二区视频免费看| 边亲边吃奶的免费视频| 精品午夜福利在线看| 国产精品久久久久久久久免| 不卡视频在线观看欧美| 亚洲欧美日韩无卡精品| 91在线精品国自产拍蜜月| 极品教师在线视频| 黄色配什么色好看| 免费播放大片免费观看视频在线观看 | 午夜日本视频在线| kizo精华| 久久久午夜欧美精品| 欧美成人a在线观看| 久久久久久久久中文| 亚洲国产精品专区欧美| 男人和女人高潮做爰伦理| 伦理电影大哥的女人| 免费看美女性在线毛片视频| 亚洲天堂国产精品一区在线| 久久久午夜欧美精品| 久久综合国产亚洲精品| 亚洲久久久久久中文字幕| 麻豆久久精品国产亚洲av| 色综合站精品国产| 嘟嘟电影网在线观看| 精品国产一区二区三区久久久樱花 | 久久精品久久久久久久性| 亚洲中文字幕日韩| 搡老妇女老女人老熟妇| 亚洲图色成人| 一区二区三区高清视频在线| 成人漫画全彩无遮挡| 成人亚洲精品av一区二区| 永久网站在线| 在线观看美女被高潮喷水网站| 天堂影院成人在线观看| 亚洲精品乱码久久久久久按摩| 91久久精品国产一区二区成人| 99久久精品一区二区三区| eeuss影院久久| 国产精品久久久久久精品电影小说 | 校园人妻丝袜中文字幕| 国产久久久一区二区三区| 成年女人看的毛片在线观看| 18禁在线播放成人免费| 又爽又黄a免费视频| 久久久久久九九精品二区国产| 午夜福利在线观看吧| 免费看光身美女| 亚洲av电影在线观看一区二区三区 | 中国美白少妇内射xxxbb| 亚洲精品乱久久久久久| 午夜精品国产一区二区电影 | 欧美日韩精品成人综合77777| 春色校园在线视频观看| 欧美+日韩+精品| 99在线人妻在线中文字幕| 亚洲不卡免费看| 亚洲真实伦在线观看| 午夜免费激情av| 伊人久久精品亚洲午夜| 在线观看一区二区三区| 久久精品国产亚洲网站| 免费观看在线日韩| av在线播放精品| 久久久a久久爽久久v久久| 国产精品一区二区性色av| 日本猛色少妇xxxxx猛交久久| 美女内射精品一级片tv| 91精品国产九色| 免费av不卡在线播放| 亚洲av福利一区| 国产成人aa在线观看| 亚洲性久久影院| 国内揄拍国产精品人妻在线| 18禁裸乳无遮挡免费网站照片| 欧美97在线视频| 2021天堂中文幕一二区在线观| 一夜夜www| 国产69精品久久久久777片| 麻豆一二三区av精品| 一区二区三区高清视频在线| 国产成人免费观看mmmm| 国产精品熟女久久久久浪| 国产午夜精品一二区理论片| 亚洲美女视频黄频| 女人被狂操c到高潮| 欧美三级亚洲精品| 我要搜黄色片| 波多野结衣巨乳人妻| 草草在线视频免费看| 国产精品一二三区在线看| 国产午夜福利久久久久久| 搡老妇女老女人老熟妇| 中文资源天堂在线| 爱豆传媒免费全集在线观看| 久久久久久久国产电影| 欧美又色又爽又黄视频| 国产人妻一区二区三区在| 一区二区三区高清视频在线| 午夜福利视频1000在线观看| 三级男女做爰猛烈吃奶摸视频| 免费观看性生交大片5| 国产午夜精品论理片| 狂野欧美激情性xxxx在线观看| 国产欧美日韩精品一区二区| 亚洲av免费在线观看| 精品人妻一区二区三区麻豆| 亚洲美女视频黄频| 午夜精品国产一区二区电影 | 国国产精品蜜臀av免费| 可以在线观看毛片的网站| 高清av免费在线| 男女那种视频在线观看| 国产精品一及| 精品99又大又爽又粗少妇毛片| 久久99蜜桃精品久久| 久久精品国产亚洲网站| 在线观看av片永久免费下载| 2021少妇久久久久久久久久久| 黄色日韩在线| 欧美zozozo另类| 国产探花极品一区二区| av在线天堂中文字幕| 久久久欧美国产精品| 一边亲一边摸免费视频| 国产精品麻豆人妻色哟哟久久 | 国产一区二区三区av在线| 久久久久性生活片| 亚洲精品日韩在线中文字幕| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| 精品一区二区三区人妻视频| 国产精品久久视频播放| 亚洲最大成人中文| 狠狠狠狠99中文字幕| 国产毛片a区久久久久| 一二三四中文在线观看免费高清| 伊人久久精品亚洲午夜| 成人午夜高清在线视频| 国产精品.久久久| 一区二区三区免费毛片| 小说图片视频综合网站| 国产在视频线精品| 天堂中文最新版在线下载 | 又黄又爽又刺激的免费视频.| 高清日韩中文字幕在线| 国产高清三级在线| 高清av免费在线| 国产在线男女| 嫩草影院精品99| 伊人久久精品亚洲午夜| 日韩av在线大香蕉| 人人妻人人澡人人爽人人夜夜 | 亚洲精品成人久久久久久| 特大巨黑吊av在线直播| 热99在线观看视频| 亚洲精品一区蜜桃| 色综合亚洲欧美另类图片| 中文字幕人妻熟人妻熟丝袜美| av又黄又爽大尺度在线免费看 | 免费观看在线日韩| 日韩高清综合在线| 欧美区成人在线视频| 99久久精品国产国产毛片| 国产高潮美女av| 建设人人有责人人尽责人人享有的 | 午夜激情欧美在线| 亚洲激情五月婷婷啪啪| 狂野欧美白嫩少妇大欣赏| 亚洲精品日韩在线中文字幕| 午夜激情福利司机影院| 国产免费一级a男人的天堂| 成人av在线播放网站| av在线观看视频网站免费| 久久精品久久久久久久性| 免费观看a级毛片全部| 老司机福利观看| 亚洲欧洲国产日韩|