王曉慧 丁智 劉寶權(quán) 王軍生 張巖

(1鞍鋼股份技術(shù)中心，鞍山 114009)(2鞍鋼集團(tuán)自動化公司，鞍山 114009)

在追求產(chǎn)品小型化、集成化的今天，出現(xiàn)了以光學(xué)為理論基礎(chǔ)，以微機(jī)械加工為手段制作的微光機(jī)電一體化系統(tǒng)，并被廣泛應(yīng)用于信息處理、空間技術(shù)、光纖通訊等諸多領(lǐng)域［1］.而具有微結(jié)構(gòu)表面的光學(xué)元件由于其光學(xué)特性、摩擦性、耐磨損性等特定性能，已經(jīng)成為微光機(jī)電產(chǎn)品的關(guān)鍵零部件［2-3］.由于傳統(tǒng)機(jī)床很難滿足微結(jié)構(gòu)表面車削加工過程中對系統(tǒng)高頻響、高分辨率、高次軌跡聯(lián)動功能的要求，近年來，利用快速伺服刀架(FTS)作為精密微位移加工模塊來車削加工微結(jié)構(gòu)表面的方法已經(jīng)成為微結(jié)構(gòu)表面切削加工的一種主流技術(shù)［4-6］.壓電陶瓷由于具有分辨率高、響應(yīng)快、剛度大等優(yōu)點(diǎn)，成為常用的FTS驅(qū)動元件，但是壓電陶瓷本身所固有的遲滯、蠕變等非線性特點(diǎn)不但會降低FTS系統(tǒng)的控制精度，而且可能造成系統(tǒng)失穩(wěn)［7］.為了減小遲滯非線性的影響，提高 FTS的跟蹤控制精度，必須對其進(jìn)行合理建模，從而設(shè)計(jì)相應(yīng)的控制策略.近年來，針對精密制造系統(tǒng)中遲滯非線性系統(tǒng)的建模日益成為各國學(xué)者關(guān)注的焦點(diǎn).文獻(xiàn)［8］在經(jīng)典Preisach模型的基礎(chǔ)上建立了廣義Preisach模型，大幅度提高了壓電陶瓷的跟蹤精度.文獻(xiàn)［9-10］將神經(jīng)網(wǎng)絡(luò)引入到遲滯非線性的建模中，通過引入Preisach型遲滯算子的概念，利用神經(jīng)網(wǎng)絡(luò)來離散化Preisach模型，簡化模型的辨識過程.但是由于Preisach模型屬于靜態(tài)模型，即遲滯輸出僅僅受到輸入極值的影響，而與輸入信號的變化頻率無關(guān)，即速率無關(guān)(rate-independent)，很難精確地對快速伺服刀架(FTS)的動態(tài)遲滯特性進(jìn)行建模.因此，針對FTS的動態(tài)遲滯特性，本文建立了FTS遲滯系統(tǒng)的RBF神經(jīng)網(wǎng)絡(luò)模型.由于理論上神經(jīng)網(wǎng)絡(luò)不能辨識壓電陶瓷驅(qū)動器的遲滯特性這類多值映射的非線性現(xiàn)象［11］，因此本文采用拓展輸入空間的方法通過引入指數(shù)型遲滯算子，將遲滯算子的輸出與FTS系統(tǒng)的輸入一起作為RBF神經(jīng)網(wǎng)絡(luò)的輸入向量，實(shí)現(xiàn)FTS遲滯系統(tǒng)由多值映射到單值映射的轉(zhuǎn)換，進(jìn)而利用神經(jīng)網(wǎng)絡(luò)對其進(jìn)行建模.

1 快速伺服刀架動態(tài)遲滯特性

快速伺服刀架作為微結(jié)構(gòu)表面車削加工的關(guān)鍵部件，其具體結(jié)構(gòu)如圖1所示.系統(tǒng)選用德國PI公司的P-841.20型壓電陶瓷作為FTS的驅(qū)動元件，其最大伸長量為30 μm，分辨率可達(dá)0.3 nm，空載頻響為14 kHz，采用德國PI公司的D-050型電容測微儀作為其實(shí)時位移檢測裝置，其量程為50.8 μm，測量分辨率為0.1 nm.

圖1 快速伺服刀架結(jié)構(gòu)示意圖

