请教压电悬臂梁问题

唐重鑫

本人在用ANSYS做压电悬臂梁能量采集器仿真，但发现网格划分有些问题，就是用了VGLUE命令后，网格划分有困难，请大家帮我看看以下的命令流（我想用来做模态分析的，大家帮我看看整个程序那些地方有问题，谢谢啦）
/prep7
! PZT-5H 材料特性参数
mp,DENS,1,7500                 ! 密度, kg/m**3
mp,perx,1,1100                 ! 介电常数
mp,pery,1,1100
mp,perz,1,827
tb,ANEL,1                      ! 弹性劲度系数, N/m^2
tbdata,1,13.2E10,7.1E10,7.31E10  ! c11,c12,c13
tbdata,7,13.2E10,7.3E10       ! c11,c33
tbdata,12,11.5E10              ! c33
tbdata,16,3.0E10              ! c44
tbdata,19,2.6E10              ! c44
tbdata,21,2.6E10              ! c66
tb,PIEZ,1                      ! 压电（应力）常数, C/m^2
tbdata,3,-5.2                  ! e13
tbdata,6,-5.2                 ! e23
tbdata,9,15.1                 ! e33
tbdata,14,12.7                ! e52
tbdata,16,12.7                ! e61
!定义主结构的材料参数
mp,dens,2,2300       !定义SI层材料参数  
MP,EX,2,130e9
MP,nuxy,2,0.27
mp,dens,3,8900       !定义NI层材料参数  
MP,EX,3,180e9
MP,nuxy,3,0.3

! 定义几何模型尺寸
  !压电层几何模型 (L1=3000um W1=600um H1 = 2um)
L1=3000e-6
W1=600e-6
H1=2e-6
  !SI层几何模型 (L2=3000um W2=600um H2 = 10um)
   L2=3000e-6
   W2=600e-6
   H2=10e-6
  !NI质量块尺寸（L3=600um W3=600um H3=500um)
L3=600e-6
W3=600e-6
H3=500e-6   
! 建立几何模型
block,0,L2,0,W2,0,H2                  !SI层
block,0,L2,0,W1,H2,H2+H1              !PZT层
block,L2-L3,L2,0,W3,H2+H1,H3+H2+H1    !NI层
vglue,all                             !将梁同压电片粘结
!定义有限元结构
et,1,solid5,3              !定义压电单元
et,2,solid95                  !定义主结构单元

!进行网格划分
LESIZE,ALL,10e-5, , , ,1, , ,1,
mat,1 $ type,1   !压电片网格划分
Vmesh,5
mat,2 $ type,2   !SI层网格划分
Vmesh,6
mat,3 $ type,2   !质量块网格划分
vmesh,4
nsel,s,loc,z,H2           !定义下层电极
cp,1,volt,all
*get,n_bot,node,0,num,min
nsel,s,loc,z,H2+H1         !定义上层电极
cp,2,volt,all
*get,n_top,node,0,num,min
fini
/solu
antype,modal
modopt,lanb,3
mxpand,3
nsel,s,loc,x,0
d,all,ux,0,,,,uy,uz
   !下层电极短路
d,n_bot,volt,0
nsel,all
solve
fini

ansyslong

不知道你解决了没有，帮你顶一下，同关注。

阿政狼

lhr799656

你这个我则么没看懂啊