c‐Type method of unified CAMG and FEA. Part 1: Beam and arch mega‐elements—3D linear and 2D non‐linear

Computer‐aided mesh generation (CAMG) dictated solely by the minimal key set of requirements of geometry, material, loading and support condition can produce ‘mega‐sized’, arbitrary‐shaped distorted elements. However, this may result in substantial cost saving and reduced bookkeeping for the subsequent finite element analysis (FEA) and reduced engineering manpower requirement for final quality assurance. A method, denoted as c‐type, has been proposed by constructively defining a finite element space whereby the above hurdles may be overcome with a minimal number of hyper‐sized elements. Bezier (and de Boor) control vectors are used as the generalized displacements and the Bernstein polynomials (and B‐splines) as the elemental basis functions. A concomitant idea of coerced parametry and inter‐element continuity on demand unifies modelling and finite element method. The c‐type method may introduce additional control, namely, an inter‐element continuity condition to the existing h‐type and p‐type methods. Adaptation of the c‐type method to existing commercial and general‐purpose computer programs based on a conventional displacement‐based finite element method is straightforward. The c‐type method with associated subdivision technique can be easily made into a hierarchic adaptive computer method with a suitable a posteriori error analysis. In this context, a summary of a geometrically exact non‐linear formulation for the two‐dimensional curved beams/arches is presented. Several beam problems ranging from truly three‐dimensional tortuous linear curved beams to geometrically extremely non‐linear two‐dimensional arches are solved to establish numerical efficiency of the method. Incremental Lagrangian curvilinear formulation may be extended to overcome rotational singularity in 3D geometric non‐linearity and to treat general material non‐linearity. Copyright © 2003 John Wiley & Sons, Ltd.     

Finite‐element model sensitivity in the vibration of partially embedded beams

Beams partially embedded in a Winkler foundation can exhibit the phenomenon of modal clustering. In this paper, a conventional finite‐element package is used in the comparison of the natural frequencies of such beams as calculated using a discretized finite‐element model and as computed from a closed‐form exact solution. The sensitivity of finite‐element models of varying degrees of refinement in predicting the incidence of modal clustering is discussed. It is shown that, in crude models, the clustering of modes is predicted accurately but other modes, even the lowest ones, can be unacceptably distorted. In very crude models, the finite‐element analysis may fail altogether in identifying the modal‐clustering phenomenon. Copyright © 1999 John Wiley & Sons, Ltd.     

Exact static element stiffness matrices of non‐symmetric thin‐walled curved beams

An evaluation procedure of exact static stiffness matrices for curved beams with non‐symmetric thin‐walled cross section are rigorously presented for the static analysis. Higher‐order differential equations for a uniform curved beam element are first transformed into a set of the first‐order simultaneous ordinary differential equations by introducing 14 displacement parameters where displacement modes corresponding to zero eigenvalues are suitably taken into account. This numerical technique is then accomplished via a generalized linear eigenvalue problem with non‐symmetric matrices. Next, the displacement functions of displacement parameters are exactly calculated by determining general solutions of simultaneous non‐homogeneous differential equations. Finally an exact stiffness matrix is evaluated using force–deformation relationships. In order to demonstrate the validity and effectiveness of this method, displacements and normal stresses of cantilever thin‐walled curved beams subjected to tip loads are evaluated and compared with those by thin‐walled curved beam elements as well as shell elements. Copyright © 2004 John Wiley & Sons, Ltd.     

A spatially curved‐beam element with warping and Wagner effects

This paper presents a new spatially curved‐beam element with warping and Wagner effects that can be used for the non‐linear large displacement analysis of members that are curved in space. The non‐linear behaviour of members curved in space shows that the Wagner effects are substantial in the large twist rotation analysis. Most existing finite beam element models, such as ABAQUS and ANSYS cannot predict the non‐linear large displacement response of members curved in space correctly because the Wagner effects, viz. the Wagner moment and the corresponding finite strain terms, have not been considered in these finite beam elements. As a consequence, these finite beam elements do not provide correct predictions for the out‐of‐plane buckling and postbuckling behaviour of arches as well. In this paper, the symmetric tangent stiffness matrix has been derived based on the finite rotations parameterized by the conventional displacements. The warping and Wagner effects: both the Wagner moment and the corresponding finite strain terms and their constitutive relationship, are included in the spatially curved‐beam element. Two components of the initial curvature, the initial twist and their interactions with the displacements are also considered in the spatially curved‐beam element. This ensures that the large twist rotation analysis for the members curved in space is accurate. Comparisons with existing experimental, analytical and numerical results show that the spatially curved‐beam element is accurate and efficient for the non‐linear elastic analysis of curved members, buckling and postbuckling analysis of arches, and in its ability to predict large deflections and twist rotations in more arbitrarily curved members. Copyright © 2005 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.     

