2016, Vol. 37
当今基础设施建设中对于施工设备的新要求成为了工程机械整机厂商的市场竞争关键，工程机械产品需要更高的挖掘速率、更低的能耗等.系统仿真技术能够在设计阶段对整机动态性能进行预测，这为工程机械开发提供了新的工具[1-2].其中平面机构建模是进行系统动态仿真研究的重要环节.目前最有效的系统动态建模方法是功率键合图法.
键合图法表达机械结构时，由于是多刚体的多维问题，需要进行全局局部坐标转换，这使得键合图模型变得复杂，因为刚体间存在铰接，通常会产生微积分代数环[3]问题.商用仿真软件如AMESim中对于机构的求解是用虚拟弹性铰接点[4]或是拉格朗日乘子[5-7]来保持铰接点的因果关系链接，这样引入新的动态环节来打破原有的因果链以防止微积分代数环的出现.这会导致整机系统状态方程出现严重的刚性问题，在连接液压、电控整机系统之后的大系统仿真运行变得十分困难.要想根本性解决该问题，需要重新确定因果关系[8].近年来对于键合图方法的研究主要集中在针对具体工程对象的建模过程[9]，而缺乏对于因果关系确定的研究.文献[10]是采用运动学解析方法针对装载机工作装置的建模实例，虽然实现了模型的高效仿真，但模型各个接口点势流关系混乱，无法规范化和模块化.
可见寻求一种统一的建模思路十分必要.本文提出了一种针对该机构的键合图因果关系分类整理方法，通过元件接口势流因果分类及部件间连接规则整理两个方面，使得刚体间连接键的因果关系得以延续，并且建立了元件间连接的规范，该分类整理可使得建模过程高效不易出错.将能够适用于任何形态结构的工作装置建模.
1 直线驱动铰链机构的定义以典型反铲挖掘机为例，其工作装置如图 1所示.
|
图 1 典型反铲挖掘机工作装置机构 Fig.1 Typical mechanism of backhoe excavator 1—动臂；2—斗杆；3—铲斗；4—动臂油缸；5—斗杆油缸；6—铲斗油缸；7—摇杆；8—摇杆；A~K—机构各铰接点；Q—铲斗尖点. |
可以看到，工程机械机构主要由两类元件组成：
1) 转动副连接的刚体部件.如动臂、摇杆、铲斗等.
2) 直线驱动器.其特点为2个转动副连接其他部分，自身内包含一个滑动副，并作为动力输入.在工程机械中通常为液压缸，如挖掘机的动臂缸.
这里将由以上两类元件组成的机构称为“直线驱动铰链机构”.本文针对此类机构进行研究.
2 机构键合图模型的元件分类整理机构键合图模型的建模详细过程可参考文献[10].这里简单介绍建模过程，并对其进行分类.
对于动臂、斗杆等刚体元件的建模，为了规范通用表达，以重心位置为局部坐标系的原点，则所有铰接点均有模数L和角度α两个参数.例如动臂1的机构如图 2a所示，其有4个铰接点，则有4个模数L和4个初始角度α.
设在其运动过程中转过的角度为θ，根据坐标系转换，以C点为例，其速度方程为
|
(1) |
根据方程可确定键合图变换器MTF的系数，其余铰接点同理计算，则其键合图模型如图 2b所示.
|
图 2 动臂键合图模型 Fig.2 Bond graph model of boom (a)—动臂元件参数定义； (b)—键合图模型. |
对液压缸建模，以动臂油缸为例，设其长度为m，伸出量为n，转过的角度为θ，如图 3a所示.
根据坐标系转换，以A点为例，其局部坐标速度应满足：
|
(2) |
同样根据方程可确定键合图变换器MTF的系数.则其键合图模型如图 3b所示.
|
图 3 液压缸键合图模型 Fig.3 Bond graph model of cylinder (a)—液压缸参数定义； (b)—键合图模型. |
对于连杆，其铰接点仅2个，第2个与刚体部件相连接.其规则与液压缸的接口类似，但模型结构不同.因此这里省略其推导过程，键合图如图 4所示.
|
图 4 连杆键合图模型 Fig.4 Bond graph model of rod |
综上，一套工作装置的组成可分为3类：液压缸，刚体部件和连杆.这3种部件可以组合出工程机械的所有直线驱动铰链机构.
3 元件接口因果关系的分类整理 3.1 接口类型的分类势流因果关系在本问题上的直接体现，即力和速度谁为输入量，谁为输出量的问题.以下分别针对刚体部件、液压缸和连杆的接口进行分类整理.
从图 4可以看到，模型对于每一个部件来说，其他变量的速度量均由一个基准点出发求解，其他铰点的平移、转动均在该基准点基础上进行.这个点可理解为工装部件的“根”，故以下称其为根点，用Rt表示.以挖掘机动臂为例，其Rt点为C点.
刚体部件中，除了根点以外，可以看到模型的其他铰接点中有一个点为x，y分解的势流源类型不同，其余均为速度输出(势源)点.而且可以看到，势流源类型不同的点为连接液压缸(或连杆)的点.从模型求解过程看，这个点需要获取外模块的力的信息(液压缸给的作用力)，也要获取速度信息(两个速度点才能确定当前姿态).以下称其为外力顶点，用Pd表示.以挖掘机动臂为例，其Pd点为B点.
除了Rt点和Pd点外，其他铰接点均为速度输出，力输入的点.以下称其为负载点，考虑到其可能有多个点，用Li(i=1,2,…)表示.以挖掘机动臂为例，其负载点为D，F点.
可以看到，刚体部件的一般特征：一个Rt点，Rt点的x,y端均为速度输入(流源)；一个Pd点，Pd点的x端为力输入，y端为速度输入(x为势源，y为流源)；若干个Li点，x,y端均为力输入.
对于液压缸来说，其同样具有Rt点，例如动臂油缸的Rt点为A点.而油缸是作为力源输出给工装部件的，因此其另一个点为与工装部件的外力顶点连接的点，由于其x，y的势流关系与刚体部件相反，为表示区别，以下称其为力输出点，用Pr表示.同时，油缸能够伸缩，其伸缩的接口为流源接口，与工作机构以外的系统连接(如液压系统)，以下用I表示.
液压缸模型的一般特征:一个Rt点，点的x，y端均为速度输入(流源)；一个Pr点,x端为速度输入，y端为力输入(x为流源，y为势源)；一个I接口，为速度输入(流源).
对于连杆，可看作不能伸缩的液压缸，因此定义与液压缸相同.
势流因果关系统一定义后的元件模型如图 5所示.
3.2 部件间的连接规则定义当各个部件如此分类后，根据键合图模型的因果关系传导，部件需符合以下规则：一个刚体部件可以有多个Li点，但必须有且只能有一个Rt点，一个Pd点；液压缸有且只有一个Rt点，一个Pr点，一个I接口；连杆有且只有一个Rt点，一个Pr点.
部件间连接可按以下规则进行：连接支架(车体)的点均为Rt点；Li点可以空闲(不连接其他点)；未连接支架的Rt点需连接于Li点上，Pd点必须连接Pr点上.
|
图 5 工作装置元件键合图因果关系统一分类 Fig.5 Classifications of element bond graph models (a)—刚体部件； (b)—液压缸； (c)—连杆. |
例如：
1) 如刚体部件2连接在刚体部件1上，则铰接点对于1为Rt点，对于2为Lx点；
2) 如刚体部件1和刚体部件2通过液压缸或连杆3连接并动作，则3连接2的铰接点对于2为Pd点，对于3为Pr点；3连接1的铰接点对于3为Lx点；对于1为Rt点.如图 6所示.
|
图 6 键合图模型因果关系分类的连接示例 Fig.6 Demonstration of classified element connections |
按以上规则进行连接，即可实现工程机械中几乎所有形态的工作机构，同时可有效避免模型错误.例如挖掘机模型可表示为图 7所示.
|
图 7 挖掘机工作装置键合图模型连接规则示例 Fig.7 Example of excavator bond graph element model connections |
本建模方法最大优势在于建模过程带来的便利性.为证明建模方法正确有效，这里用本文提出的方法对一种反铲挖掘机进行建模，同时采用AMESim软件进行建模，并用两种仿真模型进行仿真对比分析.
应用本方法对反铲挖掘机建模通过MATLAB/Simulink进行，其模型中动臂的主要参数见表 1.
|
表 1 反铲挖掘机动臂主要参数 Table 1 Main parameters of backhoe excavator boom |
以下是一个典型工况的仿真示例，是对一种工况的仿真结果，其仿真时长为10s，工作装置做挖掘—举升—卸料的动作.其初始和结束时的工作装置机构姿态用OpenGL绘图，如图 8所示.
|
图 8 挖掘机工作装置仿真工况绘图 Fig.8 Excavator geometry graph in simulation status (a)—初始姿态； (b)—结束姿态. |
分别采用本文方法的模型以及AMESim软件进行仿真,得出的挖掘机铲斗尖点Q的运动轨迹如图 9所示.
|
图 9 两种模型得出的仿真关键点空间的运动轨迹 Fig.9 Trajectories of key joints in excavator mechanism with different model |
可以看到两个仿真模型得出的结果十分吻合，这证明了本文建模方法的有效性.由于AMESim采用虚拟弹性铰接点方法求解刚体间连接关系，其铰接点被设定为通过高刚度弹簧阻尼连接，因此存在微量的位移差异来产生接触力.
5 结论本文提出了一种针对直线驱动铰链机构的因果关系分类整理，形成了标准化的建模思路.机构可分为液压缸、主部件和连杆三种元件，并分别分类了统一势流因果关系接口.按该规则可定义部件间连接规则，使得工作机构的建模中势流关系简单明了，适用于直线驱动铰链机构的系统建模.
| [1] |
Paolo C, Alvin A. Gray box modeling of an excavator′s variable displacement hydraulic pump for fast simulation of excavation cycles[J].
Control Engineering Practice, 2013, 21 (1) : 483 –494.
( 0)
|
| [2] |
Manfred H. Modelling simulation and control design for large and heavy manipulators[J].
Robotics and Autonomous Systems, 1996, 19 (1) : 167 –177.
(0)
|
| [3] |
Cacho R, Felez J, Vera C. Deriving simulation models from bond graphs with algebraic loops:the extension to multibond graph systems[J].
Journal of the Franklin Institute, 2000, 337 (1) : 579 –600.
(0)
|
| [4] |
Karnopp D C, Margolis D L. Analysis and simulation of planar mechanism systems using bond-graph[J].
Journal of Mechanical Design, 1993, 101 (2) : 187 –191.
(0)
|
| [5] |
Gawthrop P J, Smith L S. Causal augmentation of bond graphs with algebraic loops[J].
Journal of Franklin Institute, 1992, 329 (2) : 291 –303.
(0)
|
| [6] |
Wolfgang B,Cellier F.Tearing in bond graphs with dependent storage elements[C]//Proceedings of the CESA′96 Conference.Lille,1996:1113-1119. (0)
|
| [7] |
Felez J. A bond graph based simulation program[J].
Transaction ASME Journal Dynamic System Measurement Control, 1990, 112 (1) : 717 –727.
(0)
|
| [8] |
Saeed B, Amir K. Causality in vector bond graphs and its application to modeling of multi-body dynamic systems[J].
Simulation Modelling Practice and Theory, 2006, 14 (1) : 279 –295.
(0)
|
| [9] |
Wolfgang B. Bond graph modelling and simulation of multidisciplinary systems:an introduction[J].
Simulation Modelling Practice and Theory, 2009, 17 (1) : 3 –21.
(0)
|
| [10] |
Romero G, Felez J, Maroto J. Efficient simulation of mechanism kinematics using bond graphs[J].
Simulation Modelling Practice and Theory, 2009, 17 (1) : 293 –308.
(0)
|