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

    垂直平面內(nèi)二維超聲振動(dòng)銑削系統(tǒng)穩(wěn)定性研究

    2017-07-01 23:15:00趙斌斌范凱洋張躍敏
    關(guān)鍵詞:葉瓣切削力刀具

    趙 波, 趙斌斌, 范凱洋, 張躍敏

    (河南理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院 焦作, 454000)

    ?

    垂直平面內(nèi)二維超聲振動(dòng)銑削系統(tǒng)穩(wěn)定性研究

    趙 波, 趙斌斌, 范凱洋, 張躍敏

    (河南理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院 焦作, 454000)

    為研究垂直平面內(nèi)二維超聲振動(dòng)銑削系統(tǒng)的穩(wěn)定性,建立了垂直平面內(nèi)二維超聲振動(dòng)銑削穩(wěn)定性模型,采用全離散法對(duì)二維超聲振動(dòng)銑削穩(wěn)定性進(jìn)行了研究,利用Matlab軟件進(jìn)行了數(shù)值仿真,獲得了系統(tǒng)的穩(wěn)定性葉瓣圖。銑削鈦合金材料顫振實(shí)驗(yàn)結(jié)果表明,主軸轉(zhuǎn)速在1 000~3 500 r/min范圍內(nèi),數(shù)值仿真結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好,驗(yàn)證了垂直平面內(nèi)二維超聲振動(dòng)銑削穩(wěn)定性模型的正確性。在垂直平面內(nèi)施加二維超聲振動(dòng)能提高系統(tǒng)的穩(wěn)定性,軸向切深極限最大提高了約13.6%。

    二維超聲振動(dòng)銑削;全離散法;穩(wěn)定性;葉瓣圖

    引 言

    鈦合金是一種具有強(qiáng)度高、密度小、耐腐蝕性好及耐熱性高等特點(diǎn)的金屬,被廣泛應(yīng)用于航空航天等領(lǐng)域。雖然鈦合金具有上述優(yōu)點(diǎn),但其存在工藝性能差、切削加工困難等缺點(diǎn)。為了解決鈦合金等難加工材料的加工問題,業(yè)內(nèi)普遍采用的做法是合理選擇合適的加工設(shè)備和加工參數(shù),來提高難加工材料的加工質(zhì)量和加工效率。

    高速銑削具有高精度、高效率等優(yōu)點(diǎn),但在加工過程中,如果加工參數(shù)選擇不當(dāng),則存在切削顫振。切削顫振不僅嚴(yán)重影響著零件的加工質(zhì)量,還可能會(huì)破壞機(jī)床設(shè)備及刀具。因此,許多學(xué)者都對(duì)切削顫振進(jìn)行了研究。Budak等[1-2]運(yùn)用頻率法預(yù)測(cè)銑削穩(wěn)定性。Altintas等[3]考慮了用時(shí)域周期系數(shù)矩陣的時(shí)滯微分方程組描述再生效應(yīng)的動(dòng)態(tài)銑削過程,利用頻率法得出銑削系統(tǒng)穩(wěn)定性。Insperger[4]提出在時(shí)域內(nèi)運(yùn)用半離散法,并對(duì)一自由度和兩自由度銑削模型進(jìn)行了銑削穩(wěn)定性分析。Ding等[5]提出了能有效提升計(jì)算效率的全離散法,并分析了銑削系統(tǒng)穩(wěn)定性。姜燕等[6]提出了一種時(shí)間有限元預(yù)測(cè)法預(yù)測(cè)銑削系統(tǒng)的穩(wěn)定性。上述都是基于Floquet理論,將無窮維時(shí)滯系統(tǒng)離散為有限維系統(tǒng),利用系統(tǒng)傳遞矩陣特征值的模小于1時(shí)系統(tǒng)穩(wěn)定性的原理,但都需要大量實(shí)驗(yàn)驗(yàn)證其理論的正確性,因此存在經(jīng)濟(jì)性差的問題。

    超聲振動(dòng)加工是一種能夠提高加工系統(tǒng)穩(wěn)定性的方法[7],被廣泛應(yīng)用在精密加工領(lǐng)域中。速度系數(shù)k是超聲振動(dòng)加工的重要參數(shù)之一,其將超聲振動(dòng)加工分為分離型(k<1)和不分離型(k≥1)。文獻(xiàn)[8-9]研究了不分離型超聲橢圓振動(dòng)切削力,研究表明,超聲橢圓振動(dòng)切削在不分離區(qū)仍然能夠有效降低切削力。唐軍等[10]研究了分離型超聲縱扭復(fù)合銑削系統(tǒng)穩(wěn)定性的影響,并通過銑削加工碳纖維復(fù)合材料(C/C)進(jìn)行實(shí)驗(yàn)驗(yàn)證。超聲振動(dòng)加工具有經(jīng)濟(jì)效益好、裝置簡(jiǎn)單等優(yōu)點(diǎn),但目前對(duì)不分離型超聲振動(dòng)銑削系統(tǒng)穩(wěn)定性的研究還不是很多。

    考慮到目前還沒有利用全離散法對(duì)二維超聲振動(dòng)銑削系統(tǒng)穩(wěn)定性進(jìn)行研究,因此,筆者綜合利用全離散法和超聲振動(dòng)加工法,提出垂直平面內(nèi)二維超聲振動(dòng)銑削模型,建立垂直平面內(nèi)二維超聲振動(dòng)銑削系統(tǒng)的穩(wěn)定性模型。利用全離散法分析二維超聲振動(dòng)銑削穩(wěn)定性,應(yīng)用Matlab軟件進(jìn)行數(shù)值分析并得出二維超聲銑削穩(wěn)定性葉瓣圖。最后,通過二維超聲振動(dòng)銑削鈦合金實(shí)驗(yàn)驗(yàn)證二維超聲復(fù)合銑削穩(wěn)定性模型的正確性。

    1 二維超聲振動(dòng)銑削穩(wěn)定性模型

    考慮到超聲振動(dòng)分別施加在工件和刀具上,如圖1(a)所,在xOy平面上由于超聲振動(dòng)的存在,因此當(dāng)速度系數(shù)k≥1時(shí),屬于不分離型超聲振動(dòng)加工;k<1時(shí),則屬于分離型超聲振動(dòng)加工。速度系數(shù)表示為

    (1)

    其中:v為刀尖相對(duì)于工件的瞬時(shí)線速度;vc為臨界切削速度;A為超聲振幅;f為超聲振動(dòng)頻率。

    但分離型振動(dòng)加工還要滿足必要條件[11]:λ為奇數(shù),A>fz,其中λ為wc/wz,wc和wz分別表示周轉(zhuǎn)轉(zhuǎn)角速度和超聲波發(fā)生器角頻率,fz為超聲銑削系統(tǒng)進(jìn)給量。根據(jù)圖1,若不考慮刀具齒位角的變化和沿軸向超聲振動(dòng)切削厚度的影響,則考慮再生效應(yīng)的二維超聲振動(dòng)切削的切厚表達(dá)式為

    hj(t)=g(t)[sc]T(hjs(t)+hjd(t))

    (2)

    其中:g(t)為超聲振動(dòng)切削分離判定系數(shù),當(dāng)銑刀在切削工件時(shí)且滿足分離型振動(dòng)加工的必要條件時(shí),g(t)=1,否則為零;s=sin(Φj(t));c=cos(Φj(t));Φj為刀尖轉(zhuǎn)角;hjs(t)為靜態(tài)切厚,hjs(t)=[fz+Asin(2πf1t+φ1) 0]T;f1為進(jìn)給方向上超聲振動(dòng)頻率;φ1為初始相位;hjd(t)為動(dòng)態(tài)切厚,hjd(t)=[x(t)-x(t-τ)y(t)-y(t-τ)]T;τ=2π/(zwc)。

    根據(jù)式(2),采用指數(shù)型的二維超聲振動(dòng)銑削力學(xué)模型[12-13]可簡(jiǎn)化表示為

    (3)

    其中:M,C和K分別為模態(tài)質(zhì)量、模態(tài)阻尼和模態(tài)剛度;q(t)=[x(t)y(t)]T為刀齒的動(dòng)態(tài)位移;B,f1和φ2分別為軸向超聲振動(dòng)的振幅、超聲振動(dòng)頻率和初始相位。

    當(dāng)Φst(t)<Φj(t)<Φet(t)時(shí),w(Φj(t))=1;否則w(Φj(t))=0。Φet(t)和Φet(t)分別為刀具切入角和切出角。

    dv/dt=A0v(t)+A(t)v(t)+B(t)v(t-T)

    (4)

    其中:A0為常數(shù)矩陣;A(t)和B(t)分別代表周期函數(shù)矩陣,且A(t)=-B(t);T=60/(ZΩ),Ω為主軸轉(zhuǎn)速。

    在主軸旋轉(zhuǎn)的一個(gè)刀齒周期T內(nèi),把T離散為m等分,把矩陣A(t),v(t)和v(t-T)分別在每個(gè)離散的小區(qū)間上進(jìn)行拉格朗日插值法代替,根據(jù)Floquet理論,使傳遞函數(shù)特征值的模小于1,則可得到二維超聲銑削系統(tǒng)的穩(wěn)定性葉瓣圖。由于筆者研究的是超聲不分離型振動(dòng)銑削,即此時(shí)g(t)=1。

    圖1 二維超聲振動(dòng)銑削Fig.1 The ultrasonic vibration milling of two dimensional

    2 穩(wěn)定性模型實(shí)驗(yàn)驗(yàn)證

    2.1 二維超聲振動(dòng)銑削系統(tǒng)穩(wěn)定性預(yù)測(cè)圖

    圖2 實(shí)驗(yàn)裝置Fig.2 The devices of experiment

    圖3 穩(wěn)定性葉瓣圖Fig.3 The diagram of the stability lobes

    圖3表明,與普通銑削相比,垂直平面內(nèi)的二維超聲振動(dòng)銑削系統(tǒng),其最大提高銑削系統(tǒng)切削深度的極限值約為13.6%,同時(shí)也增大了穩(wěn)定性的區(qū)域,即采用全離散法分析的二維超聲振動(dòng)銑削系統(tǒng)不僅可以預(yù)測(cè)銑削系統(tǒng)的穩(wěn)定性,同時(shí)還可以提高銑削系統(tǒng)的穩(wěn)定性。這是因?yàn)樵阢娤鬟^程中,由于刀具沿主軸方向和工件沿刀具進(jìn)給方向的高頻振動(dòng)(兩個(gè)方向的超聲振動(dòng)頻率遠(yuǎn)大于主軸轉(zhuǎn)速),破環(huán)了前后刀齒留下的加工振紋,從而降低了顫振的幾率。此外,在兩個(gè)方向超聲振動(dòng)的作用下,刀具與工件為不分離切削,同時(shí)刀尖在工件上呈現(xiàn)多維運(yùn)動(dòng),這就使得加工區(qū)域更能充分地接觸到切削液,從而改善了工況,更有利于切削加工。由于刀具在主軸方向上具有超聲振動(dòng),造成加工表面形成點(diǎn)狀的超聲振動(dòng)凹坑,這些凹坑的存在影響著零件加工表面的應(yīng)力分布。

    2.2 穩(wěn)定性實(shí)驗(yàn)驗(yàn)證

    在不同主軸轉(zhuǎn)速下驗(yàn)證垂直平面內(nèi)的二維超聲振動(dòng)銑削穩(wěn)定性葉瓣圖,如圖4所示,其中“0”,“?”和“×?”等分別表示系統(tǒng)穩(wěn)定點(diǎn)、顫振點(diǎn)和無法判斷銑削穩(wěn)定性的點(diǎn)。圖4中A,B和D點(diǎn)在普通銑削中為不穩(wěn)定銑削點(diǎn);而在二維超聲振動(dòng)型銑削下A點(diǎn)轉(zhuǎn)為臨界點(diǎn),B和D點(diǎn)則變成穩(wěn)定銑削點(diǎn);C點(diǎn)在兩種加工條件下均為穩(wěn)定性點(diǎn)。

    圖4 銑削穩(wěn)定性預(yù)測(cè)圖及實(shí)驗(yàn)驗(yàn)證參數(shù)點(diǎn)Fig.4 The milling stability prediction diagram and experimental verification parameters

    為了比較二維超聲振動(dòng)銑削與普通銑削的區(qū)別,從圖4中選擇C(1 750r/min,0.15 mm)點(diǎn)進(jìn)行研究,由于y方向沒有施加超聲振動(dòng),因此只有普通銑削和二維超聲振動(dòng)銑削條件下的x,z向切削力時(shí)域信號(hào)、頻譜圖以及加工后的零件形貌圖,如圖5和圖6所示。

    由圖5和圖6可以看到,普通銑削下的切削力和二維超聲振動(dòng)銑削下的切削力大小基本相等。但從頻譜圖中可以明顯看出,二維超聲振動(dòng)銑削下的諧波分量明顯被抑制,因此提高了銑削系統(tǒng)的穩(wěn)定性,驗(yàn)證了筆者所提出的不分離型超聲復(fù)合銑削穩(wěn)定性模型的正確性。

    3 結(jié)束語

    筆者利用全離散法和超聲振動(dòng)加工法,提出垂直平面內(nèi)二維超聲振動(dòng)銑削模型,建立了垂直平面內(nèi)二維超聲振動(dòng)銑削穩(wěn)定性模型。運(yùn)用全離散法分析的垂直平面內(nèi)二維超聲振動(dòng)銑削系統(tǒng)能夠準(zhǔn)確地預(yù)測(cè)銑削系統(tǒng)穩(wěn)定性葉瓣圖,通過在數(shù)控加工中心VMC850E銑削鈦合金實(shí)驗(yàn),驗(yàn)證了垂直平面內(nèi)二維超聲振動(dòng)銑削穩(wěn)定性模型和穩(wěn)定性葉瓣圖的正確性,同時(shí)還能夠提高普通銑削系統(tǒng)的穩(wěn)定性。該實(shí)驗(yàn)裝置是在自行研制的聲學(xué)設(shè)備的基礎(chǔ)上改造而來,因此,大大提高了設(shè)備的利用率和加工的經(jīng)濟(jì)效益。

    圖5 C點(diǎn)普通銑削力信號(hào)及工件表面形貌Fig.5 The signals of traditional milling forces and the workpiece appearance of C

    圖6 C點(diǎn)二維超聲銑削力信號(hào)及工件表面形貌Fig.6 The signals of two-demension ultrasonic milling forces and the workpiece appearance of C

    [1] Budak E, Altintas Y. Analytical prediction of chatter stability in milling—part I:general formulation[J]. Journal of Dynamic System, Measurement, and Control, 1988,120:22-30.

    [2] Budak E. An analytical design method for milling cutters with nonconstant pitch to increase stability, part 2: application[J]. Journal of Manufacturing Science and Engineering, 2003, 125:35-38.

    [3] Altintas Y, Stepan G , Merdol D, et al.Chatter stability of milling in frequency and discrete time domain [J].CIRP Journal of Manufacturing Science and Technology,2008, 1(1): 35-44.

    [4] Insperger T. Updated semi-discretization method for periodic delay-differential equation with discrete delay[J]. International Journal for Numerical Mehods in Engineering,2004,61: 117-141.

    [5] Ding Ye, Zhu Limin, Zhang Xiaojian, et al. A full-discretization method for prediction of milling stability[J]. International Journal of Machine Tools & Manufacture, 2010,50:502-509.

    [6] 姜燕,郭強(qiáng),趙波.銑削穩(wěn)定性預(yù)測(cè)的時(shí)間有限元法[J].河南理工大學(xué)學(xué)報(bào):自然科學(xué)版,2016,35(5):672-676.

    Jiang Yan,Guo Qiang,Zhao Bo. The method on the stability limit prediction for milling process based on time-finite-element theory [J]. Journal of Henan Polytechnic University: Natural Science, 2016,35(5):672-676.(in Chinese)

    [7] 于勁,周小勤. 基于高頻變速特征的不分離型超聲波振動(dòng)車削抑制顫振機(jī)理[J].兵工學(xué)報(bào),1993,14 (1) : 52-57.

    Yu Jin,Zhou Xiaoqin. On the mechanism of chatter suppression with high frequency and vari-speed unseparated type ultrasonic vibration turning[J].Acta Armamentarii,1993,14(1) : 52-57. (in Chinese)

    [8] 李勛,張德遠(yuǎn). 不分離型超聲橢圓振動(dòng)切削實(shí)驗(yàn)研究[J]. 機(jī)械工程學(xué)報(bào),2010,46(19):177-182.

    Li Xun, Zhang Deyuan. Experimental study on the unseparated ultrasonic elliptical vibration cutting [J]. Journal of Mechanical Engineering, 2010,46(19):177-182.(in Chinese)

    [9] 李文,尹礁,呂壘平,等. 不分離型超聲橢圓振動(dòng)切削力特性研究[J]. 航空學(xué)報(bào),2013,34(9): 2241-2248.

    Li Wen, Yin Jiao, Lü Leiping, et al. Study on the unseparated ultrasonic elliptical vibration cutting force[J]. Acta Aeronautica et Astronnutica Sinica, 2013,34(9):2241-2248.(in Chinese)

    [10] 唐軍,趙波. 分離型縱-扭復(fù)合超聲銑削的穩(wěn)定性分析[J]. 兵工學(xué)報(bào),2015,36(7):1318-1325.

    Tang Jun, Zhao Bo. Stability analysis of the separated longitudinal-torsional composite ultrasonic milling [J]. Acta Armamentarii, 2015, 36(7): 1318-1325.(in Chinese)

    [11] 張建華. 超聲振動(dòng)輔助微細(xì)銑削運(yùn)動(dòng)學(xué)分析[EB/OL]. (2012-12-28) [2016-09-10].http:∥www.paper.edu.cn/releasepaper/content/201212-1173.

    [12] Faassen R P H, van de Wouw N, Oosterling J A J, et al. Prediction of regenerative chatter by modelingand analysis of high-speed milling[J].International Journal of Machine Tools & Manufacture,2003,43:1437-1446.

    [13] 宋清華,艾興,萬熠,等. 小徑向切深下的進(jìn)給量對(duì)銑削穩(wěn)定性的影響[J].中國(guó)機(jī)械工程,2008,19(10):1148-1152.

    Song Qinghua, Ai Xing, Wan Yi, et al. Stability prediction for high-speed milling including feed rate in low radial immersion [J]. Chinese Journal of Mechanical Engineering, 2008, 19(10):1148-1152.(in Chinese)

    Gear Crack Damage Localization of the Planetary Gearbox Based on Ensemble Empirical Mode Decomposition

    LIUHaohua1,LIFangyi1,LIGuoyan1,WANGYifan2,ZHANGShanshan1,DONGDehao1

    (1.Key Laboratory of High Efficiency Clean Mechanical Manufacture,Ministry of Education, Shandong University Jinan, 250061, China)(2.Shanghai Aerospace Equipments Manufacturer Shanghai, 200245, China)

    Abstract In order to solve the problem that it is difficult to extract the characteristic frequency and identify the location of the gear crack damage in planetary gearboxes, a gear local damage diagnosis method based on ensemble empirical mode decomposition (EEMD) and frequency demodulation analysis is put forward. Based on the vibration signal model of typical gear local damages, the method can process vibration signals of the sun, ring and planet gear with crack damage by EEMD and frequency demodulation analysis. The local damage characteristic frequency of the gear is extracted from the frequency spectrum, and the location of the gear with crack damage in the gearbox is identified. Simulation analysis and experimental results show that the method based on EEMD and frequency demodulation analysis can effectively extract the characteristic frequency from the sun, ring and planet gear with crack damage. Besides, the damage localization of the gear with crack damage in the planetary gearbox is completed.

    Keywords gearbox; frequency demodulation; damage localization; ensemble empirical mode decomposition

    Vibration Control Analysis of Trailing Edge Flap Smart Rotor with Limited Deflection Angle

    LIUShiming,YANGWeidong,YUZhihao,DONGLinghua,LUKaihua
    (National Key Laboratory of Rotorcraft Aeromechanics, Nanjing University of Aeronautics and Astronautics Nanjing, 210016, China)

    Abstract An optimization method for vibration reduction of small-scaled smart rotor with trailing edge flaps is presented. Both the inertial forces and aerodynamic forces due to the deflection of trailing edge flaps are concerned in this model. A surrogate model is developed to calculate the aerodynamic forces of flapped airfoils. The aeroelastic dynamic equations are solved with the implicit trapezoid method to get the elastic response of blade, and the vibratory blade loads and hub loads are predicted with a force integration method. The flap deflection harmonics are the design variables and the amplitutes of vibratory hub load are chosen as the objective function. The best flap deflection law for hub vibration control is found with the steepest descent method. Results show that both structural and aerodynamic loads of rotor can be precisely calculated with the current model. Hub vertical vibratory load can be effectively reduced with properly controlled flaps at different advance ratios. The lack of deflection ability can be simulated with the direct constraint method or the objective weight method. The deflection ability of trailing edge flaps significantly influences the vibration reduction effect. Despite the limited deflection angle due to the ability of actuator, vibratory loads can still be reduced with actively controlled flaps.

    Keywords helicopter; rotor; vibration; trailing edge flap; optimization method

    Research on Ultrasonic Guided Wave-Based Damage Localization for Pipeline Structure

    WANGGuofeng,LIFucai,LIUZhiqiang,MENGGuang
    (State Key Laboratory of Mechanical System and Vibration, Shanghai Jiaotong University Shanghai, 200240, China)

    Abstract Structure damage will affect the propagation of ultrasonic guided waves. Hence, non-destructive testing technology based on ultrasonic guided waves can realize structural health monitoring (SHM). A pipe with 174mm inner-diameter and 194mm outer-diameter, made of 20#carbon steel is used to investigate the damage localization in this study. Based on dispersion equations, numerical method is applied to obtain longitudinal and circumferential guided waves in the pipe structure. Central frequency of the incident signal is selected 80 kHz by simultaneously considering the dispersion curves and wave structure. Finite element method (FEM) is used to verify characteristics of guided wave propagation. According to the FEM analyses, characteristics and wave structures of longitudinal modes are close to those of circumferential waves, moreover, both of them are similar to those of Lamb waves in plate structure. Elliptical damage localization method is therefore proposed with one excitation and multiple captures for damage localization in the structure. Notch and hole are introduced to the pipe to verify effectiveness of the proposed method, and then analyzed influencing factors of damage localization deviation.

    Keywords guided wave; pipe; dispersion; damage localization; structural health monitoring (SHM)

    A Forecasting Method of Positioning Accuracy for CNC Machine Tools Feed System Based on BP Neural Network

    DENGChao1,QIANYousheng1,WUJun2,XIONGYao3,DUANChaoqun1
    (1.State Key Lab of Digital Manufacturing Equipment & Technology,

    Huazhong University of Science & Technology Wuhan, 430074, China)

    (2.Institute of Ship and Ocean Engineering, Huazhong University of Science and Technology Wuhan, 430074, China)

    (3.Wuhan Second Ship Design and Research Institute Wuhan, 430070, China)

    Abstract Aimed at the difficulty in forecasting the positioning accuracy of machine tool feed system,a new model has been established, analyzing the positioning accuracy degradation of the mechanical transmission system of feed system. Firstly, a computer numerical control(CNC) machine tool feed system model is made using a dynamic simulation software Adams to obtain positioning accuracy values for different initial state. Then, a mapping model is set up between the positioning accuracy and the gap, the ball tilt and the workpiece loading by BP neural network. According to the mapping model, the prediction method of the positioning accuracy is discussed in detail. Finaly, the method is proved to be correct and effective by conducting the test experiment of positioning accuracy on the test platform of precision motion reliability.

    Keywords feed system; BP neural network; mapping model; Adams; positioning accuracy prediction

    Research on Impact Location by Using Fiber Bragg Grating Sensor Network

    CAOLiang1,WANGJinglin1,HEZhaohua1,LIANGDakai2,SHANTianmin1,LINZeli1
    (1.Aviation Key Laboratory of Science and Technology on Fault Diagnosis and Health Management, Shanghai Aero Measurement & Control Technology Research Institute Shanghai, 201601, China)

    (2.State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics Nanjing, 210016, China)

    Abstract In view of the demand of monitoring the impact and vibration of laminated composite plates, a new method of low speed impact identification based on wavelet packet decomposition and distributed fiber Bragg grating sensing network is proposed. The sensor network layout is designed according to the load form of clamped plate structure and the sensing characteristics of fiber Bragg grating sensor. Then, the impact response signals monitored by fiber Bragg grating sensor network are analyzed by the fast Fourier transformation and wavelet packet decomposition to obtain the time domain characteristics of impact characteristics. On this basis, the mutual correlation coefficient between each characteristic decomposition signal and its corresponding time domain original signal is calculated to be the weight of the similarity. The characteristics of the impact response signals of all samples are decomposed, and the sample database is constructed. Finally, the location coordinates of impact points are determined according to the similarity calculated based on the Haudorff distance between the test signal and the sample information. The research shows that this method can realize the identification of the low speed impact position on the structure of the laminated plates.

    Keywords fiber Bragg grating sensor network; laminated plate; wavelet packet decomposition; impact location

    The Coupling Dynamic Model and Vibration Response of Straddle Type Monorail Vehicle

    WENXiaoxia1,DUZixue2,XUZhouzhou1,YINYanli1,WEIHanbing1
    (1.School of Electrical and Vehicle Engineering, Chongqing Jiaotong University Chongqing, 400074, China)

    (2.Institute of Urban Rail, Chongqing Jiaotong University Chongqing, 400074, China)

    Abstract In order to investigate the vibration response characteristics and evaluate the stability and comfort of the straddle type monorail vehicle within the designed speed, the topology relationships among the traveling wheels, guide wheels, stabilizing wheels and bogie central traction devices are analyzed. A monorail vehicle spatial coupling dynamic model, which consists of three wheel-rail contact modes, is established based on Hamilton equation. The contacts to be studied are the one between the traveling wheel and the top surface of beam, and the contact of the side surface of the rail beam with the guide wheel and stabilizing wheel respectively. Several prominent characteristics such as the vibration dynamic responses of monorail lateral and vertical acceleration are obtained based on the coupling dynamic model when the roughness is designated to be the excitation source given the vechicle moving with a constant speed. Finally, the coupling dynamic model is validated with the comparison of the related experimental results. In this way, vibration response characteristics of monorail vehicle for different design speed are obtained under three rail surface roughness excitation. According to the simulated vibration results, the monorail vehicle stability and comfort performance through straight line is evaluated. The results infer that the monorail vehicle had outstanding comfort performance and excellent running stability.

    Keywords monorail vehicle; coupling dynamic model; contact model; rail uneven excitation; vibration; stability performance

    Laplacian Eigenmaps-Support Vector Domain Description Method for Complex Electromechanical System

    YASENJIANG·Jiarula1,2,GAOJianmin1,GAOZhiyong1,JIANGHongquan1.CHENZisheng1
    (1.State Key Laboratory for Manufacturing Systems Engineering, Xi′an Jiaotong University Xi′an, 710049, China)

    (2.School of Mechanical Engineering, Xinjiang University Urumqi, 830046, China)

    Abstract The monitoring data of the complex electromechanical system has obvious high-dimensional nonlinear and complex distribution characteristics. In order to meet the requirements of complex system anomaly identification which are difficult to be satisfied by the traditional method, a kind of Laplacian eigenmaps-support vector domain description method (LE-SVDD) is proposed. On account of the fact that the points which are close in the high-dimensional feature space should also be close after being projected to the low-dimensional feature space, the improved LE method uses a weighted undirected graph to describe a popularity and find the low-dimensional embedment with an embedded method, thereby status popular structures can be found in the high-dimensional data. In the simulation experiments based on the standard Tennessee - Eastman process (TE process) test and training data, the accurate results of nonlinear feature extraction and anomaly identification at different time are given. Respectively, the average false negative rate and false alarm rate are 6.063, 6 and 5.625, 3.125, which are relatively low. It shows that LE-SVDD method has excellent non-linearity and high-dimensional data processing capability in condition monitoring, which is suitable for the monitoring and diagnosis of engineering system.

    Keywords complex electromechanical systems; anomaly detection method; feature extraction; Laplacian eigenmaps-support vector domain description method (LE-SVDD); Tennessee - Eastman proces

    Rolling Bearing Reliability Assessment and Life Prediction Based on KPCA and WPHM

    WANGFengtao,CHENXutao,LIUChenxi,LIHongkun,HANQingkai,ZHUHong
    (School of Mechanical Engineering, Dalian University of Technology Dalian, 116024, China)

    Abstract The remaining useful life (RUL) prediction of rolling bearing is significant for proactive maintenance of equipment, and selecting the features which can accurately reflect the performance degradation process as the inputs of the life prediction model is the premise of accurate RUL prediction. A novel method based on kernel principal component analysis (KPCA) and Weibull proportional hazard model (WPHM), is proposed to assess the reliability and predict the RUL of the rolling bearing. High relative feature set is constructed by selecting the effective features through extracting the time domain, frequency domain and time-frequency domain features of lifetime bearing. The kernel principal components (KPCs) which can accurately reflect the performance degradation process are obtained by KPCA. Then the KPCs are used as the covariates of WPHM to assess the reliability and predict the RUL. An example of bearing test is provided to demonstrate that this method can accurately assess the reliability and predict the RUL to provide timely maintenance resolution. Meanwhile, as the relative features are extracted, the differences in manufacturing, installation and working condition of the same type bearings are reduced, which enhances the practicability and stability of the method.

    Keywords rolling bearing; life prediction; kernel principal component analysis; Weibull proportional hazard model; relative feature

    Fault Diagnosis of Wind Turbine Gearbox Based on KFCM Optimized by Particle Swarm Optimization

    LIZhuang1,LIUYibing1,TENGWei1,LINYang1,2
    (1.School of Energy, Power and Mechanical Engineering, North China Electric Power University Beijing,102206, China)

    (2.Beijing Huaneng Xinrui Control Technique Co.,Ltd Beijing, 102209, China)

    Abstract A method based on kernel fuzzy c-means clustering (KFCM) optimized by particle swarm optimization is proposed for fault diagnosis of wind turbine gearbox. Firstly, the clustering model is built based on wrong classification rate of training samples. The training samples are classified by kernel fuzzy c-means clustering. Then particle swarm optimization is introduced for solving the clustering model while the initial clustering center and parameter of kernel function are chosen as optimization variables. The class centers of optimal clustering result are acquired. Finally, the similarity parameters in kernel space between new data samples and the class centers are calculated for diagnosing whether the new data sample belongs to knows faults. The results show that the proposed method can diagnose both the known faults and unknown faults effectively compared to traditional neural network based on supervised learning.

    Keywords kernel fuzzy c-means clustering; particle swarm optimization; wind turbine; gearbox; fault diagnosis

    A Single-Mode-Drive and Tower-Shaped Ultrasonic Motor with an Asymmetric Stator

    CHENQianwei,JUQuanyong,GAOSumeiYANGChuansen
    (School of Mechatronic Engineering, Jinling Institute of Technology Nanjing, 211169, China)

    Abstract In order to solve the problem that the existing single-mode-drive USMs either move unidirectionally or wear and tear seriously, a single-mode-drive and tower-shaped USM that can move bi-directionally is presented. The USM is composed of an asymmetric stator and a mover, and the stator is designed as an asymmetric Langevin vibrator. A low-order and a higher-order asymmetric working mode is utilized to the stator, and the scheme of PZTs′ polarization and arrangement location are designed accordingly. Switching between the low-order and the higher-order mode can drive the motor bi-directionally by the single mode. Firstly, the working principle of the motor is analyzed. Then the prototype of the motor is fabricated. Finally, the tests of vibration mode and mechanical characteristics on the prototype are accomplished. Experimental results show that under the condition of single phase excitation through phase A, the motor is working under the low-order mode, and the mover is driven forward with the maximal velocity of 112 mm/s and the maximal driving force of 2N; under the condition of single phase excitation through phase B, the motor is working under the higher-order mode, and the mover is driven backward with the maximal velocity of 94 mm/s and the maximal driving force of 3 N.

    Keywords single mode; asymmetric; ultrasonic motor; piezoelectric

    Damage Identification of Transmission Tower Based on Cross Correlation Functions Amplitude and Support Vector Machine

    HUOLinsheng1,LIXu1,LIHongnan1,ZHANGZhuoqun2
    (1.Faculty of Infrastructure Engineering, Dalian University of Technology Dalian, 116023, China)

    (2.State Nuclear Electric Power Planning Design & Research Institute Beijing, 100095, China)

    Abstract A damage detection method for transmission tower is presented based on integrating the cross correlation function amplitude and the support vector machine (SVM). The proposed method consists of two stages. Firstly, the data features, which are defined as the difference of the cross correlative function amplitude between initial and current states, is calculated from the approximate signal of single mode responses. The data features are then used to train the SVM classifier to turn the damage identification problem into classification problem. The feasibility of the method is verified by the vibration test of the transmission tower model. The results show that the proposed method only uses the dynamic responses from small number of sensors, and can identify the damage status with strong noise robustness during arbitrary excitations.

    Keywords the cross correlative function; support vector machine; damage identification; transmission tower

    The Yaw Error Compensation Method of Large Stroke Precision Positioning Platform

    LIUJizhu1,2,LIJian1,2,ZHANGWenwen1,2,WANGYangjun1,2,PANMinqiang1,2,CHENLiguo1,2
    (1.School of Mechanical and Electric Engineering, Soochow University Suzhou, 215021, China)

    (2.Collaborative Innovation Center of Suzhou Nano Science and Technology, Soochow University Suzhou, 215021, China)

    Abstract In order to reduce the impact of high precision positioning platform yaw error and improve the positioning accuracy. This paper studies the yaw error compensation method of the platform based on linear motor drive, air support and guide rails. According to the dynamic characteristics of the positioning platform's yaw errors, a highly frequent responding and non-contact yaw error detection system combined flat feet and micro displacement sensor is designed, a detection algorithm model is established, and anx,ytwo-dimensional micro-displacement compensation mechanism based on piezoelectric ceramic driver is designed. Using spring mass damper system simplified by flotation rail slider deputy, positioning platform deflection of two degrees of freedom vibration model iss established. The deflection error compensation experiments based on the error compensation system shows that positioning accuracy of precision positioning platform has been greatly improved, the final positioning accuracy of the positioning platform is better than 2 μm.

    Keywords yaw error; vibration model; pb-based lanthanumdoped zirconate titanates (PZT); error compensation

    Research on Lamb Wave and Linear PZT Array Scanning Based on Directional Damage Imaging and Evaluation

    WANGQiang
    (College of Automation,Nanjing University of Posts and Telecommunications Nanjing, 210023, China)

    Abstract Damage evaluation method is one of the key techniques in structural health monitoring. Research efforts are concentrated on Lamb wave based directional damage quantitative monitoring and evaluation method. Based on linear pb-based lanthanumdoped zirconate titanates (PZT) sensor array and phase information in different sensing signals of PZT elements, phase delay of directional damage reflection could be scanned and detected to decide the orientation and central position of damage. Time reversal damage imaging method is also improved based on linear sensor array tuning after the orientation of the damage is decided to realize the directional damage imaging and evaluation. Experiments on aluminum plate indicate the effect of the new method on damage orientation and dimension quantitative evaluation and the better anti-interference and accuracy than existing damage imaging methods.

    Keywords structural health monitoring; damage evaluation; imaging; directional damage; Lamb waves

    Modeling and Simulating of a Two-Stage Pressure Hydro-pneumatic Spring for Off-Road Vehicle

    LIZhongxing1,GUOZiquan1,WANGChuanjian2,LIMei3,MAZili1
    (1.School of Automotive and Traffic Engineering, Jiangsu University Zhenjiang, 212013, China)

    (2.Jingjiang College, Jiangsu University Zhenjiang, 212013, China)

    (3.Mechanical and Electrical Engineering College, Hainan University Haikou, 570228, China)

    Abstract In order to relieve the contradiction of dynamic performance between different loading conditions, a two-stage pressure hydro-pneumatic spring which is designed for off-load vehicle is proposed. Taking the compressibility of oil into consideration, a nonlinear mathematical model is established. A test bench for hydro-pneumatic spring is built and the accuracy of the model is verified. The one-quarter model of vehicle is built up to compare the simulations of vehicles equipped with single-accumulator hydro-pneumatic suspension to that with two-stage pressure hydro-pneumatic suspension under random road excitation. The results show that, compared with the former situation, the body acceleration root mean square (RMS) value under the latter condition decreases by 20.1% (paved road/full-load) and 10.7% (off-road/empty-load) respectively, and the tire dynamic load decreases by 36.8% and 10.4% respectively. Although the RMS value of dynamic travel of suspension increases by 11.8% and 1.9%, the probability of collision to set blocks is less than 0.1%, remaining in a reasonable range.

    Keywords off-road vehicle; two-stage pressure hydro-pneumatic; ride comfort; modeling and simulation

    Hybrid Modeling Input and Output for a Rhombic Micro-displacement Amplifier

    ZHANGChunlin,HEGuojing,YIJing
    (College of Civil Engineering and Mechanics, Central South University of Forestry and Technology Changsha, 410004, China)

    Abstract This paper proposes a rhombic micro-displacement amplifier (RMDA) for piezoelectric actuator. The amplifier consists of three parts: piezoelectric actuator (PA), flexible hinge and rhombic amplification part. However, the key drive part of RMDA is PA. PA has hysteresis effects, it is mean that the RMDA also has. To ensure the accuracy for amplifiers, the hybrid model which is based on discrete classical Preisach model (CPM) and support vector machine (SVM), is introduced to accounts for hysteresis effects. Finally, the hybrid model and its inversion model of the amplifier are modeled and successfully validated via experimental results. The results demonstrate that hybrid modeling have high accuracy for express the input/output relations of RMDA which have hysteresis effects.

    Keywords displacement amplifier; piezoelectric actuator (PA); hysteresis effects; hybrid modeling

    Acoustic Emission Mechanism and Testing of Leakage in Natural Gas Pipeline Ball Valve

    LIZhenlin1,2,ZHANGHaifeng3,HAOYibo4,ZHANGNing1,2,LEIHongxiang1,2,CHENXin1,2,LIUZhichao1,2
    (1.School of Mechanical and Transportation Engineering, China University of Petroleum-Beijing Beijing, 102249, China)

    (2.Beijing Key Laboratory of Process Fluid Filtration and Separation Beijing,102249, China)

    (3.PetroChina Pipeline R & D Center Langfang, 065000, China)

    (4.PetroChina Eastern Pipeline Co., Ltd Shanghai, 200120, China)

    Abstract The jet gas from the internal leakage in a ball valve, which is an indispensable equipment in a high-pressure natural gas transmission pipeline, can produce acoustic emission (AE) signals. Thus, the research on the characteristic rule of the AE signals will be helpful for the quantitative detection of the valve leakage flow. In dealing with this problem, the AE mechanism and testing of internal leakage in a ball valve are investigated. An AE testing system is used to measure the leakage in ball valves of three different sizes. In the analysis process of the testing data, the spectrum distribution characteristics of the AE signals on leakage levels are described, and the characteristic parameters of signals (information entropy, root mean square (RMS) and frequency peak) are extracted using the wavelet packet analysis method. On this basis, fitted curves about the characteristic parameters and the leakage flow are established, and the correlation is evaluated using the R-square (determination coefficient). The evaluation results show that the RMS value with the highest correlation (R-square =0.979) can be used for the quantitative detection of leakage in a ball valve in the natural gas pipeline.

    Keywords natural gas pipeline; ball valve; acoustic emission mechanism; testing experiment

    Sensitivity Analysis of Wheel Quality and Location on Rotor Critical Speed

    PANHonggang1,2,YUANHuiqun1,ZHAOTianyu1,YANGWenjun1
    (1.School of Mechanical Engineering & Automation, Northeastern University Shenyang, 110819, China)

    (2.School of Energy and Power Engineering, Shenyang Institute of Engineering Shenyang, 110136, China)

    Abstract To improve the efficiency of thermal power unit and reduce the pollutant emissions, the steam turbine unit in series of transformation on flow is changed, the quality and center position of mass of the rotor is changed, which leads to the change of rotor critical speed and a certain effects on the run. In this paper, using the experiment to measure the different wheel quality and position of the rotor critical speed, the results are compared with the theoretical calculation result analysis, and the sensitivity analysis methods is introduced to analyze the influence of the quality of wheel and position on changing of steam turbine rotor critical speed. Experimental study agrees well with the theoretical analysis: the rotor critical speed is decreased with the increasing of the wheel quality and is increased with the increasing of the wheel offset. The sensitivity coefficient of wheel offset on the critical speed is between 0.25 and 2.4, the greater the offset, the higher the sensitivity coefficient; the sensitivity coefficient of the quality on the critical speed is between -0.35 and -0.001, the greater the increasing ratio of the quality, the smaller the sensitivity coefficient; the influence of the offset on the critical speed is far greater than the influence of the quality, which is about 7~10 times. When the quality offset of the same wheel is greater than 40%, and the rotor critical speed changes obviously. In the same offset location, the increasing of the mass is less than 50%,, and the rotor critical speed changes obviously. The method proposed can be used to solve the problem of the adjustment of the steam turbine rotor critical speed.

    Keywords critical speed; sensitivity; offset; Bode diagram

    The Calculating Method and Application of Evidence Weight Based on Fault Sensitivity

    HUJinhai1,2,GAOXingwei1,ZHANGYu1,RENLitong1,PENGJingbo1
    (1.Aeronautics and Astronautics Engineering Institute, Air Force Engineering University Xi′an, 710038, China)

    (2.Co-Innovation Center for Advanced Aero-Engine Beijing, 100191, China)

    Abstract With a focus on the problem of the typical evidence weight calculating method that it cannot define the weight and complete the decision fusion in the high conflict evidence D-S fusion problem where only a few sensors make the judgments correctly while the most incorrectly, the paper proposes an evidence weight calculating method based on fault sensitivity. Firstly, kernel principal component analysis (KPCA) is used to obtain the nonlinear sensitive feature; then, calculating the sensitivity based on such feature and obtaining sensor decision weight based the fault detection sensitivity; applying the weight obtained above and that obtained from equal weight method and decision-making distance based method to the fusion diagnosis of the simulated fault in the rotor where three sensor are installed. The results show that: the weight obtained through the method proposed in the paper can reflect the sensitivity of different sensors when detecting the faults. High weight is given to the sensors which contain much fault information and sensitive to fault and low weight will be given to the sensors of the opposite kind. In the method proposed, evidence weights play the role as a “regulator”, which makes it possible to obtain better decision fusion result whether in the cases where only a few sensors manage to find the fault and give the correct diagnosis or in the cases where few or no conflict exists.

    Keywords multi-sensor information fusion; D-S evidence theory; fault diagnosis; fault sensitivity; evidence weight

    Optimization of Support Vector Machine and Its Application in Intelligent Fault Diagnosis

    WANGBaojian1,ZHANGXiaoli2,FUYANGAoxiao1,CHENXuefeng1
    (1.State Key Laboratory for Manufacturing and Systems Engineering, Xi′an Jiaotong University Xi′an, 710049, China)

    (2.Key Laboratory of Road Construction Technology and Equipment, Ministry of Education, Chang′an University Xi′an, 710064, China)

    Abstract Single support vector machine has low precision in fault diagnosis of bearing and gear system, the sample feature extraction method of support vector machine and the method of parameter optimization of support vector machine are studied to improve the accuracy of the support vector machine in the fault diagnosis of bearing gear system. The input samples of support vector machines are constructed by the kernel principal component analysis to reduce data redundancy, extract high dimension information of the data, then particle swarm optimization algorithm is used to optimize the kernel function parameter and penalty factor of SVM, finally, the optimized support vector machine model is used for fault diagnosis. A comparative experiment on the fault diagnosis of bearing gear is carried out in order to verify the effectiveness of the proposed method, the results show that the proposed method improves the diagnostic accuracy significantly in comparison with the general support vector machines,the effectiveness and advantages of the intelligent diagnosis method are verified.

    Keywords support vector machine; kernel principal component; particle swarm optimization algorithm; fault diagnosis

    Dynamic Testing and Modeling of a Magnetorheological Damper

    MEIZhen1,2,GAOYichao1,2,GUOZixiong1,2
    (1.College of Civil Engineering, Huaqiao University Xiamen, 361021, China)

    (2.Key Laboratory of Structural Engineering and Disaster Prevention of Fujian Province Xiamen, 361021, China)

    Abstract Establishing precise mechanical models of magnetorheological dampers is an important prerequisite not only for the response analysis and design of structures with the dampers, but also for obtaining good vibration control effect. In this paper, the dynamic performance test of a magnetorheological damper with the nominated maximum damping force of 10 kN is first carried out. Based on the test results, parametric and non-parametric dynamic models of the damper are developed and the effectiveness of the proposed models is verified. Finally, the results of the two different modeling approaches are compared. Results show that the formulated hyperbolic tangent hysteresis model (parametric model) could well describe the dynamic behavior of the magnetorheological damper. Besides, both the forward and reverse back propagation (BP) neural network models (non-parametric models) have good performance in fitting training data, generalization ability and noise immunity. Furthermore, the BP neural network models fit with a higher accuracy than the hyperbolic tangent hysteresis models, while the later enjoy a simpler expression of the damping force which is relatively easier to implement in software.

    Keywords magnetorheological damper; dynamic performance test; hyperbolic tangent hysteresis model; BP neural network model

    Multidimensional Performance Limit States for Fragility Analyze of Plane Irregular Structure

    HUANGXiaoning1,DUYongfeng1,2,LIHui1,2
    (1.Institute of Earthquake Protection and Disaster Mitigation, Lanzhou University of Technology Lanzhou, 730050, China)

    (2.Western Center of Disaster Mitigation in Civil Engineering of Ministry of Education, Lanzhou University of Technology Lanzhou, 730050, China)

    Abstract Based on the introduction of the critical angle, a new method for multidimensional performance limit states is developed to analyze the structural fragility of plane irregular frame-shear-wall structures. First, the critical angle of seismic wave is determined using the wavelet transform in Matlab. Then, concerning the dependency of the limit states, the probability of exceedance can be calculated by adopting the inter-story drift and inter-story torsion angle as the quantitative indexes. The fragility curves are achieved under two dimensional performance limit states. Using this method to analyze the aforementioned structure, the fragility curves under four performance levels are obtained, which are normal operation, immediate occupancy, life safety, and collapse prevention. The results indicated that the critical angle of seismic wave has a significant impact on the performance of the plane irregular structures. In order to avoid overestimating seismic performance of plane irregular structures, inter-story drift and inter-story torsion angle should be simultaneously considered. Therefore, for plane irregular estimate, the analysis method of fragility under multidimensional performance limit states is more safe and reliable.

    Keywords plane irregular structures; fragility analyze; multidimensional performance limit states; inter-story drift; inter-story torsion angle

    Investigation on Gas Effect of Displacer for Free-Piston Stirling Engine

    LIWei1,2,MOUJian1,2,HONGGuotong1
    (1.Technical Institute of Physics and Chemistry, Chinese Academy of Sciences Beijing, 100190, China)

    (2.University of Chinese Academy of Sciences Beijing, 100049, China)

    Abstract In order to have a deep understanding of the kinematical characteristics of the displacer of free-piston stirling engine, the effects of gas force, which is driven by the displacer only, is analyzed. And the running experiment of displacer individually is designed. During the course of the experiment a new calculation formula for natural frequency of displacer is proposed. In this paper, the gas force acting on the displacer is analyzed in the case of leading or lagging displacement of piston. And the acting mechanism of gas force is stated by means of rotation vector decomposition. The results show that a part of gas force is regarded as gas spring. The gas spring and the mechanical spring are parallelly arranged and the natural frequency of the system increases in the case of the pressure wave in advance or delaying the displacement of piston less than 90°. But the gas force is regarded as gas inertia force and the natural frequency of the system decreases in the case of the pressure wave in advance or delaying the displacement of piston greater than 90° but less than 180°. And the higher the temperature of the heat source, the smaller the natural frequency of the system. The accuracy of the effects of gas force and calculation formula for natural frequency are verified by an experimental text on a free-piston stirling engine which is designed by our laboratory.

    Keywords stirling engine; free-piston; displacer; the effects of gas; natural frequency

    Comprehensive Position Error of Delta Robot and Analysis of Its Coupling Characteristics

    ZHENGKunming1,2,ZHANGQiuju1,2
    (1.School of Mechanical Engineering, Jiangnan University Wuxi, 214122, China)

    (2.Jiangsu Key Laboratory of Advanced Food Manufacturing Equipment & Technology Wuxi, 214122, China)

    Abstract Taking Delta robot as the object of analysis, studying the position error model of the moving platform, and the coupling characteristic of the error source is analyzed. Using the position characteristic of driven arms, based on the method of geometry space vector, the mechanism error model of Delta robot is established; Using the principle of mathematical statistics and space vector as the basic, the joint clearance error model of Delta robot is derived; Based on finite element theory, on the basic of elastic dynamic model of Delta robot, the flexible error model is eatablished; Then, considering the three kinds of error source, the comprehensive position error model of Delta robot is established. Finally, by use of simulation software of Adams and Workbench, numerical calculation of Matlab and field experiment of FARO laser tracker, the correctness of comprehensive position error model is verified, and coupling characteristics of the error sources are analyzed, the relationship between position error and the axis's azimuth of system coordinate is elaborated. The results show that the error sources that affect the position error of the Delta robot's moving platform are not simply superimposed, but with obvious coupling properties, and the directional position error of the moving platform changes with the azimuth of the coordinate axis.

    Keywords Delta robot; mechanism error model; clearance error model; flexible error model; comprehensive position error; coupling characteristics

    Investigation on the Blade Tip Clearance Monitoring of Turbomachinery Based on the Pulse-Trigger Method

    WANGWeimin,SHAOHuajin,CHENLifang,QUWei
    (School of Mechanical and Electrical Engineering, Beijing University of Chemical Technology Beijing, 100029, China)

    Abstract Blade tip clearance (BTC) and blade vibration measuring are the principal tools for blade health monitoring (BHM), which can improve the efficiency and reliability of turbomachinery. In this paper, a BTC monitoring method based on pulse-trigger of eddy current sensor (ECS) is presented, in the light of the sub-Nyquist sampling caused by the limit of narrow-band, especially for the high linear velocity. The proposed method optimizes the static radial and circumferential calibration technology to obtain the sensitivity of ECS in the different relative location against the tip of blade and to obtain the function of impedance of the ECS that relates to the distance and the contact ratio of the ECS from the tip of blade. Combining blade tip timing (BTT) technology into BTC monitoring, this method calculated the contact ratio of the ECS from the blade, in order to obtain the accurate clearance depends on the function obtained before for all sampling points. Also, a BTC monitoring rig is established to validate the feasibility of this method under varying rotating speed. The results show that this method can solve the problem of sub-Nyquist sampling caused by the limit of narrow-band for ECS in the condition of high linear velocity. Compared with the location obtained by peak picking, the method presented in this paper can significantly improve the accuracy of BTC monitoring. This research can provide the basis for BTC active control and BHM based on ECS.

    Keywords turbine blades; blade tip clearance(BTC) monitoring; pulse-trigger; optimize static calibration technology; eddy current sensor (ECS)

    Fault Diagnosis and Simulation of Aircraft Air Conditioning System Based on FMEA

    LIBingyue1,SUNJianhong1,LIUHaigang2,SUNZhi1,CHENQiang1
    (1.College of Aerospace Engineering, Nanjing University of Aeronautics & Astronautics Nanjing, 210016, China)

    (2.Shenyang Aircraft Design and Research Institute Shenyang, 110035, China)

    Abstract In this paper, in order to ensure the safety and reliability of the aircraft conditioning system, the failure mode and effects analysis(FMEA) of the main parts of aircraft air conditioning system is carried out. The whole air conditioning is simulated by Matlab/Simulink, and the failure diagnosis criteria of the key parts of an air conditioning system are regulated. Then, the abnormal cabin temperature of the aircraft is taken as an example to demonstrate the fault tree analysis method (FTA). The results show that the failure of the turbine has the largest influence on the abnormal cabin temperature, which should be took into serious consideration when it comes to the work of design and maintenance.

    Keywords fault simulation; fault tree analysis; failure mode and effects analysis (FMEA); aircraft air conditioning system; severity

    Numerical Simulation and Study on Flow-Induced Noise During Automotive Throttle Quick-Opening Process

    YANGShuai1,2,WEIYanan1,2,XUELiangjun3,CHANGGuofeng1,2
    (1.New Energy Automotive Engineering Center, Tongji University Shanghai, 201804, China)

    (2.School of Automotive Studies, Tongji University Shanghai, 201804, China)

    (3.Qishuyan Locomotive Co., Ltd Qishuyan, 213011, China)

    Abstract The influence of the throttle motion on the air flow-induced noise process is studied. By applying the computational fluid dynamics (CFD) and computational aeroacoustics (CAA) coupling approach and the moving mesh technique, a 3-dimensional air transient flow and flow-induced noise process is simulated to realize the throttle quick-opening from the closed position to the fully open position. The principle of the air flow and flow-induced noise is analyzed when the throttle is opening at different rotation angle. At the beginning period of the throttle rotation, the vortex appears in the flow field behind the throttle, the pressure drops distinctly between the front and the back of the throttle, two flow-induced noise fields appear near the top and bottom position of the throttle and move towards the downstream and then merge gradually, the maximum value of the acoustic power level appears when the throttle rotation angle is opening near 40 degree, the flow-induced noise acoustic power level firstly increases and then decreases. With the opening angle increasing, the vortex is weaken and the value of pressure drop decreases between the two sides of the throttle. The time-frequency domain numerical simulation results fo the flow-induced noise show that during the throttle quick-opening process the flow-induced noise is mainly a kind of middle and low frequency broadband noise, the main part of the noise is frequency noise below about 100 Hz, the acoustic pressure is higher and does not reduce significantly with the increase of the distance between the measurement point and the throttle center. Reducing the low frequency noise plays a key role in the flow-induced noise control of the throttle.

    Keywords throttle; flow-induced noise; transient flow; numerical simulation

    Limit Cycle Oscillation of Subspan Oscillation for Two Bundled Conductors by Using Incremental Harmonic Balance Method

    YUYangyang1,2,3,GUOHulun1,3,CAOShuqian1,3,LIUBin4,CHENYushu1,3,5
    (1.Department of Mechanics, Tianjin University Tianjin, 300072, China)

    (2.Tianjin University Renai College Tianjin, 301636, China)

    (3.Tianjin Key Laboratory of Nonlinear Dynamics and Chaos Control Tianjin, 300072, China)

    (4.China Electric Power Research Institute Beijing, 100192, China)

    (5.School of Astronautics, Harbin Institute of Technology Harbin, 150001, China)

    Abstract Subspan oscillation on bundled conductors is one of the most important failures of overhead transmission lines, which is due to the leeward conductor lying in the wake of the windward conductor. In this paper, considering the aerodynamic nonlinearities, a 2-DOF subspan oscillation dynamic model of the leeward conductor lying in the wake of the windward conductor is given firstly. And then, the incremental harmonic balance method (IHB) is used to derive high-order limit cycle oscillation of subspan oscillation system, and the first to the third harmonic responses of the limit cycle oscillation are obtained. Results show that subspan oscillation only exists in a certain range of wind velocity. Moreover, the higher order of the harmonic the weaker impact it is, and the first harmonic response agree with numerical results of Runge-Kutta well. Lastly, the influence of span and original position of the leeward conductor on subspan oscillations has been analyzed, which can provide the technical support for suppression or prevention subspan oscillation.

    Keywords bundled conductors; subspan oscillation; limit cycle oscillation; incremental harmonic balance method

    Fault Diagnosis Under Variable Conditions Based on Parameter Optimized Variational Mode Decomposition and Envelope Order Spectrum

    JIANGZhanwei,ZHENGJinde,PANHaiyang,PANZiwei
    (School of Mechanical Engineering, Anhui University of Technology Maanshan, 243032, China)

    Abstract Based on the parameter optimized variational mode decomposition (POVMD) and envelope order spectrum, a new fault diagnosis method is proposed to extract the fault features of rolling bearing under variable speed condition. First, the vibration signal of rolling bearing is decomposed into several intrinsic mode functions (IMFs) by POVMD. Second, each IMF is resampled in the angular domain and transformed into stationary signals. Then, the Hilbert transform is used to estimate the envelope function of the resampled signals. The obtained envelope functionsare are analyzed using the order tracking technology and the fault feature information is read from the order spectrum. The comparisons show the superiority of POVMD over empirical mode decomposition by analyzing the simulation signals. Finally, the proposed fault diagnosis method for rolling bearing with variable speed is applied on the experimental data analysis and the results show that the propose method can effectively achieve the fault diagnosis of rolling bearing in variable speed.

    Keywords variational mode decomposition; variable conditions; envelope order spectrum; rolling bearing; fault diagnosis

    Study on the Stability of Two Dimensional Ultrasonic Vibration Milling System in Vertical Plane

    ZHAOBo,ZHAOBinbin,FANKaiyang,ZHANGYuemin
    (School of Mechanical and Power Engineering, Henan Polytechnic University Jiaozuo, 454000, China)

    Abstract In order to study the stability of the two dimension ultrasonic vibration milling system in vertical plane, the stability model of a two-dimensional ultrasonic vibration milling in vertical plane is established. The stability of the two dimensional ultrasonic vibration milling is studied using the full discrete method. The stability of the system is obtained by numerical simulation using Matlab software. The results of chatter test for milling titanium alloy material show that the numerical simulation results are in good agreement with the experimental results when the speed was 1 000~3 500 r/min. The correctness of the stability model of two dimensional ultrasonic vibration milling in vertical plane is verified. The stability of the system can be improved by applying two dimensional ultrasonic vibration in the vertical plane. The axial depth limit maximum of cut is increased by about 13.6%. The research fills up the blank of the two-dimensional ultrasonic vibration assisted milling stability study.

    Keywords two dimensional ultrasonic vibration milling; full-discretization method; stability; lobe diagram

    The vibration monitoring method is the main approach of condition monitoring for wind turbine transmission system. The vibration monitoring strategy and the process of each component′s vibration feature extraction of wind turbine transmission system are analyzed firstly, and the trend indicators of vibration feature such as sideband power factor (SBPF) and sideband energy ratio (SER) are introduced especially. Then, it is pointed out that the key to retrieve the huge economic losses caused by the fault of wind turbine transmission system is incipient fault prediction for wind turbine transmission system, and the generalization manifold learning-based incipient fault prediction method for wind turbine transmission system is introduced especially. Finally, the function and characters of the existing vibration monitoring system for wind turbine transmission system are analyzed including system architecture, data acquisition configuration and monitoring analysis methods. Furthermore, it is pointed out that big data-based predictive analytics and intelligent maintenance based on multi-source information fusion technique would be the important development trend of health management for wind turbine transmission system.

    wind turbine transmission system; vibration monitoring; feature extraction; incipient fault prognosis; big data

    10.16450/j.cnki.issn.1004-6801.2017.03.030

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51475148)

    2016-06-29;

    2016-11-05

    TB115; TB552; TH113

    趙波,男,1956年8月生,博士、教授、博士生導(dǎo)師。主要研究方向?yàn)橄冗M(jìn)制造技術(shù)、硬脆材料精密加工理論與裝備等。曾發(fā)表《分離型縱-扭復(fù)合超聲銑削的穩(wěn)定性分析》(《兵工學(xué)報(bào)》2015年第36卷第7期)等論文。 E-mail:zhaob@hpu.edu.cn

    Research Progress of Vibration Monitoring for Wind Turbine Transmission System
    TANGBaoping,LUOLei,DENGLei,HANYan(The State Key Laboratory of Mechanical Transmission, Chongqing University Chongqing, 400030, China)

    猜你喜歡
    葉瓣切削力刀具
    我家的熊掌
    多肉植物盆栽
    多肉植物盆栽
    無織構(gòu)刀具與織構(gòu)刀具銑削性能對(duì)比研究
    30CrMnSiNi2A插銑加工的切削力研究
    圓刀片切削力計(jì)算方法
    切削刀具刃口形貌對(duì)刀具使用壽命的影響
    多功能刀具
    圣母的玫瑰
    文學(xué)港(2016年7期)2016-07-06 15:48:50
    山高刀具推陳出新迎馬年——山高刀具2013秋季新聞發(fā)布會(huì)
    亚洲丝袜综合中文字幕| 亚洲欧美日韩东京热| 国产高清有码在线观看视频| 日本-黄色视频高清免费观看| 欧美人与善性xxx| 亚洲av中文字字幕乱码综合| 国产又色又爽无遮挡免| 国产精品99久久99久久久不卡 | 欧美激情极品国产一区二区三区 | 国产精品麻豆人妻色哟哟久久| 深爱激情五月婷婷| 国产黄频视频在线观看| 国产伦精品一区二区三区四那| 99久国产av精品国产电影| 久久精品人妻少妇| 80岁老熟妇乱子伦牲交| 下体分泌物呈黄色| 97在线视频观看| 国产永久视频网站| 亚洲人成网站在线观看播放| 国产成人aa在线观看| 一级毛片电影观看| 大陆偷拍与自拍| 中文在线观看免费www的网站| 99久久精品国产国产毛片| 亚州av有码| 男女边摸边吃奶| 国产成人freesex在线| 亚洲美女搞黄在线观看| 中文字幕精品免费在线观看视频 | 2018国产大陆天天弄谢| 午夜日本视频在线| 搡女人真爽免费视频火全软件| 国产一区二区在线观看日韩| 精品久久久久久久久亚洲| 亚洲经典国产精华液单| 一本久久精品| 狂野欧美激情性bbbbbb| 亚洲精品中文字幕在线视频 | 能在线免费看毛片的网站| av不卡在线播放| 欧美xxxx黑人xx丫x性爽| 欧美变态另类bdsm刘玥| 色婷婷久久久亚洲欧美| 91久久精品国产一区二区三区| 亚洲精品久久久久久婷婷小说| 熟女av电影| 丝瓜视频免费看黄片| 欧美日韩国产mv在线观看视频 | 色视频www国产| 黄色怎么调成土黄色| 日日撸夜夜添| 在线观看免费日韩欧美大片 | 久久韩国三级中文字幕| 免费播放大片免费观看视频在线观看| 在线看a的网站| 日本欧美视频一区| 99热6这里只有精品| 亚洲国产精品专区欧美| 亚洲精品自拍成人| 五月开心婷婷网| 免费人成在线观看视频色| 91精品国产九色| 久久人人爽av亚洲精品天堂 | 欧美丝袜亚洲另类| 久久午夜福利片| 亚洲av福利一区| 日本av免费视频播放| 精品国产一区二区三区久久久樱花 | 男女下面进入的视频免费午夜| 交换朋友夫妻互换小说| 成人高潮视频无遮挡免费网站| 久久精品国产鲁丝片午夜精品| 狂野欧美激情性bbbbbb| 18禁在线无遮挡免费观看视频| 久久ye,这里只有精品| 亚洲怡红院男人天堂| 国产精品国产三级专区第一集| 亚洲av二区三区四区| 免费看不卡的av| 国产av码专区亚洲av| av在线蜜桃| 黄色一级大片看看| 身体一侧抽搐| 中文欧美无线码| 少妇丰满av| av一本久久久久| 国产美女午夜福利| 亚洲精品一区蜜桃| 久久精品国产a三级三级三级| 人妻少妇偷人精品九色| 国产男女内射视频| 少妇猛男粗大的猛烈进出视频| 国产精品一区二区在线不卡| 国产综合精华液| 在线观看免费高清a一片| 老司机影院成人| videossex国产| 91精品国产国语对白视频| 综合色丁香网| 精品一区在线观看国产| 成人国产av品久久久| 国产一区二区三区综合在线观看 | 亚洲成人av在线免费| 欧美老熟妇乱子伦牲交| a级一级毛片免费在线观看| 国产精品麻豆人妻色哟哟久久| 国产亚洲av片在线观看秒播厂| 熟妇人妻不卡中文字幕| 欧美日本视频| 六月丁香七月| 青春草国产在线视频| 欧美成人精品欧美一级黄| 联通29元200g的流量卡| 国产欧美日韩一区二区三区在线 | 国产亚洲5aaaaa淫片| 久久久亚洲精品成人影院| 欧美三级亚洲精品| 亚洲色图av天堂| 人妻 亚洲 视频| 免费人成在线观看视频色| 男女无遮挡免费网站观看| 精品亚洲乱码少妇综合久久| 99热这里只有精品一区| 国产精品蜜桃在线观看| 青春草国产在线视频| 狠狠精品人妻久久久久久综合| 又黄又爽又刺激的免费视频.| 国产熟女欧美一区二区| www.色视频.com| 777米奇影视久久| 亚洲av免费高清在线观看| 在线 av 中文字幕| 黄色欧美视频在线观看| 夜夜看夜夜爽夜夜摸| 日韩中文字幕视频在线看片 | 欧美区成人在线视频| 国产精品国产三级专区第一集| 欧美一区二区亚洲| 热re99久久精品国产66热6| kizo精华| 午夜福利网站1000一区二区三区| 欧美日韩视频精品一区| 国产精品爽爽va在线观看网站| 亚洲欧洲日产国产| 丰满少妇做爰视频| 三级国产精品片| 老师上课跳d突然被开到最大视频| 日韩在线高清观看一区二区三区| av天堂中文字幕网| 亚洲激情五月婷婷啪啪| 国产成人精品婷婷| 亚洲国产精品一区三区| 国产精品一区二区在线不卡| 一二三四中文在线观看免费高清| 久久午夜福利片| 免费观看性生交大片5| 成人18禁高潮啪啪吃奶动态图 | 五月天丁香电影| av卡一久久| 免费观看性生交大片5| 啦啦啦中文免费视频观看日本| 精品少妇久久久久久888优播| 中国国产av一级| 成人漫画全彩无遮挡| 亚洲精品国产成人久久av| 精品一区在线观看国产| 精品一区在线观看国产| 国产视频首页在线观看| 男女国产视频网站| 我要看日韩黄色一级片| 亚洲精品国产成人久久av| 亚洲国产欧美人成| 国产亚洲91精品色在线| 国产国拍精品亚洲av在线观看| 日日撸夜夜添| 日韩三级伦理在线观看| 五月玫瑰六月丁香| 久久久成人免费电影| 亚洲欧美成人综合另类久久久| 亚洲精品日韩av片在线观看| 亚洲天堂av无毛| 亚洲精品日韩av片在线观看| 黄片wwwwww| 免费播放大片免费观看视频在线观看| 日韩不卡一区二区三区视频在线| 99久久精品热视频| 又大又黄又爽视频免费| 精品国产三级普通话版| 日韩不卡一区二区三区视频在线| 亚洲美女视频黄频| 99久久精品热视频| 18禁裸乳无遮挡动漫免费视频| 国产精品国产三级国产av玫瑰| 色吧在线观看| 97超视频在线观看视频| 五月玫瑰六月丁香| 欧美国产精品一级二级三级 | 黑人猛操日本美女一级片| 亚洲欧洲日产国产| 成人18禁高潮啪啪吃奶动态图 | 人妻系列 视频| 91精品国产九色| 亚洲图色成人| 国产日韩欧美在线精品| 精品久久久噜噜| 久久婷婷青草| 网址你懂的国产日韩在线| 国产男女超爽视频在线观看| av.在线天堂| 天天躁夜夜躁狠狠久久av| 免费高清在线观看视频在线观看| 日本av免费视频播放| 搡女人真爽免费视频火全软件| 国产伦精品一区二区三区视频9| 国产精品一及| 啦啦啦在线观看免费高清www| 少妇的逼好多水| 成人漫画全彩无遮挡| 高清不卡的av网站| 王馨瑶露胸无遮挡在线观看| 男人和女人高潮做爰伦理| av女优亚洲男人天堂| 天堂俺去俺来也www色官网| 中文字幕制服av| 尤物成人国产欧美一区二区三区| 极品少妇高潮喷水抽搐| 蜜桃久久精品国产亚洲av| 精品人妻一区二区三区麻豆| 国产精品麻豆人妻色哟哟久久| 成人午夜精彩视频在线观看| 一个人看视频在线观看www免费| 少妇高潮的动态图| 午夜福利在线观看免费完整高清在| 亚洲精品乱码久久久久久按摩| 亚洲精品,欧美精品| 日本av免费视频播放| 美女国产视频在线观看| 各种免费的搞黄视频| 日韩强制内射视频| 国产黄片视频在线免费观看| 久久99热这里只频精品6学生| 久久久久国产精品人妻一区二区| 99热这里只有是精品50| 亚洲欧美精品自产自拍| 国产男女内射视频| 国产中年淑女户外野战色| 少妇的逼水好多| 肉色欧美久久久久久久蜜桃| 亚洲精品亚洲一区二区| 色婷婷av一区二区三区视频| 国产精品一区二区在线不卡| 国产精品国产三级国产专区5o| 日韩视频在线欧美| 在线精品无人区一区二区三 | 中国国产av一级| 好男人视频免费观看在线| 久久久久久久精品精品| 午夜福利影视在线免费观看| 妹子高潮喷水视频| 精品99又大又爽又粗少妇毛片| 久热这里只有精品99| 精品酒店卫生间| 男女啪啪激烈高潮av片| 精华霜和精华液先用哪个| 免费久久久久久久精品成人欧美视频 | 人妻制服诱惑在线中文字幕| 欧美xxⅹ黑人| 国内揄拍国产精品人妻在线| 日韩强制内射视频| 日韩 亚洲 欧美在线| 亚洲精品中文字幕在线视频 | 男人和女人高潮做爰伦理| 干丝袜人妻中文字幕| 久久国产精品男人的天堂亚洲 | 91久久精品电影网| 亚洲伊人久久精品综合| 少妇的逼好多水| av国产久精品久网站免费入址| 免费观看无遮挡的男女| 日本黄色片子视频| 精品视频人人做人人爽| 国产白丝娇喘喷水9色精品| 边亲边吃奶的免费视频| 国产精品成人在线| 啦啦啦在线观看免费高清www| 国产成人a∨麻豆精品| 免费在线观看成人毛片| 99久久人妻综合| 91狼人影院| 中文字幕av成人在线电影| 97超视频在线观看视频| 女人十人毛片免费观看3o分钟| 久久人人爽人人片av| 不卡视频在线观看欧美| 亚洲精品自拍成人| 欧美精品国产亚洲| 性色avwww在线观看| 国产69精品久久久久777片| 国产亚洲av片在线观看秒播厂| 国产在线免费精品| 欧美精品一区二区免费开放| 97在线视频观看| 女性被躁到高潮视频| 中文在线观看免费www的网站| 国产在线免费精品| 国产黄频视频在线观看| 欧美另类一区| 午夜福利影视在线免费观看| 国产精品国产av在线观看| 麻豆成人av视频| 五月开心婷婷网| 成人毛片a级毛片在线播放| 18+在线观看网站| 蜜桃久久精品国产亚洲av| 国产 一区 欧美 日韩| 色网站视频免费| 日韩不卡一区二区三区视频在线| 男男h啪啪无遮挡| 极品教师在线视频| 国产精品不卡视频一区二区| 色综合色国产| 亚洲精品日韩av片在线观看| 伦理电影免费视频| 下体分泌物呈黄色| 赤兔流量卡办理| 夜夜看夜夜爽夜夜摸| 五月玫瑰六月丁香| 春色校园在线视频观看| 99热6这里只有精品| 亚洲人成网站在线播| 联通29元200g的流量卡| 久久国产精品大桥未久av | 成人毛片60女人毛片免费| 视频区图区小说| 成人漫画全彩无遮挡| 乱系列少妇在线播放| 亚洲丝袜综合中文字幕| 纵有疾风起免费观看全集完整版| 亚洲av男天堂| 大片免费播放器 马上看| 秋霞在线观看毛片| 三级经典国产精品| 国产极品天堂在线| 大码成人一级视频| 好男人视频免费观看在线| 国产精品久久久久久久久免| 日韩一区二区视频免费看| 免费观看在线日韩| 婷婷色综合大香蕉| 国产视频首页在线观看| 五月天丁香电影| 少妇精品久久久久久久| 99精国产麻豆久久婷婷| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av电影在线观看一区二区三区| 最近的中文字幕免费完整| av专区在线播放| 国产成人精品久久久久久| 日韩大片免费观看网站| 91精品伊人久久大香线蕉| 欧美日韩亚洲高清精品| 久久精品国产亚洲av涩爱| 国产黄色免费在线视频| 久久久精品94久久精品| 国产成人a∨麻豆精品| 国产精品久久久久久精品古装| 日日摸夜夜添夜夜添av毛片| 精品国产露脸久久av麻豆| 日本wwww免费看| 天天躁日日操中文字幕| 97精品久久久久久久久久精品| 成年女人在线观看亚洲视频| 免费看光身美女| 精品人妻视频免费看| 99九九线精品视频在线观看视频| 欧美精品亚洲一区二区| 精品视频人人做人人爽| 久久热精品热| 成人特级av手机在线观看| a级毛色黄片| 国产在线视频一区二区| 搡女人真爽免费视频火全软件| 国产免费视频播放在线视频| 亚洲不卡免费看| 蜜臀久久99精品久久宅男| 综合色丁香网| 国产色婷婷99| 哪个播放器可以免费观看大片| 国产精品一区二区在线不卡| 国产日韩欧美亚洲二区| 日日啪夜夜撸| 大片免费播放器 马上看| 午夜视频国产福利| 亚洲av成人精品一二三区| 国产精品99久久99久久久不卡 | 少妇的逼好多水| 麻豆精品久久久久久蜜桃| 日本黄色日本黄色录像| 亚洲一级一片aⅴ在线观看| 一级黄片播放器| 性高湖久久久久久久久免费观看| 久久久久视频综合| 一级黄片播放器| 久久久久久久久久人人人人人人| 97在线视频观看| 欧美高清性xxxxhd video| 成人毛片a级毛片在线播放| 嫩草影院新地址| 精品99又大又爽又粗少妇毛片| 日本黄大片高清| 国产精品女同一区二区软件| 亚洲不卡免费看| 交换朋友夫妻互换小说| 久久99热这里只频精品6学生| 国产亚洲av片在线观看秒播厂| 日韩一本色道免费dvd| 久久影院123| 晚上一个人看的免费电影| 午夜福利在线观看免费完整高清在| 美女xxoo啪啪120秒动态图| 日本欧美视频一区| 日韩精品有码人妻一区| 卡戴珊不雅视频在线播放| 青春草国产在线视频| 亚洲av二区三区四区| 男女下面进入的视频免费午夜| 国产精品嫩草影院av在线观看| 赤兔流量卡办理| 精品国产三级普通话版| 色婷婷久久久亚洲欧美| 久久久久国产网址| 综合色丁香网| 国产日韩欧美在线精品| 肉色欧美久久久久久久蜜桃| 国产国拍精品亚洲av在线观看| 五月天丁香电影| 十八禁网站网址无遮挡 | 一区在线观看完整版| 久久 成人 亚洲| 在线观看人妻少妇| 亚洲精品久久久久久婷婷小说| 99久久综合免费| 午夜福利网站1000一区二区三区| 搡老乐熟女国产| 日日啪夜夜爽| 国产黄色免费在线视频| 99久久综合免费| 国产av一区二区精品久久 | 色视频在线一区二区三区| 老女人水多毛片| 在线观看av片永久免费下载| 久久ye,这里只有精品| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久精品精品| 久久久色成人| 日韩一区二区视频免费看| 色婷婷av一区二区三区视频| 极品教师在线视频| 男人舔奶头视频| 最近手机中文字幕大全| 成人亚洲精品一区在线观看 | 欧美日韩在线观看h| 春色校园在线视频观看| av福利片在线观看| 国精品久久久久久国模美| 尾随美女入室| 啦啦啦在线观看免费高清www| 亚洲不卡免费看| 国产极品天堂在线| 国产黄色视频一区二区在线观看| 在线天堂最新版资源| 精品久久久久久久末码| 97超碰精品成人国产| 精品一区在线观看国产| 亚洲精品一二三| 噜噜噜噜噜久久久久久91| 成人午夜精彩视频在线观看| 麻豆国产97在线/欧美| 欧美一级a爱片免费观看看| 精品久久久久久电影网| 久久精品国产a三级三级三级| 五月伊人婷婷丁香| 蜜桃亚洲精品一区二区三区| 免费高清在线观看视频在线观看| 国产精品99久久99久久久不卡 | 亚洲成人手机| 伦理电影大哥的女人| 日韩强制内射视频| 国产欧美日韩精品一区二区| 精品久久国产蜜桃| 99久久精品一区二区三区| 欧美另类一区| 精品视频人人做人人爽| 亚洲国产欧美人成| 18+在线观看网站| 激情五月婷婷亚洲| 久久国产精品男人的天堂亚洲 | 男人爽女人下面视频在线观看| 高清黄色对白视频在线免费看 | 国产精品久久久久久久久免| 亚洲熟女精品中文字幕| 狂野欧美激情性bbbbbb| 亚洲av欧美aⅴ国产| 国产成人aa在线观看| 这个男人来自地球电影免费观看 | 插逼视频在线观看| 99久久精品国产国产毛片| 精品国产一区二区三区久久久樱花 | 少妇裸体淫交视频免费看高清| 内射极品少妇av片p| 日日摸夜夜添夜夜爱| 免费av不卡在线播放| 久久人人爽av亚洲精品天堂 | 亚洲av男天堂| 汤姆久久久久久久影院中文字幕| 一区二区三区免费毛片| 天堂中文最新版在线下载| 最近的中文字幕免费完整| 天天躁夜夜躁狠狠久久av| 色视频www国产| 欧美bdsm另类| 日日啪夜夜爽| 国产国拍精品亚洲av在线观看| a级一级毛片免费在线观看| 狂野欧美激情性xxxx在线观看| 免费看不卡的av| 一级片'在线观看视频| 精品久久久噜噜| a级毛片免费高清观看在线播放| 联通29元200g的流量卡| 91久久精品国产一区二区三区| 成人午夜精彩视频在线观看| 国产高清国产精品国产三级 | 色婷婷av一区二区三区视频| videos熟女内射| 最近最新中文字幕免费大全7| 精品亚洲乱码少妇综合久久| 国产精品av视频在线免费观看| 丝瓜视频免费看黄片| 亚洲欧洲国产日韩| 黄片无遮挡物在线观看| 成人二区视频| 又大又黄又爽视频免费| 亚洲国产精品国产精品| 国产一区亚洲一区在线观看| 欧美97在线视频| 亚洲欧洲日产国产| 在线播放无遮挡| av国产精品久久久久影院| 免费av不卡在线播放| 2021少妇久久久久久久久久久| 午夜福利网站1000一区二区三区| 在线免费观看不下载黄p国产| 亚洲美女黄色视频免费看| 美女视频免费永久观看网站| 又大又黄又爽视频免费| 新久久久久国产一级毛片| 九色成人免费人妻av| 欧美高清成人免费视频www| 能在线免费看毛片的网站| 国产精品偷伦视频观看了| 国产成人a区在线观看| av在线播放精品| 国产乱人视频| 中文在线观看免费www的网站| 在线观看美女被高潮喷水网站| 亚洲av日韩在线播放| 国产爽快片一区二区三区| 亚洲精品国产色婷婷电影| 老司机影院毛片| 国产乱来视频区| 五月玫瑰六月丁香| 国产精品一区二区三区四区免费观看| 一本一本综合久久| 日韩中字成人| 久热这里只有精品99| 日韩欧美精品免费久久| 欧美日韩视频精品一区| 国产伦精品一区二区三区四那| 亚洲国产最新在线播放| 久久久午夜欧美精品| 精品亚洲成a人片在线观看 | 免费观看av网站的网址| 超碰97精品在线观看| 国产成人精品婷婷| 黄色日韩在线| 国产淫片久久久久久久久| 欧美日韩在线观看h| 直男gayav资源| 国产乱人视频| 久久久精品94久久精品| 国产精品国产av在线观看| 国产亚洲最大av| 欧美激情国产日韩精品一区| 久久韩国三级中文字幕| 亚洲av欧美aⅴ国产| 女的被弄到高潮叫床怎么办| 成人国产麻豆网| 一级爰片在线观看| 色哟哟·www| 99re6热这里在线精品视频| 王馨瑶露胸无遮挡在线观看| av网站免费在线观看视频| 特大巨黑吊av在线直播| 日韩制服骚丝袜av| 亚洲人成网站高清观看| 黑人猛操日本美女一级片| 色5月婷婷丁香| 久久人人爽av亚洲精品天堂 | 有码 亚洲区| www.av在线官网国产| 国产精品嫩草影院av在线观看| 中文欧美无线码| 丝袜喷水一区|