导数场边界积分方程分析近边界应力分布

研究二维弹性力学问题边界积分方程,通过分部积分变换消除了常规导数边界积分方程中的超奇异积分,获得仅含强奇异积分的应力自然边界积分方程.对于近边界应力的计算,进一步运用正则化算法解析计算其中的几乎强奇异积分.较常规边界元法相比,应力自然边界积分方程可以求解离边界更加接近的内点应力值.算例证明了文中方法的可应用性和有效性.     

弹性理论几类导数边界积分方程之间的变换关系

导数场边界积分方程通常难以应用，因为存在着超奇异主值积分的计算障碍。弹性理论中有几类不同的位移导数边界积分方程，本文采用算子δij和∈ij(排列张量)作用于这些导数边界积分方程，做一系列变换，原有的超奇异积分被正则化为强奇异积分获解。从而建立了这些位移导数边界积分方程之间的转换关系，它们均可以归结为自然边界积分方程。自然边界积分方程仅存在容易计算的Cauchy主值积分。自然边界积分方程分析可直接获得边界应力和位移导数。     

含裂纹粘接结构的边界元应力分析

本文采用子域法和直接应力奇异单元法求解二维粘接结构中的裂纹问题。子域法把粘接结构划分为三个子域,根据每个子域的边界积分方程和子域间的界面条件,可以建立粘接结构的边界积分方程组。直接应力奇异单元能够在整个单元长度上反映裂纹端部的1/r~(1/2)奇异性,在计算时可以通过坐标变换消除奇异单元积分中的奇异性,直接计算出应力强度因子。含裂纹多层结构的数值示例和粘接补强单边裂纹板的应力测试和疲劳试验结果证实了本文方法的有效性。     

关于时域边界元法的若干研究

本文从三维弹性动力学方程的基本奇异解着手，导出了适于计算机计算的求解三弹性动力学问题的边界积分方程（ＢＩＥ），并在理论上提出了ＩＢＥ前缘系数矩阵（５）具有（１）准对角特性，（２）其各元素不随时间而变化。据此，本文给出了用时域边界元法求解弹性动力学问题的新方法。最后，数值算例验证了本文方法的正确性。     

边界元法中计算几乎奇异积分的一种无奇异算法

边界元法中存在几乎奇异积分的计算困难。引起边界单元上几乎奇异积分的因素是源点到其邻近单元的最小距离δ。本文拓展文[1]的思想,进一步采用分部积分将δ移出奇异积分式中积分核之外,转换后积分核是δ的正则函数。所以几乎强奇异和超奇异积分被化为无奇异的规则积分与解析积分的和,可由通常的Gauss数值积分解。文中应用此正则化技术求解了弹性力学平面问题的近边界点位移和应力。     

关于边界元法中奇异积分的处理

关于边界元法中奇异积分的处理臧跃龙，嵇醒（西安交通大学，７１００４９）（上海同济大学，上海２０００９２）关键词边界元法，积分奇异性，奇异性的消除１引言与有限元法不同，边界元法的数值积分通常带有对数或一、二阶奇异性．如二维问题，基本解带有对数奇异性，其...     

边界元法分析狭长体结构

针对边界元法分析狭长结构时遇到的几乎奇异积分难以计算的困难，将几乎奇异积分划分为两种类型，分别通过分部积分交换把引起积分几乎奇异的参量移至积分号之外，从而建立了一个新的正则化算法，解决了边界积分方程中几乎奇异积分的计算难题。文中用边界元法计算了弹性力学平面问题的狭长结构，算例证明了本法的有效性。     

薄体位势问题边界元法中的解析积分算法

薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法，用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元，边界积分方程离散计算的积分可归纳为三种形式。对薄体问题，源点与积分单元距离通常相距很近，这些积分产生显著几乎奇异性，直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零，公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分，首先精确计算出薄体问题边界未知位势和法向位势梯度，然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1．E-08的薄体问题，显示了边界元法分析薄体问题具有独特的优势。     

二维边界元奇异积分和多域缩聚法分析

基于基本解的一种新的表达式，对二维边界元分析中奇异积分的精确求解进行了讨论，从几何方面对基本解的奇异性进行了分析，给出了超参非连续元离散位势和弹性力学问题边界积分方程时奇异积分计算的精确式，从而为判断各种近似方法的优劣和间接方法的精度提供了依据，也为精确地分析了大规模问题提供了一条有效的途径。     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

Boundary element method to solve torsion problem of cracked cylinder

Separating the discontinuous solution by use of the single crack solution, together with the regular solution of harmonic function, the torsion problem of a cracked cylinder is reduced to solving a set of mixed-type integral equations and its numerical technique is then proposed by combining the numerical method of singular integral equation with the boundary element method. Several numerical examples are calculated which will be useful to engineering practice. The method proposed is characterized by its fine accuracy and convenience for using, which can be extended to the cases of multiple crack.The project supported by National Natural Science Foundation of China.     

Stress intensity factors of the eccentric and edge cracked cylinders

Using the single crack solution and the regular solution of harmonic function,thetorsion problem of a cracked cylinder is reduced to solving a set of mixed-type integralequations which can be solved by combining the numerical method of singular integralequation with the boundary element method.Several numerical examples arecalculated and the stress intensity factors are obtained.     

基于时间依赖基本解的奇异边界法模拟二维狄利克雷边界标量波方程

