车辆系统动力学计算机仿真.ppt

第四章车辆系统动力学数值分析方法,车辆系统动力学计算机仿真,第四章车辆系统动力学数值分析方法,一、模型,第四章车辆系统动力学数值分析方法,二、参数,第四章车辆系统动力学数值分析方法,三、激扰,第四章车辆系统动力学数值分析方法,implicitdoubleprecision(a-h,o-z)parameter(nn=40000,nd=2,gg=9.81,pi=3.14159)parameter(th=1000.)COMMON/EXPL/step,F1KC,FKC,FI1,FIdimensionzb1(nd),vb1(nd),ab1(nd)dimensiondzb1(nd),dvb1(nd),dab1(nd),d1ab1(nd)open(11,file=rusta.dat,status=unknown)open(12,file=rustd.dat,status=unknown),第四章车辆系统动力学数值分析方法,cmk=1.0e+6;dmp=2.0e+4cbk=1.0e+7;dbp=2.0e+4bm1=3250.0bm2=2100.0bim=3020.0bl=2.5,第四章车辆系统动力学数值分析方法,F1KC=1.0FKC=0.5fi1=1.5fi=0.5step=1.0e-4,第四章车辆系统动力学数值分析方法,doi=1,nn,1tt=dble(i)*stepzt2=0.05*sin(4.*pi*tt)zt3=0.05*sin(4.*pi*tt)CALLEXP1(dzb1,dvb1,dab1,d1ab1,zb1,vb1,nd)!积分CALLEXP2(dzb2,dvb2,dab2,d1ab2,zb2,vb2)!积分CALLEXP2(dzb3,dvb3,dab3,d1ab3,zb3,vb3)!积分,第四章车辆系统动力学数值分析方法,!求弹簧悬挂点质量1的位移与速度dm2=zb1(1)+0.5*bl*zb1(2)vm2=vb1(1)+0.5*bl*vb1(2)dm3=zb1(1)-0.5*bl*zb1(2)vm3=vb1(1)-0.5*bl*vb1(2),第四章车辆系统动力学数值分析方法,fm2=cmk*(dm2-zb2)+dmp*(vm2-vb2)fm3=cmk*(dm3-zb3)+dmp*(vm3-vb3)fb2=cbk*(zb2-zt2)+dbp*vb2fb3=cbk*(zb3-zt3)+dbp*vb3,第四章车辆系统动力学数值分析方法,ab1(1)=(-fm2-fm3)/bm1-ggab1(2)=0.5*bl*(fm3-fm2)/bimab2=(-fb2+fm2)/bm2-ggab3=(-fb3+fm3)/bm2-gg,第四章车辆系统动力学数值分析方法,if(mod(nn,800).eq.0)thenwrite(11,100)tt,ab1(1),ab2,ab3,zt2write(12,100)tt,zb1(1)*th,zb2*th,zb3*th,zt3endifCALLEXCH1(zb1,vb1,ab1,dzb1,dvb1,dab1,d1ab1,ND)CALLEXCH2(zb2,vb2,ab2,dzb2,dvb2,dab2,d1ab2)CALLEXCH2(zb3,vb3,ab3,dzb3,dvb3,dab3,d1ab3)enddo100format(1x,f12.5,2x,f12.5,2x,f12.5,2x,f12.5,2x,f12.5)stopend,第四章车辆系统动力学数值分析方法,SUBROUTINEEXP1(DY,DV,DA,D1A,CY,CV,ND)IMPLICITDOUBLEPRECISION(A-H,O-Z)COMMON/EXPL/VST,F1KC,FKC,FI1,FIDIMENSIONDY(ND),DV(ND),DA(ND),D1A(ND),CY(ND),CV(ND)VST2=VST*VSTDO700I=1,ND,1CY(I)=DY(I)+DV(I)*VST+(F1KC*DA(I)-FKC*D1A(I)*VST2CV(I)=DV(I)+(FI1*DA(I)-FI*D1A(I)*VST700CONTINUERETURNEND,第四章车辆系统动力学数值分析方法,SUBROUTINEEXP2(DY,DV,DA,D1A,CY,CV)IMPLICITDOUBLEPRECISION(A-H,O-Z)COMMON/EXPL/VST,F1KC,FKC,FI1,FIVST2=VST*VSTCY=DY+DV*VST+(F1KC*DA-FKC*D1A)*VST2CV=DV+(FI1*DA-FI*D1A)*VSTRETURNEND,第四章车辆系统动力学数值分析方法,SUBROUTINEEXCH1(ZZC,VZC,AZC,DZZC,DVZC,DAZC,D1AZC,ND)IMPLICITDOUBLEPRECISION(A-H,O-Z)DIMENSIONZZC(ND),VZC(ND),AZC(ND)DIMENSIONDZZC(ND),DVZC(ND),DAZC(ND),D1AZC(ND)DO5200I=1,ND,1DZZC(I)=ZZC(I)DVZC(I)=VZC(I)D1AZC(I)=DAZC(I)DAZC(I)=AZC(I)ZZC(I)=0.VZC(I)=0.AZC(I)=0.5200CONTINUERETURNEND,第四章车辆系统动力学数值分析方法,SUBROUTINEEXCH2(ZZC,VZC,AZC,DZZC,DVZC,DAZC,D1AZC)IMPLICITDOUBLEPRECISION(A-H,O-Z)DZZC=ZZCDVZC=VZCD1AZC=DAZCDAZC=AZCZZC=0.VZC=0.AZC=0.RETURNEND,第四章车辆系统动力学数值分析方法,一、无外加激扰（自由衰减振动）,第四章车辆系统动力学数值分析方法,M1加速度,第四章车辆系统动力学数值分析方法,M1速度,第四章车辆系统动力学数值分析方法,M1位移,第四章车辆系统动力学数值分析方法,M2位移,第四章车辆系统动力学数值分析方法,二、外加激扰（相同频率）,第四章车辆系统动力学数值分析方法,激扰（1）,第四章车辆系统动力学数值分析方法,加速度（M1）,第四章车辆系统动力学数值分析方法,加速度（M2）,第四章车辆系统动力学数值分析方法,M1位移,第四章车辆系统动力学数值分析方法,位移（质量2）,第四章车辆系统动力学数值分析方法,三、外加激扰（不同频率）,第四章车辆系统动力学数值分析方法,激扰,第四章车辆系统动力学数值分析方法,激扰（8）,第四章车辆系统动力学数值分析方法,加速度（质量1）,第四章车辆系统动力学数值分析方法,加速度（质量2）,第四章车辆系统动力学数值分析方法,加速度（质量3）,第四章车辆系统动力学数值分析方法,位移（质量1）,第四章车辆系统动力学数值分析方法,位移（质量2）,第四章车辆系统动力学数值分析方法,位移（质量3）,第四章车辆系统动力学数值分析方法,第四章车辆系统动力学数值分析方法,加速度（质量1）,第四章车辆系统动力学数值分析方法,角加速度（质量1）,第四章车辆系统动力学数值分析方法,加速度（质量2）,第四章车辆系统动力学数值分析方法,加速度（质量3）,第四章车辆系统动力学数值分析方法,位移（质量1）,第四章车辆系统动力学数值分析方法,角位移（质量1）,第四章车辆系统动力学数值分析方法,位移（质量2）,第四章车辆系统动力学数值分析方法,位移（质量3）,第四章车辆系统动力学数值分析方法,一、整车垂向动力学模型,课程报告,第四章车辆系统动力学数值分析方法,二、系统参数,第四章车辆系统动力学数值分析方法,三、系统激扰,第四章车辆