2.5维自愈合C/SiC复合材料的压缩力学行为

试验研究了2.5维自愈合C/SiC复合材料的压缩力学行为,根据材料的细观结构特点,建立了压缩载荷下的损伤力学模型,得到了经纬向压缩的非线性应力应变关系,预测结果与试验值吻合较好。结果表明,经向和纬向的力学行为不同,纬向的压应力逐渐增大时,切线模量逐渐增大,压缩强度为270.05MPa,而经向压应力逐渐增大时,层间损伤逐渐发生,经纱承受的弯矩越来越大,切线模量逐渐降低,破坏强度为128.66MPa。     

In-plane fracture of laminated fiber reinforced composites with varying fracture resistance: Experimental observations and numerical crack propagation simulations

A series of experimental results on the in-plane fracture of a fiber reinforced laminated composite panel is analyzed using the variational multiscale cohesive method (VMCM). The VMCM results demonstrate the influence of specimen geometry and load distribution on the propagation of large scale bridging cracks in the fiber reinforced panel. Experimentally observed variation in fracture resistance is substantiated numerically by comparing the experimental and VMCM load–displacement responses of geometrically scaled single edge-notch three point bend (SETB) specimens. The results elucidate the size dependence of the traction-separation relationship for this class of materials even in moderately large specimens, contrary to the conventional understanding of it being a material property. The existence of a “free bridging zone” (different from the conventional “full bridging zone”) is recognized, and its influence on the evolving fracture resistance is discussed. The numerical simulations and ensuing bridging zone evolution analysis demonstrates the versatility of VMCM in objectively simulating progressive crack propagation, compared against conventional numerical schemes like traditional cohesive zone modeling, which require a priori knowledge of the crack path.     

A simple cohesive zone model that generates a mode-mixity dependent toughness

A simple, mode-mixity dependent toughness cohesive zone model (MDG_c CZM) is described. This phenomenological cohesive zone model has two elements. Mode I energy dissipation is defined by a traction–separation relationship that depends only on normal separation. Mode II (III) dissipation is generated by shear yielding and slip in the cohesive surface elements that lie in front of the region where mode I separation (softening) occurs. The nature of predictions made by analyses that use the MDG_c CZM is illustrated by considering the classic problem of an elastic layer loaded by rigid grips. This geometry, which models a thin adhesive bond with a long interfacial edge crack, is similar to that which has been used to measure the dependence of interfacial toughness on crack-tip mode-mixity. The calculated effective toughness vs. applied mode-mixity relationships all display a strong dependence on applied mode-mixity with the effective toughness increasing rapidly with the magnitude of the mode-mixity. The calculated relationships also show a pronounced asymmetry with respect to the applied mode-mixity. This dependence is similar to that observed experimentally, and calculated results for a glass/epoxy interface are in good agreement with published data that was generated using a test specimen of the same type as analyzed here.     

Adhesion between elastic cylinders based on the double-Hertz model

A cohesive zone model for two-dimensional adhesive contact between elastic cylinders is developed by extending the double-Hertz model of Greenwood and Johnson (1998). In this model, the adhesive force within the cohesive zone is described by the difference between two Hertzian pressure distributions of different contact widths. Closed-form analytical solutions are obtained for the interfacial traction, deformation field and the equilibrium relation among applied load, contact half-width and the size of cohesive zone. Based on these results, a complete transition between the JKR and the Hertz type contact models is captured by defining a dimensionless transition parameter μ, which governs the range of applicability of different models. The proposed model and the corresponding analytical results can serve as an alternative cohesive zone solution to the two-dimensional adhesive cylindrical contact.     

2.5维自愈合C/SiC复合材料弹性性能预测

试验研究了2.5维自愈合C/SiC复合材料（2.5D-C/SiC）的面内弹性性能。根据复合材料的细观结构特点,建立了2.5D-C/SiC的弹性模量预测模型。计算结果与试验值吻合较好,证明了模型及计算方法的正确性。通过改变经纱编织角、纤维束中纤维数以及经纱穿过纬纱的层数,得到了经纱体积含量、纬纱体积含量、纤维总体积含量以及弹性模量随上述参数的变化规律。保持碳化硅与碳化硼总体积分数不变,改变碳化硼与碳化硅的体积比α,得到了弹性模量随α的变化关系。     

