【正文】 AREA39。VOLU39。 例如：雙疊片換能器在水中的諧響應(yīng)分析命令流文件如下：/BATCH ! /COM,ANSYS RELEASE UP20041104 23:25:18 06/06/2006 /TITLE,head!material1金屬材料mat1dens=2790mat1yang=mat1person=!material3壓電體材料PZT4mat2dens=7550c11=c12=c13=c22=c23=c33=c44=c55=c66=e13=e23=e33=e52=e61=fx=fy=fz=!material3近場水mat3dens=1000 ！水密度sonc3=1480 ！水中聲速mu3=0 ！吸收系數(shù)!!material4遠場水mat4dens=1000sonc4=1480mu4=0!!structura5邊界水mat5dens=1000sonc5=1480mu5=1!! 結(jié)構(gòu)參數(shù)h1=h2=h3=h4=h5=h6=r1=r2=r3=r4=r5=r6=h7=r6*(h1+h2+h3+h4h6)/r2h8=r6*(h1+h2+h3+h4+r5h6)/r2h9=r6*(h1+h2+h3+h4+r5h6)/r1h10=r6*(h1+h2+h3+h4+r5h6)/r4h11=r6*(h1+h2+h3+h4+r5h6)/r3h12=r6*(h1h6)/r2h13=r6*h6/r2h14=r6*h6/r1theta=45 ！ 繞軸旋轉(zhuǎn)45度!!電壓控制參數(shù)voltvalue=1!!網(wǎng)格劃分控制參數(shù)tcbj=20ouco=1sdpd=3dbj=20bhw=4wqs=18wqx=12hudu=20dihw=6sdps=12sdpx=2tchw=4ys=20js=12!!單元類型/PREP7 ET,1,SOLID45ET,2,SOLID5 ET,3,FLUID30 KEYOPT,3,2,0 ！單元類型號，任意值，近場水0遠場水1ET,4,FLUID30KEYOPT,4,2,1ET,5,FLUID130 !!實參量 R,1,20e6, ！近場水實常數(shù)，任意號，參考壓力R,2,20e6, ！遠場水實常數(shù)，任意號，參考壓力R,3,r6,0,0,0, ！邊界水實常數(shù)，半徑，中心坐標(biāo)X，Y，Z值! SAVE, definerealconstants,db,D:\ANSYSFILES\1030WATERNEW\ !!金屬材料 MPTEMP, MPTEMP,1,0 MPDATA,DENS,1,mat1dens MPTEMP, MPTEMP,1,0 MPDATA,EX,1,mat1yangMPDATA,PRXY,1,mat1person!!壓電片MPTEMP, MPTEMP,1,0 MPDATA,DENS,2,mat2dens TB,ANEL,2,1,21,0TBTEMP,0TBDATA,c11,c12,c13, TBDATA,c22,c23,c33 TBDATA,c44, TBDATA,c55,c66, MPTEMP, MPTEMP,1,0 MPDATA,PERX,2,fx MPDATA,PERY,2,fy MPDATA,PERZ,2,fzTB,PIEZ,2,0 TBMODIF,1,1,TBMODIF,1,2,TBMODIF,1,3,e13TBMODIF,2,1,TBMODIF,2,2,TBMODIF,2,3,e23TBMODIF,3,1,TBMODIF,3,2,TBMODIF,3,3,e33TBMODIF,4,1,TBMODIF,4,2,TBMODIF,4,3,TBMODIF,5,1,TBMODIF,5,2,e52TBMODIF,5,3,TBMODIF,6,1,e61TBMODIF,6,2,TBMODIF,6,3,!!近場水MPTEMP, MPTEMP,1,0 MPDATA,DENS,3,mat3dens MPTEMP, MPTEMP,1,0 MPDATA,SONC,3,sonc3 MPTEMP, MPTEMP,1,0 MPDATA,MU,3,mu3!!遠場水MPTEMP, MPTEMP,1,0 MPDATA,DENS,4,mat4dens MPTEMP, MPTEMP,1,0 MPDATA,SONC,4,sonc4MPTEMP, MPTEMP,1,0 MPDATA,MU,4,mu4!!邊界水MPTEMP, MPTEMP,1,0 MPDATA,DENS,5,mat5dens MPTEMP, MPTEMP,1,0 MPDATA,SONC,5,sonc5 MPTEMP, MPTEMP,1,0 MPDATA,MU,5,mu5 ！！K,1,0,0,h6, K,2,0,0,h1h6, K,3,0,0,h1+h2h6, K,4,0,0,h1+h2+h3h6, K,5,0,0,h1+h2+h3+h4h6, K,6,0,0,h1+h2+h3+h4+h5h6, K,7,r3,0,h1+h2h6, K,8,r3,0,h1+h2+h3h6,K,9,r3,0,h1+h2+h3+h4h6, K,10,r3,0,h1+h2+h3+h4+h5h6, K,11,r4,0,h1+h2+h3h6,K,12,r4,0,h1+h2+h3+h4h6, K,13,r4,0,h1+h2+h3+h4+h5h6, K,14,r1,0,h6,K,15,r1,0,h1h6K,16,r1,0,h1+h2+h3+h4h6,K,17,r1,0,h1+h2+h3+h4+h5h6, K,18,r2,0,h6,K,19,r2,0,h1h6, K,20,r2,0,h1+h2+h3+h4h6, K,21,r2,0,h1+h2+h3+h4+h5h6, K,22,0,0,r5,K,23,0,0,r6,K,24,r5,0,0,K,25,r6,0,0,K,26,r6,0,h7,K,27,r6,0,h8,K,28,r6,0,h9,K,29,r6,0,h10,K,30,r6,0,h11,K,31,0,0,r5,K,32,0,0,r6,K,33,0,0,0,K,34,r1,0,0,K,35,r2,0,0,K,36,r6,0,h12,K,37,r6,0,h13,K,38,r6,0,h14,!!LSTR, 1, 2LSTR, 2, 33LSTR, 33, 3LSTR, 3, 4 LSTR, 4, 5 LSTR, 5, 6 LSTR, 7, 8 LSTR, 8, 9 LSTR, 9, 10LSTR, 11, 12LSTR, 12, 13LSTR, 14, 15LSTR, 15, 34LSTR, 34, 16LSTR, 16, 17LSTR, 18, 19LSTR, 19, 35LSTR, 35, 20LSTR, 20, 21LSTR, 1, 14LSTR, 14, 18LSTR, 2, 15LSTR, 15, 19LSTR, 3, 7LSTR, 4, 8LSTR, 5, 9LSTR, 6, 10LSTR, 8, 11LSTR, 9, 12LSTR, 10, 13LSTR, 12, 16LSTR, 13, 17LSTR, 16, 20LSTR, 17, 21LSTR, 6, 22LSTR, 22, 23 LSTR, 33, 34LSTR, 34, 35LSTR, 35, 24LSTR, 24, 25LARC,22,24,34,r5 ！圓弧線，起始關(guān)鍵點號，終點關(guān)鍵點號，圓心號，半徑。VOLU39。（用藍色字表示注釋，紅色標(biāo)注提請注意）2 編輯命令流文件注意：所有以“！”開始的語句為注釋語句，并不參加執(zhí)行，可以去掉。因此我們在‘圖形操作’建模之后生成命令流文件，然后對命令流文件進行編輯整理，分塊，‘具體參數(shù)變量化’(如軸對稱模型只分析圓柱的一定角度，電極給定電壓，劃分線網(wǎng)格數(shù))等等操作之后，再由軟件讀入該文件就可以得到我們所要求解的東西。 也可選擇Main MenuGeneral PostprocPlot ResultsDeformed Shape 命令，彈出對話框選擇[Def+undef edge] ,單擊[OK]得到變形圖六 生成命令流文件以上步驟都是在圖形操作下完成的，但由于采用‘圖形操作’來建模效率太低，而且不易于改變參數(shù)（結(jié)構(gòu)參數(shù)，材料參數(shù)和網(wǎng)格參數(shù)等），因此我們有必要采用文本編程的方法來建模和求解。2. 動態(tài)演示模態(tài)結(jié)果(1） 讀取模態(tài)結(jié)果讀取第1階模態(tài)結(jié)果：主菜單－》General Postproc－》Read Results－》First Set 實際上，由于模型未在Z方向加位移限制條件，故第階1模態(tài)頻率為零，只作平動。程序進行求解，求解完成之后出現(xiàn)求解完畢提示，表示求解完成。主菜單欄－》Preprocessor－》Loads－》Define Loads－》Apply－》Electric－》Boundary－》Voltage－》On Areas選擇陶瓷下底面，如下出現(xiàn)對話框，要求輸入電壓值，電壓值填0V用同樣的方法選擇陶瓷上底面，加1V的電壓。出現(xiàn)如下對話框，選擇UX，表示固定其X方面，即在徑向上不移動。實用菜單－》WorkPlane－》Change Active CS to －》Global Cylindrical保存文件：12．節(jié)點轉(zhuǎn)換 把劃分網(wǎng)格后生成的節(jié)點轉(zhuǎn)換到柱坐標(biāo)系下，使得后續(xù)的求解過程計算速度更快。（為了采用6面體劃分體積單元，注意幾個劃分10格的線段劃分時要求單元數(shù)目一致，當(dāng)然這些線段也可以劃分成10格，但要一致）主菜單－》Preprocessor－》Meshing－》Mesh－》Volumes－》Mapped－》4 to 6 sided我們采用6面體單元劃分體，這樣我們劃分的網(wǎng)格大小就比較‘整齊’出現(xiàn)對話框，點擊Pick All，表示所有體積都按6面體單元劃分。（3）保存文件：8．指定各部分單元類型以及材料屬性（1）指定金屬殼的單元及材料屬性：主菜單欄－》Preprocessor－》Meshing－》Mesh Attributes－》Picked Volumes 用鼠標(biāo)點中金屬殼的體積單元出現(xiàn)對話框，材料參數(shù)選擇1，單元類型選擇1 SOLID45（2）指定壓電材料的單元及材料屬性：主菜單欄－》Preprocessor－》Meshing－》Mesh Attributes－》Picked Volumes 用鼠標(biāo)點中壓電材料的體積單元出現(xiàn)對話框，材料參數(shù)選擇2，單元類型選擇2 SOLID5（3） 點擊ok ，9 .劃分線段為了網(wǎng)格單元數(shù)可以控制，我們采用手動劃分網(wǎng)格，而不采用自動網(wǎng)格劃分。 生成如下圖的體。之