基于高阶双渐近透射边界的重力坝-层状地基动力相互作用分析

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

     

墩-水耦合计算模式及深水桥墩动力响应研究

为了研究地震作用下动水压力对深水桥墩的影响,以声波动理论为基础,通过对流体域边界条件的简化和流体域范围的探讨,建立了一种改进的墩-水耦合有限元模型,并提出了常用深水桥墩的理想流体域范围图。分别讨论了深水桥墩在简谐荷载、地震荷载作用下的动力响应,结果表明:动水压力增大了桥墩的动力响应,增幅在20%到50%之间,不可忽视;当桥墩基频越接近简谐荷载激励频率或地震波主频有效带宽区间时,动水压力的影响越大。     

考虑水流-结构强耦合作用的水流动力效应分析

在深水桥墩及桩基础等结构地震反应中,地震激励下的水流动水压力对水中结构的作用以及结构位形变化对水体的反作用,属于典型的流固强界面耦合问题。以一顶端伸出水面的圆柱式结构为研究对象,基于任意拉格朗日-欧拉描述的Navier-stokes方程,建立了考虑水流-结构强耦合效应的水流-结构三维有限元模型。以正弦位移波输入,考察结构材料模量、水流流速和水位、激振频率和位移幅值等多种因素,分析了结构表面作用的动水压力反应特征以及水流动力效应,探讨了水流动力效应的主要影响因素。结果表明：考虑流固强耦合作用时,结构表面作用的动水压力及其分布具有强烈的频率依赖性,高频激励可显著增强动水压力作用;由于结构周围流体具有一定的粘滞性及动载作用下具有较强的辐射阻尼效应,水流动力效应对结构的位移、内力反应均有一定的抑制作用;激振频率、水流流速和水位以及结构材料模量等因素,对水流动力效应均有一定的影响。     

Two-dimensional development of the dynamic coupled consolidation scaled boundary finite-element method for fully saturated soils

The scaled boundary finite-element method is a powerful tool used to analyse far-field boundary soil–structure interaction problems. In this paper, the method is extended to include Biot's coupled consolidation in order to deal with fully saturated soil as a two-phase medium. The advantages of this method are explained in this paper. The detailed formulation considers the general two-dimensional (2D) analysis case, accounting for body forces and surface tractions in both bounded and unbounded media.     

基于吸收边界条件的瞬变电磁法三维矢量有限元快速正演

瞬变电磁法的三维有限元正演通常采用齐次边界条件,为满足该边界条件,需要构建较大尺寸的模型,这降低了正演问题的求解速度.针对该问题,本文采用吸收边界条件代替齐次边界条件,以缩小模型体积,加快正演速度：首先,从时间域麦克斯韦方程组出发,推导了基于库伦规范的矢量势的微分控制方程,结合一阶吸收边界条件推导了相应的的弱形式方程;在此基础上采用一阶四面体矢量单元进行单元分析、Newmark法进行时间离散,实现了瞬变电磁法的快速三维正演.通过均匀半空间模型的解析解,H型地电断面的CR1Dmod解和相应模型有限元解的对比,验证了本文算法的正确性.均匀半空间模型分别采用吸收边界条件和齐次边界条件的正演结果对比表明：吸收边界条件确实可以提高三维正演的精度或者缩小模型尺寸、加快计算速度.

     

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

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

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography.The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.     

A storm event watershed model for surface runoff based on 2D fully dynamic wave equations

The main purpose of this work concerns the development and testing of an overland flow model based on the two‐dimensional fully dynamic shallow water equations. Three key aspects, fundamental to get accurate, efficient and robust computation of surface runoff at basin scale, are discussed by transferring the main findings obtained by the recent research on the topic of dam‐break wave and flood propagation in the context of rainfall–runoff modelling. In particular, attention is focused on the numerical flux and bottom slope source terms computation, on a numerical treatment of friction slope terms and on an algorithm for dealing with wetting/drying fronts. The performances of the numerical model have been preliminarily evaluated using experimental or ideal tests characterized by very critical conditions for the stability of a numerical model. Then, attention was focused on a real event occurred in a sub‐basin of Reno river in Italy to analyse the suitability of the model in simulating real flood situations. The numerical results highlight the good performances of the model in all the simulations discussed in the paper. Copyright © 2012 John Wiley & Sons, Ltd.