多轴载荷作用下马氏体耐热钢中微孔洞演化行为的模拟研究

微孔洞演化包括孔洞成核、生长和聚合三个阶段,是影响金属材料韧性断裂的重要因素.为了分析P91马氏体耐热钢中的塑性滑移对微观孔洞扩展的影响,论文提出了一种基于晶体塑性有限元的微观力学计算模型,量化了应力三轴度、洛德参数和晶体取向对微孔洞演化行为的影响.结果表明,在相对较高应力三轴度条件下,随着应力三轴度的增大,含孔洞马氏体块的等效应力-应变响应呈现出快速软化的特征,同时孔洞体积分数随着等效应变的增加而快速增加.对于给定的应力三轴度,胞元的聚合应变在[111]取向时最小,在[110]取向时最大.孔洞聚合开始时,低应力三轴度下孔洞形状趋向于椭球状,而较高应力三轴度时孔洞横向略鼓.在一定的应力三轴度和洛德参数范围内,在孔洞聚合和孔洞坍塌两种胞元失效状态之间存在着一个条带状过渡状态,在[100]晶体取向时,当洛德参数L=-1时条带最宽.论文揭示了P91马氏体耐热钢中微孔洞演化的基本机制,分析了晶体取向、应力三轴度和洛德参数对微孔洞演化的影响.这些发现为P91马氏体耐热钢的韧性损伤模型的进一步发展提供了微观理论基础.     

三轴压缩岩石细观损伤扩展特性CT实时检测

利用作者最新研制的与CT（computerized tomography)机配套的专用加载设备,完成了三轴压缩荷载作用下岩石破坏全过程的细观损伤扩展规律的实时CT检测试验。得到了在不同荷载作用下岩石中微孔洞被压密、微裂纹萌生、分叉、发展、断裂、破坏、卸载等各个阶段清晰的CT图像。对得到的CT数、CT图像等数据进行了分析,引入了初始损伤影响因子和闭合影响系数,定义了一个基于CT数的损伤变量,基于岩石细观损伤演化CT试验的结果,给出了岩石应力损伤门槛值,将岩石应力应变全过程曲线分为5段,得到了岩石损伤扩展的初步规律。     

半结晶聚合物损伤演化的实验表征与数值模拟

损伤本构模型对研究材料的断裂失效行为有重要意义, 但聚合物材料损伤演化的定量表征实验研究相对匮乏. 通过4种高密度聚乙烯(high density polythylene, HDPE)缺口圆棒试样的单轴拉伸实验获得了各类试样的载荷-位移曲线和真应力-应变曲线, 采用实验和有限元模拟相结合的方法确定了HDPE材料不同应力状态下的本构关系, 并建立了缺口半径与应力三轴度之间的关系;采用两阶段实验法定量描述了4种HDPE试样单轴拉伸过程中的弹性模量变化, 并建立了基于弹性模量衰减的损伤演化方程, 结合中断实验和扫描电子显微镜分析了应力状态对HDPE材料微观结构演化的影响. 结果表明缺口半径越小, 应力三轴度越大, 损伤起始越早、演化越快; 微观表现为: 高应力三轴度促进孔洞的萌生和发展, 但抑制纤维状结构的产生;基于实验和有限元模拟获得的断裂应变、应力三轴度、损伤演化方程等信息提出了一种适用于聚合物的损伤模型参数确定方法, 最后将本文获得的本构关系和损伤模型用于HDPE平板的冲压成形模拟, 模拟结果与实验结果吻合良好.      

含裂隙材料的空洞化损伤

本文研究了含有微小裂隙的韧性材料中细观损伤的演化。随外加应力的增大,在裂隙周围的基体中微小空洞不断地萌生并扩展。在一些情况下,由此而形成的内部损伤与仅有大小空洞的损伤有显著的不同特点。结果还表明,具有细观尺度的短裂纹,其损伤作用不宜用裂纹长度作标志。文中最后提出一个材料韧性断裂的判据。     

三轴应力场中不同形状孔洞的长大及其新模型

对不同形状孔洞在从光滑试样到裂纹试样这样广泛三轴应力场中的长大规律,本文通过控制体胞宏观应力三维度的方法进行了精确的有限元分析,计算结果表明：（1）孔洞的体积改变和形状变化是孔洞演化的两种基本机制,在不同的三轴应力场中,这两种机制的作用不同;（2）现有模型对孔洞长大规律的描述是不准确的,由它们得到的临界孔洞扩张比参数HGC与临界孔洞体积分数fc不具备一一对应关系,因此不以很好地反映也洞的实际扩张。在此基础上,提出了一个描述孔洞长大的新模型,与四种常用的现有模型相比,该模型不仅能更好地描述不同三轴应力场中孔洞的长大,而且能反映不同应力三维度水平下材料破坏模式的变化。     

