首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 156 毫秒
1.
时域局部人工边界是以行波概念为基础的,原则上不能用于模拟零频分量(消减运动)。Wolf提出了能同时模拟零频分量和行波的双渐近-多向透射边界。本文指出,在这一组合中的双渐近边界的高频渐近项是不合理的,还查能降低模拟精度。据此,本文提出一种概念清楚、简便实用的组合方案-多次射透-静力边界,并通过数值试验对此方案进行检验。  相似文献   

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

3.
基于外行波局部法向透射概念,从多项式外推的角度,结合一维和广义二维几何解释,推导建立了外行波为平面波时,人工边界节点位移解的计算公式。对该位移解的精度进行了时域分析,论证了该位移解收敛于数值精确解的必要条件和线性外推位移解的精度与数值精确解匹配的充要条件。在此基础上提出了一次“精确”透射的概念和相应的数学方法,实现了该位移解逼近于数值精确解的优化计算。  相似文献   

4.
李小军  廖振鹏 《地震学报》1995,17(3):362-369
推导出了分析二维粘及弹性场地地形对地震动影响的显式有限元-有限差分方法.这一方法中,首先利用人工边界及有限元离散方法,给出问题分析的有限元离散网格计算力学模型,并利用一种类似于差分方法的有限元方法,建立局部网格节点的动力方程,而后利用笔者提出的有阻尼体系动力方程求解的显式差分格式,及推广的多次透射边界公式,给出网格节点运动量计算的时域显式逐步积分公式.利用计算机程序实现这一方法的计算具有所需计算机内存量小及计算时间量小的优势,而且,这一方法适用于任意地形情况,具有较高的计算精度及较好的计算稳定性.   相似文献   

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

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

7.
局部解耦的时域波分析方法   总被引:8,自引:1,他引:7  
从模拟散射波的观点出发,推导了波动问题局部解耦的人工边界条件,与有限元内点的时域显式方法相结合建立了波动分析的局部解耦的直接分析方法,文中还从波动分解的概念,建立了外源问题(如地震作用)输入方法。  相似文献   

8.
本文针对显式有限元-有限差分法进行了系统的研究,包括理论分析、算法改进以及实际应用.通过理论推导定量地研究了采用这种数值模拟方法对波动传播规律的影响,以及其在结合局部透射人工边界进行计算时的稳定条件.就采用该方法结合局部透射人工边界求解波源问题和散射问题的情况分别提出了改进计算的措施,并通过数值试验验证其有效性.  相似文献   

9.
结合均匀场地模型和分层均匀场地模型,采用数值试验的方法研究了场地地震反应分析中侧向人工边界距离的选取问题。分析了各主要因素对人工边界距离选取的影响,得到了如下结论:单侧人工边界距离与土层深度之比L/h是影响人工边界距离选取的主要参数;随着侧向人工边界距离的增加,有限元模型的计算误差并不一直减小,在局部范围内,会出现计算误差随着边界距离的增加而增大的现象;场地硬度、阻尼比和激振频率等因素不会影响计算误差随边界距离的变化规律,但采用自由场边界时受输入地震动频谱特性的影响较为明显;对水平激振下的一维场地模型来说,采用结合自由度边界时所需边界距离最短,自由场边界次之,采用粘性边界达到收敛时所需边界距离最长;无论哪种边界条件,当边界距离大于某一范围后,都能得到稳定的计算结果;相对于均质场地模型来说,实际的分层均匀场地建模时可适当减小侧向人工边界距离。  相似文献   

10.
地震波数值模拟的非规则网格PML吸收边界   总被引:14,自引:4,他引:10       下载免费PDF全文
徐义  张剑锋 《地球物理学报》2008,51(5):1520-1526
以格子法为基础,以声波方程为例研究非规则网格PML(Perfectly Matched Layer)方法.本方法的核心是建立局部坐标系下的分裂方程和基于积分近似的微分方程弱形式.该非规则网格模拟方法允许在计算域内设置任意形状的人工边界.对于二维半空间问题,与采用矩形人工边界相比,采用半圆形人工边界可减少计算量20%以上.采用光滑的曲边界,不仅可减少计算区域,还可避免常规的PML吸收边界在吸收带角点区域的特殊处理.本方法事先计算和存储边界单元的局部几何参数,在计算的每一时间步查表调用这些参数,与常规的直边界PML方法相比,不增加任何计算量.  相似文献   

