On the existence of indeterminate solutions to the equations of motion under non-associated flow

Under certain conditions, an indeterminate solution exists to the equations of motion for dynamic elastic–plastic deformation of materials using constitutive laws based on non-associated flow that suggests that an initially unbounded dynamic perturbation in the stress can develop from a quiescent state on the yield surface. The existence of this indeterminate solution has been alleged to discourage use of non-associated flow rules for both dynamic and quasi-static analysis theoretically. It is shown in this paper that the indeterminate solution that may solve the equations of motion is intrinsically dynamic, and it determinately goes to zero in the quasi-static limit regardless of other indeterminate parameters. Consequently, the existence of this unstable dynamic solution has no impact on stability and use of non-associated flow rules for analysis of the quasi-static problem. More importantly, for dynamic applications, it is also shown that the indeterminate solution solves the equations of motion only if critical restrictions are applied to the constitutive equations such that the effective modulus during loading is constant and the direction of the perturbation is unidirectional over a finite time interval. It is shown that common components of the constitutive laws used in metal forming and deformation analysis are inconsistent with these restrictions. So, these common models can be generalized to include non-associated flow for analysis of the dynamic problem without concern that the solution will become indeterminate.     

A bipotential-based limit analysis and homogenization of ductile porous materials with non-associated Drucker–Prager matrix

In Gurson's footsteps, different authors have proposed macroscopic plastic models for porous solid with pressure-sensitive dilatant matrix obeying the normality law (associated materials). The main objective of the present paper is to extend this class of models to porous materials in the context of non-associated plasticity. This is the case of Drucker–Prager matrix for which the dilatancy angle is different from the friction one, and classical limit analysis theory cannot be applied. For such materials, the second last author has proposed a relevant modeling approach based on the concept of bipotential, a function of both dual variables, the plastic strain rate and stress tensors. On this ground, after recalling the basic elements of the Drucker–Prager model, we present the corresponding variational principles and the extended limit analysis theorems. Then, we formulate a new variational approach for the homogenization of porous materials with a non-associated matrix. This is implemented by considering the hollow sphere model with a non-associated Drucker–Prager matrix. The proposed procedure delivers a closed-form expression of the macroscopic bifunctional from which the criterion and a non-associated flow rule are readily obtained for the porous material. It is shown that these general results recover several available models as particular cases. Finally, the established results are assessed and validated by comparing their predictions to those obtained from finite element computations carried out on a cell representing the considered class of materials.     

Solution existence conditions for elastoplastic constitutive models of granular materials

In this paper, the conditions of solution existence for stress rates under given strain rates are investigated. The focus of the solution existence investigation is on the non-associated flow rule and elastic stress–strain relationship. Granular materials characterized with strong non-associated plastic flows are used as a particular example for analysis. Various flow rules for granular materials are analyzed, including Rowe’s, Roscoe’s flow rules and their modified versions. In the elastic stress–strain relationships of materials, the effects of Poisson’s ratio on solution existence are investigated. Both isotropic and anisotropic elasticity are considered. Given a granular material and its states, it is found that there exists a critical Poisson’s ratio for a particular non-associated flow rule. When the Poisson’s ratio of a material is above this critical Poisson’s ratio, its constitutive model is susceptible to solution non-existence. It is suggested that special attentions should be paid to the selection of material Poisson’s ratio and non-associated flow rule to ensure the existence of elastoplastic solutions.     

Review of Drucker’s postulate and the issue of plastic stability in metal forming

