田艷兵， 王濤， 王美玲， 張棟

(1.北京理工大學 自動化學院，北京 100083;2.青島理工大學自動化工程學院，山東青島 266520)

0 引言

精密伺服是超精密加工中的關(guān)鍵技術(shù)，是衡量一個國家超精密加工水平的重要指標。目前，國內(nèi)外超精密伺服控制精度已經(jīng)達到納米級。精密伺服定位系統(tǒng)多以電驅(qū)動，為避免摩擦對控制效果的影響，一般裝配有氣浮導軌。電機驅(qū)動分為直接驅(qū)動和間接驅(qū)動兩種模式。間接驅(qū)動借助滾軸絲杠等傳動機構(gòu)將電機轉(zhuǎn)動轉(zhuǎn)換為直線運動，其傳動效率較低，往返行程間隙大，限制了控制精度的提高;直接驅(qū)動采用直線電機驅(qū)動［1］，如音圈電機等［2］，傳動效率高，但是輸出力較小。2007年，東京工業(yè)大學的H.Shino等采用8個音圈電機驅(qū)動［3］，完成了5nm精度等級的二維精確定位，定位精度很高，行程受到限制。特別是對于二維定位多采用“H”型疊壓結(jié)構(gòu)，機械結(jié)構(gòu)復雜。有學者嘗試在一個平面機械結(jié)構(gòu)內(nèi)的二維定位，天津大學的馮曉梅等設(shè)計了基于音圈電機的氣浮平臺，精度能達到微米級［4］。

目前，氣動伺服系統(tǒng)多采用氣缸驅(qū)動，電磁干擾小，設(shè)備維護方便。近年來，氣缸驅(qū)動或者氣缸復合驅(qū)動的超精密定位系統(tǒng)在國內(nèi)外得到廣泛研究和運用，取得了一些成果。臺灣科技大學的Mao－Hsiung Chiang等采用氣缸和壓電驅(qū)動結(jié)合的模式［5］，較好的解決了控制精度和運動行程的問題，其控制精度達到微米級，但兩級進給模式使得機械結(jié)構(gòu)復雜，不利于控制精度的提高。

基于氣缸的定位伺服系統(tǒng)精度不高，往往需要配備雙級驅(qū)動機構(gòu)，在二維定位中采用的“H”型結(jié)構(gòu)后，使得機械平臺更復雜，各種承載平臺之間存在的摩擦不利于精密控制［6］。本文設(shè)計了一種基于新型金屬波紋管驅(qū)動的單級進給模式的氣動定位平臺，以氣浮平臺為支撐，采用高精密光柵尺進行位移反饋，簡化了系統(tǒng)機械結(jié)構(gòu)，降低了成本。針對系統(tǒng)在運動過程中存在的遲滯現(xiàn)象，建立基于PI的遲滯數(shù)學模型，為了改善系統(tǒng)的控制效果，采用了PI逆模型控制和反饋PID自整定控制，在2 500 μm的行程范圍內(nèi)，取得了較好的控制精度和響應速度。

1 硬件平臺

系統(tǒng)硬件主要由3部分組成。基于波紋管的運動驅(qū)動機構(gòu);氣浮導軌產(chǎn)生的靜壓軸承支撐的工作平臺以及位置檢測反饋部分。由單片機調(diào)節(jié)高精度壓力比例閥，控制波紋管內(nèi)部壓力，調(diào)整其伸縮量，控制平臺位移。為增加平臺剛度和穩(wěn)定性，在運動平臺另一側(cè)安裝有預緊彈簧。整個系統(tǒng)的結(jié)構(gòu)框圖和實物圖如圖1所示。

圖1 系統(tǒng)結(jié)構(gòu)框圖和實物圖Fig.1 Block diagram and actual picture of the system

1.1 運動驅(qū)動機構(gòu)

