Symmetric coupling of multi‐zone curved Galerkin boundary elements with finite elements in elasticity

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.     

Viscoplastic plane‐strain analysis of infinite solids using elastic supports

A highly efficient novel Finite Element Boundary Element Method (FEBEM) is proposed for the elasto‐viscoplastic plane‐strain analysis of displacements and stresses in infinite solids. The proposed method takes advantage of both the Finite Element Method (FEM) and the Boundary Element Method (BEM) to achieve higher efficiency and accuracy by using the concept of elastic supports to simulate the effects of unbounded solid mass surrounding the region of interest. The BEM is used to compute the stiffnesses of elastic supports and to estimate the location of the truncation boundary for the finite element model. As compared to the conventional coupled FEBEM, the proposed method has three main computational advantages. Firstly, the symmetrical and highly banded form of the standard finite element stiffness matrix is not disturbed. Secondly, the proposed technique may be implemented simply by using standard codes for elasto‐viscoplastic finite element analysis and elastic boundary element analysis. Thirdly, the yielded zone is approximately located in advance by using the BEM and hence, an unnecessarily large extent of the domain does not have to be discretized for the finite element modelling. The efficiency and accuracy of the proposed method are demonstrated by computing elastic and elasto‐plastic displacements and stresses around ‘deep’ underground openings in rock mass subject to hydrostatic and non‐hydrostatic in situ stresses. Results obtained by the proposed method are compared with ‘exact’ solutions and with those obtained by using a BEM and a coupled FEBEM. Copyright © 1999 John Wiley & Sons, Ltd.     

The Discontinuity‐Enriched Finite Element Method

We introduce a new methodology for modeling problems with both weak and strong discontinuities independently of the finite element discretization. At variance with the eXtended/Generalized Finite Element Method (X/GFEM), the new method, named the Discontinuity‐Enriched Finite Element Method (DE‐FEM), adds enriched degrees of freedom only to nodes created at the intersection between a discontinuity and edges of elements in the mesh. Although general, the method is demonstrated in the context of fracture mechanics, and its versatility is illustrated with a set of traction‐free and cohesive crack examples. We show that DE‐FEM recovers the same rate of convergence as the standard FEM with matching meshes, and we also compare the new approach to X/GFEM.     

A numerical algorithm for simulation of the electric corona discharge in the triode system

A numerical algorithm is described to calculate the charge density, electric field and corona current distribution in the corona triode. The algorithm employs a hybrid technique based on the Boundary and Finite Element Methods (FEM). FEM is used to determine the electric field because of free space charge produced by the corona discharge. The Boundary Element Method (BEM) is applied for calculating the other component of electric filed as a result of the voltage applied to the electrodes. The Method of Characteristics (MOC) is used to update the space charge density distribution. The characteristic lines are traced backwards from points of the analysed domain to the corona wire. The current density, electric field and space charge density distributions can be controlled by changing the configuration of the system. Results of calculations in a few different cases show the influence of different parameters on the work of the corona triode.     

A combined space–time extended finite element method

The Newmark method for the numerical integration of second order equations has been extensively used and studied along the past fifty years for structural dynamics and various fields of mechanical engineering. Easy implementation and nice properties of this method and its derivatives for linear problems are appreciated but the main drawback is the treatment of discontinuities. Zienkiewicz proposed an approach using finite element concept in time, which allows a new look at the Newmark method. The idea of this paper is to propose, thanks to this approach, the use of a time partition of the unity method denoted Time Extended Finite Element Method (TX‐FEM) for improved numerical simulations of time discontinuities. An enriched basis of shape functions in time is used to capture with a good accuracy the non‐polynomial part of the solution. This formulation allows a suitable form of the time‐stepping formulae to study stability and energy conservation. The case of an enrichment with the Heaviside function is developed and can be seen as an alternative approach to time discontinuous Galerkin method (T‐DGM), stability and accuracy properties of which can be derived from those of the TX‐FEM. Then Space and Time X‐FEM (STX‐FEM) are combined to obtain a unified space–time discretization. This combined STX‐FEM appears to be a suitable technique for space–time discontinuous problems like dynamic crack propagation or other applications involving moving discontinuities. Copyright © 2005 John Wiley & Sons, Ltd.     

