平面2自由度并联机构动力学性能优化

本资料包含pdf文件1个，下载需要1积分

Dynamic Performance Optimi zation for Planar 2-DOF Parallel M echanismGA0 Xiulan。LU Kaijiang，SHI Gemeng 。

(Baoi University of Arts and Sciences，Bali Shaanxi 72 1007，China)Abstract：The transitive relation of velocity an d acceleration of planar 2-DOF paralel mechanism between input and output wasobtained by using influence coeficients.Driving torque expression of active joints was acquired through establishing dynamic equationswhich was based on selecting mechanismS output as generalized coordinates.Then the maximum of the joint driving torque in the aldesign space was taken as the mechanisms dynamic performance evaluation index.The sensitivity of the performance to structural pa·ram eters was analyzed. In order to get the optimal structural param eters which satisfied both design space and dynamic perform ance，while to get the optimal dynamics perform ance in the scheduled design space，the constraints of structural paramete were establishedaccording to the requirement for the design space of the mechanism。

Keywords：Parallel mechanism；Infl uence coeficient；Dynamic perform ance；Optimal design近年来，并联机器人的应用领域在不断扩大，随着科学技术的发展，对机构的性能提出了越来越高的要求。工作空间是机器人最重要的性能指标，合理地定义工作空间是机器人结构设计的首要环节。最大限度地提高机构的运动速度，是现代机器人技术发展的主要标志，为了适应高速和高加速度的要求，许多学者从运动学和动力学方面L 提出了多种动态性能指标来量化机器人机构的动态性能。然而并联机构的工作空间与动态性能之间往往是矛盾的，为了获得良好的动态性能，只有放弃使用部分工作空间。目前并联机构的结构设计主要是基 于运动学 和动力学 J性能指标进行的，而兼顾机器人工作空间 的设计较少，其主要原因是并联机构工作空间的定义与描述本身就是-个难题。机构的工作空间与动态性能在很大程度上撒于机构的结构参数，对机构进行优化设计将是-种可行的方法，它可以兼顾工作空间和动态性能这两个方面，并使它们之间达到均衡，保证机构的总体性能最优。作者以关节驱动力矩在预定设计空间的最大值作为评定机构动力学性能的指标，将机构对设计空间的要求看作是对结构参数的约束，从而将机构的优化设计归结为约束优化问题，得到了兼顾设计空间和动力学性能的最优结构参数。

1 平面2自由度并联机构逆运动学分析平面2自由度并联机构的末端点C通过两条支链与基座相连，如图1所示。取参考坐标系Oxy。

J，IC/ / 0 A2/ J图1 平面2自由度并联机构运动简图收稿日期：2012-01-16基金项目：陕西省自然科学研究计划项目 (2011JM7015)；陕西侍育厅专项科研计划项目 (09JK327)作者简介：高秀兰 (1965-)，女，工学学士，教授，主要从事曲面数控加工编程技术及并联机器人的研究。E-mail：wgyq92wl### 163.com。

