双材料界面断裂力学模型与实验方法

纤维增强聚合物(FRP)质轻、高强, 可提高结构的刚度、强度、抗震性能和耐久性, 近年来在结构加固及工程改造中得到广泛应用. FRP与传统复合材料之间形成双材料黏结界面, 界面断裂特性是决定双材料结构性能的关键因素. 对双材料界面裂纹尖端应力场理论、界面裂纹模型、黏结界面I型、II型及混合型断裂试验及理论研究现状进行综合评述和分析. 界面模型主要有经典梁/板理论和刚性节点模型、考虑剪切变形的双亚层理论和半刚性节点模型、基于双亚层理论的柔性节点模型、考虑剪切变形的多层亚层理论和多亚层柔性节点模型、弹性地基梁模型以及黏聚模型. 还介绍了双材料界面断裂力学在FRP-混凝土研究中的应用.      

双材料界面Ⅰ型扩展裂纹尖端的弹黏塑性场

双材料界面中存在材料黏性效应, 对界面裂纹尖端场的分布和界面本身性能 的变化起着重要的影响. 考虑裂纹尖端的奇异性, 建立了双材料界面扩展裂纹尖端的弹黏塑 性控制方程. 引入界面裂纹尖端的位移势函数和边界条件, 对刚性-弹黏塑性界面I型界面 裂纹进行了数值分析, 求得了界面裂纹尖端应力应变场, 并讨论了界面裂纹尖端场随各影响 参数的变化规律. 计算结果表明, 黏性效应是研究界面扩展裂纹尖端场时的一个主要因素, 界面裂纹尖端为弹黏塑性场, 其场受材料的黏性系数、马赫数和奇异性指数控制.     

双材料梁界面力学及其应用:综述

双材料结构在工程中得到了广泛应用,如薄膜涂层、压电材料、复合材料层合板和夹层板、粘结接头、FRP加固混凝土结构等;然而,这些结构的破坏通常是从界面及其附近开始的.通常利用损伤力学的方法(材料强度的方法)来预测裂纹的萌生,利用断裂力学的方法来预测裂纹的扩展;因此,开展双材料结构的界面应力分析和具分层双材料结构的断裂以及相关分析是至关重要的.首先介绍双材料梁粘结界面应力分析的基本模型,从而为预测裂纹的萌生提供了有力工具.然后综述双材料梁界面断裂力学分析的基本方法,并详细介绍解析解求解的裂纹尖端法及其相关模型,重点强调裂纹尖端变形对分析结果的影响.最后介绍界面分层对双材料结构其它力学特性的影响,例如屈曲和振动特性.     

有限尺寸下双材料界面附近裂纹位置对裂纹尖端的影响

随着复合材料的应用和发展,不同材料组成的界面结构越来越受到人们的重视。界面层两侧材料的性能相异会引起材料界面端奇异性,同时界面和界面附近存在裂纹会引起裂尖处的应力奇异性。因此双材料界面附近的力学分析是比较复杂的。本文建立双材料直角界面模型,在材料界面附近预设初始裂纹,计算了有限材料尺寸对界面应力场及其附近裂纹应力强度因子的影响。运用弹性力学中的 Goursat 公式求得直角界面端在有限尺寸下的应力场以及其应力强度系数。通过叠加原理和格林函数法进一步得到在直角界面端附近的裂纹尖端应力强度因子。计算结果表明,在适当范围内改变材料内裂纹与界面之间的距离,界面附近裂纹尖端的应力强度因子随着裂纹与界面距离的增加而减少,并且逐渐趋于稳定。分析结果可以为预测双材料结构复合材料界面失效位置提供参考。     

考虑夹杂相互作用的复合陶瓷夹杂界面的断裂分析

