行波型桿式超聲電機定子的參數(shù)化 有限元法優(yōu)化設計

1、行波型桿式超聲電機定子的參數(shù)化有限元法優(yōu)化設計張健滔,朱華,趙淳生(南京航空航天大學精密驅(qū)動研究所,江蘇省 南京 210016摘要 :利用參數(shù)化有限元優(yōu)化方法,對行波型桿式超聲電機定子進行優(yōu)化設計。首先,在確定電機定子初始結(jié)構(gòu)的基礎上, 建立其參數(shù)化有限元模型。其次,對定子有限元模型進行模態(tài)分析,求解工作模態(tài)頻率對各結(jié)構(gòu)參數(shù)的靈敏度,選取靈敏度 高的結(jié)構(gòu)參數(shù)為設計變量,并以反映電機輸出性能的重要參數(shù)作為目標函數(shù)。同時,設計了定子結(jié)構(gòu)的優(yōu)化方案,采用了零 階優(yōu)化方法,對其結(jié)構(gòu)進行優(yōu)化設計。最后,根據(jù)優(yōu)化結(jié)果,制作了定子樣機。試驗表明:定子工作模態(tài)和端面質(zhì)點的振幅 都滿足了預期的設計要求, 試驗結(jié)

2、果與優(yōu)化設計結(jié)果相符。 研究表明, 利用該優(yōu)化設計方法能有效地縮短超聲電機設計周期。 關(guān)鍵詞:超聲電機;定子;有限元;模態(tài)分析;優(yōu)化設計中圖分類號:TM356 文獻標識碼:A引 言超聲電機具有低速大扭矩、 響應快、 扭矩 /質(zhì)量 比大、 結(jié)構(gòu)緊湊、 電磁兼容性和控制性能好等優(yōu)點, 已在精密儀器儀表、醫(yī)療器械、機器人、航空航天 以及新型武器裝備等領(lǐng)域得到廣泛的應用 1-4。 其中 圓環(huán)形行波電機的成本較高,加工工藝復雜,而行 波型桿式超聲電機(以下簡稱桿式電機由于采用 蘭杰文振子,結(jié)構(gòu)簡單,制造成本低,體積和重量 都大大減小,工藝要求相對低,易于實現(xiàn)產(chǎn)業(yè)化, 所以這種電機備受研究者的關(guān)注 5-9

3、。桿式電機定子利用兩個同頻、空間正交的彎曲 振動模態(tài)來合成端面質(zhì)點的橢圓運動,并在定子端 面產(chǎn)生一行波,推動轉(zhuǎn)子的轉(zhuǎn)動,實現(xiàn)定子的微幅 振動到轉(zhuǎn)子宏觀轉(zhuǎn)動的轉(zhuǎn)化 1,9。 因此, 在電機定子 設計時,使其具有合適的工作振型、工作頻率和足 夠大的振幅至關(guān)重要。通常設計者采用試湊的辦法 來滿足設計的要求,這樣不但設計效率低、耗費時 間多,而且有可能找不到滿意的結(jié)果。本文根據(jù)行波型桿式電機對其定子工作模態(tài) 和振幅的要求,探討了一種基于 ANSYS 有限元分 析的桿式電機定子優(yōu)化設計方法。首先,建立了定 子參數(shù)化有限元模型;其次,對其進行模態(tài)分析, 求解工作模態(tài)頻率對各結(jié)構(gòu)參數(shù)的靈敏度;根據(jù)靈 敏度分

4、析,選取設計變量;利用 ANSYS 參數(shù)化設 計語言(APDL 編程進行相應的優(yōu)化計算;根據(jù) 優(yōu)化結(jié)果制作了定子樣機,其試驗結(jié)果達到了預期 的設計要求。1 桿式電機定子優(yōu)化設計思路本文所研究的桿式電機由一個定子和二個轉(zhuǎn) 子組成。定子采用蘭杰文型振子結(jié)構(gòu),由預緊螺釘 把上、下振動塊和二組八片壓電陶瓷元件以及電極 片連接在一起。當對兩組壓電陶瓷元件分別施加正 弦、余弦交變電壓(其頻率接近定子一階彎曲共振 頻率時,將激勵出定子的兩個正交的彎曲振動模 態(tài),定子端面的質(zhì)點產(chǎn)生橢圓運動,通過定、轉(zhuǎn)子 間的摩擦力來驅(qū)動轉(zhuǎn)子轉(zhuǎn)動。 圖 1 桿式電機定子結(jié)構(gòu)及其參數(shù)基金項目:國家自然科學基金項目(5073500

5、2, 50575103在定子設計的初始階段,往往根據(jù)設計經(jīng)驗來 確定定子的基本結(jié)構(gòu)和初始結(jié)構(gòu)參數(shù)。圖 1所示是 電機的定子結(jié)構(gòu)以及相應的結(jié)構(gòu)參數(shù)。確定定子的 初始結(jié)構(gòu)后,建立定子參數(shù)化有限元模型,并計算 出工作模態(tài)頻率對各結(jié)構(gòu)參數(shù)的靈敏度。選取靈敏 度高的結(jié)構(gòu)參數(shù)為優(yōu)化設計變量,確定這些參數(shù)的 變化范圍, 以及確定優(yōu)化目標函數(shù)。 選擇優(yōu)化方法, 通過一系列的循環(huán)過程來優(yōu)化各個設計變量,直到 滿足所有設計要求為止,并輸出迭代過程和優(yōu)化結(jié) 果。定子優(yōu)化設計的思路可用圖 2表示。 圖 2 電機定子優(yōu)化設計思路2 有限元數(shù)值分析2.1 定子有限元模型桿式電機定子是一個機電耦合的動力學系統(tǒng), 同時具有力

6、場和電場的耦合作用。本文利用有限元 法建立定子的機電耦合模型。有限元法將連續(xù)的定 子結(jié)構(gòu)離散成有限個單元,并在每個單元中設定有 限個節(jié)點。對于壓電單元體,節(jié)點除了有位移自由 度外,還增加電勢自由度。假設單元內(nèi)的電場服從 線性分布,壓電單元的機電耦合作用采用線性壓電 本構(gòu)方程描述。定子系統(tǒng)的動力學微分方程可以寫 為 10,11mm m mm m me e fTme m ee e q M K K q F K K q Q +=&& (1式中:、 mm M mm K 、 me K 和 ee K 分別為定子的質(zhì)量陣、結(jié)構(gòu)剛度陣、 壓電耦合陣和介電傳導陣, f F 為轉(zhuǎn)子對定子驅(qū)動端面的廣

7、義作用力, 為壓電陶瓷片表 面電極上的電荷量, q Q m 為整個定子節(jié)點位移列向 量, 為電自由度向量。e q 式(1可以通過有限元法來求解。定子有限 元模型采用三維模型,充分考慮了壓電耦合、機械 耦合等效應。分析單元采用八節(jié)點六面體等參單 元,邊界條件為兩端自由。上、下振動塊采用的是 不銹鋼,預緊螺 釘和壓電陶瓷的材料分別是 11不銹鋼和 213Cr 89Cr Ni Ti 8PZT 。 不把壓電陶瓷結(jié)構(gòu)參數(shù)作為設計參數(shù), 用 1P 11個基本參數(shù)來確定定子的結(jié)構(gòu),初始 結(jié)構(gòu)尺寸如表 1所示。建立有定子限元模型(共 16089個節(jié)點, 12080個單元 ,如圖 3所示。11P 共 表 1 桿

8、式電機定子初始結(jié)構(gòu)尺寸參 數(shù) 1P 2P 3P 4P 5P 6P 數(shù)值 / mm 3.9 3.5 2.4 35 5 3 參 數(shù) 7P 8P 9P 10P 11P 數(shù)值 /mm 640°5.5 5.4 8 圖 3 桿式電機定子有限元模型2.2 定子結(jié)構(gòu)參數(shù)的靈敏度分析當定子的結(jié)構(gòu)形式和選用的材料確定之后,定子的振動模態(tài)就取決于定子的結(jié)構(gòu)參數(shù)。因此分析 工作模態(tài)頻率對各結(jié)構(gòu)參數(shù)的敏感程度(即靈敏度 分析 ,確定靈敏度高的結(jié)構(gòu)參數(shù),定義其為優(yōu)化 設計中的設計變量,這樣可以增強結(jié)構(gòu)修改的目的 性,提高定子優(yōu)化設計的效率。定子無阻尼自由振動方程為0mm m mm mM K &&

9、+= (2求解式 (2 廣義特征值得到 n 個特征值和對應的特征向量,也就可以得到定子的 階模態(tài)對應的 頻率和振型,其中 表示定子的自由度。n n 定子的第 階模態(tài)的固有振動方程為r2( mm rmm r K M =0 (3式中: r 為按質(zhì)量歸一化第 階的振型列陣; r r 為 第 階模態(tài)頻率。r 假設有 個結(jié)構(gòu)參數(shù) ( ,模態(tài)頻率 對結(jié)構(gòu)參數(shù) 的靈敏度表示為m i P 1, , i =m i P 212T mmmm r r r i r ii K M P P P r =(4采用差分近似得212T mmmm r r r i r i i K M P P r P =V V V V(5式中:為結(jié)構(gòu)參

10、數(shù) 的變化量; i P V i P mm K V 為定子剛度陣的變化量; mm M V 為定子質(zhì)量陣的變化量。給 一個攝動量 ,通過有限元求解得出相應的 i P i P V mm K V 和 mm M V ,并由式 (5 求出模態(tài)頻率對 的靈敏度。i P 根據(jù)前面建立的定子有限元模型,通過改變各 結(jié)構(gòu)參數(shù)(施加一攝動量 ,求解出工作模態(tài)頻率 對各結(jié)構(gòu)參數(shù)的靈敏度,如圖 4所示。 圖 4 工作模態(tài)頻率對結(jié)構(gòu)參數(shù)的靈敏度 從圖 4可以看出,工作模態(tài)頻率對結(jié)構(gòu)參數(shù) 、 、 、 、 、 的靈敏度較高。當改變這些 參數(shù)的尺寸時,可以很容易改變工作模態(tài)頻率。由 于受電機其他零部件尺寸的限制, 結(jié)構(gòu)參數(shù) 、

11、 變動范圍較小, 所以在優(yōu)化設計中把 、 固定不 變,只把參數(shù) 、 、 、 定義為設計變量。 1P 2P 3P 4P 9P 11P 9P 11P 9P 11P 1P 2P 3P 4P 3 有限元優(yōu)化設計 3.1 優(yōu)化目標函數(shù)的建立 根據(jù)行波型桿式電機運動機理,在預壓力不大 的情況下,如果認為定、轉(zhuǎn)子之間點接觸但無相對滑動,理論上電機的最大輸出轉(zhuǎn)速為 12,131060(min W f v r R=(6式中:為一階彎振橫向振幅, 0W f 為激勵頻率 (即 一階彎振共振頻率 , R 為定子的半徑。式(6可 以看出電機的轉(zhuǎn)速與一階彎振橫向振幅和共振頻 率成正比。為了使電機獲得更高的轉(zhuǎn)速和更大的轉(zhuǎn)

12、矩,必須提高定子一階彎振橫向振幅和共振頻率。 文獻 13把 稱為電機的轉(zhuǎn)速因子,它是反映電 機輸出性能的一個重要參數(shù)之一。在定子的結(jié)構(gòu)設 計中,可以通過在上、下振動塊上開槽的方法,提 高一階彎振橫向振幅,但這同時也會降低定子的一 階彎振共振頻率。所以這里存在一個優(yōu)化的問題, 如何合理地選取定子參數(shù)尺寸,以獲得比較大的轉(zhuǎn) 速因子。因此,設定轉(zhuǎn)速因子為優(yōu)化設計的目標函 數(shù)。電機定子優(yōu)化設計數(shù)學模型可表述為:0W f 01234min max (, , , . . (1, 2, 3, 4 i i i MAX W f P P P P S T P P P i =(式中:分別為變量 的極小值和極大值。 A

13、NSY 可將 7min max , i i P P i P 由于 S 程序總是最小化目標函數(shù), 因此, 問題轉(zhuǎn)化為求下式的最小值101x W f = (8這樣在優(yōu)化設計中 1x 被設定為目標函數(shù),其中 它是一種通 用的 數(shù)化有 限元 計變量參數(shù)的攝動。 有限元模 型, 與參考振型進行模態(tài)置信度 (M 1是 大于 0W f 的數(shù)值。 3.2 優(yōu)化算法在優(yōu)化過程中采用零階優(yōu)化方法, 方法,可以有效地處理絕大多數(shù)的工程問題。 其主要是通過對目標函數(shù)添加罰函數(shù),將約束問題 轉(zhuǎn)化為非約束優(yōu)化問題,再用曲線擬合來建立目標 函數(shù)和設計變量之間關(guān)系來實現(xiàn)逼近的。進行優(yōu)化分析之前,首先建立定子的參 模型,并進行

14、模態(tài)分析,提取所需的工作模態(tài) (一階彎振模態(tài)振型作為參考振型。優(yōu)化分析的 迭代過程如下。(1 進行設 (2 將攝動參數(shù)代入定子參數(shù)化 進行模態(tài)分析。(3 將計算振型 odal Assurance Criteria, MAC 計算,對一階 彎曲振動模態(tài)進行識別, 并提取其模態(tài)頻率。 MAC 計算可以用來判別待識別振型與參考振型之間的 相似度, MAC 值可以表示為 2( ( (T i w 1, , (i T Ti i w w MAC i =m = (9 式中:i 為使用有限元法計算所得的待 別的 振型, 識 定子 w 為定子一階彎曲振動模態(tài)的參考振型。 變化 度(4 進行諧響應分析。諧響應分析是

15、用來確定電機定子結(jié)構(gòu)在承受隨時間按正弦(簡諧規(guī)律 的載荷時的穩(wěn)態(tài)響應。進行諧響應分析時,將 對電機定子上的壓電陶瓷元件施加正弦、余弦兩相 交變電壓(其頻率為(3所提取的定子一階彎振 模態(tài)頻率 。諧響應分析的目的就是求解出在兩相 交變電壓的激勵下,定子端面質(zhì)點的橫向振幅。(5 計算轉(zhuǎn)速因子和目標函數(shù),滿足收斂精 11(x k (1 (x k x k +1,則停止迭代;否則返回(1評估 。 在優(yōu)化設計的過程中,程序執(zhí)行一系列的分析 修正的循環(huán)過程,重復直到滿足設計 要求 按照表 1中的定子初始結(jié)構(gòu)尺寸進行優(yōu)化計 范圍如表 2所示,優(yōu)化的收斂 誤差 為止。 3.3 優(yōu)化結(jié)果分析 算,設計變量的變化

16、取為 ×011變化范 -6。 表 2 設計變量的變化范圍 設計 變量 圍 1mm P / 3 ; 5 2P / mm P/圖 5 目標函數(shù)隨迭代步數(shù)的變化曲線圖 5給出優(yōu)化目標函數(shù)隨迭代步數(shù)的變化曲 線。由計算結(jié)果可知,經(jīng) 22次迭代后求解收斂。 表 3是設計變量以及目標函數(shù)的優(yōu)化前后變化情 況,頻率、振幅和轉(zhuǎn)速因子的變化情況也在表 3中 列出。 從表 3可見, 優(yōu)化后的轉(zhuǎn)速因子為 0.075280, 比優(yōu)化前增加了 44.8%,即理論上電機的最大輸出 轉(zhuǎn)速提高了 44.8%,說明優(yōu)化的效果很明顯。表 3 優(yōu)化前后設計變量及目標函數(shù)變化 類 別參 數(shù) 初 值 優(yōu) 化 后 1P / m

17、 m 3. 900 3. 0032P / m m 3. 500 3. 9073P / m m 2. 400 2. 998設 計 變 量3 ; 434 ; 37.53mm P /1 ; 34mm4P / m m35. 000 34. 004頻 率 / k H z39. 46236. 691 振 幅 / µm1. 317 2. 051轉(zhuǎn) 速 因 子0. 051973 0. 075280目 標 函 數(shù)0. 9480270. 9247204 電機試驗分析 根據(jù)優(yōu)化設計結(jié)果, 制作了樣機。 利用 PSV-300和 PCV-Z-040-F 多普勒激光測振儀測試了定子的振 動模態(tài),其幅頻響應曲線和

18、一階彎振振型如圖 6所 示。圖 6(a 、 (b分別為定子兩個空間正交的一 (a A 相 一 階 彎 振 振 型(b B 相 一 階 彎 振 振型(c 定子頻響曲線 圖 6 桿式電機定子頻響曲線與實測振型階彎振模態(tài)振型。從圖 6(c 可以看出,這兩個彎振模態(tài)頻率分別為 35.73kHz 和 35.67kHz ,這個結(jié)果 與優(yōu)化設計的頻率值 36.691kHz 很接近。 同時, 還利 用了多普勒激光測振儀對電機進行了定頻測試,測 得在 80V p-p 激勵電壓下,定子端面質(zhì)點一階彎振橫 向振幅為 2µm,與優(yōu)化結(jié)果振幅為 2.051µm很接近。圖 7 桿式電機實物照片圖 7所

19、示為裝配好的桿式電機,利用電機機械 特性自動測試系統(tǒng)對電機的輸出特性進行了測試, 測試結(jié)果如圖 8所示。電機通過專用的超聲電機驅(qū) 動器供電,驅(qū)動器的輸入電壓為直流 15V ,定、轉(zhuǎn) 子之間的預壓力為 9kgf 。 從圖 8可以看出, 電機轉(zhuǎn)速 隨輸出轉(zhuǎn)矩的增加而線性減小。電機的空載轉(zhuǎn)速為2061r min , 堵轉(zhuǎn)力矩為 0.258。 當電機轉(zhuǎn)速為 99N m 1r min ,輸出轉(zhuǎn)矩為 0.126時,電機的 輸出功率最大為 1.3。電機的輸出性能較文獻 5的電機有較大的提高。N m W 100 90 80 70 60 50 40 30 20 10 0 2.0 1.75 1.5 1.25 1.

20、0 0.75 0.5 0.25 0 0 0.05 要求的最優(yōu)解，并能大大的減少設計工作量，其可 為超聲電機的優(yōu)化設計提供一條新的途徑。 N Pout 參考文獻： 1 2 趙淳生超聲電機技術(shù)與應用M北京：科學出版社，2007. S. Ueha, Y. Tomikawa. Ultrasonic Motors: Theory and Applications M. Oxford, U.K.: Clarendon, 1993. Kenji Uchino. Piezoelectric Actuators and Ultrasonic Motors M. Norwell, MA: Kluwer, 1997

21、. Chunsheng Zhao, Jiantao Zhang, Jianhui Zhang, Jiamei Jin. Development and application prospects of piezoelectric precision diving technology J. Frontiers of Mechanical Engineering in China, 2008, 3(2: 119-132. M. Kurosawa, K. Nakamura, et al. An ultrasonic motor using bending vibrations of a short

22、 cylinder J. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 1989, 36: 517-521. Takashi Maeno. Recent progress of ultrasonic motor in Japan C. The First International Workshop on Ultrasonic Motors and Actuators. Yokohama, Japan, 2005: 15-17. Ichiro Okumura. A designing metho

23、d of a bar-type ultrasonic motor for autofocus lenses C. International Symposium on Theory of Machines and Mechanism, Nagoya, Japan, 1992: 836-841. Shuxiang Dong, Shuxin Wang, et al. A miniature piezoelectric ultrasonic motor based on circular bending vibration mode J. IEEE Transactions on Mechatron

24、ics, 2000, 5(4: 325-330. 朱華，陳超，趙淳生. 一種微型柱體超聲電機 J. 中國電機工 程學報，2006，26（12） ：128-133. H Allik, Thomas J R Hughes. Finite element method for piezoelectric vibration J. International Journal for Numerical Methods in Engineering, 1970, 2: 151-157. 朱華. 行波型桿式超聲電機的關(guān)鍵技術(shù)研究 D. 南京：南京航 空航天大學，2006. Pin Lu, Kwok Hon

25、g Lee, Siak Piang Lim, Wu Zhong Lin. A kinematic analysis of cylindrical ultrasonic micromotors J. Sensors and Actuators, 2001, A 87:194-197. Zhirong Li, Chunsheng Zhao, Weiqing Huang. Several key issues in developing of cylinder type 3-DOF ultrasonic motor J. Sensors and Actuators, 2007, A 136: 704

26、-709. 0.1 0.15 3 4 轉(zhuǎn)矩 T 5 圖 8 桿式電機機械特性試驗結(jié) 果曲線 6 5 結(jié)論 （1）定子的實測振動模態(tài)和端面質(zhì)點橫向振 幅，與優(yōu)化計算結(jié)果相符，表明了本文所建立的電 機參數(shù)化有限元模型和優(yōu)化設計方法是正確的。 （2）通過優(yōu)化設計，定子的轉(zhuǎn)速因子得到了 較大的提高，定子工作模態(tài)和端面質(zhì)點的振幅都滿 足了電機設計的要求。 （3）優(yōu)化后的電機，其輸出性能得到了有效 地提高。 與以往的超聲電機設計方法相比，參數(shù)化有限 元優(yōu)化設計方法可以快速地搜索到滿足電機設計 7 8 9 10 11 12 13 FEM-based parametric optimum design of

27、a traveling-wave type bar-like ultrasonic motor stator Zhang Jian-tao, Zhu Hua, Zhao Chun-sheng (Precision Driving Laboratory, Nanjing University of Aeronautics & Astronautics, Nanjing 210016, China Abstract: This paper presents the optimum design of a bar-like ultrasonic motor (USM stator using finite element (FE parametric optimum method.