Direct and Inverse Problems in Elastic Multilayers with Reflection Data

The paper deals with the reflection-transmission process generated by a multilayer with homogeneous, isotropic elastic slabs. The direct problem, namely the determination of reflected and transmitted (and internal) waves in terms of the incident one, is investigated and the solution is proved to exist and be unique. The proof is based on the validity of a first integral. Next the inverse problem is considered to determine the material parameters of the slabs through the measure of the (frequency-dependent) reflection coefficient at the boundary of the multilayer. As a result, the material parameters are found in closed form, though for two slabs only.     

弹性动力学的双互易杂交边界点法

将双互易法同杂交边界点法相结合,提出了求解弹性动力问题的新型数值方法------双互易杂交边界点方法. 该算法在求解弹性动力问题时,将控制方程非齐次项的域内积分转化为边界积分. 该方法将问题的解分为通解和特解两部分,通解使用杂交边界点法求得,特解则使用局部径向基函数插值得到,从而实现了使用静力问题的基本解来求解动力问题. 计算时仅仅需要边界上离散点的信息,无论积分还是插值都不需要网格,域内节点仅用来插值非齐次项,因此该算法仍是一种边界类型的无网格方法. 数值算例表明,该方法后处理简单,计算精度高,适合于求解弹性动力问题.      

关于线性热释电弹性介质的互等功定理及应用

根据线弹性理论以及电介质理论推出了类似于Ｂｅｔｔｉ互等功定理的热电弹性耦合体─—热释电弹性体的互等功方程；作为一种应用推导了在边界积分方法中被广泛采纳的热电弹性Ｓｏｍｉｇｌｉａｎａ方程．     

MESHLESS METHOD OF DUAL RECIPROCITY HYBRID RADIAL BOUNDARY NODE METHOD FOR ELASTICITY

Combining the radial point interpolation method （RPIM）, the dual reciprocity method （DRM） and the hybrid boundary node method （HBNM）, a dual reciprocity hybrid radial boundary node method （DHRBNM） is proposed for linear elasticity. Compared to DHBNM, RPIM is exploited to replace the moving least square （MLS） in DHRBNM, and it gets rid of the deficiency of MLS approximation, in which shape functions lack the delta function property, the boundary condition can not be applied easily and directly and it＇s computational expense is high. Besides, different approximate functions are discussed in DRM to get the interpolation property, in which the accuracy and efficiency for different basis functions are compared. Then RPIM is also applied in DRM to replace the conical function interpolation, which can greatly improve the accuracy of the present method. To demonstrate the effectiveness of the present method, DHBNM is applied for comparison, and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.     

Dual reciprocity hybrid radial boundary node method for Winkler and Pasternak foundation thin plate

An efficient dual reciprocity hybrid radial boundary node method is developed for the analysis of Winkler and Pasternak foundation thin plate, in which a hybrid displacement variational principle, radial point interpolation method (RPIM) and dual reciprocity method (DRM) are combined. Firstly, the hybrid displacement variational principle is developed, in which the domain variables are interpolated by two groups of symmetric fundamental solutions, while the boundary variables are interpolated by RPIM instead of the traditional moving least square, and the shape function obtained by RPIM satisfies the delta function property, so boundary conditions can be applied directly. Besides, DRM is exploited to evaluate the particular solutions of inhomogeneous terms, which can be used to transform the domain integrals arising from the inhomogeneous term into equivalent boundary integrals. Finally, some additional equations based on the DRM theory are proposed to overcome the problem that the boundary integral equations are not enough to solve all variables. This method has the advantages of both no element mesh of meshless method and dimensionality reduction of boundary element method. Numerical examples of Winkler and Pasternak foundation plates are given to illustrate that the present method is effective, accurate and it can be further expanded into practical engineering.     

Wave propagation and reflection-transmission in a stratified viscoelastic solid

The paper investigates time-harmonic wave propagation in continuously stratified solids and provides the results of a reflection-transmission process generated by a layer sandwiched between homogeneous half-spaces. The layer is continuously stratified and allows for jump discontinuities at a finite number of planes. The dissipative effects are accounted for through the classical Boltzmann law of viscoelasticity. By using displacement and traction as convenient vector variables, the governing equations are considered in a vector Volterra integral equation and the solution is determined by means of a matricant. Next the matricant is applied to determine the reflection and transmission coefficients of a layer, with a generic piecewise continuous profile of the material properties. The reflection-transmission process produced by an obliquely incident wave, is considered for horizontally-polarized waves. The low-frequency approximation is derived for the reflection and transmission coefficients. Next, the high-frequency approximation is investigated by a WKB-like procedure which involves a complex valued frequency-dependent shear modulus. The displacement solution is obtained for the forward- and the backward-propagating waves in the layer along with the reflection and transmission coefficients.     

