界面裂纹萌生与扩展的分子动力学模拟

运用分子动力学模拟方法研究了裂纹在界面端处萌生与沿界面扩展的临界条件. 模拟考虑了一双相材料的3种模型，即构成90°/90°和 90°/180°夹角的两个界面端和一个界面裂纹. 模拟采用了包含原子区域与连续区域的并发型多尺度模型，即在界面端尖端和裂纹尖端附近 采用分子动力学(MD)方法，MD区域之外则按照线弹性有限元方法分析. 结果表明，在断裂启动时刻，3个模型沿界面的最大应力均达到界面理想强度；而且，其界 面能恰好足以克服界面材料的本征内聚能. 因此，界面端裂纹萌生与沿界面扩展的断裂条件可以通过界面理想强度和内聚能联系起来. 并基于模拟计算结果提出了界面断裂启动的统一准则.     

瓷修复体界面断裂行为的模拟实验研究

本文利用云纹干涉法和云纹干涉－－有限元混合法，对瓷修复体的模拟双材料模型界面断裂问题进行了实验研究。用云纹干涉和数字错位云纹干涉法测量带边裂纹的双材料四点简支梁在剪切作用下界面表面的剪应变分布及界面两侧局部表面的位移场，实验表明，由于界面两两侧材料力学性质不同，表现出界面剪切断裂问题的非称性和裂尖附近复合型断裂的特点；用云纹干涉法和有限元法相结合的混合法对粘接界面角点应力奇异性进行研究，并对角点附近应力应变场作了分析，得到了应力奇异指数与边界楔角，载荷的关系，证明了用界面应力强度因子Kf来描述界面端部区域应力分布的公式，并得到了双材料界面端部区域的应力应变分布情况。本文的实验结果为进一步研究口腔金瓷修复体界面的优化设计提供了基础，同时也说明云纹干涉法对于双材料界面断裂行为的研究是有效的。     

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

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

两相材料界面裂纹的起裂方向及断裂参量研究

采用FRANC2D软件研究了两相材料含界面裂纹的断裂特性。通过在材料界面利用CASCA手动生成界面裂纹并在裂纹尖端附近设置1/4奇异等参元,得到了界面裂纹的复应力强度因子;数值模拟得到了界面裂纹的起裂方向,并分析了界面裂纹的破坏特征。计算结果表明:1/4奇异等参元很好地描述了裂纹尖端场的1/r(1/2)奇异性,FRANC2D软件能够模拟界面裂纹的扩展方向,可以得到界面裂纹尖端的应力场和复应力强度因子,为界面裂纹的断裂特性的进一步分析提供依据。     

热释电材料问题的通解与界面裂纹

该文讨论了热释电材料中的热弹性问题的一般解，进而求解了共线界面裂纹问题．利用Ｓｔｒｏｈ方法，把热释电材料的热弹性界面裂纹问题化为一向量形式的Ｈｉｌｂｅｒｔ问题，求出这一Ｈｉｌｂｅｒｔ问题的通解，进而求得了热释电材料热弹性界面裂纹的闭合解，得到了温度、热流、位移、电势、应力和电位移的全场解，得到了裂纹张开位移及电势差的精确表达式．在此基础上，还求得了均匀热释电体中单个热弹性裂纹裂尖场，单个界面裂纹裂尖场以及点热源与界面裂纹的作用．此外，该文还对界面裂纹顶点附近的端部场作了渐近分析．     

界面裂纹的路径选择与数值模拟

利用界面断裂能和荷载混合度的概念研究界面裂纹的扩展路径，利用有限元数值方法模拟界面裂纹的扩展过程，再现界面断裂的各种几何构型．研究表明，界面的断裂能和混合度完全控制了裂纹在界面附近的扩展过程．数值预测的裂纹路径与理论结果一致     

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

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

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

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

幂硬化材料复合型界面裂纹问题

