The determination of an equivalent elastic modulus for bodies with cracks by boundary element method

The equivalent elastic modulus of cracked bodies with orderly distributed cracks was computed with the boundary element method. A practical self-consistent scheme has been proposed in consideration of the mutual interaction effects of the cracks. The influence of friction coefficients and orientation of cracks has been investigated. Some computational examples have been given, and the results show that the proposed method is adequate and the scheme is efficient.This project is supported by the National Natural Science Foundation of China.     

Effective properties of elastic periodic composite media with fibers

A series solution to obtain the effective properties of some elastic composites media having periodically located heterogeneities is described. The method uses the classical expansion along Neuman series of the solution of the periodic elasticity problem in Fourier space, based on the Green's tensor, and exact expressions of factors depending on the shape of the inclusions. Some properties of convergence of the solution are presented, more specifically concerning the elasticity tensor of the reference medium, showing that the convergence occurs even for empty fibers. The solution is extended for rigid inclusions. A comparison is made with previous exact solutions for a fiber composite made of cylindrical fibers with circular cross-sections and with previous estimates. Different examples are presented for new situations concerning the study of fiber composites: composites with elliptic cross-sections and multi-phase fibrous composites.     

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.     

A boundary element formulation using the simple method

A new boundary element method is described for calculation of the steady incompressible laminar flows. The method is based on the well-known SIMPLE algorithm. The new boundary element method allows one to find the fields of the pressure and velocity corrections without inner iterations, thus reducing the computational time drastically. This makes it different from the method developed by Patankar and Spalding.³² However, the new method demands a much larger computer strorage. The boundary integral equations are discretized with the help of constant boundary elements and constant cells. The values of the integrals along the boundary elements and the cells for the two-dimensional domain are found analytically. To preserve the stability in the iteration process, under-relaxation for the convection terms is used. This paper gives the results of calculations of the flows between two plane parallel plates at Re = 20 and Re = 200, the flows in a square cavity with a moving upper lid at Re = 1 and Re = 100 and the flow in a plane channel with sudden symmetric expansion at Re =46·6.     

The 3-D boundary element method of roller bearing by plate element analogue

The mill roller bearing is made up of an internal ring, middle rolls and an external ring, the analysis of which is a multi-bodies contact problem. In this paper, based on the three-dimensional elastic contact BEM without friction, and using the structural characteristics of roller bearings, middle rolls are described by elastic plate units of different shapes, which is placed on the internal ring. The discontinuous traction problem of the contact element can be dealt with by the traction sub-element method. Therefore the contact problem can be changed from one of multi-bodies to one of two-bodies. As an elastic deformation of the plate element, it can be calculated by an elastic contact formula and is regarded as an equivalent gap to be put into the total matrix equality. The face-to-face contact condition of displacement and traction on the element will serve as the convergence judging condition. So convergence solution can be obtained before any geometry inconsistent on the contact element may take place. Therefore, the three-dimensional distribution of load of the mill roller bearing can be achieved. It has been proved that the calculating model of this method is visualized and simple, has fast convergence solution and higher accuracy, making it an effective numerical method of designing and analyzing load characteristics of the mill roller bearing. Project supported by the National Natural Science Foundation of China (50075075).     

不同材料契合弹性力学问题的虚边界元最小二乘法

本文针对不同材料契合弹性力学问题,由虚边界元方法出发,建立了引入拉氏乘子的最小二科解法,对该类问题避免了采用Ｈｅｔｅｎｙｉ’ｓ基本解的麻烦和限制。数值结果表明本文算法的有效性和计算精度高。     

Inverse doubly periodic problem of the theory of bending of a plate with elastic inclusions

Based on the balanced strength principle, a problem of determining the optimal interference for fitting elastic inclusions into holes of an isotropic elastic plate weakened by a doubly periodic system of circular holes is solved. A closed system of algebraic equations is derived, which allows solving this problem. The resultant interference increases the load-carrying capacity of the composite plate being bent. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 4, pp. 153–161, July–August, 2006.     

