韓道權(quán)，韓 超，丁寶林，葉衛(wèi)東

（1.東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院，黑龍江大慶163318；2.大慶油田裝備制造集團(tuán)，黑龍江大慶163255） ①

基于Matlab的單螺桿泵線型打扣分析

韓道權(quán)1，韓 超1，丁寶林2，葉衛(wèi)東1

（1.東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院，黑龍江大慶163318；2.大慶油田裝備制造集團(tuán)，黑龍江大慶163255）①

根據(jù)單螺桿泵定子骨線及其等距曲線的形成原理，利用Matlab軟件對(duì)單螺桿泵定子骨線打扣現(xiàn)象進(jìn)行了深入剖析，并利用數(shù)值法繪出了單螺桿泵定子頭數(shù)、變幅系數(shù)和不打扣條件下的最大等距半徑之間的關(guān)系圖形。結(jié)果表明：當(dāng)定子頭數(shù)確定時(shí)，變幅系數(shù)與不打扣最大等距半徑的關(guān)系曲線為遞減規(guī)律；當(dāng)變幅系數(shù)確定時(shí)，定子頭數(shù)越小，允許的不打扣半徑越大。

螺桿泵；曲線；干涉；數(shù)學(xué)分析

1 單螺桿泵定子骨線方程

定子骨線的研究主要是在復(fù)平柱面坐標(biāo)系下，它與常規(guī)的空間直角坐標(biāo)系O－x－y－z的不同之處在于以復(fù)平面1－O－j取代實(shí)平面x－O－y，并取向下方向?yàn)閦軸正方向。復(fù)平柱面坐標(biāo)系中，O－1稱(chēng)為實(shí)軸，單位是1；O－j稱(chēng)為虛軸，單位是，如圖1所示。單螺桿泵定子骨線是將半徑為r的滾圓沿半徑為R的導(dǎo)圓做相對(duì)純滾動(dòng)時(shí)，滾圓平面上任意點(diǎn)A在導(dǎo)圓平面上形成的軌跡。3頭定子骨線如圖2所示，圓O1為導(dǎo)圓，圓O2為滾圓。定子骨線在復(fù)平柱面坐標(biāo)系下的方程［1－2］為

式中，n為轉(zhuǎn)子的頭數(shù)；K為變幅系數(shù)，即發(fā)生點(diǎn)A到滾圓圓心O2的距離與r的比值。

圖1 復(fù)平柱面坐標(biāo)系

圖2 3頭定子骨線

2 定子輪廓線打扣現(xiàn)象分析

單螺桿泵定子輪廓線是圓心在骨線上、半徑為ra的一系列圓外包絡(luò)（即等距曲線）得到的。因此，當(dāng)?shù)染鄨A上外包絡(luò)點(diǎn)出現(xiàn)回行時(shí)，就會(huì)產(chǎn)生打扣現(xiàn)象［3－4］，如圖3所示。

圖3 等距曲線打扣

點(diǎn)A在等距曲線上的等距點(diǎn)為C，點(diǎn)B的等距點(diǎn)為D，D點(diǎn)出現(xiàn)了回行，這時(shí)就產(chǎn)生了打扣現(xiàn)象，即打扣現(xiàn)象產(chǎn)生的根本原因就是骨線上任意2點(diǎn)相鄰的法線在等距半徑范圍內(nèi)有交點(diǎn)。這時(shí)生產(chǎn)出來(lái)的定子并不是嚴(yán)格的擺線等距曲線。如果打扣較小，唯一的解決辦法是通過(guò)橡膠的過(guò)盈來(lái)進(jìn)行補(bǔ)償，但是如果打扣較大，則不能用這種方法解決。

由式（1）可得到骨線在復(fù)平柱面坐標(biāo)系下的切線方程為［5－8］

由式（2）可得直角坐標(biāo)系下坐標(biāo)方程為

由式（3）可推出，骨線在直角坐標(biāo)系下任意點(diǎn)的切線方程為

由式（3）～（4）可得骨線任意點(diǎn)上的法線方程為

