考虑断层的软化特性,建立垂直于龙门山断裂带并包含四川盆地和川西高原在内的汶川地震不稳定性地震力学模型.采用有限元方法计算得到了描述整个岩石力学系统稳定性状态和过程的平衡路径曲线.在此基础上,采用稳定性理论研究了汶川地震从孕育到发生的过程,讨论了断层倾角和断层材料参数对地震失稳的影响.数值模拟结果显示,系统只存在稳定的断层无震滑动和不稳定的地震失稳2种状态.断层倾角、初始内摩擦系数、初始黏聚力和强度曲线形状参数的增加会导致系统趋向不稳定的地震失稳状态.而强度曲线胖度参数的增加有助于系统进入稳定的断层缓慢无震滑动状态.地震失稳前,在平衡路径曲线的应力峰值点和失稳点之间,断层错动加速,应变能开始释放并且应力开始减小,是失稳的前兆.最后在失稳点发生应力突跳,地震发生,其间伴随应变能的急剧释放、应力降和断层突然错动.无论是稳定的断层无震滑动还是不稳定的地震失稳,系统重新进入应力和应变能增加状态后应力和应变能的增速由远场加载速度、岩石力学系统的结构和围岩材料属性决定,与断层软化特性参数无关.