蔡炯炯,曹 陽(yáng),程楷一,祝 亮,王子輝,蔣惠忠
(浙江科技學(xué)院 自動(dòng)化與電氣工程學(xué)院,浙江 杭州 310023)
其容積率高、精密可靠的優(yōu)點(diǎn),渦旋壓縮機(jī)因被廣泛應(yīng)用于制冷、航空航天與新能源等領(lǐng)域[1]。但其渦盤在高速運(yùn)行下,在平動(dòng)方向和軸向的高精度的動(dòng)態(tài)密封要求,對(duì)設(shè)計(jì)和加工帶來(lái)了極大挑戰(zhàn),阻礙了其往高壓力、大容量、長(zhǎng)壽命和無(wú)油方向的發(fā)展。
就軸向而言,其軸向氣體分離力隨主軸轉(zhuǎn)角和壓縮機(jī)出口氣體壓力的變化而大幅度變化,因此,要實(shí)現(xiàn)高性能軸向動(dòng)態(tài)密封[2],對(duì)渦旋壓縮機(jī)動(dòng)態(tài)的軸向氣體分離力進(jìn)行動(dòng)態(tài)精密平衡,是軸向密封的關(guān)鍵技術(shù)之一。目前,渦旋壓縮機(jī)軸向力的平衡方法中,有采用彈簧背壓式[3]、推力軸承式、氣體背壓式[4-5]等,這些方案在一定程度上能實(shí)現(xiàn)軸向力的平衡,但也存在難以實(shí)現(xiàn)高頻響高性能動(dòng)態(tài)調(diào)節(jié),或性價(jià)比不高的問(wèn)題。對(duì)此,筆者提出一種采用電磁力動(dòng)態(tài)平衡渦旋壓縮機(jī)軸向氣體分離力的方案,研究其動(dòng)態(tài)跟蹤控制方法。
渦旋壓縮機(jī)的軸向氣體分離力是周期性動(dòng)態(tài)變化的,平衡所需的電磁力在高動(dòng)態(tài)變化下其非線性特性,特別是在磁滯因素的影響下,會(huì)隨著主軸轉(zhuǎn)動(dòng)頻率的提高逐漸變得不可忽略,帶來(lái)明顯的滯后性。
對(duì)典型的大慣性、大滯后、非線性系統(tǒng)的跟蹤控制,已有較多的研究成果。動(dòng)態(tài)矩陣預(yù)測(cè)控制基于取測(cè)對(duì)象的階躍響應(yīng)為模型,采用反饋矯正滾動(dòng)實(shí)施優(yōu)化控制[6]。史密斯(Smith)預(yù)估控制,通過(guò)引入一個(gè)與被控對(duì)象并聯(lián)的補(bǔ)償器,對(duì)純滯后問(wèn)題進(jìn)行削弱與補(bǔ)償[7];通過(guò)對(duì)電磁系統(tǒng)遲滯特性研究,建立Preisach模型[8-9]、Prandtl-Ishlinskii(PI)模型[10]和Krasnoselskii-Pokrovskii(KP)模型[11]等唯象模型,和基于磁滯原理J-A模型[12]等,可較好提高系統(tǒng)動(dòng)態(tài)控制響應(yīng)能力。上述控制方法或者模型,均需要對(duì)系統(tǒng)預(yù)先進(jìn)行詳細(xì)復(fù)雜分析,確定系統(tǒng)模型結(jié)構(gòu),并以實(shí)驗(yàn)數(shù)據(jù)進(jìn)行參數(shù)識(shí)別。
在傳統(tǒng)PID控制中,本文采用引入提前跟蹤的思想,通過(guò)提前閉環(huán)跟蹤對(duì)沖對(duì)象和控制環(huán)節(jié)引入的滯后量;并在提前量的確定上,提出自適應(yīng)提前步數(shù)量尋找的方法,通過(guò)逐次迭代改進(jìn),找出最優(yōu)的步數(shù)提前量,并且在不同頻率時(shí)可以實(shí)現(xiàn)自行調(diào)整優(yōu)化。
渦旋壓縮機(jī)動(dòng)、靜渦盤之間的內(nèi)部壓縮腔在工作時(shí),會(huì)使動(dòng)、靜渦盤受到軸向氣體力的作用,該氣體力會(huì)使動(dòng)、靜渦盤分離,有使兩者之間的軸向間隙增大的趨勢(shì),會(huì)導(dǎo)致制冷劑泄漏增加,降低壓縮效率。
渦旋壓縮機(jī)渦旋盤所受的氣體軸向分離力是渦旋壓縮機(jī)固有機(jī)理所致,無(wú)法消除。筆者針對(duì)目前渦旋壓縮機(jī)動(dòng)靜渦盤間軸向力平衡或平衡方法性能上的不足,提出了一種電磁力平衡方案。
渦旋壓縮機(jī)軸向氣體力電磁平衡方案示意圖如圖1所示。
圖1 渦旋壓縮機(jī)軸向氣體力電磁平衡方案示意圖1-靜渦盤底板;2-吸力線圈鐵芯;3-吸力銅線圈;4-動(dòng)渦盤上被吸鐵芯;5-壓縮腔;6-靜渦盤渦卷;7-動(dòng)渦盤渦卷;8-壓力傳感器
渦旋壓縮機(jī)中,工質(zhì)的壓縮是通過(guò)動(dòng)、靜渦盤上的渦卷互相嚙合運(yùn)動(dòng)來(lái)實(shí)現(xiàn)的。多個(gè)壓縮腔同時(shí)進(jìn)行不同階段的壓縮循環(huán),壓縮腔內(nèi)部壓力對(duì)動(dòng)渦旋盤產(chǎn)生與角度有關(guān)的動(dòng)態(tài)變化的軸向氣體分離力:Fa=Fa1+Fa2+Fa3;在該分離力作用下,兩個(gè)渦盤之間有分離的趨勢(shì),通過(guò)采用安裝于動(dòng)靜渦盤上的環(huán)形電磁鐵,通過(guò)高性能控制可以快速高動(dòng)態(tài)平衡氣體分離力,又不至于過(guò)壓導(dǎo)致強(qiáng)烈摩擦,實(shí)現(xiàn)近零摩擦力的動(dòng)態(tài)接觸密封。
軸向氣體分離力與渦旋盤上渦卷的結(jié)構(gòu)參數(shù)和運(yùn)動(dòng)位置相關(guān),各個(gè)壓縮腔內(nèi)壓力對(duì)浮動(dòng)動(dòng)渦盤的推力之和即為動(dòng)渦盤所受軸向氣體力。由于壓縮腔面積與壓力隨著主軸轉(zhuǎn)角轉(zhuǎn)動(dòng)實(shí)時(shí)發(fā)生變化,軸向氣體力也會(huì)隨之發(fā)生變化。
一般的軸向氣體力可表示為:
Fa(θ)=
(1)
式中:ps—吸氣口壓力,MPa;P—漸開(kāi)線節(jié)距,mm;A1—中心壓縮腔軸向氣體力的作用面積,mm2;N—漸開(kāi)線的圈數(shù);ρi—第i個(gè)壓縮腔的壓力比;θ*—不同排氣壓力對(duì)應(yīng)轉(zhuǎn)角。
將模型樣機(jī)參數(shù)代入公式,可以計(jì)算得到軸向氣體力Fa隨主軸轉(zhuǎn)角θ的變化關(guān)系。此處取:N=3,ps=2 MPa,P=12 mm。
軸向力轉(zhuǎn)角關(guān)系圖如圖2所示。
圖2 軸向力-轉(zhuǎn)角關(guān)系圖
由圖2可以看出:軸向氣體分離力隨著主軸轉(zhuǎn)角呈現(xiàn)周期的變化規(guī)律,該氣體力在一個(gè)周期內(nèi)變化較大,用恒力的彈簧平衡效果很差,采用氣體背壓在高速運(yùn)轉(zhuǎn)時(shí)動(dòng)態(tài)性能極難匹配。
使用環(huán)形電磁鐵產(chǎn)生合適的動(dòng)態(tài)電磁力軌跡,用于平衡如圖2所示的氣體分離力時(shí),除了要設(shè)計(jì)合適的拓?fù)浜徒Y(jié)構(gòu)尺寸外,還需要合適的控制策略,特別需要解決在動(dòng)態(tài)跟蹤過(guò)程中,因滯后而帶來(lái)的控制難點(diǎn)。滯后的原因主要有以下二點(diǎn):
(1)鐵芯磁滯帶來(lái)的滯后。
采用圖1所示的環(huán)形電磁鐵可以產(chǎn)生高動(dòng)態(tài)性能電磁力,可用于平衡渦旋機(jī)中氣體分離力。但是電磁系統(tǒng)中控制電流與電磁力之間存在相位上的滯后。
電流-電磁力動(dòng)態(tài)關(guān)系曲線圖如圖3所示。
圖3 電流-電磁力動(dòng)態(tài)關(guān)系曲線圖
這種電磁力與勵(lì)磁電流、氣隙之間非線性關(guān)系,并非是遵循著簡(jiǎn)化的電磁力計(jì)算公式給出的理想模式,而是需要采用基于磁化原理的物理模型來(lái)表述,需要測(cè)定較多參數(shù)。不采用模型進(jìn)行動(dòng)態(tài)控制,難以對(duì)其進(jìn)行有效的跟蹤。
(2)離散控制系統(tǒng)導(dǎo)致的滯后。
軸向力控制系統(tǒng)采用力反饋控制,對(duì)于電磁壓力的采樣是周期采樣,即該系統(tǒng)具有典型的離散系統(tǒng)特點(diǎn),系統(tǒng)的運(yùn)行軌跡總是會(huì)產(chǎn)生多個(gè)采樣步數(shù)的滯后誤差,系統(tǒng)會(huì)出現(xiàn)一個(gè)長(zhǎng)度為h的時(shí)滯。針對(duì)離散系統(tǒng)不同的采樣頻率,時(shí)域上的時(shí)滯h在離散系統(tǒng)中的表達(dá)步數(shù)N也是不盡相同的。
基于上述原因,帶鐵芯電磁力動(dòng)態(tài)系統(tǒng)具有不可忽視的滯后特征。若采用在常規(guī)閉環(huán)中引入誤差微分項(xiàng)來(lái)解決,雖然微分項(xiàng)能夠在誤差到來(lái)之前,針對(duì)即將到來(lái)的誤差發(fā)出一個(gè)提前的修正信號(hào),但是微分存在著將高頻的干擾信號(hào)同時(shí)放大的缺點(diǎn);若在渦旋壓縮機(jī)軸向力平衡控制器中大比重使用微分,會(huì)造成渦旋壓縮機(jī)的抖動(dòng),嚴(yán)重影響渦旋壓縮機(jī)軸向力平衡的效果。
采用提前補(bǔ)償算法提高系統(tǒng)快速響應(yīng)性能時(shí),通過(guò)曲線提前跟蹤,對(duì)控制環(huán)節(jié)和對(duì)象本身引起的系統(tǒng)電壓-電流、電流-力之間的滯后補(bǔ)償對(duì)沖,將針對(duì)滯后進(jìn)行補(bǔ)償?shù)沫h(huán)節(jié)轉(zhuǎn)移到了閉環(huán)控制回路之外,因而,其不會(huì)對(duì)系統(tǒng)穩(wěn)定性產(chǎn)生不利影響。
通過(guò)對(duì)提前動(dòng)態(tài)跟蹤PID方法時(shí)域模型進(jìn)行分析,筆者研究了提前跟蹤效果的評(píng)價(jià)指標(biāo),并由此提出了自適應(yīng)提前步數(shù)尋優(yōu)的跟蹤算法。
基于傳統(tǒng)伺服控制系統(tǒng)中常規(guī)的PID能夠?qū)崿F(xiàn)基本的定位精度,但是由于滯后問(wèn)題的存在,輸入與輸出之間存在延時(shí)。若在PID反饋控制的基礎(chǔ)上,增加提前跟蹤環(huán)節(jié),在保證系統(tǒng)穩(wěn)定性的前提下,可對(duì)沖控制環(huán)節(jié)和對(duì)象引起的大慣性,提高系統(tǒng)動(dòng)態(tài)性能。
系統(tǒng)的控制器框圖如圖4所示。
圖4 提前跟蹤PID方法控制框圖
圖4中:在電磁鐵電磁吸力動(dòng)態(tài)控制中,PID控制環(huán)節(jié)為G1(s),驅(qū)動(dòng)電路可以簡(jiǎn)化等效成慣性環(huán)節(jié)G2(s),控制對(duì)象的電流與電磁力之間的關(guān)系等效成帶延遲特性的G3(s)。
可將其與系統(tǒng)的控制目標(biāo)R(s),反饋量C(s)組成一個(gè)閉環(huán)系統(tǒng),并在系統(tǒng)輸入中增加提前跟蹤環(huán)節(jié),來(lái)進(jìn)行分析。
圖4所示的閉環(huán)傳遞函數(shù)為:
(2)
若不加步數(shù)提前函數(shù),則系統(tǒng)的閉環(huán)傳遞函數(shù)為:
(3)
式中:G1(s)—PID調(diào)節(jié)環(huán)節(jié);G2(s)—直流斬波電路傳遞函數(shù);G3(s)—電磁吸力傳遞函數(shù)。
比較H(s)和H1(s),加入了超前環(huán)節(jié)之后系統(tǒng)的閉環(huán)特征方程根為:
D(s)=1+G1(s)G2(s)G3(s)=0
(4)
系統(tǒng)的閉環(huán)特征方程沒(méi)有改變,對(duì)于系統(tǒng)穩(wěn)定性能不會(huì)有影響,即:
(5)
當(dāng)H(s)=1時(shí),系統(tǒng)輸出完全跟隨輸入信號(hào),伺服將到達(dá)理想狀態(tài)。為使系統(tǒng)達(dá)到期望性能,需找到系統(tǒng)提前跟蹤時(shí)間τ最佳值。在嵌入式硬件執(zhí)行的離散控制中,具體為找到提前跟蹤步數(shù)N的最優(yōu)值。
對(duì)離散系統(tǒng)進(jìn)行分析設(shè)計(jì)時(shí),以一階的預(yù)測(cè)方法通過(guò)迭代的方式進(jìn)行。在用離散系統(tǒng)實(shí)現(xiàn)如2.1所示的控制策略時(shí),提前跟蹤是通過(guò)對(duì)被跟蹤目標(biāo)的信號(hào)進(jìn)行步數(shù)提前來(lái)實(shí)現(xiàn)的,即采用K+N時(shí)刻的目標(biāo)位置所需的力來(lái)決策K時(shí)刻的控制量。其中,N的值取系統(tǒng)滯后的步數(shù)τ值與系統(tǒng)采樣周期的比例,通過(guò)系統(tǒng)不斷迭代優(yōu)化,尋找系統(tǒng)能夠達(dá)到最佳跟蹤性能的最優(yōu)解,以實(shí)現(xiàn)系統(tǒng)穩(wěn)定運(yùn)行。
例如,將渦旋壓縮機(jī)工作時(shí)軸向力的每個(gè)周期按時(shí)間離散成360個(gè)點(diǎn),即式(1)中渦旋壓縮機(jī)軸向力可表示為Fα(K),K=0,1,2,3,…,359;則提前跟蹤PID方法時(shí)域的誤差的值取為:
e(K)=Fα(K+N)-F(K)
(6)
其PID控制環(huán)節(jié)的輸出為:
(7)
在上述公式中,需要設(shè)置調(diào)節(jié)的參數(shù)有PID反饋環(huán)節(jié)Kp,Ki,Kd的值。K的取值范圍為0~359,呈現(xiàn)周期性;系統(tǒng)提前步數(shù)N的選取需要根據(jù)實(shí)際跟蹤性能來(lái)優(yōu)化選取。
渦旋壓縮機(jī)軸向力電磁平衡系統(tǒng)中,存在一定程度的延時(shí),常規(guī)PID無(wú)法消除該延時(shí)量,僅僅由系統(tǒng)本身決定,針對(duì)系統(tǒng)提前步數(shù)量N的選擇需要一個(gè)評(píng)判標(biāo)準(zhǔn),以完成系統(tǒng)提前步數(shù)的最佳值確定。
筆者采用系統(tǒng)最小方差的性能評(píng)價(jià)方法,在渦旋壓縮機(jī)一個(gè)運(yùn)算周期內(nèi),系統(tǒng)不斷采樣當(dāng)前電磁吸力F(k)大小,并且與理想平衡力Fa(K)進(jìn)行比較。
PID超前預(yù)測(cè)平衡優(yōu)化性能指標(biāo)為:
(8)
系統(tǒng)超前預(yù)測(cè)步數(shù)N的值在閾值中不斷迭代尋優(yōu),以minErr值作為依據(jù),確定系統(tǒng)最適合的超前預(yù)測(cè)步數(shù)N的值,并將該參數(shù)保存至系統(tǒng)內(nèi)存中。當(dāng)系統(tǒng)啟動(dòng)時(shí),首先從處理器內(nèi)存中尋找是否有與當(dāng)前工況配合的N值,當(dāng)N存在時(shí),即可直接調(diào)用當(dāng)前N值進(jìn)行超前預(yù)測(cè)PID運(yùn)算,從而提高系統(tǒng)的運(yùn)算速度。
最佳提前量N自動(dòng)尋找的提前跟蹤PID算法工作流程圖如圖5所示。
圖5 自適應(yīng)提前步數(shù)量尋優(yōu)算法流程圖
基于前期做出相關(guān)分析工作之后,筆者確定在一定工況下系統(tǒng)的最大提前補(bǔ)償量的閾值M(目標(biāo)壓縮機(jī)最大延時(shí)時(shí)間),當(dāng)系統(tǒng)開(kāi)始運(yùn)行時(shí),系統(tǒng)依據(jù)當(dāng)前反饋量與理想值的誤差,進(jìn)行PID運(yùn)算并保存渦旋機(jī)主軸轉(zhuǎn)角一周期內(nèi)的誤差minErr。
依據(jù)上述式(8)評(píng)價(jià)系統(tǒng)控制效果,由最近3次系統(tǒng)誤差值判斷系統(tǒng)是否為穩(wěn)定狀態(tài)。當(dāng)系統(tǒng)穩(wěn)定時(shí)保存當(dāng)前指標(biāo),系統(tǒng)提前步數(shù)N迭代加一;當(dāng)系統(tǒng)迭代完成,搜索全部評(píng)價(jià)指標(biāo)值,找到誤差最小值,確定該值對(duì)應(yīng)的N值為最優(yōu)提前步數(shù)。
模擬實(shí)驗(yàn)平臺(tái)的關(guān)鍵參數(shù)如表1所示。
表1 模擬實(shí)驗(yàn)平臺(tái)關(guān)鍵參數(shù)
筆者分別采用常規(guī)PID控制器與帶自適應(yīng)提前步數(shù)的PID控制器,對(duì)渦旋壓縮機(jī)動(dòng)靜渦盤之間軸向平衡力曲線進(jìn)行跟蹤。
在控制系統(tǒng)調(diào)節(jié)過(guò)程中,采用試湊法,對(duì)PID參數(shù)進(jìn)行整定,具體為:Kp=10,Ki=0.01,Kd=0.05;控制的PWM周期T=50 μs。
在渦旋機(jī)軸向力平衡電磁伺服控制實(shí)驗(yàn)中,其動(dòng)態(tài)力跟蹤實(shí)驗(yàn)效果如圖6所示。
圖6 動(dòng)態(tài)力系統(tǒng)跟蹤控制效果
由圖6可知:在目標(biāo)曲線頻率50 Hz時(shí),在常規(guī)PID控制時(shí),系統(tǒng)存在明顯的相位滯后的現(xiàn)象;當(dāng)系統(tǒng)步數(shù)提前跟蹤算法控制介入時(shí),系統(tǒng)尋優(yōu)至最佳位置;
當(dāng)系統(tǒng)旋轉(zhuǎn)頻率為50 Hz時(shí),控制量最優(yōu)提前步數(shù)為20,目標(biāo)曲線頻率不同時(shí),所需的最優(yōu)提前步數(shù)也有所改變;如當(dāng)目標(biāo)曲線頻率為30 Hz時(shí),系統(tǒng)最優(yōu)提前步數(shù)變?yōu)?7。
除了最優(yōu)提前量的差異,當(dāng)頻率發(fā)生變化時(shí),系統(tǒng)控制誤差評(píng)價(jià)占比也有所變化。
因此,采用提前跟蹤控制可以有效地提高系統(tǒng)曲線的跟蹤效果。
在實(shí)驗(yàn)過(guò)程中,筆者記錄下每個(gè)控制周期中對(duì)應(yīng)的每個(gè)角度值的壓力的大小,將讀回的運(yùn)行數(shù)據(jù)的絕對(duì)誤差積分指標(biāo)IAE(integrated absolute error)進(jìn)行對(duì)比。
渦旋壓縮機(jī)軸向力平衡控制算法實(shí)驗(yàn)結(jié)果如表2所示。
表2 渦旋壓縮機(jī)軸向力平衡控制算法實(shí)驗(yàn)結(jié)果
根據(jù)表2可知:基于使用自適應(yīng)步數(shù)提前跟蹤PID控制跟蹤目標(biāo)曲線,可以提高系統(tǒng)的跟蹤運(yùn)行精度,通過(guò)不斷的迭代優(yōu)化尋優(yōu),找到最適合系統(tǒng)的提前步數(shù),最終系統(tǒng)的相位滯后的情況有了明顯的好轉(zhuǎn),改善了依靠單純的PID調(diào)節(jié),減少了電壓電流以及電流電磁力之間產(chǎn)生的系統(tǒng)相位差,提升了渦旋壓縮機(jī)渦旋盤之間的軸向力平衡的性能。
為實(shí)現(xiàn)渦旋壓縮機(jī)軸向動(dòng)態(tài)密封,需動(dòng)態(tài)平衡動(dòng)靜渦盤間周期性變化的軸向氣體分離力,在研究了渦旋壓縮機(jī)軸向氣體分離力變化規(guī)律后,筆者提出了采用電磁機(jī)構(gòu)產(chǎn)生電磁力快響應(yīng)動(dòng)態(tài)平衡的方案;針對(duì)高力密度帶鐵芯電磁力機(jī)構(gòu)動(dòng)態(tài)跟蹤控制時(shí)的滯后現(xiàn)象及其原因,提出了一種自適應(yīng)提前跟蹤PID閉環(huán)跟蹤控制方法,并搭建了相應(yīng)的實(shí)驗(yàn)系統(tǒng),對(duì)所提出的控制方法進(jìn)行了實(shí)驗(yàn)驗(yàn)證。
研究結(jié)果表明:
(1)采用電磁伺服機(jī)構(gòu),可快速地跟蹤目標(biāo)曲線,實(shí)現(xiàn)渦旋機(jī)軸向力的高性能動(dòng)態(tài)平衡;
(2)采用提前跟蹤的PID控制方法可以克服滯后問(wèn)題,較好地跟蹤周期性的動(dòng)態(tài)曲線;與傳統(tǒng)PID控制相比,可以提高系統(tǒng)的控制效果,準(zhǔn)確性得到了顯著的改善;
(3)對(duì)不同頻率或者周期,采用自適應(yīng)調(diào)節(jié)方法可以較好地按照相位滯后的大小,自動(dòng)調(diào)節(jié)系統(tǒng)的提前步數(shù)。
參考文獻(xiàn)(References):
[1] 楊啟超,李連生,趙遠(yuǎn)揚(yáng),等.渦旋壓縮機(jī)軸向動(dòng)態(tài)間隙的實(shí)驗(yàn)測(cè)量[J].西安交通大學(xué)學(xué)報(bào),2008,42(7):795-798,847.
[2] 王建吉,劉 濤,孫旖彤,等.渦旋壓縮機(jī)齒頂密封條摩擦磨損性能的研究[J].流體機(jī)械,2018,46(4):1-5.
[3] 余 洋.渦旋壓縮機(jī)動(dòng)力特性及仿真模擬研究[D].蘭州:蘭州理工大學(xué)能源動(dòng)力工程學(xué)院,2014.
[4] 曹 斌.基于分形理論的渦旋壓縮機(jī)渦旋盤端面摩擦特性研究[D].蘭州:蘭州理工大學(xué)石油化工學(xué)院,2018.
[5] 劉興旺,王 博,曾 強(qiáng),等.電動(dòng)汽車渦旋壓縮機(jī)變工況背壓平衡系統(tǒng)研究[J].流體機(jī)械,2019,47(8):16-22.
[6] 丁寶蒼.預(yù)測(cè)控制的理論與方法[M].北京:機(jī)械工業(yè),2008.
[7] 尹 超,李茂軍,張 靜.變頻恒壓供水系統(tǒng)的Smith預(yù)估復(fù)合控制[J].電機(jī)與控制應(yīng)用,2013,40(4)58-61.
[8] ADLY AA,MAYERGOYA I,BERGQVIST A. Preisachmodelinofmagnetostrictive hysteresis[J]. Journal of Applied Physics,1991,69(8):5777-5779.
[9] GE P,JOUANEH M. Generalized preisach model for hysteresis nonlinearity of piezoceramic actuator[J]. Precision Engineering,1997,20(2):99-111.
[10] Al J M, RAKHEJA S,SU Chun-yi.An analytical generalized Prandtl-Ishlinskii model inversion for hysteresiscompensation in micropositioning control [J].IEEE/ASME Transactions on Mechatronics,2011,16(4):734-744.
[11] 王湘江,王興松.基于KP模型的GMA遲滯系統(tǒng)辨識(shí)與補(bǔ)償[J].中國(guó)機(jī)械工程,2008,19(10):1167-1173.
[12] JILES D, ATHERTON D. Ferromagnetic hysteresis[J]. IEEE Transactions on Magnetics,1983,19(5):2183-2185.