节理岩体本构模型数值模拟及其验证

为了深入研究节理裂隙岩体的破坏机理,基于Mohr-Coulomb准则,利用(PFC)离散元程序,引入颗粒接触连接本构关系,建立了含孔洞岩样Burgers蠕变模型,以数值模拟的应力-应变曲线与试验曲线吻合作为PFC细观力学参数选取的准则,并利用获得的力学参数对含孔洞节理岩体破坏过程进行数值再现,模拟了不同颗粒单元参数试样...     

混凝土断裂的三维旋转裂缝模型研究

作为一种准脆性材料,混凝土在发生开裂后将表现一定的应变软化特性。基于Bazant提出的钝断裂带力学模型,将其扩展为三维旋转裂缝模型,并推导出相应的三维软化本构,应用于复杂应力路径下混凝土开裂扩展的三维数值模拟。通过对比四点剪切梁试验的数值模拟结果与试验结果,验证了该模型能较好地模拟混凝土混合型开裂行为。     

三维四向编织复合材料界面损伤机理数值分析

考虑界面脱粘表面压应力下摩擦力对材料界面力学性能的影响,建立损伤-摩擦相结合的界面本构模型,编写用户材料子程序VUMAT,实现其在有限元软件ABAQUS中的嵌入。基于周期性胞元分析思想,在单胞模型中纤维束/基体、纤维束/纤维束分界面引入界面单元,结合损伤-摩擦相结合的界面本构模型,建立含界面相三维四向编织复合材料的细观有限元模型。模拟典型载荷下界面损伤的起始和扩展过程,分析界面应力传递和界面破坏机理,研究界面性能对复合材料宏细观力学性能的影响规律,为实现三维四向编织复合材料界面性能优化设计和控制提供参考。      

大幅面板材渐进折弯成形的本构建模及模拟

为了提高大幅面板材成形的模拟精度,在板材折弯平面应变假设条件下,推导出基于Hill各向异性屈服准则的弹塑性本构方程.借助ABAQUS有限元软件本构模块用户子程序接口,通过编程将上述推导的应力-应变本构关系显示表达式嵌入ABAQUS分析平台.以超长大开口半椭圆形工件成形为例,建立了大幅面钢板渐进折弯的三维弹塑性有限元模型,并数值模拟了多道次渐进折弯成形及回弹全过程.模拟效果和工程应用结果表明,与传统的基于平面应力假设的本构关系模型相比,采用平面应变假设的本构关系模型的模拟结果更接近实验值.     

试论镁铝合金高应变率单轴压缩拟合本构关系的代入校核

从实验与数值模拟的结合上,确定与考察了镁铝合金在一维应力高应变率压缩状态下的本构关系。把由一维准静态及动态压缩试验结果所拟合的经验性本构关系代入镁铝合金的SHPB试验的全过程数值模拟,可再现实验测得的反射波波形及透射波波形。数值模拟也表明,采用SHPB试验经典分析所估算的应力、应变、应变率与试件典型微元在试验过程中所经受的应力应变、应变率也基本一致,然而,应变率并非常值。     

混凝土应力三轴比依赖的塑性损伤模型

提出了一种新的损伤塑性本构模型,定义了一个用应力不变量表示的应力三轴比,强调了应力三轴比对塑性屈服的影响。采用本研究设计的专门用于本构校验的计算机软件,对上述模型的应力应变加载曲线进行了数值模拟,数值结果与实验结果吻合得比较好。     

基于能量的弹塑性损伤实用本构模型

