4. 音圈电机数学模型

音圈电机(Voice Coil Motor, VCM)本质上是一台无铁芯或准无铁芯的直线直流执行器,其电磁推力与电流近似成正比,机械运动近似为单自由度直线运动。要把它接入位置环、速度环或力控环,必须先建立从输入电压到动子位移的数学模型。本节围绕电压平衡方程、力平衡方程、电气时间常数、机械时间常数和传递函数推导展开,形成后续控制器设计的“被控对象基准”。

核心结论:VCM 的电气部分是一阶 RL 电路,机械部分是一阶质量-阻尼系统;两者通过“电流产生推力”和“运动产生反电动势”耦合,最终得到从电压到位移的三阶传递函数。

一、电压平衡方程:电气侧的一阶 RL 模型

音圈电机的线圈可等效为电阻 R 与电感 L 串联,动子运动时切割磁感线产生反电动势 e = Kev,其中 Ke 为反电动势常数,v 为动子速度。根据基尔霍夫电压定律,电枢回路电压平衡方程为:

u(t) = R·i(t) + L·di(t)/dt + K_e·v(t)

其中:
u(t)   —— 驱动器施加的线圈端电压(V)
i(t)   —— 线圈电流(A)
R      —— 线圈电阻(Ω)
L      —— 线圈电感(H)
K_e    —— 反电动势常数(V·s/m)
v(t)   —— 动子速度(m/s)

若忽略反电动势,电气回路退化为标准一阶 RL 电路,其电气时间常数为 Te = L / R。该常数决定电流对电压指令的响应快慢,是电流环带宽设计的关键参数。

二、力平衡方程:机械侧的质量-阻尼模型

线圈通电后在磁场中受到洛伦兹力 F = Kfi,其中 Kf 为推力常数(N/A)。动子受到推力、粘性阻尼力、弹性回复力和外部负载力共同作用。根据牛顿第二定律:

K_f·i(t) = m·dv(t)/dt + c·v(t) + k·x(t) + F_L(t)

其中:
m      —— 动子与负载等效质量(kg)
c      —— 粘性阻尼系数(N·s/m)
k      —— 弹性回复刚度(N/m,无弹簧时取 0)
x(t)   —— 动子位移(m)
F_L(t) —— 外部扰动力(N)

当系统无弹簧约束(k = 0)时,机械部分表现为“质量 + 阻尼”的一阶速度系统,其机械时间常数为 Tm = m / c。若存在弹簧或柔性连接,则机械部分变为二阶振荡环节,需额外关注谐振频率 ωn = √(k/m)。

三、电气时间常数与机械时间常数的物理意义

两个时间常数分别刻画电气与机械子系统的“惯性”。电气时间常数 Te 反映电流建立速度,典型 VCM 的 Te 在 0.1~2 ms 量级;机械时间常数 Tm 反映动子速度响应速度,典型值在 5~100 ms 量级。二者相差一到两个数量级,这为“电流环快、位置环慢”的级联控制提供了物理依据。

参数 符号 表达式 典型量级 工程含义
电气时间常数 Te L / R 0.1~2 ms 电流环响应速度
机械时间常数 Tm m / c 5~100 ms 速度环响应速度
反电动势常数 Ke 由磁路与匝数决定 1~30 V·s/m 速度反馈强度
推力常数 Kf 由磁路与匝数决定 1~30 N/A 电流到推力增益
机电时间常数 Tem mR / (KeKf) 1~50 ms 整体机电响应

四、传递函数推导:从电压到位移

对电压平衡方程和力平衡方程分别进行拉普拉斯变换,令初始条件为零:

U(s) = R·I(s) + L·s·I(s) + K_e·s·X(s)
K_f·I(s) = m·s²·X(s) + c·s·X(s) + k·X(s) + F_L(s)

由电气方程解出电流:

I(s) = [U(s) - K_e·s·X(s)] / (R + L·s)

代入机械方程,并令 FL(s) = 0,整理可得从电压到位移的传递函数:

G(s) = X(s) / U(s)
     = K_f / [ (R + L·s)(m·s² + c·s + k) + K_e·K_f·s ]

