矢量變頻調(diào)速驅(qū)動下風機軸系脈動扭矩分析

江哲帆, 王青華, 張禮亮, 楊建剛

(1.東南大學 火電機組振動國家工程研究中心，江蘇 南京 210096 2.上海市明華電力技術(shù)工程有限公司,上海 200082)

矢量變頻調(diào)速驅(qū)動下風機軸系脈動扭矩分析

江哲帆1, 王青華2, 張禮亮1, 楊建剛1

(1.東南大學 火電機組振動國家工程研究中心，江蘇 南京 210096 2.上海市明華電力技術(shù)工程有限公司,上海 200082)

針對某型變頻運行風機上發(fā)生的軸系故障現(xiàn)象，在MATLAB/Simulink環(huán)境下建立了矢量變頻調(diào)速方式下機電耦合模型，分析了升速過程中轉(zhuǎn)速和傳遞扭矩脈動現(xiàn)象。研究表明，風機升速到一定轉(zhuǎn)速時，電磁轉(zhuǎn)矩中出現(xiàn)了與軸系扭振固有頻率相近的諧波，使得反饋到電氣系統(tǒng)的速度信號中包含部分諧振分量，該分量被控制器放大后使電磁轉(zhuǎn)矩中的諧振分量增大，從而加劇了機械系統(tǒng)諧振，導致軸系傳遞扭矩出現(xiàn)大幅脈動。在速度反饋回路中加入以軸系扭振固有頻率為中心的陷波器可以取得較好的扭振抑制效果。

變頻調(diào)速； 矢量控制； 風機； 扭振

0 引 言

在火電廠的各類輔機中，絕大部分是風機和水泵類的機械設(shè)備。送風機、引風機、一次風機、煙道風機，這四大風機總耗電量約占機組發(fā)電量的2%[1]。因此，對火電廠風機系統(tǒng)進行節(jié)能改造具有極大經(jīng)濟價值。目前國內(nèi)許多火電廠對風機系統(tǒng)進行變頻調(diào)速改造以提高經(jīng)濟效益。

在進行變頻改造之后部分火電廠出現(xiàn)風機軸系扭振增大、不穩(wěn)定導致聯(lián)軸器失效和軸段裂紋的問題，嚴重影響設(shè)備安全性。文獻[2]介紹了某電廠5 A引風機進行變頻改造后在運行中出現(xiàn)斷軸故障的案例，研究人員計算了軸系扭振固有頻率，對大軸斷裂部分進行了應(yīng)力分析，并研究了大軸的疲勞受損情況，研究結(jié)果表明本次故障主要是由扭轉(zhuǎn)共振導致的。文獻[3]研究了傳遞軸上的脈動扭矩和電磁轉(zhuǎn)矩信號中與軸系扭振固有頻率相近的諧波的相關(guān)性，并從材料力學角度進行了壽命分析。文獻[4-5]從三角載波頻率和調(diào)制波頻率的組合出發(fā)，預測了逆變過程中可能產(chǎn)生的諧波頻率。本文在文獻[1-5]等研究的基礎(chǔ)上，在MATLAB/Simulink環(huán)境下建立了風機升轉(zhuǎn)速過程中電機側(cè)和風機側(cè)轉(zhuǎn)速脈動，以及軸系傳遞扭矩脈動的機電耦合計算分析模型，從諧波角度對扭振現(xiàn)象進行了分析，并研究了基于陷波器的扭振抑制方法的可行性。

1 分析模型

1. 1 電氣模型

分析模型包括電氣和機械兩部分，電氣模型又包括矢量控制系統(tǒng)和異步電機。

矢量控制系統(tǒng)通過坐標變化，將定子電流的磁場和轉(zhuǎn)矩分量解耦，分別加以控制?？刂葡到y(tǒng)中包含一個速度環(huán)和兩個電流環(huán)，速度環(huán)中給定轉(zhuǎn)速與反饋轉(zhuǎn)速經(jīng)PI調(diào)節(jié)形成閉環(huán)控制，兩個電流環(huán)分別控制磁場和電磁轉(zhuǎn)矩的形成。控制系統(tǒng)輸出的脈沖寬度調(diào)制頻率信號通過逆變器形成頻率可調(diào)的交流電。

通過坐標變換，將異步電機a、b、c三相坐標系上的各物理量變換到兩相旋轉(zhuǎn)dq坐標系中，各物理量之間有如下關(guān)系[6]。

電壓平衡方程：

(1)