建立一个实用的弹塑性损伤本构模型。在有效应力空间采用经验公式计算塑性变形,能将模型塑性部分与损伤部分解耦,降低模型的数值处理复杂性,同时大大简化模型塑性应变的计算。结合不可逆热力学理论,基于损伤能量释放率建立损伤准则,损伤能量释放率由修正后的弹性Helmholtz自由能导出,模型中将弹性Helmholtz自由能分解为应力球量部分和应力偏量部分,将其应力球量部分产生的损伤取为零,同时根据应力状态引入折减系数对其应力偏量部分进行修正,使得模型能较为准确的模拟混凝土材料在双轴加载下的本构行为。将应力张量谱分解为正、负两部分以分别定义材料受拉、受压损伤演化,并采用受拉损伤变量、受压损伤变量分别模拟混凝土材料在拉、压加载下的本构特性。引入一个加权损伤变量使得模型能较准确的反映混凝土材料的“拉-压软化效应”。最后该文给出初步试验验证,证明了该文模型的有效性。     

二维编织C/SiC复合材料非线性损伤本构模型与应用

基于二维编织C/SiC复合材料的基本力学性能试验, 建立了该材料的宏观正交各向异性非线性损伤本构模型。模型中以可检测的应变作为变量, 采用形式简单的函数分别描述了单轴拉伸和剪切加载下的材料损伤演变下的应力-应变关系, 以及卸载状态的刚度变化规律。同时, 考虑了材料的单边效应以及拉压应力状态转换时的损伤钝化行为。将此本构模型编写成UMAT子程序并引入ABAQUS有限元软件, 可以完整描述该材料的加载非线性和卸载线性的应力-应变关系特征, 及其加卸载历史。通过对带孔板的拉伸模拟, 孔边应变分布与试验结果吻合较好, 验证了本构模型的有效性。     

修正GTN模型及其在预测硅钢冷轧边裂中的应用

考虑剪应变对微孔洞损伤演化的影响, 对GTN损伤模型的损伤演化机制进行修正, 建立了适用于不同应力三轴度水平的损伤模型. 结合隐式应力更新算法和显式有限元计算, 采用VUMAT子程序实现了修正GTN模型在有限元软件ABAQUS中的数值计算. 通过模拟纯剪切和剪切-拉伸两组试样的损伤演化和断裂行为, 验证了修正GTN模型在不同应力三轴度承载条件下的有效性. 运用修正GTN损伤模型模拟含边部缺口的带钢在轧制过程中裂纹的萌生和扩展行为, 模拟结果与实验相一致, 表明该模型可有效地用于带钢缺陷在轧制过程中扩展行为的分析和预测. 模拟和实验结果表明, 带钢边部缺口在轧制过程中, 缺口前沿和后沿均会萌生裂纹, 且后沿裂纹扩展更为明显.     

砂土统一本构模型研究及其三维数值实现

基于广义塑性理论与临界状态概念,研究提出了一个统一三维砂土本构模型,通过一组参数实现了砂土由压缩至剪切过程中状态参量的统一表述。基于ABAQUS提供的用户自定义材料子程序UMAT接口,利用Fortran语言编程实现了该三维弹塑性本构关系模型在软件中的二次开发。分别利用Toyoura砂、Fuji River砂以及Tokachi砂的剪切试验数据与数值模拟结果进行对比,结果表明:提出的有限元计算模型可以有效反映加载过程中不同围压和砂土初始密度对应力-应变曲线的影响,能够准确描述密砂的剪胀特性与应变软化特性以及松砂的剪缩特性与应变硬化特性,从而更加真实地反映三维应力状态下土的变形和强度特性。研究成果进一步扩展了ABAQUS在岩土工程中的应用范围,能够为岩土工程领域的数值分析计算提供更加快捷的解决方案。     

Growth and coalescence of penny-shaped voids in metallic alloys

