汽车空气悬架建模与仿真分析

#l6l6reducevehiclevibrationfrequencyandwheeldynamicload,makeitswingoodridecomfortandsex,handlingstabilityandhigh-speeddrivingsafety,reducingvehicleofautomobidamage.Thispaperautomobileairspringsuspensionandairsuspensionsystemcharacteristic,theinfluenceofvehiclethroughasimplecalculation,andthenthecalculationresultsintothecomputerandthegraphicbycomputersimulation,wewillgetresults.Keywords:theairsuspension,airspring,modeling,computersimulation附录编程模块grktl.c#include"stdlib.h"voidgrktl(t,y,n,h,k,z)intn,k;doublet,h,y[],z[];{externvoidgrktlf();inti,j,l;doublea[4],tt,*b,*d;b=malloc(n*sizeof(double));d=malloc(n*sizeof(double));a[0]=h/2.0;a[l]=a[0];a[2]=h;a[3]=h;for(i=0;i<=n-l;i++)z[i*k]=y[i];for(l=l;l<=k-l;l++){grktlf(t,y,n,d);for(i=0;i<=n-l;i++)b[i]=y[i];for(j=0;j<=2;j++){for(i=0;i<=n-l;i++){y[i]=z[i*k+l-l]+a[j]*d[i];b[i]=b[i]+a[j+l]*d[i]/3.0;}tt=t+a[j];grktlf(tt,y,n,d);for(i=0;i<=n-1;i++)y[i]=b[i]+h*d[i]/6.0;for(i=0;i<=n-1;i++)z[i*k+l]=y[i];t=t+h;}free(b);free(d);return;}汽车车身速度的编程程序：#include"stdio.h"#include"grkt1.c"#include"math.h"main(){inti,j;FILE*fp;doublet,h,y[4],z[4][101];y[0]=0.0;y[1]=0.0;y[2]=0.0;y[3]=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf("\n");if((fp=fopen("test3.txt","w"))==NULL){printf("cannotopenfile\n");exit(0);}for(i=0;i<=100;i++){t=i*h;printf("t=%5.2f\n",t);fprintf(fp,"%6.2f\n",z[3][i]);for(j=0;j<=3;j++)printf("y(%d)=%e",j,z[j][i]);printf("\n");}}voidgrkt1f(t,y,n,d)intn;doublet,y[],d[];{t=t;n=n;d[0]=y[1];d[1]=-200*(y[1]-y[3])-150*(y[0]-y[2]);d[2]=y[3];d[3]=-100*(y[3]-y[1])+250*(y[2]-y[0])-4000*(y[2]-sin(t));return;}汽车的车身加速度的编程：#include"stdio.h"#include"grkt1.c"#include"math.h"main(){inti,j;FILE*fp;doublet,h,y[4],z[4][101];y[0]=0.0;y[1]=0.0;y[2]=0.0;y[3]=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf("\n");if((fp=fopen("test5.txt","w"))==NULL){printf("cannotopenfile\n");exit(0);}for(i=0;i<=100;i++){t=i*h;printf("t=%5.2f\n",t);fprintf(fp,"%6.2f\n",(z[3][i]-z[3][I-1])/0.01);for(j=0;j<=3;j++)printf("y(%d)=%e",j,z[j][i]);printf("\n");}}voidgrkt1f(t,y,n,d)intn;doublet,y[],d[];{t=t;n=n;d[0]=y[1];d[1]=-200*(y[1]-y[3])-150*(y[0]-y[2]);d[2]=y[3];d[3]=-100*(y[3]-y[1])+250*(y[2]-y[0])-4000*(y[2]-sin(t));return;}汽车车轮速度的编程#include"stdio.h"#include"grkt1.c"#include"math.h"main(){inti,j;FILE*fp;doublet,h,y[4],z[4][101];y[0]=0.0;y[1]=0.0;y[2]=0.0;y[3]=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf("\n");if((fp=fopen("test2.txt","w"))==NULL){printf("cannotopenfile\n");exit(0);}for(i=0;i<=100;i++){t=i*h;printf("t=%5.2f\n",t);fprintf(fp,"%6.2f\n",z[1][i]);for(j=0;j<=3;j++)printf("y(%d)=%e",j,z[j][i]);printf("\n");}}voidgrkt1f(t,y,n,d)intn;doublet,y[],d[];{t=t;n=n;d[0]=y[1];d[1]=-200*(y[1]-y[3])-150*(y[0]-y[2]);d[2]=y[3];d[3]=-100*(y[3]-y[1])+250*(y[2]-y[0])-4000*(y[2]-sin(t));return;}汽车车轮的加速度编程#include"stdio.h"#include"grkt1.c"#include"math.h"main(){inti,j;FILE*fp;doublet,h,y[4],z[4][101];y[0]=0.0;y[1]=0.0;y[2]=0.0;y[3]=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf("\n");if((fp=fopen("test6.txt","w"))==NULL){printf("cannotopenfile\n");exit(0);}for(i=0;i<=100;i++){t=i*h;printf("t=%5.2f\n",t);fprintf(fp,"%6.2f\n",(z[1][i]-z[1][i-1])/0.01);for(j=0;j<=3;j++)printf("y(%d)=%e",j,z[j][i]);printf("\n");}}voidgrkt1f(t,y,n,d)intn;doublet,y[],d[];{t=t;n=n;d[0]=y[1];d[1]=-200*(y[1]-y[3])-150*(y[0]-y[2]);d[2]=y[3];d[3]=-100*(y[3]-y[1])+250*(y[2]-y[0])-4000*(y[2]-sin(t));return;}汽车车身的位移的编程#include"stdio.h"#include"grkt1.c"#include"math.h"main(){inti,j;FILE*fp;doublet,h,y[4],z[4][101];y[0]=0.0;y[1]=0.0;y[2]=0.0;y[3]=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf("\n");if((fp=fopen("test3.txt","w"))==NULL){printf("cannotopenfile\n");exit(0);}for(i=0;i<=100;i++){t=i*h;printf("t=%5.2f\n",t);fprintf(fp,"%6.2f\n",z[2][i]);for(j=0;j<=3;j++)printf("y(%d)=%e",j,z[j][i]);printf("\n");}}voidgrkt1f(t,y,n,d)intn;doublet,y[],d[];{t=t;n=n;d[0]=y[1];d[1]=-200*(y[1]-y[3])-150*(y[0]-y[2]);d[2]=y[3];d[3]=-100*(y[3]-y[1])+250*(y[2]-y[0])-4000*(y[2]-sin(t));return;}汽车车轮的位移的编程#include"stdio.h"#include"grkt1.c"#include"math.h"main(){inti,j;FILE*fp;doublet,h,y[4],z[4][101];y[0]=0.0;y[1]=0.0;y[2]=0.0;y[3]=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf("\n");if((fp=fopen("test1.txt","w"))==NULL){printf("cannotopenfile\n");exit(0);}for(i=0;i<=100;i++){t=i*h;printf("t=%5.2f\n",t);fprintf(fp,"%6.2f\n",z[0][i]);for(j=0;j<=3;j++)printf("y(%d)=%e",j,z[j][i]);printf("\n");}}voidgrkt1f(t,y,n,d)intn;doublet,y[],d[];{t=t;n=n;d[0]=y[1];d[1]=-200*(y[1]-y[3])-150*(y[0]-y[2]);d[2]=y[3];d[3]=-100*(y[3]-y[1])+250*(y[2]-y[0])-4000*(y[2]-sin(t));return;}Mabtlab车轮位移绘图程序：fip=fopen(‘test1.txt','r');y=fscanf(fip,'%g');fclose(fip);t=0:0.01:1.00;plot(t,y);gridon;车轮速度绘图程序fip=fopen(‘test2.txt','r');y=fscanf(fip,'%g');fclose(fip);t=0:0.01:1.00;plot(t,y);gridon;汽车车身位移绘图程序fip=fopen(‘test3.txt','r');y=fscanf(fip,'%g');fclose(fip);t=0:0.01:1.00;plot(t,y);gridon;汽车车身速度绘图程序fip=fopen(‘test4.txt','r');y=fscanf(