汽车理论第六章 修改
- 格式:ppt
- 大小:4.46 MB
- 文档页数:10
第六章作业:clear;clc;%初始化题干参数如下:f0=1.5;s=0.25;%阻尼比δgama=9;%刚度比γmu=10;%质量比μfs=3;ypss=0.25;g=9.8;ua=20;Gqn0=2.56*10^(-8);%路面不平度系数a0=10^(-6);n0=0.1;detaf=0.2;N=180;f=detaf*(0:N);lambda=f/f0;%车身频率比λlambdas=f/fs;%车轮频率比λW_f=0*f;deta=((1-lambda.^2).*(1+gama-1/mu*lambda.^2)-1).^2+4*s^2*lambda.^2.*(gama-(1/ mu+1)*lambda.^2).^2;z1_q=gama*sqrt(((1-lambda.^2).^2+4*s^2*lambda.^2)./deta); %幅频响应|z1/q|z2_z1=sqrt((1+4*s^2*lambda.^2)./((1-lambda.^2).^2+4*s^2*lambda.^2));%幅频响应|z2/z1|p_z2=sqrt((1+(2*ypss*lambdas).^2)./((1-lambdas.^2).^2+(2*ypss*lambdas).^2));%幅频响应|p/z2|z2_q=gama*sqrt((1+4*s^2*lambda.^2)./deta);p_q=p_z2.*z2_q; %幅频响应|z2/q|sqrt_Gq_f=4*pi^2*f*sqrt(Gqn0*n0^2*ua);sqrt_Gz1_f=z1_q.*sqrt_Gq_f;sqrt_Gz2_f=z2_q.*sqrt_Gq_f;sqrt_Ga_f=p_q.*sqrt_Gq_f;sigmaq=sqrt(trapz(f,sqrt_Gq_f.^2));%路面不平度加速度均方根值sigmaz1=sqrt(trapz(f,sqrt_Gz1_f.^2));%车轮加速度均方根值sigmaz2=sqrt(trapz(f,sqrt_Gz2_f.^2));%车身加速度均方根值sigma_a=sqrt(trapz(f,sqrt_Ga_f.^2));%人体加速度均方根值for i=1:(N+1)if f(i)<=2W_f(i)=0.5;elseif f(i)<=4W_f(i)=f(i)/4;elseif f(i)<=12.5W_f(i)=1;elseW_f(i)=12.5/f(i);endendaww=W_f.^2.*sqrt_Ga_f.^2;aww=sqrt(trapz(f,aww));%加权加速度均方根值Law1=20*log10(aww/a0);%加权振级disp('路面不平度加速度均方根值为');disp(sigmaq);disp('车轮加速度均方根值为');disp(sigmaz1);disp('车身加速度均方根值为');disp(sigmaz2);disp('人体加速度均方根值为');disp(sigma_a);disp('加权加速度均方根值为');disp(aww);disp('加权振级');disp(Law1);figure(1)plot(f,z1_q);title('幅频特性|z1/q|, (f=1.5Hz, δ=0.25,γ=9,μ=10)'),xlabel('激振频率f/Hz'),ylabel('|z1/q|');figure(2)plot(f,z2_z1);title('幅频特性|z2/z1|,(f=1.5Hz, δ=0.25,γ=9,μ=10)'),xlabel('激振频率f/Hz'),ylabel('|z2/z1|');figure(3)plot(f,p_z2);title('幅频特性|p/z2|,(fs=1.5Hz, δs=0.25)'),xlabel('激振频率f/Hz'),ylabel('|p/z2|'); figure(4)plot(f,sqrt_Gz1_f);title('车轮加速度均方根值√Gz1(f)谱图'),xlabel('激振频率f/Hz'),ylabel('√Gz1(f)'); figure(5)plot(f,sqrt_Gz2_f);title('车身加速度均方根值√Gz2(f)谱图'),xlabel('激振频率f/Hz'),ylabel('√Gz2(f)'); figure(6)plot(f,sqrt_Ga_f);title('人体加速度均方根值√Ga(f)谱图'),xlabel('激振频率f/Hz'),ylabel('√Ga(f)'); %第二问:ypss0=(0.125:0.005:0.5);a=0*ypss0;La=0*ypss0;M=length(ypss0);for i=1:Myps=ypss0(i);lambdas=f/fs;deta=((1-lambda.^2).*(1+gama-1/mu*lambda.^2)-1).^2+4*yps^2*lambda.^2.*(gama-(1/mu+1)*lambda.^2).^2;p_z2=sqrt((1+(2*ypss*lambdas).^2)./((1-lambdas.^2).^2+(2*ypss*lambdas).^2));z2_q=gama*sqrt((1+4*yps^2*lambda.^2)./deta);p_q=p_z2.*z2_q;jfg_Gqddf=4*pi^2*sqrt(Gqn0*n0^2*ua)*f;jfg_Gaf=p_q.*jfg_Gqddf;kk=W_f.^2.*jfg_Gaf.^2;aww(i)=sqrt(trapz(f,kk));endLaw1=20*log10(aww/a0);figure(7)plot(ypss0,aww);title('aw随δs的变化'),xlabel('“人体—座椅”系统的阻尼比δs'),ylabel('aw/m*s^-2');figure(8)plot(ypss0,Law1);title('Law随δs的变化'),xlabel('“人体—座椅”系统的阻尼比δs'),ylabel('Law/dB'); fs=1.5:0.025:6;k=length(fs);for i=1:kfs0=fs(i);lambdas1=f/fs0;deta=((1-lambda.^2).*(1+gama-1/mu*lambda.^2)-1).^2+4*yps^2*lambda.^2.*(gama-(1/mu+1)*lambda.^2).^2;p_z2=sqrt((1+(2*ypss*lambdas1).^2)./((1-lambdas1.^2).^2+(2*ypss*lambdas1).^2));z2_q=gama*sqrt((1+4*yps^2*lambda.^2)./deta);p_q=p_z2.*z2_q;jfg_Gqddf1=4*pi^2*sqrt(Gqn0*n0^2*ua)*f;jfg_Ga1f=p_q.*jfg_Gqddf1;kk1=W_f.^2.*jfg_Ga1f.^2;aww(i)=sqrt(trapz(f,kk1));endLaw2=20*log10(aww/a0);figure(9)plot(fs,aww);title('aw随fs的变化'),xlabel('“人体—座椅”系统的固有频率fs'),ylabel('aw/m*s^-2');figure(10)plot(fs,Law2);title('Law随fs的变化'),xlabel('“人体—座椅”系统的固有频率fs'),ylabel('Law/dB'); %第三问:%三个响应量均方根值随f0变化的曲线ff0=0.25:0.05:3;sigmaz2=0*ff0;sigmafd=0*ff0;sigmaFd_G=0*ff0;M=length(ff0);for i=1:Mf0=ff0(i);f=detaf*(0:N);lamta=f/f0;lamtas=f/fs;deta=((1-lamta.^2).*(1+gama-1/mu*lamta.^2)-1).^2+4*yps^2*lamta.^2.*(gama-(1/m u+1)*lamta.^2).^2;z2_qdot=2*pi*f*gama.*sqrt((1+4*yps^2*lamta.^2)./deta);fd_qdot=gama*lamta.^2./(2*pi*f+eps)./sqrt(deta);Fd_Gqdot=2*pi*f*gama/g.*sqrt(((lamta.^2/(mu+1)-1).^2+4*yps^2*lamta.^2)./deta);Gq_dotf=4*pi^2*Gqn0*n0^2*ua;Gz2f=(z2_qdot).^2*Gq_dotf;Gfd_qf=(fd_qdot).^2*Gq_dotf;GFd_Gf=(Fd_Gqdot).^2*Gq_dotf;sigmaz2(i)=sqrt(trapz(f,Gz2f));sigmafd(i)=sqrt(trapz(f,Gfd_qf));sigmaFd_G(i)=sqrt(trapz(f,GFd_Gf));if f0==1.5sgmz2=sigmaz2(i);sgmfd=sigmafd(i);sgmFd_G=sigmaFd_G(i);endendsz2=20*log10(sigmaz2/sgmz2);sfd=20*log10(sigmafd/sgmfd);sFd_G=20*log10(sigmaFd_G/sgmFd_G);plot(ff0,sz2,'r-',ff0,sfd,'b-.',ff0,sFd_G,'k--');l egend('σz2','σfd','σFd/G');gtext('σz2');gtext('σfd');gtext('σFd/G');axis([0.25 3 -25 15]);title('三个响应量均方根值随f0变化的曲线'),xlabel('车身部分固有频率f0/Hz'),ylabel('σz2/dB,σfd/dB,σFd/G/dB');%三个响应量均方根值随δ变化的曲线c=(0.5-0.125)/180;yps0=0.125:c:0.5;sigmaz2=0*yps0;sigmafd=0*yps0;sigmaFd_G=0*yps0;M=length(yps0);for i=1:Myps=yps0(i);f=detaf*(0:N);lamta=f/f0;lamtas=f/fs;deta=((1-lamta.^2).*(1+gama-1/mu*lamta.^2)-1).^2+4*yps^2*lamta.^2.*(gama-(1/m u+1)*lamta.^2).^2;z2_qdot=2*pi*f*gama.*sqrt((1+4*yps^2*lamta.^2)./deta);fd_qdot=gama*lamta.^2./(2*pi*f+eps)./sqrt(deta);Fd_Gqdot=2*pi*f*gama/g.*sqrt(((lamta.^2/(mu+1)-1).^2+4*yps^2*lamta.^2)./deta);Gq_dotf=4*pi^2*Gqn0*n0^2*ua;Gz2f=(z2_qdot).^2*Gq_dotf;Gfd_qf=(fd_qdot).^2*Gq_dotf;GFd_Gf=(Fd_Gqdot).^2*Gq_dotf;sigmaz2(i)=sqrt(trapz(f,Gz2f));sigmafd(i)=sqrt(trapz(f,Gfd_qf));sigmaFd_G(i)=sqrt(trapz(f,GFd_Gf));if yps==0.25sgmz2=sigmaz2(i);sgmfd=sigmafd(i);sgmFd_G=sigmaFd_G(i);endendsz2=20*log10(sigmaz2/sgmz2);sfd=20*log10(sigmafd/sgmfd);sFd_G=20*log10(sigmaFd_G/sgmFd_G);plot(yps0,sz2,'r-',yps0,sfd,'b-.',yps0,sFd_G,'k--');axis([0.125 0.5 -4 4]);title('三个响应量均方根值随δ变化的曲线'),xlabel('车身部分阻尼比δ'),ylabel('σz2/dB,σfd/dB,σFd/G/dB');legend('σz2','σfd','σFd/G');gtext('σz2');gtext('σfd');gtext('σFd/G');%三个响应量均方根值随γ变化的曲线gama0=4:0.1:19;sigmaz2=0*gama0;sigmafd=0*gama0;sigmaFd_G=0*gama0;M=length(gama0);for i=1:Mgama=gama0(i);f=detaf*(0:N);lamta=f/f0;lamtas=f/fs;deta=((1-lamta.^2).*(1+gama-1/mu*lamta.^2)-1).^2+4*yps^2*lamta.^2.*(gama-(1/m u+1)*lamta.^2).^2;z2_qdot=2*pi*f*gama.*sqrt((1+4*yps^2*lamta.^2)./deta);fd_qdot=gama*lamta.^2./(2*pi*f+eps)./sqrt(deta);Fd_Gqdot=2*pi*f*gama/g.*sqrt(((lamta.^2/(mu+1)-1).^2+4*yps^2*lamta.^2)./deta);Gq_dotf=4*pi^2*Gqn0*n0^2*ua;Gz2f=(z2_qdot).^2*Gq_dotf;Gfd_qf=(fd_qdot).^2*Gq_dotf;GFd_Gf=(Fd_Gqdot).^2*Gq_dotf;sigmaz2(i)=sqrt(trapz(f,Gz2f));sigmafd(i)=sqrt(trapz(f,Gfd_qf));sigmaFd_G(i)=sqrt(trapz(f,GFd_Gf));if gama==9sgmz2=sigmaz2(i);sgmfd=sigmafd(i);sgmFd_G=sigmaFd_G(i);endendsz2=20*log10(sigmaz2/sgmz2);sfd=20*log10(sigmafd/sgmfd);sFd_G=20*log10(sigmaFd_G/sgmFd_G);figure(13)plot(gama0,sz2,'r-',gama0,sfd,'b-.',gama0,sFd_G,'k--');axis([4 18 -5 6]);title('三个响应量均方根值随γ变化的曲线'),xlabel('悬架与轮胎的刚度比γ'),ylabel('σz2/dB,σfd/dB,σFd/G/dB');legend('σz2','σfd','σFd/G');gtext('σz2');gtext('σfd');gtext('σFd/G');mu0=5:0.1:20;sigmaz2=0*mu0;sigmafd=0*mu0;sigmaFd_G=0*mu0;M=length(mu0);for i=1:Mmu=mu0(i);f=detaf*(0:N);lamta=f/f0;lamtas=f/fs;deta=((1-lamta.^2).*(1+gama-1/mu*lamta.^2)-1).^2+4*yps^2*lamta.^2.*(gama-(1/m u+1)*lamta.^2).^2;z2_qdot=2*pi*f*gama.*sqrt((1+4*yps^2*lamta.^2)./deta);fd_qdot=gama*lamta.^2./(2*pi*f+eps)./sqrt(deta);Fd_Gqdot=2*pi*f*gama/g.*sqrt(((lamta.^2/(mu+1)-1).^2+4*yps^2*lamta.^2)./deta);Gq_dotf=4*pi^2*Gqn0*n0^2*ua;Gz2f=(z2_qdot).^2*Gq_dotf;Gfd_qf=(fd_qdot).^2*Gq_dotf;GFd_Gf=(Fd_Gqdot).^2*Gq_dotf;sigmaz2(i)=sqrt(trapz(f,Gz2f));sigmafd(i)=sqrt(trapz(f,Gfd_qf));sigmaFd_G(i)=sqrt(trapz(f,GFd_Gf));if mu==10sgmz2=sigmaz2(i);sgmfd=sigmafd(i);sgmFd_G=sigmaFd_G(i);endendsz2=20*log10(sigmaz2/sgmz2);sfd=20*log10(sigmafd/sgmfd);sFd_G=20*log10(sigmaFd_G/sgmFd_G);figure(14)plot(mu0,sz2,'r-',mu0,sfd,'b-.',mu0,sFd_G,'k--'); axis([5 20 -2 2]);title('三个响应量均方根值随μ变化的曲线'),xlabel('车身与车轮部分质量比μ'),ylabel('σz2/dB,σfd/dB,σFd/G/dB'); legend('σz2','σfd','σFd/G');gtext('σz2');gtext('σfd');gtext('σFd/G');计算结果:第一问结果:路面不平度加速度均方根值为0.3523车轮加速度均方根值为0.2391车身加速度均方根值为0.0168人体加速度均方根值为0.0161加权加速度均方根值为0.0100加权振级80.0291第二问的答案:第三问答案:。
余志生汽车理论第六章课后习题答案6.l 、设通过座椅支承面传至人体垂直加速度的谱密度为一白噪声,Ga ( f )=0.132m -⋅s 。
求在0.5~80H Z 频率范围加权加速度均方根值a w 和加权振级L aw ,并由表6-2查出相应人的主观感觉。
答:21805.02])()([df f G f W a a w ⎰⋅=805.125.1244225.05.121.011.041.0*5.0[dff df df f df ⎰⎰⎰⎰+⋅⋅+⋅⋅+⋅=28.24=⇒)(200a a Lg L waw=70.147)1028.24(206==-Lg查173P 图知:人的主观感觉为极不舒适。
6.2、设车速u =20m /s ,路面不平度系380q 10*56.2)(G m n -=,参考空间频率n o =0.1-1m 。
画出路面垂直位移、速度和加速度)(G q f 、)(G q f 、)(G q f 的谱图。
画图时要求用双对数坐标,选好坐标刻度值,并注明单位。
解:228220q 20*1.0*10*56.2)()(G f f u n n G f q -==29110*12.5f-= 20*1.0*10*56.2*4)(4)(G 282202q -==ππu n n G f q-710*2.02=22842204q *1.0*10*56.2*16)(16)(G f uf n n G f q -==ππ 2-710*99.3f =画出图形为:6.3、设车身-车轮二自由度汽车模型,其车身部分固有频率f o =2Hz 。
它行驶在波长λ=5m 的水泥接缝路上,求引起车身部分共振时的车速u n (km/h)。
该汽车车轮部分的固有频率f t =10Hz ,在砂石路上常用车速为30km/h 。
问由于车轮部分共振时,车轮对路面作用的动载所形成的搓板路的波长λ=?答:①当激振力等于车辆固有频率时,发生共振,所以发生共振时的车速为:2*5u 0a =⋅=f λs m /10=②搓板路的波长 :m 65106.3/30==λ6.4、设车身单质量系统的幅频 |z /q | 用双对数坐标表示时如习题图6所示。