CN114396416B - 阀后压力补偿系统参数灵敏度分析方法 - Google Patents
阀后压力补偿系统参数灵敏度分析方法 Download PDFInfo
- Publication number
- CN114396416B CN114396416B CN202210015137.9A CN202210015137A CN114396416B CN 114396416 B CN114396416 B CN 114396416B CN 202210015137 A CN202210015137 A CN 202210015137A CN 114396416 B CN114396416 B CN 114396416B
- Authority
- CN
- China
- Prior art keywords
- valve
- power
- pressure compensation
- variable
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F15—FLUID-PRESSURE ACTUATORS; HYDRAULICS OR PNEUMATICS IN GENERAL
- F15B—SYSTEMS ACTING BY MEANS OF FLUIDS IN GENERAL; FLUID-PRESSURE ACTUATORS, e.g. SERVOMOTORS; DETAILS OF FLUID-PRESSURE SYSTEMS, NOT OTHERWISE PROVIDED FOR
- F15B19/00—Testing; Calibrating; Fault detection or monitoring; Simulation or modelling of fluid-pressure systems or apparatus not otherwise provided for
- F15B19/007—Simulation or modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Fluid-Pressure Circuits (AREA)
Abstract
本发明公开一种阀后压力补偿系统参数灵敏度分析方法,其包括以下步骤,步骤一:构建阀后压力补偿系统的功率键合图模型;步骤二:推导阀后压力补偿系统的状态方程;步骤三:推导阀后压力补偿系统的参数灵敏度微分方程,确定参数灵敏度。本发明提出一种阀后压力补偿系统参数灵敏度分析方法,该方法基于功率键合图理论建立阀后压力补偿系统数学模型,以解决阀后压力补偿系统的灵敏度分析问题,目的在于提供一种针对阀后压力补偿系统的参数灵敏度分析方法,为分析压力补偿阀中影响负载敏感比例多路阀动态特性的关键参数以及为优化负载敏感比例多路阀动态特性和压力补偿阀设计参数奠定研究基础和理论依据。
Description
技术领域
本申请涉及压力补偿阀技术领域,特别涉及一种阀后压力补偿系统参数灵敏度分析方法。
背景技术
压力补偿阀作为负载敏感比例多路阀中的核心元件,保证各个执行机构动作所需的流量按比例进行分配,实现各个执行机构的运动速度不受负载的影响,对整个阀的功能乃至整个系统的工作性能都起着举足轻重的作用。采用灵敏度分析方法探究压力补偿阀中影响负载敏感比例多路阀动态特性的关键参数,为优化负载敏感比例多路阀动态特性以及压力补偿阀设计参数奠定研究基础和理论依据。
目前,为了研究阀后压力补偿系统中的各项参数对负载敏感比例多路阀的影响,主要有两种方法,一是采用AMEsim建立仿真模型对其进行分析。二是对其中某一参数建立相关动态数学模型,通过Matlab对其进行分析。而本发明基于理论知识,可以对系统中所有参数进行分析,更加全面、方便、快捷,同时为优化压力补偿阀设计参数以及负载敏感比例多路阀动态特性奠定研究基础和理论依据。
发明内容
为了克服现有技术的不足,本发明提出一种阀后压力补偿系统参数灵敏度分析方法,该方法基于功率键合图建立阀后压力补偿系统数学模型,以解决阀后压力补偿系统的灵敏度分析问题,目的在于提供一种针对阀后压力补偿系统的灵敏度分析方法,为分析压力补偿阀中影响负载敏感比例多路阀动态特性的关键参数以及为优化负载敏感比例多路阀动态特性和压力补偿阀设计参数奠定研究基础和理论依据。
为实现上述目的,本发明所采用的解决方案为提供一种阀后压力补偿系统参数灵敏度分析方法,其包括以下步骤:
步骤1:构建阀后压力补偿系统的功率键合图模型;
根据负载敏感比例多路阀的工作原理和压力补偿阀的功能,基于功率键合图理论,搭建阀后压力补偿系统的功率键合图模型;所述的压力补偿系统包括:泵-主阀节流口模块、主阀节流口-压力补偿阀模块、压力补偿阀-LS管路模块、压力补偿阀-主阀换向模块和主阀换向-负载溢流阀模块五个模块;
步骤2:推导阀后压力补偿系统的状态方程;
根据步骤1建立的阀后压力补偿系统功率键合图模型,推导系统状态方程,则阀后压力补偿系统状态方程如下式所示:
式中:为液压泵至主阀节流口管路压缩的油液体积的一阶导数;为压力补偿阀阀芯动量的一阶导数;为压力补偿阀入口管路压缩的油液体积的一阶导数;为压力补偿阀出口至LS管路受到压缩的油液体积的一阶导数;为压力补偿阀出口至主阀换向管路压缩的油液体积的一阶导数;为主阀出口至负载溢流阀管路压缩的油液体积的一阶导数;Sf为液压泵恒定输出流量;Se1为PLS压力;Se2为负载溢流阀调定压力;R1为泵的泄漏液阻;R2为主阀芯节流口液阻;R3为压力补偿阀阻尼孔液阻;R4为LS管路液阻;R5为压力补偿阀阀口液阻;R6为主阀换向阀口液阻;R7为主阀出口至负载管路液阻;C1为泵至主阀节流口管路液容;C2为压力补偿阀入口管路液容;C3为压力补偿阀出口至LS管路液容;C4为压力补偿阀出口至主阀换向部分管路液容;C5为主阀出口至负载溢流阀管路液容;F1为压力补偿阀阀芯所受液动力;F2为压力补偿阀阀芯所受摩擦力;I1为压力补偿阀阀芯等效质量;A为压力补偿阀阀芯承压面积;A1为主阀节流口过流面积;A2为压力补偿阀阀口过流面积;A3为压力补偿阀阻尼孔过流面积;A4为主阀换向部分阀口过流面积;P10为压力补偿阀阀芯动量;V2为液压泵至主阀节流口管路压缩的油液体积;V13为压力补偿阀入口管路压缩的油液体积;V17为压力补偿阀出口至LS管路受到压缩的油液体积;V24为压力补偿阀出口至主阀换向管路压缩的油液体积;V28为主阀出口至负载溢流阀管路压缩的油液体积;Cd为流量系数;ρ为油液密度;
步骤3:推导阀后压力补偿系统的参数灵敏度微分方程,确定参数灵敏度;
步骤31:简化系统的状态方程;
由于参数灵敏度微分方程是针对状态方程研究系统参数变化引起的状态矢量变化,获取步骤2求得的系统的状态方程,将其简化成下式所示的形式:
步骤32:求解状态方程,确定参数灵敏度微分方程;
状态方程的解为:
其中,已知稳态液动力和摩擦力对阀动态特性的影响程度较小,主要分析参数矢量α中其余各参数灵敏度,状态矢量x对参数α的灵敏度函数定义为:
式中:i的取值为1,2,3…,d;λi表示第i个参数矢量灵敏度;αi表示第i个参数变量;
状态矢量x是关于输入矢量u与参数矢量α的函数,在u与α相互独立的情况下,在状态方程的等式两边同时对参数α求偏导,可得:
由于u与α相互独立,则可得:
最终得到系统参数灵敏度微分方程如下式所示:
步骤33:确定阀后压力补偿系统的参数灵敏度;
一阶轨迹参数灵敏度函数的初始条件为:
初始时刻,阀后压力补偿系统状态矢量x的初值x0=0,根据式可知:
根据步骤32计算得到的一阶轨迹灵敏度微分方程,将一阶轨迹灵敏度函数初值代入,即可确定阀后压力补偿系统参数灵敏度。
可优选的是,所述步骤1中的阀后压力补偿系统具体包括以下五个模块;
所述的泵-主阀节流口模块:根据功率键合图基本原理及规则,分析高负载联泵至主阀节流口部分的动态过程中的功率流程;恒流源代表负载所需要的流量,输出的液流功率在相同压力下分为三条支流:泵泄露的损失功率,其泄漏液阻为R1、泵至主阀管路油液压缩的损失功率,其管道液容为C1、其余流向主阀节流口;流向主阀节流口输出的功率在相同流量情况下分为两条支流:主阀节流口的损失功率,其主阀节流口液阻为R2、其余流向压力补偿阀;
所述的主阀节流口-压力补偿阀模块:根据功率键合图基本原理及规则,分析高负载联主阀口至压力补偿阀的动态过程中的功率流程;流向压力补偿阀输出的功率由于压力补偿阀开启需进行液压能与机械能转换,输出的功率可以分为四个支流:压力补偿阀阀芯所受液动力为F1、压力补偿阀阀芯所受惯性力,阀芯等效质量为I、压力补偿阀阀芯所受摩擦力为F2、其余机械能再次转化成液压能;压力补偿阀阀芯开启后,输入的功率流可以分为三个支流:压力补偿阀入口管路油液压缩的损失功率,管道液容为C2、流向压力补偿阀阻尼孔、流向压力补偿阀阀口;
所述的压力补偿阀-LS管路模块:根据功率键合图基本原理及规则,分析高负载联压力补偿阀至LS管路的动态过程中的功率流程;流经压力补偿阀入口输出的一部分功率可以分为两个支流:阻尼孔损失的功率,阻尼孔液阻为R3、其余流向LS管路;流向LS管路的功率可以分为两个支流:阻尼孔至LS口管路油液压缩的损失功率,管路液容为C3、其余流向LS口;流向LS口的功率可以分为两个支流:LS管路的损失功率,LS管路液阻为R4、LS口压力Se1;
所述的压力补偿阀-主阀换向模块:根据功率键合图基本原理及规则,分析高负载联压力补偿阀至主阀换向部分的动态过程中的功率流程;流经压力补偿阀入口输出的一部分功率可以分为两个支流:压力补偿阀阀口的损失功率,压力补偿阀阀口液阻为R5、其余流向主阀换向部分管路;流向主阀换向部分管路的功率可以分为两个支流:压力补偿阀阀口至主阀换向管路油液压缩损失的功率,管路液容为C4、其余流向主阀阀口;
所述的主阀换向-负载溢流阀模块:根据功率键合图基本原理及规则,分析高负载联主阀至负载溢流阀的动态过程中的功率流程;流向主阀阀口的功率可以分为两个支流:主阀阀口的损失功率,主阀口液阻为R6、其余流向负载管路;流向负载管路的功率可以分为两个支流:主阀阀口至负载管路油液压缩的损失功率,管路液容为C5、其余流向负载;流向负载的功率可以分为两个支流:主阀阀口至负载管路的损失功率,管路液阻为R7、负载压力Se2;
所述阀后压力补偿系统的功率键合图模型是由上述各模块得到的功率键合图模型连接组成。
可优选的是,所述步骤2中推导系统状态方程具体过程如下:
所述的压力补偿系统包括容性元和感性元;取容性元C的流变量q和感性元I的力变量p作为状态变量;为了便于建立状态方程,取容性元C和感性元I功率键上自变量对时间的积分为状态变量,即取液体体积V,固体或液体的动量P为状态变量;
液压系统中,对于容性元C,其力变量p与流变量q之间的关系如下式所示:
式中:p是容性元的力变量;q是容性元的流变量;dt表示时间t是q的积分变量;V为液体的油液体积;C表示容性元;
机械系统中,对于感性元I,其流变量与力变量之间的关系如下式所示:
式中:v表示感性元的流变量;F表示感性元的力变量;I表示感性元;
则根据阀后压力补偿系统功率键合图容性元件上力变量p与状态变量V、感性元件上流变量v与状态变量P的关系如下式所示:
式中:p2为功率键合图中第2个功率键上的力变量;v10为功率键合图中第10个功率键上的流变量;p13为功率键合图中第13个功率键上的力变量;p17为功率键合图中第17个功率键上的力变量;p24为功率键合图中第24个功率键上的力变量;p28为功率键合图中第28个功率键上的力变量;
阀后压力补偿系统功率键合图阻性元件上力变量与流变量的关系如下:
p3=R1q3
式中:p3为阻性元件R1上的力变量,q3为阻性元件R1上的流变量;
阀后压力补偿系统功率键合图变换器TF上力变量与力变量、流变量与流变量的关系如下:
p6·A=p7
式中:p6为功率键合图中第6个功率键上的力变量;p7为功率键合图中第7个功率键上的力变量;q6为功率键合图中第6个功率键上的流变量;q7为功率键合图中第7个功率键上的流变量;
阀后压力补偿系统功率键合图共势结,即0结,上力变量与力变量、流变量与流变量的关系如下:
p1=p2=p3=p4
q1=q2+q3+q4
式中:p1为功率键合图中第1个功率键上的力变量;p2为功率键合图中第2个功率键上的力变量;p3为功率键合图中第3个功率键上的力变量;p4为功率键合图中第4个功率键上的力变量;q1为功率键合图中第1个功率键上的流变量;q2为功率键合图中第2个功率键上的流变量;q3为功率键合图中第3个功率键上的流变量;q4为功率键合图中第4个功率键上的流变量;
阀后压力补偿系统功率键合图共流结,即1结,上力变量与力变量、流变量与流变量的关系如下:
q4=q5=q6
p4=p5+p6
式中:q4为功率键合图中第4个功率键上的流变量;q5为功率键合图中第5个功率键上的流变量;q6为功率键合图中第6个功率键上的流变量;p4为功率键合图中第4个功率键上的力变量;p5为功率键合图中第5个功率键上的力变量;p6为功率键合图中第6个功率键上的力变量;
根据功率键合图中各变量之间的映射关系及0-1节点的意义,即可推导得到状态方程。
可优选的是,所述步骤31中的状态矢量、输入矢量、参数矢量具体如下所示:
所述的状态矢量具体包括内容如下式所示:
x=[x1,x2,x3,x4,x5,x6]T=[V2,P10,V13,V17,V24,V28]T
式中:x1,x2,…x6分别表示第1、2…6个状态变量;
所述的输入矢量具体包括内容如下式所示:
u=[u1,u2,u3]T=[Sf,Se1,Se2]T
式中:u1,u2,u3分别表示第1、2、3个输入变量;
所述的参数矢量具体包括内容如下式所示:
α=[α1,α2,α3,α4,α5,α6,α7,α8,α9,α10,α11,α12,α13,α14]T
=[R1 -1,R4 -1,R7 -1,A,A1 -1,A2,A3,A4,I1 -1,C1 -1,C2 -1,C3 -1,C4 -1,C5 -1]T
式中:α1,α2,…α14分别表示第1、2…14个参数变量。
可优选的是,所述步骤32中的一阶轨迹灵敏度方程的系数项矩阵如下所示:
式中:a11,…a66分别表示一阶轨迹灵敏度方程的系数项矩阵的36个元素;
所述系数矩阵中各元素如下所示:
式中:b1,1,…b6,14分别表示一阶轨迹灵敏度方程的自由项矩阵的84个元素;
所述自由项矩阵中各元素如下式所示:
与现有技术相比,本发明的有益效果在于:
(1)本发明基于功率键合图理论建立了多路阀阀后压力补偿系统的数学模型,解决了阀后压力补偿系统的灵敏度参数计算的问题;
(2)本发明为分析压力补偿阀中影响负载敏感比例多路阀动态特性的关键参数以及为优化负载敏感比例多路阀动态特性和压力补偿阀设计参数奠定研究基础和理论依据。
附图说明
图1为本发明实施例阀后压力补偿系统参数灵敏度分析方法的控制框图;
图2为本发明实施例负载敏感比例多路阀第一联工作原理图;
图3为本发明实施例阀后压力补偿系统功率键合图;
图4a~4n为灵敏度参数分别为λ1…λ14时实施例阀后压力补偿系统参数一阶轨迹参数灵敏度的计算结果;
具体实施方式
以下,参照附图对本发明的实施方式进行说明。
本发明提出一种阀后压力补偿系统参数灵敏度分析方法,为分析压力补偿阀中影响负载敏感比例多路阀动态特性的关键参数以及为优化负载敏感比例多路阀动态特性和压力补偿阀设计参数奠定研究基础和理论依据。为了证明本发明的适用性,将其用于某阀后压力补偿系统参数灵敏度计算的实际计算过程,如图1所示为阀后压力补偿系统参数灵敏度分析方法的控制框图;其具体包括以下步骤:
S1:构建阀后压力补偿系统的功率键合图模型;
根据如图2所示的某负载敏感比例多路阀的工作原理、图3所示的多路阀第一联工作原理图和压力补偿阀的功能,基于功率键合图理论,搭建阀后压力补偿系统的功率键合图模型;压力补偿系统包括:泵-主阀节流口模块、主阀节流口-压力补偿阀模块、压力补偿阀-LS管路模块、压力补偿阀-主阀换向模块和主阀换向-负载溢流阀模块五个模块。
泵-主阀节流口模块:根据功率键合图基本原理及规则,分析高负载联泵至主阀节流口部分的动态过程中的功率流程;恒流源代表负载所需要的流量,输出的液流功率在相同压力下分为三条支流:泵泄露的损失功率,其泄漏液阻为R1、泵至主阀管路油液压缩的损失功率,其管道液容为C1、其余流向主阀节流口。流向主阀节流口输出的功率在相同流量情况下分为两条支流:主阀节流口的损失功率,其主阀节流口液阻为R2、其余流向压力补偿阀。
主阀节流口-压力补偿阀模块:根据功率键合图基本原理及规则,分析高负载联主阀口至压力补偿阀的动态过程中的功率流程:流向压力补偿阀输出的功率由于压力补偿阀开启需进行液压能与机械能转换,输出的功率可以分为四个支流:压力补偿阀阀芯所受液动力为F1、压力补偿阀阀芯所受惯性力,阀芯等效质量为I、压力补偿阀阀芯所受摩擦力为F2、其余机械能再次转化成液压能。压力补偿阀阀芯开启后,输入的功率流可以分为三个支流:压力补偿阀入口管路油液压缩的损失功率,管道液容为C2、流向压力补偿阀阻尼孔、流向压力补偿阀阀口。
压力补偿阀-LS管路模块:根据功率键合图基本原理及规则,分析高负载联压力补偿阀至LS管路的动态过程中的功率流程;流经压力补偿阀入口输出的一部分功率可以分为两个支流:阻尼孔损失的功率,阻尼孔液阻为R3、其余流向LS管路;流向LS管路的功率可以分为两个支流:阻尼孔至LS口管路油液压缩的损失功率,管路液容为C3、其余流向LS口;流向LS口的功率可以分为两个支流:LS管路的损失功率,LS管路液阻为R4、LS口压力Se1。
压力补偿阀-主阀换向模块:根据功率键合图基本原理及规则,分析高负载联压力补偿阀至主阀换向部分的动态过程中的功率流程;流经压力补偿阀入口输出的一部分功率可以分为两个支流:压力补偿阀阀口的损失功率,压力补偿阀阀口液阻为r5、其余流向主阀换向部分管路;流向主阀换向部分管路的功率可以分为两个支流:压力补偿阀阀口至主阀换向管路油液压缩损失的功率,管路液容为C4、其余流向主阀阀口。
主阀换向-负载溢流阀模块:根据功率键合图基本原理及规则,分析高负载联主阀至负载溢流阀的动态过程中的功率流程;流向主阀阀口的功率可以分为两个支流:主阀阀口的损失功率,主阀口液阻为R6、其余流向负载管路;流向负载管路的功率可以分为两个支流:主阀阀口至负载管路油液压缩的损失功率,管路液容为C5、其余流向负载;流向负载的功率可以分为两个支流:主阀阀口至负载管路的损失功率,管路液阻为R7、负载压力Se2。
如图4所示为本实施例建立的阀后压力补偿系统功率键合图;如表1所示为该阀后压力补偿系统的功率键合图模型中各个参数的实际值。
表1阀后压力补偿系统的功率键合图模型;
S2:推导阀后压力补偿系统的状态方程;
根据步骤S1建立的阀后压力补偿系统功率键合图模型,推导系统状态方程,则阀后压力补偿系统状态方程如下式所示:
式中:为液压泵至主阀节流口管路压缩的油液体积的一阶导数;为压力补偿阀阀芯动量的一阶导数;为压力补偿阀入口管路压缩的油液体积的一阶导数;为压力补偿阀出口至LS管路受到压缩的油液体积的一阶导数;为压力补偿阀出口至主阀换向管路压缩的油液体积的一阶导数;为主阀出口至负载溢流阀管路压缩的油液体积的一阶导数;Sf为液压泵恒定输出流量;Se1为PLS压力;Se2为负载溢流阀调定压力;r1为泵的泄漏液阻;R2为主阀芯节流口液阻;R3为压力补偿阀阻尼孔液阻;R4为LS管路液阻;r5为压力补偿阀阀口液阻;r6为主阀换向阀口液阻;r7为主阀出口至负载管路液阻;c1为泵至主阀节流口管路液容;c2为压力补偿阀入口管路液容;C3为压力补偿阀出口至LS管路液容;C4为压力补偿阀出口至主阀换向部分管路液容;C5为主阀出口至负载溢流阀管路液容;F1为压力补偿阀阀芯所受液动力;F2为压力补偿阀阀芯所受摩擦力;I1为压力补偿阀阀芯等效质量;A为压力补偿阀阀芯承压面积;A1为主阀节流口过流面积;A2为压力补偿阀阀口过流面积;A3为压力补偿阀阻尼孔过流面积;A4为主阀换向部分阀口过流面积;P10为压力补偿阀阀芯动量;V2为液压泵至主阀节流口管路压缩的油液体积;V13为压力补偿阀入口管路压缩的油液体积;V17为压力补偿阀出口至LS管路受到压缩的油液体积;V24为压力补偿阀出口至主阀换向管路压缩的油液体积;V28为主阀出口至负载溢流阀管路压缩的油液体积;Cd为流量系数;ρ为油液密度;
阀后压力补偿系统状态方程的具体推倒过程如下所示:
压力补偿系统的功率键合图中包括容性元和感性元;取容性元C的流变量q和感性元I的力变量p作为状态变量;为了便于建立状态方程,取容性元C和感性元I功率键上自变量对时间的积分为状态变量,即取液体体积V,固体或液体的动量P为状态变量;
液压系统中,对于容性元C,其力变量p与流变量q之间的关系如下式所示:
式中:p是容性元的力变量;q是容性元的流变量;dt表示时间t是q的积分变量;V为液体的油液体积;C表示容性元;
机械系统中,对于感性元I,其流变量与力变量之间的关系如下式所示:
式中:v表示感性元的流变量;F表示感性元的力变量;I表示感性元;
则根据阀后压力补偿系统功率键合图容性元件上力变量p与状态变量V、感性元件上流变量v与状态变量P的关系如下式所示:
式中:p2为功率键合图中第2个功率键上的力变量;v10为功率键合图中第10个功率键上的流变量;p13为功率键合图中第13个功率键上的力变量;p17为功率键合图中第17个功率键上的力变量;p24为功率键合图中第24个功率键上的力变量;p28为功率键合图中第28个功率键上的力变量;
阀后压力补偿系统功率键合图阻性元件上力变量与流变量的关系如下:
p3=R1q3
式中:p3为阻性元件R1上的力变量,q3为阻性元件R1上的流变量;
阀后压力补偿系统功率键合图变换器TF上力变量与力变量、流变量与流变量的关系如下:
p6·A=p7
式中:p6为功率键合图中第6个功率键上的力变量;p7为功率键合图中第7个功率键上的力变量;q6为功率键合图中第6个功率键上的流变量;q7为功率键合图中第7个功率键上的流变量;
阀后压力补偿系统功率键合图共势结,即0结,上力变量与力变量、流变量与流变量的关系如下:
p1=p2=p3=p4
q1=q2+q3+q4
式中:p1为功率键合图中第1个功率键上的力变量;p2为功率键合图中第2个功率键上的力变量;p3为功率键合图中第3个功率键上的力变量;p4为功率键合图中第4个功率键上的力变量;q1为功率键合图中第1个功率键上的流变量;q2为功率键合图中第2个功率键上的流变量;q3为功率键合图中第3个功率键上的流变量;q4为功率键合图中第4个功率键上的流变量;
阀后压力补偿系统功率键合图共流结,即1结,上力变量与力变量、流变量与流变量的关系如下:
q4=q5=q6
p4=p5+p6
式中:q4为功率键合图中第4个功率键上的流变量;q5为功率键合图中第5个功率键上的流变量;q6为功率键合图中第6个功率键上的流变量;p4为功率键合图中第4个功率键上的力变量;p5为功率键合图中第5个功率键上的力变量;p6为功率键合图中第6个功率键上的力变量;
根据功率键合图中各变量之间的映射关系及0-1节点的意义,即可推导得到状态方程。
S3:推导阀后压力补偿系统的参数灵敏度微分方程,确定参数灵敏度;
S31:简化系统的状态方程;
由于参数灵敏度微分方程是针对状态方程研究系统参数变化引起的状态矢量变化,获取步骤S2求得的系统的状态方程,将其简化成下式所示的形式:
状态矢量具体包括内容如下式所示:
x=[x1,x2,x3,x4,x5,x6]T=[V2,P10,V13,V17,V24,V28]T
式中:x1,x2,…x6分别表示第1、2…6个状态变量;
输入矢量具体包括内容如下式所示:
u=[u1,u2,u3]T=[Sf,Se1,Se2]T
式中:u1,u2,u3分别表示第1、2、3个输入变量;
参数矢量具体包括内容如下式所示:
α=[α1,α2,α3,α4,α5,α6,α7,α8,α9,α10,α11,α12,α13,α14]T
=[R1 -1,R4 -1,R7 -1,A,A1 -1,A2,A3,A4,I1 -1,C1 -1,C2 -1,C3 -1,C4 -1,C5 -1]T
式中:α1,α2,…α14分别表示第1、2…14个参数变量。
S32:求解状态方程,确定参数灵敏度微分方程;
状态方程的解为:
其中,已知稳态液动力和摩擦力对阀动态特性的影响程度较小,主要分析参数矢量α中其余各参数灵敏度,状态矢量x对参数α的灵敏度函数定义为:
式中:i的取值为1,2,3…,d;λi表示第i个参数矢量灵敏度;αi表示第i个参数变量;
状态矢量x是关于输入矢量u与参数矢量α的函数,在u与α相互独立的情况下,在状态方程的等式两边同时对参数α求偏导,可得:
由于u与α相互独立,则可得:
最终得到系统参数灵敏度微分方程如下式所示:
系数项矩阵如下所示:
式中:a11,…a66分别表示一阶轨迹灵敏度方程的系数项矩阵的36个元素;
系数矩阵中各元素如下所示:
自由项矩阵如下所示:
式中:b1,1,…b6,14分别表示一阶轨迹灵敏度方程的自由项矩阵的84个元素;自由项矩阵中各元素如下式所示:
S33:确定阀后压力补偿系统的参数灵敏度;
一阶轨迹参数灵敏度函数的初始条件为:
初始时刻,阀后压力补偿系统状态矢量x的初值x0=0,根据式可知:
根据步骤S32计算得到的一阶轨迹灵敏度微分方程,将一阶轨迹灵敏度函数初值代入,即可确定阀后压力补偿系统参数灵敏度。如图4a~4n为参数灵敏度分别为λ1…λ14时实施例阀后压力补偿系统参数一阶轨迹参数灵敏度的计算结果。
通过本实施例可以进一步证明本发明的有益效果如下所示:
(1)本发明基于功率键合图理论建立了多路阀阀后压力补偿系统的数学模型,解决了阀后压力补偿系统的灵敏度参数计算的问题;
(2)本发明为分析压力补偿阀中影响负载敏感比例多路阀动态特性的关键参数以及为优化负载敏感比例多路阀动态特性和压力补偿阀设计参数奠定研究基础和理论依据。
以上所述的实施例仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (6)
1.一种阀后压力补偿系统参数灵敏度分析方法,其特征在于,其包括以下步骤:
步骤1:构建阀后压力补偿系统的功率键合图模型;
根据负载敏感比例多路阀的工作原理和压力补偿阀的功能,基于功率键合图理论,搭建阀后压力补偿系统的功率键合图模型;所述的压力补偿系统包括:泵-主阀节流口模块、主阀节流口-压力补偿阀模块、压力补偿阀-LS管路模块、压力补偿阀-主阀换向模块和主阀换向-负载溢流阀模块五个模块;
步骤2:推导阀后压力补偿系统的状态方程;
根据步骤1建立的阀后压力补偿系统功率键合图模型,推导系统状态方程,则阀后压力补偿系统状态方程如下式所示:
式中:为液压泵至主阀节流口管路压缩的油液体积的一阶导数;为压力补偿阀阀芯动量的一阶导数;为压力补偿阀入口管路压缩的油液体积的一阶导数;为压力补偿阀出口至LS管路受到压缩的油液体积的一阶导数;为压力补偿阀出口至主阀换向管路压缩的油液体积的一阶导数;为主阀出口至负载溢流阀管路压缩的油液体积的一阶导数;Sf为液压泵恒定输出流量;Se1为PLS压力;Se2为负载溢流阀调定压力;R1为泵的泄漏液阻;R2为主阀芯节流口液阻;R3为压力补偿阀阻尼孔液阻;R4为LS管路液阻;R5为压力补偿阀阀口液阻;R6为主阀换向阀口液阻;R7为主阀出口至负载管路液阻;C1为泵至主阀节流口管路液容;C2为压力补偿阀入口管路液容;C3为压力补偿阀出口至LS管路液容;C4为压力补偿阀出口至主阀换向部分管路液容;C5为主阀出口至负载溢流阀管路液容;F1为压力补偿阀阀芯所受液动力;F2为压力补偿阀阀芯所受摩擦力;I1为压力补偿阀阀芯等效质量;A为压力补偿阀阀芯承压面积;A1为主阀节流口过流面积;A2为压力补偿阀阀口过流面积;A3为压力补偿阀阻尼孔过流面积;A4为主阀换向部分阀口过流面积;P10为压力补偿阀阀芯动量;V2为液压泵至主阀节流口管路压缩的油液体积;V13为压力补偿阀入口管路压缩的油液体积;V17为压力补偿阀出口至LS管路受到压缩的油液体积;V24为压力补偿阀出口至主阀换向管路压缩的油液体积;V28为主阀出口至负载溢流阀管路压缩的油液体积;Cd为流量系数;ρ为油液密度;
步骤3:推导阀后压力补偿系统的参数灵敏度微分方程,确定参数灵敏度;
步骤31:简化系统的状态方程;
由于参数灵敏度微分方程是针对状态方程研究系统参数变化引起的状态矢量变化,获取步骤2求得的系统的状态方程,将其简化成下式所示的形式:
步骤32:求解状态方程,确定参数灵敏度微分方程;
状态方程的解为:
其中,已知稳态液动力和摩擦力对阀动态特性的影响程度较小,主要分析参数矢量α中其余各参数灵敏度,状态矢量x对参数α的灵敏度函数定义为:
式中:i的取值为1,2,3…,d;λi表示第i个参数矢量灵敏度;αi表示第i个参数变量;
状态矢量x是关于输入矢量u与参数矢量α的函数,在u与α相互独立的情况下,在状态方程的等式两边同时对参数α求偏导,可得:
由于u与α相互独立,则可得:
最终得到系统参数灵敏度微分方程如下式所示:
步骤33:确定阀后压力补偿系统的参数灵敏度;
一阶轨迹参数灵敏度函数的初始条件为:
初始时刻,阀后压力补偿系统状态矢量x的初值x0=0,根据式可知:
根据步骤32计算得到的一阶轨迹灵敏度微分方程,将一阶轨迹灵敏度函数初值代入,即可确定阀后压力补偿系统参数灵敏度。
2.根据权利要求1所述的阀后压力补偿系统参数灵敏度分析方法,其特征在于,所述步骤1中的阀后压力补偿系统具体包括以下五个模块;
所述的泵-主阀节流口模块:根据功率键合图基本原理及规则,分析高负载联泵至主阀节流口部分的动态过程中的功率流程;恒流源代表负载所需要的流量,输出的液流功率在相同压力下分为三条支流:泵泄露的损失功率,其泄漏液阻为R1、泵至主阀管路油液压缩的损失功率,其管道液容为C1、其余流向主阀节流口;流向主阀节流口输出的功率在相同流量情况下分为两条支流:主阀节流口的损失功率,其主阀节流口液阻为R2、其余流向压力补偿阀;
所述的主阀节流口-压力补偿阀模块:根据功率键合图基本原理及规则,分析高负载联主阀口至压力补偿阀的动态过程中的功率流程;流向压力补偿阀输出的功率由于压力补偿阀开启需进行液压能与机械能转换,输出的功率可以分为四个支流:压力补偿阀阀芯所受液动力为F1、压力补偿阀阀芯所受惯性力,阀芯等效质量为I、压力补偿阀阀芯所受摩擦力为F2、其余机械能再次转化成液压能;压力补偿阀阀芯开启后,输入的功率流可以分为三个支流:压力补偿阀入口管路油液压缩的损失功率,管道液容为C2、流向压力补偿阀阻尼孔、流向压力补偿阀阀口;
所述的压力补偿阀-LS管路模块:根据功率键合图基本原理及规则,分析高负载联压力补偿阀至LS管路的动态过程中的功率流程;流经压力补偿阀入口输出的一部分功率可以分为两个支流:阻尼孔损失的功率,阻尼孔液阻为R3、其余流向LS管路;流向LS管路的功率可以分为两个支流:阻尼孔至LS口管路油液压缩的损失功率,管路液容为C3、其余流向LS口;流向LS口的功率可以分为两个支流:LS管路的损失功率,LS管路液阻为R4、LS口压力Se1;
所述的压力补偿阀-主阀换向模块:根据功率键合图基本原理及规则,分析高负载联压力补偿阀至主阀换向部分的动态过程中的功率流程;流经压力补偿阀入口输出的一部分功率可以分为两个支流:压力补偿阀阀口的损失功率,压力补偿阀阀口液阻为R5、其余流向主阀换向部分管路;流向主阀换向部分管路的功率可以分为两个支流:压力补偿阀阀口至主阀换向管路油液压缩损失的功率,管路液容为C4、其余流向主阀阀口;
所述的主阀换向-负载溢流阀模块:根据功率键合图基本原理及规则,分析高负载联主阀至负载溢流阀的动态过程中的功率流程;流向主阀阀口的功率可以分为两个支流:主阀阀口的损失功率,主阀口液阻为R6、其余流向负载管路;流向负载管路的功率可以分为两个支流:主阀阀口至负载管路油液压缩的损失功率,管路液容为C5、其余流向负载;流向负载的功率可以分为两个支流:主阀阀口至负载管路的损失功率,管路液阻为R7、负载压力Se2;
所述阀后压力补偿系统的功率键合图模型是由上述各模块得到的功率键合图模型连接组成。
3.根据权利要求1所述的阀后压力补偿系统参数灵敏度分析方法,其特征在于,所述步骤2中推导系统状态方程具体过程如下:
所述的压力补偿系统包括容性元和感性元;取容性元C的流变量q和感性元I的力变量p作为状态变量;为了便于建立状态方程,取容性元C和感性元I功率键上自变量对时间的积分为状态变量,即取液体体积V,固体或液体的动量P为状态变量;
液压系统中,对于容性元C,其力变量p与流变量q之间的关系如下式所示:
式中:p是容性元的力变量;q是容性元的流变量;dt表示时间t是q的积分变量;V为液体的油液体积;C表示容性元;
机械系统中,对于感性元I,其流变量与力变量之间的关系如下式所示:
式中:v表示感性元的流变量;F表示感性元的力变量;I表示感性元;
则根据阀后压力补偿系统功率键合图容性元件上力变量p与状态变量V、感性元件上流变量v与状态变量P的关系如下式所示:
式中:p2为功率键合图中第2个功率键上的力变量;v10为功率键合图中第10个功率键上的流变量;p13为功率键合图中第13个功率键上的力变量;p17为功率键合图中第17个功率键上的力变量;p24为功率键合图中第24个功率键上的力变量;p28为功率键合图中第28个功率键上的力变量;
阀后压力补偿系统功率键合图阻性元件上力变量与流变量的关系如下:
p3=R1q3
式中:p3为阻性元件R1上的力变量,q3为阻性元件R1上的流变量;
阀后压力补偿系统功率键合图变换器TF上力变量与力变量、流变量与流变量的关系如下:
p6·A=p7
式中:p6为功率键合图中第6个功率键上的力变量;p7为功率键合图中第7个功率键上的力变量;q6为功率键合图中第6个功率键上的流变量;q7为功率键合图中第7个功率键上的流变量;
阀后压力补偿系统功率键合图共势结,即0结,上力变量与力变量、流变量与流变量的关系如下:
p1=p2=p3=p4
q1=q2+q3+q4
式中:p1为功率键合图中第1个功率键上的力变量;p2为功率键合图中第2个功率键上的力变量;p3为功率键合图中第3个功率键上的力变量;p4为功率键合图中第4个功率键上的力变量;q1为功率键合图中第1个功率键上的流变量;q2为功率键合图中第2个功率键上的流变量;q3为功率键合图中第3个功率键上的流变量;q4为功率键合图中第4个功率键上的流变量;
阀后压力补偿系统功率键合图共流结,即1结,上力变量与力变量、流变量与流变量的关系如下:
q4=q5=q6
p4=p5+p6
式中:q4为功率键合图中第4个功率键上的流变量;q5为功率键合图中第5个功率键上的流变量;q6为功率键合图中第6个功率键上的流变量;p4为功率键合图中第4个功率键上的力变量;p5为功率键合图中第5个功率键上的力变量;p6为功率键合图中第6个功率键上的力变量;
根据功率键合图中各变量之间的映射关系及0-1节点的意义,即可推导得到状态方程。
4.根据权利要求1所述的阀后压力补偿系统参数灵敏度分析方法,其特征在于,所述步骤31中的状态矢量、输入矢量、参数矢量具体如下所示:
所述的状态矢量具体包括内容如下式所示:
x=[x1,x2,x3,x4,x5,x6]T=[V2,P10,V13,V17,V24,V28]T
式中:x1,x2,…x6分别表示第1、2…6个状态变量;
所述的输入矢量具体包括内容如下式所示:
u=[u1,u2,u3]T=[Sf,Se1,Se2]T
式中:u1,u2,u3分别表示第1、2、3个输入变量;
所述的参数矢量具体包括内容如下式所示:
α=[α1,α2,α3,α4,α5,α6,α7,α8,α9,α10,α11,α12,α13,α14]T
=[R1 -1,R4 -1,R7 -1,A,A1 -1,A2,A3,A4,I1 -1,C1 -1,C2 -1,C3 -1,C4 -1,C5 -1]T
式中:α1,α2,…α14分别表示第1、2…14个参数变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210015137.9A CN114396416B (zh) | 2022-01-07 | 2022-01-07 | 阀后压力补偿系统参数灵敏度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210015137.9A CN114396416B (zh) | 2022-01-07 | 2022-01-07 | 阀后压力补偿系统参数灵敏度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114396416A CN114396416A (zh) | 2022-04-26 |
CN114396416B true CN114396416B (zh) | 2022-10-25 |
Family
ID=81228462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210015137.9A Active CN114396416B (zh) | 2022-01-07 | 2022-01-07 | 阀后压力补偿系统参数灵敏度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114396416B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105443471A (zh) * | 2015-12-04 | 2016-03-30 | 湖南三一快而居住宅工业有限公司 | 一种多路阀和多路阀的流量补偿控制系统及方法 |
CN113074155A (zh) * | 2021-04-02 | 2021-07-06 | 同济大学 | 一种负压开启插装式补油阀 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6334308B1 (en) * | 1998-03-04 | 2002-01-01 | Komatsu Ltd. | Pressure compensating valve, unloading pressure control valve and hydraulically operated device |
DE102007055377A1 (de) * | 2007-11-19 | 2009-05-20 | Robert Bosch Gmbh | Hydraulische Steueranordnung und Wegeventilsektion |
WO2009100089A1 (en) * | 2008-02-04 | 2009-08-13 | Hendrickson Usa, L.L.C. | Multi-stage height control valve including position sensitive pilot signal and pressure boost for vehicle air springs |
FR2981414B1 (fr) * | 2011-10-12 | 2013-10-18 | In Lhc | Servovalve a deux etages et etage de pilotage adapte a une telle servovalve. |
CN102799705B (zh) * | 2012-06-02 | 2015-07-22 | 福州大学 | 负载敏感控制轴向柱塞泵虚拟样机的建立方法 |
CN103573731B (zh) * | 2012-07-24 | 2015-12-16 | 徐工集团工程机械股份有限公司 | 流量控制装置以及液压油散热系统 |
CN105570491B (zh) * | 2014-10-17 | 2018-03-06 | 徐工集团工程机械股份有限公司 | 负载敏感多路阀、负载敏感系统及挖掘机 |
DE102017210823A1 (de) * | 2017-06-27 | 2018-12-27 | Robert Bosch Gmbh | Ventilblockanordnung und Verfahren für eine Ventilblockanordnung |
DE102018209793A1 (de) * | 2018-06-18 | 2019-12-19 | Zf Friedrichshafen Ag | Verfahren und Steuerungseinrichtung zur Ermittlung von Leckage an einem pneumatischen Stellsystem |
IT201800009591A1 (it) * | 2018-10-18 | 2020-04-18 | Walvoil Spa | Sistema idraulico di tipo load sensing con dispositivo idraulico di regolazione |
CN109441905B (zh) * | 2018-12-26 | 2020-01-07 | 太原理工大学 | 一种变压差负载敏感多路阀 |
CN109654009B (zh) * | 2019-02-21 | 2020-04-17 | 燕山大学 | 一种液压泵马达性能测试系统及其测试方法 |
CN110953205B (zh) * | 2019-12-16 | 2022-03-22 | 厦门理工学院 | 一种陶瓷压砖机电液比例阀位置控制方法、装置、设备及系统 |
CN214698546U (zh) * | 2021-02-22 | 2021-11-12 | 四川工程职业技术学院 | 一种负载敏感液压系统 |
-
2022
- 2022-01-07 CN CN202210015137.9A patent/CN114396416B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105443471A (zh) * | 2015-12-04 | 2016-03-30 | 湖南三一快而居住宅工业有限公司 | 一种多路阀和多路阀的流量补偿控制系统及方法 |
CN113074155A (zh) * | 2021-04-02 | 2021-07-06 | 同济大学 | 一种负压开启插装式补油阀 |
Also Published As
Publication number | Publication date |
---|---|
CN114396416A (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ramos et al. | Surge damping analysis in pipe systems: Modelling and experiments | |
Shabayek et al. | Dynamic model for subcritical combining flows in channel junctions | |
CN112415891B (zh) | 一种电液伺服系统自适应输出反馈渐近控制方法 | |
CN105843043B (zh) | 一种电液负载模拟装置自适应鲁棒力控制方法 | |
CN104698844A (zh) | 液压位置伺服系统的不确定性补偿的滑模控制方法 | |
CN114396416B (zh) | 阀后压力补偿系统参数灵敏度分析方法 | |
Dihovicni et al. | Mathematical modelling and simulation of pneumatic systems | |
Ramalingam et al. | Using the WCM for transient modeling of water distribution networks | |
Jaeger et al. | An extended κ–ϵ finite element model | |
Dulhoste et al. | Non-linear control of water flow dynamics by input-output linearization based on a collocation method model | |
Lin et al. | A novel control strategy for an energy saving hydraulic system with near-zero overflowing energy-loss | |
Boulos et al. | An explicit algorithm for calculating operating parameters for water networks | |
Guinot | Boundary condition treatment in 2× 2 systems of propagation equations | |
CN210178681U (zh) | 一种阻尼可调的液压马达的执行器 | |
Pfeiffer et al. | New hydraulic system modelling | |
Fu et al. | Research on optimal control of excavator negative control system based on secondary controllable main valve | |
Burt et al. | Irrigation-canal-simulation model usage | |
JPH11351045A (ja) | エンジンの状態を示す量の推定方法 | |
Abdalla et al. | A computational study of the actuation speed of the hydraulic cylinder under different ports’ sizes and configurations | |
Najm et al. | Hydraulic ram analysis: a new look at an old problem | |
Yang et al. | Pressure control of hydraulic servo system using proportional control valve | |
Sciancalepore et al. | Lumped parameter modeling of counterbalance valves considering the effect of flow forces | |
Marmolejo Gutierrez | Stability of tank level control systems in hydropower plants | |
Badillo-Olvera et al. | Leak isolation in pressurized pipelines using an interpolation function to approximate the fitting losses | |
James et al. | An approach to modeling real-time control of dynamic and static radial and sluice gates within EXTRAN |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |