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

    深孔鉆削刀具系統(tǒng)磁流變制振器設(shè)計(jì)及其動(dòng)態(tài)性能研究

    2016-08-04 06:16:01孔令飛孟維昌侯曉麗
    振動(dòng)與沖擊 2016年12期
    關(guān)鍵詞:磁極磁感應(yīng)阻尼

    孔令飛, 孟維昌, 侯曉麗, 王 杰, 李 超

    (西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院,西安 710048)

    深孔鉆削刀具系統(tǒng)磁流變制振器設(shè)計(jì)及其動(dòng)態(tài)性能研究

    孔令飛, 孟維昌, 侯曉麗, 王杰, 李超

    (西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院,西安710048)

    依據(jù)實(shí)際深孔鉆削刀具系統(tǒng)的結(jié)構(gòu)布局及其振動(dòng)特征,設(shè)計(jì)了一種新型刀具系統(tǒng)制振器的結(jié)構(gòu)布局形式。以制振器理想輸出阻尼力和內(nèi)磁感應(yīng)強(qiáng)度達(dá)到飽和為設(shè)計(jì)目標(biāo),給出了磁流變制振器結(jié)構(gòu)的優(yōu)化設(shè)計(jì)方法及其具體結(jié)構(gòu)參數(shù);運(yùn)用有限元法,分析了新型制振器的磁場分布特征及其輸出阻尼特征。以此為基礎(chǔ),通過深孔鉆削實(shí)驗(yàn)研究,獲得了不同勵(lì)磁電流和切削參數(shù)條件下刀具振動(dòng)的抑制效果,驗(yàn)證了所提出的新型磁流變制振構(gòu)型的可行性和有效性。

    深孔鉆削;刀具振動(dòng);制振器設(shè)計(jì);動(dòng)態(tài)特性

    近年來,隨著冶金、核電及兵器工業(yè)的發(fā)展和市場需求的多樣性,新型難加工材料的深孔加工、異型零件的深孔加工等越來越多,加工難度也越來越大[1]。實(shí)際鉆削中,由于深孔鉆削刀具處于封閉狀態(tài)、加工孔長徑比大(L/d>5)、刀具系統(tǒng)動(dòng)態(tài)剛性差等特點(diǎn),致使刀具系統(tǒng)的振動(dòng)行為極其復(fù)雜,易在工件表面形成“波浪”、“多角”、“過切”或“欠切”孔型。那么,能否有效、精準(zhǔn)地控制深孔加工刀具系統(tǒng)的動(dòng)態(tài)行為已成為人們?nèi)找骊P(guān)注的熱點(diǎn)問題[2-3]。

    為提高切削效率和加工品質(zhì),人們在鉆削加工刀具系統(tǒng)的振動(dòng)控制方面開展了大量研究工作。Mehrabadi等[4]構(gòu)建了考慮鉆削過程阻尼和質(zhì)量偏心影響的刀具系統(tǒng)動(dòng)力學(xué)模型,討論了刀具動(dòng)態(tài)運(yùn)行軌跡的特征及其穩(wěn)定性,并結(jié)合實(shí)驗(yàn)結(jié)果證實(shí)了該模型的準(zhǔn)確性。Roukema等[5-6]首次給出了包含有非線性特征、刀具幾何尺寸及全部振動(dòng)模式的鉆削過程耦合動(dòng)力學(xué)模型,通過對其進(jìn)行全面的時(shí)域仿真分析,獲得了鉆削過程刀具動(dòng)態(tài)特性的穩(wěn)定域。在此基礎(chǔ)之上,Ahmadi等[7]又提出了動(dòng)態(tài)鉆削過程的廣義穩(wěn)定性模型,模型中考慮了刀具渦動(dòng)和扭轉(zhuǎn)振動(dòng)對切厚再生效應(yīng)的影響,實(shí)現(xiàn)了在半離散時(shí)域內(nèi)鉆削刀具的穩(wěn)定性預(yù)測??傮w來說,上述這些方法可理解為從全耦合動(dòng)力學(xué)機(jī)理建模的角度來預(yù)測與控制刀具的動(dòng)態(tài)行為,它們均是通過調(diào)整切削轉(zhuǎn)速或進(jìn)給量來避開不穩(wěn)定切削區(qū)域,進(jìn)而避免再生振動(dòng)或自激振動(dòng)的出現(xiàn)??墒牵?dāng)由刀具-工件所組成的切削系統(tǒng)中有一者改變時(shí),其穩(wěn)定域也會(huì)隨之發(fā)生變化,為實(shí)際應(yīng)用帶來了困難。

    然而,上述研究給了我們一個(gè)重要的啟示:在實(shí)際鉆削過程中,若能夠抑制某些振動(dòng)模態(tài)使其不被激發(fā)出來,或降低其振動(dòng)的能量,則該振動(dòng)模態(tài)就不會(huì)對鉆削過程加工孔品質(zhì)產(chǎn)生太大影響,進(jìn)而有效提高加工孔品質(zhì)?;诖耍罁?jù)深孔鉆削刀具系統(tǒng)的實(shí)際工作特點(diǎn),設(shè)計(jì)了一種新型刀具系統(tǒng)制振器的結(jié)構(gòu)布局形式。以制振器理想輸出阻尼力和內(nèi)磁感應(yīng)強(qiáng)度達(dá)到飽和為設(shè)計(jì)目標(biāo),給出了磁流變制振器結(jié)構(gòu)的優(yōu)化設(shè)計(jì)方法及其具體結(jié)構(gòu)參數(shù),并結(jié)合有限元法,分析了新型制振器的磁場分布特征及其輸出阻尼特征。以此為基礎(chǔ),通過鉆削實(shí)驗(yàn)研究,獲得了勵(lì)磁電流、制振位置及進(jìn)給量對制振器抑振效果的影響,同時(shí)也驗(yàn)證了本文所提出的新型磁流變制振構(gòu)型的可行性和有效性。

    1深孔加工刀具振動(dòng)抑制構(gòu)型的結(jié)構(gòu)及其工作原理

    本文所提出的制振器結(jié)構(gòu)設(shè)計(jì)方案見圖1,該結(jié)構(gòu)主要由磁流變液體、殼體、活塞、活塞盤、磁極和線圈等部分組成。實(shí)際鉆削中,將鉆桿穿過導(dǎo)套,通過調(diào)整施加于制振器環(huán)形線圈的勵(lì)磁電流大小,即可改變活塞盤和磁流變液體接觸面的阻尼值和剛度特征。該結(jié)構(gòu)的最大特點(diǎn)是通過調(diào)整線圈電流數(shù)值,以期對某些有害于加工精度和表面質(zhì)量的振動(dòng)模態(tài)實(shí)現(xiàn)摧毀,或抑制其不被激發(fā)出來,最終鉆削形成預(yù)定的零件孔品質(zhì),具體深孔鉆削系統(tǒng)布局形式見圖2。

    圖1 磁流變制振器結(jié)構(gòu)示意圖Fig.1 The structure schematic diagram of magnetorheological vibration suppression device

    圖2 深孔鉆削加工原理示意圖Fig.2 The schematic diagram of deep-hole drilling process

    2制振器的設(shè)計(jì)方法

    2.1制振器設(shè)計(jì)的基本原則

    在磁流變液制振器的設(shè)計(jì)中,制振器輸出的振動(dòng)控制力及其響應(yīng)速度不僅與磁流變效應(yīng)的強(qiáng)弱及速度有關(guān),并與制振器各部件尺寸和阻尼間隙也密切相關(guān),而這些因素歸根結(jié)底又與如何構(gòu)建磁流變液黏性模型存在著緊密聯(lián)系。本文基于Bingham黏性模型[8-9],若制振器線圈未通電時(shí),制振器腔體內(nèi)的磁流變液表現(xiàn)為牛頓流體特性,其剪應(yīng)力與剪應(yīng)變成正比,而當(dāng)施加勵(lì)磁電流后,環(huán)形線圈產(chǎn)生外加磁場,相應(yīng)的阻尼通道內(nèi)磁流變液體的剪切力函數(shù)可描述為如下形式:

    (1)

    式中:τ為總剪切應(yīng)力;τy(field)與磁場大小相關(guān)的剪切應(yīng)力;η為磁流變液黏性系數(shù)(即零磁場黏度),由磁流變液性質(zhì)及剪切應(yīng)變率大小決定;γ為剪切應(yīng)變率。

    據(jù)此,制振器輸出的振動(dòng)控制力由可調(diào)庫倫阻尼力Fτ和黏性阻尼力Fη組成,其力學(xué)關(guān)系表達(dá)式為:

    F=Fτ+Fη

    (2)

    其中

    (3)

    式中:S為磁極截面積;v為活塞盤移動(dòng)速;g為阻尼通道間隙,其計(jì)算方法為:

    (4)

    式中:K為控制器可調(diào)系數(shù),即最大可調(diào)庫倫阻尼力Fτ與最大黏性阻尼力Fη的比值。

    對比式(1)、式(2)和式(4)可知,制振器產(chǎn)生的剪切應(yīng)力τy隨著間隙距離g的增加明顯下降,這主要是由于線圈在阻尼通道內(nèi)產(chǎn)生的磁感應(yīng)強(qiáng)度隨著間隙g的增加明顯下降,因而為了增大磁場強(qiáng)度則需要減小間隙g。但是,g的減小會(huì)使黏性阻尼力Fη增大,致使可調(diào)系數(shù)K變小,進(jìn)而制振力的實(shí)際輸出調(diào)節(jié)范圍變小。因此,實(shí)際設(shè)計(jì)中,制振器的基本設(shè)計(jì)原則就是需要合理設(shè)計(jì)間隙g,使得制振器輸出的振動(dòng)控制力及其調(diào)節(jié)范圍均可滿足實(shí)際使用的需求。一般說來,間隙g取值范圍為0.8~2 mm較為合理[10]。

    2.2制振器的材料選擇

    依據(jù)深孔加工制振器的實(shí)際工作特點(diǎn),并綜合考慮材料的磁學(xué)性能、結(jié)構(gòu)強(qiáng)度以及性價(jià)比等因素,新型制振器的殼體和活塞盤均選擇20#低碳鋼,而制振器的繞線盤和活塞則選擇鋁材料,這主要是基于鋁質(zhì)材料的不導(dǎo)磁性(相對磁導(dǎo)率為1),可使磁場集中在制振器的工作回路中,從而減少漏磁,提高磁場利用效率。對于磁流變液的選擇,本文選用MRF-132DG磁流變液,該材料為碳?xì)浠帕饕后w,它的剪切應(yīng)力與剪切應(yīng)變率關(guān)系、黏性系數(shù)與剪切應(yīng)變率的關(guān)系見圖3[11],屈服應(yīng)力τ和磁場強(qiáng)度H關(guān)系曲線見圖4[12]。

    圖3 MRF- 132DG剪切應(yīng)變率性質(zhì)Fig.3 Shear strain rate properties of MRF- 132DG

    圖4 MRF- 132DG屈服應(yīng)力與磁場強(qiáng)度關(guān)系Fig.4 The relationship between yield stress and magnetic field intensity of MRF- 132DG

    2.3磁路系統(tǒng)的設(shè)計(jì)

    在磁路系統(tǒng)設(shè)計(jì)中,間隙通道內(nèi)磁流變液的磁感應(yīng)強(qiáng)度B是關(guān)鍵參數(shù)。首先,根據(jù)磁流變液B-H曲線和磁特性曲線τ-H曲線得到磁流變液材料的工作點(diǎn)。本文選擇磁流變液達(dá)到磁飽和時(shí)(屈服應(yīng)力τy最大),對應(yīng)工作點(diǎn)的磁場強(qiáng)度Hf=250 kA/m,磁感應(yīng)強(qiáng)度Bf=0.9 T,磁流變液B-H曲線見圖5[12]。

    圖5 MRF-132DG B-H曲線Fig.5 B-H curve of MRF-132DG

    為了便于勵(lì)磁線圈匝數(shù)的計(jì)算,可將磁流變制振器的結(jié)構(gòu)參數(shù)作如下定義(見圖6):g為阻尼間隙;h為活塞盤厚度;e磁極厚度;L1側(cè)邊壁厚度;L2側(cè)邊壁長度;R1和R2分別為磁極內(nèi)半徑和外半徑;R3和R4分別為殼體內(nèi)半徑和外半徑。

    圖6 磁流變制振器的結(jié)構(gòu)參數(shù)Fig.6 The design of magnetorheological vibration suppression device structure

    由安培環(huán)路定理,可得整個(gè)回路的磁動(dòng)勢,

    ε=NI

    (5)

    式中:N為線圈匝數(shù);I為線圈內(nèi)流通的電流。

    同時(shí),相應(yīng)的磁路中所需總磁通量:

    Φ=BfSf

    (6)

    式中:Bf為磁流變液飽和磁感應(yīng)強(qiáng)度;Sf為阻尼間隙有效截面積。

    在不考慮漏磁的情況下,根據(jù)基爾霍夫第一定律(磁通守恒定律),磁路中各組成分磁通量關(guān)系為:

    Φ=BfSf=BASA=BESE=…

    (7)

    式中:BA、BE等為分布在磁路中各部分的磁感應(yīng)強(qiáng)度,SA、SB等為各部分的有效截面積。由于磁回路各部分材料均為20#低碳鋼,且磁感應(yīng)強(qiáng)度B呈現(xiàn)非線性變化,故相應(yīng)磁感應(yīng)強(qiáng)度Bi(Bi=BA、BE,…)可描述為:

    (8)

    由基爾霍夫第二定律和安培環(huán)路定理,可得:

    (9)

    式中:Hi為各磁路部分磁場強(qiáng)度,結(jié)合式(5)~式(7)可得:

    (10)

    式中:μf為磁流變液磁導(dǎo)率,μi為磁路各部分材料磁導(dǎo)率。當(dāng)額定工作電流Im確定后,結(jié)合圖6制振器各部分尺寸,可由下式計(jì)算所需線圈匝數(shù):

    NI=Hf2g+H1h+H22(e+L1/2)+

    H32[L2-(R2-R1)/2]+H4(L3-L1)

    (11)

    式中:H1、H2、H3和H4可由式(10)計(jì)算得到。

    3制振器結(jié)構(gòu)的優(yōu)化設(shè)計(jì)

    為了實(shí)際鉆削時(shí)新型制振器安裝便利,以及獲得足夠大的制振力輸出范圍和結(jié)構(gòu)強(qiáng)度,可先利用式(1)~式(4)進(jìn)行初步計(jì)算,得到制振器的初步設(shè)計(jì)結(jié)構(gòu)參數(shù),見表1。然后,依據(jù)制振器設(shè)計(jì)原則,將制振器具有較大的可調(diào)系數(shù)K及阻尼通道內(nèi)磁感應(yīng)強(qiáng)度達(dá)到0.9 T作為優(yōu)化設(shè)計(jì)目標(biāo),對初步結(jié)構(gòu)設(shè)計(jì)參數(shù)進(jìn)行優(yōu)化,以期實(shí)現(xiàn)磁流變液使用效率最高,且降低制振器整體體積和質(zhì)量。制振器的結(jié)構(gòu)優(yōu)化可分為兩個(gè)部分,即徑向結(jié)構(gòu)優(yōu)化和軸向結(jié)構(gòu)優(yōu)化。

    表1 磁流變制振器初步設(shè)計(jì)結(jié)構(gòu)參數(shù)

    3.1制振器徑向結(jié)構(gòu)優(yōu)化

    制振器徑向結(jié)構(gòu)的優(yōu)化應(yīng)根據(jù)制振器理想輸出控制力F,并結(jié)合制振控制力計(jì)算式(2),即可求解得到徑向結(jié)構(gòu)尺寸,其中活塞盤和磁極截面積是影響制振控制力大小和間隙通道內(nèi)磁流變液能否達(dá)到磁飽和的關(guān)鍵。另外,由于黏性阻尼力Fη比可調(diào)庫倫阻尼力Fτ小很多,且在低速區(qū)時(shí)Fη與鉆桿擾動(dòng)速度有關(guān),因而制振控制力F幾乎完全由Fτ提供。若可調(diào)庫倫阻尼力Fτ為最大值時(shí),則制振器的輸出控制力F即為制振器的理想目標(biāo)值。基于此,當(dāng)制振器輸出控制力F、磁流變液飽和剪切應(yīng)力τy確定之后,依據(jù)式(2)和式(3),就可獲得磁極截面積SE,具體計(jì)算公式為:

    (12a)

    結(jié)合磁流變制振器的結(jié)構(gòu)參數(shù)圖6,磁極截面積SE又可以寫為:

    (12b)

    式中:R1和R2分別為磁極內(nèi)半徑和外半徑。

    在制振器的徑向尺寸設(shè)計(jì)過程中,應(yīng)遵循在滿足強(qiáng)度要求的條件下盡量減小制振器體積和質(zhì)量。本文所采用鉆桿直徑為17 mm,鉆桿套厚度2~4 mm,軸承選取深溝球軸承。因而,考慮到磁極內(nèi)半徑R1受深孔鉆桿直徑影響,依據(jù)式(12),可計(jì)算獲得內(nèi)半徑R1=36.5 mm,外半徑R2=49 mm,理想輸出阻尼力Fn=300 N,τymax=44.7 kPa。

    3.2制振器軸向結(jié)構(gòu)優(yōu)化

    在制振器軸向結(jié)構(gòu)設(shè)計(jì)中,相關(guān)設(shè)計(jì)變量和約束條件如下:

    3.2.1設(shè)計(jì)變量

    軸向結(jié)構(gòu)優(yōu)化設(shè)計(jì)中,有5個(gè)設(shè)計(jì)變量,分別為阻尼間隙g、活塞盤厚度h、磁極厚度e、側(cè)邊壁厚度L1和側(cè)邊壁長度L2。

    3.2.2約束條件

    ① 非線性不等式約束

    根據(jù)設(shè)計(jì)要求, 可調(diào)庫倫阻尼力Fτ的最大值需達(dá)到制振控制力F,即

    Fτ=2Sτy≥Fn

    (13)

    ② 非線性等式約束

    磁感應(yīng)強(qiáng)度的計(jì)算必須滿足基爾霍夫第一定律,即

    Φ=BfSf=BESE

    (14)

    式中:Bf為間隙處磁感應(yīng)強(qiáng)度;Sf為阻尼間隙截面積;BE為磁極處磁感應(yīng)強(qiáng)度;SE為磁極截面積;

    ③ 線性不等式約束

    實(shí)際結(jié)構(gòu)設(shè)計(jì)中,磁極內(nèi)半徑R1應(yīng)小于磁極外半徑R2,且阻尼間隙g則應(yīng)小于活塞盤厚度h,即:

    (15)

    g-h<0

    (16)

    ④ 線性等式約束

    制振器的總厚度L3應(yīng)該等于阻尼間隙、活塞盤厚度、磁極厚度和側(cè)邊壁厚度之和:

    L3=h+2g+2e+2L1

    (17)

    3.2.3優(yōu)化目標(biāo)

    由于可調(diào)系數(shù)K反映了制振器輸出控制力的可調(diào)節(jié)能力,為增大阻尼器可調(diào)范圍,在優(yōu)化設(shè)計(jì)中,以獲得最大可調(diào)系數(shù)為優(yōu)化目標(biāo),即:

    (18)

    3.2.4優(yōu)化求解

    約束條件確定后,借助Matlab中fmincon函數(shù)來實(shí)現(xiàn)優(yōu)化計(jì)算,分別取理想控制力Fn=300 N及最大擾動(dòng)速度Vmax=3.5 m/s進(jìn)行計(jì)算,并將優(yōu)化計(jì)算結(jié)果取整,相應(yīng)優(yōu)化后制振器的結(jié)構(gòu)參數(shù)見表2。此外,考慮到勵(lì)磁線圈所占面積將會(huì)影響側(cè)邊壁厚度L1和側(cè)邊壁長度L2的大小,因而實(shí)際設(shè)計(jì)時(shí)L1和L2兩個(gè)參數(shù)還要根據(jù)線圈線徑和匝數(shù)做出具體調(diào)整。

    表2 磁流變制振器優(yōu)化后結(jié)構(gòu)參數(shù)

    由表2可知,相對于表1的初步結(jié)構(gòu)參數(shù)而言,優(yōu)化后的活塞盤厚度、磁極厚度和阻尼間隙并未改變,僅是側(cè)邊壁厚度有所增加,而側(cè)邊壁長度、磁極內(nèi)半徑和外半徑尺寸則明顯降低,因而減少了制振器的整體體積和質(zhì)量。

    4制振器磁路系統(tǒng)的數(shù)值分析

    為了驗(yàn)證新型制振器的結(jié)構(gòu)參數(shù)是否滿足實(shí)際需求,運(yùn)用ANSYS 有限元軟件中的PLANE53單元對制振器阻尼間隙、磁極和殼體回路部分的磁場分布情況進(jìn)行分析與討論。同時(shí)考慮到新型磁流變制振器是一個(gè)軸對稱結(jié)構(gòu)(見圖7所示),因而取對稱平面模型的一半分析即可,其網(wǎng)格模型見圖8。

    圖7 磁流變制振器幾何模型Fig.7Geometricmodelofmagnetorheologicalvibrationsuppressiondevice圖8 磁流變制振器網(wǎng)格劃分的模型Fig.8Themeshingmodelofmagnetorheologicalvibrationsuppressiondevice

    取表1所示的制振器初步設(shè)計(jì)尺寸進(jìn)行計(jì)算分析,求解獲得的磁感應(yīng)強(qiáng)度、磁場強(qiáng)度分布見圖9。由圖9(a)或圖9(b)可知,磁場分布集中于制振器回路結(jié)構(gòu)中,殼體處磁感應(yīng)強(qiáng)度最大,約為1.6 T,殼體材料已經(jīng)達(dá)到飽和值;阻尼間隙內(nèi)約為0.2 ~0.6 T,故磁流變液未達(dá)到飽和值,且分布不均勻。此外,從圖9(c)和圖9(d)中還可知,殼體處磁場強(qiáng)度約為3 100 A/m,間隙內(nèi)磁場強(qiáng)度約為79 000 A/m,阻尼通道內(nèi)磁場強(qiáng)度方向都垂直于活塞盤的運(yùn)動(dòng)方向,但是隨著與勵(lì)磁線圈之間距離的增加,通道內(nèi)磁場強(qiáng)度明顯減小。

    圖10為間隙通道內(nèi)磁感應(yīng)強(qiáng)度分布曲線,其中L為以制振器活塞盤頂端為起點(diǎn),豎直方向的阻尼通道長度。從圖10可知,磁力線回路主要集中于磁極處,磁力線走向符合設(shè)計(jì)預(yù)想,驗(yàn)證了制振器設(shè)計(jì)原理的正確性。但是,間隙內(nèi)磁感應(yīng)強(qiáng)度最大值約為0.62 T,未達(dá)到磁流變液飽和值。隨著阻尼通道長度L的增加,磁感應(yīng)強(qiáng)度不斷減小,且分布不均勻。這說明初步設(shè)計(jì)的制振器結(jié)構(gòu)尺寸存在冗余,制振器性能并未得到充分發(fā)揮。

    圖9 初步設(shè)計(jì)條件下磁流變制振器有限元分析結(jié)果Fig.9 The finite element analysis results of preliminary design magnetorheological vibration suppression device

    依據(jù)表2所示優(yōu)化后的制振器設(shè)計(jì)尺寸進(jìn)行有限元分析,求解獲得的磁感應(yīng)強(qiáng)度、磁場強(qiáng)度分布見圖11。由圖11(a)或圖11(b)可知,殼體處磁感應(yīng)強(qiáng)度約1.65 T,達(dá)到殼體材料的飽和值,殼體頂壁處磁感應(yīng)強(qiáng)度略大,約為1.7 T,說明頂壁尺寸略小,后續(xù)設(shè)計(jì)中,頂壁尺寸需要調(diào)整。阻尼通道內(nèi)磁感應(yīng)強(qiáng)度約為0.82 T,接近預(yù)想值0.9 T,且磁場強(qiáng)度分布均勻。此外,從圖11(c)和圖11(d)可知,殼體磁場強(qiáng)度約為3 200 A/m,間隙內(nèi)磁場強(qiáng)度最大,約為210 000 A/m,這說明阻尼通道內(nèi)磁場強(qiáng)度方向都垂直于活塞盤的運(yùn)動(dòng)方向,使得磁場、磁極和活塞盤面積得到有效利用。

    圖10 初始設(shè)計(jì)條件下磁流變制振器磁感應(yīng)強(qiáng)度分析Fig.10 The analysis of the magnetic intensity of preliminary design magnetorheological vibration suppression device

    圖11 制振器優(yōu)化設(shè)計(jì)磁流變有限元分析結(jié)果Fig.11 The finite element analysis results of optimization designmagnetorheological vibration suppression device

    圖12給出了間隙通道內(nèi)磁感應(yīng)強(qiáng)度的分布曲線,從圖中可見間隙內(nèi)磁感應(yīng)強(qiáng)度約為0.82 T。隨著阻尼通道長度L的增加,間隙內(nèi)磁感應(yīng)強(qiáng)度變化不大,分布均勻。當(dāng)距離L大于間隙通道長度后,磁感應(yīng)強(qiáng)度降至最低。這說明優(yōu)化后的制振器結(jié)構(gòu)尺寸達(dá)到了設(shè)計(jì)要求,磁場分布主要集中于磁極處,磁極作用明顯,整個(gè)制振器磁場回路都經(jīng)由磁極處通過,磁極外漏磁較少,磁場利用率明顯提高。

    圖12 優(yōu)化設(shè)計(jì)阻尼間隙長度方向磁感應(yīng)強(qiáng)度B分布曲線Fig.12 The analysis of magnetic intensity of optimization design magnetorheological vibration suppression device

    5新型制振器的抑振效果實(shí)驗(yàn)

    本節(jié)將通過實(shí)驗(yàn)研究,獲得優(yōu)化設(shè)計(jì)后的制振器構(gòu)型對刀具振動(dòng)的抑制效果,以期驗(yàn)證本文所提出的新型磁流變液制振器設(shè)計(jì)方案的可行性和有效性。

    5.1深孔鉆削抑振實(shí)驗(yàn)平臺

    本次刀具振動(dòng)抑制實(shí)驗(yàn)是在本所自制的工件回轉(zhuǎn)型深孔鉆削機(jī)床上進(jìn)行,相關(guān)實(shí)驗(yàn)平臺見圖13。實(shí)驗(yàn)過程中,利用刀具輔助支撐將制振器固定于機(jī)床導(dǎo)軌上,并通過兩個(gè)非接觸式位移傳感器分別采集實(shí)際鉆削時(shí)刀具在垂直和水平方向上的振動(dòng)信號。加工孔的圓度誤差用CONTURA-G2型三坐標(biāo)測量機(jī)來獲得。相關(guān)加工參數(shù)如下:刀桿外徑Φ17 mm,桿長1 100 mm,刀頭為錯(cuò)齒刀具,外徑Φ19 mm。材料為高強(qiáng)度合金鋼,密度為7.87×103kg/m3;供油壓力pin=2×106Pa,切削液動(dòng)力黏度μ=0.026 Pa·s;工件材質(zhì)為45#鋼,長度為300 mm。

    5.2實(shí)驗(yàn)結(jié)果分析

    5.2.1勵(lì)磁電流對抑振效果的影響

    實(shí)驗(yàn)過程中,將制振器固定于距離授油器550 mm的位置處,設(shè)定工件轉(zhuǎn)速和進(jìn)給量分別為710 r/min和0.084 0 mm/r,取0.0 A、0.2 A、0.5 A三種不同的勵(lì)磁電流進(jìn)行鉆削實(shí)驗(yàn),獲得的刀具振動(dòng)位移的時(shí)間歷程響應(yīng)見圖14。

    圖14 不同勵(lì)磁電流條件下刀具振動(dòng)的時(shí)間歷程響應(yīng)Fig.14 The vibration responses of drill tool system under different currents

    由圖14可知,當(dāng)制振器位置和進(jìn)給量不變時(shí),隨著勵(lì)磁電流強(qiáng)度的增大,刀具振動(dòng)幅值明顯減小。這說明隨著制振器勵(lì)磁線圈電流的改變,引起刀具系統(tǒng)在制振器位置處的剛度和阻尼特性也發(fā)生了變化,從而使刀具的振動(dòng)幅值得到有效控制。相應(yīng)加工條件下測量獲得的加工孔圓度誤差見表3,隨著施加的勵(lì)磁電流增大,加工孔的圓度誤差得到了改善,且均優(yōu)于施加電流為0.0 A時(shí),僅是改善的程度有所區(qū)別而已。

    表3 不同勵(lì)磁電流條件下獲得的加工孔圓度誤差

    5.2.2制振位置對抑振效果的影響

    實(shí)驗(yàn)中設(shè)定主軸轉(zhuǎn)速為710 r/min,進(jìn)給量為0.084 0 mm/r,勵(lì)磁線圈電流為0.3 A,制振器分別固定于位置1(距授油器400 mm)、位置2(距授油器550 mm)和位置3(距授油器750 mm)處,測得實(shí)際加工時(shí)刀具的振動(dòng)信號見圖15。

    從圖15可知,當(dāng)制振器從位置3變換到位置1時(shí),隨著與刀頭位置間距離的減小,刀具振動(dòng)幅值也呈減小趨勢。此外,當(dāng)制振器放置于位置3時(shí),刀具振動(dòng)幅值呈現(xiàn)出較大的波動(dòng),且振動(dòng)不平穩(wěn),這說明該位置對刀具振動(dòng)的控制能力較差。與此同時(shí),由表3可知,當(dāng)制振器放置于位置1和位置2時(shí),相應(yīng)加工孔圓度誤差值也都優(yōu)于放置于位置3處,因而實(shí)際鉆削中應(yīng)避免選用位置3進(jìn)行振動(dòng)抑制。

    圖15 不同制振位置條件下刀具振動(dòng)的時(shí)間歷程響應(yīng)Fig.15 The vibration responses of drill tool system under different suppression positions

    支撐位置/mm400550750圓度誤差/mm0.002420.002350.00562

    5.2.3進(jìn)給量對抑振效果的影響

    圖16 不同進(jìn)給量下鉆桿振動(dòng)位移曲線Fig.16 The vibration responses of drill tool system under different feed rates

    鉆削實(shí)驗(yàn)中,設(shè)定主軸轉(zhuǎn)速為710 r/min,勵(lì)磁線圈電流為0.5 A,分別選取進(jìn)給量為0.076 5 mm/r、0.084 0 mm/r和0.091 5 mm/r,測取信號繪制刀具振動(dòng)的時(shí)間歷程響應(yīng)見圖16。從圖16可知,隨著進(jìn)給量的增加,刀具振動(dòng)的幅值均明顯增大。究其原因主要是由于隨著進(jìn)給量的增加,刀具系統(tǒng)所承受的切削力不斷增大,致使刀具的振動(dòng)幅值也變大。同時(shí),刀具振動(dòng)幅值的增大,最終導(dǎo)致被加工孔圓度誤差的增加(見表5)。

    表5 不同進(jìn)給量條件下獲得的加工孔圓度誤差

    6結(jié)論

    針對實(shí)際深孔鉆削過程中刀具系統(tǒng)易發(fā)生再生振動(dòng)或自激振動(dòng)的問題,設(shè)計(jì)了一種基于磁流變液材料的新型振動(dòng)抑制構(gòu)型。依據(jù)深孔刀具的實(shí)際工作特點(diǎn),并結(jié)合有限元法,以制振器理想輸出阻尼力和阻尼間隙內(nèi)磁感應(yīng)強(qiáng)度達(dá)到磁飽和為設(shè)計(jì)目標(biāo),給出了深孔鉆削刀具制振器的優(yōu)化設(shè)計(jì)方法及相關(guān)具體結(jié)構(gòu)參數(shù)。通過實(shí)際鉆削試驗(yàn),討論了勵(lì)磁電流、進(jìn)給量及制振位置對刀具振動(dòng)抑制效果的影響。從實(shí)驗(yàn)結(jié)果可以看出,隨著勵(lì)磁線圈電流的增大,刀具振動(dòng)幅值明顯減??;當(dāng)選用不同的進(jìn)給量和制振位置時(shí),新型制振器也對刀具振動(dòng)起到了有效抑制,僅是抑制程度所不同。由此可以看出,本文所提出的新型磁流變制振構(gòu)型及其優(yōu)化設(shè)計(jì)方法,將為實(shí)現(xiàn)深孔鉆削過程中刀具系統(tǒng)動(dòng)態(tài)特征及加工孔品質(zhì)的實(shí)時(shí)控制奠定理論與應(yīng)用基礎(chǔ)。

    [1] 李言,孔令飛.振動(dòng)切削深孔加工初始偏差對孔直線度誤差的影響[J].機(jī)械工程學(xué)報(bào),2012,48(13):167-172.

    LI Yan, KONG Ling-fei. Influence of initial deviation in vibrating deep hole drilling on hole straightness error [J]. Journal of Mechanical Engineering,2012,48(13):167-173.

    [2] Yao Z H, Mei D Q, Chen Z C. Chatter suppression by parametric excitation: model and experiments[J]. Journal of Sound and Vibration, 2011, 330: 2995-3005.

    [3] Quintana G, Ciurana J. Chatter in machining processes: a review[J]. International Journal of Machine Tools and Manufacture, 2011, 51(5): 363-376.

    [4] Mehrabadi I M, Nouri M, Madoliat R. Investigating chatter vibration in deep drilling, including process damping and the gyroscopic effect [J]. International Journal of Machine Tools and Manufacture, 2009, 49: 939-946.

    [5] Roukema J C, Altintas Y. Generalized modeling of drilling vibrations. Part I: Time domain model of drilling kinematics, dynamics and hole formation [J]. International Journal of Machine Tools and Manufacture, 2007, 47: 1455-1473.

    [6] Roukema J C, Altintas Y. Generalized modeling of drilling vibrations. Part II: Chatter stability in frequency domain [J]. International Journal of Machine Tools and Manufacture, 2007, 47: 1474-1485.

    [7] Ahmadi K, Altintas Y. Stability of lateral, torsional and axial vibrations in drilling [J]. International Journal of Machine Tools and Manufacture, 2013, 68: 63-74.

    [8] Forte P, Paterno M, Rustighi E. A magnetorheological fluid damper for rotor applications[J]. International Journal of Rotating Machinery, 2004, 10(3):175-182.

    [9] Ghasemi A H, Ohadi A R, Ghaffari Saadat M H. Vibration control of a rotor-bearing system with smart bearings unsing magnetorheologic fluid [C]//15th International Congress on Sound and Vibration. Daejeon, Korea:ICSV15, 2008.

    [10] Huang J, Zhang J Q, Yang Y, et al. Analysis and design of a cylindrical magneto-rheological fluid brake [J]. Journal of Materials Processing Technology, 2002, 129:559-562.

    [11] 王修勇,孫洪鑫,陳政清. 旋轉(zhuǎn)剪切式磁流變液阻尼器設(shè)計(jì)及力學(xué)模型[J]. 振動(dòng)與沖擊,2010,29(10):77-81.

    WANG Xiu-yong, SUN Hong-xin, CHEN Zheng-qing. Design and mechanical model of rotary-mode shear MR damper [J]. Journal of Vibration and Shock,2010,29(10): 77-81.

    [12] Lord Company. MRF-132DG magneto-rheological fluid [R]. Lord Technical Data, 2008.

    Design and dynamic characteristics of magnetorheological vibration suppression device in the deep-hole drilling process

    KONG Ling-fei, MENG Wei-chang, HOU Xiao-li, WANG Jie, LI Chao

    (School of Mechanical and Precision Instrument Engineering, Xi’an University of Technology, Xi’an 710048, China)

    Considering the structure and vibration characteristics of a drilling-tools system in deep-hole drilling, a novel structural layout of the vibration suppression device is proposed. With the goal of obtaining a desired damping force output and controlling the magnetic saturation strength, the structure optimization of a novel magnetorheological vibration suppression device and its specific structure parameters are given. Combining this with the finite element method, the magnetic field distribution of the proposed vibration suppression device and its damping characteristics output are analyzed. Based on the processes mentioned above, the vibration control experiments of the tool system in deep-hole drilling is carried out so as to investigate the relevant experimental results of tool vibration control under the conditions of different excitation currents and cutting parameters. It can be seen from the experimental results that the new model of magnetorheological damping configuration proposed in this paper satisfies the demands of the practical deep-hole machining process.

    deep-hole drilling; drilling tool vibration; vibration suppression device design; dynamic characteristics

    10.13465/j.cnki.jvs.2016.12.018

    國家自然科學(xué)基金(51475367);國家科技重大專項(xiàng)(2014ZX04001-191)資助項(xiàng)目

    2015-05-20修改稿收到日期:2015-11-13

    孔令飛 男,副教授,1977年生

    TH113.1

    A

    猜你喜歡
    磁極磁感應(yīng)阻尼
    跨空海界面磁感應(yīng)通信特性分析及應(yīng)用
    同步電機(jī)轉(zhuǎn)子磁極結(jié)構(gòu)
    電磁感應(yīng)中的“知三求三”
    固定同步電機(jī)磁極用螺栓的受力分析
    淺析芳綸紙?jiān)谒啺l(fā)電機(jī)磁極中的應(yīng)用
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項(xiàng)的擴(kuò)散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    寧波漢浦工具有限公司
    具阻尼項(xiàng)的Boussinesq型方程的長時(shí)間行為
    99热6这里只有精品| 国产又色又爽无遮挡免| 成年人午夜在线观看视频| 一区在线观看完整版| 国产成人精品福利久久| 18禁观看日本| 久久国产亚洲av麻豆专区| 欧美精品国产亚洲| 国产精品偷伦视频观看了| 在线精品无人区一区二区三| 国产精品久久久av美女十八| 亚洲欧洲日产国产| 26uuu在线亚洲综合色| 国产av码专区亚洲av| 久热这里只有精品99| 久久人人97超碰香蕉20202| 国产精品嫩草影院av在线观看| 亚洲,一卡二卡三卡| 99香蕉大伊视频| 色网站视频免费| 亚洲,欧美精品.| 高清毛片免费看| 久久精品久久久久久久性| 日韩欧美精品免费久久| 久久久精品94久久精品| 国产黄色免费在线视频| 全区人妻精品视频| 99国产精品免费福利视频| 亚洲三级黄色毛片| 国产欧美日韩一区二区三区在线| 在线观看免费高清a一片| 国产日韩一区二区三区精品不卡| 黑人高潮一二区| 如日韩欧美国产精品一区二区三区| 亚洲av日韩在线播放| kizo精华| 韩国精品一区二区三区 | 春色校园在线视频观看| 午夜激情久久久久久久| 亚洲欧美成人综合另类久久久| 亚洲精品aⅴ在线观看| 免费高清在线观看日韩| 午夜免费观看性视频| 日本91视频免费播放| 欧美 日韩 精品 国产| 美国免费a级毛片| 啦啦啦中文免费视频观看日本| 男女高潮啪啪啪动态图| 日韩三级伦理在线观看| 草草在线视频免费看| 永久免费av网站大全| 91精品三级在线观看| 精品一品国产午夜福利视频| 国产乱人偷精品视频| 亚洲欧美日韩另类电影网站| 国产在视频线精品| 高清av免费在线| 18禁裸乳无遮挡动漫免费视频| 草草在线视频免费看| 婷婷色综合大香蕉| videosex国产| 人妻一区二区av| 国产又爽黄色视频| 男女午夜视频在线观看 | 日韩熟女老妇一区二区性免费视频| 国产免费又黄又爽又色| 黄色 视频免费看| 大片免费播放器 马上看| 日日爽夜夜爽网站| 成人18禁高潮啪啪吃奶动态图| 亚洲精品色激情综合| 成人国产麻豆网| 大香蕉久久成人网| 人人妻人人澡人人爽人人夜夜| 色5月婷婷丁香| 精品一区二区三卡| 乱码一卡2卡4卡精品| 国产成人欧美| 人妻少妇偷人精品九色| 一本大道久久a久久精品| 久久韩国三级中文字幕| 黄片播放在线免费| 2018国产大陆天天弄谢| 久久久久久久国产电影| 久久久精品区二区三区| 少妇 在线观看| 中文乱码字字幕精品一区二区三区| 久久午夜福利片| 国精品久久久久久国模美| 亚洲内射少妇av| xxx大片免费视频| 777米奇影视久久| 中文字幕人妻熟女乱码| 午夜av观看不卡| 欧美日韩综合久久久久久| 精品一区在线观看国产| 欧美精品国产亚洲| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产自在天天线| 精品国产一区二区久久| 精品国产一区二区三区久久久樱花| 免费在线观看完整版高清| 男女午夜视频在线观看 | 大香蕉久久网| 一区二区日韩欧美中文字幕 | av在线老鸭窝| 国产爽快片一区二区三区| 婷婷成人精品国产| 一级黄片播放器| 亚洲人与动物交配视频| 亚洲 欧美一区二区三区| www.熟女人妻精品国产 | 成人影院久久| 国产视频首页在线观看| 久久久国产一区二区| 高清av免费在线| 菩萨蛮人人尽说江南好唐韦庄| 国产熟女午夜一区二区三区| 国产欧美日韩综合在线一区二区| 伊人亚洲综合成人网| 国产成人精品久久久久久| 人人妻人人澡人人爽人人夜夜| 国产白丝娇喘喷水9色精品| 精品人妻熟女毛片av久久网站| 97超碰精品成人国产| 亚洲精品一二三| 高清视频免费观看一区二区| 色5月婷婷丁香| 在线天堂中文资源库| 日韩成人伦理影院| 日韩视频在线欧美| 久久国产精品男人的天堂亚洲 | 国产黄色免费在线视频| 一级片免费观看大全| 高清黄色对白视频在线免费看| 男女啪啪激烈高潮av片| 精品人妻偷拍中文字幕| 国产黄色免费在线视频| 午夜av观看不卡| 国产麻豆69| 最近最新中文字幕免费大全7| 桃花免费在线播放| 国产成人精品婷婷| 中文精品一卡2卡3卡4更新| 爱豆传媒免费全集在线观看| 精品午夜福利在线看| 男女边吃奶边做爰视频| 少妇人妻 视频| 亚洲国产精品专区欧美| 亚洲精品乱久久久久久| 国产精品人妻久久久影院| 亚洲av.av天堂| 国产欧美另类精品又又久久亚洲欧美| 久久久久视频综合| 成人国产麻豆网| 最近最新中文字幕免费大全7| 2018国产大陆天天弄谢| 精品第一国产精品| 啦啦啦中文免费视频观看日本| 9191精品国产免费久久| 丝袜美足系列| xxxhd国产人妻xxx| 黄色毛片三级朝国网站| 久久精品国产综合久久久 | 咕卡用的链子| 91在线精品国自产拍蜜月| 亚洲高清免费不卡视频| 日本欧美国产在线视频| 2021少妇久久久久久久久久久| 亚洲精品国产av成人精品| 欧美xxxx性猛交bbbb| 国产一级毛片在线| 在线观看人妻少妇| 亚洲人成网站在线观看播放| 亚洲久久久国产精品| 国产精品一区二区在线观看99| 九九在线视频观看精品| 久久女婷五月综合色啪小说| 成人亚洲欧美一区二区av| 亚洲av日韩在线播放| 性高湖久久久久久久久免费观看| 国产精品一区二区在线观看99| 肉色欧美久久久久久久蜜桃| 欧美精品人与动牲交sv欧美| 蜜臀久久99精品久久宅男| 26uuu在线亚洲综合色| 亚洲av福利一区| av一本久久久久| 黄色怎么调成土黄色| 18禁在线无遮挡免费观看视频| 亚洲四区av| 女的被弄到高潮叫床怎么办| 国产日韩欧美在线精品| 岛国毛片在线播放| 哪个播放器可以免费观看大片| 国产激情久久老熟女| 国产免费一区二区三区四区乱码| 亚洲图色成人| 黑丝袜美女国产一区| 亚洲av中文av极速乱| √禁漫天堂资源中文www| 欧美xxxx性猛交bbbb| 午夜免费男女啪啪视频观看| 国产一区二区三区av在线| 亚洲情色 制服丝袜| a级毛色黄片| 色5月婷婷丁香| 如日韩欧美国产精品一区二区三区| 综合色丁香网| 亚洲精品一区蜜桃| 亚洲av电影在线进入| 少妇人妻久久综合中文| 少妇猛男粗大的猛烈进出视频| 晚上一个人看的免费电影| 成年人免费黄色播放视频| 五月伊人婷婷丁香| 人妻人人澡人人爽人人| 日本欧美视频一区| 满18在线观看网站| 久久午夜综合久久蜜桃| 久久狼人影院| 99热网站在线观看| 欧美 日韩 精品 国产| 亚洲精品aⅴ在线观看| 啦啦啦视频在线资源免费观看| 日韩一区二区视频免费看| 国产极品粉嫩免费观看在线| 99热6这里只有精品| 大陆偷拍与自拍| 久久影院123| 国产精品麻豆人妻色哟哟久久| 日本色播在线视频| 中文乱码字字幕精品一区二区三区| 看免费成人av毛片| 中国国产av一级| 少妇精品久久久久久久| 最新中文字幕久久久久| 国产精品久久久久久精品古装| 天堂中文最新版在线下载| 日本av手机在线免费观看| 欧美成人精品欧美一级黄| 久久久国产精品麻豆| av免费在线看不卡| 成年人午夜在线观看视频| 久久狼人影院| 欧美亚洲 丝袜 人妻 在线| 七月丁香在线播放| 丝袜美足系列| a级毛色黄片| 在线观看一区二区三区激情| 亚洲av国产av综合av卡| 午夜福利视频精品| 成年女人在线观看亚洲视频| 一级片免费观看大全| 成人漫画全彩无遮挡| 黄色 视频免费看| 国产精品一区二区在线观看99| 亚洲一级一片aⅴ在线观看| 色婷婷久久久亚洲欧美| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产有黄有色有爽视频| 亚洲色图综合在线观看| 亚洲,欧美精品.| 国产男人的电影天堂91| 婷婷色av中文字幕| 亚洲欧美一区二区三区黑人 | 亚洲欧洲精品一区二区精品久久久 | 十分钟在线观看高清视频www| 欧美 亚洲 国产 日韩一| 久久精品aⅴ一区二区三区四区 | 永久网站在线| 国产男女内射视频| 欧美成人精品欧美一级黄| 亚洲经典国产精华液单| 啦啦啦在线观看免费高清www| 一本大道久久a久久精品| 日韩精品免费视频一区二区三区 | 日本猛色少妇xxxxx猛交久久| 最近的中文字幕免费完整| 国产探花极品一区二区| 国产成人免费观看mmmm| 侵犯人妻中文字幕一二三四区| 国产探花极品一区二区| 日韩三级伦理在线观看| 久久精品久久久久久久性| 国产成人精品在线电影| 男女边吃奶边做爰视频| 另类精品久久| 免费观看av网站的网址| 中国国产av一级| 熟妇人妻不卡中文字幕| 成人毛片60女人毛片免费| 搡老乐熟女国产| 亚洲美女搞黄在线观看| 一级,二级,三级黄色视频| 狂野欧美激情性xxxx在线观看| 在线亚洲精品国产二区图片欧美| 国产精品久久久久久久电影| 免费大片18禁| 在线观看人妻少妇| 一级,二级,三级黄色视频| 黄片播放在线免费| 天天影视国产精品| 91午夜精品亚洲一区二区三区| 久久影院123| 老熟女久久久| 蜜臀久久99精品久久宅男| 毛片一级片免费看久久久久| av在线播放精品| 黑人巨大精品欧美一区二区蜜桃 | 日韩成人伦理影院| 91久久精品国产一区二区三区| 丝袜喷水一区| 最近最新中文字幕免费大全7| 亚洲精品国产av蜜桃| 国产极品天堂在线| 免费日韩欧美在线观看| 你懂的网址亚洲精品在线观看| 免费观看无遮挡的男女| 免费人成在线观看视频色| 亚洲天堂av无毛| 夜夜骑夜夜射夜夜干| 巨乳人妻的诱惑在线观看| 中文欧美无线码| av国产久精品久网站免费入址| 午夜影院在线不卡| 欧美bdsm另类| 国产无遮挡羞羞视频在线观看| 日韩一区二区视频免费看| 观看av在线不卡| 国产精品熟女久久久久浪| 熟女av电影| 国产成人91sexporn| 水蜜桃什么品种好| 国产免费现黄频在线看| 亚洲内射少妇av| 晚上一个人看的免费电影| 全区人妻精品视频| 免费观看性生交大片5| 99热全是精品| 免费大片黄手机在线观看| 如日韩欧美国产精品一区二区三区| 国产精品一区二区在线不卡| 欧美3d第一页| 免费看光身美女| 国产极品天堂在线| 国产片内射在线| 在线免费观看不下载黄p国产| 国产成人一区二区在线| 成人黄色视频免费在线看| 精品久久久精品久久久| 亚洲精品一二三| 最近最新中文字幕免费大全7| 黄网站色视频无遮挡免费观看| 免费在线观看完整版高清| 一区二区av电影网| 精品久久久精品久久久| 欧美激情 高清一区二区三区| 亚洲欧洲国产日韩| 老司机影院成人| 国产成人一区二区在线| 久久久久久伊人网av| 黄网站色视频无遮挡免费观看| 人人妻人人澡人人爽人人夜夜| 97在线视频观看| 国产爽快片一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 精品一区在线观看国产| 高清毛片免费看| 国产精品 国内视频| 嫩草影院入口| 亚洲精品久久成人aⅴ小说| 免费观看a级毛片全部| 久久久久网色| 看免费成人av毛片| 国产精品三级大全| videos熟女内射| 欧美日韩视频精品一区| 人妻少妇偷人精品九色| 日本黄色日本黄色录像| 男女高潮啪啪啪动态图| 18禁动态无遮挡网站| 婷婷色麻豆天堂久久| 99热6这里只有精品| 2022亚洲国产成人精品| 国产成人欧美| 精品国产一区二区久久| av在线老鸭窝| 色哟哟·www| 久久毛片免费看一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 国产精品秋霞免费鲁丝片| 乱人伦中国视频| 精品一区二区三区四区五区乱码 | 看非洲黑人一级黄片| 精品99又大又爽又粗少妇毛片| 又大又黄又爽视频免费| 久久精品国产鲁丝片午夜精品| 女人精品久久久久毛片| 一级爰片在线观看| 色哟哟·www| 免费看不卡的av| 国产精品女同一区二区软件| 精品酒店卫生间| 精品卡一卡二卡四卡免费| 精品国产露脸久久av麻豆| 久久国内精品自在自线图片| a级毛片在线看网站| 欧美成人午夜精品| 日韩中字成人| av在线app专区| 国产精品人妻久久久影院| 日韩电影二区| 精品第一国产精品| 国产精品.久久久| 日韩欧美精品免费久久| 亚洲精品aⅴ在线观看| 久久精品夜色国产| 一区在线观看完整版| 内地一区二区视频在线| 日韩中文字幕视频在线看片| 少妇的逼水好多| 中国国产av一级| 国产探花极品一区二区| 精品卡一卡二卡四卡免费| 肉色欧美久久久久久久蜜桃| 亚洲人与动物交配视频| 日本色播在线视频| 新久久久久国产一级毛片| 久久久久视频综合| 国产又色又爽无遮挡免| av国产精品久久久久影院| 国产精品一区二区在线观看99| 婷婷色综合www| 中文天堂在线官网| 青春草国产在线视频| 波多野结衣一区麻豆| 狠狠精品人妻久久久久久综合| 中文字幕人妻熟女乱码| 免费黄频网站在线观看国产| 欧美少妇被猛烈插入视频| 伦精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 香蕉丝袜av| 人妻少妇偷人精品九色| 欧美日韩视频精品一区| av在线app专区| 国产男女内射视频| 亚洲在久久综合| 18禁在线无遮挡免费观看视频| 成年美女黄网站色视频大全免费| 亚洲精品日本国产第一区| 日韩av不卡免费在线播放| 欧美丝袜亚洲另类| www.色视频.com| 美女中出高潮动态图| 久久ye,这里只有精品| 欧美人与善性xxx| 黄色毛片三级朝国网站| 亚洲国产精品一区二区三区在线| 免费日韩欧美在线观看| 91国产中文字幕| 九九爱精品视频在线观看| 国产一区亚洲一区在线观看| 宅男免费午夜| 久久韩国三级中文字幕| 黄色视频在线播放观看不卡| 美女大奶头黄色视频| 青春草国产在线视频| 日韩精品有码人妻一区| 欧美97在线视频| 国产日韩一区二区三区精品不卡| 亚洲中文av在线| 春色校园在线视频观看| 超色免费av| 熟妇人妻不卡中文字幕| 大片电影免费在线观看免费| 亚洲精品一二三| 卡戴珊不雅视频在线播放| 国产熟女欧美一区二区| 日韩 亚洲 欧美在线| 成人黄色视频免费在线看| 精品亚洲成a人片在线观看| 一级a做视频免费观看| 99国产精品免费福利视频| 咕卡用的链子| 这个男人来自地球电影免费观看 | a 毛片基地| 制服诱惑二区| av在线老鸭窝| av天堂久久9| 久久精品人人爽人人爽视色| 亚洲欧美中文字幕日韩二区| 观看av在线不卡| 久久久久久久亚洲中文字幕| 久久热在线av| 99re6热这里在线精品视频| 青青草视频在线视频观看| 老熟女久久久| 青春草视频在线免费观看| 少妇人妻久久综合中文| 丝瓜视频免费看黄片| 国产黄频视频在线观看| 国产精品一区www在线观看| 国产av精品麻豆| av国产精品久久久久影院| 国产午夜精品一二区理论片| 国产国拍精品亚洲av在线观看| 亚洲伊人色综图| 亚洲国产av新网站| 丝袜在线中文字幕| 久久99蜜桃精品久久| 三级国产精品片| 亚洲精品视频女| 日本色播在线视频| 久久久久精品性色| 中文天堂在线官网| 黄色视频在线播放观看不卡| 国产免费福利视频在线观看| 高清欧美精品videossex| 熟女电影av网| 春色校园在线视频观看| 香蕉丝袜av| 最近手机中文字幕大全| 熟女电影av网| 国产在视频线精品| 国产免费福利视频在线观看| 黄色 视频免费看| 在线观看三级黄色| 在线精品无人区一区二区三| 国产精品久久久久成人av| 在线观看免费日韩欧美大片| 免费黄频网站在线观看国产| 综合色丁香网| 在线看a的网站| 丰满饥渴人妻一区二区三| 少妇人妻久久综合中文| 国产精品不卡视频一区二区| 国产在视频线精品| 如日韩欧美国产精品一区二区三区| 丝袜喷水一区| 久久ye,这里只有精品| 亚洲av电影在线进入| 国产精品99久久99久久久不卡 | 在线精品无人区一区二区三| 一本大道久久a久久精品| 中文欧美无线码| 丁香六月天网| 又粗又硬又长又爽又黄的视频| 午夜福利视频在线观看免费| 中文精品一卡2卡3卡4更新| 久久精品国产亚洲av天美| 搡女人真爽免费视频火全软件| 高清在线视频一区二区三区| 五月伊人婷婷丁香| 久久久久久久大尺度免费视频| av在线观看视频网站免费| 国产乱来视频区| 久久久亚洲精品成人影院| 亚洲精品中文字幕在线视频| 欧美日韩成人在线一区二区| 国产老妇伦熟女老妇高清| 国产成人精品久久久久久| 亚洲四区av| 日韩成人伦理影院| 国产精品三级大全| 精品一区二区三卡| 欧美另类一区| 欧美激情 高清一区二区三区| xxxhd国产人妻xxx| 性高湖久久久久久久久免费观看| 一本久久精品| 精品福利永久在线观看| 麻豆乱淫一区二区| 五月伊人婷婷丁香| 极品少妇高潮喷水抽搐| 国产成人精品婷婷| 国产精品不卡视频一区二区| 看免费av毛片| 99精国产麻豆久久婷婷| 男人爽女人下面视频在线观看| 18在线观看网站| 美女视频免费永久观看网站| 日韩制服骚丝袜av| 巨乳人妻的诱惑在线观看| 91久久精品国产一区二区三区| 久久精品国产亚洲av涩爱| 18禁在线无遮挡免费观看视频| 欧美日韩一区二区视频在线观看视频在线| 麻豆精品久久久久久蜜桃| 精品久久蜜臀av无| 国产白丝娇喘喷水9色精品| 青春草国产在线视频| 亚洲av福利一区| 国产免费现黄频在线看| 国产亚洲一区二区精品| 国产精品人妻久久久影院| 国产精品.久久久| 国产乱来视频区| 国产成人精品婷婷| 欧美最新免费一区二区三区| 韩国高清视频一区二区三区| 国产国拍精品亚洲av在线观看| 久久女婷五月综合色啪小说| 妹子高潮喷水视频| 男女边吃奶边做爰视频| 国产精品久久久久久久电影| 久久久久久久大尺度免费视频| 亚洲欧洲精品一区二区精品久久久 | 2021少妇久久久久久久久久久| 26uuu在线亚洲综合色| 欧美老熟妇乱子伦牲交| 国产老妇伦熟女老妇高清| 日韩制服骚丝袜av| 国产探花极品一区二区|