模型是一个长方体,三向尺寸分别为(x,y,z)=(600,5, 110)
,均质饱和可液化砂,流体动力耦合分析,监测孔压变化。流固耦合动力分析-孔隙水压力变化,命令流如下,运行之后发现,如图,
不知道为什么速度时程曲线有漂移?0.1s
之前,输入时程为0
,速度和孔压时程0.1S
之前的那段波动不知道从何而来?请教一下命令流是否哪里有误?或提供给我一个成功运行的超静孔压变化实例。谢谢哈
谢谢哈
new
impgrid Basin45.Flac3D
config dyn fluid
; *************************************************************************************
; ********************************************Static************************************
set dyn off fluid off
; ***************************************Mechanical property
model mechanical elastic
property density 2000 shear 2e9 bulk 3e9
; ***************************************Fluid property
model fluid fl_iso
property poro 0.5 perm 1e-8
;ini fmod 2e9
ini fmod 0
ini fdensity 1e3
ini sat 1
; ***************************************Gravity
set grav 0 0 -10
; ***************************************Boundary
fix x range x -0.1 0.1
fix x range x 599.9 600.1
fix y range y -0.1 0.1
fix y range y 4.99 5.01
fix z range z -0.1 0.1
hist id=900 unbalance
plot hist 900
hist id=901 gp yvel 300,-5,110
plot hist 901
plot zcontour sxx
plot zcontour syy
plot zcontour szz
plot add boundary
plot add axes
plot contour gpppressure
plot velocity
plot displacement
; ***************************************Initial distribution must be consistent with the gravitational gradient
ini pp 1.1e6 gradient 0 0 -1e4 range z 0 110
ini szz -2.75e6 grad 0 0 2.5e4 range z 0 110
ini sxx -8.1e5 grad 0 0 7.4e3 range z 0 110
ini syy -8.1e5 grad 0 0 7.4e3 range z 0 110
set mech rat 1e-6
solve
pause
initial xvelocity 0
initial yvelocity 0
initial zvelocity 0
initial xdis 0
initial ydis 0
initial zdis 0
;save static.sav
; *************************************************************************************
; ********************************************Dynamic************************************
set dyn on fluid on
ini fmod 2e9
model mechanical finn
property density 2000 shear 2e9 bulk 3e9 friction 35
property ff_latency=50 ff_switch=0 ff_c1=0.8 ff_c2=0.795 ff_c3=0.45 ff_c4=0.73
; ***************************************Dynamic boudary and loading
fix pp 0 range z 109.9 110.1
;fix pp range x -0.1 0.1
;fix pp range x 599.9 601.1
fix x z
free x y z range z -0.1 0.1
apply dquiet squiet range z -0.1 0.1
table 1 read ricker-7.3.dat
apply szy -2e6 hist table 1 szx 0.0 szz 0.0 range z -0.1 0.1
apply ff
set dyn damp local 0.0628
hist add id=81 gp ppressure 300,-5,90
hist add id=82 gp ppressure 300,-5,60
hist add id=83 gp ppressure 300,-5,10
hist id=1000 dytime
set mech ratio 1.0e-5
set dyn multi on
set large
plot hist 83 vs 1000
plot hist 14 vs 1000
solve age 40.0
return