新手求教flac 3D流固耦合的一些问题

qixing0211

;问题描述：本问题为模拟土体上部加载，土体为软土，在荷载作用下，发生大变形，土体侧向限制变形，只允许竖向变形，土体模型为摩尔库伦模型。拟采用流固耦合模型解决问题，求得土体中各点的孔压、应力及变形量。
new
title
hydrulic filling model
config fluid
; --- mode geometry ---
gen zone brick size 20 1 10
; --- mechnical mode ---
model mohr
set fluid biot on
prop bulk=5e8, shear=1e8, friction=35, cohesion=2e10 &

tension=1e10
fix x y range x -.1 .1 x 19.9 20.1
fix x y z range z -.1 .1
fix y
ini den 2000 ;土颗粒密度2000，手册中很多例题都没有考虑土体密度，是不是因为忽略了重力的固结因素？
apply szz 2.e6 range z 19.9 20.1 ; 土体上部加载
; --- fluid flow model ---
model fl_iso
prop perm 2.3e-6 biot_c 1. ;渗透系数应该是指水在土体中的渗透系数吧？为什么手册例题上的系数取值都那么小呢？
ini biot_mod 4.e9 ;如果没有这一行命令，计算很快，但是没有结果，有了这一行之后，计算又不收敛，这一行命令到底是做什么的呢？还有根据手册上说，比奥模量值不能太大，我改小了之后还是一样的不收敛。
ini pp 0
ini den 1000 ;水密度1000
; --- setting ---
set gravity 0 0 -9.8
set large ;大变形计算
; --- fish functions ---
def point
pnt = gp_near(10.,0.,5.)
zpnt = gp_znear(10.,0.,5.)
end
point
def pp10
pp10 = gp_pp(pnt)
ft = fltime
c_szz = z_szz(zpnt)
end
; --- drained response ---
fix pp 0
range z 19.9 20.1
set fluid on
set mech force 0 ratio 1e-4
set hist_rep 200
set fluid substep 1
set mech substep 1 auto ;slave
solve age 500
; --- histories ---
hist pp10
hist c_szz
hist ft
; --- plot ---
plot create pp
plot set title text
Vertical displacement versus time at column mid-height
plot add his 1 vs 3
plot add his 2 vs 3
plot show
;运行不收敛！哪里出了问题呢？

lmm

这么小的字怎么看的清楚？？？？？？

qixing0211

字很小吗？我再发大点

qixing0211

;问题描述：本问题为模拟土体上部加载，土体为软土，在荷载作用下，发生大变形，土体侧向限制变形，只允许竖向变形，土体模型为摩尔库伦模型。拟采用流固耦合模型解决问题，求得土体中各点的孔压、应力及变形量。
new
title
hydrulic filling model
config fluid
; --- mode geometry ---
gen zone brick size 20 1 10
; --- mechnical mode ---
model mohr
set fluid biot on
prop bulk=5e8, shear=1e8, friction=35, cohesion=2e10 &

tension=1e10
fix x y range x -.1 .1 x 19.9 20.1
fix x y z range z -.1 .1
fix y
ini den 2000 ;土颗粒密度2000，手册中很多例题都没有考虑土体密度，是不是因为忽略了重力的固结因素？
apply szz 2.e6 range z 19.9 20.1 ; 土体上部加载
; --- fluid flow model ---
model fl_iso
prop perm 2.3e-6 biot_c 1. ;渗透系数应该是指水在土体中的渗透系数吧？为什么手册例题上的系数取值都那么小呢？
ini biot_mod 4.e9 ;如果没有这一行命令，计算很快，但是没有结果，有了这一行之后，计算又不收敛，这一行命令到底是做什么的呢？还有根据手册上说，比奥模量值不能太大，我改小了之后还是一样的不收敛。
ini pp 0
ini den 1000 ;水密度1000
; --- setting ---
set gravity 0 0 -9.8
set large ;大变形计算
; --- fish functions ---
def point
pnt = gp_near(10.,0.,5.)
zpnt = gp_znear(10.,0.,5.)
end
point
def pp10
pp10 = gp_pp(pnt)
ft = fltime
c_szz = z_szz(zpnt)
end
; --- drained response ---
fix pp 0
range z 19.9 20.1
set fluid on
set mech force 0 ratio 1e-4
set hist_rep 200
set fluid substep 1
set mech substep 1 auto ;slave
solve age 500
; --- histories ---
hist pp10
hist c_szz
hist ft
; --- plot ---
plot create pp
plot set title text
Vertical displacement versus time at column mid-height
plot add his 1 vs 3
plot add his 2 vs 3
plot show
;运行不收敛！哪里出了问题呢？

qixing0211

没变化啊，只能这么大了

jfchai

土颗粒密度2000，手册中很多例题都没有考虑土体密度，是不是因为忽略了重力的固结因素？* m; l3 p5
============
手册中是为例简化。

jfchai

手册中是为了简化。

yxf4171

建议你换成弹性模型试试，不要用大变形
我现在对流固耦合也不太懂，正在研究
呵呵，帮忙顶起啊……

s20070882

看来得大家一起努力，流固耦合比较难

dizidaiwang

土工实验得出的渗透系数k cm/s还要乘以10e-6才能在FLAC3D中用，
所以手册上的渗透系数小

chensha

建议回去再看MAUNAL。

chb168

我算的船闸闸首的固结沉降，沉降曲线最后趋于水平之后，又反升上去了，不知道为什么？郁闷

cumtsafe303

正在学习中，非常感谢

zlxiaoxi3

被流固耦合给整了的感觉

fsigh

流固耦合，我的痛

cxcqjj2003

是啊 流固耦合 整死人了  杯具

abaqus532

流固耦合杯具中的杯具

chenpf09

流固耦合确实很难

云起听雨

渗透系数要乘以10-6,

云起听雨

流固耦合确实很难，正在研究！