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.     

Modelling inelastic behaviour of orthotropic metals in a unique alignment of deviatoric plane within the stress space

A finite strain constitutive model to predict the deformation behaviour of orthotropic metals is developed in this paper. The important features of this constitutive model are the multiplicative decomposition of the deformation gradient and a new Mandel stress tensor combined with the new stress tensor decomposition generalized into deviatoric and spherical parts. The elastic free energy function and the yield function are defined within an invariant theory by means of the structural tensors. The Hill’s yield criterion is adopted to characterize plastic orthotropy, and the thermally micromechanical-based model, Mechanical Threshold Model (MTS) is used as a referential curve to control the yield surface expansion using an isotropic plastic hardening assumption. The model complexity is further extended by coupling the formulation with the shock equation of state (EOS). The proposed formulation is integrated in the isoclinic configuration and allows for a unique treatment for elastic and plastic anisotropy. The effects of elastic anisotropy are taken into account through the stress tensor decomposition and plastic anisotropy through yield surface defined in the generalized deviatoric plane perpendicular to the generalized pressure. The proposed formulation of this work is implemented into the Lawrence Livermore National Laboratory-DYNA3D code by the modification of several subroutines in the code. The capability of the new constitutive model to capture strain rate and temperature sensitivity is then validated. The final part of this process is a comparison of the results generated by the proposed constitutive model against the available experimental data from both the Plate Impact test and Taylor Cylinder Impact test. A good agreement between experimental and simulation is obtained in each test.     

混凝土材料冲击特性的研究

基于混凝土材料强冲击加载下的试验研究，提出了两种损伤型动态本构模型： 损伤型黏弹性本构模型和损伤与塑性耦合的本构模型. 通过模型计算结果与冲击试验结果的 比较可发现，随着冲击速度的提高，混凝土材料内部产生了显著的塑性变形，由此 损伤型黏弹性本构模型的应用就存在一些不足. 而损伤与塑性耦合的本构模型由于考虑了裂 纹扩展引起的材料强度和刚度的弱化，以及微空洞缺陷塌陷引起的塑性变形，因而能更好地 用于模拟强冲击载荷作用下混凝土材料的冲击响应特性.     

Effects of stress invariants and reverse loading on ductile fracture initiation

Recent research studies on ductile fracture of metals have shown that the ductile fracture initiation is significantly affected by the stress state. In this study, the effects of the stress invariants as well as the effect of the reverse loading on ductile fracture are considered. To estimate the reduction of load carrying capacity and ductile fracture initiation, a scalar damage expression is proposed. This scalar damage is a function of the accumulated plastic strain, the first stress invariant and the Lode angle. To incorporate the effect of the reverse loading, the accumulated plastic strain is divided into the tension and compression components and each component has a different weight coefficient. For evaluating the plastic deformation until fracture initiation, the proposed damage function is coupled with the cyclic plasticity model which is affected by all of the stress invariants and pervious plastic deformation history.For verification and evaluation of this damage-plasticity constitutive equation a series of experimental tests are conducted on high-strength steel, DIN 1.6959. In addition finite element simulations are carried out including the integration of the constitutive equations using the modified return mapping algorithm. The modeling results show good agreement with experimental results.     

复杂加载下混凝土的弹塑性本构模型

混凝土材料在不同应力路径下或复杂加载条件下会表现出差异性显著的应力应变关系,在小幅循环加载条件下,其应力应变关系会表现出类似于弹性变形的滞回曲线.在不同应力水平下,混凝土的应力应变关系以及破坏特性都具有静水压力相关特点,即随着静水压力增大,各向异性强度特性弱化.此外,混凝土受压及受拉破坏机理不同,因而对应于混凝土硬化损伤亦有不同,即可分为受压硬化损伤,受拉硬化损伤及两者的混合硬化损伤类型.基于Hsieh模型,对该模型进行了三点改进.（1）针对小幅循环加载下混凝土无塑性变形的试验规律,而模型中在应力水平较低的循环加载条件下始终存在塑性变形的预测问题,采用在边界面模型框架下,设置了应力空间的弹性域,初始屈服面与后续临界状态屈服面几何相似的假定.（2）基于广义非线性强度准则将原模型采用变换应力方法将其推广为三维弹塑性本构模型,采用变换后模型可合理的考虑不同应力路径对于子午面以及偏平面上静水压力效应形成的影响,并避免了边界面应力点奇异问题.（3）分别对拉压两种加载损伤模式建议了相应的硬化参数表达式,可分别用于描述上述加载中产生的应变软化及强度退化行为.基于多种加载路径模拟表明:所建立的三维弹塑性本构模型可合理地用于描述混凝土的一般应力应变关系特性.      

