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.     

Stress recovery and error estimation for the scaled boundary finite‐element method

The scaled boundary finite‐element method is a novel semi‐analytical technique, combining the advantages of the finite element and the boundary element methods with unique properties of its own. This paper develops a stress recovery procedure based on a modal interpretation of the scaled boundary finite‐element method solution process, using the superconvergent patch recovery technique. The recovered stresses are superconvergent, and are used to calculate a recovery‐type error estimator. A key feature of the procedure is the compatibility of the error estimator with the standard recovery‐type finite element estimator, allowing the scaled boundary finite‐element method to be compared directly with the finite element method for the first time. A plane strain problem for which an exact solution is available is presented, both to establish the accuracy of the proposed procedures, and to demonstrate the effectiveness of the scaled boundary finite‐element method. The scaled boundary finite‐element estimator is shown to predict the true error more closely than the equivalent finite element error estimator. Unlike their finite element counterparts, the stress recovery and error estimation techniques work well with unbounded domains and stress singularities. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Interval spectral stochastic finite element analysis of structures with aggregation of random field and bounded parameters

This paper presents the study on non‐deterministic problems of structures with a mixture of random field and interval material properties under uncertain‐but‐bounded forces. Probabilistic framework is extended to handle the mixed uncertainties from structural parameters and loads by incorporating interval algorithms into spectral stochastic finite element method. Random interval formulations are developed based on K–L expansion and polynomial chaos accommodating the random field Young's modulus, interval Poisson's ratios and bounded applied forces. Numerical characteristics including mean value and standard deviation of the interval random structural responses are consequently obtained as intervals rather than deterministic values. The randomised low‐discrepancy sequences initialized particles and high‐order nonlinear inertia weight with multi‐dimensional parameters are employed to determine the change ranges of statistical moments of the random interval structural responses. The bounded probability density and cumulative distribution of the interval random response are then visualised. The feasibility, efficiency and usefulness of the proposed interval spectral stochastic finite element method are illustrated by three numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of singular stress fields at multi‐material corners under thermal loading

The scaled boundary finite‐element method is extended to the modelling of thermal stresses. The particular solution for the non‐homogeneous term caused by thermal loading is expressed as integrals in the radial direction, which are evaluated analytically for temperature changes varying as power functions of the radial coordinate. When applied to model a multi‐material corner, only the boundary of the problem domain is discretized. The boundary conditions on the straight material interfaces and the side‐faces forming the corner are satisfied analytically without discretization. The stress field is expressed semi‐analytically as a series solution. The stress distribution along the radial direction, including both the real and complex power singularity and the power‐logarithmic singularity, is represented analytically. The stress intensity factors are determined directly from their definitions in stresses. No knowledge on asymptotic expansions is required. Numerical examples are calculated to evaluate the accuracy of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd.     

A perturbation method for stochastic meshless analysis in elastostatics

A stochastic meshless method is presented for solving boundary‐value problems in linear elasticity that involves random material properties. The material property was modelled as a homogeneous random field. A meshless formulation was developed to predict stochastic structural response. Unlike the finite element method, the meshless method requires no structured mesh, since only a scattered set of nodal points is required in the domain of interest. There is no need for fixed connectivities between nodes. In conjunction with the meshless equations, classical perturbation expansions were derived to predict second‐moment characteristics of response. Numerical examples based on one‐ and two‐dimensional problems are presented to examine the accuracy and convergence of the stochastic meshless method. A good agreement is obtained between the results of the proposed method and Monte Carlo simulation. Since mesh generation of complex structures can be a far more time‐consuming and costly effort than the solution of a discrete set of equations, the meshless method provides an attractive alternative to finite element method for solving stochastic mechanics problems. Copyright © 2001 John Wiley & Sons, Ltd.     

A p‐hierarchical adaptive procedure for the scaled boundary finite element method

This paper describes a p‐hierarchical adaptive procedure based on minimizing the classical energy norm for the scaled boundary finite element method. The reference solution, which is the solution of the fine mesh formed by uniformly refining the current mesh element‐wise one order higher, is used to represent the unknown exact solution. The optimum mesh is assumed to be obtained when each element contributes equally to the global error. The refinement criteria and the energy norm‐based error estimator are described and formulated for the scaled boundary finite element method. The effectivity index is derived and used to examine quality of the proposed error estimator. An algorithm for implementing the proposed p‐hierarchical adaptive procedure is developed. Numerical studies are performed on various bounded domain and unbounded domain problems. The results reflect a number of key points. Higher‐order elements are shown to be highly efficient. The effectivity index indicates that the proposed error estimator based on the classical energy norm works effectively and that the reference solution employed is a high‐quality approximation of the exact solution. The proposed p‐hierarchical adaptive strategy works efficiently. Copyright © 2007 John Wiley & Sons, Ltd.     

