考虑塑性应变记忆恢复的循环棘轮本构模型研究

在统一粘塑性循环本构理论框架下,以Ohno-Abdel-Karim非线性随动硬化模型为基础,建立了一个循环本构模型。模型通过引入塑性应变幅值记忆效应,并在塑性应变记忆项中加入恢复系数,提高了对循环硬化材料单轴棘轮行为的预言能力。将模型应用于316L不锈钢单轴棘轮行为的描述中,模拟不同平均应力、应力幅值下的棘轮应变,均与实验数据吻合较好,证明本文改进的本构模型能合理地描述循环硬化材料的单轴棘轮行为。     

SS304不锈钢高温非比例多轴循环变形行为的粘塑性本构描述

在350℃和700℃下SS304不锈钢的非比例多轴循环变形行为进行了系统的实验研究。在此基础上,在统一粘塑性本构理论的框架下改进和发展了一个新的循环本构模型。该模型给出了新的背应力演化方程,引入了非比例度参量,考虑了温度效应和最大塑性应变幅值记忆效应,能够对材料的高温非比例多轴应变循环变形行为和棘轮行为进行统一描述。模型的模拟结果与实验结果比较表明:该模型对SS304不锈钢高温非比例多轴循环变形行为的描述比较合理。     

考虑细观变形特征的镁合金宏观本构模型及数值模拟

为了能有效描述镁合金宏观各向异性塑性行为,考虑了滑移、孪生、去孪生三种细观变形模式的特点,给出了相应的硬化函数;根据VonMises屈服准则,发展了一种镁合金宏观本构模型及其迭代算法。模型将变形模式的开启与晶粒取向相关联,同时针对镁合金孪生变形时引起的晶粒重新定向问题,描述了一种晶向偏转的方法。在此基础上编写了ABAQUS/UMAT材料用户子程序;利用开发的本构模型,开展了单轴拉伸、单轴压缩、单轴循环拉压加载条件下镁合金塑性行为的数值模拟,并对随机织构下的镁合金板材轧制过程进行了有限元仿真实验。模拟结果表明:单轴拉伸、单轴压缩和循环加载情形下的镁合金宏观硬化行为与实验结果基本吻合;轧制后镁合金板材表现出了应力-应变不均匀特性,多晶织构演化结果与实验结果基本一致。说明文中所提出的宏观本构模型、晶向偏转模型能够有效描述镁合金的宏观塑性行为和织构演化。     

一种新的高应变率复合压剪实验技术

提出了一种新的基于Hopkinson杆实验技术的在102~103s-1高应变率下实现压剪复合加载的实验装置,并给出了相应的理论分析和数值模拟。为了获取材料在复杂应力下的本构关系,借助斜飞片冲击实验的思想,对Hopkinson杆进行改造,将入射杆的末端改进为截锥形,以便在试样中同时产生压缩和剪切应力。利用有限元分析软件LS-DYNA对试样中的应力波传播进行模拟计算,并利用改进装置进行了初步实验。计算和分析结果表明,利用所设计的装置可以实现对试样的动态压剪复合加载,获得材料在高应变率复杂应力加载下的本构响应,进而建立材料在复杂应力状态下本构行为的描述。     

退火和调质42CrMo合金钢单轴循环塑性行为的实验研究

在室温下对退火和调质42CrMo合金钢进行了单轴应变控制和应力控制的系统循环实验,并对它们的应变循环和应力循环特性进行比较.揭示和分析了应变幅值、平均应变及其历史对材料应变循环特性的影响以及应力幅值、平均应力及其历史对棘轮行为的影响.讨论了应变循环和应力循环间的交互作用以及不同热处理工艺下材料循环变形行为间的区别.研究发现材料的热处理工艺、平均应力和应力幅值及其历史对材料的棘轮行为都有很大的影响.得到了一些有助于进行合理本构描述的结果.     

304不锈钢室温单轴循环棘轮行为的粘塑性本构描述