DUAL RECIPROCITY BOUNDARY ELEMENT METHOD FOR FLEXURAL WAVES IN THIN PLATE WITH CUTOUT

The theoretical analysis and numerical calculation of scattering of elastic waves and dynamic stress concentrations in the thin plate with the cutout was studied using dual reciprocity boundary element method (DRM). Based on the work equivalent law, the dual reciprocity boundary integral equations for flexural waves in the thin plate were established using static fundamental solution. As illustration, numerical results for the dynamic stress concentration factors in the thin plate with a circular hole are given. The results obtained demonstrate good agreement with other reported results and show high accuracy.     

A study of boundary integral method for a type of nonlinear problem based on generalized quasilinearlization theory

非线性泊松问题在热传导和多孔催化粒子的扩散反应等问题中是非常常见的,为此,利用广 义拟线性化迭代理论,提出了一种非线性泊松问题的新的数值迭代方法. 该方法将非线性方 程转化成一序列线性方程的迭代,其优点是初始值的选取具有一定的理论基础,并且在一定 的初始值条件下,迭代结果将单调地收敛于非线性问题的解. 将此迭代方法与边界元和双互 易杂交边界点方法结合,并用于非线性泊松问题的求解,比较了两种方法的结果精度,收敛 速度及不同初始值下的稳定性. 结果显示,基于拟线性化的双互易杂交边界点法具有较高的 稳定性和计算效率,并且收敛速度为平方阶.     

Natural convection in porous media—dual reciprocity boundary element method solution of the Darcy model

This paper describes the solution of a steady state natural convection problem in porous media by the dual reciprocity boundary element method (DRBEM). The boundary element method (BEM) for the coupled set of mass, momentum, and energy equations in two dimensions is structured by the fundamental solution of the Laplace equation. The dual reciprocity method is based on augmented scaled thin plate splines. Numerical examples include convergence studies with different mesh size, uniform and non‐uniform mesh arrangement, and constant and linear boundary field discretizations for differentially heated rectangular cavity problems at filtration with Rayleigh numbers of Ra*=25, 50, and 100 and aspect ratios of A=1/2, 1, and 2. The solution is assessed by comparison with reference results of the fine mesh finite volume method (FVM). Copyright © 2000 John Wiley & Sons, Ltd.     

Analysis of shear deformable plates with combined geometric and material nonlinearities by boundary element method

In this paper a boundary element formulation for analysis of shear deformable plates with combined geometric and material nonlinearities by boundary element method is presented. The dual reciprocity method is used in dealing with the geometric nonlinearity and domain discretization is implemented in dealing with material nonlinearity. The material is assumed to undergo large deflection with small strains. The von Mises criteria is used to evaluate the plastic zone and an elastic perfectly plastic material behaviour is assumed. An initial stress formulation is used to formulate the boundary integral equations. A total incremental method is applied to solve the nonlinear boundary integral equations. Numerical examples are presented to demonstrate the validity and the accuracy of the proposed method.     

SANDERS薄壳方程的自共轭性

本文对Sanders薄壳方程的自共轭性作了讨论.证明了以下三点: 1.通常的齐次边界条件是简单自共轭边界条件; 2.在简单自共轭边界条件下,Sanders薄壳方程是自共轭的、其蜕化(元矩)方程也是自共轭的; 3.任何薄壳理论,其满足功的互等定理与具有自共轭性所需条件是相同的: 作为6个变形分量的正定二次型的应变能函数存在. 由于Sanders薄壳理论在任意曲线坐标系中成立,故以上结论亦适用于任意曲线坐标系. 本文的讨论为采用Sanders理论对薄壳进行动力分析提供了理论准备.     

APPLICATION OF THE DUAL RECIPROCITY BOUNDARY ELEMENT METHOD TO SOLUTION OF NONLINEAR DIFFERENTIAL EQUATION

In this paper the dual reciprocity boundary element method is employed to solve nonlinear differential equation ∇² u+u+ɛu ³ =b. Results obtained in an example have a good agreement with those by FEM and show the applicability and simplicity of dual reciprocity method(DRM)in solving nonlinear differential equations.     

Multiple reciprocity method with two series of sequences of high-order fundamental solution for thin plate bending

IntroductionIt’swell_knownthatthecomplicatedfundamentalsolution[1,2 ]forHelmholtzequationΔu(x) +k2 u(x) =0　　(x∈Ω:boundedopenregioninR2 )isu (x,y) =-iH(2 )0 (k x-y ) 4,thusit’snotconvenientfornumericalcomputation .IfapplyingthesimplefundamentalsolutionofLaplaceequationu 0 (x ,y) =-ln｜x-y｜(2π) ,theexpressionforthesolutionofequationintheclosedregion Ωisc(y)u(y) + ∫Γu(x) u 0 (x,y) nx -u 0 (x ,y) u(x) n dsx =-k2∫Ωu(x)u 0 (x,y)dΩx.Astherightsideappearstheregionalintegrationinclu…     

