基于剛?cè)狁詈夏Ｐ偷谋密嚤奂軇恿W(xué)仿真與優(yōu)化

王 瑋，李 紅，駱 群，周明華，李新松

（揚州大學(xué)機械工程學(xué)院，江蘇 揚州225127）

基于剛?cè)狁詈夏Ｐ偷谋密嚤奂軇恿W(xué)仿真與優(yōu)化

王 瑋，李 紅*，駱 群，周明華，李新松

（揚州大學(xué)機械工程學(xué)院，江蘇 揚州225127）

為使混凝土泵車臂架系統(tǒng)的計算機仿真模型更好地體現(xiàn)實際情況，根據(jù)機器人運動學(xué)理論和拉格朗日方程，通過遞推形式確定了混凝土泵車柔性臂架系統(tǒng)的動力學(xué)方程，并運用計算機聯(lián)合仿真技術(shù)對混凝土泵車臂架系統(tǒng)進(jìn)行柔性仿真，同時優(yōu)化了臂架系統(tǒng)上的驅(qū)動力，結(jié)果顯示不斷改變鉸點位置能明顯降低變幅機構(gòu)中油缸的最大受力數(shù)值.仿真結(jié)果表明，柔性臂架模型更符合實際情況，為混凝土泵車臂架改進(jìn)設(shè)計和機器人化提供了理論依據(jù).

混凝土泵車；臂架系統(tǒng)；動力學(xué)仿真；剛?cè)狁詈希挥透昨?qū)動力

混凝土泵車是一種通過液壓系統(tǒng)帶動多節(jié)可折疊臂架旋轉(zhuǎn)，同時臂架帶動輸料管運動且可連續(xù)、均勻泵送混凝土到指定位置進(jìn)行澆注的專用車輛［1］.使用混凝土泵車能有效縮短建筑工期，減少勞動消耗.混凝土泵車臂架系統(tǒng)采用多冗余度懸臂梁結(jié)構(gòu)，在其工作過程中，臂架不可避免地產(chǎn)生振動.通過建立臂架系統(tǒng)合理模型，并研究臂架末端的運動軌跡或通過軟件仿真臂架的運動過程［2-3］，可獲得降低臂架振動的方法.本文通過研究臂架動力學(xué)方程來揭示臂架的各項動力學(xué)特性，分析模型運動學(xué)和動力學(xué)，優(yōu)化了末節(jié)油缸驅(qū)動力，為今后改進(jìn)臂架系統(tǒng)設(shè)計提供參考.

1　混凝土泵車柔性臂架的動力學(xué)方程

將機器人運動學(xué)理論應(yīng)用于泵車臂架系統(tǒng)中，把液壓缸視為滑塊，帶動各節(jié)臂架實現(xiàn)回轉(zhuǎn).回轉(zhuǎn)臺坐標(biāo)系逆向動力學(xué)分析的解具有唯一性，所以在建立整節(jié)臂架系統(tǒng)方程時，只須計算臂架系統(tǒng)各桿件彈性變形的運動學(xué)方程及各桿件之間的夾角，而無須考慮回轉(zhuǎn)臺運動學(xué)的影響，故臂架系統(tǒng)的運動可轉(zhuǎn)化為平面運動.運用第二類拉格朗日方程及相關(guān)理論，通過建立遞推式可得到該系統(tǒng)的動力學(xué)方程.首先研究2節(jié)臂架系統(tǒng)［4］，建立臂架系統(tǒng)動力學(xué)方程，然后將此方程推廣到多節(jié)臂架系統(tǒng)，為泵車臂架系統(tǒng)建模提供理論依據(jù).

1.1臂架上任一點的位置、速度及加速度

在臂架系統(tǒng)所在平面上建立一個全局坐標(biāo)系x Oy（見圖1），采用隨柔性體形變而變化的浮動坐標(biāo)系.若一柔性體從某一個地點運動到另一地點，則其整體運動可分解為剛性移動、剛性轉(zhuǎn)動和變形運動.像剛性體臂架一樣，假設(shè)P為第2節(jié)柔性臂上任意一點，其位置向量為rP=r0+A（sP+uP），其中r0為浮動坐標(biāo)系的原點在慣性坐標(biāo)系中的向量，rP為P點在慣性坐標(biāo)系中的向量，sP為未變形時柔性第2節(jié)臂架上P點在浮動坐標(biāo)系中的向量，uP為柔性體的相對變形向量，A為旋轉(zhuǎn)變換矩陣.

