方 遠，劉樹勇，楊理華，蘇 攀

（1.海軍工程大學(xué) 動力工程學(xué)院，武漢 430033；2.海軍潛艇學(xué)院 動力操縱系，山東 青島 266199）

壓電陶瓷材料具有機-電能量轉(zhuǎn)化效率高、穩(wěn)定性好及靈敏度高等優(yōu)良特性，既可以制成傳感器，也可以制成作動器[1-2]。將壓電陶瓷材料制成作動器時，其類型可分為疊堆型和片狀型。

疊堆型壓電陶瓷作動器主要用于超精密加工、微位移器控制等定位精度要求高的場合。為提高壓電作動器在振動主動控制中的精度，實現(xiàn)作動器電壓與位移高線性度要求，就必須控制壓電材料遲滯特性的影響。對此，Goldfarb和Celannovic利用Maxwell模型對疊堆型作動器遲滯特性進行建模，通過該模型對控制過程進行設(shè)計和分析，實現(xiàn)了較好的軌跡跟蹤效果[3]。王代華等結(jié)合Bouc-Wen模型和參數(shù)識別方法，對疊堆型作動器的輸入電壓和輸出位移之間的遲滯關(guān)系進行線性化處理，使作動器的非線性度下降10%左右[4]。李巍基于疊堆型作動器的精密運動平臺，提出非對稱Bouc-Wen模型參數(shù)辨識方法，該方法比傳統(tǒng)Bouc-Wen模型的建模精度更高[5]。

片狀壓電陶瓷作動器主要應(yīng)用于對空間體積、質(zhì)量以及受力性質(zhì)有特殊要求的場合。Tarazaga等采用PZT片，對太空環(huán)境中的剛性充氣臂振動進行主動控制，使其末端振幅減小93%[6]。黃建用壓電纖維復(fù)合材料制成作動器，對負泊松比機翼結(jié)構(gòu)進行扭轉(zhuǎn)控制，從而提升了固定翼飛行器的氣動性能[7]。孫龍飛等針對船舶、航天器中的錐形圓筒結(jié)構(gòu)，利用PZT片對一端固定的錐形圓筒進行振動主動控制，采用模糊滑?？刂剖箞A筒自由端振幅減小約50%[8]。

目前，將遲滯補償理論和疊堆型作動器控制結(jié)合的研究成果較多，然而將遲滯補償理論和片狀型作動器控制結(jié)合的應(yīng)用較為鮮見。本文提出一種通過PZT片與懸臂梁形變關(guān)系測算PZT片位移的方法，由此得到PZT片電壓-位移關(guān)系，并基于多項式擬合逆遲滯補償算法對PZT片的控制電壓進行補償，提高壓電懸臂梁的主動控制效果。

1 PZT片遲滯特性分析

1.1 懸臂梁偏轉(zhuǎn)方程

為了研究PZT片遲滯特性，需得到其電壓-位移關(guān)系。由于實驗條件限制，無法對該關(guān)系進行精確測量。因此，本文利用PZT片與懸臂梁形變關(guān)系來測算PZT片在不同電壓下的位移。

試驗中PZT片只接受正電壓作用，通電時發(fā)生縮短形變，導(dǎo)致懸臂梁產(chǎn)生偏轉(zhuǎn)或擺動，自由端產(chǎn)生位移量。PZT片與懸臂梁形變關(guān)系如圖1所示。

圖中r是懸臂梁自由端位移量，l是壓電PZT片中心到自由端距離，h是懸臂梁厚度，s為PZT片形變量。測得PZT片l=148.4 mm，h=0.4 mm。

在進行PZT片電壓-位移測試時，為了便于研究和計算，需作如下假設(shè)：

(1)懸臂梁表面長度變化量等于PZT片長度變化量s，且僅在長度方向上發(fā)生形變；

(2)懸臂梁發(fā)生彎曲形變起點位置為PZT片中心線；

(3)梁模型為Euler-Bernoulli梁；

(4)忽略粘接PZT片給懸臂梁帶來的影響。

基于以上假設(shè)，當Euler-Bernoulli梁在發(fā)生彎曲形變時，形變前垂直于梁軸線的平面在形變后保持為平面，并且仍然垂直于梁軸線，具體如圖2所示。

圖2 歐拉梁的形變

變形前相距dx的兩個橫截面，在變形后產(chǎn)生了夾角dθ，在橫截面垂直方向上與中性層的距離為Z處有軸向纖維則的長度變量為cc′

因此可算得PZT片形變量為

在圖2中，令PZT片中心到位移觀測點的距離為l，懸臂梁上位移觀測點位移量為r，可得

