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.     

A cohesive law for carbon nanotube/polymer interfaces based on the van der Waals force

We have established the cohesive law for interfaces between a carbon nanotube (CNT) and polymer that are not well bonded and are characterized by the van der Waals force. The tensile cohesive strength and cohesive energy are given in terms of the area density of carbon nanotube and volume density of polymer, as well as the parameters in the van der Waals force. For a CNT in an infinite polymer, the shear cohesive stress vanishes, and the tensile cohesive stress depends only on the opening displacement. For a CNT in a finite polymer matrix, the tensile cohesive stress remains the same, but the shear cohesive stress depends on both opening and sliding displacements, i.e., the tension/shear coupling. The simple, analytical expressions of the cohesive law are useful to study the interaction between CNT and polymer, such as in CNT-reinforced composites. The effect of polymer surface roughness on the cohesive law is also studied.     

Subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone

A recent experimental study has demonstrated the attainability of intersonic shear crack growth along weak planes in otherwise homogeneous, isotropic, linear elastic solids subjected to remote loading conditions (Rosakis et al., Science 284 (5418) (1999) 1337). The relevant experimental observations are summarized briefly here and the conditions governing the attainment of intersonic crack speeds are examined. Motivated by experimental observations, subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone is subsequently analyzed. A cohesive law is assumed, wherein the cohesive shear traction is either a constant or varies linearly with the local sliding rate. Complete decohesion is assumed to occur when the crack tip sliding displacement reaches a material-specific critical value. Closed form expressions are obtained for the near-tip fields. With a cohesive zone of finite size, it is found that the dynamic energy release rate is finite through out the intersonic regime. Crack tip stability issues are addressed and favorable speed regimes are identified. The influence of shear strength of the crack plane and of a rate parameter on crack propagation behavior is also investigated. The isochromatic fringe patterns predicted by the analytical solution are compared with the experimental observations of Rosakis et al. (1999) and comments are made on the validity of the proposed model.     

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.     

A thermomechanically coupled viscoelastic cohesive zone model at large deformation

A new viscoelastic cohesive zone model is formulated for large deformation conditions and within a fully coupled thermomechanical framework. The model is suitable for the simulation of a wide range of problems especially for polymeric materials. It can capture viscoelastic crack propagation as well as energy dissipation due to this process. Starting from the principles of thermodynamics, a 3D finite element formulation is derived for a fully coupled simultaneous solution of the thermal field and the deformation field. The viscoelastic model is constructed by extending an elastic exponential traction separation law using a simple rheology. The viscous part of the tractions is postulated to have the same characteristic length as the elastic part and that they are related by a single material parameter. A Newtonian dashpot is used to describe the evolution of the viscous separation. Furthermore, thermal effects are accounted for using temperature expressions in both the traction laws and the viscosity of the dashpot, and using a heat conduction law across the interface. The model is implemented within an implicit finite element code and the internal variable is calculated using an internal iteration. Different numerical examples are used to verify the model and a comparison with experimental data shows a satisfactory agreement.     

The cohesive law for the particle/matrix interfaces in high explosives

The debonding of particle/matrix interfaces has an important effect on the macroscopic behavior of composite materials. There are extensive analytical and numerical studies on interface debonding in composite materials based on cohesive zone models which assume a phenomenological relation between the normal (and shear) traction(s) and opening (and sliding) displacement(s) across the particle/matrix interface. However, there are little or no experiments to determine the cohesive law for particle/matrix interfaces in composite materials. In this paper, we develop a method to determine the cohesive law for particle/matrix interfaces in the high explosive PBX 9501. We use the digital image correlation technique to obtain the stress and displacement around a macroscopic crack tip in the modified compact tension experiment of PBX 9501. We use the extended Mori-Tanaka method (which accounts for the effect of interface debonding) and the equivalence of cohesive energy on the macroscale and microscale to link the macroscale compact tension experiment to the microscale cohesive law for particle/matrix interfaces. Such an approach enables us to quantitatively determine key parameters in the microscale cohesive law, namely the linear modulus, cohesive strength, and softening modulus of particle/matrix interfaces in the high explosive PBX 9501. The present study shows that Ferrante et al.'s [1982 Universal binding energy relations in metallic adhesion. In: J.M. Georges (Ed.), Microscopic Aspects of Adhesion and Lubrication, Elsevier, Amsterdam, pp. 19-30.] cohesive law, which is established primarily for bimetallic interfaces, is not suitable to the high explosive PBX 9501.     