在统一粘塑性循环本构模型的框架下对循环硬化的304不锈钢的单轴棘轮行为进行了本构描述。模型中通过随动硬化背应力演化和各向同性变形阻力演化对304不锈钢在非对称应力循环下的循环附加硬化和循环流动特性进行了分析，同时考虑了加载历史对循环棘轮行为的影响。将模型应用于304不锈钢室温单轴循环棘轮行为及其对加载历史依赖性的描述中，预言结果与实验结果吻合较好。     

金属材料在复杂应力状态下的塑性流动特性及本构模型

工程应用中,金属材料和结构往往处于复杂应力状态。材料的塑性行为会受到应力状态的影响,要精确描述材料在复杂应力状态下的塑性流动行为,必须在本构模型中考虑应力状态效应的影响。然而,由于在动态加载下材料的应变率效应和应力状态效应相互耦合、难以分离,给应力状态效应的研究和模型的建立造成很大困难。通过对Ti-6Al-4V钛合金材料开展不同加载条件下的力学性能测试,提出了一个包含应力三轴度和罗德角参数影响的新型本构模型,并通过VUMAT用户子程序嵌入ABAQUS/Explicit软件。分别采用新提出的塑性模型和Johnson-Cook模型对压剪复合试样的动态实验进行了数值模拟。结果表明,新模型不仅在对材料本构曲线的拟合方面具有较强的优势,而且由该模型所得到的透射脉冲和载荷-位移曲线均更加准确。因此,该模型能够更精确地描述和预测金属材料在复杂应力状态下的塑性流变行为。     

各向同性超弹性本构模型数值计算及验证

现有多种形式的橡胶本构模型试图预测橡胶力学性质,其中部分模型已写入有限元软件中用于仿真计算,还存在较多拟合性较好的模型无法在有限元材料库中直接获得。本文详述了由不变量和主伸长率描写的各向同性超弹性本构模型的数值实现方法,并结合最新的本构模型开发了UHYPER和UMAT子程序。将UHYPER用于有限元实现对多孔橡胶板的拉伸仿真,对比仿真和试验结果,验证子程序的正确性以及评估本构模型预测复杂应变场的准确性;将UMAT用于单轴、等双轴和剪切拉伸的有限元仿真,对比仿真和本构模型理论结果,验证子程序的可靠性。结果表明,有限元仿真结果与理论结果拟合较好,子程序能够契合本构模型的力学描述,所述方法可以用于超弹性材料的数值计算。     

304不锈钢室温和高温单轴循环塑性的实验研究

对304不锈钢进行了室温和高温单轴应变控制和应力控制下的系统循环试验。揭示和分析了循环应变幅值、平均应变及其历史和温度历史对材料应变循环特性的影响以及应力幅值、平均应力及其历史以及温度对循环棘轮行为的影响。也讨论了应变循环和应力循环间交互作用对材料循环塑性行为的影响。研究表明,无益单轴应变循环特性还是非对称单轴应力循环下的棘轮效应不仅取决于当前温度和加载状态,而且强烈依赖于其加载历史。研究得到了一些有助于304不锈钢室温和高温单轴循环行为本构描述的结果。     

饱和和非饱和介质的弹塑性损伤模型

针对饱和和非饱和工程材料变形和强度特点,给出了一个具体的弹塑性损伤本构方程,该本构方程可以描述材料性质的劣化过程（即损伤和软化）以及饱和程度对材料变形强度特性的影响。此外,还讨论了该弹塑性损伤本构方程的有限元实施方法,编制了专用的有限元程序。对常规三轴试伯进行了有限元分析,得出了一些具体的结果。     

Uniaxial and non-proportionally multiaxial ratcheting of SS304 stainless steel at room temperature: experiments and simulations

