MATLAB六自由度齿轮弯扭耦合动力学代码（考虑时变啮合刚度、齿侧间隙），根据集中质量法建模（含数学方程建立和公式推导）并在MATLAB中采用ODE45进行数值计算。 输出齿轮水平和竖直方向的振动位移、振动速度、振动加速度、轮齿间动态啮合力、相图、庞加莱图、分岔图、频谱图。

代码概述

本MATLAB代码包实现了一个完整的直齿轮传动系统动力学仿真分析工具，包含三个主要文件：Straight_Gear.m（动力学方程）、Bifurcation.m（分岔分析）和Solve.m（时域和频域分析）。代码采用集中质量法建立了六自由度弯扭耦合动力学模型，全面考虑了齿轮传动的关键非线性因素。

核心文件功能详解

1. Straight_Gear.m - 动力学方程核心

功能：定义12状态变量的齿轮系统动力学微分方程

关键模型要素：

• 时变啮合刚度计算：
  matlab
ksingle = 0.8/(0.04723+0.15551/zp+0.25791/zg)1e9; % 单齿刚度
kmean = ksingle(0.75epsilonalpha+0.25); % 平均刚度
kdouble = (kmean-(2-epsilonalpha)ksingle)/(epsilon_alpha-1); % 双齿刚度

  通过傅里叶级数展开实现刚度周期性变化

• 齿侧间隙非线性函数：
  matlab
if delta>=b
fdelta=delta-b;
elseif delta<=-b
fdelta=delta+b;
else
f_delta=0;
end

  采用三段式分段线性函数描述间隙冲击效应

• 动态啮合力计算：
  matlab
Fm = 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 rpm

2. Bifurcation.m - 参数扫描分岔分析

功能：研究齿侧间隙对系统动力学行为的影响

算法流程：

1. 设置齿侧间隙扫描范围：0-100μm，步长1e-6
2. 对每个间隙值：
   - 瞬态过程模拟（2000个啮合周期）
   - 稳态数据采集（100个啮合周期）
   - 庞加莱截面采样（每周期512点）
3. 生成分岔图，识别周期运动向混沌运动的转变

关键技术：

• 采用庞加莱截面法提取周期解
• 通过长时间瞬态模拟确保系统进入稳态
• 使用ODE45求解器保证计算精度

3. Solve.m - 综合动力学分析

功能：在固定参数下进行全面的动力学特性分析

MATLAB六自由度齿轮弯扭耦合动力学代码（考虑时变啮合刚度、齿侧间隙），根据集中质量法建模（含数学方程建立和公式推导）并在MATLAB中采用ODE45进行数值计算。 输出齿轮水平和竖直方向的振动位移、振动速度、振动加速度、轮齿间动态啮合力、相图、庞加莱图、分岔图、频谱图。

分析内容：

时域分析：

• 振动位移（xp, yp, xg, yg）时程曲线
• 振动速度时程曲线
• 振动加速度（通过梯度法数值微分）

相图分析：

• 位移-速度相平面轨迹
• 截取稳态响应段（t>0.1s）

庞加莱映射：

for i=1:cycles
    m1(i,1)=y(find(t>(0.1+Tm*i),1),1);  % 周期采样
    n1(i,1)=y(find(t>(0.1+Tm*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;

综合考虑了平移振动、扭转振动和制造误差

输出结果类型

1. 分岔图：齿侧间隙-振动位移关系，识别动力学状态转变
2. 时域波形：各自由度振动响应时间历程
3. 相图：稳态运动的相轨迹，反映系统稳定性
4. 庞加莱图：区分周期、拟周期和混沌运动
5. 频谱图：频率成分分析，识别啮合频率及其谐波

工程应用价值

本代码适用于：

• 齿轮系统动态设计参数优化
• 齿侧间隙对振动特性影响评估
• 非线性动力学行为研究
• 故障诊断特征提取
• 齿轮箱减振降噪设计

通过改变系统参数（转速、扭矩、间隙等），可全面分析直齿轮传动系统的动力学特性，为工程实践提供理论指导。