本文研究了两相幂硬化材料平面应变复合型界面裂纹问题,通过对裂端附近主奇异解的分析,得到硬化指数为 n_Ⅰ、n_Ⅱ(n_Ⅰ相似文献   

界面断裂力学简介与展望

界面断裂力学是目前国际力学界的前沿研究课题.本文对这门新兴分支学科的发展背景,近期进展及未来展望作了简要介绍.主要分析了各向同性及各向异性弹性材料间的界面裂纹;界面裂纹的力学模型;界面裂纹的弹塑性场及界面裂纹的断裂观念.讨论了界面断裂力学未来热门课题。     

三维基体裂纹与界面垂直接触的线奇异分析

使用三维断裂力学的超奇异积分方程方法,对两相材料中出现的三维基体裂纹与材料分界面垂直接触的应力线奇异问题作了理论分析,求得了线奇异界面应力的表达式     

Failure instability of a debonded interface in the presence of a main crack

The problem of fracture initiating from an edge crack in a nonhomogeneous beam made of two dissimilar linear elastic materials that are partially bonded along a common interface is studied by the strain energy density theory. The beam is subjected to three-point bending and the unbonded part of the interface is symmetrically located with regard to the applied loading. The applied load acts on the stiffer material, while the edge crack lies in the softer material. Fracture initiation from the tip of the edge crack and global instability of the composite beam are studied by considering both the local and global stationary values of the strain energy density function, dW/dV. A length parameter l defined by the relative distance between the maximum of the local and global minima of dW/dV is determined for evaluating the stability of failure initiation by fracture. Predictions on critical loads for fracture initiation from the tip of the edge crack, crack trajectories and fracture instability are made. In the analysis the load, the length of the edge crack and the length and position of the interfacial crack remained unchanged. The influence of the ratio of the moduli of elasticity of the two materials, the position of the edge crack and the width of the stiffer material on the local and global instability of the beam was examined. A general trend is that the critical load for crack initiation and fracture instability is enhanced as the width and the modulus of elasticity of the stiffer material increase. Thus, the stiffer material acts as a barrier in load transfer.     

On the path of a crack near a graded interface under large scale yielding

The trajectory of a crack lying parallel to a thin graded layer between two plastically dissimilar materials is studied using the exclusion region (ER) theory of fracture. The ER theory is a theoretical framework for surface separation within which a broad range of fracture phenomenologies can be represented. In the present study, the direction of crack advance is determined by maximizing the resolved normal-opening force on the near-tip region, whereas separation itself is governed by the intensity of plastic deformation near the tip. A computational study was undertaken using the ER theory. The special-purpose finite element analysis platform accommodates arbitrary––and a priori unknown––crack trajectories. The model problem considered herein involves two plastically dissimilar, but elastically identical, materials joined by a thin, graded interface layer. The initial crack lies parallel to the interface layer, and crack advance occurs under conditions of extensive plastic flow. It is found that the position of the initial crack relative to the interface layer has a strong influence on the fracture behavior. In general, the crack trajectories tend to curve toward the less-ductile material. Also, the presence of the interface layer leads to fracture toughnesses that significantly exceed those of either material individually for the configurations studied.     

A solution to the thermo-elastic interface crack branching in dissimilar anisotropic bi-material media

An interface crack or delamination may often branch out of the interface in a laminated composite due to thermal stresses developing around the delamination/crack tip when the media is exposed to heat flow induced by environmental events such as a sudden short-duration fire. In this paper, the thermo-elastic problem of interface crack branching in dissimilar anisotropic bi-media is studied by using the theory of Stroh’s dislocation formalism, extended to thermo-elasticity in matrix notation. Based on the complex variable method and the analytical continuation principle, the thermo-elastic interface crack/delamination problem is examined and a general solution in compact form is derived for dissimilar anisotropic bi-media. A set of Green’s functions is proposed for the dislocations (conventional dislocation and thermal dislocation/heat vortex) in anisotropic bi-media. These functions may be more suitable than those which have appeared in the literature on addressing thermo-elastic interface crack branching in dissimilar anisotropic bi-materials. Using the contour integral method, a closed form solution to the interaction between the dislocations and the interface crack is obtained. Within the scope of linear fracture mechanics, the thermo-elastic problem of interface crack branching is then solved by modelling the branched portion as a continuous distribution of dislocations. The influence of thermal loading and thermal properties on the branching behavior is examined, and criteria for predicting interface crack branching are suggested, based on the extensive numerical results from the study of various cases.     

Mode III effects on interface delamination

For crack growth along an interface between dissimilar materials the effect of combined modes I, II and III at the crack-tip is investigated. First, in order to highlight situations where crack growth is affected by a mode III contribution, examples of material configurations are discussed where mode III has an effect. Subsequently, the focus is on crack growth along an interface between an elastic-plastic solid and an elastic substrate. The analyses are carried out for conditions of small-scale yielding, with the fracture process at the interface represented by a cohesive zone model. Due to the mismatch of elastic properties across the interface the corresponding elastic solution has an oscillating stress singularity, and this solution is applied as boundary conditions on the outer edge of the region analyzed. For several combinations of modes I, II and III crack growth resistance curves are calculated numerically in order to determine the steady-state fracture toughness. For given values of K_I and K_II the minimum fracture toughness corresponds to K_III=0 in most of the range analyzed, but there is a range where the minimum occurs for a nonzero value of K_III.     

Stress intensity factor for an edge crack in two bonded dissimilar materials under convective cooling

The transient thermal stress crack problem for two bonded dissimilar materials subjected to a convective cooling on the surface containing an edge crack perpendicular to the interface is considered. The problem is solved using the principle of superposition and the uncoupled quasi-static thermoelasticity. The crack problem is formulated by applying the transient thermal stresses obtained from the uncracked medium with opposite sign on the crack surfaces to be the only external loads. Fourier integral transform is used to solve the perturbation problem resulting in a singular integral equation of Cauchy type in which the derivative of the crack surface displacement is the unknown function. The numerical results of the stress intensity factors are calculated for both the edge crack and the crack terminating at the interface using two different composite materials and illustrated as a function of time, crack length, coefficient of heat transfer, and the thickness ratio.     

Crack initiation direction from interface of bonded dissimilar media

Debonded region of an interface between two dissimilar materials are modeled as a line crack that tends to enhance the initiation of failure by fracture. Depending on the load that interacts with dissimilar materials, no a priori knowledge of how failure would initiate from an existing interface crack is assumed. By application of the strain energy density criterion, potential crack initiation sites are obtained for different biaxial loading states and materials with dissimilar properties.Numerical results are obtained for an epoxy/aluminum medium. In each case, a finite line segment of debonding is assumed. Uniform stresses are applied normal and parallel to the interface so that a biaxial load factor k determines the relative magnitude of biaxiality. Positive and negative k correspond, respectively, to applied tension and compression parallel to the interface. For a fixed ratio of the elastic moduli, crack initiation angles measured from the interface would increase with positive k and decrease with an increase of negative k. These findings are presented for different values of k. The direction of maximum yield initiation could also be determined from the stationary values of the strain density function. These locations are identified with elements that undergo excessive distortion while the possible fracture sites are assumed to coincide with regions where dilatational effects would dominate.     

Interaction between an interface crack and subinterface microcracks in metal/piezoelectric bimaterials

The J-integral analysis is presented for the interaction problem between a semi-infinite interface crack and subinterface matrix microcracks in dissimilar anisotropic materials. After deriving the fundamental solutions for an interface crack subjected to different loads and the fundamental solutions for an edge dislocation beneath the interface, the interaction problem is deduced to a system of singular integral equations with the aid of a superimposing technique. The integral equations are then solved numerically and a conservation law among three values of the J-integral is presented, which are induced from the interface crack tip, the microcracks and the remote field, respectively. The conservation law not only provides a necessary condition to confirm the numerical results derived, but also reveals that the microcrack shielding effect in such materials could be considered as a redistribution of the remote J-integral. It is this redistribution that does lead to the phenomenological shielding effect.     

与界面相交的裂纹尖端的应力奇异性分析

为了确定与结合材料的界面相交的裂纹尖端附近的应力奇异性次数，提出了一种基于最小势能原理的一维特殊有限元法，以奇异点为原点半径r0的扇形奇异区域，可以简化为一维线性领域，即一条以代表结合材料的两个自由表面为端点的线段。对该一维线性领域作网格划分，采用三节点一维等参数二次单元。数值计算结果与已有理论解的比较表明，该方法具有很高的精度和效率，最后，利用文中给出的方法，得到了各向异性结合材料中与界面以任意角相交的裂纹尖端的奇异性次数随裂纹的变化规律。     

Mode III stress intensity factors for edge-cracked circular shafts,bonded wedges,bonded half planes and DCB’s

Antiplane shear deformation of several edge-cracked geometries is considered. Analytical expressions are derived for the mode III stress intensity factor (SIF) of circular shafts with edge cracks, bonded half planes containing an interfacial edge crack, bonded wedges with an interfacial edge crack and also DCB’s. The results are extracted for simple isotropic materials as well as anisotropic materials and also bonded dissimilar materials and it is shown that the same expressions are obtained for the SIF under the same geometries but with different above-mentioned material properties. Different boundary conditions are assumed and the SIF relations are derived in each case. As the special cases, the SIF’s of the two bonded quarter planes containing an edge crack at the interface and infinite strip with a semi-infinite edge crack are extracted which coincide with the results cited in the literature.