于棟 杜斌 嚴有琪

摘 要：該文運用ANSYS對導(dǎo)波檢測時所使用的矩形pzt-5H型壓電片陶瓷振子進行了諧響應(yīng)分析和瞬態(tài)分析，瞬態(tài)分析時加載漢寧窗調(diào)制的音頻信號。首先通過諧響應(yīng)分析確定了壓電陶瓷振子橫向效應(yīng)的諧振頻率，然后通過瞬態(tài)分析分析了壓電陶瓷振子的瞬態(tài)響應(yīng)能力。結(jié)果表明，漢寧窗調(diào)制的音頻信號能夠驅(qū)動壓電陶瓷振子振動，但表現(xiàn)出遲滯現(xiàn)象。

關(guān)鍵字：壓電陶瓷 ANSYS 瞬態(tài)分析 漢寧窗

中圖分類號：TM359 文獻標識碼：A 文章編號：1672-3791（2015）08（a）-0035-03

Transverse Effect Analysis on Piezoelectric Ceramic Vibrator Based on ANSYS

Dong Yu1 Bin Du1，2 Youqi Yan1

（ 1. Special Equipment Safety Supervision Inspection Institute of Jiangsu Province，Zhenjiang 212013，China；2. School of Mechanical Engineering，Jiangsu University，Zhenjiang，Jiang Province，212013，China）

Abstract：In this paper， harmonic and transient response of Rectangular PZT-5H Piezoelectric ceramic vibrator， which were used to excite guided wave， were researched based on ANSYS， hanning window modulated audio signal were loaded during transient analysis.First，the resonant frequency of the piezoelectric vibrator transverse effect were researched by harmonic response analysis. Then transient response ability of piezoelectric ceramic vibrator were analyzed by utilizing transient analysis. The results show that the hanning window modulated audio signal were able to drive the piezoelectric ceramic vibration， but show the hysteresis phenomenon.

Key Words：Piezoelectric ceramics； ANSYS； Transient analysis； Hanning window

超聲導(dǎo)波在長距離檢測例如管道檢測中有較明顯的優(yōu)勢，在結(jié)構(gòu)的一點激勵超聲導(dǎo)波，由于導(dǎo)波衰減小，傳播的距離很遠，因此一次檢測的距離較大；同時由于超聲導(dǎo)波在被檢工件整個壁厚范圍內(nèi)存在，因此是對工件進行的整體檢測。

在激勵縱向模態(tài)導(dǎo)波時，利用了壓電陶瓷振子的橫向效應(yīng)，即長度方向伸縮振動。試驗前需要對壓電振子進行振動分析，傳統(tǒng)的解決方法是建立振子的自由振動方程，用解析方法解答。然而，這些方法數(shù)據(jù)處理起來比較繁瑣，而且處理的結(jié)果不夠直觀[1]。借助ANSYS可以直觀的對壓電振子進行振動分析，ANSYS中ANSYS Multiphysics或者ANSYS mechanical軟件包提供了壓電分析，可以進行穩(wěn)態(tài)、模態(tài)、預(yù)應(yīng)力模態(tài)、諧波、預(yù)應(yīng)力諧波和瞬態(tài)分析。

1 壓電陶瓷常用振動模式

當(dāng)施加不同方向的極化方向和電場方向時，壓電振子會表現(xiàn)出不同的振動模式，如圖1所示為四種類型的壓電振動模式[2]：圖1（a）、（b）為厚度方向極化，圖1（c）、（d）為寬度方向極化；全部為厚度方向施加電場。壓電陶瓷常用的振動模式分別為：

（a） 垂直于電軸方向伸縮（長度伸縮），簡稱為LE；

（b） 平行于電軸方向伸縮（厚度伸縮），簡稱為TE；

（c） 垂直于電軸的面內(nèi)切變（面切變，又稱輪廓切變），簡稱為FS；

（d） 平行于電軸的面內(nèi)切變（厚度切變），簡稱為TS。

