你好,我試了一下還是不行,用GUI和命令均是出現(xiàn)如下錯誤,下面是全部命令流,仿真目的是想實現(xiàn)在加載的情況下中間板向下隨意運動,向上運動受到上擋板的阻擋,但不能穿透擋板,只是發(fā)生形變,麻煩你給看看,謝謝你。
FINISH
/clear
/FILNAME,concentrated load
/TITLE,concentrated load
/prep7
/units,MKS
!********************************************************************
!輸入壓電片結構尺寸及其它參數(shù)
th=0.8e-3
wd=20e-3 !矩形板寬度
lb=80e-3 !鈹青銅基板長度
lp=40 e-3 !壓電陶瓷長度
a1=0.5 !厚度比=陶瓷厚度/整體厚度(0.05,0.1,0.15……0.95)
thb=(1-a1)*th !基板厚度
thp=a1*th !壓電陶瓷厚度
esz=0.5e-3 !網(wǎng)格劃分大小
xa=(lb-lp)/2
block,0,0.012,0,thb,0,wd
block,0.012,xa,0,thb,0,wd
block,xa,xa+lp,0,thb,0,wd
block,xa+ lp,xa+lp+0.014,0,thb,0,wd
block,0.074, lb,0,thb,0,wd
block,xa,xa+lp,thb,thb+thp,0,wd
block,xa,xa+lp,0,-thp,0,wd
block,0.074,lb,thb,thb+15*thp,0,wd
block,0.074, lb,0,-15*thp,0,wd
block,0,0.012, thb,thb+0.004,0,wd
vglue,1,2,3,4,5,6,7,8,9
!BLOCK,X1,X2,Y1,Y2,Z1,Z2建立一個長方體,以對頂角的坐標為參數(shù)。X1,X2 為X 向最小及最大坐標值,Y1,Y2 為Y 向最小及最大坐標值, Z1,Z2 為X 向最小及最大坐標值
!*********************************************************************
!定義單元類型、材料屬性
ET,1,SOLID5,3 !對SOLID5通過設置KEYOPT(1)=3可用大撓度及應力剛化功能,激活至少一個帶有壓電自由度(位移和VOLT)的單元
mp,DENS,1,7500 ! Density, kg/m**3
mp,perx,1,728.5 !對于SOLID5,PLANE13或SOLID98,如果輸入的介電系數(shù)值小于1,程序會將該值認為是絕對介電系數(shù)。即實驗室原來程序中輸入也正確
mp,pery,1,634.7
mp,perz,1,728.5
tb,ANEL,1 ! Anisotropic elastic stiffness, N/m^2,彈性常數(shù)矩陣
tbdata,1,13.9E10,7.43E10,7.78E10 ! c11,c13,c12
tbdata,7,11.5E10,7.43E10 ! c33,c13
tbdata,12,13.9E10 ! c11
tbdata,16,2.56E10 ! c44
tbdata,19,2.56E10 ! c44
tbdata,21,3.06E10 ! c66
tb,PIEZ,1 ! Piezoelectric stress coefficients, C/m^2,u、壓電常數(shù)矩陣
tbdata,2,-5.2 ! e31
tbdata,5,15.1 ! e33
tbdata,8,-5.2 ! e31
tbdata,10,12.7 ! e15
tbdata,15,12.7 ! e15
et,2,solid185
mp,ex,2,1.3e11
mp,nuxy,2,0.3
mp,dens,2,8300
et,3,solid45
mp,ex,3,1.07e11
mp,nuxy,3,0.3
mp,dens,3,8900
et,4,solid185
mp,ex,4,1.07e11
mp,nuxy,4,0.3
mp,dens,4,8900
!*********************************************************************
!網(wǎng)格劃分
esize,esz
type,1
mat,1
vmesh,15
vmesh,16
type,3
mat,3
vmesh,18
vmesh,17
esize, esz
type,4
mat,4
vmesh,10
!基板網(wǎng)格劃分
type,2
mat,2
vmesh,1
vmesh,11
vmesh,12
vmesh,13
vmesh,14
vsel,S,LOC,Y, thb,thb+0.004
vsel,R,LOC,X, 0,0.012
CM,DANGBAN,VOLU
ALLS
vsel,S,LOC,Y,0,thb
vsel,R,LOC,X, 0,0.012
CM, JIBAN,VOLU
ALLS
CMSEL,R,DANGBAN
ASLV
ASEL,R,LOC,Y, thb
CM,CONTACT,AREA
ALLS
CMSEL,R,JIBAN
ASLV
ASEL,R,LOC,Y, thb
CM, TARGET,AREA
ALLS
mp,mu,1,0.2
et,5,173
et,6,170
r,7
real,7
KEYOPT,6,9,0
CMSEL,S,contact
TYPE,5
NSLA,S,1
ESLN,S,0
ESURF,ALL
ALLS
CMSEL,S,target
TYPE,6
NSLA,S,1
ESLN,S,0
ESURF,ALL
ALLSEL
nsel,s,loc,x,0
D,all,ux,0,,,,uy,uz
ALLSEL,ALL
SAVE
FINISH
/SOLU
csys,0
ANTYPE, TRANS ! Transient analysis
TRNOPT,full ! Full solution method
Tintp,,0.25,0.5,0.5
autot,on
alphad,10
Outpr,all,all
OUTRES,ALL,ALL
nsel,all
nsel,s,loc,y, 0 !Define bottom electrode定義基板上頂部電極,如果還是用csys,5,命令為:nsel,s,loc,z,th,具體參見ANSYS幫助CSYS命令
nsel,r,loc,x,0.012,0.074
cm,cm_1,node
nsel,s,loc,y, thb
nsel,r,loc,x,0.012,0.074
cm,cm_2,node
cmsel,s, cm_1, node
cmsel,a, cm_2, node
cp,1,volt,all
d,all,volt,0
nsel,all
nsel,s,loc,y,th ! Define top electrode定義頂部電極
nsel,r,loc,x,0.02,0.06
CM,n_2,NODE
cp,2,volt,all
nsel,all
nsel,s,loc,y,-thp
nsel,r,loc,x,0.02,0.06
CM,n_3,NODE
cp,3,volt,all
allsel,all
*do,i,0.0000001,1,0.05
TIME, i
KBC,1
v=100*(-1)**(nint(mod(i,1) )+2)
D,n_2,volt,v
D,n_3,volt,v
SOLVE
*enddo
finish
![]()
QQ截圖20130302122827.png
|