圖1　兩節(jié)臂架模型示意圖Fig.1 Schematic of two booms

采用模態(tài)坐標(biāo)對點P的單元變形可描述為uP=ΦPq，式中q為變形的廣義坐標(biāo)，ΦP為點P滿足里茲基向量要求的假設(shè)變形模態(tài)矩陣.依據(jù)位置向量對時間一階和二階導(dǎo)數(shù)求得柔性臂架上P點的速度和加速度為

1.2臂架系統(tǒng)的能量

考慮到第2節(jié)臂架末端具有集中質(zhì)量m2，取其等效密度ρ2′=ρ2+m2δ（x2-l2），ρ2為均質(zhì)臂架單位長度的質(zhì)量，m2δ（x2-l2）表示臂架末端質(zhì)量的分布密度，對體積V積分可得第2節(jié)柔性臂架動能，式中mP和IP分別是點P的節(jié)點質(zhì)量和節(jié)點慣性張量；ωP為點P相對于全局坐標(biāo)基的角速度在局部坐標(biāo)基中的斜方陣表示；勢能為W=Wg（ξ）+2-1ξTHξ，其中H 為彈性勢能中結(jié)構(gòu)件的廣義剛度矩陣，ξ為廣義坐標(biāo)，因而，彈性應(yīng)變能與臂架整體剛性運動無關(guān)，通常為常量.

1.3臂架系統(tǒng)的動力學(xué)方程

式中K，M，D分別為模態(tài)剛度矩陣、質(zhì)量矩陣和模態(tài)阻尼矩陣；，Kξ分別為柔性體內(nèi)因阻尼和彈性變形引起的廣義力；fg為廣義重力；Q為廣義外力.

2　混凝土泵車剛?cè)狁詈夏Ｐ徒⑴c驗證

2.1混凝土泵車剛?cè)狁詈夏Ｐ徒?/p>

利用ANSYS軟件生成臂架系統(tǒng)模態(tài)中性文件，并在ADAMS泵車臂架剛性模型中依次讀入各節(jié)臂架的MNF文件.為了便于分析，在選定要替換的剛性體同時載入對應(yīng)的柔性體，ADAMS軟件便會對柔性體自動定位.替換成功后，原剛性體的運動形式、負(fù)荷及驅(qū)動力等都轉(zhuǎn)換到柔性體上，且柔性體會保持原剛性體的名稱和初始速度等特征［5］，最終得到整個臂架系統(tǒng)的剛?cè)狁詈夏Ｐ?，見圖2.

圖2　臂架系統(tǒng)剛?cè)狁詈夏Ｐ虵ig.2 Rigid and flexible model of boom system

2.2剛?cè)狁詈夏Ｐ偷尿炞C

剛?cè)狁詈夏Ｐ偷慕⑦^程中存在ADAMS與ANSYS軟件間數(shù)據(jù)轉(zhuǎn)換，進(jìn)行動力學(xué)仿真之前，必須保證信息傳遞的準(zhǔn)確性，否則仿真后處理結(jié)果毫無意義.模態(tài)中性文件的驗證可從以下3個方面分析：①使用ADAMS中的Info功能檢查柔性體的物理特性，包括質(zhì)量、質(zhì)心位置以及轉(zhuǎn)動慣量；②比較模型的模態(tài)信息，利用ADAMS/Linear模塊與有限元模態(tài)分析結(jié)果進(jìn)行比較；③對柔性體進(jìn)行仿真，并測量節(jié)點的靜力變形量或應(yīng)力值，比較其最大值與相同約束條件下有限元分析的結(jié)果［6］.

