渗流入门例子，从入门到精通
;帮助的 第ftd138.pdf;------------------------------参数部分
bulk modulus, K 390 MPa体积模量（土质比较硬）
shear modulus, G 280 MPa 剪切模量
soil dry density, ρd 1200 kg/m3土的密度
water density, ρw 1000 kg/m3水的密度
wall density, ρwal 1500 kg/m3档土墙的密度
permeability, k 10−12 m2/Pa-s渗透系数
porosity, n 0.3 孔隙率
fluid bulk modulus, Kf 2.0 GPa流体的体积模量
;-------------------------------------
;-----以下是命令流
;------------------------------------------------------
; Excavation in a saturated soil
;------------------------------------------------------
config fluid;设置流体算法
; --- geometrical model --- 建立模型
gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1
group soil ;设置土 组
group excavate range x 0 4 y 0 4 z 0 5 ;设置开挖范围组
group wal1 range x 4 5 y 0 5 z 0 7 ; 设置档土墙 1 组
group wal2 range x 0 4 y 4 5 z 0 7 ;设置档土墙 2 组
group wall range group wal1 any group wal2 any ;档土墙1 和2 合并为 一个档土墙组
; --- fluid flow model ---流体模型
model fl_iso ;各项同性流体
prop perm 1e-12 poro 0.3 ;设置渗透系数和？？比 or ？？率
ini fdensity 1e3 ; 初始化 水密度
ini fmod 2.0e9 ftens -1e-3 ; 初始化水的体积模量 及 ？？
model fl_null range group excavate; 开挖部分没有水流
model fl_null range group wall ;档土墙 不渗水
ini pp 0 grad 0 0 1e4 ;初始化孔隙水压力
fix pp range z -0.1 0.1 ;约束 应该是整个模型的底部吧？
fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;开挖部分也给约束住了，而其他的地方不管
; --- mechanical model ---力学模型
model elas ; 弹塑性
prop bul 3.9e6 shea 2.8e6;设置体积模量和 剪切模量，这两个必须，有一次我只设置了体积模量，就错的一塌糊涂
model null range group excavate ;力学模型的开挖
ini density 1.2e3 ;初始化 密度 ，应该是土的密度
ini density 1.5e3 range group wall ;初始化墙（档土墙）的密度，
fix x range x -.1 .1 ;施加约束，据说差分法约束0 的时候，必须取一个-0.1<0<0.1,范围
fix x range x 11.9 12.1 ;施加约束，或者叫做边界条件
fix y range y -.1 .1 ;施加约束
fix y range y 11.9 12.1
fix z range z 11.9 12.1
; initial total stresses出世总应力
ini szz 0 grad 0 0 -1.5e4 ;z方向的应力大一些
ini sxx 0 grad 0 0 -1.2e4 ;水平方向的，小一些
ini syy 0 grad 0 0 -1.2e4 ;
apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0
apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0
apply nstress -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1
; --- setting ---设置
set gravity 0 0 10 ;设置重力加速度
; --- initial state --- ;初始化状态
solve force 1 ; check initial equilibrium求解平衡
; --- histories --- 记录
set hist_rep 40;记录步为40
hist fltime ;记录时间
hist gp pp 0 0 6;记录某点的孔隙水压力
hist gp xdis 4 0 0 ;水平位移
hist gp xdis 4 0 2 ;记录水平位移
hist gp xdis 4 2 0
hist gp zdis 0 0 5 ;记录竖向位移
hist gp zdis 2 0 5
hist gp zdis 4 0 5
hist gp zdis 2 2 5
hist gp zdis 4 2 5
hist gp zdis 4 4 5
hist gp zdis 10 0 1
hist gp zdis 10 0 2
;
; --- excavation ---开挖
set fluid off;关闭水
; apply pore pressure at walls在墙上施加孔隙水压力
apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0
apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0
apply nstress -5.e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1
solve ;force 1 ;求解
save exc1.sav ; 保存
;
; --- drainage --- 排水
apply remove nstress ;撤掉刚才的压力
def relaxsetup;定义一个 释放 函数的参数
step0 = step
end
relaxsetup ;调用
def relax ;定义一个释放函数
rstep = step - step0
if rstep < ncyc then
relax=1.0-(float(rstep)/float(ncyc))
else
relax = 0.0
endif
end ;结束
set ncyc = 1000 ;先赋值，随后调用
apply nstress 0 grad 0 0 -1.e4 hist relax &
range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0
apply nstress 0 grad 0 0 -1.e4 hist relax &
range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0
apply nstress -5.e4 hist relax &
range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1
cyc ncyc
solve ;求解
save exc2.sav ;保存
;
; --- percolation --- 渗透
set fluid on ;设置水算法
fix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;施加孔隙水压力
cyc 9000 ;循环
save exc3.sav ;保存
plot create excav ;建立一个显示的 视图 一下内容都是显示内容
plot set rot 200 0 195
plot set cent 6 6 6
plot set dist 39.18
plot set magn 0.8
plot set plane ori 0 0 0 normal 0 -1 0
plot add cont pp out on
plot add block group lgra range group wall
plot add cont pp int 10000 max 110000 range x 0 4 y 0 4 z 5 7
plot add flow plane
plot show
ret
Re:渗流入门例子，从入门到精通
特定区域的渗流量如何计算呢？做过的给一个思路吧，1. 计算6年之后的渗流状况：solve age 189216000 ; six years after excavation 6年之后渗流状况
2 请问各位大侠，在进行隧道涌水量的渗流分析时，水流速度是用gp_flow( )，还是用z_qx( ),z _qy( ), z_qz( )量来计算呢，
有人这样回复：I consider we should use z_qx( ),z _qy( ), z_qz( ) , because they are specific discharge vector.
3 有关渗透率的讨论
搞清楚了，渗透性的定义是渗透系数（m/sec）除以水的容重（Pa/m），所以其单位就成了：m2/Pa*sec
比如换算成day，简单做个乘法就可以了么？
比如1m2/pa-sec = 1*24*3600 m2/pa-day？
好像是这样：flac3D的“本质渗透系数”单位为m^2/(Pa/sec)，按计算公式：
K＝Kh(cm/sec)/Y(Pa/m)=Kh*10^-2(m/sec)/9.81*10^3(Pa/m)=Kh*10.2*10^-6
4 以上都是别人讨论的，能在这个实例上进行讨论吗，比如 计算1年之后，该漏水截面的渗水量？？？
Re:渗流入门例子，从入门到精通
模型图片，见附件Re:渗流入门例子，从入门到精通
这种教学方式好, 新手可以逐句学习.Re:渗流入门例子，从入门到精通
版主苦心一片，新手甘之如饴！好帖，得益匪浅！
Re:渗流入门例子，从入门到精通
to:cjiao2000你写的真好，对于我们学渗流的肯定帮助很大，支持你提出的观点，
愿给你一块学习
Re:渗流入门例子，从入门到精通
上面这个模型图应该倒过来才对,哪有基坑开挖在下面的啊?Re:渗流入门例子，从入门到精通
这学期正好开始学习渗流力学，谢谢版主:):):)Re:渗流入门例子，从入门到精通
这学期正好开始学习渗流力学，谢谢版主:):):)Re:渗流入门例子，从入门到精通
hao a 谢谢Re:渗流入门例子，从入门到精通
我算涌水量是用z_qx( ),z _qy( ), z_qz( )量来计算的，结果还算满意，但是对于不规则断面的涌水量计算还存在问题，有人说，手册里有FISH，可是我没找到，望知道的高手讨论！Re:渗流入门例子，从入门到精通
set gravity 0 0 10 ;设置重力加速度所以基坑在下面吧
Re:渗流入门例子，从入门到精通
读了 斑竹的 例题,我对于我的论文好象有一点开窍了flac 啊,不难学,
Re:渗流入门例子，从入门到精通
我运行了一下例题，有一点不懂：在set fluid off后的那个solve,我同步显示了flow vector，发现居然有变化，按说明书上的意思用了set fluid off后好像应该关闭了水流流动，那么，flow vector为什么还在变化呢？
请高手帮指点，谢谢！
Re:渗流入门例子，从入门到精通
apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0斑竹能详细解读这个命令各参数的含义么，能在曲面上加水压力么
Re:渗流入门例子，从入门到精通
这个太好了 谢谢！Re:渗流入门例子，从入门到精通
cjiao2000 斑竹：你好，我现在作基坑渗流方面的flac分析，希望和您能够联系，给一些指导，我的电子信箱是liilinn@126.com,liilinnqd@126.com，谢谢！