24小时热门版块排行榜

返回列表

老币2016

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 53.3
帖子: 5
在线: 11.9小时
虫号: 4362745
注册: 2016-01-17
性别: GG
专业: 机械工程

[求助] 本人想做一个压电陶瓷片的瞬态分析，apdl命令流如下，请教各位大神错误在哪里？已有1人参与

finish
/clear
/PREP7
et,1,solid5

!!!!!!!!!!!!!!!!!!!!!!material-1: Keramik PIC 181(+z 极化)
mp,dens,1,7850
eps0=8.854e-12
eps11=740
eps33=624
mp,perx,1,eps11*eps0
mp,pery,1,eps11*eps0
mp,perz,1,eps33*eps0

tb,anel,1
cep11=15.230e10
cep12=8.909e10
cep13=8.547e10
cep33=13.140e10
cep44=2.830e10
cep55=2.830e10
cep66=3.161e10

tbdata,1,cep11,cep12,cep13
tbdata,7,cep11,cep13
tbdata,12,cep33
tbdata,16,cep66
tbdata,19,cep44
tbdata,21,cep55

ep31=-4.5
ep33=14.7
ep15=11.0

tb,piez,1
tbdata,3,ep31
tbdata,6,ep31
tbdata,9,ep33
tbdata,14,ep15
tbdata,16,ep15

mm=0.001
r1=25*mm
h=1.2*mm

cylind,0,r1,0,h,0,90
vatt,1,,1

cylind,0,r1,0,h,90,180
vatt,1,,1

cylind,0,r1,0,h,180,270
vatt,1,,1

cylind,0,r1,0,h,270,360
vatt,1,,1

/Replot

lsel,s,loc,z,h/2
lesize,all,,,2,,1!lesize
lplot

nummrg,kp,1.0e-6
vsel,all
mshkey,1
esize,0.002
vmesh,all
!!!nsel,s,loc,z,0 !!选择陶瓷片底面
!d,all,ux,0       !x方向约束
!d,all,uy,0       !y方向约束

finish
vsel,all
/Replot

/solu                      !瞬态分析
antype,trans
trnopt,full
tintp,,0.25,0.5,0.5
outpr,all,all
outres,all,all
!kbc,1  !载荷形式阶跃载荷

nsel,s,loc,z,0             !耦合电压
cp,1,volt,all

nsel,s,loc,z,h
cp,2,volt,all

nsel,s,loc,y,0
d,all,volt,0
allsel,all

*do,i,1,50,1
t=i/1700000
time,t
nsubst,1
v=48*sin(2*3.14159*1700000*t)

nsel,s,loc,z,h
d,all,volt,v

allsel,all
solve
*enddo

finish

回复此楼