共查询到20条相似文献,搜索用时 15 毫秒
1.
S. Ganguly J. B. Layton C. Balakrishna 《International journal for numerical methods in engineering》2000,48(5):633-654
The coupling of Finite Element Method (FEM) with a Boundary Element Method (BEM) is a desirable result that exploits the advantages of each. This paper examines the efficient symmetric coupling of a Symmetric Galerkin Multi‐zone Curved Boundary Element Analysis method with a Finite Element Method for 2‐D elastic problems. Existing collocation based multi‐zone boundary element methods are not symmetric. Thus, when they are coupled with FEM, it is very difficult to achieve symmetry, increasing the computational work to solve the problem. This paper uses a fully Symmetric curved Multi‐zone Galerkin Boundary Element Approach that is coupled to an FEM in a completely symmetric fashion. The symmetry is achieved by symmetrically converting the boundary zones into equivalent ‘macro finite elements’, that are symmetric, so that symmetry in the coupling is retained. This computationally efficient and fast approach can be used to solve a wide range of problems, although only 2‐D elastic problems are shown. Three elasticity problems, including one from the FEM‐BEM literature that explore the efficacy of the approach are presented. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
2.
S. Ganguly J. B. Layton C. Balakrishna 《International journal for numerical methods in engineering》2004,59(8):1021-1038
When the different parts of a structure are modelled independently by BEM or FEM methods, it is sometimes necessary to put the parts together without remeshing of the nodes along the part interfaces. Frequently the nodes do not match along the interface. In this work, the symmetric Galerkin multi‐zone curved boundary element is a fully symmetric formulation and is the method used for the boundary element part. For BEM–FEM coupling it is then necessary to interpolate the tractions in‐between the non‐matching nodes for the FEM part. Finally, the coupling is achieved by transforming the finite element domains to equivalent boundary element domains in a block symmetric formulation. This system is then coupled with a boundary element domain with non‐matching nodes in‐between. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
3.
J. B. Layton S. Ganguly C. Balakrishna J. H. Kane 《International journal for numerical methods in engineering》1997,40(16):2913-2931
The recent development of the symmetric Galerkin approach to boundary element analysis (BEA) has been demonstrated to be superior to the collocation method for medium to large problems. This fact has been shown in both heat conduction and elasticity. Accounts of collocation multi-zone analysis techniques have also been prevalent in the literature, dealing with multiple boundary integral relations associated with portions of overall objects. This technique results in overall system matrices with a blocked, sparse, but unsymmetric character. It has been shown that multi-zone techniques can produce smaller solution times than a single zone analysis for large problems. These techniques are useful for multi-material problems as well. This paper presents an approach for combining the benefits of both techniques resulting in a Galerkin multi-zone method, that is overall unsymmetric but contains a significant amount of block symmetry. A condensation technique in the multi-zone solver is shown to exploit the symmetry of the Galerkin formulation as well as the blocked sparsity of the multi-zone technique. This method is compared to collocation multi-zone on two elasticity problems from the literature. It is concluded that an appropriate implementation of the symmetric Galerkin multi-zone BEA indeed has the potential to be superior to the collocation based multi-zone BEA, for medium to large-scale elasticity problems. © 1997 John Wiley & Sons, Ltd. 相似文献
4.
J. J. Perez-Gavilan M. H. Aliabadi 《International journal for numerical methods in engineering》2003,57(12):1661-1693
A symmetric Galerkin boundary element formulation is developed for shear deformable plates. A mixed strategy is used for the integration process, i.e. partial regularization using simple solutions followed by a singularity subtraction technique. For the shear equation, full regularization is achieved using new kernel relationships found through a constant shear mode of deformation. Some of the strong singular integrals are avoided altogether by using a modified traction obtained through a very simple variable change; appropriate boundary conditions are defined. Details of the implementation are given and several example problems solved to verify the accuracy of the proposed formulation. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
5.
Attilio Frangi Marc Bonnet 《International journal for numerical methods in engineering》1998,41(2):337-369
A variational Boundary Element formulation is proposed for the solution of the elastic Kirchhoff plate bending problem. The stationarity conditions of an augmented potential energy functional are first discussed. After addressing the topic of the choice of the test functions, a regularization process based on integrations by parts is developed, which allows to express the formulation in terms of double integrals, the inner being at most weakly singular and the outer regular. Standard integration procedures may then be applied for their numerical evaluation in the presence of both straight and curved boundaries. The normal slope and the vertical displacement must be C0 and C1 continuous, respectively. Numerical examples show, through comparisons with analytical solutions, that a high accuracy is achieved. © 1998 John Wiley & Sons, Ltd. 相似文献
6.
R. Vodička V. Mantič F. París 《International journal for numerical methods in engineering》2010,83(1):91-128
An original approach to the solution of linear elastic domain decomposition problems by the symmetric Galerkin boundary element method is developed. The approach is based on searching for the saddle‐point of a new potential energy functional with Lagrange multipliers. The interfaces can be either straight or curved, open or closed. The two coupling conditions, equilibrium and compatibility, along an interface are fulfilled in a weak sense by means of Lagrange multipliers (interface displacements and tractions), which enables non‐matching meshes to be used at both sides of interfaces between subdomains. The accuracy and robustness of the method is tested by several numerical examples, where the numerical results are compared with the analytical solution of the solved problems, and the convergence rates of two error norms are evaluated for h‐refinements of matching and non‐matching boundary element meshes. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
7.
J. J. Prez‐Gaviln M. H. Aliabadi 《International journal for numerical methods in engineering》2004,61(7):1093-1106
A Galerkin boundary element formulation for shear deformable plate bending dynamics is developed. The formulation makes use of the static fundamental solutions for the weighted residual integral equations. The domain integrals carrying the inertia terms and generic static loads are considered as body forces and approximated with boundary values using the dual reciprocity method. The load is modelled as a series of impact loads of time varying intensity and moving in space in a predetermined path. The formulation was implemented and tested solving a benchmark problem. The results are compared with finite element solutions. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
8.
Delfim Soares Jr 《International journal for numerical methods in engineering》2009,78(9):1076-1093
The numerical modelling of interacting acoustic media by boundary element method–finite element method (BEM–FEM) coupling procedures is discussed here, taking into account time‐domain approaches. In this study, the global model is divided into different sub‐domains and each sub‐domain is analysed independently (considering BEM or FEM discretizations): the interaction between the different sub‐domains of the global model is accomplished by interface procedures. Numerical formulations based on FEM explicit and implicit time‐marching schemes are discussed, resulting in direct and optimized iterative BEM–FEM coupling techniques. A multi‐level time‐step algorithm is considered in order to improve the flexibility, accuracy and stability (especially when conditionally stable time‐marching procedures are employed) of the coupled analysis. At the end of the paper, numerical examples are presented, illustrating the potentialities and robustness of the proposed methodologies. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
9.
A.‐V. Phan J. A. L. Napier L. J. Gray T. Kaplan 《International journal for numerical methods in engineering》2003,57(6):835-851
A symmetric‐Galerkin boundary element framework for fracture analysis with frictional contact (crack friction) on the crack surfaces is presented. The algorithm employs a continuous interpolation on the crack surface (utilizing quadratic boundary elements) and enables the determination of two important quantities for the problem, namely the local normal tractions and sliding displacements on the crack surfaces. An effective iterative scheme for solving this non‐linear boundary value problem is proposed. The results of test examples are compared with available analytical solutions or with those obtained from the displacement discontinuity method (DDM) using linear elements and internal collocation. The results demonstrate that the method works well for difficult kinked/junction crack problems. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
10.
G. Bonnet 《International journal for numerical methods in engineering》2009,80(8):1110-1123
The symmetric Galerkin boundary element method is used to solve boundary value problems by keeping the symmetric nature of the matrix obtained after discretization. The matrix elements are obtained from a double integral involving the double derivative of Green's operator, which is highly singular. The paper presents a regularization of the hypersingular integrals which depend only on the properties of Green's tensor. The method is presented in the case of Laplace's operator, with an example of application. The case of elasticity is finally addressed theoretically, showing an easy extension to any case of anisotropy. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
11.
Jean‐François Remacle Nicolas Chevaugeon Émilie Marchandise Christophe Geuzaine 《International journal for numerical methods in engineering》2007,69(4):750-771
A general method for the post‐processing treatment of high‐order finite element fields is presented. The method applies to general polynomial fields, including discontinuous finite element fields. The technique uses error estimation and h‐refinement to provide an optimal visualization grid. Some filtering is added to the algorithm in order to focus the refinement on a visualization plane or on the computation of one single iso‐zero surface. 2D and 3D examples are provided that illustrate the power of the technique. In addition, schemes and algorithms that are discussed in the paper are readily available as part of an open source project that is developed by the authors, namely Gmsh. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
12.
Christian Cordes Mario Putti 《International journal for numerical methods in engineering》2001,52(4):371-387
In standard finite element simulations of groundwater flow the correspondence between hydraulic head gradients and groundwater fluxes is represented by the stiffness matrix. In two‐dimensional problems the use of linear triangular elements on Delaunay triangulations guarantees a stiffness matrix of type M. This implies that the local numerical fluxes are physically consistent with Darcy's law. This condition is fundamental to avoid the occurrence of local maxima or minima, and is of crucial importance when the calculated flow field is used in contaminant transport simulations or pathline evaluation. In three spatial dimensions, the linear Galerkin approach on tetrahedra does not lead to M‐matrices even on Delaunay meshes. By interpretation of the Galerkin approach as a subdomain collocation scheme, we develop a new approach (OSC, orthogonal subdomain collocation) that is shown to produce M‐matrices in three‐dimensional Delaunay triangulations. In case of heterogeneous and anisotropic coefficients, extra mesh properties required for M‐stiffness matrices will also be discussed. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
13.
Marc Bonnet 《International journal for numerical methods in engineering》2003,57(8):1053-1083
Procedures based on group representation theory, allowing the exploitation of geometrical symmetry in symmetric Galerkin BEM formulations, are investigated. In particular, this investigation is based on the weaker assumption of partial geometrical symmetry, where the boundary has two disconnected components, one of which is symmetric; e.g. this can be very useful for defect identification problems. The main development is expounded in the context of 3D Neumann elastostatic problems, considered as model problems; and then extended to SGBIE formulations for Dirichlet and/or scalar problems. Both Abelian and non‐Abelian finite symmetry groups are considered. The effectiveness of the present approach is demonstrated through numerical examples, where both partial and complete symmetry are considered, in connection with both Abelian and non‐Abelian symmetry groups. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
14.
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. 相似文献
15.
K. Davey M. T. Alonso Rasgado I. Rosindale 《International journal for numerical methods in engineering》1999,44(8):1031-1054
A semi‐analytical integration scheme is described in this paper which is designed to reduce the errors incurred when integrals with singular integrands are evaluated numerically. This new scheme can be applied to linear triangular elements for use in steady‐state elastodynamic BEM problems and is particularly useful for predicting displacement to high accuracy, close to surfaces for a spectrum of frequencies. The scheme involves the application of Taylor expansions to formulate the integrals into two parts. One part is regular and is evaluated numerically and the other part is singular but sufficiently simple to enable its transforma tion into a line integral. The line integral is solved numerically using Gauss–Legendre quadrature. This approach caters for all the integral types that appear in steady‐state elastodynamic boundary elements but, in particular, no special treatment is required for the evaluation of the Cauchy principal value singular integrals. Numerical tests are performed on a simple test‐problem for which a known analytical solution exists. The results obtained using the semi‐analytical approach are shown to be considerably more accurate than those obtained using standard quadrature methods. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
16.
Mark Gates Karel Matouš Michael T. Heath 《International journal for numerical methods in engineering》2008,76(9):1353-1378
We develop an asynchronous time integration and coupling method with domain decomposition for linear and non‐linear problems in mechanics. To ensure stability in the time integration and in coupling between domains, we use variational integrators with local Lagrange multipliers to enforce continuity at the domain interfaces. The asynchronous integrator lets each domain step with its own time step, using a smaller time step where required by stability and accuracy constraints and a larger time step where allowed. We show that in practice the time step is limited by accuracy requirements rather than by stability requirements. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
17.
Alok Sutradhar Glaucio H. Paulino L. J. Gray 《International journal for numerical methods in engineering》2005,62(1):122-157
A symmetric Galerkin formulation and implementation for heat conduction in a three‐dimensional functionally graded material is presented. The Green's function of the graded problem, in which the thermal conductivity varies exponentially in one co‐ordinate, is used to develop a boundary‐only formulation without any domain discretization. The main task is the evaluation of hypersingular and singular integrals, which is carried out using a direct ‘limit to the boundary’ approach. However, due to complexity of the Green's function for graded materials, the usual direct limit procedures have to be modified, incorporating Taylor expansions to obtain expressions that can be integrated analytically. Several test examples are provided to verify the numerical implementation. The results of test calculations are in good agreement with exact solutions and corresponding finite element method simulations. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
18.
L. Bergamaschi G. Pini F. Sartoretto 《International journal for numerical methods in engineering》2005,63(15):2069-2085
Finite element discretizations of flow problems involving multiaquifer systems deliver large, sparse, unstructured matrices, whose partial eigenanalysis is important for both solving the flow problem and analysing its main characteristics. We studied and implemented an effective preconditioning of the Jacobi–Davidson algorithm by FSAI‐type preconditioners. We developed efficient parallelization strategies in order to solve very large problems, which could not fit into the storage available to a single processor. We report our results about the solution of multiaquifer flow problems on an SP4 machine and a Linux Cluster. We analyse the sequential and parallel efficiency of our algorithm, also compared with standard packages. Questions regarding the parallel solution of finite element eigenproblems are addressed and discussed. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
19.
This paper is concerned with the minimization of functionals of the form ∫Γ(b) f( h ,T( b, h )) dΓ( b ) where variation of the vector b modifies the shape of the domain Ω on which the potential problem, ?2T=0, is defined. The vector h is dependent on non‐linear boundary conditions that are defined on the boundary Γ. The method proposed is founded on the material derivative adjoint variable method traditionally used for shape optimization. Attention is restricted to problems where the shape of Γ is described by a boundary element mesh, where nodal co‐ordinates are used in the definition of b . Propositions are presented to show how design sensitivities for the modified functional ∫Γ(b) f( h ,T ( b, h )) dΓ( b ) +∫Ω(b) λ( b, h ) ?2T( b, h ) dΩ( b ) can be derived more readily with knowledge of the form of the adjoint function λ determined via non‐shape variations. The methods developed in the paper are applied to a problem in pressure die casting, where the objective is the determination of cooling channel shapes for optimum cooling. The results of the method are shown to be highly convergent. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
20.
Horst Lanzerath Heinz Waller 《International journal for numerical methods in engineering》1999,45(7):841-864
The boundary element method in combination with modal techniques is used to calculate the response of transient excited structures in the time domain numerically. If the system matrices of a structure are evaluated with a fundamental solution in the frequency domain these matrices become functions of frequency which normally cannot be expressed analytically. The associated eigenvalue problem therefore is non‐linear and difficult to solve. For simplification a series expansion formula for the fundamental solution is used in different frequency ranges. Then the eigenvalue problem can be linearized and solved by direct or iterative methods. By using the orthogonal properties of the eigenfunctions, the normal modes of the dynamic problem can be uncoupled as is well known in vibration analysis. That way the transient response of a dynamic excited system in the time domain can be determined without difficulties. Displacements and stresses at different points of the structure are the result. Difficulties in the formulation of time‐dependent problems using the boundary element method can be avoided. There is no problem in considering modal damping factors, for general damping characteristics the associated fundamental solutions have to be found. Several examples are studied in the paper to illustrate how the new method can be applied. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献