地震波斜入射情形下水平成层半空间自由场地时域计算.pdfVIP

第33卷 第3期 吉林大学学报(地球科学版) V01_33No．3 OF UNIVERSITY(EARTHSCIENCEEDITION) 2003年7月 JOURNALJII．IN Ju1．2003 文章编号：1671—5888(2003)03—0372—05 地震波斜入射情形下水平成层 半空间自由场的时域计算 李山有1，王学良2，周正华1 】．中国地震局工程力学研究所，黑龙江哈尔滨】50080；2，云南省楚雄师范学院地理系，Z-南楚雄5750001 摘要：在应用有限元模拟地震工程中的波动散射问题时．为了计算边界节点的运动和实现波动能 量的输入，必须首先确定边界节点的内行波场即背景水平成层半空间的自由场。利用显式有限差分 法的内节点位移计算公式和基于水平成层介质波动传播特点建立的相邻节点问自由场运动的关系 式．给出了计算水平成层半空间自由场的时域方程(咀入射倒边界节点为例)。本方涛是显式有限元 意义上的数值精确解，可应用于近场波动数值模拟中内行波场与散射波场的时域同步模拟。 关键词：水平成层半空间I地震波倾斜八射；自由场；时域算法 中图分类号：P315．9 文献标识码：A 条件时需从总波场中扣除沿边界传播或向计算区内 0 引言 传播的非外行波场以得到外行波。对于时空解耦的 在应用有限元或有限差分法处理地震工程中的 人工边界条件来说，可以针对底边界及侧边界分90 波动散射问题时，必须从半无限背景介质中切取有 定义非外行波，例如一种具有普遍意义的定义是用 限的计算区，并在计算区边界上引人台适的人工边 人射地震波作为底边界的非外行渡，侧边界背景介 质的自由场作为侧边界的非外行波“J。目前，远离 散射体的背景介质通常假设为水平成层，而入射地 震波也假定为垂直向上传播的剪切波或压缩波，此 时自由场的求取是一个简单的一维问题，可咀方便 地在时域实现。而事实上，当震源距场地较近时，地 震波并不是垂直向上入射的，为了确保重大工程在 地震作用下的安全，工程师们已经认识到考虑地震 波入射角度不确定性的必要性。此时，为了实现波 动散射问题的数值模拟，通常采用“两步法”，即先在 频域内一次性计算好人射侧边界节点的自由场运动 时程，然后以数组或文件形式传送到时域的边界处 理中oo。这一实施方案不仅操作较繁，而且还要求 凰l 波动有限元模拟力学模型 计算机有更大的内存和存储空间。为了实现近场波 of Fig．1Themechanicalmodelfinitedement 动数值模拟中内行波场与散射波场的时域同步模 界条件”‘(图1)。由