博博博博 发表于 2010-11-6 10:05:59

蠕变分析过程

最近做了一个蠕变分析的模型，怎么总是出错，望高手指点，总是显示有非法的网格
new
config creep
gen zone cylinder p0 0 0 0 p1 0.025 0 0p2 0 0 0.1 p3 0 0.025 0 size5 15 8
gen zone reflect norm -1 0 0
gen zone reflect norm 0 -1 0
def setup
    bu = 2.8e9
    g_m =74.72e6
    g_k =56.03e6
    vis_m = 324.93e6
    vis_k = 342.03e6
    P = 0.2e6
    cons = 1./(3.*bu) + 1./g_m
    alpha = g_k/vis_k
end
setup
model cvisc
prop coh 1.2e6 tension 0.6e4
prop bulk=bu fric=20
prop mshear=g_m mvisc=vis_m
prop kshear=g_k kvisc=vis_k
fix z range z 0   
app szz=prange z 0.1
def strn_m
    strn_m = P*(crtime/vis_m + cons)/3.
    strn_k = P*((1.-exp(-alpha*crtime))/g_k + cons)/3.
    strn_b = P*(crtime/vis_m + (1.-exp(-alpha*crtime))/g_k + cons )/3.
end
strn_m
hist n=10
hist crtime
hist unbal
hist gp zdisp 0 0 10
hist strn_b
hist dt

博博博博 发表于 2010-11-12 22:45:40

怎么没有人回答啊

海底冬天 发表于 2010-11-21 12:49:16

海底冬天 发表于 2010-11-21 12:54:37

hnlgdxnyckwl 发表于 2010-11-26 00:54:50

:lol四楼很强大，学习了，顶一个。。。

baoyun 发表于 2010-11-27 14:25:39

3# 海底冬天 fish函数def strn_m
在这里是什么意思，可否解释一下呢？

lwjshyu 发表于 2011-4-10 18:47:05

给楼主留个qq409947501，请教两个问题

zhaoshuaicae 发表于 2011-8-7 09:42:23

def strn_m
; g2 F, _2 `5 z3 C- L    strn_m = P*(crtime/vis_m + cons)/34 y7 p! e$ R7 l" ^" I3 y
    strn_k = P*((1.-exp(-alpha*crtime))/g_k + cons)/3
! `, g, N& S/ k( |# G    strn_b = P*(crtime/vis_m + (1.-exp(-alpha*crtime))/g_k + cons )/3h5 G" ~+ U! l$ I& h4 R' _
end; P. P# W% O0 Y

dangwg 发表于 2011-8-11 16:42:56

楼主有没有做过流变的模拟啊？蠕变和流变应该没有本质的区别吧。

hjm460727715 发表于 2011-11-23 22:04:25

用一个不准确的等式来表示：流变=蠕变+松弛

hanmeigong2004 发表于 2012-5-26 22:13:46

海底冬天 发表于 2010-11-21 12:54 static/image/common/back.gif
经过对你的命令流稍作修改可以得到以下图形：[

得到个什么图形呢，为什么我call进去算了基本没得什么应力呢？请教!
页: [1]
查看完整版本: 蠕变分析过程