Simulation of damage evolution in ductile metals undergoing dynamic loading conditions

A set of constitutive equations for large rate-dependent elastic-plastic-damage materials at elevated temperatures is presented to be able to analyze adiabatic high strain rate deformation processes for a wide range of stress triaxialities. The model is based on the concepts of continuum damage mechanics. Since the material macroscopic thermo-mechanical response under large strain and high strain rate deformation loading is governed by different physical mechanisms, a multi-dissipative approach is proposed. It incorporates thermo-mechanical coupling effects as well as internal dissipative mechanisms through rate-dependent constitutive relations with a set of internal variables. In addition, the effect of stress triaxiality on the onset and evolution of plastic flow, damage and failure is discussed.Furthermore, the algorithm for numerical integration of the coupled constitutive rate equations is presented. It relies on operator split methodology resulting in an inelastic predictor-elastic corrector technique. The explicit finite element program LS-DYNA augmented by an user-defined material subroutine is used to approximate boundary-value problems under dynamic loading conditions. Numerical simulations of dynamic experiments with different specimens are performed and good correlation of numerical results and published experimental data is achieved. Based on numerical studies modified specimens geometries are proposed to be able to detect complex damage and failure mechanisms in Hopkinson-Bar experiments.     

Thermomechanical constitutive equations for the dynamic response of ceramics

The behavior and failure of brittle materials is significantly influenced by the existence of inhomogeneities such as pores and cracks. The proposed constitutive equations model the coupled micro-mechanical response of these inhomogeneities through evolution equations for scalar measures of porosity, and a “density” function of randomly oriented penny-shaped cracks. A specific form for the Helmholtz free energy is proposed which incorporates the known Mie–Grüneisen constitutive equation for the nonporous solid. The resulting thermomechanical constitutive equations are valid for large deformations and the elastic response is hyperelastic in the sense that the stress is related to a derivative of the Helmholtz free energy. These equations allow for the simulation of the following physical phenomena exhibited by brittle materials: (1) high compressive strength compared with much lower tensile strength; (2) inelastic deformation due to growth and nucleation of cracks and pores instead of due to dislocation dynamics associated with metal plasticity; and (3) loss of integrity (degradation of elastic moduli) due to damage accumulation. The main features of the model are demonstrated by examples of cyclic loading in homogeneous deformation and by a simulation of a dynamic plate-impact experiment on AD85 ceramic. The theoretical predictions of the model are in excellent agreement with the dynamic experimental data.     

An anisotropic elastoplastic constitutive formulation generalised for orthotropic materials

This paper presents a finite strain constitutive model to predict a complex elastoplastic deformation behaviour that involves very high pressures and shockwaves in orthotropic materials using an anisotropic Hill’s yield criterion by means of the evolving structural tensors. The yield surface of this hyperelastic–plastic constitutive model is aligned uniquely within the principal stress space due to the combination of Mandel stress tensor and a new generalised orthotropic pressure. The formulation is developed in the isoclinic configuration and allows for a unique treatment for elastic and plastic orthotropy. An isotropic hardening is adopted to define the evolution of plastic orthotropy. The important feature of the proposed hyperelastic–plastic constitutive model is the introduction of anisotropic effect in the Mie–Gruneisen equation of state (EOS). The formulation is further combined with Grady spall failure model to predict spall failure in the materials. The proposed constitutive model is implemented as a new material model in the Lawrence Livermore National Laboratory (LLNL)-DYNA3D code of UTHM’s version, named Material Type 92 (Mat92). The combination of the proposed stress tensor decomposition and the Mie–Gruneisen EOS requires some modifications in the code to reflect the formulation of the generalised orthotropic pressure. The validation approach is also presented in this paper for guidance purpose. The \({\varvec{\psi }}\) tensor used to define the alignment of the adopted yield surface is first validated. This is continued with an internal validation related to elastic isotropic, elastic orthotropic and elastic–plastic orthotropic of the proposed formulation before a comparison against range of plate impact test data at 234, 450 and \({\mathrm {895\,ms}}^{\mathrm {-1}}\) impact velocities is performed. A good agreement is obtained in each test.     

