Search results

The purpose of this paper is to present an efficient return mapping algorithm for elastoplastic constitutive problems of ductile metals with an exact closed form solution of the local constitutive problem in the small strain regime. A Newton Raphson iterative method is adopted for the solution of the boundary value problem.

An efficient return mapping algorithm is illustrated which is based on an elastic predictor and a plastic corrector scheme resulting in an implicit and accurate numerical integration method. Nonlinear kinematic hardening rules and linear isotropic hardening rules are used to describe the components of the hardening variables. In the adopted algorithmic approach the solution of the local constitutive equations reduces to only one straightforward nonlinear scalar equation.

The presented algorithmic scheme naturally leads to a particularly simple form of the nonlinear scalar equation which ultimately scales down to an algebraic (polynomial) equation with a single variable. The straightforwardness of the present approach allows to find the analytical solution of the algebraic equation in a closed form. Further, the consistent tangent operator is derived as associated with the proposed algorithmic scheme and it is shown that the proposed computational procedure ensures a quadratic rate of asymptotic convergence when used with a Newton Raphson iterative method for the global solution procedure.

In the present approach the solution of the algebraic nonlinear equation is found in a closed form and accordingly no iterative method is required to solve the problem of the local constitutive equations. The computational procedure ensures a quadratic rate of asymptotic convergence for the global solution procedure typical of computationally efficient solution schemes. In the paper it is shown that the proposed algorithmic scheme provides an efficient and robust computational solution procedure for elastoplasticity boundary value problems. Numerical examples and computational results are reported which illustrate the effectiveness and robustness of the adopted integration algorithm for the finite element analysis of elastoplastic structures also under elaborate loading conditions.

Many practical problems in engineering require fast, accurate numerical results. In particular, in cyclic plasticity or fatigue simulations, the high number of loading cycles increases the computation effort and time. The purpose of this study is to show that the return mapping technique in the framework of unconventional plasticity theories is a good compromise between efficiency and accuracy in finite element analyses.

The accuracy of the closest point projection method and the cutting plane method implementations for the subloading surface model are discussed under different loading conditions by analyzing the error as a function of the input step size and the efficiency of the algorithms.

Monotonic tests show that the two different implicit integration schemes have the same accuracy and are in good agreement with the solution obtained using an explicit forward Euler scheme, even for large input steps. However, the closest point projection method seems to describe better the evolution of the similarity centre in the cyclic loading analyses.

The purpose of this work is to show two alternative implicit integration schemes of the extended subloading surface method for metallic materials. The backward Euler integrations can guarantee a good description of the material behaviour and, at the same time, reduce the computational cost. This aspect is particularly important in the field of low or high cycle fatigue, because of the large number of cycles involved.

A detailed description of both the cutting plane and closest point projection methods is offered in this work. In particular, the two integrations schemes are compared in terms of accuracy and computation time for monotonic and cyclic loading tests.

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.

In Simo and Taylor, the classical radial return algorithm of Wilkins and Krieg and Key for plane strain and three‐dimensional J2‐flow theory, is extended to the case of plane stress. In three dimensions (or plane strain), enforcement of the discrete consistency condition reduces to a simple radial scaling of the trial stress onto the yield surface; i.e., the return map is radial. In plane stress, on the other hand, the return map, that restores the trial stress back to the yield surface, is constrained to remain in the plane stress subspace, and thus no longer reduces to a simple radial scaling. The determination of the final stress point from the trial stress now involves the solution by Newton's method of a non‐linear scalar equation, referred to as the discrete consistency equation in what follows, that yields the discrete consistency parameter λn+>0. The requirement that λn+>1 be positive is a direct consequence of the discrete Kuhn‐Tucker optimality conditions.

Recent progress in element technology in large scale explicit finite element codes has opened the way for the solution of elastoplastic shell problems of unprecedented complexity. This new capability has focused attention on the numerical issues involved in the implementation of elastoplastic material models for shells, particularly when vectorizable algorithms are required for supercomputer applications. This paper reviews four algorithms currently in the literature for plane stress and shell plasticity. First, each of the four methods is described in detail. Next, an accuracy analysis is presented for each algorithm for perfectly plastic, linear kinematic hardening, and linear isotropic hardening cases. Finally, a comparison is made of the relative computational efficiency of the four algorithms, and the importance of vectorization is illustrated.