快速伺服刀架的遲滯特性主要體現(xiàn)在其電壓位移曲線.由圖2可見，F(xiàn)TS在開環(huán)狀態(tài)下其升壓位移曲線與降壓位移曲線并不一致，同樣的輸入電壓因?yàn)榧訅簹v史的不同對應(yīng)不同的輸出位移，同樣的輸出位移對應(yīng)不同的輸入輸出曲線，而且其遲滯曲線寬度隨著輸入電壓幅值的增加而增大.

FTS的動態(tài)遲滯非線性特點(diǎn)體現(xiàn)在其驅(qū)動電壓頻率依賴性上，即速率相關(guān)(rate-dependent)性.改變輸入信號頻率，幅值60 V，輸入頻率分別為0.1，1.0，2.5，5 Hz，測試其輸入輸出特性，如圖3所示，隨著輸入電壓頻率的改變，F(xiàn)TS遲滯曲線的形狀和寬度也隨之改變.隨著輸入電壓信號頻率的增加，其最大輸出位移值也隨之變小，頻率為0.1 Hz輸入信號的最大輸出位移值比頻率為5 Hz輸入信號的最大輸出位移值大0.2 μm左右；同時，隨著輸入信號頻率的增加，其遲滯曲線的寬度也略有增加.

圖2 FTS遲滯位移曲線

圖3 FTS輸出頻率特性

由此可見，F(xiàn)TS的輸出位移不僅與當(dāng)前輸入電壓和歷史輸入電壓有關(guān)，還與輸入電壓信號的頻率有關(guān)，其動態(tài)遲滯特性是具有速率相關(guān)性的典型多對多非線性映射問題.因此，建立能夠描述其動態(tài)遲滯特性，并具有常規(guī)逼近技術(shù)優(yōu)點(diǎn)的遲滯模型是本文所做的研究工作.

2 基于RBF神經(jīng)網(wǎng)絡(luò)的FTS遲滯模型

2.1 拓展輸入空間

FTS的遲滯特性是一個多值映射問題，假設(shè)輸入輸出信號分別滿足u≤C1，y≤C2其中C1，C2為常數(shù)，多值遲滯映射定義如下:

定義1 y=f(u(t))為多值映射遲滯函數(shù)，如果滿足以下條件:

式中，t1，t2表示在輸入值u大小相同時，正、逆程對應(yīng)的時間；y(t1)表示輸入為u(t1)時的y值.對于極值時對應(yīng)的特殊點(diǎn)(見圖4)滿足:

t1=t2，u(t1)=u(t2)，y(t1)=y(t2)=y(t*)

圖4 多值映射遲滯特性

通過分析FTS系統(tǒng)的遲滯位移曲線發(fā)現(xiàn)，雖然遲滯曲線的形狀受到輸入信號幅值及頻率的影響，但是曲線上動點(diǎn)的運(yùn)動規(guī)律卻是非常的接近，在輸入到達(dá)某一極值點(diǎn)的前后分別沿著類似的兩條曲線運(yùn)動.根據(jù)遲滯曲線的這一特點(diǎn)，可以用形狀相似的數(shù)學(xué)曲線，即構(gòu)造遲滯算子f(x)來描述系統(tǒng)的遲滯動態(tài)特性.這里利用指數(shù)函數(shù)來構(gòu)造遲滯算子，其數(shù)學(xué)表達(dá)式為

式中，x為系統(tǒng)輸入；f(x)為遲滯算子輸出；x0為距離x最近的前一個輸入極值；f(x0)為輸入為x0時遲滯算子的輸出.

為了更好地跟蹤不同形狀遲滯曲線的變化趨勢，對遲滯算子進(jìn)行了如下改進(jìn):

在式(4)中，由于σ的值恒為1，很難對不同形狀的遲滯曲線都進(jìn)行高精度的辨識，因此增加一個調(diào)整遲滯算子曲線形狀的系數(shù)σ，以便于能夠更精確的跟蹤遲滯曲線的變化趨勢.由式(5)可知

將FTS系統(tǒng)的輸入x(t)與遲滯算子的輸出f［x(t)］共同作為神經(jīng)網(wǎng)絡(luò)的輸入向量，神經(jīng)網(wǎng)絡(luò)的輸入空間就由原來的一維增加到了兩維，稱之為輸入空間的拓展.這樣，就可以在多維空間中唯一地確定遲滯特性的每一個狀態(tài)，建立FTS的遲滯特性模型.

2.2 基于RBF神經(jīng)網(wǎng)絡(luò)的遲滯模型

