Analysis and design of open trench barriers in screening steady-state surface vibrations

The problem of vibration isolation by rectangular open trenches in a plane strain context is numerically studied using a finite element code, PLAXIS. The soil media is assumed to be linear elastic, isotropic, and homogeneous subjected to a vertical harmonic load producing steady-state vibration. The present model is validated by comparing it with previously published works. The key geometrical features of a trench, i.e., its depth, width, and distance from the source of excitation, are normalized with respect to the Rayleigh wavelength. The attenuation of vertical and horizontal components of vibration is studied for various trench dimensions against trench locations varied from an active to a passive case. Results are depicted in non-dimensional forms and conclusions are drawn regarding the effects of geometrical parameters in attenuating vertical and horizontal vibration components. The screening efficiency is primarily governed by the normalized depth of the barrier. The effect of width has little significance except in some specific cases. Simplified regression models are developed to estimate average amplitude reduction factors. The models applicable to vertical vibration cases are found to be in excellent agreement with previously published results.     

层状地基空沟被动隔振分析

基于薄层法在研究层状介质中波的传播问题的高效性,边界单元法处理无限域问题的精确性,结合两者提出层状地基二次形函数薄层法基本解答, 建立了基于层状半空间薄层法位移基本解答的半解析动力边界元法,可以有效的分析多层场地的动力问题,解决土-结构动力相互作用问题.同时对上软下硬的典型层状地基中的空沟被动隔振效果进行了详细的分析.研究表明:对于上软下硬地基中空沟被动隔振,为取得与匀质弹性半空间中连续屏障隔振同样的效果,需要更深的深度,同时振幅衰减曲线会出现明显的振荡现象;为使隔振效果满足Richart要求,当表层土归一化深度H＜1.0时,空沟归一化深度T应大于1.3,而1.0≤H≤2.0时,T应当大于1.8.     

三维层状地基空沟主动隔振分析

基于薄层法在研究层状介质中波的传播问题的高效性、边界单元法处理无限域问题的精确性,结合二者的优点提出三维层状地基薄层法基本解答,建立了基于三维层状半空间薄层法位移基本解答的半解析动力边界元法。该方法可有效的分析多层场地的动力问题,解决土-结构动力相互作用问题。同时分别对粘弹性上软下硬地基及上硬下软地基的三维空沟主动隔振进行了详细的分析。结果表明,两种情况下采用空沟屏障隔振均可以取得一定的隔振效果;同时,地基分层参数对空沟隔振体系的隔振效果影响显著。     

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

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

比例边界有限元子结构法研究

比例边界有限元法最初应用于土-结构的相互作用分析，经过近几年的完善和发展，如今已经能够应用到其他很多领域。但是因为比例边界有限元理论是基于相似性要求的，使得其在处理几何形状复杂的结构时，会有很大的局限性，从而在某些领域的应用仍旧受到限制。同时由于其全时空耦合，导致大量计算量和工作量，也是其应用受限的一个原因。采用子结构法，打破这些局限性，并且分别针对有限域、无限域的问题，对比例边界有限元子结构法进行了研究，得出了有利于比例边界有限元法在工程实践中应用的结论，为其在实际工程应用中提供了可靠的依据和规律。     

Elastodynamic boundary element formulation employing discontinuous curved elements

A boundary element formulation having discontinuous curved quadratic elements is presented for 2D elastodynamics. The first fundamental solution for static case is subtracted from and added to the first fundamental solution for dynamic case. As both kernels have the same order of singularity, the integral involving the regular expression arising from the subtraction can be calculated. matrix is calculated by employing the well-known rigid-body motion technique. The formulation is performed in Fourier transform space. Based on the formulation presented in this study, a general purpose computer program is developed for elastic or visco-elastic 2D elastodynamic problems. The program performs the analysis in Fourier transform space and can also be used for static analysis by assigning a very small value close to zero for the frequency. The results of some elastodynamic and dynamic soil–structure interaction problems obtained using the present study are compared with those in the literature.     

Acoustic viscoelastic modeling by frequency-domain boundary element method

Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green’s function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.     

有限元模拟中边界条件对计算结果的影响

本文从总电位法和异常电位法两个方面,分别就直流电法有限元中所采用的三类边界条件对计算结果的影响进行了分析,为采用合适的边界条件及选取最佳边界距离来提高正演精度提供了基础.结果表明,混合边界条件精度最高,可大大缩小求解区域而不影响计算精度,其次是Dirichlet边界条件精度较高,但测点越靠边界误差会越大,必须取足够大的边界区域,齐次边界条件的误差比较大,但如果采用非二极装置,通过电位差计算得到的视电阻率,由于无穷远边界对电位差的影响基本消除,视电阻率计算误差与混合边界条件下的接近,在反演中,为了节省计算时间,经常使用齐次边界条件进行有限元正演.     

基于MLS形函数处理边界的电阻率法有限单元正演模拟

