Some Remarks on the Method of Fundamental Solutions for Two- Dimensional Elasticity

In this paper,some remarks for more efficient analysis of two-dimensional elastostatic problems using the method of fundamental solutions are made.First,the effects of the distance between pseudo and main boundaries on the solution are investigated and by a numerical study a lower bound for the distance of each source point to the main boundary is suggested.In some cases,the resulting system of equations becomes ill-conditioned for which,the truncated singular value decomposition with a criterion based on the accuracy of the imposition of boundary conditions is used.Moreover,a procedure for normalizing the shear modulus is presented that significantly reduces the condition number of the system of equations.By solving two example problems with stress concentration,the effectiveness of the proposed methods is demonstrated.

SOME PROBLEMS WITH THE METHOD OF FUNDAMENTAL SOLUTION USING RADIAL BASIS FUNCTIONS

The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson＇s equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition （SVD） technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric （MQ） are fully analyzed to provide useful reference.

TWO-DIMENSIONAL ELECTROELASTIC FUNDAMENTAL SOLUTIONS FOR GENERAL ANISOTROPIC PIEZOELECTRIC MEDIA

Explicit fomulas for 2-D electroelastic fundamental solutions in general anisotropic piezoelectric media subjected to a line force and a line charge are obtained by using the plane wave decomposition method and a subsequent application of the residue calculus. 'Anisotropic' means that any material symmetry restrictions are not assumed. 'Two dimensional' includes not only in-plane problems but also anti-plane problems and problems in which in-plane and anti-plane deformations couple each other. As a special case, the solutions for transversely isotropic piezoelectric media are given.

Fundamental solution method for inverse source problem of plate equation

The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, which is referred to be an inverse source problem of a plate equation. The uniqueness theorem for this problem is stated, and the fundamental solution to the plate equation is derived. In the case that the plate is driven by the harmonic load, the fundamental solution method （FSM） and the Tikhonov regularization technique axe used to calculate the source term. Numerical experiments of the Euler-Bernoulli beam and the Kirchhoff-Love plate show that the FSM can work well for practical use, no matter the source term is smooth or piecewise.

The 2D Fundamental Solutions in BEM Applied for Piezoelectric Materials

Ｔｈｅ２ＤＦｕｎｄａｍｅｎｔａｌＳｏｌｕｔｉｏｎｓｉｎＢＥＭＡｐｐｌｉｅｄｆｏｒＰｉｅｚｏｅｌｅｃｔｒｉｃＭａｔｅｒｉａｌｓ￥（孟庆元）（杜善义）ＭＥＮＧＱｉｎｇｙｕａｎ；ＤＵＳｈａｎｙｉ（Ｄｅｐｔ．ｏｆＡｓｔｒｏｎａｕｔｉｃｓａｎｄＭｅｃｈａｎｉｃ...

MULTIPLE RECIPROCITY METHOD WITH TWO SERIES OF SEQUENCES OF HIGH-ORDER FUNDAMENTAL SOLUTION FOR THIN PLATE BENDING

The boundary value problem of plate bending problem on two_parameter foundation was discussed.Using two series of the high_order fundamental solution sequences, namely, the fundamental solution sequences for the multi_harmonic operator and Laplace operator, applying the multiple reciprocity method(MRM), the MRM boundary integral equation for plate bending problem was constructed. It proves that the boundary integral equation derived from MRM is essentially identical to the conventional boundary integral equation. Hence the convergence analysis of MRM for plate bending problem can be obtained by the error estimation for the conventional boundary integral equation. In addition, this method can extend to the case of more series of the high_order fundamental solution sequences.

The Localized Method of Fundamental Solution for Two Dimensional Signorini Problems

In this work,the localized method of fundamental solution(LMFS)is extended to Signorini problem.Unlike the traditional fundamental solution(MFS),the LMFS approximates the field quantity at each node by using the field quantities at the adjacent nodes.The idea of the LMFS is similar to the localized domain type method.The fictitious boundary nodes are proposed to impose the boundary condition and governing equations at each node to formulate a sparse matrix.The inequality boundary condition of Signorini problem is solved indirectly by introducing nonlinear complementarity problem function(NCP-function).Numerical examples are carried out to validate the reliability and effectiveness of the LMFS in solving Signorini problems.