展开分母:

G(s) = K_f / [ L·m·s³ + (L·c + R·m)·s²
              + (R·c + L·k + K_e·K_f)·s + R·k ]

当 k = 0(无弹簧)时,分母退化为:

G(s) = K_f / [ L·m·s³ + (L·c + R·m)·s²
              + (R·c + K_e·K_f)·s ]

进一步忽略电感 L(电气时间常数远小于机械时间常数)时,可降阶为经典二阶模型:

G(s) ≈ K_f / [ R·m·s² + (R·c + K_e·K_f)·s ]
     = (K_f / R) / [ m·s² + (c + K_e·K_f / R)·s ]

该式表明:反电动势项 KeKf/R 等效为附加粘性阻尼,称为“电磁阻尼”。这是 VCM 天然具备速度阻尼特性的根本原因,也是电流环闭合后系统更稳定的物理来源。

音圈电机机电耦合模型与传递函数推导 电气侧:电压平衡方程 R L 反电动势 K_e·v u = R·i + L·di/dt + K_e·v 电气时间常数 T_e = L / R 机械侧:力平衡方程 m c k K_f·i = m·dv/dt + c·v + k·x + F_L 机械时间常数 T_m = m / c K_f·i K_e·v 传递函数推导:从电压到位移 ① 拉氏变换: U(s) = (R + L·s)·I(s) + K_e·s·X(s) ② 消去电流: I(s) = [U(s) - K_e·s·X(s)] / (R + L·s) ③ 代入机械方程: K_f·I(s) = (m·s² + c·s + k)·X(s) ④ 最终传递函数: G(s) = K_f / [ (R + L·s)(m·s² + c·s + k) + K_e·K_f·s ] 无弹簧、忽略电感时降阶为二阶:G(s) ≈ (K_f / R) / [ m·s² + (c + K_e·K_f / R)·s ]

五、案例:某 50 mm 行程 VCM 的参数辨识与模型验证

某精密贴装设备使用 50 mm 行程 VCM,实测参数为:R = 4.2 Ω,L = 0.8 mH,Ke = 8.5 V·s/m,Kf = 8.5 N/A,m = 0.35 kg,c = 12 N·s/m,无弹簧。计算得:

T_e = L / R = 0.8e-3 / 4.2 ≈ 0.19 ms
T_m = m / c = 0.35 / 12 ≈ 29.2 ms
电磁阻尼 = K_e·K_f / R = 8.5 × 8.5 / 4.2 ≈ 17.2 N·s/m
等效总阻尼 = c + K_e·K_f / R ≈ 12 + 17.2 = 29.2 N·s/m

可见电磁阻尼已超过机械阻尼,系统实际阻尼比远高于仅按机械参数估算的值。若忽略反电动势,会严重低估阻尼、误判系统为欠阻尼,导致位置环增益设置过大而振荡。

常见误区:把 VCM 当作纯“力发生器”建模,忽略反电动势 Kev。在低速大力矩场景下误差尚可接受,但在高速往返扫描或高带宽位置环中,反电动势会显著改变系统阻尼和相位,必须纳入模型。

六、复盘与工程要点

  • 电气时间常数决定电流环带宽上限,通常要求电流环带宽 ωci ≥ 5 / Te。
  • 机械时间常数决定速度环和位置环的响应节奏,位置环带宽一般不超过 1 / (5Tm)。
  • 反电动势提供天然电磁阻尼,是 VCM 比同尺寸直线电机更易稳定的重要原因。
  • 电感降阶仅在 Te ≪ Tm 时成立;若驱动器 PWM 频率低或线圈电感大,需保留三阶模型。
  • 弹簧刚度 k不可忽略时,系统出现谐振峰,需在控制器中加入陷波滤波器。

综上,VCM 的数学模型由电压平衡方程和力平衡方程耦合而成,电气时间常数与机械时间常数分别刻画两个子系统的响应速度,最终传递函数是理解电流环、速度环和位置环设计的基础。掌握该模型后,才能进一步开展控制器参数整定、带宽分配和抗扰设计。