考虑三轴约束时孔洞的聚合机理及有效能量准则

通过体胞分析方法,对不同状孔洞在从光滑试样到裂纹试样的三轴应力场中的聚合机理进行了较精解的有限元分析,计算结果表明：（１）孔洞的相互靠近和横向扩展是导致相邻孔洞发生内颈缩聚合的两种基本机制,在应力三维度Ｒσ等于１．２５附近,这两种机制发生较明显的变化。（２）单纯以孔洞体积分数ｆＣ概念为基础的材料破坏参数一般敏感于应力三维度,不能很好地预报不同三轴应力场中材料的破坏,在此基础上,提出了描述孔洞聚合的     

混凝土三维细观模型的建模方法与力学特性分析

根据混凝土材料的细观组成和结构特点，基于三维Voronoi图形提出了一种简单高效的混凝土细观模型生成方法，利用塑性损伤模型对该细观模型进行了单、多轴应力状态下的准静态分析以及SHPB动态有限元分析。结果表明，数值模拟得到的应力应变曲线和破坏模式与实验结果基本吻合，本文中提出的混凝土三维细观模型可较好地模拟混凝土的静、动态力学特性，为进一步从细观力学角度研究混凝土损伤演化规律和破坏机理提供了模型基础。     

颗粒缺陷相互作用下复合材料的细观损伤模型

含尖角的非椭球颗粒附近应力集中较大,诱导缺陷形成裂纹是材料损伤的重要来源.对于强界面颗粒,大刚度颗粒诱导裂纹向基体中扩展形成近似平面片状裂纹,认为诱导裂纹受颗粒应力附近应力场控制,基于有效自洽理论建立了材料细观损伤模型,得到了单向拉伸下的损伤演化,并分析了颗粒形状、尺寸、颗粒性能以及颗粒与初始缺陷相对位置等因素对材料损伤的影响.结果表明,非椭球颗粒更易诱发裂纹,同样外载应力下,损伤程度更大,含非椭球颗粒材料强度更低;含扁平型的颗粒材料裂纹损伤过程更加明显并且材料强度更大;提高颗粒刚度和含量能够增大材料强度.材料中存在尺寸过大或过小的初始裂纹时材料损伤过程不明显.     

层裂损伤的细观数值分析

采用体胞模型的分析方法推导了材料在弹塑性变形阶段的孔洞增长方程。假设所有孔洞的内外半径之比相同,用数值方法定性分析了材料在层裂损伤过程中孔洞数密度分布的变化。通过分析孔洞数密度分布、孔洞体积累积百分比、不同大小孔洞所占体积份额的计算结果,指出初始损伤对损伤演化有直接影响。     

复杂应力下脆性岩石材料的微裂纹损伤特性

基于无限大体深埋平片状椭圆形裂纹的变形场及其扩展条件,分别推导了脆性岩石材料内部具有任意空间取向的单个张开型或闭合型椭圆形微裂纹及其扩展引起的附加柔度张量;考虑微裂纹系统对材料变形的影响,引入了概率密度函数,得到了任意应力状态下脆性岩石材料的宏细观损伤模型;分析了椭圆形微裂纹的短长轴比率对材料损伤的影响.计算结果显示:材料损伤随着短长轴比率的增大而增大,短长轴比率对材料损伤的影响会随着载荷的增大而提高.将本文模型应用于混凝土的单轴拉伸和花岗岩的单轴压缩,结果表明本文模型能够对实验现象给予很好地解释.     

A molecular dynamics study of void growth and coalescence in single crystal nickel

Molecular dynamics simulations using Modified Embedded Atom Method (MEAM) potentials were performed to analyze material length scale influences on damage progression of single crystal nickel. Damage evolution by void growth and coalescence was simulated at very high strain rates (10⁸–10¹⁰/s) involving four specimen sizes ranging from ≈5000 to 170,000 atoms with the same initial void volume fraction. 3D rectangular specimens with uniform thickness were provided with one and two embedded cylindrical voids and were subjected to remote uniaxial tension at a constant strain rate. Void volume fraction evolution and the corresponding stress–strain responses were monitored as the voids grew under the increasing applied tractions.The results showed that the specimen length scale changes the dislocation pattern, the evolving void aspect ratio, and the stress–strain response. At small strain levels (0–20%), a damage evolution size scale effect can be observed from the damage-strain and stress–strain curves, which is consistent with dislocation nucleation argument of Horstemeyer et al. [Horstemeyer, M.F., Baskes, M.I., Plimpton, S.J., 2001a. Length scale and time scale effects on the plastic flow of FCC metals. Acta Mater. 49, pp. 4363–4374] playing a dominant role. However, when the void volume fraction evolution is plotted versus the applied true strain at large plastic strains (>20%), minimal size scale differences were observed, even with very different dislocation patterns occurring in the specimen. At this larger strain level, the size scale differences cease to be relevant, because the effects of dislocation nucleation were overcome by dislocation interaction.This study provides fodder for bridging material length scales from the nanoscale to the larger scales by examining plasticity and damage quantities from a continuum perspective that were generated from atomistic results.     

