一种泵体阀杆的减振设计方法

1.本发明属于机械结构减振设计技术领域，涉及了一种泵体阀杆的减振设计方法。


背景技术：

2.泵是机械工业的主要产品也是生产生活的一种通用机械。以我国东北输油管网系统为例，该系统利用太阳生泵站和管道系统对石油能源进行输送。随着设备服役年份增加以及改造工程的影响，这一能源输送系统出现了一系列影响生产安全和工作性能的因素。振动就是其中之一，而泵体阀杆的振动是最激烈的，这是由于阀杆顶端的控制电机和阀杆立柱构成的振动体把阀盖体上盖传来的振动放大造成的。阀杆振动会使得阀出口单向阀螺栓松动、漏油甚至造成管道破裂等严重后果。因此要保证能源运输的安全和效益必须对振动进行减振。而不同的工作状态受到的激振频率也不尽相同，对阀杆进行减振时需要考虑到频率多样的影响因素，针对特定频率的减振效果变得尤为重要。声子晶体是一种人工复合材料，对特定频率的弹性波具有优良的抑制作用。在泵站处于工作状态时，某一频率的弹性波传播到阀杆时使得阀杆产生共振等放大振动的现象。实际生活中，振动无法避免，将振动的影响危害降低到一定程度是大多数工程应用的目标。可应用声子晶体结构的减振原理对阀杆进行设计，使得阀杆在工作过程中受到某一频率振动时不是放大振动而是能起到抑制振动的关键作用，并对减振的程度达到满意的效果即不影响泵的工作效果。因此有必要对一种泵体阀杆的减振设计提供一种方法。


技术实现要素：

3.发明目的：本发明的目的是为了解决现有技术中存在的缺陷，而提出了一种泵体阀杆的减振设计方法。
4.技术方案：为了实现上述目的，本发明采用了如下技术方案：
5.本发明所述的一种泵体阀杆的减振设计方法，其具体操作步骤如下：
6.步骤(1)、对阀杆主体以等截面二组元声子晶体杆为目标进行离散处理，将阀杆主体离散成弹簧-质点系统；
7.步骤(2)、根据建立好的阀杆主体的弹簧-质点系统，计算出弹簧-质点系统的振动方程；
8.步骤(3)、将计算出的弹簧-质点系统的振动方程与bloch定理结合，应用弹性波在周期性结构中传播的原理建立周期性阀杆减振频率及衰减因子计算的色散多项式；
9.步骤(4)、明确色散多项式和减振频率以及色散多项式和衰减因子的对应关系；
10.步骤(5)、建立满足指定减振频率及衰减因子的逆向设计存在性判据；
11.步骤(6)、利用sturm定理结合判据求解符合目标减振频率及衰减因子要求的阀杆体参数。
12.进一步的，所述步骤(2)中，在计算弹簧-质点系统的振动方程中，需对阀杆主体部分的刚度进行求解，即求解弹簧的刚度；
13.其中，弹簧的刚度近似等于截面所受应力除以应变，
14.则同种材料之间的弹簧等效刚度z表示为：
[0015][0016]
式中，e表示该等效弹簧所对应材料的杨氏模量，r表示杆的半径，c表示离散单元轴向长度；
[0017]
而两种不同种材料之间的弹簧等效刚度z则表示为：
[0018][0019]
式中，e1表示第一种材料的杨氏模量；e2表示第二种材料的杨氏模量；
[0020]
另外，弹簧-质点系统中的第j个质点质量mj表示
[0021]
为：mj＝πρjcr2(3)式中，ρ表示离散单元密度，ρj表示第j个质点对应的离散单元的单元密度；
[0022]
其次，根据牛顿第二定律可知弹簧-质点系统中第j个质点的运动方程：
[0023]mjaj
＝zj(v
j 1-vj)-z
j-1
(v
j-v
j-1
)(j＝1,2,...,n)(4)
[0024]
式中，aj表示第j个质点的加速度，等于位移的二阶导，v表示位移。
[0025]
进一步的，在所述步骤(3)中，将弹簧-质点系统的振动方程与bloch定理结合如下式所示：
[0026][0027]
式中，u表示振幅，i表示虚数的单位，q表示波矢，l表示杆的长度，n表示离散单元数，ω表示频率，t表示时间；
[0028]
上述得到的振动方程(4)可变换得：
[0029][0030]
将振动方程(6)同周期性边界条件结合可得一线性方程组，展开为矩阵形式表示为：
[0031]
[ω2i-x(q)]u＝0(7)
[0032]
其中，i表示单位矩阵，x(q)表示关于材料参数和结构参数以及波矢和频率有关的矩阵，u表示是由各质点振幅组成的列向量；
[0033]
令矩阵x(q)-ω2i行列式为零，可得阀杆主体部分单胞的色散关系方程：
[0034]an
ω
2n
a
n-1
ω
2(n-1)

