CN105574329A - 基于周期分段函数移相叠加算法的脉动概率密度获取方法 - Google Patents

基于周期分段函数移相叠加算法的脉动概率密度获取方法 Download PDF

Info

Publication number
CN105574329A
CN105574329A CN201510924468.4A CN201510924468A CN105574329A CN 105574329 A CN105574329 A CN 105574329A CN 201510924468 A CN201510924468 A CN 201510924468A CN 105574329 A CN105574329 A CN 105574329A
Authority
CN
China
Prior art keywords
function
pulsation
waypoint
phase shift
probability density
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.)
Pending
Application number
CN201510924468.4A
Other languages
English (en)
Inventor
李晶
廖攀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN201510924468.4A priority Critical patent/CN105574329A/zh
Publication of CN105574329A publication Critical patent/CN105574329A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Complex Calculations (AREA)

Abstract

本发明涉及一种基于周期分段函数移相叠加算法的脉动概率密度获取方法,用于实现双泵合流流量脉动概率密度函数的获取,该方法具体包括以下步骤:1)分别获取两个泵的单泵流量函数F(x)和 为双泵间的移相角,其中,F为周期分段函数,周期为T,2)采用周期分段函数移相叠加算法计算以获得的表达式作为双泵合流流量函数3)根据步骤2)获得双泵合流流量函数计算双泵合流流量脉动概率密度函数。与现有技术相比,本发明具有计算效率高、简单方便等优点。

Description

