三相驱动下非光滑振动驱动系统的动力学分析1)

朱诗慧 周震 吕敬 王琪

(北京航空航天大学航空科学与工程学院,北京 100083)

引言

近年来可移动式机器人的发展受到了广泛关注,其运动方式仍然是以轮式或腿式为主[1-2],但是在管道清理、灾后救援等特殊情况下,这些运动方式并不能很好的满足人们的需求.Chernous’ko[3-4]提出了振动驱动移动系统,这类系统依靠内部质量块的周期性振动来驱动外部箱体在有阻力环境中移动,具有结构简单、易于密封化和微型化等优点.

学者们在振动驱动系统的直线平移运动特性方面做了大量研究,其中Sergey 等[5]指出振动驱动系统的定向运动依赖于“非对称的摩擦力”或“非对称的内部驱动方式”; Chernous’ko[6-7]引入非对称驱动方式:两相驱动和三相驱动,分别对内部质量块施加速度控制模式和加速度控制模式,分析了系统实现周期性稳态运动的条件.Li 等[8]和Sobolev 等[9]分别通过倒摆车和含多个转子的振动驱动实验模型验证了直线平移的相关理论.方虹斌等[10-12]运用平均法得到了三相驱动下振动驱动系统在黏性摩擦环境中运动速度的近似解,并从滑移分岔的角度分析了其在库伦干摩擦作用下的运动情况.Du 等[13]通过压电悬臂梁实验对简谐振动驱动系统的黏滑运动进行了实验验证,观察到了明显的黏滞−滑移现象.陈祺等[14-15]基于一维库伦干摩擦模型(为便于计算假设动、静摩擦系数相等),首次分析了系统不同黏滑运动的分岔条件.

目前,对于作平面一般运动振动驱动系统的运动特性研究得较少[16-21],其中Volkva[16]等研究了这样一类振动驱动系统,它借助外部4 个倾斜的支撑足与地面保持刚性接触,内部有两个做正弦运动的质量块,在各向异性(前后运动的摩擦系数不相同)的黏性摩擦力作用下可以实现S 形的折线运动.殷蓬勃等[17]在此基础上分析了三相驱动下系统的平面运动特性,通过调节两个质量块的振动周期比和相位差实现了外部箱体上某点的运动轨迹在直线和圆弧间的切换.占雄等[18-19]引入各向异性的库伦干摩擦模型,研究了一类振动驱动系统(滑道垂直,正弦激励),其外部箱体与地面间保持面接触; 并假设箱体在运动过程中,其纵轴始终与箱体形心的运动轨迹相切;揭示了箱体平动和转动中出现的黏滞−滑移切换现象.

这类振动驱动系统箱体的运动是通过外部的摩擦力实现的.干摩擦的不连续性带来的黏滞−滑移效应的存在[22-23],使得系统的运动情况更为复杂,但充分利用其特性可以有效提高系统运动效率.以往的研究大多采用库伦干摩擦模型分析系统直线平移运动中的黏滞−滑移切换现象[13].在研究箱体作平面一般运动问题时,若采用库伦模型处理多点摩擦问题,数值求解难度较大[24-25].

近年来,LuGre 摩擦模型也得到了广泛的应用[26-27].它可反映库伦干摩擦、黏性摩擦和Stribeck摩擦等多种摩擦特性[28],在降低数值计算难度的同时较为全面地揭示摩擦特性.该模型的参数较多且不易选取,在工程手册上也难以查到,但可借助一些经验公式[29-30]进行选取.若模型参数选取不合适时,会导致计算仿真结果失真.

本文将借鉴文献[14-18]的研究方法研究三相振动驱动系统的动力学特性; 为弥补原有研究中所采用的摩擦模型的不足,文中将采用LuGre 摩擦模型描述箱体与地面接触点的二维摩擦,能更好地反映静摩擦特性对黏滞−滑移运动的影响,且便于计算;利用Lagrange 方程,给出振动驱动系统作平面一般运动的动力学模型和算法; 分析三相驱动参数的变化对系统运动特性的影响,为振动驱动型机器人的设计和控制提供理论支撑.

1 振动驱动系统

1.1 力学模型

设振动驱动系统由质量为M的均质刚性箱体和质量分别为m1,m2的质量块(视为质点)组成,其中质量块在箱体内的两条平行滑道中相对箱体作周期运动,箱体通过三个刚性支撑足Ci(i=1,2,3)始终与地面保持接触,如图1 所示.箱体的长和宽分别为2a和2b,高度不计.长度均为2L的两个平行滑道平行于箱体的对称轴,到对称轴的距离均为l.支撑足到箱体形心O1的距离O1Ci=a,且tan φ=b/a.

