利用谱元法计算弹性波场的若干理论问题 |
| |
作者姓名: | 王秀明 G.SERIANI 林伟军 |
| |
作者单位: | 1. 中国科学院声学研究所,北京,100080;CSIRO Petroleum, ARRC, POBOX 1130, Technology Park, Bentley, WA 6102, Australia 2. Instituto Nazonnale di Oceangrafiae di Geofisica Sperimentale Borga Grotta Gigante 42/c-Sgonico, I-34016, Trieste, Italy 3. 中国科学院声学研究所,北京,100080 |
| |
基金项目: | 国家自然科学基金;UNESCO国际理论物理中心基金、澳大利亚教育和科技部国际科学合作计划以及中国科学院百人计划联合基金 |
| |
摘 要: | 近来在解决弹性动力学问题中新发展起来的谱元法是一种基于组成整个研究区域的单元Galerkin算法, 并利用弹性波动力学方程的弱形式进行数值求解的方法. 针对求解弹性波动方程的谱元法中的若干数学问题以及相关的算法原理进行了系统研究, 介绍了Legendre 和Chebyshev多项式构造的基函数及其用在参考元上数值积分的Gauss-Lobbatto配置点, 推导了在利用Legendre和Chebyshev多项式展开时的单元积分的具体表达式. 在声阻抗差别很大的非均匀介质中, 通过引入空间域中在预先条件下的共轭梯度的元到元算法和时间域中时间积分的交错网格的预期/多次校正算法, 谱元法不需要形成有限元中的全局矩阵和有效载荷矢量, 能够在很大程度上提高计算精度和计算效率. 另外, 在某些情况下, 如果再利用单元积分的解析式, 也不需要形成单元矩阵, 由于采用解析式, 计算中可以节约大量内存而不会损失太多计算效率. 元到元算法中使用了最优的张量乘积技术, 使得该方法比有限元法在内存需求量和计算时间等方面更为有效.
|
关 键 词: | 弹性波模拟 谱元 元到元算法 弱形式 |
收稿时间: | 2006-01-10 |
修稿时间: | 2006-01-102006-08-28 |
本文献已被 CNKI 维普 万方数据 等数据库收录! |
| 点击此处可从《中国科学(G辑)》浏览原始摘要信息 |
|
点击此处可从《中国科学(G辑)》下载全文 |
|