基于周期分段函数移相叠加算法的脉动概率密度获取方法
技术领域
本发明涉及双泵合流流量特性获取方法,尤其是涉及一种基于周期分段函数移相叠加算法的脉动概率密度获取方法。
背景技术
液压系统常采用双泵合流来提高能量利用率。双泵合流的流量为两个泵流量之和。单泵流量为关于泵转角的周期函数,且多表现为周期分段函数。双泵合流流量即为两个周期分段函数叠加。由于双泵初始转角相位不同,两个周期分段函数的初始相位不同,叠加的函数表达式也不同。在分段较多的情况下,叠加过程非常复杂,直接进行叠加计算费时费力且容易出错,因此需要设计一种通用的周期分段函数移相叠加方法,以方便实现双泵合流流量特性的获取。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种计算效率高、简单方便的基于周期分段函数移相叠加算法的脉动概率密度获取方法。
本发明的目的可以通过以下技术方案来实现:
一种基于周期分段函数移相叠加算法的脉动概率密度获取方法,用于实现双泵合流流量脉动概率密度函数的获取,该方法具体包括以下步骤:
1)分别获取两个泵的单泵流量函数F(x)和为双泵间的移相角,其中,F为周期分段函数,周期为T,
2)采用周期分段函数移相叠加算法计算以获得的表达式作为双泵合流流量函数
3)根据步骤2)获得双泵合流流量函数计算双泵合流流量脉动概率密度函数。
所述周期分段函数F的分段数为n,每段长度为li在区间[0,T]中,每个分段对应的函数为fi(x),i=1,2,3…,n。
所述周期分段函数移相叠加算法具体为:
201)根据周期分段函数中的li写出数列A和B:
A = Σ u = k i l u , k = 1 , 2 , ... , i ( i = 1 , 2 , ... , n )
B = Σ r = 1 i l r + Σ s = k n l s , k = i + 2 , i + 3 , ... , n ( i = 1 , 2 , ... , n - 1 )
202)将数列A和B中的数从小到大排序,删除重复的数,得到数列:
k1,k2,k3…,km
其中km=T,m表示排序后的数列元素个数,m≤n2-n+1,利用该数列将的分段区间设置为[0,k1),[k1,k2),[k2,k3),…,[km-1,km);
203)对的任一分段区间[kr,kr+1),1≤r≤m-1,任意取该分段区间中的某个值,计算移相叠加函数分段点位置,找出范围在(0,T)之间的分段点,并由小到大排序得到数列:
s1,s2,s3…,st
其中st=T,t表示找出的分段点个数,设置叠加函数自变量x的分段区间为[0,s1),[s1,s2),[s2,s3),…,[st-1,st);
204)对于叠加函数自变量x的每一分段区间,根据分段点所在位置对应的函数计算叠加函数各个分段函数,获得移相叠加函数的表达式。
所述步骤203)中,移相叠加函数分段点位置由以下3个部分的集合中获取:
第一个周期分段点
第二个周期分段点
F(x)分段点0,l1,l1+l2,l1+l2+l3,…,l1+l2+l3+…+lm
所述步骤204)中,根据分段点所在位置对应的函数计算叠加函数各个分段函数具体为:
401)设定两个索引值p、q,其中,令p的初始值为0,q的初始值通过以下方式确定:
若s1源自第一个周期分段点,则令q=m+n,若s1源自第二个周期分段点,则令q=m+2n;
402)计算移相叠加函数在区间[0,s1)上的函数表达式:
其中,j=p或q;
403)遍历si,i=1,2…,t-1,若si源自F(x)分段点,则令p=p+1,若源自于分段点,则令q=q+1,其对应区间x∈[si,si+1)上的函数表达式为:
404)改变r的值,求出的不同分段区间上,的表达式。
与现有技术相比,本发明具有以下优点:
1)本发明设计了一种新型周期分段函数移相叠加算法,实现了两个不同初始相角的周期函数的自动叠加计算,计算效率高,计算结果准确。
2)本发明的周期分段函数移相叠加算法中,函数分段数n和每段长度l均为可变量,亦即在对应双泵合流的流量脉动概率密度函数的计算中,涉及单泵结构与流量概率密度函数对应的柱塞数和柱塞孔间距及缓冲槽等的数值均可变化,在实际应用中非常方便。
3)本发明具有可拓展性,可以进一步变形,拓展至三个及以上周期函数的叠加,实现三泵或者更多泵组合的流量密度函数的计算。
附图说明
图1为本发明的流程示意图;
图2为实施例中的分段周期函数F(x)示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本实施例提供一种基于周期分段函数移相叠加算法的脉动概率密度获取方法,用于实现双泵合流流量脉动概率密度函数的获取,该方法具体包括以下步骤:
1)分别获取两个泵的单泵流量函数F(x)和为双泵间的移相角,其中,F为周期分段函数,周期为T,周期分段函数F的分段数为n,每段长度为li在区间[0,T]中,每个分段对应的函数为fi(x),i=1,2,3…,n;
2)采用周期分段函数移相叠加算法计算以获得的表达式作为双泵合流流量函数
3)根据步骤2)获得双泵合流流量函数计算双泵合流流量脉动概率密度函数。
如图1所示,周期分段函数移相叠加算法具体为:
201)根据周期分段函数中的li写出数列A和B(仅在n≥3时才存在B数列):
A = Σ u = k i l u , k = 1 , 2 , ... , i ( i = 1 , 2 , ... , n )
B = Σ r = 1 i l r + Σ s = k n l s , k = i + 2 , i + 3 , ... , n ( i = 1 , 2 , ... , n - 1 ) .
202)将数列A和B中的数从小到大排序,得到数列:
k1,k2,k3…,km
其中km=T,m的值与分段数n以及每个分段的长度li有关,如果A,B数列中没有重复的数,那么此时m最大,为n2-n+1。由于A,B数列中可能存在重复的数,从小到大排列后去掉重复的数,因此有m≤n2-n+1。利用该数列将的分段区间设置为[0,k1),[k1,k2),[k2,k3),0,[km-1,km)。
203)对的任一分段区间[kr,kr+1),1≤r≤m-1,任意取该分段区间中的某个值,计算移相叠加函数分段点位置,找出范围在(0,T)之间的分段点,并由小到大排序得到数列:
s1,s2,s3…,st
其中st=T,t表示找出的分段点个数,考虑到可能有重复元素,t的值在不同情况下有所不同,具体数值与F(x)分段情况有关,分段点详细计算在下面步骤中会给出。设置叠加函数自变量x的分段区间为[0,s1),[s1,s2),[s2,s3),…,[st-1,st)。
204)对于叠加函数自变量x的每一分段区间,只要知道分段点在原函数上还是在移相函数上,然后根据分段点所在位置对应的函数(分段点前面或后面的函数),就可计算叠加函数各个分段函数,从而获得移相叠加函数的表达式。
根据分段点所在位置对应的函数计算叠加函数各个分段函数具体为:
401)设定两个索引值p、q,其中,令p的初始值为0,q的初始值通过以下方式确定:
若s1源自第一个周期分段点,则令q=m+n,若s1源自第二个周期分段点,则令q=m+2n。
402)计算移相叠加函数在区间[0,s1)上的函数表达式:
其中,j=p或q;
j的取值范围为1~3n,其原因在于步骤203)涉及的分段点包含3个部分,而每个部分的分段点区间对应函数不同,对于F(x)分段点,相应区间上的函数用fi(x)来表示;对于第一个周期分段点,相应区间上的函数用表示;对于第二个周期分段点,相应区间上的函数用表示。均可以用fi函数表示,为了方便编程计算,人为地将这三种函数从1开始编号,每个编号对应一个函数,由于i=1,2,3…,n,即F(x)分段点区间上的函数有n个,3个部分就是索引值就为1~3n,每个编号就对应一个函数,方便计算机运算表达。
403)遍历si,i=1,2…,t-1,若si源自F(x)分段点,则令p=p+1,若源自于分段点,则令q=q+1,其对应区间x∈[si,si+1)上的函数表达式为:
据此可以求出当时,在区间[0,T)上的全部表达式。
取不同的r值,即可求出的不同分段区间上,的表达式,从而求出全部的表达式,以获得的表达式作为双泵合流流量函数可计算双泵合流流量脉动概率密度函数。
计算双泵合流流量脉动概率密度函数具体过程为:
a、首先根据计算流量脉动函数
其中分别表示双泵合流流量函数最大值和最小值,表示双泵合流流量函数的平均值。
b、根据流量脉动函数计算流量脉动分布函数F(δ)及其导数f(δ),f(δ)表示流量脉动为δ时的概率密度,即为双泵合流流量脉动特性,其值表达了出现流量脉动为δ的可能性大小。其中,流量脉动分布函数F(δ)为
F ( &delta; ) = 0 , &delta; < &delta; min L ( &delta; ) / T , &delta; min &le; &delta; &le; &delta; max 1 , &delta; > &delta; max
其中δmin和δmax分别表示流量脉动函数的最小值和最大值,L(δ)表示满足所对应的区间长度。
如果双泵合流流量脉动函数在区间单调递减;或在某个区间内单调递减,且该区间内计算的分布函数F(δ)与在区间上计算结果一致,那么所对应的区间长度L(δ)为
L(δ)=T-g-1(δ)
其中g-1(δ)为流量脉动函数的反函数。
双泵合流流量脉动函数在区间单调递减,易得满足取值范围为那么有区间的长度即为L(δ)=T-g-1(δ)。
以7柱塞轴向柱塞泵为例,其单泵的流量函数为
该函数周期为2β,其中:
K为与柱塞泵自身结构及转速有关的系数;
β为相邻两柱塞夹角的一半,即β=π/Z,Z为柱塞数目,在这里等于7;
为柱塞泵的转角;
——闭死角,这里为π/12;
——错配角,这里为π/24;
显然,该周期分段函数分为3段,根据发明,采用计算机编程计算,得到的周期分段函数移相叠加函数表达式为
可以看出,原函数仅仅分了三段,但计算出来的结果很复杂,而采用本发明方法可以很方便地计算周期分段函数移相叠加函数。
以上所述,仅为本发明的较佳实施实例,不能以之限制本发明权利要求的范围。即凡依本发明权利要求所做的均等变化及修饰,仍将不失发明的要义所在,亦不脱离本发明的精神和权利要求范围,故都应视为本发明的进一步实施状况。