以第1節(jié)臂架為例進(jìn)行模態(tài)中性文件的驗證.打開第1節(jié)臂架的柔性體信息，與相應(yīng)剛性體比較發(fā)現(xiàn)，質(zhì)心位置重合，質(zhì)量和轉(zhuǎn)動慣量相差很小.采用ADAMS/Linear模塊，在不考慮阻尼狀態(tài)下，比較第1節(jié)臂架柔性體的固有頻率與ANSYS中的模態(tài)分析結(jié)果（見表1）.通過驗證發(fā)現(xiàn)，在相同約束條件下，臂架模型的ANSYS模態(tài)分析結(jié)果與ADAMS線性分析結(jié)果近似相等，說明模態(tài)中性文件是正確的，可利用臂架系統(tǒng)剛?cè)狁詈夏Ｐ妥鬟M(jìn)一步的仿真研究.

表1　第1節(jié)臂架固有頻率Tab.1 Natural frequency of boom 1

3　混凝土泵車剛?cè)狁詈夏Ｐ蛢?yōu)化

對驗證后的臂架系統(tǒng)剛?cè)狁詈夏Ｐ瓦M(jìn)行運動學(xué)和動力學(xué)仿真分析［7］.圖3，4為模型第5節(jié)臂架位移仿真結(jié)果.由圖可見，大范圍剛體運動和局部彈性變形引起的耦合效應(yīng)可使臂架系統(tǒng)出現(xiàn)一定程度的彈性變形，且臂架末端受到疊加影響，振動現(xiàn)象尤為明顯.這種現(xiàn)象對于輕質(zhì)長臂桿類泵車臂架系統(tǒng)來說是無法避免的，只能采取一些減振措施盡量減小振動帶來的負(fù)面影響.

圖3　第5節(jié)臂架末端點縱向位移曲線Fig.3 Longitudinal displacement curve of boom 5

圖4　第5節(jié)臂架末端點橫向位移曲線Fig.4 Transverse displacement curve of boom 5

對臂架系統(tǒng)剛?cè)狁詈夏Ｐ头抡婧蟀l(fā)現(xiàn)，臂架末端振動最明顯，油缸受力也比較復(fù)雜.選取油缸5作為研究對象，對油缸受力進(jìn)行優(yōu)化.首先，對臂架油缸5處的鉸接點、連桿連接點等位置進(jìn)行參數(shù)化處理，如圖5所示.

圖5　臂架4與臂架5間各鉸點Fig.5 Parametric name of hinge points between boom 4 and 5

創(chuàng)建設(shè)計變量v1，v2，…，v10代替鉸鏈位置的坐標(biāo)，各點坐標(biāo)為A（v1，v2），B（v3，v4），C（v5，v6），D（v7，v8），E（v9，v10）.理論上，應(yīng)把10個參數(shù)全部作為設(shè)計變量進(jìn)行優(yōu)化設(shè)計，但在此空間中尋求最優(yōu)解計算量巨大，效率極低；因此，在優(yōu)化迭代前，先選用靈敏度對設(shè)計變量進(jìn)行分析，過濾靈敏度值低的變量［8-9］.結(jié)果顯示，v3，v6，v7，v8在初始位置時的靈敏度值較其他變量高，意味著它們對臂架系統(tǒng)變幅機構(gòu)影響較大；因此，選用v3，v6，v7，v8對油缸5受力影響較大的變量進(jìn)行優(yōu)化設(shè)計，從而找到使油缸5最大受力減小時各鉸點的坐標(biāo)值.

圖6為第5節(jié)臂架油缸受力的優(yōu)化曲線，可見曲線在橫軸上下不規(guī)則地跳動，說明油缸受力非常復(fù)雜.在第5節(jié)臂架回收過程中，由于輕質(zhì)長桿存在振動，必然導(dǎo)致油缸受力不穩(wěn)定，在臂架回收接近終點時圖像振蕩幅度減弱，說明油缸受力已趨于穩(wěn)定.這是由于第5節(jié)臂架在接近收攏時，油缸5與第5節(jié)臂架處于坐標(biāo)系原點同側(cè)，力臂減小，因而第5節(jié)臂架振動幅度減弱，油缸受力趨于平緩.優(yōu)化后的油缸受力基本小于優(yōu)化前的受力，而且通過多次改變變幅機構(gòu)中鉸點坐標(biāo)值，使得油缸5在最大受力工況下的最大受力值逐漸降低［10］，最終由438 kN減小到286 kN，減幅達(dá)34.6%.通過多次改變變幅機構(gòu)中鉸點的坐標(biāo)值，可達(dá)到減小油缸最大受力的效果，且為了避免因改變鉸點位置而使系統(tǒng)內(nèi)發(fā)生干涉現(xiàn)象，鉸點坐標(biāo)的優(yōu)化幅度為0.3%～2.3%（見表2）.

