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

第四章 车辆系统动力学数值分析方法,车辆系统动力学计算机仿真,第四章 车辆系统动力学数值分析方法,一、模型,第四章 车辆系统动力学数值分析方法,二、参数,第四章 车辆系统动力学数值分析方法,三、激扰,第四章 车辆系统动力学数值分析方法,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,fi dimension zb1(nd),vb1(nd),ab1(nd) dimension dzb1(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+4 cbk=1.0e+7;dbp=2.0e+4 bm1=3250.0 bm2=2100.0 bim=3020.0 bl=2.5,第四章 车辆系统动力学数值分析方法,f1kc=1.0 fkc=0.5 fi1=1.5 fi=0.5 step=1.0e-4,第四章 车辆系统动力学数值分析方法,do i=1,nn,1 tt=dble(i)*step 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*bl*vb1(2),第四章 车辆系统动力学数值分析方法,fm2=cmk*(dm2-zb2)+dmp*(vm2-vb2) fm3=cmk*(dm3-zb3)+dmp*(vm3-vb3) fb2=cbk*(zb2-zt2)+dbp*vb2 fb3=cbk*(zb3-zt3)+dbp*vb3,第四章 车辆系统动力学数值分析方法,ab1(1)=(-fm2-fm3)/bm1-gg ab1(2)=0.5*bl*(fm3-fm2)/bim ab2=(-fb2+fm2)/bm2-gg ab3=(-fb3+fm3)/bm2-gg,第四章 车辆系统动力学数值分析方法,if(mod(nn,800).eq.0)then write(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) stop end,第四章 车辆系统动力学数值分析方法,subroutine exp1(dy,dv,da,d1a,cy,cv,nd) implicit double precision(a-h,o-z) common /expl/ vst,f1kc,fkc,fi1,fi dimension dy(nd),dv(nd),da(nd),d1a(nd),cy(nd),cv(nd) vst2=vst*vst do 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)*vst 700 continue return end,第四章 车辆系统动力学数值分析方法,subroutine exp2(dy,dv,da,d1a,cy,cv) implicit double precision(a-h,o-z) common /expl/ vst,f1kc,fkc,fi1,fi vst2=vst*vst cy=dy+dv*vst+(f1kc*da-fkc*d1a)*vst2 cv=dv+(fi1*da-fi*d1a)*vst return end,第四章 车辆系统动力学数值分析方法,subroutine exch1(zzc,vzc,azc,dzzc,dvzc,dazc,d1azc,nd) implicit double precision(a-h,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 continue return end,第四章 车辆系统动力学数值分析方法,subroutine exch2(zzc,vzc,azc,dzzc,dvzc,dazc,d1azc) implicit double precision(a-h,o-z) dzzc=zzc dvzc=vzc d1azc=dazc dazc=azc zzc=0. vzc=0. azc=0. return end,第四章 车辆系统动力学数值分析方法,一、无外加激扰（自由衰减振动）,第四章 车辆系统动力学数值分析方法,m1加速度,第四章 车辆系统动力学数值分析方法,m1速度,第四章 车辆系统动力学数值分析方法,m1位移,第四章 车辆系统动力学数值分析方法,m2位移,第四章 车辆系统动力学数值分析方法,二、外加激扰（相同频率）,第四章 车辆系统动力学数值分析方法,激扰（1）,第四章 车辆系统动力学数值分析方法,加速度（m1）,第四章 车辆系统动力学数值分析方法,加速度（m2）,第四章 车辆系统动力学数值分析方法,m1位移,第四章 车辆系统动力学数值分析方法,位移（质量2）,第四章 车辆系统动力学数值分析方法,三、外加激扰（不同频率）,第四章 车辆系统动力学数值分析方法,激扰,第四章 车辆系统动力学数值分析方法,激扰（8）,第四章 车辆系统动力学数值分析方法,加速度（质量1）,第四章 车辆系统动力学数值分析方法,加速度（质量2）,第四章 车辆系统动力学数值分析方法,加速度（质量3）,第四章 车辆系统动力学数值分析方法,位移（质量1）,第四章 车辆系统动力学数值分析方法,位移（质量2）,第四章 车辆系统动力学数值分析方法,位移（质量3）,第四章 车辆系统动力学数值分析方法,第四章 车辆系统动力学数值分析方法,加速度（质量1）,第四章 车辆系统动力学数值分析方法,角加速度（质量1）,第四章 车辆系统动力学数值分析方法,加速度（质量2）,第四章 车辆系统动力学数值分析方法,加速度（质量3）,第四章 车辆系统动力学数值分析方法,位移（质量1）,第四章 车辆系统动力学数值分析方法,角位移（质量1）,第四章 车辆系统动力学数值分析方法,位移（质量2）,第四章 车辆系统动力学数值分析方法,位移（质量3）,第四章 车辆系统动力学数值分析方法,一、整车垂向动力学模型,课程报告,第四章 车辆系统动力学数值分析方法,二、系统参数,第四章 车辆系统动力学数值分析方法,三、系统激扰,第四章 车