时域局部透射边界的计算飘移失稳1).pdf

第 28 卷 第 5 期 V o l. 28, N o. 5 1996 年 9 月 A CTA M ECHAN ICA S IN ICA Sep t, 1996 时域局部透射边界的计算飘移失稳1) 李小军 ( 国家地震局地球物理研究所, 北京 100081) 廖振鹏 ( 国家地震局工程力学研究所, 哈尔滨 150080) 摘要　本文通过对时域局部透射边界特性及应用结果分析, 发现了多次透射边界的一种 非高频振荡失稳的失稳性, 并就此给出了人工边界的计算飘移失稳的概念文中对多次透 射边界的计算飘移失稳现象及失稳机制作了分析, 探讨了飘移失稳发生的条件, 提出了 一种简单有效的消除飘移失稳的方法 关键词　多次透射边界, 飘移失稳, 振荡失稳 引 　 言 地震及其它动力荷载下的地表地球介质及建筑物的动力问题均属于近场波动问题有 限元方法与时域逐步积分方法在近场波动问题分析中的应用导致了时域局部人工边界的出 现时域局部人工边界中, 显式时域多次透射边界与有限元及有限差分方法配合使用成为建 立近场波动问题分析的显式数值方法[ 1, 2 ] 的唯一可用的高精度人工边界方法 [3 ] 多次透射边界 是精度可控的高精度外推人工边界, 但它存在着有待进一步研究计算 [4 ] 稳定性问题廖振鹏等 其高频振荡失稳问题作了系统的研究, 给出了抑制高频失稳现象的 措施时域局部透射边界还存在另一失稳性问题, 即计算飘移失稳, 这正是本文所要论述的 问题 1　计算飘移失稳 ( ) 对于有限元离散模型 如图1所示 , 时域局部多次透射边界的表示式为 N uP + 1 = (- 1) n+ 1CN uP + 1 ( 1) R J n n , J n= 1 图1 空间离散模型中的透射边界 F ig 1 T ran sm it t ing boundary in th e d iscrete m odel in sp ace 1) ( ) 地震科学联合基金资助项 目 900002 1 , 国家地震局地球物理研究所论著: 96 0032 A 收到第一稿, 19960 129收到修改稿 628 力　　　学　　　学　　　报　　　　　　1996 年 第 28 卷 其中 P ′ 1