Spallation analysis with a closed trans-scale formulation of damage evolution

A closed, trans-scale formulation of damage evolution based on the statistical microdamage mechanics is summarized in this paper. The dynamic function of damage bridges the mesoscopic and macroscopic evolution of damage. The spallation in an aluminium plate is studied with this formulation. It is found that the damage evolution is governed by several dimensionless parameters, i.e., imposed Deborah numbersDe ^* andDe, Mach numberM and damage numberS. In particular, the most critical mode of the macroscopic damage evolution, i.e., the damage localization, is determined by Deborah numberDe ^*. Deborah numberDe ^* reflects the coupling and competition between the macroscopic loading and the microdamage growth. Therefore, our results reveal the multi-scale nature of spallation. In fact, the damage localization results from the nonlinearity of the microdamage growth. In addition, the dependence of the damage rate on imposed Deborah numbersDe ^* andDe, Mach numberM and damage numberS is discussed. The project supported by the National Natural Science Foundation of China (10172084, 10232040, 10232050, 10372012, 10302029) and the Special Funds for Major State Research Project (G200077305)     

固体破坏理论的若干问题

本文列举和初步探讨了位于当代固体破坏理论的前沿上的十五个问题。它们包括：固体破坏理论的封闭、宏微观非线性力学、尺度效应、多层次计算、细微观实验力学技术、断裂过程区的描述、异质体的强度与韧性、材料的多层次抗破坏设计、高速破坏过程、非均匀介质的动态断裂与前兆、疲劳的起源、破坏元技术、内损伤破坏与微结构形貌演化失稳的识别、材料的微结构演化、电致失效力学。     

斜波压缩下锡的相变和层裂行为

获取不同热力学路径下锡的动态响应实验数据,是深入研究其相变和损伤物理过程的基础.利用小型磁驱装置CQ-4完成了金属锡的斜波加载实验,获取了锡含有相变和层裂损伤物理信息的实验数据.实验结果显示,在加载段锡依次经历了弹塑性转变和β-γ相变两种物理过程,屈服强度约0.194 GPa,相变压力随着锡厚度的增加从7.54 GPa减小到7.14 GPa.在卸载段出现了明显的层裂损伤,层裂强度约1.1 GPa,与相同加载压力下冲击实验结果有巨大差异,层裂片厚度约0.38 mm.结合由锡的多相Helmholtz自由能计算的多相状态方程、Hayes相变动力学方程和损伤度理论,对斜波压缩实验过程进行一维流体动力学数值模拟,计算结果可以很好描述锡的弹塑性转变、相变和层裂三个物理过程.     

改进SPH方法在陶瓷材料层裂数值模拟中的应用

为了探讨铝飞片撞击陶瓷材料时的层裂现象,采用改进SPH方法模拟应力波在陶瓷材料中的传播。结果表明,当离散粒子分布不均匀时,数值模拟计算的自由面速度时程曲线与实测曲线吻合良好。对比CSPM方法,改进SPH方法的精度更高。提出适用于数值模拟的陶瓷材料损伤演化方程,对脉冲载荷下陶瓷/钢层合板层裂的破坏过程进行数值模拟,结果表明,由于陶瓷的波阻抗高于钢的,且抗压强度远高于抗拉强度,因此拉应力引起的层裂破坏是主要的。即使在材料内部传播的只是弹性压缩波,当弹性波到达材料界面时,由界面反射引起的卸载波也能导致陶瓷发生层裂破坏。     

From Gradient Damage Laws to Griffith’s Theory of Crack Propagation

This paper is devoted to the comparison of the evolution of damage governed by a gradient damage model with the evolution of a crack predicted by Griffith’s law. The analysis is made in a two-dimensional setting, assuming that damage is concentrated inside thin bands whose width is proportional to the internal length of the material. Taking advantage of the variational formulation based on the three principles of irreversibility, stability and energy balance, one introduces a generalized Rice path integral which contains terms involving the gradient of damage. Assuming that the internal length of the material is small by comparison with the dimension of the body, a separation of scales is achieved. Owing to the energy balance and the stability condition, one first proves some properties of this path integral with respect to the path. Then, one shows that the evolution of the damage zone is governed by Griffith’s law, the dissipated surface energy being given by the energy dissipated in the damage process zone.     

延性金属层裂自由面速度曲线特征多尺度模拟研究

