非均匀椭球地球地震和地表载荷库伦应力计算方法

摘要

本发明公开了一种非均匀椭球地球地震和地表载荷库伦应力计算方法，利用将地震剪切位错和地表载荷分别等效为有限元计算中的体力项和面力项，再通过局部网格自适应加密技术保证计算精度，此方法可计算地震包括震间、同震和震后，以及地表载荷包括地形剥蚀、冰川消融和水库抽蓄水的加卸载引起的库伦应力变化；具体步骤是：地震剪切位错与地表荷载的等效替代、网格自适应加密、粘弹性maxwell体积分型本构递推实现和横向非均匀椭球形地球的实现，最后通过程序实现自动计算。本发明的有益效果是：其计算更为全面和准确的体现了地球浅表圈层地震和火山库伦应力变化的全过程。