企业官网建设流程全解析

简介：本资源是一份面向机械工程、动力学建模与传动系统分析领域的专业MATLAB建模工具，聚焦行星齿轮传动系统的扭振特性研究，适用于高校研究生、传动系统设计工程师及振动噪声分析从业者。资源解决的核心问题是行星齿轮机构在动态载荷下的扭转振动建模、模态识别与共振风险评估，尤其适用于汽车变速器、风电增速箱、机器人关节等高精度动力传输场景的动力学仿真需求。压缩包仅含1个核心文件NewMyModel.m（MATLAB脚本），大小3KB，完整实现了行星齿轮系（含太阳轮、行星轮、行星架、齿圈）的集中参数化建模、线性化动力学方程构建、质量-刚度-阻尼矩阵组装及特征值求解流程，可直接输出扭振固有频率与模态振型。已有367人学习下载，读者可快速复现行星齿轮传动扭振模态分析全过程，获取可修改的底层建模逻辑、清晰的物理量映射关系及标准化特征值求解框架，为故障诊断、结构优化与减振设计提供可靠理论支撑。

1. 动笔之前：齿轮扭振问题为什么值得单独建模

1.1 一个真实案例引发的思考

前段时间我接手了一套行星齿轮减速机构的振动问题分析。设备在试验台上跑起来之后，壳体上某个测点在中高速段出现明显的加速度峰，频谱上对应的是齿轮啮合频率附近的边带带。现场工程师第一反应是“轴承坏了”，换了轴承之后问题依旧。后来把齿轮箱拆开检查，齿面没有明显点蚀，间隙也合格，但壳体振动始终压不下去。

这个案子最后兜兜转转回到了动力学仿真上。我搭了一个以纯扭转为主的齿轮传动系统模型，把太阳轮、行星轮、行星架、齿圈完整地放进去，算它的动力学方程和扭振模态，然后和实测频谱对照，问题才真正定位出来：不是轴承失效，而是啮合频率的某倍频和系统某阶扭振固有频率对上了，形成了扭转共振。齿轮箱本身没坏，但整个系统在这一转速下被持续激励，振动自然压不住。

这就是齿轮传动扭振分析的典型价值。工程上很多异常振动、异响、断齿、齿面胶合，追根溯源都跟传动系统的扭振特性有关。摸清楚系统的固有频率和振型，搞清楚激励从哪来、频率值是多少，才能预判共振转速，而不是等问题在台架上暴露了才去猜。

1.2 扭振到底是个什么东西

扭振，简单说就是旋转体系绕轴线的扭转振动。别把它想得太玄，你可以把传动轴想象成一根橡皮筋，一端输入扭矩，一端拖着负载。扭矩不是瞬间传过去的，中间会有弹性扭转、回弹、再扭转的过程，表现出来就是各旋转部件之间的角速度此起彼伏，角位移彼此波动，这就是扭振。

对于齿轮传动系统，扭振分析最直接的办法是建立动力学方程：把每个齿轮看成具有转动惯量的刚体，把轮齿啮合看成弹簧和阻尼元件，把齿轮轴看成扭转弹簧，然后列出一组二阶常微分方程。解这组方程的特征值，就得到系统的扭振模态。模态这词听着高级，说白了就是系统在不同频率下最容易被“晃起来”的那种振动形态。频率和振型知道了，激励频率再一对，共振条件一目了然。

这篇内容适合谁看？我自己觉得，凡是做齿轮箱设计、传动系统NVH、旋转机械故障诊断的工程师，以及机械专业在读研究生，都可以参考。文章里我会把建模思路、方程怎么搭、参数怎么取、模态怎么解读、常见坑怎么避开，完整走一遍。用的是一套行星齿轮减速机构作为演示对象，但方法同样适用于定轴齿轮箱和其他旋转机械。

2. 行星齿轮和定轴齿轮在动力学上到底差在哪

2.1 行星轮系的结构特殊性

定轴齿轮箱，轴都是固定的，齿轮只是绕自己的轴转，建模相对直接：一对啮合齿轮给一个啮合刚度就完事。行星齿轮则完全不同，它有一个行星架带着好几个行星轮一边自转一边公转，太阳轮在中间，齿圈在最外面。这套结构把一条动力路径分成了多条并联功率流：太阳轮把扭矩分给所有行星轮，行星轮再合力把扭矩汇集到行星架或齿圈上。

