Miroslav Halilovič

NICE-An explicit numerical scheme for efficient integration of nonlinear constitutive equations

The paper presents a simple but efficient new numerical scheme for the integration of nonlinear constitutive equations. Although it can be used for the integration of a system of algebraic and differential equations in general, the scheme is primarily developed for use with the direct solution methods for solving boundary value problems, e.g. explicit dynamic analysis in ABAQUS/Explicit. In the developed explicit scheme, where no iteration is required, the implementation simplicity of the forward-Euler scheme and the accuracy of the backward-Euler scheme are successfully combined. The properties of the proposed NICE scheme, which was also implemented into ABAQUS/Explicit via User Material Subroutine (VUMAT) interface platform, are compared with the properties of the classical forward-Euler scheme and backward-Euler scheme. For this purpose two highly nonlinear examples, with the von Mises and GTN material model considered, have been studied. The accuracy of the new scheme is demonstrated to be at least of the same level as experienced by the backward-Euler scheme, if we compare them on the condition of the same CPU time consumption. Besides, the simplicity of the NICE scheme, which is due to implementation similarity with the classical forward-Euler scheme, is its great Advantage.

Earing Prediction in Cup Drawing using the BBC2008 Yield Criterion

The paper deals with constitutive modelling of highly anisotropic sheet metals. It presents FEM based earing predictions in cup drawing simulation of highly anisotropic aluminium alloys where more than four ears occur. For that purpose the BBC2008 yield criterion, which is a plane‐stress yield criterion formulated in the form of a finite series, is used. Thus defined criterion can be expanded to retain more or less terms, depending on the amount of given experimental data. In order to use the model in sheet metal forming simulations we have implemented it in a general purpose finite element code ABAQUS/Explicit via VUMAT subroutine, considering alternatively eight or sixteen parameters (8p and 16p version). For the integration of the constitutive model the explicit NICE (Next Increment Corrects Error) integration scheme has been used. Due to the scheme effectiveness the CPU time consumption for a simulation is comparable to the time consumption of built‐in constitutive models. Two aluminium alloys, namely...

Impact Of Elastic Modulus Degradation On Springback In Sheet Metal Forming

Strain recovery after removal of forming loads, commonly defined as springback, is of great concern in sheet metal forming, in particular with regard to proper prediction of the final shape of the part. To control the problem a lot of work has been done, either by minimizing the springback on the material side or by increasing the estimation precision in corresponding process simulations. Unfortunately, by currently available software springback still cannot be adequately predicted, because most analyses of springback are using linear, isotropic and constant Young’s modulus and Poisson’s ratio. But, as it was measured and reported, none of it is true. The aim of this work is to propose an upgraded mechanical model which takes evolution of damage and related orthotropic stiffness degradation into account. Damage is considered by inclusion of ellipsoidal cavities, and their influence on the stiffness degradation is taken in accordance with the Mori‐Tanaka theory, adopting the GTN model for plastic flow. In ...

NICE h : a higher-order explicit numerical scheme for integration of constitutive models in plasticity

The article introduces, as a result of further development of the first-order scheme NICE, a simple and efficient higher-order explicit numerical scheme for the integration of a system of ordinary differential equations which is constrained by an algebraic condition (DAE). The scheme is based on the truncated Taylor expansion of the constraint equation with order h of the scheme being determined by the highest exponent in the truncated Taylor series. The integration scheme thus conceived will be named NICEh, considering both principal premises of its construction. In conjunction with a direct solution technique used to solve the boundary value problem, the NICEh scheme is very convenient for integrating constitutive models in plasticity. The plasticity models are defined mostly by a system of algebraic and differential equations in which the yield criterion represents the constraint condition. To study the properties of the new integration scheme, which, like the forward-Euler scheme, is characterised by its implementation simplicity due to the explicitness of its formulations, a damage constitutive model (Gurson–Tvergaard–Needleman model) is considered. The general opinion that the implicit backward-Euler scheme is much more accurate than the thus-far known explicit schemes is challenged by the introduction of the NICEh scheme. The accuracy of the higher-order explicit scheme in the studied cases is significantly higher than the accuracy of the classical backward-Euler scheme, if we compare them under the condition of a similar CPU time consumption.