令懸臂梁厚度為h，則PZT片位移量s為

1.2 PZT片穩(wěn)態(tài)電壓控制試驗

基于懸臂梁的運動特性，文中測算PZT位移的方法不利于描述PZT片在不同頻率下的遲滯特性，因此僅對PZT片靜態(tài)遲滯特性進行分析。在PZT片上施加電壓值為0 V、10 V、20 V、…、300 V、290 V、280 V、…、0 V，當懸臂梁穩(wěn)定時得到自由端位移r。將測得的數(shù)據(jù)r代入式（4），得到PZT片長度變化量與控制電壓關(guān)系如圖3所示。

圖3 PZT片電壓-位移曲線

由圖3可知，PZT片電壓-位移關(guān)系具有明顯的遲滯特性，該特性會影響PZT片的控制精度，因此需要對遲滯進行抑制以消除該特性的不利影響。

1.3 利用最小二乘法進行多項式擬合

為了對PZT片電壓-位移特性進行定量分析，需建立PZT片遲滯特性模型。針對PZT片靜態(tài)時壓電-位移數(shù)據(jù)分布特點，采取最小二乘法對離散點進行多項式擬合，獲得PZT片遲滯模型。假設(shè)n個點的擬合多項式為

每個點與曲線的距離平方和為

將式（6）兩邊同時對aj(j=0,1,2,…,k)求偏導(dǎo)，可得

化簡后可得

進一步寫成矩陣形式，可得

將式（9）中的矩陣化簡后可得

上式為X·A=Y的形式，通過xi、yi(i=1，2，…，n)求得系數(shù)矩陣A后，能夠進一步求得擬合函數(shù)。采取以上算法，可以得到當x最高階數(shù)為5時，兩條曲線擬合多項式分別為

f1(x)和f2(x)分別對應(yīng)圖4中電壓遞增時和電壓遞減時的擬合曲線。兩條擬合曲線帶來的誤差較小且易于計算，可作為PZT片的遲滯模型。

圖4 遲滯特性擬合曲線

2 逆遲滯補償算法

2.1 逆遲滯模型

在所建立的PZT片遲滯模型基礎(chǔ)上，進行逆模型計算。已知f1(x)和f2(x)在變化范圍內(nèi)單調(diào)遞增，因此存在逆遲滯模型。由于五次函數(shù)求逆過程較復(fù)雜，故采取擬合方法求得式（12），該式即為逆遲滯模型，如圖5所示。其輸入為PZT片位移，輸出為補償后電壓。

圖5 逆遲滯擬合曲線

2.2 逆遲滯補償

利用PZT片靜態(tài)遲滯特性對控制電壓進行逆遲滯補償，實現(xiàn)控制電壓與PZT片位移間的線性關(guān)系，完成懸臂梁振幅控制系統(tǒng)的閉環(huán)PID控制。具體流程如圖6所示。

圖6 控制流程

對控制系統(tǒng)輸入一個期望值，PID控制器會把激光位移傳感器采集所得數(shù)據(jù)與該期望值作比較，并根據(jù)誤差值來輸出控制電壓u。隨后通過逆遲滯補償模型對u進行逆遲滯補償，實現(xiàn)PZT片電壓-位移關(guān)系線性化。

逆遲滯補償前，對控制電壓u進行升降判斷。以u上升時的逆遲滯補償為例，當PID輸出電壓u時，期望位移為s，然而u對應(yīng)實際位移為sreal，如圖7所示。為了對PZT片壓電-位移進行線性化，使補償后的電壓與s對應(yīng)，可以算得s滿足

將式（13）代入式（12）中，可以求得補償后的電壓

圖7 逆遲滯補償原理

3 試驗過程與結(jié)果

3.1 硬件設(shè)置和軟件調(diào)試

PZT片粘貼在懸臂梁表面，在梁的寬度方向上居中，在梁的長度方向上靠近懸臂梁根部，如圖8所示。

將懸臂梁和激光位移傳感器固定在鋁制框架上，連接框架與激振器，試驗系統(tǒng)如圖9所示。

圖8 PZT片在懸臂梁上的位置

圖9 實驗設(shè)備

數(shù)據(jù)采集和處理流程如圖10所示，上位機通過數(shù)據(jù)輸出板卡使激振器產(chǎn)生正弦信號，從而使框架上的懸臂梁發(fā)生簡諧振動；激光位移傳感器采集懸臂梁自由端位移信號，通過數(shù)據(jù)采集板卡傳輸給上位機；上位機進行控制處理后，通過數(shù)據(jù)輸出板卡輸出控制電壓，壓電陶瓷驅(qū)動電源將控制電壓進行放大后作用在PZT片上，從而實現(xiàn)懸臂梁振動閉環(huán)控制。

