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

    含間隙RU-RPR解耦并聯(lián)機(jī)構(gòu)混沌與沖擊現(xiàn)象

    2017-02-15 00:48:19侯雨雷井國(guó)寧曾達(dá)幸邱雪松李慧劍
    振動(dòng)與沖擊 2017年1期
    關(guān)鍵詞:并聯(lián)間隙沖擊

    侯雨雷, 汪 毅, 井國(guó)寧, 曾達(dá)幸, 邱雪松, 李慧劍

    (1. 燕山大學(xué) 機(jī)械工程學(xué)院,河北 秦皇島 066004;2.燕山大學(xué) 河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗(yàn)室,河北 秦皇島 066004)

    含間隙RU-RPR解耦并聯(lián)機(jī)構(gòu)混沌與沖擊現(xiàn)象

    侯雨雷1, 汪 毅1, 井國(guó)寧1, 曾達(dá)幸1, 邱雪松1, 李慧劍2

    (1. 燕山大學(xué) 機(jī)械工程學(xué)院,河北 秦皇島 066004;2.燕山大學(xué) 河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗(yàn)室,河北 秦皇島 066004)

    針對(duì)含間隙機(jī)構(gòu)中的混沌與沖擊現(xiàn)象,以自主提出的RU-RPR(R為轉(zhuǎn)動(dòng)副,U為虎克鉸,P為移動(dòng)副)兩轉(zhuǎn)動(dòng)解耦并聯(lián)機(jī)構(gòu)為研究對(duì)象,考慮運(yùn)動(dòng)副間隙,結(jié)合Lankarani-Nikravesh接觸力模型與拉格朗日方程建立間隙機(jī)構(gòu)動(dòng)力學(xué)模型,分析不同的驅(qū)動(dòng)速度和摩擦因數(shù)對(duì)機(jī)構(gòu)中混沌現(xiàn)象與沖擊的影響,同時(shí)探討機(jī)構(gòu)的穩(wěn)定性與沖擊現(xiàn)象的關(guān)系,并分析增設(shè)彈簧對(duì)機(jī)構(gòu)動(dòng)力學(xué)特性的影響。研究結(jié)果表明,改變驅(qū)動(dòng)速度與摩擦因數(shù)時(shí),沖擊現(xiàn)象仍然存在,而機(jī)構(gòu)可以由混沌運(yùn)動(dòng)變?yōu)橹芷谶\(yùn)動(dòng),穩(wěn)定性增強(qiáng);機(jī)構(gòu)的穩(wěn)定性與沖擊現(xiàn)象并沒(méi)有必然的聯(lián)系,且增加彈簧可以顯著的減弱機(jī)構(gòu)中的沖擊程度。

    并聯(lián)機(jī)構(gòu);間隙;動(dòng)力學(xué);混沌;沖擊

    隨著機(jī)械產(chǎn)品日益朝著高速、重載、高精度的目標(biāo)發(fā)展,對(duì)機(jī)構(gòu)動(dòng)態(tài)性能的要求也越來(lái)越高,而不可避免的運(yùn)動(dòng)副間隙的影響愈發(fā)引起眾多學(xué)者和工程界的重視。并聯(lián)機(jī)構(gòu)的獨(dú)特優(yōu)勢(shì)及其廣泛應(yīng)用已為業(yè)界認(rèn)可,作為多維的強(qiáng)非線性系統(tǒng),并聯(lián)機(jī)構(gòu)包含著豐富的動(dòng)力學(xué)特性,其能否適應(yīng)更高的作業(yè)條件要求而繼續(xù)得以應(yīng)用并發(fā)展,將取決于動(dòng)力學(xué)這一關(guān)鍵問(wèn)題的解決程度。

    目前,很多學(xué)者對(duì)含間隙機(jī)構(gòu)的動(dòng)力學(xué)建模[1-2]、動(dòng)態(tài)特性[3-6]、混沌現(xiàn)象[7-10]、實(shí)驗(yàn)分析[11]等,已經(jīng)進(jìn)行了深入而廣泛的研究。時(shí)兵等[12]針對(duì)大型重載機(jī)構(gòu)虛擬樣機(jī)中的間隙旋轉(zhuǎn)鉸,基于間隙矢量模型對(duì)一汽輪機(jī)閥門機(jī)構(gòu)進(jìn)行仿真計(jì)算,分析鉸間隙、制造和裝配誤差等因素對(duì)機(jī)構(gòu)動(dòng)力學(xué)特性的影響。BAI等[13]采用一種混合接觸力模型分析含間隙機(jī)構(gòu)的磨損。RAHMANIAN等[8]繪制了含間隙曲柄滑塊機(jī)構(gòu)隨間隙尺寸變化的分岔圖。魏曉輝等[14]分析了裝配間隙對(duì)起落架收放機(jī)構(gòu)的影響,并提出了改進(jìn)設(shè)計(jì)。余躍慶等[15]針對(duì)運(yùn)動(dòng)副間隙引起的并聯(lián)機(jī)器人定位誤差,探討了運(yùn)動(dòng)副間隙誤差引起的系統(tǒng)結(jié)構(gòu)誤差的變化規(guī)律,并提出優(yōu)化補(bǔ)償方法。趙剛練等[16]研究了含空間圓柱副間隙的動(dòng)力學(xué)建模方法,并分析了圓柱副間隙對(duì)其他構(gòu)件的動(dòng)態(tài)特性的影響。王庚祥等[17]對(duì)含空間球面副間隙對(duì)新型減振平臺(tái)4-SPS/CU并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)性能影響進(jìn)行研究,分析球面副間隙對(duì)該并聯(lián)機(jī)構(gòu)運(yùn)動(dòng)學(xué)特性及驅(qū)動(dòng)力的影響。趙寬等[18]分析了含間隙的曲柄滑塊機(jī)構(gòu)在考慮運(yùn)動(dòng)副間隙、摩擦力及系統(tǒng)參數(shù)具有隨機(jī)性時(shí)的動(dòng)力響應(yīng)問(wèn)題。

    混沌是非線性動(dòng)力系統(tǒng)所特有的且普遍存在的現(xiàn)象,混沌現(xiàn)象在有些情況下于機(jī)構(gòu)性能有利,而在有些情況下將產(chǎn)生不良的影響,導(dǎo)致本來(lái)認(rèn)為無(wú)論如何是安全和穩(wěn)定的系統(tǒng)失靈,為工程系統(tǒng)的響應(yīng)和可靠性增加了不確定因素。對(duì)并聯(lián)機(jī)構(gòu)進(jìn)行混沌現(xiàn)象辨識(shí),進(jìn)而體現(xiàn)于動(dòng)態(tài)設(shè)計(jì)中并實(shí)現(xiàn)有效的控制,對(duì)機(jī)構(gòu)動(dòng)態(tài)特性的優(yōu)劣具有關(guān)鍵性的作用。

    在間隙機(jī)構(gòu)中,運(yùn)動(dòng)副元素之間的分離碰撞將產(chǎn)生沖擊現(xiàn)象,嚴(yán)重的沖擊會(huì)導(dǎo)致機(jī)構(gòu)的運(yùn)動(dòng)精度降低,加劇磨損,對(duì)機(jī)構(gòu)產(chǎn)生破壞,降低機(jī)構(gòu)的使用壽命?,F(xiàn)有的文獻(xiàn)中,含間隙機(jī)構(gòu)動(dòng)態(tài)性能的分析已經(jīng)比較普遍,而針對(duì)引起間隙機(jī)構(gòu)動(dòng)態(tài)性能惡化的沖擊現(xiàn)象卻鮮有學(xué)者進(jìn)行深入的研究。

    機(jī)構(gòu)的穩(wěn)定性是指機(jī)構(gòu)保持原有運(yùn)動(dòng)狀態(tài)的能力。當(dāng)給定機(jī)構(gòu)的驅(qū)動(dòng)為周期性時(shí),如果機(jī)構(gòu)處于周期運(yùn)動(dòng)狀態(tài),則機(jī)構(gòu)是穩(wěn)定的,而如果機(jī)構(gòu)處于混沌運(yùn)動(dòng)狀態(tài),則機(jī)構(gòu)是不穩(wěn)定的。運(yùn)動(dòng)副的間隙會(huì)造成機(jī)構(gòu)運(yùn)動(dòng)的不確定性,機(jī)構(gòu)的動(dòng)態(tài)特性在某種條件下會(huì)劇烈地變化,而且呈現(xiàn)出混沌特性[19]。OLYAEI等[19]通過(guò)某種控制方法使機(jī)構(gòu)由混沌狀態(tài)變成不混沌狀態(tài),從而使機(jī)構(gòu)由不穩(wěn)定狀態(tài)變成穩(wěn)定狀態(tài)。ZHANG等[4]對(duì)含多間隙的平面并聯(lián)機(jī)構(gòu)進(jìn)行動(dòng)力學(xué)建模,分析了載荷、運(yùn)動(dòng)速度和運(yùn)動(dòng)軌跡對(duì)機(jī)構(gòu)穩(wěn)定性的影響。王威等[20]以平面四桿機(jī)構(gòu)為例研究了間隙與干摩擦對(duì)機(jī)構(gòu)運(yùn)行穩(wěn)定的影響。

    本文即在此背景下,以含間隙RU-RPR解耦并聯(lián)機(jī)構(gòu)[21]為研究對(duì)象,利用Lankarani-Nikravesh接觸力模型計(jì)算法向接觸力,將修正的Coulomb摩擦力法則嵌入到切向接觸力的計(jì)算過(guò)程中,應(yīng)用拉格朗日方程建立間隙機(jī)構(gòu)動(dòng)力學(xué)模型,采用變步長(zhǎng)四階龍格庫(kù)塔法予以數(shù)值求解,辨識(shí)解耦并聯(lián)機(jī)構(gòu)非線性動(dòng)力學(xué)中的混沌現(xiàn)象,分析不同的因素對(duì)機(jī)構(gòu)中混沌與沖擊現(xiàn)象的影響,探討機(jī)構(gòu)的穩(wěn)定性與混沌現(xiàn)象間關(guān)系,并分析增設(shè)彈簧對(duì)機(jī)構(gòu)沖擊現(xiàn)象的影響,以期對(duì)并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)研究有所深化,為研制高效、高速和穩(wěn)定的并聯(lián)裝備提供理論依據(jù)與技術(shù)支撐。

    1 轉(zhuǎn)動(dòng)副間隙中的接觸力分析

    (1)

    其單位矢量為

    n=e/e

    (2)

    圖1 轉(zhuǎn)動(dòng)副間隙模型Fig.1 Gap model of the rotation pair

    (3)

    式中:Ri、Rj分別為軸承與軸的半徑。

    式(3)兩邊對(duì)時(shí)間求導(dǎo),可得碰撞點(diǎn)速度為

    (4)

    將碰撞點(diǎn)的相對(duì)速度投影到碰撞平面的切向和法向平面,可得相對(duì)的法向速度vn和切向速度vt,其標(biāo)量形式為

    (5)

    式中:切向矢量t可以通過(guò)將n逆時(shí)針旋轉(zhuǎn)90°獲得。

    軸承和軸的穿透深度為

    δ=e-c

    (6)

    式中:c為半徑間隙,且c=Ri-Rj。當(dāng)δ為正值時(shí),軸承與軸接觸,反之,軸承與軸沒(méi)有接觸。

    碰撞過(guò)程中必然產(chǎn)生碰撞力,導(dǎo)致一定的沖擊、振動(dòng),從而對(duì)機(jī)構(gòu)的性能造成不利影響。因此,采用一種合適的接觸力模型來(lái)估算碰撞過(guò)程中的碰撞力非常重要。本文選用Lankarani-Nikravesh接觸力模型[22],該模型既考慮了碰撞過(guò)程中的能量損失,又較全面地考慮了碰撞體的材料屬性、局部彈性變形以及碰撞速度等信息[23],廣泛用于含間隙機(jī)構(gòu)的動(dòng)力學(xué)分析。軸對(duì)軸承的法向接觸力的表達(dá)式為

    (8)

    (9)

    式中:cf為摩擦因數(shù);v0、v1為給定的速度界限。

    該法則可以防止切向接觸力的方向突變,從而有利于數(shù)值積分的穩(wěn)定。

    軸對(duì)軸承的碰撞力可表示為

    (10)

    式中: fix、fiy分別為軸對(duì)軸承的碰撞力在X方向和Y方向上的分量。

    2 含間隙RU-RPR解耦并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)建模

    2.1 RU-RPR解耦并聯(lián)機(jī)構(gòu)組成

    RU-RPR兩轉(zhuǎn)動(dòng)解耦并聯(lián)機(jī)構(gòu)如圖2所示,該機(jī)構(gòu)由定平臺(tái)、動(dòng)平臺(tái)以及連接兩平臺(tái)的兩條支鏈組成,從定平臺(tái)到動(dòng)平臺(tái),第一條分支鏈依次由轉(zhuǎn)動(dòng)副和虎克鉸相連,第二條分支鏈依次由轉(zhuǎn)動(dòng)副、移動(dòng)副和轉(zhuǎn)動(dòng)副相連。其中,第二條分支中移動(dòng)副的移動(dòng)方向垂直于該分支兩個(gè)轉(zhuǎn)動(dòng)副的軸線,兩分支中與固定平臺(tái)相連的兩轉(zhuǎn)動(dòng)副軸線相互平行且與虎克鉸的一條軸線平行,虎克鉸的另一條軸線與第二分支中與動(dòng)平臺(tái)相連的轉(zhuǎn)動(dòng)副軸線重合。

    圖2 RU-RPR兩轉(zhuǎn)動(dòng)并聯(lián)機(jī)構(gòu)組成示意圖Fig.2 Composition schematic of RU-RPR 2-DOF rotational parallel mechanism

    在定平臺(tái)上固連定坐標(biāo)系O-XYZ,其中,坐標(biāo)原點(diǎn)O為與連線中點(diǎn)在定平臺(tái)平面的投影,Z軸豎直向上,Y軸與和定平臺(tái)相連的兩轉(zhuǎn)動(dòng)副A、D軸線垂直并指向A;在動(dòng)平臺(tái)上固連動(dòng)坐標(biāo)系o-xyz,其中,坐標(biāo)原點(diǎn)o建立在動(dòng)平臺(tái)幾何中心,y軸和分支2中與動(dòng)平臺(tái)相連的轉(zhuǎn)動(dòng)副的軸線共線,方向指向虎克鉸,z軸與動(dòng)平臺(tái)平面垂直,方向向上。機(jī)構(gòu)處于初始位置時(shí),AB與BC共線。

    以RU-RPR并聯(lián)機(jī)構(gòu)第一分支的轉(zhuǎn)動(dòng)副A和第二分支與動(dòng)平臺(tái)相連的轉(zhuǎn)動(dòng)副C作為主動(dòng)副,θ與γ分別表示轉(zhuǎn)動(dòng)副A、轉(zhuǎn)動(dòng)副C的旋轉(zhuǎn)角度;α與β分別表示動(dòng)平臺(tái)繞動(dòng)坐標(biāo)系的z軸和y旋轉(zhuǎn)角度。

    θ、γ為輸入角度,而α和β為輸出角度。顯見(jiàn),θ與α一一對(duì)應(yīng),而γ與β也是一一對(duì)應(yīng)關(guān)系,所以,RU-RPR并聯(lián)機(jī)構(gòu)繞y軸和z軸的兩轉(zhuǎn)動(dòng)自由度即兩轉(zhuǎn)動(dòng)互不影響。因此,在研究轉(zhuǎn)動(dòng)副D處的間隙時(shí),可暫不考慮機(jī)構(gòu)繞y軸的轉(zhuǎn)動(dòng),而只考慮動(dòng)平臺(tái)繞z軸(鉛垂方向)的轉(zhuǎn)動(dòng),故可將RU-RPR解耦并聯(lián)機(jī)構(gòu)視作一個(gè)平面四桿機(jī)構(gòu),其等效機(jī)構(gòu)簡(jiǎn)圖如圖3所示。

    圖3 含間隙RU-RPR并聯(lián)機(jī)構(gòu)等效機(jī)構(gòu)簡(jiǎn)圖Fig.3 Schematic diagram of equivalent agencies of RU-RPR parallel mechanism with clearance

    2.2 含間隙RU-RPR并聯(lián)機(jī)構(gòu)運(yùn)動(dòng)分析

    如圖3所示,針對(duì)等效的含間隙RU-RPR解耦并聯(lián)機(jī)構(gòu)建立平面坐標(biāo)系X′O′Y′,并對(duì)后續(xù)動(dòng)力學(xué)分析所需相關(guān)參數(shù)定義如下:l1、θ1、J1分別表示桿AB的長(zhǎng)度、轉(zhuǎn)角以及對(duì)轉(zhuǎn)動(dòng)連接副A的轉(zhuǎn)動(dòng)慣量;l2、ls2、θ2、m2、Js2分別表示桿BC的長(zhǎng)度、質(zhì)心B0距運(yùn)動(dòng)副B的距離、轉(zhuǎn)角(與X′軸正向的夾角)、質(zhì)量以及對(duì)自身質(zhì)心的轉(zhuǎn)動(dòng)慣量;m3、ls3、Js3分別表示桿CE的質(zhì)量、質(zhì)心C0與C之間的距離以及對(duì)自身質(zhì)心的轉(zhuǎn)動(dòng)慣量;m4、ls4、Js4分別表示桿DE的質(zhì)量、質(zhì)心D0距運(yùn)動(dòng)副E的距離以及對(duì)自身質(zhì)心的轉(zhuǎn)動(dòng)慣量;l表示CD之間的距離,l5表示機(jī)架的長(zhǎng)度;x,y分別表示間隙(圖中小圓表示軸,大圓表示軸承)的橫向、縱向分量,即軸和軸承間的橫向、縱向相對(duì)位移分量。

    由于間隙的存在,此時(shí)該四桿機(jī)構(gòu)有三個(gè)自由度。設(shè)定機(jī)構(gòu)的廣義坐標(biāo)為:桿AB的轉(zhuǎn)動(dòng)角度θ1,軸銷與軸承橫向相對(duì)位移x,軸銷與軸承的縱向相對(duì)位移y。通過(guò)引進(jìn)擾動(dòng)坐標(biāo),建立含間隙解耦并聯(lián)機(jī)構(gòu)的運(yùn)動(dòng)學(xué)模型。

    對(duì)無(wú)間隙四桿機(jī)構(gòu),有如下的封閉方程

    (11)

    式中:θ20為無(wú)間隙機(jī)構(gòu)中BC桿轉(zhuǎn)角,l0為無(wú)間隙機(jī)構(gòu)中CD桿的長(zhǎng)度,cθ1=cosθ1,sθ1=sinθ1,其余類似。

    對(duì)D處含一個(gè)間隙的四桿機(jī)構(gòu),有如下封閉方程

    (12)

    利用擾動(dòng)坐標(biāo)法,含間隙機(jī)構(gòu)的轉(zhuǎn)角均由無(wú)間隙機(jī)構(gòu)的轉(zhuǎn)角加上一個(gè)小擾動(dòng)角構(gòu)成,桿長(zhǎng)度由無(wú)間隙機(jī)構(gòu)的長(zhǎng)度加上一個(gè)小的擾動(dòng)長(zhǎng)度構(gòu)成。于是有

    (13)

    將式(13)代入式(12)中,根據(jù)三角公式和等價(jià)無(wú)窮小概念化簡(jiǎn),可得

    (14)

    聯(lián)立式(11)和式(14),得

    (15)

    則式(15)可簡(jiǎn)化為

    (16)

    整個(gè)RU-RPR并聯(lián)機(jī)構(gòu)的動(dòng)能T及勢(shì)能V分別為

    (17)

    式中:F=Js2+m2ls22+Js3+m3l22+Js4+m3ls32+m4ls42,G=m3l1l2+m2l1ls2。

    拉格朗日函數(shù)為

    L=T-V

    (18)

    進(jìn)而,拉格朗日運(yùn)動(dòng)方程可表示為

    (19)

    式中:Qj表示非保守系統(tǒng)中對(duì)應(yīng)廣義坐標(biāo)qj的廣義力,可以表示為

    (20)

    式中:Fi和Mi分別為作用在物體i質(zhì)心的外力、外力矩,ri、φi分別為物體i質(zhì)心的移動(dòng)位移和轉(zhuǎn)動(dòng)角度。

    在小學(xué)數(shù)學(xué)教學(xué)中培養(yǎng)學(xué)生的自主學(xué)習(xí)能力,需要加強(qiáng)對(duì)數(shù)學(xué)課程的策略改革,創(chuàng)新教學(xué)方式,激發(fā)學(xué)生的學(xué)習(xí)興趣,促使學(xué)生積極進(jìn)行學(xué)習(xí)活動(dòng),鼓勵(lì)學(xué)生進(jìn)行自主學(xué)習(xí),提高學(xué)生的自主學(xué)習(xí)能力,促進(jìn)學(xué)生素質(zhì)能力的培養(yǎng)。

    由式(20),廣義力可以表示為

    (21)

    將式(18)、式(21)代入式(19),可得

    (22)

    (23)

    式中:

    3 含間隙RU-RPR并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)數(shù)值仿真

    機(jī)構(gòu)中軸承與軸之間存在碰撞與非碰撞兩種狀態(tài),而在分析每次碰撞的初始碰撞速度及碰撞面的方向時(shí),精確的檢測(cè)出碰撞點(diǎn)至關(guān)重要。已知軸承與軸發(fā)生碰撞的條件為

    δ(tn)δ(tn+1)≤0

    (24)

    則在時(shí)間間隔tn與tn+1之間至少發(fā)生一次碰撞。碰撞時(shí)間非常短暫,故為保證計(jì)算效率與精度,采用變步長(zhǎng)四階龍格庫(kù)塔法進(jìn)行數(shù)值積分,當(dāng)碰撞深度大于積分誤差時(shí),積分步長(zhǎng)設(shè)定為前一個(gè)步長(zhǎng)的一半,重新計(jì)算直至滿足誤差要求,本文借鑒文獻(xiàn)[25]中的方法來(lái)控制步長(zhǎng)。

    式(23)包含兩個(gè)二階非線性方程,采用變步長(zhǎng)四階龍格庫(kù)塔法予以求解。機(jī)構(gòu)桿件形狀參見(jiàn)圖2,多選用長(zhǎng)方形和圓柱形,其材料為45鋼,則可得到RU-RPR機(jī)構(gòu)相關(guān)參數(shù)如下:

    l1=130 mm,l2=308 mm,ls2=154 mm,ls3=80 mm,ls4=54 mm,l5=524 mm,J1=7.351×10-3kg·mm3,Js2=6.160×10-2kg·mm3,Js3=5.772×10-3kg·mm3,Js4=3.008×10-3kg·mm3,m2=8.516 kg,m3=1.161 kg,m4=1.405 kg。

    機(jī)構(gòu)動(dòng)力學(xué)仿真參數(shù)如表1所示,機(jī)構(gòu)初始時(shí)桿OB與桿BC共線,機(jī)構(gòu)驅(qū)動(dòng)角速度為ω=10π rad/s。選取計(jì)算了2個(gè)周期的仿真結(jié)果,機(jī)構(gòu)碰撞深度、接觸力和軸心軌跡分別如圖4(a)、圖4(b)和圖4(c)所示;圖4(d)為機(jī)構(gòu)Poincare映射圖,計(jì)算了500個(gè)映射點(diǎn)。Poincare映射[8]可以方便地進(jìn)行混沌辨識(shí)。Poincare映射上的一個(gè)點(diǎn)及少數(shù)離散點(diǎn)、閉合曲線、成片的具有分形結(jié)構(gòu)的密集點(diǎn)分別表示系統(tǒng)的周期運(yùn)動(dòng)、擬周期運(yùn)動(dòng)和混沌運(yùn)動(dòng)。

    圖4(a)為軸承與軸之間的穿透深度,由圖4(a)可知,軸與軸承在分離和碰撞之間不斷變換,并非處于連續(xù)接觸的狀態(tài)。圖4(c)中軸心軌跡圖也可體現(xiàn)該特點(diǎn)(其中虛線圓表示軸承邊界,該圓以軸承為中心為圓心,以軸承半徑為半徑)。軸與軸承在不斷的分離與碰撞過(guò)程中,必然產(chǎn)生沖擊現(xiàn)象,而圖4(b)中高頻振蕩的接觸力曲線也印證了此現(xiàn)象,且圖4(a)中穿透深度振蕩的時(shí)間點(diǎn)與圖4(b)中接觸力振蕩的時(shí)間點(diǎn)是一致的。圖4(d)中間隙處相對(duì)位移x與相對(duì)位移y的Poincare映射為具有一定的規(guī)律的密集點(diǎn),表現(xiàn)為混沌狀態(tài)。

    表1 RU-RPR并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)仿真參數(shù)

    圖4 含間隙RU-RPR并聯(lián)機(jī)構(gòu)的動(dòng)力學(xué)響應(yīng)Fig.4 Dynamic responses of RU-RPR parallel mechanism with clearance

    4 RU-RPR并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)特性參數(shù)影響分析

    驅(qū)動(dòng)速度以及摩擦因數(shù)是影響機(jī)構(gòu)動(dòng)態(tài)性能的重要因素。如下將根據(jù)接觸力、軸心軌跡以及Poincare映射分析不同的驅(qū)動(dòng)速度及不同的摩擦因數(shù)對(duì)機(jī)構(gòu)中沖擊現(xiàn)象的影響,同時(shí)探討機(jī)構(gòu)穩(wěn)定性與沖擊現(xiàn)象的相互關(guān)系。

    首先依據(jù)表1所列參數(shù),參照前述仿真方法,依次取驅(qū)動(dòng)速度為2π rad/s、8π rad/s、20π rad/s、40π rad/s,所得相應(yīng)的機(jī)構(gòu)接觸力、軸心軌跡及Poincare映射分別如圖5~圖7所示。

    圖5 不同驅(qū)動(dòng)速度下RU-RPR機(jī)構(gòu)的接觸力Fig.5 Contact forces of RU-RPR mechanism at different driving speeds

    由圖5可看出,隨著驅(qū)動(dòng)速度的增加,接觸力曲線的峰值增大,意味著隨著驅(qū)動(dòng)速度的增大,軸與軸承之間的穿透深度更大,軸與軸承之間的沖擊現(xiàn)象更嚴(yán)重,圖6也體現(xiàn)了該特點(diǎn)。

    圖6 不同驅(qū)動(dòng)速度下RU-RPR機(jī)構(gòu)的軸心軌跡圖Fig.6 Axle center trajectory diagrams of RU-RPR mechanism at different driving speeds

    圖7表明,隨著驅(qū)動(dòng)速度的增大,機(jī)構(gòu)中的混沌現(xiàn)象減弱,甚至形成周期運(yùn)動(dòng),機(jī)構(gòu)的穩(wěn)定性增強(qiáng)。而圖6顯示出在各種驅(qū)動(dòng)速度下,軸與軸承之間始終無(wú)法在整個(gè)運(yùn)動(dòng)階段連續(xù)接觸,即沖擊現(xiàn)象依然存在。由此說(shuō)明機(jī)構(gòu)的穩(wěn)定性與沖擊現(xiàn)象并沒(méi)有必然的聯(lián)系。

    其次,令驅(qū)動(dòng)速度為5πrad/s,按摩擦系數(shù)依次取0.01、0.02、0.05、0.1四種情況,分析機(jī)構(gòu)接觸力、軸心軌跡及Poincare映射分別如圖8~圖10所示。

    圖7 不同驅(qū)動(dòng)速度下RU-RPR機(jī)構(gòu)的Poincare映射Fig.7 Poincare mapping diagrams of RU-RPR mechanism at different driving speeds

    圖8 不同摩擦系數(shù)下RU-RPR機(jī)構(gòu)的接觸力Fig.8 Contact forces of RU-RPR mechanism at different friction coefficients

    由圖8可知,當(dāng)摩擦系數(shù)比較小時(shí),接觸力的峰值隨著摩擦系數(shù)的增大而增大,沖擊現(xiàn)象更為嚴(yán)重;而當(dāng)摩擦系數(shù)較大時(shí),隨著摩擦系數(shù)增大,接觸力的峰值幾乎不變,沖擊現(xiàn)象的程度基本保持穩(wěn)定。

    圖10表明,摩擦系數(shù)越大,機(jī)構(gòu)中的混沌現(xiàn)象越弱,機(jī)構(gòu)的穩(wěn)定性增強(qiáng);但是軸與軸承之間的沖擊碰撞依然存在(參見(jiàn)圖9),也再次說(shuō)明機(jī)構(gòu)的穩(wěn)定性與沖擊現(xiàn)象間并不存在必然的聯(lián)系。

    圖9 不同摩擦系數(shù)下RU-RPR機(jī)構(gòu)的軸心軌跡圖Fig.9 Axle center trajectory diagrams of RU-RPR mechanism at different friction coefficients

    圖10 不同摩擦系數(shù)下RU-RPR機(jī)構(gòu)的Poincare映射Fig.10 Poincare mapping diagrams of RU-RPR mechanism at different friction coefficients

    5 RU-RPR并聯(lián)機(jī)構(gòu)間隙影響的控制

    根據(jù)RU-RPR并聯(lián)機(jī)構(gòu)動(dòng)態(tài)特性參數(shù)影響分析,隨著驅(qū)動(dòng)速度或摩擦系數(shù)的增大,機(jī)構(gòu)中的混沌現(xiàn)象減弱,機(jī)構(gòu)的穩(wěn)定性增強(qiáng),但軸與軸承之間始終無(wú)法連續(xù)接觸,沖擊現(xiàn)象依然存在,甚至更加嚴(yán)重,而這必將降低機(jī)構(gòu)的精度、壽命等,影響其正常工作。

    故此,為減少軸與軸承之間的沖擊,考慮在桿DE的質(zhì)心與機(jī)架OD的幾何中心之間增加一個(gè)彈簧,如圖11所示。設(shè)彈簧原長(zhǎng)為L(zhǎng)s0,剛度系數(shù)為ks,伸長(zhǎng)后彈簧長(zhǎng)度為L(zhǎng)s,則增加彈簧后RU-RPR機(jī)構(gòu)的勢(shì)能為

    (25)

    圖11 增加彈簧后的RU-RPR機(jī)構(gòu)簡(jiǎn)圖Fig.11 RU-RPR mechanism diagram after adding the spring

    為了達(dá)到預(yù)期的效果,使彈簧始終處于拉伸狀態(tài),且設(shè)定彈簧的原長(zhǎng)為0.15 m,剛度系數(shù)為1 000 N/m,驅(qū)動(dòng)速度ω=10π rad/s,其它動(dòng)力學(xué)參數(shù)參照表1,仿真結(jié)果如圖12所示。由圖12(a)可知,增加彈簧后,軸與軸承之間始終連續(xù)接觸,沒(méi)有發(fā)生沖擊現(xiàn)象,圖12(b)、(c)也說(shuō)明該特點(diǎn)。此外,對(duì)比圖12(b)與圖4(b)可知,增加彈簧后,接觸力波動(dòng)趨于平緩,且最大值由6 778降為1 629,降低了76%。圖12(d)則表明增加彈簧后,機(jī)構(gòu)中的混沌現(xiàn)象消失。

    圖12 增設(shè)彈簧后RU-RPR并聯(lián)機(jī)構(gòu)的動(dòng)力學(xué)響應(yīng)Fig.12 Dynamic responses of RU-RPR parallel mechanism after adding spring

    6 結(jié) 論

    (1)利用Lankarani-Nikravesh接觸力模型、修正的Coulomb摩擦力法則和拉格朗日方程建立了含間隙RU-RPR解耦并聯(lián)機(jī)構(gòu)的動(dòng)力學(xué)模型,借助變步長(zhǎng)四階龍格庫(kù)塔法予以求解,動(dòng)力學(xué)仿真結(jié)果表明RU-RPR機(jī)構(gòu)轉(zhuǎn)動(dòng)副中間隙的存在使得機(jī)構(gòu)在某種條件下存在混沌和沖擊現(xiàn)象。

    (2)當(dāng)含間隙RU-RPR機(jī)構(gòu)中驅(qū)動(dòng)速度、摩擦系數(shù)增大時(shí),機(jī)構(gòu)由混沌運(yùn)動(dòng)變?yōu)橹芷谶\(yùn)動(dòng),即機(jī)構(gòu)的穩(wěn)定性提高,而沖擊現(xiàn)象依然存在,說(shuō)明機(jī)構(gòu)的穩(wěn)定性與沖擊現(xiàn)象并沒(méi)有必然的聯(lián)系,且通過(guò)增加彈簧并預(yù)緊能夠有效消除機(jī)構(gòu)中的沖擊現(xiàn)象,對(duì)高可靠性并聯(lián)裝備的設(shè)計(jì)與控制具有一定的理論參考意義和實(shí)踐應(yīng)用價(jià)值。

    [ 1 ] 白爭(zhēng)鋒. 考慮鉸間間隙的機(jī)構(gòu)動(dòng)力學(xué)特性研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué),2011.

    [ 2 ] FLORES P. Dynamic analysis of mechanical systems with imperfect kinematic joints[D].Portugal: Universidade Do Minho, 2004.

    [ 3 ] 王庚祥,劉宏昭. 考慮球面副間隙的4-SPS/CU并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)分析[J]. 機(jī)械工程學(xué)報(bào),2015, 51(1): 43-51. WANG Gengxiang, LIU Hongzhao. Dynamics analysis of 4-SPS/CU parallel mechanism with spherical joint clearance[J]. Journal of Mechanical Engineering, 2015, 51(1): 43-51.

    [ 4 ] ZHANG Xuchong, ZHANG Xianmin, CHEN Zhong. Dynamic analysis of a 3-RRR parallel mechanism with multiple clearance joints[J]. Mechanism and Machine Theory, 2014, 78(8): 105-115.

    [ 5 ] 白爭(zhēng)鋒, 趙陽(yáng), 趙志剛. 考慮運(yùn)動(dòng)副間隙的機(jī)構(gòu)動(dòng)態(tài)特性研究[J]. 振動(dòng)與沖擊, 2011, 30(11): 17-20. BAI Zhengfeng, ZHAO Yang, ZHAO Zhigang. Dynamic characteristics of mechanisms with joint clearance[J]. Journal of Vibration and Shock, 2011, 30(11): 17-20.

    [ 6 ] 張?jiān)隼? 巫世晶, 趙文強(qiáng), 等. 含間隙高速多連桿傳動(dòng)機(jī)構(gòu)動(dòng)力學(xué)特性研究[J]. 振動(dòng)與沖擊, 2014, 33(14): 66-71. ZHANG Zenglei, WU Shijing, ZHAO Wenqiang, et al. Dynamic characteristics of high-speed multi-link transmission mechanisms with clearance[J]. Journal of Vibration and Shock, 2014, 33(14): 66-71.

    [ 7 ] 侯雨雷,汪毅,李明洋,等.一種含間隙兩轉(zhuǎn)動(dòng)一移動(dòng)解耦并聯(lián)機(jī)構(gòu)混沌現(xiàn)象辨識(shí)[J]. 中國(guó)機(jī)械工程,2015, 26(13): 1759-1766. HOU Yulei, WANG Yi, LI Mingyang, et al. Chaos identification of two rotational and one translational decoupled parallel mechanism concerning clearance[J]. China Mechanical Engineering, 2015, 26(13): 1759-1766.

    [ 8 ] RAHMANIAN S, GHAZAVI M R. Bifurcation in planar slider-crank mechanism with revolute clearance joint[J]. Mechanism and Machine Theory, 2015, 91(5): 86-101.

    [ 9 ] TANG Yuanguang, CHANG Zongyu, DONG Xiaogang, et al. Nonlinear dynamics and analysis of a four-bar linkage with clearance[J]. Frontiers of Mechanical Engineering, 2013, 8(2): 160-168.

    [10] 何文運(yùn), 張俊紅, 馬梁. 剪式齒輪傳動(dòng)機(jī)構(gòu)非線性動(dòng)力學(xué)特性研究[J]. 振動(dòng)與沖擊, 2015, 34(15): 115-121. HE Wenyun, ZHANG Junhong, MA Liang. Nonlinear dynamic characteristics of a scissor gear transmission system[J]. Journal of Vibration and Shock, 2015, 34(15): 115-121.

    [11] ERKAYA S, DOGAN S, ULUS S. Effects of joint clearance on the dynamics of a partly compliant mechanism:Numerical and experimental studies[J]. Mechanism and Machine Theory, 2015, 88(2): 125-140.

    [12] 時(shí)兵,金燁. 面向虛擬樣機(jī)的機(jī)構(gòu)間隙旋轉(zhuǎn)鉸建模與動(dòng)力學(xué)仿真[J]. 機(jī)械工程學(xué)報(bào),2009, 45(4): 299-303. SHI Bing, JIN Ye. Dynamic simulation and modeling of revolute clearance joint for virtual prototyping[J]. Journal of Mechanical Engineering, 2009, 45(4): 299-303.

    [13] BAI Zhengfeng, ZHAO Yang, CHEN Jun. Dynamics analysis of planar mechanical system considering revolute clearance joint wear[J]. Tribology International, 2013, 64(3): 85-95.

    [14] 魏小輝, 王鈺龍, 印寅, 等. 某起落架收放聯(lián)動(dòng)機(jī)構(gòu)故障分析及改進(jìn)設(shè)計(jì)[J]. 中國(guó)機(jī)械工程,2014, 25(3): 399-403. WEI Xiaohui, WANG Yulong, YIN Yin, et al. Failure analysis and improved design of a landing gear retraction mechanism[J]. China Mechanical Engineering, 2014, 25(3): 399-403.

    [15] 余躍慶, 田浩. 運(yùn)動(dòng)副間隙引起的并聯(lián)機(jī)器人誤差及其補(bǔ)償[J]. 光學(xué)精密工程,2015, 23(5): 1331-1339. YU Yueqing, TIAN Hao. Error and compensation of parallel robot with joint clearances[J]. Optics and Precision Engineering, 2015, 23(5): 1331-1339.

    [16] 趙剛練, 姜毅, 郝繼光, 等. 考慮圓柱鉸鏈間隙的多剛體系統(tǒng)動(dòng)力學(xué)計(jì)算方法[J]. 振動(dòng)與沖擊, 2013, 32(17): 171-176. ZHAO Ganglian, JING Yi, HAO Jiguang. Computational method of rigid mllltibody system dynamics considering cylindrical joint clearance[J]. Journal of Vibration and Shock, 2013, 32(17): 171-176.

    [17] 王庚祥, 劉宏昭, 鄧培生. 考慮球面副間隙的并聯(lián)機(jī)構(gòu)動(dòng)力學(xué)模型[J]. 振動(dòng)與沖擊, 2014, 33(10): 43-49. WANG Gengxiang, LIU Hongzhao, DENG Peisheng. Dynamic modeling for a parallel mechanism considering spherical joint clearance[J]. Journal of Vibration and Shock, 2014, 33(10): 43-49.

    [18] 趙寬, 陳建軍, 曹鴻鈞, 等. 含間隙曲柄滑塊系統(tǒng)動(dòng)力響應(yīng)的隨機(jī)性分析[J]. 振動(dòng)與沖擊, 2014, 33(11): 130-135. ZHAO Kuan, CHEN Jianjun, CAO Hongjun, et al.Randomness analysis of dynamic responses of a slider-crank system with clearances[J]. Journal of Vibration and Shock, 2014, 33(11): 130-135.

    [19] OLYAEI A A, GHAZAVI M R. Stabilizing slider-crank mechanism with clearance joints[J]. Mechanism and Machine Theory, 2012, 53(3): 17-29.

    [20] 王威, 沈政, 宋玉玲, 等. 含間隙和干摩擦的連桿機(jī)構(gòu)系統(tǒng)動(dòng)力學(xué)研究[J]. 振動(dòng)與沖擊, 2015, 34(18): 210-214. WANG Wei, SHEN Zheng, SONG Yuling. System dynamics of linkage mechanism with clearance and dry friction[J]. Journal of Vibration and Shock, 2015, 34(18): 210-214.

    [21] HOU Yulei, ZENG Daxing, ZHANG Zhanye, et al. A novel two degrees of freedom rotational decoupled parallel mechanism[J]. Applied Mechanics and Materials, 2012, 215/216: 293-296.

    [22] LANKARANI H M, NIKRAVESH P E. Continuous contact force models for impact analysis in multibody systems[J]. Nonlinear Dynamics, 1994, 5(2): 193-207.

    [23] FLORES P, AMBRSIO J, CLARO J C P, et al. Influence of the contact—Impact force model on the dynamic response of multi-body systems[J]. Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, 2006, 220(1): 21-34.

    [25] FLORES P, AMBRSIO J. On the contact detection for contact-impact analysis in multibody systems[J]. Multibody System Dynamics, 2010, 24(5): 103-122.

    Chaos andimpact phenomena of a RU-RPR decoupled parallel mechanism containing clearance

    HOU Yulei1, WANG Yi1, JING Guoning1, ZENG Daxing1, QIU Xuesong1, LI Huijian2

    (1. School of Mechanical Engineering, Yanshan University, Qinhuangdao 066004, China;2. Key Laboratory of Mechanical Reliability for Heavy Equipments and Large Structures of Hebei Province, Yanshan University, Qinhuangdao 066004, China)

    For chaos and impact phenomena in a mechanism with clearance, taking a two-rotation decoupled parallel mechanism RU-RPR (R stands for a revolute pair, U stands for a Hooke joint, P stands for a translation pair) proposed independently as a study object, considering clearance existing in a kinematic pair, combining with Lankarani-Nikravesh contact force model and Lagrange equation, the dynamic equation of the mechanism was established, and its chaos and impact phenomena were studied. Then the influences of different driving speeds and friction coefficients on chaos and impact phenomena were analyzed, and the relationship between the stability of the mechanism and impact phenomenon was investigated at the same time. The influences of an additional spring on dynamic characteristics of the mechanism were analyzed. The results showed that when changing drive speed or friction coefficient, the impact phenomenon still exists, while the motion of the mechanism can be changed from chaotic motion to periodic one, the stability of the mechanism is enhanced; there is no necessary relation between the stability of the mechanism and its impact phenomenon, and the impact level can be significantly weakened by adding springs.

    parallel mechanism; clearance; dynamics; chaos; impact

    國(guó)家自然科學(xué)基金(51205339;51305384);中國(guó)博士后科學(xué)基金(2013M541199)

    2015-08-16 修改稿收到日期:2015-12-30

    侯雨雷 男,博士后,教授,博士生導(dǎo)師,1980年生

    李慧劍 男,博士,教授,博士生導(dǎo)師,1980年生

    TH113

    A

    10.13465/j.cnki.jvs.2017.01.032

    猜你喜歡
    并聯(lián)間隙沖擊
    間隙
    識(shí)別串、并聯(lián)電路的方法
    飛行過(guò)載及安裝間隙對(duì)主安裝節(jié)推力測(cè)量的影響
    緊流形上的Schr?dinger算子的譜間隙估計(jì)
    審批由“串聯(lián)”改“并聯(lián)”好在哪里?
    淺談保護(hù)間隙的利弊與應(yīng)用
    廣西電力(2016年4期)2016-07-10 10:23:38
    并聯(lián)型APF中SVPWM的零矢量分配
    一種軟開關(guān)的交錯(cuò)并聯(lián)Buck/Boost雙向DC/DC變換器
    奧迪Q5換擋沖擊
    奧迪A8L換擋沖擊
    亚洲精品一区蜜桃| 激情 狠狠 欧美| 熟女电影av网| 亚洲人成网站在线观看播放| 小蜜桃在线观看免费完整版高清| 日日摸夜夜添夜夜添av毛片| 免费人成在线观看视频色| 亚洲av成人精品一区久久| 久热久热在线精品观看| 亚洲精品第二区| 黄色视频在线播放观看不卡| 天堂俺去俺来也www色官网| 下体分泌物呈黄色| 免费在线观看成人毛片| 天堂俺去俺来也www色官网| 久久久久久久国产电影| 欧美成人午夜免费资源| 一级爰片在线观看| 国产大屁股一区二区在线视频| 日日啪夜夜爽| 欧美一级a爱片免费观看看| 极品少妇高潮喷水抽搐| 新久久久久国产一级毛片| 一级毛片aaaaaa免费看小| 国产欧美日韩精品一区二区| 人人妻人人爽人人添夜夜欢视频 | 王馨瑶露胸无遮挡在线观看| 一区二区三区四区激情视频| 久久久久久伊人网av| 美女主播在线视频| 成人国产av品久久久| 国产欧美亚洲国产| 网址你懂的国产日韩在线| 国产毛片在线视频| 国产成人精品福利久久| 国内少妇人妻偷人精品xxx网站| 黄片无遮挡物在线观看| 一级毛片 在线播放| 啦啦啦在线观看免费高清www| 亚洲国产av新网站| 日韩中文字幕视频在线看片 | 欧美3d第一页| av线在线观看网站| 尤物成人国产欧美一区二区三区| 精品视频人人做人人爽| 午夜视频国产福利| 人妻一区二区av| 波野结衣二区三区在线| 熟女电影av网| 亚洲av国产av综合av卡| 观看av在线不卡| 黑人高潮一二区| 国产精品女同一区二区软件| av国产精品久久久久影院| 欧美成人一区二区免费高清观看| 丝袜脚勾引网站| 韩国高清视频一区二区三区| 国产日韩欧美亚洲二区| 亚洲欧美成人精品一区二区| h日本视频在线播放| 日本黄大片高清| 国产一区二区在线观看日韩| 亚洲精品乱久久久久久| 久久综合国产亚洲精品| av国产久精品久网站免费入址| 亚洲国产色片| 亚洲人成网站在线播| 久久精品熟女亚洲av麻豆精品| 十八禁网站网址无遮挡 | 精品国产一区二区三区久久久樱花 | tube8黄色片| 国产精品一区www在线观看| 国产黄片视频在线免费观看| a级毛片免费高清观看在线播放| 国产精品精品国产色婷婷| 欧美另类一区| 日本猛色少妇xxxxx猛交久久| 色吧在线观看| 国产成人免费观看mmmm| 国产综合精华液| 中文字幕精品免费在线观看视频 | 高清在线视频一区二区三区| 大又大粗又爽又黄少妇毛片口| 人人妻人人爽人人添夜夜欢视频 | 久久久久国产网址| 免费观看在线日韩| 国产亚洲5aaaaa淫片| 十分钟在线观看高清视频www | 最新中文字幕久久久久| 国产成人免费观看mmmm| 亚洲国产精品一区三区| 午夜老司机福利剧场| 亚洲欧美一区二区三区国产| 欧美高清性xxxxhd video| 欧美变态另类bdsm刘玥| 人妻制服诱惑在线中文字幕| 亚洲国产高清在线一区二区三| 日产精品乱码卡一卡2卡三| 99热网站在线观看| 激情五月婷婷亚洲| 啦啦啦中文免费视频观看日本| 夜夜骑夜夜射夜夜干| 少妇人妻精品综合一区二区| 国产精品久久久久久久久免| tube8黄色片| 日韩,欧美,国产一区二区三区| 高清午夜精品一区二区三区| 91aial.com中文字幕在线观看| 久久久久久久久大av| 国产综合精华液| 水蜜桃什么品种好| 中文字幕免费在线视频6| 国产在线男女| 精品国产一区二区三区久久久樱花 | 中文天堂在线官网| 亚洲成人中文字幕在线播放| 国产精品人妻久久久影院| 久久99热这里只频精品6学生| 一区二区三区乱码不卡18| 男男h啪啪无遮挡| 亚洲国产精品成人久久小说| 久热久热在线精品观看| 午夜免费观看性视频| 九九在线视频观看精品| 午夜日本视频在线| 成人国产av品久久久| 国内揄拍国产精品人妻在线| 又粗又硬又长又爽又黄的视频| 综合色丁香网| av不卡在线播放| 卡戴珊不雅视频在线播放| 春色校园在线视频观看| 成年免费大片在线观看| 熟妇人妻不卡中文字幕| 麻豆成人午夜福利视频| 啦啦啦啦在线视频资源| 亚洲国产精品成人久久小说| 久久精品国产亚洲av天美| 国产精品人妻久久久久久| 国产有黄有色有爽视频| 欧美亚洲 丝袜 人妻 在线| 国产精品国产三级国产av玫瑰| 18禁在线播放成人免费| 毛片女人毛片| 久久99热6这里只有精品| 久久久久久久亚洲中文字幕| 九草在线视频观看| 国产乱人视频| 精品一区二区免费观看| 久久久精品94久久精品| 亚洲va在线va天堂va国产| 熟妇人妻不卡中文字幕| 亚洲精品一二三| 国产成人精品久久久久久| 亚洲欧美日韩无卡精品| 18禁裸乳无遮挡动漫免费视频| 久久久精品94久久精品| 日韩免费高清中文字幕av| 嘟嘟电影网在线观看| 久久久欧美国产精品| 大话2 男鬼变身卡| 精品久久久久久电影网| 日韩成人av中文字幕在线观看| 热99国产精品久久久久久7| 国产v大片淫在线免费观看| 国产亚洲av片在线观看秒播厂| 网址你懂的国产日韩在线| 成人毛片a级毛片在线播放| 精品人妻视频免费看| 国产精品三级大全| 高清欧美精品videossex| 欧美高清性xxxxhd video| 熟女av电影| 永久网站在线| 99久久精品热视频| 免费观看无遮挡的男女| a级毛色黄片| 一级爰片在线观看| 亚洲性久久影院| 亚洲婷婷狠狠爱综合网| 午夜激情福利司机影院| 亚洲欧美成人精品一区二区| 亚洲av国产av综合av卡| 久久ye,这里只有精品| 少妇 在线观看| 小蜜桃在线观看免费完整版高清| 网址你懂的国产日韩在线| 国产成人aa在线观看| 在线观看人妻少妇| 国产av码专区亚洲av| 三级国产精品欧美在线观看| 久久99热这里只有精品18| 免费大片18禁| 51国产日韩欧美| 最新中文字幕久久久久| 精品久久久久久久久亚洲| 超碰av人人做人人爽久久| 水蜜桃什么品种好| 亚洲丝袜综合中文字幕| 亚洲av不卡在线观看| 久久精品国产a三级三级三级| 欧美人与善性xxx| 国产精品偷伦视频观看了| 欧美国产精品一级二级三级 | 久久久久久久精品精品| 婷婷色综合www| 亚洲精品乱码久久久v下载方式| 国产探花极品一区二区| 免费人妻精品一区二区三区视频| 成人18禁高潮啪啪吃奶动态图 | 久久人妻熟女aⅴ| 日日撸夜夜添| 91精品伊人久久大香线蕉| 日本黄色片子视频| 久久久精品免费免费高清| 夜夜看夜夜爽夜夜摸| 高清日韩中文字幕在线| 国产成人a∨麻豆精品| 免费看av在线观看网站| 国产午夜精品一二区理论片| 国产高潮美女av| 午夜激情久久久久久久| 国语对白做爰xxxⅹ性视频网站| 国产一级毛片在线| 蜜桃在线观看..| 精品人妻偷拍中文字幕| 九色成人免费人妻av| 国产精品久久久久成人av| 一级爰片在线观看| 夜夜看夜夜爽夜夜摸| 国产男人的电影天堂91| 婷婷色综合www| 久久久a久久爽久久v久久| 熟女电影av网| 大码成人一级视频| 中文字幕精品免费在线观看视频 | 黄色一级大片看看| 少妇人妻一区二区三区视频| 日韩成人av中文字幕在线观看| 国产69精品久久久久777片| 91狼人影院| 亚洲成人中文字幕在线播放| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品古装| 麻豆成人av视频| 亚洲美女黄色视频免费看| 深爱激情五月婷婷| 香蕉精品网在线| 日韩,欧美,国产一区二区三区| 日韩不卡一区二区三区视频在线| 国产成人a∨麻豆精品| 丝瓜视频免费看黄片| h日本视频在线播放| 久久久成人免费电影| 99久久精品国产国产毛片| 又黄又爽又刺激的免费视频.| 老司机影院毛片| 亚洲美女视频黄频| 天美传媒精品一区二区| 国产美女午夜福利| 国产精品人妻久久久久久| 亚洲精品aⅴ在线观看| 免费av中文字幕在线| 国产成人aa在线观看| kizo精华| 亚洲欧美精品专区久久| 校园人妻丝袜中文字幕| 色视频www国产| 免费看光身美女| 3wmmmm亚洲av在线观看| 美女福利国产在线 | 国产熟女欧美一区二区| 亚洲精品视频女| 美女主播在线视频| 一本—道久久a久久精品蜜桃钙片| 欧美国产精品一级二级三级 | 日本vs欧美在线观看视频 | 午夜视频国产福利| 久久亚洲国产成人精品v| 少妇高潮的动态图| 国产男女超爽视频在线观看| 肉色欧美久久久久久久蜜桃| 亚洲成人av在线免费| 亚洲精品日韩在线中文字幕| 国产真实伦视频高清在线观看| 亚州av有码| 成人影院久久| 美女内射精品一级片tv| 老师上课跳d突然被开到最大视频| 美女脱内裤让男人舔精品视频| 日日撸夜夜添| 高清欧美精品videossex| 成人综合一区亚洲| 久久久亚洲精品成人影院| 我要看日韩黄色一级片| 国产精品久久久久久精品古装| 日日摸夜夜添夜夜添av毛片| 只有这里有精品99| 一级二级三级毛片免费看| 人妻系列 视频| 日韩国内少妇激情av| 国产精品秋霞免费鲁丝片| 2021少妇久久久久久久久久久| 联通29元200g的流量卡| 夜夜看夜夜爽夜夜摸| 国产精品三级大全| 韩国高清视频一区二区三区| 岛国毛片在线播放| 亚洲欧美成人综合另类久久久| 天天躁日日操中文字幕| 精品亚洲成国产av| 日韩av不卡免费在线播放| 欧美日韩视频高清一区二区三区二| 成年av动漫网址| 大话2 男鬼变身卡| 欧美性感艳星| 蜜桃亚洲精品一区二区三区| 黄色配什么色好看| 中国三级夫妇交换| 熟女人妻精品中文字幕| 亚洲精华国产精华液的使用体验| 日韩人妻高清精品专区| 欧美zozozo另类| 一级毛片aaaaaa免费看小| 只有这里有精品99| 久久 成人 亚洲| 联通29元200g的流量卡| av播播在线观看一区| 少妇猛男粗大的猛烈进出视频| 久热这里只有精品99| 亚洲欧洲国产日韩| 欧美xxxx性猛交bbbb| 人妻制服诱惑在线中文字幕| 自拍偷自拍亚洲精品老妇| 91精品伊人久久大香线蕉| 久久久精品免费免费高清| 国产真实伦视频高清在线观看| 99久久精品国产国产毛片| 另类亚洲欧美激情| 久久久久精品久久久久真实原创| 亚洲精品国产av成人精品| 成人午夜精彩视频在线观看| 一区在线观看完整版| 亚洲va在线va天堂va国产| 精品一品国产午夜福利视频| 中文欧美无线码| 好男人视频免费观看在线| 丝袜喷水一区| 欧美亚洲 丝袜 人妻 在线| 又大又黄又爽视频免费| 亚洲av福利一区| 国产成人a∨麻豆精品| 性色av一级| 在线观看免费视频网站a站| 精品酒店卫生间| 街头女战士在线观看网站| 男女国产视频网站| 国内精品宾馆在线| 亚洲美女搞黄在线观看| 国产av一区二区精品久久 | 国产免费一级a男人的天堂| 国产成人91sexporn| 亚洲国产精品一区三区| 日本一二三区视频观看| 欧美 日韩 精品 国产| 久热这里只有精品99| av天堂中文字幕网| 久久99精品国语久久久| 亚州av有码| 91久久精品国产一区二区成人| 亚洲人与动物交配视频| 亚洲精品一区蜜桃| 国产探花极品一区二区| 久久久亚洲精品成人影院| 在线观看美女被高潮喷水网站| 丝袜喷水一区| 国产 一区精品| 我要看黄色一级片免费的| 久久这里有精品视频免费| 哪个播放器可以免费观看大片| 美女cb高潮喷水在线观看| 人妻一区二区av| 99热网站在线观看| 最近2019中文字幕mv第一页| .国产精品久久| 麻豆国产97在线/欧美| 91精品国产九色| 久久精品人妻少妇| 观看美女的网站| 丝袜喷水一区| 欧美日韩在线观看h| 女性被躁到高潮视频| 99re6热这里在线精品视频| 亚洲欧洲日产国产| 高清毛片免费看| 男女边摸边吃奶| 春色校园在线视频观看| 伦精品一区二区三区| xxx大片免费视频| 老熟女久久久| 一区二区三区精品91| 2021少妇久久久久久久久久久| 男女国产视频网站| 亚洲精品第二区| 欧美日韩一区二区视频在线观看视频在线| 成人漫画全彩无遮挡| 国产高清三级在线| 日本色播在线视频| 日韩一本色道免费dvd| 亚洲国产精品一区三区| av在线app专区| 欧美精品一区二区免费开放| videossex国产| 麻豆成人av视频| 国产精品三级大全| 蜜桃久久精品国产亚洲av| 午夜福利在线在线| 国产精品久久久久成人av| 国产在线免费精品| 国产精品欧美亚洲77777| 高清在线视频一区二区三区| 男女边摸边吃奶| 高清欧美精品videossex| 人妻一区二区av| 久久久国产一区二区| 久久久久久久国产电影| 国产乱人视频| 美女内射精品一级片tv| 久久亚洲国产成人精品v| 我的女老师完整版在线观看| 男人狂女人下面高潮的视频| 高清黄色对白视频在线免费看 | 国产综合精华液| 青春草视频在线免费观看| 婷婷色综合大香蕉| 免费黄色在线免费观看| 美女脱内裤让男人舔精品视频| 成人影院久久| 亚洲av成人精品一区久久| 久久精品久久久久久噜噜老黄| 国产成人aa在线观看| 亚洲国产高清在线一区二区三| 校园人妻丝袜中文字幕| 国产av一区二区精品久久 | 亚洲国产最新在线播放| 九草在线视频观看| 18禁裸乳无遮挡动漫免费视频| 日本一二三区视频观看| 高清毛片免费看| 色综合色国产| 汤姆久久久久久久影院中文字幕| 日韩三级伦理在线观看| 一区二区三区免费毛片| 97在线人人人人妻| 欧美区成人在线视频| 欧美日韩精品成人综合77777| 免费播放大片免费观看视频在线观看| 国产精品av视频在线免费观看| 日韩亚洲欧美综合| 免费观看性生交大片5| 久久久久久久亚洲中文字幕| 草草在线视频免费看| 亚洲人与动物交配视频| 国产高潮美女av| 欧美成人一区二区免费高清观看| 亚洲精华国产精华液的使用体验| 亚洲一级一片aⅴ在线观看| 婷婷色av中文字幕| 亚洲最大成人中文| 激情 狠狠 欧美| 精品一区二区三卡| 亚洲国产av新网站| 男人狂女人下面高潮的视频| 一区二区三区精品91| 三级国产精品片| 亚洲精品视频女| 国产淫片久久久久久久久| 日本-黄色视频高清免费观看| av卡一久久| 岛国毛片在线播放| 欧美一区二区亚洲| 美女中出高潮动态图| 七月丁香在线播放| 国产在线一区二区三区精| 久久久久精品久久久久真实原创| 成人二区视频| 中文欧美无线码| 性色av一级| 中文欧美无线码| 美女视频免费永久观看网站| av播播在线观看一区| 啦啦啦啦在线视频资源| 国产精品国产三级专区第一集| 水蜜桃什么品种好| 日韩一本色道免费dvd| 丝袜脚勾引网站| 亚洲欧美精品自产自拍| 一级毛片我不卡| 亚洲成色77777| 91精品伊人久久大香线蕉| 日韩,欧美,国产一区二区三区| 国产亚洲最大av| 在线观看人妻少妇| 亚洲综合色惰| 中文在线观看免费www的网站| 亚洲av.av天堂| 日日啪夜夜爽| 99热网站在线观看| 97在线人人人人妻| 七月丁香在线播放| 涩涩av久久男人的天堂| 久久久久久久久大av| 内地一区二区视频在线| 26uuu在线亚洲综合色| 啦啦啦啦在线视频资源| 中文欧美无线码| 天堂俺去俺来也www色官网| 日日啪夜夜撸| 青春草国产在线视频| 日日啪夜夜撸| 久久人人爽人人片av| 国产一区二区三区av在线| 久久精品国产亚洲网站| 午夜福利在线在线| 色视频www国产| 日本av免费视频播放| 国产视频内射| 免费高清在线观看视频在线观看| 99热6这里只有精品| 日产精品乱码卡一卡2卡三| 狂野欧美激情性bbbbbb| 天天躁夜夜躁狠狠久久av| 51国产日韩欧美| 国产黄色免费在线视频| 麻豆乱淫一区二区| 国产精品成人在线| 美女高潮的动态| 91aial.com中文字幕在线观看| 乱系列少妇在线播放| 国产高潮美女av| 亚洲人与动物交配视频| 亚洲av中文av极速乱| 久久精品久久久久久久性| 欧美日韩在线观看h| 中文字幕制服av| 日韩不卡一区二区三区视频在线| 内地一区二区视频在线| 日本vs欧美在线观看视频 | 亚洲图色成人| 亚洲av男天堂| 亚洲欧美精品专区久久| 国产免费一区二区三区四区乱码| h日本视频在线播放| av国产久精品久网站免费入址| 人体艺术视频欧美日本| 美女cb高潮喷水在线观看| 九九在线视频观看精品| 丰满少妇做爰视频| 高清午夜精品一区二区三区| 一区二区三区乱码不卡18| 欧美日韩一区二区视频在线观看视频在线| av不卡在线播放| 精品人妻偷拍中文字幕| 免费黄色在线免费观看| av播播在线观看一区| 特大巨黑吊av在线直播| 亚洲欧美成人综合另类久久久| 久久精品国产a三级三级三级| 十分钟在线观看高清视频www | 久久久久国产精品人妻一区二区| 午夜福利在线在线| 干丝袜人妻中文字幕| 欧美日韩视频精品一区| 欧美日韩亚洲高清精品| 久久精品久久久久久噜噜老黄| 国产乱人偷精品视频| 国产欧美另类精品又又久久亚洲欧美| 国产淫语在线视频| 我要看日韩黄色一级片| 久久久久国产网址| 午夜福利视频精品| 欧美区成人在线视频| 女性被躁到高潮视频| 国产乱人视频| 妹子高潮喷水视频| 日韩伦理黄色片| 国产 一区 欧美 日韩| 欧美另类一区| 国产探花极品一区二区| 一级毛片我不卡| 成人毛片60女人毛片免费| 国产中年淑女户外野战色| 国产午夜精品久久久久久一区二区三区| 免费看日本二区| 久久99热这里只有精品18| 亚洲av成人精品一区久久| 多毛熟女@视频| 国产精品99久久久久久久久| 免费观看性生交大片5| 精品国产一区二区三区久久久樱花 | 欧美另类一区| 久久精品国产亚洲av涩爱| 18禁在线无遮挡免费观看视频| .国产精品久久| 国产片特级美女逼逼视频| 欧美亚洲 丝袜 人妻 在线| 在线亚洲精品国产二区图片欧美 | 欧美日本视频| 国精品久久久久久国模美| 九九在线视频观看精品| 下体分泌物呈黄色| 汤姆久久久久久久影院中文字幕| 晚上一个人看的免费电影| 亚洲av男天堂| 成年女人在线观看亚洲视频|