复合材料中夹杂含量较高时,夹杂间的相互作用能显著改变材料细观应力应变场分布,基体和夹杂中的平均应力应变水平也会发生较大变化,导致复合材料强度等力学性能发生显著变化. 为修正单一夹杂模型运用在实际材料中的误差,基于相互作用直推估计法,建立一种考虑含夹杂相互作用的夹杂界面裂纹开裂模型. 首先根据相互作用直推估计法,得到残余应力和外载应力共同作用下夹杂中的平均应力,再计算无限大基体中相同的夹杂达到相同应力场时的等效加载应力,将此加载应力作为含界面裂纹夹杂的等效应力边界条件,在此边界条件下求得界面裂纹尖端的应力强度因子,进而得到界面裂纹开裂的极限加载条件,并分析了夹杂弹性性能、含量、热残余应力、夹杂尺寸等因素对界面裂纹开裂条件的影响. 结果表明,方法能够有效修正单夹杂模型运用在实际材料中的误差,较大的残余应力对界面裂纹开裂有重要的影响,夹杂刚度的影响并非单调且比较复杂;在残余应力较小时,降低柔性夹杂刚度或者增大刚性夹杂刚度都有利于提高材料强度;扩大夹杂尺寸将导致裂纹开裂极限应力显著降低,从而降低材料强度.      

压剪载荷作用下界面裂纹尖端场的研究

建立了弹性-幂律蠕变双材料界面裂纹准静态扩展的力学模型，求得了裂纹尖端应力、应变和位移场分离变量形式的解及其数值结果；讨论了材料性能参数对裂纹尖端场的影响；计算和分析了界面裂纹的摩擦效应，并且得出了给定条件下裂尖场的轮廓图形．     

双材料自由边界面端部问题的数值分析

对各向异性双材料自由边界面端部奇异性场问题进行了研究,利用有限元分析法所得到的各向异性双材料自由边界面端部的应力奇异性指数以及角分布函数,构造了一个自由边界面端部单元,据此建立了自由边界面端部奇异性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型.与四节点单元相结合,提出一种求解自由边界面端部广义应力强度因子的杂交元法.考核例结果表明:本文方法的数值解精度高,可应用于各向异性材料双材料自由边界面端部问题.     

光弹法实测正交异性双材料界面裂纹断裂参数

采用光弹贴片法实测正交异性双材料界面裂纹尖端区域的应力应变场, 进而求出界面裂纹的断裂力学参量. 将正交异性双材料板加工成拉伸试件,在聚碳酸酯贴片 的单侧表面镀金属铝膜,以提高贴片的反射效率. 沿贴片后的双材料界面预制裂缝,逐渐加 大载荷,得到一系列清晰的等差线条纹图. 利用正交异性双材料界面裂纹尖端应力分量表达 式计算出应力强度因子. 实验表明,光弹贴片法可有效地分析正交异性双材料界面裂纹问题.     

刚性-粘弹性材料界面Ⅰ型动态扩展裂纹的尖端场

裂纹尖端渐近场的研究是断裂力学研究的重要课题之一。为了研究粘性效应作用下的界面动态扩展裂纹尖端渐近场,建立了刚性.粘弹性材料界面Ⅰ型动态扩展裂纹的力学模型;在稳态蠕变阶段,弹性变形和粘性变形同时在裂纹尖端场中占主导地位,应力和应变具有相同的奇异量级,即(σ,ε)∝r-1/(n-1)。当n→∞,幂硬化粘弹性材料动态扩展裂纹尖端场与Freund给出的理想塑性材料动态扩展裂纹尖端场具有相近的奇异量级;结合运动和协调方程,推导出粘弹性材料动态扩展裂尖场的控制方程。根据问题的边界条件和连续条件,通过数值计算,得到了裂纹尖端连续的分离变量形式的应力、应变和位移场。数值计算表明,裂纹尖端场主要受材料的蠕变指数n和马赫数M的控制,这为解决工程实践中所遇到的相应的问题和建立材料的破坏准则提供理论的参考。     

幂硬化弹塑性-刚性界面Ⅲ型定常扩展裂纹尖端的渐近场

对幂硬化弹塑性材料－刚性材料界面上裂纹以定常方式扩展的Ⅲ型问题进行弹塑性渐近分析，给出裂纹尖端的应力，应变和位移场解。通过数值计算，考察了不同Mach数以及裂纹尖端混合参数对场解的构造以及应力，应变分布的影响，为给出合理的断裂准则提供理论依据。     

Mechanics and fracture of crack tip deformable bi-material interface

