共查询到20条相似文献，搜索用时 609 毫秒
1.
流体饱和多孔介质黏弹性动力人工边界 总被引：1，自引：0，他引：1
基于Biot流体饱和多孔介质本构方程，采用平面波和远场散射波经验叠加来反映外行波传播，以经验参数反映人工边界外行波动的衰减和多角度透射特性。在人工边界处分别施加反映固相和液相介质传播效应的弹簧及阻尼来模拟人工边界以外的无限域介质对来自有限域的外行波的能量的吸收作用。从而形成一种流体饱和多孔介质的黏弹性动力人工边界。数值算例表明：边界的精度和稳定性高于现有的黏性边界、黏弹性人工边界及一阶透射边界。 相似文献
2.
数值模拟是研究半无限域中结构地震反应的重要手段,而边界条件的选取则直接影响到数值模拟结果的可靠性.本文应用有限元软件ABAQUS,分别采用固定边界条件、黏性边界条件和远置边界条件对简单弹性半无限介质的波源问题和散射问题进行了数值模拟分析,结果显示:使用固定边界时,波动不能穿过边界面进入半无限介质,这与波动在半无限介质中... 相似文献
3.
P-SV波斜入射时成层半空间自由场的时域算法 总被引：3，自引：0，他引：3
刘晶波和王艳提出了一种弹性水平成层半空间中平面P-SV波斜入射时平面内自由波场时域计算的一维化有限元方法,该方法采用黏性人工边界条件近似地模拟下部基岩半空间的辐射阻尼,可导致平面P-SV波以大角度入射时自由场计算精度降低.本文提出一种精确模拟基岩半空间辐射阻尼的人工边界条件.由于基岩半空间中外行波是传播方向已知的平面P波和SV波,利用弹性介质的应力-位移本构关系建立了人工边界处应力与速度的阻抗边界条件;采用该人工边界条件替代黏性边界条件,改进了P-SV波斜入射时成层半空间自由场的时域算法.数值试验表明,与采用黏性边界条件的自由场算法相比,改进方法具有更高的计算精度,其计算结果与理论解吻合更好. 相似文献
4.
《地震工程与工程振动》2016,(2)
在波动方程有限差分波场数值模拟中,为了使计算得到的边界波场值更接近于真实的边界波场值,需要使用边界条件以减少来自计算区域边界的人为反射能量。传统的分裂式完全匹配层(SPML)吸收边界条件不能有效吸收掠射波,后来发展的不分裂卷积完全匹配层(CPML)能够较好地处理这个问题,并且CPML在处理边界问题时无需对波场进行非物理分裂。本文针对Kelvin-Voigt模型的黏弹性介质,采用高阶交错网格有限差分方法进行数值模拟,同时利用CPML吸收边界条件进行边界处理。数值模拟结果表明:与弹性波相比,大角度入射产生的低频掠射波对黏弹性波的影响更大;CPML吸收边界条件对低频掠射波的吸收效果比传统的SPML吸收边界条件更好。 相似文献
5.
传统的人工边界是在声波或弹性波分解为内行和外行分量的基础上设计的.对于较复杂的介质模型（如粘弹介质、弹塑介质、塑性介质模型等）,地震波不能简单地分解为各自满足其运动方程的纵波、横波等成分,且地震波在传播过程中还有衰减现象,因此传统的人工边界已不能有效地用于非弹性的较复杂介质的波场模拟.本文提出了一种适合于更一般的地震模型的波场计算的透射边界条件,并以粘弹介质的一个较简单的二维几何模型为例,通过有限元数值求解这种透射边界和刚性边界下粘弹波的激发和传播,得到地面合成地震记录,证明了这种透射边界用于较复杂介质波场计算是可行的. 相似文献
6.
针对地下工程抗震分析中边界条件和地震动输入方式选择问题,对一地下半无限体模型在一剪切波垂直入射条件下进行不同边界条件和不同地震动输入方法的算例试算。结果发现:两侧采用竖直向位移约束、底部采用黏弹性边界、并在模型底部用等效节点力输入地震动的方法,与三边均采用黏弹性边界条件和三边均采用等效节点荷载输入地震动的方法可以获得同样的结果。该方法易于在ANSYS中通过APDL编程实现,可为地下工程抗震分析中黏弹性边界的使用提供便利。 相似文献
7.
采用有限模型数值方法求解能量开放系统中的波动问题时,虚拟人工边界的处理方式对计算结果的准确性和精度具有重要的影响.本文针对无限域饱和多孔介质中波传播问题的人工边界处理方式进行了研究,提出了饱和多孔介质近场波动分析的一种黏弹性人工边界处理方法.在考虑多孔介质中固相和液相的相互作用的情况下,通过在人工边界处分别施加反映固相和液相介质波传播效应的弹簧及阻尼来模拟饱和多孔介质中波的能量辐射效应影响.算例表明,本文建议的黏弹性人工边界具有较好的模拟效果. 相似文献
8.
无限域地基有限元分析的简化粘弹性边界 总被引：4，自引：0，他引：4
本文提出一种能用于无限域地基动、静力有限元分析的简单人工边界——粘弹边界，采用一层边界单元来实现。该边界弹簧部分采用弹性半空间Mindlin解计算，边界单元采用具有粘性阻尼的粘弹材料，通过确定合理的弹性模量、泊松比和粘性阻尼系数，不仅可以使边界单元能够模拟无限域远场的弹性恢复力，同时可以模拟粘性边界。通过算例分析，验证了该边界的有效性，它可用于土-结构动力相互作用分析。 相似文献
9.
利用有限元软件Ansys二次开发语言apdl进行二次开发,将多次透射边界(MTF边界)和弹阻边界(V-S边界)添加到软件中,从而实现对土-结构动力相互作用问题中无限域场地的模拟。以野外大比例(1/2)土-箱基-框架结构试验模型的牵引释放自由振动试验(波源问题)和爆破震动试验(散射问题)数据为基础,通过有限元数值计算结果和试验结果的对比可知,施加了人工边界的土体模型可以较好地模拟波在无限域土体中的传播,有效地减少有限元计算模拟时的计算规模。同时多次透射边界在计算精度以及散射问题中地震波输入的便捷性方面要优于弹阻边界。 相似文献
10.
11.
基于饱和两相介质弹性波动方程分析SV波在饱和土体自由表面的反射问题,引入波动方程的势函数解答,求解出二维问题中SV波入射情况下饱和土体自由场的位移、速度、加速度和应力响应。在饱和土体自由场响应解析解基础上,建立SV波入射下饱和土体自由场静、动力有限元模型。建模中考虑了如下几方面因素:(1)在不同分析步,对土体单元赋予不同材料本构。通过*model change命令进行单元生死设定,从而实现在初始应力场平衡的静力状态下采用DuncanChang本构模型,而地震波动输入时采用Davidenkov动力本构模型;(2)采用多孔介质黏弹性人工边界条件,在人工边界上分别施加固相和液相介质的弹簧和阻尼来模拟饱和土体中能量的传播;(3)将地震波转化为作用在人工边界上的等效地震荷载,施加到人工边界节点上;(4)土体单元采用4结点平面应变孔压单元(CPE4P)。有限元计算与解析解比较结果表明:SV波在垂直入射和掠入射时,竖向位移响应为零;在45°左右入射时,水平位移响应最大;60°左右入射时,竖向位移响应最大。这些结论与解析解吻合较好,本文模型为建立土-结构动力相互作用模型打下良好的基础。 相似文献
12.
在地震波传播数值模拟的过程中,需要使用吸收边界条件从而达到衰减人为边界 反射的目的. 本文基于傍轴近似法提出了计算三维弹性波方程的吸收边界条件公式,表示了 各边界面、边棱和角点处波场所满足的单程波方程,并在三维弹性波数值模拟中进行了应用 . 理论模型及三维盐丘地质模型波场切片快照试算结果表明,该吸收边界条件可以有效地吸 收人为边界反射,适用于较大入射角情况,从而消除了边界有效波信息的干扰. 由于采用四 阶近似方程,在保证计算精度的前提下,该方法具有节省计算工作量和易于实现的特点. 相似文献
13.
14.
基于薄层法在研究层状介质中波的传播问题的高效性,边界单元法处理无限域问题的精确性,结合两者提出层状地基二次形函数薄层法基本解答, 建立了基于层状半空间薄层法位移基本解答的半解析动力边界元法,可以有效的分析多层场地的动力问题,解决土-结构动力相互作用问题.同时对上软下硬的典型层状地基中的空沟被动隔振效果进行了详细的分析.研究表明:对于上软下硬地基中空沟被动隔振,为取得与匀质弹性半空间中连续屏障隔振同样的效果,需要更深的深度,同时振幅衰减曲线会出现明显的振荡现象;为使隔振效果满足Richart要求,当表层土归一化深度H＜1.0时,空沟归一化深度T应大于1.3,而1.0≤H≤2.0时,T应当大于1.8. 相似文献
15.
《应用地球物理》2019,(4)
无网格有限差分法能有效提高数值模拟的几何灵活性,且无需网格映射或复杂的网格生成过程。RBF-FD (基于径向基函数的有限差分)是最常用的无网格有限差分法之一,可以准确模拟地震波在非矩形计算域中的传播。本文提出适于弹性波方程无网格有限差分数值解的PML (完全匹配层)吸收边界条件,可以应用于非矩形速度模型的边界。但是PML吸收边界对近掠射波、低频波的吸收效果不好。为此,我们继续提出适于弹性波方程无网格有限差分数值解的CFS-PML(复频移完全匹配层)吸收边界条件。本文所提两种边界条件均是通过构造辅助微分方程,得到不分裂时域表达式,具有存储量小、便于编程实现的特点。模拟结果表明,两种PML吸收边界条件都能有效地消除无网格有限差分数值模拟的人工边界反射。此外,本文所提CFS-PML相较PML对近掠射波和损耗波的吸收效果更好。 相似文献
16.
基于单相介质中地震波理论的高频面波法已广泛应用于求取浅地表S波的速度.然而水文地质条件表明,普遍的浅地表地球介质富含孔隙.孔隙中充填的流体会显著地影响面波在浅地表的传播,进而造成频散和衰减的变化.本文研究了地震勘探频段内针对含流体孔隙介质边界条件的面波的传播特性.孔隙流体在自由表面存在完全疏通、完全闭合以及部分疏通的情况.孔隙单一流体饱和时,任何流体边界条件下存在R1模式波,与弹性介质中的Rayleigh波类似,相速度稍小于S波并在地震记录中显示强振幅.由于介质的内在衰减,R1在均匀半空间中也存在频散,相速度和衰减在不同流体边界下存在差异.Biot固流耦合系数(孔隙流体黏滞度与骨架渗透率之比)控制频散的特征频率,高耦合系数会在地震勘探频带内明显消除这种差异.介质的迂曲度等其他物性参数对不同流体边界下的R1波的影响也有不同的敏感度.完全闭合和部分疏通流体边界下存在R2模式波,相速度略低于慢P波.在多数条件下,如慢P波在时频响应中难以观察到.但是在耦合系数较低时会显现,一定条件下甚至会以非物理波形式接收R1波的辐射,显示强振幅.浅表风化层低速带存在,震源激发时的运动会显著影响面波的传播.对于接收点径向运动会造成面波的Doppler频移,横向运动会造成面波的时频畸变.孔隙存在多相流体时,中观尺度下不均匀斑块饱和能很好地解释体波在地震频带内的衰减.快P波受到斑块饱和显著影响,R1波与快P波有更明显关联,与完全饱和模型中不同,也更易于等效模型建立.频散特征频率受孔隙空间不同流体成分比例变化的控制,为面波方法探测浅地表流体分布与迁移提供可能性.通常情况孔隙介质频散特征频率较高,标准线性黏弹性固体可以在相对低频的地震勘探频带内等效表征孔隙介质中R1波的传播特征,特别在时域,可在面波成像反演建模中应用. 相似文献
17.
基于Kelvin空间半无限体黏弹性解,以单个轮轴荷载为例,研究地铁列车荷载下地基土动应力状态、应力路径及主应力轴旋转等应力变化规律,并与相应弹性解下的各动应力特征进行对比分析,探讨更适用于工程实际的地基附加动应力特征解。结果表明:弹性地基与黏弹性地基中的各应力状态变化总体趋势、主应力轴旋转变化趋势相同,但极限应力状态各应力分量大小及应力路径形状有明显差别,考虑黏弹性解更符合软土区地铁长期变形实际工况;并对单轮荷载下动应力特征的黏弹性解进行参数分析,发现特征系数对应力分量影响具有不同特性。 相似文献
18.
三维各向异性介质弹性波模拟的透射边界条件 总被引：1，自引：1，他引：0
在弹性波数值模拟中,为了压制人工边界反射,需要使用边界条件.本文从三维各向异性介质中弹性渡方程出发,在边界入射波为平面波的假设前提下,利用特征值分解以及坐标系旋转,导出人工边界处质点位移矢量计算公式,实现了在边界处为任意倾斜入射时的透射边界条件.并显示了边界条件的使用效果. 相似文献
19.
基于薄层法在研究层状介质中波的传播问题的高效性、边界单元法处理无限域问题的精确性,结合二者的优点提出三维层状地基薄层法基本解答,建立了基于三维层状半空间薄层法位移基本解答的半解析动力边界元法。该方法可有效的分析多层场地的动力问题,解决土-结构动力相互作用问题。同时分别对粘弹性上软下硬地基及上硬下软地基的三维空沟主动隔振进行了详细的分析。结果表明,两种情况下采用空沟屏障隔振均可以取得一定的隔振效果;同时,地基分层参数对空沟隔振体系的隔振效果影响显著。 相似文献
20.
将土体视为固-液两相介质,基于饱和土体有效应力原理,建立饱和土体-地下综合管廊结构体系相互作用动力模型：在地应力平衡的静力状态下,采用Duncan-Chang非线性弹性本构模型,在地震波作用的动力状态下,采用Davidenkov非线性黏弹性本构模型;考虑饱和土体黏弹性动力人工边界条件,并将地震动作用转化为作用在人工边界节点上的动力荷载。模型考察不同土体材料、结构特性以及土-结构接触摩擦对结构地震响应的影响,得出如下结论：（1）地震波的卓越周期与场地卓越周期相近时,引起结构上的变形最大;（2）综合管廊结构管廊壁厚越薄,埋深越深,结构尺寸越大,结构刚度越小,结构变形越大;（3）不考虑土-结构接触面的状态非线性将会增大结构变形。 相似文献