论坛 › G0201：FLAC3D/FLAC2D › 一个自编蠕变的程序,希望大家能给点意见,更希望斑竹和各位大虾莅临指导

alloyer007 发表于 2006-8-1 01:09:26

一个自编蠕变的程序,希望大家能给点意见,更希望斑竹和各位大虾莅临指导

config creep
gen zone cshell size 5 100 20 p0 0,0,0 p1 0,0.090,0 p2 0,0,50 p3 0.090,0,0 p4 0,0.090,50 p5 0.090,0,50 p8 0,0.079,0 p9 0.079,0,0 p10 0,0.079,50 p11 0.079,0,50 group taoguan

gen zone cshell size3 100 20 p0 0,0,0 p1 0,0.1095,0 p2 0,0,50 p3 0.1095,0,0 p4 0,0.1095,50 p5 0.1095,0,50 p8 0,0.090,0 p9 0.090,0,0 p10 0,0.090,50 p11 0.090,0,50 group shuinihuan

gen zone radcyl size 5 100 20 20 p0 0,0,0 p1 0,25,0 p2 0,0,50 p3 25,0,0 p4 0,25,50 p8 0,0.1095,0 p9 0.1095,0,0 p10 0,0.1095,50 p11 0.1095,0,50 group weiyan

model power range group weiyan
prop bulk=1.002e3 shear=0.301e3 a_1=3.9e-7 n_1=3 dens=0.0023 range model power
model elastic range group shuinihuan
prop bulk=1.051e4 shear=8.9974e3 dens=0.002 range group shuinihuan
model elastic range group taoguan
prop bulk=17.5e4 shear=8.0769e4 dens=0.0078 range group taoguan

ini szz=-110.4 grad 0,0,-0.02254 range group weiyan
ini sxx=-63.08 grad 0,0,-0.01287 range group weiyan
ini syy=-63.08 grad 0,0,-0.01287 range group weiyan
ini density 0.0023range group weiyan
apply szz=-110.4 range z 49.99 50.01

set gravity 0 0 -9.81
fix x range x -0.001 0.001
fix x range x 24.999 25.001
fix y range y -0.001 0.001
fix y range y 24.999 25.001
fix z range z -0.001 0.001
set creep off

hist unbal
hist z sxx 0.09 0 25
hist z syy 0.09 0 25
hist z sxy0.09 0 25
hist z szz0.09 0 25
hist z syz0.09 0 25
hist z sxz0.09 0 25
hist z smax 0.09 0 25
hist z smin 0.09 0 25
hist z smid 0.09 0 25

hist z sxx 0 0.09 25
hist z syy 0 0.09 25
hist z sxy0 0.09 25
hist z szz0 0.09 25
hist z syz0 0.09 25
hist z sxz0 0.09 25
hist z smax 0 0.09 25
hist z smin 0 0.09 25
hist z smid 0 0.09 25

hist z sxx 0.06364 0.06364 25
hist z syy 0.06364 0.06364 25
hist z sxy0.06364 0.06364 25
hist z szz0.06364 0.06364 25
hist z syz0.06364 0.06364 25
hist z sxz0.06364 0.06364 25
hist z smax 0.06364 0.06364 25
hist z smin 0.06364 0.06364 25
hist z smid 0.06364 0.06364 25

step 4000

set creep on
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
set cr dt 1.141553e-4
set cr mindt 1.0e-4
set cr maxdt 8.0e-4
set cr dt auto on

hist reset
hist unbal
hist crtime
hist dt   


hist z sxx 0.09 0 25
hist z syy 0.09 0 25
hist z sxy0.09 0 25
hist z szz0.09 0 25
hist z syz0.09 0 25
hist z sxz0.09 0 25
hist z smax 0.09 0 25
hist z smin 0.09 0 25
hist z smid 0.09 0 25

hist z sxx 0 0.09 25
hist z syy 0 0.09 25
hist z sxy0 0.09 25
hist z szz0 0.09 25
hist z syz0 0.09 25
hist z sxz0 0.09 25
hist z smax 0 0.09 25
hist z smin 0 0.09 25
hist z smid 0 0.09 25

hist z sxx 0.06364 0.06364 25
hist z syy 0.06364 0.06364 25
hist z sxy0.06364 0.06364 25
hist z szz0.06364 0.06364 25
hist z syz0.06364 0.06364 25
hist z sxz0.06364 0.06364 25
hist z smax 0.06364 0.06364 25
hist z smin 0.06364 0.06364 25
hist z smid 0.06364 0.06364 25

solve age 5
这是我写的一个程序,不知道对不对?望大家批评指正~
谢谢~

freddie_1999 发表于 2006-8-1 09:21:28

蠕变还没接触过，帮你顶一下

alloyer007 发表于 2006-8-1 13:05:59

高手快来看看~

who_is 发表于 2006-8-1 14:30:49

蠕变参数都没输入，怎么算？
粘性系数。

alloyer007 发表于 2006-8-1 21:36:02

power模型没有粘性系数这个参数吧?

lakewater 发表于 2006-8-3 19:43:16

给大家一些前后处理的图片看看啊？

alloyer007 发表于 2006-8-5 16:46:02

我得到的结果不是最大不平衡力不收敛,就是考察点处的力很不规则,与实际和理论分析相差很大.(按照理论分析及其他资料说明考察点处的力是随时间逐渐增大的,且最后趋近于一个最大值)
我分析了一下:
问题可能出现在以下几个方面:

1\模型及其属性定义是否正确.我建的模型分3层,各层分别定义属性,但不知道属性是否都正确的加进去了?
2\初始和边界条件.本例中,下表面受垂直该面的竖向简支约束,前后左右四表面均作用有发向约束.在上表面作用一个下压力,在两外表面分别施加应力载荷,考察最里层结构的外围边界受力情况.
如图:

http://cn.f11.yahoofs.com/users/44a3b77bz6e43b206/d28b/__sr_/3e34cnb.jpg?phAyF1EBbdy_9ssc


3\应力加载是否正确?

alloyer007 发表于 2006-8-6 14:52:44

顶一下
请高手看看

elegant_2005 发表于 2006-8-6 19:29:46

顶!!!

alloyer007 发表于 2006-8-7 21:16:51

自顶~

sy1884 发表于 2008-12-9 14:52:27

怎么没看见有高手指点啊！同样问题！

heibing1984 发表于 2010-1-17 12:33:16

订我也遇到同样的问题 不平衡力好像趋于一个很大的值

baoyun 发表于 2010-8-4 16:58:45

power模型不需要粘性系数

查看完整版本: 一个自编蠕变的程序,希望大家能给点意见,更希望斑竹和各位大虾莅临指导