垂向动力学作业3和4

垂向动力学作业3和4垂向动力学作业3和4垂向动力学作业3和4第二次作业课程作业三解:G∵上式及下文用到的简化表达式：A∴解：G∵∴解：G∵=∴解：G∵=ω=ω=ω∴解：G∵=∴课程作业四解：车身质心垂向加速度增益：Z∴Z∵ZZH车身俯仰角增益:ϕ=ϕ解:车身质心:G前轮：G∵后轮：G∵解：悬架动行程（前轮）:G∵同理得后轮悬架动行程：G轮胎动载荷（前轮)：G∵∴同理得后轮轮胎：G源程序:课程作业三clearm_A=256；m_B=31；d_A=1150;c_A=20.2e3;c_R=128e3;％车身—车轮双质量系统参数f0=1。5;zeta=d_A。/（2*sqrt（c_A。＊m_A））;%ζgamma=c_R。/c_A;％γmu=m_A./m_B；％μw0=sqrt（c_A./m_A）;％ω0％路面和车速u=20;Gqn0=16e—6;n0=.1;％定义频率域f=linspace（0，16，1000)；w=2＊pi*f;%ωlam=w/w0;%λdelta=((1-(w/w0)。^2）。*(1+gamma-(w/w0）。^2/mu）—1)。^2+4*zeta^2＊(w/w0）。^2.＊(gamma-(1/mu+1）*（w/w0).^2）。^2；%%A1=1i*w＊d_A+c_A;A2=-w。^2*m_A+1i*w＊d_A+c_A；A3=—w。^2＊m_B+1i*w*d_A+c_A+c_R；N=A3.＊A2-A1.^2；%|zA/zh|zA_zh=gamma＊（（1+4*zeta.^2。＊lam。^2)./delta).^.5;%%＃1G1=4*pi^2＊Gqn0＊n0^2*u。＊w。^2.*zA_zh.^2;figure（1);plot（f，G1，’linewidth’，2）；title('车身加速度Z^’'^’'_A'，'fontsize'，12）xlabel('/Hz’，'fontsize',12)ylabel（'/(m^2/s^3)’，'fontsize',12）％%#2Hw_2=abs((A1.＊c_R-N）。/N）；G2=4*pi^2*Gqn0*n0^2＊u。＊w.^-2.＊(Hw_2)。^2；figure(2);plot（f,G2，'linewidth',2);title('车身离地面距离Z_A—Z_h'，'fontsize',12）xlabel（'/Hz’,’fontsize',12）ylabel(’/m^2\timess'，’fontsize’，12）%%#3Hw_3=abs(c_A＊c_R*(A1—A2）./N);G3=4*pi^2*Gqn0＊n0^2＊u。*w.^—2.*(Hw_3).^2；figure(3)；plot（f,G3,’linewidth’，2);title（’弹簧力动载荷C_A（Z_A—Z_R）'，'fontsize'，12)xlabel（'/Hz'，'fontsize’，12)ylabel（’/N^2\timess’，'fontsize’，12)%%＃4Hw_4=abs(d_A＊c_R*(A1-A2)。/N)；G4=4*pi^2＊Gqn0*n0^2*u*(Hw_4).^2;figure(4)；plot(f，G4,’linewidth'，2)；title（’减震器力动载荷d_A（Z^'’_A-Z^'’_R)’,'fontsize',12）xlabel（'/Hz’,'fontsize’,12）ylabel(’/N^2/s','fontsize’,12)％%#5Hw_5=abs(c_R＊(c_R*A2./N-1）);G5=4＊pi^2*Gqn0＊n0^2*u。*w.^-2。*(Hw_5)。^2;figure（5）;plot(f,G5,’linewidth’,2);title('车轮动载荷C_R(Z_R—Z_h）’,'fontsize’，12)xlabel(’/Hz'，’fontsize’，10）ylabel(’/N^2\timess’，’fontsize’，10）课程作业四clearclcm1f=40。5;m1r=45.4;m2f=377。5;m2r=312.5；a=1。25;b=1.51;L=2。76；k1f=192e3；k1r=192e3;k2f=17e3;k2r=22e3；c2f=1500;c2r=1500;u=20；Gqn0=16e—6;n0=.1；f=linspace(0,16，1000）；w=2*pi＊f;%ω%%frontA1_f=1i＊w*c2f+k2f；A2_f=—w。^2*m2f+1i*w*c2f+k2f；A3_f=-w。^2*m1f+1i＊w*c2f+k2f+k1f;N_f=A3_f。＊A2_f—A1_f.^2;%%rearA1_r=1i*w＊c2r+k2r；A2_r=-w.^2＊m2r+1i＊w*c2r+k2r;A3_r=-w。^2*m1r+1i*w*c2r+k2r+k1r;N_r=A3_r.*A2_r-A1_r。^2;%%质心垂向位移Hw_z2=w。^2。＊abs（b/L*（A1_f＊k1f./N_f）+a/L*(A1_r＊k1r。/N_r)。*exp(-1i＊L*w/u));G_z2=4＊pi^2＊Gqn0*n0^2*u.＊w。^—2.*Hw_z2。^2；％%俯仰角Hw_phi_y=abs（1/L*((A1_r*k1r。/N_r）。*exp(—1i＊L*w/u）—（A1_f＊k1f。/N_f))）;％G_phi_y=Hw_phi_y.*4＊pi^2*Gqn0*n0^2＊u。*w。^-2;%％（1)图figure（1）plot(f，Hw_z2，’linewidth'，2);title('车身质心处垂向加速度增益’，’fontsize',12)xlabel（'/Hz’，'fontsize’，12)ylabel(’/s^-^2’，’fontsize’，12)figure(2)plot(f，Hw_phi_y,’linewidth’，2);title（'车身质心处俯仰角增益'，’fontsize',12）xlabel(’/Hz’，'fontsize’，12）ylabel('/m^-^1’，'fontsize’,12)%％前轮Hw_z1f=w。^2.*abs(A2_f*k1f。/N_f);G_z1f=4＊pi^2*Gqn0*n0^2*u。＊w.^—2.*Hw_z1f。^2;%％后轮Hw_z1r=w.^2。＊abs(A2_r＊k1r。/N_r)；G_z1r=4*pi^2*Gqn0*n0^2＊u.*w。^-2。*Hw_z1r。^2;figure（3)plot(f，G_z2,’linewidth’，2）title（'车身质心垂向加速度的功率谱密度’，'fontsize’,12）xlabel（’/Hz’，'fontsize’，12)ylabel(’/（m^2\timess^3)'，'fontsize'，12）figure（4)plot(f，G_z1f，f，G_z1r,'——'，'linewidth',2）title(’前轮和后轮垂向加速度的功率谱密度’，'fontsize'，12)xlabel（’/Hz'，'fontsize’，12）ylabel（'/（m^2\timess^3)’，’fontsize'，12)legend（’前轮’，'后轮’）％％（3）悬架动行程％悬架动行程Hw_fd_f=abs（k1f./N_f。*(A1_f-A2_f））；%前悬G_fd_f=4＊pi^2*Gqn0＊n0^2*u。*w。^—2。＊Hw_fd_f。^2；Hw_fd_r=abs（k1r。/N_r。＊(A1_r—A2_r）)；%后悬G_fd_r=4＊pi^2＊Gqn0*n0^2＊u.＊w。^—2。＊Hw_fd_r.^2;figure（5)plot(f，G_fd_f,f，G_fd_r,'——’,’linewidth',2)title（'悬架动行程功率谱密度'，’fontsize’，12）xlabel(’/Hz’，’fontsize’,12）ylabel（’/m^2\timess’，'fontsize’，12)legend（'前轮’,'后轮')%轮胎动载荷Hw_Fd_f=abs(k1f*（A2_f*k1f。/N_f-1)）;%前悬G_Fd_f=4＊pi^2＊Gqn0＊n0^2*u.＊w。^-2.*Hw_Fd_f.^2;Hw_Fd_r=abs(k1r*(A2_r＊k1r。/N_r-1));%后悬G_Fd_r=