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

《车辆系统动力学计算机仿真.ppt》由会员分享，可在线阅读，更多相关《车辆系统动力学计算机仿真.ppt（56页珍藏版）》请在三一办公上搜索。

1、第四章 车辆系统动力学数值分析方法,车辆系统动力学计算机仿真,第四章 车辆系统动力学数值分析方法,一、模型,第四章 车辆系统动力学数值分析方法,二、参数,第四章 车辆系统动力学数值分析方法,三、激扰,第四章 车辆系统动力学数值分析方法,implicit double precision(a-h,o-z)parameter(nn=40000,nd=2,gg=9.81,pi=3.14159)parameter(th=1000.)COMMON/EXPL/step,F1KC,FKC,FI1,FIdimension zb1(nd),vb1(nd),ab1(nd)dimension dzb1(nd),dv

2、b1(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,第四章 车辆系统动力学数值分析方法,do i=1,nn,1 tt=dble(i)*st

3、ep zt2=0.05*sin(4.*pi*tt)zt3=0.05*sin(4.*pi*tt)CALL EXP1(dzb1,dvb1,dab1,d1ab1,zb1,vb1,nd)!积分 CALL EXP2(dzb2,dvb2,dab2,d1ab2,zb2,vb2)!积分 CALL EXP2(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

4、*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)then wri

5、te(11,100)tt,ab1(1),ab2,ab3,zt2 write(12,100)tt,zb1(1)*th,zb2*th,zb3*th,zt3 end if CALL EXCH1(zb1,vb1,ab1,dzb1,dvb1,dab1,d1ab1,ND)CALL EXCH2(zb2,vb2,ab2,dzb2,dvb2,dab2,d1ab2)CALL EXCH2(zb3,vb3,ab3,dzb3,dvb3,dab3,d1ab3)end do 100 format(1x,f12.5,2x,f12.5,2x,f12.5,2x,f12.5,2x,f12.5)stopend,第四章 车辆系统动力学

6、数值分析方法,SUBROUTINE EXP1(DY,DV,DA,D1A,CY,CV,ND)IMPLICIT DOUBLE PRECISION(A-H,O-Z)COMMON/EXPL/VST,F1KC,FKC,FI1,FIDIMENSION DY(ND),DV(ND),DA(ND),D1A(ND),CY(ND),CV(ND)VST2=VST*VSTDO 700 I=1,ND,1 CY(I)=DY(I)+DV(I)*VST+(F1KC*DA(I)-FKC*D1A(I)*VST2 CV(I)=DV(I)+(FI1*DA(I)-FI*D1A(I)*VST700 CONTINUERETURNEND,第四

7、章 车辆系统动力学数值分析方法,SUBROUTINE EXP2(DY,DV,DA,D1A,CY,CV)IMPLICIT DOUBLE PRECISION(A-H,O-Z)COMMON/EXPL/VST,F1KC,FKC,FI1,FIVST2=VST*VST CY=DY+DV*VST+(F1KC*DA-FKC*D1A)*VST2 CV=DV+(FI1*DA-FI*D1A)*VSTRETURNEND,第四章 车辆系统动力学数值分析方法,SUBROUTINE EXCH1(ZZC,VZC,AZC,DZZC,DVZC,DAZC,D1AZC,ND)IMPLICIT DOUBLE PRECISION(A-H

8、,O-Z)DIMENSION ZZC(ND),VZC(ND),AZC(ND)DIMENSION DZZC(ND),DVZC(ND),DAZC(ND),D1AZC(ND)DO 5200 I=1,ND,1 DZZC(I)=ZZC(I)DVZC(I)=VZC(I)D1AZC(I)=DAZC(I)DAZC(I)=AZC(I)ZZC(I)=0.VZC(I)=0.AZC(I)=0.5200 CONTINUERETURNEND,第四章 车辆系统动力学数值分析方法,SUBROUTINE EXCH2(ZZC,VZC,AZC,DZZC,DVZC,DAZC,D1AZC)IMPLICIT DOUBLE PRECISI

9、ON(A-H,O-Z)DZZC=ZZC DVZC=VZC D1AZC=DAZC DAZC=AZC ZZC=0.VZC=0.AZC=0.RETURNEND,第四章 车辆系统动力学数值分析方法,一、无外加激扰（自由衰减振动）,第四章 车辆系统动力学数值分析方法,M1加速度,第四章 车辆系统动力学数值分析方法,M1速度,第四章 车辆系统动力学数值分析方法,M1位移,第四章 车辆系统动力学数值分析方法,M2位移,第四章 车辆系统动力学数值分析方法,二、外加激扰（相同频率）,第四章 车辆系统动力学数值分析方法,激扰（1）,第四章 车辆系统动力学数值分析方法,加速度（M1）,第四章 车辆系统动力学数值分析

10、方法,加速度（M2）,第四章 车辆系统动力学数值分析方法,M1位移,第四章 车辆系统动力学数值分析方法,位移（质量2）,第四章 车辆系统动力学数值分析方法,三、外加激扰（不同频率）,第四章 车辆系统动力学数值分析方法,激扰,第四章 车辆系统动力学数值分析方法,激扰（8）,第四章 车辆系统动力学数值分析方法,加速度（质量1）,第四章 车辆系统动力学数值分析方法,加速度（质量2）,第四章 车辆系统动力学数值分析方法,加速度（质量3）,第四章 车辆系统动力学数值分析方法,位移（质量1）,第四章 车辆系统动力学数值分析方法,位移（质量2）,第四章 车辆系统动力学数值分析方法,位移（质量3）,第四章 车

11、辆系统动力学数值分析方法,第四章 车辆系统动力学数值分析方法,加速度（质量1）,第四章 车辆系统动力学数值分析方法,角加速度（质量1）,第四章 车辆系统动力学数值分析方法,加速度（质量2）,第四章 车辆系统动力学数值分析方法,加速度（质量3）,第四章 车辆系统动力学数值分析方法,位移（质量1）,第四章 车辆系统动力学数值分析方法,角位移（质量1）,第四章 车辆系统动力学数值分析方法,位移（质量2）,第四章 车辆系统动力学数值分析方法,位移（质量3）,第四章 车辆系统动力学数值分析方法,一、整车垂向动力学模型,课程报告,第四章 车辆系统动力学数值分析方法,二、系统参数,第四章 车辆系统动力学数值分析方法,三、系统激扰,第四章 车辆系统动力学数值分析方法,第四章 车辆系统动力学数值分析方法,车体:,第四章 车辆系统动力学数值分析方法,构架:,第四章 车辆系统动力学数值分析方法,二系力,一系力,第四章 车辆系统动力学数值分析方法,轮轨力：,车体：垂向、侧滚(2)构架：垂向、侧滚、点头（3）轮对：垂向、侧滚（4）,系统自由度：,第四章 车辆系统动力学数值分析方法,二、系统参数,