爆炸和冲击载荷下金属材料及结构的动态失效仿真

通过数值模拟研究爆炸冲击载荷下金属材料和结构的动态失效规律对表征爆炸冲击毁伤效应及设计新型抗冲击结构具有重要意义.强动载下金属材料失效涉及材料大变形、热力耦合、材料状态变化等多个复杂物理过程,给数值仿真带来了极大挑战,其中包括裂纹、剪切带等复杂失效模式的几何描述、动态失效准则的确定、塑性与损伤耦合演化的描述等问题.针对这些挑战性问题,基于能量变分建立描述金属动态失效过程的热弹塑性相场理论和计算模型,实现了断裂与剪切带失效模式的统一描述,并提出了其显式有限元高效求解策略.进一步将该模型应用于爆炸冲击载荷下金属脆韧失效模式转变、绝热剪切带（ASBs）自组织及冲击波作用下薄壁圆盘失效形式转变三个典型金属动态失效问题,验证了理论模型的准确性及计算模型的稳健性.该工作为后续开展基于仿真的爆炸冲击毁伤评估及防护结构设计研究奠定了基础.     

混凝土材料裂纹尖端损伤带尺寸特性*

本文采用作者研究各向异性材料裂纹尖端塑性区特性的分析方法,研究了混凝土材料裂纹尖端损伤带尺寸特性.基于单轴应力假定,考虑了混凝土材料的非线性软化特性等影响,提出用两种分析模型:即能量模型和断裂模型来分别讨论.最后,还对这些模型作了比较.     

弹性-应变软化粘塑性材料反平面剪切动态扩展裂纹尖端的渐近解*

本文首先给出了一种用于描述材料软化,并存在有粘塑性的材料模型.用这种模型对反平面剪切型动态扩展状态下,裂纹尖端的弹粘塑性场进行了渐近分析,给出了弹性－应变软化粘塑性材料反平面剪切动态扩展裂纹尖端的渐近解方程.分析结果表明,在裂纹尖端应变具有(ln(R/r))1/(n＋1)的奇异性,应力具有(ln(R/r))-n/(n＋1)的奇异性.从而本文揭示了应变软化粘塑性材料反平面剪切动态扩展裂纹尖端的渐近行为.     

韧性材料中的一个新的动态孔洞增长模型

本文提出了一个可用来描述韧性材料动态断裂过程的新模型.基体是由具有过应力形式的弹-粘塑性本构关系的材料构成.文中给出了在常等效应变率下的含孔洞材料的动态加载面及其相应的近似解析表达式.当材料单元经受球对称载荷作用时,本文的模型将退化为Carroll-Holt-Johnson模型.当材料的应变率敏感参数趋于零时,本文模型将退化为Gurson模型.本文的研究还表明,在一般情况下正交性法则是不成立的.最后,将本文模型与近年来Pan,Saje,Needleman和Perzyna等人所提出的含孔洞材料的粘塑性本构模型进行了比较.     

高温合金GH4169的失效特性及失效模型研究

对GH4169高温合金开展了不同应力三轴度（-0.33~0.33）、不同应变率（0.001~5 000 s-1）、不同温度（293~1 073 K）条件下的材料性能试验.基于Johnson-Cook失效模型的框架,研究了Johnson-Cook(JC)失效模型及已有文献提出的修改形式中应力三轴度项拟合结果的不确定性及应变率对失效应变的线性关系描述局限性问题,通过提出的特定参数确定方法与耦合应力三轴度的应变率效应指数函数,建立了一种唯象修正的失效模型.基于GH4169高温合金的试验结果,标定了修正的失效模型与JC模型中各个参数.结果表明:在不同应力三轴度下,GH4169的失效应变表现出不同的应变率效应;与传统的JC模型相比,修正的失效模型更能够较好地描述GH4169的失效行为;同时能够保证失效应变的非负性.     

简单剪切振荡现象及弹塑性本构的限制条件

基于在真应力空间刻划弹塑性物质的强化、软化和理想塑性特性,本文以刚塑性随动强化模型为例说明产生简单剪切振荡现象,是与模型在简单剪切变形下,其强化和软化特性发生交替变化的现象有关。为使弹塑性本构模型更符合实际,要求它们必须满足如下条件:即对任意弹塑性加载变形过程,本构模型所给出的应力应该是非振荡的,所描述的强化和软化特性,不存在从软化阶段至强化阶段的过渡。     

纵波在率相关塑性软化材料中的传播

