Atomistic/continuum simulation of interfacial fracture Part II: Atomistic/dislocation/continuum simulation

Coupled atomistic/dislocation/continuum simulation of interfacial fracture is performed in this paper. The model consists of a nanoscopic core made by atomistic assembly and a surrounding elastic continuum with discrete dislocations. Atomistic dislocations nucleate from the crack tip and move to the continuum layer where they glide according to the dislocation dynamics curve. An atoms/continuum averlapping belt is devised to facilitate the transition between the two scales. The continuum constraint on the atomic assembly is imposed through the mechanics atmosphere along the overlapping belt. Transmissions of mechanics parameters such as displacements, stresses, masses and momenta across the belt are realized. The present model allows us to explore interfacial fracture processes under different mode mixity. The effect of atomistic zigzag interface on the fracture process is revealed: it hinders dislocation emission from the crack tip, especially under high mode mixity. The project supported by the National Natural Science Foundation of China     

A piezoelectric biaxial loading device for interfacial fracture experiments

In this paper we describe the development of a new biaxial loading device for investigating mixed-mode fracture at bimaterial interfaces. The new device makes use of piezoelectric actuators and specially arranged flexures to provide independent displacements normal and tangential to the interface. Capacitive probes and special washer load cells were used for measuring the displacements and reactive loads, respectively. A closed-loop circuit was formed with a personal computer to control the applied displacements to within 10 nm. Preliminary experiments with quartz/epoxy/aluminum sandwich specimens with cracks growing between the quartz and the epoxy found that the intrinsic toughness of this interface was 30% lower than the value for a glass/epoxy interface. Crack opening interferometry measurements having a resolution of 30 nm revealed the presence of a cohesive zone whose size was about 0.5 μm.     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

压电复合材料粘接界面断裂有限元模拟

根据数字化FRMM(Fix-Ratio Mix-Mode)断裂试验,得到了压电复合材料试件的断裂韧性和位移及应变场。本文在试验的基础上,通过非线性有限元软件ABAQUS及用户子程序UMAT进行了模拟分析,采用基于损伤力学的粘聚区模型(CZM)对压电复合材料界面的起裂和脱胶扩展进行了分析,并与VCCT方法进行了比较。计算得到的荷载位移曲线更接近于试验结果,但在裂纹扩展路径上的吻合需要对粘聚区法则进一步修正。通过进一步对CZM参数进行分析,表明界面粘结强度和界面刚度对计算结果的影响很大。研究结果表明,粘聚区模型可以很好地表征压电复合材料弱粘接界面脱胶断裂问题。     

高应变率下断裂韧性实验的数值模拟

采用有限元软件ANSYS/LS-DYNA程序对静态和冲击荷载作用下的含裂纹半圆弯曲(SCB)实验进行了数值模拟。根据静态实验的模拟结果,提出了适合复合型加载的Ⅰ型应力强度因子拟合公式,采用该公式计算应力强度因子的最大误差不超过10%。动态实验的模拟结果表明:对于纯Ⅰ型加载的SCB实验,动态应力强度因子随着试样半径、支座间距以及相对裂纹长度的变化呈现规律性变化;当试样半径小于60mm、相对支座间距为1.2、相对裂纹长度在0.1~0.4范围内时,惯性效应的影响较小,采用静态拟合公式计算裂尖的动态应力强度因子的误差约10%;对于复合型加载的SCB实验,当相对裂纹长度为0.2~0.4、裂纹倾角在10°~40°范围内时,采用静态拟合公式计算裂尖的动态应力强度因子的误差小于10%。     

断裂过程的有限元模拟

讨论了材料断裂过程的有限元模拟技术。基于自适应有限元的一般原理，并针对多相材料的裂纹扩展的特点，提出了一种简化的高精度和高效率有限元网格的动态重新划分策略。裂纹被假设沿着单元之间的路径连续扩展，利用节点力释放技术生成新的裂纹自由表面，发展了一种可随裂尖连续移动的网格动态加密和释放方法。这种方法已在各种裂纹问题中得以实现与应用。     

