摩擦模型及其磁滞特性的研究

摩擦模型及其磁滞特性的研究1

张晓巧，王兴松

东南大学机械工程学院，江苏南京 （210096）

E-mail: linlangqiao820926@

摘 要：摩擦是存在于两接触面间的一种非常复杂的现象，直接影响运动控制系统的定位精度和速度。根据运动控制发展的需要，本文提出将机械伺服系统中出现的摩擦现象与磁滞现象综合考虑的理念。首先介绍了摩擦非线性环节典型的静力学模型和动力学模型，以及各种模型的特点及其对摩擦特征的描述能力，重点剖析了用偏微分方程描述的动力学模型，在速度反向时摩擦力的滞后现象，以及低速时的Stribeck 现象。最后在模型分析的基础上，通过仿真和实验研究，验证了摩擦具有磁滞特性。研究工作为摩擦、磁滞的综合建模和统一补偿控制，提供了理论基础。

关键词：摩擦，静力学模型，动力学模型，摩擦力滞后，磁滞特性

1. 引言

摩擦存在于所有具有相对运动或相对运动趋势的两接触面之间，是一种复杂的、非线性的、具有不确定性的自然现象。有时我们需要利用摩擦，如刹车、离合器等，但更多的时候我们要减小摩擦影响。在机械伺服运动系统中，由摩擦非线性引起的运动特性将直接影响加工定位的精度和速度。

随着先进装备运动精度要求的提高和数字控制的广泛应用，对摩擦补偿的研究也更具有其经济实用性。许多学者根据摩擦现象的动力学特性提出了相关的数学模型，并在此基础上运用控制算法实现摩擦的补偿，取得一定成效。

但在多数实际机械系统中，摩擦和其它非线性现象往往同时存在、同时作用（如间隙、饱和、磁滞等），比如用压电微驱动系统进行高精度误差补偿时，会引入磁滞问题。如果对它们进行分别建模、补偿，由于方法不统一，实际运用困难。基于摩擦本身就具有滞后效应（friction lag ）的特性，作者提出对摩擦、磁滞进行统一建模补偿的理念，以提高精密伺服运动系统的精度。本文主要讨论摩擦模型的描述特点及其磁滞特性。

2. 摩擦数学模型分析

建立准确的摩擦非线性数学模型，对深入认识摩擦现象、进行精确补偿具有重要的意义。因此，有关摩擦建模的研究一直是该领域学者关注的热点。到目前为止，提出的摩擦模型共有30余种，它们为摩擦非线性的补偿研究奠定了

基础。

2.1 摩擦静力学模型

2.1.1 古典摩擦模型

古典摩擦模型[1]首次由达芬奇（1519）提出，

该模型认为摩擦力只与法向负载成比例关系，与

运动方向相反，而与运动速度和接触面无关。后

由库仑(1785)等科学家将其发展完善。随后莫因

1本课题得到国家自然科学基金（50475076）资助

(Morin,1833)引入静摩擦、雷诺(Reynolds,1866)引入粘性流体摩擦，继而形成了在工程中广泛应用的摩擦模型，即库仑摩擦+粘滞摩擦+静摩擦。如图1(c)所示。

库仑模型，又可称sgn 函数模型，即

sgn()c F F v = (1)

该模型完美诠释了达芬奇提出的摩擦假设，并以其简洁的表达式曾一度被广泛用于补偿技术中。但是在运动速度过零点处的摩擦力具有不确定性，可在[],c c F F −区间取任意值。

19世纪流体动力学的发展引出力与润滑粘度相关的描述。粘性摩擦（viscous friction ）模型的提出首次将速度变量引入到摩擦力表达式中，即

sgn()v v F F v v δ= (2)

粘性模型通常与库仑模型结合形成如图1(b)所示的摩擦特性。其中v δ是几何学参数，当1v δ=时，摩擦力与速度在单向区间成线性关系。

静摩擦（stiction ）描述了物体在静止状态下的摩擦特性，通常比库仑力大。很显然在没有速度的状态下，静摩擦力不能由速度表达式来描述，而是外界力的函数，即

0&sgn()0&e e s s e e s

F if v F F F F F if v F F =<⎧=⎨=≥⎩ (3) 由以上三种古典模型结合形成的摩擦模型如图1(c)所示，它由常数和速度函数组成，描述的摩擦力是不连续的。然而在实际应用中，速度和摩擦力都是连续变化的，如图1(d)所示。所以在以往模型的基础上，Stribeck 曲线模型的提出，具有更普遍的意义，即

