OVERHAUSER TRIANGULAR ELEMENTS FOR THREE-DIMENSIONAL POTENTIAL PROBLEMS USING BOUNDARY ELEMENT METHODS

Interpolation to boundary data and one-dimensional Overhauser parabola blending methods are used to derive Overhauser triangular elements. The elements are C¹-continuous at inter-element nodes and no functional derivatives are required as nodal parameters. These efficient parametric representation elements are used to solve three-dimensional potential problems using the Boundary Element Method (BEM). Results obtained are generally as accurate as those obtained using Overhauser quadrilateral elements.     

C2-CONTINUOUS ELEMENTS FOR BOUNDARY ELEMENT ANALYSIS

Overhauser's original idea of linearly blending two sets of quadratic C⁰-continuous basis functions to produce a set of C¹-continuous basis functions is employed by linearly blending two sets of quadratic C¹-continuous basis functions. The result is a set of eight basis functions which are C²-continuous from element to element and can be used for boundary element analysis where post-processing of the solution is required. Solutions to Laplace's equation in simple geometries are used to demonstrate the accuracy of the solutions obtained with the new elements. These new elements also provide more accurate values for the first and second derivatives of the solution in the tangential direction at all points on a smooth boundary and good approximations at corner points. © 1997 by John Wiley & Sons, Ltd.     

AN INDIRECT BOUNDARY ELEMENT METHOD FOR PLATE BENDING ANALYSIS

An indirect Boundary Element Method is employed for the static analysis of homogeneous isotropic and linear elastic Kirchhoff plates of an arbitrary geometry. The objectives of this paper consists of a construction and a study of the resulting boundary integral equations as well as a development of stable powerful algorithms for their numerical approximation. These equations involve integrals with high-order kernel singularities. The treatment of singular and hypersingular integrals and a construction of solutions in the neighborhood of the irregular points on the boundary are discussed. Numerical examples illustrate the procedure and demonstrate its advantages. © 1997 by John Wiley & Sons, Ltd.     

R-ADAPTIVE BOUNDARY ELEMENT METHOD FOR UNSTEADY FREE-SURFACE FLOW ANALYSIS

An adaptive meshing based on the r-method is developed for two-dimensional unsteady non-linear flows with a free surface. Coupling of a boundary element equation and a weighted residual formulation of the pressure equation on the free surface is employed in solving the wave problems. A mesh optimization scheme is constructed for these two matrix equations. The final mesh distribution on the free surface is determined from the weighted average of these two adaptive meshes. Through numerical analyses of a non-linear sloshing problem and solitary waves in a tank, the influence of time interval of the remeshing and the weight factor for two adaptive meshes on the accuracy and efficiency of the proposed approach is investigated.     

二维正交各向异性结构弹塑性问题的边界元分析

给出了二维正交各向异性结构弹塑性问题的边界元分析方法, 包括相应边界积分方程、内点应力公式、边界元求解格式以及弹塑性应力计算方法。在弹塑性分析中, 引入了Hill-Tsai 屈服准则, 采用初应力法和切向预测径向返回法确定实际应力状态。通过具体算例分析了二维正交各向异性结构的弹塑性应力和塑性区分布情况, 部分数值结果与已有结果进行了比较, 两者基本吻合。结果表明, 本文中给出的边界元法可以有效地用于求解二维正交各向异性结构的弹塑性问题。      

复杂边界条件功能梯度板三维分析的细观元法

发展一种细观分析和宏观计算相结合的计算方法—细观元法。细观元法是使功能梯度构件宏观响应和材料组分的几何、物理、构造参数直接发生关联的分析方法,实现材料细观结构到构件宏观响应的直接过度分析,为解决功能梯度构件宏、细观跨尺度分析提供了一种有力工具。细观元方法不增加结点自由度,却使得功能梯度板件的任意功能梯度变化、各种复杂边界条件得到反映。用细观元法得到了几种具有不同复杂边界条件的功能梯度板件的力学量三维分布形态。     

复合材料边界元素法研究

边界元法是一种迅速发展的数值解法.本文研究了在复合材料层合板中运用平面二次边界元法的有关问题.首先,讨论了将层合板等效为一般正交各向异性单层板的方法;其次,给出了偏轴情况下的基本解,研究了直线型二次边界元偏轴影响系数的解析计算以及域内应力表达式;最后对算例进行了计算,所得结果与相应的解析解和有限元相符甚好.     

SENSITIVITY ANALYSIS FORMULATION FOR THREE-DIMENSIONAL CONDUCTION HEAT TRANSFER WITH COMPLEX GEOMETRIES USING A BOUNDARY ELEMENT METHOD