A locking‐free selective smoothed finite element method using tetrahedral and triangular elements with adaptive mesh rezoning for large deformation problems

A novel finite element (FE) formulation with adaptive mesh rezoning for large deformation problems is proposed. The proposed method takes the advantage of the selective smoothed FE method (S‐FEM), which has been recently developed as a locking‐free FE formulation with strain smoothing technique. We adopt the selective face‐based smoothed/node‐based smoothed FEM (FS/NS‐FEM‐T4) and edge‐based smoothed/node‐based smoothed FEM (ES/NS‐FEM‐T3) basically but modify them partly so that our method can handle any kind of material constitutive models other than elastic models. We also present an adaptive mesh rezoning method specialized for our S‐FEM formulation with material constitutive models in total form. Because of the modification of the selective S‐FEMs and specialization of adaptive mesh rezoning, our method is locking‐free for severely large deformation problems even with the use of tetrahedral and triangular meshes. The formulation details for static implicit analysis and several examples of analysis of the proposed method are presented in this paper to demonstrate its efficiency. Copyright © 2014 John Wiley & Sons, Ltd.     

The DRM‐MD integral equation method: an efficient approach for the numerical solution of domain dominant problems

This work presents a multi‐domain decomposition integral equation method for the numerical solution of domain dominant problems, for which it is known that the standard Boundary Element Method (BEM) is in disadvantage in comparison with classical domain schemes, such as Finite Difference (FDM) and Finite Element (FEM) methods. As in the recently developed Green Element Method (GEM), in the present approach the original domain is divided into several subdomains. In each of them the corresponding Green's integral representational formula is applied, and on the interfaces of the adjacent subregions the full matching conditions are imposed. In contrast with the GEM, where in each subregion the domain integrals are computed by the use of cell integration, here those integrals are transformed into surface integrals at the contour of each subregion via the Dual Reciprocity Method (DRM), using some of the most efficient radial basis functions known in the literature on mathematical interpolation. In the numerical examples presented in the paper, the contour elements are defined in terms of isoparametric linear elements, for which the analytical integrations of the kernels of the integral representation formula are known. As in the FEM and GEM the obtained global matrix system possesses a banded structure. However in contrast with these two methods (GEM and non‐Hermitian FEM), here one is able to solve the system for the complete internal nodal variables, i.e. the field variables and their derivatives, without any additional interpolation. Finally, some examples showing the accuracy, the efficiency, and the flexibility of the method for the solution of the linear and non‐linear convection–diffusion equation are presented. Copyright © 1999 John Wiley & Sons, Ltd.     

A formulation for electrostatic‐mechanical contact and its numerical solution

The progress in advanced technology fields requires more and more sophisticated formulations to consider contact problems properly. This paper is devoted to the development of a new constitutive model for electrostatic‐mechanical contacts, based on a micro–macro approach to describe the contact behaviour. The electric‐mechanical contact constitutive law is obtained considering the real microscopic shape of the contacting surfaces, the microscopic behaviour of force transmission and current flow. Some thermo‐mechanical macroscopic models based on microscopic characterizations have already been developed to compute the normal and tangential contact stiffness and the thermal contact resistance. On the basis of such macroscopic models, a similar model, suitable for the electric‐mechanical field, is developed. With reference to the thermal constriction resistance the electric contact resistance is studied, assuming a flux tube around each contacting asperity, and choosing a suitable geometry for its narrowing at the contact zone. The contact element geometry is based on well known theoretical and experimental micro‐mechanical laws, suitably adapted for the FEM formulation. The macroscopic stiffness matrix is calculated on the basis of the microscopic laws and it is continuously updated as a function of the changes in the mechanical and electric significant parameters. A consistent linearization of the set of equations is developed to improve the computational speed, within the framework of implicit methods. Copyright © 2005 John Wiley & Sons, Ltd.     

Discrete differential operators on irregular nodes (DDIN)