基于数值模拟的混凝土系列有效断裂判据

基于ANSYS平台,对某系列混凝土断裂破坏试验的加载过程进行有限元数值模拟,获得了荷载-加力 点位移曲线.利用获得的曲线,在线性渐进叠加假定的基础上,确定了试件在失稳断裂时的临界有效裂缝长度.将模拟峰值荷载和临界有效裂缝长度代入由LEFM推导出的公式计算了考虑断裂过程区影响的有效断裂韧度.最后,通过对系列试验有效断裂韧...     

The dynamical complexity of work-hardening： a large-scale molecular dynamics simulation

We analyze a large-scale molecular dynamics simulation of work hardening in a model system of a ductile solid. With tensile loading, we observe emission of thousands of dislocations from two sharp cracks. The dislocations interact in a complex way, revealing three fundamental mechanisms of work-hardening in this ductile material. These are (1) dislocation cutting processes, jog formation and generation of trails of point defects; (2) activation of secondary slip systems by Frank-Read and cross-slip mechanisms; and (3) formation of sessile dislocations such as Lomer-Cottrell locks. We report the discovery of a new class of point defects referred to as trail of partial point defects, which could play an important role in situations when partial dislocations dominate plasticity. Another important result of the present work is the rediscovery of the Fleischer-mechanism of cross-slip of partial dislocations that was theoretically proposed more than 50 years ago, and is now, for the first time, confirmed by atomistic simulation. On the typical time scale of molecular dynamics simulations, the dislocations self-organize into a complex sessile defect topology. Our analysis illustrates numerous mechanisms formerly only conjectured in textbooks and observed indirectly in experiments. It is the first time that such a rich set of fundamental phenomena have been revealed in a single computer simulation, and its dynamical evolution has been studied. The present study exemplifies the simulation and analysis of the complex nonlinear dynamics of a many-particle system during failure using ultra-large scale computing.     

功能梯度双材料弱/微间断界面的冲击断裂分析

提出强间断、弱间断、微间断和全连续界面的概念与分类,建立功能梯度弹性双材料弱间断 界面冲击断裂问题的力学模型,采用积分变换法推导问题的Cauchy奇异积分方程,并用配 点法求得数值解. 分析表明,弱/微间断性对于FGMs界面裂纹应力强度因子有着重要影响, 而且微间断性是优于弱间断性的一种界面力学性能连接关系. 以FGMs界面某一侧 的力学性能函数在界面处的Taylor展开式的低阶项作为界面另一侧的力学性能函数,便可 以使FGMs界面成为``微间断'界面. 界面的一阶微间断对应力强度因子的减小作用较为明 显,而高阶(二阶及以上)微间断对应力强度因子的影响较小. 减小界面的弱间断程度或使 FGMs界面具备``微间断性',都将利于提高功能梯度双材料界面抗冲击断裂能力,在一定 程度上达到界面增韧的目的.     

Application of split Hopkinson tension bar technique to the study of dynamic fracture properties of materials

A novel approach is proposed in determining dynamic fracture toughness(DFT) of high strength steel,using the split Hopkinson tension bar(SHTB) apparatus,combined with a hybrid experimental-numerical method.The center-cracked tension specimen is connected between the bars with a specially designed fixture device.The fracture initiation time is measured by the strain gage method,and dynamic stress intensity factors(DSIF) are obtained with the aid of 3D finite element analysis(FEA).In this approach,the dimensions of the specimen are not restricted by the connection strength or the stress-state equilibrium conditions,and hence plane strain state can be attained conveniently at the crack tip.Through comparison between the obtained results and those in open publication,it is concluded that the experimental data are valid,and the method proposed here is reliable.The validity of the obtained DFT is checked with the ASTM criteria,and fracture surfaces are examined at the end of paper.     

缝合复合材料层间开裂的试验和数值研究

