企业官网建设流程全解析

1. 齿轮动力学纯扭转模型：先搞清楚它在解决什么问题

做齿轮传动系统分析这么多年，我遇到最多的情况是：拿着一个减速器或者增速箱的设计图纸，客户一句话就是"帮我算算这齿轮组稳不稳定，振动大不大"。但真要动手分析，你不可能一上来就把齿轮副、轴、轴承、箱体全做成三维弹性体——那个计算量直接起飞，而且大量参数你根本凑不齐。所以工程上最常用的做法，就是先用纯扭转模型把整个传动链的动力学底子打出来。这也是标题里"纯扭转"三个字的核心价值所在。

纯扭转模型说直白点，就是忽略轴的弯曲、箱体的弹性变形、轴承的柔性这些因素，把每一根轴只保留一个转动自由度，把齿轮、联轴器、负载这些质量体看成具有转动惯量的刚性圆盘，把齿面啮合、轴的扭转变形看成扭转弹簧和阻尼器。这样一来，整个齿轮传动系统就被简化成"多圆盘-弹簧-阻尼"的振动系统。虽然做了大幅简化，但在这个模型下，你能非常清楚地回答几个关键问题：系统的固有频率在哪个频段、啮合激励会引起多大的动态响应、动载系数大概是多少、转速在什么范围内会踩到共振点。这些结论在工程上非常有用，它决定了后面要不要做更精细的有限元分析，也决定了前期选型和避振方案的方向。

我个人的习惯是，拿到一个齿轮传动系统，第一版分析永远是纯扭转模型，基本一天之内就能在Simulink里搭完并跑出结果。等这个模型验证了基本设计合理，再考虑要不要上弯曲-扭转耦合模型，甚至是整机多体动力学。这件事背后的逻辑很简单：纯扭转模型是整个动力学分析的最小可行方案，但它的输出足够支撑大部分工程决策。

这篇博文适合谁看？我觉得三类人最需要：第一类是机械专业的研究生，做齿轮动力学课题起步时卡在建模思路上；第二类是搞传动系统设计的工程师，想快速评估一个设计方案，但又不想一上来就啃又贵又慢的商业软件；第三类是刚接触Simulink的仿真新人，想找一个不太难、又有明确物理意义的入门项目来练手。无论你是哪类，跟着这篇文章把模型搭起来、把结果解读出来，你就能建立一套从物理模型到仿真验证的完整思维方式。

2. 建模前的数学准备：公式推导与参数意义

2.1 集中质量法下的运动方程

纯扭转模型的数学基础是集中质量法。你想象把一根传动轴上的所有零件都"压扁"到轴线上，每个零件变成一个只有转动惯量的圆盘，零件之间通过轴的扭转刚度和啮合刚度连接。这样，整个系统就可以用一组二阶常微分方程描述，每个方程对应一个自由度的角位移或扭转角。

拿最常见的单级圆柱齿轮减速器来说，它简化以后就是两自由度系统：电机侧（主动轴）的转动惯量和负载侧（从动轴）的转动惯量，中间通过齿轮啮合刚度连起来。运动方程长这样：

[ J_p \ddot{\theta}_p + c_m r_p (r_p \dot{\theta}_p - r_g \dot{\theta}_g) + k_m(t) r_p (r_p \theta_p - r_g \theta_g) = T_p ]

[ J_g \ddot{\theta}_g + c_m r_g (r_g \dot{\theta}_g - r_p \dot{\theta}_p) + k_m(t) r_g (r_g \theta_g - r_p \theta_p) = -T_g ]

其中 (J_p)、(J_g) 分别是主动轮和从动轮（包含各自轴系）的转动惯量，(r_p)、(r_g) 是分度圆半径，(k_m(t)) 是时变啮合刚度，(c_m) 是啮合阻尼，(T_p)、(T_g) 是驱动转矩和负载转矩。

这两个方程本质上在描述什么？你可以把它类比成两个人中间拉着一根"弹性绳"在转身。你转得快了，绳子就绷紧，力量就传过去；你转得慢了，绳子就松，力量就减弱。齿轮啮合就是这根弹性绳，而绳子的松紧在我们这里被啮合刚度 (k_m) 控制。这个类比虽然粗，但特别有助于你记住模型的结构。

2.2 时变啮合刚度的两种处理方式

