Search results

The paper deals with the rheological‐dynamical analogy in which the three‐dimensional stress‐strain relations are defined under cyclic variation of stress for Hencky’s total strain theory. In many practical visco‐elasto‐plastic problems, like as multiaxial fatigue under loading at constant stress amplitude and constant stress ratio, the load‐carrying members are subjected to proportional loading. The classical Hencky’s theory has the advantage of mathematical convenience but its disadvantage is that the deformations predicted for the volume element are independent of the loading path. The existing formulations of the constitutive models for metals are mainly based on the Prandtl‐Reuss incremental theory of elasto‐plasticity, slip theory of plasticity or continuum damage mechanics. They have been shown capable of reproducing satisfactorily most experimental results available for metallic specimens. However, from the theoretical viewpoint little has been said about how these formulations relate to realistic predicting many different inelastic and time dependent problems of two‐ or threedimensional solids, such as fatigue, discontinuous plastic deformation etc. In this paper, fundamentally new aspect of isochronous constitutive relations for Hencky’s theory, which are dependent of the each loading path, is achieved by systematically introducing RDA concept into the continuum framework. Specific inelastic and fatigue formulation of triaxial state of stress is developed and discussed within the new theoretical tool and related to von Mises plasticity..

This paper is concerned with a new proposal regarding the analysis of visco‐elastoplasticity and fatigue and is based on rheological‐dynamical theory. Due to the analogy between rheological model and dynamical model with viscous damping, it becomes obvious that inelastic response of members is essentially a dynamical problem. An analytical rheological‐dynamical viscoelasto‐ plastic solution of one‐dimensional longitudinal continuous vibration under loading and solution for the stress relaxation as unloading have been developed and used to obtain the fatigue limit of thin long bars. Rheologic behavior of the bar can be characterized by one parameter, like in a single‐degree‐of‐freedom spring mass system. In all inelastic strains time rate effects are always present to some degree. Whether or not their exclusion has a significant influence on the prediction of the material fatigue behavior depends upon several factors like: maximum absolute stress in the cycle, coefficient of asymmetry of cycle, creep coefficient, slope of the strain hardening portion of the stress‐strain curve, relative frequency and uniaxial yield stress. This paper provides description of dynamic magnification factor, relaxation of stress, stress concentration and the fatigue limit of thin long symmetrical bars.

Various stress return algorithms in elastoplastic analyses using the finite element method require the evaluation of the contact (or penetration) stress state (Figure 1), defining the transition from elastic to elastoplastic behaviour. Various iterative schemes are commonly used to evaluate contact stress state with a great degree of precision, as subsequent analysis process (forward Euler, mid‐point rule stress return scheme) is greatly affected by the evaluation of the contact stress state, as has been stressed by several authors.

This paper presents a general procedure for solving 3D contact problems with implicit finite element codes. Emphasis is put on generality and robustness. Bodies in contact can be 3D solids or shells. Material and geometric nonlinearities can be dealt with (elasto‐plasticity, elasto‐visco‐plasticity, nonlinear elasticity, large displacements, strains and rotations). Different kinds of interaction are supported (tied, slip, friction). Advantage is taken of the solution history in order to improve the efficiency of the search algorithm. Numerical examples illustrate the general character of the proposed algorithm.

A variety of meshless methods have been developed in the last 20 years with an intention to solve practical engineering problems, but are limited to small academic problems due to associated high computational cost as compared to the standard finite element methods (FEM). The purpose of this paper is to develop an efficient and accurate algorithms based on meshless methods for the solution of problems involving both material and geometrical nonlinearities.

A parallel two-dimensional linear elastic computer code is presented for a maximum entropy basis functions based meshless method. The two-dimensional algorithm is subsequently extended to three-dimensional adaptive nonlinear and three-dimensional parallel nonlinear adaptively coupled finite element, meshless method cases. The Prandtl-Reuss constitutive model is used to model elasto-plasticity and total Lagrangian formulations are used to model finite deformation. Furthermore, Zienkiewicz and Zhu and Chung and Belytschko error estimation procedure are used in the FE and meshless regions of the problem domain, respectively. The message passing interface library and open-source software packages, METIS and MUltifrontal Massively Parallel Solver are used for the high performance computation.

