Bilateral bounds for structures under dynamic shakedown conditions

Bounds to functions of the plastic strains produced during the elastic-plastic stage of a dynamic loading process leading to shakedown can be constructed on the basis of the solution of a single, static, holonomic elastic-plastic problem. This bounding procedure is described in detail with respect to the case of periodic external loading, and the generalization to an arbitrary dynamic loading condition is provided. Simple examples illustrate the bounding effectiveness of the proposed technique.

---

sommario Si mostra come sia possibile produrre delimitazioni di funzioni delle deformazioni plastiche che avvengono durante la fase iniziale di un processo dinamico che si concluda con l'adattamento in regime elastico (shakedown). Si illustra in dettaglio la procedure di delimitazione con riferimento al caso di carichi esterni periodici, e si fornisce l'estensione concettuale al caso di una condizione di carico generica. Due esempi numerici mostrano l'accuratezza e l'efficienza numerica della tecnica proposta.

     

弹塑性结构安定下限分析的无网格局部

将基于Voronoi结构的无网格局部Petrov-Galerkin法与减缩基技术相结合,建立了一种安定下限分析的新方法．为了克服移动最小二乘近似难以准确施加本质边界条件的缺点,采用了自然邻近插值构造试函数.通过引入基准载荷域上载荷角点的概念,消除了安定下限分析中由时间参数所引起的求解困难．利用减缩基技术,将安定分析问题化为一系列未知变量较少的非线性规划子问题．在每个非线性规划子问题中,自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模拟,而这些自平衡应力场基矢量可应用弹塑性增量分析中的平衡迭代结果得到．算例结果证明了提出的分析方法的有效性．      

Initiation and propagation of localized deformation in elasto-plastic strips under uniaxial tension

This paper deals with the evolution of inhomogeneous deformation in shape memory alloy strips and mild steel strips under uniaxial tension. New experiments on NiTi strips, which initially are in an austenitic phase, show that at a critical stress level martensite nucleates in sharp bands inclined at 55° to the axis of loading. Under prescribed end displacement martensite subsequently spreads either by steady-state propagation of inclined transition fronts or via a criss-cross pattern of finger-like features. Similar events have been reported in the literature regarding the evolution of Lüders bands in fine grained steel strips and wires. The similarity of macroscopic events, despite the different mechanisms of instability at the micro-level, prompted us to approximate the material behavior as a finitely deforming elasto-plastic solid with a trilinear up-down-up nominal stress-strain response. Two such stress-strain responses were used in finite element simulations of strip tension tests. In the first the true stress-strain response maintains its stability and in the second the intermediate branch has a negative slope. While both material models produced inhomogeneous deformations with features similar to those of the experiments, the larger initiation peak associated with the second gave results which closely resembled specific experiments. The numerical simulations confirmed that the evolution of events seen in experiments on SMAs and mild steels is strongly influenced by overall geometric (structural) effects. Furthermore, the success of this simple continuum constitutive model strongly suggests that continuum level events remain dominant players in such fine grained materials.     

Endochronic analysis for finite elasto-plastic deformation and application to metal tube under torsion and metal rectangular block under biaxial compression

The ordinary differential constitutive equations of endochronic theory are extended to simulate elasto-plastic deformation in the range of finite strain using the concept of corotational rate. Different corotational stress rates (Jaumann, Cotter-Rivlin, Truesdell, Dienes and Mandel) are incorporated into the theory. In addition, a new formulation of the plastic spin, which can be used in the Mandel stress rate, is derived. Theoretical simulations of the axial effects for various materials subjected to simple and pure torsional loading cases are discussed in this study. It is shown that the endochronic theory incorporated with the Mandel stress rate yields the most satisfactory result, as indicated from comparison with the experimental data found in literature.

Finally, theoretical investigation of the deformation subjected to finite proportional and non-proportional biaxial compression is presented. The true relationship between stress and strain can be converted to a nominal stress-strain relationship for biaxial loading through the explicit transformation equations derived in this paper. Experimental data tested by Khan and Wang [1990] (“An Experimental study of Large finite Plastic Deformation in Annealed 1100 Aluminum During Proportional and Non-proportional Biaxial Compression” Int. J. Plasticity, 6, 485) are suitably described by the theory demonstrated from a comparison with the theoretical prediction according to rigid-plastic and elastic-plastic models employed by Huang and Khan [1991]. “An Analysis of Finite Elastic-Plastic Deformation under Biaxial Compression”, Int. J. Plasticity, 7, 219).     

Finite deformation crack-line fields in a thin elasto-plastic sheet

Finite deformation in the crack-tip zone of plastic deformation is investigated for Mode-I opening of a crack in a thin sheet of elasto-plastic material. The material obeys the von Mises yield criterion in the true stresses, and the stretching tensor satisfies a flow law of the Prandtl-Reuss type. Incompressibility and a state of generalized plane stress are assumed. It is assumed that linearized elasticity applies outside the zone of plastic deformation. On the crack-line between the crack-tip and the elastic—plastic boundary, two distinct regions have been recognized: the near-tip zone and the intermediate region. In the near-tip zone the fields are controlled by the radius of curvature of the blunted crack-tip. Here the stress field has been approximated by classical plane stress results. It has been assumed that the crack-line stresses may be taken as uniform in the intermediate region. In each region, deformation variables have been determined by the use of the constitutive relations, and the results have been matched to the corresponding quantities in the neighboring region(s). In this manner expressions have been constructed for the deformation gradients on the crack-line, in terms of the distance to the crack-tip in the deformed configuration, the yield stress in shear, and the stress intensity factor of linear elastic fracture mechanics.     

