基于黏弹性边界的地下隧道非线性地震响应分析模型

采用黏弹性边界方法和地震动输入等效节点力方法模拟半无限场地的地震激励,结合通用有限元软件ANSYS,建立地下隧道地震响应分析模型。模型采用等效线性化方法模拟土体非线性,采用瑞利阻尼和材料阻尼相结合方法模拟土体阻尼,分别考虑衬砌中混凝土和钢筋的非线性,以及衬砌和土体之间的接触。通过与文献比较,对整体模型的精度进行验证,并以天津市深厚覆盖层场地中某地铁双隧道为例进行非线性地震响应分析。研究表明,隧道间距较小时隧道之间存在显著相互作用,非线性地震动应力显著增大。本文模型可用于地下隧道一类结构的抗震验算分析。     

考虑土体非线性及倾角影响的盆地地震响应研究

基于二维沉积盆地模型,采用D-P弹塑性模型模拟盆地的非线性特征。利用显式有限元与黏弹性边界结合的方法,通过改变盆地边缘倾角,在时域和频域内分析盆地地表的地震动响应,对比线性与非线性盆地地震反应的差异。结果表明:(1)土体非线性对整个盆地范围内地震动的影响都较显著。考虑非线性时地震动放大系数明显降低,降低幅度在30%~50%。同时,考虑非线性和倾角影响时最强烈放大区域的范围和位置变化,且很小倾角下的分布特征显著不同。(2)两分量的放大系数都有随边缘倾角的增大而增强的趋势,但均是盆地边缘区域受非线性的影响最为显著。此外,真实地震波输入下显著放大区域的范围及线性与非线性结果的差异程度相对更大。(3)考虑非线性时,对于不同频率地震波的放大系数差别明显,但都表现出从低频到高频谱比分布越来越复杂的现象,同时盆地倾角的影响程度随频率的增大更加明显。(4)考虑土体非线性并未改变地震波传播的总体特征,但各震相强度相对降低。     

基于黏弹性模型的苍山地震库仑应力变化研究

采用多层黏弹性模型,计算郯城地震和苍山地震引起的同震及震后黏弹松弛库仑应力变化,讨论这两个地震之间的触发关系、苍山地震对余震分布及周围地区应力变化的影响。结果表明:(1)郯城地震对苍山地震的发生有一定的抑制作用;(2)苍山地震的大部分余震分布在库仑应力增加区域,库仑应力变化值为0.08~0.12bar,且余震主要分布在8~14km处的断层面上;(3)安丘—莒县断裂的东北部、苍尼断裂的中段、凫山断裂的东段以及韩庄断裂、峄城断裂位于库仑应力增加区,均存在较强的地震危险性。     

流体饱和多孔介质黏弹性动力人工边界

基于Biot流体饱和多孔介质本构方程，采用平面波和远场散射波经验叠加来反映外行波传播，以经验参数反映人工边界外行波动的衰减和多角度透射特性。在人工边界处分别施加反映固相和液相介质传播效应的弹簧及阻尼来模拟人工边界以外的无限域介质对来自有限域的外行波的能量的吸收作用。从而形成一种流体饱和多孔介质的黏弹性动力人工边界。数值算例表明：边界的精度和稳定性高于现有的黏性边界、黏弹性人工边界及一阶透射边界。     

基于黏弹性边界的地震动输入方法和边界条件选择研究

针对地下工程抗震分析中边界条件和地震动输入方式选择问题,对一地下半无限体模型在一剪切波垂直入射条件下进行不同边界条件和不同地震动输入方法的算例试算。结果发现:两侧采用竖直向位移约束、底部采用黏弹性边界、并在模型底部用等效节点力输入地震动的方法,与三边均采用黏弹性边界条件和三边均采用等效节点荷载输入地震动的方法可以获得同样的结果。该方法易于在ANSYS中通过APDL编程实现,可为地下工程抗震分析中黏弹性边界的使用提供便利。     

黏弹性人工边界在ADINA中的应用

探讨利用ADINA 的弹簧阻尼单元模拟黏弹性人工边界的方法,并进行了均匀半空间算例与成层半空间的计算,验正了采用的黏弹性人工边界具有较高的精度,可以方便的应用于波动问题的模拟解析,从而可以有效地分析结构-地基的动力相互作用问题.     

饱和岩石的各向异性及非线性黏弹性响应

