A SPACE-TIME COUPLED p-VERSION LEAST SQUARES FINITE ELEMENT FORMULATION FOR UNSTEADY TWO-DIMENSIONAL NAVIER–STOKES EQUATIONS

This paper presents a p-version least squares finite element formulation for two-dimensional unsteady fluid flow described by Navier–Stokes equations where the effects of space and time are coupled. The dimensionless form of the Navier–Stokes equations are first cast into a set of first-order differential equations by introducing auxiliary variables. This permits the use of C⁰ element approximation. The element properties are derived by utilizing the p-version approximation functions in both space and time and then minimizing the error functional given by the space–time integral of the sum of squares of the errors resulting from the set of first-order differential equations. This results in a true space–time coupled least squares minimization procedure. The application of least squares minimization to the set of coupled first-order partial differential equations results in finding a solution vector {δ} which makes gradient of error functional with respect to {δ} a null vector. This is accomplished by using Newton's method with a line search. A time marching procedure is developed in which the solution for the current time step provides the initial conditions for the next time step. Equilibrium iterations are carried out for each time step until the error functional and each component of the gradient of the error functional with respect to nodal degrees of freedom are below a certain prespecified tolerance. The space–time coupled p-version approximation functions provide the ability to control truncation error which, in turn, permits very large time steps. What literally requires hundreds of time steps in uncoupled conventional time marching procedures can be accomplished in a single time step using the present space–time coupled approach. The generality, success and superiority of the present formulation procedure is demonstrated by presenting specific numerical examples for transient couette flow and transient lid driven cavity. The results are compared with the analytical solutions and those reported in the literature. The formulation presented here is ideally suited for space–time adaptive procedures. The element error functional values provide a mechanism for adaptive h, p or hp refinements.     

Formulation and implementation of stress‐driven and/or strain‐driven computational homogenization for finite strain

In this paper, we present a homogenization approach that can be used in the geometrically nonlinear regime for stress‐driven and strain‐driven homogenization and even a combination of both. Special attention is paid to the straightforward implementation in combination with the finite‐element method. The formulation follows directly from the principle of virtual work, the periodic boundary conditions, and the Hill–Mandel principle of macro‐homogeneity. The periodic boundary conditions are implemented using the Lagrange multiplier method to link macroscopic strain to the boundary displacements of the computational model of a representative volume element. We include the macroscopic strain as a set of additional degrees of freedom in the formulation. Via the Lagrange multipliers, the macroscopic stress naturally arises as the associated ‘forces’ that are conjugate to the macroscopic strain ‘displacements’. In contrast to most homogenization schemes, the second Piola–Kirchhoff stress and Green–Lagrange strain have been chosen for the macroscopic stress and strain measures in this formulation. The usage of other stress and strain measures such as the first Piola–Kirchhoff stress and the deformation gradient is discussed in the Appendix. Copyright © 2015 John Wiley & Sons, Ltd.     

A finite‐strain quadrilateral shell element based on discrete Kirchhoff–Love constraints

This paper improves the 16 degrees‐of‐freedom quadrilateral shell element based on pointwise Kirchhoff–Love constraints and introduces a consistent large strain formulation for this element. The model is based on classical shell kinematics combined with continuum constitutive laws. The resulting element is valid for large rotations and displacements. The degrees‐of‐freedom are the displacements at the corner nodes and one rotation at each mid‐side node. The formulation is free of enhancements, it is almost fully integrated and is found to be immune to locking or unstable modes. The patch test is satisfied. In addition, the formulation is simple and amenable to efficient incorporation in large‐scale codes as no internal degrees‐of‐freedom are employed, and the overall calculations are very efficient. Results are presented for linear and non‐linear problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Hybrid equilibrium element with interelement interface for the analysis of delamination and crack propagation problems

This article proposes a formulation for the analysis of delamination and fracture propagation problems at the interelement interface, with perfect adhesion at the pre-failure condition and with linear softening at the post-failure regime. The proposed formulation is based on the hybrid equilibrium element (HEE) model, with stress fields which strongly verify the homogeneous equilibrium equations and interelement equilibrium equations. The HEE can easily model high-order stress fields and can implicitly model the initially rigid behavior of an extrinsic interface at the element sides. The interface model is defined as a function of the same degrees of freedom of the HEE (generalized stresses) and the pre- and post-failure behavior of the interface can be modeled without any additional degree of freedom. The proposed formulation is developed for a nine-node triangular HEE with quadratic stress fields. This article also proposes an extrinsic cohesive model, rigorously developed in the damage mechanics framework.     

Structural damage identification: influence of model incompleteness and errors