以延性金属钽为研究对象,对钽在平板撞击下的层裂行为进行了多尺度下的数值模拟研究,从微观视角对自由面速度曲线上的典型特征进行了新的解读。在宏观尺度,对比分析了光滑粒子流体动力学法（smootfied particle hydrodynamics, SPH）与Lagrange网格法以及几种本构模型的模拟结果及其适用性。通过与实验数据的对比表明,Steinberg-Cochran-Guinan本构模型在层裂模拟中与实验数据吻合较好,通过改变加载条件获得了不同应变率下的自由面速度曲线,分析了不同应变率下的自由面速度曲线中的典型特征。在微观尺度,采用分子动力学方法获得层裂区域内损伤演化情况,揭示了宏观尺度自由面速度曲线典型特征所蕴含的物理内涵。分析表明,层裂表现为材料内部微孔洞形核、长大和聚集的损伤演化过程,自由面速度曲线上的典型特征与层裂区域的损伤演化过程存在密切关联。Pullback信号是层裂区域内微孔洞形核的宏观表征;自由面速度曲线的下降幅值在一定程度上反映了微孔洞的形核条件,由此计算得到的层裂强度实际上是微孔洞的形核强度。此外,Pullback信号后的速度回跳速率反映了微损伤演化的速率。     

45#钢的损伤演化方程和层裂准则研究

通过唯象分析和细观物理统计相结合的方法建立了一种韧性材料的损伤演化方程。在试验结果和内变量理论的基础上得到了45#钢的含损伤热—粘塑性本构关系。用有限差分方法计算了45#钢平板撞击所产生的应力波传播规律、损伤演化规律及层裂进程。通过自由面速度历史的数值模拟,并基于计算结果与试验结果间的最佳一致性,得到了损伤演化方程中的材料参数和极限损伤数值,并以此为依据建立了材料的应力率层裂准则。     

Multi-scale analysis in elasto-viscoplasticity coupled with damage

The purpose of this study is to present a micromechanical approach, based on the transformation field analysis (TFA), proposed by Dvorak, which has been generalized at Onera in order to analyze the nonlinear behavior of heterogeneous materials in elasto-viscoplasticity coupled with damage. In such analysis, the macroscopic constitutive equations are not purely phenomenological but are built up from multi-scale approaches starting from the knowledge of the properties of the constituents at the microscopic or mesoscopic scales. The model can take into account some local characteristics that can evolve during the thermo-mechanical applied loads or the manufacturing process, like the grain size for metallic alloys or the fiber volume fraction for composites.The determination of some specific tensors which are present in this formulation is closely linked to the microstructure morphology of heterogeneous materials constituting the macroscopic structure. For example, an Eshelby’s based approach is more appropriate to characterize polycrystalline materials with a random microstructure, while the homogenization of periodic media technique can be used for composite materials with a sufficiently regular microstructure. The proposed methodologies allowing to perform this nonlinear analysis across the scales are illustrated with examples based on the behavior of structures reinforced with a long fiber unidirectional metal matrix composite.     

固体的统计细观力学——-连接多个耦合的时空尺度

从固体力学所面临的新的挑战------多物理、多尺度耦合及其现状的描述开始, 以层裂 过程为例, 说明了这些多尺度非平衡问题的基本困难在于, 在固体中不同尺度上有不同的微 结构层次及不同的演化物理和速率. 接下来, 概述了一些针对这一困难的独特的思路及 其成果. 第3部分强调了一些统计平均方法的范式, 以及处理包含多个时间和空间尺度的问 题的新思路, 特别是非平衡损伤演化导致宏观失效的问题. 在第4部分, 简要评述了一些连 接多个空间和时间尺度的细观力学框架, 如位错理论, 物理细观力学, Weibull理论, 随机 理论等, 并且阐述了其中蕴含的跨尺度耦合的机理. 然后, 在第5部分, 回到了描述损 伤演化过程的框架, 也就是统计细观损伤力学以及它的跨尺度封闭近似. 基于这些跨尺度框 架, 在第6部分, 对控制跨尺度耦合的可能机理进行了评述和比较. 由于对失效时灾变 的洞察与跨尺度强耦合紧密相关, 一些非平衡和强相互作用的新概念在第7部分进行了讨 论. 最后, 以一个简短的总结和一些建议结束.     

Synchronous multi-scale observations on rock damage and rupture