The growth and coalescence of penny-shaped voids resulting from particle fracture is a common damage process for many metallic alloys. A three steps modeling strategy has been followed to investigate this specific failure process. Finite element cell calculations involving very flat voids shielded or not by a particle have been performed in order to enlighten the specific features of a damage mechanism starting with initially flat voids with respect to more rounded voids. An extended Gurson-type constitutive model supplemented by micromechanics-based criteria for both void nucleation and void coalescence is assessed for the limit of very flat voids towards the FE calculations. The constitutive model is then used to generate a parametric study of the effects of the stress state, the microstructure and the mechanical properties on the ductility. Based on these results, a simple closed-form model for the ductility is finally proposed. The main outcomes of this study are that (i) the ductility of metal alloys involving penny-shaped voids is primarily controlled by the relative void spacing; (ii) the definition of an effective porosity in terms of an equivalent population of spherical voids is valid for low particle volume fraction; (iii) the presence of a particle shielding the void does not significantly affect the void growth rates and void aspect evolution; (iv) early fracture by void coalescence can occur under very low stress triaxiality conditions if the particle volume fraction is large enough, explaining that some alloys and composites can fail through a transgranular ductile fracture mode under uniaxial tension condition before the onset of necking; (v) the fracture mechanism moves from void growth controlled to void nucleation controlled when increasing the void nucleation stress, lowering the stress triaxiality, and increasing the initial void aspect ratio.     

Three-dimensional cell model analyses of void growth in ductile materials

Three-dimensional micromechanical models were developed to study the damage by void growth in ductile materials. Special emphasis is given to the influence of the spatial arrangement of the voids. Therefore, periodical void arrays of cubic primitive, body centered cubic and hexagonal structure are investigated by analyzing representative unit cells. The isotropic behaviour of the matrix material is modelled using either v. Mises plasticity or the modified Gurson-Tvergaard constitutive law. The cell models are analyzed by the large strain finite element method under monotonic loading while keeping the stress triaxiality constant. The obtained mesoscopic deformation response and the void growth of the unit cells show a high dependence on the value of triaxiality. The spatial arrangement has only a weak influence on the deformation behaviour, whereas the type and onset of the plastic collapse behaviour are strongly affected. The parameters of the Gurson-Tvergaard model can be calibrated to the cell model results even for large porosity, emphasizing its usefulness and justifying its broad applicability.     

Cell model simulation of void growth in nodular cast iron under cyclic loading

The failure of cast iron under high plastic cyclic strains is controlled by the mechanisms of formation, growth and coalescence of voids. A cell model approach is used to simulate nodular cast iron as a periodic array of loosely bonded spherical inclusions in the matrix material. The models are analyzed by the finite element method under cyclic loading while keeping the stress triaxiality constant. Different types of matrix hardening are used: isotropic, kinematic and combined hardening. The graphite inclusions are simulated by a rigid body. Deformation and void growth are studied in dependence on stress triaxiality and strain range. In most cases after a few cycles a non-symmetric stationary mesoscopic cyclic stress–strain curve is established. The deformation response and the development of the void volume fraction are strongly affected by the value of triaxiality. The void volume is incrementally increasing with each load cycle in a ratcheting manner. The void growth rate depends on the chosen hardening type and is smallest for kinematic hardening. The comparison with simulations in absence of graphite inclusions revealed that void evolution is favored by the inclusions.     

Stress triaxiality effect on void nucleation in ductile metals

The stress triaxiality effect on the strain required for void nucleation by particle‐matrix debonding has been investigated by means of micromechanical modelling. A unit‐cell model considering an elastic spherical particle embedded in an elastic‐plastic matrix was developed to the purpose. Particle‐matrix decohesion was simulated through the progressive failure of a cohesive interface. It has been shown that the parameters of matrix‐particle cohesive interface are correlated with macroscopic material properties. Here, a simple relationship for the maximum cohesive opening at interface failure as a function of material fracture toughness and yield stress has been derived. Results seem to confirm that, increasing stress triaxiality, the strain at which void nucleation is predicted to occur decreases exponentially in a similar way as for fracture strain. This result has substantial implications in modelling of ductile damage because it indicates that if the stress triaxiality is high enough, ductile fracture can occur at plastic strain lower than that necessary to nucleate damage for moderate or low stress triaxiality regime.     

