A least-squares variational method for full-field reconstruction of elastic deformations in shear-deformable plates and shells

A variational principle is formulated for the inverse problem of full-field reconstruction of three-dimensional plate/shell deformations from experimentally measured surface strains. The formulation is based upon the minimization of a least-squares functional that uses the complete set of strain measures consistent with linear, first-order shear-deformation theory. The formulation, which accommodates for transverse shear-deformation, is applicable for the analysis of thin and moderately thick plate and shell structures. The main benefit of the variational principle is that it is well-suited for C⁰-continuous displacement finite element discretizations, thus enabling the development of robust algorithms for application to complex civil and aeronautical structures. The methodology is especially aimed at the next generation of aerospace vehicles for use in real-time structural health monitoring systems.     

Mixed least-squares finite element model for the static analysis of laminated composite plates

This paper presents a mixed finite element model for the static analysis of laminated composite plates. The formulation is based on the least-squares variational principle, which is an alternative approach to the mixed weak form finite element models. The mixed least-squares finite element model considers the first-order shear deformation theory with generalized displacements and stress resultants as independent variables. Specifically, the mixed model is developed using equal-order C⁰ Lagrange interpolation functions of high p-levels along with full integration. This mixed least-squares-based discrete model yields a symmetric and positive-definite system of algebraic equations. The predictive capability of the proposed model is demonstrated by numerical examples of the static analysis of four laminated composite plates, with different boundary conditions and various side-to-thickness ratios. Particularly, the mixed least-squares model with high-order interpolation functions is shown to be insensitive to shear-locking.     

Variational formulation of linear problems with nonhomogeneous boundary conditions and internal discontinuities

A variational formulation applicable to linear operators with nonhomogeneous boundary conditions and jump discontinuities is presented. For the formulation to be applicable, the boundary condition and discontinuities have to be consistent with the operator governing the field problem. The problem is set up in a space of suitable continuous bilinear mapping. Thus, operators on inner product spaces, convolution spaces and energy spaces are included as specializations. The basic construction can be used to generate dual-complementary variational principles. Implementation is illustrated by examples. The role of boundary terms in finite element discretizations based on interpolants of limited smoothness is discussed.     

Insight into finite element shell discretizations by use of the “basic shell mathematical model”

The objective of this paper is to gain insight into finite element discretizations of shells using the basic shell mathematical model and, in particular, regarding the sources of “locking”. We briefly review the “basic shell mathematical model” and present a formulation of shell finite elements based on this model. These shell finite elements are equivalent to the widely-used continuum mechanics based shell finite elements. We consider a free hyperboloid shell problem, which is known to be difficult to solve accurately. Using a fine mesh of MITC9 elements based on the basic shell mathematical model, a detailed analysis is performed giving the distributions of all strain terms. A similar analysis using the MITC6 shell element shows why this element locks when the shell thickness is very small.     

Discontinuous Petrov–Galerkin method with optimal test functions for thin-body problems in solid mechanics

We study the applicability of the discontinuous Petrov–Galerkin (DPG) variational framework for thin-body problems in structural mechanics. Our numerical approach is based on discontinuous piecewise polynomial finite element spaces for the trial functions and approximate, local computation of the corresponding ‘optimal’ test functions. In the Timoshenko beam problem, the proposed method is shown to provide the best approximation in an energy-type norm which is equivalent to the L₂-norm for all the unknowns, uniformly with respect to the thickness parameter. The same formulation remains valid also for the asymptotic Euler–Bernoulli solution. As another one-dimensional model problem we consider the modelling of the so called basic edge effect in shell deformations. In particular, we derive a special norm for the test space which leads to a robust method in terms of the shell thickness. Finally, we demonstrate how a posteriori error estimator arising directly from the discontinuous variational framework can be utilized to generate an optimal hp-mesh for resolving the boundary layer.     

A stabilized MITC element for accurate wave response in Reissner–Mindlin plates

Residual based finite element methods are developed for accurate time-harmonic wave response of the Reissner–Mindlin plate model. The methods are obtained by appending a generalized least-squares term to the mixed variational form for the finite element approximation. Through judicious selection of the design parameters inherent in the least-squares modification, this formulation provides a consistent and general framework for enhancing the wave accuracy of mixed plate elements. In this paper, the mixed interpolation technique of the well-established MITC4 element is used to develop a new mixed least-squares (MLS4) four-node quadrilateral plate element with improved wave accuracy. Complex wave number dispersion analysis is used to design optimal mesh parameters, which for a given wave angle, match both propagating and evanescent analytical wave numbers for Reissner–Mindlin plates. Numerical results demonstrates the significantly improved accuracy of the new MLS4 plate element compared to the underlying MITC4 element.     

