发明名称 基于单叶双曲面各向同性并联机构全局优化设计方法
摘要 本发明公开了一种基于单叶双曲面各向同性并联机构全局优化设计方法,采用复合单叶双曲面型式的并联机构,消除了负载惯性参数的约束,实现了完全各向同性。在此基础上提出用于衡量全域工作空间模态变化程度的全局模态灵敏度指标,采用该指标作为优化目标,通过优化结构参数实现了并联机构的具体参数设计。采用该发明设计的并联机构,不仅实现了控制中心的完全解耦及动态各向同性,而且保证了全域工作空间内的最优性能。本方法消除了传统六自由度并联机构无法实现完全各向同性的缺陷,实现了并联机构的完全解耦及优越的动态特性,使得其在工业应用中无需为了提升控制性能而研究复杂的控制策略,降低了工业应用成本。
申请公布号 CN103605893B 申请公布日期 2016.03.23
申请号 CN201310608105.0 申请日期 2013.11.19
申请人 哈尔滨工业大学 发明人 田体先;姜洪洲;佟志忠;何景峰
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 基于单叶双曲面各向同性并联机构全局优化方法,其特征在于,方法如下:步骤1:确定系统构型根据负载特性M<sub>t</sub>,计算满足完全各向同性的中位雅可比矩阵J<sub>lx0</sub>:<img file="FSB0000149103110000011.GIF" wi="2097" he="182" />式(1)中:p<sub>1</sub>(α)=[‑k<sub>a1</sub>sinα k<sub>a1</sub>cosα k<sub>c1</sub> ‑a<sub>1z</sub>k<sub>a1</sub>cosα+r<sub>1</sub>k<sub>c1</sub>sinα ‑a<sub>1z</sub>k<sub>a1</sub>sinα‑r<sub>1</sub>k<sub>c1</sub>cosα r<sub>1</sub>k<sub>a1</sub>]<sup>T</sup>p<sub>2</sub>(α)=[‑k<sub>a2</sub>sinα ‑k<sub>a2</sub>cosα k<sub>c2</sub> ‑a<sub>2z</sub>k<sub>a2</sub>cosα‑r<sub>2</sub>k<sub>c2</sub>sinα a<sub>2z</sub>k<sub>a2</sub>sinα‑r<sub>2</sub>k<sub>c2</sub>cosα ‑r<sub>2</sub>k<sub>a2</sub>]<sup>T</sup><img file="FSB0000149103110000012.GIF" wi="1268" he="164" />复合单叶双曲面的特征参数包括:喉部半径r<sub>1</sub>及r<sub>2</sub>,双曲面中心距a<sub>1z</sub>及a<sub>2z</sub>,系数c<sub>1</sub>及c<sub>2</sub>;负载特性M<sub>t</sub>=[m<sub>x</sub> m<sub>y</sub> m<sub>z</sub> I<sub>xx</sub> I<sub>yy</sub> I<sub>zz</sub>]应满足:m<sub>x</sub>=m<sub>y</sub>=m<sub>z</sub>=m,I<sub>xx</sub>=I<sub>yy</sub>;各参数计算过程如下:定义喉部半径比<img file="FSB0000149103110000013.GIF" wi="155" he="144" />(1)当n=1时:<img file="FSB0000149103110000014.GIF" wi="1534" he="174" />(2)当n≠1时:<img file="FSB0000149103110000015.GIF" wi="1039" he="326" /><img file="FSB0000149103110000016.GIF" wi="581" he="112" /><img file="FSB0000149103110000017.GIF" wi="629" he="175" /><img file="FSB0000149103110000018.GIF" wi="1150" he="191" /><img file="FSB0000149103110000021.GIF" wi="496" he="186" />根据设计要求选取喉部半径比n及角度α,运用公式(1)计算中位雅可比矩阵J<sub>lx0</sub>=J<sub>lx</sub>(α),其中0≤n≤1,<img file="FSB0000149103110000022.GIF" wi="238" he="122" />步骤2:选择优化参数工程设计中固定支腿长度l,并通过调节支腿布置方式来获得具体的并联机构结构设计参数;支腿长度l:<img file="FSB0000149103110000023.GIF" wi="377" he="162" />l<sub>scale</sub>为特征尺度;可优化参数包括:平台中心距d、上平台高度h;各参数选取原则为:平台中心距d:位于两组双曲面中的上平台中心与解耦中心距离,当n=1时,d=0;平移高度h:上下平台沿z向平移的高度;本发明中采用单参数优化,每次优化取上述各参数其中之一;步骤3:结构参数计算对并联机构的结构参数设计实质为根据步骤1中获得的中位雅可比矩阵J<sub>lx0</sub>求取上下平台铰点空间阵A和B;上平台铰点空间阵:A=[a<sub>1</sub> a<sub>2</sub> a<sub>3</sub> a<sub>4</sub> a<sub>5</sub> a<sub>6</sub>] 下平台铰点空间阵:B=[b<sub>1</sub> b<sub>2</sub> b<sub>3</sub> b<sub>4</sub> b<sub>5</sub> b<sub>6</sub>] a<sub>i</sub>为上平台各铰点空间矢量,b<sub>i</sub>为下平台各铰点空间矢量i=1,2…6;具体设计过程如下:(1)采用步骤1求得的中位雅可比矩阵J<sub>lx0</sub>,代入公式(2)中,取出中位各支腿单位矢量l<sub>ni0</sub>和矢量矩V<sub>i0</sub>;<img file="FSB0000149103110000031.GIF" wi="1602" he="494" />(2)计算上下平台铰点空间阵A和B:<img file="FSB0000149103110000032.GIF" wi="720" he="113" /><img file="FSB0000149103110000033.GIF" wi="593" he="206" />当l<sub>ni0</sub>(3)=0时,k<sub>1i</sub>=0;<img file="FSB0000149103110000034.GIF" wi="317" he="169" />当l<sub>ni0</sub>(3)=0时,k<sub>2i</sub>=0;a<sub>i</sub>=p<sub>i</sub>+(k<sub>1i</sub>+k<sub>2i</sub>)l<sub>ni0</sub>L<sub>0</sub>=l·[1<sub>n10</sub> l<sub>n20</sub> l<sub>n30</sub> l<sub>n40</sub> l<sub>n50</sub> l<sub>n60</sub>] B=A‑L<sub>0</sub>步骤4:建立模态灵敏度函数<img file="FSB0000149103110000035.GIF" wi="1704" he="244" />模态灵敏度矩阵:<img file="FSB0000149103110000036.GIF" wi="1446" he="399" />其中:<img file="FSB0000149103110000041.GIF" wi="1543" he="179" /><img file="FSB0000149103110000042.GIF" wi="1552" he="184" /><img file="FSB0000149103110000043.GIF" wi="1666" he="174" /><img file="FSB0000149103110000044.GIF" wi="1628" he="187" /><img file="FSB0000149103110000045.GIF" wi="1604" he="169" /><img file="FSB0000149103110000046.GIF" wi="1606" he="188" />ε为一微小摄动量,取ε=10<sup>‑5</sup>;<img file="FSB0000149103110000047.GIF" wi="715" he="100" /><img file="FSB0000149103110000048.GIF" wi="730" he="108" />G<sub>T</sub>(sx)为广义频率阵,U(sx)<sub>i</sub>为模态解耦阵U(sx)的第i列,U(sx)为广义频率阵G<sub>T</sub>(sx)的模态解耦阵,雅可比矩阵J<sub>lx</sub>(sx):<img file="FSB0000149103110000049.GIF" wi="1608" he="408" />其中:<img file="FSB00001491031100000410.GIF" wi="560" he="88" />l<sub>ni</sub>为各支腿空间单位矢量i=1,2…6;<img file="FSB00001491031100000411.GIF" wi="745" he="235" />c=[x y z]<sup>T</sup><img file="FSB0000149103110000051.GIF" wi="1496" he="264" />步骤5:计算最小灵敏度运用公知的黄金分割法寻找最小值:ρ<sub>min</sub>=minρ(d,h);步骤6:生成优化曲线根据优化变量的不同,绘制优化曲线<img file="FSB0000149103110000052.GIF" wi="346" he="149" />确定优化目标阈值,0≤f≤3db,选取符合目的的优化参数;步骤7:校核’检验设计的结构参数是否存在干涉,若存在,返回步骤2重新修改设计参数进行优化;步骤8:结束。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号哈尔滨工业大学机电工程学院