A refined cohesive zone model that accounts for inertia of cohesive zone of a moving crack

Existing cohesive zone models assume that actual fracture zone of non-zero mass can be modeled by a line segment (cohesive zone) with no mass and inertia. In the present work, a simplified mass-spring model is presented to study inertia effect of cohesive zone on a mode-I steady-state moving crack. It is showed that fracture energy predicted by the present model increases dramatically when a finite limiting crack speed is approached. Reasonable agreement with known experiments indicates that the present model has the potential to catch the inertia effect of cohesive zone which has been ignored in existing cohesive zone models and better simulate dynamic fracture at high crack speed.     

纤维复合材料界面横向拉伸分析

以单纤维十字型横向拉伸试验为研究对象,对纤维/基体界面采用弹性-软化双线性内聚力模型,建立了纤维复合材料在横向拉伸作用下界面法向失效过程的解析模型。得到了沿纤维/基体圆周界面的法向应力分布,纤维/基体界面的状态与界面承载力和单纤维复合材料承载力的关系,以及内聚力参数和试件几何尺寸对它们的影响。结果表明:纤维/基体圆周界面在脱粘前经历全部弹性及弹性+软化两种状态;当界面为弹性状态时,界面法向应力随界面强度线性增加;当界面为弹性+软化状态时,界面软化范围随界面裂纹萌生位移的增加而增大;界面初始脱粘位置与拉伸荷载方向重合;界面初始脱粘时的界面承载力随界面强度及界面裂纹萌生位移的增加而增加,随界面裂纹生成位移的增加而降低;单纤维复合材料的脱粘荷载受基体截面尺寸的影响,当纤维体积含量相同时,沿荷载方向截面尺寸的增大对提高脱粘荷载更显著。     

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.     

Cohesive fracture simulation and failure modes of FRP–concrete bonded interfaces

A work-of-fracture method using three-point bend beam (3PBB) specimen, commonly employed to determine the fracture energy of concrete, is adapted to evaluate the mode-I cohesive fracture of fiber reinforced plastic (FRP) composite–concrete adhesively bonded interfaces. In this study, a bilinear damage cohesive zone model (CZM) is used to simulate cohesive fracture of FRP–concrete bonded interfaces. The interface cohesive process damage model is proposed to simulate the adhesive–concrete interface debonding; while a tensile plastic damage model is used to account for the cohesive cracking of concrete near the bond line. The influences of the important interface parameters, such as the interface cohesive strength, concrete tensile strength, critical interface energy, and concrete fracture energy, on the interface failure modes and load-carrying capacity are discussed in detail through a numerical finite element parametric study. The results of numerical simulations indicate that there is a transition of the failure modes controlling the interface fracture process. Three failure modes in the mode-I fracture of FRP–concrete interface bond are identified: (1) complete adhesive–concrete interface debonding (a weak bond), (2) complete concrete cohesive cracking near the bond line (a strong bond), and (3) a combined failure of interface debonding and concrete cohesive cracking. With the change of interface parameters, the transition of failure modes from interface debonding to concrete cohesive cracking is captured, and such a transition cannot be revealed by using a conventional fracture mechanics-based approach, in which only an energy criterion for fracture is employed. The proposed cohesive damage models for the interface and concrete combined with the numerical finite element simulation can be used to analyze the interface fracture process, predict the load-carrying capacity and ductility, and optimize the interface design, and they can further shed new light on the interface failure modes and transition mechanism which emulate the practical application.     

内聚力模型裂纹问题分析的解析奇异单元

内聚力模型已经被广泛应用于需要考虑断裂过程区的裂纹问题当中,然而常用的数值方法应用于分析内聚力模型裂纹问题时还存在着一些不足,比如不能准确的给出断裂过程区的长度、需要网格加密等。为了克服这些缺点,论文构造了一个新型的解析奇异单元,并将之应用于基于内聚力模型的裂纹分析当中。首先将虚拟裂纹表面处的内聚力用拉格拉日插值的方法近似表示为多项式的形式,而多项式表示的内聚力所对应的特解可以被解析地给出。然后利用一个简单的迭代分析,基于内聚力模型的裂纹问题就可以被模拟出来了。最后,给出二个数值算例来证明本文方法的有效性。     