系統(tǒng)驅(qū)動機構(gòu)采用電氣比例閥控制的驅(qū)動波紋管。波紋管是一類子午線呈波紋狀的旋轉(zhuǎn)體，依靠內(nèi)部氣體壓力的變化產(chǎn)生伸縮運動，產(chǎn)生位移。按形狀可分為“U”型、“Ω”型、“S”型等［7］。其中“U”型波紋管承受壓力大，靜態(tài)特性較好，輸出位移與輸入信號線性關(guān)系好，無摩擦，模型簡單。通過單片機的12位D/A數(shù)據(jù)口控制高速比例閥，實現(xiàn)波紋管的位置控制［8］。

1.2 支撐體結(jié)構(gòu)

支撐體結(jié)構(gòu)主要由氣浮導軌和工作平臺組成，工作平臺由高剛度的多孔質(zhì)氣浮導軌產(chǎn)生靜壓空氣軸承支撐。

氣浮導軌為工作臺定向并提供無摩擦支撐。氣浮導軌的氣孔噴射壓縮空氣到工作臺表面產(chǎn)生懸浮力，當懸浮力與工作平臺重量平衡時，平臺處于穩(wěn)定懸浮狀態(tài)。支撐臺負載發(fā)生變化或有外界擾動時，氣浮平臺表面和載物臺之間的間隙改變，懸浮力重新達到平衡［9］，具有自平衡能力。該系統(tǒng)氣浮導軌為封閉矩形，氣浮導軌的上下、左右均有氣孔。上下兩側(cè)氣浮導軌會在上、下兩個面產(chǎn)生靜壓面，以提高靜壓軸承的剛性，左、右兩個面靜壓面起到導向作用。對于空氣靜壓軸承設(shè)計可以采用有限元分析方［10－11］，在此不再贅述，多孔質(zhì)空氣靜壓軸承的原理示意圖如圖2所示。

圖2 氣浮導軌Fig.2 Aerostatic slider

1.3 檢測元件

系統(tǒng)采用MicroE公司生產(chǎn)的Mercury3500系列的光柵尺進行位置反饋，其檢測分辨率可以達到5nm［12］。

2 系統(tǒng)遲滯模型

波紋管作為驅(qū)動機構(gòu)，近年來已有學者對其進行了一定的研究，其工作壽命和機械特性已經(jīng)在實驗中得到驗證［13］。作為執(zhí)行機構(gòu)，波紋管在低速或者靜態(tài)下，有著很好的線性，但是其采用氣體驅(qū)動，由于氣體本身的可壓縮性，在動態(tài)運動的情況下，波紋管驅(qū)動系統(tǒng)輸入和輸出之間存在遲滯現(xiàn)象。遲滯是指系統(tǒng)的運行軌跡在上行和下行之間存在間隙，也即系統(tǒng)的輸入和輸出之間并不是一一對應的關(guān)系，具體到波紋管來說，就是其電磁比例閥的輸入電壓所對應的波紋管產(chǎn)生的輸出位移不是一一對應的，這對于超精密定位系統(tǒng)有著較大的影響。

目前關(guān)于波紋管的遲滯特性大多基于金屬材料特性進行改進，來減小遲滯的影響，但是效果不好。本文對其遲滯特性進行研究，通過建立系統(tǒng)遲滯模型，找出針對性的控制策略，降低遲滯對系統(tǒng)定位精度的影響［14］。

系統(tǒng)遲滯特性由其初載曲線確定。在該系統(tǒng)中，其初載曲線可以通過實驗獲取，將驅(qū)動電磁壓力比例閥從無源狀態(tài)下開始，逐步施加驅(qū)動電壓到最大值過程中，采集位移參數(shù)，得到波紋管驅(qū)動平臺的運行軌跡。為了驗證波紋管的遲滯特性，平臺加載后，進行初載特性實驗，其實驗數(shù)據(jù)見表1。

波紋管驅(qū)動平臺的初載特性實驗在低頻輸入下進行，電壓的變化速率較慢。實驗中系統(tǒng)輸入電壓信號變化速率為0.5 V/s，系統(tǒng)的輸出響應可以很好的抑制波紋管的動態(tài)特性，增加系統(tǒng)模型的精度。

表1 實驗數(shù)據(jù)Table 1 Experimental data