Drucker’s postulate defines a class of stable work hardening materials that are classified as non-energetic and is equivalent to the associated flow rule (AFR). The postulate has been shown to be a sufficient condition for plastic stability. However, experiments indicate that plastic deformation of aluminum and steel alloys does not adhere to the constraints of the AFR. Therefore, the requirement for accuracy suggests that the metal forming industry should also consider material models that are based on non-associated flow. But Drucker’s work raises the issue of stability when considering the use of non-associated flow in material models. While this concern is merited and many types of instability arises from certain types of non-associated flow, this has led to a widely accepted view that Drucker’s postulate is a necessary condition for stability. This perception is inhibiting the acceptance or consideration of more accurate material models that are suggested from the experimental observations about violations of the AFR. This paper proposes a specific class of material models based on non-associated flow and derives the constraints on this class of models to ensure stability. The existence of this class of non-AFR models proves that Drucker’s postulate is a sufficient but not necessary condition for stability. Furthermore, the class of models described in this paper is quite general and provides a framework for consideration of potentially more accurate material models while guaranteeing the same level of stability as typically associated with materials that satisfy Drucker’s postulate.     

Three-dimensional limit analysis of rigid blocks assemblages. Part II: Load-path following solution procedure and validation

A novel solution procedure for the non-associated limit analysis of rigid blocks assemblages is proposed. This proposal produces better solutions than previously proposed procedures and it is also able to provide an insight into the structural behaviour prior to failure. The limit analysis model proposed in Part I of this paper and the solution procedure are validated through illustrative examples in three-dimensional masonry piers and walls. The use of limit analysis for three-dimensional problems incorporating non-associated flow rules and a coupled yield surface is novel in the literature.     

Plastic collapse in non-associated hardening materials with application to Cam-clay

This paper presents a variational formulation for the analysis of plastic collapse conditions for a class of hardening materials that accounts for some non-associated flow laws such as the modified Cam-clay model of soils. In this framework, classical statical and kinematical principles of limit analysis do not hold. The variational principle is formulated for the general class of materials whose flow equations are derived from a kind of generalized potentials named bipotentials by de Saxcé.The plastic collapse phenomenon for hardening materials is considered first and formulated as a system of equations. In particular, the case of the usual modified Cam-clay model is analyzed. The paper follows with the proposal of a minimization principle whose solution is then related to the solution of the plastic collapse problem. We demonstrate the use of this minimum principle in a simple example of triaxial compression of a modified Cam-clay material. Finally, we discuss the particular form of the proposed variational formulation for the case of associated plasticity.     

A micromechanics-based elastoplastic damage model for granular materials at low confining pressure

This paper is devoted to the formulation of a micromechanics-based constitutive model for granular materials under relatively low confining pressure. The constitutive formulation is performed within the general framework of homogenization for granular materials. However, new rigorous stress localization laws are proposed. Some local constitutive relations are established under the consideration of irreversible thermodynamics. Macroscopic plastic deformation is obtained by considering local plastic sliding in a limit number of families of contact planes. The plastic sliding at each contact plane is described by a non-associated plastic flow rule, taking into account pressure sensitivity and normal dilatancy. Nonlinear elastic deformation related to progressive compaction of contacts is also taken into account. Material softening is described by involving damage process related to degradation of microstructure fabric. The proposed model is applied to some typical granular materials (sands). The numerical predictions are compared with experimental data.     

A plasticity model for powder compaction processes incorporating particle deformation and rearrangement

This paper develops a mechanistic model of granular materials that can be used with a commercial finite element package (ABAQUS). The model draws on the ideas of critical state soil mechanics and combines them with the theory of envelopes to develop an elasto-plastic model with a non-associated flow rule. The model incorporates both local deformation at the granule contacts, and rearrangement of the granules so that jointly they account for any bulk deformation. The mechanics of the model closely reflect the physicality of the material behaviour and the model parameters are closely linked (although not simplistically identical) to the characteristics of the granules. This not only gives an insight into the material behaviour, but also enables the model to be used to facilitate design of the material, its processing properties and, hence, component development. The model is used to simulate drained triaxial tests, settlement of a powder in a bin, and some examples of die pressing. Simulations are compared with experimental data and with predictions obtained using other models.     

Nonlocal continuum theory of anisotropically damaged metals

