考虑地下水的隧道实例（含子网格技术）

CQMaxwell

;考虑耦合
;###1、建立力学模型定义参数###
config gwflow ats extra 5
gr 93 70
gen 0,0 0,40 40,40 40,0 i=1,14 j=1,21
m e i=1,13 j=1,20
gen 40,0 40,40 72,40 72,0 i=14,30 j=1,21
m e i=14,29 j=1,20
gen 72,0 72,40 112,40 112,0 i=30,45 j=1,21
m e i=30,44 j=1,20
gen 0,40 0,72 40,72 40,40 i=1,14 j=21,41
m e i=1,13 j=21,40
gen 40,40 40,72 72,72 72,40 i=46,94 j=1,41
m e i=46,93 j=1,40
gen 72,40 72,72 112,72 112,40 i=30,45 j=21,41
m e i=30,44 j=21,40
gen 0,72 0,160 40,152.86 40,72 i=1,14 j=41,71
m e i=1,13 i=41,70
gen 40,72 40,152.86 72,147.14 72,72 i=14,30  j=41,71
m e i=14,29 j=41,70
gen 72,72 72,147.14 112,140 112,72 i=30,45  j=41,71
m e i=30,44 j=41 70
;粘结节点
attach aside from 46 1 to 46 41 bside from 14 21 to 14 41
attach aside from 46 41 to 94 41 bside from 14 41 to 30 41
attach aside from 30 41 to 30 21 bside from 94 41 to 94 1
attach aside from 14 21 to 30 21 bside from 46 1 to 94 1
;定义参数
group 'ground' notnull
model mohr notnull group 'ground'
prop dens=2500 shear=5.64e8 bulk=1.47e9 fric=23 coh=1.0e5
prop por=0.3 perm=4.5E-10 i=46,94 j=1,41;渗流参数
prop por=0.3 perm=4.5E-10 i=1,14 j=21,41
prop por=0.3 perm=4.5E-10 i=30,45 j=21,41
;开挖界线
gen arc 56.00,51.60 62.20,51.60 180
gen arc 54.74,51.74 49.81,52.01 44
gen arc 56.00,62.12 50.76,48.83 44
gen arc 56.00,52.31 61.09,48.77 44
;开挖线
gen arc 56.00,62.40 51.18,55.52 70
;位移边界
fix y i=1,45 j=1
fix x i=1
fix x i=45
;关闭渗流模式只进行力学模式
set grav=9.81
set flow=off
set large
;计算初始位移
solve ela
save initial.sav  
;########2\饱和流固耦合分析#########
rest initial.sav
set flow=on ;开启渗流模式
water bulk=10000.0
water density=1000.0
initial pp 320000.0 var 0.0,-320000.0 i=46,94 j=1,41
initial pp 320000.0 var 0.0,-320000.0 i=1,14 j=21,41
initial pp 320000.0 var 0.0,-320000.0 i=30,45 j=21,41
fix pp i=46,94 j=41
fix pp i=1,14 j=41
fix pp i=30,45 j=41
ini pp 0 i=46,94 j=41
ini pp 0 i=1,14 j=41
ini pp 0 i=30,45 j=41
fix pp mark
ini pp 0 mark
;########\进行力学计算微台阶法#######
rest ini2.save
ini xdisp=0 ydisp=0
;隧道监测拱顶沉降和周边收敛
hist unbal ;hist 1
hist xdisp i=63 j=21 ;hist 2
hist xdisp i=77 j=21 ;hist 3
hist xdisp i=61 j=15 ;hist 4
hist xdisp i=79 j=15 ;hist 5
hist ydisp i=70 j=23 ;hist 6
;1#侧线
def closure1
    closure1=100.0*(xdisp(63,21)-xdisp(77,21))/(x(77,21)-x(63,21))
end
hist closure1 ;hist 7
;2#测线
def closure2
    closure2=100.0*(xdisp(61,15)-xdisp(79,15))/(x(79,15)-x(61,15))
end
hist closure2 ;hist 8
;3#测线
def closure3
    closure3=100.0*(ydisp(70,23))/(y(70,23)-y(70,11))
end
hist closure3 ;hist 9
;隧道上台阶开挖
mod null reg=69 20
step 60
;隧道上台阶初期支护
struct prop 1 dens=2400 e=2.4e10 thick=0.22 area=0.22 pr=0.25 i=8.87e-4
              & shape=0.3 syield=1.1e6 sycomp=1.2e7

