a deformation dependent stabilization technique exemplified by eas elements at large strains abstract stabilized finite element methods have been developed mainly in the context of computational fluid dynamics cfd and have shown to be able to add stability to previously unstable formulations in a consistent way in this contribution a deformation dependent stabilization technique conceptually based on the above mentioned developments in the cfd area is developed for solid mechanics to cure the well known enhanced assumed strain eas method from artificial instabilities hourglass modes that have been observed in the range of large compressive strains in investigating the defect of the original formulation the dominating role of the kinematic equation as cause for the instabilities is revealed this observation serves as key ingredient for the design of the stabilizing term introduced on the level of the variational equation a proper design for the stabilization parameter is given based on a mechanical interpretation of the underlying defect as well as of the stabilizing action this stabilizing action can be thought of an additional constraint introduced into the reparametrized hu washizu functional in a least square form together with a deformation dependent stabilization parameter numerical examples show the capability of this approach to effectively eliminate spurious hourglass modes which otherwise may appear in the presence of large compressive strains while preserving the advantageous features of the eas method namely the reduction of the stiffness for an in plane bending mode i e when plane stress elements are used in a bending situation

multiscale phenomena green s functions the dirichlet to neumann formulation subgrid scale models bubbles and the origins of stabilized methods abstract an approach is developed for deriving variational methods capable of representing multiscale phenomena the ideas are first illustrated on the exterior problem for the helmholtz equation this leads to the well known dirichlet to neumann formulation next a class of subgrid scale models is developed and the relationships to bubble function methods and stabilized methods are established it is shown that both the latter methods are approximate subgrid scale models the identification for stabilized methods leads to an analytical formula for the intrinsic time scale whose origins have been a mystery heretofore

a 4 node finite shell element for the implementation of general hyperelastic 3d elasticity at finite strains abstract in this paper a finite shell element for large deformations is presented based on extensible director kinematics the essential feature is an interface to arbitrary three dimensional material laws the non linear lagrangian formulation is based on the three field variational principle parametrized with the displacement vector enhanced green lagrangian strain tensor and second piola kirchhoff stress tensor the developed quadrilateral shell element is characterized by a course mesh accuracy and distortion insensitivity compared with bilinear displacement approaches furthermore plane stress response is approximately recovered in the asymptotic case of vanishing thickness a number of example problems investigating large deformation as well as finite strain applications are presented compressible and incompressible hyperelastic materials of the st venant kirchhoff neo hookean and mooney rivlin type are particularly used

bubble stabilization of finite element methods for the linearized incompressible navier stokes equations abstract in this paper we discuss the stabilization via bubble functions of a finite element method for stationary linearized incompressible navier stokes equations it is shown that residual free bubbles can reproduce supg

assumed strain stabilization procedure for the 9 node lagrange shell element an assumed strain strain interpolation method is used to construct a stabilization matrix for the 9 node shell element the stabilization procedure can be justified based on the hellinger reissner variational method it involves a projection vector which is orthogonal to both linear and quadratic fields in the local co ordinate system of each quadrature point all terms in the development involve 2 2 quadrature in the 9 node element example problems show good accuracy and an almost optimal rate of convergence

geometrically non linear enhanced strain mixed methods and the method of incompatible modes a class of assumed strain mixed finite element methods for fully non linear problems in solid mechanics is presented which when restricted to geometrically linear problems encompasses the classical method of incompatible modes as a particular case the method relies crucially on a local multiplicative decomposition of the deformation gradient into a conforming and an enhanced part formulated in the context of a three field variational formulation the resulting class of mixed methods provides a possible extension to the non linear regime of well known incompatible mode formulations in addition this class of methods includes non linear generalizations of recently proposed enhanced strain interpolations for axisymmetric problems which cannot be interpreted as incompatible modes elements the good performance of the proposed methodology is illustrated in a number of simulations including 2 d 3 d and axisymmetric finite deformation problems in elasticity and elastoplasticity remarkably these methods appear to be specially well suited for problems involving localization of the deformation as illustrated in several numerical examples

a new 3 d finite element for nonlinear elasticity using the theory of a cosserat point abstract the theory of a cosserat point has been used to formulate a new 3 d finite element for the numerical analysis of dynamic problems in nonlinear elasticity the kinematics of this element are consistent with the standard tri linear approximation in an eight node brick element specifically the cosserat point is characterized by eight director vectors which are determined by balance laws and constitutive equations for hyperelastic response the constitutive equations for the director couples are determined by derivatives of a strain energy function restrictions are imposed on the strain energy function which ensure that the element satisfies a nonlinear version of the patch test it is shown that the cosserat balance laws are in one to one correspondence with those obtained using a bubnov galerkin formulation nevertheless there is an essential difference between the two approaches in the procedure for obtaining the strain energy function specifically the cosserat approach determines the constitutive coefficients for inhomogeneous deformations by comparison with exact solutions or experimental data in contrast the bubnov galerkin approach determines these constitutive coefficients by integrating the 3 d strain energy function using the kinematic approximation it is shown that the resulting cosserat equations eliminate unphysical locking and hourglassing in large compression without the need for using assumed enhanced strains or special weighting functions