The paper deals with a consistent and systematic general framework for the development of anisotropic continuum damage in ductile metals based on thermodynamic laws and nonlocal theories. The proposed model relies on finite strain kinematics based on the consideration of damaged as well as fictitious undamaged configurations related via metric transformation tensors which allow for the interpretation of damage tensors. The formulation is accomplished by rate-independent plasticity using a nonlocal yield condition of Drucker–Prager type, anisotropic damage based on a nonlocal damage growth criterion as well as non-associated flow and damage rules. The nonlocal theory of inelastic continua is established to be able to take into account long-range microstructural interaction. The approach incorporates macroscopic interstate variables and their higher-order gradients which properly describe the change in the internal structure and investigate the size effect of statistical inhomogeneity of the heterogeneous material. The idea of bridging length-scales is made by using higher-order gradients in the evolution equations of the equivalent inelastic strain measures which leads to a system of elliptic partial differential equations which is solved using the finite difference method at each iteration of the loading step and the displacement-based finite element procedure is governed by the standard principle of virtual work. Numerical simulations of the elastic–plastic deformation behavior of damaged solids demonstrate the efficiency of the formulation. Tension tests undergoing large strains are used to investigate the damage growth in high strength steel. The influence of various model parameters on the prediction of the deformation and localization of ductile metals is discussed.     

Micromechanical approach to the strength properties of frictional geomaterials

The present paper describes a micromechanics-based approach to the strength properties of composite materials with a Drucker–Prager matrix in the situation of non-associated plasticity. The concept of limit stress states for such materials is first extended to the context of homogenization. It is shown that the macroscopic limit stress states can theoretically be obtained from the solution to a sequence of viscoplastic problems stated on the representative elementary volume. The strategy of resolution implements a non-linear homogenization technique based on the modified secant method. This procedure is applied to the determination of the macroscopic strength properties and plastic flow rule of materials reinforced by rigid inclusions, as well as for porous media. The role of the matrix dilatancy coefficient is in particular discussed in both cases. Finally, finite element solutions are derived for a porous medium and compared to the micromechanical predictions.     

Anisotropic parameter identification using inhomogeneous tensile test

In this contribution, an inverse identification strategy of constitutive laws for elastoplastic behaviour is presented. The proposed inverse algorithm is composed on an appropriate finite element calculation combined with an optimisation procedure. It is applied to identify material anisotropic coefficients using a set up of easy performed laboratory tests. The used experimental data are the plane tensile test and the off axes tensile tests. The identified behaviour models are mainly based on Hill's quadratic yield criterion. Two cases of this yield criterion have been considered: the transverse isotropic and the orthotropic one under an associated and non-associated flow rule assumptions for each case. The yield surface has been assumed to expand isotropically (isotropic strain hardening law) as a function of the plastic work.In order to better describe anisotropic plastic properties of the studied materials, a recently planar anisotropic yield function is used. It is a non-quadratic yield criterion which takes account of anisotropic yield stresses as well as anisotropic strain ratios. It is subsequently shown that the agreement between inverse identification results and experimental measurements were improved.We prove also that the presented strategy is a good alternative to the simplified homogeneous tests assumption, especially for the plane tensile test.     

UH模型切线刚度矩阵对称化及其应用

分析了目前常用的几种土的本构模型三维化方法以及变换应力方法的优势. 详细分析了用变换应力三维化本构模型非关联流动的本质, 针对因采用非关联流动法则而带来的本构模型切线刚度矩阵的非对称性, 建立了本构模型切线刚度矩阵的对称化方法. 将此方法应用于UH模型并进行了有限元分析, 进行了三轴伸长试验有限元分析. 通过以上分析比较发现, 切线刚度矩阵对称化方法能够明显改善有限元分析的收敛性, 缩短分析时间, 有着非常好的应用前景.      

Evaluation of advanced anisotropic models with mixed hardening for general associated and non-associated flow metal plasticity

