FLAC3D蠕变分析QQ群

lg14507

在论坛里逛了很久，发现很多关于蠕变的帖子都没有人回，本人也在用FLAC3D做一个隧道的蠕变，也遇到很多问题，希望懂的朋友和学习的朋友能一起进来交流下。QQ群号：373217143

mr2014

我们能够联系啊下吗

yangxin546

可以问我，我做蠕变本构及其二次开发的。

lg14507

谢谢朋友，能加你QQ吗？我QQ1079004766。

lg14507

yangxin546 发表于 2014-2-19 16:56
可以问我，我做蠕变本构及其二次开发的。

谢谢朋友，能加你QQ吗？我QQ1079004766。

yangxin546

就在论坛里面说吧

angeljunior

yangxin546 发表于 2014-2-19 21:26
就在论坛里面说吧

你做过蠕变损伤的二次开发吗

yangxin546

做过，就这个问题？

lg14507

yangxin546 发表于 2014-2-19 21:26
就在论坛里面说吧

我在做的是软岩隧道的蠕变，你看下我出的一个拱顶和一个拱腰的位移图。这是算了10天的蠕变，拱顶的是持续增大，后面趋于稳定，而拱腰的怎么后面又变小了，而且算到最后还变成0了。

yangxin546

可能有很多原因，
1 塑性区发生了变化，导致该位置发生了破坏
2 你计算的时候是先静力平衡，再算蠕变，还是一起计算，如果是一起计算，可能会出现这种情况
3 蠕变太大，导致网格出现畸形

lg14507

yangxin546 发表于 2014-2-21 16:12
可能有很多原因，
1 塑性区发生了变化，导致该位置发生了破坏
2 你计算的时候是先静力平衡，再算蠕变，还是 ...

我是先用弹塑性模型算自重平衡，然后开挖也是用弹塑性模型算的，是算到平衡了再进行支护用蠕变计算的。

lg14507

yangxin546 发表于 2014-2-21 16:12
可能有很多原因，
1 塑性区发生了变化，导致该位置发生了破坏
2 你计算的时候是先静力平衡，再算蠕变，还是 ...

应该蠕变也不算太大的，才算到1.几个CM就成这样了，拱顶都算那么大了也还是正常的。

yangxin546

还有网格大小，最好截张图
塑性区大小，截张图

lg14507

yangxin546 发表于 2014-2-22 23:51
还有网格大小，最好截张图
塑性区大小，截张图

我在论坛里看到有人说他是拱腰的位移趋势是争取的，拱顶的位移趋势确实变小的。

lg14507

yangxin546 发表于 2014-2-22 23:51
还有网格大小，最好截张图
塑性区大小，截张图

还是希望能加你QQ问下，这样问起来很麻烦，难说清楚。。

yangxin546

从塑性区上来看，估计是破坏了，，，你得把参数C，fai，调大点再算。

lg14507

yangxin546 发表于 2014-2-22 23:51
还有网格大小，最好截张图
塑性区大小，截张图

问下关于FLAC3D手册中一个BURGER模型的蠕变例子中，最后P=-1后面这三个参数都表示什么吗？

白鹿中的赢家

请问一下，大神，我这个命令问题出现在哪里啊，总是出现math errorver flow invalid  file:\cycle.cpp line :466或者math error:over flow invalid  file:\cycle.cpp line :1276的错误，谢谢啦
; --- prop ---
model cpower range group 1
model cpower range group 2
model cpower range group 3

;material properties
prop bulk 10.87e9 shear 5.91e9 range group 1
prop cohesion 1e6 friction 35 tension 1e6 range group 1
prop a_1=1.4e-23 n_1=4.37 range group 1
prop bulk 15.42e9 shear 7.12e9 range group 2
prop cohesion 1.25e6 friction 40 tension 1e6 range group 2
prop a_1=6e-24 n_1=3.82 range group 2
prop bulk 8.33e9 shear 3.85e9 range group 3
prop cohesion 0.5e6 friction 30 tension 0.5e6 range group 3
prop a_1=1.2e-23 n_1=3.91 range group 3

; --- settings ---
ini density 2600 range group 1
ini density 2200 range group 2
ini density 2300 range group 3
set gravity 0 -9.81 0
set large

; --- 初始条件 ---
ini syy -16.8e6  szz -16.8e6  sxx -16.8e6

; --- 边界条件 ---
fix x ran x -250
fix x ran x 250
fix y ran y 0
fix z ran z 250
fix z ran z 0
apply syy -17.1e6 range y 799.9 800.1
; --- 监控 ---
save 初始模型.sav


set hist_rep 100

hist id=1 unbal
hist id=2 gp ydisp 0 480 0
hist id=3 gp ydisp 0 320 0
hist id=4 gp ydisp -50 400 0
hist id=5 gp ydisp  50 400 0

; --- small time elastic response ---
;pause
save 初始模型1.sav
solve
save 初始应力-1.sav

plot history 1
plot history 2
plot history 3
plot history 4
plot history 5
plot con disp ou on
plot con syy ou on
pause
; --- small time elastic response ---
solve
pause

;--- 施加椭球洞穴的内压---
ini xv=0 yv=0 zv=0 xdisp 0 ydisp 0 zdisp 0
apply nstress -3.0e6 range x -50 50 y 320 480 z 0 50
;pause
solve  
save 初始内压-1.sav
;pause
plot con syy ou on
plot con ydisp ou on
plot block state


; --- creep response ---
config creep
set large
print creep
set creep on
set creep mindt 1.0e-4 maxdt 1
set creep dt auto on
solve age 1
save 一年蠕变.sav
pause
plot con syy ou on
plot con ydisp ou on
plot block state



set creep mindt 1e-4 maxdt 1
set creep dt auto on
solve age 4
save 五年蠕变.sav
pause
plot con syy ou on
plot con ydisp ou on
plot block state