岩石膨胀的三维本构关系和1976年唐山地震前垂直形变的有限元模拟

本文研究了岩石在应力作用下的膨胀问题。依据实验资料和参考Nur的岩石膨胀二维本构关系,提出了岩石介质膨胀的三维非线性本构关系的一般表达式。並将这组非线性关系式推导成增量形式的有限元矩阵表达式。 考虑到导出的膨胀的本构关系矩阵是非对称矩阵,因此,本文编制了非对称全带宽存储及直接求解的三维非线性有限元程序。为了模拟实际断层的情况,考虑断层的蠕滑,在程序中还专门设计一种非线性空间断层单元,以便在研究真实地震问题或有关岩石膨胀的问题中作数值模拟。 用上述带有弹性-膨胀矩阵三维非线性有限单元法,假设不同模型,对唐山大震前出现的明显隆起进行三维数值模拟研究,並结合小震分布等震前现象,对1976年7月28日唐山地震的孕育过程提出看法。