Trefftz boundary elements—multi‐region formulations

This paper is concerned with an effective numerical implementation of the Trefftz boundary element method, for the analysis of two‐dimensional potential problems, defined in arbitrarily shaped domains. The domain is first discretized into multiple subdomains or regions. Each region is treated as a single domain, either finite or infinite, for which a complete set of solutions of the problem is known in the form of an expansion with unknown coefficients. Through the use of weighted residuals, this solution expansion is then forced to satisfy the boundary conditions of the actual domain of the problem, leading thus to a system of equations, from which the unknowns can be readily determined. When this basic procedure is adopted, in the analysis of multiple‐region problems, proper boundary integral equations must be used, along common region interfaces, in order to couple to each other the unknowns of the solution expansions relative to the neighbouring regions. These boundary integrals are obtained from weighted residuals of the coupling conditions which allow the implementation of any order of continuity of the potential field, across the interface boundary, between neighbouring regions. The technique used in the formulation of the region‐coupling conditions drives the performance of the Trefftz boundary element method. While both of the collocation and Galerkin techniques do not generate new unknowns in the problem, the technique of Galerkin presents an additional and unique feature: the size of the matrix of the final algebraic system of equations which is always square and symmetric, does not depend on the number of boundary elements used in the discretization of both the actual and region‐interface boundaries. This feature which is not shared by other numerical methods, allows the Galerkin technique of the Trefftz boundary element method to be effectively applied to problems with multiple regions, as a simple, economic and accurate solution technique. A very difficult example is analysed with this procedure. The accuracy and efficiency of the implementations described herein make the Trefftz boundary element method ideal for the study of potential problems in general arbitrarily‐shaped two‐dimensional domains. Copyright © 1999 John Wiley & Sons, Ltd.     

Solution of non-linear boundary value problems by discrete least squares

This paper presents a new discrete (point-matching) least squares method for solving general, non-linear boundary value problems. The technique uses a mixed formulation and minimizes the sum of squared residuals by Powell's (1965) algorithm. The present method is compared with finite difference and weighted residual finite element methods in the examples and discussion. It shows considerable promise as a general differential equation solver for both large and small computer systems.     

The scaled boundary finite element method—alias consistent infinitesimal finite element cell method—for diffusion

The scaled boundary finite element method, alias the consistent infinitesimal finite element cell method, is developed starting from the diffusion equation. Only the boundary of the medium is discretized with surface finite elements yielding a reduction of the spatial dimension by one. No fundamental solution is necessary, and thus no singular integrals need to be evaluated. Essential and natural boundary conditions on surfaces and conditions on interfaces between different materials are enforced exactly without any discretization. The solution of the function in the radial direction is analytical. This method is thus exact in the radial direction and converges to the exact solution in the finite element sense in the circumferential directions. The semi‐analytical solution inside the domain leads to an efficient procedure to calculate singularities accurately without discretization in the vicinity of the singular point. For a bounded medium symmetric steady‐state stiffness and mass matrices with respect to the degrees of freedom on the boundary result without any additional assumption. Copyright © 1999 John Wiley & Sons, Ltd.     

The simple boundary element method for multiple cracks in functionally graded media governed by potential theory: a three‐dimensional Galerkin approach

The simple boundary element method consists of recycling existing codes for homogeneous media to solve problems in non‐homogeneous media while maintaining a purely boundary‐only formulation. Within this scope, this paper presents a ‘simple’ Galerkin boundary element method for multiple cracks in problems governed by potential theory in functionally graded media. Steady‐state heat conduction is investigated for thermal conductivity varying either parabolically, exponentially, or trigonometrically in one or more co‐ordinates. A three‐dimensional implementation which merges the dual boundary integral equation technique with the Galerkin approach is presented. Special emphasis is given to the treatment of crack surfaces and boundary conditions. The test examples simulated with the present method are verified with finite element results using graded finite elements. The numerical examples demonstrate the accuracy and efficiency of the present method especially when multiple interacting cracks are involved. Copyright © 2005 John Wiley & Sons, Ltd.     

Assumed‐deformation gradient finite elements with nodal integration for nearly incompressible large deformation analysis

