第３３卷第３期
２００３年７月
吉林大学学报（地球科学版）
ＪＯＵＲＮＡＬＯＦＪＩＩ．ＩＮＵＮＩＶＥＲＳＩＴＹ（ＥＡＲＴＨ
ＳＣＩＥＮＣＥＥＤＩＴＩＯＮ）
Ｖ０１＿３３
Ｊ
Ｎｏ．３
ｕ１．２００３
文章编号：１６７１—５８８８（２００３）０３—０３７２—０５
地震波斜入射情形下水平成层半空间自由场的时域计算
李山有１，王学良２，周正华１
】．中国地震局工程力学研究所，黑龙江哈尔滨】５００８０；２，云南省楚雄师范学院地理系，Ｚ－南楚雄５７５０００１
摘要：在应用有限元模拟地震工程中的波动散射问题时．为了计算边界节点的运动和实现波动能量的输入，必须首先确定边界节点的内行波场即背景水平成层半空间的自由场。利用显式有限差分
法的内节点位移计算公式和基于水平成层介质波动传播特点建立的相邻节点问自由场运动的关系
式．给出了计算水平成层半空间自由场的时域方程（咀入射倒边界节点为例）。本方涛是显式有限元意义上的数值精确解，可应用于近场波动数值模拟中内行波场与散射波场的时域同步模拟。
关键词：水平成层半空间Ｉ地震波倾斜八射；自由场；时域算法中图分类号：Ｐ３１５．９
文献标识码：Ａ
０
引言
在应用有限元或有限差分法处理地震工程中的
条件时需从总波场中扣除沿边界传播或向计算区内
传播的非外行波场以得到外行波。对于时空解耦的
人工边界条件来说，可以针对底边界及侧边界分９０
定义非外行波，例如一种具有普遍意义的定义是用
波动散射问题时，必须从半无限背景介质中切取有限的计算区，并在计算区边界上引人台适的人工边
人射地震波作为底边界的非外行渡，侧边界背景介
质的自由场作为侧边界的非外行波“Ｊ。目前，远离散射体的背景介质通常假设为水平成层，而入射地震波也假定为垂直向上传播的剪切波或压缩波，此
时自由场的求取是一个简单的一维问题，可咀方便
地在时域实现。而事实上，当震源距场地较近时，地震波并不是垂直向上入射的，为了确保重大工程在
地震作用下的安全，工程师们已经认识到考虑地震波入射角度不确定性的必要性。此时，为了实现波动散射问题的数值模拟，通常采用“两步法”，即先在
频域内一次性计算好人射侧边界节点的自由场运动时程，然后以数组或文件形式传送到时域的边界处
凰ｌ
Ｆｉｇ．１
Ｔｈｅ
波动有限元模拟力学模型
ｍｅｃｈａｎｉｃａｌｍｏｄｅｌｏｆｆｉｎｉｔｅｄｅｍｅｎｔ
理中ｏｏ。这一实施方案不仅操作较繁，而且还要求
计算机有更大的内存和存储空间。为了实现近场波动数值模拟中内行波场与散射波场的时域同步模
拟，本文基于水平成层介质在地震波斜入射情形下波动传播水平视波速不变且已知的特点，建立了人射侧边界节点自由场计算的内插公式，并将其与内
界条件”‘（图１）。由于人工边界条件只能模拟向
计算区外传播的外行波，这样在数值实现这些边界
收稿日期：２００２—０８
２８
基金项目：国家自然科学基金资助项目（５０１７８０６５）；地震科学联合基金资助课题（ｉ０１０６２）
作者简介：李山有（１９６５一），男，河北碰县人，高级工程师．在站博士后．主要从事强震动数据分析与近场波动数值模拟
研究，Ｅ—ｍａｉｌ：ｌｉ
ｓｈａｎｙｏｕ＠１６３
ｃｏｍ．