针对这一问题,岩土力学与工程国家重点实验室计算岩石力学组首次采用虚单元法对于含节理大坝中的接触问题进行了模拟。与传统的数值方法相比,虚单元法不需要对形函数进行显示的表达,计算高效及处理网格灵活,构造精度可达任意阶数,能克服复杂性的块石几何形状带来的畸形网格问题,适用于任意复杂的单元形状。基于虚单元法的优点,推导了虚单元法在不同接触状态下的刚度矩阵。
为了验证虚单元法求解接触问题的准确性,首先算例验证了虚单元法对于网格扭曲度敏感性,随后通过一系列的数值算例验证了虚单元法求解接触问题的可行性(图1),最后通过含节理大坝的模拟,说明了节理对于大坝稳定性有重要的影响(图2)。
相关研究成果以“The Virtual Element Method for the Dam Foundation With Joint”为题,发表在Frontiers in Earth Science。研究工作获得了国家自然科学基金面上项目(11972043)的资助。
论文链接:https://www.frontiersin.org/articles/10.3389/feart.2022.875561/full
图1 均匀压缩下的水平接触界面模拟(大主应力云图)
图2 带节理的混泥土大坝 (x方向位移云图)