本文所推导出的率相关软化材料在塑性软化状态下所满足的运动方程是具有双曲型微分方程性质,从而证实了由于应变率效应,应力波在软化状态下仍将传播.这与率无关理论有着质的不同.从上述方程求解出半无限长杆端部受冲击载荷的瞬态解,结果表明,软化状态下的塑性波波速、软化区界面传播的距离与材料的率敏感性和软化特征有关.     

双势积分算法在非关联材料中的应用

在双势理论的框架下,根据材料自由能形式,材料可以被划分为显式标准材料和隐式标准材料.以经典的非关联D-P模型为例,对其本构锥体进行了描述,并引入了一对对偶锥体.证明了在对偶锥体的描述下,不仅能满足非关联D-P模型自身本构关系,其应力和塑性应变也能满足隐式流动表达.结合双势理论和D-P模型自身的本构特点,推导出了非关联D-P模型率形式弹性状态下、率形式塑性状态下、增量形式弹性状态下、增量形式塑性状态下和增量形式弹塑性状态下的双势函数,从而得到了非关联D-P模型的双势积分算法.通过数值模拟算例验证了双势积分算法的准确性和稳定性.     

求解层间界面反平面剪切破坏的剪切梁模型(Ⅰ)——基本特性

在反平面剪切载荷及侧压力共同作用下引起的裂纹及裂纹扩展导致的层间界面失效,是岩土工程层间界面及砌体结构中界面层上典型的失效方式.运用弹性力学和断裂力学的理论原理,提出了能够反映上述层间界面断裂失效问题力学特性的剪切梁模型.文中采用具有应力软化特性的“粘性裂纹”(内聚力裂纹)模型来表述层间裂纹前方损伤过程区的本构行为.对通过粘性层结合在一起的两个弹性板,在反平面剪切载荷及侧压力共同作用下的力学行为作了解析分析计算,研究了层间界面裂纹扩展规律.     

基于有限元的金属基短纤维复合材料MMC的一种统计蠕变模型

在有限元分析的基础上建立了一个单向应力状态下金属基短纤维复合材料(MMC)的统计蠕变模型.首先建立细胞模型并进行有限元分析,得到了单向应力状态下材料细观尺寸及载荷方向对宏观蠕变响应的影响规律.通过在细胞模型中增加一界面层(考虑材料特性和厚度)来研究基体和纤维的界面对MMC宏观蠕变响应的影响.基于细胞模型的数值结果,提出了一适用于纤维平面随机分布的随机统计模型,该模型考虑了纤维的断裂.通过试验获得纤维的统计分布规律.分析结果表明随机统计模型可以满意地描述试验结果.进一步讨论了材料细观尺寸,纤维的断裂特性以及界面层的材料特性和厚度对MMC宏观蠕变响应的影响.     

Simulating quasi-brittle failures including damage-induced softening based on the mechanism of stress redistribution

There have been extensive analytical and numerical studies on crack propagation and softening behavior in quasi-brittle materials, such as concrete, rock and bone, based on the smeared crack model or the discrete crack model, each of which has advantages and disadvantages. It is well known that stress redistribution has an important effect on cracking behavior and may be taken as the physics mechanism. Using the principle of superposition, the theoretical explanation is provided for stress redistribution during the progressive fracture process. Although this explanation deals with problems of linear elastic fracture mechanics, it is a clue of modeling softening behavior. Then, it is shown that the above stress redistribution based method satisfies the necessary condition that the specimen returns to the unstressed state whenever all the external forces are removed. Finally, a new model for simulating quasi-brittle softening is established via a kind of "sub-element" methodology. Numerical examples show the effectiveness of the above subelement method in modeling the behavior of block media with quasi-brittle block-block interfaces.     

Interlaminar Crack Propagation Analysis of ENF Specimens Based on the Cohesive Zone Model北大核心CSCD

Based on the classical laminated plate theory and the cohesive zone model, a theoretical model for general delamination cracked laminates was established for crack propagation of pure mode Ⅱ ENF specimens. Compared with the conventional beam theory, the proposed model fully considered the softening process of the cohesive zone and introduced the nonlinear behavior of ENF specimens before failure. The predicted failure load is smaller than that under the beam theory and closer to the experimental data in literatures. Compared with the beam theory with only fracture toughness considered, the proposed model can simultaneously analyze the influences of the interface strength, the fracture toughness and the initial interface stiffness on the load-displacement curves in ENF tests. The results show that, the interface strength mainly affects the mechanical behavior of specimens before failure, but has no influence on crack propagation. The fracture toughness is the main parameter affecting crack propagation, and the initial interface stiffness only affects the linear elastic loading stage. The cohesive zone length increases with the fracture toughness and decreases with the interface strength. The effect of the interface strength on the cohesive zone length is more obvious than that of the fracture toughness. When the adhesive zone tip reaches the half length of the specimen, the adhesive zone length will decrease to a certain extent. Copyright ©2022 Applied Mathematics and Mechanics. All rights reserved.     

