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

概述

本代码实现了一个完整的六自由度直齿轮传动系统弯扭耦合动力学模型，采用集中质量法建立数学模型，并通过MATLAB的ODE45求解器进行数值计算。该模型综合考虑了时变啮合刚度、齿侧间隙非线性、误差激励等关键因素，能够全面分析齿轮系统的动态特性。

核心数学模型

系统动力学方程

根据集中质量法，系统建立以下12个一阶微分方程：

主动轮（小齿轮）方程：

• 水平方向：mp·ẍp + cpx·ẋp + kpx·xp = -Fm·cos(α)
• 垂直方向：mp·ÿp + cpy·ẏp + kpy·yp = Fm·sin(α) - mp·g
• 扭转方向：Jp·θ̈p = Tp - Fm·Rp

从动轮（大齿轮）方程：

• 水平方向：mg·ẍg + cgx·ẋg + kgx·xg = Fm·cos(α)
• 垂直方向：mg·ÿg + cgy·ẏg + kgy·yg = -Fm·sin(α) - mg·g
• 扭转方向：Jg·θ̈g = -Tg + Fm·Rg

关键物理模型

1. 时变啮合刚度

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

采用傅里叶级数展开模拟时变刚度：

km = (a₀/2 + Σ[aₙ·cos(n·(ωm·t + π·(2-εα)))])·B

其中单双齿啮合刚度通过重合度εα进行加权计算。

2. 齿侧间隙非线性

采用分段线性函数描述：

       { δ - b,  δ ≥ b
f(δ) = { 0,     |δ| < b
       { δ + b,  δ ≤ -b

3. 动态啮合力

Fm = cm·δ̇ + km·f(δ)

其中δ为啮合变形量，包含齿轮相对位移和误差激励。

代码结构

1. 主求解文件 (`Solve.m`)

• 参数初始化（齿数、转速、间隙等）
• 调用ODE45求解动力学方程
• 生成全面的结果分析图表

2. 分岔分析文件 (`Bifurcation.m`)

• 齿侧间隙参数扫描（0-100μm）
• 采用庞加莱截面法提取周期解
• 生成分岔图揭示系统非线性行为

3. 动力学方程文件 (`Straight_Gear.m`)

• 实现完整的12自由度动力学方程
• 计算时变啮合刚度、啮合阻尼、误差激励
• 处理齿侧间隙非线性

输出结果分析

时域分析

• 振动位移：各方向位移响应时程曲线
• 振动速度：系统运动状态变化率
• 振动加速度：通过数值微分计算，反映动态载荷

频域分析

• 频谱图：通过FFT变换识别系统主要频率成分
• 可观察到啮合频率及其谐波分量

非线性特性分析

• 相图：位移-速度相平面轨迹，揭示系统稳定性
• 庞加莱截面：每隔激励周期采样，区分周期/混沌运动
• 分岔图：展示系统随齿侧间隙变化的动力学行为转变

物理意义与工程应用

1. 参数影响分析：通过改变齿侧间隙，研究其对系统稳定性的影响
2. 故障诊断：异常的振动特征可指示齿轮磨损、断齿等故障
3. 设计优化：为齿轮参数选择提供理论依据，避免共振和混沌
4. 动态载荷预测：为强度设计和寿命评估提供载荷数据

技术特点

• 高精度数值计算：采用变步长ODE45算法，保证计算精度
• 完整非线性考虑：同时包含几何非线性和参数激励
• 多维度分析：从时域、频域、相空间多角度全面分析
• 工程实用性：输出结果可直接用于工程设计和故障诊断

该代码为齿轮系统动力学研究提供了完整的仿真平台，能够深入揭示齿轮传动系统的复杂动力学行为，对高性能齿轮传动设计具有重要指导意义。