consistent tangent operators for rate independent elastoplasticity abstract it is shown that for problems involving rate constitutive equations such as rate independent elastoplasticity the notion of consistency between the tangent stiffness operator and the integration algorithm employed in the solution of the incremental problem plays a crucial role in preserving the quadratic rate of asymptotic convergence of iterative solution schemes based upon newton s method within the framework of closest point projection algorithms a methodology is presented whereby tangent operators consistent with this class of algorithms may be systematically developed to wit associative j 2 flow rules with general nonlinear kinematic and isotropic hardening rules as well as a class of non associative flow rules are considered the resulting iterative solution scheme preserves the asymptotic quadratic convergence characteristic of newton s method whereas use of the socalled elastoplastic tangent in conjunction with a radial return integration algorithm a procedure often employed results in newton type of algorithms with suboptimal rate of convergence application is made to a set of numerical examples which include saturation hardening laws of exponential type

an adaptive stabilization strategy for enhanced strain methods in non linear elasticity this paper proposes and analyzes an adaptive stabilization strategy for enhanced strain es methods applied to quasistatic non linear elasticity problems the approach is formulated for any type of enhancements or material models and it is distinguished by the fact that the stabilization term is solution dependent the stabilization strategy is first constructed for general linearized elasticity problems and then extended to the non linear elastic regime via an incremental variational principle a heuristic choice of the stabilization parameters is proposed which in the numerical examples proved to provide stable approximations for a large range of deformations different problems and material models we also provide explicit lower bounds for the stabilization parameters that guarantee that the method will be stable these are not advocated since they are generally larger than the ones based on heuristics and hence prone to deteriorate the locking free behavior of es methods numerical examples with two different non linear elastic models in thin geometries and incompressible situations show that the method remains stable and locking free over a large range of deformations finally the method is strongly based on earlier developments for discontinuous galerkin methods and hence throughout the paper we offer a perspective about the similarities between the two

design of galerkin generalized least squares methods for timoshenko beams abstract a class of finite element methods the galerkin generalized least squares methods are developed and applied to model the steady state response of timoshenko beams an optimal method is designed using a linear interpolation of the response such that there is zero finite element dispersion error the classical method of selective reduced integration and a modified version of selective reduced integration with mass lumping are shown to fall under the galerkin generalized least squares framework numerical experiments in wave propagation demonstrate the dramatic superiority of the new optimal method over the standard approaches the goal of the new methods is to decrease the computational burden required to achieve a desired accurcy level at a particular frequency thereby enabling larger scale higher frequency computations for a given platform

on inverse scattering for then body schrodinger equation abstract we consider the generalized n body schrodinger operator h q x q x a a v a x a x r d x a x a r d x a x a with short range regular interactions we obtain in particular the following results 1 we give new formulas and proofs relating high energy asymptotics in a weak sense of the scattering operators with the same cluster decomposition a for and and the x ray transform pi defined on the set of all lines in x a with nonsingular directions of the effective potentials i x a x a x a these results significantly clarify some of those given in the literature 2 we describe completely ker p and give a method for reconstruction of i mod ker p from pi 3 we prove pointwise high energy asymptotics for the two cluster two cluster scattering amplitudes f with the fourier transform i in the leading term for the case when the cluster decomposition for and is the same 4 we give several additional results for the case of perturbed stratified medium when q x v a x a v b x b x a x 1 x b x 1 x d x and each v c x c rapidly decreases as x c c a b our results including proofs in some cases significantly simplify methods of high energy inverse scattering for the n body schrodinger operator given earlier by wang and by enss and weder

finite element methods for linear hyperbolic problems abstract we give a survey of some recent work by the authors on finite element methods for convectiondiffusion problems and first order linear hyperbolic problems

error analysis of some galerkin least squares methods for the elasticity equations the recent technique of stabilizing mixed finite element methods by augmenting the galerkin formulation with least squares terms calculated separately on each element is considered the error analysis is performed in a unified manner yielding improved results for some methods introduced earlier in addition a new formulation is introduced and analyzed method is used within two different contexts the first is in connection with the plain strain or three dimensional problems near the incompressible limit for which the poisson ratio equals one half it is well known that for these problems a standard lower order displacement finite element method will break down by now it is also well known how a working method can be obtained a new unknown the pressure is introduced and discretized with a finite element space different from that used for the displacements the finite element method obtained falls into the class of saddlepoint problems for which there is a general theory developed by brezzi and babuska 5 7 the application of this theory shows that the method converges optimally if the discrete spaces for the displacement and pressure satisfy a stability condition for a particular choice of displacement and pressure approximations it can be a difficult task to verify this babuska brezzi or inf sup condition during the last decade this problem has been intensively studied and the performances of many combinations are now known cf 18 and 23 26 for some recent extensions in addition there exists a general technique utilizing bubble functions by which a stable method can

