车辆系统动力学仿真大作业(带程序).docx

《车辆系统动力学仿真大作业(带程序).docx》由会员分享，可在线阅读，更多相关《车辆系统动力学仿真大作业(带程序).docx（17页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、车辆系统动力学仿真大作业(带程序)文档视界车辆系统动力学仿真大作业(带程序)车辆系统动力学仿真大作业(带程序)文档视界车辆系统动力学仿真大作业(带程序)车辆系统动力学仿真大作业(带程序)AccordingtotheNewssecondLaw,wecangettheequation:2)()(221211mgzzczzkzm-+-=?221212)()(zkmgzzczzkzmw+-=?0)()()()(222111222111=-+-+-+-+?wwwwzLzkzLzkzLzczLzczm0)()()()(2222111122221111=-+-+-?wwwwzLzLkzLzLkzLzLczL

2、zLcJdwwwwQzLzkzLzczm,111111111)()(-=-?dwwwwQzLzkzLzczm,222222222)()(-=-+-+-?Whenthereisnoexcitationwecangettheequation:2)()(221211mgzzczzkzm-+-=?221212)()(zkmgzzczzkzmw+-=?Thenwesubstitudethedataintotheequation,wewriteaproceduretosimulatethesystem:Date:?=?=?=MN/m0.10km25.1s/mkN0.20MN/m0.1mkg3020kg21

3、00kg3250w2lckImmbywbProgram:1.Forthenoexcitationsituation.Buidafilenamedrigid1.mfunctiondy=rigid1(t,y)dy=zeros(4,1);dy(1)=y(2);dy(3)=y(4);dy(2)=-(615.38*y(1)-615.38*y(3)+12.31*dy(1)-12.31*dy(3)-9.81);dy(4)=-(5238.1*y(3)-476.19*y(1)+9.52*dy(3)-9.52*dy(1)-9.81);Buidanotherfilenamedtest1.mT,Y=ode45(rig

4、id1,01.5,0000);figure(1)plot(T,Y(:,1)gridon,xlabel(time(sec),ylabel(displayment(m),title(displaymentofboge)figure(2)plot(T,Y(:,2)gridon,xlabel(time(sec),ylabel(velocity(m/s),title(velocityofboge)figure(3)plot(T,Y(:,3)gridon,xlabel(time(sec),ylabel(displayment(m),title(displaymentofwheel)figure(4)plo

5、t(T,Y(:,4)gridon,xlabel(time(sec),ylabel(velocity(m/s),title(velocityofwheel)文档视界车辆系统动力学仿真大作业(带程序)车辆系统动力学仿真大作业(带程序)文档视界车辆系统动力学仿真大作业(带程序)车辆系统动力学仿真大作业(带程序)Whenthereareexcitations,accordingtoNewtonssecondLawandmomentumequationwecanget0)()()()(32321=-+-+-+-+?zlzkzlzkzlzczlzczm0)()()()(3232=-+-+-?zlzklzl

6、zklzlzclzlzlcJdQzlzkzlzczm,222)()(22-=-?dQzlzkzlzczm,33333)()(-=-+-+-?)(22,2xzkQwd-=)(33,3xzkQwd-=Therearefourexcitations,eachwheelsetgetoneexcitation.X1=0.005sin(4t)X2=0.005sin(8t)X1andX2actonboge1X3=0.005sin(4t)X4=0.005sin(6t)X1andX2actonboge2Program:Mw=2100;%wheelmassMb=3250;%bogemassMc=0;%carbod

7、ymassKw=1.0e7;%wheelelasticcoeffetionKs=0e6;%secendarysuspensioncoeffetionKp=1.0e6;%primerysuspensioncoeffetionCs=1;%secendarydumpcoeffetionCp=5.0e4;%primerysuspensioncoeffetionlc=0;%carbodylenthlt=1.25;%bogelenthic=1;%carbodyinertiaib=3020;%bogeinertiaA=-(Kp+Cp+Kw+Mw),0,0,0,(Kp+Cp),0,0,-lt*(Kp+Cp),

8、0,0;0,-(Kp+Cp+Kw+Mw),0,0,(Kp+Cp),0,0,-lt*(Kp+Cp),0,0;0,0,-(Kp+Cp+Kw+Mw),0,0,(Kp+Cp),0,0,-lt*(Kp+Cp),0;0,0,0,-(Kp+Cp+Kw+Mw),0,(Kp+Cp),0,0,-lt*(Kp+Cp),0;(Kp+Cp),(Kp+Cp),0,0,-(Ks+Cs+2*Kp+2*Cp+Mb),0,(Ks+Cs),0,0,-lc*(Ks+Cs);0,0,(Kp+Cp),(Kp+Cp),0,-(Ks+Cs+2*Kp+2*Cp+Mb),(Ks+Cs),0,0,-lc*(Ks+Cs);0,0,0,0,(Ks+C

9、s),(Ks+Cs),-(2*Ks+2*Cs+Mc),0,0,0;lt*(Kp+Cp),-lt*(Kp+Cp),0,0,0,0,0,(2*lt2*(Kp+Cp)-ib),0,0;0,0,(lt*Kp+lt*Cp),-(lt*Kp+lt*Cp),0,0,0,0,(2*lt2*(Kp+Cp)-ib),0;0,0,0,0,(lc*Ks+lc*Cs),-(lc*Ks+lc*Cs),0,0,0,(2*lc2*(Ks+Cs)-ic);B=-(2*Mw+Cp),0,0,0,Cp,0,0,-Cp*lt,0,0;0,-(2*Mw+Cp),0,0,Cp,0,0,Cp*lt,0,0;0,0,-(2*Mw+Cp),0

10、,0,Cp,0,0,-Cp*lt,0;0,0,0,-(2*Mw+Cp),0,Cp,0,0,Cp*lt,0;Cp,Cp,0,0,-(Cs+2*Cp+2*Mb),0,Cs,0,0,-Cs*lc;0,0,Cp,Cp,0,-(Cs+2*Cp+2*Mb),Cs,0,0,-Cs*lc;0,0,0,0,Cs,Cs,-(2*Mc+2*Cs),0,0,0;lt*Cp,-lt*Cp,0,0,0,0,0,(2*lt2*Cp-2*ib),0,0;0,0,lt*Cp,-lt*Cp,0,0,0,0,(2*lt2*Cp-2*ib),0;0,0,0,0,lc*Cs,-lc*Cs,0,0,0,(2*lc2*Cs-2*ic);C

11、=diag(Mw,Mw,Mw,Mw,Mb,Mb,Mc,ib,ib,ic);D=9.8*Mw,Mw,Mw,Mw,Mb,Mb,Mc,0,0,0;Ts=0.001;%seperate10sinto10000partst=0:Ts:10;x1=0.05*sin(4*pi*t);%theexcitationofwheel1x2=0.05*sin(8*pi*t);%theexcitationofwheel2x3=0.05*sin(4*pi*t);%theexcitationofwheel3x4=0.05*sin(6*pi*t);%theexcitationofwheel4N=4000;YY=zeros(1

12、0,2);VV=zeros(10,2);AA=zeros(10,2);forj=3:NXn=Kw*x1(j-2);x2(j-2);x3(j-2);x4(j-2);0;0;0;0;0;0;YY(:,j)=AB*YY(:,j-1)+AC*YY(:,j-2)-AXn-AD;%displaymentmatrixVV(:,j)=YY(:,j)-YY(:,j-1);%speedmatrixAA(:,j)=VV(:,j)-VV(:,j-1);%accelerationmatrixendi=3:N;figure(1)plot(i,YY(1,i)gridon,xlabel(time(Ts),ylabel(dis

13、playment(m),title(displaymentofwheel1);figure(2)plot(i,YY(2,i)gridon,xlabel(time(Ts),ylabel(displayment(m),title(displaymentofwheel2);figure(3)plot(i,YY(3,i)gridon,xlabel(time(Ts),ylabel(displayment(m),title(displaymentofwheel3);figure(4)plot(i,YY(4,i)gridon,xlabel(time(Ts),ylabel(displayment(m),tit

14、le(displaymentofwheel4);figure(5)plot(i,YY(5,i)gridon,xlabel(time(Ts),ylabel(displayment(m),title(displaymentofboge1);figure(6)plot(i,YY(6,i)gridon,xlabel(time(Ts),ylabel(displayment(m),title(displaymentofboge2);figure(7)plot(i,YY(8,i)gridon,xlabel(time(Ts),ylabel(angulardisplayment(m),title(angular

15、displaymentofboge1);figure(8)plot(i,YY(9,i)gridon,xlabel(time(Ts),ylabel(angulardisplayment(m),title(angulardisplaymentofboge2);figure(9)plot(i,VV(1,i)gridon,xlabel(time(Ts),ylabel(velocity(m/s),title(velocityofwheel1);figure(10)plot(i,VV(2,i)gridon,xlabel(time(Ts),ylabel(velocity(m/s),title(velocit

16、yofwheel2);figure(11)plot(i,VV(3,i)gridon,xlabel(time(Ts),ylabel(velocity(m/s),title(velocityofwheel3);figure(12)plot(i,VV(4,i)gridon,xlabel(time(Ts),ylabel(velocity(m/s),title(velocityofwheel4);figure(13)plot(i,VV(5,i)gridon,xlabel(time(Ts),ylabel(velocity(m/s),title(velocityofboge1);figure(14)plot

17、(i,VV(6,i)gridon,xlabel(time(Ts),ylabel(velocity(m/s),title(velocityofboge2);figure(15)plot(i,VV(8,i)gridon,xlabel(time(Ts),ylabel(angulervelocity(rad/s),title(angularvelocityofboge1);figure(16)plot(i,VV(9,i)gridon,xlabel(time(Ts),ylabel(angulervelocity(rad/s),title(angularvelocityofboge2);figure(17

18、)plot(i,AA(1,i)gridon,xlabel(time(Ts),ylabel(acceleration(m/s2),title(accelerationofwheel1);figure(18)plot(i,AA(2,i)gridon,xlabel(time(Ts),ylabel(acceleration(m/s2),title(accelerationofwheel2);figure(19)plot(i,AA(3,i)gridon,xlabel(time(Ts),ylabel(acceleration(m/s2),title(accelerationofwheel3);figure(20)plot(i,AA(4,i)gridon,xlabel(time(Ts),ylabel(acceleration(m/s2),title(accelerationofwheel4);figure(21)plot(i,AA(5,i)gridon,xlabel(time(Ts),ylabel(acceleration(m/s2),title(accelerationofboge1);figure(22)plot(i,AA(6,i)gridon,xlabel(time(Ts),ylabel(acceleration(m/s2),title(accelerationofboge2);