边界积分方程中超奇异积分的解法

本文对边界积分方程中所存在的超奇异积分的数值解法作了综述，并介绍了它的一些应用。     

精细积分的非线性动力学积分方程及其解法

给出了非线性动力学积分方程的表达式,针对该方程提出了一个显式预测-校正的单步四阶精度的精细积分算法,适用于多自由度、强非线性,非保守系统.算例表明该方法精度高、计算量较少.     

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

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

非线性动力学积分方程分块积分解法

对于非线性动力学方程组分块地应用精细积分算法,使其化成积分方程表达式,求解的表达式中具有相对低阶的转换矩阵,从而使精细积分更适用于多自由度、强非线性、变系数、非保守系统,针对积分方程提出了一个显示预测-校正的单步四阶精度自起步的精细积分算法。算例表明本方法是有效的。     

三维有限体中平片裂纹的超奇异积分方程方法—Ⅱ.数值方法

超奇异积分方程方法的理论分析已在本文的第Ｉ部分中给出，这一部分是经的数值方法，及用此方法求解的若干典型的平片裂纹问题。     

一类奇异积分计算方法及其在断裂力学中的应用

一类奇异积分计算方法及其在断裂力学中的应用任传波,云大真（大连理工大学,大连１１６０２３）在力学及其它的工程计算中,常常遇到奇异积分,不同奇异程度的积分可以采用不同的方法来解决。本文提出的方法可以解决如下一类的奇异积分问题其中１求解方法对于式（１）的...     

等离子非线性扩散问题的积分方程数值分析

为了对等离子体密度非线性扩散方程进行数值分析,首先以差商代替微商,然后用对应方程的基本解,使之化为积分方程、再用边界单元法求解.     

三维断裂力学的超奇异积分方程方法

本文利用有限部积分的概念和方法,严格地证明了三维弹性体中受任意载荷作用的平片裂纹问题的超奇异积分方程组,并对未知解的性态作了理论分析,得到了性态指数,在此基础上通过主部分析,精确地求得了裂纹前沿光滑点附近的奇性应力场,从而找到了以裂纹面位移间断(位错)表示的应力强度因子表达式,最后对所得的超奇异积分方程组建立了数值法,并用此计算了若干典型的平片裂纹问题,数值结果令人满意。     

弹性力学平面问题中一类无奇异边界积分方程

从理论上提出一种新的方法,归化出间接变量无奇异边界积分方程. 采用Lagrange二次单元,建立一个数值求解框架系统. 此外,基于问题的计算区域的特殊性,给出一种边界近似方法. 数值算例表明该方法所取得的数值结果与精确解相当接近,特别是边界量的数值结果. 此外,该方法容易被推广到三维问题.和已有的直接变量的情形相比较,具有优点：1)无需处理HFP积分. 大大降低处理问题的复杂性,并提高了计算效率和解的精度;2)摆脱了问题的具体形式,进入纯代数操作.这样做的好处是从理论上建立一种普遍适用的方法,不仅适用于弹性力学问题,同样可应用于其它问题,如位势问题, Stokes问题等. 3)提供了一种计算CPV积分的方法.     

一类非线性奇异振动方程的渐近解

本文讨论了一类型典的非一性奇异振动方程－ＷＣＭ方程的斩近解法、分析了该方程解的稳定性问题，并将结果与摄动法求解的结果以及数值积分的结果相比较以说明所用方法能用于此类非线性奇异方程的分析。     

A numerical method for fractional integral with applications

IntroductionThefractionalcalculushasalonghistoryandthereareamassofworkstodiscussthefractionalderivativesandfractionalintegralswitharbitrary (realorcomplex)order[1- 3 ].Thefractionalcalculushasawideapplicationbackground ,especiallyinthefieldsofchemistry ,electromagnetics,materialscienceandmechanics.Forexample,Gement[4 ]proposedthefractionalderivativeconstitutivemodelsofaviscoelasticmaterialatfirst.Themodelshavereceivedincreasingattention[5 - 7].Onlyafewparametersarecontainedinthemodelsandthemo…     

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method (LBIEM),singularities in the local boundary integrals need to be treated specially. In the current paper,local integral equations are adopted for the nodes inside the domain trod moving least square approximation (MLSA) for the nodes on the global boundary,thus singularities will not occur in the new al- gorithm.At the same time,approximation errors of boundary integrals are reduced significantly.As applications and numerical tests,Laplace equation and Helmholtz equa- tion problems are considered and excellent numerical results are obtained.Furthermore, when solving the Hehnholtz problems,the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions.Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.     