表2　優(yōu)化前后設(shè)計變量的比較Tab.2 Comparison of design variables during optimization

圖6　油缸5受力優(yōu)化曲線Fig.6 Optimization curve of oil cylinder 5

［1］LIN Peichun，WANG Jenhun，HUANG Shanhun，et al. Dispatching ready mixed concrete trucks under demand postponement and weight limit regulation［J］.Automat Constr，2010，19（6）：798-807.

［2］陳棟，王剛，謝秀芬.混凝土泵車臂架振動的動態(tài)特性 ［J］.機械設(shè)計與研究，2011，27（4）：92-95.

［3］戴麗.基于柔性多體動力學(xué)的混凝土泵車臂架系統(tǒng)的建模與仿真 ［D］.沈陽：東北大學(xué)，2007.

［4］WONG K V，BARIN E.Oil spill containment by a flexible boom system［J］.Spill Sci Technol B，2003，8（5）：509-520.

［5］馮波.混凝土泵車臂架系統(tǒng)的仿真分析 ［D］.長沙：長沙理工大學(xué)，2012.

［6］朱大煒.基于模態(tài)綜合法的雷達(dá)天線機構(gòu)的柔體動力學(xué)仿真 ［D］.上海：上海交通大學(xué)，2007.

［7］CAZZULANI G，GHIELMETTI C，GIBERTI H，et al.A test rig and numerical model for investigating truck mounted concrete pumps［J］.Automat Constr，2011，20（8）：1133-1142.

［8］YAN Shangyao，LAI Weishen，CHEN Maonan.Production scheduling and truck dispatching of ready mixed concrete［J］.Transport Res E，2008，44（1）：164-179.

［9］姚菁琳，鄭翔，張新國.八連桿壓力機傳動系統(tǒng)的設(shè)計及優(yōu)化 ［J］.揚州大學(xué)學(xué)報：自然科學(xué)版，2014，17（1）：46-49.

［10］SONG Ping，LI Qingzhou，LI Keiie.Sensitivity optimization of a monolithic high-shock three-axis piezoresistive accelerometer with single sensing element［J］.Chin Sci Bull，2009，54（19）：3600-3607.

Dynamics simulation and optimization of the boom system of the concrete pump truck based on rigid and flexible model

WANG Wei，LI Hong*，LUO Qun，ZHOU Minghua，LI Xinsong
（Sch of Mech Engin，Yangzhou Univ，Yangzhou 225127，China）

Combining with the theories of robotics and Lagrange’s equation，a dynamical function for flexible boom of concrete pump truck is established by the recurrence methods in order to make the computer model better reflect the actual situation.Based on the dynamical function，the rigid and flexible model of concrete pump truck is built and simulated through the ioint application of ANSYSand ADAMS.Using the rigid and flexible mode of boom system，the cylinder stress of luffing mechanism are studied by using optimization module after verifying the accuracy of dynamic model.It is found from the simulation results that the rigid and flexible model provides a theoretical basis for truck improvement and robotizing.

concrete pump truck；boom system；dynamics simulation；rigid and flexible；cylinder stress

TU 646

A

1007-824X（2015）01-0061-04

（責(zé)任編輯 賈慧鳴，秋 實）

2013-10-16.*聯(lián)系人，E-mail：lhchina@163.com.

國家科技部火炬計劃資助項目（2012GH530716）.

王瑋，李紅，駱群，等.基于剛?cè)狁詈夏Ｐ偷谋密嚤奂軇恿W(xué)仿真與優(yōu)化［ J］.揚州大學(xué)學(xué)報：自然科學(xué)版，2015，18（1）：61-64.