发明名称 基于关节力的空间六维力测量方法
摘要 本发明提供的一种基于关节力的空间六维力测量方法,包括上平台和下平台,所述下平台由六个电机分别独立驱动的六个驱动杆与上平台连接,六个电机形成六个关节,所述上平台和下平台分别具有六个连接点,通过建立坐标系、运动学逆解、运动学正解以及Jacobi矩阵的求解;能够在极端环境中对机器人末端的六维力和力矩进行准确测量计算,并且能够适应恶劣的工况环境,可靠性强。
申请公布号 CN104634498A 申请公布日期 2015.05.20
申请号 CN201510057991.1 申请日期 2015.02.04
申请人 重庆大学 发明人 刘飞;赵云花;谢更新;熊辉;詹涵菁;张元勋
分类号 G01L5/16(2006.01)I 主分类号 G01L5/16(2006.01)I
代理机构 北京海虹嘉诚知识产权代理有限公司 11129 代理人 谢殿武
主权项 一种基于关节力的空间六维力测量方法,其特征在于:包括上平台和下平台,所述下平台由六个电机分别独立驱动的六个驱动杆与上平台连接形成并联结构,六个电机形成六个关节,所述上平台和下平台分别具有六个连接点;包括如下步骤:S1.建立坐标系在上平台的中心建立工件坐标系O1‑x1y1z1,在下平台的中心建立基坐标系O‑xyz,Bi和Ai分别为上下平台对应的六个连接点,其中:下平台的各连接点在基坐标系中的向量为(Aix,Aiy,Aiz),上平台的各连接点在基坐标系中的向量表示为(Bix,Biy,Biz),在工件坐标系中的向量表示为(bix,biy,biz),工件坐标系在基坐标系中的位置矢量为O1,li为上平台和下平台对应点的长度,上平台相对于下平台的位姿x(t)、y(t)、z(t)、θ<sub>x</sub>(t)、θ<sub>y</sub>(t)、θ<sub>z</sub>(t)表示为x、y、z、θ<sub>x</sub>、θ<sub>y</sub>、θ<sub>z</sub>;S2.运动学逆解当上平台的位姿改变时,根据平面与平面上点的关系求出此时新点的坐标值Bi,B<sub>i</sub>=Rb<sub>i</sub>+O<sub>1</sub>     (1.1)其中,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>y</mi></msub><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>x</mi></msub></mtd><mtd><mo>-</mo><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>z</mi></msub><mo>+</mo><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><mi>sin</mi><msub><mi>&theta;</mi><mi>y</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>z</mi></msub></mtd><mtd><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>z</mi></msub><mo>+</mo><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>y</mi></msub><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>z</mi></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>y</mi></msub></mtd><mtd><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>y</mi></msub><mi>sin</mi><msub><mi>&theta;</mi><mi>z</mi></msub><mo>+</mo><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>z</mi></msub></mtd><mtd><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>y</mi></msub><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>z</mi></msub><mo>-</mo><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>z</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>y</mi></msub></mtd><mtd><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>y</mi></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>z</mi></msub></mtd><mtd><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>y</mi></msub><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>z</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000667143930000011.GIF" wi="1948" he="298" /></maths>R为上平台姿态的旋转矩阵,O<sub>1</sub>为上平台上的工件坐标系相对于下平台基坐标系中的位移矢量,b<sub>i</sub>为B<sub>i</sub>点在工件坐标系中的位置矢量;驱动杆在基坐标系o‑xyz中的位置矢量为L<sub>i</sub>=B<sub>i</sub>‑A<sub>i</sub>=Rb<sub>i</sub>+O<sub>1</sub>‑A<sub>i</sub>,其中,(i=1,2,...,6),驱动杆的实时长度可表示称被测运动问题位姿参数的函数:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>l</mi><mi>i</mi></msub><mo>=</mo><mo>|</mo><mo>|</mo><msub><mi>L</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo><mo>=</mo><mo>|</mo><msub><mi>B</mi><mi>i</mi></msub><mo>-</mo><msub><mi>A</mi><mi>i</mi></msub><mo>|</mo><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>Rb</mi><mi>i</mi></msub><mo>+</mo><msub><mi>O</mi><mn>1</mn></msub><mo>-</mo><msub><mi>A</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>Rb</mi><mi>i</mi></msub><mo>+</mo><msub><mi>O</mi><mn>1</mn></msub><mo>-</mo><msub><mi>A</mi><mi>i</mi></msub><mo>)</mo></mrow></msqrt><mo>=</mo><msub><mi>h</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>,</mo><msub><mi>&theta;</mi><mi>x</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>y</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>z</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1.2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000667143930000021.GIF" wi="1749" he="135" /></maths>B<sub>i</sub>、b<sub>i</sub>分别为上平台顶点在基坐标系、工件坐标系中的位置向量,A<sub>i</sub>为下平台顶点在基坐标系中的位置向量,将被测物体的六个位姿参数代入到(1.2)式中,可求得驱动杆的长度;S3.运动学正解正解方程由逆解方程变化而来:l<sub>i</sub><sup>2</sup>=(B<sub>i</sub>‑A<sub>i</sub>)(B<sub>i</sub>‑A<sub>i</sub>)<sup>T</sup>,其中B<sub>i</sub>,A<sub>i</sub>分别表示上平台和下平台的连接点在基坐标系中的坐标;f<sub>i</sub>(B<sub>i</sub>)=f<sub>i</sub>(x,y,z,θ<sub>x</sub>,θ<sub>y</sub>,θ<sub>z</sub>)=(Rb<sub>i</sub>+O<sub>1</sub>‑A<sub>i</sub>)(Rb<sub>i</sub>+O<sub>1</sub>‑A<sub>i</sub>)<sup>T</sup>‑l<sub>i</sub><sup>2</sup>=0,i=1,2,...,6   (1.3);首先令x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,x<sub>4</sub>,x<sub>5</sub>,x<sub>6</sub>=(x,y,z,θ<sub>x</sub>,θ<sub>y</sub>,θ<sub>z</sub>)且初始点(x,y,z,θ<sub>x</sub>,θ<sub>y</sub>,θ<sub>z</sub>)=(0,0,0,0,0,0),然后将f<sub>i</sub>(B<sub>i</sub>)(i=1,2,...,6)在B<sub>i</sub>附近进行Taylor展开,取其一阶线性部分可得<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>f</mi><mi>i</mi></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>-</mo><msub><mi>x</mi><msub><mi>k</mi><mn>0</mn></msub></msub><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><msub><mi>f</mi><mi>i</mi></msub></mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>k</mi></msub></mfrac><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1.4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000667143930000022.GIF" wi="1225" he="252" /></maths>式(1.4)中为x<sub>i</sub>的线性方程组,对应的系数矩阵为J<sub>1</sub>且系数矩阵对应的第i行第j列元素为<img file="FDA0000667143930000023.GIF" wi="358" he="259" />通过求系数矩阵J<sub>1</sub>的逆矩阵以及解方程组即可求得上平台的位姿;S4.力Jacobi矩阵的求解机构的速度Jacobi矩阵J为:J=[e<sup>T</sup> (Rb×e)<sup>T</sup>]<sup>‑1</sup>     (1.5)对于六个关节驱动力(或力矩)组成的关节矢量为τ=[τ<sub>1</sub>,τ<sub>2</sub>,τ<sub>3</sub>,τ<sub>4</sub>,τ<sub>5</sub>,τ<sub>6</sub>]<sup>T</sup>关节矢量与运动平台的广义操作力矢量F=[f<sub>1</sub>,f<sub>2</sub>,f<sub>3</sub>,f<sub>4</sub>,f<sub>5</sub>,f<sub>6</sub>]<sup>T</sup>具有如下关系:τ=J<sup>T</sup>(q)F     (1.6)力Jacobi矩阵J<sub>F</sub>是静平衡状态下,关节力向操作力映射的线性关系,即F=J<sub>F</sub>τ,则<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mrow><msub><mi>J</mi><mi>F</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>J</mi><mi>T</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>J</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mi>T</mi></msup><mo>=</mo><mo>[</mo><msup><mi>e</mi><mi>T</mi></msup><msup><mrow><mo>(</mo><mi>Rb</mi><mo>&times;</mo><mi>e</mi><mo>)</mo></mrow><mi>T</mi></msup><msup><mo>]</mo><mi>T</mi></msup><mo>=</mo></mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>e</mi></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mi>Rb</mi><mo>&times;</mo><mi>e</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1.7</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000667143930000032.GIF" wi="1841" he="214" /></maths>
地址 400044 重庆市沙坪坝区沙坪坝正街174号