2.1 PI模型

從表1的實驗數(shù)據(jù)中可以看出，波紋管驅(qū)動系統(tǒng)遲滯現(xiàn)象明顯，需要建立針對性的遲滯模型。目前常見的遲滯模型有Preisach模型、KP模型、PI模型等。PI模型是從Preisach模型發(fā)展而來的一種模型［15］，由Krasnosel’skii和Pokrovskii建立，由許多加權(quán)重疊的間隙算子(Backlash Operator)組成的，算子的特性由閾值r和權(quán)值ω決定，不同閾值的遲滯算子通過加權(quán)ω疊加，形成遲滯特性曲線。閾值r決定了基本遲滯算子的寬度，權(quán)值ω決定了基本遲滯算子的斜度［16］。

在構(gòu)建遲滯模型時，間隙算子的數(shù)學形式為

其中，T為采樣周期。式(1)是一個基于y(0)的遞歸公式，有

y0為系統(tǒng)在沒有激勵源的情況下的初值，一般情況下可以取0。

加權(quán)ω之后形成加權(quán)遲滯算子，即

從而PI遲滯模型可以通過n個加權(quán)遲滯算子累加得到，有

PI模型的一個顯著優(yōu)點是模型的解析逆容易得到，從而可以通過PI逆模型對控制對象線性化，線性化后的模型便于控制和分析，在線性化的基礎(chǔ)上進行控制方案的研究，其線性化原理如圖3所示［17］。

圖3 基于PI逆矩陣的控制器設(shè)計原理Fig.3 Controller design based on PI inverse matrix

2.2 模型參數(shù)辨識

從式(4)可以看出，決定系統(tǒng)PI模型的2個主要參數(shù)為r和ω。在實際操作過程中，PI模型參數(shù)的辨識可以利用系統(tǒng)帶載實驗數(shù)據(jù)形成的初載曲線進行，根據(jù)實驗數(shù)據(jù)，對式(5)所描述的數(shù)學模型，用Matlab的ployfit函數(shù)進行擬合，并用最小方差優(yōu)化，可以得到辨識后的系統(tǒng)參數(shù)如表2所示。

表2 參數(shù)辨識結(jié)果Table 2 Result of parameters identification

結(jié)合PI模型，其解析逆模型參數(shù)也容易求得，可以建立基于逆模型的前饋控制器，能夠較好的完成軌跡跟蹤。

2.3 模型驗證

為了驗證建立模型的精度，在實驗室平臺上進行實驗［18］。仿真實驗中將系統(tǒng)模型放在SIMULINK中進行，壓力比例閥輸入信號和壓力輸出按線性關(guān)系處理。而對于實驗平臺則是按照一定規(guī)律輸入電壓信號，通過光柵尺檢測位移輸出。上位機軟件采用Labwindows進行設(shè)計，下位機采用MSP430單片機系統(tǒng)對高精度壓力比例閥進行控制，調(diào)整壓力，控制波紋管伸縮量，并和上位機進行通信。在實驗中，電磁比例閥控制電壓信號變化軌跡為(0.0v－1.5v－0.0v－1.4v－0.2v－1.0v－0.5v－0.8v－0.6v)。

變化速率為0.5V/s，實驗輸出和模型輸出如圖4所示。

圖4 模型輸出和實驗輸出Fig.4 Model output and experimental output

從圖中可以看出，通過初載曲線擬合建立的PI遲滯模型可以很好地逼近系統(tǒng)的實際輸出，所建立的遲滯模型是有效的，精度較高。

3 控制策略研究

實驗表明，在波紋管驅(qū)動特性中，遲滯特性是影響其控制精度的主要因素，因此在控制器設(shè)計中，重點考慮遲滯特性的線性化，根據(jù)PI模型采用PI逆模型進行處理。

3.1 前饋開環(huán)控制