· 66· 机床与液压 第41卷选鳃构-个分支，构成开式运动链 ，该分支关节速度到动平台的运动速度 的映射关系[ [- 三， - [耋 ]垒[G：](r) 式中： 为末端点C的位置矢量， [ Yr，上标 r表示分支序号， 中的下标表示分支中运动副的序号。

当已知动平台的速度，可求得该分支关节相对速度 [G：] [G (1)同样，可以得到每个分支关节运动加速度 到动平台的运动加速度矢量 的映射关系 [茹 梦] [ ] [G：] 西式中：[ ] R ，称为二阶影响系数矩阵，其每个元素都是-个2维矢量f，n-X yY Y - m >式中： 、)， 、 、y 的下标与该分支关节序号对应。

当已知动平台的加速度矢量 ，可求得该分支关节相对加速度 [G ] - [ ] (2)2 并联机构凯恩方程2.I 输入构件广义主动力和广义惯性力计算选取末端点 C的坐标 、Y为广义坐标。由式(1)得(，)[G [G输入杆质心G，的速度J：。， (r)[-sin咖 cos ：口 何(3)式中：独立速度 、 的系数为质心G，对各独立速度的偏速度。

， (-sin， fcosl )[Gr断 (-sinl icosl .，)[G ]：输入构件A，B，对各独立速度的偏角速度分别为∞ [G ]：Ik ∞ [G：]：作用在输入构件上的主动力有重力及驱动力矩，取质心为简化中心，其主矢和主矩分别为Rc-mI gJMc ，输入构件对各独立速度的广义主动力和广义惯性力fK印RG 唧MG。 。yi ： . 。yLc . 。y y ，y式中：足 -m≮ ，口 为质心的加速度，将式 (3)对时间t求导得nG：L1 -.1(r (-sinl fc。s咖： 0)- L1 2l(r 。

(cosl isinl ) - .f 粤r ×口 dm-a12l-L2 (kr 表示质量微元 dm相对简化中心 (质心)的矢径，口 为微元 dm相对简化中心的加速度。

口 f 1 (-sin1 icoscb1 )-f咖 ·(cosl isinl )2.2 从动构件广义主动力和广义惯性力计算由式 (1)得关节的相对角速度f )[G [G [G： [G 2，2Y从动构件角速度矢量 (咖 )k从动构件对各独立速度 、 的偏角速度依次为∞ [G [G k∞似[G [G k质心 Q的速度，。 (r)(-sin4，1 icosl ) ( r]咖 )[-sin($l 十咖 )icos(4，1 咖 ).，] (4)按独立速度整理 lQyy式中：独立速度 、 的系数为质心 Q对各独立速度的偏速度。

H ：L (-sin icol .j)[G [-sin(咖： 咖 )fcos(l ).J]([G ] [G ] ： )L (-sin，l icosl J)[G [-sin(咖： 咖 )icos(咖： 咖 n) ]([ ]： [G ] )机构在末端点 C处作用有集中力F 和F (见图2)，它使得从动杆产生附加力，由于杆的重力比F和 要小得多，因而忽略杆的重力而将从动杆看作二力杆，认为从动杆对末端点的作用力沿杆轴线方向，对末端点进行静力分析求得附加力第3期 高秀兰 等：平面2自由度并联机构动力学性能优化 ·67·式中：r和s表示分支序号，6 ： ，6 币 。

图2 末端点C的受力分析作用在从动构件上的主动力有重力和附加力。取质心为简化中心，其主矢和主矩分别为RQ sin6 -F cos8sin(6 -艿 )F sin6 -FycoS8sin(6 -艿 ) i rJ-m2g JMQ：0从动构件对各独立速度的广义主动力和广义惯性力fK RQ。H竹MQ。∞i ：露 . Qy . 竹式中：R -m：口Q，口Q为质心加速度，将式 (4)对时间t求导得口p厶咖。(-sinl icosl ) ( ： )·[-sin( ： )icos(4，1 咖 ).，]- o ·(COs4，1 isin4，1 )- /'2-( ) [cos( ： )isin(bl ).J]争 .争三 -f生， xa dm-f生z( )qdlk- )将所有构件上的广义主动力和广义惯性力分别求和，得到机构相对于各个独立速度广义主动力和广义惯性力2 2 ∑ ∑K2 2 ，Y ∑ ∑式中： 、 分别为分支 r上输入构件对独立速度 的广义主动力和广义惯性力；、 分别为分支 r上中间连杆对独立速度的广义主动力和广义惯性力。

基于凯恩方法，可以得到机构的动力学方程Ky0 ，Y写成矩阵形式7 ( ) 曰( ) P(x) (5)式中： 为主动关节驱动力矩矢量，下[r。r：]T；y( )为操作空间的惯性矩阵，B(x)∈ ；P(x)是克服外界负载和重力所需的关节驱动力矩，仅与机构的位形有关。

3 动力学性能优化3.1 机构的工作空间与设计空间机构每个分支的工作空间是以 ，为中心的环形区域，如图3所示，机构的工作空间是各支链子空间的交集，它由图中4段圆弧围成。如果要求工作空间能包容-个边长为 b的正方形，中心在 P( ，Yo)，其顶点为 C[ [式中： 为 与 轴的夹角， "IT(i-1) ，i：1-4。

图3 机构工作空间解析如果给定设计空间，则机构结构尺寸必须满足下面的几何约束条件：(1)过 A。向 C C4引垂线，如果垂足 D 位于c c 之间，则要求l z，-f，l≤A。Dl如果垂足位于C 之外，则要求I z1-f3 I≤A1Cl同样对于分支 2，如果垂足位于C C：之间，则要求· 68· 机床与液压 第 4l卷l f -f4 1≤A2 D2如果垂足位于C c 之外，则要求l z2-z4 l≤A2c2(2)在三角形 A G A 中l z -z，ll z -z l≥20(3)各支链子空间外边界应包含 c。、c：、c 、C ，即fl i>max/( c.0) Yc ，i1，2，3，4z2z4≥max/( c -口) y2c ，i：1，2，3，4(4)在△A。A：G 和△A A G 中f2z4-I z1- I≤20f z3-I z：-z4 I≤2口3.2 机构动力学性能指标和动力学优化设计对机器人进行轨迹规划或实际控制时，必须考虑在整个设计空间以及机构运行速度和加速度的范围内，关节驱动力矩的最大值丁7- max ( ， ， ，b)，i1，2， ∈D，l l≤X ， l I≤ (6)式中：b为设计变量，包括机构的结构参数和动力参数；i为分支序号；D为设计空间； 、 分别为允许运行的最大速度和加速度；r 为主动关节的驱动力矩。

r可以作为评定机构动力学品质优劣的性能指标，它代表机器人在整个工作空间的总体性能。

并联机构由多条分支组成，其工作空间是各分支子空间的交集，与传统的串联机器人相比，并联机构最大的不足就是工作空间相对较校另外在工作空间中心部位，动力学性能良好，由里向外，机构的动力特性变差。为了获得良好的动力学性能，就要放弃使用靠近边界的这部分工作空间。对机构进行动力学优化设计，就是要兼顾设计空间与动力学性能这两个矛盾的方面，将机构对设计空间的要求看作是对结构尺寸的约束，从而将动力学优化描述成：给定机构的设计空间，确定-组结构参数和动力学参数，在结构参数满足工作空间要求的几何约束条件下，使机构的某种全局性能指标在整个设计空间最优。有实用意义的优化模型有：(1)对于仅在低速运行的机构，为了使机构结构最紧凑，占用的空间最小，应使其结构尺寸的总和最小，即在结构尺寸满足第3.1节中的几何约束条件下，使f (b )minf(6) ，(b)0LlL2L3L(7)(2)对于加减速频繁和高速运行的机构，在结构尺寸满足第3.1节中的几何约束条件下，应使机构的动力性能最优，即r (b )mint (b) (8)式 (6)的性能指标的含义是 ：在给定机构的结构参数和动力学参数以后，使机构的末端点以允许范围的速度和加速度在整个设计空问运行时，各个主动关节驱动力矩的最大值。

4 应用实例根据使用要求，机构工作空问包含-个边长b200 mm、中心坐标P (50 mm，300 mm)的正方形区域，将其视为机构预定的设计空间，并且要求机构工作时，末端点在此空间能达到的最大速度和加速度 (±2.5，±2.5)m/s， (±250， ±250)m/s 。杆AB和BC为均质杆，直径分别为32 mm和25 mm。

让性能指标 分别取-系列指定值，可以得到关节驱动力矩在工作空间的等位线，如图4所示∩见：在工作空间的中部，关节驱动力矩小，机构动力学性能良好，由里向外，机构的动力学性能逐渐变差。如果将预定设计空间也表示在图中，将不难看出机构在设计空间中的最大驱动力矩。

. 0.3 -0.2 -0.1 0 0.1 0.2 0.3xlrn图4 工作空间边界及关节驱动力矩等位线图以机构结构参数为设计变量 b [ ， ， ，L ，n]，求解优化模型 (7)和 (8)，得到这两种情形下的最优结构参数，如表1所示。按尺度最优的设计方案绘制机构的工作空间，如图5所示，它与设计空间的关系表明，所得的结构参数是满足设计空间要求的最紧凑的构型要求。

第3期 高秀兰 等：平面2自由度并联机构动力学性能优化 ·69·表 1 两种情形下的最优结构参数图5 尺度最优方案下的工作空间与设计空间为了测试优化结果对动力学性能影响的程度，用动力学仿真的方法，计算末端点在整个设计空间运动时的关节驱动力矩，结果如图6所示，关节驱动力矩在整个预定的设计空间峰值小，变化平缓，变化范围校这预示着，按最优结构参数设计的机构，在预定的设计空间具有较好的动力学性能。