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.     

Director‐based beam finite elements relying on the geometrically exact beam theory formulated in skew coordinates

In the present work, a new director‐based finite element formulation for geometrically exact beams is proposed. The new beam finite element exhibits drastically improved numerical performance when compared with the previously developed director‐based formulations. This improvement is accomplished by adjusting the underlying variational beam formulation to the specific features of the director interpolation. In particular, the present approach does not rely on the assumption of an orthonormal director frame. The excellent performance of the new approach is illustrated with representative numerical examples. Copyright © 2013 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 geometrically non‐linear piezoelectric solid shell element based on a mixed multi‐field variational formulation

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

Frame‐indifferent beam finite elements based upon the geometrically exact beam theory

This paper describes a new beam finite element formulation based upon the geometrically exact beam theory. In contrast to many previously proposed beam finite element formulations the present discretization approach retains the frame‐indifference (or objectivity) of the underlying beam theory. The space interpolation of rotational degrees‐of‐freedom is circumvented by the introduction of a reparameterization of the weak form corresponding to the equations of motion of the geometrically exact beam theory. Copyright © 2002 John Wiley & Sons, Ltd.     

Theory and numerics for finite deformation fracture modelling using strong discontinuities

A general finite element approach for the modelling of fracture is presented for the geometrically non‐linear case. The kinematical representation is based on a strong discontinuity formulation in line with the concept of partition of unity for finite elements. Thus, the deformation map is defined in terms of one continuous and one discontinuous portion, considered as mutually independent, giving rise to a weak formulation of the equilibrium consisting of two coupled equations. In addition, two different fracture criteria are considered. Firstly, a principle stress criterion in terms of the material Mandel stress in conjunction with a material cohesive zone law, relating the cohesive Mandel traction to a material displacement ‘jump’ associated with the direct discontinuity. Secondly, a criterion of Griffith type is formulated in terms of the material‐crack‐driving force (MCDF) with the crack propagation direction determined by the direction of the force, corresponding to the direction of maximum energy release. Apart from the material modelling, the numerical treatment and aspects of computational implementation of the proposed approach is also thoroughly discussed and the paper is concluded with a few numerical examples illustrating the capabilities of the proposed approach and the connection between the two fracture criteria. Copyright © 2005 John Wiley & Sons, Ltd.     

Three‐dimensional numerical modelling by XFEM of spring‐layer imperfect curved interfaces with applications to linearly elastic composite materials

The spring‐layer interface model is widely used in describing some imperfect interfaces frequently involved in materials and structures. Typically, it is appropriate for modelling a thin soft interphase layer between two relatively stiff bulk media. According to the spring‐layer interface model, the displacement vector suffers a jump across an interface whereas the traction vector is continuous across the same interface and is, in the linear case, proportional to the displacement vector jump. In the present work, an efficient three‐dimensional numerical approach based on the extended finite element method is first proposed to model linear spring‐layer curved imperfect interfaces and then applied to predict the effective elastic moduli of composites in which such imperfect interfaces intervene. In particular, a rigorous derivation of the linear spring‐layer interface model is provided to clarify its domain of validity. The accuracy and convergence rate of the elaborated numerical approach are assessed via benchmark tests for which exact analytical solutions are available. The computated effective elastic moduli of composites are compared with the relevant analytical lower and upper bounds. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of approximations in analyses of beams of open thin‐walled cross‐section—part I: Flexural–torsional stability

