非齐次弹性力学问题双互易边界元方法研究北大核心CSCD

基于弹性力学边界元方法理论,将边界元法与双互易法结合,采用指数型基函数对非齐次项进行插值得到双互易边界积分方程.将边界积分方程离散为代数方程组,利用已知边界条件和方程特解求解方程组,得出域内位移和边界面力.指数型基函数的形状参数是由插值点最近距离的最小值决定,采用这种形状参数变化方案,分析径向基函数(RBF)插值精度以及插值稳定性.再次将指数型基函数应用到双互易边界元法中,分析双互易边界元方法下计算精度及稳定性,验证了指数型插值函数作为双互易边界元方法的径向基函数解决弹性力学域内体力项问题的有效性.     

一类非线性传输问题的区域分解方法（英文）

针对一类非线性传输问题提出了有限元与边界元的耦合方法并设计了基于耦合法的区域分解算法.该算法避免了求解边界积分方程,从而计算量大大减少.算法的收敛性分析和数值算例验证了该算法的合理和有效性.     

焊接过程瞬态温度场与残余应力场

本文根据焊接热传导及热弹塑性的研究现状,对如下两个问题进行了探讨.首先,采用非线性问题线性化的方法对非线性非定常温度场问题的边界元法做了改进,并将其用于焊接热传导分析:其次,提出“等放线膨胀系数法”考虑相变对于应力场的影响,并将其用于焊接热弹塑性分析.实例的数值计算结果与实测数据的对比分析表明本文方法行之有效.     

求解双材料裂纹结构全域应力场的扩展边界元法

在线弹性理论中,复合材料裂纹尖端具有多重应力奇异性,常规数值方法不易求解.该文建立的扩展边界元法(XBEM)对围绕尖端区域位移函数采用自尖端径向距离r的渐近级数展开式表达,其幅值系数作为基本未知量,而尖端外部区域采用常规边界元法离散方程.两方程联立求解可获得裂纹结构完整的位移和应力场.对两相材料裂纹结构尖端的两个材料域分别采用合理的应力特征对,然后对其进行计算,通过计算结果的对比分析,表明了扩展边界元法求解两相材料裂纹结构全域应力场的准确性和有效性.     

一个电镀模型的正则边界元分析—Signorini问题的正则边界元近似

本文应用正则边界元方法研究了由一个电镀模型问题导出的Signorini问题的数值解,给出了近似解的误差分析,计算例子表明,用正则边界元方法求解Signorini问题是行之有效的,并具有计算简便、节省计算时间与内存等优点。     

Ⅱ型平面应力裂纹线场的弹塑性精确解

本文采用线场分析方法对理想弹塑性Ⅱ型平面应力裂纹裂纹线附近的应力场及弹塑性边界进行了精确分析。本文完全放弃了小范围屈服条件,探讨了弹塑性边界上弹塑性应力场匹配条件的正确提法,通过将裂纹线附近塑性区应力场的通解(而不是过去采用的特解)与弹性应力场的精确解(而不是通常的裂尖应力强度因子K场)在裂纹线附近的弹塑性边界上匹配,本文得出了塑性区应力场,塑性区长度及弹塑性边界的单位法向量在裂纹线附近的足够精确的表达式。     

Ⅱ型平面动力裂纹线场的弹塑性精确解

本采用线场分析方法对理想弹塑性Ⅱ型平面应力裂纹裂纹线附近的应力场及弹塑性边界进行了精确分析，本完全放弃了小范围屈服条件，探讨了弹塑性边界上弹塑性应力场匹配条件的正确提法，通过将裂纹线附近塑性区应力场的通解（而不是过去采用的特解）与弹性应力场的精确解（而不是通常的裂尖应力强度因子Ｋ场）在裂纹线附近的弹塑性边界上匹配，本得出了塑性区应力场，塑性区长度及弹塑性边界的单位法向量在裂纹线附近的足够精确     

和式极限的积分方法

计算极限是极限理论的重要内容,大多数函数的极限运算问题可用常规的算法及运算法则解决.而无限多项的和式的极限是极限论当中很难求解的,具有一定难度.本文给出了积分在和式极限求解中的若干命题及计算方法.     

和式极限的积分方法

计算极限是极限理论的重要内容，大多数函数的极限运算问题可用常规的算法及运算法则解决.而无限多项的和式的极限是极限论当中很难求解的，具有一定难度.本文给出了积分在和式极限求解中的若干命题及计算方法.     

多边形应力杂交单元的接触算法研究

针对单元较小的情况下,现有的非均匀模型难以得到精确应力场和应力集中现象的问题,采用直接约束法,提出了一种基于多边形应力杂交单元的优化接触算法.多边形应力杂交单元在应力函数的构造以及积分区域划分上的优势,使其能适应复杂的模型边界与材料边界,更易划分网格.根据上述理论研究成果编制了完整的计算程序,算例结果发现该方法能够得到粉末压制过程的宏观非线性力学响应、高精度的应力场以及明显的应力集中现象,为复杂优化问题的求解提供了有效手段.     

固体弹性三维问题统一解

