论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 6617阅读
  • 15回复

求助:希望问一下有关flac设置水位变动时边界条件的问题 [复制链接]

上一主题 下一主题
离线hbwzh
 

发帖
47
土币
60
威望
303
原创币
0
只看楼主 正序阅读 使用道具 楼主  发表于: 2010-12-22
在做的事有关水位变动对边坡稳定的影响,所以要用非稳定流来解决问题。
用flac3d在做的时候把水位变动的过程分为三个步:第一步,水位上升,第二步,高水位饱和,第三步,水位下降。在水位上升时通过不断释放更新固定边界孔压值来实现水位的上升过程,并在这个过程中更新边界静水压力;高水位饱和就是固定孔压,设定求解时间,没什么好说的;在水位下降时,采取了和水位上升时一样的方法,就是不断随时间更新孔压值。
在模拟水位上升时,做出来的孔压云图效果还行,至少像那么回事儿(由于模型尺寸偏小,时步很短,计算时间很长,不得已调整了流体模量,做出来的时间值和真实时间值有偏差),和seep计算结果的云图很像。在水位饱和过程中也可以,看得出来孔压随时间变化。
但是,水位下降就不是那么回事儿了:水位下降过程中,显示孔压只是在边界有改变,并没有反应到边坡内部。也就是说,只在边坡边界孔压有变化,而边坡内部,孔压还是很高,几乎不变。
我分析是flac3d默认的是不透水边界,当边界孔压上升时,与坡体内部有孔压差,而边坡是透水的,因而可以渗近来,反应到坡体内部的孔压值;但在水位下降时,虽然边坡边界的孔压值改变了,边坡内比边界及边界以外空压高,但是由于默认的是不透水的,因而水无法渗出,所以坡体内部的孔压无法改变。
不知道我的思路对不对,大家在做的时候遇到过这种问题没?希望求教
离线datong2012

发帖
385
土币
5406
威望
11
原创币
0
只看该作者 15楼 发表于: 2013-01-09
学习一下怎么玩渗流
离线hbwzh

发帖
47
土币
60
威望
303
原创币
0
只看该作者 14楼 发表于: 2010-12-27
回 13楼(haiqyantu) 的帖子
你看一下,只在边界的孔压改变了,坡内的还是没有改变,计算时间也很长。
def waterdown
fltime1=fltime-259200
pp_2=7.6e3-fltime1*7.6e3/10800
pp_2_2=pp_2-2e3
p3=-1*pp_2
p4=-1*pp_2_2
line=pp_2/1e4
if line>0.2
command
apply remove nstress
free pp ran z 0 x 3.46 4.02
apply pp pp_1 ran z 0 x 3.46 4.02
apply nstress p1 ran z 0 x 3.46 4.02
free pp ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 0.2
apply pp pp_2 grad 0 0 -1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 0.2
free pp ran x 3.1 3.16 z 0.1999 0.2001
apply pp pp_2_2  ran x 3.1 3.16 z 0.1999 0.2001
free pp ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
apply pp pp_2 grad 0 0 -1e4  ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 line
apply nstress p3 grad 0 0 1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 0.2
apply nstress p4 ran x 3.1 3.16 z 0.1999 0.2001
apply nstress p3 grad 0 0 1e4 ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 line
endcommand
else
command
apply remove nstress
free pp ran x 3.1 3.16 z 0.1999 0.2001
apply pp 0  ran x 3.1 3.16 z 0.1999 0.2001
fix pp ran x 3.1 3.16 z 0.1999 0.2001
free pp ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
apply pp 0 ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
fix pp  ran plane norm 1 0 2 ori 3.1 0 0.2 z 0.2 0.76
free pp ran plane norm 1 0 1.5 ori 0 0 3.46 z line 0.2
apply pp 0  ran plane norm 1 0 1.5 ori 0 0 3.46 z  line 0.2
fix pp  ran plane norm 1 0 1.5 ori 0 0 3.46 z  line 0.2
free pp ran plane norm 1 0 1.5 ori 0 0 3.46 z line 0.2
apply pp pp_2 grad 0 0 -1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
fix pp ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
apply nstress p3 grad 0 0 1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
free pp ran z 0 x 3.46 4.02
apply pp pp_1 ran z 0 x 3.46 4.02
apply nstress p1 ran z 0 x 3.46 4.02
endcommand
endif
end

set fishcall 10 waterdown
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 13楼 发表于: 2010-12-26
是透水边界。

apply pp 0 range
就可以设置透水边界
半生闲吟今终止,
一步江湖无尽期。
离线hbwzh

发帖
47
土币
60
威望
303
原创币
0
只看该作者 12楼 发表于: 2010-12-26
回 11楼(haiqyantu) 的帖子
当fix pp 0的时候还是透水边界么?
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 11楼 发表于: 2010-12-25
陈博士的书中p233、234给出一个真空预压的例子,其中对竖向排水通道的模拟就用fix pp设置的排水边界,
你上面说的“更新孔压值”,设定了孔压值,应该就是排水边界吧
半生闲吟今终止,
一步江湖无尽期。
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 10楼 发表于: 2010-12-25
陈育民博士的书里p216里提到了,手册中肯定也有了,
设置透水边界可采用FIX pp或APPLY pp命令,
那么反过来,用这两个命令设置孔压时,是否就可认为是透水边界了?

设置sat得用ini或prop吧

我没怎么算过渗流,只是刚好看到,多讨论啊
半生闲吟今终止,
一步江湖无尽期。
离线hbwzh

发帖
47
土币
60
威望
303
原创币
0
只看该作者 9楼 发表于: 2010-12-25
回 7楼(haiqyantu) 的帖子
我看到某论文里面说可以固定边界饱和度,我想如果把边界饱和度固定设置为0,应该也可以持续排水,但是这个固定用fix sat 系统提示有问题,fix后似乎不能跟sat,我用set sat如何?可以用么?
离线hbwzh

发帖
47
土币
60
威望
303
原创币
0
只看该作者 8楼 发表于: 2010-12-25
回 7楼(haiqyantu) 的帖子
用fix pp或apply pp设置的时候,就是透水边界么,哪本书里说了?
离线haiqyantu

发帖
186
土币
394
威望
226
原创币
0
只看该作者 7楼 发表于: 2010-12-25
向楼主学习一下啊。
不是说用fix pp或apply pp设置的时候,就是透水边界么?
在某一特别小的时段内,孔压能否认为不变呢,此时到底应认为透水还是不透水呢?
半生闲吟今终止,
一步江湖无尽期。
离线love36524

发帖
3160
土币
4687
威望
9899
原创币
0
只看该作者 6楼 发表于: 2010-12-23
回 5楼(hbwzh) 的帖子
昨天就在的啊,你没看到吗

不为失败找借口,只为成功找方法。
离线hbwzh

发帖
47
土币
60
威望
303
原创币
0
只看该作者 5楼 发表于: 2010-12-23
回 1楼(love36524) 的帖子
能不能放到紧急求助版块里面啊?我之前发了一个类似的,但是一直没人回复
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      浙公网安备 33010602003799号 浙ICP备14021682号-1

      工业和信息化部备案管理系统网站