()00&sgn()e s

e F v i

f v e s F F if v F F F otherwise ⋅≠⎧⎪=⋅=⎨⎪⎩F < (4) 其中()F v 是速度相关曲线函数，用来描述摩擦力与速度的关系特性。目前已提出多种参数化函数形式[2]，比较常用的表达式为

/()()s s v v c s c v F v F F F e −=+−+F v (5)

其中s v 为Stribeck 速度，如图2所示，即摩擦力负斜率段的初始

切线与库仑力的交点值，该值也直接影响低速段摩擦特性曲线

的陡度。选取指数衰减函数作为静摩擦力到粘滞摩擦的过渡，

能够很好地描述摩擦力的连续性。通过实验测得实际的摩擦力

曲线通常不对称，这主要因为反向时接触面的润滑性质不同。

古典摩擦模型虽然有模型简单，参数容易确定的优点，但它只描述了摩擦与速度的静态关系，不能解释摩擦的磁滞效应，也不能描述随实验条件变化临

界摩擦力变化的现象，以及静摩擦状态下接触面间小的预偏移

(Dahl 效应)现象。这也是这类模型在高精度、低速跟踪控制时，

不能达到理想效果的原因。

2.1.2 Karnopp 摩擦模型

上述模型在不同层次上描述了摩擦力随速度变化的特性，但

都存在一个基本问题，即速度零点的摩擦力描述不确定性。由于

摩擦力在速度零点的不连续性使得其很难检测。

这对于通过由位置

信号求导估算速度的系统，更为复杂。由于测量噪声的存在，离散位置值计算出的速度，可能根本没有零值。为解决这个问题，Karnopp 提出的模型对其进行了修补[3]，巧妙地克服了零速度检测的困难，并消除了静止和滑移时不同状态方程之间的切换问题。如图3所示，该模型通过在较小的速度附近定义一个微小的“粘性”区域士Dv ，当v Dv <时认定（速度死区），摩擦力与作用在系统上的其它合力平衡，直到超过最大静摩擦力时，进入滑动阶段。

0v =karnopp 模型是介于静态模型和动态模型之间的一种简化模型，通过规定静止区间，降低了对速度信号质量的要求，更适合于工程实践。但它与系统其它受力有很大的关系，需要根据不同的系统结构，来剪裁合适的模型。

2.1.3 Armstrong 摩擦模型

Armstrong 模型[1]在修改了古典模型的基础上，考虑了随时间变化的静态阻力和Stribeck 效应，及预滑移运动。模型由三部分组成，通过机械开关或逻辑状态量进行切换，并在切换时初始化参数。

2,,,()1(,)((,)

)sgn()1(()/)(,)()f t c s d v l d s d s a s s a d

x k x v t F F t v F v v t vs t t F F F t γτγγ∞=−=+++−=+−+ （6） 表达式共包含7个参数，所以又称为七参数模型。为速度过零时的驻留时间。(6)式中的滑移摩擦力方程采用速度的延时形式d t (l v t )τ−来替代Stribeck 摩擦力相对应的瞬时速度，以静力学方程的形式体现其动力学特性。

Armstrong 模型虽然能够反映摩擦的动、静态特性，但其实质是将静态模型与动态特性生硬地组合，缺乏明确的物理意义，且存在冗余参数。Armstrong 还应用量纲和摄动分析方法建立了具有5个参数的模型。

2.2 摩擦动力学模型

以上介绍的典型摩擦静力学模型，已经基本描述了摩擦在静力作用下的表现特性。但是随着精密伺服系统和高端硬件技术的发展，基于静摩擦模型的控制精度已达极限，人们对复杂摩擦动力学模型的研究兴趣日趋浓厚。

2.2.1 Dahl 摩擦模型

Dahl 以滚珠轴承伺服控制实验为出发点[4]，发现了类似于固体摩擦的摩擦特性，并由此提出了简单的摩擦数学模型。Dahl 模型是重要的动态模型，在其基础上发展的LuGre 模型和Bliman-Sorine 模型也代表了这一领域的最新成果。

Dahl 从应力-应变曲线出发，用偏微分方程将应力作用过程描述出来，即

(1sgn )c

dF F v dx F

ασ=

− （7）

其中x 代表变形量，F 即摩擦力，c F 为库仑摩擦力，σ

是刚度系数，α参数决定了曲线的形状。当1α=时，

可以得到如图4所示的曲线。曲线在库仑力范围内变

化，当速度反向时，由于sgn 符号函数的作用使得微增量