將波紋管驅(qū)動系統(tǒng)辨識得到的遲滯模型進行求逆，逆模型直接串聯(lián)在被控對象前面，將系統(tǒng)線性化，逆模型作為前饋控制器。這種方案稱為前饋控制方案，盡管缺少反饋環(huán)節(jié)，但是在某些精度要求不高的場所，該種方案可以減少位置檢測反饋裝置，從而從很大程度上降低成本。實驗研究表明，采用直接前饋開環(huán)控制，控制精度可以達到十幾個微米的等級，因此也有一定的實用價值。

考慮系統(tǒng)在運動過程中存在的動態(tài)特性影響，可以在前饋控制器前加入針對動態(tài)特性的控制環(huán)節(jié)，作為總的前饋控制器。波紋管驅(qū)動精密定位平臺前饋控制框圖如圖5所示。

圖5 系統(tǒng)遲滯逆模型開環(huán)控制Fig.5 Open-loop control based on inverse hysteresis model

3.2 基于PSO自整定的復合控制方案

為了進一步提高系統(tǒng)的控制精度，在定位系統(tǒng)中采用了前饋開環(huán)控制加反饋環(huán)節(jié)的控制方案。前饋控制用來校正波紋管驅(qū)動執(zhí)行器的遲滯非線性，提高對參考位移信號的跟蹤能力;反饋調(diào)節(jié)用來進一步校正前饋補償沒有消除的偏差以及由模型的不確定性帶來的誤差，從而彌補建模帶來的定位誤差，使得工作平臺跟蹤精度進一步提高［19］。

反饋環(huán)節(jié)采用數(shù)字PID控制，PID參數(shù)對于控制效果非常重要，但參數(shù)整定較為繁瑣，且效果不一定最佳，特別是系統(tǒng)模型參數(shù)變化的情況下，PID參數(shù)需要針對不同情況進行整定。對于定位平臺，由于其承載質(zhì)量的不同，不同的PID參數(shù)會影響控制效果［20］。為了改善參數(shù)整定效果，考慮到運動控制系統(tǒng)調(diào)整周期短，在復合控制中，通過粒子群算法對PID參數(shù)進行優(yōu)化，確定系統(tǒng)最優(yōu)參數(shù)，可以將不同負載時的控制參數(shù)預先整定儲存，以改善系統(tǒng)的實時性能，提高系統(tǒng)的控制精度。

粒子群算法(Particle Swarm Optimization，PSO)是遺傳算法的一種，和遺傳算法相似，它也是從一定范圍內(nèi)的隨機解出發(fā)，通過迭代尋找最優(yōu)解，通過適應度函數(shù)－fitness來評價辨識參數(shù)的品質(zhì)，但它比遺傳算法規(guī)則更為簡單，它沒有遺傳算法的“交叉”(Crossover)和“變異”(Mutation)操作，主要是通過模擬鳥群尋找食物的過程，提出的優(yōu)化算法［21］。該算法簡單，編程較容易實現(xiàn)。PSO算法的關(guān)鍵在于適宜度函數(shù)的設(shè)計，根據(jù)運動控制系統(tǒng)的特點和控制效果，適宜度函數(shù)可以選擇為

其中，W1，W2，W3，W4為權(quán)重函數(shù)，在尋優(yōu)過程中，充分考慮運動定位控制的特點，對系統(tǒng)的穩(wěn)態(tài)誤差和調(diào)節(jié)時間作為2個主要衡量指標，權(quán)值分別取值為50，30，2，10。

對于該適宜度函數(shù)，采用PSO算法進行尋優(yōu)［22］，目的是使適宜度函數(shù)達到極小值。PID控制效果取決于3個參數(shù)Kp、Ki、Kd，因此訓練的粒子維數(shù)為3，粒子數(shù)目取30個，最大迭代次數(shù)2 000次。PSO算法參數(shù)整定的具體實現(xiàn)步驟為:

1)系統(tǒng)初始化，設(shè)定Kp、Ki、Kd初值，確定Kp、Ki、Kd的范圍，賦予粒子初始位置和速度以及pbest和gbest，確定系統(tǒng)的合法的Vmax和Xmax。

2)對每個粒子Kp、Ki、Kd階躍響應，實驗得到時域內(nèi)的性能指標e(∞)、ts、tr和σ%。