…
a1ω2 a0＝0(8)
[0035]
式中，a0～an表示与q、r、e、ρ、l、n相关的待定系数；
[0036]
其中，当色散关系方程(8)不等于0时，则定义其为色散多项式：
[0037]
q(ω2)＝anω
2n
a
n-1
ω
2(n-1)

…
a1ω2 a0(9)
[0038]
设式(7)中的波矢q由实部和虚部来表示：
[0039]
q＝g bi(10)
[0040]
式中，g和b均表示为实数，i表示虚数符号；当b不为零时则表示衰减存在，b的数值
大小表示该结构衰减程度的大小，故b被称为衰减因子。
[0041]
进一步的，在步骤(4)中，所述明确色散多项式和减振频率以及色散多项式和衰减因子的对应关系是：
[0042]
首先，将数值带入到色散多项式并进行减振频率的求解计算；
[0043]
接着，通过控制变量法分析色散多项式和指定减振频率、衰减因子的对应关系，从而明确色散多项式和减振频率以及色散多项式和衰减因子的对应关系。
[0044]
进一步的，在步骤(5)中，所述建立满足指定减振频率及衰减因子的逆向设计存在性判据如下式所示：
[0045][0046]
其中，b
max
表示是需要满足的衰减因子大小；
[0047]
为满足波矢q在取值范围内的取值要求，波矢取值后色散多项式应满足如下式要求：
[0048][0049]
有益效果：本发明与现有技术相比，本发明的特点是：本发明将声子晶体原理引入阀杆的结构设计中，通过变换阀杆构造来实现减振；可根据阀杆工作时受到的振动频率快速设计出能够对目标减振频率及减振程度进行减振的阀杆；根据本发明能对不同工作环境下的不同泵站提供可行的减振阀杆设计方法。
附图说明
[0050]
图1是本发明的操作流程图；
[0051]
图2是本发明中一种泵体阀杆的减振设计方法的阀杆模型图，
[0052]
图中：1是螺纹连接部、2是泵体阀盖体、3是杆肩、4是阀杆主体部、5是螺纹连接部、6是垫圈；
[0053]
图3为本发明一种泵体阀杆的减振设计的杆主体部模型离散示意图；
[0054]
图4是本发明一种泵体阀杆的减振设计中求解结果数据点集图；
[0055]
其中，纵坐标表示杨氏模量横坐标表示密度；
[0056]
图5是本发明中一个编程计算流程图；
[0057]
图6是本发明实施例中求解结果参数绘制得衰减因子曲线图；
[0058]
图7是本发明实施例中非求解结果参数绘制得衰减因子曲线图；
[0059]
图8是本发明一种泵体阀杆的减振设计方法设计出的元胞结构的仿真结果模态图；
[0060]
图9是本发明一种泵体阀杆的减振设计方法设计出的元胞结构的仿真结果模态图；
[0061]
图10是本发明一种泵体阀杆的减振设计方法设计出的元胞结构的仿真结果模态图；
[0062]
图11是本发明一种泵体阀杆的减振设计方法设计出的元胞结构的仿真结果模态图。
具体实施方式
[0063]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
[0064]
如图所述，本发明所述的本发明所述的一种泵体阀杆的减振设计方法，其具体操作步骤如下：
[0065]
步骤(1)、对阀杆主体以等截面二组元声子晶体杆为目标进行离散处理，将阀杆主体离散成弹簧-质点系统；
[0066]
步骤(2)、根据建立好的阀杆主体的弹簧-质点系统，计算出弹簧-质点系统的振动方程；
[0067]
步骤(3)、将计算出的弹簧-质点系统的振动方程与bloch定理结合，应用弹性波在周期性结构中传播的原理建立周期性阀杆减振频率及衰减因子计算的色散多项式；
[0068]
步骤(4)、明确色散多项式和减振频率以及色散多项式和衰减因子的对应关系；
[0069]
步骤(5)、建立满足指定减振频率及衰减因子的逆向设计存在性判据；
[0070]
步骤(6)、利用sturm定理结合判据求解符合目标减振频率及衰减因子要求的阀杆体参数。
[0071]
进一步的，所述步骤(2)中，在计算弹簧-质点系统的振动方程中，需对阀杆主体部分的刚度进行求解，即求解弹簧的刚度；
[0072]
其中，弹簧的刚度近似等于截面所受应力除以应变，
[0073]
则同种材料之间的弹簧等效刚度z表示为：
[0074][0075]
式中，e表示该等效弹簧所对应材料的杨氏模量，r表示杆的半径，c表示离散单元轴向长度；
[0076]
而两种不同种材料之间的弹簧等效刚度z则表示为：
[0077][0078]
式中，e1表示第一种材料的杨氏模量；e2表示第二种材料的杨氏模量；
[0079]
另外，弹簧-质点系统中的第j个质点质量mj表示
[0080]
为：mj＝πρjcr2(3)
[0081]
式中，ρ表示离散单元密度，ρj表示第j个质点对应的离散单元的单元密度；
[0082]
其次，根据牛顿第二定律可知弹簧-质点系统中第j个质点的运动方程：
[0083]mjaj
＝zj(v
j 1-vj)-z
j-1
(v
j-v
j-1
)(j＝1,2,...,n)(4)
[0084]
式中，aj表示第j个质点的加速度，等于位移的二阶导，v表示位移。
[0085]
进一步的，在所述步骤(3)中，将弹簧-质点系统的振动方程与bloch定理结合如下
式所示：
[0086][0087]
式中，u表示振幅，i表示虚数的单位，q表示波矢，l表示杆的长度，n表示离散单元数，ω表示频率，t表示时间；
[0088]
上述得到的振动方程(4)可变换得：
[0089][0090]
单个泵体阀杆元胞假设离散成n个质点，在周期性边界条件下，两元胞接合面处的弹簧等效刚度可表示为一元胞的最后一个弹簧等效刚度zn或者表示为另一元胞的第0个等效弹簧刚度z0，z0＝zm；同理z1＝z
n 1
；同样应用到质点质量m0＝mm，m1＝m
n 1
；对于振幅也同样适用这一周期性边界条件，则u0＝um，u1＝u
n 1
；
[0091]
将振动方程(6)同周期性边界条件结合可得一线性方程组，展开为矩阵形式表示为：
[0092]
[ω2i-x(q)]u＝0(7)
[0093]
其中，i表示单位矩阵，x(q)表示关于材料参数和结构参数以及波矢和频率有关的矩阵，u表示是由各质点振幅组成的列向量；
[0094]
令矩阵x(q)-ω2i行列式为零，可得阀杆主体部分单胞的色散关系方程：
[0095]an
ω
2n
a
n-1
ω
2(n-1)