The elastic solution of concentrated force acting in orthogonal anisotropic half-plane and constant element fundamental formulae of boundary element method

In this paper, the elastic solutions of concentrated force acting in orthogonal anisotropic half-plane are derived by imaginal method and the formulae of coefficient matrix for constant element are put forward. To solve half-plane problems numerically by BEM, this paper provides the necessary formulae. Because the expressions of fundamental solutions are very simple, the object functions could be obtained for every integral of constant element and higher order element of indirect BEM. Thus, the procedure of integration could be avoided in calculation program.     

A boundary element method for anisotropic coupled thermoelasticity

Summary A boundary element formulation is presented for the solution of the equations of fully coupled thermoelasticity for materials of arbitrary degree of anisotropy. By employing the fundamental solutions of anisotropic elastostatics and stationary heat conduction, a system of equations with time-independent matrices is obtained. Since the fundamental solutions are uncoupled and time-independent, a domain integral remains in the representation formula which contains the time-dependence as well as the thermoelastic coupling. This domain integral is transformed to the boundary by means of the dual reciprocity method. By taking this approach, the use of dynamic fundamental solutions is avoided, which enables an efficient calculation of system matrices. In addition, the solution of transient processes as well as, free and forced vibration analysis becomes straightforward and can be carried out with standard time-stepping schemes and eigensystem solvers. Another important advantage of the present formulation is its versatility, since it includes a number of simplified thermoelastic theories, viz. the theory of thermal stresses, coupled and uncoupled quasi-static thermoelasticity, and stationary thermoelasticity. The accuracy of the new thermoelastic boundary element method is demonstrated by a number of example problems. Support by the Deutsche Forschungsgemeinschaft (DFG) of the Graduate Collegium Modelling and discretization methods for continua and fluids (GKKS) at the University of Stuttgart is gratefully acknowledged.     

虚边界元法中的旋转周期对称结构

本文证明具有旋转周期对称性的结构,在对称适应的指标架下,其虚边界元方程的系数矩阵具有块循环的形式,给出一种分解算法（即将原问题分解成一系列相互独立的子问题的求解法）,适于任意形式的载荷分布。     

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.     

基于边界元法求解三维弹性摩擦接触问题

边界元法求解三维摩擦接触问题,其中一个关键点在于如何确定滑移方向。即当出现相对滑移时,滑移方向如何确定。当前常采用的方法是,粘结点利用切向面力得到滑移方向,滑移点利用切向相对位移得到滑移方向。不过该方法难以保证收敛性。针对这一问题,本文采用滑移方向预测技术得到滑移方向。即以后出现相对滑移时,滑移方向采用预测技术中得到的滑移方向。由于摩擦接触问题和历史加载相关,本文采用增量法求解。不同摩擦系数下的数值结果都证明了本文算法的有效性和收敛性及滑移方向预测技术的有效性。     

A boundary element method for calculating elastic waves scattered by axisymmetric inclusion

A Boundary Element Method (BEM) is described to compute the scattering of elastic waves by an axisymmetric inclusion in an infinite elastic medium. The boundary loads applied to the inclusion is expanded in terms of Fourier series in an infinite space. The boundary integral equation is solved in the general direction of the axisymmetric inclusion by BEM. The problem of the 3-D scattering of elastic waves is reduced to a 1-Done. According to the geometric features of the axisymmetric in clusion the ring shell elements are adopted in this method. A comparison is made with other BEM methods. The numerical results show this method can reduce the amount of calculation and enhance the speed of convergence. Supported by Foundation of Ph. D Program of State Education Commission of China     

The numerical solution of Stokes flow in a domain with re-entrant boundaries by the boundary element method

