曲 磊， 王建國， 丁根芳， 覃 艷

(合肥工業(yè)大學(xué)土木與水利工程學(xué)院，安徽合肥 230009)

0 引 言

壓電層合結(jié)構(gòu)作為傳感器與致動器，在航天、精密控制與測量及微機(jī)電系統(tǒng)中廣為采用[1，2]。壓電層合結(jié)構(gòu)已成為近20年來國內(nèi)外學(xué)者研究的熱點，近年來用于壓電多層結(jié)構(gòu)分析的有限元數(shù)值方法得到了飛速發(fā)展。

文獻(xiàn)[3，4]最早開展了這方面的研究工作;文獻(xiàn)[5]總結(jié)了壓電三明治層合結(jié)構(gòu)的有限元理論;文獻(xiàn)[6]基于離散層理論提出了板殼的有限元分析模型;文獻(xiàn)[7]提出了一種能用于拉伸致動與剪切致動3層壓電耦合結(jié)構(gòu)分析的通用有限元模型;文獻(xiàn)[8]運用雙線性四節(jié)點板單元研究了熱壓電層合板;文獻(xiàn)[9]采用改進(jìn)的Reissner-Mindlin模型用于壓電層合板的有限元分析;文獻(xiàn)[10，11]在壓電有限元方面也做了很多工作;文獻(xiàn)[12]基于一階剪切變形理論，采用一種新的具有多自由度的四結(jié)點四邊形有限單元CTMQE對壓電復(fù)合板進(jìn)行了研究;文獻(xiàn)[13]基于狀態(tài)空間法系統(tǒng)地研究了多層壓電介質(zhì)空間軸對稱與非軸對稱問題的解。

本文在文獻(xiàn)[14，15]的基礎(chǔ)上，建立了含有壓電傳感器與致動器的多層耦合結(jié)構(gòu)的有限元方程，采用二維四結(jié)點四邊形壓電耦合單元，借助ANSYS/APDL語言，編制了適用任意組壓電片通過面層材料粘貼于智能結(jié)構(gòu)表面的力-電多場通用有限元分析程序(M PFEMP)。

本文首先運用該程序驗算了疊層梁[16]，與文獻(xiàn)[16]比較證明了方法的可行性與有效性;計算了由壓電片、黏結(jié)層和主體材料組成的5層復(fù)合梁，算例中提供了1組壓電片和2組壓電片在不同邊界條件下的計算結(jié)果;討論了壓電片長度和位置變化對梁的變形影響，還提供了閉環(huán)情況下前5階自振頻率和C-F情況下的振型圖。

本文結(jié)論對智能結(jié)構(gòu)的設(shè)計具有一定的參考價值，本文方法可以作為復(fù)雜不規(guī)則壓電復(fù)合結(jié)構(gòu)有限元分析的計算基礎(chǔ)。

1 線彈性壓電體有限元方程

壓電層合結(jié)構(gòu)作如下假定:①結(jié)構(gòu)各層之間理想黏結(jié);②遵循線性彈性壓電體理論;③各個壓電片之間無相互作用的電磁場干擾;④壓電片上下表面完全覆蓋電極。線性壓電結(jié)構(gòu)有限元方程[14]為:

由方程(1)的第2式解出φi=Kφφ-1(g i-Kφu u i)，代入方程(1)的第1式可得:

方程(2)是以結(jié)點位移表示的線性壓電結(jié)構(gòu)的有限元法動力方程。

2 MPFEMP程序驗算

力-電有限元分析程序(MPFEMP)主要采用參數(shù)化設(shè)計語言APDL，利用*DO、*ENDDO循環(huán)語句、*IF、*ELSE、*ENDIF等判斷語句和*DIM矩陣定義語句來實現(xiàn)在壓電片上模擬任意數(shù)量的電極，在任意電極上加載任意變化的電壓載荷。為驗證程序有效性，本文先對中間層全為壓電層的疊層梁[16]展開對比驗證分析。幾何參數(shù)、材料常數(shù)、載荷和邊界條件均同文獻(xiàn)[16]，本文計算所得梁端撓度如圖1所示。

由圖1可看出，本文結(jié)果與文獻(xiàn)[16]解析解(1.19E-7 m)吻合良好，誤差率為1.3%，證明了本文計算方法的有效性。同時，本文的退化程序在文獻(xiàn)[14]中也得到了驗證。

