- 积分
- 0
- 注册时间
- 2006-8-10
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-12-10 10:19:42
|
显示全部楼层
来自 上海
哈哈,我这里有相应的3D的命令流,我加了备注,里面也有我的疑惑通通用?表示了,如果你也会一点儿3D的话我想你应该看能看懂,其实会2D也是相通的哈。。。也希望有大虾来指点。。。并且还有一个渗流速度的问题?我一直想把这个渗流速度怎样记录下来,可惜还是不知道怎么弄,希望大虾能给个意见!!!谢谢亲们咯
title
Steady state flow through a vertical embankment - case 1 ;一个竖向坝基上的稳态渗流之情况一
def setup
c_perm = 1e-10 ; 渗透系数1e-4cm/s
c_poro = 0.3 ; 孔隙率
c_kw = 1e3 ; 流体模量
c_L = 9. ; 模型长度
c_h1 = 6. ; 模型高度
c_h2 = 1.2 ; 含水层厚度
end
setup
config fl
gen zone brick size 30 1 20 p1 c_L 0 0 p2 0 0.15 0 p3 0 0 c_h1 ; 长10宽0.15高6的长方体几何模型
; --- fluid flow model ---
model fl_iso
prop perm c_perm poros c_poro ; default value: poros 0.5 biot_c 1 ;赋渗透系数和孔压值
set fl biot off
ini fmodulus c_kw
ini sat 0.0
ini sat 1 range z -.1 1.1999 ; 含水层饱和度赋1
ini sat 1 range x -.0001 .0001 ; 左边界饱和度赋1???不解何意(应该是外边界的水头高,应保持饱和)
ini pp 1.2e4 grad 0 0 -1e4 range z -.1 1.1999 ; 含水层孔压值由上至下线性增加
fix pp range x 8.99 9.01 ; 固定右侧孔压值
ini pp 6e4 grad 0 0 -1e4 range x -.0001 .0001 ; 左侧孔压由上至下线性增加
fix pp range x -.0001 .0001 ; 并将左侧孔压固定
; --- settings ---
ini fdensity 1e3
ini ftens 0.0
set gravity 0 0 -10
set mech off ; 力学计算关闭
set fl on ; 渗流计算开启
; --- test ---
solve ratio 1.e-3 ; 以最大不平衡力比率为1e-3控制计算
save ch2a.sav
def qflac ; 定义流量记录
qval = 0.0 ; 置初值归0
pnt = gp_head ; 节点指针
loop while pnt # null ; 节点指针非零即循环
fval = gp_xpos(pnt) - 0.001 ; 指针节点的x坐标减0.001函数
if fval < 0.0 then ; 若函数小于0
qval =qval + gp_flow(pnt) ; 那么将节点不平衡流量计入总流量
end_if
pnt = gp_next(pnt) ; 指针跳转
end_loop
qflac = qval /0.15 ; scale for unit thickness ; 单位长度(0.15)的流量
qsol = 1e-10*1e3*10.*(6.*6.-1.2*1.2)/(2.*9.) ; 单位换算系数m/s???后面的高度换算代表什么不解何意——这个应该是理论值-但是我还是不知道为什么要这么计算?(应该是按达西定律来的,q=kiA,其中1e-10*1e3*10.是k的单位换算从FLAC单位变成了m/s,然后(6.*6.-1.2*1.2)/(2.*9.)其实是i与A的合并,i=h/L,h=6-1.2,L=9,由于长度方向为1所以A=a=6+1.2,然后那个2是一半模型 ),上一个是数值解的值
end
print qflac qsol ; 输出以便查看理论和数值解的差别
pau
def seepage_face
array h2i(3),h2f(3)
h2i(1) = c_L
h2i(2) = 0.
h2i(3) = 1.2
h2f(1) = c_L
h2f(2) = 0.
h2f(3) = 1.8
status = set_line_width(2)
status = draw_line(h2i,h2f) ; 一条线由首未点构成?但这个status是什么东东还不明白
end
plot create view_case1
plot set plane origin 0 1e-5 0 normal 0 1 0
plot add flow plane
plot add grid plane
plot add fish seepage_face ; 注意这个命令流?
plot show
|
|