Contact position controlling for two-dimensional motion bodies by the boundary element method

An algorithm is presented for controlling two-dimensional motion contact bodies with conforming discretization. Since a kind of special boundary element is utilized in the algorithm, the displacement compatibility and traction equilibrium conditions at nodes can be satisfied simultaneously in arbitrary locations of the contact interface. In addition, a method is also proposed in which the contact boundary location can be moved flexibly on the possible contact boundary. This method is effective to deal with moving and rolling contact problems on a possible larger moving or rolling contact region. Numerical examples show effectiveness of the presented scheme.     

Multipole BEM for 3-D Elasto-Plastic Contact with Friction

The analysis of 3-D elasto-plastic contact with friction is a highly nonlinear problem. The elements in the contact and plastic zones should be refined to obtain accurate information about the real size,displacement, and traction in the contact zone. However, the increase in the number of degrees of freedom is limited when traditional boundary element method (BEM) is used with the larger memory size and long CPU time required for the solution procedure. This paper describes the additional mathematical friction model to the 3-D elastic multipole BEM to develop a 3-D elasto-plastic contact multipole BEM. Numerical tests show that with this new method, the needed computer memory size is only 2% of the traditional BEM model with friction, which erases large-scale computing with refined meshes and improves the computational accuracy.     

FM-BEM Evaluation for Effective Elastic Moduli of Microcracked Solids

A fast multipole boundary element method (FM-BEM) was applied for the analysis of microcracked solids. Both the computational complexity and memory requirement are reduced to O(N), where N is the number of degrees of freedom. The effective elastic moduli of a 2-D solid containing thousands of randomly distributed microcracks were evaluated using the FM-BEM. The results prove that both the differential method and the method proposed by Feng and Yu provide satisfactory estimates to such problems. The effect of a non-uniform distribution of microcracks has been studied using a novel model. The numerical results show that the non-uniform distribution induces a small increase in the global stiffness.     

Mathematical Programming Solution for the Frictional Contact Muitipole BEM

This paper presents a new mathematical model for the highly nonlinear problem of frictional contact. A programming model, multipole boundary element method (BEM), was developed for 3-D elastic contact with friction to replace the Monte Carlo method. A numerical example shows that the optimization programming model for the point-to-surface contact with friction and the fast optimization generalized minimal residual algorithm (GMRES(m)) significantly improve the analysis of such problems relative to the conventional BEM.     

Shakedown Analysis of 3-D Structures Using the Boundary Element Method Based on the Static Theorem

The static shakedown theorem was reformulated for the boundary element method (BEM) rather than the finite element method with Melan‘s theorem, then used to develop a numerical solution procedure for shakedown analysis. The self-equilibrium stress field was constructed by a linear combination of several basis self-equilibrium stress fields with undetermined parameters. These basis self-equilibrium stress fields were expressed as elastic responses of the body to imposed permanent strains obtained using a 3-D BEM elastic-plastic incremental analysis. The lower bound for the shakedown load was obtained from a series of nonlinear mathematical programming problems solved using the Complex method. Numerical examples verified the precision of the present method.     

Micromechanics-BEM Analysis for Piezoelectric Composites

The effective material properties of piezoelectric composites are predicted using micromechanics models of the composite structure combined with a boundary element method (BEM) solution of the governing equation. The composites consist of inclusion and matrix phases. The micromechanics method gives formulae for the overall material constants as functions of the concentration matrix, while the boundary element simulation gives numerical solutions of the boundary displacement and electric potential equations for inclusion or hole problems. Numerical results for a piezoelectric plate with circular inclusions are presented to illustrate applications of the proposed micromechanics-BEM formulation.     

Boundary Integral Equations and A Posteriori Error Estimates

Adaptive methods have been rapidly developed and applied in many fields of scientific and engineering computing, Reliable and efficient a posteriori error estimates play key roles for both adaptive finite element and boundary element methods. The aim of this paper is to develop a posteriori error estimates for boundary element methods. The standard a posteriori error estimates for boundary element methods are obtained from the classical boundary integral equations. This paper presents hyper-singular a posteriori error estimates based on the hyper-singular integral equations, Three kinds of residuals are used as the estimates for boundary element errors. The theoretical analysis and numerical examples show that the hypersingular residuals are good a posteriori error indicators in many adaptive boundary element computations.     

An element by element spectral element method for elastic wave modeling

The spectral element method which combines the advantages of spectral method with those of finite element method, provides an efficient tool in simulating elastic wave equation in complex medium. Based on weak form of elastodynamic equations, mathematical formulations for Legendre spectral element method are presented. The wave field on an element is discretized using high?order Lagrange interpolation, and integration over the element is accomplished based upon the Gauss?Lobatto?Legendre integration rule. This results in a diagonal mass matrix which leads to a greatly simplified algorithm. In addition, the element by element technique is introduced in our method to reduce the memory sizes and improve the computation efficiency. Finally, some numerical examples are presented to demonstrate the spectral accuracy and the efficiency. Because of combinations of the finite element scheme and spectral algorithms, this method can be used for complex models, including free surface boundaries and strong heterogeneity.     

A fast SVM training algorithm based on the set segmentation and k-means clustering

