以下是我的初始平衡命令流：
group '右线1台阶初支' range group 128 any group 18 any group 63 any group 125 any group 54 any
group '右线2台阶初支' range group 17 any group 16 any group 53 any group 52 any
group '右线3台阶初支' range group 13 any group 19 any group 31 any group 49 any group 37 any
group '右线1台阶土体' range group 62 any group 127 any group 30 any group 61 any group 36 any group 42 any group 123 any
group '右线2台阶土体' range group 29 any group 35 any group 41 any group 28 any group 34 any group 40 any
group '右线3台阶土体' range group 27 any group 26 any group 32 any group 39 any group 38 any
group '右线1台阶二衬' range group 48 any group 24 any group 124 any
group '右线2台阶二衬' range group 23 any group 22 any group 47 any group 46 any
group '右线3台阶二衬' range group 20 any group 44 any
group '右线拱部环形加固层' range group 64
;******下穿段左线设置
group '左线1台阶初支' range group 136 any group 112 any group 121 any group 132 any group 76 any
group '左线2台阶初支' range group 111 any group 110 any group 75 any group 74 any
group '左线3台阶初支' range group 107 any group 95 any group 89 any group 71 any group 77 any
group '左线1台阶土体' range group 134 any group 119 any group 130 any group 100 any group 94 any group 88 any
group '左线2台阶土体' range group 99 any group 93 any group 87 any group 98 any group 92 any group 86 any
group '左线3台阶土体' range group 97 any group 96 any group 90 any group 85 any group 84 any
group '左线1台阶二衬' range group 135 any group 120 any group 106 any group 131 any group 82 any
group '左线2台阶二衬' range group 105 any group 104 any group 81 any group 80 any
group '左线3台阶二衬' range group 102 any group 78 any
group '左线拱部环形加固层' range group 122
;******1号线主体设置
group '1号线上部主体' range group 6
group '1号线左底角主体' range group 144
group '1号线中底角主体' range group 145
group '1号线右底角主体' range group 146
group '1号线左底部主体' range group 142
group '1号线右底部主体' range group 143
;******围岩设置
group '杂填土层' range group 1
group '粉土填土层' range group 2
group '砺层' range group 3
group '砾岩A层' range group 138 any group 4 any group 139 any
group '砾岩B层' range group 140 any group 141 any group 7 any group 55 any group 65 any group 113 any
group '下穿段上方加固层' range group 5
;;;;;;制定摩尔库伦模型
mo mohr;;;定义材料属性
;杂填土层def mat9
s_mod9=y_mod9/(2.0*(1.0+p_ratio9))
b_mod9=y_mod9/(3.0*(1.0-2.0*p_ratio9))
end
set y_mod9=67e6 p_ratio9=0.4
mat9
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)
pro bulk b_mod9 she s_mod9 fric 10 coh 60e3 range group '杂填土层'
;粉土填土层
def mat10
s_mod10=y_mod10/(2.0*(1.0+p_ratio10))
b_mod10=y_mod10/(3.0*(1.0-2.0*p_ratio10))
end
set y_mod10=52e7 p_ratio10=0.29
mat10
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)pro bulk b_mod10 she s_mod10 fric 30 coh 65e3 range group '粉土填土层'
;砺层def mat11
s_mod11=y_mod11/(2.0*(1.0+p_ratio11))
b_mod11=y_mod11/(3.0*(1.0-2.0*p_ratio11))
end
set y_mod11=104e7 p_ratio11=0.2
mat11
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)
pro bulk b_mod11 she s_mod11 fric 38 coh 68e3 range group '砺层'
;砾岩A层
def mat12
s_mod12=y_mod12/(2.0*(1.0+p_ratio12))
b_mod12=y_mod12/(3.0*(1.0-2.0*p_ratio12))
end
set y_mod12=2.2e9 p_ratio12=0.23
mat12
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)pro bulk b_mod12 she s_mod12 fric 40 coh 70e3 range group '砾岩A层'
;砾岩B层和下穿段def mat13
s_mod13=y_mod13/(2.0*(1.0+p_ratio13))
b_mod13=y_mod13/(3.0*(1.0-2.0*p_ratio13))
end
set y_mod13=2.8e9 p_ratio13=0.23
mat13
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)
pro bulk b_mod13 she s_mod13 fric 40 coh 70e3 range group '砾岩B层' any group '右线1台阶初支' any group '右线2台阶初支' any group '右线3台阶初支' any group '右线1台阶土体' any group '右线2台阶土体' any group '右线3台阶土体' any group '右线1台阶二衬' any group '右线2台阶二衬' any group '右线3台阶二衬' any group '左线1台阶初支' any group '左线2台阶初支' any group '左线3台阶初支' any group '左线1台阶土体' any group '左线2台阶土体' any group '左线3台阶土体' any group '左线1台阶二衬' any group '左线2台阶二衬' any group '左线3台阶二衬' any
;下穿段上方加固层
def mat14
s_mod14=y_mod14/(2.0*(1.0+p_ratio14))
b_mod14=y_mod14/(3.0*(1.0-2.0*p_ratio14))
end
set y_mod14=4.2e9 p_ratio14=0.21
mat14
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)pro bulk b_mod14 she s_mod14 range group '下穿段上方加固层'
;一号线def mat15
s_mod15=y_mod15/(2.0*(1.0+p_ratio15))
b_mod15=y_mod15/(3.0*(1.0-2.0*p_ratio15))
end
set y_mod15=2.5e10 p_ratio15=0.2
mat15
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)
pro bulk b_mod15 she s_mod15 range group '1号线上部主体' any group '1号线左底角主体' any group '1号线中底角主体' any group '1号线右底角主体' any group '1号线左底部主体' any group '1号线右底部主体' any
;下穿段左右线上方两环形加固层
def mat16
s_mod16=y_mod16/(2.0*(1.0+p_ratio16))
b_mod16=y_mod16/(3.0*(1.0-2.0*p_ratio16))
end
set y_mod16=4.2e9 p_ratio16=0.21
mat16
; 围岩：E=1.5e9,v=0.4,K=E/3(1-2v),G=E/2(1+v)pro bulk b_mod16 she s_mod16 range group '右线拱部环形加固层' any group '左线拱部环形加固层' any
;;;;;;;限制边界条件
fix x range x -30.01 -29.99
fix x range x 29.99 30.01
fix z range z -10.01 -9.99
fix y range y -0.01 0.01
fix y range y 23.19 23.21
apply szz -2e4 range z 34.99 35.01;;;;;;;加载地面均布荷载
apply szz -2.8e4 range z 24.537 24.539 group '1号线左底部主体' any group '1号线右底部主体' any
;;;;;;;设置土层和1号线主体的密度
ini density 1650 range group '杂填土层'
ini density 1950 range group '粉土填土层'
ini density 2050 range group '砺层'
ini density 2130 range group '砾岩A层'
ini density 2130 range group '砾岩B层' any group '右线1台阶初支' any group '右线2台阶初支' any group '右线3台阶初支' any group '右线1台阶土体' any group '右线2台阶土体' any group '右线3台阶土体' any group '右线1台阶二衬' any group '右线2台阶二衬' any group '右线3台阶二衬' any group '右线1台阶初支' any group '左线2台阶初支' any group '左线3台阶初支' any group '左线1台阶土体' any group '左线2台阶土体' any group '左线3台阶土体' any group '左线1台阶二衬' any group '左线2台阶二衬' any group '左线3台阶二衬' any
ini density 2400 range group '下穿段上方加固层'
ini density 2400 range group '右线拱部环形加固层' any group '左线拱部环形加固层' any
ini density 2500 range group '1号线上部主体' any group '1号线左底角主体' any group '1号线中底角主体' any group '1号线右底角主体' any group '1号线左底部主体' any group '1号线右底部主体' any
set gravity 0 0 -9.8
solve
save 3T-ini.sav
