2. 扭振数学模型：多惯量系统建模、扭转刚度与阻尼、运动方程推导
各位工程师朋友，咱们今天聊聊扭振的数学模型。说实话，很多人一听到「数学模型」四个字就头大，觉得那是搞理论的人干的事。但我在现场摸爬滚打这么多年，越来越觉得——不懂模型，你连故障在哪儿都说不清楚。
扭振分析，说白了就是研究「转轴怎么拧着抖」。你想想看，一根轴连着电机和负载，电机一发力，轴就扭。扭得均匀还好，要是扭得忽快忽慢，那就是扭振。咱们今天就把这个事儿用数学说清楚。
核心思路：把复杂的传动链拆成几个「惯性块」，中间用「弹簧」和「阻尼器」连着。这就是多惯量系统建模的基本思想。
2.1 多惯量系统建模：把传动链拆开看
我在做轧机传动系统诊断时，遇到过一台主传动电机，振动频谱上全是扭振的边频。当时我就在想，这轴系到底是怎么抖起来的？后来我把它拆成三个部分：电机转子、联轴器、轧辊。每个部分都有自己的转动惯量，中间用扭转刚度和阻尼连接。
为什么要拆成多惯量？因为实际传动链不是一根刚体。你想想看：
- 电机转子：质量大，惯性大，像个飞轮
- 齿轮箱：内部有多级齿轮，每级都有自己的惯量和刚度
- 联轴器：弹性元件，主要贡献刚度和阻尼
- 负载：比如轧辊、风机叶片，惯量往往最大
每个部件都可以看作一个「惯量块」。惯量块之间通过轴段连接，轴段提供扭转刚度和结构阻尼。这就是多惯量系统的物理基础。
我的经验：实际建模时，不是惯量块越多越好。我曾经把一个风机传动链拆成7个惯量块，结果模型复杂到算不动，而且参数也测不准。后来简化成3个惯量块，反而抓住了主要问题。记住：模型是工具，不是目的。
2.2 扭转刚度与阻尼：轴是怎么「软」的
扭转刚度，说白了就是轴抵抗扭转变形的能力。你拿一根筷子，两头一拧，筷子就扭了。拧到一定程度，「啪」断了。轴也是一样，只不过它不会断，但会扭。
扭转刚度的计算公式很简单：
K = G * J / L其中：
- G：剪切模量，材料属性（钢大约是80 GPa）
- J：极惯性矩，跟轴径的四次方成正比（轴粗一点，刚度大很多）
- L：轴段长度（轴越长，越软）
这里有个坑，我踩过。有一次我算一根长轴的扭转刚度，用了名义直径。结果现场测试发现，实际刚度比计算值低了30%。后来一查，轴上有键槽和退刀槽，这些局部削弱会显著降低刚度。所以，计算刚度时一定要考虑应力集中和截面变化。
再说阻尼。扭转阻尼比刚度难搞多了。它来自：
- 材料内阻尼：轴材料本身的耗能特性
- 结构阻尼：联轴器、花键、键连接处的摩擦
- 外部阻尼：比如齿轮啮合时的油膜阻尼
阻尼系数通常用 C 表示，单位是 N·m·s/rad。实际工程中，阻尼很难精确计算，我一般用经验公式或者模态测试反推。
注意：阻尼对扭振的影响非常大。阻尼太小，共振时振幅会无限大（理论上）。阻尼足够大，共振就被抑制了。我见过一个案例，风机传动链在启动过程中扭振剧烈，后来在联轴器上加了粘滞阻尼器，问题就解决了。所以，别小看阻尼。
2.3 运动方程推导：从物理到数学
好，现在咱们把物理模型变成数学方程。以最简单的两惯量系统为例：
一个电机（惯量 J₁）通过一根轴（刚度 K，阻尼 C）驱动负载（惯量 J₂）。
对每个惯量块，应用牛顿第二定律（旋转形式）：
J₁ * α₁ = T₁ - T_shaft
J₂ * α₂ = T_shaft - T_load其中：
- α₁、α₂：角加速度
- T₁：电机电磁转矩
- T_load：负载转矩
- T_shaft：轴传递的转矩
轴转矩由两部分组成：弹性恢复力（刚度项）和阻尼力（阻尼项）：
T_shaft = K * (θ₁ - θ₂) + C * (ω₁ - ω₂)其中 θ₁、θ₂ 是角位移，ω₁、ω₂ 是角速度。
把上面三个方程联立，就得到两惯量系统的运动方程：
J₁ * d²θ₁/dt² + C * (dθ₁/dt - dθ₂/dt) + K * (θ₁ - θ₂) = T₁
J₂ * d²θ₂/dt² + C * (dθ₂/dt - dθ₁/dt) + K * (θ₂ - θ₁) = -T_load写成矩阵形式更清爽：
[J₁ 0 ] [d²θ₁/dt²] [ C -C] [dθ₁/dt] [ K -K] [θ₁] [ T₁ ]
[0 J₂] [d²θ₂/dt²] + [-C C] [dθ₂/dt] + [-K K] [θ₂] = [-T_load]这就是两惯量扭振系统的标准数学模型。你看着可能觉得复杂，但它的物理意义很清楚：
- 左边第一项：惯性力（跟加速度有关）
- 左边第二项：阻尼力（跟速度差有关）
- 左边第三项：弹性力（跟位移差有关）
- 右边：外部激励（电机转矩和负载转矩）
关键点：扭振的本质就是惯性力、弹性力和阻尼力之间的动态平衡。当外部激励的频率接近系统的固有频率时，就会发生共振。
2.4 多惯量系统的扩展
实际传动链往往不止两个惯量块。比如风电齿轮箱，有高速轴、中间轴、低速轴，每根轴都有自己的惯量和刚度。这时候就需要建立多惯量模型。
对于 n 个惯量块的系统，运动方程可以写成通用形式：
[J] * {d²θ/dt²} + [C] * {dθ/dt} + [K] * {θ} = {T}其中：
- [J]：惯量矩阵（对角矩阵）
- [C]：阻尼矩阵（三对角矩阵，相邻惯量块之间有耦合）
- [K]：刚度矩阵（也是三对角矩阵）
- {T}：外部转矩向量
这个方程的解，就是系统的扭振响应。求解方法有：
- 时域法：直接数值积分（比如龙格-库塔法），得到角位移、角速度随时间的变化
- 频域法：做拉普拉斯变换或傅里叶变换，得到系统的频率响应函数
我个人习惯先用频域法算固有频率和振型，看看哪些频率容易出问题。然后再用时域法模拟具体工况下的响应。
避坑指南：我曾经用四惯量模型算一个造纸机传动系统，结果固有频率算出来跟实测差了15%。后来发现，我忽略了纸卷的惯量是随着卷径变化的。所以，对于变惯量系统，一定要考虑惯量的时变性。否则模型就是错的。
2.5 知识体系总览
下面这张图把本章的核心逻辑串起来了。你可以看到，从物理系统到数学模型，再到求解方法，最后到工程应用，是一条完整的链路。
这张图把咱们今天讲的内容串起来了。从物理系统出发，经过多惯量建模和参数获取，得到运动方程，然后选择时域或频域方法求解，最后应用到工程实际中。每一步都有坑，每一步也都有技巧。
好了，关于扭振数学模型，咱们就聊到这儿。记住：模型是工具，不是真理。好的模型能帮你抓住主要矛盾，坏的模型只会让你在错误的路上越走越远。下次你在现场遇到扭振问题，不妨先画个多惯量模型，算算固有频率，再决定怎么下手。