group 支护结构 range group 隧道 cyl end1 0 0 13.516 end2 0 50 13.516 rad 13.52 cyl end1 0 0 13.516 end2 0 50 13.516 rad 13.37 not
group 支护结构 range group 隧道 cyl end1 3.3684 0 4.1354 end2 3.3684 50 4.1354 rad 3.5554 cyl end1 3.3684 0 4.1354 end2 3.3684 50 4.1354 rad 3.3900 not x 5.49 6.53 z 1.285 2.495
group 支护结构 range group 隧道 cyl end1 -3.3684 0 4.1354 end2 -3.3684 50 4.1354 rad 3.5554 cyl end1 -3.3684 0 4.1354 end2 -3.3684 50 4.1354 rad 3.3900 not x -6.53 -5.49 z 1.285 2.495
;创建显示图
plot create view_inter plot add surface
plot add interface red plot show
;清零位移 ini xdis 0 ini ydis 0 ini zdis 0
隧道埋深的影响
apply szz -5.9e6 range z 52 52.2
开挖一个group中的一段
model null range group 1 y 0 90 ;
导出指定位置的数据
set logfile 挖前节点位置.log set log on
print gp position range x -0.1 0.1 z 7.5 7.7 set log off
;set logfile 挖前位移.log ;set log on
;print gp displacement range x -0.1 0.1 z 7.5 7.7 ;set log off
1.3建模过程
1、在CAD中建模
a，模型中的所有线条都须改为直线段，不能用多段线，线段与线段之间不能重叠。线段数量不宜太多； b，将模型的高程与地形图高程对应。将模型端点坐标放置恰当位置（端点坐标最好为整数）。确保模型在第一象限。
c，将模型保存为“DXF”格式；
2、将DXF格式的模型导入algor中
a，打开algor，点击“SD2”－file－import－autocad－F10－enclose b，此时，algor中将出现模型的轮廓 c，将模型导出，file－export－iges d，此时得到igs格式的模型
新方法：
1、在CAD中将点数据导出：若地形线为多段线，则选中多段线，输入list命令，即可将多段线各点导出
2、用excel编辑点，用ansys创建点 3、用spline线创建线
3、在ansys中处理模型
a，打开ansys软件－file－import－iges；将模型导入ansys中
删除line only，创建spline线，创建大面，用divide将大面切割成小面。面被线divide时，线消失。
创建面出现问题时可能是点出现重复，merge一下即可。
b，设置单元类型：element type－add－solid－quad 4node 42－
c，设置材料属性：material props－material models;模型中有多少个group有建多少个材料属性，材料属性可选择structural-density,赋值应赋不同的值，以便merger。 d，创建面单元：create-areas-arbitrary-by lines
e，给每个面赋材料属性：meshtool-给每个面赋材料属性－element attributes-areas-set；面单元编号与材料属性编号应对应
d，网分：smart size自动网分-Tri三角形－Free自由映射－mesh（大小一致）设置line的长度-size controls-lines-set
选择等边长网分时，注意用clear命令清楚等边长的间隔，如未clear，这线为虚线。
e，合并节点和单元，重新编号：numbering ctrls-merge items-notes/elements
numbering ctrls-compress numbers-notes/elements 网分后不能合并所有的对象，只能notes/elements，如果all则所有材料合并为一个 当材料属性所附的参数不相等时，可以在网分后merge all。
f，导出节点和单元数据:create-nodes-write node file-node.dat-ok
create-elements-write elemet file-elment.dat-ok
4、利用ansys与flac的转换程序，将模型转入flac 平面模型转换-solid 42 最新转换程序：
平面模型转换-Grid文件-flac模型文件，导出的是一个模型，而非命令流，但速度很快。 打开flac，file-imput grid-已转换好的grid模型文件。
group 冲积物 range group 2
group 支护 range group 隧道
group 开挖体 range group \已开挖\
5、在flac中计算：
参数：剪切模量G=E/2(1+u),体积模量K=E/3(1-2u) 各岩石的参数： model mohr
prop bulk 28.33e7 shear 21.5e6 coh 5000 fric 30.4 dens 1950 range group 1 ;粉细砂，弹模34Mpa，泊松比0.3
prop bulk 27.778e7 shear 15.8736e7 coh 6000 fric 33.2 dens 2050 range group 2 ;粗砂，弹模40Mpa，泊松比0.26
prop bulk 2.57e9 shear 1.77e9 coh 1.2e6 fric 40 dens 2320 range group 3 ;泥岩，弹模4320Mpa，泊松比0.22
prop bulk 2.57e9 shear 1.77e9 coh 1.2e6 fric 40 dens 2320 range group 4 ;泥岩，弹模4320Mpa，泊松比0.22
prop bulk 16.67e9 shear 12.5e9 dens 2400 range group pile;C30号混凝土，弹模30Gpa，泊松比0.2
prop bulk 16.96e9 shear 11.934e9 coh 1.78e6 fric 47 dens 2530 range group 1 ;低品位石灰岩覆盖层泊松比0.215
prop bulk 0.67e9 shear 0.22e9 coh 0.03e6 fric 21 dens 2250 range group 2 ;粘土夹层泊松比0.35
prop bulk 17e9 shear 12e9 coh 1.8e6 fric 47 dens 2580 range group 3 ;石灰岩矿石层泊松比0.22
prop bulk 6.67e9 shear 2.22e9 coh 0.7e6 fric 39 dens 2450 range group 2 ;四类围岩，变模6000Mpa，泊松比0.35，内摩擦角39，粘聚力0.7Mpa，密度2450kg/m3
property she=0.09e9 bulk=0.2e9 coh=0.04e6 fric =27 ten=0.02e6 den=2020 range gro m-1;河谷堆积物
property she=0.5e9 bulk=1.4e9 coh=0.2e6 fric =27 ten=0.2e6 den=2250 range gro m-3;强风化带，5类围岩
property she=1.4e9 bulk=3.5e9 coh=0.7e6 fric =33 ten=0.45e6 den=2325 range gro m-5;弱风化带，4类围岩
property she=5.1e9 bulk=9.6e9 coh=1.1e6 fric =44.5 ten=7e6 den=2550 range gro m-7;基岩，3类围岩
property she=0.09e9 bulk=0.2e9 coh=0.04e6 fric =40 ten=0.02e6 den=2020 range gro m-1;河床冲积物
prop bulk 16.67e9 shear 12.5e9 dens 2400 range group m-1;C30号混凝土，弹模30Gpa，泊松比0.2
2-1定义一个FISH函数
new def abc
abc = 25 * 3 + 5 End print abc
2-2使用一个变量
new def abc hh = 25
abc = hh * 3 + 5 End Print hh Print abc
2-3对变量和函数的理解
new def abc hh = 25
abc = hh * 3 + 5 End
set abc=0 hh=0 print hh print abc print hh new def abc
abc = hh * 3 + 5 end
set hh=25 print abc
set abc=0 hh=0 print hh print abc print hh
2-4获取变量的历史记录
new
gen zone brick size 1 2 1
model mohr
prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10 fix x y z range y -0.1 0.1
apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot block group def get_ad
ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load
load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load
hist gp ydis 0,2,0 step 1000
plot his 1 vs -2
2-5用FISH函数计算体积模量和剪砌模量
new
def derive
s_mod = y_mod / (2.0 * (1.0 + p_ratio))
b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end
set y_mod = 5e8 p_ratio = 0.25 derive
print b_mod print s_mod
2-6 在FLAC输入中使用符号变量
New
def derive
s_mod = y_mod / (2.0 * (1.0 + p_ratio))
b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end
set y_mod = 5e8 p_ratio = 0.25 derive
gen zone brick size 2,2,2 model elastic