The main objective of this paper is to develop a generalized finite element formulation of stress integration method for non-quadratic yield functions and potentials with mixed nonlinear hardening under non-associated flow rule. Different approaches to analyze the anisotropic behavior of sheet materials were compared in this paper. The first model was based on a non-associated formulation with both quadratic yield and potential functions in the form of Hill’s (1948). The anisotropy coefficients in the yield and potential functions were determined from the yield stresses and r-values in different orientations, respectively. The second model was an associated non-quadratic model (Yld2000-2d) proposed by Barlat et al. (2003). The anisotropy in this model was introduced by using two linear transformations on the stress tensor. The third model was a non-quadratic non-associated model in which the yield function was defined based on Yld91 proposed by Barlat et al. (1991) and the potential function was defined based on Yld89 proposed by Barlat and Lian (1989). Anisotropy coefficients of Yld91 and Yld89 functions were determined by yield stresses and r-values, respectively. The formulations for the three models were derived for the mixed isotropic-nonlinear kinematic hardening framework that is more suitable for cyclic loadings (though it can easily be derived for pure isotropic hardening). After developing a general non-associated mixed hardening numerical stress integration algorithm based on backward-Euler method, all models were implemented in the commercial finite element code ABAQUS as user-defined material subroutines. Different sheet metal forming simulations were performed with these anisotropic models: cup drawing processes and springback of channel draw processes with different drawbead penetrations. The earing profiles and the springback results obtained from simulations with the three different models were compared with experimental results, while the computational costs were compared. Also, in-plane cyclic tension–compression tests for the extraction of the mixed hardening parameters used in the springback simulations were performed for two sheet materials.     

A finite element formulation based on non-associated plasticity for sheet metal forming

In the present paper, a finite element formulation based on non-associated plasticity is developed. In the constitutive formulation, isotropic hardening is assumed and an evolution equation for the hardening parameter consistent with the principle of plastic work equivalence is introduced. The yield function and plastic potential function are considered as two different functions with functional form as the yield function of Hill [Hill, R., 1948. Theory of yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. A 193, 281–297] or Karafillis–Boyce associated model [Karafillis, A.P. Boyce, M., 1993. A general anisotropic yield criterion using bounds and a transformation weighting tensor. J. Mech. Phys. Solids 41, 1859–1886]. Algorithmic formulations of constitutive models that utilize associated or non-associated flow rule coupled with Hill or Karafillis–Boyce stress functions are derived by application of implicit return mapping procedure. Capabilities in predicting planar anisotropy of the Hill and Karafillis–Boyce stress functions are investigated considering material data of Al2008-T4 and Al2090-T3 sheet samples. The accuracy of the derived stress integration procedures is investigated by calculating iso-error maps.     

Mathematical and numerical modeling of the non-associated plasticity of soils—Part 1: The boundary value problem

The soil is characterized by the influence of the hydrostatic stress, which leads to a yield surface with a shape of a pyramid for Mohr–Coulomb criteria and a shape of a cone for Drucker–Prager one. These materials are also characterized by a non-associated plasticity where the plastic yielding rule does not follow the normality rule. The usual mechanical models use two independent functions to describe this particular collapse. Unfortunately, this manner broke the model formulation. The purpose of this work is to present a consistent formulation of the non-associated plasticity of soil. The frame of the mathematical analysis is the concept of the implicit standard material. The cornerstone of this new idea is the construction of a single function called the bipotential playing in the same time the roles of the yield surface and the plastic potential. The bipotential concept is then intended to involve the constitutive law, cover the normality rule even for the non-associated soil and the proof of the solution existence. The formulation was initially performed for the case of a regular point out of the cone apex and in present, it is extended to the irregular point located at the apex. The paper presents firstly the implicit standard material method. Then, the methodology to build a full model for the boundary value problem is detailed. Particular expressions and relations are sufficiently explained and discussed. Attention is made to the evolution problem and the variational principles related to the elastic–plastic behavior.     

基于非关联流动规律的Gotoh屈服准则的参数确定方法