Bilinear softening parameters and equivalent LEFM R-curve in quasibrittle failure

For composites and adhesive joints, the determination of the cohesive zone parameters from Double Cantilever Beam specimens loaded with pure moments is now well established and documented. However, for quasibrittle materials used in Civil Engineering such as concrete or wood, the difficulty to apply a pure bending moment lies inappropriated the method used for composites. Nevertheless, the one-to-one correspondence which exists between the R-curve and the softening curve is here revisited and adapted for any kind of specimen geometry and for the bilinear approximation of the softening function, well-known to successfully describe the failure of a wide group of quasibrittle materials. It is shown that even though the connections between the cohesive parameters and the ‘equivalent LEFM’ R-curve are geometry and material dependent, their trends are preserved whatever the specimen geometry and the material are. The outline of a general estimation procedure of the cohesive zone parameters funded on the equivalent LEFM R-curve is proposed.     

Fracture of bicrystal metal/ceramic interfaces: A study via the mechanism-based strain gradient crystal plasticity theory

Two continuum mechanical models of crystal plasticity theory namely, conventional crystal plasticity theory and mechanism-based crystal plasticity theory, are used to perform a comparative study of stresses that are reached at and ahead of the crack tip of a bicrystal niobium/alumina specimen. Finite element analyses are done for a stationary crack tip and growing cracks using a cohesive modelling approach. Using mechanism-based strain gradient crystal plasticity theory the stresses reached ahead of the crack tip are found to be two times larger than the stresses obtained from conventional crystal plasticity theory. Results also show that strain gradient effects strongly depend on the intrinsic material length to the size of plastic zone ratio (l/R₀). It is found that the larger the (l/R₀) ratio, the higher the stresses reached using mechanism-based strain gradient crystal plasticity theory. An insight into the role of cohesive strength and work of adhesion in macroscopic fracture is also presented which can be used by experimentalists to design better bimaterials by varying cohesive strength and work of adhesion.     

建筑织物膜材双轴剪切试验与分析

为了研究建筑织物膜材的剪切力学性能,提出一种新的剪切测试方法。采用中心区域宽度和四臂长度均为16cm的十字形试件,试件纱线的经纬向与加载方向呈45°角。根据膜材变形和应力关系,推导了剪应力和剪应变的计算方法。定义了使试件中心区域产生三个循环剪应力场的加载谱,循环产生正负交替的剪切应力。试验采用表面抛光的(Polyvinylidene Fluoride,聚偏氟乙烯)涂层膜材,测量x、y两个方向的应力和应变,通过计算得到剪切应力应变曲线,并对试验结果进行分析。结果表明,新的剪切测试方法能够反映建筑织物膜材剪切力学性能。最后通过有限元方法模拟材料受剪状态下的应力和应变,与试验得到的应力值和应变值相近。     

Quasi-Static Torsion Characterization of Micro-diameter Copper Wires

A method to measure very small torques that subject micro-diameter copper wires to plasticity is developed for quasi-static torsion experiment. Following the concept in the work by Fleck et al. (Acta Metall. Mater. 42:2, 1994), we employed a glass fiber filament as the torque cell. To calculate the small torques applied on the micro-diameter copper wires, an additional rotation sensor is required to measure the rotation of the glass fiber torque cell. The rotation sensor system is attached between the glass fiber and the copper wire specimen. It uses a laser extensometer to gauge the distance between a helical and a horizontal reflection tapes on a foam cylinder, which is calibrated with the angle of rotation. A new set of torsional experimental data for the copper wires with four different diameters, from 16–180 μm, are presented. All copper wires exhibit a typical elastic-plastic response. The torsional properties of these copper wires were not found to be significantly different. The uncertainties of the measurement and analysis are discussed.     

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

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

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

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

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.     

内聚力界面单元与复合材料的界面损伤分析

推导了一种基于内聚力模型无厚的界面单元,用来模拟复合材料纤维与基体之间的界面层．研究了纤维周期分布的复合材料受横向荷载时,在界面不同的强韧性条件下其界面损伤演化的规律和对复合材料整体性质的影响     

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.