This paper presents a study of the influence of model incompleteness and errors in a structural damage identification technique based on modal data, leading to a set of linear equations. To deal with the incompleteness three methods are used: (1) dynamic and (2) static expansions of measured degrees of freedom and (3) replacement of unmeasured degrees of freedom by their undamaged counterparts. These methods are applied on damage identification of a laminated rectangular plate. The numerical tests show that the best accuracy is obtained with the dynamic expansion, followed by the static expansion and finally, the replacement of unmeasured degrees of freedom. It is also demonstrated that for small damage the errors are the main influence, whereas for large damage the incompleteness becomes the most important factor in the results. The identification of small and large damages is performed by deleting non-reliable equations, i.e., equations that contain large errors.     

Thin plate elements with relaxed Kirchhoff constraints

An analysis of thin plates by C⁰ elements is presented. The fundamental part of it is the formulation of a set of constraints that relate the rotational degrees of freedom to the translational degrees of freedom and which are weaker than the usual Kirchhoff constraints. A penalty term associated with these constraints is then introduced instead of the usual expression for the shear strain energy. It is shown that the approach is effective for both quadrilateral and triangular elements.     

Phase‐field modeling of ductile fracture at finite strains: A robust variational‐based numerical implementation of a gradient‐extended theory by micromorphic regularization

This work provides a robust variational‐based numerical implementation of a phase field model of ductile fracture in elastic–plastic solids undergoing large strains. This covers a computationally efficient micromorphic regularization of the coupled gradient plasticity‐damage formulation. The phase field approach regularizes sharp crack surfaces within a pure continuum setting by a specific gradient damage modeling with geometric features rooted in fracture mechanics. It has proven immensely successful with regard to the analysis of complex crack topologies without the need for fracture‐specific computational structures such as finite element design of crack discontinuities or intricate crack‐tracking algorithms. The proposed gradient‐extended plasticity‐damage formulation includes two independent length scales that regularize both the plastic response as well as the crack discontinuities. This ensures that the damage zones of ductile fracture are inside of plastic zones or vice versa and guarantees on the computational side a mesh objectivity in post‐critical ranges. The proposed setting is rooted in a canonical variational principle. The coupling of gradient plasticity to gradient damage is realized by a constitutive work density function that includes the stored elastic energy and the dissipated work due to plasticity and fracture. The latter represents a coupled resistance to plasticity and damage, depending on the gradient‐extended internal variables that enter plastic yield functions and fracture threshold functions. With this viewpoint on the generalized internal variables at hand, the thermodynamic formulation is outlined for gradient‐extended dissipative solids with generalized internal variables that are passive in nature. It is specified for a conceptual model of von Mises‐type elasto‐plasticity at finite strains coupled with fracture. The canonical theory proposed is shown to be governed by a rate‐type minimization principle, which fully determines the coupled multi‐field evolution problem. This is exploited on the numerical side by a fully symmetric monolithic finite element implementation. An important aspect of this work is the regularization towards a micromorphic gradient plasticity‐damage setting by taking into account additional internal variable fields linked to the original ones by penalty terms. This enhances the robustness of the finite element implementation, in particular, on the side of gradient plasticity. The performance of the formulation is demonstrated by means of some representative examples. Copyright © 2016 John Wiley & Sons, Ltd.     

A space–time coupled p-version least-squares finite element formulation for unsteady fluid dynamics problems

This paper presents a p-version least-squares finite element formulation for unsteady fluid dynamics problems where the effects of space and time are coupled. The dimensionless form of the differential equations describing the problem are first cast into a set of first-order differential equations by introducing auxiliary variables. This permits the use of C° element approximation. The element properties are derived by utilizing p-version approximation functions in both space and time and then minimizing the error functional given by the space–time integral of the sum of squares of the errors resulting from the set of first-order differential equations. This results in a true space–time coupled least-squares minimization procedure. A time marching procedure is developed in which the solution for the current time step provides the initial conditions for the next time step. The space–time coupled p-version approximation functions provide the ability to control truncation error which, in turn, permits very large time steps. What literally requires hundreds of time steps in uncoupled conventional time marching procedures can be accomplished in a single time step using the present space–time coupled approach. For non-linear problems the non-linear algebraic equations resulting from the least-squares process are solved using Newton's method with a line search. This procedure results in a symmetric Hessian matrix. Equilibrium iterations are carried out for each time step until the error functional and each component of the gradient of the error functional with respect to nodal degrees of freedom are below a certain prespecified tolerance. The generality, success and superiority of the present formulation procedure is demonstrated by presenting specific formulations and examples for the advection–diffusion and Burgers equations. The results are compared with the analytical solutions and those reported in the literature. The formulation presented here is ideally suited for space–time adaptive procedures. The element error functional values provide a mechanism for adaptive h, p or hp refinements. The work presented in this paper provides the basis for the extension of the space–time coupled least-squares minimization concept to two- and three-dimensional unsteady fluid flow.     

