螺栓和法兰连接的接触分析(带预紧)

jtxqhm

从<ANSYS 高级工程有限元分析范例精选>这本书中打了这个命令流,比较匆忙,没写注释,大家帮忙补充一下吧,也许有用
*set,rad_i,1250/2   
*set,rad_pani,705   
*set,rad_o,2450/2   
*set,rad_pano,1830  
*set,rad_b,180/2
*set,m_bolt,160
*set,m_bh1,350+170+50   
*set,m_bh2,350+170+250  
*set,m_nut,280/2
*set,h_nut,170  
*set,rad_drill,2140/2   
*set,rad_m,750+150  
*set,h_m,18
*set,rad_bi,1950/2
*set,rad_bo,2330/2  
*set,rad_bh,5   
*set,f_rad,45   
*set,fh_rad,1855/2  
*set,r_fillt,125
*set,n,15   
*set,key_w,290  
*set,key_h,130  
*set,th,200
*set,flang_h1,350   
*set,flang_h2,350   
*set,h_shaft,350+490
*set,h_sh,H_shaft+500   
*set,pi,acos(-1)
*set,elemsize,60
/graphics,power
*set,disp_b,0.75
*set,zmax,h_sh  
*set,f_ext,22.9215e+5*12
*set,f_eef,30.827e+5*12
*set,r_out,rad_pano
*set,t_w,375
*set,u_r,r_out-t_w  
/prep7  
et,1,185
mp,ex,1,2.06e5  
mp,nuxy,1,0.29  
mp,dens,1,7.85e-9   
mp,ex,2,2.06e5  
mp,nuxy,2,0.29  
mp,dens,2,7.85e-9   
mp,ex,3,2.06e5  
mp,nuxy,3,0.29  
mp,dens,3,7.85e-9   
mp,mu,4,0.15
k,1,rad_i,0,h_m
k,2,rad_m,0,0   
k,3,rad_i+th
k,4,rad_o,  
k,5,rad_o,0,flang_h1
k,6,rad_m,0,flang_h1
k,7,rad_i+th,0,flang_h1
k,8,rad_i+th,0,h_shaft  
k,9,rad_i,0,h_shaft
k,10,rad_i,0,flang_h1   
k,11,rad_m,0,H_M
k,12,rad_i+th,0,H_sh
k,13,rad_i,0,h_sh   
k,14,rad_bi,0,rad_bh
k,15,rad_bo,0,rad_bh
k,16,rad_bi
k,17,rad_bo
k,18,rad_bi,0,flang_h1  
k,19,rad_bo,0,flang_h1  
k,20,rad_i+th+f_rad,0,flang_h1  
k,21,fh_rad
k,22,fh_rad,0,flang_h1+h_shaft  
k,23,fh_rad,0,flang_h1  
l,8,20  
l,20,5  
l,21,22
lfillt,1,2,125  
lstr,1,10   
lstr,10,9   
lstr,9,13   
lstr,13,12  
lstr,12,8   
lstr,8,9
lstr,8,24   
lstr,1,11   
lstr,11,2   
lstr,2,21   
lstr,21,4   
lstr,4,5
lstr,10,23  
lstr,5,23   
lptn,all
al,7,8,9,10
al,1,22,21,16,6,10  
al,5,11,12,13,18,16
a,5,4,21,23
alls
lsla
lsel,invert
ldele,all   
ksll
ksel,invert
kdele,all   
APLOT   
alls
csys,1  
kgen,2,1,,,,n,,,0   
l,1,3   
vdrag,all,,,,,,2
csys,0  
cyl4,,,rad_bi,,rad_bo,,rad_bh   
v***v,4,5
cyl4,rad_drill,0,rad_b,,,,flang_h1+100  
csys,1  
vgen,,4,,,,n,,,,1   
v***v,6,4
wpro,n,,
csys,4  
vsymm,y,all,,,,0,0  
wpro,-n,,   
csys,0  
rectng,0,rad_o+100,0,key_w/2
vsel,all
aslv
asel,invert
vext,all,,,0,0,key_h
alls
vptn,all
vdele,11,,,1
vdele,12,,,1
vdele,15,,,1
vdele,16,,,1
vdele,18,,,1
wpoff,0,key_w/2,0   
wpro,,90
v***w,all
wpcsys,-1,0
wpoff,0,0,key_h
v***w,all
wpcsys,-1,0
csys,0  
cyl4,,,rad_bi,,rad_bo,,flang_h1
vptn,all
vdele,33,,,1
wpoff,0,0,rad_bh
v***w,all
wpoff,0,0,-rad_bh   
cylind,rad_pani,rad_pano,0,-flang_h2,0,2*n  
cyl4,rad_drill,0,rad_b,,,,flang_h2+100  
csys,1  
vgen,,16,,,,n,-flang_h2,,,1
v***v,15,16  
vsel,s,loc,z,0,h_sh+100
cm,shaft,volu   
alls
vsel,s,loc,z,0,-flang_h2-100
cm,pan,volu
alls
esize,elemsize,0,   
mat,1   
cmsel,,shaft
vsweep,all  
mat,2   
alls
cmsel,,pan  
vsweep,all  
alls
cmsel,,shaft
aslv
lsla
ksll
wpro,2*n,,  
csys,1  
cmsel,u,pan
vsel,r,loc,y,n,2*n  
csys,4  
vsymm,y,all,,,,0,0  
wpro,n,,
vsel,s,loc,z,0,h_sh+100
vsymm,y,all,,,,0,0  
eslv
nsle
nummrg,kp   
nummrg,node
alls
wpcsys,-1,0
cmsel,,pan  
aslv
lsla
ksll
wpro,2*n
csys,4  
vsymm,y,all,,,,0,0  
wpro,2*n
csys,4  
csys,1  
vsel,r,loc,y,2*n,4*n
csys,4  
vsymm,y,all,,,,0,0  
alls
wpcsys,-1,0
vsel,s,loc,z,0,-flang_h2-100
aslv
lsla
ksll
eslv
nsle
nummrg,kp   
nummrg,node
vsel,s,loc,z,0,h_sh+100
cm,shaft,volu   
alls
vsel,s,loc,z,0,-flang_h2-100
cm,pan,volu
alls
csys,0  
cylind,m_bolt/2,,m_bh1,-m_bh2,0,360
rpr4,6,,,m_nut,2*n,h_nut
vgen,,95,,,,,flang_h1,,,1   
vgen,2,95,,,,,-flang_h2-flang_h1-h_nut,,0   
cmsel,u,pan
cmsel,u,shaft   
vptn,all
cm,b_temp,volu  
vgen,,all,,,rad_drill,,,,,1
csys,1  
vgen,,all,,,,n,,,,1
mat,3   
esize,elemsize/2,0,
vsweep,all  
vgen,4,all,,,,n*2,,,0   
vclear,b_temp   
vdele,b_temp
vgen,,all,,,,-n*2,,,,1  
cm,bvolu,volu   
alls
psmesh,1,b100,,volu,95,0,z,0,,,,b_100   
psmesh,2,b200,,volu,116,0,z,0,,,,b_200  
psmesh,3,b300,,volu,109,0,z,0,,,,b_300  
cmsel,r,shaft   
aslv
asel,r,loc,z,0  
cm,as_con,area  
alls
cmsel,r,pan
aslv
asel,r,loc,z,0  
cm,ap_tar,area  
alls
/com,contact pair creation -start   
mat,4   
r,3
real,3  
et,3,170
et,4,174
keyopt,4,9,0
cmsel,s,ap_tar  
type,3  
nsla,s,1
esln,s,0
esurf,all   
alls
cmsel,s,as_con  
type,4  
nsla,s,1
esln,s,0
esurf,all   
allsel  
cmsel,r,bvolu   
aslv
asel,r,loc,z,flang_h1   
cm,ab1_con,area
alls
cmsel,r,shaft   
aslv
asel,r,loc,z,flang_h1   
csys,1  
asel,r,loc,x,fh_rad,rad_o   
csys,0  
cm,abs_tar,area
alls
/com, contact pair creation - start
mat,4   
r,4
real,4  
et,5,170
et,6,174
keyopt,5,9,0
cmsel,s,abs_tar
type,5  
nsla,s,1
esln,s,0
esurf,all   
alls
cmsel,s,ab1_con
type,6  
nsla,s,1
esln,s,0
esurf,all   
allsel  
/com, contact pair creation - end   
cmsel,r,bvolu   
aslv
asel,r,loc,z,-flang_h2  
cm,ab2_con,area
alls
cmsel,r,pan
aslv
asel,r,loc,z,-flang_h2  
cm,abp_tar,area
alls
/com, contact pair creation - start
mat,4   
r,5
real,5  
et,7,170
et,8,174
keyopt,7,9,0
cmsel,s,abp_tar
type,7  
nsla,s,1
esln,s,0
esurf,all   
alls
cmsel,s,ab2_con
type,8  
nsla,s,1
esln,s,0
esurf,all   
allsel  
/com, contact pair creation - end   