Multiscale temporal integration

We study the application of the variational multiscale method to the problem of temporal integration, with the final goal of designing integration schemes that go beyond the classical notion of upwinding.We develop a formulation based on a mixed hybrid finite element method, where the fine scale mode problems automatically decouple at the element level without the need to resort to a localization assumption. We give general orthogonality conditions for the trial and test spaces that allow to construct hierarchical p methods.We test some simple ideas for the modeling of the unresolved scales. The resulting algorithms are analyzed using classical analytical measures.     

A penalty duality method for the budiansky-sanders shell model

Starting from the Budiansky-Sanders shell model, we introduce a perturbed formulation based on a penalty-duality argument for which only C ° finite elements are required.An error estimate and numerical procedure for solving the finite element model are given.     

A thin cylindrical shell finite element based on a mixed formulation

A specialized functional for thin cylindrical shells derived from the Washizu-Hu variational principle using considerations of relaxed continuity requirements is presented. A mixed formulation for a cylindrical thin shell finite element is developed from this functional. The assumed fields for displacements and stress resultants are bilinear functions in the longitudinal and circunferential directions.The agreement between the present results and those obtained in previous formulations is good if the comparison is based on the precision related to the number of variables involved in the problem.     

An 8-node assumed stress hybrid element for analysis of shells

In this paper an 8-node quadrilateral assumed-stress hybrid shell element is presented. The formulation is based on Hellinger–Reissner variational principle. The element is developed by flat shell approach by combining a membrane element with a Mindlin plate element. The proposed element has six degrees of freedom per node and permits an easy connection to other types of finite elements. Numerical examples are presented to show that the validity and efficiency of the present element for static and free vibration analysis.     

Nonlinear preconditioned conjugate gradient and least-squares finite elements

A least-squares variational formulation for first-order systems is extended to a class of nonlinear problems. A finite element analysis and associated preconditioners for iterative solution are developed. The preconditioned system is not sensitive to problem (mesh) size as is demonstrated in numerical comparison studies.     

Solid-to-shell transition elements for the computation of interlaminar stresses

This paper presents an accurate and practical technique for coupling shell element models to three-dimensional continuum finite element models. The compatibility between these two types of formulations is enforced by degenerating a continuum element through kinematic constraints compatible with shell deformations. Two formulations of two-dimensional/three-dimensional transition elements are presented. The first and simplest formulation is based on the Mindlin-Reissner plate assumptions, and is found to perform well in a variety of problems involving the analysis of geometrically linear/non-linear laminated structures. The second formulation is based on a higher-order shell theory that allows stretching in the through-the-thickness direction. This additional freedom virtually eliminates the interlaminar normal stress boundary layer that can form in lower-order transition elements. Finally, the coupling of two-dimensional to three-dimensional subdomains is enriched with the use of an interface element, which can be used in conjunction with either transition formulation. The interface element improves the efficiency of the solid-to-shell transition modeling scheme by allowing the independent selection of optimal mesh sizes in the shell and the three-dimensional regions of the model.     

Stability of cylindrical shell panels subjected to follower forces based on a mixed finite element formulation

A mixed finite element formulation is developed from a weak variational priniciple. This formulation is applied to stability analysis of cylindrical shell structures subjected to follower loading. Bilinear trial functions are used for all field variables. The rectangular curved elements presented here satisfy the continuity requirements for the field variables at the element interface. Two examples of a cantilevered cylindrical shell panel under different kinds of loading are solved.     

A robust non-linear solid shell element based on a mixed variational formulation

The paper is concerned with a geometrically non-linear solid shell finite element formulation, which is based on the Hu-Washizu variational principle. For the approximation of the independent displacement, stress and strain fields, the strain field is additively decomposed into two parts. Due to the fact that one part of the strain field is interpolated in the same manner as proposed by the enhanced assumed strain (EAS) method, it is denoted as EAS field. The other strain field is approximated with the same interpolation functions as the stress field. In contrast to the EAS concept the approximation spaces of the stresses and the enhanced assumed strains are not orthogonal. Consequently the stress field is not eliminated from the finite element equations. For the displacements tri-linear shape functions are considered. Shear locking and curvature thickness locking are treated using assumed natural strain interpolations. A static condensation leads to a simple low order hexahedral solid shell element. Numerical tests show that the present model is very robust and allows larger load steps than an EAS solid shell element.     

