地震波动数值模拟及稳定性研究

论文摘要

本文针对显式有限元.有限差分法进行了系统的研究,包括理论分析,算法改进以及实际应用。通过理论推导定量地研究了采用这种数值模拟方法对波动传播规律的影响,以及其在结合局部透射人工边界进行计算时的稳定条件。分别就采用该方法结合局部透射人工边界求解波源问题和散射问题的情况分别提出了改进计算的措施,并通过数值试验验证其有效性。提出了一种有限元和显式有限元.有限差分法耦合方法,并采用其对可控震源激震过程进行了数值模拟。利用在汶川地震中获得的数据,进行场地效应分析。主要工作如下：1.根据由李小军等人提出的一种求解动力方程的显式积分格式应用于显式有限元-有限差分法进行波动数值模拟时的递推形式,给出了波动传递函数的形式；以一维均匀离散体系为模型,推导了该情况下数值波动的传递函数,并进行数值试验,验证了传递函数求解的正确性；并通过分析求解所得的传递函数,探讨了这种显式积分格式对离散网格中波动传播规律包括截止频率、频散现象和能量耗散等方面的影响。2.针对显式有限元一有限差分法结合局部透射边界进行波动数值模拟的情况,推导并计算了循环系数,给出了计算过程的稳定条件,并通过数值试验对该判定准则进行验证；在给出的循环系数的基础上,定量地分析了积分格式的能耗特性对由边界条件引入的高频失稳的抑制和消除作用；并探讨了如何更有效地利用能耗特性保证计算稳定性,即实际计算模型中时间离散步距以及空间模型大小等参数如何取值既能确保计算稳定又能最大程度减小计算量。3.分别针对采用显式有限元一有限差分法结合局部透射人工边界求解波源问题和散射问题的情况分别提出了改进措施。对于波源问题,建议针对每个边界节点可以采用不同的人工波速,该人工波速可以通过试算一定时间长度,然后通过相关处理获得。通过数值实验证明了,在波源问题的求解中采用这种方法选择人工波速可以在不增加计算失稳风险的同时提高数值计算精度。对于散射问题,则给出了一种关于波动输入方法的改进措施,即在分区分离波场中入射波不再采用解析解而是采用离散数值解代替,数值实验证明了采用该方法可以有效地抑制漂移失稳。4.给出了一种基于有限元和显式有限元.有限差分法相耦合的区域分裂法。将计算模型分为两个区域,对于广义结构及其附近区域采用有限元法模拟,而土体其它部分采用显式有限元.有限差分法,两个区域的模拟通过搭接区域实现数据的交换。分别代表波源问题和散射问题的数值试验都验证了该方法的有效性。采用提出的耦合方法针对可控震源激振过程建立的数值计算模型,并对不同的工况进行计算,研究分析了各因素包括接触土体密度、接触土体刚度、接触土层厚度、软、硬夹层以及底平板半径对可控震源能量输入大地体系的影响,为如何能高效地将激震能量输入土体提供了参考。5.利用自贡地形影响台阵在2008年汶川地震中获得的主震强震记录,采用水平/垂直谱比法（HVSR）对该山脊地形场地对地震动的放大效应进行了分析,并探讨了所获记录结果的可靠性；采用显式有限元一有限差分法结合局部透射人工边界的二维场地影响分析模型与方法,模拟了该场地的地震动响应。对模拟结果与HVSR法分析结果进行对比分析,表明数值模拟结果与观测记录分析结果具有较好的一致性,利用数值计算方法建立的简单二维分析模型基本上可以反映山脊地形对地震动的影响。

论文目录

• 摘要
• Abstract
• 目录
• 第一章 绪论
• 1.1 研究背景
• 1.2 研究现状
• 1.2.1 数值模拟方法
• 1.2.2 边界条件
• 1.2.3 波动输入或加载方式
• 1.3 存在的问题
• 1.4 论文研究内容及大纲
• 第二章 显式积分格式对波动传播特性的影响
• 2.1 引言
• 2.2 有限元-有限差分法概述
• 2.2.1 空间离散对波动传播规律影响
• 2.2.2 时间离散方法
• 2.3 显式积分格式对波动传播规律影响的研究
• 2.3.1 传递函数的推导
• 2.3.2 数值实验
• 2.3.4 传递函数分析研究
• 2.4 显式积分格式与精细时程积分格式的对比
• 2.5 小结
• 第三章 显式积分格式对局部透射边界高频失稳的抑制和消除作用
• 3.1 引言
• 3.2 全局人工边界条件
• 3.2.1 一致边界
• 3.2.2 边界元法
• 3.2.3 其它一些边界条件
• 3.3 局部人工边界
• 3.3.1 旁轴近似边界
• 3.3.2 higdon边界
• 3.3.3 粘性边界
• 3.3.4 粘弹性边界
• 3.3.5 完美匹配层
• 3.3.6 透射边界
• 3.4 多次透射公式及其稳定性问题
• 3.4.1 多次透射公式及其数值实现
• 3.4.2 循环系数
• 3.4.3 数值试验
• 3.5 小结
• 第四章 显式有限元-有限差分法结合局部透射人工边界波动模拟应用改进措施
• 4.1 引言
• 4.2 二维波动模拟数值实现
• 4.2.1 二维显式有限元-有限差分法
• 4.2.2 局部透射人工边界条件
• 4.2.3 波动输入方法
• 4.3 波源问题
• 4.3.1 问题的提出和改进措施
• 4.3.2 数值实验
• 4.4 散射问题
• 4.4.1 问题的提出以及改进措施
• 4.4.2 数值实验
• 4.5 小结
• 第五章 有限元和显式有限元-有限差分法区域分裂耦合法及对可控震源激震过程数值模拟
• 5.1 引言
• 5.2 土-结构相互作用的混合方法
• 5.2.1 有限元与无穷元耦合法
• 5.2.2 有限元与边界元耦合法
• 5.2.3 边界元与无限边界元耦合法
• 5.2.4 有限元与有限差分耦合法
• 5.2.5 其它耦合法
• 5.3 有限元与显式有限元-有限差分耦合法
• 5.4 数值验证及讨论分析
• 5.4.1 算例一
• 5.4.2 算例二
• 5.5 可控震源的原理及研究现状
• 5.6 数值试验
• 5.7 小结
• 第六章 自贡西山公园山脊地形场地效应分析
• 6.1 引言
• 6.2 场地效应分析方法
• 6.2.1 经验法
• 6.2.2 理论法
• 6.3 自贡地形影响台阵及数据
• 6.4 数据分析及数值模拟计算结果
• 6.5 小结
• 第七章 总结与展望
• 7.1 本文主要成果及结论
• 7.2 进一步研究工作展望
• 参考文献
• 作者简介
• 攻读博士期间主要参与的课题
• 攻读博士期间发表文章
• 致谢