/SOL
alls
lsclear,all
antype,0
autots,1
nsubst,25,100,1
outres,all,all  
eqslv,pcg,1e-8  
time,15
csys,1  
asel,s,loc,y,0  
da,all,symm
alls
asel,s,loc,y,90
da,all,symm
alls
sload,1,9,lock,disp,disp_b,1,2  
sload,2,9,lock,disp,disp_b,1,2  
sload,3,9,lock,disp,disp_b,1,2  
cmsel,s,pan
nslv,r,1
nsel,r,loc,z,0  
csys,1  
nsel,r,loc,x,u_r,r_out  
d,all,all   
alls
***ctran
lswrite,1   
lswrite,2   
nsel,s,loc,z,zmax-1,zmax+10
*get,nload,node,,count  
f,all,fz,f_ext/nload/4  
alls
lswrite,3   
lssolve,1,3,1
/post
set,last
plnsol,s,eqv,0,1

yjwung9

楼主辛苦了。
建议给积分。

yjwung9

楼主辛苦了。
建议给积分。

jtxqhm

^_^

jinqinghai

it should be corrected seriously

jtxqhm

我在Ansys9.0上可以运行的
这个是网格图和mise应力云图

jtxqhm

mise云图

jinqinghai

!!! checked apdl , all "***" should be "S B".

