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

    基于流固耦合的磁流變減振器特性分析

    2020-05-15 10:32:28侯鎖軍趙向陽
    河南工學(xué)院學(xué)報 2020年1期
    關(guān)鍵詞:減振器氮?dú)?/a>活塞

    侯鎖軍,趙向陽

    基于流固耦合的磁流變減振器特性分析

    侯鎖軍,趙向陽

    (河南工學(xué)院 車輛與交通工程學(xué)院,河南 新鄉(xiāng) 453003)

    針對磁流變減振器示功特性出現(xiàn)畸形的問題,首先采用Hypermesh軟件建立磁流變減振器的液體模型和固體模型,將模型導(dǎo)入ADIAN液固耦合軟件并建立流固耦合模型,并通過減振器臺架示功特性試驗(yàn)驗(yàn)證了流固耦合模型的正確性;采用該流固耦合模型分析了不同速度、不同磁場及不同氮?dú)鈮毫ο碌膬?nèi)部速度分布和壓力分布,通過對不同氮?dú)馐覊毫Φ臏p振器內(nèi)部壓力場仿真,合理選擇了磁流變減振器內(nèi)的氮?dú)鈮毫?,避免了減振器示功特性出現(xiàn)畸形。這一磁流變減振器特性分析方法為避免磁流變減振器示功特性出現(xiàn)畸形提供了理論依據(jù)。

    流固耦合;磁流變減振器;臺架試驗(yàn);內(nèi)特性

    磁流變減振器的內(nèi)特性是指其內(nèi)部壓力及流速特性等,外特性指減振器的示功特性以及速度特性。磁流變減振器的性能設(shè)計需要工程師對減振器的結(jié)構(gòu)進(jìn)行反復(fù)調(diào)整才能獲得合理的減振器外特性曲線,因?yàn)橥馓匦曰兊臏p振器必然會影響其與懸架匹配的工作質(zhì)量,因此合理選擇減振器結(jié)構(gòu)參數(shù)來避免減振器外特性畸變是減振器設(shè)計的重要內(nèi)容[1]。有學(xué)者對普通減振器外特性的各種畸變及原因做了分析[2],這對磁流變減振器的設(shè)計有一定的參考意義。

    磁流變減振器是一種新型的阻尼可調(diào)減振器,利用磁流變液的粘度可控特性,能實(shí)現(xiàn)減振器阻尼力連續(xù)可調(diào),以滿足汽車行駛時不同工況下的減振要求[3-6]。磁流變減振器相對于普通減振器,其內(nèi)外腔灌裝磁流變液體粘度大,根據(jù)臺架試驗(yàn)可知,如果其氮?dú)馐覊毫x擇不合理,就會出現(xiàn)減振器外特性畸形現(xiàn)象,影響其與懸架的匹配質(zhì)量。

    本文擬通過流固耦合有限元理論建立磁流變減振器的流固耦合模型,研究磁流變減振器的內(nèi)部流場及壓力分布,為合理選擇減振器的氮?dú)鈮毫?,充分調(diào)試磁流變減振器的特性曲線提供理論參考。

    1 流固耦合有限元理論

    自動動態(tài)增量非線性分析(Automatic Dynamic Incremental Nonlinear Analysis,ADINA)軟件[7]是美國ADINA R&D公司開發(fā)的有限元通用分析仿真軟件,目前在各個工業(yè)領(lǐng)域應(yīng)用非常廣泛。ADINA軟件可解決多種物理場耦合的問題,如流固耦合、熱力耦合、流固熱耦合、聲流體和結(jié)構(gòu)耦合等。ADINA提供的gap邊界條件、參數(shù)化動網(wǎng)格、自適應(yīng)動網(wǎng)格、網(wǎng)格重畫分、強(qiáng)耦合等技術(shù)為減振器的流固耦合仿真提供了可靠的技術(shù)保障。同時,ADINA還具有穩(wěn)定的非線性問題求解專利技術(shù)。

    減振器的工作頻率一般低于50Hz,其內(nèi)部油液與結(jié)構(gòu)之間的耦合屬于低頻動態(tài)邊界條件耦合問題,流固耦合邊界上的基本運(yùn)動學(xué)和動力學(xué)條件是:

    上式對應(yīng)的條件為非滑移邊界,對應(yīng)滑移邊界,上式變?yōu)椋?/p>

    液固耦合計算過程中,液體作用于耦合邊界的力為:

    一般來說,在流固耦合分析中,為了得到更好的網(wǎng)格質(zhì)量和模擬精度,常常使用動網(wǎng)格技術(shù)(ALE),對于不可壓縮流體,其N-S方程為:

    上式為其連續(xù)性方程,其動量方程為

    如圖1所示,在流固耦合有限元分析中,液體單元和固體單元在流固耦合邊界相互作用,該作用導(dǎo)致流固耦合面是可變的,因此網(wǎng)格必須具有變形能力,否則網(wǎng)格易出現(xiàn)扭曲現(xiàn)象而導(dǎo)致計算終止,通??梢岳脛泳W(wǎng)格技術(shù)解決這一問題。

    圖1 流固耦合模型簡圖

    在利用有限元軟件處理流固耦合問題時,分別構(gòu)建液體和固體模型,其耦合關(guān)系通過定義耦合邊界條件實(shí)現(xiàn),在軟件中分別定義固體模型和液體模型不同的材料、邊界條件、流固耦合面、載荷及時間步長等。任何可能與流體發(fā)生相互作用的邊界都需要定義為FSI邊界,有些邊界在計算開始時可能沒有與液體接觸,但在一段時間以后會和液體接觸,這些在流固耦合計算過程中所有可能與液體接觸的面都要設(shè)置為流固耦合面。

    2 磁流變減振器流固耦合模型建立

    2.1 磁流變減振器結(jié)構(gòu)模型

    在建立流固耦合有限元模型之前,首先建立磁流變減振器的三維CATIA模型,如圖2所示。該模型在原有CAD模型的基礎(chǔ)上做了一些簡化,將減振器外筒和上下端蓋合成一體,去掉了螺栓孔、密封圈以及緊定螺釘?shù)取?/p>

    圖2 磁流變減振器三維模型

    2.2 磁流變減振器流固耦合模型建立

    磁流變減振器流固耦合仿真過程中殼體的變形對減振器的運(yùn)動特性可以忽略不計,只考慮內(nèi)筒活塞和氮?dú)馐一钊倪\(yùn)動特性對減振器特性的影響。磁流變減振器氮?dú)馐覂?nèi)氣體對氮?dú)馐一钊膲毫烧J(rèn)為是線性的,在0~10MPa壓力范圍和0~200℃溫度范圍內(nèi),氮?dú)饪烧J(rèn)為是理想氣體,氮?dú)庠诖帕髯儨p振器中的工作過程可認(rèn)為是絕熱過程,因此在流固耦合模型中,氮?dú)饪捎脧椈蓡卧妗8鶕?jù)公式(6)可計算出氮?dú)馐一钊幱谥虚g行程壓力為1MPa時的氮?dú)鈴椈商匦?,如圖3所示。氮?dú)馐业膲毫?i>P可表示為:

    式中,Pn表示氮?dú)馐一钊幱谥虚g行程時的壓力(Pa),Vn表示氮?dú)馐一钊幱谥虚g行程時的氮?dú)馐殷w積(m3),VN表示活塞處于任意位置時的氮?dú)馐殷w積(m3),xp表示氮?dú)馐一钊灰疲╩)。

    將磁流變減振器建立的三維CATIA模型以IGS的格式導(dǎo)入有限元前處理軟件Hypermesh中,在Hypermesh軟件中分別建立固體網(wǎng)格模型和液體網(wǎng)格模型。

    固體模型僅僅需要活塞桿、活塞及氮?dú)饣钊哪P?,如圖4所示;液體模型是基于磁流變減振器的三維CATIA固體模型而建立,在Hypermesh軟件中畫出固體模型所包圍的液體部分,如圖5所示,其中上通道液體模型和下通道液體模型為受磁場控制的部分。固體和液體模型共102 475個六面體網(wǎng)格,劃分網(wǎng)格后給每個零件組賦予材料屬性,然后存儲成nas格式導(dǎo)入流固耦合ADINA軟件中,在ADINA軟件中通過Node Set定義壁面、流固耦合面、約束和載荷[8]。

    圖4 固體網(wǎng)格模型

    圖5 液體網(wǎng)格模型

    在固體模型中,活塞桿及活塞材料都為鋼,密度為7.8×10-9t/mm3,泊松比為0.29,彈性模量為2.07×105N/mm2;定義載荷時間函數(shù),激勵幅值設(shè)定為25mm,時間函數(shù)按照激勵頻率設(shè)定,如圖6所示。約束活塞桿及氮?dú)饣钊倪\(yùn)動方向,此模型設(shè)定只能在方向運(yùn)動。設(shè)定活塞上端面、活塞下端面和氮?dú)馐一钊露嗣鏋轳詈厦?,分別設(shè)定為Interface1、Interface2和Interface3,這三個流固耦合面和液體模型中設(shè)定的耦合面要一一對應(yīng)。定義彈簧單元并設(shè)定彈簧剛度,如圖7所示。

    在液體模型中依次定義磁流變液材料參數(shù)、流固耦合面、壁面及位移約束面等,如圖8所示。磁流變液動力粘度為2.9×10-7MPa.s,密度為2.8×10-9t/mm3。在ADINA流固耦合計算過程中載荷在固體模型中設(shè)置,計算時間步長取決于液體模型中的時間步長,因此時間步長在液體模型中設(shè)置。

    圖6 活塞桿上端面位移載荷

    圖7 流固耦合面及氮?dú)鈴椈?/p>

    圖8 液體模型Wall和FSI邊界條件

    3 磁流變減振器流固耦合模型驗(yàn)證

    流固耦合模型的幾點(diǎn)假設(shè):從磁流變液的試驗(yàn)中發(fā)現(xiàn),溫度對磁流變液的粘度和剪切應(yīng)力的影響可以忽略不計,因此模型中不考慮溫度對減振器的影響,為了避免橡膠襯套的彈性影響示功曲線的形狀,在流固耦合模型中去掉了減振器兩端吊耳內(nèi)的橡膠襯套,工作間隙定為1.5mm,上下端阻尼通道單側(cè)長度均為30mm,激勵振幅選擇為25mm,激勵頻率選擇0.33Hz、0.83Hz、1.66Hz、2.48Hz和3.33Hz,對應(yīng)的減振器最大速度分別為0.052m/s、0.13m/s、0.26m/s、0.39m/s和0.52m/s。

    圖9所示為激勵頻率1.66Hz不同電流下的示功特性曲線,圖10所示為不同激勵頻率不同電流下的速度特性曲線,可以看出,理論值與試驗(yàn)值大小、趨勢一致。由于臺架自身存在不穩(wěn)定性、臺架傳感器測試精度存在誤差以及磁流變減振器制造和安裝工藝的誤差,使得數(shù)據(jù)存在一定波動,因此仿真與試驗(yàn)存在一定的誤差,不同工況下的仿真值和試驗(yàn)值阻尼力最大差值小于125N,最大誤差為4.5%。

    圖9 磁流變減振器示功特性仿真和試驗(yàn)結(jié)果

    圖10 磁流變減振器速度特性仿真和試驗(yàn)結(jié)果

    4 磁流變減振器內(nèi)特性分析

    4.1 不同速度下流場及壓力場分布

    對比磁流變減振器工作頻率分別為0.33Hz、1.66Hz和3.33Hz下減振器內(nèi)部液體流動速度分布及壓力分布,其他設(shè)定條件為氮?dú)鈮毫?MPa,阻尼通道上下端粘度設(shè)定為100Pa.s,相當(dāng)于1A的電流,其他部分粘度設(shè)置為0.29Pa.s,圖11~16所示為活塞處于拉伸行程中間位置時三種頻率下的速度分布圖。

    圖11 0.33Hz速度三維分布

    圖12 0.33Hz速度截面分布

    圖13 1.66Hz速度三維分布

    圖14 1.66Hz速度截面分布

    圖15 3.33Hz速度三維分布

    圖16 3.33Hz速度截面分布

    由圖11~16可以看出,活塞帶動其上方磁流變液流入頂端環(huán)形通道,在環(huán)形通道入口處,通道面積減小,形成擠流,流速增加,不同活塞速度下,頂端環(huán)形通道的流動速度不同,活塞運(yùn)動頻率越大,環(huán)形通道液體流動速度越大。仿真結(jié)果也符合這一規(guī)律。

    圖17為0.33Hz、1.66Hz和3.33Hz三種頻率下減振器不同部分的粘度值。圖18~20所示為0.33Hz、1.66Hz和3.33Hz三種頻率下活塞處于拉伸中間行程的壓力分布圖??梢钥闯龌钊隙说江h(huán)形通道之間的液體壓力均勻分布,大小變化不明顯;而在上下端環(huán)形阻尼通道處,由于流通面積減小,其液體流動速度增大,液體壓力迅速減小,并且在環(huán)形通道處液體壓力呈梯度減小。隨著激勵頻率的增加,液體流經(jīng)環(huán)形通道后壓力下降幅值逐漸增大,三種頻率下,上端通道壓力降分別為0.234MPa、0.400MPa和0.800MPa;下端通道壓力降分別為0.233MPa、0.400MPa和0.760MPa,頻率越大阻尼通道壓差越大,活塞桿拉伸的阻力就越大。

    圖17 粘度分布

    圖18 0.33Hz壓力分布

    圖19 1.66Hz壓力分布

    圖20 3.33Hz壓力分布

    4.2 不同磁場下流場及壓力場分布

    對比磁流變減振器上下端線圈施加磁場(上下端阻尼通道液體粘度設(shè)定為100Pa.s,相當(dāng)于1A電流,其余部分為0.29Pa.s)和不施加磁場下(液體粘度全部設(shè)定為0.29Pa.s)的內(nèi)部液體流場及壓力分布,其他設(shè)定條件為減振器激勵頻率為3.33Hz,氮?dú)鈮毫?MPa。有圖21和22所示,零磁場下和施加磁場下速度分布一致;由圖23~26可以看出,零磁場與施加磁場兩種情況下磁流變液流過上端阻尼通道后壓力降分別為0.24MPa和0.750MPa,零磁場與施加磁場兩種情況下磁流變液流過下端通道后壓力降分別為0.23MPa和0.748MPa。施加磁場后液體流經(jīng)環(huán)形通道后壓力下降值明顯增大,因此通過磁場控制磁流變液黏度可以實(shí)現(xiàn)對磁流變減振器阻尼力的控制。

    圖21 不加磁場流場分布

    圖22 施加磁場流場分布

    圖23 不加磁場上端通道壓力截面分布

    圖24 不加磁場下端通道壓力截面分布

    圖25 施加磁場上端通道壓力截面分布

    圖26 施加磁場下端通道壓力截面分布

    4.3 不同氮?dú)鈮毫ψ饔孟麓帕髯儨p振器壓力場及流速分布

    對比激振頻率為1.66Hz,上下端阻尼通道粘度為100Pa.s,相當(dāng)于磁場線圈施加1A電流,其他部分液體粘度為0.29Pa.s,氮?dú)鈮毫Ψ謩e為0、0.5和1MPa時,磁流變減振器的流速及壓力分布如圖27和28所示。從圖27速度分布圖看,活塞處于拉伸行程中間位置時,活塞上方液體在上端阻尼通道入口處形成擠流,由上端阻尼通道處流向內(nèi)外缸筒之間,然后液體由下端阻尼通道流入活塞下端空腔,并在下端阻尼通道出口處形成旋流。從速度分布圖可以看出,其速度分布合理,并且不同氮?dú)鈮毫ο滤俣确植夹螒B(tài)相似。

    圖27 速度分布

    由壓力分布圖28可以看出,在氮?dú)鈮毫?時,活塞下方的壓力為負(fù)值;氮?dú)鈮毫?.5MPa時,下方壓力為0.0667MPa;氮?dú)鈮毫υ黾拥?MPa時,下方壓力增加到大約0.6MPa。結(jié)果表明,沒有氮?dú)鈮毫虻獨(dú)鈮毫^低時減振器的示功圖都會出現(xiàn)畸形,該仿真結(jié)果為選擇合理的氮?dú)鈮毫μ峁┝死碚撘罁?jù)。

    圖28 壓力分布

    圖29為磁流變減振器在不同氮?dú)鈮毫ο碌脑囼?yàn)示功曲線,可以看出在氮?dú)鈮毫?時,示功曲線出現(xiàn)畸形,曲線左上端和右下端出現(xiàn)凹陷現(xiàn)象,其原因是在拉伸行程中磁流變液體不能迅速流回活塞下腔,在壓縮行程中磁流變液體不能迅速流回活塞上腔,導(dǎo)致在拉伸和壓縮行程初期,出現(xiàn)空行程狀態(tài),從而使示功曲線出現(xiàn)畸形現(xiàn)象;當(dāng)?shù)獨(dú)鈮毫?.5MPa時,凹陷減小;當(dāng)?shù)獨(dú)鈮毫μ岣叩?MPa時,凹陷現(xiàn)象逐漸消失。此試驗(yàn)結(jié)果也驗(yàn)證了流固耦合有限元仿真結(jié)果。根據(jù)不同工況下反復(fù)試驗(yàn)的結(jié)果,最終確定磁流變減振器的氮?dú)鈮毫?MPa。

    圖29 不同氮?dú)鈮毫Υ帕髯儨p振器示功特性曲線

    5 結(jié)論

    (1)本文通過流固耦合有限元分析軟件ADINA建立了磁流變減振器的流固耦合模型,仿真分析了不同速度、不同磁場、不同氮?dú)鈮毫ο碌臏p振器速度場和壓力場分布,仿真結(jié)果符合減振器實(shí)際內(nèi)部液體運(yùn)動規(guī)律。

    (2)通過對不同氮?dú)馐覊毫ψ饔孟碌臏p振器內(nèi)部壓力場仿真分析,合理選擇了磁流變減振器內(nèi)的氮?dú)鈮毫?,避免了減振器示功特性出現(xiàn)畸形現(xiàn)象,為磁流變減振器氮?dú)鈮毫Φ倪x擇提供了理論依據(jù)。

    [1] 孫勝利.位移相關(guān)減振器動力學(xué)建模及對車輛性能影響的研究[D].長春:吉林大學(xué),2008.

    [2] 俞德孚.車輛懸架減振器的理論和實(shí)踐[M].北京:兵器工業(yè)出版社,2003.

    [3] 劉韶慶. 磁流變可調(diào)阻尼減振器的特性研究[D]. 鎮(zhèn)江: 江蘇大學(xué), 2007.

    [4] WERELEY N M, CHO J U, CHOI Y T, et al. Magnetorheological dampers in shear mode [J]. Smart materials and structures, 2008, 17(1): 1-11.

    [5] HONG S R, CHOI S B, CHOI Y T, et al. Non-dimensional analysis and design of a magnetorheological damper[J]. Journal of sound and vibration, 2005, 288(4-5): 847- 863.

    [6] 陳杰平.基于磁流變減振器的汽車半主動懸架設(shè)計與控制研究[D].合肥:合肥工業(yè)大學(xué),2010.

    [7] 岳戈,梁宇白,陳晨,等. ADINA流體與流固耦合功能的高級應(yīng)用[M].北京:人民交通出版社, 2010.

    [8] 馬野, 袁志丹,曹金鳳. ADINA有限元經(jīng)典實(shí)例分析[M]. 北京: 機(jī)械工業(yè)出版社,2011.

    Characteristics Analysis of Magneto-rheological Damper Based on Fluid-solid Coupling

    HOU Suo-jun, ZHAO Xiang-yang

    (School of Vehicle and Traffic Engineering, Henan Institute of Technology, Xinxiang 453003, China)

    In order to solve the problem of abnormal force-displacement characteristics of magneto-rheological(MR) damper, firstly, the fluid model and solid model of MR damper were established by Hyper mesh software, the model was imported into ADINA fluid-solid coupling software to establish the fluid-solid coupling model, and the correctness of the fluid-solid coupling model was verified by the bench test of the force-displacement characteristics of the damper. The fluid-solid coupling model was used to analyze the internal velocity distribution and pressure distribution under different velocity, magnetic field and nitrogen pressure. Through the simulation of the internal pressure field of the MR damper with different nitrogen chamber pressure, the nitrogen pressure in the MR damper was reasonably selected to avoid the deformity of the force-displacement characteristics of the MR damper. The characteristic analysis method of MR damper studied in this paper provides a theoretical basis for avoiding the deformity of the force-displacement characteristic of MR damper.

    fluid-solid coupling; MR damper; bench test; internal characteristics

    U461.4

    A

    2096–7772(2020)01–0039–07

    2019-12-27

    河南省科技攻關(guān)項(xiàng)目(182102210034);河南省高等學(xué)校青年骨干教師培養(yǎng)計劃項(xiàng)目(2016GGJS-196,2018GCJS172)

    侯鎖軍(1979―),男,河南輝縣人,副教授,博士,主要從事汽車系統(tǒng)動力學(xué)及懸架控制研究。

    (責(zé)任編輯呂春紅)

    猜你喜歡
    減振器氮?dú)?/a>活塞
    一種活塞擠壓式室內(nèi)通風(fēng)設(shè)備的研發(fā)與應(yīng)用
    海上平臺氮?dú)庀到y(tǒng)流程工藝
    汽車減振器與磁流變材料分析與運(yùn)算
    氮?dú)馀菽瓑毫岩盒阅芗皯?yīng)用評價
    小天使·六年級語數(shù)英綜合(2014年10期)2014-10-20 07:03:37
    饋能式磁流變減振器自供電特性研究
    汽車科技(2014年6期)2014-03-11 17:45:36
    裝多級扭轉(zhuǎn)減振器的發(fā)動機(jī)曲軸系統(tǒng)扭振分析
    氮?dú)狻?xì)水霧滅火關(guān)鍵技術(shù)研究
    河南科技(2014年19期)2014-02-27 14:15:31
    KS Kolbenschmidt公司的新型鋼活塞
    低噪聲活塞——降低活塞銷的拍擊噪聲
    99久久综合精品五月天人人| 国产视频一区二区在线看| 国产精品影院久久| 长腿黑丝高跟| 麻豆国产av国片精品| 亚洲av不卡在线观看| 小蜜桃在线观看免费完整版高清| 麻豆久久精品国产亚洲av| 国内揄拍国产精品人妻在线| 夜夜爽天天搞| 日本在线视频免费播放| 成人国产综合亚洲| 国产不卡一卡二| 免费大片18禁| 18禁黄网站禁片午夜丰满| 亚洲欧美日韩高清在线视频| 欧美日韩黄片免| 香蕉av资源在线| 99久久无色码亚洲精品果冻| 久久伊人香网站| 色老头精品视频在线观看| 男女午夜视频在线观看| 国产伦人伦偷精品视频| 欧美丝袜亚洲另类 | 久久伊人香网站| 日韩精品中文字幕看吧| 午夜福利18| or卡值多少钱| 欧美丝袜亚洲另类 | 欧美乱码精品一区二区三区| 97碰自拍视频| 国产精品自产拍在线观看55亚洲| e午夜精品久久久久久久| 精品久久久久久久人妻蜜臀av| 我的老师免费观看完整版| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美一区二区三区黑人| 校园春色视频在线观看| av天堂中文字幕网| 国产私拍福利视频在线观看| 青草久久国产| 可以在线观看毛片的网站| 手机成人av网站| 在线观看av片永久免费下载| 大型黄色视频在线免费观看| 精品人妻偷拍中文字幕| 午夜免费成人在线视频| av天堂在线播放| 中文亚洲av片在线观看爽| 麻豆久久精品国产亚洲av| 国内少妇人妻偷人精品xxx网站| 我的老师免费观看完整版| 人妻夜夜爽99麻豆av| 免费在线观看影片大全网站| 999久久久精品免费观看国产| 日韩精品青青久久久久久| 日本 av在线| 色播亚洲综合网| 一卡2卡三卡四卡精品乱码亚洲| 99国产精品一区二区三区| 少妇人妻一区二区三区视频| 欧美黑人巨大hd| 国产成人啪精品午夜网站| 国产精品国产高清国产av| 亚洲精品成人久久久久久| 日本免费a在线| 欧美黄色淫秽网站| 老司机在亚洲福利影院| 黄片小视频在线播放| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 久久精品人妻少妇| 香蕉丝袜av| 欧美中文日本在线观看视频| 网址你懂的国产日韩在线| 男人舔女人下体高潮全视频| 男女视频在线观看网站免费| 亚洲av成人精品一区久久| 国产一级毛片七仙女欲春2| 亚洲av免费在线观看| 99riav亚洲国产免费| 天堂影院成人在线观看| 一进一出抽搐gif免费好疼| www.999成人在线观看| 99在线视频只有这里精品首页| 国产探花在线观看一区二区| 91字幕亚洲| 露出奶头的视频| 国产日本99.免费观看| 国产熟女xx| 亚洲精品在线观看二区| 欧美性猛交╳xxx乱大交人| 免费在线观看亚洲国产| 日本熟妇午夜| 免费电影在线观看免费观看| av女优亚洲男人天堂| 久久久久久久亚洲中文字幕 | eeuss影院久久| 日韩av在线大香蕉| 精品一区二区三区av网在线观看| 男人舔奶头视频| 中文字幕熟女人妻在线| 老司机午夜福利在线观看视频| 日韩精品中文字幕看吧| 天堂动漫精品| 最新在线观看一区二区三区| 免费看光身美女| 男女下面进入的视频免费午夜| 欧美成人a在线观看| 国产高清videossex| 亚洲国产高清在线一区二区三| 日本精品一区二区三区蜜桃| 亚洲人与动物交配视频| 激情在线观看视频在线高清| 亚洲精品色激情综合| 一级作爱视频免费观看| 欧美+日韩+精品| www日本黄色视频网| 国产精品久久久久久人妻精品电影| 欧美一级毛片孕妇| 国产老妇女一区| 亚洲av成人av| 老司机午夜十八禁免费视频| 少妇高潮的动态图| 岛国视频午夜一区免费看| 欧美成狂野欧美在线观看| 夜夜躁狠狠躁天天躁| 日韩免费av在线播放| 黄色视频,在线免费观看| 欧美成人免费av一区二区三区| 黄色日韩在线| 国产亚洲欧美在线一区二区| 色综合欧美亚洲国产小说| 三级男女做爰猛烈吃奶摸视频| 岛国视频午夜一区免费看| 女警被强在线播放| 亚洲av五月六月丁香网| 19禁男女啪啪无遮挡网站| 无人区码免费观看不卡| av天堂中文字幕网| 国产一区二区三区在线臀色熟女| 国产av一区在线观看免费| 悠悠久久av| 极品教师在线免费播放| 亚洲男人的天堂狠狠| 一区二区三区激情视频| av在线天堂中文字幕| 亚洲av电影不卡..在线观看| 一区二区三区高清视频在线| 久久久久久久精品吃奶| 怎么达到女性高潮| 美女大奶头视频| 国产精品久久久久久人妻精品电影| 久久国产精品人妻蜜桃| 国产午夜精品久久久久久一区二区三区 | www.色视频.com| 禁无遮挡网站| 国产探花极品一区二区| 亚洲精品美女久久久久99蜜臀| 中文字幕精品亚洲无线码一区| 亚洲va日本ⅴa欧美va伊人久久| 九九久久精品国产亚洲av麻豆| 亚洲专区国产一区二区| АⅤ资源中文在线天堂| 国产淫片久久久久久久久 | 小说图片视频综合网站| 久久香蕉国产精品| 欧美中文综合在线视频| 欧美一区二区精品小视频在线| 最近视频中文字幕2019在线8| 国产97色在线日韩免费| 色哟哟哟哟哟哟| 黄色视频,在线免费观看| 亚洲精品乱码久久久v下载方式 | 99热这里只有精品一区| 少妇高潮的动态图| av国产免费在线观看| 搡女人真爽免费视频火全软件 | 成年女人毛片免费观看观看9| 小说图片视频综合网站| 制服丝袜大香蕉在线| 欧美另类亚洲清纯唯美| 男人舔奶头视频| 日韩欧美免费精品| 淫妇啪啪啪对白视频| av天堂在线播放| 中亚洲国语对白在线视频| 国产激情欧美一区二区| 午夜a级毛片| 国产成人影院久久av| 真人做人爱边吃奶动态| 国产蜜桃级精品一区二区三区| 99精品久久久久人妻精品| 老司机深夜福利视频在线观看| 色吧在线观看| 久久草成人影院| 亚洲 欧美 日韩 在线 免费| 搡老妇女老女人老熟妇| 午夜视频国产福利| 国产91精品成人一区二区三区| 少妇人妻精品综合一区二区 | 色哟哟哟哟哟哟| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 亚洲激情在线av| 三级毛片av免费| 亚洲色图av天堂| 亚洲内射少妇av| 特大巨黑吊av在线直播| 嫩草影视91久久| 亚洲精华国产精华精| 精品久久久久久久人妻蜜臀av| 亚洲片人在线观看| av片东京热男人的天堂| 两个人的视频大全免费| 中国美女看黄片| 国产激情偷乱视频一区二区| 91麻豆精品激情在线观看国产| 非洲黑人性xxxx精品又粗又长| 一个人看的www免费观看视频| 国产三级黄色录像| 国产在视频线在精品| 婷婷精品国产亚洲av在线| 欧美xxxx黑人xx丫x性爽| 波多野结衣巨乳人妻| 国产真人三级小视频在线观看| 国产成人影院久久av| 一本综合久久免费| 琪琪午夜伦伦电影理论片6080| 亚洲不卡免费看| 国内少妇人妻偷人精品xxx网站| 成人午夜高清在线视频| 一区二区三区国产精品乱码| 色吧在线观看| 亚洲精品456在线播放app | 欧美精品啪啪一区二区三区| 亚洲无线在线观看| 欧洲精品卡2卡3卡4卡5卡区| 日本五十路高清| 俺也久久电影网| 18+在线观看网站| 国产91精品成人一区二区三区| 美女高潮的动态| 精品一区二区三区人妻视频| 中国美女看黄片| 亚洲无线观看免费| 女生性感内裤真人,穿戴方法视频| 日本三级黄在线观看| 97碰自拍视频| 一级作爱视频免费观看| 精品一区二区三区视频在线观看免费| 国模一区二区三区四区视频| 欧美精品啪啪一区二区三区| 欧美一区二区亚洲| 色综合婷婷激情| 黄片小视频在线播放| 亚洲欧美精品综合久久99| 国产精品综合久久久久久久免费| 黄片大片在线免费观看| 国产精品一区二区三区四区免费观看 | 亚洲久久久久久中文字幕| 好男人电影高清在线观看| 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 国产av不卡久久| 亚洲精品乱码久久久v下载方式 | 欧美一级a爱片免费观看看| 一级a爱片免费观看的视频| 一进一出抽搐动态| 一区二区三区高清视频在线| 欧美成狂野欧美在线观看| 又黄又爽又免费观看的视频| 偷拍熟女少妇极品色| 免费观看精品视频网站| 精品乱码久久久久久99久播| 久久伊人香网站| 无遮挡黄片免费观看| 麻豆成人午夜福利视频| 日本黄大片高清| 欧美黑人巨大hd| 国产精品亚洲美女久久久| 久久中文看片网| 色综合欧美亚洲国产小说| 久久欧美精品欧美久久欧美| 在线观看免费午夜福利视频| 国内精品美女久久久久久| 夜夜爽天天搞| 手机成人av网站| 18禁在线播放成人免费| 91在线观看av| 噜噜噜噜噜久久久久久91| 亚洲在线观看片| 九九久久精品国产亚洲av麻豆| 久久精品国产99精品国产亚洲性色| 在线观看av片永久免费下载| 国产淫片久久久久久久久 | 此物有八面人人有两片| 国产真实乱freesex| 午夜两性在线视频| a在线观看视频网站| 老熟妇乱子伦视频在线观看| 岛国在线免费视频观看| 亚洲欧美日韩东京热| 久久久精品欧美日韩精品| 国产色爽女视频免费观看| 最新中文字幕久久久久| 色av中文字幕| av国产免费在线观看| av片东京热男人的天堂| 2021天堂中文幕一二区在线观| 在线观看免费午夜福利视频| 成人av一区二区三区在线看| 久久久色成人| 亚洲一区二区三区色噜噜| 亚洲精品久久国产高清桃花| 国产精品99久久99久久久不卡| 国内揄拍国产精品人妻在线| 成人亚洲精品av一区二区| 免费看a级黄色片| bbb黄色大片| 亚洲一区二区三区色噜噜| 99riav亚洲国产免费| 亚洲av二区三区四区| 男女之事视频高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品自产拍在线观看55亚洲| 欧美日本视频| 亚洲性夜色夜夜综合| 国产伦精品一区二区三区四那| 美女高潮喷水抽搐中文字幕| 久久久久久久亚洲中文字幕 | 男人的好看免费观看在线视频| 一个人看的www免费观看视频| 69人妻影院| 99热精品在线国产| 97人妻精品一区二区三区麻豆| 色精品久久人妻99蜜桃| 精品一区二区三区av网在线观看| 国内精品美女久久久久久| 90打野战视频偷拍视频| 成人鲁丝片一二三区免费| 特级一级黄色大片| 男女下面进入的视频免费午夜| 桃红色精品国产亚洲av| 日本三级黄在线观看| 精品福利观看| 国产私拍福利视频在线观看| 久久久国产成人免费| 亚洲中文日韩欧美视频| 欧美在线黄色| 欧美成人一区二区免费高清观看| 国产高清三级在线| 亚洲精品一区av在线观看| 午夜精品在线福利| 99久久久亚洲精品蜜臀av| 久久午夜亚洲精品久久| 免费在线观看影片大全网站| 精品熟女少妇八av免费久了| 嫩草影视91久久| 狂野欧美白嫩少妇大欣赏| 久久久久久人人人人人| av女优亚洲男人天堂| 每晚都被弄得嗷嗷叫到高潮| 成人av在线播放网站| 悠悠久久av| 麻豆成人午夜福利视频| 伊人久久大香线蕉亚洲五| 一进一出抽搐动态| 午夜免费成人在线视频| 久久精品国产99精品国产亚洲性色| 久久人人精品亚洲av| 色视频www国产| 欧美日韩黄片免| 51午夜福利影视在线观看| 在线国产一区二区在线| 色综合站精品国产| 亚洲天堂国产精品一区在线| 中文字幕人成人乱码亚洲影| 日本三级黄在线观看| 国产精品久久久人人做人人爽| 国产高清videossex| 久久精品国产清高在天天线| 男女那种视频在线观看| 哪里可以看免费的av片| 中文字幕高清在线视频| 偷拍熟女少妇极品色| 国产精品99久久久久久久久| 欧美激情在线99| 手机成人av网站| 欧美黑人巨大hd| 麻豆成人av在线观看| 一边摸一边抽搐一进一小说| 动漫黄色视频在线观看| 国产欧美日韩一区二区精品| 嫩草影视91久久| 免费在线观看亚洲国产| 精华霜和精华液先用哪个| 午夜福利成人在线免费观看| 久久久成人免费电影| 制服人妻中文乱码| 免费高清视频大片| 国产精品野战在线观看| 欧美性猛交╳xxx乱大交人| 欧美乱色亚洲激情| 亚洲七黄色美女视频| 亚洲最大成人中文| 欧美成人一区二区免费高清观看| 一区二区三区激情视频| 99精品久久久久人妻精品| 又粗又爽又猛毛片免费看| 午夜激情欧美在线| 成人三级黄色视频| 欧美国产日韩亚洲一区| 99精品欧美一区二区三区四区| 看免费av毛片| 露出奶头的视频| 国产伦在线观看视频一区| 精品免费久久久久久久清纯| 国产精品99久久99久久久不卡| www.www免费av| 亚洲 欧美 日韩 在线 免费| 一进一出抽搐动态| 少妇人妻一区二区三区视频| 国产精品 欧美亚洲| 97超视频在线观看视频| 美女被艹到高潮喷水动态| 国产一区二区三区在线臀色熟女| 伊人久久精品亚洲午夜| 中亚洲国语对白在线视频| 日本黄色视频三级网站网址| 国产真实乱freesex| 此物有八面人人有两片| 欧美成人免费av一区二区三区| 久久婷婷人人爽人人干人人爱| 国产综合懂色| 琪琪午夜伦伦电影理论片6080| 91九色精品人成在线观看| 国产高清有码在线观看视频| 小说图片视频综合网站| 欧美日韩福利视频一区二区| 国产午夜精品久久久久久一区二区三区 | 亚洲成av人片在线播放无| 性色av乱码一区二区三区2| 日韩精品中文字幕看吧| 国产亚洲av嫩草精品影院| 国产一级毛片七仙女欲春2| 男人的好看免费观看在线视频| 在线天堂最新版资源| 久9热在线精品视频| 又粗又爽又猛毛片免费看| 国产精品亚洲av一区麻豆| 精华霜和精华液先用哪个| 久久久久精品国产欧美久久久| 波多野结衣高清作品| 中文字幕高清在线视频| 女人十人毛片免费观看3o分钟| 欧美区成人在线视频| 麻豆国产av国片精品| 亚洲真实伦在线观看| 欧美极品一区二区三区四区| 欧美黄色片欧美黄色片| 一级毛片高清免费大全| 亚洲七黄色美女视频| 国产视频一区二区在线看| 成人性生交大片免费视频hd| 91av网一区二区| 丰满人妻熟妇乱又伦精品不卡| 亚洲 国产 在线| 桃红色精品国产亚洲av| 久久欧美精品欧美久久欧美| 99热这里只有精品一区| 国产欧美日韩精品亚洲av| 中文资源天堂在线| 日韩av在线大香蕉| 亚洲av五月六月丁香网| 日韩av在线大香蕉| 成人一区二区视频在线观看| 又黄又爽又免费观看的视频| 欧美性感艳星| 在线观看舔阴道视频| 国产一区二区亚洲精品在线观看| 99热精品在线国产| 久久精品影院6| 午夜精品一区二区三区免费看| 国产高清视频在线观看网站| 欧美一级a爱片免费观看看| 亚洲18禁久久av| 亚洲在线自拍视频| 宅男免费午夜| 亚洲不卡免费看| 久久久成人免费电影| 日韩欧美国产在线观看| 丰满人妻一区二区三区视频av | 精品久久久久久久人妻蜜臀av| 国产高清有码在线观看视频| 色尼玛亚洲综合影院| 少妇人妻一区二区三区视频| 夜夜夜夜夜久久久久| 在线观看免费视频日本深夜| 亚洲黑人精品在线| 亚洲欧美日韩高清在线视频| 精品国产三级普通话版| 午夜福利成人在线免费观看| 色在线成人网| 很黄的视频免费| 嫩草影视91久久| 国产亚洲av嫩草精品影院| 国产精品亚洲一级av第二区| 嫩草影院入口| 黄色日韩在线| 俺也久久电影网| 成人特级黄色片久久久久久久| 日本 欧美在线| 波野结衣二区三区在线 | 精品电影一区二区在线| 一本久久中文字幕| 中文亚洲av片在线观看爽| a在线观看视频网站| 国产老妇女一区| 99国产精品一区二区蜜桃av| 久久久久国产精品人妻aⅴ院| 亚洲av一区综合| 精品久久久久久久人妻蜜臀av| 免费高清视频大片| 日韩大尺度精品在线看网址| 亚洲av成人av| 免费在线观看日本一区| 18禁国产床啪视频网站| 观看免费一级毛片| 在线免费观看的www视频| 国产精品久久电影中文字幕| 亚洲片人在线观看| svipshipincom国产片| 日韩成人在线观看一区二区三区| 国产69精品久久久久777片| 日韩免费av在线播放| 别揉我奶头~嗯~啊~动态视频| 国产爱豆传媒在线观看| 亚洲人成伊人成综合网2020| 久久久成人免费电影| 日韩欧美在线乱码| 99riav亚洲国产免费| 最好的美女福利视频网| 久久国产精品人妻蜜桃| 观看美女的网站| 国产一区二区三区视频了| 国产精品久久久久久久久免 | 香蕉丝袜av| 国产午夜福利久久久久久| www国产在线视频色| www.999成人在线观看| 欧美精品啪啪一区二区三区| 欧美日韩福利视频一区二区| 一a级毛片在线观看| 窝窝影院91人妻| 欧美成人a在线观看| 国产精品久久久久久亚洲av鲁大| 午夜久久久久精精品| 日韩欧美免费精品| 欧美日韩黄片免| 精品日产1卡2卡| 麻豆国产97在线/欧美| 一个人看的www免费观看视频| 非洲黑人性xxxx精品又粗又长| 午夜免费观看网址| 又黄又粗又硬又大视频| 亚洲国产色片| 婷婷六月久久综合丁香| 9191精品国产免费久久| 国产亚洲精品久久久com| 91av网一区二区| 窝窝影院91人妻| 精品久久久久久久末码| 免费看a级黄色片| 亚洲在线自拍视频| 99国产精品一区二区蜜桃av| 国产三级中文精品| 国产aⅴ精品一区二区三区波| 一个人免费在线观看电影| 99久久久亚洲精品蜜臀av| 国产久久久一区二区三区| 亚洲性夜色夜夜综合| 午夜福利在线观看免费完整高清在 | www国产在线视频色| 九色成人免费人妻av| 国产精品1区2区在线观看.| 国产视频内射| 波多野结衣高清作品| 老司机在亚洲福利影院| 欧美色欧美亚洲另类二区| 欧美黑人欧美精品刺激| 白带黄色成豆腐渣| 亚洲成a人片在线一区二区| 波多野结衣高清作品| h日本视频在线播放| 夜夜夜夜夜久久久久| 天堂√8在线中文| av福利片在线观看| 亚洲一区高清亚洲精品| 床上黄色一级片| av天堂中文字幕网| 久久精品国产99精品国产亚洲性色| 亚洲内射少妇av| 欧美区成人在线视频| 99久国产av精品| 露出奶头的视频| 亚洲中文字幕一区二区三区有码在线看| 精品99又大又爽又粗少妇毛片 | 老汉色av国产亚洲站长工具| 国语自产精品视频在线第100页| 成年免费大片在线观看| 变态另类丝袜制服| 国产精品亚洲一级av第二区| 婷婷丁香在线五月| 嫩草影院精品99|