3)根據(jù)階躍響應參數(shù)，據(jù)式(6)計算種群中每個粒子適應度值。

4)對于每一個粒子，根據(jù)計算的適宜度函數(shù)，尋求更優(yōu)值，確定是否更新其pbest和gbest值。

5)更新粒子的速度和位置。

6)判斷粒子位置和速度的合法性。

7)如果達到最大迭代次數(shù)轉(zhuǎn)至8)，否則轉(zhuǎn)至2)。

8)尋優(yōu)結(jié)束，得到PSO尋優(yōu)后PID的參數(shù)，即Kp、Ki、Kd。

PSO尋優(yōu)的PID復合控制方案框圖如圖6所示。

圖6 基于PSO的復合控制方案Fig.6 Compound control based on PSO

4 實驗和分析

分別利用傳統(tǒng)PID控制和自尋優(yōu)PID復合控制方案控制，在波紋管驅(qū)動平臺上進行實驗，系統(tǒng)主要參數(shù)為

m為移動平臺的自重和負重;A為波紋管有效受力面積;ps為氣源壓力;p0為靜壓軸承氣源壓力。

4.1 PID控制方案

對于傳統(tǒng)的PID控制，進行實驗研究，取Kp=0.46，Ki=0.06，Kd=0.01，其階躍穩(wěn)態(tài)誤差為0.16 μm、超調(diào)量為3.7%，如圖7所示。

圖7 PID控制階躍響應Fig.7 Step response

三角波跟蹤和正弦跟蹤平均誤差分別為0.41 μm和0.67 μm，見表3。

表3 控制效果比較Table 3 Positioning accuracy of different control schemes

4.2 基于PSO的PID自整定方案

實驗按照定位實驗和軌跡跟蹤實驗進行，階躍響應選取30 μm階躍信號，其實驗響應曲線如圖8所示，系統(tǒng)的調(diào)整時間0.06 s左右，超調(diào)在2%左右，穩(wěn)定誤差小于0.05 μm。

圖8 PSO尋優(yōu)PID階躍響應Fig.8 Step response based on PSO

跟蹤實驗采用三角波和正弦波輸入。三角波實驗時，波紋管預加壓伸縮到1 000 μm，誤差曲線以1 000 μm為基準。正弦波輸入頻率為2.5 Hz，避免負值輸入，輸入從零值開始，數(shù)值偏向時間軸一側(cè)。其軌跡跟蹤效果如圖9和圖10所示。

圖9 鋸齒波跟蹤實驗Fig.9 Sawtooth tracking experiment

從圖中可以看出在PSO優(yōu)化參數(shù)后的復合控制方案下，系統(tǒng)跟蹤效果很好。三角波輸入，跟蹤最大誤差為0.32 μm，平均為0.12 μm;正弦波輸入，跟蹤誤差平均為0.18 μm，達到亞微米級或納米級。滿足超精密定位平臺的需求。與傳統(tǒng)PID控制相比，其平均誤差和最大誤差均有較大改進，見表3。

圖10 正弦跟蹤實驗Fig.10 Sine tracking experiment

5 結(jié)論

本文設(shè)計了一種基于波紋管驅(qū)動的超精密運動平臺，實驗表明，實驗平臺定位精度達到了超精密定位的要求。波紋管驅(qū)動的定位平臺在國內(nèi)外研究和應用較少，其主要優(yōu)點有:

1)波紋管驅(qū)動相比較電動驅(qū)動而言，波紋管采用氣動作為驅(qū)動能量，環(huán)境干擾小，維護方便。

2)波紋管本身在較大行程中有著很高的定位精度，這一點是其他超精密驅(qū)動裝置，如壓電等所部具備的，因此在一些較大行程的超精密驅(qū)動中可以避免兩級驅(qū)動的模式。采用PSO進行參數(shù)優(yōu)化之后的復合控制方案，控制精度較高。

3)波紋管造價低廉，有著明顯的價格和成本優(yōu)勢。

