FLAC3D做动力响应分析命令流

title 库木吐喇动力响应分析
config dyn
set dyn off

model elastic

prop bulk=1.275e9 shear=1.121e9

set gra=0,0,-9.81
ini density 2500

;边界条件
fix x y z range z 0
fix x range x 0
fix x range x 65
fix y range y 0
fix y range y 60

hist id=1 unbal
solve

save 初始地应力.sav

ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0

model mohr

;定义材料属性灰岩
prop bulk=8.3e7 shear=3.85e7 ten=1.17e6 coh=27.2e6 fric=27.8

solve
save 静力分析摩尔库伦模型.sav


;动力学计算
restore 初始地应力.sav

set dyn on
set dyn multi on
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0

;去掉模型底部原有的静力边界

;地震波输入
table 1 read acc1.txt
table 2 read acc2.txt
table 3 read acc3.txt

free x y z range z -0.1 0.1
apply nquiet squiet dquiet ran z -0.1 0.1
apply yacc 1.0 hist table 1 range z -0.1 0.1
apply zacc 1.0 hist table 3 range z -0.1 0.1

apply ff
group 一维自由场网格
group 二维自由场网格 ran x 0 65
group 二维自由场网格 ran y 0 60
group 主体网格 ran x 0 65 y 0 60

;设置阻尼
set dyn damp rayleigh 0.05 9


set dyn time=0
hist reset ;清空里面所有的监测点
hist id=1 unbal

hist id=2 gp disp 0 30 24.6
hist id=3 gp disp 8.1089 30 25
hist id=4 gp disp 14.2921 30 28
hist id=5 gp disp 22.3255 30 29
hist id=6 gp disp 24.6548 30 30
hist id=7 gp disp 28.0801 30 34
hist id=8 gp disp 30.0305 30 40.2326
hist id=9 gp disp 31.655 30 64
hist id=10 gp disp 34.7475 30 69
hist id=11 gp disp 37.9218 30 72
hist id=12 gp disp 43.846 30 72.6824
hist id=13 gp disp 46.6029 30 73
hist id=14 gp disp 55.3267 30 74
hist id=15 gp disp 59.6752 30 75
hist id=16 gp disp 62.3912 30 75.2851
hist id=17 gp disp 65 30 75.6
hist id=18 gp yacc 0 30 24.6
hist id=19 gp yacc 8.1089 30 25
hist id=20 gp yacc 14.2921 30 28
hist id=21 gp yacc 22.3255 30 29
hist id=22 gp yacc 24.6548 30 30
hist id=23 gp yacc 28.0801 30 34
hist id=24 gp yacc 30.0305 30 40.2326
hist id=25 gp yacc 31.655 30 64
hist id=26 gp yacc 34.7475 30 69
hist id=27 gp yacc 37.9218 30 72
hist id=28 gp yacc 43.846 30 72.6824
hist id=29 gp yacc 46.6029 30 73
hist id=30 gp yacc 55.3267 30 74
hist id=31 gp yacc 59.6752 30 75
hist id=32 gp yacc 62.3912 30 75.2851
hist id=33 gp yacc 65 30 75.6
hist id=34 gp yvel 0 30 24.6
hist id=35 gp yvel 8.1089 30 25
hist id=36 gp yvel 14.2921 30 28
hist id=37 gp yvel 22.3255 30 29
hist id=38 gp yvel 24.6548 30 30
hist id=39 gp yvel 28.0801 30 34
hist id=40 gp yvel 30.0305 30 40.2326
hist id=41 gp yvel 31.655 30 64
hist id=42 gp yvel 34.7475 30 69
hist id=43 gp yvel 37.9218 30 72
hist id=44 gp yvel 43.846 30 72.6824
hist id=45 gp yvel 46.6029 30 73
hist id=46 gp yvel 55.3267 30 74
hist id=47 gp yvel 59.6752 30 75
hist id=48 gp yvel 62.3912 30 75.2851
hist id=49 gp yvel 65 30 75.6
hist id=50 gp za

