各位老师,以下是我的命令流,运行报错,烦请各位老师帮帮我,不甚感激!
/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 !定义Cu层材料参数
MP,EX,3,110e9
MP,nuxy,3,0.324
mp,dens,4,8900 !定义NI层材料参数
MP,EX,4,180e9
MP,nuxy,4,0.3
! 定义关键点
K,1,0,0,0 !SIO2层
k,2,4e-3,3e-4,0
K,3,4e-3,9e-4,0
k,4,0,12e-4,0
K,5,0,0,3e-6
k,6,4e-3,3e-4,3e-6
K,7,4e-3,9e-4,3e-6
k,8,0,12e-4,3e-6
K,9,0,0,13e-6 !Cu层
k,10,4e-3,3e-4,13e-6
K,11,4e-3,9e-4,13e-6
k,12,0,12e-4,13e-6
K,13,0,0,28e-6 !PZT层
k,14,4e-3,3e-4,28e-6
K,15,4e-3,9e-4,28e-6
k,16,0,12e-4,28e-6
K,17,3.5e-3,3e-4,28e-6 !NI层
k,18,4e-3,3e-4,28e-6
K,19,4e-3,9e-4,28e-6
k,20,3.5e-3,9e-4,28e-6
K,21,3.5e-3,3e-4,628e-6
k,22,4e-3,3e-4,628e-6
K,23,4e-3,9e-4,628e-6
k,24,3.5e-3,9e-4,628e-6
! 定义体
V,1,2,3,4,5,6,7,8
V,5,6,7,8,9,10,11,12
V,9,10,11,12,13,14,15,16
V,17,18,19,20,21,22,23,24
vglue,all !将梁同压电片粘结
!定义有限元结构
et,1,solid5,3 !定义压电单元
et,2,solid45 !定义主结构单元
!采用映射划分网格连接相邻面
asel,s,loc,z,28e-6
cm,CM_2,area
accat,CM_2
!进行网格划分
LESIZE,ALL,10e-5, , , ,1, , ,1,
mat,1 $ type,1 !压电片网格划分
Vmesh,6
mat,2 $ type,2 !SI层网格划分
Vmesh,1
mat,3 $ type,2 !Cu网格划分
vmesh,2
mat,4 $ type,2 !Ni网格划分
vmesh,5
nsel,s,loc,z,3e-6 !定义下层电极
cp,1,volt,all
*get,n_bot,node,0,num,min
d,n_bot,volt,0
nsel,s,loc,z,28e-6 !定义上层电极
cp,2,volt,all
*get,n_top,node,0,num,min
nsel,s,loc,x,0
d,all,ux,0,,,,uy,uz
!下层电极短路
nsel,all
fini
/solu
antype,modal
modopt,lanb,3
mxpand,3
solve
fini |