TW201323902A - 感應機參數估測方法 - Google Patents
感應機參數估測方法 Download PDFInfo
- Publication number
- TW201323902A TW201323902A TW100145775A TW100145775A TW201323902A TW 201323902 A TW201323902 A TW 201323902A TW 100145775 A TW100145775 A TW 100145775A TW 100145775 A TW100145775 A TW 100145775A TW 201323902 A TW201323902 A TW 201323902A
- Authority
- TW
- Taiwan
- Prior art keywords
- parameter
- impedance
- induction machine
- speed
- steady
- Prior art date
Links
Landscapes
- Control Of Ac Motors In General (AREA)
Abstract
一種感應機參數估測方法,係包含:一轉速上升率降低步驟;一阻抗轉差率曲線求取步驟;一初值求取步驟;一阻抗參數最佳化步驟;及一機械參數計算步驟;該方法係於低壓無載啟動下,估測該感應機之參數。
Description
本發明係關於一種感應機參數估測方法,尤其是一種低壓無載啟動狀態之感應機參數估測方法。
感應電動機廣泛應用於工業領域,其應用層面已由單純的動力提供,逐漸延伸至精確控制(轉速、轉矩及定位),為了提高該感應電動機之控制精確度,該感應電動機之一相關參數是不可或缺的資訊。而該相關參數之取得方法,可利用該感應電動機之等效模型進行估測,其中,常用的等效模型可分成穩態等效模型與動態等效模型二種。
在穩態等效模型中,該相關參數中之阻抗參數可由標準測試取得,例如透過堵轉試驗、無載啟動試驗及定子側的直流測試,即可獲得該穩態等效模型之阻抗參數。在該穩態等效模型架構下,使用者雖可分析穩態時之運轉狀況,然而卻因未加入負載的變化與電源端信號變動的暫態考量,使得穩態等效模型在暫態特性上的描述稍嫌不足。
在動態等效模型中,係考慮電壓、電流、定子頻率及負載變動等暫態效應。常見的估測方式,有頻譜分析技術、參考模型適應技術及最佳化技術等。相關技術所獲得的動態等效模型雖能更完整的描述該暫態特性,然而在該相關參數之估測過程中,須使用變流器等額外之控制電路,藉以激勵出響應,造成估測所需之實施硬體增加。
此外,不論是穩態等效模型或動態等效模型,在分析估測的過程中,皆僅能取得該相關參數中之阻抗參數,無法得到該相關參數中之一機械參數(轉動慣量與摩擦係數)。有鑑於此,必須要有一套兼顧低硬體成本、高估測準度與完整參數求取之感應機參數估測方法,以期符合不同需要。
本發明之主要目的係提供一種感應機參數估測方法,該參數估測方法具有較低之硬體成本負擔與簡易的架構。
本發明之另一目的係提供一種感應機參數估測方法,該參數估測方法可估測感應機較準確之一相關參數。
本發明之另一目的係提供一種感應機參數估測方法,該參數估測方法可取得較完整之相關參數。
為達到前述發明目的,本發明之感應機參數估測方法,係包含:一轉速上升率降低步驟,係降低無載啟動之一感應機之輸入電壓,使該感應機啟動之轉速上升率降低;一阻抗轉差率曲線求取步驟,係利用一量測儀器量測並記錄該感應機由啟動至穩態轉速期間之一時變電壓、一時變電流及一時變轉速,並藉由一時域分析求得一阻抗轉差率曲線;一初值求取步驟,係利用定子側之直流測試搭配阻抗轉差率曲線,以最小平方法來計算阻抗參數之初值;一阻抗參數最佳化步驟,係藉由阻抗參數的初值為基礎,以一最小均方根法結合一粒子群演算法進行該阻抗參數之最佳化處理;及一機械參數計算步驟,係將該阻抗參數代入穩態等效電路中,並搭配由啟動至穩態轉速期間的時變電壓、轉速信號與機械方程式求得一機械參數。
本發明之感應機參數估測方法,其中該輸入電壓係降低至待測之該感應機的最低啟動電壓。
本發明之感應機參數估測方法,其中該最小均方根法係用以定義一目標函數,且該目標函數具有數梯度量。
本發明之感應機參數估測方法,其中該數梯度量係作為該粒子群演算法之一修正量的修正分量。
為讓本發明之上述及其他目的、特徵及優點能更明顯易懂,下文特舉本發明之較佳實施例,並配合所附圖式,作詳細說明如下:本發明所述之「梯度量」,係指目標函數的一次微分。
本發明所述之「無載」,係指該感應機之動力輸出端無連接任何傳動構件,為無負載狀態之運轉。
本發明所述之「轉差率」,係指一轉差值與定子旋轉磁場之一同步轉速的比值,其中,該轉差值為該同步轉速和實際轉速之差值。
請參照第1圖所示,本發明之感應機參數估測方法係包含:一轉速上升率降低步驟S1、一阻抗轉差率曲線求取步驟S2、一初值求取步驟S3、一阻抗參數最佳化步驟S4及一機械參數計算步驟S5,並以一量測儀器記錄一感應機於低壓無載啟動時,由啟動至穩態轉速期間的時變電壓、時變電流與時變轉速,再透過一電腦等運算處理器進行該感應機之相關參數估測。
該轉速上升率降低步驟S1,係降低無載啟動之該感應機之輸入電壓,使該感應機啟動之轉速上升率降低。在該感應機正常啟動狀態下,該感應機之動態等效模型與穩態等效模型具有不同特性。更詳言之,在定子參考模型下,該動態等效模型可表示為:
其中i qs ,i ds 為定子電流、i qr ,i dr 為轉子電流、v qs ,v ds 為定子電壓、v qr ,v dr 為轉子電壓、R s 為定子電阻、R r 為轉子電阻、L m 為互感、L s 為定子電感、L r 為轉子電感、ω r 為轉子轉速及p為微分因子。
在方程式(1)中阻抗矩陣的部分項與轉速相依,為轉速的函數,在轉速變動的情況下,方程式組為非線性。若轉速為定值或是零,則方程式組為線性。當降低該感應機啟動之輸入電壓,使該感應機處於低壓無載啟動狀態時,可使該感應機之轉速上升率在單位時間內變化的幅度較小,即dω r /dt 0,使轉速變動所產生之暫態效應降至最低,便可忽略該暫態效應,即可利用穩態等效模型來進行該感應機之阻抗參數估測,且不需額外之控制電路。該動態等效模型所近似之該穩態等效模型可如第2圖所示,且該穩態等效電路可表示為轉速的函數,並具有一等效電阻與一等效電抗,該等效電阻與等效電抗分別為:
其中,該R s 為定子電阻、R r 為轉子等效電阻、X m 為磁激電抗、X s 為定子等效電抗、X r 為轉子等效電抗、s(n)為轉差率、ω s 為同步轉速及ω r (n)為轉速,其中n=0,1,2…N-1。
在本實施例中,該輸入電壓欲降低之底限在此並不設限,可根據各個不同規格之感應機進行調整,較佳可將該輸入電壓降低至待測之該感應機的最低啟動電壓,使該感應機之轉速可於啟動時緩慢上升。
該阻抗轉差率曲線求取步驟S2,係利用該量測儀器量測並記錄該感應機由啟動至穩態轉速期間之一時變電壓、一時變電流及一時變轉速,並藉由一時域分析求得一阻抗轉差率曲線。
當該感應機已藉由降低啟動電壓,致使該感應機之轉速可於啟動時緩慢上升後,該感應機之阻抗參數便可利用穩態等效模型描述,並具有該等效電阻與該等效電抗。為準確求得該等效電阻與等效電抗與轉差率之關係,可利用該量測儀器記錄該感應機由啟動至穩態轉速期間之時變電壓、時變電流及時變轉速,且該時變電壓、時變電流及時變轉速具有相同之取樣刻度。其中,該時變電壓與時變電流可經由時域分析解析出每一取樣刻度之振幅與相位,並藉此獲得各取樣點相對應之一次側等效電阻與等效電抗;該轉差率可藉由各個取樣點之該轉速求得,並藉此獲得各取樣點相對應之轉差率。在具有相同取樣刻度之條件下,可得到該阻抗轉差率曲線,在本實施例中,該阻抗轉差率曲線如第3圖所示。
該初值求取步驟S3,係利用定子側之直流測試搭配阻抗轉差率曲線,以最小平方法來計算阻抗參數之初值。阻抗參數的最佳化必須以一組初值作為搜尋的出發點,而適當合理的初值將有助於提升收斂效率。在本實施例中,係於該感應機操作於一額定電流時,利用該量測儀器量測該感應機之操作電壓與操作電流,並求得該定子電阻,若定子繞組為Y接型式,則該方程式為:
其中,V DC 為該操作電壓,I DC 為該操作電流。
當該定子電阻為已知時,可運用該操作電壓、操作電流及轉速計算得該s(n)與該R(s(n)),透過將該定子電阻代入方程式(2)整理可得:
其中
在最小平方法中,令一誤差平方和為:
其中最小化的必要條件為:
經整理可得:
當轉差率s(n)與Q(s(n))為已知的情況,代入方程式(8)即可計算得到K 1 與K 2 。在該阻抗參數中,由於磁激電抗X m >>X s 且X m >>X r ,故在方程式(3)中假設X r =X s ,並代入方程式(3):
其中,
同理,令一誤差平方和為:
其中,K(n)=(1+(s(n))2 K 1),同理可得:
在此定義K 5=K 4/K 3=,則其他阻抗參數之初值為:
X s =X r =K 3-X m (14)
透過上述計算過程可決定該阻抗參數之初值,該阻抗參數之初值可作為後續之阻抗參數最佳化處理之基礎,以取得較佳之該阻抗參數。
該阻抗參數最佳化步驟S4,係藉由阻抗參數的初值為基礎,以一最小均方根法(Least Mean Square,LMS)結合一粒子群演算法(Particle Swarm Optimization,PSO)進行該阻抗參數之最佳化處理,以搜尋該阻抗參數的廣域最佳解。首先詳述最小均方根法的部份,最小均方根法可處理多變數非線性函數的最小化搜尋。將步驟S2所得之阻抗轉差率曲線表示為R(s(n))、X(s(n))及s(n),n=0,1,2…N-1。在LMS可將其目標函數ER與EX表示為:
其中,
M(n)=(R r /s(n))2+(X m +X r )2 (20)
該目標函數ER與EX皆為具有一極小值之函數,在該極小值時,該二函數對每個阻抗參數的梯度量均為0。本發明以R s 、R r 搜尋目標函數ER的最小值,以X s 、X r 及X m 搜尋目標函數EX的最小值,該電阻的誤差量可表示為:
因此該目標函數ER與目標函數EX最小值的數梯度量可表示為:
最小均方根法在本發明中可強化局部搜尋能力,再配合粒子群演算法,即可完成廣域搜尋的最佳化。在此為達到整體最佳化的目的,其適應函數可定義表示為:
FIT=ER+EX (28)
因本發明於步驟S3中提供一組合理的初值,在粒子群演算法中的族群觀念下,僅需將粒子以亂數形式散佈於初值附近的參數空間上,以便針對適當區域進行廣域搜尋,每個粒子在在參數空間可表示為:
該粒子群演算法之粒子位置更新方式如下所示:
其中,代表更新後之粒子位置,代表更新前之粒子位置,代表粒子位置更新之一修正量。其中,該修正量具有梯度量,且該梯度量係作為該修正量之修正分量,該修正量如下所示:
其中,C 1 為粒子自身經驗之認知參數,C 2 為粒子群體經驗之認知參數,為粒子本身在k次疊代過程中之最佳位置,為族群中k次疊代過程中之最佳位置,▽E i 為該目標函數ER與EX之數梯度量,即為方程式(23)至(27),α k 為加速因子,iter max為設定之疊代次數,αmax為設定之最大加速因子,αmin為設定之最小加速因子。
在本實施例中,係設定該C 1 =C 2 =1,設定該αmax=0.5,設定該αmin=0.3,且當該粒子位置之更新次數到達該疊代次數時,便結束該阻抗參數最佳化步驟S4。更詳言之,在阻抗參數搜尋的過程中,由於粒子群演算法的特性,使參數最佳化的過程避開收斂性的問題。對粒子群演算法而言,其本身具有於廣域中搜尋最佳解之特性,但當其粒子身處最佳解附近時,因其本身鄰近搜尋能力較差,有時候可能會遺漏掉最佳解,故在該粒子群演算法中加入該梯度量,以該梯度量做為該粒子群演算法之修正量的修正分量,可使該最佳化演算法同時兼具廣域與鄰近搜尋的能力,以求得最佳之該阻抗參數。
該機械參數計算步驟S5,係將該阻抗參數代入穩態等效電路中,並搭配由啟動至穩態轉速期間的時變電壓、轉速信號與機械方程式求得一機械參數。請參照第2圖所示,更詳言之,將阻抗參數代入穩態等效電路中,並搭配由啟動至穩態轉速期間的時變電壓與轉速,可計算一時變的氣隙功率P AG ,該氣隙功率可表示為:
一機械功率P m 可表示為:
其中該P RCL 為轉子電路的等效銅損。
一電磁轉矩T e 可表示為:
在動態等效模型中,該感應機之該機械方程式可表示為:
其中ω m 為轉子的機械轉速,即ω m =2ω r /P。P為該感應機的極數,J為轉動慣量,B為摩擦係數,在無載啟動的條件下,令T Load =0,整理上式可得:
在無載啟動的條件下,當轉子轉速由靜止狀態逐漸加速至穩態轉速時,考慮在穩態的狀況下,轉速為固定,即dω m /dt=0。將該條件代入方程式(36),可得該摩擦係數B為:
當該摩擦係數為已知時,可藉由下列方程式求得該轉動慣量J:
其中,該轉動慣量J之平均值可由轉速自靜止狀態加速至穩態轉速之區間求得。
透過將該感應機操作於低壓無載啟動狀態,可在單位時間內,使該感應機之動態特性近似於穩態特性,令該阻抗參數之取得不需透過額外的控制電路。藉由量測該感應機啟動至穩態轉速期間之時變電壓、時變電流及時變轉速,可求出該感應機之阻抗轉差率曲線,並藉此以最小平方法計算出一組合理的阻抗參數的初值。接著以該初值為基礎,利用該粒子群演算法結合最小均方根法進行最佳解之搜索,且該粒子群演算法之修正量中,含有該梯度量為其修正量之修正分量,使該阻抗參數得以取得最佳解。且在低壓無載啟動狀態下,配合最佳化之該阻抗參數,可逐步獲得較準確之該機械參數。
本發明之感應機參數估測方法,係降低啟動電壓使該感應機處於低壓無載啟動狀態,以進行後續之相關參數估測,該方式不需設置複雜的控制電路,具有節省額外硬體成本之功效。
本發明之感應機參數估測方法,係利用一梯度量結合一粒子群演算法進行阻抗參數之最佳化搜索,具有最佳化阻抗參數之功效。
本發明之感應機參數估測方法,係降低啟動電壓使該感應機處於低壓無載啟動狀態,有助於取得一機械參數,具有取得較完整之相關參數功效。
雖然本發明已利用上述較佳實施例揭示,然其並非用以限定本發明,任何熟習此技藝者在不脫離本發明之精神和範圍之內,相對上述實施例進行各種更動與修改仍屬本發明所保護之技術範疇,因此本發明之保護範圍當視後附之申請專利範圍所界定者為準。
[本發明]
S1...轉速上升率降低步驟
S2...阻抗轉差率曲線求取步驟
S3...初值求取步驟
S4...阻抗參數最佳化步驟
S5...機械參數計算步驟
第1圖:本發明感應機參數估測方法之流程圖。
第2圖:本發明感應機參數估測方法之穩態等效電路圖。
第3圖:本發明感應機參數估測方法之阻抗轉差率曲線。
S1...轉速上升率降低步驟
S2...阻抗轉差率曲線求取步驟
S3...初值求取步驟
S4...阻抗參數最佳化步驟
S5...機械參數計算步驟
Claims (4)
- 一種感應機參數估測方法,係包含:一轉速上升率降低步驟,係降低無載啟動之一感應機之輸入電壓,使該感應機啟動之轉速上升率降低;一阻抗轉差率曲線求取步驟,係利用一量測儀器量測並記錄該感應機由啟動至穩態轉速期間之一時變電壓、一時變電流及一時變轉速,並藉由一時域分析求得一阻抗轉差率曲線;一初值求取步驟,係利用定子側之直流測試搭配阻抗轉差率曲線,以最小平方法計算阻抗參數之初值;一阻抗參數最佳化步驟,係藉由阻抗參數的初值為基礎該梯度量作為一最佳化演算法之修正分量,以一最小均方根法結合一粒子群演算法為基礎,利用該最佳化演算法進行該阻抗參數之最佳化處理;及一機械參數計算步驟,係將該阻抗參數代入穩態等效電路中,並搭配由啟動至穩態轉速期間的時變電壓、轉速信號與機械方程式求得一機械參數。
- 如申請專利範圍第1項所述之感應機參數估測方法,其中該輸入電壓係降低至待測之該感應機的最低啟動電壓。
- 如申請專利範圍第1項所述之感應機參數估測方法,其中該最小均方根法係用以定義一目標函數,且該目標函數具有數梯度量。
- 如申請專利範圍第3項所述之感應機參數估測方法,其中該數梯度量係作為該粒子群演算法之一修正量的修正分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW100145775A TWI431300B (zh) | 2011-12-12 | 2011-12-12 | 感應機參數估測方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW100145775A TWI431300B (zh) | 2011-12-12 | 2011-12-12 | 感應機參數估測方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201323902A true TW201323902A (zh) | 2013-06-16 |
TWI431300B TWI431300B (zh) | 2014-03-21 |
Family
ID=49032903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW100145775A TWI431300B (zh) | 2011-12-12 | 2011-12-12 | 感應機參數估測方法 |
Country Status (1)
Country | Link |
---|---|
TW (1) | TWI431300B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106452258A (zh) * | 2016-11-11 | 2017-02-22 | 福建睿能科技股份有限公司 | 一种三相感应电机参数检测方法及装置 |
CN115001353A (zh) * | 2022-05-09 | 2022-09-02 | 上海达坦能源科技股份有限公司 | 抽油机井变频调速智能控制方法及控制系统 |
-
2011
- 2011-12-12 TW TW100145775A patent/TWI431300B/zh not_active IP Right Cessation
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106452258A (zh) * | 2016-11-11 | 2017-02-22 | 福建睿能科技股份有限公司 | 一种三相感应电机参数检测方法及装置 |
CN106452258B (zh) * | 2016-11-11 | 2019-06-11 | 福建睿能科技股份有限公司 | 一种三相感应电机参数检测方法及装置 |
CN115001353A (zh) * | 2022-05-09 | 2022-09-02 | 上海达坦能源科技股份有限公司 | 抽油机井变频调速智能控制方法及控制系统 |
Also Published As
Publication number | Publication date |
---|---|
TWI431300B (zh) | 2014-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lin et al. | Parameter identification of induction machine with a starting no-load low-voltage test | |
Raca et al. | Robust magnet polarity estimation for initialization of PM synchronous machines with near-zero saliency | |
Feng et al. | Particle-filter-based magnet flux linkage estimation for PMSM magnet condition monitoring using harmonics in machine speed | |
Lai et al. | Torque ripple minimization for interior PMSM with consideration of magnetic saturation incorporating online parameter identification | |
CA2802448A1 (en) | System and method of speed detection in an ac induction machine | |
CN110535392B (zh) | 一种基于lm算法的永磁同步电机参数辨识方法 | |
JP2012519464A (ja) | 永久磁石励磁形同期機械の回転子シャフトの位置を、電流信号と電圧信号に基づいて角度センサなしで決定する方法および装置 | |
JP2016073200A (ja) | 永久磁石同期モーターにおける鉄損を決定するための方法及びシステム | |
Odhano et al. | Parameter identification and self-commissioning of AC permanent magnet machines-A review | |
Dlamini et al. | A non-intrusive method for estimating motor efficiency using vibration signature analysis | |
TW201710924A (zh) | 以多項式估測感應機參數之方法 | |
CN210111901U (zh) | 用于检测绕组温度的装置以及感应电机 | |
Grantham et al. | Rapid parameter determination for induction motor analysis and control | |
CN103718454B (zh) | 永磁同步电动机的电动机常数计算方法以及电动机常数计算装置 | |
CN106602953B (zh) | 基于磁场定向准确性的感应电机转子时间常数的验证方法 | |
TWI431300B (zh) | 感應機參數估測方法 | |
CN105699896B (zh) | 一种适用于感应电机转子故障诊断方法 | |
JP2008286779A (ja) | Ipmモータのためのトルク推定器 | |
Zhang et al. | Prediction and experiment of DC-bias iron loss in radial magnetic bearing for a small scale turbomolecular pump | |
KR101591198B1 (ko) | 영구자석 동기 전동기의 인덕턴스 추정기 및 영구자석 동기 전동기의 인덕턴스 추정방법, 그 방법을 수행하기 위한 프로그램이 기록된 기록매체 | |
CN103155400A (zh) | Pm电动机的电动机常数计算方法以及电动机常数计算装置 | |
CN109932648B (zh) | 一种同步电机q轴电感饱和特性试验测量方法 | |
WO2023109403A1 (zh) | 永磁同步电机的转子温度的估算方法、装置和计算机设备 | |
Babel et al. | Parametric sensitivity in the analysis and control of permanent magnet synchronous machines | |
CN111726050B (zh) | 确定永磁同步电机的铁损的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Annulment or lapse of patent due to non-payment of fees |