齿轮啮合刚度和普通弹簧最大的区别在于：它是时变的。因为齿轮在啮合过程中，参与啮合的齿对数是变化的——单齿啮合区和双齿啮合区交替出现，综合啮合刚度自然呈现周期性波动。这个波动正是齿轮系统振动的主要激励源之一，你要是把它当成常数，模型就算白建了。

我通常的处理思路有两种。第一种是查标准公式计算平均啮合刚度，然后在平均刚度基础上叠加一个按啮合频率变化的简谐波动，即：

[ k_m(t) = k_{avg} + k_a \cos(2\pi f_m t + \phi) ]

其中 (f_m = z \cdot n / 60) 是啮合频率（(z) 是齿数，(n) 是转速），(k_a) 是波动幅值，一般取平均刚度的20%~30%。这种方法的优点是简单直接，适合快速评估共振风险；缺点是没有捕捉到单双齿交替带来的非谐波特征。

第二种是更接近真实情况的方波/分段函数法：把一个啮合周期分成单齿啮合区和双齿啮合区，分别赋予不同的刚度值。比如单齿区刚度是 (k_1)，双齿区刚度是 (k_2 \approx 1.5k_1)，啮合周期的比例关系按重合度 (\varepsilon) 来分配。这样得到的刚度激励包含了更多频谱分量，仿真结果也更可信。

我建议做课题或者写论文的读者直接用第二种，因为审稿人看到你用时变刚度，第一反应就会问你"你的激励形式有没有包含单双齿交替特征"，用方波法至少这个点站得住。

2.3 阻尼的估算与取值陷阱

阻尼是动力学模型里最玄学的参数。它不像质量、刚度那样可以直接测量或精确计算，但它的取值直接影响共振峰的高度和系统衰减快慢。在啮合阻尼的估算上，我一般用下面这个经验公式：

[ c_m = 2\zeta \sqrt{\frac{k_m \cdot J_p \cdot J_g}{(r_g^2 J_p + r_p^2 J_g)}} ]

其中 (\zeta) 是阻尼比，齿轮啮合通常取0.03~0.1。这里有个陷阱：如果你分别给主动轴和从动轴的转动惯量设成相差很大的数值，直接套公式很容易算出负阻尼或者数值异常的阻尼系数。我实际测试过，当从动轮转动惯量比主动轮大一个数量级以上时，公式里的分母会非常大，算出来的阻尼等于几乎没有阻尼，共振峰值会异常尖。解决办法是适当调高阻尼比到0.1，或者改用基于能量等效的阻尼模型。这块我建议你在仿真脚本里把阻尼单独做成一个参数，方便反复调。

3. Simulink建模实操：从零搭一个纯扭转齿轮模型

3.1 整体架构设计：模块怎么划分

Simulink建模这件事，最怕的就是一上来就拖模块，拖到最后连线都理不清。我的习惯是先把系统的物理结构画成框图，然后在Simulink里按"一个物理对象一个子系统"的原则来组织。

对于纯扭转齿轮模型，我通常分成三个子模块：电机输入模块（提供驱动转矩）、齿轮啮合动力学模块（包含啮合刚度、阻尼、传动比关系）、负载模块（提供阻力矩）。其中，齿轮啮合动力学模块是整个模型的关键。

为什么要强调按物理对象划分子系统？因为齿轮系统最后通常要扩展到多级传动、轴系扭振、甚至和电机控制联合仿真。如果你把模型全堆在一个画布里，后续每一个扩展都会让你痛不欲生。而按物理对象划分后，每个子系统内部就是独立的"质量-刚度-阻尼"单元，新增一级齿轮传动，就相当于再复制粘贴一组子系统，方法没有任何变化。这一点在做完整个多级齿轮箱模型后体会尤其深。

一个典型的单级齿轮纯扭转模型在Simulink中的信号流是这样的：电机的驱动转矩作用于主动轮，主动轮的角速度通过齿轮啮合关系传递到从动轮，啮合刚度按时间或转角信号更新，从动轮的负载转矩作为反向激励参与计算，最终输出的信号是主动轮和从动轮的角速度、角位移，以及动态啮合力。

3.2 核心模块搭建细节

我在搭建时，具体操作步骤如下：