In formulating a finite element model for the flexural–torsional stability and 3‐D non‐linear analyses of thin‐walled beams, a rotation matrix is usually used to obtain the non‐linear strain–displacement relationships. Because of the coupling between displacements, twist rotations and their derivatives, the components of the rotation matrix are both lengthy and complicated. To facilitate the formulation, approximations have been used to simplify the rotation matrix. A simplified small rotation matrix is often used in the formulation of finite element models for the flexural–torsional stability analysis of thin‐walled beams of open cross‐section. However, the approximations in the small rotation matrix may lead to the loss of some significant terms in the stability stiffness matrix. Without these terms, a finite element line model may predict the incorrect flexural–torsional buckling load of a beam. This paper investigates the effects of approximations in the elastic flexural–torsional stability analysis of thin‐walled beams, while a companion paper investigates the effects of approximations in the 3‐D non‐linear analysis. It is found that a finite element line model based on a small rotation matrix may predict incorrect elastic flexural–torsional buckling loads of beams. To perform a correct flexural–torsional stability analysis of thin‐walled beams, modification of the model is needed, or a finite element model based on a second‐order rotation matrix can be used. Copyright © 2001 John Wiley & Sons, Ltd.     

Finite rotation quadrilateral element for multi‐layer beams

The paper presents the finite rotations beam equations derived on use of the generalized Reissner hypothesis with a scalar parameter for the transverse extension. The beam strain and change of curvature measures are obtained from the right stretch strain, and the virtual work is given for Biot‐type stress and couple resultants. The strain energy for the first‐order isotropic elastic material is assumed in terms of the right stretch strain, and constitutive equations for the beam stress and couple resultants are derived. Two finite rotation elements are developed from the derived beam equations: a beam element with the transverse stretch and a quadrilateral element. First, the beam element with the uniformly under‐integrated tangent operator is developed. Next, the formula linking the middle‐line variables and the interface variables of the beam is introduced consistently with the generalized Reissner kinematics. Linearization of this formula is performed, and the derived tangent operator is used to convert the two‐node beam element to a four‐node quadrilateral. Both the finite elements have been tested on several numerical examples, some of highly non‐linear characteristics, and their accuracy is very good. It has been established that the quadrilateral element, which is intended for applications to multi‐layer beams, performs very well for high elemental aspect ratios, and can therefore be applied to modelling of very thin layers. Copyright © 1999 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.     

Effects of approximations in analyses of beams of open thin‐walled cross‐section—part II: 3‐D non‐linear behaviour

In a companion paper, the effects of approximations in the flexural‐torsional stability analysis of beams was studied, and it was shown that a second‐order rotation matrix was sufficiently accurate for a flexural‐torsional stability analysis. However, the second‐order rotation matrix is not necessarily accurate in formulating finite element model for a 3‐D non‐linear analysis of thin‐walled beams of open cross‐section. The approximations in the second‐order rotation matrix may introduce ‘self‐straining’ due to superimposed rigid‐body motions, which may lead to physically incorrect predictions of the 3‐D non‐linear behaviour of beams. In a 3‐D non‐linear elastic–plastic analysis, numerical integration over the cross‐section is usually used to check the yield criterion and to calculate the stress increments, the stress resultants, the elastic–plastic stress–strain matrix and the tangent modulus matrix. A scheme of the arrangement of sampling points over the cross‐section that is not consistent with the strain distributions may lead to incorrect predictions of the 3‐D non‐linear elastic–plastic behaviour of beams. This paper investigates the effects of approximations on the 3‐D non‐linear analysis of beams. It is found that a finite element model for 3‐D non‐linear analysis based on the second‐order rotation matrix leads to over‐stiff predictions of the flexural‐torsional buckling and postbuckling response and to an overestimate of the maximum load‐carrying capacities of beams in some cases. To perform a correct 3‐D non‐linear analysis of beams, an accurate model of the rotations must be used. A scheme of the arrangement of sampling points over the cross‐section that is consistent with both the longitudinal normal and shear strain distributions is needed to predict the correct 3‐D non‐linear elastic–plastic behaviour of beams. Copyright © 2001 John Wiley & Sons, Ltd.     

Buckling of core wall structures coupled with connecting beams

A rational and accurate numerical analysis is presented for the buckling of core walls coupled with connecting beams by using spline finite member element method (SFMEM) in which the effects of torsion, warping and, especially, the shearing strains in the middle surface of the walls are taken into account. The core wall structure is discretized and modelled as an equivalent thin‐walled closed section, while the connecting beams between openings are replacement by an equivalent shear diaphragm of thickness. The numerical method combines the advantages of B₃‐spline, the finite member element method and Vlasov's thin‐walled bar theory. The simplicity and accuracy of the proposed scheme are demonstrated by a numerical example. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical integration of the stiff dynamics of geometrically exact shells: an energy‐dissipative momentum‐conserving scheme