stabilized finite element method for incompressible flows with high reynolds number in the following paper we discuss the exhaustive use and implementation of stabilization finite element methods for the resolution of the 3d time dependent incompressible navier stokes equations the proposed method starts by the use of a finite element variational multiscale vms method which consists in here of a decomposition for both the velocity and the pressure fields into coarse resolved scales and fine unresolved scales this choice of decomposition is shown to be favorable for simulating flows at high reynolds number we explore the behaviour and accuracy of the proposed approximation on three test cases first the lid driven square cavity at reynolds number up to 50 000 is compared with the highly resolved numerical simulations and second the lid driven cubic cavity up to re 12 000 is compared with the experimental data finally we study the flow over a 2d backward facing step at re 42 000 results show that the present implementation is able to exhibit good stability and accuracy properties for high reynolds number flows with unstructured meshes

hourglass control in linear and nonlinear problems abstract mesh stabilization techniques for controlling the hourglass modes in under integrated hexahedral and quadrilateral elements are described it is shown that the orthogonal hourglass techniques previously developed can be obtained from simple requirements that insure the consistency of the finite element equations in the sense that the gradients of linear fields are evaluated correctly it is also shown that this leads to an hourglass control that satisfies the patch test the nature of the parameters which relate the generalized stresses and strains for controlling hourglass modes is examined by means of a mixed variational principle and some guidelines for their selection are discussed finally effective means of implementing these hourglass procedure in computer codes are described applications to both the laplace equation and the equations of solid mechanics in 2 and 3 dimensions are considered

exponential bounds and absence of positive eigenvalues for n body schrodinger operators for a large class ofn body potentialsv we prove that if is an eigenfunction of v with eigenvaluee then sup 2 e 0 exp x l2 is either a threshold or consequences of this result are the absence of positive eigenvalues and optimal l2 exponential lower bounds

a class of mixed assumed strain methods and the method of incompatible modes a three field mixed formulation in terms of displacements stresses and an enhanced strain field is presented which encompasses as a particular case the classical method of incompatible modes within this frame work incompatible elements arise as particular compatible mixed approximations of the enhanced strain field the conditions that the stress interpolation contain piece wise constant functions and be l2 ortho gonal to the enhanced strain interpolation ensure satisfaction of the patch test and allow the elimination of the stress field from the formulation the preceding conditions are formulated in a form particularly convenient for element design as an illustration of the methodology three new elements are developed and shown to exhibit good performance a plane 3d elastic plastic quad an axisymmetric element and a thick plate bending quad the formulation described herein is suitable for non linear analysis

three dimensional extension of non linear shell formulation based on the enhanced assumed strain concept conventional shell formulations such as 3 or 5 parameter theories or even 6 parameter theories including the thickness change as extra parameter require a condensation of the constitutive law in order to avoid a significant error due to the assumption of a linear displacement field across the thickness this means that the normal stress in thickness direction has to either vanish or be constant in general these extra constraints cannot be satisfied explicitly or they lead to elaborate strain expressions the main objective of the present study is to introduce directly a complete 3 d constitutive law without modification therefore a 7 parameter theory is utilized which includes a linear varying thickness stretch as extra variable allowing also large strain effects

a note on enhanced strain methods for large deformations abstract enhanced strain elements have been applied to simulate geometrically and materially non linear problems their main advantage is due to the fact that they perform very good in the incompressible limit as well as in bending situations however in compressive deformation states these elements depict stability modes which are associated with hour glass forms and thus denote a rank deficiency for such deformations in this paper the appearance of these modes will be investigated analytically by means of a simple representative example

streamline upwind petrov galerkin formulations for convection dominated flows with particular emphasis on the incompressible navier stokes equations abstract a new finite element formulation for convection dominated flows is developed the basis of the formulation is the streamline upwind concept which provides an accurate multidimensional generalization of optimal one dimensional upwind schemes when implemented as a consistent petrov galerkin weighted residual method it is shown that the new formulation is not subject to the artificial diffusion criticisms associated with many classical upwind methods the accuracy of the streamline upwind petrov galerkin formulation for the linear advection diffusion equation is demonstrated on several numerical examples the formulation is extended to the incompressible navier stokes equations an efficient implicit pressure explicit velocity transient algorithm is developed which accomodates several treatments of the incompressibility constraint and allows for multiple iterations within a time step the effectiveness of the algorithm is demonstrated on the problem of vortex shedding from a circular cylinder at a reynolds number of 100