Assumed strain nodally integrated hexahedral finite element formulation for elastoplastic applications

In this work, a linear hexahedral element based on an assumed strain finite element technique is presented for the solution of plasticity problems. The element stems from the Nodally Integrated Continuum Element (NICE) formulation and its extensions. Assumed gradient operators are derived via nodal integration from the kinematic‐weighted residual; the degrees of freedom are only the displacements at the nodes. The adopted constitutive model is the classical associative von Mises plasticity model with isotropic and kinematic hardening; in particular, a double‐step midpoint integration algorithm is adopted for the integration and solution of the relevant nonlinear evolution equations. Efficiency of the proposed method is assessed through simple benchmark problems and comparison with reference solutions. Copyright © 2014 John Wiley & Sons, Ltd.     

Dynamic analysis of rail vehicle axle

In this paper, in order to obtain the dynamic forces on the passenger coach axle, a full rail vehicle model with 19-dof (degrees of freedom) has been considered. For a specific example, the variations of these dynamic forces with velocity of the passenger coach, suspension characteristics and way conditions have been examined. Dynamic forces found in the resonance regions at the range of 2–5 m/s (7.2–18 km/h) has been discussed. Theoretical results obtained for the dynamic forces have been successfully compared with the experimental results of German Railways (Deutsche Bahn-DB).     

Influence of porosity on the flexural and vibration response of gradient plate using nonpolynomial higher-order shear and normal deformation theory

The present article deals with flexural and vibration response of functionally graded plates with porosity. The basic formulation is based on the recently developed non-polynomial higher-order shear and normal deformation theory by the authors’. The present theory contains only four unknowns and also accommodate the thickness stretching effect. The effective material properties at each point are determined by two micromechanics models (Voigt and Mori–Tanaka scheme). The governing equations for FGM plates are derived using variational approach. Results have been obtained by employing a C⁰ continuous isoparametric Lagrangian finite element with eight degrees of freedom per node. Convergence and comparative study with the reported results in the literature, confirm the accuracy and efficiency of the present model and finite element formulation. The influence of the porosity, various boundary conditions, geometrical configuration and micromechanics models on the flexural and vibration behavior of FGM plates is examined.     

Proper orthogonal decomposition and modal analysis for acceleration of transient FEM thermal analysis

A method of reducing the number of degrees of freedom and the overall computing times in finite element method (FEM) has been devised. The technique is valid for linear problems and arbitrary temporal variation of boundary conditions. At the first stage of the method standard FEM time stepping procedure is invoked. The temperature fields obtained for the first few time steps undergo statistical analysis yielding an optimal set of globally defined trial and weighting functions for the Galerkin solution of the problem at hand. Simple matrix manipulations applied to the original FEM system produce a set of ordinary differential equations of a dimensionality greatly reduced when compared with the original FEM formulation. Using the concept of modal analysis the set is then solved analytically. Treatment of non‐homogeneous initial conditions, time‐dependent boundary conditions and controlling the error introduced by the reduction of the degrees of freedom are discussed. Several numerical examples are included for validation of the approach. Copyright © 2004 John Wiley & Sons, Ltd.     

Non‐linear exact geometry 12‐node solid‐shell element with three translational degrees of freedom per node

This paper presents the finite rotation exact geometry (EG) 12‐node solid‐shell element with 36 displacement degrees of freedom. The term ‘EG’ reflects the fact that coefficients of the first and second fundamental forms of the reference surface and Christoffel symbols are taken exactly at each element node. The finite element formulation developed is based on the 9‐parameter shell model by employing a new concept of sampling surfaces (S‐surfaces) inside the shell body. We introduce three S‐surfaces, namely, bottom, middle and top, and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows one to represent the finite rotation higher order EG solid‐shell element formulation in a very compact form and to derive the strain–displacement relationships, which are objective, that is, invariant under arbitrarily large rigid‐body shell motions in convected curvilinear coordinates. The tangent stiffness matrix is evaluated by using 3D analytical integration and the explicit presentation of this matrix is given. The latter is unusual for the non‐linear EG shell element formulation. Copyright © 2011 John Wiley & Sons, Ltd.     

Non-local yield limit degradation

