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

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

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

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

三种介质耦联系统动力反应分析的显式有限元法及其应用

在笔者已建立的流体饱和多孔介质动力分析的显式有限元法的基础上，提出了可分析任意形状的流体饱和多孔介质-单相弹性固体介质-理想流体介质耦联的复杂系统的动力响应的显式有限元方法．该方法建立的有限元方程列式具有解耦特征，不需求解联立方程组，因而极大地提高了计算效率．将这一方法用于分析了考虑库水、坝、淤泥层和基岩这一复杂系统动力相互作用问题的斜坝面的地震响应，并给出了一些计算结果．      

波在饱和多孔介质与弹性固体介质交界面上的界面效应

本利用Ｂｉｏｔ饱和多孔介质动力学理论及饱和多孔介质和固体介质交界面的连续条件,推导了平面ＰＩ波以任意角饱和多孔介质向弹性固体介质入射时,反射快Ｐ１波,慢ＰＨ波和ＳＶ波和透射Ｐ波、ＳＶ波系数的计算公式,并给出一算例,本不考虑饱和多孔介质中液体渗流的能量耗散。     

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

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

基于流固耦合动力模型的饱和土体-隧道体系地震反应研究

基于ABAQUS软件平台,应用自行开发的流固耦合动力模型孔压单元模拟场地土体,并通过黏弹性人工边界方法实现地震动的输入,对饱和土体场地中的双孔隧道结构在地震荷载作用下的动力反应进行研究。计算结果表明:在地震反应结束时刻,场地土体位移幅值在两隧道之间以及两隧道的附近区域较大,而远离隧道的区域则较小;场地底部区域土体的孔压幅值较大,而场地顶部区域土体则较小;隧道左右两侧拱腰部位的衬砌的应力较大,而拱顶部位则较小。计算结果同时表明了流固耦合动力模型孔压单元在饱和土体-隧道体系地震反应研究中的适用性。     

流体饱和多孔隙介质弹性波方程边界元解法研究

基于流体饱和多孔隙各向同性介质模型,本文首先推导了流体饱和多孔隙介质中弹性波传播的频率域系统动力方程及边界积分方程,然后给出了流体饱和多孔隙介质弹性波方程的基本解,最后,利用本文给出的边界元方法对流体饱和多孔隙各向同性介质中的弹性波传播进行了数值模拟.结果表明：不论是从固相位移,还是液相位移的地震合成记录都能看到明显的慢速P波,本文提出的流体饱和多孔隙介质弹性波边界元法是有效可行的.     

平面S波在饱和多孔热弹性介质边界的反射问题研究

针对饱和多孔介质中热弹性波的传播特性问题,基于多孔介质理论和广义的热弹性模型,研究平面S波在饱和多孔热弹性介质边界上的反射问题。以考虑流-固耦合的饱和多孔介质波动方程和热-弹耦合的广义热弹性基本方程出发,建立饱和多孔介质的热-流-固耦合弹性波动模型。通过引入势函数并考虑自由透水和绝热的边界条件,经过理论推导最终给出在饱和多孔热弹性介质边界上的四种反射波的振幅反射率的理论表达式。在此基础上进行数值计算,分别讨论平面S波的入射频率、入射角和热膨胀系数等参数对四种反射波的振幅反射率的影响情况。结果表明:各反射波的振幅反射率分别随频率和热膨胀系数的增大而增大,同时也受到平面S波入射角变化的影响。该结论对于土动力学的理论研究及其相关的工程勘探具有一定的指导意义。     

黏弹性双相介质参数反演的同伦方法

本文以介质的动态时域响应（位移、速度、加速度）作为反演的依据,对黏弹性双相介质的材料参数进行了反分析研究.根据附加条件,并考虑到由于测量误差和测量数据的不充分造成的反演问题的不适定性,将黏弹性双相介质的参数反演问题转化为非线性泛函的极小值问题,然后应用大范围收敛的同伦方法求出使泛函极小的根作为反问题的解.最后,以二维半空间黏弹性双相介质模型为例进行了数值反演分析,数值结果表明本文方法在参数的反演过程中目标函数稳定收敛,且具有一定的抗“噪声”干扰性能;与利用速度和加速度信息反演的结果相比,利用位移信息反演的结果具有较高的精度和抗噪声能力.     

