24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

fzyling

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 69.3
帖子: 9
在线: 3小时
虫号: 3744095
注册: 2015-03-17
专业: 机械动力学

[求助] ansys接触不能计算已有4人参与

用ansys做销轴在耳板内接触问题，总是出现类似图片里出现的问题，不能计算，不知道是不是接触的设置问题，请求帮助，下面附上命令流。

fini
/clear
/filname,k40_60z2

/prep7
csys,4
d1_cg=198.7
d2_cg=180.8
h1_cg=365
k,,d1_cg,0,0
k1_cg=kp(d1_cg,0,0)
k,,d2_cg,0,h1_cg
k2_cg=kp(d2_cg,0,h1_cg)
k,,d1_cg,d1_cg,0
k1_temp=kp(d1_cg,d1_cg,0)
kwplan,-1,k1_cg,k2_cg,k1_temp
l1_cg=100
l2_cg=96
t1_cg=8.5
t2_cg=5.3
wprota,,,90
blc4,-l2_cg/2,-l1_cg/2,l2_cg,l1_cg,h1_cg
blc4,-(l2_cg-t2_cg*2)/2,-(l1_cg-t1_cg*2)/2,l2_cg-t2_cg*2,l1_cg-t1_cg*2,h1_cg
vsbv,1,2,,delete,delete
wpcsys,-1,0
h1_csb=350+l2_cg/2
d1_csb=476
l1_csb=100
l2_csb=96
t1_csb=8.5
t2_csb=5.3
wpoffs,,,h1_csb
wprota,,,90
blc4,-l2_csb/2,-l1_csb/2,l2_csb,l1_csb,d1_csb/2
blc4,-(l2_csb-t2_csb*2)/2,-(l1_csb-t1_csb*2)/2,l2_csb-t2_csb*2,l1_csb-t1_csb*2,d1_csb/2
vsbv,1,2,,delete,delete
wpoffs,l2_csb/2
wprota,,,90
vsbw,3,,1
vdele,2,,,1
kwplan,-1,k1_cg,k2_cg,k1_temp
kwpave,39
vsbw,4,,1
vdele,3,,,1
nummrg,all

wpcsys,-1,0
l1_eb=150
l2_eb=140
l3_eb=150
r1_eb=65
r2_eb=60/2
t1_eb=20
t2_eb=10
kwpave,9
wprota,,,90
wpoffs,,-(l1_eb-l1_csb)/2
blc4,-l2_eb,0,l2_eb,l1_eb,t1_eb
wpoffs,-l2_eb,l1_eb/2
cyl4,0,0,l1_eb/2,90,l1_eb/2,270,t1_eb
vadd,3,4
cyl4,0,0,r2_eb,,,,t1_eb
vsbv,5,3,,delete,delete
wpoffs,,,t1_eb
cyl4,0,0,r2_eb,,r1_eb,,t2_eb
kwplan,-1,k1_cg,k2_cg,k1_temp
kwpave,27
k,,-l3_eb,0,0
k,,-l3_eb,l1_eb,0
k1_eb=kp(0,0,0)
k2_eb=kp(-l3_eb,0,0)
k3_eb=kp(-l3_eb,l1_eb,0)
k4_eb=kp(0,l1_eb,0)
lstr,k1_eb,k2_eb
lstr,k2_eb,k3_eb
lstr,k3_eb,k4_eb
lsla,u
lsel,a,,,31
al,all
alls
vdrag,37,,,,,,42
nummrg,all

wpcsys,-1,0
kwpave,21
wprota,,,-90
e1_neb=31
e2_neb=64
l1_neb=70
wpoffs,,,-e1_neb
wpoffs,,-(l1_eb-l1_csb)/2
blc4,0,0,l2_eb,l1_eb,-t1_eb
blc4,-l1_neb,0,l1_neb,(l1_eb-l1_csb)/2,-t1_eb
wpoffs,l2_eb,l1_eb/2
cyl4,0,0,l1_eb/2,-90,l1_eb/2,90,-t1_eb
vadd,6,8
cyl4,0,0,r2_eb,,,,-t1_eb
vsbv,9,6,,delete,delete
cyl4,0,0,r2_eb,,r1_eb,,t2_eb
vsymm,y,7
nummrg,all
vsel,s,loc,z,-t1_eb,0
vgen,2,all,,,,,-(e2_neb+t1_eb)
alls
wpoffs,,,-(e2_neb+t1_eb)
cyl4,0,0,r2_eb,,r1_eb,,e2_neb

