不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

加权残值法对Stokes流的应用

本文通过综述加权残值法对流体力学一个分支——Stokes流动的应用,试图说明这一方法对于解决流体力学问题具有广阔的前景。      

无奇异性边界积分方程法

边界元法中如何有效地处理奇异积分,一直是人们极为关心的课题。本文提出了建立由“线源”产生的边界积分方程,其优点是可任意阶地降低奇异性阶次。计算实例表明,本文所提出的方法优于常规的边界元法。尤其是提高了位移的计算精度。对边界单元网格不等长划分,更显示出该法的优越性。     

积分方程的加权残数配点法

对于用积分方程式表示的物理问题,本文应用加权残数的配点法进行求解,使得求解方法简便、可靠、正确。与常规的求解微分方程的加权残数配点法相比较,本方法不需假定任何试函数。文中给出较多的算例,论证了本文提出的积分方程的加权残数配点法的有效性,可靠性和精确性。     

积分单元法计算双柱绕流

本文利用解不可压纳维尔-斯托克斯方程的积分单元法计算了双柱绕流的流场.结果表明,积分单元法是计算柱群绕流的一种方法.     

各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程，数值计算结果表明在退化尺度附近，充要的随机边界积分程较习用的随机边界积分方程有较大的优越性。     

各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程。数值计算结果表明在退化尺度附近，充要的随机边界积分方程较习用的随机边界积分方程有较大的优越性。     

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

边界积分方程的若干基本理论问题

本文提出了边界积分方程理论体系框图,对其中研究得不多的和容易引起误解的环节,如基本解的定义、解的边界积分表示、超定问题有解的充要条件等,作了专门的讨论.     

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

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

样条积分方程法分析弹塑性板弯曲

基于形变理论和Mises准则本文用虚载法分别导出Reissner型和Kirchhoff型板弹塑性弯曲方程,对它们间在多边形简支和轴对称弯曲下的相通性给出论证,并用弹性板样条积分方程法来求解,对诸如塑性域的范围和深度以及各点的弹塑性内力和位移等即使在稀疏剖分下也能有良好的计算精度。     

广义Stokes方程的完全边界积分表示式及其在求解N-S方程中的应用

为了分解N-S方程组各变量相互偶合,本文采用Peaceman-Rachford算子分裂法,将时间相依的N-S方程组分解成不存在上述偶合特性的线性和非线性的子问题。线性子问题具有广义Stokes方程类型。本文采用多重互易法,即采用多阶拉普拉斯算子基本解逐步变换,将其解表示成完全边界积分形式,从而使问题的计算维数降低一维。广义Stokes方程的算例以及二维圆柱在剪切流中的Stokes绕流解,都表明多重互易算法具有高效特点,而且后者与文[3]解析解吻合得非常好。     

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

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

AN EXTENDED BOUNDARY INTEGRAL EQUATION METHOD FOR THE REMOVAL OF IRREGULAR FREQUENCY EFFECTS

Numerical techniques for the analysis of wave–body interactions are developed by the combined use of two boundary integral equation formulations. The velocity potential, which is expressed in a perturbation expansion, is obtained directly from the application of Green's theorem (the ‘potential formulation’), while the fluid velocity is obtained from the gradient of the alternative form where the potential is represented by a source distribution (the ‘source formulation’). In both formulations, the integral equations are modified to remove the effect of the irregular frequencies. It is well known from earlier works that if the normal velocity is prescribed on the interior free surface, inside the body, an extended boundary integral equation can be derived which is free of the irregular frequency effects. It is shown here that the value of the normal velocity on the interior free surface must be continuous with that outside the body, to avoid a logarithmic singularity in the source strength at the waterline. Thus the analysis must be carried out sequentially in order to evaluate the fluid velocity correctly: first for the velocity potential and then for the source strength. Computations are made to demonstrate the effectiveness of the extended boundary integral euations in the potential and source formulations. Results are shown which include the added-mass and damping coefficients and the first-order wave-exciting forces for simple three-dimensional bodies and the second-order forces on a tension-leg-platform. The latter example illustrates the importance of removing irregular frequency effects in the context of second-order wave loads.     

边界积分方程法在动态断裂力学中的数值计算研究

本文将D.Nardini和C.A.Brebbia所提出的一种动态边界积分方程新解法应用于动态断裂力学数值计算，对数值实现问题，尤其对数值稳定性及精度问题进行了详细研究，得到了能保证数值稳定性的数值解法，给出了动态断裂力学计算实例，同己有的数值结果比较，表明本文的计算是成功的。     

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years .It uses the moving least square (MLS) approximation as a shape function . The smoothness of the MLS approximation is determined by that of the basic function and of the weight function, and is mainly determined by that of the weight function. Therefore, the weight function greatly affects the accuracy of results obtained. Different kinds of weight functions, such as the spline function, the Gauss function and so on, are proposed recently by many researchers. In the present work, the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method. The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed. Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and a in Gauss and exponential weight functions are in the range of reasonable values, respectively, and the higher the smoothness of the weight function, the better the features of the solutions.     

时程积分过程中结构运动参数的协调

通过在时程积分过程中补充一个能量方程,并借助原时程积分方法的过程数据,提出一种协调时间步长内结构计算位移、速度和加速度的方法.该方法蕴涵了对原时程积分方法计算假设条件的修正,使之更符合实际.以wilson-θ法为例,加入协调过程后,提高了计算精度和减弱计算结果的超越现象.算例结果表明时间步长内结构的计算运动参数不协调,是时程积分结果误差的一个主要原因.