抱歉楼主，是我没有太理解你的本意。
ANSYS遇到接触非线性分析不收敛，可以结合output中的信息对模型进行检查、调试，通常的原因有以下几点：1.模型缺少足够的约束条件；2.接触刚度过大；3.接触对设置不正确；4.求解子步设置过小；5.未定义摩擦系数。另外，施加位移输入相比力输入要更加容易收敛。
按照以上经验，对楼主的代码进行了适当修改。1.将两端力拉伸改为一端约束，另一端位移输入；2.默认的接触刚度系数减小。在自家电脑上运行收敛通过，至少证明接触对的建立是没有问题。若一定要求输入条件是力，可能还需要楼主进一步调试。修改后收敛的代码如下：
/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
r, 4, 0, 0, 0.001, 0.1, 0, 0 ! contact pair 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
r, 5, 0, 0, 0.001, 0.1, 0, 0 ! contact pair 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
r, 6, 0, 0, 0.001, 0.1, 0, 0 ! contact pair 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
r, 7, 0, 0, 0.001, 0.1, 0, 0 ! contact pair 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
r, 8, 0, 0, 0.001, 0.1, 0, 0 ! contact pair 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, all, uy, 0
d, all, uz, 0
d, all, ux, 1.0e-5
! 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,2,ux,all
d, all, uy, 0
d, all, uz, 0
d, all, ux, 0
! nsel,s,loc,x,-76/1000,-76/1000
! nsel,r,loc,y
! nsel,r,loc,z
! f,all,fx,pload2
allsel
solve
finish |
回复此楼
5+3
物理学调剂
+4