请教
本人在做 耳片分析工程中，需要在销钉与耳片之间加一个干涉衬套，这样就形成了两个接触对。可是在分析中出现了错误，我实在无法解决，请各位高手帮忙解答，谢谢！。以下是APDL。r1_r2=2.5 ! 耳片半径 比 销钉空半径
e1_e2=0.36 ! 耳片材料 E 比 销钉材料 E
mu1_mu2=1.18 ! 耳片材料 v 比 销钉材料 v
r2=10 ! 销钉空半径
e2=200000 ! 销钉材料 弹性模量
mu2=0.28 ! 销钉材料 泊松比
mu1=mu1_mu2*mu2 ! 耳片材料 泊松比
e1=e1_e2*e2 ! 耳片材料 弹性模量
r1=r1_r2*r2 ! 耳片半径
!建模
/prep7
k,1,0,r1,0
k,2,0,-r1,0
k,3,-50,-r1,0
k,4,-50,r1,0
k,5,0,0,0
k,6,r1,0,0
l,4,1
larc,1,6,5,r1
larc,6,2,5,r1
l,2,3
l,3,4
al,1,2,3,4,5
cyl4,kx(5),ky(5),0,0,r2,360
asba,1,2
numcmp,all
hd_xd=0.2 ! 衬套的相对厚度
gsl=0.002 ! 干涉量
ct_ri=r2/(1-gsl+hd_xd) ! 衬套内径
ct_ro=ct_ri+hd_xd*ct_ri ! 衬套外径
cyl4,kx(5),ky(5),ct_ro,0,ct_ri,360 !创建干涉衬套
cyl4,kx(5),ky(5),0,0,ct_ri-(ct_ro-r2),360 !创建销钉
et,1,183
mp,ex,1,e1 ! 耳片
mp,prxy,1,mu1
mp,mu,1,0.2
mp,ex,2,e2 ! 销钉
mp,prxy,2,mu2
mp,ex,3,e2 !衬套
mp,prxy,3,mu2
! 网格划分
mshape,1,2d ! 指定三角形单元分网
mshkey,0 ! 指定自由网格划分方式
type,1
mat,1
esize,1.5
amesh,1
esize,0.5
type,1
mat,3
amesh,2
esize,0.5
type,1
mat,2
amesh,3
et,2,169 ! 2D目标面接触单元
et,3,172 ! 2D接触面接触单元
r,3,,,1,1 ! 定义接触对实常数
nropt,unsym !牛顿—拉普森 矩阵 为非对称矩阵
lsel,s,,,10
lsel,a,,,11
lsel,a,,,12
lsel,a,,,13
nsll,s,1
esln,s,0
type,3
mat,2
real,3
esurf,all,top ! 耳片与衬套接触对
allsel
lsel,s,,,6
lsel,a,,,7
lsel,a,,,8
lsel,a,,,9
nsll,s,1
type,2
mat,1
real,3
esurf,all,top
! 衬套与销钉接触对
lsel,s,,,18
lsel,a,,,19
lsel,a,,,20
lsel,a,,,21
nsll,s,1
esln,s,0
type,3
mat,2
real,3
esurf,all,top
allsel
lsel,s,,,14
lsel,a,,,15
lsel,a,,,16
lsel,a,,,17
nsll,s,1
type,2
mat,2
real,3
esurf,all,top
allsel
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1
allsel
lsel,s,,,19
lsel,a,,,20
sfl,all,pres,150
allsel
lsel,s,,,5
nsll,s,1
d,all,all
! 定义求解
/solu
allsel
antype,0 !指定分析类型为 一次静力分析
nlgeom,on !开启大变形效应
autots,on !关闭自动时间步长跟踪
time,3
nsubst,100,1000,3 ! 设置子步
solve 那位大侠会呀，帮忙看看，小弟很急呀、、、、、、、、:(
页:
[1]