24小时热门版块排行榜

返回列表

samooning

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 43
帖子: 10
在线: 4.3小时
虫号: 2815107
注册: 2013-11-20

[求助] 各位大神帮忙看下ansys模拟螺栓钢管接触问题分析为什么没法算！已有1人参与

奉上命令流先
*******************************************
fini
/clear
/filename, full creating contact
*SET,fchd,0.005 !辅材厚度
*set,Uhou,0.01 !U型板厚
*set,llll,0.5    !管长
*set,ls,0.02    !螺栓直径
*set,lsk,0.022 !实际螺栓孔径
*SET,cbdj,lsk*2 !插板端距
*SET,cbzxj,lsk*3  !插板中心距
*SET,br,0.102/2 !支3半径
*SET,xj,0.1    !插板与支管相交
*SET,Uh,0.06    !U型板高
ld=0.081       !板短边宽
lc=0.11          !板长边宽
la=0.1081       !板长
h=0.01          !板高
jj=0.015       !板管间距

/prep7
CYLIND,br-fchd,br,0,llll,-90,90,
!建立腹杆插板
wpro,,,90
wpro,,90
wpoffs,,,0.015
csys,4
BLOCK,-xj,0.015+2*cbdj,,br+0.02,-uhou,,
BLOCK,-xj,0.015+2*cbdj,br+0.01,br+0.02,-uhou,uh-uhou
VOVLAP,ALL

WPRO,,,-90
VSBW,ALL
WPOF,,,xj
VSBW,ALL
WPRO,,,90
VSBW,ALL
wpof,,,-uhou
vsbw,all
ALLS

WPCSYS,-1,0
wpof,,,-jj
VSBW,ALL

wpro,,-90
wpof,,,-h/2

nummrg,all
numcmp,all

k,100,
k,101,ld
k,102,lc,la
k,103,,la
k,104,,,h
a,100,101,102,103
l,100,104
vdrag,76,,,,,,119

wpof,cbzxj/2,0.015+cbdj-jj,-0.02
CYLIND, ,ls/2,-0.02,0.1,0,360,
vsel,s,,,15,17,2
cm,aa,volu
alls
vsbv,aa,18,,,1

kwpave,49
wpro,,,-90
vsbw,all
kwpave,35,
vsbw,all
kwpave,22,
wpro,,-90
vsbw,all
save

!板材切割
wpcsys,-1,0
vsel,s,loc,y,-0.011,0
vsel,a,,,14
wpof,cbzxj/2,,-0.015-cbdj
vsbw,all
wpro,,,90
vsbw,all
wpro,,,90
wpof,,,br+0.01-cbzxj/2
vsbw,all
wpof,,,-2*(br+0.01-cbzxj/2)
vsbw,all
wpro,,,90
wpof,,,br+0.01-cbzxj/2
vsbw,all
wpcsys,-1,0
wpof,cbzxj/2,,-0.015-cbdj
vsel,r,loc,x,-0.03,0.03
vsel,r,loc,z,-0.03,0.03
k,1000
kwplan,-1,1000,139,131
vsbw,all
kwplan,-1,1000,159,151
vsbw,all
alls
save

!建立螺栓模型
wpcsys,-1,0
wpof,cbzxj/2,,-0.015-cbdj
wpro,,-90
wpof,,,-0.025
cyl4,,,ls/2,360,,,0.04
cyl4,,,ls,,,,ls/2
wpoff,0,0,0.03
cyl4,,,ls,,,,ls/2
numcmp,all