并联带来的直接后果就是模态变得密集。同一套结构里因为多个行星轮完全对称，系统里会出现很多频率非常接近甚至成对出现的模态。做模态分析时你会发现，有些固有频率是成双成对的，振型形式上互为镜像；有些模态对应的所有行星轮振动相位完全一致，有些模态则是某些行星轮运动、其它行星轮静止。这些对称特性是定轴齿轮箱里不太会遇到的，也是行星齿轮扭振建模最需要注意的地方。

另外，行星齿轮系统里同时存在外啮合（太阳轮-行星轮）和内啮合（行星轮-齿圈）。外啮合和内啮合的激励特性不一样，轮齿接触应力分布、啮合刚度变化规律也有差异。建模的时候如果只用“一对齿轮啮合”的思路去套，很容易丢掉内啮合这一半的激励信息。

2.2 动力传递路径上的激励源不止一个

定轴齿轮箱的激励源主要是齿轮副啮合产生的传递误差，再叠加轴的弯曲、轴承游隙等因素。行星齿轮箱的激励源要多得多：太阳轮和每个行星轮之间有一组啮合激励，每个行星轮和齿圈之间又有一组啮合激励。行星轮数量一多，激励源数量直接翻倍。

更麻烦的是，各个行星轮与太阳轮的啮合并不是同步的。因为行星轮在圆周方向等间距布置，它们与太阳轮啮合的相位存在固定差值。这个相位差直接影响激励叠加效果。同相位叠加可能让振动增强，反相位叠加可能相互抵消。设计时行星轮个数、齿数选择，其实已经默默决定了这些相位关系。这也是为什么有些齿轮箱某些阶次特别安静，某些阶次怎么改进都压不下去——根子就在相位组合上。

还有一点经常被忽略：行星轮在公转，啮合点位置也在旋转。站在绝对坐标系里看，行星轮与太阳轮之间的啮合力方向是随时间变化的。纯扭转模型可以避开这个复杂问题，但如果要做弯扭耦合或壳体振动分析，这个方向变化就必须考虑进去。模型复杂度的取舍，要在动手前想清楚。

2.3 建模前先想清楚三个问题

我在搭模型之前通常会先问自己三个问题。第一，要解决什么问题？如果只关心系统的扭转固有频率和共振转速，纯扭转模型就够了，自由度少，参数少，出结果快。如果要分析壳体振动偏大、异响和齿轮箱结构耦合，就得至少加入轴和轴承的横向自由度。第二，精度要求是什么？前期方案设计阶段用平均啮合刚度完全够用，后期做故障诊断或优化降噪就需要时变啮合刚度加上传递误差激励。第三，能拿到的参数有哪些？有没有齿轮参数表、轴系尺寸、材料属性、实测频率数据？参数不全，模型再精细也是空中楼阁。

这三个问题想清楚了，再决定模型自由度和参数精度，就不容易走弯路。

3. 从零搭动力学方程：自由度、刚度、阻尼一个都不能少

3.1 建模基本假设与自由度选择

我这次演示用的是纯扭转行星齿轮模型。基本假设如下：

• 所有齿轮和行星架视为刚体，不考虑齿体弹性变形和轴的横向弯曲
• 轮齿啮合简化为沿着啮合线方向作用的弹簧-阻尼单元
• 忽略轴承游隙、齿侧间隙和轮齿摩擦
• 材料为线弹性，系统做小变形振动

这套假设对应一个 N+3 自由度的系统（齿圈固定时 N+2 个）：太阳轮转角 θ_s、行星架转角 θ_c、齿圈转角 θ_r、每个行星轮转角 θ_pi（i = 1, 2, ..., N）。我这里的演示模型取 N = 3，也就是三个行星轮。齿圈在实际结构中经常做成固定，也就是 θ_r 直接约束为 0，自由度进一步减少。

在列方程之前有个非常关键的操作：把所有角度量都转换成啮合线方向的位移。为什么要这么做？因为齿轮啮合的实质是两个轮齿在啮合线方向上的弹性接触，用角位移写方程，啮合刚度单位会变成 N·m/rad，还得换算；用啮合线位移写方程，所有啮合刚度的单位统一为 N/m，方程装配更直观。换算很简单，角位移乘以基圆半径就得到啮合线位移。

太阳轮和行星轮 i 之间的啮合线相对位移，我习惯写作：

δ_spi = r_bs · (θ_s - θ_c) - r_bp · (θ_pi - θ_c) + e_spi(t)

行星轮 i 和齿圈之间的啮合线相对位移写作：

δ_rpi = r_br · (θ_r - θ_c) + r_bp · (θ_pi - θ_c) + e_rpi(t)

