24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2241  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fzyling

新虫 (初入文坛)

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

用ansys做销轴在耳板内接触问题,总是出现类似图片里出现的问题,不能计算,不知道是不是接触的设置问题,请求帮助,下面附上命令流。
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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haoyanglijun

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by shgao20 at 2015-03-25 08:14:12
抱歉,漏了一点:另外还需要将 nsubst 加大,比如:
      nsubst,  1e4,1e8,10
求解时仍有可能出现你遇到的错误提示,点击 preceed,可以继续计算。

您好!请问这样计算的结果可靠吗?能满足精度要求吗?
12楼2016-02-15 13:42:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
从提示框里面里面 就可以 看出哪里有问题了 ,24335号节点在Z轴方向缺少约束。。。仔细检查一下。。。
数值计算
2楼2015-03-24 15:33:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzyling

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 独孤神宇 at 2015-03-24 15:33:51
从提示框里面里面 就可以 看出哪里有问题了 ,24335号节点在Z轴方向缺少约束。。。仔细检查一下。。。

轴在耳板里是接触关系,还用对其进行约束吗?
3楼2015-03-24 16:09:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
这是因为没有发生接触的原因,因而轴发生了刚体运动。
  首先检查一下接触对中,接触单元和目标单元的法线方向是否互相指向对方,如不是,则需要颠倒错误的法线方向,使它们互相指向对方。
  其次可以调整一下接触属性的设置,一个常用的设置是将接触单元的 keyopt(5) 设置为 3 (Close gap/reduce penetration with auto CNOF)。
  另外,还可以对求解参数做一些设置,比如:打开自动载荷步选项;设置较大的子步数,如: NSUBST, 1000,1e8,10
  等。
4楼2015-03-24 19:21:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见