11.
两相介质近场波动模拟的一种解耦有限元方法   总被引:6,自引:1,他引:5       下载免费PDF全文
本文将求解近场波动问题的一种解耦技术推广到两相介质,得到了一种求解两相介质近场波动问题的直接解耦方法,包括集中质量有限元模型、时域显式积分格式和局部人工边界条件. 首先应用加权残数法,并依据波动模拟的精度要求,得到了两相介质集中质量有限元模型. 然后,结合两相介质中波动的衰减特性,实现了透射边界在两相介质近场波动中的运用. 最后,通过数值实验,并与解析解对比,验证了本文方法的有效性.  相似文献   

12.
波动数值模拟是力学、地球物理和多个工程学科共同关注的领域,研究者在基础和应用研究方面皆取得了丰硕成果.然而,由边界引入局部失稳问题,即使在线性范围亦尚未彻底解决.这里“边界”泛指人工边界、物理边界或不同介质的分界面.  相似文献   

13.
本文将Newmark法的时间差分格式的有限元与边界元方程相结合,建立了动力问题的边界元与有限元耦合模型。这种耦合模型在求解时域范围内基础动力响应问题中,可充分发挥有限元求解局域问题和边界元求解无限域问题各自的优点。  相似文献   

14.
本文采用一维模型分析了时域离散化后有限元离散网格中波动的基本特征。文中讨论了时域离散化对波动的频散、截止频率及寄生振荡的影响,分析了由时域离散化引起的在频率域及由空间离散化引起的在波数域波动解的多重分支现象,并指出在研究中应取其基本分支。文中也给出了考虑时间步长影响的有限元离散化准则。研究结果表明,从模拟连续模型的精度上看,时域离散化使集中质最有限元法优于一致质量法,前者可取较大时间步长。从计算精度、容易程度及经济性各方面衡量,集中质量法均比一致质量法可取。本文的工作不但对波动有限元模拟的参数选取有指导意义,同时也是分析时域人工边界稳定性的基础工作之一。  相似文献   

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

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

17.
P-SV波斜入射时成层半空间自由场的时域算法   总被引:3,自引:0,他引:3  
刘晶波和王艳提出了一种弹性水平成层半空间中平面P-SV波斜入射时平面内自由波场时域计算的一维化有限元方法,该方法采用黏性人工边界条件近似地模拟下部基岩半空间的辐射阻尼,可导致平面P-SV波以大角度入射时自由场计算精度降低.本文提出一种精确模拟基岩半空间辐射阻尼的人工边界条件.由于基岩半空间中外行波是传播方向已知的平面P波和SV波,利用弹性介质的应力-位移本构关系建立了人工边界处应力与速度的阻抗边界条件;采用该人工边界条件替代黏性边界条件,改进了P-SV波斜入射时成层半空间自由场的时域算法.数值试验表明,与采用黏性边界条件的自由场算法相比,改进方法具有更高的计算精度,其计算结果与理论解吻合更好.  相似文献   

18.
二维场地地震反应有限元分析的问题探讨   总被引:2,自引:2,他引:0  
本文通过数值计算,在线弹性范围内初步分析了土层地震反应随不同侧向人工边界、有限土域取值范围、土介质的阻尼系数的变化情况;采用滞后阻尼假定,探讨了在土层时域分析中如何由滞后阻尼系数形成阻尼矩阵,给出了一个阻尼系数转换频率取值的经验公式。在此研究基础上编制了二维有限元SR2D计算程序。  相似文献   

19.
时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析.  相似文献   

20.
鉴于核电厂选址的愈发复杂多样性,合理有效地模拟复杂非均质地基条件下抗震I类物项结构的静动力耦合响应是进行抗震安全性评价的关键技术问题。以某核电厂联合泵房的实际厂址地基为背景,针对静动力分析中的边界处理问题,基于Newmark隐式积分特点和黏弹性人工边界理论,提出了一种适于静动力耦合分析求解的人工边界模型,并通过创建无厚度紧支黏弹性人工边界单元实现了所开发模型在ANSYS平台的嵌入。进而,通过联合泵房结构静力分析与极大初始时间步分析的对比,验证所提出模型的合理可靠性。最后,从结构应力、响应及地基稳定性等方面对联合泵房进行了综合评价,并就地基局部非均质的影响进行探讨。研究结果表明:实际厂址中局部捕虏体的存在对联合泵房结构的地震响应及地基安全性影响较小,且所建立模型具有较好地解决非线性动力相互作用问题的潜力。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号