第一步，建立驱动输入。用一个恒定阶跃信号或者一个斜坡信号来模拟电机启动。很多人在这里直接用一个常数模块，但我建议用Step模块并设置合理的阶跃时间和终值。因为齿轮系统在启动瞬间会出现冲击响应，如果直接加一个阶跃转矩而不做任何平滑处理，仿真初期的数值振荡会非常明显。你可以把阶跃时间设置成0.1~0.2秒，让转矩在一个小时间段内爬升，这样既不会丢失启动冲击信息，又能让求解器稳定过渡。

第二步，搭建惯性环节。每个转动惯量 (J) 对应的传递函数是 (1/(Js))，也就是角加速度积分得到角速度，角速度再积分得到角位移。在Simulink中直接用Integrator模块级联两个积分器即可。这个环节的关键在于：主动轮和从动轮各有这么一组积分器，而且两组积分器之间必须通过啮合力互相耦合，不能独立运行。我见过很多人把两个轮子建成了两个互不相干的单自由度系统，那就是完全错误的模型。

第三步，实现啮合刚度。这是模型最核心的部分。如果采用简谐波动法，用Clock模块获取当前仿真时间，再通过一个Fcn模块或者MATLAB Function模块计算啮合刚度值。但要注意，啮合刚度应该以齿轮转角为自变量，而不是直接以时间为自变量。因为齿轮转速可能发生变化，如果直接按时间周期性变化，当转速波动时刚度变化频率会失真。更正的做法是：用主动轮的角位移和齿轮的齿数计算出啮合周期的相位，再套入刚度公式。

如果采用方波法，可以用MATLAB Function模块写一段判断逻辑：根据重合度 (\varepsilon) 判断当前时刻是单齿啮合还是双齿啮合，返回对应的刚度值。这段逻辑非常简单，但效果很真实。贴一段我常用的代码：

function k = mesh_stiffness(phase, eps, k1, k2) % phase: 归一化啮合相位 [0, 1) % eps: 重合度 % k1: 单齿啮合刚度, k2: 双齿啮合刚度 double_frac = eps - floor(eps); % 双齿啮合占一个周期的比例 if phase < double_frac k = k2; else k = k1; end end

这里有个细节需要注意：(k_2) 并不是 (k_1) 的两倍。双齿啮合时载荷由两对齿分担，但每对齿的刚度并不相同，且计入齿基变形后总刚度提升一般在1.3~1.7倍之间。我见过不少初学者上来就写 (k_2 = 2k_1)，仿真结果直接失真。推荐工程上取 (k_2 = 1.5k_1) 起步，然后根据具体齿轮参数修正。

第四步，搭建啮合力计算。动态啮合力可以写作：

[ F_m = k_m(t) \cdot (r_p\theta_p - r_g\theta_g) + c_m \cdot (r_p\dot{\theta}_p - r_g\dot{\theta}_g) ]

在Simulink中，用减法模块计算主动轮和从动轮在啮合线方向上的位移差和速度差，然后分别乘上刚度和阻尼，最后相加。这个力就是连接两个惯性环节的"桥梁"。我建议在力信号输出端接一个Scope或者To Workspace模块，把动态啮合力记录下来，这是后面分析动载系数的主要数据来源。

第五步，加载负载转矩。负载转矩可以是一个恒定值，也可以按照实际工况设置为转速的函数或时间函数。如果是恒定负载，直接加一个Constant模块取负号即可。如果想模拟变负载工况，比如碎石机、压缩机这类冲击负载，可以把负载设置成阶跃变化的序列。这一步看似简单，但对系统动态响应影响很大。恒定负载下系统稳态运行平稳，动态特性主要由啮合刚度激励主导；变负载下则会叠加外部激励的响应，分析的时候需要把两种激励分开来开。

3.3 仿真参数设置与求解器选择

很多初学者在Simulink里直接默认是用变步长求解器ode45，然后仿真结果一顿乱抖，就以为模型出问题了。实际上，对于齿轮动力学这类存在高频周期性激励的系统，求解器设置是有讲究的。

我的一般做法是：先用变步长求解器ode45跑一遍，步长上限设置为啮合频率对应周期的1/50左右。比如啮合频率是1000Hz，那Minimum step size和Max step size至少应该控制在2e-5秒的量级。如果你不限制最大步长，求解器会自动把步长拉大，导致高频响应丢失，你最后看到的曲线是严重失真的"平滑"曲线。

