基于颗粒离散元法的连接键应变软化模型及其应用

基于颗粒间的有限接触假设,提出了可表述颗粒间力、力矩传递的连接键模型. 为了表征连接键的塑性、损伤及断裂过程,在连接键中引入了考虑应变软化效应的Mohr-Coulomb 准则及最大拉应力准则. 单一连接键的单向拉伸测试及直剪测试表明了上述连接键应变软化模型的计算精度. 研究了颗粒体系的宏观应变能与颗粒平均配位数的对应关系. 通过计算发现,对于二维颗粒体系,当平均配位数为5 时,颗粒体系的宏观应变能与相同参数下连续介质方法(如有限元等) 计算获得的应变能基本一致. 利用上述连接键应变软化模型对岩石的单轴压缩过程进行了模拟,计算结果表明:岩石单轴压缩的应力应变曲线经历了线性上升段、非线性上升段、非线性下降段及缓变段等4 个阶段,并给出了上述4 个阶段与岩石内部损伤破裂状态的内在联系. 计算结果还表明,随着断裂应变的增大,岩石的破裂模式逐渐由拉剪复合型破裂向单一压剪型破裂转化;随着断裂应变的增大,峰值应力及达到峰值应力时的应变均逐渐增大,但峰值时的破裂度及终态时的破裂度将逐渐减小.      

Tension–torsion fracture experiments – Part II: Simulations with the extended Gurson model and a ductile fracture criterion based on plastic strain

An extension of the Gurson model that incorporates damage development in shear is used to simulate the tension–torsion test fracture data presented in Faleskog and Barsoum (2013) (Part I) for two steels, Weldox 420 and 960. Two parameters characterize damage in the constitutive model: the effective void volume fraction and a shear damage coefficient. For each of the steels, the initial effective void volume fraction is calibrated against data for fracture of notched round tensile bars and the shear damage coefficient is calibrated against fracture in shear. The calibrated constitutive model reproduces the full range of data in the tension–torsion tests thereby providing a convincing demonstration of the effectiveness of the extended Gurson model. The model reinforces the experiments by highlighting that for ductile alloys the effective plastic strain at fracture cannot be based solely on stress triaxiality. For nominally isotropic alloys, a ductile fracture criterion is proposed for engineering purposes that depends on stress triaxiality and a second stress invariant that discriminates between axisymmetric stressing and shear dominated stressing.     

基于J-C模型的镁合金MB2动静态拉伸破坏行为

为了研究不同应力状态和应变率条件下镁合金MB2的拉伸破坏行为,利用材料试验机和分离式Hopkinson拉杆(SHTB),对镁合金MB2的光滑及缺口圆柱试件进行了动静态拉伸加载;拟合得到了镁合金MB2的动静态拉伸本构关系,建立了其修正的Johnson-Cook失效破坏准则,并对不同试件的拉伸破坏行为进行了数值模拟;利用SEM对宏观破坏模式对应的微观损伤机理进行了分析。结果表明,随着应力三轴度的增加,镁合金MB2的等效破坏应变先增大后减小,宏观破坏模式由剪切转为正拉断,微观损伤机制由混合断裂转变为韧窝断裂;而随着应变率的增加,等效破坏应变不断减小,破坏模式不发生改变。Johnson-Cook本构关系和修正后的Johnson-Cook失效破坏准则能较好地拟合动态静态拉伸实验结果并预测不同试件的杯锥形破坏特征。     

Ductile fracture: Experiments and computations

