- 积分
- 0
- 注册时间
- 2013-10-14
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-3-6 16:05:17
|
显示全部楼层
来自 江苏南京
程序代码如下：
n
impgrid Flacmodel_haitang.Flac3D
pl bl gr
def modu_parameter
array y_mod(3),p_ratio(3),shea_mod(3),bul_mod(3)
y_mod(1)=6.74e8 ;E强~全风化弹性模量(Pa)
y_mod(2)=1.67e9 ;E弱风化弹性模量(Pa)
y_mod(3)=2.15e9 ;E微风化弹性模量(Pa)
p_ratio(1)=0.35 ;E强~全风化泊松比
p_ratio(2)=0.28 ;E弱风化泊松比
p_ratio(3)=0.23 ;E微风化泊松比
loop ii (1,3)
shea_mod(ii)=y_mod(ii)/(2.0*(1+p_ratio(ii)))
bul_mod(ii)=y_mod(ii)/(3.0*(1-2.0*p_ratio(ii)))
endloop
shear1=shea_mod(1)
bulk1=bul_mod(1)
shear2=shea_mod(2)
bulk2=bul_mod(2)
shear3=shea_mod(3)
bulk3=bul_mod(3)
end
modu_parameter
model mohr
pro bulk=bulk1 she=shear1 co 5.0e4 fric 20 ten 0 dilation=0 range group 6 any group 21 any group 23 any group 12 any group 20 any group 2 any group 11 any group 16 any group 1 any group 3 any group 26 any
pro bulk=bulk2 she=shear2 co 7.0e4 fric 28 ten 8.75e3 dilation=0 range group 4 any group 8 any group 22 any group 24 any group 15 any group 18 any group 14 any group 13 any group 27 any group 5 any
pro bulk=bulk3 she=shear3 co 3.55e5 fric 30 ten 2e4 dilation=0 range group 7 any group 10 any group 17 any group 19 any group 9 any group 25 any
ini density=1970 range group 6 any group 21 any group 23 any group 12 any group 20 any group 2 any group 11 any group 16 any group 1 any group 3 any group 26 any
ini density=2200 range group 4 any group 8 any group 22 any group 24 any group 15 any group 18 any group 14 any group 13 any group 27 any group 5 any
ini density=2650 range group 7 any group 10 any group 17 any group 19 any group 9 any group 25 any
set grav=0,0,-9.81
;初始化材料密度: 饱和的和非饱和的
;设置边界约束条件
fix x y z range z -50.1 -49.9
fix x range x -20.1 -19.9
fix x range x 224.9 225.1
fix y range y -0.1 0.1
fix y range y 435.1 435.3
hist unbal
solve e
save gravtybalance1.sav
;位移和速度归零
ini xdisp 0
ini ydisp 0
ini zdisp 0
ini xvel 0
ini yvel 0
ini zvel 0
ini state=0
save dispvelto01.sav
hist id 2 gp xdisp 1 , 320 , 194
hist id 3 gp xdisp 1 , 340 , 194
hist id 4 gp xdisp 21 , 300 , 174
hist id 5 gp xdisp 21 , 320 , 174
hist id 6 gp xdisp 21 , 340 , 174
hist id 7 gp xdisp 51 , 280 , 154
hist id 8 gp xdisp 51 , 300 , 154
hist id 9 gp xdisp 51 , 320 , 154
hist id 10 gp xdisp 51 , 340 , 154
hist id 11 gp xdisp 69 , 220 , 134
hist id 12 gp xdisp 69 , 240 , 134
hist id 13 gp xdisp 69 , 260 , 134
hist id 14 gp xdisp 69 , 280 , 134
hist id 15 gp xdisp 69 , 300 , 134
hist id 16 gp xdisp 69 , 320 , 134
hist id 17 gp xdisp 69 , 340 , 134
hist id 18 gp xdisp 93 , 160 , 114
hist id 19 gp xdisp 93 , 180 , 114
hist id 20 gp xdisp 93 , 200 , 114
hist id 21 gp xdisp 93 , 220 , 114
hist id 22 gp xdisp 93 , 240 , 114
hist id 23 gp xdisp 93 , 260 , 114
hist id 24 gp xdisp 93 , 280 , 114
hist id 25 gp xdisp 93 , 300 , 114
hist id 26 gp xdisp 93 , 320 , 114
hist id 27 gp xdisp 93 , 340 , 114
hist id 28 gp xdisp 93 , 360 , 114
hist id 29 gp xdisp 115 , 160 , 94
hist id 30 gp xdisp 115 , 180 , 94
hist id 31 gp xdisp 115 , 200 , 94
hist id 32 gp xdisp 115 , 220 , 94
hist id 33 gp xdisp 115 , 240 , 94
hist id 34 gp xdisp 115 , 260 , 94
hist id 35 gp xdisp 115 , 280 , 94
hist id 36 gp xdisp 115 , 300 , 94
hist id 37 gp xdisp 115 , 320 , 94
hist id 38 gp xdisp 115 , 340 , 94
hist id 39 gp xdisp 115 , 360 , 94
hist id 40 gp xdisp 139 , 160 , 74
hist id 41 gp xdisp 139 , 180 , 74
hist id 42 gp xdisp 139 , 200 , 74
hist id 43 gp xdisp 139 , 220 , 74
hist id 44 gp xdisp 139 , 240 , 74
hist id 45 gp xdisp 139 , 260 , 74
hist id 46 gp xdisp 139 , 280 , 74
hist id 47 gp xdisp 139 , 300 , 74
hist id 48 gp xdisp 139 , 320 , 74
hist id 49 gp xdisp 139 , 340 , 74
hist id 50 gp xdisp 139 , 360 , 74
hist id 51 gp xdisp 161 , 160 , 54
hist id 52 gp xdisp 161 , 180 , 54
hist id 53 gp xdisp 161 , 200 , 54
hist id 54 gp xdisp 161 , 220 , 54
hist id 55 gp xdisp 161 , 240 , 54
hist id 56 gp xdisp 161 , 260 , 54
hist id 57 gp xdisp 161 , 280 , 54
hist id 58 gp xdisp 161 , 300 , 54
hist id 59 gp xdisp 161 , 320 , 54
hist id 60 gp xdisp 161 , 340 , 54
hist id 61 gp xdisp 161 , 360 , 54
hist id 62 gp xdisp 1 , 320 , 170
hist id 63 gp xdisp 21 , 320 , 150
hist id 64 gp xdisp 51 , 320 , 131
hist id 65 gp xdisp 69 , 220 , 88
hist id 66 gp xdisp 69 , 340 , 88
hist id 67 gp xdisp 115 , 160 , 48
hist id 68 gp xdisp 115 , 340 , 48
hist id 69 gp xdisp 115 , 360 , 69
hist id 70 gp xdisp 139 , 320 , 58
hist id 71 gp xdisp 161 , 160 , 14
hist id 72 gp xdisp 161 , 300 , 5
hist id 73 gp xdisp 161 , 360 , 12
hist id 74 gp zdisp 1 , 320 , 194
hist id 75 gp zdisp 1 , 340 , 194
hist id 76 gp zdisp 21 , 300 , 174
hist id 77 gp zdisp 21 , 320 , 174
hist id 78 gp zdisp 21 , 340 , 174
hist id 79 gp zdisp 51 , 280 , 154
hist id 80 gp zdisp 51 , 300 , 154
hist id 81 gp zdisp 51 , 320 , 154
hist id 82 gp zdisp 51 , 340 , 154
hist id 83 gp zdisp 69 , 220 , 134
hist id 84 gp zdisp 69 , 240 , 134
hist id 85 gp zdisp 69 , 260 , 134
hist id 86 gp zdisp 69 , 280 , 134
hist id 87 gp zdisp 69 , 300 , 134
hist id 88 gp zdisp 69 , 320 , 134
hist id 89 gp zdisp 69 , 340 , 134
hist id 90 gp zdisp 93 , 160 , 114
hist id 91 gp zdisp 93 , 180 , 114
hist id 92 gp zdisp 93 , 200 , 114
hist id 93 gp zdisp 93 , 220 , 114
hist id 94 gp zdisp 93 , 240 , 114
hist id 95 gp zdisp 93 , 260 , 114
hist id 96 gp zdisp 93 , 280 , 114
hist id 97 gp zdisp 93 , 300 , 114
hist id 98 gp zdisp 93 , 320 , 114
hist id 99 gp zdisp 93 , 340 , 114
hist id 100 gp zdisp 93 , 360 , 114
hist id 101 gp zdisp 115 , 160 , 94
hist id 102 gp zdisp 115 , 180 , 94
hist id 103 gp zdisp 115 , 200 , 94
hist id 104 gp zdisp 115 , 220 , 94
hist id 105 gp zdisp 115 , 240 , 94
hist id 106 gp zdisp 115 , 260 , 94
hist id 107 gp zdisp 115 , 280 , 94
hist id 108 gp zdisp 115 , 300 , 94
hist id 109 gp zdisp 115 , 320 , 94
hist id 110 gp zdisp 115 , 340 , 94
hist id 111 gp zdisp 115 , 360 , 94
hist id 112 gp zdisp 139 , 160 , 74
hist id 113 gp zdisp 139 , 180 , 74
hist id 114 gp zdisp 139 , 200 , 74
hist id 115 gp zdisp 139 , 220 , 74
hist id 116 gp zdisp 139 , 240 , 74
hist id 117 gp zdisp 139 , 260 , 74
hist id 118 gp zdisp 139 , 280 , 74
hist id 119 gp zdisp 139 , 300 , 74
hist id 120 gp zdisp 139 , 320 , 74
hist id 121 gp zdisp 139 , 340 , 74
hist id 122 gp zdisp 139 , 360 , 74
hist id 123 gp zdisp 161 , 160 , 54
hist id 124 gp zdisp 161 , 180 , 54
hist id 125 gp zdisp 161 , 200 , 54
hist id 126 gp zdisp 161 , 220 , 54
hist id 127 gp zdisp 161 , 240 , 54
hist id 128 gp zdisp 161 , 260 , 54
hist id 129 gp zdisp 161 , 280 , 54
hist id 130 gp zdisp 161 , 300 , 54
hist id 131 gp zdisp 161 , 320 , 54
hist id 132 gp zdisp 161 , 340 , 54
hist id 133 gp zdisp 161 , 360 , 54
;
;第1台阶（18-20）开挖完
m null range gro 4 any gro 6 any
step 300
save kw1-1.sav
;第2台阶（15-17）开挖完
m null range gro 8 any gro 21 any
step 300
save kw2-1.sav
;第3台阶（12-14）开挖完
m null range gro 22 any gro 23 any gro 7 any
step 300
save kw3-1.sav
;第4台阶（9-11）开挖完
m null range gro 10 any gro 20 any gro 12 any gro 24 any
step 300
save kw4-1.sav
;第5台阶（6-8）开挖完
m null range gro 2 any gro 17 any gro 18 any gro 11 any gro 15 any
step 300
save kw5-1.sav
;第6台阶（3-5）开挖完
m null range gro 5 any gro 16 any gro 19 any gro 14 any
step 300
save kw6-1.sav
;第7台阶（1-2）开挖完
m null range gro 1 any gro 13 any gro 9 any
step 300
save kw7-1.sav
solve
save kw7.sav
hist write 1 file hist2.txt
set log on ;通过日志文件输出单元和节点结果数据
set logfile &weiyi.log
print gp dis
hist write 2 vs 3 file 194.dat
hist write 4 5 vs 6 file 174.dat
hist write 7 8 9 vs 10 file 154.dat
hist write 11 12 13 14 15 16 vs 17 file 134.dat
hist write 18 19 20 21 22 23 24 25 26 27 vs 28 file 114.dat
hist write 29 30 31 32 33 34 35 36 37 38 vs 39 file 94.dat
hist write 40 41 42 43 44 45 46 47 48 49 vs 50 file 74.dat
hist write 51 52 53 54 55 56 57 58 59 60 vs 61 file 54.dat
set log off
|
|
IP卡
狗仔卡
发表于 2015-3-6 16:02:33
提升卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
~~
:):)