用粘性边界有限元法分析弹性半无限地基中的动力问题

本文使用设置粘性边界单元的有限元方法，分析了简谐集中激振力产生的地表位移反应，刚性基础及桩基础的阻抗函数等半无限地基中的动力问题。计算结果与其他数值分析方法结果的比较表明，粘性边界单元的有限元方法适用于分析弹性半无限地基中的动力问题。本文还讨论了有限元网格尺寸及模型大小对计算结果的影响。     

结构-地基系统静-动力联合分析模型

以有限区域模拟无限区域的动力计算必须考虑边界上波的传播效应,应不使波在边界上发生反射再返回到计算域中,并且应考虑初始静力场对其动力反应的影响。因此,基于无限元的静、动力特性提出静-动力统一无限元,通过算例验证无限元静-动力统一人工边界在静-动力联合分析模型中的精度和稳定性。结果表明:无限元静-动力统一人工边界在静力、动力及静-动力联合分析模型中具有很高的精度和良好的稳定性,突破了传统动力人工边界不能兼顾静力效应的局限。     

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

基于Biot流体饱和多孔介质本构方程，分别考察具有辐射阻尼性质的外行柱面波和球 面波在圆柱面和球面人工边界上引起的法向、切向应力的表达式. 在应力表达形式上，固相 介质和孔隙流体的法向和切向应力都是由两项组成，它们分别与质点的位移和速度成正比， 因此，可在人工边界的法向和切向设置连续分布的并联弹簧------黏滞阻尼器，用来模拟人工边 界以外的无限域介质对来自有限计算域的外行波动的能量吸收作用，从而形成了流体饱和多 孔介质的黏弹性动力人工边界. 流体饱和多孔介质的黏弹性动力人工边界可方便地与大型通 用软件结合，用于分析饱和土中复杂的结构-地基动力相互作用问题. 算例表明流体饱和多 孔介质黏弹性动力人工边界具有较好的精度和稳定性.     

成层半空间出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性水平成层半空间中自由波场时域计算的一维化有 限元方法. 在进行有限元网格划分时，竖向单元取满足有限元模拟精度的任意尺寸，水平向 网格尺寸由时间离散步长和水平视波速确定，并自动进行虚拟网格划分. 基底设置人工边界， 并将波动输入转化为等效荷载施加在边界节点上. 然后将集中质量有限元法和中心差分法相 结合建立节点运动方程，并将水平方向相邻节点的运动用该节点相邻时刻的运动表示，从而 将求解节点运动的二维方程组转化为一维方程组. 求解此方程组，即得到自由场中竖向一列 节点的运动. 最后根据行波传播的特点，可方便地确定全部自由波场. 理论分析和数值算例 表明，该方法具有较高的精度和良好的稳定性.     

黏弹性人工边界等效荷载计算的改进方法

黏弹性人工边界在场地地震反应和结构-地基动力相互作用等问题的计算中已得到了广泛的应用.地震波在黏弹性人工边界中的输入是通过将地震波转化为作用于人工边界处的等效载荷来实现的.计算等效节点载荷的常规方法默认边界节点对应区域内的应力为均布力,但实际上该节点对应区域内的应力分布通常是不均匀的.本文在有限元方法结合黏弹性局部人工边界的显式时域波动方法的基础上,建立了无限域散射问题地震波等效载荷计算的一种改进方法.该方法采用细化网格与应力积分相结合的方法计算人工边界等效节点力,有效地降低了人工边界上等效节点力的计算误差.以不同角度入射地震波的二维算例为例,算例给出的波场位移云图和节点位移时程曲线验证了本文方法的有效性,其计算精度与网格尺寸和地震波入射角度密切相关,且网格越小、入射角度越小,计算精度越高.对于相同的网格尺寸,本文采用方法的计算精度明显高于常规方法,尤其是对于斜入射问题优势更为明显.     

三维裂纹分析的主值型面力边界积分方程法

给出了一组只包含Ｃａｕｃｈｙ主值积分、不含有强奇异积分的三维静动力边界积分方程及其应用于裂纹问题的具体列式，并给出了几何轴对称问题的相应半解析边界元求解方法，将三维问题降阶为一维数值问题．文中分析了无限、半无限介质中圆裂纹、平行圆裂纹系、球面裂纹等在静载及应力波作用下的静力或瞬态动力响应问题，求得了相应的应力强度因子．     

数值流形方法的粘性边界问题初探

在实际工程数值流形方法分析中,采用固定约束边界的方法处理无限域或者半无限域的情况,边界处应力波的反射造成模拟结果与实际情况不符.本文基于Lysmer等人提出的粘性边界理论,在边界上设置阻尼器,推导相应粘性边界条件下流形单元刚度矩阵的数值计算格式,经岩石长条中弹性波传播算例,并与有限元结果对比,验证了该粘性边界的有效性,有利于数值流形方法的工程中推广应用.     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

