4. 音圈电机数学模型
音圈电机(Voice Coil Motor, VCM)本质上是一台无铁芯或准无铁芯的直线直流执行器,其电磁推力与电流近似成正比,机械运动近似为单自由度直线运动。要把它接入位置环、速度环或力控环,必须先建立从输入电压到动子位移的数学模型。本节围绕电压平衡方程、力平衡方程、电气时间常数、机械时间常数和传递函数推导展开,形成后续控制器设计的“被控对象基准”。
一、电压平衡方程:电气侧的一阶 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 天然具备速度阻尼特性的根本原因,也是电流环闭合后系统更稳定的物理来源。
五、案例:某 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
可见电磁阻尼已超过机械阻尼,系统实际阻尼比远高于仅按机械参数估算的值。若忽略反电动势,会严重低估阻尼、误判系统为欠阻尼,导致位置环增益设置过大而振荡。
Kev。在低速大力矩场景下误差尚可接受,但在高速往返扫描或高带宽位置环中,反电动势会显著改变系统阻尼和相位,必须纳入模型。
六、复盘与工程要点
- 电气时间常数决定电流环带宽上限,通常要求电流环带宽
ωci ≥ 5 / Te。 - 机械时间常数决定速度环和位置环的响应节奏,位置环带宽一般不超过
1 / (5Tm)。 - 反电动势提供天然电磁阻尼,是 VCM 比同尺寸直线电机更易稳定的重要原因。
- 电感降阶仅在
Te ≪ Tm时成立;若驱动器 PWM 频率低或线圈电感大,需保留三阶模型。 - 弹簧刚度 k不可忽略时,系统出现谐振峰,需在控制器中加入陷波滤波器。
综上,VCM 的数学模型由电压平衡方程和力平衡方程耦合而成,电气时间常数与机械时间常数分别刻画两个子系统的响应速度,最终传递函数是理解电流环、速度环和位置环设计的基础。掌握该模型后,才能进一步开展控制器参数整定、带宽分配和抗扰设计。