你好,我试了一下还是不行,用GUI和命令均是出现如下错误,下面是全部命令流,仿真目的是想实现在加载的情况下中间板向下随意运动,向上运动受到上挡板的阻挡,但不能穿透挡板,只是发生形变,麻烦你给看看,谢谢你。
FINISH
/clear
/FILNAME,concentrated load
/TITLE,concentrated load
/prep7
/units,MKS
!********************************************************************
!输入压电片结构尺寸及其它参数
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 !网格划分大小
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建立一个长方体,以对顶角的坐标为参数。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,如果输入的介电系数值小于1,程序会将该值认为是绝对介电系数。即实验室原来程序中输入也正确
mp,pery,1,634.7
mp,perz,1,728.5
tb,ANEL,1 ! Anisotropic elastic stiffness, N/m^2,弹性常数矩阵
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、压电常数矩阵
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
!*********************************************************************
!网格划分
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
!基板网格划分
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
|