2. 传动系统刚度建模基础

在精密传动系统中,刚度描述的是抵抗弹性变形的能力。若把传动链看作一条从驱动端到负载端的“弹性路径”,那么每一个齿轮啮合点、每一段轴、每一个轴承、每一圈丝杠滚道,都会在这条路径上贡献自己的柔度。理解刚度建模,首先要建立两个基本量的清晰定义:刚度 k 与柔度 c。刚度定义为产生单位变形所需的力或力矩,柔度则是刚度的倒数。在直线方向,k = F / δ,单位为 N/μm 或 N/m;在扭转方向,kθ = T / θ,单位为 N·m/rad。柔度 c = 1 / k,它更直观地表达了“在给定载荷下会产生多大变形”。工程上常把柔度作为串联计算的中间量,因为串联柔度可以直接相加,而串联刚度不能直接相加。

串联与并联是传动系统刚度合成中最基本的两种拓扑关系。串联指载荷依次经过多个弹性环节,例如“电机轴—联轴器—齿轮轴—齿轮啮合—输出轴”这条路径,各环节承受相同力或力矩,总变形为各环节变形之和,因此总柔度等于各环节柔度之和:c总 = c1 + c2 + … + cn,总刚度 k总 = 1 / (1/k1 + 1/k2 + … + 1/kn)。串联系统的总刚度永远小于其中最弱的环节刚度,这就是“短板效应”。并联指多个弹性环节共同承担同一载荷、变形相同,例如双电机驱动同一齿轮、双螺母预紧丝杠、多根立柱支撑同一工作台,此时总刚度等于各并联刚度之和:k总 = k1 + k2 + … + kn。并联能提高整体刚度,但前提是各支路变形协调、载荷分配合理。

串联与并联的判别口诀:载荷“依次经过”的是串联,柔度相加;载荷“同时分担”的是并联,刚度相加。实际传动链往往是串并联混合结构,需要先拆分子路径,再逐级合成。

传动系统中常见的刚度类型有三种:扭转刚度、弯曲刚度和轴向刚度。扭转刚度描述轴类零件抵抗扭转变形的能力,kθ = G·Ip / L,其中 G 为剪切模量,Ip 为极惯性矩,L 为轴段长度。它直接决定齿轮、丝杠、谐波减速器在传递扭矩时的角位移滞后。弯曲刚度描述轴或梁抵抗弯曲变形的能力,kb = 3E·I / L3(简支梁中点加载情形),其中 E 为弹性模量,I 为截面惯性矩。弯曲刚度不足会导致齿轮啮合错位、丝杠下垂、轴承偏载。轴向刚度描述零件沿轴线方向抵抗拉压变形的能力,ka = E·A / L,A 为截面积。丝杠的轴向刚度、轴承的轴向刚度、谐波减速器柔轮的轴向刚度都属于这一类。三种刚度并非孤立,实际变形往往是耦合的,但在初步建模中可先按主导方向分别计算,再通过坐标变换或等效折算统一到同一参考轴。

传动系统刚度建模:串联/并联合成与三类刚度 串联路径(柔度相加) 电机轴 联轴器 齿轮啮合 输出轴 c总 = c1 + c2 + c3 + c4 k总 = 1 / (1/k1 + 1/k2 + 1/k3 + 1/k4) 总刚度 ≤ 最弱环节刚度 并联路径(刚度相加) 支路 A 支路 B 共同负载 k总 = kA + kB 变形协调,载荷按刚度分配 三类基本刚度及其物理含义 扭转刚度 kθ kθ = G·Ip / L 抵抗扭矩引起的角变形 影响:齿轮/丝杠/谐波 角位移滞后与传动误差 单位:N·m/rad 弯曲刚度 kb kb = 3E·I / L³ 抵抗横向力引起的挠度 影响:啮合错位/丝杠下垂 轴承偏载与振动 单位:N/μm 轴向刚度 ka ka = E·A / L 抵抗轴向力引起的伸缩 影响:丝杠定位精度 轴承预紧与热变形 单位:N/μm

集中参数建模方法的核心思想是:把连续分布的弹性体离散为有限个“质量—弹簧—阻尼”单元,每个单元用集中刚度、集中质量和集中阻尼描述。对于齿轮传动,常把一对齿轮的啮合刚度等效为沿啮合线方向的时变弹簧 km(t),同时考虑齿轮轴的扭转刚度 kθ1、kθ2 和轴承的径向支撑刚度 kb。对于丝杠进给系统,常把丝杠离散为若干段,每段用轴向刚度 ka,i 和扭转刚度 kθ,i 表示,螺母副用轴向接触刚度 knut 表示,两端轴承用轴向和径向刚度表示。对于谐波减速器,柔轮是核心弹性环节,其扭转刚度 kflex 通常呈非线性,需用分段线性或多项式拟合;波发生器轴承刚度和刚轮固定刚度也需纳入模型。

下面给出一个集中参数建模的简化示例,用于说明参数如何组织与合成。假设某单级齿轮传动系统,电机轴扭转刚度 kθ1 = 2.0×104 N·m/rad,联轴器扭转刚度 kθ2 = 1.5×104 N·m/rad,齿轮啮合等效扭转刚度(折算到电机轴)kθ3 = 8.0×103 N·m/rad,输出轴扭转刚度 kθ4 = 3.0×104 N·m/rad。由于这些环节沿扭矩传递路径依次串联,总柔度为:

c总 = 1/20000 + 1/15000 + 1/8000 + 1/30000
    = 5.00e-5 + 6.67e-5 + 1.25e-4 + 3.33e-5
    = 2.75e-4 rad/(N·m)

k总 = 1 / c总 ≈ 3636 N·m/rad

可以看到,齿轮啮合环节的柔度最大(1.25×10-4),是系统的薄弱环节,总刚度约 3636 N·m/rad,远低于任一单独环节的刚度。若把齿轮啮合刚度提高到 1.6×104 N·m/rad,则总柔度降为 2.125×10-4,总刚度提升到约 4706 N·m/rad,提升约 29%。这说明在串联系统中,改善最弱环节对总刚度提升最有效。

建模提示:集中参数建模时,先画传动链拓扑图,标出每个弹性环节的刚度类型和方向;再判断串并联关系;最后统一折算到同一参考轴(通常选电机轴或负载轴)。折算时注意传动比平方关系:从动轴刚度折算到主动轴时乘以传动比平方,反之除以传动比平方。

在齿轮、丝杠、谐波减速器中的初步应用可以这样理解:齿轮系统重点关注啮合刚度的时变特性和轴系扭转刚度,啮合刚度随啮合位置周期变化,是齿轮振动的主要激励源;丝杠系统重点关注轴向刚度沿行程的变化,因为丝杠一端固定一端支撑时,不同位置的轴向刚度不同,导致定位精度随位置波动;谐波减速器重点关注柔轮扭转刚度的非线性,其刚度随负载扭矩增大而降低,在小扭矩区尤为明显,直接影响机器人关节的动态响应和定位精度。把这三类对象的刚度模型统一到“集中参数—串并联合成—折算到参考轴”的框架下,就完成了传动系统刚度建模的第一步。

实际建模时还需注意:刚度参数应尽量通过实验辨识或供应商数据获取,避免仅凭几何公式估算;对于预紧轴承、预压丝杠螺母,其刚度随预紧力变化,需在模型中体现;温度变化会引起材料模量和预紧状态改变,高精度场合应引入热—刚度耦合修正。完成刚度建模后,下一步才能进行固有频率计算、动态响应分析和振动控制策略设计。