A 2.5-D dynamic model for a saturated porous medium. Part II: Boundary element method

The 3-D boundary integral equation is derived in terms of the reciprocal work theorem and used along with the 2.5-D Green’s function developed in Part I [Lu, J.F., Jeng, D.S., Williams, S., submitted for publication. A 2.5-D dynamic model for a saturated porous medium: Part I. Green’s function. Int. J. Solids Struct.] to develop the 2.5-D boundary integral equation for a saturated porous medium. The 2.5-D boundary integral equations for the wave scattering problem and the moving load problem are established. The Cauchy type singularity of the 2.5-D boundary integral equation is eliminated through introduction of an auxiliary problem and the treatment of the weakly singular kernel is also addressed. Discretisation of the 2.5-D boundary integral equation is achieved using boundary iso-parametric elements. The discrete wavenumber domain solution is obtained via the 2.5-D boundary element method, and the space domain solution is recovered using the inverse Fourier transform. To validate the new methodology, numerical results of this paper are compared with those obtained using an analytical approach; also, some numerical results and corresponding analysis are presented.     

Elastic wave scattering from a penny-shaped interface crack in coated materials

This paper is concerned with the elastic wave scattering induced by a penny-shaped interface crack in coated materials. Using the integral transform, the problem of wave scattering is reduced to a set of singular integral equations in matrix form. The singular integral equations are solved by the asymptotic analysis and contour integral technique, and the expressions for the stress and displacement as well as the dynamic stress intensity factors (SIFs) are obtained. Using numerical analysis, this approach is verified by the finite element method (FEM), and the numerical results agree well with the theoretical results. For various crack sizes and material combinations, the relations between the SIFs and the incident frequency are analyzed, and the amplitudes of the crack opening displacements (CODs) are plotted versus incident wavenumber. The investigation provides a theoretical basis for the dynamic failure analysis and nondestructive evaluation of coated materials.     

Multiple reciprocity boundary element analysis of two-dimensional anisotropic thermoelasticity involving an internal arbitrary non-uniform volume heat source

In the direct boundary element method (BEM) formulation of anisotropic thermoelasticity, thermal loads manifest themselves as additional volume integral terms in the boundary integral equation (BIE). Conventionally, this requires internal cell discretisation throughout the whole domain. In this paper, the multiple reciprocity method in BEM analysis is employed to treat the general 2D thermoelasticity problem when the thermal loading is due to an internal non-uniform volume heat source. By successively performing the “volume-to-surface” integral transformation, the general formulation of the associated BIE for the problem is derived. The successful implementation of such a scheme is illustrated by three numerical examples.     

非均匀介质散射问题的体积分方程数值解法

将非均匀介质视为某一均匀背景介质中的扰动，可建立用均匀背景介质格林函数作基本解的体积分方程．给出了配置法求解体积分方程的数值方法，首先解得扰动域内各点以速度扰动为权的波场函数，然后回代计算得到观测面上各接收点的散射波场．与边界元法和Ｂｏｒｎ近似法计算结果比较表明该方法具有很高的精度，可得到穿过非     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程，结合广义Betti-Rayleigh互易等式与时域下的边界积分方程，推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解，将经过主部分析的积分核函数分解为静态和动态部分，其中动态积分核不具有奇异性。在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换，采用配置点法计算超奇异积分，获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序，得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律，数值结果稳定且收敛速度快。     

Solution of multiple crack problem in a finite plate using coupled integral equations

This paper studies a numerical solution of multiple crack problem in a finite plate using coupled integral equations. After using the principle of superposition, the multiple crack problem in a finite plate can be converted into two problems: (a) the multiple crack problem in an infinite plate and (b) a usual boundary value problem for the finite plate. For the former problem, the Fredholm integral equation is used. For the latter problem, a BIE based on complex variable is suggested in which a Cauchy singular kernel exists. For the proposed BIE, after using the inverse matrix technique, the dependence of the traction at a domain point from the boundary tractions is formulated indirectly. This is a particular advantage of the present study. Several numerical examples are provided and the computed results for stress intensity factor and T-stress at crack tips are given.     

Finite element-boundary integral formulation for electromagnetic scattering

A new formulation is described which combines the most robust attributes of the volume finite element and surface integral equation approaches to electromagnetic boundary value solutions. The result is a numerical technique which may be applied to scattering problems involving configurations having metallic surfaces and inhomogeneous penetrable material situated in open spatial regions. This is accomplished by way of coupling internal region finite element modal field solutions to equivalent currents on the surrounding boundary surface through an appropriate surface integral equation. The method is demonstrated for the special case of scattering by axisymmetric inhomogeneous penetrable objects. Example numerical calculations are presented for validation of the procedure and potential problem areas are discussed.     

