Analysis of crack growth and crack-tip plasticity in ductile materials using cohesive zone models

Crack initiation and crack growth resistance in elastic plastic materials, dominated by crack-tip plasticity are analyzed with the crack modeled as a cohesive zone. Two different types (exponential and bilinear) of cohesive zone models (CZMs) have been used to represent the mechanical behavior of the cohesive zones. In this work, it is suggested that different forms of CZMs (e.g., exponential, bilinear) are the manifestations of different micromechanisms-based inelastic processes that participate in dissipating energy during the fracture process and each form is specific to each material system. It is postulated that the total energy release rate comprises the plastic dissipation rate in the bounding material and the separation energy rate within the fracture process zone, the latter is determined by CZMs. The total energy release rate then becomes a function of the material properties (e.g., yield strength, strain hardening exponent) and cohesive properties of the fracture process zone (e.g., cohesive strength and cohesive energy), and the form of cohesive zone model (CZM) that determines the rate of energy dissipation in the forward and wake regions of the crack. The effects of material parameters, cohesive zone parameters as well as the form/shape of CZMs in predicting the crack growth resistance and the size of plastic zone (SPZ) surrounding the crack tip are systematically examined. It is found that in addition to the cohesive strength and cohesive energy, the form (shape) of the traction–separation law of CZM plays a very critical role in determining the crack growth resistance (R-curve) of a given material. It is further observed that the shape of the CZM corresponds to inelastic processes active in the forward and wake regions of the crack, and has a profound influence on the R-curve and SPZ.     

数学网格和物理网格分离的有限单元法(I):基本理论

常规有限单元法在复杂边界问题的网格剖分、可移动边界和非连续变形问题的数值模拟等方面存在困难.本文将常规的有限单元分离为几何上相互独立的数学单元和物理单元,基于数学单元构造近似函数,引入位移模式关联法则以确定物理单元的位移模式,提出了在现有有限单元法框架内、基于数学网格和物理网格分离的强化有限单元法(FEM++).与常规有限单元法(SFEM)比较表明,强化有限单元法不仅很好地克服了常规有限单元法网格剖分上的困难,而且提供了一条更简便、更自然的分析移动边界问题和非连续变形问题的新途径.最后,通过数值算例验证了强化有限单元法的适用性和有效性.     

数学网格和物理网格分离的有限单元法(II):粘聚裂纹扩展问题中的应用

强化有限单元法将物理网格与数学网格分离开来,可以方便地描述非连续变形;粘聚区域模型是模拟断裂过程区作用最简单有效的方法,且可以避免裂纹尖端的应力奇异性.本文以平面问题为例,将强化有限单元法与粘聚区域模型相结合,利用富集数学节点描述任意粘聚裂纹扩展过程中的非连续变形问题,提出了裂纹扩展过程中数学节点富集和数学单元定义的方法.本文还导出了与平面4～8节点平面等参单元对应的8～16节点粘聚裂纹单元列武.最后,通过三点弯梁的裂纹扩展过程模拟验证了本文提出的粘聚裂纹扩展模拟方法的有效性.     

Cohesive zone modeling of interfacial delamination in PZT thin films

Previous experimental investigations [Shang, F., Kitamura, T., Hirakata, H., Kanno, I., Kotera, H., Terada, K., 2005. Experimental and theoretical investigations of delamination at free edge of interface between piezoelectric thin films on a substrate. International Journal of Solids and Structures 42 (5–6) 1729–1741] have demonstrated that multilayered Cr/PZT/PLT/Pt/Ti thin films deposited on single-crystal silicon substrates are delaminated along the interface between Cr and PZT layers in a brittle manner. This study starts with a model based on the cohesive zone concept and carries out numerical simulations to check the fracture behavior of this interfacial delamination. Three types of cohesive zone models (CZMs) are adopted, including the exponential, bilinear, and trapezoidal models. Characteristic CZM parameters are extracted through comparisons with experimental results. The simulation results show that (i) cohesive strength and work of separation are the dominating parameters in the CZMs; (ii) the bilinear CZM more suitably describes this brittle interfacial delamination; and (iii) in comparison with typical several mm-thick film/coating materials, the fracture energy of this weak Cr/PZT interface is quite low. Our study demonstrates the applicability of CZM in characterizing the interface fracture behavior of film materials with micrometer thicknesses.     

Some issues in the application of cohesive zone models for metal–ceramic interfaces

