flac3d边坡地震模型。 采用自由场边界，使用瑞利阻尼，输入地震波采用阪神地震波。 鲁甸波。 。

最近在折腾FLAC3D边坡地震模型，发现自由场边界配瑞利阻尼这个组合确实有意思。特别是输入阪神地震波的时候，边坡动力响应特别能反映真实场景。今儿就带大伙儿看看这模型怎么搭，中间顺带聊聊代码里藏着的小心思。

先说边界条件这块，自由场边界真是个好东西。它能模拟无限域的地震动传递，比固定边界真实多了。代码里用这个：

boundary ff

看起来平平无奇对吧？其实背后是FLAC3D用等效单元模拟远场效果。特别要注意的是模型底部必须设置成黏性边界，否则地震波会在模型底部反射搞事情：

boundary yviscous range pos-z -10 0

这里-z方向是模型底部，yviscous参数就是黏性边界的开关。

瑞利阻尼的设置就更有讲究了，直接上硬核参数：

rayleigh 0.02 0.001 frequency 2.5

第一个0.02是质量阻尼系数，后面0.001是刚度阻尼系数。重点是这个2.5Hz的中心频率，得根据边坡前几阶固有频率来定。建议先用modal分析摸个底，不然阻尼参数设了也白设。

地震波输入这块，阪神波的处理流程要特别注意基线校正。原始加速度积分后位移漂到姥姥家了，得这么搞：

table 1 import 'kobe.acc' factor 9.81
wave generate table 1 type acceleration damp 0.02
wave apply x range union name 'Base'

factor 9.81是把g单位转换成m/s²，damp参数用来消除高频震荡。这里有个坑——FLAC3D的wave apply默认是速度边界，要是在动力分析设置里没选对类型，边坡直接给你表演原地升天。

flac3d边坡地震模型。 采用自由场边界，使用瑞利阻尼，输入地震波采用阪神地震波。 鲁甸波。 。

说到结果处理，边坡关键点的时程曲线必须存下来：

history interval=50
hist gp zdisp 0,0,0
hist gp xvel 10,0,5

interval=50这个采样间隔别瞎改，太密了数据量爆炸，太疏了错过峰值。建议先试算个5秒看看数据量，再调整采样频率。

最后来个暴力提醒：动力计算时把

model large-strain on

关了！大变形和动力分析同时开，轻则计算结果漂移，重则数值爆炸。别问我是怎么知道的，说多了都是泪。

模型跑起来之后，用这个命令看边坡塑性区发展：

plot block state plane strain on

要是看到剪切带像树枝分叉一样延伸，恭喜你——边坡确实在地震作用下发生渐进破坏。这时候记得检查节点速度时程，阪神波的脉冲特性会让边坡出现明显的"回弹-滑动"交替现象。

总之这套组合拳打下来，边坡在地震中的动力响应基本能摸清。不过要真做论文的话，记得拿实际地质参数校核，别直接用默认参数——毕竟现实中的边坡可比模型调皮多了。