缝合复合材料层合板中贯穿厚度方向的缝线,能有效增强层合板的抗分层能力。本文对一种碳纤/环氧缝合复合材料层板进行了短梁三点弯试验,测得了压头的载荷-位移曲线,并观察了层间裂纹的扩展,证实了缝线对层间裂纹的阻滞作用。建立了三维有限元模型模拟了上述试验,模型中相邻的铺层之间布置了一层初始无厚度的界面单元,界面单元的失效自然模拟层间开裂,而缝线简化为面积等效的梁单元,数值结果与试验观测吻合。     

岩石、混凝土类材料断裂破坏有限元数值模拟中的有效方法

岩石、混凝土类材料断裂破坏有限元数值模拟中的网格重划,依据单元畸变和裂缝介质间的单元干涉作为网格重划判据,采用几何体重构技术把几何实体分解成能在ANSYS上实现六面体网格划分的几个部分,利用体积判断法确定新结点在旧单元的单元编号,在场量传递上采用基于解析性质的等参有限元逆变换,把旧网格场量信息传递到新网格中。本文对ANSYS进行二次开发,实现了三维网格重划,网格重划采用单元畸变和界面干涉两个判据,在网格再划分前进行几何体重构,提取变形后的点线面信息重新生成实体,充分利用AN-SYS的函数和体积判断法找到新结点在旧网格中的位置,在新旧网格间的场量传递中采用基于解析逆等参单元法。在平台上实现了三维有限元网格重划技术,最后利用方料的单轴压缩断裂模拟计算检验了传递前后等效塑性应变分布用载荷信息的变化,证明了所开发系统的正确性。     

无钴合金钢的层裂断裂及数值模拟研究

利用一级轻气炮作为加载手段,对无钴合金钢的层裂特性进行了研究,获得了Hugoniot弹性极限、层裂强度、层裂片厚度以及塑性应变率等动力学参数。对回收的样品进行了断口分析和金相分析,从宏微观角度分析了无钴合金在不同应变率下的断裂机制。利用LS-DYNA有限元程序对层裂现象进行了数值模拟,计算结果与实验结果基本吻合。     

加载速率对船用钢断裂韧性的影响

在加载速率为100～106MPam1/2/s的范围内,分别采用准静态、示波冲击、Hopkinson杆型试验装置对某船用钢进行了断裂韧性测试。试验结果表明,此钢的断裂韧性对加载速率敏感,即随着加载速率的升高,断裂韧性下降;应用位错动力学对实验现象作出了解释。     

Brittle fracture: Variation of fracture toughness with constraint and crack curving under mode I conditions

The effect of constraint on brittle fracture of solids under predominantly elastic deformation and mode I loading conditions is studied. Using different cracked specimen geometry, the variation of constraint is achieved in this work. Fracture tests of polymethyl methacrylate were performed using single edge notch, compact tension and double cantilever beam specimens to cover a bread range of constraint. The test data demonstrate that the apparent fracture toughness of the material varies with the specimen geometry or the constraint level. Theory is developed using the critical stress (strain) as the fracture criterion to show that this variation can be interpreted using the critical stress intensity factorK _Cand a second parameterT orA ₃,whereT andA ₃are the amplitudes of the second and the third term in the Williams series solution, respectively. The implication of this constraint effect to the ASTM fracture toughness value, crack tip opening displacement fracture criterion and energy release rateG _Cis discussed. Using the same critical stress (strain) as the fracture criterion, the theory further predicts crack curving or instability under mode I loading conditions. Experimental data are presented and compared with the theory.     

纤维复合材料损伤过程的数值模拟

利用界面断裂力学和有限元法数值模拟纤维增强复合材料的细观损伤过程，研究各种主要破坏模式之间的相互转变和影响，指出以断裂能和混合度表示的界面性能是控制复合材料损伤过程的主要细观参数。分析了界面韧度对破坏性能的影响，探讨了基于破坏模式控制的复合材料韧度设计的新途径。     

Three-dimensional interfacial fracture analysis of a one-dimensional hexagonal quasicrystal coating