A previous research made an integral mathematical contribution for obtaining local function interpolation using neighboring nodal values of the solution function. Subsequent researchers developed mesh‐free methods for Finite Element Method (FEM). This principle can also be used to obtain discrete differential operators on irregular nodes. They may be successfully applied to Finite Difference method, Moving Particle Semi‐implicit (MPS) method and Random Collocation Method (RCM). In this paper, we obtain discrete differential operators on irregular nodes and successfully apply them to solve differential equations using the RCM. We also discuss mathematical aspects of the MPS method. Copyright © 2011 John Wiley & Sons, Ltd.     

Stress analysis around crack tips in finite strain problems using the eXtended finite element method

Fracture of rubber‐like materials is still an open problem. Indeed, it deals with modelling issues (crack growth law, bulk behaviour) and computational issues (robust crack growth in 2D and 3D, incompressibility). The present study focuses on the application of the eXtended Finite Element Method (X‐FEM) to large strain fracture mechanics for plane stress problems. Two important issues are investigated: the choice of the formulation used to solve the problem and the determination of suitable enrichment functions. It is demonstrated that the results obtained with the method are in good agreement with previously published works. Copyright © 2005 John Wiley & Sons, Ltd.     

A geometrically non‐linear piezoelectric solid shell element based on a mixed multi‐field variational formulation

This paper is concerned with a geometrically non‐linear solid shell element to analyse piezoelectric structures. The finite element formulation is based on a variational principle of the Hu–Washizu type and includes six independent fields: displacements, electric potential, strains, electric field, mechanical stresses and dielectric displacements. The element has eight nodes with four nodal degrees of freedoms, three displacements and the electric potential. A bilinear distribution through the thickness of the independent electric field is assumed to fulfill the electric charge conservation law in bending dominated situations exactly. The presented finite shell element is able to model arbitrary curved shell structures and incorporates a 3D‐material law. A geometrically non‐linear theory allows large deformations and includes stability problems. Linear and non‐linear numerical examples demonstrate the ability of the proposed model to analyse piezoelectric devices. Copyright © 2005 John Wiley & Sons, Ltd.     

A semi‐Lagrangean time‐integration approach for extended finite element methods

Many computational problems incorporate discontinuities that evolve in time. The eXtendend Finite Element Method (XFEM) is able to represent discontinuities sharply on fixed arbitrary meshes, but numerical difficulties arise if these discontinuities move in time. We point out that this issue is crucial for interface problems with strongly discontinuous fields on fixed grids. A method using semi‐Lagrangean techniques is proposed to adequately handle time integration based on finite difference schemes in the context of the XFEM. The basic idea is to adapt previous numerical solutions to the current interface position by tracking back virtual Lagrangean particles to their previous positions, where an appropriate solution can be extrapolated from a smooth field. Convergence properties of the proposed method in time and space are thoroughly studied for two one‐dimensional model problems. Finally, the method is applied to the particularly challenging problem of premixed combustion, where the discontinuity appears at the flame front separating the burnt from the unburnt gases. A two‐dimensional and a three‐dimensional expanding flame demonstrates that the method is sufficiently accurate to retain the properties of the overall Nitsche‐type formulation for interface problems with embedded strong discontinuities. Copyright © 2014 John Wiley & Sons, Ltd.     

Nodal sensitivities as error estimates in computational mechanics

Summary This paper proposes the use of special sensitivities, called nodal sensitivities, as error indicators and estimators for numerical analysis in mechanics. Nodal sensitivities are defined as rates of change of response quantities with respect to nodal positions. Direct analytical differentiation is used to obtain the sensitivities, and the infinitesimal perturbations of the nodes are forced to lie along the elements. The idea proposed here can be used in conjunction with general purpose computational methods such as the Finite Element Method (FEM), the Boundary Element Method (BEM) or the Finite Difference Method (FDM); however, the BEM is the method of choice in this paper. The performance of the error indicators is evaluated through two numerical examples in linear elasticity.     

Embedded interfaces by polytope FEM