*set,rad_i,1250/2
*set,rad_pani,705
*set,rad_o,2450/2
*set,rad_pano,1830
*set,rad_b,180/2
*set,m_bolt,160
*set,m_bh1,350+170+50
*set,m_bh2,350+170+250
*set,m_nut,280/2
*set,h_nut,170
*set,rad_drill,2140/2
*set,rad_m,750+150
*set,h_m,18
*set,rad_bi,1950/2
*set,rad_bo,2330/2
*set,rad_bh,5
*set,f_rad,45
*set,fh_rad,1855/2
*set,r_fillt,125
*set,n,15
*set,key_w,290
*set,key_h,130
*set,th,200
*set,flang_h1,350
*set,flang_h2,350
*set,h_shaft,350+490
*set,h_sh,H_shaft+500
*set,pi,acos(-1)
*set,elemsize,60
/graphics,power
*set,disp_b,0.75
*set,zmax,h_sh
*set,f_ext,22.9215e+5*12
*set,f_eef,30.827e+5*12
*set,r_out,rad_pano
*set,t_w,375
*set,u_r,r_out-t_w
/prep7
et,1,185
mp,ex,1,2.06e5
mp,nuxy,1,0.29
mp,dens,1,7.85e-9
mp,ex,2,2.06e5
mp,nuxy,2,0.29
mp,dens,2,7.85e-9
mp,ex,3,2.06e5
mp,nuxy,3,0.29
mp,dens,3,7.85e-9
mp,mu,4,0.15
k,1,rad_i,0,h_m
k,2,rad_m,0,0
k,3,rad_i+th
k,4,rad_o, 0,
k,5,rad_o,0,flang_h1
k,6,rad_m,0,flang_h1
k,7,rad_i+th,0,flang_h1
k,8,rad_i+th,0,h_shaft
k,9,rad_i,0,h_shaft
k,10,rad_i,0,flang_h1
k,11,rad_m,0,H_M
k,12,rad_i+th,0,H_sh
k,13,rad_i,0,h_sh
k,14,rad_bi,0,rad_bh
k,15,rad_bo,0,rad_bh
k,16,rad_bi
k,17,rad_bo
k,18,rad_bi,0,flang_h1
k,19,rad_bo,0,flang_h1
k,20,rad_i+th+f_rad,0,flang_h1
k,21,fh_rad ,
k,22,fh_rad,0,flang_h1+h_shaft
k,23,fh_rad,0,flang_h1
l,8,20
l,20,5
l,21,22
lfillt,1,2,125
lstr,1,10
lstr,10,9
lstr,9,13
lstr,13,12
lstr,12,8
lstr,8,9
lstr,8,24
lstr,1,11
lstr,11,2
lstr,2,21
lstr,21,4
lstr,4,5
lstr,10,23
lstr,5,23
lptn,all
al,7,8,9,10
al,1,22,21,16,6,10
al,5,11,12,13,18,16
a,5,4,21,23
alls
lsla
lsel,invert
ldele,all
ksll
ksel,invert
kdele,all
APLOT
alls
csys,1
kgen,2,1,,,,n,,,0
l,1,3
vdrag,all,,,,,,2
csys,0
cyl4,,,rad_bi,,rad_bo,,rad_bh
v***v,4,5
cyl4,rad_drill,0,rad_b,,,,flang_h1+100
csys,1
vgen,,4,,,,n,,,,1
v***v,6,4
wpro,n,,
csys,4
vsymm,y,all,,,,0,0
wpro,-n,,
csys,0
rectng,0,rad_o+100,0,key_w/2
vsel,all
aslv
asel,invert
vext,all,,,0,0,key_h
alls
vptn,all
vdele,11,,,1
vdele,12,,,1
vdele,15,,,1
vdele,16,,,1
vdele,18,,,1
wpoff,0,key_w/2,0
wpro,,90
v***w,all
wpcsys,-1,0
wpoff,0,0,key_h
v***w,all
wpcsys,-1,0
csys,0
cyl4,,,rad_bi,,rad_bo,,flang_h1
vptn,all
vdele,33,,,1
wpoff,0,0,rad_bh
v***w,all
wpoff,0,0,-rad_bh

