合肥工業(yè)大學(xué)-直動(dòng)式溢流閥仿真.doc

《合肥工業(yè)大學(xué)-直動(dòng)式溢流閥仿真.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、計(jì)算機(jī)仿真直動(dòng)式溢流閥仿真實(shí)驗(yàn)班級(jí)：機(jī)設(shè)13-4班學(xué)號(hào)：2013210360姓名：楊尚武授課教師：翟華日期：2016年4月13日一、實(shí)驗(yàn)?zāi)康谋緦?shí)驗(yàn)要求學(xué)生能掌握連續(xù)系統(tǒng)仿真的一般過程，狀態(tài)變量法的一般過程，鍵合圖法仿真的一般過程，以四階定步長(zhǎng)龍格——庫塔法的計(jì)算機(jī)程序的編寫。二、實(shí)驗(yàn)要求本實(shí)驗(yàn)要求學(xué)生能掌握連續(xù)系統(tǒng)仿真的一般過程，狀態(tài)變量法的一般過程，鍵合圖法仿真的一般過程，以及四階定步長(zhǎng)龍格—庫塔法的計(jì)算機(jī)程序的編號(hào)要求每個(gè)學(xué)生參考上述源程序，獨(dú)立編寫C語言源程序或其他高級(jí)語言程序，正確計(jì)算仿真結(jié)果，并繪制壓力P（t）和時(shí)間

2、t的關(guān)系圖，以及閥芯位移x（t）和時(shí)間t的關(guān)系。三、參考C語言程序#include"stdio.h"#include"math.h"voidmain(){FILE*fp;inti,j,e,g,l,s,n1,n2,b[3][2];floatd1,r1,r2,i1,c1,c2,x1,h1,a1,c,p0,pp,q1,x2,t;floata[3][3],y[3],u[2],k[3][5],h[5],p[3][5],z[3],d[3][5];if((fp=fopen("fz1","wb"))==NULL){printf("cannot

3、openfile");}fprintf(fp,"digitalsimulationofhydraulicreliefvalve");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;a[0][0]=-a1*a1*r1/i1;a[0][1]=-1.0/c2;a[0][2]=a1/

4、c1;a[1][0]=1/i1;a[1][1]=0.0;a[1][2]=0.0;a[2][0]=-a1/i1;a[2][1]=0.0;a[2][2]=-1/(r2*c1);for(i=0;i<=2;i++){for(j=0;j<=1;j++)b[i][j]=0;}b[0][0]=-1;b[2][1]=1;y[0]=0.0;y[1]=0.0;y[2]=p0*c1;x2=q1/(0.7*3.14*d1*sqrt(2.0*pp/900.0));u[0]=pp*a1-(x1+x2)/c2;u[1]=q1;c=-0.7*3.14*d1*

5、sqrt(2.0/900.0);for(i=0;i<=2;i++){k[i][0]=0.0;}h[0]=0.0;h[1]=h1/2.0;h[2]=h1/2.0;h[3]=h1;fprintf(fp,"TP(T)X(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++){p[i][j]=h[j-1]*k[i][j-1];z[i]=y[i];z[i]=z[i]+p[i][j];}if(z[1]<0.0)z[1]=

6、0.0;if(y[1]==0.0&&z[0]<0.0)z[0]=0.0;if(z[2]<0.0)z[2]=0.0;for(i=0;i<=2;i++){d[i][j]=0.0;for(l=0;l<=2;l++){d[i][j]=d[i][j]+a[i][l]*z[l];}k[i][j]=d[i][j];for(s=0;s<=1;s++){k[i][j]=k[i][j]+b[i][s]*u[s];}}if(y[1]==0.0&&k[0][j]<0.0)k[0][j]=0.0;if(y[1]==0.0&&k[1][j]<0.0)k[

7、1][j]=0.0;if(y[2]==0.0&&k[2][j]<0.0)k[2][j]=0.0;if(y[1]>x1){k[2][j]=k[2][j]+c*(z[1]-x1)*sqrt(z[2]/c1);if(y[2]==0.0&&k[2][j]<0.0)k[2][j]=0.0;}}for(i=0;i<=2;i++){y[i]=y[i]+h1*(k[i][1]+2*k[i][2]+2*k[i][3]+k[i][4])/6.0;}if(y[1]<0.0)y[1]=0.0;if(y[1]==0.0&&y[0]<0.0)y[0]=0

8、.0;if(y[2]<0.0)y[2]=0.0;}fprintf(fp,"%7.5f%e%e",t,y[2]/c1,y[1]);}}四、輸出結(jié)果digitalsimulationofhydraulicreliefvalveTP(T)X(T)0.001501.449418e