A robust explicit integration of elasto-plastic constitutive models, based on simple subincrement size estimation

Purpose The purpose of this study, which is designed for the implementation of models in the implicit finite element framework, is to propose a robust, stable and efficient explicit integration algorithm for rate-independent elasto-plastic constitutive models. Design/methodology/approach The proposed automatic substepping algorithm is founded on an explicit integration scheme. The estimation of the maximal subincrement size is based on the stability analysis. Findings In contrast to other explicit substepping schemes, the algorithm is self-correcting by definition and generates no cumulative drift. Although the integration proceeds with maximal possible subincrements, high level of accuracy is attained. Algorithmic tangent stiffness is calculated in explicit form and optionally no analytical second-order derivatives are needed. Research limitations/implications The algorithm is convenient for elasto-plastic constitutive models, described with an algebraic constraint and a set of differential equations. This covers a large family of materials in the field of metal plasticity, damage mechanics, etc. However, it cannot be directly used for a general material model, because the presented algorithm is convenient for solving a set of equations of a particular type. Practical implications The estimation of the maximal stable subincrement size is computationally cheap. All expressions in the algorithm are in explicit form, thus the implementation is simple and straightforward. The overall performance of the approach (i.e. accuracy, time consumption) is fully comparable with a default (built-in) ABAQUS/Standard algorithm. Originality/value The estimated maximal subincrement size enables the algorithm to be stable by definition. Subincrements are much larger than those in conventional substepping algorithms. No error control, error correction or local iterations are required even in the case of large increments.

An approach to consider the arterial residual stresses in modelling of a patient-specific artery

In this work, the residual stress state of a human common carotid artery is predicted using the so-called thermomechanical analogy approach. The purpose of the approach is to enable consistent mapping of residual stresses and the respective configuration from a circular arterial segment to a patient-specific arterial geometry. This is achieved by applying proper volumetric dilatations to the actual arterial stress-free in vivo geometry, which makes use of the analogy that states that the bending stresses can be obtained on an equivalent manner by applying proper thermal dilatations. The common carotid artery data are obtained in vivo from a healthy 28-year-old man using non-invasive methods. The predicted residual stresses of the common carotid artery are in good quantitative agreement with the data from prior work in this field. The approach is validated by predicting the common carotid artery zero-stress state configuration, where a sector-like (cut-open) state is obtained. With this approach, it is thus possible to predict the residual stresses as well as the configuration of patient-specific arterial geometry without the need to model its cut-open zero-stress configuration.

Characterization of Mechanical Behaviour of Sheet Metal: Hardening at Large Strains and Anisotropy Identification

The paper focuses on a reliable identification of the hardening curve of sheet metal after the onset of necking by means of a purposely developed shear test. The main comparative advantage of the test is its simple implementation needing no further development of the experimental equipment used in standard tensile testing. In addition, in order to achieve more flexible hardening curve description from those given by the classical analytical hardening laws, the paper presents the hardening curve approximation based on the cubic spline technique. The identified hardening curves are validated and compared to the analytically extrapolated ones obtained by a continuation of classical hardening laws for large plastic strain. It is shown, that due to large difference in response the analytical extrapolation is seriously disputable.

Approaches To Modelling Of Elastic Modulus Degradation In Sheet Metal Forming

Strain recovery after removal of forming loads, commonly defined as springback, is of great concern in sheet metal forming, in particular with regard to proper prediction of the final shape of the part. To control the problem a lot of work has been done, either by minimizing the springback on the material side or by increasing the estimation precision in corresponding process simulations. Unfortunately, by currently available software springback still cannot be adequately predicted, because most analyses of springback are using linear, isotropic and constant Young’s modulus and Poisson’s ratio. But, as it was measured and reported, none of it is true. The aim of this work is to propose an upgraded mechanical model which takes evolution of damage and related orthotropic stiffness degradation into account. Damage is considered by inclusion of ellipsoidal cavities, and their influence on the stiffness degradation is taken in accordance with the Mori‐Tanaka theory, adopting the GTN model for plastic flow. Wit...