多次透射公式的一种稳定实现措施

基于对双曲偏微分议程数值解稳定性的GKS（Gustafsson,Kreiss and Sundstrom）准则的物理解释，本文提出了在波动模拟中多次透射公式飘移失稳的一个简单措施，并用波源问题和散射问题的数值试验检验了建议的稳定措施的有效性。     

多次透射公式实用形式稳定性分析

本文以二维出平面波动方程为例,利用数值方法研究了多次透射公式的实用形式的稳定问题,结果表明：多次透射公式的实用形式与波动方程的集中质量显式有限元结合,不但存在零频漂移,而且同样存在高频振荡失稳,而与中心差分方法结合仅存在零频漂移失稳。     

求解散射问题中消除多次透射边界飘移失稳措施的效果分析

针对应用多次透射人工边界公式求解散射问题时可能出现的飘移失稳现象,对比分析了实时降阶法、修正算子法、附加黏弹器法和改进输入法的消飘效果和对计算精度的影响,结果表明:4种方法均可有效消除或抑制比较缓慢的飘移失稳现象;飘移趋势强烈时,修正算子法和附加黏弹器法可能无法有效地抑制失稳;当阻尼较小时,采用实时降阶法可能引入额外的低频误差;在能够确定适当参数值的前提下,修正算子法和附加黏弹器法能够提高低频计算精度;改进输入法则无需人为设置参数,能够在抑制飘移失稳的同时提高计算精度。     

多次透射边界模拟成层介质中瞬态波动的精度分析

本文将对多次透射边界在瞬态波动模拟中的精度分析推广到介质中存在多个物理波速的情形，阐明了它模拟成层介质中瞬态波动的适用性，指出多次透射边界的误差随着透射阶数的增加而迅速减小，二阶透射边界已具有足够高的模拟精度。     

多次透射公式的一种高频失稳机制

本文利用一维双曲型偏微分方程组初边值问题数值稳定性的GKS定理的物理解释和推广，分析了多次透射公式在数值实现中可能出现的一种高频失稳机制，即在多维离散网格，沿某一空间方向的外行简谐波，由于与其他空间方向节点运动的耦合效应可以使能量传播方向反向。文中通过数值试验对这一失稳机制作了初步验证。     

散射问题中消除多次透射边界高频振荡失稳措施比较分析

多次透射边界是一种应用广泛的人工边界,类似于其它局部人工边界,它也存在计算稳定性问题。本文分析多次透射边界高频振荡失稳的原因,通过数值试验对已有的消除高频失稳措施进行比较分析,并给出在散射问题中已有消除高频失稳措施的适用性。研究表明,对于完全弹性问题,建议用滤波方法消除高频振荡失稳,但滤波系数的取值需要依据一定的经验;对于黏弹性问题,建议适当控制阻尼取值,利用差分格式本身的高能耗特性来抑制高频振荡失稳,该方法不存在参数选取的问题,更便于应用实施。     

散射问题中消除多次透射边界高频振荡失稳措施比较分析

多次透射边界是一种应用广泛的人工边界,类似于其它局部人工边界,它也存在计算稳定性问题。本文分析多次透射边界高频振荡失稳的原因,通过数值试验对已有的消除高频失稳措施进行比较分析,并给出在散射问题中已有消除高频失稳措施的适用性。研究表明,对于完全弹性问题,建议用滤波方法消除高频振荡失稳,但滤波系数的取值需要依据一定的经验;对于黏弹性问题,建议适当控制阻尼取值,利用差分格式本身的高能耗特性来抑制高频振荡失稳,该方法不存在参数选取的问题,更便于应用实施。     

局部透射边界的精度

离散局部透射边界和旁轴近似的精度分析表明,它们所以缺乏透射大角度入射波的能力乃是源于采用单一的透射方向。本文提出将边界节点附近的入射波表示为一系列预定入射角的平面波的叠加以消除这一缺点。在这一思想上建立的局部透射边界能够在全部入射角范围内减小反射。通过模拟稳态SH波动简要地阐明了上述思想及其在建立高精度、稳定的局部透射边界中的意义。     