cc 0 30 24.6
hist id=51 gp zacc 8.1089 30 25
hist id=52 gp zacc 14.2921 30 28
hist id=53 gp zacc 22.3255 30 29
hist id=54 gp zacc 24.6548 30 30
hist id=55 gp zacc 28.0801 30 34
hist id=56 gp zacc 30.0305 30 40.2326
hist id=57 gp zacc 31.655 30 64
hist id=58 gp zacc 34.7475 30 69
hist id=59 gp zacc 37.9218 30 72
hist id=60 gp zacc 43.846 30 72.6824
hist id=61 gp zacc 46.6029 30 73
hist id=62 gp zacc 55.3267 30 74
hist id=63 gp zacc 59.6752 30 75
hist id=64 gp zacc 62.3912 30 75.2851
hist id=65 gp zacc 65 30 75.6
hist id=66 gp zvel 0 30 24.6
hist id=67 gp zvel 8.1089 30 25
hist id=68 gp zvel 14.2921 30 28
hist id=69 gp zvel 22.3255 30 29
hist id=70 gp zvel 24.6548 30 30
hist id=71 gp zvel 28.0801 30 34
hist id=72 gp zvel 30.0305 30 40.2326
hist id=73 gp zvel 31.655 30 64
hist id=74 gp zvel 34.7475 30 69
hist id=75 gp zvel 37.9218 30 72
hist id=76 gp zvel 43.846 30 72.6824
hist id=77 gp zvel 46.6029 30 73
hist id=78 gp zvel 55.3267 30 74
hist id=79 gp zvel 59.6752 30 75
hist id=80 gp zvel 62.3912 30 75.2851
hist id=81 gp zvel 65 30 75.6
hist id=82 gp xdisp 0 30 24.6
hist id=83 gp xdisp 8.1089 30 25
hist id=84 gp xdisp 14.2921 30 28
hist id=85 gp xdisp 22.3255 30 29
hist id=86 gp xdisp 24.6548 30 30
hist id=87 gp xdisp 28.0801 30 34
hist id=88 gp xdisp 30.0305 30 40.2326
hist id=89 gp xdisp 31.655 30 64
hist id=90 gp xdisp 34.7475 30 69
hist id=91 gp xdisp 37.9218 30 72
hist id=92 gp xdisp 43.846 30 72.6824
hist id=93 gp xdisp 46.6029 30 73
hist id=94 gp xdisp 55.3267 30 74
hist id=95 gp xdisp 59.6752 30 75
hist id=96 gp xdisp 62.3912 30 75.2851
hist id=97 gp xdisp 65 30 75.6
hist id=98 gp ydisp 0 30 24.6
hist id=99 gp ydisp 8.1089 30 25
hist id=100 gp ydisp 14.2921 30 28
hist id=101 gp ydisp 22.3255 30 29
hist id=102 gp ydisp 24.6548 30 30
hist id=103 gp ydisp 28.0801 30 34
hist id=104 gp ydisp 30.0305 30 40.2326
hist id=105 gp ydisp 31.655 30 64
hist id=106 gp ydisp 34.7475 30 69
hist id=107 gp ydisp 37.9218 30 72
hist id=108 gp ydisp 43.846 30 72.6824
hist id=109 gp ydisp 46.6029 30 73
hist id=110 gp ydisp 55.3267 30 74
hist id=111 gp ydisp 59.6752 30 75
hist id=112 gp ydisp 62.3912 30 75.2851
hist id=113 gp ydisp 65 30 75.6
hist id=114 gp zdisp 0 30 24.6
hist id=115 gp zdisp 8.1089 30 25
hist id=116 gp zdisp 14.2921 30 28
hist id=117 gp zdisp 22.3255 30 29
hist id=118 gp zdisp 24.6548 30 30
hist id=119 gp zdisp 28.0801 30 34
hist id=120 gp zdisp 30.0305 30 40.2326
hist id=121 gp zdisp 31.655 30 64
hist id=122 gp zdisp 34.7475 30 69
hist id=123 gp zdisp 37.9218 30 72
hist id=124 gp zdisp 43.846 30 72.6824
hist id=125 gp zdisp 46.6029 30 73
hist id=126 gp zdisp 55.3267 30 74
hist id=127 gp zdisp 59.6752 30 75
hist id=128 gp zdisp 62.3912 30 75.2851
hist id=129 gp zdisp 65 30 75.6
hist id=130 gp force 0 30 24.6
hist id=131 gp force 8.1089 30 25
hist id=132 gp force 14.2921 30 28
hist id=133 gp force 22.3255 30 29
hist id=134 gp force 24.6548 30 30
hist id=135 gp for

