Search results

This paper proposes a unified original general framework, designed to theoretically develop and to extremely easily implement elastoplastic constitutive laws defined in the so called two-invariants space, both in small and finite strain regime.

A general return mapping algorithm is proposed, and particularly a standard procedure is developed to compute the two algorithmic tangent operators, required to solve the Newton–Raphson scheme at the local and global level and thus cast the elastoplastic algorithm within a FEM code.

This work demonstrates that the proposed procedure is fully general and can be applied whatever is the elastic law, the yield surface, the plastic potential function and the hardening law. Several numerical examples are reported, not only to demonstrate the accuracy and robustness of the algorithm, but also explain how to use this general algorithm also in other applications.

The proposed algorithm and its numerical implementation into a FEM code is new and original. The usefulness and the value of the algorithm is twofold: (1) it can be implemented in a small and finite strain simulation FEM code, in order to handle different types of constitutive laws in the same modular way, thus fully leveraging on modern object-oriented coding approach; (2) it can be used as a framework to develop (and then to implement) new constitutive models, since the researcher can simply define the relevant functions (and its main derivatives) and automatically get the numerical algorithm.

The objective of this work is to develop an element technology to recover the plane stress response without any plane stress specific modifications in the large strain regime. Therefore, the essential feature of the proposed element formulation is an interface to arbitrary three‐dimensional constitutive laws. The easily implemented and computational cheap four‐noded element is characterized by coarse mesh accuracy and the satisfaction of the plane stress constraint in a weak sense. A number of example problems involving arbitrary small and large strain constitutive models demonstrate the excellent performance of the concept pursued in this work.

Under restriction of an isotropic elastic response of deformed lattice, develops a covariant theory of finite elastoplasticity in principal axes of a pair of deformation tensors. In material description, the tensor pair consists of the plastic deformation tensor and the total deformation Cauchy‐Green tensor. Applies the proposed theory to elastoplastic membrane shells, whose references and current configurations can be arbitrary space‐curved surfaces. Pressure‐insensitive von Mises yield criterion with isotropic hardening and a quadratic form of the strain energy function given in terms of elastic principal stretches are considered as a model problem. Through an explicit enforcement of the plane stress condition we arrive at a reduced two‐dimensional problem representation, which is set in the membrane tangent plane. Numerical implementation of the presented theory relies crucially on the operator split methodology to simplify the state update computation. Presents a set of numerical examples in order to illustrate the performance of the presented methodology and indicate possible applications in the area of sheet metal forming.

Many industrial analyses require the resolution of complex nonlinear problems. For such calculations, error‐controlled adaptive strategies must be used to improve the quality of the results. In this paper, adaptive strategies for nonlinear calculations in plasticity based on an enhanced error on the constitutive relation are presented. We focus on the adaptivity of the mesh and of the time discretization.

An efficient boundary element method (BEM) formulation for three‐dimensional elastoplasticity is presented in this paper. The BEM formulation for nonlinear problems requires discretization of the surface as well as part of the volume. In this paper nine‐noded quadrilateral elements are used for modelling the surface and 27‐noded brick elements for the volume. Particular attention is paid to the accurate evaluation of the Cauchy principal value volume integrals appearing in the interior stress calculations. An explicit initial strain formulation is used to satisfy the non‐linearity. The accuracy of the proposed method is demonstrated by solving a number of benchmark problems.

Plane strain constitutive behaviour of von Mises and isotropic Hoffman materials is examined using single element tests. Two kinds of tests are conducted – (a) prescribed displacement tests; and (b) tests with a mixture of displacements and boundary tractions prescribed. While (a) are used to understand the manner of stress traversal on the yield surface in principal stress space, (b) are employed to study the load displacement response and the possibility of ensuing localization. Associated plasticity is assumed throughout. The tests are conducted using perfect and strain softening plasticity. It is found that for the von Mises criterion limited exact solutions can be evolved even under softening (or hardening) conditions. For isotropic Hoffman materials the nature of the stress traversal, load deflection response and the satisfaction of the localization conditions are strongly influenced by the ratio and difference of uniaxial yield strengths, in tension and compression, as well as by the softening parameters.

The numerical simulation of metal forming processes follows a highly non‐linear analysis where general aspects such as elastoplasticity, finite deformation and contact mechanics are combined. Approximated solutions obtained by finite element techniques require strong computational effort, that contradicts the need of interactive industrial applications. The first part of the work deals with the description of the main elements of the formulation, with attention to mathematical modelling and the approximating algorithms in the incremental iterative frame of non‐linear analysis, ending with the results obtained in hot rolling simulation. The second part is dedicated to computational efficiency analysis and the presentation of the related methods and results obtained in this work.

The finite element quasi‐static analysis of elastoplastic systems is studied by making use of a generalized variable approach for the spatial discretization and a generalized mid‐point rule for the time integration. Both the classical form of the constitutive law and the convex analysis formulation are presented. The relation between the mid‐point time integration and the extremal path theory is discussed. Extremal properties for the finite‐step solution are formulated.

A numerical formulation for solving homogeneous anisotropic heat conduction problems based on the use of an isotropic fundamental solution is presented in detail. The analysis is carried out assuming a generic position of the coordinate axes, which may not coincide with the principal directions of orthotropy of the material. The two primary integral equations of the method are derived from the governing differential equation of the problem. Then, the numerical procedure is developed by rewriting the internal degrees of freedom that arise from the domain discretization in terms of the boundary nodes and solving the resulting system of linear equations for the boundary unknowns only. Special attention is given to the differentiation of singular integrals which yields additional terms as well as to the evaluation of the resulting Cauchy principal value integral. The main feature of the proposed formulation is its generality, which makes possible its direct extension to solve the problem of three‐dimensional heat conduction in anisotropic media and, foremost, to three‐dimensional orthotropic and anisotropic elasticity or elastoplasticity.

The boundary element method (BEM) and the finite element method (FEM) may be computationally expensive if complex problems are to be solved; thus there is the need of implementing them on fast computer architectures, especially parallel computers. Because these methods are complementary to each other, the coupling of FEM and BEM is widely used. In this paper, the coupling of displacement‐based FEM and collocation BEM and its implementation on a distributed memory system (Parsytec MultiCluster2) is described. The parallelization is performed by data partitioning which leads to a very high efficiency. As model problems, we assume linear elasticity for the boundary element method and elastoplasticity for the finite element method. The efficiency of our implementation is shown by various test examples. By numerical examples we show that a multiplicative Schwarz method for coupling BEM with FEM is very well suited for parallel implementation.