Novel interface deformable bi-layer beam theory is developed to account for local effects at crack tip of bi-material interface by modeling a bi-layer composite beam as two separate shear deformable sub-layers with consideration of crack tip deformation. Unlike the sub-layer model in the literature in which the crack tip deformations under the interface peel and shear stresses are ignored and thus a “rigid” joint is used, the present study introduces two interface compliances to account for the effect of interface stresses on the crack tip deformation which is referred to as the elastic foundation effect; thus a flexible condition along the interface is considered. Closed-form solutions of resultant forces, deformations, and interface stresses are obtained for each sub-layer in the bi-layer beam, of which the local effects at the crack tip are demonstrated. In this study, an elastic deformable crack tip model is presented for the first time which can improve the split beam solution. The present model is in excellent agreements with analytical 2-D continuum solutions and finite element analyses. The resulting crack tip rotation is then used to calculate the energy release rate (ERR) and stress intensity factor (SIF) of interface fracture in bi-layer materials. Explicit closed-form solutions for ERR and SIF are obtained for which both the transverse shear and crack tip deformation effects are accounted. Compared to the full continuum elasticity analysis, such as finite element analysis, the present solutions are much explicit, more applicable, while comparable in accuracy. Further, the concept of deformable crack tip model can be applied to other bi-layer beam analyses (e.g., delamination buckling and vibration, etc.).     

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.     

Buckling of delaminated bi-layer beam-columns

An improved analytical model is presented to analyze the delamination buckling of a bi-layer beam-column with a through-the-width delamination. Both the transverse shear deformation and local delamination tip deformations are taken into consideration, and two delaminated sub-layers as well as two substrates in the intact (un-delaminated) regions are modeled as individual Timoshenko beams. A deformable interface is introduced to establish the continuity condition between the two substrates in the intact regions. Consequently, a flexible joint is formed at the delamination tip, and it is different from the conventional rigid joint given in most of studies in the literature, in which the local delamination tip deformations are completely ignored. In contrast to the local delamination buckling in our previous study (Qiao et al., 2010), the present model accounts for the global deformations of the intact region in the delaminated composite beam-column, thus capable of capturing the buckling mode shape transitions from the global, to global–local coexistent, and to local buckling for asymmetric delamination as the interface delamination increases. Good agreement of the present analytical solutions with the full 2-D elastic finite element analysis demonstrates the local deformation effects around the delamination tip and verifies the accuracy of the present model. Parametric studies are conducted to investigate the effects of loading eccentricity, delaminated sub-layer thickness ratio, and interface compliance on the critical buckling load for the delaminated composite beam-column. Transitions of buckling modes from the global to local delamination buckling are also disclosed as the thickness of one sub-layer reduces from the thick sub-layer to a thin film. The developed delamination buckling solution facilitates the design analysis and optimization of laminated composite structures, and it can be used with confidence in buckling analysis of delaminated composite structures.     

A new simple third-order shear deformation theory of plates

An improved simple third-order shear deformation theory for the analysis of shear flexible plates is presented in this paper. This new plate theory is composed of three parts: the simple third-order kinematics of displacements reduced from the higher-order displacement field derived previously by the author; a system of 10th-order differential equilibrium equations in terms of the three generalized displacements of bending plates; five boundary conditions at each edge of plate boundaries. Although the resulting displacement field is the same as that proposed by Murthy, the variational consistent governing equations and the associated proper boundary conditions are derived and identified in this work for the first time in the literature. The applications and accuracy of the present shear deformation theory of plates are demonstrated by analytically solving the differential governing equations of a twisting plate, a bending beam and two bending plates to which the 3-D elasticity solutions are available, and excellent agreements are achieved even for the torsion of a plate with square cross-section as well the local effects of stresses at plate boundaries can be characterized accurately. These analytical solutions clearly show that the simple third-order shear deformation theory developed in this work indeed gives better results than the first-order shear deformation theories and other simple higher-order shear deformation theories, since the present third-order shear flexible theory is based on a more rigorous kinematics of displacements and consists of not only a system of variational consistent differential equations, but also a group of consistent boundary conditions associated with the differential equations. The present simple third-order shear deformation theory can easily be applied to the static and dynamic finite element analysis of laminated plates just like the applications of other popular shear flexible plate theories, and improved results could be obtained from the present simple third-order shear deformable theories of plates.     

