空隙水压力不对，如何修改，附命令
zone import 'zhujiangquan.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 elastic
zone property bulk @bulk shear @shear
;
; --------------fluid flow model----------------
zone fluid cmodel assign isotropic
zone fluid property permeability 3.8e-10
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
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
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.49e6 gradient 0 0 -10000
model mechanical active on
model fluid active off
model solve unbalanced-maximum 0 or ratio 1e-3
为何空隙水压力是反的，而且还是负的
页:
[1]