!!! set up the lower circle plate

cylind,rad_pani,rad_pano,0,-flang_h2,0,2*n
cyl4,rad_drill,0,rad_b,,,,flang_h2+100
csys,1
vgen,,16,,,,n,-flang_h2,,,1
v***v,15,16

!!! set up components

vsel,s,loc,z,0,h_sh+100
cm,shaft,volu
alls
vsel,s,loc,z,0,-flang_h2-100
cm,pan,volu
alls
esize,elemsize,0,
mat,1
cmsel,,shaft
vsweep,all
mat,2
alls
cmsel,,pan
vsweep,all
alls
cmsel,,shaft
aslv
lsla
ksll
wpro,2*n,,
csys,1
cmsel,u,pan
vsel,r,loc,y,n,2*n
csys,4
vsymm,y,all,,,,0,0
wpro,n,,
vsel,s,loc,z,0,h_sh+100
vsymm,y,all,,,,0,0
eslv
nsle
nummrg,kp
nummrg,node
alls
wpcsys,-1,0
cmsel,,pan
aslv
lsla
ksll
wpro,2*n
csys,4
vsymm,y,all,,,,0,0
wpro,2*n
csys,4
csys,1
vsel,r,loc,y,2*n,4*n
csys,4
vsymm,y,all,,,,0,0
alls
wpcsys,-1,0
vsel,s,loc,z,0,-flang_h2-100
aslv
lsla
ksll
eslv
nsle
nummrg,kp
nummrg,node
vsel,s,loc,z,0,h_sh+100
cm,shaft,volu
alls
vsel,s,loc,z,0,-flang_h2-100
cm,pan,volu
alls
csys,0
cylind,m_bolt/2,,m_bh1,-m_bh2,0,360
rpr4,6,,,m_nut,2*n,h_nut
vgen,,95,,,,,flang_h1,,,1
vgen,2,95,,,,,-flang_h2-flang_h1-h_nut,,0
cmsel,u,pan
cmsel,u,shaft
vptn,all
cm,b_temp,volu
vgen,,all,,,rad_drill,,,,,1
csys,1
vgen,,all,,,,n,,,,1
mat,3
esize,elemsize/2,0,
vsweep,all
vgen,4,all,,,,n*2,,,0
vclear,b_temp
vdele,b_temp
vgen,,all,,,,-n*2,,,,1
cm,bvolu,volu
alls
psmesh,1,b100,,volu,95,0,z,0,,,,b_100
psmesh,2,b200,,volu,116,0,z,0,,,,b_200
psmesh,3,b300,,volu,109,0,z,0,,,,b_300
cmsel,r,shaft
aslv
asel,r,loc,z,0
cm,as_con,area
alls
cmsel,r,pan
aslv
asel,r,loc,z,0
cm,ap_tar,area
alls
/com,contact pair creation -start
mat,4
r,3
real,3
et,3,170
et,4,174
keyopt,4,9,0
cmsel,s,ap_tar
type,3
nsla,s,1
esln,s,0
esurf,all
alls
cmsel,s,as_con
type,4
nsla,s,1
esln,s,0
esurf,all
allsel
cmsel,r,bvolu
aslv
asel,r,loc,z,flang_h1
cm,ab1_con,area
alls
cmsel,r,shaft
aslv
asel,r,loc,z,flang_h1
csys,1
asel,r,loc,x,fh_rad,rad_o
csys,0
cm,abs_tar,area
alls
/com, contact pair creation - start
mat,4
r,4
real,4
et,5,170
et,6,174
keyopt,5,9,0
cmsel,s,abs_tar
type,5
nsla,s,1
esln,s,0
esurf,all
alls
cmsel,s,ab1_con
type,6
nsla,s,1
esln,s,0
esurf,all
allsel
/com, contact pair creation - end
cmsel,r,bvolu
aslv
asel,r,loc,z,-flang_h2
cm,ab2_con,area
alls
cmsel,r,pan
aslv
asel,r,loc,z,-flang_h2
cm,abp_tar,area
alls
/com, contact pair creation - start
mat,4
r,5
real,5
et,7,170
et,8,174
keyopt,7,9,0
cmsel,s,abp_tar
type,7
nsla,s,1
esln,s,0
esurf,all
alls
cmsel,s,ab2_con
type,8
nsla,s,1
esln,s,0
esurf,all
allsel
/com, contact pair creation - end

