最近做一个接触模拟，锚杆拉拔，但是始终不收敛，请求各位帮我看下指导一下我谢谢了
浏览：6012 回答：1
帮我看下吧，我努力了很久，闭门造车什么都没弄出来
finish$/clear$/prep7
k,6,10,1000
k,7,30,1000
k,8,32,1000
k,9,32,2000
k,10,30,2000
k,11,10,2000
l,6,7
l,7,8
l,8,9
l,9,10
l,10,11
l,11,6
l,7,10
al,1,2,3,4,5,6
asbl,1,7,sepo
k,15,0,0
k,16,1000,0
k,17,1000,2000
k,18,0,2000
k,19,0,1000
k,20,632,2000
k,21,32,2014
k,22,632,2014
k,23,0,999.5
k,24,32,999.5
l,15,16
l,16,17
l,17,20
l,20,22
l,22,21
l,21,13
l,20,13
l,4,18
l,18,19
l,1,19
l,12,24
l,23,24
l,23,15
ldele,4,6,1
al,1,2,3,7,13,19,20,21
al,16,17,18,10
!土体参数
E_1=30E6
U_1=0.3
DENS_1=1900
C_1=30E3
FI_1=25
DFI_1=19
MP,EX,1,E_1
MP,PRXY,1,U_1
MP,DENS,1,DENS_1
TB,DP,1
TBDATA,1,C_1,FI_1,DFI_1
!锚固体参数
E_2=26E9
U_2=0.2
DENS_2=2500
C_2=1.2E6
FI_2=32
DFI_2=8
MP,EX,2,E_2
MP,PRXY,2,U_2
MP,DENS,2,DENS_2
TB,DP,2
TBMODIF,1,1,C_2
TBMODIF,1,2,FI_2
TBMODIF,1,3,DFI_2
!锚杆参数
E_4=17E10
U_4=0.3
DENS_4=7900
MP,EX,4,E_4
MP,PRXY,4,U_4
MP,DENS,4,DENS_4
!界面单元
E_3=10E7
U_3=0.2
DENS_3=2200
C_3=40E3
FI_3=25
DFI_3=23
MP,EX,3,E_3
MP,PRXY,3,U_3
MP,DENS,3,DENS_3
TB,DP,3
TBMODIF,1,1,C_3
TBMODIF,1,2,FI_3
TBMODIF,1,3,DFI_3
ET,1,42
keyopt,1,3,1
ET,2,169
ET,3,171
mp,mu,5,0.1
r,1,,,,,,,,,,,,,35e2,,,,,,,,1.5,8
keyopt,3,1,0
keyopt,3,2,0
keyopt,3,4,0
keyopt,3,5,0
keyopt,3,7,0
keyopt,3,8,0
keyopt,3,9,0
keyopt,3,10,0
keyopt,3,11,0
keyopt,3,12,2
lsel,s,line,,11,15,2
LESIZE,ALL,,,250
lsel,s,line,,10,17,7
LESIZE,ALL,,,250
lsel,s,line,,12,14,2
LESIZE,ALL,,,1
ASEL,S,AREA,,3
AATT,3,,1
mshkey,1
mshape,0
AMESH,ALL
lsel,s,line,,8,9,1
LESIZE,ALL,,,10
lsel,s,line,,16,18,2
LESIZE,ALL,,,4
ASEL,S,AREA,,2
AATT,2,,1
mshkey,0
mshape,0
AMESH,ALL
ASEL,S,AREA,,4
AATT,4,,1
mshkey,0
mshape,0
AMESH,ALL
lsel,s,line,,20
LESIZE,ALL,,,6
lsel,s,line,,19
LESIZE,ALL,,,1
lsel,s,line,,1,2,1
LESIZE,ALL,,,30
lsel,s,line,,21
LESIZE,ALL,,,50,0.5
lsel,s,line,,3
LESIZE,ALL,,,10
lsel,s,line,,7
LESIZE,ALL,,,180,0.3
ASEL,S,AREA,,1
AATT,1,,1
mshkey,0
mshape,0
AMESH,ALL
!创建目标面（土钉/锚杆）
lsel,s,line,,11
nsll,s,1
esln,s,0
type,2
real,1
mat,5
esurf
allsel,all,all
!创建接触面（界面层）
lsel,s,line,,15
nsll,s,1
esln,s,0
type,3
real,1
mat,5
esurf
esurf,all,reverse
allsel,all,all
!施加约束及荷载
dl,1,,uy
dl,2,,ux
dl,21,,ux
dl,17,,ux
allsel,all
LSEL,S,LINE,,16
SFL,ALL,PRES,-10
allsel,all,all
lsel,s,line,,15
nsll,s,1
esln,s,0
esurf,all,reverse
allsel,all,all
/solu
antype,static
nlgeom,on
autots,on
time,1
NSUBST,100
OUTRES,ALL,ALL
solve
finish$/clear$/prep7
k,6,10,1000
k,7,30,1000
k,8,32,1000
k,9,32,2000
k,10,30,2000
k,11,10,2000
l,6,7
l,7,8
l,8,9
l,9,10
l,10,11
l,11,6
l,7,10
al,1,2,3,4,5,6
asbl,1,7,sepo
k,15,0,0
k,16,1000,0
k,17,1000,2000
k,18,0,2000
k,19,0,1000
k,20,632,2000
k,21,32,2014
k,22,632,2014
k,23,0,999.5
k,24,32,999.5
l,15,16
l,16,17
l,17,20
l,20,22
l,22,21
l,21,13
l,20,13
l,4,18
l,18,19
l,1,19
l,12,24
l,23,24
l,23,15
ldele,4,6,1
al,1,2,3,7,13,19,20,21
al,16,17,18,10
!土体参数
E_1=30E6
U_1=0.3
DENS_1=1900
C_1=30E3
FI_1=25
DFI_1=19
MP,EX,1,E_1
MP,PRXY,1,U_1
MP,DENS,1,DENS_1
TB,DP,1
TBDATA,1,C_1,FI_1,DFI_1
!锚固体参数
E_2=26E9
U_2=0.2
DENS_2=2500
C_2=1.2E6
FI_2=32
DFI_2=8
MP,EX,2,E_2
MP,PRXY,2,U_2
MP,DENS,2,DENS_2
TB,DP,2
TBMODIF,1,1,C_2
TBMODIF,1,2,FI_2
TBMODIF,1,3,DFI_2
!锚杆参数
E_4=17E10
U_4=0.3
DENS_4=7900
MP,EX,4,E_4
MP,PRXY,4,U_4
MP,DENS,4,DENS_4
!界面单元
E_3=10E7
U_3=0.2
DENS_3=2200
C_3=40E3
FI_3=25
DFI_3=23
MP,EX,3,E_3
MP,PRXY,3,U_3
MP,DENS,3,DENS_3
TB,DP,3
TBMODIF,1,1,C_3
TBMODIF,1,2,FI_3
TBMODIF,1,3,DFI_3
ET,1,42
keyopt,1,3,1
ET,2,169
ET,3,171
mp,mu,5,0.1
r,1,,,,,,,,,,,,,35e2,,,,,,,,1.5,8
keyopt,3,1,0
keyopt,3,2,0
keyopt,3,4,0
keyopt,3,5,0
keyopt,3,7,0
keyopt,3,8,0
keyopt,3,9,0
keyopt,3,10,0
keyopt,3,11,0
keyopt,3,12,2
lsel,s,line,,11,15,2
LESIZE,ALL,,,250
lsel,s,line,,10,17,7
LESIZE,ALL,,,250
lsel,s,line,,12,14,2
LESIZE,ALL,,,1
ASEL,S,AREA,,3
AATT,3,,1
mshkey,1
mshape,0
AMESH,ALL
lsel,s,line,,8,9,1
LESIZE,ALL,,,10
lsel,s,line,,16,18,2
LESIZE,ALL,,,4
ASEL,S,AREA,,2
AATT,2,,1
mshkey,0
mshape,0
AMESH,ALL
ASEL,S,AREA,,4
AATT,4,,1
mshkey,0
mshape,0
AMESH,ALL
lsel,s,line,,20
LESIZE,ALL,,,6
lsel,s,line,,19
LESIZE,ALL,,,1
lsel,s,line,,1,2,1
LESIZE,ALL,,,30
lsel,s,line,,21
LESIZE,ALL,,,50,0.5
lsel,s,line,,3
LESIZE,ALL,,,10
lsel,s,line,,7
LESIZE,ALL,,,180,0.3
ASEL,S,AREA,,1
AATT,1,,1
mshkey,0
mshape,0
AMESH,ALL
!创建目标面（土钉/锚杆）
lsel,s,line,,11
nsll,s,1
esln,s,0
type,2
real,1
mat,5
esurf
allsel,all,all
!创建接触面（界面层）
lsel,s,line,,15
nsll,s,1
esln,s,0
type,3
real,1
mat,5
esurf
esurf,all,reverse
allsel,all,all
!施加约束及荷载
dl,1,,uy
dl,2,,ux
dl,21,,ux
dl,17,,ux
allsel,all
LSEL,S,LINE,,16
SFL,ALL,PRES,-10
allsel,all,all
lsel,s,line,,15
nsll,s,1
esln,s,0
esurf,all,reverse
allsel,all,all
/solu
antype,static
nlgeom,on
autots,on
time,1
NSUBST,100
OUTRES,ALL,ALL
solve
kakaqun 