式中：u——電壓；r——電阻；λ——磁通； 下標q、d——q軸和d軸上的物理量分量； 下標s、r——定子和轉(zhuǎn)子上的物理量分量；w、wr——兩相旋轉(zhuǎn)d、q坐標系旋轉(zhuǎn)速度和電動機轉(zhuǎn)子旋轉(zhuǎn)速度。

磁通平衡方程：

(2)

式中：L——電抗； 下標m、l——等效繞組之間的互感、漏電感。

電磁轉(zhuǎn)矩計算方程：

(3)

式中：Te——電磁轉(zhuǎn)矩；p——電機磁極對數(shù)。

當轉(zhuǎn)子磁鏈始終與d軸保持同向且為數(shù)值時，式(3)可以通過變換簡化為

(4)

式中：C——與電機參數(shù)有關(guān)的常數(shù)。

電磁轉(zhuǎn)矩完全由q軸上的定子電流分量控制，與d軸上的定子電流分量無關(guān)[7]。

1. 2 機械模型

機械模型包括電機、風機、傳遞軸、聯(lián)軸器、軸承等部分。本文將機械模型簡化為雙質(zhì)量塊模型[8]，將電機和風機兩個質(zhì)量塊通過傳遞軸相連接，如圖1所示。使用矩陣方式描述機械系統(tǒng)。

圖1 雙質(zhì)量塊模型

雙質(zhì)量塊模型的轉(zhuǎn)動慣量矩陣可以表示為

(5)

式中：Jm、Jl——電機和風機的轉(zhuǎn)動慣量。

抗扭剛度矩陣可以表示為

(6)

式中：kc——傳遞軸抗扭剛度。

傳遞軸上阻尼矩陣可以表示為

(7)

式中：cc——傳遞軸上的阻尼。

軸承阻尼矩陣可以表示為

(8)

式中：cm、cl——集中到電機和風機質(zhì)量塊上的軸承阻尼。

轉(zhuǎn)矩矩陣可以表示為

(9)

式中：Te、TL——電磁轉(zhuǎn)矩和負載轉(zhuǎn)矩。

(10)

完整的機械模型可以用矩陣形式表示為

(11)

1. 3 機電耦合模型

利用狀態(tài)空間法建立機電耦合模型，狀態(tài)空間方程可以表示為

(12)

在狀態(tài)空間中Te可表示為

(13)

機電耦合模型將電磁轉(zhuǎn)矩作為機械系統(tǒng)輸入，并采用轉(zhuǎn)速反饋方式。電氣模型和機械模型部分如圖2中所示。

1. 4 加入陷波器的機電耦合模型

電氣系統(tǒng)輸出的電磁轉(zhuǎn)矩信號帶有頻率豐富的諧波分量，其中頻率與機械系統(tǒng)扭振固有頻率重合的諧波分量將引起機械系統(tǒng)諧振。機械系統(tǒng)諧振導致反饋到電氣系統(tǒng)的速度信號帶有諧振成分。該成分經(jīng)過電氣系統(tǒng)控制器放大后使得電磁轉(zhuǎn)矩中諧振分量進一步增大，從而加劇了機械系統(tǒng)諧振。這樣的惡性循環(huán)使得軸系傳遞扭矩脈動過大，最終導致軸系損壞。

圖2 機電耦合模型

由上述分析可得，如果能在速度反饋環(huán)節(jié)中消除諧振成分，使得電磁轉(zhuǎn)矩中的諧振分量不再增大，將起到抑制扭振的作用。在測試轉(zhuǎn)速的場合中，通常將轉(zhuǎn)速信號轉(zhuǎn)換成電信號進行測量[9]，對電信號可以通過陷波器消除特定頻率的諧波。因此，考慮在速度反饋回路中加入陷波器。陷波器是一種特殊的濾波器，其在對信號中特定頻率的信號起到較好的抑制減小作用的同時，做到不影響其他頻率的信號。速度反饋回路中加入的陷波器如圖3所示。回路中以17.2 Hz為中心的陷波器特性如圖4所示。

圖3 串聯(lián)陷波器的機電耦合模型

2 計算分析實例

2. 1 系統(tǒng)參數(shù)

以某電機-風機系統(tǒng)參數(shù)為例進行分析，系統(tǒng)主要參數(shù)如表1所示。計算步長為1×10-5s。

表1 系統(tǒng)主要參數(shù)

2. 2 機電耦合模型扭振現(xiàn)象

圖5～圖7顯示了風機升速過程中電機側(cè)和風機側(cè)轉(zhuǎn)速以及軸系傳遞扭矩脈動的情況。

