三维地震断层动力破裂的显式并行有限元解法 - 中国地震信息网.PDF

维普资讯 第23卷 第 4期 地 震 工 程 与 工 程 振 动 Vo1．23．No．4 2003年 8月 EARTHQUAKEENGINEERINGANDENGINEERINGVIBRATION Aug．，2003 文章编号 ：1000—1301(2003)04—0022—05 三维地震断层动力破裂 的显式并行有 限元解法 刘启方 ，丁海平 ，袁一凡 ，金 星 ，黄 勇 (1．中国地震局工程力学研究所 ，黑龙江 哈尔滨 150080；2．哈尔滨工业大学 黑龙江 哈尔滨 150o90) 摘要 ：本文在近场波动有限元的基础上，根据断层动力破裂的摩擦模型 ，详细地推导出了断层面节 点 在破裂的不 同时刻 的计算公式 ，提 出了一种模拟断层动力破裂的显式并行有限元方法 ，利用此方法 ， 可以模拟断层 的破裂过程 、地表破裂及由此产生 的地面运动。 关键词 ：断层动力破裂 ；摩擦模型 ；显式有 限元 ；并行计算 中图分类号 ：P315．01 文献标识码 ：A An explicitparallelfinite-elementmethodofdynamic rupturein 3D earthquakefault LIUQi—fang。DINGHai—ping一，YUANYi—fan ，JINXing，HUANGYong (1．InstituteofEngineeringMechanics，ChinaSeismologicalBureau，Harbin150080，China； 2．HarbinInstituteofTechnology，Harbin150090，China) Abstract：Basedonthefinite—elementmethodofReal一fieldwavemotion，wepresentan explicitparallelfinite—ele— mentmethod tosimulatethedynamicruptureof3D earthquakefaultaccordingtothefrictionalmode1．Theformula— tionofdifferenttypeofthenodeinthefaulthasbeenobtained．Theurptureandslipprocessofhtefault，htestrong groundmotionandthesurfaceurpturecan besimulatedbythismethod ． Keywords：dynamicurptureoffault；frictionalmodel；explicitfinite—element；praallelcompute 1 前言 近年来的几次大地震 (1994，Northridge；1995，Kobe；1999，Chi—Chi)表明，断层 的破裂过程 以及地震波与 断层的动力相互作用对近场地面运动有重要影响，许多学者利用数值方法对断层的动力破裂过程及 由此产 生的地面运动作了研究 。通过对断层动力破裂过程的研究可以得到断层面的位错过程、断层破裂引起 的地面运动以及断层能否破裂至地表和地表破裂的尺度 ，同时可以进一步加深对震源物理 的认识。断层的 动力破裂的数值模拟方法主要有 3种 ，即：有限差分法 (FDM)¨ 、边界积分方程法 (BIEM) -5]和有限元方 法 (FEM) 。有限差分法不适宜处理复杂的边界 问题，边界积分方程法是建立在弹性互易定理的基础上 的，所 以不适宜处理非线性问题 ，有限元方法则既能处理复杂的边界又能用于计算非线性问题 ，它是处理断 层动力破裂的最有效的方法 。Aagaard 利用三维有限元方法成功地模拟 了不同初始条件下断层的动力破 裂过程。由于断层的动力破裂是一个三