Cohesive zone models (CZMs) are being increasingly used to simulate discrete fracture processes in a number of homogeneous and inhomogeneous material systems. The models are typically expressed as a function of normal and tangential tractions in terms of separation distances. The forms of the functions and parameters vary from model to model. In this work, two different forms of CZMs (exponential and bilinear) are used to evaluate the response of interfaces in titanium matrix composites reinforced by silicon carbide (SCS-6) fibers. The computational results are then compared to thin slice push-out experimental data. It is observed that the bilinear CZM reproduces the macroscopic mechanical response and the failure process while the exponential form fails to do so. From the numerical simulations, the parameters that describe the bilinear CZM are determined. The sensitivity of the various cohesive zone parameters in predicting the overall interfacial mechanical response (as observed in the thin-slice push out test) is carefully examined. Many researchers have suggested that two independent parameters (the cohesive energy, and either of the cohesive strength or the separation displacement) are sufficient to model cohesive zones implying that the form (shape) of the traction–separation equations is unimportant. However, it is shown in this work that in addition to the two independent parameters, the form of the traction–separation equations for CZMs plays a very critical role in determining the macroscopic mechanical response of the composite system.     

Nonlinear interface shear fracture of end notched flexure specimens

The nonlinear analytical solutions of an end notched flexure adhesive joint or fracture test specimen with identical or dissimilar adherends are investigated. In the current study, a cohesive zone model (with arbitrary nonlinear cohesive laws) based analytical solution is obtained for the interface shear fracture of an end notched flexure (ENF) specimen with sufficiently long bond length. It is found that the scatter and inconsistency in calculating Mode II toughness may be significantly reduced by this model. The present work indicates that the Mode II toughness G_IIc under pure shear cracking condition is indeed very weakly dependent on the initial crack length. And this conclusion is well supported by the experimental results found in the literature. The parametric studies show that the interface shear strength is the most dominant parameter on the critical load. It is also interesting to note that with very short initial crack length and identical interface shear strength, higher Mode II toughness indeed cannot increase the critical load. Unlike the high insensitivity of critical load to the detailed shape of the cohesive law for Mode I peel fracture, the shape of the cohesive law becomes relatively important for the critical load of joints under pure Mode II fracture conditions, especially for joints with short initial crack length. The current study may help researchers deepen the understanding of interface shear fracture and clarify some previous concepts on this fracture mode.     

Separation work analysis of cohesive law and consistently coupled cohesive law

An appropriate coupled cohesive law for predicting the mixed mode failure is established by combining normal separation and tangential separation of surfaces in the cohesive zone model(CZM) and the cohesive element method.The Xu-Needleman exponential cohesive law with the fully shear failure mechanism is one of the most popular models.Based on the proposed consistently coupled rule/principle,the Xu-Needleman law with the fully shear failure mechanism is proved to be a non-consistently coupled cohesive la...     

Local Interface Shear Fracture of Bonded Steel Joints with Various Bondline Thicknesses

Interfacial fracture is a critical issue for extensive applications of adhesively bonded structures to a variety of modern industries. Extensive global experimental tests have been conducted to measure the global behavior of adhesively bonded joint, such as ultimate load capacity and toughness. Recently, several studies have also been employed to characterize the local interfacial traction–separation laws. However, very few tests have investigated the dependency of the local interfacial constitutive laws on the adhesive thickness, particularly, under Mode-II loading conditions. In this work, six typical adhesive thicknesses (from 0.1 mm to 1.0 mm) are prepared for the bonded joints with a configuration of end notched flexure (ENF) specimen to realize the Mode-II fracture loading (shear fracture). With a recently developed analytical model, the global energy release rates of the ENF specimens are experimentally measured. Meanwhile, with the image analysis technique, the local slips between the two adherends are obtained. Finally, based on the J-integral theory, the local interfacial constitutive laws at different bondline thicknesses are obtained. Several experimental findings are reported in this work. This work may provide valuable baseline experimental data for the input in cohesive zone model (CZM) based analytical and numerical simulations.     

Separation work analysis of cohesive law and consistently coupled cohesive law

An appropriate coupled cohesive law for predicting the mixed mode failure is established by combining normal separation and tangential separation of surfaces in the cohesive zone model (CZM) and the cohesive element method. The Xu-Needleman exponential cohesive law with the fully shear failure mechanism is one of the most popular models. Based on the proposed consistently coupled rule/principle, the Xu-Needleman law with the fully shear failure mechanism is proved to be a non-consistently coupled cohesive law by analyzing the surface separation work. It is shown that the Xu-Needleman law is only valid in the mixed mode fracture when the normal separation work equals the tangential separation work. Based on the consistently coupled principle and the modification of the Xu-Needleman law, a consistently coupled cohesive (CCC) law is given. It is shown that the proposed CCC law has already overcome the non-consistency defect of the Xu-Needleman law with great promise in mixed mode analyses.     

