A hierarchical finite element for geometrically non‐linear vibration of doubly curved,moderately thick isotropic shallow shells

A p‐version, hierarchical finite element for doubly curved, moderately thick, isotropic shallow shells is derived and geometrically non‐linear free vibrations of panels with rectangular planform are investigated. The geometrical non‐linearity is due to large displacements, and the effects of the rotatory inertia and transverse shear are considered. The time domain equations of motion are obtained by applying the principle of virtual work and the d'Alembert's principle. These equations are mapped to the frequency domain by the harmonic balance method, and are finally solved by a predictor–corrector method. The convergence properties of the element proposed and the influence of several parameters on the dynamic response are studied. These parameters are the shell's thickness, the width‐to‐length ratio, the curvature‐to‐width ratio and the ratio between curvature radii. The first and higher order modes are analysed. Some results are compared with results published or calculated using a commercial finite element package. It is demonstrated that with the proposed element low‐dimensional, accurate models are obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

An accurate two-node finite element for shear deformable curved beams

An accurate two-node (three degrees of freedom per node) finite element is developed for curved shear deformable beams. The element formulation is based on shape functions that satisfy the homogeneous form of the partial differential equations of motion which renders it free of shear and membrane locking. The element is demonstrated to converge to the results obtained from a shear deformable straight beam when the beam becomes shallower. Numerical examples were performed to demonstrate the accuracy and efficiency with respect to previously published formulations. © 1998 John Wiley & Sons, Ltd.     

Geometrically non‐linear thermoelastic analysis of functionally graded shells using finite element method

A finite element formulation governing the geometrically non‐linear thermoelastic behaviour of plates and shells made of functionally graded materials is derived in this paper using the updated Lagrangian approach. Derivation of the formulation is based on rewriting the Green–Lagrange strain as well as the 2nd Piola–Kirchhoff stress as two second‐order functions in terms of a through‐the‐thickness parameter. Material properties are assumed to vary through the thickness according to the commonly used power law distribution of the volume fraction of the constituents. Within a non‐linear finite element analysis framework, the main focus of the paper is the proposal of a formulation to account for non‐linear stress distribution in FG plates and shells, particularly, near the inner and outer surfaces for small and large values of the grading index parameter. The non‐linear heat transfer equation is also solved for thermal distribution through the thickness by the Rayleigh–Ritz method. Advantages of the proposed approach are assessed and comparisons with available solutions are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

The p‐version of the finite element method for three‐dimensional curved thin walled structures

In this paper we present an implementation of a three‐dimensional p‐version for structural problems of solids with almost arbitrarily curved surfaces. Applying the blending function method, complex structures can often be modelled by a few p‐elements, being the basis for a higher order approximation. Numerical examples will demonstrate, that the p‐version with anisotropic Ansatz spaces allows to predict the structural behaviour of three‐dimensional plates and shells with approximately the same amount of degrees of freedom as in the two‐dimensional case, yet significantly more accurate due to the three‐dimensional model. Furthermore, it is advantageous to compute complex structures exclusively with three‐dimensional discretizations as no special elements are needed to model the transition from dimensionally reduced formulations like plates or shells to fully three‐dimensional solid elements. Using the p‐version with anisotropic Ansatz spaces the whole structure can be efficiently discretized with solid elements, even if the aspect ratio of the elements becomes very large. Copyright © 2001 John Wiley Sons, Ltd.     

Vibration analysis of beams with non‐local foundations using the finite element method

In this paper, a non‐local viscoelastic foundation model is proposed and used to analyse the dynamics of beams with different boundary conditions using the finite element method. Unlike local foundation models the reaction of the non‐local model is obtained as a weighted average of state variables over a spatial domain via convolution integrals with spatial kernel functions that depend on a distance measure. In the finite element analysis, the interpolating shape functions of the element displacement field are identical to those of standard two‐node beam elements. However, for non‐local elasticity or damping, nodes remote from the element do have an effect on the energy expressions, and hence the damping and stiffness matrices. The expressions of these direct and cross‐matrices for stiffness and damping may be obtained explicitly for some common spatial kernel functions. Alternatively numerical integration may be applied to obtain solutions. Numerical results for eigenvalues and associated eigenmodes of Euler–Bernoulli beams are presented and compared (where possible) with results in literature using exact solutions and Galerkin approximations. The examples demonstrate that the finite element technique is efficient for the dynamic analysis of beams with non‐local viscoelastic foundations. Copyright © 2007 John Wiley & Sons, Ltd.     

A refined non‐linear non‐conforming triangular plate/shell element

