合肥工业大学-直动式溢流阀仿真.doc

计算机仿真直动式溢流阀仿真实验班 级： 机设13-4班 学 号： 2013210360 姓 名： 杨 尚 武 授课教师： 翟 华 日 期： 2016年4月13日 1、 实验目的 本实验要求学生能掌握连续系统仿真的一般过程，状态变量法的一般过程，键合图法仿真的一般过程，以四阶定步长龙格库塔法的计算机程序的编写。二、 实验要求本实验要求学生能掌握连续系统仿真的一般过程，状态变量法的一般过程，键合图法仿真的一般过程，以及四阶定步长龙格库塔法的计算机程序的编号要求每个学生参考上述源程序，独立编写C语言源程序或其他高级语言程序，正确计算仿真结果，并绘制压力P（t）和时间t的关系图，以及阀芯位移x（t）和时间t的关系。三、参考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.39e+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;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;fprintf(fp, T P(T) X(T)n);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) 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(y20.0)y2=0.0;fprintf(fp,%7.5f %e %en,t,y2/c1,y1);四、输出结果digital simulation of hydraulic relief valveT P(T) X(T)0.00150 1.449418e+06 0.000000e+000.00300 2.288071e+06 0.000000e+000.00450 3.099895e+06 1.151061e-040.00600 3.868414e+06 4.654126e-040.00750 4.600404e+06 1.006565e-030.00900 4.988608e+06 1.692020e-030.01050 3.161974e+06 2.095751e-030.01200 1.816280e+06 1.854696e-030.01350 1.982722e+06 1.433021e-030.01500 2.839527e+06 1.252615e-030.01650 3.644874e+06 1.364071e-030.01800 3.882667e+06 1.673995e-030.01950 3.016941e+06 1.839893e-030.02100 2.377060e+06 1.706981e-030.02250 2.531030e+06 1.497600e-030.02400 3.156087e+06 1.446740e-030.02550 3.501333e+06 1.583777e-030.02700 3.186811e+06 1.717583e-030.02850 2.736281e+06 1.695083e-030.03000 2.673098e+06 1.579036e-030.03150 2.970536e+06 1.515308e-030.03300 3.243553e+06 1.563448e-030.03450 3.179934e+06 1.648301e-030.03600 2.923158e+06 1.665780e-030.03750 2.805473e+06 1.611564e-030.03900 2.920318e+06 1.561015e-030.04050 3.096681e+06 1.568833e-030.04200 3.124609e+06 1.614838e-030.04350 3.001740e+06 1.639431e-030.04500 2.899607e+06 1.619705e-030.04650 2.925323e+06 1.587159e-030.04800 3.023885e+06 1.580377e-030.04950 3.071301e+06 1.601658e-030.05100 3.024275e+06 1.621452e-030.05250 2.955756e+06 1.617937e-030.05400 2.946312e+06 1.600020e-030.05550 2.993725e+06 1.590586e-030.05700 3.033829e+06 1.598347e-030.05850 3.023208e+06 1.611109e-030.06000 2.984370e+06 1.613578e-030.06150 2.966569e+06 1.605188e-030.06300 2.985009e+06 1.597548e-030.06450 3.011714e+06 1.598983e-030.06600 3.015217e+06 1.606013e-030.06750 2.996430e+06 1.609634e-030.06900 2.980998e+06 1.606533e-030.07050 2.985335e+06 1.601573e-030.07200 3.000486e+06 1.600654e-030.07350 3.007349e+06 1.603957e-030.07500 2.999997e+06 1.606921e-030.07650 2.989607e+06 1.606318e-030.07800 2.988369e+06 1.