流固耦合中为什么---主从进程法---只能运行从进程

jlinz2008

new
config fluid
;-------------------------------------------------------------------------------------------------
gen zone radcylinder p0 0 0 0 p1 2.5 0 0 p2 0 0 -0.7 p3 0 2.5 0 &
         dim 0.25 0.25 0.25 0.25 size 3 3 4 10 ratio 1 1 1 1.1
gen zone radcylinder p0 0 0 -.7 p1 2.5 0 -.7 p2 0 0 -2.5 p3 0 2.5 -.7 &
         dim 0.25 0.25 0.25 0.25 size 3 6 4 10 ratio 1 1 1 1.1
gen zone radcylinder p0 0 0 -2.5 p1 2.5 0 -2.5 p2 0 0 -3.6 p3 0 2.5 -2.5 &
        dim 0.25 0.25 0.25 0.25 size 3 3 4 10 ratio 1 1 1 1.1  
gen zone radcylinder p0 0 0 -3.6 p1 2.5 0 -3.6 p2 0 0 -6.2 p3 0 2.5 -3.6 &
        dim 0.25 0.25 0.25 0.25 size 3 6 4 10 ratio 1 1 1 1.1

plot add axes black
;---------------------------generate stonepile------------------------------------------
gen zone cyl p0 0 0 0    p1 0.25 0 0    p2 0 0 -.7  p3 0 .25 0     size 3 3 4 ratio 1 1 1
gen zone cyl p0 0 0 -.7  p1 0.25 0 -.7  p2 0 0 -2.5 p3 0 .25 -.7   size 3 6 4 ratio 1 1 1
gen zone cyl p0 0 0 -2.5 p1 0.25 0 -2.5 p2 0 0 -3.6 p3 0 .25 -2.5  size 3 3 4 ratio 1 1 1
gen zone cyl p0 0 0 -3.6 p1 0.25 0 -3.6 p2 0 0 -6.2 p3 0 .25 -3.6  size 3 6 4 ratio 1 1 1
group stone range cylinder end1 0 0 0 end 0 0 -6.2 radius 0.25
;-----------------------bengou----------------------------------------------------
ini den 1.56e3 range z 0 -0.7
ini den 1.85e3 range z -0.7 -2.5
ini den 1.25e3 range z -2.5 -3.6
ini den 1.89e3 range z -3.6 -6.2
set gravity 0 0 -10

model mohr
prop bulk 7.43e7 shear 3.43e7 tension 0 coh 4e3 fric 25 dilation 5 range z 0 -0.7
prop bulk 0.93e7 shear 0.43e7 tension 0 coh 4e3 fric 28 dilation 5 range z -0.7 -2.5
prop bulk 0.49e7 shear 0.23e7 tension 0 coh 0   fric 25 dilation 0 range z -2.5 -3.6
prop bulk 1.55e7 shear 0.72e7 tension 0 coh 2e3 fric 35 dilation 25 range z -3.6 -6.2
;
;
;---------------------------------------------------------------------------------
;moder mohr range group stone
;ini den 2.07e3 range group stone
;prop bulk 2.78e7 shear 1.28e7 coh 5e3 fric 45 tension 0 di 10 range group stone
;
;
;---------------------displacement -----------------------------------------------
fix x range x -0.1 0.1
fix x range x 2.4 2.6
fix y range y -0.1 0.1
fix y range y 2.4 2.6
fix z range z -6.1 -6.2
;---------------------stress ---------------------------------------------------------
;
ini sxx 0 grad 0 0 7.8e3 range z 0 -0.7
ini syy 0 grad 0 0 7.8e3 range z 0 -0.7
ini szz 0 grad 0 0 15.6e3 range z 0 -0.7
ini szz 2.03e3 grad 0 0 18.5e3 range z -0.7 -2.5
ini sxx 1.015e3 grad 0 0 9.25e3 range z -0.7 -2.5
ini syy 1.015e3 grad 0 0 9.25e3 range z -0.7 -2.5
ini szz -12.97e3 grad 0 0 12.5e3 range z -2.5 -3.6
ini sxx -6.485e3 grad 0 0 6.25e3 range z -2.5 -3.6
ini syy -6.485e3 grad 0 0 6.25e3 range z -2.5 -3.6
ini szz 10.07e3 grad 0 0 18.9e3 range z -3.6 -6.2
ini sxx 5.035e3 grad 0 0 9.45e3 range z -3.6 -6.2
ini syy 5.035e3 grad 0 0 9.45e3 range z -3.6 -6.2

;------------------------pore  pressure------------------------------------------
;
water density 1000
water table ori 0 0 -0.7 norm 0 0 1
;
;
;-----------------------model fluid---------------------------------------------

model fl_iso
ini fden 1000 ften -1e10
ini fmod 2e9
ini sat 1.0 range z -0.7 -6.2
ini sat 0.0 range z 0 -0.7
prop por 0.444 perm 1.2e-5 range z 0 -.7
prop por 0.5 perm 5.8e-6  range z -.7 -2.5
prop por 0.55 perm 1.7e-9 range z -2.5 -3.6
prop por 0.41 perm 5.0e-5 range z -3.6 -6.2
prop por 0.23 perm 1.2e-4 range group stone
;---------------------hydraulic confinement-------------------------------------
fix pp 0 range z -0.7
fix pp 55e3 range z -6.1 -6.3
apply pp -7e3 grad 0 0 -10e3 range x -0.1 0.1 z -0.7 -6.2
apply pp -7e3 grad 0 0 -10e3 range x 2.4 2.6 z -0.7 -6.2
apply pp -7e3 grad 0 0 -10e3 range y -0.1 0.1 z -0.7 -6.2
apply pp -7e3 grad 0 0 -10e3 range y 2.4 2.6 z -0.7 -6.2

;solve
;save ppbalance.sav
;ini xdis 0 ydis 0 zdis 0 xvel 0 yvel 0 zvel 0
;set mech on fluid off

;
;--------------------running form---------------------------------------------------
set fluid substep 50000
set mech substep 6000 auto


hist gp pp 0 0 -2.5
plot create histview
plot hist

plot show
solve age 20e6
save balance1.sav

望高手指点下--后面的主从模式有问题么？这个只能运行力学模式，为什么渗透模式运行不了？谢谢

mazhongyuan9049

渗流模型开启后 本构只能是弹性的（比奥模型） 谈塑性（mohr-coulomb）模型加上渗流模式就关闭了！

Andy2947

主从进程步数设置的不合理！

yanqisha01

mazhongyuan9049 发表于 2009-10-29 09:15
渗流模型开启后 本构只能是弹性的（比奥模型） 谈塑性（mohr-coulomb）模型加上渗流模式就关闭了！ ...

我用mohr模型还是可以计算的