Kinematic limit analysis of frictional materials using nonlinear programming

In this paper, a nonlinear numerical technique is developed to calculate the plastic limit loads and failure modes of frictional materials by means of mathematical programming, limit analysis and the conventional displacement-based finite element method. The analysis is based on a general yield function which can take the form of the Mohr–Coulomb or Drucker–Prager criterion. By using an associated flow rule, a general nonlinear yield criterion can be directly introduced into the kinematic theorem of limit analysis without linearization. The plastic dissipation power can then be expressed in terms of kinematically admissible velocity fields and a nonlinear optimization formulation is obtained. The nonlinear formulation only has one constraint and requires considerably less computational effort than a linear programming formulation. The calculation is based entirely on kinematically admissible velocities without calculation of the stress field. The finite element formulation of kinematic limit analysis is developed and solved as a nonlinear mathematical programming problem subject to a single equality constraint. The objective function corresponds to the plastic dissipation power which is then minimized to give an upper bound to the true limit load. An effective, direct iterative algorithm for kinematic limit analysis is proposed in this paper to solve the resulting nonlinear mathematical programming problem. The effectiveness and efficiency of the proposed method have been illustrated through a number of numerical examples.     

Limit analysis of composite materials based on an ellipsoid yield criterion

Employing repeating unit cell (RUC) to represent the microstructure of periodic composite materials, this paper develops a numerical technique to calculate the plastic limit loads and failure modes of composites by means of homogenization technique and limit analysis in conjunction with the displacement-based finite element method. With the aid of homogenization theory, the classical kinematic limit theorem is generalized to incorporate the microstructure of composites. Using an associated flow rule, the plastic dissipation power for an ellipsoid yield criterion is expressed in terms of the kinematically admissible velocity. Based on nonlinear mathematical programming techniques, the finite element modelling of kinematic limit analysis is then 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 limit load of a composite is then obtained. The nonlinear formulation has a very small number of constraints and requires much less computational effort than a linear formulation. An effective, direct iterative algorithm is proposed to solve the resulting nonlinear programming problem. The effectiveness and efficiency of the proposed method have been validated by several numerical examples. The proposed method can provide theoretical foundation and serve as a powerful numerical tool for the engineering design of composite materials.     

??????????????????????

将有限元方法引入到塑性极限分析中,采用刚体有限元离散挡土墙后土体计算区域, 同时构造运动许可速度场,在满足屈服条件、流动法则、虚功方程以及相应的边界条件的基 础上,建立约束方程,引入数学规划方法求解挡土墙在不同变位模式下极限土压力分布. 算 例说明了该方法的正确性和有效性.     

On kinematic method in shakedown theory: I. Duality of extremum problems

A kinematic method for determining the safety factor in shakedown problems is developed. An upper bound kinematic functional is defined on a set of kinematically admissible time-independent velocity fields. Every value of the functional is an upper bound for the safety factor. Using convex analysis methods, conditions are established under which the infimum of the kinematic upper bounds equals the safety factor, in particular, conditions under which it is sufficient to consider only smooth velocity fields for the safety factor calculation. The method generalizes that recently proposed for the case of spherical yield surfaces by Kamenjarzh and Weichert. The extension covers a wide class of yield surfaces and inhomogeneous bodies. A shakedown problem for a beam subjected to a concentrated load is considered as an example.     

STRENGTH PROPERTIES OF JOINTED ROCK MASSES BASED ON THE HOMOGENIZATION METHOD

