初始状态

zhaoge

model new
zone import 'shuangtaijiewangge.f3grid'
model title 'Excavation and Support for a Shallow Tunnel'

model configure fluid

fish define parameter
l=50
r0=11
r1=r0-0.001
r2=r0+0.001
h1=17
E=9e6
V=0.27
bulk=E/(3-6*V)
shear=E/(2+2*V)
end
@parameter
zone cmodel assign mohr-coulomb
zone property bulk @bulk shear @shear friction 63 cohesion 2e5 tension 1.0e5 dilation 0 range position-z 0 49
;
; --------------fluid flow model----------------
zone fluid cmodel assign isotropic
zone fluid property permeability 0.322
zone fluid property porosity 0.5
zone gridpoint initialize fluid-modulus 2e9
zone gridpoint initialize fluid-tension -1e15
zone initialize fluid-density 1000
;
; assign boundary conditions               ;閿?璐稿崵???鐭?t閿無璐??鑼呰矊??鍗璐窪T???a??鐭娓閿歽鑴
zone gridpoint fix velocity-x range position-x -40
zone gridpoint fix velocity-x range position-x 40
zone gridpoint fix velocity-z range position-z 0
zone gridpoint fix velocity-y range position-y 0
zone gridpoint fix velocity-y range position-y 50
;
; assign initial stress state             ;3?閿璐竱璋﹟鑴抽垾瑁
model gravity 0 0 -10
zone initialize density 2300
zone initialize stress-zz -1.23e7 gradient 0 0 23000
zone initialize stress-xx -4.55e6 gradient 0 0 8507
zone initialize stress-yy -4.55e6 gradient 0 0 8507
zone gridpoint initialize pore-pressure 4.37e6 gradient 0 0 -10000
;
model mechanical active on
model fluid active off
model history mechanical unbalanced-maximum
model solve unbalanced-maximum 0 or ratio 1e-5
model save 'chushi'
为什么出来的应力图是不太对

zhaoge

请问为何出来的孔隙水压力是负的