图1 三足支撑的振动驱动系统模型Fig.1 The model of vibration-driven system with tripod support

设Oxy为惯性坐标系,O1ξη 为固连在箱体上的动坐标系,θ 表示轴O1ξ 与轴Ox间的夹角,以逆时针为正; (xM,yM) 是箱体形心O1在惯性坐标系中的坐标,则系统的广义坐标可表示为

1.2 三相驱动

三相驱动是指将质量块相对箱体的运动周期T∗分为τ1,τ2,τ3三段,在每个时间段,质量块相对箱体作匀变速直线运动,其相对运动的数学表达式为[3]

其中,wk> 0 (k=1,2,3) 称为驱动参数,T∗=τ1+τ2+τ3.图2 给出了该三相驱动方式下质量块的运动特征.

图2 三相驱动下质量块的运动特征Fig.2 Movement characteristics of a mass under three-phase control

设t时刻质量块相对于箱体的横坐标为ξj(t)(j=1,2),且满足

以质量块在区间[0,L]上的正向运动为例,分析其运动时的参数特征.为避免质量块与箱体间发生碰撞,并满足周期性运动条件,可得

且存在某个时刻t0j有

对于质量块mj(j=1,2),式(2)中的w1,w2,w3分别定义为w1j,w2j,w3j,由式(2)∼式(5)可以导出[3,7]

2 动力学方程及其算法

2.1 箱体接触足的摩擦力

动坐标系O1ξη 与惯性坐标系Oxy的转换矩阵可表示为

由图1 可知,支撑足Ci(i=1,2,3)相对箱体的位置向量可表示为

支撑足在惯性系中的位置向量可表示为

其中rO1=[xM,yM]T,将上式对时间求导可得其速度vCi=.

用二维LuGre 摩擦模型[31-32]描述支撑足与地面间的摩擦力,其表达式为

其中,FTi=[FTix,FTiy]T为支撑足Ci的摩擦力向量,FNi为支撑力.µi=[µix,µiy]T为LuGre 摩擦系数矩阵,其表达式为

其中,zi=[zix,ziy]T(i=1,2,3)是鬃毛的平均变形量列向量;σ0和σ1为鬃毛刚度和阻尼系数,σ2为黏性系数,µc,µs分别为动/静摩擦系数,vs为Stribeck 速度.

2.2 振动驱动系统的动力学方程

设质量块m1,m2在动坐标系中的相对位置向量为

则质量块mj(j=1,2)的绝对位置向量为

由第二类拉格朗日方程可得该系统的动力学方程

其中,Qq为系统广义力列向量.T为系统动能,且T=TM+T1+T2,其具体表达式为

其中JM=M(a2+b2)/3.

系统的广义力Qq可表示为

其中,FT x和FTy是3 个支撑足的摩擦力向点O1简化的主矢FT在x轴和y轴上的投影,MT是这些摩擦力向点O1简化的主矩.

当不考虑箱体高度时,支撑力可由下列方程求得

将求出的支撑力代入式(10) 得到摩擦力,并代入式(16),利用动力学方程(14),通过数值仿真可分析该系统的动力学特性.

3 算例分析

设箱体内的质量块按照图3 给出的4 种方式进行驱动,分别就箱体直线平移和平面一般运动进行动力学分析,算例中系统的基本参数如下[14,20]:a=0.05 m,b=0.03 m,L=0.01 m,l=0.02 m,M=0.1 kg,m1=m2=0.03 kg,g=10 m/s2,µs=0.30,µc=0.25,vs=0.004 m/s,σ0=106/m,σ1=103s/m,σ2=0 s/m.

图3 质量块的不同驱动方式Fig.3 Different methods of mass-driven

3.1 直线平移运动的仿真与分析

首先,设两个质量块的驱动参数相同,即wk1=wk2(k=1,2,3),ξ1=ξ2=ξ,如图3(a)所示.由于驱动的对称性,以及质量分布与受力的对称性,系统沿x轴作直线平移运动.

当wk1=wk2(k=1,2,3),其值取表1 给出的值.由于摩擦的存在,箱体稳态的周期性运动有4 种类型[14,33]:擦滑、穿滑、回滑和不黏(grazingsliding、crossing-sliding、switching-sliding,no-stick).

表1 不同运动类型的驱动参数Table 1 Different drive parameters

图4 给出了箱体4 种运动类型下的速度时间历程图,其中黑虚线是基于LuGre 摩擦模型的仿真结果,红线是基于库伦摩擦模型的仿真结果,篮圈是由文献[15] 得到的理论解析结果.3 种结果完全一致,验证了本文方法的有效性.