Modeling of crack growth in ductile solids: a three-dimensional analysis

A population of several spherical voids is included in a three-dimensional, small scale yielding model. Two distinct void growth mechanisms, put forth by [Int. J. Solids Struct. 39 (2002) 3581] for the case of a two-dimensional model containing cylindrical voids, are well contained in the model developed in this study for spherical voids. A material failure criterion, based on the occurrence of void coalescence in the unit cell model, is established. The critical ligament reduction ratio, which varies with stress triaxiality and initial porosity, is used to determine ligament failure between the crack tip and the nearest void. A comparison of crack initiation toughness of the model containing cylindrical voids with the model containing spherical voids reveals that the material having a sizeable fraction of spherical voids is tougher than the material having cylindrical voids. The proposed material failure determination method is then used to establish the fracture resistance curve (J–R curve) of the material. For a ductile material containing a small volume fraction of microscopic voids initially, the void by void growth mechanism prevails, which results in a J–R curve having steep slope. On the other hand, for a ductile material containing a large volume fraction of initial voids, the multiple voids interaction mechanism prevails, which results in a flat J–R curve. Next, the effect of T-stress on fracture resistance is examined. Finally, nucleation and growth of secondary microvoids and their effects on void coalescence are briefly discussed.     

Formulation of a damage internal state variable model for amorphous glassy polymers

The following article proposes a damage model that is implemented into a glassy, amorphous thermoplastic thermomechanical inelastic internal state variable framework. Internal state variable evolution equations are defined through thermodynamics, kinematics, and kinetics for isotropic damage arising from two different inclusion types: pores and particles. The damage arising from the particles and crazing is accounted for by three processes of damage: nucleation, growth, and coalescence. Nucleation is defined as the number density of voids/crazes with an associated internal state variable rate equation and is a function of stress state, molecular weight, fracture toughness, particle size, particle volume fraction, temperature, and strain rate. The damage growth is based upon a single void growing as an internal state variable rate equation that is a function of stress state, rate sensitivity, and strain rate. The coalescence internal state variable rate equation is an interactive term between voids and crazes and is a function of the nearest neighbor distance of voids/crazes and size of voids/crazes, temperature, and strain rate. The damage arising from the pre-existing voids employs the Cocks–Ashby void growth rule. The total damage progression is a summation of the damage volume fraction arising from particles and pores and subsequent crazing. The modeling results compare well to experimental findings garnered from the literature. Finally, this formulation can be readily implemented into a finite element analysis.     

Multiscale modeling of the damaged plastic behavior and failure of Al/SiCp composites

We studied the tensile behavior and damage of an aluminium X2080 reinforced with different volume fractions of silicon carbide particles. The main damage mechanism is particle failure. Regions of the matrix adjacent to broken particles are sites with high hydrostatic tension and hence the nucleation of cavities is expected. Using J integral and HRR modified solution it is possible to calculate the growth of these voids. Macroscopic failure is governed by a critical volume fraction of voids. The originality of this work lies in the modeling of the composite using a micromechanical approach.     

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.     

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.     

单晶与纳米多晶锡层裂的分子动力学研究

低熔点金属的层裂是目前延性金属动态断裂的基础科学问题之一。采用非平衡态分子动力学方法模拟了冲击压力在13.5～61.0 GPa下单晶和纳米多晶锡的经典层裂和微层裂过程。研究结果表明：在加载阶段,冲击速度不影响单晶模型中的波形演化规律,但影响纳米多晶模型中的波形演化规律,其中经典层裂中晶界滑移是影响应力波前沿宽度的重要因素;在单晶模型中,经典层裂和微层裂中孔洞成核位置位于高势能处;在纳米多晶模型中,经典层裂中的孔洞多在晶界（含三晶界交界处）处成核,并沿晶定向长大,产生沿晶断裂,而微层裂中孔洞在晶界和晶粒内部成核,导致沿晶断裂、晶内断裂和穿晶断裂;孔洞体积分数呈现指数增长,相同冲击速度下单晶和纳米多晶Sn孔洞体积分数变化规律一致;经典层裂中孔洞体积分数曲线的两个转折点分别表示孔洞成核与长大的过渡和材料从损伤到断裂的灾变性转变。