圖5 升速過程中電機側(cè)轉(zhuǎn)速

圖6 升速過程中風機側(cè)轉(zhuǎn)速

圖7 升速過程中軸系傳遞扭矩

計算分析結(jié)果表明，電機升速到510 r/min時，電機側(cè)和風機側(cè)轉(zhuǎn)速出現(xiàn)脈動，電機側(cè)轉(zhuǎn)速脈動幅度分別為13 rad/s，風機側(cè)轉(zhuǎn)速脈動幅度為1.6 rad/s。因為電機側(cè)和風機側(cè)轉(zhuǎn)速脈動幅度相差很大，軸系傳遞扭矩出現(xiàn)大幅脈動，脈動幅度高達2.3×106N·m。

2. 3 升速過程頻譜分析

對風機升速過程中電機側(cè)和風機側(cè)轉(zhuǎn)速、電氣系統(tǒng)輸出電磁轉(zhuǎn)矩和軸系傳遞扭矩的交流量進行頻譜分析。

2. 3. 1 轉(zhuǎn)速頻譜

圖8～圖9顯示了電機側(cè)和風機側(cè)轉(zhuǎn)速脈動后的頻譜分析。

電機側(cè)和風機側(cè)的轉(zhuǎn)速發(fā)生脈動前，信號中對應(yīng)軸系扭振固有頻率17.2 Hz的諧波分量很少，轉(zhuǎn)速發(fā)生脈動后，信號中17.2 Hz的諧波分量增多，電機側(cè)轉(zhuǎn)速信號中分量較大，幅值為47.8 r/min，風機側(cè)轉(zhuǎn)速信號中17.2 Hz分量較小，幅值為6.05 r/min，因此電機側(cè)轉(zhuǎn)速脈動幅度較風機側(cè)要大。

圖8 轉(zhuǎn)速脈動后電機側(cè)轉(zhuǎn)速頻譜

圖9 轉(zhuǎn)速脈動后風機側(cè)轉(zhuǎn)速頻譜

2. 3. 2 電磁轉(zhuǎn)矩頻譜

圖10為電機側(cè)和風機側(cè)轉(zhuǎn)速出現(xiàn)脈動后的電磁轉(zhuǎn)矩的頻譜分析。電機側(cè)和風機側(cè)轉(zhuǎn)速發(fā)生脈動前，電氣系統(tǒng)輸出的電磁轉(zhuǎn)矩信號中17.2 Hz的分量很少。轉(zhuǎn)速發(fā)生脈動后，電磁轉(zhuǎn)矩信號中出現(xiàn)大量頻率與軸系扭振固有頻率的各倍頻相近的諧波分量，其中，一至四倍頻諧波的分量幅值較大。電磁轉(zhuǎn)矩中17.2 Hz的諧波分量將導致機械系統(tǒng)諧振。

圖10 轉(zhuǎn)速脈動后電磁轉(zhuǎn)矩頻譜

2. 3. 3 軸系傳遞扭矩頻譜

圖11～圖12為電機側(cè)和風機側(cè)轉(zhuǎn)速出現(xiàn)脈動前后軸系傳遞扭矩的頻譜分析。轉(zhuǎn)速發(fā)生脈動前軸系傳遞扭矩信號中僅有少量頻率靠近軸系扭振固有頻率的諧波分量。轉(zhuǎn)速發(fā)生脈動后對應(yīng)軸系扭振固有頻率的17.2 Hz的諧波分量顯著增大，幅值為之前的數(shù)百倍。軸系傳遞扭矩中巨大的交流量將對機組安全運行產(chǎn)生影響。

圖11 轉(zhuǎn)速脈動前軸系傳遞扭矩頻譜

圖12 轉(zhuǎn)速脈動后軸系傳遞扭矩頻譜

2. 4 加入陷波器后扭振現(xiàn)象

圖13～圖15顯示了在速度反饋回路中加入陷波器后，電機升速過程中的轉(zhuǎn)速和傳遞扭矩。

圖13 加入陷波器后升速過程中電機側(cè)轉(zhuǎn)速

圖14 加入陷波器后升速過程中風機側(cè)轉(zhuǎn)速

圖15 加入陷波器后升速過程中軸系傳遞扭矩

結(jié)果表明，在陷波器投入使用后電機側(cè)轉(zhuǎn)速脈動幅度下降到0.84 rad/s，風機側(cè)轉(zhuǎn)速脈動基本消失，軸系傳遞扭矩脈動幅度為1.3×105N·m，僅為之前的5.7%，傳遞扭矩脈動現(xiàn)象得到抑制。