An assumed‐strain finite element technique for non‐linear finite deformation is presented. The weighted‐residual method enforces weakly the balance equation with the natural boundary condition and also the kinematic equation that links the elementwise and the assumed‐deformation gradient. Assumed gradient operators are derived via nodal integration from the kinematic‐weighted residual. A variety of finite element shapes fits the derived framework: four‐node tetrahedra, eight‐, 27‐, and 64‐node hexahedra are presented here. Since the assumed‐deformation gradients are expressed entirely in terms of the nodal displacements, the degrees of freedom are only the primitive variables (displacements at the nodes). The formulation allows for general anisotropic materials and no volumetric/deviatoric split is required. The consistent tangent operator is inexpensive and symmetric. Furthermore, the material update and the tangent moduli computation are carried out exactly as for classical displacement‐based models; the only deviation is the consistent use of the assumed‐deformation gradient in place of the displacement‐derived deformation gradient. Examples illustrate the performance with respect to the ability of the present technique to resist volumetric locking. A constraint count can partially explain the insensitivity of the resulting finite element models to locking in the incompressible limit. Copyright © 2008 John Wiley & Sons, Ltd.     

Mixed plate bending elements based on least‐squares formulation

A finite element formulation for the bending of thin and thick plates based on least‐squares variational principles is presented. Finite element models for both the classical plate theory and the first‐order shear deformation plate theory (also known as the Kirchhoff and Mindlin plate theories, respectively) are considered. High‐order nodal expansions are used to construct the discrete finite element model based on the least‐squares formulation. Exponentially fast decay of the least‐squares functional, which is constructed using the L₂ norms of the equations residuals, is verified for increasing order of the nodal expansions. Numerical examples for the bending of circular, rectangular and skew plates with various boundary conditions and plate thickness are presented to demonstrate the predictive capability and robustness of the new plate bending elements. Plate bending elements based on this formulation are shown to be insensitive to both shear‐locking and geometric distortions. Copyright © 2004 John Wiley & Sons, Ltd.     

DQEM for free vibration analysis of Timoshenko beams on elastic foundations

 A differential quadrature element method (DQEM) based on first order shear deformation theory is developed for free vibration analysis of non-uniform beams on elastic foundations. By decomposing the system into a series of sub-domains or elements, any discontinuity in loading, geometry, material properties, and even elastic foundations can be considered conveniently. Using this method, the vibration analysis of general beam-like structures is to be studied. The governing equations of each element, natural compatibility conditions at the interface of two adjacent elements and the external boundary conditions are developed in a systematic manner, using Hamilton's principle. The present DQEM is to be implemented to Timoshenko beams resting on partially supported elastic foundations with various types of boundary conditions under the action of axial loading. The general versality, accuracy, and efficiency of the presented DQEM are demonstrated having solved different examples and compared to the exact or other numerical procedure solutions. Received: 11 October 2002/Accepted: 26 November 2002     

Finite elements for materials with strain gradient effects

A finite element implementation is reported of the Fleck–Hutchinson phenomenological strain gradient theory. This theory fits within the Toupin–Mindlin framework and deals with first‐order strain gradients and the associated work‐conjugate higher‐order stresses. In conventional displacement‐based approaches, the interpolation of displacement requires C¹‐continuity in order to ensure convergence of the finite element procedure for higher‐order theories. Mixed‐type finite elements are developed herein for the Fleck–Hutchinson theory; these elements use standard C⁰‐continuous shape functions and can achieve the same convergence as C¹ elements. These C⁰ elements use displacements and displacement gradients as nodal degrees of freedom. Kinematic constraints between displacement gradients are enforced via the Lagrange multiplier method. The elements developed all pass a patch test. The resulting finite element scheme is used to solve some representative linear elastic boundary value problems and the comparative accuracy of various types of element is evaluated. Copyright © 1999 John Wiley & Sons, Ltd.     

Isoparametric Hermite elements

