一个双联拱隧道的例子,渗流的一塌糊涂
newconfig 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调用
Re:一个双联拱隧道的例子,渗流的一塌糊涂
见图片,这个程序仅提供参考,希望大家讨论哦Re:一个双联拱隧道的例子,渗流的一塌糊涂
cjiao兄,你的隧道的围岩是什么岩石?怎么渗透系数那么小?一般的量级是e-11cm/s.下面是我做的流固耦合,请大侠指点,我们多多交流!
Re:一个双联拱隧道的例子,渗流的一塌糊涂
矢量图Re:一个双联拱隧道的例子,渗流的一塌糊涂
难得糊涂~Re:一个双联拱隧道的例子,渗流的一塌糊涂
cjiao兄,你的隧道的围岩是什么岩石?怎么渗透系数那么小?一般的量级是e-11cm/s.下面是我做的流固耦合,请大侠指点,我们多多交流!
小还漏的一塌糊涂呢,大的话,就漏的没有水了哦,
这个不是渗透系数,好像是渗透性,大家都这么称呼,不知道有没有道理哦,
主要是二次衬砌和围岩之间,有的是放水布、有的是放水钢板,估计比较小小。
我的外面的围岩 确实 稍微小了点哦,只是一个 实例吗?!
Re:一个双联拱隧道的例子,渗流的一塌糊涂
哦!我明白了。cjiao兄在数值计算和编程方面很有造诣,希望多多交流!Re:一个双联拱隧道的例子,渗流的一塌糊涂
cjiao2000 wrote:这个不是渗透系数,好像是渗透性,大家都这么称呼,不知道有没有道理哦,
量纲为[长度/时间]的这个东东就叫渗透系数,也叫水力传导系数!它不仅取决于多孔介质的性质,而且和渗透液体的物理性质有关!
在渗流力学中还有一个与渗透性有关的东东,它的量纲是[长度的平方],称之为渗透率,它只与多孔介质的骨架性质有关.
Re:一个双联拱隧道的例子,渗流的一塌糊涂
在flac中,【渗透性】这个数值好像和传统得渗透系数有一个10e-6的关系,他好像是渗透系数除水的密度,这个见本版其他的讨论,量纲 也所有变化。Re:一个双联拱隧道的例子,渗流的一塌糊涂
今天一天没有来,因为停电,晚上来电,赶紧报道一下哦。Re:一个双联拱隧道的例子,渗流的一塌糊涂
cj2000兄:做得很好啊,以后关于水的问题要向你多多请教 了。
Re:一个双联拱隧道的例子,渗流的一塌糊涂
谦虚哦,这个源代码直接就可以用,稍微修改就可以哦。其实我是看了大家的帖子才有这样子的感受的,希望多多向你学习哦。Re:一个双联拱隧道的例子,渗流的一塌糊涂
只要放的源代码多了哦,那么大家都会了哦,我建议,所有的源代码,都奖励3~5分,不知道这个提议,总斑竹是不是支持啊,看看ansys版,看看人家的源代码,那才叫做爽。Re:一个双联拱隧道的例子,渗流的一塌糊涂
渗流着一块成立一个小集体,真是好,不知cjiao20004愿不愿意带我们学好渗流这一块?Re:一个双联拱隧道的例子,渗流的一塌糊涂
cjiao2000不知你想过没有这样一个问题,只有一层含水层,上部和下部都是隔水层,
怎么模拟这种现象,我模拟的结果很不好
Re:一个双联拱隧道的例子,渗流的一塌糊涂
我刚学习了渗流的手册,感觉几个物理概念还不是很清楚。手册上的咚咚,和我们教材上的不同啊。
Re:一个双联拱隧道的例子,渗流的一塌糊涂
是不一样,在flac中找不到弹性模量的身影,而且渗透系数也不一样。Re:一个双联拱隧道的例子,渗流的一塌糊涂
谢谢版主无私奉献!Re:一个双联拱隧道的例子,渗流的一塌糊涂
是不一样,在flac中找不到弹性模量的身影,而且渗透系数也不一样。if you select anisotropic model, you will meet Young's modulus.
Re:一个双联拱隧道的例子,渗流的一塌糊涂
我现在正在搞一个联拱隧道渗流实例,现出差在外,回来把成果整理汇总向各位高手学习一下!我在现场做了一个试验,结果发现裂隙渗透得一塌糊涂啊!爆破对渗流影响非常严重。
对了,flac能做裂隙渗流模拟吗?急。