Perturbation finite volume method for convective-diffusion integral equation

A perturbation finite volume (PFV) method for the convective-diffusion integral equation is developed in this paper. The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order integration approximations, with the least nodes similar to the standard three-point schemes, that is, the number of the nodes needed is equal to unity plus the face-number of the control volume. For instance, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D linear and nonlinear problems, 2-D and 3-D flow model equations. Comparing with other standard three-point schemes, the PFV scheme has much smaller numerical diffusion than the first-order upwind scheme (UDS). Its numerical accuracies are also higher than the second-order central scheme (CDS), the power-law scheme (PLS) and QUICK scheme. The project supported by the National Natural Science Foundation of China (10272106, 10372106)     

New numerical method for volterra integral equation of the second kind in piezoelastic dynamic problems

The elastodynamic problems of piezoelectric hollow cylinders and spheres under radial deformation can be transformed into a second kind Volterra integral equation about a function with respect to time, which greatly simplifies the solving procedure for such elastodynamic problems. Meanwhile, it becomes very important to find a way to solve the second kind Volterra integral equation effectively and quickly. By using an interpolation function to approximate the unknown function, two new recursive formulae were derived, based on which numerical solution can be obtained step by step. The present method can provide accurate numerical results efficiently. It is also very stable for long time calculating.     

A numerical method for Cauchy principal value of singular integral under parametric transformation in BEM

A numerical method for the Cauchy principal value of the singular integral in BEM is developed using the concept of finite part integration under integral variable transformation. It is applied to the numerical integration on isoparametric element successfully, as shown in the examples in this paper.The project supported by National Natural Foundation of China.     

A perturbation boundary integral equation method for analysis of interaction between foundation and nonlinear rock and soil medium

A scheme is developed for analysing the interaction between a foundation and a nonlinear rock and soil medium, in which the foundation is considered as a linear elastic body and a typical boundary integral equation method (BIEM) is employed. On the basis of taking the nonlinear properties of the medium into account, a perturbation BIEM is developed. The fundamental equations for the nonlinear coupling analysis are formulated, and typical problems are solved and discussed by the present method.     

The first integral method for solving some important nonlinear partial differential equations

Exact solutions of some important nonlinear partial differential equations are obtained by using the first integral method. The efficiency of the method is demonstrated by applying it for two selected equations.     

Impact of singularity of Navier-Stokes equation upon atmospheric motion equations

Some conclusions about the smooth function classes stability for the basic system of equations of atmospheric motion and instability for Navier-Stokes equation are summarized. On the basis of this, by taking the basic system of equations of atmospheric motion via Boussinesq approximation as example to explain in detail that the instability about some simplified models of the basic system of equations for atmospheric motion is caused by the instability of Navier-Stokes equation, thereby, a principle to guarantee the stability of simplified equation is drawn in simplifying the basic system of equations.     

Isotropicalized spline integral equation method for the analysis of anisotropic plates

In the view of Reissner’s and Kirchhoff’s theories,respectively,we formulate theisotropicalized governing equations for the anisotropic plates,and give the proof of theequivalence relation between these two plate-models for the simply-supported rectangularorthotropic plates.The well-known fundamental solutions of the isotropic plates are utlizedfor the spline integral equation analysis of anisotropic plates.Even with sparse meshes thesatisfactory results can be obtained.The analysis of plates on two-parameter elasticfoundation is so simple as the common case that only a few terms should be added to theformulas of fictitious loads.     

Two-dimensional polynomial eigenstrain formulation of boundary integral equation with numerical verification

The low-order polynomial-distributed eigenstrain formulation of the boundary integral equation (BIE) and the corresponding definition of the Eshelby tensors are proposed for the elliptical inhomogeneities in two-dimensional elastic media. Taking the results of the traditional subdomain boundary element method (BEM) as the control, the effectiveness of the present algorithm is verified for the elastic media with a single elliptical inhomogeneity. With the present computational model and algorithm, significant improvements are achieved in terms of the efficiency as compared with the traditional BEM and the accuracy as compared with the constant eigenstrain formulation of the BIE.