MECHANLCAL FAILURE IN HAVERSIAN SYSTEMS ON THE BASIS OF MICROSTRUCTURE

The aim of this paper is to observe variations in the mechanical properties of Haversian systems including Haversian lamellae and canals,by means of the compressive test of five bone specimens subjected to different loads corresponding to 2KN,3KN,dKN,5KN and 6KN respectively. From a series of related microscopic pictures it is clearly seen that:(1) in the lower level of loading (2KN) ,the Haversian canals much the same as the one unsubjected to loads,and conditions of Harversian lamellae are also the same,but,the individual Harversian canal was slightly bent; (2) in the case of the higher levels of loading,it is first deformed in the weaked portion of the Haversian lamellae; (3) with increasing loads,the fracture shapes of the Haversian system are very complex,but the fractures always take place in the cement line between osteons.At the same time,deformed dimensions of the Harversian lamellae and canals were measured. The observed phenomena were qualitatively interpreted by the theory of the linear viscoelas-ticity.     

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.     

Numerical simulation of quasi-brittle fracture using damaging cohesive surfaces

The cohesive surface methodology is used in a numerical study of fracture of concrete. The traction vs. separation response is governed by an isotropic damage law in which damage evolves according to a prescribed one-dimensional linear or exponential softening law. Cohesive surfaces are immersed in the continuum to allow for a maximum freedom of crack path selection. The single edge notched four point shear beam and the double edge notched tensile bar are used to study: (i) the influence of the tangential cohesive response on the development of the fracture path and (ii) the mesh alignment sensitivity. It is shown that in the present formulation, the tangential cohesive response has a minor influence on both crack path and global characteristics. Mesh alignment does have a significant influence on the outcome of the numerical analysis.     

CRACK PROPAGATION IN POLYCRYSTALLINE ELASTIC-VISCOPLASTIC MATERIALS USING COHESIVE ZONE MODELS

Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior's material properties (e.g. strain rate sensitivity) characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.     

用于三维裂纹扩展分析的有限变形不可逆内聚元

为了能准确和高效的跟踪动态裂纹扩展,我们发展了三维有限变形的内聚元和一系列不可逆内聚力关系.该内聚元通过采用不可逆内聚力关系来控制裂纹两侧物质的逐渐分离和形成自由表面,这一点可类比于传统有限元对块体材料的离散化.为了展示该方法的预测能力及便于灵活使用的特点,我们模拟了Zehnder和Rosakis$^{[1]}$所做的重物落下动态断裂实验,值得注意的是该方法可以近似模拟出裂尖的轨迹.      

考虑粘结面厚度的局部粘结单元法

传统无厚度粘结单元法CFEM (Cohesive finite element method)在模拟脆性材料断裂方面具有很强的优势,但也存在很大问题.一是单元尺寸增大,收敛性变差;二是单元尺寸变小,模型刚度发生折减.为了克服这两个问题,发展了考虑厚度的局部粘结单元法,即在裂纹可能扩展区插入具有一定厚度的粘结面单元.粘结面单元采用拓展虚内键本构(Augmented virtual internal bond)描述.由于考虑了厚度,粘结面交叉处会形成多边形空缺.为了弥补这一空缺,将其看作多边形键元胞,采用离散虚内键模型(Discretized virtual internal bond)对其建模,保证了模型的几何完整性.模拟结果表明,本文方法有效,克服了传统CFEM方法的刚度折减问题,提高了计算稳定性和收敛性.     

Bifurcation phenomena in the rigid inclusion power law matrix composite sphere

Bifurcation of interface separation related to cavity nucleation is analyzed for a radially loaded composite sphere consisting of a rigid inclusion separated from a power law matrix by a uniform, non-linear cohesive zone. Equations for the spherically symmetric and non-symmetric problems are obtained from a hyperelastic finite strain theory by a limiting process that preserves non-linear matrix and interface response at infinitesimal strain. A complete solution to the symmetric problem is presented including bifurcation load, stresses, and evolution of elasto-plastic boundary and interface separation. An analysis of non-symmetric bifurcation, under symmetric conditions of geometry and loading, yields the bifurcation load and first non-symmetric mode shape associated with rigid inclusion displacement. An energy analysis is carried out for both symmetric and non-symmetric problems in order to assess stability of spherically symmetric states to spherically symmetric and non-symmetric “rigid body mode” perturbations.Results are provided for an interface force law that captures interface failure in normal mode and linear response in shear mode. For the symmetric problem, (i) there are threshold parameter values above which bifurcation will generally not occur, (ii) threshold values below which there do not exist equilibria in the post bifurcation regime, (iii) bifurcation occurs after attainment of the maximum interface strength. For the non-symmetric problem, (i) bifurcation always occurs, although it can be delayed by interfacial shear, (ii) for the smooth interface, non-symmetric bifurcation occurs after attainment of the maximum interface strength and always precedes symmetric bifurcation.     