Claims (5)

1.一种基于周期分段函数移相叠加算法的脉动概率密度获取方法,用于实现双泵合流流量脉动概率密度函数的获取,其特征在于,该方法具体包括以下步骤:
1)分别获取两个泵的单泵流量函数F(x)和 为双泵间的移相角,其中,F为周期分段函数,周期为T,
2)采用周期分段函数移相叠加算法计算以获得的表达式作为双泵合流流量函数
3)根据步骤2)获得双泵合流流量函数计算双泵合流流量脉动概率密度函数。
2.根据权利要求1所述的基于周期分段函数移相叠加算法的脉动概率密度获取方法,其特征在于,所述周期分段函数F的分段数为n,每段长度为li在区间[0,T]中,每个分段对应的函数为fi(x),i=1,2,3...,n。
3.根据权利要求2所述的基于周期分段函数移相叠加算法的脉动概率密度获取方法,其特征在于,所述周期分段函数移相叠加算法具体为:
201)根据周期分段函数中的li写出数列A和B:
A = &Sigma; u = k i l u , k = 1 , 2 , ... , i ( i = 1 , 2 , ... , n )
B = &Sigma; r = 1 i l r + &Sigma; s = k n l s , k = i + 2 , i + 3 , ... , n ( i = 1 , 2 , ... , n - 1 )
202)将数列A和B中的数从小到大排序,删除重复的数,得到数列:
k1,k2,k3...,km
其中km=T,m表示排序后的数列元素个数,m≤n2-n+1,利用该数列将的分段区间设置为[0,k1),[k1,k2),[k2,k3),...,[km-1,km);
203)对的任一分段区间[kr,r+1),1≤r≤m-1,任意取该分段区间中的某个值,计算移相叠加函数分段点位置,找出范围在(0,T)之间的分段点,并由小到大排序得到数列:
s1,s2,s3...,st
其中st=T,t表示找出的分段点个数,设置叠加函数自变量x的分段区间为[0,s1),[s1,2),[s2,s3),...,[st-1,st);
204)对于叠加函数自变量x的每一分段区间,根据分段点所在位置对应的函数计算叠加函数各个分段函数,获得移相叠加函数的表达式。
4.根据权利要求3所述的基于周期分段函数移相叠加算法的脉动概率密度获取方法,其特征在于,所述步骤203)中,移相叠加函数分段点位置由以下3个部分的集合中获取:
第一个周期分段点
第二个周期分段点
F(x)分段点0,l1,l1+l2,l1+l2+l3,...,l1+l2+l3+…+ln
5.根据权利要求4所述的基于周期分段函数移相叠加算法的脉动概率密度获取方法,其特征在于,所述步骤204)中,根据分段点所在位置对应的函数计算叠加函数各个分段函数具体为:
401)设定两个索引值p、q,其中,令p的初始值为0,q的初始值通过以下方式确定:
若s1源自第一个周期分段点,则令q=m+n,若s1源自第二个周期分段点,则令q=m+2n;
402)计算移相叠加函数在区间[0,s1)上的函数表达式:
其中,j=p或q;
403)遍历si,i=1,2...,t-1,若si源自F(x)分段点,则令p=p+1,若源自于分段点,则令q=q+1,其对应区间x∈[si,si+1)上的函数表达式为:
404)改变r的值,求出的不同分段区间上,的表达式。
CN201510924468.4A 2015-12-11 2015-12-11 基于周期分段函数移相叠加算法的脉动概率密度获取方法 Pending CN105574329A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510924468.4A CN105574329A (zh) 2015-12-11 2015-12-11 基于周期分段函数移相叠加算法的脉动概率密度获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510924468.4A CN105574329A (zh) 2015-12-11 2015-12-11 基于周期分段函数移相叠加算法的脉动概率密度获取方法