4)波紋管能量來源為氣源，相對于電動而言，其驅(qū)動力的大小取決于管腔內(nèi)的壓力和波紋管截面積，因此在同等體積下，其輸出推力相對于電機推力較大。

5)由于波紋管本身的柔性特性，可以考慮在一個平面內(nèi)實現(xiàn)二維定位，從而避免的傳統(tǒng)的“H”型結(jié)構(gòu)，降低了系統(tǒng)的機械結(jié)構(gòu)復雜程度，便于設(shè)計和控制，也降低了平臺的硬件成本。

波紋管作為一種新型的驅(qū)動器件，研究表明其在超精密控制中有很好的應用前景，但是也存在諸多問題，目前還沒有專門的企業(yè)生產(chǎn)驅(qū)動用波紋管，本實驗中用到的波紋管為專門設(shè)計，其機械性能還有待進一步研究，其機械壽命也還需要進一步考慮和驗證，希望更多地專家和學者給予關(guān)注。

［1］張從鵬，劉強.直線電機氣浮精密定位平臺設(shè)計與控制［J］.北京航空航天大學學報，2008，34(2):224－228.

ZHANG Congpeng，LIU Qiang.Design and control of air bearing precision positioning stage driven by linear motors［J］.Journal of Beijing University of Aeronautics and Astronautics，2008，34(2):224－228.

［2］趙興玉，張勝泉，張大衛(wèi).基于音圈電機精密定位平臺的控制系統(tǒng)設(shè)計與仿真［J］，天津大學學報，2007，40(2):127－131.

ZHAO Xingyu，ZHANG Shengquan，ZHANG Dawei.Design and simulation of the control system of precision positioning table driven by voice coil actuator［J］.Journal of Tianjin University，2007，40(2):127－131.

［3］SHINNO H，YOSHIKA H，TANIGUCHI K.A newly developed linear motor-driven aerostatic x-y planar motion table system for nano-machining［J］.Annals of the CIRP，2007，56(1):369－372.

［4］馮曉梅，張大衛(wèi)，趙興玉，等.基于音圈電機的新型高速精密定位系統(tǒng)設(shè)計方法［J］.中國機械工程，2005，16(16):1414－1418.

FENG Xiaomei，ZHANG Dawei，ZHAO Xingyu，et al.The new high-speed precision positioning system design methodology based on voice coil motor［J］.Chinese Mechanical Engineering，2005，16(16):1414－1418.

［5］CHIANG M，CHEN C，TSOU T.Large stroke and high precision pneumatic-piezoelectric hybrid positioning control using adaptive discrete variable structure control［J］.Mechtronics，2005(15):523－545.

［6］張建明.現(xiàn)代超精密加工技術(shù)和裝備的研究與發(fā)展［J］.航空精密制造技術(shù)，2008，44(1):1－7.

ZHANG Jiangming.Research and development of modern ultraprecision processing technology and equipment［J］.Aviation Precision Manufacturing Technology，2008，44(1):1－7.

［7］王占彬，周浩洋，石朝鋒，等.一種U型波紋管強度及穩(wěn)定性分析方法研究［J］.強度與環(huán)境，2012，39(4):19－24.

WANG Zhanbin，ZHOU Haoyang，SHI Chaofeng，et al.Research on the analytical method of U-shape bellows strength and Stabilization［J］.Structure and Environment Engineering，2012，39(4):19－24.

［8］樊大均.波紋管設(shè)計學［M］.北京:北京理工大學出版社，1988:60－93.

［9］劉強，張從鵬.直線電機驅(qū)動的H型氣浮導軌運動平臺［J］.光學精密工程，2007，15(10):1540－1545.

LIU Qiang，ZHANG Congpeng.H-type air-bearing motion stage driven by liner motors［J］.Optics and Precision Engineering，2007，15(10):1540－1545.

［10］陶家生.閉式雙排節(jié)流孔矩形氣浮導軌的計算方法［J］.光學精密工程，1996，4(2):31－36.

TAO Jiasheng.Calculation of close double row orifice rectangular aerostatic slideways［J］.Optics and Precision Engineering，1996，4(2):31－36.

