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

alloyer007

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 size  3 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.0023  range 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 sxy  0.09 0 25
hist z szz  0.09 0 25
hist z syz  0.09 0 25
hist z sxz  0.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 sxy  0 0.09 25
hist z szz  0 0.09 25
hist z syz  0 0.09 25
hist z sxz  0 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 sxy  0.06364 0.06364 25
hist z szz  0.06364 0.06364 25
hist z syz  0.06364 0.06364 25
hist z sxz  0.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 sxy  0.09 0 25
hist z szz  0.09 0 25
hist z syz  0.09 0 25
hist z sxz  0.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 sxy  0 0.09 25
hist z szz  0 0.09 25
hist z syz  0 0.09 25
hist z sxz  0 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 sxy  0.06364 0.06364 25
hist z szz  0.06364 0.06364 25
hist z syz  0.06364 0.06364 25
hist z sxz  0.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

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

alloyer007

高手快来看看~

who_is

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

alloyer007

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

lakewater

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

alloyer007

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

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




3\应力加载是否正确?

alloyer007

顶一下
请高手看看

elegant_2005

顶!!!

alloyer007

自顶~

sy1884

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

heibing1984

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

baoyun

power模型不需要粘性系数