sel node id=124 36 -12 19 sel beamsel node 117 120 sel beamsel node 118 120 sel beamsel node 119 120 sel beamsel node 120 121 sel beamsel node 122 121 sel beamsel node 123 121 sel beamsel node 124 121 sel link id 113 117 target zone sel link id 114 118 target zone sel link id 115 119 target zone sel link id 116 122 target zone sel link id 117 123 target zone sel link id 118 124 target zone ;内撑赋属性

Sel beam prop emod=32e9 nu=0.167

sel beam prop XCArea=7.2e-1 XCIz=4.86e-2 XCIy=3.48e-2 XCJ=0.0 ;监测内撑轴力--------------------------------- h id 1 sel beams force Fx cid=1 h id 2 sel beams force Fx cid=2 h id 3 sel beams force Fx cid=3 h id 4 sel beams force Fx cid=7 h id 5 sel beams force Fx cid=8 h id 6 sel beams force Fx cid=9 ;显示结构单元 plot con zdisp

plot add sel geom node off plot add sel beam force fx plot set rot 40 0 280 Solve

Save stage2.sav ;保存工况二，方便下一工况调用 运行上面命令流，得到位移云图，如下所示：

FLAC3D 3.00Step 3625 Model Perspective10:31:45 Wed Apr 16 2014Center: X: -1.845e+001 Y: 9.395e+000 Z: -5.000e-001Dist: 3.114e+002Rotation: X: 40.000 Y: 0.000 Z: 280.000Mag.: 0.512Ang.: 22.500Contour of Y-Displacement Magfac = 0.000e+000-7.3768e-003 to -6.0000e-003-6.0000e-003 to -4.0000e-003-4.0000e-003 to -2.0000e-003-2.0000e-003 to 0.0000e+000 0.0000e+000 to 2.0000e-003 2.0000e-003 to 4.0000e-003 4.0000e-003 to 6.0000e-003 6.0000e-003 to 7.3931e-003 Interval = 2.0e-003SEL Geometry Magfac = 0.000e+000beam Force Fx Magfac = 0.000e+000 positive wrt SEL systemnegative wrt SEL system Itasca Consulting Group, Inc.Minneapolis, MN USA 图 工况二水平位移云图 ③工况三命令流：

restore stage2.sav ;第三步开挖

model null range x 0 36 y -12 12 z 6.2 12.8

;设置第二道钢 beam内撑------------------------------------- ;设置斜撑

sel beam begin 3,12,13.3 end 0,9,13.3 sel beam begin 6,12,13.3 end 0,6,13.3 sel beam begin 9,12,13.3 end 0,3,13.3 sel beam begin 3,-12,13.3 end 0,-9,13.3 sel beam begin 6,-12,13.3 end 0,-6,13.3 sel beam begin 9,-12,13.3 end 0,-3,13.3 ;设置对撑

sel beam begin 12,12,13.3 end 12,-12,13.3 sel beam begin 15,12,13.3 end 15,-12,13.3 sel beam begin 18,12,13.3 end 18,-12,13.3 sel beam begin 21,12,13.3 end 21,-12,13.3 sel beam begin 24,12,13.3 end 24,-12,13.3 sel beam begin 27,12,13.3 end 27,-12,13.3 sel beam begin 30,12,13.3 end 30,-12,13.3 sel beam begin 33,12,13.3 end 33,-12,13.3 sel beam begin 36,12,13.3 end 36,-12,13.3 ;内撑赋属性---------------------------------------

Sel beam prop emod=200e9 nu=0.3 range z 13.2 13.4

sel beam prop XCArea=3.05e-3 XCIz=1.31e-3 XCIy=1.31e-3 XCJ=0.0 range z 13.2 13.4

;显示结构单元 plot con zdisp

plot add sel geom node off plot add sel beam force fx plot set rot 40 0 280 Solve

Save stage3.sav ;保存工况三，方便下一工况调用

运行上面命令流，得到位移云图，如下所示：