消除多次透射公式高频振荡失稳的一种措施

多次透射公式（MTF）物理概念简单,便于在计算机上实现时空解藕的高精度波动数值模拟。然而,MTF与其它局部人工边界条件类似,存在数值模拟失稳问题,如高频振荡便是可能出现的失稳现象。本文在分析MTF高频振荡失稳机理的基础上,提出了在波动有限元数值模拟中消除MTF高频振荡失稳的一种措施,即在整个有限元数值模拟区内施加与应变速率成正比的较小粘性阻尼;同时,讨论了这一稳定措施的有效性及其对数值计算精度的影响,并通过数值试验检验了这一稳定措施的可行性。结果表明,消除高频振荡失稳的措施行之有效,且只对波动有限元数值模拟中无意义的高频分量具有抑制作用,而对有意义的较低频段内的波动有限元数值模拟精度影响较小。     

局部透射边界在稳态波动有限元模拟中的应用

本文将局部透射边界和集中质量有限元技术结合起来形成一个模拟稳态近场波动的力学模型。这一模型中任一节点的运动仅与其邻近的节点耦联。这一特点导致稳态波动有限元模拟的一种十分简便和有效的算法。本文将以基底固定的弹性层内的出平面波动为例阐明这一算法,并用一种新的局部透射边界模拟行进波动。     

较复杂介质波场计算中的一种透射边界

传统的人工边界是在声波或弹性波分解为内行和外行分量的基础上设计的.对于较复杂的介质模型（如粘弹介质、弹塑介质、塑性介质模型等）,地震波不能简单地分解为各自满足其运动方程的纵波、横波等成分,且地震波在传播过程中还有衰减现象,因此传统的人工边界已不能有效地用于非弹性的较复杂介质的波场模拟.本文提出了一种适合于更一般的地震模型的波场计算的透射边界条件,并以粘弹介质的一个较简单的二维几何模型为例,通过有限元数值求解这种透射边界和刚性边界下粘弹波的激发和传播,得到地面合成地震记录,证明了这种透射边界用于较复杂介质波场计算是可行的.     

较复杂介质波场计算中的一种透射边界

传统的人工边界是在声波或弹性波分解为内行和外行分量的基础上设计的.对于较复杂的介质模型（如粘弹介质、弹塑介质、塑性介质模型等）,地震波不能简单地分解为各自满足其运动方程的纵波、横波等成分,且地震波在传播过程中还有衰减现象,因此传统的人工边界已不能有效地用于非弹性的较复杂介质的波场模拟.本文提出了一种适合于更一般的地震模型的波场计算的透射边界条件,并以粘弹介质的一个较简单的二维几何模型为例,通过有限元数值求解这种透射边界和刚性边界下粘弹波的激发和传播,得到地面合成地震记录,证明了这种透射边界用于较复杂介质波场计算是可行的.     

一种时域复阻尼本构方程

为了将复阻尼理论应用于时步积分计算,本文根据文献［１］提供的建立时域复阻尼本构的思路,首先讨论了通过组合标准线性体构成与复阻尼理论在具有工程意义的频段内等效的时域本构方程的方法,即通过适当选取本构参数使滞回阻尼比和弹性模量与复阻尼理论相应值基本保持一致。然后,用一维波动数值模拟实例,对比了本文建议的本构方程、频域复阻尼理论及粘性阻尼本构对计算结果的影响,论证了本文建议的本构方程应用于地震工程的可行     

基于混合有限元格式的完美匹配层与多次透射公式人工边界比较研究

介绍了完美匹配层(PML)人工边界可以吸收不同频率和任意角度入射波的原理以及PML人工边界的构造方法. 在此基础上,将PML人工边界应用于地震波动数值模拟的速度应力混合有限元格式中,探讨了PML应用的可行性,并通过数值试验研究了PML人工边界的反射率,比较了PML人工边界与多次透射公式(MTF)人工边界应用于体波和面波模拟中数值反射的差异,对两种边界的透射效果进行了分析. 结果表明, 尽管数值离散后PML人工边界不再保持完美匹配特性,但PML人工边界在近场波动数值模拟中可获得比MTF人工边界更为理想的吸收效果,在角点透射、大角度掠射情形下尤为明显;PML人工边界在混合有限元格式的数值算法中,未见失稳等不良反应,比MTF人工边界有更好的稳定性;在合理选择参数的情况下,PML人工边界的运算量可接受.      

