魔神 发表于 2009-11-22 13:49:48

并联式钢坯修磨机器人动力学及轨迹规划研究

&nbsp;&nbsp;&nbsp;&nbsp;<P align=left><FONT face=宋体><STRONG>1 引言</STRONG></FONT></P><P align=left><FONT face=宋体>  本文介绍了一种新型的三自由度并联式钢坯修磨机器人机构,它是由运动平台、固定平台、空间平动机构及连接两个平台的三个分支组成,每个分支之间有一个作为驱动的移动副,其他关节均为虎克铰(图1).它的工作原理是由三个驱动轴控制运动平台沿<EM>X、Y、Z</EM>三个方向移动,使安装在运动平台上的砂轮对钢坯进行修磨.其作用等效于一台三轴联动的数控磨床,但它的结构却比传统磨床简单许多.另外,该机构具有运动学和动力学正反计算简单、作业空间大和无位置奇异等优点,这些优点给其设计和使用带来许多方便.本文首先提出这种并联机构运动学的正逆解,然后应用拉格朗日方法讨论了该机构在操作空间和关节空间的动力学方程,再根据钢坯修磨的作业特点进行了修磨轨迹规划,尤其是对拐角轨道进行了专门设计,最后对运动学和动力学进行计算机仿真.这些工作将对钢坯修磨机器人的控制和结构优化设计提供有益的依据.</FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_oqyzpi20083181515.gif"></P><P align=center><FONT face=宋体><STRONG>图1 三自由度度联式钢坯修磨机器人机构</STRONG></FONT></P><P><STRONG><FONT face=宋体>2 运动学方程</FONT></STRONG></P><P><FONT face=宋体>  在图1所示的机构中上平台为固定平台,下平台为运动平台,它们均为正三角形.将坐标系<EM>XYZ</EM>放置在固定平台中心<EM>O<SUB>b</SUB></EM>,<EM>Z</EM>轴向下.设运动平台中心<EM>O<SUB>p</SUB></EM>在坐标系中的坐标为{<EM>x<SUB>p</SUB> y<SUB>p</SUB> z<SUB>p</SUB></EM>}<SUP>T</SUP>.因为运动平台只有三个移动自由度,而位姿并无变化,所以机构的位置逆解<SUP>[1]</SUP>为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_1otg77200831815123.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(1)</FONT><P align=left><FONT face=宋体>式中<EM>d</EM>=(<EM>R</EM>-<EM>r</EM>)/2,R和r分别是上下平台外接圆的半径.<BR>  从上式中可求出运动平台中心的坐标,即机构位置正解为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_6gfjvb200831815145.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(2)</FONT><P align=left><FONT face=宋体>由此可以求出运动平台的雅可比矩阵为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_l4ybni20083181525.gif"><FONT face=宋体>(3)</FONT><P align=left><STRONG><FONT face=宋体>3 动力学方程</FONT></STRONG></P><P align=left><FONT face=宋体>  下面应用拉格朗日方法对机构的动力学进行分析,对于任何机械系统,操作空间内的主动力为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_ww6wkg200831815227.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(4)</FONT><P><FONT face=宋体>式中:<EM>L</EM>称为拉格朗日函数,表示机构的动能K和势能<EM>P</EM>之差,<EM>q<SUB>j</SUB></EM>为机构的广义坐标,<EM>q<SUB>j</SUB></EM>={<EM>x<SUB>p</SUB></EM>,<EM>y<SUB>p</SUB></EM>,<EM>z<SUB>p</SUB></EM>}<SUP>T</SUP>.<BR>  将<EM>L</EM>=<EM>K</EM>-<EM>P</EM>代入式(4)中,求得机构在操作空间内的主动力<EM>F<SUB>j</SUB></EM>为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_24kd6s2008318151147.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(5)</FONT><P align=left><FONT face=宋体>在本机构中驱动杆的结构如图2所示,它的质心距固定端的距离为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_vdw4wl200831815124.gif">&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(6)</FONT><P align=left><FONT face=宋体>式中<EM>m</EM>为杆<EM>i</EM>的质量;<BR>  <EM>m</EM><SUB>1</SUB>为杆<EM>i</EM>中摆动杆的质量;<BR>  <EM>l</EM><SUB>01</SUB>为摆动杆质心到固定点B<SUB>i</SUB>的距离;<BR>  <EM>l</EM><SUB>02</SUB>为伸缩杆的质心到活动点A<SUB>i</SUB>的距离;<BR>  <EM>l<SUB>i</SUB></EM>为杆<EM>i</EM>的总长,它是随时间变化的.<BR>因此,机构的动能</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_ffkjll2008318151221.gif">&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(7)</FONT><P align=left><FONT face=宋体>式中<EM>M</EM>为运动平台的质量;<BR>   <EM>θ<SUB>i</SUB></EM>是杆<EM>i</EM>与固定平台的夹角</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_jarc9z2008318151237.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(8)</FONT><P align=left><FONT face=宋体>机构的势能</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_hilq772008318151252.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(9)</FONT><P align=left><FONT face=宋体>将式(7)和(9)代入式(5)得</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_y3xbls200831815138.gif"><FONT face=宋体>(10)</FONT><P align=left><FONT face=宋体>在式(10)中,忽略了各关节中摩擦力和空间平动机构对驱动力的影响.借助于雅可比矩阵<EM>J</EM>,可将操作空间的主动力转换成关节空间的关节驱动力</FONT></P><P align=left><FONT face=宋体><EM>τ<SUB>j</SUB></EM>=<EM>J<SUP>T</SUP>F<SUB>j</SUB></EM>   <EM>j</EM>=1,2,3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </FONT><FONT face=宋体>(11)</FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_e5s4kw2008318151330.gif"></P><P align=center><FONT face=宋体><STRONG>图2 驱动杆结构示意图</STRONG></FONT></P><P align=left><STRONG><FONT face=宋体>4 修磨轨迹规划</FONT></STRONG></P><P align=left><FONT face=宋体>  钢坯修磨机器人是为了解决钢坯局部修磨自动化而设计的专用机器人,具体修磨轨迹如图3所示.</FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_vh261x200831815149.gif"></P><P align=center><FONT face=宋体><STRONG>图3 钢坯修磨轨迹</STRONG></FONT></P><P align=left><FONT face=宋体>  本文采用直角坐标空间的规划方法进行轨迹规划.由于钢坯修磨采用恒压力,为避免产生磨削烧伤,保证修磨质量,作业时要求运动平台携带砂轮匀速通过全部作业轨迹(两端除外).这就需要对拐角轨道进行特殊的光滑设计,以保证砂轮能够高速平稳地通过这一区域.如采用普通圆弧接续的拐角轨道将能够保证运动速度连续,但在直线与圆弧接点处,由于曲率的突然变化而造成了加速度的突变,从而导致了修磨运动的不稳定.为了解决这个问题,这里采用了回旋曲线接续的拐角轨道<SUP>[3]</SUP>.回旋曲线具有在接续点处曲率连续变化、一阶和二阶导数连续、拐角轨道光滑等特点,可使砂轮在不减速的情况下平稳地通过拐角轨道,特别适合机器人高速度和高精度的控制要求.<BR>  在图3中,修磨作业从A点开始到H点结束,AC、DE和FH为直线段,CD和EF分别为两段以对偶形式连接的回旋曲线.<BR><STRONG>4.1 直线轨迹规划</STRONG><BR>  在轨迹曲线AH中,AB和GH分别为加速和减速阶段,其它部分均为匀速运动.加速和减速阶段采用四次多项式曲线与匀速运动阶段进行拟合,四次多项式曲线方程如下:</FONT></P><P align=left><FONT face=宋体><EM>P</EM>=<EM>K</EM><SUB>1</SUB><EM>t</EM><SUP>4</SUP>+<EM>K</EM><SUB>2</SUB><EM>t</EM><SUP>3</SUP>+<EM>K</EM><SUB>3</SUB><EM>t</EM><SUP>2</SUP>+<EM>K</EM><SUB>4</SUB><EM>t</EM>+<EM>K</EM><SUB>0&nbsp;&nbsp;&nbsp;&nbsp; </SUB></FONT><FONT face=宋体>(12)</FONT></P><P align=left><FONT face=宋体>式中<EM>K</EM><SUB>0</SUB>是每段起始点的坐标.设加(减)速过程所用时间均为<EM>t<SUB>f</SUB></EM>,边值条件为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_ulybsn2008318151437.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(13)</FONT><P align=left><FONT face=宋体>由上述边值条件可求出式(12)的各个系数并代回原式中得</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_p8ixzy2008318151443.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(14)</FONT><P align=left><FONT face=宋体>当<EM>V</EM><SUB>1</SUB>=0,<EM>V</EM><SUB>2</SUB>=<EM>V</EM><SUB>匀</SUB>时,为启动加速的<EM>AB</EM>阶段,运动方程为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_satgmc2008318151534.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(15)</FONT><P align=left><FONT face=宋体>当<EM>V</EM><SUB>1</SUB>=<EM>V</EM><SUB>匀</SUB>,<EM>V</EM><SUB>2</SUB>=0时,为减速停止的GH阶段,运动方程为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_zwbawg2008318151554.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(16)</FONT><P align=left><FONT face=宋体>当<EM>V</EM><SUB>1</SUB>=<EM>V</EM><SUB>2</SUB>=<EM>V</EM><SUB>匀</SUB>时,为匀速直线运动阶段,其中BC段运动方程为</FONT></P><P align=left><FONT face=宋体><EM>x</EM>=<EM>V</EM><SUB>匀</SUB><EM>t</EM>+<EM>x<SUB>B</SUB></EM>, <EM>y</EM>=<EM>y<SUB>B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </SUB></EM></FONT><FONT face=宋体>(17)</FONT></P><P align=left><FONT face=宋体>DE段</FONT></P><P align=left><FONT face=宋体><EM>x</EM>=-<EM>V</EM><SUB>匀</SUB><EM>t</EM>+<EM>x<SUB>D</SUB></EM>, <EM>y</EM>=<EM>y<SUB>D&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </SUB></EM></FONT><FONT face=宋体>(18)</FONT></P><P align=left><FONT face=宋体>FG段</FONT></P><P align=left><FONT face=宋体><EM>x</EM>=<EM>V</EM><SUB>匀</SUB><EM>t</EM>+<EM>x<SUB>F</SUB></EM>, <EM>y</EM>=<EM>y<SUB>F&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </SUB></EM></FONT><FONT face=宋体>(19)</FONT></P><P align=left><FONT face=宋体><STRONG>4. 2 拐角轨道规划</STRONG><BR>  拐角轨道<EM>CM</EM><SUB>1</SUB><EM>D</EM>和<EM>EM</EM><SUB>2</SUB><EM>F</EM>是两段分别以<EM>M</EM><SUB>1</SUB>和<EM>M</EM><SUB>2</SUB>为对称中心的对偶回旋曲线,其中<EM>CM</EM><SUB>1</SUB>段方程为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_lhthhj2008318151627.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(20)</FONT><P align=left><FONT face=宋体>式中: <EM>a</EM>为长度系数,<SUB><IMG src="http://www.chmcw.com/upload_files/article/20/1_xekny72008318151637.gif"></SUB>;<EM>b</EM>为比例系数,<EM>b</EM>=<EM>φ<SUB>m</SUB></EM>/<EM>t</EM><SUP>2</SUP><EM><SUB>m</SUB></EM>;<EM>φ<SUB>m</SUB></EM>是对接点<EM>M</EM><SUB>1</SUB>处的切线角;<EM>t<SUB>m</SUB></EM>是对接点处的参变量;<EM>φ</EM><SUB>0</SUB>是曲线与直线对接点处的倾角.同理<EM>DM</EM><SUB>1</SUB>段方程为</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_i4nons200831815177.gif">&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(21)</FONT><P align=left><FONT face=宋体>EM<SUB>2</SUB>段</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_uzd5vp2008318151713.gif">&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(22)</FONT><P align=left><FONT face=宋体>FM<SUB>2</SUB>段</FONT></P><IMG src="http://www.chmcw.com/upload_files/article/20/1_hwn57n2008318151727.gif">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT face=宋体>(23)</FONT><P align=left><STRONG><FONT face=宋体>5 计算机仿真</FONT></STRONG></P><P align=left><FONT face=宋体>  为分析在修磨过程中各关节和运动平台的速度、加速度及驱动力的变化,应用MATLAB软件对这些因素进行了计算机仿真.设整个机构及修磨轨迹方程的有关参数如下表.</FONT><STRONG></P><P align=center><FONT face=宋体>仿真参数表</FONT></STRONG></P><CENTER><TABLE border=1><TBODY><TR><TD align=middle><FONT face=宋体 size=2>平台参数</FONT></TD><TD align=middle><FONT face=宋体 size=2>运动方程参数</FONT></TD><TD align=middle colSpan=2><FONT face=宋体 size=2>曲线节点坐标(m)</FONT></TD></TR><TR><TD align=middle><FONT face=宋体><FONT size=2><EM>M</EM>=50kg</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>t<SUB>f</SUB></EM>=0.3s</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>x<SUB>A</SUB></EM>=<EM>x<SUB>E</SUB></EM>=<EM>x<SUB>F</SUB></EM>=0.1</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>y<SUB>A</SUB></EM>=<EM>y<SUB>B</SUB></EM>=y<SUB>C</SUB>=0.1</FONT></FONT></TD></TR><TR><TD align=middle><FONT face=宋体><FONT size=2><EM>m</EM>=25kg</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>t<SUB>m</SUB></EM>=0.19s</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>x<SUB>C</SUB></EM>=<EM>x<SUB>D</SUB></EM>=<EM>x<SUB>H</SUB></EM>=0.199</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>y<SUB>D</SUB></EM>=<EM>y<SUB>E</SUB></EM>=0.12</FONT></FONT></TD></TR><TR><TD align=middle><FONT face=宋体><FONT size=2><EM>m</EM><SUB>1</SUB>=20kg</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>φ<SUB>m</SUB></EM>=π/2</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>x<SUB>B</SUB></EM>=0.109</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>y<SUB>F</SUB></EM>=<EM>y<SUB>G</SUB></EM>=<EM>y<SUB>H</SUB></EM>=0.13</FONT></FONT></TD></TR><TR><TD align=middle><FONT face=宋体><FONT size=2><EM>R</EM>=0.65m</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>φ</EM><SUB>0</SUB>=0</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>x<SUB>G</SUB></EM>=0.190</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>z</EM>=1.0</FONT></FONT></TD></TR><TR><TD align=middle><FONT face=宋体><FONT size=2><EM>r</EM>=0.2m</FONT></FONT></TD><TD align=middle><FONT face=宋体><FONT size=2><EM>V</EM><SUB>匀</SUB>=0.06 m/s</FONT></FONT></TD><TD align=middle><FONT size=2> </FONT></TD><TD align=middle><FONT size=2> </FONT></TD></TR></TBODY></TABLE></CENTER><TABLE width="90%" border=0><TBODY><TR><TD><FONT face=宋体><BR><FONT size=2>  首先对式(1)求导,可求得驱动关节的速度和加速度方程,其仿真结果如图4所示.另外,通过对式(15)~(23)的求导,也可求得运动平台在操作空间的速度和加速度方程,仿真结果如图5所示.从图4~图5中可见在全部修磨过程中,驱动关节和运动平台的速度和加速度都是连续的,可以满足高速修磨作业要求.</FONT></FONT></TD></TR></TBODY></TABLE><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_98n7fk2008318151853.gif"></P><P align=center><FONT face=宋体>a. 各杆的速度由线       b. 各杆加速度曲线<BR><STRONG>图4 驱动关节仿真结果</STRONG></FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_lvjczt2008318151923.gif"></P><P align=center><FONT face=宋体>a. 运动平台沿XY方向的速度   b. 运动平台沿XY方向的加速度<BR><STRONG>图5 仿真结果</STRONG></FONT></P><P align=left><FONT face=宋体>  将表1中参数代入(12)式,可得三杆各移动副的驱动力变化如图6.从图6中可见,各杆驱动力在加、减速和拐角轨迹(回旋曲线)阶段有较大的变动,而在直线轨迹阶段,它的变化比较平稳的.并且在全部修磨过程中,驱动力一直是连续变化的.</FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_ei1ppm2008318151957.gif"></P><P align=center><FONT face=宋体><STRONG>图6 各杆的驱动力</STRONG></FONT></P><P><STRONG><FONT face=宋体>6 结论</FONT></STRONG></P><P><FONT face=宋体>  (1) 通过对一种三自由度并联式钢坯修磨机器人机构的研究,得到了相对简单的运动学正逆解,并应用拉格朗日方法求得了机构的动力学方程的显示解.<BR>  (2) 在钢坯修磨轨迹中,加减速及匀速直线运动阶段采用了四次曲线拟合,而在拐角轨道处应用了对偶形式的回旋曲线.<BR>  (3) 仿真结果表明,各杆和运动平台的速度、加速度及杆的关节驱动力均为连续函数,无突变,能够满足高速修磨作业的要求.</FONT></P>
页: [1]
查看完整版本: 并联式钢坯修磨机器人动力学及轨迹规划研究

中国磨床技术论坛
论 坛 声 明 郑重声明:本论坛属技术交流,非盈利性论坛。本论坛言论纯属发表者个人意见,与“中国磨削技术论坛”立场无关。 涉及政治言论一律删除,请所有会员注意.论坛资源由会员从网上收集整理所得,版权属于原作者. 论坛所有资源是进行学习和科研测试之用,请在下载后24小时删除, 本站出于学习和科研的目的进行交流和讨论,如有侵犯原作者的版权, 请来信告知,我们将立即做出整改,并给予相应的答复,谢谢合作!

中国磨削网