稀疏存储的显式有限元三角网格地震波数值模拟及其犘犕犔吸收.pdfVIP

第 卷 第 期 地 球 物 理 学 报 5 6 9 , V o l . 5 6 N o . 9 年 月 2 0 1 3 9 , C H I N E S E J O U R N A L O F G E O P H Y S I C S S e p . 2 0 1 3 刘 有 山 , 滕 吉 文 , 刘 少 林 等 稀 疏 存 储 的 显 式 有 限 元 三 角 网 格 地 震 波 数 值 模 拟 及 其 吸 收 边 界 条 件 地 球 物 理 学 报 , , . P M L . 2 0 1 3 ( ) : , : / 5 6 9 3 0 8 5 ? 3 0 9 9 d o i 1 0 . 6 0 3 8 c j g 2 0 1 3 0 9 2 1 . , , , L i u Y S T e n g J W L i u S L e t a l . E x p l i c i t f i n i t e e l e m e n t m e t h o d w i t h t r i a n g l e m e s h e s s t o r e d b y s p a r s e f o r m a t a n d i t s ( ) , , ( ) : , : e r f e c t l m a t c h e d l a e r s a b s o r b i n b o u n d a r c o n d i t i o n . ? ? ? ? ? ? ? ? . ? ? ? ? ? . i n C h i n e s e 2 0 1 3 5 6 9 3 0 8 5 ? 3 0 9 9 d o i 1 0 . p y y g y ? ? / 6 0 3 8 c j g 2 0 1 3 0 9 2 1 . 稀 疏 存 储 的 显 式 有 限 元 三 角 网 格 地 震 波 数 值 模 拟 及 其 吸 收 边 界 条 件 ? ? ? , 1 2 1 1 1 刘 有 山 , 滕 吉 文 , 刘 少 林 , 徐 涛 中 国 科 学 院 地 质 与 地 球 物 理 研 究 所 , 北 京 1 1 0 0 0 2 9 中 国 科 学 院 大 学 , 北 京 2 1 0 0 0 4 9 摘 要 有 限 元 法 是 复 杂 介 质 地 震 模 拟 的 有 力 工 具 , 它 能 比 较 客 观 地 反 映 地 震 波 的 传 播 , 比 较 细 致 地 再 现 地 震 图 像 但 是 , 为 了 获 得 较 精 确 的 结 果 , 有 限 元 法 模 拟 地 震 波 的 传 播 需 要 的 网 格 点 数 多 , 具 有 计 算 量 大 和 消 耗 内 存 多 的 . 缺 点 针 对 上 述 缺 点 , 本 文 对 刚 度 矩 阵 采 用 压 缩 存 储 行 ( ) 格 式 , 以 减 少 计 算 量 并 节 省 内 存 ; 采 用 集 中 质 量 矩 阵 . C S R 得 到 对 角 的 质 量 矩 阵 以 提 高 有 限 元 法 ( 显 式 有 限 元 ) 的 计 算 效 率 ; 时 间 离 散 采 用 保 能 量 的 算 法 以 提 高 有 N e w m a r k 限 元 法 的 计 算 精 度 ; 采 用 变 分 形 式 ( 弱 形 式 ) 的 吸 收 边 界 条 件 对 人 工 截 断 边 界 进 行 处 理 通 过 与 高 精 度 的 数 值 P M L . 方 法 ? ? ? 谱 元 法 的 数 值 试 验 的 对 比 表 明 , 上 述 方 法 的 引 入 可 使 有 限 元 法 在 计 算 精 度 和 计 算 效 率 方 面 均 可 取 得 比 较 显 著 的 改 进 为 了 获 得 相 当 的 计 算 精 度 , 相 比 于 阶 谱 元 法 , 显 式 有 限 元 法 需 要 更 精 细 的 网 格 然 而 , 显 式 有 限 元 法 . 7 . 的 计 算 速 度 比 前 者 快 近 倍 , 而 内 存 需 求 仅 为 谱 元 法 的 / / 2 1 4 ~ 1 6 . 关 键 词 显 式 有 限 元 法 , 压 缩 存 储 行 ( ) 格 式 , 集 中 质 量 矩 阵 , 算 法 , 吸 收 边 界 条 件 , 谱 元 法 C S R N e w m a r k P M L : / 中 图 分 类 号 收 稿 日 期 , 收 修 定 稿 d o i 1 0 . 6 0 3 8 c j g 2 0 1 3 0 9 2 1 P 6 3 1 2 0 1 2 ? 1 1 ? 0 6 2 0 1 3 ? 0 6 ? 2 5 ? ? ? ? ? ? ? ? ? ?