应用数学进展  >> Vol. 10 No. 4 (April 2021)

平面库伦摩擦系统的动力学分析
Dynamic Analysis of Planar Coulomb Friction System

DOI: 10.12677/AAM.2021.104143, PDF, HTML, XML, 下载: 8  浏览: 23 

作者: 陈会敏, 冷 轩:长沙理工大学数学与统计学院,湖南 长沙

关键词: 库伦摩擦系统滑模域Poincaré映射全局吸引子收敛Coulomb Friction System Sliding Mode Domain Poincaré Mapping Global Attractor Convergence

摘要: 机械系统库伦摩擦的存在,导致了系统的不连续性的存在。本文研究了在水平面上运动的物体在库伦摩擦力作用下的全局动力学行为。在这种情况下,不连续系统的滑模域上的点都是系统的伪平衡点。在此基础上,通过构造Poincaré映射,给出了滑模域中的点构成的集合均为系统的全局吸引子这一结论,最后通过数值模拟验证了该吸引子在有限时间内是收敛的这一结论。
Abstract: The existence of Coulomb friction in mechanical system leads to the existence of discontinuity. In this paper, the global dynamic behavior of an object moving on a horizontal plane under Coulomb friction is studied. In this case, the points on the sliding mode domain of the discontinuous system are all pseudo-equilibrium points of the system. On this basis, by constructing the Poincaré mapping, the conclusion that the set of points in the sliding mode domain is the global attractor of the system is given. Finally, the conclusion that the attractor is convergent in finite time is proved by numerical simulation.

文章引用: 陈会敏, 冷轩. 平面库伦摩擦系统的动力学分析[J]. 应用数学进展, 2021, 10(4): 1336-1341. https://doi.org/10.12677/AAM.2021.104143

1. 引言

摩擦是表面接触时产生的一种非常复杂的现象,摩擦现象无处不在,其存在于两个物体之间,并不是单一存在的 [1] [2]。而库伦摩擦模型是较早存在的对摩擦现象进行描述的一种模型。许多学者对其进行了研究并取得了一些实质性的进展。Edward Berger [3] 综述许多不同工程领域的大量文献,并批判性地考察了不同情况下各种摩擦模型的使用,得到系统模型和摩擦模型从根本上是耦合的,不能独立选择,摩擦模型的有用性和系统动力学模型的成功是相互依赖的,摩擦性质和表面相互作用科学的基本问题将是未来许多领域的重要研究方向。Andrea Bacciotti [4] 利用局部Lipschitz条件和Lyapunov函数研究了Filippov意义下的具有右端不连续系统的稳定性和正定性。得到了一般微分包含情况下的稳定性,并给出了Jurdjevic-Quinn型反馈律稳定的一些充分条件。Leine R.I..和Campen D. [5] 用广义雅可比矩阵分析了两个平面非光滑连续系统的平衡分支,研究了一个属于微分包含类的非自治Filippov系统的力学例子,给出了周期解的一些显著的不连续分支,为后续相关研究提供了参考思路。Alvarez J. [6] 基于右端不连续常微分方程的n不变性原理,研究了具有库伦摩擦的一自由度机械振子的渐近稳定性,并进行了仿真验证,实现了控制器和观测器的有限时间收敛。以往关于摩擦动力学的研究大多忽视其实际的动力学模型,将其近似线性或弱化线性系统简化为线性系统来分析研究 [7] [8]。而本文在前人研究的基础上,重点讨论了在水平面上运动的物体所受的库伦摩擦力模型,并进一步研究系统运动的动力学行为。

2. 系统数学模型的建立

考虑形式较简单的库伦摩擦系统,如图1所示:

Figure 1. Coulomb friction system

图1. 库伦摩擦系统

假设一个质量为m的物体在水平面上运动,物体只受到弹簧的弹力和地面对它的摩擦力,在运动过程中物体所受的摩擦力的大小是恒定的,其方向与运动方向相反,则根据牛顿第一运动定律,系统的动力学行为可描述为

m x ¨ + p ( x ˙ ) + k x = 0 (1.1)

其中:x表示物体的位移,k为弹簧的弹力系数,且

