楼主你好,解决不了你的问题,我倒是有个问题请教。
你的情况是 "由于接触脱开造成了不收敛。在梁端施加位移后,体系开始计算,但是当两个接触面发生滑移后,发生错误,认为接触对从接触成为了不接触",我想要的正是这种状态,因为我想看接触面参数降到多低时才出现滑移。但是我建了一个小模型发现无论参数多低都不会不收敛,就算摩擦系数降为零也一样。。不知道楼主是否清楚。下面是我的小模型,就是为了验证这个做的很小的模型,麻烦过目一下,谢谢!
- /prep7
- block,0,100,0,100,0,100
- k,,0,50,100
- k,,100,50,100
- k,,100,100,50
- k,,0,100,50
- a,9,10,11,12
- vsba,1,7
- vdele,2
- a,9,10,11,12
- va,2,8,9,10,11
- et,1,solid45
- et,2,170
- et,3,173
- mp,ex,1,2.6e10
- mp,prxy,1,0.167
- mp,dens,1,2000
- mp,mu,3,1.2 !摩擦系数
- R,1,,,, ,,
- RMORE,,,, ,,
- RMORE,2e2,,, ,, !粘聚力
- TYPE,1
- MAT,1
- REAL,1
- ESYS,0
- MSHAPE,1,3D
- MSHKEY,0
- CM,_Y,VOLU
- VSEL,s,,,1
- vsel,a,,,3
- CM,_Y1,VOLU
- CHKMSH,'VOLU'
- CMSEL,S,_Y
- VMESH,_Y1
- CMDELE,_Y
- CMDELE,_Y1
- CMDELE,_Y2
- asel,s,,,7
- cm,target,area
- type,2
- nsla,s,1
- esurf,all
- esel,s,type,,2
- eplot
- allsel
- asel,s,,,2
- aplot
- cm,contact,area
- type,3
- nsla,s,1
- esurf,all
- esel,s,type,,3
- eplot
- allsel
- NSEL,S,,,140
- NSEL,A,,,141
- NSEL,A,,,103
- NSEL,A,,,93
- D,all, , , , , ,UX,UY,UZ, , ,
- acel,,9.8,
- finish
- /solu
- NSUBST,5,20,1
- nlgeom,1
- allsel
- solve
复制代码 |