eXtended Stochastic Finite Element Method for the numerical simulation of heterogeneous materials with random material interfaces

An eXtended Stochastic Finite Element Method has been recently proposed for the numerical solution of partial differential equations defined on random domains. This method is based on a marriage between the eXtended Finite Element Method and spectral stochastic methods. In this article, we propose an extension of this method for the numerical simulation of random multi‐phased materials. The random geometry of material interfaces is described implicitly by using random level set functions. A fixed deterministic finite element mesh, which is not conforming to the random interfaces, is then introduced in order to approximate the geometry and the solution. Classical spectral stochastic finite element approximation spaces are not able to capture the irregularities of the solution field with respect to spatial and stochastic variables, which leads to a deterioration of the accuracy and convergence properties of the approximate solution. In order to recover optimal convergence properties of the approximation, we propose an extension of the partition of unity method to the spectral stochastic framework. This technique allows the enrichment of approximation spaces with suitable functions based on an a priori knowledge of the irregularities in the solution. Numerical examples illustrate the efficiency of the proposed method and demonstrate the relevance of the enrichment procedure. Copyright © 2010 John Wiley & Sons, Ltd.     

Use of higher‐order shape functions in the scaled boundary finite element method

The scaled boundary finite element method is a novel semi‐analytical technique, whose versatility, accuracy and efficiency are not only equal to, but potentially better than the finite element method and the boundary element method for certain problems. This paper investigates the possibility of using higher‐order polynomial functions for the shape functions. Two techniques for generating the higher‐order shape functions are investigated. In the first, the spectral element approach is used with Lagrange interpolation functions. In the second, hierarchical polynomial shape functions are employed to add new degrees of freedom into the domain without changing the existing ones, as in the p‐version of the finite element method. To check the accuracy of the proposed procedures, a plane strain problem for which an exact solution is available is employed. A more complex example involving three scaled boundary subdomains is also addressed. The rates of convergence of these examples under p‐refinement are compared with the corresponding rates of convergence achieved when uniform h‐refinement is used, allowing direct comparison of the computational cost of the two approaches. The results show that it is advantageous to use higher‐order elements, and that higher rates of convergence can be obtained using p‐refinement instead of h‐refinement. Copyright © 2005 John Wiley & Sons, Ltd.     

Hybrid stress quadrilateral finite element approximation for stochastic plane elasticity equations

This paper considers stochastic hybrid stress quadrilateral finite element analysis of plane elasticity equations with stochastic Young's modulus and stochastic loads. Firstly, we apply Karhunen–Loève expansion to stochastic Young's modulus and stochastic loads so as to turn the original problem into a system containing a finite number of deterministic parameters. Then we deal with the stochastic field and the space field by k ?version/p ?version finite element methods and a hybrid stress quadrilateral finite element method, respectively. We derive a priori error estimates, which are uniform with respect to the Lamè constant λ ∈(0,+∞ ). Finally, we provide some numerical results. Copyright © 2016 John Wiley & Sons, Ltd.     

An adaptive remeshing procedure for discontinuous finite element limit analysis

An adaptive remeshing procedure is proposed for discontinuous finite element limit analysis. The procedure proceeds by iteratively adjusting the element sizes in the mesh to distribute local errors uniformly over the domain. To facilitate the redefinition of element sizes in the new mesh, the interelements discontinuous field of elemental bound gaps is converted into a continuous field, ie, the intensity of bound gap, using a patch‐based approximation technique. An analogous technique is subsequently used for the approximation of element sizes in the old mesh. With these information, an optimized distribution of element sizes in the new mesh is defined and then scaled to match the total number of elements specified for each iteration in the adaptive remeshing process. Finally, a new mesh is generated using the advancing front technique. This adaptive remeshing procedure is repeated several times until an optimal mesh is found. Additionally, for problems involving discontinuous boundary loads, a novel algorithm for the generation of fan‐type meshes around singular points is proposed explicitly and incorporated into the main adaptive remeshing procedure. To demonstrate the feasibility of our proposed method, some classical examples extracted from the existing literary works are studied in detail.     

A Frobenius solution to the scaled boundary finite element equations in frequency domain for bounded media