This paper aims to determine the strength properties of jointed rock masses by means of the homogenization method.To reflect the microstructure of jointed rock masses,a representative element volume (R...     

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.     

Shakedown analysis of shell structures of kinematic hardening materials

It is of great practical importance to analyze the shakedown of shell structures under cyclic loading, especially of those made of strain hardening materials.In this paper, some further understanding of the shakedown theorem for kinematic hardening materials has been made, and it is applied to analyze the shakedown of shell structures. Though the residual stress of a real state is related to plastic strain, the time-independent residual stress field as we will show in the theorem may be unrelated to the time-independent kinematically admissible plastic strain field . For the engineering application, it will be much more convenient to point this out clearly and definitely, otherwise it will be very difficult. Also we have proposed a new method of proving this theorem. The above theorem is applied to the shakedown analysis of a cylindrical shell with hemispherical ends. According to the elastic solution, various possible residual stress and plastic strain fields, the shakedown analysis of the structure can be reduced to a mathematical programming problem. The results of calculation show that the shakedown load of strain hardening materials is about 30–40% higher than that of ideal plastic materials. So it is very important to consider the hardening of materials in the shakedown analysis, for it can greatly increase the structure design capacity, and meanwhile provide a scientific basis to improve the design of shell structures.     

多孔材料塑性极限载荷及其破坏模式分析

运用塑性力学中的机动极限分析理论，研究韧性基体多孔材料的塑性极限承载能力和破坏模式。以多孔材料的细观结构为研究对象，将细观力学中的均匀化理论引入到塑性极限分析中，并结合有限元技术，建立细观结构极限载荷的一般计算格式，并提出相应的求解算法。数值算例表明：细观孔洞对材料的宏观强度影响明显；在单向拉伸作用下，孔洞呈现膨胀扩大规律；多孔材料破坏源于基体塑性区的贯通。     

基于均匀化理论韧性复合材料塑性极限分析

运用细观力学中的均匀化方法分析了韧性复合材料的塑性极限承载能力.从反映复合材料细观结构的代表性胞元入手,将均匀化理论运用到塑性极限分析中,计算由理想刚塑性、Mises组分材料构成的复合材料的极限承载能力.运用机动极限方法和有限元技术,最终将上述问题归结为求解一组带等式约束的非线性数学规划问题,并采用一种无搜索直接迭代算法求解.为复合材料的强度分析提供了一个有效手段.     

On shakedown of three-dimensional elastoplastic strain-hardening structures

The present article considers the shakedown problem of structures made of either kinematic or mixed strain-hardening materials. Some basic and useful shakedown properties of elastoplastic strain-hardening structures are proved mathematically. It is impossible for a kinematic strain-hardening structure to be involved in incremental plastic collapse, and so its only possible failure mode is that of alternating plasticity. A time-independent self-equilibrium stress field has no influence on the shakedown of a kinematic strain-hardening structure although it contributes to the magnitude of plastic deformation. The sufficient shakedown conditions for either kinematic or mixed strain-hardening structures are deduced, from which the lower bound of shakedown load domain can be obtained via a mathematical programming problem. It should be pointed out that, to guarantee the safety of an elastoplastic strain-hardening structure, the damage analysis is also necessary to determine the maximum load factor the structure can bear. The shakedown analysis of strain-hardening structures can be simplified by the conclusions obtained in this article, as is illustrated by two simple examples.     

Performance of the MLPG method for static shakedown analysis for bounded kinematic hardening structures

Shakedown analysis is an extension of plastic limit analysis to the case of variable repeated loads and plays a significant role in safety assessment and structural design. This paper presents a solution procedure based on the meshless local Petrov–Galerkin (MLPG) method for lower-bound shakedown analysis of bounded kinematic hardening structures. The numerical implementation is very simple and convenient because it is only necessary to construct an array of nodes in the targeted domain. Moreover, the natural neighbour interpolation (NNI) is employed to construct trial functions for simplifying the imposition of essential boundary conditions. The kinematic hardening behaviour is simulated by an overlay model and the numerical difficulties caused by the time parameter are overcome by introducing the conception of load corner. The reduced-basis technique is applied to solve the mathematical programming iteratively through a sequence of reduced residual stress subspaces with very low dimensions and the resulting non-linear programming sub-problems are solved via the Complex method. Numerical examples demonstrate that the proposed solution procedure is feasible and effective to determine the shakedown loads of bounded kinematic hardening structures as well as unbounded kinematic hardening structures.     

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.     

Shakedown reduced kinematic formulation,separated collapse modes,and numerical implementation

Our shakedown reduced kinematic formulation is developed to solve some typical plane stress problems, using finite element method. Whenever the comparisons are available, our results agree with the available ones in the literature. The advantage of our approach is its simplicity, computational effectiveness, and the separation of collapse modes for possible different treatments. Second-order cone programming developed for kinematic plastic limit analysis is effectively implemented to study the incremental plasticity collapse mode. The approach is ready to be used to solve general shakedown problems, including those for elastic–plastic kinematic hardening materials and under dynamic loading.     

Elastoplastic Analysis with Limited Plastic Deformations and Displacements*

ABSTRACT

This paper presents solution methods for elastoplastic and shakedown analysis of linearly elastic, perfectly plastic bodies for which the conventional classical formulations of these problems are completed by constraints on overall plastic deformation and elastoplastic displacement. The methods are described in terms of nonlinear mathematical programming and provide solutions when the plastic reserves of the body are not fully exhausted, and the plastic performance and the plastic deformations are controlled. Application of the method is illustrated by an example.     

Theoretical and computational aspects in the shakedown analysis of finite elastoplasticity

A fully nonlinear shakedown analysis is considered for structures undergoing large elastic-plastic strains. The underlying kinematics of finite elastoplasticity are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. It is Shown that the notion of a fictitious, self-equilibrated residual stress field of Melan's linear shakesdown theorem has to be replaced by the notion of real, self-equilibrated residual state. Path-dependent and path-independent shakedown theorems are presented that can be realized in an incremental step-by-step procedure using Finite Element codes. The numerical implementation is considered for highly nonlinear truss structures undergoing large cyclic deformations with ideal-plastic, isotropic and kinematic hardening material behavior. Path-dependency of the residual states in the case of non-adaptation and path-independency in the case of shakedown are shown, and the shakedown domain is determined taking into account also the stability boundaries of the structure.     

Shakedown kinematic theorem for elastic–perfectly plastic bodies

The classical shakedown kinematic theorem due to Koiter for elastic–perfectly plastic bodies is re-examined and divided into separated shakedown and nonshakedown theorems. While the shakedown theorem is based on the set of Koiter's plastic strain rate cycles, the non-shakedown one involves a broader set of admissible plastic strain rate cycles, the end-cycle accumulated strains of which are deviatoric parts of compatible strain fields. For certain broad classes of practical problems the two statements are unified to yield the unique theorem in Koiter's sense.     

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.     

Reduction of Kinematic: Inequality and Optimization of Perfectly Plastic Frameworks

Abstract

In the kinematic theory of structures consisting of perfectly plastic elements, an inequality between the plastic dissipation work and the load work is used. This inequality, which we will term “the kinematic inequality,” must hold for all kinematically admissible mechanisms. These mechanisms are generated by certain parameters which usually remain in the kinematic inequality and which thereby preclude the general application of the kinematic approach. In this paper we overcome this difficulty in the case of frames and provide various applications of the method. By using new theorems we eliminate the parameters and reduce the kinematic inequality to a finite system of inequalities which depend only on frame geometry and on loads. Based on these theorems, a procedure is offered for deriving a system of independent inequalities for general multistory multibay frames. New theorems are then obtained regarding the existence and the rotation of certain plastic hinges in collapse mechanisms. The overall theory is illustrated by a specific example. Finally, the formulations obtained following our method are used to minimize the mass of a fixed-base rectangular portal frame for any length, height, and system of loads.     

Micro/macromechanical plastic limit analyses of composite materials and structures

The load-bearing capacities of ductile composite materials and structures are studied by means of a combined micro/macromechanics approach. Firstly, on the microscopic scale, the aim is to get the macroscopic strength domains by means of the homogenization theory of micromechanics. A representative volume element (RVE) is selected to reflect the microstructures of the composite materials. By introducing the homogenization theory into the kinematic limit theorem of plastic limit analysis, an optimization format to directly calculate the limit loads of the RVE is obtained. And the macroscopic yield criterion can be determined according to the relation between macroscopic and microscopic fields. Secondly, on the macroscopic scale, by introducing the Hill's yield criterion into the kinematic limit theorem, the limit loads of orthotropic structures such as unidirectional fiber-reinforced composite structures are worked out. The finite element modeling of the kinematic limit analysis is deduced into a nonlinear mathematical programming with equality-constraint conditions that can be solved by means of a direct iterative algorithm. Finally, some examples are illustrated to show the application of the present approach. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctoral Dissertation of China (No. 200025), the Fund of the Ministry of Education of China for Returned Oversea Scholars and the Basic Research Foundation of Tsinghua University.     

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

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