…
a1ω2 a0＝0(8)
[0096]
式中，a0～an表示与q、r、e、ρ、l、n相关的待定系数；
[0097]
其中，当色散关系方程(8)不等于0时，则定义其为色散多项式：
[0098]
q(ω2)＝anω
2n
a
n-1
ω
2(n-1)

…
a1ω2 a0(9)
[0099]
设式(7)中的波矢q由实部和虚部来表示：
[0100]
q＝g bi(10)
[0101]
式中，g和b均表示为实数，i表示虚数符号；当b不为零时则表示衰减存在，b的数值大小表示该结构衰减程度的大小，故b被称为衰减因子。
[0102]
进一步的，在步骤(4)中，所述明确色散多项式和减振频率以及色散多项式和衰减因子的对应关系是：
[0103]
首先，将数值带入到色散多项式并进行减振频率的求解计算；
[0104]
接着，通过控制变量法分析色散多项式和指定减振频率、衰减因子的对应关系，从而明确色散多项式和减振频率以及色散多项式和衰减因子的对应关系。
[0105]
进一步的，在步骤(5)中，所述建立满足指定减振频率及衰减因子的逆向设计存在性判据如下式所示：
[0106][0107]
其中，b
max
表示是需要满足的衰减因子大小；
[0108]
为满足波矢q在取值范围内的取值要求，波矢取值后色散多项式应满足如下式要
求：
[0109][0110]
进一步的，在步骤(6)中，引用sturm定理结合存在性判据(11)和(12)进行求解符合目标减振频率及衰减因子要求的阀杆体结构或材料参数；所述的sturm定理是一种确定实数函数根的计算方法，首先满足判据(11)和(12)需将结构参数和材料参数以及ω2＝x带入色散多项式，色散多项式q(x)要满足q∈[0,π/l]以及衰减区间范围内sturm序列变号数等于0。相对应sturm序列变号数等于0说明色散多项式q在指定的频率范围内无根，则相应的参数即为满足要求的参数；
[0111]
具体的，为确保指定衰减频率在波矢q的取值范围内没有根，应用sturm定理进行计算，可有效的节省计算时间和确保计算结果的准确性；将ω2＝x带入色散多项式，有q(ω2)＝q(x)，根据sturm定理，可以给出一个色散多项式的sturm序列，取q0(x)＝q(x)，q1(x)＝q`(x)，对q1(x)和q0(x)应用带余除法，得q0(x)＝q1(x)q1(x)-q2(x)，再对q2(x)和q1(x)应用带余除法，得q1(x)＝q2(x)q2(x)-q3(x)，依次类推，直至有qi(x)＝0；计算时将取值区间的最小值和最大值带入色散多项式q，分别求得最小值和最大值相对应的q的sturm序列，即q0(min)、...、qn(min)和q0(max)、...、qm(max)；若相邻两个序列正负异号则为变号，序列q0(a)，q1(a)，
…
，qm(a)的变号数v(a)与序列q0(b)，q1(b)，
…
，qm(b)的变号数v(b)的差v(a)-v(b)恰是f(x)在区间(a，b)内实根的个数；为确保在区间内色散多项式q无根，相对应的q(min)和q(max)的变号数要相等即v(min)-v(max)＝0。
[0112]
设本实施实例中需要满足减振频率为15000hz、衰减因子大于等于20；以一维等截面二组元杆为例，已知结构参数和b组元的材料参数，求a组元的材料参数即杨氏模量和密度。
[0113]
如图2所示，本发明提供的一种泵体阀杆的减振设计，
[0114]
包括以下步骤：
[0115]
步骤一、对阀杆主体以等截面二组元声子晶体杆结构进行离散，元胞离散为存在十个质量点的弹簧质点系统，设阀杆主体长为1m；
[0116]
步骤二、根据步骤3建立色散多项式方程；色散多项式方程的推导由下列公式表示：
[0117]
表1材料参数
[0118][0119]
杆单胞长度l＝0.1m，r＝0.005m；根据等效原理，可得每段杆的等效质量：
[0120][0121]
式中，ρa表示为待求材质的密度，ρb表示为钢材质的密度；
[0122]
相邻质点间在轴向上的刚度系数为：
[0123][0124]
式中，ei、e
i 1
分别表示为第i个与第i 1个离散单元的杨氏模量；
[0125]
如图2所示，图2是弹簧质点链振动方程：
[0126]mjaj
＝zj(v
j 1-vj)-z
j-1
(v
j-v
j-1
)(j＝1,2,...,10)(15)
[0127]
步骤三、根据上述振动方程结合bloch定理有：
[0128][0129]
其中：n＝10；将位移方程和式(16)结合后，可变换得：
[0130][0131]
展开得矩阵形式为：
[0132]
(ω2i-x(q))u
t
＝0(18)
[0133][0134]
式中，第一种材料结构的离散质点质量m1～m5等于第二种材料结构的离散质点质量m6～m
10
等于相对的左右质点属于同种材料的等效弹簧刚度k1～k4等于k6～k9等于左右质点不属于同种材料的等效弹簧刚度k5、k
10
等于
[0135]
u＝[u1…u10
](20)
[0136]
步骤四、由矩阵ω2i-x(q)行列式为零，可得色散关系方程：
[0137]a10
ω
20
a9ω
18

…
a1ω2 a0＝0(21)
[0138]
式中：a0～a
10
为与q、r、e、ρ、l、n相关的待定系数；
[0139]
色散多项式为：
[0140]
q(ω2)＝a
10
ω
20
a9ω
18

…
a1ω2 a0(22)
[0141]
步骤五、建立逆向设计的存在性判据。根据步骤二得到的色散多项式建立逆向设计的存在性判据如下：
[0142]
且
[0143]
式中，ω0＝2πf；f＝15000，b
max
＝20。
[0144]
下表为计算结果：
[0145]
表2
[0146][0147][0148]
结果体现，通过本发明计算得出的参数作为元胞的材料参数可以对指定频率15000hz产生衰减，并且由图4和图5可知衰减程度耶同样满足即衰减因子数值足够大。
[0149]
步骤六、结合sturm定理进行编程计算，求解出第一种材料的杨氏模量ea和密度ρa；最后选取结果点集中的数据进行带隙计算验证。
[0150]
将求解所得参数中的一个数据点(杨氏模量为1.5
×
109pa、密度为1000kg/m3)作为a段的材料参数，然后计算得到杆元胞结构的衰减因子图，如图6所示，纵坐标表示频率，横坐标表示衰减因子b；为体现求解数据的可靠性，再将不是求解所得参数的一个数据点(杨氏模量为1.5
×
109pa、密度为2500kg/m3)作为a段的材料参数，然后计算得到杆元胞结构的衰减因子图，如图7所示，纵坐标表示频率，横坐标表示衰减因子b；利用comsol有限元仿真软件对元胞在带隙内的某一频率处的振动模态进行仿真验证。令元胞a段的杨氏模量为1.5
×
109pa；密度为1000kg/m3；在5000hz处的仿真结果如图8所示、在15000hz处的仿真结果如图9所示，可以看出对在5000hz和15000hz频率的元胞一端施加激励，而元胞末端没有出现明显位移，说明在5000hz处减振杆减振效果好。令元胞a段的杨氏模量为1.5
×
109pa；密度为2500kg/m3；在5000hz处的仿真结果如图10所示、在12000hz处的仿真结果如图11所示，可以看出对在5000hz和12000hz频率的元胞一端施加激励，而元胞末端没有出现明显位移，说明在5000hz和12000hz处减振杆减振效果好。
[0151]
以上所述，仅为本发明较佳的具体实施方式，但其保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。