ansys接触分析不收敛，希望高手指点指点

qyt913

/filename,jiechufenxi                !定义文件名
/title, aim to verify solid model    !定义标题
/prep7
bosongbi=0.48                         !设置泊松比为全局变量
thick=0.2e-3                             !设置表皮层厚度为全局变量
ET,1,PLANE183                         !定义单元
KEYOPT,1,3,2                          !单元选项为平面应变
MP,EX,1,1.36e5   
MP,PRXY,1,bosongbi   
MP,EX,2,8e4
MP,PRXY,2,bosongbi
MP,EX,3,3.4e4  
MP,PRXY,3,bosongbi
MP,EX,4,2e11
MP,PRXY,4,0.3
                        !四种材料定义完毕
!建立模型
CSWPLA,11,1,0.75,1,                   !创建局部坐标系
K,1,8e-3,,,                           !创建关键点
K,2,-8e-3,,,  
L,       2,       1
K,3,thick-8e-3,,,
K,4,-thick+8e-3,,,
K,5,-thick+7e-3,,,
K,6,thick-7e-3,,,
L,       3,       4
L,       6,       5
CSYS,0  
lsymm,y,all
CSWPLA,12,1,4/5.33,1,                !创建局部坐标系
k,101,2.65e-3,0
k,102,-2.65e-3,0
l,101,102
csys,0
lsymm,y,7                            !建立指骨   
lgen,,7,8,,,1e-3,,,,1                !向上平移1mm
nummrg,kp
WPCSYS,-1,0
RECTNG,-7e-3,7e-3,-6e-3,-7e-3
al,1,4,2,5
al,3,2,5,6
al,3,6
al,7,8
asba,4,5
aglue,2,3,6
wpoff,0,1e-3,0
wpro,,,-90  
asbw,all
wpro,,-90
asbw,all
wpcsys,-1,0
alls
asel,s,loc,x,-10,0
adele,all,,,1
alls
numcmp,all
lcomb,16,22
lcomb,11,24
lcomb,10,23               !模型创建完毕
!划分网格  
asel,s,area,,1,6,5
aatt,1,,1,0
asel,s,area,,3,7,4
aatt,2,,1,0
asel,s,area,,4,5,1
aatt,3,,1,0
asel,s,area,,2
aatt,4,,1,0             !赋予材料属性
lesize,12,,,4
lesize,13,,,4
lesize,3,,,4
lesize,1,,,40
lesize,2,,,40
lesize,10,,,60
lesize,11,,,60
lesize,15,,,5
lesize,6,,,5
lesize,14,,,5
lesize,16,,,60
lesize,4,,,40
lesize,19,,,15
lesize,21,,,15
lesize,20,,,15
lesize,17,,,40
lesize,18,,,60
lesize,5,,,4
lesize,9,,,4
lesize,7,,,20
lesize,8,,,20
alls
mshape,0
mshkey,1
amesh,all            !网格划分完毕
!设置接触对
/COM, CONTACT PAIR CREATION - START
  
MAT,1   
R,3                     !定义接触实常数
REAL,3  
ET,2,169                !定义接触单元类型
ET,3,172
KEYOPT,3,5,1
KEYOPT,3,10,2   
R,3,0,0,0.1, 0,1e-008,0
! Generate the target surface   
LSEL,S,,,8  
CM,_TARGET,LINE
TYPE,2  
NSLL,S,1
ESLN,S,0
ESURF,all   
! Generate the contact surface  
LSEL,S,,,10
CM,_CONTACT,LINE
TYPE,3  
NSLL,S,1
ESLN,S,0
ESURF,all
cmdel,_TARGET
cmdel,_CONTACT
ALLSEL,all  
Eplot
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT
finish                       !接触对创建完毕
!进入求解器
/solu
antype,0
nlgeom,1
nsubst,100,200,20
outres,all,all
eqslv,spar
lnsrch,1
pred,on,,on
time,2
NROPT,FULL, ,OFF
!施加约束
allsel,all
lplot
dl,9,,symm
dl,12,,symm
dl,13,,symm
dl,14,,symm
dl,15,,symm
dl,19,,symm
dl,20,,symm
dl,17,,all,0
dl,18,,all,0
dl,7,,uy,2e-3
solve                     !求解

qyt913

我从网格质量上以及接触对属性上进行修改，但总是没有成功，结果不理想。
问题主要集中在：
1.一个单元或多个单元过度扭曲。
2.结果不收敛，到1.63s左右就迭代不下去了。
希望高手指点指点，困惑我很久了！

ferriswanglei

位移太大，单元扭曲不可避免，将网格再细一点试试看，如果不行，可采取Advanced Analysis Techniques Guide中的方法，见
http://forum.simwe.com/thread-967491-1-1.html

yuqianjun123

我表示关注啊...

lixiaodong163

可以改变一下接触刚度和初始穿透度试试

白云青溪

好东西，感谢楼主

haozhiqiang888

学习了  非常感谢！

cc332255

同样关注着这个问题