图4 四种运动类型箱体移动速度的时间历程图Fig.4 Time history diagram of box moving-velocity in four motion types

图4 四种运动类型箱体移动速度的时间历程图(续)Fig.4 Time history diagram of box moving-velocity in four motion types(continued)

由图4 可以看出,当系统是擦滑运动类型时,箱体时而向前运动,时而保持静止,但不存在倒退现象.而其他运动类型,箱体均存在着倒退现象.

图5 驱动参数与 的关系图Fig.5 The relationship between and the drive parameters

从图5 可以发现,存在明显的脊线.当w1,w2取脊线上的值时,其对应的值大于两侧的值;且有x˙随w2的增大而增大.因此当w2的值给定时(如图5(c)所示),在脊线的交界处,取得最大值.

图6(a)和图6(b)分别展示了w1和w3对的影响,以及箱体的位置和速度随驱动参数变化的区别.在转折点处所对应的驱动参数是箱体由擦滑运动类型转为其他运动类型的临界点,这时系统的移动效率最优,避免了后退运动的出现,减少了能量的消耗.

图6 驱动参数w1,w3 对的影响Fig.6 The effect of drive parameters w1,w3 on the

图6 驱动参数w1,w3 对 的影响(续)Fig.6 The effect of drive parameters w1,w3 on the(continued)

3.2 平面一般运动的仿真与分析

3.2.1 定轴转动

当按照图3(b)所示的模式驱动时,即wk1=−wk2(k=1,2,3),ξ1=−ξ2,由于驱动的反对称性,若初始系统静止,3 个支撑足处摩擦力的主矢FT=0,向箱体形心简化的主矩MT0,箱体绕其形心作定轴转动.同样其转动运动类型也有4 种,即擦滑、穿滑、回滑和不黏,如图7 所示.

图7 四种运动类型箱体转动速度的时间历程图Fig.7 Time history diagram of box rotation-velocity in four motion types

图8 驱动参数与的关系图Fig.8 The relationship betweenand the drive parameters

在图8 中,同样可以看到脊线,当驱动参数取脊线对应的数值时,取得最大值; 脊线对应的驱动参数也是擦滑运动类型与其他运动类型的分界值.

3.2.2 平面一般运动

当按照图3(c) 所示的模式驱动时,即|wk1||wk2|(k=1,2,3),令n=wk2/wk1,箱体形心运动轨迹为圆周.箱体作平面一般运动.

若w11=w31=8 m/s2,w21=30 m/s2,取不同的比值n,箱体形心运动轨迹为半径不同的圆周,如图9 所示.并描绘了n=1.2 及n=−3 时,形心的y坐标与箱体转动角度的关系图.可见,箱体形心运动一周,箱体转动一周.

图10 给出了不同驱动参数和不同摩擦系数情况下,箱体形心运动轨迹的曲率半径r.

图9 不同n 值的形心运动轨迹Fig.9 The centroid-trajectory at different values of n

图10 曲率半径r 与参数比n 的关系图Fig.10 The relationship between radius of curvature r and parameters ratio n

通过上述算例的分析可知,根据不同的摩擦系数,通过调整驱动参数可改变箱体的运动形式和形心运动轨迹.

4 结论

本文研究了一类内部受两个相互平行的三相驱动激励的振动驱动系统,在各向同性的二维干摩擦环境中的平面运动规律.首先,建立了具有3 个足支撑的振动驱动系统的力学模型; 采用二维LuGre 摩擦模型描述支撑足与地面间的摩擦,用拉格朗日方程建立了该系统的动力学方程.然后,通过数值仿真对该系统的动力学特性进行了分析:(1)利用系统做直线平移运动的理论解和基于库伦摩擦的数值仿真解,校验了基于LuGre 摩擦模型的数值仿真结果,表明了本文的建模方法与数值算法的合理性; (2)通过调整驱动参数,实现了箱体的平移运动、定轴转动和平面一般运动,揭示了箱体平动和转动过程中均存在擦滑、穿滑、回滑和不黏4 种运动类型; (3) 不仅分析了驱动参数对箱体移动速度的影响,还分析了其对箱体转动快慢和箱体形心运动轨迹曲率半径的影响.分析结果表明,这类振动驱动系统拥有丰富的黏滞−滑移动力学行为和形式多样的平面运动能力.本文给出的方法可为该类型的移动机器人的设计和研究提供一种方法.后续可以进一步研究振动驱动型机器人在三维空间的建模方法、数值算法和动力学特性.