怎么载荷变大就脱离接触了？有命令流

haiyo

各位大侠帮忙给看看，下面是我计算成功的一个轴对称接触模型，但是把倒数第6行的载荷变大（比如为30或者200），就脱离接触了
不知该怎么调参数，请不吝赐教，在下谢谢啦！
finish
/clear
/prep7
k, 1,643.0906,689.5541         !建立几何模型，单位为mm
k, 2,644.7560,668.9149
k, 3,643.4967,671.9369
k, 4,643.0906,677.9331
k, 5,645.0506,635.6731
k, 6,643.4706,677.5531
k, 7,649.6456,689.5541
k, 8,649.6456,675.6731
k, 9,649.9256,675.6731
k,10,649.6900,635.6731
k,11,645.0377,635.3584
k,12,642.4506,635.2339
k,13,641.2456,635.2339
k,14,641.2456,636.8339
k,15,644.1009,636.8339
k,16,644.4735,634.1508
k,17,643.3856,676.9339
k,18,643.0556,671.8960
k,19,644.7456,668.7589
k,20,644.3871,668.7589
k,21,643.3856,671.9631
k,22,643.0556,676.9339
k,23,644.4656,633.3739
k,24,645.0456,633.3739
k,25,645.0456,632.2339
k,26,642.4506,632.2339
a,1,4,6,3,2,5,10,9,8,7
a,15,14,13,12,26,25,24,23,16,11,19,21,17,22,18,20
lsel,s,,,1,5
cm,guan_nei,line    !定义刚性接触面组件
lsel,s,,,20,22
cm,ke_wai,line      !定义柔性接触面组件
lsel,s,,,24,26
cm,ke_nei,line      !定义施加内压组件
et,1,plane42    !定义单元类型，轴对称
mp,ex,1,105E+3      !定义材料参数，单位为Mpa
mp,prxy,1,0.3
mp,ex,2,210E+3
mp,prxy,2,0.3
mp,mu,3,0.2
mat,2               !划分网格
esize,1
amesh,1
mat,1
esize,0.5
amesh,2
et,3,conta171       !定义接触单元类型
et,4,targe169
r,1,,,0.05
mat,3
type,3              !划分接触单元
lsel,s,,,ke_wai
nsll,s,1
esurf,all
allsel,all
type,4
lsel,s,,,guan_nei
nsll,s,1
esurf,all
allsel,all
keyopt,1,3,1
keyopt,3,5,1
/solu
lsel,s,,,10        
nsll,s,1
d,all,uy,0
allsel
lsel,s,,,15
nsll,s,1
d,all,uy,0
allsel
lsel,s,,,ke_nei   
sfl,all,pres,3! 若是将此处载荷变为30，或者200，就脱离接触了
allsel
solve
/post1
/DSCALE,1,1.0
plnsol,s,eqv

newton

你的约束好象不完全吧?  X方向没有约束,力大了当然跑了



压力为3的时候变形已经很大了,压力200的变形没法看了.

[ 本帖最后由 newton 于 2007-4-8 13:18 编辑 ]