求助快开结构的接触分析问题？

tiexie1986

finish
/clear
!!!!!!!!!!环境设置!!!!!!!!!!!!
*afun,deg  
!设定结构基本参数,单位N,mm,MPa
*set,P,51                       
!结构尺寸参数
*set,r1,100
*set,r2,120
*set,r3,140
*set,r4,170
*set,r5,143
*set,r6,30
*set,r7,173
!
*set,t1,140
*set,t2,120
*set,t3,140
*set,t4,160
*set,t5,180
*set,t6,185
*set,t7,190
*set,ang,20
*set,ang1,4
*set,ang2,8
*set,ang3,12
*set,ang4,16
!
exx1=2.0e5                   !封头
nuxy1=0.3
exx2=2.0e5                   !壳体及上法兰
nuxy2=0.3
exx3=2.0e5                   !卡箍
nuxy3=0.3
*set,M1,2E4
*set,M2,5
*set,P1,60
*set,P3,65
*set,N1,1e1
*set,N2,2e1
*set,N3,5e1
*set,N4,1e2
*set,N5,2e2
*set,N6,5e2
*set,N7,1e3
*set,N8,2e3
*set,N9,5e3
*set,N10,1e4
*set,N11,2e4
*set,N12,5e4
*set,N13,1e5
*set,N14,2e5
*set,N15,5e5
*set,N16,1e6
*set,S1,4000
*set,S2,2828
*set,S3,1897
*set,S4,1414
*set,S5,1069
*set,S6,724
*set,S7,572
*set,S8,441
*set,S9,331
*set,S10,262
*set,S11,214
*set,S12,159
*set,S13,138
*set,S14,114
*set,S15,93.1
*set,S16,86.2
/prep7                           
!!!!!!定义单元及材料属性!!!!!!
et,1,solid92                  
et,2,TARGE170                  
et,3,CONTA174
keyopt,3,9,1
keyopt,3,10,1
keyopt,3,12,0
keyopt,3,7,1
keyopt,3,2,0
r,1                              
!定义封头材料1
mp,ex,1,exx1                  
mp,nuxy,1,nuxy1                 
!定义壳体及上法兰材料2
mp,ex,2,exx2                 
mp,nuxy,2,nuxy2
!定义卡箍材料3
mp,ex,3,exx3               
mp,nuxy,3,nuxy3   
!!!!!!创建模型!!!!!!
!define kp
k,1,0,t7+t5-t4+t3-t2
k,2,r3,t7+t5-t4+t3-t2
k,3,r3,t7+t5-t4
k,4,r3,t7
k,5,r3,t6
k,6,r1,t6
k,7,r1,t7
k,8,0,t7
k,9,r1,t6
k,10,r3,t6
k,11,r3,t5
k,12,r3,t4
k,13,r3,t1
k,14,r2,t2
k,15,r2,0
k,16,r1,0
k,17,r5+r4-r3,t7+t5-t4+t3-t2
k,18,r5+r4-r3+r6,t7+t5-t4+t3-t2
k,19,r5+r4-r3+r6,t3
k,20,r5+r4-r3,t3
k,21,r7,t4
k,22,r7,t7+t5-t4
l,1,2
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7
l,7,8
l,1,8
l,9,10
l,10,11
l,11,12
l,12,13
l,13,14
l,14,15
l,15,16
l,16,9
l,17,18
l,18,19
l,19,20
l,21,20
l,21,22
l,22,17
a,1,2,3,4,5,6,7,8
a,9,10,11,12,13,14,15,16
a,17,18,19,20,21,22
allsel
vrotat,all,,,,,,1,8,ang