A physically-based constitutive model for anisotropic damage in rubber-toughened glassy polymers during finite deformation

The present work focuses on the development of a physically-based model for large deformation stress-strain response and anisotropic damage in rubber-toughened glassy polymers. The main features leading to a microstructural evolution (regarding cavitation, void aspect ratio, matrix plastic anisotropy and rubbery phase deformation) in rubber-toughened glassy polymers are introduced in the proposed constitutive model. The constitutive response of the glassy polymer matrix is modelled using the hyperelastic-viscoplastic model of [Boyce et al., 1988] and [Boyce et al., 2000]. The deformation mechanisms of the matrix material are accounted for by two resistances: an elastic-viscoplastic isotropic intermolecular resistance acting in parallel with a visco-hyperelastic anisotropic network resistance, each resistance being modified to account for damage effects by void growth with a variation of the void aspect ratio. The effective contribution of the hyperelastic particles to the overall composite behaviour is taken into account by treating the overall system in a composite scheme framework. The capabilities of the proposed constitutive model are checked by comparing experimental data with numerical simulations. The deformation behaviour of rubber-toughened poly(methyl methacrylate) was investigated experimentally in tension at a temperature of 80 °C and for different constant true strain rates monitored by a video-controlled technique. The reinforcing phase is of the soft core-hard shell type and its diameter is of the order of one hundred nanometers. The particle volume fraction was adjusted from 15% to 45% by increments of 5%. The stress-strain response and the inelastic volumetric strain are found to depend markedly on particle volume fraction. For a wide range of rubber volume fractions, the model simulations are in good agreement with the experimental results. Finally, a parametric analysis demonstrates the importance of accounting for void shape, matrix plastic anisotropy and rubber content.     

Finite deformation of a polymer: experiments and modeling

The response of a polymer (polytetrafluoroethylene) to quasi-static and dynamic loading is determined and modeled. The polytetrafluoroethylene is extremely ductile and highly nonlinear in elastic as well as plastic behaviors including elastic unloading. Constitutive model developed earlier by Khan, Huang and Liang (KHL) is extended to include the responses of polymeric materials. The strain rate hardening, creep, and relaxation behaviors of polytetrafluoroethylene were determined through extensive experimental study. Based on the observation that both viscoelastic and viscoplastic deformation of polytetrafluoroethylene are time dependent and nonlinear, a phenomenalogical viscoelasto–plastic constitutive model is presented by a series connection of a viscoelastic deformation module (represented by three elements standard solid spring dashpot model), and a viscoplastic deformation module represented by KHL model. The KHL module is affected only when the stress exceeds the initial yield stress. The comparison between the predictions from the extended model and experimental data for uniaxial static and dynamic compression, creep and relaxation demonstrate that the proposed constitutive model is able to represent the observed time dependent mechanical behavior of polytetrafluoroethylene polytetrafluoroethylene qualitatively and quantitatively.     

A quasi-flow constitutive model with strain-rate dependence

In this paper, the proposed is a quasi-flow constitutive model with strain-rate sensitivity for elastic plastic large deformation. The model is based on the Quasi-flow Corner theory, and is suitable for the sheet metal forming process simulation with a variable punch machine velocity. Uniaxial tensile tests and deep-drawing tests of a circular blank with square punch are carried out and numerically simulated. The consistency between the experimental and the numerically simulated results shows the validity of the present new constitutive model. The project supported by the Scientific Foundation of National Outstanding Youth of China (10125208), the National Natural Science Foundation of China (19832020), and the National Education Committee of China     

Finite viscoplasticity of amorphous glassy polymers in the logarithmic strain space

