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.

Exact solutions for Mohr Coulomb elastoplasticity are developed. Using these solutions an exact stress increment for a given finite strain increment can be computed. The developed solutions are valid for perfect and linear hardening/softening plasticity using isotropic work hardening hypotheses. The solutions can be used to check computer codes and assess their ability to handle multiple active yield surfaces. Illustrative examples are included.

This paper aims to present a direct analysis to demonstrate why markedly different tensile and compressive behaviors of concretes could not be simulated with the Drucker–Prager yield criterion.

This study proposed an extended form of the latter for establishing a new elastoplasticity model with evolving yield strengths.

Explicit closed-form solutions to non-symmetric tensile and compressive responses of uniaxial specimens at finite strain are for the first time obtained from hardening to softening.

With such exact solutions, the yield strengths in tension and compression can be explicitly prescribed by uniaxial tensile and compressive stress-strain functions. Then, the latter two are further provided in explicit forms toward accurately simulating tensile and compressive behaviors. Numerical examples are supplied for meso-scale heterogeneous concrete (MSHC) and high-performance concrete (HPC), etc. Model predictions are in good agreement with test data.

This paper addresses the loading/unloading conditions of thediscrete initial—value problem of plastic flow at infinitesimal deformations. As in the continuum problem, it is established that the strain—space formulation of the loading conditions is primary. Generalized trapezoidal and mid‐point rules are discussed. The loading conditions established for the general non‐associated flow problem are shown to naturally reduce to well‐known inequalities for flow rules obeying normality.

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.

An implicit integration algorithm for elastoplastic constitutive equations in plane stress analysis is presented. The error associated with this algorithm is of the same order as the one reached in three‐dimensional analysis with the radial return algorithm. No subincrementation is needed. Moreover, the exact elastoplastic stress—strain matrix related to this algorithm is derived.

The purpose of this study is to show how gradual strength degradation of metal beams under cyclic bending up to fatigue failure is simulated based on a new elastoplasticity model free of any yield criterion.

A new approach is proposed toward accurately and explicitly prescribing evolution of non-uniform stress distribution on beam cross-section under cyclic bending and, as such, gradual degradation of the bending strength can be directly determined.

Explicit results for the bending response in a whole cyclic process up to fatigue failure are obtained and the fatigue characteristic curve is for the first time simulated directly between the curvature amplitude and the cycle number to failure.

First, explicit and accurate determination of the non-uniform stress distribution on beam cross-section is achieved with asymptotic softening effects. Second, degradation of the bending strength can be directly deduced cycle by cycle. Finally, the relationship between the bending moment and the curvature is calculated using new and efficient numerical algorithms, thus bypassing usual time-consuming calculations with finite element procedures. Numerical results are presented and in good agreement with experimental data.

This work addresses the computational aspects of a model for elastoplastic damage at finite strains. The model is a modification of a previously established model for large strain elastoplasticity described by Perić et al. which is here extended to include isotropic damage and kinematic hardening. Within the computational scheme, the constitutive equations are numerically integrated by an algorithm based on operator split methodology (elastic predictor—plastic corrector). The Newton—Raphson method is used to solve the discretized evolution equations in the plastic corrector stage. A numerical assessment of accuracy and stability of the integration algorithm is carried out based on iso‐error maps. To improve the stability of the local N—R scheme, the standard elastic predictor is replaced by improvedinitial estimates ensuring convergence for large increments. Several possibilities are explored and their effect on the stability of the N—R scheme is investigated. The finite element method is used in the approximation of the incremental equilibrium problem and the resulting equations are solved by the standard Newton—Raphson procedure. Two numerical examples are presented. The results are compared with those obtained by the original elastoplastic model.

Stress computation in finite element materially non‐linear analysis is an important problem that has perhaps been receiving less attention than it deserves. Not only does it consume a significant share of total computer time, but also inaccuracies and ‘savings’ thereupon may well jeopardize the gains aimed at by sophisticating elsewhere the numerical strategy. A well established algorithm for stress computation is reviewed in detail, illustrating a number of computational hazards and proposing simple solutions.

The purpose of this paper is to establish a stiffness model of fixed joint considering self-affinity and elastoplasticity of asperities.

The proposed model considers that asperities of different scales are interrelated rather than independent. For elastoplastic contact, a spring-damper model and an elastic deformation ratio function were proposed to calculate the contact stiffness of asperities.

A revised fractal asperity model was proposed to calculate the contact stiffness of fixed joint, the impacts of the fractal dimension, the fractal roughness parameter and the Meyer index on the contact stiffness were discussed, and the present experimental results and the Jiang’s experimental results showed that the stiffness can be well predicted by proposed model.

The contradiction between the Majumdar and Bhushan model and the Morag and Etsion model can be well explained by considering the interaction among asperities of different scales. For elastoplastic contact, elastic deformation ratio should be considered, and the stiffness of asperities increases first and then decreases with the increasing of interference.