离散单元法PFC中动力人工边界的设置方法

数值计算是研究岩土体地震动力响应规律及灾变机制的有效方法，计算工作的关键之一就是设置合理动力人工边界条件。动力人工边界理论上应当实现对原介质应力应变场的精确模拟，保证波在人工边界处的传播特性与原介质一致，使外行散射波通过人工边界时无反射效应，发生完全的透射或被人工边界完全吸收。

颗粒离散单元法（DEM）由于对每个颗粒单元应用牛顿第二定律，在接触位置应用力-位移定律并不断更新，颗粒单元运动不受变形量的限制，在研究岩土体失稳大变形问题时优越性显著。

为了能够采用颗粒离散单元法（DEM）正确有效地模拟斜坡等岩土体震裂破坏及运动大变形的全过程，建立合理的DEM动力人工边界显得尤为关键。

本案例计算了在离散元PFC中进行人工边界的两种方法：（1）高阻尼法，（2）黏性人工边界。同时还对固定人工边界的情况进行了对比

其中高阻尼法的原理是为边界处的颗粒设置局部阻尼系数为1.0，以避免边界处的波动反射。

黏性人工边界的方法则是基于人工边界的理论实现的，（参考了这篇论文，表以敬意。周兴涛,盛谦,崔臻,冷先伦,付晓东,马亚丽娜.颗粒离散单元法动力人工边界设置方法[J].岩土力学,2018,39(07):2671-2680+2690.）

颗粒离散单元数值模型的黏性边界条件设置条件如下： 

本算例实现了在颗粒流PFC对黏性人工边界的模拟。本算例进行了一维波的测试，并模拟对比了两种人工边界的模拟方法。

首先建立模型，在边界墙体的伺服功能下平衡模型：

删除边界墙体，对右侧边界附近的颗粒运动进行约束模拟边界，对于高阻尼的方法模拟是为右侧边界处的颗粒设置局部阻尼系数为1.0。黏性人工边界的方法则是对边界颗粒施加荷载，吸收掉入射的波动能量，以模拟无限介质。为左侧边界附近的颗粒施加脉冲波。

固定边界模型中部处球颗粒的水平速度时程曲线如下：

高阻尼边界模型中部处球颗粒的水平速度时程曲线如下：

黏性边界模型中部处球颗粒的水平速度时程曲线如下：

付费内容为算例的完整代码。