The Gelfand-Levitan, the Marchenko, and the Gopinath-Sondhi integral equations of inverse scattering theory, regarded in the context of inverse impulse-response problems

In this paper, we are concerned with interpreting some classical inverse-scattering theories so that they are relevant to the one-dimensional inverse impulse-response problem which arises in reflection seismology

First the Gelfand-Levitan integral equations (which arise in the inverse scattering theory for the Schrödinger equation) are derived strictly in the time domain. Originally these celebrated equations were derived as a means of solving an inverse spectral problem, which is naturally posed in the frequency domain.

We show not only that these equations have a time-dependent formulation, but that their derivation is actually simpler here than in its original context. Next we give a similar derivation for the Marchenko integral equation.

We then obtain, by an independent method, the integral equation of Gopinath and Sondhi, which is not unlike the Gelfand-Levitan linear equation. It was used by them as a means of solving a time-dependent inverse problem arising in speech synthesis. A new integral equation, similar to the Marchenko equation is also derived.

Finally the integral equations are related to each other by direct transformation independent of the corresponding differential equations.

The present paper opens with a section in which the equation of one-dimensional elastic waves and the corresponding seismic inverse impulse-response problem are transformed into a form to which the Gelfand-Levitan theory applies, and then into the equations which arose in Gopinath and Sondhi's work.

We regard the integral equation of Gopinath and Sondhi as being directly applicable to the interpretation of seismic reflection data. The Gelfand-Levitan theory is also applicable but only after considerable transformation.

Our calculations are entirely in the time domain. The resulting equations have the merit that in order to recover the unknown coefficient on a finite interval (0, L) we need to use the boundary data also only on a finite segment of the reflection time series. The length of the record is just the two-way travel time associated with length L. By contrast, frequency-domain approaches require that long records be Fourier transformed. We distinguish results which are true only when the unknown coefficient is continuous from those true when it is bounded but only piecewise continuous.

In this work we have not addressed the important problem of deconvolution.     

二维介质参数的大扰动反演方法

对非均匀介质参数反演问题进行了研究，并提出了用于反演二维介质参数的广义射线近似方法．利用参考场量和扰动变量对声波方程中的介质参数进行处理，并利用Ｇｒｅｅｎ函数理论得到扰动参数比的积分方程．基于非均匀介质中波函数的局部理论和射线理论，引入了全波场的广义射线近似形式，通过定义介质参数函数，把反演目标归结为其第一类Ｆｒｅｄｈｏｌｍ积分方程．利用积分变换方法得到二维介质的介质参数函数，从而得到介质参数，在Ｂｏｒｎ近似方法中，反演的介质参数扰动不能超过２０％，但是在本文中介绍的方法能够有效地反演其扰动比不超过５０％的变化情况     

High‐order filtering for control volume flow simulation

A general methodology is presented in order to obtain a hierarchy of high‐order filter functions, starting from the standard top‐hat filter, naturally linked to control volumes flow simulations. The goal is to have a new filtered variable better represented in its high resolved wavenumber components by using a suitable deconvolution. The proposed formulation is applied to the integral momentum equation, that is the evolution equation for the top‐hat filtered variable, by performing a spatial reconstruction based on the approximate inversion of the averaging operator. A theoretical analysis for the Burgers' model equation is presented, demonstrating that the local de‐averaging is an effective tool to obtain a higher‐order accuracy. It is also shown that the subgrid‐scale term, to be modeled in the deconvolved balance equation, has a smaller absolute importance in the resolved wavenumber range for increasing deconvolution order. A numerical analysis of the procedure is presented, based on high‐order upwind and central fluxes reconstruction, leading to congruent control volume schemes. Finally, the features of the present high‐order conservative formulation are tested in the numerical simulation of a sample turbulent flow: the flow behind a backward‐facing step. Copyright © 2001 John Wiley & Sons, Ltd.     

求解三维裂纹前缘SIF分布时间历程的通用权函数法和有限变分法

将三维热权函数法扩展为适用于表面力、体积力和温度载荷的通用权函数法(UWF).推导出以变分型积分方程表达的UWF法基本方程,从变分的角度,将求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE)改造为可以适用于一般的变分型积分方程的一类新型数值方法--有限变分法(FVM).在FVM中可以引入无穷多种线性无关的局部变分模式,可以根据计算要求在求解域中插入任意多个计算节点,单一型裂纹问题FVM所得到的最终方程组的系数矩阵总是一个对称的窄带矩阵,而且对角元总是大数,具有良好的数值计算性能.FVM对于SIF沿裂纹前缘急剧变化的复杂情况具有较好的数值模拟能力和较高的计算精度,利用自身一致性,可以求得三维裂纹前缘SIF的高精度解.     

