2014年10月 农业机械学报第45卷 第10期 d o i : 1 0 . 6 0 4 1 / j . i s s n . 1 0 0 0 1 2 9 8 . 2 0 1 4 . 1 0 . 0 5 1 基于流体力学泄漏模型的螺杆泵泄漏机理分析 张元勋1 唐倩2 李忠华2 严迪2 ( 1 . 重庆大学航空航天学院,重庆 4 0 0 0 4 4 ;2 . 重庆大学机械传动国家重点实验室,重庆 4 0 0 0 4 4 ) 摘要:针对螺杆泵筒壁间隙、 啮合间隙等决定螺杆泵性能的关键间隙设计问题, 在研究螺杆各项间隙构成原理的基 础上, 运用流体力学间隙泄漏理论, 建立泵腔间隙泄漏压差流和剪切流模型, 分别分析螺杆泵筒壁间隙和啮合间隙 的泄漏量, 得到螺杆泵不同间隙泄漏的表达式.以主动螺杆 2头、 从动螺杆 3头的双螺杆泵为例, 通过理论计算与 试验对比分析了螺杆转速、 压差对泵泄露的作用机理, 揭示了螺杆转速、 泵进出口压差对螺杆泵流量特性及容积效 率特性的影响规律, 从而完善了螺杆泵转子设计理论, 提升了螺杆泵输送性能. 关键词:螺杆泵 泄漏 间隙 容积效率 中图分类号:T H 3 2 7 文献标识码:A 文章编号:1 0 0 0 1 2 9 8 ( 2 0 1 4 ) 1 0 0 3 2 6 0 7 收稿日期:2 0 1 3 0 9 2 6 修回日期:2 0 1 3 1 0 2 7 国家自然科学基金 资助项目(51175522)、教育部创新团队发展计划资助项目(IRT1196)和重庆市自然科学基金(重点) 资助项目(cstc2013jjB70002)作者简介:张元勋, 博士生, 主要从事螺杆泵和水泵等泵类产品基础理论研究, E m a i l :y u a n x u n . z h a n g @y a h o o . c o m 通讯作者:唐倩, 教授, 博士生导师, 主要从事机械传动、 流体机械及计算机辅助仿真研究, E m a i l :t q c q u @c q u . e d u . c n 引言 螺杆泵作为容积式泵系列中的典型产品, 具有 结构简单、 输送介 质粘度范围 广、 液 力脉 动小 等优 点[ 1- 2 ] .多年来, 螺杆泵的研究者分别在螺杆齿形 曲线设计优化[ 3- 5 ] 、 螺杆加工刀具设计[ 6- 8 ] 及螺杆 泵密封性能[ 9- 1 2 ] 等方面开展了诸多研究.然而, 由 于螺杆齿面构型的复杂性使得螺杆泵的间隙控制相 比其他类型容积式泵具有更大的不确定性和不均匀 性[ 1 3-1 4 ] , 严重影响螺杆泵的容积效率和使用寿命. 螺杆间隙的波动性直接决定螺杆泵的性能及寿 命, 间隙过大导致 泵的内泄漏量增 加, 容积 效率降 低; 间隙过小则导致运转部件间的摩擦增加, 使用寿 命降低[ 1 5 ] , 因此, 螺杆间隙的合理设计备受学者们 重视[ 1 6- 1 8 ] .但是, 对于间隙的构成及其影响没有形 成系统的理论; 普遍采用样机试验的方法, 对不同间 隙条件下泵的性能进行预测, 研发效率低, 资源浪费 严重.因此, 深入研究螺杆泵间隙对泵输送性能的 作用机理, 获得螺杆泵各项间隙的构成及其对泵性 能的影响规律, 是螺杆泵性能提高的关键. 本文在研究螺杆各项间隙构成原理的基础上, 对螺杆间隙进行分类, 针对螺杆泵筒壁间隙、 啮合间 隙等决定螺杆泵性能的关键间隙的设计问题, 运用 流体力学间隙泄漏理论, 分别建立泵腔间隙泄漏压 差流和剪切流模型, 分析螺杆泵筒壁间隙和啮合间 隙的泄漏量, 得到螺 杆泵不同 间隙泄漏 的表 达式. 以主动螺杆 2头、 从动螺杆 3头双螺杆泵为例, 分别 分析螺杆转速、 泵进出口压差对螺杆泵流量特性及 容积效率特性的影响, 并与试验结果进行对比, 验证 模型的正确性. 1 螺杆泵间隙构成 螺杆泵工作过程中泄漏间隙主要有: ①螺杆转 子齿顶与泵缸筒壁之间形成的筒壁 间隙 δ 1 .② 啮 合区主动螺杆齿顶与从动螺杆齿根或者主动螺杆齿 根与从动螺杆齿顶所形成的径向间 隙δ2.③ 啮合区螺杆齿面之间沿接触线均匀分布的法向间隙, 将 其沿圆周投影到轴截面即为法向间隙 δ 3 ( 图1).其 中, 由螺杆啮合线及筒壁所形成的泄漏三角形, 由于 其与法向间隙及筒壁间隙相连, 且靠近啮合区, 因此 将泄漏三角形间隙合并到法向间隙. 1 1 筒壁间隙构成 螺杆泵靠相互啮合的螺杆绕各自的轴线旋转产 生容积变化将液体从一端吸入, 输送到另一端排出. 研究发现, 螺杆的筒壁泄漏间隙有效泄漏长度与泵 缸缸筒长度并不相等.如图 2所示, 在低压区, 由于 筒壁间隙内液体与进口液体相连, 筒壁间隙内液体 压力等于进口压力; 同理, 在高压区, 齿槽内液体与 出口压力相等.因此, 在螺杆进口端与出口端的某 一范围内, 筒壁间隙不发生泄漏. 图1螺杆泵泄漏间隙示意图 F i g . 1 D i a g r a m m a t i cs k e t c ho f c l e a r a n c e o f p u m pl e a k a g e 图2筒壁间隙泄漏长度 F i g . 2 L e n g t ho f t o pc l e a r a n c eo f t o o t h 过渡区即为有效泄漏区.随着螺杆的转动, 液 体将沿着轴向流向泵的排出端, 过渡区液体压力从 进口端到出口端逐渐增加.螺杆每转一周, 过渡区 将会重复一次上述运动.由于主动螺杆是双头的, 因此过渡区长度即为螺杆螺距的 2倍, 即主动螺杆 导程, 称过渡区为有效泄漏区域. 1 2 法向间隙构成 法向间隙又称为啮合间隙, 是指螺杆齿面在啮 合过程中接触线所形成的间隙.螺杆在相互啮合运 动中, 主动螺杆与从动螺杆齿面相接触所形成的空 间曲线即为螺杆的理论接触线.为了避免螺杆的接 触磨损造成功率损失, 实际加工过程往往要留出一 定的间隙, 该间隙即为沿着接触线分布的主、 从螺杆 齿面的法向距离, 因此又称螺杆法向间隙.径向间 隙理论上是螺杆的齿顶与齿根在运动中接触线所形 成的间隙, 严格意义上属于法向间隙.为了计算方 便, 本文将径向间隙归为法向间隙处理( 图3).泄漏三角形理论上是由主动螺杆、 从动螺杆接 触线的顶点( 该点不在齿顶圆上) 与其分别对应的 筒壁的接触点( 该点为齿顶圆与筒壁的交点) 所形 成的一个空间曲边三角形( 图3).将该空间曲边三 角形沿圆周方向投影到轴截面即为图 4 a所示三角 形MO N . 螺杆泵在 运转过程中, 理 论啮 合间隙为零, 此时, 泵的主要泄漏发生在泄漏三角形区域( 图4a中 MN O区域) .为避免螺杆啮合时功率损失, 在加工 螺杆时预留一定的间隙, 使得齿面在满足理论齿形 图3螺杆接触线示意图 F i g . 3 D i a g r a m o f c o n t a c t l i n eo f s c r e wr o t o r s 图4泄漏三角形演化示意图 F i g . 4 D i a g r a m o f e v o l u t i o no f l e a k a g et r i a n g l e ( a )泄漏三角形理论廓形 ( b )泄漏三角形实际廓形 共轭的同时, 保证其实际运转中不接触, 螺杆齿面接 触线也将沿着各自螺杆齿面法向向内偏移成为两条 空间的等距曲线, 两空间曲线间的距离即为螺杆齿 面的法向间隙.进一步, 螺杆齿顶接触点 O也将分 为主动螺杆齿顶接触点 O m 和从动螺杆齿顶接触点 O f , 原来的泄漏三角形也变为图 4所示曲边四边形 N MO mO f .由于啮 合间 隙的存在, 原来封闭的泄漏三角形与螺杆齿面法向间隙连通, 使其联合构成螺 杆齿面法向间隙泄漏. 2 螺杆泵内泄漏模型 2 1 筒壁间隙泄漏计算模型 由11节分析可知, 螺杆筒壁间隙泄漏只发生 在螺杆有效泄漏区域, 而且, 该区域液体也随着螺杆 的转动, 由进口端面移动到出口端面.其有效泄漏 长度为螺杆导程 S .进一步分析螺杆及泵缸缸筒几 何关系( 图5)可以得到, 螺杆的筒壁间隙泄漏可以 理解为如下情形: 一个∞形物体在一个∞形筒壁内 沿轴向匀速移动, 其形成的间隙是均匀的, 且该物体 在筒内的运动方向是沿轴向由低压方向运动到高压 方向, 该物体在筒内运动过程中所形成的间隙泄漏 7 2 3 第10期 张元勋 等:基于流体力学泄漏模型的螺杆泵泄漏机理分析 量即为螺杆筒壁间隙导致的泄漏量. 图5筒壁间隙示意图 F i g . 5 D i a g r a m m a t i cs k e t c ho f t o pc l e a r a n c eo f t o o t h 如图 5所示筒壁间隙为 δ 1 , 将该间隙沿某一方 向展开到平面上, 那么该筒壁间隙泄漏量的计算即 可转换为两平行平板间间隙的泄漏量.平板间隙中 的流动一般为层流, 分 3种类型: 压差造成的流动称 压差流动、 相对运动造成的流动称剪切流动和压差 与剪切同时作用 下的流动.螺杆 泵 是一 种容 积式 泵, 将液体由低压区输送到高压区, 结合前面分析, 螺杆泵的筒壁间隙泄漏是包含压差流和剪切流的混 合流动, 其中, 压差流来自泵进出口的压差; 剪切流 则由螺杆的旋转带动螺杆齿顶沿 壁面 的轴 向移动 ( 速度 u 1 ) 引起. 图6平行平板间隙流动示意图 F i g . 6 D i a g r a m o f p a r a l l e l p l a t ef l o w 由上述分析可知, 假定图 6中平板 A为螺杆齿 顶有效泄漏区移动构成的柱面所展开的平板, 平板 B为对应筒壁所展开的平板, 其中平板 B相对静止, 平板 A相对平板 B有相对移动速度 u 1 , 平板宽度均 为b,平板间隙 h=δ 1 , 有效长度 l 为主动螺杆的导 程.假定液体不可压缩, 质量力忽略不计, 粘度不 变.在液体中取一个微元体 d x 、 d y ( 宽度方向取单 位长度) , 作用在它与液流相垂直的 2个表面上的 压力为 p 和p+d p , 作用在它与液流相平行的上下 两个表面上的切应力为 τ和τ+d τ , 因此微元体的 受力平衡方程为 p d y +( τ+d τ ) d x =( p +d p ) d y +τ d x ( 1 ) 式中 p — — —压力 τ — — —切应力 结合牛顿内摩擦定律并积分可得 u= y 2 2 μ d p d x +C 1 y +C 2 ( 2 ) 式中 C 1 、 C 2 — — —积分常数 μ— — —动力粘度 由平板 间隙流动方式并结合以上分析可知, 图 6所示间隙流动是压差流和剪切流的叠加.因此 对式( 2 ) 分以下两种情况讨论. ( 1 ) 压差流 平板 A相对平板 B静止, 在间隙两端只有压差 存在.根据边界条件 y= 0时, u=0 ; y=h时, u=0 可以得到 C 1 =- h 2 μ d p d x C 2 = 0 于是由压差引起的平板间隙泄漏量为 q p = b h 3 1 2 μ l Δ p ( 3 ) 从式( 3 ) 可以看出, 通过间隙的流量与间隙的 三次方呈正比. 在压差作用下, 流过固定平行平板缝隙的流量 与缝隙宽度 b 、 压力差 Δ p 、 缝隙厚度 h的三次方呈正 比, 而与液体的动力粘度、 缝隙长度 l 呈反比. ( 2 ) 剪切流 两平行平板有相对运动, 速度为 u i , 无压差, 这 种流动称 为纯剪切流动.根据边界条件, 当y=0 时, u=0 ; 当y=h时, u=u i , 且dp/dx=0 , 代入式( 2 ) 得C1=uihC2=0因此, 可以得到螺杆泵在工作过程中由于筒壁 间隙而导致的泄漏量为 q w = b h 3 1 2 μ l Δ p ± b h 2 u i ( 4 ) 其中 b = R 1 +R 2 2 =R M h=δ 1 u i= 2 π n R Mt a n λ 式( 4 ) 中正、 负号的确定: 当长平板 B相对于短 平板 A的运动方向和压差流动方向一致时, 取"+" 号; 反之取"-" 号.本文中泵缸缸筒相对螺杆齿顶 的移动方向与压差方向相同, 结合螺杆泵几何参数 可以得到螺杆泵筒壁间隙泄漏计算模型 q w = R Mδ 3 1 1 2 μ S ( p 1 -p 2 )+π n R 2 Mδ 1 t a n λ ( 5 ) 2 2 法向间隙泄漏计算模型 由12节分析螺杆法向间隙的形成原理, 结合 螺杆运动过程中啮合区的几何特性, 可以将螺杆法 向间隙泄漏看作是局部管道截面突然缩小的管道流 8 2 3 农业机械学报2014年动.当流体从大直径管( 即螺杆高压腔) 流入小直 径管( 即螺杆法向间隙) 过程中, 流线弯曲, 流束截 面收缩; 进入小直径管后, 由于惯性, 流束继续收缩 至最小截面, 即为本文所述螺杆法向间隙, 而后逐渐 增大流入螺杆低 压腔中, 构成 螺杆 法向间 隙泄 漏. 由于螺杆法向间隙是由均匀的等距空间曲线构成, 其法向啮合间隙将构成一个空间的曲边矩形. 综上所述, 螺杆法向间隙泄漏可以看成是下述 情形: 主、 从螺杆啮合区的法向间隙 δ 3 随着螺杆的 轴向移动由低压腔逐渐移动到高压腔, 如图 7所示, 物体 M 即为螺杆啮合齿面, 物体 N为泵缸缸筒, M 相对 N有一轴向速度 u 1 , p 1 为高压区压力, p 2 为低 压区压力, 由于螺杆的啮合运动, 结合泵缸缸筒的密 封作用, 将高低压腔液体分割开, 并由物体 M 的轴 向移动将高压腔液体输送到出口处.在此过程中, 螺杆法向间隙的泄漏假定为矩形薄壁小孔, 其中小 孔的高度即为法向啮合间隙 δ 3 , 小孔的长度即为螺杆 齿面接触线长度与泄露曲边四边形的径向长度之和. 图7小孔泄漏模型示意图 F i g . 7 D i a g r a m o f l e a k a g em o d e l f o r s m a l l h o l e s 根据以上分析, 对于图 7所示的通过薄壁小孔 的液流, 取截面 1 1和22为计算截面, 设截面 1 1处的压力和平均速度分别为 p 1 、 u 1 , 截面 2 2 处的压力和平均速度分别为 p 2 、 u 2 .选取同一高度 处液体, 伯努利方程为 p 1 +ρ g h 1 + 1 2 ρ α 1 u 2 1 =p 2 +ρ g h 2 + 1 2 ρ α 2 u 2 2 +p W ( 6 ) 其中 p W =K ρ u 2 2 2 ( 7 ) 式中 ρ — — —液体密度 p W— — —局部压力损失 K — — —修圆进口局部损失因数 α 1 、 α 2 — — —截面 1 1 、 2 2液体动能修正系 数, 即截面上实际动能与平均流动 计算的动能之比值 由于液体在流经小孔时, 其压力能损失主要是 液体断面的突然变小而导致的局部损失, 液体粘度 在该情况下可以忽略.因此, 各截面动能修正系数 在此处取值为 1 . 关于相对运动速度的计算方法可参考齿轮啮合 原理[ 1 9 ] , 此处不再赘述.将式( 7 ) 代入式( 6 ) 并化 简得 u 2 = p 1 -p 2 + 1 2 ρ u 2 1 1 2 ρ ( 1+K 槡)(8)结合表 1所列修圆进口局部损失因数, 以及螺 杆齿面及法向间隙几何特性, 可以确定本文中修圆 进口局部损失因数为 0 0 4 . 表1修圆进口的局部损失因数 T a b . 1 P a r t i a l l o s s f a c t o ro f t i pr o u n d i n gi mp o r t r / D 0 0 2 0 0 6 ≥0 1 5 K 0 2 8 0 1 5 0 0 4 通过以上分析, 可以得出螺杆泵法向间隙泄漏 量为 q n =( l c+l t ) δ 3 p 1 -p 2 + 1 2 ρ u 2 1 1 2 ρ ( 1+K 槡)(9)式中 l c — — —螺杆接触线长度[ 2 0 ] l t — — —泄漏曲边四边形径向方向长度 2 3 实例计算模型 对于螺杆泵法向间隙, 由于螺杆啮合区几何特 性复杂, 且多为空间分布, 不能直观判断螺杆法向间 隙的构成.这里, 将图 8 a 中从动螺杆沿着半径为 r ′ 的圆展开到一个平面上, 其中半径为 r ′ 的圆是一个 与AB线 相切的圆. 其剖面图如图8b所 示. 在图8c中, 空白区域为从动螺杆齿槽部分, M、 N分别 对应主动螺杆的两头. 以从动螺杆为相对静止状态, 考察主动螺杆的 运动情况.主动螺杆的旋转运动在图 8 c中可以看 作是主动螺杆齿顶 M、 N的x与z方向运动的合成, 即主动螺杆转动一个角度, 在图 8 c中体现为 M、 N 由位置Ⅰ Ⅰ运动到Ⅱ Ⅱ.其中包含着螺杆的圆周 运动( x 方向移动) 和轴向运动( z 方向移动) .进一 步可以得到, 螺杆在相互运动过程中, 泄漏三角形的 大小及其相对主动螺杆齿顶的位置没有发生变化; 主动螺杆与从动螺杆在半径 r ′ 的圆上接触点位置相 对从动螺杆( 地面坐标系) 来说, 则是沿着 z 方向逐 渐移动到出口.同理, 在螺杆齿槽不同半径圆上的 接触点, 其相对地面坐标系仅做轴向运动( z 方向移 动) , 彼此相邻的接触点则处于相对静止, 这样就构 成隔离螺杆高、 低压腔的分割线 ( 即螺杆齿面接触 线) . 结合图 3和8c可以看出, 当螺杆泵工作时, 主、 从动螺杆啮合区将形成稳定的接触线, 其有效作用 9 2 3 第10期 张元勋 等:基于流体力学泄漏模型的螺杆泵泄漏机理分析 图8泵腔内流场特性示意图 F i g . 8 D i a g r a m o f f l o wc h a r a c t e r i s t i c si np u m pc a v i t y 的长度即为主动螺杆 M、 N作用的轴向区域.在该 区域内, M、 N分别 分割着一个螺 距内从动螺杆的高、 低压腔, 且由于主动螺杆是两头, 因此, 该螺杆泵 接触线的有效长度即为主动螺杆的导程. 图 8所示泄漏三角形分别分布在螺杆齿面的上 下两侧, 且与泵缸缸筒相连.从图中也可以明显看 出, 分割螺杆泵齿槽高、 低压腔的泄漏三角形个数为 主动螺杆头数的两倍, 即在该型螺杆泵中有效泄漏 三角形个数为 4个. 在获得螺杆啮合区几何特性及密封特性后, 即 可用式( 9 ) 计算螺杆泵法向间隙泄漏.以某一双螺 杆泵为例, 计算分析压差、 转速、 导程及间隙等对螺杆 泵流量、 容积效率的影响关系.各参数如表 2 所示. 表2螺杆几何参数 T a b . 2 Ge o me t r yp a r a me t e r so f s c r e w r o t o r s 参数 数值 主动螺杆 从动螺杆 头数 2 3 节圆半径 / m m 4 2 6 3 根圆半径 / m m 3 5 3 5 顶圆半径 / m m 7 0 7 0 螺距 / m m 6 1 6 1 中心距 / m m 1 0 5 接触线长度 / m m 1 5 5 泄漏三角形等效长度 / m m 7 5 3 试验及讨论 如图 9所示, 由电动机经扭矩、 转速仪将扭矩传 递到待试验螺杆泵上, 油路经过滤器进入到螺杆泵, 泵出口依次连接流量计、 过滤器并经冷却器回流油 箱.其中, 阀门控制泵进油口流量, 各仪器仪表测量精 度依次为: 转速表 ± 0 5 %、 压力表不低于 1 5级、 温度 表±1℃、 流量计不大于 0 5 %、 转矩仪小于 ± 1 5 %. 图9螺杆泵试验装置简图 F i g . 9 D i a g r a m m a t i cs k e t c ho f s c r e wp u m pt e s t u n i t 1 . 泵2.扭矩、 转速仪 3 . 电动机 4 . 过滤器 5 . 冷却器 6 . 油箱 7 . 过滤器 8 、 1 1 . 阀门 9 . 真空表 1 0 . 压力表 1 2 . 流量计 试验用油 为CD40润滑 油, 密度889k g / m 3 , 油温50~ 7 0 ℃, 动力粘度 5 2 5* 1 0 - 2 P a ·s . 在计算螺杆泵理论流量及容积效率时, 螺杆法 向间隙与筒壁间隙均为均匀变化, 各间隙宽度不随 螺杆的转动而变化, 是一种理想状况.然而, 由于螺 杆加工误差的存在, 不能保证各项间隙均匀变化, 导 致螺杆在啮合运动过程中法向间隙与筒壁间隙随螺 杆的运转而产生周期性的变化.由于螺杆间隙周期 性变化的存在, 使得螺杆在运转过程中泄漏量也随 着产生变化, 且与计算结果存在一定差距. 图10所 示为初始进口压力50k P a , 转速为2100r/min,导程 1 2 2 m m , 螺杆齿面法向间隙 0 1 2 m m , 筒壁间隙 0 2 4 m m , 理论流量为 1 6 4 1 5 m 3 / h时试验 与计算的结果对比.当出口压力低, 即压差小时, 计 算结果与试验结果偏差较小.随着压差的增加, 在 间隙一定的条件下, 螺杆泵泄漏量将随之增加.说 明压力对螺杆的容积效率影响很大, 进一步可以得 到, 压差流在螺杆泵泄漏中占主导地位. 图10泵流量 压力 容积效率特性曲线 F i g . 1 0 C h a r a c t e r i s t i cc u r v eo f r a t eo f f l o w p r e s s u r e v o l u m e t r i ce f f i c i e n c y 由图 1 0可知, 随着压差的增加, 计算结果与试 0 3 3 农业机械学报2014年 验结果的偏差也越来越大, 正如前面分析, 计算是在 间隙均匀的条件下完成的, 而试验所得结果是在螺 杆间隙不均匀条件下完成的, 这就说明, 螺杆加工间 隙的误差在低压条件下表现不明显, 当压差增加时, 其间隙的不均匀性也将凸显出来. 螺杆泵理论流量为 q i d e a l=S f l u i d L m a i n n m a i n 式中 q i d e a l — — —螺杆泵在转速为 n ma i n 时的理论流量 S f l u i d — — —泵腔流道截面面积 L ma i n — — —主动螺杆导程 n ma i n — — —主动螺杆转速 由上述分析可知, 螺杆啮合间隙的不均匀性在 转速恒定的条件下, 随着压力的增加将表现的愈来 愈明显; 当螺杆压差恒定时, 螺杆间隙对泵流量特性 的影响也将随转速的变化而变化( 图11,试验条件 为初始进口压力为 5 0k P a , 出口压力 8 5 0k P a , 导程 1 2 2m m , 螺杆齿面法向间隙012m m , 筒壁间隙024m m) .由图中流量曲线可以得知, 在低转速(小于 1 9 0 0r / m i n ) 时, 计算结果与试验结果差别较 小, 即在低转速条件下, 以润滑油作为输送介质的螺 杆泵而言, 其对螺杆由于加工而引起的间隙的不均 匀性表现不是很敏感, 当转速超过 19 0 0r / m i n时, 试验所得螺杆泵的泄漏量相比计算值将越来越大, 其容积效率也随之急剧下降. 图11泵流量 转速 容积效率特性曲线 F i g . 1 1 C h a r a c t e r i s t i cc u r v eo f r a t eo f f l o w r o t a t i o n s p e e d v o l u m e t r i ce f f i c i e n c y 同时, 随着转速的增加, 螺杆相对筒壁的轴向速 度也随之增加, 由于液体粘性的作用, 泵体对泵腔间 隙处液体的剪切作用也将越来越明显, 由于泵体对 液体的剪切作用而形成的剪切流也将随着螺杆转速 的增加而增加.如图 1 1所示容积效率曲线, 当转速 小于 1 0 0 0r / m i n时, 试验与计算得到的容积效率都 迅速增加, 此时螺杆由于加工导致的不均匀间隙对 螺杆的 泄漏量影响不大; 当转速大于10 0 0r / m i n 时, 计算结果显示容积效率增加速度低于前者, 而试 验结果偏离计算结果的幅度越来越大, 导致容积效 率在转速大于 19 0 0r / m i n时出现明显的拐点而急 剧下降.因此, 对于特定的输送介质而言, 过高的转 速并不利于螺杆泵容积效率的提高. 为了更加直观判定所建泄漏模型的合理性, 通 过分析图 1 0与图 1 1数据计算得理论值与实测值的 误差( 图12).从图 1 3 ( 试验条件为初始进口压力 为50k P a , 导程为122m m , 螺杆齿面法向间隙012m m , 筒壁间隙 0 2 4m m) 中可以看出, 在螺 杆 泵额定转 速范围内, 理论值与实测值误差不超过6%, 验证了本文所建模型的合理性. 图12理论计算值与实际测量值误差 F i g . 1 2 E r r o r c u r v eo f t h e o r ya n dm e a s u r e m e n t 图13泵转速 压力 容积效率特性曲线 F i g . 1 3 C h a r a c t e r i s t i cc u r v eo f r o t a t i o n s p e e d p r e s s u r e v o l u m e t r i ce f f i c i e n c y 对比图 1 3不同转速均匀间隙条件下螺杆泵容 积效率随压差的变化关系, 可以得知, 螺杆泵容积效 率降低的速度将随着压差的增加而增大, 尤其表现 在低转速条件下.泵容积效率降低的速度随着压差 的增加而增加, 随着转速的降低而增加.说明在均 匀间隙条件下, 过低的转速将会使螺杆泵的做功有 一定的损耗, 降低了螺杆泵的总效率. 4 结论 ( 1 ) 分析了螺杆泵各项间隙的构成原理, 对螺 杆泵主要泄漏形式, 运用流体力学间隙泄漏及小孔 1 3 3 第10期 张元勋 等:基于流体力学泄漏模型的螺杆泵泄漏机理分析 泄漏计算方法, 分别建立了螺杆泵筒壁间隙泄漏模 型和法向啮合间隙泄漏模型, 得到螺杆泵理论泄漏 量数值计算方法. ( 2 ) 搭建了螺杆泵流量特性试验台, 理论与试 验对比分析了螺杆泵螺杆转速、 压差对容积效率的 影响, 获得了不同工况条件下螺杆间隙对螺杆泵流 量特性的作用规律. ( 3 ) 比较了螺杆泵主要工况下理论计算值与实 测值的大小, 二者结果较为一致, 且理论计算值最大 误差小于 6 %, 说明本方法是切实可行的. 参考文献1NakashimaC e l s oY ,J r S i l v i od eO l i v e i r a ,C a e t a n oEF .H e a t t r a n s f e r i nat w i n s c r e wm u l t i p h a s ep u m p :t h e r m a l m o d e l i n ga n do n e a p p l i c a t i o ni nt h ep e t r o l e u m i n d u s t r y [ J ] .E n e r g y ,2 0 0 6 ,3 1 ( 1 5 ) :3 4 1 5- 3 4 2 5 . 2 李福天.螺杆泵[ M] .北京:机械工业出版社,2 0 1 0 . 3 张波,魏益民,Wo l f g a n gS i e t z e , 等.双螺杆挤压机螺杆元件类型对扭矩和压力的影响[ J ] .农业机械学报,2 0 0 7 ,3 8 ( 9 ) : 7 1- 7 4 . Z h a n gB o ,We i Y i m i n ,Wo l f g a n gS i e t z e ,e t a l .I n f l u e n c eo f t w i n s c r e we x t r u d e rs c r e we l e m e n t t y p eo nt h et o r q u ea n dp r e s s u r e s d u r i n ge x t r u s i o n[ J ] .T r a n s a c t i o n so f t h eC h i n e s eS o c i e t yf o r A g r i c u l t u r a l Ma c h i n e r y ,2 0 0 7 ,3 8 ( 9 ) :7 1- 7 4 .( i nC h i n e s e ) 4 P e n n a c c h i P ,S e x t oLF .D e s i g ni m p r o v e m e n t o f s c r e wp u m pp o w e r s o u r c e s f o r h y d r a u l i ce l e v a t o r s t or e d u c en o i s ee m i s s i o n s [ J ] . N o i s eC o n t r o l E n g i n e e r i n gJ o u r n a l ,2 0 0 7 ,5 5 ( 2 ) :1 6 4- 1 7 1 . 5 潘中永,潘希伟,李晓俊,等.离心泵泄漏流转子动力学特性分析[ J ] . 农业机械学报,2 0 1 3 ,4 4 ( 1 1 ) :6 7- 7 1 . P a nZ h o n g y o n g ,P a nX i w e i ,L i X i a o j u n ,e t a l .R o t o r d y n a m i ca n a l y s i s o f l e a k a g ef l o wi nc e n t r i f u g a l p u m p[ J ] .T r a n s a c t i o n s o f t h e C h i n e s eS o c i e t yf o r A g r i c u l t u r a l Ma c h i n e r y ,2 0 1 3 ,4 4 ( 1 1 ) :6 7- 7 1 .( i nC h i n e s e ) 6 WuY u r e n ,F o n gZ h a n g h u a .I m p r o v e dr o t o rp r o f i l i n gb a s e do nt h ea r b i t r a r ys e a l i n gl i n ef o rt w i n s c r e w c o m p r e s s o r s[ J ] . Me c h a n i s m a n dMa c h i n eT h e o r y ,2 0 0 8 ,4 3 ( 6 ) :6 9 5- 7 1 1 . 7 王可,赵文珍,唐宗军, 等.异性螺杆无瞬心包络铣削技术研究 [ J ] .中国机械工程,2 0 0 1 ,1 2 ( 3 ) :2 9 4- 2 9 6 . 8 叶晓琰,石海峡,胡敬宁, 等.外环流转子泵转子与泵体间间隙的优化选择[ J ] .农业机械学报,2 0 0 7 , 3 8 ( 4 ) :8 1- 8 4 . Y eX i a o y a n ,S h i H a i x i a ,H uJ i n g n i n g ,e t a l .O p t i m i z a t i o ns e l e c t i o no f g a pb e t w e e nr o t o ra n dp u m pb o d yo f e x t e r n a l c i r c u l a t i o n r o t o r p u m p[ J ] .T r a n s a c t i o n so f t h eC h i n e s eS o c i e t yf o r A g r i c u l t u r a l Ma c h i n e r y ,2 0 0 7 ,3 8 ( 4 ) :8 1- 8 4 .( i nC h i n e s e ) 9 D z h a n a k h m e d o vAK h .I n v e s t i g a t i o n s o f t h ed e s t r u c t i v ee f f e c t o f t w o p h a s el i q u i do ns c r e wp u m pl i ps e a l s [ J ] .J o u r n a l o f F r i c t i o n a n dWe a r ,2 0 0 8 ,2 9 ( 4 ) :3 1 0- 3 1 3 . 1 0 李晓俊,袁寿其,潘中永,等. 基于结构化网格的离心泵全流场数值模拟[ J ] . 农业机械学报,2 0 1 3 ,4 4 ( 7 ) :4 9- 5 4 . L i X i a o j u n ,Y u a nS h o u q i ,P a nZ h o n g y o n g ,e t a l .N u m e r i c a l s i m u l a t i o no f w h o l ef l o wf i e l df o r c e n t r i f u g a l p u m pw i t hs t r u c t u r e d g r i d[ J ] .T r a n s a c t i o n so f t h eC h i n e s eS o c i e t yf o r A g r i c u l t u r a l Ma c h i n e r y ,2 0 1 3 ,4 4 ( 7 ) :4 9- 5 4 .( i nC h i n e s e ) 1 1 H s i e hC h i u f a n ,H w a n gY i i w e n ,F o n gZ h a n g h u a .S t u d yo nt h et o o t hp r o f i l ef o r t h es c r e wc l a w t y p ep u m p[ J ] .Me c h a n i s m a n d Ma c h i n eT h e o r y ,2 0 0 8 ,4 3 ( 7 ) :8 1 2- 8 2 8 . 1 2 熊伟, 冯全科. 双螺杆压缩机齿间间隙分布的计算[ J ] .西安交通大学学报,2 0 0 4 ,3 8 ( 7 ) :6 8 2- 6 8 5 . X i o n gWe i ,F e n gQ u a n k e .C a l c u l a t i o nt oi n t e r l o b ec l e a r a n c ed i s t r i b u t i o no f t w i n s c r e wc o m p r e s s o r[ J ] .J o u r n a l o f X i a nJ i a o T o n gU n i v e r s i t y ,2 0 0 4 ,3 8 ( 7 ) :6 8 2- 6 8 5 .( i nC h i n e s e ) 1 3 屈文涛, 胥掌世, 张弘, 等.潜油双螺杆泵泄露模型的理论研究[ J ] .石油钻探技术,2 0 0 7 ,3 5 ( 6 ) :7 6- 7 8 . Q uWe n t a o ,X uZ h a n g s h i ,Z h a n gH o n g ,e ta l .T h e o r e t i c a ls t u d yo nl e a k a g em o d e lo fs u b m e r s i b l et w i n s c r e w p u m p[ J ] . P e t r o l e u m D r i l l i n gT e c h n i q u e s ,2 0 0 7 ,3 5 ( 6 ) :7 6- 7 8 .( i nC h i n e s e ) 1 4 R y a z a n t s e vVM,P l y a s o vVV .D e t e r m i n i n gt h ef o r c e so nt h es c r e wi nt w o a n dt h r e e b e a r i n gt w o s c r e wp u m p s [ J ] .R u s s i a n E n g i n e e r i n gR e s e a r c h ,2 0 1 0 ,3 0 ( 9 ) :8 7 7- 8 8 5 . 1 5 刘在伦,王东伟,梁森.离心泵叶轮平衡孔液体泄漏量特性试验[ J ] . 农业机械学报,2 0 1 2 ,4 3 ( 7 ) :8 4- 8 8 . L i uZ a i l u n ,Wa n gD o n g w e i ,L i a n gS e n .F l u i dl e a k a g ec h a r a c t e r i s t i ct e s t o nb a l a n c ea p e r t u r eo f c e n t r i f u g a l p u m pi m p e l l e r [ J ] . T r a n s a c t i o n s o f t h eC h i n e s eS o c i e t yf o r A g r i c u l t u r a l Ma c h i n e r y ,2 0 1 2 ,4 3 ( 7 ) :8 4- 8 8 .( i nC h i n e s e ) 1 6 黎耀军,沈金峰,刘竹青,等. 轴流泵轮缘间隙非定常流动的大涡模拟[ J ] . 农业机械学报,2 0 1 3 ,4 4 ( 增刊) :1 1 3- 1 1 8 . L i Y a o j u n ,S h e nJ i n f e n g ,L i uZ h u q i n g ,e t a l .L a r g ee d d ys i m u l a t i o no f u n s t e a d yf l o wi nt i pr e g i o no f a x i a l f l o w p u m p[ J ] . T r a n s a c t i o n s o f t h eC h i n e s eS o c i e t yf o r A g r i c u l t u r a l Ma c h i n e r y ,2 0 1 3 ,4 4 ( S u p p . ) :1 1 3- 1 1 8 .( i nC h i n e s e ) 1 7 F o n gZ h a n g h u a ,H u a n gFC .E v a l u a t i n gt h ei n t e r l o b ec l e a r a n c ea n dd e t e r m i n i n gt h es i z e s a n ds h a p e so f a l l t h el e a k a g ep a t h s f o r t w i n s c r e wv a c u u m p u m p[ J ] .P r o c .I Me c h E , P a r t C :J o u r n a l o f Me c h a n i c a l E n g i n e e r i n gS c i e n c e ,2 0 0 6 ,2 2 0 ( 4 ) :4 9 9- 5 0 6 . 1 8 F o n gZ h a n g h u a ,H u a n gFC ,F a n gH S .E v a l u a t i n gt h ei n t e r l o b ec l e a r a n c eo f t h et w i n s c r e wc o m p r e s s o r b yt h eI S O c l e a r a n c e c o n t o u r d i a g r a m ( I C C D )[ J ] .Me c h a n i s m a n dMa c h i n eT h e o r y ,2 0 0 1 ,3 6 ( 6 ) :7 2 5- 7 4 2 . 1 9 F a y d o r LL i t v i n ,A l f o n s oF u e n t e s .G e a r g e o m e t r ya n da p p l i e dt h e o r y[ M] .L o n d o n :C a m b r i d g eU n i v e r s i t yP r e s s ,2 0 0 4 . 2 0 X uJ i a n ,F e n gQ u a n k e ,WuWe i f e n g .G e o m e t r i c a l d e s i g na n di n v e s t i g a t i o no f an e wp r o f i l eo f t h et h r e es c r e wp u m p[ J ] .A S ME J o u r n a l o f Me c h a n i c a l D e s i g n ,2 0 1 1 ,1 3 3 ( 9 ) :0 9 4 5 0 1 1- 5 . ( 下转第 3 3 9页) 2 3 3 农业机械学报2014年AC o mp o u n dC o n t r o l S t r a t e g yC o mb i n i n gL o a dC o mp e n s a t i o nw i t h A D R Co f E l e c t r o h y d r a u l i cP o s i t i o nC o n t r o l S y s t e m G a oB i n g w e i S h a oJ u n p e n g L i J i a n y i n g J i Y a j u a n ( C o l l e g eo f Me c h a n i c a l a n dP o w e r E n g i n e e r i n g ,H a r b i nU n i v e r s i t yo f S c i e n c e a n dT e c h n o l o g y ,H a r b i n1 5 0 0 8 0 ,C h i n a ) A b s t r a c t :T h ee l e c t r o h y d r a u l i cp o s i t i o nc o n t r o l s y s t e m ,w h i c hh a s t h es h o r t c o m i n g s o f h i g h l yn o n l i n e a r , t i m e v a r y i n go f t h ei n t e r n a l p a r a m e t e r sa n dd i s t u r b a n c eo f t h ee x t e r n a l l o a d ,s e r i o u s l ya f f e c t st h es t a t i c a n dd y n a m i cc o n t r o l .I n o r d e rt os o l v et h i sp r o b l e m ,ac o m p o u n d c o n t r o ls t r a t e g yc o m b i n i n gl o a d c o m p e n s a t i o nw i t hA D R C ( a u t od i s t u r b a n c er e j e c t i o nc o n t r o l l e r ) w a sp u tf o r w a r d ,a n dt h ew o r k i n g p r i n c i p l eo f t h ec o m p o s i t ec o n t r o l s t r a t e g yw a sg i v e ni nt h i sp a p e r .A D R Cc o n t r o l l e r w a sd e s i g n e d ,a n d t h e e x t e n d e d s t a t e o b s e r v e rw a su s e d f o ro b s e r v i n g i n t e r n a lp a r a m e t e r su n c e r t a i n t i e sa n d e x t e r n a l d i s t u r b a n c e s ,s o t h a tt h e d i s t u r b a n c e o ft h e s y s t e m i ss u p p r e s s e d e f f e c t i v e l y . L o a d c o m p e n s a t i o n c o n t r o l l e r w a sd e s i g n e da n dt h ec o m p e n s a t i o nm o d e lw a sd e r i v e d . T h e yf u r t h e rw e a k e n e dt h ea d v e r s e e f f e c t s o nt h es y s t e m d u et ot h ec h a n g e s i nt h ee x t e r n a l l o a d ,a n di m p r o v e dt h ep o s i t i o nc o n t r o l a c c u r a c y o f t h es y s t e m .T h ec o m p o s i t ec o n t r o l s t r a t e g yw a s v e r i f i e do nt h es i m u l a t i o ni nMa t l a ba n ds e m i p h y s i c a l s i m u l a t i o np l a t f o r m ,r e s p e c t i v e l y .A n dt h es i m u l a t i o na n de x p e r i m e n t a lr e s u l t ss h o w t h a tt h ee l e c t r o h y d r a u l i cp o s i t i o nc o n t r o l s y s t e mw i t hA D R Cc o n t r o l l e r c o u l de f f e c t i v e l yi n h i b i t t h ei n t e r f e r e n c e o f e x t e r n a l d i s t u r b a n c e s ,a n dt h ep r e c i s ep o s i t i o n i n gc o n t r o lw a sr e a l i z e da f t e ri n t r o d u c i n gt h el o a dc o m p e n s a t i o n c o n t r o l l e r .A n dt h e nt h ee f f i c i e n c yo f t h em e t h o dp r e s e n t e dw a s v e r i f i e dt h r o u g hs i m u l a t i o n s . Ke yw o r d s :E l e c t r o h y d r a u l i cp o s i t i o nc o n t r o l s y s t e m L o a dc o m p e n s a t i o n A D R Cc o n t r o l l e r ( 上接第 3 3 2页) L e a k a g eMe c h a n i s m o f S c r e wP u mpB a s e do n L e a k a g eMo d e l i nF l u i dMe c h a n i c s Z h a n gY u a n x u n 1 T a n gQ i a n 2 L i Z h o n g h u a 2 Y a nD i 2 ( 1 . C o l l e g eo f A e r o s p a c eE n g i n e e r i n g ,C h o n g q i n gU n i v e r s i t y ,C h o n g q i n g4 0 0 0 3 0 ,C h i n a 2 . S t a t eK e yL a b o r a t o r yo f Me c h a n i c a l T r a n s m i s s i o n ,C h o n g q i n gU n i v e r s i t y ,C h o n g q i n g4 0 0 0 3 0 ,C h i n a ) A b s t r a c t :C l e a r a n c el e a k a g ei s t h em a i nc a u s eo f p u m pl e a k a g e .A i m i n ga t t h ed e s i g na n dc o n t r o l o f t h e k e yc l e a r a n c e ,w h i c hd e t e r m i n e st h ep e r f o r m a n c eo f s c r e wp u m ps u c ha sc l e a r a n c eo f b a r r e l w a l l a n d m e s h i n ga r e a ,t h em o d e l o f d i f f e r e n t i a l p r e s s u r ef l o wa n ds h e a r f l o wo f c l e a r a n c el e a k a g ei np u m pc a v i t y i s e s t a b l i s h e d .T h em a t h e m a t i c a l e x p r e s s i o n s o f l e a k a g ef r o md i f f e r e n t c l e a r a n c ea r eo b t a i n e db ya n a l y z i n g t h ec l e a r a n c el e a k a g eo f b a r r e l w a l l a n dm e s h i n ga r e a .T h ei n f l u e n c eo f t h er o t a t i o ns p e e da n dp r e s s u r e d i f f e r e n c eb e t w e e ni n l e ta n do u t l e to nf l o w c h a r a c t e r i s t i c sa n dv o l u m e t r i ce f f i c i e n c yc h a r a c t e r i s t i c so f s c r e wp u m pi s r e v e a l e db yt h e o r e t i c a l c a l c u l a t i o na n de x p e r i m e n t a l a n a l y s i s o f t w i ns c r e wp u m p .B a s e do n t h e s er e s u l t s ,t h er o t o r d e s i g nt h e o r yo f t h es c r e wp u m pa s w e l l a s t h ep e r f o r m a n c eo f t h es c r e wp u m pc a n b ei m p r o v e d . Ke yw o r d s :S c r e wp u m p L e a k a g e C l e a r a n c e V o l u m e t r i ce f f i c i e n c y 9 3 3 第10期 高炳微 等:基于负载力补偿的自抗扰复合电液位置控制方法