当变步长跑通之后，建议换成固定步长求解器ode4（四阶龙格库塔）验证一遍。固定步长取啮合频率对应周期的1/100，既能保证精度，又能在后续做实时仿真或硬件在环时保持一致的行为。这里有个实操经验：如果你的仿真耗时过长，先检查是不是最大步长设置得太小了。我曾经帮同事排查一个仿真卡死问题，发现他把最大步长设成了1e-7秒，明明只需要1e-5秒就能满足精度要求，白白慢了三个数量级。

仿真时间我建议设置为不少于20个啮合周期的时长。如果啮合频率1000Hz，那仿真时间至少0.02秒。但这只是最低要求，如果你想提取稳态响应，仿真时长要多留一些余量，一般取0.1秒以上。启动阶段的瞬态响应需要足够的衰减时间才能进入稳态，你要是只跑了0.02秒，看到的全是启动冲击的暂态成分，频谱分析出来什么结论都不可靠。

4. 仿真结果分析方法：从波形到力学结论

4.1 传动误差与动态啮合力解读

纯扭转模型搭好之后，第一步先看什么？我的习惯是先看传动误差曲线。传动误差定义为：

[ TE = r_p\theta_p - r_g\theta_g ]

这个量的物理意义非常直观：当着两个值完全相等时，说明齿轮副像刚性体一样同步转动，没有任何弹性变形；当它们出现偏差时，说明齿面发生了弹性变形和相对滑动，这正是振动和噪声的来源。

在仿真结果中，你会看到传动误差包含几个成分：一个恒定偏移（平均弹性变形量）、一个按啮合频率周期变化的波动（由时变刚度引起）、以及可能的低频成分（由系统固有振动引起）。如果你在Simulink里把传动误差的信号接入Scope，直接观察曲线的形状，就能初步判断模型行为是否合理。正常情况下，传动误差曲线应该呈现"锯齿状"或者"正弦波动状"，如果曲线是一条完美直线，那说明时变刚度没有被正确引入，模型退化成了定刚度系统。

动态啮合力是另一个核心输出。把 (F_m) 曲线拉出来看，它的均值应该接近名义载荷 (T_p / r_p)，波动幅值则反映了系统的动态放大效应。这里有一个非常实用的指标——动载系数 (K_v)，定义为动态啮合力的峰值与名义载荷的比值：

[ K_v = \frac{F_{m,max}}{T_p / r_p} ]

这个系数直接告诉你齿轮在实际运转中承受的最大载荷比理论静载荷高出多少。如果 (K_v) 超过了1.3，说明系统动态响应偏大，需要检查是否存在共振问题。按照齿轮设计标准，(K_v) 通常要求控制在1.05~1.3范围内，超出这个范围就要改参数。

4.2 频谱分析的具体操作

波形只能定性看，频谱才能定量分析。把动态啮合力或传动误差信号送到MATLAB工作空间，然后用fft函数做频谱分析。这是我的标准操作流程：

% 假设 simout 是 To Workspace 输出的结构体 t = simout.time; F = simout.signals.values; % 去除均值 F = F - mean(F); % FFT分析 Fs = 1 / (t(2) - t(1)); L = length(F); Y = fft(F); P2 = abs(Y / L); P1 = P2(1:L/2+1); P1(2:end-1) = 2*P1(2:end-1); f = Fs * (0:(L/2)) / L; % 绘制频谱 plot(f, P1); xlabel('频率 (Hz)'); ylabel('幅值'); xlim([0, max(f)/2]); grid on;

跑完这段代码，你会看到频谱上在啮合频率 (f_m)、其倍频 (2f_m)、(3f_m) 处有明显的峰值。这个现象本身就能告诉你：激励源就是啮合刚度的周期性波动。接下来要关注的是系统中是否在某个频率点出现异常大的峰值——特别是当啮合频率接近系统固有频率时，峰值会被大幅放大。如果你发现频谱在固有频率附近出现了比啮合频率基频还要高的峰，那基本可以确认系统运行在共振区或者接近共振区。

4.3 扫频分析找出共振转速

确认共振风险最直接的手段就是做扫频分析。方法很简单：保持负载不变，把输入转速从一个较低值逐步增加到较高值，在每个转速下仿真到稳态，然后记录该转速下动态啮合力的RMS值或峰值。最后把"转速-响应幅值"曲线画出来，出现明显峰值的转速就是临界转速。

