SIMPLE SHAKEDOWN OF STRUCTURES UNDER VARIABLE MULTI-LOADINGS

The shakedown analysis of structures under variable multi-loadings is considered, and the corresponding simple shakedown condition is presented in this paper. Distribution of fixed stresses field is given, and the self-equilibrium of fixed stresses field is analyzed. Elastic shakedown and plastic shakedown conditions are presented based on the fixed stresses field. The theorem is convenient to evaluate the shakedown limit of structures under cyclical variable multiloadings through solving positive scalar fields and fixed stresses field factors at a series of dangerous positions of the structure, and tedious computations are avoided. Finally the theorem is applied to a thick-walled cylindrical tube under variable pressure and temperature, and the rolling contact problem. The results are in good agreement with some computational results.     

A modified basis reduction method for limited kinematic hardening shakedown analysis under complex loads

A novel extension of the basis reduction method for kinematic hardening shakedown problem is presented. Firstly, the basis reduction method is implemented based on the modified Newton–Raphson (N-R) method. Then a new technique for the construction of back stress field is introduced, where the simultaneous influence of multiple load corners in shakedown is taken into consideration. Finally, two typical numerical examples are investigated. The results compared with previous works in literatures demonstrated that the proposed method is accurate and the performance in reducing of computation time is significant.     

Lower bound shakedown analysis by the symmetric Galerkin boundary element method

In this paper, the static shakedown theorem is reformulated making use of the symmetric Galerkin boundary element method (SGBEM) rather than of finite element method. Based on the classical Melan’s theorem, a numerical solution procedure is presented for shakedown analysis of structures made of elastic-perfectly plastic material. The self-equilibrium stress field is constructed by linear combination of several basis self-equilibrium stress fields with parameters to be determined. These basis self-equilibrium stress fields are expressed as elastic responses of the body to imposed permanent strains obtained through elastic–plastic incremental analysis. The lower bound of shakedown load is obtained via a non-linear mathematical programming problem solved by the Complex method. Numerical examples show that it is feasible and efficient to solve the problems of shakedown analysis by using the SGBEM.     

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.     

Variational methods for the steady state response of elastic–plastic solids subjected to cyclic loads

Solids (or structures) of elastic–plastic internal variable material models and subjected to cyclic loads are considered. A minimum net resistant power theorem, direct consequence of the classical maximum intrinsic dissipation theorem of plasticity theory, is envisioned which describes the material behavior by determining the plastic flow mechanism (if any) corresponding to a given stress/hardening state. A maximum principle is provided which characterizes the optimal initial stress/hardening state of a cyclically loaded structure as the one such that the plastic strain and kinematic internal variable increments produced over a cycle are kinematically admissible. A steady cycle minimum principle, integrated form of the aforementioned minimum net resistant power theorem, is provided, which characterizes the structure’s steady state response (steady cycle) and proves to be an extension to the present context of known principles of perfect plasticity. The optimality equations of this minimum principle are studied and two particular cases are considered: (i) loads not exceeding the shakedown limit (so recovering known results of shakedown theory) and (ii) specimen under uniform cyclic stress (or strain). Criteria to assess the structure’s ratchet limit loads are given. These, together with some insensitivity features of the structure’s alternating plasticity state, provide the basis to the ratchet limit load analysis problem, for which solution procedures are discussed.     

Influence of geometrical nonlinearities on the shakedown of damaged structures

A generalization of Melan’s shakedown theorem is presented taking into account geometrical effects and plastic ductile damage. Numerical results illustrate the proposed method.     

Dynamic shakedown of structures with variable appended masses and subjected to repeated excitations

Elastic shakedown for discrete, or finite-element discretized, structures subjected to combinations of static and time-variable loads is addressed in the hypothesis of elastic-perfectly plastic material behavior. The static load is conceived as the weight of an additional mass appended to the structure, whereas the time-variable load is conceived as an unknown sequence of excitations belonging to a specified domain, with intervals between subsequent excitations during which the structure is considered as being motionless. It is shown that, in the plane of the static and time-variable load parameters, the structure's dynamic shakedown domain is nonconvex and that its boundary curve generally exhibits local minima and maxima at those static load values at which a resonant and anti-resonant structural behavior, respectively occurs. It is also shown that, for static loads close to the resonant behavior values, the shakedown limit load can be sensibly smaller than the value computed without taking the appended mass into account. The problem of evaluating the shakedown limit load is discussed and a numerical example presented.     