allsel
vgen,2,all,,,0,-t7,0,,,1                 !移动体到一定位置
CSYS,5
k,51,r4,-ang3,0
k,52,r3+0.95*(r4-r3),-ang3-0.125*(ang4-ang3),0
k,53,r3+0.88*(r4-r3),-ang3-0.25*(ang4-ang3),0
k,54,r3+0.76*(r4-r3),-ang3-0.375*(ang4-ang3),0
k,55,r3+0.58*(r4-r3),-ang3-0.5*(ang4-ang3),0
k,56,r3+0.35*(r4-r3),-ang3-0.625*(ang4-ang3),0
k,57,r3+0.16*(r4-r3),-ang3-0.75*(ang4-ang3),0
k,58,r3+0.08*(r4-r3),-ang3-0.875*(ang4-ang3),0
k,59,r3,-ang4,0
k,61,r4,-ang2,0
k,62,r3+0.95*(r4-r3),-ang2-(ang1-ang2)/8,0
k,63,r3+0.88*(r4-r3),-ang2-(ang1-ang2)/4,0
k,64,r3+0.76*(r4-r3),-ang2-3*(ang1-ang2)/8,0
k,65,r3+0.58*(r4-r3),-ang2-(ang1-ang2)/2,0
k,66,r3+0.35*(r4-r3),-ang2-5*(ang1-ang2)/8,0
k,67,r3+0.16*(r4-r3),-ang2-6*(ang1-ang2)/8,0
k,68,r3+0.08*(r4-r3),-ang2-7*(ang1-ang2)/8,0
k,69,r3,-ang1,0
ksel,s,,,51,59
bsplin,all
allsel
ksel,s,,,61,69
bsplin,all
allsel
l,51,61
allsel
l,59,69
allsel
al,65,67,64,66
VEXT, 28,,, 0,0,t4-t3                          !拉伸面
allsel
vgen,2,4,,,0,0,-t7+t4,,,0
wprot,0,0,-ang/2
CSYS,4
LSYMM, x, 68, 70,2,,,, 0
LGEN, 2, 88, 89, , r4+r5, 0, 0,,, 1                          
LPTN, 88,89,63
LDELE, 93, 95,2,,,  1
wprot,0,0,ang/2
CSYS,5
l,73,75,2
al,92,88,94,96
CSYS,0
VEXT, 40,,,0,t3-t2,0
vgen,2,6,,,0,-t7+t4-t5+t3,0,,,0
VGLUE, 1,4
VGLUE, 2,5
VGLUE, 3,6,7
allsel
!划分网格
mshkey,0
mshape,1,3d
VMESH,8,9,1
allsel
mat,2
VMESH,1,4,3
allsel
mat,3
VMESH,2,10,8
vmesh,5
!创建接触
!创建封头与卡箍的接触
CSYS,0
allsel
asel,s,area,,40
type,2
real,1
r,,,,0.01,0.1
rmore,,,250/1.732,,,1
nsla,s,1
esln,s,0
esurf,all
asel,s,area,,53
type,3
nsla,s,1
esln,s,0
esurf,all
allsel
!创建筒体与卡箍的接触
et,4,TARGE170                  
et,5,CONTA174
r,2
allsel
asel,s,area,,47
type,4
real,2
r,,,,0.01,0.1
rmore,,,250/1.732,,,1
nsla,s,1
esln,s,0
esurf,all
asel,s,area,,5
type,5
nsla,s,1
esln,s,0
esurf,all
allsel
!创建筒体与封头的接触
et,6,TARGE170                  
et,7,CONTA174
r,3
allsel
asel,s,area,,12
type,6
real,3
nsla,s,1
esln,s,0
esurf,all
asel,s,area,,8
type,7
nsla,s,1
esln,s,0
esurf,all
!!!!!!!!接触创建完毕!!!!!!!!
allsel
finish
/sol
lsclear,all
nropt,unsym
pred,on
lnsrch,on
autots,1
antype,static
autots,1
nsubst,25,100,1
outres,all,all
wpsty,defa
time,10
allsel
lswrite,1
!time,20
nropt,unsym
asel,s,area,,18
da,all,uy
allsel
asel,s,area,,1,3,1
da,all,symm
allsel
asel,s,area,,11,20,9
asel,a,area,,27
da,all,symm
allsel
asel,s,area,,9,10,1
asel,a,area,,19
sfa,all,1,pres,p
allsel
lswrite,2
lssolve,1,2,1
finish
以上是代码。
分析结果SINT应力值达到1700MPa，超许用应力太多，但是接触应力只有200MPa。这倒符合实际。现在就是不知道此模型的接触是不是有参数设置错误了，还请各位大虾帮忙看看。