MATLAB六自由度齿轮弯扭耦合动力学代码考虑时变啮合刚度、齿侧间隙根据集中质量法建模含数学方程建立和公式推导并在MATLAB中采用ODE45进行数值计算。 输出齿轮水平和竖直方向的振动位移、振动速度、振动加速度、轮齿间动态啮合力、相图、庞加莱图、分岔图、频谱图。代码概述本MATLAB代码包实现了一个完整的直齿轮传动系统动力学仿真分析工具包含三个主要文件Straight_Gear.m动力学方程、Bifurcation.m分岔分析和Solve.m时域和频域分析。代码采用集中质量法建立了六自由度弯扭耦合动力学模型全面考虑了齿轮传动的关键非线性因素。核心文件功能详解1. Straight_Gear.m - 动力学方程核心功能定义12状态变量的齿轮系统动力学微分方程关键模型要素时变啮合刚度计算matlabksingle 0.8/(0.047230.15551/zp0.25791/zg)1e9; % 单齿刚度kmean ksingle(0.75epsilonalpha0.25); % 平均刚度kdouble (kmean-(2-epsilonalpha)ksingle)/(epsilon_alpha-1); % 双齿刚度通过傅里叶级数展开实现刚度周期性变化齿侧间隙非线性函数matlabif deltabfdeltadelta-b;elseif delta-bfdeltadeltab;elsef_delta0;end采用三段式分段线性函数描述间隙冲击效应动态啮合力计算matlabFm cmdeltat kmf_delta; % 阻尼力 弹性力系统参数小齿轮55齿质量1.53kg转动惯量0.041kg·m²大齿轮75齿质量3.01kg转动惯量0.079kg·m²支撑刚度2.0e8 N/m支撑阻尼1.8e3 N·s/m输入扭矩200 N·m转速1000 rpm2. Bifurcation.m - 参数扫描分岔分析功能研究齿侧间隙对系统动力学行为的影响算法流程设置齿侧间隙扫描范围0-100μm步长1e-6对每个间隙值- 瞬态过程模拟2000个啮合周期- 稳态数据采集100个啮合周期- 庞加莱截面采样每周期512点生成分岔图识别周期运动向混沌运动的转变关键技术采用庞加莱截面法提取周期解通过长时间瞬态模拟确保系统进入稳态使用ODE45求解器保证计算精度3. Solve.m - 综合动力学分析功能在固定参数下进行全面的动力学特性分析MATLAB六自由度齿轮弯扭耦合动力学代码考虑时变啮合刚度、齿侧间隙根据集中质量法建模含数学方程建立和公式推导并在MATLAB中采用ODE45进行数值计算。 输出齿轮水平和竖直方向的振动位移、振动速度、振动加速度、轮齿间动态啮合力、相图、庞加莱图、分岔图、频谱图。分析内容时域分析振动位移xp, yp, xg, yg时程曲线振动速度时程曲线振动加速度通过梯度法数值微分相图分析位移-速度相平面轨迹截取稳态响应段t0.1s庞加莱映射for i1:cycles m1(i,1)y(find(t(0.1Tm*i),1),1); % 周期采样 n1(i,1)y(find(t(0.1Tm*i),1),2); end每啮合周期采样一点形成庞加莱截面频域分析FFT频谱分析频率范围0-1000Hz采样频率100kHz满足采样定理去除直流分量准确识别频率成分数学模型基础运动微分方程系统包含12个一阶微分方程对应6个自由度的位移和速度主动轮水平方向dy(1) y(2); dy(2) (-Fm*cos(alpha) - cpx*y(2) - kpx*y(1))/mp;从动轮扭转方向dy(11) y(12); dy(12) (-Tg Fm*Rg)/Jg;啮合变形计算delta (y(1)-y(7))*cos(alpha) (y(3)-y(9))*sin(alpha) (Rp*y(5)-Rg*y(11)) - et;综合考虑了平移振动、扭转振动和制造误差输出结果类型分岔图齿侧间隙-振动位移关系识别动力学状态转变时域波形各自由度振动响应时间历程相图稳态运动的相轨迹反映系统稳定性庞加莱图区分周期、拟周期和混沌运动频谱图频率成分分析识别啮合频率及其谐波工程应用价值本代码适用于齿轮系统动态设计参数优化齿侧间隙对振动特性影响评估非线性动力学行为研究故障诊断特征提取齿轮箱减振降噪设计通过改变系统参数转速、扭矩、间隙等可全面分析直齿轮传动系统的动力学特性为工程实践提供理论指导。