电阻率法有限单元正演模拟中,采用第三类边界条件时为保证精度仍要求较大范围的计算域.无单元法为地球物理领域的新兴正演模拟方法,其计算效率低,但其中采用的移动最小二乘（MLS）形函数相比于有限单元法形函数具有良好的连续性,模拟精度高.本文将MLS形函数应用于电阻率法有限单元2.5维正演的第三类边界条件处理,提出电阻率法有限单元-移动最小二乘（FEM-MLS）耦合正演方法.通过不同正演方法的模型算例模拟结果对比,验证了本文算法的有效性,并讨论了各个参数选择对模拟结果的影响.本文数值模拟结果表明采用第三类边界条件时,在同等计算精度前提下,FEM-MLS耦合法相比于有限单元法可进一步缩小计算域并提高了计算效率,相比于采用较大计算域满足边界条件的有限单元法计算效率提高了约一倍,相比于采用相同小范围计算域的有限单元法平均精度提高了约一倍.

     

Rayleigh波入射下相邻山体地震响应边界元法模拟

采用间接边界元法（IBEM）,对Rayleigh波入射下两邻近山体的地震响应进行了定量分析。结果表明:与单个山体在Rayleigh波入射下的地震响应相比,两山地形的地震反应规律更为复杂,反应特征受控于入射波频率和山体间距等因素。总体上看:入射Rayleigh波频率较低时,两山间距对双山地形地震反应影响较大。低频波入射,受邻近山体影响,迎波面山体水平位移峰值约为入射波水平位移的6.3倍,放大效应可达单山作用的1.5倍。且山体山脚处竖向位移的频谱振荡更为剧烈,特定频率下山脚的竖向位移反应可达单山放大效应的1.6倍。较高频波入射时,迎波面一侧山体受邻近山体影响较小,且对波表现出明显的屏障效应,背波面山体地震反应强度被削弱。     

基于混合边界条件的有限单元法GPR正演模拟

从Maxwell方程组出发,推导了探地雷达(GPR)有限元波动方程.阐述了透射边界条件和Sarma边界条件的原理,推导了这两种边界条件的理论公式;通过在衰减层内加入过渡带优化了Sarma边界条件的加载方法,压制了介质区和衰减层交界面处的人为反射.考虑到透射边界条件与Sarma边界条件不同的理论机制,提出了一种结合透射边界条件和Sarma边界条件的混合边界条件,它利用Sarma边界条件对到达边界区域的GPR波能量衰减功能和透射边界对GPR波能量的透射功能,使GPR波经过Sarma边界条件的衰减吸收后,再通过透射边界条件将剩余能量透射出去,集成了二者的优势.并以二维均匀模型中的中心脉冲激励源方式为例,通过Matlab程序实现,以GPR的全波场快照的直观方式,对比了有、无边界条件及不同边界条件对人工截断边界的处理效果,说明了该混合边界条件对到达截断边界处的GPR波的处理优于单一边界条件.最后,以基于混合边界条件的有限单元法对两个典型的GPR地电模型进行了正演模拟,指导了GPR数据处理与工程实践.     

多源条件下直流电阻率法有限元三维数值模拟中一种近似边界条件

在多源直流电阻率法有限元三维数值模拟中,传统混合边界条件由于刚度矩阵与源点位置相关,无法实现线性方程组的快速回代求解,目前常用齐次边界条件或无限元边界进行替代,虽然实现了快速回代求解,但同时也降低了数值模拟的精度.为了实现快速回代求解,并确保数值模拟的计算精度,本文提出了一种近似边界条件方法.其核心思想是将与场源位置相关的边界系数矩阵从刚度矩阵中分离出来,使得分离后的刚度矩阵与场源位置无关.并将边界系数矩阵与边界处一次电场向量的乘积移到线性方程组右端源项中,当场源位置发生改变时,只需要重新计算右端源项就可以实现快速回代求解.理论模型数值计算表明,在水平地形条件下,本文边界条件数值精度优于混合边界条件;在起伏地形条件下,与齐次边界条件相比,本文边界条件数值结果与混合边界条件吻合度更高.     

基于比例边界有限元方法的拱坝子结构分区形式研究

比例边界有限元方法是一种半解析的数值计算方法,具有降维、网格灵活、严格模拟无限域和无需基本解等特点。比例边界有限元方法的基本理论是在整体坐标与局部坐标的比例边界转换基础之上建立的,相似中心的选取是否合理对分析计算具有重要的影响,导致在模拟拱坝这种不规则的空间壳体结构时,具有一定的局限性。采用子结构方法,将坝体分为若干满足相似性要求的区域可解决上述问题,以某拱坝为例给出了合理的坝体子结构分区形式,验证了子结构方法的精确性,为建立基于比例边界有限元方法的坝体-库水-地基系统的计算模型奠定了基础。     

单元畸变对动力有限元方法计算精度的影响

在有限元分析中,由于求解区域或边界的不规则,有限单元的划分会产生畸变单元。本文通过数值试验分析了有限单元畸变对动力有限元计算精度的影响,结果表明:长宽比畸变对动力有限元计算精度没有影响;当斜交角不小于30°时,斜交角对动力有限元计算精度的影响可忽略;锥度对动力有限元计算精度有影响,锥度越大,锥度对动力有限元计算精度的影响越大。     