A special boundary integral formulation had been proposed to analyse many engineering problems of conduction heat transfer in complex three-dimensional geometries (closely spaced surface and circular hole in infinite domain or simple modification of it) by Rezayat and Burton. One example of such geometries is the mold sets in the injection molding process. In this paper, an efficient and accurate approach for the design sensitivity analysis (DSA) is presented for these kinds of problems in the similar complex geometries using the direct differentiation approach (DDA) based on the above special boundary integral formulation. The present approach utilizes the implicit differentiation of the boundary integral equations with respect to the design variables (radii and locations of circular holes) to yield the sensitivity equations. A sample problem (heat transfer of injection molding cooling system) is solved to demonstrate the accuracy of the present sensitivity analysis formulation. Although the techniques introduced here are applied to a particular problem in heat transfer of injection molding cooling system, their potential application is quite broad.     

能量平衡结合有限元数值计算分析天然气管道裂纹稳定扩展问题

本文提出了在天然气管道裂纹稳定扩展问题中，应用能量平衡方法结合有限元数值计算结果来分析计算裂纹驱动力。能量数值计算采用了可模拟动裂纹在管道上扩展的有限元程序ＰＦＲＡＣ，它包括了对未开裂管道和裂纹扩展管道的能量计算。通过分析外力作功和内部能量在裂纹扩展时的变化率，应用能量平衡方法计算了裂纹驱动力，并与在ＰＦＲＡＣ程序中应用节点力释放方法计算的裂纹驱动力的结果进行了比较     

随机边界元法在大坝可靠性分析中的应用

本文综合考虑重力坝受自重，水压力及温度场的影响，应用随机边界元法，结合一阶二次矩法，分析了一大坝的强度可靠度的稳定可靠度，所得结果与随机有限元法所得相符。     

AUTOMATIC ADAPTIVE 3-D FINITE ELEMENT REFINEMENT USING DIFFERENT-ORDER TETRAHEDRAL ELEMENTS

Automatic refinement finite element analyses were carried out employing three different-order tetrahedral solid elements for the solution of 3-D stress analysis problems. Numerical results indicated that the adaptive refinement procedure could eliminate effectively the effect of singularities and the optimal convergence rate was achieved in all the examples tested. The preconditioned conjugate gradient technique was used for the solution of the large system of simultaneous equations. By interpolating the initial guess of the iteration solver from the previous converged solution, the number of iterations needed for the solution is lower than expected. Furthermore, when the mesh density distribution pattern has converged, it became even more efficient and independent of the number of degrees of freedom in the finite element mesh. The relative efficiency of the three different-order tetrahedral elements has also been compared in terms of storage and computational cost needed for achieving a certain accuracy. It is found that although the cubic T20 element can achieve the highest convergence rate, the T10 element is the most competitive and effective element in terms of storage and computational cost needed. © 1997 John Wiley & Sons, Ltd.     

用虚裂纹模型研究混凝土裂缝扩展的边界单元法

本文将虚裂纹模型引入边界元法中,从而推出了用虚裂纹模型研究混凝土裂缝扩展过程的边界元法。由于边界元法在解决断裂问题时的优势,使得该法成为分析混凝土裂缝的有效方法之一。文中利用自编程序,给出了算例。     

横观各向同性土介质与结构物共同作用的有限元与边界元耦合分析

本文采用有限元与边界元耦合方法对上与结构物进行共同作用的分析．其中上部结构采用有限元子结构法，地基上考虑为横观各向同性特性，应用边界元法，通过基础界面上位移连续、面力平衡的条件进行耦合分析．编制了相应的计算程序，对某高层空间剪力培与地基的共同作用进行了计算，与有限元法和实测结果进行了对比分析，得到了满意的结果．说明：有限元和边界元的耦合方法能充分发挥两者的优越性，不失为共同作用分析的一种有效途径．     

STABILITY ANALYSIS AND DESIGN OF TIME-STEPPING SCHEMES FOR GENERAL ELASTODYNAMIC BOUNDARY ELEMENT MODELS

In the literature there is growing evidence of instabilities in standard time-stepping schemes to solve boundary integral elastodynamic models. However, there has been no theory to support scientists and engineers in assessing the stability of their boundary element algorithms or to help them with the design of new, more stable algorithms. In this paper we present a general framework for the analysis of the stability of any time-domain boundary element model. We illustrate how the stability theory can be used to assess the stability of existing boundary element models and how the insight gained from this analysis can be used to design more stable time-stepping schemes. In particular, we describe a new time-stepping procedure that we have developed, which has substantially enhanced stability characteristics and greater accuracy for the same computational effort. The new scheme, which we have called ‘the half-step scheme’, is shown to have substantially improved performance for the displacement discontinuity boundary element method commonly used to model dynamic fracture interaction and propagation. © 1997 by John Wiley & Sons, Ltd.     