奇异边界法是一个半解析边界配点强格式方法,具有无数值积分和无网格、编程容易以及数学简单等优点。本文首次将时间依赖基本解运用于奇异边界法,计算模拟二维标量波方程;结合确定源点强度因子的反插值技术,提出了二维狄利克雷边界标量波方程源点强度因子的一个经验公式;引进了解决波方程基本解G奇异性的一种无奇异积分处理方法。数值实验证明,基于时间依赖基本解的奇异边界法可精确高效地模拟二维狄利克雷边界标量波方程,在计算效率、精度、稳定性和适应性等方面有明显优势。     

Singular integral equations and boundary element method of cracks in thermally stressed planar solids

Ｉｔｉｓｋｎｏｗｎｔｈａｔｍｏｓｔａｇｒｉｃｕｌｔｕｒａｌｐｒｏｄｕｃｔｓａｎｄｆｏｏｄｓａｒｅｐｒｏｃｅｓｓｅｄａｎｄｔｒａｎｓｐｏｒｔｅｄｕｎｄｅｒｃｅｒｔａｉｎｔｅｍｐｅｒａｔｕｒｅｃｏｎｄｉｔｉｏｎｓ,ａｎｄｔｈｅｓｔｒｕｃｔｕｒａｌｃｏｍｐｏｎｅｎｔｓａｌｓｏｗｏｒｋｕｎｄｅｒａｔｈｅｒｍａｌｅｎｖｉｒｏｎｍｅｎｔ.Ｔｅｍｐｅｒａｔｕｒｅｉｎｄｕｃｅｄｓｔｒｅｓｓｅｓｕｓｕａｌｌｙｌｅａｄｔｏｄａｍａｇｅｏｆｆｌａｗｅｄｓｏｌｉｄｓ.Ｔｈｕｓ,ｔｈｅｉｎｖｅｓｔｉｇａｔｉｏｎｏｆｔｈｅｃｒ…     

压电材料平面问题的虚边界元-等额配点解法

利用压电材料平面问题的基本解和弹性力学虚边界元方法的基本思想,提出了压电材料平面问题的虚边界元-等额配点解法。该解法继承了传统边界元方法的优点,而避免了传统边界元方法遇到的边界积分奇异性问题。最后给出了压电材料平面问题的一些具体算例,并与解析解作了比较。结果表明本文的方法有很高的精度,是该问题一个十分有效的数值求解方法。     

双材料基本解弹塑性边界元法

提出并研究采用双材料基本解的弹塑性边界元法,得到了内点应力公式中有关奇点塑性应变自由项的完整表达式,并利用非连续边界单元和非连续区域单元解决了当奇点位于界面上时该自由项难于确定,以及计算区域Cauchy主值积分的常塑性应变场法在与界面相连的奇异区域单元上无法实施的困难．采用双材料基本解的弹塑性边界无法针对双材料的结构特点,特别适于分析有关弹塑性双材料界面及界面裂纹问题．     

二维位势边界元法高阶单元几乎奇异积分半解析算法

准确计算几乎奇异积分是边界元法难题之一。目前,对于一般的高阶单元的几乎奇异积分尚缺乏通用高效的计算方法。本文在单元局部坐标系中表征了二维高阶单元的几何特征,提出了源点相对高阶单元的接近度概念。针对二维位势边界元法的3节点二次等参单元,构造出与单元积分核具有相同几乎奇异性的近似奇异核函数。从二维位势几乎奇异积分单元积分核中扣除近似奇异核函数,把几乎奇异积分项转换为规则积分和奇异积分两部分之和,规则积分部分用常规Gauss数值积分计算,奇异积分部分由导出的解析公式计算,从而建立了二维位势问题高阶单元几乎强奇异和超奇异积分的半解析算法。算例结果表明了本文半解析算法的有效性和计算精度。     

层体边界元法在二维层体断裂分析中的应用

在刚度矩阵法的基础上建立了用于进行二维多层体结构断裂分析的边界单元法（ＢＥＭＬＭ）由于ＢＥＭＬＭ的基本方程中已经包含了层体表面和裂纹缝面的边界条件，因而不需要对这些边界进行单元离散，从而其断裂分析可望有较好的精度通过与柯西积分方程法进行结合，算例表明ＢＥ ＭＬＭ是可靠并有效的     

一种改进格林元方法及在渗流问题中的应用

采用格林公式和基本解推导出直接边界积分方程来求解渗流问题.边界积分方程数值离散基于格林元方法(Green element methond),改进了原方法中压力和压力导数的求解方法,命名为混合边界元方法(Mixed boundary element method).相较于格林元类方法,该方法显式考虑了求解节点的外法向流量值和压力值,并使求得的数值解在求解区域上能够连续,符合实际的物理过程,在不增加额外未知数的情况下提高了计算精度.分析了不同网格类型对模拟计算结果的影响,并对稳定渗流问题、非稳定(瞬态)渗流问题和非稳态问题进行了实例计算,结果显示改进方法提高了计算精度,并对各类渗流问题有较好的适应性.     

THE MESHLESS VIRTUAL BOUNDARY METHOD AND ITS APPLICATIONS TO 2D ELASTICITY PROBLEMS

A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network (IRBFN) constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations (BIE) can be established while the fundamental solutions to the problems are introduced. The singular value decomposition (SVD) method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.