Development of a continuum plasticity model for the commercial finite element code ABAQUS

The present work relates to the development of computational material models for sheet metal forming simulations. In this specific study, an implicit scheme with consistent Jacobian is used for integration of large deformation formulation and plane stress elements. As a privilege to the explicit scheme, the implicit integration scheme is unconditionally stable. The backward Euler method is used to update trial stress values lying outside the yield surface by correcting them back to the yield surface at every time increment. In this study, the implicit integration of isotropic hardening with the von Mises yield criterion is discussed in detail. In future work it will be implemented into the commercial finite element code ABAQUS by means of a user material subroutine

High order unfitted finite element methods on level set domains using isoparametric mappings

We introduce a new class of unfitted finite element methods with high order accurate numerical integration over curved surfaces and volumes which are only implicitly defined by level set functions. An unfitted finite element method which is suitable for the case of piecewise planar interfaces is combined with a parametric mapping of the underlying mesh resulting in an isoparametric unfitted finite element method. The parametric mapping is constructed in a way such that the quality of the piecewise planar interface reconstruction is significantly improved allowing for high order accurate computations of (unfitted) domain and surface integrals. This approach is new. We present the method, discuss implementational aspects and present numerical examples which demonstrate the quality and potential of this method.Comment: 18 pages, 8 figure

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme

High-Order Unstructured Lagrangian One-Step WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems: Applications to Compressible Multi-Phase Flows

In this article we present the first better than second order accurate unstructured Lagrangian-type one-step WENO finite volume scheme for the solution of hyperbolic partial differential equations with non-conservative products. The method achieves high order of accuracy in space together with essentially non-oscillatory behavior using a nonlinear WENO reconstruction operator on unstructured triangular meshes. High order accuracy in time is obtained via a local Lagrangian space-time Galerkin predictor method that evolves the spatial reconstruction polynomials in time within each element. The final one-step finite volume scheme is derived by integration over a moving space-time control volume, where the non-conservative products are treated by a path-conservative approach that defines the jump terms on the element boundaries. The entire method is formulated as an Arbitrary-Lagrangian-Eulerian (ALE) method, where the mesh velocity can be chosen independently of the fluid velocity. The new scheme is applied to the full seven-equation Baer-Nunziato model of compressible multi-phase flows in two space dimensions. The use of a Lagrangian approach allows an excellent resolution of the solid contact and the resolution of jumps in the volume fraction. The high order of accuracy of the scheme in space and time is confirmed via a numerical convergence study. Finally, the proposed method is also applied to a reduced version of the compressible Baer-Nunziato model for the simulation of free surface water waves in moving domains. In particular, the phenomenon of sloshing is studied in a moving water tank and comparisons with experimental data are provided

Numerical integration of the incrementally non-linear, zero elastic range, bounding surface plasticity model for sand

SANISAND-Z is a recently developed plasticity model for sands with zero purely elastic range in stress space within the framework of Bounding Surface (BS) plasticity. As a consequence of zero elastic range the plastic strain increment direction, and consequently the elastic-plastic moduli fourth order tensor depends on the direction of the stress increment, rendering the model incrementally non-linear and intrinsically implicit. An iterative algorithm based on the Backward Euler method is presented to solve the non-linear system of ordinary differential equations. A non-traditional consistency condition based on the plastic multiplier is introduced as a core element of the system. A thorough analysis of the stability and accuracy of the algorithm is presented based on error estimation. The proposed integration scheme allows the use of SANISAND-Z framework in Finite Element Analysis