Fundamental Solution of Dirichlet Boundary Value Problem of Axisymmetric Helmholtz Equation

Fundamental solution of Dirichlet boundary value problem of axisymmetric Helmholtz equation is constructed via modi?ed Bessel function of the second kind, which uni?ed the formulas of fundamental solution of Helmholtz equation, elliptic type Euler-Poisson-Darboux equation and Laplace equation in any dimensional space.

THE THREE-DIMENSIONAL FUNDAMENTAL SOLUTION TO STOKES FLOW IN THE OBLATE SPHEROIDAL COORDINATES WITH APPLICATIONS TO MULTIPLES SPHEROID PROBLEMS

A new three-dimensional fundamental solution to the Stokes flow was proposed by transforming the solid harmonic functions in Lamb's solution into expressions in terms Of the oblate spheroidal coordinates. These fundamental solutions are advantageous in treating flows past an arbitrary number of arbitrarily positioned and oriented oblate spheroids. The least squares technique was adopted herein so that the convergence difficulties often encountered in solving three-dimensional problems were completely avoided. The examples demonstrate that present approach is highly accurate, consistently stable and computationally efficient. The oblate spheroid may be used to model a variety of particle shapes between a circular disk and a sphere. For the first time, the effect of various geometric factors on the forces and torques exerted on two oblate spheroids were systematically studied by using the proposed fundamental solutions. The generality of this approach was illustrated by two problems of three spheroids.

Fundamental Solution for Welding Problem by Two Dissimilar Isotropic Semi－Planes

A fundamental solution was obtained for an infinite plane bonded by two dissimilar isotropic semi-planes by employing plane elastic complex variable method and theory of boundary value problems for analytic functions.Fundamental solution was prepared for solving these types of problems with boundary element method.

A New Approximate Fundamental Solution for Orthotropic Plate

A weight double trigonometric series is presented as an approximate fundamental solution for orthotropic plate.Integral equation of orthotropic plate bending is solved by a new method, which only needs one basic boundary integral Eq., puts one fictitious boundary outside plate domain. Examples show that the approximate fundamental solution and solving method proposed in this paper are simple, reliable and quite precise. And they are applicable for various boundary conditions.

Electroelastic Analysis of Two-Dimensional Piezoelectric Structures by the Localized Method of Fundamental Solutions

Accurate and efficient analysis of the coupled electroelastic behavior of piezoelectric structures is a challenging task in the community of computational mechanics.During the past few decades,the method of fundamental solutions(MFS)has emerged as a popular and well-established meshless boundary collocation method for the numerical solution of many engineering applications.The classical MFS formulation,however,leads to dense and non-symmetric coefficient matrices which will be computationally expensive for large-scale engineering simulations.In this paper,a localized version of the MFS(LMFS)is devised for electroelastic analysis of twodimensional(2D)piezoelectric structures.In the LMFS,the entire computational domain is divided into a set of overlapping small sub-domains where the MFS-based approximation and the moving least square(MLS)technique are employed.Different to the classical MFS,the LMFS will produce banded and sparse coefficient matrices which makes the method very attractive for large-scale simulations.Preliminary numerical experiments illustrate that the present LMFM is very promising for coupled electroelastic analysis of piezoelectric materials.

Localized Method of Fundamental Solutions for Acoustic Analysis Inside a Car Cavity with Sound-Absorbing Material

This paper documents the first attempt to apply a localized method of fundamental solutions(LMFS)to the acoustic analysis of car cavity containing soundabsorbing materials.The LMFS is a recently developed meshless approach with the merits of being mathematically simple,numerically accurate,and requiring less computer time and storage.Compared with the traditional method of fundamental solutions(MFS)with a full interpolation matrix,the LMFS can obtain a sparse banded linear algebraic system,and can circumvent the perplexing issue of fictitious boundary encountered in the MFS for complex solution domains.In the LMFS,only circular or spherical fictitious boundary is involved.Based on these advantages,the method can be regarded as a competitive alternative to the standard method,especially for high-dimensional and large-scale problems.Three benchmark numerical examples are provided to verify the effectiveness and performance of the present method for the solution of car cavity acoustic problems with impedance conditions.

流体饱和半空间中不均质体对地震波散射的MFS求解

