两种晶体塑性损伤模型在有限变形条件下的比较

对建立在连续损伤力学和内变量理论基础上的两种晶体塑性损伤模型进行了比较,旨在比较两种模型在描述材料物理性能方面的适用性以及由加载而引起变形响应的不同之处.模拟结果显示,两种模型均能反映出塑性各向异性和损伤的演化;由加载而导致有限变形的响应不仅依赖于变形,而且也依赖于晶格取向;尽管两种模型在揭示单晶体的物理性能方面是不同的,但是在预测材料力学性能方面有着相同的预测趋势.     

延性单晶非均匀变形的三维有限元模拟

对延性单晶在拉伸载荷作用下的应变局域化和颈缩等非均匀变形过程进行了三维有限元数值模拟。将相关晶体塑性本构模型及一种新的数值积分方法补充到ABAQUS6.1商用有限元软件中。该方法的特点是，利用晶体塑性的动力学方程，获得一个关于晶体弹性变形梯度的演化方程，采用半隐式积分方案进行求解。本文推导出一种新的应力变本构矩阵。按此方式更新本构矩阵，计算速度和计算稳定性大大提高。加载方式，边界条件和变形程度等因素影响着滑移系的启动状况，这是平面模型所不能预测的。本文利用三维有限元方法模拟了不同取向下滑移系的启动状况，全面地考虑了FCC单晶材料12个可能滑移系在变形过程中的启动状况，合理地模拟了FCC面心立方单晶沿不同取向加载时晶轴旋转导致的应变局域化和颈缩等非均匀变形过程。     

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

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

铜单晶拉伸试样表面滑移带痕迹的晶体塑性分析

采用建立在晶体塑性理论基础上的晶体塑性有限变形计算方法,针对铜单晶试样单轴拉伸过程中晶体滑移在试样表面留下的滑移带痕迹进行了数值研究.作者利用三维有限元模拟不同取向铜单晶试样的拉伸变形,通过晶体塑性滑移面与试样表面交线的几何分析,得到了试样在不同晶向拉伸下不同滑移系启动造成的试样表面滑移痕迹,并对数值计算的试样表面滑移...     

材料三维微结构表征及其晶体塑性有限元模拟

材料的力学性能，尤其是在有限变形下所呈现的宏观各向异性，是材料结构设计和服役寿命考虑的关键因素。由于宏观模型不能较好地反映材料微观结构（晶粒的形貌和取向等）对宏观塑性各向异性的影响，因此，本文建立了能实际反映晶粒形貌的三维Voronoi模型，并基于晶体塑性理论对铝合金在有限变形下的响应进行计算。首先，建立反映材料微结构的代表性体积单元RVE模型进行计算，并与实验结果进行对比验证。然后，以单向拉伸为例，分析了有限变形过程中试件的晶粒形貌和取向分布等微观因素对宏观各向异性演化的影响，并从材料和结构两个层面讨论了微观结构对宏观力学性能的影响。结果表明，本文模型能够反映微观结构对宏观力学性能的影响，为实际生产制造领域构件的力学性能提供可靠的预测。     

晶粒数量对多晶集合体初始各向异性的影响

Taylor类多晶晶体粘塑性模型被用于研究晶粒数量对随机分布多晶体拉伸塑性各向异性的影响。分别沿包含不同晶粒数量的多晶集合体的各方向进行单向拉伸数值模拟实验，得到多晶集合体各方向在一定等效应变下的等效应力，并用云图和等高线表示在多晶体的参考球面上。定义了描述多晶集合体各向异性程度的参考指标。讨论了三种确定晶体随机取向的方法。计算结果表明：晶粒数量有限的多晶集合体的应力应变响应仍有一定的各向异性，且随着晶粒数量增多，多晶集合体的各向异性程度降低；就所包含晶粒数相同的多晶集合体来说，在确定晶粒随机取向时，选取不同的方法对它的各向异性程度也有一定的影响。     

镁合金晶体塑性本构模型与非均匀孪生变形分析