其中 r_bs、r_bp、r_br 分别是太阳轮、行星轮、齿圈的基圆半径；e_spi、e_rpi 是综合啮合误差，作为激励源。这里正负号以我的坐标系和啮合方向定义为准，你建模时可以不同，但必须保证相对位移公式能正确反映“压缩为正”的物理事实。我见过很多建模错误都出在这一步符号上，所以建议列完公式后先做一个静态检查：给太阳轮一个微小正转角，看啮合位移是增大还是减小，是否符合实际。

3.2 参数怎么取：从设计参数表到自己算

参数是建模里最考验耐心的一步。先看齿轮基本参数。我演示用的行星轮系几何参数如下表，这个表按行业惯例可以从设计图纸或齿轮参数表中直接获得，单位为国际单位制。

基圆半径根据齿数和模数计算，标准齿轮的基圆半径公式为：

r_b = m · z · cos(α) / 2

太阳轮基圆半径 r_bs = 2 × 24 × cos20° / 2 ≈ 22.55 mm = 0.02255 m；行星轮 r_bp = 2 × 16 × cos20° / 2 ≈ 15.04 mm = 0.01504 m；齿圈 r_br = 2 × 56 × cos20° / 2 ≈ 52.63 mm = 0.05263 m。这些值会在装配方程时反复用到，建议一开始就算好，不要等写到一半再查。

啮合刚度是另一个关键参数。工程上常用 ISO 6336 或 GB/T 3480 的公式估算单齿啮合刚度，也可以用有限元算一对轮齿从啮入到啮出过程中的刚度变化曲线。我这里为了演示，取平均啮合刚度 k = 3 × 10^8 N/m，啮合阻尼比取 0.05。这个刚度值对钢制齿轮、20mm 齿宽来说处于合理量级，但你的具体模型一定要用自己的参数重算，不能照抄。

3.3 方程装配的实操细节

有了相对位移，运动方程就可以写了。以太阳轮为例，它的转动方程包含所有行星轮与它的啮合力力矩和输入扭矩：