Numerous criteria have been developed for ductile fracture (DF) prediction in metal plastic deformation. Finding a way to select these DF criteria (DFCs) and identify their applicability and reliability, however, is a non-trivial issue that still needs to be addressed in greater depth. In this study, several criteria under the categories of ‘uncoupled damage criterion’ and the ‘coupled damage criterion’, including the continuum damage mechanics (CDM)-based Lemaitre model and the Gurson-Tvergaard-Needleman (GTN) model, are investigated to determine their reliability in ductile failure prediction. To create diverse stress and strain states and fracture modes, different deformation scenarios are generated using tensile and compression tests of Al-alloy 6061 (T6) with different sample geometries and dimensions. The two categories of criteria are coded into finite element (FE) models based on the unconditional stress integration algorithm in the VUMAT/ABAQUS platform. Through physical experiments, computations and three industrial case studies, the entire correlation panorama of the DFCs, deformation modes and DF mechanisms is established and articulated. The experimental and simulation results show the following. (1) The mixed DF mode exists in every deformation of concern in this study, even in the tensile test of the round bar sample with the smallest notch radius. A decrease of stress triaxiality (η-value) leads to a reduction in the accuracy of DF prediction by the two DFC categories of DFCs, due to the interplay between the principal stress dominant fracture and the shear-stress dominant factor. (2) For deformations with a higher η-value, both categories of DFCs predict the fracture location reasonably well. For those with a lower or even negative η-value, the GTN and CDM-based criteria and some of the uncoupled criteria, including the C&L, Ayada and Oyane models, provide relatively better predictions. Only the Tresca and Freudenthal models can properly predict the shear dominant fracture. The reliability sequence of fracture moment prediction is thus the GTN model, followed by the CDM-based model and the uncoupled models. (3) The applicability of the DFCs depends on the use of suitable damage evolution rules (void nucleation/growth/coalescence and shear band) and consideration of several influential factors, including pressure stress, stress triaxiality, the Lode parameter, and the equivalent plastic strain or shear stress. These parameters determine the deformation mode (shear dominant or maximum principal stress dominant deformation) and, further, the DF mechanism (dimple fracture/shear fracture/mixed fracture).     

Effect of temperature on the fracture modes of BHW-35 steel under mode-II loading

Finite element method (FEM) has been used to analyze the stress and strain fields and the stress tri-axial levels around the tip of the crack under mode- II loading. The results show that: under mode- II loading, the direction of the maximum tensile stress and that of the maximum tri-axial levels (R _o ) exist at an angle of –75. 3° from the original crack plane; the maximum shear stress andR _o = 0 exist along the original crack plane.Mode- II loading experiment using BHW-35 steel at different temperatures show that there are two kinds of fracture mode, opening mode (or tensile mode) and sliding mode (or shear mode). A decrease in temperature causes the fracture mode to change from shear mode to tensile mode. For BHW-35 steel, this critical temperature is about –90 C. Actually, under any kind of loading mode (mode I . mode II , mode III or mixed mode), there always exist several kinds of potenital fracture modes (for example, opening mode, sliding mode, tearing mode or mixed mode). The effect of temperature under mode- II loading is actually related to the change of the elastic-plastic properties of the material.     

电磁加载下7075铝环的膨胀断裂模式转变研究

利用电磁膨胀环实验技术,研究了7075铝环在2700～8700 s?1拉伸加载应变率下的断裂模式转变现象。实验结果表明:在低应变率下,铝环的断裂受最大正应力控制,发生拉伸断裂;在高应变率下,铝环的断裂受最大剪应力控制,发生剪切断裂;在中应变率下,铝环的断裂同时受最大正应力和最大剪应力控制,为拉伸和剪切同时存在的拉剪混合断裂模式;随着应变率的增加,铝环的破片数量呈先增加、再减小、最后增加的趋势,破片数量变化拐点与断裂模式转变点基本吻合。     

An interacting micro-crack damage model for failure of brittle materials under compression

A model is developed for brittle failure under compressive loading with an explicit accounting of micro-crack interactions. The model incorporates a pre-existing flaw distribution in the material. The macroscopic inelastic deformation is assumed to be due to the nucleation and growth of tensile “wing” micro-cracks associated with frictional sliding on these flaws. Interactions among the cracks are modeled by means of a crack-matrix-effective-medium approach in which each crack experiences a stress field different from that acting on isolated cracks. This yields an effective stress intensity factor at the crack tips which is utilized in the formulation of the crack growth dynamics. Load-induced damage in the material is defined in terms of a scalar crack density parameter, the evolution of which is a function of the existing flaw distribution and the crack growth dynamics. This methodology is applied for the case of uniaxial compression under constant strain rate loading. The model provides a natural prediction of a peak stress (defined as the compressive strength of the material) and also of a transition strain rate, beyond which the compressive strength increases dramatically with the imposed strain rate. The influences of the crack growth dynamics, the initial flaw distribution, and the imposed strain rate on the constitutive response and the damage evolution are studied. It is shown that different characteristics of the flaw distribution are dominant at different imposed strain rates: at low rates the spread of the distribution is critical, while at high strain rates the total flaw density is critical.     