vsel,s,,,57
vsel,a,,,59,70,11
vovlap,all
wpro,,-90
vsbw,all
wpro,,,90
vsbw,all
cm,bolt,volu
vsel,inve
!定义组件
wpcsys,-1,0
vsel,s,loc,y,0
cm,board,volu
vsel,inve
cmse,u,bolt
cm,cha,volu
alls
numcmp,all
!划分单元
et,1,solid185
mp,ex,1,2.06e11  !悬臂梁和柱
mp,prxy,1,0.3
tb,bkin,1
tbdata,1,235e6
mp,ex,2,2.2e5  !螺栓组
mp,prxy,2,0.25
tb,bkin,2
tbdata,1,400e6 !假设的
mp,ex,3,2.06e11  !端板
mp,prxy,3,0.3
tb,bkin,3
tbdata,1,235e6
r,1
r,2
r,3
save
!螺栓划分单元
cmsel,s,bolt
vatt,2,2,1
mshkey,1
esize,0.005
vmesh,all
!连接板划分单元
cmse,s,bolt
vsel,inve
kwplan,-1,165,159,127
vsel,r,loc,x,-0.03,0.03
vsel,r,loc,y,,0.06
vatt,1,1,1
mshkey,1
esize,0.007
vmesh,all

cmsel,s,board
vatt,1,1,1
mshkey,1
esize,0.007
vmesh,all
!插板划分单元
cmsel,s,cha
vatt,3,3,1
lsel,s,length,,0.4
lesize,all,0.03
wpcsys,-1,0
lsel,s,loc,y,-h
lesize,all,0.005
esize,0.007
vsel,r,loc,z,-jj,-1
vmesh,all
vsel,s,loc,z,0,1
vmesh,all
vsel,s,loc,z,-jj,0
SMRT,6
vsweep,all
save
!生成接触和目标单元
et,2,conta174
et,3,targe170
dfkn=1.0
dfto=0.01
!插板与连接板
r,4,,,dfkn,dfto
cmsel,s,board
nslv,s,1
nsel,r,loc,y,-h/2
real,4
type,2
esurf
cmsel,s,cha
nslv,s,1
nsel,r,loc,y,-h/2
nsel,r,loc,z,-jj,-1
type,3
esurf
!螺母与端板之间(标准接触)
r,5,,,dfkn,dfto
asel,,,,375,379,4
asel,a,,,383,386,3
nsla,s,1
real,5
type,2
esurf
cmse,s,board
nslv,s,1
nsel,r,loc,y,h/2
kwplan,-1,189,193,159
nsel,r,loc,y,,2*(br+0.01-cbzxj/2)
nsel,r,loc,x,br+0.01-cbzxj/2,-(br+0.01-cbzxj/2)
type,3
esurf
!螺栓和孔壁
r,6,,,dfkn,dfto
wpcsys,-1,0
wpof,cbzxj/2,,-0.015-cbdj
cmse,s,bolt
nslv,s,1
wpro,,-90
cswpla,11,1,1,1
nsel,s,loc,x,ls/2
nsel,r,loc,z,-0.015,0.005
real,6
type,3
esurf
cmse,s,bolt
vsel,inve
nslv,s,1
nsel,s,loc,x,ls/2
nsel,r,loc,z,-0.015,0.005
type,2
esurf
alls
!螺母与端板之间(标准接触)
r,7,,,dfkn,dfto
asel,s,,,361,365,4
asel,a,,,369,372,3
nsla,s,1
real,7
type,2
esurf
alls
cmse,s,cha
nslv,s,1
aslv,s
csys,4
asel,r,loc,z,-0.015
asel,r,loc,x,-0.03,0.03
asel,r,loc,y,-0.03,0.03
nsla,s,1
type,3
esurf

!施加约束和荷载
wpcsys,-1,0
csys,4
nsel,s,loc,z,-jj-la
D,all, , , , , ,ALL, , , , ,
asel,s,loc,x,0
da,all,symm
nsel,s,loc,z,0.5
F,all,FZ,10
alls
!求解
/solu
cncheck
nlgeom,on
outres,all,all
time,1
nsubst,50
pred,off
solve

******************
结构总是没法算，好像是约束不够，发生结构整体位移，但是实在找不出来问题出在那里了，求大神指导指导，谢谢！！！！

2JS7YTT]52(G~6F0~B0{S@W.png

回复此楼