Numerical solutions are presented for two-dimensional low Reynolds number flow in a rotating tank with stationary barriers. The boundary element method is employed, assuming straight panels and quadratic source distribution. The feasibility of repositioning the nodes as a way to minimize the error is explored. A stretching parameter places smaller elements near the re-entrant regions. Elementary error analysis shows uniform improvement in the solution with stretching. The changing eddy pattern for different numbers and sizes of the barriers is compared with experimental results.     

Effective elastic moduli of periodic granular composite with transversely isotropic phases

We find a rigorous solution describing the macroscopically uniform stress state of a periodic granular composite with transversely isotropic phases. The structure of the composite is modeled by a cube containing a finite number of arbitrarily arranged and oriented, transversely isotropic spherical inclusions. This provides the model with a flexible means of describing the microstructure. Applying periodic vector solutions and local expansion formulas reduces the initial boundary-value problem to a system of linear algebraic equations. By averaging the solution over the unit cell, we derived exact finite expressions for the components of the effective stiffness tensor. The numerical data presented help to evaluate the efficiency of the method and the limits of applicability of available approximate theories.Translated from Prikladnaya Mekhanika, Vol. 40, No. 9, pp. 123–130, September 2004.     

Analysis of convective diffusion problem with first-order chemical reaction by boundary element method

A boundary element method is presented for a steady state convective diffusion problem with a first-order chemical reaction. In the chemical reaction system the usual conservation law for calculation of the diagonal components of the matrix H cannot be satisfied. Thus in the boundary element formulation the emphasis is laid on how to compute these diagonal components. The main aim of this paper is to present the method of numerical quadrature to obtain the principal values, using the splitting technique to treat the singular integration.     

A boundary element modeling of fatigue crack growth in a plane elastic plate

This paper presents an extension of a boundary element method to fatigue growth analysis of mixed-mode cracked plane elastic bodies. The method consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity element due to the author. In the boundary element implementation the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the modified maximum strain energy density criterion. In numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the boundary element method. Crack growth is simulated by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characters of some related elements are adjusted according to the manner in which the boundary element method is implemented. Some numerical results of fatigue growth in a plane elastic plate with a center-inclined crack under uniaxial cyclic loading are given.     

Numerical analysis of effective elastic properties of geomaterials containing voids using 3D random fields and finite elements

The purpose of the study is to investigate the influence of porosity and void size on effective elastic geotechnical engineering properties with a 3D model of random fields and finite element. The random field theory is used to generate models of geomaterials containing spatially random voids with controlled porosity and void size. A “tied freedom” analysis is developed to evaluate the effective Young’s modulus and Poisson’s ratio in an ideal block material of finite elements. To deliver a mean and standard deviation of the elastic parameters, this approach uses Monte-Carlo simulations and finite elements, where each simulation leads to an effective value of the property under investigation. The results are extended to investigate an influence of representative volume element (RVE). A comparison of the effective elastic stiffness of 2D and 3D models is also discussed.     

谱元法研究桁架结构的振动带隙特性

采用谱元法研究了桁架周期结构的带隙特性.从杆和梁的运动方程出发,推导出与频率相关的插值函数,得到了杆单元和梁单元的动力学刚度矩阵.在频域下将铰结构考虑为一个谱单元,并推导出铰单元的动力学刚度矩阵.将杆/梁单元和铰单元加以整合得到整体结构的动力学刚度阵,进而建立整体结构的运动方程.通过求解整体结构的动力学方程,获得结构的频域响应,进而研究结构的带隙特性.本文将谱元法求解得到的固有频率结果与有限元法进行了对比,分析了单胞数量和材料的变化对结构带隙特性的影响,拓展了谱元法的应用领域.     

Numerical solution of axisymmetric cavity flows using the boundary element method

This paper presents a formulation of the boundary element method (BEM) for solution of axisymmetric cavity flow problems. The governing equation is written in terms of Stokes' stream function, requiring a new fundamental solution to be found. The iterative procedure for adjusting the free-surface position is similar to that used for planar cavity flows. Numerical results are compared with finite difference and finite element solutions, showing the robustness of the BEM model.