Plasticity aspects of interfacial fracture for polymer/glass specimens

Experimental results suggest that the interfacial fracture resistance is minimal for approximate near tip Mode I accompanied by positive and negative near tip Mode II. Finite-strain FE analysis is made for an elastic–plastic medium bonded to an ideally elastic medium with an interface crack. Small-scale plasticity conditions are invoked and examined in relation to the elastic–plastic stress distribution along the bond line. Plasticity engenders a tendency to turn near tip biaxiality towards pure Mode I regardless of the mixed-mode loading. High levels of hydrostatic stress are attained. For different mode mixities of the applied load, the dependence of the elastic–plastic normal bond stress on load level is examined. It is found that under positive Mode II loading, the normal bond stress σ_yy tends to saturate as the load level rises. This does not occur for Mode I and negative Mode II loading. In addition, deformation patterns inside the plastic zone are examined for mixed-mode situations. A displacement criterion based on the normal bond crack opening suggests a dependence of the critical load level on the extent of mixed mode. Under positive mode II fracture, traces of the ductile material are found at the top of the elastic substrate. Some of these conclusions appear to be consistent with the fracture patterns observed for LD-polyethylene/glass interfacial mixed-mode fracture.     

金属材料层裂破坏的内聚力模型

本文把内聚力单元嵌入到连续介质有限元单元之间,构建了一个层裂破坏的内聚力模型,以计及层裂破坏过程中能量耗散行为。采用该模型对平板撞击条件下的20号钢层裂实验进行了数值模拟研究,重点讨论了内聚力模型特征参数对计算结果的影响规律。研究结果表明采用指数型损伤演化行为的内聚力模型可以较好地描述弹塑性材料层裂破坏过程中的非线性能量耗散行为。利用一发实测自由面速度波剖面对计算结果进行对比校准,可确定内聚力模型特征参数。该特征参数可同样成功地预示不同撞击速度下的层裂实验,获得的模拟曲线与实验曲线之间符合程度很好,特别是自由面速度“回跳”后波形振荡周期和幅值与实验结果非常接近。这表明了内聚力模型在描述层裂过程中能量耗散行为方面具有较好适用性,并且不难由简单实验标定相关的特征参数。     

Constraint effects in adhesive joint fracture

Constraint effects in adhesive joint fracture are investigated by modelling the adherents as well as a finite thickness adhesive layer in which a single row of cohesive zone elements representing the fracture process is embedded. Both the adhesive and the adherents are elastic-plastic with strain hardening. The bond toughness Γ (work per unit area) is equal to Γ₀+Γ_p, where Γ₀ is the intrinsic work of fracture associated with the embedded cohesive zone response and Γ_p is the extra contribution to the bond toughness arising from plastic dissipation and stored elastic energy within the adhesive layer. The parameters of the model are identified from experiments on two different adhesives exhibiting very different fracture properties. Most of the tests were performed using the wedge-peel test method for a variety of adhesives, adherents and wedge thicknesses. The model captures the constraint effects resulting from the change in Γ_p: (i) the plastic dissipation increases with increasing bond line thickness in the fully plastic regime and then decreases to reach a constant value for very thick adhesive layers; (ii) the plastic dissipation in the fully plastic regime increases drastically as the thickness of the adherent decreases. Finally, this model is used to assess a simpler approach which consists of simulating the full adhesive layer as a single row of cohesive elements.     

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

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

Thermal de-cohesion model for time of craze failure in cohesive zone

The cohesive parameter corresponding to craze failure time is predicted for thermoplastics material. A craze failure separation criterion is proposed for a cohesive zone subjected to a melt layer formed and thickened by adiabatic deformation heat from a craze drawing. The numerical simulation of cohesive zone separation is based on non-linear thermal conduction and convection in the craze region and bulk region around the active layer, associated with a mechanical craze fibrils drawing in an uniaxial direction. The craze failure time is predicted with the assumption of the constant craze thickening rate and cohesive stress for a pipe-grade polyethylene. The numerically computed model reveals the inverse power law decay of the craze failure time, t_f, with increasing in craze thickening rate, v_c, (almost, t_f∝V_c⁻¹) for the thermoplastics. The full notch impact test experimental results are consistent with the analysis prediction. It is concluded that the craze failure time can be theoretically predicted using the numerical modeling.     