Isoparametric Hermite elements are created using Bogner–Fox–Schmit rectangles on a reference domain and mapping these numerically onto the computational domain. The difficulties involved in devising explicit C¹ shape functions for isoparametric elements are thus avoided, and the resulting elements have all the benefits of full C¹ continuity, the simplicity of the Bogner–Fox–Schmit element and the geometrical flexibility expected from higher-order isoparametric elements. The numerical mapping consists in the finite element solution of a linear boundary value problem, which is inexpensive and is carried out as a preprocessing operation—the required derivatives of the mapping then being supplied to the main analysis as data. Some care is required in defining the differential boundary conditions, and guidance on this is provided. Examples are given showing the success of the mapping procedure, and the use of the resulting elements in the solution of some boundary value problems. The numerical results confirm a convergence analysis provided for the new isoparametric Hermite element.     

A super‐element for crack analysis in the time domain

A super‐element for the dynamic analysis of two‐dimensional crack problems is developed based on the scaled boundary finite‐element method. The boundary of the super‐element containing a crack tip is discretized with line elements. The governing partial differential equations formulated in the scaled boundary co‐ordinates are transformed to ordinary differential equations in the frequency domain by applying the Galerkin's weighted residual technique. The displacements in the radial direction from the crack tip to a point on the boundary are solved analytically without any a priori assumption. The scaled boundary finite‐element formulation leads to symmetric static stiffness and mass matrices. The super‐element can be coupled seamlessly with standard finite elements. The transient response is evaluated directly in the time domain using a standard time‐integration scheme. The stress field, including the singularity around the crack tip, is expressed semi‐analytically. The stress intensity factors are evaluated without directly addressing singular functions, as the limit in their definitions is performed analytically. Copyright © 2004 John Wiley & Sons, Ltd.     

Efficient non‐linear solid–fluid interaction analysis by an iterative BEM/FEM coupling

An iterative coupling of finite element and boundary element methods for the time domain modelling of coupled fluid–solid systems is presented. While finite elements are used to model the solid, the adjacent fluid is represented by boundary elements. In order to perform the coupling of the two numerical methods, a successive renewal of the variables on the interface between the two subdomains is performed through an iterative procedure until the final convergence is achieved. In the case of local non‐linearities within the finite element subdomain, it is straightforward to perform the iterative coupling together with the iterations needed to solve the non‐linear system. In particular a more efficient and a more stable performance of the new coupling procedure is achieved by a special formulation that allows to use different time steps in each subdomain. Copyright © 2005 John Wiley & Sons, Ltd.     

The boundary element method applied to the analysis of two-dimensional nonlinear sloshing problems

This paper deals with an application of the boundary element method to the analysis of nonlinear sloshing problems, namely nonlinear oscillations of a liquid in a container subjected to forced oscillations. First, the problem is formulated mathematically as a nonlinear initial-boundary value problem by the use of a governing differential equation and boundary conditions, assuming the fluid to be inviscid and incompressible and the flow to be irrotational. Next, the governing equation (Laplace equation) and boundary conditions, except the dynamic boundary condition on the free surface, are transformed into an integral equation by employing the Galerkin method. Two dynamic boundary condition is reduced to a weighted residual equation by employing the Galerkin method. Two equations thus obtained are discretized by the use of the finite element method spacewise and the finite difference method timewise. Collocation method is employed for the discretization of the integral equation. Due to the nonlinearity of the problem, the incremental method is used for the numerical analysis. Numerical results obtained by the present boundary element method are compared with those obtained by the conventional finite element method and also with existing analytical solutions of the nonlinear theory. Good agreements are obtained, and this indicates the availability of the boundary element method as a numerical technique for nonlinear free surface fluid problems.     

Assumed strain finite element methods for conserving temporal integrations in non‐linear solid dynamics