At present, studies on training algorithms for support vector machines (SVM) are important issues in the field of machine learning. It is a challenging task to improve the efficiency of the algorithm without reducing the generalization performance of SVM. To face this challenge, a new SVM training algorithm based on the set segmentation and k-means clustering is presented in this paper. The new idea is to divide all the original training data into many subsets, followed by clustering each subset using k-means clustering and finally train SVM using the new data set obtained from clustering centroids. Considering that the decomposition algorithm such as SVMlight is one of the major methods for solving support vector machines, the SVMlight is used in our experiments. Simulations on different types of problems show that the proposed method can solve efficiently not only large linear classification problems but also large nonlinear ones.     

An Adaptive Boundary Collocation Method for Plate Bending Problems

A boundary collocation method based on the least-square technique and a corresponding adaptive computation process have been developed for the plate bending problem. The trial functions are constructed using a series of the biharmonic polynomials, and the local error indicators are given by the residu- als of the energy density on the boundary. In comparison with the conventional collocation methods, the solution accuracy in the present method can be improved in an economical and efficient way. In order to dem- onstrate the efficiency and advantages of the adaptive boundary collocation method proposed in this paper, two numerical examples are presented for circular plates subjected to uniform loads and restrained by mixed boundary conditions. The numerical results for the examples show good agreement with ones presented in the literature.     

显式有限元侵蚀接触搜索算法研究

侵蚀接触是弹体穿甲/侵彻等冲击动力学问题数值模拟与显式有限元程序的核心算法.文中基于接触界面的点-面离散模型,发展了一种有效、精确的侵蚀接触搜索算法.新算法通过引入历史状态量跟踪从点对主片的穿透历史,并结合接触强度方法,克服从点对主片穿透定义的歧义性,精确确定出从点-主片接触对及其相应的接触方向;对于接触点计算,建立了一种非迭代的计算方法,提高了接触计算的稳定性.文中算法基于自主研发的大型并行显式有限元程序PANDA-Impact实现,并利用典型算例进行了验证.数值模拟结果表明,文中建立的侵蚀接触搜索算法可以用于复杂问题模拟并具有较高的计算精度.      

三维热传导问题的间接边界元法

以三维无限域内单位强度点热源作用的温度场为基本解，用二维虚拟热源法求解三维热传导问题的稳态温度场。由于间接边界元方法的数值离散特性，应用二次非协调单元简便可行地解决了数值计算中的角域问题和边界点法向取值不确定问题。对于含ｒ－１的奇异积分取其柯西主值积分，含ｒ－２的奇异积分由间接方法求得。算例表明此方法有很高的精度。为碾压混凝土拱坝仿真计算提供了有力工具。     

应用边界元法解决接触界面附近带有孔洞的2维移动接触问题

为利用边界元法在解决边界接触这类固有的非线性问题上的优势 ,给出了一种处理 2维移动接触问题的方法。由于在该方法中引入了一种特殊的边界单元 ,使得当接触面移动时接触物体之间的位移协调性和面力平衡条件能够同时得到满足。另外还给出了一种在可能接触边界上可灵活地移动接触边界位置的方法 ,该方法对于处理较大范围内的移动接触问题是有效的     

直齿圆柱齿轮三维接触应力边界元分析

提供了一种有摩擦的三维弹性接触问题边界元分析方法及其前处理技术，并将此法用于一对处于啮合状态的直齿轮，得到了满意的结果。该程序使用方便，只需输入齿轮的基本参数，即可得到齿轮不同啮合位置的网格、应力分布。     

三维Helmholtz方程外边值问题的虚边界元法

基于位势的延拓,推导出三维虚边界积分方程.通过选择不同的虚边界,避免相应内问题的特征值与波数重合,从而保证解的唯一性.数值算例验证了该方法求解任意波数三维Helmholtz方程外边值问题的有效性.     

三维Helmholtz方程的边界点方法

将移动最小二乘近似和边界积分方程相结合,提出了求解三维Helmholtz方程内外边值问题的无网格边界点方法.该方法用单层位势理论将Helmholtz方程转化为间接边界积分方程,并用边界点法离散间接边界积分方程.由于边界积分方程中含有基本解的积分计算时会出现弱奇异,详细推导了弱奇异积分的计算方式.数值算例表明了间接边界点法求解三维Helmholtz方程的有效性.     

三维快速多极边界元高性能并行计算

该文实现了快速多极边界元法的一种高性能并行计算。其并行求解器基于自适应新版本快速多极边界元算法,采用三维二次等参元和等精度积分格式,并通过实测的任务量进行分布式并行环境下的合理负载划分。数值算例表明,该求解器在保持高次边界元高精度优点的基础上,对于几何形状不规则的结构仍能保持较好的并行效率,和传统边界元法相比使解题规模有了数量级的提高。这种并行计算为边界元法在大规模复杂工程问题中的应用提供了有效方案。     

混合有限元法在弹性接触问题中的应用

讨论了利用Hellinger-Reissner二类变量广义变分原理解决定边界接触问题的混合有限元法,并且给出这种计算方法的数值实例。数值结果表明该方法在定边界弹性接触问题中的应用行之有效。     

焊接熔池的边界积分方程计算模型

本文提出了准稳态焊接熔池的三维边界积分方程计算模型,介绍了熔池作用力和边界条件的处理方法及计算模型的求解步骤。该计算模型是用边界单元法进行熔池现象计算机模拟的基础。