车辆和转向架蛇行运动模拟实验——实验报告——林硕、陈俊龙、陈虹、谢天、徐育颋 参数选择如下Rr1(m) Rr2(m) Rw1(m) Rw2(m) 200*10^-360*10^-3 100*10^-3 inf 200*9,8 Rr2 的选择依据为:根据 60 轨的尺寸,其接触点的曲率半径为 300mm,因为模型为 1:5, 所以为60mm。 认为模型轮轨的弹性模量和泊松比为一样的,E=200GPa,μ=0.3 计算程序如下:%单位制可以为m-kg-N-Pa clear; clc; 初始参数输入E1=200*10^9;%Pa E2=200*10^9;%Pa u1=0.3; u2=0.3; Rr1=200*10^-3;%m Rr2=300*10^-3/5;%m 1:5 的模型 Rw1=100*10^-3;%m Rw2=inf;%infinte P=200*9.8;%N 开始计算AplusB=0.5*(1/Rr1+1/Rr2+1/Rw1+1/Rw2); BminusA=0.5*(1/Rr1-1/Rr2+1/Rw1-1/Rw2); thelt=acos(abs(BminusA/AplusB)); thelt=thelt/pi*180; theltTemp=[30,40,50,60,70,80,90]; mTemp=[2.73,2.14,1.75,1.49,1.28,1.13,1.00]; nTemp=[0.49,0.57,0.64,0.72,0.80,0.89,1.00]; m=interp1(theltTemp,mTemp,thelt);%进行线性差值 n=interp1(theltTemp,nTemp,thelt); k1=(1-u1^2)/(pi*E1); k2=(1-u2^2)/(pi*E2); a=m*(3*pi*P*(k1+k2)/(4*AplusB))^(1/3); b=n*(3*pi*P*(k1+k2)/(4*AplusB))^(1/3); BminusA/AplusB