MATLAB六自由度齿轮弯扭耦合动力学代码（含时变啮合刚度、齿侧间隙及集中质量法建模的数值计算分析）
MATLAB六自由度齿轮弯扭耦合动力学代码（考虑时变啮合刚度、齿侧间隙），根据集中质量法建模（含数学方程建立和公式推导）并在MATLAB中采用ODE45进行数值计算。 输出齿轮水平和竖直方向的振动位移、振动速度、振动加速度、轮齿间动态啮合力、相图、庞加莱图、分岔图、频谱图。
一、程序概述
本MATLAB程序基于集中质量法构建了六自由度齿轮弯扭扭耦合动力学模型，重点用于分析齿轮传动系统在考虑时变啮合刚度和齿侧间隙等非线性因素影响下的动态特性。程序通过数值计算方法求解复杂的动力学方程，能够输出齿轮系统的振动响应、动态啮合力及多种种非线性动力学特征图谱，为齿轮系统的设计优化、故障诊断和动态性能评估提供理论依据和仿真工具。
二、模型构建理论基础
1. 动力学模型自由度定义
程序采用集中质量法将齿轮系统简化为六自由度动力学模型，具体自由度包括：
- 主动轮：水平方向位移xₚ、竖直方向位移yₚ、扭转角位移θₚ
- 从动轮：水平方向位移xᵍ、竖直方向位移yᵍ、扭转角位移θᵍ
对应每个位移自由度，包含相应的速度变量，形成12阶状态向量用于微分方程求解。
2. 关键力学模型与公式
（1）啮合变形计算
齿轮啮合点的综合变形量δ考虑了平移、转动和啮合误差的影响：
δ = (xₚ - xᵍ)cosα + (yₚ - yᵍ)sinα + (Rₚθₚ - Rᵍθᵍ) - e(t)
其中：α为压力角，Rₚ、Rᵍ分别为主动轮和从动轮的分度圆半径，e(t)为啮合误差函数（包含均值和余弦波动分量）
（2）齿侧间隙非线性模型
考虑齿侧间隙b的分段非线性特性，啮合变形有效量f(δ)定义为：
- 当δ ≥ b时，f(δ) = δ - b
- 当δ ≤ -b时，f(δ) = δ + b
- 当|δ| < b时，f(δ) = 0（轮齿脱离接触状态）
（3）时变啮合刚度模型
啮合刚度kₘ(t)采用傅里叶级数形式描述其周期性变化：
kₘ(t) = [a₀/2 + Σaₙcos(nωₘt + φₙ)] × B
MATLAB六自由度齿轮弯扭耦合动力学代码（考虑时变啮合刚度、齿侧间隙），根据集中质量法建模（含数学方程建立和公式推导）并在MATLAB中采用ODE45进行数值计算。 输出齿轮水平和竖直方向的振动位移、振动速度、振动加速度、轮齿间动态啮合力、相图、庞加莱图、分岔图、频谱图。
其中：a₀、aₙ为傅里叶系数，ωₘ为啮合频率，B为齿宽，同时考虑单齿啮合与双齿啮合阶段的刚度差异
（4）动态啮合力计算
轮齿间动态啮合力Fₘ由阻尼力和弹性力组成：
Fₘ = cₘ×dδ/dt + kₘ(t)×f(δ)
其中：cₘ为啮合阻尼系数，dδ/dt为啮合变形速率
（5）运动微分方程
系统动力学方程基于牛顿第二定律和动量矩定理建立，以主动轮水平方向运动为例：
mₚ×ẍₚ = -Fₘcosα - cₚₓ×ẋₚ - kₚₓ×xₚ
其他五个自由度的运动方程形式类似，分别考虑各自方向的惯性力、阻尼力、弹性力和啮合力分量。
三、程序结构与功能模块
1. 核心文件组成
| 文件名 | 主要功能 |
|---|---|
| Straight_Gear.m | 定义系统物理参数，构建动力学微分方程，计算时变啮合刚度和动态啮合力 |
| Solve.m | 调用ODE45数值求解器求解微分方程，处理计算结果并生成各类动力学特性图表 |
| Bifurcation.m | 分析系统随关键参数（齿侧间隙）变化的非线性动力学行为，绘制分岔图 |
2. 主要功能模块
（1）参数设置模块
在Straight_Gear.m中定义齿轮系统的基本参数，包括：
- 齿轮几何参数：齿数、模数、压力角、分度圆半径等
- 物理参数：齿轮质量、转动惯量、支撑刚度和阻尼系数
- 工况参数：输入转速、齿侧间隙、啮合误差幅值等
（2）动力学方程求解模块
在Solve.m中实现：
- 采用ODE45数值方法求解12阶非线性微分方程组
- 设置积分时间范围和初始条件（微小扰动避免零初始值奇点）
- 提取并处理系统响应数据，分离暂态响应和稳态响应
（3）结果分析与可视化模块
生成多种动力学特性图表：
- 时域响应：水平/竖直方向振动位移、速度、加速度随时间变化曲线
- 动态啮合力：轮齿间接触力的时域变化曲线
- 相图：位移-速度平面轨迹，反映系统运动稳定性
- 庞加莱图：特定周期时刻的状态点分布，识别运动周期特性
- 频谱图：通过FFT变换获得振动信号的频率成分和幅值分布
- 分岔图：随齿侧间隙变化的系统响应分岔特性，揭示非线性动力学行为转变
四、关键计算流程
- 参数初始化：设定齿轮基本参数和仿真工况，计算啮合频率、转动角速度等派生参数
- 方程构建：在Straight_Gear.m中定义状态向量导数函数，实现动力学方程的数值表达
- 数值求解：调用ODE45求解器对微分方程组进行积分计算，获得系统状态随时间的变化
- 结果后处理：
- 振动加速度通过速度信号的数值微分获得
- 稳态响应提取：选取足够长时间后的响应数据（如t>0.1s）
- 频谱分析：对稳态振动信号进行快速傅里叶变换
- 相图与庞加莱图：提取特定状态变量组合并绘图 - 分岔分析：在Bifurcation.m中通过改变齿侧间隙参数，重复求解过程，绘制分岔图谱
五、程序应用与扩展
1. 主要应用场景
- 齿轮系统动态特性预测与分析
- 齿侧间隙、时变刚度等参数对系统动态性能的影响研究
- 齿轮非线性振动与动力学行为分析
- 齿轮故障机理研究与故障特征提取
2. 程序扩展建议
- 增加更多故障模型（如齿面磨损、断齿等）
- 引入更复杂的啮合刚度计算模型
- 扩展为齿轮传动系统多体动力学模型
- 增加参数敏感性分析和优化模块
本程序通过系统化的建模方法和完整的求解流程，为齿轮弯扭耦合动力学分析提供了可靠的数值仿真工具，可有效支持齿轮传动系统的设计改进和性能评估。