24小时热门版块排行榜

返回列表

11050209

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 938.2
散金: 100
帖子: 321
在线: 68.3小时
虫号: 3896498
注册: 2015-05-28
专业: 金属材料的凝固与结晶学

[求助] 接触问题求助已有1人参与

两块板通过螺栓连接，拉伸问题的求助，搞不懂什么原因导致算出来的应力特别小！由于第一次接触到接触分析，所以不是很懂，我觉得问题可能出在接触这一块，不知道大神能不能帮看下，还有一个可能就是我施加的约束不够?如果是约束的问题又该怎么加约束呢
附上命令：
/prep7
/TITLE, lashen
!定义单元类型及材料属性
et,1,solid185
et,2,plane182
mp,ex,1,190e6
mp,ex,2,184e6
mp,ex,3,69e9
mp,prxy,1,0.34
mp,prxy,2,0.3
mp,prxy,3,0.33
!定义参数
l1=100/1000 !试样长度
l2=34/1000 !夹头长度
w1=24/1000 !试样宽度
w2=33.8/1000 !夹头宽度
h=1.5/1000 !试样厚度
l3=48/1000 !试样交叠部分长度
r=4.9/1000    !连接圆弧半径
r1=6.25/2/1000    !圆台下表面半径
r2=8.75/2/1000    !圆台上表面半径
r3=4/1000    !螺母半径
pload1=100000000*2*h*w2
pload2=-pload1
!建立几何模型(TA32)
block,-l3/2,-l3/4,0,w1/2,0,h
block,-l3/4,0,0,w1/2,0,h
block,-l3/4,0,-w1/2,0,0,h
block,-l3/2,-l3/4,-w1/2,0,0,h
block,0,l3/4,0,w1/2,0,h
block,l3/4,l3/2,0,w1/2,0,h
block,l3/4,l3/2,-w1/2,0,0,h
block,0,l3/4,-w1/2,0,0,h
wpoff,-12/1000,0,0
cone,r1,r2,0,h
wpoff,24/1000,0,0
cone,r1,r2,0,h
vsel,s,,,1,4,1
allsel
vsbv,all,9
allsel
vsel,s,,,5,8,1
allsel
vsbv,all,10
numcmp,all
CSYS,0
WPAVE,0,0,0
CSYS,0
block,-37.1/1000,-24/1000,-w1/2,w1/2,0,h
k,73,-42/1000,w1/2+r,h
k,74,-76/1000,w1/2+r,h
k,75,-76/1000,-(w1/2+r),h
k,76,-42/1000,-(w1/2+r),h
larc,72,73,54,r
l,73,74
l,74,75
l,75,76
larc,76,69,54,r
lsel,s,,,125,129,1
lsel,a,,,120
al,all
allsel
vdrag,63,,,,,,121
block,-76/1000,-42/1000,-w2/2,w2/2,-h,0
numcmp,all
nummrg,all
a,73,82,79,76
vsba,10,8
numcmp,all
allsel
vglue,all
!建立几何模型(7050)
block,-l3/2,-l3/4,0,w1/2,-h,0
block,-l3/4,0,0,w1/2,-h,0
block,-l3/4,0,-w1/2,0,-h,0
block,-l3/2,-l3/4,-w1/2,0,-h,0
block,0,l3/4,0,w1/2,-h,0
block,l3/4,l3/2,0,w1/2,-h,0
block,l3/4,l3/2,-w1/2,0,-h,0
block,0,l3/4,-w1/2,0,-h,0
wpoff,-12/1000,0,0
cylind,0,r1,0,-h
wpoff,24/1000,0,0
cylind,0,r1,0,-h
vsel,s,,,14,16,1
vsel,a,,,9
allsel
vsbv,all,21
vsel,s,,,17,20,1
allsel
vsbv,all,22
CSYS,0
WPAVE,0,0,0
CSYS,0
block,24/1000,37.1/1000,-w1/2,w1/2,-h,0

k,139,42/1000,w1/2+r,0
k,140,76/1000,w1/2+r,0
k,141,76/1000,-(w1/2+r),0
k,142,42/1000,-(w1/2+r),0

