摘要：运用三维非线性动态有限元方法（Non-linear Dynamic FEM）仿真模拟研究了菱形模型的动态应力作用，分析了1976年唐山7．8级地震．结果表明：（1）唐山地震前增强的中强地震产生的冲击力形成的应力波反射后发生半波损失，成为动态的拉伸应力，通过减小内摩擦而引起岩石产生裂纹、局部破裂并发生串通，一方面触发地震，另一方面为孕育唐山地震积累能量，表现为正反馈过程．（2）唐山地震产生的冲击力加卸载存在时间差，造成加卸载冲击波相位不同，结果产生的动应力和围压叠加致使围压波动和下降，导致抗剪切强度降低，使初始破裂分别向北东、西南方向传播．（3）冲击力的卸载，相对于产生拉伸性的应力波，反射后发生半波损失，成为动态的压缩应力，致使围压增大并导致破裂停止；而动态应力在传播过程中衰减．因此，余震的发生是一个负反馈过程．

关键词：唐山地震 三维非线性动态有限元方法 动态应力作用 低值波动 

单位：中国地震应急搜救中心 北京100049 华中科技大学 武汉430074 中国地震局地壳应力研究所 北京100085 中国地震台网中心 北京100045