Free vibrations of solid and hollow wedge beams with rectangular or circular cross‐sections and carrying any number of point masses

Since the literature relating to the natural frequencies and mode shapes of the double‐tapered wedge beams carrying multiple point masses is rare, the object of this paper is to present some information in this aspect. First of all, the closed‐form solutions in terms of the Bessel functions for the natural frequencies and normal mode shapes of the ‘bare’ wedge beams (without carrying any point masses) were determined. Next, the partial differential equation of motion for the ‘loading’ wedge beams (carrying any number of point masses) were transformed into the matrix equation by using the expansion theorem and the foregoing natural frequencies and normal mode shapes of the ‘bare’ wedge beam. Finally, the eigenvalue equation associated with the last matrix equation was solved to give the natural frequencies and the mode shapes of the ‘loading’ wedge beams. The formulation of this paper is available for the solid and hollow wedge beams with square, rectangular or circular cross sections. In other words, the taper ratio for the width and that for the depth may be equal or unequal. All the numerical results were compared with the existing literature or the conventional finite element method results and good agreement was achieved. Copyright © 2004 John Wiley & Sons, Ltd.     

Theory and numerics of three‐dimensional beams with elastoplastic material behaviour

A theory of space curved beams with arbitrary cross‐sections and an associated finite element formulation is presented. Within the present beam theory the reference point, the centroid, the centre of shear and the loading point are arbitrary points of the cross‐section. The beam strains are based on a kinematic assumption where torsion‐warping deformation is included. Each node of the derived finite element possesses seven degrees of freedom. The update of the rotational parameters at the finite element nodes is achieved in an additive way. Applying the isoparametric concept the kinematic quantities are approximated using Lagrangian interpolation functions. Since the reference curve lies arbitrarily with respect to the centroid the developed element can be used to discretize eccentric stiffener of shells. Due to the implemented constitutive equations for elastoplastic material behaviour the element can be used to evaluate the load‐carrying capacity of beam structures. Copyright © 2000 John Wiley & Sons, Ltd.     

A rational elasto‐plastic spatially curved thin‐walled beam element

Torsion is one of the primary actions in members curved in space, and so an accurate spatially curved‐beam element needs to be able to predict the elasto‐plastic torsional behaviour of such members correctly. However, there are two major difficulties in most existing finite thin‐walled beam elements, such as in ABAQUS and ANSYS, which may lead to incorrect predictions of the elasto‐plastic behaviour of members curved in space. Firstly, the integration sample point scheme cannot capture the shear strain and stress information resulting from uniform torsion. Secondly, the higher‐order twists are ignored which leads to loss of the significant effects of Wagner moments on the large twist torsional behaviour. In addition, the initial geometric imperfections and residual stresses are significant for the elasto‐plastic behaviour of members curved in space. Many existing finite thin‐walled beam element models do not provide facilities to deal with initial geometric imperfections. Although ABAQUS and ANSYS have facilities for the input of residual stresses as initial stresses, they cannot describe the complicated distribution patterns of residual stresses in thin‐walled members. Furthermore, external loads and elastic restraints may be applied remote from shear centres or centroids. The effects of the load (and restraint) positions are important, but are not considered in many beam elements. This paper presents an elasto‐plastic spatially curved element with arbitrary thin‐walled cross‐sections that can correctly capture the uniform shear strain and stress information for integration, and includes initial geometric imperfections, residual stresses and the effects of the load and restraint positions. The element also includes elastic restraints and supports, which have to be modelled separately as spring elements in some other finite thin‐walled beam elements. Comparisons with existing experimental and analytical results show that the elasto‐plastic spatially curved‐beam element is accurate and efficient. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Finite element formulation for anisotropic coupled piezoelectro‐hygro‐thermo‐viscoelasto‐dynamic problems