Publications (1)

Publication Number Publication Date
CN105574329A true CN105574329A (zh) 2016-05-11

Family

ID=55884455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510924468.4A Pending CN105574329A (zh) 2015-12-11 2015-12-11 基于周期分段函数移相叠加算法的脉动概率密度获取方法

Country Status (1)

Country Link
CN (1) CN105574329A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109900422A (zh) * 2019-01-02 2019-06-18 同济大学 一种基于图形的多泵合流流量脉动特性检测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3506068B2 (ja) * 1999-09-29 2004-03-15 日本電気株式会社 外れ値度計算装置
CN104005924A (zh) * 2013-02-25 2014-08-27 白巨章 大功率径向柱塞泵

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3506068B2 (ja) * 1999-09-29 2004-03-15 日本電気株式会社 外れ値度計算装置
CN104005924A (zh) * 2013-02-25 2014-08-27 白巨章 大功率径向柱塞泵

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JING LI ET AL.: "Analysis of the flow ripple stochastic characteristics in dual-pump converging system", 《2015 INTERNATIONAL CONFERENCE ON FLUID POWER AND MECHATRONICS》 *
梁晓东: "液压系统双泵合流技术探究", 《有色设备》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109900422A (zh) * 2019-01-02 2019-06-18 同济大学 一种基于图形的多泵合流流量脉动特性检测方法

Similar Documents

Publication Publication Date Title
CN103942425B (zh) 一种数据处理方法和装置
CN103276686A (zh) 一种梯级水库下游设计洪水的推求方法
CN104408241B (zh) 一种修形圆柱齿轮的有限元网格自动生成方法
CN102141064A (zh) 空间过滤法建立湍流模型的构建方法
CN102982250A (zh) 基于随机响应面估计参数不确定性的随机模型修正方法
CN106875047A (zh) 水库流域径流预报方法及系统
CN104182634B (zh) 一种梯级水电站联合运行水位控制断面优选方法
CN105574329A (zh) 基于周期分段函数移相叠加算法的脉动概率密度获取方法
CN106320256A (zh) 一种多沙河流抽水蓄能电站过机泥沙通量的确定方法
CN103177166B (zh) 一种基于多项式拟合的模锻件毛坯设计方法
CN101748703A (zh) 洪水再分配方法
CN106320255A (zh) 一种漫滩水流滩槽过流量的计算方法
CN104675713A (zh) 一种基于数据样本的离心泵无过载设计方法
CN103970610A (zh) 一种供水管网中节点流量的监控方法
CN104021288A (zh) 用于导管架平台频谱疲劳分析的基本波确定方法
CN105302520A (zh) 一种倒数运算的求解方法及系统
CN105512352A (zh) 一种面向气动载荷的分析方法
Burova et al. On construction third order approximation using values of integrals
CN106066912A (zh) 一种多分区结构化网格的生成方法
Huynh An upwind moment scheme for conservation laws
CN104296748B (zh) 基于双边补偿结构的捷联惯导姿态圆锥补偿算法设计方法
CN114239434A (zh) 一种针对树网格格子玻尔兹曼方法中虚拟分层边界的高精度处理方法
CN113268710A (zh) 一种干旱极值信息输出方法、装置、电子设备和介质
CN105550510A (zh) 一种双泵合流流量脉动特性提取方法
CN113421341A (zh) 一种水库冲刷淤积计算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160511

WD01 Invention patent application deemed withdrawn after publication