在MTS伺服压机上对不同饱和状态的砂岩、大理岩标本进行了垂直层理和平行层理两个方向的正弦波加载试验,研究了饱和岩石的各向异性及非线性黏弹性行为. 在最大载荷低于岩石的屈服点时,获得的衰减、杨氏模量、泊松比、波速都表现出显著的各向异性,同时它们还具有较明显的应变振幅效应和频率效应. 随着应变振幅的增大,衰减近似呈线性增长,杨氏模量和泊松比近似呈线性下降. 在0.005-4Hz频段内除衰减不明显依赖频率以外,杨氏模量、泊松比、波速对频率的依赖性都较强; 当频段扩大到15Hz时,这4个参数都表现出较强的频率效应. 饱和岩石的衰减、杨氏模量、泊松比、波速均随饱和液体的黏滞系数增大而增大.     

基于黏弹性边界的P波斜入射波场分解方法精度比较

黏弹性人工边界以其概念清晰,稳定性高,精度良好的优点,在地下结构抗震的研究中已经得到了较为广泛的应用。在使用黏弹性边界进行地震动的输入时,需要把地震波转化为模型网格节点上的等效节点力,这个过程就要涉及到波场分解的问题。波场分解的方式可以不同,但不同的波场分解方式会得出精度不同的计算结果,而学界目前尚未就此达成一致。本文在前人研究的基础上,总结了3种地震波斜入射的波场分解方法,编制了相应的Python脚本来实现黏弹性边界和等效节点力的批量施加,并且以不同角度的地震波斜入射为例在同一有限元模型中计算以进行比较研究。结果表明计算精度与波场分解的方式密切相关:在模型边界上考虑的自由场波动叠加越完整,计算结果越精确。     

饱和多孔介质近场波动分析的一种黏弹性人工边界

采用有限模型数值方法求解能量开放系统中的波动问题时,虚拟人工边界的处理方式对计算结果的准确性和精度具有重要的影响.本文针对无限域饱和多孔介质中波传播问题的人工边界处理方式进行了研究,提出了饱和多孔介质近场波动分析的一种黏弹性人工边界处理方法.在考虑多孔介质中固相和液相的相互作用的情况下,通过在人工边界处分别施加反映固相和液相介质波传播效应的弹簧及阻尼来模拟饱和多孔介质中波的能量辐射效应影响.算例表明,本文建议的黏弹性人工边界具有较好的模拟效果.     

饱和岩石的各向异性及非线性黏弹性响应

在MTS伺服压机上对不同饱和状态的砂岩、大理岩标本进行了垂直层理和平行层理两个方向的正弦波加载试验,研究了饱和岩石的各向异性及非线性黏弹性行为. 在最大载荷低于岩石的屈服点时,获得的衰减、杨氏模量、泊松比、波速都表现出显著的各向异性,同时它们还具有较明显的应变振幅效应和频率效应. 随着应变振幅的增大,衰减近似呈线性增长,杨氏模量和泊松比近似呈线性下降. 在0.005-4Hz频段内除衰减不明显依赖频率以外,杨氏模量、泊松比、波速对频率的依赖性都较强; 当频段扩大到15Hz时,这4个参数都表现出较强的频率效应. 饱和岩石的衰减、杨氏模量、泊松比、波速均随饱和液体的黏滞系数增大而增大.     

基于有限元方法研究岫岩罗圈里陨石坑地形对地震反应的影响

基于有限元方法对岫岩罗圈里陨石坑地形对地震反应的影响进行了探索。结合文献资料、岩土勘察工作建立了岫岩罗圈里陨石坑二维有限元分析模型。选用粘弹性人工边界和等效节点力输入方法解决了有限元模拟半空间无限域时的边界问题和地震动输入问题。为了突出地形因素的影响,采用弹性材料介质进行数值模拟分析,并与水平场地模型的分析结果进行了对比。数值模拟结果和对比表明:罗圈里陨石坑地形对于陨石坑中部的地震动有显著的放大作用,远高于水平场地模型。此外,盆地中部短周期部分的反应谱也明显高于盆地边缘和水平一维场地模型,这是导致海城地震中盆地中部自振周期较短的单层房屋破坏明显严重的主要原因。     

弹性波边界元法正演模拟

弹性波边界元地震模型方法（ＢＥＥＳＭ）,实现了二维和三维问题的纵、横波及转换波的同时模拟,并且能模拟任意复杂构造的地震声波正演模型．根据地震模型的特点,本文发展了数值积分计算与矩阵消元同步进行的块状高斯消元法;用解析法处理奇异积分;用无限元法处理边界吸收问题;采用单元长度随介质速度和计算频率变化的变单元算法,及自动剖分单元等技术,提高了计算精度,节省了内存,缩短了计算时间．     

