采用势能法编写的行星齿轮内啮合齿轮副时变啮合刚度程序(健康齿)程序中考虑了精确的渐开线齿形以及齿轮变位内齿圈固定行星架旋转同时考虑了考虑各啮合齿轮副之间的相位差可提供相位差计算小程序 假设内齿轮基体为刚性一、程序概述本程序套件基于势能法开发专注于计算健康状态下行星齿轮内啮合副的时变啮合刚度TVMS核心针对“行星轮-内齿圈”啮合场景。程序严格遵循齿轮啮合力学原理考虑了齿轮传动的关键几何与物理特性同时适配行星轮系的运动特性可用于齿轮动力学分析、振动噪声预测及故障诊断的基础刚度建模。1.1 核心功能定位计算对象行星齿轮系统中“行星轮主动轮-内齿圈从动轮固定”的内啮合副。核心输出随行星架旋转角度变化的啮合刚度曲线含单齿刚度、综合啮合刚度。适用场景健康齿轮的动力学建模、有限元如ANSYS结果验证、齿轮系统设计阶段的刚度特性分析。1.2 程序文件构成程序套件包含3个核心MATLAB脚本功能各有侧重可单独运行或配合使用具体分工如下表所示文件名核心功能关键输出适用场景innerplanetarygear.m计算内啮合副的综合时变啮合刚度区分单齿/双齿啮合区行星轮旋转角度-综合刚度曲线含赫兹接触刚度整体啮合刚度特性分析、系统动力学建模输入single_tooth.m拆分输出单齿对刚度前齿对i、后齿对j含主动轮/从动轮单独刚度行星轮旋转角度-单齿对刚度曲线主动轮/从动轮分别输出单齿刚度验证、啮合区刚度贡献分析singletoothR1.m优化单齿刚度计算逻辑聚焦单个齿的独立刚度剔除冗余啮合对行星轮旋转角度-单齿刚度曲线行星轮/内齿圈明确标注单齿刚度精准验证、有限元对比实验二、关键理论基础与假设2.1 核心计算方法势能法程序以势能守恒原理为核心将齿轮啮合过程中的弹性变形能分解为4个部分通过“总势能各部分势能之和”推导刚度刚度为势能对变形的二阶导数具体分解如下弯曲变形能轮齿受载时产生的弯曲形变对应的能量需考虑齿根过渡曲线与渐开线的连续几何特性。剪切变形能轮齿剪切形变对应的能量采用1.2修正系数符合齿轮力学经典假设。压缩变形能轮齿接触区域压缩形变对应的能量基于横截面面积与弹性模量计算。赫兹接触变形能两齿面接触点的局部变形能量采用赫兹接触理论公式计算接触刚度。2.2 关键假设与边界条件程序通过合理假设简化计算复杂度同时保证工程精度核心假设如下内齿轮基体刚性不考虑内齿圈基体的弹性变形因内齿轮基体刚度成熟计算公式缺失且工程中内齿圈通常固定且刚度远大于轮齿。运动约束内齿圈固定不动行星架带动行星轮绕内齿圈旋转符合行星轮系典型运动模式。啮合相位差通过等分行星轮旋转角度、区分啮合周期内的齿对接触顺序间接体现各啮合副的相位差异。健康齿假设不考虑轮齿裂纹、磨损、点蚀等故障仅计算健康轮齿的啮合刚度。三、程序核心特性与优势3.1 精准的几何建模精确渐开线齿形通过渐开线方程基圆、压力角、展角关联计算齿面坐标同时考虑齿根过渡曲线基于齿条刀具加工原理计算刀顶圆角对应的过渡曲线坐标避免传统“简化齿形”带来的刚度计算误差。齿轮变位适配支持变位齿轮计算通过输入变位系数xp为行星轮变位系数xg为内齿圈变位系数由程序自动推导适配非标准齿轮如为避免根切、提高承载能力设计的变位齿轮的刚度计算。3.2 全面的物理特性考虑程序在参数设置中涵盖齿轮传动的关键物理与几何参数确保计算的工程实用性核心参数如下表所示参数类别参数符号物理意义示例值程序默认齿轮几何参数Zs/Zp/Zr太阳轮/行星轮/内齿圈齿数21/33/87| |mn| 法向模数 | 2.5mm程序中单位为m || |alfa0| 标准齿形角 | 20°程序中转换为弧度 || |hax/cx| 齿顶高系数/顶隙系数 | 1.0/0.25 || 材料参数 |E| 弹性模量 | 2.068×10¹¹Pa钢 |采用势能法编写的行星齿轮内啮合齿轮副时变啮合刚度程序(健康齿)程序中考虑了精确的渐开线齿形以及齿轮变位内齿圈固定行星架旋转同时考虑了考虑各啮合齿轮副之间的相位差可提供相位差计算小程序 假设内齿轮基体为刚性| |v| 泊松比 | 0.3钢 || |G| 剪切模量 | 由E和v自动推导GE/[2(1v)] || 运动参数 |ns| 太阳轮转速 | 1500r/min || |Np| 行星轮个数 | 4个 |3.3 啮合过程精细化模拟啮合区划分自动区分“双齿啮合区”与“单齿啮合区”- 双齿啮合区前齿对即将退出啮合标记为i与后齿对刚进入啮合标记为j同时接触刚度为两齿对刚度叠加。- 单齿啮合区仅后齿对j接触刚度为单齿对刚度。相位差体现通过n_div360*zp每转等分刻度数将行星轮旋转一周划分为等角度步长每个步长对应不同齿对的接触状态间接体现多行星轮啮合时的相位差异。四、程序结构与关键模块解析4.1 通用程序框架以inner_planetary_gear.m为例所有程序均遵循“参数初始化→几何/运动参数计算→刚度积分求解→结果绘图”的逻辑框架核心模块分工明确便于修改与扩展模块1参数初始化齿轮基本参数定义齿轮几何参数齿数、模数、齿形角、材料参数弹性模量、泊松比、运动参数转速、行星轮个数代码片段如下matlab% 齿轮基本参数hax1.0; cx0.25; % 齿顶高系数、顶隙系数Zs21; Zp33; Zr87; % 太阳轮/行星轮/内齿圈齿数mn2.5/1000; L30/1000; % 模数m、齿宽malfa0deg2rad(20); % 标准齿形角radE2.068e11; v0.3; GE/(1v)/2; % 弹性模量、泊松比、剪切模量模块2几何与运动参数计算基于齿轮啮合原理自动推导分度圆半径、基圆半径、实际啮合角、变位系数、齿顶/根圆半径等关键参数示例如下matlabrpmnzp/2; rgmnzg/2; % 行星轮/内齿圈分度圆半径alfaacos(d0cos(alfa0)/d); % 实际啮合角考虑中心距变动rbprpcos(alfa0); rbgrgcos(alfa0); % 基圆半径raprp(haxxpdelty)mn; ragrg-(hax-xgdelty)*mn; % 齿顶圆半径模块3啮合起始/终止位置计算确定啮合起始点P1与分离点P2的坐标与角度参数明确啮合长度与重合度CR代码逻辑如下matlabPP1rbg(tan(alfa)-tan(alfaag)); % 节点P到起始点P1的距离PP2rbp(tan(alfaap)-tan(alfa)); % 节点P到分离点P2的距离P1P2PP1PP2; CRP1P2/Pb; % 实际啮合长度、重合度模块4刚度积分求解核心模块通过梯形数值积分法trapz函数计算轮齿的弯曲、剪切、压缩刚度均为刚度的倒数积分再求总刚度同时叠加赫兹接触刚度示例如下matlab% 弯曲刚度积分渐开线部分kbj2ptrapz(tau,fkbj2p); % 渐开线部分弯曲刚度倒数积分kbj1ptrapz(gamma,fkbj1p); % 过渡曲线部分弯曲刚度倒数积分kbjpkbj1pkbj2p; % 总弯曲刚度倒数% 综合刚度计算串联刚度叠加kt(m)1./(1/khkbjpksjpkajpkbjgksjgkajg); % 单齿啮合区总刚度模块5结果可视化绘制“行星轮旋转角度-刚度”曲线标注坐标轴、图例支持导出高清图符合学术绘图规范代码如下matlabtheta1thetap/pi180; % 角度转换rad→°plot(theta1,kt10^(-8),LineWidth,1.5); grid off;xlabel(Rotation angle of planet gear \thetaP(°),FontName,Times New Roman);ylabel(TVMS kt (\times10^8 N/m),FontName,Times New Roman);4.2 single_tooth_R1.m的优化点作为修改后的单齿刚度计算程序singletoothR1.m在通用框架基础上做了以下优化更适配单齿刚度验证场景简化啮合对逻辑剔除“前齿对i”冗余计算仅聚焦单个啮合对j的单齿刚度避免多齿对干扰。明确刚度输出标注图例直接区分“行星轮单齿刚度planet gear”与“内齿圈单齿刚度ring gear”且关闭图例边框set(h, Box, off)提升图面简洁度。调整啮合周期范围基于重合度CR确定啮合周期thetap[0:2pi/zp/1000:2pi/zp*CR]确保覆盖完整啮合过程无多余角度数据。五、程序使用说明5.1 环境要求MATLAB版本兼容R2016b及以上版本需支持deg2rad、trapz等基础函数。运行依赖无额外工具箱依赖仅需MATLAB基础环境。5.2 参数修改指南用户可根据实际齿轮系统参数在程序“参数初始化”模块修改以下关键参数其余参数由程序自动推导无需手动修改参数符号修改说明注意事项Zs/Zp/Zr太阳轮/行星轮/内齿圈齿数需满足行星轮系传动比关系ZrZs2*Zp默认参数已满足mn法向模数单位m默认2.5mm需转换为0.0025m需与实际齿轮模数一致单位转换需准确L齿宽单位m默认30mm转换为0.03m齿宽直接影响刚度大小需准确输入xp行星轮变位系数默认0即标准齿轮若为变位齿轮需输入设计的变位系数ns太阳轮转速单位r/min默认1500影响啮合频率计算需与实际工况一致5.3 结果解读程序运行后输出的核心结果为“旋转角度-刚度”曲线需结合啮合区特性解读双齿啮合区刚度值较高且相对平稳两齿对刚度叠加对应曲线的“平台段”。单齿啮合区刚度值较低且随角度变化单齿对刚度受齿面接触位置影响对应曲线的“波动段”。单齿刚度验证singletoothR1.m输出的行星轮/内齿圈单齿刚度可与有限元如ANSYS计算结果对比验证势能法的准确性程序作者提示因假设与简化结果可能存在小幅差异需结合具体工况分析。六、程序局限性与注意事项内齿轮基体刚性假设程序未考虑内齿圈基体的弹性变形因缺乏成熟计算公式若需更精准结果需结合有限元补充基体刚度修正。健康齿假设仅适用于无故障的健康齿轮若需计算裂纹、磨损等故障状态的刚度需扩展故障齿形建模模块。有限元对比提示程序作者明确说明“不保证与有限元结果完全吻合”影响因素包括基体刚度简化、赫兹接触假设、积分步长等对比时需考虑这些简化带来的误差。相位差体现方式程序通过等分旋转角度间接体现相位差未直接计算多行星轮的相位角若需精准相位差分析需在“运动参数计算”模块补充相位角公式。七、参考文献与理论依据程序开发严格参考齿轮力学与行星轮系领域的经典文献确保理论正确性核心参考文献如下Xihui Liang, Ming J. Zuo.Analytically evaluating the influence of crack on the mesh stiffness of a planetary gear set(2014)马辉, 左明键等.Time-varying mesh stiffness calculation of cracked spur gears势能法计算时变啮合刚度的经典文献罗阳, 马辉等. 行星齿轮系统动力学建模与刚度计算相关文献程序备注中明确提及八、总结本程序套件是行星齿轮内啮合副时变啮合刚度计算的高效工具具有几何建模精准、物理特性全面、使用便捷的优势可满足健康齿轮动力学分析的基础需求。singletoothR1.m作为优化后的单齿刚度计算程序更适用于刚度验证与对比实验innerplanetarygear.m则适用于系统级动力学建模的综合刚度输入。用户在使用时需注意参数单位转换的准确性、实际工况与程序假设的匹配度若需扩展至故障齿轮或更复杂场景可基于现有框架补充故障建模、基体刚度修正等模块进一步提升程序的工程适用性。