依据三维弹性力学问题的Kelvin解,用三维虚边界元法来建立积分方程,从而使三维实体和各类板、壳等问题的求解思想得到统一.对各类三维问题采用统一的思想建立数学模型,更有利于程序模块化,增强了程序的通用性.另外,建立积分方程时直接引用Kelvin解,而未引入任何其它假设,使该方法的解更偏于实际,且使应用范围拓宽.再者,与边界元直接法相比,该方法的优点在于无需处理奇异积分,且系数阵是对称的.最后给出部分算例,以证明方法的有效性和计算精度.     

Boundary element method for dynamic inelastic analysis

The paper shows formulation and application of the boundary element method (BEM) for dynamic analysis of elastoplastic materials. The initial stress approach is used in the elastoplastic analysis. The mass matrices are computed by the dual reciprocity method (DRBEM). Displacements and stresses are computed by the iterative procedure in each time step. The time dependent problem is solved by the direct integration Houbolt method. The methods presented in the paper are applied to compute displacements and stresses in a crank loaded by dynamic forces. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

薄板的局部Petrov-Galerkin方法

利用薄板控制微分方程的等效积分对称弱形式和对变量(挠度)采用移动最小二乘近似函数进行插值,研究了薄板弯曲问题的无网格局部Petrov-Galerkin方法．这是一种真正的无网格方法,它不需要任何有限元或边界元网格,不管这种网格是用于能量积分还是进行插值的目的．所有的积分都在规则形状的子域及其边界上进行,并用罚因子法施加本质边界条件．数值例子表明,无网格局部Petrov-Galerkin法不但能够求解二阶微分方程的边值问题,而且求解四阶微分方程的边值问题也很有效,也具有收敛快、稳定性好、对挠度和内力都具有精度高的特点．     

利用边界元法求解一类重调和方程

边界元法(BEM)和多重互易法(MRM)相结合求解一类重调和方程.通过重调和基本解序列给出的MRM-方法和BEM, 推导出该类问题的MRM-边界变分方程, 用边界元法求解该变分方程, 从而得到重调和方程的近似解, 并给出了解的存在唯一性证明.通过数值算例说明了MRM-方法具有收敛速度快、计算精度高, 易编程等优点, 为使用边界元法数值求解重调和方程提供了方法和理论依据.适合于工程中的实际运算.     

Efficient crack growth analyzes by combining fast methods for BEM

The combination of fast methods for the boundary element method (BEM) for efficient crack growth analyzes is presented. Due to the nonlinearity of fatigue crack growth an incremental procedure has to be applied. Within each increment a stress analysis is needed. Based on the asymptotic stress field the stress intensity factors (SIFs) are calculated by an extrapolation method. Then, a new crack front is determined by a reliable 3D crack growth criterion. Finally, the numerical model has to be updated for the next increment. The time dominant factor in each increment is the computation of the stress field. Due to the stress concentration problem the BEM is utilized. To speed-up the calculation several independent fast methods are exploited. An algebraic technique is the adaptive cross approximation (ACA) method which is acting on the system matrix itself. The application of the substructure technique leads to a blockwise band matrix and therefore to reduced memory requirements. Further savings in memory and computation time are reached by modelling cracks with the dual discontinuity method (DDM) and using the ACA method in each substructure. The efficiency of the combined methods is shown by a complex industrial example. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

New Applications of the Boundary Element Method in Dynamics of Solids

The boundary element method (BEM) is developed and applied in new fields of dynamic fracture mechanics, dynamics of composite, elasto–plastic and piezoelectric materials. The BEM results are compared with solutions computed by the finite element method (FEM) showing high accuracy of the BEM. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical solution of the generalized Laplace equation by coupling the boundary element method and the perturbation method

The boundary element method (BEM) has, in general, some advantages with respect to domain methods, in so far as no internal discretization of the domain is required. Such an advantage, however, is lost if the BEM is used directly for the numerical solution of the generalized Laplace equation (GLE). This paper demonstrates that the GLE can also be dealt with advantageously by coupling the BEM with the perturbation method (PM). The technique involves transforming the starting equation (GLE) into an equation without partial first derivatives and then solving a few problems by BEM. The method requires an internal network but the unknowns are only on the boundary. The procedure is applied to numerically solve two test problems with known analytical solutions.     

A coupled BEM and FEM for the interior transmission problem in acoustics

The interior transmission problem (ITP) is a boundary value problem arising in inverse scattering theory, and it has important applications in qualitative methods. In this paper, we propose a coupled boundary element method (BEM) and a finite element method (FEM) for the ITP in two dimensions. The coupling procedure is realized by applying the direct boundary integral equation method to define the so-called Dirichlet-to-Neumann (DtN) mappings. We show the existence of the solution to the ITP for the anisotropic medium. Numerical results are provided to illustrate the accuracy of the coupling method.     

Solution of Poisson’s equation by analytical boundary element integration

The solution of Poisson’s equation is essential for many branches of science and engineering such as fluid-mechanics, geosciences, and electrostatics. Solution of two-dimensional Poisson’s equations is carried out by BEM based on Galerkin Vector Method in which the integrals that appear in the boundary element method are expressed by analytical integration. In this paper, the Galerkin vector method is developed for more general case than presented in the previous works. The integrals are computed for constant and linear elements in BEM. By employing analytical integration in BEM computation, the numerical schemes and coordinate transformations can be avoided. The presented method can also be used for the multiple domain case. The results of the analytical integration are employed in BEM code and the obtained analytical expression will be applied to several examples where the exact solution exists. The produced results are in good agreement with the exact solution.