裂尖大应变细观断裂研究

本文用反映空穴形核成长的Gurson本构方程来描述裂尖区域材料在大应变情形下的力学特性，并进一步考虑了空穴演变对材料杨氏模量的影响。文中用上述本构方程分别结合弹塑性大应变有限元方法对平面应变I型裂纹问题作了计算，分析了裂尖应力分布、裂尖形状变化和裂尖空穴演变过程，并与用Prandtl-Reuss本构方程教育处的结果作了比较。     

金属韧性断裂的细观研究

对三点弯曲试件裂尖及断裂过程分别用 Prandtl-Reuss本构、Gurson本构方程作了大变形弹塑性有限元分析 ,并比较了二者的结果。对裂尖应力分布 ,应变分布和裂尖处空穴演化作了初步研究。计算模拟结果表明裂尖空穴化区域是一个很小量级的区域 ,采用 Gurson本构方程来研究裂尖“过程区”比采用 Prandtl-Reuss本构方程分析“过程区”来的合理 ,更接近材料实际断裂过程。     

脆性固体损伤和局部软化的有限元分析

本文研究混凝土、岩石一类工程中常用的应变软化材料的有限元分析方法。在作者以往有关粘塑性损伤本构模型的工作基础上，给出了一组便于有限元计算的本构方程表达式。包括损伤弹性矩阵和局部损伤软化矩阵，分别运用于计算硬化和软化阶段的有限元刚度矩阵；对所提出的本构方程的实验验证计算和一些算例的有限元数值分析，表明文中给出的本构方程是可行的，相应的有限元算法能较好地对损伤固体的局部软化效应进行数值分析，并可成功地追踪应力应变响应的软化曲线     

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

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

可压超弹性-塑性材料中的空穴生成

本文研究了材料的弹塑性性质对球体中空穴生成问题的影响，材料的弹性用一种可压超弹性材料的本构关系来描述，材料的塑性用满足材料的不可压条件和Tresca屈服条件的理想塑性材料的本构关系来描述。这类超弹性．塑性材料中可以发生空穴的生成现象，得到了在表面拉伸作用下球体中空穴生成时空穴半径与临界拉伸之间的关系式和临界拉伸。球体的变形可分为弹-塑性变形阶段和完全塑性变形阶段，球体中心首先形成塑性变形区域，并有空穴的突然生成；塑性变形区域能够快速增长，并且使球体很快进入完全塑性变形阶段；空穴在弹-塑性变形阶段迅速增长，但进入完全塑性变形阶段后增长较慢。同时给出了不同变形阶段球体中的应力分布。数值结果表明材料的塑性性质对材料中的空穴生成有明显的影响。     

粘塑性损伤模型模拟准超塑性单轴拉伸行为

发展了Chaboche粘塑性本构模型的大变形隐式算法，用损伤(DM)和无损伤(NDM)模型模拟准超塑性单轴拉伸。发现变形过程可分为三个阶段：均匀变形、颈缩发展、断裂破坏阶段。DM可准确模拟前两个阶段变形，NDM只能较好地模拟均匀变形阶段，表明DM可以较精确地描述稳定发展的动态过程。由于有限元方法只能描述连续介质，因此对于断裂破坏阶段，NDM模拟载荷大于试验结果，DM的载荷小于试验结果，这是由高应变速率敏感性造成。DM能够描述试验中出现地多处颈缩现象，局部应变速率分布随时间演化反映了颈缩发展程度。严重颈缩部位的距离代表着超塑性变形能力，距离越大，抗颈缩能力越好。     

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

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

基于试验与有限元耦合技术的延性材料全程单轴本构关系获取方法

摘 要: 材料拉伸直至断裂的全程单轴本构关系对材料大变形和断裂机理研究具有重要意义。传统拉伸试验获取的材料真应力-真应变曲线在试样颈缩后不可测。借助可以精确测量三维变形的DIC(Digital image correlate) 技术和有限元分析技术(Finite element analysis),本文提出了基于漏斗试样拉伸试验获取材料全程单轴本构关系的新方法,即TF(Test and FEA)方法。该方法将TF方法获取的材料全程单轴应力应变关系曲线作为有限元软件中的材料本构关系对漏斗试样拉伸变形过程进行模拟,其模拟载荷-位移曲线、漏斗根部直径-位移曲线和漏斗变形轮廓线等均与试验结果吻合良好,试样表面模拟应变也与DIC测试结果吻合, 根据不同半径漏斗试样模拟获得的全程真应力-真应变曲线保持良好一致性。最后,还对试样颈缩断面的力学行为进行了讨论,并给出了304不锈钢、汽轮机叶片材料2Cr12Ni4Mo3VNBN和 1Gr12Ni3Mo2VN、汽轮机转子材料30Cr2Ni4MoV的全程单轴本构关系模型参数、破断应力和破断应变。     

延性金属膨胀环的运动及断裂的三维数值模拟

利用LS-DYNA三维动力有限元软件对延性金属环的膨胀运动与断裂进行数值模拟。在膨胀环圆周加入泊松随机分布断裂成核点,利用J-C本构模型,研究诸如颈缩形成时间、颈缩区与均匀变形区的温度、应力、应变的对比等颈缩形成机理,以及讨论了环圆周上断裂成核点的泊松随机分布碎裂特性的影响。数值计算结果与实验结果、理论分析结果吻合较好,表明施加泊松随机断裂成核点的数值模拟方法是合理的。     

粗晶LY-12合金准超塑性单轴拉伸模拟