ce 28.0801 30 34
hist id=136 gp force 30.0305 30 40.2326
hist id=137 gp force 31.655 30 64
hist id=138 gp force 34.7475 30 69
hist id=139 gp force 37.9218 30 72
hist id=140 gp force 43.846 30 72.6824
hist id=141 gp force 46.6029 30 73
hist id=142 gp force 55.3267 30 74
hist id=143 gp force 59.6752 30 75
hist id=144 gp force 62.3912 30 75.2851
hist id=145 gp force 65 30 75.6
hist id=146 gp xforce 0 30 24.6
hist id=147 gp xforce 8.1089 30 25
hist id=148 gp xforce 14.2921 30 28
hist id=149 gp xforce 22.3255 30 29
hist id=150 gp xforce 24.6548 30 30
hist id=151 gp xforce 28.0801 30 34
hist id=152 gp xforce 30.0305 30 40.2326
hist id=153 gp xforce 31.655 30 64
hist id=154 gp xforce 34.7475 30 69
hist id=155 gp xforce 37.9218 30 72
hist id=156 gp xforce 43.846 30 72.6824
hist id=157 gp xforce 46.6029 30 73
hist id=158 gp xforce 55.3267 30 74
hist id=159 gp xforce 59.6752 30 75
hist id=160 gp xforce 62.3912 30 75.2851
hist id=161 gp xforce 65 30 75.6
hist id=162 gp yforce 0 30 24.6
hist id=163 gp yforce 8.1089 30 25
hist id=164 gp yforce 14.2921 30 28
hist id=165 gp yforce 22.3255 30 29
hist id=166 gp yforce 24.6548 30 30
hist id=167 gp yforce 28.0801 30 34
hist id=168 gp yforce 30.0305 30 40.2326
hist id=169 gp yforce 31.655 30 64
hist id=170 gp yforce 34.7475 30 69
hist id=171 gp yforce 37.9218 30 72
hist id=172 gp yforce 43.846 30 72.6824
hist id=173 gp yforce 46.6029 30 73
hist id=174 gp yforce 55.3267 30 74
hist id=175 gp yforce 59.6752 30 75
hist id=176 gp yforce 62.3912 30 75.2851
hist id=177 gp yforce 65 30 75.6
hist id=178 gp zforce 0 30 24.6
hist id=179 gp zforce 8.1089 30 25
hist id=180 gp zforce 14.2921 30 28
hist id=181 gp zforce 22.3255 30 29
hist id=182 gp zforce 24.6548 30 30
hist id=183 gp zforce 28.0801 30 34
hist id=184 gp zforce 30.0305 30 40.2326
hist id=185 gp zforce 31.655 30 64
hist id=186 gp zforce 34.7475 30 69
hist id=187 gp zforce 37.9218 30 72
hist id=188 gp zforce 43.846 30 72.6824
hist id=189 gp zforce 46.6029 30 73
hist id=190 gp zforce 55.3267 30 74
hist id=191 gp zforce 59.6752 30 75
hist id=192 gp zforce 62.3912 30 75.2851
hist id=193 gp zforce 65 30 75.6
hist id=194 zone smax 0 30 24.6
hist id=195 zone smax 8.1089 30 25
hist id=196 zone smax 14.2921 30 28
hist id=197 zone smax 22.3255 30 29
hist id=198 zone smax 24.6548 30 30
hist id=199 zone smax 28.0801 30 34
hist id=200 zone smax 30.0305 30 40.2326
hist id=201 zone smax 31.655 30 64
hist id=202 zone smax 34.7475 30 69
hist id=203 zone smax 37.9218 30 72
hist id=204 zone smax 43.846 30 72.6824
hist id=205 zone smax 46.6029 30 73
hist id=206 zone smax 55.3267 30 74
hist id=207 zone smax 59.6752 30 75
hist id=208 zone smax 62.3912 30 75.2851
hist id=209 zone smax 65 30 75.6
hist id=210 zone smin 0 30 24.6
hist id=211 zone smin 8.1089 30 25
hist id=212 zone smin 14.2921 30 28
hist id=213 zone smin 22.3255 30 29
hist id=214 zone smin 24.6548 30 30
hist id=215 zone smin 28.0801 30 34
hist id

=216 zone smin 30.0305 30 40.2326
hist id=217 zone smin 31.655 30 64
hist id=218 zone smin 34.7475 30 69
hist id=219 zone smin 37.9218 30 72
hist id=220 zone smin 43.846 30 72.6824
hist id=221 zone smin 46.6029 30 73
hist id=222 zone smin 55.3267 30 74
hist id=223 zone smin 59.6752 30 75
hist id=224 zone smin 62.3912 30 75.2851
hist id=225 zone smin 65 30 75.6
hist id=226 zone sxx 0 30 24.6
hist id=227 zone sxx 8.1089 30 25
hist id=228 zone sxx 14.2921 30 28
hist id=229 zone sxx 22.3255 30 29
hist id=230 zone sxx 24.6548 30 30
hist id=231 zone sxx 28.0801 30 34
hist id=232 zone sxx 30.0305 30 40.2326
hist id=233 zone sxx 31.655 30 64
hist id=234 zone sxx 34.7475 30 69
hist id=235 zone sxx 37.9218 30 72
hist id=236 zone sxx 43.846 30 72.6824
hist id=237 zone sxx 46.6029 30 73
hist id=238 zone sxx 55.3267 30 74
hist id=239 zone sxx 59.6752 30 75
hist id=240 zone sxx 62.3912 30 75.2851
hist id=241 zone sxx 65 30 75.6
hist id=242 zone szz 0 30 24.6
hist id=243 zone szz 8.1089 30 25
hist id=244 zone szz 14.2921 30 28
hist id=245 zone szz 22.3255 30 29
hist id=246 zone szz 24.6548 30 30
hist id=247 zone szz 28.0801 30 34
hist id=248 zone szz 30.0305 30 40.2326
hist id=249 zone szz 31.655 30 64
hist id=250 zone szz 34.7475 30 69
hist id=251 zone szz 37.9218 30 72
hist id=252 zone szz 43.846 30 72.6824
hist id=253 zone szz 46.6029 30 73
hist id=254 zone szz 55.3267 30 74
hist id=255 zone szz 59.6752 30 75
hist id=256 zone szz 62.3912 30 75.2851
hist id=257 zone szz 65 30 75.6


hist id=258 dytime

set hist_rep 100 ;设置历史记录时间间隔
solve age 10 ;进行10s的动力计算
save 动力分析.sav
plot hist 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 vs 258
hist write 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 vs 137 file dynamic11.txt