《合肥工业大学》由会员分享，可在线阅读，更多相关《合肥工业大学（7页珍藏版）》请在金锄头文库上搜索。
1、计算机仿真直动式溢流阀仿真实验班 级：机设13-4班学 号：2013210360姓 名：杨尚武授课教师：翟华日 期：2016年4月13日一、实验目的本实验要求学生能掌握连续系统仿真的一般过程，状态变量法的一般过程，键合 图法仿真的一般过程，以四阶定步长龙格库塔法的计算机程序的编写。二、 实验要求本实验要求学生能掌握连续系统仿真的一般过程，状态变量法的一般过程， 键合图法仿真的一般过程，以及四阶定步长龙格库塔法的计算机程序的编号 要求每个学生参考上述源程序，独立编写 C 语言源程序或其他高级语言程序，正 确计算仿真结果，并绘制压力P (t)和时间t的关系图，以及阀芯位移x (t) 和时间 t 的
2、关系。三、参考C语言程序#include stdio.h#include math.hvoid main()FILE * fp;int i,j,e,g,l,s,n1,n2,b32;float d1,r1,r2,i1,c1,c2,x1,h1,a1,c,p0,pp,q1,x2,t;float a33,y3,u2,k35,h5,p35,z3,d35; if(fp=fopen(fz1,wb)=NULL)printf(cannot open filen);fprintf(fp,digital simulation of hydraulic relief valven); d1=0.012;r1=0.39
3、e+11;r2=0.147e+12;i1=0.0614;c1=0.8e-12;c2=0.2e-4;x1=0.0014;p0=0.6e+6;pp=0.3e+7;q1=0.46e-3;h1=0.0003;t=0.0;n1=5;n2=55;a1=0.785*d1*d1;a00=-a1*a1*r1/i1; a01=-1.0/c2;a02=a1/c1; a10=1/i1;a11=0.0;a12=0.0; a20=-a1/i1;a21=0.0; a22=-1/(r2*c1);for (i=0;i=2;i+)for (j=0;j=1;j+) bij=0;b00=-1;b21=1;y0=0.0;y1=0.0;
4、y2=p0*c1;x2=q1/(0.7*3.14*d1*sqrt(2.0*pp/900.0); u0=pp*a1-(x1+x2)/c2;u1=q1;c=-0.7*3.14*d1*sqrt(2.0/900.0); for(i=0;i=2;i+)ki0=0.0;h0=0.0;h1=h1/2.0;h2=h1/2.0;h3=h1;X(T)n);fprintf(fp, T P(T)for (e=1;e=n2;e+)for (g=1;g=n1;g+)t=t+h1;for(j=1;j=4;j+) for(i=0;i=2;i+) pij=hj-1*kij-1; zi=yi;zi=zi+pij;if(z10.0
5、) z1=0.0; if(y1=0.0&z00.0)z0=0.0;if(z20.0) z2=0.0;for (i=0;i=2;i+)dij=0.0;for(l=0;l=2;l+)dij=dij+ail*zl; kij=dij;for(s=0;s=1;s+) kij=kij+bis*us; if(y1=0.0&k0j0.0)k0j=0.0;if(y1=0.0&k1j0.0)k1j=0.0;if(y2=0.0&k2jx1)k2j=k2j+c*(z1-x1)*sqrt(z2/c1);if(y2=0.0&k2j0.0) k2j=0.0;for(i=0;i=2;i+)yi=yi+h1*(ki1+2*ki2+2*ki3+ki4)/6.0; if(y10.0)y1=0.0;if(y1=0.0&y00.0)y0=0.0;if(y2- -C9 3 -M 弋业壬型图_匕屯4.阀芯位移随时间变化的曲线为：:1虫丄$ 、h -爭玄丘 d El田 33