The scaled boundary finite element method (FEM) is a recently developed semi‐analytical numerical approach combining advantages of the FEM and the boundary element method. Although for elastostatics, the governing homogeneous differential equations in the radial co‐ordinate can be solved analytically without much effort, an analytical solution to the non‐homogeneous differential equations in frequency domain for elastodynamics has so far only been obtained by a rather tedious series‐expansion procedure. This paper develops a much simpler procedure to obtain such an analytical solution by increasing the number of power series in the solution until the required accuracy is achieved. The procedure is applied to an extensive study of the steady‐state frequency response of a square plate subjected to harmonic excitation. Comparison of the results with those obtained using ABAQUS shows that the new method is as accurate as a detailed finite element model in calculating steady‐state responses for a wide range of frequencies using only a fraction of the degrees of freedom required in the latter. Copyright © 2006 John Wiley & Sons, Ltd.     

Dynamic soil–structure interaction analysis of layered unbounded media via a coupled finite element/boundary element/scaled boundary finite element model

In this paper, a coupled model based on finite element method (FEM), boundary element method (BEM) and scaled boundary FEM (SBFEM) (also referred to as the consistent infinitesimal finite element cell method) for dynamic response of 2D structures resting on layered soil media is presented. The SBFEM proposed by Wolf and Song (Finite‐element Modelling of Unbounded Media. Wiley: England, 1996) and BEM are used for modelling the dynamic response of the unbounded media (far‐field). The standard FEM is used for modelling the finite region (near‐field) and the structure. In SBFEM, which is a semi‐analytical technique, the radiation condition at infinity is satisfied exactly without requiring the fundamental solution. This method, also eliminates the need for the discretization of interfaces between different layers. In both SBFEM and BEM, the spatial dimension is decreased by one. The objective of the development of this coupled model is to combine advantages of above‐mentioned three numerical models to solve various soil–structure interaction (SSI) problems efficiently and effectively. These three methods are coupled (FE–BE–SBFEM) via substructuring method, and a computer programme is developed for the harmonic analyses of SSI systems. The coupled model is established in such a way that, depending upon the problem and far‐field properties, one can choose BEM and/or SBFEM in modelling related far‐field region(s). Thus, BEM and/or SBFEM can be used efficiently in modelling the far‐field. The proposed model is applied to investigate dynamic response of rigid and elastic structures resting on layered soil media. To assess the proposed SSI model, several problems existing in the literature are chosen and analysed. The results of the proposed model agree with the results presented in the literature for the chosen problems. The advantages of the model are demonstrated through these comparisons. Copyright © 2004 John Wiley & Sons, Ltd.     

Stochastic nonlinear free vibration of laminated composite plates resting on elastic foundation in thermal environments

This paper presents the nonlinear free vibration analysis of laminated composite plates resting on elastic foundation with random system properties in thermal environments. System parameters are modeled as basic random variables for accurate prediction of system behavior. A C ⁰ nonlinear finite element based on HSDT in von Karman sense is used to descretize the laminate. A direct iterative method in conjunction with first-order perturbation technique is outlined and applied to solve the stochastic nonlinear generalized eigenvalue problem. The developed stochastic procedure is successfully used for thermally induced nonlinear free vibration problem with a reasonable accuracy. Numerical results for various combinations of boundary conditions, geometric parameters, amplitude ratios, foundation parameters and thermal loading have been compared with those available in literature and an independent MCS. Some new results are also presented which clearly demonstrate the importance of the randomness in the system parameters on the response of the structures.     

Computation of 3-D stress singularities for multiple cracks and crack intersections by the scaled boundary finite element method

Among the various possible ways of dealing with notch and crack situations, the scaled boundary finite element method [SBFEM, (Wolf and Song in Finite element modelling of unbounded structures. Wiley, Chichester, 1996; Wolf in The scaled boundary finite element method. Wiley, Chichester, 2003)] has been adopted in this work. This method has been proved to be versatile, much less time consuming than the finite element method and generates highly accurate numerical predictions in cases of structures with notches and cracks. The SBFEM gives the advantage of boundary element method by reducing one dimension in modelling the structures but the mathematical formulations are more related to conventional displacement based finite element method. This method requires a certain scalability of the given structure with respect to a point called similarity center. Like in the case of the boundary element method, the structure needs to be discretized only at the surface where standard displacement based isoparametric finite element formulations are adequate. Unlike in the boundary element method, however, no fundamental solution is required by the scaled boundary finite element method. The similarity or scalability of the method requires separation of coordinates such that in the radial direction (i.e. scaling direction) it yields simple differential equations that can be solved analytically. So this approach can be considered as a semi-analytical method. Several two-dimensional examples have been analysed for crack and notch situations that are well known cases in fracture mechanics. A number of three-dimensional cases have been considered for different crack configurations that yield high order of singularity. The results, according to the authors’ knowledge are up to now unpublished in the open literature. Parametric studies are conducted for structures with bi-material interfaces.     

