发明名称 肘关节驱动安装结构及其优化设计方法
摘要 本发明公开了一种基于贯通式直线电机驱动的经过优化的肘关节驱动安装结构及其优化设计方法。包括以下内容:确定贯通式直线电机驱动的肘关节的安装结构形式;确定设计中待优化变量;建立肘关节力学参数的动态数学模型;建立约束函数与目标函数;采用有约束优化算法求得本结构的综合最优解;输出优化设计结果及其运动仿真图形。本发明从电机安装空间,电机传动效率,手臂运动范围,需要的最大输出力矩等多方面的因素考虑,求得一个综合最优解。
申请公布号 CN101785717B 申请公布日期 2011.09.28
申请号 CN201010109902.0 申请日期 2010.02.06
申请人 山东科技大学 发明人 樊炳辉;孙高祚;孙爱芹;江浩;张志献;王传江;徐文尚;邹吉祥;樊东哲
分类号 A61F2/70(2006.01)I;A61F2/58(2006.01)I;B25J17/02(2006.01)I 主分类号 A61F2/70(2006.01)I
代理机构 代理人
主权项 1.一种肘关节驱动安装结构的优化设计方法,所述的肘关节驱动安装结构是一种演化的摆动导杆机构,它由大臂、小臂和贯通式直线步进电机组成;其中,大臂和小臂用肘关节铰链铰接在一起,大臂上设有上铰链支座,小臂上设有下铰链支座,贯通式直线步进电机的机体铰接在大臂的上铰链支座上,贯通式直线步进电机的螺杆前端铰接在小臂的下铰链支座上;其特征在于,设计步骤如下:第一步:确定设计变量及其它参数的选择设上铰链支座的中心点为O<sub>1</sub>、下铰链支座的中心点为O<sub>2</sub>、肘关节铰链的中心点为O<sub>3</sub>、小臂的重心在P点处、当小臂处于初始铅垂位置时过O<sub>3</sub>点的垂线与过P点的水平线交与E点、当小臂处于初始铅垂位置时过O<sub>3</sub>点的水平线与过O<sub>1</sub>点的垂线交与C点、当小臂处于初始铅垂位置时过O<sub>3</sub>点的水平线与过O<sub>2</sub>点的垂线交与D点、当小臂处于初始铅垂位置时过O<sub>1</sub>点的垂线与过O<sub>2</sub>点的水平线交与B点、小臂处于初始铅垂位置时P点在水平方向上偏离O<sub>3</sub>一定值距离EP、小臂处于初始铅垂位置时P点在高度方向上偏离O<sub>3</sub>一定值距离EO<sub>3</sub>;取C、O<sub>1</sub>两点间的距离CO<sub>1</sub>长度为设计变量x<sub>1</sub>(单位:mm);取B、C两点间的距离BC长度为设计变量x<sub>2</sub>(单位:mm);取C、O<sub>3</sub>两点间的距离CO<sub>3</sub>长度为设计变量x<sub>3</sub>(单位:mm);取D、O<sub>3</sub>两点间的距离DO<sub>3</sub>长度为设计变量x<sub>4</sub>(单位:mm);第二步:根据具体设计要求确定x<sub>1</sub>、x<sub>2</sub>、x<sub>3</sub>、x<sub>4</sub>的取值范围并进行初始化赋值;第三步:在小臂运动范围内,建立其重力矩M<sub>z</sub>(单位:Nm)、驱动力矩M<sub>d</sub>(单位:Nm)力学参数的动态数学模型:<maths num="0001"><![CDATA[<math><mrow><msub><mi>&phi;</mi><mn>0</mn></msub><mo>=</mo><mi>arctg</mi><mrow><mo>(</mo><mfrac><mi>EP</mi><msub><mi>EO</mi><mn>3</mn></msub></mfrac><mo>)</mo></mrow></mrow></math>]]></maths><img file="FSB00000535536200022.GIF" wi="204" he="61" /><maths num="0002"><![CDATA[<math><mrow><msup><mi>a</mi><mo>,</mo></msup><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mi>x</mi><mn>4</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mi>b</mi><mo>=</mo><msqrt><msubsup><mi>x</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>3</mn><mn>2</mn></msubsup></msqrt></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>c</mi><mo>=</mo><msqrt><msubsup><mi>x</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>4</mn><mn>2</mn></msubsup></msqrt></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mn>0</mn></msub><mo>=</mo><mi>arccos</mi><mrow><mo>(</mo><mfrac><mrow><msup><mi>b</mi><mn>2</mn></msup><mo>+</mo><msup><mi>c</mi><mn>2</mn></msup><mo>-</mo><msup><mi>a</mi><mrow><mo>,</mo><mn>2</mn></mrow></msup></mrow><mrow><mn>2</mn><mi>bc</mi></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>a</mi><mo>=</mo><msqrt><msup><mi>b</mi><mn>2</mn></msup><mo>+</mo><msup><mi>c</mi><mn>2</mn></msup><mo>-</mo><mn>2</mn><mo>&CenterDot;</mo><mi>b</mi><mo>&CenterDot;</mo><mi>c</mi><mo>&CenterDot;</mo><mi>cos</mi><mi>&alpha;</mi></msqrt></mrow></math>]]></maths>p=(a+b+c)/2<maths num="0007"><![CDATA[<math><mrow><msub><mi>L</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mo>&CenterDot;</mo><msqrt><mi>p</mi><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mi>a</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mi>b</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mi>c</mi><mo>)</mo></mrow></msqrt></mrow><mrow><mn>1000</mn><mo>&CenterDot;</mo><mi>a</mi></mrow></mfrac></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>L</mi><mn>2</mn></msub><mo>=</mo><mi>sin</mi><mi>&phi;</mi><mo>&CenterDot;</mo><msqrt><msubsup><mi>EO</mi><mn>3</mn><mn>2</mn></msubsup><mo>+</mo><msup><mi>EP</mi><mn>2</mn></msup></msqrt><mo>/</mo><mn>1000</mn></mrow></math>]]></maths>M<sub>d</sub>=F·L<sub>1</sub>M<sub>z</sub>=m·g·L<sub>2</sub>其中<img file="FSB000005355362000210.GIF" wi="34" he="42" />为自变量(单位:度),是小臂相对于肘关节转角为0度初始位置时的运动增量数值,取值范围是0°-135°;φ是PO<sub>3</sub>连线与过O<sub>3</sub>点垂线间的夹角(单位:度),数值为负时,夹角在垂线左侧,数值为正时,夹角在垂线右侧,运动中是变量;c=O<sub>2</sub>O<sub>3</sub>,是定值(单位:mm);b=O<sub>1</sub>O<sub>3</sub>,是定值(单位:mm);a=O<sub>1</sub>O<sub>2</sub>,是变量(单位:mm),等价于直线步进电机螺杆伸出长度;<img file="FSB000005355362000211.GIF" wi="485" he="79" />运动中是变量(单位:度),α<sub>0</sub>是小臂在初始铅垂位置时∠O<sub>1</sub>O<sub>3</sub>O<sub>2</sub>的角度,可根据各参数利用几何和三角函数知识编程求得;L<sub>1</sub>是动力臂(单位:m);L<sub>2</sub>是重力臂(单位:m);F是直线步进电机额定推力(单位:N);m为小臂重量(单位:kg);g为重力加速度;第四步:根据直线电机安装空间、各个连杆的安装空间、小臂运动范围、不能发生运动干涉、传动效率各方面要求,设计变量的约束函数如下:g<sub>1</sub>(x)=10-x<sub>1</sub>≤0g<sub>2</sub>(x)=x<sub>1</sub>-55≤0g<sub>3</sub>(x)=10-x<sub>2</sub>≤0g<sub>4</sub>(x)=x<sub>2</sub>-40≤0g<sub>5</sub>(x)=40-x<sub>3</sub>≤0g<sub>6</sub>(x)=x<sub>3</sub>-120≤0g<sub>7</sub>(x)=40-x<sub>4</sub>≤0g<sub>8</sub>(x)=x<sub>4</sub>-x<sub>3</sub>≤0g<sub>9</sub>(x)=20°-α<sub>0</sub>≤0g<sub>10</sub>(x)=α<sub>0</sub>-50°≤0第五步:建立以直线电机驱动力在小臂活动范围内产生的驱动力矩与小臂重力矩与之差绝对值的最大值达到最小的目标函数:ΔM<sub>i</sub>=abs(M<sub>di</sub>-M<sub>zi</sub>)f(x)=max(ΔM<sub>i</sub>)其中:M<sub>di</sub>表示肘关节在运动过程中某直线电机推动力产生的驱动力矩在小臂处于不同位置时的数值,i=1......n;M<sub>zi</sub>表示肘关节在运动过程中小臂重力矩在小臂处于不同位置时的数值,i=1......n;ΔM<sub>i</sub>表示肘关节在运动过程中驱动力矩与小臂重力矩之差的绝对值在小臂处于不同位置时的数值,i=1......n,一般驱动力矩未必在任何工作位置都要大于小臂重力矩,残余的重力矩可采用某种重力矩平衡装置来解决;f(x)是ΔM<sub>i</sub>中的最大值,是目标函数,优化设计的结果是使得f(x)达到最小,即,使得小臂在要求的状态和活动范围内,其驱动力矩与小臂的重力矩之差绝对值的最大值达到最小;第六步:根据设计变量、约束函数、小臂力学参数的动态数学模型及目标函数编制优化设计的计算机程序,并输入计算机进行运行,采用有约束的优化设计算法对各个设计变量x<sub>1</sub>、x<sub>2</sub>、x<sub>3</sub>、x<sub>4</sub>进行优化计算,直至达到期望的优化值;第七步:输出优化设计计算结果及其运动仿真图形。
地址 266510 山东省青岛经济技术开发区前湾港路579号山东科技大学