FLAC3D 3.00Step 5653 Model Perspective10:35:55 Wed Apr 16 2014Center:Rotation: X: -2.040e+001 X: 40.000 Y: -1.663e+000 Y: 0.000 Z: -5.000e-001 Z: 280.000Dist: 3.114e+002Mag.: 0.512Ang.: 22.500Contour of Y-Displacement Magfac = 0.000e+000-1.6873e-002 to -1.5000e-002-1.5000e-002 to -1.0000e-002-1.0000e-002 to -5.0000e-003-5.0000e-003 to 0.0000e+000 0.0000e+000 to 5.0000e-003 5.0000e-003 to 1.0000e-002 1.0000e-002 to 1.5000e-002 1.5000e-002 to 1.6582e-002 Interval = 5.0e-003SEL Geometry Magfac = 0.000e+000beam Force Fx Magfac = 0.000e+000 positive wrt SEL systemnegative wrt SEL system Itasca Consulting Group, Inc.Minneapolis, MN USA 图 工况三水平位移云图

④工况四命令流：

restore stage3.sav ;第四步开挖

model null range x 0 36 y -12 12 z 0 6.2

;设置第三道钢 beam内撑------------------------------------- ;设置斜撑

sel beam begin 3,12,6.7 end 0,9,6.7 sel beam begin 6,12,6.7 end 0,6,6.7 sel beam begin 9,12,6.7 end 0,3,6.7 sel beam begin 3,-12,6.7 end 0,-9,6.7 sel beam begin 6,-12,6.7 end 0,-6,6.7

sel beam begin 9,-12,6.7 end 0,-3,6.7 ;设置对撑

sel beam begin 12,12,6.7 end 12,-12,6.7 sel beam begin 15,12,6.7 end 15,-12,6.7 sel beam begin 18,12,6.7 end 18,-12,6.7 sel beam begin 21,12,6.7 end 21,-12,6.7 sel beam begin 24,12,6.7 end 24,-12,6.7 sel beam begin 27,12,6.7 end 27,-12,6.7 sel beam begin 30,12,6.7 end 30,-12,6.7 sel beam begin 33,12,6.7 end 33,-12,6.7 sel beam begin 36,12,6.7 end 36,-12,6.7

;内撑赋属性---------------------------------------

Sel beam prop emod=200e9 nu=0.3 range z 13.2 13.4

sel beam prop XCArea=3.05e-3 XCIz=1.31e-3 XCIy=1.31e-3 XCJ=0.0 range z 13.2 13.4

;显示结构单元 plot con zdisp

plot add sel geom node off plot add sel beam force fx plot set rot 40 0 280 Solve

Save result.sav ;保存最后结果，通过print write n file.文件名.dat 提取数据

运行上面命令流，得到位移云图，如下所示：

FLAC3D 3.00Step 7179 Model Perspective10:39:46 Wed Apr 16 2014Center: X: -1.505e+001 Y: -6.417e-001 Z: -6.574e+000Dist: 3.114e+002Rotation: X: 40.000 Y: 0.000 Z: 280.000Mag.: 0.512Ang.: 22.500Contour of Y-Displacement Magfac = 0.000e+000-2.0750e-002 to -2.0000e-002-2.0000e-002 to -1.5000e-002-1.5000e-002 to -1.0000e-002-1.0000e-002 to -5.0000e-003-5.0000e-003 to 0.0000e+000 0.0000e+000 to 5.0000e-003 5.0000e-003 to 1.0000e-002 1.0000e-002 to 1.5000e-002 1.5000e-002 to 2.0000e-002 2.0000e-002 to 2.0678e-002 Interval = 5.0e-003SEL Geometry Magfac = 0.000e+000beam Force Fx Magfac = 0.000e+000 positive wrt SEL systemItasca Consulting Group, Inc.Minneapolis, MN USA 图 工况三水平位移云图

四、计算结果

为了更清晰的反映计算结果，寻找内在规律，可将五个工况的地下连续墙侧移绘制在一张图中，通过Origin 作图等绘图软件。很明显，只有对数据进行更深入的处理，才能发现计算结果的规律，反映问题的本质。同时也能够让读者以最直观、最清晰的方式解读计算成果。

提取数据命令流如下：

hist write 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 file 中段水平位移.dat hist write 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 file 边角水平位移.dat hist write 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 file 地表沉降.dat hist write 1 2 3 4 5 6 file 轴力.dat