e1_cgtb=70
e2_cgtb=10
l1_cgtb=110
l2_cgtb=150
l3_cgtb=30
t1_cgtb=10
wpcsys,-1,0
wpoffs,,,h1_csb
wpoffs,e1_cgtb,-l1_csb/2,-(l2_csb/2-e2_cgtb)
k,,0,0,0
k,,l1_cgtb,0,0
k,,l1_cgtb,0,-l2_cgtb
k,,l1_cgtb-l3_cgtb,0,-l2_cgtb
k,,0,0,-l3_cgtb
k1_cgtb=kp(0,0,0)
k2_cgtb=kp(l1_cgtb,0,0)
k3_cgtb=kp(l1_cgtb,0,-l2_cgtb)
k4_cgtb=kp(l1_cgtb-l3_cgtb,0,-l2_cgtb)
k5_cgtb=kp(0,0,-l3_cgtb)
lstr,k1_cgtb,k2_cgtb
lstr,k2_cgtb,k3_cgtb
lstr,k3_cgtb,k4_cgtb
lstr,k4_cgtb,k5_cgtb
lstr,k5_cgtb,k1_cgtb
lsla,u
al,all
vext,98,,,,-t1_cgtb
vgen,2,14,,,,l1_csb+t1_cgtb
alls
vglue,all
nummrg,all

wpcsys,-1,0
wpoffs,0,0,h1_csb+l2_csb/2+l2_eb
wprota,,,90
l1_xz=265
r1_xz=60/2
cylind,r1_xz,,0,l1_xz
l1_g=24/2
l2_g=9
l3_g=24
wpoffs,,,l1_g
vsel,s,,,1
vsbw,all,,1
wpoffs,,,e1_neb-l1_g+t1_eb+e2_neb+t1_eb
vsbw,all,,1
wpoffs,,,24
vsbw,all,,1
wpoffs,,,l2_g
vsbw,all,,1
wpoffs,,,24
vsbw,all,,1
wpoffs,,,l2_g
vsbw,all,,1
kwpave,25
vsbw,all,,1
cm,xz_vo,volu
alls

mp,ex,1,2.06e5
mp,nuxy,1,0.3
mp,dens,1,7.85e-6
mp,mu,1,0.1
mp,ex,2,2.06e5
mp,nuxy,2,0.3
mp,dens,2,7.85e-6

mat,2
et,1,186
cmsel,u,xz_vo
vsel,u,,,25,28,3
mshape,1,3d
mshkey,0
esize,15
vmesh,all
alls
cmsel,u,xz_vo
mshape,1,3d
mshkey,0
esize,30
vmesh,all
cmsel,s,xz_vo
mshape,1,3d
mshkey,0
esize,10
vsweep,all
alls

wpcsys,-1,0
wpoffs,,,h1_csb+l2_csb/2+l2_eb
wprota,,,90
!!第一个接触对
mat,1
et,3,targe170
et,4,174
keyopt,4,5,3
r,3,,,0.1
real,3
!目标面
vsel,s,,,20,26,6
vsel,a,,,17,27,10
allsel,below,volume
wpstyle,0.05,0.1,-1,1,0.003,1,2,,5
asel,r,loc,x,r2_eb
nsla,s,1
esln,s
type,3
esurf
alls
!接触面
vsel,s,,,1
allsel,below,volume
asel,r,loc,x,r2_eb
nsla,s,1
esln,s
type,4
esurf
alls
!!第二个接触对
et,5,170
et,6,174
keyopt,6,5,3
r,4,,,0.1
real,4
!目标面
vsel,s,,,16,24,8
  allsel,below,volume
asel,r,loc,x,r2_eb
nsla,s,1
esln,s
type,5
esurf
alls
  !接触面
vsel,s,,,11
allsel,below,volume
asel,r,loc,x,r2_eb
nsla,s,1
esln,s
type,6
esurf
alls

wpcsys,-1,0
asel,s,,,13
da,all,all
alls
asel,s,,,1,23,22
da,all,symm
alls

load_x1=23362*9.8
load_y1=27655*9.8
load_x2=10390*9.8
load_y2=178*9.8
load_x3=12971*9.8
load_y3=3170*9.8
wpcsys,-1,0
asel,s,,,55
nsla,s,1
*get,ncount1,node,,count
f,all,fy,load_x3/ncount1
f,all,fz,-load_y3/ncount1
alls
asel,s,,,6
nsla,s,1
*get,ncount2,node,,count
f,all,fy,load_x2/ncount2
f,all,fz,-load_y2/ncount2
alls
asel,s,,,30,57,27
nsla,s,1
*get,ncount3,node,,count
f,all,fy,-load_x1/ncount3
f,all,fz,-load_y1/ncount3
alls

/solu
outres,basic,all
autots,on
  solc,on
  pred,on
  lnsrch,1
  nropt,unsym
  eqslv,pcg,1e-4
  time,1
  nsubst,20,200,10

solve

回复此楼