Numerical examples are given to demonstrate the correct implementation and performance of the parallel algorithms. The agreement between the numerical and analytical results in the case of linear elastic example is excellent. For the nonlinear problems load-displacement curve are compared with the reference FEM and found in a very good agreement. As compared to the FEM, no volumetric locking was observed in the case of meshless method. Furthermore, it is shown that increasing the number of processors up to a given number improve the performance of parallel algorithms in term of simulation time, speedup and efficiency.

Problems involving both material and geometrical nonlinearities are of practical importance in many engineering applications, e.g. geomechanics, metal forming and biomechanics. A family of parallel algorithms has been developed in this paper for these problems using adaptively coupled finite element, meshless method (based on maximum entropy basis functions) for distributed memory computer architectures.

The purpose of this paper is to present eight local elasto‐plastic beam element formulations incorporated into the corotational framework for two‐noded three‐dimensional beams. These formulations capture the warping torsional effects of open cross‐sections and are suitable for the analysis of the nonlinear buckling and post‐buckling of thin‐walled frames with generic cross‐sections. The paper highlights the similarities and discrepancies between the different local element formulations. The primary goal of this study is to compare all the local element formulations in terms of accuracy, efficiency and CPU‐running time.

The definition of the corotational framework for a two‐noded three‐dimensional beam element is presented, based upon the works of Battini .The definitions of the local element kinematics and displacements shape functions are developed based on both Timoshenko and Bernoulli assumptions, and considering low‐order as well as higher‐order terms in the second‐order approximation of the Green‐Lagrange strains. Element forces interpolations and generalized stress resultant vectors are then presented for both mixed‐based Timoshenko and Bernoulli formulations. Subsequently, the local internal force vector and tangent stiffness matrix are derived using the principle of virtual work for displacement‐based elements and the two‐field Hellinger‐Reissner assumed stress variational principle for mixed‐based formulations, respectively. A full comparison and assessment of the different local element models are performed by means of several numerical examples.

In this study, it is shown that the higher order elements are more accurate than the low‐order ones, and that the use of the higher order mixed‐based Bernoulli element seems to require the least number of FEs to accurately model the structural behavior, and therefore allows some reduction of the CPU time compared to the other converged solutions; where a larger number of elements are needed to efficiently discretize the structure.

The paper reports computation times for each model in order to assess their relative efficiency. The effect of the numbers of Gauss points along the element length and within the cross‐section are also investigated.

A consistent formulation for unilateral contact problems including frictional work hardening or softening is proposed. The approach is based on an augmented Lagrangian approach coupled to an implicit quasi‐static Finite Element Method. Analogous to classical work hardening theory in elasto‐plasticity, the frictional work is chosen as the internal variable for formulating the evolution of the friction convex. In order to facilitate the implementation of a wide range of phenomenological models, the friction coefficient is defined in a parametrised form in terms of Bernstein polynomials. Numerical simulation of a 3D deep‐drawing operation demonstrates the performance of the methods for predicting frictional contact phenomena in the case of large sliding paths including high curvatures.

Liquefaction phenomenon and its catastrophic nature can be analysed as a particular material behaviour of granular media under certain loading paths. Proposes a definition of liquefaction and its modelling by constitutive relations. Discusses this modelling in relation to the questions of stability and uniqueness. Considers the signs of three scalar quantities: the work of second order, the determinant of the symmetric part of the tangent constitutive tensor and the determinant of the tensor itself. Concludes that the liquefaction path is situated inside a potentially unstable domain and that in some cases this path reaches some states of loss of uniqueness, which are essentially bifurcation points.

Focuses on a finite element computational model for the Timoshenko beam which is idealized as an elasto‐plastic‐damage medium governed by Lemaitre’s continuum damage mechanics (CDM) model for ductile fracture. Response under monotonically increasing loading does not show any deviation from elasto‐plastic simulation. However, a marked difference in the residual stress field is noted by virtue of the unloading phase, in which the CDM model allows for stiffness degradation in contrast to classical elasto‐plasticity which requires unloading at the (frozen) initial stiffness of the material.

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from pathological mesh‐dependence when strain‐softening models are employed in failure analyses. In this contribution the governing field equations are regularized by adding rotational degrees‐of‐freedom to the conventional translational degrees‐of‐freedom. This so‐called elasto‐plastic Cosserat continuum model, for which an efficient and accurate integration algorithm and a consistent tangent operator are also derived in this contribution, warrants convergence of the load—deflection curve to a unique solution upon mesh refinement and a finite width of the localization zone. This is demonstrated for an infinitely long shear layer and a biaxial test of a strain‐softening elasto‐plastic von Mises material.