透射边界真实反射系数的研究

本文在全部入射角范围内系统地讨论了时域局部透射边界真实的反射系数及其与理想稳态和理想暂态反射系数的关系。文中指出，后两者可分别用作前者上，下限的估计；同时，就时域逐步积分计算而言理想暂态较之理想稳态更接近真实情形，简便的理想暂态反射系数可用于在全部入射角范围内度量边界是的误差。     

波动数值模拟中透射边界的稳定实现

从波动能量在计算区内累积增大的观点出发,通过简单的一维弹性波模型,系统地阐明了在近场波动数值模拟中透射边界两类数值失稳--"高频振荡"和"零频飘移"的机理:前者源于对波动数值模拟无意义的高频波动在人工边界上的放大和波在有限计算区内多次反射产生的反复放大;后者则源于透射边界允许零频和接近零频的分量不断进入计算区. 由此提出了稳定实现透射边界的完整方案包括两项简单措施:第一,在全部计算区内按文中建议的方法注入小阻尼,以消除高频振荡;第二, 给出一种具有明确物理意义的消除零频飘移的算子算法. 最后,提供了三维波源问题和散射问题的详细数值试验结果.     

波动数值模拟中透射边界的稳定实现

从波动能量在计算区内累积增大的观点出发,通过简单的一维弹性波模型,系统地阐明了在近场波动数值模拟中透射边界两类数值失稳--"高频振荡"和"零频飘移"的机理:前者源于对波动数值模拟无意义的高频波动在人工边界上的放大和波在有限计算区内多次反射产生的反复放大;后者则源于透射边界允许零频和接近零频的分量不断进入计算区. 由此提出了稳定实现透射边界的完整方案包括两项简单措施:第一,在全部计算区内按文中建议的方法注入小阻尼,以消除高频振荡;第二, 给出一种具有明确物理意义的消除零频飘移的算子算法. 最后,提供了三维波源问题和散射问题的详细数值试验结果.     

场地地震波动模拟中透射边界稳定性问题研究

本文对场地地震波动模拟中透射边界的稳定性问题进行了一系列研究。对已有的透射边界稳定措施做了比较分析,并给出了不同措施的适用性。结合粘弹性边界的物理模型提出了一个新的消除透射边界飘移失稳的措施,参考高频滤波的思路探讨了一种新的消除透射边界葛频振荡失稳的措施,并用数值试验验证了这两种新的透射边界稳定方法的有效性。将显式行限元一有限差分方法结合透射边界的数值解法的计算结果与解析结果作对比分析,验证这弹一套数值解法的计算精度。利用二维模型模拟了自贡土层和山脊地形的场地效应,利用显戈有限元一有限差分方法结合透射边界这套数值解法进行数值模拟计算,应用汶川地震的实示强震观测记录来验证这套数值解法对于实际地震波动模拟的有效性。     

局部平面波假定与优化傍轴透射边界

本文研究了透射边界及相关问题,指出在近源波动的一般情况下,沿人工边界法向传播的外行波动应为曲面波,但在人工边界点的足够小的时空邻域内,可以将曲平面波近似处理为平面波,亦即采用了局部平面波假定。认为对不同的人工边界点,局部平面波模型中的基本变量是时空变量,是影响透射边界近源精度和稳定性的决定性因素。建立了对物理视波速进行优化求解的数学模型,以及相应的1阶和2阶优化透射边界。     

针对场地地震反应分析的ANSYS二次开发

以往的场地地震反应分析程序往往缺乏很好的前后处理,难以应用于实际工程。本文利用通用有限元分析程序ANSYS进行二次开发,引入多次透射边界以适用于场地地震反应分析。