4. 音圈电机数学模型

音圈电机（Voice Coil Motor, VCM）本质上是一台无铁芯或准无铁芯的直线永磁同步执行器，其电磁与机械过程耦合紧密。要设计电流环、位置环乃至前馈补偿，必须先把物理对象抽象成可计算的数学模型。本节从电气侧、电磁力侧、机械侧三条主线出发，建立电压方程、力方程与运动方程，并进一步整理成状态空间形式，为后续控制器设计提供统一基础。

一、三条基本方程：从电路到运动

音圈电机的动子线圈处于永磁体产生的气隙磁场中。当线圈通电时，受到安培力驱动；同时线圈在磁场中运动又会产生反电动势。因此，电气与机械两个子系统通过“电流—力”和“速度—反电动势”两条通道双向耦合。

1. 电压方程（电气侧）

设线圈端电压为 u(t)，回路电流为 i(t)，线圈电阻为 R，电感为 L，反电动势系数为 Ke，动子速度为 v(t)，则电压平衡关系为：

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

其中 Ke·v(t) 是反电动势项，它随速度增大而增大，直接削弱有效驱动电压。这是音圈电机高速运行时电流难以维持的根源。

2. 力方程（电磁侧）

设电磁推力为 Fm(t)，推力常数（力常数）为 Kf，则：

Fm(t) = Kf·i(t)

在理想线性区内，Kf 近似为常数，单位常用 N/A。实际中由于磁路边缘效应、行程端部磁场不均，Kf 会随位置轻微波动，高精度场合需做位置相关修正。

3. 运动方程（机械侧）

设动子及负载总质量为 m，粘滞阻尼系数为 c，弹簧或柔性支撑刚度为 k，外部扰动力为 Fd(t)，位移为 x(t)，则：

m·d²x(t)/dt² = Kf·i(t) − c·dx(t)/dt − k·x(t) − Fd(t)

若为无弹簧直驱结构，可令 k = 0；若存在导轨摩擦、线缆拖拽力，则统一归入 Fd(t) 或等效阻尼。

二、从微分方程到状态空间模型

为了便于用现代控制方法分析，选取状态变量：电流 i、位移 x、速度 v，即状态向量 X = [i, x, v]ᵀ。由前述三式可整理为：

di/dt = (1/L)·u − (R/L)·i − (Ke/L)·v
dx/dt = v
dv/dt = (Kf/m)·i − (c/m)·v − (k/m)·x − (1/m)·Fd

写成标准状态空间形式 Ẋ = A·X + B·u + E·Fd，输出方程可取 y = C·X：

A = [ −R/L    0      −Ke/L
        0      0        1
       Kf/m  −k/m    −c/m  ]

B = [ 1/L
       0
       0 ]

E = [  0
       0
     −1/m ]

C = [ 0  1  0 ]   （以位移为输出）

该模型清晰揭示了三个关键极点来源：电气时间常数 L/R、机械谐振频率 √(k/m) 以及阻尼比 c/(2√(km))。当 k=0 时，机械侧退化为一个积分环节加一阶惯性，系统呈现“双积分 + 电气惯性”的特征。

三、参数关系与典型取值

下表汇总了音圈电机建模中常用参数及其对系统行为的影响，便于在辨识与调试时快速定位敏感项。

参数	符号	典型单位	主要影响
线圈电阻	R	Ω	决定稳态电流与铜耗，影响电气时间常数
线圈电感	L	mH	限制电流变化率，影响电流环带宽
反电动势系数	Ke	V·s/m	高速时削弱有效电压，限制最高速度
推力常数	Kf	N/A	决定同等电流下的输出推力
动子质量	m	kg	影响加速度能力与机械谐振频率
阻尼系数	c	N·s/m	影响超调与稳定时间
支撑刚度	k	N/m	决定机械谐振点，柔性支撑需重点关注

四、典型场景：位置阶跃响应中的模型表现

假设某音圈电机参数为 R=8 Ω、L=2 mH、Ke=12 V·s/m、Kf=12 N/A、m=0.35 kg、c=6 N·s/m、k=0。当给定电压阶跃时，电流先快速上升，随后因动子加速产生反电动势而回落；位移则经历“加速—匀速—减速”过程。若忽略反电动势，仿真会明显高估高速段推力，导致位置环增益设计偏大，实际出现振荡或过冲。

五、常见误区与复盘

• 误区一：把 Kf 与 Ke 当作无关参数。在理想单位制下二者数值近似相等（SI 单位下 Kf = Ke），辨识时应交叉验证，避免量纲混乱。
• 误区二：忽略电感动态。当电流环带宽接近 R/L 时，若仍按纯电阻负载整定 PI 参数，电流会出现明显滞后与超调。
• 误区三：把线缆拖拽力当作恒定扰动。线缆力随位置变化，若只用常值前馈补偿，行程两端仍会有残差。
• 误区四：状态空间建模时遗漏输出矩阵。位置控制需要 x 作为输出，速度前馈需要 v，输出矩阵选择错误会导致观测器设计失效。

掌握上述三方程与状态空间表达后，电流环可基于电气一阶模型整定，位置环可基于机械二阶模型设计，前馈项则可直接从 Ke·v 与 m·a 推导。模型越贴近真实机电耦合关系，后续控制器的鲁棒性与定位精度就越有保障。