通過(guò)上述公式利用軟件繪制定子骨線及其法線圖形。由于定子骨線是軸對(duì)稱(chēng)圖形，因此，只研究其1個(gè)對(duì)稱(chēng)區(qū)間的圖形即可。當(dāng)K＝1時(shí)，定子骨線及其法線圖形如圖4所示。當(dāng)θ＝0時(shí)，法線是垂直的，它與其他法線的焦點(diǎn)隨著θ的減小而越來(lái)越接近骨線，這時(shí)無(wú)論等距半徑取多大，等距曲線都會(huì)產(chǎn)生打扣現(xiàn)象。這也是由于骨線在θ＝0處為尖點(diǎn)所致。

圖4 K＝1時(shí)3頭骨線及其法線

K＝0.8時(shí)，骨線及其法線圖形如圖5所示。可以看出，在θ＝0處，骨線為弧形，因此法線是水平的，而后法線斜率先遞增，后遞減。但法線在一定范圍內(nèi)彼此間仍有交點(diǎn)。因此，等距曲線的打扣現(xiàn)象是在等距半徑增大到一定范圍內(nèi)產(chǎn)生的。

圖5 K＝0.8時(shí)3頭骨線及其法線

K＝0.5時(shí)，骨線及其法線圖形如圖6所示?？梢钥闯觯讦龋?處法線仍是水平的，而后法線呈遞增規(guī)律。法線在骨線外側(cè)是發(fā)散的，沒(méi)有交點(diǎn)，這時(shí)無(wú)論等距半徑取多大，等距曲線都不會(huì)產(chǎn)生打扣現(xiàn)象。

圖6 K＝0.5時(shí)3頭骨線及其法線

以3頭單螺桿泵骨線為例，通過(guò)Matlab繪制定子骨線法線斜率圖形，進(jìn)一步直觀判斷定子打扣現(xiàn)象。圖7自下而上分別是K為0.1～0.9（步長(zhǎng)0.1）時(shí)法線斜率；K＝1時(shí)法線斜率如圖8所示。

圖7 K＝0.1～0.9時(shí)3頭骨線法線斜率

圖8 K＝1時(shí)3頭骨線法線斜率

如圖7所示，當(dāng)K趨近于0時(shí)，斜率曲線向右下方下凹，并且單調(diào)遞增。隨著K的增加，曲線下凹趨勢(shì)不斷變緩，直至上凸。如圖8所示，當(dāng)K＝1時(shí)，上凸曲線左側(cè)消失，只剩下右側(cè)，變成唯一一條向左下方下凹的曲線，此時(shí)，在θ＝0處法線是垂直的。當(dāng)K＜1時(shí)，骨線在θ＝0處為弧形，因此法線為水平方向。當(dāng)法線斜率單調(diào)遞增時(shí)，彼此之間就不會(huì)有交點(diǎn)；當(dāng)法線斜率先遞增后遞減時(shí)，打扣現(xiàn)象是在一定范圍內(nèi)產(chǎn)生的；當(dāng)法線遞減（K＝1）時(shí)，定子輪廓線永遠(yuǎn)打扣。

3 數(shù)值法求解不打扣最大等距半徑

通過(guò)上述分析可知，定子輪廓線不打扣的最大等距半徑即為骨線外側(cè)法線交點(diǎn)到骨線距離中的最小值。基于這個(gè)條件，用Matlab進(jìn)行編程，求出這個(gè)最小值，過(guò)程如下：

通過(guò)上述程序可以找到在n、K條件下定子不打扣的最大等距半徑，并得到1組頭數(shù)不變、K變化的數(shù)據(jù)。通過(guò)曲線擬合的形式可以做出不同頭數(shù)、變幅系數(shù)與最大等距半徑之間的關(guān)系圖形，如圖9所示。圖9中，取ra＝1，即單位擺線；當(dāng)ra的值不為1時(shí)，最大等距半徑數(shù)值上等于單位擺線的等距半徑乘以滾圓半徑。

圖9 3～6頭骨線最大等距半徑曲線

由圖9可以看出，定子不打扣的最大等距半徑是隨著變幅系數(shù)的增加而逐漸減小的；圖中清晰地反應(yīng)出了定子頭數(shù)、最大不打扣等距半徑和變幅系數(shù)K三者之間的關(guān)系。例如，當(dāng)K＝0.85時(shí)，從圖中可得3～6頭定子不打扣的最大等距半徑分別為9.4、5.9、4.7和4.2mm，即頭數(shù)越少，允許的不打扣等距半徑越大。