struct liner beg gr 70 23 end gr 69 23 seg 3 prop 1
struct liner beg gr 69 23 end gr 68 23 seg 3 prop 1
struct liner beg gr 68 23 end gr 67 23 seg 3 prop 1
struct liner beg gr 67 23 end gr 66 23 seg 3 prop 1
struct liner beg gr 66 23 end gr 66 22 seg 3 prop 1
struct liner beg gr 66 22 end gr 65 22 seg 3 prop 1
struct liner beg gr 65 22 end gr 64 22 seg 3 prop 1
struct liner beg gr 64 22 end gr 64 21 seg 3 prop 1
struct liner beg gr 64 21 end gr 63 21 seg 3 prop 1
struct liner beg gr 63 21 end gr 63 20 seg 3 prop 1
struct liner beg gr 77 20 end gr 77 21 seg 3 prop 1
struct liner beg gr 77 21 end gr 76 21 seg 3 prop 1
struct liner beg gr 76 21 end gr 76 22 seg 3 prop 1
struct liner beg gr 76 22 end gr 75 22 seg 3 prop 1
struct liner beg gr 75 22 end gr 74 22 seg 3 prop 1
struct liner beg gr 74 22 end gr 74 23 seg 3 prop 1
struct liner beg gr 74 23 end gr 73 23 seg 3 prop 1
struct liner beg gr 73 23 end gr 72 23 seg 3 prop 1
struct liner beg gr 72 23 end gr 71 23 seg 3 prop 1
struct liner beg gr 71 23 end gr 70 23 seg 3 prop 1
;接触面
int 1 as fr node 31 to node 32 bs long fr 63 20 to 77 20
interface 1 kn=1e8 ks=1e8 fric=14 coh=1.0e5
;隧道锚杆3m
;锚杆参数
struct prop 2 e=2.1e11 a=4.91e-4 yield=3.06e9 kb=4.94e8 sb=2.64e4 spac 1.0
              & sfric=34.77 ycomp=3.06e9 yield=30.06e8 perimeter=7.85e-2
