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.

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.

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.

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.

Axisymmetric Finite Element Method for Analysis of Plate on Layered Soil

To obtain the fundamental solution of soil has become the key problem for the semi-analytical and semi-numerical (SASN) method in analyzing plate on layered soil. By applying axisymmetric finite element method (FEM),an expression relating the surface settlement and the reaction of the layered soil can be obtained. Such a reaction can be treated as load acting on the applied external load. Having the plate modelled by four-node elements,the governing equation of the plate can be formed and solved. In this case, the fundamental solution can be introduced into the global soil stiffness matrix and five-node or nine-node element soil stiffness matrix.The existing commercial FEM software can be used to solve the fundamental solution of soil, which can bypass the complicated formula derivation and boasts high computational efficiency as well.

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.

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.

The 2D Fundamental Solutions in BEM Applied for Piezoelectric Materials

It is expected to solve the piezoelectric problems by using numerical method,e.g.,the boundary element method,since analytical solutions can hardly be obtained due to the deformation,electric field coupling effects of the material and anisotropy of the piezoelectric properties.However,the appropriate fundamental solutions used in boundary element method have not been found so far.In this paper,the solutions of displacement,potential,stress and dielectric displacement field effected by a unit mechanical force and a unit electric charge applied respectively at a point in a two-dimensional infinite domain are derived based on the coupling theory of electromagnetics and continuum mechanics.Hence,the eight indispensable fundamental solutions in BEM are provided.

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.

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.

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

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

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

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

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

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

双调和方程奇异边界法的改进及在Kirchhoff板弯曲中的应用

奇异边界法(SBM)是一种基于边界离散的无网格数值方法,在很多科学计算和工程领域中得到广泛的应用.该方法在处理复杂几何区域或者多连通区域时比基本解方法(MFS)数值计算更为稳定,具有易于实施、精度高等优点.SBM数值计算的关键之处在于源强度因子的计算,特别是相对于Laplace方程更为复杂的双调和方程的边界条件下源强度因子的计算.在高阶导数边界条件下,采用反插或者“加减项”原理计算源强度因子相对繁琐.本文对双调和方程的SBM进行了改进,将其中一个插值基函数改进为非奇异基函数形式,避免计算该基函数的源强度因子,极大简化了SBM的数值计算.本文改进对MFS同样有效,可以作为对传统MFS数值算法的补充.数值算例结果表明,本文提出的改进均能得到误差很小的数值解,且算法稳定,计算效率较高.

Application of Method of Fundamental Solutions in Solving Potential Flow Problems for Ship Motion Prediction

A novel panel-free approach based on the method of fundamental solutions (MFS) is proposed to solve the potential flow for predicting ship motion responses in the frequency domain according to strip theory. Compared with the conventional boundary element method (BEM), MFS is a desingularized, panel-free and integration-free approach. As a result, it is mathematically simple and easy for programming. The velocity potential is described by radial basis function (RBF) approximations and any degree of continuity of the velocity potential gradient can be obtained. Desingularization is achieved through collating singularities on a pseudo boundary outside the real fluid domain. Practical implementation and numerical characteristics of the MFS for solving the potential flow problem concerning ship hydrodynamics are elaborated through the computation of a 2D rectangular section. Then, the current method is further integrated with frequency domain strip theory to predict the heave and pitch responses of a containership and a very large crude carrier (VLCC) in regular head waves. The results of both ships agree well with the 3D frequency domain panel method and experimental data. Thus, the correctness and usefulness of the proposed approach are proved. We hope that this paper will serve as a motivation for other researchers to apply the MFS to various challenging problems in the field of ship hydrodynamics.

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.