［11］張新宇，陳忠基，姚瑤，等.基于有限元法的氣浮支撐系統(tǒng)的數(shù)值模擬與實驗研究［J］.遼寧科技大學學報，2008，31(3/4):285－287.

ZHANG Xinyu，CHEN Zhongji，YAO Yao，et al.Numerical simulation based on finite element method and experiment of aerostatic bearing system［J］.Journal of University of Science and Technology Liaoning，2008，31(3/4):285－287.

［12］王文，李欣欣.超精密定位平臺的測量系統(tǒng)研究［J］.機電工程，2006，23(4):13－16.

WANG Wen，LI Xinxin.Study on measurement system of ultraprecision positioner［J］.Journal of Mechanical and Electrical Engineering，2006，23(4):13－16.

［13］KAWASHIMA K，ARAIA T，TADANO K，et al.Development of coarse/fine dual stage using pneumatically driven bellows actuator and cylinder with air bearings［J］.Precision Engineering，2010，34(3):526－533.

［14］張棟，張承進，魏強.壓電微動工作臺的動態(tài)遲滯模型［J］.光學精密工程，2009，17(3):549－556.

ZHANG Dong，ZHANG Chenjin，WEI Qiang.Dynamic hysteresis model of piezo-positioning stage［J］.Optics and Precision Engineering，2009，17(3):549－556.

［15］賴志林，劉向東，陳振，等.壓電陶瓷執(zhí)行器多模時變滑模逆補償控制［J］.電機與控制學報，2012，16(1):92－103.

LAI Zhilin，LIU Xiangdong，CHEN Zhen，et al.Multi-mode time-varying sliding mode control of Piezoceramic actuator based on inverse compensation［J］.Electric Machines and Control，2012，16(1):92－103.

［16］曹金鳳，黨選舉，張凱.基于PI和神經(jīng)網(wǎng)絡混合模型的音圈電機遲滯建模［J］.系統(tǒng)仿真學報，201l，123(2):386－389.

CAO Jinfeng，DANG Xuanju，ZHANG Kai.Hybrid hysteresis modeling of voice coil motor based on PI and neural network［J］.Journal of System Simulation，201l，123(2):386－389.

［17］ANG Weitech，KHOSLA P K，RIVIERE C N.Feed forward controller with inverse rate-dependent model for piezoelectric actuators in trajectory-tracking applications［J］.IEEE/ASME Transactions on Mechatronics，2007，12(2):134－142.

［18］黨選舉，姜輝，楊青，等.高頻響音圈電機的復雜遲滯建模與驗證［J］.微電機，2012，45(6):1－5.

DANG Xuanju，JIANG Hui，YANG Qing，et al.Modeling and verify of complex hysteresis for high frequency voice coil motor［J］.Micromotors，2012，45(6):1－5.

［19］王希花，郭書祥，葉秀芬，等.壓電陶瓷遲滯特性的建模及復合控制［J］.電機與控制學報，2009，13(5):766－771.

WANG Xihua，GUO Shuxiang，YE Xiufen，et al.Modeling and feed-forward control based on piezoelectric ceramic hysteretic［J］.Electric Machines and Control，2009，13(5):766－771.

［20］崔玉國，董維杰，孫寶元，等.壓電微動工作臺的位移復合控制［J］.機械工程學報，2006，42(3):156－161.

CUI Yuguo，DONG Weijie，SUN Baoyuan，et al.Displacement compounding control of piezoelectric micro-motion worktable［J］.Optics and Precision Engineering，2006，42(3):156－161.

［21］KENNEDY J，EBERHART R.Particle swarm optimization［C］//Proceedings of IEEE international conference on neural networks，Nov 27－Dec 1，1995，Perth，Australia.1995，4(2):1942－1948.

［22］田艷兵.基于PSO的自整定PID溫度控制研究［J］.化工儀表自動化，2010，37(11):31－33.

TIAN Yanbing.The design of self-tuning PID temperature controller based on PSO［J］.Control and Instruments in Chemical Industry，2010，37(11):31－33.