动态断裂力学的无限相似边界元法

对弹性动力学的相似边界元法进行了进一步研究,推导了相应的计算公式,并在此基础上提出了动态断裂力学的无限相似边界元法.与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,大大减少了计算量.对动态断裂力学问题,无限相似边界元法由于在裂纹尖端的边界上设置了逼近于裂纹尖端的无限个相似边界单元,可直接得到裂纹尖端具有奇异性的应力,而不需要设置奇异单元,从而突破了奇异单元对应力奇异性阶次的局限.另外,还讨论了无限相似边界元法得到的无限阶的线性代数方程组的求解方法.     

反平面断裂问题的无单元伽辽金比例边界法

将比例边界法与无单元伽辽金法相结合,建立了反平面断裂分析的无单元伽辽金比例边界法。这是一种边界型无网格法,在环向方向上采用无单元伽辽金法进行离散,因此计算时仅需要边界上的节点信息,不需要边界元所要求的基本解。为了便于施加本质边界条件,通过建立节点值和虚拟节点值之间的关系给出了修正的移动最小二乘形函数。在径向方向上,该方法利用解析的方法求解,因此是一种半解析的数值方法。最后,给出了数值算例,并验证了所提方法后处理简单和计算精度高的特点,适合于求解反平面断裂问题。     

A HYBRID METHOD FOR ANALYZING THE DYNAMIC RESPONSES OF CAVITIES OR SHELLS BURIED IN AN ELASTIC HALF-PLANE

In this paper, based on a variational formalism which originally proposed by Mei [1] for infinite elastic medium and extended by Yeh, et al. [2,3] for elastic half-plane, a hybrid method which combines the finite element and series expansion method is implemented to solve the diffraction of plane waves by a cavity buried in an elastic half-plane. The finite domain which encloses all inhomogeneities including the cavity can be easily formulated by finite element methods. The unknown boundary data obtained by subtracting the known free fields from the total fields which include the boundary nodal displacements and tractions at the interface between the finite domain and the surrounding elastic half-plane are not independent of each other and can be correlated through aseries repre sentation. Due to the continuity condition at the interface, the same series representation is still valid for the exterior elastic half-plane to represents the scattered wave. The unknown coefficients of this series are treated as generalized coordinates and can be easily formulated by the same variational principle. The expansion function of the series is composed of basis function. Each basis function is constructed from the basis function for an infinite plane by superimposing an additional homogeneous reflective term to satisfy both traction free conditions at ground surface and radiation conditions at infinity. The numerical results are made against those obtained by boundary element methods, and good agreements are found.     

Linear viscoelastic analysis of a semi-infinite porous medium

This paper considers the problem of a semi-infinite, isotropic, linear viscoelastic half-plane containing multiple, non-overlapping circular holes. The sizes and the locations of the holes are arbitrary. Constant or time dependent far-field stress acts parallel to the boundary of the half-plane and the boundaries of the holes are subjected to uniform pressure. Three types of loading conditions are assumed at the boundary of the half-plane: a point force, a force uniformly distributed over a segment, a force uniformly distributed over the whole boundary of the half-plane. The solution of the problem is based on the use of the correspondence principle. The direct boundary integral method is applied to obtain the governing equation in the Laplace domain. The unknown transformed displacements on the boundaries of the holes are approximated by a truncated complex Fourier series. A system of linear equations is obtained by using a Taylor series expansion. The viscoelastic stresses and displacements at any point of the half-plane are found by using the viscoelastic analogs of Kolosov–Muskhelishvili’s potentials. The solution in the time domain is obtained by the application of the inverse Laplace transform. All the operations of space integration, the Laplace transform and its inversion are performed analytically. The method described in the paper allows one to adopt a variety of viscoelastic models. For the sake of illustration only one model in which the material responds as the standard solid in shear and elastically in bulk is considered. The accuracy and efficiency of the method are demonstrated by the comparison of selected results with the solutions obtained by using finite element software ANSYS.     

弹性波与单侧界裂纹相互作用问题的边界元法

措助边界元法设计了一种迭代修正方法来求解单侧界面裂纹模型与弹性波的相互作用问题,作为对算法的检验,用这种方法我们具体地分析了平面简谐弹性波对一个则界面裂纹的入射,给出了裂纹面的接触形态及应力场。     

Diffraction of acoustic and elastic waves on a half-plane for boundary conditions of various types