内聚力模型的形状对胶接结构断裂过程的影响

内聚力模型被广泛应用于粘接结构的断裂数值模拟过程中,为深入分析不同形状内聚力模型与胶黏剂性质和粘接结构断裂之间的关系,本文分别采用脆性和延展性两种类型胶黏剂,对其粘接的对接试件进行了单轴拉伸、剪切实验,以及其粘接的双臂梁试件进行了断裂实验.3种类型的内聚力模型(抛物线型、双线型和三线型)分别模拟了以上粘接结构的断裂过程,并与实验结果进行对比.结果发现:双线型的内聚力模型适用计算脆性胶黏剂的拉伸与剪切的断裂过程;指数型内聚力模型较适合计算延展性胶黏剂的拉伸和剪切的断裂过程,临界应力、断裂能和模型的形状参数是分析拉伸和剪切的重要参数;双臂梁试件的断裂过程模拟结果发现,断裂曲线与胶黏剂性质有关,内聚力模型形状参数也有影响.通过实验与计算结果分析,双线型内聚力模型更适合脆性胶黏剂粘接的双臂梁断裂计算,而三线型更适合计算延展性胶黏剂粘接的双臂梁断裂过程,此研究结果对胶黏剂的使用和粘接结构的断裂分析有很重要意义.     

An improved atomistic simulation based mixed-mode cohesive zone law considering non-planar crack growth

A novel and improved atomistic simulation based cohesive zone law characterizing interfacial debonding is developed which explicitly accounts for the non-planarity of the crack propagation. Group of atoms in the simulation constituting cohesive zones which are used to obtain local stress and crack opening displacement data are determined dynamically during the non-planar crack growth as they cannot be determined apriori. The methodology is used to study the debonding of Σ5 (2 1 0)/[0 0 1] symmetric tilt grain boundary interface in a Cu bicrystal under several mixed mode loading conditions. Simulations show that such bicrystalline specimen exhibits three types of energy dissipative mechanisms – shear coupled GB migration (SCM) away from the crack-tips, change in spacial orientation of GB structural units rendering highly disordered grain boundary near the crack tips and brittle intergranular fracture. Which combination of these three deformation mechanism will be active influencing the degree of non-planarity of the crack propagation at various stages of loading depends on the loading mode-mixity. As the ratio of shear component of the loading parallel to the GB plane and normal to the tilt axis with respect to the normal loading increases (thereby increasing the mode-mixity), overall strain-to-failure also increases and SCM tends to become the dominant deformation mechanism. Through this framework, analytical functional forms and parameters describing cohesive laws for both normal and shear traction as a function of the mode-mixity of the loading and crack opening displacement are predicted.     

Effects of interfacial properties on the ductility of polymer-supported metal films for flexible electronics

Polymer-supported metal films as interconnects for flexible, large area electronics may rupture when they are stretched, and the rupture strain is strongly dependent upon the film/substrate interfacial properties. This paper investigates the influence of interfacial properties on the ductility of polymer-supported metal films by modeling the microstructure of the metal film as well as the film/substrate interface using the method of finite elements and the cohesive zone model (CZM). The influence of various system parameters including substrate thickness, Young’s modulus of substrate material, film/substrate interfacial stiffness, strength and interfacial fracture energy on the ductility of polymer-supported metal films is systematically studied. Obtained results demonstrate that the ductility of polymer-supported metal films increases as the interfacial strength increases, but the increasing trend is affected distinctly by the interfacial stiffness.     

DETERMINATION OF INTERFACIAL MECHANICAL PARAMETERS FOR AN Al/EPOXY/Al2O3 SYSTEM BY USING PEEL TEST SIMULATIONS

Peel test measurements and simulations of the interfacial mechanical parameters for the Al/Epoxy/Al2O3 system are performed in the present investigation. A series of Al film thicknesses between 20 and 250 microns and three peel angles of 90, 135 and 180 degrees are considered. Two types of epoxy adhesives are adopted to obtain both strong and weak interface adhesions. A finite element model with cohesive zone elements is used to identify the interfacial parameters and simulate the peel test process. By simulating and recording normal stress near the crack tip, the separation strength is obtained, Furthermore, the cohesive energy is identified by comparing the simulated steady-state peel force and the experimental result. It is found from the research that both the cohesive energy and the separation strength can be taken as the intrinsic interfacial parameters which are dependent on the thickness of the adhesive layer and independent of the film thickness and peel angle.