Bounds for shakedown of cohesive-frictional materials under moving surface loads

In this paper, shakedown of a cohesive-frictional half space subjected to moving surface loads is investigated using Melan’s static shakedown theorem. The material in the half space is modelled as a Mohr–Coulomb medium. The sliding and rolling contact between a roller and the half space is assumed to be plane strain and can be approximated by a trapezoidal as well as a Hertzian load distribution. A closed form solution to the elastic stress field for the trapezoidal contact is derived, and is then used for the shakedown analysis. It is demonstrated that, by relaxing either the equilibrium or the yield constraints (or both) on the residual stress field, the shakedown analysis leads to various bounds for the elastic shakedown limit. The differences among the various shakedown load factors are quantitatively compared, and the influence of both Hertzian and trapezoidal contacts for the half space under moving surface loads is studied. The various bounds and shakedown limits obtained in the paper serve as useful benchmarks for future numerical shakedown analysis, and also provide a valuable reference for the safe design of pavements.     

A return-map algorithm for general associative isotropic elasto-plastic materials in large deformation regimes

The present paper addresses a flexible solution algorithm for associative isotropic elasto-plastic materials, i.e. for materials whose elastic and plastic behaviors are described through an isotropic free-energy function, an isotropic yield function and an associative flow rule. The discussion is relative to a large deformation regime, while no hardening mechanisms are included. The algorithm is based on a combination of the operator split method and a return map scheme. Both the algorithm linearization and the requirements for the yield criterion convexity are discussed in detail. Finally, to show the algorithm flexibility and performance, the discussion is specialized to three yield criteria and some test problems are studied.     

Thermocouple measurement of metal temperature under pulsed deformation conditions

Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 137–143, January–February, 1991.     

Necking of elasto-plastic rods under tension

This work is concerned with the modeling of the necking phenomenon in metallic tension members using a non-linear theory of elasto-plastic rods with deformable cross section. The study of the coupling between the axial deformation and the cross section deformation in tensile specimens is a basic step to understand the strain-softening and localization behaviors in more complex metallic structures. The main features of the model are illustrated through the simulation of the necking of an ASTM 6351 aluminium bar under tension.     

Limit and shakedown analysis under hydrogen embrittlement condition

A modified shakedown theorem and its solving technique are presented to involve hydrogen embrittlement of steel into limit and shakedown analysis. Firstly, the shakedown theorem for hydrogen embrittled material is derived from a limited kinematic hardening shakedown theorem and hydrogen enhanced localized plasticity mechanism of hydrogen embrittlement. In the presented theorem, hydrogen’s effect is taken into account by the synergistic action of both strength reduction and stress redistribution. Secondly, a novel solving technique is developed based on the basis reduction method, in which the complicated constraints in the resulting nonlinear mathematical programming are released. At last, three numerical examples are carried out to verify the performance of the proposed method and to reveal hydrogen’s effect on the limit and shakedown load of structure. The numerical results are discussed and compared with those from literatures, which proves the accuracy and high efficiency of the introduced solving technique. It is concluded that the proposed theorem can predict the limit and shakedown load of hydrogen embrittled structure reasonably.     

Shape optimization for a composite crack-length sensor

Conductive polymeric sensors can be used as crack-length gages. They have the advantages of being inexpensive, versatile in size and require inexpensive instrumentation. They also measure crack length continuously. The shape of these gages can be varied in several ways to improve their sensitivity. The gages with the 'optimized' or tapered shape overcome the principal disadvantages of rectangular gages, namely poor sensitivity at small crack lengths.     

Procedure for describing the deformation of materials to fracture under conditions of near-superplasticity

Novosibirsk. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 1, pp. 149–157, January–February, 1995.     

Computational formulation of shakedown analysis as a conic quadratic optimization problem

In the analysis of shakedown problems that have been discretized by means of the finite element method, large and sparse optimization problems are generated. The purpose of this paper is to provide the details of how such a problem can be cast in the form of a conic quadratic optimization problem, making use of Melan’s static theorem. An effective algorithm for the solution of the optimization problem is proposed.     

Analyses on micro-fracture of crack tip under large deformation conditions

In the present paper, Gurson's constitutive equation, which takes into account the development of voids, is used to study the behaviour of the material in the region near crack tip. Furthermore, the effect of void development on Young's modulus, which was not considered by Gurson, is taken into consideration. The analyses on void development, on stress distribution near crack tip, and on the variance of COD for the plane strain mode I problem are carried out with the large elastic-plastic deformation finite element method. The results are compared with those estimated from the Prandtl-Reuss constitutive equation.