A Study of Large Plastic Deformations in Dual Phase Steel Using Digital Image Correlation and FE Analysis

Large plastic deformation in sheets made of dual phase steel DP800 is studied experimentally and numerically. Shear testing is applied to obtain large plastic strains in sheet metals without strain localisation. In the experiments, full-field displacement measurements are carried out by means of digital image correlation, and based on these measurements the strain field of the deformed specimen is calculated. In the numerical analyses, an elastoplastic constitutive model with isotropic hardening and the Cockcroft–Latham fracture criterion is adopted to predict the observed behaviour. The strain hardening parameters are obtained from a standard uniaxial tensile test for small and moderate strains, while the shear test is used to determine the strain hardening for large strains and to calibrate the fracture criterion. Finite Element (FE) calculations with shell and brick elements are performed using the non-linear FE code LS–DYNA. The local strains in the shear zone and the nominal shear stress-elongation characteristics obtained by experiments and FE simulations are compared, and, in general, good agreement is obtained. It is demonstrated how the strain hardening at large strains and the Cockcroft–Latham fracture criterion can be calibrated from the in-plane shear test with the aid of non-linear FE analyses. An erratum to this article can be found at     

中子辐照金属材料的脆化模型研究

金属材料的辐照脆化问题一直以来都是核能安全领域亟待解决的关键问题之一.为了更准确地预测金属材料的辐照脆化行为,基于Johnson-Cook本构模型,将未辐照金属材料的断裂真应力取作辐照材料的断裂真应力,建立了通过辐照退火态金属材料屈服强度就能够预测其整个真应力$\!$-$\!$-$\!$应变曲线,以及断裂真应变的辐照脆化模型.实验研究了不同中子剂量辐照退火态高纯铝的准静态拉伸真应力$\!$-$\!$-$\!$应变曲线、断裂真应力和断裂真应变随辐照剂量的变化规律.结果表明,辐照剂量越高,高纯铝的屈服强度越高,断裂真应变越低,但断裂真应力几乎不变.通过TEM显微分析获得了高纯铝内部辐照缺陷的尺寸和数密度随辐照剂量的变化规律,结果表明,辐照剂量越高,孔洞的尺寸和数密度越高,但位错环尺寸和数密度始终很小,难以准确统计.由辐照高纯铝实验数据拟合得到了辐照脆化模型所需参数,并检验了该模型的预测效果.结果表明,无论是通过实验还是显微分析得到辐照高纯铝的屈服强度,模型的预测结果均能够与实验结果较好地吻合,且模型对退火态高纯铝临界中子剂量的预测值也与文献结果一致.      

褶皱与晶界耦合作用对石墨烯断裂行为的影响

由CVD方法制备的石墨烯含有大量的晶界,通常还带有许多褶皱,本文通过分子动力学方法研究了具有褶皱和晶界的石墨烯 平面拉伸断裂行为,结果显示,在垂直晶界方向,褶皱能够显著提高小角度双晶石墨烯的断裂应力,断裂应力增幅最 大约为50%,褶皱对断裂应力的影响随晶界角的增大减弱,导致双晶石墨烯断裂应力对晶界角不敏感,只略低于单晶石墨 烯,和实验结果完全吻合;在沿晶界方向,褶皱对双晶石墨烯断裂应力影响不明显. 另外,褶皱可以显著提高双晶石墨烯的断 裂应变,增幅最大约为100%. 增强机制归纳主要如下:通过面外变形,褶皱可以部分释放晶界5-7环中C---C键的预拉伸变形, 提高双晶石墨烯的断裂应力;褶皱可以降低相邻5-7环之间相互作用,导致断裂应力对晶界角不敏感;在拉伸作用下,褶皱被部分 拉平,这可以显著降低C---C键面内拉伸变形,导致断裂应变显著增大. 本研究为准确理解多晶石墨烯断裂行为提供重要帮助.      