;隧道上台阶锚杆关键点
struct node 62  56.1516,60.8143  
struct node 63 54.7860,60.7353
struct node 64 53.3367,60.4223
struct node 65  52.1592,59.9771
struct node 66 50.9474,59.3071
struct node 67 49.9123,58.5186
struct node 68 49.0478,57.6493
struct node 69 57.3663,60.7137
struct node 70 58.6675,60.4211
struct node 71 59.8500,59.9729
struct node 72 61.0525,59.3071
struct node 73 62.0910,58.5157
struct node 74 62.9587,57.6419
;隧道上台阶锚杆
struct cable beg node 1 end node 62 seg 4 prop 2
struct cable beg node 5 end node 63 seg 4 prop 2
struct cable beg node 9 end node 64 seg 4 prop 2
struct cable beg node 14 end node 65 seg 4 prop 2
struct cable beg node 19 end node 66 seg 4 prop 2
struct cable beg node 24 end node 67 seg 4 prop 2
struct cable beg node 29 end node 68 seg 4 prop 2
struct cable beg node 58 end node 69 seg 4 prop 2
struct cable beg node 54 end node 70 seg 4 prop 2
struct cable beg node 49 end node 71 seg 4 prop 2
struct cable beg node 44 end node 72 seg 4 prop 2
struct cable beg node 39 end node 73 seg 4 prop 2
struct cable beg node 34 end node 74 seg 4 prop 2
step 1000
;隧道下台阶开挖
mod null reg=70 13
;隧道下台阶初期支护
struct liner beg gr 63 20 end gr 62 20 seg 3 prop 1
struct liner beg gr 62 20 end gr 62 19 seg 3 prop 1
struct liner beg gr 62 19 end gr 61 19 seg 3 prop 1
struct liner beg gr 61 19 end gr 61 18 seg 3 prop 1
struct liner beg gr 61 18 end gr 61 17 seg 3 prop 1
struct liner beg gr 61 17 end gr 61 16 seg 3 prop 1
struct liner beg gr 61 16 end gr 61 15 seg 3 prop 1
struct liner beg gr 61 15 end gr 61 14 seg 3 prop 1
struct liner beg gr 61 14 end gr 61 13 seg 3 prop 1
struct liner beg gr 61 13 end gr 62 13 seg 3 prop 1
struct liner beg gr 62 13 end gr 62 12 seg 3 prop 1
struct liner beg gr 62 12 end gr 63 12 seg 3 prop 1
struct liner beg gr 63 12 end gr 64 12 seg 3 prop 1
struct liner beg gr 64 12 end gr 64 11 seg 3 prop 1
struct liner beg gr 64 11 end gr 65 11 seg 3 prop 1
struct liner beg gr 65 11 end gr 66 11 seg 3 prop 1
struct liner beg gr 66 11 end gr 67 11 seg 3 prop 1
struct liner beg gr 67 11 end gr 68 11 seg 3 prop 1
struct liner beg gr 68 11 end gr 69 11 seg 3 prop 1
struct liner beg gr 69 11 end gr 70 11 seg 3 prop 1
struct liner beg gr 70 11 end gr 71 11 seg 3 prop 1
struct liner beg gr 71 11 end gr 72 11 seg 3 prop 1
struct liner beg gr 72 11 end gr 73 11 seg 3 prop 1
struct liner beg gr 73 11 end gr 74 11 seg 3 prop 1
struct liner beg gr 74 11 end gr 75 11 seg 3 prop 1
struct liner beg gr 75 11 end gr 76 11 seg 3 prop 1
struct liner beg gr 76 11 end gr 76 12 seg 3 prop 1
struct liner beg gr 76 12 end gr 77 12 seg 3 prop 1
struct liner beg gr 77 12 end gr 78 12 seg 3 prop 1
struct liner beg gr 78 12 end gr 78 13 seg 3 prop 1
struct liner beg gr 78 13 end gr 79 13 seg 3 prop 1
struct liner beg gr 79 13 end gr 79 14 seg 3 prop 1
struct liner beg gr 79 14 end gr 79 15 seg 3 prop 1
struct liner beg gr 79 15 end gr 79 16 seg 3 prop 1
struct liner beg gr 79 16 end gr 79 17 seg 3 prop 1
struct liner beg gr 79 17 end gr 79 18 seg 3 prop 1
struct liner beg gr 79 18 end gr 79 19 seg 3 prop 1
struct liner beg gr 79 19 end gr 78 19 seg 3 prop 1
struct liner beg gr 78 19 end gr 78 20 seg 3 prop 1
struct liner beg gr 78 20 end gr 77 20 seg 3 prop 1
;接触面
int 2 as fr node 32 to node 114 bs long fr 77 20 to 63 20
interface 2 kn=1e8 ks=1e8 fric=14 coh=1.0e5
;隧道下台阶锚杆关键点
struct node 233 48.2920,56.6511
struct node 234 47.6636 , 55.5283
struct node 235 47.2888 ,54.6070
struct node 236 46.9338, 53.2525
struct node 237 47.0192 ,51.3537
struct node 238 47.5788 ,49.7201
struct node 239 63.7144 ,56.6414
struct node 240 64.3360, 55.5292
struct node 241 64.7120,54.6048
struct node 242 65.0660 ,53.2538
struct node 243 64.9802 ,51.3051
struct node 244 64.4658 ,49.7375  
struct node 245 48.5139 ,48.2242
struct node 246 63.7440,48.4610
;下台阶锚杆
struct cable beg node 116 end node 233 seg 4 prop 2
struct cable beg node 120 end node 234 seg 4 prop 2
struct cable beg node 124 end node 235 seg 4 prop 2
struct cable beg node 128 end node 236 seg 4 prop 2
struct cable beg node 131 end node 237 seg 4 prop 2
struct cable beg node 135 end node 238 seg 4 prop 2
struct cable beg node 230 end node 239 seg 4 prop 2
struct cable beg node 226 end node 240 seg 4 prop 2
struct cable beg node 222 end node 241 seg 4 prop 2
struct cable beg node 218 end node 242 seg 4 prop 2
struct cable beg node 215 end node 243 seg 4 prop 2
struct cable beg node 211 end node 244 seg 4 prop 2
struct cable beg node 140 end node 245 seg 4 prop 2
struct cable beg node 206 end node 246 seg 4 prop 2
solve

几何模型：
   
子网格（subgrid）：

xautzhao

谢谢啦，呵呵

tangjincai28

非常感谢无私的楼主

深夜来访电话

好，认真研究一下
谢谢楼主的大公无私

zoubaoping

谢谢，会好好研究的

yang200909

谢谢楼主！！

lsssl

好东西啊，多谢

huhai

1# CQMaxwell


非常感谢

abblcn

好好学习一下。

bluesea7910

非常感谢无私的楼主

nuaaer

好东西，感谢感谢~~

tjudamer

谢谢楼主 好好学习一下

amandcat

楼主真是一个有心人啊！

gzlhai

楼主真是一个有心人啊！

liujia1985546

看看，值得学习

江beyond

非常感谢楼主的无私奉献！！！

zhuhaha

相当的经典啊

luxiang

好多命令都看不明白的

sjzdh

不错，谢谢提供分享

牙密

flac不识别第三行的ats，怎么回事