根據(jù)彈性波傳播方向與電軸方向之間的關(guān)系，可將壓電陶瓷振動方式分為兩種：傳播方向垂直于電軸方向，稱為橫向效應(yīng)包括LE模式和FS模式；傳播方向平行于電軸方向，稱為縱向效應(yīng)包括TE模式和TS模式。根據(jù)縱向?qū)Р詈徒邮盏睦碚?，激勵和接收縱向模態(tài)導(dǎo)波時所使用的是壓電陶瓷振子長度方向伸縮模式，即圖1（a）所示LE模式。

2 ANSYS模型建立

2.1 壓電陶瓷的輸入?yún)?shù)

壓電分析采用第二類壓電方程：

（1）

其中，T，S 表示應(yīng)力和應(yīng)變矢量；D，E 分別表示電位移矢量和電場矢量；e，cE，S 分別表示壓電應(yīng)力系數(shù)、壓電彈性系數(shù)和壓電介電系數(shù)。

壓電陶瓷片選用pzt-5A型進行分析，其模型的輸入?yún)?shù)分別為[3]：相對介電常數(shù)矩陣[]=，壓電應(yīng)變系數(shù)矩[d]=（10-12m/v），壓電剛度矩陣[c]=（1010N/m2）。

分析時選用solid226單元，keyopt（1）=1 001。solid226為20節(jié)點6面體單元，采用高階單元（中間節(jié)點形式二階單元），與solid5相比，分析精度提高，改善了分析準確性。

2.2 壓電陶瓷模型建立

在ANSYS前處理模塊，采用自上而下的建模方式，直接生成長方形壓電陶瓷振子模型。壓電陶瓷振子的尺寸為20 mmX4 mmX1 mm。采用solid226單元進行并且運用掃略的方式進行網(wǎng)格劃分，控制網(wǎng)格單元尺寸為0.5mm，最終生成的模型如圖2所示，模型中包含640個單元，3 909個節(jié)點。

3 諧響應(yīng)分析

3.1 諧響應(yīng)分析簡介

諧響應(yīng)分析是用于確定線性結(jié)構(gòu)在承受隨時間按正弦（簡諧）規(guī)律變化的載荷時的穩(wěn)態(tài)響應(yīng)的一種技術(shù)。分析的目的是計算出結(jié)構(gòu)在幾種頻率下的響應(yīng)并得到一些響應(yīng)值（通常是位移）對頻率的曲線。諧響應(yīng)分析可采用三種方法：完全法（Full）、縮減法（Reduced）、模態(tài)疊加法（Mode？Superposition）。其中完全法使用簡單，并且使用完全矩陣，不涉及質(zhì)量矩陣的近似，求解準確，該文將使用完全法進行壓電陶瓷的諧響應(yīng)分析。

3.2 邊界條件與加載

該文主要分析壓電陶瓷在自由條件下的橫向效應(yīng)，因此不對模型施加位移約束。為了模擬實際中的電極面，分別耦合z軸上下表面電壓，耦合電壓后的表面在以后的分析過程中始終保持電勢相同。諧響應(yīng)分析時需要設(shè)置載荷的幅值和相位，因此加載條件為：z軸下表面接地，上表面加載電壓幅值為20 V，相位為0，即正弦信號。

3.3 結(jié)果分析

提取1242節(jié)點（0，0.002，0.001）x軸方向的位移Ux隨頻率變化的曲線圖，如圖3所示，放大70 KHz至85 KHz的位移圖，可以觀察出1242節(jié)點在77 KHz時位移最大，即在77 KHz時壓電陶瓷產(chǎn)生諧振。如圖4所示，模型77KHz下的變形圖，實體為變形前模型，虛線為變形后模型?？梢缘贸鲈谕饧?7 KHz正弦函數(shù)激勵時，模型主要振動方式為長度方向伸縮振動，即LE模式。

