qixing0211 发表于 2009-1-10 14:16:24

求教一个简单的流固耦合问题

很简单的一个模型,模拟土体在自重应力下一维固结。
new
config fluid
; --- geometry mode ---
gen zone brick size 20 1 10
; --- mechnical mode ---
model mohr
prop bulk=5e8, shear=1e8, friction=35, cohesion=2e10 &
   tension=1e10
fix x y range x -.1 .1 x 19.9 20.1 ;模型左右两边x、y固定
fix x y z range z -.1 .1;模型底部固定,且不透水
fix y;平面问题
ini dens 1500;土体干密度,随意取的数值
; --- fluid flow model ---
model fl_iso
prop perm 2.3e-8 poro 0.3
ini fmod 2e9
ini fden 1000
ini pp 0
; --- setting ---
set gravity 0 0 -9.8
set large;软土,会发生大变形
; --- first establish undrained response --- 建立自重应力平衡
set fluid off
solve ratio 1e-4
save con_und.sav
; --- drained response ---开始流固耦合计算
fix pp 0range z 19.9 20.1;模型顶部为排水面
set fluid on
set mech force 0 ratio 1e-4
set fluid substep 1
set mech substep 10 auto ;slave substep该怎么取值?
step 1000
; --- histories ---
hist fltime
hist gp zdis 10 1 5
hist gp pp 10 1 6
hist zone szz 10 1 5
plot his 2 vs 1
plot show
;问题:图中无曲线出现,但是等值线图显示有计算结果,而且发现曲线图坐标轴的数值非常大,达10e20,附图。这是怎么回事呢?到底是显示的原因,还是命令的问题?求高手帮忙,谢谢!

qixing0211 发表于 2009-1-10 14:28:55

呵呵,我有些傻了,我应该把监测的命令放在step前面的

lookcity 发表于 2009-1-10 14:52:51

是的 :) 呵呵

disciple_ 发表于 2009-8-12 23:39:22

:lol
学习下!
页: [1]
查看完整版本: 求教一个简单的流固耦合问题