圖10 數(shù)據(jù)采集和處理流程

試驗控制程序框圖如圖11所示。在程序框圖中，使用位移寄存器來判斷懸臂梁位移值是增大或減小，兩種情況分別對應(yīng)條件語句中不同的補償函數(shù)。此外，為消除干擾信號的影響，在對u進行判斷前采用巴特沃斯濾波器進行濾波，通過頻率為5 Hz～8 Hz。如果位移值增大，則以圖4中電壓遞增曲線對PID輸出進行逆遲滯補償；如果位移值下降，則以圖4中電壓遞減曲線對u進行補償。

圖11 labview軟件程序框圖

驅(qū)動電源的功率放大倍數(shù)為250倍，輸出電壓為0至300 V，因此PID控制器輸出電壓范圍設(shè)定為0至1.2 V。實驗中PZT片只能接受正電壓作用，產(chǎn)生縮短形變，因此懸臂梁與激光傳感器之間距離期望值選定為2 mm。試驗中得到PID控制器參數(shù)最優(yōu)為Kc=0.230，Ti=0.003，Td=0.000。

3.2 試驗結(jié)果與分析

試驗過程中測得壓電懸臂梁固有頻率為6.952 Hz。為了便于分析，設(shè)定激振器輸出信號為頻率6.952 Hz、幅值0.02 V的正弦波，在此基礎(chǔ)上驗證PID及PID逆遲滯補償?shù)恼駝又鲃涌刂菩Ч?。截取穩(wěn)定后的懸臂梁振幅曲線如圖12所示。

由圖12可知，在進行振動主動控制前，懸臂梁自由端振幅達5.219 mm，經(jīng)PID控制后幅值減小為1.557 mm，經(jīng)PID逆遲滯補償控制后振動幅值降為1.400 mm，比補償前振動控制效果提高了10.083%。因此，逆遲滯補償PID控制器優(yōu)于傳統(tǒng)PID控制器對壓電懸臂梁的振動主動控制效果。

圖12 不同工況下的懸臂梁振幅

4 結(jié)語

本文對PZT片遲滯特性進行了研究，并利用PZT片對懸臂梁進行振動主動控制試驗研究，得到如下結(jié)論：

(1)提出一種根據(jù)PZT片與懸臂梁形變關(guān)系來測算PZT片電壓-位移關(guān)系的方法，并通過試驗驗證了PZT片的遲滯特性；

(2)應(yīng)用最小二乘法對PZT片遲滯特性進行建模并求得逆遲滯補償模型，并用于對PZT片控制電壓進行補償；

(3)基于最小二乘逆遲滯補償方法，開展了壓電懸臂梁振動主動控制試驗研究，結(jié)果表明：經(jīng)過PID逆遲滯補償控制后的懸臂梁振幅比傳統(tǒng)PID控制振幅減小了10.083%，本文所提補償方法能有效增強壓電懸臂梁振動控制效果，對于壓電懸臂梁振動主動控制具有重要的參考價值。

參考文獻：

[1]楊曉明，李宗津.新型水泥基壓電傳感器的基本性能研究[J]. 傳感技術(shù)學(xué)報，2012，25(3)：349-354.

[2]張婷.基于遲滯特性與時滯現(xiàn)象的智能懸臂梁自適應(yīng)振動控制研究[D].上海：上海交通大學(xué)，2014.

[3]GOLDFARB M,CELANOVIC N.Modeling piezoelectric stack actuators forcontrolofmicromanipulation[J].Control Systems IEEE,1997,17(3):69-79.

[4]王代華，嚴松林，朱煒.基于Bouc-Wen模型的壓電執(zhí)行器的前饋線性化控制器[J].儀器儀表學(xué)報，2015，36(7)：1514-1521.

[5]李巍.壓電作動器遲滯非線性建模與補償控制研究[D].武漢：華中科技大學(xué)，2013.

[6]TARAZAGA PABLO,INMAN D,WILKIE W.Control of a space rigidizable-inflatable boom using macro-fiber composite actuators[J].JournalofVibration and Control,2007,13(7):935-950.

[7]黃建.基于壓電纖維復(fù)合材料的可扭轉(zhuǎn)機翼結(jié)構(gòu)研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2013.

[8]孫龍飛，李維嘉，吳耀中.壓電陶瓷在錐形圓筒振動主動控制中的應(yīng)用[J].噪聲與振動控制，2016，36(4)：188-192.