在试验的基础上，利用Chaboche本构模型对LY-12准超塑性单轴拉伸变形进行数值模拟，对其能够实现大延伸率的力学现象进行分析，由于材料具有高应变速率敏感指数，颈缩不发生在标距中间，而是有多处出现颈缩，不同颈缩相互的松弛作用对变形起重要作用，模拟和试验都表明，温度一定，不同应变速率下，颈缩发生的位置和颈缩程度不同，从而最大破坏延伸率也不同，利用该模型可以预测粗晶LY—12材料的超塑性变形能力和可能破坏的部位。     

伴随微孔洞生长的裂纹弹塑性定常扩展

裂纹在韧性材料中扩展时,将们随着微孔洞的萌生和生长,孔洞的萌生和深化将直接影响着材料的总体断裂韧性和强度,以往的研究主要集中在将裂纹的扩展刻划为微孔洞的萌生、生长和汇合这样一个过程。从传统的断裂过程区模型出发研究微孔洞的萌生和生长对材料总体断裂韧性的影响,通过采用Ｇｕｒｓｏｎ模型,建立塑性增量本构关系,然后针对定常扩展情况直接进行分析,孔洞对材料断裂韧性的影响由本构关系刻划,而在孔洞汇合模型中,上     

Analyses on micro-fracture of crack tip under large deformation conditions

In the present paper, Gurson's constitutive equation, which takes into account the development of voids, is used to study the behaviour of the material in the region near crack tip. Furthermore, the effect of void development on Young's modulus, which was not considered by Gurson, is taken into consideration. The analyses on void development, on stress distribution near crack tip, and on the variance of COD for the plane strain mode I problem are carried out with the large elastic-plastic deformation finite element method. The results are compared with those estimated from the Prandtl-Reuss constitutive equation.     

Continuum modeling of a porous solid with pressure-sensitive dilatant matrix

The pressure-sensitive plastic response of a material has been studied in terms of the intrinsic sensitivity of its yield stress to pressure and the presence and growth of cavities. This work focuses on the interplay between these two distinctly different mechanisms and the attendant material behavior. To this end, a constitutive model is proposed taking both mechanisms into account. Using Gurson's homogenization, an upper bound model is developed for a voided solid with a plastically dilatant matrix material. This model is built around a three-parameter axisymmetric velocity field for a unit sphere containing a spherical void. The void is also subjected to internal pressure; this can be relevant for polymeric adhesives permeated by moisture that vaporizes at elevated temperatures. The plastic response of the matrix material is described by Drucker–Prager's yield criterion and an associated flow rule. The resulting yield surface and porosity evolution law of the homogenized constitutive model are presented in parametric form. Using the solutions to special cases as building blocks, approximate models with explicit forms are proposed. The parametric form and an approximate explicit form are compared against full-field solutions obtained from finite element analysis. They are also studied for loading under generalized tension conditions. These computational simulations shed light on the interplay between the two mechanisms and its enhanced effect on yield strength and plastic flow. Among other things, the tensile yield strength of the porous solid is greatly reduced by the internal void pressure, particularly when a liquid/vapor phase is the source of the internal pressure.     

On the deformation and fracture of porous materials

The constitutive behavior of porous materials (including the yield loci, the void growth rate, the macro stress-strain relation and the strain to localization instability) is examined based on the lower bound approach proposed by the present authors. These results are then compared with the experimental and the finite element results as well as those predicted by Gurson's equations. Emphasis is placed on approaching the real behavior from the upper and the lower bound analysis. Calculation is also made on the influence of void nucleation on the critical strain to instability and a modified strain-controlled nucleation criterion is proposed. Finally the instability and fracture of AISI4340 steel in plane strain tension is examined and comparison is made between theoretical and experimental results.     

The void-size effect on plastic flow localization in the Gurson model

Recent studies have shown that the size of microvoids has a significant effect on the void growth rate. The purpose of this paper is to explore whether the void size effect can influence the plastic flow localization in ductile materials. We have used the extended Gurson‘s dilatational plasticity theory, which accounts for the void size effect, to study the plastic flow localization in porous solids with long cylindrical voids. The localization model of Rice is adopted, in which the material inside the band may display a different response from that outside the band at the incipient plastic flow localization. The present study shows that it has little effect on the shear band angle.     

Atomistic insights into dislocation-based mechanisms of void growth and coalescence

One of the low-temperature failure mechanisms in ductile metallic alloys is the growth of voids and their coalescence. In the present work we attempt to obtain atomistic insights into the mechanisms underpinning cavitation in a representative metal, namely Aluminum. Often the pre-existing voids in metallic alloys such as Al have complex shapes (e.g. corrosion pits) and the defromation/damage mechanisms exhibit a rich size-dependent behavior across various material length scales. We focus on these two issues in this paper through large-scale calculations on specimens of sizes ranging from 18 thousand to 1.08 million atoms. In addition to the elucidation of the dislocation propagation based void growth mechanism we highlight the observed length scale effect reflected in the effective stress-strain response, stress triaxiality and void fraction evolution. Furthermore, as expected, the conventionally used Gurson's model fails to capture the observed size-effects calling for a mechanistic modification that incorporates the mechanisms observed in our (and other researchers') simulation. Finally, in our multi-void simulations, we find that, the splitting of a big void into a distribution of small ones increases the load-carrying capacity of specimens. However, no obvious dependence of the void fraction evolution on void coalescence is observed.     

Analysis of shear localization in porous materials

A set of constitutive equations are derived based on the authors' lower bound yield loci for porous materials. By using these equations, the conditions for shear localization in porous materials are then investigated and the results are compared with those of Gurson's equations and the finite element analysis. The advantages of the present constitutive equations are fully illustrated.