共查询到20条相似文献,搜索用时 15 毫秒
1.
Rafael Gallego Javier Surez 《International journal for numerical methods in engineering》2000,49(4):501-518
In this paper a procedure to solve the identification inverse problems for two‐dimensional potential fields is presented. The procedure relies on a boundary integral equation (BIE) for the variations of the potential, flux, and geometry. This equation is a linearization of the regular BIE for small changes in the geometry. The aim in the identification inverse problems is to find an unknown part of the boundary of the domain, usually an internal flaw, using experimental measurements as additional information. In this paper this problem is solved without resorting to a minimization of a functional, but by an iterative algorithm which alternately solves the regular BIE and the variation BIE. The variation of the geometry of the flaw is modelled by a virtual strainfield, which allows for greater flexibility in the shape of the assumed flaw. Several numerical examples demonstrate the effectiveness and reliability of the proposed approach. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
2.
P. D. SHAH C. L. TAN X. WANG 《Fatigue & Fracture of Engineering Materials & Structures》2006,29(5):343-356
The importance of a two‐parameter approach in the fracture mechanics analysis of many cracked components is increasingly being recognized in engineering industry. In addition to the stress intensity factor, the T stress is the second parameter considered in fracture assessments. In this paper, the path‐independent mutual M‐integral method to evaluate the T stress is extended to treat plane, generally anisotropic cracked bodies. It is implemented into the boundary element method for two‐dimensional elasticity. Examples are presented to demonstrate the veracity of the formulations developed and its applicability. The numerical solutions obtained show that material anisotropy can have a significant effect on the T stress for a given cracked geometry. 相似文献
3.
A simple demonstration of the existence of the Cauchy principal value (CPV) of the strongly singular surface integral in the Somigliana Identity at a non-smooth boundary point is presented. First a regularization of the strongly singular integral by analytical integration of the singular term in the radial direction in pre-image planes of smooth surface patches is carried out. Then it is shown that the sum of the angular integrals of the characteristic of the tractions of the Kelvin fundamental solution is zero, a formula for the transformation of angles between the tangent plane of a suface patch and the pre-image plane at smooth mapping of the surface patch being derived for this purpose. 相似文献
4.
A dual integral formulation for the interior problem of the Laplace equation with a smooth boundary is extended to the exterior problem. Two regularized versions are proposed and compared with the interior problem. It is found that an additional free term is present in the second regularized version of the exterior problem. An analytical solution for a benchmark example in ISBE is derived by two methods, conformal mapping and the Poisson integral formula using symbolic software. The potential gradient on the boundary is calculated by using the hypersingular integral equation except on the two singular points where the potential is discontinuous instead of failure in ISBE benchmarks. Based on the matrix relations between the interior and exterior problems, the BEPO2D program for the interior problem can be easily reintegrated. This benchmark example was used to check the validity of the dual integral formulation, and the numerical results match the exact solution well. 相似文献
5.
Marc Bonnet Massimo Guiggiani 《International journal for numerical methods in engineering》1998,41(7):1255-1275
This paper investigates the evaluation of the sensitivity, with respect to tangential perturbations of the singular point, of boundary integrals having either weak or strong singularity. Both scalar potential and elastic problems are considered. A proper definition of the derivative of a strongly singular integral with respect to singular point perturbations should accommodate the concomitant perturbation of the vanishing exclusion neighbourhood involved in the limiting process used in the definition of the integral itself. This is done here by esorting to a shape sensitivity approach, considering a particular class of infinitesimal domain perturbations that ‘move’ individual points, and especially the singular point, but leave the initial domain globally unchanged. This somewhat indirect strategy provides a proper mathematical setting for the analysis. Moreover, the resulting sensitivity expressions apply to arbitrary potential-type integrals with densities only subjected to some regularity requirements at the singular point, and thus are applicable to approximate as well as exact BEM solutions. Quite remarkable is the fact that the analysis is applicable when the singular point is located on an edge and simply continuous elements are used. The hypersingular BIE residual function is found to be equal to the derivative of the strongly singular BIE residual when the same values of the boundary variables are substituted in both SBIE and HBIE formulations, with interesting consequences for some error indicator computation strategies. © 1998 John Wiley & Sons, Ltd. 相似文献
6.
Taku Tada Raúl Madariaga 《International journal for numerical methods in engineering》2001,50(1):227-251
We present an efficient numerical method for solving indirect boundary integral equations that describe the dynamics of a flat two‐dimensional (2‐D) crack in all modes of fracture. The method is based on a piecewise‐constant interpolation, both in space and time, of the slip‐rate function, by which the original equation is reduced to a discrete convolution, in space and time, of the slip‐rate and a set of analytically obtained coefficients. If the time‐step interval is set sufficiently small with respect to the spatial grid size, the discrete equations decouple and can be solved explicitly. This semi‐analytic scheme can be extended to the calculation of the wave field off the crack plane. A necessary condition for the numerical stability of this scheme is investigated by way of an exhaustive set of trial runs for a kinematic problem. For the case investigated, our scheme is very stable for a fairly wide range of control parameters in modes III and I, whereas, in mode II, it is unstable except for some limited ranges of the parameters. The use of Peirce and Siebrits' ε‐scheme in time collocation is found helpful in stabilizing the numerical calculation. Our scheme also allows for variable time steps. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
7.
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. 相似文献
8.
Vitaly E. Bulgakov Marina V. Bulgakova 《International journal for numerical methods in engineering》1998,43(3):533-548
A finite element constructed on the basis of boundary integral equations is proposed. This element has a flexible shape and arbitrary number of nodes. It also has good approximation properties. A procedure of constructing an element stiffness matrix is demonstrated first for one-dimensional case and then for two-dimensional steady-state heat conduction problem. Numerical examples demonstrate applicability and advantages of the method. © 1998 John Wiley & Sons, Ltd. 相似文献
9.
Yan Gu Wen Chen Hongwei Gao Chuanzeng Zhang 《International journal for numerical methods in engineering》2016,107(2):109-126
This study documents the first attempt to extend the singular boundary method, a novel meshless boundary collocation method, for the solution of 3D elasticity problems. The singular boundary method involves a coupling between the regularized BEM and the method of fundamental solutions. The main idea here is to fully inherit the dimensionality and stability advantages of the former and the meshless and integration‐free attributes of the later. This makes it particularly attractive for problems in complex geometries and three dimensions. Four benchmark 3D problems in linear elasticity are well studied to demonstrate the feasibility and accuracy of the proposed method. The advantages, disadvantages, and potential applications of the proposed method, as compared with the FEM, BEM, and method of fundamental solutions, are also examined and discussed. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
10.
Haojie Lian Pierre Kerfriden Stéphane P. A. Bordas 《International journal for numerical methods in engineering》2016,106(12):972-1017
The present work addresses shape sensitivity analysis and optimization in two‐dimensional elasticity with a regularized isogeometric boundary element method (IGABEM). Non‐uniform rational B‐splines are used both for the geometry and the basis functions to discretize the regularized boundary integral equations. With the advantage of tight integration of design and analysis, the application of IGABEM in shape optimization reduces the mesh generation/regeneration burden greatly. The work is distinct from the previous literatures in IGABEM shape optimization mainly in two aspects: (1) the structural and sensitivity analysis takes advantage of the regularized form of the boundary integral equations, eliminating completely the need of evaluating strongly singular integrals and jump terms and their shape derivatives, which were the main implementation difficulty in IGABEM, and (2) although based on the same Computer Aided Design (CAD) model, the mesh for structural and shape sensitivity analysis is separated from the geometrical design mesh, thus achieving a balance between less design variables for efficiency and refined mesh for accuracy. This technique was initially used in isogeometric finite element method and was incorporated into the present IGABEM implementation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
11.
T. A. CRUSE J. D. RICHARDSON 《International journal for numerical methods in engineering》1996,39(19):3273-3304
The paper presents two fully equivalent and regular forms of the hypersingular Somigliana stress identity in elasticity that are appropriate for problems in which the displacement field (and resulting stresses) is C1,α continuous. Each form is found as the result of a single decomposition process on the kernels of the Somigliana stress identity in three dimensions. The results show that the use of a simple stress state for regularization arises in a direct manner from the Somigliana stress identity, just as the use of a constant displacement state regularization arose naturally for the Somigliana displacement identity. The results also show that the same construction leads naturally to a finite part form of the same identity. While various indirect constructions of the equivalents to these findings are published, none of the earlier forms address the fundamental issue of the usual discontinuities of boundary data in the hypersingular Somigliana stress identity that arise at corners and edges. These new findings specifically focus on the corner problem and establish that the previous requirements for continuity on the densities in the hypersingular Somigliana stress identity are replaced by a sole requirement on displacement field continuity. The resulting regularized and finite part forms of the Somigliana stress identity leads to a regularized form of the stress boundary integral equation (stress-BIE). The regularized stress-BIE is shown to properly allow piecewise discontinuity of the boundary data subject only to C1,α continuity of the underlying displacement field. The importance of the findings is in their application to boundary element modeling of the hypersingular problem. The piecewise discontinuity derivation for corners is found to provide a rigorous and non-singular basis for collocation of the discontinuous boundary data for both the regularized and finite part forms of the stress-BIE. The boundary stress solution for both forms is found to be an average of the computed stresses at collocation points at the vertices of boundary element meshes. Collocation at these points is shown to be without any unbounded terms in the formulation thereby eliminating the use of non-conforming elements for the hypersingular equations. The analytical findings in this paper confirm the correct use of both regularized and finite part forms of the stress-BIE that have been the basis of boundary element analysis previously published by the first author of the current paper. 相似文献
12.
Cheng-Hung Huang Chih-Wei Chen 《International journal for numerical methods in engineering》1998,42(5):943-965
A Boundary Element Method (BEM)-based inverse algorithm utilizing the iterative regularization method, i.e. the conjugate gradient method (CGM), is used to solve the Inverse Heat Conduction Problem (IHCP) of estimating the unknown transient boundary temperatures in a multi-dimensional domain with arbitrary geometry. The results obtained by the CGM are compared with that obtained by the standard Regularization Method (RM). The error estimation based on the statistical analysis is derived from the formulation of the RM. A 99 per cent confidence bound is thus obtained. Finally, the effects of the measurement errors to the inverse solutions are discussed. Results show that the advantages of applying the CGM in the inverse calculations lie in that (i) the major difficulties in choosing a suitable quadratic norm, determining a proper regularization order and determining the optimal smoothing (or regularization) coefficient in the RM are avoided and (ii) it is less sensitive to the measurement errors, i.e. more accurate solutions are obtained. © 1998 John Wiley & Sons, Ltd. 相似文献
13.
Jeng-Tzong Chen Wen-Sheng Huang Ying-Te Lee Shyh-Rong Kuo Shing-Kai Kao 《先进材料力学与结构力学》2017,24(1):1-15
The boundary integral equation method in conjunction with the degenerate kernel, the direct searching technique (singular value decomposition), and the only two-trials technique (2 × 2 matrix eigenvalue problem) are analytically and numerically used to find the degenerate scales, respectively. In the continuous system of boundary integral equation, the degenerate kernel for the 2D Kelvin solution in the polar coordinates is reviewed and the degenerate kernel in the elliptical coordinates is derived. Using the degenerate kernel, an analytical solution of the degenerate scales for the elasticity problem of circular and elliptical cases is obtained and compared with the numerical result. Further, the triangular case and square case were also numerically demonstrated. 相似文献
14.
L. Távara J. E. Ortiz V. Mantič F. París 《International journal for numerical methods in engineering》2008,74(5):776-798
A general, efficient and robust boundary element method (BEM) formulation for the numerical solution of three‐dimensional linear elastic problems in transversely isotropic solids is developed in the present work. The BEM formulation is based on the closed‐form real‐variable expressions of the fundamental solution in displacements Uik and in tractions Tik, originated by a unit point force, valid for any combination of material properties and for any orientation of the radius vector between the source and field points. A compact expression of this kind for Uik was introduced by Ting and Lee (Q. J. Mech. Appl. Math. 1997; 50 :407–426) in terms of the Stroh eigenvalues on the oblique plane normal to the radius vector. Working from this expression of Uik, and after a revision of their final formula, a new approach (based on the application of the rotational symmetry of the material) for deducing the derivative kernel Uik, j and the corresponding stress kernel Σijk and traction kernel Tik has been developed in the present work. These expressions of Uik, Uik, j, Σijk and Tik do not suffer from the difficulties of some previous expressions, obtained by other authors in different ways, with complex‐valued functions appearing for some combinations of material parameters and/or with division by zero for the radius vector at the rotational‐symmetry axis. The expressions of Uik, Uik, j, Σijk and Tik have been presented in a form suitable for an efficient computational implementation. The correctness of these expressions and of their implementation in a three‐dimensional collocational BEM code has been tested numerically by solving problems with known analytical solutions for different classes of transversely isotropic materials. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
15.
Yuqing Niu Marijan Dravinski 《International journal for numerical methods in engineering》2003,58(7):979-998
Based on the full‐space Green's functions, a three‐dimensional time‐harmonic boundary element method is presented for the scattering of elastic waves in a triclinic full space. The boundary integral equations for incident, scattered and total wave fields are given. An efficient numerical method is proposed to calculate the free terms for any geometry. The discretization of the boundary integral equation is achieved by using a linear triangular element. Applications are discussed for scattering of elastic waves by a spherical cavity in a 3D triclinic medium. The method has been tested by comparing the numerical results with the existing analytical solutions for an isotropic problem. The results show that, in addition to the frequency of the incident waves, the scattered waves strongly depend on the anisotropy of the media. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
16.
T. S. Fisher K. E. Torrance 《International journal for numerical methods in engineering》1999,45(11):1631-1655
The complex‐variable boundary element method (CVBEM) is used to analyse forced convection in cooling passages with general, convex cross‐sections. Quadratic spline interpolation is employed for the modelling of coupled velocity and temperature fields. The method is well‐suited for ducts with curved surfaces and high geometric aspect ratios. The method is illustrated for ducts with rounded rectangular, elliptical and rounded diamond cross‐sections. General correlations are presented for the fully developed Nusselt number and Moody friction factor, and the Hagenbach entrance region factor. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
17.
Víctor D. Fachinotti Alberto Cardona Philippe Jetteur 《International journal for numerical methods in engineering》2008,74(6):894-910
This paper introduces a finite element model for the inverse design of pieces with large displacements in the elastic range. The problem consists in determining the initial shape of the piece, such that it attains the designed shape under the effect of service loads. The model is particularly focused on the design of parts with a markedly anisotropic behavior, like laminated turbine blades. Although the formulation expresses equilibrium on the distorted configuration, it uses a standard constitutive equations library that is expressed as usual for measures attached to the undistorted configuration. In this way, the modifications to a standard finite elements code are limited to the routines for the computation of the finite element internal forces and tangent matrix. Two examples are given, the first one for validation purposes, while the second is an application which has industrial interest for the design of turbine blades. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
18.
Nikolaos I. Ioakimidis 《International journal for numerical methods in engineering》2000,47(11):1865-1885
Complex boundary integral equations (Fredholm‐type regular or Cauchy‐type singular or even Hadamard–Mangler‐type hypersingular) have been used for the numerical solution of general plane isotropic elasticity problems. The related Muskhelishvili and, particularly, Lauricella–Sherman equations are famous in the literature, but several more extensions of the Lauricella–Sherman equations have also been proposed. In this paper it is just mentioned that the stress and displacement components can be very accurately computed near either external or internal simple closed boundaries (for anyone of the above equations: regular or singular or hypersingular, but with a prerequisite their actual numerical solution) through the appropriate use of the even more classical elementary Cauchy theorem in complex analysis. This approach has been already used for the accurate numerical computation of analytic functions and their derivatives by Ioakimidis, Papadakis and Perdios (BIT 1991; 31 : 276–285), without applications to elasticity problems, but here the much more complicated case of the elastic complex potentials is studied even when just an appropriate non‐analytic complex density function (such as an edge dislocation/loading distribution density) is numerically available on the boundary. The present results are also directly applicable to inclusion problems, anisotropic elasticity, antiplane elasticity and classical two‐dimensional fluid dynamics, but, unfortunately, not to crack problems in fracture mechanics. Brief numerical results (for the complex potentials), showing the dramatic increase of the computational accuracy, are also displayed and few generalizations proposed. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
19.
K. M. Liew G. P. Zou S. Rajendran 《International journal for numerical methods in engineering》2003,57(5):599-616
In this paper we present a novel spline strip kernel particle method (SSKPM) that has been developed for solving a class of two‐dimensional (2D) elasticity problems. This new approach combines the concepts of the mesh‐free methods and the spline strip method. For the interpolation of the assumed displacement field, we employed the kernel particle shape functions in the transverse direction, and the B3‐spline function in the longitudinal direction. The formulation is validated on several beam and semi‐infinite plate problems. The numerical results of these test problems are then compared with the existing solutions obtained by the exact or numerical methods. From this study we conclude that the SSKPM is a potential alternative to the classical finite strip method (FSM). Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
20.
L. Marin L. Elliott P. J. Heggs D. B. Ingham D. Lesnic X. Wen 《International journal for numerical methods in engineering》2004,60(11):1933-1947
In this paper, several boundary element regularization methods, such as iterative, conjugate gradient, Tikhonov regularization and singular value decomposition methods, for solving the Cauchy problem associated to the Helmholtz equation are developed and compared. Regularizing stopping criteria are developed and the convergence, as well as the stability, of the numerical methods proposed are analysed. The Cauchy problem for the Helmholtz equation can be regularized by various methods, such as the general regularization methods presented in this paper, but more accurate results are obtained by classical methods, such as the singular value decomposition and the Tikhonov regularization methods. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献