Shakedown analysis of engineering structures with limited kinematical hardening

We present a numerical method for the computation of shakedown loads of engineering structures with limited kinematical hardening under thermo-mechanical loading. The method is based on Melan’s statical shakedown theorem, which results in a nonlinear convex optimization problem. This is solved by an interior-point algorithm recently developed by the authors, specially designed for lower bound shakedown analysis of large-scale problems. Limited kinematical hardening is taken into account by use of a two-surface model, such that both alternating plasticity and incremental collapse can be captured. For the yield surface as well as for the bounding surface the von Mises criterion is used. The proposed method is validated by two examples, where numerical results are compared to those of literature where available.     

A method for the evaluation of a ratchet limit and the amplitude of plastic strain for bodies subjected to cyclic loading

An extension of the upper bound shakedown theorem to load histories in excess of shakedown has been presented elsewhere in this issue. Here the minimisation process described therein is applied to the solutions of the ratchet limit as well as shakedown and limit load for a range of simple problems. The solutions provide an estimate of the maxima of the varying plastic strain magnitudes, which is compared with the Neuber approximate values. The position of the ratchet boundary is confirmed by comparison step-by-step analysis.     

应变强化模型的安定准则研究

基于Melan经典的安定理论和von Mises屈服准则,建立了塑性应变强化条件下结构安定的数学模型,根据与时间无关的应力场的特性,对结构中与时间无关的应力场进行了合理的数学变换,将其与载荷变化系数联系起来,推导出与其对应的结构安定极限范围的表达式,给出塑性应变强化模型安定性存在的简化条件.该结论有利于简化应变强化条件下结构的安定分析.     

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

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

Structural shakedown for elastic–plastic materials with hardening saturation surface

An elastic–plastic material model with internal variables and thermodynamic potential, not admitting hardening states out of a saturation surface, is assumed as a basis to formulate a statical Melan-type shakedown theorem. Grounding on the optimality conditions relative to the shakedown load multiplier problem for a structure subjected to cyclic loads, the impending inadaptation collapse mechanism at the shakedown limit state is analyzed and discussed. It is shown that the adopted model is able to catch ratchetting collapse mode at a structural level. Numerical results for a simple structure are finally reported.     

A minimum theorem for cyclic load in excess of shakedown,with application to the evaluation of a ratchet limit

Although the shakedown theorems for perfect plasticity have been known since Koiter's 1960 review paper, extensions of the theory to situations where ratchetting or reverse plasticity occurs in excess of shakedown have not appeared in the literature. In this paper a generalisation of the upper bound theorem is derived which reduces to the upper bound shakedown theorem in the limiting case when the load point approaches the shakedown boundary. The new theory is used to develop a method for identifying the ratchet limit for a class of loading histories through the sequential minimisation of two functionals. A programming method, based on the Elastic Compensation method for shakedown is then derived and convergence proven. Numerical examples of the application of the method to practical problems are discussed by us in an accompanying paper.     

Shape optimization for elasto-plastic deformation under shakedown conditions

An integrated approach for all necessary variations within direct analysis, variational design sensitivity analysis and shakedown analysis based on Melan’s static shakedown theorem for linear unlimited kinematic hardening material behavior is formulated. Using an adequate formulation of the optimization problem of shakedown analysis the necessary variations of residuals, objectives and constraints can be derived easily. Subsequent discretizations w.r.t. displacements and geometry using e.g. isoparametric finite elements yield the well known ‘tangent stiffness matrix’ and ‘tangent sensitivity matrix’, as well as the corresponding matrices for the variation of the Lagrangian-functional which are discussed in detail. Remarks on the computer implementation and numerical examples show the efficiency of the proposed formulation. Important effects of shakedown conditions in shape optimization with elasto-plastic deformations are highlighted in a comparison with elastic and elasto-plastic material behavior and the necessity of applying shakedown conditions when optimizing structures with elasto-plastic deformations is concluded.     

A SOLUTION PROCEDURE FOR LOWER BOUND LIMIT AND SHAKEDOWN ANALYSIS BY SGBEM