Interfacial mixed-mode fracture of adhesive bonds undergoing large deformation

Interfacial fracture of adhesive bonds undergoing large-scale yielding is studied using a combined experimental/finite-element approach. The full range of in-plane mode mixity is produced over bond thickness ranging from 30 to 500 μm using the scarf and the ENF joint geometries. Novel techniques for introducing pre-cracks and surface decoration, together with in situ observations, facilitate accurate determination of the bond-average and the local shear strains at the crack tip during the onset as well as the rest of the crack propagation event. The crack generally grew along one of the two interfaces of the bond, although the failure was always fully cohesive. The local shear strain at the crack tip is independent of the bond thickness, and, under quasi-static conditions, it remains constant throughout the growth, which make it a viable fracture parameter. This quantity strongly depends on the mode mixity, the sign of the phase angle (i.e., shearing direction) and the crack speed, however.A finite-element analysis is used to obtain the crack tip deformation field for an interface crack in adhesively bonded scarf and ENF joints. Large-strain and quasi-static conditions are assumed. A distinct material model in the fracture process zone that allows for volume change in the post-yield regime is incorporated into the analysis. The local deformation is characterized by a pair of bond-normal and tangential displacements corresponding to the nodal points adjacent to the crack tip. The critical values of these quantities are obtained when the FEM bond-average shear strain at the crack tip becomes equal to its experimental counterpart. The so defined critical local displacements, after an appropriate normalization, seem to conform to a single-valued, linear type interrelationship over the entire range of mode mixity. The fact that this relationship is independent of the bond thickness, and furthermore it encompasses both cases of positive and negative phase angles, makes it a viable candidate for characterizing mixed-mode interfacial fracture under large-deformation conditions.     

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.     

The effect of the irregular interface geometry in deformation and fracture of a steel substrate–boride coating composite

Presented in this paper is a computational analysis of the mechanisms involved in plastic deformation and fracture of a composite with coating under compressive and tensile loading. Using a steel specimen surface-hardened by diffusion borating, a role of the irregular geometry of the interface between the base material and hardened surface layer is investigated. In order to describe the mechanical behavior of the steel substrate and brittle coating, use is made of a plastic flow model including isotropic strain hardening and a fracture model, respectively. Using the Huber fracture criterion, the model takes into account the difference in the critical strength values for different types of local compressive and tensile states. It is shown that the irregular, serrated shape of the substrate–coating interface retards propagation of a longitudinal crack into this coating and prevents it from spalling under external compression of this composite. It is found out that even in the case of a simple uniaxial compression of the mesovolumes of this composite the boride “teeth” are subjected to tensile stresses, whose values are comparable with those of the external compressive load, and the direction of crack propagation and the general fracture behavior largely depend on the external loading conditions.     

Debonding of FRP-plated reinforced concrete beam,a bond-slip analysis. I. Theoretical formulation

External bonding of FRP plates or sheets has emerged as a popular method for strengthening reinforced concrete. Debonding along the FRP–concrete interface can lead to premature failure of the structure. In this study, a bond-slip model is established to study the interface debonding induced by a flexural crack in a FRP-plated concrete beam. The reinforced concrete beam and FRP plate are modeled as two linearly elastic Euler–Bernoulli beams bonded together through a thin layer of FRP–concrete interface. The interface layer is essentially modeled as a large fracture processing zone of which the stress–deformation relationship is described by a nonlinear bond-slip model. Three different bond-slip models (bi-linear, triangular and linear-damaging) are used. By dividing the debonding process into several stages, governing equations of interfacial shear and normal stresses are obtained. Closed-form solutions are then obtained for the interfacial shear and normal stresses and the deflection of the beam in each stage of debonding. In such a way, the proposed model unifies the whole debonding process, including elastic deformation, debonding initiation and growth, into one model. With such a superior feature, the proposed model provides an efficient and effective analytical tool to study FRP–concrete interface debonding.