Numerical analysis of electromagnetic wave scattering and interaction based on frequency-domain integral equation and method of moments techniques

This paper discusses numerical modeling of electromagnetic wave scattering and interaction by general arbitrary-shaped two- and three-dimensional material objects in free space based on the frequency-domain integral equation method. Depending on the type of material object modeled, straightforward electromagnetic boundary value problems can be set up in terms of a set of coupled integro-differential equations with unknown surface equivalent current distributions. These unknown surface equivalent currents are determined by a numerical technique based on the method of moments (MOM) which involves reducing the exact coupled integro-differential equations to a corresponding partitioned dense matrix equation. Various electromagnetic fields (penetrating, near scattered, and far scattered) as well as the radar cross section can be calculated once the unknown surface current distributions are determined.

Canonical two- and three-dimensional scattering and interaction examples are reported which illustrate the capabilities of this method for analyzing arbitrary-shaped objects comprised of conducting, dielectric, and anisotropic materials. For specific geometries, both the induced surface current distribution and the corresponding radar cross section are presented with either analytical or code-to-code validations. The paper concludes with a discussion of means to extend the applicability of the integral equation modeling technique to the case of electrically-large material objects.     

A wideband fast multipole boundary element method for half-space/plane-symmetric acoustic wave problems

This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.     

Dynamic stress concentration for multiple multilayered inclusions embedded in an elastic half-space subjected to SH-waves

The dynamic stress concentration factor (DSCF) is evaluated along the interfaces of multiple multilayered inclusions embedded in a half-space when subjected to a plane harmonic SH-wave. A weak form of Helmholtz equation is utilized to derive a non-hypersingular boundary integral equations to compute the stresses. Eliminating the need to rely on hypersingular integrals, greatly simplifies the procedure. The numerical results obtained by the proposed method, are validated against analytical solutions.Various contributing factors that can influence the DSCF are investigated, including multiple scattering, layering, stiffness of the adjacent inclusions, and impedance contrast of the layers. The DSCF is found to be highly prone to these changes, particularly with the soft materials. Therefore, accurate analysis of stresses requires models that consider multiple scattering and layering. The presented result could be used for predicting the seismic failure of pipes and underground tunnels and for estimating the stress failure in strong ground motion seismology due to subsurface irregularities.     

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

A HIGH ORDER KERNEL INDEPENDENT FAST MULTIPOLE BOUNDARY ELEMENT METHOD FOR ELASTODYNAMICS

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nyström 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nyström 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nyström 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.     

弹性动力学高阶核无关快速多极边界元法

基于核无关的快速多极方法, 发展了一种弹性动力学问题的快速、高精度边界元分析方法. 采用基于二次曲面单元的Nystr?m 离散, 将边界积分方程转化为求和形式, 可以方便地进行加速计算;由于采用二次元, 边界元分析精度很高. 将一种新型快速多极方法用于Nystr?m 边界元法的加速计算, 该方法的数值实现简便、不依赖于积分方程基本解的表达式, 因此通用性很好;该方法还具有最优的计算量和存储量、精度高且可以控制. 结合Nystr?m 边界元系数矩阵和快速多极方法转换矩阵的特点, 提出一种大幅度降低边界元内存消耗的策略. 数值结果表明, 该方法无论在分析精度, 还是计算速度和内存消耗上, 都大大优于同类方法, 是一种快速、通用的工程弹性动力学问题大规模数值分析方法.      

A coupling method of difference with high order accuracy and boundary integral equation for evolutionary equation and its error estimates

In the present paper,a new numerical method for solving initial-boundary valueproblems of evolutionary equations is proposed and studied,combining difference methodwith high accuracy with boundary integral equation method.The numerical approximateschemes for both problems on a bounded or unbounded domain in R~3 are proposed and theirprior error estimates are obtained.     

大延伸非均匀介质中地震波全弹性散射理论II-弹性波多次散射理论

推出大延伸非均匀连续介质的弹性波能量传输表达式,并以此为基础建立并发展了相应的弹性动力学能量传递理论.构造了一个基于非均匀薄层或非均匀相屏单次散射迭代法的多次散射模型.该模型既适用于弱散射,也适用于强散驰 既适用于普通散射,也适用于转换散射;尽管高频情况下只考虑普通散射及前向散射.应用该模型计算了弹性多次散射的能通量,处理了散射衰减问题。数值实验的结果表明,短周期地震图上的尾波主要来源于S波散射。