level 4
模拟边坡的流固耦合,不知道哪里出问题了,一直不收敛,求大佬看看
model config fluid
zone cmodel assign mohr-coulomb
zone property bulk 3e9 shear 3.1e8 friction 45 range group 'chuqiba'
zone initialize density 2250 range group 'chuqiba'
zone property bulk 3.43e7 shear 6.32e6 friction 16 cohesion 1.078e3 range group 'jiyan'
zone initialize density 1950 range group 'jiyan'
zone property bulk 3.148e7 shear 6.03e6 friction 28 cohesion 9.8e3 range group 'weikuang1'
zone initialize density 1900 range group 'weikuang1'
zone property bulk 2.667e7 shear 5.71e6 friction 30 cohesion 9.8e3 range group 'weikuang2'
zone initialize density 1850 range group 'weikuang2'
zone property bulk 2.273e7 shear 5.4e6 friction 33 cohesion 7.84e3 range group 'weikuang3'
zone initialize density 1800 range group 'weikuang3'
zone fluid cmodel assign isotropic
zone fluid property permeability 5e-6 poros 1.5 range group 'chuqiba'
zone fluid property permeability 3.06e-12 poros 1.0 range group 'jiyan'
zone fluid property permeability 1.275e-10 poros 0.95 range group 'weikuang1'
zone fluid property permeability 3.825e-10 poros 0.9 range group 'weikuang2'
zone fluid property permeability 1.428e-10 poros 0.85 range group 'weikuang3'
zone gridpoint initialize fluid-modulus 2.2e9
zone gridpoint initialize fluid-tension 0
zone initialize fluid-density 1000
zone gridpoint fix velocity-x range position-x 0
zone gridpoint fix velocity-x range position-x 200
zone gridpoint fix velocity range position-z 0
zone gridpoint fix velocity-y
model gravity 0,0,-10
model fluid active off
model history mechanical ratio-maximum
model solve
model save 'h20210506001'
zone gridpoint initialize displacement (0,0,0)
zone gridpoint initialize velocity (0,0,0)
model fluid active on
model fluid substep 100
model mechanical substep 100
model mechanical slave on
model history name 'fluid-alltime' fluid time-total
zone history name 'chuqiba-pp1' pore-pressure position 26.5, 0, 5.0
zone history name 'chuqiba-pp2' pore-pressure position 26.5, 30, 5.0
zone history name 'chuqiba-pp3' pore-pressure position 26.5, 60, 5.0
zone history name 'chuqiba-pp4' pore-pressure position 21.58, 0, 7.76
zone history name 'chuqiba-pp5' pore-pressure position 21.58, 30, 7.76
zone history name 'chuqiba-pp6' pore-pressure position 30, 0, 2.5
zone history name 'chuqiba-pp7' pore-pressure position 30, 0, 2.5
zone history name 'pp8' pore-pressure position 70, 0, 0
zone history name 'pp9' pore-pressure position 135, 0, 0
zone history name 'pp10' pore-pressure position 195, 0, 0
zone history name 'pp11' pore-pressure position 135, 0, 20
zone history name 'pp12' pore-pressure position 180, 0, 35
zone history name 'pp13' pore-pressure position 70, 0, 20
zone history name 'pp14' pore-pressure position 115, 0, 35
zone history name 'pp15' pore-pressure position 165, 0, 50
zone history name 'pp16' pore-pressure position 50, 0, 20
zone history name 'pp17' pore-pressure position 140, 0, 55
zone history name 'dis1' displacement-x position 24,0,12
zone history name 'dis2' displacement-x position 50,0,22
zone history name 'dis3' displacement-x position 70,0,28
zone history name 'dis4' displacement-x position 95,0,40
zone history name 'dis5' displacement-x position 125,0,50
zone history name 'dis6' displacement-z position 24,0,12
zone history name 'dis7' displacement-z position 50,0,22
zone history name 'dis8' displacement-z position 70,0,28
zone history name 'dis9' displacement-z position 95,0,40
zone history name 'dis10' displacement-z position 125,0,
model solve fluid time-total 8e7 or mechanical ratio 1e-4
model save 'h20210506002'
2021年05月06日 07点05分
1
model config fluid
zone cmodel assign mohr-coulomb
zone property bulk 3e9 shear 3.1e8 friction 45 range group 'chuqiba'
zone initialize density 2250 range group 'chuqiba'
zone property bulk 3.43e7 shear 6.32e6 friction 16 cohesion 1.078e3 range group 'jiyan'
zone initialize density 1950 range group 'jiyan'
zone property bulk 3.148e7 shear 6.03e6 friction 28 cohesion 9.8e3 range group 'weikuang1'
zone initialize density 1900 range group 'weikuang1'
zone property bulk 2.667e7 shear 5.71e6 friction 30 cohesion 9.8e3 range group 'weikuang2'
zone initialize density 1850 range group 'weikuang2'
zone property bulk 2.273e7 shear 5.4e6 friction 33 cohesion 7.84e3 range group 'weikuang3'
zone initialize density 1800 range group 'weikuang3'
zone fluid cmodel assign isotropic
zone fluid property permeability 5e-6 poros 1.5 range group 'chuqiba'
zone fluid property permeability 3.06e-12 poros 1.0 range group 'jiyan'
zone fluid property permeability 1.275e-10 poros 0.95 range group 'weikuang1'
zone fluid property permeability 3.825e-10 poros 0.9 range group 'weikuang2'
zone fluid property permeability 1.428e-10 poros 0.85 range group 'weikuang3'
zone gridpoint initialize fluid-modulus 2.2e9
zone gridpoint initialize fluid-tension 0
zone initialize fluid-density 1000
zone gridpoint fix velocity-x range position-x 0
zone gridpoint fix velocity-x range position-x 200
zone gridpoint fix velocity range position-z 0
zone gridpoint fix velocity-y
model gravity 0,0,-10
model fluid active off
model history mechanical ratio-maximum
model solve
model save 'h20210506001'
zone gridpoint initialize displacement (0,0,0)
zone gridpoint initialize velocity (0,0,0)
model fluid active on
model fluid substep 100
model mechanical substep 100
model mechanical slave on
model history name 'fluid-alltime' fluid time-total
zone history name 'chuqiba-pp1' pore-pressure position 26.5, 0, 5.0
zone history name 'chuqiba-pp2' pore-pressure position 26.5, 30, 5.0
zone history name 'chuqiba-pp3' pore-pressure position 26.5, 60, 5.0
zone history name 'chuqiba-pp4' pore-pressure position 21.58, 0, 7.76
zone history name 'chuqiba-pp5' pore-pressure position 21.58, 30, 7.76
zone history name 'chuqiba-pp6' pore-pressure position 30, 0, 2.5
zone history name 'chuqiba-pp7' pore-pressure position 30, 0, 2.5
zone history name 'pp8' pore-pressure position 70, 0, 0
zone history name 'pp9' pore-pressure position 135, 0, 0
zone history name 'pp10' pore-pressure position 195, 0, 0
zone history name 'pp11' pore-pressure position 135, 0, 20
zone history name 'pp12' pore-pressure position 180, 0, 35
zone history name 'pp13' pore-pressure position 70, 0, 20
zone history name 'pp14' pore-pressure position 115, 0, 35
zone history name 'pp15' pore-pressure position 165, 0, 50
zone history name 'pp16' pore-pressure position 50, 0, 20
zone history name 'pp17' pore-pressure position 140, 0, 55
zone history name 'dis1' displacement-x position 24,0,12
zone history name 'dis2' displacement-x position 50,0,22
zone history name 'dis3' displacement-x position 70,0,28
zone history name 'dis4' displacement-x position 95,0,40
zone history name 'dis5' displacement-x position 125,0,50
zone history name 'dis6' displacement-z position 24,0,12
zone history name 'dis7' displacement-z position 50,0,22
zone history name 'dis8' displacement-z position 70,0,28
zone history name 'dis9' displacement-z position 95,0,40
zone history name 'dis10' displacement-z position 125,0,
model solve fluid time-total 8e7 or mechanical ratio 1e-4
model save 'h20210506002'