分子动力学模拟含不同初始缺陷单晶镍的断裂行为与应力演化

断裂是一个跨尺度复杂的物理过程,对宏观尺度的断裂行为已有深入的研究和发展,然而对微观尺度的断裂行为及断裂过程中应力场的变化缺乏深入的理解。本文通过分子动力学模拟,研究了具有不同初始缺陷(尖锐裂纹、钝裂纹和孔洞)的单晶镍的断裂行为和应力分布特征。结果表明,不同的初始缺陷导致了不同的断裂机制、断裂强度和抗断裂性能。含初始孔洞的单晶镍样品有最高的断裂强度和最强的抗断裂性能,这与孔洞扩展过程中堆积层错的形成密切相关。其次是含初始钝裂纹的样品,在裂纹扩展过程中出现由[100]超位错发射引起的裂尖钝化;含尖锐裂纹的样品表现为脆性断裂,裂尖原子没有出现微结构的变化,其强度和抗断裂性能最低。此外,不同的初始缺陷也会导致断裂过程中应力分布的变化,对含有尖锐裂纹的脆性断裂试样,高应力(拉伸应力、平均应力和米塞斯应力)总是出现在扩展裂纹的裂尖。而对于含有钝裂纹或孔洞的韧性断裂试样,高应力不仅分布在裂尖,也分布在位错发射和堆积层错形成的区域,在裂纹/孔洞扩展之前,应力随着加载时间的增加而迅速增加,而一旦裂纹或孔洞开始扩展,应力增加非常缓慢或几乎不增加,但拉伸应力值始终大于平均应力和米塞斯应力值。这表明,在I型...     

几种金属材料断裂条件的试验研究

考察了几种金属材料在常规破坏试验与拉-扭双轴破坏试验中，随受力条件不同，材料不同形式断裂破坏变化规律和相应的断裂条件；利用几种韧性材料复合型断裂试验结果，分析了随应力三维度变化，材料中孔洞成核形状与聚合方式的变化规律，分析了不同断裂形式时启裂点、启裂方向变化规律及主要影响因素。研究表明，从材料断裂物理机制来看，裂纹体与无裂纹体断裂破坏实质是相同的，随材料塑性因素与应力三维度变化的影响。材料内孔洞成核形状由椭球形逐渐变化为细长形，材料断裂由正拉断转向剪断，裂纹体与无裂体材料的断裂形式、断裂危险点、危险点上断裂方向等宏观量也有着相同的变化规律；区分不同物理机制建立断裂条件，可能适合解决金属材料不同形式的断裂力学问题。     

含氦泡辐照老化材料层裂损伤计算方法分析

辐照条件下,一些材料内部产生大量的氦泡等微缺陷,氦泡的大小和数密度随着辐照年限的增长而增长。氦泡分布特征的变化不仅影响材料本身的物理、力学性质,而且直接影响材料层裂损伤演化后期材料破坏颗粒度的分布特征。延性材料的层裂损伤演化过程一般包括孔洞的成核、增长和汇合,但因已有孔洞对新成核孔洞存在抑制作用,当初始孔洞数密度达到一定临界值时,材料内部没有新的孔洞成核,因此,层裂损伤的计算可以不考虑新孔洞成核的影响。本文中基于损伤早期演化的特征,给出了这一临界值的计算方法,并进一步探讨了含氦泡辐照老化钚材料层裂损伤的计算方法。同时,在完善孔洞增长(void growth, VG)层裂损伤模型中参数的确定方法的基础上,借助含氦泡常规铝材料的层裂实验结果,对此问题进行了定性的分析：在氦泡尺寸变化不大的情况下,当氦泡浓度低于临界氦泡浓度时,需要考虑初始氦泡以及新增孔洞的综合影响;反之,可以采用简单的层裂损伤模型,不需要计算孔洞成核,但由于增长孔洞之间的相互影响,损伤模型的初始损伤参数需要重新确定。     

影响三维编织CMC断裂性能的孔洞问题

