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

    基于剛?cè)峤佑|動力學(xué)的打結(jié)器虛擬打結(jié)方法研究

    2016-10-27 02:03:51尹建軍陳亞明
    農(nóng)業(yè)機械學(xué)報 2016年9期
    關(guān)鍵詞:模型

    尹建軍 高 強 陳亞明

    (江蘇大學(xué)現(xiàn)代農(nóng)業(yè)裝備與技術(shù)教育部重點實驗室, 鎮(zhèn)江 212013)

    ?

    基于剛?cè)峤佑|動力學(xué)的打結(jié)器虛擬打結(jié)方法研究

    尹建軍高強陳亞明

    (江蘇大學(xué)現(xiàn)代農(nóng)業(yè)裝備與技術(shù)教育部重點實驗室, 鎮(zhèn)江 212013)

    針對打結(jié)器結(jié)構(gòu)改進或參數(shù)調(diào)整后在設(shè)計階段難以評判優(yōu)劣而增加制造投入成本的問題,以打結(jié)器能否打出繩結(jié)為評判準則,提出了一種基于剛?cè)峤佑|動力學(xué)的打結(jié)器虛擬打結(jié)方法。該方法利用ADAMS的Bushing連接建立大變形柔性捆繩的動力學(xué)模型,通過對捆繩和剛性構(gòu)件之間設(shè)置多接觸約束和載荷施加,實現(xiàn)了完整的打結(jié)器虛擬打結(jié)仿真過程,研究了捆繩在拉力作用下成結(jié)的剛?cè)峤换チW(xué)行為。通過對打結(jié)器樣機的成結(jié)動作與虛擬打結(jié)的仿真結(jié)果進行對比分析,驗證了捆繩動力學(xué)模型與剛?cè)峤佑|動力學(xué)仿真模型的正確性,可作為評價設(shè)計優(yōu)劣、保證打結(jié)器一次性制造成功的有效方法。利用該虛擬打結(jié)方法分析了打結(jié)嘴與夾繩盤動作時序差φ對打結(jié)器成結(jié)的影響,對判定打結(jié)器參數(shù)匹配后能否成結(jié)或成結(jié)優(yōu)劣起到了直觀驗證作用。

    打捆機; 打結(jié)器; 接觸動力學(xué); 虛擬打結(jié); ADAMS

    引言

    D型打結(jié)器是方捆機捆扎裝置的核心部件,但我國方捆機的生產(chǎn)長期依賴進口打結(jié)器。為了簡化打結(jié)器結(jié)構(gòu),對打結(jié)器進行結(jié)構(gòu)改進或參數(shù)調(diào)整[1-3],一般采用運動分析結(jié)合虛擬裝配后即進入樣機試制階段,但改進后的打結(jié)器能否打出理想繩結(jié),缺乏直觀可視化的實現(xiàn)方法來驗證設(shè)計優(yōu)劣,未經(jīng)虛擬打結(jié)驗證設(shè)計優(yōu)劣的打結(jié)器試制,可能導(dǎo)致試制失敗而使制造投入面臨風(fēng)險。因此,有必要研究打結(jié)器虛擬打結(jié)實現(xiàn)方法來完善打結(jié)器設(shè)計理論。

    與打結(jié)器剛性執(zhí)行機構(gòu)的動作仿真[4-6]相比,由于捆繩的力學(xué)特性差異和柔性大變形問題,對打結(jié)過程的剛?cè)峤佑|動力學(xué)研究十分困難。國外對繩索類動力學(xué)研究主要基于質(zhì)子-彈簧模型對水下漁網(wǎng)運動過程進行了數(shù)值模擬[7-8]。國內(nèi)較多采用商用軟件對復(fù)雜繩索體運動過程進行模擬,比如:文獻[9]采用LS-DYNA 安全帶單元研究復(fù)雜繩索體動力學(xué)仿真建模方法,采用“節(jié)點-表面”接觸類型實現(xiàn)柔性網(wǎng)體與對象的接觸模擬。文獻[10]采用ANSYS和ADAMS軟件對多操作機排牙機器人系統(tǒng)中的牙弓曲線發(fā)生器進行聯(lián)合仿真。文獻[11]研究了虛擬手術(shù)中繩子打結(jié)的仿真實現(xiàn)。上述研究方法為打結(jié)器剛?cè)峤佑|動力學(xué)行為研究提供了思路和有效途徑。

    本文針對打結(jié)器結(jié)構(gòu)改進或參數(shù)調(diào)整后在設(shè)計階段難以評判優(yōu)劣而增加制造投入風(fēng)險的問題,基于ADAMS的Bushing模型建立反映捆繩柔性和力學(xué)特性的動力學(xué)模型,以打結(jié)器能否打出繩結(jié)為評判準則,提出打結(jié)器虛擬打結(jié)的實現(xiàn)方法,研究捆繩在拉力作用下打結(jié)過程的剛?cè)峤换チW(xué)行為,為打結(jié)器一次性制造成功提供評價設(shè)計優(yōu)劣的有效方法。

    1 捆繩動力學(xué)模型

    捆繩為大變形柔性體,而柔性體采用模態(tài)來表示物體的彈性,相對捆繩打結(jié)的大變形而言,捆繩的彈性是小變形。假設(shè)捆繩由若干微段圓柱組成,每個微段均視作剛體,相鄰兩個微段圓柱通過ADAMS軟件中的Bushing力連接,當微段圓柱相比捆繩總長度很小時,該微段剛體組合模型可視作捆繩的柔性體模型。在兩微段圓柱間施加Bushing力并使其相互作用,通過在兩微段圓柱中心連接處定義沿坐標系軸的3個分力和3個分力矩(Fx、Fy、Fz、Tx、Ty、Tz)獲得柔性力連接受力模型,如圖1所示。

    圖1 基于Bushing力連接的捆繩動力學(xué)模型Fig.1 Dynamic model of the rope based on Bushing connection

    根據(jù)力學(xué)理論,Bushing力[12]可表示為

    (1)

    式中K——捆繩剛度矩陣,為6×6對角矩陣

    C——捆繩阻尼系數(shù)矩陣,為6×6對角矩陣

    (Fx0,F(xiàn)y0,F(xiàn)z0)——沿坐標軸方向分力初始值(Tx0,Ty0,Tz0)——沿坐標軸方向分力矩初始值

    (rx,ry,rz)——第1微段圓柱中心處坐標系相對于第2微段圓柱中心處坐標系的相對位移

    (θx,θy,θz)——第1微段圓柱中心處坐標系相對于第2微段圓柱中心處坐標系的轉(zhuǎn)角

    (vx,vy,vz)——第1微段圓柱中心處坐標系相對于第2微段圓柱中心處坐標系的速度

    (ωx,ωy,ωz)——第1微段圓柱中心處坐標系相對于第2微段圓柱中心處坐標系的角速度

    (rx,ry,rz)、(θx,θy,θz)、(vx,vy,vz)、(ωx,ωy,ωz)由求解器自動確定。

    捆繩剛度矩陣對角元素計算式為

    式中K11——拉伸剛度

    K22、K33——剪切剛度

    K44——扭轉(zhuǎn)剛度

    K55、K66——彎曲剛度

    E——捆繩彈性模量

    G——捆繩剪切模量μ——泊松比

    I——捆繩中心慣性矩

    A、D、L——捆繩截面積、直徑和長度

    由式(1)可看出,相鄰兩微段圓柱間的柔性力大小與捆繩材料特性、相對位移、轉(zhuǎn)角、速度、角速度有關(guān)。在此,捆繩選直徑D=2.5 mm的聚丙烯繩,長度L=0.48 m。根據(jù)聚丙烯的力學(xué)特性[13],E=1.32 GPa,μ=0.35,經(jīng)計算得出G=0.49 GPa,I=1.92 mm4,K11=13 502.5 N/m,K22=K33=5 012.3 N/m,K44=3 920.0 N/m,K55=K66=5 280.0 N/m。根據(jù)聚丙烯的蠕變特性[13],阻尼系數(shù)矩陣C的對角線元素C11= C22= C33=0.71N·s/m,C44= C55=C66=0.1N·s/m。

    按上述參數(shù)在ADAMS中設(shè)定相鄰兩微段捆繩模型的剛度及阻尼系數(shù),使整段捆繩模型的力學(xué)性能接近捆繩真實力學(xué)性能,從而可模擬整段捆繩的變形及動力學(xué)性能。

    2 打結(jié)器剛?cè)峤佑|動力學(xué)建模

    2.1打結(jié)器成結(jié)裝置的虛擬樣機模型

    采用SolidWorks軟件建立打結(jié)器成結(jié)裝置模型,如圖2a所示??紤]到打結(jié)過程中撥繩盤的撥繩動作迫使捆繩搭接在打結(jié)嘴上,故將撥繩機構(gòu)一并建立到虛擬打結(jié)模型中。撥繩機構(gòu)由撥繩凸輪機構(gòu)與空間曲柄搖桿機構(gòu)串聯(lián)而成,其功能是撥繩盤將捆繩撥向打結(jié)嘴使捆繩搭接在打結(jié)嘴表面上,為打結(jié)嘴繞扣做準備。成結(jié)裝置模型中采用雙齒盤驅(qū)動打結(jié)器[14],其結(jié)構(gòu)組成如圖2b所示。將SolidWorks軟件下的打結(jié)器成結(jié)裝置模型轉(zhuǎn)換成parasolid圖形核心的 *.x_t格式,再將該中間文件導(dǎo)入到ADAMS軟件中,對導(dǎo)入的模型修改質(zhì)量特性以保持與實際材料特性相符,并刪除局部不影響仿真分析的零件,以簡化仿真模型。

    圖2 打結(jié)器成結(jié)裝置虛擬樣機模型Fig.2 Virtual prototyping model of knotting device of knotter1.雙齒盤驅(qū)動打結(jié)器 2.連桿 3.撥繩盤 4.滾輪 5.捆繩 6.機架 7.撥繩曲柄 8.擺桿 9.驅(qū)動主軸 10.撥繩凸輪 11.大齒盤 12.打結(jié)器支架 13.小齒盤 14.打結(jié)嘴錐齒輪 15.壓簧 16.壓板 17.割繩刀 18.打結(jié)嘴 19.鉤鉗 20.刀臂 21.夾繩盤 22.夾繩靜片 23.蝸桿軸錐齒輪 24.蝸桿 25.斜齒輪 26.夾繩動片 27.夾繩壓板

    為了檢驗設(shè)計的雙齒盤驅(qū)動打結(jié)器是否可以實現(xiàn)捆繩打結(jié),在圖2所示的虛擬樣機模型中建立捆繩模型。在此,采用240個直徑2.5mm、長度2mm的小圓柱體通過Bushing力相互連接起來形成整根捆繩,按照送繩結(jié)束之后捆繩的形態(tài)進行裝配,即捆繩一端搭接在夾繩盤槽口,捆繩另一端與滑輪形成接觸。

    2.2打結(jié)器成結(jié)裝置仿真模型設(shè)置

    2.2.1驅(qū)動及約束設(shè)置

    打結(jié)器驅(qū)動主軸的額定轉(zhuǎn)速為90r/min,故對驅(qū)動主軸施加旋轉(zhuǎn)驅(qū)動為540°/s。根據(jù)打結(jié)器成結(jié)裝置構(gòu)件之間的實際連接關(guān)系在ADAMS下添加相應(yīng)的運動副,如表1所示。對于錐齒輪傳動和蝸桿傳動、捆繩和剛性構(gòu)件均存在接觸關(guān)系,將其運動副均設(shè)置為“solidtosolid”接觸連接。捆繩的240個小圓柱之間均添加Bushing力連接,Bushing設(shè)置分為兩部分,平移運動部分按捆繩動力學(xué)模型中的K11、K22、K33、C11、C22、C33數(shù)值進行設(shè)定,旋轉(zhuǎn)運動部分按捆繩動力學(xué)模型中的K44、K55、K66、C44、C55、C66數(shù)值進行設(shè)定。

    對于多接觸約束機械系統(tǒng)的虛擬樣機模型,ADAMS求解器解算一組微分-代數(shù)方程構(gòu)成的動力學(xué)方程組[15],通常采用數(shù)值迭代法。根據(jù)ADAMS采用的接觸判斷準則和接觸動力學(xué)模型,求解動力學(xué)方程組時,ADAMS求解器將接觸力并入廣義力矩陣中求解。對不同構(gòu)件之間的接觸作用,需要通過定義接觸表面,設(shè)定接觸類型來實現(xiàn)。為此,為了避免捆繩與剛性構(gòu)件之間發(fā)生穿透現(xiàn)象,須將捆繩和與打結(jié)器中有接觸關(guān)系的構(gòu)件之間添加接觸副,在本模型中將打結(jié)嘴、勾鉗、刀臂、夾繩動片、夾繩盤、

    表1 構(gòu)件之間的約束關(guān)系設(shè)置Tab.1 Constraint relation configuration among components

    撥繩盤與捆繩模型的小圓柱之間添加接觸連接。ADAMS接觸連接對話框中提供了接觸零件的動摩擦因數(shù)和靜摩擦因數(shù)的參數(shù)設(shè)置,為了提高仿真的真實性,必須對需要接觸的零件的動摩擦因數(shù)和靜摩擦因數(shù)進行設(shè)置。由于實際打結(jié)中采用聚丙烯繩,添加捆繩與打結(jié)嘴、勾鉗、夾繩盤、夾繩刀、割繩刀臂的靜摩擦因數(shù)為0.37、動摩擦因數(shù)為0.18[16]。同時,捆繩各段圓柱體之間也須添加接觸連接,避免打結(jié)過程中各段圓柱之間發(fā)生穿透而無法形成繩結(jié)。由于有接觸關(guān)系的各段圓柱之間只占整根捆繩的一小段,故只需將捆繩模型中有接觸關(guān)系的各段圓柱之間添加接觸,可以大大提高ADAMS求解器的計算效率。

    2.2.2載荷施加

    如圖3所示,利用室內(nèi)臺架試驗測量打結(jié)過程中捆繩拉力F變化曲線,測量方法如下:

    在初級捆繩和次級捆繩上分別串聯(lián)一個拉力傳感器,手動驅(qū)動飛輪使打結(jié)器主軸轉(zhuǎn)動進行打空結(jié),通過拉力采集計算機系統(tǒng)得到成結(jié)過程的捆繩拉力,如圖4所示。

    圖3 捆繩拉力原理與測量Fig.3 Pulling force measurement and principle of rope1.打結(jié)嘴 2.鉤鉗 3.初級捆繩拉力傳感器 4.初級捆繩 5.次級捆繩 6.次級捆繩拉力傳感器 7.次級捆繩拉力傳感器顯示器 8.初級捆繩拉力傳感器顯示器 9.拉力采集計算機系統(tǒng)

    圖4 捆繩拉力變化曲線Fig.4 Pulling force curves of rope

    由于虛擬打結(jié)過程是單股繩,而實際打結(jié)過程是雙股繩。在仿真時取初級捆繩和次級捆繩拉力之和的峰值模擬捆繩的實際受力,由圖4拉力變化曲線可知該峰值取為200 N。據(jù)此,在捆繩模型末端沿端面法向添加2個角度固定、大小為200 N的拉力。為了模擬打結(jié)器在成結(jié)后捆繩被刀臂上的刀片切斷的情況,在仿真模型的夾繩盤處捆繩末端添加一個STEP函數(shù)驅(qū)動的拉力,STEP函數(shù)表示如下:

    F=STEP(time,0,200,0.48,200)+

    STEP(time,0.48,0,0.5,-200)

    設(shè)置該力的意義是:在割繩刀接觸到打結(jié)嘴與夾繩盤之間的捆繩時,夾繩盤處捆繩末端拉力消失,可近似模擬捆繩被割繩刀片割斷。通過多次仿真嘗試設(shè)定時間t=0.5 s,割繩刀接觸到打結(jié)嘴與夾繩盤之間的捆繩。

    對式(1)中捆繩各段圓柱之間的分力和分力矩初始值Fx0、Fy0、Fz0、Tx0、Ty0、Tz0分別取為200 N、0 N、0 N和0 N·m、0 N·m、0 N·m。

    3 打結(jié)器虛擬打結(jié)與試驗驗證

    3.1虛擬打結(jié)仿真

    經(jīng)過上述驅(qū)動及約束設(shè)置和載荷施加后,利用ADAMS求解器對打結(jié)器成結(jié)裝置虛擬樣機模型進行仿真,其仿真結(jié)果如圖5所示。

    圖5 雙齒盤驅(qū)動打結(jié)器打結(jié)仿真過程圖Fig.5 Knotting simulation process drawings of the knotter driven by double gear-discs

    為降低仿真計算量和難度,仿真時捆繩為單股,而實際打結(jié)是雙股,這并不影響虛擬打結(jié)的正確性。從圖5a到圖5i中的成結(jié)仿真過程可見,成結(jié)裝置的虛擬樣機模型中多接觸約束設(shè)置合適,捆繩的動力學(xué)模型建模方法可行,參數(shù)設(shè)置正確,可以模擬捆繩的大變形及動力學(xué)行為,能再現(xiàn)捆繩在拉力作用下與剛性構(gòu)件之間的接觸約束過程。利用虛擬打結(jié)可實時觀察捆繩在拉力作用下形成繩結(jié)的剛?cè)峤换ミ^程,包括撥繩盤將捆繩撥向打結(jié)嘴搭接位置的動態(tài)過程、捆繩被刀臂的彎曲弧面導(dǎo)入弧形槽口、捆繩被夾繩器夾持旋轉(zhuǎn)、捆繩被打結(jié)嘴鉸扭成環(huán)、鉤鉗張開鉗咬捆繩、鉤鉗閉合鉤出繩頭和刀臂強制將繩扣從打結(jié)嘴上推脫。對于捆繩模型,文獻[4]僅考慮了捆繩剛度系數(shù),且剛度系數(shù)矩陣為定值,未考慮捆繩阻尼系數(shù)矩陣,而本文建立的捆繩動力學(xué)模型接近捆繩真實力學(xué)性能,可獲得較為準確的整段捆繩的變形及動力學(xué)性能。

    在虛擬打結(jié)仿真中,若夾繩盤處捆繩末端被固定,滾輪側(cè)捆繩末端沿端面法向施加拉力,仿真結(jié)果如圖6所示,這與實際打結(jié)過程不符,其原因是ADAMS求解器難以仿真捆繩被刀片割斷的過程,捆繩一端被夾繩器夾持使得絞扭成環(huán)的繩扣纏繞在打結(jié)嘴上不能被刀臂強制推脫,而是繩扣繞打結(jié)嘴表面反卷到鉤鉗上,直至繩扣散開,打結(jié)仿真失敗。這也是在夾繩盤處捆繩末端添加一個STEP函數(shù)驅(qū)動拉力的原因,來近似模擬刀片割繩動作。文獻[4]通過ADAMS二次開發(fā)程序控制微段圓柱運動來近似模擬割繩動作,對ADAMS二次開發(fā)要求高。而本文模擬割繩動作的方法更具有通用性,便于實現(xiàn)。

    圖6 虛擬打結(jié)的割繩動作存在的問題Fig.6 Problem of cutting rope motion of virtual knotting

    3.2虛擬打結(jié)的試驗驗證

    將設(shè)計的雙齒盤驅(qū)動打結(jié)器安裝到室內(nèi)打結(jié)臺架上,試驗時調(diào)整捆繩張力、鉤鉗壓簧和夾繩壓板的壓力達到合適值,在不上秸稈物料的情況下,慢速驅(qū)動打結(jié)器做打結(jié)動作,并用數(shù)碼相機記錄打結(jié)器成結(jié)過程,與虛擬打結(jié)仿真結(jié)果進行對比,如圖7所示,其中左側(cè)圖為虛擬樣機打結(jié),右側(cè)圖為打結(jié)器樣機打結(jié)。

    圖7 虛擬樣機與打結(jié)器樣機打結(jié)過程對比Fig.7 Contrast of knotting process between virtual prototype and knotter prototype

    由圖7可見,虛擬打結(jié)與打結(jié)器物理樣機的打結(jié)過程非常吻合,表明該虛擬打結(jié)實現(xiàn)方法可以用于設(shè)計階段檢驗打結(jié)器的結(jié)構(gòu)改進和參數(shù)調(diào)整后的成結(jié)效果,保證打結(jié)器物理樣機一次性制造的成功,大大降低了打結(jié)器制造的投入成本。

    3.3虛擬打結(jié)應(yīng)用

    打結(jié)嘴和夾繩盤動作時序差φ和打結(jié)嘴軸線與主軸軸線夾角α是打結(jié)器設(shè)計的兩個基本參數(shù)[1]。利用上述方法可對φ和α在不同組合下的打結(jié)器模型進行虛擬打結(jié),以驗證參數(shù)匹配是否合適。在此,α取98°,φ取2.5°、7.5°、12.5°、17.5°、22.5°共5組值。如圖8所示,蝸桿軸錐齒輪與打結(jié)嘴錐齒輪安裝相位差ψ取為30°,小齒盤上的不完全錐齒輪與大齒盤上的不完全錐齒輪的安裝相位差按ε=φ+ψ計算[14],當φ取2.5°、7.5°、12.5°、17.5°、22.5°時,ε取32.5°、37.5°、42.5°、47.5°、52.5°。

    圖8 齒盤錐齒輪安裝相位差示意圖Fig.8 Installation phase difference between bevel gears on the gear-discs1.小齒盤 2.大齒盤 3.打結(jié)嘴軸心線 4.夾繩盤蝸桿軸心線

    修改打結(jié)器虛擬樣機模型使小齒盤上的不完全錐齒輪與大齒盤上的不完全錐齒輪的安裝相位差ε分別為32.5°、37.5°、47.5°、52.5°,利用前述方法對其虛擬打結(jié),并與φ取12.5°時的虛擬打結(jié)過程對比,如圖9~13所示。

    圖9 φ取2.5°時的虛擬打結(jié)Fig.9 Virtual knotting when φ is 2.5°

    圖10 φ取7.5°時的虛擬打結(jié)Fig.10 Virtual knotting when φ is 7.5°

    圖11 φ取12.5°時的虛擬打結(jié)Fig.11 Virtual knotting when φ is 12.5°

    圖12 φ取17.5°時的虛擬打結(jié)Fig.12 Virtual knotting when φ is 17.5°

    由圖9~13的圖a可見,當φ為2.5°、7.5°、12.5°、17.5°、22.5°時,繞扣開始時夾繩盤相對打結(jié)嘴的轉(zhuǎn)動逐漸超前,隨著打結(jié)嘴與夾繩盤時序差φ增大,繞扣開始時捆繩沿夾繩盤逆時針方向捆繩夾持點向夾繩動片偏移,逐漸處于夾繩盤與夾繩動片形成的夾持閉合空間。當φ=2.5°時,繞扣開始時夾繩盤未將捆繩送入夾繩動片夾緊,開始繞扣后捆繩被繩針帶回,虛擬打結(jié)失敗,表明小齒盤上的不完全錐齒輪與蝸桿軸錐齒輪嚙合滯后大齒盤上的不完全錐齒輪與打結(jié)嘴軸錐齒輪嚙合太多,即φ取2.5°過小。當φ=7.5°時,繞扣開始時夾繩盤并未將捆繩完全送入夾繩動片,但在繩針夾持捆繩的輔助作用下,隨著打結(jié)嘴旋轉(zhuǎn)繞扣的同時捆繩被逐漸夾緊,不影響虛擬打結(jié)。當φ=12.5°時,繞扣開始時捆繩恰好被夾繩盤完全送入夾繩動片夾緊。當φ=17.5°和φ=22.5°時,繞扣開始時捆繩已被夾繩盤與夾繩動片夾緊。

    由圖10c、10d可見,當φ=7.5°時,捆繩處于鉤鉗末端與打結(jié)嘴末端連線的中間偏上位置,與勾鉗距離較近,鉤鉗可咬住捆繩;由圖12c~12e可見,當φ=17.5°時,捆繩處于鉤鉗末端與打結(jié)嘴末端連線的中間偏下位置,貼近打結(jié)嘴上表面,鉤鉗可咬住捆繩;由圖11c~11e可見,當φ=12.5°時,捆繩處于鉤鉗末端與打結(jié)嘴末端連線的中間位置,鉤鉗可靠咬住捆繩,表明φ=12.5°可獲得鉗咬捆繩最佳位置;而當φ=22.5°時,捆繩在夾繩盤上的夾持點過低,鉗咬時捆繩滑動到打結(jié)嘴側(cè)方,導(dǎo)致鉤鉗無法鉗咬住捆繩,鉗咬動作失敗而未能成結(jié),如圖13c、13d所示,表明φ取22.5°過大。

    由上述分析,當α=98°、φ為7.5°~17.5°時,可以實現(xiàn)虛擬打結(jié),表明小齒盤上的不完全錐齒輪

    圖13 φ取22.5°時的虛擬打結(jié)Fig.13 Virtual knotting when φ is 22.5°

    與大齒盤上的不完全錐齒輪的安裝相位差允許在(42.5±5)°間變動,約為半個齒的角度。

    從圖9~13可清楚看出打結(jié)嘴和夾繩盤動作時序差φ取不同值時,虛擬打結(jié)效果比解析計算后的數(shù)值模擬結(jié)果更加直觀,對判定打結(jié)器參數(shù)匹配是否合適起到直觀驗證作用,有助于形成新的打結(jié)器設(shè)計方案。利用該方法能解決打結(jié)器設(shè)計后難以直觀驗證成結(jié)效果的問題,在設(shè)計階段可發(fā)現(xiàn)設(shè)計缺陷,降低打結(jié)器制造的投入成本。

    4 結(jié)論

    (1)建立了完整的打結(jié)器剛?cè)峤佑|虛擬打結(jié)實現(xiàn)方法,通過對打結(jié)器樣機的成結(jié)動作與虛擬打結(jié)的仿真結(jié)果進行對比分析,驗證了捆繩動力學(xué)模型與剛?cè)峤佑|動力學(xué)仿真的正確性,可作為評價設(shè)計優(yōu)劣、保證打結(jié)器一次性制造成功的有效方法。

    (2)通過對雙齒盤打結(jié)器打結(jié)過程的剛?cè)峤佑|動力學(xué)仿真,驗證了ADAMS的Bushing連接可解決大變形柔性繩索類模型的建立與仿真問題,對剛?cè)峤佑|的機構(gòu)設(shè)計與分析具有借鑒價值。

    1尹建軍, 張萬慶, 陳亞明, 等. 打結(jié)器夾繩-繞扣-鉗咬動作參數(shù)分析與打結(jié)試驗[J]. 農(nóng)業(yè)機械學(xué)報, 2015,46(9):135-143.

    YIN Jianjun, ZHANG Wanqing, CHEN Yaming, et al. Parameters analysis of rope-holding motion, knot-winding motion, rope-biting motion of knotter and knotting tests[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(9):135-143. (in Chinese)

    2LI Hui, HE Jin, WANG Qingjie, et al. Design and experiment of φ-type-knots knotters on Chinese small square balers[J]. Chinese Journal of Mechanical Engineering: English Edition, 2014, 27(1): 154-164.

    3張紹英, 李海濤, 曹慶和, 等. 雙α結(jié)打結(jié)器關(guān)鍵傳動機構(gòu)設(shè)計[J]. 農(nóng)業(yè)機械學(xué)報, 2013, 44(12): 74-79.

    ZHANG Shaoying, LI Haitao, CAO Qinghe, et al. Design of key transmission mechanism of double-α-knot knotter[J]. Transactions of the Chinese Society for Agricultural Machinery,2013, 44(12): 74-79. (in Chinese)

    4尹建軍, 李雙, 李耀明. D 型打結(jié)器及其輔助機構(gòu)運動仿真與時序分析[J]. 農(nóng)業(yè)機械學(xué)報, 2011, 42(6): 103-107.

    YIN Jianjun, LI Shuang, LI Yaoming. Kinematic simulation and time series analysis of D-knotter and its ancillary mechanisms[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(6): 103-107. (in Chinese)

    5陳龍健, 李誠, 張安琪, 等. 基于逆向工程的D 型打結(jié)器重構(gòu)與運動仿真[J]. 農(nóng)業(yè)機械學(xué)報, 2014,45(12):104-108.

    CHEN Longjian, LI Cheng, ZHANG Anqi, et al. Reconstitution and motion simulation of D-bale knotter based on reverse engineering[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014,45(12):104-108. (in Chinese)

    6黃偉,張志毅,劉志強,等. D型打結(jié)器Pro/Engineer 仿真及動作功能分析[J]. 農(nóng)機化研究, 2009,31(9):58-60.

    HUANG Wei, ZHANG Zhiyi, LIU Zhiqiang, et al. Pro/Engineer imitation and the analysis on movements and functions of twine knott[J]. Journal of Agricultural Mechanization Research, 2009,31(9):58-60. (in Chinese)

    7BESSONNEAU J S, MARICHAL D. Study of the dynamics of submerged supple nets (applications to trawls) [J]. Ocean Engineering, 1998, 25(7): 563-583.

    8SHARF Inna, ZHANG Yunning. A contact force solutionfor non-colliding contact dynamics simulation[C]∥ASME 2006 International Design & Engineering Technical Conferences and Computers & Information in Engineering Conference, 2006: 677-687.

    9趙國偉, 朱鷗寧, 徐云飛, 等. LS-DYNA 帶單元在繩索體仿真中的應(yīng)用研究[J]. 機械科學(xué)與技術(shù), 2013, 32(9): 1249-1253.

    ZHAO Guowei, ZHU Ouning, XU Yunfei, et al. The application of LS-DYNA seatbelt element in net simulation[J]. Mechanical Science and Technology, 2013, 32(9): 1249-1253. (in Chinese)

    10張永德, 汪洋濤, 王沫楠,等. 基于ANSYS 與ADAMS 的柔性體聯(lián)合仿真[J]. 系統(tǒng)仿真學(xué)報, 2008, 20(17): 4501-4504.

    ZHANG Yongde, WANG Yangtao, WANG Monan, et al.Co-simulation of flexible body based on ANSYS and ADAMS [J]. Journal of System Simulation, 2008, 20(17): 4501-4504. (in Chinese)

    11賈世宇, 潘振寬. 具備觸覺反饋的實時繩子打結(jié)仿真[J]. 系統(tǒng)仿真學(xué)報, 2009, 21(20): 6519-6523.

    JIA Shiyu, PAN Zhenkuan. Real-time knot tying simulation of ropes with haptics[J]. Journal of System Simulation,2009, 21(20): 6519-6523. (in Chinese)

    12KLISCH T. Contact mechanics in multibody systems [J]. Mechanism and Machine Theory,1999, 34(5): 665-675.

    13周魯豫. UHMWPE纖維繩索的編織工藝和力學(xué)性能研究[D]. 青島: 青島大學(xué), 2013.

    ZHOU Luyu. Study on the knitting process and mechanical properties of UHMWPE fiber rope[D]. Qingdao: Qingdao University, 2013. (in Chinese)

    14尹建軍,張萬慶,陳亞明, 等. 雙齒盤驅(qū)動打結(jié)器設(shè)計與成結(jié)試驗分析[J]. 農(nóng)業(yè)機械學(xué)報,2016,47(3):98-105.

    YIN Jianjun, ZHANG Wanqing, CHEN Yaming, et al.Design and knotting test analysis of knotter driven by double gear-discs [J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3):98-105.(in Chinese)

    15石明全, 張鵬, 范樹遷,等. 基于多接觸約束的動力學(xué)仿真計算效率[J]. 機械工程學(xué)報, 2011, 46(7): 108-113.

    SHI Mingquan, ZHANG Peng, FAN Shuqian,et al. Computational efficiency for dynamics simulation based on multi-contact mechanical system[J]. Journal of Mechanical Engineering, 2011, 46(7): 108-113. (in Chinese)

    16陳龍健, 李誠, 張安琪, 等. 秸稈捆扎過程中打結(jié)鉗嘴載荷試驗分析[J]. 農(nóng)業(yè)機械學(xué)報, 2015, 46(9): 128-134.

    CHEN Longjian, LI Cheng, ZHANG Anqi, et al. Load experimental analysis of bill hook during straw baling[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(9): 128-134. (in Chinese)

    Virtual Knotting Method of Knotter Based on Rigid-flexible Contact Dynamics

    Yin JianjunGao QiangChen Yaming

    (KeyLaboratoryofModernAgriculturalEquipmentandTechnology,MinistryofEducation,JiangsuUniversity,Zhenjiang212013,China)

    Aiming at the problem of increasing investment cost of knotter production due to difficulty in evaluation of good or bad design after structure improvement or parameters adjustment of the knotter, a kind of virtual knotting method of knotter based on rigid-flexible contact dynamics was presented according to the judgment rule of whether the knotter could tie a knot during design phase. The method used the Bushing connection tool under ADAMS to build up dynamic model of flexible rope with large deformation. A complete virtual knotting simulation process of the knotter was implemented by setting up multiple contact constraint between the rope and rigid component of the knotter and applying load to the rope. Rigid-flexible interactive mechanical behaviors of tying a knot of the knotter were studied under the action of pulling force of the rope. The comparison analyses between knotting movements of the knotter prototype and virtual knotting results validated the correctness of dynamic model of the rope and simulation model of rigid-flexible contact dynamics. It may be regarded as an effective method to evaluate good or bad design and guarantee disposable successful production of the knotter. The effects of the timing sequence difference valueφbetween the motion of knotting hook and the motion of rope-holding plate on knotting motion of the knotter were analyzed by using the proposed virtual knotting method, which can play visual validation role to judge whether the knotter can knot or tie a good and bad knot after parameters match of the knotter.

    baler; knotter; contact dynamics; virtual knotting; ADAMS

    10.6041/j.issn.1000-1298.2016.09.013

    2016-02-19

    2016-03-18

    國家自然科學(xué)基金項目(51375215)、鎮(zhèn)江市重點研發(fā)計劃項目(NY2015009)和江蘇高校優(yōu)勢學(xué)科建設(shè)工程項目(蘇政辦發(fā)(2014)37號)

    尹建軍(1973—),男,研究員,博士,主要從事收獲機械集成設(shè)計和農(nóng)業(yè)機器人技術(shù)研究,E-mail: yinjianjun@ujs.edu.cn

    S491

    A

    1000-1298(2016)09-0085-08

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    黄片小视频在线播放| 啦啦啦在线免费观看视频4| 七月丁香在线播放| 少妇人妻 视频| 香蕉国产在线看| 免费女性裸体啪啪无遮挡网站| 最近中文字幕2019免费版| 你懂的网址亚洲精品在线观看| 亚洲熟女精品中文字幕| 免费在线观看日本一区| 国产日韩欧美视频二区| 久久毛片免费看一区二区三区| 在线看a的网站| 一区二区av电影网| 国产av国产精品国产| 午夜两性在线视频| 人人妻人人爽人人添夜夜欢视频| 欧美精品一区二区免费开放| 色视频在线一区二区三区| 亚洲欧美日韩另类电影网站| 欧美激情 高清一区二区三区| av网站免费在线观看视频| 十分钟在线观看高清视频www| 中文字幕av电影在线播放| 国产无遮挡羞羞视频在线观看| 国产午夜精品一二区理论片| 欧美亚洲 丝袜 人妻 在线| 欧美黄色淫秽网站| 欧美日韩国产mv在线观看视频| 国产成人欧美| 久久鲁丝午夜福利片| 在线观看免费日韩欧美大片| 黑人欧美特级aaaaaa片| 美女视频免费永久观看网站| 大香蕉久久网| 亚洲欧美成人综合另类久久久| 亚洲精品久久午夜乱码| 一级黄片播放器| 欧美97在线视频| 黄频高清免费视频| av网站在线播放免费| 高清欧美精品videossex| 精品人妻1区二区| 女人久久www免费人成看片| 国产精品久久久久久人妻精品电影 | 久久久久久久久免费视频了| 国产成人欧美| 成人亚洲精品一区在线观看| e午夜精品久久久久久久| 国产成人a∨麻豆精品| 纯流量卡能插随身wifi吗| 午夜免费男女啪啪视频观看| 90打野战视频偷拍视频| 久久久国产精品麻豆| 亚洲国产中文字幕在线视频| 少妇人妻 视频| 99热国产这里只有精品6| 在线av久久热| 免费久久久久久久精品成人欧美视频| 又粗又硬又长又爽又黄的视频| 男女边吃奶边做爰视频| 久久影院123| 在线观看免费高清a一片| 国产在视频线精品| 99热全是精品| 中文字幕亚洲精品专区| 巨乳人妻的诱惑在线观看| 1024香蕉在线观看| 99精国产麻豆久久婷婷| 成人国产av品久久久| 久久久久久久久久久久大奶| 国产片内射在线| 亚洲五月色婷婷综合| 在线观看免费视频网站a站| 久久久久久久久久久久大奶| 久久99热这里只频精品6学生| 国产精品久久久久久精品古装| 亚洲,欧美,日韩| 两性夫妻黄色片| 啦啦啦中文免费视频观看日本| 18禁观看日本| 中文字幕人妻丝袜制服| 激情五月婷婷亚洲| 美女视频免费永久观看网站| 久久中文字幕一级| 大码成人一级视频| 国产精品久久久久久精品电影小说| 美女中出高潮动态图| 男女下面插进去视频免费观看| 国产精品免费大片| 热re99久久精品国产66热6| 人妻一区二区av| 国产精品一国产av| 久久精品国产综合久久久| 男女午夜视频在线观看| 黄频高清免费视频| 国产精品一区二区免费欧美 | 亚洲男人天堂网一区| 久久九九热精品免费| 高潮久久久久久久久久久不卡| 亚洲国产欧美日韩在线播放| 免费观看av网站的网址| 久久久久精品国产欧美久久久 | 伊人久久大香线蕉亚洲五| 在线天堂中文资源库| 69精品国产乱码久久久| 亚洲精品久久成人aⅴ小说| 欧美大码av| 日日夜夜操网爽| 亚洲专区国产一区二区| 韩国精品一区二区三区| 天天影视国产精品| 精品一品国产午夜福利视频| 中文字幕人妻丝袜一区二区| 精品国产一区二区三区久久久樱花| 美女大奶头黄色视频| 日本午夜av视频| av片东京热男人的天堂| 国产亚洲欧美精品永久| 国产精品一区二区在线不卡| 国产免费视频播放在线视频| 性色av一级| 少妇精品久久久久久久| 国产91精品成人一区二区三区 | 伦理电影免费视频| 久久天躁狠狠躁夜夜2o2o | 久久国产亚洲av麻豆专区| 婷婷色综合www| 日韩中文字幕欧美一区二区 | 亚洲av男天堂| 激情五月婷婷亚洲| 亚洲精品日本国产第一区| 亚洲成人免费av在线播放| 人妻 亚洲 视频| 亚洲成av片中文字幕在线观看| av有码第一页| 天天添夜夜摸| 免费高清在线观看视频在线观看| 一区二区av电影网| 国产精品 欧美亚洲| 国产99久久九九免费精品| 欧美性长视频在线观看| 午夜精品国产一区二区电影| 一区福利在线观看| 久久久国产精品麻豆| 超色免费av| 精品久久蜜臀av无| 50天的宝宝边吃奶边哭怎么回事| 精品福利观看| 少妇精品久久久久久久| 九草在线视频观看| 在线观看免费午夜福利视频| 老司机在亚洲福利影院| 国产亚洲欧美在线一区二区| 2018国产大陆天天弄谢| 在线亚洲精品国产二区图片欧美| 亚洲av日韩在线播放| 国产人伦9x9x在线观看| 女人高潮潮喷娇喘18禁视频| 国产成人影院久久av| 国产av国产精品国产| tube8黄色片| 久久久久久久久久久久大奶| 精品卡一卡二卡四卡免费| 国产亚洲一区二区精品| 91麻豆精品激情在线观看国产 | 国产男女超爽视频在线观看| 黄色一级大片看看| 久久精品国产亚洲av高清一级| 黄色片一级片一级黄色片| 国产精品 欧美亚洲| 免费在线观看视频国产中文字幕亚洲 | 日日爽夜夜爽网站| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲熟女毛片儿| 国产真人三级小视频在线观看| 老司机影院毛片| 美女主播在线视频| 亚洲精品久久午夜乱码| 极品少妇高潮喷水抽搐| 久久精品成人免费网站| 一级毛片我不卡| 午夜免费观看性视频| 国产精品成人在线| 不卡av一区二区三区| 深夜精品福利| 国产亚洲午夜精品一区二区久久| 乱人伦中国视频| 国产免费福利视频在线观看| 精品一区在线观看国产| 十八禁高潮呻吟视频| 男女下面插进去视频免费观看| 一级毛片女人18水好多 | 成人影院久久| 久久久国产欧美日韩av| 欧美精品一区二区免费开放| 9热在线视频观看99| 男女无遮挡免费网站观看| av有码第一页| 亚洲专区中文字幕在线| 嫩草影视91久久| 欧美变态另类bdsm刘玥| 久久久国产精品麻豆| 国产黄频视频在线观看| www.精华液| 亚洲欧美清纯卡通| 精品人妻一区二区三区麻豆| 亚洲中文av在线| 免费看十八禁软件| 亚洲精品国产一区二区精华液| 午夜激情av网站| 欧美精品高潮呻吟av久久| 日韩,欧美,国产一区二区三区| 国产成人影院久久av| 多毛熟女@视频| 欧美日韩精品网址| videos熟女内射| 亚洲av电影在线进入| 九色亚洲精品在线播放| 人妻 亚洲 视频| 亚洲av成人精品一二三区| 亚洲欧美一区二区三区国产| 国产有黄有色有爽视频| 一区二区三区四区激情视频| 80岁老熟妇乱子伦牲交| 久久性视频一级片| 咕卡用的链子| 捣出白浆h1v1| 久久久精品国产亚洲av高清涩受| 我要看黄色一级片免费的| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 久久久久久免费高清国产稀缺| 大话2 男鬼变身卡| 午夜91福利影院| 亚洲av美国av| www.精华液| 亚洲av成人不卡在线观看播放网 | 亚洲人成电影免费在线| 久热爱精品视频在线9| 另类亚洲欧美激情| 99热网站在线观看| 午夜视频精品福利| 日韩一本色道免费dvd| 精品少妇内射三级| www.999成人在线观看| 91麻豆av在线| 无遮挡黄片免费观看| 亚洲天堂av无毛| 国产老妇伦熟女老妇高清| 亚洲视频免费观看视频| 国产精品三级大全| 国产又爽黄色视频| 亚洲免费av在线视频| 国产成人欧美| 中文欧美无线码| 日本欧美国产在线视频| 国产熟女欧美一区二区| 国产片内射在线| av天堂在线播放| 青草久久国产| 久久久久国产精品人妻一区二区| 天天操日日干夜夜撸| 国产成人a∨麻豆精品| 1024香蕉在线观看| 久久久精品94久久精品| 国产真人三级小视频在线观看| 少妇人妻 视频| 丰满饥渴人妻一区二区三| 51午夜福利影视在线观看| 国产免费视频播放在线视频| 欧美日本中文国产一区发布| 在线观看国产h片| 国产精品麻豆人妻色哟哟久久| 日本色播在线视频| 亚洲精品第二区| 少妇裸体淫交视频免费看高清 | 9热在线视频观看99| 黑丝袜美女国产一区| 看免费av毛片| 欧美激情 高清一区二区三区| 啦啦啦在线观看免费高清www| 91字幕亚洲| 超色免费av| 丝袜脚勾引网站| 成人亚洲欧美一区二区av| 亚洲情色 制服丝袜| 亚洲国产看品久久| 中文精品一卡2卡3卡4更新| 亚洲av在线观看美女高潮| 免费高清在线观看视频在线观看| 久久久久国产精品人妻一区二区| 久久免费观看电影| 欧美国产精品va在线观看不卡| 女人久久www免费人成看片| 久久久久久久久免费视频了| 亚洲,欧美精品.| 又粗又硬又长又爽又黄的视频| 精品人妻一区二区三区麻豆| 观看av在线不卡| 建设人人有责人人尽责人人享有的| 高清黄色对白视频在线免费看| 精品国产国语对白av| www.999成人在线观看| 男女国产视频网站| 亚洲人成网站在线观看播放| 欧美另类一区| 精品人妻在线不人妻| 国产精品人妻久久久影院| 亚洲伊人色综图| 亚洲熟女精品中文字幕| 婷婷成人精品国产| 国产色视频综合| 老司机在亚洲福利影院| av网站免费在线观看视频| 99国产精品一区二区三区| 成人亚洲欧美一区二区av| 日日摸夜夜添夜夜爱| a级片在线免费高清观看视频| 成人三级做爰电影| 国产日韩欧美在线精品| 两人在一起打扑克的视频| 黄色 视频免费看| 亚洲精品日韩在线中文字幕| 一级黄片播放器| 国产女主播在线喷水免费视频网站| 夜夜骑夜夜射夜夜干| 性高湖久久久久久久久免费观看| 亚洲 欧美一区二区三区| 宅男免费午夜| 黑人欧美特级aaaaaa片| 国产精品二区激情视频| 欧美大码av| 欧美变态另类bdsm刘玥| 欧美精品av麻豆av| av视频免费观看在线观看| 91麻豆av在线| 日韩精品免费视频一区二区三区| 成人三级做爰电影| 午夜免费观看性视频| 欧美日韩亚洲国产一区二区在线观看 | 久久精品人人爽人人爽视色| 99久久综合免费| 欧美日本中文国产一区发布| 久久99热这里只频精品6学生| 男女床上黄色一级片免费看| 一级毛片女人18水好多 | av电影中文网址| 成在线人永久免费视频| 日韩免费高清中文字幕av| 人体艺术视频欧美日本| 国语对白做爰xxxⅹ性视频网站| 大陆偷拍与自拍| 亚洲中文日韩欧美视频| 夫妻性生交免费视频一级片| 亚洲av片天天在线观看| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 精品免费久久久久久久清纯 | 只有这里有精品99| 国产精品久久久久久精品古装| 国产色视频综合| 在线观看www视频免费| 最新的欧美精品一区二区| 亚洲欧美一区二区三区国产| 最近中文字幕2019免费版| 久久ye,这里只有精品| 嫩草影视91久久| 久久免费观看电影| 人人妻人人添人人爽欧美一区卜| 新久久久久国产一级毛片| 亚洲av国产av综合av卡| av天堂久久9| 亚洲成人免费av在线播放| 日韩大码丰满熟妇| 黑人猛操日本美女一级片| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久人妻精品电影 | 精品一区二区三卡| 国产成人91sexporn| 又粗又硬又长又爽又黄的视频| 久久久久久人人人人人| 美女大奶头黄色视频| av国产久精品久网站免费入址| 国产欧美日韩一区二区三 | 国产成人影院久久av| 免费在线观看视频国产中文字幕亚洲 | 成年人免费黄色播放视频| 精品第一国产精品| 岛国毛片在线播放| 丰满少妇做爰视频| 免费女性裸体啪啪无遮挡网站| 亚洲精品国产av成人精品| 天堂中文最新版在线下载| 免费在线观看视频国产中文字幕亚洲 | 又大又爽又粗| 欧美亚洲日本最大视频资源| 久久精品成人免费网站| 国产欧美日韩一区二区三区在线| 国产一区有黄有色的免费视频| 中文字幕人妻熟女乱码| 亚洲欧美精品综合一区二区三区| 一级毛片黄色毛片免费观看视频| 在线观看国产h片| 国产亚洲av片在线观看秒播厂| 国产成人影院久久av| 在线观看免费高清a一片| 亚洲国产精品999| 午夜激情久久久久久久| 乱人伦中国视频| 十八禁人妻一区二区| 啦啦啦视频在线资源免费观看| 赤兔流量卡办理| 亚洲成国产人片在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品 国内视频| 精品久久久精品久久久| 国产在线视频一区二区| 欧美少妇被猛烈插入视频| 天天躁狠狠躁夜夜躁狠狠躁| 欧美人与性动交α欧美精品济南到| 亚洲欧美中文字幕日韩二区| 成人亚洲精品一区在线观看| 亚洲国产精品一区三区| 在线观看免费视频网站a站| 桃花免费在线播放| 纯流量卡能插随身wifi吗| 成人手机av| 精品少妇一区二区三区视频日本电影| 国产精品国产三级专区第一集| 婷婷成人精品国产| 日韩一卡2卡3卡4卡2021年| 丁香六月天网| 国产精品免费大片| 男女床上黄色一级片免费看| 午夜日韩欧美国产| 国产福利在线免费观看视频| 黄网站色视频无遮挡免费观看| 如日韩欧美国产精品一区二区三区| 成年av动漫网址| 一本—道久久a久久精品蜜桃钙片| 亚洲第一av免费看| 久久久久精品国产欧美久久久 | 美女中出高潮动态图| 成人影院久久| 五月开心婷婷网| 一本大道久久a久久精品| 美女大奶头黄色视频| 久久久精品区二区三区| 日韩制服骚丝袜av| 欧美在线黄色| 青春草亚洲视频在线观看| 国产一区二区激情短视频 | 免费黄频网站在线观看国产| av一本久久久久| 黑人欧美特级aaaaaa片| 国产在线观看jvid| 欧美大码av| 大香蕉久久网| 国产精品熟女久久久久浪| 老鸭窝网址在线观看| 蜜桃国产av成人99| 99热网站在线观看| 80岁老熟妇乱子伦牲交| 日日夜夜操网爽| 香蕉丝袜av| 在现免费观看毛片| 久久九九热精品免费| 亚洲av在线观看美女高潮| 大码成人一级视频| 亚洲第一青青草原| av国产久精品久网站免费入址| 永久免费av网站大全| 日韩一区二区三区影片| 久久精品久久久久久噜噜老黄| 亚洲国产精品国产精品| 国产亚洲欧美精品永久| 久久精品国产亚洲av涩爱| 天天躁日日躁夜夜躁夜夜| 男女下面插进去视频免费观看| 国产亚洲一区二区精品| av片东京热男人的天堂| 国产欧美日韩一区二区三 | 曰老女人黄片| 日本午夜av视频| av在线老鸭窝| 国产人伦9x9x在线观看| 90打野战视频偷拍视频| 国产精品99久久99久久久不卡| 两性夫妻黄色片| av在线app专区| 夫妻午夜视频| 国产一区二区三区综合在线观看| 国产精品久久久久成人av| 丝袜人妻中文字幕| av天堂在线播放| 丰满人妻熟妇乱又伦精品不卡| 国产在线视频一区二区| 国产精品熟女久久久久浪| 国产精品成人在线| 丝袜喷水一区| 制服人妻中文乱码| 国产av精品麻豆| 侵犯人妻中文字幕一二三四区| 日本一区二区免费在线视频| 一区福利在线观看| 国产1区2区3区精品| 男女之事视频高清在线观看 | 久久久久久久久久久久大奶| 中文字幕制服av| 777米奇影视久久| 国产伦理片在线播放av一区| 亚洲精品国产av成人精品| 十八禁人妻一区二区| av天堂在线播放| 久久久久精品国产欧美久久久 | a级毛片在线看网站| 在线亚洲精品国产二区图片欧美| 久久久国产一区二区| 捣出白浆h1v1| 成年人午夜在线观看视频| 观看av在线不卡| 美女午夜性视频免费| 午夜激情久久久久久久| 亚洲成av片中文字幕在线观看| 美国免费a级毛片| 超色免费av| 精品亚洲成国产av| 人人澡人人妻人| 99国产综合亚洲精品| 亚洲欧美激情在线| 色婷婷av一区二区三区视频| 宅男免费午夜| 咕卡用的链子| 在线观看免费日韩欧美大片| 亚洲国产毛片av蜜桃av| 久久久久久久久免费视频了| 精品免费久久久久久久清纯 | av片东京热男人的天堂| 欧美激情 高清一区二区三区| 美国免费a级毛片| 久久性视频一级片| 日本欧美视频一区| 啦啦啦在线免费观看视频4| 精品福利永久在线观看| 亚洲精品国产av蜜桃| √禁漫天堂资源中文www| 色94色欧美一区二区| 一级黄片播放器| 最新在线观看一区二区三区 | 操美女的视频在线观看| 91麻豆精品激情在线观看国产 | 晚上一个人看的免费电影| 美女大奶头黄色视频| 曰老女人黄片| 日韩一本色道免费dvd| 亚洲欧美中文字幕日韩二区| 老鸭窝网址在线观看| 丰满少妇做爰视频| 大陆偷拍与自拍| 欧美国产精品va在线观看不卡| 国产免费福利视频在线观看| 又紧又爽又黄一区二区| 久久九九热精品免费| 在线亚洲精品国产二区图片欧美| 在线观看免费午夜福利视频| 丝瓜视频免费看黄片| av片东京热男人的天堂| 美女午夜性视频免费| 国产视频首页在线观看| 亚洲精品美女久久久久99蜜臀 | 高清视频免费观看一区二区| 丰满迷人的少妇在线观看| 亚洲久久久国产精品| 久久精品久久精品一区二区三区| 国产成人精品在线电影| 99国产精品一区二区三区| 1024视频免费在线观看| 老司机靠b影院| 99热网站在线观看| 亚洲国产看品久久| 国产日韩一区二区三区精品不卡| 国产三级黄色录像| 女人爽到高潮嗷嗷叫在线视频| 欧美在线一区亚洲| 建设人人有责人人尽责人人享有的| 99精品久久久久人妻精品| 91成人精品电影| 国产亚洲午夜精品一区二区久久| 一本一本久久a久久精品综合妖精| av国产精品久久久久影院| 午夜福利免费观看在线| 亚洲国产精品成人久久小说| 日韩制服骚丝袜av| 久久精品aⅴ一区二区三区四区| 欧美日韩国产mv在线观看视频| 亚洲五月婷婷丁香| 亚洲第一av免费看| 国产成人精品久久二区二区免费| 又大又爽又粗| 国产免费又黄又爽又色| 丁香六月欧美| 中国美女看黄片| 在线观看一区二区三区激情| 国产99久久九九免费精品| 麻豆av在线久日| 在线观看人妻少妇| 中文字幕av电影在线播放| 亚洲人成电影观看| 国产在线观看jvid| 精品少妇内射三级| 国产精品一国产av| 赤兔流量卡办理| 欧美性长视频在线观看| 国产麻豆69|