A finite element algorithm has been developed for the efficient analysis of smart composite structures with piezoelectric polymer sensors or/and actuators based on piezoelectro‐hygro‐thermo‐viscoelasticity. Variational principles for anisotropic coupled piezoelectro‐hygro‐thermo‐viscoelasto‐dynamic problems have also been proposed in this study. As illustrative studies, dynamic responses in laminated composite beams and plates with PVDF sensors and actuators are obtained as functions of time using the present finite element procedures. The voltage feedback control scheme is utilized. The proposed numerical method can be used for analysing problems in the design of smart structures as well as smart sensors and actuators. Copyright © 1999 John Wiley & Sons, Ltd.     

Free vibrations of shear‐flexible and compressible arches by FEM

The purpose of this paper is to analyse free vibrations of arches with influence of shear and axial forces taken into account. Arches with various depth of cross‐section and various types of supports are considered. In the calculations, the curved finite element elaborated by the authors is adopted. It is the plane two‐node, six‐degree‐of‐freedom arch element with constant curvature. Its application to the static analysis yields the exact results, coinciding with the analytical ones. This feature results from the use of the exact shape functions in derivation of the element stiffness matrix. In the free vibration analysis the consistent mass matrix is used. It is obtained on the base of the same functions. Their coefficients contain the influences of shear flexibility and compressibility of the arch. The numerical results are compared with the results obtained for the simple diagonal mass matrix representing the lumped mass model. The natural frequencies are also compared with the ones for the continuous arches for which the analytically determined frequencies are known. The advantage of the paper is a thorough analysis of selected examples, where the influences of shear forces, axial forces as well as the rotary and tangential inertia on the natural frequencies are examined. Copyright © 2001 John Wiley & Sons, Ltd.     

Modelling of reinforced‐concrete structures providing crack‐spacing based on X‐FEM,ED‐FEM and novel operator split solution procedure

In this work, we present a novel approach to the finite element modelling of reinforced‐concrete (RC) structures that provides the details of the constitutive behavior of each constituent (concrete, steel and bond‐slip), while keeping formally the same appearance as the classical finite element model. Each component constitutive behavior can be brought to fully non‐linear range, where we can consider cracking (or localized failure) of concrete, the plastic yielding and failure of steel bars and bond‐slip at concrete steel interface accounting for confining pressure effects. The standard finite element code architecture is preserved by using embedded discontinuity (ED‐FEM) and extended (X‐FEM) finite element strain representation for concrete and slip, respectively, along with the operator split solution method that separates the problem into computing the deformations of RC (with frozen slip) and the current value of the bond‐slip. Several numerical examples are presented in order to illustrate very satisfying performance of the proposed methodology. Copyright © 2010 John Wiley & Sons, Ltd.     

hp‐Generalized FEM and crack surface representation for non‐planar 3‐D cracks

A high‐order generalized finite element method (GFEM) for non‐planar three‐dimensional crack surfaces is presented. Discontinuous p‐hierarchical enrichment functions are applied to strongly graded tetrahedral meshes automatically created around crack fronts. The GFEM is able to model a crack arbitrarily located within a finite element (FE) mesh and thus the proposed method allows fully automated fracture analysis using an existing FE discretization without cracks. We also propose a crack surface representation that is independent of the underlying GFEM discretization and controlled only by the physics of the problem. The representation preserves continuity of the crack surface while being able to represent non‐planar, non‐smooth, crack surfaces inside of elements of any size. The proposed representation also provides support for the implementation of accurate, robust, and computationally efficient numerical integration of the weak form over elements cut by the crack surface. Numerical simulations using the proposed GFEM show high convergence rates of extracted stress intensity factors along non‐planar curved crack fronts and the robustness of the method. Copyright © 2008 John Wiley & Sons, Ltd.     

Study of warping torsion of thin‐walled beams with open cross‐section using macro‐elements

