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

    混合流動式磁流變閥結(jié)構(gòu)設(shè)計與壓降性能分析

    2016-10-27 02:04:22胡國良廖明科龔美慧
    農(nóng)業(yè)機械學報 2016年9期
    關(guān)鍵詞:圓盤圓環(huán)阻尼

    胡國良 鐘 芳 廖明科 龔美慧

    (華東交通大學載運工具與裝備教育部重點實驗室, 南昌 330013)

    ?

    混合流動式磁流變閥結(jié)構(gòu)設(shè)計與壓降性能分析

    胡國良鐘芳廖明科龔美慧

    (華東交通大學載運工具與裝備教育部重點實驗室, 南昌 330013)

    針對目前磁流變閥結(jié)構(gòu)單一且體積大的不足,設(shè)計了一種結(jié)構(gòu)緊湊的混合流動式磁流變閥,該磁流變閥阻尼間隙液流通道由單個軸向圓環(huán)流動、單個徑向圓盤流動和單個中心小孔流動串聯(lián)組合而成。闡述了混合流動式磁流變閥結(jié)構(gòu)及工作原理;分析了勵磁線圈作用下有效阻尼間隙處的磁路分布,同時建立了其壓降數(shù)學模型。采用有限元法對混合流動式磁流變閥電磁場進行了建模仿真,對阻尼間隙處的磁場強度和剪切屈服應力的分布規(guī)律進行了分析;仿真結(jié)果表明軸向圓環(huán)阻尼間隙厚度為1.0 mm,徑向圓盤阻尼間隙厚度為0.5 mm,加載電流為1.2 A時磁流變閥壓降最大,為3 342 kPa。對混合流動式磁流變閥壓降性能進行了試驗測試,具體分析了加載電流以及徑向圓盤阻尼間隙厚度對磁流變閥壓降的影響,試驗結(jié)果表明軸向圓環(huán)阻尼間隙厚度為1.0 mm,徑向圓盤阻尼間隙厚度為0.5 mm,加載電流為1.2 A時磁流變閥壓降最大,為2 650 kPa,試驗結(jié)果與仿真結(jié)果變化趨勢基本一致。

    磁流變閥; 混合流動; 壓降; 有限元分析

    引言

    在磁場作用下,磁流變液能在毫秒級時間內(nèi)從流動狀態(tài)瞬變?yōu)榘牍虘B(tài),并且這種變化可逆[1]。借助于磁流變液的流變特性,以磁流變液為工作介質(zhì)的磁流變智能器件目前已在汽車交通、橋梁建筑等行業(yè)得到廣泛應用[2-3]。

    磁流變閥是一種典型的磁流變智能器件,它可作為旁通閥用來控制阻尼器,應用于汽車半主動懸架系統(tǒng)[4];也可內(nèi)置于磁流變懸置里面,用來控制液流通道的流速,應用于海洋平臺的減振系統(tǒng)[5-6]。由于磁流變閥的閥芯和閥體相對固定,只需控制激勵電流,改變阻尼間隙處的磁感應強度,就可獲得可控的進出口壓降,故操作方便且響應速度快。相比傳統(tǒng)液壓控制閥,磁流變閥機械結(jié)構(gòu)簡單、響應速度快、運動平穩(wěn)且系統(tǒng)綜合成本低[7-8]。

    國內(nèi)外學者對磁流變閥的研究主要集中在結(jié)構(gòu)設(shè)計及參數(shù)優(yōu)化兩方面。GORODKIN等[9]首次提出了一種用于被動阻尼系統(tǒng)中的磁流變節(jié)流閥,并進行了試驗分析。SALLOOM等[10-11]提出了一種由4個圓環(huán)形單線圈磁流變閥單元組合而成的三位四通方向控制閥,該磁流變閥通過控制流體流動來控制系統(tǒng)速度及位置。ICHWAN等[12]提出了一種模塊式磁流變閥,仿真和試驗結(jié)果表明,通過把2個或多個相同的磁流變閥單元串聯(lián)組合在一起,可產(chǎn)生較大的壓降。趙燦等[13]在磁流變閥入口及出口處將液流通道設(shè)計為非直線型,這種結(jié)構(gòu)設(shè)計可有效增大磁流變閥壓降調(diào)節(jié)范圍。本課題組也先后設(shè)計了一種采用雙線圈進行激勵的圓環(huán)型磁流變閥以及一種阻尼間隙可在1~2 mm之間機械無級可調(diào)的磁流變閥[14-15];同時對比分析了軸向圓環(huán)流磁流變閥和徑向圓盤流磁流變閥的壓降變化規(guī)律[16-17]。此外,NGUYEN等[18-20]為了得到最佳磁流變閥結(jié)構(gòu),采用ANSYS有限元分析方法分別對單線圈、雙線圈、三線圈以及圓環(huán)圓盤組合型4種不同結(jié)構(gòu)的磁流變閥參數(shù)進行了仿真優(yōu)化分析。

    上述研究中磁流變閥液流通道結(jié)構(gòu)主要是單一軸向圓環(huán)流動式或單一徑向圓盤流動式;在不增加閥外形結(jié)構(gòu)尺寸前提下,這兩種液流通道結(jié)構(gòu)的磁流變閥進出口壓降均較小,壓降調(diào)節(jié)范圍并沒有取得理想效果?;诖?,國內(nèi)外學者提出把軸向圓環(huán)流和徑向圓盤流兩種液流通道組合在一起,達到延長有限阻尼間隙長度,增加進出口壓降的目的[21-24]。WANG等[21-22]首先提出并設(shè)計了一種同時具有軸向圓環(huán)流和徑向圓盤流液流通道的磁流變閥,同等磁場強度下能有效提高磁流變閥壓降。但該磁流變閥液流通道由2個軸向圓環(huán)流、2個徑向圓盤流以及1個中心小孔流組成,并且設(shè)計時假設(shè)軸向圓環(huán)阻尼間隙和徑向圓盤阻尼間隙處的磁場強度相等,導致壓降沒有得到優(yōu)化。在此基礎(chǔ)上,IMADUDDIN等[23-24]提出了一種集成多個圓環(huán)流和圓盤流的蜿蜒式磁流變閥,該磁流變閥液流通道由5個圓環(huán)流、6個圓盤流和2個中心小孔流按順序串聯(lián)組合而成,仿真和試驗表明磁流變閥能產(chǎn)生大于2.5 MPa的壓降。但該閥液流通道結(jié)構(gòu)復雜,很容易導致閥工作時堵塞;另外壓降主要由徑向圓盤阻尼間隙提供,實際上是一個徑向流磁流變閥。

    本文提出并設(shè)計一種結(jié)構(gòu)緊湊的混合流動式磁流變閥,該磁流變閥阻尼間隙液流通道由單個軸向圓環(huán)流動、單個徑向圓盤流動和單個中心小孔流動共同組成。在保持外形結(jié)構(gòu)尺寸不變時,這種結(jié)構(gòu)設(shè)計可有效提高進出口壓降。同時,通過理論計算、電磁場仿真及試驗測試分析混合流動式磁流變閥的壓降變化規(guī)律,并研究徑向阻尼間隙對混合流動式磁流變閥壓降的影響。

    1 混合流動式磁流變閥結(jié)構(gòu)設(shè)計

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

    磁流變閥性能與其液流通道結(jié)構(gòu)密切相關(guān)。設(shè)計時,首先應使外加磁場產(chǎn)生的磁力線方向與磁流變液流經(jīng)的液流通道相垂直,這樣磁流變液所受剪切應力最大,磁流變效應最好;其次,應盡可能延長有效阻尼間隙長度,從而增大磁流變閥壓降可調(diào)范圍;最后應選擇合適的阻尼間隙厚度以提高進出口壓降。基于此,在課題組前期研究基礎(chǔ)上[14-17],設(shè)計了如圖1所示的混合流動式磁流變閥。該磁流變閥液流通道由單個軸向圓環(huán)流、單個徑向圓盤流和單個中心小孔流串聯(lián)組合而成,構(gòu)成混合流動式液流通道。其中,導磁圓盤與閥體之間圍成的圓環(huán)形阻尼間隙構(gòu)成軸向圓環(huán)流液流通道;導磁圓盤和閥芯之間圍成的圓盤形阻尼間隙構(gòu)成徑向圓盤流液流通道;閥芯中間的通孔形成中心小孔流液流通道。在保證外形結(jié)構(gòu)尺寸的前提下,這種結(jié)構(gòu)設(shè)計有效增加了阻尼長度,從而增大了閥的壓降。工作時,磁流變液首先經(jīng)過左端蓋中心通孔,流經(jīng)導流盤,然后通過軸向圓環(huán)形阻尼間隙,再流經(jīng)徑向圓盤形阻尼間隙,接著到達閥芯中心小孔液流通道,最后從右端蓋中心通孔處流出。

    如圖1所示,定位圓盤和閥芯通過螺紋固定連接,兩者圍成圓環(huán)形繞線槽,勵磁線圈纏繞在圓環(huán)形繞線槽內(nèi)。導磁圓盤、定位圓盤和墊圈通過螺釘固定連接,此時徑向圓盤阻尼間隙厚度可通過墊圈厚度來保證。試驗時,為了分析徑向圓盤阻尼間隙厚度對所設(shè)計的磁流變閥壓降的影響,設(shè)計了0.5、0.8、1.0、1.5 mm 4種規(guī)格的墊圈,從而形成4種徑向圓盤阻尼間隙。另外,為防止磁流變閥出現(xiàn)堵塞問題但又不影響磁流變效應,該磁流變閥軸向圓環(huán)流阻尼間隙厚度設(shè)計為1.0 mm,中心小孔直徑設(shè)計為4.0 mm。整個混合流動式磁流變閥的外形直徑為62 mm,軸向長度為80 mm。

    給勵磁線圈通入一定大小的直流電,在磁場作用下,磁力線經(jīng)過閥芯垂直穿過徑向圓盤阻尼間隙,到達導磁圓盤,再垂直穿過軸向圓環(huán)阻尼間隙,經(jīng)過閥體回到閥芯,形成一個完整的閉合回路。此時軸向圓環(huán)阻尼間隙和徑向圓盤阻尼間隙處的液流方向和磁力線走向相垂直。通過改變電流大小來控制軸向圓環(huán)和徑向圓盤阻尼間隙處的磁感應強度,進而改變剪切應力,可實現(xiàn)對磁流變閥壓降的無級可調(diào)。

    圖1 混合流動磁流變閥結(jié)構(gòu)原理與三維模型圖Fig.1 Schematic diagram of hybrid fluid flow MR valve1.左端蓋 2.導流盤 3.螺釘 4.導磁圓盤 5.墊圈 6. 定位圓盤 7.閥體 8.勵磁線圈 9.閥芯 10.右端蓋

    1.2磁路分析

    圖2所示為混合流動式磁流變閥的磁路簡圖,左右磁路軸對稱。假設(shè)磁力線均勻分布,且不考慮漏磁,則整個閉合回路中有

    ΦMR,r=ΦMR,a=Φsteel=Φ

    (1)

    式中Φ——勵磁線圈產(chǎn)生的磁通量

    ΦMR,a——通過軸向圓環(huán)阻尼間隙處磁通量

    ΦMR,r——通過徑向圓盤阻尼間隙處磁通量

    Φsteel——閥芯、導磁圓盤及閥套等導磁材料的磁通量

    圖2 混合流動磁流變閥磁路簡圖Fig.2 Simplified magnetic circuit of hybrid fluid flow MR valve

    采用基爾霍夫定律對磁路進行分析,可表示為

    (2)

    式中Nc——勵磁線圈匝數(shù)

    I——勵磁線圈中電流

    Hi——磁路中第i部分的磁場強度

    li——磁路中第i部分的有效長度

    通過勵磁線圈的磁通量可表示為

    Φ=∮cBdS=BiSi

    (3)

    式中Bi——磁路中第i部分的磁感應強度

    Si——磁路中第i部分垂直于磁力線的橫截面面積

    圖2中磁路每一部分的有效長度li為

    (4)

    式中L——閥芯長度

    td——定位圓盤厚度

    tc——閥芯右側(cè)厚度

    gr——徑向圓盤阻尼間隙厚度

    La——導磁圓盤厚度

    R——磁流變閥半徑

    th——閥芯厚度

    Rd——閥芯左側(cè)帶螺紋處半徑

    R0——孔口半徑

    ga——軸向圓環(huán)阻尼間隙厚度

    Rc——閥芯左側(cè)不帶螺紋處半徑

    每一部分垂直于磁力線的橫截面積可表示為

    (5)

    式中SMR,r——徑向圓盤阻尼間隙處橫截面積

    SMR,a——軸向圓環(huán)阻尼間隙處橫截面積

    根據(jù)電磁理論,磁感應強度B與磁場強度H兩者之間的關(guān)系為

    Bi=μ0μiHi

    (6)

    式中μ0——真空絕對磁導率,取4π×10-7H/m

    μi——磁流變材料的相對磁導率

    磁路中每一部分的磁阻Ri可表示為

    (7)

    因此,公式(2)可進一步表示為

    (8)

    磁路中各部分的磁感應強度Bi可表示為

    (9)

    式中Bjsat——所有磁性材料最小飽和磁感應強度

    綜合公式(7)~(9),軸向圓環(huán)阻尼間隙處的磁感應強度BMR,a和徑向圓盤阻尼間隙處的磁感應強度BMR,r可分別表示為

    (10)

    (11)

    1.3壓降數(shù)學模型

    忽略磁流變液流經(jīng)左、右端蓋及導流盤中的粘滯壓降,由圖1可知,混合流動式磁流變閥的壓降Δp可表示為

    Δp=Δpa+Δpr+Δpo

    (12)

    式中Δpa——軸向圓環(huán)阻尼間隙處產(chǎn)生的壓降

    Δpr——徑向圓盤阻尼間隙處產(chǎn)生的壓降

    Δpo——閥芯中心小孔處產(chǎn)生的壓降

    軸向圓環(huán)阻尼間隙處產(chǎn)生的壓降Δpa為

    (13)

    式中Δpa,η——軸向圓環(huán)阻尼間隙處產(chǎn)生的粘滯壓降

    Δpa,τ——軸向圓環(huán)阻尼間隙處產(chǎn)生的隨磁場強度變化的屈服壓降

    η——磁流變液零場粘度

    q——系統(tǒng)流量

    τy,a——軸向圓環(huán)阻尼間隙處動態(tài)剪切應力

    c——修正系數(shù),取值范圍為2~3

    徑向圓盤阻尼間隙處產(chǎn)生的壓降Δpr為

    Δpr=Δpr,η+Δpr,τ=

    (14)

    式中Δpr,η——徑向圓盤阻尼間隙處產(chǎn)生的粘滯壓降

    Δpr,τ——徑向圓盤阻尼間隙處產(chǎn)生的隨磁場強度變化的屈服壓降

    τy,r——徑向圓盤阻尼間隙處產(chǎn)生的動態(tài)剪切應力

    Lr——徑向圓盤阻尼間隙處有效阻尼長度

    閥芯中心小孔處產(chǎn)生的壓降Δpo為

    (15)

    把式(13)~(15)代入到式(12)中,可得到混合流動式磁流變閥的壓降Δp為

    (16)

    壓降可調(diào)范圍λ也是一個判斷磁流變閥性能好壞的重要指標,此處定義壓降可調(diào)范圍為Δpτ與Δpη的比值,可表示為

    (17)

    由式(16)可知,影響混合流動式磁流變閥壓降的主要因素包括:軸向圓環(huán)阻尼間隙長度和厚度、徑向圓盤阻尼間隙長度和厚度以及這兩個阻尼間隙處的剪切應力。其中,軸向圓環(huán)阻尼間隙和徑向圓盤阻尼間隙處產(chǎn)生的屈服壓降均隨著阻尼間隙厚度的減小而增大,并隨著阻尼間隙長度的增加而增大;另外,兩個阻尼間隙處的剪切應力的增大也會導致壓降的增大。因此設(shè)計時必須對這些結(jié)構(gòu)參數(shù)和磁場參數(shù)進行綜合考慮。

    2 混合流動式磁流變閥電磁場仿真分析

    混合流動式磁流變閥工作介質(zhì)磁流變液選用重慶材料研究所生產(chǎn)的MRF-J01T型磁流變液。該磁流變液相關(guān)性能參數(shù)如表1所示。由磁流變液的工作模式可知,磁流變液在發(fā)生流變效應時,外界環(huán)境中存在較強的磁場。因此磁流變液的磁場特性對于磁流變閥的設(shè)計至關(guān)重要。應用最小二乘法對MRF-J01T型磁流變液τy-B曲線進行三次多項式擬合[16],可得到軸向圓環(huán)形阻尼間隙及徑向圓盤形阻尼間隙處的動態(tài)剪切屈服應力與磁感應強度的關(guān)系為

    (18)

    (19)

    式中,a3=-984.274 2 kPa/T3,a2=865.390 1 kPa/T2,a1=-48.464 4 kPa/T,a0=0.018 2 kPa。

    表1 MRF-J01T型磁流變液相關(guān)性能參數(shù)Tab.1 Performance index of MR fluid with MRF-J01T

    圖3 混合流動式磁流變閥實體模型Fig.3 Entity model of hybrid fluid flow MR valve1.導磁圓盤 2.液流通道 3. 定位圓盤 4.勵磁線圈 5.閥體 6.閥芯

    為了直觀了解混合流動式磁流變閥的壓降性能,采用有限元仿真軟件ANSYS電磁場模塊建立了圖3所示的有限元實體模型。考慮到混合流動式磁流變閥截面均為規(guī)則的軸對稱圖形,在不影響仿真精度的前提下,取截面的1/2作為仿真對象,以減少計算量;同時選擇PLANE53作為磁場仿真單元。

    圖4所示是加載電流為1.2 A,軸向圓環(huán)形阻尼間隙厚度為1.0 mm,徑向圓盤形阻尼間隙厚度為0.5 mm時的混合流動式磁流變閥磁力線分布圖。由圖可觀察到,該磁流變閥的磁力線絕大部分都垂直穿過徑向圓盤阻尼間隙和軸向圓環(huán)形阻尼間隙,但也存在少部分漏磁,這是因為理想狀況下閥芯繞線圈部分與閥體假定為無窮遠,因此不存在漏磁狀況;但實際情況中兩者之間距離一定。另外,圖4給出了磁流變液流經(jīng)的液流通道長度,其中有效軸向圓環(huán)阻尼長度La為4.2 mm,有效徑向圓盤阻尼長度Lr為11 mm。

    圖4 磁力線分布及液流通道Fig.4 Distribution of magnetic flux line and fluid flow paths

    圖5所示為液流通道處磁感應強度變化曲線。由圖可知,軸向圓環(huán)和徑向圓盤有效阻尼間隙處的磁感應強度均隨著加載電流的增加而增大;另外,相同加載電流下,徑向圓盤有效阻尼間隙處的磁感應強度要大于軸向圓環(huán)有效阻尼間隙處的磁感應強度,但差值不大,說明所設(shè)計的混合流動式磁流變閥磁場利用率高。

    圖5 磁感應強度隨路徑變化曲線Fig.5 Magnetic flux density along the fluid flow paths

    圖6所示為軸向圓環(huán)和徑向圓盤有效阻尼間隙處磁感應強度隨加載電流的變化曲線。由圖可知,同一電流下徑向圓盤有效阻尼間隙處的磁感應強度要略大于軸向圓環(huán)有效阻尼間隙處的磁感應強度,但差值不大,這也進一步表明混合流動式磁流變閥磁場利用率高。

    圖6 磁感應強度隨電流變化曲線Fig.6 Magnetic flux density under applied currents

    圖7 剪切應力隨路徑變化曲線Fig.7 Yield stress along the fluid flow paths

    圖7所示為不同徑向圓盤阻尼間隙厚度時剪切應力隨液流通道路徑的變化曲線。從圖中可觀察到,當加載電流為0.8 A時,徑向圓盤阻尼間隙厚度越小,剪切應力越大,這是由于徑向圓盤阻尼間隙厚度小的磁流變閥獲取的磁感應強度要大于徑向圓盤阻尼間隙厚度大的磁流變閥。另外,當加載電流為2.0 A時,剪切應力基本上保持在70 kPa左右不變,這是因為2.0 A下的磁流變液剪切應力已經(jīng)飽和。

    圖8所示為軸向圓環(huán)阻尼間隙厚度為1.0 mm,徑向圓盤阻尼間隙厚度為0.5 mm時剪切應力及壓降隨電流的變化曲線。從圖中可看出,當加載電流為1.2 A時,徑向圓盤阻尼間隙處的剪切應力與軸向圓環(huán)阻尼間隙處的剪切應力基本上同時達到飽和,此時混合流動式磁流變閥壓降也達到最大值,最大值為3 342 kPa。

    圖8 剪切應力及壓降隨電流變化曲線Fig.8 Yield stress and estimated pressure drop under different applied currents

    圖9所示為軸向圓環(huán)阻尼間隙厚度為1.0 mm,徑向圓盤阻尼間隙厚度為0.5 mm時理論壓降和仿真壓降對比曲線,從圖中可觀察到仿真壓降和理論壓降趨勢基本相同。當加載電流小于0.6 A時,磁滯壓降未完全發(fā)揮作用,壓降模型中的粘滯壓降占整個壓降的比例較大,導致仿真結(jié)果比理論值要大;當加載電流等于0.6 A時,仿真結(jié)果與理論計算結(jié)果一致;但隨著加載電流的增加,仿真結(jié)果比理論值略低,主要是因為加載電流增加到一定值時,磁致壓降發(fā)揮作用并逐漸趨于最大值,由于仿真模型的近似處理,導致仿真計算結(jié)果偏低,但總體趨勢一致。

    圖9 仿真和理論壓降隨電流變化對比曲線Fig.9 Comparison of simulated pressure drop and estimated pressure drop under different applied currents

    圖10所示為軸向圓環(huán)阻尼間隙厚度為1.0 mm,徑向圓盤阻尼間隙厚度分別為0.5、0.8、1.0、1.5 mm時的混合流動式磁流變閥壓降隨電流變化曲線。從圖中可看到,在磁流變液飽和之前,混合磁流變閥壓降隨電流增加而增大;另外,加載電流一定時,徑向圓盤阻尼間隙厚度越小,所對應的壓降越大;徑向圓盤阻尼間隙厚度為0.5 mm,軸向圓環(huán)阻尼間隙厚度為1.0 mm時所對應的壓降最大,最大壓降為3 342 kPa。

    圖10 不同徑向圓盤阻尼間隙厚度時壓降隨電流變化曲線Fig.10 Estimated pressure drop under different applied currents and different radical resistance gaps

    3 混合流動式磁流變閥壓降性能分析

    3.1混合流動式磁流變閥樣機制作

    在理論分析、電磁場仿真及結(jié)構(gòu)設(shè)計基礎(chǔ)上,加工了如圖11所示的混合流動式磁流變閥。該閥的主要結(jié)構(gòu)參數(shù)如表2所示。

    圖11 混合流動式磁流變閥實物圖Fig.11 Prototyping of hybrid fluid flow MR valve

    參數(shù)數(shù)值軸向圓環(huán)阻尼間隙厚度ga/mm1徑向圓盤阻尼間隙厚度gr/mm0.5~1.5閥體厚度th/mm10磁流變閥半徑R/mm31閥芯長度L/mm41繞線槽厚度Wc/mm7孔口半徑R0/mm2閥芯左側(cè)帶螺紋處半徑Rd/mm13閥芯左側(cè)不帶螺紋處半徑Rc/mm14閥芯右側(cè)厚度tc/mm8軸向圓環(huán)阻尼間隙長度La/mm4.2徑向圓盤阻尼間隙長度Lr/mm11定位圓盤厚度td/mm8

    3.2壓降性能試驗

    在自行搭建的磁流變閥性能測試試驗臺上對所設(shè)計的混合流動式磁流變閥壓降性能進行試驗測試分析[16]。為了研究徑向圓盤阻尼間隙對磁流變閥壓降的影響,試驗中加工了如圖12所示的大小相同,厚度分別為0.5、0.8、1.0、1.5 mm的4種不同尺寸的墊圈,用于調(diào)節(jié)徑向圓盤阻尼間隙厚度。

    圖12 4種不同尺寸的墊圈Fig.12 Four sizes of washer used in the experiment

    圖13所示為4種徑向圓盤阻尼間隙厚度下的磁流變閥壓力隨電流變化曲線。由圖可知,磁流變閥的入口壓力和壓降均隨著加載電流的增加而增大,而出口壓力基本保持在700 kPa左右。這是因為壓降是由入口壓力減去出口壓力所得,而出口壓力是通過調(diào)節(jié)試驗系統(tǒng)中的溢流閥旋鈕來實現(xiàn)的。旋鈕固定時,出口壓力也一定。從圖中還可觀察到出口壓力存在少量的波動,這是因為電動機在運轉(zhuǎn)過程中會引起振動,同時齒輪泵由于長期工作,內(nèi)嚙合齒輪會發(fā)生磨損,使得系統(tǒng)流量產(chǎn)生波動。

    圖13 4種徑向圓盤阻尼間隙厚度下壓力隨電流變化曲線Fig.13 Pressure change under different applied currents and radial resistance gaps

    圖14所示為4種徑向圓盤阻尼間隙厚度下混合流動式磁流變閥壓降隨電流變化曲線。由圖可知,壓降隨加載電流的增加而增大,并逐漸趨于最大值并達到穩(wěn)定;徑向圓盤阻尼間隙厚度越小,混合流動式磁流變閥所對應的飽和電流也越??;另外,同一加載電流下,徑向圓盤阻尼間隙厚度越小,系統(tǒng)產(chǎn)生的壓降也越大,說明該磁流變閥可以在體積不變的情況下實現(xiàn)增大壓降的目的。當徑向圓盤阻尼間隙厚度為0.5 mm時,混合流動式磁流變閥所對應的壓降最大,為2 650 kPa。

    圖14 4種徑向圓盤阻尼間隙厚度下壓降隨電流變化曲線Fig.14 Pressure drop change under different applied currents and radial resistance gaps

    通過公式(17)以及圖14,可以得到混合流動式磁流變閥的壓降可調(diào)范圍λ。當徑向圓盤阻尼間隙厚度分別為0.5、0.8、1.0、1.5 mm時,壓降可調(diào)范圍λ分別為7.31、5.21、5.03和4.20。由數(shù)據(jù)可知,徑向圓盤阻尼間隙厚度越小,壓降可調(diào)范圍就越大。這也進一步說明所設(shè)計的混合流動式磁流變閥在保持外形結(jié)構(gòu)尺寸不變的情況下,通過更換墊圈就可實現(xiàn)增大壓降可調(diào)范圍的目的。

    圖15所示為軸向圓環(huán)阻尼間隙厚度為1.0 mm,徑向圓盤阻尼間隙厚度為0.5 mm時試驗壓降和仿真壓降對比曲線。由圖可知,加載電流較小時,試驗壓降與仿真壓降結(jié)果基本相等;但隨著電流的增加,試驗壓降比仿真壓降要低,穩(wěn)定后的差值為692 kPa。產(chǎn)生誤差的主要原因在于磁場作用下,液流通道里面的磁流變液鐵粉顆粒由單鏈狀聚集成柱狀,使得仿真時的飽和剪切應力要大于實際試驗中的飽和剪切應力。另外,試驗過程中的磁場存在一定程度的漏磁,也會導致實際壓降測量結(jié)果偏低,但試驗壓降和仿真壓降總體趨勢一致。

    圖15 仿真與試驗壓降對比曲線Fig.15 Pressure drop comparison of simulated and experimental results

    4 結(jié)論

    (1)提出并設(shè)計了一種結(jié)構(gòu)緊湊的混合流動式磁流變閥,該磁流變閥阻尼間隙處液流通道由單個軸向圓環(huán)流動、單個徑向圓盤流動和單個中心小孔流動串聯(lián)組合而成,整個混合流動式磁流變閥的外形直徑為62 mm,軸向長度為80 mm。在保證外形尺寸的前提下,這種結(jié)構(gòu)設(shè)計增加了有效阻尼長度,閥的壓降也相應得到增加。

    (2)對混合流動式磁流變閥的磁路進行了理論分析,同時推導了其壓降數(shù)學模型。采用有限元法對其電磁場進行了建模仿真,仿真結(jié)果表明,當軸向圓環(huán)阻尼間隙厚度為1.0 mm,徑向圓盤阻尼間隙厚度為0.5 mm,加載電流為1.2 A時壓降最大,為3 342 kPa。

    (3)對混合流動式磁流變閥壓降性能進行了試驗測試,試驗結(jié)果表明,當軸向圓環(huán)阻尼間隙厚度為1.0 mm、徑向圓盤阻尼間隙厚度為0.5 mm時,加載電流為1.2 A時的壓降最大,為2 650 kPa。試驗壓降和仿真壓降總體趨勢一致。

    (4)仿真和試驗結(jié)果表明徑向圓盤阻尼間隙厚度越小,所對應的壓降和壓降可調(diào)范圍就越大,這進一步說明所設(shè)計的混合流動式磁流變閥在保持體積一定的情況下,能夠?qū)崿F(xiàn)增大壓降和壓降可調(diào)范圍的目的。

    1GONCALVES F D, KOO J H, AHMADIAN M. A review of the state of the art in magnetorheological fluid technologies—Part I: MR fluid and MR fluid models[J]. The Shock and Vibration Digest, 2006, 38(3): 203-219.

    2ZHU X, JING X, CHENG L. Magnetorheological fluid dampers: a review on structure design and analysis[J]. Journal of Intelligent Material Systems and Structures, 2012, 23(8): 839-873.

    3SUN S, DENG H, LI W, et al. Improving the critical speeds of high-speed trains using magnetorheological technology[J]. Smart Materials and Structures, 2013, 22(11): 115012.

    4SAHIN H, WANG X, GORDANINEJAD F. Magneto-rheological fluid flow through complex valve geometries[J]. International Journal of Vehicle Design, 2013, 63(2-3): 241-255.

    5NGUYEN Q H, CHOI S B, LEE Y S, et al. Optimal design of high damping force engine mount featuring MR valve structure with both annular and radial flow paths[J]. Smart Materials and Structures, 2013, 22(11): 115024.

    6SHAH K, CHOI S B. A new magnetorheological mount featured by changeable damping gaps using a moved-plate valve structure[J]. Smart Materials and Structures, 2014, 23(12): 125022.

    7LI W H, DU H, GUO N Q. Finite element analysis and simulation evaluation of a magnetorheological valve[J]. The International Journal of Advanced Manufacturing Technology, 2003, 21(6): 438-445.

    8ABD FATAH A Y, MAZLAN S A, KOGA T, et al. A review of design and modeling of magnetorheological valve[J]. International Journal of Modern Physics B, 2015, 29(4): 1530004.

    9GORODKIN S, LUKIANOVICH A. Magnetorheological throttle valve in passive damping system[J]. Journal of Intelligent Material Systems and Structures, 1998, 9(8): 637-641.

    10SALLOOM M Y, SAMAD Z. Magneto-rheological directional control valve[J]. The International Journal of Advanced Manufacturing Technology, 2012, 58(1-4): 279-292.

    11SALLOOM M Y, SAMAD Z. Design and modeling magneto rheological directional control valve[J]. Journal of Intelligent Material Systems and Structures, 2012, 23(2): 155-167.

    12ICHWAN B, MAZLAN S A, IMADUDDIN F, et al. Development of a modular MR valve using meandering flow path structure[J]. Smart Materials and Structures, 2016, 25(3): 037001.

    13趙燦, 劉丹丹, 湯春瑞. 一種非直線液流通道的磁流變閥: 中國,200810065763.9[P]. 2010-06-02.

    ZHAO Can, LIU Dandan, TANG Chunrui. A magnetorheological valve with a non-linear flow path: ZL 200810065763.9[P]. 2010-06-02.(in Chinese)

    14HU G, LONG M, HUANG M, et al. Design, analysis, prototyping, and experimental evaluation of an efficient double coil magnetorheological valve[J]. Advances in Mechanical Engineering, 2014, 6: 403410.

    15HU G, LONG M, YU L, et al. Design and performance evaluation of a novel magnetorheological valve with a tunable resistance gap[J]. Smart Materials and Structures, 2014, 23(12): 127001.

    16胡國良, 李海燕, 張海云. 圓環(huán)流磁流變閥壓降性能分析與試驗[J]. 農(nóng)業(yè)機械學報, 2016, 47(3): 381-388.

    HU Guoliang, LI Haiyan, ZHANG Haiyun. Performance analysis and experimental tests of pressure drop of annular type magnetorheological valve[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3): 381-388. (in Chinese)

    17胡國良, 李海燕, 李衛(wèi)華. 徑向流和圓環(huán)流磁流變閥壓降性能對比與試驗[J]. 農(nóng)業(yè)機械學報, 2016, 47(4): 364-371, 405.

    HU Guoliang, LI Haiyan, LI Weihua. Comparison and experiment of pressure drop of radial and annular type magnetorheological valves[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(4): 364-371, 405. (in Chinese)

    18NGUYEN Q H, HAN Y M, CHOI S B, et al. Geometry optimization of MR valves constrained in a specific volume using the finite element method[J]. Smart Materials and Structures, 2007, 16(6): 2242-2252.

    19NGUYEN Q H, CHOI S B, WERELEY N M. Optimal design of magnetorheological valves via a finite element method considering control energy and a time constant[J]. Smart Materials and Structures, 2008, 17(2): 025024.

    20NGUYEN Q H, CHOI S B, LEE Y S, et al. An analytical method for optimal design of MR valve structures[J]. Smart Materials and Structures, 2009, 18(9): 095032.

    21AI H X, WANG D H, LIAO W H. Design and modeling of a magnetorheological valve with both annular and radial flow paths[J]. Journal of Intelligent Material Systems and Structures, 2006, 17(4): 327-334.

    22WANG D H, AI H X, LIAO W H. A magnetorheological valve with both annular and radial fluid flow resistance gaps[J]. Smart Materials and Structures, 2009, 18(11): 115001.

    23IMADUDDIN F, MAZLAN S A, RAHMAN M A A, et al. A high performance magnetorheological valve with a meandering flow path[J]. Smart Materials and Structures, 2014, 23(6): 065017.

    24IMADUDDIN F, MAZLAN S A, ZAMZURI H, et al. Design and performance analysis of a compact magnetorheological valve with multiple annular and radial gaps[J]. Journal of Intelligent Material Systems and Structures, 2015, 26(9): 1038-1049.

    Design and Pressure Drop Analysis of Hybrid Fluid Flow Magnetorheological Valve

    Hu GuoliangZhong FangLiao MingkeGong Meihui

    (KeyLaboratoryofConveyanceandEquipment,MinistryofEducation,EastChinaJiaotongUniversity,Nanchang330013,China)

    Magnetorheological (MR) valve is a kind of smart control mechanism that using magnetorheologcial fluid as the working fluid. A compact hybrid fluid flow MR valve was developed and prototyped, the fluid flow paths of the proposed MR valve were consisted of a single annular flow channel, a single radial flow channel and an orifice flow channel through structural design. In order to investigate the changes of the magnetic flux density and yield stress along the fluid flow paths under the four different radial resistance gaps, the magnetic circuit design and finite element analysis were carried out. The simulation results showed that this design methodology can provide a larger pressure drop of 3 342 kPa at 0.5 mm radial resistance gap and 1.0 mm annular resistance gap, which constrained in a volume space of MR valve with an outer diameter of 62 mm and overall length of 80 mm. Meanwhile, the experimental pressure drop under the four different radial resistance gaps was also carried out on the test rig. The results showed that the pressure drop was increased as the applied current increased, and the pressure drop also increased as the radial resistance gap decreased from 1.5 mm to 0.5 mm, the maximum pressure drop was approximately 2 650 kPa at the applied current of 1.2 A when the radial resistance gap was 0.5 mm, which was well accorded with the simulation results. This proposed MR valve significantly improved its efficiency through compact design and changeable radial resistance gap, which is beneficial to control the MR damper as a bypass control valve for its good pressure regulating capability.

    MR valve; hybrid fluid flow; pressure drop; finite element analysis

    10.6041/j.issn.1000-1298.2016.09.052

    2016-06-09

    2016-07-06

    國家自然科學基金項目(51165005、51475165、11462004)和江西省自然科學基金項目(20151BAB206035)

    胡國良(1973—),男,教授,博士生導師,主要從事磁流變智能器件及結(jié)構(gòu)、流體傳動及控制研究,E-mail: glhu@ecjtu.edu.cn

    TH137.5

    A

    1000-1298(2016)09-0389-09

    猜你喜歡
    圓盤圓環(huán)阻尼
    加權(quán)全能量最小的圓環(huán)形變
    豬圓環(huán)病毒病的發(fā)生、診斷和防治
    一例鴨圓環(huán)病毒病的診斷
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項的擴散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    圓盤鋸刀頭的一種改進工藝
    石材(2020年6期)2020-08-24 08:27:00
    圓環(huán)上的覆蓋曲面不等式及其應用
    單位圓盤上全純映照模的精細Schwarz引理
    奇怪的大圓盤
    国产亚洲精品久久久com| 看片在线看免费视频| 人妻制服诱惑在线中文字幕| 亚洲一区二区三区色噜噜| 日韩欧美精品v在线| 免费高清视频大片| 男插女下体视频免费在线播放| 欧美3d第一页| 成人精品一区二区免费| 国产高清视频在线播放一区| 国产男人的电影天堂91| 好男人在线观看高清免费视频| 欧美高清性xxxxhd video| 午夜免费男女啪啪视频观看 | 亚洲最大成人av| 在现免费观看毛片| 久久精品国产鲁丝片午夜精品| 在线播放无遮挡| 欧美三级亚洲精品| 免费电影在线观看免费观看| 亚洲精品国产av成人精品 | 久久中文看片网| 99国产精品一区二区蜜桃av| 午夜福利成人在线免费观看| 亚洲国产精品合色在线| 最新中文字幕久久久久| 亚洲欧美日韩卡通动漫| 成人鲁丝片一二三区免费| 精品久久久久久久人妻蜜臀av| 国产老妇女一区| 在线观看一区二区三区| 欧美性感艳星| 欧美中文日本在线观看视频| 天天一区二区日本电影三级| 一级a爱片免费观看的视频| 舔av片在线| 51国产日韩欧美| 一级毛片aaaaaa免费看小| 夜夜看夜夜爽夜夜摸| 久久精品夜夜夜夜夜久久蜜豆| 嫩草影院新地址| 啦啦啦啦在线视频资源| 婷婷亚洲欧美| 国产蜜桃级精品一区二区三区| 丝袜美腿在线中文| 国内久久婷婷六月综合欲色啪| 天天一区二区日本电影三级| 亚洲美女视频黄频| 性色avwww在线观看| 亚洲精品乱码久久久v下载方式| 日韩国内少妇激情av| 无遮挡黄片免费观看| 狂野欧美激情性xxxx在线观看| 成人性生交大片免费视频hd| 少妇的逼水好多| 午夜福利高清视频| 久久精品国产清高在天天线| 女生性感内裤真人,穿戴方法视频| 欧美日韩精品成人综合77777| 美女大奶头视频| 国产真实伦视频高清在线观看| 国产午夜精品论理片| 午夜免费激情av| 99热6这里只有精品| 特大巨黑吊av在线直播| 国产 一区 欧美 日韩| 久久久久久久久久成人| 变态另类成人亚洲欧美熟女| 高清午夜精品一区二区三区 | 97超碰精品成人国产| 免费无遮挡裸体视频| av在线蜜桃| 国产 一区精品| 欧美精品国产亚洲| 亚洲性夜色夜夜综合| 一边摸一边抽搐一进一小说| 男人和女人高潮做爰伦理| 免费看av在线观看网站| 极品教师在线视频| 禁无遮挡网站| 超碰av人人做人人爽久久| or卡值多少钱| 欧美性感艳星| 国产视频内射| 国产成人aa在线观看| 久久99热6这里只有精品| 亚洲国产精品国产精品| 美女被艹到高潮喷水动态| 一个人看视频在线观看www免费| 亚洲精品一卡2卡三卡4卡5卡| 99久久九九国产精品国产免费| 国产aⅴ精品一区二区三区波| 亚洲天堂国产精品一区在线| 在线播放无遮挡| 高清午夜精品一区二区三区 | av国产免费在线观看| 国产av不卡久久| 两个人视频免费观看高清| 亚洲国产色片| 成熟少妇高潮喷水视频| 亚洲欧美日韩高清在线视频| 日本熟妇午夜| 两个人的视频大全免费| 国内揄拍国产精品人妻在线| 91久久精品国产一区二区三区| 一个人观看的视频www高清免费观看| 亚洲中文日韩欧美视频| 18禁黄网站禁片免费观看直播| 69av精品久久久久久| 校园春色视频在线观看| 国产成人精品久久久久久| 日本精品一区二区三区蜜桃| 麻豆成人午夜福利视频| 久久精品国产亚洲av涩爱 | 一夜夜www| 午夜精品国产一区二区电影 | 日韩精品青青久久久久久| 国产精品爽爽va在线观看网站| aaaaa片日本免费| 国产成人一区二区在线| 精品久久久久久久久亚洲| av视频在线观看入口| 日本三级黄在线观看| 欧美区成人在线视频| 五月玫瑰六月丁香| 免费看av在线观看网站| 夜夜看夜夜爽夜夜摸| 国产欧美日韩精品亚洲av| 乱系列少妇在线播放| 亚洲人成网站在线观看播放| 国内揄拍国产精品人妻在线| 日本欧美国产在线视频| 特级一级黄色大片| 色吧在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲国产日韩欧美精品在线观看| 免费电影在线观看免费观看| 欧美丝袜亚洲另类| 国产高清不卡午夜福利| 亚洲性夜色夜夜综合| 一卡2卡三卡四卡精品乱码亚洲| 在线a可以看的网站| 精品午夜福利在线看| 亚洲人与动物交配视频| 免费av不卡在线播放| 深夜精品福利| 精品午夜福利视频在线观看一区| 三级国产精品欧美在线观看| 一级毛片久久久久久久久女| 欧美一级a爱片免费观看看| 日本免费一区二区三区高清不卡| 69av精品久久久久久| 国产综合懂色| 亚洲真实伦在线观看| 91久久精品国产一区二区成人| 国产综合懂色| 欧美日韩在线观看h| 嫩草影院入口| 赤兔流量卡办理| 亚洲成人精品中文字幕电影| 男人的好看免费观看在线视频| 日韩大尺度精品在线看网址| 国产精品久久久久久亚洲av鲁大| 中国美女看黄片| 中文在线观看免费www的网站| 日韩一本色道免费dvd| 偷拍熟女少妇极品色| 色av中文字幕| 身体一侧抽搐| 狠狠狠狠99中文字幕| 国产av一区在线观看免费| 精品少妇黑人巨大在线播放 | 亚洲精品一区av在线观看| 日产精品乱码卡一卡2卡三| 亚州av有码| 校园人妻丝袜中文字幕| 久久久国产成人免费| 亚洲国产精品久久男人天堂| 成年版毛片免费区| 少妇丰满av| 日本免费a在线| 久久99热6这里只有精品| av专区在线播放| 亚洲成人久久性| 黑人高潮一二区| 免费观看的影片在线观看| 色av中文字幕| 日韩,欧美,国产一区二区三区 | 久久久久久久午夜电影| 女人十人毛片免费观看3o分钟| 国产美女午夜福利| 国产成人a区在线观看| 欧美不卡视频在线免费观看| 日日啪夜夜撸| 亚洲av成人av| 午夜免费激情av| 久久精品国产亚洲av天美| 欧美人与善性xxx| 亚洲精品一区av在线观看| 欧美一区二区国产精品久久精品| 精品欧美国产一区二区三| 成年av动漫网址| 全区人妻精品视频| 久久综合国产亚洲精品| 久久久久久久亚洲中文字幕| 欧美日韩在线观看h| 国产精品久久久久久av不卡| 搡老熟女国产l中国老女人| 少妇被粗大猛烈的视频| 亚洲精品国产av成人精品 | 91久久精品国产一区二区成人| 国产色爽女视频免费观看| 波多野结衣高清无吗| 两个人视频免费观看高清| 久久久精品大字幕| 美女大奶头视频| 高清毛片免费观看视频网站| 亚洲精华国产精华液的使用体验 | 嫩草影院新地址| 久久人妻av系列| 日韩人妻高清精品专区| 国产精品久久久久久av不卡| 中文资源天堂在线| 亚洲国产欧洲综合997久久,| 一卡2卡三卡四卡精品乱码亚洲| 高清午夜精品一区二区三区 | 久久99热6这里只有精品| 婷婷精品国产亚洲av| 床上黄色一级片| 久久久久性生活片| 高清午夜精品一区二区三区 | 干丝袜人妻中文字幕| 国产精品一区二区三区四区免费观看 | 俺也久久电影网| 国产探花在线观看一区二区| 色噜噜av男人的天堂激情| 国产乱人偷精品视频| 国产精品一区二区性色av| 欧美一区二区亚洲| 黄色欧美视频在线观看| 最新中文字幕久久久久| 亚洲一区高清亚洲精品| av.在线天堂| 免费观看人在逋| 亚洲人成网站高清观看| 一区福利在线观看| 乱系列少妇在线播放| 亚洲av二区三区四区| 露出奶头的视频| 日本在线视频免费播放| 波野结衣二区三区在线| 一级a爱片免费观看的视频| 成人永久免费在线观看视频| 亚洲欧美精品自产自拍| av中文乱码字幕在线| 国产高清激情床上av| 免费人成视频x8x8入口观看| 一卡2卡三卡四卡精品乱码亚洲| 人妻久久中文字幕网| 免费一级毛片在线播放高清视频| 国内少妇人妻偷人精品xxx网站| 天天躁日日操中文字幕| 成人高潮视频无遮挡免费网站| 精品午夜福利在线看| 人妻久久中文字幕网| 嫩草影视91久久| 国模一区二区三区四区视频| 少妇熟女aⅴ在线视频| 精品久久国产蜜桃| 久久久精品大字幕| 我的女老师完整版在线观看| 亚洲欧美精品综合久久99| 亚洲成人久久性| 两个人视频免费观看高清| 99精品在免费线老司机午夜| 亚州av有码| 校园春色视频在线观看| 成人三级黄色视频| 一级黄片播放器| 全区人妻精品视频| 午夜福利在线观看吧| 欧美高清成人免费视频www| 色综合亚洲欧美另类图片| 成年免费大片在线观看| 国产午夜福利久久久久久| 两个人视频免费观看高清| 日日摸夜夜添夜夜添小说| 精品国产三级普通话版| 99热6这里只有精品| 国产精品久久久久久久电影| 最近中文字幕高清免费大全6| 国产在线精品亚洲第一网站| 日韩制服骚丝袜av| 女的被弄到高潮叫床怎么办| 自拍偷自拍亚洲精品老妇| 男人的好看免费观看在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久国产蜜桃| 精品久久久噜噜| 国产蜜桃级精品一区二区三区| 日本黄色片子视频| 中文字幕av成人在线电影| 亚洲五月天丁香| 97超碰精品成人国产| АⅤ资源中文在线天堂| 丝袜喷水一区| 久久亚洲精品不卡| 一卡2卡三卡四卡精品乱码亚洲| 国产精品福利在线免费观看| 色视频www国产| 久久久久久九九精品二区国产| 国产亚洲精品久久久com| 免费av不卡在线播放| 亚洲欧美日韩卡通动漫| 国产毛片a区久久久久| 18禁在线播放成人免费| 色5月婷婷丁香| 淫秽高清视频在线观看| 亚洲精华国产精华液的使用体验 | av卡一久久| 日本熟妇午夜| 亚洲av熟女| 熟妇人妻久久中文字幕3abv| 99在线人妻在线中文字幕| 桃色一区二区三区在线观看| 精品人妻偷拍中文字幕| 一本一本综合久久| 亚洲欧美日韩东京热| 床上黄色一级片| 亚洲在线观看片| 香蕉av资源在线| 蜜桃久久精品国产亚洲av| 18禁黄网站禁片免费观看直播| 亚洲一区高清亚洲精品| 日本色播在线视频| 极品教师在线视频| 一a级毛片在线观看| 国产v大片淫在线免费观看| 国产午夜精品久久久久久一区二区三区 | 淫妇啪啪啪对白视频| 99久久中文字幕三级久久日本| 国产精品久久电影中文字幕| 国内揄拍国产精品人妻在线| 午夜精品国产一区二区电影 | 国产亚洲91精品色在线| 中文字幕久久专区| 床上黄色一级片| 精品国内亚洲2022精品成人| 国产欧美日韩精品亚洲av| 黄色欧美视频在线观看| 亚洲av成人av| 天堂av国产一区二区熟女人妻| 麻豆国产97在线/欧美| 啦啦啦观看免费观看视频高清| 久久这里只有精品中国| 中文字幕人妻熟人妻熟丝袜美| 搡女人真爽免费视频火全软件 | 99热只有精品国产| 高清毛片免费看| 一级a爱片免费观看的视频| 最新在线观看一区二区三区| 国产精品无大码| 干丝袜人妻中文字幕| 免费看a级黄色片| 校园人妻丝袜中文字幕| av在线亚洲专区| 搡女人真爽免费视频火全软件 | 中文亚洲av片在线观看爽| 久久午夜福利片| 大又大粗又爽又黄少妇毛片口| 亚州av有码| 亚洲人成网站高清观看| 亚洲av美国av| 国产单亲对白刺激| 黄片wwwwww| 精品久久久久久成人av| 午夜福利18| 欧美丝袜亚洲另类| 久久人人精品亚洲av| 成人鲁丝片一二三区免费| 国产真实伦视频高清在线观看| 久久精品国产99精品国产亚洲性色| 日韩欧美精品v在线| 一级av片app| 成年女人看的毛片在线观看| 女人被狂操c到高潮| 成人亚洲欧美一区二区av| 性色avwww在线观看| 你懂的网址亚洲精品在线观看 | 久久久久久久久大av| 亚洲四区av| 国产一级毛片七仙女欲春2| 国产亚洲av嫩草精品影院| 日韩精品中文字幕看吧| 久久久精品94久久精品| 精品人妻一区二区三区麻豆 | 久久久久国内视频| 搞女人的毛片| 69人妻影院| 不卡视频在线观看欧美| 中文字幕熟女人妻在线| 亚洲欧美日韩东京热| 亚洲aⅴ乱码一区二区在线播放| 人人妻人人看人人澡| 淫秽高清视频在线观看| 美女高潮的动态| 精品人妻偷拍中文字幕| 色视频www国产| 12—13女人毛片做爰片一| 俺也久久电影网| 日韩国内少妇激情av| 久久精品国产亚洲av香蕉五月| 一级黄片播放器| 网址你懂的国产日韩在线| 欧美日韩乱码在线| 99热这里只有是精品50| 1000部很黄的大片| 床上黄色一级片| 女人被狂操c到高潮| 少妇高潮的动态图| 91久久精品国产一区二区三区| 国产 一区 欧美 日韩| 免费高清视频大片| 国产久久久一区二区三区| 久久精品国产自在天天线| 99久国产av精品| av.在线天堂| 我要搜黄色片| 国产69精品久久久久777片| 国产精品1区2区在线观看.| 女人十人毛片免费观看3o分钟| 在线播放国产精品三级| 日本熟妇午夜| 一级毛片电影观看 | 国产一区二区亚洲精品在线观看| 看黄色毛片网站| 色吧在线观看| 色av中文字幕| 欧美日韩精品成人综合77777| 18禁裸乳无遮挡免费网站照片| 精品一区二区三区人妻视频| 性色avwww在线观看| 亚洲成人久久性| 午夜老司机福利剧场| 美女被艹到高潮喷水动态| 成人漫画全彩无遮挡| 亚洲av中文字字幕乱码综合| 一个人免费在线观看电影| 精品国内亚洲2022精品成人| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 欧美精品国产亚洲| 别揉我奶头 嗯啊视频| 国产精品不卡视频一区二区| 久久精品国产亚洲av香蕉五月| 久久午夜福利片| 亚洲中文日韩欧美视频| 99久久久亚洲精品蜜臀av| 女生性感内裤真人,穿戴方法视频| 中文字幕人妻熟人妻熟丝袜美| 一区福利在线观看| 精品久久久噜噜| 久久久久久伊人网av| 成人一区二区视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 韩国av在线不卡| 亚洲欧美日韩无卡精品| 99热全是精品| av在线观看视频网站免费| 日韩亚洲欧美综合| 欧美日本亚洲视频在线播放| 可以在线观看的亚洲视频| 午夜福利在线在线| 国产黄色视频一区二区在线观看 | 久久精品国产鲁丝片午夜精品| 搡老岳熟女国产| 一本精品99久久精品77| 男女下面进入的视频免费午夜| av在线老鸭窝| 亚洲精品成人久久久久久| 日本欧美国产在线视频| 99精品在免费线老司机午夜| 久久精品影院6| 在线播放国产精品三级| 国产精品一及| 国产精品一二三区在线看| 一区福利在线观看| 成人鲁丝片一二三区免费| 色综合亚洲欧美另类图片| 亚洲国产精品久久男人天堂| 精品久久国产蜜桃| 天堂av国产一区二区熟女人妻| 亚洲av第一区精品v没综合| 国产精品电影一区二区三区| 一个人免费在线观看电影| 中文资源天堂在线| 久久午夜福利片| 波野结衣二区三区在线| av在线播放精品| 日本-黄色视频高清免费观看| 久久韩国三级中文字幕| 国产成人精品久久久久久| 午夜免费激情av| 亚洲精品在线观看二区| 国产黄色视频一区二区在线观看 | 黄色一级大片看看| 狠狠狠狠99中文字幕| 亚洲欧美精品自产自拍| 久久久色成人| 日本与韩国留学比较| 又爽又黄a免费视频| 亚洲无线在线观看| 成年版毛片免费区| 国产成人aa在线观看| 男人舔女人下体高潮全视频| 一级毛片aaaaaa免费看小| 亚洲七黄色美女视频| 寂寞人妻少妇视频99o| 国产单亲对白刺激| 国产一区亚洲一区在线观看| av在线亚洲专区| 俺也久久电影网| 国产av在哪里看| 国产视频内射| 美女黄网站色视频| 天堂网av新在线| 免费大片18禁| 少妇熟女aⅴ在线视频| 日韩亚洲欧美综合| 三级国产精品欧美在线观看| 99久国产av精品| 长腿黑丝高跟| 麻豆一二三区av精品| 亚洲aⅴ乱码一区二区在线播放| 国产免费男女视频| 国产成人a∨麻豆精品| 日韩一本色道免费dvd| 国模一区二区三区四区视频| 日日啪夜夜撸| 精品人妻一区二区三区麻豆 | 国国产精品蜜臀av免费| 日韩,欧美,国产一区二区三区 | 少妇的逼好多水| 成人综合一区亚洲| а√天堂www在线а√下载| 精品一区二区三区视频在线观看免费| 欧美色视频一区免费| 日本色播在线视频| 成人鲁丝片一二三区免费| 少妇丰满av| 国产激情偷乱视频一区二区| 日本爱情动作片www.在线观看 | 亚洲av中文av极速乱| 在线播放国产精品三级| 国产精品国产三级国产av玫瑰| 成人美女网站在线观看视频| 欧美一区二区国产精品久久精品| av天堂中文字幕网| 99热全是精品| 亚洲精品色激情综合| 国产精品99久久久久久久久| 成人午夜高清在线视频| 不卡视频在线观看欧美| 亚洲高清免费不卡视频| 国产探花极品一区二区| 一区二区三区四区激情视频 | 免费电影在线观看免费观看| 精品99又大又爽又粗少妇毛片| 国产伦精品一区二区三区视频9| 麻豆精品久久久久久蜜桃| 少妇人妻一区二区三区视频| 大又大粗又爽又黄少妇毛片口| 天堂动漫精品| 成年版毛片免费区| 久久久色成人| 日本五十路高清| 成年女人毛片免费观看观看9| 亚洲欧美日韩无卡精品| 亚洲人成网站在线播| 97在线视频观看| 国产精品亚洲美女久久久| 麻豆国产av国片精品| 欧美最黄视频在线播放免费| 丝袜美腿在线中文| 99国产极品粉嫩在线观看| 国产精品久久久久久久久免| 亚洲内射少妇av| 青春草视频在线免费观看| 日韩强制内射视频| 精品一区二区三区视频在线观看免费| 91精品国产九色| 久久久成人免费电影| 国产成人一区二区在线| 麻豆成人午夜福利视频| 菩萨蛮人人尽说江南好唐韦庄 | 男女视频在线观看网站免费| 波多野结衣高清无吗| 久久久久久久亚洲中文字幕| 亚洲内射少妇av| 熟妇人妻久久中文字幕3abv| 女人十人毛片免费观看3o分钟| 亚洲精华国产精华液的使用体验 | 欧美xxxx性猛交bbbb| 国产午夜福利久久久久久| 一个人观看的视频www高清免费观看| 看十八女毛片水多多多| av视频在线观看入口| 91在线观看av| 精品久久国产蜜桃| 联通29元200g的流量卡| 国产精品爽爽va在线观看网站| 日韩欧美在线乱码| 亚洲七黄色美女视频| 久久久久久久久久成人| 一区二区三区免费毛片|