抱歉樓主，是我沒(méi)有太理解你的本意。
ANSYS遇到接觸非線性分析不收斂，可以結(jié)合output中的信息對(duì)模型進(jìn)行檢查、調(diào)試，通常的原因有以下幾點(diǎn)：1.模型缺少足夠的約束條件；2.接觸剛度過(guò)大；3.接觸對(duì)設(shè)置不正確；4.求解子步設(shè)置過(guò)小；5.未定義摩擦系數(shù)。另外，施加位移輸入相比力輸入要更加容易收斂。
按照以上經(jīng)驗(yàn)，對(duì)樓主的代碼進(jìn)行了適當(dāng)修改。1.將兩端力拉伸改為一端約束，另一端位移輸入；2.默認(rèn)的接觸剛度系數(shù)減小。在自家電腦上運(yùn)行收斂通過(guò)，至少證明接觸對(duì)的建立是沒(méi)有問(wèn)題。若一定要求輸入條件是力，可能還需要樓主進(jìn)一步調(diào)試。修改后收斂的代碼如下：
/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
!定義參數(shù)
l1=100/1000 !試樣長(zhǎng)度
l2=34/1000 !夾頭長(zhǎng)度
w1=24/1000 !試樣寬度
w2=33.8/1000 !夾頭寬度
h=1.5/1000 !試樣厚度
l3=48/1000 !試樣交疊部分長(zhǎng)度
r=4.9/1000 !連接圓弧半徑
r1=6.25/2/1000 !圓臺(tái)下表面半徑
r2=8.75/2/1000 !圓臺(tái)上表面半徑
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
!劃分網(wǎng)格(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
!劃分網(wǎng)格(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
!劃分網(wǎng)格(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
!建立接觸對(duì)
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 |
回復(fù)此樓