Thin‐walled beams with open cross‐section under torsion or complex load are studied based on the hypotheses of the classical theory (Vlasov). Different from previous techniques presented in the literature, the concept of a strip‐plate is introduced. This concept is used to accurately model the effect of bending induced by torsion and to define an alternate finite element called macro‐element. The macro‐elements are shown to model more accurately the thin‐walled beams under warping torsion or complex load therefore giving better results than the classical theory. Copyright © 1999 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.     

Coupled mechanics and finite element for non‐linear laminated piezoelectric shallow shells undergoing large displacements and rotations

A theoretical framework is presented for analysing the coupled non‐linear response of shallow doubly curved adaptive laminated piezoelectric shells undergoing large displacements and rotations. The formulated mechanics incorporate coupling between in‐plane and flexural stiffness terms due to geometric curvature, coupling between mechanical and electric fields, and encompass geometric non‐linearity effects due to large displacements and rotations. The governing equations are formulated explicitly in orthogonal curvilinear co‐ordinates and are combined with the kinematic assumptions of a mixed‐field shear‐layerwise shell laminate theory. Based on the above formulation, a finite element methodology together with an incremental‐iterative technique, based on Newton–Raphson method is formulated. An eight‐node coupled non‐linear shell element is also developed. Various evaluation cases on laminated curved beams and cylindrical panels illustrate the capability of the shell finite element to predict the complex non‐linear behaviour of active shell structures including buckling, which is not captured by linear shell models. The numerical results also show the inherent capability of piezoelectric shell structures to actively induce large displacements through piezoelectric actuators, by jumping between multiple equilibrium states. Copyright © 2005 John Wiley & Sons, Ltd.     

On the dynamic behaviour of the Timoshenko beam finite elements

First, various finite element models of the Timoshenko beam theory for static analysis are reviewed, and a novel derivation of the 4 × 4 stiffness matrix (for the pure bending case) of the superconvergent finite element model for static problems is presented using two alternative approaches: (1) assumed-strain finite element model of the conventional Timoshenko beam theory, and (2) assumed-displacement finite element model of a modified Timoshenko beam theory. Next, dynamic versions of various finite element models are discussed. Numerical results for natural frequencies of simply supported beams are presented to evaluate various Timoshenko beam finite elements. It is found that the reduced integration element predicts the natural frequencies accurately, provided a sufficient number of elements is used. The research reported herein is supported by theOscar S. Wyatt Endowed Chair.     

A flexibility‐based finite element for linear analysis of arbitrarily curved arches

A new finite element for linear analysis of arbitrarily curved plane beams is developed. Based on a flexibility formulation the exact stiffness matrix and nodal force vector are obtained. In this way the effectiveness of the arch element depends solely on the capability to describe accurately the arch geometry. To this purpose a cubic parametric representation is adopted, which is local to the element and capable to model curved beams of any shape, while preserving C¹ continuity. In addition the proposed representation offers two inner parameters which can be used to locally adapt the arch shape. The resultant arch element is readily implementable into existing finite element codes and its effectiveness is validated through various benchmark problems, well established in the literature. Copyright © 2005 John Wiley & Sons, Ltd.     

Distortion‐immune nine‐node displacement‐based quadrilateral thick plate finite elements that satisfy constant‐bending patch test

We employ the linked interpolation concept to develop two higher‐order nine‐node quadrilateral plate finite elements with curved sides that pass the constant bending patch test for arbitrary node positions. The linked interpolation for the plate displacements is expanded with three bubble parameters to get polynomial completeness necessary to satisfy the patch test. In contrast to some other techniques, the elements developed in this way retain a symmetric stiffness matrix at a marginal computational expense at the element level. The new elements generated using this concept are tested on several examples with curved sides or some other kind of geometric distortion. Copyright © 2014 John Wiley & Sons, Ltd.     

Buckling Study of Thin‐walled Channel Beams with Double‐box Flanges in Pure Bending

Abstract: This paper provides the results of numerical and experimental investigations of buckling problems of cold‐formed, thin‐walled channel beams with double‐box flanges under pure bending. A local and global buckling analysis is realised numerically with the use of the finite strip method. A local buckling has been experimentally studied and also numerically with the use of the finite element method. Experimental tests of beams subjected to pure bending are conducted. The results of numerical and experimental investigations are presented and compared. A fundamental influence of double‐box flanges on the critical load is shown.