发明名称 防止模态跃迁的六自由度并联机构全局模态空间控制方法
摘要 本发明公开了一种防止模态跃迁的六自由度并联机构全局模态空间控制方法。本发明采用一种新的矩阵特征值迭代求解算法获得并联机构在全局工作空间内的实时模态解耦矩阵。该迭代求解算法利用模态连续性原理,基于数值迭代思想,采用从中位模态向全局工作空间逐步迭代逼近的方法,沿着位姿变化方向,通过若干步迭代完成模态频率及模态解耦阵的求解,消除了采用传统方法求解造成的模态跃迁现象。本发明具有在六自由度并联机构全局工作空间内实现完全模态解耦的优点,能够使模态空间控制器的适用性增强。在此基础上构造的全局模态空间控制器,扩展了模态空间控制的适用范围,使之适用于六自由度并联机构大范围全局工作空间。
申请公布号 CN103170972B 申请公布日期 2015.01.28
申请号 CN201310054212.3 申请日期 2013.01.31
申请人 哈尔滨工业大学 发明人 田体先;姜洪洲;何景峰;佟志忠
分类号 B25J9/16(2006.01)I 主分类号 B25J9/16(2006.01)I
代理机构 代理人
主权项 一种防止模态跃迁的六自由度并联机构全局模态空间控制方法,并联机构的广义频率矩阵Ω=f(ε)是并联机构位姿ε的函数,当并联机构平台由中位ε<sub>0</sub>处移动到空间某一位姿ε<sub>T</sub>处时,广义频率矩阵需由中位条件下计算得到的Ω<sub>0</sub>=f(ε<sub>0</sub>)实时更新为以新的位姿ε<sub>T</sub>为变量计算得到的广义频率矩阵Ω<sub>T</sub>=f(ε<sub>T</sub>),矩阵Ω<sub>T</sub>的特征值分解Ω<sub>T</sub>=U<sup>T</sup>∑U,其特征在于,方法如下:步骤1:设置初值:包括初始步长Δ<sub>h0</sub>、初始模态解耦阵U<sub>0</sub>及迭代精度e;初始步长Δ<sub>h0</sub>的选取:将Ω<sub>T</sub>分解为Ω<sub>T</sub>=Ω<sub>0</sub>+Δ<sub>T</sub>,取<img file="FSB0000133223650000011.GIF" wi="211" he="118" />N为正整数;初始模态解耦阵U<sub>0</sub>的选取:1)计算广义频率矩阵Ω<sub>0</sub>,Ω<sub>0</sub>=f(ε<sub>0</sub>);2)使用公知的常规数值解法算法SVD求得Ω<sub>0</sub>的特征值分解:<img file="FSB0000133223650000012.GIF" wi="305" he="62" />迭代精度e的选取:依据所需计算精度而定;初次迭代计算中Δ<sub>h</sub>=Δ<sub>h0</sub>,U=U<sub>0</sub>;步骤2:完成矩阵迭代:Ω<sub>i</sub>=Ω<sub>i‑1</sub>+Δ<sub>h</sub>,Δ<sub>t</sub>=Δ<sub>t‑1</sub>+Δ<sub>h</sub>;步骤3:对矩阵Ω<sub>i</sub>进行相似变换H=U<sup>T</sup>Ω<sub>i</sub>U;步骤4:采用公知的雅可比Jacobi方法对矩阵H进行特征值分解H=V<sup>T</sup>∑V;步骤5:计算相似性判别指标<img file="FSB0000133223650000013.GIF" wi="352" he="150" />k为矩阵V及E的列数,k=1…n,n为矩阵Ω<sub>T</sub>的阶次;步骤6:根据判别指标γ判断是否发生模态跃迁,若γ≠E,则说明此时发生了模态跃迁,应减小步长,取<img file="FSB0000133223650000014.GIF" wi="168" he="110" />返回步骤2重新计算;步骤7:若γ=E,说明计算正确,进行模态解耦阵迭代:U=UV;步骤8:判断迭代是否应该结束;e<sub>max</sub>的计算式为:<img file="FSB0000133223650000015.GIF" wi="331" he="92" />Φ=U<sup>T</sup>Ω<sub>T</sub>U;若e<sub>max</sub>≤e则程序结束,U即为矩阵Ω<sub>T</sub>的模态解耦阵,∑为相应的模态频率阵;若e<sub>max</sub>>e且||Δ<sub>t</sub>||>||Δ<sub>T</sub>||,取<img file="FSB0000133223650000016.GIF" wi="229" he="128" />返回步骤2重新计算;若e<sub>max</sub>>e且||Δ<sub>t</sub>||<||Δ<sub>T</sub>||,重置步长Δ<sub>h</sub>=Δ<sub>h0</sub>返回步骤2重新计算。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号哈尔滨工业大学机电工程学院