matlab:行星齿轮动力学集中质量参数模型行星齿轮箱的振动特性总是让人着迷——那种行星轮、太阳轮和齿圈之间复杂的相互作用像极了机械系统里的芭蕾舞剧。今天咱们用MATLAB搞点实战试试用集中质量模型来抓取这出舞蹈的动力学特征。先给模型定个基调每个齿轮都看作刚性圆盘用弹簧阻尼单元连接。假设有3个行星轮工程常见配置系统总自由度包括太阳轮旋转、行星架旋转、行星轮自转。敲黑板重点参数来了% 系统基本参数 Ns 20; % 太阳轮齿数 Np 30; % 行星轮齿数 Nr 80; % 齿圈齿数 m_p 0.8; % 行星轮质量/kg k_mesh 5e7; % 啮合刚度/(N/m) c_mesh 500; % 啮合阻尼/(N·s/m)这里有个坑要注意啮合刚度不是定值实际应该用时变刚度但咱们先简化为平均值方便计算。接下来构建运动方程核心在于处理行星轮的相位差% 行星轮位置角度均布 theta_pos linspace(0, 2*pi*(1-1/3), 3); % 构建质量矩阵M M diag([Js, Jc, Jp*ones(1,3)]); % Js太阳轮惯量, Jc行星架惯量 % 刚度矩阵K初始化 K zeros(5,5); for i 1:3 phase theta_pos(i); % 太阳轮-行星轮啮合刚度 K(1, 2i) -k_mesh * sin(phase); % 行星轮-齿圈啮合项 K(2i, 2i) k_mesh*(1 (Nr/Np)^2); end这段代码暗藏玄机刚度矩阵的非对角项反映了太阳轮与行星轮的动态耦合行星轮位置的正弦项体现了空间相位的影响。齿圈固定时其作用转化为对行星轮的等效刚度增强。matlab:行星齿轮动力学集中质量参数模型上主菜——求解微分方程组。用ode45解算前需要把二阶方程转换成一阶function dx gear_ode(t, x, M, K, C) % 状态变量分解 q x(1:5); % 位移 dq x(6:10); % 速度 % 构造加速度 ddq M \ (-C*dq - K*q external_force(t)); dx [dq; ddq]; end注意这里的external_force函数可以模拟负载波动比如设置成sin函数来模拟扭矩波动。跑仿真时初始条件别乱设比如给太阳轮一个微小扰动x0 [zeros(5,1); [0; 0.1; zeros(3,1)]]; [t, X] ode45((t,x)gear_ode(t,x,M,K,C), [0 0.5], x0);后处理阶段才是精华所在。看时域响应可能平平无奇但频谱分析会暴露玄机% 太阳轮振动频谱 [pxx, f] pwelch(X(:,1), 512, [], [], 1/0.001); plot(f, 10*log10(pxx)) xlabel(频率 (Hz)) ylabel(功率/dB)通常在频谱上能看到行星轮过轴频率行星架转频×行星轮数量及其边带还有啮合频率的调制现象。比如当啮合频率为500Hz时出现以该频率为中心的边带群间距正是行星架转频。仿真时遇到过数值发散试试调整阻尼系数c_mesh。刚度矩阵的条件数过大时可能需要隐式算法如ode15s。这个模型虽然简化了轮齿接触细节但足以捕捉系统的主要共振特征比多体动力学仿真快两个数量级特别适合参数优化初期使用。下次试试引入时变刚度和齿侧间隙那会让频谱出现超有趣的次谐波成分。不过那是另一个故事了... MATLAB中的行星齿轮动力学：集中质量参数模型
matlab:行星齿轮动力学集中质量参数模型行星齿轮箱的振动特性总是让人着迷——那种行星轮、太阳轮和齿圈之间复杂的相互作用像极了机械系统里的芭蕾舞剧。今天咱们用MATLAB搞点实战试试用集中质量模型来抓取这出舞蹈的动力学特征。先给模型定个基调每个齿轮都看作刚性圆盘用弹簧阻尼单元连接。假设有3个行星轮工程常见配置系统总自由度包括太阳轮旋转、行星架旋转、行星轮自转。敲黑板重点参数来了% 系统基本参数 Ns 20; % 太阳轮齿数 Np 30; % 行星轮齿数 Nr 80; % 齿圈齿数 m_p 0.8; % 行星轮质量/kg k_mesh 5e7; % 啮合刚度/(N/m) c_mesh 500; % 啮合阻尼/(N·s/m)这里有个坑要注意啮合刚度不是定值实际应该用时变刚度但咱们先简化为平均值方便计算。接下来构建运动方程核心在于处理行星轮的相位差% 行星轮位置角度均布 theta_pos linspace(0, 2*pi*(1-1/3), 3); % 构建质量矩阵M M diag([Js, Jc, Jp*ones(1,3)]); % Js太阳轮惯量, Jc行星架惯量 % 刚度矩阵K初始化 K zeros(5,5); for i 1:3 phase theta_pos(i); % 太阳轮-行星轮啮合刚度 K(1, 2i) -k_mesh * sin(phase); % 行星轮-齿圈啮合项 K(2i, 2i) k_mesh*(1 (Nr/Np)^2); end这段代码暗藏玄机刚度矩阵的非对角项反映了太阳轮与行星轮的动态耦合行星轮位置的正弦项体现了空间相位的影响。齿圈固定时其作用转化为对行星轮的等效刚度增强。matlab:行星齿轮动力学集中质量参数模型上主菜——求解微分方程组。用ode45解算前需要把二阶方程转换成一阶function dx gear_ode(t, x, M, K, C) % 状态变量分解 q x(1:5); % 位移 dq x(6:10); % 速度 % 构造加速度 ddq M \ (-C*dq - K*q external_force(t)); dx [dq; ddq]; end注意这里的external_force函数可以模拟负载波动比如设置成sin函数来模拟扭矩波动。跑仿真时初始条件别乱设比如给太阳轮一个微小扰动x0 [zeros(5,1); [0; 0.1; zeros(3,1)]]; [t, X] ode45((t,x)gear_ode(t,x,M,K,C), [0 0.5], x0);后处理阶段才是精华所在。看时域响应可能平平无奇但频谱分析会暴露玄机% 太阳轮振动频谱 [pxx, f] pwelch(X(:,1), 512, [], [], 1/0.001); plot(f, 10*log10(pxx)) xlabel(频率 (Hz)) ylabel(功率/dB)通常在频谱上能看到行星轮过轴频率行星架转频×行星轮数量及其边带还有啮合频率的调制现象。比如当啮合频率为500Hz时出现以该频率为中心的边带群间距正是行星架转频。仿真时遇到过数值发散试试调整阻尼系数c_mesh。刚度矩阵的条件数过大时可能需要隐式算法如ode15s。这个模型虽然简化了轮齿接触细节但足以捕捉系统的主要共振特征比多体动力学仿真快两个数量级特别适合参数优化初期使用。下次试试引入时变刚度和齿侧间隙那会让频谱出现超有趣的次谐波成分。不过那是另一个故事了...