流固耦合-2

1、渗流时间步与渗流时长

①渗流时长Timestep={（zone尺寸^2+孔隙率）/（渗透系数k+比奥模量M）} * 步数

注：：自定义时间步的时间不长约≥1.25 倍渗流特征长度（时间步）时，软件会报错

Timestep too large for zone fluid calculation.

model new
Program threads 1
model title "2-1 timestep in fluid by chengxiaojie QQ845767883"
fish automatic-create off
model configure fluid
;model fluid timestep fix 1.9999999e-15
;model mechanical time-total 0.4e-2
zone create brick size 10 10 10 point 1 (10,0,0) point 2 (0,10,0) ...
                              point 3 (0,0,10)
; --- fluid flow model ---
zone fluid cmodel assign isotropic

zone fluid biot on

zone fluid property permeability 1 porosity 0.2 biot 0.6 ; 
zone gridpoint initialize saturation 0.9
zone gridpoint initialize biot 1e4 ;不是biot-modulus
;zone gridpoint initialize fluid-modulus 2e3 ;error in flac3d
zone water density 1000
zone face apply pore-pressure 1 range position-z 0
; --- settings ---


zone face skin
zone cmodel assign elastic  
zone property bulk 4e3 shear 2e3 density 3000
zone face apply velocity-x 0 range group 'West' or 'East'  ;or 'North' or 'South' or 'Bottom'
zone face apply velocity-y 0 range group 'North' or 'South'
zone face apply velocity-z 0 range group 'Bottom'
model mechanical active on
model fluid active on
zone fluid implicit off ;默认是显式积分法explicit；
;隐式分析法implicit计算模式只适用于介质保持完全饱和；
model gravity 10
model solve  fluid time-total 0.0019
model save '2-1'
fish define ts
     global fixtime=zone.fluid.timestep
     global shenliushijian=zone.fluid.time.total
     global bushu=global.step 
     global shijianbu=0
            shijianbu=shenliushijian/bushu
    global lixueshijianbu=mech.time.total/mech.step

end
@ts

list @fixtime @shenliushijian @bushu @shijianbu @lixueshijianbu

2、流固耦合的方式

①前后计算

model configure fluid
zone fluid active on
zone mechanical active off

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★then ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ 

model configure fluid
zone fluid active off
zone mechanical active on

②主从进程

model configure fluid
zone fluid active on
zone mechanical substep fm
zone mechanical slave on
model fluid substep fc

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★or ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ 

model configure fluid
zone fluid active on
zone fluid substep fc
zone fluid slave on
model mechanical substep fm

③同步求解

model configure fluid
zone fluid active on
zone mechanical active on

3、流固耦合的选择依据

help位置：flac理论背景→fluid-mechanical …→Solving Flow-Only …→Selection of …

4、渗流场参数关联力学参数的自定义耦合方法

问：渗流场状态的改变关联到力学参数上，最难解决的问题是什么？

答：渗流场很多流体相关的参数都是在节点上的，但是力学模型参数基本上都是在zone上，如何实现历遍后关联到zone上！

model new

model restore '2-4-1'
zone Property density 2850
fish define by_chengxiaojie

     t=0

    loop foreach local r gp.list
       t=gp.sat(r)+0.0001
       if t>0.5 then
         zzhizhen=gp.link.zone(r)
           v=gp.pos(r)
           wz=zone.near.all(v)
           zone.prop(zzhizhen,'friction')=25
           zone.prop(wz,'friction')=25
        endif
       if t >1 then
          zzhizhen=gp.link.zone(r)
           v=gp.pos(r)
           wz=zone.near.all(v)
           zone.prop(zzhizhen,'friction')=38
           zone.prop(wz,'friction')=38
        endif
        
    endloop

end
@by_chengxiaojie
model save '2-4-1sav'

5、真三维潜水面自定义耦合方法

1

1

6、流固耦合的控制和章节总结

1

1

1