Gauss–Newton Multilevel Methods for Least-Squares Finite Element Computations of Variably Saturated Subsurface Flow

We apply the least-squares mixed finite element framework to the nonlinear elliptic problems arising in each time-step of an implicit Euler discretization for variably saturated flow. This approach allows the combination of standard piecewise linear H ¹-conforming finite elements for the hydraulic potential with the H(div)-conforming Raviart–Thomas spaces for the flux. It also provides an a posteriori error estimator which may be used in an adaptive mesh refinement strategy. The resulting nonlinear algebraic least-squares problems are solved by an inexact Gauss–Newton method using a stopping criterion for the inner iteration which is based on the change of the linearized least-squares functional relative to the nonlinear least-squares functional. The inner iteration is carried out using an adaptive multilevel method with a block Gauss–Seidel smoothing iteration. For a realistic water table recharge problem, the results of computational experiments are presented. Received January 4, 1999; revised July 19, 1999     

Transition finite elements with temperature and temperature gradients as primary variables for axisymmetric heat conduction

This paper presents finite element formulation for a special class of elements referred to as “transition finite elements” for axisymmetric heat conduction. The transition elements are necessary in applications requiring the use of both axisymmetric solid elements and axisymmetric shell elements. The elements permit transition from the solid portion of the structure to the shell portion of the structure. A novel feature of the formulation presented here is that nodal temperatures as well as nodal temperature gradients are retained as primary variables. The weak formulation of the Fourier heat conduction equation is constructed in the cylindrical co-ordinate system (r, z). The element geometry is defined in terms of the co-ordinates of the nodes as well as the nodal point normals for the nodes lying on the middle surface of the element. The element temperature field is approximated in terms of element approximation functions, nodal temperatures and the nodal temperature gradients. The properties of the transition elements are then derived using the weak formulation and the element temperature approximation. The formulation presented here permits linear temperature distribution through the element thickness. Convective boundaries as well as distributed heat flux is permitted on all four faces of the element. Furthermore, the element formulation also permits distributed heat flux and orthotropic material behaviour. Numerical examples are presented, first to illustrate the accuracy of the formulation and second to demonstrate its usefulness in practical applications. Numerical results are also compared with the theoretical solutions.     

Discrete energy method for the analysis of cylindrical shells

This paper presents an investigation into the use of the closely associated finite difference technique for the analysis of shell structures as a feasible alternative to the finite element method. The method discretises the total energy of the structure into energy due to extension and bending and that due to shear and twisting, contributed by two separate sets of rectangular elements formed by a suitable finite difference network. The derivatives in the corresponding energy expressions are replaced by their finite difference forms and the nodal displacements then constitute the undetermined parameters in the variational formulation. The formulation is also extended to a cylindrical shell element of rectangular planform. The results obtained by DEM are compared with existing results and they show excellent agreement.     

Hybrid strain based three node flat triangular shell elements—II. Numerical investigation of nonlinear problems

In a companion paper [M. L. Liu and C. W. S. To, Comput. Struct. 54, 1031–1056 (1995)] theories and incremental formulation of nonlinear shell structures discretized by the finite element method are discussed. The updated Lagrangian formulation and the incremental Hellinger-Reissner variational principle are adopted. The independently assumed fields employed are the incremental displacements and incremental strains. Based on the theory and incremental formulation explicit element stiffness and mass matrices of three node flat triangular shell finite elements are derived. In the present paper the derived element matrices are applied to nine examples. The latter include static and dynamic response analysis of shell structures with geometrical, material, and geometrical and material nonlinearities. The formulation adopted and element matrices derived are found to be accurate, flexible and applicable to various types of shell structures with geometrical and material nonlinearities.     

Analysis of multibody systems with flexible plates using variational graph-theoretic methods

This work treats the problem of modelling multibody systems with structural flexibility. By combining linear graph theory with the principle of virtual work and finite elements, a dynamic formulation is obtained that extends graph-theoretic (GT) modelling methods to the analysis of thin flexible plates for multibody systems. The system is represented by a linear graph, in which nodes represent reference frames on flexible plates, and edges represent components that connect these frames. To generate the equations of motion with elastic deformations, the flexible plates are discretized using a triangular thin shell finite element based on the discrete Kirchhoff criterion and can be used to discretize bidirectional bodies such as satellite panels, flatbed trailers, and mechanisms with plates. Three flexible systems with plates are analyzed to illustrate the performance of this new variational graph-theoretic formulation and its ability to generate directly a set of motion equations for flexible multibody systems (FMS) without additional user input.