01. 案例介绍

水坝发生地震时，由于地震惯性力突然增加，水库大坝在内部（坝体结构本身）与外部（库水的动水压力）各种作用力的作用下，大坝坝体将发生变形、坝基也将产生相对位移或不同期位移。土坝可能产生的主要损坏有：裂缝、滑坡、防浪墙与混凝土建筑物的破坏等。在水库水位较高的情况下，如果不及时抢护，短时间内就会发生溃坝，造成巨大损失。因此，掌握地震时水坝的压力分布，对于评估水坝安全性、防止塌坝等灾害的发生具有重大意义。

02. 模拟过程

常规的模拟方法是将流体压力计算入水坝的等效质量中，同时将与水坝接触的土地视为刚性物体，不考虑大地在地震时的位移。这样计算方法简单，但结果往往较实际偏差较大。在本文中，我们利用SALOME_MECA，采取边界元与有限元耦合法，考虑流体、水坝、土地之间的相互作用，模拟地震时三者的动态变化、水坝压力分布情况，并用实验取得的数据对比模拟结果，可以看出新型算法的模拟曲线很好的贴合了实际曲线，也再次检验了SALOME_MECA强大而准确的数值模拟能力。

建模几何示意图

建模网格示意图

03. 模拟结果展示

模拟结果如下图所示，蓝色为模拟曲线，红色为实际曲线。可以看到，常规算法的曲线在4到6赫兹区间时偏差很大，而新型算法的曲线始终保持与实际曲线较好的贴合度，结果更为真实、准确。

图一 常规算法与实际曲线对比

图二 新型算法与实际曲线对比

将模拟的震时水坝最大承压值（图中红色区域）与水坝极限承压值相比较，即可预测地震时水坝的安全性，并为水坝的设计、改造提供理论依据。由下图可知，使用常规算法，地震时水坝承压值最大增加1.4兆帕，而新型算法算得只会增加0.57兆帕。因此使用常规法模拟并设计水坝时，对水坝的承压能力要求过高，虽然保证水坝的安全性，但无疑增加了不必要的材料费用。使用结果更为准确的新型算法则可以在安全的前提下节省成本，提供更为合理的设计依据。

图三 常规算法 地震时水坝的压力分布

图四 有限元法 地震时水坝的压力分布

根据实验结果，SALOME_MECA能够很好地进行抗震分析。在实际使用中，它还可以处理其他多种线性、非线性问题（接触，摩擦，断裂力学等），用户可以根据自身需求编写程序，改进算法，不断开发出新型物理模块，从而更好地解决问题。

更多资讯可登录格物CAE官方网站

https://cae.yuansuan.cn/

或关注公众号【远算云学院】

远算在bilibili、头条、知乎、技术邻定期发布课程视频等内容

敬请关注