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

    單線圈徑向流磁流變閥優(yōu)化設(shè)計(jì)與性能分析

    2019-08-13 01:42:52胡國良李林森胡安琪喻理梵
    關(guān)鍵詞:磁路磁感應(yīng)阻尼

    胡國良 李林森 胡安琪 喻理梵

    (華東交通大學(xué)載運(yùn)工具與裝備教育部重點(diǎn)實(shí)驗(yàn)室, 南昌 330013)

    0 引言

    液壓系統(tǒng)在處理生產(chǎn)活動(dòng)中大慣量與重負(fù)載扭矩方面的良好表現(xiàn)使其在工業(yè)應(yīng)用領(lǐng)域占有重要地位。作為現(xiàn)代液壓系統(tǒng)中的核心部件,電液伺服閥可以實(shí)現(xiàn)精確、快速的控制,在精確定位與速度控制領(lǐng)域獲得了廣泛的應(yīng)用。但電液伺服閥系統(tǒng)本身結(jié)構(gòu)較為復(fù)雜,且存在響應(yīng)慢、成本高、體積大等缺點(diǎn),因此替代液壓系統(tǒng)傳統(tǒng)驅(qū)動(dòng)機(jī)制的研究備受關(guān)注。新興智能流體材料——磁流變液可通過對其流體域施加磁場作用從而改變其流變特性,以磁流變液為工作介質(zhì)的閥、減振器以及制動(dòng)器等器件可以在不同程度上改善相應(yīng)領(lǐng)域傳動(dòng)性能。磁流變閥通過改變工作電流、調(diào)節(jié)內(nèi)部阻尼通道處的磁感應(yīng)強(qiáng)度,即可調(diào)節(jié)磁流變液流變程度,從而控制磁流變閥內(nèi)部壓降。磁流變閥無復(fù)雜的活動(dòng)部件,結(jié)構(gòu)更加緊湊,響應(yīng)速度更快,壓降調(diào)節(jié)也更加平穩(wěn)[1-9]。

    磁流變閥工作模式一般為流動(dòng)模式,通過增加有效阻尼間隙內(nèi)的磁感應(yīng)強(qiáng)度,可增大磁流變液的剪切屈服應(yīng)力,進(jìn)而有效提升磁流變閥的壓降性能。根據(jù)磁路歐姆定律可知,通過增加磁動(dòng)勢或降低磁阻可以增加磁路中的磁場強(qiáng)度,因此主要通過增加勵(lì)磁線圈數(shù)量和減小有效阻尼間隙厚度來提升磁流變閥的壓降性能。在普通單線圈磁流變閥基礎(chǔ)上增加勵(lì)磁線圈數(shù)量,可在一定程度上增大壓降可控范圍,降低磁路響應(yīng)時(shí)間,但線圈個(gè)數(shù)增加會導(dǎo)致閥體積與能耗增加[10-11]。HU等[12]利用調(diào)節(jié)阻尼間隙厚度能改變磁流變閥壓降性能的特性,設(shè)計(jì)了一種阻尼間隙可調(diào)式磁流變閥,通過閥芯在圓錐環(huán)形通道內(nèi)的軸向移動(dòng),可實(shí)現(xiàn)阻尼間隙在1~2 mm間的無級切換,進(jìn)一步提高了磁流變閥壓降可調(diào)范圍。但活動(dòng)部件的增加導(dǎo)致結(jié)構(gòu)較復(fù)雜,且較小的阻尼通道易堵塞液流通道,從而影響磁流變閥工作穩(wěn)定性。

    增大磁流變閥有效阻尼間隙是另一種提高壓降性能的方法,單純增加軸向或徑向液流通道長度會導(dǎo)致閥外形尺寸的增加。近年來,國內(nèi)外學(xué)者主要采用在磁流變閥內(nèi)設(shè)置蜿蜒式路徑的方式,間接增加其有效液流通道長度[13-17]。

    改變磁流變閥內(nèi)部結(jié)構(gòu)提高壓降性能的方法會在不同程度上增加磁流變閥結(jié)構(gòu)的復(fù)雜性以及工作的不穩(wěn)定性。基于此,在不改變傳統(tǒng)閥結(jié)構(gòu)基礎(chǔ)上,部分學(xué)者提出多種對尺寸參數(shù)進(jìn)行優(yōu)化方法,以提高其壓降性能[18-21]。

    上述關(guān)于磁流變閥性能優(yōu)化的研究中,只對特定應(yīng)用場合以及特定種類的磁流變閥設(shè)計(jì)具有一定的參考價(jià)值,且磁流變閥體積一定時(shí),徑向流結(jié)構(gòu)較軸向流結(jié)構(gòu)能達(dá)到的可調(diào)壓降范圍更大。

    本文利用ANSYS有限元仿真軟件中的優(yōu)化工具,在分析尺寸參數(shù)對關(guān)鍵性能影響的基礎(chǔ)上,對單線圈徑向流磁流變閥可控壓降性能進(jìn)行優(yōu)化。不改變普通單線圈徑向流磁流變閥內(nèi)部基本結(jié)構(gòu)及外觀尺寸,通過增加隔磁材料、引導(dǎo)磁力線走向改善磁流變閥的磁路設(shè)計(jì),以期降低磁路磁漏,增加間隙處的磁感應(yīng)強(qiáng)度。利用有限元分析軟件的Design Opt模塊,結(jié)合零階與一階優(yōu)化工具,通過對目標(biāo)函數(shù)設(shè)立適當(dāng)?shù)臋?quán)重系數(shù),對壓降與壓降可調(diào)系數(shù)進(jìn)行多目標(biāo)優(yōu)化,以實(shí)現(xiàn)對磁流變閥最優(yōu)結(jié)構(gòu)尺寸的快速準(zhǔn)確求解。對優(yōu)化前后的單線圈徑向流磁流變閥進(jìn)行仿真對比,以期得到不同電流下有效阻尼通道的平均磁感應(yīng)強(qiáng)度與壓降變化規(guī)律,并在磁流變閥性能測試平臺上對兩種磁流變閥的壓降及響應(yīng)性能進(jìn)行測試。

    1 工作原理及結(jié)構(gòu)

    1.1 工作原理

    圖1 單線圈徑向流磁流變閥結(jié)構(gòu)示意圖Fig.1 Structure diagram of single coil radial flow MR valve1.左端蓋 2.側(cè)翼套筒 3.閥體 4.隔磁套筒 5.勵(lì)磁線圈 6.閥芯 7.定位盤 8.右端蓋

    磁流變液的工作模式?jīng)Q定了垂直于磁流變液流動(dòng)方向的磁場分量最有利于其流變效應(yīng)的充分發(fā)揮,因此磁力線更均勻集中的垂直穿過有效液流通道能最大程度改善磁流變閥可控壓降性能?;诖?,在傳統(tǒng)磁流變閥的研究基礎(chǔ)上,設(shè)計(jì)了單線圈徑向流磁流變閥(圖1)。該磁流變閥液流通道由兩段中心小孔流、兩段徑向圓盤流與一段軸向圓環(huán)流共同組成,其中左右兩端定位盤中間通孔構(gòu)成中心小孔流通道,閥芯與左右兩端定位盤圍成徑向圓盤流通道,閥芯與繞線架圍成軸向圓環(huán)液流通道。工作時(shí),磁流變液首先通過左端蓋中心通孔,流經(jīng)左定位盤中心小孔流道,隨后通過左端徑向圓盤阻尼間隙,進(jìn)入軸向環(huán)形阻尼間隙,再通過右端徑向圓盤阻尼間隙,到達(dá)右定位盤中心小孔流道,最后從右端蓋中心通孔流出。

    如圖1所示,單線圈徑向流磁流變閥左右兩側(cè)翼套筒與隔磁套筒過盈配合構(gòu)成線圈繞線架,閥芯兩側(cè)端面與圓周表面分別設(shè)計(jì)有4個(gè)定位凸臺,約束其軸向與徑向位移。導(dǎo)磁材料較磁流變液擁有更低的磁阻,磁力線更容易在磁阻低的材料中形成閉合回路,為了能夠降低磁路中對有效阻尼間隙磁場強(qiáng)度貢獻(xiàn)較大的垂直分量的損耗,該磁流變閥在線圈繞線架中部設(shè)置了隔磁性零件隔磁套筒,同時(shí)左右兩端側(cè)翼套筒設(shè)置為導(dǎo)磁材料。利用隔磁材料的低導(dǎo)磁特性,引導(dǎo)磁力線更集中均勻地通過徑向圓盤阻尼間隙,有效提高了磁流變閥磁場利用率,在不增加外形尺寸的情況下,可進(jìn)一步提升磁流變閥的壓降性能。

    1.2 磁路設(shè)計(jì)

    圖2為單線圈徑向流磁流變閥磁路示意圖,磁力線垂直穿過的零件與區(qū)域包括閥體、側(cè)翼套筒、定位盤、閥芯以及兩段徑向阻尼間隙。假定磁流變閥磁路中磁力線均勻分布,無磁漏現(xiàn)象,則回路中有

    Φr=Φl=Φsteel=Φ

    (1)

    式中Φ——通電線圈總磁通量

    Φl——經(jīng)過左端徑向圓盤阻尼間隙磁通量

    Φr——經(jīng)過右端徑向圓盤阻尼間隙磁通量

    Φsteel——閥體、閥芯及定位盤等導(dǎo)磁零件磁通量

    磁流變閥磁路共分為8段有效部分,分別為兩段徑向圓盤阻尼間隙、閥體、繞線架左右兩端側(cè)翼套筒、左右兩端定位盤以及閥芯?;诖耍捎没鶢柣舴蚨蓪Υ怕愤M(jìn)行分析,可表示為

    (2)

    式中Nc——通電線圈匝數(shù)

    I——通電線圈電流

    Hi——磁路中第i部分的磁場強(qiáng)度

    li——磁力線經(jīng)過第i部分的有效長度

    圖2 單線圈徑向流磁流變閥磁路示意圖Fig.2 Magnetic circuit diagram of single coil radial flow MR valve

    通電線圈產(chǎn)生的磁通量可表示為

    Φ=∮cBdS=BiSi

    (3)

    式中Bi——磁路中第i部分的磁感應(yīng)強(qiáng)度

    Si——磁力線垂直經(jīng)過區(qū)域第i部分的橫截面面積

    磁感應(yīng)強(qiáng)度Bi與磁場強(qiáng)度Hi兩者之間的關(guān)系為

    Bi=μ0μiHi

    (4)

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

    μi——各區(qū)域材料的相對磁導(dǎo)率

    閉合磁路中有效區(qū)域磁阻Ri表示為

    (5)

    閥體磁阻R1為

    (6)

    側(cè)翼套筒及定位盤與其豎直相接部分磁阻R2和R8為

    (7)

    定位盤磁阻R3和R7為

    (8)

    徑向圓盤阻尼間隙磁阻R4和R6為

    (9)

    閥芯磁阻R5為

    (10)

    式中r——閥體半徑H——閥體長度

    th——閥體厚度

    b——側(cè)翼套筒寬度

    a——定位盤寬度

    gr——徑向圓盤阻尼間隙

    ga——軸向圓盤阻尼間隙

    Ha——閥芯寬度

    r1——閥芯半徑

    r0——小孔流道半徑

    μs——10號鋼的相對磁導(dǎo)率,取1 000

    μf——磁流變液的相對磁導(dǎo)率,取2.5

    因此,單線圈徑向流磁流變閥總磁阻Rm可表示為

    Rm=R1+2R2+2R3+2R4+R5

    (11)

    將式(3)、(4)代入式(2),可得

    (12)

    故磁路各組成部分的磁感應(yīng)強(qiáng)度Bi可表示為

    (13)

    1.3 壓降數(shù)學(xué)模型

    圖3為單線圈徑向流磁流變閥液流通道結(jié)構(gòu)圖。單線圈徑向流磁流變閥包含小孔流道、徑向流道與軸向流道3種不同類型的液流通道(圖2),忽略左右端蓋處的粘滯壓降,劃分為5處液流區(qū)域,其中Ⅰ、Ⅴ為小孔流道,Ⅱ、Ⅳ為徑向流道,Ⅲ為軸向流道。小孔流道Ⅰ、Ⅴ與軸向流道Ⅲ液流方向與磁力線平行,其壓降Δp1、Δp5、Δp3僅包含粘滯壓降,徑向流道Ⅱ、Ⅳ的液流方向與磁力線垂直,壓降Δp2、Δp4由粘滯壓降和磁致壓降兩部分構(gòu)成。單線圈徑向流磁流變閥壓降Δp可表示為

    Δp=Δp1+Δp2+Δp3+Δp4+Δp5=Δp1η+
    Δp2η+Δp2τ+Δp3η+Δp4η+Δp4τ+Δp5η

    (14)

    式中 Δp——單線圈徑向流磁流變閥總壓降

    Δp1ηΔp5η——液流區(qū)域Ⅰ~Ⅴ粘滯壓降

    Δp2τ、Δp4τ——區(qū)域Ⅱ、Ⅳ磁致壓降

    圖3 單線圈徑向流磁流變閥液流通道結(jié)構(gòu)圖Fig.3 Fluid channel structure diagram of single coil radial flow MR valve

    液流區(qū)域Ⅰ與液流區(qū)域Ⅴ為等徑等長小孔流道,其壓降Δp1和Δp5可表示為

    (15)

    式中q——系統(tǒng)流量,L/min

    η——磁流變液零場粘度,取0.8 Pa·s

    區(qū)域Ⅲ處軸向圓環(huán)阻尼間隙產(chǎn)生的壓降Δp3為

    (16)

    區(qū)域Ⅱ與區(qū)域Ⅳ處徑向圓盤阻尼間隙處產(chǎn)生的壓降Δp2和Δp4分為磁致壓降與粘滯壓降兩部分,其粘滯壓降為

    (17)

    磁致壓降可表示為

    (18)

    式中τy——?jiǎng)討B(tài)剪切屈服應(yīng)力

    c——修正系數(shù),無量綱,取2~3

    將式(15)~(18)代入式(14),得到單線圈徑向流磁流變閥總壓降Δp為

    (19)

    壓降可調(diào)范圍是衡量磁流變閥應(yīng)用范圍的重要指標(biāo),定義壓降可調(diào)系數(shù)K為Δpτ與Δpη的比值,可表示為

    (20)

    2 單線圈徑向流磁流變閥磁場與尺寸參數(shù)仿真分析

    2.1 仿真模型建立

    磁流變閥內(nèi)部磁路中的磁力線走向很大程度上取決于材料屬性的設(shè)置,通過隔磁材料與導(dǎo)磁材料的合理配合,可以減少磁路中的磁漏現(xiàn)象,引導(dǎo)磁力線均勻密集地穿過有效阻尼間隙。為了驗(yàn)證設(shè)計(jì)的單線圈徑向流磁流變閥結(jié)構(gòu)的合理性,了解關(guān)鍵尺寸參數(shù)對于其可控壓降性能的影響,采用有限元仿真軟件ANSYS的電磁場模塊建立單線圈徑向流磁流變閥1/2二維仿真實(shí)體模型(圖4)。為了提高建模與運(yùn)算的簡便快捷性,將單線圈徑向流磁流變閥三維結(jié)構(gòu)簡化成二維平面問題,考慮其截面為規(guī)則軸對稱圖形,在不影響模型求解精度的情況下,選取1/2截面為仿真對象。選擇1等級的SmartSize智能網(wǎng)格劃分以及四邊形網(wǎng)格結(jié)構(gòu)單元,施加無磁漏的磁力線平行邊界,并在線圈單元施加電流密度輸入勵(lì)磁載荷。

    圖4 單線圈徑向流磁流變閥2D仿真模型Fig.4 2D model of single coil radial flow MR valve1.液流通道 2.定位盤 3.閥芯 4.側(cè)翼套筒 5.隔磁套筒 6.勵(lì)磁線圈 7.閥體

    根據(jù)材料導(dǎo)磁特性的不同,二維模型各單元被劃分為隔磁材料、導(dǎo)磁材料、勵(lì)磁線圈與磁流變液4類屬性不同的區(qū)域,其中隔磁材料區(qū)域?yàn)檫x用不銹鋼材料的隔磁套筒,導(dǎo)磁材料區(qū)域由選用10號低碳鋼材料的定位盤、閥芯、側(cè)翼套筒以及閥體共同組成。采用可視為均勻材料,且相對磁導(dǎo)率為1的均布銅線作為勵(lì)磁線圈,采用重慶材料研究所開發(fā)的MRF-J01T型磁流變液作為磁流變閥的工作介質(zhì),使用最小二乘法對MRF-J01T型磁流變液τy-B曲線的非線性關(guān)系進(jìn)行三次多項(xiàng)式擬合[15-17],得到徑向圓盤阻尼間隙處剪切屈服應(yīng)力與磁感應(yīng)強(qiáng)度的關(guān)系為

    τy=a3B3+a2B2+a1B+a0

    (21)

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

    2.2 壓降性能仿真分析

    圖5為單線圈徑向流磁流變閥在加載電流為1 A時(shí)磁力線走向及路徑分布,圖中定義了徑向圓盤液流通道S1、S3與軸向圓環(huán)液流通道S2共3條路徑。由圖5可知,通過在繞線架中部設(shè)置不導(dǎo)磁隔磁套筒,絕大部分磁力線都能夠垂直穿過徑向阻尼間隙,避免了磁力線直接在繞線架處形成閉合回路現(xiàn)象的發(fā)生,提高了磁場利用率。同時(shí)還可觀察到,越靠近閥芯中部的磁力線其分布越集中均勻,走向也越垂直于徑向阻尼間隙,這也表明了磁流變閥磁路設(shè)計(jì)的合理性。

    圖5 單線圈徑向流磁流變閥磁力線分布Fig.5 Distribution of magnetic flux density

    圖6 徑向阻尼間隙內(nèi)磁感應(yīng)強(qiáng)度隨路徑分布曲線Fig.6 Distribution of magnetic flux density in radial damping gap along defined paths

    圖6為加載電流為1A時(shí)路徑S1、S3處磁場強(qiáng)度隨距離的變化情況。從圖6可看出,路徑S1、S3處的磁感應(yīng)強(qiáng)度大小及分布規(guī)律基本相同,其路徑兩端的磁感應(yīng)強(qiáng)度較低,沿路徑逐漸增大,位于路徑中間區(qū)段的磁感應(yīng)強(qiáng)度集中且穩(wěn)定。同時(shí)2條路徑的磁感應(yīng)強(qiáng)度均在0.453~0.495 T之間,接近MRF-J01T型磁流變液的飽和磁場強(qiáng)度,基本滿足磁流變閥對磁場參數(shù)的要求,為后續(xù)磁流變閥尺寸參數(shù)優(yōu)化提供了數(shù)據(jù)參考。

    (22)

    式中l(wèi)——液流通道長度

    B⊥——垂直阻尼間隙方向的磁感應(yīng)強(qiáng)度

    由式(22)可知,提高磁感應(yīng)強(qiáng)度垂直阻尼間隙方向的分量可增加相應(yīng)阻尼間隙內(nèi)的平均磁感應(yīng)強(qiáng)度。

    圖7為單線圈徑向流磁流變閥壓降以及阻尼間隙內(nèi)平均磁感應(yīng)強(qiáng)度隨加載電流的變化曲線。從圖7可看出,平均磁感應(yīng)強(qiáng)度與壓降都隨輸入電流的增加而增大,因?yàn)檩斎腚娏髟龃髸鸫怕分写艅?dòng)勢增加,因此間隙內(nèi)平均磁感應(yīng)強(qiáng)度增大;磁感應(yīng)強(qiáng)度的增加引起了磁流變液剪切屈服應(yīng)力的提高,從而導(dǎo)致磁流變閥壓降提升。此外,在輸入電流為1 A時(shí),磁流變閥可產(chǎn)生約2 MPa壓降。

    圖7 不同電流下的平均磁感應(yīng)強(qiáng)度及壓降變化曲線Fig.7 Changing curves of average magnetic flux density and pressure drop under different currents

    2.3 尺寸參數(shù)對性能的影響

    關(guān)鍵部件尺寸參數(shù)對提升磁流變閥可控壓降性能具有重要影響,通過改變徑向圓盤阻尼間隙gr和定位盤寬度a來探究磁流變閥壓降與可調(diào)系數(shù)變化狀況。選取約束條件為

    Ha=H-2gr-2a

    (23)

    圖8為在約束條件基礎(chǔ)上,只改變徑向圓盤阻尼間隙gr得到的磁感應(yīng)強(qiáng)度與壓降可調(diào)系數(shù)變化曲線。從圖8可看出,隨著徑向圓盤阻尼間隙增加,磁流變閥有效阻尼間隙內(nèi)的磁感應(yīng)強(qiáng)度隨之減小,這主要是因?yàn)殚g隙的增大必然會帶來磁路磁阻增加。此外,間隙的增加導(dǎo)致了壓降可調(diào)系數(shù)的增加,由式(19)可知,由于磁致壓降與圓盤阻尼間隙一次方呈反比,而粘滯壓降與阻尼間隙三次方呈反比,粘滯壓降的下降速度較磁致壓降更快。

    圖8 不同阻尼間隙下的磁感應(yīng)強(qiáng)度及壓降可調(diào)系數(shù)Fig.8 Magnetic flux density and pressure drop adjustable coefficient under different damping gaps

    圖9為磁感應(yīng)強(qiáng)度與壓降可調(diào)系數(shù)隨定位盤寬度a的變化曲線。由圖9可知,隨著定位盤寬度的增加,阻尼間隙內(nèi)的磁感應(yīng)強(qiáng)度基本不變,但壓降可調(diào)系數(shù)明顯上升。因?yàn)樵诩s束條件下,定位盤寬度的增加減少了軸向圓環(huán)液流通道的長度,粘滯壓降會有所下降,而整體磁路磁阻基本不變,因此磁致壓降不變,從而導(dǎo)致磁流變閥壓降可調(diào)系數(shù)增加。

    圖9 不同定位盤寬度下的磁感應(yīng)強(qiáng)度及壓降可調(diào)系數(shù)Fig.9 Magnetic flux density and pressure drop adjustable coefficient under different location width

    3 單線圈徑向流磁流變閥優(yōu)化

    3.1 可控壓降性能優(yōu)化

    利用ANSYS軟件中的優(yōu)化工具改善相關(guān)關(guān)鍵部件尺寸,進(jìn)一步提高單線圈徑向流磁流變閥性能效果。僅使用一階優(yōu)化方法進(jìn)行仿真求解易產(chǎn)生局部最優(yōu)解,而零階優(yōu)化方法求解精度不高,因此綜合運(yùn)用上述兩種優(yōu)化工具進(jìn)行優(yōu)化仿真。首先利用零階方法確定最優(yōu)解近似區(qū)域,然后通過一階方法精確求解得到最優(yōu)參數(shù),該方法能夠在得到準(zhǔn)確結(jié)果的同時(shí),減少優(yōu)化算法的迭代時(shí)間,更加快速地得到最佳技術(shù)方案,其優(yōu)化流程如圖10所示。在初步優(yōu)化設(shè)計(jì)中,為了避免在非合理設(shè)計(jì)空間上得到局部最小值,首先使用零階方法進(jìn)行曲線擬合逼近合理設(shè)計(jì)序列,循環(huán)結(jié)束后以收斂值為起點(diǎn)運(yùn)行一階方法,添加懲罰函數(shù)得到無約束目標(biāo)函數(shù)后,搜尋設(shè)計(jì)變量方向矢量計(jì)算得到新的設(shè)計(jì)變量,最終經(jīng)過多次迭代后得到最優(yōu)解。

    圖10 單線圈徑向流磁流變閥優(yōu)化設(shè)計(jì)流程圖Fig.10 Flow chart of optimal design

    性能優(yōu)異的磁流變閥應(yīng)在擁有較高壓降的同時(shí)兼具較寬的壓降可調(diào)范圍,為了實(shí)現(xiàn)相應(yīng)性能指標(biāo)的多目標(biāo)優(yōu)化,采用線性加權(quán)法設(shè)立指標(biāo)權(quán)重系數(shù)構(gòu)建目標(biāo)函數(shù)為

    (24)

    式中p′τ——優(yōu)化前磁致壓降

    pτ——優(yōu)化后磁致壓降

    K′——優(yōu)化前壓降可調(diào)系數(shù)

    K——優(yōu)化后壓降可調(diào)系數(shù)

    λ1——磁致壓降權(quán)重系數(shù)

    λ2——壓降可調(diào)系數(shù)的權(quán)重系數(shù)

    表1 磁流變閥優(yōu)化設(shè)計(jì)變量

    Tab.1 Optimization design variable of MR valvemm

    參數(shù)數(shù)值小孔流道半徑r03~6閥芯半徑r112~20定位盤寬度a8~16閥體厚度th6~10繞線槽寬T25~40徑向圓盤阻尼間隙gr1~2

    考慮磁流變閥綜合性能對于磁致壓降與壓降可調(diào)系數(shù)兩項(xiàng)指標(biāo)的需求,設(shè)置了3種不同權(quán)重系數(shù)λ1=1、λ2=0;λ1=0、λ2=1以及λ1=0.75、λ2=0.25的目標(biāo)函數(shù)分別進(jìn)行優(yōu)化。圖11為目標(biāo)函數(shù)權(quán)重系數(shù)λ1=0.75、λ2=0.25的磁流變閥相關(guān)優(yōu)化設(shè)計(jì)參數(shù)迭代圖,從圖11b可看出,優(yōu)化過程共經(jīng)歷12次迭代,目標(biāo)函數(shù)在第12次迭代中獲得最小值,即得到最優(yōu)結(jié)果。圖11a中設(shè)計(jì)變量與狀態(tài)變量參數(shù)在經(jīng)歷最初較大幅值變化后逐漸趨于平穩(wěn),因?yàn)榍捌诹汶A優(yōu)化在初始變量設(shè)計(jì)上下限內(nèi)得到小范圍最優(yōu)解近似區(qū)域,一階優(yōu)化在此基礎(chǔ)上進(jìn)而獲得準(zhǔn)確最優(yōu)結(jié)果。最終優(yōu)化結(jié)果如表2所示。

    圖11 單線圈徑向流磁流變閥多目標(biāo)優(yōu)化結(jié)果Fig.11 Optimization result of proposed MR valve

    參數(shù)優(yōu)化前優(yōu)化后λ1=1λ2=0λ1=0λ2=1λ1=0.75λ2=0.25目標(biāo)函數(shù)f1.000.3980.5670.685磁致壓降pτ/MPa1.9094.7981.5153.801粘滯壓降pη/MPa0.0920.0400.0410.227壓降可調(diào)系數(shù)K20.7112.1036.5116.77平均磁感應(yīng)強(qiáng)度B/T0.4890.5260.4780.518小孔流道半徑r0/mm4.53.04.94.2閥芯半徑r1/mm15.020.016.219.6定位盤寬度a/mm1213.951616閥體厚度th/mm8666繞線槽寬T/mm34404040徑向圓盤阻尼間隙gr/mm1.51.02.01.1

    由表2可知,當(dāng)權(quán)重系數(shù)λ1=1、λ2=0時(shí),磁致壓降pτ由1.909 MPa增加到4.798 MPa,提高了151.3%,壓降可調(diào)系數(shù)K由20.71下降至12.10,下降41.2%。當(dāng)權(quán)重系數(shù)λ1=0、λ2=1時(shí),磁致壓降pτ由1.909 MPa下降至1.515 MPa,下降20.7%,壓降可調(diào)系數(shù)K由20.71增加到36.51,提高了76.3%。當(dāng)權(quán)重系數(shù)λ1=0.75、λ2=0.25時(shí),磁致壓降pτ由1.909 MPa增加到3.801 MPa,提高了99.1%,壓降可調(diào)系數(shù)K由20.71下降至16.77,下降19.0%。從上述優(yōu)化結(jié)果可看出,權(quán)重系數(shù)為λ1=0.75、λ2=0.25,輸入電流為1 A時(shí),磁流變閥犧牲較少的壓降可調(diào)系數(shù)性能,其磁致壓降pτ增加約2.9 MPa,性能提升顯著,避免了單性能指標(biāo)優(yōu)化導(dǎo)致另一性能大幅度下降的問題。

    3.2 優(yōu)化結(jié)果分析

    圖12 優(yōu)化前后平均磁感應(yīng)強(qiáng)度隨電流的變化曲線Fig.12 Changing curves of average magnetic flux density before and after optimization

    圖13 優(yōu)化前后壓降隨電流的變化曲線Fig.13 Changing curves of pressure drop before and after optimization

    圖13為優(yōu)化前后單線圈徑向流磁流變閥壓降隨加載電流的變化仿真曲線。從圖13可看出,優(yōu)化前后磁流變閥壓降隨輸入電流的增加而上升。同一輸入電流下,磁流變閥優(yōu)化前后阻尼間隙內(nèi)平均磁感應(yīng)強(qiáng)度差值較小,優(yōu)化后磁流變閥壓降相比優(yōu)化前壓降有明顯提升。因?yàn)榇胖聣航凳芗羟袘?yīng)力與尺寸共同影響,由式(18)、(21)可知,優(yōu)化前后磁流變液動(dòng)態(tài)剪切屈服應(yīng)力變化較小,但磁致壓降式(18)的系數(shù)從優(yōu)化前的7增加到優(yōu)化后的14,因此優(yōu)化后磁流變閥壓降提升較為顯著。仿真結(jié)果顯示,當(dāng)加載電流為1 A時(shí),優(yōu)化后的磁流變閥壓降可達(dá)4 MPa,相比優(yōu)化前的2 MPa,性能提升100%,優(yōu)化效果明顯。

    4 壓降性能實(shí)驗(yàn)

    4.1 磁流變閥樣機(jī)及性能測試系統(tǒng)

    在仿真分析及優(yōu)化設(shè)計(jì)基礎(chǔ)上,加工了圖14所示的優(yōu)化前和優(yōu)化后2種內(nèi)部尺寸不同的單線圈徑向流磁流變閥。設(shè)計(jì)優(yōu)化前后磁流變閥性能測試方法,搭建了圖15所示單線圈徑向流磁流變閥性能測試系統(tǒng)。為了更好地將具有一定壓力和流量的磁流變液輸送到磁流變閥中,液壓驅(qū)動(dòng)系統(tǒng)采用了由三相異步電動(dòng)機(jī)以及齒輪式定量泵組成的動(dòng)力單元,其中定量泵系統(tǒng)流量為4 L/min;常閉狀態(tài)的溢流閥Ⅰ作為安全閥安裝在齒輪泵與磁流變液缸之間,安全閥值為6.3 MPa;溢流閥Ⅱ作為背壓閥安裝在磁流變閥與磁流變液缸回油口之間,通過調(diào)整閥口開度來模擬系統(tǒng)負(fù)載變化。控制檢測系統(tǒng)中的壓力傳感器Ⅰ和Ⅱ分別檢測磁流變閥進(jìn)口和出口壓力,并將壓力信號轉(zhuǎn)換為模擬信號,經(jīng)由數(shù)據(jù)采集卡傳送到計(jì)算機(jī)上的圖形化編程語言LabVIEW中進(jìn)行計(jì)算處理。

    圖14 單線圈徑向流磁流變閥實(shí)物圖Fig.14 Prototyping of single coil radial flow MR valve1.繞線架 2.閥體 3.左端蓋 4.右端蓋 5.閥芯 6.定位盤 7.優(yōu)化前磁流變閥 8.優(yōu)化后磁流變閥

    4.2 優(yōu)化后單線圈徑向流磁流變閥性能

    手動(dòng)旋轉(zhuǎn)溢流閥Ⅱ的手柄可以改變磁流變閥出口負(fù)載壓力,可以模擬系統(tǒng)負(fù)載工況的變化。定義溢流閥Ⅱ手柄未旋轉(zhuǎn),閥口處于最大開度狀態(tài)時(shí)為負(fù)載工況1;順時(shí)針旋轉(zhuǎn)溢流閥Ⅱ手柄2圈為負(fù)載工況2,順時(shí)針旋轉(zhuǎn)溢流閥Ⅱ手柄4圈為負(fù)載工況3。

    圖16為優(yōu)化后單線圈徑向流磁流變閥在負(fù)載工況1下進(jìn)口壓力、出口壓力以及壓降隨電流的變化曲線。從圖16可看出,隨著加載電流的增加,進(jìn)口壓力與壓降均以相同的上升趨勢不斷增加,與仿真結(jié)果趨勢一致。當(dāng)輸入電流達(dá)到1.8 A時(shí),可得到約2.58 MPa的最大壓降。溢流閥Ⅱ作為背壓閥可保證出口壓力基本保持恒定,但出口壓力仍有小幅波動(dòng),因?yàn)橄到y(tǒng)在工作過程中磁流變閥中磁流變液流速會因齒輪泵磨損以及電動(dòng)機(jī)振動(dòng)而不穩(wěn)定。

    圖16 優(yōu)化后磁流變閥壓力在負(fù)載工況1下變化趨勢Fig.16 Pressure of MR valve after optimization at load condition 1

    圖17為負(fù)載變化對優(yōu)化后單線圈徑向流磁流變閥壓降性能的影響。由圖17可知,不同負(fù)載工況下的壓降曲線基本重合,這說明優(yōu)化后磁流變閥壓降性能不受負(fù)載變化的影響,這種負(fù)載變化不改變壓降的特性可以進(jìn)一步拓寬磁流變閥的使用范圍,尤其是將磁流變閥作為旁通閥并聯(lián)在阻尼器旁,可改善阻尼器的阻尼力調(diào)節(jié)范圍,提升磁流變閥控阻尼器設(shè)計(jì)的靈活性,使其能夠應(yīng)用于不同阻尼工況下。

    圖17 優(yōu)化后磁流變閥不同負(fù)載工況下壓降比較Fig.17 Pressure drop comparison of MR valve under different load conditions

    4.3 磁流變閥壓降性能對比

    圖18為優(yōu)化前后單線圈徑向流磁流變閥壓降隨電流變化的對比曲線。與圖13優(yōu)化前后磁流變閥壓降仿真曲線相比,實(shí)驗(yàn)壓降與仿真壓降具有相同的變化規(guī)律。從圖18可看出,優(yōu)化后磁流變閥比優(yōu)化前擁有更優(yōu)異的壓降性能,在輸入電流為0.3 A時(shí),優(yōu)化后磁流變閥壓降相比優(yōu)化前開始明顯增大,在加載電流1.8 A時(shí),優(yōu)化前磁流變閥總壓降為1.84 MPa,優(yōu)化后磁流變閥總壓降可達(dá)2.58 MPa,提升40.22%。可控壓降性能優(yōu)化效果顯著。

    圖18 優(yōu)化前后磁流變閥壓降性能對比Fig.18 Comparison of pressure drop of MR valve before and after optimization

    表3為優(yōu)化前后單線圈徑向流磁流變閥壓降與壓降可調(diào)系數(shù)。從表3可看出,優(yōu)化后磁流變閥壓降性能有較大改善,壓降可調(diào)系數(shù)也有較大提升。優(yōu)化前磁流變閥壓降可調(diào)系數(shù)為7.94,優(yōu)化后壓降可調(diào)系數(shù)可達(dá)10.07,提升26.83%。與優(yōu)化仿真結(jié)果相比有一定差異,原因?yàn)椋孩僭诖艌龇抡嬷胁捎美硐脒吔鐥l件,未考慮實(shí)際磁流變閥存在的漏磁現(xiàn)象。②采用Bingham模型推導(dǎo)壓降數(shù)學(xué)模型未考慮磁流變液剪切變稠現(xiàn)象的發(fā)生。③由于實(shí)驗(yàn)條件的限制,未對實(shí)際材料特性參數(shù)進(jìn)行測量,導(dǎo)致仿真設(shè)定的參數(shù)與實(shí)際情況有差異。

    表3 優(yōu)化前后磁流變閥壓降及壓降可調(diào)系數(shù)Tab.3 Pressure drop and adjustable coefficient of MR valve before and after optimization

    5 結(jié)論

    (1)基于傳統(tǒng)徑向流磁流變閥,設(shè)計(jì)了一種在線圈繞線架中部設(shè)置隔磁套筒的單線圈徑向流磁流變閥,合理改善了磁路設(shè)計(jì)。結(jié)合推導(dǎo)的壓降數(shù)學(xué)模型,對單線圈徑向流磁流變閥進(jìn)行了建模仿真。仿真結(jié)果表明,磁流變閥初始結(jié)構(gòu)設(shè)計(jì)及其基本參數(shù)設(shè)定合理,基于此進(jìn)行優(yōu)化設(shè)計(jì)獲得的最優(yōu)解正確合理。

    (2)運(yùn)用零階與一階優(yōu)化工具,以權(quán)重系數(shù)λ1=0.75、λ2=0.25對磁流變閥建立的多目標(biāo)函數(shù)進(jìn)行優(yōu)化分析。結(jié)果表明,在加載電流為1 A時(shí),磁流變閥壓降由優(yōu)化前的2 MPa增至優(yōu)化后的4 MPa,性能提高了100%;優(yōu)化后磁流變閥在犧牲較少壓降可調(diào)系數(shù)性能的情況下,大幅提升了磁致壓降,性能提升顯著。

    (3)對優(yōu)化后單線圈徑向流磁流變閥壓降性能進(jìn)行實(shí)驗(yàn)測試對比,實(shí)驗(yàn)結(jié)果表明,在勵(lì)磁電流1.8 A時(shí),優(yōu)化后磁流變閥總壓降為2.58 MPa,較優(yōu)化前提升了40.22%;優(yōu)化后磁流變閥壓降可調(diào)系數(shù)為10.07,較優(yōu)化前提升了26.83%。

    猜你喜歡
    磁路磁感應(yīng)阻尼
    基于分布磁路法的非晶合金高速磁浮直線電機(jī)懸浮力計(jì)算
    跨空海界面磁感應(yīng)通信特性分析及應(yīng)用
    電磁感應(yīng)中的“知三求三”
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項(xiàng)的擴(kuò)散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    基于等效磁路法的永磁同步電機(jī)特性分析
    具阻尼項(xiàng)的Boussinesq型方程的長時(shí)間行為
    基于ANSYS Workbench的微型揚(yáng)聲器磁路優(yōu)化分析
    永磁魔環(huán)磁感應(yīng)強(qiáng)度的仿真模擬
    物理與工程(2013年3期)2013-03-11 16:04:35
    精品电影一区二区在线| 亚洲av日韩精品久久久久久密| 成人国产综合亚洲| 不卡av一区二区三区| 婷婷精品国产亚洲av| 亚洲熟女毛片儿| 全区人妻精品视频| 日韩国内少妇激情av| 亚洲成av人片在线播放无| 又黄又爽又免费观看的视频| 国产精品98久久久久久宅男小说| 18禁观看日本| 91av网一区二区| 日韩欧美三级三区| 国产精品精品国产色婷婷| 叶爱在线成人免费视频播放| 亚洲九九香蕉| 国产精品香港三级国产av潘金莲| 午夜影院日韩av| 一区福利在线观看| 99国产极品粉嫩在线观看| 两个人的视频大全免费| 12—13女人毛片做爰片一| 国产高潮美女av| 男人舔女人的私密视频| 午夜福利18| 国产不卡一卡二| 九色成人免费人妻av| 国产成人精品无人区| 精品免费久久久久久久清纯| 亚洲国产欧美一区二区综合| 亚洲国产色片| 一进一出抽搐动态| 狂野欧美白嫩少妇大欣赏| 国产成人精品久久二区二区91| 国产成人精品无人区| 日韩av在线大香蕉| 别揉我奶头~嗯~啊~动态视频| 欧美日韩一级在线毛片| 国产成人影院久久av| 在线播放国产精品三级| 男女那种视频在线观看| 亚洲国产欧美网| 一区二区三区激情视频| 女人高潮潮喷娇喘18禁视频| 久久亚洲精品不卡| 99re在线观看精品视频| 国产精品久久久人人做人人爽| 五月伊人婷婷丁香| 日韩欧美一区二区三区在线观看| 美女免费视频网站| 麻豆av在线久日| 最近最新中文字幕大全免费视频| a在线观看视频网站| 久久久色成人| 美女 人体艺术 gogo| 午夜精品在线福利| 国产精品乱码一区二三区的特点| 日韩欧美在线乱码| 国产精品永久免费网站| 国产高清videossex| 亚洲av中文字字幕乱码综合| 9191精品国产免费久久| 每晚都被弄得嗷嗷叫到高潮| 国产成年人精品一区二区| 国产成人精品无人区| 亚洲成人久久爱视频| 啦啦啦观看免费观看视频高清| 色老头精品视频在线观看| 亚洲最大成人中文| 国产精品久久久久久久电影 | 欧美另类亚洲清纯唯美| 亚洲一区高清亚洲精品| x7x7x7水蜜桃| 18禁黄网站禁片午夜丰满| 日本黄色片子视频| 国内久久婷婷六月综合欲色啪| 一级毛片女人18水好多| 在线观看舔阴道视频| 曰老女人黄片| 日韩免费av在线播放| 色综合亚洲欧美另类图片| 琪琪午夜伦伦电影理论片6080| 一夜夜www| tocl精华| 性欧美人与动物交配| 欧美乱色亚洲激情| 99re在线观看精品视频| a级毛片在线看网站| www.熟女人妻精品国产| 亚洲中文字幕日韩| 国产成年人精品一区二区| 久久久国产成人免费| 国产欧美日韩精品一区二区| 日本与韩国留学比较| 亚洲av第一区精品v没综合| 熟女少妇亚洲综合色aaa.| 白带黄色成豆腐渣| 精品乱码久久久久久99久播| 久久久色成人| 在线国产一区二区在线| 成人三级黄色视频| 天天躁日日操中文字幕| 最新美女视频免费是黄的| 国内精品美女久久久久久| 成年女人永久免费观看视频| 国产97色在线日韩免费| 国产乱人伦免费视频| 香蕉av资源在线| 欧美色视频一区免费| 中文字幕高清在线视频| 啦啦啦免费观看视频1| 日本 欧美在线| 国产精品一区二区免费欧美| 97碰自拍视频| 狠狠狠狠99中文字幕| 亚洲第一电影网av| 老汉色∧v一级毛片| 久久久久久国产a免费观看| aaaaa片日本免费| 黑人操中国人逼视频| 美女免费视频网站| 嫁个100分男人电影在线观看| 人妻丰满熟妇av一区二区三区| 成人国产一区最新在线观看| 免费高清视频大片| 听说在线观看完整版免费高清| 精品熟女少妇八av免费久了| 亚洲欧美精品综合久久99| 天堂网av新在线| 在线观看免费视频日本深夜| 在线视频色国产色| 亚洲九九香蕉| x7x7x7水蜜桃| 亚洲av成人精品一区久久| 亚洲va日本ⅴa欧美va伊人久久| 成人特级黄色片久久久久久久| 成人三级黄色视频| 国产成人福利小说| 欧美xxxx黑人xx丫x性爽| 成人鲁丝片一二三区免费| 国产探花在线观看一区二区| 亚洲色图av天堂| 女人高潮潮喷娇喘18禁视频| 亚洲精品色激情综合| 又大又爽又粗| 成人鲁丝片一二三区免费| 日日夜夜操网爽| 99视频精品全部免费 在线 | 欧美中文日本在线观看视频| 香蕉久久夜色| 一级黄色大片毛片| 啦啦啦观看免费观看视频高清| 婷婷精品国产亚洲av| 国产一区在线观看成人免费| 国产成人精品久久二区二区91| 91麻豆av在线| 少妇熟女aⅴ在线视频| 搡老岳熟女国产| 性色av乱码一区二区三区2| 成年免费大片在线观看| 日韩国内少妇激情av| 国产精品香港三级国产av潘金莲| 色播亚洲综合网| 亚洲欧洲精品一区二区精品久久久| 国产毛片a区久久久久| 国产主播在线观看一区二区| 欧美zozozo另类| 欧美日韩一级在线毛片| 女人被狂操c到高潮| 99精品在免费线老司机午夜| 99久久精品国产亚洲精品| 伦理电影免费视频| 欧美不卡视频在线免费观看| 此物有八面人人有两片| 国产成年人精品一区二区| 老鸭窝网址在线观看| 亚洲人与动物交配视频| 欧洲精品卡2卡3卡4卡5卡区| 国产伦在线观看视频一区| 日韩人妻高清精品专区| 亚洲av美国av| 在线观看午夜福利视频| 99热这里只有精品一区 | 小说图片视频综合网站| 99re在线观看精品视频| 香蕉丝袜av| 国产成人一区二区三区免费视频网站| 日日夜夜操网爽| 999久久久国产精品视频| 国产三级黄色录像| 女人被狂操c到高潮| 欧美乱色亚洲激情| 国产精品,欧美在线| 精品无人区乱码1区二区| 激情在线观看视频在线高清| 琪琪午夜伦伦电影理论片6080| 午夜福利成人在线免费观看| 老汉色∧v一级毛片| 搡老妇女老女人老熟妇| 精品国内亚洲2022精品成人| 一区福利在线观看| 中文字幕精品亚洲无线码一区| 亚洲av五月六月丁香网| 一个人免费在线观看的高清视频| 99久久成人亚洲精品观看| 成人国产一区最新在线观看| 又大又爽又粗| 久久精品亚洲精品国产色婷小说| 亚洲精品色激情综合| 欧美中文日本在线观看视频| 99久久久亚洲精品蜜臀av| 男女之事视频高清在线观看| 国产激情欧美一区二区| 日本 av在线| 真人一进一出gif抽搐免费| 欧美绝顶高潮抽搐喷水| 波多野结衣高清作品| 精品久久蜜臀av无| 久久久久久人人人人人| 母亲3免费完整高清在线观看| 天堂av国产一区二区熟女人妻| 最新中文字幕久久久久 | 国产精品久久久久久久电影 | 狂野欧美激情性xxxx| 不卡一级毛片| 久久亚洲真实| 好看av亚洲va欧美ⅴa在| 亚洲av美国av| 欧美中文综合在线视频| 男人舔奶头视频| 丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品九九99| www.精华液| 欧美+亚洲+日韩+国产| 中国美女看黄片| 国产熟女xx| 男女床上黄色一级片免费看| 人妻久久中文字幕网| 成人av一区二区三区在线看| 最近最新免费中文字幕在线| 国产免费男女视频| 视频区欧美日本亚洲| 国产高清videossex| 女人高潮潮喷娇喘18禁视频| 亚洲色图 男人天堂 中文字幕| 色老头精品视频在线观看| 欧美一区二区国产精品久久精品| 级片在线观看| www日本黄色视频网| 国产麻豆成人av免费视频| 三级男女做爰猛烈吃奶摸视频| 叶爱在线成人免费视频播放| 日日摸夜夜添夜夜添小说| 国产伦精品一区二区三区视频9 | 中文在线观看免费www的网站| 色在线成人网| 99久久国产精品久久久| 国产一级毛片七仙女欲春2| 欧美在线一区亚洲| 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 男女那种视频在线观看| 久久婷婷人人爽人人干人人爱| 成年免费大片在线观看| 国产乱人视频| 久久久色成人| 国产亚洲精品久久久久久毛片| 久久中文字幕人妻熟女| 一级毛片高清免费大全| 天堂动漫精品| 国产又黄又爽又无遮挡在线| 国产成人系列免费观看| 97超级碰碰碰精品色视频在线观看| 亚洲精品在线美女| 一边摸一边抽搐一进一小说| 激情在线观看视频在线高清| 一个人免费在线观看电影 | 亚洲av熟女| 亚洲午夜理论影院| 女生性感内裤真人,穿戴方法视频| 中文字幕最新亚洲高清| 国产高清激情床上av| 性欧美人与动物交配| 亚洲 国产 在线| 美女cb高潮喷水在线观看 | 精品欧美国产一区二区三| 网址你懂的国产日韩在线| 偷拍熟女少妇极品色| 免费观看精品视频网站| 深夜精品福利| 国产人伦9x9x在线观看| 亚洲片人在线观看| 国产v大片淫在线免费观看| 嫩草影视91久久| 国产欧美日韩一区二区精品| 亚洲成人免费电影在线观看| 在线观看美女被高潮喷水网站 | 国产午夜福利久久久久久| 波多野结衣高清作品| 性欧美人与动物交配| 亚洲午夜精品一区,二区,三区| 国产精品一区二区三区四区久久| 日本撒尿小便嘘嘘汇集6| 国产伦在线观看视频一区| 男人舔女人下体高潮全视频| 男人的好看免费观看在线视频| 可以在线观看毛片的网站| 久久人人精品亚洲av| 中国美女看黄片| 成人欧美大片| 18禁黄网站禁片免费观看直播| 这个男人来自地球电影免费观看| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区视频在线观看免费| АⅤ资源中文在线天堂| 成人无遮挡网站| 国产激情偷乱视频一区二区| 免费观看人在逋| 蜜桃久久精品国产亚洲av| 亚洲一区二区三区色噜噜| 噜噜噜噜噜久久久久久91| 欧美极品一区二区三区四区| 一区二区三区激情视频| or卡值多少钱| 国产久久久一区二区三区| 亚洲 国产 在线| 久久亚洲真实| 少妇的逼水好多| 日韩欧美一区二区三区在线观看| 操出白浆在线播放| 亚洲av五月六月丁香网| 久久中文看片网| 少妇的丰满在线观看| 999精品在线视频| 一夜夜www| 看免费av毛片| 97人妻精品一区二区三区麻豆| 中亚洲国语对白在线视频| 国产成人精品久久二区二区免费| 啦啦啦免费观看视频1| 欧美乱妇无乱码| 久久久久免费精品人妻一区二区| 在线永久观看黄色视频| 熟女人妻精品中文字幕| 成人国产综合亚洲| 亚洲精品国产精品久久久不卡| 午夜福利成人在线免费观看| 悠悠久久av| 最近最新中文字幕大全免费视频| 日韩大尺度精品在线看网址| 五月伊人婷婷丁香| 久久久国产成人免费| 国产毛片a区久久久久| 久久精品综合一区二区三区| 国产激情久久老熟女| ponron亚洲| 久久精品91无色码中文字幕| 91在线观看av| 久久这里只有精品中国| 亚洲成人久久爱视频| 欧美日韩瑟瑟在线播放| 久久久成人免费电影| 操出白浆在线播放| 激情在线观看视频在线高清| 九色成人免费人妻av| 欧美日韩亚洲国产一区二区在线观看| 老汉色∧v一级毛片| 狠狠狠狠99中文字幕| 一级毛片女人18水好多| 男人舔女人下体高潮全视频| 女警被强在线播放| 小说图片视频综合网站| 黄色视频,在线免费观看| 久久久久久久久免费视频了| 亚洲国产精品sss在线观看| 精品一区二区三区视频在线 | 亚洲欧美日韩东京热| 久久香蕉国产精品| 夜夜躁狠狠躁天天躁| 丰满人妻一区二区三区视频av | 黑人操中国人逼视频| 亚洲狠狠婷婷综合久久图片| 搞女人的毛片| 国产精品久久久久久人妻精品电影| 欧美不卡视频在线免费观看| 又大又爽又粗| 我的老师免费观看完整版| 国产淫片久久久久久久久 | 婷婷丁香在线五月| 人人妻人人看人人澡| 九色国产91popny在线| 变态另类成人亚洲欧美熟女| 国产91精品成人一区二区三区| 女警被强在线播放| 99视频精品全部免费 在线 | 久久久久久大精品| 国产人伦9x9x在线观看| 国产成人系列免费观看| 一进一出抽搐动态| 亚洲人成网站在线播放欧美日韩| 日本a在线网址| 99国产精品99久久久久| 国产欧美日韩一区二区三| 老司机午夜福利在线观看视频| 波多野结衣高清无吗| 久久久国产成人免费| 三级男女做爰猛烈吃奶摸视频| 狠狠狠狠99中文字幕| 亚洲av片天天在线观看| 国产成年人精品一区二区| 国产精品久久电影中文字幕| 亚洲精品一区av在线观看| 国产淫片久久久久久久久 | 夜夜爽天天搞| 国产午夜精品论理片| 久久久久久国产a免费观看| 淫妇啪啪啪对白视频| 男女做爰动态图高潮gif福利片| 国产成人欧美在线观看| 村上凉子中文字幕在线| 亚洲人与动物交配视频| 日韩免费av在线播放| 老司机午夜十八禁免费视频| 欧美在线黄色| av黄色大香蕉| 国产一区二区在线av高清观看| 中文字幕av在线有码专区| 国产精品女同一区二区软件 | 两性夫妻黄色片| 国产精品综合久久久久久久免费| 99久久精品国产亚洲精品| av国产免费在线观看| 欧美国产日韩亚洲一区| 久久精品影院6| 99久久99久久久精品蜜桃| 性色av乱码一区二区三区2| h日本视频在线播放| 国产精华一区二区三区| 97碰自拍视频| 欧美极品一区二区三区四区| 亚洲av免费在线观看| 美女被艹到高潮喷水动态| 极品教师在线免费播放| 757午夜福利合集在线观看| 欧美xxxx黑人xx丫x性爽| 身体一侧抽搐| 亚洲五月婷婷丁香| 国产精品99久久99久久久不卡| 国产综合懂色| 国产精品1区2区在线观看.| 久久久久久大精品| 真实男女啪啪啪动态图| 一区福利在线观看| 他把我摸到了高潮在线观看| 又爽又黄无遮挡网站| www.自偷自拍.com| 麻豆成人午夜福利视频| 欧美三级亚洲精品| 18美女黄网站色大片免费观看| 又黄又爽又免费观看的视频| 极品教师在线免费播放| 精品久久久久久久人妻蜜臀av| 免费高清视频大片| 亚洲一区二区三区不卡视频| 久久久久久国产a免费观看| svipshipincom国产片| 一区福利在线观看| 国产激情久久老熟女| 99久久精品热视频| www.熟女人妻精品国产| 久久99热这里只有精品18| 色av中文字幕| 亚洲国产精品sss在线观看| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 亚洲一区高清亚洲精品| 欧美日本亚洲视频在线播放| 国产成人一区二区三区免费视频网站| 日韩人妻高清精品专区| 很黄的视频免费| 亚洲美女黄片视频| 成人鲁丝片一二三区免费| 午夜福利18| 色视频www国产| 国产高清视频在线播放一区| 操出白浆在线播放| 日韩人妻高清精品专区| 久久久久性生活片| 欧美日本亚洲视频在线播放| 欧美乱妇无乱码| 日本黄色视频三级网站网址| 青草久久国产| 成人特级黄色片久久久久久久| 久久性视频一级片| 1000部很黄的大片| 禁无遮挡网站| 日本黄色视频三级网站网址| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 国产精品自产拍在线观看55亚洲| 亚洲中文日韩欧美视频| 哪里可以看免费的av片| 一级毛片精品| 免费看日本二区| 午夜福利在线在线| 久久精品国产99精品国产亚洲性色| 18禁黄网站禁片免费观看直播| 特大巨黑吊av在线直播| 一级a爱片免费观看的视频| 成人三级做爰电影| 男女下面进入的视频免费午夜| 国内久久婷婷六月综合欲色啪| 国产欧美日韩一区二区三| 九色国产91popny在线| 国产av不卡久久| 一级毛片高清免费大全| 麻豆国产97在线/欧美| 美女被艹到高潮喷水动态| 欧美一区二区国产精品久久精品| 看片在线看免费视频| 亚洲中文av在线| 国产亚洲精品综合一区在线观看| 欧美性猛交╳xxx乱大交人| 91九色精品人成在线观看| 一夜夜www| 国产精品1区2区在线观看.| 午夜福利在线观看吧| 亚洲欧美日韩东京热| 天天添夜夜摸| 一区二区三区国产精品乱码| 久久久久久大精品| 在线观看舔阴道视频| 在线观看美女被高潮喷水网站 | 国产又色又爽无遮挡免费看| 国产日本99.免费观看| 日本三级黄在线观看| 国产三级中文精品| 99热6这里只有精品| 国内揄拍国产精品人妻在线| 精品99又大又爽又粗少妇毛片 | 亚洲国产看品久久| 91麻豆av在线| 久久欧美精品欧美久久欧美| 日本撒尿小便嘘嘘汇集6| 人人妻,人人澡人人爽秒播| 国内精品久久久久久久电影| 精品久久久久久,| 黄频高清免费视频| 成年免费大片在线观看| 日本五十路高清| 国产日本99.免费观看| 亚洲美女黄片视频| 性欧美人与动物交配| 一级毛片高清免费大全| 最近视频中文字幕2019在线8| 91九色精品人成在线观看| 亚洲av电影在线进入| 欧美日韩瑟瑟在线播放| 美女黄网站色视频| 少妇的逼水好多| 亚洲七黄色美女视频| 麻豆成人av在线观看| 久9热在线精品视频| 免费看a级黄色片| 99久国产av精品| 久久这里只有精品中国| 国产又黄又爽又无遮挡在线| 国产黄a三级三级三级人| 一本综合久久免费| cao死你这个sao货| 国产精品亚洲一级av第二区| av女优亚洲男人天堂 | 一个人观看的视频www高清免费观看 | 国产成人精品久久二区二区91| 变态另类成人亚洲欧美熟女| 欧美成狂野欧美在线观看| 亚洲专区中文字幕在线| 欧美3d第一页| 在线观看一区二区三区| 欧美大码av| 国产三级黄色录像| 熟妇人妻久久中文字幕3abv| 久久亚洲真实| 色哟哟哟哟哟哟| 又爽又黄无遮挡网站| 国产私拍福利视频在线观看| 精品久久久久久久人妻蜜臀av| 免费观看人在逋| av天堂中文字幕网| 欧美色视频一区免费| 日本五十路高清| 欧美中文综合在线视频| 国产精品一区二区三区四区久久| 欧美高清成人免费视频www| 狠狠狠狠99中文字幕| 国产视频一区二区在线看| 一级作爱视频免费观看| 亚洲成人精品中文字幕电影| 色尼玛亚洲综合影院| 国产欧美日韩精品亚洲av| 欧美乱码精品一区二区三区| 夜夜爽天天搞| 一进一出抽搐gif免费好疼| 又爽又黄无遮挡网站| 国产91精品成人一区二区三区| 亚洲av中文字字幕乱码综合| 精品免费久久久久久久清纯| 国产精品亚洲一级av第二区| 亚洲成人久久爱视频| 精华霜和精华液先用哪个| 两个人看的免费小视频| 中文字幕熟女人妻在线| 老司机在亚洲福利影院|