颗粒离散单元法地震动力时程计算黏性人工边界及其应用

采用颗粒离散单元法(PFC)进行地震动力时程计算时,人工截断边界上需设置吸收边界条件,以防止波的反射。鉴于PFC数值计算模型的人工边界上颗粒单元半径大小不一、边界面凸凹不平,在连续介质的黏性吸收边界条件方程基础之上,推导出适用于离散介质的等效方程。在离散介质的等效方程中引入微调系数,建立确定其最优值的方法,以实现对波的最佳吸收。采用合成岩体方法(SRM)建立顺层岩质边坡和逆层岩质边坡结构面网络模型。基于提出的方法,在边坡颗粒流数值模型的人工边界上设置黏性吸收边界,对含有正交次级节理的顺层岩质边坡和逆层岩质边坡的地震破坏过程进行模拟,对其地震破坏机制进行系统分析,相关研究结论对边坡地震破坏机制研究具有重要意义。

ABAQUS中动力问题边界条件的选取

数值模拟是研究半无限域中结构地震反应的重要手段,而边界条件的选取则直接影响到数值模拟结果的可靠性。本文应用有限元软件ABAQUS,分别采用固定边界条件、黏性边界条件和远置边界条件对简单弹性半无限介质的波源问题和散射问题进行了数值模拟分析,结果显示:使用固定边界时,波动不能穿过边界面进入半无限介质,这与波动在半无限介质中传播的规律不符。通过与远置边界对比可知:使用黏性边界条件可以很好地模拟地震波穿过人工边界向半无限介质传播的真实情况,可以得到与使用远置边界相同精度的结果。

基于FLAC3D的斜入射地震波作用的数值模拟方法研究

地震动输入是进行土层场地或土-结构相互作用体系地震反应分析的重要一步。通过波场分解,将人工边界上的波动场分解为自由波场和散射波场,其中人工边界上的自由波场根据波在自由表面的反射理论求得,散射波场由黏性人工边界直接吸收,进而在人工边界上施加等效荷载实现斜入射地震波作用的输入。借助MATLAB辅助计算程序和FLAC3D的前后处理功能,进行不同角度(0°和30°)入射条件下SV波作用的时域数值模拟计算,数值解和理论解吻合较好。算例表明,采用黏性人工边界和等效荷载输入方法具有较高的计算精度,建立的数值模拟方法为土层场地和土-结构动力相互作用在地震波斜入射下的反应分析提供了简单易行的计算途径。