一个双联拱隧道的例子，渗流的一塌糊涂

cjiao2000

new
config fluid

call Anstofl.dat
model fl_isotropic

prop perm 6.4e-20 poro 0.01

prop perm 6.4e-19 poro 0.015 range group 9  ;;
prop perm 6.4e-19 poro 0.01 range group 10
prop perm 6.4e-19 poro 0.01 range group 15
prop perm 6.4e-19 poro 0.015 range group 3    ;;
prop perm 6.4e-19 poro 0.01 range group 4
prop perm 6.4e-19 poro 0.01 range group 12

prop perm 6.4e-19 poro 0.025 range group 13  ;;
prop perm 6.4e-19 poro 0.02 range group 11
prop perm 6.4e-19 poro 0.02 range group 21
prop perm 6.4e-19 poro 0.025 range group 6    ;;
prop perm 6.4e-19 poro 0.02 range group 5
prop perm 6.4e-19 poro 0.02 range group 22

prop perm 6.4e-28 poro 0.0001 range group 18
prop perm 6.4e-28 poro 0.0001 range group 20
prop perm 6.4e-28 poro 0.0001 range group 23
;prop perm 6.4e-28 poro 0.0001 range group 27

;prop perm 6.4e-28 poro 0.0001 range group 126
prop perm 6.4e-28 poro 0.0001 range group 19
prop perm 6.4e-28 poro 0.0001 range group 16
prop perm 6.4e-28 poro 0.0001 range group 24

prop perm 6.4e-28 poro 0.0001 range x -11.51 11.51 z -45.13 -35.5

ini fdensity 1e3
ini fmod 2.0e9 ftens -1e-3
model fl_null range group 26
model fl_null range group 25

model fl_null range group 27
model fl_null range group 126

;model fl_null range group 23
;model fl_null range group 18
;model fl_null range group 20
;model fl_null range group 27
;model fl_null range group 126
;model fl_null range group 19
;model fl_null range group 16
;model fl_null range group 24
ini pp 0 grad 0 0 -1e4
fix pp 0 range z -0.1 0.1
;fix pp 0 range z -80.1 79.9
;fix pp range y -0.1 0.1
;fix pp range y 0.9 1.1
;fix pp range X -75.1 -74.9
;fix pp range X 74.9 75.1

fix pp 0 range group 25
fix pp 0 range group 26
;fix pp 0 range group block
plot cont pp
;----------------------------------------
mo mohr
pro bulk 15e9 she 6e9 fric 20 coh 200e3 ten 5e3 dil 3 range z -80 0
model null range group 25
model null range group 26

; assign boundary conditions
fix x range x -75.1 -74.9
fix x range x 74.9 75.1
fix z range z -80.1 -79.9
fix y range y -.1 .1
fix y range y 0.9 1.1

ini density 2650
ini density 2650 range group 5
ini density 2650 range group 6
ini density 2650 range group 22

ini density 2650 range group 21
ini density 2650 range group 11

ini szz 0 grad 0 0 1.5e4
ini sxx 0 grad 0 0 1.2e4
ini syy 0 grad 0 0 1.2e4

set grav 0 0 -10

set hist_rep 40
hist fltime

;set fluid off
;solve force 1000
;step 1000

step 1500
save cq003.sav

;该实例计算了一个双联拱隧道渗流场，有待进一步完善，此程序另存为text.txt，附件为模型文件解压缩到同一个目录，即可用flac调用

cjiao2000

见图片，这个程序仅提供参考，希望大家讨论哦

我为FLAC3D狂

cjiao兄，你的隧道的围岩是什么岩石？怎么渗透系数那么小？一般的量级是e-11cm/s.
下面是我做的流固耦合，请大侠指点，我们多多交流！

我为FLAC3D狂

矢量图

tmandysha

难得糊涂~

cjiao2000

cjiao兄，你的隧道的围岩是什么岩石？怎么渗透系数那么小？一般的量级是e-11cm/s.
下面是我做的流固耦合，请大侠指点，我们多多交流！

小还漏的一塌糊涂呢，大的话，就漏的没有水了哦，
这个不是渗透系数，好像是渗透性，大家都这么称呼，不知道有没有道理哦，

主要是二次衬砌和围岩之间，有的是放水布、有的是放水钢板，估计比较小小。
我的外面的围岩 确实 稍微小了点哦，只是一个 实例吗？！

我为FLAC3D狂

哦！我明白了。cjiao兄在数值计算和编程方面很有造诣，希望多多交流！

youxianyuanyang

cjiao2000 wrote:
这个不是渗透系数，好像是渗透性，大家都这么称呼，不知道有没有道理哦，

量纲为［长度／时间］的这个东东就叫渗透系数，也叫水力传导系数！它不仅取决于多孔介质的性质，而且和渗透液体的物理性质有关！
在渗流力学中还有一个与渗透性有关的东东，它的量纲是［长度的平方］，称之为渗透率，它只与多孔介质的骨架性质有关．

cjiao2000

在flac中，【渗透性】这个数值好像和传统得渗透系数有一个10e-6的关系，他好像是渗透系数除水的密度，这个见本版其他的讨论，量纲 也所有变化。

cjiao2000

今天一天没有来，因为停电，晚上来电，赶紧报道一下哦。

starsmoon

cj2000兄：
做得很好啊，以后关于水的问题要向你多多请教　了。

cjiao2000

谦虚哦，这个源代码直接就可以用，稍微修改就可以哦。其实我是看了大家的帖子才有这样子的感受的，希望多多向你学习哦。

cjiao2000

只要放的源代码多了哦，那么大家都会了哦，我建议，所有的源代码，都奖励3～5分，不知道这个提议，总斑竹是不是支持啊，看看ansys版，看看人家的源代码，那才叫做爽。

colorfish

渗流着一块成立一个小集体，真是好，不知cjiao20004愿不愿意带我们学好渗流这一块？

colorfish

cjiao2000
不知你想过没有这样一个问题，只有一层含水层，上部和下部都是隔水层，
怎么模拟这种现象，我模拟的结果很不好

starsmoon

我刚学习了渗流的手册，感觉几个物理概念还不是很清楚。
手册上的咚咚，和我们教材上的不同啊。

cjiao2000

是不一样，在flac中找不到弹性模量的身影,而且渗透系数也不一样。

ying111

谢谢版主无私奉献！

hh199422

是不一样，在flac中找不到弹性模量的身影,而且渗透系数也不一样。

if you select anisotropic model, you will meet Young's modulus.

lidiyuan123

我现在正在搞一个联拱隧道渗流实例，现出差在外，回来把成果整理汇总向各位高手学习一下！
我在现场做了一个试验，结果发现裂隙渗透得一塌糊涂啊！爆破对渗流影响非常严重。
对了，flac能做裂隙渗流模拟吗？急。