膜结构有限元分析中的平面单元与曲面单元的比较

在对膜结构进行有限元模拟时,通常采用两种面单元形式——平面常应变三角形单元和六结点曲面三角形等参元,究竟哪一种单元形式更优一直是存在争议的问题。首先就这两种单元的有限元格式进行探讨,比较了它们在位移模式、几何方程及刚度矩阵等方面的力学本质差异;然后,通过一些典型算例,比较了它们在计算精度和计算效率等方面的优劣。分析结果表明:在有限元格式方面,曲面三角形等参元较之平面常应变三角形单元更为完备,更加符合结构的实际受力特征;在计算精度和计算效率方面,曲面三角形等参元也要优于平面常应变三角形单元;因此,从结构精细化分析的角度来看,宜采用曲面三角形等参元来进行膜结构有限元分析。     

SOLVING TRANSIENT DYNAMIC CRACK PROBLEMS BY THE HYPERSINGULAR BOUNDARY ELEMENT METHOD

Abstract— The subject of hypersingular boundary integral equations is a rapidly developing topic due to the advantages which this kind of formulation offers compared to the standard boundary integral method. The hypersingular formulation is particularly well suited for fracture mechanics problems, where there are important gradients of the stress field and singularities. This formulation for time domain antiplane problems has been recently addressed by the authors and in the present paper, the formulation for time domain plane problems is presented and applied for the first time. A mixed Boundary Element approach based on the standard integral equation and the hypersingular integral equation is developed. The mixed formulation allows for a very simple discretization of the problem, where no subregion is needed. Conforming quadratic elements are used for the crack and the external boundaries. The hypersingular integral equation is used for collocation points within the crack elements, while the standard integral representation is used for the external boundaries. Several examples with different crack geometries are studied to illustrate the possibilities of the method. The Stress Intensity Factor (S.I.F.) is very accurately computed from the crack tip opening displacements along the crack tip element. The results show that the proposed approach for S.I.F. evaluation is simple and produces accurate solutions.     

ELASTOHYDRODYNAMIC LUBRICATION ANALYSIS OF LAYERED LINE CONTACT BY THE BOUNDARY ELEMENT METHOD

This paper presents a numerical routine to compute the contact characteristics of elastomer layered cylinders lubricated by isoviscous liquids. The indentation of the elastic layer is calculated from boundary integral equations which are solved by linear and quadratic boundary element methods for a finite plane model and a circular representation of the junction. The hydrodynamic equation is also transformed into a boundary integral equation and solved by Simpson's rule. Some factors which possibly affect numerical accuracy are examined. Examples for finite plane and circular layer are analysed with reference to parameters for printing press roller contact, in which results are obtained for the indentation, film thickness and liquid pressure, as well as internal stresses through the simultaneous solution of the elasticity and hydrodynamic equations. The results show that high precision is easily achieved and the method is efficient for such layered problems.     

求解三维物体波浪荷载的边界元模型

建立了求解规则波中任意形状三维物体绕射问题的边界元模型。设置辐射面将流场分为内、外场,外场中的绕射势用特征函数展开式解析表达,内场则由边界积分方程求解,所取的Green 函数为简单的1/r 形式。内场边界采用四边形面元和轴对称面元混合布置的划分方式,轴对称面元上奇点强度沿周向为Fourier 级数展开的高次分布,提高了求解三维绕射问题的精度和效率。经与直立圆柱绕射的线性理论解比较,验证了数值方法的可靠性。利用谱分析原理,将规则波中的结果拓展,得到多向不规则波中绕射物体所受波浪力的统计特性。     

激光干涉仪反射镜三维温度场的快速多极边界元分析

温差造成光刻机激光干涉仪反射镜热变形,从而影响光刻机的精度。该文将基于二次单元的快速边界元法用于激光干涉仪反射镜的大规模温度场模拟。不连续单元的引入可以有效处理角点问题;新型快速多极算法用于边界元法的加速求解。建立统一的二次单元多极展开格式以处理混合边界。数值算例分析了快速多极边界元法的计算精度和效率,并和常规算法比较;使用该算法对激光干涉仪反射镜进行了大规模温度场计算,并和有限元法比较。结果表明：基于二次单元的快速多极边界元法可以高精度求解大规模三维传热问题。     

弹性地基板广义边值问题的边界元法

本文利用Ｈａｎｋｅｌ变换导出了弹性地基板弯曲问题的基本解，该基本解对于Ｗｉｎｋｌｅｒ地基、Ｐａｓｔｅｒｎａｋ地基和弹性半空间地基模型具有统一的表达形式。在此基础上，建立了适用于弹性地基板广义边值问题的边界积分方程组，最后文中给出了若干数值算例。