在气相沉积三维编织陶瓷基复合材料(CMC)的编织结构中，孔洞的存在在所难免。研究孔洞对断裂性能的影响必须面对孔洞带来的问题。过去研究孔洞问题的工作主要集中在等效模量和强度上，并取得了一些理论成果。然而，除了微裂纹、纤维和偏转增韧等的研究外，很少有关于孔洞对断裂性能的文献报道。因此这里探讨一下孔洞对断裂的影响问题。通常研究者认为孔洞能够通过两种途径增韧材料的断裂韧性：即裂尖的钝化和材料柔性增加。前者只能够定性地认识和分析裂尖形貌的增韧效果(本文不涉及此内容)；后者可以做一些定量的分析。基于此，首先研究了孔洞的形态与分布，并给出了孔洞模型。其次用两次等效的工程方法均匀化材料的有效弹性模量，并且和实验相比较，其结果令人满意。接着又分析了应力强度因子K及K控制区随孔洞体积分量的变化规律。最后引用能量释放率G用图的形式表述了孔洞的增韧效果。显然图中的结论表明，孔洞的增韧效果并不理想。     

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

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

三点弯曲试样动态应力强度因子计算研究

利用Hopkinson压杆对三点弯曲试样进行冲击加载,采集了垂直裂纹面距裂尖2mm和与裂纹面成60°距裂尖5mm处的应变信号。根据裂尖附近测试的应变信号计算试样的动态应力强度因子,并与有限元计算结果进行比较,结果表明由于裂尖有一段疲劳裂纹区,通过裂尖附近应变信号来计算动态应力强度因子时,如果裂尖位置确定不准及粘贴应变片位置不够准确对计算结果将带来很大影响。因此利用应变片法计算动态应力强度因子时,为了获得更准确的计算结果,在实验后应对试件裂纹面进行分析测量,重新确定裂尖位置,必要时需对应变片至裂尖距离进行修正后再计算动态应力强度因子值。     

混凝土的断裂参数及断裂过程区

由于混凝土断裂破坏时，裂尖具有微裂缝区及主缝的亚临界扩展这一特征，采用实验室小尺寸试件和线弹性断裂力学确定的断裂参数，往往有明显的尺寸效应，即ｋＩＣ随试件的尺寸的增大而增大，净截面名义应力随试件尺寸的增大而减小。     

韧性断裂的微观模型及其在约束断裂中的应用

详细分析了不同形状断裂试件及小范围屈服模型裂纹端部的损伤演化，提出了韧性断裂的宏观起裂相当于裂尖前方一特征位置处的损伤达到一临界值。利用此模型获得了与实验相一致的宏观断裂韧性及与约束无关的理论断裂韧性。     

高强铝合金的动态拉伸断裂行为实验研究

在分离式霍普金森拉杆、三点弯曲和平板撞击加载下对棒材铝合金（2024-T4、7075-T6）进行动态拉伸断裂实验研究。实验结果表明：1）一维应力动态加载下 7075-T6铝合金的初始屈服应力与断裂应变明显高于2024-T4铝合金,但三点弯曲和平板撞击层裂实验中发现2024-T4铝合金相比于7075-T6铝合金具有更好的抗裂纹扩展与层裂失效能力,这表明应力状态对两种铝合金拉伸断裂行为有明显的影响; 2）断口的光学与扫描电镜分析发现：2024-T4铝合金主要表现出脆性断裂行为,起因于孔洞或裂纹主要成核于晶内强化相形成穿晶断裂;而7075-T6铝合金则展现出韧性和脆性混合断裂特征,原因是部分孔洞或裂纹在晶界成核增长发生沿晶断裂,部分在晶内强化相周围形成孔洞从而造成穿晶断裂。, 在分离式霍普金森拉杆、三点弯曲和平板撞击加载下对棒材铝合金（2024-T4、7075-T6）进行动态拉伸断裂实验研究。实验结果表明：1）一维应力动态加载下 7075-T6铝合金的初始屈服应力与断裂应变明显高于2024-T4铝合金,但三点弯曲和平板撞击层裂实验中发现2024-T4铝合金相比于7075-T6铝合金具有更好的抗裂纹扩展与层裂失效能力,这表明应力状态对两种铝合金拉伸断裂行为有明显的影响; 2）断口的光学与扫描电镜分析发现：2024-T4铝合金主要表现出脆性断裂行为,起因于孔洞或裂纹主要成核于晶内强化相形成穿晶断裂;而7075-T6铝合金则展现出韧性和脆性混合断裂特征,原因是部分孔洞或裂纹在晶界成核增长发生沿晶断裂,部分在晶内强化相周围形成孔洞从而造成穿晶断裂。     