The classical problem of wave diffraction on a half-plane with boundary conditions of different types and its generalizations to elastic media are considered. As a solution method it is proposed to combine the Fourier method of separation of variables and the series summation technique based on the use of integral representations of Bessel functions. The analytic solutions thus obtained are equally efficient in the near- and far-field diffraction regions. The two-term singularity at a corner point (in stresses for elastic media and in the velocity for acoustic media) was discovered for the first time. The knowledge of singularities in the scalar problem allowed one to construct the solution of the vector problem of elastic longitudinal wave diffraction. It is investigated how different types of boundary conditions on both sides of the half-plane affect the solution behavior in the far-field region. Possible physical interpretations of the obtained results are given.     

Dynamic response of underground structures by time domain SBEM and SFEM

The dynamic interaction problems of three-dimensional linear elastic structures witharbitrary shaped section embedded in a homogeneous,isotropic and linear elastic half spaceunder dynamic disturbances are numerically solved.The numerical method employed is acombination of the time domain semi-analytical boundary element method(SBEM)usedfor the semi-infinite soil medium and the semi-analytical finite element method(SFEM)used for the three-dimensional structure.The two methods are combined throughequilibrium and compatibility conditions at the soil-structure interface.Displacements,velocities,accelerations and interaction forces at the interface between undergroundstructure and soil medium produced by the diffraction of wave by an underground structurefor every time step are obtained.In dynamic soil-structure interaction problems,it isadvantageous to combine the SBEM and the SFEM in an effort to produce an optimumnumerical hybrid scheme which is characterized by the main advantages of the two methods.The     

黏弹性多层介质中SH波动的一种吸收边界条件

提出一种高精度时域吸收边界条件,与有限元法结合用于模拟瞬态标量SH波在达朗贝尔黏弹性多层介质中传播问题.建立时域吸收边界条件的过程是:首先将半无限域沿着竖向半离散得到半离散的位移方程以及人工边界处的力-位移关系,再通过引入模态分解, 将物理空间下的量转化到模态空间,从而获得半无限域模态空间下的频域动力刚度;其次采用一种在全频范围内收敛的连分式精确逼近单层介质模态空间下标量形式的频域动力刚度,将标量连分式扩展为矩阵形式用来表示多层介质的频域动力刚度;最后通过引入辅助变量技术,将模态空间下基于连分式的频域动力刚度关系转化为时域吸收边界条件,进一步转换到物理空间后得到物理空间下的时域吸收边界条件.单层介质和五层介质的数值算例表明,建立的高精度时域吸收边界条件对于达朗贝尔黏弹性单层介质是精确且稳定的;对于达朗贝尔黏弹性多层介质, 为了保证其高精度特性,需要将人工边界放置在距离感兴趣区域约为0.5倍无限域高度的位置处.      

Diffraction of a shear plane wave in a compound elastic space with a half-infinite crack parallel to the inhomogeneity line

Problems of stress wave propagation and diffraction in elastic inhomogeneous media are undoubtedly of interest to scientists from the viewpoint of investigation of fundamental laws of dynamic processes and of the use of the results in technical and technological applications. The paper deals with the dynamic contact problem of shear plane wave diffraction at the edge of a semi-infinite crack in a compound space consisting of two elastic half-spaces. The questions related to the onset of surface waves and the wave field behavior in far-field regions are also considered.     

Mixed Boundary Value Problem in a Non-homogeneous Medium Under Steady Distribution of Temperature

This paper is concerned with a mixed boundary value problem of a non-homogeneous medium under steady distribution of temperature whose elastic constants are exponential functions of y. By using Fourier cosine transforms the mixed boundary value problem of heat conduction is reduced to a Fredholm integral equation of the second kind. Then the elastic problem of the non-homogeneous semi-infinite half-plane under distribution of load over a plane face is solved. The influence of the non-homogeneity of the material on the thermal stress distribution is illustrated graphically.     

多层波导中矢量波动的时域人工边界条件

本文提出了一种近似的时域人工边界条件(artificial boundary condition, ABC)用来模拟含有瑞利阻尼的线弹性多层波导平面内的矢量波动,该ABC是时域稳定的, 且能与有限元法无缝耦合. 建立ABC的思路是,首先将多层波导的矢量波动方程简化为$x$方向和$y$方向解耦的两个标量波动方程;其次基于比例边界有限元法得到无限域$x$方向和$y$方向模态空间下半离散的频域动力刚度,再用矩阵连分式近似表示$x$方向和$y$方向的频域动力刚度;最后通过辅助变量技术将连分式时域化,从而分别得到人工边界上$x$方向和$y$方向的时域ABC.方法中影响计算精度和计算效率的参数有无限域的模态数$n$、连分式阶数$J$和人工边界远离兴趣域的距离$L$. 数值算例表明,仅需将被载荷激起的无限域的模态数$n$参与计算, 一般可以取$J$=3,$L$的取值基本与地下结构尺寸无关, 它与土层的总层高$H$成正比关系,关系系数与土层的材料参数有关.      

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.