Transition Matrix Method for Determining Stress Distribution of an Inhomogeneity in Elastostatic Medium

The transition matrix method has been extensively utilized to solve scattering in elastodynamic media. It is based on the reciprocity theorem, continuity of the interface boundary conditions, and applicable to arbitrary shape of inhomogeneity in systematic matrix multiplication. However, the transition matrix method has never been applied to determine stress distribution in elastostatic media. One important reason is the problem of the shortage of the basis functions of the elastostatic media that must be used to develop the transition matrix. This study investigates the required basis functions, and finds a set thereof that include Love??s special solutions of three dimensional elastostatics and three vector functions that are applicable to elastic waves. The proposed basis functions also can be adopted to derive the three significant orthogonality conditions for reciprocity at the surface of the inhomogeneity, which are useful in developing the transition matrix. The novel basis functions make the process of derivation of the T-matrix in elastostatics similar to that in elastodynamics. This process is illustrated for a spherical inhomogeneity that embedded in an elastic medium and stress patterns are compared with Goodier??s solutions, demonstrating high accuracy.     

Scattering and excitation of SH-surface waves by a protrusion at the mass-loaded boundary of an elastic half-space

A rigorous theory of the scattering and excitation of SH-surface waves by a protrusion at the mass-loaded boundary of an elastic half-space is presented. The boundary value problem (which is of the third kind) is solved by employing two suitably chosen Green functions. One of them is represented as a Fourier type of integral, the other is taken to be the Bessel function of the second kind and order zero. The procedure leads to a system of three, coupled, integral equations. This system is solved numerically. In case of an incident bulk wave, the amplitude of the launched surface wave is computed; in case of an incident surface wave, its transmission and reflection factor are computed. For both cases, an expression for the far-field radiation pattern of the scattered bulk wave is derived. A reciprocity relation is shown to exist between the amplitude of the launched surface wave and the far-field bulk wave radiation pattern. Numerical results are presented for a triangularly-prismatic protrusion; they are compared with the results pertaining to a corresponding indentation in the mass-loaded boundary, that have been obtained in a previous paper.     

Scattering of elastic waves by elastically transparent obstacles (integral-equation method)

A formulation of elastodynamic diffraction problems for sinusoidally in time varying disturbances in a linearly elastic medium is presented. Starting with the elastodynamic reciprocity relation, an integral representation for the particle displacement is derived. In it, the particle displacement and the traction at the boundary of the obstacle occur. From the integral representation, an associated integral equation is obtained by letting the point of observation approach the boundary of the obstacle. The “obstacle” may be either a rigid body, a void, or a body with elastic properties differing from those of its environment, or a combination of these. The integral equation thus obtained is well-suited for numerical treatment, when obstacles up to a few wavelengths in maximum diameter are considered.     

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.     

????????????????????????е????

提出了将杂交边界点法和双重互易法结合求解势问题的一种新的算法. 将势问题的解分为通解和特解两部分，通解使用 杂交边界点方法求解，特解则利用局部径向基函数近似. 该方法输入数据只是求解域上离散 的点，不需要额外的方程来计算域内物理量，后处理十分简便. 数值算例表明了该方法的稳 定性和有效性.     

Dual Reciprocity Boundary Element Method for Magnetohydrodynamic Flow Using Radial Basis Functions

A dual reciprocity boundary element method is given to obtain the solution in terms of velocity and induced magnetic field for the study of MHD (magnetohydrodynamic) flow through a rectangular duct having insulating walls. The equations are transformed to two types of nonlinear Poisson equations and the right-hand sides in these equations are approximated using combinations of two classes of radial basis functions (the value of the function and its normal derivatives are utilized for approximation). Computations are carried out for several values of the Hartman number (0 h M h 10) by using constant boundary elements. Comparisons are made for two types of formulations and for traditional and osculatory type approximations of the right-hand side functions. It is found that osculatory interpolation gives better results than traditional interpolation and the type of the Poisson equation, which contains derivative of the unknown function, is better than the other type, which contains unknown function only. The results for velocity and induced magnetic field are illustrated by some selected graphs.     

Écart à la réciprocité et identification de fissures en thermoélasticité isotrope transitoireReciprocity gap and identification of cracks in transitory isotropic thermoelasticity

We consider a three-dimensional homogeneous isotropic elastic solid containing a family of planar cracks submitted to a time-dependent thermal loading. The displacement and surface traction fields are measured over the whole external boundary of the solid. We propose in this Note to define and exploit a reciprocity gap, based only on the mechanical quantities available on the boundary, and which enables us to derive explicit formula for the location of the plane where the cracks are lying. Boundary conditions on the cracks can be of any nature provided they ensure that the normal heat flux and surface traction vector are continuous across the crack surfaces. To cite this article: S. Andrieux, H.D. Bui, C. R. Mecanique 334 (2006).