基于有限断裂法和比例边界有限元法的裂纹扩展模拟

基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。     

裂纹面摩擦接触引起的断裂韧性增长的研究

采用弹黏塑性的材料本构关系, 建立了压、剪混合型裂纹常速准静 态扩展的力学模型, 求得了裂纹面摩擦接触条件下裂纹尖端场的数值解, 并基于数 值结果讨论了扩展裂纹的摩擦效应. 计算和分析表明, 裂纹面的摩擦效应主要表现 在两个方面. 第一方面是摩擦会导致裂纹尖端区材料的断裂韧性增高, 并且裂纹面间的摩擦作用越强, 增韧效果越显著. 摩擦增韧的机制可以解释为裂纹 面间的摩擦作用导致裂纹尖端塑性区尺寸变大, 使裂纹尖端场的塑性变形能增加, 从而使得裂纹尖端区材料增韧. 摩擦生热并不是导致材料断裂韧性增长的根本机制. 第二方面是摩擦会导致``断裂延缓'. 利用裂纹面的摩擦来提高构件的承载能力和延长构件的服役寿命具有较大的工程实用价值.     

Dynamics of crack penetration vs. branching at a weak interface: An experimental study

In this paper, the dynamic crack-interface interactions and the related mechanics of crack penetration vs. branching at a weak interface are studied experimentally. The interface is oriented perpendicular to the incoming mode-I crack in an otherwise homogeneous bilayer. The focus of this investigation is on the effect of interface location and the associated crack-tip parameters within the bilayer on the mechanics of the ensuing fracture behavior based on the optical methodologies laid down in Ref. Sundaram and Tippur (2016). Time-resolved optical measurement of crack-tip deformations, velocity and stress intensity factor histories in different bilayer configurations is performed using Digital Gradient Sensing (DGS) technique in conjunction with high-speed photography. The results show that the crack path selection at the interface and subsequently the second layer are greatly affected by the location of the interface within the geometry. Using optically measured fracture parameters, the mechanics of crack penetration and branching are explained. Counter to the intuition, a dynamically growing mode-I approaching a weak interface at a lower velocity and stress intensity factor penetrates the interface whereas a higher velocity and stress intensity factor counterpart gets trapped by the interface producing branched daughter cracks until they kink out into the next layer. An interesting empirical observation based on measured crack-tip parameters for crack penetration and branching is also made.     

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.     

A general solution method for an anti-plane shear crack dynamically accelerating along a bimaterial interface

We develop a general solution method for a dynamically accelerating crack under anti-plane shear conditions along the interface between two different homogeneous isotropic elastic materials. The crack is initially at rest, and after loading is applied the crack-tip speed which may accelerate up to the shear wave speed of the more compliant material. The analysis includes an exact, closed-form expression for the stress intensity factor for an arbitrary time-dependent crack-face traction, as well as expressions for computing the crack-face displacements and the stress in front of the crack. We also present some numerical examples for fixed loads and for loads moving with the crack tip, using a stress intensity factor fracture criterion, in order to examine the predicted effect of material mismatch on interfacial fracture.     

BRITTLE-DAMAGE ANALYSIS OF CRACKS UNDER CONDITIONS OF PLANE-STRAIN TENSILE LOADING