This paper presents a new class of assumed strain finite elements to use in combination with general energy‐momentum‐conserving time‐stepping algorithms so that these conservation properties in time are preserved by the fully discretized system in space and time. The case of interest corresponds to nearly incompressible material responses, in the fully non‐linear finite strain elastic and elastoplastic ranges. The new elements consider the classical scaling of the deformation gradient with an assumed Jacobian (its determinant) defined locally through a weighted averaging procedure at the element level. The key aspect of the newly proposed formulation is the definition of the associated linearized strain operator or B‐bar operator. The developments presented here start by identifying the conditions that this discrete operator must satisfy for the fully discrete system in time and space to inherit exactly the conservation laws of linear and angular momenta, and the conservation/dissipation law of energy for elastic and inelastic problems, respectively. Care is also taken of the preservation of the relative equilibria and the corresponding group motions associated with the momentum conservation laws, and characterized by purely rotational and translational motions superimposed to the equilibrium deformed configuration. With these developments at hand, a new general B‐bar operator is introduced that satisfies these conditions. The new operator not only accounts for the spatial interpolations (e.g. bilinear displacements with piece‐wise constant volume) but also depends on the discrete structure of the equations in time. The aforementioned conservation/dissipation properties of energy and momenta are then proven to hold rigorously for the final numerical schemes, unconditionally of the time step size and the material model (elastic or elastoplastic). Different finite elements are considered in this framework, including quadrilateral and triangular elements for plane problems and brick elements for three‐dimensional problems. Several representative numerical simulations are presented involving, in particular, the use of energy‐dissipating momentum‐conserving time‐stepping schemes recently developed by the author and co‐workers for general finite strain elastoplasticity in order to illustrate the properties of the new finite elements, including these conservation/dissipation properties in time and their locking‐free response in the quasi‐incompressible case. Copyright © 2007 John Wiley & Sons, Ltd.     

Finite element response sensitivity analysis using force‐based frame models

This paper presents a method to compute consistent response sensitivities of force‐based finite element models of structural frame systems to both material constitutive and discrete loading parameters. It has been shown that force‐based frame elements are superior to classical displacement‐based elements in the sense that they enable, at no significant additional costs, a drastic reduction in the number of elements required for a given level of accuracy in the computed response of the finite element model. This advantage of force‐based elements is of even more interest in structural reliability analysis, which requires accurate and efficient computation of structural response and structural response sensitivities. This paper focuses on material non‐linearities in the context of both static and dynamic response analysis. The formulation presented herein assumes the use of a general‐purpose non‐linear finite element analysis program based on the direct stiffness method. It is based on the general so‐called direct differentiation method (DDM) for computing response sensitivities. The complete analytical formulation is presented at the element level and details are provided about its implementation in a general‐purpose finite element analysis program. The new formulation and its implementation are validated through some application examples, in which analytical response sensitivities are compared with their counterparts obtained using forward finite difference (FFD) analysis. The force‐based finite element methodology augmented with the developed procedure for analytical response sensitivity computation offers a powerful general tool for structural response sensitivity analysis. Copyright © 2004 John Wiley & Sons, Ltd.     

An h‐hierarchical adaptive procedure for the scaled boundary finite‐element method

The scaled boundary finite‐element method (a novel semi‐analytical method for solving linear partial differential equations) involves the solution of a quadratic eigenproblem, the computational expense of which rises rapidly as the number of degrees of freedom increases. Consequently, it is desirable to use the minimum number of degrees of freedom necessary to achieve the accuracy desired. Stress recovery and error estimation techniques for the method have recently been developed. This paper describes an h‐hierarchical adaptive procedure for the scaled boundary finite‐element method. To allow full advantage to be taken of the ability of the scaled boundary finite‐element method to model stress singularities at the scaling centre, and to avoid discretization of certain adjacent segments of the boundary, a sub‐structuring technique is used. The effectiveness of the procedure is demonstrated through a set of examples. The procedure is compared with a similar h‐hierarchical finite element procedure. Since the error estimators in both cases evaluate the energy norm of the stress error, the computational cost of solutions of similar overall accuracy can be compared directly. The examples include the first reported direct comparison of the computational efficiency of the scaled boundary finite‐element method and the finite element method. The scaled boundary finite‐element method is found to reduce the computational effort considerably. Copyright © 2002 John Wiley & Sons, Ltd.     

Consistent coupling of beam and shell models for thermo‐elastic analysis