圖1 3層壓電梁的變形曲線

3 n組壓電片粘貼的5層壓電梁

5層壓電梁如圖2所示，中間層主體材料為鋁，黏粘層為剛性泡沫，壓電層粘貼在黏粘面層上。壓電材料極化方向平行于x軸。

圖2 5層壓電耦合梁示意圖

本文作如下假定:①電極相對于壓電材料和鋁層不計質(zhì)量及剛度;②)壓電層、電極、剛性泡沫和鋁層的黏結(jié)是理想的，不存在相對位移;③壓電材料不存在相互的電磁效應(yīng);④材料沿Z向變形均勻。梁總長L為100 mm，鋁層h3為8 mm，黏粘面層為0.25mm，上下壓電層等厚均為2mm，第i組壓電片的長度ai和位置bi表達(dá)式為:

鋁、剛性泡沫的材料參數(shù)見文獻(xiàn)[7]，壓電材料參數(shù)見文獻(xiàn)[17]。在每個壓電片上施加100 V電壓荷載。分析中考慮梁兩端固支(C-C)、兩端均鉸結(jié)(H-H)、一端固支一端自由(C-F)、一端固支一端鉸結(jié)(C-H)4種邊界條件。有限元分析采用四結(jié)點四邊形單元(PLANE13)，5層梁結(jié)構(gòu)劃分成550個網(wǎng)格，兼顧了計算效率和計算精度。

3.1 粘貼1組壓電片對梁的變形分析

壓電片長度a=0.01m，位置b變化時，梁的變形情況如圖3所示。

壓電片位置b=0.05m，長度a變化時，梁的變形情況如圖4所示。

圖3 壓電片b變化時梁的變形情況

圖4 壓電片a變化時梁的變形情況

3.2 粘貼2組壓電片對梁的變形分析

對于粘貼2組壓電片對的梁，當(dāng)壓電片長度a1=a2=0.01 m，位置b1、b2變化時，梁的變形情況如圖5所示。2組壓電片的位置b1=0.025m，b2=0.075 m，長度a1、a2變化時，梁的變形情況如圖6所示。

圖5 壓電片b1、b2變化時梁的變形情況

圖6 壓電片a1、a2變化時梁的變形情況

4 壓電梁的自由振動

研究結(jié)構(gòu)的自振頻率與振型是結(jié)構(gòu)形狀控制和主動控制的基礎(chǔ)。本文采用APDL語言編制層合梁的自由振動分析程序，壓電單元采用二維四結(jié)點應(yīng)力耦合單元，普通材料采用平面等參結(jié)構(gòu)單元，計算:①單壓電片對粘貼于結(jié)構(gòu)表面的自由振動，a=0.01 m，b=0.05 m;②2組壓電片對粘貼于結(jié)構(gòu)表面的自由振動，a1=a2= 0.01 m，b1=0.025m，b2=0.075 m。材料參數(shù)同上，閉環(huán)前5階自振頻率見表1所列。

粘貼1組壓電片對和粘貼2組壓電片對懸臂梁的前3階振型圖如圖7所示。

由圖7可知，2種不同的粘貼方式，變形趨勢基本相同，Benjeddou[7]的梁幾何尺寸與本文不同，但是對比發(fā)現(xiàn)，梁整體的變形趨勢基本相同。

表1 單壓電片對及2組壓電片對黏結(jié)梁的前5階自振頻率 H z

圖7 單壓電片對與2組壓電片對粘貼梁的前3階閉環(huán)振型圖

5 結(jié) 論

本文利用APDL語言編制的有限元程序?qū)Ρ砻嬲迟N壓電片的層合梁進(jìn)行了計算分析，考察了壓電層位置和長度對梁變形的影響，分析了其自振頻率和C-F情況下的振型圖。結(jié)論如下:

(1)壓電片相同的長度，不同的黏結(jié)位置和不同的邊界條件，所產(chǎn)生的梁的最大位移值均不相同。H-H、C-H和C-C 3種情況下，梁的最大位移發(fā)生在壓電片粘貼的中心處，隨著壓電片位置的變化，無論對于黏結(jié)1組壓電片和2組壓電片的梁，最大變形均受影響較大。所以，壓電片在梁中的位置存在著一個最優(yōu)值。同時在壓電片長度固定的情況下，壓電片位置的變化對于懸臂梁的最大變形影響不大。

