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

    基于磁場(chǎng)FE和CFD的磁流變阻尼器力學(xué)性能分析

    2020-11-25 03:54:18劉昕運(yùn)吳大林馬吉?jiǎng)?/span>
    關(guān)鍵詞:磁通阻尼器活塞

    劉昕運(yùn) 吳大林 馬吉?jiǎng)?/p>

    (1.陸軍工程大學(xué)火炮工程系, 石家莊 050001; 2.西京大學(xué)機(jī)械工程學(xué)院, 西安 710123)

    0 引言

    磁流變(MR)阻尼器的內(nèi)部通常填充有特殊的磁流變液體,該液體由非導(dǎo)電液體、微小的鐵磁性粉末和少量添加劑混合而成[1]。這種特殊的流體在磁場(chǎng)作用下表現(xiàn)出粘塑性非牛頓流體的特性[2]。流體的屈服應(yīng)力由磁感應(yīng)強(qiáng)度決定,即為磁流變阻尼器的阻尼力可控的原因[3]。相比于膠質(zhì)阻尼器和其他阻尼器[4-5],MR阻尼器的滯回環(huán)范圍更大[6]。目前,MR阻尼器發(fā)展迅速,已在車輛工程[7-8]、土木工程[9]、醫(yī)療器械[10]、飛行器[11]和武器工業(yè)[12]等領(lǐng)域被廣泛使用。

    隨著計(jì)算機(jī)技術(shù)的發(fā)展,有限元分析(Finite element analysis,F(xiàn)EA)和計(jì)算流體動(dòng)力學(xué)(Computational fluid dynamics,CFD)仿真技術(shù)已成為分析磁流變液體和磁流變阻尼器的主要手段。與傳統(tǒng)的一維解析模型相比,F(xiàn)EA和CFD數(shù)值模擬更加準(zhǔn)確,且能夠獲得更多可視化參數(shù)[13]。文獻(xiàn)[14]在COMSOL中耦合了磁場(chǎng)和流體,并進(jìn)行了實(shí)驗(yàn)和仿真。文獻(xiàn)[15-16]在ANSYS/CFX中用CLL語言耦合了磁流變阻尼器模型。文獻(xiàn)[17]模擬了磁流變阻尼器的磁場(chǎng)分布,并考慮了溫度的影響。以上研究屬于邊界已知的“顯式”主動(dòng)運(yùn)動(dòng)問題,而本文研究的沖擊運(yùn)動(dòng)屬于邊界條件不可知的“隱式”被動(dòng)運(yùn)動(dòng)問題。目前,此類問題在國內(nèi)外尚未見相關(guān)研究和報(bào)道。

    本文提出單向耦合數(shù)值模擬方法。首先在COMSOL中對(duì)磁場(chǎng)進(jìn)行有限元分析,然后將磁場(chǎng)結(jié)果導(dǎo)入Fluent中,結(jié)合Six DOF模型和UDF文件對(duì)沖擊過程進(jìn)行分析,通過實(shí)驗(yàn)進(jìn)行驗(yàn)證,以期準(zhǔn)確獲取磁流變阻尼器的力學(xué)特性。

    1 數(shù)值模型

    1.1 結(jié)構(gòu)原理

    圖1 磁流變阻尼器CAD設(shè)計(jì)圖Fig.1 CAD design of MR damper

    近年來,磁流變阻尼器出現(xiàn)了多種結(jié)構(gòu),單出桿式[18]和雙出桿式[19],有氣體腔[20]和無氣體腔[21],特殊式[22]等。本文涉及的MR阻尼器用于緩沖艦炮系統(tǒng)中的高速運(yùn)動(dòng)的實(shí)驗(yàn)彈,CAD設(shè)計(jì)如圖1所示。由于外部幾何限制,對(duì)阻尼器結(jié)構(gòu)進(jìn)行了特殊設(shè)計(jì),如圖2所示,主要由活塞桿、磁流體活塞、空氣彈簧活塞、缸體、限位塊、線圈、壓縮空氣、磁流變液組成?;钊麠U兩端直徑相等,彈性回復(fù)力主要來自空氣壓縮?;钊欣@有兩組線圈,其被外殼密封。導(dǎo)線從活塞桿內(nèi)部通過并連接外部直流電源。當(dāng)線圈通電時(shí),將產(chǎn)生大致如圖2所示的磁場(chǎng)分布。磁場(chǎng)導(dǎo)致磁流變液中的鐵磁固體顆粒按磁感線形成鏈狀,導(dǎo)致環(huán)形間隙中的部分流體表現(xiàn)出極大的屈服應(yīng)力,流體通過間隙時(shí)會(huì)產(chǎn)生較大的阻尼力。阻尼器主要尺寸如表1所示。

    1.2 靜磁場(chǎng)FE模型

    有限元模型在多物理場(chǎng)仿真軟件COMSOL/Multiphysics中建立。依據(jù)麥克斯韋方程組[10]計(jì)算電磁場(chǎng),包括恒定電荷密度的連續(xù)性方程、安培定律和法拉第定律,即

    (1)

    其中

    B=μoμrH

    (2)

    式中H——磁場(chǎng)強(qiáng)度

    J——電流密度

    A——磁矢量勢(shì)

    B——磁通密度

    μr、μo——材料磁導(dǎo)率、真空磁導(dǎo)率

    根據(jù)阻尼器的對(duì)稱性,生成二維軸對(duì)稱網(wǎng)格,如圖3所示。計(jì)算域總共有70 375個(gè)四邊形網(wǎng)格單元。為了評(píng)估網(wǎng)格的相關(guān)性,生成另一組更加細(xì)密的網(wǎng)格,共281 500個(gè)四邊形單元。比較兩組網(wǎng)格的磁通強(qiáng)度計(jì)算結(jié)果,輸出曲線幾乎相同。這說明生成的網(wǎng)格已經(jīng)非常精細(xì),結(jié)果誤差與網(wǎng)格尺寸基本無關(guān)。

    圖3 磁流變阻尼器的有限元網(wǎng)格Fig.3 FE mesh of MR damper

    采用LORD公司的MRF-140CG型磁流變液,圖4為該磁流變液的特性曲線?;钊^和缸體采用碳素結(jié)構(gòu)鋼1045,活塞桿采用低磁導(dǎo)率金屬。線圈導(dǎo)線采用美國線規(guī)AWG 24標(biāo)準(zhǔn)銅線。部分材料參數(shù)如表2所示。

    圖4 MRF-140CG型磁流變液的特性曲線Fig.4 Characteristics of MRF-140CG MR fluid

    1.3 磁流變液CFD模型

    MR阻尼器的CFD分析需考慮湍流和能量傳遞,其中流體控制方程包括連續(xù)性方程、Navier-Stokes方程和能量方程[23]

    表2 材料屬性參數(shù)Tab.2 Material property parameters

    (3)

    式中ρ——流體密度V——流體速度

    fb——體積力p——差壓

    μ——?jiǎng)恿φ扯萮——比焓

    λ——熱導(dǎo)率Φ——粘性耗散項(xiàng)

    Sh——源項(xiàng)T——溫度

    1.3.1流體本構(gòu)模型

    根據(jù)圖4中磁流變液材料屬性定義磁通密度B和屈服應(yīng)力τy的關(guān)系曲線。使用Hill sigmoidal函數(shù)擬合為

    (4)

    式中,ε=117 024.8,β=0.967 66,k=1.693 63。

    目前,針對(duì)粘塑性流體提出了多種本構(gòu)模型[24-26]。本文采用由雙曲正切函數(shù)定義的流變粘度模型[13-14],即

    (5)

    式中α——粘度增長(zhǎng)控制因子,取0.03

    ζ——粘度增長(zhǎng)控制因子,取0.1

    對(duì)速度沖擊過程進(jìn)行CFD分析時(shí),必須考慮MR流體的可壓縮性。因?yàn)樽枘崞髟跊_擊瞬間壓力變化非常大,如果不考慮流體的壓縮性將導(dǎo)致壓力計(jì)算出現(xiàn)異常。液體壓力用Tait狀態(tài)方程描述,即

    (6)

    式中m——密度指數(shù)ρ0——參考密度

    E0——參考體積模量

    1.3.2空氣彈簧

    磁流變阻尼器有空氣彈簧結(jié)構(gòu),空氣腔存在一定初始?jí)毫Γ渲饕饔檬翘峁┗貜?fù)力,使被壓縮的阻尼器自動(dòng)復(fù)位??諝鈴椈闪s一般依據(jù)氣體多變過程方程來描述,即

    (7)

    式中p0——初始?jí)毫?,?.5 MPa

    Aair——活塞面積,取1 756.15 mm2

    W0——初始體積,取52 684.5 mm3

    n——?dú)怏w多變指數(shù),取1.3

    x——活塞位移

    1.3.3流體網(wǎng)格定義

    流體計(jì)算域仍然采用二維軸對(duì)稱四邊形網(wǎng)格。動(dòng)網(wǎng)格使用Layering層鋪技術(shù),網(wǎng)格被壓縮至一定程度時(shí),貼近壁面的網(wǎng)格會(huì)坍縮與合并。而網(wǎng)格被拉伸至一定程度時(shí),貼近壁面的網(wǎng)格會(huì)分裂出新層。為了順利實(shí)現(xiàn)網(wǎng)格層鋪法,在生成網(wǎng)格時(shí)需考慮Mesh Interface。如圖5所示,活塞兩側(cè)的壁面被定義為剛體移動(dòng)壁面。僅生成流體域網(wǎng)格,間隙處的網(wǎng)格被加密,總共包含了95 403個(gè)四邊形單元。

    圖5 流體網(wǎng)格定義示意圖Fig.5 Schematic of fluid mesh definition

    1.3.4用戶自定義函數(shù)

    磁流變阻尼器的活塞移動(dòng)時(shí),磁場(chǎng)的影響區(qū)域會(huì)跟著變化。所以,對(duì)于非牛頓區(qū)域的跟蹤和捕捉是必要的。與顯式邊界模型不同,隱式邊界模型的“磁場(chǎng)作用區(qū)域”不能直接定義,其需要依賴移動(dòng)壁面的實(shí)時(shí)坐標(biāo)。因此,基于C語言編寫了用戶自定義函數(shù)(UDF文件),并在求解器中進(jìn)行編譯。

    使用“DEFINE_PROPERTY”宏定義計(jì)算域的動(dòng)力粘度屬性。首先,使用“Get_Domain”宏和“Lookup_Thread”宏獲取移動(dòng)壁面的指針,在此之上用“F_CENTROID”宏訪問移動(dòng)壁面的坐標(biāo)信息并存入數(shù)組;然后寫入式(4)、(5)和磁場(chǎng)FE分析得到的磁通密度。以上函數(shù)中涉及2個(gè)變量,其中坐標(biāo)信息從數(shù)組中提取,剪切速率由“C_STRAIN_ RATE_MAG”宏訪問得到;最后用“C_ CENTROID”宏訪問計(jì)算域的所有網(wǎng)格單元,并把非牛頓粘度函數(shù)賦值于受磁場(chǎng)影響的單元。

    活塞只在x方向移動(dòng),而且存在空氣彈簧力和限位塊,故還需要定義6自由度模塊的UDF文件。根據(jù)頭文件six_dof.h,使用“DEFINE_SDOF_PROPERTIES”宏定義剛體屬性,代碼可以與屬性宏編寫到一個(gè)文件中。使用“SDOF_MASS”宏定義剛體的質(zhì)量,其值取阻尼器負(fù)載的質(zhì)量;使用“SDOFO_1DOF_T_P”宏定義自由度為1個(gè)平移自由度;使用“SDOFO_LOC, SDOFO_MIN,SDOFO_MAX”宏定義剛體的位移限制;最后用“SDOF_LOAD_F_X”宏寫入空氣彈簧力式(7)。

    2 數(shù)值結(jié)果

    對(duì)于單向耦合模型,式(4)、(5)是耦合磁場(chǎng)FEA和CFD的關(guān)鍵。通過靜磁場(chǎng)分析得到磁通密度B在流體域的分布,然后將其轉(zhuǎn)換為動(dòng)力粘度在流體域的分布,并導(dǎo)入CFD模型。

    2.1 靜磁場(chǎng)FE分析

    當(dāng)輸入電流I=4.0 A時(shí),計(jì)算域內(nèi)的磁通密度等值線如圖6所示。大部分磁感線從活塞的兩端穿過流體進(jìn)入阻尼器的缸體,由于流體和活塞桿的磁導(dǎo)率都非常低,磁通量更多地分布在磁性材料中。輸出環(huán)形間隙中間線上的磁通密度分布曲線如圖7所示,磁通密度存在兩個(gè)峰值,分別在環(huán)形間隙的兩端,最大值達(dá)到1.96 T。而在兩個(gè)峰值的中部以及外側(cè),磁通密度迅速降低至接近0。

    圖6 流體域中的磁通密度分布Fig.6 Magnetic flux density distribution in computational domain

    圖7 環(huán)形間隙中線上的磁通密度分布Fig.7 Magnetic flux density distribution on middle line of annular gap

    2.2 模型驗(yàn)證

    模型采用瞬態(tài)壓力基求解器(PBS)。使用Realizablek-ε模型,激活能量方程和粘性熱。設(shè)置Coupled算法并采用二階離散格式,求解控制Courant數(shù)為50,最大求解時(shí)間步長(zhǎng)為1×10-5s。

    磁通強(qiáng)度在流體中的主要作用區(qū)域如圖8所示。受磁場(chǎng)影響的主要區(qū)域不全在環(huán)形間隙中,間隙的兩端存在溢出??赡軙?huì)影響阻尼器的力學(xué)特性計(jì)算結(jié)果。因此,在UDF代碼中編寫環(huán)形間隙兩端的“圓弧”區(qū)域,其結(jié)合圖7的磁通密度分布曲線能夠更準(zhǔn)確地描述非牛頓粘度。

    圖8 環(huán)形間隙兩端的圓函數(shù)區(qū)Fig.8 Circular function area at both ends of annular gap

    在CSS-55100型萬能試驗(yàn)機(jī)上對(duì)阻尼器進(jìn)行勻速準(zhǔn)靜態(tài)運(yùn)動(dòng)實(shí)驗(yàn),如圖9所示。固定缸體,對(duì)活塞桿進(jìn)行勻速的加載和卸載。利用壓力傳感器和位移傳感器采集阻力和位移。圖10、11分別為勻速準(zhǔn)靜態(tài)運(yùn)動(dòng)的阻抗力-時(shí)間曲線和滯回曲線。對(duì)比阻尼器的實(shí)驗(yàn)和仿真曲線,驗(yàn)證了模型的準(zhǔn)確性。

    圖9 勻速準(zhǔn)靜態(tài)實(shí)驗(yàn)設(shè)備Fig.9 Uniform quasi-static experimental equipment

    圖10 勻速準(zhǔn)靜態(tài)運(yùn)動(dòng)的阻抗力-時(shí)間曲線Fig.10 Resistance-time curves of uniform quasi-static motion

    圖11 勻速準(zhǔn)靜態(tài)運(yùn)動(dòng)的滯回曲線Fig.11 Hysteresis loop of uniform quasi-static motion

    2.3 沖擊CFD分析

    當(dāng)剛體初速度為7.4 m/s,壁面邊界的運(yùn)動(dòng)狀態(tài)由流固雙向耦合的6自由度模塊自動(dòng)解算,總仿真時(shí)間為80 ms時(shí),整個(gè)沖擊過程的等值線和流線圖如圖12所示。沖擊初期,流體在環(huán)形間隙中瞬時(shí)最大速度達(dá)到了75 m/s以上,但隨著動(dòng)能的耗散,流體速度會(huì)很快降低,達(dá)到最大位移僅需4 ms。磁流變阻尼器在運(yùn)動(dòng)過程中,流體產(chǎn)生了許多不停變化的渦,且流線極為不規(guī)則。但較大的渦都出現(xiàn)在活塞運(yùn)動(dòng)方向的后面。

    阻尼器吸收的動(dòng)能主要轉(zhuǎn)化為流體熱能。由于液體具有粘性,內(nèi)部摩擦將產(chǎn)生熱量,初始溫度為300 K時(shí),流體溫度分布如圖13所示。高溫液體主要產(chǎn)生于環(huán)形間隙,因?yàn)榇颂幖羟兴俾屎蛣?dòng)力粘度最大。沖擊運(yùn)動(dòng)結(jié)束時(shí)溫度分布并不均勻,活塞一側(cè)的平均溫度暫時(shí)高于另一側(cè)。

    輸出的磁流變阻尼器運(yùn)動(dòng)學(xué)特性曲線如圖14所示。MR阻尼器的壓縮過程持續(xù)4 ms,拉伸過程持續(xù)大約70 ms。圖14a顯示,MR阻尼器最終無法完全復(fù)位,在位移為1.3 mm處幾乎停止運(yùn)動(dòng),這是由Bingham流體性質(zhì)決定的。阻尼器在復(fù)原運(yùn)動(dòng)時(shí),流體域的剪切速率逐漸降低。由式(5)可知,在剪切速率很低時(shí),表觀粘度會(huì)逐漸大幅升高,從而導(dǎo)致流體剪切應(yīng)力變大。而隨著位移行程的減小,空氣彈簧力逐漸降低。當(dāng)空氣彈簧力小于流體阻力時(shí),活塞開始減速直到停止運(yùn)動(dòng)。阻尼器在復(fù)原運(yùn)動(dòng)階段先加速后減速的特征如圖14b所示。在實(shí)際應(yīng)用中,這種現(xiàn)象并不會(huì)影響磁流體阻尼器的使用。

    圖12 磁流變阻尼器速度等值線和流線Fig.12 Velocity contours and streamtraces of MR damper

    圖13 磁流變阻尼器溫度等值線Fig.13 Temperature contours of MR damper

    圖14 沖擊運(yùn)動(dòng)的運(yùn)動(dòng)學(xué)特性曲線Fig.14 Kinematic characteristic of impact motion

    這是由于當(dāng)控制電源斷開后,磁流變液會(huì)恢復(fù)牛頓流體特性,阻尼器將很快完全復(fù)位。圖15為沖擊運(yùn)動(dòng)的滯回曲線,阻尼器的阻抗力曲線變化非常不規(guī)則,最大阻抗力出現(xiàn)在沖擊初期,在行程2 mm處達(dá)到最大阻抗力110 kN,而后迅速降低并逐漸趨于平穩(wěn)。復(fù)原階段中的阻抗力保持在200~400 N極低的水平。該磁流體阻尼器在沖擊過程中幾乎吸收了所有的初始動(dòng)能,其緩沖效能較為優(yōu)異。

    圖15 沖擊運(yùn)動(dòng)的滯回曲線Fig.15 Hysteresis loop of impact motion

    3 結(jié)論

    (1)針對(duì)某磁流變阻尼器提出了結(jié)合磁路有限元分析和計(jì)算流體動(dòng)力學(xué)分析的單向耦合仿真研究方法。在實(shí)驗(yàn)驗(yàn)證勻速準(zhǔn)靜態(tài)運(yùn)動(dòng)模型的基礎(chǔ)上,對(duì)模型進(jìn)行了特定沖擊載荷的流固雙向耦合模擬,獲得了其力學(xué)特性規(guī)律。

    (2)在隱式邊界條件的CFD分析中,解決了粘塑性流體本構(gòu)定義、非牛頓粘度區(qū)的跟蹤和捕獲、以及被動(dòng)運(yùn)動(dòng)的動(dòng)網(wǎng)格處理等難點(diǎn)。在磁路FE分析得到雙峰分布的磁通密度的基礎(chǔ)上,用雙曲正切函數(shù)定義了粘塑性磁流變液的非牛頓粘度,并用“圓函數(shù)”準(zhǔn)確捕獲非牛頓粘度區(qū)。模型也考慮了液體的可壓縮性和多變方程描述的空氣彈簧力。

    (3)在沖擊運(yùn)動(dòng)過程中,流體域產(chǎn)生大小不一的多個(gè)渦,其中最明顯的渦產(chǎn)生于活塞運(yùn)動(dòng)方向的后側(cè)。磁流變液中的熱量主要產(chǎn)生于環(huán)形間隙,最終在活塞一側(cè)的平均溫度高于另一側(cè)。在持續(xù)供電的情況下,阻尼器將無法完全復(fù)位,但電路控制可以解決這一問題。

    猜你喜歡
    磁通阻尼器活塞
    核電廠機(jī)械式阻尼器故障分析及處理
    軸向磁通電勵(lì)磁雙凸極電機(jī)及容錯(cuò)運(yùn)行控制策略
    一種活塞擠壓式室內(nèi)通風(fēng)設(shè)備的研發(fā)與應(yīng)用
    幾種軟鋼阻尼器對(duì)建筑安全的影響
    連梁阻尼器的初步設(shè)計(jì)方法
    四川建筑(2018年4期)2018-09-14 00:16:52
    基于LabVIEW的微型磁通門磁強(qiáng)計(jì)測(cè)試系統(tǒng)搭建
    基于磁通門原理的零磁通交直流電流傳感器
    基于FPGA的數(shù)字磁通計(jì)設(shè)計(jì)
    面向600MW汽輪發(fā)電機(jī)的顆粒阻尼器減振試驗(yàn)
    KS Kolbenschmidt公司的新型鋼活塞
    国产中年淑女户外野战色| 性欧美人与动物交配| 丰满人妻熟妇乱又伦精品不卡| 欧美乱码精品一区二区三区| 好男人电影高清在线观看| 黄色女人牲交| 18美女黄网站色大片免费观看| 97超级碰碰碰精品色视频在线观看| 黄色丝袜av网址大全| 亚洲国产日韩欧美精品在线观看 | 亚洲av一区综合| 亚洲欧美日韩高清专用| 精品久久久久久久久久久久久| 国产午夜精品论理片| 欧美乱码精品一区二区三区| 在线十欧美十亚洲十日本专区| h日本视频在线播放| 成人特级黄色片久久久久久久| 日韩中文字幕欧美一区二区| 高清日韩中文字幕在线| 一个人观看的视频www高清免费观看| 国内精品久久久久久久电影| 国产亚洲精品av在线| ponron亚洲| 美女黄网站色视频| 国产精品嫩草影院av在线观看 | 国产精品爽爽va在线观看网站| 久久精品91无色码中文字幕| 国产精品永久免费网站| 在线a可以看的网站| 欧美bdsm另类| 亚洲av一区综合| 老司机午夜十八禁免费视频| 丁香六月欧美| 久久久久久人人人人人| 欧美性猛交黑人性爽| 一级作爱视频免费观看| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 老汉色av国产亚洲站长工具| 亚洲人与动物交配视频| 日日摸夜夜添夜夜添小说| 国内精品美女久久久久久| 国产精品免费一区二区三区在线| 禁无遮挡网站| 日日摸夜夜添夜夜添小说| 国产精品综合久久久久久久免费| 十八禁网站免费在线| 久久性视频一级片| 国产av麻豆久久久久久久| 一区二区三区高清视频在线| 90打野战视频偷拍视频| 在线看三级毛片| 免费av观看视频| 悠悠久久av| 久久香蕉国产精品| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品综合一区在线观看| 美女大奶头视频| 观看美女的网站| 国产69精品久久久久777片| 久久精品国产综合久久久| 亚洲 国产 在线| 国产三级黄色录像| 麻豆成人午夜福利视频| 真人做人爱边吃奶动态| 色在线成人网| 国产精品久久久久久人妻精品电影| 成熟少妇高潮喷水视频| 一个人看的www免费观看视频| 午夜福利欧美成人| 一边摸一边抽搐一进一小说| 亚洲天堂国产精品一区在线| 悠悠久久av| 亚洲精品粉嫩美女一区| 色视频www国产| 中文资源天堂在线| 一区福利在线观看| 给我免费播放毛片高清在线观看| 日日摸夜夜添夜夜添小说| 欧美+亚洲+日韩+国产| 国产亚洲欧美在线一区二区| 欧美中文日本在线观看视频| 免费av不卡在线播放| 香蕉av资源在线| 日本一二三区视频观看| 久久久色成人| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 精品久久久久久久毛片微露脸| 国内精品久久久久久久电影| 免费无遮挡裸体视频| 婷婷六月久久综合丁香| 欧美zozozo另类| 亚洲avbb在线观看| 亚洲专区国产一区二区| 麻豆国产97在线/欧美| 熟女人妻精品中文字幕| 欧美乱色亚洲激情| 中文字幕人成人乱码亚洲影| 给我免费播放毛片高清在线观看| 久久久久久久久中文| 国产精品一区二区三区四区久久| 很黄的视频免费| 欧美绝顶高潮抽搐喷水| eeuss影院久久| 嫩草影院入口| 免费人成在线观看视频色| 大型黄色视频在线免费观看| 日日干狠狠操夜夜爽| www.熟女人妻精品国产| 国产精品综合久久久久久久免费| 观看美女的网站| 精品久久久久久成人av| 亚洲成人免费电影在线观看| АⅤ资源中文在线天堂| xxxwww97欧美| 亚洲一区二区三区不卡视频| 免费av毛片视频| 禁无遮挡网站| 国产精品99久久久久久久久| 欧美黄色片欧美黄色片| 国产欧美日韩一区二区精品| 成人永久免费在线观看视频| 免费电影在线观看免费观看| 午夜激情欧美在线| 亚洲最大成人中文| 久久精品国产清高在天天线| 精品熟女少妇八av免费久了| 日韩精品中文字幕看吧| 变态另类丝袜制服| 国产亚洲欧美98| 亚洲色图av天堂| 精品乱码久久久久久99久播| 禁无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 丰满人妻一区二区三区视频av | 99久久九九国产精品国产免费| 18禁国产床啪视频网站| 在线免费观看不下载黄p国产 | 麻豆一二三区av精品| 午夜影院日韩av| 男女视频在线观看网站免费| 最新美女视频免费是黄的| 日本免费一区二区三区高清不卡| 欧美zozozo另类| xxx96com| 欧美日本亚洲视频在线播放| 禁无遮挡网站| 中文资源天堂在线| 国产精品一及| 免费在线观看成人毛片| 精品不卡国产一区二区三区| 成人特级av手机在线观看| 性欧美人与动物交配| 波野结衣二区三区在线 | 欧美国产日韩亚洲一区| 悠悠久久av| 香蕉丝袜av| 在线播放国产精品三级| 99国产极品粉嫩在线观看| 岛国在线免费视频观看| 真实男女啪啪啪动态图| 成人三级黄色视频| 1024手机看黄色片| 色老头精品视频在线观看| 亚洲,欧美精品.| 免费在线观看亚洲国产| 亚洲 欧美 日韩 在线 免费| 日韩欧美免费精品| 久久久精品大字幕| 美女大奶头视频| 亚洲五月婷婷丁香| 亚洲一区二区三区不卡视频| 久久久久久久久大av| 黄色片一级片一级黄色片| 在线十欧美十亚洲十日本专区| 久久中文看片网| 国产中年淑女户外野战色| 国产乱人伦免费视频| 九色成人免费人妻av| 性欧美人与动物交配| 熟女电影av网| 舔av片在线| 成人18禁在线播放| 日韩免费av在线播放| 51午夜福利影视在线观看| 最新美女视频免费是黄的| 亚洲精品在线观看二区| 免费搜索国产男女视频| netflix在线观看网站| 丁香欧美五月| 亚洲人与动物交配视频| 色老头精品视频在线观看| 2021天堂中文幕一二区在线观| 久久亚洲真实| 久久国产精品人妻蜜桃| 一区二区三区激情视频| 最近最新中文字幕大全免费视频| 99久久精品一区二区三区| 国产精品香港三级国产av潘金莲| 可以在线观看的亚洲视频| 欧美日韩瑟瑟在线播放| 中文字幕熟女人妻在线| 国产高清videossex| 国产三级中文精品| 久久午夜亚洲精品久久| 亚洲av二区三区四区| 丰满乱子伦码专区| 久久久久精品国产欧美久久久| 成人特级黄色片久久久久久久| 99国产综合亚洲精品| xxx96com| 久久婷婷人人爽人人干人人爱| 中文字幕精品亚洲无线码一区| 免费无遮挡裸体视频| 国产精品日韩av在线免费观看| 亚洲人成网站在线播| 成人av一区二区三区在线看| 亚洲精品久久国产高清桃花| 久久精品国产亚洲av涩爱 | 欧美乱妇无乱码| 国产伦精品一区二区三区四那| 日本五十路高清| 亚洲aⅴ乱码一区二区在线播放| 国产探花在线观看一区二区| 国产精品久久久人人做人人爽| x7x7x7水蜜桃| 免费在线观看影片大全网站| 亚洲精品久久国产高清桃花| 欧美激情久久久久久爽电影| av在线天堂中文字幕| 真实男女啪啪啪动态图| 操出白浆在线播放| 国产熟女xx| 久久久成人免费电影| 亚洲人成伊人成综合网2020| 日韩欧美 国产精品| av视频在线观看入口| 窝窝影院91人妻| 国产免费av片在线观看野外av| x7x7x7水蜜桃| 老熟妇仑乱视频hdxx| 中亚洲国语对白在线视频| 欧美乱码精品一区二区三区| 国产国拍精品亚洲av在线观看 | 神马国产精品三级电影在线观看| 日韩欧美精品免费久久 | 成人三级黄色视频| 免费在线观看成人毛片| 日韩免费av在线播放| 成人无遮挡网站| 91九色精品人成在线观看| 五月玫瑰六月丁香| 中文在线观看免费www的网站| 午夜福利18| 午夜精品久久久久久毛片777| 婷婷亚洲欧美| 精品免费久久久久久久清纯| 欧美日本视频| 床上黄色一级片| 国产在视频线在精品| 无遮挡黄片免费观看| 99在线视频只有这里精品首页| 老熟妇乱子伦视频在线观看| 18禁美女被吸乳视频| 99久久精品一区二区三区| 国产91精品成人一区二区三区| aaaaa片日本免费| 亚洲 欧美 日韩 在线 免费| 在线观看午夜福利视频| 一本综合久久免费| 九九热线精品视视频播放| 一本精品99久久精品77| 久久这里只有精品中国| 国产97色在线日韩免费| 97碰自拍视频| 亚洲人与动物交配视频| 久久久久久九九精品二区国产| 欧美成人一区二区免费高清观看| 又黄又爽又免费观看的视频| 亚洲av第一区精品v没综合| 国产av麻豆久久久久久久| www日本黄色视频网| 国产精品美女特级片免费视频播放器| 欧美三级亚洲精品| 久久久久国内视频| 国产黄片美女视频| 夜夜爽天天搞| 国产精品一区二区三区四区免费观看 | 亚洲av熟女| 成人精品一区二区免费| 成年女人永久免费观看视频| 天堂网av新在线| 丰满的人妻完整版| 久久久久精品国产欧美久久久| 97人妻精品一区二区三区麻豆| 精品无人区乱码1区二区| 老汉色∧v一级毛片| 中亚洲国语对白在线视频| 伊人久久精品亚洲午夜| 一个人免费在线观看电影| 麻豆一二三区av精品| www.999成人在线观看| 国产成人啪精品午夜网站| 成人三级黄色视频| 国产精品久久久久久久久免 | 亚洲美女视频黄频| 日韩欧美在线二视频| 色在线成人网| 久久精品国产99精品国产亚洲性色| 日韩免费av在线播放| 国产精品嫩草影院av在线观看 | 国产色爽女视频免费观看| avwww免费| 色视频www国产| 嫩草影视91久久| 真人一进一出gif抽搐免费| 99视频精品全部免费 在线| 老鸭窝网址在线观看| 亚洲欧美日韩东京热| 舔av片在线| 真人做人爱边吃奶动态| 亚洲av日韩精品久久久久久密| 国产精品一区二区免费欧美| 宅男免费午夜| 中文在线观看免费www的网站| 搡女人真爽免费视频火全软件 | av中文乱码字幕在线| 99精品久久久久人妻精品| 成年女人看的毛片在线观看| www.www免费av| 黄色日韩在线| 亚洲精品影视一区二区三区av| 欧美绝顶高潮抽搐喷水| 日韩欧美免费精品| 日本熟妇午夜| 69av精品久久久久久| 一个人看视频在线观看www免费 | 国产精品久久久久久久电影 | 国产高潮美女av| 禁无遮挡网站| 精品国产亚洲在线| 国产欧美日韩一区二区精品| 欧美一区二区精品小视频在线| 成年人黄色毛片网站| 亚洲电影在线观看av| www.www免费av| 又紧又爽又黄一区二区| 日本一二三区视频观看| 麻豆久久精品国产亚洲av| 久久性视频一级片| 两个人的视频大全免费| 一进一出抽搐动态| 97碰自拍视频| 国产69精品久久久久777片| 男女午夜视频在线观看| 老司机午夜福利在线观看视频| 久久人人精品亚洲av| 免费在线观看成人毛片| 精品久久久久久,| 首页视频小说图片口味搜索| 一级a爱片免费观看的视频| 狂野欧美激情性xxxx| 亚洲乱码一区二区免费版| 两人在一起打扑克的视频| 国产男靠女视频免费网站| 亚洲欧美激情综合另类| 日本与韩国留学比较| 国产97色在线日韩免费| www日本黄色视频网| 精品日产1卡2卡| 91av网一区二区| 久久精品国产99精品国产亚洲性色| 亚洲成人精品中文字幕电影| 精品日产1卡2卡| 在线观看一区二区三区| 久久伊人香网站| 久久久成人免费电影| 在线播放国产精品三级| 美女高潮喷水抽搐中文字幕| 日日夜夜操网爽| 琪琪午夜伦伦电影理论片6080| 丁香六月欧美| av在线天堂中文字幕| 香蕉久久夜色| 欧美在线一区亚洲| 久久香蕉精品热| 高清毛片免费观看视频网站| 国产免费一级a男人的天堂| 黄片大片在线免费观看| 国产免费一级a男人的天堂| 国产真实乱freesex| 亚洲专区中文字幕在线| 色视频www国产| 成人三级黄色视频| 久久伊人香网站| 亚洲激情在线av| 最近最新中文字幕大全电影3| 宅男免费午夜| 窝窝影院91人妻| 淫秽高清视频在线观看| 国产探花极品一区二区| 亚洲精华国产精华精| 国产视频一区二区在线看| 禁无遮挡网站| 亚洲精品亚洲一区二区| xxxwww97欧美| 在线观看一区二区三区| 国产又黄又爽又无遮挡在线| 波多野结衣高清作品| 美女被艹到高潮喷水动态| 90打野战视频偷拍视频| 天天一区二区日本电影三级| 欧美xxxx黑人xx丫x性爽| 全区人妻精品视频| 国产淫片久久久久久久久 | www日本黄色视频网| 国产精品 欧美亚洲| 久久中文看片网| 一级毛片女人18水好多| 亚洲性夜色夜夜综合| 18+在线观看网站| 观看美女的网站| 老熟妇乱子伦视频在线观看| 美女cb高潮喷水在线观看| 1000部很黄的大片| 亚洲欧美激情综合另类| 国产不卡一卡二| 国产久久久一区二区三区| 一进一出好大好爽视频| 美女被艹到高潮喷水动态| 中文字幕熟女人妻在线| 少妇高潮的动态图| 青草久久国产| 亚洲 国产 在线| 首页视频小说图片口味搜索| 久久这里只有精品中国| 怎么达到女性高潮| 国产一区二区三区在线臀色熟女| 成人高潮视频无遮挡免费网站| 色精品久久人妻99蜜桃| 99久久久亚洲精品蜜臀av| 亚洲人成伊人成综合网2020| 少妇丰满av| 黄色丝袜av网址大全| 亚洲精品美女久久久久99蜜臀| 丁香六月欧美| 日本a在线网址| 国产乱人伦免费视频| 禁无遮挡网站| av天堂在线播放| 中文字幕精品亚洲无线码一区| 欧美av亚洲av综合av国产av| 国产av一区在线观看免费| 十八禁网站免费在线| 国产亚洲精品综合一区在线观看| 欧美一区二区亚洲| 国产乱人伦免费视频| 亚洲成av人片免费观看| 久久久久亚洲av毛片大全| 蜜桃亚洲精品一区二区三区| 亚洲激情在线av| 国产av麻豆久久久久久久| 日韩欧美免费精品| 国产视频一区二区在线看| 精品熟女少妇八av免费久了| 亚洲精品日韩av片在线观看 | 久久国产乱子伦精品免费另类| 成熟少妇高潮喷水视频| av国产免费在线观看| 日韩人妻高清精品专区| 国产精品美女特级片免费视频播放器| 18禁黄网站禁片免费观看直播| 黄色丝袜av网址大全| 床上黄色一级片| 搞女人的毛片| 俺也久久电影网| 又粗又爽又猛毛片免费看| 特大巨黑吊av在线直播| 国产精品免费一区二区三区在线| 小蜜桃在线观看免费完整版高清| 亚洲aⅴ乱码一区二区在线播放| 久久精品91蜜桃| 狂野欧美白嫩少妇大欣赏| 91av网一区二区| 制服人妻中文乱码| 久久久国产精品麻豆| 又黄又爽又免费观看的视频| 真人一进一出gif抽搐免费| 欧美日本亚洲视频在线播放| 国产在视频线在精品| 欧美一区二区国产精品久久精品| 欧美不卡视频在线免费观看| 久久午夜亚洲精品久久| 欧美av亚洲av综合av国产av| 国产久久久一区二区三区| 九色成人免费人妻av| 日韩国内少妇激情av| 日本一二三区视频观看| 国产野战对白在线观看| 尤物成人国产欧美一区二区三区| 综合色av麻豆| av福利片在线观看| bbb黄色大片| avwww免费| 啦啦啦韩国在线观看视频| 一夜夜www| 成人特级黄色片久久久久久久| 国内精品久久久久久久电影| 丰满的人妻完整版| 亚洲国产欧美网| 久久天躁狠狠躁夜夜2o2o| 久久久久久九九精品二区国产| 欧美成人免费av一区二区三区| 久久久久久久久大av| 丁香六月欧美| 亚洲黑人精品在线| 亚洲精品在线美女| 欧美3d第一页| 亚洲av免费在线观看| 国产亚洲精品av在线| 免费看十八禁软件| 国产av在哪里看| av黄色大香蕉| 天天添夜夜摸| 夜夜看夜夜爽夜夜摸| xxxwww97欧美| 97超视频在线观看视频| 91九色精品人成在线观看| 18禁裸乳无遮挡免费网站照片| av天堂在线播放| 国产在视频线在精品| 久久久精品欧美日韩精品| 欧美日韩综合久久久久久 | 亚洲18禁久久av| 欧美最黄视频在线播放免费| 久久久久久久精品吃奶| 在线观看av片永久免费下载| 日本在线视频免费播放| 中文字幕熟女人妻在线| 夜夜看夜夜爽夜夜摸| 性欧美人与动物交配| 国产高清视频在线观看网站| 日本熟妇午夜| 精品熟女少妇八av免费久了| 国产主播在线观看一区二区| 中文字幕高清在线视频| 午夜老司机福利剧场| 日本与韩国留学比较| 久久久久久久亚洲中文字幕 | 午夜福利成人在线免费观看| 国产精品三级大全| 欧美激情久久久久久爽电影| 亚洲成人久久爱视频| 九九久久精品国产亚洲av麻豆| 日日摸夜夜添夜夜添小说| 国产精品自产拍在线观看55亚洲| 亚洲最大成人中文| 久久性视频一级片| av国产免费在线观看| 国产成人a区在线观看| 亚洲国产欧美人成| 无限看片的www在线观看| 高清毛片免费观看视频网站| 在线视频色国产色| 夜夜躁狠狠躁天天躁| 丝袜美腿在线中文| 此物有八面人人有两片| 久久国产精品人妻蜜桃| 精品免费久久久久久久清纯| 久久久久久九九精品二区国产| 亚洲国产精品久久男人天堂| 国产精品自产拍在线观看55亚洲| 精品国产亚洲在线| 精品日产1卡2卡| 欧美一区二区国产精品久久精品| 欧美日韩综合久久久久久 | 深夜精品福利| 少妇的逼水好多| 一区福利在线观看| 日韩欧美免费精品| 网址你懂的国产日韩在线| 日本一二三区视频观看| 国产乱人视频| 中出人妻视频一区二区| 国产欧美日韩一区二区三| 日韩欧美国产在线观看| 国产精品一区二区三区四区久久| 少妇的逼水好多| 18禁美女被吸乳视频| 国产精品1区2区在线观看.| 久久午夜亚洲精品久久| 亚洲人成网站在线播放欧美日韩| 国产亚洲精品久久久com| ponron亚洲| 桃色一区二区三区在线观看| 村上凉子中文字幕在线| 国产精品久久久久久久久免 | 久久久久免费精品人妻一区二区| 男女之事视频高清在线观看| 国产熟女xx| 欧美激情在线99| 女同久久另类99精品国产91| 欧美色视频一区免费| 91久久精品电影网| 国产成人av教育| 搡老岳熟女国产| 高清在线国产一区| 哪里可以看免费的av片| 亚洲成av人片在线播放无| 久久久久久久精品吃奶| 国产精品久久久久久久电影 | 日本三级黄在线观看| 亚洲熟妇熟女久久|