共查询到20条相似文献，搜索用时 15 毫秒
1.
J. Dominguez M. P. Ariza R. Gallego 《International journal for numerical methods in engineering》2000,48(1):111-135
The present paper deals with a boundary element formulation based on the traction elasticity boundary integral equation (potential derivative for Laplace's problem). The hypersingular and strongly singular integrals appearing in the formulation are analytically transformed to yield line and surface integrals which are at most weakly singular. Regularization and analytical transformation of the boundary integrals is done prior to any boundary discretization. The integration process does not require any change of co‐ordinates and the resulting integrals can be numerically evaluated in a simple and efficient way. The formulation presented is completely general and valid for arbitrary shaped open or closed boundaries. Analytical expressions for all the required hypersingular or strongly singular integrals are given in the paper. To fulfil the continuity requirement over the primary density a simple BE discretization strategy is adopted. Continuous elements are used whereas the collocation points are shifted towards the interior of the elements. This paper pretends to contribute to the transformation of hypersingular boundary element formulations as something clear, general and easy to handle similar to in the classical formulation. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
2.
Rafael Gallego Lucia Comino Alberto Ruiz‐Cabello 《International journal for numerical methods in engineering》2006,66(12):1913-1933
In this paper, the material constant sensitivity boundary integral equation is presented, and its numerical solution proposed, based on boundary element techniques. The formulation deals with plane problems with general rectilinear anisotropy. Expressions for the computation of sensitivities for displacements, tractions, strains and stresses are derived, both for boundary and interior points. The sensitivities can be computed with respect to the bulk material properties or to the properties of part of the domain (inclusions, coatings, etc.). To assess the accuracy of the proposed approach, the computed results are compared to analytical ones derived from exact solutions obtained by complex potential theory, when possible, or finite difference derivatives otherwise. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
3.
T. Takahashi A. Kawai T. Ebisuzaki 《International journal for numerical methods in engineering》2006,66(3):529-548
We present a new solution to accelerate the boundary integral equation method (BIEM). The calculation time of the BIEM is dominated by the evaluation of the layer potential in the boundary integral equation. We performed this task using MDGRAPE‐2, a special‐purpose computer designed for molecular dynamics simulations. MDGRAPE‐2 calculates pairwise interactions among particles (e.g. atoms and ions) using hardwired‐pipeline processors. We combined this hardware with an iterative solver. During the iteration process, MDGRAPE‐2 evaluates the layer potential. The rest of the calculation is performed on a conventional PC connected to MDGRAPE‐2. We applied this solution to the Laplace and Helmholtz equations in three dimensions. Numerical tests showed that BIEM is accelerated by a factor of 10–100. Our rather naive solution has a calculation cost of O(N2 × Niter), where N is the number of unknowns and Niter is the number of iterations. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
4.
K. M. Liew Yuzhou Sun S. Kitipornchai 《International journal for numerical methods in engineering》2007,69(4):729-749
This paper considers a 2‐D fracture analysis of anisotropic piezoelectric solids by a boundary element‐free method. A traction boundary integral equation (BIE) that only involves the singular terms of order 1/r is first derived using integration by parts. New variables, namely, the tangential derivative of the extended displacement (the extended displacement density) for the general boundary and the tangential derivative of the extended crack opening displacement (the extended displacement dislocation density), are introduced to the equation so that solution to curved crack problems is possible. This resulted equation can be directly applied to general boundary and crack surface, and no separate treatments are necessary for the upper and lower surfaces of the crack. The extended displacement dislocation densities on the crack surface are expressed as the product of the characteristic terms and unknown weight functions, and the unknown weight functions are modelled using the moving least‐squares (MLS) approximation. The numerical scheme of the boundary element‐free method is established, and an effective numerical procedure is adopted to evaluate the singular integrals. The extended ‘stress intensity factors’ (SIFs) are computed for some selected example problems that contain straight or curved cracks, and good numerical results are obtained. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
5.
W. S. Hwang 《International journal for numerical methods in engineering》1999,44(12):1775-1783
A boundary spectral method is developed to solve acoustical problems with arbitrary boundary conditions. A formulation, originally derived by Burton and Miller, is used to overcome the non‐uniqueness problem in the high wave number range. This formulation is further modified into a globally non‐singular form to simplify the procedure of numerical quadrature when spectral methods are applied. In the present approach, generalized Fourier coefficients are determined instead of local variables at nodes as in conventional methods. The convergence of solutions is estimated through the decay of magnitude of the generalized Fourier coefficients. Several scattering and radiation problems from a sphere are demonstrated with high wave numbers in the present paper. Copyright © 1999 John Wiey & Sons, Ltd. 相似文献
6.
论Helmholtz方程的一类边界积分方程的合理性 总被引：5，自引：0，他引：5
本文导出了Helmholtz 方程超定边值问题有解的一个充要条件,和用非解析开拓法证明了文[1]中的Helmholtz 方程在外域中的解的边界积分表示式的合理性,并将此类边界积分表示式推广用于带空洞的有限域。这样就比较严密而又浅近地证明了基于该表示式建立起来的间接变量和直接变量边界积分方程的合理性。 相似文献
7.
Davoud Mirzaei Mehdi Dehghan 《International journal for numerical methods in engineering》2009,79(13):1662-1682
In this paper the meshless local boundary integral equation (LBIE) method for numerically solving the non‐linear two‐dimensional sine‐Gordon (SG) equation is developed. The method is based on the LBIE with moving least‐squares (MLS) approximation. For the MLS, nodal points spread over the analyzed domain are utilized to approximate the interior and boundary variables. The approximation functions are constructed entirely using a set of scattered nodes, and no element or connectivity of the nodes is needed for either the interpolation or the integration purposes. A time‐stepping method is employed to deal with the time derivative and a simple predictor–corrector scheme is performed to eliminate the non‐linearity. A brief discussion is outlined for numerical integrations in the proposed algorithm. Some examples involving line and ring solitons are demonstrated and the conservation of energy in undamped SG equation is investigated. The final numerical results confirm the ability of method to deal with the unsteady non‐linear problems in large domains. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
8.
边界元法是边界积分方程的数值解法 ,是随着计算机技术的发展而出现的。建立声学边界积分方程分两种方法 :直接法与间接法。本文介绍了边界元法在环境声学中的应用 ,如声屏障和不同情况下道路周围的声场分布、复杂气象条件对声传播的影响的问题等。由于边界元法是半解析半数值解法。在解边界积分方程时会遇到解的存在与唯一性问题。 相似文献
9.
10.
Barbara M. Johnston Peter R. Johnston 《International journal for numerical methods in engineering》2003,56(4):589-607
Accurate numerical evaluation of integrals arising in the boundary element method is fundamental to achieving useful results via this solution technique. In this paper, a number of techniques are considered to evaluate the weakly singular integrals which arise in the solution of Laplace's equation in three dimensions and Poisson's equation in two dimensions. Both are two‐dimensional weakly singular integrals and are evaluated using (in a product fashion) methods which have recently been used for evaluating one‐dimensional weakly singular integrals arising in the boundary element method. The methods used are based on various polynomial transformations of conventional Gaussian quadrature points where the transformation polynomial has zero Jacobian at the singular point. Methods which split the region of integration into sub‐regions are considered as well as non‐splitting methods. In particular, the newly introduced and highly accurate generalized composite subtraction of singularity and non‐linear transformation approach (GSSNT) is applied to various two‐dimensional weakly singular integrals. A study of the different methods reveals complex relationships between transformation orders, position of the singular point, integration kernel and basis function. It is concluded that the GSSNT method gives the best overall results for the two‐dimensional weakly singular integrals studied. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
11.
A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected
to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method
combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations
of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are
applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor
solutions to several typical planar crack problems in a finite body are evaluated.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献
12.
L. S. Miers J. C. F. Telles 《International journal for numerical methods in engineering》2008,76(7):1090-1107
A meshless procedure, based on boundary integral equations, is proposed to analyze elastoplastic problems. To cope with non‐linear problems, the usual boundary element method introduces domain discretization cells, often considered a ‘drawback’ of the method. Here, to get rid of the standard element and cell, i.e. boundary and domain discretization, the orthogonal moving least squares (also known as improved moving least squares) method is used. The algorithm adopted to solve these particular inelastic non‐linear problems is a well‐established, criterion‐independent implicit procedure, previously developed by the authors. Comparative results are presented at the end to illustrate the effectiveness of the proposed techniques. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
13.
W. S. Hwang 《International journal for numerical methods in engineering》1998,41(6):1077-1085
The boundary spectral method for solving three-dimensional non-lifting potential problems is developed. This method combines spectral approximations and the direct numerical integration such as Gaussian quadrature or trapezoidal rules successfully. The singularities of the integral equation are completely removed by subtracting known solutions from the Laplace equation. After discretization, every element of the resultant matrix only contains integrals with non-singular kernels. Therefore, all the integrals can be implemented easily and efficiently. By spectral approximations, the unknown variable is expressed as a truncated series of basis functions, which are orthogonal usually. Instead of solving the variables at collocation points in the conventional methods, the coefficients of basis functions are determined in the spectral approach. It is shown that the new method reduces a lot of number of unknowns, storage of matrix elements, and computer time for solving the algebraic equations. © 1998 John Wiley & Sons, Ltd. 相似文献
14.
EDUARD S. VENTSEL 《International journal for numerical methods in engineering》1997,40(9):1597-1610
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. 相似文献
15.
车内声场的数学模型建立 总被引：1，自引：1，他引：1
本文首先利用Ｈｅｌｍｈｏｌｔｚ方程和Ｇｒｅｅｎ定理推导出适合多种边界条件的车内声场边界积分方程，然后利用边界元数值分析技术离散方程，得到已知某一封闭空间边界的振动特性求解其内部声压的边界元数学模型。作为验证，本文还对两个实例进行了试验，结果表明边界元计算值与理论值和试验实测值吻合良好。 相似文献
16.
A Zener-Stroh curved crack is defined such that the crack undergoes an initial displacement discontinuity. A singular integral equation is suggested to solve the Zener-Stroh curved crack problem. General formulation for evaluating the stress intensity factors and the T-stresses at the crack tips of a Zener-Stroh curved crack is carried out. For the Zener-Stroh arc crack, T-stresses at the crack tips can be evaluated in a closed form. 相似文献
17.
M. Eskandari 《Engineering Fracture Mechanics》2010,77(5):819-3172
The analytical treatment of an energetically consistent annular crack in a piezoelectric solid subjected to remote opening electromechanical loading is addressed. Potential functions and Hankel transform in combination with a robust technique are employed to reduce the solution of the mixed boundary value problem into a Fredholm integral equation of the second kind. The limiting case of a penny-shaped crack in a piezoelectric medium with energetically consistent boundary conditions over the crack faces is extracted for the first time. The electrical discharge phenomenon within the crack gap is modeled utilizing a non-linear constitutive law and the effects of the breakdown field on the energy release rate are delineated. The energy release rate, the electric displacement inside the crack gap, and the closing traction on crack faces are plotted for all possible geometries of a non-discharging annular crack. 相似文献
18.
Mixed-mode crack analysis in unidirectionally and bidirectionally functionally graded materials is performed by using a boundary integral equation method. To make the analysis tractable, the Young's modulus of the functionally graded materials is assumed to be exponentially dependent on spatial variables, while the Poisson's ratio is assumed to be constant. The corresponding boundary value problem is formulated as a set of hypersingular traction boundary integral equations, which are solved numerically by using a Galerkin method. The present method is especially suited for straight cracks in infinite FGMs. Numerical results for the elastostatic stress intensity factors are presented and discussed. Special attention of the analysis is devoted to investigate the effects of the material gradients and the crack orientation on the elastostatic stress intensity factors. 相似文献
19.
A numerical solution for half-ellipse-shaped crack is presented. In the solution, the hypersingular integral equation in conjunction with the curve length technique is used. Particular attention is paid to the case when b/c (b and c the minor and major axis of an ellipse) is a small value, for example b/c=0.02. A reasonable accuracy requires a sufficient discretization. 相似文献
20.
J. Englund 《International journal for numerical methods in engineering》2006,66(11):1791-1816
The stress field in a finite, edge cracked specimen under load is computed using algorithms based on two slightly different integral equations of the second kind. These integral equations are obtained through left regularizations of a first kind integral equation. In numerical experiments it is demonstrated that the stress field can be accurately computed. Highly accurate stress intensity factors and T‐stresses are presented for several setups and extensive comparisons with results from the literature are made. For simple geometries the algorithms presented here achieve relative errors of less than 10?10. It is also shown that the present algorithms can accurately handle both geometries with arbitrarily shaped edge cracks and geometries containing several hundred edge cracks. All computations were performed on an ordinary workstation. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献