A refined non‐conforming triangular plate/shell element for geometric non‐linear analysis of plates/shells using the total Lagrangian/updated Lagrangian approach is constructed in this paper based on the refined non‐conforming element method for geometric non‐linear analysis. The Allman's triangular plane element with vertex degrees of freedom and the refined triangular plate‐bending element RT9 are used to construct the present element. Numerical examples demonstrate that the accuracy of the new element is quite high in the geometric non‐linear analysis of plates/shells. Copyright © 2003 John Wiley & Sons, Ltd.     

Element‐wise a posteriori estimates based on hierarchical bases for non‐linear parabolic problems

We show that the issue of a posteriori estimate the errors in the numerical simulation of non‐linear parabolic equations can be reduced to a posteriori estimate the errors in the approximation of an elliptic problem with the right‐hand side depending on known data of the problem and the computed numerical solution. A procedure to obtain local error estimates for the p version of the finite element method by solving small discrete elliptic problems with right‐hand side the residual of the p‐FEM solution is introduced. The boundary conditions are inherited by those of the space of hierarchical bases to which the error estimator belongs. We prove that the error in the numerical solution can be reduced by adding the estimators that behave as a locally defined correction to the computed approximation. When the error being estimated is that of a elliptic problem constant free local lower bounds are obtained. The local error estimation procedure is applied to non‐linear parabolic differential equations in several space dimensions. Some numerical experiments for both the elliptic and the non‐linear parabolic cases are provided. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite element of slender beams in finite transformations: a geometrically exact approach

This article is devoted to the modelling of thin beams undergoing finite deformations essentially due to bending and torsion and to their numerical resolution by the finite element method. The solution proposed here differs from the approaches usually implemented to treat thin beams, as it can be qualified as ‘geometrically exact’. Two numerical models are proposed. The first one is a non‐linear Euler–Bernoulli model while the second one is a non‐linear Rayleigh model. The finite element method is tested on several numerical examples in statics and dynamics, and validated through comparison with analytical solutions, experimental observations and the geometrically exact approach of the Reissner beam theory initiated by Simo. The numerical result shows that this approach is a good alternative to the modelling of non‐linear beams, especially in statics. Copyright © 2003 John Wiley & Sons, Ltd.     

A semi‐analytical curved element for linear elasticity based on the scaled boundary finite element method

This work introduces a semi‐analytical formulation for the simulation and modeling of curved structures based on the scaled boundary finite element method (SBFEM). This approach adapts the fundamental idea of the SBFEM concept to scale a boundary to describe a geometry. Until now, scaling in SBFEM has exclusively been performed along a straight coordinate that enlarges, shrinks, or shifts a given boundary. In this novel approach, scaling is based on a polar or cylindrical coordinate system such that a boundary is shifted along a curved scaling direction. The derived formulations are used to compute the static and dynamic stiffness matrices of homogeneous curved structures. The resulting elements can be coupled to general SBFEM or FEM domains. For elastodynamic problems, computations are performed in the frequency domain. Results of this work are validated using the global matrix method and standard finite element analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

Comparison of high‐order curved finite elements

Several finite element techniques used in domains with curved boundaries are discussed and compared, with particular emphasis in two issues: the exact boundary representation of the domain and the consistency of the approximation. The influence of the number of integration points on the accuracy of the computation is also studied. Two‐dimensional numerical examples, solved with continuous and discontinuous Galerkin formulations, are used to test and compare all these methodologies. In every example shown, the recently proposed NURBS‐enhanced finite element method (NEFEM) provides the maximum accuracy for a given spatial discretization, at least one order of magnitude more accurate than classical isoparametric finite element method (FEM). Moreover, NEFEM outperforms Cartesian FEM and p‐FEM, stressing the importance of the geometrical model as well as the relevance of a consistent approximation in finite element simulations. Copyright © 2011 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.     

Numerical finite element formulation of the Schapery non‐linear viscoelastic material model

This study presents a numerical integration method for the non‐linear viscoelastic behaviour of isotropic materials and structures. The Schapery's three‐dimensional (3D) non‐linear viscoelastic material model is integrated within a displacement‐based finite element (FE) environment. The deviatoric and volumetric responses are decoupled and the strain vector is decomposed into instantaneous and hereditary parts. The hereditary strains are updated at the end of each time increment using a recursive formulation. The constitutive equations are expressed in an incremental form for each time step, assuming a constant incremental strain rate. A new iterative procedure with predictor–corrector type steps is combined with the recursive integration method. A general polynomial form for the parameters of the non‐linear Schapery model is proposed. The consistent algorithmic tangent stiffness matrix is realized and used to enhance convergence and help achieve a correct convergent state. Verifications of the proposed numerical formulation are performed and compared with a previous work using experimental data for a glassy amorphous polymer PMMA. Copyright © 2003 John Wiley & Sons, Ltd.     

Energy and momentum conserving elastodynamics of a non‐linear brick‐type mixed finite shell element

