你好,这是我写的建立接触对的程序,加入整体运行后,下面板在向上运动过程中仍能穿透上面的挡板,希望你帮着看看如何改进程序避免此现象。
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,solid45
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
mp,mu,1,0.2
mat,1
r,3
real,3
et,5,170
et,6,174
r,3,,,0.1,0.1,,
NROPT,UNSYM
ASEL,S,,,4
CM,_TARGET,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
ESURF,ALL
ASEL,S,,,57
CM,_CONTACT,AREA
TYPE,6
NSLA,S,1
ESLN,S,0
ESURF,ALL
CMDEL,_TARGET
CMDEL,_ CONTACT
ALLSEL,ALL
![]()
QQ截图20130227163509.png
|