This paper reports a multi-scale study on damage evolution process and rupture of gabbro under uniaxial compression with several experimental techniques, including MTS810 testing machine, white digital speckle correlation method, and acoustic emission technique. In particular, the synchronization of the three experimental systems is realized for the study of relationship of deformation and damage at multiple scales. It is found that there are significant correlation between damage evolution at small and large length scales, and rupture at sample scale, especially it displays critical sensitivity at multiple scales and trans-scale fluctuations.     

脆性岩石破裂过程损伤与渗流耦合数值模型研究

大量的实验结果表明，脆性岩石的渗透性不是一个常量，而是应力和应力诱发损伤破裂的函数．建立了一个描述非均匀岩石渗流-应力-损伤耦合数学模型(FSD Model)，开发出岩石破裂过程渗流-应力-损伤耦合分析计算系统(F-RFPA^2D)．在该系统中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流-损伤耦合迭代来实现．该系统能够对岩石试件在孔隙水压力和双轴荷载作用下裂纹的萌生、扩展过程中渗透率演化规律及其渗流-应力耦合机制进行模拟分析．最后，给出两个算例：算例1模拟载荷作用下岩石应力应变-渗透率全过程．模拟结果表明，非均匀性对岩石的应力峰值强度、峰值前后其渗透性演化规律及其破裂机制影响十分明显，模拟结果和实验结果较为一致；算例2模拟孔隙水压力作用下岩石拉伸断裂过程，通过和物理实验对比验证，验证了模型计算结果的可靠性。     

金属材料层裂再压实的模拟研究

激波在自由面卸载后金属内部经常出现层裂现象。若金属内层裂区再次受到冲击加载,则处于拉伸稀疏状态下的金属会逐渐被再次压实为密实介质,直至层裂区消失、再压实过程完成。由于金属层裂区初始拉伸状态的复杂性及再压实后物质状态的不确定性,复杂加载情况下宏观模拟该问题的可靠性验证存在困难。目前,在实验诊断难以准确给出金属层裂区进入再压实过程的初始状态及再压实状态的情况下,具有层裂区内部细节描述能力的直接数值模拟成为了验证宏观模拟可靠性的一种有效手段。首先,在直接数值模拟建模中将金属层裂区初始拉伸状态建模为仅含层裂片、仅含孔洞、同时含有孔洞与层裂片3类情况。然后,通过不同孔隙度、再压实速率、层裂片数及孔洞数下的直接数值模拟,统计得到了对应工况下金属层裂区的再压实状态。最后,在保证直接模拟与宏观模拟具有良好可比性的情况下,对层裂再压实过程进行了宏观建模及模拟分析。分析认为:在宏观网格断裂后处理算法使用全应力置零和温度不变的情况下,宏观模拟能够较好地模拟稀疏区内含层裂片情况下的金属层裂再压实过程及再压实状态;若金属层裂区内部以仅含孔洞的初始状态进入再压实过程,则无论孔洞塌缩是否形成界面喷射,宏观模拟均无法较好模拟该层裂再压实过程及再压实状态。     

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

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

CHARACTERISTIC DIMENSIONLESS NUMBERS IN MULTl-SCALE AND RATE-DEPENDENT PROCESSES

Multi-scale modeling of materials properties and chemical processes has drawn great attention from science and engineering. For these multi-scale and rate-dependent processes, how to characterize their trans-scale for-mulation is a key point. Three questions should be addressed:How do multi-sizes affect the problems?How are length scales coupled with time scales?How to identify emergence of new structure in process and its effect?For this sake, the macroscopic equations of mechanics and the kinetic equations of the microstructural transforma-tions should form a unified set that be solved simultaneously.As a case study of coupling length and time scales, the trans-scale formulation of wave-induced damage evolution due to mesoscopic nucleation and growth is discussed. In this problem, the trans-scaling could be reduced to two inde-pendent dimensionless numbers: the imposed Deborah number De=(ac)/(LV) and the intrinsic Deborah num-ber D = (nNc5)/V* ,where a. L, c, V and nN are wave speed, sample size, micr     

An Experimental Technique for Spalling of Concrete

The spalling strength of concrete is measured by examining the strain wave profiles in a polymer buffer bar behind the slender concrete bar specimen placed between a large diameter (Φ100 mm) Hopkinson bar and the buffer bar. The experimental results indicate that the spalling strength is related to not only the compressive strength of concrete but also the impact velocities (the loading rates). The rate effect of spalling strength mainly results from the different cracking paths in concrete under different impact velocities. However when the input compressive stress to specimen exceeds the threshold required to trigger the compressive damage, the spalling strength decreases due to the evolution and cumulation of compressive damage in concretes. The repeated impact loading experiments indicate that damage plays an important role in the spallation process of concrete. The high speed video of the spalling fracture process shows that multiple spalling fractures may occur in the scab and damage accumulation resulting from stress wave propagation in scab is the main reason for the producing of multiple spallations.     