自适应非结构有限元MT二维起伏地形正反演研究

在山区进行MT勘探时,用规则网格有限元方法模拟起伏地形会受到限制.本文采用非结构三角网格可以有效地模拟任意二维地质结构,如起伏地形、倾斜岩层和多尺度构造等.正演引入自适应有限元方法,其在网格剖分过程中能根据单元误差自动细化网格,保证了正演结果的精度.将自适应有限元与Occam算法结合,且引用并行处理技术提高正反演计算速度.通过对比两个理论模型,讨论了地形对MT正演响应的影响;其次进行了不同地电模型带地形反演展示了本文算法的正确性和适用性;最后将该方法应用于实测MT数据处理,证明了自适应非结构有限元方法是复杂地形下处理MT数据的有力工具.     

复电阻率法二维有限元数值模拟

伴随着复电阻率法的广泛应用,发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度,对无穷远边界进行了近似处理.整个正演计算过程分为两步,首先采用有限单元法计算四个不同频率的视复电阻率数据,然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比,验证了本文方法的正确性.设计了两个二维极化模型,数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映.     

Uniform half-plane elastodynamic problems by an approximate boundary element method

In this paper the applicability of an approximate Boundary Element Method to uniform half-plane elastodynamic problems is investigated. This method employs the concept of images to construct approximate fundamental solutions for the half-plane and does not require any half-plane surface discretization. The method is formulated in the frequency domain for the case of harmonic disturbances or the Laplace transform domain for the case of transient disturbances. In the latter case a numerical inversion of the transformed solution is necessary to obtain the time domain response. The proposed method can be used as an alternative to boundary element methods that either utilize the infinite plane fundamental solution and thus require a half-plane surface discretization, or employ the exact half-plane fundamental solution, which even though leads to no surface discretization, is of a very lengthy and complicated form. Two characteristics numerical examples are used to illustrate the proposed method and study its advantages and disadvantages.     

三维CSAMT法非结构化网格有限元数值模拟

考虑到可控源音频大地电磁法（CSAMT）电偶极发射源与地下介质的三维结构特点,本文采用非结构化网格剖分技术,开展了三维CSAMT方法有限元数值模拟研究,将三维电磁场的背景场和异常场分别求解,避免了电偶极发射源的奇异性问题,并减小了计算区域.推导了三维异常电场遵循的有限元方程,加入散度条件进行约束以消除电场伪解;对非结构化网格单元采用高斯加权平均算法,得到了精度较高的异常磁场.针对层状介质模型,与积分方程法对比,验证了有限元算法的正确性;计算分析了典型三维地质模型的电磁响应,异常体反映明显.结果表明本文算法正确、可靠,适用于三维地质模型的CSAMT方法正反演研究.

     

Modelling of pile wave barriers by effective trenches and their screening effectiveness

The three-dimensional problem of isolation of vibration by a row of piles is studied numerically on the basis of a model replacing the row of piles by an effective trench in order to reduce the modelling complexity. The analysis is accomplished with the aid of an advanced frequency domain boundary element method, which is used for both the infilled trench and the soil medium in conjunction with a coupling procedure based on enforcement of equilibrium and compatibility at the trench–soil interface. Linear elastic or viscoelastic material behaviour is assumed for both the piles and the soil. The piles can be tubular or solid and have circular or square cross-section. The vibration source is a vertical force, harmonically varying with time, and the row of piles acts as a passive wave barrier. The effective trench model is constructed by invoking well known homogenization techniques used in the mechanics of fibre-reinforced composite materials, and its accuracy is compared against a rigorous boundary element analysis modelling each pile separately in full contact with the soil medium. On the basis of the effective trench model, the screening effectiveness of a row of piles is studied through parametric studies.     

Modeling granular soils liquefaction using coupled lattice Boltzmann method and discrete element method

In this paper, a novel coupled pore-scale model of pore-fluid interacting with discrete particles is presented for modeling liquefaction of saturated granular soil. A microscale idealization of the solid phase is achieved using the discrete element method (DEM) while the fluid phase is modeled at a pore-scale using the lattice Boltzmann method (LBM). The fluid forces applied on the particles are calculated based on the momentum exchange between the fluid and particles. The presented model is based on a first principles formulation in which pore-pressure develops due to actual changes in pore space as particles׳ rearrangement occurs during shaking. The proposed approach is used to model the response of a saturated soil deposit subjected to low and large amplitude seismic excitations. Results of conducted simulations show that at low amplitude shaking, the input motion propagates following the theory of wave propagation in elastic solids. The deposit response to the strong input motion indicates that liquefaction took place and it was due to reduction in void space during shaking that led to buildup in pore-fluid pressure. Soil liquefaction was associated with soil stiffness degradation and significant loss of interparticle contacts. Simulation results also indicate that the level of shaking-induced shear strains and associated volumetric strains play a major role in the onset of liquefaction and the rate of pore-pressure buildup.