p ( x ˙ ) = { p 0 , x ˙ > 0 p 0 , x ˙ < 0 (1.2)

由此可知,系统(1.1)~(1.2)为不连续系统。

p 0 0 时,由于系统(1.1)~(1.2)是一个右端不连续的微分系统,当考虑到系统(1.1)在Filippov意义下的解,不难知道系统(1.1)的解是存在的,且满足初值条件 x ( 0 ) = x 0 , y ( 0 ) = y 0 , 存在区间为 [ 0 , + ) ,并且系统的解是右唯一的。为了后续的方便讨论,将系统(1.1)~(1.2)化为等价的平面微分系统,令 x ˙ = y

{ x ˙ = y y ˙ = k x m p ( y ) m (1.3)

其中

p ( y ) = { p 0 , y > 0 p 0 , y < 0 (1.4)

由(1.3)~(1.4)可得系统(1.3)的上半线性子系统为

{ x ˙ = y y ˙ = k x m p 0 m , y > 0 (1.5)

下半线性子系统为

{ x ˙ = y y ˙ = k x m + p 0 m , y < 0 (1.6)

对于上半子系统 e U = ( p 0 k , 0 ) 为其中心,且通解为

m y 2 + k ( x + p 0 k ) 2 = c 1 (1.7)

其中 c 1 为非负常数。

同理,下半子系统以 e L = ( p 0 k , 0 ) 为中心,其通解为

m y 2 + k ( x p 0 k ) 2 = c 2 (1.8)

其中 c 2 为非负常数。

3. 切换线上滑模动力学行为

由(1.3)式可知,系统的切换线是 Σ = { x 2 : y = 0 } ,紧接着着重分析系统(1.3)在切换线 Σ = { x 2 : y = 0 } 上和它周围附近的动力学行为。则系统(1.3)在 y = 0 附近的向量场分析如图2所示:

Figure 2. Vector field distribution

图2. 向量场分布

一般地,当弹簧被拉直时,物体将不会再运动,不妨设弹簧被拉直时所处的位置范围是 { x | A 1 < x < A 2 } ,为了后面的方便讨论,给出了 A 1 > 3 p 0 k A 2 > 3 p 0 k ,但 A 1 A 2 不一定相等,则 Σ 可分为三个部分

Σ C 1 = { ( x , y ) | A 1 < x < p 0 k , y = 0 }

Σ s = { ( x , y ) | p 0 k < x < p 0 k , y = 0 }

Σ C 2 = { ( x , y ) | p 0 k < x < A 2 , y = 0 }

其中 Σ C 1 Σ C 2 为穿越集,且在 Σ C 1 上系统(1.3)的解从下往上穿越,而在 Σ C 2 上系统(1.3)的解从上往下穿越, Σ s 为滑模域,并且该滑模域是吸引的。为了分析滑模域上的动力学行为,其滑模方程起着很重要的作用,根据Filippov凸方法可以可得,在滑模域 Σ s 上的滑模方程为

{ x ˙ = y , y ˙ = 0 , (2.1)

由此可知, Σ s 上的所有点均为系统(1.3)的伪平衡点,即 Σ s 为系统(1.3)的伪平衡点集,显然 Σ s 中的伪平衡点均不孤立。

4. 系统的全局动力学行为

下面将用Poincaré映射来分析系统(1.3)的全局动力学行为。

注意到从 ( x 0 , 0 ) Σ C 2 出发的解将进入 S 2 ,并经过有限的时间再次与 Σ 相交不妨设交点为 ( x 1 L , 0 ) ,且 x 1 L < p 0 k ,则下半Poincaré映射 P L ( ) 定义如下:

P L ( x 0 ) = { x 1 L , x 0 ( p 0 k , A 2 ] x 0 , x 0 = [ A 1 , p 0 k ]

另一方面,若 ( x 0 , 0 ) Σ C 1 ,即 x 0 < p 0 k 出发的解将进入 S 1 ,并经过有限时间再次与 Σ 相交。不妨设交点为 ( x 1 U , 0 ) ,且 x 1 U > p 0 k 。同理,我们可定义上半Poincaré映射 P U ( ) 如下所示:

P U ( x 0 ) = { x 1 U , x 0 [ A 1 , p 0 k ) x 0 , x 0 = [ p 0 k , A 2 ]

由上述通解公式可知,上半子系统的通解为椭圆,椭圆的中心为 e U = ( p 0 k , 0 ) ,从而

x 1 U = P U ( x 0 ) = x 0 2 p 0 k , x 0 [ A 1 , p 0 k ) (3.1)

同理可得

x 1 L = P L ( x 0 ) = x 0 + 2 p 0 k , x 0 ( p 0 k , A 2 ] (3.2)

进一步,根据(3.1)和(3.2)可定义Poincaré映射如下

P ( x 0 ) = P U ( P L ( x 0 ) ) = x 0 4 p 0 k , x 0 [ 3 p 0 k , A 2 ] (3.3)

定理3.1 由滑模域中的点构成的集合 Σ s = { ( x , 0 ) | | x | p 0 k } 为系统(1.3)的全局吸引子,且该吸引子在有限时间内是收敛的。

证明:由子系统的通解公式(1.7)和(1.8)可知,系统(1.3)的任意一个解都将在有限时间内到达切换线 Σ 上,下面只需证明从切换线 Σ 出发的解在有限时间内收敛到平衡点集 Σ S 上,从而只需证明对于任意 x 0 Σ S ,存在正整数 N 使得

P ( x 0 ) Σ S , P N = P P P N

事实上,根据Poincaré映射的表达式(3.3),当 x 0 [ 3 p 0 k , 3 p 0 k ] 时,则 P ( x 0 ) Σ S ,即此时取 N = 1 ,当 x 0 ( 3 p 0 k , A 2 ] 时,由于 | P ( x 0 ) x 0 | = 4 p 0 k ,因此取 N = [ k x 0 p 0 4 p 0 ] + 1 ,则 P N ( x 0 ) S ;而当 x 0 [ A 1 , 3 p 0 k ) 时, P ( x 0 ) ( 3 p 0 k , A 2 ] ,从而取 N = [ k x 0 3 p 0 4 p 0 ] + 2 可得 P N ( x 0 ) S 。证毕!

5. 数值仿真及应用分析

k = 1 , m = 1 , p 0 = 1 2 时,系统的滑模域区间为 [ 0.5 , 0.5 ] ,从图3中可以看出,在有限时间内从任何点出发的轨线,最后都收敛到 [ 0.5 , 0.5 ] 区间的滑模域上。

Figure 3. Trajectory phase diagram of plane Coulomb friction system

图3. 平面库伦摩擦系统的轨线相图

6. 结论

线性运动理论是运动力学研究中应用最为广泛的工具。以往研究人员都将近似线性或弱化线性系统简化为线性系统来分析研究,忽视所研究的实际动力学系统模型。本章首先建立了在水平面上带有库伦摩擦系统的动力学模型,进而研究了其运动行为,然后讨论了具有两个子系统且子系统的切换线为直线的分段线性微分系统的动力学行为,并通过构造Poincaré映射得到了滑模域中的点构成的集合为系统的全局吸引子,且该吸引子在有限时间内是收敛的这一结论。最后通过数值模拟,证实了这一结论的正确性。

参考文献

[1] Urbakh, M., Klafter, J., Gourdon, D., et al. (2004) The Nonlinear Nature of Friction. Nature, 430, 525.
https://doi.org/10.1038/nature02750
[2] 王松年. 摩擦学原理及应用[M]. 北京: 中国铁道出版社, 1990.
[3] Berger, E. (2002) Friction Modeling for Dynamic System Simulation. Applied Mechanics Reviews, 55, 25-32.
https://doi.org/10.1115/1.1501080
[4] Bacciotti, A. and Ceragioli. F. (1999) Stability and Stabilization of Discontinuous Systems and Nonsmooth Lyapunov Functions. ESAIM Control Optimisation and Calculus of Variations, 4, 361-376.
https://doi.org/10.1051/cocv:1999113
[5] Leine, R.I. and Campen, D. (2006) Bifurcation Phenomena in Non-Smooth Dynamical Systems. European Journal of Mechanics—A/Solids, 25, 595-616.
https://doi.org/10.1016/j.euromechsol.2006.04.004
[6] Alvarez, J., Orlov, I. and Acho, L. (2000) An Invariance Principle for Discontinuous Dynamic Systems with Application to a Coulomb Friction Oscillator. Journal of Dynamic Systems Measurement and Control, 122, 317-322.
https://doi.org/10.1115/1.1317229
[7] 黄立宏, 郭振远, 王佳伏. 右端不连续微分方程理论与应用[M]. 北京: 科学出版社, 2011.
[8] 丁千, 翟红梅. 机械系统摩擦动力学研究进展[J]. 力学进展, 2013, 43(1): 112-131.