larc,85,139,19,r
l,139,140
l,140,141
l,141,142
larc,82,142,24,r
lsel,s,,,186,187,1
lsel,a,,,192,193,1
lsel,a,,,167,197,30
al,all
vdrag,99,,,,,,178
block,42/1000,76/1000,-w2/2,w2/2,0,h
a,117,121,114,100
vsba,18,143
vsel,s,,,14,17,1
vsel,a,,,23,26,1
vsel,a,,,19,21,1
vsel,a,,,9
vglue,all
numcmp,all
!建立几何模型(1s627)
wpoff,-12/1000,0,0
cone,r1,r2,0,h
cylind,0,r1,0,-(h+h)
cylind,r1,r3,-h,-(h+h)
allsel
vsel,s,,,25,27,1
vglue,all
wpoff,24/1000,0,0
cone,r1,r2,0,h
cylind,0,r1,0,-(h+h)
cylind,r1,r3,-h,-(h+h)
allsel
vsel,s,,,27,30,3
vsel,a,,,26
vglue,all
CSYS,0
WPAVE,0,0,0
CSYS,0
allsel
!划分网格(TA32)
MSHAPE,0,2D
MSHKEY,1
mat,1
lsel,s,,,7,8
lesize,all,,,6
lsel,s,,,74,76,2
lesize,all,,,4
ALLSEL,ALL
AMAP,39,4,33,32,5
lsel,s,,,17,18,1
lesize,all,,,6
AMAP,43,5,32,31,9
allsel
AMAP,47,9,31,34,13
AMAP,50,13,34,33,4
lsel,s,,,4,6,2
lesize,all,,,4
ALLSEL,ALL
AMAP,7,9,41,40,18
AMAP,14,18,40,39,21
AMAP,22,21,39,42,25
AMAP,34,25,42,41,9
lsel,s,,,35
lesize,all,,,3
vsel,s,,,1,8,1
vsweep,all
asel,s,,,39,43,4
asel,a,,,47,50,3
asel,a,,,7,14,7
asel,a,,,22,34,12
aclear,all
allsel
vsweep,13
lsel,s,,,93
lesize,all,,,2,2
allsel
lsel,s,,,97
lesize,all,,,2,0.5
vsweep,11
lsel,s,,,94,96,2
lesize,all,,,17
vsweep,12
vsweep,10
allsel
!划分网格(7050)
mat,3
AMAP,96,62,101,100,63
AMAP,121,63,100,99,76
AMAP,132,76,99,102,68
AMAP,129,68,102,101,62
AMAP,73,76,109,108,77
AMAP,114,77,108,107,83
AMAP,117,83,107,110,88
AMAP,111,88,110,109,76
lsel,s,,,125
lesize,all,,,3
vsel,s,,,17,21,1
vsel,a,,,23,24,1
vsel,a,,,9
vsweep,all
asel,s,,,73,96,23
asel,a,,,111,117,3
asel,a,,,121,129,8
asel,a,,,132
aclear,all
allsel
lsel,s,,,218,220,2
lesize,all,,,7
vsweep,22
lsel,s,,,199,200,1
lesize,all,,,2,2
vsweep,14
lsel,s,,,185,187,2
lesize,all,,,17
vsweep,16
vsweep,15
!划分网格(1s627)
mat,2
lsel,s,,,231,234,1
lesize,all,,,12
vsweep,25
lsel,s,,,259,262,1
lesize,all,,,12
vsweep,28
lsel,s,,,255,258,1
lesize,all,,,12
vsweep,29
lsel,s,,,241,244,1
LESIZE,all, , ,12
vsweep,26
lsel,s,,,291,294,1
lesize,all,,,12
vsweep,31
lsel,s,,,287,290,1
lesize,all,,,12
vsweep,32
allsel
!建立接触对
mp,mu,4,0.3
mat,4
r,4
real,4
et,3,targe170
et,4,conta173
!generate the target surface
asel,s,,,135,136,1
cm,left_up_bolt,area
cmsel,s,left_up_bolt
type,3
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,37,45,4
asel,a,,,48
cm,left_up_plate,area
cmsel,s,left_up_plate
type,4
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,5
real,5
et,5,targe170
et,6,conta173
asel,s,,,147,148,1
cm,left_down_bolt,area
cmsel,s,left_down_bolt
type,5
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,94,97,3
asel,a,,,99,100,1
cm,left_down_plate,area
cmsel,s,left_down_plate,area
type,6
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,6
real,6
et,7,targe170
et,8,conta173
asel,s,,,140,141,1
cm,right_up_bolt,area
cmsel,s,right_up_bolt
type,7
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,3,9,6
asel,a,,,19,32,13
cm,right_up_plate,area
cmsel,s,right_up_plate
type,8
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,7
real,7
et,9,targe170
et,10,conta173
asel,s,,,160,161,1
cm,right_down_bolt,area
cmsel,s,right_down_bolt
type,9
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,56,74,18
asel,a,,,77,87,10
cm,right_down_plate,area
cmsel,s,right_down_plate
type,10
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,8
real,8
et,11,targe170
et,12,conta173
asel,s,,,6,13,7
asel,a,,,20,33,13
asel,a,,,38,46,4
asel,a,,,49
cm,left_plate,area
cmsel,s,left_plate
type,11
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,129,132,3
asel,a,,,96,121,25
asel,a,,,111,117,6
asel,a,,,73,114,41
cm,right_plate,area
cmsel,s,right_plate
type,12
nsla,s,1
esln,s,0
esurf,all
allsel
/sol
antype,static
nlgeom,on
autots,on
lnsrcch,on
pred,on
outres,all,all
time,1
nsubst,80,1000,20
nsel,s,loc,x,76/1000,76/1000
cp,1,ux,all
d,7377,uy,0
d,7377,uz,0
nsel,s,loc,x,76/1000,76/1000
nsel,r,loc,y
nsel,r,loc,z
f,all,fx,pload1
allsel
nsel,s,loc,x,-76/1000,-76/1000
cp,1,ux,all
nsel,s,loc,x,-76/1000,-76/1000
nsel,r,loc,y
nsel,r,loc,z
f,all,fx,pload2
allsel
solve
finish

回复此楼