The uniaxial and non-proportionally multiaxial ratcheting behaviors of SS304 stainless steel at room temperature were initially researched by experiment and then were theoretically described by a cyclic constitutive model in the framework of unified visco-plasticity. The effects of cyclic stress amplitude, mean stress, and their histories on the ratcheting were experimentally investigated under uniaxial and different multiaxial loading paths. The shapes of non-proportional loading paths were linear, circular, elliptical and rhombic, respectively. In the constitutive model, the rate-dependent behavior of the material was reflected by a viscous term; the cyclic flow and cyclic hardening behaviors of the material under asymmetrical stress-controlled cycling were reflected by the evolution rules of kinematic hardening back stress and isotropic deforming resistance, respectively. The effect of loading history on the ratcheting was also considered by introducing two fading memorization functions for maximum inelastic strain amplitude and isotropic deformation resistance, respectively, into the constitutive model. The effect of multiaxial loading path on the ratcheting was reflected by a non-proportional factor defined in this work. The predicting ability of the developed model was proved to be good by comparing the simulations with corresponding experiments.     

非比例循环塑性内时本构关系研究

在Valanis的内时本构理论的框架中,引入内结构张量以反映由于非比例加载而引起金属材料的附加等向强化及异向强化效应,同时提出材料强化程度的度量采用沿路径法线方向的塑性应变分量来描述.这些考虑的有效性已经通过用所建模型对304不锈钢材料在一些典型非比例循环加载路径下的响应进行的理论预测得到了验证;将该模型应用于U71Mn材料室温单轴棘轮行为描述中,结果显示内结构张量的引入不仅能较好地反映应变控制下的非比例附加效应,而且也能较好地反映应力控制下塑性应变的累积及变化率.     

Reconcile the perfectly elastoplastic model to simulate the cyclic behavior and ratcheting

Perfectly elastoplastic constitutive model is modified through a smoothing factor introduced by Liu [Liu, C.-S., 2003. Smoothing elastoplastic stress–strain curves obtained by a critical modification of conventional models. Int. J. Solids Struct. 40, 2121–2145]. The new model allows plasticity to happen in a non-zero-measure yield volume in stress space, rather than that of conventional zero-measure yield surface, and within the yield volume the plastic modulus is varying continuously. It endows a specific strain-hardening rule of flow stress and is able to describe the phenomena of strain hardening, cyclic hardening, the Bauschinger effect, mean-stress relaxation, strain ratcheting, out-of-phase hardening, as well as erasure-of-memory. In order to suppress the over prediction of ratcheting we consider a scalar function of smoothing factor, which can simulate the saturation behavior of uniaxial/multiaxial strain ratcheting. These effects are demonstrated through numerical examples. The existence of stress equilibrium point and limiting surface is a natural result without requiring an extra design. Moreover, the non-linear constitutive equations can be converted into a linear system for augmented stress in the Minkowski space, of which the symmetry group is a proper orthochronous Lorentz group SO_o(5, 1). The augmented stress is a time-like vector, moving on hyperboloids inside the cone. When taking the Prager kinematic hardening rule into account we can simulate some cyclic behaviors of SAE 4340 and grade 60 steels within a certain accuracy through the use of only three material constants and a fixed smoothing factor. To simulate the ratcheting behaviors of SS304 stainless steel we allow the smoothing factor to be an exponential decaying function of λ.     

考虑尺寸效应的多晶金属循环塑性分析

本文采用多晶塑性分析方法,设材料点包含一定数量的各向异性单晶晶粒并考虑晶粒尺寸的影响,计算材料点的应力和应变时利用了Taylor假设。模型引入考虑尺寸效应的晶体滑移硬化函数,同时针对晶体滑移引入背应力及其方向性硬化的描述,以反映不同晶粒尺寸材料在循环加载条件下的力学行为。利用该模型,本文第一作者采用显式格式编制了与ABAQUS商用有限元软件接口的用户材料子程序（VUMAT）,实例计算证实该模型可以反映和描述多晶金属材料在材料反复加载条件下的循环塑性行为与尺寸效应。     

SUPER-ELASTIC CONSTITUTIVE MODEL CONSIDERING PLASTICITY AND ITS FINITE ELEMENT IMPLEMENTATION