(2)壓電片長度越大，梁的最大位移越大。無論粘貼1組壓電片還是2組壓電片的層合梁，壓電片長度增加，最大變形均增大，其中C-F情況下，最大位移發(fā)生在梁的自由端。而在其它的邊界條件下，配置1組壓電片的結(jié)構(gòu)大部分發(fā)生在梁中間處;而配置了2組壓電片的結(jié)構(gòu)最終都在壓電片所在的區(qū)域發(fā)生最大變形，形成非常有趣的雙波峰現(xiàn)象。所以，利用多組壓電片粘貼于結(jié)構(gòu)表面，可以按照需要設(shè)計結(jié)構(gòu)的變形。

(3)對比埋置1組壓電片和埋置2組壓電片的層合結(jié)構(gòu)，發(fā)現(xiàn)壓電片數(shù)量的增加對梁的最大位移影響不大，但是對于梁任何一處獲得足夠理想的致動效果影響較大。

(4)梁固支后自振頻率明顯增大。增加壓電片的數(shù)量后，結(jié)構(gòu)的自振頻率略有增大，但是總體振型趨勢變化不大，說明表面黏結(jié)形式的壓電梁，壓電片的數(shù)量和位置對梁的自由振動影響不大。

[1] Dem eis R B.Piezoelectric dam ping for space structures[J]. Aerospace AM，1985，23(9):32-33.

[2] Liao W H，W ang D H，H uang S L.W irelessmonitoring of cab le tension of cable-stayed b ridges using PVDF piezoelectric films[J].J Intel Mat Syst Str，2001，12(5):331-339.

[3] A llik H，Hughes T J.Finite element method for piezoelectric vib ration[J].Int J Num er Meth Eng，1970，2: 151-157.

[4] C raw ley E F，Deluis J.Use of piezoelectric actuators as elements of intelligent structu res[J].AIAA Journal，1987，25(10):1373-1385.

[5] H a K H.Finite-element analy sis of sandw ich plates:an overview[J].Compu t Struct，1990，37(4):397-403.

[6] Tzou H S，Ye R.Analysis of piezoelectric structu res with lam inated piezoelectric triangle shell elements[J].AA IA Journal，1996，34:110-115.

[7] Benjeddou A，Trindade M A，Ohayon R.A unified beam finite elem en tmodel for extension and shear piezoelectric actuationm echanisms[J].J Intel M at Sy st Str，1997，8(12): 1012-1025.

[8] Lee H J，Saravanos D A.Generalized finite element formulation for smart multilayered thermal piezoelectric composite plates[J].Int J Solids Structures，1997，34(26): 3355-3371.

[9] Carrera E.An improved Reissner-M indlin-type m odel for the electromechanical analysis of m ultilayered plates including piezo-layers[J].J Intel Mat Syst Str，1997，8(3): 232-248.

[10] Tian X G，Shen Y P.Finite element analysis of thermoelastic behavior of piezoelectric structuresunder finite deformation[J].Acta M ech Solida Sin，2002，15(4):312-322.

[11] 丁皓江，池毓蔚，國鳳林，等.壓電材料軸對稱有限元分析[J].計算力學(xué)學(xué)報，2000，17(1):1-7.

[12] Cen Song，Soh A K，Long Y Q，et al.A new 4-node quad rilateral FE model w ith variable electrical deg ree of freedom for the analysis of piezoelectric lam inated composite p lates [J].Composite Stru ctures，2002，58(4):583-599.

[13] W ang JG，F(xiàn)ang SS，Chen L F.The state vector methods for space axisymmetric p rob lem s in multilayered piezoelectric media[J].Int J Solids Structu res，2002，39: 3959-3970.

[14] 丁根芳，王建國，覃 艷.自感知壓電層合梁的數(shù)值分析與數(shù)值模擬[J].工程力學(xué)，2006，23(2):131-136.

[15] 王建國，覃 艷，丁根芳.單點配置壓電層合結(jié)構(gòu)的優(yōu)化設(shè)計[J].合肥工業(yè)大學(xué)學(xué)報:自然科學(xué)版，2007，30(5): 594-597.

[16] Zhang X D，Sun C T.Formu lation of an adap tive sandw ich beam[J].Smart Mater Struct，1996，5(6):814-823.

[17] Sm its J G，Dalke S I，Coony T K.The constituent equationsof piezoelectric bim orphs[J].Sensor and Actuators A:Physical，1991，28:41-61.