简介齿轮传动系统的振动分析是机械动力学领域的核心问题其中扭振由轮齿啮合弹性与转动惯量耦合产生直接决定齿轮箱运行平稳性与寿命。通过集中参数模型将太阳轮、行星轮等构件简化为惯量圆盘以弹簧-阻尼模拟啮合刚度可快速求解系统的固有频率与模态振型实现共振风险预测。该方法在风电、工程机械等行星齿轮传动装置的早期设计阶段具有重要价值能高效筛选临界转速并评估参数敏感性。本文以单级行星齿轮为例完整演示了从几何参数计算、刚度矩阵组装到MATLAB模态分析的全流程为齿轮箱扭振校核提供可复现代码与工程参考。 齿轮箱跑起来之后在某些转速区间突然噪声变大、振动骤增拆开一看齿面已经出现微点蚀——这种场景对做传动系统的人一定不陌生。我前阵子用自建的NewMyModel模型对一台单级行星齿轮传动系统做扭振动力学建模和模态分析从推导动力学方程到跑出完整的扭振模态参数整个过程踩了不少坑。这篇文章把完整思路、为什么这样建模、以及可直接复现的代码全部分享出来适合研究齿轮传动动力学的研究生、刚接触行星齿轮振动分析的工程师以及正在做齿轮箱扭振校核的机械设计人员参考。看完你不仅能自己搭建一个纯扭转集中参数模型还能系统地分析齿轮传动系统的共振风险。这套方法的核心思路是行星齿轮传动系统在扭矩传递过程中由于轮齿啮合弹性、各构件转动惯量和周期性激励的共同作用会表现出明显的扭转振动特性。用集中质量法把太阳轮、行星轮、齿圈、行星架简化为具有转动惯量的刚性圆盘把轮齿啮合简化为弹簧-阻尼元件建立动力学方程后求解特征值问题就能得到系统的固有频率和扭振模态。这套流程在工程上属于性价比极高的方案——比有限元建模快得多又能回答该不该避开某个转速这种实际问题。1. 行星齿轮传动为什么会扭振从结构到工程事故1.1 行星齿轮的基本结构和运动关系单级行星齿轮传动由一个太阳轮、若干个行星轮、一个内齿圈和一个行星架组成。典型的构型是太阳轮作为输入行星架作为输出齿圈固定不动。行星轮均匀分布在太阳轮和齿圈之间同时与两者啮合实现功率分流。在齿数设计上有一个基本约束叫同心条件Z_r Z_s 2 Z_p也就是齿圈齿数等于太阳轮齿数加上两倍行星轮齿数。例如太阳轮20齿、行星轮34齿、齿圈88齿正好满足202×3488。这个条件保证了三个齿轮的中心线能够重合行星轮才能同时与太阳轮和齿圈正确啮合。传动比公式也很关键。对于太阳轮输入、行星架输出、齿圈固定的情况i 1 Z_r / Z_s以上面这套齿数计算传动比就是188/205.4。行星齿轮之所以在风电、工程机械、汽车自动变速器里被广泛使用就是因为在很小的体积内能实现大传动比而且多行星轮分流让每个啮合点承受的载荷大大降低。1.2 扭振的物理本质弹性系统加周期性激励要理解扭振先把机械系统抽象一下。齿轮箱里的每一个构件都不是绝对的刚体轮齿啮合处有接触弹性变形轴有扭转柔度这些弹性元件把各个转动惯量串在一起构成了一个典型的多自由度扭转振动系统。拿生活中的例子类比你在秋千上坐着别人每荡到最低点推一把如果推的频率恰好等于秋千的固有频率秋千就会越荡越高。齿轮箱里的情况本质上完全相同——啮合激励就是那只手系统固有频率就是秋千的频率两者一旦对上振幅就会被放大轮齿动载荷成倍增加。行星齿轮系统和普通平行轴齿轮副的一个关键差异在于多个行星轮同时参与啮合太阳轮、行星架、齿圈之间的刚度耦合是高度交互的。一方面功率分流让系统更紧凑另一方面也意味着任何一个构件的振动都可能通过多路啮合路径传递给其他构件。这就让行星齿轮的扭转模态比普通圆柱齿轮箱复杂得多。1.3 激励源有哪些以及共振的破坏形式齿轮传动的振动激励主要分三类时变啮合刚度激励齿轮啮合过程中单齿对啮合和双齿对啮合交替进行啮合刚度周期性变化这是最主要的内部激励。它的基频就是啮合频率。传动误差激励轮齿的制造误差、装配误差、齿面变形、修形量都会让实际啮合位置偏离理论位置形成位移激励。外部负载波动负载扭矩的变化、原动机的输出波动、启动冲击也会激发系统振动。这些激励一旦和系统某一阶固有频率重合就会引发共振。共振带来的直接后果包括动载系数急剧增大、齿面接触应力超标导致点蚀或胶合、轮齿弯曲疲劳断裂、轴承过早失效以及齿轮箱辐射噪声增大。在工程现场最常见的表现就是某个转速区间噪声异常明显出问题往往就发生在这些转速点上。2. 建模思路与动力学方程推导从物理到数学2.1 为什么最终选择了纯扭转集中参数模型行星齿轮系统的动力学建模大致可以分三个台阶纯扭转集中参数模型每个构件只保留一个旋转自由度建模最快适合做固有频率筛选和参数敏感性分析。平面-扭转集中参数模型每个构件增加两个平移自由度可以分析支撑刚度、偏心、不对中等因素的影响是学术研究和工程校核中很常用的模型。有限元模型完整考虑齿轮体弹性、轮齿变形、箱体刚度精度最高但建模成本也最高通常只用于最终方案的精细校核。我这次选择纯扭转模型原因很实际首先本次关注的核心问题就是传动系统的扭转共振重点在于各构件之间的扭转耦合关系其次概念设计阶段往往需要快速比较多组齿数、惯量参数每改一次参数就重跑一次有限元不现实。纯扭转模型的自由度数量级只有几个到十几个求解特征值问题几乎不花时间参数扫描效率极高。当然要到后期确认具体结构的应力水平还是得回到有限元。但先说清楚做扭振模态筛选和临界转速评估纯扭转模型在工程上是足够靠谱的前提是建模时物理参数给得准。2.2 自由度定义与啮合变形表达先定义广义坐标。对于Np个行星轮的单级行星轮系齿圈固定取广义坐标向量q [θ_s, θ_p1, θ_p2, ..., θ_pNp, θ_c]ᵀ其中θ_s为太阳轮转角θ_pi为第i个行星轮转角θ_c为行星架转角。齿圈固定所以没有θ_r自由度。为了把旋转自由度映射到啮合线上的弹性变形需要把转角乘以基圆半径换算成切向线位移u_s r_bs · θ_su_p r_bp · θ_pu_c r_c · θ_c这里r_bs、r_bp是太阳轮和行星轮的基圆半径r_c是行星架回转中心到行星轮中心的距离即中心距。接下来是关键一步定义啮合线上的相对位移以齿面受压方向为正。在工程上常用的近似表达为δ_spi r_bs · θ_s - r_bp · θ_pi - r_c · θ_cδ_rpi r_bp · θ_pi - r_c · θ_c第二条式子表达的是行星轮与固定齿圈之间的啮合变形齿圈固定所以没有θ_r项行星架转动会带着行星轮整体离开或靠近啮合点这就是-r_c·θ_c项的物理意义。这里要特别提醒不同文献里符号约定差异很大有的把受压方向定义为正有的定义为负。这直接影响动力学方程中力的方向进而影响刚度矩阵的对称性。在做推导之前一定要自己把几何关系画一遍明确每个变形的正方向。这个坑我后面会专门展开讲。2.3 用拉格朗日方程建立运动方程有了广义坐标和啮合变形表达式就可以写系统的动能和势能。系统动能T 0.5·J_s·θ̇_s² 0.5·ΣJ_p·θ̇_pi² 0.5·J_c·θ̇_c²系统弹性势能忽略重力势能只考虑啮合弹簧的变形能V 0.5·Σk_sp·δ_spi² 0.5·Σk_rp·δ_rpi²其中k_sp、k_rp分别是太阳轮-行星轮啮合副和行星轮-齿圈啮合副的等效啮合刚度单位是N/m。把T和V代入拉格朗日方程d/dt(∂T/∂q̇) - ∂T/∂q ∂V/∂q Q就可以得到系统的运动微分方程整理成矩阵形式为M·q̈ K·q T(t)这里没有显式写出阻尼矩阵C因为模态分析求解固有频率和振型时阻尼的影响可以忽略到做谐响应或时域响应时才需要给阻尼。工程上常用比例阻尼处理即C α·M β·K两个系数用实测或经验模态阻尼比反推。2.4 质量、刚度和阻尼矩阵怎么组装根据拉格朗日方程推导刚度矩阵K各元素就是势能对广义坐标的二阶偏导K_ij ∂²V / ∂q_i ∂q_j以Np3为例系统共有5个自由度质量矩阵是对角阵M diag(J_s, J_p, J_p, J_p, J_c)刚度矩阵是一个5×5对称矩阵每个位置的具体表达式如下矩阵位置刚度表达式K_θs,θsNp·k_sp·r_bs²K_θs,θpi-k_sp·r_bs·r_bpK_θs,θc-Np·k_sp·r_bs·r_cK_θpi,θpir_bp²·(k_sp k_rp)K_θpi,θc-r_bp·r_c·(k_sp k_rp)K_θc,θcNp·r_c²·(k_sp k_rp)非对角线行星轮之间0从矩阵结构可以直观看出行星齿轮扭振模型的几个特点太阳轮和行星架的刚度项都带有Np的倍数因子意味着行星轮个数增多会直接提高中心构件的刚度贡献行星轮之间在刚度矩阵里没有直接耦合它们通过太阳轮和行星架间接联系。阻尼矩阵如果不采用比例阻尼可以在每个啮合副上单独加阻尼力。啮合阻尼系数c_sp、c_rp与啮合刚度之间有关系c 2·ζ·sqrt(k·m_eq)ζ一般取0.01到0.1。注意这个公式里的m_eq是等效质量需要用两个构件的等效惯量折算。3. NewMyModel模型落地参数计算与MATLAB实现3.1 几何参数计算齿数、基圆半径、中心距与传动比模型的几何参数是一切动力学计算的基础。这一步最容易出错我把完整的计算流程写在代码里直接用SI单位制避免单位混乱。%% NewMyModel_GeomAndModal.m % 单级行星齿轮纯扭转动力学模型 % 所有单位采用SI米、千克、牛、秒、赫兹 clear; clc; close all; % ---- 1. 基本齿数参数 ---- Z_s 20; % 太阳轮齿数 Z_p 34; % 行星轮齿数 Z_r 88; % 齿圈齿数 Np 3; % 行星轮个数 % 同心条件校验 if Z_r ~ Z_s 2*Z_p error(齿数不满足同心条件); end % ---- 2. 尺寸参数 ---- m 2e-3; % 模数注意单位是米 alpha 20; % 压力角度 b 30e-3; % 齿宽米 rho 7850; % 材料密度钢 kg/m^3 % 分度圆半径 r_s m * Z_s / 2; r_p m * Z_p / 2; r_r m * Z_r / 2; % 基圆半径 r_bs r_s * cosd(alpha); r_bp r_p * cosd(alpha); r_br r_r * cosd(alpha); % 中心距行星架回转半径 r_c r_s r_p; % 传动比太阳轮输入行星架输出齿圈固定 i_ratio 1 Z_r / Z_s; fprintf(传动比 i %.2f\n, i_ratio);这里有一个非常值得注意的细节模数m直接使用2e-3米而不是2毫米。许多人在写代码时习惯用毫米后面算惯量时半径、宽度、密度混在一起单位一旦不统一计算结果会错一个甚至几个数量级而且很难排查。我的习惯是一次性把长度全部折算成米。3.2 转动惯量估算与啮合刚度选取转动惯量如果齿轮是实心结构可以直接按圆柱体公式估算J 0.5·ρ·π·b·r⁴行星架则不能简单按实心圆盘估算因为它带着行星轮轴、轴承座和支撑筋结构实际惯量通常比实心盘小但也要考虑行星轮绕自身轴线的惯量传递。工程上常用等效系数处理我这次取1.3倍的圆盘惯量作为行星架等效惯量因为行星架上还有三个行星轮轴和额外的结构质量。更精确的做法是把行星架分解成几个圆柱和方块的惯量叠加但评估阶段没必要那么细。啮合刚度选取上根据ISO 6336的经验数据钢制齿轮副沿啮合线的啮合刚度大约在10~20 N/(mm·μm)每毫米齿宽。取中间值14乘以齿宽30mm再换算到SI单位k 14 × 30 420 N/μm 4.2e8 N/m这就是单对齿啮合的等效刚度。太阳轮-行星轮啮合副和行星轮-齿圈啮合副在工程上常用同一个值近似。% ---- 3. 惯量计算 ---- J_s 0.5 * pi * rho * b * r_s^4; J_p 0.5 * pi * rho * b * r_p^4; J_c 1.3 * (0.5 * pi * rho * b * r_c^4); % 经验系数1.3 % ---- 4. 啮合刚度估算 ---- % 经验值钢制齿轮副啮合线刚度约14 N/(mm*um) k_sp 14 * (b * 1e3) * 1e6; % 转换为 N/m k_rp k_sp; fprintf(太阳轮惯量 J_s %.3e kg·m^2\n, J_s); fprintf(行星轮惯量 J_p %.3e kg·m^2\n, J_p); fprintf(行星架惯量 J_c %.3e kg·m^2\n, J_c); fprintf(啮合刚度 k %.3e N/m\n, k_sp);3.3本文还有配套的精品资源点击获取 行星齿轮扭振动力学建模与模态分析实战
简介齿轮传动系统的振动分析是机械动力学领域的核心问题其中扭振由轮齿啮合弹性与转动惯量耦合产生直接决定齿轮箱运行平稳性与寿命。通过集中参数模型将太阳轮、行星轮等构件简化为惯量圆盘以弹簧-阻尼模拟啮合刚度可快速求解系统的固有频率与模态振型实现共振风险预测。该方法在风电、工程机械等行星齿轮传动装置的早期设计阶段具有重要价值能高效筛选临界转速并评估参数敏感性。本文以单级行星齿轮为例完整演示了从几何参数计算、刚度矩阵组装到MATLAB模态分析的全流程为齿轮箱扭振校核提供可复现代码与工程参考。 齿轮箱跑起来之后在某些转速区间突然噪声变大、振动骤增拆开一看齿面已经出现微点蚀——这种场景对做传动系统的人一定不陌生。我前阵子用自建的NewMyModel模型对一台单级行星齿轮传动系统做扭振动力学建模和模态分析从推导动力学方程到跑出完整的扭振模态参数整个过程踩了不少坑。这篇文章把完整思路、为什么这样建模、以及可直接复现的代码全部分享出来适合研究齿轮传动动力学的研究生、刚接触行星齿轮振动分析的工程师以及正在做齿轮箱扭振校核的机械设计人员参考。看完你不仅能自己搭建一个纯扭转集中参数模型还能系统地分析齿轮传动系统的共振风险。这套方法的核心思路是行星齿轮传动系统在扭矩传递过程中由于轮齿啮合弹性、各构件转动惯量和周期性激励的共同作用会表现出明显的扭转振动特性。用集中质量法把太阳轮、行星轮、齿圈、行星架简化为具有转动惯量的刚性圆盘把轮齿啮合简化为弹簧-阻尼元件建立动力学方程后求解特征值问题就能得到系统的固有频率和扭振模态。这套流程在工程上属于性价比极高的方案——比有限元建模快得多又能回答该不该避开某个转速这种实际问题。1. 行星齿轮传动为什么会扭振从结构到工程事故1.1 行星齿轮的基本结构和运动关系单级行星齿轮传动由一个太阳轮、若干个行星轮、一个内齿圈和一个行星架组成。典型的构型是太阳轮作为输入行星架作为输出齿圈固定不动。行星轮均匀分布在太阳轮和齿圈之间同时与两者啮合实现功率分流。在齿数设计上有一个基本约束叫同心条件Z_r Z_s 2 Z_p也就是齿圈齿数等于太阳轮齿数加上两倍行星轮齿数。例如太阳轮20齿、行星轮34齿、齿圈88齿正好满足202×3488。这个条件保证了三个齿轮的中心线能够重合行星轮才能同时与太阳轮和齿圈正确啮合。传动比公式也很关键。对于太阳轮输入、行星架输出、齿圈固定的情况i 1 Z_r / Z_s以上面这套齿数计算传动比就是188/205.4。行星齿轮之所以在风电、工程机械、汽车自动变速器里被广泛使用就是因为在很小的体积内能实现大传动比而且多行星轮分流让每个啮合点承受的载荷大大降低。1.2 扭振的物理本质弹性系统加周期性激励要理解扭振先把机械系统抽象一下。齿轮箱里的每一个构件都不是绝对的刚体轮齿啮合处有接触弹性变形轴有扭转柔度这些弹性元件把各个转动惯量串在一起构成了一个典型的多自由度扭转振动系统。拿生活中的例子类比你在秋千上坐着别人每荡到最低点推一把如果推的频率恰好等于秋千的固有频率秋千就会越荡越高。齿轮箱里的情况本质上完全相同——啮合激励就是那只手系统固有频率就是秋千的频率两者一旦对上振幅就会被放大轮齿动载荷成倍增加。行星齿轮系统和普通平行轴齿轮副的一个关键差异在于多个行星轮同时参与啮合太阳轮、行星架、齿圈之间的刚度耦合是高度交互的。一方面功率分流让系统更紧凑另一方面也意味着任何一个构件的振动都可能通过多路啮合路径传递给其他构件。这就让行星齿轮的扭转模态比普通圆柱齿轮箱复杂得多。1.3 激励源有哪些以及共振的破坏形式齿轮传动的振动激励主要分三类时变啮合刚度激励齿轮啮合过程中单齿对啮合和双齿对啮合交替进行啮合刚度周期性变化这是最主要的内部激励。它的基频就是啮合频率。传动误差激励轮齿的制造误差、装配误差、齿面变形、修形量都会让实际啮合位置偏离理论位置形成位移激励。外部负载波动负载扭矩的变化、原动机的输出波动、启动冲击也会激发系统振动。这些激励一旦和系统某一阶固有频率重合就会引发共振。共振带来的直接后果包括动载系数急剧增大、齿面接触应力超标导致点蚀或胶合、轮齿弯曲疲劳断裂、轴承过早失效以及齿轮箱辐射噪声增大。在工程现场最常见的表现就是某个转速区间噪声异常明显出问题往往就发生在这些转速点上。2. 建模思路与动力学方程推导从物理到数学2.1 为什么最终选择了纯扭转集中参数模型行星齿轮系统的动力学建模大致可以分三个台阶纯扭转集中参数模型每个构件只保留一个旋转自由度建模最快适合做固有频率筛选和参数敏感性分析。平面-扭转集中参数模型每个构件增加两个平移自由度可以分析支撑刚度、偏心、不对中等因素的影响是学术研究和工程校核中很常用的模型。有限元模型完整考虑齿轮体弹性、轮齿变形、箱体刚度精度最高但建模成本也最高通常只用于最终方案的精细校核。我这次选择纯扭转模型原因很实际首先本次关注的核心问题就是传动系统的扭转共振重点在于各构件之间的扭转耦合关系其次概念设计阶段往往需要快速比较多组齿数、惯量参数每改一次参数就重跑一次有限元不现实。纯扭转模型的自由度数量级只有几个到十几个求解特征值问题几乎不花时间参数扫描效率极高。当然要到后期确认具体结构的应力水平还是得回到有限元。但先说清楚做扭振模态筛选和临界转速评估纯扭转模型在工程上是足够靠谱的前提是建模时物理参数给得准。2.2 自由度定义与啮合变形表达先定义广义坐标。对于Np个行星轮的单级行星轮系齿圈固定取广义坐标向量q [θ_s, θ_p1, θ_p2, ..., θ_pNp, θ_c]ᵀ其中θ_s为太阳轮转角θ_pi为第i个行星轮转角θ_c为行星架转角。齿圈固定所以没有θ_r自由度。为了把旋转自由度映射到啮合线上的弹性变形需要把转角乘以基圆半径换算成切向线位移u_s r_bs · θ_su_p r_bp · θ_pu_c r_c · θ_c这里r_bs、r_bp是太阳轮和行星轮的基圆半径r_c是行星架回转中心到行星轮中心的距离即中心距。接下来是关键一步定义啮合线上的相对位移以齿面受压方向为正。在工程上常用的近似表达为δ_spi r_bs · θ_s - r_bp · θ_pi - r_c · θ_cδ_rpi r_bp · θ_pi - r_c · θ_c第二条式子表达的是行星轮与固定齿圈之间的啮合变形齿圈固定所以没有θ_r项行星架转动会带着行星轮整体离开或靠近啮合点这就是-r_c·θ_c项的物理意义。这里要特别提醒不同文献里符号约定差异很大有的把受压方向定义为正有的定义为负。这直接影响动力学方程中力的方向进而影响刚度矩阵的对称性。在做推导之前一定要自己把几何关系画一遍明确每个变形的正方向。这个坑我后面会专门展开讲。2.3 用拉格朗日方程建立运动方程有了广义坐标和啮合变形表达式就可以写系统的动能和势能。系统动能T 0.5·J_s·θ̇_s² 0.5·ΣJ_p·θ̇_pi² 0.5·J_c·θ̇_c²系统弹性势能忽略重力势能只考虑啮合弹簧的变形能V 0.5·Σk_sp·δ_spi² 0.5·Σk_rp·δ_rpi²其中k_sp、k_rp分别是太阳轮-行星轮啮合副和行星轮-齿圈啮合副的等效啮合刚度单位是N/m。把T和V代入拉格朗日方程d/dt(∂T/∂q̇) - ∂T/∂q ∂V/∂q Q就可以得到系统的运动微分方程整理成矩阵形式为M·q̈ K·q T(t)这里没有显式写出阻尼矩阵C因为模态分析求解固有频率和振型时阻尼的影响可以忽略到做谐响应或时域响应时才需要给阻尼。工程上常用比例阻尼处理即C α·M β·K两个系数用实测或经验模态阻尼比反推。2.4 质量、刚度和阻尼矩阵怎么组装根据拉格朗日方程推导刚度矩阵K各元素就是势能对广义坐标的二阶偏导K_ij ∂²V / ∂q_i ∂q_j以Np3为例系统共有5个自由度质量矩阵是对角阵M diag(J_s, J_p, J_p, J_p, J_c)刚度矩阵是一个5×5对称矩阵每个位置的具体表达式如下矩阵位置刚度表达式K_θs,θsNp·k_sp·r_bs²K_θs,θpi-k_sp·r_bs·r_bpK_θs,θc-Np·k_sp·r_bs·r_cK_θpi,θpir_bp²·(k_sp k_rp)K_θpi,θc-r_bp·r_c·(k_sp k_rp)K_θc,θcNp·r_c²·(k_sp k_rp)非对角线行星轮之间0从矩阵结构可以直观看出行星齿轮扭振模型的几个特点太阳轮和行星架的刚度项都带有Np的倍数因子意味着行星轮个数增多会直接提高中心构件的刚度贡献行星轮之间在刚度矩阵里没有直接耦合它们通过太阳轮和行星架间接联系。阻尼矩阵如果不采用比例阻尼可以在每个啮合副上单独加阻尼力。啮合阻尼系数c_sp、c_rp与啮合刚度之间有关系c 2·ζ·sqrt(k·m_eq)ζ一般取0.01到0.1。注意这个公式里的m_eq是等效质量需要用两个构件的等效惯量折算。3. NewMyModel模型落地参数计算与MATLAB实现3.1 几何参数计算齿数、基圆半径、中心距与传动比模型的几何参数是一切动力学计算的基础。这一步最容易出错我把完整的计算流程写在代码里直接用SI单位制避免单位混乱。%% NewMyModel_GeomAndModal.m % 单级行星齿轮纯扭转动力学模型 % 所有单位采用SI米、千克、牛、秒、赫兹 clear; clc; close all; % ---- 1. 基本齿数参数 ---- Z_s 20; % 太阳轮齿数 Z_p 34; % 行星轮齿数 Z_r 88; % 齿圈齿数 Np 3; % 行星轮个数 % 同心条件校验 if Z_r ~ Z_s 2*Z_p error(齿数不满足同心条件); end % ---- 2. 尺寸参数 ---- m 2e-3; % 模数注意单位是米 alpha 20; % 压力角度 b 30e-3; % 齿宽米 rho 7850; % 材料密度钢 kg/m^3 % 分度圆半径 r_s m * Z_s / 2; r_p m * Z_p / 2; r_r m * Z_r / 2; % 基圆半径 r_bs r_s * cosd(alpha); r_bp r_p * cosd(alpha); r_br r_r * cosd(alpha); % 中心距行星架回转半径 r_c r_s r_p; % 传动比太阳轮输入行星架输出齿圈固定 i_ratio 1 Z_r / Z_s; fprintf(传动比 i %.2f\n, i_ratio);这里有一个非常值得注意的细节模数m直接使用2e-3米而不是2毫米。许多人在写代码时习惯用毫米后面算惯量时半径、宽度、密度混在一起单位一旦不统一计算结果会错一个甚至几个数量级而且很难排查。我的习惯是一次性把长度全部折算成米。3.2 转动惯量估算与啮合刚度选取转动惯量如果齿轮是实心结构可以直接按圆柱体公式估算J 0.5·ρ·π·b·r⁴行星架则不能简单按实心圆盘估算因为它带着行星轮轴、轴承座和支撑筋结构实际惯量通常比实心盘小但也要考虑行星轮绕自身轴线的惯量传递。工程上常用等效系数处理我这次取1.3倍的圆盘惯量作为行星架等效惯量因为行星架上还有三个行星轮轴和额外的结构质量。更精确的做法是把行星架分解成几个圆柱和方块的惯量叠加但评估阶段没必要那么细。啮合刚度选取上根据ISO 6336的经验数据钢制齿轮副沿啮合线的啮合刚度大约在10~20 N/(mm·μm)每毫米齿宽。取中间值14乘以齿宽30mm再换算到SI单位k 14 × 30 420 N/μm 4.2e8 N/m这就是单对齿啮合的等效刚度。太阳轮-行星轮啮合副和行星轮-齿圈啮合副在工程上常用同一个值近似。% ---- 3. 惯量计算 ---- J_s 0.5 * pi * rho * b * r_s^4; J_p 0.5 * pi * rho * b * r_p^4; J_c 1.3 * (0.5 * pi * rho * b * r_c^4); % 经验系数1.3 % ---- 4. 啮合刚度估算 ---- % 经验值钢制齿轮副啮合线刚度约14 N/(mm*um) k_sp 14 * (b * 1e3) * 1e6; % 转换为 N/m k_rp k_sp; fprintf(太阳轮惯量 J_s %.3e kg·m^2\n, J_s); fprintf(行星轮惯量 J_p %.3e kg·m^2\n, J_p); fprintf(行星架惯量 J_c %.3e kg·m^2\n, J_c); fprintf(啮合刚度 k %.3e N/m\n, k_sp);3.3本文还有配套的精品资源点击获取