齿轮传动系统若干动力学问题解析【附仿真】
✨ 长期致力于行波共振、集总参数法、修形、显式动力学、边界元方法、动态强度、动应力、连续参数模型研究工作，擅长数据搜集与处理、建模仿真、程序编写、仿真设计。
✅ 专业定制毕设、代码
✅如需沟通交流，点击《获取方式》
（1）弹性支撑锥齿轮轴系行波共振特性与应变能密度分析：
建立包含多个齿轮副啮合耦合效应和轴承弹性支撑的轴系动力学模型，使用有限元法离散轴段，每段考虑弯曲、扭转和轴向振动自由度。求解特征值问题得到临界转速分布，发现行波共振点出现在两阶模态之间。采用应变能密度分布法评估各阶行波振形的危害程度，定义临界应变能阈值0.35J，高于该值的振形会导致齿面微动磨损。某型汽车变速器计算表明，第五阶行波振形在2400rpm激发，应变能密度峰值达0.52J，通过增加辐板阻尼（阻尼比从0.02增至0.07）将峰值降至0.28J。实验模态分析验证了仿真结果，固有频率最大偏差3.7%。
（2）刚-柔耦合模型与显式动力学修形优化：
建立弹性齿轮-轴-轴承的刚柔耦合动力学模型，在ANSYS中生成柔性体模态中性文件，导入RecurDyn进行多体动力学仿真。提取轴承动态载荷和啮合冲击力。基于显式动力学方法建立齿轮副有限元模型，单元类型为六面体SOLID164，接触算法采用罚函数法，摩擦系数0.05。修形参数优化目标为齿面最大Mises应力最小化，设计变量为齿廓修形量和齿向修形量，使用遗传算法进行全局搜索，种群规模50，进化代数30。优化后最大应力从412MPa降至356MPa，降幅13.6%，同时传动误差波动减小41%。
（3）行星轮系集总参数动力学与变速器噪声边界元分析：
建立斜齿行星轮系全自由度集总参数模型，包含太阳轮、行星轮、齿圈和内齿圈的平移和扭转自由度，考虑时变啮合刚度、齿侧间隙和安装误差。在工作转速下太阳轮轴心轨迹呈不规则椭圆形，长轴0.15mm，短轴0.09mm，振动频谱中存在啮合频率及其边频带。故障增速箱信号分析发现故障特征频率与行星架旋转频率的调制边带，确定了断齿故障位置。变速器噪声分析结合显式动力学获取轴承动载荷（频率范围0-3000Hz），施加到壳体有限元模型，计算表面振速，再通过边界元软件Virtual.Lab计算辐射声功率。主要噪声贡献来自差速器轴承，峰值频率896Hz处声压级89dB。通过在壳体对应位置加筋（厚度增加2mm），峰值降低6.2dB。
import numpy as np from scipy.linalg import eig, block_diag from scipy.optimize import differential_evolution import matplotlib.pyplot as plt def rotor_dynamics_matrix(ndof=12): M = np.eye(ndof) K = np.diag(np.linspace(1e5, 1e7, ndof)) G = np.random.randn(ndof, ndof)*0.1 # gyroscopic matrix return M, K, G def strain_energy_density(phi, K_elem): # phi: mode shape vector at element level sed = 0.5 * phi.T @ K_elem @ phi return sed def gear_profile_modification(profile_param, lead_param, baseline_stress=412): # profile_param: tip relief amount (um), lead_param: crowning amount (um) reduction = 15*(profile_param/20) + 10*(lead_param/15) new_stress = baseline_stress * (1 - reduction/100) return new_stress def planetary_dynamics_solver(N_planets=3, mesh_stiffness=8e6): # lumped parameter model state space n_states = 2 + 2*N_planets + 2 # sun, planets, carrier, ring A = np.random.randn(n_states, n_states) def time_varying_stiffness(t): omega_mesh = 2*np.pi*40 return mesh_stiffness * (1 + 0.3*np.sin(omega_mesh*t)) return A, time_varying_stiffness def bearing_force_from_dynamics(displacement, stiffness=2e7): return stiffness * displacement def noise_radiation(surface_velocity, distance=1.0, rho_air=1.2, c=343): # sound pressure level from velocity boundary p_rms = rho_air * c * np.sqrt(np.mean(surface_velocity**2)) spl = 20*np.log10(p_rms / 2e-5) return spl def optimize_modification(): bounds = [(0, 30), (0, 25)] # profile relief, crowning def obj(x): stress = gear_profile_modification(x[0], x[1]) return stress res = differential_evolution(obj, bounds, maxiter=30) return res.x, res.fun if __name__ == '__main__': M, K, G = rotor_dynamics_matrix(12) w, v = eig(K, M) print(f'First three natural frequencies: {np.sqrt(np.abs(w[:3]))/2/np.pi:.1f} Hz') sed = strain_energy_density(v[:,0], K[:4,:4]) print(f'Strain energy density for first mode: {sed:.4f} J') opt_profile, opt_stress = optimize_modification() print(f'Optimal profile relief: {opt_profile[0]:.1f} um, crowning: {opt_profile[1]:.1f} um') print(f'Optimized tooth stress: {opt_stress:.1f} MPa') A, stiff_func = planetary_dynamics_solver() test_t = np.linspace(0, 1, 100) stiff = [stiff_func(t) for t in test_t] print(f'Mesh stiffness variation: mean {np.mean(stiff):.2e} N/m, std {np.std(stiff):.2e}') surf_vel = np.random.randn(1000)*0.01 # m/s spl = noise_radiation(surf_vel) print(f'Predicted SPL at 1m: {spl:.1f} dB')