This paper presents a new family of time‐stepping algorithms for the integration of the dynamics of non‐linear shells. We consider the geometrically exact shell theory involving an inextensible director field (the so‐called five‐parameter shell model). The main characteristic of this model is the presence of the group of finite rotations in the configuration manifold describing the deformation of the solid. In this context, we develop time‐stepping algorithms whose discrete solutions exhibit the same conservation laws of linear and angular momenta as the underlying physical system, and allow the introduction of a controllable non‐negative energy dissipation to handle the high numerical stiffness characteristic of these problems. A series of algorithmic parameters for the different components of the deformation of the shell (i.e. membrane, bending and transverse shear) fully control this numerical dissipation, recovering existing energy‐momentum schemes as a particular choice of these algorithmic parameters. We present rigorous proofs of the numerical properties of the resulting algorithms in the full non‐linear range. Furthermore, it is argued that the numerical dissipation is introduced in the high‐frequency range by considering the proposed algorithm in the context of a linear problem. The finite element implementation of the resulting methods is described in detail as well as considered in the final arguments proving the aforementioned conservation/dissipation properties. We present several representative numerical simulations illustrating the performance of the newly proposed methods. The robustness gained over existing methods in these stiff problems is confirmed in particular. Copyright © 2002 John Wiley & Sons, Ltd.     

An extension of Hill's three‐component model to include different fibre types in finite element modelling of muscle

Most of the proposed versions of the Hill's model use a sliding‐element theory, considering a single sarcomere. However, a muscle represents a collection of different fibre types with a large range in contractile properties among them. An extension of Hill's three‐component model is proposed here to take into account different fibre types. We present a model consisting of a number of sarcomeras of different types coupled in parallel with the connective tissue. Each sarcomere is modelled by one non‐linear elastic element connected in series with one non‐linear contractile element. Using the finite element method, in an incremental‐iterative scheme of calculating equilibrium configurations of a muscle, the key step is the determination of stresses corresponding to strain increments. The stress calculation procedure for the extended Hill's model is reduced to the solution of a number of independent non‐linear equations with respect to the stretch increments of the serial elastic elements in each sarcomere. Since the distribution of the specific fibre type is non‐uniform over the muscle volume, we have material heterogeneity which we modelled by using the so‐called ‘Generalized Isoparametric Element Formulation’ for functionally graded materials (FGMs). The proposed computational scheme is built in our FE package PAK, so that muscles of complex three‐dimensional shapes can be modelled. In numerical examples, we illustrate the main characteristics of the developed numerical model and some possibilities of realistic modelling of muscle functioning. Copyright © 2006 John Wiley & Sons, Ltd.     

On evaluation of shape sensitivities of non‐linear critical loads

The present paper focuses on the evaluation of the shape sensitivities of the limit and bifurcation loads of geometrically non‐linear structures. The analytical approach is applied for isoparametric elements, leading to exact results for a given mesh. Since this approach is difficult to apply to other element types, the semi‐analytical method has been widely used for shape sensitivity computation. This method combines ease of implementation with computational efficiency, but presents severe accuracy problems. Thus, a general procedure to improve the semi‐analytical sensitivities of the non‐linear critical loads is presented. The numerical examples show that this procedure leads to sensitivities with sufficient accuracy for shape optimization applications. Copyright © 2002 John Wiley & Sons, Ltd.     

The enriched space–time finite element method (EST) for simultaneous solution of fluid–structure interaction