表2顯示了加入陷波器前后電機側(cè)和風機側(cè)轉(zhuǎn)速信號中對應(yīng)軸系扭振固有頻率17.2 Hz的諧波分量的絕對值，結(jié)果表明速度信號中該分量顯著減少，電機-風機系統(tǒng)的穩(wěn)定性得到改善。在速度反饋回路中增加陷波器來抑制扭振故障具有可行性。

表2 加入陷波器前后速度信號中17.2 Hz諧波分量

3 結(jié) 語

(1) 在電機-風機系統(tǒng)投變頻運行時，由于電機側(cè)和風機側(cè)轉(zhuǎn)速脈動幅度相差很大，導致軸系承受巨大的脈動扭矩，容易導致軸系斷裂等故障。

(2) 導致扭振故障的原因是電磁轉(zhuǎn)矩中頻率與軸系扭振固有頻率17.2 Hz相近的諧波分量引起了機械系統(tǒng)諧振，使反饋到電氣系統(tǒng)的速度信號帶有了諧振成分，經(jīng)過電氣系統(tǒng)控制器放大后電磁轉(zhuǎn)矩中諧振分量進一步增大，從而加劇了機械系統(tǒng)速度諧振,最終軸系傳遞扭矩脈動過大。

(3)在系統(tǒng)速度反饋回路中加入以軸系扭振固有頻率為陷波頻率的陷波器可以取得較好的扭振抑制效果。

[1] 周培建.火電廠引風機變頻節(jié)能改造的研究[D].青島: 青島理工大學,2014.

[2] 張磊,張俊杰,葛鴻昌,等.某型引風機變頻改造后電機斷軸原因分析[J].華北電力技術(shù),2013(12): 50-52.

[3] 石順梅,朱勝.增壓風機傳動軸斷裂失效分析[J].廣西電力,2012,35(4): 77-79.

[4] MANGUELLE J S, SCHRODER S, GEYER T, et al. Prediction of mechanical shaft failures due to pulsating torques of variable-frequency drives[J]. IEEE Transactions on Industry Applications, 2010, 46(5): 1979-1988.

[5] HAN X, PALAZZOLO A B. VFD machinery vibration fatigue life and multilevel inverter effect[J]. IEEE Transactions on Industry Applications, 2013, 49(6): 2562-2575.

[6] HAN X. Analysis and simulation of mechanical trains driven by variable frequency drivesystems[D]. Texas A&M University, 2010.

[7] 楊耕,羅應(yīng)立.電機與運動控制系統(tǒng)[M].北京:清華大學出版社,2014.

[8] 蔡昆,許鎮(zhèn)琳,張弢,等.電氣傳動中的扭振現(xiàn)象及其抑制[J].電氣自動化,2000(5): 19-21.

[9] 梁毓明,汪鋒鎖,劉志剛,等.JCZ型轉(zhuǎn)矩轉(zhuǎn)速傳感器輸出信號的處理方法[J].計算機測量與控制,2005,13(11): 1296-1298.

Analysis of Shaft Fluctuating Torque for Vector Variable Frequency Speed Regulation Driven Fan

JIANGZhefan1,WANGQinghua2,ZHANGLiliang1,YANGJiangang1

(1. National Engineering Research Center of Turbo Generator Vibration, Southeast University, Nanjing 210096, China 2. Shanghai Minghua Electric Power Technology Engineering Co., Ltd., Shanghai 200082, China)

On account of the shaft fault of a vector variable frequency speed regulation driven fan, mechanical-electromagnetic coupling model of the vector variable frequency speed regulation driven fan with the field-oriented control method was built in the MATLAB/Simulink. The phenomenon that the speed and transmission torque on the shaft fluctuate during the process of speed rise was analyzed. The research demonstrates harmonic with certain frequency close to the torsion vibration nature frequency of the shaft appears in the signal of electromagnetic torque after reaching a certain speed. Then the signal of speed which was the feedback of electrical system contains harmonic with the same frequency which aggravated the resonance after the being amplified by the controller. Finally the transmission torque on the shafting fluctuates wildly. The torsional vibration was well controlled with the notch filter whose trap frequency was the torsion vibration nature frequency in the speed feedback loop.

variable frequency speed regulation; vector control; fan; torsional vibration

江哲帆(1993—)，男，碩士研究生，研究方向為旋轉(zhuǎn)機械振動監(jiān)測與故障診斷。

TM 301.2

A

1673-6540(2017)04- 0101- 06

2016 -08 -31