仿蝗虫跳跃机器人的后腿跳跃机械结构设计_徐成宇
第39卷第4期长春理工大学学报(自然科学版)
Journal of Changchun University of Science and Technology (Natural Science Edition )
Vol.39No.4仿蝗虫跳跃机器人的后腿跳跃机械结构设计
徐成宇,朱人杰,张宝庆,曹晓捷,张连青,袁慧静,熊梦强
(长春理工大学
摘
机电工程学院,长春130022)
要:依据蝗虫身体结构及其运动特点,采用齿轮、凸轮传动,弹簧储能等机构部件,构建仿蝗虫跳跃机器人三维模
型。用摆动从动件凸轮的柔性膝关节模拟半月板的能力释放,建立了以基圆半径最小作为优化目标,按压力角优化的数学模型,确定了仿蝗虫机器人后腿跳跃机构尺寸。仿真结果表明:通过MATLAB 绘制凸轮的理论与实际廓线,使体积在优化后降低了22%,通过ADAMS 的接触应力分析,确定了凸轮副最大承载力不应低于15KN 。此设计研究确定了仿蝗虫机器人的后腿跳跃机械结构尺寸,减小蓄力凸轮体积,为其下一步系统结构优化设计奠定基础。关键词:仿蝗虫机器人;凸轮廓线优化;柔性膝关节;ADAMS ;MATLAB 中图分类号:
TP242
文献标识码:A
文章编号:1672-9870(2016)04-0067-06
Research on the Driving Characteristics of the Take-off
Procedure and Modelling of Locust-like Hopping Robot
XU Chengyu ,ZHU Renjie ,ZHANG Baoqing ,CAO Xiaojie ,ZHANG Lianqing ,YUAN Huijing ,XIONG Mengqiang
(School of Mechatronic Engineering ,Changchun University of Science and Technology ,Changchun 130022)
Abstract :According to the body structure and movement characteristics of locust ,adopting gears ,cam
transmissions
and institutions such as the spring energy storage components ,a 3D prototype of locust-like hopping robot is present-ed. In order to simulating the meniscus to release the power by using oscillating follower cam flexible knee joint ,a mathematical model for optimizing pressure angle is established to minimize the radius of base circle. Finally the jump-ing mechanism size of locust-like robot ’s hind legs is determined. The simulation results shows that :drawing the theo-retic contour line and actual contour line through MATLAB reduces the volume of cam by 22%,the cam ’s maximum bearing capacity no less than 15KN is determined by the ADAMS analysis of the contact stress. This design research has identified the jumping mechanism size of locust-like robot ’s hind legs ,decreased the volume of collected-power cam ,and provided basis for system construction optimal design in next step.
Key words :locust-like robot ;contour line optimization ;flexible knee ;ADAMS ;MATLAB
仿生跳跃机器人涉及了诸如仿生学、机械学和控制学等多个学科,通过从生物体灵巧的运动机构和机敏的运动模式上获取灵感,起到提升跳跃机器人性能的作用,故逐渐成为了热门研究方向。随着研究的深入,这种技术现在大多被用来在外星探索领域,再加之月球与火星表面重力加速度大大小于地球这一因素,跳跃机器人可以轻松越过数倍甚至数十倍于自身尺寸的障碍物,因此越障这一功能在跳跃机器人上广泛运用。近年来,世界各知名科研
收稿日期:2016-02-24
院所,大学的科学家、学者、教授们正试图研制能模
仿动物运动机理的机构模型,例如:美国宇航局(NASA )下属的喷气动力学实验室(JPL )与加利福尼亚技术学院联合研制的仿青蛙弹跳机,在扩大漫游车在星际探索中崎岖地形地貌下的活动范围起到
[1]
作用;Hyon S H 和Mita T 研究设计出的’KenKen ’机器人[2]意在通过模仿肌腱在跑和跳过程中起到的作用,证明了此结构对腿型机器人小跑步态的有效性。国内的西北工业大学葛文杰教授开展了对于仿
基金项目:吉林省科技发展计划资助项目(KYC-JC-XM-2014-025);吉林省大学生创新创业训练计划项目(2014S032)作者简介:徐成宇(1955-),男,博士,讲师,E-mail :[email protected]
长春理工大学学报(自然科学版)2016年
袋鼠跳跃机器人的动力学、运动学特性的深入研究
[3]
。Heitler W J [4]通过动物身体结构角度分析指出,在蝗虫跳跃时,骨骼结构为杠杆,关节结构为支点,而肌肉、肌腱则提供动力;Burrows M 等人[5]则
腿部肌肉与半月板结构在起跳过程中发生变形,因
此引入了柔性机构,利用它的这一特点来模拟蝗虫肌肉储能作用,分析蝗虫起跳阶段中柔性体对起跳力、力矩的影响。其不仅能减小对地面的振动和冲击,还能合理有效的储存和释放能量,从而达到低能耗的目的。
通过前述分析,所设计的仿蝗虫跳跃机器人,如图2所示。其主体由后腿弹跳机构、前腿爬行减震机构、翅膀结构所组成。跳跃机器人采用后腿弹跳工作原理,其机构由凸轮机构、六连杆机构和弹簧装置构成,图3为总体机构运动简图,通过蜗轮蜗杆机构、齿轮机构传动,使凸轮机构的蓄力,依靠弹簧的瞬间释放给地面很大的作用力,从而获得很大加速度带动机体腾空而起。前腿爬行减震机构为多自由度连杆机构,配有减震弹簧,减轻蝗虫机器人落地重
量。
对蝗虫后腿的关节组织进行深入研究,采用解剖学方法揭示了后腿关节在起跳阶段的运动变化。生物学研究表明:蝗虫后腿的膝关节是跳跃运动核心,其中的半月板结构是跳跃的重要保证,一次成功的起跳需要后腿足够大的蹬地力和与之相对应的速度
[6]
。目前,仿生跳跃机器人的研制大多都将关注点放在了跳跃运动机理和机构实现这两个方面,其研究的模型与地面多为点接触,缺乏对跳跃系统中单个部件能量变化的研究,忽略其变化对整个跳跃系统的影响,本文着眼于仿蝗虫机器人的关键部位——膝关节,提出了采用弹簧作为储能元件、凸轮作为蓄力和释放装置的仿生关节的设计方案,建立柔性膝关节,并对凸轮进行了优化设计。
1模型建立
蝗虫的跳跃运动过程分为着地阶段和腾空阶段,着地阶段可分为起跳与落地两个阶段。蝗虫的腿部肌肉的能量变化主要反映在着地阶段,而跳跃能力的高低主要体现在起跳阶段。腾空过程中身体质心变化较小,运动相对平稳。因此,本文着重对于仿蝗虫跳跃机器人着地阶段中的起跳阶段进行研
究。
图2
仿蝗虫机器人总体结构图
图1蝗虫外部结构形态图
蝗虫的主要肌肉,其一为胫骨伸展肌肉,引起腿的伸开,另一块为曲肌,使得腿收起。如图1所示,其腿节与胫节相连,可在一定方向内转动。腿节上肌肉一端连接在腿节,另一端与腿节的腱筋相连,腱筋则通过与腿节和胫节的关节连接到达胫节和跗节。研究者发现,在蝗虫后足腿节与胫节的连接处有一个半月板装置,类似于弹簧,蝗虫起跳前通过腿节的肌肉拉动半月板,起跳时突然释放其中的能量,使胫节快速后摆,完成蹬地动作。
为了更加真实模拟蝗虫的上述特性,特别是其
图3总体机构运动简图
2柔性膝关节设计
如前所述,蝗虫的后腿的膝关节受到腿节处的
伸肌、屈肌即关节内半月板的共同作用。起跳时,屈肌突然放松,立即释放出半月板积聚的能量,使腿部快速伸展。根据机器人膝关节和腿部的运动特点,
第4期徐成宇,等:仿蝗虫跳跃机器人的后腿跳跃机械结构设计
蝗虫的半月板结构使用线性弹簧机构来模拟,从而提出柔性膝关节模型,用于获得较大的收缩角度,同时产生较大的蹬地力。如图4所示,通过凸轮轮廓的特殊设计,在凸轮装置带动从动件摆过一定角度之后,凸轮与从动件突然解除约束,使得从动件在弹簧的带动下快速回摆,从而实现了收腿与摆腿的动
作。
轮基圆半径、摆杆长度、凸轮轴和摆杆回转轴的中心距等基本尺寸。其中需重点考虑压力角对于传动效率的影响,由于此处凸轮机构的升程起蓄力作用,载荷较大,若压力角过大,会导致蓄力过程中的摩擦损耗增大,
有卡死的情况发生,故优化设计的目的是在给定许用压力角条件下,优化各参数,使整个机构的尺寸尽可能小,结构上更加紧凑。
图4弹簧储能元件
图6摆动从动件平面凸轮
2.1蓄力凸轮运动规律确定
图6为摆动从动件平面凸轮的基本构型,令凸轮转动方向同从动件的转动方向相反。凸轮理论轮摆杆初始位置角廓的基圆半径r 0与滚子半径r r 、
摆杆长度L 、凸轮轴和摆杆回转轴的中心距C β0、有关,即
(r 0+r r ) 2=L 2+C 2-2LC cos β0
相较于一般凸轮机构的推程、回程以及休止期,
柔性关节中只需休止期和推程。推程时,从动件带动胫节朝内回收,回程时从动件和凸轮廓线脱离,使从动件完成一个类似急回的动作,而休止期的设置是为了急回之后到下一次的推程过渡地更加平稳。
凸轮从动件的运动规律,选用等加速等减速运动。
推程段从动件运动方程为:
2m 2æϕ÷
φç0≤φ≤ö=βèøϕ推程段从动件运动方程为:2m 2ϕ-ϕ-φ(≤φ≤ϕ) =()m ββ
ϕ(3)
此处取C =55mm ;r 0=15mm ;L =42mm 。则凸轮理论廓线的直角坐标方程为:
ìx =C sin φ-L sin(φ+β0+β1) í=C cos φ-L cos(φ++)
β0β1îy
(4)
(1)
由于理论廓线与实际廓线在法线方向上的距离故可得到实际廓线处处相等,且等于滚子半径r r ,的表达式为:
ì
ï
ïx ' =x +r
ïïïíïï'
ïy =y -r
ïïî
d (2)
此处取βm =80°,ϕ=180°,另取凸轮休止角
即可得到摆动从动件休止期和推程的运ϕs =180°,
动规律线图(图5):
(摆动的最大角度),ϕ为凸式中:βm 为从动件升程
轮推程角,β为从动件摆角。
(5)
其中:
d x =C cos φ-L cos(φ++) æ1+d 1ö
÷β0β1çøè
d 1öd æ
=-C sin φ+L sin(φ+β0+β1) ç1+÷øè
图5摆动件运动规律线图
2.2凸轮轮廓线设计
根据从动件的运动规律,采用优化设计确定凸
取r r =10mm ,通过MATLAB 绘图,可得凸轮廓线,如图7所示:
长春理工大学学报(自然科学版)2016
年
选取多目标优化函数FMINCON 函数[8]求解具有不等式约束的非线性规划问题,其标准形式为:
{
min F (X )
s . t . G (X ) ≤0
(14)
4后腿跳跃机构尺寸确定
r 0=10mm ;L =36.9mm ,此时的凸轮廓线如图8所
图7蓄力凸轮廓线
经过计算的最优组合为:r r =5mm ,βm =110°,[α]=45°,C =48.8mm ,d s =2mm
,
示和优化前后各量的对比见表1。
3约束条件的确定及优化
根据文献[7-8]介绍的摆动从动件凸轮按压力
角优化设计的方法[7],优化目标为基圆半径最小要求,可得到基圆半径的目标函数为:
F (φ)=(R 0-R a )
2
=L 2+1-2L cos β0-2R a L 2+1-2L cos β0+R a 2
(6)
1.5d s +12
,d s 为凸轮轴直径,此处取d s =2mm ,这样是为了保证优化后最小基圆半径不
式中:R a =
图8表1
对比量
L /mm
优化后凸轮廓线
优化前后各量值的对比
C /mm
βm /°
r r /mm
小于按结构需要的最小半径。
L 2=
2éæd 2d öæ2αêcos []ç1-÷+çêèøçèë
d 2ùd öú÷÷úøû
2
(7)
优化前优化后
4236.9
5548.8
80110
105
-[α]-β1β0=öæ
ç1-÷èø(8)
约束条件为:
(1)保证从动杆在初始位置和终了位置的压力
角不得超过[α]:
g 1(x )=cos [α]L 2+1-2L cos β0-sin β0≤0
由表中各量对比可知机构体积比原来减小了=22%,适应了小型化的要求。
根据上述的分析,得出仿蝗虫机器人电子样机的整机尺寸:体长1100mm ,体高260mm ,体宽200mm ,前、中足腿长200mm ,后腿胫节300mm 。重
弹簧量约为
5kg 。其中,弹簧刚度系数k =25N/mm,
原长125mm ,最大伸长量100mm 。后腿弹跳机构示意图如图9所示。
(9)
g 2(x )=cos [α]L 2+1-2L cos(β0+βm ) -sin(β0+βm ) ≤0
(10)
(2)凸轮转角限制:
g 3(x ) =-φ≤0
g 4(x ) =φ-2π≤0
(11)(12)(13)
(3)凸轮最小曲率半径限制:
g 5(x ) =3-r r -1≤0
min
图9弹跳机构结构示意图
第4期徐成宇,等:仿蝗虫跳跃机器人的后腿跳跃机械结构设计
5接触应力分析
ABS 材料,设定材料参数(表2)后,
进行运动仿真,得到接触应力曲线如图9所示:
表2
名称ABS
E /GPa
滚子从动件盘形凸轮机构中,滚子与凸轮轮廓表面理论上为线接触,但受到材料的弹性变形的作用,滚子与凸轮表面实际情况是微小的面接触,由弹性力学中赫兹公式可知,接触应力发生在接触区中线上,计算公式[9]为:
σH =
凸轮副材料参数表
μ
硬度65HRC
σb /MPa
0.20.39432
(15)
式中:F n ——凸轮滚子法向载荷(N );b ——滚子与凸轮接触线长度(mm );ρ1,ρ2——凸轮实际轮廓线曲率半径、滚子半径(mm );μ1,μ2——凸轮与E 1,E 2——凸轮与滚子材料的滚子材料的泊松比;
图10接触应力曲线
弹性模量(MPa )。
如图
9所示,F n 为凸轮滚子法向载荷,F q 为滚α为图示点处的压力角,有几何关系可子所受载荷,
知:
F n =F q cos α(16)
从图10中可看出最大应力不超过12KN ,取许
用系数
1.2后取整,故设定的最大承载力不应低于15KN
。
6跳跃结果
图11描述了样机系统在跳跃腾空阶段中的某一瞬间的姿态,从图可看到,样机在仿真过程中姿态基本平稳,与水平面夹角都在±8°以内,如图12所示,它的跳跃曲线的最高点高度接近1750mm 。
图9凸轮受力示意图
凸轮廓线上各点压力角与曲率半径均不相同,
为求解接触应力,需求出各点的压力角和曲率半径。当凸轮转动到φ时,从动件摆角为β1,此点距离转动中心为C ,该点的压力角与理论廓线曲率半径为:
d öæL ç1-÷øè
tan α=cot(β1+β0) -(17)
β1β0æd ö32
(C +ç÷)
èø
2
2
图11仿蝗虫机器人跳跃姿态
æd öæd öC 2+2ç÷+C ç÷
èøèø
凸轮实际廓线曲率半径为
ρ=ρ1-ρ2
ρ=
(18)
图12仿蝗虫机器人跳跃曲线
7
(19)
结论
将凸轮导入ADAMS 中,凸轮与滚子均选用
基于仿生学原理,建立仿蝗虫跳跃机器人的机
构模型,用摆动从动件凸轮的柔性膝关节模拟半月
长春理工大学学报(自然科学版)2016年
板的能力释放,使腿节在弹簧的带动下快速回摆;以[4]Li X Y. Kinetic analysis of flexible robot kangaroo
hopping mechanism [D ]. Northwestern Polytechnic 基圆半径最小作为优化目标,按压力角优化的数学
University ,Xian. 2006.3. 模型,通过MATLAB 绘制凸轮的理论与实际廓线,
[5]Heitler W J.How locust jump [EB/OL][2012-07-21].
使体积在优化后降低了22%,设计了仿蝗虫机器人
Standrews :The University of St Andrews ,2007.
后腿跳跃机构尺寸;通过ADAMS 的接触应力分析,
[6]Burrows M ,Sutton P. Locusts use a composite of
确定了凸轮副最大承载力不应低于15kN ,跳跃机器
resilin and hard cuticle as an energy store for jump-人的跳跃高度为1750mm 。ing and kicking [J ].The Journal of Experimental Bi-参考文献
[1]Hale E ,Schara N ,Burdick J. A minimally actuated
hopping rover for exploration of celestial bodies [C ]. Proceedings of the IEEE International Conference on Robotics and Automation. San Francisco :IEEE Press ,2000:420-427.
[2]Hyon S H ,Mita T. Development of biologically in-spired hopping robot [C ]. Proceedings of the 2002IEEE international conference on robotics automa-tion Washington DC ,2002.
[3]李向阳. 仿袋鼠机器人柔性跳跃机构的动力学分析
[D ]. 西安:西北工业大学,2006.3.
ology ,2012,215(19):3501-3512. [7]
Chen Diansheng ,Yin Junmao ,Zhao Kai. Bionic mechanism and kinematics analysis of hopping robot inspired by locust jumping [J ].Journal of Bionic En-gineering ,2011,8(4):429-439.
[8]吕庸厚,周榴明. 摆动从动杆盘状凸轮机构按压力角的
优化设计[J ]. 无锡轻工业学院学报,1986,5(1):38-54.
[9]杨绘宇,王石刚,梁庆华,等. 仿蝗虫跳跃机器人起跳过
程运动学建模及分析[J ]. 机械设计,2010,27(4):54-57.
[10]黄尚兵,金光,安源,等. 滚子摆动从动件盘形凸轮的
设计及优化[J ]. 工程设计学报,2012,12(6):449-453.
(上接第62页)
还原倾向实验证明,可有效提高产品透过率。
(3)控制光学玻璃熔炼时的氧化还原气氛可以改善玻璃的透过率,通过气氛抑制变价元素易着色价态,使变价元素的吸收带远离使用光谱波长,达到提高透过率的效果。改变熔炼气氛的方法亦有两种,一是向玻璃中引入氧化剂或还原剂,靠氧化剂或还原剂来抑制或促成元素有利价态的形成,二是使用各类气体如O 2、N 2、CO 2、CO 、H 2等制造有利的化学反应条件。
参考文献
[1]李维民. 稀土光学玻璃现状与发展[J ]. 稀土信息,2012,
5(11):19-22.
[2]江朗. 稀土光学玻璃研制总结[J ]. 仪表材料,1980,6
(22):44-45.
[3]冯宏,李维民. 光学玻璃工艺学[M ]. 重庆:中国兵器装
备集团公司西南兵工局教育中心,2007:23-236.
[4]衣德学,翟守元,何文. 电熔锆刚玉砖侵蚀特征的研究
[J ]. 山东建材学院学报,1994,8(4):66-71.
[5]西北轻工业学院. 玻璃工艺学[M ]. 北京:中国轻工业出
版社,2004:216-217.
[6]贺文燕,谷帮新. 氧化态和氧化铁对镧系玻璃光吸收的
影响[J ]. 光电工程,1980.2(2):42-48.
[7]张晓伟,李梅,柳召刚,等. 稀土在玻璃中的作色和紫外
吸收作用的研究进展[J ]. 硅酸盐通报,2009,28(06):1208-1212.
[8]林玉珍,吕桂兰. 稀土光学玻璃着色的研究[J ]. 光学精
密工程,2011,8(22):6-12.
[9]陈新政,王海风. 透远紫外玻璃中痕量铁对紫外透过率
的影响[J ]. 玻璃与搪瓷,2011,39(4):13-15.
[10]贺文燕. 稀土杂质对镧系玻璃光吸收的影响及其测定
方法的研究[J ]. 硅酸盐通报,1983,4(28):39-44.
[11]曹振亚,张泽田. 有色玻璃与特种玻璃[M ]. 四川:四川
科学技术出版社,1987:25-105.
[12]田英良,孙诗兵. 新编玻璃工艺学[M ]. 北京:中国轻工
业出版社,2011:154-173.