找回密码
 注册
Simdroid-非首页
楼主: hh199422

[流固耦合] 边坡降雨入渗分析_含fish

[复制链接]
发表于 2016-4-13 14:24:02 | 显示全部楼层 来自 湖北武汉
感谢版主的无私奉献!
回复 不支持

使用道具 举报

发表于 2016-4-14 23:58:43 | 显示全部楼层 来自 湖北武汉
Simdroid开发平台
本帖最后由 dongxx0221 于 2016-4-15 12:14 编辑
hh199422 发表于 2016-4-12 21:19
你说的问题请你见我的论文就清楚了。这里一两句话说不明白

谢谢蒋老师,已拜读大作。但是贴出的代码不完整,只有调节渗透系数和饱和度的fish。故请教以下问题。

①初始非饱和土的孔压,里面的ini_surf_val是什么,怎么得到的?
  1. ;if gpn_pp <0 then
  2.    command
  3.      fix pp inipp_surf_val range id nm_gp
  4.    endcommand
  5. ;endif
复制代码

这个初始孔压场,水面以上的负压区是怎么生成的?

“地面产流前,在降雨入渗过程中,当降雨强度小于边坡表层土体最大入渗率时,入渗量取降雨强度;当降雨强度大于边坡表层土体最大入渗率时,入渗量取最大入渗率”。
按水文地质学,在特定土层、地貌条件下,入渗量=降雨量*入渗系数,入渗系数一般在0.15-0.20左右。
按文中表述,降雨强度不足以形成表面流时,全部入渗;当强度足以形成表面流时,按最大入渗强度。最大入渗量是由饱和渗透系数决定,入渗量和渗透系数有何计算关系?才能去判断按哪种强度入渗。

③"按照表面interface,通过interface 单元信息,寻找地表表层每一个实体单元(zone)的外表面,并在该表面上施加流量边界,边界上的流量数值可以随计算时间实施调整,以模拟不同时期降雨强度的变化;在每一个计算时间步,对比降雨强度值与边坡岩土体的入渗率(本研究取边坡表层实体单元饱和渗透系数)之间的大小关系,取两者之间的小值作为边界上的流量输入值qs 。"
"通过interface 节点信息,定位地表单元的节点位置,并在每一个计算步中,判断地表单元节点的孔隙水压力ps是否大于0。地表节点孔隙水压力大于0,表明该节点达到饱和状态,地表将出现积水,因此,在后一计算时间步中,将该节点的流量边界修改为压力边界(可取压力为0)。"
1类边界用pp,2类边界应该是指用discharge,1类边界和2类边界如何根据降雨量自动识别并转换的?比如案例中用到的降雨强度3×10-6 m/s (130mm/12h,特大暴雨)和表土的饱和渗透系数,对比饱和渗透系数ks与降雨量qs?这两个好像没法直接对比吧,不是一个概念。

“降雨停止后,由于降雨入渗进入斜坡体内的水体在重力作用下,继续保持向边坡体内和下部渗流的趋势,在边坡斜坡表面,坡体内的水分既可向边
界外溢出,也可能在重力作用下继续下渗。当斜坡上点因坡体内水分溢出而达到正孔隙水压力时,边坡表层节点的边界调整为压力为0 的压力边界。降
雨停止后边坡表面上的孔隙水压力在重力作用下变为负值时,渗流边界改为零流量边界,即坡体内的水分不再溢出。”

“降雨停止后,移除地表单元表面的流量边界。对地表节点的孔隙水压力进行实时监控,当地表节点孔隙水压力大于0 时,将该节点孔隙水压力
强制置零。节点孔隙水压力为0的部位即为坡体内的地下水溢出点。”

从降雨后1d、3d、5d图,发现-50kpa区没变化,上面入渗的水(50kpa区)只沿pp=50和pp=-50分割线向左下流,没有向下入渗。按常理向-50kpa负压区渗流比侧向渗流更容易。

⑤楼主这大作只是纯渗流没有流固耦合吗?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 6 不支持 0

使用道具 举报

发表于 2016-5-31 15:55:41 | 显示全部楼层 来自 湖南湘潭
正在学习当中,谢谢楼主的分享
回复 不支持

使用道具 举报

发表于 2016-8-4 10:02:09 | 显示全部楼层 来自 北京
还是有些没看懂啊 楼主还在么
回复 不支持

使用道具 举报

 楼主| 发表于 2016-8-7 20:45:08 | 显示全部楼层 来自 湖南长沙
现在只是偶尔上论坛,可以Email联系,相互探讨
回复 不支持

使用道具 举报

发表于 2016-9-11 10:38:06 | 显示全部楼层 来自 北京
研究了大约3个月了 还是没能做出来 哎~~~~
回复 不支持

使用道具 举报

发表于 2016-9-21 17:05:56 | 显示全部楼层 来自 广东广州
谢谢楼主分享啊
回复 不支持

使用道具 举报

发表于 2017-1-19 19:13:34 | 显示全部楼层 来自 湖北武汉
渗流真真是很高端啊,需要前人指路。
回复 不支持

使用道具 举报

发表于 2017-2-28 20:58:18 | 显示全部楼层 来自 甘肃兰州
请问楼主给出的是全部的fish语言吗?要是全部的话就太好了
回复 不支持

使用道具 举报

发表于 2017-6-3 17:35:48 | 显示全部楼层 来自 湖北武汉
原文中提到式(8)计算相对渗透系数,而全文公式最大编号只到式(7)。。。。并没有式8.。。

点评

VG模型是式6  发表于 2017-9-9 22:02
回复 不支持

使用道具 举报

 楼主| 发表于 2017-6-11 15:53:58 | 显示全部楼层 来自 湖南长沙
这是个蚕 发表于 2017-6-3 17:35
原文中提到式(8)计算相对渗透系数,而全文公式最大编号只到式(7)。。。。并没有式8.。。 ...

kr=s^2(3-2s)
回复 1 不支持 0

使用道具 举报

发表于 2017-9-9 11:00:57 | 显示全部楼层 来自 湖北武汉
本帖最后由 3Dgeosciences 于 2017-9-28 16:02 编辑

谢谢蒋教授的奉献~
回复 不支持

使用道具 举报

发表于 2018-6-21 15:23:34 | 显示全部楼层 来自 奥地利
老师,您好,您的初始孔压的分布,是怎么计算得到的
回复 不支持

使用道具 举报

发表于 2018-7-3 17:46:52 | 显示全部楼层 来自 陕西西安
感谢蒋老师无私的精神
回复 不支持

使用道具 举报

发表于 2018-10-19 21:50:34 | 显示全部楼层 来自 陕西西安
你好,我最近刚开始学习边坡入渗,对模拟降雨的过程还是不太熟悉,因为是1级会员,也不能发消息,我的邮箱是644479616@qq.com,可以留下QQ或者其他联系方式吗?我有些问题想请教,万分感谢!
回复 不支持

使用道具 举报

发表于 2018-11-22 15:15:37 | 显示全部楼层 来自 贵州贵阳
楼主能一起学习学习吗
回复 不支持

使用道具 举报

发表于 2019-1-24 15:33:43 来自手机 | 显示全部楼层 来自 四川凉山州西昌
请问降雨边界如何施加,请大神讲解一下
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2026-1-3 12:06 , Processed in 0.039693 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表