直齿行星齿轮纯扭转非线性动力学建模(含3个或4个行星轮)包括动力学方程推导过程考虑了各啮合齿轮副之间的啮合相位可出相图频谱图分岔图庞加莱映射。行星齿轮系统的非线性动力学分析总带着点机械工程的浪漫——既要数学的严谨又得兼顾物理的直观。今天咱们直接上手搞个3行星轮的直齿系统扭转振动模型手撕动力学方程的同时顺带搞出相图频谱这些可视化效果。先看齿轮副的时变啮合刚度怎么处理。假设太阳轮与行星轮间的啮合刚度为k_sp1.5e8(10.2sin(ωtφ))这里φ就是相位差的关键。当行星轮数量n3时相邻行星轮的相位差Δφ2π/3。这个相位差直接影响着系统共振点的分布好比合唱团里声部的错位让和声更稳定。直齿行星齿轮纯扭转非线性动力学建模(含3个或4个行星轮)包括动力学方程推导过程考虑了各啮合齿轮副之间的啮合相位可出相图频谱图分岔图庞加莱映射。建立纯扭转模型时每个构件只有绕轴线的转动自由度。定义太阳轮转角θs齿圈θr三个行星轮θp1~θp3。系统动能T的表达式得考虑各部件转动惯量J_s 0.15 # 太阳轮 J_p 0.02 # 单个行星轮 J_r 0.3 # 齿圈 def kinetic_energy(theta_s, theta_p, theta_r): return 0.5*(J_s*theta_s**2 J_r*theta_r**2 sum([J_p*p**2 for p in theta_p]))势能计算才是重头戏得处理各齿轮副的弹性变形。以太阳轮-行星轮1的啮合为例相对位移δsp1 (θs - θp1)rs - (θp1 - θc)rp其中r是基圆半径θc是行星架转角这里假设行星架固定。这式子里的非线性项可能来自齿侧间隙咱们暂且用三次多项式近似def mesh_force(delta, k1e8, c5e3, gap1e-4): if abs(delta) gap: return 0 else: return k*(delta - gap*np.sign(delta)) c*delta_dot # 含迟滞阻尼组装完所有齿轮副的作用力后用拉格朗日方程导出微分方程组时会蹦出耦合项。特别是行星轮间的相位差会让方程出现时变系数处理这类项建议用状态空间法。最终得到的方程大概长这样J_s*θ_s sum(k_spi*(θ_s - θ_pi)) T_input J_p*θ_pi k_spi*(θ_pi - θ_s) k_rpi*(θ_pi - θ_r) 0 J_r*θ_r sum(k_rpi*(θ_r - θ_pi)) -T_load数值求解时推荐用Runge-Kutta法Python里可以用scipy.integrate.odeint。注意步长要小于最快变化的刚度周期比如当转速3000rpm时啮合频率约500Hz步长至少得小于1/(10*500)0.0002秒。跑完仿真后该出图了。相图绘制只需简单plot位移和速度plt.plot(theta_s_history, theta_s_velocity_history, ,b, alpha0.3) plt.xlabel(Angular displacement [rad]) plt.ylabel(Angular velocity [rad/s])频谱分析用FFT要注意加窗不然会出现频谱泄露。试试Blackman-Harris窗from scipy.signal import get_window window get_window(blackmanharris, N) fft_result np.fft.rfft(theta_s_history * window)分岔图需要参数扫描。比如固定输入扭矩逐渐改变转速并记录每个转速下的位移极值。用numpy的广播技巧能加速计算omega_range np.linspace(200, 800, 500) peaks [] for w in omega_range: # 运行仿真并提取极值 peaks.append(find_local_maxima(theta_s_history[-1000:]))庞加莱截面取的是系统状态在特定相位点的采样相当于动力学的快照。取行星轮1啮合相位为0的时刻poincare_samples [] for i in range(len(time)): if abs((omega_mesh*time[i]) % (2*np.pi)) 1e-3: # 相位判断 poincare_samples.append([theta_s_history[i], theta_s_velocity_history[i]])调试这种模型常会遇到数值发散这时候回头检查阻尼系数是否合理。经验之谈阻尼比一般取0.05-0.15对应到Rayleigh阻尼中的β系数得根据系统固有频率调整。另外初始条件别设成全零给个微小扰动才能激发非线性特性。最终得到的各种图中分岔图上的倍周期现象会提示系统进入混沌状态庞加莱点集从有限个点变成连续分布时就得警惕了。这些非线性特征对行星齿轮的故障预警可是重要指标——动力学模型的价值就在于此。