The symmetric Galerkin boundary element method (SGBEM) instead of the finite element method is used to perform lower bound limit and shakedown analysis of structures. The self-equilibrium stress fields are constructed by a linear combination of several basic self-equilibrium stress fields with parameters to be determined. These basic self-equilibrium stress fields are expressed as elastic responses of the body to imposed permanent strains and obtained through elastic-plastic incremental analysis. The complex method is used to solve nonlinear programming and determine the maximal load amplifier. The limit analysis is treated as a special case of shakedown analysis in which only the proportional loading is considered. The numerical results show that SGBEM is efficient and accurate for solving limit and shakedown analysis problems. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctorial Dissertation of China (No. 200025) and the Basic Research Foundation of Tsinghua University.     

Inelastic responses of a two-bar system with temperature-dependent elastic modulus under cyclic thermomechanical loadings

This paper is concerned with the elastic plastic response of a two-bar system with temperature-dependent elastic coefficients under cyclic thermomechanical loadings. Such materials are characterized by lack of results concerning the asymptotic behaviors and conditions for shakedown occurrence. This study shows that the considered simple structure is sufficiently complex to experience different periodic long-term behaviors as in classical elastoplasticity. In order to understand how Melan–Koiter method works for such materials, the evolution of the structure’s response until the stabilization of the plastic strain (‘shakedown’) or the asymptotic dissipative behavior (‘alternating plasticity’ or ‘ratcheting’) is analytically addressed and the Bree diagram is then constructed. The main result of this work is that the residual stress and strain fields are time-dependent even when shakedown occurs. Besides, we proved that Halphen’s conjecture (Halphen, 2005) giving a sufficient condition for shakedown occurrence is not a necessary condition. Finally, numerical results performed by an incremental finite element procedure are presented.     

A nonlinear programming approach to kinematic shakedown analysis of frictional materials

This paper develops a novel nonlinear numerical method to perform shakedown analysis of structures subjected to variable loads by means of nonlinear programming techniques and the displacement-based finite element method. The analysis is based on a general yield function which can take the form of most soil yield criteria (e.g. the Mohr–Coulomb or Drucker–Prager criterion). Using an associated flow rule, a general yield criterion can be directly introduced into the kinematic theorem of shakedown analysis without linearization. The plastic dissipation power can then be expressed in terms of the kinematically admissible velocity and a nonlinear formulation is obtained. By means of nonlinear mathematical programming techniques and the finite element method, a numerical model for kinematic shakedown analysis is developed as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the plastic dissipation power which is to be minimized and an upper bound to the shakedown load can be calculated. An effective, direct iterative algorithm is then proposed to solve the resulting nonlinear programming problem. The calculation is based on the kinematically admissible velocity with one-step calculation of the elastic stress field. Only a small number of equality constraints are introduced and the computational effort is very modest. The effectiveness and efficiency of the proposed numerical method have been validated by several numerical examples.     

含缺陷轴对称体的安定与极限分析

利用静力安全定理得到了计算轴对称体安定与极限载荷的统一格式，采用温度参数法构造安定分析所需的残余应力场，为了克服对工程实际问题进行安定分析时解题规模与计算精度的矛盾，针对轴对称体的特点，采用两种线性化方案对屈服面进行线性化处理，即直接内接法和在降维应力偏量空间中对屈服面的线性化处理，使安定分析转化为一线性规划问题，在简化过程中合理选择线性化方案以便使应力校核点接近精确的屈服面；为了减小计算量，在求     

Interfacial slip and creep in rolling contact incorporating a cylinder with an elastic layer

A systematic approach for investigating the interfacial behaviour of tyred systems is proposed. A two-dimensional contact model of an elastic strip, shrink-fitted onto a wheel, and subjected to different rolling contact conditions, has been adopted to illustrate the method. The model combines existing techniques to explore individual elastic contact problems and it enables us to characterise the behaviour at the strip/substrate interface caused by loads induced by a quasi-static application of stationary and moving loads on the surface of the layer. The solution is compared to the stationary load case and regimes of local slip, full stick, separation and frictional creep are identified and collated for a variety of loading conditions, materials and geometries. Further, this article presents an investigation of frictional shakedown for layered systems subjected to periodic contact loading. The term shakedown is here referred to as the possibility of developing interfacial residual stresses at the layer/substrate interface such that frictional slip, originally activated by the applied external contact load, ceases after a few loading cycles. The possible applicability of the Melan’s theorem for elastic frictional system is investigated and preliminary results presented.