Semi-implicit type cutting plane method (CPM) and fully implicit type closest point projection method (CPPM) are the two most widely used frameworks for numerical stress integration. CPM is simple, easy to implement and accurate up to first order. CPPM is unconditionally stable and accurate up to second order though the formulation is complex. Therefore, this study aims to develop a less complex and accurate stress integration method for complex constitutive models.

Two integration techniques are formulated using the midpoint and Romberg method by modifying CPM. The algorithms are implemented for three different classes of soil constitutive model. The efficiency of the algorithms is judged via stress point analysis and solving a boundary value problem.

Stress point analysis indicates that the proposed algorithms are stable even with a large step size. In addition, numerical analysis for solving boundary value problem demonstrates a significant reduction in central processing unit (CPU) time with the use of the semi-implicit-type midpoint algorithm.

Traditionally, midpoint and Romberg algorithms are formulated from explicit integration techniques, whereas the present study uses a semi-implicit approach to enhance stability. In addition, the proposed stress integration algorithms provide an efficient means to solve boundary value problems pertaining to geotechnical engineering.

Presents a robust and unconditionally stable return‐mapping algorithm based on the discrete counterpart of the principle of maximum plastic dissipation. Develops the explicit expression for the consistent elasto‐plastic tangent modulus. All expressions are derived via tensor formulation showing the advantage over the classical matrix notation. The integration algorithm is implemented in the formulation of the four‐node isoparametric assumed‐strain finite‐rotation shell element employing the Mindlin‐Reissner‐type shell model. By applying the layered model, plastic zones can be displayed through the shell thickness. Material non‐linearity described by the von Mises yield criterion and isotropic hardening is combined with a geometrically non‐linear response assuming finite rotations. Numerical examples illustrate the efficiency of the present formulation in conjunction with the standard Newton iteration approach, in which no line search procedures are required. Demonstrates the excellent performance of the algorithm for large time respective load steps.

Numerical solutions in computational plasticity are severely challenged when concrete and geomaterials are considered with non‐regular yield surfaces, strain‐softening and non‐associated flow. There are two aspects that are of immediate concern within load steps which are truly finite: first, the iterative corrector must assure that the equilibrium stress state and the plastic process variables do satisfy multiple yield conditions with corners, Fi(σ, q) = 0, at discrete stages of the solution process. To this end, a reliable return mapping algorithm is required which minimizes the error of the plastic return step. Second, the solution of non‐linear equations of motion on the global structural level must account for limit points and premature bifurcation of the equilibrium path. The current paper is mainly concerned with the implicit integration of elasto‐plastic hardening/softening relations considering non‐associated flow and the presence of composite yield conditions with corners.

The purpose of this paper is to suggest a robust hybrid method for updating the stress and plastic internal variables in plasticity considering damage mechanics.

By benefiting the properties of the well-known explicit and implicit integrations, a new mixed method is derived. In fact, the advantages of the mentioned techniques are used to achieve an efficient integration.

The numerical studies demonstrate the high precision and robustness of the suggested algorithm.

The perfect von-Mises plasticity together with Lemaitre damage model is considered within the realm of small deformations.

Updating stress and plastic internal variables are of utmost importance in elastoplastic analyses of structures. The accuracy and efficiency of stress-updating methods significantly affect the final outcomes of nonlinear analyses.

The idea which is used to derive the hybrid method leads to an efficient integration method for updating the constitutive equations of the damage mechanics.

This paper aims at the problem of surrounding rock excavation damage zone of tunneling in the rich water region, this paper aims to propose a new seepage-stress-damage coupling model and studied the numerical algorithm. This reflects the characteristics of rock damage evolution, accompanied by plastic flow deformation and multi-field interaction.

First of all, rock elastoplastic damage constitutive model based on the Drucker–Prager criterion is established, the fully implicit return mapping algorithm is adopted to realize the numerical solution. Second, based on the relation between damage variation and permeability coefficient, the rock stress-seepage-damage model and multi-field coupling solving iterative method are presented. Finally, using the C++ language compiled the corresponding programs and simulated tunnel engineering in the rich water region.

Results show that difference evolution-based back analysis inversed damage parameters well, at the same time the established coupling model and calculating program have more advantages than general conventional methods. Multiple field coupling effects should be more considered for the design of tunnel support.

The proposed method provides an effective numerical simulation method for the construction of the tunnel and other geotechnical engineering involved underground water problems.