Die design by a boundary element method

A simple technique to design extrusion dies for three dimensional profiles based on Boundary Element Methods is reported. The technique is applied to design a few dies for triangular and square Newtonian extrudates. The results obtained compare well with the available data on common design practice.     

边界元理论在复杂外边界油藏水平井渗流中的应用

受构造作用的影响,实际油藏的外边界往往是复杂多样的.本文从渗流理论出发建立了复杂外边界油藏水平井渗流数学模型,并采用Lord Kelvin点源解、贝塞尔函数积分和泊松叠加公式等方法求解了复杂外边界油藏水平井的边界元基本解,利用边界元的理论建立了复杂外边界油藏水平井井底压力响应数学模型.通过计算得到了无因次压力和压力导数双对数理论图版,并在其基础上分析了复杂外边界油藏水平井渗流特征及其影响因素.     

Optimum shape and topology design using the boundary element method

A procedure is developed for simultaneous shape and topology design optimization of linear elastic two-dimensional continuum structures. An intuitive approach is presented to treat such topological problems whereby material is eliminated from within the structure (by introducing holes at regions of low stress) through a sequence of shape optimization processes. A mathematical programming technique coupled with the boundary element (BE) method of response and sensitivity analyses enables the optimal positioning of these holes plus optimization of the overall structural shape. The analytical derivative BE formulation is explained together with the use of appropriate design velocity fields, and example problems are solved to demonstrate the optimization procedure.     

Application of the Dorodnitsyn finite element method to swirling boundary layer flow

The Dorodnitsyn finite element method for turbulent boundary layer flow with surface mass transfer is extended to include axisymmetric swirling internal boundary layer flow. Turbulence effects are represented by the two-layer eddy viscosity model of Cebeci and Smith¹ with extensions to allow for the effect of swirl. The method is applied to duct entry flow and a 10 degree included-angle conical diffuser, and produces results in close agreement with experimental measurements with only 11 grid points across the boundary layer. The introduction of swirl (w_e/u_e = 0.4) is found to have little effect on the axial skin friction in either a slightly favourable or adverse pressure gradient, but does cause an increase in the displacement area for an adverse pressure gradient. Surface mass transfer (blowing or suction) causes a substantial reduction (blowing) in axial skin friction and an increase in the displacement area. Both suction and the adverse pressure gradient have little influence on the circumferential velocity and shear stress components. Consequently in an adverse pressure gradient the flow direction adjacent to the wall is expected to approach the circumferential direction at some downstream location.     

Computer-aided spillway design using the boundary element method and non-linear programming

The general context of this paper is to support the design of spillways by a direct mathematical approach instead of trial-and-error methods. First, a two-dimensional model is formulated to determine the free surface and the discharge for a stationary, incompressible, homogeneous, non-viscous and irrotational flow over a fixed spillway. The flow satisfies the Laplace equation and the Bernoulli equation (potential flow). An important feature of the model is that it can be extended to design the spillway structure when the spillway is not fixed but the pressure on the spillway is described by a cavitation criterion. Next, the continuous model is discretized by the boundary element method (BEM). We use a non-linear programming algorithm to calculate the pressures and the shape of the spillway. A computer-aided design package is developed on a PC using the equations describing the free surface, the BEM and standard optimization techniques. The input and output of the model are realized using graphical routines. Finally, we discuss the convergence and the computation time of the algorithms.     

滚动体几何相似条件在轴承边界元法中的应用

滚动轴承作为多物体接触系统,具有高度的非线性。为了求解滚动轴承的三维载荷分布,考虑所有滚动体的形状和边界条件相同,利用几何相似条件通过初始位置滚动体的离散模型,得到其他滚动体的离散数据。将所有滚动体看作是一个物体,轴承系统可以用3个物体无摩擦接触边界元法进行模拟计算。考虑全部滚动体作为一个整体的特殊性,叙述了边界积分方程和耦合的矩阵方程的建立,编制了基于滚动体几何相似条件的轴承边界元法Fortran源程序。利用该程序对轧机四列圆锥滚子轴承进行数值模拟,得到轴承的接触压力和载荷,滚动体接触宽度和接触滚动体个数,并说明了方法的有效性。     

Site of ductile fracture initiation in cold forging: A finite element model

Ductile fracture occurs due to micro-void nucleation, growth and finally coalescence into micro-crack. In this study a new ductile fracture condition that based on the microscopic phenomena of void nucleation, growth and coalescence was proposed. Using this condition and combining with finite element model to predict the fracture locations in bulk metal forming, the results show that it is in close accordance with observations of some experimental specimens. However, in order to obtaining the high trustiness many experiments have to be carried out.     

Wideband fast multipole boundary element method: Application to acoustic scattering from aerodynamic bodies

