第1部分 命令流按照顺序进行
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
t abc=0 hh=0
print hh
print abc
print hh
new
def abc
abc = hh * 3 + 5
end
t hh=25
print abc
t 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 t 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
t 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
t y_mod = 5e8 p_ratio = 0.25
derive
gen zone brick size 2,2,2
陈大发
model elastic
prop bulk=b_mod shear=s_mod
print zone prop bulk
print zone prop shear
2-7 控制循环
New
行车记录仪十大品牌def xxx
sum = 0
prod = 1
loop n (1,10)
sum = sum + n
prod = prod * n
end_loop
end
xxx
print sum, prod
new
gen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10)
model elas
plot t rotation 0 0 45
plot block group
def install
pnt = zone_head
loop while pnt #null
z_depth = -z_zcen(pnt)
y_mod = y_zero + cc * sqrt(z_depth)
z_prop(pnt, ’shear’) = y_mod / (2.0*(1.0+p_ratio))
z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio))
pnt = zone_next(pnt)
end_loop
end
set p_ratio=0.25 y_zero=1e7 cc=1e8
install
2-8 拆分命令行
new ;example of a sum of many things
def long_sum
猝然的意思temp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10
long_sum = temp + v11 + v12 + v13 + v14 + v15
end
2-9 变量类型
new
def haveone
aa = 2
bb = 3.4
cc = ’Have a nice day’
气概的意思dd = aa * bb
ee = cc + ’, old chap’
end
haveone
print fish
2-10 IF条件语句
new
def abc
if xx > 0 then
abc = 1000
el
abc = -1000
end_if
钻石怎么画end
t xx = 10
print abc
t xx = 0
print abc
2-11 索单元自动生成
new
gen zone brick size 10 3 5
plot t rotation 0 0 45
plot block group
def place_cables
loop n (1,5)
z_d = float(n) - 0.5
command
l cable beg 0.0,1.5,z_d end 7.0,1.5,z_d ng 7
end_command
end_loop
end
place_cables
plot grid l geom red
new
gen zone brick size 10 3 5
plot t rotation 15 0 60
plot block group
mod mohr
prop bulk 1e8 shear .3e8 fric 35
prop coh 1e3 tens 1e3
ini dens 1000
t grav 0,0,-10
fix x y z range z -.1 .1
fix y range y -.1 .1
fix y range y 2.9 3.1
fix x range x -.1 .1
fix x range x 9.9 10.1
t large
hist unbal
solve
save cab_str.sav
ini xdis 0 ydis 0 zdis 0
hist gp xdisp 0,1,5
def place_cables
loop n (1,5)
z_d = 5.5 - float(n)
z_t = z_d + 0.5
z_b = z_d - 0.5
command
free x range x -.1,.1 z z_b z_t
solve
l cable beg 0.0,0.5,z_d end 7.0,0.5,z_d ng 7
l cable beg 0.0,1.5,z_d end 7.0,1.5,z_d ng 7
l cable beg 0.0,2.5,z_d end 7.0,2.5,z_d ng 7
l cable prop emod 2e10 ytension 1e8 xcarea 1.0 &
gr_k 2e10 gr_coh 1e10 gr_per 1.0
end_command
end_loop
end
place_cables
save cab_end.sav
plot sketch l cable force red
2-12圆形隧道开挖模拟计算
;建立模型
gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 &
size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderock
gen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &
size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock
gen zon reflect dip 90 dd 90 orig 0 0 0
gen zon reflect dip 0 dd 0 ori 0 0 0
gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1
gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2
gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3
gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3
gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4
gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5
gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1
gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1