Finite polycrystalline elastoplasticity and damage: multiscale kinematics

A physically realistic macroscopic decomposition of the deformation gradient for metallic polycrystals should explicitly account for all relevant sub-macroscopic kinematic processes, including lattice deformation, plastic flow, and evolution of damage, that significantly contribute to the homogenized deformation at the macroscale. The present work suggests such a decomposition, based on principles of volume averaging and focusing upon elastoplasticity and a variety of damage modes including intergranular fracture, void growth and coalescence, and shear discontinuities. This decomposition, of hybrid additive–multiplicative form, captures precisely the kinematics of arbitrarily anisotropic damage and also offers insight into mesoscopic distributions of residual elastic lattice strain attributed to heterogeneity of local deformation occurring at both intergranular and intragranular scales.     

Multiscale Dynamics of Heterogeneous Media in?the?Peridynamic Formulation

A methodology is presented for investigating the dynamics of heterogeneous media using the nonlocal continuum model given by the peridynamic formulation. The approach presented here provides the ability to model the macroscopic dynamics while at the same time resolving the dynamics at the length scales of the microstructure. Central to the methodology is a novel two-scale evolution equation. The rescaled solution of this equation is shown to provide a strong approximation to the actual deformation inside the peridynamic material. The two scale evolution can be split into a microscopic component tracking the dynamics at the length scale of the heterogeneities and a macroscopic component tracking the volume averaged (homogenized) dynamics. The interplay between the microscopic and macroscopic dynamics is given by a coupled system of evolution equations. The equations show that the forces generated by the homogenized deformation inside the medium are related to the homogenized deformation through a history dependent constitutive relation.     

滑移爆轰驱动钢管的层裂

将VG损伤模型推广到二维情况,考虑了最大主应力方向对损伤演化的影响,并使用显式断裂算法对20钢管在GI-920炸药滑移爆轰驱动下层裂的问题进行二维数值模拟。分析了一维内爆和滑移爆轰两种加载方式下作用于钢管外表面的压力及钢管内部受力状态的区别,考察了滑移爆轰加载方式下钢管外表面的受力随炸药厚度变化的规律,进而研究了钢管内损伤的分布和演化,以及裂纹的产生和扩展现象。计算得到的层裂片初始厚度随炸药厚度的变化规律与实验结果符合较好。     

Analysis of the spallation mechanism suppression in plasma-sprayed TBCs through the use of heterogeneous bond coat architectures

This paper critically examines the use of heterogeneous bond coats to increase the durability of plasma-sprayed thermal barrier coatings under spatially-uniform cyclic thermal loading. A major failure mechanism in these types of coatings involves spallation of the top coat caused by the top/bond coat thermal expansion mismatch concomitant with deposition-induced top/bond coat interfacial roughness, oxide film growth and creep-induced normal stress reversal at the rough interface’s peaks. The reduction of the top/bond coat thermal expansion mismatch aimed at increasing coating durability can be achieved by embedding alumina particles in the bond coat. Herein, we analyze the evolution of local stress and inelastic strain fields in the vicinity of the rough top/bond coat interface during thermal cycling, and how these fields are influenced by the presence of spatially uniform and non-uniform (graded) distributions of alumina particles in the metallic bond coat. The analysis is conducted using the higher-order theory for functionally graded materials which accounts for the high-temperature creep/relaxation effects within the individual TBC constituents. In the presence of two-phase bond coat microstructures, both the actual and homogenized bond coat properties are employed in the analysis in order to highlight the limitations of the prevalent homogenization-based approach applied to graded materials. The results reveal that the use of heterogeneous, two-phase bond coats, with spatially uniform or graded microstructures, while slightly suppressing the normal stress component evolution in the interfacial peak region, increases the magnitude of the shear stress component as well as the inelastic strain evolution in this region, thereby potentially promoting delamination initiation. The analysis based on homogenized bond coat microstructure produces misleading results relative to how the bond coat heterogeneity affects the magnitude of the normal and shear stress, and inelastic strain, components.