通過諧響應(yīng)分析可以分析模型在外加不同頻率的正弦變化的載荷時的振動模式。該文通過分析20 ～200 KHz之間模型的諧響應(yīng)的出了模型在77KHz下進行長度方向伸縮諧振。

4 瞬態(tài)分析

4.1 瞬態(tài)分析簡介

瞬態(tài)分析又稱時間歷程分析，是用于確定承受任意的隨時間變化的載荷的結(jié)構(gòu)動力學(xué)響應(yīng)的一種方法。對于壓電瞬態(tài)分析，為使求解穩(wěn)定應(yīng)用TINTP命令設(shè)置積分參數(shù)ALPHA=0.25、DELTA=0.5和THETA=0.5[4]。

4.2 加載設(shè)置

導(dǎo)波激勵時通常使用經(jīng)過漢寧窗調(diào)制的音頻信號。與正弦音頻信號相比，經(jīng)過漢寧窗調(diào)制的音頻信號的能量高并且能保持較窄的頻帶，這有利于長距離檢測和激勵單一模態(tài)的導(dǎo)波[5]。該文中加載的信號為經(jīng)過漢寧窗調(diào)制的10周期音頻信號，加載信號如圖5所示。加載條件為z軸下表面接地，上表面加載圖5所示的信號。

4.3 結(jié)果分析

為研究模型在加載圖5所示信號下的響應(yīng)，提取了節(jié)點1242在x軸方向隨時間變化的位移圖，如圖6所示。由圖6可知，在加載調(diào)制的音頻信號后，壓電陶瓷按照調(diào)制音頻信號的變化規(guī)律變化，說明調(diào)制的音頻信號能夠驅(qū)動壓電陶瓷片發(fā)生振動。對比圖5和圖6可知，在加載10個周期的音頻信號后，在對應(yīng)的時間內(nèi)壓電陶瓷振子在x軸方向的振動并不是完整的10個周期，這是由于壓電陶瓷的遲滯性所致的[6]。

5 結(jié)語

該文運用ANSYS軟件對矩形壓電陶瓷振子的橫向效應(yīng)進行了分析，分別進行了諧響應(yīng)分析和瞬態(tài)分析。通過諧響應(yīng)分析確定了壓電陶瓷振子橫向效應(yīng)的諧振頻率。瞬態(tài)分析時加載了頻率和諧振頻率一致并經(jīng)過漢寧窗調(diào)制的音頻信號，結(jié)果表明經(jīng)漢寧窗調(diào)制的音頻信號能夠驅(qū)動壓電陶瓷振子振動，但是表現(xiàn)出一定的遲滯性。該文沒有研究壓電陶瓷遲滯性對導(dǎo)波激勵的影響，還有待進一步探究

參考文獻

[1] 羅平.基于ANSYS 模態(tài)分析的壓電陶瓷參數(shù)數(shù)據(jù)的轉(zhuǎn)換及其輸入[J].廣東技術(shù)師范學(xué)院學(xué)報，2009（9）：11-15.

[2] 韓寶坤，孫超，呂航.基于ANSYS 的純壓電振子瞬態(tài)特性分析[J].微特電機.2012，40（4）：35-37.

[3] 魏勝.基于ANSYS 的壓電換能器有效機電耦合系數(shù)的仿真分析[J].深圳職業(yè)技術(shù)學(xué)院學(xué)報，2012，11（3）：7-12.

[4] 呂航，范明慶，孫超，等.基于ANSYS 的矩形壓電振子瞬態(tài)特性分析[J].微特電機，2012，40（3）：34-37.

[5] Jennifer E. Michaels， Sang Jun Lee， Anthony J. Croxford， et al. Chirp excitation of ultrasonic guided waves[J].Ultrasonics，2013，53（1）：265-270.

[6] 王社良，曹文華，樊禹江.基于ANSYS的新型壓電陶瓷管驅(qū)動器的性能分析與研究[J].材料導(dǎo)報B：研究篇，2012，26（12）：154-157.