基于微孔贯通细观损伤模型的金属韧性断裂分析

韧性材料断裂过程通常可看作是材料内部微孔洞的形核、扩展及相互贯通的积累。经典的Gurson- Tvergaard (GT)模型能够很好地模拟具有变形均匀、各向同性的孔洞的萌生及扩展过程;但无法模拟由孔洞贯通而引起的局部变形过程,因此需要对其修正,引入相应的孔洞贯通准则。该文采用两种贯通准则对经典GT模型进行修正,即Thomason的塑性极限载荷准则和临界等效塑性应变准则。借助用户自定义程序UMAT将采用这两种贯通准则修正的GT本构关系嵌入至商用有限元软件ABAQUS中,从而可通过对金属材料应力状态和断裂机理的分析控制孔洞的贯通。以一组含有不同缺口根半径的圆棒拉伸试验件为例,分析了该类金属构件自孔洞萌生至最终断裂的整个损伤演化过程,并与试验数据进行了对比,验证了该模型的有效性。该文还讨论了金属断裂过程中应力三轴度对微裂纹萌生与扩展的影响。     

Use of a modified Gurson model approach for the simulation of ductile fracture by growth and coalescence of microvoids under low, medium and high stress triaxiality loadings

In the paper the modified Gurson model is developed for the simulation of damage growth and ductile fracture under low, medium and high stress triaxiality loadings. A new coalescence criterion is introduced based on a simple assumption that singular value of the effective stress triggers the coalescence of microvoids in materials. According to the introduced approach the void coalescence described by means of the modified Gurson model is not only determined by the so-called critical, constant void volume fraction but also by the stress triaxiality ratio. Computational simulations have been carried out for Al 2024–T351 aluminum specimens. In order to find some improvements of micromechanical damage models, two different approaches have been compared for modeling the shear driven microvoid coalescence under low stress triaxiality loadings.     

An exponential ductile continuum damage model for metals

In the frame of continuum damage mechanics an isotropic ductile plastic damage model is derived. The model is based on void damage variable, defined using effective stress concept and thermodynamics. The damage evolution from this model is exponential with equivalent plastic strain as experienced in some low carbon steel like AISI 1015. The damage model is sensitive to stress triaxiality and emulates the damage evolution as recorded in experiments conducted on such metals. The model is validated by comparison with the Rice-Tracey model and other experimental results published in the literature. This model can be used to study the growth and coalescence of micro voids, influence of stress triaxiality on strain to rupture and crack initiation phenomena.     

混凝土变参数等效应变损伤模型

根据混凝土4种不同受力状态下应力-应变曲线的试验数据,建立了以等效应变作为参变量的变四参数等效应变损伤模型.该模型可以较好地描述复杂空间应力状态下混凝土的损伤本构关系,运用该模型计算得到的4种不同受力状态下的全过程应力-应变曲线与试验结果具有较好的一致性.     

Effect of multiaxial stress state on creep damage in ASME T92 welded joint

This study reports the type IV fracture process and the influence of multiaxial stress state in ASME T92 welded joints during creep. The type IV fracture occurs at the fine‐grained heat‐affected zone (ie, FGHAZ), involving void initiation, growth, and coalescence, microcrack occurrence, propagation and extension, and eventual macrocrack with consequent joint failure. The creep damage is not uniformly distributed along the thickness direction in the FGHAZ, and the central part of the welded joint is the most seriously damaged region. The equivalent creep strain is higher at the external surface, but the stress triaxiality is larger in the centre section. Large equivalent creep strain could promote creep void initiation, whereas high hydrostatic pressure and stress triaxiality factor accelerate void growth in the FGHAZ of T92 joints. Besides, reducing groove angle and HAZ width of the joints is recommended to delay the occurrence of type IV cracking because of lower equivalent creep strain and stress triaxiality factor.