A crack on the grain boundary — An analysis based on crystal plasticity theory

The influence of the mismatch of the lattice orientation on the deformation and stress fields of a crack located on the grain boundary is studied by means of the finite-element analysis taking account of finite deformatio and finite lattice rotation. The plane strain calculations for an fcc crystal subjected to mode I loading are performed on the basis of the crystalline plasticity described by a planar three-slip model. For the crack-tip shapes and the dominant deformation modes on slip systems, results of all the cases analysed here are in qualitative agreement with the earlier analytical and numerical solutions. Our results indicate that the lattice orientation difference may greatly influence the shear stress along the grain boundary which is related to grain-boundary sliding, while the normal stress along the grain boundary, which may induce cleavage fracture, is virtually insensitive to it. The influence of the lattice orientations on the crack-tip fields is also investigated under small-scale-yielding conditions and the comparison with the results of finite deformation is made.     

Effect of stress-state and spacing on voids in a shear-field

Unit cell model analyses are carried out for a material with a periodic array of voids, subject to shear loading. Thus the focus is on ductile fracture in conditions of low stress triaxiality. It has been shown recently that voids in shear are flattened out to micro-cracks, which rotate and elongate until interaction with neighboring micro-cracks gives coalescence, so that the failure mechanism is very different from that under tensile loading. In the present studies the plane strain unit cell has fully periodic boundary conditions, so that any combination of the stress components in the overall average stress state can be prescribed. This also allows for studies of the effect of different initial void spacing in the two in-plane coordinate directions. The stress states considered are essentially simple shear, with various levels of tensile stresses or compressive stresses superposed, i.e. low positive stress triaxiality or even negative stress triaxiality. For high aspect ratio unit cells a clear localization band is found inside the cell, which actually represents several parallel bands, due to periodicity. In the materials represented by a low aspect ratio unit cell localization would also occur after that the maximum shear stress has been passed, but this is not shown when periodicity is enforced. The effect of superposed tensile or compressive stresses is found to be bigger for high aspect ratio unit cells than for low aspect ratios.     

45钢柱壳爆炸膨胀断裂的SPH模拟分析

金属柱壳爆炸膨胀断裂存在拉伸、剪切及拉剪混合等多种断裂模式,目前其物理机制及影响因素还不清晰。本文中采用光滑粒子流体动力学方法（smoothed particle hydrodynamics, SPH）对45钢柱壳在JOB-9003及RHT-901不同装药条件下的外爆实验进行了数值模拟,探讨柱壳在不同装药条件下发生的剪切断裂、拉剪混合断裂模式及其演化过程,模拟结果与实验结果一致。SPH数值模拟结果表明:在爆炸加载阶段,随着冲击波在柱壳内、外壁间来回反射形成二次塑性区,沿柱壳壁厚等效塑性应变演化呈凸形分布,壁厚中部区域等效塑性应变较内、外壁大;在较高爆炸压力（JOB-9003）作用下,柱壳断裂发生在爆轰波加载阶段,损伤裂纹从塑性应变积累较大的壁厚中部开始沿剪切方向向内、外壁扩展,形成剪切型断裂模式;而在RHT-901空心炸药加载下,虽然裂纹仍从壁厚中部开始沿剪切方向扩展,但随后柱壳进入自由膨胀阶段,未断区域处于拉伸应力状态,柱壳局部发生结构失稳,形成类似“颈缩”现象,裂纹从剪切方向转向沿颈缩区向外扩展,呈现拉剪混合断裂模式。拉伸裂纹占截面的比例与柱壳结构失稳时刻相关。可见,柱壳断裂演化是一个爆炸冲击波与柱壳结构相互作用的过程,不能简单将其作为一系列膨胀拉伸环处理。     

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

NONLINEAR UNIFIED STRENGTH MODEL OF GEOMATERIALS