A thermodynamic method for the construction of a cohesive law from a nonlocal damage model

Several published papers deal with the possibility of replacing a damage finite element model by a combination of cohesive zones and finite elements. The focus of the paper is to show under which conditions this change of model can be done in an energy-wise manner.The objective is to build a cohesive model based on a known damage model, without making any assumption on the shape of the cohesive law. The method is characterized, on the one hand, by the use of a well-defined thermodynamic framework for the cohesive model and, on the other hand, by the idea that the main quantity which must be maintained through the change of model is the energy dissipated by the structure. An analysis of the stability criteria enables us to determine the domains of validity of the different models. Thus, we show that it is consistent to derive the cohesive law from a given nonlocal damage model because the occurrence of a discontinuity can be viewed as an alternative way to limit localization. The method is illustrated on one-dimensional examples and a numerical resolution method for the problem with a cohesive zone is presented.     

Constitutive behaviour of mixed mode loaded adhesive layer

Mixed mode testing of adhesive layer is performed with the Mixed mode double Cantilever Beam specimen. During the experiments, the specimens are loaded by transversal and/or shear forces; seven different mode mixities are tested. The J-integral is used to evaluate the energy dissipation in the failure process zone. The constitutive behaviour of the adhesive layer is obtained by a so called inverse method and fitting an existing mixed mode cohesive model, which uses a coupled formulation to describe a mode dependent constitutive behaviour. The cohesive parameters are determined by optimizing the parameters of the cohesive model to the experimental data. A comparison is made with the results of two fitting procedures. It is concluded that the constitutive properties are coupled, i.e. the peel and shear stress depend on both the peel and shear deformations. Moreover, the experiments show that the critical deformation in the peel direction is virtually independent of the mode mixity.     

粘结单元在模拟FRP层合板低速冲击响应中的应用

纤维增强树脂基复合材料层合板(fibre reinforced plastic composites,FRP)在航空、航天、交通、造船等诸多工程中得到了日益广泛的应用,而其在冲击载荷下的响应和破坏特别是分层一直为学术界所关注。本文中对FRP层合板在冲击载荷下的响应和破坏进行数值模拟,并通过引入粘结层重点研究其分层破坏。首先,介绍一种基于改进的粘结区域方法的粘结层损伤模型;其次,详细介绍了有限元模型建模过程和建模细节;最后,对有限元模型进行验证,并分析分层损伤发生的原因。模拟结果表明,该模型不仅能准确预测FRP层合板在低速冲击载荷下的载荷-时间曲线和载荷-位移曲线,还能成功地预测其分层破坏。     

Multiscale cohesive failure modeling of heterogeneous adhesives

A novel multiscale cohesive approach that enables prediction of the macroscopic properties of heterogeneous thin layers is presented. The proposed multiscale model relies on the Hill's energy equivalence lemma, implemented in the computational homogenization scheme, to couple the micro- and macro-scales and allows to relate the homogenized cohesive law used to model the failure of the adhesive layer at the macro-scale to the complex damage evolution taking place at the micro-scale. A simple isotropic damage model is used to describe the failure processes at the micro-scale. We establish the upper and lower bounds on the multiscale model and solve several examples to demonstrate the ability of the method to extract physically based macroscopic properties.     

裂纹扩展过程中线性内聚力模型计算的半解析有限元法

提出了求解基于线性内聚力模型的平面裂纹扩展问题的半解析有限元法,利用弹性平面扇形域哈密顿体系的方程,通过分离变量法及共轭辛本征函数向量展开法,推导了一个环形和一个圆形奇异超级解析单元列式,组装这两个超级单元能准确地描述裂纹表面作用有双线性内聚力的平面裂纹尖端场。将该解析元与有限元相结合,构成半解析的有限元法,可求解任意几何形状和载荷的基于线性内聚力模型的平面裂纹扩展问题。典型算例的计算结果表明本文方法简单有效,具有令人满意的精度。     

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

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