屈服准则对板料成形过程的理论解析、工艺优化和有限元模拟有着重要的影响. 通过提高屈服准则的各向异性表征能力, 可以确保成形过程的可靠性及实际预测的准确性. 本文基于非关联流动法则, 给出了Gotoh屈服准则一套全新的参数求解方法. 在结合常用屈服准则并考虑流动规律的基础上, 分别以5754O铝合金、DP980先进高强钢和SAPH440结构钢作为研究对象, 进行了不同加载路径下各向异性变形行为的预测. 根据Gotoh屈服准则推导的屈服函数、塑性势函数以及基于关联流动的理论函数计算出屈服应力和各向异性指数$r$值随加载角度的分布趋势, 进而针对平面应力状态的屈服轨迹展开分析, 验证了不同屈服准则和流动规律对各向异性屈服行为的预测精度. 理论与实验数据对比结果表明: 不同屈服准则针对同种板料在流动规律一致的情形下其表征各向异性的能力有显著差异; 相同屈服准则基于不同流动规律其表征能力也具有明显差别. 基于非关联流动的屈服准则能极大地提高精度, 各向异性表征能力显著加强. 相关结果能够为各向异性屈服准则在塑性成形领域的实际应用方案提供重要参考.      

Three-dimensional limit analysis of rigid blocks assemblages. Part I: Torsion failure on frictional interfaces and limit analysis formulation

The governing equations of plastic torsion of, arbitrarily shaped, frictional interfaces are established and a general solution is proposed. The case of rectangular interfaces is used to study the interactions of torsion strength with bending moments and shear forces. In order to solve limit analysis problems, a piecewise linear approximation of the yield function for rectangular interfaces is proposed. A model for the limit analysis of three-dimensional block assemblages interacting through no-tension frictional interfaces is presented including the proposal for the torsion failure mode. This model takes into account non-associated flow rules and limited compressive stresses at the interfaces.     

Unified elastoplastic finite difference and its application

Two elastoplastic constitutive models based on the unified strength theory (UST) are established and implemented in an explicit finite difference code, fast Lagrangian analysis of continua (FLAC/FLAC3D), which includes an associated/non-associated flow rule, strain-hardening/softening, and solutions of singularities. Those two constitutive models are appropriate for metallic and strength-different (SD) materials, respectively. Two verification examples are used to compare the computation results and test data using the two-dimensional finite difference code FLAC and the finite element code ANSYS, and the two constitutive models proposed in this paper are verified. Two application examples, the large deformation of a prismatic bar and the strain-softening behavior of soft rock under a complex stress state, are analyzed using the three-dimensional code FLAC3D. The two new elastoplastic constitutive models proposed in this paper can be used in bearing capacity evaluation or stability analysis of structures built of metallic or SD materials. The effect of the intermediate principal stress on metallic or SD material structures under complex stress states, including large deformation, three-dimensional and non-association problems, can be analyzed easily using the two constitutive models proposed in this paper.     

Closed set of the uniqueness conditions and bifurcation criteria in generalized coupled thermoplasticity for small deformations

This paper reports the results of a study into global and local conditions of uniqueness and the criteria excluding the possibility of bifurcation of the equilibrium state for small strains. The conditions and criteria are derived on the basis of an analysis of the problem of uniqueness of a solution involving the basic incremental boundary problem of coupled generalized thermo-elasto-plasticity. This work forms a follow-up of previous research (?loderbach in Bifurcations criteria for equilibrium states in generalized thermoplasticity, IFTR Reports, 1980, Arch Mech 3(35):337–349, 351–367, 1983), but contains a new derivation of global and local criteria excluding a possibility of bifurcation of an equilibrium state regarding a comparison body dependent on the admissible fields of stress rate. The thermal elasto-plastic coupling effects, non-associated laws of plastic flow and influence of plastic strains on thermoplastic properties of a body were taken into account in this work. Thus, the mathematical problem considered here is not a self-conjugated problem.