Presented is a new type of a non-local continuum model which avoids problems of convergence at mesh refinement and spurious mesh sensitivity in a softening continuum characterized by degradation of the yield limit. The key idea, which has recently been proposed in a general context and has already been applied to softening damage due to stiffness degradation, is to apply the non-local concept only to those parameters which cause the degradation while keeping the definition of the strains local. Compared to the previously advanced fully non-local continuum formulation, the new approach has the advantage that the stresses are subjected to the standard differential equations of equilibrium and standard boundary or interface conditions. The new formulation exhibits no zero-energy periodic modes, imbrication of finite elements is unnecessary and finite elements with standard continuity requirements are sufficient. Two-dimensional finite element solutions with up to 3248 degrees of freedom are presented to document convergence and efficacy. The formulation is applied to tunnel excavation in a soil stabilized by cement grouting, with the objective of preventing cave-in (burst) of the tunnel sides due to compression softening.     

A hybrid element formulation for the three‐dimensional penalty finite element analysis of incompressible fluids

This work presents a hybrid element formulation for the three‐dimensional penalty finite element analysis of incompressible Newtonian fluids. The formulation is based on a mixed variational statement in which velocity and stresses are treated as independent field variables. The main advantage of this formulation is that it bypasses the use of ad hoc techniques such as selective reduced integration that are commonly used in penalty‐based finite element formulations, and directly yields high accuracy for the velocity and stress fields without the need to carry out smoothing. In addition, since the stress degrees of freedom are condensed out at an element level, the cost of solving for the global degrees of freedom is the same as in a standard penalty finite element method, although the gain in accuracy for both the velocity and stress (including the pressure) fields is quite significant. Copyright © 2007 John Wiley & Sons, Ltd.     

A penalty finite element approach for couple stress elasticity

There is mounting evidence for size dependent elastic deformation at micron and submicron length scales. Material formulations incorporating higher order gradients in displacements have been successful in modeling such size dependent deformation behavior. A couple stress theory without micro-rotation is considered here as micro-rotations increase complexity and necessitate parameters that are difficult to determine. Higher order gradient theories require continuity in both displacements and their derivatives and direct approaches with both displacements and their derivatives as nodal variables results in a large number of degrees of freedom. Here nodal rotations are applied along with nodal displacements to obtain a simpler finite element formulation with fewer degrees of freedom. The difference in rotation gradients determined with nodal displacements and rotations are minimized by a penalty term. To assess the suggested approach simulations are presented and discussed, where the material parameters have been obtained from experiments of epoxy microbeams in the literature.     

Structural-Damage Detection by Distributed Piezoelectric Transducers and Tuned Electric Circuits

ABSTRACT

A novel technique for damage detection of structures is introduced and discussed. It is based on purely electric measurements of the state variables of an electric network coupled to the main structure through a distributed set of piezoelectric patches. The constitutive parameters of this auxiliary network are optimized to increase the sensitivity of global measurements—as the frequency, response functions relative to selected electric degrees of freedom—with respect to a given class of variations in the structural–mechanical properties. Because the proposed method is based on purely electric input and output measurements, it allows for accurate results in the identification and localization of damages. Use of the electric frequency-response function to identify the mechanical damage leads to nonconvex optimization problems; therefore the proposed sensitivity-enhanced identification procedure becomes computationally efficient if an a priori knowledge about the damage is available.     

Improved assumed-stress hybrid shell element with drilling degrees of freedom for linear stress,buckling and free vibration analyses

An improved 4-node quadrilateral assumed-stress hybrid shell element with drilling degrees of freedom is presented. The formulation is based on Hellinger–Reissner variational principle and the shape functions are formulated directly for the 4-node element. The element has 12 membrane degrees of freedom and 12 bending degrees of freedom. It has 9 independent stress parameters to describe the membrane stress resultant field and 13 independent stress parameters to describe the moment and transverse shear stress resultant field. The formulation encompasses linear stress, linear buckling and linear free vibration problems. The element is validated with standard test cases and is shown to be robust. Numerical results are presented for linear stress, buckling, and free vibration analyses.     

Ductile damage analysis of elasto‐plastic shells at large inelastic strains

The objective of this contribution is to model ductile damage phenomena under consideration of large inelastic strains, to couple the corresponding constitutive law with a multi‐layer shell kinematics and to give finally an adequate finite element formulation. An elastic–plastic constitutive law is formulated by using a spatial hyperelasto‐plastic formulation based on the multiplicative decomposition of the deformation gradient. To include isotropic ductile damage the continuum damage model of Rousselier is modified so as to consider large strains and additionally extended by various void nucleation and macro‐crack criteria. In order to achieve numerical efficiency, elastic strains are supposed to be sufficiently small providing a numerical effective integration based on the backward Euler rule. Finite element formulation is enriched by means of the enhanced strain concept. Thus the well‐known deficiencies due to incompressible deformations and the inclusion of transverse strains are avoided. Several examples are given to demonstrate the performance of the algorithms developed concerning large inelastic strains of shells and ductile damage phenomena. Copyright © 2000 John Wiley & Sons, Ltd.