基于MATLAB Simulink的定轴齿轮与行星齿轮仿真模型

一、定轴齿轮传动系统建模

1. 模型架构设计

核心模块 ：使用Simscape Mechanical库中的Rotational Gear模块构建齿轮副，结合Transfer Function或State-Space模块模拟传动特性。
参数设置 ：
- 传动比：根据齿轮齿数比计算（如主动轮齿数Z1=20，从动轮Z2=60，传动比i=3）。
- 转动惯量 ：齿轮质量转换为等效转动惯量（J=m⋅r2J=m⋅r^2J=m⋅r2）。
- 啮合刚度 ：通过Rotational Damper模块模拟时变刚度（公式：k(t)=k0⋅(1+α⋅sin(ωt))k(t)=k_0⋅(1+α⋅sin(ωt))k(t)=k0⋅(1+α⋅sin(ωt))）。

2. Simulink实现代码

matlab 复制代码

%% 定轴齿轮模型搭建
% 齿轮参数
Z1 = 20; Z2 = 60; % 主动轮/从动轮齿数
i = Z2/Z1;        % 传动比
J_gear1 = 0.1;    % 齿轮1等效转动惯量 (kg·m²)
J_gear2 = 0.4;    % 齿轮2等效转动惯量 (kg·m²)
k_mesh = 1e6;     % 啮合刚度 (N·m/rad)

% 模块搭建
gear1 = simscape.mechanical.rotational.RotationalGear('GearRatio', 1, ...
    'Inertia', J_gear1, 'Damping', 0.1);
gear2 = simscape.mechanical.rotational.RotationalGear('GearRatio', i, ...
    'Inertia', J_gear2, 'Damping', 0.1);
damper = simscape.mechanical.rotational.RotationalDamper('DampingCoefficient', 0.05);

% 连接关系
gear1.Outputs(1).Connect(gear2.Inputs(1));
damper.Inputs(1).Connect(gear1.Outputs(1));

3. 关键仿真设置

输入信号 ：使用Constant Torque模块施加负载扭矩（如50 N·m）。
求解器 ：选择ode45，步长0.001秒。
输出监测 ：通过Scope观察齿轮转速差（Δω）和扭矩波动。

二、行星齿轮传动系统建模

1... Simulink实现代码

matlab 复制代码

%% 行星齿轮模型搭建（集中质量法）
% 参数定义
m_s = 1.2; m_p = 0.8; m_r = 3.0; % 太阳轮/行星轮/齿圈质量 (kg)
k_mesh = 1.5e6; c_mesh = 1500;   % 啮合刚度与阻尼
r_s = 0.1; r_p = 0.08; r_r = 0.3; % 节圆半径 (m)

% 状态变量定义（θ_s, θ_p, θ_r, ω_s, ω_p, ω_r）
states = @(t,y) [y(4); y(5); y(6); ...
    (k_mesh*(y(3)-y(1)) - c_mesh*y(4))/m_s; ...
    (-k_mesh*(y(2)-y(3)) - c_mesh*y(5))/m_p; ...
    (k_mesh*(y(1)-y(2)) - c_mesh*y(6))/m_r];

% ODE45求解器设置
tspan = [0 10]; % 仿真时间
y0 = [0; 0; 0; 100; 0; 0]; % 初始条件（角度/角速度）

% 仿真求解
[t,y] = ode45(states, tspan, y0);

% 结果可视化
figure;
subplot(2,1,1);
plot(t, y(:,1:3)*180/pi); % 转角曲线
xlabel('时间 (s)'); ylabel('转角 (°)');
legend('太阳轮', '行星轮', '齿圈');

subplot(2,1,2);
plot(t, y(:,4:6)); % 角速度曲线
xlabel('时间 (s)'); ylabel('角速度 (rad/s)');
legend('ω_s', 'ω_p', 'ω_r');

2. 动态特性分析

收敛性验证：通过相图（θ vs. ω）检查自由度是否收敛（参考）。
频谱分析：使用FFT分析振动信号，识别啮合频率及其倍频（公式：fmesh=z⋅n/60，z为齿数）。

三、多自由度耦合仿真优化

1. 参数敏感性分析

刚度影响：调整kmesh观察共振频率变化（刚度增大→共振频率升高）。
阻尼优化 ：通过fminsearch寻找最优cmesh抑制振动幅值。

2. 故障模拟扩展

齿面剥落：在啮合力矩中叠加周期性冲击（代码示例）：

matlab 复制代码

% 齿面剥落故障模型
fault_gain = 0.5; % 故障程度 (0-1)
T_fault = fault_gain * 1000 * sin(2*pi*5*t); % 5Hz冲击
y(4) = y(4) + T_fault;

3. 控制策略集成

PID调速：在Simulink中添加PID控制器模块，调节行星架转速：

matlab 复制代码

% PID控制器参数
Kp = 0.8; Ki = 0.1; Kd = 0.05;
controller = pid(Kp, Ki, Kd);
output = controller(y(6) - ref_speed);

参考代码定轴齿轮和行星齿轮的仿真模型 www.youwenfan.com/contentcsq/52490.html

四、仿真结果验证与调试

1. 验证方法

理论对比：将仿真结果与集中质量法解析解对比（误差应<5%）。
实验数据匹配：导入实际齿轮箱振动数据，调整刚度/阻尼参数使仿真曲线吻合。

2. 常见问题调试

发散问题：检查刚度矩阵是否正定，或降低积分步长。
非物理振荡：引入瑞利阻尼（c=2ξωm）抑制高频振荡。

注：实际建模时需根据具体齿轮参数（如模数、压力角）调整刚度计算公式，并通过实验数据标定模型精度。