由于徑向基函數(shù)(RBF)神經(jīng)網(wǎng)絡(luò)具有計(jì)算量小、收斂速度快、無局部極小等優(yōu)點(diǎn)，因此本文利用RBF神經(jīng)網(wǎng)絡(luò)作為遲滯非線性系統(tǒng)的建模工具.利用拓展輸入空間法建立的FTS遲滯模型基本結(jié)構(gòu)如圖5所示，將FTS系統(tǒng)輸入x(t)與遲滯算子輸出f［x(t)］一同作為RBF神經(jīng)網(wǎng)絡(luò)的輸入向量.

圖5 拓展輸入空間法建立FTS系統(tǒng)遲滯模型

3 模型的仿真研究

為了驗(yàn)證模型的有效性，對FTS系統(tǒng)施以頻率為1 Hz，幅值為0～60 V的正弦周期激勵，采樣頻率為1 kHz，采集2000個實(shí)測數(shù)據(jù)對，隨機(jī)挑選其中的1500個數(shù)據(jù)對用來訓(xùn)練神經(jīng)網(wǎng)絡(luò)，其余500個數(shù)據(jù)對用來檢驗(yàn)?zāi)Ｐ托Ч?由于采集的數(shù)據(jù)單位不一致且均為正值，因而在訓(xùn)練之前必須對數(shù)據(jù)進(jìn)行歸一化處理，歸一化后的數(shù)據(jù)將分布在［0，1］區(qū)間內(nèi).利用拓展輸入空間法進(jìn)行系統(tǒng)辨識的結(jié)果如圖6(a)所示.圖6(b)為實(shí)測數(shù)據(jù)的辨識誤差，模型的驗(yàn)證均方差MSE=5.1633×10－6.

圖6 RBF神經(jīng)網(wǎng)絡(luò)模型的辨識結(jié)果及誤差

4 結(jié)語

本文針對快速伺服刀架(FTS)的動態(tài)遲滯特性，通過引入指數(shù)型遲滯算子，采用拓展神經(jīng)網(wǎng)絡(luò)輸入空間的方法成功實(shí)現(xiàn)了FTS遲滯系統(tǒng)由多值映射到單值映射的轉(zhuǎn)換，并利用RBF神經(jīng)網(wǎng)絡(luò)對其進(jìn)行建模.實(shí)驗(yàn)結(jié)果表明，該遲滯模型可以很好地預(yù)測快速伺服刀架的遲滯特性，模型的驗(yàn)證均方差MSE為5.1633×10－6.

References)

［1］王仕璠，朱自強(qiáng).現(xiàn)代光學(xué)原理［M］.成都:電子科技大學(xué)出版社，1998:2-12.

［2］Asoh H，Oide A，Ono S.Formation of microstructured silicon surfaces by electrochemical etching using colloidal crystal as mask［J］.Electrochemistry Communications，2006，8(12):1817-1820.

［3］Huang Ying，Liu Ran，Lai Jianjun，et al.Design and fabrication of negative microlens array［J］.Optics＆Laser Technology，2008，40(8):1047-1050.

［4］Sze-Weia Gan，Han-Seok Lim，Rahman M，et al.A fine tool servo system for global position error compensation for a miniature ultra-precision lathe［J］.Machine Tools and Manufacture，2007，47(7/8):1302-1310.

［5］Gao W，Tano M，Takeshi A，et al.Measurement and compensation of error motions of a diamond turning machine［J］.Precision Engineering，2007，31(3):310-316.

［6］楊元華.基于FTS的微結(jié)構(gòu)功能表面超精密切削加工關(guān)鍵技術(shù)［D］.哈爾濱:哈爾濱工業(yè)大學(xué)，2007.

［7］Li Chuntao，Tan Yonghong.A neural networks model for hysteresis nonlinearity［J］.Sensors and Actuators，2004，7112(3):49-54.

［8］Ge P，Jouaneh Musa.Generalized preisach model for hysteresis nonlinearity of piezoceramic actuators［J］.Precision Engineering，1997，20(2):99-111.

［9］Adly A A，Abd-El-Hafiz S K.Using neural networks in the identification of preisach-type hysteresis models［J］.IEEE Trans Magnetics，1998，34(3):629-635.

［10］Claudio S，Ciro V.Magnetic hysteresis modeling via feed-forward neural networks［J］.IEEE Trans Magnetics，1998，34(3):623-628.

［11］Wei J D，Sun C T.Constructing hysteretic memory in neural networks［J］.IEEE Trans Systems，Man and Cybernetics Part B:Cybernetics，2000，30(4):601-609.