基于MATLAB的直动式减压阀的P-Q特性曲线的仿真.docx

基于MATLAB的直动式减压阀的P-Q特性曲线的仿真一概述减压阀的作用是将系统压力减压、稳压的一种控制元件，其调节方式粉直动式、先导式两种。如图51 所示，直动式减压阀由主阀芯、膜片、弹簧、调节手柄、主阀体组成，其核心部件是主阀芯、膜片和调节弹簧。P1 为气源压力，P2 为减压输出压力，FS 为弹簧压紧力，当P2AFS 时，主阀芯向下移动，主阀口打开，P2 上升；当P2AFS 时，主阀芯上移，将主阀口关闭，膜片继续上移，气体会从膜片的泄压口溢流到大气中，使P2 下降；当输出压力降到调定压力FS/A 时，膜片上的受力会保持平衡状态。（1） 图51 直动式减压阀二减压阀的特性曲线减压阀的特性曲线主要为压力特性曲线、流量特性、溢流特性。（只对流量特性进行仿真）压力特性曲线是指在流量不变时，输入压力变化引起输出压力变化的特性曲线；流量特性曲线是指在输入压力不变时，输出流量变化引起输出压力变化的特性曲线；溢流特性曲线是指溢流量与输出压力间的关系，一般减压阀的溢流能力很小，不会在系统中当溢流阀使用，在绝大多数系统中，一般不用考虑减压阀的溢流特性曲线。理想的特性曲线如图52 所示。当P1 小于设定压力时，主阀芯全开，P2等于P1，当P1 大于设定压力时，理想的状态是P2 等于设定压力，并保持不变；理想的流量特性曲线为不管流经减压阀的流量如何变化，P2 应等于设定压力，并保持不变。实际上，由于减压阀各参数的相互影响，实际特性与理想特性相差甚远。减压阀的理想流量特性曲线三减压阀的工作原理和受力分析作用在主阀芯及膜片受力分析：“向上的力”有P1 作用在主阀芯底部的力，P2 作用在膜片上的向上的力；“向下的力”有弹簧力FS，P2 作用在主阀芯的向下的力假设作用在主阀芯的面积为A2，弹簧预紧压缩量为X0，弹簧刚度为K，阀口开度为X，则力学平衡方程如下：P1A2+P2A=KX0-X+P2A2当流量增加时，阀口开度增加，P2 会减小；当P1 增加时，为了保证流量不变，阀口开度会减小，P2 会增加。四减压阀工作时的动态性能和分析由于减压阀在工作时还有稳态液动力，因此有：P1A1+P2A2=KX0-X+P2A1+2CdxRcos（P1-P2）流经小孔的的流量为：q=Acvc=CdA02(P1-P2)小孔面积： A0=12（2r12cos-2r1tan-x2cossin）五通过matlab对减压阀进行性能仿真经过计算发现稳态液动力很小，可以忽略，于是有：P1A1+P2A2=KX0-X+P2A1q=Acvc=Cd2*2r12cos-2r1tan-x2cossin*2(P1-P2)其中，r1为A1的半径。设r1=0.01m，A2=0.005m2，P1=10MPa，=45，K=100kNm，=5kgm3通过matlab仿真后得出曲线如下图：Matlab计算程序：r1=0.01;%主阀芯底部半径A1=3.14*(r1)2;%主阀芯底部作用面积A2=0.005;%P2作用膜片面积（上部）K=1e5;%弹簧刚度P1=1e7;%入口处压力P2=1000:1000:1e7;%将P2从10000变到1000000，其中每次变化的数值为1000q=zeros(1,10000);%建立一个1维矩阵a=45;%主阀芯底部倾斜角Cd=0.7;for i=1:10000 x=(P1*A1+P2(i)*A2-P2(i)*A1)/K; q(i)=Cd/2*(2*3.14*r12/cos(a)-2*3.14*(r1*tan(a)-x)*cos(a)*sin(a)*sqrt(2*(P1-P2(i)/5); endplot(q,P2)axis(0 100 9000000 11000000)xlabel qylabel P2通过仿真可以看出，减压阀的流量特性曲线基本接近理论曲线。不过由于不知道一般减压阀的参数，不知道流量q和压力P是不是会太大。六作业感想气动这门课平时没有作业，只有这么一个算不小的作业，对我来说感触很深。首先我没有用过matlab，因此为了做这个作业，让我了解了基本的matlab的操作方法，其间花了不少时间和经历（单单安装