In this paper, the three-dimensional (3D) interfacial fracture is analyzed in a one-dimensional (1D) hexagonal quasicrystal (QC) coating structure under mechanical loading. A planar interface crack with arbitrary shape is studied by a displacement discontinuity method. Fundamental solutions of interfacial concentrated displacement discontinuities are obtained by the Hankel transform technique, and the corresponding boundary integral-differential equations are constructed with the superposition principle. Green’s functions of constant interfacial displacement discontinuities within a rectangular element are derived, and a boundary element method is proposed for numerical simulation. The singularity of stresses near the crack front is investigated, and the stress intensity factors (SIFs) as well as energy release rates (ERRs) are determined. Finally, relevant influencing factors on the fracture behavior are discussed.     

基于离散元的杆件断裂全过程数值模拟

杆件的断裂会涉及到大变形、非线性以及不连续等问题,通常的数值计算方法模拟这种复杂力学行为具有局限性。本文基于颗粒离散元法DEM,将接触粘结处的分布式弹簧用梁纤维进行等效,提出了一种适于结构弹塑性分析的DEM纤维梁模型,然后在此基础上构建了构件断裂模拟算法以及纤维破环准则。将该模型应用于悬臂梁结构,模拟了悬臂梁从弹性到弹塑性阶段,再到断裂破坏的全过程,数值模拟得到的结构响应和截面开裂破坏形态均较合理。最后将该方法应用于单层网壳倒塌破坏模拟,并与网壳振动台倒塌试验进行对比,结果表明,数值模拟得到的杆件断裂过程及结构倒塌模式与试验现象一致,验证了该模型的正确性和适用性。     

Enhanced fracture energy and effects of temperature at spalling in ceramics: continuum damage modelling

Summary Enhanced fracture energy losses at spalling and the temperature dependence of the spalling strength of alumina ceramic bars are analysed on the basis of the experimental tests conducted both in room temperature and within the temperature range up to 1500°C at strain rates of some 500 s⁻¹. The experimental method and the measurements are first shortly outlined. The mechanical response of ceramic bars is modelled then as a heterogeneous distribution of brittle-elastic mesoelements undergoing continuum damage at the known strain history, corresponding to that registered in the experiments. The mesoelements are characterised by the values of initial damage randomly fluctuating within a given band-width superposed on a deterministic distribution, which corresponds to the fabrication conditions of the ceramic bars. The model has been tested in the evaluation of room-temperature experiments, its parameters: the average value of the initial damage, Young's modulus of the undamaged material and the energy absorption capacity in continuum damage are taken from the calibration fitting the experimental data. The registered energy losses at spalling, which exceed the static values of fracture energy by almost an order of magnitude, can be explained thus by the enhancement of the dissipation due to bulk damage, which is computed on the basis of the above parameters. The temperature change of the Young's modulus of the matrix material is taken as corresponding to the measured change of the uniaxial wave velocity in the bar, and corrected by the temperature change of the mass density. The analysis of the model shows that the drop in the spalling strength of the specimens with the increase of the temperature is phenomenologically related to the falling energy absorption capacity within the continuum damage mechanism. An explanation of this phenomenon is attempted, based on the grain-size-related mechanisms of the microfracture from pre-existing intergranular flaws distributed over the bulk of ceramics. Received 7 May 1999; accepted for publication 14 June 1999     

3-D numerical simulation of fracture processes in heterogeneous brittle materials

By using the lattice model combined with finite element methods and statistical techniques, a numerical approach is developed to establish mechanical models of three-dimensional heterogeneous brittle materials. A special numerical code is introduced, in which a lattice model and statistical approaches are used to simulate the initial heterogeneity of material properties. The size of displacement-load step is adaptively determined so that only few elements would fail in each load step. When the tensile principal strain in an element exceeds the ultimate strain of this element, the element is considered broken and its Young's modulus is set to be very low. Some important behaviors of heterogeneous brittle materials are indicated using this code. Load-displacement curves and figures of three-dimensional fracture patterns are also numerically obtained, which are similar to those observed in laboratory tests.