A viscous-spring transmitting boundary for cylindrical wave propagation in saturated poroelastic media

Based on the u–U formulation of Biot equation and the assumption of zero permeability coefficient, a viscous-spring transmitting boundary which is frequency independent is derived to simulate the cylindrical elastic wave propagation in unbounded saturated porous media. By this viscous-spring boundary the effective stress and pore fluid pressure on the truncated boundary of the numerical model are replaced by a set of spring, dashpot and mass elements, and its simplified form is also given. A u–U formulation FEA program is compiled and the proposed transmitting boundaries are incorporated therein. Numerical examples show that the proposed viscous-spring boundary and its simplified form can provide accurate results for cylindrical elastic wave propagation problems with low or intermediate values of permeability or frequency content. For general two dimensional wave propagation problems, spuriously reflected waves can be greatly suppressed and acceptable accuracy can still be achieved by placing the simplified boundary at relatively large distance from the wave source.     

Coupling multiscale finite element method for consolidation analysis of heterogeneous saturated porous media

The multiscale finite element method is developed for solving the coupling problems of consolidation of heterogeneous saturated porous media under external loading conditions. Two sets of multiscale base functions are constructed, respectively, for the pressure field of fluid flow and the displacement field of solid skeleton. The coupling problems are then solved with a multiscale numerical procedure in space and time domain. The heterogeneities induced by permeabilities and mechanical parameters of the saturated porous media are both taken into account. Numerical experiments are carried out for different cases in comparison with the standard finite element method. The numerical results show that the coupling multiscale finite element method can be successfully used for solving the complicated coupling problems. It reduces greatly the computing effort in both memory and time for transient problems.     

三维土-结构动力相互作用体系分析的两步时域显式波动有限元过程

为减少直接分析三维大尺度复杂土-结构动力相互作用问题的计算量,提高计算效率,本文直接从波动方程出发,提出了较常规子结构法更简单的两步简化计算过程,即第一步简化上部复杂结构体系为集中质量杆系模型,并求基础处等效输入,第二步通过等效输入求上部结构各位置的动力反应.其中第一步计算主要采用集中质量显式有限单元法结合局部透射人工...     

Analysis of Rayleigh waves in saturated porous elastic media by finite element method

A numerical procedure for the analysis of Rayleigh waves in saturated porous elastic media is proposed by use of the finite element method. The layer stiffness matrix, the layer mass matrix and the layer damping matrix in a layered system are presented for the discretized form of the solid-fluid equilibrium equation proposed by Biot. In order to consider the influence of the permeability coefficient on the behavior of Rayleigh waves, attention is focused on the following states: ‘drained’ state, ‘undrained’ state and the states between two extremes of ‘drained’ and ‘undrained’ states. It is found from computed results that the permeability coefficient exerts a significant effect on dispersion curves and displacement distributions of Rayleigh waves in saturated porous media.     

3D viscous-spring artificial boundary in time domain

After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.     

An explicit finite element method for dy—namic analysis in three—medium coupling system and its application

In this paper,an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium,elastic single-phase medium and ideal fluid medium.This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling,and which does not need to solve sys-tem of linear equations.The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water,dam,sediment and basement rock.The vertical displacement at the top point of the dam is calculated and some conclusions are given.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性, 文中选用二次场算法, 将背景电阻率设置为水平层状且各向异性, 场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组, 通过Galerkin加权余量法形成大型稀疏有限元方程, 采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位; 进而, 利用滑动平均方法得到二次矢量位和标量位在空间的导数, 由此得到二次电磁场; 通过一维模型对算法的可靠性进行验证, 与此同时, 针对实际复杂海洋电磁模型, 比较有限元模拟结果与积分方程模拟结果, 进一步验证算法精度.若干计算结果均表明, 文中算法具有良好的通用性, 适用于井中电磁、航空电磁, 环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.