摘要
大量的震害表明，地表地形对地面运动有很强的放大和缩小作用，因而直接影响到震害的分布。对地形条件进行研究的目的有两个，一是解释地震时地面运动的放大作用和空间变化，再就是预测在将来受地震时地面运动的特性。因此地形对地面运动的影响分析是地震工程学中非常重要且具吸引力的领域。通过大量研究，局部不规则地形对地震动峰值加速度的影响我们已有相当的了解，但是地形对地震动谱特性影响的研究相对较少。本文就以探讨地形对地震动谱特性的影响规律为研究目标，通过数值计算结果分析，着重开展了以下几个方面的研究工作，并初步得到了一些有意义的结果。
1．数值计算方法的选择与分析
本文选择的数值计算方法为李小军等提出的粘弹性场地地震反应分析的一种时域高精度显式有限元-有限差分方法。第二章系统回顾了显式有限元有限差分方法的理论框架。另外，这种时域高精度的显式有限元有限差分方法的优点在于采用了李小军等的显式逐步积分方法，而逐步积分方法的基本特性除了计算格式的收敛性、稳定性、精度以外还必须考虑数值计算的高频能耗特性(dissipation)和计算起步的幅值超越性(overshoot)。因此，作为前人工作的补充，本文还在第二章中分析了中心差分格式、中心差分结合单边差分的格式以及李小军等的积分格式这三种不同显式格式的能耗特性(对高频成分的抑制和消除作用)与幅值超越性(初始计算数据的放大特性)，并对三种显式格式显示出来的这两种特性做一综合评判。研究结果表明：李小军等的显式差分格式在物理阻尼比为0到1的范围内均具有较强的能耗特性，这种特性优于另外两种显式积分格式，能有效地克服透射边界的高频失稳，另外，这种格式也具有较小的速度幅值超越性，因此，在处理复杂场地地震反应方面是一种比较理想的显式积分格式。
2．地形对地震动谱特性的影响分析
由于实际工程遇到的场地地形条件是复杂多变的，地形对输入地震波的谱特性的影响也与很多因素有关，这些因素大体上包括场地特性和输入地震波特性两个部分。场地特性包括地形的特征尺寸如坡度和高度、介质的性质等，输入地震动特性包括输入地震波的频谱特性、入射角度等。为了便于总结规律，本文控制了若干影响因素以简化问题的分析。首先研究了坡地地形在高宽比不同、体波入射角度不同时出平面运动谱特性的差异，接着进一步研究了简单平台地形及凹陷地形对入射体波谱特性的影响，给出了若干定性及定量的结果。这些研究结果表明：地形的特征尺寸和入射角度不仅对地震动的峰值有影响，而且对地震动的谱特性有较大的影响，其中入射角度的影响尤为显著。另外，研究还表明基于脉冲输入研究粘弹性地形谱特性的影响仍不失为一种有效的方法。
3．地形对输入地震动反应谱加速度峰值及特征周期影响规律的分析
若干震害调查结果表明地形影响地震动的峰值及频谱特性，即影响反应谱的形状。现有的场地反应谱对地形的考虑是以地震力的形式给出的，实际工程中碰到局部不规则地形的问题时，工程师们无法从反应谱中直接获得特征周期或峰值加速度的调整方案，本文以一个具体的平台地形为算例，通过拟合的人造地震动输入，计算获得有地形影响的场地观测点的加速度时程，再利用美国抗震设计样本规范ATC3-06建议的方法以及双参数标定反应谱的方法确定观测点加速度时程所对应的场地反应谱特征周期值，进而简要分析了地形对输入地震动反应谱峰值加速度及特征周期的影响规律，并尝试着给出了特定条件下地形对观测点反应谱的调整方案。
Many investigations of earthquake damages manifest that surfacetopography has significant amplification or reduction effects on groundmotion, and, therefore, it can influence the distribution of earthquakedamages directly. The purpose of studying on the influence of surfacetopography is to explain the reasons for the large amplifications andspatial variations of the ground motion and to predict some of the featuresof the ground motion during future earthquakes. So analysis of the effectof topography is one of the most important and attractive field inearthquake engineering. After a lot of investigation, we have know quite afew about the effect of surface topography on peak ground acceleration,but comparatively, the effect of surface topography on spectrum propertyof earthquake is still obscure for us. The objective of this thesis is todiscuss effects of topography on spectrum property of earthquakes,through numerical experiments and result analysis, the followingproblems are analyzed, and some applicable results are obtained.
1. the choice and analysis of numerical integration methods
A numerical integration method executed in this thesis is the explicitfinite difference method suggested by Li Xiaojun et al. for analyzing theearthquake response of two-dimensional visco-elastic site. The systematicreview of this method is made in the second chapter. Besides, the virtue ofthis method is using the explicit integration formula also suggested by LiXiaojun. We know for a explicit integration formula, basic characteristicsincluding astringency、stability、precision、dissipation and overshootcharacteristic. So, as complementarity of former people's work, dissipationand overshoot characteristics of three different integration formulas arediscussed in detail. The conclusions show that the explicit integrationformula suggested by Li Xiaojun and others has stronger numericaldissipation characteristic than the other 2 formulas while the real dampingvaries from 0 to 1, and the larger the real damping is, the phenomenon ismore obvious, this kind of formula can depress or eliminate the high-frequency instability induced by Local Transmitting Boundary, and,under the condition of the same precision, this kind of formula has lesservelocity overshoot characteristic, so it is an ideal integration formula forsolving the earthquake response of complex site.
2. effect analysis of topography on spectrum property of earthquake
Surface topography is complex and various, actually, effect oftopography on spectrum property refers to many factors, such ascharacteristic dimensions of topography (including height and slope angleet al.)、characters of soil or rock、spectrum properties of input wavelets、angle of incidence et al. These factors can be classified into two sorts,characters of site and characters of earthquake. In order to summarizeconveniently, firstly, we discussed the difference of spectrum property ofanti-plane movement under the conditions of different topographic factorsand direction of the incident waves, then research expands to mesa andvalley. The conclusions show that the presence of topographic factors andthe direction of the incident waves not only significantly affect the peakground acceleration of site, but also affect the spectrum property of theincoming waves, and the effect of direction of the incident waves isespecially remarkable. The study also find that it is reasonable to analysespectrum properties with the input of impulse wavelets in the finiteelement simulation of wave motion.
3. the effects analysis of topography on peak-acceleration andcharacteristic periods
Investigations after earthquakes show that surface topography affectsthe peak-accleration and spectrum property of earthquake, so, it affectsthe shape of design spectra. Considerations of topographic effects in theseismic design code for buildings are suggested in form of earthquakeforce, while it is hard for engineers to get adjusting project throughseismic design code if they have to face topographic problems. In thisthesis, a model of mesa is established, a lot of artificial ground motions aregenerated and have been used as input movements, time histories of observation points are gained, finally, two methods (one is suggested byACT3-06 and another is suggested by Liao Zhenpeng and Li Dahua)areused to get the characteristic periods of sites. With these results, theeffects analysis of topography on peak-acceleration and characteristicperiods is executed and corresponding adjusting project for design spectrais suggested.
引文
陈丙午，1982，不规则地形对地震动及震害的影响，地震工程与工程振动，2(1)，12～20。
程民宪，1988，负屈服刚度条件下数值积分的收敛性和稳定性，力学学报，2(1)，31～40。
窦立军、杨柏坡，2001，场地条件对传递函数的影响，建筑科学，17(3)，40～46。
杜修力、熊建国等，1993，平面SH波散射问题的边界积分方程分析法，中国地震，03，30～40。
杜修力、王进延，2000，阻尼弹性结构动力计算的显式差分法，工程力学，17(5)，36～43。
耿淑伟，2005，抗震设计规范中地震作用的规定，哈尔滨：中国地震局工程力学研究所博士学位论文。
郭自强，1982，固体中的波，北京：地震出版社。
戈革，1980，地震波动力学基础，石油工业出版社。
胡聿贤，1988，地震工程学，北京：地震出版社。
景立平，2004，波动有限元方程显式逐步积分格式稳定性分析，地震工程与工程振动，24(5)，20～26。
刘洪兵等，1999，地震中地形放大效应的观测和研究进展，世界地震工程，15(3)，20～25。
刘晶波，1989，波动有限元模拟及复杂场地对地震动的影响，哈尔滨：中国地震局工程力学研究所博士学位论文。
刘晶波、廖振鹏，1992，有限元离散模型中的出平面波动，力学学报．，24(2)，207～215。
刘晶波、廖振鹏，1989，离散网格中的弹性波动(2)，地震工程与工程振动，9(2)，1～11。
刘晶波、廖振鹏，1990，离散网格中的弹性波动(3)，地震工程与工程振动，10(2)，1～10。
刘晶波，1996，局部不规则地形对地震地面运动的影响，地震学报，18(2)，239～245。
刘曾武等，1982，地震波在工程中的应用，地震出版社
梁建文、张郁山、顾晓鲁等，2003，圆弧形层状凹陷地形对平面SH波的散射，振动工程学报，16(2)：158～165。
李大华，1992，振动反应数值分析的连锁公式法，地震学报，14(2)，236～242。
李大华，1989，工程抗震标准设计反应谱的统计误差，工程抗震，4，35～38。
李小军，1996，地震工程中动力方程求解的逐步积分方法，工程力学，13(2)，110～118。
李小军，1993，非线性场地地震反应分析方法的研究，哈尔滨：中国地震局工程力学研究所博士学位论文。
李小军、刘爱文，2000，动力方程求解的显式积分格式及其稳定性与适用性，世界地震工程，16(2)，8～12。
李小军、廖振鹏，1993，非线性结构动力方程求解的显式差分格式的特性分析，工程力学，10(3)，141～148。
李小军、廖振鹏、关慧敏，1995，粘弹性场地地形对地震动影响分析的显式有限元有限差分方法，地震学报，17(3)，362～369。
李小军、廖振鹏、杜修力，1992，有阻尼体系动力问题的一种显式差公解法，地震工程与工程振动，12(4)，74～80。
李方杰，2005，局部不规则场地对Rayleigh波的散射[D]，天津：天津大学
廖振鹏，2002，工程波动理论导论(第二版)，北京：科技出版社。
廖振鹏，1984，近场波动的有限元模拟，地震工程与工程振动，4(2)，1～14。
廖振鹏、刘晶波，1986，离散网格中的弹性波动(1)，地震工程与工程振动，6(2)，1～16。
梁建文、张郁山等，2003，圆弧形层状凹陷地形对平面SH波的散射，16(2)，158～165
梁建文、严林隽，2001，圆弧形层状沉积谷地对入射平面SV波散射解析解，24(2)，235～243
彭承光、李运贵，2004，《场地地震效应工程勘察基础》，北京：地震出版社，189～217。
彭国伦，2002，FORTRAN95程序设计，北京：中国电力出版社。
孙焕纯，1981，非线性结构动力分析的改进的θ法，力学学报，02，31～40。
唐晖，2005，显式积分格式及其在波动有限元模拟中的应用研究，北京：中国地震局地球物理研究所硕士学位论文。
王进廷、杜修力、赵成刚等，2003，分层介质中波动传播分析的显式有限元法，岩石力学与工程学报，22(1)，75～80。
徐士良，1982，FORTRAN常用算法，北京：清华大学出版社。
杨柏坡、陈庆彬，1992，显式中心差分有限单元法在复杂场地地震反应分析中的应用，地震研究，15(1)，70～78。
杨柏坡、陈庆彬、袁一凡，1991，地震小区划中复杂场地影响的修正方法，地震工程与工程振动，11(4)，19～27。
杨柏坡、寞立军，2000，地震安全性评价和高层建筑的地震动输入，地震工程与工程振动，20(1)，35～41。
杨柏坡、杨笑梅，1997，复杂场地上结构地震反应的研究，地震工程与工程振动，17(2)，10～15。
杨柏坡、杨笑梅，1996，二维不均匀场地地震地而运动的估算，华南地震，16(2)，1～9。
袁晓铭、廖振鹏，1993，圆弧形凹陷地形对平面SH波散射问题的级数解答，地震工程与工程振动，13(2)，1～11。
袁晓铭、廖振鹏，1995，圆弧型沉积盆地对平面SH波的散射，华南地震，15(2)，1～8。
朱镜清，1983，论结构动力分析中的数值稳定性，力学学报，1(4)，388～395。
赵成刚、李伟华等，2005，考虑耦合质量影响的饱和多孔介质动力响应分析的显式有限元法，计算力学学报，22(5)，555～561。
赵成刚、王进延、史培新等，2001，流体饱和两相多孔介质动力反应分析的显式有限元法，岩土工程学报，23(2)，178～182。
周根寿、王志良、韩清宇，1984，场地土中软弱夹层对地面地震反应的影响，《地基与工业建筑抗震》，地震出版社，20～30。
Assimaki, D., Gazetas, G., 2004, Soil and topographic amplification on canyon banks and the 1999 Athens Earthquake, JEE, 8 (1), 1～43。
Boore, D.M., 1973, The Effect of Simple Topography on Seismic Waves; Implicate on for the Aeccleration Recorded at Pacoima Dam, San Fernando Valley California, BSSA, 63 (5), 1603～1609。
Bouchon, M., 1973, Effect of topography on surface motion, Bulletin of the Seismological Society of America, 63, 615～632。
De Barros,F.C.P. and Luco,J.E., 1995, Amplication of obliquely incident waves by a cylindrical valley embedded in a layered half-space, Soil Dynamics and Earthquake Enginceling, 14, 163～174。
Geli, L., Bard, P.Y., Jullien, B., 1998, The effect of topography on earthquake ground motion: A Review and New Results, BSSA, 78 (1), 42～63。
Hilber,H.M, Hughes, T.T.R, 1997, Numerical Dissipation for Time Integration Algorithms in Structure Dynamics, Earthquake Engineering and Structure Dynamics, 5 (3), 283～292。
Lee, V.W., 1982, A note on scattering of elastic plane waves by a hemispherical canyon, Soil Dynalnics and Earthquake Engineering, 1,122～129。
Lawrence, L.D. and Lewis, R.W., 1973, Observed effects of topography on ground motion, Bulletin of the Seismological Society of America, 63, 283～298。
Liu, D.K and Han, F., 1990, Scattering of plane SH-wave by a cylindrical canyon of arbitrary shape in anisotopic media, ACTA Mechanica Sinica, 6 (3), 256～266。
Sanchez-Sesma, F.J., 1990, Elementary solutions for the response of a wedge-shaped mediuln to incident SH and SV waves, BSSA., 80 (3), 737～742。
Trifunac, M.D., 1971, Surface motion of a semi-cylindrical valley for incident plane SH waves, Bulletin of the Seismological Society of America, 61, 1755～1770。
Todorovska M, Lee V W.Surface, 1991, motion of shallow circular alluvial valleys for incident plane SH waves: analytical solution. Soil Dynamic and Earthquake Engineering, 4:192～200。
Wong, H.L. and Trifunac, M.D., 1974, Surface motion of a semi-elliptical alluvial valley for incident plane SH waves, BSSA, 64, 1389～1408。
Wong, H.L. and Jennings P. C., 1975, Effects of canyon topography on strong ground motion, BSSA, 65, 1239～1257。
William, Z.S., 2004, An Exact Solution for Effects of Topography on Free Rayleigh Waves, Bulletin of the Seismological Society of America, 94, 1706～1727。