这一步在Simulink里可以通过写一个for循环脚本来批量仿真，每次修改电机转速参数，然后批量跑仿真并收集结果。我以前是手动改一次参数跑一次，后来发现对于参数扫描这活完全可以自动化。写一个脚本：

speed_list = 500:50:3000; % 转速范围 r/min Kv_list = zeros(size(speed_list)); for i = 1:length(speed_list) % 更新 Simulink 工作区参数 in = Simulink.SimulationInput('gear_torsional_model'); in = in.setVariable('n_input', speed_list(i)); simOut = sim(in); % 提取动态啮合力并计算动载系数 Fm = simOut.DynamicMeshForce.signals.values; Kv_list(i) = max(Fm) / (T_p / r_p); end plot(speed_list, Kv_list); xlabel('输入转速 (r/min)'); ylabel('动载系数 Kv'); grid on;

当你画出这条速度-动载系数曲线时，你会发现存在几个明显的"山峰"。这些山峰对应的转速就是系统的临界转速。在这个转速附近运行，齿轮副的动应力会显著增大，严重时甚至引起齿面点蚀或断齿。知道这些临界转速，你就可以在设计中刻意避开它们，或者通过改变转动惯量、刚度来调整固有频率，把临界转速挪到工作转速范围之外。

5. 工程实践中的常见问题与排查技巧

5.1 仿真发散：最大步长和刚度波动是首要嫌疑

仿真一发散，新手第一反应是"模型错了"。但实际上，对于纯扭转齿轮模型，发散最常见的两个原因跟模型逻辑关系不大。第一个原因是最大步长设置过大，导致求解器在啮合刚度的突变点处跳过了关键信息，产生数值振荡甚至发散。因为时变刚度在单双齿交替的边界处有跳变，如果求解器在这个跳变附近的步长远大于跳变持续的时间尺度，数值误差就会被放大，最终导致发散。

排查方法是：把最大步长强制设置为啮合频率周期的1/100，再跑一遍。如果发散消失，说明问题就是步长不够小。第二个原因是阻尼系数设置错误。如果啮合阻尼 (c_m) 是负值或者数量级不对，系统的总阻尼可能出现负值，能量不断被注入系统，振幅指数增长，看起来就是发散。检查一下你公式中的各个参数单位是否统一，尤其是转动惯量用 (kg \cdot m^2)、刚度用 (N/m) 时，阻尼的单位必须是 (N \cdot s/m)，别把这里搞混了。

如果你已经排除了以上两点仍然发散，那就要检查初始条件了。我建议给主动轮的初始角速度设置为额定转速对应的值，不要从零开始。这样可以让仿真在接近稳态工作点的位置启动，极大降低启动冲击引发发散的几率。这一招我屡试不爽。

5.2 结果平滑得不像话：误以为模型正确

发散是极端问题，还有一种问题是"结果太完美"——传动误差曲线平滑得跟正弦波一样，每一周期完全重复，毫无噪声和波动特征。这种情况恰恰说明模型可能有问题。

最可能的原因是：时变刚度根本没被正确引入，系统实际上是在恒定刚度下运行。你在Simulink里检查一下刚度信号，如果刚度值是一条水平直线，说明你写的MATLAB Function模块里相位计算出了问题。最常见的错误是：你直接用了Clock模块的时间信号而没有除以啮合周期，导致 (phase) 变量一直超出 [0,1) 的范围，代码里判断单双齿啮合的逻辑永远走同一个分支，刚度自然就变成常数了。

另一个可能的原因是，输出信号的采样率不够，Scope的显示抽稀把高频波动隐藏掉了。这种问题看起来是结果光滑完美，实际上是数据采集的假象。解决方法是把数据通过To Workspace输出到MATLAB，然后用完整分辨率画图，不要只看Scope上的示意图。

5.3 参数单位与数量级校验清单

我想最后整理一个我在仿真之前必做的单位与数量级检查清单，帮助你避免很多莫名其妙的错误：

每次建完模型，先把这些参数过一遍，看看数量级是否合理。比如齿轮的啮合刚度不会到10⁵这个数量级以下，也不会超过10¹⁰这个数量级。如果你算出来的啮合刚度落在了异常区间，优先检查齿宽、模数、弹性模量这些基础参数有没有输错。