The paper introduces a weighted residual‐based approach for the numerical investigation of the interaction of fluid flow and thin flexible structures. The presented method enables one to treat strongly coupled systems involving large structural motion and deformation of multiple‐flow‐immersed solid objects. The fluid flow is described by the incompressible Navier–Stokes equations. The current configuration of the thin structure of linear elastic material with non‐linear kinematics is mapped to the flow using the zero iso‐contour of an updated level set function. The formulation of fluid, structure and coupling conditions uniformly uses velocities as unknowns. The integration of the weak form is performed on a space–time finite element discretization of the domain. Interfacial constraints of the multi‐field problem are ensured by distributed Lagrange multipliers. The proposed formulation and discretization techniques lead to a monolithic algebraic system, well suited for strongly coupled fluid–structure systems. Embedding a thin structure into a flow results in non‐smooth fields for the fluid. Based on the concept of the extended finite element method, the space–time approximations of fluid pressure and velocity are properly enriched to capture weakly and strongly discontinuous solutions. This leads to the present enriched space–time (EST) method. Numerical examples of fluid–structure interaction show the eligibility of the developed numerical approach in order to describe the behavior of such coupled systems. The test cases demonstrate the application of the proposed technique to problems where mesh moving strategies often fail. Copyright © 2007 John Wiley & Sons, Ltd.     

A dissipation‐based arc‐length method for robust simulation of brittle and ductile failure

A robust method to trace the equilibrium path in non‐linear solid mechanics problems is proposed. A general arc‐length constraint based on the energy release rate is developed. Constraints have been derived for the cases of geometrically linear damage, geometrically linear plasticity and geometrically non‐linear damage. All three constraints can efficiently be applied in a finite element context. Numerical simulations demonstrate that the proposed framework gives robust results for these cases. Applicability of the proposed framework to other types of constitutive and/or kinematic behaviour is predicted. Copyright © 2008 John Wiley & Sons, Ltd.     

Aeroelastic forces and dynamic response of long‐span bridges

In this paper a time domain approach for predicting the non‐linear dynamic response of long‐span bridges is presented. In particular the method that leads to the formulation of aeroelastic and buffeting forces in the time domain is illustrated in detail, where a recursive algorithm for the memory term's integration is properly developed. Moreover in such an approach the forces' expressions, usually formulated according to quasi‐static theory, have been substituted by expressions including the frequency‐dependent characteristics. Such expressions of aeroelastic and buffeting forces are made explicit in the time domain by means of the convolution integral that involves the impulse functions and the structural motion or the fluctuating velocities. A finite element model (FEM) has been developed within the framework of geometrically non linear analysis, by using 3‐d degenerated finite element. The proposed procedure can be used to analyze both the flutter instability phenomenon and buffeting response. Moreover, working in the geometrically non‐linearity range, it verifies the possibility of strongly flexible structures of actively resisting the wind loading. Copyright © 2004 John Wiley & Sons, Ltd.     

Masonry infilled reinforced concrete frames under horizontal loading / Stahlbetonrahmen mit Ausfachungen aus Mauerwerk unter horizontalen Belastungen

The behaviour of infilled reinforced concrete frames under horizontal load has been widely investigated, both experimentally and numerically. Since experimental tests represent large investments, numerical simulations offer an efficient approach for a more comprehensive analysis. When RC frames with masonry infill walls are subjected to horizontal loading, their behaviour is highly non‐linear after a certain limit, which makes their analysis quite difficult. The non‐linear behaviour results from the complex inelastic material properties of the concrete, infill wall and conditions at the wall‐frame interface. In order to investigate this non‐linear behaviour in detail, a finite element model using a micro modelling approach is developed, which is able to predict the complex non‐linear behaviour resulting from the different materials and their interaction. Concrete and bricks are represented by a non‐linear material model, while each reinforcement bar is represented as an individual part installed in the concrete part and behaving elasto‐plastically. Each brick is modelled individually and connected taking into account the non‐linearity of a brick mortar interface. The same approach is followed using two finite element software packages and the results are compared with the experimental results. The numerical models show a good agreement with the experiments in predicting the overall behaviour, but also very good matching for strength capacity and drift. The results emphasize the quality and the valuable contribution of the numerical models for use in parametric studies, which are needed for the derivation of design recommendations for infilled frame structures.