Maximum likelihood estimation of stochastic chaos representations from experimental data

This paper deals with the identification of probabilistic models of the random coefficients in stochastic boundary value problems (SBVP). The data used in the identification correspond to measurements of the displacement field along the boundary of domains subjected to specified external forcing. Starting with a particular mathematical model for the mechanical behaviour of the specimen, the unknown field to be identified is projected on an adapted functional basis such as that provided by a finite element discretization. For each set of measurements of the displacement field along the boundary, an inverse problem is formulated to calculate the corresponding optimal realization of the coefficients of the unknown random field on the adapted basis. Realizations of these coefficients are then used, in conjunction with the maximum likelihood principle, to set‐up and solve an optimization problem for the estimation of the coefficients in a polynomial chaos representation of the parameters of the SBVP. Copyright © 2005 John Wiley & Sons, Ltd.     

A definition and evaluation procedure of generalized stress intensity factors at cracks and multi-material wedges

This paper proposes a definition of generalized stress intensity factors that includes classical definitions for crack problems as special cases. Based on the semi-analytical solution obtained from the scaled boundary finite-element method, the singular stress field is expressed as a matrix power function with its dimension equal to the number of singular terms. Not only real and complex power singularities but also power-logarithmic singularities are represented in a unified expression without explicitly determining the type of singularity. The generalized stress intensity factors are evaluated directly from the scaled boundary finite-element solution for the singular stress field by following standard stress recovery procedures in the finite element method. The definition and evaluation procedure are valid to multi-material wedges composed of any number of isotropic and anisotropic materials. Numerical examples, including a cracked homogeneous plate, a bimaterial plate with an interfacial crack, a V-notched bimaterial plate and a crack terminating at a material interface, are analyzed. Features of this unified definition are discussed.     

Lower‐bound limit analysis by using the EFG method and non‐linear programming

Intended to avoid the complicated computations of elasto‐plastic incremental analysis, limit analysis is an appealing direct method for determining the load‐carrying capacity of structures. On the basis of the static limit analysis theorem, a solution procedure for lower‐bound limit analysis is presented firstly, making use of the element‐free Galerkin (EFG) method rather than traditional numerical methods such as the finite element method and boundary element method. The numerical implementation is very simple and convenient because it is only necessary to construct an array of nodes in the domain under consideration. The reduced‐basis technique is adopted to solve the mathematical programming iteratively in a sequence of reduced self‐equilibrium stress subspaces with very low dimensions. The self‐equilibrium stress field is expressed by a linear combination of several self‐equilibrium stress basis vectors with parameters to be determined. These self‐equilibrium stress basis vectors are generated by performing an equilibrium iteration procedure during elasto‐plastic incremental analysis. The Complex method is used to solve these non‐linear programming sub‐problems and determine the maximal load amplifier. Numerical examples show that it is feasible and effective to solve the problems of limit analysis by using the EFG method and non‐linear programming. Copyright © 2007 John Wiley & Sons, Ltd.     

A modified scaled boundary finite element method for three‐dimensional dynamic soil‐structure interaction in layered soil

This paper is devoted to the analysis of elastodynamic problems in 3D‐layered systems which are unbounded in the horizontal direction. For this purpose, a finite element model of the near field is coupled to a scaled boundary finite element model (SBFEM) of the far field. The SBFEM is originally based on describing the geometry of a half‐space or full‐space domain by scaling the geometry of the near field / far field interface using a radial coordinate. A modified form of the SBFEM for waves in a 2D layer is also available. None of these existing formulations can be used to describe a 3D‐layered medium. In this paper, a modified SBFEM for the analysis of 3D‐layered continua is derived. Based on the use of a scaling line instead of a scaling centre, a suitable scaled boundary transformation is proposed. The derivation of the corresponding scaled boundary finite element (SBFE) equations in displacement and stiffness is presented in detail. The latter is a nonlinear differential equation with respect to the radial coordinate, which has to be solved numerically for each excitation frequency considered in the analysis. Various numerical examples demonstrate the accuracy of the new method and its correct implementation. These include rigid circular and square foundations embedded in or resting on the surface of layered homogeneous or inhomogeneous 3D soil deposits over rigid bedrock. Hysteretic damping is assumed in some cases. The dynamic stiffness coefficients calculated using the proposed method are compared with analytical solutions or existing highly accurate numerical results. Copyright © 2011 John Wiley & Sons, Ltd.