An abrupt damage model, taking full account of finite geometrychanges, is used to study both the shapes of damage zones and the stress strain fieldsfor a plane-strain tensile crack under the small-scale yielding condition. Two typicalcrack-tip damage profiles are simulated by the element vanish technique. The fracturetoughness increment due to damage dissipation is evaluated in terms of the energyrelease rate.     

Analyzing the Laws of Stable Subcritical Growth of Cracks in Polymeric Materials on the Basis of Fracture Mesomechanics Models. Theory and Experiment

The laws of stable crack growth are analyzed using fracture mesomechanics models for polymeric materials under long-term subcritical tension. A review is given of experimental and theoretical studies of crack-tip process zones. The studies were conducted using physical and mechanical methods and fracture mesomechanics models, allowing for the structural and rheological features of process zones. Theoretical and experimental results on the behavior of process zones are generalized and the theory of stable crack growth in viscoelastic polymers, which assumes the autonomy of the process zone during crack development, is justified     

基于扩展有限元法的混凝土细观断裂破坏过程模拟

扩展有限元法（XFEM）是分析不连续力学问题（特别是断裂问题）的一种有效的数值方法。在常规的有限元位移模式中,基于单位分解的思想加入一个跳跃函数和渐进缝尖位移场来对不连续体附近的节点自由度进行局部加强,从而反映了位移的不连续性。介绍了扩展有限元的基本原理,给出了扩展有限元进行混凝土开裂及裂纹扩展的分析方法,最后采用扩展有限元法模拟了湿筛混凝土单轴拉伸作用下及WinklerL-型混凝土板的细观断裂破坏过程。分析了混凝土裂纹萌生、扩展的过程及破坏形态,数值结果与实验结果吻合良好。研究表明：扩展有限元法通过特定的位移模式,使裂纹两侧不连续位移场的表达独立于网格划分,能有效地模拟混凝土材料细观断裂破坏过程。     

Subcritical Growth of an Internal Circular Crack in an Aging Viscoelastic Laminated Composite

Delayed fracture of a laminated composite under tensile loads applied at infinity is studied. The composite consists of alternating elastic and aging viscoelastic layers and contains an internal penny-shaped mode I macrocrack located in parallel to the layers. A modified Leonov–Panasyuk–Dugdale crack model and the critical crack-tip opening criterion constitute a fracture model. The subcritical crack growth equations are derived using the Volterra principle and the method of operator continued fractions. The laws governing delayed fracture are studied for a specific composite material     

Evaluation of fracture behavior of iron aluminides

Comparative fracture tests of three Fe-28%Al iron aluminides showed that alloys with Zr and C addition (FA-187) or with B, Zr, and C addition (FA-189) are extrinsically more susceptible to environmental embrittlement than the base ternary alloy (FA-186) under constant tensile loading condition. This may be caused by the variations of grain boundary morphology (i.e. changes of grain size and grain boundary cohesive strength) caused by the alloy addition. The effect of grain boundary size and cohesive strength are further investigated with reference to the susceptibility of hydrogen embrittlement. Finite element simulation of initial intergranular fracture of two iron aluminides (FA-186 and FA-189) are made. The computational scheme involves coupling the stress and mass diffusion analyses to determine crack-tip stress state and the crack tip hydrogen diffusion. Maximum strain failure criteria was adopted to simulate intergranular fracture. The numerical modeling results correlated well with the experimental data. The result further confirmed that the grain boundary morphology is important as it appears to control the intrinsic and extrinsic fracture behavior of iron aluminides.     

The influence of a rate sensitive effect of material for a propagating crack-tip field

The elastic-viscoplastic model proposed in [1] is used in this paper to analyze the quasi-statically growing crack field and dynamic propagating crack-tip field of mode III. By analysis, we obtained more reasonable results than those of an elastic-plastic material. When the effect of rate sensitivity of a material is considered, it is found that only the quasi-statically growing crack-tip field is the special case of dynamic propagating crack-tip field when Mach numberM approaches zero.