论文部分内容阅读
针对有限元法(FEM)处理不连续位移场所存在的困难,结合离散元法(DEM)的优点,提出了一种DEM/FEM耦合计算模型。基于罚函数法,通过建立虚功方程及变分原理,得到了界面力的计算公式,实现了两种计算域之间的平滑过渡。在小变形区域用FEM模拟,在弹塑性大变形区域采用DEM进行模拟。将耦合模型应用于单层球面网壳倒塌数值模拟分析,并与课题组前期开展的球壳振动台试验模型进行对比。计算结果表明,网壳第一、第二圈斜杆首先进入塑性并发生弯曲,逐级提高地震波峰值加速度,进入塑性杆件不断增多直至结构最终坍塌,与试验过程及试验现象均较为吻合,说明了本文DEM/FEM耦合计算模型的有效性。