In this paper, the Polytope Finite Element Method is employed to model an embedded interface through the body, independent of the background FEM mesh. The elements that are crossed by the embedded interface are decomposed into new polytope elements which have some nodes on the interface line. The interface introduces discontinuity into the primary variable (strong) or into its derivatives (weak). Both strong and weak discontinuities are studied by the proposed method through different numerical examples including fracture problems with traction‐free and cohesive cracks, and heat conduction problems with Dirichlet and Dirichlet–Neumann types of boundary conditions on the embedded interface. For traction‐free cracks which have tip singularity, the nodes near the crack tip are enriched with the singular functions through the eXtended Finite Element Method. The concept of Natural Element Coordinates (NECs) is invoked to drive shape functions for the produced polytopes. A simple treatment is proposed for concave polytopes produced by a kinked interface and also for locating crack tip inside an element prior to using the singularity enrichment. The proposed method pursues some implementational details of eXtended/Generalized Finite Element Methods for interfaces. But here the additional DOFs are constructed on the interface lines in contrast to X/G‐FEM, which attach enriched DOFs to the previously existed nodes. Copyright © 2011 John Wiley & Sons, Ltd.     

Overcoming element quality dependence of finite elements with adaptive extended stencil FEM (AES‐FEM)

The finite element methods (FEMs) are important techniques in engineering for solving partial differential equations, but they depend heavily on element shape quality for stability and good performance. In this paper, we introduce the Adaptive Extended Stencil Finite Element Method (AES‐FEM) as a means for overcoming this dependence on element shape quality. Our method replaces the traditional basis functions with a set of generalized Lagrange polynomial basis functions, which we construct using local weighted least‐squares approximations. The method preserves the theoretical framework of FEM and allows imposing essential boundary conditions and integrating the stiffness matrix in the same way as the classical FEM. In addition, AES‐FEM can use higher‐degree polynomial basis functions than the classical FEM, while virtually preserving the sparsity pattern of the stiffness matrix. We describe the formulation and implementation of AES‐FEM and analyze its consistency and stability. We present numerical experiments in both 2D and 3D for the Poisson equation and a time‐independent convection–diffusion equation. The numerical results demonstrate that AES‐FEM is more accurate than linear FEM, is also more efficient than linear FEM in terms of error versus runtime, and enables much better stability and faster convergence of iterative solvers than linear FEM over poor‐quality meshes. Copyright © 2016 John Wiley & Sons, Ltd.     

Issues in convergence improvement for non-linear finite element programs

Systems of non-linear equations as they arise when analysing various physical phenomena and technological processes by the implicit Finite Element Method (FEM) are commonly solved by the Newton–Raphson method. The modelling of sheet metal forming processes is one example of highly non-linear problems where the iterative solution procedure can become very slow or diverge. This paper focuses on techniques to overcome these numerical difficulties. Several methods to generate initial guesses within the radius of convergence are proposed. Appropriate stopping criteria for the iterative procedure are discussed. A combination of various line search methods with the continuation method is proposed. The efficiency and robustness of these numerical procedures are compared based on a set of test examples. A particular form of line search was identified which allows the stable and efficient solution of highly non-linear sheet metal forming problems. Even though the present investigations were motivated by the application of the implicit FEM to the simulation of sheet metal forming processes, the findings are general enough to be applicable to a wide spectrum of non-linear FEM applications. © 1997 John Wiley & Sons, Ltd.     

楔形工字梁抗剪极限承载力试验研究

该文介绍了楔形变截面工字钢短梁的抗剪承载力试验,与有限元结果进行了对比,提出了将楔形变截面梁的承载力按等截面梁的承载力乘以楔率折减系数的计算方法,与试验和扩充的有限元计算结果对比表明:公式有很好的精度.     

Four-parameter functionally graded cracked plates of arbitrary shape: A GDQFEM solution for free vibrations

The dynamic behavior of moderately thick FGM plates with geometric discontinuities and arbitrarily curved boundaries is investigated. The Generalized Differential Quadrature Finite Element Method (GDQFEM) is proposed as a numerical approach. The irregular physical domain in Cartesian coordinates is transformed into a regular domain in natural coordinates. Several types of cracked FGM plates are investigated. It appears that GDQFEM is analogous to the well-known Finite Element Method (FEM). With reference to the proposed technique the governing FSDT equations are solved in their strong form and the connections between the elements are imposed with the inter-element compatibility conditions. The results show excellent agreement with other numerical solutions obtained by FEM.