The paper presents aspects of the finite element formulation of momentum and energy conserving algorithms for the non‐linear dynamic analysis of shell‐like structures. The key contribution is a detailed analysis of the implementation of a Simó–Tarnow‐type conservation scheme in a recently developed new mixed finite shell element. This continuum‐based shell element provides a well‐defined interface to strain‐driven constitutive stress updates algorithms. It is based on the classic brick‐type trilinear displacement element and is equipped with specific gradient‐type enhanced strain modes and shell‐typical assumed strain modifications. The excellent performance of the proposed dynamic shell formulation with respect to conservation properties and numerical stability behaviour is demonstrated by means of three representative numerical examples of elastodynamics which exhibit complex free motions of flexible structures undergoing large strains and large rigid‐body motions. 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.     

Calculation of strict error bounds for finite element approximations of non‐linear pointwise quantities of interest

This paper deals with the verification of simulations performed using the finite element method. More specifically, it addresses the calculation of strict bounds on the discretization errors affecting pointwise outputs of interest which may be non‐linear with respect to the displacement field. The method is based on classical tools, such as the constitutive relation error and extraction techniques associated with the solution of an adjoint problem. However, it uses two specific and innovative techniques: the enrichment of the adjoint solution using a partition of unity method, which enables one to consider truly pointwise quantities of interest, and the decomposition of the non‐linear quantities of interest by means of projection properties in order to take into account higher‐order terms in establishing the bounds. Thus, no linearization is performed and the property that the local error bounds are guaranteed is preserved. The effectiveness of the approach and the quality of the bounds are illustrated with two‐dimensional applications in the context of elastic fatigue problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Adaptive superposition of finite element meshes in non‐linear transient solid mechanics problems

An s‐adaptive finite element procedure is developed for the transient analysis of 2‐D solid mechanics problems with material non‐linearity due to progressive damage. The resulting adaptive method simultaneously estimates and controls both the spatial error and temporal error within user‐specified tolerances. The spatial error is quantified by the Zienkiewicz–Zhu error estimator and computed via superconvergent patch recovery, while the estimation of temporal error is based on the assumption of a linearly varying third‐order time derivatives of the displacement field in conjunction with direct numerical time integration. The distinguishing characteristic of the s‐adaptive procedure is the use of finite element mesh superposition (s‐refinement) to provide spatial adaptivity. Mesh superposition proves to be particularly advantageous in computationally demanding non‐linear transient problems since it is faster, simpler and more efficient than traditional h‐refinement schemes. Numerical examples are provided to demonstrate the performance characteristics of the s‐adaptive method for quasi‐static and transient problems with material non‐linearity. Copyright © 2007 John Wiley & Sons, Ltd.     

An a posteriori error estimator for the p‐ and hp‐versions of the finite element method

An a posteriori error estimator is proposed in this paper for the p‐ and hp‐versions of the finite element method in two‐dimensional linear elastostatic problems. The local error estimator consists in an enhancement of an error indicator proposed by Bertóti and Szabó (Int. J. Numer. Meth. Engng. 1998; 42 :561–587), which is based on the minimum complementary energy principle. In order to obtain the local error estimate, this error indicator is corrected by a factor which depends only on the polynomial degree of the element. The proposed error estimator shows a good effectivity index in meshes with uniform and non‐uniform polynomial distributions, especially when the global error is estimated. Furthermore, the local error estimator is reliable enough to guide p‐ and hp‐adaptive refinement strategies. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical p‐version refinement studies for the regularized stress‐BEM

In the development of the boundary element method (BEM) and the finite element method (FEM) researchers have typically selected similar basis functions. That is, both methods typically employ low‐order interpolations such as piece‐wise linear or piece‐wise quadratic and rely on h‐version refinement to increase accuracy as required. In the case of the FEM, the decision to use low‐order elements is made for computational efficiency as an attractive compromise between local modeling accuracy and sparseness of the resulting linear system. However, in many BEM formulations, low‐order elements may be the only practical choice given the complexity of using analytic integration formulae in conjunction with special integral interpretations. Unlike their efficient use in the FEM, fine meshes of low‐order elements in the BEM are highly inefficient from a computational standpoint given the dense nature of BEM systems. Moreover, owing to singularities in the kernel functions, the BEM should be expected to benefit more so than the FEM from very high levels of local accuracy. Through the use of regularized algorithms which only require numerical integration, p‐version refinement in the BEM is easily extended to include any set of basis functions with no significant increase in programming complexity. Numerical results show that by using interpolations as high as 12th and 16th order, one can expect reductions in error by as many as five orders of magnitude over comparable algorithms based on similar system size. For two‐dimensional problems, it is also shown that, for a given level of error, one can expect reductions in system size by an order of magnitude, thus leading to a reduction in computational expense for conventional algorithms by three orders of magnitude. Copyright © 2003 John Wiley & Sons, Ltd.