微结构演化对镁合金材料力学性能有着显著的影响,为了揭示镁合金宏观塑性各向异性特性与非均匀孪生变形的关系,开展了不同路径下的单轴加载试验以及采用含滑移、孪生机制的晶体塑性本构模型对试验条件下的镁合金变形行为进行数值模拟研究。文中本构模型描述了滑移与孪生变形机制以及晶格转动的机制,同时研究采用三维微结构代表性有限元模型,其包含晶粒尺寸、晶向和晶界倾角等微结构参数。研究结果表明,轧制镁合金具有强烈的宏观塑性各向异性行为,并对这种镁合金塑性各向异性行为的模拟结果以及多晶织构的模拟演化结果与试验测量进行对比,结果都基本吻合。对孪生非均匀变形模拟分析表明,镁合金宏观塑性各向异性行为与滑移、孪生变形机制的不同启动组合紧密相关,同时多晶体内应力的非均匀分布受到孪生变形的严重影响。而不同晶粒尺寸的晶粒所发生的孪生变形有比较大的差异,造成孪晶变体在晶粒内的分布极不均匀。本研究可为通过微结构的合理配置来设计和控制材料的力学性能提供理论依据.     

基于晶体塑性理论研究铝材料高压高应变率下的强度特性

为了了解金属材料在极端加载下复杂动态响应过程中的多种机制和效应,重点针对Al材料在高压、高应变率加载下的塑性变形机制,在经典晶体塑性模型的基础上,对其中的非线性弹性、位错动力学和硬化形式进行改进,建立适用于高压、高应变率加载下的热弹-黏塑性晶体塑性模型。该模型可以较好地描述单晶铝和多晶铝材料屈服强度随压力的变化过程,相比宏观模型,用该模型还获得了多晶Al材料在冲击加载下的织构演化规律,揭示了织构择优取向行为和压力的关系。     

考虑运动硬化的DD3镍基合金力学行为研究

详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,引入了运动硬化规律,考虑了温度和应变率对晶体塑性变形的影响,通过针对每个滑移系考虑屈服准则和流动规律建立了晶体塑性模型. 对积分过程进行了推导,通过编写ABAQUS材料用户子程序(UMAT), 实现本构模型的有限元积分算法. 在此基础上模拟了DD3镍基单晶合金在单轴拉伸和循环载荷下的响应,并与实验数据进行了对比. 利用该模型可以很好地模拟镍基单晶所具有的各向异性特性,体现了镍基单晶在循环载荷作用下的拉-压不对称性.      

复杂应力状态镍基单晶合金低周疲劳损伤模型

根据连续介质损伤力学理论,采用应变能释放率作为热力学广义力描述正交异性材料的疲劳损伤过程,引入取向函数考虑镍基单晶合金晶体取向对疲劳损伤的非线性影响,提出了一个各向异性疲劳损伤模型。应用多元线性回归分析方法,拟合疲劳试验数据可确定模型参数。从应变能释放率的应变空间表达式出发,导出了含有3个弹性常数的单晶合金应变三轴性因子,它既反映了材料性能的晶体取向相关性,又反映了正应力和剪应力的相互作用,并可退化为各向同性材料的应变三轴性因子。利用该模型对CMSX-2镍基单晶合金在应力控制对称循环拉-扭载荷作用下的低周疲劳寿命进行预测,预测值与试验值吻合的相当好,试验所得数据均落在2.2倍偏差的分布带内。     

Interpretation of the behavior of metals under large plastic shear deformations: comparison of macroscopic predictions to physically based predictions

A large plastic shear problem is analyzed by application of a macroscopic anisotropic plasticity model (Kuroda, M., 1997. Interpretation of the behavior of metals under large plastic shear deformations: a macroscopic approach. Int. J. Plasticity 13, 359–383), and the results are compared to predictions based on crystal plasticity with the Taylor assumption. It is found that these two different-scale models provide very similar predictions. The interpretations for such similarities are pursued in detail. The present macroscopic model reproduces quite well the change in orientation of anisotropy, which is directly predicted in the crystal plasticity analyses as a macroscopic manifestation of texture development. Consequently, the predictions for the rotation of the yield surface by the different-scale models become very similar. It is clearly shown that, in a macroscopic sense, the rotation of the anisotropic yield surface is a main cause of the axial effects in large plastic shear deformation.     