将材料的破坏归结为剪切破坏,每种材料对应于特定的剪切滑动面,抗剪强度为滑动面上正应力的函数,基于不同材料的强度特性将一系列的剪切滑动面统一起来,建立了岩土材料的非线性统一强度模型.非线性统一强度模型的滑动面为β应力空间内的等倾面,在β应力空间内的强度面为圆锥面;在普通应力空间内的强度面为一系列连续光滑、外凸的锥面,在偏平面上强度曲线涵盖了从下限Matsuoka-Nakai曲线到上限Drucker-Prager圆之间的所有区域,子午面上强度线为直线.非线性统一强度模型只有3个材料参数,参数都具有明确的物理意义,通过与国内外学者已取得的岩土类材料真三轴强度试验结果的比较,表明模型可适用于多种类型的材料,并合理描述其非线性强度特性.     

Modification of the fracture criterion for V-shaped notches (plane problem). Relationship between toughness and strength and structural parameters

A fracture criterion of the type of the Neuber-Novozhilov criterion is proposed to describe the fracture in the vicinity of the tip of a V-shaped notch under tensile and shear loading. In the proposed criterion, the limits of averaging of the stresses along the notch axis depend on the presence, location, and size of the initial defects in the material. The crystal lattice parameter of the initial material is chosen for the characteristic linear size. For a V-shaped notch subjected to tension and shear, simple equations are obtained that relate the stress intensity factors for the modified singularity coefficients, the singularity coefficients themselves, and the theoretical tensile and shear strengths of a single crystal of the material taking into account the damage to the material in the vicinity of the notch tip. The equations obtained allow a passage to the limit from a notch to a crack. It is shown that the classical critical stress intensity factor used in the strength analysis of cracked solids is not a material constant.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 106–115, January–February, 2005.     

Interaction between hydraulic and natural fractures

The results of a numerical simulation of the interaction between a hydraulic fracture propagating under plane-strain conditions in an infinite impermeable elastic medium and an already existing natural fault are presented. The fracture is created by an incompressible Newtonian fluid injected at a constant flow rate from a source located at its center; the behavior of the fault is described by a dry friction model of the Mohr-Coulomb type. The results obtained are in agreement with the available studies in this field of research.     

岩土材料的非线性统一强度模型

将材料的破坏归结为剪切破坏,每种材料对应于特定的剪切滑动面,抗剪强度为滑动面上正应力的函数,基于不同材料的强度特性将一系列的剪切滑动面统一起来,建立了岩土材料的非线性统一强度模型.非线性统一强度模型的滑动面为β应力空间内的等倾面,在β应力空间内的强度面为圆锥面;在普通应力空间内的强度面为一系列连续光滑、外凸的锥面,在偏平面上强度曲线涵盖了从下限Matsuoka-Nakai曲线到上限Drucker-Prager圆之间的所有区域,子午面上强度线为直线.非线性统一强度模型只有3个材料参数,参数都具有明确的物理意义,通过与国内外学者已取得的岩土类材料真三轴强度试验结果的比较,表明模型可适用于多种类型的材料,并合理描述其非线性强度特性.      

PMMA膨胀环动态拉伸碎裂实验研究

在强动载作用下, 脆性材料的碎裂问题是一个重要的研究课题, 而脆性材料在冲击拉伸载荷下的力学行为的实验研究相对较匮乏. 提出了一种动态拉伸断(碎)裂的液压膨胀环实验技术, 可用于准脆性/脆性材料的动态拉伸. 利用该技术对有机玻璃(PMMA)圆环试件进行了不同膨胀速度下的动态碎裂实验研究. 从回收碎片的断口形貌和碎片内部残余裂纹观察可知试件的破碎由环向拉伸应力造成, 碎片断口处发出的稀疏波会将周围的拉伸应力卸载, 从而抑制其他裂纹的进一步发展. 利用超高速相机记录了试件的膨胀碎裂过程, 利用DISAR激光速度干涉仪获得了试件外表面粒子的径向膨胀速度历史, 通过试件上的应变片获得了试件的应变历史和断裂应变. 实验结果表明: 在拉伸应变率150\mathrm{~}500s^{-1}范围, 材料的动态断裂应变低于准静态加载下的断裂应变, 体现出“动脆”现象; 随着加载应变率的提高, PMMA 材料的碎片尺寸减小; 无量纲化的PMMA圆环的平均碎片尺寸介于韧性碎裂模型和脆性碎裂模型的预测数值之间, 反映出材料的准脆性特性.