Based on the experimental results of super-elastic NiTi alloy, a three-dimensional super-elastic constitutive model including both of stress-induced martensite transformation and plasticity is constructed in a framework of general inelasticity. In the proposed model, transformation hardening, reverse transformation of stress-induced martensite, elastic mismatch between the austenite and martensite phases, and temperature-dependence of transformation stress and elastic modulus of each phase are considered. The plastic yielding of martensite occurred under high stress is addressed by a bilinear isotropic hardening rule. Drucker-Prager-typed transformation surfaces are employed to describe the asymmetric behavior of NiTi alloy in tension and compression. The prediction capability of the proposed model is verified by comparing the simulated results with the correspondent experimental ones. Based on backward Euler＇s integration, a new expression of consistent tangent modulus is derived. The proposed model is then implemented into a finite element package ABAQUS by user-subroutine UMAT. Finally, the validity of such implementation was verified by some numerical samples.     

Cyclic plasticity modeling with the distribution of non-linear relaxations approach

Motivated by the distribution of non-linear relaxation (DNLR) approach, a phenomenological model is proposed in order to describe the cyclic plasticity behavior of metals under proportional and non-proportional loading paths with strain-controlled conditions. Such a model is based on the generalization of the Gibbs's relationship outside the equilibrium of uniform system and the use of the fluctuation theory to analyze the material dissipation due to its internal reorganization. The non-linear cyclic stress–strain behavior of metals notably under complex loading is of particular interest in this study. Since the hardening effects are described appropriately and implicitly by the model, thus, a host of inelastic behavior of metals under uniaxial and multiaxial cyclic loading paths are successfully predicted such as, Bauschinger, strain memory effects as well as additional hardening. After calibrating the model parameters for two metallic materials, the model has demonstrated obviously its ability to describe the cyclic elastic-inelastic behavior of the nickel base alloy Waspaloy and the stainless steel 316L. The model is then implemented in a commercial finite element code simulating the cyclic stress–strain response of a thin-walled tube specimen. The numerical responses are in good agreement with experimental results.     

Computational modeling of inelastic large ratcheting strains

A framework for phenomenological hyperelasto-plasticity with initial anisotropy, kinematic hardening as well as anisotropic damage is presented in [Menzel et al., Int. J. Plasticity (2004), in press]. In this contribution, we exploit and extend this framework to include several back-stresses in order to capture the ratcheting response of polycrystalline metals subjected to cyclic stress with non-zero mid-value. The evolution equations for kinematic hardening resemble a linear combination of the multiple-Armstrong–Frederick and the Burlet–Cailletaud models, which are extended to the large strain setting. The capability of the model to capture various phenomenological characteristics, in particular multi-axial ratcheting, is illustrated by numerical examples. Comparisons with uni-axial and bi-axial experimental ratcheting results for carbon steel are given. Finally, the finite element analysis of a simplified railway turnout component subjected to cyclic loading is presented.     

On the implementation of a multi-surface kinematic hardening plasticity and its applications

This paper is concerned with an application of the multi-surface plasticity in solid mechanics and geotechnical problems. The model is of a von-Mises type with associated flow rule, originally proposed by Montans. The Mroz translation rule is implemented to the movements of the yield surfaces and the fully implicit scheme with radial mapping method is applied in numerical computations. Algorithmic consistent tangent modulus with numerical integration algorithm of constitutive equations is extracted. The model is developed in the class of kinematic hardening models, so the ‘Masing’ rule is preserved. The model is able to consider the plastic strain accumulation in constant axial stress state, such as ratcheting. The implementation is validated by means of a simple deformation path of combined extension and compression test, a pure shear test with pseudo-random loading, a test which demonstrates the capabilities of the model in simulation of cyclic loading and ratcheting, a cyclic shear test in saturated undrained sand and finally, the analysis of a plate with holes, which presents the shear band using the multi-surface plasticity model.