另外要特别提醒一个细节：Simulink模型中如果用了角度信号，默认单位是弧度而不是度。很多人把角度从Scope里导出来一看是几千的值，以为仿真出了问题，其实是因为角度累积超过了 (2\pi) 的周期。在后续信号处理中要注意对角度取模 (2\pi)，避免数值太大导致FFT分析时出现精度损失。传动误差的计算用的是弧度的线性组合，单位必须是弧度，否则算出来的刚度力就会差57.3倍，结果完全不可信。

5.4 如何验证模型可信度

仿真模型不是搭完跑出结果就完事，你还需要说服自己（以及以后可能的审稿人或同事）这个模型是可信的。我的验证思路分三步。

第一步是静力学验证：把时变刚度改成恒定刚度，阻尼设为零，输入恒定转矩，仿真到稳态后检查传动误差的均值是否等于 (T_p / (k_m r_p^2)) 的解析解。如果吻合，说明模型的基础受力逻辑正确。这一步看着简单，但能过滤掉一大半建模错误。

第二步是能量验证：在稳态下检查系统的输入功率是否等于输出功率加上阻尼耗散功率。在Simulink里，你可以分别计算驱动转矩做的功率 (T_p \cdot \omega_p) 和负载吸收的功率 (T_g \cdot \omega_g)，以及阻尼元件消耗的功率，做一个能量平衡表。如果三者不守恒，说明模型有能量泄漏或者莫名引入能量。这一步很有效，能揪出连接错误或者符号方向错误的问题。

第三步是对比验证：如果你有条件，用商业软件或者已有的实验数据做对比。没有条件也不要紧，至少可以把固有频率的结果用手算的扭转振动公式校核一遍。单级齿轮副的扭转固有频率可以用简化公式估算：

[ f_n = \frac{1}{2\pi} \sqrt{\frac{k_m}{J_e}} ]

其中等效转动惯量 (J_e = \frac{J_p J_g r_p^2}{J_p r_g^2 + J_g r_p^2})，这个值和Simulink模型频响中的峰值频率对得上，就说明模型的基本动力学行为是可信的。

6. 扩展思路：这个模型还能往哪里走

纯扭转模型是起点，不是终点。当你把单级齿轮副的纯扭转模型跑通之后，可以进行一系列扩展，不断逼近真实系统的复杂度。

最简单的扩展是多级齿轮传动。把两级或三级齿轮副的模型串联起来，每一个啮合点都有一对刚度-阻尼-传动比关系，级与级之间通过轴段扭转刚度连接。这时系统就有更多自由度，也会出现更多固有频率和共振转速。在Simulink里做这个扩展非常顺畅，因为你已经按物理对象划分子系统了，直接复制粘贴再连线就行。改参数的事情就交给数组和循环脚本处理。

再进一步，可以引入时变啮合位置和误差激励。真实的齿轮副不可避免存在齿距误差、齿形误差和齿面磨损，这些都会产生额外的位移激励。你可以把误差建模为叠加在传动误差信号上的周期函数，通过测量或者公差标准确定幅值，然后观察它对动态响应的影响。这一块对分析齿轮噪声、传动精度衰减非常有价值。

还有一个很实用的方向是和电机控制模型联合仿真。纯扭转模型加上电机电磁转矩模型和控制器模型，就能看到系统启动加速、突然加载、转速波动等瞬态过程中的齿轮动态响应。很多电动驱动系统的扭转振动问题（比如抖动、冲击噪声）必须在联合仿真中才能复现。在Simulink里，电机模型和齿轮模型天然兼容，我见过不少混动变速器标定团队就是用这套方法做换挡品质优化前期的初步分析。

不过我在实际中得到的最大体会是：纯扭转模型给你提供的不是"精确值"，而是"趋势判断"和"风险预警"。它能告诉你共振区间在哪里、动载系数在什么量级、整条传动链哪一级最脆弱，但它无法告诉你齿根应力峰值精确到多少兆帕，更无法告诉你齿面接触温度会不会超限。那些问题必须交给更精细的有限元分析和热分析。但反过来，如果你没有纯扭转模型提供的宏观动力学信息就直接上有限元，那你连边界条件都说不清楚，算出来的局部应力也不可信。合理的工程路线永远是：先全局、后局部；先扭转模型、后精细分析。

按这个思路，你能从一个相对简单的单级齿轮副模型出发，逐步构建起完整的传动系统仿真能力，并且每一步都有清晰的物理意义和工程价值。这也是我推荐每个做齿轮动力学的人从纯扭转模型入手的原因——它简单，但不简陋。