clc;clear all;
V=48;%车速
Rw=0.38;%车轮半径
Mw=300;%车轮质量和1/3轴重
Tw=0.025;%车轮辐板厚度
m=63.4;%单位长度踏面质量
M1=2*pi*Rw*m;%轮重
M2=106;%轴套质量加1/3轴重
Iwr=3.26e-6;%车轮踏面径向惯性矩
E=2.06e11;%车轮弹性模量
K=pi*E*Tw;%车轮刚度
wl=sqrt(K/M2);%低阶固有频率
wh=sqrt(K*Mw/M1*M2);%高阶固有频率
n=1;
% K1=K2=E*Tw.*n/(2*Rw)
K1=E*Tw.*n./(2*Rw);
K2=E*Tw.*n./(2*Rw);
% w0=[125,250,500,1000,2000,4000,8000];
% w0=[200,250,315,400,500,600,800,1000,1250,1600,2000,2500,3150,4000,5000];
w0=[200,250,315,350,400,450,500,550,600,650,700,750,800,850,900,950,1000,1100,1200,1250,1600,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3150,3200,3300,3400,3500,3600,3700,3800,3900,4000,4200,4400,4600,4800,5000];
% w0=100:1:5000;
w=w0.*pi.*2;%定义频率
ee=(1-w.^2./wl^2)./(-j.*w.*Mw.*(1-w.^2./wh.^2));%一阶阻抗
y1(1,:)=(-j.*w.*Rw*n.^2/pi)./((E*Iwr/Rw^2).*n^2.*(n^2-1).^2+Rw^2*(K1.*n^2+K2)-m*Rw^2.*w.^2*(n.^2+1));
for n=2:100
y1(n,:)=(-j.*w.*Rw*(n+1).^2/pi)./((E*Iwr/Rw^2).*(n+1)^2.*((n+1)^2-1).^2+Rw^2*(K1.*(n+1)^2+K2)-m*Rw^2.*w.^2*((n+1).^2+1));
y1(n,:)=y1(n-1,:)+y1(n,:);
end
y2(1,:)=y1(100,:);
y3=ee+y2;
y4=1./y3;
y5=log(y4);
z=abs(y4);
%绘制车轮阻抗图
loglog(w0,z,'linewidth',2);
title('车轮阻抗');
xlabel('频率/ HZ');
ylabel('车轮阻抗 Ns/m')
% grid on
Zw=z;
%计算车轮沿环向均速度
p=66000;%轮静载荷,沿车轮径向分布力
n=1;
An(1,:)=(-n.^2*Rw*p)./(m.*Rw.*w.^2.*(n.^2+1)-E.*Iwr./Rw^3.*n.^2.*(1-n.^2)*(n.^2-1)-Rw.*(K1.*n.^2+K2));
dd(1,:)=norm(An.^2);
% dd=norm(An.^2);
dd_1(1,:)=dd;
ff=An;
ff(1,:)=dd(1,:);
for n=2:100;
An(n,:)=(-n.^2*Rw*p)./(m.*Rw.*w.^2.*(n.^2+1)-E.*Iwr./Rw^3.*n.^2.*(1-n.^2)*(n.^2-1)-Rw.*(K1.*n.^2+K2));
dd(n,:)=norm(An(n,:).^2);
dd_1(n,:)=dd(n,:)+dd_1(n-1,:);
% dd_1=dd(n,:)+dd_1(n-1,:);
ff(n,:)=dd(n,:)+ff(n-1,:);
end
dd_2(1,:)=dd_1(end,:);
% dd_2=dd_1(end,:);
ff_1=ff(end,:);
ff_2=norm(ff_1).^2;
gg=dd_2./(2*ff_2);
save gg;
% save Zw;