The paper outlines a new constitutive model and experimental results of rate-dependent finite elastic–plastic behavior of amorphous glassy polymers. In contrast to existing kinematical approaches to finite viscoplasticity of glassy polymers, the formulation proposed is constructed in the logarithmic strain space and related to a six-dimensional plastic metric. Therefore, it a priori avoids difficulties concerning with the uniqueness of a plastic rotation. The constitutive framework consists of three major steps: (i) A geometric pre-processing defines a total and a plastic logarithmic strain measures determined from the current and plastic metrics, respectively. (ii) The constitutive model describes the stresses and the consistent moduli work-conjugate to the logarithmic strain measures in an analogous structure to the geometrically linear theory. (iii) A geometric post-processing maps the stresses and the algorithmic tangent moduli computed in the logarithmic strain space to their nominal, material or spatial counterparts in the finite deformation space. The analogy between the formulation of finite plasticity in the logarithmic strain space and the geometrically linear theory of plasticity makes this framework very attractive, in particular regarding the algorithmic implementation. The flow rule for viscoplastic strains in the logarithmic strain space is adopted from the celebrated double-kink theory. The post-yield kinematic hardening is modeled by different network models. Here, we compare the response of the eight chain model with the newly proposed non-affine micro-sphere model. Apart from the constitutive model, experimental results obtained from both the homogeneous compression and inhomogeneous tension tests on polycarbonate are presented. Besides the load–displacement data acquired from inhomogeneous experiments, quantitative three-dimensional optical measurements of the surface strain fields are carried out. With regard to these experimental data, the excellent predictive quality of the theory proposed is demonstrated by means of representative numerical simulations.     

An indentation technique for estimating the energy density as fracture toughness with Berkovich indenter for ductile bulk materials

A technique is proposed to estimate the energy density as fracture toughness for ductile bulk materials with an indentation system equipped with a Berkovich indenter based on the theory of plastic deformation energy transforming into the indentation energy of fracture. With progressive increase of penetration loads, the material damage is exhibited on the effective elastic modulus. A quadratic polynomial relationship between the plastic penetration depth and penetration load, and an approximate linear relationship between logarithmic plastic penetration depth and logarithmic effective elastic modulus are exhibited by indentation investigation with Berkovich indenter. The parameter of damage variable is proposed to determine the critical effective elastic modulus at the fracture point. And the strain energy density factor is calculated according to the equations of penetration load, plastic penetration depth and effective elastic modulus. The fracture toughness of aluminum alloy and stainless steel are evaluated by both indentation tests and K_IC fracture toughness tests. The predicted S_cr values of indentation tests are in good agreement with experimental results of CT tests.     

气泡在液体中运动过程的数值模拟

本文用数值方法预测气泡在液体中的百定常运动。运用位标函数进行界面的隐含跟踪并且与有限体积法相结合构成一种可行的计算方法。本文方法允许在界面处存在很大的物性差，而且较容易将表面张力引入控制方程。我们对气液两相流中单个气泡的运动进行了计算，得到了与实验结果符合很好的数值结果。     

含热冲击预损伤的陶瓷基复合材料损伤本构模型

陶瓷基复合材料结构在服役过程中不可避免地经受热冲击(较高的热应力梯度)而产生热机械损伤, 因此, 建立含循环热冲击预损伤材料的损伤本构模型, 以描述材料在热机械载荷作用下的力学行为, 对材料结构损伤容限设计与结构完整性评估非常重要. 本文首先对经历了循环热冲击的材料进行单调拉伸损伤实验, 发现对于含循环热冲击预损伤的材料, 其弹性模量的下降与所施加的应变直接相关. 然后在连续介质损伤力学的框架下, 基于平面应力假设, 建立了含循环热冲击预损伤材料的损伤演化模型, 该模型所涉及的参数可通过一个偏轴(45$^\circ$)以及两个正轴(平行于两个主方向)的单调拉伸试验获得. 最后, 采用经典塑性理论对由基体损伤引起的非弹性应变进行了描述. 本文所提出的应变损伤宏观模型可以描述陶瓷基复合材料在热机械载荷作用下的损伤演化, 同时弥补了含预损伤的陶瓷基复合材料在机械载荷下损伤本构模型在理论及实验研究方面的不足.      

Constitutive relations for isotropic or kinematic hardening at finite elastic–plastic deformations

The rate-type constitutive relations of rate-independent metals with isotropic or kinematic hardening at finite elastic–plastic deformations were presented through a phenomenological approach. This approach includes the decomposition of finite deformation into elastic and plastic parts, which is different from both the elastic–plastic additive decomposition of deformation rate and Lee’s elastic–plastic multiplicative decomposition of deformation gradient. The objectivity of the constitutive relations was dealt with in integrating the constitutive equations. A new objective derivative of back stress was proposed for kinematic hardening. In addition, the loading criteria were discussed. Finally, the stress for simple shear elastic–plastic deformation was worked out.