The phase-field model with an auto-calibrated degradation function based on general softening laws for cohesive fracture

Phase-field models have become popular to simulate cohesive failure problems because of their capability of predicting crack initiation and propagation without additional criteria. In this paper, a new phase-field damage model coupled with general softening laws for cohesive fracture is proposed based on the unified phase-field theory. The commonly used quadratic geometric function in the classical phase-field model is implemented in the proposed model. The modified degradation function related to the failure strength and length scale is used to obtain the length scale insensitive model. Based on the analytical solution of a 1-D case, general softening laws in cohesive zone models can be considered. Parameters in the degradation function can be calibrated according to different softening curves and material properties. Numerical examples show that the results obtained by the proposed model have a good agreement with experimental results and the length scale has a negligible influence on the load-displacement curves in most cases, which cannot be observed in classical phase-field model.     

Damage and fracture behavior in dynamic tension tests: Modelling and numerical simulations

The continuum damage model is based on a general thermodynamic framework for the modeling of rate and temperature dependent behavior of anisotropically damaged elastic-plastic materials subjected to fast deformation. The introduction of damaged and fictitious undamaged configurations allows the definition of damage tensors and the corresponding free energy functions lead to material laws affected by damage and temperature. The damage condition and the corresponding damage rule strongly depend on stress triaxiality. Furthermore, the rate and temperature dependence is reflected in a multiplicative decomposition of the plastic hardening and damage softening functions. The macro crack behavior is characterized by a triaxiality dependent fracture criterion. The continuum damage model is implemented into LS-DYNA as user defined material model. Corresponding numerical simulations of unnotched and notched tension tests with high strain rates demonstrate the plastic and damage processes during the deformation leading to final fracture numerically predicted by an element erosion technique. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical approach for the simulation of ductile fracture with embedded discontinuities

In the present study, a computational approach for the numerical simulation of ductile fracture within the framework of the finite element method is proposed. In the developed macroscopic formulation, the inelastic behavior in the bulk of the material is described by the finite elasto‐plastic material model proposed in [4]. The failure process is modeled by introducing discontinuities when a special local fracture criterion is satisfied. The discontinuities are incorporated via special triangular finite elements with embedded interfaces following the line of [2]. Finally, the numerical procedure is evaluated for a twodimensional representative test problem. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

半无穷大裂纹端部粘聚力分析

准脆性材料裂纹端部断裂过程区粘聚力是导致非线性断裂特性的重要原因,根据准脆性材料的断裂特性,对存在粘聚力分布的半无穷大裂纹力学分析模型,由变形叠加原理得到以该粘聚应力分布为未知函数的积分方程,通过对积分方程的分析推证,得到了该分布函数解的数学结构和级数型表达式;提出了由实际裂纹张开位移,确定裂纹端部粘聚力分布函数的两种方法:其一由连续的裂纹张开位移通过积分变换求解未知函数级数展开项的系数,其二是由离散的裂纹张开位移数据通过最小二乘法确定该函数;推导出了相应方法求解未知量的代数方程,并且给出了适当的算例和讨论。     

Constitutive modeling of cyclic stress softening in filled elastomers

Quasi-static cyclic loading tests on filled rubber-like materials reveal a significant stress softening in the first cycle. The magnitude of this softening, widely known as the Mullins effect, reduces in the next cycles until it reaches a stabilized value referred to as hysteresis. In this contribution, we associate the hysteresis with the fracture of carbon bonds in the filler network. In order to calculate the cyclic energy dissipation, we modify the classic concept of network decomposition [1, 2] and add a new network considered to be responsible for the breakage and re-formation of carbon black aggregates during loading and unloading. The proposed model is in-line with a wide range of experimental observations. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A gradient model for finite strain elastoplasticity coupled with damage

This paper describes the formulation of an implicit gradient damage model for finite strain elastoplasticity problems including strain softening. The strain softening behavior is modeled through a variant of Lemaitre's damage evolution law. The resulting constitutive equations are intimately coupled with the finite element formulation, in contrast with standard local material models. A 3D finite element including enhanced strains is used with this material model and coupling peculiarities are fully described. The proposed formulation results in an element which possesses spatial position variables, nonlocal damage variables and also enhanced strain variables. Emphasis is put on the exact consistent linearization of the arising discretized equations.

A numerical set of examples comparing the results of local and the gradient formulations relative to the mesh size influence is presented and some examples comparing results from other authors are also presented, illustrating the capabilities of the present proposal.