In this paper, the finite element formulation of a transition element for consistent coupling between shell and beam finite element models of thin‐walled beam‐like structures in thermo‐elastic problems is presented. Thin‐walled beam‐like structures modelled only with beam elements cannot be used to study local stress concentrations or to provide local mechanical or thermal boundary conditions. For this purpose, the structure has to be modelled using shell elements. However, computations using shell elements are a lot more expensive as compared to beam elements. The finite element model can be more efficient when the shell elements are used only in regions where the local effects are to be studied or local boundary conditions have to be provided. The remaining part of the structure can be modelled with beam elements. To couple these two models (i.e. shell and beam models) at transitional cross‐sections, transition elements are derived here for thermo‐elastic problems. The formulation encloses large displacement and rotational behaviour, which is important in case of thin‐walled beam‐like structures. Copyright © 2004 John Wiley & Sons, Ltd.     

Stress analysis around holes or notches by special finite elements

In this paper, three special elements for stress concentration analysis around holes or notches are presented. One is based on a variational principle involving displacements, unsymmetric stresses and rotations as independent variables, the other two on the classical Hellinger–Reissner principle. The new elements are designed to be used in a single layer along the cavity rim, where enhanced properties are needed to get a satisfactory level of accuracy. Special properties are obtained by enforcing the stress approximation to satisfy a priori the traction‐free boundary condition on the element side positioned along the cavity rim and by properly enriching the displacement interpolation. The proposed elements are: applicable to geometries of arbitrary shape, easy‐to‐use, low‐cost and readily implementable into existing finite‐element codes. Copyright © 2005 John Wiley & Sons, Ltd.     

A coupled boundary/finite element method for the computation of magnetically and electrostatically levitated droplet shapes

A coupled finite element and boundary element method is developed to predict the magnetic vector and scalar potential distributions in the droplets levitated in an alternating magnetic or electrostatic field. The computational algorithm entails the application of boundary elements in the region of free space and finite elements in the droplet region, the two being coupled along the droplet–air interface. The coupled boundary and finite element scheme is further integrated with a WRM‐based algorithm to predict the free surface deformation of magnetically and electrostatically levitated droplets. Several corner treatments for the boundary and finite element coupling and their implications to free surface calculations are discussed. Detailed formulation and numerical implementation are given. Numerical results are compared with available analytical solutions whenever available. A selection of computed results is presented for mag‐ netically or electrostatically levitated droplets under both terrestrial and microgravity conditions. Copyright © 1999 John Wiley & Sons, Ltd.     

An effective hybrid displacement function element method for solving the edge effect of Mindlin–Reissner plate

In bending problems of Mindlin–Reissner plate, the resultant forces often vary dramatically within a narrow range near free and soft simply‐supported (SS1) boundaries. This is so‐called the edge effect or the boundary layer effect, a challenging problem for conventional finite element method. In this paper, an effective finite element method for analysis of such edge effect is developed. The construction procedure is based on the hybrid displacement function (HDF) element method [1], a simple hybrid‐Trefftz stress element method proposed recently. What is different is that an additional displacement function f related to the edge effect is considered, and its analytical solutions are employed as the additional trial functions for the first time. Furthermore, the free and the SS1 boundary conditions are also applied to modify the element assumed resultant fields. Then, two new special elements, HDF‐P4‐Free and HDF‐P4‐SS1, are successfully constructed. These new elements are allocated along the corresponding boundaries of the plate, while the other region is modeled by the usual HDF plate element HDF‐P4‐11 β [1]. Numerical tests demonstrate that the present method can effectively capture the edge effects and exactly satisfy the corresponding boundary conditions by only using relatively coarse meshes. Copyright © 2015 John Wiley & Sons, Ltd.     

Mesh refinement and iterative solution methods for finite element computations

Global and element residuals are introduced to determine a posteriori, computable, error bounds for finite element computations on a given mesh. The element residuals provide a criterion for determining where a finite element mesh requires refinement. This indicator is implemented in an algorithm in a finite element research program. There it is utilized to automatically refine the mesh for sample two-point problems exhibiting boundary layer and interior layer solutions. Results for both linear and nonlinear problems are presented. An important aspect of this investigation concerns the use of adaptive refinement in conjunction with iterative methods for system solution. As the mesh is being enriched through the refinement process, the solution on a given mesh provides an accurate starting iterate for the next mesh, and so on. A wide range of iterative methods are examined in a feasibility study and strategies for interweaving refinement and iteration are compared.