关于contact的奇怪现象

lossxie

这个本来是手册上的一条例题，我在整理contact的计算模式时偶然发现下述情况，当使用应力边界时结果就会莫名其妙，跪求斑竹和大侠解答
;create shear test model
new
round 0.001
bl (-0.05,-0.1) (-0.05,0.1) (0.25,0.1) (0.25,-0.1)
crack -1 0 1 0
crack 0 0.1 0 0
crack 0.2 0.1 0.2 0
del range -0.05 0 0 0.1
del range 0.2 0.25 0 0.1

gen quad 1 1 range 0 1 -1 0
gen quad 1 1 range 0 1 0 1

prop mat=1 d=2.6e-3 k=45000 g=30000      ;单位全为Mpa

;Coulomb joint model
change jcons=2
set jcondf=2
;set add_dil on
;
prop jmat=1 jkn=40000 jks=40000 &
            jfric=30 jdil=6 zdil=4e-4 jcoh=2.0 &
            jtension=10e10
;
; apply boundary conditions

; apply normal load
bound stress (0,0,6.98) range -1,1 0.09,0.11      
;                                                   ;如果使用的是位移边界则删除contact的时间主要看jtension和btol的关系如果
;                                                   ;jtension/jkn>btol则删除时间由jtension决定，当contact的内力>jtension时，contact
;                                                   ;被删除，反之由btol控制，当contact位移>btolcontact才会被删除。
;                                                   ;但是如果使用的是应力边界则结果就会莫名其妙，如本例只要stress的拉力大于6.99
;                                                   ;contact就会完全失效，无论jtension为多大,是我的d版问题吗？我用的是UDEC3.10
;bou yvel=0.1 range -1,1 0.09,0.11
bound xvel=0 range -0.06,-0.04 -1,1
bound xvel=0 range 0.24,0.26 -1,1
bound yvel=0 range -1,1 -0.11,-0.09
hist unbal nc=1
hist nstr 1.0e-3 0.0
hist ndis 1.0e-3 0.0
step 10
print contact force
print contact disp

lossxie

不知道，大家有没有试过，如果大家试过后出现了同一样的情况，不妨留言。
如果上述情况大家都有碰到的话，强烈建议斑竹，以加分征解形式把这个问题置顶，请求高手回答。
这个问题虽然看起来好像不大，但是实际上在很多地方我门对contact都是要jtension的，如果它无法正常工作，这样对我们的结果会有很大的影响。
希望斑竹大家努力回帖，多多支持

lossxie

怎么没有人响应，斑竹也一直不献身指点迷津，只好自己顶一下

jl19830927

我也是在做关于岩体节理模拟的。对接触问题比较关心，如果接触的抗拉强度都不起作用了，那结果就没有意义了。希望能听到高手的指点。

lossxie

是啊，不知道你是不是在计算上面命令时也出现了同样的问题？？大家一起顶

lossxie

我找到了4.0版得程序，计算同样发现上述问题，强烈建议斑竹置顶加分征解

最近这段忙毕业和工作，实在是比较忙，根本没时间上网。请见谅。

[ 本帖最后由 glj264 于 2007-3-13 22:22 编辑 ]

whohuxins

首先的问题是你要干什么？
你做这个例子的目的是什么？
因为UDEC是一个非常灵活的数值分析软件
一个数字的变化或者是一个命令的变化都会导致计算结果有很大的出入
烦请楼主告诉我们此算例的目的
我也好帮你解决

cundall

几个发现：
第一，不是正版和盗版的问题，使用了udec4.0的正版，还是出现上述问题
第二，如果在施加应力边界的同时，施加初始应力场
比如
insitu stress 0 0 7.00
bound stress (0,0,7.00) range -1,1 0.09,0.11
计算是没有问题的。

有可能的原因是7.0MPa的突加应力太大了，块体产生了很大的加速度，导致btol太大，从而接触删除。
但是无法解释为什么6.98可以，7.0就不可以，同样期待指点

cundall

恩，想了一下，给大家一个解释，看看是否正确
第一，udec中块体之间的接触采用cell的方式判断临近块体。当施加的初始应力过大时，块体在一个时步内才是了过大的加速度和位移，而块体的变形则在下一个时步才能得到反应，这就出现了原本相邻cell的块体不在相邻cell的误判，从而导致接触判断误认为这两个块体已经不相邻，从而检索不到接触关系。
比如，块体1在5号cell，块体2在4号cell，原本两者接触，由于块体2受到了过大的加速度，使得在一个时步内，其型心已经移动到3号cell内，从而触发接触检索，这时块体的变形还没有完成（下一步完成），于是出现了块体2“刚性”运动到3号cell的假想，由于其已经不再是相邻cell，所以块体接触检索检索不到这两者的接触。这是由于离散元本身时步步进的显式算法和接触关系算法共同造成的，一般情况下不会出现（没有那么大的初始加速度），大家可以放心使用。
－－－－－－－－－－－－
|  5  |    4   |  3    |
－－－－－－－－－－－－
|  6  |   7   |   8    |
－－－－－－－－－－－－

第二、关于6.99，这个数应该与块体的尺寸和密度有关，更改块体尺寸和密度，这个数应该会发生变化。我懒得试验了，大家可以试试。

第三、如果万不得已，碰到类似情况，建议大家施加初始应力场，这可以有效降低初始加速度，比如上述例子，只要施加初始应力场，不会出现这类问题，详细见本人楼上的帖子。

cundall

第二点解释不对，我再琢磨一下

liuxl04

这可能是一个bug

一夜

bound stress (0,0,6.98) range -1,1 0.09,0.11  
在计算中写成6.98e６，即可，一直在用

cumtsafe303

增加见识了。。。 :lol

yangliuren

我在编程序时也用到了contact的思想

liyongming1569

应该将6.98，改为6.98e6，因为这里的单位是Pa！

wuhai111

看样子学习3dec的人是在太少了啊

孤独

3DEC还没学呢

suihansanyou

谢谢您的帮助，呵呵

suihansanyou

谢谢您的帮助，呵呵

dead5845065

等待高手中。。。。