24小时热门版块排行榜    

查看: 729  |  回复: 1

奇璇璇

新虫 (初入文坛)

[求助] 做压电悬臂梁分析,程序编好后查看电场老是提示错误,老师要求得到电压值

这是程序
finish
/prep7       !进入前处理器
!pzt下
L1=0.04
W1=0.008
H1=0.0003
!CU层
L2=0.05
W2=0.008
H2=0.0005
!PZT上
L3=0.04
W3=0.008
H3=0.0003
!NI 1
L4=0.005
W4=0.008
H4=0.008
BLOCK,L2-L1,L2,0,W1,0,H1
!PZT下
BLOCK,0,L2,0,W2,H1,H1+H2
!cu层
BLOCK,L2-L3,L2,0,W3,H1+H2,H1+H2+H3
!pzt上
BLOCK,L2-L4,L2,0,W4,H1+H2+H3,H1+H2+H3+H4
!NI 1
vglue,all          !将所有的体粘接
!PZT下材料特性参数
mp,dens,1,7500   !介电常数矩阵mp,perx,1,1100
mp,pery,1,1100
mp,perz,1,827
tb,anel,1        !刚度矩阵??还是弹性劲度系数??tbdata,1,13.2e10,7.1e10,7.3e10
tbdata,7,13.2e10,7.3e10
tbdata,12,11.5e10
tbdata,16,3.0e10
tbdata,19,2.6e10
tbdata,21,2.6e10
tb,piez,1        !压电应力常数,压电常数矩阵的输入
tbdata,3,-5.2
tbdata,6,-5.2
tbdata,9,15.1
tbdata,14,12.7
tbdata,16,12.7!PZT上材料特性参数
mp,dens,3,7500     !介电常数矩阵mp,perx,3,1100
mp,pery,3,1100
mp,perz,3,827tb,anel,3       !刚度矩阵??还是弹性劲度系数??tbdata,1,13.2e10,7.1e10,7.3e10
tbdata,7,13.2e10,7.3e10
tbdata,12,11.5e10
tbdata,16,3.0e10
tbdata,19,2.6e10
tbdata,21,2.6e10
tb,piez,3      !压电应力常数,压电常数矩阵的输入
tbdata,3,-5.2
tbdata,6,-5.2
tbdata,9,15.1
tbdata,14,12.7
tbdata,16,12.7!定义结构的材料参数
mp,dens,1,7500   !前面已经定义了pzt的弹性矩阵(刚度矩阵)所以此处就不要再
mp,nuxy,1,0.287
mp,dens,2,8300
mp,ex,2,110e9
mp,nuxy,2,0.324
mp,dens,3,7500
mp,nuxy,3,0.287
mp,dens,4,8900
mp,ex,4,180e9
mp,nuxy,4,0.3!网格划分,定义有限元结构
et,1,solid5,3et,2,solid45!进行网格划分
lesize,all,5e-4,,,,1,,,1,mat,1 $type,1  !下压电片网格划分
mat,3 $type,1  !上压电片网格划分
mat,2$type,2  !cu层网格划分
mat,4$type,2  !ni质量块网格划分
vsweep,all!边界条件、载荷及求解
nsel,s,loc,z,0!选中z坐标为0的所有点,定义下压电片的下层电极
cp,1,volt,all
*get,n_bot,node,0,num,mind,n_bot,volt,0!下压电片下层电极短路
nsel,s,loc,z,h1!定义下压电片上层电极
cp,2,volt,all
*get,n_top,node,0,num,min
nsel,s,loc,z,h1+h2!定义上压电片下层电极
cp,3,volt,all
*get,n_bot,node,0,num,min
d,n_bot,volt,0!上压电片下层电极短路
nsel,s,loc,z,h1+h2+h3!定义上压电片上层电极
cp,4,volt,all
*get,n_top,node,0,num,min!施加固定约束
da,35,ux,0
da,35,uy,0
da,35,uz,0
da,36,ux,0
da,36,uy,0
da,36,uz,0
!考虑重力载荷
acel,0,0,9.8
nsel,all
!静态分析
/solu
antype,static!静态分析
solve
fini
/post1
pldisp,1!显示变形和未变形形状
plnsol,u,sum!绘制位移等值线图
finish
!模态分析
/solu
antype,modal
modopt,lanb,4
mxpand,4
solve
finish
!下压电片输出电压
!谐响应分析
/soluan
type,harmic
harfrq,0,3000!频率范围
nsubst,20!子步数量
kbc,1
solve
finish
分析完以后查看电场模块就会提示 The requested EF data is not available.  The PLNSOL command is ignored求高手指教!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

发错版了,去相关专业版吧
2楼2013-01-06 16:06:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 奇璇璇 的主题更新
信息提示
请填处理意见