J_s · θ_s'' + Σ_{i=1..N} r_bs · [k_sp · δ_spi + c_sp · δ_spi'] = T_s

行星架方程包含行星轴的支反力等效项和负载扭矩，行星轮方程包含它与太阳轮、齿圈两侧的啮合力力矩。把所有方程合在一起，写成矩阵形式就是：

M · x'' + C · x' + K · x = F(t)

其中 M 是质量矩阵（对角矩阵，元素是各构件的转动惯量），K 是刚度矩阵，由啮合刚度投影到各自由度后装配得到，C 是阻尼矩阵，F(t) 是包括输入扭矩、负载扭矩和传递误差激励的合力向量。如果齿圈固定，直接删掉 θ_r 对应的行列，同时把啮合位移公式里的 θ_r 置为 0。

方程装配完成之后，做模态分析之前，我一定会做一个静态合理性检查：把刚度矩阵乘以一个单位位移向量，看得到的力方向是否符合物理直觉；再对所有元素做量纲检查，确保每一项单位是 N·m/rad 而不是别的。这一步花十分钟，但能省掉后面排查问题的几个小时。

4. 扭振模态分析实操：特征值算完还要会解读

4.1 求解方法与行星齿轮模态的三大分类

忽略激励项，只研究自由振动，把位移设成 x = X·e^(ωt) 的形式代回方程，就得到广义特征值问题：

det(K - ω² · M) = 0

解这个方程可以得到系统的固有角频率 ω 和对应的振型向量 X。固有频率 f = ω / (2π)。实际求解可以直接用 MATLAB 的 eig 函数，或者用商业有限元软件。自由度不多时，MATLAB 手写也就几十行代码的事。

行星齿轮系统有个非常经典的理论结论，值得在这里提一下。按文献里对标准行星齿轮模型的分析，系统的模态可以分成三大类：

• 扭转模态：中心构件（太阳轮、齿圈、行星架）只发生扭转运动，所有行星轮的振动呈相同模式。这类模态在固有频率谱上往往是孤立的一支。
• 平移模态：中心构件发生横向平移，多个行星轮参与振动，而且频率成对出现，也就是重数为二。
• 行星轮模态：主要是行星轮之间的相对振动，中心构件几乎不动，频率的重数跟行星轮数量相关，当行星轮数为 N 时，有一部分模态重数表现为 N-3 的形式。

这里我不展开推导了，但记住这个分类对解读结果非常有帮助。比如你在扫频实验里看到频率靠得很近的一对峰，很可能就是平移模态；看到某个模态下中心轴振动很弱、只有行星轮在动，那多半是行星轮模态。

4.2 算一组结果的工程解读

按上面那组参数，齿圈固定，系统有 5 个自由度（太阳轮、行星架、3 个行星轮）。用平均啮合刚度 3×10^8 N/m、啮合阻尼比 0.05 做模态分析，第一阶扭转固有频率大致落在 1200 Hz 到 1500 Hz 之间，具体值取决于行星架的支撑刚度折算。

这个频率高不高？直接看数值没概念，要和激励频率对比。行星齿轮啮合频率计算公式为：

f_m = (n_s - n_c) · z_s / 60 = n_c · z_r / 60

太阳轮输入 1500 r/min，齿圈固定，传动比 i = 1 + z_r / z_s = 1 + 56/24 ≈ 3.333，行星架输出转速 n_c = 1500 / 3.333 ≈ 450 r/min。代入公式：

f_m = (1500 - 450) × 24 / 60 = 420 Hz

也就是说，在 1500 r/min 时，太阳轮-行星轮和行星轮-齿圈的啮合频率都是 420 Hz，其整数倍频是 840 Hz、1260 Hz、1680 Hz 等。如果第一阶扭转固有频率落在 1260 Hz 附近，那正好和三倍啮合频率重合，这个转速下就会出现明显的扭振放大。这就是我开头那个案例里实测频谱看到的现象。

每次算完固有频率，我都会把激励倍频列一张表对照，而不是只看孤立数字。固有频率和激励频率“撞上”才是问题的起点，单独讨论固有频率多高没有工程意义。

4.3 Campbell 图与共振转速判断

把固有频率和激励频率画在同一张图上，就成了工程上常用的 Campbell 图。横轴是输入转速，纵轴是频率。固有频率水平线若干条，啮合频率及其倍频是从原点出发的射线，转速变化时射线斜率固定。射线与水平线相交的点，就是潜在共振点。

判断共振条件时有几个坑。首先，啮合频率的倍频不只到三倍，十阶以内都要看，尤其在高转速下齿频激励的谐波分量可能很丰富。其次，不要只关注齿轮啮合激励，电机的电磁激励、负载侧的周期性扰动也要考虑。第三，行星齿轮由于多个行星轮分担载荷，实际激励中某些谐波会被相位条件削弱甚至消除，这在精确建模时可以用传递误差相位来体现。

图纸上看到交点后，还要估算共振裕度。工程上一般要求共振转速避开工作转速范围，或者至少留出 10% 以上的裕度。如果躲不开，就要从结构上调整刚度（比如改轴径、优化齿宽），或者通过增加阻尼来抑制共振峰值。

5. 建模和仿真中常见的坑，我替你踩过了

5.1 刚体模态和约束问题

最容易坑到新人的是系统里出现 0 频刚体模态。行星齿轮在纯扭转模型里如果没有施加任何接地约束，整个传动链可以自由旋转，刚度矩阵奇异，求特征值时前几阶接近 0。做模态分析之前一定要想清楚哪些构件被约束了。齿圈固定就把 θ_r 删掉；输入轴和负载轴如果没有接地约束，系统会多出刚体自由度，这在自由振动分析里是合理的，但如果你只关心弹性模态，就需要先把刚体模态剔除。

我自己习惯的做法是：在输入轴端和负载端分别加一个微小的接地扭转弹簧（比如相对于系统刚度小两三个数量级），这样既能抑制刚体模态，又不会显著影响弹性模态频率。但这只是处理手段，不是物理真实的模拟，报告里要写清楚。

5.2 单位转换和方向符号

单位错误是另一个高频翻车点。转速、角速度经常遇到 r/min 和 rad/s 混用。比如计算啮合频率时，速度单位必须用 r/s 或 r/min 配 60；算激励频率时，角频率单位必须用 rad/s。我刚做这套模型时就在这儿翻过车，算出来的啮合频率差了 9.55 倍，结果 Campbell 图上共振点全错位。

方向符号问题更隐蔽。行星轮同时和太阳轮、齿圈啮合，齿圈固定后，行星轮两侧的啮合力方向是相反的。列方程时如果不把符号统一，方程装配后刚度矩阵不对称，求出的固有频率会有虚部，结果一眼看上去没毛病，实际振型全不对。检查方法很简单：给某个自由度一个正向单位位移，手动算一遍所有啮合位移的变化方向，看是否和公式一致。

5.3 平均刚度 vs 时变刚度

纯扭转模型最常用平均啮合刚度，计算简单，适合做模态分析。但实际问题里啮合刚度是随啮合位置周期变化的，时变刚度会引入参数激励，可能产生组合共振现象。我实际测试过同一个模型：用平均刚度算出的固有频率和用时变刚度算出的主共振点基本一致，但时变刚度模型还能预测到组合共振（比如某一阶固有频率等于啮合频率的某一半或者差值），这在实际台架上有时候能看到，平均刚度模型就完全预测不到。

所以我的建议是：前期方案筛选用平均刚度，快速出结果；后期做共振校核、振动预测，至少要引入一阶谐波变化的时变啮合刚度。这个升级在 MATLAB 里实现并不难，把 K 矩阵里的恒定刚度 k 换成 k(t) = k_m + k_1 · cos(ω_m · t + φ) 即可，代价只是求解从特征值问题变成参数激励稳定性问题，计算时间变长。

5.4 模型验证：仿真结果和实测对不上怎么办

模型建得再精细，最终都要面对实测数据的检验。我通常按三条路径验证：一是对比固有频率，用锤击法或停机测试激励齿轮箱，看加速度频谱上的峰和模型固有频率是否对得上；二是对比稳态响应，在已知转速下实测壳体或轴端的振动频谱，看啮合频率及其边带的幅值趋势是否一致；三是对比振型，这个最难，通常需要多个测点才能判断。

如果对不上，先不要怀疑模型，先检查参数。转动惯量算没算错？齿轮宽度是按多少算的？支撑刚度折进去了吗？啮合刚度是不是取自手册但没有修正？这些参数经常能解释 20% 以上的偏差。确认参数无误之后再看模型自由度，很多时候是忽略了负载惯量或联轴器刚度。总的来说，模型和实测误差在 10% 以内就算可用，在 5% 以内已经是很好的结果。

6. 工具链选型与后续扩展建议

6.1 用 MATLAB 手写方程还是用商业软件

平心而论，工具选择取决于你的具体场景。系统自由度少、阶段是快速方案比较，我强烈建议直接用 MATLAB 或 Python 手写运动方程。自由度几十个以内，手写矩阵装配非常快，而且每行代码都对应明确的物理意义，调参方便，查错容易，你永远知道模型里有什么、没有什么。这也是我从一开始就推荐先建模再仿真的原因。

如果项目涉及复杂的轴系、轴承、箱体耦合，手写工作量就上来了。这时候可以考虑 Romax、Dyn或者类似软件，它们内置了齿轮建模、轴承刚度、传递误差计算，省去大量重复劳动。也可以先用多体动力学软件如 ADAMS 搭刚柔耦合模型，更直观，但参数标定周期长，前期学习成本高。

我的经验是：先用手写纯扭转模型把系统的频率特性摸清楚，再用商业软件做详细校核。两个模型结果一致，说明理解是正确的；结果不一致，往往意味着手写模型遗漏了关键自由度，或者商业软件的某些默认参数和你的假设不一致，这时候反而能发现更深层的问题。

6.2 沿着这个模型还能往哪走

这套行星齿轮扭振模型的扩展空间很大。一个方向是加入弯扭轴耦合，把轴的弯曲自由度、轴承刚度和扭转变形耦合起来，这样能分析更接近实际的齿轮箱振动，也能解释壳体振动问题。另一个方向是加入齿轮误差激励的统计特性，用实测传递误差数据代替简谐假设，预测不同制造精度下的振动水平，这直接和齿轮精度等级挂钩。

还有一个很实用的方向是参数灵敏度分析。改变太阳轮齿数、行星轮数量、齿宽、输入力矩，观察固有频率和响应幅值如何变化，找到最敏感的几何参数，再做设计优化。这种分析的价值在于：等台架出了问题再改齿轮已经来不及了，设计阶段就知道哪些参数对 NVH 影响大，项目成本会省很多。

我在实际项目里还有个体会：仿真结果最好和试验数据沉淀成一套企业内部的数据库。每次新项目做完，记录下齿轮参数、啮合刚度、实测固有频率和异常振动案例，积累两三个项目之后，建模的初始参数估计会越来越准，很多问题还没仿真就能提前预判。对做传动设计的团队来说，这套数据资产比单个仿真模型值钱得多。

我做齿轮箱扭振这些年，最大的心得是：模型永远是为工程问题服务的，不要为了追求高保真而把模型搞得越来越复杂。纯扭转模型不适合做噪声预测，但如果只关心共振转速和扭振模态，它已经是最好用的工具。先把简单模型吃透，再把自由度一个个加进去，每一步都和实测对比，这才是做动力学分析最稳妥的路子。

本文还有配套的精品资源，点击获取