1
案例背景
此例描述的是在旋轉座標系下帶有科里奧利效應的壓電耦合分析。幾何模型是由兩個齒和底座組成的石英音叉組成的,兩個齒由底座相連線。模型示意圖如圖1 ,其中,音叉模型寬度W=1250 um,齒寬W_t = 450um ;總長L=4800um; 齒長L_t= 3200um :厚度T = 350um。
石英材料的主要屬性有彈性係數、壓電係數、介電常數和密度。由於資料較大, 在這裡不做列舉,將在前處理中一一介紹。
圖1 模型示意圖
2
操作演示
window。DATA。videoArr。push({“title”:“ANSYS科里奧利效應的壓電耦合分析”,“vid”:“j3313orvjy1”,“img”:“http://inews。gtimg。com/newsapp_ls/0/14415436039_960495/0”,“desc”:“”})
3
模型結果圖片
4
APDL程式碼
!定義標題
/title, Coriolis Effect in a Vibrating Quartz Tuning Fork
/nopr
pi = 4*atan(1)
/VIEW,1,1,1,1
/TRIAD,lbot
!預處理
/PREP7
! 材料引數
c11 = 86。74e3
c12 = 6。99e3
c13 = 11。91e3
c14 = 17。91e3
c33 = 107。2e3
c44 = 57。94e3
tb,ANEL,1
tbdata, 1, c11, c12, c13, 0, c14, 0
tbdata, 7, c11, c13, 0,-c14, 0
tbdata,12, c33, 0, 0, 0
tbdata,16, (c11-c12)/2, 0, c14
tbdata,19, c44, 0
tbdata,21, c44
e11 = 0。171
e14 =-0。0406
tb,PIEZ,1
tbdata, 1, e11, 0, 0
tbdata, 4, -e11, 0, 0
tbdata, 7, 0, 0, 0
tbdata,10, 0, -e11, 0
tbdata,13, e14, 0, 0
tbdata,16, 0, -e14, 0
! -絕緣體常量
emunit,EPZRO,8。854e-6
mp,PERx,1, 4。43
mp,PERy,1, 4。43
mp,PERz,1, 4。63
! 密度
mp,DENS,1,2649e-18
! 尺寸引數
thick = 350
leng_TF = 4800
leng_tin = 3200
dist_t = 350
width_t = 450
x_t_in = dist_t/2
x_t_out = dist_t/2 + width_t
! 定義單元屬性
et,1,SOLID226,1001
!建模關鍵點
k, 1, 0, 0 , -thick/2
k, 2, 0, leng_TF-leng_tin , -thick/2
k, 3, x_t_in, 0 , -thick/2
k, 4, x_t_in, leng_TF-leng_tin , -thick/2
k, 5, x_t_in, leng_TF , -thick/2
k, 6, x_t_out, 0 , -thick/2
k, 7, x_t_out, leng_TF-leng_tin , -thick/2
k, 8, x_t_out, leng_TF , -thick/2
! 建立面
a,1,3,4,2
a,3,6,7,4
a,4,7,8,5
! 建立線
lesize, 5,,, 4,
*repeat,3,2
lesize, 1,,, 2,
lesize, 3,,, 2,
lesize, 8,,, 14, 3
lesize,10,,, 14, 1/3
lesize, 2,,, 8, -2
*repeat,3,2
*get,n_lin,LINE,,count
lgen,2,1, n_lin, 1,,, thick,20
l,1,21, 4,
*repeat,8,1,1
lsymm,X,all,,,100
! -建立體
v, 1, 3, 4, 2, 21, 23, 24, 22
v, 3, 6, 7, 4, 23, 26, 27, 24
v, 4, 7, 8, 5, 24, 27, 28, 25
v,101,103,104,102, 121,123,124,122
v,103,106,107,104, 123,126,127,124
v,104,107,108,105, 124,127,128,125
vplot
nummrg,kp
! 劃分網格
type,1
vmesh,all
! 定義電極
delta=20
! 電極載荷
nsel,s,loc,x, x_t_in+delta, x_t_out-delta
nsel,u,loc,z, -thick/2+1, thick/2-1
nsel,a,loc,x, -x_t_out-1, -x_t_out+1
nsel,a,loc,x, -x_t_in-1, -x_t_in+1
nsel,r,loc,y, leng_TF-leng_tin-1, leng_TF-leng_tin*0。45
cp,1,volt,all
n_load=ndnext(0)
! 電極接地
nsel,s,loc,x, -x_t_out+delta, -x_t_in-delta
nsel,u,loc,z, -thick/2+1, thick/2-1
nsel,a,loc,x, x_t_out-1, x_t_out+1
nsel,a,loc,x, x_t_in-1, x_t_in+1
nsel,r,loc,y, leng_TF-leng_tin-1, leng_TF-leng_tin*0。45 ! select tine-nodes
cp,2,volt,all
n_ground=ndnext(0)
nsel,all
! 求解設定
/SOLU
! —— Structural constraints
nsel,s,loc,y
d,all,ux,0,,,,uy,uz
nsel,all
! —— Ground electrode
d,n_ground,volt,0
! —— Loaded electrode
d,n_load,volt,1
! —— Add Coriolis effect
coriolis,on,,,off
omega,,1。e4
fini
!計算
/SOLU
antype,modal
modopt,QRDAMP,4
solve
fini
! 諧分析
/SOLU
antype,harm
dmprat,0。02
outres,all,all
solve
fini
!後處理
/POST1
set,1,1
/dscale,1,6
plns,uz
anharm
fini
案例來源:ANSYS 15。0多物理耦合場有限元分析從入門到精通