A wideband adaptive multi‐level fast multipole method (MLFMM) is used to accelerate the matrix–vector products arising from a boundary element method (BEM) formulation which solves the Burton–Miller boundary integral equation (BIE). The wideband MLFMM presented here applies a plane wave expansion formulation with fast interpolation and filtering for calculations in the high‐frequency regime and a partial wave expansion formulation with rotation‐coaxial translation in the low‐frequency regime. The iterative solvers GMRES, Bi‐CGSTAB and CGS are tested and compared and a block diagonal preconditioner is used to improve the condition number of the BEM matrices and to accelerate the convergence of the iterative solvers. Details on the implementation of the formulations are described, including the treatment of singular integrals. Results for acoustic scattering from a wing plus engine nacelle configuration for a prescribed source in a subsonic uniform flow are presented for a broad range of frequencies in order to assess the implemented capability. Copyright © 2010 John Wiley & Sons, Ltd.     

一种新型的求解扭转问题的边界元方法

证明了柱体自由扭转的边界积分方程被积函数的散度等于零,将翘曲函数表示为翘曲势函数在边界点的数值计算,避免求解奇异的数值积分。实例计算表明,该表度较高。     

Stochastic boundary element method in elasticity

The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and randomly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations concerning the mean and deviation of the displacements are derived, respectively. It is found that the randomness of material parameters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared with those obtained by theoretical analysis or other numerical methods. The project supported by the National Natural Science Foundation of China and the State Education Commission Foundation of China     

A symmetrical nonsingular boundary element method

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｏｕｇｈｓｉｎｇｕｌａｒＢｏｕｎｄａｒｙｅｌｅｍｅｎｔｉｓｒｅｇａｒｄｅｄａｓａｋｉｎｄｏｆｗｅｉｇｈｔｅｄｒｅｓｉｄｕａｌｓ,ｔｈｅｓｉｎｇｕｌａｒｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｈａｓａｕｎｉｑｕｅｃｈａｒａｃｔｅｒｉｓ...     

Damage analysis by boundary element method

The response of cracked bodies subjected to loading was investigated by the boundary element method in this paper. The two-law elastic-cohesive-softening model was used for crack propagation analysis. The interface conditions for uncracked, craze, open crack, adhesive crack and slid crack parts were discussed and the corresponding incremental iteration algorithm was given. A simplified damage propagation model was presented. The technique has been applied to some specific examples which give the evidence that the method is satisfactory and efficient.     

Application of an improved boundary element method on the problem of elastic torsion

An improved boundary clement method has been used in analyzing and calculating the problems of the torsion of a prismatic bar with elliptical cross-section. In this paper the calculated results correspond with the values of boundary element method. However, the quantity of data required by the improved boundary element method is much less than that required by boundary element method, and the calculating time will be greatly reduced. Therefore, the procedure of this paper is an economical and efficient numerical computational way for solving Poisson equation problem.     

基于单元子分法的结构多尺度边界单元法

建立在基于单元子分法的一种有效自适应格式以及多区域边界元三步求解技术基础上提出了一种计算结构多尺度问题的多区域边界元法。首先,通过高斯积分误差分析公式确定边界单元在满足精度要求下所需要的高斯点数,当所需高斯点数超过规定数目时该单元就被自动划分成一定数量的子单元,从而消除结构多尺度所引起的近奇异性。在单元子分技术的基础上采用多区域边界元三步求解技术来处理材料非均质问题:第一步消除各子域的内部未知量,第二步消除各子域独自拥有的边界未知量,第三步根据位移相容性条件和面力平衡条件建立系统方程组并求解公共界面节点位移以及每个子域的其他未知量。数值算例结果表明本方法可以用较少的计算时间得到满意的结果,是处理结构多尺度问题的一种有效方法。     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

Finite-part integral and boundary element method to solve flat crack problems

Using the Somigliana formula and the concepts of finite-part integral, a set of hypersingular integral equations to solve the arbitrary fiat crack in three-dimensional elasticity is derived and its numerical method is then proposed by combining the finitepart integral method with boundary element method. In order to verify the method, several numerical examples are carried out. The results of the displacement discontinuities of the crack surface and the stress intensity factors at the crack front are in good agrernent with the' theoretical solutions.     

Infinite spline boundary element method and its application to structure-foundation interaction

In view of the infinity behaviors of3-D Kelvin solution,we constructed an infinite spline boundary element which has fine precision in the analysis of the half space foundation subjected to uniform pressure on the circular domain.We also analysed a square plate resting on elastic half space foundation.The results indicate that this model not only fits for the coupled analysis of foundation and structures but also has the advantage of fewer degrees of freedom and fine precision.     

The boundary element method applied to elastic contact problems with friction

In the current paper the boundary integral equations (BIE) for elastic contact problems with friction are derived from the incremental virtual work principle. After introducing contact conditions of adhesion and slip into BIE all variants on boundary are made to discretize by quadratic isoparametric boundary element. In the current paper not only an auto-increment loading law is presented but also the iterative calculation laws for open, slip and adhesion condition are given. The results of numerical examples are satisfactory.