基于正交试验法的随机地震响应的参数分析

考虑地震的随机性和土体的非线性,引入粘弹性边界和正交试验法,提出了一种在谐响应分析中直接输入加速度的加载方式,建立了基于虚拟激励法、等效线性化法和确定性动力有限元法相结合的非线性桩-土体系的随机地震响应的参数敏感性分析方法.研究了桩土模量比、地震动输入水平、桩的长径比等无量纲参数对随机响应的影响.结果表明:桩土模量比是最敏感的因素.     

用线性波传播模型再现视非线性地震响应

在地震工程中,地震荷载的预测了解沉积层地面运动的放大作用,它定义为沉积 与岩石场地地面运动大小的比值。在１９９４年北岭地震中,观测到的沉积层场地余震的（级为３－４）弱地面运动放大作用是主震（震级为６．７）放大作用的２倍。这种现象的一种解释是沉积层振幅的放大（非线性）作用。然后,利用经验脉冲响应进行的地震模拟和用弱非生,随机三维地壳速度弹性有限差分法的计算表明,线性波传播模型可以解释所观测到的（表现     

饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法

本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.     

衰减补偿的黏弹性最小二乘逆时偏移

地震勘探是发现和认识油气藏的重要方法,随着国内油气勘探领域、勘探对象的日趋复杂化,油气勘探迫切需求地震成像的精准度、分辨率和保真度的大幅度提升.相比于常规逆时偏移方法,最小二乘逆时偏移成像有较少的伪像,成像振幅均衡,成像分辨率高.目前常用弹性最小二乘逆时偏移方法来对地下介质进行成像.但由于地下介质通常表现出黏弹性性质,...     

弹性波边界元法正演模拟

弹性波边界元地震模型方法（ＢＥＥＳＭ），实现了二维和三维问题的纵、横波及转换波的同时模拟，并且能模拟任意复杂构造的地震声波正演模型．根据地震模型的特点，本文发展了数值积分计算与矩阵消元同步进行的块状高斯消元法；用解析法处理奇异积分；用无限元法处理边界吸收问题；采用单元长度随介质速度和计算频率变化的变单元算法，及自动剖分单元等技术，提高了计算精度，节省了内存，缩短了计算时间．     

二维复杂黏弹性各向异性介质中的地震射线追踪

实际地层中传播的地震波普遍存在速度各向异性和能量衰减现象,因此基于弹性介质假设条件下的地震射线正反演算法具有一定的局限性,而研究黏弹性各向异性介质中的地震波传播规律可为揭示地下结构提供更加可靠的理论依据.射线追踪技术是揭示高频地震波传播特性的有效手段之一,然而绝大多数研究仅限于弹性介质.针对黏弹性各向异性介质,本文首先给出了一种射线速度和振幅衰减的计算方法,然后结合改进型最短路径算法,在实空间内实现了计算复杂介质模型中地震波(qP,qSV,qSH)的射线路径和传播走时以及能量衰减(虚走时).该算法适用于复杂黏弹性各向异性介质.误差分析结果显示,实走时的最大相对误差小于0.13％,虚走时的最大相对误差小于0.55％,表明该算法具有较高的计算精度.     

城市地下变电站非线性地震反应分析

为研究城市地下变电站结构及其内部电气设备的动力特性及地震反应规律,以西安市快速轨道交通二号线张家堡地下变电站为工程背景,依据黏弹性人工边界理论及有限单元法,建立考虑土体-结构-设备动力相互作用的三维有限元计算模型,首先通过半空间自由场模型验证了地震动输入方法的有效性,继而开展地下变电站系统的模态分析和动力时程分析,得到了该类结构的动力特性及地震响应规律。研究结果表明,结构上覆土体及电气设备在地震作用下产生的惯性力是结构构件应力增长的主要原因,8度罕遇地震作用下地下变电站内电气设备的加速度放大系数明显超出规范要求,需要采取相应保护措施。     

2010年马乌莱MW8.8地震震后形变三维黏弹性数值模拟

2010年智利马乌莱MW8.8地震发生在纳斯卡板块与南美板块的板块边界处,引起了显著的同震和震后效应.GPS台网数据显示记录到的同震海向位移最大约5 m,垂向沉降最大约50 cm.在经过对俯冲效应、季节变化等效应的校正后,震后6年的海向最大位移约68 cm,垂向抬升最大约20 cm.马乌莱地震显著的震后形变对该区域的地...