普通內(nèi)擺線是當(dāng)式（1）中K＝1時(shí)的擺線，共軛曲線副是定轉(zhuǎn)子頭數(shù)相差1得到的。短副內(nèi)擺線型單螺桿泵是短副內(nèi)擺線做定子骨線，用被定子骨線內(nèi)包絡(luò)的共軛曲線作為轉(zhuǎn)子骨線。因此，如果定子是不打扣的，那么用這種包絡(luò)曲線的方法得到的轉(zhuǎn)子骨線及其等距曲線也是不打扣的［5］。因此，只要保證基礎(chǔ)骨線R（θ）的等距曲線不打扣，就可以得到嚴(yán)格的定轉(zhuǎn)子共軛曲線。基礎(chǔ)骨線的等距曲線打扣條件遵從圖9的規(guī)律。

4 結(jié)論

1） 當(dāng)變幅系數(shù)K＝1時(shí)，等距曲線永遠(yuǎn)打扣；當(dāng)0＜K＜1時(shí)，等距曲線永遠(yuǎn)不打扣。

2） 當(dāng)變幅系數(shù)K趨近于1時(shí)，不打扣等距半徑迅速下降，越接近1，可用等距半徑越小。

3） 不同頭數(shù)、相同變幅系數(shù)、頭數(shù)越少時(shí)，允許的不打扣等距半徑越大。

4） 相同頭數(shù)、變幅系數(shù)越小時(shí)，允許的不打扣等距半徑越大。

［1］ 蘇義腦，謝竹莊.螺桿鉆具馬達(dá)線型分析基礎(chǔ)及研究方法［J］.石油機(jī)械，1985，13（6）：10－20.

［2］ 杜秀華，任 彬，韓國(guó)有.雙頭單螺桿泵的線型設(shè)計(jì)及虛擬建模［J］.石油礦場(chǎng)機(jī)械，2007，36（1）：33－35.

［3］ 萬(wàn)邦烈.單螺桿泵式水利機(jī)械的研究與開(kāi)發(fā)［J］.石油礦場(chǎng)機(jī)械，1995，24（3）：14－18.

［4］ 熊建新，陳次昌，魏存祥.雙頭單螺桿泵襯套線型分析［J］.石油機(jī)械，2008，36（12）：31－33.

［5］ 宋玉杰，溫后珍，孟碧霞.單螺桿式水力機(jī)械變幅系數(shù)的選擇［J］.中國(guó)機(jī)械工程，2010，21（18）：69－71.

［6］ 王國(guó)慶，師國(guó)臣，馬志權(quán)，等.大慶油田螺桿泵機(jī)采井系統(tǒng)效率現(xiàn)狀及對(duì)策［J］.石油礦場(chǎng)機(jī)械，2011，40（7）：25－28.

［7］ 韓國(guó)有，姜長(zhǎng)鑫，杜秀華.普通內(nèi)擺線型雙頭螺桿泵建模及運(yùn)動(dòng)仿真［J］.石油礦場(chǎng)機(jī)械，2011，40（1）：1－4.

［8］ 王春陽(yáng)，夏秀文，李良軍，等.預(yù)輪廓定子螺桿鉆具性能研究［J］.石油礦場(chǎng)機(jī)械，2011，40（5）：90－92.

Analysis of Linear Dozen Buckle of Single Screw Pump Based on Matlab

HAN Dao－quan1，HAN Chao1，DING Bao－lin2，YE Wei－dong1
（1.College of Mechanical Science and Engineering，Northeast Petroleum University，Daqing163318，China；2.Daqing Petroleum Equipment Group，Daqing163255，China）

According to formation mechanism of stator skeleton lines and its isometric profiles of screw pump，the phenomenon of Dozen buckle of bone line of single screw pump stator has been analyzed by mathematical software Matlab，the relationship between the different number of single screw pump stator，the radius variation ratio and the radius of the maximum isometric under the condition of no dozen buckle has been painted out.The results show that：when the stator head is determined，the radius variation ratio and the radius of the maximum isometric show a decreasing trend；when the radius variation ratio is determined，the fewer the number of stator head，the larger the radius.

screw pump；curve；interaction；mathematic analysis

1001－3482（2011）12－0038－04

TE933.3

A

2011－06－13

黑龍江省教育廳科技攻關(guān)項(xiàng)目（11551014）

韓道權(quán)（1966－），男，黑龍江大慶人，副教授，主要從事機(jī)械設(shè)計(jì)制造及其自動(dòng)化研究，E－mail：handq2003＠163.com。