织构可调塑性本构模型标定及其应用

提出了利用率相关晶体塑性模型标定织相可调本构模型的求解步骤,得出了一组依赖于晶粒间相互作用假设而独立于具体板材织构的本构相关系数.以此为基础再结合板材织构系数所得出的本构模型系数可避免出现屈服面非外凸的情形.利用所提求解步骤对在不同热处理条件下产生不同织构的AL5052铝合金板的深拉成形过程进行了有限元模拟.结果再现了典型织构在板材成形过程中所出现的塑性各向异性,从而表明求解步骤的可行性.     

Plastic potentials for anisotropic porous solids

The aim of this paper is to incorporate plastic anisotropy into constitutive equations of porous ductile metals. It is shown that plastic anisotropy of the matrix surrounding the voids in a ductile material could have an influence on both effective stress–strain relation and damage evolution. Two theoretical frameworks are envisageable to study the influence of plastic flow anisotropy: continuum thermodynamics and micromechanics. By going through the Rousselier thermodynamical formulation, one can account for the overall plastic anisotropy, in a very simple manner. However, since this model is based on a weak coupling between plasticity and damage dissipative processes, it does not predict any influence of plastic anisotropy on cavity growth, unless a more suitable choice of the thermodynamical force associated with the damage parameter is made. Micromechanically-based models are then proposed. They consist of extending the famous Gurson model for spherical and cylindrical voids to the case of an orthotropic material. We derive an upper bound of the yield surface of a hollow sphere, or a hollow cylinder, made of a perfectly plastic matrix obeying the Hill criterion. The main findings are related to the so-called ‘scalar effect’ and ‘directional effect’. First, the effect of plastic flow anisotropy on the spherical term of the plastic potential is quantified. This allows a classification of sheet materials with regard to the anisotropy factor h; this is the scalar effect. A second feature of the model is the plasticity-induced damage anisotropy. This results in directionality of fracture properties (‘directional effect’). The latter is mainly due to the principal Hill coefficients whilst the scalar effect is enhanced by ‘shear’ Hill coefficients. Results are compared to some micromechanical calculations using the finite element method.     

塑性各向异性板料本构关系及冲压成形的数值模拟

推导了具有一般屈服函数形式的弹塑性速率型本构关系;给出了用于板料成形的Hill塑性各向异性屈服模型下本构关系的具体形式;用有限元动力显式计算程序MSC/DYTRAN模拟了金属板料的冲压成形;通过算例分析,考察了塑性各向异性对凸耳形成和大小以及对成形模拟结果准确性的影响;数值结果和实验结果表明：各向（厚向）异性本构模型比各向同性本构模型更真实地反映了板料的成形性。     

Unit cell debonding analyses for arbitrary orientation of plastic anisotropy

Debonding of rigid inclusions embedded in the elastic–plastic aluminum alloy Al 2090-T3 is analyzed numerically using a unit cell model taking full account of finite strains. The cell is subjected to overall biaxial plane strain tension and periodical boundary conditions are applied to represent arbitrary orientations of plastic anisotropy. Plastic anisotropy is considered using two phenomenological anisotropic yield criteria, namely Hill [Proceedings of the Royal Society of London A 193 (1948) 281] and Barlat et al. [International Journal of Plasticity 7 (1991) 693]. For this material plastic anisotropy delays debonding compared to plastic isotropy except for the case of Hill’s yield function when the tensile directions coincided with the principal axes of anisotropy. For some inclinations of the principal axes of anisotropy relative to the tensile directions, the stress strain responses are identical but the deformation modes are mirror images of each other.     

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.