尹建軍　高　強　陳亞明

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

?

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

尹建軍高強陳亞明

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

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

打捆機; 打結(jié)器; 接觸動力學; 虛擬打結(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)方法來驗證設計優(yōu)劣，未經(jīng)虛擬打結(jié)驗證設計優(yōu)劣的打結(jié)器試制，可能導致試制失敗而使制造投入面臨風險。因此，有必要研究打結(jié)器虛擬打結(jié)實現(xiàn)方法來完善打結(jié)器設計理論。

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

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

1　捆繩動力學模型

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

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

根據(jù)力學理論，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)角、速度、角速度有關。在此，捆繩選直徑D=2.5 mm的聚丙烯繩，長度L=0.48 m。根據(jù)聚丙烯的力學特性[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ù)，使整段捆繩模型的力學性能接近捆繩真實力學性能，從而可模擬整段捆繩的變形及動力學性能。

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

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格式，再將該中間文件導入到ADAMS軟件中，對導入的模型修改質(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.夾繩壓板

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

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

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

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

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

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

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

設置該力的意義是：在割繩刀接觸到打結(jié)嘴與夾繩盤之間的捆繩時，夾繩盤處捆繩末端拉力消失，可近似模擬捆繩被割繩刀片割斷。通過多次仿真嘗試設定時間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ū)動及約束設置和載荷施加后，利用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ù)設置正確，可以模擬捆繩的大變形及動力學行為，能再現(xiàn)捆繩在拉力作用下與剛性構(gòu)件之間的接觸約束過程。利用虛擬打結(jié)可實時觀察捆繩在拉力作用下形成繩結(jié)的剛?cè)峤换ミ^程，包括撥繩盤將捆繩撥向打結(jié)嘴搭接位置的動態(tài)過程、捆繩被刀臂的彎曲弧面導入弧形槽口、捆繩被夾繩器夾持旋轉(zhuǎn)、捆繩被打結(jié)嘴鉸扭成環(huán)、鉤鉗張開鉗咬捆繩、鉤鉗閉合鉤出繩頭和刀臂強制將繩扣從打結(jié)嘴上推脫。對于捆繩模型，文獻[4]僅考慮了捆繩剛度系數(shù)，且剛度系數(shù)矩陣為定值，未考慮捆繩阻尼系數(shù)矩陣，而本文建立的捆繩動力學模型接近捆繩真實力學性能，可獲得較為準確的整段捆繩的變形及動力學性能。

在虛擬打結(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é)的試驗驗證

將設計的雙齒盤驅(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)方法可以用于設計階段檢驗打結(jié)器的結(jié)構(gòu)改進和參數(shù)調(diào)整后的成結(jié)效果，保證打結(jié)器物理樣機一次性制造的成功，大大降低了打結(jié)器制造的投入成本。

3.3虛擬打結(jié)應用

打結(jié)嘴和夾繩盤動作時序差φ和打結(jié)嘴軸線與主軸軸線夾角α是打結(jié)器設計的兩個基本參數(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è)方，導致鉤鉗無法鉗咬住捆繩，鉗咬動作失敗而未能成結(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é)器設計方案。利用該方法能解決打結(jié)器設計后難以直觀驗證成結(jié)效果的問題，在設計階段可發(fā)現(xiàn)設計缺陷，降低打結(jié)器制造的投入成本。

4　結(jié)論

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

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

1尹建軍, 張萬慶, 陳亞明, 等. 打結(jié)器夾繩-繞扣-鉗咬動作參數(shù)分析與打結(jié)試驗[J]. 農(nóng)業(yè)機械學報, 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é)器關鍵傳動機構(gòu)設計[J]. 農(nóng)業(yè)機械學報, 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è)機械學報, 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è)機械學報, 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 帶單元在繩索體仿真中的應用研究[J]. 機械科學與技術(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)仿真學報, 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)仿真學報, 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纖維繩索的編織工藝和力學性能研究[D]. 青島: 青島大學, 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é)器設計與成結(jié)試驗分析[J]. 農(nóng)業(yè)機械學報,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石明全, 張鵬, 范樹遷,等. 基于多接觸約束的動力學仿真計算效率[J]. 機械工程學報, 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è)機械學報, 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

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

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

S491

A

1000-1298(2016)09-0085-08