针对饱和半空间中任意形状不均质体对地震波的散射问题,采用一种新型基本解方法(MFS)进行求解分析。该方法结合饱和半空间中膨胀波源和剪切波源的格林函数,首先由分布在不均质体和半空间交界面附近两虚拟波源面上的PI、PII、SV波源分别构造了不均质体内外的散射波场,然后由交界面连续性条件建立方程并求解确定了虚拟波源密度,总波场反应由自由波场和散射波场叠加而得,最后在精度检验的基础上,通过一组典型算例研究了平面P波在饱和半空间中不均质体周围散射的基本规律。研究结果表明:两相介质波动的MFS模拟具有极高的数值精度和高频稳定性;波在饱和不均质体中的散射特征取决于入射波的频率和方向、边界渗透条件、介质孔隙率等参数,与在单相介质中的情况具有很大差别;地表位移特征反映了波的散射和饱和不均质体的自振特性,随着介质孔隙率增大,位移幅值谱振荡更为剧烈,且不均质体刚度越小,地表位移放大效应越显著,但最大不超过自由场反应的2倍。

基于快速多极子基本解方法(FMM-MFS)的弹性波二维散射模拟研究

针对弹性波二维散射问题,发展一种新的快速多极子基本解方法(FMM-MFS)。方法基于单层位势理论,通过在虚边界上设置膨胀波线源和剪切波线源以构造散射波场,从而避免了奇异性的处理和边界单元离散;结合快速多极子展开技术(FMM),大幅度降低了计算量和存储量,突破了传统方法难以处理大规模散射问题的瓶颈。以全空间孔洞对P、SV波的二维散射为例,给出了具体求解步骤,并在个人计算机上实现了上百万自由度问题的快速精确计算。在方法效率和精度检验基础上,分别以单孔洞和随机孔洞群对平面波(P、SV波)的散射为例进行计算模拟,揭示了孔洞(群)周围弹性波散射的若干重要规律。

薄板大挠度弯曲问题的DRM-MFS无网格方法

在介绍薄板大挠度问题的控制方程及其推导过程的基础上,根据渐近迭代方法,分别用DRM方法和MFS方法近似特解和齐次解,得到求解薄板大挠度弯曲问题的DRM-MFS无网格近似方法,再通过数值算例验证方法的有效性及准确性.

Integral Representations for the Solutions of the Generalized Schroedinger Equation in a Finite Interval

We reduce the initial value problem for the generalized Schroedinger equation with piecewise-constant leading coefficient to the system of Volterra type integral equations and construct new useful integral representations for the fundamental solutions of the Schroedinger equation. We also investigate some significant properties of the kernels of these integral representations. The integral representations of fundamental solutions enable to obtain the basic integral equations, which are a powerful tool for solving inverse spectral problems.

基于蒙特卡洛法的Euler-Bernoulli梁基频和振型求解方法

将Rayleigh法和蒙特卡洛法相结合,在Euler-Bernoulli梁理论假设下求解了均匀梁、变截面梁和附带集中质量的变截面梁自由振动问题.对原本连续的梁结构模型进行离散化处理,利用蒙特卡洛法给出梁结构的假设振型.将假设得到的梁结构振型函数代入Rayleigh法,多次计算过程中,将历次基频所得值与计算所得最小值进行比较,根据其相对误差判断是否满足收敛条件,进而求得基频及对应的振型.结果表明,不同计算模型中基频最大误差不超过10%,能够满足工程需求,且精度和时间的控制参数调整灵活,使用者可根据自身需要自行调节.该方法理论简明,适用范围广泛,能够快速准确地求解诸多类型的梁结构基频和振型.

Localized Method of Fundamental Solutions for Three-Dimensional Elasticity Problems: Theory

A localized version of the method of fundamental solution(LMFS)is devised in this paper for the numerical solutions of three-dimensional(3D)elasticity problems.The present method combines the advantages of high computational efficiency of localized discretization schemes and the pseudo-spectral convergence rate of the classical MFS formulation.Such a combination will be an important improvement to the classical MFS for complicated and large-scale engineering simulations.Numerical examples with up to 100,000 unknowns can be solved without any difficulty on a personal computer using the developed methodologies.The advantages,disadvantages and potential applications of the proposed method,as compared with the classical MFS and boundary element method(BEM),are discussed.