/SOL
alls
lsclear,all
antype,0
autots,1
nsubst,25,100,1
outres,all,all
eqslv,pcg,1e-8
time,15
csys,1
asel,s,loc,y,0
da,all,symm
alls
asel,s,loc,y,90
da,all,symm
alls
sload,1,9,lock,disp,disp_b,1,2
sload,2,9,lock,disp,disp_b,1,2
sload,3,9,lock,disp,disp_b,1,2
cmsel,s,pan
nslv,r,1
nsel,r,loc,z,0
csys,1
nsel,r,loc,x,u_r,r_out
d,all,all
alls
***ctran
lswrite,1
lswrite,2
nsel,s,loc,z,zmax-1,zmax+10
*get,nload,node,,count
f,all,fz,f_ext/nload/4
alls
lswrite,3
lssolve,1,3,1
/post
set,last
plnsol,s,eqv,0,1

[ 本帖最后由 jinqinghai 于 2006-10-2 14:09 编辑 ]

jinqinghai

image of stress of SEQv

sunflower

GOOD ^-^

jtxqhm

版主辛苦了!
后来忙于找工作,写论文一直没动,真是不好意思!

wuhaiyu9713

我的版本是10.0，我输入你的命令流，结果运行不起来
！！！！！！！！
不知楼主是怎么弄的

mikehe

呵呵 我的是9.0的进去错误也是无穷 晕倒
也不知道是怎么回事？

ultrastar

其实命令流本身是没错误的

错就错在 论坛的关键字屏蔽掉了，我下载了份PDF的和此命令流对照过了

论坛屏蔽了 V*S*B*W等等命令，所以会出现错误。

稍后我会上传份更正过的TXT供大家阅览

heminbox

我对加载的部分不是太明白，nsubst,25,100,1；
time,15
asel,s,loc,y,0
da,all,symm

sload,1,9,lock,disp,disp_b,1,2
sload,2,9,lock,disp,disp_b,1,2
sload,3,9,lock,disp,disp_b,1,2
能否解释下？为何要这样处理

yangmeijun

在第一步载荷步中对螺栓预紧，即b_100预紧单元中，施加预紧位移为disp_bmm

第二步载荷 步锁紧这个位移最终产生螺栓预紧作用。

sunnygu

换了***,怎么还是错误阿。什么时候当附件穿一下啊

liuliukan

我上传了这个螺栓分析命令了，请斑竹加点分，附件里已完整修改原书中觉得有些不妥的地方，并全部添加了注释。在ansys9， ansys10中测试通过。
注意：因论坛限制，附件中已将文件后缀名改为rar了，请下载后将文件后缀名改回dat或者txt.

hehe369

我做过
现在说边界条件没定义好
怎么办?

月陨星沉

谢谢了，若有问题再请教啊……